THEMIS Particle Moments (Ground): Difference between revisions

From SPEDAS Wiki
Jump to navigation Jump to search
mNo edit summary
(various updates)
Line 1: Line 1:
[[File:peef_moments_example.png|200px|thumb|right|Density, velocity, and temperature for ESA electrons.]]
[[File:peef_moments_example.png|200px|thumb|right|Density, velocity, and temperature for ESA electrons.]]


Moments of the particle distribution may be calculated manually.  This allows for manual calibration and specification of input parameters.  Moments may be produced after selecting a spacecraft, particle data type, and time range.  Some products also require magnetic field data to be present.
Moments of the particle distribution may be calculated on the ground from raw ESA and/or SST particle data.  This allows for manual calibration and specification of input parameters.  Moments can be generated with command line routines or through the [[GUI|SPEDAS GUI]].  Some products also require magnetic field data to be present.


For command line usage see:
For command line usage see:
Line 13: Line 13:


==Overview==
==Overview==
Particle moments are produced by thm_part_products from [[Level 0 particle data|raw particle data]].  All data products are calculated simultaneously for a particular data type.
Moments are produced with ''thm_part_products'' from [[Level 0 particle data|raw particle data]].  All data products are calculated simultaneously for a particular data type.  Particle data must be loaded explicitly beforehand.


<div style="border: 1px solid LightGray; background-color:#F9F9F9; padding: 2px 10px; margin-left: 5px; overflow:hidden">
<div style="border: 1px solid LightGray; background-color:#F9F9F9; padding: 2px 10px; margin-left: 5px; overflow:hidden">
SPEDAS>  thm_part_load, probe='a', datatype='peif', trange=[t0,t1]
SPEDAS>  thm_part_products, probe='a', datatype='peif', trange=[t0,t1], '''outputs = 'moments' '''
SPEDAS>  thm_part_products, probe='a', datatype='peif', trange=[t0,t1], '''outputs = 'moments' '''
</div>
</div>
Line 28: Line 30:
Products marked with a * require magnetic field data to be loaded.
Products marked with a * require magnetic field data to be loaded.


All "vector" data types are 3-vectors.  Symmetric tensors are stored in six dimensions: xx, yy, zz, xy, xz, yz.
All "vector" data types are 3-vectors: <x,y,z>.  Symmetric tensors are stored in six dimensions: <xx, yy, zz, xy, xz, yz>.


{| class="wikitable" style="border: 1px solid LightGray; background-color:#F9F9F9; margin-left:5px;"
{| class="wikitable" style="border: 1px solid LightGray; background-color:#F9F9F9; margin-left:5px;"
Line 42: Line 44:
|scalar
|scalar
|-
|-
|current
|sc_current
|Current
|Total number flux
|
|1/(s&middot;cm<sup>2</sup>)
|scalar
|scalar
|-
|-
Line 117: Line 119:
====Spacecraft potential====
====Spacecraft potential====
[[File:pot_dens_compare.png|200px|thumb|right|Spacecraft potential plotted with corrected and uncorrected ESA density data]]
[[File:pot_dens_compare.png|200px|thumb|right|Spacecraft potential plotted with corrected and uncorrected ESA density data]]
Increased spacecraft potential may affect the lower energy bins for ESA data.  A tplot variable containing spacecraft potential can be specified when generating particle data products.  If no variable is specified then '''"th?_pxxm_pot"''' will be used, if present.
Accumulated charge on the spacecraft will affect low energy ESA data.  A tplot variable containing spacecraft potential can be used to correct for this when generating ground moments.  If no variable is specified then '''"th?_pxxm_pot"''' will be used automatically, if present. Spacecraft potential can be loaded with the [[THEMIS Particle Moments (Spacecraft)|on-board moments]].
Spacecraft potential must be loaded separately:
 
<div style="border: 1px solid LightGray; background-color:#F9F9F9; padding: 2px 10px; margin-left: 5px; overflow:hidden;">
<div style="border: 1px solid LightGray; background-color:#F9F9F9; padding: 2px 10px; margin-left: 5px; overflow:hidden;">
SPEDAS> thm_load_mom, probe='a', trange=[t0,t1]
SPEDAS> thm_load_mom, probe='a', trange=[t0,t1]
Line 127: Line 129:
====ESA background subtraction====
====ESA background subtraction====
Background subtraction for ESA data is enabled by default when producing spectra.  For background removal options see [[ESA Background Removal]].  For example usage see:
Background subtraction for ESA data is enabled by default when producing spectra.  For background removal options see [[ESA Background Removal]].  For example usage see:
<div style="border: 1px solid LightGray; background-color:#F9F9F9; padding: 2px 10px; margin-left: 5px; overflow:hidden;">
<div style="border: 1px solid LightGray; background-color:#F9F9F9; padding: 2px 10px; margin-left: 5px; overflow:hidden;">
/projects/THEMIS/examples/advanced/thm_crib_esa_bgnd_remove.pro
/projects/THEMIS/examples/advanced/thm_crib_esa_bgnd_remove.pro
Line 133: Line 136:
====SST contamination====
====SST contamination====
In most cases, contaminated SST bins should be masked and interpolated over automatically.  If the default set of bins is inadequate then a list may be specified instead.  For more information see [[SST Contamination Removal]].
In most cases, contaminated SST bins should be masked and interpolated over automatically.  If the default set of bins is inadequate then a list may be specified instead.  For more information see [[SST Contamination Removal]].
<div style="border: 1px solid LightGray; background-color:#F9F9F9; padding: 2px 10px; margin-left: 5px; overflow:hidden;">
<div style="border: 1px solid LightGray; background-color:#F9F9F9; padding: 2px 10px; margin-left: 5px; overflow:hidden;">
SPEDAS> thm_part_products, probe='a', datatype='psif', trange=[t0,t1], '''sst_sun_bins=[0,8,16,24,32,33,34,40,47,48,49,50,55,56,57]'''
SPEDAS> thm_part_products, probe='a', datatype='psif', trange=[t0,t1], '''sst_sun_bins=[0,8,16,24,32,33,34,40,47,48,49,50,55,56,57]'''
Line 139: Line 143:


==Uncertainty estimates==
==Uncertainty estimates==
The statistical uncertainty based on the original measurement in counts can be propagated through the moment calculation and output as tplot variables.  The error is taken to be the standard deviation of the measurement assuming they are Poisson distributed.
<div style="border: 1px solid LightGray; background-color:#F9F9F9; padding: 2px 10px; margin-left: 5px; overflow:hidden;">
SPEDAS> thm_part_products, probe='a', datatype='psif', trange=[t0,t1], '''/get_error'''
</div>





Revision as of 01:15, 28 May 2015

Density, velocity, and temperature for ESA electrons.

Moments of the particle distribution may be calculated on the ground from raw ESA and/or SST particle data. This allows for manual calibration and specification of input parameters. Moments can be generated with command line routines or through the SPEDAS GUI. Some products also require magnetic field data to be present.

For command line usage see:

/projects/THEMIS/examples/basic/thm_crib_part_products.pro


For spacecraft processed moments see on-board moments.


Overview

Moments are produced with thm_part_products from raw particle data. All data products are calculated simultaneously for a particular data type. Particle data must be loaded explicitly beforehand.

SPEDAS> thm_part_load, probe='a', datatype='peif', trange=[t0,t1]

SPEDAS> thm_part_products, probe='a', datatype='peif', trange=[t0,t1], outputs = 'moments'

Support data will be used automatically if it is present. By default the software will use "th?_fgs" as the magnetic field and "th?_pxxm_pot" as the spacecraft potential ("?" denoting probe designation). If no support data is found it will be ignored; any derived quantities will not be calculated. Keywords may be used to specify support data:

SPEDAS> thm_part_products, probe='a', datatype='peif', trange=[t0,t1], outputs = 'moments', mag_name='tha_fgl', sc_pot_name='tha_peef_sc_pot'

Data Products

Products marked with a * require magnetic field data to be loaded.

All "vector" data types are 3-vectors: <x,y,z>. Symmetric tensors are stored in six dimensions: <xx, yy, zz, xy, xz, yz>.

Variable Name Description Units Data type
density Density 1/cm3 scalar
sc_current Total number flux 1/(s·cm2) scalar
flux Number flux 1/(s·cm2) vector
eflux Energy flux eV/(s·cm2) vector
velocity Velocity km/s vector
t3 Diagonalized temperature eV vector
magt3 Magnetic field aligned temperature* eV vector
avgtemp Average temperature eV scalar
vthermal Thermal velocity km/s vector
mftens Momentum flux tensor eV/cm3 symmetric tensor
ptens Pressure tensor eV/cm3 symmetric tensor
symm Direction of pressure tensor symmetry n/a vector
symm_phi Phi (longitude) direction of symmetry vector ° scalar
symm_theta Theta (latitude) direction of symmetry vector ° scalar
symm_ang Angle between symmetry vector and magnetic field* ° scalar


Contamination

Spacecraft potential

Spacecraft potential plotted with corrected and uncorrected ESA density data

Accumulated charge on the spacecraft will affect low energy ESA data. A tplot variable containing spacecraft potential can be used to correct for this when generating ground moments. If no variable is specified then "th?_pxxm_pot" will be used automatically, if present. Spacecraft potential can be loaded with the on-board moments.

SPEDAS> thm_load_mom, probe='a', trange=[t0,t1]

SPEDAS> thm_part_products, probe='a', datatype='peef', trange=[t0,t1], output='moments', sc_pot_name='tha_pxxm_pot'

ESA background subtraction

Background subtraction for ESA data is enabled by default when producing spectra. For background removal options see ESA Background Removal. For example usage see:

/projects/THEMIS/examples/advanced/thm_crib_esa_bgnd_remove.pro

SST contamination

In most cases, contaminated SST bins should be masked and interpolated over automatically. If the default set of bins is inadequate then a list may be specified instead. For more information see SST Contamination Removal.

SPEDAS> thm_part_products, probe='a', datatype='psif', trange=[t0,t1], sst_sun_bins=[0,8,16,24,32,33,34,40,47,48,49,50,55,56,57]


Uncertainty estimates

The statistical uncertainty based on the original measurement in counts can be propagated through the moment calculation and output as tplot variables. The error is taken to be the standard deviation of the measurement assuming they are Poisson distributed.

SPEDAS> thm_part_products, probe='a', datatype='psif', trange=[t0,t1], /get_error


Resources

ESA instrument page

SST instrument page