42 |
uVel(i,j,k,bi,bj)=uVel(i,j,k,bi,bj) |
uVel(i,j,k,bi,bj)=uVel(i,j,k,bi,bj) |
43 |
& +deltaTmom*(ab15*gU(i,j,k,bi,bj)+ab05*gUNm1(i,j,k,bi,bj) |
& +deltaTmom*(ab15*gU(i,j,k,bi,bj)+ab05*gUNm1(i,j,k,bi,bj) |
44 |
& -hxFac*pSurfX(i,j)/rhonil |
& -hxFac*pSurfX(i,j)/rhonil |
45 |
& )*maskW(i,j,k,bi,bj) |
& )*_maskW(i,j,k,bi,bj) |
46 |
gUNm1(i,j,k,bi,bj)=gU(i,j,k,bi,bj) |
gUNm1(i,j,k,bi,bj)=gU(i,j,k,bi,bj) |
47 |
ENDDO |
ENDDO |
48 |
ENDDO |
ENDDO |
52 |
vVel(i,j,k,bi,bj)=vVel(i,j,k,bi,bj) |
vVel(i,j,k,bi,bj)=vVel(i,j,k,bi,bj) |
53 |
& +deltaTmom*(ab15*gV(i,j,k,bi,bj)+ab05*gVNm1(i,j,k,bi,bj) |
& +deltaTmom*(ab15*gV(i,j,k,bi,bj)+ab05*gVNm1(i,j,k,bi,bj) |
54 |
& -hyFac*pSurfY(i,j)/rhonil |
& -hyFac*pSurfY(i,j)/rhonil |
55 |
& )*maskS(i,j,k,bi,bj) |
& )*_maskS(i,j,k,bi,bj) |
56 |
gVNm1(i,j,k,bi,bj)=gV(i,j,k,bi,bj) |
gVNm1(i,j,k,bi,bj)=gV(i,j,k,bi,bj) |
57 |
ENDDO |
ENDDO |
58 |
ENDDO |
ENDDO |