87 |
DO i=iMin,iMax |
DO i=iMin,iMax |
88 |
IF ( hFacW(i,j,kLev,bi,bj) .NE. 0. ) THEN |
IF ( hFacW(i,j,kLev,bi,bj) .NE. 0. ) THEN |
89 |
gU(i,j,kLev,bi,bj) = gU(i,j,kLev,bi,bj) |
gU(i,j,kLev,bi,bj) = gU(i,j,kLev,bi,bj) |
90 |
& -( (1.-maskC(i-1,j,kLev-1,bi,bj)) |
& -( (1.-maskC(i-1,j,kLev-1,bi,bj))*aim_drag(i-1,j,bi,bj) |
91 |
& *aim_drag(i-1,j,bi,bj)*0.0 |
& +(1.-maskC( i ,j,kLev-1,bi,bj))*aim_drag( i ,j,bi,bj) |
92 |
& +(1.-maskC( i ,j,kLev-1,bi,bj)) |
& )* 0.5 _d 0 * uVel(i,j,kLev,bi,bj) |
|
& *aim_drag( i ,j,bi,bj)*0.0 |
|
|
& )* 0.5 _d 0 * uVel(i,j,kLev,bi,bj) |
|
93 |
ENDIF |
ENDIF |
94 |
ENDDO |
ENDDO |
95 |
ENDDO |
ENDDO |