c **************************************************************************
c ** subroutine to do B-spline Interpolation **
c **************************************************************************
c ** **
c ** NOTE: link with option '/IMSL'. **
c ** **
c ** FDATA(NXDATA,NYDATA) -- input array to do interpolation **
c ** VALUE(NXVEC,NYVEC) -- output array to do interpolation **
c ** BGX,BGY -- begin value of x and y grid (in agreement with FDATA) **
c ** VGX,VGY -- begin value of x and y grid (in agreement with VALUE) **
c ** G1 -- the grid interval of x & y dimension (in agreement with FDATA)**
c ** G2 -- the grid interval of x & y dimension (in agreement with FDATA)**
c ** **
c **************************************************************************
SUBROUTINE INTERLIB(FDATA,NXDATA,NYDATA,BGX,BGY,
+ VALUE,NXVEC, NYVEC, VGX,VGY, G1,G2)
DIMENSION FDATA(NXDATA,NYDATA), VALUE(NXVEC,NYVEC)
DIMENSION BSCOEF(NXDATA,NYDATA), XDATA(NXDATA), YDATA(NYDATA)
DIMENSION XKNOT(NXDATA+NXVEC), YKNOT(NYDATA+NYVEC)
DIMENSION XVEC(NXVEC), YVEC(NYVEC)
EXTERNAL BS2GD, BS2IN, BSNAK