calculation of the rigidity cutoff in selected number of points is improved

the model is initialized as in 2022-11-05--Rigidity-cutoff--SEP-transport-in-geospace--and--rigidity-cutoff

Setting the SEP model:

                gunzip data/T05/2017_OMNI_5m_with_TS05_variables.dat.gz

    Transport of SEPs in geospace:

                ./Config.pl -application=earth-geospace_boundary_injection_sep_t96
                ./Config.pl -spice-path=/nobackupp17/vtenishe/SPICE/cspice -spice-kernels=/nobackupp17/vtenishe/SPICE/Kernels
                cp input/geospace-rigidity.in .
                cp input/sep-spectrum-09102017-event.dat .
                ./Config.pl -input=geospace-rigidity.in

                make -j

    Calculation of the rigidity cutoff:

                ./Config.pl -application=earth-cutoff-rigidity_t05
                make -j

IMF is in the background geomagnetic field model

Setting up the SEPs in geospace model is fixed. The model is initialized the same way as in case of 2022-10-25--Rigidity-cutoff--SEP-transport-in-geospace--and--rigidity-cutoff

            git checkout earth_cutoff_rigidity_model

            gunzip data/T05/2017_OMNI_5m_with_TS05_variables.dat.gz

Transport of SEPs in geospace:

            ./Config.pl -application=earth-geospace_boundary_injection_sep_t96
            ./Config.pl -spice-path=/nobackupp17/vtenishe/SPICE/cspice -spice-kernels=/nobackupp17/vtenishe/SPICE/Kernels
            cp input/geospace-rigidity.in .
            cp input/sep-spectrum-09102017-event.dat .
            ./Config.pl -input=geospace-rigidity.in

            make -j

Calculation of the rigidity cutoff:

            ./Config.pl -application=earth-cutoff-rigidity_t05
            make -j

it seems that the model of SEP transport in geospace does something meaninful. The fields are calculated with T05 and recalculated during a model run. T05 is set with the parameters from data/T05/2017_OMNI_5m_with_TS05_variables.dat

Setting the model:
        git checkout earth_cutoff_rigidity_model

        gunzip data/T05/2017_OMNI_5m_with_TS05_variables.dat.gz

        ./Config.pl -application=earth-geospace_boundary_injection_sep_t96
        ./Config.pl -spice-path=/nobackupp17/vtenishe/SPICE/cspice -spice-kernels=/nobackupp17/vtenishe/SPICE/Kernels
        cp input/geospace-rigidity.in .
        ./Config.pl -input=geospace-rigidity.in
        make -j

a test of the ridigity cutoff calculates as a surface os a sphere. Sphere vs individual poins are defined in input/geospace-rigidity.in
Initialization of the model:
    git checkout earth_cutoff_rigidity_model
    ./Config.pl -application=earth-cutoff-rigidity_t96
    ./Config.pl -spice-path=/nobackupp17/vtenishe/SPICE/cspice -spice-kernels=/nobackupp17/vtenishe/SPICE/Kernels
    cp input/geospace-rigidity.in .
    ./Config.pl -input=geospace-rigidity.in
    make -j

2022-10-14--Rigidity-cutoff--Forward-model--SW-injected-in-geospace rigidity cutoff mode: added a capability for swittching to the forward modeling. That can be used for modeling transient effects (capturing SEPs), etc.

      Example: injects SW in the domain
      git checkout earth_cutoff_rigidity_model
      ./Config.pl -application=earth-geospace_boundary_injection_sw_t96
      ./Config.pl -spice-path=/nobackupp17/vtenishe/SPICE/cspice -spice-kernels=/nobackupp17/vtenishe/SPICE/Kernels
      cp input/geospace-rigidity.in .
      ./Config.pl -input=geospace-rigidity.in
      cp input/points.dat .
      make -j

     Most of the settings are as in 2022-04-20--rigidity-cutoff....

This version of AMPS was used in modeling SEPs in the heliosphere that has been completed for CCMC 2022 meeting. The SWMF PARAM.in is in srcSEP. The AMPS' input file is input/test/sep_focused_transport-ih--fl-He__diffusion.input

test test_SEP-SC-IH--FL--FL-attachment--Droge_diffusion runs. PARAM.in that is used for the test: SWMF_data/PT/AMPS/data/input/SEP/SC-IH--FL--PT/PARAM.in (commit id e23ae0705587579477452847abd54de73c1e11ce)

example of modeling SEPs moving along magnetic field lines extracted from M-FLAMPA
to run test:
from AMPS directory:
1. make test_SEP-SC-IH--FL--FL-attachment--Droge_diffusion
2. cp data/input/SEP/SC-IH--FL--PT/PARAM_locate_shock.in SWMF-IH/SWMF/run_test/PARAM.in

data/input/SEP/SC-IH--FL--PT/PARAM_locate_shock.in: commit id 6b81fe50a5d0af6b165ef36f70e9aff9bda0283d of SWMF_data

implemented capability todefine the range of the injected particle velocities in terms of either energy or rigidity range (in input/geospace-rigidity.in)

added: post-compile input file
example of the model run: input - test/sep_parker_spiral__field_line.input, post-compile input - test/sep_parker_spiral__field_line.in
important: ./Config.pl -spice-path=nospice
important: a field line completely belongs to a given cpu

SEP model: the example is a stand-alone model with Parker spiral and T96
Input file: test/sep_parker_spiral__field_line.input
Mover: particle simply transported along the magnetic field line without any change of its velocity; no scattering implemented
Important: runs only -np=1. Like in the SWMF SC/IH/OH/AMPS runs, a field line entirelly belongs to a particular CPU. If run with -np!=1, the code does not work. Need to be fixed in the future.

Important: when confic the code need to execute ./Config.pl -spice-path=nospice

1. sampling can be turned off. The results in strange valued of the cutoff rigidity
2. defining of the individual sampling points is fixed

calculation of the rigidity for a set of points as well as for a pherical map works (tested with T96)

added a capability of using a post-compile input file to set the parameters of T96, T05, and location of the sphere where the cutoff is calculated.

Example:
./Config.pl -application=earth-cutoff-rigidity_t96
./Config.pl -spice-path=/nobackupp17/vtenishe/SPICE/cspice -spice-kernels=/nobackupp17/vtenishe/SPICE/Kernels
cp input/geospace-rigidity.in
./Config.pl -input=geospace-rigidity.in
make -j

Geospace rigidity cutoff model: the models works now with T96 (earth-cutoff-rigidity_t96) and T05 (earth-cutoff-rigidity_t05) background magnetic field models. In both cases rigidity cutoff plots are produced. Need to verify the calculated values.

Rigidity cutoff model:
1. earth.input.pl is adopted from the madter branch
2. the model for T96 is set in earth-cutoff-rigidity_t96.input