rbsp_efw_cal_spec.pro (Documentation for /home/davin/idl/socware/)

;+
; NAME: rbsp_efw_cal_spec
;
; PURPOSE: Calibrate RBSP spectral data
;
; NOTES: Intended to be called from rbsp_load_efw_spec.pro
;
; TO DO: 
;		-No accounting for possible midday source change
;		-not using boom shorting factor
;
;		There are 7 returned spec channels. The possibilities are
;			Select 7 of: E12dc,E34dc,E56dc
;						 E12ac,E34ac,E56ac
;						 Edcpar,Edcprp
;						 Eacpar,Eacprp
;						 V1ac,V2ac,V3ac,V4ac,V5ac,V6ac
;						 SCMU,SCMV,SCMW
;						 SCMpar,SCMprp,
;						 (V1dc+V2dc+V3dc+V4dc)/4,
;						 Edcprp2, Eacprp2, SCMprp2
;		
;
; Calibration procedure
;
;	subtract offset
;		vals = (double(vals) - offset[j])
;   Compensate for 4 bit dynamic range compression
;       vals /= 16.
;	Divide the 0th bin by 4 (from weirdness doc)
;		vals[*,0] /= 4.
;	Compensate for using Hanning window
;		vals /= 0.375^2
;	Apply ADC conversion
;		vals *= adc_factor^2
;	Apply channel gain
;		vals *= gain[j]^2   (V^2/Hz at this point)
;	Apply boom length to change from V^2/Hz to mV/m^2/Hz				
;		vals *=  conversion[j]^2   ;mV/m ^2 /Hz
;
;	Apply freq dependent gain
;			
;
;
;
; HISTORY:
;   2012-09-xx: Created by Aaron W Breneman, University of Minnesota
;	2013-04-xx: added automated channel selection
;	2013-04-23: added freq-dependent gain
;	2013-05-15: freq-dependent gain curves now obtained from rbsp_efw_get_gain_results.pro and
;				are no longer "hard-wired" into this code. Also added the full 
;				freq-dependent correction as well as the mu-metal square can (nT/v)
;				curves to correct the SCM data.
;
;
; VERSION:
; $LastChangedBy:  $
; $LastChangedDate:  $
; $LastChangedRevision:  $
; $URL: svn+ssh://thmsvn@ambrosia.ssl.berkeley.edu/repos/ssl_general/trunk/missions/rbsp/efw/rbsp_efw_cal_fbk.pro $
;
;-


pro rbsp_efw_cal_spec, probe=probe, trange=trange, datatype=datatype, pT=pT


;Spec sources
get_data,'rbsp'+probe+'_efw_spec_64_ccsds_data_DFB_config',data=dfbc
dfbconfig = rbsp_get_efw_dfb_config(dfbc.y)
src = dfbconfig[0].spec_config.spec_src





;***********GET THE CALIBRATION FILE*****************
;path = '~/Desktop/code/Aaron/RBSP/TDAS_trunk_svn/ssl_general/missions/rbsp/efw/calibration_files/'
;restore,path + 'gain_curves.idl'
   
;plot,tmp.freqs,tmp.gain


;
	;rbsp_efw_init
;
;	rbsp_burst=!rbsp_efw
;	rbsp_burst.remote_data_dir='http://tetra.space.umn.edu/rbspdata/'
;	
;	dt=1.e-6
;	
;	if keyword_set(probe) then p_var=probe else p_var='*'
;	vprobes = ['a','b']
;	p_var = strfilter(vprobes, p_var ,delimiter=' ',/string)
;	
;	for p=0,size(p_var,/n_elements)-1 do begin
;	
;		rbspx = 'rbsp'+ p_var[p]
;	
;		for b=1,2 do begin
;	
;			bid='vb'+string(b,format='(I0)')
;		
;			format = rbspx + '/'+bid+'_playback/YYYY/'+ $
;					rbspx+'_efw_'+bid+'_playback_YYYYMMDD_v*.txt'
;
;			relpathnames = file_dailynames(file_format=format,trange=trange, $
;					addmaster=addmaster)
;
;			files=file_retrieve(relpathnames,/last_version,_extra=rbsp_burst)
;

;***********GET THE CALIBRATION FILE*****************





if datatype eq 'spec' then begin

	compile_opt idl2
	
	;dprint,verbose=verbose,dlevel=4,'$Id: rbsp_efw_cal_filterbank.pro 10920 2012-09-17 17:08:00Z aaron_breneman $'
	
	if ~keyword_set(trange) then trange = timerange()
	
	cp0 = rbsp_efw_get_cal_params(trange[0])
	;Get gain(f) parameters
	fcals = rbsp_efw_get_gain_results()
	;Get boom lengths
	bl = rbsp_efw_boom_deploy_history(trange[0])
	
	
	; first see if we have 36, 64, or 112 bins
	names=''
	nn=''
	names=tnames('rbsp'+probe+'*spec?')
	if strlen(names[0]) ne 0 then nn=strmid(names[0],10,1)

	case nn of
		'3': bins = '36'
		'6': bins = '64'
		'1': bins = '112'
		else: dprint,'NOT CORRECT NUMBER OF BINS...CALIBRATION ABORTED**************'
	endcase

	case strlowcase(probe) of
	  'a': cp = cp0.a
	  'b': cp = cp0.b
	  else: dprint, 'INVALID PROBE NAME. CALIBRATION ABORTED*************'
	endcase
	

	adc_factor = 2.5d / 32767.5d
	
;	boom_shorting_factor = cp.boom_shorting_factor
	
	rbspx = 'rbsp' + probe[0]
	
	
	tvar = rbspx + ['_efw_64_spec0',$
					'_efw_64_spec1',$
					'_efw_64_spec2',$
					'_efw_64_spec3',$
					'_efw_64_spec4',$
					'_efw_64_spec5',$
					'_efw_64_spec6']


	gain = replicate(0.,n_elements(src))
	offset = replicate(0.,n_elements(src))
	conversion = replicate(0.,n_elements(src))
	units = replicate('',n_elements(src))
	zlims = replicate(0d,n_elements(src),2)


	if ~keyword_set(pT) then pTf = 1. else pTf = 1000.
	if ~keyword_set(pT) then pTnT = 'nT' else pTnT = 'pT'
	if ~keyword_set(pT) then pTnT_zscale = 1. else pTnT_zscale=1d3


	;frequency bins
	fbins_36L = fcals.cal_spec.freq_spec36L
	fbins_36H = fcals.cal_spec.freq_spec36H
	fbins_36C = fcals.cal_spec.freq_spec36C
	fbins_64L = fcals.cal_spec.freq_spec64L
	fbins_64H = fcals.cal_spec.freq_spec64H
	fbins_64C = fcals.cal_spec.freq_spec64C
	fbins_112L = fcals.cal_spec.freq_spec112L
	fbins_112H = fcals.cal_spec.freq_spec112H
	fbins_112C = fcals.cal_spec.freq_spec112C



	if bins eq '36' then gainf = replicate(0.,36,7)
	if bins eq '64' then gainf = replicate(0.,64,7)
	if bins eq '112' then gainf = replicate(0.,112,7)
	

	;Determine parameters for each channel
	for bb=0,6 do begin
		case src[bb] of
	
			'E12AC': begin
						gain[bb] = cp.adc_gain_eac[0]
						if bins eq '36' then gainf[*,bb] = replicate(1.,36)
						if bins eq '64' then gainf[*,bb] = replicate(1.,64)
						if bins eq '112' then gainf[*,bb] = replicate(1.,112)
						offset[bb] = cp.adc_offset_eac[0]
;						conversion[bb] = 1000./cp.boom_length[0]
						if probe eq 'a' then conversion[bb] = 1000./bl.A12
						if probe eq 'b' then conversion[bb] = 1000./bl.B12
						units[bb] = '(mV/m)!U2!N'
						zlims[bb,*] = [1d-3^2,1d-1^2]						
					 end
			'E34AC': begin
						gain[bb] = cp.adc_gain_eac[1]
						if bins eq '36' then gainf[*,bb] = replicate(1.,36)
						if bins eq '64' then gainf[*,bb] = replicate(1.,64)
						if bins eq '112' then gainf[*,bb] = replicate(1.,112)
						offset[bb] = cp.adc_offset_eac[1]
;						conversion[bb] = 1000./cp.boom_length[1]
						if probe eq 'a' then conversion[bb] = 1000./bl.A34
						if probe eq 'b' then conversion[bb] = 1000./bl.B34
						units[bb] = '(mV/m)!U2!N'
						zlims[bb,*] = [1d-3^2,1d-1^2]
					 end
			'E56AC': begin
						gain[bb] = cp.adc_gain_eac[2]
						if bins eq '36' then gainf[*,bb] = replicate(1.,36)
						if bins eq '64' then gainf[*,bb] = replicate(1.,64)
						if bins eq '112' then gainf[*,bb] = replicate(1.,112)
						offset[bb] = cp.adc_offset_eac[2]
;						conversion[bb] = 1000./cp.boom_length[2]
						if probe eq 'a' then conversion[bb] = 1000./bl.A56
						if probe eq 'b' then conversion[bb] = 1000./bl.B56
						units[bb] = '(mV/m)!U2!N'
						zlims[bb,*] = [1d-3^2,1d-1^2]
					 end
			'E12DC': begin
						gain[bb] = cp.adc_gain_edc[0]
						if probe eq 'a' then begin
							if bins eq '36' then gainf[*,bb] = fcals.cal_spec.e12DC_a_spec36.gain_vs_freq
							if bins eq '64' then gainf[*,bb] = fcals.cal_spec.e12DC_a_spec64.gain_vs_freq
							if bins eq '112' then gainf[*,bb] = fcals.cal_spec.e12DC_a_spec112.gain_vs_freq
						endif else begin
							if bins eq '36' then gainf[*,bb] = fcals.cal_spec.e12DC_b_spec36.gain_vs_freq
							if bins eq '64' then gainf[*,bb] = fcals.cal_spec.e12DC_b_spec64.gain_vs_freq
							if bins eq '112' then gainf[*,bb] = fcals.cal_spec.e12DC_b_spec112.gain_vs_freq
						endelse							
						offset[bb] = cp.adc_offset_edc[0]
;						conversion[bb] = 1000./cp.boom_length[0]
						if probe eq 'a' then conversion[bb] = 1000./bl.A12
						if probe eq 'b' then conversion[bb] = 1000./bl.B12
						units[bb] = '(mV/m)!U2!N'
						zlims[bb,*] = [1d-3^2,1d-1^2]
					 end
			'E34DC': begin
						gain[bb] = cp.adc_gain_edc[1]
						if probe eq 'a' then begin
							if bins eq '36' then gainf[*,bb] = fcals.cal_spec.e34DC_a_spec36.gain_vs_freq
							if bins eq '64' then gainf[*,bb] = fcals.cal_spec.e34DC_a_spec64.gain_vs_freq
							if bins eq '112' then gainf[*,bb] = fcals.cal_spec.e34DC_a_spec112.gain_vs_freq
						endif else begin
							if bins eq '36' then gainf[*,bb] = fcals.cal_spec.e34DC_b_spec36.gain_vs_freq
							if bins eq '64' then gainf[*,bb] = fcals.cal_spec.e34DC_b_spec64.gain_vs_freq
							if bins eq '112' then gainf[*,bb] = fcals.cal_spec.e34DC_b_spec112.gain_vs_freq
						endelse
						offset[bb] = cp.adc_offset_edc[1]
;						conversion[bb] = 1000./cp.boom_length[1]
						if probe eq 'a' then conversion[bb] = 1000./bl.A34
						if probe eq 'b' then conversion[bb] = 1000./bl.B34
						units[bb] = '(mV/m)!U2!N'
						zlims[bb,*] = [1d-3^2,1d-1^2]
					 end
			'E56DC': begin
						gain[bb] = cp.adc_gain_edc[2]
						if probe eq 'a' then begin
							if bins eq '36' then gainf[*,bb] = fcals.cal_spec.e56DC_a_spec36.gain_vs_freq
							if bins eq '64' then gainf[*,bb] = fcals.cal_spec.e56DC_a_spec64.gain_vs_freq
							if bins eq '112' then gainf[*,bb] = fcals.cal_spec.e56DC_a_spec112.gain_vs_freq
						endif else begin
							if bins eq '36' then gainf[*,bb] = fcals.cal_spec.e56DC_b_spec36.gain_vs_freq
							if bins eq '64' then gainf[*,bb] = fcals.cal_spec.e56DC_b_spec64.gain_vs_freq
							if bins eq '112' then gainf[*,bb] = fcals.cal_spec.e56DC_b_spec112.gain_vs_freq
						endelse
						offset[bb] = cp.adc_offset_edc[2]
;						conversion[bb] = 1000./cp.boom_length[2]
						if probe eq 'a' then conversion[bb] = 1000./bl.A56
						if probe eq 'b' then conversion[bb] = 1000./bl.B56
						units[bb] = '(mV/m)!U2!N'
						zlims[bb,*] = [1d-3^2,1d-1^2]
					 end
			'V1DC': begin
						gain[bb] = cp.adc_gain_vdc[0]
						if bins eq '36' then gainf[*,bb] = replicate(1.,36)
						if bins eq '64' then gainf[*,bb] = replicate(1.,64)
						if bins eq '112' then gainf[*,bb] = replicate(1.,112)
						offset[bb] = cp.adc_offset_vdc[0]
						conversion[bb] = 1.
						units[bb] = 'volts!U2!N'
						zlims[bb,*] = [1d-4^2,1d-2^2]
					end
			'V2DC': begin
						gain[bb] = cp.adc_gain_vdc[1]
						if bins eq '36' then gainf[*,bb] = replicate(1.,36)
						if bins eq '64' then gainf[*,bb] = replicate(1.,64)
						if bins eq '112' then gainf[*,bb] = replicate(1.,112)
						offset[bb] = cp.adc_offset_vdc[1]
						conversion[bb] = 1.
						units[bb] = 'volts!U2!N'
						zlims[bb,*] = [1d-4^2,1d-2^2]
					end
			'V3DC': begin
						gain[bb] = cp.adc_gain_vdc[2]
						if bins eq '36' then gainf[*,bb] = replicate(1.,36)
						if bins eq '64' then gainf[*,bb] = replicate(1.,64)
						if bins eq '112' then gainf[*,bb] = replicate(1.,112)
						offset[bb] = cp.adc_offset_vdc[2]
						conversion[bb] = 1.
						units[bb] = 'volts!U2!N'
						zlims[bb,*] = [1d-4^2,1d-2^2]
					end
			'V4DC': begin
						gain[bb] = cp.adc_gain_vdc[3]
						if bins eq '36' then gainf[*,bb] = replicate(1.,36)
						if bins eq '64' then gainf[*,bb] = replicate(1.,64)
						if bins eq '112' then gainf[*,bb] = replicate(1.,112)
						offset[bb] = cp.adc_offset_vdc[3]
						conversion[bb] = 1.
						units[bb] = 'volts!U2!N'
						zlims[bb,*] = [1d-4^2,1d-2^2]
					end
			'V5DC': begin
						gain[bb] = cp.adc_gain_vdc[4]
						if bins eq '36' then gainf[*,bb] = replicate(1.,36)
						if bins eq '64' then gainf[*,bb] = replicate(1.,64)
						if bins eq '112' then gainf[*,bb] = replicate(1.,112)
						offset[bb] = cp.adc_offset_vdc[4]
						conversion[bb] = 1.
						units[bb] = 'volts!U2!N'
						zlims[bb,*] = [1d-4^2,1d-2^2]
					end
			'V6DC': begin
						gain[bb] = cp.adc_gain_vdc[5]
						if bins eq '36' then gainf[*,bb] = replicate(1.,36)
						if bins eq '64' then gainf[*,bb] = replicate(1.,64)
						if bins eq '112' then gainf[*,bb] = replicate(1.,112)
						offset[bb] = cp.adc_offset_vdc[5]
						conversion[bb] = 1.
						units[bb] = 'volts!U2!N'
						zlims[bb,*] = [1d-4^2,1d-2^2]
					end

			'V1AC': begin
						gain[bb] = cp.adc_gain_vac[0]
						if bins eq '36' then gainf[*,bb] = replicate(1.,36)
						if bins eq '64' then gainf[*,bb] = replicate(1.,64)
						if bins eq '112' then gainf[*,bb] = replicate(1.,112)
						offset[bb] = cp.adc_offset_vac[0]
						conversion[bb] = 1.
						units[bb] = 'volts!U2!N'
						zlims[bb,*] = [1d-4^2,1d-2^2]
					end
			'V2AC': begin
						gain[bb] = cp.adc_gain_vac[1]
						if bins eq '36' then gainf[*,bb] = replicate(1.,36)
						if bins eq '64' then gainf[*,bb] = replicate(1.,64)
						if bins eq '112' then gainf[*,bb] = replicate(1.,112)
						offset[bb] = cp.adc_offset_vac[1]
						conversion[bb] = 1.
						units[bb] = 'volts!U2!N'
						zlims[bb,*] = [1d-4^2,1d-2^2]
					end
			'V3AC': begin
						gain[bb] = cp.adc_gain_vac[2]
						if bins eq '36' then gainf[*,bb] = replicate(1.,36)
						if bins eq '64' then gainf[*,bb] = replicate(1.,64)
						if bins eq '112' then gainf[*,bb] = replicate(1.,112)
						offset[bb] = cp.adc_offset_vac[2]
						conversion[bb] = 1.
						units[bb] = 'volts!U2!N'
						zlims[bb,*] = [1d-4^2,1d-2^2]
					end
			'V4AC': begin
						gain[bb] = cp.adc_gain_vac[3]
						if bins eq '36' then gainf[*,bb] = replicate(1.,36)
						if bins eq '64' then gainf[*,bb] = replicate(1.,64)
						if bins eq '112' then gainf[*,bb] = replicate(1.,112)
						offset[bb] = cp.adc_offset_vac[3]
						conversion[bb] = 1.
						units[bb] = 'volts!U2!N'
						zlims[bb,*] = [1d-4^2,1d-2^2]
					end
			'V5AC': begin
						gain[bb] = cp.adc_gain_vac[4]
						if bins eq '36' then gainf[*,bb] = replicate(1.,36)
						if bins eq '64' then gainf[*,bb] = replicate(1.,64)
						if bins eq '112' then gainf[*,bb] = replicate(1.,112)
						offset[bb] = cp.adc_offset_vac[4]
						conversion[bb] = 1.
						units[bb] = 'volts!U2!N'
						zlims[bb,*] = [1d-4^2,1d-2^2]
					end
			'V6AC': begin
						gain[bb] = cp.adc_gain_vac[5]
						if bins eq '36' then gainf[*,bb] = replicate(1.,36)
						if bins eq '64' then gainf[*,bb] = replicate(1.,64)
						if bins eq '112' then gainf[*,bb] = replicate(1.,112)
						offset[bb] = cp.adc_offset_vac[5]
						conversion[bb] = 1.
						units[bb] = 'volts!U2!N'
						zlims[bb,*] = [1d-4^2,1d-2^2]
					end

			'SCMU': begin
						gain[bb] = cp.adc_gain_msc[0]
						;*******************
						;***TEMPORARY - rbsp_efw_get_cal_params HAS A ROUGH SQ-CAN FACTOR
						;BUILT IN. I GET THE CORRECT CURVES FROM rbsp_efw_get_gain_results. SO, I NEED
						;TO MANUALLY REMOVE THE CORRECTION HERE. 
						;if probe eq 'a' then gain[bb] /= 2.35 else gain[bb] /= 2.33
						;*******************

						if probe eq 'a' then begin
							if bins eq '36' then gain_tmp = fcals.cal_spec.scmu_a_spec36.gain_vs_freq
							if bins eq '64' then gain_tmp = fcals.cal_spec.scmu_a_spec64.gain_vs_freq
							if bins eq '112' then gain_tmp = fcals.cal_spec.scmu_a_spec112.gain_vs_freq
			  				;nT per volt square can conversion curve
							if bins eq '36' then nT_volt = fcals.cal_spec.scmu_a_spec36.stimcoil_nt2v
							if bins eq '64' then nT_volt = fcals.cal_spec.scmu_a_spec64.stimcoil_nt2v
							if bins eq '112' then nT_volt = fcals.cal_spec.scmu_a_spec112.stimcoil_nt2v
						endif else begin
							if bins eq '36' then gain_tmp = fcals.cal_spec.scmu_b_spec36.gain_vs_freq
							if bins eq '64' then gain_tmp = fcals.cal_spec.scmu_b_spec64.gain_vs_freq
							if bins eq '112' then gain_tmp = fcals.cal_spec.scmu_b_spec112.gain_vs_freq
			  				;nT per volt square can conversion curve
							if bins eq '36' then nT_volt = fcals.cal_spec.scmu_b_spec36.stimcoil_nt2v
							if bins eq '64' then nT_volt = fcals.cal_spec.scmu_b_spec64.stimcoil_nt2v
							if bins eq '112' then nT_volt = fcals.cal_spec.scmu_b_spec112.stimcoil_nt2v							
						endelse

						gainf[*,bb] = 1/(gain_tmp * nT_volt)
						offset[bb] = cp.adc_offset_msc[0]
						conversion[bb] = pTf
						units[bb] = pTnT+'!U2!N'
						zlims[bb,*] = [(pTnT_zscale*1d-4)^2,(pTnT_zscale*1d-2)^2]
					end
			'SCMV': begin
						gain[bb] = cp.adc_gain_msc[1]
						;*******************
						;***TEMPORARY - rbsp_efw_get_cal_params HAS A ROUGH SQ-CAN FACTOR
						;BUILT IN. I GET THE CORRECT CURVES FROM rbsp_efw_get_gain_results. SO, I NEED
						;TO MANUALLY REMOVE THE CORRECTION HERE. 
						;if probe eq 'a' then gain[bb] /= 2.22 else gain[bb] /= 2.10
						;*******************
						if probe eq 'a' then begin
							if bins eq '36' then gain_tmp = fcals.cal_spec.scmv_a_spec36.gain_vs_freq
							if bins eq '64' then gain_tmp = fcals.cal_spec.scmv_a_spec64.gain_vs_freq
							if bins eq '112' then gain_tmp = fcals.cal_spec.scmv_a_spec112.gain_vs_freq
			  				;nT per volt square can conversion curve
							if bins eq '36' then nT_volt = fcals.cal_spec.scmv_a_spec36.stimcoil_nt2v
							if bins eq '64' then nT_volt = fcals.cal_spec.scmv_a_spec64.stimcoil_nt2v
							if bins eq '112' then nT_volt = fcals.cal_spec.scmv_a_spec112.stimcoil_nt2v
						endif else begin
							if bins eq '36' then gain_tmp = fcals.cal_spec.scmv_b_spec36.gain_vs_freq
							if bins eq '64' then gain_tmp = fcals.cal_spec.scmv_b_spec64.gain_vs_freq
							if bins eq '112' then gain_tmp = fcals.cal_spec.scm_b_spec112.gain_vs_freq
			  				;nT per volt square can conversion curve
							if bins eq '36' then nT_volt = fcals.cal_spec.scmv_b_spec36.stimcoil_nt2v
							if bins eq '64' then nT_volt = fcals.cal_spec.scmv_b_spec64.stimcoil_nt2v
							if bins eq '112' then nT_volt = fcals.cal_spec.scmv_b_spec112.stimcoil_nt2v							
						endelse


						gainf[*,bb] = 1/(gain_tmp * nT_volt)
						offset[bb] = cp.adc_offset_msc[1]
						conversion[bb] = pTf
						units[bb] = pTnT+'!U2!N'
						zlims[bb,*] = [(pTnT_zscale*1d-4)^2,(pTnT_zscale*1d-2)^2]
					end
			'SCMW': begin
						gain[bb] = cp.adc_gain_msc[2]
						;*******************
						;***TEMPORARY - rbsp_efw_get_cal_params HAS A ROUGH SQ-CAN FACTOR
						;BUILT IN. I GET THE CORRECT CURVES FROM rbsp_efw_get_gain_results. SO, I NEED
						;TO MANUALLY REMOVE THE CORRECTION HERE. 
						;if probe eq 'a' then gain[bb] /= 2.22 else gain[bb] /= 2.22
						;*******************
						if probe eq 'a' then begin
							if bins eq '36' then gain_tmp = fcals.cal_spec.scmw_a_spec36.gain_vs_freq
							if bins eq '64' then gain_tmp = fcals.cal_spec.scmw_a_spec64.gain_vs_freq
							if bins eq '112' then gain_tmp = fcals.cal_spec.scmw_a_spec112.gain_vs_freq
			  				;nT per volt square can conversion curve
							if bins eq '36' then nT_volt = fcals.cal_spec.scmw_a_spec36.stimcoil_nt2v
							if bins eq '64' then nT_volt = fcals.cal_spec.scmw_a_spec64.stimcoil_nt2v
							if bins eq '112' then nT_volt = fcals.cal_spec.scmw_a_spec112.stimcoil_nt2v
						endif else begin
							if bins eq '36' then gain_tmp = fcals.cal_spec.scmw_b_spec36.gain_vs_freq
							if bins eq '64' then gain_tmp = fcals.cal_spec.scmw_b_spec64.gain_vs_freq
							if bins eq '112' then gain_tmp = fcals.cal_spec.scmw_b_spec112.gain_vs_freq
			  				;nT per volt square can conversion curve
							if bins eq '36' then nT_volt = fcals.cal_spec.scmw_b_spec36.stimcoil_nt2v
							if bins eq '64' then nT_volt = fcals.cal_spec.scmw_b_spec64.stimcoil_nt2v
							if bins eq '112' then nT_volt = fcals.cal_spec.scmw_b_spec112.stimcoil_nt2v							
						endelse
						gainf[*,bb] = 1/(gain_tmp * nT_volt)
						offset[bb] = cp.adc_offset_msc[2]
						conversion[bb] = pTf
						units[bb] = pTnT+'!U2!N'
						zlims[bb,*] = [(pTnT_zscale*1d-4)^2,(pTnT_zscale*1d-2)^2]
					end
		endcase
	endfor


		gain = gain/adc_factor  
	
		ylabel = src
		zlabel = 'RBSP'+probe+' EFW spec'
	
	
		ytitle = ylabel + '!C[Hz]'	
		ztitle = strarr(7)
		for j=0,6 do ztitle[j] = '!C' + units[j] + '/[Hz]!C!C' + zlabel
	
		for j=0,6 do begin
	
	
			get_data,tvar[j],data=data,dlimits=dlim
	
			if dlim.data_att.units eq 'ADC' then begin
	
			tst = size(data)
			
				;make sure there's data
				if tst[0] ne 0 then begin
		

					;subtract offset
					new_y = (double(data.y) - offset[j])

					;Compensate for 4 bit dynamic range compression
					new_y /= 4.^2
				
					;Divide the 0th bin by 4 (from weirdness doc)
					new_y[*,0] /= 4.

					;Compensate for using Hanning window
					new_y /= 0.375^2
				
					;Apply ADC conversion
					new_y *= adc_factor^2
				
					;Apply channel gain
					new_y *= gain[j]^2

					;Apply boom length to change from V to mV/m or from V to nT (or pT)		
					new_y *=  conversion[j]^2 
	
	
					;Apply freq dependent gain (and square can conversion of nT/v for SCM)
					gtmp = reform(gainf[*,j])^2
					for qq=0,n_elements(new_y[*,0])-1 do new_y[qq,*] = new_y[qq,*]*gtmp
	
				
			
					new_data = {x:data.x,y:new_y,v:data.v}
					dlim.data_att.units = units[j]
					dlim.data_att.channel = src[j]
					
					colors = [1]
				
					str_element, dlim, 'cdf', /delete
					str_element, dlim, 'code_id', /delete
				
					newname = tvar[j]
					store_data, newname,data=new_data,dlim=dlim
				
					options,newname,labflag=1
					options,newname,'ytitle',ytitle[j]
					options,newname,'ztitle',ztitle[j]
					ylim,newname,3,10000,1

		
				endif	
			endif else print,tvar[j], ' is already calibrated'
		endfor



	endif else begin
		print, ''
		dprint, 'Invalid datatype. Calibration aborted.'
		print, ''
		return
	endelse




	zlim,rbspx +'_efw_'+bins+'_spec0',zlims[0,0],zlims[0,1],1
	zlim,rbspx +'_efw_'+bins+'_spec1',zlims[1,0],zlims[1,1],1
	zlim,rbspx +'_efw_'+bins+'_spec2',zlims[2,0],zlims[2,1],1
	zlim,rbspx +'_efw_'+bins+'_spec3',zlims[3,0],zlims[3,1],1
	zlim,rbspx +'_efw_'+bins+'_spec4',zlims[4,0],zlims[4,1],1
	zlim,rbspx +'_efw_'+bins+'_spec5',zlims[5,0],zlims[5,1],1
	zlim,rbspx +'_efw_'+bins+'_spec6',zlims[6,0],zlims[6,1],1



	;--------------------------------------------
	;Fix the bin labels in the tplot structure
	;--------------------------------------------

	if bins eq '36' then begin
		tplot_names,'*efw_36_spec*',names=tnames
		fbin_label_36 = strtrim(fbins_36L,2) + '-' + strtrim(fbins_36H,2) + ' Hz'
		for i=0,size(tnames,/n_elements)-1 do begin
			get_data,tnames[i],data=d,limits=l,dlimits=dl
			if is_struct(d) then d.v=fbins_36c[0:35]
			store_data,tnames[i],data=d,limits=l,dlimits=dl
		endfor
	endif

	if bins eq '64' then begin
		tplot_names,'*efw_64_spec*',names=tnames
		fbin_label_64 = strtrim(fbins_64L,2) + '-' + strtrim(fbins_64H,2) + ' Hz'
		for i=0,size(tnames,/n_elements)-1 do begin
			get_data,tnames[i],data=d,limits=l,dlimits=dl
			if is_struct(d) then d.v=fbins_64c[0:63]
			store_data,tnames[i],data=d,limits=l,dlimits=dl
		endfor
	endif


	if bins eq '112' then begin
		tplot_names,'*efw_112_spec*',names=tnames
		fbin_label_112 = strtrim(fbins_112L,2) + '-' + strtrim(fbins_112H,2) + ' Hz'
		for i=0,size(tnames,/n_elements)-1 do begin
			get_data,tnames[i],data=d,limits=l,dlimits=dl
			if is_struct(d) then d.v=fbins_112c[0:111]
			store_data,tnames[i],data=d,limits=l,dlimits=dl
		endfor
	endif

end