;+
; PRO: das2dlm_load_cassini_rpws_waveform, ...
;
; Description: Cassini data from Radio and Plasma Wave Science, waveformw (PSD level 3 filed), das2 dataset: /Cassini/RPWS/...;
; Available datasets:
; 'HiRes_LoFreq_Waveform' - Collection: 100 Hz and 7.14 kHz sample rate, correlated 5-Component waveforms from the WFR
; 'HiRes_MidFreq_Waveform' - Collection: 27.8 kHz and 222 kHz sample rate waveforms from the WBR (PDS level 3 files)
; Data may be retured from diffrent datasets. The tplot variables will indicate the dataset name (e.g. '_Ex_05')
;
; Keywords:
; trange: Sets the time tange
; source (optional): String that defines dataset: 'MidFreq', 'LoFreq' (default: 'LoFreq')
; output (optional): String (Case sensetive!), or array of strings of output data from this collection
; Collection HiRes_LoFreq_Waveform:
; 'Bx' - Waveform from the tri-axial search coil magnetic antenna Bx, detect magnetic components of electromagnetic waves
; 'By' - Waveform from the By antenna used as a monopole, detect magnetic components of electromagnetic waves
; 'Bz' - Waveform from the Bz antenna used as a monopole, detect magnetic components of electromagnetic waves
; 'Ew' - Waveform from the Ew electric monopole antenna
; 'Ex' - Waveform from the Eu and Ev electric dipole antennas, aligned along the x axis of the spacecraft
; '7kHz' - Output 140 mks sample interval 7.143 kHz rolloff data, by default 10 ms sample data 100 Hz rolloff are sent
; Collection HiRes_MidFreq_Waveform (at this moment, onlye one parameter: Ew or Ex can be selected):
; 'Ew' - Waveform from the Ew electric monopole antenna
; 'Ex' (Currently, does not work) - Waveform from the Eu and Ev electric dipole antennas, aligned along the X axis of the spacecraft
; '80khz' - chage band rolloff to 80khz (default is 10khz)
; Examples:
; output='Ex'
; output=['Ex','Ew','7kHz']
; output=['Ew','80khz']
; resolution (optional): string of the resolution, e.g. '.21' (default, '')
; parameter (optional): string of optional das2 parameters
;
; Note:
; Please reffer to the das2.org catalog for additional information or examples of das2 use:
; https://das2.org/browse/uiowa/cassini/rpws
;
; CREATED BY:
; Alexander Drozdov (adrozdov@ucla.edu)
;
; $LastChangedBy: adrozdov $
; $Date: 2021-01-25 20:29:41 -0800 (Mon, 25 Jan 2021) $
; $Revision: 29621 $
; $URL: svn+ssh://thmsvn@ambrosia.ssl.berkeley.edu/repos/spdsoft/tags/spedas_6_1/projects/cassini/das2dlm_load_cassini_rpws_waveform.pro $
;-
pro das2dlm_load_cassini_rpws_waveform, trange=trange, source=source, resolution=resolution, output=output, parameter=parameter
das2dlm_cassini_init
if ~undefined(trange) && n_elements(trange) eq 2 $
then tr = timerange(trange) $
else tr = timerange()
if undefined(source) then source = 'LoFreq'
if undefined(output) then output = ''
if undefined(parameter) then parameter = ''
case strlowcase(source) of
'lofreq': begin
source = 'HiRes_LoFreq_Waveform'
t_name = 'time'
; filter output and leave only acceptable parameters
v_arr = strfilter(output,['Ex','Ew','Bx','By','Bz','7kHz'])
; We parameter is not empty, merge with the request
if parameter ne '' then parameter = strjoin([v_arr, parameter], '+') $
else parameter = strjoin(v_arr, '+')
; Remove '7khz' from v_arr, v_arr is a list of variables from now
if array_contains(v_arr, '7kHz') then v_arr = v_arr[where(v_arr ne '7kHz')]
; Check default (empty v_arrarr)
if not keyword_set(v_arr) then v_arr = ['Ex'] ; Default variable
; NOTE: v_arr is not in use anymore, since we need to determine the variable name dynamically
end
'midfreq': begin
source = 'HiRes_MidFreq_Waveform'
t_name = 'time'
v_arr = ['WBR']
; This collection uses diffrent way to process the parameters. Currently parameters Ew and Ex do not return data
output = strfilter(output,['Ex','Ew','80khz'])
; Set additional parameters
aparam = []
if array_contains(output, '80khz') then aparam = [aparam, '--80khz']
if array_contains(output, 'Ex') then aparam = [aparam, '--Ex']
if array_contains(output, 'Ew') then aparam = [aparam, '--Ew']
if parameter ne '' then parameter = strjoin([aparam, parameter], '+') $
else if keyword_set(aparam) then parameter = strjoin(aparam, '+')
; NOTE: v_arr is not in use anymore, since we need to determine the variable name dynamically
end
else: begin
dprint, dlevel = 0, 'Unknown source. Accepatable sources are: LoFreq or MidFreq'
return
end
endcase
if parameter ne '' then parameter = '¶ms=' + parameter.Replace(' ','+')
if undefined(resolution) $
then resolution = ''
if resolution ne '' $
then resolution = '&resolution=' + resolution
time_format = 'YYYY-MM-DDThh:mm:ss'
url = 'http://planet.physics.uiowa.edu/das/das2Server?server=dataset'
dataset = 'dataset=Cassini/RPWS/' + source
time1 = 'start_time=' + time_string( tr[0] , tformat=time_format)
time2 = 'end_time=' + time_string( tr[1] , tformat=time_format)
requestUrl = url + '&' + dataset + '&' + time1 + '&' + time2 + resolution + parameter
print, requestUrl
query = das2c_readhttp(requestUrl)
if undefined(query) then begin
dprint, dlevel = 0, 'Das2 returned invalid or empty query.'
return
endif
; Get dataset
;for i=0,size(v_arr,/n_elem)-1 do begin
for i=0,query.n_dsets-1 do begin
nset = i
ds = das2c_datasets(query, nset)
; Get time
das2dlm_get_ds_var, ds[0], 'time', 'center', p=pt, v=vt, m=mt, d=dt
; Exit on empty data
;if undefined(dt) then begin
; dprint, dlevel = 0, 'Dataset has no data for the selected period.'
; return
;endif
; Convert time
dt = das2dlm_time_to_unixtime(dt, vt.units)
; Manually fix the dimentions according to variable's rank
dt=dt[*]
; v_name = v_arr[i]
; We have to determine v_name dynamically, since it is not nessesary in the order of parameters
das2dlm_get_ds_var_name, ds[0], vnames=vnames, exclude=['time']
v_name = vnames[0]
; Get variables
das2dlm_get_ds_var, ds[0], v_name, 'center', p=pf, v=vf, m=mf, d=df
; Manually fix the dimentions according to variable's rank
df=df[*]
;dt = transpose(dt[0, *],[1, 0])
;df = df[*, 0]
;da = transpose(da, [1, 0])
tvarname = 'cassini_rpws_' + strlowcase(source) + '_' + ds[0].name ; + v_name ; ds[0].name also contains the dataset number
store_data, tvarname, data={x:dt, y:df}, $
dlimits={spec:0, ylog:0, zlog:0}
; Metadata
das2dlm_get_ds_meta, ds[0], meta=mds, title=das2name
str_element, DAS2, 'url', requestUrl, /add
str_element, DAS2, 'name', das2name, /add
str_element, DAS2, 'propds', mds, /add ; add data set property
das2dlm_add_metadata, DAS2, p=pt, v=vt, m=mt, add='t'
das2dlm_add_metadata, DAS2, p=pf, v=vf, m=mf, add='y'
options, /default, tvarname, 'DAS2', DAS2 ; Store metadata (this should not affect graphics)
options, /default, tvarname, 'title', DAS2.name
; Data Label
ytitle = DAS2.namey + ', ' + DAS2.unitsy
options, /default, tvarname, 'ytitle', ytitle ;
; Clear ds (some strange crashes happen if this is not included'
ds = !null
endfor
; Cleaning up
res = das2c_free(query)
end