Directory: | ./ |
---|---|
File: | rad/suphy0.f90 |
Date: | 2022-01-11 19:19:34 |
Exec | Total | Coverage | |
---|---|---|---|
Lines: | 276 | 293 | 94.2% |
Branches: | 9 | 44 | 20.5% |
Line | Branch | Exec | Source |
---|---|---|---|
1 | |||
2 | !OPTIONS XOPT(NOEVAL) | ||
3 | 1 | SUBROUTINE SUPHY0(KULOUT) | |
4 | |||
5 | !**** *SUPHY0* - Initialize common YOMPHY0 physics atmospheric | ||
6 | ! constants | ||
7 | |||
8 | ! Purpose. | ||
9 | ! -------- | ||
10 | ! Initialize YOMPHY0, the common that contains the parameters | ||
11 | ! for the atmospheric part of the physics of the model. | ||
12 | |||
13 | !** Interface. | ||
14 | ! ---------- | ||
15 | ! *CALL* *SUPHY0(KULOUT) | ||
16 | |||
17 | ! Explicit arguments : | ||
18 | ! -------------------- | ||
19 | ! KULOUT : Logical unit for the output | ||
20 | |||
21 | ! Implicit arguments : | ||
22 | ! -------------------- | ||
23 | ! COMMON YOMPHY0 | ||
24 | |||
25 | ! Method. | ||
26 | ! ------- | ||
27 | ! See documentation | ||
28 | |||
29 | ! Externals. | ||
30 | ! ---------- | ||
31 | |||
32 | ! Reference. | ||
33 | ! ---------- | ||
34 | ! Documentation ARPEGE | ||
35 | |||
36 | ! Author. | ||
37 | ! ------- | ||
38 | ! J.-F. Geleyn . | ||
39 | |||
40 | ! Modifications. | ||
41 | ! -------------- | ||
42 | ! Original : 90-9-1 | ||
43 | ! Ajout de GWDAMP (J.-F. Geleyn) : 91-2-2 | ||
44 | ! Modified by Michel Deque 91-04-01 (param. for convect. clouds) | ||
45 | ! Ajout de HOBST, NPCLO1/2, XNBMAX et REVGSL, remplacement de | ||
46 | ! GWDCOE par GWDSE (J.-F. Geleyn, M. Deque, L. Labbe) : 92-4-8 | ||
47 | ! Ajout de VZ0CM (E. Bazile) : 92-3-27 | ||
48 | ! Modified by C. Castejon and E. Gerard 92-02-28 (stat. clouds) | ||
49 | ! Modified by R. EL Khatib : 93-04-02 Set-up defaults controled by LECMWF | ||
50 | ! Modified by M. Deque : 95-03-20 (USURIC and USUPRC) | ||
51 | ! Modified by J.F. Geleyn C. Bossuet : 95-12-20 (GWDBC) | ||
52 | ! Modified by Luc Gerard: 97-02-26 (entrainement de la qte de mvt horizontale). | ||
53 | ! Modified by J.M. Piriou: 97-02-28 (schema de nebulosite ACNEBN). | ||
54 | ! Modified by J.M. Piriou: 97-04-17 (valeurs par defaut). | ||
55 | ! Modified by J.L. Ricard : seuil (SCO) sur les precip. conv. | ||
56 | ! modif de la turbulence residuelle | ||
57 | ! Modified by J.M. Piriou: 97-08-21 (introduce Xu-Randall cloudiness). | ||
58 | ! Modified by M. Deque : 97-10-21 (introduce GWD lift). | ||
59 | ! Modified by M. Deque : 98-02-21 (new mixing length profile). | ||
60 | ! Modified by V. Lorant : 98-08-05 (new mixing length profile). | ||
61 | ! Modified by J.M. Piriou: 98-02-11 (introduce downdrafts tuning parameters). | ||
62 | ! Modified by J.M. Piriou: 98-03-10 (introduce GRCVPP) | ||
63 | ! Modified by L. Gerard : 98-11-30 (TUDGP, TDDGP, GCOMOD) | ||
64 | ! Modified by R. EL Khatib :98-12-14 Remove LRDSPIL | ||
65 | ! Modified by J.M. Piriou: 99-01-04 (introduce GCVADS, GCVBETA) | ||
66 | ! Modified by J.M. Piriou: 99-06-18 (introduce GCVPSI, GCVALFA, USURICL AND USURICE. | ||
67 | ! Change default value for GDDEVA). | ||
68 | ! Modified by J.M. Piriou: 2000-08-23 (new use of the Richardson critical number (USURID)). | ||
69 | ! Modified by J.M. Piriou: 2000-08-23 (cloud core buoyancy as a fraction of an undilute plume (GCVNU)). | ||
70 | ! Modified by J.M. Piriou: 2000-10-06 (exponent USURIDE). | ||
71 | ! Modified by E. Bazile : 2000-12-12 CYCORA's default value. | ||
72 | ! Modified by J.M. Piriou: 2001-04-05 (introduce GCVMLT). | ||
73 | ! Modified by Y. Bouteloup:2002-06-14 (introduce RCVEVAP). | ||
74 | ! Modified by F. Bouyssel: 2002-06-25 (introduce UTILGUST, RRGAMMA, RRSCALE). | ||
75 | ! Modified by J.M. Piriou: 2002-08-19 (introduce GPBLH*). | ||
76 | ! Modified by D. Banciu: 2002-12-09 (introduce GCVPSIE). | ||
77 | ! R. El Khatib : 2001-08-07 Pruning options | ||
78 | ! J.M. Piriou : 2002-01-10 set default values to operational ones. | ||
79 | ! 03-06, move rnlcurv rnegat into yomphy0 (F. Bouyssel, C. Fischer) | ||
80 | ! E. Bazile : 2004-02-24 (introduce EDK). | ||
81 | ! M.Hamrud 01-Oct-2003 CY28 Cleaning | ||
82 | ! E. Bazile : 2004-06-30 (inroduces XKLM). | ||
83 | ! R. Brozkova : 2004-11 modifs for Xue-Randall cloud. scheme | ||
84 | ! P. Marquet & F.Bouyssel : 2004-08-18 (Lopez) | ||
85 | ! F.Bouyssel : 2005-07-18 (new defaults and new parameters for Lopez) | ||
86 | ! F.Bouyssel : 2006-01-25 (new defaults for Lopez!) | ||
87 | ! F. Vana : 2006-01-30 tunables for pTKE. | ||
88 | ! R.Brozkova : 2006-03-03 : tuning constants for mixing lengths and Charnock formulae | ||
89 | ! A.Alias : 2006-03-10 renaming KDN in RKDN | ||
90 | ! Modified by GMGEC/EAC : 2006-05 list of modif. | ||
91 | ! V. Lorant : 99-01-05 (new mixing length profile). | ||
92 | ! M. Deque : 00-03-21 (new entrainment rate VL ). | ||
93 | ! P. Marquet : 2002-11-05 TRENTRVL=0. if not ACPBLH | ||
94 | ! P. Marquet : 2004-05-27 TFVR and TFVS for ADVPRC. | ||
95 | ! P. Marquet : 2004-10-13 RAUTEFR for snow (ACMICRO). | ||
96 | ! P. Marquet : 2004-10-14 RAUTSBET for snow (ACMICRO). | ||
97 | ! A. Alias : 2005-06-23 (param. for ACCVIMPGDY) | ||
98 | ! E. Bazile & P. Marquet : 2006-04-11 AGRE1,AGRERICR,AJBUMIN,RCOFLM | ||
99 | ! for LPBLE. | ||
100 | ! F.Bouyssel : 2006-10-30 RQLCV,RQICVMAX,RQICVMIN,RHEVAP | ||
101 | ! M. Bellus : 03-Oct-2006 : ALARO-0 phasing (defaults for prognostic | ||
102 | ! convection, pTKE and PIL microphysics) | ||
103 | ! A.Alias : 2006-07 FEVAPC added for ACCVIMPGY (JF Gueremy) | ||
104 | ! ------------------------------------------------------------------------ | ||
105 | |||
106 | USE PARKIND1 ,ONLY : JPIM ,JPRB | ||
107 | USE YOMHOOK ,ONLY : LHOOK, DR_HOOK | ||
108 | |||
109 | ! Ce qui concerne NULNAM commente par MPL le 15.04.09 | ||
110 | !USE YOMLUN , ONLY : NULNAM | ||
111 | USE YOMCT0B , ONLY : LECMWF | ||
112 | USE YOMSIMPHL, ONLY : LSIMPH ,LSMOOTHD ,LSMOOTHA | ||
113 | USE YOMPHY , ONLY : LCVCAS ,LCVLIS ,LCVRA ,LMPHYS ,LCONDWT ,& | ||
114 | & LPROCLD ,LNEBCO ,LNEBR ,LECT ,LNEBGR ,LCVRAV3 ,LPIL | ||
115 | USE YOMPHY0 , ONLY : TYM ,NPCLO1 ,NPCLO2 ,AERCS1 ,AERCS3 ,& | ||
116 | & AERCS5 ,ALMAV ,BEDIFV ,ECMNP ,EDB ,EDC ,& | ||
117 | & EDD ,EVAP ,FONT ,GALP ,GCCSV ,GCOMOD ,& | ||
118 | & GCVADS ,GCVALFA ,GCVBETA ,GCVPSI ,GDDEVA ,GDDSDE ,& | ||
119 | & GRCVPP ,GWDAMP ,GWDBC ,GWDCCO ,GWDCD ,GWDLT ,& | ||
120 | & AHCLPV ,GWDSE ,HOBST ,HUCOE ,HUCOE2 ,HUTIL ,& | ||
121 | & HUTIL1 ,HUTIL2 ,& | ||
122 | & QSSC ,QSMIN ,QSNEBC ,QSNEBS ,QSSUSC ,QSSUSS ,& | ||
123 | & QXRAL ,QXRDEL ,QXRR ,REVGSL ,RTCAPE ,SCO ,& | ||
124 | & SENSL ,SNNBCO ,SPNBCO ,SXNBCO ,TCA ,TCT ,& | ||
125 | & TCW ,TENTR ,TENTRX ,TDDGP ,TUDGP ,TURB ,& | ||
126 | & TVF ,UHDIFV ,USDMLT ,USUPRC ,USURIC ,USURID ,& | ||
127 | & USURIDE ,USURICE ,USURICL ,VCHRNK ,VKARMN ,VZ0CM ,& | ||
128 | & XNBMAX ,RICRLM ,XBLM ,XMINLM ,XMAXLM ,XWSALM ,& | ||
129 | & XWSBLM ,GCVNU ,GCVMLT ,RCVEVAP ,GPBLHK0 ,GPBLHRA ,& | ||
130 | & UTILGUST ,RRGAMMA ,RRSCALE ,GCVPSIE ,QSUSXC ,QSUSXS ,& | ||
131 | & ETACUT ,RNEGAT ,RNLCURV ,QSSUSV ,QXRHX ,GCISMIN ,& | ||
132 | & GWDPROF ,GWDVALI ,RCIN ,EDK ,XKLM ,RPHI0 ,& | ||
133 | & RPHIR ,QXRTGH ,ADISE ,ADISI ,AECLS3 ,AECLS4 ,& | ||
134 | & AKN ,ALD ,ALPHAE ,ALPHAT ,ECTMIN ,UCWSTAR ,& | ||
135 | & UDECT ,USHEARM ,UPRETMIN ,UPRETMAX ,ARSCH ,ARSCQ ,& | ||
136 | & ARSC1 ,ARSB2 ,ACBRPHIM ,ALMAVE ,RICRET ,STTBMIN ,& | ||
137 | & AGREKE ,RDTFAC ,RAUTEFS ,RNINTR ,RNINTS ,RQLCR ,& | ||
138 | & RQICRMAX ,RQICRMIN ,RACCEF ,RRIMEF ,RHCRIT1 ,RHCRIT2 ,& | ||
139 | & RETAMIN ,TFVR ,TFVS ,RAUTEFR ,RAUTSBET ,GRHCMOD ,& | ||
140 | & RQICRT1 ,RQICRT2 ,RQICRSN ,RQCRNS ,RFACNSM ,RAGGEF ,& | ||
141 | & RQLCV ,RQICVMAX ,RQICVMIN ,RHEVAP ,& | ||
142 | & AGRE1 ,AGRERICR ,AJBUMIN ,RCOFLM ,& | ||
143 | & A0ML_AU ,A0ML_AT ,A0ML_BU ,A0ML_BT ,VZIUSTAR0,& | ||
144 | & TENTRVL ,TRENTRV ,UETEPS ,UPRECLP ,& | ||
145 | & ARSC2 ,ARSCT ,AGRE2 ,AGREF ,& | ||
146 | & AJ1PEPS ,AJ1MEPS ,NAJITER ,& | ||
147 | & ALFX ,TCTC ,TVFC ,GAMAP1 ,RKDN ,& | ||
148 | & VVN ,VVX ,FENTRT ,HCMIN ,FQLIC ,FNEBC ,& | ||
149 | & NUPTKE ,GAMTKE ,RCOLL ,RFALLL ,TDDBU ,& | ||
150 | & TDDFR ,TUDBU ,TUDFR ,GCVACHI ,GCVALMX ,GCVADMW ,& | ||
151 | & GCVBEE ,GCVEEX ,ECMNPI ,GFRIC ,GCVSQDN ,GCVSQDR ,& | ||
152 | & GCVSQDCX ,GRRINTE ,GRRMINA ,GDDBETA ,GDDEVF ,GDDWPF ,& | ||
153 | & TENTRD ,RDPHIC ,GWBFAUT ,RWBF1 ,RWBF2 ,RAUITN ,& | ||
154 | & RAUITX ,RAUIUSTE ,RSMDNEBX ,RSMDTX ,NSMTPA ,NSMTPB ,& | ||
155 | & FEVAPC | ||
156 | |||
157 | |||
158 | IMPLICIT NONE | ||
159 | |||
160 | INTEGER(KIND=JPIM),INTENT(IN) :: KULOUT | ||
161 | REAL(KIND=JPRB) :: ZHOOK_HANDLE | ||
162 | |||
163 | INTERFACE | ||
164 | SUBROUTINE ABOR1(CDTEXT) | ||
165 | CHARACTER(LEN=*) :: CDTEXT | ||
166 | END SUBROUTINE ABOR1 | ||
167 | END INTERFACE | ||
168 | INTERFACE | ||
169 | SUBROUTINE POSNAM(KULNAM,CDNAML) | ||
170 | USE PARKIND1 ,ONLY : JPIM ,JPRB | ||
171 | INTEGER(KIND=JPIM),INTENT(IN) :: KULNAM | ||
172 | CHARACTER(LEN=*) ,INTENT(IN) :: CDNAML | ||
173 | END SUBROUTINE POSNAM | ||
174 | END INTERFACE | ||
175 | INTERFACE | ||
176 | SUBROUTINE SURHCRI(KULOUT) | ||
177 | USE PARKIND1 ,ONLY : JPIM ,JPRB | ||
178 | INTEGER (KIND=JPIM),INTENT(IN) :: KULOUT | ||
179 | END SUBROUTINE SURHCRI | ||
180 | END INTERFACE | ||
181 | |||
182 | ! ------------------------------------------------------------------ | ||
183 | NAMELIST/NAMPHY0/& | ||
184 | &AERCS1,AERCS3,AERCS5,ALMAV,BEDIFV,ECMNP,ETACUT,AHCLPV & | ||
185 | &,EDB,EDC,EDD,EVAP,FONT,GALP,GCISMIN,GCCSV,GCVADS,GCVALFA,GCVBETA,GCVMLT & | ||
186 | &,GCVNU,GCVPSI,GCVPSIE,GCOMOD & | ||
187 | &,GDDEVA,GDDSDE,GPBLHK0,GPBLHRA,GRCVPP,GWDAMP,GWDBC,GWDCCO,GWDCD & | ||
188 | &,GWDLT,GWDSE,GWDPROF,GWDVALI,HOBST,HUCOE,HUTIL,HUTIL1,HUTIL2,NPCLO1 & | ||
189 | &,NPCLO2,QSSC,QSMIN,QSNEBC,QSNEBS,QSSUSC,QSSUSS,QSSUSV,QSUSXC,QSUSXS & | ||
190 | &,QXRAL,QXRDEL,QXRHX,QXRR,QXRTGH,RPHI0,RPHIR & | ||
191 | &,RCIN,RCVEVAP,RDTFAC,REVGSL,RTCAPE,RICRLM,RNEGAT,RNLCURV,RRGAMMA,RRSCALE & | ||
192 | &,SCO,SENSL & | ||
193 | &,SNNBCO,SPNBCO & | ||
194 | &,SXNBCO,TCA,TCT,TCW,TENTR,TENTRX,TDDGP,TUDGP,TURB,TVF & | ||
195 | &,TYM,UHDIFV,USDMLT,USUPRC,USURIC,USURID,USURIDE,USURICE,USURICL,UTILGUST & | ||
196 | &,VCHRNK,VKARMN,VZ0CM,VZIUSTAR0,XNBMAX,XBLM,XMAXLM,XMINLM,XNBMAX,EDK,XKLM & | ||
197 | &,ADISE,ADISI,AECLS3,AECLS4,AKN,ALD,ALPHAE,ALPHAT,ECTMIN & | ||
198 | &,UCWSTAR,UDECT,USHEARM,UPRETMIN,UPRETMAX & | ||
199 | &,ARSCH,ARSCQ,ARSC1,ARSB2,ACBRPHIM,ALMAVE,RICRET,STTBMIN,AGREKE & | ||
200 | &,RAUTEFS,RNINTR,RNINTS,RQLCR,RQICRMAX,RQICRMIN,RACCEF,RRIMEF,RAGGEF & | ||
201 | &,RHCRIT1,RHCRIT2,RETAMIN,RFACNSM,TFVR,TFVS,RAUTEFR,RAUTSBET,GRHCMOD & | ||
202 | &,RQICRT1,RQICRT2,RQICRSN,RQCRNS,RQICVMAX,RQICVMIN,RQLCV,RHEVAP & | ||
203 | &,AGRE1,AGRERICR,AJBUMIN,RCOFLM,A0ML_AU,A0ML_AT,A0ML_BU,A0ML_BT & | ||
204 | &,NUPTKE,GAMTKE,RCOLL,RFALLL & | ||
205 | &,TDDBU,TDDFR,TUDBU,TUDFR,GCVACHI,GCVALMX,GCVADMW,GCVBEE,GCVEEX,ECMNPI & | ||
206 | &,GFRIC,GCVSQDN,GCVSQDR,GCVSQDCX,GRRINTE,GRRMINA,GDDBETA,GDDEVF,GDDWPF & | ||
207 | &,TENTRD & | ||
208 | ! For old ACCVIMP | ||
209 | &,TENTRVL,TRENTRV & | ||
210 | ! For TKE scheme | ||
211 | &,UETEPS,UPRECLP & | ||
212 | &,ARSC2,ARSCT & | ||
213 | ! For Grenier (2000) top-PBL entrainment scheme | ||
214 | &,AGRE2,AGREF & | ||
215 | ! For dry conv. adjustment scheme | ||
216 | &,AJ1MEPS,AJ1PEPS,NAJITER & | ||
217 | ! For ACCVIMPGY | ||
218 | &,ALFX,TCTC,TVFC,GAMAP1,RKDN,VVN,VVX,FENTRT,HCMIN,FQLIC,FNEBC,FEVAPC & | ||
219 | ! | ||
220 | &,RDPHIC,GWBFAUT,RWBF1,RWBF2,RAUITN,RAUITX,RAUIUSTE & | ||
221 | &,RSMDNEBX,RSMDTX,NSMTPA,NSMTPB | ||
222 | ! ------------------------------------------------------------------ | ||
223 | |||
224 | ! ------------------------------------------------------------------ | ||
225 | |||
226 | !* 1. Set default values. | ||
227 | ! ------------------- | ||
228 | |||
229 | ! 1.1 Set implicit default values | ||
230 | |||
231 |
1/2✓ Branch 0 taken 1 times.
✗ Branch 1 not taken.
|
1 | IF (LHOOK) CALL DR_HOOK('SUPHY0',0,ZHOOK_HANDLE) |
232 | 1 | A0ML_AU=4.5_JPRB | |
233 | 1 | A0ML_AT=5.0_JPRB | |
234 | 1 | A0ML_BU=3.0_JPRB | |
235 | 1 | A0ML_BT=0.8_JPRB | |
236 | 1 | AERCS1=0.2726_JPRB | |
237 | 1 | AERCS3=-0.4239_JPRB | |
238 | 1 | AERCS5=0.3595_JPRB | |
239 | 1 | ALMAV=300._JPRB | |
240 | 1 | BEDIFV=0.05_JPRB | |
241 | 1 | ECMNP=3000._JPRB | |
242 | 1 | EDB=5._JPRB | |
243 | 1 | EDC=5._JPRB | |
244 | 1 | EDD=5._JPRB | |
245 | 1 | EDK=1.0_JPRB | |
246 | 1 | ETACUT=1.0_JPRB | |
247 | 1 | EVAP=4.8E+06_JPRB | |
248 | 1 | FONT=2.4E+04_JPRB | |
249 | 1 | GALP=0.53_JPRB | |
250 | 1 | GCISMIN=6.7E-05_JPRB | |
251 | 1 | GCCSV=0._JPRB | |
252 | 1 | GCOMOD=1._JPRB | |
253 | 1 | GCVADS=0.8_JPRB | |
254 | 1 | GCVALFA=4.5E-05_JPRB | |
255 | 1 | GCVBETA=0.2_JPRB | |
256 | 1 | GCVMLT=0.00016_JPRB | |
257 | 1 | GCVNU=2.5E-05_JPRB | |
258 | 1 | GCVPSI=0.5_JPRB | |
259 | 1 | GCVPSIE=0.0_JPRB | |
260 | 1 | GDDEVA=0.25_JPRB | |
261 | 1 | GDDSDE=0.5_JPRB | |
262 | 1 | GPBLHK0=0.25_JPRB | |
263 | 1 | GPBLHRA=4._JPRB | |
264 | 1 | GRCVPP=1._JPRB | |
265 | 1 | GWDAMP=0.6_JPRB | |
266 | 1 | GWDBC=2._JPRB | |
267 | 1 | GWDCCO=1._JPRB | |
268 | 1 | GWDCD=6._JPRB | |
269 | 1 | GWDLT=0._JPRB | |
270 | 1 | GWDSE=3.5E-03_JPRB | |
271 | 1 | GWDPROF=1._JPRB | |
272 | 1 | GWDVALI=0._JPRB | |
273 | 1 | HOBST=3._JPRB | |
274 | 1 | HUCOE=2._JPRB | |
275 | 1 | HUCOE2=0.4_JPRB | |
276 | 1 | HUTIL=1.8_JPRB | |
277 | 1 | HUTIL1=-0.6_JPRB | |
278 | 1 | HUTIL2=1.1_JPRB | |
279 | 1 | NPCLO1=1 | |
280 | 1 | NPCLO2=1 | |
281 | 1 | QSSC=1600._JPRB | |
282 | 1 | QSMIN=1.E-4_JPRB | |
283 | 1 | QSNEBC=26000._JPRB | |
284 | 1 | QSNEBS=0.7_JPRB | |
285 | 1 | QSSUSC=1._JPRB | |
286 | 1 | QSSUSS=0.25_JPRB | |
287 | 1 | QSSUSV=0._JPRB | |
288 | 1 | QSUSXC=3.3E-05_JPRB | |
289 | 1 | QSUSXS=3.3E-05_JPRB | |
290 | 1 | RPHI0=0._JPRB | |
291 | 1 | RPHIR=1750._JPRB | |
292 | 1 | QXRAL=10000._JPRB | |
293 | 1 | QXRDEL=0._JPRB | |
294 | 1 | QXRHX=1._JPRB | |
295 | 1 | QXRR=0.5_JPRB | |
296 | 1 | QXRTGH=3.5_JPRB | |
297 | 1 | RCIN=0._JPRB | |
298 | 1 | RCVEVAP=0._JPRB | |
299 | 1 | REVGSL=80._JPRB | |
300 | 1 | RTCAPE=10800._JPRB | |
301 | 1 | SCO=-20._JPRB | |
302 | 1 | SENSL=1._JPRB | |
303 | 1 | SNNBCO=0._JPRB | |
304 | 1 | SPNBCO=3000._JPRB | |
305 | 1 | SXNBCO=0.5_JPRB | |
306 | 1 | TCA=1._JPRB | |
307 | 1 | TCT=1.E-4_JPRB | |
308 | 1 | TCW=8.E-4_JPRB | |
309 | 1 | TENTR=2.5E-06_JPRB | |
310 | 1 | TENTRX=8.E-05_JPRB | |
311 | 1 | TDDGP=0.8_JPRB | |
312 | 1 | TUDGP=0.8_JPRB | |
313 | 1 | TURB=1._JPRB | |
314 | 1 | TVF=1._JPRB | |
315 | 1 | TYM(1)=0.92_JPRB | |
316 | 1 | TYM(2)=0.74_JPRB | |
317 | 1 | TYM(3)=16.6_JPRB | |
318 | 1 | TYM(4)=10.1_JPRB | |
319 | 1 | TYM(5)=0.08_JPRB | |
320 | 1 | UHDIFV=8.E-04_JPRB | |
321 | 1 | USDMLT=1.25E+04_JPRB | |
322 | 1 | USUPRC=0.0_JPRB | |
323 | 1 | USURIC=1.0_JPRB | |
324 | 1 | USURID=0.035_JPRB | |
325 | 1 | USURIDE=1.0_JPRB | |
326 | 1 | USURICE=0.5_JPRB | |
327 | 1 | USURICL=4._JPRB | |
328 | 1 | UTILGUST=0.125_JPRB | |
329 | 1 | VCHRNK=0.021_JPRB | |
330 | 1 | VKARMN=0.4_JPRB | |
331 | 1 | VZ0CM=1.5E-04_JPRB | |
332 | 1 | VZIUSTAR0=0._JPRB | |
333 | 1 | XNBMAX=1._JPRB | |
334 | 1 | AHCLPV=1000._JPRB | |
335 | 1 | RICRLM=0.5_JPRB | |
336 | 1 | RRGAMMA=0.8_JPRB | |
337 | 1 | RRSCALE=1.15E-4_JPRB | |
338 | 1 | XBLM=6.5_JPRB | |
339 | 1 | XKLM=1.0_JPRB | |
340 | 1 | XMAXLM=3000._JPRB | |
341 | 1 | XMINLM=500._JPRB | |
342 | 1 | XWSALM=0.1_JPRB | |
343 | 1 | XWSBLM=7.0_JPRB | |
344 | 1 | RNEGAT = -7.E-05_JPRB | |
345 | 1 | RNLCURV = 7.E+04_JPRB | |
346 | |||
347 | ! - - - - - - - - - - - - - - | ||
348 | ! The old Convective scheme : | ||
349 | ! - - - - - - - - - - - - - - | ||
350 | 1 | TRENTRV=1._JPRB | |
351 | 1 | TENTRVL=-1.0_JPRB | |
352 | |||
353 |
1/4✗ Branch 0 not taken.
✓ Branch 1 taken 1 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
|
1 | IF ( LCONDWT.AND.LPROCLD ) THEN |
354 | ✗ | RDTFAC=0.5_JPRB | |
355 | ELSE | ||
356 | 1 | RDTFAC=1.0_JPRB | |
357 | ENDIF | ||
358 | |||
359 | ! - - - - - - - - - - | ||
360 | ! Lopez Microphysics : | ||
361 | ! - - - - - - - - - - | ||
362 | 1 | RAUTEFR=1.E-03_JPRB | |
363 | 1 | RAUTEFS=1.E-03_JPRB | |
364 | 1 | RAUTSBET=0.025_JPRB | |
365 | 1 | RNINTR=8.E+06_JPRB | |
366 | 1 | RNINTS=2.E+06_JPRB | |
367 | 1 | RQLCR=2.E-04_JPRB | |
368 | 1 | RQICRMAX=0.3E-04_JPRB | |
369 | 1 | RQICRMIN=0.2E-06_JPRB | |
370 | 1 | RQLCV=2.E-04_JPRB | |
371 | 1 | RQICVMAX=0.3E-04_JPRB | |
372 | 1 | RQICVMIN=0.2E-06_JPRB | |
373 | 1 | RQICRT1=-80._JPRB | |
374 | 1 | RQICRT2=30._JPRB | |
375 | 1 | RQICRSN=0.5_JPRB | |
376 | 1 | RQCRNS=0.03_JPRB | |
377 | 1 | RACCEF=1._JPRB | |
378 | 1 | RAGGEF=0.2_JPRB | |
379 | 1 | RRIMEF=1._JPRB | |
380 | 1 | RHEVAP=0.0_JPRB | |
381 | 1 | RHCRIT1=0.5_JPRB | |
382 | 1 | RHCRIT2=0.91_JPRB | |
383 | 1 | RETAMIN=0.4_JPRB | |
384 | 1 | RFACNSM=1.4_JPRB | |
385 | 1 | TFVR=5.0_JPRB | |
386 | 1 | TFVS=0.6_JPRB | |
387 | 1 | GRHCMOD=0.3_JPRB | |
388 | |||
389 | |||
390 | !-------------------------- | ||
391 | ! PIL MICROPHYSICS | ||
392 | !-------------------------- | ||
393 | 1 | RDPHIC=10000._JPRB | |
394 | 1 | GWBFAUT=15._JPRB | |
395 | ! two constants for ACPLUIE_PROG: | ||
396 | 1 | RWBF1=300._JPRB | |
397 | 1 | RWBF2=4._JPRB | |
398 | ! | ||
399 | 1 | RAUITN=233.15_JPRB | |
400 | 1 | RAUITX=263.15_JPRB | |
401 | 1 | RAUIUSTE=0.025_JPRB | |
402 | 1 | RSMDNEBX=0.2_JPRB | |
403 | 1 | RSMDTX=1.0_JPRB | |
404 | 1 | NSMTPA=2 | |
405 | 1 | NSMTPB=3 | |
406 | 1 | RCOLL=6.9E-03_JPRB | |
407 | 1 | RFALLL=1.0_JPRB | |
408 | !-------------------------- | ||
409 | ! Prognostic convection physical parameters | ||
410 | !-------------------------- | ||
411 | 1 | TUDBU=0.5_JPRB | |
412 | 1 | TDDBU=0.5_JPRB | |
413 | 1 | TUDFR=0.0012_JPRB | |
414 | 1 | TDDFR=0.0006_JPRB | |
415 | 1 | GCVALMX=0.95_JPRB | |
416 | ! Pas d'activite historique: defaut tres haut, | ||
417 | ! pour que KUO joue tout seul: | ||
418 | 1 | GCVACHI=1.E9_JPRB | |
419 | 1 | GCVADMW =0 | |
420 | ! Explicit entrainment: | ||
421 | 1 | GCVBEE=0._JPRB ! 0.2_JPRB | |
422 | 1 | GCVEEX=1._JPRB | |
423 | 1 | ECMNPI=3000._JPRB | |
424 | 1 | GFRIC=-1.0_JPRB !1.E-3_JPRB | |
425 | ! Sqeezing: | ||
426 | 1 | GCVSQDN=0.01_JPRB | |
427 | 1 | GCVSQDR=0.8_JPRB | |
428 | 1 | GCVSQDCX=1.0_JPRB ! No squeezing | |
429 | ! PRECIPITATING AREA PARAMETERS (aplpar): | ||
430 | 1 | GRRINTE=2._JPRB | |
431 | 1 | GRRMINA=1.E-5_JPRB | |
432 | ! DD explicit detrainment: | ||
433 | 1 | GDDBETA=0.2_JPRB | |
434 | ! ACMODO DD PARAMETERS: | ||
435 | 1 | GDDEVF=0.5_JPRB | |
436 | 1 | GDDWPF=0._JPRB | |
437 | ! DD ENTRAINMENT RATE: | ||
438 | 1 | TENTRD=1.E-4_JPRB | |
439 | |||
440 | ! - - - - - - - - | ||
441 | ! For TKE scheme : | ||
442 | ! - - - - - - - - | ||
443 | 1 | ADISE=-0.5_JPRB | |
444 | 1 | ADISI=1.5_JPRB | |
445 | 1 | AECLS3=3.75_JPRB | |
446 | 1 | AECLS4=0.3_JPRB | |
447 | 1 | AKN=0.2_JPRB | |
448 | 1 | ALD=1.4_JPRB | |
449 | 1 | ALPHAE=1.0_JPRB | |
450 | 1 | ALPHAT=1.0_JPRB | |
451 | 1 | ECTMIN=1.E-10_JPRB | |
452 | |||
453 | 1 | UCWSTAR=1.0_JPRB/3._JPRB | |
454 | 1 | UDECT=5._JPRB | |
455 | 1 | USHEARM=1.E-04_JPRB | |
456 | 1 | UPRETMIN=60000._JPRB | |
457 | 1 | UPRETMAX=97500._JPRB | |
458 | 1 | ARSCH=4._JPRB | |
459 | 1 | ARSCQ=1.2_JPRB | |
460 | 1 | ARSC1=2.0_JPRB/(3._JPRB*ARSCH*ARSCQ) | |
461 | 1 | ARSB2=3._JPRB*ARSCH*ARSC1/2.0_JPRB | |
462 | 1 | ACBRPHIM=2.2_JPRB | |
463 | 1 | ALMAVE=0._JPRB | |
464 | 1 | RICRET=0.195_JPRB | |
465 | |||
466 | 1 | STTBMIN=SQRT(3._JPRB) | |
467 | ! A minimum value for the adimentional jump in boyancy : d(Theta)/Theta | ||
468 | 1 | AJBUMIN=0.005_JPRB | |
469 | ! The "Master Length" is equal to "RCOFLM*Z_PBL" | ||
470 | 1 | RCOFLM=0.085_JPRB | |
471 | |||
472 | |||
473 | ! TKE (P.Marquet) | ||
474 | 1 | UETEPS=1.0_JPRB | |
475 | |||
476 | ! - - - - - - - - - - - - - - - - | ||
477 | ! For dry conv. adjustment scheme : | ||
478 | ! - - - - - - - - - - - - - - - - | ||
479 | 1 | AJ1MEPS=0.99_JPRB | |
480 | 1 | AJ1PEPS=10.0_JPRB | |
481 | 1 | NAJITER=30 | |
482 | |||
483 | ! - - - - - - - - - - - - - - - - - - - - - - - | ||
484 | ! For Grenier (2000) top-PBL entrainment scheme : | ||
485 | ! - - - - - - - - - - - - - - - - - - - - - - - | ||
486 | 1 | AGRE1=0.16_JPRB | |
487 | 1 | AGRE2=15._JPRB | |
488 | 1 | AGREF=0.8_JPRB | |
489 | 1 | AGREKE=5.0_JPRB | |
490 | 1 | AGRERICR=50._JPRB | |
491 | |||
492 | ! Pseudo prognostic TKE scheme | ||
493 | 1 | NUPTKE=0.52_JPRB | |
494 | 1 | GAMTKE=0.5_JPRB | |
495 | |||
496 | |||
497 | |||
498 | ! - - - - - - - | ||
499 | ! For ACCVIMPGY : | ||
500 | ! - - - - - - - | ||
501 | 1 | ALFX=10.E-02_JPRB | |
502 | 1 | TCTC=1.60E-04_JPRB | |
503 | 1 | TVFC=1._JPRB | |
504 | 1 | GAMAP1=1.5_JPRB | |
505 | 1 | RKDN=30.E-06_JPRB | |
506 | 1 | VVN=0.0_JPRB | |
507 | 1 | VVX=-45._JPRB | |
508 | 1 | FENTRT=2.5_JPRB | |
509 | 1 | HCMIN=0.0_JPRB | |
510 | 1 | FQLIC=2.5_JPRB | |
511 | 1 | FNEBC=25.0_JPRB | |
512 | 1 | FEVAPC=3.5_JPRB | |
513 | |||
514 | ! 1.2 Modify default values according to LECMWF | ||
515 | |||
516 | IF (LECMWF) THEN | ||
517 | ELSE | ||
518 | ENDIF | ||
519 | |||
520 | ! ------------------------------------------------------------------ | ||
521 | |||
522 | !* 2. Modify default values. | ||
523 | ! ---------------------- | ||
524 | |||
525 | ! Ce qui concerne NAMPHY0 commente par MPL le 15.04.09 | ||
526 | !CALL POSNAM(NULNAM,'NAMPHY0') | ||
527 | !READ(NULNAM,NAMPHY0) | ||
528 | |||
529 |
2/10✗ Branch 0 not taken.
✓ Branch 1 taken 1 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✓ Branch 9 taken 1 times.
|
1 | IF(.NOT.(LNEBCO.AND.(LNEBR.OR.LNEBGR.OR.LECT)) .AND.LCVRAV3) THEN |
530 | ✗ | WRITE(UNIT=KULOUT,FMT='(A)') ' ' | |
531 | ✗ | WRITE(UNIT=KULOUT,FMT='(A)') ' !' | |
532 | ✗ | WRITE(UNIT=KULOUT,FMT='(A)') ' ! TENTRVL is set to 0. in SUPHY0 !!' | |
533 | ✗ | WRITE(UNIT=KULOUT,FMT='(A)') ' !' | |
534 | ✗ | WRITE(UNIT=KULOUT,FMT='(A)') ' ' | |
535 | ✗ | TENTRVL = 0.0_JPRB | |
536 | ENDIF | ||
537 | ! ------------------------------------------------------------------ | ||
538 | |||
539 | !* 3. Print final values. | ||
540 | ! ------------------- | ||
541 | |||
542 | 1 | WRITE(UNIT=KULOUT,FMT='('' COMMON YOMPHY0 '')') | |
543 | WRITE(UNIT=KULOUT,FMT='('' AERCS1 = '',E11.4,'' AERCS3 = '',E11.4 & | ||
544 | & ,'' AERCS5 = '',E11.4,'' ALMAV = '',E11.4,'' ECMNP = '',E11.4 & | ||
545 | & ,'' EDB = '',E11.4,/,'' EDC = '',E11.4,'' EDD = '',E11.4 & | ||
546 | & ,'' EDK = '',E11.4,'' ETACUT = '',E11.4 & | ||
547 | & ,'' EVAP = '',E11.4,'' FONT = '',E11.4,'' GWDAMP = '',E11.4 & | ||
548 | & ,'' GWDSE = '',E11.4,'' GWDBC = '',E11.4,/& | ||
549 | & ,'' GWDCD = '',E11.4,/& | ||
550 | & ,'' GWDPROF = '',E11.4,'' GWDVALI = '',E11.4,/& | ||
551 | & ,'' HUCOE = '',E11.4,'' HUTIL = '',E11.4 ,'' HUTIL1 = '',E11.4 ,'' HUTIL2 = '',E11.4 & | ||
552 | & ,'' VCHRNK = '',E11.4,'' VKARMN = '',E11.4 & | ||
553 | & ,'' SNNBCO = '',E11.4,/,'' SPNBCO = '',E11.4 & | ||
554 | & ,'' SXNBCO = '',E11.4,'' HOBST = '',E11.4,'' NPCLO1 = '',I4 & | ||
555 | & ,'' NPCLO2 = '',I4,'' RCIN = '',E11.4,'' RCVEVAP = '',E11.4 & | ||
556 | & ,'' REVGSL = '',E11.4,/& | ||
557 | & ,'' RTCAPE= '',E11.4,'' GCOMOD= '',E11.4,/& | ||
558 | & ,'' UHDIFV = '',E11.4,'' VZ0CM = '',E11.4,'' VZIUSTAR0 = '',E11.4 & | ||
559 | & ,'' XNBMAX = '',E11.4,'' GALP = '',E11.4,'' SENSL = '',E11.4 & | ||
560 | & ,'' TCA = '',E11.4,/,'' TCT = '',E11.4,'' TCW = '',E11.4 & | ||
561 | & ,'' TURB = '',E11.4,'' TVF = '',E11.4,/& | ||
562 | & ,'' TYM = '',5E11.4,'' QSMIN = '',E11.4,/& | ||
563 | & ,'' QSSC = '',E11.4,'' RPHI0 = '',E11.4,'' RPHIR = '',E11.4 & | ||
564 | & ,'' BEDIFV = '',E11.4,'' SCO = '',E11.4 & | ||
565 | & ,'' USDMLT = '',E11.4,/& | ||
566 | & ,'' GDDEVA = '',E11.4,'' GDDSDE = '',E11.4,/& | ||
567 | & ,'' GWDCCO = '',E11.4 & | ||
568 | & ,'' HUCOE2 = '',E11.4,/& | ||
569 | & ,'' TENTR = '',E11.4,'' TENTRX = '',E11.4 & | ||
570 | & ,'' TUDGP = '',E11.4,'' TDDGP = '',E11.4 & | ||
571 | & ,'' GRCVPP = '',E11.4 & | ||
572 | & )')& | ||
573 | 1 | & AERCS1,AERCS3,AERCS5,ALMAV,ECMNP,EDB,EDC,EDD,EDK,ETACUT,EVAP,FONT,& | |
574 | 1 | & GWDAMP,GWDSE,GWDBC,GWDCD,GWDPROF,GWDVALI,& | |
575 | 1 | & HUCOE,HUTIL,HUTIL1,HUTIL2,VCHRNK,VKARMN,SNNBCO,SPNBCO,SXNBCO,HOBST,& | |
576 | 1 | & NPCLO1,NPCLO2,RCIN,RCVEVAP,REVGSL,RTCAPE,GCOMOD,& | |
577 | 1 | & UHDIFV,VZ0CM,VZIUSTAR0,XNBMAX,GALP,SENSL,TCA,TCT,TCW,& | |
578 | 1 | & TURB,TVF,TYM,QSMIN,QSSC,RPHI0,RPHIR,BEDIFV,SCO,USDMLT,& | |
579 | 1 | & GDDEVA,GDDSDE,GWDCCO,& | |
580 | 1 | & HUCOE2,& | |
581 | 2 | & TENTR,TENTRX,TUDGP,TDDGP,GRCVPP | |
582 | WRITE(UNIT=KULOUT,FMT='('' USUPRC = '',E11.4,'' USURIC = '',E11.4 & | ||
583 | & ,'' QSNEBC = '',E11.4,'' QSNEBS = '',E11.4 & | ||
584 | & ,'' QSSUSC = '',E11.4,'' QSSUSS = '',E11.4,'' QSSUSV = '',E11.4 & | ||
585 | & ,'' QSUSXC = '',E11.4,'' QSUSXS = '',E11.4 & | ||
586 | & ,'' GCCSV = '',E11.4,/& | ||
587 | & ,'' QXRAL = '',E11.4,'' QXRDEL = '',E11.4 & | ||
588 | & ,'' QXRHX = '',E11.4,'' QXRR = '',E11.4,'' QXRTGH = '',E11.4,/& | ||
589 | & ,'' GWDLT = '',E11.4,'' AHCLPV = '',E11.4,/& | ||
590 | & ,'' GCVADS = '',E11.4,'' GCVBETA= '',E11.4 & | ||
591 | & ,'' RICRLM = '',E11.4,'' XBLM = '',E11.4,'' XKLM = '',E11.4 & | ||
592 | & ,'' XMINLM = '',E11.4,'' XMAXLM = '',E11.4 & | ||
593 | & )')& | ||
594 | 1 | & USUPRC,USURIC,QSNEBC,QSNEBS,QSSUSC,QSSUSS,QSSUSV,QSUSXC,QSUSXS,GCCSV,QXRAL & | |
595 | 1 | & ,QXRDEL,QXRHX,QXRR,QXRTGH,GWDLT,AHCLPV,GCVADS,GCVBETA,RICRLM,XBLM,XKLM,XMINLM & | |
596 | 2 | & ,XMAXLM | |
597 | WRITE(UNIT=KULOUT,FMT='('' XWSALM = '',E11.4,'' XWSBLM = '',E11.4 & | ||
598 | & ,'' GCVALFA= '',E11.4,'' GCVPSI = '',E11.4,'' GCVPSIE = '',E11.4 & | ||
599 | & ,'' USURICL= '',E11.4,'' USURICE= '',E11.4,/ & | ||
600 | & ,'' USURID= '',E11.4,'' USURIDE= '',E11.4,'' GCVNU= '',E11.4 & | ||
601 | & ,'' GCVMLT= '',E11.4,'' GPBLHK0= '',E11.4,'' GPBLHRA = '',E11.4 & | ||
602 | & ,'' UTILGUST= '',E11.4,'' RRGAMMA= '',E11.4,'' RRSCALE = '',E11.4 & | ||
603 | & )')& | ||
604 | 1 | & XWSALM,XWSBLM,GCVALFA,GCVPSI,GCVPSIE,USURICL,USURICE,USURID,USURIDE,GCVNU,GCVMLT,GPBLHK0 & | |
605 | 2 | & ,GPBLHRA,UTILGUST,RRGAMMA,RRSCALE | |
606 | WRITE(UNIT=KULOUT,FMT='('' RNEGAT == '',E11.4,'' RNLCURV = '',E11.4 & | ||
607 | 1 | & )') RNEGAT,RNLCURV | |
608 | 1 | WRITE(UNIT=KULOUT,FMT='('' RDTFAC == '',E11.4 )') RDTFAC | |
609 | 1 | WRITE(UNIT=KULOUT,FMT='(9(5(A,E11.4),/))') 'GCISMIN=',GCISMIN | |
610 | WRITE(UNIT=KULOUT,FMT='('' A0ML_AU == '',E11.4,'' A0ML_AT == '',E11.4 & | ||
611 | 1 | &,'' A0ML_BU == '',E11.4,'' A0ML_BT == '',E11.4 )') A0ML_AU,A0ML_AT& | |
612 | 2 | &,A0ML_BU,A0ML_BT | |
613 | |||
614 | ! - - - - - - - - | ||
615 | ! For TKE scheme : | ||
616 | ! - - - - - - - - | ||
617 | ! WRITE(UNIT=KULOUT,FMT='( '' ADISE = '',E11.4,'' ADISI = '',E11.4,/ & | ||
618 | ! &,'' AECLS3 = '',E11.4,'' AECLS4 = '',E11.4,'' AKN = '',E11.4,/ & | ||
619 | ! &,'' ALD = '',E11.4,'' ALPHAE = '',E11.4,'' ALPHAT = '',E11.4,/ & | ||
620 | ! &,'' ECTMIN = '',E11.4,'' UCWSTAR = '',E11.4,'' UDECT = '',E11.4,/ & | ||
621 | ! &,'' USHEARM = '',E11.4,/ & | ||
622 | ! &,'' UPRETMIN= '',E11.4,'' UPRETMAX= '',E11.4,'' ARSCH = '',E11.4,/ & | ||
623 | ! &,'' ARSCQ = '',E11.4,'' ARSC1 = '',E11.4,/ & | ||
624 | ! &,'' ARSB2 = '',E11.4,'' STTBMIN = '',E11.4,/ & | ||
625 | ! &,'' ACBRPHIM= '',E11.4,'' ALMAVE = '',E11.4,'' RICRET = '',E11.4,/ & | ||
626 | ! &,'' UETEPS = '',E11.4,/ & | ||
627 | ! &,'' ABJUMIN = '',E11.4,'' RCOFLM = '',E11.4,& | ||
628 | ! &)')& | ||
629 | WRITE(UNIT=KULOUT,FMT=*) & | ||
630 | 1 | &ADISE,ADISI,AECLS3,AECLS4,AKN,ALD,ALPHAE,ALPHAT,ECTMIN,UCWSTAR,UDECT,& | |
631 | 1 | &USHEARM,UPRETMIN,UPRETMAX,ARSCH,ARSCQ,ARSC1,ARSB2,STTBMIN,& | |
632 | 2 | &ACBRPHIM,ALMAVE,RICRET,UETEPS,AJBUMIN,RCOFLM | |
633 | |||
634 | |||
635 | ! - - - - - - - - - - | ||
636 | ! Lopez Microphysics : | ||
637 | ! - - - - - - - - - - | ||
638 | 1 | WRITE(UNIT=KULOUT,FMT='('' - - - - - - - - - - - '')') | |
639 | 1 | WRITE(UNIT=KULOUT,FMT='('' Microphysics scheme '')') | |
640 | 1 | WRITE(UNIT=KULOUT,FMT='('' - - - - - - - - - - - '')') | |
641 | WRITE(UNIT=KULOUT,FMT='( '' RAUTEFR = '',E11.4 & | ||
642 | &,'' RAUTEFS = '',E11.4,'' RAUTSBET= '',E11.4,/ & | ||
643 | &,'' RNINTR = '',E11.4,'' RNINTS = '',E11.4,/ & | ||
644 | &,'' RQLCR = '',E11.4,'' RQCRNS = '',E11.4,/ & | ||
645 | &,'' RQICRMIN= '',E11.4,'' RQICRMAX= '',E11.4,/ & | ||
646 | &,'' RQLCV = '',E11.4 & | ||
647 | &,'' RQICVMIN= '',E11.4,'' RQICVMAX= '',E11.4,/ & | ||
648 | &,'' RQICRT1 = '',E11.4,'' RQICRT2 = '',E11.4,/ & | ||
649 | &,'' RQICRSN = '',E11.4,'' RACCEF = '',F11.4,/ & | ||
650 | &,'' RAGGEF = '',F11.4,'' RRIMEF = '',F11.4,/ & | ||
651 | &,'' RHCRIT1 = '',F11.6,'' RHCRIT2 = '',F11.6,/ & | ||
652 | &,'' RETAMIN = '',F11.6,'' RFACNSM = '',F11.6,/ & | ||
653 | &,'' TFVR = '',F11.6,'' TFVS = '',F11.6,/ & | ||
654 | &,'' GRHCMOD = '',F11.6,'' RHEVAP = '',F11.6,/ & | ||
655 | & )') & | ||
656 | 1 | &RAUTEFR,RAUTEFS,RAUTSBET,RNINTR,RNINTS,RQLCR,RQCRNS,RQICRMIN,RQICRMAX, & | |
657 | 1 | &RQLCV,RQICVMIN,RQICVMAX,RQICRT1,RQICRT2,RQICRSN,RACCEF,RAGGEF,RRIMEF, & | |
658 | 2 | &RHCRIT1,RHCRIT2,RETAMIN,RFACNSM,TFVR,TFVS,GRHCMOD,RHEVAP | |
659 | |||
660 | ! - - - - - - - - - - - | ||
661 | ! For Grenier scheme : | ||
662 | ! - - - - - - - - - - - | ||
663 | !WRITE(UNIT=KULOUT,FMT='('' AGRE1 = '',E11.4,'' AGRE2 = '',E11.4,'' AGREF = '',E11.4,/ & | ||
664 | ! & ,'' AGRERICR= '',E11.4,'' AGREKE = '',E11.4, & | ||
665 | ! & )')& | ||
666 | WRITE(UNIT=KULOUT,FMT=*) & | ||
667 | 1 | & AGRE1, AGRE2, AGREF, AGRERICR, AGREKE | |
668 | |||
669 | ! - - - - - - - - - - - - - - - - | ||
670 | ! For dry conv. adjustment scheme : - | ||
671 | ! - - - - - - - - - - - - - - - - | ||
672 | WRITE(UNIT=KULOUT,FMT='(& | ||
673 | & '' AJ1PEPS = '',E11.4,'' AJ1MEPS = '',E11.4,'' NAJITER = '',I3,8X & | ||
674 | & )')& | ||
675 | 1 | & AJ1PEPS, AJ1MEPS, NAJITER | |
676 | |||
677 | ! - - - - - - - | ||
678 | ! For ACCVIMPGY : | ||
679 | ! - - - - - - - | ||
680 | WRITE(UNIT=KULOUT,FMT='( '' ALFX = '',E11.4 & | ||
681 | &,'' TCTC = '',E11.4,'' TVFC = '',F11.6,/ & | ||
682 | &,'' GAMAP1 = '',F11.6,'' RKDN = '',E11.4,/ & | ||
683 | &,'' VVN = '',F11.6,'' VVX = '',F11.6,/ & | ||
684 | &,'' FENTRT = '',F11.6,'' HCMIN = '',F11.6 & | ||
685 | &,'' FQLIC = '',F11.6,'' FNEBC = '',F11.6,/ & | ||
686 | &,'' FEVAPC = '',F11.6)')& | ||
687 | 1 | & ALFX, TCTC, TVFC, GAMAP1, RKDN, VVN, VVX, FENTRT, & | |
688 | 2 | & HCMIN, FQLIC, FNEBC, FEVAPC | |
689 | |||
690 | |||
691 | 1 | WRITE(UNIT=KULOUT,FMT='('' Pseudo prognostic TKE scheme '')') | |
692 | WRITE(UNIT=KULOUT,FMT='( '' NUPTKE = '',E11.4 & | ||
693 | 1 | &'' GAMTKE = '',E11.4)') NUPTKE,GAMTKE | |
694 | |||
695 | 1 | WRITE(KULOUT,'(/'' PIL MICROPHYSICS : '')') | |
696 | WRITE(UNIT=KULOUT,FMT='('' RAUIUSTE = '',E10.4 & | ||
697 | &,'' RAUITN = '',F8.3,'' RAUITX = '',F8.3 & | ||
698 | &,'' RDPHIC = '',F8.2 & | ||
699 | &,'' GWBFAUT = '',F5.2 & | ||
700 | &,'' RWBF1 = '',F6.2,'' RWBF2 = '',F5.2 & | ||
701 | &,'' RSMDNEBX = '',F5.2 & | ||
702 | &,'' RSMDTX = '',F5.2 & | ||
703 | &,'' NSMTPA = '',I2 & | ||
704 | &,'' NSMTPB = '',I2 & | ||
705 | &,'' RCOLL = '',E11.4 & | ||
706 | &,'' RFALLL = '',E11.4 & | ||
707 | & )') & | ||
708 | 1 | &RAUIUSTE,RAUITN,RAUITX,RDPHIC,GWBFAUT,RWBF1,RWBF2,& | |
709 | 2 | &RSMDNEBX, RSMDTX, NSMTPA, NSMTPB, RCOLL, RFALLL | |
710 | !---------------------------------------- | ||
711 | 1 | WRITE(KULOUT,'(/'' PROGNOSTIC CONVECTION '')') | |
712 | if (1==0) then !!!!! A REVOIR (MPL) | ||
713 | WRITE(KULOUT,& | ||
714 | &'('' TUDBU='',E16.6,'' TUDFR='',E16.6,'' TDDBU='',E16.6,'' TDDFR='', & | ||
715 | & E16.6,& | ||
716 | & '' GCVALMX='',G10.4,'' GCVACHI='',E16.6,& | ||
717 | & '' GCVADMW='',G10.4,& | ||
718 | & '' GCVEEX ='',G10.4,& | ||
719 | & '' GCVBEE ='',G10.4,& | ||
720 | & '' GCVSQDN ='',G10.4,& | ||
721 | & '' GCVSQDR ='',G10.4,& | ||
722 | & '' GCVSQDCX ='',G10.4,& | ||
723 | & '' ECMNPI ='',G10.4,& | ||
724 | & '' GFRIC ='',G10.4& | ||
725 | & )')TUDBU,TUDFR,TDDBU,TDDFR,GCVALMX,GCVACHI,& | ||
726 | & GCVADMW,GCVEEX,GCVBEE, & | ||
727 | & GCVSQDN, GCVSQDR, GCVSQDCX, ECMNPI,GFRIC | ||
728 | else | ||
729 | 1 | print*,'>>>>> TOTO1 ', TUDBU,TUDFR,TDDBU,TDDFR,GCVALMX,GCVACHI,& | |
730 | 1 | & GCVADMW,GCVEEX,GCVBEE, & | |
731 | 2 | & GCVSQDN, GCVSQDR, GCVSQDCX, ECMNPI,GFRIC | |
732 | endif | ||
733 | WRITE(KULOUT,& | ||
734 | &'('' GDDEVF='',G10.4,'' GDDWPF='',G10.4,'' GDDBETA='',G10.4)')& | ||
735 | 1 | & GDDEVF, GDDWPF, GDDBETA | |
736 | WRITE(KULOUT,& | ||
737 | & '('' GRRINTE ='',G10.4,'' GRRMINA ='',G10.4& | ||
738 | 1 | & )') GRRINTE, GRRMINA | |
739 | ! ------------------------------------------------------------------ | ||
740 | |||
741 | !* 4. Check consistency between logical and real namelist inputs. | ||
742 | ! ------------------- | ||
743 | |||
744 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 1 times.
|
1 | IF (LMPHYS) THEN |
745 | |||
746 | ✗ | IF(LCVRA.AND.((TUDGP /= 0.0_JPRB.OR.TDDGP /= 0.0_JPRB).AND..NOT.LCVCAS)) THEN | |
747 | WRITE(UNIT=KULOUT,FMT='(A)') & | ||
748 | ✗ | & 'INCONSISTENCY BETWEEN TUDGP, TDDGP AND LCVCAS!...' | |
749 | ✗ | CALL ABOR1('TUDP<>0. OR TDDGP<>0. IMPLIES LCVCAS=T!...') | |
750 | ENDIF | ||
751 | |||
752 | ✗ | IF(.NOT.LCVLIS.AND.GCVPSI /= 0.0_JPRB) THEN | |
753 | ✗ | WRITE(UNIT=KULOUT,FMT='(A)') 'INCONSISTENCY BETWEEN LCVLIS AND GCVPSI!...' | |
754 | ✗ | CALL ABOR1('LCVLIS=F IMPLIES GCVPSI=0.!...') | |
755 | ENDIF | ||
756 | |||
757 | ENDIF | ||
758 | |||
759 | !* Consistency check for simplified physics keys | ||
760 | |||
761 |
2/8✓ Branch 0 taken 1 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 1 times.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✗ Branch 7 not taken.
|
1 | IF ((LMPHYS.OR.LSIMPH).AND.(LSMOOTHD.OR.LSMOOTHA)) THEN |
762 | |||
763 | ✗ | IF(RNLCURV == 0.0_JPRB) THEN | |
764 | ✗ | WRITE(UNIT=KULOUT,FMT='(A)') 'SMOOTHING IN SIM. PH. BUT RNLCURV = ZERO' | |
765 | ✗ | CALL ABOR1('PHYSICS AND SMOOTHING IMPLY RNLCURV /= ZERO') | |
766 | ENDIF | ||
767 | |||
768 | ENDIF | ||
769 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 1 times.
|
1 | IF (LPIL) THEN |
770 | ! Setup RHCRI profile(s) | ||
771 | ! (for LAM case... GAW set up by suecuv, called before suphy) | ||
772 | ! (For Not LAM, SULEG called by sugem called before suphy) | ||
773 | ✗ | CALL SURHCRI(KULOUT) | |
774 | ! CALL SULOCST(KULOUT) | ||
775 | ENDIF | ||
776 | |||
777 |
1/2✓ Branch 0 taken 1 times.
✗ Branch 1 not taken.
|
1 | IF (LHOOK) CALL DR_HOOK('SUPHY0',1,ZHOOK_HANDLE) |
778 | 1 | END SUBROUTINE SUPHY0 | |
779 |