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. xmn = 0. & xmx = 339. lfilenm = 'cmp4list' filenm=' ' close,13 & openr,13,lfilenm imgB=intarr(339,339) imgC=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')) 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))) img = img - bias cavg = total ( img (100:2033-100 , 1017 ) ) / (2034. - 200. ) print,filenm,' ',expt, bias , cavg opal = strmid(filenm,0,2) if ( opal eq '5B' ) then begin ;{ imgB = rebin(img, 339,339) wset,0 & & erase & tv,bytscl(imgB,min=22000,max=24000) endif if ( opal eq '5C' ) then begin ;{ imgC = rebin(img, 339,339) imgC = float(imgC) wset,0 & tv,bytscl(imgC,min=2200,max=2500),0, 339 endif wset,2 plot,[xmn,xmx],[pmn,pmx],/nodata,xstyle=1,ystyle=1 oplot,imgB(*,339/2), thick=2 ;; oplot,imgC(*,339/2) * ( 7.76e-7 / 8.00e-8 )/3.91,linestyle=2,thick=2 oplot,imgC(*,339/2) * ( 7.76e-7 / 7.94e-8 )/3.91,linestyle=2,thick=2 endwhile ;} close,13 end