/[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.12 by cnh, Sun Feb 4 14:38:45 2001 UTC revision 1.13 by adcroft, Tue May 29 14:01:36 2001 UTC
# Line 6  C $Name$ Line 6  C $Name$
6  CStartOfInterFace  CStartOfInterFace
7        SUBROUTINE CALC_DIFFUSIVITY(        SUBROUTINE CALC_DIFFUSIVITY(
8       I        bi,bj,iMin,iMax,jMin,jMax,K,       I        bi,bj,iMin,iMax,jMin,jMax,K,
9       I        maskC,maskUp,       I        maskUp,
10       O        KappaRT,KappaRS,KappaRU,KappaRV,       O        KappaRT,KappaRS,KappaRU,KappaRV,
11       I        myThid)       I        myThid)
12    
# Line 25  C     == GLobal variables == Line 25  C     == GLobal variables ==
25    
26  C     == Routine arguments ==  C     == Routine arguments ==
27  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     maskC   - land/water mask for tracer points  
28  C     maskUp  - land/water mask for Wvel points (above tracer level)  C     maskUp  - land/water mask for Wvel points (above tracer level)
29  C     myThid - Instance number for this innvocation of CALC_COMMON_FACTORS  C     myThid - Instance number for this innvocation of CALC_COMMON_FACTORS
30  C  C
31        INTEGER bi,bj,iMin,iMax,jMin,jMax,K        INTEGER bi,bj,iMin,iMax,jMin,jMax,K
       _RS maskC(1-Olx:sNx+Olx,1-Oly:sNy+Oly)  
32        _RS maskUp(1-Olx:sNx+Olx,1-Oly:sNy+Oly)        _RS maskUp(1-Olx:sNx+Olx,1-Oly:sNy+Oly)
33        _RL KappaRT(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)        _RL KappaRT(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)
34        _RL KappaRS(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)        _RL KappaRS(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)
# Line 45  C     I, J, K - Loop counters Line 43  C     I, J, K - Loop counters
43    
44        DO j=jMin,jMax        DO j=jMin,jMax
45         DO i=iMin,iMax         DO i=iMin,iMax
46          KappaRT(i,j,k) = maskC(i,j)*maskUp(i,j)*          KappaRT(i,j,k) = maskUp(i,j)*
47       &      (       &      (
48       &         KappaRT(i,j,k)       &         KappaRT(i,j,k)
49       &       + diffKrT       &       + diffKrT
# Line 55  C     I, J, K - Loop counters Line 53  C     I, J, K - Loop counters
53    
54        DO j=jMin,jMax        DO j=jMin,jMax
55         DO i=iMin,iMax         DO i=iMin,iMax
56          KappaRS(i,j,k) = maskC(i,j)*maskUp(i,j)*          KappaRS(i,j,k) = maskUp(i,j)*
57       &      (       &      (
58       &         KappaRS(i,j,k)       &         KappaRS(i,j,k)
59       &       + diffKrS       &       + diffKrS
Legend:
Removed from v.1.12  
changed lines
  Added in v.1.13
  ViewVC Help
Powered by ViewVC 1.1.22