70 |
A4tmp=min(viscA4+viscA4Grid*(rAs(i,j,bi,bj)**2)/deltaTmom, |
A4tmp=min(viscA4+viscA4Grid*(rAs(i,j,bi,bj)**2)/deltaTmom, |
71 |
& viscA4Max) |
& viscA4Max) |
72 |
IF (viscA4GridMax.GT.0.) THEN |
IF (viscA4GridMax.GT.0.) THEN |
73 |
A4tmp=min(A4tmp,viscA4GridMax*(rAw(i,j,bi,bj)**2)/deltaTmom) |
A4tmp=min(A4tmp,viscA4GridMax*(rAs(i,j,bi,bj)**2)/deltaTmom) |
74 |
ENDIF |
ENDIF |
75 |
A4tmp=max(A4tmp,viscA4GridMin*(rAw(i,j,bi,bj)**2)/deltaTmom) |
A4tmp=max(A4tmp,viscA4GridMin*(rAs(i,j,bi,bj)**2)/deltaTmom) |
76 |
vDragTerms(i,j) = |
vDragTerms(i,j) = |
77 |
& -_recip_hFacS(i,j,k,bi,bj) |
& -_recip_hFacS(i,j,k,bi,bj) |
78 |
& *recip_drF(k)*recip_rAs(i,j,bi,bj) |
& *recip_drF(k)*recip_rAs(i,j,bi,bj) |