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"

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

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

!*       1.    Set default values.

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

!*       1.1.1 Set default values for Meteo-France physics

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

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

! ALARO-0 (cloud)

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

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

! ALARO-0 (prognostic convection)

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

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

! Cloud and precipitation prognostic scheme (Lopez) :

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

!AROME physics

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

! Module YOMPHY :

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

!*    1.1.2  Set default values for simplified physical parametrization

!            of Meteo-France

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

!*    1.2.1  Set default values for ECMWF physics package

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

LE4ALB=.FALSE.

!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

!     NSTPW=30

!  Setup for 3 degree resolution wave model

!     RNORTW= 72.0

!     RSOUTW=-63.0

!     RDEGREW=3.0

!  Setup for 1.5 degree resolution wave model

!        1.3 Modify default values according to LECMWF

  LEOCWA=.FALSE.

  LEOCCO=.FALSE.

  LECOND=.FALSE.

  LEEVAP=.TRUE.

  LESICE=.TRUE.

  ENDIF

!   LERADI2 = .TRUE.

!   LERADS2 = .TRUE.

!   LEVDIF2 = .TRUE.

!   LEGWDG2 = .TRUE.

!   LECUMF2 = .TRUE.

!   LECOND2 = .TRUE.

!    IF(L131TL) THEN

!      LEPHYS=.TRUE.

!      LECOND=.TRUE.

!      LEPCLD=.FALSE.

!    ENDIF

    LSDRLC=.TRUE.

    LZMCON=.TRUE.

    LKEXP =.TRUE.

  ENDIF

!   LETRAJP = .TRUE.

!   LERADI2 = .TRUE.

!   LERADS2 = .TRUE.

!   LEVDIF2 = .TRUE.

!   LEGWDG2 = .TRUE.

!   LECUMF2 = .TRUE.

!   LECOND2 = .TRUE.

  ENDIF

!      Full physics when computing the trajectory

!      Simple scheme for TL and ADJ

    LVDFDS  = .TRUE.

    LCZDEB  = .FALSE.

    ALPHA   = 3._JPRB

    AH0     = 1000.0_JPRB

    USTARL  = 0.5_JPRB

    USTARS  = 0.2_JPRB

    ALANDZ0 = 0.05_JPRB

    ASEAZ0  = 0.0005_JPRB

  ENDIF

ELSE

  IF ((NCONF == 1).OR.(NCONF == 131).OR.(NCONF == 401)&

     & .OR.(NCONF == 501).OR.(NCONF == 601)&

    LSRCONT=.FALSE.

    LNEBT=.FALSE.

    LNEBNXR=.FALSE.

    LNEIGE=.TRUE.

      !extinction of MF phsic's keys for arome's run

      !lfpos is added for 927 confs for example.

      LGWDC=.FALSE.

      LNEBCO=.FALSE.

      LNEBR=.FALSE.

      LNEBT=.FALSE.

      LOZONE=.FALSE.

      LREWS=.FALSE.

      LRAYFM=.FALSE.

      LRAYFM15=.FALSE.

      LSNV=.FALSE.

      LSTRAS=.FALSE.

      !initialisation of AROME's physic ones

      LMICRO=.FALSE.

      LTURB=.FALSE.

      LMSE=.FALSE.

      LKFBCONV=.FALSE.

      LUSECHEM=.FALSE.

      LKFBD=.FALSE.

      LKFBS=.FALSE.

      LUSECHEM=.FALSE.

      LORILAM=.FALSE.

      LRDUST=.FALSE.

    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

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

ENDIF

!IF(.NOT.(LEPHYS.OR.LMPHYS)) THEN

!  LAGPHY=.FALSE.

!  LSLPHY=.FALSE.

!  WRITE(NULOUT,'('' INFO - LSLPHY RESET TO .FALSE. '')')

!ENDIF

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 .

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

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

!*       3.    Do tests of consistency.

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

! * Test that adequate physics is activated when at least one of the

!   options LRETCFOU, LWRTCFOU, LRCOEF,is set to .TRUE.

  ENDIF

ENDIF

ENDIF

  ENDIF

ENDIF

! * Test that options LRETCFOU, LWRTCFOU are set to .FALSE.

!   when at least LRCOEF=.TRUE.

ENDIF

! * Test that option LSRCONT is set to .FALSE. when LSRCON=.FALSE.

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)

ENDIF

! TESTS FOR ALARO-0 (prognostic convection)

  ENDIF

ENDIF

  ! * For the time being, LPHSPSH requires L3MT.OR.LSTRAPRO (missing code

  !   in CPTEND in the other cases, to be coded later).

ENDIF

! * Use simplified physics with LMPA=T?

  ! * Simplified physics is not consistent with AROME physics

  !   concerning most of the physical parameterizations.

ENDIF

ENDIF

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

!*       4.    Print final values.

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

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

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 &

 & )')&

WRITE(UNIT=KULOUT,FMT='('' NBITER = '',I2,'' NDPSFI = '',I2,'' NPHYREP = '',I2)') &

!WRITE(UNIT=KULOUT,FMT='('' NPHYINT = '',I2)')NPHYINT

WRITE(UNIT=KULOUT,FMT='(&

 & '' LECT    = '',L5,'' LAJUCV  = '',L5,/, &

 & '' LCVRAV3 = '',L5, &

 & '' LNEBGR  = '',L5,'' LNEBGY  = '',L5,'' LBCCOND = '',L5 &

WRITE(UNIT=KULOUT,FMT='(  '' LPROCLD = '',L5 &

 &,'' LEVAPP  = '',L5,'' LCOLLEC = '',L5,'' LNSMLIS  = '',L5 &

 &,'' LADJCLD  = '',L5,'' LAUTONEB  = '',L5,'' LSSD = '',L5)')&

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 &

WRITE(UNIT=KULOUT,&

  &FMT='('' LCVPRO ='',L1,'' LCDDPRO ='',L1,'' LSCMF ='',L1)') &

WRITE(UNIT=KULOUT,FMT='('' LSIMPH= '',L5,'' LTRAJPS = '',L5 &

 & ,'' LTRAJPST = '',L5 &

 & ,'' LSMOOTHD = '',L5,'' LSMOOTHA = '',L5,'' LSMOOTHB = '',L5 &

 & ,'' LCLOUDS = '',L5 )')&

WRITE(UNIT=KULOUT,FMT='('' LCVRASP = '',L5,'' LGWDSP = '',L5 &

 & ,'' LRAYSP = '',L5,'' LSTRASP = '',L5,'' LVDIFSP = '',L5 &

 & ,'' LRRMESSP = '',L5)')&

WRITE(UNIT=KULOUT,FMT='('' LRCOEF= '',L5 &

 & ,'' LTLADDIA = '',L5,'' LGLOBRAD = '',L5)')&