trend

- 4

Trapped Radiation ENvironment model Development

UD329

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

Name

    Subroutine UD329 - interpolate between two or three points

Synopsis

 
      RECORD /zseg/ mele(*)
      INTEGER*4     ja, jb, jtmp, ifail, npt      
      RECORD /zpnt/ mgeo(npt)

      CALL UD329 (mele, ja, jb, jtmp, mgeo, npt,ifail)
   

Arguments

    mele
    field line segment [in]
    ja
    index of the first point [in]
    jb
    index of the last point [in]
    jtmp
    index of another point [in]
    mgeo
    interpolated points (including the points indiced by ja and jb) [out]
    npt
    number of points [in]
    ifail
    Error flag (see diagnostics) [out]

Description

The subroutine UD329 interpolates npt-2 points along a magnetic field line segment. The indices of the points are passed by the arguments ja, jb and jtmp. If the argument jtmp is set to a value less than zero, the subroutine makes a linear interpolation between the two points mele(ja) and mele(jb). Otherwise, a parabolic interpolation is applied. The magnetic field vector is evaluated at the interpolated points.

History

    The subroutine UD329 does not exist in version 1.05 and earlier.

Diagnostics

Common Blocks

    None

Dependencies

Called by

    UD321, evaluate a time integral over a magnetic field line

Calls

    UM530, evaluate the magnetic field vector
    UT541, convert spherical coordinates to cartesian coordinates
    UT546, convert cartesian coordinates to spherical coordinates
    UT999, radius of curvature in a dipolar magnetic field

See also

    None

Reported Bugs

    None

Examples

    None