/[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.5 by jmc, Mon Nov 5 18:48:04 2007 UTC revision 1.6 by dfer, Tue Feb 12 23:02:47 2008 UTC
# Line 33  C     == GLobal variables == Line 33  C     == GLobal variables ==
33  #include "PTRACERS_SIZE.h"  #include "PTRACERS_SIZE.h"
34  #include "PTRACERS_PARAMS.h"  #include "PTRACERS_PARAMS.h"
35  c #include "PTRACERS_FIELDS.h"  c #include "PTRACERS_FIELDS.h"
36  c #include "GAD.h"  #include "GAD.h"
37    
38  C     !INPUT/OUTPUT PARAMETERS:  C     !INPUT/OUTPUT PARAMETERS:
39  C     == Routine arguments ==  C     == Routine arguments ==
# Line 52  C     KappaRS   :: Net diffusivity for s Line 52  C     KappaRS   :: Net diffusivity for s
52  C     !LOCAL VARIABLES:  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,trIdentity
56        _RL KbryanLewis79        _RL KbryanLewis79
57  #ifdef ALLOW_BL79_LAT_VARY  #ifdef ALLOW_BL79_LAT_VARY
58        _RL KbryanLewisEQ        _RL KbryanLewisEQ
# Line 69  CEOP Line 69  CEOP
69  C Loop over tracers  C Loop over tracers
70        DO iTr=1,PTRACERS_numInUse        DO iTr=1,PTRACERS_numInUse
71    
72           trIdentity = iTr + GAD_TR1 - 1
73    
74         IF ( .NOT. PTRACERS_useKPP(iTr) ) THEN         IF ( .NOT. PTRACERS_useKPP(iTr) ) THEN
75          DO j = 1-Oly, sNy+Oly          DO j = 1-Oly, sNy+Oly
76           DO i = 1-Olx, sNx+Olx           DO i = 1-Olx, sNx+Olx
# Line 103  C--   Set vertical diffusivity contribut Line 105  C--   Set vertical diffusivity contribut
105           CALL GMREDI_CALC_DIFF(           CALL GMREDI_CALC_DIFF(
106       I        bi,bj,iMin,iMax,jMin,jMax,k,1,       I        bi,bj,iMin,iMax,jMin,jMax,k,1,
107       U        KappaRtr(1-Olx,1-Oly,iTr),       U        KappaRtr(1-Olx,1-Oly,iTr),
108       I        myThid)       I        trIdentity,myThid)
109          ENDIF          ENDIF
110  #endif  #endif
111    
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22