c c----------------------------------------------------------------------- c Fourth order horizontal pressure gradient (uniform grid). c----------------------------------------------------------------------- c #ifdef cubspl call splintslab(tdepth(1,1,jrsm1),rho(1,1,jrsm1), * d2rhodz2(1,1,jrsm1),imt,km,i-1,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-1,klow(i,3), * khi(i,3),vdepth(i,k,jrn),r(3),xtrp(3)) call splintslab(tdepth(1,1,jrnp1),rho(1,1,jrnp1), * d2rhodz2(1,1,jrnp1),imt,km,i-1,klow(i,4), * khi(i,4),vdepth(i,k,jrn),r(4),xtrp(4)) c call splintslab(tdepth(1,1,jrsm1),rho(1,1,jrsm1), * d2rhodz2(1,1,jrsm1),imt,km,i,klow(i,5), * khi(i,5),vdepth(i,k,jrn),r(5),xtrp(5)) call splintslab(tdepth(1,1,jrs),rho(1,1,jrs), * d2rhodz2(1,1,jrs),imt,km,i,klow(i,6), * khi(i,6),vdepth(i,k,jrn),r(6),xtrp(6)) call splintslab(tdepth(1,1,jrn),rho(1,1,jrn), * d2rhodz2(1,1,jrn),imt,km,i,klow(i,7), * khi(i,7),vdepth(i,k,jrn),r(7),xtrp(7)) call splintslab(tdepth(1,1,jrnp1),rho(1,1,jrnp1), * d2rhodz2(1,1,jrnp1),imt,km,i,klow(i,8), * khi(i,8),vdepth(i,k,jrn),r(8),xtrp(8)) c call splintslab(tdepth(1,1,jrsm1),rho(1,1,jrsm1), * d2rhodz2(1,1,jrsm1),imt,km,i+1,klow(i,9), * khi(i,9),vdepth(i,k,jrn),r(9),xtrp(9)) call splintslab(tdepth(1,1,jrs),rho(1,1,jrs), * d2rhodz2(1,1,jrs),imt,km,i+1,klow(i,10), * khi(i,10),vdepth(i,k,jrn),r(10),xtrp(10)) call splintslab(tdepth(1,1,jrn),rho(1,1,jrn), * d2rhodz2(1,1,jrn),imt,km,i+1,klow(i,11), * khi(i,11),vdepth(i,k,jrn),r(11),xtrp(11)) call splintslab(tdepth(1,1,jrnp1),rho(1,1,jrnp1), * d2rhodz2(1,1,jrnp1),imt,km,i+1,klow(i,12), * khi(i,12),vdepth(i,k,jrn),r(12),xtrp(12)) c call splintslab(tdepth(1,1,jrsm1),rho(1,1,jrsm1), * d2rhodz2(1,1,jrsm1),imt,km,i+2,klow(i,13), * khi(i,13),vdepth(i,k,jrn),r(13),xtrp(13)) call splintslab(tdepth(1,1,jrs),rho(1,1,jrs), * d2rhodz2(1,1,jrs),imt,km,i+2,klow(i,14), * khi(i,14),vdepth(i,k,jrn),r(14),xtrp(14)) call splintslab(tdepth(1,1,jrn),rho(1,1,jrn), * d2rhodz2(1,1,jrn),imt,km,i+2,klow(i,15), * khi(i,15),vdepth(i,k,jrn),r(15),xtrp(15)) call splintslab(tdepth(1,1,jrnp1),rho(1,1,jrnp1), * d2rhodz2(1,1,jrnp1),imt,km,i+2,klow(i,16), * khi(i,16),vdepth(i,k,jrn),r(16),xtrp(16)) #else call lintslab(tdepth(1,1,jrsm1),rho(1,1,jrsm1),imt,km,i-1, * 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-1, * klow(i,3),khi(i,3),vdepth(i,k,jrn),r(3),xtrp(3)) call lintslab(tdepth(1,1,jrnp1),rho(1,1,jrnp1),imt,km,i-1, * klow(i,4),khi(i,4),vdepth(i,k,jrn),r(4),xtrp(4)) c call lintslab(tdepth(1,1,jrsm1),rho(1,1,jrsm1),imt,km,i, * klow(i,5),khi(i,5),vdepth(i,k,jrn),r(5),xtrp(5)) call lintslab(tdepth(1,1,jrs),rho(1,1,jrs),imt,km,i, * klow(i,6),khi(i,6),vdepth(i,k,jrn),r(6),xtrp(6)) call lintslab(tdepth(1,1,jrn),rho(1,1,jrn),imt,km,i, * klow(i,7),khi(i,7),vdepth(i,k,jrn),r(7),xtrp(7)) call lintslab(tdepth(1,1,jrnp1),rho(1,1,jrnp1),imt,km,i, * klow(i,8),khi(i,8),vdepth(i,k,jrn),r(8),xtrp(8)) c call lintslab(tdepth(1,1,jrsm1),rho(1,1,jrsm1),imt,km,i+1, * klow(i,9),khi(i,9),vdepth(i,k,jrn),r(9),xtrp(9)) call lintslab(tdepth(1,1,jrs),rho(1,1,jrs),imt,km,i+1, * klow(i,10),khi(i,10),vdepth(i,k,jrn),r(10), * xtrp(10)) call lintslab(tdepth(1,1,jrn),rho(1,1,jrn),imt,km,i+1, * klow(i,11),khi(i,11),vdepth(i,k,jrn),r(11), * xtrp(11)) call lintslab(tdepth(1,1,jrnp1),rho(1,1,jrnp1),imt,km,i+1, * klow(i,12),khi(i,12),vdepth(i,k,jrn),r(12), * xtrp(12)) c call lintslab(tdepth(1,1,jrsm1),rho(1,1,jrsm1),imt,km,i+2, * klow(i,13),khi(i,13),vdepth(i,k,jrn),r(13), * xtrp(13)) call lintslab(tdepth(1,1,jrs),rho(1,1,jrs),imt,km,i+2, * klow(i,14),khi(i,14),vdepth(i,k,jrn),r(14), * xtrp(14)) call lintslab(tdepth(1,1,jrn),rho(1,1,jrn),imt,km,i+2, * klow(i,15),khi(i,15),vdepth(i,k,jrn),r(15), * xtrp(15)) call lintslab(tdepth(1,1,jrnp1),rho(1,1,jrnp1),imt,km,i+2, * klow(i,16),khi(i,16),vdepth(i,k,jrn),r(16), * xtrp(16)) #endif c kwt(i,k)=0 call extrap(16,vdepth(i,k,jrn),r,xtrp) c tr1=r24*(c27*(r(9)-r(5))+(r(13)-r(1)))*cstr(j-1) tr2=r24*(c27*(r(10)-r(6))+(r(14)-r(2)))*cstr(j) tr3=r24*(c27*(r(11)-r(7))+(r(15)-r(3)))*cstr(j+1) tr4=r24*(c27*(r(12)-r(8))+(r(16)-r(4)))*cstr(j+2) dpx(i,k)=r16*(c9*(tr3+tr2)-tr4-tr1)*dxur(i) c tr1=r24*(c27*(r(3)-r(2))+(r(4)-r(1))) tr2=r24*(c27*(r(7)-r(6))+(r(8)-r(5))) tr3=r24*(c27*(r(11)-r(10))+(r(12)-r(9))) tr4=r24*(c27*(r(15)-r(14))+(r(16)-r(13))) dpy(i,k)=r16*(c9*(tr3+tr2)-tr4-tr1)*dyur(j) c