dome = 1 if ( dome ) then begin ;{ window,xsize=339,ysize=678 window,2,xsize=600,ysize=400 endif ;} close,1 & close,13 pmn = 0. & pmx = 30000. dmn = 22000. & dmx = 26000. xmn = 0. & xmx = 339. lfilenm = 'tlist' filenm=' ' close,13 & openr,13,lfilenm imgt1=intarr(339,339) imgt2=intarr(339,339) while ( not eof(13) ) do begin ;{ readf,13,filenm img = readfits(filenm,hdr) xsz = fix(sxpar (hdr , 'NAXIS1')) ysz = fix(sxpar (hdr , 'NAXIS2')) expt = float(sxpar (hdr , 'SHUTTER')) lcvrv = fix(sxpar (hdr , 'LCVR-V')) lcvrt = float(sxpar (hdr , 'LCVR-T')) print,filenm,' ',expt byteorder,img img = long(img) img = (img lt 0 ) * ( img + 32768 + 32768) + (img ge 0)*(img) ;; Use a column --> column (y-row) at x=4 bias_line = img(4,0+100:2033-100) bias = total(bias_line) / (float(n_elements(bias_line))) print,'bias ',bias img = img - bias opal = strmid(filenm,0,2) if ( expt eq 0.250 ) then begin ;{ imgt1 = rebin(img, 339,339) imgt1 = float(imgt1) *4. wset,0 & tv,bytscl(imgt1,min=dmn,max=dmx) endif if ( expt eq 1.0 ) then begin ;}{ imgt2 = rebin(img, 339,339) imgt2 = float(imgt2) wset,0 & tv,bytscl(imgt2,min=dmn,max=dmx),0,339 endif wset,2 plot,[xmn,xmx],[pmn,pmx],/nodata,xstyle=1,ystyle=1 oplot,imgt1(*,339/2), linestyle=2 oplot,imgt2(*,339/2) endwhile ;} close,13 end