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

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

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

revision 1.13 by baylor, Mon Sep 26 15:27:11 2005 UTC revision 1.14 by jmc, Tue Sep 27 00:18:20 2005 UTC
# Line 1  Line 1 
1    C $Header$
2    C $Name$
3    
4  #include "MOM_COMMON_OPTIONS.h"  #include "MOM_COMMON_OPTIONS.h"
5    
# Line 70  C     == Global variables == Line 72  C     == Global variables ==
72  #include "GRID.h"  #include "GRID.h"
73  #include "EEPARAMS.h"  #include "EEPARAMS.h"
74  #include "PARAMS.h"  #include "PARAMS.h"
 #ifdef ALLOW_EXCH2  
 #include "W2_EXCH2_TOPOLOGY.h"  
 #include "W2_EXCH2_PARAMS.h"  
 #endif /* ALLOW_EXCH2 */  
75    
76  C     == Routine arguments ==  C     == Routine arguments ==
77        INTEGER bi,bj,k        INTEGER bi,bj,k
# Line 118  C     == Local variables == Line 116  C     == Local variables ==
116        _RL viscA4_ZSmg(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL viscA4_ZSmg(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
117        _RL viscA4_DSmg(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL viscA4_DSmg(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
118        LOGICAL calcLeith,calcSmag        LOGICAL calcLeith,calcSmag
       LOGICAL northWestCorner, northEastCorner,  
      &        southWestCorner, southEastCorner  
 #ifdef ALLOW_EXCH2  
       INTEGER myTile  
 #endif /* ALLOW_EXCH2 */  
   
 C     Special stuff for Cubed Sphere  
       southWestCorner = .FALSE.  
       southEastCorner = .FALSE.  
       northWestCorner = .FALSE.  
       northEastCorner = .FALSE.  
       IF (useCubedSphereExchange) THEN  
 #ifdef ALLOW_EXCH2  
        myTile = W2_myTileList(bi)  
        IF ( exch2_isWedge(myTile) .EQ. 1 .AND.  
      &      exch2_isSedge(myTile) .EQ. 1 ) THEN  
         southWestCorner = .TRUE.  
        ENDIF  
        IF ( exch2_isEedge(myTile) .EQ. 1 .AND.  
      &      exch2_isSedge(myTile) .EQ. 1 ) THEN  
         southEastCorner = .TRUE.  
        ENDIF  
        IF ( exch2_isEedge(myTile) .EQ. 1 .AND.  
      &      exch2_isNedge(myTile) .EQ. 1 ) THEN  
         northEastCorner = .TRUE.  
        ENDIF  
        IF ( exch2_isWedge(myTile) .EQ. 1 .AND.  
      &      exch2_isNedge(myTile) .EQ. 1 ) THEN  
         northWestCorner = .TRUE.  
        ENDIF  
 #else  
        southWestCorner = .TRUE.  
        southEastCorner = .TRUE.  
        northWestCorner = .TRUE.  
        northEastCorner = .TRUE.  
 #endif /* ALLOW_EXCH2 */  
       ENDIF  
119    
120        useVariableViscosity=        useVariableViscosity=
121       &      (viscAhGrid.NE.0.)       &      (viscAhGrid.NE.0.)
# Line 354  C These are (powers of) length scales Line 315  C These are (powers of) length scales
315    
316           L2rdt=0.25 _d 0*recip_dt*L2           L2rdt=0.25 _d 0*recip_dt*L2
317           IF (useAreaViscLength) THEN           IF (useAreaViscLength) THEN
318            L4rdt=0.125 _d 0*recip_dt*RaZ(i,j,bi,bj)**2            L4rdt=0.125 _d 0*recip_dt*rAz(i,j,bi,bj)**2
319           ELSE           ELSE
320            L4rdt=recip_dt/            L4rdt=recip_dt/
321       &     ( 6. _d 0*(recip_DXV(I,J,bi,bj)**4+recip_DYU(I,J,bi,bj)**4)       &     ( 6. _d 0*(recip_DXV(I,J,bi,bj)**4+recip_DYU(I,J,bi,bj)**4)
Legend:
Removed from v.1.13  
changed lines
  Added in v.1.14
  ViewVC Help
Powered by ViewVC 1.1.22