GCC Code Coverage Report


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