/[MITgcm]/MITgcm/model/src/calc_div_ghat.F
ViewVC logotype

Diff of /MITgcm/model/src/calc_div_ghat.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.24 by jmc, Sun Aug 2 19:54:50 2009 UTC revision 1.25 by jmc, Sun Sep 27 23:18:53 2009 UTC
# Line 95  C      => Filtering of uVel,vVel is nece Line 95  C      => Filtering of uVel,vVel is nece
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)
Legend:
Removed from v.1.24  
changed lines
  Added in v.1.25
  ViewVC Help
Powered by ViewVC 1.1.22