pro cg_Get_FUV_99_Lab_Calibration,f_wavelength,f_calibration,f_calibration_error,slit_width,input_spectrum,bin_definition
if n_params() le 0 then begin
;Print,'Pro Get_FUV_99_Lab_Calibration,f_wavelength,f_calibration,f_calibration_error,slit_width,input_spectrum,bin_definition'
;Print,'1999 FUV Channel Laboratory Calibration including adjustment to fit H2 model electron lamp spectrum'
;Print,' '
;Print,'f_wavelength is the wavelength vector
;Print,' '
;Print,'f_calibration is the data MULTIPLIER that converts couts/sec to kilorayleighs/pixel (or kilorayleighs/Angstrom)'
;Print,' '
;Print,'f_calibration_error is the 1-sigma error in f_calibration'
;Print,' '
;print,'slit_width selects ocultation (0), low_resolution (1), or high_resolution (2)'
;print,'
;print,'input_spectrum selects calibration for continuous spectra (1) or for discrete lines (2) '
;print,' '
;print,'bin_definition is a two element array that defines the spectralxspatial format'
;print,'for the wavelength scale and the calibration vector'
return
end
;print,' '
;print,'*********************************************************************************************************'
;Print,'** 1999 FUV Channel Laboratory Calibration including adjustment to fit H2 model electron lamp spectrum **'
;*********************************************
;set up calibration matrix parameters and format
; slit width
if n_params() lt 4 then slit_width=1 ;default is low resolution slit
if (slit_width ne 1) and (slit_width ne 2) and (slit_width ne 0) then read,'*******  Select occultation (0), low resolution (1), or high resolution entrance slit ',slit_width
;if slit_width eq 1 then print,'** Calibration values for low_resolution (0.15mm wide) spectrograph slit                               **'
;if slit_width eq 2 then print,'** Calibration values for hi_resolution (0.075mm wide) spectrograph slit                               **'
;if slit_width eq 0 then print,'** Calibration values for occultation (0.80mm wide) spectrograph slit                               **'
if slit_width eq 0 then slit_width=1.5/8
; input spectrum
if n_params() lt 5 then input_spectrum=1 ;default is continuous
if (input_spectrum ne 1) and (input_spectrum ne 2) then read,'*******  Select either a continus input spectrum (1) or discrete line (2) ',input_spectrum
;if input_spectrum eq 2 then Print,'** Multiply the Data (counts/pixel/sec) by f_calibration to obtain kiloRayleighs/pixel                 **'
;if input_spectrum eq 1 then Print,'** Multiply the Data (counts/pixel/sec) by f_calibration to obtain kiloRayleighs/Angstrom              **'
; bin_definition
if n_params() lt 6 then bin_definition=[1,1]  ;defalut is spectral & spatial binning=1
if n_elements(bin_definition) ne 2 then begin
	bin_definition=intarr(2)
	read,'*******  input spectral binning parameter (1=no binning) ' ,temp
	bin_definition(0)=temp
	read,'*******  input spatial binning parameter ((1=no binning) ' ,temp
	bin_definition(1)=temp
end
spectral_bin=bin_definition(0)
spatial_bin=bin_definition(1)
;if (n_params() eq 6) and (n_elements(bin_definition) eq 2) then spectral_bin=bin_definition(0)
;if (n_params() eq 6) and (n_elements(bin_definition) eq 2) then spatial_bin=bin_definition(1)
;Print,'** Spectral x Spatial Binning is ', strtrim(string(spectral_bin),2),' x ',strtrim(string(spatial_bin),2),'                                                                 **'
;print,'*********************************************************************************************************'
;print,' '
;**************************
;read calibration data file
get_lun,uu
openr,uu,'FUV_1999_Lab_Cal.dat'
x=' '
readf,uu,x
dat=fltarr(13*6)
readf,uu,dat
free_lun,uu
wcal=fltarr(26)
fcalerror=wcal
fcal=wcal
for k=0,12 do wcal(k)=dat(k*6)
for k=0,12 do fcal(k)=dat(k*6+1)
for k=0,12 do fcalerror(k)=dat(k*6+2)
for k=13,25 do wcal(k)=dat((k-13)*6+3)
for k=13,25 do fcal(k)=dat((k-13)*6+4)
for k=13,25 do fcalerror(k)=dat((k-13)*6+5)
close,1
;*********************
;generate calibration matrix
cg_f_flight_wavelength,wf
f_cal=fltarr(1024,64)
f_cal_error=fltarr(1024,64)
fcal=interpol(fcal,wcal,wf)/60.
fcalerror=interpol(fcalerror,wcal,wf)/60
for k=2,61 do f_cal(*,k)=fcal            ; calibration per spatial pixel at the detector center row
for k=2,61 do f_cal_error(*,k)=fcalerror  ; error per spatial pixel at the detector center row
;apply slit width correction
f_cal=f_cal/slit_width
f_cal_error=f_cal_error/slit_width
;apply factor for continuous input spectrum
if input_spectrum eq 1 then begin
	f_cal=f_cal*.78
	f_cal_error=f_cal_error*.78
end
;perform spatial binning
f_calibration_temp=fltarr(1024,64/spatial_bin)
f_calibration_error_temp=fltarr(1024,64/spatial_bin)
;check for spatial binning
if spatial_bin eq 1 then  begin  ; set up tempory binning vectors
	f_calibration_temp=f_cal
	f_calibration_error_temp=f_cal_error
end
if spatial_bin eq 64 then  begin   ; bin everything
	f_calibration_temp(*,0)=total(f_cal,2)
	f_calibration_error_temp(*,0)=total(f_cal_error,2)
end
if (spatial_bin eq 64) or (spatial_bin eq 1) then goto,spect_bin
for k=0,64/spatial_bin-1 do begin
	f_calibration_temp(*,k)=total(f_cal(*,k*spatial_bin:(k+1)*spatial_bin-1),2)
	f_calibration_error_temp(*,k)=total(f_cal_error(*,k*spatial_bin:(k+1)*spatial_bin-1),2)
end
;perform spectral binning
spect_bin: f_wavelength=wf ;define final vectors for  binning in spatial dimension is 1
f_calibration=f_calibration_temp
f_calibration_error=f_calibration_error_temp
if spectral_bin eq 1 then goto,fin
f_wavelength=fltarr(1024/spectral_bin)  ;set up vectors for spectral binning
f_calibration=fltarr(1024/spectral_bin,64/spatial_bin)
f_calibration_error=fltarr(1024/spectral_bin,64/spatial_bin)
if spectral_bin eq 1024 then begin  ; bin everything in spectral dimension
	f_calibration(0,*)=total(f_calibration_temp,1)
	f_calibration_error(0,*)=total(f_calibration_error_temp,1)
	f_wavelength=total(wf)/1024
end
if (spectral_bin eq 1024) and (spatial_bin eq 64) then begin  ;bin everything in both spatial and spectral dimensions
	f_calibration(0,0)=total(f_calibration_temp,1)
	f_calibration_error(0,0)=total(f_calibration_error_temp,1)
	f_wavelength=total(wf)/1024
end
if spectral_bin eq 1024 then goto,fin
num=spectral_bin
if spatial_bin eq 64 then begin  ;bin in spectral when binning in spatial=64
	for k=0,1024/num-1 do begin
		f_wavelength(k)=total(wf(k*num:(k+1)*num-1))/num
		f_calibration(k)=total(f_calibration_temp(k*num:(k+1)*num-1))
		f_calibration_error(k)=total(f_calibration_error_temp(k*num:(k+1)*num-1))
	end
end
if spatial_bin eq 64 then goto,fin
for k=0,1024/num-1 do f_wavelength(k)=total(wf(k*num:(k+1)*num-1))/num	;bin in spectral when binning in spatial ne 64 or 1
for k=0,1024/num-1 do begin
	f_calibration(k,*)=total(f_calibration_temp(k*num:(k+1)*num-1,*),1)
	f_calibration_error(k,*)=total(f_calibration_error_temp(k*num:(k+1)*num-1,*),1)
end
fin: if n_elements(f_calibration eq 1) then begin
	f_calibration=1./f_calibration
	f_calibration_error=1./f_calibration_error
	return
end
cal_temp=(f_calibration eq 0)*1.e10+f_calibration
err_temp=(f_calibration eq 0)*1.e10+f_calibration_error
f_calibration=1./cal_temp
f_calibration_error=err_temp/cal_temp/cal_temp
f_calibration=(f_calibration gt 1.e-9)*f_calibration
f_calibration_error=(f_calibration_error gt 1.e-9)*f_calibration_error
end
return
end