/[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.7 by mlosch, Thu Sep 16 09:35:11 2004 UTC revision 1.8 by jmc, Mon Oct 20 23:51:39 2008 UTC
# Line 7  C $Name$ Line 7  C $Name$
7  CBOP  CBOP
8  C     !ROUTINE: CALC_VISCOSITY  C     !ROUTINE: CALC_VISCOSITY
9  C     !INTERFACE:  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,
12       O        KappaRU,KappaRV,       O           KappaRU, KappaRV,
13       I        myThid)       I           myThid )
14    
15  C     !DESCRIPTION: \bv  C     !DESCRIPTION: \bv
16  C     *==========================================================*  C     *==========================================================*
17  C     | SUBROUTINE CALC_VISCOSITY                                  C     | SUBROUTINE CALC_VISCOSITY
18  C     | o Calculate net vertical viscosity                          C     | o Calculate net vertical viscosity
19  C     *==========================================================*  C     *==========================================================*
20  C     \ev  C     \ev
21    
# Line 30  C     == GLobal variables == Line 30  C     == GLobal variables ==
30    
31  C     !INPUT/OUTPUT PARAMETERS:  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     iMin,iMax,jMin,jMax :: Range of points for which calculation
34  C     myThid - Instance number for this innvocation of CALC_COMMON_FACTORS  C     bi,bj   :: current tile indices
35  C     KappaRU :: Total vertical viscosity for zonal flow.  C     KappaRU :: Total vertical viscosity for zonal flow.
36  C     KappaRV :: Total vertical viscosity for meridional flow.  C     KappaRV :: Total vertical viscosity for meridional flow.
37        INTEGER bi,bj,iMin,iMax,jMin,jMax,K  C     myThid  :: my Thread Id number
38          INTEGER iMin,iMax,jMin,jMax
39          INTEGER bi,bj
40        _RL KappaRU(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)        _RL KappaRU(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)
41        _RL KappaRV(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)        _RL KappaRV(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)
42        INTEGER myThid        INTEGER myThid
43    
44  C     !LOCAL VARIABLES:  C     !LOCAL VARIABLES:
45  C     == Local variables ==  C     == Local variables ==
46  C     I, J, K - Loop counters  C     i, j, k :: Loop counters
47        INTEGER i,j        INTEGER i,j,k
48  CEOP  CEOP
49    
50        DO j = 1-Oly, sNy+Oly        DO k = 1,Nr
        DO i = 1-Olx, sNx+Olx  
         KappaRU(i,j,k) = viscAr  
        ENDDO  
       ENDDO  
51    
52        DO j = 1-Oly, sNy+Oly         DO j = 1-Oly, sNy+Oly
53         DO i = 1-Olx, sNx+Olx          DO i = 1-Olx, sNx+Olx
54          KappaRV(i,j,k) = viscAr           KappaRU(i,j,k) = viscAr
55             KappaRV(i,j,k) = viscAr
56            ENDDO
57         ENDDO         ENDDO
       ENDDO  
58    
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,iMax,jMin,jMax, k,
62       O        KappaRU,KappaRV,       O        KappaRU,KappaRV,
63       I        myThid)       I        myThid)
64  #endif  #endif
65    
66  #ifdef ALLOW_PP81  #ifdef ALLOW_PP81
67        IF (usePP81) CALL PP81_CALC_VISC(         IF (usePP81) CALL PP81_CALC_VISC(
68       I        bi,bj,iMin+1,iMax,jMin+1,jMax,K,       I        bi,bj, iMin,iMax,jMin,jMax, k,
69       O        KappaRU,KappaRV,       O        KappaRU,KappaRV,
70       I        myThid)       I        myThid)
71  #endif  #endif
72    
73  #ifdef ALLOW_MY82  #ifdef ALLOW_MY82
74        IF (useMY82) CALL MY82_CALC_VISC(         IF (useMY82) CALL MY82_CALC_VISC(
75       I        bi,bj,iMin+1,iMax,jMin+1,jMax,K,       I        bi,bj, iMin,iMax,jMin,jMax, k,
76       O        KappaRU,KappaRV,       O        KappaRU,KappaRV,
77       I        myThid)       I        myThid)
78  #endif  #endif
79    
80  #ifdef ALLOW_GGL90  #ifdef ALLOW_GGL90
81        IF (useGGL90) CALL GGL90_CALC_VISC(         IF (useGGL90) CALL GGL90_CALC_VISC(
82       I        bi,bj,iMin+1,iMax,jMin+1,jMax,K,       I        bi,bj, iMin,iMax,jMin,jMax, k,
83       O        KappaRU,KappaRV,       O        KappaRU,KappaRV,
84       I        myThid)       I        myThid)
85  #endif  #endif
86    
87    C--   end of k loop
88          ENDDO
89    
90        RETURN        RETURN
91        END        END
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22