subroutine akim2(n) c This routine calls a monotonic spline fit to the observations and then c calculates the profile vcr for all values of r from the spline coefficients. c Note that this fit is performed on the raw observations WITHOUT the Holland c profile subtracted out. c Ray Richardson, 7/18/02 parameter (imx=250, np=26, incfd=1) common /gdinf/ ngd,ngr,ntr,dt,js,jn,ie,iw,iimax,imax,jjmax, * jmax,nstflg,icx,icy,ihx,ihy,dftx,dfty common /sm/ crw(np),rc(np),vc(np),r(imx),vcr(imx),wc(imx) real rc2(n),crw2(n),d(n) integer ierr logical skip skip =.true. do i= 1 , imx vcr(i) = 0.0 enddo do nn=1,n rc2(nn)=rc(nn+1) crw2(nn)=crw(nn+1) end do call pchim(n, rc2, crw2, d, incfd, ierr) call pchfe(n, rc2, crw2, d, incfd, skip, imx, r, vcr, ierr) return end