173 |
ENDIF |
ENDIF |
174 |
#endif |
#endif |
175 |
|
|
176 |
|
#ifdef ALLOW_ADDFLUID |
177 |
|
IF ( selectAddFluid.GE.1 ) THEN |
178 |
|
DO j=1,sNy |
179 |
|
DO i=1,sNx |
180 |
|
cg2d_b(i,j,bi,bj) = cg2d_b(i,j,bi,bj) |
181 |
|
& - addMass(i,j,k,bi,bj)*mass2rUnit/deltaTmom |
182 |
|
ENDDO |
183 |
|
ENDDO |
184 |
|
#ifdef ALLOW_NONHYDROSTATIC |
185 |
|
IF (use3Dsolver) THEN |
186 |
|
DO j=1,sNy |
187 |
|
DO i=1,sNx |
188 |
|
cg3d_b(i,j,k,bi,bj) = cg3d_b(i,j,k,bi,bj) |
189 |
|
& - addMass(i,j,k,bi,bj)*mass2rUnit/deltaTmom |
190 |
|
ENDDO |
191 |
|
ENDDO |
192 |
|
ENDIF |
193 |
|
#endif |
194 |
|
ENDIF |
195 |
|
#endif /* ALLOW_ADDFLUID */ |
196 |
|
|
197 |
RETURN |
RETURN |
198 |
END |
END |