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

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

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

revision 1.1 by jmc, Fri Oct 22 16:02:51 2004 UTC revision 1.2 by jmc, Thu Oct 28 00:30:57 2004 UTC
# Line 64  C     msgBuf     :: message buffer Line 64  C     msgBuf     :: message buffer
64        CHARACTER*(MAX_LEN_MBUF) msgBuf        CHARACTER*(MAX_LEN_MBUF) msgBuf
65  CEOP  CEOP
66    
       KbryanLewis79=diffKrBL79surf+(diffKrBL79deep-diffKrBL79surf)  
      & *( atan( -( rF(k)-diffKrBL79Ho )/diffKrBL79scl )/PI+0.5 _d 0)  
   
67        IF ( trIdentity.EQ.GAD_TEMPERATURE) THEN        IF ( trIdentity.EQ.GAD_TEMPERATURE) THEN
68    
69         DO k = 1,Nr         DO k = 1,Nr
70            KbryanLewis79=diffKrBL79surf+(diffKrBL79deep-diffKrBL79surf)
71         &  *( atan( -( rF(k)-diffKrBL79Ho )/diffKrBL79scl )/PI+0.5 _d 0)
72          DO j = 1-Oly,sNy+Oly          DO j = 1-Oly,sNy+Oly
73           DO i = 1-Olx,sNx+Olx           DO i = 1-Olx,sNx+Olx
74            KappaRTr(i,j,k) =            KappaRTr(i,j,k) =
# Line 87  CEOP Line 86  CEOP
86        ELSEIF ( trIdentity.EQ.GAD_SALINITY) THEN        ELSEIF ( trIdentity.EQ.GAD_SALINITY) THEN
87    
88         DO k = 1,Nr         DO k = 1,Nr
89            KbryanLewis79=diffKrBL79surf+(diffKrBL79deep-diffKrBL79surf)
90         &  *( atan( -( rF(k)-diffKrBL79Ho )/diffKrBL79scl )/PI+0.5 _d 0)
91          DO j = 1-Oly, sNy+Oly          DO j = 1-Oly, sNy+Oly
92           DO i = 1-Olx, sNx+Olx           DO i = 1-Olx, sNx+Olx
93            KappaRTr(i,j,k) =            KappaRTr(i,j,k) =
# Line 107  CEOP Line 108  CEOP
108    
109         iTr = trIdentity - GAD_TR1 + 1         iTr = trIdentity - GAD_TR1 + 1
110         DO k = 1,Nr         DO k = 1,Nr
111            KbryanLewis79=diffKrBL79surf+(diffKrBL79deep-diffKrBL79surf)
112         &  *( atan( -( rF(k)-diffKrBL79Ho )/diffKrBL79scl )/PI+0.5 _d 0)
113          DO j = 1-Oly, sNy+Oly          DO j = 1-Oly, sNy+Oly
114           DO i = 1-Olx, sNx+Olx           DO i = 1-Olx, sNx+Olx
115            KappaRTr(i,j,k) =            KappaRTr(i,j,k) =
# Line 144  C--   Add physical pacakge contributions Line 147  C--   Add physical pacakge contributions
147        IF (trUseKPP) THEN        IF (trUseKPP) THEN
148         IF (trIdentity.EQ.GAD_TEMPERATURE) THEN         IF (trIdentity.EQ.GAD_TEMPERATURE) THEN
149           CALL KPP_CALC_DIFF_T(           CALL KPP_CALC_DIFF_T(
150       I        bi,bj,iMin+1,iMax,jMin+1,jMax,0,Nr,       I        bi,bj,iMin,iMax,jMin,jMax,0,Nr,
151       U        KappaRTr,       U        KappaRTr,
152       I        myThid)       I        myThid)
153         ELSE         ELSE
154           CALL KPP_CALC_DIFF_S(           CALL KPP_CALC_DIFF_S(
155       I        bi,bj,iMin+1,iMax,jMin+1,jMax,0,Nr,       I        bi,bj,iMin,iMax,jMin,jMax,0,Nr,
156       U        KappaRTr,       U        KappaRTr,
157       I        myThid)       I        myThid)
158         ENDIF         ENDIF
# Line 159  C--   Add physical pacakge contributions Line 162  C--   Add physical pacakge contributions
162  #ifdef ALLOW_PP81  #ifdef ALLOW_PP81
163        IF (usePP81) THEN        IF (usePP81) THEN
164           CALL PP81_CALC_DIFF(           CALL PP81_CALC_DIFF(
165       I        bi,bj,iMin+1,iMax,jMin+1,jMax,0,Nr,       I        bi,bj,iMin,iMax,jMin,jMax,0,Nr,
166       U        KappaRTr,       U        KappaRTr,
167       I        myThid)       I        myThid)
168        ENDIF        ENDIF
# Line 168  C--   Add physical pacakge contributions Line 171  C--   Add physical pacakge contributions
171  #ifdef ALLOW_MY82  #ifdef ALLOW_MY82
172        IF (useMY82) THEN        IF (useMY82) THEN
173           CALL MY82_CALC_DIFF(           CALL MY82_CALC_DIFF(
174       I        bi,bj,iMin+1,iMax,jMin+1,jMax,0,Nr,       I        bi,bj,iMin,iMax,jMin,jMax,0,Nr,
175       U        KappaRTr,       U        KappaRTr,
176       I        myThid)       I        myThid)
177        ENDIF        ENDIF
# Line 177  C--   Add physical pacakge contributions Line 180  C--   Add physical pacakge contributions
180  #ifdef ALLOW_GGL90  #ifdef ALLOW_GGL90
181        IF (useGGL90) THEN        IF (useGGL90) THEN
182           CALL GGL90_CALC_DIFF(           CALL GGL90_CALC_DIFF(
183       I        bi,bj,iMin+1,iMax,jMin+1,jMax,0,Nr,       I        bi,bj,iMin,iMax,jMin,jMax,0,Nr,
184       O        KappaRTr,       O        KappaRTr,
185       I        myThid)       I        myThid)
186        ENDIF        ENDIF
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22