57 |
DO j=1-Oly,sNy+Oly |
DO j=1-Oly,sNy+Oly |
58 |
DO i=1-Olx,sNx+Olx |
DO i=1-Olx,sNx+Olx |
59 |
Rjp=(tracer(i,j,k,bi,bj)-tracer(i,j,kp1,bi,bj)) |
Rjp=(tracer(i,j,k,bi,bj)-tracer(i,j,kp1,bi,bj)) |
60 |
& *maskC(i,j,kp1,bi,bj) |
& *maskC(i,j,kp1,bi_arg,bj_arg) |
61 |
Rj =(tracer(i,j,km1,bi,bj)-tracer(i,j,k,bi,bj)) |
Rj =(tracer(i,j,km1,bi,bj)-tracer(i,j,k,bi,bj)) |
62 |
& *maskC(i,j,k,bi,bj)*maskC(i,j,km1,bi,bj) |
& *maskC(i,j,k,bi_arg,bj_arg)*maskC(i,j,km1,bi_arg,bj_arg) |
63 |
Rjm=(tracer(i,j,km2,bi,bj)-tracer(i,j,km1,bi,bj)) |
Rjm=(tracer(i,j,km2,bi,bj)-tracer(i,j,km1,bi,bj)) |
64 |
& *maskC(i,j,km1,bi,bj) |
& *maskC(i,j,km1,bi_arg,bj_arg) |
65 |
|
|
66 |
cfl=abs(wVel(i,j,k,bi,bj)*dTarg*recip_drc(k)) |
cfl=abs(wVel(i,j,k,bi_arg,bj_arg)*dTarg*recip_drc(k)) |
67 |
d0=(2.-cfl)*(1.-cfl)*oneSixth |
d0=(2.-cfl)*(1.-cfl)*oneSixth |
68 |
d1=(1.-cfl*cfl)*oneSixth |
d1=(1.-cfl*cfl)*oneSixth |
69 |
c thetaP=0. |
c thetaP=0. |