277 |
DO K=1,Nr |
DO K=1,Nr |
278 |
DO j = 1,sNy |
DO j = 1,sNy |
279 |
DO i = 1,sNx+1 |
DO i = 1,sNx+1 |
280 |
tmpMk(i,j,K,bi,bj) = uVel(i,j,K,bi,bj)*0.5 _d 0 |
tmpMk(i,j,K,bi,bj) = uVel(i,j,K,bi,bj)*hFacW(i,j,k,bi,bj) |
281 |
& *(totPhiHyd(i,j,K,bi,bj)+totPhiHyd(i-1,j,K,bi,bj)) |
& *0.5 _d 0*(totPhiHyd(i,j,K,bi,bj)+totPhiHyd(i-1,j,K,bi,bj)) |
282 |
ENDDO |
ENDDO |
283 |
ENDDO |
ENDDO |
284 |
ENDDO |
ENDDO |
293 |
DO K=1,Nr |
DO K=1,Nr |
294 |
DO j = 1,sNy+1 |
DO j = 1,sNy+1 |
295 |
DO i = 1,sNx |
DO i = 1,sNx |
296 |
tmpMk(i,j,K,bi,bj) = vVel(i,j,K,bi,bj)*0.5 _d 0 |
tmpMk(i,j,K,bi,bj) = vVel(i,j,K,bi,bj)*hFacS(i,j,k,bi,bj) |
297 |
& *(totPhiHyd(i,j,K,bi,bj)+totPhiHyd(i,j-1,K,bi,bj)) |
& *0.5 _d 0*(totPhiHyd(i,j,K,bi,bj)+totPhiHyd(i,j-1,K,bi,bj)) |
298 |
ENDDO |
ENDDO |
299 |
ENDDO |
ENDDO |
300 |
ENDDO |
ENDDO |