.ipynb .pdf

repository open issue

JupyterHub

Contents

MAGxLR_1B (Magnetic field 1Hz)

Abstract: Access to the low rate (1Hz) magnetic data (level 1b product), together with geomagnetic model evaluations (level 2 products).

%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.11.0
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()

Product information

This is one of the main products from Swarm - the 1Hz measurements of the magnetic field vector (B_NEC) and total intensity (F). These are derived from the Vector Field Magnetometer (VFM) and Absolute Scalar Magnetomer (ASM).

Documentation:

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

Measurements are available through VirES as part of collections with names containing MAGx_LR, for each Swarm spacecraft:

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

{'MAG': ['SW_OPER_MAGA_LR_1B', 'SW_OPER_MAGB_LR_1B', 'SW_OPER_MAGC_LR_1B']}

The measurements can be used together with geomagnetic model evaluations as shall be shown below.

Check what “MAG” data variables are available

request.available_measurements("MAG")

['F',
 'dF_AOCS',
 'dF_other',
 'F_error',
 'B_VFM',
 'B_NEC',
 'dB_Sun',
 'dB_AOCS',
 'dB_other',
 'B_error',
 'q_NEC_CRF',
 'Att_error',
 'Flags_F',
 'Flags_B',
 'Flags_q',
 'Flags_Platform',
 'ASM_Freq_Dev']

Check the names of available models

request.available_models(details=False)

['IGRF',
 'LCS-1',
 'MF7',
 'CHAOS-Core',
 'CHAOS-Static',
 'CHAOS-MMA-Primary',
 'CHAOS-MMA-Secondary',
 'MCO_SHA_2C',
 'MCO_SHA_2D',
 'MLI_SHA_2C',
 'MLI_SHA_2D',
 'MLI_SHA_2E',
 'MMA_SHA_2C-Primary',
 'MMA_SHA_2C-Secondary',
 'MMA_SHA_2F-Primary',
 'MMA_SHA_2F-Secondary',
 'MIO_SHA_2C-Primary',
 'MIO_SHA_2C-Secondary',
 'MIO_SHA_2D-Primary',
 'MIO_SHA_2D-Secondary',
 'AMPS',
 'MCO_SHA_2X',
 'CHAOS',
 'CHAOS-MMA',
 'MMA_SHA_2C',
 'MMA_SHA_2F',
 'MIO_SHA_2C',
 'MIO_SHA_2D',
 'SwarmCI']

Fetch some MAG data and models

We can fetch the data and the model predictions (evaluated on demand) at the same time. We can also subsample the data - here we subsample it to 10-seconds by specifying the “PT10S” sampling_step.

request.set_collection("SW_OPER_MAGA_LR_1B")
request.set_products(
    measurements=["F", "B_NEC"],
    models=["CHAOS-Core", "MCO_SHA_2D"],
    sampling_step="PT10S"
)
data = request.get_between(
    # 2014-01-01 00:00:00
    start_time = dt.datetime(2014,1,1, 0),
    # 2014-01-01 01:00:00
    end_time = dt.datetime(2014,1,1, 1)
)

See a list of the source files

data.sources

['SW_OPER_MAGA_LR_1B_20140101T000000_20140101T235959_0505_MDR_MAG_LR',
 'SW_OPER_MCO_SHA_2D_20131126T000000_20180101T000000_0401',
 'SW_OPER_MCO_SHA_2X_19970101T000000_20220807T235959_0711']

Load as a pandas dataframe

Use expand=True to extract vectors (B_NEC…) as separate columns (…_N, …_E, …_C)

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

	Radius	Latitude	Longitude	F_CHAOS-Core	F	F_MCO_SHA_2D	Spacecraft	B_NEC_CHAOS-Core_N	B_NEC_CHAOS-Core_E	B_NEC_CHAOS-Core_C	B_NEC_N	B_NEC_E	B_NEC_C	B_NEC_MCO_SHA_2D_N	B_NEC_MCO_SHA_2D_E	B_NEC_MCO_SHA_2D_C
Timestamp
2014-01-01 00:00:00	6878309.22	-1.228938	-14.116674	22874.436246	22867.5503	22874.211509	A	20113.226153	-4127.051332	-10082.926887	20103.5246	-4126.2621	-10086.9888	20113.623921	-4127.463956	-10081.454567
2014-01-01 00:00:10	6878381.17	-1.862521	-14.131424	22821.171065	22814.5656	22820.941425	A	19824.732246	-4162.697190	-10509.890211	19815.0914	-4160.9933	-10514.4074	19825.161844	-4163.127549	-10508.410652
2014-01-01 00:00:20	6878452.05	-2.496090	-14.146155	22769.599713	22763.2585	22769.369161	A	19533.089786	-4197.081116	-10922.343366	19523.4946	-4195.1968	-10926.9664	19533.553905	-4197.529054	-10920.860481
2014-01-01 00:00:30	6878521.87	-3.129644	-14.160861	22719.465912	22713.3703	22719.238240	A	19238.842861	-4230.136455	-11320.203278	19229.2386	-4228.4747	-11324.8335	19239.343572	-4230.601819	-11318.721366
2014-01-01 00:00:40	6878590.61	-3.763184	-14.175534	22670.525965	22664.7202	22670.304681	A	18942.536458	-4261.797994	-11703.424196	18932.8807	-4260.8424	-11708.0897	18943.075144	-4262.280634	-11701.947796

… or as an xarray dataset:

ds = data.as_xarray()
ds

<xarray.Dataset>
Dimensions:           (Timestamp: 360, NEC: 3)
Coordinates:
  * Timestamp         (Timestamp) datetime64[ns] 2014-01-01 ... 2014-01-01T00...
  * NEC               (NEC) <U1 'N' 'E' 'C'
Data variables:
    Spacecraft        (Timestamp) object 'A' 'A' 'A' 'A' 'A' ... 'A' 'A' 'A' 'A'
    Radius            (Timestamp) float64 6.878e+06 6.878e+06 ... 6.868e+06
    Latitude          (Timestamp) float64 -1.229 -1.863 -2.496 ... 48.14 48.77
    Longitude         (Timestamp) float64 -14.12 -14.13 -14.15 ... 153.6 153.6
    F_CHAOS-Core      (Timestamp) float64 2.287e+04 2.282e+04 ... 4.02e+04
    F                 (Timestamp) float64 2.287e+04 2.281e+04 ... 4.021e+04
    F_MCO_SHA_2D      (Timestamp) float64 2.287e+04 2.282e+04 ... 4.021e+04
    B_NEC_CHAOS-Core  (Timestamp, NEC) float64 2.011e+04 ... 3.557e+04
    B_NEC_MCO_SHA_2D  (Timestamp, NEC) float64 2.011e+04 ... 3.557e+04
    B_NEC             (Timestamp, NEC) float64 2.01e+04 -4.126e+03 ... 3.558e+04
Attributes:
    Sources:         ['SW_OPER_MAGA_LR_1B_20140101T000000_20140101T235959_050...
    MagneticModels:  ["CHAOS-Core = 'CHAOS-Core'(max_degree=20,min_degree=1)"...
    AppliedFilters:  []

xarray.Dataset

• Dimensions:
  • Timestamp: 360
  • NEC: 3
• Coordinates: (2)
  • Timestamp
    (Timestamp)
    datetime64[ns]
    2014-01-01 ... 2014-01-01T00:59:50

    description :

    Time stamp

    array(['2014-01-01T00:00:00.000000000', '2014-01-01T00:00:10.000000000',
           '2014-01-01T00:00:20.000000000', ..., '2014-01-01T00:59:30.000000000',
           '2014-01-01T00:59:40.000000000', '2014-01-01T00:59:50.000000000'],
          dtype='datetime64[ns]')

  • NEC
    (NEC)
    <U1
    'N' 'E' 'C'

    units :

    description :

    NEC frame - North, East, Centre (down)

    array(['N', 'E', 'C'], dtype='<U1')

• Data variables: (10)
  • 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', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A'], dtype=object)

  • Radius
    (Timestamp)
    float64
    6.878e+06 6.878e+06 ... 6.868e+06

    units :

    m

    description :

    Position in ITRF – Radius

    array([6878309.22, 6878381.17, 6878452.05, 6878521.87, 6878590.61,
           6878658.27, 6878724.84, 6878790.32, 6878854.7 , 6878917.97,
           6878980.13, 6879041.18, 6879101.12, 6879159.94, 6879217.63,
           6879274.2 , 6879329.64, 6879383.95, 6879437.14, 6879489.19,
           6879540.12, 6879589.91, 6879638.58, 6879686.11, 6879732.52,
           6879777.8 , 6879821.95, 6879864.98, 6879906.89, 6879947.68,
           6879987.36, 6880025.93, 6880063.39, 6880099.74, 6880135.  ,
           6880169.17, 6880202.24, 6880234.24, 6880265.15, 6880295.  ,
           6880323.77, 6880351.49, 6880378.17, 6880403.79, 6880428.38,
           6880451.94, 6880474.48, 6880496.01, 6880516.53, 6880536.06,
           6880554.61, 6880572.17, 6880588.77, 6880604.42, 6880619.11,
           6880632.87, 6880645.71, 6880657.62, 6880668.64, 6880678.75,
           6880687.99, 6880696.35, 6880703.85, 6880710.51, 6880716.32,
           6880721.31, 6880725.49, 6880728.87, 6880731.46, 6880733.28,
           6880734.33, 6880734.64, 6880734.2 , 6880733.05, 6880731.18,
           6880728.61, 6880725.36, 6880721.44, 6880716.86, 6880711.63,
           6880705.77, 6880699.3 , 6880692.22, 6880684.55, 6880676.31,
           6880667.5 , 6880658.13, 6880648.24, 6880637.81, 6880626.88,
           6880615.45, 6880603.53, 6880591.15, 6880578.3 , 6880565.01,
           6880551.29, 6880537.15, 6880522.61, 6880507.67, 6880492.35,
    ...
           6875648.21, 6875587.21, 6875525.54, 6875463.22, 6875400.22,
           6875336.57, 6875272.26, 6875207.3 , 6875141.69, 6875075.42,
           6875008.52, 6874940.98, 6874872.8 , 6874804.  , 6874734.58,
           6874664.55, 6874593.9 , 6874522.66, 6874450.83, 6874378.4 ,
           6874305.41, 6874231.84, 6874157.71, 6874083.03, 6874007.82,
           6873932.07, 6873855.8 , 6873779.02, 6873701.74, 6873623.97,
           6873545.73, 6873467.02, 6873387.86, 6873308.26, 6873228.24,
           6873147.8 , 6873066.96, 6872985.74, 6872904.15, 6872822.2 ,
           6872739.91, 6872657.3 , 6872574.37, 6872491.16, 6872407.66,
           6872323.91, 6872239.9 , 6872155.68, 6872071.23, 6871986.6 ,
           6871901.79, 6871816.82, 6871731.71, 6871646.48, 6871561.14,
           6871475.72, 6871390.23, 6871304.7 , 6871219.15, 6871133.59,
           6871048.04, 6870962.54, 6870877.08, 6870791.71, 6870706.43,
           6870621.28, 6870536.27, 6870451.42, 6870366.75, 6870282.3 ,
           6870198.07, 6870114.09, 6870030.39, 6869946.98, 6869863.89,
           6869781.14, 6869698.75, 6869616.74, 6869535.15, 6869453.98,
           6869373.26, 6869293.01, 6869213.26, 6869134.03, 6869055.34,
           6868977.21, 6868899.67, 6868822.72, 6868746.41, 6868670.75,
           6868595.76, 6868521.47, 6868447.89, 6868375.05, 6868302.96,
           6868231.66, 6868161.16, 6868091.48, 6868022.65, 6867954.68])

  • Latitude
    (Timestamp)
    float64
    -1.229 -1.863 ... 48.14 48.77

    units :

    deg

    description :

    Position in ITRF – Latitude

    array([ -1.2289377,  -1.8625209,  -2.4960899,  -3.1296444,  -3.7631844,
            -4.3967096,  -5.03022  ,  -5.6637153,  -6.2971954,  -6.9306601,
            -7.5641093,  -8.1975429,  -8.8309606,  -9.4643624, -10.097748 ,
           -10.7311174, -11.3644703, -11.9978066, -12.6311262, -13.2644288,
           -13.8977145, -14.5309829, -15.1642339, -15.7974674, -16.4306832,
           -17.0638812, -17.6970611, -18.3302228, -18.9633662, -19.5964911,
           -20.2295973, -20.8626845, -21.4957527, -22.1288017, -22.7618312,
           -23.394841 , -24.027831 , -24.6608009, -25.2937506, -25.9266797,
           -26.5595882, -27.1924756, -27.8253419, -28.4581867, -29.0910098,
           -29.7238109, -30.3565897, -30.9893459, -31.6220792, -32.2547894,
           -32.887476 , -33.5201387, -34.1527771, -34.785391 , -35.4179798,
           -36.0505431, -36.6830806, -37.3155917, -37.9480759, -38.5805328,
           -39.2129618, -39.8453623, -40.4777338, -41.1100756, -41.742387 ,
           -42.3746673, -43.0069157, -43.6391316, -44.271314 , -44.903462 ,
           -45.5355748, -46.1676512, -46.7996903, -47.4316908, -48.0636517,
           -48.6955715, -49.3274489, -49.9592825, -50.5910706, -51.2228117,
           -51.8545038, -52.4861451, -53.1177336, -53.7492669, -54.3807427,
           -55.0121585, -55.6435114, -56.2747986, -56.9060167, -57.5371624,
           -58.1682318, -58.7992209, -59.4301252, -60.06094  , -60.69166  ,
           -61.3222795, -61.9527923, -62.5831914, -63.2134696, -63.8436185,
    ...
           -14.0397279, -13.4057211, -12.7716953, -12.1376508, -11.5035876,
           -10.8695063, -10.2354069,  -9.6012897,  -8.967155 ,  -8.333003 ,
            -7.6988341,  -7.0646483,  -6.4304461,  -5.7962276,  -5.161993 ,
            -4.5277428,  -3.8934769,  -3.2591959,  -2.6248998,  -1.990589 ,
            -1.3562637,  -0.7219241,  -0.0875706,   0.5467966,   1.1811773,
             1.8155712,   2.449978 ,   3.0843975,   3.7188294,   4.3532734,
             4.9877293,   5.6221967,   6.2566755,   6.8911652,   7.5256657,
             8.1601767,   8.7946978,   9.4292288,  10.0637693,  10.6983191,
            11.3328779,  11.9674454,  12.6020212,  13.2366051,  13.8711967,
            14.5057958,  15.1404019,  15.7750149,  16.4096342,  17.0442597,
            17.6788909,  18.3135275,  18.9481691,  19.5828155,  20.2174661,
            20.8521207,  21.4867788,  22.1214401,  22.7561042,  23.3907706,
            24.025439 ,  24.6601089,  25.2947799,  25.9294516,  26.5641235,
            27.1987952,  27.8334661,  28.4681359,  29.102804 ,  29.7374699,
            30.3721332,  31.0067932,  31.6414495,  32.2761015,  32.9107485,
            33.5453901,  34.1800256,  34.8146544,  35.4492757,  36.0838891,
            36.7184936,  37.3530887,  37.9876736,  38.6222475,  39.2568096,
            39.891359 ,  40.525895 ,  41.1604165,  41.7949226,  42.4294124,
            43.0638848,  43.6983388,  44.3327731,  44.9671867,  45.6015782,
            46.2359464,  46.87029  ,  47.5046074,  48.1388972,  48.7731577])

  • Longitude
    (Timestamp)
    float64
    -14.12 -14.13 ... 153.6 153.6

    units :

    deg

    description :

    Position in ITRF – Longitude

    array([-14.1166743, -14.1314243, -14.1461554, -14.1608609, -14.1755342,
           -14.1901685, -14.2047571, -14.2192932, -14.2337699, -14.2481804,
           -14.2625177, -14.2767747, -14.2909444, -14.3050195, -14.3189927,
           -14.3328566, -14.3466037, -14.3602264, -14.3737168, -14.3870671,
           -14.4002693, -14.413315 , -14.4261961, -14.4389038, -14.4514295,
           -14.4637642, -14.4758987, -14.4878238, -14.4995298, -14.5110068,
           -14.5222447, -14.5332332, -14.5439615, -14.5544186, -14.5645933,
           -14.5744739, -14.5840483, -14.5933041, -14.6022285, -14.6108083,
           -14.6190299, -14.6268789, -14.6343409, -14.6414006, -14.6480422,
           -14.6542495, -14.6600055, -14.6652927, -14.6700926, -14.6743864,
           -14.6781542, -14.6813754, -14.6840285, -14.6860912, -14.6875401,
           -14.6883507, -14.6884977, -14.6879543, -14.6866927, -14.6846837,
           -14.6818967, -14.6782997, -14.673859 , -14.6685394, -14.6623038,
           -14.6551131, -14.6469263, -14.6377003, -14.6273894, -14.6159457,
           -14.6033184, -14.589454 , -14.5742958, -14.557784 , -14.539855 ,
           -14.5204415, -14.4994721, -14.476871 , -14.4525575, -14.4264457,
           -14.3984442, -14.3684554, -14.3363752, -14.3020922, -14.2654873,
           -14.2264328, -14.1847917, -14.1404169, -14.0931501, -14.0428207,
           -13.989245 , -13.9322242, -13.8715433, -13.8069696, -13.73825  ,
           -13.6651097, -13.5872493, -13.5043418, -13.41603  , -13.3219222,
    ...
           154.3456616, 154.3325292, 154.3192466, 154.3058219, 154.292263 ,
           154.2785777, 154.2647738, 154.2508587, 154.23684  , 154.222725 ,
           154.208521 , 154.194235 , 154.1798741, 154.1654455, 154.1509559,
           154.1364123, 154.1218215, 154.1071903, 154.0925253, 154.0778333,
           154.063121 , 154.0483951, 154.0336621, 154.0189287, 154.0042016,
           153.9894873, 153.9747926, 153.9601241, 153.9454886, 153.9308927,
           153.9163432, 153.901847 , 153.8874109, 153.8730419, 153.8587469,
           153.8445331, 153.8304075, 153.8163775, 153.8024503, 153.7886334,
           153.7749344, 153.761361 , 153.7479209, 153.7346222, 153.7214729,
           153.7084812, 153.6956557, 153.6830049, 153.6705376, 153.6582629,
           153.6461898, 153.634328 , 153.622687 , 153.6112768, 153.6001076,
           153.5891898, 153.5785343, 153.5681521, 153.5580547, 153.5482538,
           153.5387615, 153.5295904, 153.5207534, 153.5122639, 153.5041357,
           153.496383 , 153.4890208, 153.4820643, 153.4755295, 153.4694328,
           153.4637914, 153.4586231, 153.4539464, 153.4497804, 153.4461454,
           153.443062 , 153.4405521, 153.4386384, 153.4373444, 153.4366949,
           153.4367158, 153.437434 , 153.4388778, 153.4410767, 153.4440619,
           153.4478658, 153.4525225, 153.4580678, 153.4645395, 153.4719772,
           153.4804224, 153.4899192, 153.5005138, 153.5122551, 153.5251946,
           153.5393869, 153.5548894, 153.5717634, 153.5900734, 153.6098879])

  • F_CHAOS-Core
    (Timestamp)
    float64
    2.287e+04 2.282e+04 ... 4.02e+04

    units :

    nT

    description :

    Magnetic field intensity, calculated by the CHAOS-Core spherical harmonic model

    array([22874.43624621, 22821.17106481, 22769.59971337, 22719.46591166,
           22670.52596496, 22622.54967897, 22575.32183564, 22528.64304193,
           22482.33072681, 22436.21963862, 22390.16217088, 22344.02870303,
           22297.70765105, 22251.10548042, 22204.14635522, 22156.77151398,
           22108.93924075, 22060.62391357, 22011.8153655 , 21962.51852238,
           21912.75187454, 21862.54755391, 21811.94953231, 21761.01369558,
           21709.80609701, 21658.40285374, 21606.8889278 , 21555.35728442,
           21503.9083919 , 21452.64938365, 21401.69328829, 21351.15860276,
           21301.1686272 , 21251.85096256, 21203.33687346, 21155.76125822,
           21109.2622065 , 21063.98030763, 21020.059182  , 20977.64450599,
           20936.8847003 , 20897.930137  , 20860.93360176, 20826.0505287 ,
           20793.43838892, 20763.25746472, 20735.67070275, 20710.84394774,
           20688.94622584, 20670.14966433, 20654.62994294, 20642.56647779,
           20634.14212553, 20629.54369107, 20628.9619615 , 20632.59123883,
           20640.6296496 , 20653.2789663 , 20670.74384286, 20693.2322135 ,
           20720.95401614, 20754.12138269, 20792.9474765 , 20837.64593927,
           20888.43028204, 20945.51253631, 21009.10266318, 21079.40742639,
           21156.62925036, 21240.96507525, 21332.60540642, 21431.7327481 ,
           21538.52097179, 21653.13346233, 21775.72271998, 21906.42858222,
           22045.3774666 , 22192.68135902, 22348.436745  , 22512.72385296,
    ...
           29752.85399359, 29566.4729752 , 29389.14415891, 29221.09019256,
           29062.51904143, 28913.62417727, 28774.58336394, 28645.55891291,
           28526.69728625, 28418.12887124, 28319.9675907 , 28232.31116402,
           28155.24038831, 28088.81946616, 28033.09558997, 27988.09929761,
           27953.84397622, 27930.32606055, 27917.5253019 , 27915.404695  ,
           27923.91055485, 27942.97286473, 27972.5058195 , 28012.40742056,
           28062.56104774, 28122.83472932, 28193.08289309, 28273.14578651,
           28362.8516702 , 28462.01601292, 28570.44375987, 28687.92938891,
           28814.25817235, 28949.20716149, 29092.54624155, 29244.03894887,
           29403.44380381, 29570.51495296, 29745.00337567, 29926.65779969,
           30115.22531558, 30310.45183827, 30512.08375035, 30719.86693143,
           30933.54857976, 31152.87638438, 31377.59947975, 31607.46824194,
           31842.23411475, 32081.64942323, 32325.4679786 , 32573.44377055,
           32825.33133684, 33080.88560234, 33339.86100795, 33602.01198176,
           33867.09222095, 34134.85477914, 34405.05156725, 34677.43441721,
           34951.75372603, 35227.75967768, 35505.2017655 , 35783.82953864,
           36063.39272352, 36343.64163347, 36624.32747863, 36905.20315001,
           37186.02257686, 37466.54247645, 37746.52179388, 38025.72200997,
           38303.90774086, 38580.84639745, 38856.30861568, 39130.06759317,
           39401.89991915, 39671.58451681, 39938.90274171, 40203.63834787])

  • F
    (Timestamp)
    float64
    2.287e+04 2.281e+04 ... 4.021e+04

    units :

    nT

    description :

    Magnetic field intensity

    array([22867.5503, 22814.5656, 22763.2585, 22713.3703, 22664.7202,
           22617.0636, 22570.1752, 22523.8285, 22477.9168, 22431.9755,
           22386.2023, 22340.3553, 22294.3361, 22248.0264, 22201.3642,
           22154.302 , 22106.7829, 22058.7928, 22010.306 , 21961.3227,
           21911.8445, 21861.9265, 21811.5971, 21760.9304, 21709.9609,
           21658.785 , 21607.4943, 21556.1812, 21504.9585, 21453.9223,
           21403.2124, 21352.9377, 21303.2112, 21254.1668, 21205.9354,
           21158.6381, 21112.4242, 21067.4182, 21023.7546, 20981.5974,
           20941.0772, 20902.3547, 20865.5785, 20830.913 , 20798.4836,
           20768.4778, 20741.0709, 20716.4062, 20694.6578, 20676.0056,
           20660.6319, 20648.7177, 20640.4698, 20636.0623, 20635.6721,
           20639.4882, 20647.7111, 20660.5247, 20678.1302, 20700.7427,
           20728.5752, 20761.8515, 20800.79  , 20845.6239, 20896.5678,
           20953.8253, 21017.6249, 21088.1532, 21165.6066, 21250.1829,
           21342.0523, 21441.4069, 21548.4016, 21663.1958, 21785.954 ,
           21916.7992, 22055.8738, 22203.2925, 22359.1802, 22523.585 ,
           22696.5783, 22878.2053, 23068.4931, 23267.3885, 23474.9191,
           23691.0173, 23915.6834, 24148.8139, 24390.3318, 24640.1631,
           24898.1484, 25164.1641, 25438.1193, 25719.9361, 26009.3985,
           26306.3882, 26610.6352, 26921.9652, 27240.0891, 27564.5525,
    ...
           34918.1532, 34614.2141, 34313.0318, 34014.9108, 33720.193 ,
           33429.1474, 33142.0841, 32859.283 , 32581.112 , 32307.9308,
           32040.1152, 31778.0859, 31522.3264, 31273.2564, 31031.2569,
           30796.66  , 30569.6687, 30350.4844, 30139.266 , 29936.2242,
           29741.52  , 29555.4189, 29378.0553, 29209.599 , 29050.2582,
           28900.2357, 28759.6507, 28628.7647, 28507.836 , 28397.152 ,
           28296.9839, 28207.6563, 28129.5151, 28062.8582, 28007.8342,
           27964.4422, 27932.5326, 27911.8126, 27901.9587, 27902.6501,
           27913.6336, 27934.7316, 27965.7887, 28006.7474, 28057.5472,
           28118.1506, 28188.4974, 28268.4893, 28358.0448, 28457.0458,
           28565.3181, 28682.6517, 28808.8122, 28943.5375, 29086.6526,
           29237.9713, 29397.2705, 29564.2635, 29738.709 , 29920.3739,
           30109.0413, 30304.4248, 30506.2857, 30714.3751, 30928.3665,
           31147.9944, 31372.9382, 31602.9716, 31837.7928, 32077.1767,
           32320.8922, 32568.7267, 32820.418 , 33075.8155, 33334.7625,
           33597.0478, 33862.3655, 34130.4558, 34401.05  , 34673.9007,
           34948.7001, 35225.2151, 35503.1946, 35782.3497, 36062.4684,
           36343.3197, 36624.6573, 36906.1629, 37187.5748, 37468.5924,
           37748.9907, 38028.5159, 38306.9906, 38584.2225, 38860.0354,
           39134.1607, 39406.3744, 39676.4155, 39943.9243, 40208.6581])

  • F_MCO_SHA_2D
    (Timestamp)
    float64
    2.287e+04 2.282e+04 ... 4.021e+04

    units :

    nT

    description :

    Magnetic field intensity, calculated by the MCO_SHA_2D spherical harmonic model

    array([22874.21150876, 22820.9414253 , 22769.3691607 , 22719.23823997,
           22670.30468073, 22622.33792   , 22575.12230553, 22528.4579632 ,
           22482.16181488, 22436.06809821, 22390.02871531, 22343.91359596,
           22297.61076902, 22251.02639378, 22204.08442392, 22156.72599394,
           22108.90939598, 22060.60912872, 22011.81525336, 21962.53302081,
           21912.7813279 , 21862.59277419, 21812.01183649, 21761.09491627,
           21709.90856487, 21658.52935136, 21607.04261608, 21555.54160377,
           21504.12694187, 21452.90578566, 21401.99103716, 21351.50091329,
           21301.55828264, 21252.29017168, 21203.82714295, 21156.30328653,
           21109.85580438, 21064.62434957, 21020.75159148, 20978.38227606,
           20937.66394936, 20898.74619803, 20861.78114204, 20826.92369534,
           20794.33097448, 20764.16308639, 20736.58298782, 20711.75671783,
           20689.85367248, 20671.0465089 , 20655.51157484, 20643.4290656 ,
           20634.98269735, 20630.36018174, 20629.75322598, 20633.35703247,
           20641.37057788, 20653.99640776, 20671.43985026, 20693.90939912,
           20721.61542764, 20754.77037638, 20793.58759432, 20838.28079566,
           20889.06346574, 20946.14753118, 21009.74279054, 21080.05581012,
           21157.28880438, 21241.63851139, 21333.29525952, 21432.44141394,
           21539.25075267, 21653.88661456, 21776.50149756, 21907.2352712 ,
           22046.21440529, 22193.55094038, 22349.34139861, 22513.66600569,
    ...
           29753.73149158, 29567.29857266, 29389.91737916, 29221.81070888,
           29063.18671068, 28914.23907109, 28775.14579608, 28646.06946135,
           28527.1568092 , 28418.53851618, 28320.32879594, 28232.62565262,
           28155.51015437, 28089.04675398, 28033.28286664, 27988.24922   ,
           27953.9593535 , 27930.40981354, 27917.58042056, 27915.43419568,
           27923.91743863, 27942.96007831, 27972.47621988, 28012.36374472,
           28062.5058881 , 28122.77051576, 28193.01188183, 28273.07005607,
           28362.77312409, 28461.93638897, 28570.36464514, 28687.85223957,
           28814.18433744, 28949.13790944, 29092.48278814, 29243.9824857 ,
           29403.39552589, 29570.4760846 , 29744.9751937 , 29926.64165333,
           30115.2226426 , 30310.46417701, 30512.11274767, 30719.91434675,
           30933.61628543, 31152.96636296, 31377.71381857, 31607.60912575,
           31842.40381692, 32081.850296  , 32325.70244284, 32573.71430537,
           32825.640469  , 33081.23589554, 33340.25505164, 33602.45237909,
           33867.58157525, 34135.3956783 , 34405.6465653 , 34678.08601239,
           34952.46433484, 35228.53160452, 35506.03716763, 35784.73038678,
           36064.36075853, 36344.67831995, 36625.43395701, 36906.38018846,
           37187.27052549, 37467.86122592, 37747.91074114, 38027.18003385,
           38305.43319013, 38582.43709407, 38857.96187529, 39131.78026476,
           39403.66844511, 39673.40501375, 39940.77110175, 40205.55035918])

  • B_NEC_CHAOS-Core
    (Timestamp, NEC)
    float64
    2.011e+04 -4.127e+03 ... 3.557e+04

    units :

    nT

    description :

    Magnetic field vector, NEC frame, calculated by the CHAOS-Core spherical harmonic model

    array([[ 20113.22615297,  -4127.05133197, -10082.92688681],
           [ 19824.73224609,  -4162.69719018, -10509.89021079],
           [ 19533.08978642,  -4197.08111622, -10922.34336608],
           ...,
           [ 19127.53601852,  -1776.09964288,  34710.4804633 ],
           [ 18887.30169639,  -1793.46048933,  35144.9752641 ],
           [ 18643.03672966,  -1808.26071658,  35573.86556285]])

  • B_NEC_MCO_SHA_2D
    (Timestamp, NEC)
    float64
    2.011e+04 -4.127e+03 ... 3.557e+04

    units :

    nT

    description :

    Magnetic field vector, NEC frame, calculated by the MCO_SHA_2D spherical harmonic model

    array([[ 20113.62392147,  -4127.46395613, -10081.45456687],
           [ 19825.16184358,  -4163.12754932, -10508.41065152],
           [ 19533.55390549,  -4197.52905375, -10920.86048078],
           ...,
           [ 19129.93319584,  -1775.95410044,  34711.24757684],
           [ 18889.6896485 ,  -1793.32967817,  35145.82179792],
           [ 18645.41210832,  -1808.14125488,  35574.78759174]])

  • B_NEC
    (Timestamp, NEC)
    float64
    2.01e+04 -4.126e+03 ... 3.558e+04

    units :

    nT

    description :

    Magnetic field vector, NEC frame

    array([[ 20103.5246,  -4126.2621, -10086.9888],
           [ 19815.0914,  -4160.9933, -10514.4074],
           [ 19523.4946,  -4195.1968, -10926.9664],
           ...,
           [ 19122.434 ,  -1772.6702,  34719.0503],
           [ 18879.9888,  -1790.0997,  35154.8312],
           [ 18635.3729,  -1805.3596,  35583.7714]])

• Attributes: (3)

  Sources :

  ['SW_OPER_MAGA_LR_1B_20140101T000000_20140101T235959_0505_MDR_MAG_LR', 'SW_OPER_MCO_SHA_2D_20131126T000000_20180101T000000_0401', 'SW_OPER_MCO_SHA_2X_19970101T000000_20220807T235959_0711']

  MagneticModels :

  ["CHAOS-Core = 'CHAOS-Core'(max_degree=20,min_degree=1)", 'MCO_SHA_2D = MCO_SHA_2D(max_degree=20,min_degree=1)']

  AppliedFilters :

  []

Fetch the residuals directly

Adding residuals=True to .set_products() will instead directly evaluate and return all data-model residuals

request = SwarmRequest()
request.set_collection("SW_OPER_MAGA_LR_1B")
request.set_products(
    measurements=["F", "B_NEC"],
    models=["CHAOS-Core", "MCO_SHA_2D"],
    residuals=True,
    sampling_step="PT10S"
)
data = request.get_between(
    start_time = dt.datetime(2014,1,1, 0),
    end_time = dt.datetime(2014,1,1, 1)
)
df = data.as_dataframe(expand=True)
df.head()

	Spacecraft	Radius	Latitude	Longitude	F_res_MCO_SHA_2D	F_res_CHAOS-Core	B_NEC_res_MCO_SHA_2D_N	B_NEC_res_MCO_SHA_2D_E	B_NEC_res_MCO_SHA_2D_C	B_NEC_res_CHAOS-Core_N	B_NEC_res_CHAOS-Core_E	B_NEC_res_CHAOS-Core_C
Timestamp
2014-01-01 00:00:00	A	6878309.22	-1.228938	-14.116674	-6.661209	-6.885946	-10.099321	1.201856	-5.534233	-9.701553	0.789232	-4.061913
2014-01-01 00:00:10	A	6878381.17	-1.862521	-14.131424	-6.375825	-6.605465	-10.070444	2.134249	-5.996748	-9.640846	1.703890	-4.517189
2014-01-01 00:00:20	A	6878452.05	-2.496090	-14.146155	-6.110661	-6.341213	-10.059305	2.332254	-6.105919	-9.595186	1.884316	-4.623034
2014-01-01 00:00:30	A	6878521.87	-3.129644	-14.160861	-5.867940	-6.095612	-10.104972	2.127119	-6.112134	-9.604261	1.661755	-4.630222
2014-01-01 00:00:40	A	6878590.61	-3.763184	-14.175534	-5.584481	-5.805765	-10.194444	1.438234	-6.141904	-9.655758	0.955594	-4.665504

Plot the scalar residuals

… using the pandas method:

ax = df.plot(
    y=["F_res_CHAOS-Core", "F_res_MCO_SHA_2D"],
    figsize=(15,5),
    grid=True
)
ax.set_xlabel("Timestamp")
ax.set_ylabel("[nT]");

… using matplotlib interface

NB: we are doing plt.plot(x, y) with x as df.index (the time-based index of df), and y as df[".."]

plt.figure(figsize=(15,5))
plt.plot(
    df.index,
    df["F_res_CHAOS-Core"],
    label="F_res_CHAOS-Core"
)
plt.plot(
    df.index,
    df["F_res_MCO_SHA_2D"],
    label="F_res_MCO_SHA_2D"
)
plt.xlabel("Timestamp")
plt.ylabel("[nT]")
plt.grid()
plt.legend();

… using matplotlib interface (Object Oriented style)

This is the recommended route for making more complicated figures

fig, ax = plt.subplots(figsize=(15,5))
ax.plot(
    df.index,
    df["F_res_CHAOS-Core"],
    label="F_res_CHAOS-Core"
)
ax.plot(
    df.index,
    df["F_res_MCO_SHA_2D"],
    label="F_res_MCO_SHA_2D"
)
ax.set_xlabel("Timestamp")
ax.set_ylabel("[nT]")
ax.grid()
ax.legend();

Plot the vector components

fig, axes = plt.subplots(nrows=3, ncols=1, figsize=(15,10), sharex=True)
for component, ax in zip("NEC", axes):
    for model_name in ("CHAOS-Core", "MCO_SHA_2D"):
        ax.plot(
            df.index,
            df[f"B_NEC_res_{model_name}_{component}"],
            label=model_name
        )
    ax.set_ylabel(f"{component}\n[nT]")
    ax.legend()
axes[0].set_title("Residuals to models (NEC components)")
axes[2].set_xlabel("Timestamp");

Similar plotting, using the data via xarray instead

xarray provides a more sophisticated data structure that is more suitable for the complex vector data we are accessing, together with nice stuff like unit and other metadata support. Unfortunately due to the extra complexity, this can make it difficult to use right away.

ds = data.as_xarray()
ds

<xarray.Dataset>
Dimensions:               (Timestamp: 360, NEC: 3)
Coordinates:
  * Timestamp             (Timestamp) datetime64[ns] 2014-01-01 ... 2014-01-0...
  * NEC                   (NEC) <U1 'N' 'E' 'C'
Data variables:
    Spacecraft            (Timestamp) object 'A' 'A' 'A' 'A' ... 'A' 'A' 'A' 'A'
    Radius                (Timestamp) float64 6.878e+06 6.878e+06 ... 6.868e+06
    Latitude              (Timestamp) float64 -1.229 -1.863 ... 48.14 48.77
    Longitude             (Timestamp) float64 -14.12 -14.13 ... 153.6 153.6
    B_NEC_res_MCO_SHA_2D  (Timestamp, NEC) float64 -10.1 1.202 ... 2.782 8.984
    B_NEC_res_CHAOS-Core  (Timestamp, NEC) float64 -9.702 0.7892 ... 2.901 9.906
    F_res_MCO_SHA_2D      (Timestamp) float64 -6.661 -6.376 ... 3.153 3.108
    F_res_CHAOS-Core      (Timestamp) float64 -6.886 -6.605 ... 5.022 5.02
Attributes:
    Sources:         ['SW_OPER_MAGA_LR_1B_20140101T000000_20140101T235959_050...
    MagneticModels:  ["CHAOS-Core = 'CHAOS-Core'(max_degree=20,min_degree=1)"...
    AppliedFilters:  []

xarray.Dataset

• Dimensions:
  • Timestamp: 360
  • NEC: 3
• Coordinates: (2)
  • Timestamp
    (Timestamp)
    datetime64[ns]
    2014-01-01 ... 2014-01-01T00:59:50

    description :

    Time stamp

    array(['2014-01-01T00:00:00.000000000', '2014-01-01T00:00:10.000000000',
           '2014-01-01T00:00:20.000000000', ..., '2014-01-01T00:59:30.000000000',
           '2014-01-01T00:59:40.000000000', '2014-01-01T00:59:50.000000000'],
          dtype='datetime64[ns]')

  • NEC
    (NEC)
    <U1
    'N' 'E' 'C'

    units :

    description :

    NEC frame - North, East, Centre (down)

    array(['N', 'E', 'C'], dtype='<U1')

• Data variables: (8)
  • 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', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A'], dtype=object)

  • Radius
    (Timestamp)
    float64
    6.878e+06 6.878e+06 ... 6.868e+06

    units :

    m

    description :

    Position in ITRF – Radius

    array([6878309.22, 6878381.17, 6878452.05, 6878521.87, 6878590.61,
           6878658.27, 6878724.84, 6878790.32, 6878854.7 , 6878917.97,
           6878980.13, 6879041.18, 6879101.12, 6879159.94, 6879217.63,
           6879274.2 , 6879329.64, 6879383.95, 6879437.14, 6879489.19,
           6879540.12, 6879589.91, 6879638.58, 6879686.11, 6879732.52,
           6879777.8 , 6879821.95, 6879864.98, 6879906.89, 6879947.68,
           6879987.36, 6880025.93, 6880063.39, 6880099.74, 6880135.  ,
           6880169.17, 6880202.24, 6880234.24, 6880265.15, 6880295.  ,
           6880323.77, 6880351.49, 6880378.17, 6880403.79, 6880428.38,
           6880451.94, 6880474.48, 6880496.01, 6880516.53, 6880536.06,
           6880554.61, 6880572.17, 6880588.77, 6880604.42, 6880619.11,
           6880632.87, 6880645.71, 6880657.62, 6880668.64, 6880678.75,
           6880687.99, 6880696.35, 6880703.85, 6880710.51, 6880716.32,
           6880721.31, 6880725.49, 6880728.87, 6880731.46, 6880733.28,
           6880734.33, 6880734.64, 6880734.2 , 6880733.05, 6880731.18,
           6880728.61, 6880725.36, 6880721.44, 6880716.86, 6880711.63,
           6880705.77, 6880699.3 , 6880692.22, 6880684.55, 6880676.31,
           6880667.5 , 6880658.13, 6880648.24, 6880637.81, 6880626.88,
           6880615.45, 6880603.53, 6880591.15, 6880578.3 , 6880565.01,
           6880551.29, 6880537.15, 6880522.61, 6880507.67, 6880492.35,
    ...
           6875648.21, 6875587.21, 6875525.54, 6875463.22, 6875400.22,
           6875336.57, 6875272.26, 6875207.3 , 6875141.69, 6875075.42,
           6875008.52, 6874940.98, 6874872.8 , 6874804.  , 6874734.58,
           6874664.55, 6874593.9 , 6874522.66, 6874450.83, 6874378.4 ,
           6874305.41, 6874231.84, 6874157.71, 6874083.03, 6874007.82,
           6873932.07, 6873855.8 , 6873779.02, 6873701.74, 6873623.97,
           6873545.73, 6873467.02, 6873387.86, 6873308.26, 6873228.24,
           6873147.8 , 6873066.96, 6872985.74, 6872904.15, 6872822.2 ,
           6872739.91, 6872657.3 , 6872574.37, 6872491.16, 6872407.66,
           6872323.91, 6872239.9 , 6872155.68, 6872071.23, 6871986.6 ,
           6871901.79, 6871816.82, 6871731.71, 6871646.48, 6871561.14,
           6871475.72, 6871390.23, 6871304.7 , 6871219.15, 6871133.59,
           6871048.04, 6870962.54, 6870877.08, 6870791.71, 6870706.43,
           6870621.28, 6870536.27, 6870451.42, 6870366.75, 6870282.3 ,
           6870198.07, 6870114.09, 6870030.39, 6869946.98, 6869863.89,
           6869781.14, 6869698.75, 6869616.74, 6869535.15, 6869453.98,
           6869373.26, 6869293.01, 6869213.26, 6869134.03, 6869055.34,
           6868977.21, 6868899.67, 6868822.72, 6868746.41, 6868670.75,
           6868595.76, 6868521.47, 6868447.89, 6868375.05, 6868302.96,
           6868231.66, 6868161.16, 6868091.48, 6868022.65, 6867954.68])

  • Latitude
    (Timestamp)
    float64
    -1.229 -1.863 ... 48.14 48.77

    units :

    deg

    description :

    Position in ITRF – Latitude

    array([ -1.2289377,  -1.8625209,  -2.4960899,  -3.1296444,  -3.7631844,
            -4.3967096,  -5.03022  ,  -5.6637153,  -6.2971954,  -6.9306601,
            -7.5641093,  -8.1975429,  -8.8309606,  -9.4643624, -10.097748 ,
           -10.7311174, -11.3644703, -11.9978066, -12.6311262, -13.2644288,
           -13.8977145, -14.5309829, -15.1642339, -15.7974674, -16.4306832,
           -17.0638812, -17.6970611, -18.3302228, -18.9633662, -19.5964911,
           -20.2295973, -20.8626845, -21.4957527, -22.1288017, -22.7618312,
           -23.394841 , -24.027831 , -24.6608009, -25.2937506, -25.9266797,
           -26.5595882, -27.1924756, -27.8253419, -28.4581867, -29.0910098,
           -29.7238109, -30.3565897, -30.9893459, -31.6220792, -32.2547894,
           -32.887476 , -33.5201387, -34.1527771, -34.785391 , -35.4179798,
           -36.0505431, -36.6830806, -37.3155917, -37.9480759, -38.5805328,
           -39.2129618, -39.8453623, -40.4777338, -41.1100756, -41.742387 ,
           -42.3746673, -43.0069157, -43.6391316, -44.271314 , -44.903462 ,
           -45.5355748, -46.1676512, -46.7996903, -47.4316908, -48.0636517,
           -48.6955715, -49.3274489, -49.9592825, -50.5910706, -51.2228117,
           -51.8545038, -52.4861451, -53.1177336, -53.7492669, -54.3807427,
           -55.0121585, -55.6435114, -56.2747986, -56.9060167, -57.5371624,
           -58.1682318, -58.7992209, -59.4301252, -60.06094  , -60.69166  ,
           -61.3222795, -61.9527923, -62.5831914, -63.2134696, -63.8436185,
    ...
           -14.0397279, -13.4057211, -12.7716953, -12.1376508, -11.5035876,
           -10.8695063, -10.2354069,  -9.6012897,  -8.967155 ,  -8.333003 ,
            -7.6988341,  -7.0646483,  -6.4304461,  -5.7962276,  -5.161993 ,
            -4.5277428,  -3.8934769,  -3.2591959,  -2.6248998,  -1.990589 ,
            -1.3562637,  -0.7219241,  -0.0875706,   0.5467966,   1.1811773,
             1.8155712,   2.449978 ,   3.0843975,   3.7188294,   4.3532734,
             4.9877293,   5.6221967,   6.2566755,   6.8911652,   7.5256657,
             8.1601767,   8.7946978,   9.4292288,  10.0637693,  10.6983191,
            11.3328779,  11.9674454,  12.6020212,  13.2366051,  13.8711967,
            14.5057958,  15.1404019,  15.7750149,  16.4096342,  17.0442597,
            17.6788909,  18.3135275,  18.9481691,  19.5828155,  20.2174661,
            20.8521207,  21.4867788,  22.1214401,  22.7561042,  23.3907706,
            24.025439 ,  24.6601089,  25.2947799,  25.9294516,  26.5641235,
            27.1987952,  27.8334661,  28.4681359,  29.102804 ,  29.7374699,
            30.3721332,  31.0067932,  31.6414495,  32.2761015,  32.9107485,
            33.5453901,  34.1800256,  34.8146544,  35.4492757,  36.0838891,
            36.7184936,  37.3530887,  37.9876736,  38.6222475,  39.2568096,
            39.891359 ,  40.525895 ,  41.1604165,  41.7949226,  42.4294124,
            43.0638848,  43.6983388,  44.3327731,  44.9671867,  45.6015782,
            46.2359464,  46.87029  ,  47.5046074,  48.1388972,  48.7731577])

  • Longitude
    (Timestamp)
    float64
    -14.12 -14.13 ... 153.6 153.6

    units :

    deg

    description :

    Position in ITRF – Longitude

    array([-14.1166743, -14.1314243, -14.1461554, -14.1608609, -14.1755342,
           -14.1901685, -14.2047571, -14.2192932, -14.2337699, -14.2481804,
           -14.2625177, -14.2767747, -14.2909444, -14.3050195, -14.3189927,
           -14.3328566, -14.3466037, -14.3602264, -14.3737168, -14.3870671,
           -14.4002693, -14.413315 , -14.4261961, -14.4389038, -14.4514295,
           -14.4637642, -14.4758987, -14.4878238, -14.4995298, -14.5110068,
           -14.5222447, -14.5332332, -14.5439615, -14.5544186, -14.5645933,
           -14.5744739, -14.5840483, -14.5933041, -14.6022285, -14.6108083,
           -14.6190299, -14.6268789, -14.6343409, -14.6414006, -14.6480422,
           -14.6542495, -14.6600055, -14.6652927, -14.6700926, -14.6743864,
           -14.6781542, -14.6813754, -14.6840285, -14.6860912, -14.6875401,
           -14.6883507, -14.6884977, -14.6879543, -14.6866927, -14.6846837,
           -14.6818967, -14.6782997, -14.673859 , -14.6685394, -14.6623038,
           -14.6551131, -14.6469263, -14.6377003, -14.6273894, -14.6159457,
           -14.6033184, -14.589454 , -14.5742958, -14.557784 , -14.539855 ,
           -14.5204415, -14.4994721, -14.476871 , -14.4525575, -14.4264457,
           -14.3984442, -14.3684554, -14.3363752, -14.3020922, -14.2654873,
           -14.2264328, -14.1847917, -14.1404169, -14.0931501, -14.0428207,
           -13.989245 , -13.9322242, -13.8715433, -13.8069696, -13.73825  ,
           -13.6651097, -13.5872493, -13.5043418, -13.41603  , -13.3219222,
    ...
           154.3456616, 154.3325292, 154.3192466, 154.3058219, 154.292263 ,
           154.2785777, 154.2647738, 154.2508587, 154.23684  , 154.222725 ,
           154.208521 , 154.194235 , 154.1798741, 154.1654455, 154.1509559,
           154.1364123, 154.1218215, 154.1071903, 154.0925253, 154.0778333,
           154.063121 , 154.0483951, 154.0336621, 154.0189287, 154.0042016,
           153.9894873, 153.9747926, 153.9601241, 153.9454886, 153.9308927,
           153.9163432, 153.901847 , 153.8874109, 153.8730419, 153.8587469,
           153.8445331, 153.8304075, 153.8163775, 153.8024503, 153.7886334,
           153.7749344, 153.761361 , 153.7479209, 153.7346222, 153.7214729,
           153.7084812, 153.6956557, 153.6830049, 153.6705376, 153.6582629,
           153.6461898, 153.634328 , 153.622687 , 153.6112768, 153.6001076,
           153.5891898, 153.5785343, 153.5681521, 153.5580547, 153.5482538,
           153.5387615, 153.5295904, 153.5207534, 153.5122639, 153.5041357,
           153.496383 , 153.4890208, 153.4820643, 153.4755295, 153.4694328,
           153.4637914, 153.4586231, 153.4539464, 153.4497804, 153.4461454,
           153.443062 , 153.4405521, 153.4386384, 153.4373444, 153.4366949,
           153.4367158, 153.437434 , 153.4388778, 153.4410767, 153.4440619,
           153.4478658, 153.4525225, 153.4580678, 153.4645395, 153.4719772,
           153.4804224, 153.4899192, 153.5005138, 153.5122551, 153.5251946,
           153.5393869, 153.5548894, 153.5717634, 153.5900734, 153.6098879])

  • B_NEC_res_MCO_SHA_2D
    (Timestamp, NEC)
    float64
    -10.1 1.202 -5.534 ... 2.782 8.984

    units :

    nT

    description :

    Magnetic field vector residual, NEC frame, calculated as a difference of the measurement and value of the MCO_SHA_2D spherical harmonic model

    array([[-10.09932147,   1.20185613,  -5.53423313],
           [-10.07044358,   2.13424932,  -5.99674848],
           [-10.05930549,   2.33225375,  -6.10591922],
           ...,
           [ -7.49919584,   3.28390044,   7.80272316],
           [ -9.7008485 ,   3.22997817,   9.00940208],
           [-10.03920832,   2.78165488,   8.98380826]])

  • B_NEC_res_CHAOS-Core
    (Timestamp, NEC)
    float64
    -9.702 0.7892 ... 2.901 9.906

    units :

    nT

    description :

    Magnetic field vector residual, NEC frame, calculated as a difference of the measurement and value of the CHAOS-Core spherical harmonic model

    array([[-9.70155297,  0.78923197, -4.06191319],
           [-9.64084609,  1.70389018, -4.51718921],
           [-9.59518642,  1.88431622, -4.62303392],
           ...,
           [-5.10201852,  3.42944288,  8.5698367 ],
           [-7.31289639,  3.36078933,  9.8559359 ],
           [-7.66382966,  2.90111658,  9.90583715]])

  • F_res_MCO_SHA_2D
    (Timestamp)
    float64
    -6.661 -6.376 ... 3.153 3.108

    units :

    nT

    description :

    Magnetic field intensity residual, calculated as a difference of the measurement and value of the MCO_SHA_2D spherical harmonic model

    array([-6.66120876e+00, -6.37582530e+00, -6.11066070e+00, -5.86793997e+00,
           -5.58448073e+00, -5.27432000e+00, -4.94710553e+00, -4.62946320e+00,
           -4.24501488e+00, -4.09259821e+00, -3.82641531e+00, -3.55829596e+00,
           -3.27466902e+00, -2.99999378e+00, -2.72022392e+00, -2.42399394e+00,
           -2.12649598e+00, -1.81632872e+00, -1.50925336e+00, -1.21032081e+00,
           -9.36827897e-01, -6.66274187e-01, -4.14736489e-01, -1.64516271e-01,
            5.23351295e-02,  2.55648645e-01,  4.51683915e-01,  6.39596232e-01,
            8.31558133e-01,  1.01651434e+00,  1.22136284e+00,  1.43678671e+00,
            1.65291736e+00,  1.87662832e+00,  2.10825705e+00,  2.33481347e+00,
            2.56839562e+00,  2.79385043e+00,  3.00300852e+00,  3.21512394e+00,
            3.41325064e+00,  3.60850197e+00,  3.79735796e+00,  3.98930466e+00,
            4.15262552e+00,  4.31471361e+00,  4.48791218e+00,  4.64948217e+00,
            4.80412752e+00,  4.95909110e+00,  5.12032516e+00,  5.28863440e+00,
            5.48710265e+00,  5.70211826e+00,  5.91887402e+00,  6.13116753e+00,
            6.34052212e+00,  6.52829224e+00,  6.69034974e+00,  6.83330088e+00,
            6.95977236e+00,  7.08112362e+00,  7.20240568e+00,  7.34310434e+00,
            7.50433426e+00,  7.67776882e+00,  7.88210946e+00,  8.09738988e+00,
            8.31779562e+00,  8.54438861e+00,  8.75704048e+00,  8.96548606e+00,
            9.15084733e+00,  9.30918544e+00,  9.45250244e+00,  9.56392880e+00,
            9.65939471e+00,  9.74155962e+00,  9.83880139e+00,  9.91899431e+00,
    ...
           -1.22114916e+01, -1.18796727e+01, -1.18620792e+01, -1.22117089e+01,
           -1.29285107e+01, -1.40033711e+01, -1.54950961e+01, -1.73047614e+01,
           -1.93208092e+01, -2.13865162e+01, -2.33448959e+01, -2.49693526e+01,
           -2.59950544e+01, -2.61885540e+01, -2.54486666e+01, -2.38070200e+01,
           -2.14267535e+01, -1.85972135e+01, -1.56217206e+01, -1.27840957e+01,
           -1.02838386e+01, -8.22847831e+00, -6.68751988e+00, -5.61634472e+00,
           -4.95868810e+00, -4.61991576e+00, -4.51448183e+00, -4.58075607e+00,
           -4.72832409e+00, -4.89058897e+00, -5.04654514e+00, -5.20053957e+00,
           -5.37213744e+00, -5.60040944e+00, -5.83018814e+00, -6.01118570e+00,
           -6.12502589e+00, -6.21258460e+00, -6.26619370e+00, -6.26775333e+00,
           -6.18134260e+00, -6.03937701e+00, -5.82704767e+00, -5.53924675e+00,
           -5.24978543e+00, -4.97196296e+00, -4.77561857e+00, -4.63752575e+00,
           -4.61101692e+00, -4.67359600e+00, -4.81024284e+00, -4.98760537e+00,
           -5.22246900e+00, -5.42039554e+00, -5.49255164e+00, -5.40457909e+00,
           -5.21607525e+00, -4.93987830e+00, -4.59656530e+00, -4.18531239e+00,
           -3.76423484e+00, -3.31650452e+00, -2.84256763e+00, -2.38068678e+00,
           -1.89235853e+00, -1.35861995e+00, -7.76657008e-01, -2.17288460e-01,
            3.04274509e-01,  7.31174078e-01,  1.07995886e+00,  1.33586615e+00,
            1.55740987e+00,  1.78540593e+00,  2.07352471e+00,  2.38043524e+00,
            2.70595489e+00,  3.01048625e+00,  3.15319825e+00,  3.10774082e+00])

  • F_res_CHAOS-Core
    (Timestamp)
    float64
    -6.886 -6.605 -6.341 ... 5.022 5.02

    units :

    nT

    description :

    Magnetic field intensity residual, calculated as a difference of the measurement and value of the CHAOS-Core spherical harmonic model

    array([-6.88594621e+00, -6.60546481e+00, -6.34121337e+00, -6.09561166e+00,
           -5.80576496e+00, -5.48607897e+00, -5.14663564e+00, -4.81454193e+00,
           -4.41392681e+00, -4.24413862e+00, -3.95987088e+00, -3.67340303e+00,
           -3.37155105e+00, -3.07908042e+00, -2.78215522e+00, -2.46951398e+00,
           -2.15634075e+00, -1.83111357e+00, -1.50936550e+00, -1.19582238e+00,
           -9.07374542e-01, -6.21053906e-01, -3.52432308e-01, -8.32955795e-02,
            1.54802987e-01,  3.82146263e-01,  6.05372195e-01,  8.23915580e-01,
            1.05010810e+00,  1.27291635e+00,  1.51911171e+00,  1.77909724e+00,
            2.04257280e+00,  2.31583744e+00,  2.59852654e+00,  2.87684178e+00,
            3.16199350e+00,  3.43789237e+00,  3.69541800e+00,  3.95289401e+00,
            4.19249970e+00,  4.42456300e+00,  4.64489824e+00,  4.86247130e+00,
            5.04521108e+00,  5.22033528e+00,  5.40019725e+00,  5.56225226e+00,
            5.71157416e+00,  5.85593567e+00,  6.00195706e+00,  6.15122221e+00,
            6.32767447e+00,  6.51860893e+00,  6.71013850e+00,  6.89696117e+00,
            7.08145040e+00,  7.24573370e+00,  7.38635714e+00,  7.51048650e+00,
            7.62118386e+00,  7.73011731e+00,  7.84252350e+00,  7.97796073e+00,
            8.13751796e+00,  8.31276369e+00,  8.52223682e+00,  8.74577361e+00,
            8.97734964e+00,  9.21782475e+00,  9.44689358e+00,  9.67415190e+00,
            9.88062821e+00,  1.00623377e+01,  1.02312800e+01,  1.03706178e+01,
            1.04963334e+01,  1.06111410e+01,  1.07434550e+01,  1.08611470e+01,
    ...
           -1.13339936e+01, -1.10540752e+01, -1.10888589e+01, -1.14911926e+01,
           -1.22608414e+01, -1.33884773e+01, -1.49326639e+01, -1.67942129e+01,
           -1.88612863e+01, -2.09768712e+01, -2.29836907e+01, -2.46548640e+01,
           -2.57252883e+01, -2.59612662e+01, -2.52613900e+01, -2.36570976e+01,
           -2.13113762e+01, -1.85134606e+01, -1.55666019e+01, -1.27545950e+01,
           -1.02769548e+01, -8.24126473e+00, -6.71711950e+00, -5.66002056e+00,
           -5.01384774e+00, -4.68412932e+00, -4.58549309e+00, -4.65648651e+00,
           -4.80687020e+00, -4.97021292e+00, -5.12565987e+00, -5.27768891e+00,
           -5.44597235e+00, -5.66966149e+00, -5.89364155e+00, -6.06764887e+00,
           -6.17330381e+00, -6.25145296e+00, -6.29437567e+00, -6.28389969e+00,
           -6.18401558e+00, -6.02703827e+00, -5.79805035e+00, -5.49183143e+00,
           -5.18207976e+00, -4.88198438e+00, -4.66127975e+00, -4.49664194e+00,
           -4.44131475e+00, -4.47272323e+00, -4.57577860e+00, -4.71707055e+00,
           -4.91333684e+00, -5.07010234e+00, -5.09850795e+00, -4.96418176e+00,
           -4.72672095e+00, -4.39897914e+00, -4.00156725e+00, -3.53371721e+00,
           -3.05362603e+00, -2.54457768e+00, -2.00716550e+00, -1.47983864e+00,
           -9.24323524e-01, -3.21933467e-01,  3.29821373e-01,  9.59749995e-01,
            1.55222314e+00,  2.04992355e+00,  2.46890612e+00,  2.79389003e+00,
            3.08285914e+00,  3.37610255e+00,  3.72678432e+00,  4.09310683e+00,
            4.47448085e+00,  4.83098319e+00,  5.02155829e+00,  5.01975213e+00])

• Attributes: (3)

  Sources :

  ['SW_OPER_MAGA_LR_1B_20140101T000000_20140101T235959_0505_MDR_MAG_LR', 'SW_OPER_MCO_SHA_2D_20131126T000000_20180101T000000_0401', 'SW_OPER_MCO_SHA_2X_19970101T000000_20220807T235959_0711']

  MagneticModels :

  ["CHAOS-Core = 'CHAOS-Core'(max_degree=20,min_degree=1)", 'MCO_SHA_2D = MCO_SHA_2D(max_degree=20,min_degree=1)']

  AppliedFilters :

  []

fig, axes = plt.subplots(nrows=3, ncols=1, figsize=(15,10), sharex=True)
for i, ax in enumerate(axes):
    for model_name in ("CHAOS-Core", "MCO_SHA_2D"):
        ax.plot(
            ds["Timestamp"],
            ds[f"B_NEC_res_{model_name}"][:, i],
            label=model_name
        )
    ax.set_ylabel("NEC"[i] + " [nT]")
    ax.legend()
axes[0].set_title("Residuals to models (NEC components)")
axes[2].set_xlabel("Timestamp");

Note that xarray also allows convenient direct plotting like:

ds["B_NEC_res_CHAOS-Core"].plot.line(x="Timestamp");

Access multiple MAG datasets simultaneously

It is possible to fetch data from multiple collections simultaneously. Here we fetch the measurements from Swarm Alpha and Bravo. In the returned data, you can differentiate between them using the “Spacecraft” column.

request = SwarmRequest()
request.set_collection("SW_OPER_MAGA_LR_1B", "SW_OPER_MAGC_LR_1B")
request.set_products(
    measurements=["F", "B_NEC"],
    models=["CHAOS-Core",],
    residuals=True,
    sampling_step="PT10S"
)
data = request.get_between(
    start_time = dt.datetime(2014,1,1, 0),
    end_time = dt.datetime(2014,1,1, 1)
)
df = data.as_dataframe(expand=True)
df.head()

	Radius	Latitude	Longitude	Spacecraft	F_res_CHAOS-Core	B_NEC_res_CHAOS-Core_N	B_NEC_res_CHAOS-Core_E	B_NEC_res_CHAOS-Core_C
Timestamp
2014-01-01 00:00:00	6878309.22	-1.228938	-14.116674	A	-6.885946	-9.701553	0.789232	-4.061913
2014-01-01 00:00:10	6878381.17	-1.862521	-14.131424	A	-6.605465	-9.640846	1.703890	-4.517189
2014-01-01 00:00:20	6878452.05	-2.496090	-14.146155	A	-6.341213	-9.595186	1.884316	-4.623034
2014-01-01 00:00:30	6878521.87	-3.129644	-14.160861	A	-6.095612	-9.604261	1.661755	-4.630222
2014-01-01 00:00:40	6878590.61	-3.763184	-14.175534	A	-5.805765	-9.655758	0.955594	-4.665504

df[df["Spacecraft"] == "A"].head()

	Radius	Latitude	Longitude	Spacecraft	F_res_CHAOS-Core	B_NEC_res_CHAOS-Core_N	B_NEC_res_CHAOS-Core_E	B_NEC_res_CHAOS-Core_C
Timestamp
2014-01-01 00:00:00	6878309.22	-1.228938	-14.116674	A	-6.885946	-9.701553	0.789232	-4.061913
2014-01-01 00:00:10	6878381.17	-1.862521	-14.131424	A	-6.605465	-9.640846	1.703890	-4.517189
2014-01-01 00:00:20	6878452.05	-2.496090	-14.146155	A	-6.341213	-9.595186	1.884316	-4.623034
2014-01-01 00:00:30	6878521.87	-3.129644	-14.160861	A	-6.095612	-9.604261	1.661755	-4.630222
2014-01-01 00:00:40	6878590.61	-3.763184	-14.175534	A	-5.805765	-9.655758	0.955594	-4.665504

df[df["Spacecraft"] == "C"].head()

	Radius	Latitude	Longitude	Spacecraft	F_res_CHAOS-Core	B_NEC_res_CHAOS-Core_N	B_NEC_res_CHAOS-Core_E	B_NEC_res_CHAOS-Core_C
Timestamp
2014-01-01 00:00:00	6877665.99	5.908082	-14.420068	C	-10.283541	-10.312742	1.869603	-0.402592
2014-01-01 00:00:10	6877747.67	5.274386	-14.434576	C	-9.908776	-10.106798	2.009393	-1.033176
2014-01-01 00:00:20	6877828.39	4.640702	-14.449141	C	-9.616648	-10.019548	1.917093	-1.481472
2014-01-01 00:00:30	6877908.15	4.007030	-14.463755	C	-9.401442	-10.127591	1.524755	-2.023153
2014-01-01 00:00:40	6877986.93	3.373371	-14.478412	C	-9.178214	-10.224281	1.057798	-2.427985

… or using xarray

ds = data.as_xarray()
ds.where(ds["Spacecraft"] == "A", drop=True)

<xarray.Dataset>
Dimensions:               (Timestamp: 360, NEC: 3)
Coordinates:
  * Timestamp             (Timestamp) datetime64[ns] 2014-01-01 ... 2014-01-0...
  * NEC                   (NEC) <U1 'N' 'E' 'C'
Data variables:
    Spacecraft            (Timestamp) object 'A' 'A' 'A' 'A' ... 'A' 'A' 'A' 'A'
    Radius                (Timestamp) float64 6.878e+06 6.878e+06 ... 6.868e+06
    Latitude              (Timestamp) float64 -1.229 -1.863 ... 48.14 48.77
    Longitude             (Timestamp) float64 -14.12 -14.13 ... 153.6 153.6
    B_NEC_res_CHAOS-Core  (Timestamp, NEC) float64 -9.702 0.7892 ... 2.901 9.906
    F_res_CHAOS-Core      (Timestamp) float64 -6.886 -6.605 ... 5.022 5.02
Attributes:
    Sources:         ['SW_OPER_MAGA_LR_1B_20140101T000000_20140101T235959_050...
    MagneticModels:  ["CHAOS-Core = 'CHAOS-Core'(max_degree=20,min_degree=1)"]
    AppliedFilters:  []

xarray.Dataset

• Dimensions:
  • Timestamp: 360
  • NEC: 3
• Coordinates: (2)
  • Timestamp
    (Timestamp)
    datetime64[ns]
    2014-01-01 ... 2014-01-01T00:59:50

    description :

    Time stamp

    array(['2014-01-01T00:00:00.000000000', '2014-01-01T00:00:10.000000000',
           '2014-01-01T00:00:20.000000000', ..., '2014-01-01T00:59:30.000000000',
           '2014-01-01T00:59:40.000000000', '2014-01-01T00:59:50.000000000'],
          dtype='datetime64[ns]')

  • NEC
    (NEC)
    <U1
    'N' 'E' 'C'

    units :

    description :

    NEC frame - North, East, Centre (down)

    array(['N', 'E', 'C'], dtype='<U1')

• Data variables: (6)
  • 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', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A',
           'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A', 'A'], dtype=object)

  • Radius
    (Timestamp)
    float64
    6.878e+06 6.878e+06 ... 6.868e+06

    units :

    m

    description :

    Position in ITRF – Radius

    array([6878309.22, 6878381.17, 6878452.05, 6878521.87, 6878590.61,
           6878658.27, 6878724.84, 6878790.32, 6878854.7 , 6878917.97,
           6878980.13, 6879041.18, 6879101.12, 6879159.94, 6879217.63,
           6879274.2 , 6879329.64, 6879383.95, 6879437.14, 6879489.19,
           6879540.12, 6879589.91, 6879638.58, 6879686.11, 6879732.52,
           6879777.8 , 6879821.95, 6879864.98, 6879906.89, 6879947.68,
           6879987.36, 6880025.93, 6880063.39, 6880099.74, 6880135.  ,
           6880169.17, 6880202.24, 6880234.24, 6880265.15, 6880295.  ,
           6880323.77, 6880351.49, 6880378.17, 6880403.79, 6880428.38,
           6880451.94, 6880474.48, 6880496.01, 6880516.53, 6880536.06,
           6880554.61, 6880572.17, 6880588.77, 6880604.42, 6880619.11,
           6880632.87, 6880645.71, 6880657.62, 6880668.64, 6880678.75,
           6880687.99, 6880696.35, 6880703.85, 6880710.51, 6880716.32,
           6880721.31, 6880725.49, 6880728.87, 6880731.46, 6880733.28,
           6880734.33, 6880734.64, 6880734.2 , 6880733.05, 6880731.18,
           6880728.61, 6880725.36, 6880721.44, 6880716.86, 6880711.63,
           6880705.77, 6880699.3 , 6880692.22, 6880684.55, 6880676.31,
           6880667.5 , 6880658.13, 6880648.24, 6880637.81, 6880626.88,
           6880615.45, 6880603.53, 6880591.15, 6880578.3 , 6880565.01,
           6880551.29, 6880537.15, 6880522.61, 6880507.67, 6880492.35,
    ...
           6875648.21, 6875587.21, 6875525.54, 6875463.22, 6875400.22,
           6875336.57, 6875272.26, 6875207.3 , 6875141.69, 6875075.42,
           6875008.52, 6874940.98, 6874872.8 , 6874804.  , 6874734.58,
           6874664.55, 6874593.9 , 6874522.66, 6874450.83, 6874378.4 ,
           6874305.41, 6874231.84, 6874157.71, 6874083.03, 6874007.82,
           6873932.07, 6873855.8 , 6873779.02, 6873701.74, 6873623.97,
           6873545.73, 6873467.02, 6873387.86, 6873308.26, 6873228.24,
           6873147.8 , 6873066.96, 6872985.74, 6872904.15, 6872822.2 ,
           6872739.91, 6872657.3 , 6872574.37, 6872491.16, 6872407.66,
           6872323.91, 6872239.9 , 6872155.68, 6872071.23, 6871986.6 ,
           6871901.79, 6871816.82, 6871731.71, 6871646.48, 6871561.14,
           6871475.72, 6871390.23, 6871304.7 , 6871219.15, 6871133.59,
           6871048.04, 6870962.54, 6870877.08, 6870791.71, 6870706.43,
           6870621.28, 6870536.27, 6870451.42, 6870366.75, 6870282.3 ,
           6870198.07, 6870114.09, 6870030.39, 6869946.98, 6869863.89,
           6869781.14, 6869698.75, 6869616.74, 6869535.15, 6869453.98,
           6869373.26, 6869293.01, 6869213.26, 6869134.03, 6869055.34,
           6868977.21, 6868899.67, 6868822.72, 6868746.41, 6868670.75,
           6868595.76, 6868521.47, 6868447.89, 6868375.05, 6868302.96,
           6868231.66, 6868161.16, 6868091.48, 6868022.65, 6867954.68])

  • Latitude
    (Timestamp)
    float64
    -1.229 -1.863 ... 48.14 48.77

    units :

    deg

    description :

    Position in ITRF – Latitude

    array([ -1.2289377,  -1.8625209,  -2.4960899,  -3.1296444,  -3.7631844,
            -4.3967096,  -5.03022  ,  -5.6637153,  -6.2971954,  -6.9306601,
            -7.5641093,  -8.1975429,  -8.8309606,  -9.4643624, -10.097748 ,
           -10.7311174, -11.3644703, -11.9978066, -12.6311262, -13.2644288,
           -13.8977145, -14.5309829, -15.1642339, -15.7974674, -16.4306832,
           -17.0638812, -17.6970611, -18.3302228, -18.9633662, -19.5964911,
           -20.2295973, -20.8626845, -21.4957527, -22.1288017, -22.7618312,
           -23.394841 , -24.027831 , -24.6608009, -25.2937506, -25.9266797,
           -26.5595882, -27.1924756, -27.8253419, -28.4581867, -29.0910098,
           -29.7238109, -30.3565897, -30.9893459, -31.6220792, -32.2547894,
           -32.887476 , -33.5201387, -34.1527771, -34.785391 , -35.4179798,
           -36.0505431, -36.6830806, -37.3155917, -37.9480759, -38.5805328,
           -39.2129618, -39.8453623, -40.4777338, -41.1100756, -41.742387 ,
           -42.3746673, -43.0069157, -43.6391316, -44.271314 , -44.903462 ,
           -45.5355748, -46.1676512, -46.7996903, -47.4316908, -48.0636517,
           -48.6955715, -49.3274489, -49.9592825, -50.5910706, -51.2228117,
           -51.8545038, -52.4861451, -53.1177336, -53.7492669, -54.3807427,
           -55.0121585, -55.6435114, -56.2747986, -56.9060167, -57.5371624,
           -58.1682318, -58.7992209, -59.4301252, -60.06094  , -60.69166  ,
           -61.3222795, -61.9527923, -62.5831914, -63.2134696, -63.8436185,
    ...
           -14.0397279, -13.4057211, -12.7716953, -12.1376508, -11.5035876,
           -10.8695063, -10.2354069,  -9.6012897,  -8.967155 ,  -8.333003 ,
            -7.6988341,  -7.0646483,  -6.4304461,  -5.7962276,  -5.161993 ,
            -4.5277428,  -3.8934769,  -3.2591959,  -2.6248998,  -1.990589 ,
            -1.3562637,  -0.7219241,  -0.0875706,   0.5467966,   1.1811773,
             1.8155712,   2.449978 ,   3.0843975,   3.7188294,   4.3532734,
             4.9877293,   5.6221967,   6.2566755,   6.8911652,   7.5256657,
             8.1601767,   8.7946978,   9.4292288,  10.0637693,  10.6983191,
            11.3328779,  11.9674454,  12.6020212,  13.2366051,  13.8711967,
            14.5057958,  15.1404019,  15.7750149,  16.4096342,  17.0442597,
            17.6788909,  18.3135275,  18.9481691,  19.5828155,  20.2174661,
            20.8521207,  21.4867788,  22.1214401,  22.7561042,  23.3907706,
            24.025439 ,  24.6601089,  25.2947799,  25.9294516,  26.5641235,
            27.1987952,  27.8334661,  28.4681359,  29.102804 ,  29.7374699,
            30.3721332,  31.0067932,  31.6414495,  32.2761015,  32.9107485,
            33.5453901,  34.1800256,  34.8146544,  35.4492757,  36.0838891,
            36.7184936,  37.3530887,  37.9876736,  38.6222475,  39.2568096,
            39.891359 ,  40.525895 ,  41.1604165,  41.7949226,  42.4294124,
            43.0638848,  43.6983388,  44.3327731,  44.9671867,  45.6015782,
            46.2359464,  46.87029  ,  47.5046074,  48.1388972,  48.7731577])

  • Longitude
    (Timestamp)
    float64
    -14.12 -14.13 ... 153.6 153.6

    units :

    deg

    description :

    Position in ITRF – Longitude

    array([-14.1166743, -14.1314243, -14.1461554, -14.1608609, -14.1755342,
           -14.1901685, -14.2047571, -14.2192932, -14.2337699, -14.2481804,
           -14.2625177, -14.2767747, -14.2909444, -14.3050195, -14.3189927,
           -14.3328566, -14.3466037, -14.3602264, -14.3737168, -14.3870671,
           -14.4002693, -14.413315 , -14.4261961, -14.4389038, -14.4514295,
           -14.4637642, -14.4758987, -14.4878238, -14.4995298, -14.5110068,
           -14.5222447, -14.5332332, -14.5439615, -14.5544186, -14.5645933,
           -14.5744739, -14.5840483, -14.5933041, -14.6022285, -14.6108083,
           -14.6190299, -14.6268789, -14.6343409, -14.6414006, -14.6480422,
           -14.6542495, -14.6600055, -14.6652927, -14.6700926, -14.6743864,
           -14.6781542, -14.6813754, -14.6840285, -14.6860912, -14.6875401,
           -14.6883507, -14.6884977, -14.6879543, -14.6866927, -14.6846837,
           -14.6818967, -14.6782997, -14.673859 , -14.6685394, -14.6623038,
           -14.6551131, -14.6469263, -14.6377003, -14.6273894, -14.6159457,
           -14.6033184, -14.589454 , -14.5742958, -14.557784 , -14.539855 ,
           -14.5204415, -14.4994721, -14.476871 , -14.4525575, -14.4264457,
           -14.3984442, -14.3684554, -14.3363752, -14.3020922, -14.2654873,
           -14.2264328, -14.1847917, -14.1404169, -14.0931501, -14.0428207,
           -13.989245 , -13.9322242, -13.8715433, -13.8069696, -13.73825  ,
           -13.6651097, -13.5872493, -13.5043418, -13.41603  , -13.3219222,
    ...
           154.3456616, 154.3325292, 154.3192466, 154.3058219, 154.292263 ,
           154.2785777, 154.2647738, 154.2508587, 154.23684  , 154.222725 ,
           154.208521 , 154.194235 , 154.1798741, 154.1654455, 154.1509559,
           154.1364123, 154.1218215, 154.1071903, 154.0925253, 154.0778333,
           154.063121 , 154.0483951, 154.0336621, 154.0189287, 154.0042016,
           153.9894873, 153.9747926, 153.9601241, 153.9454886, 153.9308927,
           153.9163432, 153.901847 , 153.8874109, 153.8730419, 153.8587469,
           153.8445331, 153.8304075, 153.8163775, 153.8024503, 153.7886334,
           153.7749344, 153.761361 , 153.7479209, 153.7346222, 153.7214729,
           153.7084812, 153.6956557, 153.6830049, 153.6705376, 153.6582629,
           153.6461898, 153.634328 , 153.622687 , 153.6112768, 153.6001076,
           153.5891898, 153.5785343, 153.5681521, 153.5580547, 153.5482538,
           153.5387615, 153.5295904, 153.5207534, 153.5122639, 153.5041357,
           153.496383 , 153.4890208, 153.4820643, 153.4755295, 153.4694328,
           153.4637914, 153.4586231, 153.4539464, 153.4497804, 153.4461454,
           153.443062 , 153.4405521, 153.4386384, 153.4373444, 153.4366949,
           153.4367158, 153.437434 , 153.4388778, 153.4410767, 153.4440619,
           153.4478658, 153.4525225, 153.4580678, 153.4645395, 153.4719772,
           153.4804224, 153.4899192, 153.5005138, 153.5122551, 153.5251946,
           153.5393869, 153.5548894, 153.5717634, 153.5900734, 153.6098879])

  • B_NEC_res_CHAOS-Core
    (Timestamp, NEC)
    float64
    -9.702 0.7892 ... 2.901 9.906

    units :

    nT

    description :

    Magnetic field vector residual, NEC frame, calculated as a difference of the measurement and value of the CHAOS-Core spherical harmonic model

    array([[-9.70155297,  0.78923197, -4.06191319],
           [-9.64084609,  1.70389018, -4.51718921],
           [-9.59518642,  1.88431622, -4.62303392],
           ...,
           [-5.10201852,  3.42944288,  8.5698367 ],
           [-7.31289639,  3.36078933,  9.8559359 ],
           [-7.66382966,  2.90111658,  9.90583715]])

  • F_res_CHAOS-Core
    (Timestamp)
    float64
    -6.886 -6.605 -6.341 ... 5.022 5.02

    units :

    nT

    description :

    Magnetic field intensity residual, calculated as a difference of the measurement and value of the CHAOS-Core spherical harmonic model

    array([-6.88594621e+00, -6.60546481e+00, -6.34121337e+00, -6.09561166e+00,
           -5.80576496e+00, -5.48607897e+00, -5.14663564e+00, -4.81454193e+00,
           -4.41392681e+00, -4.24413862e+00, -3.95987088e+00, -3.67340303e+00,
           -3.37155105e+00, -3.07908042e+00, -2.78215522e+00, -2.46951398e+00,
           -2.15634075e+00, -1.83111357e+00, -1.50936550e+00, -1.19582238e+00,
           -9.07374542e-01, -6.21053906e-01, -3.52432308e-01, -8.32955795e-02,
            1.54802987e-01,  3.82146263e-01,  6.05372195e-01,  8.23915580e-01,
            1.05010810e+00,  1.27291635e+00,  1.51911171e+00,  1.77909724e+00,
            2.04257280e+00,  2.31583744e+00,  2.59852654e+00,  2.87684178e+00,
            3.16199350e+00,  3.43789237e+00,  3.69541800e+00,  3.95289401e+00,
            4.19249970e+00,  4.42456300e+00,  4.64489824e+00,  4.86247130e+00,
            5.04521108e+00,  5.22033528e+00,  5.40019725e+00,  5.56225226e+00,
            5.71157416e+00,  5.85593567e+00,  6.00195706e+00,  6.15122221e+00,
            6.32767447e+00,  6.51860893e+00,  6.71013850e+00,  6.89696117e+00,
            7.08145040e+00,  7.24573370e+00,  7.38635714e+00,  7.51048650e+00,
            7.62118386e+00,  7.73011731e+00,  7.84252350e+00,  7.97796073e+00,
            8.13751796e+00,  8.31276369e+00,  8.52223682e+00,  8.74577361e+00,
            8.97734964e+00,  9.21782475e+00,  9.44689358e+00,  9.67415190e+00,
            9.88062821e+00,  1.00623377e+01,  1.02312800e+01,  1.03706178e+01,
            1.04963334e+01,  1.06111410e+01,  1.07434550e+01,  1.08611470e+01,
    ...
           -1.13339936e+01, -1.10540752e+01, -1.10888589e+01, -1.14911926e+01,
           -1.22608414e+01, -1.33884773e+01, -1.49326639e+01, -1.67942129e+01,
           -1.88612863e+01, -2.09768712e+01, -2.29836907e+01, -2.46548640e+01,
           -2.57252883e+01, -2.59612662e+01, -2.52613900e+01, -2.36570976e+01,
           -2.13113762e+01, -1.85134606e+01, -1.55666019e+01, -1.27545950e+01,
           -1.02769548e+01, -8.24126473e+00, -6.71711950e+00, -5.66002056e+00,
           -5.01384774e+00, -4.68412932e+00, -4.58549309e+00, -4.65648651e+00,
           -4.80687020e+00, -4.97021292e+00, -5.12565987e+00, -5.27768891e+00,
           -5.44597235e+00, -5.66966149e+00, -5.89364155e+00, -6.06764887e+00,
           -6.17330381e+00, -6.25145296e+00, -6.29437567e+00, -6.28389969e+00,
           -6.18401558e+00, -6.02703827e+00, -5.79805035e+00, -5.49183143e+00,
           -5.18207976e+00, -4.88198438e+00, -4.66127975e+00, -4.49664194e+00,
           -4.44131475e+00, -4.47272323e+00, -4.57577860e+00, -4.71707055e+00,
           -4.91333684e+00, -5.07010234e+00, -5.09850795e+00, -4.96418176e+00,
           -4.72672095e+00, -4.39897914e+00, -4.00156725e+00, -3.53371721e+00,
           -3.05362603e+00, -2.54457768e+00, -2.00716550e+00, -1.47983864e+00,
           -9.24323524e-01, -3.21933467e-01,  3.29821373e-01,  9.59749995e-01,
            1.55222314e+00,  2.04992355e+00,  2.46890612e+00,  2.79389003e+00,
            3.08285914e+00,  3.37610255e+00,  3.72678432e+00,  4.09310683e+00,
            4.47448085e+00,  4.83098319e+00,  5.02155829e+00,  5.01975213e+00])

• Attributes: (3)

  Sources :

  ['SW_OPER_MAGA_LR_1B_20140101T000000_20140101T235959_0505_MDR_MAG_LR', 'SW_OPER_MAGC_LR_1B_20140101T000000_20140101T235959_0505_MDR_MAG_LR', 'SW_OPER_MCO_SHA_2X_19970101T000000_20220807T235959_0711']

  MagneticModels :

  ["CHAOS-Core = 'CHAOS-Core'(max_degree=20,min_degree=1)"]

  AppliedFilters :

  []

previous

Finding conjunctions

next

MAGxHR_1B (Magnetic field 50Hz)