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