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 = 'plist' filenm=' ' close,13 & openr,13,lfilenm imgp1=intarr(339,339) imgp2=intarr(339,339) imgp3=intarr(339,339) imgp4=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 ( lcvrv eq 1770 ) then begin ;{ print,'p1' imgp1 = rebin(img, 339,339) wset,0 & erase & tv,bytscl(imgp1,min=dmn*expt,max=dmx*expt) endif if ( lcvrv eq 2210 ) then begin ;}{ print,'p2' imgp2 = rebin(img, 339,339) wset,0 & tv,bytscl(imgp2,min=dmn*expt,max=dmx*expt),0,339*0.5 endif if ( lcvrv eq 2840 ) then begin ;}{ print,'p3' imgp3 = rebin(img, 339,339) wset,0 & tv,bytscl(imgp3,min=dmn*expt,max=dmx*expt),0,339 endif if ( lcvrv eq 4450 ) then begin ;}{ print,'p4' imgp4 = rebin(img, 339,339) wset,0 & tv,bytscl(imgp4,min=dmn*expt,max=dmx*expt),0,339*1.5 endif wset,2 if( expt eq 0.25 ) then begin & pmn = 5000. & pmx = 7000. & endif if( expt eq 1.0 ) then begin & pmn = 5000.*4. & pmx = 7000.*4. & endif if( expt eq 4.0 ) then begin & pmn = 5000.*16. & pmx = 7000.*16. & endif plot,[xmn,xmx],[pmn,pmx],/nodata,xstyle=1,ystyle=1 oplot,imgp1(*,339/2), linestyle=2 oplot,imgp2(*,339/2) oplot,imgp3(*,339/2) oplot,imgp4(*,339/2) endwhile ;} close,13 end