Previous: fittpr Up: ../plot79_f.html Next: fittsr

FITTSF

       SUBROUTINE  FITTSF (ZE, MX,MY, NX,NY, X,Y, P,NP, WORK,MAXWRK,
      X                    IWORK,MAXIWK)
 C$    (Triangulate - Surface Grid)
 C$    This routine provides a convenient user interface to  Frank
 C$    Little's triangulation  routines, FITT**.   It is  intended
 C$    for general  use  when  an arbitrary  set  of  data  points
 C$    defining a surface is available, and from which is  desired
 C$    the  construction   of  an   interpolated  surface   on   a
 C$    rectangular grid.  If the  grid points are  equally-spaced,
 C$    the surface can  then be  input to  any of  the contour  or
 C$    hidden-line routines for graphing.  The arguments are:
 C$
 C$    ZE(NX,NY).......Interpolated surface values (returned).
 C$    (MX,MY).........Actual declared dimensions of ZE(*,*).
 C$    (X(NX),Y(NY))...Coordinate  arrays  defining  the  grid  at
 C$                    which  interpolation   is   desired.    For
 C$                    graphing  purposes,   these  must   contain
 C$                    equally-spaced   points,    and   can    be
 C$                    conveniently   initialized   with   utility
 C$                    routine UTRIV.
 C$    P(3,NP).........Array containing (X,Y,Z) coordinates of
 C$                    the points defining the surface.
 C$    NP..............Number of points defining the surface.
 C$    WORK(*).........REAL scratch array.
 C$    MAXWRK..........Dimension of WORK(*).  This should be at
 C$                    least 3*NP.
 C$    IWORK(*)........INTEGER  scratch  array.   NOTE:   It    is
 C$                    assumed that this  space can  be passed  to
 C$                    other routines and used as a LOGICAL array.
 C$                    This is  a violation  of Portable  FORTRAN,
 C$                    but does work on most machines.
 C$    MAXIWK..........Dimension of IWORK(*).  This should be at
 C$                    least (10.5*NP + 52*SQRT(NP)).
 C$    (19-MAR-82)