59 |
Rj =(tracer(i, j )-tracer(i,j-1))*maskLocS(i, j ) |
Rj =(tracer(i, j )-tracer(i,j-1))*maskLocS(i, j ) |
60 |
Rjm=(tracer(i,j-1)-tracer(i,j-2))*maskLocS(i,j-1) |
Rjm=(tracer(i,j-1)-tracer(i,j-2))*maskLocS(i,j-1) |
61 |
|
|
62 |
c vLoc = vFld(i,j) |
vLoc = vFld(i,j) |
63 |
vLoc = vTrans(i,j)*recip_dxG(i,j,bi,bj) |
c vLoc = vTrans(i,j)*recip_dxG(i,j,bi,bj) |
64 |
& *recip_drF(k)*_recip_hFacS(i,j,k,bi,bj) |
c & *recip_drF(k)*_recip_hFacS(i,j,k,bi,bj) |
65 |
cfl=abs(vLoc*deltaTloc*recip_dyC(i,j,bi,bj)) |
cfl=abs(vLoc*deltaTloc*recip_dyC(i,j,bi,bj)) |
66 |
d0=(2.-cfl)*(1.-cfl)*oneSixth |
d0=(2.-cfl)*(1.-cfl)*oneSixth |
67 |
d1=(1.-cfl*cfl)*oneSixth |
d1=(1.-cfl*cfl)*oneSixth |