GitHub

Wind

Overview

The Wind data loader provides access to energetic particle data from the 3DP (3-Dimensional Plasma and Energetic Particle Investigation) instrument aboard the Wind spacecraft.

References

Supported Instruments

3DP (3-Dimensional Plasma and Energetic Particle Investigation)

Electron Datasets

  • SFSP: Suprathermal electron omnidirectional fluxes (27-520 keV)
  • SFPD: Suprathermal electron energy-angle distributions (27-520 keV)
  • ELSP: Low energy electron omnidirectional (7-300 keV)
  • ELPD: Low energy electron energy-angle (7-300 keV)
  • EHSP: High energy electron omnidirectional (100-3200 keV)
  • EHPD: High energy electron energy-angle (100-3200 keV)

Proton Datasets

  • SOSP: Suprathermal proton omnidirectional fluxes (70 keV - 6.8 MeV)
  • SOPD: Suprathermal proton energy-angle distributions (70 keV - 6.8 MeV)
  • PLSP: Low energy proton omnidirectional (30 keV - 2 MeV)
  • PLPD: Low energy proton energy-angle (30 keV - 2 MeV)
  • PHSP: High energy proton omnidirectional (1-30 MeV)
  • PHPD: High energy proton energy-angle (1-30 MeV)

Basic Usage

List available Wind datasets:

using SolarEnergeticParticle

datasets = get_datasets(:WI, "3DP")
15-element Vector{String}:
 "WI_PM_3DP"
 "WI_PLSP_3DP"
 "WI_SOSP_3DP"
 "WI_SOPD_3DP"
 "WI_EM_3DP"
 "WI_ELM2_3DP"
 "WI_ELSP_3DP"
 "WI_ELPD_3DP"
 "WI_EHSP_3DP"
 "WI_EHPD_3DP"
 "WI_SFSP_3DP"
 "WI_SFPD_3DP"
 "WI_K0_3DP"
 "WI_EMFITS_E0_3DP"
 "WIND_3DP_ECHSFITS_E0-YR"

Load Wind 3DP data:

# Suprathermal electron data
dataset = "WI_SFSP_3DP"
tmin = "2021/10/28T06"
tmax = "2021/10/29T12"
data = get_data(dataset, tmin, tmax; verbose=true)
(FLUX = Float32[0.05040033 0.00024625 … 1.5791571f-6 1.1752819f-6; 0.05002167 0.00022134681 … 1.5791571f-6 1.1752819f-6; … ; 0.0704314 0.019753063 … 0.0005144044 0.000165026; 0.07213554 0.01896915 … 0.0005243757 0.00017532954], FLUX_STACKED = Float32[0.05040033 0.00024625 … 1.5791571f-6 1.1752819f-6; 0.05002167 0.00022134681 … 1.5791571f-6 1.1752819f-6; … ; 0.0704314 0.019753063 … 0.0005144044 0.000165026; 0.07213554 0.01896915 … 0.0005243757 0.00017532954], FLUX_SIMPLE = Float32[0.05040033 0.00024625 … 1.5791571f-6 1.1752819f-6; 0.05002167 0.00022134681 … 1.5791571f-6 1.1752819f-6; … ; 0.0704314 0.019753063 … 0.0005144044 0.000165026; 0.07213554 0.01896915 … 0.0005243757 0.00017532954], ENERGY_filled = Float32[26994.957 40138.668 … 309508.94 516823.4; 26994.957 40138.668 … 309508.94 516823.4; … ; 26994.957 40138.668 … 309508.94 516823.4; 26994.957 40138.668 … 309508.94 516823.4], TIME = [1.6354008019574304e9; 1.6354008144877086e9; … ; 1.6355087760780115e9; 1.6355087886082902e9;;])
using SpacePhysicsMakie, CairoMakie

tvars2plot = replace!(data.FLUX, 0 => eps())
f = tplot(tvars2plot; plottype=Stairs)
ylims!.(f.axes, 5e-6, 8e-2)
f
Example block output