pro linearize_nirc2,fitslist
;
; Procedure to linearize NIRC2 (treated as a single detector).
; Run on all images before running anything else.  
;
; fitslist	- list of FITS files in 1024x1024 format to linearize

imgsize = 1024
fitsarr = fltarr(imgsize,imgsize)
norm = fltarr(imgsize,imgsize,/nozero)
coeff = [1.001,-6.9e-6,-0.70e-10]

openr, 1, fitslist
fitsfile = ''
while not eof(1) do begin
	readf, 1, fitsfile
	fitsarr = readfits(fitsfile,header)
	coadds = sxpar(header,'COADDS')
	norm = coeff[0]+coeff[1]*fitsarr/coadds+coeff[2]*(fitsarr/coadds)^2
	fitsarr = fitsarr / norm

	time = systime()
	time = 'Linearized:  ' + time
	fxaddpar, header, 'HISTORY', time
	writefits, fitsfile, fitsarr, header
endwhile
close, 1
end