50 |
Rj =(tracer(i, j )-tracer(i,j-1))*maskLocS(i, j ) |
Rj =(tracer(i, j )-tracer(i,j-1))*maskLocS(i, j ) |
51 |
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) |
52 |
|
|
53 |
c vLoc = vFld(i,j) |
vLoc = vFld(i,j) |
54 |
vLoc = vTrans(i,j)*recip_dxG(i,j,bi,bj) |
c vLoc = vTrans(i,j)*recip_dxG(i,j,bi,bj) |
55 |
& *recip_drF(k)*_recip_hFacS(i,j,k,bi,bj) |
c & *recip_drF(k)*_recip_hFacS(i,j,k,bi,bj) |
56 |
cfl=abs(vLoc*deltaTloc*recip_dyC(i,j,bi,bj)) |
cfl=abs(vLoc*deltaTloc*recip_dyC(i,j,bi,bj)) |
57 |
d0=(2. _d 0 -cfl)*(1. _d 0 -cfl)*oneSixth |
d0=(2. _d 0 -cfl)*(1. _d 0 -cfl)*oneSixth |
58 |
d1=(1. _d 0 -cfl*cfl)*oneSixth |
d1=(1. _d 0 -cfl*cfl)*oneSixth |