/[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.3 by cnh, Sun Feb 4 14:38:46 2001 UTC revision 1.4 by jmc, Wed Feb 14 22:51:27 2001 UTC
# Line 19  C     == Global data == Line 19  C     == Global data ==
19  #include "DYNVARS.h"  #include "DYNVARS.h"
20  #include "GRID.h"  #include "GRID.h"
21    
22    #ifdef INCLUDE_DIAGNOSTICS_INTERFACE_CODE
23    #include "AVER.h"
24    #endif
25    
26        EXTERNAL DIFFERENT_MULTIPLE        EXTERNAL DIFFERENT_MULTIPLE
27        LOGICAL  DIFFERENT_MULTIPLE        LOGICAL  DIFFERENT_MULTIPLE
28    
# Line 35  C     myThid - Thread number of this ins Line 39  C     myThid - Thread number of this ins
39  #ifdef INCLUDE_CONVECT_CALL  #ifdef INCLUDE_CONVECT_CALL
40    
41  C     == Local variables ==  C     == Local variables ==
42        INTEGER K        INTEGER i, j, K
43        _RL rhoKm1(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL rhoKm1(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
44        _RL rhoK  (1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL rhoK  (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
45        _RL ConvectCount(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)        _RL ConvectCount(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
# Line 43  C     == Local variables == Line 47  C     == Local variables ==
47  C--   Check to see if should convect now  C--   Check to see if should convect now
48        IF ( DIFFERENT_MULTIPLE(cAdjFreq,myTime,myTime-deltaTClock) ) THEN        IF ( DIFFERENT_MULTIPLE(cAdjFreq,myTime,myTime-deltaTClock) ) THEN
49    
50    C-      Initialisation of Convection Counter
51            DO K=1,Nr
52             DO j=1-OLy,sNy+OLy
53              DO i=1-OLx,sNx+OLx
54               ConvectCount(i,j,k) = 0.
55              ENDDO
56             ENDDO
57            ENDDO
58    
59  #ifdef  ALLOW_KPP  #ifdef  ALLOW_KPP
60          IF (          IF (
61       &      (.NOT.useKPP) ! CONVECT not needed with KPP mixing       &      (.NOT.useKPP) ! CONVECT not needed with KPP mixing
# Line 84  C-          Check static stability with Line 97  C-          Check static stability with
97       U           ConvectCount,       U           ConvectCount,
98       I           myTime,myIter,myThid)       I           myTime,myIter,myThid)
99    
100    #ifdef INCLUDE_DIAGNOSTICS_INTERFACE_CODE
101                IF (myIter.ne.nIter0 .AND. taveFreq.GT.0.) THEN
102                  CALL TIMEAVER_1FLD_XYZ(ConvectCount, ConvectCountTave,
103         I                               deltaTclock, bi, bj, K, myThid)
104                ENDIF
105    #endif /* INCLUDE_DIAGNOSTICS_INTERFACE_CODE */
106    
107  C--       End DO K=1,Nr  C--       End DO K=1,Nr
108            ENDDO            ENDDO
109                
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22