pro redH4rg,date=date,skipsky=skipsky,skipflat=skipflat

	if not(keyword_set(date)) then date='07Apr24'
	rootdir='c:\figerdev\data\H4RG\'
	cd,rootdir+date

	if not(keyword_set(skipsky)) then skipsky=0
	if not(keyword_set(skipflat)) then skipflat=0

	if skipflat eq 0 then begin
		file='BigWhiteSpot_FlatField_H4RG_SIPIN_2_Reads_Apr25_2007_05_26_13.fits'
		dummy=readfits(file,h)
		bias=readfits('BiasFrame_H4RG_SIPIN_2_Reads_Apr25_2007_18_32_39.fits')
		writefits,'flat.fits',dummy-bias,h
	endif

	if skipsky eq 0 then begin
		files=file_search('9152*.fits')
		dummy=readfits(files(0),h)
		imsize=size(dummy)

		im=intarr(imsize(1),imsize(2),n_elements(files))
		dummy=0

		for i=0,n_elements(files)-1 do begin
			dummy=(readfits(files(i)))
			im(*,*,i)=dummy(*,*,3)+dummy(*,*,2)-(dummy(*,*,1)+dummy(*,*,0))
			dummy=0
		endfor

		medim=median(im,dimension=3)
		sxaddpar,h,'BZERO','32768'
		sxaddpar,h,'BSCALE','1'

		writefits,'medsky'+date+'.fits',medim,h

	endif else begin

		medim=readfits('medsky'+date+'.fits',hsky)

	endelse

	files=file_search('NGC_6791*.fits')
	for i=0,n_elements(files)-1 do begin
		im=readfits(files(i),h)
		diff=im(*,*,3)+im(*,*,2)-(im(*,*,1)+im(*,*,0))
		diff=diff-medim
		medvector=median(diff(0:127,*),dimension=1)
		for j=0,4095 do begin
			diff(*,j)=diff(*,j)-medvector(j)
		endfor
		writefits,strmid(files(i),0,strpos(files(i),'.fits'))+'_s.fits',diff,h
	endfor




stop
end