47 |
C i,j :: loop indices |
C i,j :: loop indices |
48 |
INTEGER I,J |
INTEGER I,J |
49 |
_RS gridScalingA4 |
_RS gridScalingA4 |
50 |
|
_RS gridScalingAh |
51 |
CEOP |
CEOP |
52 |
|
|
53 |
C - Laplacian and bi-harmonic terms |
C - Laplacian and bi-harmonic terms |
54 |
DO j=1-Oly+1,sNy+Oly |
DO j=1-Oly+1,sNy+Oly |
55 |
DO i=1-Olx,sNx+Olx |
DO i=1-Olx,sNx+Olx |
56 |
gridScalingA4 = (_dyU(i,j,bi,bj)**4)/(5000.**4) |
gridScalingA4 = (_dyU(i,j,bi,bj)**4)/(5000.**4) |
57 |
|
gridScalingAh = (_dyU(i,j,bi,bj)**2)/(5000.**2) |
58 |
yViscFluxU(i,j) = |
yViscFluxU(i,j) = |
59 |
& _dxV(i,j,bi,bj)*drF(k)*hFacZ(i,j) |
& _dxV(i,j,bi,bj)*drF(k)*hFacZ(i,j) |
60 |
& *( |
& *( |
61 |
& -viscAh*(uFld(i,j)-uFld(i,j-1)) |
& -viscAh*(uFld(i,j)-uFld(i,j-1)) |
62 |
|
& *gridScalingAh |
63 |
#ifdef ISOTROPIC_COS_SCALING |
#ifdef ISOTROPIC_COS_SCALING |
64 |
& *cosFacV(J,bi,bj) |
& *cosFacV(J,bi,bj) |
65 |
#endif |
#endif |