suphy0.F90

Go to the documentation of this file.

!OPTIONS XOPT(NOEVAL)
SUBROUTINE suphy0(KULOUT)

!**** *SUPHY0*   - Initialize common YOMPHY0 physics atmospheric
!                  constants

!     Purpose.
!     --------
!           Initialize YOMPHY0, the common that contains the parameters
!           for the atmospheric part of the physics of the model.

!**   Interface.
!     ----------
!        *CALL* *SUPHY0(KULOUT)

!        Explicit arguments :
!        --------------------
!        KULOUT : Logical unit for the output

!        Implicit arguments :
!        --------------------
!        COMMON YOMPHY0

!     Method.
!     -------
!        See documentation

!     Externals.
!     ----------

!     Reference.
!     ----------
!        Documentation ARPEGE

!     Author.
!     -------
!        J.-F. Geleyn .

!     Modifications.
!     --------------
!        Original : 90-9-1
!        Ajout de GWDAMP (J.-F. Geleyn) : 91-2-2
!        Modified by Michel Deque 91-04-01 (param. for convect. clouds)
!        Ajout de HOBST, NPCLO1/2, XNBMAX et REVGSL, remplacement de
!           GWDCOE par GWDSE (J.-F. Geleyn, M. Deque, L. Labbe) : 92-4-8
!        Ajout de VZ0CM (E. Bazile) : 92-3-27
!        Modified by C. Castejon and E. Gerard 92-02-28 (stat. clouds)
!        Modified by R. EL Khatib : 93-04-02 Set-up defaults controled by LECMWF
!        Modified by M. Deque : 95-03-20 (USURIC and USUPRC)
!        Modified by J.F. Geleyn C. Bossuet : 95-12-20 (GWDBC)
!        Modified by Luc Gerard: 97-02-26 (entrainement de la qte de mvt horizontale).
!        Modified by J.M. Piriou: 97-02-28 (schema de nebulosite ACNEBN).
!        Modified by J.M. Piriou: 97-04-17 (valeurs par defaut).
!        Modified by J.L. Ricard : seuil (SCO) sur les precip. conv.
!                                  modif de la turbulence residuelle
!        Modified by J.M. Piriou: 97-08-21 (introduce Xu-Randall cloudiness).
!        Modified by M. Deque   : 97-10-21 (introduce GWD lift).
!        Modified by M. Deque   : 98-02-21 (new mixing length profile).
!        Modified by V. Lorant  : 98-08-05 (new mixing length profile).
!        Modified by J.M. Piriou: 98-02-11 (introduce downdrafts tuning parameters).
!        Modified by J.M. Piriou: 98-03-10 (introduce GRCVPP)
!        Modified by L. Gerard  : 98-11-30 (TUDGP, TDDGP, GCOMOD)
!        Modified by R. EL Khatib :98-12-14 Remove LRDSPIL
!        Modified by J.M. Piriou: 99-01-04 (introduce GCVADS, GCVBETA)
!        Modified by J.M. Piriou: 99-06-18 (introduce GCVPSI, GCVALFA, USURICL AND USURICE.
!                                 Change default value for GDDEVA).
!        Modified by J.M. Piriou: 2000-08-23 (new use of the Richardson critical number (USURID)).
!        Modified by J.M. Piriou: 2000-08-23 (cloud core buoyancy as a fraction of an undilute plume (GCVNU)).
!        Modified by J.M. Piriou: 2000-10-06 (exponent USURIDE).
!        Modified by E. Bazile  : 2000-12-12 CYCORA's default value.
!        Modified by J.M. Piriou: 2001-04-05 (introduce GCVMLT).
!        Modified by Y. Bouteloup:2002-06-14 (introduce RCVEVAP).
!        Modified by F. Bouyssel: 2002-06-25 (introduce UTILGUST, RRGAMMA, RRSCALE).
!        Modified by J.M. Piriou: 2002-08-19 (introduce GPBLH*).
!        Modified by D. Banciu:   2002-12-09 (introduce GCVPSIE).
!        R. El Khatib : 2001-08-07 Pruning options
!        J.M. Piriou  : 2002-01-10 set default values to operational ones.
!        03-06, move rnlcurv rnegat into yomphy0 (F. Bouyssel, C. Fischer)
!        E. Bazile : 2004-02-24 (introduce EDK).
!        M.Hamrud      01-Oct-2003 CY28 Cleaning
!        E. Bazile : 2004-06-30 (inroduces XKLM).
!        R. Brozkova : 2004-11 modifs for Xue-Randall cloud. scheme
!        P. Marquet & F.Bouyssel : 2004-08-18 (Lopez)
!        F.Bouyssel : 2005-07-18 (new defaults and new parameters for Lopez)
!        F.Bouyssel : 2006-01-25 (new defaults for Lopez!)
!        F. Vana    : 2006-01-30 tunables for pTKE.
!        R.Brozkova : 2006-03-03 : tuning constants for mixing lengths and Charnock formulae
!        A.Alias    : 2006-03-10 renaming KDN in RKDN
!        Modified by GMGEC/EAC  : 2006-05 list of modif.
!                    V. Lorant  : 99-01-05 (new mixing length profile).
!                    M. Deque   : 00-03-21 (new entrainment rate VL  ).
!                    P. Marquet : 2002-11-05 TRENTRVL=0. if not ACPBLH
!                    P. Marquet : 2004-05-27 TFVR and TFVS for ADVPRC.
!                    P. Marquet : 2004-10-13 RAUTEFR for snow (ACMICRO).
!                    P. Marquet : 2004-10-14 RAUTSBET for snow (ACMICRO).
!                    A. Alias   : 2005-06-23 (param. for ACCVIMPGDY)
!        E. Bazile & P. Marquet : 2006-04-11 AGRE1,AGRERICR,AJBUMIN,RCOFLM 
!                                  for LPBLE.
!        F.Bouyssel : 2006-10-30 RQLCV,RQICVMAX,RQICVMIN,RHEVAP
!        M. Bellus  : 03-Oct-2006 : ALARO-0 phasing (defaults for prognostic
!                                   convection, pTKE and PIL microphysics)
!        A.Alias     : 2006-07   FEVAPC added for ACCVIMPGY (JF Gueremy)
! ------------------------------------------------------------------------

USE parkind1  ,ONLY : jpim     ,jprb
USE yomhook   ,ONLY : lhook,   dr_hook

! Ce qui concerne NULNAM commente par MPL le 15.04.09
!USE YOMLUN   , ONLY : NULNAM
USE yomct0b  , ONLY : lecmwf
USE yomsimphl, ONLY : lsimph   ,lsmoothd ,lsmootha
USE yomphy   , ONLY : lcvcas   ,lcvlis   ,lcvra    ,lmphys   ,lcondwt  ,&
 & lprocld  ,lnebco   ,lnebr    ,lect     ,lnebgr   ,lcvrav3 ,lpil
USE yomphy0  , ONLY : tym      ,npclo1   ,npclo2   ,aercs1   ,aercs3   ,&
 & aercs5   ,almav    ,bedifv   ,ecmnp    ,edb      ,edc      ,&
 & edd      ,evap     ,font     ,galp     ,gccsv    ,gcomod   ,&
 & gcvads   ,gcvalfa  ,gcvbeta  ,gcvpsi   ,gddeva   ,gddsde   ,&
 & grcvpp   ,gwdamp   ,gwdbc    ,gwdcco   ,gwdcd    ,gwdlt    ,&
 & ahclpv   ,gwdse    ,hobst    ,hucoe    ,hucoe2   ,hutil    ,&
 & hutil1   ,hutil2   ,&
 & qssc     ,qsmin    ,qsnebc   ,qsnebs   ,qssusc   ,qssuss   ,&
 & qxral    ,qxrdel   ,qxrr     ,revgsl   ,rtcape   ,sco      ,&
 & sensl    ,snnbco   ,spnbco   ,sxnbco   ,tca      ,tct      ,&
 & tcw      ,tentr    ,tentrx   ,tddgp    ,tudgp    ,turb     ,&
 & tvf      ,uhdifv   ,usdmlt   ,usuprc   ,usuric   ,usurid   ,&
 & usuride  ,usurice  ,usuricl  ,vchrnk   ,vkarmn   ,vz0cm    ,&
 & xnbmax   ,ricrlm   ,xblm     ,xminlm   ,xmaxlm   ,xwsalm   ,&
 & xwsblm   ,gcvnu    ,gcvmlt   ,rcvevap  ,gpblhk0  ,gpblhra  ,&
 & utilgust ,rrgamma  ,rrscale  ,gcvpsie  ,qsusxc   ,qsusxs   ,&
 & etacut   ,rnegat   ,rnlcurv  ,qssusv   ,qxrhx    ,gcismin  ,&
 & gwdprof  ,gwdvali  ,rcin     ,edk      ,xklm     ,rphi0    ,&
 & rphir    ,qxrtgh   ,adise    ,adisi    ,aecls3   ,aecls4   ,&
 & akn      ,ald      ,alphae   ,alphat   ,ectmin   ,ucwstar  ,&
 & udect    ,ushearm  ,upretmin ,upretmax ,arsch    ,arscq    ,&
 & arsc1    ,arsb2    ,acbrphim ,almave   ,ricret   ,sttbmin  ,&
 & agreke   ,rdtfac   ,rautefs  ,rnintr   ,rnints   ,rqlcr    ,&
 & rqicrmax ,rqicrmin ,raccef   ,rrimef   ,rhcrit1  ,rhcrit2  ,&
 & retamin  ,tfvr     ,tfvs     ,rautefr  ,rautsbet ,grhcmod  ,&
 & rqicrt1  ,rqicrt2  ,rqicrsn  ,rqcrns   ,rfacnsm  ,raggef   ,&
 & rqlcv    ,rqicvmax ,rqicvmin ,rhevap   ,&
 & agre1    ,agrericr ,ajbumin  ,rcoflm   ,&
 & a0ml_au  ,a0ml_at  ,a0ml_bu  ,a0ml_bt  ,vziustar0,&
 & tentrvl  ,trentrv  ,ueteps   ,upreclp  ,&
 & arsc2    ,arsct    ,agre2    ,agref    ,&
 & aj1peps  ,aj1meps  ,najiter  ,&
 & alfx     ,tctc     ,tvfc     ,gamap1   ,rkdn      ,&
 & vvn      ,vvx      ,fentrt   ,hcmin    ,fqlic    ,fnebc    ,&
 & nuptke   ,gamtke   ,rcoll    ,rfalll   ,tddbu    ,&
 & tddfr    ,tudbu    ,tudfr    ,gcvachi  ,gcvalmx  ,gcvadmw  ,&
 & gcvbee   ,gcveex   ,ecmnpi   ,gfric    ,gcvsqdn  ,gcvsqdr  ,&
 & gcvsqdcx ,grrinte  ,grrmina  ,gddbeta  ,gddevf   ,gddwpf   ,&
 & tentrd   ,rdphic   ,gwbfaut  ,rwbf1    ,rwbf2    ,rauitn   ,&
 & rauitx   ,rauiuste ,rsmdnebx ,rsmdtx   ,nsmtpa   ,nsmtpb   ,&
 & fevapc


IMPLICIT NONE

INTEGER(KIND=JPIM),INTENT(IN)    :: KULOUT 
REAL(KIND=JPRB) :: ZHOOK_HANDLE

#include "abor1.intfb.h"
#include "posnam.intfb.h"
#include "surhcri.intfb.h"

#include "namphy0.h"

!     ------------------------------------------------------------------

!*       1.    Set default values.
!              -------------------

!        1.1 Set implicit default values

IF (lhook) CALL dr_hook('SUPHY0',0,zhook_handle)
a0ml_au=4.5_jprb
a0ml_at=5.0_jprb
a0ml_bu=3.0_jprb
a0ml_bt=0.8_jprb
aercs1=0.2726_jprb
aercs3=-0.4239_jprb
aercs5=0.3595_jprb
almav=300._jprb
bedifv=0.05_jprb
ecmnp=3000._jprb
edb=5._jprb
edc=5._jprb
edd=5._jprb
edk=1.0_jprb
etacut=1.0_jprb
evap=4.8e+06_jprb
font=2.4e+04_jprb
galp=0.53_jprb
gcismin=6.7e-05_jprb
gccsv=0._jprb
gcomod=1._jprb
gcvads=0.8_jprb
gcvalfa=4.5e-05_jprb
gcvbeta=0.2_jprb
gcvmlt=0.00016_jprb
gcvnu=2.5e-05_jprb
gcvpsi=0.5_jprb
gcvpsie=0.0_jprb
gddeva=0.25_jprb
gddsde=0.5_jprb
gpblhk0=0.25_jprb
gpblhra=4._jprb
grcvpp=1._jprb
gwdamp=0.6_jprb
gwdbc=2._jprb
gwdcco=1._jprb
gwdcd=6._jprb
gwdlt=0._jprb
gwdse=3.5e-03_jprb
gwdprof=1._jprb
gwdvali=0._jprb
hobst=3._jprb
hucoe=2._jprb
hucoe2=0.4_jprb
hutil=1.8_jprb
hutil1=-0.6_jprb
hutil2=1.1_jprb
npclo1=1
npclo2=1
qssc=1600._jprb
qsmin=1.e-4_jprb
qsnebc=26000._jprb
qsnebs=0.7_jprb
qssusc=1._jprb
qssuss=0.25_jprb
qssusv=0._jprb
qsusxc=3.3e-05_jprb
qsusxs=3.3e-05_jprb
rphi0=0._jprb
rphir=1750._jprb
qxral=10000._jprb
qxrdel=0._jprb
qxrhx=1._jprb
qxrr=0.5_jprb
qxrtgh=3.5_jprb
rcin=0._jprb
rcvevap=0._jprb
revgsl=80._jprb
rtcape=10800._jprb
sco=-20._jprb
sensl=1._jprb
snnbco=0._jprb
spnbco=3000._jprb
sxnbco=0.5_jprb
tca=1._jprb
tct=1.e-4_jprb
tcw=8.e-4_jprb
tentr=2.5e-06_jprb
tentrx=8.e-05_jprb
tddgp=0.8_jprb
tudgp=0.8_jprb
turb=1._jprb
tvf=1._jprb
tym(1)=0.92_jprb
tym(2)=0.74_jprb
tym(3)=16.6_jprb
tym(4)=10.1_jprb
tym(5)=0.08_jprb
uhdifv=8.e-04_jprb
usdmlt=1.25e+04_jprb
usuprc=0.0_jprb
usuric=1.0_jprb
usurid=0.035_jprb
usuride=1.0_jprb
usurice=0.5_jprb
usuricl=4._jprb
utilgust=0.125_jprb
vchrnk=0.021_jprb
vkarmn=0.4_jprb
vz0cm=1.5e-04_jprb
vziustar0=0._jprb
xnbmax=1._jprb
ahclpv=1000._jprb
ricrlm=0.5_jprb
rrgamma=0.8_jprb
rrscale=1.15e-4_jprb
xblm=6.5_jprb
xklm=1.0_jprb
xmaxlm=3000._jprb
xminlm=500._jprb
xwsalm=0.1_jprb
xwsblm=7.0_jprb
rnegat = -7.e-05_jprb
rnlcurv = 7.e+04_jprb

!     - - - - - - - - - - - - - -
!     The old Convective scheme :
!     - - - - - - - - - - - - - -
trentrv=1._jprb
tentrvl=-1.0_jprb

IF ( lcondwt.AND.lprocld ) THEN
  rdtfac=0.5_jprb
ELSE
  rdtfac=1.0_jprb
ENDIF

!        - - - - - - - - - -
!        Lopez Microphysics :
!        - - - - - - - - - -
rautefr=1.e-03_jprb
rautefs=1.e-03_jprb
rautsbet=0.025_jprb
rnintr=8.e+06_jprb
rnints=2.e+06_jprb
rqlcr=2.e-04_jprb
rqicrmax=0.3e-04_jprb
rqicrmin=0.2e-06_jprb
rqlcv=2.e-04_jprb
rqicvmax=0.3e-04_jprb
rqicvmin=0.2e-06_jprb
rqicrt1=-80._jprb
rqicrt2=30._jprb
rqicrsn=0.5_jprb
rqcrns=0.03_jprb
raccef=1._jprb
raggef=0.2_jprb
rrimef=1._jprb
rhevap=0.0_jprb
rhcrit1=0.5_jprb
rhcrit2=0.91_jprb
retamin=0.4_jprb
rfacnsm=1.4_jprb
tfvr=5.0_jprb
tfvs=0.6_jprb
grhcmod=0.3_jprb


!--------------------------
! PIL MICROPHYSICS
!--------------------------
rdphic=10000._jprb
gwbfaut=15._jprb
! two constants for ACPLUIE_PROG:
rwbf1=300._jprb
rwbf2=4._jprb
!
rauitn=233.15_jprb
rauitx=263.15_jprb
rauiuste=0.025_jprb
rsmdnebx=0.2_jprb
rsmdtx=1.0_jprb
nsmtpa=2
nsmtpb=3
rcoll=6.9e-03_jprb
rfalll=1.0_jprb
!--------------------------
! Prognostic convection physical parameters
!--------------------------
tudbu=0.5_jprb
tddbu=0.5_jprb
tudfr=0.0012_jprb
tddfr=0.0006_jprb
gcvalmx=0.95_jprb
! Pas d'activite historique: defaut tres haut,
! pour que KUO joue tout seul:
gcvachi=1.e9_jprb
gcvadmw =0
! Explicit entrainment:
gcvbee=0._jprb ! 0.2_JPRB
gcveex=1._jprb
ecmnpi=3000._jprb
gfric=-1.0_jprb !1.E-3_JPRB
! Sqeezing:
gcvsqdn=0.01_jprb
gcvsqdr=0.8_jprb
gcvsqdcx=1.0_jprb ! No squeezing
! PRECIPITATING AREA PARAMETERS (aplpar):
grrinte=2._jprb
grrmina=1.e-5_jprb
! DD explicit detrainment:
gddbeta=0.2_jprb
! ACMODO DD PARAMETERS:
gddevf=0.5_jprb
gddwpf=0._jprb
! DD ENTRAINMENT RATE:
tentrd=1.e-4_jprb

!     - - - - - - - -
!     For TKE scheme :
!     - - - - - - - -
adise=-0.5_jprb
adisi=1.5_jprb
aecls3=3.75_jprb
aecls4=0.3_jprb
akn=0.2_jprb
ald=1.4_jprb
alphae=1.0_jprb
alphat=1.0_jprb
ectmin=1.e-10_jprb

ucwstar=1.0_jprb/3._jprb
udect=5._jprb
ushearm=1.e-04_jprb
upretmin=60000._jprb
upretmax=97500._jprb
arsch=4._jprb
arscq=1.2_jprb
arsc1=2.0_jprb/(3._jprb*arsch*arscq)
arsb2=3._jprb*arsch*arsc1/2.0_jprb
acbrphim=2.2_jprb
almave=0._jprb
ricret=0.195_jprb

sttbmin=sqrt(3._jprb)
! A minimum value for the adimentional jump in boyancy : d(Theta)/Theta
ajbumin=0.005_jprb
! The "Master Length" is equal to "RCOFLM*Z_PBL"
rcoflm=0.085_jprb


! TKE (P.Marquet)
ueteps=1.0_jprb 

!     - - - - - - - - - - - - - - - -
!     For dry conv. adjustment scheme :
!     - - - - - - - - - - - - - - - -
aj1meps=0.99_jprb
aj1peps=10.0_jprb
najiter=30

!     - - - - - - - - - - - - - - - - - - - - - - -
!     For Grenier (2000) top-PBL entrainment scheme :
!     - - - - - - - - - - - - - - - - - - - - - - -
agre1=0.16_jprb
agre2=15._jprb
agref=0.8_jprb
agreke=5.0_jprb
agrericr=50._jprb

! Pseudo prognostic TKE scheme
nuptke=0.52_jprb
gamtke=0.5_jprb



!     - - - - - - -
!     For ACCVIMPGY :
!     - - - - - - -
alfx=10.e-02_jprb
tctc=1.60e-04_jprb
tvfc=1._jprb
gamap1=1.5_jprb
rkdn=30.e-06_jprb
vvn=0.0_jprb
vvx=-45._jprb
fentrt=2.5_jprb
hcmin=0.0_jprb
fqlic=2.5_jprb
fnebc=25.0_jprb
fevapc=3.5_jprb

!        1.2 Modify default values according to LECMWF

IF (lecmwf) THEN
ELSE
ENDIF

!     ------------------------------------------------------------------

!*       2.    Modify default values.
!              ----------------------

! Ce qui concerne NAMPHY0 commente par MPL le 15.04.09
!CALL POSNAM(NULNAM,'NAMPHY0')
!READ(NULNAM,NAMPHY0)

IF(.NOT.(lnebco.AND.(lnebr.OR.lnebgr.OR.lect)) .AND.lcvrav3) THEN
  WRITE(unit=kulout,fmt='(A)') ' '
  WRITE(unit=kulout,fmt='(A)') ' !'
  WRITE(unit=kulout,fmt='(A)') ' ! TENTRVL is set to 0. in SUPHY0  !!'
  WRITE(unit=kulout,fmt='(A)') ' !'
  WRITE(unit=kulout,fmt='(A)') ' '
  tentrvl = 0.0_jprb
ENDIF
!     ------------------------------------------------------------------

!*       3.    Print final values.
!              -------------------

WRITE(unit=kulout,fmt='('' COMMON YOMPHY0 '')')
WRITE(unit=kulout,fmt='('' AERCS1 = '',E11.4,'' AERCS3 = '',E11.4 &
 & ,'' AERCS5 = '',E11.4,'' ALMAV = '',E11.4,'' ECMNP = '',E11.4 &
 & ,'' EDB = '',E11.4,/,'' EDC = '',E11.4,'' EDD = '',E11.4 &
 & ,'' EDK = '',E11.4,'' ETACUT = '',E11.4 &
 & ,'' EVAP = '',E11.4,'' FONT = '',E11.4,'' GWDAMP = '',E11.4 &
 & ,'' GWDSE = '',E11.4,'' GWDBC = '',E11.4,/&
 & ,'' GWDCD = '',E11.4,/&
 & ,'' GWDPROF = '',E11.4,'' GWDVALI = '',E11.4,/&
 & ,'' HUCOE = '',E11.4,'' HUTIL = '',E11.4 ,'' HUTIL1 = '',E11.4 ,'' HUTIL2 = '',E11.4 &
 & ,'' VCHRNK = '',E11.4,'' VKARMN = '',E11.4 &
 & ,'' SNNBCO = '',E11.4,/,'' SPNBCO = '',E11.4 &
 & ,'' SXNBCO = '',E11.4,'' HOBST = '',E11.4,'' NPCLO1 = '',I4 &
 & ,'' NPCLO2 = '',I4,'' RCIN = '',E11.4,'' RCVEVAP = '',E11.4 &
 & ,'' REVGSL = '',E11.4,/&
 & ,'' RTCAPE= '',E11.4,'' GCOMOD= '',E11.4,/&
 & ,'' UHDIFV = '',E11.4,'' VZ0CM = '',E11.4,'' VZIUSTAR0 = '',E11.4 &
 & ,'' XNBMAX = '',E11.4,'' GALP = '',E11.4,'' SENSL = '',E11.4 &
 & ,'' TCA = '',E11.4,/,'' TCT = '',E11.4,'' TCW = '',E11.4 &
 & ,'' TURB = '',E11.4,'' TVF = '',E11.4,/&
 & ,'' TYM = '',5E11.4,'' QSMIN = '',E11.4,/&
 & ,'' QSSC = '',E11.4,'' RPHI0 = '',E11.4,'' RPHIR = '',E11.4 &
 & ,'' BEDIFV = '',E11.4,'' SCO = '',E11.4 &
 & ,'' USDMLT = '',E11.4,/&
 & ,'' GDDEVA = '',E11.4,'' GDDSDE = '',E11.4,/&
 & ,'' GWDCCO = '',E11.4 &
 & ,'' HUCOE2 = '',E11.4,/&
 & ,'' TENTR = '',E11.4,'' TENTRX = '',E11.4 &
 & ,'' TUDGP = '',E11.4,'' TDDGP = '',E11.4 &
 & ,'' GRCVPP = '',E11.4 &
 & )')&
 & aercs1,aercs3,aercs5,almav,ecmnp,edb,edc,edd,edk,etacut,evap,font,&
 & gwdamp,gwdse,gwdbc,gwdcd,gwdprof,gwdvali,&
 & hucoe,hutil,hutil1,hutil2,vchrnk,vkarmn,snnbco,spnbco,sxnbco,hobst,&
 & npclo1,npclo2,rcin,rcvevap,revgsl,rtcape,gcomod,&
 & uhdifv,vz0cm,vziustar0,xnbmax,galp,sensl,tca,tct,tcw,&
 & turb,tvf,tym,qsmin,qssc,rphi0,rphir,bedifv,sco,usdmlt,&
 & gddeva,gddsde,gwdcco,&
 & hucoe2,&
 & tentr,tentrx,tudgp,tddgp,grcvpp  
WRITE(unit=kulout,fmt='('' USUPRC = '',E11.4,'' USURIC = '',E11.4 &
 & ,'' QSNEBC = '',E11.4,'' QSNEBS = '',E11.4 &
 & ,'' QSSUSC = '',E11.4,'' QSSUSS = '',E11.4,'' QSSUSV = '',E11.4 &
 & ,'' QSUSXC = '',E11.4,'' QSUSXS = '',E11.4 &
 & ,'' GCCSV  = '',E11.4,/&
 & ,'' QXRAL  = '',E11.4,'' QXRDEL = '',E11.4 &
 & ,'' QXRHX  = '',E11.4,'' QXRR   = '',E11.4,'' QXRTGH = '',E11.4,/&
 & ,'' GWDLT  = '',E11.4,'' AHCLPV = '',E11.4,/&
 & ,'' GCVADS = '',E11.4,'' GCVBETA= '',E11.4 &
 & ,'' RICRLM = '',E11.4,'' XBLM   = '',E11.4,'' XKLM = '',E11.4 &
 & ,'' XMINLM = '',E11.4,'' XMAXLM = '',E11.4 &
 & )')&
 & usuprc,usuric,qsnebc,qsnebs,qssusc,qssuss,qssusv,qsusxc,qsusxs,gccsv,qxral &
 & ,qxrdel,qxrhx,qxrr,qxrtgh,gwdlt,ahclpv,gcvads,gcvbeta,ricrlm,xblm,xklm,xminlm &
 & ,xmaxlm  
WRITE(unit=kulout,fmt='('' XWSALM = '',E11.4,'' XWSBLM = '',E11.4 &
 & ,'' GCVALFA= '',E11.4,'' GCVPSI = '',E11.4,'' GCVPSIE = '',E11.4 &
 & ,'' USURICL= '',E11.4,'' USURICE= '',E11.4,/ &
 & ,'' USURID= '',E11.4,'' USURIDE= '',E11.4,'' GCVNU= '',E11.4 &
 & ,'' GCVMLT= '',E11.4,'' GPBLHK0= '',E11.4,'' GPBLHRA = '',E11.4 &
 & ,'' UTILGUST= '',E11.4,'' RRGAMMA= '',E11.4,'' RRSCALE = '',E11.4 &
 & )')&
 & xwsalm,xwsblm,gcvalfa,gcvpsi,gcvpsie,usuricl,usurice,usurid,usuride,gcvnu,gcvmlt,gpblhk0 &
 & ,gpblhra,utilgust,rrgamma,rrscale  
WRITE(unit=kulout,fmt='('' RNEGAT == '',E11.4,'' RNLCURV = '',E11.4 &
 & )') rnegat,rnlcurv  
WRITE(unit=kulout,fmt='('' RDTFAC == '',E11.4 )') rdtfac
WRITE(unit=kulout,fmt='(9(5(A,E11.4),/))') 'GCISMIN=',gcismin
WRITE(unit=kulout,fmt='('' A0ML_AU == '',E11.4,'' A0ML_AT == '',E11.4 &
 &,'' A0ML_BU == '',E11.4,'' A0ML_BT == '',E11.4 )') a0ml_au,a0ml_at&
 &,a0ml_bu,a0ml_bt

!     - - - - - - - -
!     For TKE scheme :
!     - - - - - - - -
! WRITE(UNIT=KULOUT,FMT='( '' ADISE   = '',E11.4,'' ADISI   = '',E11.4,/ &
!  &,'' AECLS3  = '',E11.4,'' AECLS4  = '',E11.4,'' AKN     = '',E11.4,/ &
!  &,'' ALD     = '',E11.4,'' ALPHAE  = '',E11.4,'' ALPHAT  = '',E11.4,/ &
!  &,'' ECTMIN  = '',E11.4,'' UCWSTAR = '',E11.4,'' UDECT   = '',E11.4,/ &
!  &,'' USHEARM = '',E11.4,/ &
!  &,'' UPRETMIN= '',E11.4,'' UPRETMAX= '',E11.4,'' ARSCH   = '',E11.4,/ &
!  &,'' ARSCQ   = '',E11.4,'' ARSC1   = '',E11.4,/ &
!  &,'' ARSB2   = '',E11.4,'' STTBMIN = '',E11.4,/ &
!  &,'' ACBRPHIM= '',E11.4,'' ALMAVE  = '',E11.4,'' RICRET  = '',E11.4,/ &
!  &,'' UETEPS  = '',E11.4,/ &
!  &,'' ABJUMIN = '',E11.4,'' RCOFLM  = '',E11.4,&
!      &)')&
   WRITE(unit=kulout,fmt=*) &
  &adise,adisi,aecls3,aecls4,akn,ald,alphae,alphat,ectmin,ucwstar,udect,&
  &ushearm,upretmin,upretmax,arsch,arscq,arsc1,arsb2,sttbmin,&
  &acbrphim,almave,ricret,ueteps,ajbumin,rcoflm


!       - - - - - - - - - -
!       Lopez Microphysics :
!       - - - - - - - - - -
WRITE(unit=kulout,fmt='('' - - - - - - - - - - - '')')
WRITE(unit=kulout,fmt='(''  Microphysics scheme  '')')
WRITE(unit=kulout,fmt='('' - - - - - - - - - - - '')')
WRITE(unit=kulout,fmt='( '' RAUTEFR = '',E11.4   &
 &,'' RAUTEFS = '',E11.4,'' RAUTSBET= '',E11.4,/ &
 &,'' RNINTR  = '',E11.4,'' RNINTS  = '',E11.4,/ &
 &,'' RQLCR   = '',E11.4,'' RQCRNS  = '',E11.4,/ &
 &,'' RQICRMIN= '',E11.4,'' RQICRMAX= '',E11.4,/ &
 &,'' RQLCV   = '',E11.4 &
 &,'' RQICVMIN= '',E11.4,'' RQICVMAX= '',E11.4,/ &
 &,'' RQICRT1 = '',E11.4,'' RQICRT2 = '',E11.4,/ &
 &,'' RQICRSN = '',E11.4,'' RACCEF  = '',F11.4,/ &
 &,'' RAGGEF  = '',F11.4,'' RRIMEF  = '',F11.4,/ &
 &,'' RHCRIT1 = '',F11.6,'' RHCRIT2 = '',F11.6,/ &
 &,'' RETAMIN = '',F11.6,'' RFACNSM = '',F11.6,/ &
 &,'' TFVR    = '',F11.6,'' TFVS    = '',F11.6,/ &
 &,'' GRHCMOD = '',F11.6,'' RHEVAP  = '',F11.6,/ &
 &    )') &
 &rautefr,rautefs,rautsbet,rnintr,rnints,rqlcr,rqcrns,rqicrmin,rqicrmax, &
 &rqlcv,rqicvmin,rqicvmax,rqicrt1,rqicrt2,rqicrsn,raccef,raggef,rrimef, &
 &rhcrit1,rhcrit2,retamin,rfacnsm,tfvr,tfvs,grhcmod,rhevap

!     - - - - - - - - - - -
!     For Grenier scheme :
!     - - - - - - - - - - -
!WRITE(UNIT=KULOUT,FMT='('' AGRE1   = '',E11.4,'' AGRE2   = '',E11.4,'' AGREF   = '',E11.4,/ &
! & ,'' AGRERICR= '',E11.4,'' AGREKE  = '',E11.4, &
! & )')&
WRITE(unit=kulout,fmt=*) &
 & agre1,    agre2,   agref,      agrericr, agreke

!     - - - - - - - - - - - - - - - -
!     For dry conv. adjustment scheme :                                -
!     - - - - - - - - - - - - - - - -
WRITE(unit=kulout,fmt='(&
 & '' AJ1PEPS = '',E11.4,'' AJ1MEPS = '',E11.4,'' NAJITER = '',I3,8X &
 & )')&
 & aj1peps, aj1meps, najiter

!     - - - - - - -
!     For ACCVIMPGY :
!     - - - - - - -
WRITE(unit=kulout,fmt='( '' ALFX    = '',E11.4  &
 &,'' TCTC    = '',E11.4,'' TVFC    = '',F11.6,/ &
 &,'' GAMAP1  = '',F11.6,'' RKDN    = '',E11.4,/ &
 &,'' VVN     = '',F11.6,'' VVX     = '',F11.6,/ &
 &,'' FENTRT  = '',F11.6,'' HCMIN   = '',F11.6  &
 &,'' FQLIC   = '',F11.6,'' FNEBC   = '',F11.6,/ &
 &,'' FEVAPC  = '',F11.6)')&
 & alfx, tctc, tvfc, gamap1, rkdn, vvn, vvx, fentrt, &
 & hcmin, fqlic, fnebc, fevapc


WRITE(unit=kulout,fmt='(''  Pseudo prognostic TKE scheme  '')')
WRITE(unit=kulout,fmt='( '' NUPTKE = '',E11.4  &
 &'' GAMTKE = '',E11.4)') nuptke,gamtke

WRITE(kulout,'(/'' PIL MICROPHYSICS : '')')
WRITE(unit=kulout,fmt='('' RAUIUSTE = '',E10.4  &
 &,'' RAUITN = '',F8.3,'' RAUITX = '',F8.3  &
 &,'' RDPHIC = '',F8.2   &
 &,'' GWBFAUT = '',F5.2 &
 &,'' RWBF1 = '',F6.2,'' RWBF2 = '',F5.2 &
 &,'' RSMDNEBX = '',F5.2 &
 &,'' RSMDTX = '',F5.2 &
 &,'' NSMTPA = '',I2 &
 &,'' NSMTPB = '',I2 &
 &,'' RCOLL = '',E11.4 &
 &,'' RFALLL = '',E11.4 &
 &    )') &
 &rauiuste,rauitn,rauitx,rdphic,gwbfaut,rwbf1,rwbf2,&
 &rsmdnebx, rsmdtx, nsmtpa, nsmtpb, rcoll, rfalll
!----------------------------------------
WRITE(kulout,'(/'' PROGNOSTIC CONVECTION '')')
if (1==0) then     !!!!! A REVOIR (MPL)
WRITE(kulout,&
 &'('' TUDBU='',E16.6,'' TUDFR='',E16.6,'' TDDBU='',E16.6,'' TDDFR='', &
 & E16.6,&
 & '' GCVALMX='',G10.4,'' GCVACHI='',E16.6,&
 & '' GCVADMW='',G10.4,&
 & '' GCVEEX ='',G10.4,&
 & '' GCVBEE ='',G10.4,&
 & '' GCVSQDN ='',G10.4,&
 & '' GCVSQDR ='',G10.4,&
 & '' GCVSQDCX ='',G10.4,&
 & '' ECMNPI ='',G10.4,&
 & '' GFRIC ='',G10.4&
 & )')tudbu,tudfr,tddbu,tddfr,gcvalmx,gcvachi,&
 &    gcvadmw,gcveex,gcvbee, &
 &    gcvsqdn, gcvsqdr, gcvsqdcx, ecmnpi,gfric
else
  print*,'>>>>> TOTO1 ', tudbu,tudfr,tddbu,tddfr,gcvalmx,gcvachi,&
 &    gcvadmw,gcveex,gcvbee, &
 &    gcvsqdn, gcvsqdr, gcvsqdcx, ecmnpi,gfric
endif
WRITE(kulout,&
 &'('' GDDEVF='',G10.4,'' GDDWPF='',G10.4,'' GDDBETA='',G10.4)')&
 & gddevf, gddwpf, gddbeta
WRITE(kulout,&
 & '('' GRRINTE ='',G10.4,'' GRRMINA ='',G10.4&
 & )') grrinte, grrmina
!     ------------------------------------------------------------------

!*       4.    Check consistency between logical and real namelist inputs.
!              -------------------

IF (lmphys) THEN

  IF(lcvra.AND.((tudgp /= 0.0_jprb.OR.tddgp /= 0.0_jprb).AND..NOT.lcvcas)) THEN
    WRITE(unit=kulout,fmt='(A)') &
     & 'INCONSISTENCY BETWEEN TUDGP, TDDGP AND LCVCAS!...'  
    CALL abor1('TUDP<>0. OR TDDGP<>0. IMPLIES LCVCAS=T!...')
  ENDIF

  IF(.NOT.lcvlis.AND.gcvpsi /= 0.0_jprb) THEN
    WRITE(unit=kulout,fmt='(A)') 'INCONSISTENCY BETWEEN LCVLIS AND GCVPSI!...'
    CALL abor1('LCVLIS=F IMPLIES GCVPSI=0.!...')
  ENDIF

ENDIF

!*  Consistency check for simplified physics keys

IF ((lmphys.OR.lsimph).AND.(lsmoothd.OR.lsmootha)) THEN

  IF(rnlcurv == 0.0_jprb) THEN
    WRITE(unit=kulout,fmt='(A)') 'SMOOTHING IN SIM. PH. BUT RNLCURV = ZERO'
    CALL abor1('PHYSICS AND SMOOTHING IMPLY RNLCURV /= ZERO')
  ENDIF

ENDIF
IF (lpil) THEN
! Setup RHCRI profile(s)
! (for LAM case... GAW set up by suecuv, called before suphy)
! (For Not LAM, SULEG called by sugem called before suphy)
  CALL surhcri(kulout)
! CALL SULOCST(KULOUT)
ENDIF

IF (lhook) CALL dr_hook('SUPHY0',1,zhook_handle)
END SUBROUTINE suphy0