This page was created by the IDL library routine
mk_html_help2.
Last modified: Sat Apr 5 18:16:45 2025.
PRO: REE_SET_GSM, sc
PURPOSE: LOADS STATE AD SETS UP PLOTS FOR GSM
INPUT:
SC - REQUIRED. STRING. Spacecraft
EXAMPLE: ree_set_gsm, 'a'
OUTPUT: GSM Coordinates are laoded into tplot.
WARNING: BE SURE TO SET TIMESPAN FIRST!
INITIAL VERSION: REE 08-10-31
MODIFICATION HISTORY:
LASP, CU
(See projects/themis/spacecraft/fields/LASP/ree_set_gsm.pro)
PRO: THM_EFI_CLEAN_EFP
PURPOSE:
FIXES AXIAL, REMOVES SPIKES, SPINTONE, SC POTENTIAL, AND OTHER EFP ERRORS.
This routine is designed for general use. It will take EFP data in DSL as
input, clean up the data, and transform the data into FAC and GSM, and
return the cleaned EFP data in DSL, FAC, and GSM. If the data is gathered in
the dayside, set keyword /SUBSOLAR to use keyword defaults specifically for
dayside data.
The routine can work on one-day long data, but if only a single particle
burst data is interesting, use the keyword TRANGE to specify the time range
of that burst can greatly reduce run time.
CAUTION: The input keyword ENAME is optional. When it is not specified,
ENAME is defaulted to 'thx_efp', where 'x' is probe name and should be
one of these: 'a', 'b', 'c', 'd', and 'e'. In this case, if the existed
'thx_efp' is not in DSL, the output will be wrong. So, make sure the
existed 'thx_efp' is in DSL, or specify ENAME to another tplot variable
which is EFP data in DSL.
NOTE: 1) Spike removal does not always work, so don't be surprised if there
still are spikes left on the data. Hopefully the spike issue can be
resolved in the future.
2) Merge can exhibit pathological behavior if B is near the spin plane.
spike removal does not always work.
3) Make sure thx_state_spinper is available! Best to set up for gsm.
Make sure timespan is set!!!!
EXAMPLES:
For typical use, see thm_crib_cleanefp.pro and thm_crib_cleanefi.pro which
are typically located under TDAS_DIR/idl/themis/examples.
INPUT:
KEYWORDS (general):
probe NEEDED: Program does only one sc at a time!
Ename OPTIONAL: Valid TPLOT name for E, DEFAULT = 'thX_efp' (will
fetch)
Vname OPTIONAL: Valid TPLOT name for E, DEFAULT = 'thX_vap' (will
fetch)
Bdslname OPTIONAL, Valid TPLOT name for B, DEFAULT = 'thX_fgh_dsl'
(will fetch)
trange OPTIONAL. Time range.
USE TRANGE TO ISOLATE SINGLE PARTICLE BURST AND GREATLY
REDUCE RUN TIME!!!!
subsolar OPTIONAL. If set, defaults for subsolar region are used.
DEFAULT = TAIL
talk OPTIONAL. Plots diagnostics. Can be fun. Slow. DEFAULT = 0
duration_threshold: OPTIONAL. If the duration of a burst is less than
duration_threshold, it is skipped. DEFAULT = 15.
In units of seconds.
KEYWORDS (for output tplot variables):
Edslname OPTIONAL: Name for the cleaned efp in DSL.
DEFAULT = Ename + '_clean_dsl'
Egsmname OPTIONAL: Name for the cleaned efp in DSL.
DEFAULT = Ename + '_clean_gsm'
Efacname OPTIONAL: Name for the cleaned efp in DSL.
DEFAULT = Ename + '_clean_fac'
KEYWORDS (Remove_SpinTone):
SpinRemove If entered as 0, suppresses spintone removal. DEFAULT = 4
(Spintone removed). The following are the behaviors for
different values of SpinRemove.
1: Only spintone, 2: Spintone and 2nd harm. 3: Spintone and
4th harm. 4: Spintone, 2nd, and 4th harm.
SpinNsmooth Activates median smoothing. # of half periods. Must be odd.
DEFAULT = 19 for subsolar and 39 otherwise.
SpinPoly Activates polyfit smoothing. Can be unstable if >9. DEFAULT
= 0
KEYWORDS (Remove_Potential):
VscRemove If entered as 0, suppresses potential removal. DEFAULT = 1
(Potential removed.)
VscPole Frequency to median smooth Vsc before applying to Ez.
DEFAULT = 5.0 Hz
Vpoly Actives polyfit of Ez to Vsc. DEFAULT = 2 (5 /Subsolar)
Use_Electrons (OBSOLETE) Kept for backward compatibility.
TeName (OBSOLETE) Kept for backward compatibility.
KEYWORDS (Remove_Spikes):
SpikeRemove If entered as 0, suppresses spike removal. DEFAULT = 1
(Spikes are removed.)
SpikeNwin Number of points in spike search window. DFLT = 16
SpikeSig Sigma of spikes. DFLT = 5 (Sometimes Sig = 6 works better)
SpikeSigmin Minimun of sigma. DFLT = 0.01 mV/m. (Sometimes Sigmin = 0
works better)
SpikeNfit Number of points in the fit window. DFLT = 16
SpikeFit If set, will do a Gaussian fit to Spikes. DFLT = 0
KEYWORDS (Remove_Spin_Epoch):
EpochRemove If set, removes all spin-epoch signals. DEFAULT = 0 (Not
used.) Only useful for short stretches or if plasma
condtitions are constant.
KEYWORDS (Fix_Axial)
FixAxial If entered as 0, suppresses axial fix. DEFAULT=1 (Axial is
fixed.)
AxPoly Forces polynomial fit of Ez-Eder of order AxPoly. DEFAULT=2
(9 /SubSolar)
Merge Merges Ederived with Eaxial. DEFAULT = 0 No merging
(1 /Subsolar)
1 - Abrupt merging. If Eder is avaliable Eder is used
(0 - Fmerge). Else Ez.
2 - Soft merging. Bz/|B| is considered when merging.
WARNING: MERGE=1 OR 2 CAN EXHIBIT PATHOLOGICAL BEHAVIOR IF
MAGNETIC FIELD IS NEAR THE SPIN PLANE. PLEASE DOUBLE
CHECK RESULTS IF USED.
Fmerge Frequency of crossover for merging Ederived with Eaxial.
DEFAULT = 5.0 Hz (/Subsolar).
MergeRatio Mimimum value of Bz/|B| to start merge (USED IF MERGE=2).
DEFAULT = 0.05
MinBz Needed to avoid divide by zero in Ederive. DEFAULT = 0.01 nT
MinRatio Mimimum value of Bz/|B| to calculate Ederive. DEFAULT = 0.1
(0.05 /Subsolar)
KEYWORDS (B_Smooth)
BspinPoly Activates polyfit smoothing. Can be unstable if >9. DEFAULT
= 2. -1 suppresses.
Bsmooth Activates B smoothing for fac rotation. DEFAULT = 11. 1 or
less suppresses.
HISTORY:
2013-06-20: JBT. Added deconvolution of sheath response.
2011-05-21: JBT. Added keyword DURATION_THRESHOLD.
2011-05-20: JBT. Fixed a bug that happens when the spin-tone removing fails.
2010-09-13: JBT. The default of FixAxial was changed to 0 (not to fix
axial E-field component).
Updated documents.
2010-04-08: JBT. Fixed some bugs.
2009-06-04: REE/JBT. Second Release.
2009-06-04: Jianbao Tao (JBT). The metadata (dlim.data_att) was improved.
2009-05-05: REE. First Release.
(See projects/themis/spacecraft/fields/LASP/thm_efi_clean_efp.pro)
PRO: THM_EFI_CLEAN_EFW
PURPOSE:
FIXES HF, REMOVES SPIKES, AND LOW FREQUENCY.
This routine is designed for general use. It will take EFW data in
DSL as input, clean up the data, and transform the data into FAC, and return
the cleaned EFP data in DSL/FAC.
The routine can work on one-day long data, but if only a single particle
burst data is interesting, use the keyword TRANGE to specify the time range
of that burst can greatly reduce run time.
NOTE: 1) Spike removal does not always work, so don't be surprised if there
still are spikes left on the data. Hopefully the spike issue can be
resolved in the future.
2) Make sure thx_state_spinper is available! Best to set up for gsm.
Make sure timespan is set!!!!
EXAMPLES:
For typical use, see thm_crib_cleanefw.pro and thm_crib_cleanefi.pro which
are typically located under TDAS_DIR/idl/themis/examples.
INPUT:
KEYWORDS (general):
probe NEEDED: Program does only one sc at a time!
Ename OPTIONAL: Valid TPLOT name for E WAVE, DEFAULT = 'thX_efw'
(will fetch)
Bdslname OPTIONAL, Valid TPLOT name for B, DEFAULT = 'thX_fgh_dsl'
(will fetch)
trange OPTIONAL. Time range.
USE TRANGE TO ISOLATE SINGLE PARTICLE BURST AND GREATLY
REDUCE RUN TIME!!!!
talk OPTIONAL. Plots diagnostics. Can be fun. Slow. DEFAULT = 0
EfpName OPTIONAL: Valid TPLOT name for EFP, DEFAULT = 'thX_efp'
(will fetch)
EFP is used to find spike positions.
status OPTIONAL: A named variable to return the exiting status. If
exit successfully, status = 0, otherwise status = 1
EFWHEDACNAME OPTIONAL: Valid TPLOT name for the AC-coupled E-field header info
DEFAULT = 'thX_efw_hed_ac'
KEYWORDS (Filter):
FPole (OBSOLETE) Kept for backward compatibility.
KEYWORDS (for output tplot variables):
Edslname OPTIONAL: Name for the cleaned efp in DSL.
DEFAULT = Ename + '_clean_dsl'
Efacname OPTIONAL: Name for the cleaned efp in DSL.
DEFAULT = Ename + '_clean_fac'
KEYWORDS (Remove_Spikes):
SpikeRemove If entered as 0, suppresses spike removal.
DEFAULT = 1 (Spikes are removed.)
SpikeNwin Number of points in spike search window. DFLT = 16
SpikeSig Sigma of spikes. DFLT = 5
(Sometimes Sig = 6 works better)
SpikeSigmin Minimun of sigma. DFLT = 0.01 mV/m.
(Sometimes Sigmin = 0 works better)
SpikeNfit Number of points in the fit window. DFLT = 16
SpikeFit If set, will do a Gaussian fit to Spikes. DFLT = 0
HISTORY:
2009-05-06: REE. First Release.
2009-06-04: Jianbao Tao (JBT), CU/LASP.
The metadata (dlim.data_att) is modified properly.
2009-06-17: JBT, CU/LASP.
Keyword STATUS added.
2010-04-08: JBT, CU/LASP.
The high-pass filter is changed from an FFT filter to an
FIR filter. Some modification is made to be compatible with
two different sample rates 8192 Hz and 16384 Hz.
2010-09-13: JBT, CU/LASP.
Cleaned up documents.
Added keyword dt in calling thm_efi_fix_freq_and_phase.
2012-07-26: JBT, SSL, UC Berkeley.
1. Incorporated C. Cully's calibration for
1) Plasma-probe coupling
2) Anti-aliasing Bessel filter
3) ADC response
4) DFB response
2012-08-30: CMC, University of Calgary.
Fixed problems with AC-coupled data sampled at < 16 ksps
VERSION:
$LastChangedBy $
$LastChangedDate $
$LastChangedRevision $
$URL $
(See projects/themis/spacecraft/fields/LASP/thm_efi_clean_efw.pro)
PRO: THM_EFI_DERIVE_EZ, Edsl, Bdsl, new_name=new_name
PURPOSE: DERIVES EZ ASSUMING EdotB = 0
INPUT:
Edsl - REQUIRED. STRING Efield data name (tplot).
Bdsl - REQUIRED. STRING Bfield data (tplot)
CALLING: thm_efi_derive_Ez, Edsl, Bdsl
OUTPUT: Tplot
INITIAL VERSION: REE 08-11-04
MODIFICATION HISTORY:
LASP, CU
(See projects/themis/spacecraft/fields/LASP/thm_efi_derive_ez.pro)
PRO: THM_EFI_ExB, Ename, Bname, Sname=Sname, Vname=Vname, $
EdotB=EdotB, Btot=Btot
PURPOSE: Calculate Poynting Flux and Flow Velocity
INPUT:
Eanme - REQUIRED. STRING Efield data name (tplot).
Bname - REQUIRED. STRING Bfield data (tplot)
Sname - OPTIONAL. Poynting Flux name (tplot)
Vname - OPTIONAL. Velocity Flow name (tplot)
CALLING: thm_efi_exb, Ename, Bname, Sname
OUTPUT: Poynting Flux is Calculated
INITIAL VERSION: REE 08-10-31
LASP, CU
MODIFICATION HISTORY:
2010-02-21: Added a check of the coordinate systems of Ename and Bname.
-JBT, CU/LASP.
(See projects/themis/spacecraft/fields/LASP/thm_efi_exb.pro)
FUNCTION: THM_EFI_FIX_FREQ_AND_PAHSE, E, freq=freq, ghf=ghf, dt=dt, axial=axial
PURPOSE: A routine which integrates data to fix the gain and remove the
phase shift due to resistive to capacitive crossover of the
THEMIS preamps. ONLY USEFUL FOR BURST DATA.
INPUT:
data - REQUIRED. Data to be fixed. DEFAULT = SPIN PLANE
KEYWORDS:
axial - OPTIONAL. Changes freq and gh to AXIAL default values.
freq - OPTIONAL. Crossover frequency in HZ. DEFAULT = SPIN PLANE
Careful! Derived from Rsheath(Csheath+Cin)
ghf - OPTIONAL. Gain at high frequency. DEFAULT = SPIN PLANE
Csheath/(Csheath+Cin)
dt - OPTIONAL. Time between samples (s). DEFAULT = 1/8192
CALLING: Esp = thm_efi_fix_freq_and_phase(Esp) or
Eax = thm_efi_fix_freq_and_phase(Ex, /ax)
OUTPUT: Data is integrated to remove gain/phase error from preamp RC
crossover.
NOTE ON DEFAULTS: Spin Plane: The front-end network, R=100k, C = 10pF is
ignored. Instead, I fudge Csh = 14 pF and Cin = 11.5 pF to realize a close
approximation of measured gain/phase (see Bonnell et al. paper).
Axial: Same system. I use Csh = 4 pf, Cin = 13 pF to approximate.
One could improve the calculation with 2 poles, but the location of Cin
(before or after the input network) would need to be questioned.
Rsh is assumed to be 20 MOhm
BEHAVIOR:
(1) USE /AX for AXIAL SIGNALS!!!
(2) DATA AT BEGINNING OF AN ARRAY MAY NOT BE CORRECTED.
Use large arrays if possible. The program relaxes by one
e-fold every 3 ms.
(3) Electric field is improved for all plasma conditions. However,
PHASE/GAIN CORRECTION MAY NOT BE ENOUGH FOR ALL PLASMA CONDITIONS.
(4) PROBES MUST BE IN SUNLIGHT!
(5) ONLY USEFUL FOR BURST DATA!
(6) BE SURE TO SET DT = 1/16384 FOR AC BURST!!!
INITIAL VERSION: REE 08-08-26
MODIFICATION HISTORY:
LASP, CU
(See projects/themis/spacecraft/fields/LASP/thm_efi_fix_freq_and_phase.pro)
FUNCTION: THM_EFI_GET_POTENTIAL, Vname
NOT FOR GENERAL USE. CALLED BY THM_EFI...
ONLY FOR ISOLATED PARTICLE OR WAVE BURSTS
NOT FOR ENTIRE ORBIT.
PURPOSE:
Remove SC potential from axial signal.
INPUT:
Vname TPLOT name of voltages.
KEYWORDS:
HISTORY:
2009-03-30: REE.
(See projects/themis/spacecraft/fields/LASP/thm_efi_get_potential.pro)
PROCEEDURE: THM_EFI_SIN_FIT, e, t, es=es, ec=ec, per=per
Called by THM_EFI_REMOVE_OFFSET_AND_SPIN, name, ask=ask
PURPOSE: DOES A QUICK SIN AND COS FIT
INPUT:
E - REQUIRED (Electric Field)
T - REQUIRED
WARNING: USE SHORT SEGMENTS
OUTPUT:
AMP - AMPLITUDE
PHS - PHASE
Per - PERIOD (Sypply and save a lot of time!)
INITIAL VERSION: REE 08-10-31
University of Colorado
MODIFICATION HISTORY:
(See projects/themis/spacecraft/fields/LASP/thm_efi_sin_fit.pro)
NAME:
THM_JBT_GET_BTRANGE (FUNCTION)
PURPOSE:
This routine is to find the starting time and ending time of each
continuous section of a give tplot variable which essentially is specified
by a tplot name such as 'tha_efp'.
If the routine exit unsuccessfully, it will return -1. Otherwise, it will
return a 2D array as [2,number_of_total_bursts] which stores the
the starting time and the ending time of each continuous section.
CALLING SEQUENCE:
btrange = thm_jbt_get_btrange(tvar, nbursts=nbursts, tind=tind)
ARGUMENTS:
tvar: (INPUT, REQUIRED). The name of a tplot variable.
KEYWORDS:
nbursts: (OUTPUT, OPTIONAL) A named variable to return the number of
sections.
tind: (OUTPUT, OPTIONAL) A named variable to return a 2D array of the index
of starting and ending time points with structure
[[starting],[ending]]
EXAMPLES:
tvar = 'tha_efp'
btrange = thm_jbt_get_btrange(tvar, nb = nb, tind = tind)
HISTORY:
2009-05-04, written by Jianbao Tao, in CU/LASP.
2010-04-08: Updated the help information. JBT, CU/LASP.
(See projects/themis/spacecraft/fields/LASP/thm_jbt_get_btrange.pro)
FUNCTION: THM_EFI_DERIVE_EZ, E, B, minBz=minBz, minRat=minRat, ratio=ratio
PURPOSE: DERIVES EZ ASSUMING EdotB = 0
INPUT:
E - REQUIRED. Data structure in DSL
B - REQUIRED. Data structure in DSL
KEYWORDS:
minBz - OPTIONAL. Keeps from divide by zero. DEFAULT = 0.01 nT
minRat - OPTIONAL. Minimum value of Bz/|B|. DEFAULT = 0.1
ratio - OUTPUT. = Bz/|B|
CALLING: thm_efi_derive_Ez, Edsl, Bdsl
OUTPUT: Tplot
INITIAL VERSION: REE 08-11-04
MODIFICATION HISTORY:
REE 09_05_01
LASP, CU
(See projects/themis/spacecraft/fields/LASP/thm_lsp_derive_ez.pro)
FUNCTION: THM_LSP_FILTER, x, dt, flow, fhigh, db=db
PURPOSE: Filters the data COMP*.
INPUT:
x - NEEDED. Just the data. Not a structure or name.
dt - NEEDED. Time between points.
freq - NEEDED. Pole of filter. If band_pass: [f1,f2] f1<f2
- If low-pass: [0,f]. If high-pass [f,0]
KEYWORDS:
db - OPTIONAL. If convol option is taken, default = 120.
CREATED BY: REE, 97-03-17 - modified 97-10-03 REE added buf_dt
FILE: fa_fields_filter.pro
VERSION: 0.0
LAST MODIFICATION:
REE 08-11-02, rewrite for THEMIS - no call external so much slower.
REE 09-04-30, made it a simple function.
REE 09-05-05, Added Gaussian option. Changed passband to match digital_fliter
(See projects/themis/spacecraft/fields/LASP/thm_lsp_filter.pro)
NAME:
THM_LSP_FILTER_HIGHPASS (FUNCTION)
PURPOSE:
CAUTION: THIS ROUTINE ONLY WORKS FOR DATA SAMPLED AT EITHER 8192 HZ OR
16384 HZ.
In general, this routine works as a high-pass filter on the input data.
The data to be filtered must be a continuously sampled 1D array. The basic
process of this routine is the following. First, the input data will be
low-pass filtered with two Finite Impulse Response filters and 2:1
decimations. The details of this procedure is described in the paper
Cully, Space Sci Rev, 2008, V141, 343-355.
Second, the output from the first step is extracted from the input data to
generate the final output. This whole process is equivalent to a high-pass
filter.
If the routine exit unsuccessfully, -1 will be retured.
CALLING SEQUENCE:
output = thm_lsp_filter_highpass(datain, dt, [keywords])
ARGUMENTS:
datain: (INPUT, REQUIRED) The data to be filtered. It must be a 1D array.
dt: (INPUT, REQUIRED) The sample interval of DATAIN in seconds;
1/dt = sample_rate.
KEYWORDS:
freqlow: (OUTPUT, OPTIONAL) The approximate lower bound of the frequency
range of the output of the routine.
HISTORY:
2009-05-03: Created by Jianbao Tao, CU/LASP.
2009-05-17: Fixed the edge effect due to sectioning original data.
Jianbao Tao, CU/LASP
2010-01-22: Fixed the issue of dealing with a too-short input array.
Jianbao Tao, CU/LASP
VERSION:
$LastChangedBy: jianbao_tao $
$LastChangedDate: 2012-07-27 12:29:50 -0700 (Fri, 27 Jul 2012) $
$LastChangedRevision: 10753 $
$URL: svn+ssh://thmsvn@ambrosia.ssl.berkeley.edu/repos/spdsoft/trunk/projects/themis/spacecraft/fields/LASP/thm_lsp_filter_highpass.pro $
(See projects/themis/spacecraft/fields/LASP/thm_lsp_filter_highpass.pro)
PRO: THM_LSP_FIND_BURST
NOT FOR GENERAL USE.
PURPOSE:
Isolate an indivitual burst for analysis.
INPUT:
Data -NEEDED. A data structure
tfind -OPTIONAL. Time of interest.
KEYWORD:
HISTORY:
2009-05-30: REE. Broke out to run with wave burst.
(See projects/themis/spacecraft/fields/LASP/thm_lsp_find_burst.pro)
FUNCTION: THM_LSP_FIND_SPIKES
NOT FOR GENERAL USE. CALLED BY THM_EFI_REMOVE_SPIKES
ONLY FOR ISOLATED PARTICLE OR WAVE BURSTS
NOT FOR ENTIRE ORBIT.
PURPOSE:
Remove the non-physical spiky signals in the efw data.
INPUT:
t -NEEDED. Time array
xx -NEEDED. EFP data, DO NOT GIVE EFW DATA. Won't work.
per -NEEDED. Spin period.
KEYWORD:
nwin -OPTIONAL. Number of points in spike search window. DFLT = 16
spikesig -OPTIONAL. Sigma of spikes. DFLT = 5
sigmin -OPTIONAL. Minimun of sigma. DFLT = 0 mV/m.
OUTPUT:
tpks RETURN VALUE. Time of peaks (mod spin per) + t0
Amp Estimated amplitude of spikes. Helps remove_spikes.
HISTORY:
2009-05-30: REE. Broke out to run with wave burst.
(See projects/themis/spacecraft/fields/LASP/thm_lsp_find_spikes.pro)
NAME:
THM_LSP_GET_SPEC (PROCEDURE)
PURPOSE:
Get the power spectrum density (PSD) of a given tplot variable.
CALLING SEQUENCE:
thm_lsp_get_spec, tvar, units = units, prefix = prefix, fftlen = fftlen, $
yrange = yrange
ARGUMENTS:
tvar: (INPUT, REQUIRED) The name of a tplot variable to calculate the
PSD from.
KEYWORDS:
units: (INPUT, OPTIONAL) A string of the units of the data in the tvar. By
default, it is obtained from dlim.data_att.units, or 'unknown' if
dlim.data_att.units is not available.
prefix: (INPUT, OPTIONAL) The prefix for tplot variables that contain the
resulting PSD. By default, prefix = tvar.
fftlen: (INPUT, OPTIONAL) The number of data points to FFT. By default,
fftlen = 512.
yrange: (INPUT, OPTIONAL) The yrange for the spectra. By default, yrange =
[1, Nyquist] in units of Hz.
ufactor: (INPUT, OPTIONAL) A factor to convert units from one to another.
instr_resp: (INPUT, OPTIONAL) The instrument response as a function of
frequency. The number of elements of instr_resp should be the same
as (fftlen/2 + 1). By default, instr_resp = 1.
/checkenergy: (INPUT, OPTIONAL) If set, the ratios of time-domain energy to
freq-domain energy are stored into tplot variables. By default, no
such output is generated.
/phase: If set, phase will be calculated and saved into tplots.
/ft : If set, the FFT of the input data will be saved.
EXAMPLES:
See thm_crib_lsp_get_spec.pro which is typically located under
TDAS_DIR/idl/themis/examples.
HISTORY:
2010-03-07: Created by Jianbao Tao (JBT), CU/LASP.
2012-07-26: JBT, SSL, UC Berkeley.
1. Removed the spectral averaging.
2. Added keywords *phase* and *ft* to output phase and FFT data.
VERSION:
$LastChangedBy: jianbao_tao $
$LastChangedDate: 2012-07-27 12:29:50 -0700 (Fri, 27 Jul 2012) $
$LastChangedRevision: 10753 $
$URL: svn+ssh://thmsvn@ambrosia.ssl.berkeley.edu/repos/spdsoft/trunk/projects/themis/spacecraft/fields/LASP/thm_lsp_get_spec.pro $
(See projects/themis/spacecraft/fields/LASP/thm_lsp_get_spec.pro)
FUNCTION: THM_LSP_MEDIAN_SMOOTH, x, nsmooth
NOT FOR GENERAL USE. CALLED BY THM_EFI...
ONLY FOR ISOLATED PARTICLE OR WAVE BURSTS
NOT FOR ENTIRE ORBIT.
PURPOSE:
Smooth data with median rather than average - removes spikes.
INPUT:
x -NEEDED. Array
nsmooth -OPTIONAL. Number of points in smoothing. DEFAULT = 11
-MUST BE ODD.
KEYWORDS:
HISTORY:
2009-04-28: REE.
(See projects/themis/spacecraft/fields/LASP/thm_lsp_median_smooth.pro)
FUNCTION: THM_LSP_NOTCH_SPIKES
NOT FOR GENERAL USE.
ONLY FOR ISOLATED PARTICLE OR WAVE BURSTS
NOT FOR ENTIRE ORBIT.
PURPOSE:
Remove the non-physical spiky signals in the efw data.
INPUT:
t -NEEDED. Time array
x -NEEDED. Data
per -NEEDED. Spin period.
tpks -NEEDED. Phase of spikes.
KEYWORD:
nfit -OPTIONAL. Number of points in the fit window. DFLT = 16 (256 for wave)
fit -OPTIONAL. If set, will perform a Gauss fit.
Amp -OPTIONAL. Estimated amplitudes. Same number of elements as tpks.
Diagnose -OPTIONAL. Plots fits and notchs.
Talk -OPTIONAL. Indicates number of fits and notchs.
HISTORY:
2009-05-30: REE. Broke out to run with wave burst.
(See projects/themis/spacecraft/fields/LASP/thm_lsp_notch_spikes.pro)
PRO: THM_LSP_REMOVE_SPIKES
NOT FOR GENERAL USE.
ONLY FOR ISOLATED PARTICLE OR WAVE BURSTS
NOT FOR ENTIRE ORBIT.
PURPOSE:
Remove the non-physical spiky signals in the efw data.
INPUT:
t -NEEDED. Time array
Ex, Ey, Ez -NEEDED. Data. Leave blank or set to zero to skip.
per -NEEDED. Spin period.
Efp -NEEDED for wave burst. Not needed for particle burst.
KEYWORD:
nwin -OPTIONAL. Number of points in spike search window. DFLT = 16
spikesig -OPTIONAL. Sigma of spikes. DFLT = 5
sigmin -OPTIONAL. Minimun of sigma. DFLT = 0 mV/m.
nfit -OPTIONAL. Number of points in the fit window. DFLT = 16
fit -OPTIONAL. If set, will fit spikes to Gaussian. DFLT = 0
HISTORY:
2009-05-12: REE. Complete rewrite for wave burst.
VERSION:
$LastChangedBy$
$LastChangedDate$
$LastChangedRevision$
$URL$
(See projects/themis/spacecraft/fields/LASP/thm_lsp_remove_spikes.pro)
FUNCTION: THM_LSP_REMOVE_SPINTONE, t, x, per, SpinPoly=SpinPoly, $
nsmooth=nsmooth, talk=talk
NOT FOR GENERAL USE. CALLED BY THM_EFI...
ONLY FOR ISOLATED PARTICLE OR WAVE BURSTS
NOT FOR ENTIRE ORBIT.
PURPOSE: GETS RID OF SPIN TONE AND OFFSET ON MEDIUM SIZE STRETCHES
INPUT:
t -NEEDED. Time array.
x -NEEDED. Data array.
per -NEEDED. Spin period. NOTE: input per/2 for 2nd harmonic, etc.
KEYWORDS:
-DEFAULT ACTION. If neither nsmooth or SpinPoly are set, sin
and cos fits over entire period are used.
nsmooth -OPTIONAL. If set, it takes precedence. It uses qfit
and "median smooth" sin and cos over
nsmooth points.
SpinPoly -OPTIONAL. If set (and nsmooth not set), it uses qfit
and polyfits sin and cos over nsmooth points.
talk -OPTIONAL. Plots diagnostics.
fail -OPTIONAL. If the removing fails, fail = 1.
Otherwise fail = 0.
OUTPUT: TPLOT STORE
HISTORY:
INITIAL VERSION: REE 99-03-25 (ff_remove_spintone)
Space Scienes Lab, UCBerkeley
MODIFICATION HISTORY:
08-10-31 Modified for THEMIS by REE
University of Colorado
2011-05-20: Added the keyword FAIL.
Jianbao Tao (JBT), CU/LASP.
(See projects/themis/spacecraft/fields/LASP/thm_lsp_remove_spintone.pro)
FUNCTION: THM_LSP_REMOVE_SPIN_EPOCH, t, x, per, talk=talk
NOT FOR GENERAL USE. CALLED BY THM_EFI_REMOVE_SPIKES
ONLY FOR ISOLATED PARTICLE OR WAVE BURSTS
NOT FOR ENTIRE ORBIT.
PURPOSE:
Remove the non-physical spiky signals in the efw data.
INPUT:
t -NEEDED. Time array
x -NEEDED. Data
per -NEEDED. Spin period.
KEYWORDS:
talk -OPTIONAL. Provides diagnostic plots.
HISTORY:
2008-11-08: Created by Jianbao Tao at LASP@CU-Boulder.
2008-11-10: The head comment was added.
2009-03-30: REE. Rewrite to test epoch analysis.
(See projects/themis/spacecraft/fields/LASP/thm_lsp_remove_spin_epoch.pro)
FUNCTION: THM_QFIT, data, phs, phsf=phsf, es=es, ec=ec, zero=zero,
do_sigma=do_sigma, sigma=sigma,
period=period, slide=slide, n_fitpts=n_fitpts,
out=out, max_err=max_err, bad_pts=bad_pts
PURPOSE: User unfriendly fit routine.
NOTE: REQUIRES CONTINUOUS BUFFERS!
NOT FOR GENERAL USE.
INPUT:
data - REQUIRED. A DATA ARRAY - NOT A STRUCTURE!
RECOMMEND: dat_in does not have NANS. One nan
may destroy entire period.
phs - REQUIRED. An array of phases.
IMPORTANT! All phases must be valid! No NANs!
KEYWORDS:
period - INPUT. Fit period. DEFAULT = 4pi
slide - NO LONGER AN OPTION!
NOTE: Slide will be forced to be pi/2.
n_fitpts - INPUT. Number of points per fit. DEFAULT = 64.
do_sigma - OPTION. /do_sigma fills sigma.
out - OPTION. /out fills bad_pts.
NOTE: FOR OUTLYER REJECTION, SEE BELOW!
max_err - INPUT. Maximum allowable error of bad_pts. DEFAULT=25 nT
OUTPUT:
phsf - OUTPUT. Phase of fit -> Time of fit.
es - OUTPUT. Sin phase of fit.
ec - OUTPUT. Cos phase of fit.
zero - OUTPUT. Zero level of fit.
sigma - OUTPUT. Deviation in nT. Only filled if /do_sigma
bad_pts - OUTPUT. Deviation in nT. Only filled if /out
CALLING:
SEE ff_magdc for an example. MUST USE fa_fields_bufs, and
ff_zero_crossing first.
OUT-LYING POINTS REJECTION:
The program must be iterated. For example, below shows one interation:
ff_qfit,data,phs,phsf=phsf,es=es,ec=ec,zero=zero,/out, bad_pts=bad_pts
data(bad_pts) = !values.f_nan
index = where(finite(data), n_index)
IF n_index GT 0 then BEGIN
data = data(index)
phs = phs(index)
ff_qfit,data,phs,phsf=phsf,es=es,ec=ec,zero=zero
ENDIF ELSE print, "Big trouble! No valid points."
HISTORY:
INITIAL VERSION: REE 97-10-05
Space Sciences Lab, UCBerkeley
2011-05-20: Added the keyword FAIL.
Jianbao Tao (JBT), CU/LASP.
(See projects/themis/spacecraft/fields/LASP/thm_qfit.pro)
NAME:
TRANGE_CLIP (PROCEDURE)
PURPOSE:
Utility to trim the time range of a tplot variable and remove
excess data.
CALLING SEQUENCE:
trange_clip,name, t1, t2, newname=newname, data_in=data_in,$
remove_match=remove_match, BadClip=BadClip
ARGUMENTS:
name: (INPUT, REQUIRED) Either tplot variable name
or data structure (if data_in keyword specified)
t1: (INPUT, REQUIRED) Start time (double, time since 1970)
t2: (INPUT, REQUIRED) Stop time (double, time since 1970)
KEYWORDS:
data_in: (INPUT, OPTIONAL) Set to specify that the input is a data
structure, not a tplot name.
remove_match: (INPUT, OPTIONAL) Removes data between t1 and t2
(default is keep only data between t1 and t2)
BadClip: (OUTPUT, OPTIONAL) A named variable of the status of the clip. If
the clip is valid, BadClip = 0, otherwise 1.
newname: (INPUT, OPTIONAL) A new name for storing the clip into a tplot
variable if the name is also a tplot variable name.
Outputs:
Modifies the specified tplot name or data structure.
Examples:
tvar = 'thd_fgh_dsl'
newname = tvar + '_clip'
trange_clip, tvar, t1, t2, newname = newname
HISTORY:
REE. 09-05-11. Changed to work with CLEAN EFP. Added BadClip=BadClip
VERSION:
$LastChangedBy$
$LastChangedDate$
$LastChangedRevision$
$URL$
(See projects/themis/spacecraft/fields/LASP/trange_clip.pro)
NAME:
VH (short for ViewHelp)
PURPOSE:
Make a shortcut to the IDL procedure DOC_LIBRARY
CALLING SEQUENCE
vh, pro_name
ARGUMENTS:
pro_name: (INPUT, REQUIRED) The name of the routine whose help document is
to be shown. It could be a string. For example, the following IDL
command shows the help info of the routine TPLOT.
vh, 'tplot'
It could also simply be the name of the routine (without the quotes),
if no variable with the same name as the routine is defined. For
example, if no variable is named TPLOT, the following command command
does the same thing as the one above.
vh, tplot
HISTORY:
2008-05-23: Created by Jianbao Tao (JBT), CU/LASP
2010-12-12: JBT, CU/LASP.
1. The help document was improved.
2012-07-20: JBT, SSL, UC Berkeley.
1. Added the VERSION to the documentation comment.
VERSION:
$LastChangedBy: jianbao_tao $
$LastChangedDate: 2012-07-23 17:44:33 -0700 (Mon, 23 Jul 2012) $
$LastChangedRevision: 10737 $
$URL: svn+ssh://thmsvn@ambrosia.ssl.berkeley.edu/repos/spdsoft/trunk/projects/themis/spacecraft/fields/LASP/vh.pro $
(See projects/themis/spacecraft/fields/LASP/vh.pro)
NAME:
FFT_CONVOLUTION
PURPOSE:
Computes the convolution of two or many 1D to 3D arrays by using the FFT.
DESCRIPTION:
This function computes in one time the convolution of an array with
one or more convolution kernels, e.g. F*K1*K2*...Kn.
1 dimensional to 3 dimensional arrays are supported.
CATEGORY:
Mathematics.
CALLING SEQUENCE:
Result = FFT_CONVOLUTION(array, kernel, /DOUBLE)
INPUTS:
Array: 1-D to 3-D array to be convolved.
Kernel: 1-D to 3-D array to be used as convolution kernel.
Kernel Must have the SAME SIZE as Array.
For multiple convolutions Kernel must have the dimension of Array + 1 and
must have one convolution kernel for each element of the Last dimension:
e.g.: for F(x,y)*K1(x,y)*K2(x,y) where F is of LX*LY size the kernel is
an array of [LX, LY, 2] in size.
The Kernel is always centered on the Array points to make convolution.
KEYWORD PARAMETERS:
DOUBLE: if set uses the double precision for the FFT.
MODIFICATION HISTORY:
Feb 2004 - Gianluca Li Causi & Massimo De Luca, INAF - Rome Astronomical Observatory
licausi@mporzio.astro.it
http://www.mporzio.astro.it/~licausi/
(See projects/themis/spacecraft/fields/LASP/fft_convolution.pro)