46 |
CEOP |
CEOP |
47 |
|
|
48 |
C Zonal flux d/dx V |
C Zonal flux d/dx V |
49 |
DO j=1-Oly,sNy+Oly |
DO j=1-Oly+1,sNy+Oly-1 |
50 |
DO i=1-Olx+1,sNx+Olx |
DO i=1-Olx+1,sNx+Olx |
51 |
fZon(i,j) = drF(k)*hFacZ(i,j) |
fZon(i,j) = drF(k)*hFacZ(i,j) |
52 |
& *_dyU(i,j,bi,bj) |
& *_dyU(i,j,bi,bj) |
60 |
|
|
61 |
C Meridional flux d/dy V |
C Meridional flux d/dy V |
62 |
DO j=1-Oly,sNy+Oly-1 |
DO j=1-Oly,sNy+Oly-1 |
63 |
DO i=1-Olx,sNx+Olx |
DO i=1-Olx+1,sNx+Olx-1 |
64 |
fMer(i,j) = drF(k)*_hFacC(i,j,k,bi,bj) |
fMer(i,j) = drF(k)*_hFacC(i,j,k,bi,bj) |
65 |
& *_dxF(i,j,bi,bj) |
& *_dxF(i,j,bi,bj) |
66 |
& *_recip_dyF(i,j,bi,bj) |
& *_recip_dyF(i,j,bi,bj) |
68 |
ENDDO |
ENDDO |
69 |
ENDDO |
ENDDO |
70 |
|
|
71 |
C del^2 U |
C del^2 V |
72 |
DO j=0,sNy+2 |
DO j=1-Oly+1,sNy+Oly-1 |
73 |
DO i=0,sNx+1 |
DO i=1-Olx+1,sNx+Olx-1 |
74 |
del2v(i,j) = recip_drF(k)*_recip_hFacS(i,j,k,bi,bj) |
del2v(i,j) = recip_drF(k)*_recip_hFacS(i,j,k,bi,bj) |
75 |
& *recip_rAs(i,j,bi,bj) |
& *recip_rAs(i,j,bi,bj) |
76 |
& *( fZon(i+1,j) - fZon(i, j ) |
& *( fZon(i+1,j) - fZon(i, j ) |
81 |
|
|
82 |
IF (no_slip_sides) THEN |
IF (no_slip_sides) THEN |
83 |
C-- No-slip BCs impose a drag at walls... |
C-- No-slip BCs impose a drag at walls... |
84 |
DO j=0,sNy+2 |
DO j=1-Oly+1,sNy+Oly-1 |
85 |
DO i=0,sNx+1 |
DO i=1-Olx+1,sNx+Olx-1 |
86 |
hFacZClosedW = _hFacS(i,j,k,bi,bj) - hFacZ(i,j) |
hFacZClosedW = _hFacS(i,j,k,bi,bj) - hFacZ(i,j) |
87 |
hFacZClosedE = _hFacS(i,j,k,bi,bj) - hFacZ(i+1,j) |
hFacZClosedE = _hFacS(i,j,k,bi,bj) - hFacZ(i+1,j) |
88 |
del2v(i,j) = del2v(i,j) |
del2v(i,j) = del2v(i,j) |