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