;.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)
mean=fltarr(1024,1024)
pixel_check=fltarr(33)

for i=0,32 do begin

;Add up all frames to get cumulative sum
cum_sum=cum_sum+fldata(*,*,i)

;Keep track of one pixel over 33 reads to see what's happening
pixel_check(i)=fldata(270,932,i)

endfor

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


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

print,My_header
print,pixel_check

;b=data(row,*,read)

;plot,b

END