c----------------------------------------------------------------------
        subroutine DSPLINT(XA,YA,Y2A,N,X,Y)
c is a copy of NUM REC routine splint, in real*8. It gives the value of 
c the splined function at point x. DSPLINE must be called beforehand.
c
c author: Adapted from Numerical Recipes, Kieron Burke, 9.15.99
c     in: n - no of points
c	  xa - vector of x-values
c	  ya - vector of y-values
c	  y2a - vector of coefs
c	  x - chosen x
c    out: y - splined value for x
c  calls: --
      implicit real*8(a-h,o-z)
      DIMENSION XA(N),YA(N),Y2A(N)