subroutine akim(n) c @(#)akim.f 1.1 08/07/00 parameter (imx=95, jmx=1, kmax=18, 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) common /poly/ c(4,np) cc dimension break(np) real rsw do 2 i = 1 , imx vcr(i) = 0.0 2 continue cc c cc istat = 0 rsw = 1.0e31 call spline(n, rc, crw, rsw, rsw, break) if (istat .ne. 0) then write(*,'(''warning: istat ='',i5,'' from e01bee'')') end if do 10 i = 1 , imx do 5 ir = 1 , n-1 if(r(i).lt.rc(ir).or.r(i).gt.rc(ir+1))go to 5 d = r(i) - rc(ir) istat = 0 call splint(n,rc,crw,break,r(i),vcr(i)) if (istat .ne. 0) then write(*,'(''warning: istat ='',i5,'' from e01bfe'')') end if go to 10 5 continue 10 continue cc return end