SUBROUTINE INTRP (FI, ifi, jfi, f, if, jf, xi, yi, x, y, FX, FV, FXY,
                  hx, hy, kl, p, s, c, a, d, work)

ARGUMENTS

OUT

FI(IFI,JFI) interpolated field                                      real
FX(IF,JF)   partial derivative of F with respect to x               real
FY(IF,JF)   partial derivative with respect to y                    real
FXY(IF,JF)  partial second derivative with respect to x and y       real

IN

IFI         x-dimension of array FI                                 integer
JFI         y-dimension of array FI                                 integer 
F(IF,JF)    original field                                          integer
IF          x-dimension of arrays F, FX, FY and FXY                 integer
JF          y-dimension of arrays F, FX, FY and FXY                 integer
XI(IFI)     x location of the values of the interpolated field FI   real
YI(JFI)     y location of the values of the interpolated field FI   real
X(IF)       x location of the values of the original field F        real
Y(JF)       y location of the values of the original field F        real
HX(IF)      grid-lengths defined by HX(i)=X(i+1)-X(i),i=1,JF-1      real
HY(JF)      grid-lengths defined by HY(i)=Y(i+1)-Y(i) i=1,JF-1      real
KL          0 or 1 (how to compute FXY: FXY or FYX)                 integer 
P           working storage vectors of length l=max(IF,JF)          real
S           working storage vectors of length l=max(IF,JF)          real
C           working storage vectors of length l=max(IF,JF)          real
A           working storage vectors of length l=max(IF,JF)          real
D           working storage vectors of length l=max(IF,JF)          real
WORK(JFI,4) work area used by intrp only                            real
  (1,JFI)      (IFI,JFI)  
   +-------------+
   |             | 
   |             |  array FI 
   |             |
   |             |
   |             |
   +-------------+
  (1,1)       (IFI,1)
  (1,JF)      (IF,JF)  
   +-------------+
   |             | 
   |             |  arrays F, FX, FY and FXY 
   |             |
   |             |
   |             |
   +-------------+
  (1,1)       (IF,1)

DESCRIPTION

EXAMPLE

SEE ALSO