;.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)
	if i ne 32 then begin
	difference=fldata(*,*,(i+1))-fldata(*,*,i)
	poss_cos=where(difference gt 10*fldata(*,*,i))
	print,poss_cos
	endif 

;Keep track of one pixel over 33 reads to see what's happening
pixel_check(i)=fldata(910,286,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
print,poss_cos

;b=data(row,*,read)

;plot,b

END