| 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 |