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

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

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

revision 1.1 by heimbach, Mon Aug 13 18:05:26 2001 UTC revision 1.5 by edhill, Thu Oct 9 04:19:18 2003 UTC
# Line 1  Line 1 
1    C $Header$
2    C $Name$
3    
4    #include "PACKAGES_CONFIG.h"
5  #include "CPP_OPTIONS.h"  #include "CPP_OPTIONS.h"
6    
7  CStartOfInterFace  CBOP
8    C     !ROUTINE: CALC_VISCOSITY
9    C     !INTERFACE:
10        SUBROUTINE CALC_VISCOSITY(        SUBROUTINE CALC_VISCOSITY(
11       I        bi,bj,iMin,iMax,jMin,jMax,K,       I        bi,bj,iMin,iMax,jMin,jMax,K,
      I        maskUp,  
12       O        KappaRU,KappaRV,       O        KappaRU,KappaRV,
13       I        myThid)       I        myThid)
14    
15  C     /==========================================================\  C     !DESCRIPTION: \bv
16  C     | SUBROUTINE CALC_VISCOSITY                                |  C     *==========================================================*
17  C     | o Calculate net viscosity                                |  C     | SUBROUTINE CALC_VISCOSITY                                
18  C     \==========================================================/  C     | o Calculate net vertical viscosity                        
19        IMPLICIT NONE  C     *==========================================================*
20    C     \ev
21    
22    C     !USES:
23          IMPLICIT NONE
24  C     == GLobal variables ==  C     == GLobal variables ==
25  #include "SIZE.h"  #include "SIZE.h"
26  #include "EEPARAMS.h"  #include "EEPARAMS.h"
# Line 21  C     == GLobal variables == Line 28  C     == GLobal variables ==
28  #include "DYNVARS.h"  #include "DYNVARS.h"
29  #include "GRID.h"  #include "GRID.h"
30    
31    C     !INPUT/OUTPUT PARAMETERS:
32  C     == Routine arguments ==  C     == Routine arguments ==
33  C     bi, bj, iMin, iMax, jMin, jMax - Range of points for which calculation  C     bi, bj, iMin, iMax, jMin, jMax - Range of points for which calculation
 C     maskUp  - land/water mask for Wvel points (above tracer level)  
34  C     myThid - Instance number for this innvocation of CALC_COMMON_FACTORS  C     myThid - Instance number for this innvocation of CALC_COMMON_FACTORS
35  C  C     KappaRU :: Total vertical viscosity for zonal flow.
36    C     KappaRV :: Total vertical viscosity for meridional flow.
37        INTEGER bi,bj,iMin,iMax,jMin,jMax,K        INTEGER bi,bj,iMin,iMax,jMin,jMax,K
       _RS maskUp(1-Olx:sNx+Olx,1-Oly:sNy+Oly)  
38        _RL KappaRU(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)        _RL KappaRU(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)
39        _RL KappaRV(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)        _RL KappaRV(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)
40        INTEGER myThid        INTEGER myThid
 CEndOfInterface  
41    
42    C     !LOCAL VARIABLES:
43  C     == Local variables ==  C     == Local variables ==
44  C     I, J, K - Loop counters  C     I, J, K - Loop counters
45        INTEGER i,j        INTEGER i,j
46    CEOP
47    
48        DO j=jMin,jMax        DO j = 1-Oly, sNy+Oly
49         DO i=iMin,iMax         DO i = 1-Olx, sNx+Olx
50          KappaRU(i,j,k) = viscAr          KappaRU(i,j,k) = viscAr
51         ENDDO         ENDDO
52        ENDDO        ENDDO
53    
54        DO j=jMin,jMax        DO j = 1-Oly, sNy+Oly
55         DO i=iMin,iMax         DO i = 1-Olx, sNx+Olx
56          KappaRV(i,j,k) = viscAr          KappaRV(i,j,k) = viscAr
57         ENDDO         ENDDO
58        ENDDO        ENDDO
# Line 52  C     I, J, K - Loop counters Line 60  C     I, J, K - Loop counters
60  #ifdef ALLOW_KPP  #ifdef ALLOW_KPP
61        IF (useKPP) CALL KPP_CALC_VISC(        IF (useKPP) CALL KPP_CALC_VISC(
62       I        bi,bj,iMin+1,iMax,jMin+1,jMax,K,       I        bi,bj,iMin+1,iMax,jMin+1,jMax,K,
      I        maskUp,  
63       O        KappaRU,KappaRV,       O        KappaRU,KappaRV,
64       I        myThid)       I        myThid)
65  #endif  #endif
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22