ud321 - evaluate a time integral over a magnetic field line¶
- subroutine ud321(mind, mele, kpar1, rpar2, zsub, result, rnorm, ifail)¶
- Parameters:
mind [zind,in] :: pointers into
mele
mele (*) [zseg,in] :: set of the elementary field line segments
kpar1 [integer4,in] :: control parameters passed to
zsub
rpar2 [real8,in] :: control parameter passed to
zsub
zsub [EXTERNAL,in] :: external subroutine (see note of
ud320()
)result [real8,out] :: result of the integration
rnorm [real8,out] :: time integral of the unity
ifail [integer4,out] :: error flag (see Diagnostics)
Description¶
The subroutine
ud321()
averages physical quantities over a magnetic field line segment in the same way as subroutineud320()
. As input, the subroutine accepts an external subroutine and two parameters to be passed to this subroutine (see note of subroutineud320()
for synopsis and description). As result, the subroutine returns the value of the time integral\(\int f(r) \frac{ds}{\sqrt{1-B/B_m)}}\)
evaluated over the magnetic field line segment for the function provided by the subroutine zsub and for the unity (arguments
result
andrnorm
, respectively). In this equation, ds is a element of length along the field line and the denominator is proportional to the parallel velocity.
History¶
The subroutine
ud321()
does not exist in version 1.05 and earlier.The subroutine
ud321()
has been completely reviewed in version 2.01, as well as subroutineud320()
. The major effects are the addition of a special entry (ud323()
) and a modification of the synopsis of thezsub()
external subroutine (seeud320()
for more details).
Diagnostics¶
-32101, invalid field line description
-32102, wrong elementary segment
-32103, problem with an equatorial single point
Common Blocks¶
None
Dependencies¶
Called by¶
Calls¶
See also¶
None
Reported Bugs¶
None
Examples¶
None
UNILIB/tags/v3.02