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

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

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

revision 1.19 by jmc, Wed Mar 6 01:59:02 2002 UTC revision 1.20 by heimbach, Sun Mar 24 02:36:39 2002 UTC
# Line 2  C $Header$ Line 2  C $Header$
2  C $Name$  C $Name$
3    
4  #include "CPP_OPTIONS.h"  #include "CPP_OPTIONS.h"
5    #ifdef ALLOW_GMREDI
6    # include "GMREDI_OPTIONS.h"
7    #endif
8    #ifdef ALLOW_KPP
9    # include "KPP_OPTIONS.h"
10    #endif
11    
12  CBOP  CBOP
13  C     !ROUTINE: THERMODYNAMICS  C     !ROUTINE: THERMODYNAMICS
# Line 169  C     uninitialised but inert locations. Line 175  C     uninitialised but inert locations.
175          yA(i,j)      = 0. _d 0          yA(i,j)      = 0. _d 0
176          uTrans(i,j)  = 0. _d 0          uTrans(i,j)  = 0. _d 0
177          vTrans(i,j)  = 0. _d 0          vTrans(i,j)  = 0. _d 0
         DO k=1,Nr  
          phiHyd(i,j,k)  = 0. _d 0  
          sigmaX(i,j,k) = 0. _d 0  
          sigmaY(i,j,k) = 0. _d 0  
          sigmaR(i,j,k) = 0. _d 0  
         ENDDO  
         rhoKM1 (i,j) = 0. _d 0  
178          rhok   (i,j) = 0. _d 0          rhok   (i,j) = 0. _d 0
179          phiSurfX(i,j) = 0. _d 0          phiSurfX(i,j) = 0. _d 0
180          phiSurfY(i,j) = 0. _d 0          phiSurfY(i,j) = 0. _d 0
# Line 223  C--     Set up work arrays that need val Line 222  C--     Set up work arrays that need val
222            fVerS  (i,j,2) = 0. _d 0            fVerS  (i,j,2) = 0. _d 0
223            fVerTr1(i,j,1) = 0. _d 0            fVerTr1(i,j,1) = 0. _d 0
224            fVerTr1(i,j,2) = 0. _d 0            fVerTr1(i,j,2) = 0. _d 0
225              rhoKM1 (i,j)   = 0. _d 0
226           ENDDO           ENDDO
227          ENDDO          ENDDO
228    
# Line 230  C--     Set up work arrays that need val Line 230  C--     Set up work arrays that need val
230           DO j=1-OLy,sNy+OLy           DO j=1-OLy,sNy+OLy
231            DO i=1-OLx,sNx+OLx            DO i=1-OLx,sNx+OLx
232  C This is currently also used by IVDC and Diagnostics  C This is currently also used by IVDC and Diagnostics
233               phiHyd(i,j,k)  = 0. _d 0
234               sigmaX(i,j,k) = 0. _d 0
235               sigmaY(i,j,k) = 0. _d 0
236               sigmaR(i,j,k) = 0. _d 0
237             ConvectCount(i,j,k) = 0.             ConvectCount(i,j,k) = 0.
238             KappaRT(i,j,k) = 0. _d 0             KappaRT(i,j,k) = 0. _d 0
239             KappaRS(i,j,k) = 0. _d 0             KappaRS(i,j,k) = 0. _d 0
# Line 239  C This is currently also used by IVDC an Line 243  C This is currently also used by IVDC an
243  #ifdef ALLOW_PASSIVE_TRACER  #ifdef ALLOW_PASSIVE_TRACER
244             gTr1(i,j,k,bi,bj) = 0. _d 0             gTr1(i,j,k,bi,bj) = 0. _d 0
245  #endif  #endif
246    #ifdef ALLOW_GMREDI
247               Kwx(i,j,k,bi,bj)    = 0. _d 0
248               Kwy(i,j,k,bi,bj)    = 0. _d 0
249               Kwz(i,j,k,bi,bj)    = 0. _d 0
250    #ifdef GM_NON_UNITY_DIAGONAL
251               Kux(i,j,k,bi,bj)    = 0. _d 0
252               Kvy(i,j,k,bi,bj)    = 0. _d 0
253    #endif
254    #endif /* ALLOW_GMREDI */
255  #endif  #endif
256            ENDDO            ENDDO
257           ENDDO           ENDDO
# Line 364  CADJ &     = comlev1_bibj, key=ikey, byt Line 377  CADJ &     = comlev1_bibj, key=ikey, byt
377    
378  #ifdef  ALLOW_GMREDI  #ifdef  ALLOW_GMREDI
379    
 #ifdef ALLOW_AUTODIFF_TAMC  
 CADJ STORE sigmaX(:,:,:) = comlev1, key=ikey, byte=isbyte  
 CADJ STORE sigmaY(:,:,:) = comlev1, key=ikey, byte=isbyte  
 CADJ STORE sigmaR(:,:,:) = comlev1, key=ikey, byte=isbyte  
 #endif /* ALLOW_AUTODIFF_TAMC */  
380  C--     Calculate iso-neutral slopes for the GM/Redi parameterisation  C--     Calculate iso-neutral slopes for the GM/Redi parameterisation
381          IF (useGMRedi) THEN          IF (useGMRedi) THEN
382            CALL GMREDI_CALC_TENSOR(            CALL GMREDI_CALC_TENSOR(
Legend:
Removed from v.1.19  
changed lines
  Added in v.1.20
  ViewVC Help
Powered by ViewVC 1.1.22