This page was generated by nbsphinx from docs/examples/solo_epd.ipynb.
Interactive online version:

Solar Orbiter HET data

Only for Google Colab users:

%pip install --upgrade ipympl speasy

try:
    from google.colab import output
    output.enable_custom_widget_manager()
except:
    print("Not running inside Google Collab")

For all users:

import speasy as spz
%matplotlib widget
# Use this instead if you are not using jupyterlab yet
#%matplotlib notebook

import matplotlib
import matplotlib.pyplot as plt

import pandas as pd
import numpy as np

start = "2020-11-28T00:00:00"
stop = "2020-12-05T00:00:00"

Solar Orbiter electron flux

solo_eflux:spz.SpeasyVariable = spz.get_data("amda/solo_het_omni_eflux", start, stop)

print(np.unique(np.diff(solo_eflux.time), return_counts=True))

(array([1000000000, 1001000000], dtype='timedelta64[ns]'), array([603798,   1000]))

eflux_df = solo_eflux.to_dataframe()
# resample to 1s
eflux_df = eflux_df.resample("1s").ffill()
eflux_df.rolling(600).mean().plot()
plt.yscale("log")

Solar Orbiter proton flux

solo_hflux:spz.SpeasyVariable = spz.get_data("amda/solo_het_omni_hhcflux", start, stop)

print(np.unique(np.diff(solo_eflux.time), return_counts=True))

(array([1000000000, 1001000000], dtype='timedelta64[ns]'), array([603798,   1000]))

hflux_df = solo_hflux.to_dataframe()
# resample to 1s
hflux_df = hflux_df.resample("1s").ffill()
hflux_df.rolling(600).mean().plot()
plt.yscale("log")

plt.figure()
pas_omni:spz.SpeasyVariable = spz.get_data("amda/solo_het_omni_hflux", start, stop)
pas_omni.plot(cmap='rainbow', vmin=1e-3, vmax=1e2, edgecolors="face")
plt.tight_layout()