trend

- 4

Trapped Radiation ENvironment model Development

UT990

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

Name

    Subroutine UT990 - initialize the UNILIB libray

Synopsis

      INTEGER*4 kunit, kinit, iver
      
      CALL UT990 ( kunit, kinit, iver )

Arguments

    kunit
    File unit for output (see note) [in]
    kinit
    Control parameter (default: 1) [in]
    iver
    Current version [out]

Description

    The subroutine UT990 initializes the different common blocks of the UNILIB libray and displays the value of different constants. To generate a complete initialisation, the argument kinit has to be set to a positive value, e.g. 1. A single common block can be initialized by using a negative value for the argument kinit. To only initialize (or reset) the common block UCddd, the argument kinit has to be set to the value -ddd, i.e. -110, -120, -130, -140, -150, -160, -170, -190 or -192. When the argument is set to an other negative value, none of the common block is initialized.

    The values assigned to the arguments of the different common blocks are displayed in the table below (the floating values are only given with 6 significant digits).

    Common block UC110
    mlblnot init. k1st-1
    klmp-1
    Common block UC120
    nbrfl0 kurfl-1
    mfl(*)not init.
    Common block UC130
    nbrsg0 kursg-1
    mseg(*)not init.
    Common block UC140
    mintinit. by UM510 mextinit. by UM520
    msuninit. by UM520
    Common block UC150
    matminit. by UA610 ntspec-1
    nnspec-1 kspec(*)init. by UA610
    kflag(1...50) 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0
    Common block UC160
    pi 3.14159 deg0.0174533
    re6371.2 gmagmo0.311653
    eclipt23.4415 geoid(1...3)6378.16, 1.00674, 0.00673970
    uma(1...30) 1.6604 10-24, 6.6416 10-24, 33.208 10-24, 66.416 10-24, 139.47410-24, 217.512 10-24, 3.3208 10-24, 46.4912 10-24, 53.1328 10-24, 46.4912 10-24, 73.0576 10-24, 28.2268 10-24, 26.5664 10-24, 0, 0, 28.2268 10-24, 23.2456 10-24, 0, 9.1091 10-28, 1.6604 10-24, 6.6416 10-24, 19.9248 10-24, 23.2456 10-24, 26.5664 10-24, 32.208 10-24, 46.4912 10-24, 49.812 10-24, 53.1328 10-24, 0, 0
    Common block UC170
    nsg0 kgp-1
    mlabnot init. mlinnot init.
    mele(*)not init.
    Common block UC190
    prop0.2 stepx0.075
    stpmin2.0 umsq0.292893
    upsq1.70711 uk2-1.0
    uk30.757359 epskm0.2
    epsrel6.0 10-6 stplst15.0
    xclatnot init. kmflg0
    kum5331
    Common block UC192
    xrmin1000.0 xbmin0.00004
    xtmin0.984808 xbmax1000000.0
    epslon0.05 epsfl0.001
    fvet0.65 pvet3.452
    epsomeg0.008 dltlat1.0

    As output, the subroutine returns the current library version in the argument ifail.

Note

    The argument kunit is used to specify a FORTRAN file unit for the output displayed by the subroutine. When the argument is set to a value less or equal to zero, the printings inside the subroutine are not executed. Commonly when the file unit is set to 6, the output are displayed on the screen (i.e. on the standard output). The FORTRAN statement OPEN has to be used to specify a file name as illustrated below:
      OPEN( unit=3, file='result.txt', status='new' )
      CALL UT990( 3, 1, IVER )
      CLOSE( 3 )

History

    In version 1.10 and earlier, some arguments of common blocks UC190 and UC192 received other values at the initialisation. In common block UC190, the argument nxstp was defined instead of argument kum533. The old values are listed in the table below.

    Old values (version 1.10 and earlier)
    stepx0.4 stpmin5.0
    epskm0.25 stplst25.0
    nxstp400 dltlat2.0

    From version 2.01, the argument kflag(39) of common block UC150 is active and initialized to 1.

    The initialisation of a single common block (kinit = -110, -120, -130, -140, -150, -160, -170, -190 or -192) is only implemented since version 2.03

Diagnostics

    -99001, error on output device

Common Blocks

    UC110, magnetic drift shell description
    UC120, magnetic field line description
    UC130, set of magnetic field line segments
    UC160, general constants
    UC170, temporary magnetic field line
    UC190, control parameters, set 1
    UC192, control parameters, set 2

Dependencies

Called by

    None

Calls

    None

See also

    None

Reported Bugs

    None

Examples

    #1, evaluation of the magnetic field vector
    #2, evaluation of the (B, L) coordinates
    #3, evaluation of the third invariant
    #4, search the point with the lowest altitude on a magnetic drift shell
    #5, namelist to initialize UNILIB
    G.06 How to start using the UNILIB library ?
    T.03 How to customize the magnetic field model used by the library ?
    T.04 Why are Roederer's L* and McIlwain's L parameters different ?
    T.05 How to pass correctly the subroutine arguments in IDL ?