;.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_2.0_30.fits',0,My_header)
fldata=float(data)
poss_cos=0
poss_cos_oned=0
cum_sum=fltarr(1024,1024)
difference=fltarr(1024,1024)
my_mean=fltarr(1024,1024)
pixel_check=fltarr(33)

for i=5,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_oned=where(difference gt 10000)
		;poss_cos_oned=where(difference gt 10000)
		if poss_cos_oned(0) ne -1 then begin
			poss_cos=array_indices(difference,poss_cos_oned)
		;	print,poss_cos
		endif
	;print,poss_cos
	endif 

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

print,pixel_check(i)
endfor

my_mean=cum_sum/33
diff=fldata(*,*,32)-my_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