;+
;FUNCTION: mvn_lpw_prd_lp_swp_setupparam_new
; Create initial Sweep parameter set forMAVEN/LPW sweep analysis.
;
;INPUTS:
;OUTPUTS: gives LP Sweep parameter set specific for the mvn_lpw_prd_lp analys program.
;KEYWORDS:
; win: use this keyword to see the every sweeep plot during the run
;
;EXAMPLE:
; mvn_lpw_prd_lp_swp_setupparam, PP, 'l0'
;
;CREATED BY: Michiko Morooka 10-21-14
;FILE: mvn_lpw_prd_lp_n_t_fit.pro
;VERSION: 3.0
;LAST MODIFICATION:
; 2014-10-20 M. Morooka
; 2014-11-17 M. Morooka Add err information
; Add current_l0, voltage_l0
; 2015-01-13 M. Morooka Add extra dummy fields
; 2015-04-16 M. Morooka dTe added (2.3)
; 2015-05-26 M. Morooka Change for new_fit (3.0)
; 2015-06-19 M. Morooka Minor addition (3.1) (Ni2, mi2)
;-
;
;------- this_version_mvn_lpw_prd_lp_swp_setupparam -------
function this_version_mvn_lpw_prd_lp_swp_setupparam
ver = 3.1
prd_ver= 'version mvn_lpw_prd_lp_swp_setupparam: ' + string(ver,format='(F4.1)')
return, prd_ver
end
;--------- this_version_mvn_lpw_prd_lp_swp_setupparam -----
; ------ swp_setupparam -----------------------------------
function mvn_lpw_prd_lp_swp_setupparam, prd_ver_in
;------ the version number of this routine --------------------------------------------------------
t_routine=SYSTIME(0)
prd_ver= this_version_mvn_lpw_prd_lp_swp_setupparam()
;print, '------------------------------' & print, prd_ver & print, '------------------------------'
;--------------------------------------------------------------------------------------------------
PROJECT = 'MAVEN/LPW'
rp = 6.35e-3*0.5 ; Diameter of Langmuir Probe [m]
lp = 0.4 ; Length of cylindrical Langmuir Probe [m]
Ap = !pi*6.35e-3*0.4 ; Area of probe surface [m^2]
XA = 6.35e-3*0.4 ; Effective probe surface in [m^2], (lp*rp)
RXA = 1.0 ; persentage of probe surface (0-1.0)
R_sun = 1.5 ; distance from sun in [AU]
v_len = 128 ; Max number of Voltage steps
PP = create_struct( $
'ptitle', 'Sweep Current Plots' $ ; plotting title
,'xlim',[!values.F_nan, !values.F_nan] $
,'proj', PROJECT $ ; project name
,'probe', 0 $ ; Probe number
,'rp', rp $ ; Radius of Langmuir Probe [m]
,'lp', lp $ ; Length of cylindrical Langmuir Probe [m]
,'Ap', Ap $ ; Area of probe surface [m^2]
,'XA', XA $ ; Effective probe surface in [m^2], (2*lp*rp)
,'RXA', RXA $ ; persentage of probe surface (0-1.0)
,'R_sun', R_sun $ ; distance from sun in [AU]
; +++++ Operation information +++++
,'swp_mode', -1.0 $ ; sweep mode
; +++++ Spacectaft Attitude information +++++
; +++++ Sweep Info. +++++
,'voltage_l0', make_array(v_len,1,value=!values.F_nan) $
,'current_l0', make_array(v_len,1,value=!values.F_nan) $
,'time_l0', make_array(v_len,1,value=double(!values.F_nan)) $
,'voltage', make_array(v_len,1,value=!values.F_nan) $
,'time', double(!values.F_nan) $
,'current', make_array(v_len,1,value=!values.F_nan) $
,'I_photo', make_array(v_len,1,value=!values.F_nan) $
,'I_ion', make_array(v_len,1,value=!values.F_nan) $
,'I_electron1', make_array(v_len,1,value=!values.F_nan) $
,'I_electron2', make_array(v_len,1,value=!values.F_nan) $
,'I_electron3', make_array(v_len,1,value=!values.F_nan) $
,'I_tot', make_array(v_len,1,value=!values.F_nan) $
,'I_tmp', make_array(v_len,1,value=!values.F_nan) $
,'I_ion2', make_array(v_len,1,value=!values.F_nan) $
,'I_tot2', make_array(v_len,1,value=!values.F_nan) $
,'I_off', !values.F_nan $ ; Instrumental current error (set 8.78647e-09 for ver. before 2.0)
; +++++ Plasma parameters +++++
,'Vsc', !values.F_nan $ ; Spacecraft velosity respect to plasma [km]
; +++++ positive bias side +++++
,'No_e', !values.F_nan $ ; Number of electron populations (1, 2, or 3 only)
,'Ne_tot', !values.F_nan $ ; Total electron density in [cm(-3)]
,'Ne1', !values.F_nan $ ; Electron density of comp1 in [cm(-3)]
,'Ne2', !values.F_nan $ ; Electron density of comp2 in [cm(-3)]
,'Ne3', !values.F_nan $ ; Electron density of comp3 in [cm(-3)]
,'dNe_tot', !values.F_nan $ ; Total electron density in [cm(-3)]
,'dNe1', !values.F_nan $ ; error
,'dNe2', !values.F_nan $ ; error
,'dNe3', !values.F_nan $ ; error
,'Neprox', !values.F_nan $ ; Enectron density in thin plasma region using floating potential proxy
,'dNeprox', !values.F_nan $ ; error
,'U_zero', !values.F_nan $ ; Voltage that current goes to zero in [V]
,'U0', !values.F_nan $ ; Probe Floating potential [V]
,'U1', !values.F_nan $ ; Characteristic potentials between probe-to-plasma1 [V]
,'U2', !values.F_nan $ ; Characteristic potentials between probe-to-plasma2 [V]
,'dU0', !values.F_nan $ ; error
,'dU1', !values.F_nan $ ; error
,'dU2', !values.F_nan $ ; error
,'Usc', !values.F_nan $ ; Spacecraft potential (Debye length assumed with obtained density)
,'dUsc', !values.F_nan $ ; error
,'Usc_corr', !values.F_nan $ ; Spacecraft potential corrected by attitudes
,'Te', !values.F_nan $ ; Average Electron temperature in [eV]
,'Te1', !values.F_nan $ ; Electron temperature of comp1 in [eV]
,'Te2', !values.F_nan $ ; Electron temperature of comp2 in [eV]
,'Te3', !values.F_nan $ ; Electron temperature of comp3 in [eV]
,'dTe', !values.F_nan $ ; error of main electron component in [eV]
,'dTe1', !values.F_nan $ ; error
,'dTe2', !values.F_nan $ ; error
,'dTe3', !values.F_nan $ ; error
,'fnorm', !values.F_nan $ ;
; +++++ negative bias side +++++
,'Ni', !values.F_nan $ ; Ion effective density in [cm(-3)]
,'Ni2', !values.F_nan $ ; Ion effective density in [cm(-3)]
,'Ufloat', !values.F_nan $ ; Ufloat to calcurate Ni [V] (=Usc)
,'Ufloat2', !values.F_nan $ ; Alternative Ufloat to calcurate Ni2 [V]
,'Ti', !values.F_nan $ ; Ion thermal velocity in [eV]
,'Vi', !values.F_nan $ ; Ion drift velocity in [km/s]
,'m', !values.F_nan $ ; Linearised Ion current constans
,'m2', !values.F_nan $ ; Linearised Ion current constans (squared m)
,'m_corr', !values.F_nan $ ; 'M' value corrected by attitude
,'b', !values.F_nan $ ; Linearised Ion current gradient
,'b2', !values.F_nan $ ;Linearised Ion current gradient (squared b)
,'mi', !values.F_nan $ ; Average Ion mass number
,'mi2', !values.F_nan $ ; Average Ion mass number
,'dNi', !values.F_nan $ ; error
,'dTi', !values.F_nan $ ; error
,'dVi', !values.F_nan $ ; error
,'dM', !values.F_nan $ ; error
,'dB', !values.F_nan $ ; error
,'dmi', !values.F_nan $ ; error
,'mb_norm', !values.F_nan $ ; Norm value for linear fitting
,'UV', !values.F_nan $ ; UV level
,'dUV', !values.F_nan $ ; error
; +++++ fitting info +++++
,'Ii_ind', make_array(v_len,1,value=!values.F_nan) $
,'Ie_ind', make_array(v_len,1,value=!values.F_nan) $
; +++++ numbers for calibrations
,'ifit_points', 0 $ ; number of points used for ion fitting for IOSP
,'efit_points', [0, 0, 0] $ ; number of points used for electron fitting for IOSP
,'flg', 0 $
;,'flag_info', '' $
;,'flag_source', '' $
,'fit_err', make_array(v_len,2,value=!values.F_nan) $ ; fitting err value
,'fit_err2', !values.F_nan $ ; fitting err value
,'fit_err3', !values.F_nan $ ; fitting err value
,'fit_err4', make_array(v_len,2,value=!values.F_nan) $ ; fitting err value
; +++++ add prd ver in the end
,'prd_ver', prd_ver_in + ' # ' + prd_ver $
,'fit_function_name', '' $ ; fitting routine name;
,'U_input', [!values.F_nan,!values.F_nan,!values.F_nan] $ ; Initial suggestion of U if needed
; +++++ add extra field for in case
,'dummy_arr1', make_array(v_len,1,value=!values.F_nan) $ ;---
,'dummy_arr2', make_array(v_len,1,value=!values.F_nan) $ ;---
,'dummy_arr3', make_array(v_len,1,value=!values.F_nan) $ ;---
,'dummy_arr4', make_array(v_len,1,value=!values.F_nan) $ ; Ie_ind for N1 in case of two components
,'dummy1', !values.F_nan $ ; ----- Ufloat from ion fitting (20150606)
,'dummy2', !values.F_nan $ ; ----- unfixed RXA for ion (20150606)
,'dummy3', !values.F_nan $ ; ----- sq_linear fit param A
,'Beta', !values.F_nan $
,'peaks', make_array(3,1,value = !values.F_nan) $
,'WGaussian', make_array(3,1,value = !values.F_nan) $
,'HGaussian', make_array(3,1,value = !values.F_nan) $
,'U_Imin', !values.F_nan $
,'U_didvmax', !values.F_nan $
)
return, PP
end
; ------------------------------------ swp_setupparam -----