| 61 |
& viscAhMax) |
& viscAhMax) |
| 62 |
A4tmp=min(viscA4+viscA4Grid*(rAs(i,j,bi,bj)**2)/deltaTmom, |
A4tmp=min(viscA4+viscA4Grid*(rAs(i,j,bi,bj)**2)/deltaTmom, |
| 63 |
& viscA4Max) |
& viscA4Max) |
| 64 |
|
IF (viscA4GridMax.GT.0.) THEN |
| 65 |
|
A4tmp=min(A4tmp,viscA4GridMax*(rAw(i,j,bi,bj)**2)/deltaTmom) |
| 66 |
|
ENDIF |
| 67 |
|
A4tmp=max(A4tmp,viscA4GridMin*(rAw(i,j,bi,bj)**2)/deltaTmom) |
| 68 |
vDragTerms(i,j) = |
vDragTerms(i,j) = |
| 69 |
& -_recip_hFacS(i,j,k,bi,bj) |
& -_recip_hFacS(i,j,k,bi,bj) |
| 70 |
& *recip_drF(k)*recip_rAs(i,j,bi,bj) |
& *recip_drF(k)*recip_rAs(i,j,bi,bj) |