CERN homepage INFN homepage
FLUKA: 16.7} The various SDUM options available with command SPECSOUR Previous Index Next

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-IONF: 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 by the GCRIONF code 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. 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