.ipynb .pdf

repository open issue

JupyterHub

Contents

EFIx_LP_1B (Langmuir probe 2Hz)

Abstract: Access to the electric field instrument langmuir probe data (2Hz) (level 1b product).

%load_ext watermark
%watermark -i -v -p viresclient,pandas,xarray,matplotlib

Python implementation: CPython
Python version       : 3.9.7
IPython version      : 8.0.1

viresclient: 0.10.1
pandas     : 1.4.1
xarray     : 0.21.1
matplotlib : 3.5.1

from viresclient import SwarmRequest
import datetime as dt
import matplotlib.pyplot as plt

request = SwarmRequest()

EFIx_LP_1B product information

Measurements from the Langmuir Probe (LP) of the Electric Field Instrument (EFI) at 2Hz, for each Swarm spacecraft.

Documentation:

• https://earth.esa.int/web/guest/missions/esa-eo-missions/swarm/data-handbook/level-1b-product-definitions#EFIX_LP_1B_Product

Check what “EFI” data variables are available

request.available_collections("EFI", details=False)

{'EFI': ['SW_OPER_EFIA_LP_1B', 'SW_OPER_EFIB_LP_1B', 'SW_OPER_EFIC_LP_1B']}

request.available_measurements("EFI")

['U_orbit',
 'Ne',
 'Ne_error',
 'Te',
 'Te_error',
 'Vs',
 'Vs_error',
 'Flags_LP',
 'Flags_Ne',
 'Flags_Te',
 'Flags_Vs']

Fetch one day of EFI data

request.set_collection("SW_OPER_EFIA_LP_1B")
request.set_products(
    measurements=['U_orbit',
                 'Ne',
                 'Ne_error',
                 'Te',
                 'Te_error',
                 'Vs',
                 'Vs_error',
                 'Flags_LP',
                 'Flags_Ne',
                 'Flags_Te',
                 'Flags_Vs']
)
data = request.get_between(
    dt.datetime(2016,1,1),
    dt.datetime(2016,1,2)
)

data.sources

['SW_OPER_EFIA_LP_1B_20160101T000000_20160101T235959_0501_MDR_EFI_LP']

Load and plot using pandas/matplotlib

df = data.as_dataframe()
df.head()

	Radius	Latitude	U_orbit	Te_error	Spacecraft	Longitude	Te	Flags_Ne	Ne_error	Vs	Vs_error	Flags_LP	Flags_Te	Ne	Flags_Vs
Timestamp
2016-01-01 00:00:00.196999936	6833852.72	-72.511716	7604.407	9.999990e+09	A	92.799630	2945.20	20	9.999990e+09	-2.201	9.999990e+09	1	20	126188.4	20
2016-01-01 00:00:00.696000000	6833853.09	-72.543238	7604.398	9.999990e+09	A	92.813944	2891.97	20	9.999990e+09	-2.193	9.999990e+09	1	20	127792.9	20
2016-01-01 00:00:01.196999936	6833853.91	-72.574886	7604.389	9.999990e+09	A	92.828370	2921.96	20	9.999990e+09	-2.200	9.999990e+09	1	20	132515.3	20
2016-01-01 00:00:01.696000000	6833854.27	-72.606406	7604.380	9.999990e+09	A	92.842799	2936.74	20	9.999990e+09	-2.194	9.999990e+09	1	20	137933.0	20
2016-01-01 00:00:02.196999936	6833855.09	-72.638051	7604.371	9.999990e+09	A	92.857342	2870.09	20	9.999990e+09	-2.190	9.999990e+09	1	20	138913.8	20

df.plot(y=["Ne", "Te", "Vs"], subplots=True, figsize=(20,5));

df.plot(x="Latitude", y="Ne");

Load as xarray

ds = data.as_xarray()
ds

<xarray.Dataset>
Dimensions:     (Timestamp: 172776)
Coordinates:
  * Timestamp   (Timestamp) datetime64[ns] 2016-01-01T00:00:00.196999936 ... ...
Data variables: (12/15)
    Spacecraft  (Timestamp) object 'A' 'A' 'A' 'A' 'A' ... 'A' 'A' 'A' 'A' 'A'
    U_orbit     (Timestamp) float64 7.604e+03 7.604e+03 ... 7.634e+03 7.634e+03
    Longitude   (Timestamp) float64 92.8 92.81 92.83 ... -95.37 -95.37 -95.37
    Te          (Timestamp) float64 2.945e+03 2.892e+03 ... 2.527e+03 2.545e+03
    Flags_Ne    (Timestamp) uint8 20 20 20 20 20 20 20 ... 20 20 20 20 20 20 20
    Ne_error    (Timestamp) float64 1e+10 1e+10 1e+10 ... 1e+10 1e+10 1e+10
    ...          ...
    Flags_Vs    (Timestamp) uint8 20 20 20 20 20 20 20 ... 20 20 20 20 20 20 20
    Radius      (Timestamp) float64 6.834e+06 6.834e+06 ... 6.823e+06 6.823e+06
    Latitude    (Timestamp) float64 -72.51 -72.54 -72.57 ... 31.65 31.69 31.72
    Te_error    (Timestamp) float64 1e+10 1e+10 1e+10 ... 1e+10 1e+10 1e+10
    Vs          (Timestamp) float64 -2.201 -2.193 -2.2 ... -2.244 -2.239 -2.243
    Ne          (Timestamp) float64 1.262e+05 1.278e+05 ... 6.485e+04 6.43e+04
Attributes:
    Sources:         ['SW_OPER_EFIA_LP_1B_20160101T000000_20160101T235959_050...
    MagneticModels:  []
    RangeFilters:    []

xarray.Dataset

• Dimensions:
  • Timestamp: 172776
• Coordinates: (1)
  • Timestamp
    (Timestamp)
    datetime64[ns]
    2016-01-01T00:00:00.196999936 .....

    description :

    Time stamp

    array(['2016-01-01T00:00:00.196999936', '2016-01-01T00:00:00.696000000',
           '2016-01-01T00:00:01.196999936', ..., '2016-01-01T23:59:58.696000000',
           '2016-01-01T23:59:59.196999936', '2016-01-01T23:59:59.696000000'],
          dtype='datetime64[ns]')

• Data variables: (15)
  • Spacecraft
    (Timestamp)
    object
    'A' 'A' 'A' 'A' ... 'A' 'A' 'A' 'A'

    units :

    -

    description :

    Spacecraft identifier (values: 'A', 'B', 'C' or '-' if not available).

    array(['A', 'A', 'A', ..., 'A', 'A', 'A'], dtype=object)

  • U_orbit
    (Timestamp)
    float64
    7.604e+03 7.604e+03 ... 7.634e+03

    units :

    m/s

    description :

    Spacecraft velocity in the ITRF

    array([7604.407, 7604.398, 7604.389, ..., 7633.896, 7633.891, 7633.888])

  • Longitude
    (Timestamp)
    float64
    92.8 92.81 92.83 ... -95.37 -95.37

    units :

    deg

    description :

    Position in ITRF – Longitude

    array([ 92.7996305,  92.8139436,  92.8283703, ..., -95.3696257,
           -95.3696648, -95.3697037])

  • Te
    (Timestamp)
    float64
    2.945e+03 2.892e+03 ... 2.545e+03

    units :

    K

    description :

    Plasma electron temperature

    array([2945.2 , 2891.97, 2921.96, ..., 2568.07, 2526.62, 2545.17])

  • Flags_Ne
    (Timestamp)
    uint8
    20 20 20 20 20 ... 20 20 20 20 20

    units :

    -

    description :

    Flags characterizing the plasma density measurement

    array([20, 20, 20, ..., 20, 20, 20], dtype=uint8)

  • Ne_error
    (Timestamp)
    float64
    1e+10 1e+10 1e+10 ... 1e+10 1e+10

    units :

    cm-3

    description :

    Error estimate of plasma density (Ne)

    array([9.99999e+09, 9.99999e+09, 9.99999e+09, ..., 9.99999e+09,
           9.99999e+09, 9.99999e+09])

  • Flags_LP
    (Timestamp)
    uint8
    1 1 1 1 1 1 1 1 ... 1 1 1 1 1 1 1 1

    units :

    -

    description :

    Flags indicating the source/method of LP measurements (ne, Te, Vs)

    array([1, 1, 1, ..., 1, 1, 1], dtype=uint8)

  • Vs_error
    (Timestamp)
    float64
    1e+10 1e+10 1e+10 ... 1e+10 1e+10

    units :

    V

    description :

    Error estimate of spacecraft potential (U_pol)

    array([9.99999e+09, 9.99999e+09, 9.99999e+09, ..., 9.99999e+09,
           9.99999e+09, 9.99999e+09])

  • Flags_Te
    (Timestamp)
    uint8
    20 20 20 20 20 ... 20 20 20 20 20

    units :

    -

    description :

    Flags characterizing the electron temperature measurement

    array([20, 20, 20, ..., 20, 20, 20], dtype=uint8)

  • Flags_Vs
    (Timestamp)
    uint8
    20 20 20 20 20 ... 20 20 20 20 20

    units :

    -

    description :

    Flags characterizing the spacecraft potential measurement

    array([20, 20, 20, ..., 20, 20, 20], dtype=uint8)

  • Radius
    (Timestamp)
    float64
    6.834e+06 6.834e+06 ... 6.823e+06

    units :

    m

    description :

    Position in ITRF – Radius

    array([6833852.72, 6833853.09, 6833853.91, ..., 6822823.74, 6822819.66,
           6822817.73])

  • Latitude
    (Timestamp)
    float64
    -72.51 -72.54 ... 31.69 31.72

    units :

    deg

    description :

    Position in ITRF – Latitude

    array([-72.5117156, -72.5432383, -72.5748861, ...,  31.6548056,
            31.6869231,  31.7189125])

  • Te_error
    (Timestamp)
    float64
    1e+10 1e+10 1e+10 ... 1e+10 1e+10

    units :

    K

    description :

    Error estimate of plasma electron temperature (Te)

    array([9.99999e+09, 9.99999e+09, 9.99999e+09, ..., 9.99999e+09,
           9.99999e+09, 9.99999e+09])

  • Vs
    (Timestamp)
    float64
    -2.201 -2.193 ... -2.239 -2.243

    units :

    V

    description :

    Spacecraft potential

    array([-2.201, -2.193, -2.2  , ..., -2.244, -2.239, -2.243])

  • Ne
    (Timestamp)
    float64
    1.262e+05 1.278e+05 ... 6.43e+04

    units :

    cm-3

    description :

    Plasma density (electron)

    array([126188.4, 127792.9, 132515.3, ...,  64617.9,  64852.5,  64298.7])

• Attributes: (3)

  Sources :

  ['SW_OPER_EFIA_LP_1B_20160101T000000_20160101T235959_0501_MDR_EFI_LP']

  MagneticModels :

  []

  RangeFilters :

  []

fig, (ax1, ax2) = plt.subplots(figsize=(10, 5), nrows=2, sharex=True)
def subplot(_ax, da):
    """Plot a given DataArray on a given axes"""
    _ax.plot(da)
    _ax.set_ylabel(f"{da.name} [{da.units}]")
for var, ax in zip(("Ne", "Te"), (ax1, ax2)):
    subplot(ax, ds[var])

previous

MAGxHR_1B (Magnetic field 50Hz)

next

IPDxIRR_2F (Ionospheric plasma densities)