;Procedure: RBSP_LOAD_EFW_XSPEC
;
;Purpose: Loads RBSP EFW XSPEC data
;
;keywords:
; probe = Probe name. The default is 'all', i.e., load all available probes.
; This can be an array of strings, e.g., ['a', 'b'] or a
; single string delimited by spaces, e.g., 'a b'
; varformat=string
; TRANGE= (Optional) Time range of interest (2 element array), if
; this is not set, the default is to prompt the user. Note
; that if the input time range is not a full day, a full
; day's data is loaded
; VARNAMES: names of variables to load from cdf: default is all.
; /GET_SUPPORT_DATA: load support_data variables as well as data variables
; into tplot variables. (NOT IMPLEMENTED YET)
; /DOWNLOADONLY: download file but don't read it. (NOT IMPLEMENTED YET)
; /INTEGRATION: If set, load data from integration.
; /valid_names, if set, then this routine will return the valid probe, datatype
; and/or level options in named variables supplied as
; arguments to the corresponding keywords.
; files named varible for output of pathnames of local files.
; /VERBOSE set to output some useful info
; type: set to 'calibrated' to automatically convert data into physical units
;Example:
; rbsp_load_efw_xspec,/get_suppport_data,probe=['a', 'b']
;Notes:
; 1. Written by Peter Schroeder, February 2012
;
; $LastChangedBy: peters $
; $LastChangedDate: 2012-03-16 10:52:30 -0700 (Fri, 16 Mar 2012) $
; $LastChangedRevision: 10145 $
; $URL: svn+ssh://thmsvn@ambrosia.ssl.berkeley.edu/repos/ssl_general/tags/tdas_7_00/missions/rbsp/efw/rbsp_load_efw_xspec.pro $
;-
pro rbsp_load_efw_xspec,probe=probe, datatype=datatype, trange=trange, $
level=level, verbose=verbose, downloadonly=downloadonly, $
cdf_data=cdf_data,get_support_data=get_support_data, $
tplotnames=tns, make_multi_tplotvar=make_multi_tplotvar, $
varformat=varformat, valid_names = valid_names, files=files,$
type=type, integration=integration
rbsp_efw_init
dprint,verbose=verbose,dlevel=4,'$Id: rbsp_load_efw_xspec.pro 10145 2012-03-16 17:52:30Z peters $'
if(keyword_set(probe)) then $
p_var = probe
vb = keyword_set(verbose) ? verbose : 0
vb = vb > !rbsp_efw.verbose
vprobes = ['a','b']
vlevels = ['l1','l2']
vdatatypes=['xspec']
if ~keyword_set(type) then begin
type = 'raw'
endif
if keyword_set(valid_names) then begin
probe = vprobes
level = vlevels
datatype = vdatatypes
return
endif
if not keyword_set(p_var) then p_var='*'
p_var = strfilter(vprobes, p_var ,delimiter=' ',/string)
if not keyword_set(datatype) then datatype='*'
datatype = strfilter(vdatatypes, datatype ,delimiter=' ',/string)
if not keyword_set(level) then level='*'
level = strfilter(vdatatypes, level ,delimiter=' ',/string)
addmaster=0
probe_colors = ['m','b']
for s=0,n_elements(p_var)-1 do begin
rbspx = 'rbsp'+ p_var[s]
if keyword_set(integration) then rbsppref = rbspx + '/l1_int' $
else rbsppref = rbspx + '/l1'
format = rbsppref + '/xspec/YYYY/'+rbspx+'_l1_xspec_32_YYYYMMDD_v00.cdf'
relpathnames = file_dailynames(file_format=format,trange=trange,addmaster=addmaster)
; if vb ge 4 then printdat,/pgmtrace,relpathnames
dprint,dlevel=3,verbose=verbose,relpathnames,/phelp
files = file_retrieve(relpathnames, _extra=!rbsp_efw)
if keyword_set(!rbsp_efw.downloadonly) or keyword_set(downloadonly) then continue
; suf='_raw'
suf=''
; midfix='_hsk_beb_analog_'
prefix=rbspx+'_efw_xspec_32_'
; if keyword_set(get_support_data) then $
cdf2tplot,file=files,varformat=varformat,all=0,prefix=prefix,suffix=suf,verbose=vb, $
tplotnames=tns,/convert_int1_to_int2,get_support_data=get_support_data ; load data into tplot variables
if is_string(tns) then begin
xspecname = rbspx+'_efw_xspec_32_xspec_32'
get_data,xspecname+'_src1',data=src1data
get_data,xspecname+'_src2',data=src2data
get_data,xspecname+'_rc',data=rcdata
get_data,xspecname+'_ic',data=icdata
for i = 0,3 do begin
specnewname = rbspx+'_efw_32_xspec'+strcompress(i,/rem)
store_data,specnewname+'_src1',data={x: src1data.x, $
y: reform(src1data.y[*,*,i]), v: src1data.v2}, $
dlimits = {spec: 1}
store_data,specnewname+'_src2',data={x: src2data.x, $
y: reform(src2data.y[*,*,i]), v: src2data.v2}, $
dlimits = {spec: 1}
store_data,specnewname+'_rc',data={x: rcdata.x, $
y: reform(rcdata.y[*,*,i]), v: rcdata.v2}, $
dlimits = {spec: 1}
store_data,specnewname+'_ic',data={x: icdata.x, $
y: reform(icdata.y[*,*,i]), v: icdata.v2}, $
dlimits = {spec: 1}
endfor
; pn = byte(p_var[s]) - byte('a')
; options, /def, tns, colors = probe_colors[pn]
dprint, dlevel = 5, verbose = verbose, 'Setting options...'
options, /def, tns, code_id = '$Id: rbsp_load_efw_xspec.pro 10145 2012-03-16 17:52:30Z peters $'
dprint, dwait = 5., verbose = verbose, 'Flushing output'
dprint, dlevel = 4, verbose = verbose, 'XSPEC 32 data Loaded for probe: '+p_var[s]
endif else begin
dprint, dlevel = 0, verbose = verbose, 'No EFW XSPEC 32 data loaded...'+' Probe: '+p_var[s]
; dprint, dlevel = 0, verbose = verbose, 'Try using get_support_data keyword'
endelse
format = rbsppref + '/xspec/YYYY/'+rbspx+'_l1_xspec_64_YYYYMMDD_v00.cdf'
relpathnames = file_dailynames(file_format=format,trange=trange,addmaster=addmaster)
; if vb ge 4 then printdat,/pgmtrace,relpathnames
dprint,dlevel=3,verbose=verbose,relpathnames,/phelp
files = file_retrieve(relpathnames, _extra=!rbsp_efw)
if keyword_set(!rbsp_efw.downloadonly) or keyword_set(downloadonly) then continue
; suf='_raw'
suf=''
; midfix='_hsk_beb_analog_'
prefix=rbspx+'_efw_xspec_64_'
; if keyword_set(get_support_data) then $
cdf2tplot,file=files,varformat=varformat,all=0,prefix=prefix,suffix=suf,verbose=vb, $
tplotnames=tns,/convert_int1_to_int2,get_support_data=get_support_data ; load data into tplot variables
if is_string(tns) then begin
xspecname = rbspx+'_efw_xspec_64_xspec_64'
get_data,xspecname+'_src1',data=src1data
get_data,xspecname+'_src2',data=src2data
get_data,xspecname+'_rc',data=rcdata
get_data,xspecname+'_ic',data=icdata
for i = 0,3 do begin
specnewname = rbspx+'_efw_64_xspec'+strcompress(i,/rem)
store_data,specnewname+'_src1',data={x: src1data.x, $
y: reform(src1data.y[*,*,i]), v: src1data.v2}, $
dlimits = {spec: 1}
store_data,specnewname+'_src2',data={x: src2data.x, $
y: reform(src2data.y[*,*,i]), v: src2data.v2}, $
dlimits = {spec: 1}
store_data,specnewname+'_rc',data={x: rcdata.x, $
y: reform(rcdata.y[*,*,i]), v: rcdata.v2}, $
dlimits = {spec: 1}
store_data,specnewname+'_ic',data={x: icdata.x, $
y: reform(icdata.y[*,*,i]), v: icdata.v2}, $
dlimits = {spec: 1}
endfor
; pn = byte(p_var[s]) - byte('a')
; options, /def, tns, colors = probe_colors[pn]
dprint, dlevel = 5, verbose = verbose, 'Setting options...'
options, /def, tns, code_id = '$Id: rbsp_load_efw_xspec.pro 10145 2012-03-16 17:52:30Z peters $'
dprint, dwait = 5., verbose = verbose, 'Flushing output'
dprint, dlevel = 4, verbose = verbose, 'XSPEC 64 data Loaded for probe: '+p_var[s]
endif else begin
dprint, dlevel = 0, verbose = verbose, 'No EFW XSPEC 64 data loaded...'+' Probe: '+p_var[s]
; dprint, dlevel = 0, verbose = verbose, 'Try using get_support_data keyword'
endelse
format = rbsppref + '/xspec/YYYY/'+rbspx+'_l1_xspec_112_YYYYMMDD_v00.cdf'
relpathnames = file_dailynames(file_format=format,trange=trange,addmaster=addmaster)
; if vb ge 4 then printdat,/pgmtrace,relpathnames
dprint,dlevel=3,verbose=verbose,relpathnames,/phelp
files = file_retrieve(relpathnames, _extra=!rbsp_efw)
if keyword_set(!rbsp_efw.downloadonly) or keyword_set(downloadonly) then continue
; suf='_raw'
suf=''
; midfix='_hsk_beb_analog_'
prefix=rbspx+'_efw_xspec_112_'
; if keyword_set(get_support_data) then $
cdf2tplot,file=files,varformat=varformat,all=0,prefix=prefix,suffix=suf,verbose=vb, $
tplotnames=tns,/convert_int1_to_int2,get_support_data=get_support_data ; load data into tplot variables
if is_string(tns) then begin
xspecname = rbspx+'_efw_xspec_112_xspec_112'
get_data,xspecname+'_src1',data=src1data
get_data,xspecname+'_src2',data=src2data
get_data,xspecname+'_rc',data=rcdata
get_data,xspecname+'_ic',data=icdata
for i = 0,3 do begin
xspecnewname = rbspx+'_efw_112_xspec'+strcompress(i,/rem)
store_data,xspecnewname+'_src1',data={x: src1data.x, $
y: reform(src1data.y[*,*,i]), v: src1data.v2}, $
dlimits = {spec: 1}
store_data,xspecnewname+'_src2',data={x: src2data.x, $
y: reform(src2data.y[*,*,i]), v: src2data.v2}, $
dlimits = {spec: 1}
store_data,xspecnewname+'_rc',data={x: rcdata.x, $
y: reform(rcdata.y[*,*,i]), v: rcdata.v2}, $
dlimits = {spec: 1}
store_data,xspecnewname+'_ic',data={x: icdata.x, $
y: reform(icdata.y[*,*,i]), v: icdata.v2}, $
dlimits = {spec: 1}
endfor
pn = byte(p_var[s]) - byte('a')
options, /def, tns, colors = probe_colors[pn]
dprint, dlevel = 5, verbose = verbose, 'Setting options...'
options, /def, tns, code_id = '$Id: rbsp_load_efw_xspec.pro 10145 2012-03-16 17:52:30Z peters $'
dprint, dwait = 5., verbose = verbose, 'Flushing output'
dprint, dlevel = 4, verbose = verbose, 'XSPEC 112 data Loaded for probe: '+p_var[s]
endif else begin
dprint, dlevel = 0, verbose = verbose, 'No EFW XSPEC 112 data loaded...'+' Probe: '+p_var[s]
; dprint, dlevel = 0, verbose = verbose, 'Try using get_support_data keyword'
endelse
endfor
end