European Space Agency

Royal Belgian Institute for Aeronomy

ud329 - interpolate between two or three points

subroutine  ud329(mele, ja, jb, jtmp, mgeo, npt, ifail)

Parameters:

• mele (*) [zseg,in] :: field line segment

• ja [integer4,in] :: index of the first point

• jb [integer4,in] :: index of the last point

• jtmp [integer4,in] :: index of another point

• mgeo (npt) [zpnt,in] :: interpolated points (including the points indiced by ja and jb)

• npt [integer4,in] :: number of points

• ifail [integer4,out] :: error flag (see Diagnostics)

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

• -32901, wrong value of the argument npt

• -32902, -32903, not separated points

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

---

UNILIB/tags/v3.03