95 |
C-- Now the filter are applied in the_correction_step(). |
C-- Now the filter are applied in the_correction_step(). |
96 |
C We have left this code here to indicate where the filters used to be |
C We have left this code here to indicate where the filters used to be |
97 |
C in the algorithm before JMC moved them to after the pressure solver. |
C in the algorithm before JMC moved them to after the pressure solver. |
98 |
C- |
c#ifdef ALLOW_ZONAL_FILT |
99 |
C#ifdef ALLOW_ZONAL_FILT |
c IF (zonal_filt_lat.LT.90.) THEN |
100 |
C IF (zonal_filt_lat.LT.90.) THEN |
c CALL ZONAL_FILTER( |
101 |
C CALL ZONAL_FILTER( |
c U uVel( 1-OLx,1-OLy,k,bi,bj), |
102 |
C & uVel, hFacW, 1-1, sNy+1, k, k, bi, bj, 1, myThid) |
c I hFacW(1-OLx,1-OLy,k,bi,bj), |
103 |
C CALL ZONAL_FILTER( |
c I 0, sNy+1, 1, bi, bj, 1, myThid ) |
104 |
C & vVel, hFacS, 1-1, sNy+1, k, k, bi, bj, 2, myThid) |
c CALL ZONAL_FILTER( |
105 |
C ENDIF |
c U vVel( 1-OLx,1-OLy,k,bi,bj), |
106 |
C#endif |
c I hFacS(1-OLx,1-OLy,k,bi,bj), |
107 |
|
c I 0, sNy+1, 1, bi, bj, 2, myThid ) |
108 |
|
c ENDIF |
109 |
|
c#endif |
110 |
DO j=1,sNy |
DO j=1,sNy |
111 |
DO i=1,sNx+1 |
DO i=1,sNx+1 |
112 |
pf(i,j) = ( implicDiv2Dflow * gU(i,j,k,bi,bj) |
pf(i,j) = ( implicDiv2Dflow * gU(i,j,k,bi,bj) |