16.7} The various SDUM options available with command SPECSOUR

    For 
SDUM
 = GCR-ALLF: All-nucleon flux

    All-nucleon flux as explained in 16.1.2}.
    Three different options (average, maximum and minimum flux) are available.
    The program reads fluxes from a file named "allnucok.dat" in which are given
    the total energy (GeV), the fluxes (E.dN/dE) and the neutron/proton ratios.
    It is possible to give an energy interval and to choose a starting radius
    (radius of the emission sphere in case of spherical geometry) or starting
    height (the emission height in case of flat geometry). It is possible to
    activate the vertical geomagnetic cutoff and to give the cutoff value at
    the central latitude, otherwise the geomagnetic cutoff will be not taken
    into account. Ions are treated like separate nucleons, or as alphas and
    protons.

    
WHAT(1)
 = 1: central value
            = 2: minimum value
            = 3: maximum value

    
WHAT(2)
 >= 0: starting radius (cm)
             < 0: starting height (cm)

    
WHAT(3)
 = Minimum energy

    
WHAT(4)
 = Maximum energy

    
WHAT(5)
 = Spectral index for sampling (below transition energy)

    
WHAT(6)
 = Transition energy for sampling (above it, sample from 1/E)

    
SDUM
 : not used

    Continuation card:

    
WHAT(1)
 = 0: no geomagnetic cutoff
            = 1: geomagnetic cutoff is requested
            = 2: the vertical geomagnetic cutoff is read as 
WHAT(2)

    
WHAT(2)
 = vertical geomagnetic cutoff at central latitude for 
WHAT(1)
 = 2,
              no meaning otherwise

    
WHAT(3)
-
WHAT(5)
: no meaning

    
WHAT(6)
 =< 0: nucleons are transported separately
             > 0: transport as many alphas as can be built by neutrons, and the
                  remaining protons

    
SDUM
    = "&" in any position in column 71 to 78 (or in the last field if
              free format is used)


    For 
SDUM
 = GCR-AMS: AMS All-particle flux

    All-particle flux (ion flux), as explained in 16.1.1}.
    The particle composition of the flux can be modified by choosing the minimum
    and maximum atomic number (1 =< Z =< 28). The spectrum components have been
    produced applying the solar modulation with the modified Badhwar code 
    [Bad96] for various modulation parameters to the LIS spectra based
    on AMS02 data of ref. [DTL22] and written on '.spc' files (Z+A+<PhiMV>+.spc). 
    It is possible to give an energy interval and to choose a starting radius 
    (radius of the emission sphere in case of spherical geometry) or starting 
    height (the emission height in case of flat geometry). It is possible to 
    activate the geomagnetic cutoff (
WHAT(7)
 in SPECSOUR) and to input 
    optionally the vertical cutoff value at the central latitude.
    Ions are treated like real ions or can be split.
    The optimized value for spectral index for sampling (below transition
    energy) is gamma = 1.75 (
WHAT(5)
). Above transition energy, the spectrum
    will be assumed to have a 1/E shape. Proper weights are then applied
    to account for the actual spectral indeces. This procedure allows to
    sample also the high energy tails of the primary spectra. It must be
    noted that this is a form of source spectrum biasing and therefore no
    analogue runs are possible.


    For 
SDUM
 = GCR-IONF: OLD All-particle flux

    All-particle flux (ion flux), as explained in 16.1.2}.
    The particle composition of the flux can be modified by choosing the minimum
    and maximum atomic number (1 =< Z =< 28). The spectrum components have been
    produced by the modified Badhwar code [Bad96] for various modulation
    parameters and written on '.spc' files (Z+<PhiMV>+.spc). It is possible to
    give an energy interval and to choose a starting radius (radius of the
    emission sphere in case of spherical geometry) or starting height (the
    emission height in case of flat geometry). It is possible to activate the
    geomagnetic cutoff (
WHAT(7)
 in SPECSOUR) and to input optionally the
    vertical cutoff value at the central latitude.
    Ions are treated like real ions or can be splitted.
    The optimized value for spectral index for sampling (below transition
    energy) is gamma = 1.75 (
WHAT(5)
). Above transition energy, the spectrum
    will be assumed to have a 1/E shape. Proper weights are then applied
    to account for the actual spectral indeces. This procedure allows to
    sample also the high energy tails of the primary spectra. It must be
    noted that this is a form of source spectrum biasing and therefore no
    analogue runs are possible.


    For 
SDUM
 = SPE-SPEC, SPE-2003 or SPE-2005, the source is a Solar
               Particle Event.

    The input parameters are the same for the three 
SDUM
 values.

    For 
SDUM
 = SPE-SPEC or SPE-2005, the spectrum is read from a
        file sep20jan2005.spc

    For 
SDUM
 = SPE-2003, the spectrum is read from a file sep28oct2003.spc

    WHAT (1) = Z_max + 100 * Z_min (Z_min = 1 if none is defined)

    WHAT (2) = Starting radius (cm)

    WHAT (3) = Minimum energy

    WHAT (4) = Maximum energy
               If maximum and minimum energy differ by less than 5% then a fixed
               energy (= Maximum energy) is sampled

    WHAT (5) = Spectral index for sampling (below transition energy)

    WHAT (6) = Transition energy for sampling (above it, sample from 1/E)

    Continuation card:

    
WHAT(1)
 = 0: no geomagnetic cutoff
            = 1: geomagnetic cutoff is requested
            = 2: the vertical geomagnetic cutoff is read from 
WHAT(2)

    
WHAT(2)
 = vertical geomagnetic cutoff at central latitude for 
WHAT(1)
 = 2,
              no meaning otherwise

    
WHAT(3)
 = number of energy point in the spectra
              
Default
: 50

    
WHAT(4)
 = if > 0 vertical run

    
WHAT(5)
 = if > 0 probabilities 1 / (2 x Z) are used for the various ions
              (1 for Z = 1)

    
WHAT(6)
 =< 0: ions are split
             > 0: ions are treated like real ions

    
SDUM
    = "&" in any position in column 71 to 78 (or in the last field if
              free format is used)

Previous Index Next