33 |
INTEGER myThid |
INTEGER myThid |
34 |
|
|
35 |
C == Local variables == |
C == Local variables == |
36 |
|
C wFld :: velocity, vertical component |
37 |
INTEGER i,j,kp1,km1,km2,bi,bj |
INTEGER i,j,kp1,km1,km2,bi,bj |
38 |
_RL Rjm,Rj,Rjp,cfl,d0,d1 |
_RL Rjm,Rj,Rjp,cfl,d0,d1 |
39 |
_RL psiP,psiM,thetaP,thetaM |
_RL psiP,psiM,thetaP,thetaM |
40 |
|
_RL wFld |
41 |
|
|
42 |
IF (.NOT. multiDimAdvection) THEN |
IF (.NOT. multiDimAdvection) THEN |
43 |
C If using the standard time-stepping/advection schemes (ie. AB-II) |
C If using the standard time-stepping/advection schemes (ie. AB-II) |
65 |
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)) |
66 |
& *maskC(i,j,km1,bi_arg,bj_arg) |
& *maskC(i,j,km1,bi_arg,bj_arg) |
67 |
|
|
68 |
cfl=abs(wVel(i,j,k,bi_arg,bj_arg)*dTarg*recip_drc(k)) |
c wFld = wVel(i,j,k,bi_arg,bj_arg) |
69 |
|
wFld = rTrans(i,j)*recip_rA(i,j,bi_arg,bj_arg) |
70 |
|
cfl=abs(wFld*dTarg*recip_drC(k)) |
71 |
d0=(2.-cfl)*(1.-cfl)*oneSixth |
d0=(2.-cfl)*(1.-cfl)*oneSixth |
72 |
d1=(1.-cfl*cfl)*oneSixth |
d1=(1.-cfl*cfl)*oneSixth |
73 |
c thetaP=0. |
c thetaP=0. |