GitHub

Solar and Heliospheric Observatory (SOHO)

Overview

The SOHO data loader provides access to energetic particle data from multiple instruments aboard the Solar and Heliospheric Observatory spacecraft.

References

Supported Instruments

CELIAS (Charge, Element, and Isotope Analysis System)

  • SEM: Solar Extreme-ultraviolet Monitor
  • Energy range: 0.08-2.0 MeV (protons), 0.32-8.0 MeV/nucleon (He4), 0.64-16.0 MeV/nucleon (CNO)

COSTEP-EPHIN (Comprehensive Suprathermal and Energetic Particle Analyzer)

  • EPHIN: Electron Proton Helium Instrument
  • Energy range: 4.3-53.0 MeV (protons), 4.3-53.0 MeV/nucleon (helium), 0.25-3.0 MeV (electrons)

ERNE (Energetic and Relativistic Nuclei and Electron experiment)

  • LED: Low Energy Detector (1.3-13.0 MeV)
  • HED: High Energy Detector (13.0-100.0 MeV)

Basic Usage

List available SOHO datasets:

using SolarEnergeticParticle

datasets = get_datasets(:SOHO)
16-element Vector{String}:
 "SOHO_CELIAS-SEM_15S"
 "SOHO_CELIAS-SEM_1DAY"
 "SOHO_CELIAS-PM_30S"
 "SOHO_CELIAS-PM_5MIN"
 "SOHO_COSTEP-EPHIN_L3I-1MIN"
 "SOHO_COSTEP-EPHIN_L3I-5MIN"
 "SOHO_COSTEP-EPHIN_L3I-10MIN"
 "SOHO_COSTEP-EPHIN_L3I-30MIN"
 "SOHO_COSTEP-EPHIN_L3I-1HR"
 "SOHO_COSTEP-EPHIN_L3I-1DAY"
 "SOHO_ERNE-LED_L2-1MIN"
 "SOHO_ERNE-HED_L2-1MIN"
 "SOHO_ERNE_HEAVY-ION-1HR"
 "SOHO_ERNE_HEAVY-ION-5MIN"
 "SOHO_HELIO1DAY_POSITION"
 "SOHO_HELIO1HR_POSITION"

Load SOHO data:

# ERNE High Energy Detector data
dataset = "SOHO_ERNE-HED_L2-1MIN"
tmin = "2021/10/28T06"
tmax = "2021/10/29T12"
data = get_data(dataset, tmin, tmax; verbose=true)
(est = UInt32[0x00004080; 0x00004080; … ; 0x00004080; 0x00004080;;], PH = Float32[0.00022 0.0 … 0.000214 0.000208; 0.0 0.0 … 0.000218 0.000141; … ; 0.453 0.385 … 0.00502 0.000626; 0.466 0.418 … 0.00573 0.00117], AH = Float32[0.0 0.0 … 6.05f-5 0.0; 0.0 0.0 … 5.35f-5 9.02f-5; … ; 0.00106 0.000587 … 9.76f-5 0.000105; 0.00079 0.000864 … 0.000186 0.0], PHC = Int32[1 0 … 6 3; 0 0 … 6 2; … ; 1597 2000 … 109 7; 1637 2165 … 124 13], AHC = Int32[0 0 … 1 0; 0 0 … 1 2; … ; 4 3 … 3 2; 3 6 … 3 0])

Advanced Usage

using SpacePhysicsMakie, CairoMakie

begin
    tvars2plot = map([data.AH, data.PH]) do x
        replace!(select_channel(x, 1:3:10), 0 => eps())
    end
    f = tplot(tvars2plot; plottype=Stairs)
    ylims!.(f.axes, 5e-5, 8e-1)
    f
end
Example block output