/[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.16 by jmc, Tue Mar 16 00:08:27 2010 UTC revision 1.17 by jmc, Tue Sep 20 22:25:23 2011 UTC
# Line 112  CEOP Line 112  CEOP
112           DO j = 1-Oly,sNy+Oly           DO j = 1-Oly,sNy+Oly
113            DO i = 1-Olx,sNx+Olx            DO i = 1-Olx,sNx+Olx
114             KappaRTr(i,j,k) = KappaRTr(i,j,k)             KappaRTr(i,j,k) = KappaRTr(i,j,k)
115  #if (defined ALLOW_3D_DIFFKR || \  #ifdef ALLOW_3D_DIFFKR
      (defined (ALLOW_AUTODIFF_TAMC) && defined (ALLOW_DIFFKR_CONTROL)))  
116       &          + diffKr(i,j,k,bi,bj)       &          + diffKr(i,j,k,bi,bj)
117  #else  #else
118       &          + diffKrNrT(k)       &          + diffKrNrT(k)
# Line 126  CEOP Line 125  CEOP
125           DO j = 1-Oly, sNy+Oly           DO j = 1-Oly, sNy+Oly
126            DO i = 1-Olx, sNx+Olx            DO i = 1-Olx, sNx+Olx
127             KappaRTr(i,j,k) = KappaRTr(i,j,k)             KappaRTr(i,j,k) = KappaRTr(i,j,k)
128  #if (defined ALLOW_3D_DIFFKR || \  #ifdef ALLOW_3D_DIFFKR
      (defined (ALLOW_AUTODIFF_TAMC) && defined (ALLOW_DIFFKR_CONTROL)))  
129       &          + diffKr(i,j,k,bi,bj)       &          + diffKr(i,j,k,bi,bj)
130  #else  #else
131       &          + diffKrNrS(k)       &          + diffKrNrS(k)
# Line 143  CEOP Line 141  CEOP
141           DO j = 1-Oly, sNy+Oly           DO j = 1-Oly, sNy+Oly
142            DO i = 1-Olx, sNx+Olx            DO i = 1-Olx, sNx+Olx
143             KappaRTr(i,j,k) = KappaRTr(i,j,k)             KappaRTr(i,j,k) = KappaRTr(i,j,k)
144  #if (defined ALLOW_3D_DIFFKR || \  #ifdef ALLOW_3D_DIFFKR
      (defined (ALLOW_AUTODIFF_TAMC) && defined (ALLOW_DIFFKR_CONTROL)))  
145       &          + diffKr(i,j,k,bi,bj)       &          + diffKr(i,j,k,bi,bj)
146  #else  #else
147       &          + PTRACERS_diffKrNr(k,iTr)       &          + PTRACERS_diffKrNr(k,iTr)
# Line 176  C--   Set vertical diffusivity contribut Line 173  C--   Set vertical diffusivity contribut
173       I        bi,bj,iMin,iMax,jMin,jMax,0,Nr,       I        bi,bj,iMin,iMax,jMin,jMax,0,Nr,
174       O        KappaRTr,       O        KappaRTr,
175       I        myThid)       I        myThid)
176    #ifdef ALLOW_PTRACERS
177         ELSEIF ( trIdentity.GE.GAD_TR1) THEN         ELSEIF ( trIdentity.GE.GAD_TR1) THEN
178             iTr = trIdentity - GAD_TR1 + 1
179           CALL KPP_CALC_DIFF_Ptr(           CALL KPP_CALC_DIFF_Ptr(
180       I        bi,bj,iMin,iMax,jMin,jMax,0,Nr,       I        bi,bj,iMin,iMax,jMin,jMax,0,Nr,
181       O        KappaRTr,       O        KappaRTr,
182       I        myThid)       I        iTr, myThid )
183         ENDIF  #endif /* ALLOW_PTRACERS */
184  #if (defined ALLOW_PTRACERS && ! (defined ALLOW_3D_DIFFKR || \         ELSE
185       (defined (ALLOW_AUTODIFF_TAMC) && defined (ALLOW_DIFFKR_CONTROL))))          WRITE(msgBuf,'(A,I4)')
186         IF ( trIdentity.GE.GAD_TR1) THEN       &       ' CALC_3D_DIFFUSIVITY: Invalid tracer Id: ',trIdentity
187          iTr = trIdentity - GAD_TR1 + 1          CALL PRINT_ERROR( msgBuf, myThid )
188          DO k = 1,Nr          STOP 'ABNORMAL END: S/R CALC_3D_DIFFUSIVITY'
          DO j = 1-Oly, sNy+Oly  
           DO i = 1-Olx, sNx+Olx  
            KappaRTr(i,j,k) = KappaRTr(i,j,k)  
      &          - diffKrNrS(k) + PTRACERS_diffKrNr(k,iTr)  
           ENDDO  
          ENDDO  
         ENDDO  
189         ENDIF         ENDIF
 #endif  
190        ENDIF        ENDIF
191  #endif /* ALLOW_KPP */  #endif /* ALLOW_KPP */
192    
Legend:
Removed from v.1.16  
changed lines
  Added in v.1.17
  ViewVC Help
Powered by ViewVC 1.1.22