- 4

Trapped Radiation ENvironment model Development

Example 3 - evaluation of the third invariant

[ Description | Source | Results ]


    Evaluation of the third invariant of a magnetic drift shell labeled by a pair of (Bm, Lm) parameters. The magnetic field is described with the help of the DGRF/IGRF 1995 geomagnetic model.

    The program contains three parts: the initialisation of the library and the magnetic model (UT990, UM510 and UM520), the tracing of the magnetic drift shell (UD310) and the evaluation of the third invariant (UD330). Note that subroutine UD317 can be used instead of subroutine UD310 (see note).

    In this example the drift shell parameters are set to Bm = 0.19 Gauss and Lm = 2. The foot points of the magnetic drift shell in the Northern hemisphere are shown on the lefthand panel of the figure below. The crosses indicate the points used to evaluate the magnetic flux. In this case, the third invariant is equal to 3.8 10+7 Gauss km2 which corresponds to L* = 2.07. The righthand panel of the figure shows the foot points in the Southern hemisphere.

    Note that the foot points of the drift shell are stored in the common block UC120 [arguments mfl().footpn and mfl().footps ].


      PROGRAM example3
      INCLUDE 'structure.h'
      INTEGER*4    kunit, kinit, ifail, kint, kext
      CHARACTER*32 lbint, lbext
      REAL*8       year, param(10), amjd
      INTEGER*4    knfl, ktyplus
      REAL*8       fbm0, flm0, falt, phi, star
C     initialize variables
      DATA kunit, kinit/ 6, 1/
      DATA kint, kext/ 0, 0/
      DATA year, param, amjd/ 1995.d0, 10*0.0d0, 0.0d0/
C     initialize the library
      CALL UT990 (kunit, kinit, ifail)
      IF( ifail .LT. 0 )STOP
C     set the geomagnetic field model
C     (kint=0, DGRF/IGRF)
      CALL UM510 (kint, year, lbint, kunit, ifail)
      IF( ifail .LT. 0 )STOP
C     set no external magnetic field model
      CALL UM520 (kext, amjd, param,
     :           lbext, kunit, ifail)
      IF( ifail .LT. 0 )STOP
C     trace the drift shell
      fbm0          =    0.19d0
      flm0          =    2.0d0
      falt          =    0.0d0
      knfl          =   60
      ktyplus       =    3
      CALL UD310 (fbm0, flm0, falt, knfl, 
     :           ktyplus, ifail)
      IF( ifail .LT. 0 )STOP
C     evaluate the third invariant
      CALL UD330 (phi, star, ifail) 
      IF( ifail .LT. 0 )STOP
C     write the result
      WRITE(*,*) phi,star


                                UNILIB Library
                                 Version  2.01
                        Developed by BIRA/IASB for ESA


--- Main control parameters ---

     UC160 (general constants):
                     Geoid major axis = 6378.16     km
                     Geoid minor axis = 6356.77     km
               Mean Earth radius (Re) = 6371.20     km
                 Ecliptic inclination =   23.441512 deg
                McIlwain Earth dipole =    0.311653 G/Re^3

     UC190 (field line tracing):
                    Maximum step size =  477.8      km at one Re
             Maximum "last" step size =   15.0      km
                    Minimum step size =    2.0      km
                   Altitude precision =  200.       m
             Magnetic field precision =    0.0006   %

     UC192 (magnetic field evaluation):
            Max. geomagnetic latitude =   85.0      deg
                       Minimum radius =    0.157    Re
           (drift shell tracing):
                       Precision on L =    0.10000  %
           Precision on the longitude =    0.05     deg

--- Geomagnetic field model ---

     Model ( 0): IGRF 1995                          Epoch 1995.
                            Order + 1 =   11
                    Calculation epoch = 1995.0      year
        Colatitude of the dipole pole =   10.70     deg
         Longitude of the dipole pole =  -71.41     deg
                  Earth dipole moment =    0.302077 G/Re^3
         Correction for the SAA drift =    0.00     deg

--- External magnetic field model ---

     Model (0): None
          Planetary activity index Kp =    0.0      (magnetopause check)
                       Universal Time =    0.00     hrs
           Greenwich angle of equinox =  100.08     deg
          GEO direction of the Sun, X =   -0.9200
                                    Y =   -0.0131
                                    Z =   -0.3918
                    dipole tilt angle =  -25.92     deg
 -38317592.4315653        2.07441404002132