This software computes omnidirectional, integral, or differential (in energy) fluxes of electrons and protons trapped in Earth's radiation belts for specified energy, L-value, and B/B0 range. The fluxes are obtained by an interpolation procedure from the experimental data. The program allows users to specify input parameters and options from file and generates flux tables that stored to file.
This software computes omnidirectional, integral, or differential (in energy) fluxes of electrons and protons trapped in Earth's radiation belts for specified energy, L-value, and B/B0 range. The fluxes are obtained by an interpolation procedure from the experimental data. The program allows the user to specify input parameters and options from file and to save the generated flux tables to disk.
The program will first request the names of the input and output files. The working directory must not contain files with the name of the output file.
Input file should contain lines of following format:
L B/B0 MODEL FLUXTYPE E1 E2
Example:
2.15 2.5 2 1 0.5 0.7
Where: L - McIlwain parameter,
B/B0 - fraction of magnetic field strength in the position of observer and in the position of field line top,
MODEL - one of following integer numbers, flux model:
- AP8MAX, protons, solar maximum
- AP8MIN, protons, solar minimum
- AE4MAX, electrons, solar maximum
- AE4MIN, electrons, solar minimum
- AEI7HI, electrons, high estimate
- AEI7LO, electrons, low estimate
- AE8MAX, electrons, solar maximum
- AE8MIN, electrons, solar minimum
FLUXTYPE - one of following integer numbers, selection of flux type:
- integral
- differential
E1 - for integral flux - threshold of energy (MeV), to calculate Flux(>E1),for differential - lowest threshold of energy (MeV) to calculate Flux(E1<E<E2)
E2 - for integral flux program doesn't take this number into account for differential - highest threshold of energy (MeV) to calculate Flux(E1<E<E2)
Format of the resulting file:
L B/B0 MODEL FLUXTYPE E1 E2 FLUX1
2.1500001 2.5000000 2.0000000 1.0000000 0.5000000 0.7000000 1.2473444E+006
Where L, B/B0, MODEL, FLUXTYPE, E1, E2 are the same as in the input file, but integer numbers obtain a decimal part.
The last column FLUX1 means flux calculated using model of selected name.
Integral flux units are particles/(cm*cm*sec),differential flux units are particles/(cm*cm*sec*MeV)