ud330 - evaluate the third invariant¶
- subroutine ud330(phi, star, ifail)¶
- Parameters:
phi [real8,out] :: Third invariant
star [real8,out] :: Roederer’s parameter L*
ifail [integer4,out] :: error flag (see Diagnostics)
Description¶
The subroutine
ud330()
computes the third invariant of the magnetic drift shell stored in the common blockuc110
,uc120
anduc130
. The third invariant is evaluated by the integration of the magnetic flux over a spherical polar cap. Note that the subroutine needs the foot prints of the magnetic field lines contained in the drift shell. To this effect, the argumentfalt
has to be set to a positive value when the subroutineud310()
is called to build the drift shell. The integration steps are controlled by the argumentsepsomeg
anddltlat
of the common blockuc192
.Roederer’s shell parameter L* is deduced directly from the value of the third invariant and of the moment M of the magnetic dipole. The value of M depends on the argument
kmflg
of the common blockuc190
.
References¶
Roederer, J.G., Dynamics of geomagnetically trapped radiation, Springer-Verlag (1970)
Diagnostics¶
None
Common Blocks¶
Dependencies¶
Called by¶
None
Calls¶
See also¶
Reported Bugs¶
None
Examples¶
UNILIB/tags/v3.02