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

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

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

revision 1.14 by mlosch, Wed Sep 18 16:38:01 2002 UTC revision 1.15 by mlosch, Fri Nov 1 22:00:33 2002 UTC
# Line 54  C     !LOCAL VARIABLES: Line 54  C     !LOCAL VARIABLES:
54  C     == Local variables ==  C     == Local variables ==
55  C     rhoKm1, rhoK :: Density at adjacent levels (common ref. level)  C     rhoKm1, rhoK :: Density at adjacent levels (common ref. level)
56  C     ConvectCount :: Convection mixing freq. counter.  C     ConvectCount :: Convection mixing freq. counter.
57        INTEGER i, j, K        INTEGER i, j, K, kTop, kBottom, kDir, kRef
58        _RL rhoKm1(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL rhoKm1(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
59        _RL rhoK  (1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL rhoK  (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
60        _RL ConvectCount(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)        _RL ConvectCount(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
# Line 94  C--       CONVECT not needed with KPP mi Line 94  C--       CONVECT not needed with KPP mi
94       &     ) THEN       &     ) THEN
95  #endif  /* ALLOW_KPP */  #endif  /* ALLOW_KPP */
96    
97             IF ( buoyancyRelation .eq. 'OCEANIC' ) THEN
98              kTop    = 2
99              kBottom = Nr
100              kDir    = 1
101             ELSEIF ( buoyancyRelation .eq. 'OCEANICP' ) THEN
102              kTop    = Nr-1
103              kBottom =  1
104              kDir    = -1
105             ELSEIF ( buoyancyRelation .eq. 'ATMOSPHERIC' ) THEN
106              kTop    = Nr-1
107              kBottom =  1
108              kDir    = -1
109             ELSE
110              STOP 'CONVECTIVE_ADJUSTMENT: We should never reach this point'
111             ENDIF  
112    
113  C--       Loop over all *interior* layers  C--       Loop over all *interior* layers
114            DO K=2,Nr            DO K=kTop,kBottom,kDir
115    
116               IF ( buoyancyRelation .eq. 'OCEANIC' ) THEN
117                kRef = k-1
118               ELSEIF ( buoyancyRelation .eq. 'OCEANICP' ) THEN
119                kRef = k
120               ELSEIF ( buoyancyRelation .eq. 'ATMOSPHERIC' ) THEN
121                kRef = k
122               ENDIF        
123    
124  #ifdef ALLOW_AUTODIFF_TAMC  #ifdef ALLOW_AUTODIFF_TAMC
125              kkey = (ikey-1)*Nr + k              kkey = (ikey-1)*Nr + k
# Line 104  CADJ STORE salt (:,:,k-1,bi,bj) = comlev Line 128  CADJ STORE salt (:,:,k-1,bi,bj) = comlev
128  #endif /* ALLOW_AUTODIFF_TAMC */  #endif /* ALLOW_AUTODIFF_TAMC */
129  C-          Density of K-1 layer (above W(K)) reference to K-1 T-level  C-          Density of K-1 layer (above W(K)) reference to K-1 T-level
130              CALL FIND_RHO(              CALL FIND_RHO(
131       I           bi,bj,iMin,iMax,jMin,jMax,K-1,K-1,       I           bi,bj,iMin,iMax,jMin,jMax,K-1,KRef,
132       I           theta,salt,       I           theta,salt,
133       O           rhoKm1,       O           rhoKm1,
134       I           myThid )       I           myThid )
# Line 115  CADJ STORE theta(:,:,k,bi,bj) = comlev1_ Line 139  CADJ STORE theta(:,:,k,bi,bj) = comlev1_
139  CADJ STORE salt (:,:,k,bi,bj) = comlev1_bibj_k, key=kkey, byte=isbyte  CADJ STORE salt (:,:,k,bi,bj) = comlev1_bibj_k, key=kkey, byte=isbyte
140  #endif /* ALLOW_AUTODIFF_TAMC */  #endif /* ALLOW_AUTODIFF_TAMC */
141              CALL FIND_RHO(              CALL FIND_RHO(
142       I           bi,bj,iMin,iMax,jMin,jMax,K,K-1,       I           bi,bj,iMin,iMax,jMin,jMax,K,KRef,
143       I           theta,salt,       I           theta,salt,
144       O           rhoK,       O           rhoK,
145       I           myThid )       I           myThid )
Legend:
Removed from v.1.14  
changed lines
  Added in v.1.15
  ViewVC Help
Powered by ViewVC 1.1.22