PREP_IDEAL_CASE

The PREP_IDEAL_CASE program prepares a Meso-NH file, that contains all the parameters and fields necessary for the execution of the Meso-NH model. Specifically, the grid parameters, the initial fields and the geophysical fields are included in this file. It is possible using this program to generate idealized fields defined by few parameters. The generated initial conditions are produced analytically, leading to quasi-1D fields or 3D fields or a single profile build with either:

layers of constant Brunt-Vaisala frequency, shear and humidity
a Radiosounding and ideal surface fields
a Radiosounding and real physiographic fields
a Radiosounding and real and ideal surface fields at the same time

For these latter cases, the initial fields may be hydrostatically or geostrophically balanced or not. For these fields to satisfy the anelastic constraint, a final correction is applied to them. The interaction between the PREP_IDEAL_CASE program and the user is made through the PRE_IDEA1.nam file. The degrees of freedom are collected in a set of namelists, read by this program.

PREP_IDEAL_CASE program and its corresponding namelist and function
Executable	Namelist	Function
PREP_IDEAL_CASE	PRE_IDEA1.nam	Prepare atmospheric initial and boundary file

Contrary to other namelist PRE_IDEA1.nam is made of two parts :

A namelist-part with directives for the preparation of an idealized case (always present). The order of namelists is free and unset namelists can be ommited.
A free-formatted part describing a vertical profile of n layers of constant moist Brunt-Vaisala frequency or a radiosounding and sometimes the explicit list of the heights of the vertical levels. This part can be present or absent in the other cases.

To initialize a simulation with a radiosounding and real terrain conditions, it is necessary to perform the PREP_PGD program to create a Meso-NH physiographic data file. This data file contains the orography and the physiographic data fields (related to the soil scheme). It is also possible to perform a complete ideal case with ideal orography and non trivial surface conditions. The user can combine the two possibilities with flags included in the namelist NAM_REAL_PGD and initialize a simulation with a real orography and idealized homogeneous surface fields. If a PGD file is specified and if the flags in namelist NAM_REAL_PGD are set to FALSE, homogeneous values can be imposed by the user in namelists from the externalized surface facility PGD (namelists NAM_COVER and NAM_ISBA), else the PREP_PGD fields are taken into account.

The following namelists can be used in the PRE_IDEA1.nam file :

NAM_CONFIO

NAM_CONFIO content
Fortran name	Fortran type	Default value
CIO_DIR	CHARACTER(LEN=512)
LVERB_OUTLST	LOGICAL	.TRUE.
LVERB_STDOUT	LOGICAL	.FALSE.
LVERB_ALLPRC	LOGICAL	.FALSE.
NGEN_VERB	INTEGER	4
NGEN_ABORT_LEVEL	INTEGER	2
NBUD_VERB	INTEGER	4
NBUD_ABORT_LEVEL	INTEGER	2
NIO_VERB	INTEGER	4
NIO_ABORT_LEVEL	INTEGER	2
LIO_COMPRESS	LOGICAL	.TRUE.
CIO_COMPRESS_ALGO	CHARACTER(LEN=10)	‘ZSTD’
NIO_COMPRESS_LEVEL	INTEGER	4
LDIAG_REDUCE_FLOAT_PRECISION	LOGICAL	.FALSE.
LIO_ALLOW_REDUCED_PRECISION_BACKUP	LOGICAL	.FALSE.
LIO_ALLOW_NO_BACKUP	LOGICAL	.FALSE.
LIO_NO_WRITE	LOGICAL	.FALSE.
NFILE_NUM_MAX	INTEGER	999

Warning

If a file is not found in the netCDF fileformat, Meso-NH will check if it exists in the LFI format and use it if found. This could be useful if you need to mix the reading of different files with different fileformats.

CIO_DIR : directory used to write outputs, backups and diachronic files (current directory by default). It can be overridden by CBAK_DIR for backups and diachronic files and by COUT_DIR for outputs.
LVERB_OUTLST : flag to write application messages in OUTPUT_LISTINGn files (in current directory, n is for the current model)
LVERB_STDOUT : flag to write application messages on the standard output
NGEN_VERB : set the verbosity level for generic messages
- 0 : no messages
- 1 : fatal messages
- 2 : error messages (and lower values)
- 3 : warning messages (and lower values)
- 4 : info messages (and lower values)
- 5 : debug messages (and lower values)
NGEN_ABORT_LEVEL : set the minimum level of generic message to abort the application (same levels as for NGEN_VERB)
NBUD_VERB : set the verbosity level for budget messages (same levels as for NGEN_VERB)
NBUD_ABORT_LEVEL : set the minimum level of budget message to abort the application (same levels as for NGEN_VERB)
NIO_VERB : set the verbosity level for IO messages (same levels as for NGEN_VERB)
NIO_ABORT_LEVEL : set the minimum level of IO message to abort the application (same levels as for NGEN_VERB)

Warning

Not all messages use this infrastructure. Therefore, some of them are not affected by these options.

LIO_COMPRESS : enable lossless compression of data for all files. This can have a negative impact on performance. This option takes precedence over their equivalent NAM_BACKUP and NAM_OUTPUT namelists.
CIO_COMPRESS_ALGO: set the compression algorithm (only for files in netCDF format, not for LFI format). The allowed values are ‘ZSTD’ (for Zstandard compression,default value), ‘DEFLATE’ (for zlib compression) or ‘NONE’. This option takes precedence over their equivalent in NAM_BACKUP and NAM_OUTPUT namelists (only if LIO_COMPRESS=.TRUE. which is the default). If set to ‘NONE’, all compression will be disabled (that stands also for lossy compression).
LOUT_COMPRESS_LEVEL : set the compression level. The value must be in the 0 to 9 interval (0 for no compression, 9 for maximum compression). This option takes precedence over their equivalent in NAM_BACKUP and NAM_OUTPUT namelists (only if LIO_COMPRESS=.TRUE. which is the default).
LDIAG_REDUCE_FLOAT_PRECISION : force writing of floating points numbers in single precision for diagnostic files (written by the DIAG program)
LIO_ALLOW_REDUCED_PRECISION_BACKUP : flag to allow writing of backup files with a reduced precision as well as reading of reduced precision files and files written with Meso-NH compiled with a lower precision for integers or reals (ie MNH_INT=4 and MNH_REAL=4).
LIO_ALLOW_NO_BACKUP : allow to have no valid backup time (useful for some tests)
LIO_NO_WRITE : disable file writes (useful for benchs)
NFILE_NUM_MAX : maximum number for numbered files (mainly backup and output files). If less than 1000, the numbers will be on 3 digits. From 1000, they will be on the number of digits of NFILE_NUM_MAX (5 if NFILE_NUM_MAX=12345).

NAM_AERO_PRE

NAM_AERO_PRE content
Fortran name	Fortran type	Default value
LORILAM	LOGICAL	.FALSE.
LDUST	LOGICAL	.FALSE.
LSALT	LOGICAL	.FALSE.
LINITPM	LOGICAL	.FALSE.
XINIRADIUSI	REAL	0.05
XINIRADIUSJ	REAL	0.2
XINISIGI	REAL	1.8
XINISIGJ	REAL	2.0
XN0IMIN	REAL	10.0
XN0JMIN	REAL	1.0
CRGUNIT	CHARACTER	‘NUMB’
NMODE_DST	INTEGER	3
XN0MIN	REAL	1.e3, 1.e1, 1.e-2
XINIRADIUS	REAL	0.044, 0.3215, 1.575
XINISIG	REAL	2.0, 1.78, 1.85
CRGUNITD	CHARACTER	‘NUMB’
NMODE_SLT	INTEGER	3
XN0MIN_SLT	REAL	1.e4, 1.e2, 1.e-1
XINIRADIUS_SLT	REAL	0.14, 1.125, 7.64
XINISIG_SLT	REAL	1.9, 2., 2.
CRGUNITS	CHARACTER	‘MASS’
CCLOUD	CHARACTER	‘NONE’
NMOD_CCN	INTEGER	1
NMOD_IFN	INTEGER	1
LDSTCAMS	LOGICAL	.FALSE.
LSLTCAMS	LOGICAL	.FALSE.

LORILAM : flag to activate chemical aerosol initialization (only if LCH_INIT_FIELD=T in NAM_CH_MNHCN_PRE).
LDUST : flag to activate passive dust initialization.
LSALT : flag to activate passive sea salt initialization.
LINITPM : flag to activate primary aerosol initialization (Black and Organic carbon) from concentration of CO (only if LORILAM=T in NAM_CH_MNHCN_PRE).
XINIRADIUSI : initial mean radius of aitken mode in $\mu m$ (only if LORILAM=T in NAM_AERO_PRE).
XINIRADIUSJ : initial mean radius of accumulation mode in $\mu m$ (only if LORILAM=T in NAM_AERO_PRE).
XINISIGI : initial standard deviation of aitken mode (only if LORILAM=T in NAM_AERO_PRE).
XINISIGJ : initial standard deviation of accumulation mode (only if LORILAM=T in NAM_AERO_PRE).
XN0IMIN : minimum number concentration of aitken mode (only if LORILAM=T in NAM_AERO_PRE).
XN0JMIN : minimum number concentration of accumulation mode (only if LORILAM=T in NAM_AERO_PRE).
CRGUNIT : definition of XINIRADIUSI and XINIRADIUSJ : mean radius can be in mass (‘MASS’) or in number (‘NUMB’) (only if LORILAM=T in NAM_AERO_PRE).
NMODE_DST : number of DUST mode (between 1 and 3 and only if LDUST=T in NAM_AERO_PRE).
XN0MIN : minimum number concentration of the NMODE_DST in particles by $m^3$ (only if LDUST=T in NAM_AERO_PRE).
XINIRADIUS : initial mean radius of the NMODE_DST modes in $\mu m$ (only if LDUST=T in NAM_AERO_PRE).
XINISIG : initial standard deviation of the NMODE_DST modes (only if LDUST=T in NAM_AERO_PRE).
CRGUNITD : definition of XINIRADIUS : mean radius can be in mass (‘MASS’) or in number (‘NUMB’) (only if LDUST=T in NAM_AERO_PRE).
NMODE_SLT : number of SALT mode in $\mu m$ (between 1 and 3 and only if LSALT=T in NAM_AERO_PRE).
XN0MIN_SLT : minimum number concentration of the NMODE_SLT in particles by $m^3$ (only if LSALT=T in NAM_AERO_PRE).
XINIRADIUS_SLT : initial mean radius of the NMODE_SLT modes (only if LSALT=T in NAM_AERO_PRE).
XINISIG_SLT : initial standard deviation of the NMODE_SLT modes (only if LSALT=T in NAM_AERO_PRE).
CRGUNITS : definition of XINIRADIUS_SLT : mean radius can be in mass (‘MASS’) or in number (‘NUMB’) (only if LSALT=T in NAM_AERO_PRE).
CCLOUD : microphysics scheme (only ‘LIMA’ possible) to use with aerosols coupling
NMOD_CCN : number of CCN modes
NMOD_IFN : number of IFN modes
LDSTCAMS : flag to activate initialization of dust aerosols from CAMS file.
LSLTCAMS : flag to activate initialization of sea-salt aerosols from CAMS file.

NAM_BLANKn

NAM_BLANKn content
Fortran name	Fortran type	Default value
XDUMMY1	REAL	0.0
XDUMMY2	REAL	0.0
XDUMMY3	REAL	0.0
XDUMMY4	REAL	0.0
XDUMMY5	REAL	0.0
XDUMMY6	REAL	0.0
XDUMMY7	REAL	0.0
XDUMMY8	REAL	0.0
NDUMMY1	INTEGER	0
NDUMMY2	INTEGER	0
NDUMMY3	INTEGER	0
NDUMMY4	INTEGER	0
NDUMMY5	INTEGER	0
NDUMMY6	INTEGER	0
NDUMMY7	INTEGER	0
NDUMMY8	INTEGER	0
LDUMMY1	LOGICAL	TRUE
LDUMMY2	LOGICAL	TRUE
LDUMMY3	LOGICAL	TRUE
LDUMMY4	LOGICAL	TRUE
LDUMMY5	LOGICAL	TRUE
LDUMMY6	LOGICAL	TRUE
LDUMMY7	LOGICAL	TRUE
LDUMMY8	LOGICAL	TRUE
CDUMMY1	CHARACTER(LEN=80)
CDUMMY2	CHARACTER(LEN=80)
CDUMMY3	CHARACTER(LEN=80)
CDUMMY4	CHARACTER(LEN=80)
CDUMMY5	CHARACTER(LEN=80)
CDUMMY6	CHARACTER(LEN=80)
CDUMMY7	CHARACTER(LEN=80)
CDUMMY8	CHARACTER(LEN=80)

Eight dummy variables and arrays (real, integer, logical, and character of length 80) are defined for testing and debugging. They are read through the namelist but are not used by any Meso-NH routine. If a developer wants to temporarily add a parameter to a subroutine, they can include a USE MODD_BLANK_n statement in that subroutine. This allows them to access and modify these variables via the namelist input.

NAM_CH_MNHCN_PRE

NAM_CH_MNHCN_PRE content
Fortran name	Fortran type	Default value
LCH_INIT_FIELD	LOGICAL	.FALSE.
CCHEM_INPUT_FILE	CHARACTER(LEN=128)	MNHC.input

LCH_INIT_FIELD : Flag to activate initialization subroutine CH_INIT_FIELD.
CCHEM_INPUT_FILE : Name of the general purpose input file for initialization.

NAM_CONF_PRE

Warning

This namelist comes from SURFEX 9.0.0 user guide https://www.umr-cnrm.fr/surfex/IMG/pdf/surfex_tecdoc.pdf.

NAM_CONF_PRE content
Fortran name	Fortran type	Default value
LCARTESIAN	LOGICAL	.TRUE.
LPACK	LOGICAL	.TRUE.
CEQNSYS	CHARACTER(LEN=3)	‘DUR’
NVERB	INTEGER	5
CIDEAL	CHARACTER(LEN=4)	‘CSTN’
CZS	CHARACTER(LEN=4)	‘FLAT’
LBOUSS	LOGICAL	.FALSE.
LPERTURB	LOGICAL	.FALSE.
LFORCING	LOGICAL	.FALSE.
LSHIFT	LOGICAL	.FALSE.
L2D_ADV_FRC	LOGICAL	.FALSE.
L2D_REL_FRC	LOGICAL	.FALSE.
NHALO	INTEGER	1
JPHEXT	INTEGER	1
LOCEAN	LOGICAL	.FALSE.

LCARTESIAN : Flag for cartesian geometry
- .TRUE. for cartesian geometry
- .FALSE. for conformal projection
LPACK : Flag to compact 1D or 2D fields
CEQNSYS : Equation system resolved by Meso-NH
- ‘LHE’ Lipps and HEmler anelastic system
- ‘DUR’ approximated form of the DURran version of the anelastic sytem
- ‘MAE’ classical Modified Anelastic Equations but with not any approximation in the momentum equation
NVERB : verbosity level
- 0 for minimum of prints
- 5 for intermediate level of prints
- 10 for maximum of prints.
Note

If CSURF=’EXTE’ in namelist NAM_GRn_PRE, NVERB=10 prints two LaTeX files containing the initialisation of surface scheme variables for each type of surface cover.
CIDEAL : kind of idealized fields
- ‘CSTN’ : Constant moist Brunt Vaisala frequency case
- ‘RSOU’ : radiosounding case
CZS : orography selector. The formulae are given below in the description of the namelist NAM_GRIDH_PRE.
- ‘FLAT’ : constant XHMAX orography (zero by default)
- ‘SINE’ : sine-shaped orography
- ‘BELL’ : bell-shaped orography
- ‘AGNE’ : orography with $h*a^2/(x^2+a^2)$ shape
- ‘DATA’ : discretized orography. The data describing the orography are given in the free format part (Free-format part).
LBOUSS : Flag for a Boussinesq version.
- .TRUE. The reference anelastic state is $\theta _{ref} = cte = \theta _{ref} (z=0)$ and $\rho _{ref} = cte = \rho _{ref} (z=0)$. In this case, the stratification is taken into account in the Meso-NH model in the flottability term. The typical length, on which this stratification varies, is much greater than the domain heigth and the $\theta_{ref}$ variation can be therefore neglected.
- .FALSE. The reference anelastic state varies with the altitude.
LPERTURB : Flag to add a perturbation on the initially horizontally homogeneous fields. This perturbation is not balanced. 3 perturbation types are implemented in the routine set_perturb.f90:
- a spherical perturbation on the dry potential temperature and the moisture fields (typical for convection initialization).
- a perturbation on the horizontal components of the wind derived from a streamfunction (typical for large scale studies). This prevents the wind from becoming divergent.
- a perturbation on the dry potential temperature field at the first mass level near the ground, corresponding to a white noise (uniform amplitude in the spectral space) (typical for Large Eddy Simulations initialization)
Note

When set to .TRUE., the parameters for the exact definition of the perturbation can be set in the namelist NAM_PERT_PRE or sometimes can be modified directly in the subroutine set_perturb.f90.
LFORCING : Flag to specify forcing sources. When .TRUE., the precise definition of the forcing is set in the free-format part of PRE_IDEA1.nam. LFORCING must be then set to .TRUE. in EXSEG1.nam (NAM_CONF).
LSHIFT : flag to shift altitudes in boundary layer. If LGEOSBAL=TRUE, LSHIFT will be set to FALSE.
L2D_ADV_FRC : flag to activate advecting forcing (2D simulation only). When .TRUE., the precise definition of the advecting forcing is set in the free-format part of PRE_IDEA1.nam.
L2D_REL_FRC : flag to activate relaxation forcing (2D simulation only). When .TRUE., the precise definition of the relaxation forcing is set in the free-format part of PRE_IDEA1.nam.
NHALO : Size of the halo for parallel distribution. This variable is related to computer performance but has no impact on simulation results.
JPHEXT : Horizontal External points number JPHEXT must be equal to 3 for cyclic cases with WENO5.
LOCEAN : flag to activate the Ocean version of Meso-NH. Pronostic variables are: Current (U and V), Vertical velocity (W), Temperature (TH), Subgrid Turbulent Kinetic Energy (TKE). Salinity (RV) can be activated with LUSERV=T. The Z-axis is directed upward (as in the atmosphere version), i.e. top of model domain corresponds to the sea surface. The initial profile must be defined in the Free-format part.

NAM_CONFn

NAM_CONFn content
Fortran name	Fortran type	Default value
LUSERV	LOGICAL	.TRUE.
LUSECI	LOGICAL	.FALSE.
LUSERC	LOGICAL	.FALSE.
LUSERR	LOGICAL	.FALSE.
LUSERI	LOGICAL	.FALSE.
LUSERS	LOGICAL	.FALSE.
LUSERG	LOGICAL	.FALSE.
LUSERH	LOGICAL	.FALSE.
NSV_USER	INTEGER	0

LUSERV : Flag to use vapor mixing ratio ($r_v$)
LUSECI : Flag to use pristine ice ($C_i$)
LUSERC : Flag to use cloud mixing ratio ($r_c$)
LUSERR : Flag to use rain mixing ratio ($r_r$
LUSERI : Flag to use ice mixing ratio ($r_i$)
LUSERS : Flag to use snow mixing ratio ($r_s$)
LUSERG : Flag to use graupel mixing ratio ($r_g$)
LUSERH : Flag to use hail mixing ratio ($r_h$)

Note

You don’t need to fill this records : they are directly managed by CCLOUD.

NSV_USER : Number of user passive scalar variables

Warning

Scalar variables needed for the 2-moment microphysical schemes, lagrangian trajectory, passive pollutants or the chemistry options are treated automatically by the model and should not be counted here.

NAM_CONFZ

NAM_CONFZ content
Fortran name	Fortran type	Default value
NZ_VERB	INTEGER	0
NZ_PROC	INTEGER	0
NB_PROCIO_R	INTEGER	1
NB_PROCIO_W	INTEGER	1
MPI_BUFFER_SIZE	INTEGER	40
LMNH_MPI_BSEND	LOGICAL	.TRUE.
LMNH_MPI_ALLTOALLV_REMAP	LOGICAL	.FALSE.
NZ_SPLITTING	INTEGER	10
NPMAX_T1DFLAT_R	INTEGER	130

NZ_VERB: level of message for NZ solver and I/O
NZ_PROC: number of processes to use in the Z splitting. The default value (0) yields an automatic calculation of the number.
NB_PROCIO_R: number of processes to use for parallel I/O when reading file. The default value (1) yields a reading from 1 file only. If more than 1 file, the 3D field are written as several 2D slices.
NB_PROCIO_W: Number of processes to use for parallel I/O when writing file. The default value (1) yields a writing into 1 file only. If more than 1 file, the 3D field are written as several 2D slices.
MPI_BUFFER_SIZE: default size for MPI_BSEND buffer in $10^6$ bytes. MPI_BUFFER_SIZE corresponds approximately to the size of the domain, that is, $NX*NY*NZ$ for I/O in 1 file, and $NX*NY$ for I/O in N 2D-slice files.
LMNH_MPI_BSEND: during HALO exchange and FFT transposition, switch to use bufferized either MPI_BSEND routine or asynchrone MPI_ISEND routine. Depending on the computer and size of the problem, one or the other option could run faster. MPI_ISEND also uses less memory so MPI_BUFFER_SIZE should be decreased.
LMNH_MPI_ALLTOALLV_REMAP:
- FALSE: FFT remap with send/recv <=> NZ_SPLITTING=10
- TRUE: FFT remap with mpi_alltoallv <=> NZ_SPLITTING=14 (BG/MPICH optimization)
NZ_SPLITTING: setting by namelist for debugging by expert user only. The non-expert user will use LMNH_MPI_ALLTOALLV_REMAP=T/F only: IZ=1=flat_inv; IZ=2=flat_invz; IZ=1+2=the two; +8=P1/P2.
NPMAX_T1DFLAT_R: setting to determine the size of the memory buffer allocated for the GPU version of Meso-NH to store real fields. The total buffer size is this setting multiplied by the number of mesh points. This buffer is used to remove allocations and deallocations that are very costly on GPUs. This value can be increased when needed (if too small, an error will be raised at runtime), but it should not be too large to avoid wasting memory.

NAM_DIMn_PRE

NAM_DIMn_PRE content
Fortran name	Fortran type	Default value
NIMAX	INTEGER	10
NJMAX	INTEGER	10

NIMAX : number of mass points in x-direction of the initial file is NIMAX+2*JPHEXT (JPHEXT corresponds to the number of marginal points in the horizontal directions and is fixed to 1 for the present Meso-NH version ). NIMAX must be equal to $2^m \times 3^n \times 5^p$ with $(m,n,p) \ge 0$.
NJMAX : number of mass points in y-direction of the physical domain. The total size of the array written in the initial file is NJMAX+2*JPHEXT. NJMAX must be equal to $2^m \times 3^n \times 5^p$ with $(m,n,p) \ge 0$.

NAM_DYNn_PRE

NAM_DYNn_PRE content
Fortran name	Fortran type	Default value
CPRESOPT	CHARACTER(LEN=5)	‘CRESI’
NITR	INTEGER	4
XRELAX	REAL	1.0
LRES	LOGICAL	.FALSE.
XRES	REAL	1.E-07

CPRESOPT : gives the type of pressure solver used for the elliptic equation (‘RICHA’, ‘CGRAD’, ‘CRESI’, ‘ZRESI’). This equation is solved in order to ensure the anelastic constraint for the initial wind field. Note that the solver is applied even for the ‘FLAT’ case when the Earth spericity is taken into account.
NITR : number of iterations used for the resolution of the elliptic equation (solver = “CPRESOPT”).
XRELAX : relaxation factor used by the Richardson method (CPRESOPT = “RICHA”).
LRES : flag to change the residual divergence limit
XRES : Value of the residual divergence limit

NAM_GRID_PRE

NAM_GRID_PRE content
Fortran name	Fortran type	Default value
XLON0	REAL	0.0
XLAT0	REAL	60.0
XBETA	REAL	0.0
XRPK	REAL	1.0
XLONORI	REAL	350.0
XLATORI	REAL	37.0

XLON0 : reference longitude for conformal projection and cartesian plane (if LCARTESIAN =.TRUE. this value can be usefull to compute local solar time)
XLAT0 : reference latitude for conformal projection and cartesian plane
XBETA : rotation angle for conformal projection and cartesian plane
XRPK : cone factor for the projection (only if LCARTESIAN =.FALSE.):
- XRPK=1: polar stereographic projection from south pole
- 1>XRPK>0: Lambert projection from south pole
- XRPK=0: Mercator projection from earth center
- -1<XRPK<0: Lambert projection from north pole
- XRPK=-1: polar stereographic projection from north pole
XLONORI : Longitude (in degrees) of the origine point (not used if LCARTESIAN =.TRUE.). This point is the mass point of conformal coordinates (x=0,y=0) of the Meso-NH grids.
XLATORI : Latitude (in degrees) of the origine point (not used if LCARTESIAN =.TRUE.)

NAM_GRIDH_PRE

NAM_GRIDH_PRE content
Fortran name	Fortran type	Default value
XLATCEN	REAL	XUNDEF
XLONCEN	REAL	XUNDEF
XDELTAX	REAL	5000.0
XDELTAY	REAL	5000.0
XHMAX	REAL	300.0 / 0.0
NEXPX	INTEGER	3
NEXPY	INTEGER	1
XAX	REAL	10000.0
XAY	REAL	10000.0
NIZS	INTEGER	5
NJZS	INTEGER	5

XLATCEN : latitude of the center of the domain for initialization. This point is vertical vorticity point (Grid).
XLONCEN : longitude of the center of the domain for initialization. This point is vertical vorticity point (Grid).
XDELTAX : mesh length (in meters) in x-direction on the conformal or cartesian plane. It is not used if you read information in a Meso-NH constant file (PGD_FILE).
XDELTAY : mesh length (in meters) in y-direction on the conformal or cartesian plane. It is not used if you read information in a Meso-NH constant file (PGD_FILE).
XHMAX : Maximum height (in meters) $h_{max}$ for orography (case CZS $\neq$ ‘FLAT’) or ground level for flat orography. Default is 0.0 for ‘FLAT’ and 300.0 elsewhere.
NEXPX : Exponent $exp_{x}$ for orography in case of CZS=’SINE’.
NEXPY : Exponent $exp_{y}$ for orography in case of CZS=’SINE’.
XAX : Widths (in meters) $a_{x}$ along x for orography in case CZS=’BELL’.
XAY : Width (in meters) $a_{y}$ along y for orography in case CZS=’BELL’.
NIZS : Localization in x-direction in the physical domain of the mountain center in the case CZS =’BELL’. ($x_{s} = NIZS * XDELTAX$). It refers to a vertical velocity point at the ground ( NIZS, NJZS )(Grid).
NJZS : Localization in y-direction in the physical domain of the mountain center in the case CZS =’BELL’. ($y_{s} = NJZS * XDELTAY$).

Note

In case of CZS = ‘BELL’ :

in the three-dimensional case :

\[z_s \left( \hat{x} , \hat{y} \right) = { h_{max} \over \left( 1 + \left( { \hat{x} - NIZS * XDELTAX \over XAX } \right) ^2 + \left( { \hat{y} - NJZS * XDELTAX \over XAY } \right) ^2 \right) ^{1.5 } }\]

in the two-dimensional case :

\[z_s \left( \hat{x} \right) = { h_{max} \over 1 + \left( { \hat{x} - NIZS * XDELTAX \over XAX } \right) ^2 }\]

NAM_GRn_PRE

Warning

This namelist comes from SURFEX 9.0.0 user guide https://www.umr-cnrm.fr/surfex/IMG/pdf/surfex_tecdoc.pdf.

NAM_GRn_PRE content
Fortran name	Fortran type	Default value
CSURF	CHARACTER(LEN=4)	‘NONE’

CSURF : ground selector.
- ‘NONE’ : no surface scheme will be activated during the future MESONH simulation, we therefore do not need any surface parameters. All the namelists of the externalized surface will be ignored.
- ‘EXTE’ : the externalized surface is used. See the SURFEX documentation for more details.

NAM_IBM_LSF

NAM_IBM_LSF content
Fortran name	Fortran type	Default value
LIBM_LSF	LOGICAL	.FALSE.
CIBM_TYPE	CHARACTER(LEN=4)	‘NONE’
NIBM_SMOOTH	INTEGER	1
XIBM_SMOOTH	REAL	0.0001

LIBM_LSF : Flag to calculate LevelSet Function (the minimum distance to the obstacles) or not.
- .TRUE.: The LevelSet Function is calculated.
- .FALSE.: The LevelSet Function is not calculated.
CIBM_TYPE : The way the obstacles are described.
- ‘NONE’: No obstacles.
- ‘IDEA’: Idealised ellipsoide or parallelepiped obstacles that are described in the additional namelist ibm_idea.nam.
- ‘GENE’: Generic user defined obstacles in .obj format.
NIBM_SMOOTH : The number of iterations for smoothing the LevelSet Function. In the case a considerable smoothing shall be done, it is recommended to use NIBM_SMOOTH=10.
XIBM_SMOOTH : The characteristic length scale used for smoothing the LevelSet Function. It is recommended to use a value of XIBM_SMOOTH close to the grid size.

ibm_idea.nam

ibm_idea.nam content
Fortran name	Fortran type	Default value
NOBJ	INTEGER
NOBJ_TYPE	INTEGER
NTYPE	INTEGER
NOBJ_ETYPE	INTEGER
XX1	REAL
XX2	REAL
XY1	REAL
XY2	REAL
XZ1	REAL
XZ2	REAL

NOBJ : Total number of obstacles (parallelepiped and ellipsoide).
NOBJ_TYPE : Number of obstacles type (parallelepiped and/or ellipsoide).
NTYPE : Obstacle type:
- 1: Parallelepiped.
- 2: Ellipsoide.
NOBJ_ETYPE : Number of obstacles for each type.
XX1, XX2, XY1, XY2, XZ1, XZ2 : Locations of the obstacles.
- if type 1 (parallelepiped): I1/I2 = min/max in direction (X,Y,Z)
- if type 2 (ellipsoide): X1/Y1/Z1 locations of object center, X2/Y2/Z2 axe lenght in each direction

Note

Example (case of 4 parallelepipeds and 1 ellipsoide) :

5 2 # NOBJ / NOBJ_TYPE
1 4 # NTYPE / NOBJ_ETYPE
+51.00 +53.42 +51.45 +63.65 -1.00 +2.54 # XX1/XX2/XY1/XY2/XZ1/XZ2
+51.00 +53.42 +71.55 +83.75 -1.00 +2.54 #
+51.00 +53.42 +91.65 +103.85 -1.00 +2.54 #
+51.00 +53.42 +111.75 +123.95 -1.00 +2.54 #
2 1 # NTYPE / NOBJ_ETYPE
+100. +10. +100. +10. -1. 10. # XX1/XX2/XY1/XY2/XZ1/XZ2

In the case of obstacles in contact with the ground it is necessary to have negative value of XZ1.

ibm_gene.obj

The ibm_gene.obj file contains the information of the triangles constituting the faces of the obstacles. The .obj file must have a particular organization:

A line with ‘usemtl’ indicates the two materials of each side of the interface. Only the faces with their external face in contact with the outside air are read (mat2=air).
A line starting with ‘v’ indicates the location (x,y,z coordinates) of a triangle vertex.
A line starting with ‘f’ indicates the vertex constituting one triangle.

usemtl mat1:mat2
v xv1 yv1 zv1
v xv2 yv2 zv2
v xv3 yv3 zv3
v xv4 yv4 zv4
f 1 2 3
f 1 3 4
usemtl mat3:mat2
v xv5 yv5 zv5
v xv6 yv6 zv6
v xv7 yv7 zv7
v xv8 yv8 zv8
f 5 6 7
f 5 7 8

Note

In the above example, the first triangle is formed by vertex numbers 1, 2, 3 and the second by vertex numbers 1, 3, 4. The two triangles form the interface mat1:mat2. The interface mat3:mat2 is defined by two other triangles. The vertex number increment starts at the first vertex defined in the file. There are numerous other ways to write an .obj file. For instance, the normal or the texture of the faces can be defined. For the sake of simplicity only material (‘usemtl’), vertex coordinates (‘v’) and face definition (‘f’) are considered in Meso-NH. The face normal is computed directly in the code and the face texture is irrelevant in the present version of the code.

NAM_LBCn_PRE

NAM_LBCn_PRE content
Fortran name	Fortran type	Default value
CLBCX	ARRAY(2 CHARACTERS(LEN=4))	2*’CYCL’
CLBCY	ARRAY(2 CHARACTERS(LEN=4))	2*’CYCL’

CLBCX : represent the type of lateral boundary condition at the left and right boundaries along x (CLBCX(1) and CLBCX(2) respectively). Possible values are :
- ‘CYCL’ : for cyclic boundary condition
- ‘OPEN’ : for open boundary condition
- ‘WALL’ : for rigid wall boundary condition

Note

It should be note that CLBCX(1) or CLBCY(1) refers to the lowest index values (IIB, IJB for X and Y directions) and CLBCX(2) or CLBCY(2) to the highest index values (IIE and IJE).
The same boundary conditions must be used for the MESONH run itself (NAM_LBCn)
Note also that CYCLIC conditions are not possible with a PGD file (CPGD_FILE different to ‘ ‘ in NAM_REAL_PGD).

CLBCY : same as CLBCX but for the left and right boundaries along y (CLBCY(1) and CLBCY(2) respectively). They are strings of 4 characters.

NAM_LUNITn

NAM_LUNITn content
Fortran name	Fortran type	Default value
CINIFILE	CHARACTER(LEN=128)	‘INIFILE’
CINIFILEPGD	CHARACTER(LEN=128)	‘ ‘
CCPLFILE	CHARACTER(LEN=128)(:)	NONE

Warning

This namelist is shared in by PREP_IDEAL_CASE and MESONH programs but the CCPLFILE option is only relevant for MESONH.

CINIFILE : name of the initial Meso-NH file produced by PREP_IDEAL_CASE, it will then be used as initial file in a MESONH simulation.
CINIFILEPGD : name of the PGD file if CSURF $\neq$ ‘NONE’ :
- If you use an input PGD file for the step PREP_IDEAL_CASE (CPGD_FILE in NAM_REAL_PGD), you must have CINIFILEPGD=CPGD_FILE.
- If there is no input PGD, CINIFILEPGD is the name of the PGD file produced by PREP_IDEAL_CASE.
CCPLFILE : name of the files which contains the field values used for the coupling of the outermost MESONH model. No more than JPCPLFILEMAX=1000 (since MNH-V6-0-0) files can be used in a simulation. These CCPLFILE file names are only meaningful for the outermost model which finds its boundary conditions from a previously executed run of Meso-NH or another model (prepared by PREP_REAL_CASE). No constraint are imposed on the coupling file names only that they must be temporally ordered
If the coupling files are given by
CCPLFILE(1)= ’F_1’ -> t1

CCPLFILE(2)= ’F_2’ -> t2

CCPLFILE(3)= ’A_2’ -> t3

CCPLFILE(4)= ’A_5’ -> t4
then, the instants must satisfy : tsegment ≤ t1 < t2 < t3 < t4. If it is not the case, the program stops. If the coupling fields are not time dependent, no coupling files are required because the coupling fields are read from the inital MESONH file of model 1 as the Larger scale fields ( LSUM, LSVM, LSWM, LSTHM, LSRVM ). More details can be found in the scientific documentation of the model.

NAM_PERT_PRE

NAM_PERT_PRE content
Fortran name	Fortran type	Default value
CPERT_KIND	CHARACTER(LEN=2)	‘TH’
XAMPLITH	REAL	1.5
XAMPLIRV	REAL	0.0
XAMPLIUV	REAL	1.0834
XAMPLIWH	REAL	0.1
NKWH	INTEGER	2
LSET_RHU	LOGICAL	.TRUE.
XCENTERZ	REAL	2000.0
XRADX	REAL	10000.0
XRADY	REAL	10000.0
XRADZ	REAL	2000.0
LWH_LBXU	LOGICAL	.FALSE.
LWH_LBYV	LOGICAL	.FALSE.

CPERT_KIND : Defines the type of the perturbation
- ‘TH’ : thermodynamical fields perturbation ($\theta$ and $r_v$)
- ‘UV’ : horizontal wind fields perturbation (U and V)
- ‘WH’ : white noise applied to $\theta$
- ‘WW’ : white noise applied to wind components
XAMPLITH : maximum perturbation for $\theta$
XAMPLIRV : maximum perturbation for $r_v$
XAMPLIUV : maximum perturbation for U and V
XAMPLIWH : maximum perturbation for the normalized white noise (temperature or wind)
NKWH : Upper level of the layer starting from the ground where the white noise is applied
LSET_RHU : Conservation of the relative humidity
- TRUE the relative humidity is conserved in the $\theta$ perturbation
- FALSE the $r_v$ perturbation is computed with the XAMPLIRV amplitude
XCENTERZ : Height of the maximum of the $\theta$ perturbation (m)
XRADX : radius of the perturbation along X (m)
XRADY : radius of the perturbation along Y (m)
XRADZ : radius of the perturbation along Z (m)
LWH_LBXU : White noise in inflow and outflow LBC of U
LWH_LBXV : White noise in inflow and outflow LBC of V

NAM_REAL_PGD

NAM_REAL_PGD content
Fortran name	Fortran type	Default value
CPGD_FILE	CHARACTER(LEN=128)
LREAD_ZS	LOGICAL	.FALSE.
LREAD_GROUND_PARAM	LOGICAL	.FALSE.

CPGD_FILE : name of the physiographic data file containing the ground data fields. The file must be generated by the PREP_PGD program.

Note

For a purely ideal case, the CPGD_FILE variable may be deleted from the namelist or set to its default value ‘ ‘.
The horizontal grid will be read in the PGD file and therefore, the mesh increments XDELTAX and XDELTAY are no more used.

LREAD_GROUND_PARAM : Flag to use or not the surface cover types (COVERnnn) and all other physiographic fields (except orographic ones) read in the PGD file.
- .TRUE. to read the data in the PGD file
- .FALSE. to use XUNIF_COVER idealized homogeneous values given in the namelist NAM_COVER (from the externalized surface) and scratch the PGD_FILE data
LREAD_ZS : Flag to use or not the orography parameters read in the PGD file.
- .TRUE. to use the data read in the PGD_FILE
- .FALSE. to use an idealized orography given in the namelist NAM_GRIDH_PRE and scratch the PGD_FILE data

NAM_SLEVE

This namelist defines the orography file and orographic treatment to be done.

NAM_SLEVE content
Fortran name	Fortran type	Default value
NSLEVE	INTEGER	12
XSMOOTH_ZS	REAL	XUNDEF

NSLEVE : number of iteration for computation of smooth orography.
XSMOOTH_ZS : optional uniform smooth orography.

NAM_VER_GRID

There are three ways to compute the vertical grid:

constant grid mesh: only the number of levels NKMAX and the grid mesh sizes ZDZGRD and ZDZTOP are used. ZDZGRD and ZDZTOP must have the same value. The type of grid YZGRID_TYPE has to be set to ‘FUNCTN’.
two layers are defined, with different stretching in each layer: the grid mesh size is given near the ground with ZDZGRD and at top of the model with ZDZTOP and the stretching coefficients ZSTRGRD, ZSTRTOP and ZZMAX_STRGRD has to be defined. It is possible that the top grid size is never reached if the number of points is not enough for the prescribed stretchings. The type of grid YZGRID_TYPE has to be set to ‘FUNCTN’.
the vertical discretization is given by the user: the type of grid YZGRID_TYPE has to be set to ‘MANUAL’ and only the number of levels NKMAX is used.

The variables of this namelist are:

NAM_VER_GRID content
Fortran name	Fortran type	Default value
LTHINSHELL	LOGICAL	.FALSE.
NKMAX	INTEGER	10
YZGRID_TYPE	CHARACTER(LEN=6)	‘FUNCTN’
ZDZGRD	REAL	300.0
ZDZTOP	REAL	300.0
ZZMAX_STRGRD	REAL	0.0
ZSTRGRD	REAL	0.0
ZSTRTOP	REAL	0.0
LSLEVE	LOGICAL	.FALSE.
XLEN1	REAL	7500.0
XLEN2	REAL	2500.0

LTHINSHELL : Flag for the thinshell approximation
NKMAX : number of points in z-direction of the required physical domain. The total size of the array written in initial file will be NKMAX + 2*JPVEXT (JPVEXT is fixed to 1 for the present version of Meso-NH).
YZGRID_TYPE : type of vertical grid definition:
- ‘FUNCTN’: the vertical grid is given by a regular logarithmic function, whose variation is determined by the values of free parameters ZDZGRD, ZDZTOP, ZSTRGRD, ZSTRTOP and ZZMAX_STRGRD.
- ‘MANUAL’: the levels are explicitly given in the free-formatted part with the keyword ZHAT by entering the heights of the different levels from K=2 to K= KMAX + 2).
ZDZGRD : mesh length in z-direction near the ground
ZDZTOP : mesh length in z-direction near the top of the model
ZZMAX_STRGRD : Altitude separating the two constant stretching layers
ZSTRGRD : Constant imposed stretching (in %) in the lower layer (below ZZMAX_STRGRD)
ZSTRTOP : Constant imposed stretching (in %) in the upper layer (above ZZMAX_STRGRD)
LSLEVE : flag for Sleve vertical coordinate.
XLEN1 : decay scale for smooth topography (in meters)
XLEN2 : decay scale for smale-scale topography deviation (in meters)

NAM_VPROF_PRE

NAM_VPROF_PRE content
Fortran name	Fortran type	Default value
LGEOSBAL	LOGICAL	FALSE
CFUNU	CHARACTER(LEN=3)	ZZZ
CFUNV	CHARACTER(LEN=3)	ZZZ
CTYPELOC	CHARACTER(LEN=6)	IJGRID
XLATLOC	REAL	45.0
XLONLOC	REAL	0.0
XXHATLOC	REAL	20000.0
XYHATLOC	REAL	20000.0
NILOC	INTEGER	4
NJLOC	INTEGER	4

LGEOSBAL : Flag to fulfill the geostrophic balance or not
- .TRUE. the geostrophic balance is satisfied by the initial fields
- .FALSE. the geostrophic balance is not satisfied by the initial fields
CFUNU : String of 3 characters, describing the type of function, which gives the x component of the wind. Possible configurations are listed below :
- ‘ZZZ’ : U = U(z). The U(z) values are taken from the Radio-Sounding or analitycal profile given in the free-formatted part of the PRE_IDEA1.nam file.
- ‘Y*Z’ : U= F(Y)*U(Z). The U(z) values are build in the same way as the ‘ZZZ’ case and the function F(Y) is a simple function of Y, which must be adapted by modifying its definition directly in the routine FUNUY. The default function is :
\[FUNUY(\hat{y}) = {1 \over \cosh \left({ \hat{y} - \hat{y_0} \over zwidth } \right)}\]
- ‘Y,Z’ : U= G(Y,Z). The function G must also be adapted by modifying its definition directly in the routine FUNUYZ. The default function is :
\[FUNUYZ(\hat{y},z) = { 1 \over \cosh \left(\left( { \hat{y} - \hat{y_0} \over zwidthy } \right) ^2 + \left( { z - z_0 \over zwidthz } \right) ^2 \right) }\]

Note

Notice that in this case the U(z) values given by the profile are not used.
CFUNV : String of 3 characters, describing the type of function, which gives the y component of the wind. Possible configurations are listed below
- ‘ZZZ’ : V = V(z). The V(z) values are taken from the Radio-Sounding or analitycal profile given in the free-formatted part of the PRE_IDEA1.nam file.
- ‘X*Z’ : V= F(X)*V(Z). The V(z) values are build in the same way as the ‘ZZZ’ case and the function F(X) is a simple function of X, which must be adapted by modifying its definition directly in the routine FUNVX. The default function is :
\[FUNVX(\hat{x}) = { 1 \over \cosh \left({ \hat{x} - \hat{x_0} \over zwidth } \right) }\]
- ‘X,Z’ : V= G(X,Z). The function G must also be adapted by modifying its definition directly in the routine FUNVXZ. The default function is :
\[FUNVXZ(\hat{x},z) = { 1 \over \cosh \left(\left( { \hat{x} - \hat{x_0} \over zwidthx } \right) ^2 + \left( { z - z_0 \over zwidthz } \right) ^2 \right) }\]

Note

Notice that in this case the V(z) values given by the profile are not used.
CTYPELOC : Type of information used to give the localization of vertical profile :
- ‘IJGRID’ for (i,j) point on index space
- ‘XYHATM’ for (x,y) coordinates on conformal plane or cartesian plane
- ‘LATLON’ for (latitude,longitude) on spherical earth
XLATLOC : Latitude (in degrees) of the vertical profile localization (used in case CTYPELOC=’LATLON’)
XLONLOC : Longitude (in degrees) of the vertical profile localization (used in case CTYPELOC=’LATLON’)
XXHATLOC : position (in meters) x of the vertical profile localization (used in cases CTYPELOC=’XYHATM’)
XYHATLOC : position (in meters) y of the vertical profile localization (used in cases CTYPELOC=’XYHATM’)
NILOC : position i in the physical domain of the vertical profile localization (used in cases CTYPELOC=’IJGRID’). If you use a 1D model, then NILOC is reset to 1 by the program.
NJLOC : position j in the physical domain of the vertical profile localization (used in cases CTYPELOC=’IJGRID’) If you use a 1D or a 2D model, then NJLOC is reset to 1 by the program.

Note

It is important to setup NILOC and NJLOC properly especially if a topography is present in the simulation. The vertical profile is interpolated on the horizontal regarding the topography for the other grid points.

NAM_PGD_SCHEMES

Warning

This namelist comes from SURFEX 9.0.0 user guide https://www.umr-cnrm.fr/surfex/IMG/pdf/surfex_tecdoc.pdf.

This namelist defines the four schemes that will be used, one for each type of surface (sea, inland water, town, vegetation).

NAM_PGD_SCHEMES content
Fortran name	Fortran type	Default value
CNATURE	CHARACTER(LEN=6)	‘ISBA’
CSEA	CHARACTER(LEN=6)	‘SEAFLX’
CWATER	CHARACTER(LEN=6)	‘WATFLX’
CTOWN	CHARACTER(LEN=6)	‘TEB’

CNATURE : scheme used for vegegation and natural soil covers. The different possibilities are:
- “NONE”: no scheme used. No fluxes will be cmputed at the surface.
- “FLUX”: ideal fluxes are prescribed. The have to be set in the fortran routine init_ideal_flux.f90.
- “TSZ0”: In this cheme, the fluxes are computed according to the ISBA physics, but the surface characteristics (temperature, humidity, etc…) remain constant with time.
- “ISBA”: this is the full ISBA scheme [Noilhan and Planton, 1989], with all options developped since this initial paper.
CSEA : scheme used for sea and ocean. The different possibilities are:
- “NONE” : no scheme used. No fluxes will be cmputed at the surface.
- “FLUX” : ideal fluxes are prescribed. The have to be set in the fortran routine init_ideal_flux.f90.
- “SEAFLX” : this is a relatively simple scheme, using the Charnock formula.
CWATER : scheme used for inland water. The different possibilities are:
- “NONE” : no scheme used. No fluxes will be cmputed at the surface.
- “FLUX” : ideal fluxes are prescribed. The have to be set in the fortran routine init_ideal_flux.f90.
- “WATFLX” : this is a relatively simple scheme, using the Charnock formula.
- “FLAKE” : this is lake scheme from Mironov, 2005.
CTOWN : scheme used for towns. The different possibilities are:
- “NONE” : no scheme used. No fluxes will be cmputed at the surface.
- “FLUX” : ideal fluxes are prescribed. The have to be set in the fortran routine init_ideal_flux.f90.
- “TEB” : this is the Town Energy Balance scheme [Masson, 2000], with all the subsequent ameliorations of the scheme.
LGARDEN : general flag to activate TEB_GARDEN

NAM_COVER

Warning

This namelist comes from SURFEX 9.0.0 user guide https://www.umr-cnrm.fr/surfex/IMG/pdf/surfex_tecdoc.pdf.

This namelist gives the information to compute the surface cover fractions when ECOCLIMAP is used. It is possible to use an existing ECOCLIMAP map or the define the ECOCLIMAP covers for the user’s domain.

NAM_COVER content
Fortran name	Fortran type	Default value
XUNIF_COVER	REAL
YCOVER	CHARACTER(LEN=28)
YCOVERFILETYPE	CHARACTER(LEN=6)
XRM_COVER	REAL	10-6
XRM_COAST	REAL	1.0
XRM_LAKE	REAL	0.0
XRM_SEA	REAL	0.0
LORCA_GRID	LOGICAL	.FALSE.
XLAT_ANT	REAL	-77.0
LIMP_COVER	LOGICAL	.FALSE.
LRM_RIVER	LOGICAL	.FALSE.

XUNIF_COVER : specified values for uniform cover fractions. For each index i between 1 and 573, XUNIF_COVER(i) is the fraction of the ith ecosystem of ecoclimap. The same fraction of each ecosystem is set to all points of the grid. The sum of all ecosystem fractions must be equal to one. If XUNIF_COVER is set, it has priority on the use of an ecosystem file (see next item: YCOVER). In the case of grid without any reference to geographical coordinates (“CARTESIAN “ or “NONE “), XUNIF_COVER must be set.
YCOVER : ecoclimap data file name. It is used only if XUNIF_COVER is not set.
YCOVERFILETYPE : type of YCOVER file (‘DIRECT’, ‘BINLLV’, ‘BINLLF’, ‘ASCLLV’).
XRM_COVER : for each point, all fractions of ecosystems that are below XRM_COVER are removed (i.e. set to zero), and the corresponding area fractions are distributed among the remaining ecosystem fractions. Whatever the value of XRM_COVER, at least one ecosystem remains for each grid point.
XRM_COAST : limit of coast coverage under which the coast is replaced by sea or inland water.
XRM_LAKE : limit of inland lake coverage under which the water is removed.
XRM_SEA : limit of sea coverage under which the sea is removed.
LORCA_GRID : flag to ensure the compatibility between surfex and Orca grid which minimal latitude over Antarctica is 77S
XLAT_ANT : minimum Orca grid latitude over Antarctica
LIMP_COVER : reads the cover fractions in an existing PGD file to avoid their computation
LRM_RIVER : if T, rivers (cover 3) are removed

NAM_ISBA

Warning

This namelist comes from SURFEX 9.0.0 user guide https://www.umr-cnrm.fr/surfex/IMG/pdf/surfex_tecdoc.pdf.

NAM_ISBA content
Fortran name	Fortran type	Default value
NPATCH	INTEGER	1
NGROUND_LAYER	INTEGER	1E+9
CISBA	CHARACTER(LEN=3)	‘3-L’
CPEDO_FUNCTION	CHARACTER(LEN=4)	‘CH78’
CPHOTO	CHARACTER(LEN=3)	‘NON’
LTR_ML	LOGICAL	.FALSE.
CALBEDO	CHARACTER(LEN=3)	‘CM13’
XRM_PATCH	REAL	0.0
XUNIF_CLAY	REAL	0.33
YCLAY	CHARACTER(LEN=28)
YCLAYFILETYPE	CHARACTER(LEN=6)	‘ASCLLV’
XUNIF_SAND	REAL	0.33
YSAND	CHARACTER(LEN=28)
YSANDFILETYPE	CHARACTER(LEN=6)	‘ASCLLV’
XUNIF_RUNOFFB	REAL	0.5
YRUNOFFB	CHARACTER(LEN=28)
YRUNOFFBFILETYPE	CHARACTER(LEN=6)	‘ASCLLV’
XUNIF_WDRAIN	REAL	0.0
YWDRAIN	CHARACTER(LEN=28)
YWDRAINFILETYPE	CHARACTER(LEN=6)	‘ASCLLV’
YCTI	CHARACTER(LEN=28)
YCTIFILETYPE	CHARACTER(LEN=6)
XUNIF_SOC_TOP	REAL	1.E+20
XUNIF_SOC_SUB	REAL	1.E+20
YSOC_TOP	CHARACTER(LEN=28)
YSOC_SUB	CHARACTER(LEN=28)
YSOCFILETYPE	CHARACTER(LEN=6)
XUNIF_PERM	REAL	1.E+20
YPERM	CHARACTER(LEN=28)
YPERMFILETYPE	CHARACTER(LEN=28)
XUNIF_PH	REAL	1.E+20
YPH	CHARACTER(LEN=28)
YPHFILETYPE	CHARACTER(LEN=28)
XUNIF_FERT	REAL	1.E+20
YFERT	CHARACTER(LEN=28)
YFERTFILETYPE	CHARACTER(LEN=28)
LIMP_SAND	LOGICAL	.FALSE.
LIMP_CLAY	LOGICAL	.FALSE.
LIMP_CTI	LOGICAL	.FALSE.
LIMP_SOC	LOGICAL	.FALSE.
LIMP_PERM	LOGICAL	.FALSE.
XSOILGRID	REAL(150)	1.E+20
LMEB	LOGICAL	.FALSE.
LLULCC	LOGICAL	.FALSE.

NPATCH : number of patches used in ISBA. One patch corresponds to aggregated parameters. 12 patches correspond to separate energy budgets for all vegetation types present in ISBA. 3 patches correspond to bare soil types, low vegetation, trees. If CPHOTO equals ‘NON’ any number of patches between 1 and 12 is possible, for the other values of CPHOTO, 12 patches are required. The order and the signification of each patch is the following:
- 1: no vegetation (smooth) - NO
- 2: no vegetation (rocks) - ROCK
- 3: permanent snow and ice - SNOW
- 4: temperate broadleaf cold-deciduous summergreen - TEBD (TREE)
- 5: boreal needleleaf evergreen - BONE (CONI)
- 6: tropical broadleaf evergreen - EVER
- 7: C3 cultures types - C3
- 8: C4 cultures types - C4
- 9: irrigated crops - IRR
- 10: grassland (C3) - GRAS
- 11: tropical grassland (C4) - TROG
- 12: peat bogs, parks and gardens (irrigated grass) - PARK
- 13: tropical broadleaf deciduous - TRBD (TREE)
- 14: temperate broadleaf evergreen - TEBE (TREE)
- 15: temperate needleleaf evergreen - TENE (CONI)
- 16: boreal broadleaf cold-deciduous summergreen - BOBD (TREE)
- 17: boreal needleleaf cold-deciduous summergreen - BOND (CONI)
- 18: boreal grass - BOGR (GRASS)
- 19: shrub - SHRB (TREE)
NGROUND_LAYER : number of soil layer used in case of diffusion physics in the soil (CISBA = ‘DIF’):
- with CISBA = 2-L, NGROUND_LAYER default is 2
- with CISBA = 3-L, NGROUND_LAYER default is 3
- with CISBA= DIF and LECOCLIMAP, NGROUND_LAYER default is 14
CISBA : type of soil discretization and physics in ISBA:
- ‘2-L’: force-restore method with 2 layers for hydrology
- ‘3-L’: force-restore method with 3 layers for hydrology
- ‘DIF’: diffusion layer, with any number of layers
CPEDO_FUNCTION : Pedo-transfert function for DIF. The following options are currently available:
- “CH78”: Clapp and Hornberger 1978 for BC
- “C084”: Cosby et al. 1988 for BC
CPHOTO :type of photosynthesis physics. The following options are currently available:
- “NON”: none is used. Jarvis formula is used for plant transpiration.
- “AST”: ISBA-AGS with offensive/defensive stress, without evolving Leaf Area Index
- “NIT”: ISBA-AGS with nitrogen, with evolving Leaf Area Index
- “NCB”: ISBA-AGS with nitrogen, with evolving Leaf Area Index and wood, soil, roots biomass
LTR_ML : to activate new radiative transfert calculation, only if CPHOTO/=NON.
CALBEDO : type of bare soil albedo. The following options are currently available
- “DRY “: dry bare soil albedo
- “WET “: wet bare soil albedo
- “MEAN”: albedo for bare soil half wet, half dry
- “EVOL”: albedo of bare soil evolving with soil humidity
- “CM13”: albedo by cover and vegetation type processed from satellite data
XRM_PATCH : threshol to remove little fractions of patches
XUNIF_CLAY : uniform prescribed value of clay fraction.
YCLAY : clay fraction data file name.
YCLAYFILETYPE : type of clay data file (‘DIRECT’, ‘BINLLF’, ‘BINLLV’, ‘ASCLLV’)
XUNIF_SAND : uniform prescribed value of sand fraction.
YSAND : sand fraction data file name.
YSANDFILETYPE : type of sand data file (‘DIRECT’, ‘BINLLF’, ‘BINLLV’, ‘ASCLLV’)
XUNIF_RUNOFFB : uniform prescribed value of subgrid runoff coefficient.
YRUNOFFB : subgrid runoff coefficient data file name.
YRUNOFFBFILETYPE : type of subgrid runoff data file (‘DIRECT’, ‘BINLLF’, ‘BINLLV’, ‘ASCLLV’)
XUNIF_WDRAIN : uniform prescribed value of subgrid drainage.
YWDRAIN : subgrid drainage data file name.
YWDRAINFILETYPE : type of subgrid drainage data file (‘DIRECT’, ‘BINLLF’, ‘BINLLV’, ‘ASCLLV’)
YCTI : topographic indices file name.
YCTIFILETYPE : type of topographic file (‘DIRECT’, ‘BINLLF’, ‘BINLLV’, ‘ASCLLV’)
XUNIF_SOC_TOP : uniform prescribed value of topsoil organic carbon (used only in CSOC=SGH in NAM_ISBAn)
XUNIF_SOC_SUB : uniform prescribed value of subsoil organic carbon (used only in CSOC=SGH in NAM_ISBAn)
YSOC_TOP : organic carbon topsoil data file name (used only in CSOC=SGH in NAM_ISBAn).
YSOC_SUB : organic carbon subsoil data file name (used only in CSOC=SGH in NAM_ISBAn).
YSOCFILETYPE : type of organic matter data file (‘DIRECT’, ‘BINLLF’, ‘BINLLV’, ‘ASCLLV’) (used only in CSOC=SGH in NAM_ISBAn)
XUNIF_PERM : uniform value of permafrost distribution (used only if CISBA=DIF)
YPERM : file name for permafrost distribution (used only if CISBA=DIF)
YPERMFILETYPE : permafrost distribution data file type(‘DIRECT’, ‘BINLLF’, ‘BINLLV’, ‘ASCLLV’) (used only if CISBA=DIF)
XUNIF_PH : uniform value of soil pH (used only if LCH_NO_FLUX=T)
YPH : file name for soil pH (used only if LCH_NO_FLUX=T)
YPHFILETYPE : soil pH data file type (‘DIRECT’, ‘BINLLF’, ‘BINLLV’, ‘ASCLLV’) (used only if LCH_NO_FLUX=T)
XUNIF_FERT : uniform value of soil fertilization rate (kgN/ha/h) (used only if LCH_NO_FLUX=T)
YFERT : file name for soil fertilisation rate (kgN/ha/h) (used only if LCH_NO_FLUX=T)
YFERTFILETYPE : soil fertilisation rate file type (kgN/ha/h)(‘DIRECT’, ‘BINLLF’, ‘BINLLV’, ‘ASCLLV’) (used only if LCH_NO_FLUX=T)
LIMP_SAND : reads sand fraction in an existing PGD file
LIMP_CLAY : reads clay fraction in an existing PGD file
LIMP_CTI : reads topographic indices in an existing PGD file
LIMP_SOC : reads organic carbon in an existing PGD file
LIMP_PERM : reads permafrost distribution in an existing PGD file
XSOILGRID : uniform soil depth grid for CISBA=DIF. Default with CISBA=DIF and LECOCLIMAP is (/0.01,0.04,0.10,0.20,0.40,0.60,0.80,1.00,1.50,2.00,3.00,5.00,8.00,12.0/)
LMEB : Flag to activate MEB (please note that by default, MEB uses the TR_LM radiation scheme, so when LMEB=T, LTR_ML=T automatically)
LLULCC : land-use land cover change scheme activation key

NAM_CH_EMIS_PGD

This namelist is used to initialize chemistry components emissions. You can treat up to 999 such fields. These fields will be written on all the files you will use later (after prognostic fields initialization, or during and after run, etc…). Their name in the files are ‘EMIS_GRnnn’, where nnn goes from 001 to 999. During the execution of the programs, these fields are stored in the XEMIS_GR_FIELDS(:,:) (first dimension: spatial dimension, second dimension: total number of fields), in the module MODD_EMIS_GR_FIELDn. The temporal evolution, the aggregation of prescribed emissions and the link with the corresponding chemical prognostic variables are handled by the subroutine CH_EMISSION_FLUXn.f90.

NAM_CH_EMIS_PGD content
Fortran name	Fortran type	Default value
NEMIS_PGD_NBR	INTEGER	0
CEMIS_PGD_NAME(:)	CHARACTER(LEN=20)
CEMIS_PGD_FILE(:)	CHARACTER(LEN=28)
CEMIS_PGD_COMMENT(:)	CHARACTER(LEN=40)
NEMIS_PGD_TIME	INTEGER	0
CEMIS_PGD_FILETYPE(:)	CHARACTER(LEN=6)	DIRECT
CEMIS_PGD_AREA(:)	CHARACTER(LEN=3)	ALL
CEMIS_PGD_ATYPE(:)	CHARACTER(LEN=3)	ARI

Only the first NEMIS_PGD_NBR values in these arrays are meaningfull.

NEMIS_PGD_NBR : number of dummy fields.
CEMIS_PGD_NAME(:) : list of the dummy fields you want to initialize with your own data. You can give any name you want. This is a way to describe what is the field. This information is not used by the program. It is just written in the Meso-NH files.
CEMIS_PGD_FILE(:) : list of the names of the files containing the data for the fields you have specified in CEMIS_PGD_NAME(:).
CEMIS_PGD_COMMENT(:) : list of the comments associated to each emission field.
NEMIS_PGD_TIME(:) : list of the time of the files containing the data for the fields you have specified in CEMIS_PGD_NAME(:).
CEMIS_PGD_FILETYPE(:) : list of the types of the files containing the data for the fields you have specified in CEMIS_PGD_NAME(:) (‘DIRECT’, ‘BINLLF’, ‘BINLLV’, ‘ASCLLV’).
CEMIS_PGD_AREA(:) : area of meaningfullness of the fields you have specified in CEMIS_PGD_NAME(:) (‘ALL’, ‘NAT’, ‘TWN’, ‘SEA’, ‘WAT’, ‘LAN’, respectively for everywhere, natural areas, town areas, sea, inland waters, land = natural cover + town). For example, oceanic emission of DNS is relevant on ‘SEA’.
CEMIS_PGD_ATYPE(:) : type of averaging (during PGD for the fields you have specified in CEMIS_PGD_NAME(:) (‘ARI’, ‘INV’, ‘LOG’, respectively for arithmetic, inverse and logarithmic averaging).

Note

Example:

&NAM_CH_EMIS_PGD
   NEMIS_PGD_NBR = 2,
   CEMIS_PGD_NAME(1)='COE',
   NEMIS_PGD_TIME(1)=0,
   CEMIS_PGD_COMMENT(1)='CO_00h00',
   CEMIS_PGD_AREA(1)='LAN',
   CEMIS_PGD_ATYPE(1)='ARI',
   CEMIS_PGD_FILE(1)='co_00.asc',
   CEMIS_PGD_FILETYPE(1)='ASCLLV',
   CEMIS_PGD_NAME(2)='COE',
   NEMIS_PGD_TIME(2)=43200,
   CEMIS_PGD_COMMENT(2)='CO_12h00',
   CEMIS_PGD_AREA(2)='LAN',
   CEMIS_PGD_ATYPE(2)='ARI',
   CEMIS_PGD_FILE(2)='co_12.asc',
   CEMIS_PGD_FILETYPE(2)='ASCLLV',
   CEMIS_PGD_NAME(3)='DMSE',
   NEMIS_PGD_TIME(3)=0,
   CEMIS_PGD_COMMENT(3)='dms_cte',
   CEMIS_PGD_AREA(3)='SEA',
   CEMIS_PGD_ATYPE(3)='ARI',
   CEMIS_PGD_FILE(3)='dms.asc',
   CEMIS_PGD_FILETYPE(3)='ASCLLV'
/

NAM_DUMMY_PGD

Warning

This namelist comes from SURFEX 9.0.0 user guide https://www.umr-cnrm.fr/surfex/IMG/pdf/surfex_tecdoc.pdf.

This namelist allows to incorporate into the physiographic file any surface field. You can treat up to 999 such fields. These fields will be written on all the files you will use later(after prognostic fields initialization, or during and after run, etc…). Their name in the files are ‘DUMMY_GRnnn’, where nnn goes from 001 to 999. During the execution of the programs, these fields are stored in the XDUMMY_FIELDS(:,:) (first dimension: spatial dimension, second dimension: total number of fields), in the module MODD_DUMMY_SURF_FIELDn. You must modify the fortran source, where you want to use them.

NAM_DUMMY_PGD content
Fortran name	Fortran type	Default value
NDUMMY_NBR	INTEGER	0
CDUMMY_NAME(:)	CHARACTER(LEN=20)
CDUMMY_FILE(:)	CHARACTER(LEN=28)
CDUMMY_FILETYPE(:)	CHARACTER(LEN=6)
CDUMMY_AREA(:)	CHARACTER(LEN=3)	‘ALL’
CDUMMY_ATYPE(:)	CHARACTER(LEN=3)	‘ARI’

Only the first NDUMMY_NBR values in these arrays are meaningfull.

NDUMMY_NBR : number of dummy fields.
CDUMMY_NAME(:) : list of the dummy fields you want to initialize with your own data. You can give any name you want. This is a way to describe what is the field. This information is not used by the program. It is just written in the Meso-NH files.
CDUMMY_FILE(:) : list of the names of the files containing the data for the fields you have specified in CDUMMY_NAME(:).
CDUMMY_FILETYPE(:) : list of the types of the files containing the data for the fields you have specified in CDUMMY_NAME(:) (‘DIRECT’, ‘LATLON’, ‘BINLLF’, ‘BINLLV’, ‘ASCLLV’).
CDUMMY_AREA(:) : area of meaningfullness of the fields you have specified in CDUMMY_NAME(:) (‘ALL’, ‘NAT’, ‘TWN’, ‘SEA’, ‘WAT’, ‘LAN’, respectively for everywhere, natural areas, town areas, sea, inland waters, land = natural cover + town). For example, oceanic emission of DNS is relevant on ‘SEA’.
CDUMMY_ATYPE(:) : type of averaging (during PGD for the fields you have specified in CDUMMY_NAME(:) (‘ARI’, ‘INV’, ‘LOG’, respectively for arithmetic, inverse and logarithmic averaging).

NAM_PREP_SURF_ATM

Warning

This namelist comes from SURFEX 9.0.0 user guide https://www.umr-cnrm.fr/surfex/IMG/pdf/surfex_tecdoc.pdf.

This namelist is used to:

initialize the date of all surface schemes. The namelist information is used only if no input data file is used, either from namelist or by fortran code (as in MESONH program). If a file is used, the date is read in it.
define the default file in which each scheme can read the needed data (e.g. temperature). Note that, all the information given in this namelist can be erased for each scheme by the namelist corresponding to this scheme, as the information in the shceme namelists have priority on namelist NAM_PREP_SURF_ATM.

NAM_PREP_SURF_ATM content
Fortran name	Fortran type	Default value
CFILE	CHARACTER(LEN=28)
CFILETYPE	CHARACTER(LEN=6)
CFILEPGD	CHARACTER(LEN=28)
CFILEPGDTYPE	CHARACTER(LEN=6)
NYEAR	INTEGER
NMONTH	INTEGER
NDAY	INTEGER
XTIME	REAL
NHALO_PREP	INTEGER	0
LWRITE_EXTERN	LOGICAL	.FALSE.

CFILE / CFILEPGD : name of the prep / pgd file used to define the date and the file in which to read the needed data (e.g. temperature).

Note

The use of a file or prescribed value in each scheme namelist has priority on the data in CFILE / CFILEPGD file of namelist NAM_PREP_SURF_ATM.
CFILE and CFILEPGD can identify the same file.

CFILETYPE / CFILEPGDTYPE : type of the CFILE / CFILEPGD file, if the latter is provided. The following values are currently usable:
- “MESONH”: the file type is a MESONH file.
- “GRIB”: the file type is a GRIB file, coming from any of these models:
  - “ECMWF “: european center forecast model
  - “ARPEGE”: Arpege french forecast model
  - “AROME”: AROME french forecast local model
  - “MOCAGE”: Mocage french research chemistry model
- “ASCII “: ASCII Surfex PREP/PGD file
- “LFI “: LFI Surfex PREP/PGD file
NYEAR : year of surface UTC time. It is used only if no atmospheric file or no surface file is given (in those the date can be read).
NMONTH : month of surface UTC time. It is used only if no atmospheric file or no surface file is given (in those the date can be read).
NDAY : day of surface UTC time. It is used only if no atmospheric file or no surface file is given (in those the date can be read).
XTIME : time from midnight of surface UTC time. It is used only if no atmospheric file or no surface file is given (in those the date can be read). (seconds).
NHALO_PREP : halo size for the extrapolation of pronostic fields from input file
LWRITE_EXTERN : The new key LWRITE_EXTERN is added. If LWRITE_EXTERN=T, soil depths for ISBA and TEB are written in the current output PREP file.

NAM_PREP_SEAFLUX

Warning

This namelist comes from SURFEX 9.0.0 user guide https://www.umr-cnrm.fr/surfex/IMG/pdf/surfex_tecdoc.pdf.

This namelist information is used to initialize the “SEAFLX” sea scheme temperature.

NAM_PREP_SEAFLUX content
Fortran name	Fortran type	Default value
XSST_UNIF	REAL
CFILE_SEAFLX	CHARACTER(LEN=28)
CTYPE_SEAFLX	CHARACTER(LEN=6)
CFILEPGD_SEAFLX	CHARACTER(LEN=28)
CTYPEPGD	CHARACTER(LEN=6)
CFILEWAVE_SEAFLX	CHARACTER(LEN=28)
CTYPEWAVE	CHARACTER(LEN=6)
NYEAR	INTEGER
NMONTH	INTEGER
NDAY	INTEGER
XTIME	REAL
LSEA_SBL	LOGICAL	.FALSE.
LOCEAN_MERCATOR	LOGICAL	.FALSE.
LOCEAN_CURRENT	LOGICAL	.FALSE.
XTIME_REL	REAL	25920000.0
LCUR_REL	LOGICAL	.FALSE.
LTS_REL	LOGICAL	.FALSE.
LZERO_FLUX	LOGICAL	.FALSE.
LCORR_FLUX	LOGICAL	.FALSE.
XCORFLX	REAL	0.0
LDIAPYC	LOGICAL	.FALSE.
CSEAICE_SCHEME	CHARACTER(LEN=6)
XSSS_UNIF	REAL	1.E+20
XSIC_UNIF	REAL	1.E+20

XSST_UNIF : uniform prescribed value of Sea Surface Temperature. This prescribed value, if defined, has priority on the use of CFILE_SEAFLX data.
CFILE_SEAFLX / CFILEPGD_SEAFLX : name of the PREP/PGD files used to define the Sea surface Temperature. The use of a file or prescribed value XSST_UNIF has priority on the data in CFILE_SEAFLX file.
CTYPE_SEAFLX / CTYPEPGD : type of the CFILE_SEAFLX / CFILEPGD_SEAFLX files, if the latter is provided. The following values are currently usable:
- “MESONH”: the file type is a MESONH file.
- “GRIB”: the file type is a GRIB file, coming from any of these models:
  - “ECMWF “: european center forecast model
  - “ARPEGE”: Arpege french forecast model
  - “AROME”: AROME french forecast local model
  - “MOCAGE”: Mocage french research chemistry model
- “NETCDF”: the file type is a NETCDF file, coming from MERCATOR (possible only for CTYPE_SEAFLX)
- “ASCII “: PREP/PGD Surfex ASCII file
- “LFI “: PREP/PGD Surfex LFI file
CFILEWAVE_SEAFLX : name of the file used to define the significant wave height (Hs) and the peak period (Tp)
CTYPEWAVE : type of the CFILEWAVE_SEAFLX file if the latter is provided. CTYPEWAVE must be given. The ‘NETCDF’ value (if the file type is a netcdf file) is the only one usable. Other ‘PREP_SEAFLUX’ types are under development and lead now to uniform values of wave parameters.
NYEAR : year of surface UTC time. It is used only if no atmospheric file or no surface file is given (in those the date can be read).
NMONTH : month of surface UTC time. It is used only if no atmospheric file or no surface file is given (in those the date can be read).
NDAY : day of surface UTC time. It is used only if no atmospheric file or no surface file is given (in those the date can be read).
XTIME : time from midnight of surface UTC time. It is used only if no atmospheric file or no surface file is given (in those the date can be read). (seconds).
LSEA_SBL : activates surface boundary multi layer scheme over sea.
LOCEAN_MERCATOR : oceanic variables initialized from MERCATOR if T
LOCEAN_CURRENT : initial ocean state with current (if F ucur=0, vcur=0)
XTIME_REL : time of relaxation (s)
LCUR_REL : flag for relaxation on current
LTS_REL : flag for relaxation on ocean temperature
LZERO_FLUX : flag for testing zero incoming flux at the ocean surface
LCORR_FLUX : flag for flux correction
XCORFLX : correction coefficient for surface fluxes
LDIAPYC : flag for diapycnal mixing activation
XSSS_UNIF : from V8, uniform prescribed value of Sea Surface Salinity. This prescribed value, if defined, has priority on the use of CFILE_SEAFLX data.
CSEAICE_SCHEME : from V8, name of the sea-ice scheme to activate.
XSIC_UNIF : uniform sea ice covert fraction

NAM_PREP_WATFLUX

Warning

This namelist comes from SURFEX 9.0.0 user guide https://www.umr-cnrm.fr/surfex/IMG/pdf/surfex_tecdoc.pdf.

This namelist information is used to initialize the “WATFLX” sea scheme temperature.

NAM_PREP_WATFLUX content
Fortran name	Fortran type	Default value
XTS_WATER_UNIF	REAL
CFILE_WATFLX	CHARACTER(LEN=28)
CTYPE	CHARACTER(LEN=6)
CFILEPGD_WATFLX	CHARACTER(LEN=28)
CTYPEPGD	CHARACTER(LEN=6)
NYEAR	INTEGER
NMONTH	INTEGER
NDAY	INTEGER
XTIME	REAL
LWAT_SBL	LOGICAL	.FALSE.

XTS_WATER_UNIF : uniform prescribed value of water surface temperature supposed at an altitude of 0m (mean sea level altitude). The temperature is then modified for each point depending on its altitude, following a uniform vertical gradient of -6.5 K km-1. This prescribed value, if defined, has priority on the use of CFILE_WATFLX data.
CFILE_WATFLX / CFILEPGD_WATFLX : name of the PREP / PGD files used to define the Sea surface Temperature. The use of a file or prescribed value XTS_WATER_UNIF has priority on the data in CFILE_WATFLX file.
CTYPE / CTYPEPGD : type of the CFILE_WATFLX / CFILEPGD_WATFLX file, if the latter is provided. CTYPE / CTYPEPGD must then be given. The following values are currently usable:
- “MESONH”: the file type is a MESONH file.
- “GRIB”: the file type is a GRIB file, coming from any of these models:
  - “ECMWF “: european center forecast model
  - “ARPEGE”: Arpege french forecast model
  - “AROME”: AROME french forecast local model
  - “MOCAGE”: Mocage french research chemistry model
- “ASCII “: PREP / PGD Surfex ASCII file
- “LFI “: PREP/PGD Surfex LFI file
NYEAR : year of surface UTC time. It is used only if no atmospheric file or no surface file is given (in those the date can be read).
NMONTH : month of surface UTC time. It is used only if no atmospheric file or no surface file is given (in those the date can be read).
NDAY : day of surface UTC time. It is used only if no atmospheric file or no surface file is given (in those the date can be read).
XTIME : time from midnight of surface UTC time (seconds). It is used only if no atmospheric file or no surface file is given (in those the date can be read).
LWAT_SBL : activates surface boundary multi layer scheme over inland water.

NAM_PREP_GREENROOF_SNOW

NAM_PREP_GREENROOF_SNOW content
Fortran name	Fortran type	Default value
CSNOW_GR	CHARACTER(LEN=3)	‘3-L’
NSNOW_LAYER_GR	INTEGER	3
CFILE_SNOW_GR	CHARACTER(LEN=28)	CFILE_GR in
		{footnotesize NAM_PREP_TEB_GREENROOF}
CTYPE_SNOW	CHARACTER(LEN=6)	CTYPE in
		{footnotesize NAM_PREP_TEB_GREENROOF}
CFILEPGD_SNOW_GR	CHARACTER(LEN=28)	CFILEPGD_GR in
		{footnotesize NAM_PREP_TEB_GREENROOF}
CTYPEPGD_SNOW	CHARACTER(LEN=6)	CTYPEPGD in
		{footnotesize NAM_PREP_TEB_GREENROOF}
LSNOW_IDEAL_GR	LOGICAL	F
XWSNOW_GR	REAL(20)
XZSNOW_GR	REAL(20)	1.E+20
XTSNOW_GR	REAL(20)	273.16
XLWCSNOW_GR	REAL(20)
XRSNOW_GR	REAL(20)
XASNOW_GR	REAL	0.5

CSNOW_GR : type of snow scheme. Possible snow schemes are:
- ‘D95’: Douville et al (1995) snow scheme.
- ‘3-L’: Boone and Etchevers (2000) three layers snow scheme.
- ‘EBA’: Bogatchev and Bazile (2005), Arpege operational snow scheme.
NSNOW_LAYER_GR : number of snow layers
CFILE_SNOW_GR : name of the file used to define the snow profiles. The use of a file or prescribed value of XRSNOW_GR, XTSNOW_GR, XWSNOW_GR and XASNOW_GR has priority on the data in CFILE_SNOW_GR file
CTYPE_SNOW : type of the CFILE_SNOW_GR file, if the latter is provided. CTYPE_SNOW must then be given. The following values are currently usable:
- “MESONH”: the file type is a MESONH file.
- “GRIB “: the file type is a GRIB file, coming from any of these models: * “ECMWF “: european center forecast model * “ARPEGE”: Arpege french forecast model * “AROME”: Arome french forecast local model * “MOCAGE”: Mocage french research chemistry model
- “LFI “: LFI PREP file
- “ASCII”: ASCII PREP FILE
CFILEPGD_SNOW_GR : name of the associated PGD file if CFILE_SNOW_GR is a PREP files.
CTYPEPGD_SNOW : type of the CFILEPGD_SNOW file, if the latter is provided. CTYPEPGD_SNOW must then be given. The following values are currently usable:
- “MESONH”: the file type is a MESONH file.
- “LFI “: LFI PREP file
- “ASCII”: ASCII PREP FILE
LSNOW_IDEAL_GR : if LSNOW_IDEAL_GR = F , only one value can be given for following snow parameters and a vertical interpolation is processed. If LSNOW_IDEAL_GR = T, values are given for each layer and there is no vertical interpolation performed.
XWSNOW_GR : uniform value to initialize snow content, one for each layer
XZSNOW_GR : uniform value to initialize snow depth, one for each layer (m)
XTSNOW_GR : uniform value to initialize snow temperature, one for each layer
XLWCSNOW_GR : uniform value to initialize snow liquid water content, one for each layer (kg/m3)
XRSNOW_GR : uniform value to initialize snow density, one for each layer
XASNOW_GR : uniform value to initialize snow albedo

NAM_PREP_GARDEN_SNOW

NAM_PREP_GARDEN_SNOW content
Fortran name	Fortran type	Default value
CSNOW_GD	CHARACTER(LEN=3)	‘D95’
NSNOW_LAYER_GD	INTEGER	1
CFILE_SNOW_GD	CHARACTER(LEN=28)	CFILE_GD in
		NAM_PREP_TEB_GARDEN
CTYPE_SNOW	CHARACTER(LEN=6)	CTYPE in
		NAM_PREP_TEB_GARDEN
CFILEPGD_SNOW_GD	CHARACTER(LEN=28)	CFILEPGD_GD in
		NAM_PREP_TEB_GARDEN
CTYPEPGD_SNOW	CHARACTER(LEN=6)	CTYPEPGD in
		NAM_PREP_TEB_GARDEN
LSNOW_IDEAL_GD	LOGICAL	F
XWSNOW_GD	REAL(20)
XZSNOW_GD	REAL(20)	1.E+20
XTSNOW_GD	REAL(20)	273.16
XLWCSNOW_GD	REAL(20)
XRSNOW_GD	REAL(20)
XASNOW_GD	REAL	0.5

CSNOW_GD : type of snow scheme. Possible snow schemes are:
- ‘D95’: Douville et al (1995) snow scheme.
- ‘3-L’: Boone and Etchevers (2000) three layers snow scheme.
- ‘EBA’: Bogatchev and Bazile (2005), Arpege operational snow scheme.
NSNOW_LAYER_GD : number of snow layers
CFILE_SNOW_GD : name of the file used to define the snow profiles. The use of a file or prescribed value of XRSNOW_GD, XTSNOW_GD, XWSNOW_GD and XASNOW_GD has priority on the data in CFILE_SNOW_GD file
CTYPE_SNOW : type of the CFILE_SNOW_GD file, if the latter is provided. CTYPE_SNOW must then be given. The following values are currently usable:
- “MESONH”: the file type is a MESONH file.
- “GRIB “: the file type is a GRIB file, coming from any of these models: * “ECMWF “: european center forecast model * “ARPEGE”: Arpege french forecast model * “AROME”: Arome french forecast local model * “MOCAGE”: Mocage french research chemistry model
- “LFI “: LFI PREP file
- “ASCII”: ASCII PREP FILE
CFILEPGD_SNOW_GD : name of the associated PGD file if CFILE_SNOW_GD is a PREP files.
CTYPEPGD_SNOW : type of the CFILEPGD_SNOW file, if the latter is provided. CTYPEPGD_SNOW must then be given. The following values are currently usable:
- “MESONH”: the file type is a MESONH file.
- “LFI “: LFI PREP file
- “ASCII”: ASCII PREP FILE
LSNOW_IDEAL_GD : if LSNOW_IDEAL_GD = F , only one value can be given for following snow parameters and a vertical interpolation is processed. If LSNOW_IDEAL_GD = T, values are given for each layer and there is no vertical interpolation performed.
XWSNOW_GD : uniform value to initialize snow content, one for each layer (kg/m$^{2}$)
XZSNOW_GD : uniform value to initialize snow depth, one for each layer (m) (alternative to XWSNOW_GD)
XTSNOW_GD : uniform value to initialize snow temperature, one for each layer (K)
XLWCSNOW_GD : uniform value to initialize liquid snow water contents, one for each layer (kg/m$^{3}$)
XRSNOW_GD : uniform value to initialize snow density, one for each layer (kg/m$^{3}$)
XASNOW_GD : uniform value to initialize snow albedo (-)

NAM_PREP_TEB

Warning

This namelist comes from SURFEX 9.0.0 user guide https://www.umr-cnrm.fr/surfex/IMG/pdf/surfex_tecdoc.pdf.

This namelist information is used to initialize the “TEB “ urban scheme variables: road, roof and wall temperature profiles, water intercepted by roofs and roads, snow, building internal temperature.

NAM_PREP_TEB content
Fortran name	Fortran type	Default value
XWS_ROAD	REAL	‘NONE’
XWS_ROOF	REAL	‘NONE’
CFILE_WS	CHARACTER(LEN=28)
CTYPE_WS	CHARACTER(LEN=6)
XTS_ROAD	REAL	‘NONE’
XTS_ROOF	REAL	‘NONE’
XTS_WALL	REAL	‘NONE’
XTI_BLD	REAL	‘NONE’
XHUI_BLD	REAL	‘NONE’
CROAD_DIR	CHARACTER(LEN=4)
CWALL_OPT	CHARACTER(LEN=4)
CFILE_TS	CHARACTER(LEN=28)
CTYPE_TS	CHARACTER(LEN=6)
CFILE_TEB	CHARACTER(LEN=28)
CTYPE	CHARACTER(LEN=6)
CFILEPGD_TEB	CHARACTER(LEN=28)
CTYPEPGD	CHARACTER(LEN=6)
NYEAR	INTEGER	‘NONE’
NMONTH	INTEGER	‘NONE’
NDAY	INTEGER	‘NONE’
XTIME	REAL	‘NONE’
LTEB_CANOPY	LOGICAL	.FALSE.
LATM_CANOPY	LOGICAL	.FALSE.
XTDEEP_TEB	REAL	1.E+20
XTS_BLD	REAL	17+XTT

XWS_ROAD : uniform prescribed value of soil water interception for the road reservoir. This prescribed value, if defined, has priority on the use of CFILE_WS and CFILE_TEB data.
XWS_ROOF : uniform prescribed value of soil water interception for the roof reservoir. This prescribed value, if defined, has priority on the use of CFILE_WS and CFILE_TEB data.
CFILE_WS : name of the file used to define the soil water reservoirs. The use of a file or prescribed value of XWS_ROAD and XWS_ROOF has priority on the data in CFILE_WS file.
CTYPE_WS : type of the CFILE_WS file, if the latter is provided. CTYPE_WS must then be given. The following values are currently usable:
- “MESONH”: the file type is a MESONH file.
- “GRIB “: the file type is a GRIB file, coming from any of these models:
  - “ECMWF “: european center forecast model
  - “ARPEGE”: Arpege french forecast model
  - “AROME”: Arome french forecast local model
  - “MOCAGE”: Mocage french research chemistry model
XTS_ROAD : uniform prescribed value of temperature for road, supposed at an altitude of 0m (mean sea level altitude). The temperature is then modified for each point depending on its altitude, following a uniform vertical gradient of -6.5 K km-1. This prescribed value, if defined, has priority on the use of CFILE_TS and CFILE_TEB data.
XTS_ROOF : uniform prescribed value of temperature for roof, supposed at an altitude of 0m (mean sea level altitude). The temperature is then modified for each point depending on its altitude, following a uniform vertical gradient of -6.5 K km-1. This prescribed value, if defined, has priority on the use of CFILE_TS and CFILE_TEB data.
XTS_WALL : uniform prescribed value of temperature for wall, supposed at an altitude of 0m (mean sea level altitude). The temperature is then modified for each point depending on its altitude, following a uniform vertical gradient of -6.5 K km-1. This prescribed value, if defined, has priority on the use of CFILE_TS and CFILE_TEB data.
XTI_BLD : uniform prescribed value of internal building temperature. This temperature is not dependent on altitude. This prescribed value, if defined, has priority on the use of CFILE_TS and CFILE_TEB data.
XHUI_BLD : uniform bulding relative hum (between 0-1)
CROAD_DIR : TEB option for road direction:
- “UNIF” : no specific direction
- “ORIE” : many road ORIEntations (linked to NTEB_PATCH)
CWALL_OPT : TEB option for walls:
- “UNIF” : uniform walls
- “TWO” : two separated walls
CFILE_TS : name of the file used to define the soil temperature profile. The use of a file or prescribed value of XTS_ROAD, XTS_ROOF, XTS_WALL, XTI_BLD or XTI_ROAD has priority on the data in CFILE_TS file.
CTYPE_TS : type of the CFILE_TS file, if the latter is provided. CTYPE_TS must then be given. The following values are currently usable:
- “MESONH”: the file type is a MESONH file.
- “GRIB “: the file type is a GRIB file, coming from any of these models:
  - “ECMWF “: european center forecast model
  - “ARPEGE”: Arpege french forecast model
  - “AROME”: Arome french forecast local model
  - “MOCAGE”: Mocage french research chemistry model
CFILE_TEB / CFILEPGD_TEB : name of the PREP/PGD files used to define any TEB variable. The use of a file or prescribed value XWS_ROAD, XWS_ROOF, XTS_ROAD, XTS_ROOF, XTS_WALL, XTI_BLD, XTI_ROAD, CFILE_WS or CFILE_TS has priority on the data in CFILE_TEB file.
CTYPE / CTYPEPGD : type of the CFILE_TEB / CFILEPGD_TEB file, if the latter is provided. CTYPE must then be given. The following values are currently usable:
- “MESONH”: the file type is a MESONH file.
- “GRIB “: the file type is a GRIB file, coming from any of these models:
  - “ECMWF “: european center forecast model
  - “ARPEGE”: Arpege french forecast model
  - “AROME”: Arome french forecast local model
  - “MOCAGE”: Mocage french research chemistry model
- “ASCII “: PREP/PGD Surfex ASCII file
- “LFI “: PREP/PGD Surfex LFI file
NYEAR : year of surface UTC time. It is used only if no atmospheric file or no surface file is given (in those the date can be read).
NMONTH : month of surface UTC time. It is used only if no atmospheric file or no surface file is given (in those the date can be read).
NDAY : day of surface UTC time. It is used only if no atmospheric file or no surface file is given (in those the date can be read).
XTIME : time from midnight of surface UTC time. It is used only if no atmospheric file or no surface file is given (in those the date can be read). (seconds).
LTEB_CANOPY : activates surface boundary multi layer scheme over town.
LATM_CANOPY : flag to replace canopy prognostic variables by atmospheric models prognostic variables
XTDEEP_TEB : deep temperature for TEB soil (K). XTS_BLD corresponds to the old key XTI_ROAD (before V9)
XTS_BLD : soil below buildings uniform temperature (K). Default value is 17+XTT (XTT for triple point temperature = 273.16K).

NAM_PREP_TEB_GARDEN

NAM_PREP_TEB_GARDEN content
Fortran name	Fortran type	Default value
XHUG_SURF_GD	REAL	none
XHUG_ROOT_GD	REAL	none
XHUG_DEEP_GD	REAL	none
XHUGI_SURF_GD	REAL	none
XHUGI_ROOT_GD	REAL	none
XHUGI_DEEP_GD	REAL	none
CFILE_HUG_SURF_GD	CHARACTER(LEN=28)	CFILE_HUG_GD
		in this namelist
CFILE_HUG_ROOT_GD	CHARACTER(LEN=28)	CFILE_HUG_GD
		in this namelist
CFILE_HUG_DEEP_GD	CHARACTER(LEN=28)	CFILE_HUG_GD
		in this namelist
CFILE_HUG_GD	CHARACTER(LEN=28)	CFILE_GD
		in this namelist
CTYPE_HUG	CHARACTER(LEN=6)	CTYPE
		in this namelist

XTG_SURF_GD	REAL	none
XTG_ROOT_GD	REAL	none
XTG_DEEP_GD	REAL	none
CFILE_TG_SURF_GD	CHARACTER(LEN=28)	CFILE_TG_GD
		in this namelist
CFILE_TG_ROOT_GD	CHARACTER(LEN=28)	CFILE_TG_GD
		in this namelist
CFILE_TG_DEEP_GD	CHARACTER(LEN=28)	CFILE_TG_GD
		in this namelist
CFILE_TG_GD	CHARACTER(LEN=28)	CFILE_GD
		in this namelist
CTYPE_TG	CHARACTER(LEN=6)	CTYPE
		in this namelist

CFILE_GD	CHARACTER(LEN=28)	CFILE in
		NAM_PREP_SURF_ATM
CTYPE	CHARACTER(LEN=6)	CFILETYPE in
		NAM_PREP_SURF_ATM
CFILEPGD_GD	CHARACTER(LEN=28)	CFILEPGD in
		NAM_PREP_SURF_ATM
CTYPEPGD	CHARACTER(LEN=6)	CFILEPGDTYPE in
		NAM_PREP_SURF_ATM

XHUG_SURF_GD : uniform prescribed value of liquid soil water index (SWI) for the surface soil layer. This prescribed value, if defined, has priority on the use of CFILE_HUG_GD and CFILE_GD data.
XHUG_ROOT_GD : uniform prescribed value of liquid soil water index (SWI) for the root zone soil layer(s). This prescribed value, if defined, has priority on the use of CFILE_HUG_GD and CFILE_GD data.
XHUG_DEEP_GD : uniform prescribed value of liquid soil water index (SWI) for the deep soil layer(s). This prescribed value, if defined, has priority on the use of CFILE_HUG_GD and CFILE_GD data.
XHUGI_SURF_GD : uniform prescribed value of ice soil water index (SWI) for the surface soil layer. This prescribed value, if defined, has priority on the use of CFILE_HUG_GD and CFILE_GD data.
XHUGI_ROOT_GD : uniform prescribed value of ice soil water index (SWI) for the root zone soil layer(s). This prescribed value, if defined, has priority on the use of CFILE_HUG_GD and CFILE_GD data.
XHUGI_DEEP_GD : uniform prescribed value of ice soil water index (SWI) for the deep soil layer(s). This prescribed value, if defined, has priority on the use of CFILE_HUG_GD and CFILE_GD data.
CFILE_HUG_SURF_GD : name of the file used to define the liquid soil water index (SWI) for the surface soil layer.
CFILE_HUG_ROOT_GD : name of the file used to define the liquid soil water index (SWI) for the root zone soil layer(s).
CFILE_HUG_DEEP_GD : name of the file used to define the liquid soil water index (SWI) for the deep soil layer(s).
CFILE_HUG_GD : name of the file used to define the soil water profiles.\
CTYPE_HUG : type of the CFILE_HUG_GD file, if the latter is provided. CTYPE_HUG must then be given. The following values are currently usable:
- “MESONH”: the file type is a MESONH file.
- “GRIB “: the file type is a GRIB file, coming from any of these models: * “ECMWF “: european center forecast model * “ARPEGE”: Arpege french forecast model * “AROME”: Arome french forecast local model * “MOCAGE”: Mocage french research chemistry model
- “ASCII / LFI “: PREP file from Surfex
- “ASCLLV”: ASCII latlonval file (one file for each depth)
XTG_SURF_GD : uniform prescribed value of temperature for the surface soil layer, supposed at an altitude of 0m (mean sea level altitude). The temperature is then modified for each point depending on its altitude, following a uniform vertical gradient of -6.5 K km-1. This prescribed value, if defined, has priority on the use of CFILE_TG_GD and CFILE_GD data.
XTG_ROOT_GD : uniform prescribed value of temperature for the root zone soil layer(s), supposed at an altitude of 0m (mean sea level altitude). The temperature is then modified for each point depending on its altitude, following a uniform vertical gradient of -6.5 K km-1. This prescribed value, if defined, has priority on the use of CFILE_TG_GD and CFILE_GD data.
XTG_DEEP_GD : uniform prescribed value of temperature for the deep soil layer(s), supposed at an altitude of 0m (mean sea level altitude). The temperature is then modified for each point depending on its altitude, following a uniform vertical gradient of -6.5 K km-1. This prescribed value, if defined, has priority on the use of CFILE_TG_GD and CFILE_GD data.
CFILE_TG_SURF_GD : name of the file used to define the surface soil temperature profile.
CFILE_TG_ROOT_GD : name of the file used to define the root zone soil temperature profile.
CFILE_TG_DEEP_GD : name of the file used to define the deep soil temperature profile.
CFILE_TG_GD : name of the file used to define the soil temperature profile.
The use of a file or prescribed value of XTG_SURF_GD, XTG_ROOT_GD and XTG_DEEP_GD has priority on the data in CFILE_TG_GD file.
CTYPE_TG : type of the CFILE_TG_GD file, if the latter is provided. CTYPE_TG must then be given. The following values are currently usable:
- “MESONH”: the file type is a MESONH file.
- “GRIB “: the file type is a GRIB file, coming from any of these models: * “ECMWF “: european center forecast model * “ARPEGE”: Arpege french forecast model * “AROME”: Arome french forecast local model * “MOCAGE”: Mocage french research chemistry model
- “ASCII / LFI “: PREP file from Surfex
- “ASCLLV”: ASCII latlonval file (one file for each depth)
CFILE_GD / CFILEPGD_GD : name of the PREP / PGD files used to define any GARDEN variable. The use of a file or prescribed value XHUG_SURF_GD, XHUG_ROOT_GD, XHUG_DEEP_GD, XTG_SURF_GD, XTG_ROOT_GD, XTG_DEEP_GD, CFILE_WG_GD and CFILE_TG_GD has priority on the data in CFILE_GD file.
CTYPE / CTYPEPGD : type of the CFILE_GD / CFILEPGD_GD files, if the latter is provided. CTYPE / CTYPEPGD must then be given. The following values are currently usable:
- “MESONH”: the file type is a MESONH file.
- “GRIB “: the file type is a GRIB file, coming from any of these models: * “ECMWF “: european center forecast model * “ARPEGE”: Arpege french forecast model * “AROME”: Arome french forecast local model * “MOCAGE”: Mocage french research chemistry model
- “ASCII “: PREP/PGD Surfex ASCII file
- “LFI “: PREP/PGD Surfex LFI file

NAM_PREP_TEB_GREENROOF

NAM_PREP_TEB_GREENROOF content
Fortran name	Fortran type	Default value
XHUG_SURF_GR	REAL	none
XHUG_ROOT_GR	REAL	none
XHUG_DEEP_GR	REAL	none
XHUGI_SURF_GR	REAL	none
XHUGI_ROOT_GR	REAL	none
XHUGI_DEEP_GR	REAL	none
CFILE_HUG_SURF_GR	CHARACTER(LEN=28)	CFILE_HUG_GR
		in this namelist
CFILE_HUG_ROOT_GR	CHARACTER(LEN=28)	CFILE_HUG_GR
		in this namelist
CFILE_HUG_DEEP_GR	CHARACTER(LEN=28)	CFILE_HUG_GR
		in this namelist
CFILE_HUG_GR	CHARACTER(LEN=28)	CFILE_GR
		in this namelist
CTYPE_HUG	CHARACTER(LEN=6)	CTYPE
		in this namelist

XTG_SURF_GR	REAL	none
XTG_ROOT_GR	REAL	none
XTG_DEEP_GR	REAL	none
CFILE_TG_SURF_GR	CHARACTER(LEN=28)	CFILE_TG_GR
		in this namelist
CFILE_TG_ROOT_GR	CHARACTER(LEN=28)	CFILE_TG_GR
		in this namelist
CFILE_TG_DEEP_GR	CHARACTER(LEN=28)	CFILE_TG_GR
		in this namelist
CFILE_TG_GR	CHARACTER(LEN=28)	CFILE_GR
		in this namelist
CTYPE_TG	CHARACTER(LEN=6)	CTYPE
		in this namelist
		in this namelist
CFILE_GR	CHARACTER(LEN=28)	CFILE in
		NAM_PREP_SURF_ATM
CTYPE	CHARACTER(LEN=6)	CFILETYPE in
		NAM_PREP_SURF_ATM
CFILEPGD_GR	CHARACTER(LEN=28)	CFILEPGD in
		NAM_PREP_SURF_ATM
CTYPEPGD	CHARACTER(LEN=6)	CFILEPGDTYPE in
		NAM_PREP_SURF_ATM

XHUG_SURF_GR : uniform prescribed value of liquid soil water index (SWI) for the surface soil layer. This prescribed value, if defined, has priority on the use of CFILE_HUG_GR and CFILE_GR data.
XHUG_ROOT_GR : uniform prescribed value of liquid soil water index (SWI) for the root zone soil layer(s). This prescribed value, if defined, has priority on the use of CFILE_HUG_GR and CFILE_GR data.
XHUG_DEEP_GR : uniform prescribed value of liquid soil water index (SWI) for the deep soil layer(s). This prescribed value, if defined, has priority on the use of CFILE_HUG_GR and CFILE_GR data.
XHUGI_SURF_GR : uniform prescribed value of ice soil water index (SWI) for the surface soil layer. This prescribed value, if defined, has priority on the use of CFILE_HUG_GR and CFILE_GR data.
XHUGI_ROOT_GR : uniform prescribed value of ice soil water index (SWI) for the root zone soil layer(s). This prescribed value, if defined, has priority on the use of CFILE_HUG_GR and CFILE_GR data.
XHUGI_DEEP_GR : uniform prescribed value of ice soil water index (SWI) for the deep soil layer(s). This prescribed value, if defined, has priority on the use of CFILE_HUG_GR and CFILE_GR data.
CFILE_HUG_SURF_GR : name of the file used to define the liquid soil water index (SWI) for the surface soil layer.
CFILE_HUG_ROOT_GR : name of the file used to define the liquid soil water index (SWI) for the root zone soil layer(s).
CFILE_HUG_DEEP_GR : name of the file used to define the liquid soil water index (SWI) for the deep soil layer(s).
CFILE_HUG_GR : name of the file used to define the soil water profiles.\
CTYPE_HUG : type of the CFILE_HUG_GR file, if the latter is provided. CTYPE_HUG must then be given. The following values are currently usable:
- “MESONH”: the file type is a MESONH file.
- “GRIB “: the file type is a GRIB file, coming from any of these models: * “ECMWF “: european center forecast model * “ARPEGE”: Arpege french forecast model * “AROME”: Arome french forecast local model * “MOCAGE”: Mocage french research chemistry model
- “ASCII / LFI “: PREP file from Surfex
- “ASCLLV”: ASCII latlonval file (one file for each depth)
XTG_SURF_GR : uniform prescribed value of temperature for the surface soil layer, supposed at an altitude of 0m (mean sea level altitude). The temperature is then modified for each point depending on its altitude, following a uniform vertical gradient of -6.5 K km-1. This prescribed value, if defined, has priority on the use of CFILE_TG_GR and CFILE_GR data.
XTG_ROOT_GR : uniform prescribed value of temperature for the root zone soil layer(s), supposed at an altitude of 0m (mean sea level altitude). The temperature is then modified for each point depending on its altitude, following a uniform vertical gradient of -6.5 K km-1. This prescribed value, if defined, has priority on the use of CFILE_TG_GR and CFILE_GR data.
XTG_DEEP_GR : uniform prescribed value of temperature for the deep soil layer(s), supposed at an altitude of 0m (mean sea level altitude). The temperature is then modified for each point depending on its altitude, following a uniform vertical gradient of -6.5 K km-1. This prescribed value, if defined, has priority on the use of CFILE_TG_GR and CFILE_GR data.
CFILE_TG_SURF_GR : name of the file used to define the surface soil temperature profile.
CFILE_TG_ROOT_GR : name of the file used to define the root zone soil temperature profile.
CFILE_TG_DEEP_GR : name of the file used to define the deep soil temperature profile.
CFILE_TG_GR : name of the file used to define the soil temperature profile.\
CTYPE_TG : type of the CFILE_TG_GR file, if the latter is provided. CTYPE_TG must then be given. The following values are currently usable:
- “MESONH”: the file type is a MESONH file.
- “GRIB “: the file type is a GRIB file, coming from any of these models: * “ECMWF “: european center forecast model * “ARPEGE”: Arpege french forecast model * “AROME”: Arome french forecast local model * “MOCAGE”: Mocage french research chemistry model
- “ASCII / LFI “: PREP file from Surfex
- “ASCLLV”: ASCII latlonval file (one file for each depth)
CFILE_GR / CFILEPGD_GR : name of the PREP / PGD files used to define any GARDEN variable. The use of a file or prescribed value XHUG_SURF_GR, XHUG_ROOT_GR, XHUG_DEEP_GR, XTG_SURF_GR, XTG_ROOT_GR, XTG_DEEP_GR, CFILE_WG_GR and CFILE_TG_GR has priority on the data in CFILE_GR file.
CTYPE / CTYPEPGD : type of the CFILE_GR / CFILEPGD_GR files, if the latter is provided. CTYPE / CTYPEPGD must then be given. The following values are currently usable:
- “MESONH”: the file type is a MESONH file.
- “GRIB “: the file type is a GRIB file, coming from any of these models: * “ECMWF “: european center forecast model * “ARPEGE”: Arpege french forecast model * “AROME”: Arome french forecast local model * “MOCAGE”: Mocage french research chemistry model
- “ASCII “: PREP/PGD Surfex ASCII file
- “LFI “: PREP/PGD Surfex LFI file

NAM_PREP_TEB_SNOW

NAM_PREP_TEB_SNOW content
Fortran name	Fortran type	Default value
CSNOW_ROOF	CHARACTER(LEN=6)	‘1-L’
CSNOW_ROAD	CHARACTER(LEN=6)	‘1-L’
CFILE_SNOW_TEB	CHARACTER(LEN=28)	‘’
CTYPE_SNOW	CHARACTER(LEN=6)	‘’

CFILEPGD_SNOW_TEB	CHARACTER(LEN=28)	‘’
CTYPEPGD_SNOW	CHARACTER(LEN=6)	‘’
XWSNOW_ROOF	REAL	none
XWSNOW_ROAD	REAL	none
XTSNOW_ROOF	REAL	none
XTSNOW_ROAD	REAL	none
XLWCSNOW_ROOF	REAL	none
XLWCSNOW_ROAD	REAL	none
XASNOW_ROOF	REAL	none
XASNOW_ROAD	REAL	none
XRSNOW_ROOF	REAL	none
XRSNOW_ROAD	REAL	none
LSNOW_IDEAL_TEB	LOGICAL	F

CSNOW_ROAD : snow scheme used over roads
CSNOW_ROOF : snow scheme used over roofs
CFILE_SNOW_TEB : name of the file used to define the snow profiles. The use of a file or prescribed value of XRSNOW_ROOF/ROAD, XTSNOW_ROOF/ROAD, XWSNOW_ROOF/ROAD and XASNOW_ROOF/ROAD has priority on the data in CFILE_SNOW_TEB file
CTYPE_SNOW : type of the CFILE_SNOW_TEB file, if the latter is provided. CTYPE_SNOW must then be given. The following values are currently usable:
- “MESONH”: the file type is a MESONH file.
- “GRIB “: the file type is a GRIB file, coming from any of these models: * “ECMWF “: european center forecast model * “ARPEGE”: Arpege french forecast model * “AROME”: Arome french forecast local model * “MOCAGE”: Mocage french research chemistry model
- “LFI “: LFI PREP file
- “ASCII”: ASCII PREP FILE
CFILEPGD_SNOW_TEB : name of the associated PGD file if CFILE_SNOW_TEB is a PREP files.
CTYPEPGD_SNOW : type of the CFILEPGD_SNOW
XWSNOW_ROAD : snow reservoir for roads (kg/m$^{2}$)
XWSNOW_ROOF : snow reservoir for roofs (kg/m$^{2}$)
XTSNOW_ROAD : snow temperature for roads (K)
XTSNOW_ROOF : snow temperature for roofs (k)
XLWCSNOW_ROAD : snow liquid water content for roads (kg/m$^{3}$)
XLWCSNOW_ROOF : snow liquid water content for roofs (kg/m$^{3}$)
XRSNOW_ROOF : snow density for roofs (kg/m$^{3}$)
XRSNOW_ROAD : snow density for roads (kg/m$^{3}$)
XASNOW_ROAD : snow albedo for roads (-)
XASNOW_ROOF : snow albedo for roofs (-)
LSNOW_IDEAL_TEB : if LSNOW_IDEAL_TEB = F , only one value can be given for following snow parameters and a vertical interpolation is processed. If LSNOW_IDEAL_TEB = T, values are given for each layer and there is no vertical interpolation performed.

NAM_PREP_ISBA

Warning

This namelist comes from SURFEX 9.0.0 user guide https://www.umr-cnrm.fr/surfex/IMG/pdf/surfex_tecdoc.pdf.

NAM_PREP_ISBA content
Fortran name	Fortran type	Default value
CFILE_ISBA	CHARACTER(LEN=28)
CTYPE	CHARACTER(LEN=6)
CFILEPGD_ISBA	CHARACTER(LEN=28)
CTYPEPGD	CHARACTER(LEN=6)
XHUG_SURF	REAL	‘NONE’
XHUG_ROOT	REAL	‘NONE’
XHUG_DEEP	REAL	‘NONE’
XHUGI_SURF	REAL	‘NONE’
XHUGI_ROOT	REAL	‘NONE’
XHUGI_DEEP	REAL	‘NONE’
CFILE_HUG_SURF	CHARACTER(LEN=28)
CFILE_HUG_ROOT	CHARACTER(LEN=28)
CFILE_HUG_DEEP	CHARACTER(LEN=28)
CFILE_HUG	CHARACTER
CTYPE_HUG	CHARACTER(LEN=6)
XTG_SURF	REAL	‘NONE’
XTG_ROOT	REAL	‘NONE’
XTG_DEEP	REAL	‘NONE’
CFILE_TG_SURF	CHARACTER(LEN=28)
CFILE_TG_ROOT	CHARACTER(LEN=28)
CFILE_TG_DEEP	CHARACTER(LEN=28)
CFILE_TG	CHARACTER(LEN=28)
CTYPE_TG	CHARACTER(LEN=6)
NYEAR	INTEGER	‘NONE’
NMONTH	INTEGER	‘NONE’
NDAY	INTEGER	‘NONE’
XTIME	REAL	‘NONE’
LISBA_CANOPY	LOGICAL	.FALSE.
LEXTRAP_TG	LOGICAL	.FALSE.
LEXTRAP_WG	LOGICAL	.FALSE.
LEXTRAP_WGI	LOGICAL	.FALSE.
LEXTRAP_SN	LOGICAL	.FALSE.

CFILE_ISBA / CFILEPGD_ISBA : name of the PREP / PGD files used to define any ISBA variable. The use of a file or prescribed value XHUG_SURF, XHUG_ROOT, XHUG_DEEP, XTG_SURF, XTG_ROOT, XTG_DEEP, CFILE_WG and CFILE_TG has priority on the data in CFILE_ISBA file.
CTYPE / CTYPEPGD : type of the CFILE_ISBA / CFILEPGD_ISBA files, if the latter is provided. CTYPE / CTYPEPGD must then be given. The following values are currently usable:
- “MESONH”: the file type is a MESONH file.
- “GRIB “: the file type is a GRIB file, coming from any of these models:
  - “ECMWF “: european center forecast model
  - “ARPEGE”: Arpege french forecast model
  - “AROME”: AROME french forecast local model
  - “MOCAGE”: Mocage french research chemistry model
- “ASCII “: PREP/PGD Surfex ASCII file
- “LFI “: PREP/PGD Surfex LFI file
XHUG_SURF : uniform prescribed value of liquid soil water index (SWI) for the surface soil layer. This prescribed value, if defined, has priority on the use of CFILE_HUG and CFILE_ISBA data.
XHUG_ROOT : uniform prescribed value of liquid soil water index (SWI) for the root zone soil layer(s). This prescribed value, if defined, has priority on the use of CFILE_HUG and CFILE_ISBA data.
XHUG_DEEP : uniform prescribed value of liquid soil water index (SWI) for the deep soil layer(s). This prescribed value, if defined, has priority on the use of CFILE_HUG and CFILE_ISBA data.
XHUGI_SURF : uniform prescribed value of ice soil water index (SWI) for the surface soil layer. This prescribed value, if defined, has priority on the use of CFILE_HUG and CFILE_ISBA data.
XHUGI_ROOT : uniform prescribed value of ice soil water index (SWI) for the root zone soil layer(s). This prescribed value, if defined, has priority on the use of CFILE_HUG and CFILE_ISBA data.
XHUGI_DEEP : uniform prescribed value of ice soil water index (SWI) for the deep soil layer(s). This prescribed value, if defined, has priority on the use of CFILE_HUG and CFILE_ISBA data.
CFILE_HUG_SURF : name of the file used to define the liquid soil water index (SWI) for the surface soil layer.
CFILE_HUG_ROOT : name of the file used to define the liquid soil water index (SWI) for the root zone soil layer(s).
CFILE_HUG_DEEP : name of the file used to define the liquid soil water index (SWI) for the deep soil layer(s).
CFILE_HUG : name of the file used to define the soil water profiles.

Warning

The use of a file or prescribed value of XHUG_SURF, XHUG_ROOT and XHUG_DEEP has priority on the data in CFILE_HUG file.

CTYPE_HUG : type of the CFILE_HUG file, if the latter is provided. CTYPE_HUG must then be given. The following values are currently usable:
- “MESONH”: the file type is a MESONH file.
- “GRIB “: the file type is a GRIB file, coming from any of these models:
  - “ECMWF “: european center forecast model
  - “ARPEGE”: Arpege french forecast model
  - “AROME”: AROME french forecast local model
  - “MOCAGE”: Mocage french research chemistry model
- “ASCII / LFI “: PREP file from Surfex
- “ASCLLV”: ASCII latlonval file (one file for each depth)
- “NETCDF”: netcdf standard file (one variable by depth)
XTG_SURF : uniform prescribed value of temperature for the surface soil layer, supposed at an altitude of 0m (mean sea level altitude). The temperature is then modified for each point depending on its altitude, following a uniform vertical gradient of -6.5 K km-1. This prescribed value, if defined, has priority on the use of CFILE_TG and CFILE_ISBA data.
XTG_ROOT : uniform prescribed value of temperature for the root zone soil layer(s), supposed at an altitude of 0m (mean sea level altitude). The temperature is then modified for each point depending on its altitude, following a uniform vertical gradient of -6.5 K km-1. This prescribed value, if defined, has priority on the use of CFILE_TG and CFILE_ISBA data.
XTG_DEEP : uniform prescribed value of temperature for the deep soil layer(s), supposed at an altitude of 0m (mean sea level altitude). The temperature is then modified for each point depending on its altitude, following a uniform vertical gradient of -6.5 K km-1. This prescribed value, if defined, has priority on the use of CFILE_TG and CFILE_ISBA data.
CFILE_TG_SURF : name of the file used to define the surface soil temperature profile.
CFILE_TG_ROOT : name of the file used to define the root zone soil temperature profile.
CFILE_TG_DEEP : name of the file used to define the deep soil temperature profile.
CFILE_TG : name of the file used to define the soil temperature profile. The use of a file or prescribed value of XTG_SURF, XTG_ROOT and XTG_DEEP has priorityon the data in CFILE_TG file.
CTYPE_TG : type of the CFILE_TG file, if the latter is provided. CTYPE_TG must then be given. The following values are currently usable:
- “MESONH”: the file type is a MESONH file.
- “GRIB “: the file type is a GRIB file, coming from any of these models:
  - “ECMWF “: european center forecast model
  - “ARPEGE”: Arpege french forecast model
  - “AROME”: AROME french forecast local model
  - “MOCAGE”: Mocage french research chemistry model
- “ASCII / LFI “: PREP file from Surfex
- “ASCLLV”: ASCII latlonval file (one file for each depth)
- “NETCDF”: netcdf standard file (one variable by depth)
NYEAR : year of surface UTC time. It is used only if no atmospheric file or no surface file is given (in those the date can be read).
NMONTH : month of surface UTC time. It is used only if no atmospheric file or no surface file is given (in those the date can be read).
NDAY : day of surface UTC time. It is used only if no atmospheric file or no surface file is given (in those the date can be read).
XTIME : time from midnight of surface UTC time. It is used only if no atmospheric file or no surface file is given (in those the date can be read). (seconds).
LISBA_CANOPY : activates surface boundary multi layer scheme over vegetation.
LEXTRAP_TG : extrapolate TG points where LSM < 0.5 (buffer only)
LEXTRAP_WG : extrapolate WG points where LSM < 0.5 (buffer only)
LEXTRAP_WGI : extrapolate WGI points where LSM < 0.5 (buffer only)
LEXTRAP_SN : extrapolate SNOW (SWE/depth) points where LSM < 0.5 (buffer only)

NAM_PREP_ISBA_CARBON

NAM_PREP_ISBA_CARBON content
Fortran name	Fortran type	Default value
CRESPSL	CHARACTER(LEN=3)	‘DEF’
LSOILGAS	LOGICAL	F
LRESET_CSOIL	LOGICAL	F

CRESPSL : soil respiration option. Possible values are:
- ‘DEF’: no soil respiration
- ‘N92’: Ecosystem respiration from Norman et al. 1992 (odl ‘DEF’ option befoire V9)
- ‘PRM’: Rivalland 2003
- ‘CNT’: Heterotrophic respiration following CENTURY model from Gibelin et al. 2008
- ‘DIF’: activation of the carbon soil dynamics (discretization of soil carbon) from Morel et al. 2019 (JAMES)
LSOILGAS : activation of the soil gas diffusion module to simulate O2, CO2 and CH4 soil dynamics from Morel et al. 2019 (JAMES). !!! this scheme is actually a prototype !!!
LRESET_CSOIL : Flag to initialize isba physic but not soil carbon

NAM_PREP_ISBA_SNOW

NAM_PREP_ISBA_SNOW content
Fortran name	Fortran type	Default value
CSNOW	CHARACTER(LEN=3)	‘D95’

NSNOW_LAYER	INTEGER	1
CFILE_SNOW	CHARACTER(LEN=28)	CFILE_ISBA in
		NAM_PREP_ISBA
CTYPE_SNOW	CHARACTER(LEN=6)	CTYPE in
		NAM_PREP_ISBA
CFILEPGD_SNOW	CHARACTER(LEN=28)	CFILEPGD_ISBA in
		NAM_PREP_ISBA
CTYPEPGD_SNOW	CHARACTER(LEN=6)	CTYPEPGD in
		NAM_PREP_ISBA
LSNOW_IDEAL	LOGICAL	F
LSNOW_FRAC_TOT	LOGICAL	F
LSNOW_PREP_PERM	LOGICAL	T
XWSNOW	REAL(20)
XZSNOW	REAL(20)	1.E+20
XTSNOW	REAL(20)	273.16
XLWCSNOW	REAL(20)
XRSNOW	REAL(20)
XASNOW	REAL	0.5
XSG1SNOW	REAL(20)	none
XSG2SNOW	REAL(20)	none
XHISTSNOW	REAL(20)	none
XAGESNOW	REAL(20)	none
LSWEMAX	LOGICAL	F
XSWEMAX	REAL
NIMPUR	INTEGER	0

CSNOW : type of snow scheme. Possible snow schemes are:
- ‘D95’: Douville et al (1995) snow scheme.
- ‘3-L’: Boone and Etchevers (2001); Decharme et al. (2016) N-layer (default 12) snow scheme
- ‘EBA’: Bogatchev and Bazile (2005), Arpege operational snow scheme.
- ‘CRO’: Crocus model
NSNOW_LAYER : number of snow layers
CFILE_SNOW : name of the file used to define the snow profiles. The use of a file or prescribed value of XRSNOW, XTSNOW, XWSNOW and XASNOW (and XSG1SNOW, XSG2SNOW, XHISTSNOW and XAGESNOW in case of CSNOW = CROCUS) has priority on the data in CFILE_SNOW file
CTYPE_SNOW : type of the CFILE_SNOW file, if the latter is provided. CTYPE_SNOW must then be given. The following values are currently usable:
- “MESONH”: the file type is a MESONH file.
- “GRIB “: the file type is a GRIB file, coming from any of these models: * “ECMWF “: european center forecast model * “ARPEGE”: Arpege french forecast model * “AROME”: AROME french forecast local model * “MOCAGE”: Mocage french research chemistry model
- “LFI “: LFI PREP file
- “ASCII”: ASCII PREP FILE
CFILEPGD_SNOW : name of the associated PGD file if CFILE_SNOW is a PREP files.
CTYPEPGD_SNOW : type of the CFILEPGD_SNOW file, if the latter is provided. CTYPEPGD_SNOW must then be given. The following values are currently usable:
- “MESONH”: the file type is a MESONH file.
- “LFI “: LFI PREP file
- “ASCII”: ASCII PREP FILE
LSNOW_IDEAL : if LSNOW_IDEAL = F , only one value can be given for following snow parameters and a vertical interpolation is processed. If LSNOW_IDEAL = T, values are given for each layer and there is no vertical interpolation performed.
LSNOW_FRAC_TOT : if LSNOW_FRAC_TOT = T, the total snow fraction XPSN = MIN(1.0, ZSNOWSWE(:)/XWCRN_EXPL) where ZSNOWSWE is the snow liquid water content, and XWCRN_EXPL is the critical value of the equivalent water content of the snow reservoir.
LSNOW_PREP_PERM : activates or disactivates initialization over permanent ice areas.
XWSNOW : uniform value to initialize snow content, one for each layer
XZSNOW : depth of snow layers (m). Alternative to XWSNOW.
XTSNOW : uniform value to initialize snow temperature, one for each layer
XLWCSNOW : snow liquid water content (kg/m3)
XRSNOW : uniform value to initialize snow density, one for each layer
XASNOW : uniform value to initialize snow albedo
XSG1SNOW : uniform value to initialize snow layers grain feature 1 for Crocus, one for each layer
XSG2SNOW : uniform value to initialize snow layers grain feature 2 for Crocus, one for each layer
XHISTSNOW : uniform value to initialize snow layer grain historical parameter for Crocus, one for each layer
XAGESNOW : uniform value to initialize snow grain age for Crocus, one for each layer
LSWEMAX : logical switch to set an upper limit on initial snow water equivalent
XSWEMAX : upper limit of initial snow water equivalent
NIMPUR : number of impurity you want to use in your simulation. NIMPUR=1 with black carbon only and NIMPUR=2 with black carbon and dust (to run a simulation with dust only you can set NIMPUR=2 and prescribe no black carbon deposition)

Free-format part

The free-format part of the PRE_IDEA1.nam namelist corresponds to blocks of text that are necessary when using certain options, which are described in the following sections. Each section of the free format part must be introduced by its corresponding keyword (writen on a separated line).

Vertical grid

If you want to define your own vertical grid (CZGRID_TYPE = “MANUAL” in NAM_VER_GRID inside PRE_IDEA1.nam), you must give the heights of the vertical velocity levels from the surface to the top of the domain (NKMAX + 1 values). The corresponding keyword is ZHAT.

Note

Example of free part of PRE_IDEA1.nam :

&NAM_VER_GRID NKMAX       = 10,
              YZGRID_TYPE = 'MANUAL' /

ZHAT












Radiosounding case

For idealized simulations, it is possible to use a radiosounding to initialize the model (CIDEAL = “RSOU” in NAM_CONF_PRE in PRE_IDEA1.nam). The corresponding keyword is RSOU and the free format part has to be set in this order:

RSOU
year month day time

Note

year has to be an integer (example : 1994)
month has to be an integer (example : 4)
day has to be an integer (example : 22)
time has to be a real in seconds (example : 36000. for 10 h)

kind of radiosounding

Note

Nine kind of radiosounding are possible : “STANDARD”, “PUVTHVMR”, “PUVTHVHU”, “ZUVTHVMR”, “ZUVTHVHU”, “PUVTHDMR”, “PUVTHDHU”, “ZUVTHDMR”, “ZUVTHLMR”. Except for the “STANDARD”:

the first letter represents the kind of altitude variable (P for pressure (Pa) and Z for height (m)),
the second and third letters represent the kind of wind variables (U for zonal wind (m/s), V for meridian wind (m/s)),
the fourth, fifth and sixth letters represent the kind of temperature variable (THV for virtual potential temperature (K), THD for dry potential temperature (K) and THL for liquid potential temperature (K)),
the seventh and eighth letters represent the kind of moist variable (HU for relative humidity (%) and MR for vapor mixing ratio (kg/kg)).

In case of “STANDARD”, the altitude variable is the pressure (Pa), the wind variables are direction and wind force (m/s), the temperature variable is the temperature (K) and the moist variable is the dew point temperature (K).

height of the ground level (real, in meters)
pressure at the ground level (real,in Pascal)
temperature at ground level (real, in Kelvin)
humidity at ground level (real, unit depends on the kind radiosounding)

number of wind data levels (integer)
altitude_level_1 first_wind_variable second_wind_variable
altitude_level_2 first_wind_variable second_wind_variable
...
altitude_level_top first_wind_variable second_wind_variable

Note

Units of altitude_levels, first_wind_variable and second_wind_variable depends on the kind radiosounding you choose.

number of mass data levels (integer)
altitude_level_2 temperature humidity additional_cloud_variable(s)
altitude_level_3 temperature humidity additional_cloud_variable(s)
...
altitude_level_top temperature humidity additional_cloud_variable(s)

Note

Number of mass data levels includes the ground level (i.e. the first level). That is why the following list starts at level 2.
Units of altitude_levels, temperature, humidity and additional_cloud_variable(s) depends on the kind radiosounding you choose).

Warning

You should make sure that the highest level of the radiosounding is located above the highest vertical level of the model.
Additional cloud variables. For the moment, this configuration works only for kind=”PUVTHDMR” or “ZUVTHDMR” and L1D=.TRUE.. It is planned to compute radiation diagnostics with the DIAG program :
- cloud mixing ratio if LUSERC=T or LUSERI=T (real, in kg/kg)
- ice mixing ratio if LUSERI=T (real, in kg/kg)

Note

Example of free part of PRE_IDEA1.nam :

RSOU
1990 10 3 72000.
'STANDARD'
200.
100240.
287.5
276.
2
85000. 20. 10.
70000. 30. 10.
3
90000. 280. 275.
60000. 271. 269.

Constant moist Brunt-Vaisala case

For idealized simulations, it is possible to use a constant Brunt-Vaisala frequency, shear and humidity layers to initialize the model (CIDEAL = “CSTN” in NAM_CONF_PRE in PRE_IDEA1.nam). The corresponding keyword is CSTN and the free format part has to be set in this order:

CSTN
year month day time

Note

year has to be an integer (example : 1994)
month has to be an integer (example : 4)
day has to be an integer (example : 22)
time has to be a real in seconds (example : 36000. for 10 h)

number of levels (integer)
virtual potential temperature at ground level (real, in K)
pressure at ground level (real, in Pa)
height at all levels.
zonal wind component at all levels
meridian wind component at all levels
relative humidity at all levels
moist Brunt Vaisala frequency at all layers

Note

The first level is the ground level.
For moist Brunt Vaisala frequency, the number of layers is the number of levels - 1.

In this case, the level number can even be equal to 1, because the profile information is linearly interpolated on the model grid without orography (wind components, $\theta_v$ and humidity) before the application of the Laplace relation to deduce the pressure and the vapor mixing ratio. Thus, the layers’ thicknesses are never too large to invalidate the Laplace relation.

Note

Example of free part of PRE_IDEA1.nam :

CSTN
06 06 21600.
5
5
100240.
1000. 1500. 3000. 4000.
20. 25. 30. 35.
10. 12.5 11.5 15.
84. 85. 79. 87.
01 0.014 0.015 0.016

The forced version

For idealized simulations a forced mode can be useful to impose the effects of a simplified large scale environment to the model solution. This functionality works when LFORCING=.TRUE. and CIDEAL=’RSOU’ or ‘CSTN’ in NAM_CONF_PRE and when LGEOSBAL =.FALSE. in NAM_VPROF_PRE for inclusion of a geostrophic wind forcing in PRE_IDEA1.nam. All forcing fields are issued from spatial interpolation of chronological series of 1D data (provided by the user onto the model grid).

The forcing fields can be time dependent. Application of the forcing begins as soon as the date and time of the first set of forcing field given by the user, is lower or equal to the current date and time of the model run. The forcing action of the last forcing field is remanant, this is a way to impose a stationnary forcing. When the current date and time of the model run is bounded by two successive forcing fields, a simple linear interpolation in time is made. Note that an available Newtonian relaxation forcing type on [u, v] and/or [$\theta$, rv] is exclusive from the other physical forcings.

The forcing information and soundings have to be added at the end of the free-format part already written for CIDEAL=”CSTN” or “RSOU”. Depending of the altitude of the forcing data two keywords are available, ZFRC means that the altitude of the forcing data are in height scale (m) and PFRC means that the altitude of the forcing data are in pressure scale (Pa).

The free format part has to be set in this order:

ZFRC or PFRC
number of time dependent forcing (integer)
year month day time

Note

year has to be an integer (example : 1994)
month has to be an integer (example : 4)
day has to be an integer (example : 22)
time has to be a real in seconds (example : 36000. for 10 h)

Warning

The 1D forcing data are different from the one used to initialize the model because specific data have to be entered. The data used to define each forcing are given sequentially in the following order (one item per line):

ground height (real, m)
ground pressure (real, Pa)
potential temperature at ground level (real, K)

Note

It is used later in the code to compute - if asked - a time varying sea surface temperature).

humidity (real, kg/kg) at ground level
number of level (integer)
altitude of level 1 list of forcing fields
idem at level 2
...
idem at level top

Note

Unit of altitudes depends on the keyword chosen ZFRC for (m) and PFRC for (Pa)
List of forcing field is :
- $u_{frc}$ component at the corresponding level (real, m/s)
- $v_{frc}$ component at the corresponding level (real, m/s)
- $\theta_{frc}$ at the corresponding level (real, K)
- $rv_{frc}$ at the corresponding level (real, kg/kg)
- $w_{frc}$ at the corresponding level (real, m/s)
- $(\partial\theta / \partial t)_{frc}$ at the corresponding level (real, K/s)
- $(\partial rv / \partial t)_{frc}$ at the corresponding level (real, 1/s)
- $(\partial u / \partial t)_{frc}$ at the corresponding level (real, m/s2)
- $(\partial v / \partial t)_{frc}$ at the corresponding level (real, m/s2).

If PFRC is the forcing type, an additional sounding is given in order to convert the pressure levels into height levels with enough accuracy. Data are organized as follows:

number of level (integer)
pressure at level 1 (real, Pa), $\theta$ at level 1 (real, K) and rv at level 1 (real, kg/kg).

This operation is repeated until the previous number of sounding is reached.

Note

Example of free part of PRE_IDEA1.nam :

ZFRC
1
1983 07 01 0.
0
1000000
284.5
0.008
6
   5. -7.0 0.0 281.10 0.00540 -0.00000 0. 0. 0. 0.
  15. -7.0 0.0 281.10 0.00540 -0.00000 0. 0. 0. 0.
1095. -7.0 0.0 280.75 0.00540 -0.00300 0. 0. 0. 0.
1145. -7.0 0.0 290.60 0.00190 -0.00300 0. 0. 0. 0.
3000. -7.0 0.0 304.15 0.00190 -0.00300 0. 0. 0. 0.
9000. -7.0 0.0 346.15 0.00190 -0.00300 0. 0. 0. 0.

The advective forcing

For 2D idealized simulation, an advective forcing can be used to impose effects to the model solution. This functionality works when L2D_ADV_FRC=.TRUE. in NAM_CONF_PRE in PRE_IDEA1.nam. The advecting forcings mimic the latidudinal humidity and temperature advection not taken into account in a 2D model. The forcing information and soundings have to be added at the end of the free-format part already written for CIDEAL=”CSTN” or “RSOU”. The corresponding keyword is ZFRC_ADV and the free format part has to be set in this order:

ZFRC_ADV
number of forcing files
type of forcing : ZADV2D for Z levels or PADV2D for pressure levels
number vertical levels for the file
date of first forcing : year (integer) month (integer) day (integer) time (real, seconds)
name of the file with horizontal mean profile of theta, rv
name of the advective forcing file

Note

Example of free part of PRE_IDEA1.nam :

ZFRC_ADV
1
ZADV2D
52
1997 07 15 00000.
"mean_atm_07.dat"
"frc_ideal_7_70km.dat"

The relaxation forcing

For 2D idealized simulation, a relaxation forcing can be used to impose effects to the model solution. This functionality works when L2D_REL_FRC=.TRUE. in NAM_CONF_PRE in PRE_IDEA1.nam. The relaxation forcing allows the relax the model fields towards a 2D climatology for temperature and humidity. The forcing information and soundings have to be added at the end of the free-format part already written for CIDEAL=”CSTN” or “RSOU”. The corresponding keyword is ZFRC_REL and the free format part has to be set in this order:

ZFRC_REL
number of forcing files
type of forcing : ZREL2D for Z levels or PREL2D for pressure levels
number vertical levels for the file
date of first forcing : year (integer) month (integer) day (integer) time (real, seconds)
name of the file with horizontal mean profile of theta, rv
name of the advective forcing file

Note

Example of free part of PRE_IDEA1.nam :

ZFRC_REL
1
ZREL2D
52
1997 07 15 00000.
"mean_atm_07.dat"
"frc_ideal_7_70km.dat"

Discretized orography

You can prescribe your own orography when CZS = “DATA” in NAM_CONF_PRE in PRE_IDEA1.nam. Only the orography corresponding to the computational domain must be provided in the free format part. For 3D orography, data are read like if it was a map (the first line is the Northern border and the first data is the North-West corner) with one line per Y-axis increment. The corresponding keyword is ZSDATA.

Note

Example of free part of PRE_IDEA1.nam :

ZSDATA
30. 35. 50. 30. 30.
59.5 133.3 100.2 136.7 100.
89.5 183.3 200.2 299.7 170.5
112.5 193.0 210.2 206.7 120.
82.5 153.0 180.5 156.7 100.3

The ocean version

For oceanic version of Meso-NH, the initial and forcing profiles of the ocean are written in the free-format part of PRE_IDEA1.nam file, where the altitude variable is the depth. To follow usual convention for ocean data, the 1D profiles are given starting from the surface (positive value). The corresponding keyword is RSOU and the profile and forcing information are written in the following order :

RSOU
year month day time

Note

year has to be an integer (example : 1994)
month has to be an integer (example : 4)
day has to be an integer (example : 22)
time has to be a real in seconds (example : 36000. for 10 h)

kind of data used for the profile

Note

Two kind are possible :

KIND=”IDEALOCE” : data written in PRE_IDEA1.nam
KIND=”STANDOCE” : data written in a NetCDF file (See set_rsou.f90 for details)

The following format is valid for KIND=”IDEALOCE” :

atmospheric pressure at the surface (real, in Pascal)

Warning

The surface corresponds to the top domain of the oceanic model.

sea surface temperature at the surface (real, in K)
sea surface salinity at the surface (real, in g/kg)

number of sea current levels (integer)
level 1 : depth (real, m)  u-current (real, m/s)  v-current (real, m/s)
level 2 : depth (real, m)  u-current (real, m/s)  v-current (real, m/s)
level ...

number of mass data levels (integer).
level 2 : depth (real, m)  temperature (real, K), salinity (real, g/kg)
level 3 : depth (real, m)  temperature (real, K), salinity (real, g/kg)
level ...

number of time-varying forcing (integer).

Warning

The data used to define each forcing are given sequentially in the following order (one item per line):

year month day time

Note

year has to be an integer (example : 1994)
month has to be an integer (example : 4)
day has to be an integer (example : 22)
time has to be a real in seconds (example : 36000. for 10 h)

u-stress (real, m2/s2)
v-stress (real, m2/s2)
heat turbulent flux (real, W/m2)
radiative flux (real, W/m2)

Note

Surface fluxes are positive when going upward (from the ocean to the atmosphere).
In the KIND=”STANDOCE” initial 1D profiles and surface fluxes(t) are read from 2 netcdf files (See set_rsou.f90 for details)