/[MITgcm]/MITgcm/pkg/mom_fluxform/mom_v_del2v.F
ViewVC logotype

Diff of /MITgcm/pkg/mom_fluxform/mom_v_del2v.F

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

revision 1.7 by jmc, Tue Dec 5 05:30:38 2006 UTC revision 1.8 by jmc, Mon Apr 25 22:46:41 2011 UTC
# Line 22  C !USES: =============================== Line 22  C !USES: ===============================
22  #include "EEPARAMS.h"  #include "EEPARAMS.h"
23  #include "PARAMS.h"  #include "PARAMS.h"
24  #include "GRID.h"  #include "GRID.h"
25    #include "SURFACE.h"
26    
27  C !INPUT PARAMETERS: ===================================================  C !INPUT PARAMETERS: ===================================================
28  C  bi,bj                :: tile indices  C  bi,bj                :: tile indices
# Line 89  C     del^2 V Line 90  C     del^2 V
90    
91        IF (no_slip_sides) THEN        IF (no_slip_sides) THEN
92  C-- No-slip BCs impose a drag at walls...  C-- No-slip BCs impose a drag at walls...
93        DO j=1-Oly+1,sNy+Oly-1         DO j=1-Oly+1,sNy+Oly-1
94         DO i=1-Olx+1,sNx+Olx-1          DO i=1-Olx+1,sNx+Olx-1
95          hFacZClosedW = _hFacS(i,j,k,bi,bj) - hFacZ(i,j)  #ifdef NONLIN_FRSURF
96          hFacZClosedE = _hFacS(i,j,k,bi,bj) - hFacZ(i+1,j)  C-    this will not give any side-drag along thin wall.
97          del2v(i,j) = del2v(i,j)  C     (but this might just be what we want ...)
98       &  -_recip_hFacS(i,j,k,bi,bj)*recip_drF(k)           hFacZClosedW = h0FacS(i,j,k,bi,bj) - h0FacS(i-1,j,k,bi,bj)
99       &   *recip_rAs(i,j,bi,bj)*recip_deepFac2C(k)           hFacZClosedE = h0FacS(i,j,k,bi,bj) - h0FacS(i+1,j,k,bi,bj)
100       &   *( hFacZClosedW*dyU( i ,j,bi,bj)           hFacZClosedW = MAX( 0. _d 0, hFacZClosedW )
101       &      *_recip_dxV( i ,j,bi,bj)           hFacZClosedE = MAX( 0. _d 0, hFacZClosedE )
102       &     +hFacZClosedE*dyU(i+1,j,bi,bj)  #else
103       &      *_recip_dxV(i+1,j,bi,bj)           hFacZClosedW = _hFacS(i,j,k,bi,bj) - hFacZ(i,j)
104       &    )*drF(k)*2.*vFld(i,j)           hFacZClosedE = _hFacS(i,j,k,bi,bj) - hFacZ(i+1,j)
105       &     *_maskS(i,j,k,bi,bj)  #endif
106             del2v(i,j) = del2v(i,j)
107         &    -_recip_hFacS(i,j,k,bi,bj)
108         &       *recip_rAs(i,j,bi,bj)*recip_deepFac2C(k)
109         &       *( hFacZClosedW*dyU( i ,j,bi,bj)
110         &                      *_recip_dxV( i ,j,bi,bj)
111         &         +hFacZClosedE*dyU(i+1,j,bi,bj)
112         &                      *_recip_dxV(i+1,j,bi,bj)
113         &        )*vFld(i,j)*sideDragFactor
114         &         *_maskS(i,j,k,bi,bj)
115            ENDDO
116         ENDDO         ENDDO
       ENDDO  
117        ENDIF        ENDIF
118    
119        RETURN        RETURN
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22