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
[3]:
start = "2020-11-28T00:00:00"
stop = "2020-12-05T00:00:00"
Solar Orbiter electron flux
[7]:
solo_eflux:spz.SpeasyVariable = spz.get_data("amda/solo_het_omni_eflux", start, stop)
[8]:
print(np.unique(np.diff(solo_eflux.time), return_counts=True))
(array([1000000000, 1001000000], dtype='timedelta64[ns]'), array([603798, 1000]))
[9]:
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
[14]:
solo_hflux:spz.SpeasyVariable = spz.get_data("amda/solo_het_omni_hhcflux", start, stop)
[15]:
print(np.unique(np.diff(solo_eflux.time), return_counts=True))
(array([1000000000, 1001000000], dtype='timedelta64[ns]'), array([603798, 1000]))
[16]:
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")
[17]:
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()