SUBROUTINE AKIM(N) c c 6/1/99 - Subroutine modified by Tim Marchok to replace calls c to csakm and csval with a call to scsint since the c IBM scientific library (ESSL) does not contain c either CSAKM or CSVAL. scsint is a cubic spline c interpolation subroutine. c USE setparms c PARAMETER (IMX=95, JMX=1, KMAX=42, LGI=20 , NP = 26) COMMON /GDINF/ NGD,NGR,NTR,DT,JS,JN,IE,IW,IIMAX,IMAX,JJMAX, * JMAX,NSTFLG,ICX,ICY,IHX,IHY,DFTX,DFTY CC COMMON /SM/ CRW(NP),RC(NP),VC(NP),R(IMX),VCR(IMX),WC(IMX) CTM COMMON /POLY/ C(4,NP) CTM DIMENSION BREAK(NP) REAL coeff(N,4) c VCR = 0.0 initflag = 0 if (kind(RC) == real_single .and. kind(CRW) == real_single) then call scsint (RC,CRW,coeff,N,initflag,R,VCR,IMX) else if (kind(RC) == real_double .and. kind(CRW) == real_double) &then call dcsint (RC,CRW,coeff,N,initflag,R,VCR,IMX) endif CTM CALL CSAKM(N,RC,CRW,BREAK,C) CTM CTM DO 10 I = 1 , IMX CTM DO 5 IR = 1 , N-1 CTM IF(R(I).LT.RC(IR).OR.R(I).GT.RC(IR+1))GO TO 5 CTM D = R(I) - RC(IR) CTM VCR(I) = CSVAL(R(I),N-1,BREAK,C) CTM GO TO 10 CTM5 CONTINUE CTM10 CONTINUE RETURN END