;.run mrdfits

;READ, "Enter Name of File:",File_name,format=string
;READ, "Enter Row:",row,format=long
;READ, "Enter Read:",Read,format=long

data=mrdfits('../CMOS/dark_160K_33_min_1_8.0_40.fits',0,My_header)
fldata=float(data)

cum_sum=fltarr(1024,1024)

for i=0,32 do begin

cum_sum=cum_sum+fldata(*,*,i)

endfor

mean=cum_sum/33
diff=fldata(*,*,32)-mean


mwrfits,mean,'../CMOS/mean.fits',My_header
mwrfits,mean,'../CMOS/diff.fits',My_header

print,My_header


;b=data(row,*,read)

;plot,b

END