58 |
IF ( k.LE.1 ) maskM1 = 0. |
IF ( k.LE.1 ) maskM1 = 0. |
59 |
IF ( k.GE.Nr ) maskP1 = 0. |
IF ( k.GE.Nr ) maskP1 = 0. |
60 |
|
|
61 |
|
C- Initialise output arrays: |
62 |
|
DO j=1-OLy,sNy+OLy |
63 |
|
DO i=1-OLx,sNx+OLx |
64 |
|
uDissip(i,j) = 0. _d 0 |
65 |
|
vDissip(i,j) = 0. _d 0 |
66 |
|
ENDDO |
67 |
|
ENDDO |
68 |
|
|
69 |
C - Laplacian and bi-harmonic terms |
C - Laplacian and bi-harmonic terms |
70 |
c IF (harmonic) THEN |
c IF (harmonic) THEN |
71 |
|
|
98 |
|
|
99 |
DO j=2-OLy,sNy+OLy-1 |
DO j=2-OLy,sNy+OLy-1 |
100 |
DO i=2-OLx,sNx+OLx-1 |
DO i=2-OLx,sNx+OLx-1 |
101 |
vDissip(i,j) = recip_rAs(i,j,bi,bj)*( |
vDissip(i,j) = recip_rAs(i,j,bi,bj)*( |
102 |
& ( viscAh3d_12(i+1, j , k )*str12(i+1, j ,k) |
& ( viscAh3d_12(i+1, j , k )*str12(i+1, j ,k) |
103 |
& *dyU(i+1, j ,bi,bj) |
& *dyU(i+1, j ,bi,bj) |
104 |
& -viscAh3d_12( i , j , k )*str12( i , j ,k) |
& -viscAh3d_12( i , j , k )*str12( i , j ,k) |