COSMIC RAY FLUXES CALCULATION BASED ON PARMA

Calculates the fluxes (/cm²/s) of cosmic (secondary) particles as a function of altitude (angular and energy integrated), based on the PARMA program (Fortran 90). An illustrative example is included for electrons with min energy thresholds set at 0.3 and 1 MeV, and altitudes between 10 and 20 km.

• Requirements: a Fortran 90 compiler (like `gfortran`), and optionally `python`.
• Open a terminal and run `make` command to compile the Fortran 90 code and produce executable `electron_fluxes`
• the executable `electron_fluxes` takes 5 input arguments : `particle ID (0:neutron, 1-28:H-Ni, 29-30:muon+-, 31:e-, 32:e+, 33:photon)`, `min energy threshold (MeV)`, `altitude (km)`, `latitude (deg)` and `longitude`.
• the fortran code just processes one set of parameters at a time.
• the python script `run_on_grid.py` runs the compiled executable 
  `electron_fluxes` for a range of input parameters, outputs the results in `.csv` file, and plots results in a `.png` file.
• example results are showed in `results.csv` and `flux_vs_altitude_electron.png`

* use instructions in folder `run_using_docker` to compile and run the fortran code in any platform using docker.