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]

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.