c
c-----------------------------------------------------------------------
c  Second order horizontal pressure gradient.
c-----------------------------------------------------------------------
c
#ifdef cubspl
          call splintslab(tdepth(1,1,jrs),rho(1,1,jrs),
     *                    d2rhodz2(1,1,jrs),imt,km,i,klow(i,1),
     *                    khi(i,1),vdepth(i,k,jrn),r(1),xtrp(1))
          call splintslab(tdepth(1,1,jrs),rho(1,1,jrs),
     *                    d2rhodz2(1,1,jrs),imt,km,i+1,klow(i,2),
     *                    khi(i,2),vdepth(i,k,jrn),r(2),xtrp(2))
          call splintslab(tdepth(1,1,jrn),rho(1,1,jrn),
     *                    d2rhodz2(1,1,jrn),imt,km,i,klow(i,3),
     *                    khi(i,3),vdepth(i,k,jrn),r(3),xtrp(3))
          call splintslab(tdepth(1,1,jrn),rho(1,1,jrn),
     *                    d2rhodz2(1,1,jrn),imt,km,i+1,klow(i,4),
     *                    khi(i,4),vdepth(i,k,jrn),r(4),xtrp(4))
#else
          call lintslab(tdepth(1,1,jrs),rho(1,1,jrs),imt,km,i,
     *                  klow(i,1),khi(i,1),vdepth(i,k,jrn),r(1),xtrp(1))
          call lintslab(tdepth(1,1,jrs),rho(1,1,jrs),imt,km,i+1,
     *                  klow(i,2),khi(i,2),vdepth(i,k,jrn),r(2),xtrp(2))
          call lintslab(tdepth(1,1,jrn),rho(1,1,jrn),imt,km,i,
     *                  klow(i,3),khi(i,3),vdepth(i,k,jrn),r(3),xtrp(3))
          call lintslab(tdepth(1,1,jrn),rho(1,1,jrn),imt,km,i+1,
     *                  klow(i,4),khi(i,4),vdepth(i,k,jrn),r(4),xtrp(4))
#endif
c
#ifdef noxtrp
          if(xtrp(1).or.xtrp(2).or.xtrp(3).or.xtrp(4)) then
            kwt(i,k)=1
          else
            kwt(i,k)=0
          endif
# else
          kwt(i,k)=0
          call extrap(4,vdepth(i,k,jrn),r,xtrp)
#endif
c
          dpx(i,k)=((r(4)-r(3))*cstr(j+1)+(r(2)-r(1))*cstr(j))*dxu2r(i)
          dpy(i,k)=((r(4)-r(2))+(r(3)-r(1)))*dyu2r(j)
c