function raw_summary, image
;
;
;
DIM = SIZE(IMAGE) 
NX = DIM(1) 
NY = DIM(2) 
N  = DIM(3) 
NP = FLOAT(NX)*NY 
case n of
    21: begin
                dark = image(*,*,0)
                im1 = 1         ;first of 10 M1 tuning images
                im2 = 11        ;first of 10 M2 tuning images
        end
    22: begin
                dark = image(*,*,0)
                iiobs = 1
                im1 = 2         ;first of 10 M1 tuning images
                im2 = 12        ;first of 10 M2 tuning images
        end
    27: begin
		dark = image(*,*,6)
		im1 = 7		;first of 10 M1 tuning images
		im2 = 17	;first of 10 M2 tuning images
	end
    29: begin
		DARK = (IMAGE(*,*,1) + IMAGE(*,*,23))/2
		im1 = 2		;first of 10 M1 tuning images
		im2 = 12	;first of 10 M2 tuning images
	end
    ELSE: begin
		print, 'Data cube must have 27 or 29 images'
		print, 'Only found: ', n
		print, 'Returning...
		return, ' '
	end
endcase
;
avg = fix(total(dark)/n_elements(dark))
fmt = '(i5)'
out = string(avg, format=fmt)
;
out = out + '   '
for i=im1,im1+9 do begin
    avg = fix(total(image(*,*,i)-dark)/n_elements(dark))
    out = out + string(avg, format=fmt)
end
out = out + '   '
for i=im2,im2+9 do begin
    avg = fix(total(image(*,*,i)-dark)/n_elements(dark))
    out = out + string(avg, format=fmt)
end
;
return, out
end