dome = 1 if ( dome ) then begin ;{ window,xsize=678,ysize=678 window,2,xsize=600,ysize=400 endif ;} close,1 & close,13 mywin=0 fpos=lonarr(500) & ftmp=long(0) & ii=0 pmn = 0. & pmx = 40000. xmn = 0. & xmx = 2034. lfilenm = 'list' filenm=' ' close,13 & openr,13,lfilenm while ( not eof(13) ) do begin ;{ point_lun,-13,ftmp & fpos(ii)=ftmp 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) mn = min(img(200:xsz-200,200:ysz-200),max=mx) simg = rebin(img, 678,678) print,filenm,' min,max,avg ',mn,mx, $ total(float(simg(384-50:384+49,384-50:384+49))) /( 100.*100.) print, expt,lcvrv,lcvrt,' time ', ' volts ' , ' temp ' print,' ' opal = strmid(filenm,0,2) if ( expt eq 4. ) then begin tmn = 45000 & tmx = 68000 if (opal eq '5C' ) then begin tmn = 8000 & tmx = 14000 endif endif if ( expt eq 1. ) then begin tmn = 17000 & tmx = 28000 if (opal eq '5C' ) then begin tmn = 4000 & tmx = 7000 endif endif if ( expt eq 0.25 ) then begin tmn = 7000 & tmx = 10000 if (opal eq '5C' ) then begin tmn = 3000 & tmx = 4000 endif endif if ( dome ) then begin ;{ wset,0 & tv,bytscl(simg,min=tmn,max=tmx) ;; wset,2 ;; plot,[xmn,xmx],[pmn,pmx],/nodata,xstyle=1,ystyle=1 ;; oplot,img(*,2034/2), linestyle=2 ;; oplot,img(2034/2,*) ;; wset,0 ;; rdpix,simg ;;if(mywin eq 0)then begin ;{ ;;cursor,x,y,/NORMAL ;;if ( !ERR eq 2 ) then begin ;{ ;; ii=ii-1 & if (ii le 0) then ii=0 ;; point_lun,13,fpos(ii) ;;endif else begin ;}{ ;; ii=ii+1 ;;endelse ;} ;;endif;} endif;} endwhile ;} close,13 end