/[MITgcm]/MITgcm/pkg/mom_common/mom_v_sidedrag.F
ViewVC logotype

Diff of /MITgcm/pkg/mom_common/mom_v_sidedrag.F

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

revision 1.10 by jmc, Fri Sep 30 00:11:34 2005 UTC revision 1.11 by jmc, Wed Oct 12 21:07:01 2005 UTC
# Line 2  C $Header$ Line 2  C $Header$
2  C $Name$  C $Name$
3    
4  #include "MOM_COMMON_OPTIONS.h"  #include "MOM_COMMON_OPTIONS.h"
 #define SIDEDRAG_OLD_VERSION  
5    
6  CBOP  CBOP
7  C !ROUTINE: MOM_V_SIDEDRAG  C !ROUTINE: MOM_V_SIDEDRAG
# Line 54  C !LOCAL VARIABLES: ==================== Line 53  C !LOCAL VARIABLES: ====================
53  C  i,j                  :: loop indices  C  i,j                  :: loop indices
54  C  hFacZClosedE         :: fractional open water to east  C  hFacZClosedE         :: fractional open water to east
55  C  hFacZClosedW         :: fractional open water to west  C  hFacZClosedW         :: fractional open water to west
56        INTEGER I,J        INTEGER i,j
57        _RL hFacZClosedE,hFacZClosedW        _RL hFacZClosedE,hFacZClosedW
58  CEOP  CEOP
   
 #ifdef SIDEDRAG_OLD_VERSION  
59        _RL Ahtmp,A4tmp        _RL Ahtmp,A4tmp
60    
61  C     - Laplacian  and bi-harmonic terms  
62          IF ( sideDragFactor.LE.0. ) THEN
63    C--   Laplacian and bi-harmonic terms: variable-Viscosity coeff. from
64    C     MOM_CALC_VISC are not used here (corresponds to the "old" version).
65    
66        DO j=2-Oly,sNy+Oly-1        DO j=2-Oly,sNy+Oly-1
67         DO i=2-Olx,sNx+Olx-1         DO i=2-Olx,sNx+Olx-1
68  #ifdef NONLIN_FRSURF  #ifdef NONLIN_FRSURF
# Line 91  C       (but this might just be what we Line 92  C       (but this might just be what we
92       &     +hFacZClosedE*_dyU(i+1,j,bi,bj)       &     +hFacZClosedE*_dyU(i+1,j,bi,bj)
93       &      *_recip_dxV(i+1,j,bi,bj) )       &      *_recip_dxV(i+1,j,bi,bj) )
94       &   *drF(k)*2.*(       &   *drF(k)*2.*(
95       &                Ahtmp*vFld(i,j)*cosFacV(J,bi,bj)       &                Ahtmp*vFld(i,j)*cosFacV(j,bi,bj)
96  #ifdef COSINEMETH_III  #ifdef COSINEMETH_III
97       &               -viscA4*del2v(i,j)*sqcosFacV(J,bi,bj)       &               -viscA4*del2v(i,j)*sqcosFacV(j,bi,bj)
98  #else  #else
99       &               -A4tmp*del2v(i,j)*cosFacV(J,bi,bj)       &               -A4tmp*del2v(i,j)*cosFacV(j,bi,bj)
100  #endif  #endif
101       &              )       &              )
102         ENDDO         ENDDO
103        ENDDO        ENDDO
104    
105  #else /* SIDEDRAG_OLD_VERSION */        ELSE
106    C--   Laplacian and bi-harmonic terms: using variable-Viscosity coeff.
107    C     from MOM_CALC_VISC, consistent with dissipation in the interior
108    
 C     - Laplacian  and bi-harmonic terms  
109        DO j=2-Oly,sNy+Oly-1        DO j=2-Oly,sNy+Oly-1
110         DO i=2-Olx,sNx+Olx-1         DO i=2-Olx,sNx+Olx-1
111  C-      this will not give any side-drag along thin wall.  C-      this will not give any side-drag along thin wall.
# Line 121  C       (but this might just be what we Line 123  C       (but this might just be what we
123       &   -_recip_hFacS(i,j,k,bi,bj)       &   -_recip_hFacS(i,j,k,bi,bj)
124       &   *recip_drF(k)*recip_rAs(i,j,bi,bj)       &   *recip_drF(k)*recip_rAs(i,j,bi,bj)
125       &   *( hFacZClosedW*_dyU( i ,j,bi,bj)*_recip_dxV( i ,j,bi,bj)       &   *( hFacZClosedW*_dyU( i ,j,bi,bj)*_recip_dxV( i ,j,bi,bj)
126       &         *(viscAh_Z(i  ,j)*vFld(i,j)*cosFacV(J,bi,bj)       &         *( viscAh_Z(i  ,j)*vFld(i,j)*cosFacV(j,bi,bj)
127  #ifdef COSINEMETH_III  #ifdef COSINEMETH_III
128       &           -viscA4_Z(i  ,j)*del2v(i,j)*sqcosFacV(J,bi,bj))       &           -viscA4_Z(i  ,j)*del2v(i,j)*sqcosFacV(j,bi,bj) )
129  #else  #else
130       &           -viscA4_Z(i  ,j)*del2v(i,j)*cosFacV(J,bi,bj))       &           -viscA4_Z(i  ,j)*del2v(i,j)*cosFacV(j,bi,bj) )
131  #endif  #endif
132       &     +hFacZClosedE*_dyU(i+1,j,bi,bj)*_recip_dxV(i+1,j,bi,bj)       &     +hFacZClosedE*_dyU(i+1,j,bi,bj)*_recip_dxV(i+1,j,bi,bj)
133       &         *(viscAh_Z(i+1,j)*vFld(i,j)*cosFacV(J,bi,bj)       &         *( viscAh_Z(i+1,j)*vFld(i,j)*cosFacV(j,bi,bj)
134  #ifdef COSINEMETH_III  #ifdef COSINEMETH_III
135       &           -viscA4_Z(i+1,j)*del2v(i,j)*sqcosFacV(J,bi,bj))       &           -viscA4_Z(i+1,j)*del2v(i,j)*sqcosFacV(j,bi,bj) )
136  #else  #else
137       &           -viscA4_Z(i+1,j)*del2v(i,j)*cosFacV(J,bi,bj))       &           -viscA4_Z(i+1,j)*del2v(i,j)*cosFacV(j,bi,bj) )
138  #endif  #endif
139       &              )       &    )*drF(k)*sideDragFactor
      &   *drF(k)*2.  
140         ENDDO         ENDDO
141        ENDDO        ENDDO
142    
143  #endif /* (not) SIDEDRAG_OLD_VERSION */  C--   end old-version / new-version IF blocks
144          ENDIF
145    
146  #ifdef ALLOW_DIAGNOSTICS  #ifdef ALLOW_DIAGNOSTICS
147        IF (useDiagnostics) THEN        IF (useDiagnostics) THEN
Legend:
Removed from v.1.10  
changed lines
  Added in v.1.11
  ViewVC Help
Powered by ViewVC 1.1.22