su0phy.F90

Go to the documentation of this file.

SUBROUTINE su0phy(KULOUT)

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

!USE YOMLUN   , ONLY : NULOUT   ,NULNAM  MPL 14.04.09
USE yomlun   , ONLY : nulout
! Ce qui concerne LSLPHY commente par MPL 24.11.08
!USE YOMSLPHY , ONLY : LSLPHY
USE yomphy   , ONLY : nbiter   ,noir, ndpsfi, nphyrep ,lmphys   ,lreasur  ,&
 & lcape    ,lcondwt  ,lcvpp    ,lcvdd    ,lhuneg   ,&
 & lneige   ,lrnumx   ,lclsatur,l2phys    ,lcvra    ,lgwd     ,&
 & lgwdc    ,lhmto    ,lnebco   ,lnebn    ,lnebr    ,lqxrtgh, lhucn,&
 & lnebt    ,lnd2diff  ,lozone   ,lray     ,lrayfm   ,lrayfm15 ,&
 & lrrmes   ,lsfhyd   ,lsnv     ,lsolv    ,lfgel    ,&
 & lsrcon   ,lsrcont  ,lslc     ,lrrgust  ,lrelaxw  ,&
 & laerosea ,laerolan ,laerosoo ,laerodes ,laerovol ,laerosul ,lrelaxt  ,&
 & lo3abc   ,lstra    ,lstras   ,lthermo  ,lvdif    ,&
 & lraylu   ,lrews    ,lrprox   ,lrmix    ,lrstab   ,&
 & lrautoev ,lraypl   ,lcvlis   ,lcvcas   ,lvgsn    ,&
 & lnebnxr  ,lfpcor   ,lnoias   ,cgmixlen ,lprgml   ,lglt     ,&
 & lnewd    ,lrtpp    ,lrtdl    ,ldifcons ,lect     ,&
 & lcvpgy   ,lprocld  ,levapp   ,lcollec  ,lptke    ,l3mt     ,&
 & lcvpro   ,lcddpro  ,lscmf    ,lvoigt   ,lvfull   ,&
 & lnsmlis, lphcdpi   ,&
 & nphy     ,jphyec   ,jphymf   ,jphyaro  , & 
 & lajucv   ,lpble    ,lnebgr   ,lnebgy   ,&
 & lbccond  ,lcvrav3  ,lz0hsrel ,lblvar   ,&
 & ladjcld  ,lautoneb ,lssd     ,lcvppkf  ,lectfl ,&
 & lpil     ,lphspsh  ,lsmrot   ,lsmnimbt ,&
 & lsmtps   ,l1drhcri ,lgwrhcri ,nsmtbot   , &
 & nsmdneb  ,nprag    ,nprac    ,nprri    ,lstrapro, lnewstat
USE yomarphy , ONLY : lmpa    ,lmicro   ,lturb    ,&
 & lmse  ,lkfbconv ,lkfbd    ,lkfbs    ,lusechem ,&
 & lorilam  ,lrdust, lbuflux ,ccoupling 
USE yoephy   , ONLY : lephys    ,&
 & lecond   ,lecumf   ,ledcld   ,leevap   ,legwdg   ,&
 & leozoc   ,leqngt   ,leradi   ,lerads   ,&
 & leshcv   ,lesice   ,lesurf   ,levdif   ,&
 & lagphy   ,lepcld   ,leco2diu ,&
 & leo3ch   ,lbud23   ,lemethox ,lera40   ,lecurr   ,lvdftrac ,&
 & leocwa   ,leocco   ,leocsa   ,lmftrac  ,lerain   ,le4alb   ,&
 & rthrfrti ,nephys_pcfull, lemwave
USE yoewcou  , ONLY : nstpw    ,rsoutw   ,rnortw   ,rdegrew  ,&
 & lwcou    ,lwcou2w
! Ce qui concerne RCLDTOPP commente par MPL 24.11.08
!USE YOECLDP  , ONLY : RCLDTOPP  
USE yophlc   , ONLY : alpha    ,ah0      ,ustarl   ,&
 & ustars   ,alandz0  ,aseaz0   ,lsphlc   ,lvdflc   ,&
 & lsdrlc   ,lczdeb   ,lzmcon   ,lkexp    ,lvdfds   ,&
 & lsdrds  
USE yomct0   , ONLY : nconf    ,&
 & lretcfou ,lwrtcfou, larome,lfpos,lpc_full
USE yomct0b  , ONLY : lecmwf
! Tous les YOM* ci-dessous commentes par MPL 24.11.08
!USE YOMINI   , ONLY : NEINI
!USE YOMDFI   , ONLY : NEDFI
!USE YOMVRTL  , ONLY : L131TL
!USE YOPHNC   , ONLY : LETRAJP  ,LETRAJPT ,LERADI2  ,LERADS2   ,&
! & LERADSW2 ,LERADN2 ,LERADFL2 ,LEDCLD2  ,LENCLD2   ,&
! & LEVAPLS2 ,LEVDIF2 ,LEGWDG2  ,LECUMF2  ,LECUBM2   ,&
! & LECOND2  ,LEQNGT2 ,LESURF2  ,LEKPERT  ,LTRACLNPH
!USE YOEPHLI  , ONLY : LENOPERT
!USE YOMNCL   , ONLY : LNCLIN   ,LREGCL
USE yomsimphl, ONLY : lsimph   ,ltrajps  ,ltrajpst ,&
 & lsmoothd ,lsmootha ,lsmoothb ,lcvrasp  ,&
 & lgwdsp   ,lraysp   ,lstrasp  ,lvdifsp  ,lrrmessp ,lclouds  
USE yomrcoef , ONLY : lrcoef    ,ltladdia ,lglobrad
!USE YOMIOP   , ONLY : NPCKFT95 ,NEXPBT95
!USE YOMDYNA  , ONLY : LGWADV
!USE YOMCOAPHY   , ONLY : NPHYINT

!**** *SU0PHY*   - Initialize common YOxPHY controlling physics

!     Purpose.
!     --------
!           Initialize YOxPHY, the common that includes the
!           basic switches for the physics of the model.

!**   Interface.
!     ----------
!        *CALL* *SU0PHY(KULOUT) from SU0YOMA

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

!        Implicit arguments :
!        --------------------

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

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

!     Reference.
!     ----------
!        ECMWF Research Department documentation of the IFS

!        or

!        Documentation ARPEGE (depending on which physics will be used)

!     Author.
!     -------
!        J.-J. Morcrette                    *ECMWF*
!        J.-F. Geleyn for the ARPEGE rewriting.

!     Modifications.
!     --------------
!        Original : 91-11-12
!        Modified 92-02-22 by M. Deque (tests of consistency with *YOMDPHY*)
!        Modified by R. EL Khatib : 93-04-02 Set-up defaults controled by LECMWF
!        Modified 94-02-28 by M.  Deque  : Shallow convection clouds
!        Modified 93-10-28 by Ph. Dandin : FMR scheme with MF physics
!        Modified 93-08-24 by D. Giard (test of consistency with digital filter)
!        Modified by M. Hamrud    : 93-06-05 Make use of LECMWF for ECMWF
!        Modified 95-11-27 by M. Deque (2nd call to APLPAR)
!        Modified 96-01-10 by M. Janiskova (logical switches for simpl.ph.par.)
!        Modified 97-02-28 by J.M. Piriou (cloudiness scheme switch LNEBN)
!        Modified 97-04-17 by J.M. Piriou (default values)
!        Modified by F. Rabier    : 96-09-25 Full physics set-up for 801 job
!        Modified by G. Hello     : 97-07-31 Full MF physics set-up for 801
!        Modified by M. Deque     : 97-05-25 Frozen FMR
!        Modified by E.Bazile     : 97-11-18 Soil freezing (LFGEL)
!        Modified by M. Deque     : 98-01-05 Cleaning of NDPSFI
!        Modified by M. Janiskova : 98-11-18 Full MF physics set-up for 131
!                                 : 99-02-21 Set-up for radiation coef.
!        Modified by C. Jakob     : 98-04    Methane oxidation
!        Modified by T. Bergot    : 98-08  Full MF physics set-up for 601
!        Modified by E.Bazile     : 99-02-12 Superficial soil freezing (LFGELS)
!        Modified by L. Gerard    : 98-12-07 LSRCON
!        Modified by J.M. Piriou  : 99-04-19 Moon radiation
!        Modified 99-02-17 by K. YESSAD: options LRETCFOU, LWRTCFOU.
!        Modified by J.M. Piriou  : 99-06-18 LCVLIS and LCVCAS
!        Modified by J.M. Piriou  : 99-07-07 Introduce reproductibility in physics (NPHYREP).
!        Modified by D. Giard     : 2000-10-25 LVGSN (snow and vegetation)
!        Modified by E. Bazile    : 2000-11-12 CYCORA's default value.      
!        Modified by F. Bouyssel  : 2001-03-03 LRRMESSP     
!        Modified by J.M. Piriou  : 2001-11-27 LSRCONT
!        Modified by J.M. Piriou  : 2002-11-15 LNEBNXR
!        Modified by F. Bouyssel  : 2002-06-25 LCVPPLIS and LRRGUST
!        R. El Khatib : 2001-08-07 Pruning options
!        J.M. Piriou  : 2002-01-10 set default values to operational ones.
!        Modified by A.Beljaars   : 2002-11-12 LECURR (Ocean current)
!        Modified by Y. Bouteloup : 2002-03-05 LO3ABC
!        Modified 08-2002 C. Smith : use "w" as prognostic variable in the
!         semi-lag advection of vertical divergence in the NH model. 
!        Modified by F. Bouyssel  : 2002-12-18 Cleaning of LCVPPLIS
!        Modified by E. Bazile      : 2003-02-13 LFPCOR
!        Modified by E. Bazile      : 2003-02-18 LNOIAS
!        Modified by M. Janiskova : 2003-05 set-up for ECMWF stat.cloud scheme
!        M.Hamrud      01-Oct-2003 CY28 Cleaning
!        P. Viterbo   ECMWF   03-12-2004  Include user-defined RTHRFRTI
!        Modified by Y Seity      : 2004-11-16 For AROME setup, default values 
!            for namarphy keys,read namarphy, switch off arp/ald physics keys
!        Modified by R. Brozkova : 2004-11 modifs for Xu-Randall cloud. scheme
!        P. Marquet and F. Bouyssel : 2004-08-18 (Lopez) 
!        G. Hello : 2005-04-25 Lground initialization (surfex and arome)
!        Y. Seity : 2005-09-25 LRDUST, LORILAM and LUSECHEM initialisation
!                 for AROME
!        T. Kovacic : 2006-03-17 LPHCDPI, NPHY, LBUFLUX; for DDH and BUDGET
!        Modified by F. Bouyssel  : 2005-01-25 Change default of LNSMLIS
!        D.Salmond     22-Nov-2005 Mods for coarser/finer physics
!        P. Lopez      14-02-2006  Added switch LTRACLNPH for including  
!                                  tracers in linearized physics
!        R.Brozkova : pre ALARO 0 modset: mixing lengths computation
!        Modified by GMGEC/EAC   : 2006-03 list of modif.
!                 P. Marquet  : 99-01-18 Dry Conv. Adj.     (LAJUCV)
!                 P. Marquet  : 02-02-14 YOMPHY
!                 P. Marquet  : 02-06-18 LBCCOND (Becht/Chab ACCOND)
!                 P. Marquet  : 02-08-30 LCVRAV3    (old ACCVIMP_V3)
!                     and new cloud model : 2006-03-03
!        Modified by E. Bazile    : 2006-04-11 Add LPBLE in case of LECT
!        Modified by E. Bazile    : 2006-04-20 Add LADJCLD,LCVPPKF, LECTFL
!                 JF. Gueremy : LZ0HSREL initialised (used in ACHMT)
!        Modified by F. Bouyssel  : 2006-10-30 Add LAUTONEB, LSSD
!        Modified by F. Vana     : 2006-01-30 LPTKE
!        Modifed by D. Banciu    : 2006-08-31 LND2DIFF
!        M. Bellus : 28-Sep-2006 ALARO-0 phasing: L3MT,LCVPRO,LCDDPRO,LSCMF,
!                    LVOIGT,LVFULL,LPIL,LPHSPSH,LSMROT,LSMNIMBT,LSMTPS,
!                    L1DRHCRI,LGWRHCRI,NSMTBOT,NSMDNEB,NPRAG,NPRAC,NPRRI,
!                    LSTRAPRO,LNEWSTAT ; removed duplicity in LPBLE
!        JJMorcrette 20060721 MODIS albedo
!     ------------------------------------------------------------------

IMPLICIT NONE

INTEGER(KIND=JPIM),INTENT(IN)    :: KULOUT 

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

INTEGER(KIND=JPIM) :: IERR
REAL(KIND=JPRB) :: ZHOOK_HANDLE

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

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

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

#include "namphy.h"
!#include "namarphy.h"
#include "naephy.h"
#include "naphlc.h"
!#include "namtrajp.h"
!#include "namsimphl.h"
!#include "namrcoef.h"

!     ------------------------------------------------------------------
print*,'Dans SUOPHY ', kulout

IF (lhook) CALL dr_hook('SU0PHY',0,zhook_handle)

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

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

!*       1.1.1 Set default values for Meteo-France physics
!              -------------------------------------------

lmphys=.false.
lreasur=.true.

cgmixlen='Z'
lprgml=.false.
lcape=.false.
lcondwt=.false.
lcvcas=.false.
lcvlis=.false.
lcvpp=.false.
lcvppkf=.false.
lcvdd=.false.
lhuneg=.true.
lneige=.true.
lrnumx=.false.
lclsatur=.false.
lvoigt=.false.
lvfull=.false.
lrrgust=.false.
l2phys=.false.
lo3abc=.false.
laerosea=.false.
laerolan=.false.
laerosoo=.false.
laerodes=.false.
laerovol=.false.
laerosul=.false.
lrelaxt=.false.
lrelaxw=.false.
le4alb=.false.
lglt=.false.
lnewd=.false.
ldifcons=.false.
lect=.false.
lptke=.false.
lectfl=.false.
lpble=.false.
lcvpgy=.false.
l3mt=.false.

lcvra=.false.
lgwdc=.false.
lgwd=.false.
lhmto=.false.
lnebco=.false.
lnebn=.false.
lnebnxr=.false.
lnebr=.false.
lnebt=.false.
lnd2diff=.false.
lqxrtgh=.false.
lhucn=.false.
lozone=.false.
lray=.false.
lraylu=.false.
lrews=.false.
lrprox=.false.
lrmix=.false.
lrstab=.false.
lrautoev=.false.
lrtpp=.false.
lrtdl=.false.
lraypl=.false.
lrayfm=.false.
lrayfm15=.false.
lrrmes=.false.
lsfhyd=.false.
lsnv=.false.
lsolv=.false.
lfgel=.false.
lsrcon=.false.
lsrcont=.false.
lslc=.false.
lstra=.false.
lstras=.false.
lthermo=.false.
lvdif=.false.
lvgsn=.false.
lfpcor=.false.
lnoias=.false.
lphcdpi=.false.
lblvar=.false.

lz0hsrel=.false.

nbiter=2
ndpsfi=0
nphyrep=1
noir=0

! ---------------------------------------------------
! ALARO-0 (cloud)
! ---------------------------------------------------
lpil=.false.
lstrapro=.false.
lnewstat=.true.
lphspsh=.false. ! Pseudo Historic Surface Precip Sensible Heat Flux
lsmrot=.false.
lsmtps=.false.
lsmnimbt=.false.
l1drhcri=.false.
lgwrhcri=.false.
nsmtbot=0 ! interpolate
nsmdneb=2 ! gradient limitation
nprag=1
nprac=1
nprri=1
! ---------------------------------------------------
! ALARO-0 (prognostic convection)
! ---------------------------------------------------
lcvpro=.false.
lcddpro=.false. 
lscmf=.false.

! - - - - - - - - - - - - - - - - - - - - - - - - - -
! Cloud and precipitation prognostic scheme (Lopez) :
! - - - - - - - - - - - - - - - - - - - - - - - - - -

lprocld=.false.
levapp=.true.
lcollec=.true.
lnsmlis=.true.
ladjcld=.true.
lautoneb=.false.
lssd=.false.

!AROME physics
lmpa=.false.
lmicro=.false.
lturb=.false.
lmse=.false.
lkfbconv=.false.
lkfbd=.false.
lkfbs=.false.
lusechem=.false.
lorilam=.false.
lrdust=.false.
lbuflux=.true.
ccoupling='E'

! - - - - - - - - - - -
! Module YOMPHY :
! - - - - - - - - - - -

lajucv=.false.
lnebgr=.false.
lnebgy=.false.
lbccond=.false.

lcvrav3=.false.

!*    1.1.2  Set default values for simplified physical parametrization
!            of Meteo-France
!     -----------------------------------------------------------------

lsimph=.false.
ltrajps=.false.
ltrajpst=.false.
lsmoothd=.false.
lsmootha=.false.
lsmoothb=.false.
lclouds=.false.

lcvrasp=.false.
lgwdsp=.false.
lraysp=.false.
lstrasp=.false.
lvdifsp=.false.
lrrmessp=.false.

lrcoef=.false.
ltladdia=.false.
lglobrad=.false.

!*    1.2.1  Set default values for ECMWF physics package
!            --------------------------------------------

lephys=.false.
lagphy=.true.

lecond=.false.
lepcld=.false.
lecumf=.false.
ledcld=.false.
leevap=.true.
legwdg=.false.
leozoc=.false.
leqngt=.false.
leradi=.false.
lerads=.false.
leshcv=.false.
lesice=.true.
lesurf=.false.
levdif=.false.
leocwa=.false.
leocco=.false.
leocsa=.false.
lemethox=.false.
lera40=.false.
lecurr=.false.
lvdftrac=.true.
lmftrac=.true.
lerain=.false.
le4alb=.false.
rthrfrti=0.0_jprb
!NPHYINT=0
le4alb=.false.

!-------------------------------------------------------
! pressure above which cloud scheme is not called
! !!!WARNING!!! this has to be in the part of the domain 
! where the pure pressure level grid is used, otherwise
! the code is not bit-reproducible!!!
! Don't call the cloud scheme for pressures lower than 1hPa
!RCLDTOPP=100.0_JPRB

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

!*     1.2.2  Set-up linearized physical parametrization of ECMWF
!             ---------------------------------------------------

!LETRAJP = .FALSE.
!LETRAJPT= .FALSE.
!LERADI2 = .FALSE.
!LERADS2 = .FALSE.
!LERADSW2= .FALSE.
!LERADN2 = .FALSE.
!LERADFL2= .FALSE.
!LEDCLD2 = .FALSE.
!LENCLD2 = .FALSE.
!LEVAPLS2= .FALSE.
!LEVDIF2 = .FALSE.
!LEGWDG2 = .FALSE.
!LECUMF2 = .FALSE.
!LECUBM2 = .FALSE.
!LECOND2 = .FALSE.
!LEQNGT2 = .FALSE.
!LESURF2 = .FALSE.
!LEKPERT = .FALSE.
!LNCLIN  = .FALSE.
!LREGCL  = .FALSE.
!LTRACLNPH = .FALSE.

! No perturbation of surface arrays
!LENOPERT = .TRUE.

!*       Packing parameters
!        -------------------
!NPCKFT95 = 1
!NEXPBT95 = 6

!     LOGICAL FOR THE VERT DIFF SCHEME VDIFLCZ USED IN CONF 601

lsphlc  = .false.
lvdflc  = .false.
lvdfds  = .true.
lsdrlc  = .true.
lsdrds  = .false.
lczdeb  = .false.
lzmcon  = .true.
lkexp   = .true.
alpha   = 3._jprb
ah0     = 1000.0_jprb
ustarl  = 0.5_jprb
ustars  = 0.2_jprb
alandz0 = 0.05_jprb
aseaz0  = 0.0005_jprb

nstpw=2
!     NSTPW=30

lwcou=.false.
lwcou2w=.false.

!  Setup for 3 degree resolution wave model

!     RNORTW= 72.0
!     RSOUTW=-63.0
!     RDEGREW=3.0

!  Setup for 1.5 degree resolution wave model

rnortw= 81.0_jprb
rsoutw=-81.0_jprb
rdegrew=1.5_jprb

!        1.3 Modify default values according to LECMWF

IF (lecmwf) THEN
  lepcld=.true.
  levdif=.true.
  leocwa=.false.
  leocco=.false.
  leocsa=.true.
  lesurf=.true.
  lecond=.false.
  lecumf=.true.
  leevap=.true.
  legwdg=.true.
  leozoc=.true.
  leqngt=.true.
  leradi=.true.
  lerads=.true.
  lesice=.true.
  ledcld=.true.
  leo3ch=.false.
  leco2diu=.false.
  lemethox=.true.
  IF(nconf == 1) THEN
    lephys=.true.
  ENDIF
  IF(nconf == 131) THEN
!   LERADI2 = .TRUE.
!   LERADS2 = .TRUE.
!   LEVDIF2 = .TRUE.
!   LEGWDG2 = .TRUE.
!   LECUMF2 = .TRUE.
!   LECOND2 = .TRUE.
!    IF(L131TL) THEN
!      LEPHYS=.TRUE.
!      LECOND=.TRUE.
!      LEPCLD=.FALSE.
!    ENDIF
    lsphlc=.true.
    lvdflc=.true.
    lsdrlc=.true.
    lzmcon=.true.
    lkexp =.true.
  ENDIF
  IF (nconf == 401 .OR. nconf == 501) THEN
    lephys = .true.
    lecond = .true.
    lepcld = .false.

!   LETRAJP = .TRUE.
!   LERADI2 = .TRUE.
!   LERADS2 = .TRUE.
!   LEVDIF2 = .TRUE.
!   LEGWDG2 = .TRUE.
!   LECUMF2 = .TRUE.
!   LECOND2 = .TRUE.
  ENDIF

  IF(nconf == 601.OR.nconf == 801) THEN

!      Full physics when computing the trajectory

    lephys=.true.

!      Simple scheme for TL and ADJ

    lsphlc  = .true.
    lvdflc  = .true.
    lvdfds  = .true.
    lsdrlc  = .true.
    lsdrds  = .true.
    lczdeb  = .false.
    lzmcon  = .true.
    lkexp   = .true.
    alpha   = 3._jprb
    ah0     = 1000.0_jprb
    ustarl  = 0.5_jprb
    ustars  = 0.2_jprb
    alandz0 = 0.05_jprb
    aseaz0  = 0.0005_jprb

  ENDIF
  nphy = jphyec
ELSE
  IF ((nconf == 1).OR.(nconf == 131).OR.(nconf == 401)&
     & .OR.(nconf == 501).OR.(nconf == 601)&
     & .OR.(nconf == 701).OR.(nconf == 801)) THEN  
    lagphy=.false.
    lcvcas=.true.
    lcvdd=.true.
    lcvlis=.true.
    lcvpp=.true.
    lsrcont=.false.
    lcvra=.true.
    lfgel=.true.
    lgwd=.true.
    lhmto=.true.
    lmphys=.true.
    lnebn=.true.
    lnebt=.false.
    lnebnxr=.false.
    lneige=.true.
    lray=.true.
    lraylu=.true.
    lraypl=.true.
    lrrmes=.true.
    lsfhyd=.true.
    lsolv=.true.
    lsrcon=.true.
    lstra=.true.
    lthermo=.true.
    lvdif=.true.
    nphy = jphymf
    IF(larome.AND..NOT.lfpos) THEN
      !extinction of MF phsic's keys for arome's run
      !lfpos is added for 927 confs for example.
      lcvra=.false.
      lgwd=.false.
      lgwdc=.false.
      lnebco=.false.
      lnebn=.false.
      lnebr=.false.
      lnebt=.false.
      lozone=.false.
      lray=.false.
      lraylu=.false.
      lrews=.false.
      lraypl=.false.
      lrayfm=.false.
      lrayfm15=.false.
      lrrmes=.false.
      lsfhyd=.false.
      lsnv=.false.
      lsolv=.false.
      lfgel=.false.
      lstra=.false.
      lstras=.false.
      lvdif=.false.
      !initialisation of AROME's physic ones
      lmpa=.true.
      lmicro=.false.
      lturb=.false.
      lmse=.false.
      lkfbconv=.false.
      lusechem=.false.
      lkfbd=.false.
      lkfbs=.false.
      lusechem=.false.
      lorilam=.false.
      lrdust=.false.
      nphy = jphyaro 
    ENDIF
  ENDIF
ENDIF

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

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

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

! Ce qui concerne NAMARPHY commente par MPL le 24.11.08
!CALL POSNAM(NULNAM,'NAMARPHY')
!READ(NULNAM,NAMARPHY)
!IF(LMSE) THEN
!  LSOLV=.FALSE.
!  LFGEL=.FALSE.
!  LSNV=.FALSE.
!  LVGSN=.FALSE.
!  WRITE(NULOUT,'('' INFO -old isba-lsolv,lfgel,lsnv,lvgsn- reset to .F. '', &
!                 &'' when using lground i.e. surfex '')')
!ENDIF

nephys_pcfull=3

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

IF(lpc_full)THEN
  IF(nephys_pcfull < 2 .OR.nephys_pcfull > 3)THEN
    CALL abor1(' SU0PHY: NEPHYS_PCFULL')
  ENDIF
ELSE
  nephys_pcfull=3
ENDIF

IF(     nconf == 201.OR.nconf == 202 &
   & .OR.nconf == 421.OR.nconf == 422 &
   & .OR.nconf == 521.OR.nconf == 522  )THEN  
  lreasur=.false.
ENDIF
!IF(.NOT.(LEPHYS.OR.LMPHYS)) THEN
!  LAGPHY=.FALSE.
!  LSLPHY=.FALSE.
!  WRITE(NULOUT,'('' INFO - LSLPHY RESET TO .FALSE. '')')
!ENDIF
IF(lepcld) THEN
  lecond=.false.
ENDIF
IF(lagphy.AND.ndpsfi /= 0) THEN
  ndpsfi=0
  WRITE(nulout,'('' WARNING - NDPSFI RESET TO 0 '')')
ENDIF

! Commente par MPL 24.11.08
!CALL POSNAM(NULNAM,'NAPHLC')
!READ(NULNAM,NAPHLC)

! Commente par MPL 24.11.08
!CALL POSNAM (NULNAM,'NAMTRAJP')
!READ (NULNAM,NAMTRAJP)

! Commente par MPL 24.11.08
!CALL POSNAM(NULNAM,'NAMSIMPHL')
!READ(NULNAM,NAMSIMPHL)

! Commente par MPL 24.11.08
!CALL POSNAM(NULNAM,'NAMRCOEF')
!READ(NULNAM,NAMRCOEF)

!-------------------------------------------------
! Initialize profile extractions for the Single Column Model .
!-------------------------------------------------

CALL suscm(kulout)

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

!*       3.    Do tests of consistency.
!              ------------------------

ierr=0

! * Test that adequate physics is activated when at least one of the
!   options LRETCFOU, LWRTCFOU, LRCOEF,is set to .TRUE.
IF( lwrtcfou .AND. ((.NOT.lmphys).OR.lephys.OR.lagphy) ) THEN
  ierr=ierr+1
  WRITE(kulout,'(1X,A,I3,A)') ' ! SU0PHY: ERROR NR ',ierr,' !!!'
  WRITE(kulout,'(1X,A,A,A)') ' SU0PHY: options ',&
   & 'LWRTC, LWRTCFOU ',&
   & 'require the following options for physics:'  
  WRITE(kulout,'(1X,A)') '         - LMPHYS=.TRUE. '
  WRITE(kulout,'(1X,A)') '         - LEPHYS=.FALSE. '
  WRITE(kulout,'(1X,A)') '         - LAGPHY=.FALSE. '
  IF(.NOT.lraysp) THEN
    WRITE(kulout,'(1X,A,A,A)') ' SU0PHY: option ',&
     & ' LWRTCFOU ',&
     & 'require the following option for simplified physics:'  
    WRITE(kulout,'(1X,A)') '         - LRAYSP=.TRUE. '
  ENDIF
ENDIF
IF( lretcfou .AND. ((.NOT.(lmphys.OR.lsimph)).OR.lephys.OR.lagphy) ) THEN
  ierr=ierr+1
  WRITE(kulout,'(1X,A,I3,A)') ' ! SU0PHY: ERROR NR ',ierr,' !!!'
  WRITE(kulout,'(1X,A,A,A)') ' SU0PHY: option ',&
   & ' LRETCFOU ',&
   & 'require the following options for physics:'  
  WRITE(kulout,'(1X,A)') '         - LMPHYS=.TRUE. or LSIMPH=.TRUE.'
  WRITE(kulout,'(1X,A)') '         - LEPHYS=.FALSE. '
  WRITE(kulout,'(1X,A)') '         - LAGPHY=.FALSE. '
ENDIF
IF( lrcoef .AND. ((.NOT.lmphys).OR.lephys.OR.lagphy) ) THEN
  ierr=ierr+1
  WRITE(kulout,'(1X,A,I3,A)') ' ! SU0PHY: ERROR NR ',ierr,' !!!'
  WRITE(kulout,'(1X,A,A,A)') ' SU0PHY: options ',&
   & 'LRCOEF ',&
   & 'require the following options for physics:'  
  WRITE(kulout,'(1X,A)') '         - LMPHYS=.TRUE. '
  WRITE(kulout,'(1X,A)') '         - LEPHYS=.FALSE. '
  WRITE(kulout,'(1X,A)') '         - LAGPHY=.FALSE. '
  IF(.NOT.lraysp) THEN
    WRITE(kulout,'(1X,A,A,A)') ' SU0PHY: options ',&
     & 'LRCOEF ',&
     & 'require the following option for simplified physics:'  
    WRITE(kulout,'(1X,A)') '         - LRAYSP=.TRUE. '
  ENDIF
ENDIF

! * Test that options LRETCFOU, LWRTCFOU are set to .FALSE.
!   when at least LRCOEF=.TRUE. 
IF( (lrcoef) .AND. (lretcfou.OR.lwrtcfou) ) THEN
  ierr=ierr+1
  WRITE(kulout,'(1X,A,I3,A)') ' ! SU0PHY: ERROR NR ',ierr,' !!!'
  WRITE(kulout,'(1X,A,A,A)') ' SU0PHY: options ',&
   & 'LRCOEF ',&
   & 'require that:'  
  WRITE(kulout,'(1X,A)') '         - LRETCFOU=.FALSE. '
  WRITE(kulout,'(1X,A)') '         - LWRTCFOU=.FALSE. '
ENDIF

! * Test that option LSRCONT is set to .FALSE. when LSRCON=.FALSE. 
IF( (lsrcont) .AND. (.NOT.lsrcon) ) THEN
  ierr=ierr+1
  WRITE(kulout,'(1X,A,I3,A)') ' ! SU0PHY: ERROR NR ',ierr,' !!!'
  WRITE(kulout,'(1X,A,A,A)') ' SU0PHY: options ',&
   & 'LSRCONT ',&
   & 'require that:'  
  WRITE(kulout,'(1X,A)') '         - LSRCON=.TRUE. '
ENDIF

! * Use NDPSFI=1 when not available or inconsistent with some other options?
!IF(LGWADV .AND. (NDPSFI==1)) THEN
  ! * these two options actually can run together on an informatic
  !   point of view, but the current assumptions done with LGWADV=T
  !   are that a particle which is on the Earth surface remains
  !   on the Earth surface; this is equivalent to assume that "etadot_surf"
  !   is always zero. This condition is satisfied when NDPSFI=0
  !   but not when NDPSFI=1 which gives a non-zero value to "etadot_surf".
  !   So one forbids combination "LGWADV .AND. (NDPSFI==1)" which leads
  !   to inconsistencies in the model at the surface.
!  IERR=IERR+1
!  WRITE(KULOUT,'(1X,A,I3,A)') ' ! SU0PHY: ERROR NR ',IERR,' !!!'
!  WRITE(KULOUT,'(1X,A,A,A)') ' SU0PHY: options ',&
!   & 'NDPSFI=1 ',&
!   & 'require that:'  
!  WRITE(KULOUT,'(1X,A)') '         - LGWADV=.FALSE. '
!ENDIF

! TESTS FOR ALARO-0 (cloud)
IF (lpil.AND..NOT.lcondwt) THEN
  ierr=ierr+1
  WRITE(kulout,'(1X,A,I3,A)') ' ! SU0PHY: ERROR NR ',ierr,' !!!'
  WRITE(kulout,'(1X,A)') ' SU0PHY: LPIL requires LCONDWT '
ENDIF

! TESTS FOR ALARO-0 (prognostic convection)
IF (lcvpro) THEN
  lcddpro=.true.
  IF (nconf /= 1)THEN
    WRITE(kulout,'(1X, A)') 'WARNING: LCVPRO requires NCONF=1'
    lcvpro=.false.
    lcddpro=.false.
  ENDIF
ENDIF

IF (lphspsh .AND. .NOT.(l3mt.OR.lstrapro)) THEN
  ! * For the time being, LPHSPSH requires L3MT.OR.LSTRAPRO (missing code
  !   in CPTEND in the other cases, to be coded later).
  ierr=ierr+1
  WRITE(kulout,'(1X,A,I3,A)') ' ! SU0PHY: ERROR NR ',ierr,' !!!'
  WRITE(kulout,'(1X,A)') .OR.' SU0PHY: LPHSPSH requires L3MTLSTRAPRO '
ENDIF

! * Use simplified physics with LMPA=T?
IF(lmpa.AND.lsimph) THEN
  ! * Simplified physics is not consistent with AROME physics
  !   concerning most of the physical parameterizations.
  ierr=ierr+1
  WRITE(kulout,'(1X,A,I3,A)') ' ! SU0PHY: ERROR NR ',ierr,' !!!'
  WRITE(kulout,'(1X,A)') ' SU0PHY: LMPA and LSIMPH cannot be both to T.'
ENDIF


IF (ierr >= 1) THEN
  CALL flush(kulout)
  CALL abor1(' SU0PHY: ABOR1 CALLED')
ENDIF

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

!*       4.    Print final values.
!              -------------------

WRITE(unit=kulout,fmt='('' COMMON YOMPHY '')')
WRITE(unit=kulout,fmt='('' LMPHYS= '',L5,'' LREASUR = '',L5)')lmphys,lreasur
if (1==0) then       !!!!! A REVOIR (MPL)
WRITE(unit=kulout,fmt='('' LCONDWT='',L5,'' LCVPP = '',L5,&
 & '' LCVPPKF = '',L5,&
 & '' LNEIGE = '',L5,'' LRNUMX = '',L5,'' LCLSATUR = '',L5,'' LHUNEG = '',L5,&
 & '' LVOIGT = '',L5,'' LVFULL = '',L5,&
 & '' L2PHYS = '',L5,&
 & '' LCVDD = '',L5,'' LCAPE= '',L5,'' LSRCON = '',L5,&
 & '' LCVLIS = '',L5,'' LCVCAS= '',L5,'' LSRCONT= '',L5,'' LO3ABC= '',L5,  &
 & '' LAEROSEA = '',L5,'' LAEROLAN = '',L5,'' LAEROSOO = '',L5, &
 & '' LAERODES = '',L5,'' LAEROVOL = '',L5,'' LAEROSUL = '',L5, &
 & '' LRELAXT  = '',L5,'' LRELAXW  = '',L5,&
 & '' LSLC = '',L5,'' LRRGUST= '',L5, &
 & '' LNEWD = '',L5,'' LGLT= '',L5,'' LDIFCONS= '',L5, &
 & '' LECT = '',L5,'' LECTFL = '',L5,'' LPBLE = '',L5,'' LCVPGY= '',L5, &
 & '' LPTKE= '',L5,'' L3MT= '',L5 &
 & )')&
 & lcondwt,lcvpp,lcvppkf,lneige,lrnumx,lclsatur,lhuneg,lvoigt,lvfull,l2phys,lcvdd,&
 & lcape,lsrcon,lcvlis,lcvcas,lsrcont,lo3abc,laerosea,laerolan,laerosoo, &
 & laerodes,laerovol,laerosul,lrelaxt,lrelaxw,lslc,lrrgust,lnewd,lglt,&
 & ldifcons,lect,lectfl,lpble,lcvpgy,lptke,l3mt 
else
  print*,'OKL'
  print*,'LOGICS',&
 & lcondwt,lcvpp,lcvppkf,lneige,lrnumx,lclsatur,lhuneg,lvoigt,lvfull,l2phys,lcvdd,&
 & lcape,lsrcon,lcvlis,lcvcas,lsrcont,lo3abc,laerosea,laerolan,laerosoo, &
 & laerodes,laerovol,laerosul,lrelaxt,lrelaxw,lslc,lrrgust,lnewd,lglt,&
 & ldifcons,lect,lectfl,lpble,lcvpgy,lptke,l3mt 
endif
WRITE(unit=kulout,fmt='('' LCVRA = '',L5 ,'' LFPCOR = '',L5&
 & ,'' LNOIAS = '',L5 &
 & ,'' LGWD = '',L5,'' LHMTO = '',L5,'' LRAY = '',L5 &
 & ,'' LSFHYD = '',L5,'' LSTRA = '',L5,'' LTHERMO = '',L5,/&
 & ,'' LVDIF = '',L5,'' LNEBCO = '',L5,'' LNEBT = '',L5,/&
 & ,'' LRRMES = '',L5,'' LOZONE = '',L5,'' LNEBR = '',L5,/&
 & ,'' LSNV = '',L5,'' LSOLV = '',L5,'' LFGEL = '',L5,/&
 & ,'' LVGSN ='',L5,'' LND2DIFF ='',L5,/&
 & ,'' LNEBN = '',L5,'' LNEBNXR = '',L5,'' LQXRTGH = '',L5,'' LSTRAS = '',L5,/&
 & ,'' LHUCN = '',L5,'' LGWDC     = '',L5,'' LRAYFM = '',L5,'' LRAYFM15= '',L5,/&
 & ,'' LRAYLU  = '',L5,'' LREWS  = '',L5,'' LRPROX = '',L5 /&
 & ,'' LRMIX   = '',L5,'' LRSTAB = '',L5,''LRAUTOEV = '',L5 /&
 & ,'' LRTPP = '',L5,'' LRTDL = '',L5,'' LRAYPL = '',L5 &
 & )')&
 & lcvra,lfpcor,lnoias,lgwd,lhmto,lray,lsfhyd,lstra,lthermo,lvdif,lnebco,&
 & lnebt,lrrmes,lozone,lnebr,lsnv,lsolv,lfgel,lvgsn,lnd2diff,&
 & lnebn,lnebnxr,lqxrtgh,lstras,lhucn,lgwdc,&
 & lrayfm,lrayfm15,lraylu,lrews,&
 & lrprox,lrmix,lrstab,lrautoev,lrtpp,lrtdl,lraypl
WRITE(unit=kulout,fmt='('' NBITER = '',I2,'' NDPSFI = '',I2,'' NPHYREP = '',I2)') &
 & nbiter,ndpsfi,nphyrep  
WRITE(unit=kulout,fmt='('' NOIR = '',I2)')noir
!WRITE(UNIT=KULOUT,FMT='('' NPHYINT = '',I2)')NPHYINT

WRITE(unit=kulout,fmt='('' - - - - - - - - - -'')')
WRITE(unit=kulout,fmt='('' COMMON YOMPHY'')')
WRITE(unit=kulout,fmt='('' - - - - - - - - - -'')')
WRITE(unit=kulout,fmt='(&
 & '' LECT    = '',L5,'' LAJUCV  = '',L5,/, &
 & '' LCVRAV3 = '',L5, &
 & '' LNEBGR  = '',L5,'' LNEBGY  = '',L5,'' LBCCOND = '',L5 &
 & )') lect,    lajucv,  &
 & lcvrav3, lnebgr, lnebgy, lbccond

WRITE(unit=kulout,fmt='('' - - - - - - - - -'')')
WRITE(unit=kulout,fmt='('' MICROPHYSICS KEY '')')
WRITE(unit=kulout,fmt='('' - - - - - - - - -'')')
WRITE(unit=kulout,fmt='(  '' LPROCLD = '',L5 &
 &,'' LEVAPP  = '',L5,'' LCOLLEC = '',L5,'' LNSMLIS  = '',L5 &
 &,'' LADJCLD  = '',L5,'' LAUTONEB  = '',L5,'' LSSD = '',L5)')&
 &lprocld,levapp,lcollec,lnsmlis,ladjcld,lautoneb,lssd

WRITE(unit=kulout,fmt='('' ALARO-0 cloud '')')
WRITE(unit=kulout,fmt='('' LPIL = '',L5,'' LSTRAPRO = '',L5 &
  &,'' LNEWSTAT = '',L5 &
  &,'' LPHSPSH = '',L5 ,'' LSMROT = '',L5 &
  &,'' LSMTPS = '',L5 ,'' LSMNIMBT = '',L5 &
  &,'' L1DRHCRI = '',L5 ,'' LGWRHCRI = '',L5 &
  &,'' NSMTBOT = '',I2,'' NSMDNEB = '',I2 &
  &,'' NPRAG = '',I2,'' NPRAC = '',I2 &
  &,'' NPRRI = '',I2 &
  &) ')lpil,lstrapro,lnewstat,lphspsh,lsmrot,lsmtps, lsmnimbt, &
  &    l1drhcri, lgwrhcri,&
  &    nsmtbot,nsmdneb, nprag,nprac,nprri

WRITE(unit=kulout,fmt='('' ALARO-0 prognostic convection  '')')
WRITE(unit=kulout,&
  &fmt='('' LCVPRO ='',L1,'' LCDDPRO ='',L1,'' LSCMF ='',L1)') &
  &    lcvpro,lcddpro,lscmf

WRITE(unit=kulout,fmt='('' COMMON YOMSIMPHL '')')
WRITE(unit=kulout,fmt='('' LSIMPH= '',L5,'' LTRAJPS = '',L5 &
 & ,'' LTRAJPST = '',L5 &
 & ,'' LSMOOTHD = '',L5,'' LSMOOTHA = '',L5,'' LSMOOTHB = '',L5 &
 & ,'' LCLOUDS = '',L5 )')&
 & lsimph,ltrajps,ltrajpst,lsmoothd,lsmootha, &
 & lsmoothb,lclouds  
WRITE(unit=kulout,fmt='('' LCVRASP = '',L5,'' LGWDSP = '',L5 &
 & ,'' LRAYSP = '',L5,'' LSTRASP = '',L5,'' LVDIFSP = '',L5 &
 & ,'' LRRMESSP = '',L5)')&
 & lcvrasp,lgwdsp,lraysp,lstrasp,lvdifsp,lrrmessp  
WRITE(unit=kulout,fmt='('' COMMON YOMRCOEF '')')
WRITE(unit=kulout,fmt='('' LRCOEF= '',L5 &
 & ,'' LTLADDIA = '',L5,'' LGLOBRAD = '',L5)')&
 & lrcoef,ltladdia,lglobrad  

! TEST OF CONSISTENCY FOR RADIATION SCHEMES

IF(lmphys) THEN
  IF((lray.AND.lrayfm).OR.(lray.AND.lrayfm15).OR.(lrayfm15.AND.lrayfm)) THEN
    WRITE(nulout,'('' WARNING - 2 RADIATION SCHEMES... '')')
    CALL abor1('SU0PHY: ABOR1 CALLED')
  ENDIF
ENDIF

!     TESTS OF CONSISTENCY INSIDE YOMPHY

IF(lstras.AND.lrnumx)THEN
  WRITE(nulout,fmt='('' ACPLUIS AND LRNUMX ARE NOT COMPATIBLE''&
   & ,''      FOR THE MOMENT'')')  
  CALL abor1('SU0PHY: ABOR1 CALLED')
ENDIF

!     TESTS OF CONSISTENCY BETWEEN YOMPHY AND INITIALIZATION

!IF (.NOT.LMPHYS.AND..NOT.LEPHYS) THEN
! Digital Filter Initialisation
!  IF ((NEINI == 2.OR.NEINI == 4).AND.NEDFI >= 2) THEN
!    IF (NEDFI <= 5) THEN
!      WRITE(NULOUT,FMT='('' YOMPHY AND YEMDFI ARE NOT COMPATIBLE'')')
!      CALL ABOR1(' DIABATIC DFI WITHOUT PHYSICS !')
!    ELSE
!      WRITE(NULOUT,FMT='('' CAUTION : FORWARD DFI WITHOUT PHYSICS'')')
!    ENDIF
!  ENDIF
!ENDIF

WRITE(unit=kulout,fmt='('' COMMON YOMARPHY '')')
WRITE(unit=kulout,fmt='('' LMPA = '',L5, &
 & '' LMICRO = '',L5,'' LTURB = '',L5, &
 & '' LMSE = '',L5,'' LKFBCONV = '',L5,&
 & '' LKFBD = '',L5,'' LKFBS = '',L5,&
 & '' LUSECHEM = '',L5,'' LORILAM = '',L5,'' LRDUST = '',L5)')&
 & lmpa,lmicro,lturb,lmse,lkfbconv,lkfbd,lkfbs,lusechem,&
 & lorilam,lrdust

WRITE(unit=kulout,fmt='('' COMMON YOEPHY '')')
WRITE(unit=kulout,fmt='('' LEPHYS = '',L5, &
 & '' LECOND = '',L5,'' LECUMF = '',L5 &
 & )')&
 & lephys,lecond,lecumf  
!WRITE(NULOUT,'("RCLDTOPP=",F10.2)') RCLDTOPP
WRITE(unit=kulout,fmt='('' LEDCLD = '',L5,'' LEGWDG = '',L5 &
 & ,'' LEOZOC = '',L5,'' LEQNGT = '',L5 &
 & ,'' LEO3CH = '',L5)')&
 & ledcld,legwdg,leozoc,leqngt,leo3ch  
WRITE(unit=kulout,fmt='('' LERADI = '',L5 &
 & ,'' LESHCV = '',L5,'' LESURF = '',L5,'' LEVDIF = '',L5 &
 & ,'' LEOCWA = '',L5,'' LEOCCO = '',L5,'' LEOCSA = '',L5 &
 & ,'' LECURR = '',L5,'' RTHRFRTI = '',F7.2 &
 & )')&
 & leradi,leshcv,lesurf,levdif,leocwa,leocco,leocsa,lecurr,rthrfrti 
WRITE(unit=kulout,fmt='('' LEPCLD = '',L5,'' LEMETHOX= '',L5 &
 & ,'' LE4ALB = '',L5 &
 & )')&
 & lepcld,lemethox,le4alb  

! Do not allow storage of trajectory in TL physics if tracers 
! are to be included in linearized physics (temporary) 
!IF (LETRAJP .AND. LETRAJPT .AND. LTRACLNPH) THEN
!  WRITE(UNIT=KULOUT,FMT='('' CAUTION : STORAGE OF TRAJECTORY IN LINEAR. PHYSICS NOT YET '')')
!  WRITE(UNIT=KULOUT,FMT='('' IMPLEMENTED FOR TRACERS: LETRAJPT WILL BE RESET TO .FALSE. '')')
!ENDIF

!IF (LETRAJP) THEN
!  WRITE(UNIT=KULOUT,FMT='('' ======================== '')')
!  WRITE(UNIT=KULOUT,FMT='('' LINEAR PHYSICS ACTIVATED '')')
!  WRITE(UNIT=KULOUT,FMT='('' ======================== '')')
!  WRITE(UNIT=KULOUT,FMT='('' LERADI2 = '',L4 &
!   & ,'' LERADS2 = '',L4 &
!   & ,'' LERADSW2= '',L4 &
!   & ,'' LERADN2 = '',L4 &
!   & ,'' LERADFL2= '',L4 &
!   & )')&
!   & LERADI2,LERADS2,LERADSW2,LERADN2,LERADFL2  
!  WRITE(UNIT=KULOUT,FMT='('' LEDCLD2 = '',L4 &
!   & ,'' LENCLD2 = '',L4 &
!   & ,'' LEVAPLS2= '',L4 &
!   & ,'' LREGCL  = '',L4 &
!   & ,'' LNCLIN  = '',L4 &
!   & )')&
!   & LEDCLD2,LENCLD2,LEVAPLS2,LREGCL,LNCLIN
!  WRITE(UNIT=KULOUT,FMT='('' LECUMF2 = '',L4 &
!   & ,'' LECUBM2 = '',L4 &
!   & ,'' LECOND2 = '',L4    &
!   & )')&
!   & LECUMF2,LECUBM2,LECOND2  
!  WRITE(UNIT=KULOUT,FMT='('' LEVDIF2 = '',L4 &
!   & ,'' LEGWDG2 = '',L4 &
!   & ,'' LEQNGT2 = '',L4 &
!   & ,'' LESURF2 = '',L4 &
!   & )')&
!   & LEVDIF2,LEGWDG2,LEQNGT2,LESURF2  
!  WRITE(UNIT=KULOUT,FMT='('' LEKPERT = '',L4,'' LENOPERT = '',L4)')&
!   & LEKPERT,LENOPERT
!  WRITE(UNIT=KULOUT,FMT='('' ======================== '')')
!  WRITE(UNIT=KULOUT,FMT='('' LETRAJPT = '',L4)')LETRAJPT
!  WRITE(UNIT=KULOUT,FMT='('' ======================== '')')
!  WRITE(UNIT=KULOUT,FMT='('' LTRACLNPH = '',L4)')LTRACLNPH
!  WRITE(UNIT=KULOUT,FMT='('' ======================== '')')
!ELSE
!  WRITE(UNIT=KULOUT,FMT='('' ============================ '')')
!  WRITE(UNIT=KULOUT,FMT='('' LINEAR PHYSICS NOT ACTIVATED '')')
!  WRITE(UNIT=KULOUT,FMT='('' ============================ '')')
!  WRITE(UNIT=KULOUT,FMT='('' LENCLD2 = '',L4)')LENCLD2
!  WRITE(UNIT=KULOUT,FMT='('' LEVAPLS2= '',L4)')LEVAPLS2
!ENDIF

WRITE(unit=kulout,fmt='('' LSPHLC = '',L5,'' LVDFLC = '',L5 &
 & ,'' LVDFDS = '',L5,'' LSDRLC = '',L5,'' LSDRDS = '',L5)')&
 & lsphlc,lvdflc,lvdfds,lsdrlc,lsdrds  
WRITE(unit=kulout,fmt='('' LCZDEB = '',L5,'' LZMCON = '',L5 &
 & ,'' LKEXP  = '',L5)') lczdeb,lzmcon,lkexp  
WRITE(unit=kulout,fmt='('' ALPHA  = '',E16.6)') alpha
WRITE(unit=kulout,fmt='('' AH0    = '',E16.6,''m'')') ah0
WRITE(unit=kulout,fmt='('' USTARL = '',E16.6,''m s-1'')') ustarl
WRITE(unit=kulout,fmt='('' USTARS = '',E16.6,''m s-1'')') ustars
WRITE(unit=kulout,fmt='('' ALANDZ0= '',E16.6,''m'')') alandz0
WRITE(unit=kulout,fmt='('' ASEAZ0 = '',E16.6,''m'')') aseaz0

IF (lagphy) THEN
  WRITE(unit=kulout,fmt='('' LAGPHY = '',L5,'' THE PHYSICS PAC&
   & KAGE IS CALLED AFTER THE DYNAMICS IN LAGGED MODE'')')lagphy  
ELSE
  WRITE(unit=kulout,fmt='('' LAGPHY = '',L5,'' THE PHYSICS PAC&
   & KAGE IS CALLED BEFORE THE DYNAMICS '')')lagphy  
ENDIF
WRITE(unit=kulout,fmt=*) 'CGMIXLEN=',trim(cgmixlen)
WRITE(unit=kulout,fmt='('' LPRGML = '',L4)') lprgml

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

IF (lhook) CALL dr_hook('SU0PHY',1,zhook_handle)
END SUBROUTINE su0phy