infil = '/md65/log/sci160k/940422_024117
infil = '/md66/log/sci160k/940420_185642
if (n_elements(img) eq 0) then begin
    read_mdi, infil, 0, img, h
    read_mdi, infil, 1, img2, h2
end
;
for iloop=0,1 do begin
    case iloop of
	0: list = [15, 69, 512, 661]			;image 0
	1: list = [20, 362, 512, 570, 840, 954]		;image 1
    endcase
    for ilist=0,n_elements(list)-1 do begin
	ii = list(ilist)
	for ismooth=0,1 do begin
	    case iloop of
		0: begin & col0 = img(ii,*) & col1 = img(ii-1,*) & col2 = img(ii+1,*) & end
		1: begin & col0 = img2(ii,*) & col1 = img2(ii-1,*) & col2 = img2(ii+1,*) & end
	    endcase
	    tit = infil + ' Image ' + strtrim(iloop,2) + '   Column ' + strtrim(ii,2)
	    ytit = 'Column ' + strtrim(ii,2) + ' (dashed is col on either side)'
	    if (ismooth eq 1) then begin
		col0 = smooth(col0*1., 15)
		col1 = smooth(col1*1., 15)
		col2 = smooth(col2*1., 15)
		tit = tit + ' (smooth 15)'
	    end
	    plot, col0, /ynoz, tit=tit, xtit='Line Number', ytit=ytit
	    oplot, col1, linestyle=1
	    oplot, col2, linestyle=2
	    plottime, 0, 0, 'JAKE2.PRO Ver 1.0' & plottime & pause
	end
    end
end
;
end