99 |
faceArea = _dyG(I,J,bi,bj)*drF(K) |
faceArea = _dyG(I,J,bi,bj)*drF(K) |
100 |
& *_hFacW(I,J,K,bi,bj) |
& *_hFacW(I,J,K,bi,bj) |
101 |
aW2d(I,J,bi,bj) = aW2d(I,J,bi,bj) + |
aW2d(I,J,bi,bj) = aW2d(I,J,bi,bj) + |
102 |
& implicSurfPress*implicDiv2DFlow* |
& implicSurfPress*implicDiv2DFlow |
103 |
& gBaro*faceArea*recip_dxC(I,J,bi,bj) |
& *faceArea*recip_dxC(I,J,bi,bj) |
104 |
faceArea = _dxG(I,J,bi,bj)*drF(K) |
faceArea = _dxG(I,J,bi,bj)*drF(K) |
105 |
& *_hFacS(I,J,K,bi,bj) |
& *_hFacS(I,J,K,bi,bj) |
106 |
aS2d(I,J,bi,bj) = aS2d(I,J,bi,bj) + |
aS2d(I,J,bi,bj) = aS2d(I,J,bi,bj) + |
107 |
& implicSurfPress*implicDiv2DFlow* |
& implicSurfPress*implicDiv2DFlow |
108 |
& gBaro*faceArea*recip_dyC(I,J,bi,bj) |
& *faceArea*recip_dyC(I,J,bi,bj) |
109 |
ENDDO |
ENDDO |
110 |
ENDDO |
ENDDO |
111 |
ENDDO |
ENDDO |
194 |
aC = -( |
aC = -( |
195 |
& aW2d(I,J,bi,bj) + aW2d(I+1,J ,bi,bj) |
& aW2d(I,J,bi,bj) + aW2d(I+1,J ,bi,bj) |
196 |
& +aS2d(I,J,bi,bj) + aS2d(I ,J+1,bi,bj) |
& +aS2d(I,J,bi,bj) + aS2d(I ,J+1,bi,bj) |
197 |
& +freeSurfFac*myNorm* horiVertRatio* |
& +freeSurfFac*myNorm*recip_Bo(I,J,bi,bj)* |
198 |
& rA(I,J,bi,bj)/deltaTMom/deltaTMom |
& rA(I,J,bi,bj)/deltaTMom/deltaTMom |
199 |
& ) |
& ) |
200 |
aCs = -( |
aCs = -( |
201 |
& aW2d(I,J-1,bi,bj) + aW2d(I+1,J-1,bi,bj) |
& aW2d(I,J-1,bi,bj) + aW2d(I+1,J-1,bi,bj) |
202 |
& +aS2d(I,J-1,bi,bj) + aS2d(I ,J ,bi,bj) |
& +aS2d(I,J-1,bi,bj) + aS2d(I ,J ,bi,bj) |
203 |
& +freeSurfFac*myNorm* horiVertRatio* |
& +freeSurfFac*myNorm*recip_Bo(I,J-1,bi,bj)* |
204 |
& rA(I,J-1,bi,bj)/deltaTMom/deltaTMom |
& rA(I,J-1,bi,bj)/deltaTMom/deltaTMom |
205 |
& ) |
& ) |
206 |
aCw = -( |
aCw = -( |
207 |
& aW2d(I-1,J,bi,bj) + aW2d(I ,J ,bi,bj) |
& aW2d(I-1,J,bi,bj) + aW2d(I ,J ,bi,bj) |
208 |
& +aS2d(I-1,J,bi,bj) + aS2d(I-1,J+1,bi,bj) |
& +aS2d(I-1,J,bi,bj) + aS2d(I-1,J+1,bi,bj) |
209 |
& +freeSurfFac*myNorm* horiVertRatio* |
& +freeSurfFac*myNorm*recip_Bo(I-1,J,bi,bj)* |
210 |
& rA(I-1,J,bi,bj)/deltaTMom/deltaTMom |
& rA(I-1,J,bi,bj)/deltaTMom/deltaTMom |
211 |
& ) |
& ) |
212 |
IF ( aC .EQ. 0. ) THEN |
IF ( aC .EQ. 0. ) THEN |