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

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

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

revision 1.32 by dimitri, Mon Apr 23 20:46:49 2007 UTC revision 1.33 by jmc, Tue Sep 4 16:49:44 2007 UTC
# Line 6  C $Name$ Line 6  C $Name$
6    
7  CBOP  CBOP
8  C     !ROUTINE: CALC_DIFFUSIVITY  C     !ROUTINE: CALC_DIFFUSIVITY
9  C     !INTERFACE:  C     !INTERFACE:
10        SUBROUTINE CALC_DIFFUSIVITY(        SUBROUTINE CALC_DIFFUSIVITY(
11       I        bi,bj,iMin,iMax,jMin,jMax,k,       I        bi,bj,iMin,iMax,jMin,jMax,k,
12       I        maskUp,       I        maskUp,
13       O        KappaRT,KappaRS,       O        KappaRT,KappaRS,
# Line 15  C     !INTERFACE: Line 15  C     !INTERFACE:
15    
16  C     !DESCRIPTION: \bv  C     !DESCRIPTION: \bv
17  C     *==========================================================*  C     *==========================================================*
18  C     | SUBROUTINE CALC_DIFFUSIVITY                                C     | SUBROUTINE CALC_DIFFUSIVITY
19  C     | o Calculate net vertical diffusivity                                C     | o Calculate net vertical diffusivity
20  C     *==========================================================*  C     *==========================================================*
21  C     | Combines spatially varying diffusion coefficients from  C     | Combines spatially varying diffusion coefficients from
22  C     | KPP and/or GM and/or convective stability test.  C     | KPP and/or GM and/or convective stability test.
# Line 36  C     !INPUT/OUTPUT PARAMETERS: Line 36  C     !INPUT/OUTPUT PARAMETERS:
36  C     == Routine arguments ==  C     == Routine arguments ==
37  C     bi, bj, :: tile indices  C     bi, bj, :: tile indices
38  C     iMin,   :: Range of points for which calculation is performed.  C     iMin,   :: Range of points for which calculation is performed.
39  C     iMax,  C     iMax,
40  C     jMin,  C     jMin,
41  C     jMax  C     jMax
42  C     maskUp  :: land/water mask for Wvel points (above tracer level)  C     maskUp  :: land/water mask for Wvel points (above tracer level)
43  C     myThid  :: Instance number for this innvocation of CALC_DIFFUSIVITY  C     myThid  :: Instance number for this innvocation of CALC_DIFFUSIVITY
# Line 53  C     !LOCAL VARIABLES: Line 53  C     !LOCAL VARIABLES:
53  C     == Local variables ==  C     == Local variables ==
54  C     I, J :: Loop counters  C     I, J :: Loop counters
55        INTEGER i,j        INTEGER i,j
56        _RL KbryanLewis79, KbryanLewisEQ        _RL KbryanLewis79
57    #ifdef ALLOW_BL79_LAT_VARY
58          _RL KbryanLewisEQ
59    #endif
60  CEOP  CEOP
61    
62        IF ( .NOT. UseKPP ) THEN        IF ( .NOT. UseKPP ) THEN
# Line 100  CEOP Line 103  CEOP
103        ENDIF        ENDIF
104    
105  #ifdef ALLOW_GMREDI  #ifdef ALLOW_GMREDI
106        IF (useGMRedi) THEN        IF (useGMRedi) THEN
107           CALL GMREDI_CALC_DIFF(           CALL GMREDI_CALC_DIFF(
108       I        bi,bj,iMin,iMax,jMin,jMax,k,1,       I        bi,bj,iMin,iMax,jMin,jMax,k,1,
109       U        KappaRT,       U        KappaRT,
# Line 150  CEOP Line 153  CEOP
153       I        myThid)       I        myThid)
154        ENDIF        ENDIF
155  #endif  #endif
156          
157  #ifdef ALLOW_GGL90  #ifdef ALLOW_GGL90
158        IF (useGGL90) THEN        IF (useGGL90) THEN
159           CALL GGL90_CALC_DIFF(           CALL GGL90_CALC_DIFF(
# Line 163  CEOP Line 166  CEOP
166       I        myThid)       I        myThid)
167        ENDIF        ENDIF
168  #endif  #endif
169          
170  C-    Apply mask to vertical diffusivity  C-    Apply mask to vertical diffusivity
171  C jmc: don't have the impression that masking is needed  C jmc: don't have the impression that masking is needed
172  C      but could be removed later if it's the case.  C      but could be removed later if it's the case.
173        DO j = 1-Oly, sNy+Oly        DO j = 1-Oly, sNy+Oly
174         DO i = 1-Olx, sNx+Olx         DO i = 1-Olx, sNx+Olx
Legend:
Removed from v.1.32  
changed lines
  Added in v.1.33
  ViewVC Help
Powered by ViewVC 1.1.22