pro osr

; Calculates the electronics ofset (O) plus Serial Register dark current for the images

Outfile='O+Rs.out'
close,2
openw,2,outfile

cd,'', current=cdir & cd,cdir
DB="DB"
printf,2,cdir


files=findfile(cdir+'\'+DB+'\Image\*_Img',count=nfiles)
;files=t_getfil('*.*') & siz=size(files) & nfiles=siz(1)	;for lots of files
if nfiles eq 0 then print,'*** No Files Found in***  in ', DB

print,string(10b),'There are: ',sstr(nfiles),' files.'
printf,2,string(10b),'There are: ',sstr(nfiles),' files.',string(10b)

printf,2,''
printf,2,' File  Set    Time     Temp    Exp.   #    Min    Mean     Max     Lamp  Nul2 Nul3   O+SR   '
printf,2,'  #     #     (sec)    (K)     (ms)  Col   (DN)   (DN)     (DN)          (DN) (DN)   (DN)   '
printf,2,''

for i=0,nfiles-1 do begin
 if i eq nfiles then goto,exit
 d_read,files(i),h,p
 set=d_value(h,87)
 time=d_value(h,88)
 temp=d_value(h,22)
 exp=d_value(h,36)
 ncol=d_value(h,121)
 lamp=d_value(h,102)
 nul2=d_value(h,122)
 nul3=d_value(h,123)
 osr=((nul2/4.+0.125)+(nul3/4.+0.125))/2.

 type=mtype(files(i))

 ;check to see if data is saturated...
 sat=''
 maxp=max(p)	;the maximum
 meanp=mean(p)	;the average
 minp=min(p)
 if maxp ge 4095 then sat=' *'

 printf,2,i,set,time,temp,exp,ncol,minp,meanp,maxp,sat,lamp,nul2,nul3,osr,$
  format='(i4,i5,f10.2,f8.1,f7.1,i6,i7,f9.2,i7,a2,1x,a5,i5,i5,f8.2)'

endfor	;for i

exit:
print,''
print,' * = Saturated Data'

cd,cdir
;stop
close,/all
print,'Done, results are in ',cdir+'\'+outfile
end