def test_center_brst_moments_data(self): data = mms_load_hpca(trange=['2015-10-16/13:06', '2015-10-16/13:07'], data_rate='brst') centered = mms_load_hpca(trange=['2015-10-16/13:06', '2015-10-16/13:07'], data_rate='brst', center_measurement=True, suffix='_centered') t, d = get_data('mms1_hpca_hplus_ion_bulk_velocity') c, d = get_data('mms1_hpca_hplus_ion_bulk_velocity_centered') self.assertTrue(np.round(c[0]-t[0], decimals=3) == 5.0)
def test_load_ion_omni(self): del_data('*') data = mms_load_hpca(trange=['2016-10-16/5:00', '2016-10-16/6:00'], datatype='ion') mms_hpca_calc_anodes(fov=[0, 360], probe='1') mms_hpca_spin_sum() self.assertTrue(data_exists('mms1_hpca_hplus_flux_elev_0-360_spin')) self.assertTrue(data_exists('mms1_hpca_heplus_flux_elev_0-360_spin')) self.assertTrue(data_exists('mms1_hpca_heplusplus_flux_elev_0-360_spin')) self.assertTrue(data_exists('mms1_hpca_oplus_flux_elev_0-360_spin'))
def test_load_ion_omni_suffix(self): del_data('*') data = mms_load_hpca( probe=2, trange=['2016-08-09/09:10', '2016-08-09/10:10:00'], datatype='ion', data_rate='brst', suffix='_brst') mms_hpca_calc_anodes(fov=[0, 360], probe=2, suffix='_brst') mms_hpca_spin_sum(probe=2, suffix='_brst', avg=True) self.assertTrue( data_exists('mms2_hpca_hplus_flux_brst_elev_0-360_spin'))
print(time_double('2015-10-16 13:06:00.00451')) # create new tplot variables with store_data from pytplot import store_data # save the B-field vector store_data('b_vector', data={'x': times, 'y': fgm_data[:, 0:3]}) # save the B-field magnitude store_data('b_mag', data={'x': times, 'y': fgm_data[:, 3]}) # the keywords are very flexible, e.g., from pyspedas.mms import mms_load_hpca, mms_load_eis, mms_load_feeps, mms_load_aspoc # specify multiple probes as integers, and multiple datatypes mms_load_hpca(probe=[1, 2, 4], data_rate='brst', datatype=['moments', 'ion'], trange=['2015-10-16/13:06', '2015-10-16/13:07']) mms_load_eis(probe='4', data_rate='brst', datatype=['phxtof', 'extof'], trange=['2015-10-16/13:06', '2015-10-16/13:07']) mms_load_feeps(get_support_data=True, probe=2, data_rate=['srvy', 'brst'], trange=['2015-10-16/13:06', '2015-10-16/13:07']) mms_load_aspoc() ''' ========================================================================== Plotting MMS Data in Python
def test_load_spdf_data(self): data = mms_load_hpca(trange=['2015-10-16', '2015-10-16/01:00'], spdf=True) self.assertTrue(data_exists('mms1_hpca_hplus_number_density'))