def test_pad_electrons_srvy_probe(self):
mms_load_feeps(probe=4)
mms_feeps_pad(probe=4)
self.assertTrue(
data_exists(
'mms4_epd_feeps_srvy_l2_electron_intensity_70-600keV_pad_spin')
)
self.assertTrue(
data_exists(
'mms4_epd_feeps_srvy_l2_electron_intensity_70-600keV_pad'))
def test_pad_ions_srvy(self):
mms_load_feeps(probe=4,
data_rate='brst',
datatype='ion',
trange=['2015-10-01/10:48:16', '2015-10-01/10:49:16'])
mms_feeps_pad(probe=4,
data_rate='brst',
datatype='ion',
angles_from_bfield=True)
self.assertTrue(
data_exists('mms4_epd_feeps_brst_l2_ion_intensity_70-600keV_pad'))
self.assertTrue(
data_exists(
'mms4_epd_feeps_brst_l2_ion_intensity_70-600keV_pad_spin'))
'''
This script loads and prints FEEPS brst-mode ion data at several data points.
This is meant to be called from the IDL test suite for comparison with the data loaded via IDL SPEDAS
'''
from pyspedas import mms_load_feeps
from pytplot import get_data
mms_load_feeps(trange=['2015-10-16/13:06', '2015-10-16/13:07'],
data_rate='brst')
t, d, e = get_data('mms1_epd_feeps_brst_l2_electron_intensity_omni')
# times
print(t[0:10].round(6).tolist())
# energy table
print(e.round(6).tolist())
# spot-check some data values
print(d[2, :].tolist())
print(d[200, :].tolist())
print(d[300, :].tolist())
print(d[350, :].tolist())
'''
This script loads and prints FEEPS srvy-mode ion data at several data points.
This is meant to be called from the IDL test suite for comparison with the data loaded via IDL SPEDAS
'''
from pyspedas import mms_load_feeps
from pytplot import get_data
mms_load_feeps(datatype='ion')
t, d, e = get_data('mms1_epd_feeps_srvy_l2_ion_intensity_omni')
# times
print(t[0:10].round(6).tolist())
# energy table
print(e.round(6).tolist())
# spot-check some data values
print(d[2, :].tolist())
print(d[2000, :].tolist())
print(d[5000, :].tolist())
print(d[10000, :].tolist())
print(d[15000, :].tolist())
'''
This script loads and prints FEEPS brst-mode ion data at several data points.
This is meant to be called from the IDL test suite for comparison with the data loaded via IDL SPEDAS
'''
from pyspedas import mms_load_feeps
from pytplot import get_data
mms_load_feeps(datatype='ion',
trange=['2015-10-16/13:06', '2015-10-16/13:07'],
data_rate='brst')
t, d, e = get_data('mms1_epd_feeps_brst_l2_ion_intensity_omni')
# times
print(t[0:10].round(6).tolist())
# energy table
print(e.round(6).tolist())
# spot-check some data values
print(d[2, :].tolist())
print(d[200, :].tolist())
print(d[300, :].tolist())
print(d[350, :].tolist())
'''
This script loads and prints FEEPS data at several data points.
This is meant to be called from the IDL test suite for comparison with the data loaded via IDL SPEDAS
'''
from pyspedas import mms_load_feeps
from pytplot import get_data
mms_load_feeps()
t, d, e = get_data('mms1_epd_feeps_srvy_l2_electron_intensity_omni')
# times
print(t[0:10].round(6).tolist())
# energy table
print(e.round(6).tolist())
# spot-check some data values
print(d[2, :].tolist())
print(d[2000, :].tolist())
print(d[5000, :].tolist())
print(d[10000, :].tolist())
print(d[15000, :].tolist())
'''
This crib sheet shows how to load and plot FEEPS data
'''
from pyspedas import mms_load_feeps, mms_feeps_pad
from pytplot import tplot
# load electron data
mms_load_feeps(trange=['2015-10-16', '2015-10-17'], datatype='electron')
# plot the omni-directional electron intensity
tplot([
'mms1_epd_feeps_srvy_l2_electron_intensity_omni',
'mms1_epd_feeps_srvy_l2_electron_intensity_omni_spin'
])
# calculate the pitch angle distribution
mms_feeps_pad()
# plot the PAD
tplot('mms1_epd_feeps_srvy_l2_electron_intensity_70-600keV_pad')
