trend

- 4

Trapped Radiation ENvironment model Development

UA636

[ Name | Synopsis | Arguments | Description | Diagnostics | Common Blocks | Dependencies | See Also | Reported Bugs | Examples ]

Name

    Subroutine UA636 - Evaluate a weighted atmospheric mass

Synopsis

      record /zgeo/ mpos 
      record /zvec/mb
      integer*4     ktrl, ifail
      real*8        eng, rmass

      CALL UA636 (ktrl, eng, mpos, rmass, ifail)

Arguments

    ktrl
    Control parameter [in]
    eng
    Energy [in]
    mpos
    Geographic location [in]
    mb
    Magnetic field vector (not used) [in]
    rmass
    Weighted atmospheric mass [out]
    ifail
    Error flag (see Diagnostics) [out]

Description

    The subroutine UA636 provides an easy way to average atmospheric density numbers. The argument ktrl is used to control the kind of average (see Table below). When the argument ktrl is different than zero, the weighting includes particle cross sections evaluated by subroutine UA637 or UA638. The cross sections compiled by Pierrard (see subroutine UA637) are used when ktrl is set to 1 or 2, while the Hassit cross sections (see subroutine UA636) are used when ktrl is set to 11 or 12. When the argument ktrl is equal to 0, the density numbers are weighted only by the atomic mass of each component.

    ktrl Weight for the density numbers
    0Atomic mass
    1,11(Atomic mass) x (Proton cross section)
    2,12Proton cross section
    1,2Subroutine UA637 is used
    11,12Subroutine UA638 is used (Hassit)

    Note that the units of the output argument rmass are function of the paramater ktrl. The argument eng is then used to indicate the particle energy. The argument mpos indicates where the atmospheric densities have to be evaluated. The subroutine UA636 can be used as argument of subroutine UD320. Be aware that before any call to the subroutine UA636, an atmospheric model has to be selected with the help of subroutine UA610.

History

    In version 2.00 and earlier, the sysnopsis of subroutine UA636 does not include the argument mb. This argument has been added from version 2.01 to follow a modification of synopsys in subroutines UD320 and UD321.

Diagnostics

    -63601, atmospheric model incompatible with ktrl parameter
    -63602, illegal ktrl parameter

Common Blocks

    UC150, atmospheric model description

Dependencies

Called by

    None

Calls

    UA630, evaluate the atmospheric densities
    UA634, compute mass density from number density
    UA635, evaluate the Debye length
    UA637, proton cross sections
    UA638, Hassitt proton cross sections

See also

    UD320, evaluate a time integral over a drift shell
    UA610, select an atmospheric, ionospheric and/or plasmaspheric model

Reported Bugs

    None

Examples

    None