33 |
DO i=iMin,iMax |
DO i=iMin,iMax |
34 |
gUNm1(i,j,k,bi,bj)=uVel(i,j,k,bi,bj) |
gUNm1(i,j,k,bi,bj)=uVel(i,j,k,bi,bj) |
35 |
& +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) |
36 |
#ifdef ALLOW_CD |
#ifdef INCLUDE_CD_CODE |
37 |
& +guCD(i,j,k,bi,bj) |
& +guCD(i,j,k,bi,bj) |
38 |
#endif |
#endif |
39 |
& )*_maskW(i,j,k,bi,bj) |
& )*_maskW(i,j,k,bi,bj) |
45 |
DO i=iMin,iMax |
DO i=iMin,iMax |
46 |
gVNm1(i,j,k,bi,bj)=vVel(i,j,k,bi,bj) |
gVNm1(i,j,k,bi,bj)=vVel(i,j,k,bi,bj) |
47 |
& +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) |
48 |
#ifdef ALLOW_CD |
#ifdef INCLUDE_CD_CODE |
49 |
& +gvCD(i,j,k,bi,bj) |
& +gvCD(i,j,k,bi,bj) |
50 |
#endif |
#endif |
51 |
& )*_maskS(i,j,k,bi,bj) |
& )*_maskS(i,j,k,bi,bj) |