trend

- 4

Trapped Radiation ENvironment model Development

UD321

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

Name

    Subroutine UD321 - evaluate a time integral over a magnetic field line

Synopsis

    
      REAL*8    rpar2, result, rnorm
      INTEGER*4 kpar1, ifail
      EXTERNAL  zsub
      RECORD /zseg/ mele(*)
      RECORD /zind/ mind

      CALL UD321 (mind, mele, kpar1, rpar2, zsub, result, rnorm, ifail)   

Arguments

    mind
    pointers into mele [in]
    mele
    Set of the elementary field line segments [in]
    kpar1
    control parameter passed to zsub [in]
    rpar2
    control parameter passed to zsub [in]
    zsub
    external subroutine (see note of subroutine UD320) [in]
    result
    result of the integration [out]
    rnorm
    time integral of the unity [out]
    ifail
    Error flag (see diagnostics) [out]

Description

    The subroutine UD321 averages physical quantities over a magnetic field line segment in the same way as subroutine UD320. As input, the subroutine accepts an external subroutine and two parameters to be passed to this subroutine (see note of subroutine UD320 for synopsis and description). As result, the subroutine returns the value of the time integral
    [Eq]
    evaluated over the magnetic field line segment for the function provided by the subroutine zsub and for the unity (arguments result and rnorm, respectively). In this equation, ds is a element of length along the field line and the denominators 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 subroutine UD320. The major effects are the addition of a special entry (UD323) and a modification of the synopsis of the zsub external subroutine (see UD320 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

    UD320, evaluate a double time integral over a drift shell

Calls

    UD328, evaluate the second derivative of B
    UD329, interpolate between two or three points

See also

    None

Reported Bugs

    None

Examples

    None