/[MITgcm]/MITgcm/pkg/ptracers/ptracers_calc_diff.F
ViewVC logotype

Diff of /MITgcm/pkg/ptracers/ptracers_calc_diff.F

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

revision 1.4 by dimitri, Thu Sep 27 14:42:59 2007 UTC revision 1.5 by jmc, Mon Nov 5 18:48:04 2007 UTC
# Line 5  C $Name$ Line 5  C $Name$
5    
6  CBOP  CBOP
7  C     !ROUTINE: PTRACERS_CALC_DIFF  C     !ROUTINE: PTRACERS_CALC_DIFF
8  C     !INTERFACE:  C     !INTERFACE:
9        SUBROUTINE PTRACERS_CALC_DIFF(        SUBROUTINE PTRACERS_CALC_DIFF(
10       I        bi,bj,iMin,iMax,jMin,jMax,k,       I        bi,bj,iMin,iMax,jMin,jMax,k,
11       I        maskUp,       I        maskUp,
12       O        KappaRtr,       O        KappaRtr,
# Line 31  C     == GLobal variables == Line 31  C     == GLobal variables ==
31  #include "GRID.h"  #include "GRID.h"
32  #include "DYNVARS.h"  #include "DYNVARS.h"
33  #include "PTRACERS_SIZE.h"  #include "PTRACERS_SIZE.h"
34  #include "PTRACERS.h"  #include "PTRACERS_PARAMS.h"
35    c #include "PTRACERS_FIELDS.h"
36  c #include "GAD.h"  c #include "GAD.h"
37    
38  C     !INPUT/OUTPUT PARAMETERS:  C     !INPUT/OUTPUT PARAMETERS:
# Line 52  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,iTr        INTEGER i,j,iTr
56        _RL KbryanLewis79, KbryanLewisEQ        _RL KbryanLewis79
57    #ifdef ALLOW_BL79_LAT_VARY
58          _RL KbryanLewisEQ
59    #endif
60  CEOP  CEOP
61    
62        KbryanLewis79=diffKrBL79surf+(diffKrBL79deep-diffKrBL79surf)        KbryanLewis79=diffKrBL79surf+(diffKrBL79deep-diffKrBL79surf)
# Line 95  C--   Set vertical diffusivity contribut Line 99  C--   Set vertical diffusivity contribut
99  #endif  #endif
100    
101  #ifdef ALLOW_GMREDI  #ifdef ALLOW_GMREDI
102          IF ( PTRACERS_useGMRedi(iTr) ) THEN          IF ( PTRACERS_useGMRedi(iTr) ) THEN
103           CALL GMREDI_CALC_DIFF(           CALL GMREDI_CALC_DIFF(
104       I        bi,bj,iMin,iMax,jMin,jMax,k,1,       I        bi,bj,iMin,iMax,jMin,jMax,k,1,
105       U        KappaRtr(1-Olx,1-Oly,iTr),       U        KappaRtr(1-Olx,1-Oly,iTr),
# Line 120  C--   Set vertical diffusivity contribut Line 124  C--   Set vertical diffusivity contribut
124       I        myThid)       I        myThid)
125          ENDIF          ENDIF
126  #endif  #endif
127          
128  #ifdef ALLOW_GGL90  #ifdef ALLOW_GGL90
129          IF (useGGL90) THEN          IF (useGGL90) THEN
130           CALL GGL90_CALC_DIFF(           CALL GGL90_CALC_DIFF(
# Line 129  C--   Set vertical diffusivity contribut Line 133  C--   Set vertical diffusivity contribut
133       I        myThid)       I        myThid)
134          ENDIF          ENDIF
135  #endif  #endif
136          
137  C-    Apply mask to vertical diffusivity  C-    Apply mask to vertical diffusivity
138  C jmc: don't have the impression that masking is needed  C jmc: don't have the impression that masking is needed
139  C      but could be removed later if it's the case.  C      but could be removed later if it's the case.
140          DO j = 1-Oly, sNy+Oly          DO j = 1-Oly, sNy+Oly
141           DO i = 1-Olx, sNx+Olx           DO i = 1-Olx, sNx+Olx
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22