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

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

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

Revision 1.17 - (hide annotations) (download)
Fri Nov 15 03:01:21 2002 UTC (21 years, 6 months ago) by heimbach
Branch: MAIN
CVS Tags: checkpoint47e_post, checkpoint47c_post, checkpoint48e_post, checkpoint50c_pre, checkpoint48i_post, checkpoint50, checkpoint50b_pre, checkpoint48b_post, checkpoint48c_pre, checkpoint47d_pre, checkpoint47a_post, checkpoint48d_pre, checkpoint47i_post, checkpoint47d_post, checkpoint48d_post, checkpoint48f_post, checkpoint48h_post, checkpoint47g_post, checkpoint48a_post, checkpoint50a_post, checkpoint47j_post, branch-exfmods-tag, checkpoint48c_post, checkpoint47b_post, checkpoint47f_post, checkpoint47, checkpoint48, checkpoint49, checkpoint48g_post, checkpoint47h_post, checkpoint50b_post
Branch point for: branch-exfmods-curt
Changes since 1.16: +15 -14 lines 
differentiable version of checkpoint46n_post
o external_fields_load now part of differentiation list
o pressure needs multiple storing;
  would be nice to have store_pressure at beginning or
  end of forward_step, e.g. by having phiHyd global (5-dim.)
  (NB: pressure is needed for certain cases in find_rho,
  which is also invoked through convective_adjustment).
o recomputations in find_rho for cases
 'JMD95'/'UNESCO' or 'MDJWF' are OK.
o #define ATMOSPHERIC_LOADING should be differentiable
o ini_forcing shifted to begining of initialise_varia
1 heimbach 1.17 C $Header: /u/gcmpack/MITgcm/model/src/convective_adjustment.F,v 1.16 2002/11/02 20:07:39 mlosch Exp $
2 cnh 1.8 C $Name: $
3 adcroft 1.2
4     #include "CPP_OPTIONS.h"
5    
6 cnh 1.8 CBOP
7     C !ROUTINE: CONVECTIVE_ADJUSTMENT
8     C !INTERFACE:
9 adcroft 1.2 SUBROUTINE CONVECTIVE_ADJUSTMENT(
10     I bi, bj, iMin, iMax, jMin, jMax,
11     I myTime, myIter, myThid )
12 cnh 1.8 C !DESCRIPTION: \bv
13     C *==========================================================*
14     C | SUBROUTINE CONVECTIVE_ADJUSTMENT
15     C | o Driver for vertical mixing or similar parameterization
16     C *==========================================================*
17     C \ev
18    
19     C !USES:
20 adcroft 1.2 IMPLICIT NONE
21     C == Global data ==
22     #include "SIZE.h"
23     #include "EEPARAMS.h"
24     #include "PARAMS.h"
25     #include "DYNVARS.h"
26     #include "GRID.h"
27 heimbach 1.12 #ifdef ALLOW_PASSIVE_TRACER
28     #include "TR1.h"
29     #endif
30 jmc 1.5 #ifdef ALLOW_TIMEAVE
31     #include "TIMEAVE_STATV.h"
32 jmc 1.4 #endif
33 heimbach 1.6 #ifdef ALLOW_AUTODIFF_TAMC
34     #include "tamc.h"
35     #include "tamc_keys.h"
36     #endif /* ALLOW_AUTODIFF_TAMC */
37 adcroft 1.2 EXTERNAL DIFFERENT_MULTIPLE
38     LOGICAL DIFFERENT_MULTIPLE
39    
40 cnh 1.8 C !INPUT/OUTPUT PARAMETERS:
41 adcroft 1.2 C == Routine arguments ==
42     C bi,bj,iMin,iMax,jMin,jMax,K - Loop counters
43     C myTime - Current time in simulation
44     C myIter - Current iteration in simulation
45     C myThid - Thread number of this instance of S/R CONVECT
46     INTEGER bi,bj,iMin,iMax,jMin,jMax
47     _RL myTime
48     INTEGER myIter
49     INTEGER myThid
50    
51     #ifdef INCLUDE_CONVECT_CALL
52    
53 cnh 1.8 C !LOCAL VARIABLES:
54 adcroft 1.2 C == Local variables ==
55 cnh 1.8 C rhoKm1, rhoK :: Density at adjacent levels (common ref. level)
56     C ConvectCount :: Convection mixing freq. counter.
57 heimbach 1.17 INTEGER i, j, K, kTop, kBottom, kDir, deltaK
58 adcroft 1.2 _RL rhoKm1(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
59     _RL rhoK (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
60     _RL ConvectCount(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
61 adcroft 1.10 _RL weightA(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
62     _RL weightB(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
63 cnh 1.8 CEOP
64 adcroft 1.2
65     C-- Check to see if should convect now
66     IF ( DIFFERENT_MULTIPLE(cAdjFreq,myTime,myTime-deltaTClock) ) THEN
67    
68 heimbach 1.17 C-- Initialise counters
69     kTop = 0
70     kBottom = 0
71     kDir = 0
72     deltaK = 0
73    
74 jmc 1.4 C- Initialisation of Convection Counter
75     DO K=1,Nr
76     DO j=1-OLy,sNy+OLy
77     DO i=1-OLx,sNx+OLx
78     ConvectCount(i,j,k) = 0.
79     ENDDO
80     ENDDO
81     ENDDO
82    
83 heimbach 1.6 #ifdef ALLOW_AUTODIFF_TAMC
84     act1 = bi - myBxLo(myThid)
85     max1 = myBxHi(myThid) - myBxLo(myThid) + 1
86     act2 = bj - myByLo(myThid)
87     max2 = myByHi(myThid) - myByLo(myThid) + 1
88     act3 = myThid - 1
89     max3 = nTx*nTy
90     act4 = ikey_dynamics - 1
91     ikey = (act1 + 1) + act2*max1
92     & + act3*max1*max2
93     & + act4*max1*max2*max3
94     #endif /* ALLOW_AUTODIFF_TAMC */
95    
96 adcroft 1.2 #ifdef ALLOW_KPP
97 heimbach 1.13 C-- CONVECT not needed with KPP mixing
98 adcroft 1.2 IF (
99 heimbach 1.13 & (.NOT.useKPP)
100 adcroft 1.2 & ) THEN
101     #endif /* ALLOW_KPP */
102    
103 mlosch 1.15 IF ( buoyancyRelation .eq. 'OCEANIC' ) THEN
104 heimbach 1.17 kTop = 2
105 mlosch 1.15 kBottom = Nr
106 heimbach 1.17 kDir = 1
107     deltaK = -1
108 mlosch 1.15 ELSEIF ( buoyancyRelation .eq. 'OCEANICP' ) THEN
109 mlosch 1.16 kTop = Nr
110     kBottom = 2
111 mlosch 1.15 kDir = -1
112 heimbach 1.17 deltaK = 0
113 mlosch 1.15 ELSEIF ( buoyancyRelation .eq. 'ATMOSPHERIC' ) THEN
114 mlosch 1.16 kTop = Nr
115     kBottom = 2
116 mlosch 1.15 kDir = -1
117 heimbach 1.17 deltaK = 0
118 mlosch 1.15 ELSE
119     STOP 'CONVECTIVE_ADJUSTMENT: We should never reach this point'
120     ENDIF
121    
122 adcroft 1.2 C-- Loop over all *interior* layers
123 mlosch 1.15 DO K=kTop,kBottom,kDir
124    
125 adcroft 1.2 #ifdef ALLOW_AUTODIFF_TAMC
126 heimbach 1.7 kkey = (ikey-1)*Nr + k
127     CADJ STORE theta(:,:,k-1,bi,bj) = comlev1_bibj_k, key=kkey, byte=isbyte
128     CADJ STORE salt (:,:,k-1,bi,bj) = comlev1_bibj_k, key=kkey, byte=isbyte
129 adcroft 1.2 #endif /* ALLOW_AUTODIFF_TAMC */
130     C- Density of K-1 layer (above W(K)) reference to K-1 T-level
131     CALL FIND_RHO(
132 heimbach 1.17 I bi,bj,iMin,iMax,jMin,jMax,K-1,K+deltaK,
133 adcroft 1.2 I theta,salt,
134     O rhoKm1,
135     I myThid )
136    
137     C- Density of K layer (below W(K)) reference to K-1 T-level.
138     #ifdef ALLOW_AUTODIFF_TAMC
139 heimbach 1.7 CADJ STORE theta(:,:,k,bi,bj) = comlev1_bibj_k, key=kkey, byte=isbyte
140     CADJ STORE salt (:,:,k,bi,bj) = comlev1_bibj_k, key=kkey, byte=isbyte
141 adcroft 1.2 #endif /* ALLOW_AUTODIFF_TAMC */
142     CALL FIND_RHO(
143 heimbach 1.17 I bi,bj,iMin,iMax,jMin,jMax,K,K+deltaK,
144 adcroft 1.2 I theta,salt,
145     O rhoK,
146     I myThid )
147    
148     #ifdef ALLOW_AUTODIFF_TAMC
149 heimbach 1.7 CADJ STORE rhoKm1(:,:) = comlev1_bibj_k, key = kkey, byte = isbyte
150     CADJ STORE rhoK (:,:) = comlev1_bibj_k, key = kkey, byte = isbyte
151 adcroft 1.2 #endif /* ALLOW_AUTODIFF_TAMC */
152     C- Check static stability with layer below and mix as needed.
153 adcroft 1.10 c CALL CONVECT(
154     c I bi,bj,iMin,iMax,jMin,jMax,K,rhoKm1,rhoK,
155     c U ConvectCount,
156     c I myTime,myIter,myThid)
157    
158     C- Pre-calculate mixing weights for interface K
159     CALL CONVECTIVE_WEIGHTS(
160     I bi,bj,K,rhoKm1,rhoK,
161     O weightA,weightB,ConvectCount,
162     I myThid)
163    
164     C- Convectively mix heat across interface K
165     CALL CONVECTIVELY_MIXTRACER(
166     I bi,bj,k,weightA,weightB,
167     U theta,
168     I myThid)
169    
170     C- Convectively mix salt across interface K
171     CALL CONVECTIVELY_MIXTRACER(
172     I bi,bj,k,weightA,weightB,
173     U salt,
174     I myThid)
175    
176     #ifdef ALLOW_PASSIVE_TRACER
177     C- Convectively mix passive tracer across interface K
178     CALL CONVECTIVELY_MIXTRACER(
179     I bi,bj,k,weightA,weightB,
180     U Tr1,
181     I myThid)
182     #endif /* ALLOW_PASSIVE_TRACER */
183 adcroft 1.11
184     #ifdef ALLOW_PTRACERS
185     C- Convectively mix passive tracers across interface K
186     IF ( usePTRACERS ) THEN
187     CALL PTRACERS_CONVECT(
188     I bi,bj,k,weightA,weightB,myThid)
189     ENDIF
190     #endif /* ALLOW_PTRACERS */
191 jmc 1.4
192 adcroft 1.2 C-- End DO K=1,Nr
193     ENDDO
194    
195 jmc 1.5 #ifdef ALLOW_TIMEAVE
196     IF (myIter.ne.nIter0 .AND. taveFreq.GT.0.) THEN
197     CALL TIMEAVE_CUMULATE(ConvectCountTave, ConvectCount, Nr,
198     I deltaTclock, bi, bj, myThid)
199     ENDIF
200     #endif /* ALLOW_TIMEAVE */
201    
202 adcroft 1.2 #ifdef ALLOW_KPP
203     C-- End IF (.NOT.useKPP)
204     ENDIF
205     #endif /* ALLOW_KPP */
206    
207     C-- End IF (DIFFERENT_MULTIPLE)
208     ENDIF
209    
210     #endif /* INCLUDE_CONVECT_CALL */
211    
212     RETURN
213     END
  ViewVC Help
Powered by ViewVC 1.1.22