PREP_REAL_CASE

PREP_REAL_CASE program performs the change of orography and vertical grid by interpolating horizontally and vertically for a GRIB file or only vertically for a Meso-NH file. The Meso-NH output file will be used either for the beginning of the simulation or for coupling. The main hypothesis is that hydrostatism is verified. Therefore, if the input file is a Meso-NH file, there is a small loss of information.

PREP_REAL_CASE program and its corresponding namelist and function

Executable

Namelist

Function

PREP_REAL_CASE

PRE_REAL1.nam

Prepare atmospheric initial and boundary file

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

NAM_AERO_CONF

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

NAM_CH_CONF content

Fortran name

Fortran type

Default value

LUSECHAQ

LOGICAL

.FALSE.

LUSECHIC

LOGICAL

.FALSE.

LUSECHEM

LOGICAL

.FALSE.

  • LUSECHAQ : flag to activate aqueous phase chemistry.

  • LUSECHIC : flag to activate ice phase chemistry.

  • LUSECHEM : flag to activate gazeous phase chemistry whithout aqueous phase chemistry.

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

NAM_FILE_NAMES content

Fortran name

Fortran type

Default value

HATMFILE

CHARACTER(LEN=128)

‘’

HATMFILETYPE

CHARACTER(LEN=6)

‘MESONH’

HPGDFILE

CHARACTER(LEN=128)

‘’

HSURFFILE

CHARACTER(LEN=128)

‘’

HSURFFILETYPE

CHARACTER(LEN=6)

‘MESONH’

HCHEMFILE

CHARACTER(LEN=128)

‘’

HCHEMFILETYPE

CHARACTER (LEN=6)

‘MESONH’

HCAMSFILE

CHARACTER(LEN=128)

‘’

HCAMSFILETYPE

CHARACTER(LEN=6)

‘NETCDF’

CINIFILE

CHARACTER(LEN=128)

‘INIFILE’

  • HATMFILE : name of the atmospheric file.

  • HATMFILETYPE : type of atmospheric file (‘GRIBEX’, ‘MESONH’)

  • HPGDFILE : name of the Physiographic Data File.

  • HSURFFILE : optional name of the file containing the surface fields.

  • HSURFFILETYPE : type of surface file (‘GRIBEX’, ‘MESONH’)

  • HCHEMFILE : optional name of the file containing the chemical species if they are not in the HATMFILE or if the ones of the HATMFILE are not used (only if HATMFILETYPE is ‘GRIBEX’). The grids must be the same as the ones of the output file (CINIFILE).

  • HCHEMFILETYPE : type of the chemical file (‘GRIBEX’, ‘MESONH’, ‘MOZART’,’CAMSEU’). MOZART and CAMSEU must be in netcdf format

  • HCAMSFILE : optional name of the CAMS file containing the aerosols species. The file must be a NETCDF.

  • HCAMSFILETYPE : type of the CAMS file (‘NETCDF’)

  • CINIFILE : name of the Meso-NH file, used as initial or coupling file in a Meso-NH simulation

NAM_HURR_CONF

The purpose of this namelist is to replace an existing cyclone (hurricane filtering) with another one whose characteristics are defined (vortex bogussing). Each step (hurricane filtering and vortex bogussing) has to be done separately with the PREP_REAL_CASE program.

In a mono-model configuration, the first PREP_REAL_CASE job allows to remove analysed hurricane from the input GRIB fields: filtered and interpolated fields are written in a Meso-NH file (hurricane filtering). It is used as input file for the second PREP_REAL_CASE job during which the analytical vortex is added (vortex bogussing).

For a grid-nesting simulation, the hurricane filtering is first applied for the outer domain (dad model), with the program PREP_REAL_CASE. The filtered fields are then horizontally interpolated for inner domains with the program SPAWNING. Then, for each inner domain, a vortex bogussing is added with the program PREP_REAL_CASE.

Note

  • The hurricane filtering can be applied on five atmospheric GRIB fields (HATMFILETYPE=’GRIBEX’ in NAM_FILE_NAMES in PRE_REAL1.nam) : the two horizontal components of wind, the absolute temperature, the humidity and the surface pressure reduced to ground level. Each output field is decomposed into three parts: first, the BASic part is computed by the low-pass Barnes filter; then the hurricane (symmetric) disturbance is computed from the remainder disturbance part. The initial fields are then remplaced by their ENVironmental part: total field minus hurricane disturbance part.

  • The vortex bogussing consists on a symmetric vortex added to the input atmospheric Meso-NH fields (HATMFILETYPE=’MESONH’ in NAM_FILE_NAMES in PRE_REAL1.nam). The tangential wind is computed from an analytical formulation [Holland, 1980]: Mercator projection must be used to respect hypotheses of Holland. Then, the balanced mass field is deduced from the thermal wind relation. The bogus of the two horizontal components of wind and the potential temperature is added to the initial (filtered) fields.

Warning

For hurricane filtering and vortex bogussing PREP_REAL_CASE must be launch on 1 core. It’s not yet parallelized.

NAM_HURR_CONF content

Fortran name

Fortran type

Default value

LFILTERING

LOGICAL

.FALSE.

CFILTERING

CHARACTER(LEN=5)

‘UVT’

NK

INTEGER

50

XLAMBDA

REAL

0.2

XLATGUESS

REAL

XUNDEF

XLONGUESS

REAL

XUNDEF

XBOXWIND

REAL

XUNDEF

XRADGUESS

REAL

XUNDEF

NPHIL

INTEGER

24

NDIAG_FILT

INTEGER

-1

NLEVELR0

INTEGER

15

LBOGUSSING

LOGICAL

.FALSE.

XLATBOG

REAL

XUNDEF

XLONBOG

REAL

XUNDEF

XVTMAXSURF

REAL

XUNDEF

XRADWINDSURF

REAL

XUNDEF

XMAX

REAL

16000.0

XC

REAL

0.7

XRHO_Z

REAL

-0.3

XRHO_ZZ

REAL

0.9

XB_0

REAL

1.65

XBETA_Z

REAL

-0.5

XBETA_ZZ

REAL

0.35

XANGCONV0

REAL

0.0

XANGCONV1000

REAL

0.0

XANGCONV2000

REAL

0.0

CDADATMFILE

CHARACTER(LEN=128)

CDADBOGFILE

CHARACTER(LEN=128)

  • LFILTERING : Flag to filter the fields (U,V,T,Q,reduced Ps) of the atmospheric file

  • CFILTERING : to choose the fields to be filtered (U,V,T,reduced Ps).

    • ‘UVT’ : U,V,T are filtered (default),

    • ‘UVTQ’ : U,V,T,Q are filtered (default),

    • ‘UVTP’ : U,V,T and reduced PS are filtered,

    • ‘UVTPQ’ : U,V,T,Q and reduced PS are filtered,

  • NK : number of points of the half-width of the window in which the Barnes filter is applied to compute low-pass component of a given field

  • XLAMBDA : a coefficient in the exponential weighting function of the Barnes filter

  • XLATGUESS : latitude of the guessed position of the cyclone center

  • XLONGUESS : longitude of the guessed position of the cyclone center

  • XBOXWIND : half-width of the box inside which the dynamical center is searched from the guessed position (km)

  • XRADGUESS : guess of the radius of the domain in which the cyclone will be filtered (km)

  • NPHIL : number of azimuthal directions used for the cylindrical coordinates

  • NDIAG_FILT : allows storage of several components calculated from total fields. Be careful, the components are on the GRIB vertical grid. Then to visualize all the GRIB vertical levels, the number of MesoNH vertical levels must be equal or greater than the number of levels in the input GRIB file.

    • 0 : total (unfiltered) fields: UT15, VT15 for wind components; TEMPTOT, PRESTOT for absolute temperature and surface pressure, environmental (filtered) fields (total field minus hurricane disturbance component): UT16, VT16, TEMPENV, PRESENV,

    • 0,1 : basic fields (low-pass component isolated by the Barnes filter): UT17, VT17, TEMPBAS, PRESBAS,

    • 0,1,2 : total disturbance tangential wind component (XVTDIS).

  • NLEVELR0 : level used to compute R0

  • LBOGUSSING : Flag to switch on the addition of the bogus vortex (logical)

  • XLATBOG : latitude of the bogussed position of the analytical cyclone center

  • XLONBOG : longitude of the bogussed position of the analytical cyclone center

  • XVTMAXSURF : maximum tangential wind near the surface or about 500 m altitude (m/s)

  • XRADWINDSURF : radius of maximum wind near the surface or about 500 m altitude (km)

  • XMAX : altitude where the tangential wind vanishes (m)

  • The following variables are parameters describing tangential wind in Holland’s law (see formulation in routine holland_vt.f90).

    • XC : standard coefficient for maximum tangential wind,

    • XRHO_Z, XRHO_Z : standard coefficients for radius of maximum wind,

    • XB_0, XBETA_Z, XBETA_ZZ : standard coefficients for B parameter.

  • XANGCONV0, XANGCONV1000, XANGCONV2000 : convergence angle of wind near the surface, at 1000m and 2000m altitude.

  • CDADATMFILE : if LBOGUSSING=.TRUE. : name of the dad of HATMFILE

  • CDADBOGFILE : if LBOGUSSING=.TRUE. : name of the dad of CINIFILE. The program will check that CDADATMFILE and CDADBOGFILE have the same characteristics, before replacing the dad name of CINIFILE by CDADBOGFILE instead of CDADATMFILE. CDADBOGFILE must exist before running the PREP_REAL_CASE job.

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.

NAM_REAL_CONF

NAM_REAL_CONF content

Fortran name

Fortran type

Default value

CEQNSYS

CHARACTER(LEN=3)

‘DUR’

CPRESOPT

CHARACTER(LEN=5)

‘CRESI’

NVERB

INTEGER

1

LSHIFT

LOGICAL

.TRUE./.FALSE.

LDUMMY_REAL

LOGICAL

.FALSE.

LRES

LOGICAL

.FALSE.

XRES

REAL

1.E-07

NITR

INTEGER

4

LCOUPLING

LOGICAL

.FALSE.

NHALO

INTEGER

1

JPHEXT

INTEGER

1

NVERLEVEL_GRIB

INTEGER

120

  • CEQNSYS : EQuatioN SYStem

    • ‘LHE’: Lipps-HEmler 1982

    • ‘MAE’: Modified Anelastic Equations

    • ‘DUR’: following DURran 1990 derivations

Note

By default, if if HATMFILETYPE = ‘GRIBEX’, CEQNSYS=’DUR’ and if HATMFILETYPE = ‘MESONH’, CEQNSYS is fixed to the one in input Meso-NH file.

  • CPRESOPT : option for pressure solver (‘RICHA’, ‘CGRAD’, ‘CRESI’, ‘ZRESI’).

  • NVERB : verbosity level (error diagnostics are computed if NVERB>4)

  • LSHIFT : flag to shift altitudes in boundary layer

Note

By default, if if HATMFILETYPE = ‘GRIBEX’, LSHIFT=.TRUE. and if HATMFILETYPE = ‘MESONH’, LSHIFT is fixed to the one in input Meso-NH file.

  • LDUMMY_REAL : flag to read dummy fields stored in the GRIB file asking for additional fields by modifying FULLPOS namelist or {bf extractecmwf} by modifying MARS requests).You have to fill also a free formatted part as described in the exemple of AROME field.

  • LRES : flag to change the residual divergence limit.

  • XRES : Value of the residual divergence limit

  • NITR : number of iterations used for the resolution of the elliptic equation (solver = “CPRESOPT”).

  • LCOUPLING : Flag to handle the output file as coupling file instead of initial file : no surface field will be computed. (time saving and smaller file) }

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

  • NVERLEVEL_GRIB : number of vertical levels of ICON grib data files used as init or coupling (120 for ICON-global; 74 for ICON-EU)

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

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 “FLAKE” sea scheme temperature.

NAM_PREP_FLAKE content

Fortran name

Fortran type

Default value

XTS_UNIF

REAL

XUNIF_T_SNOW

REAL

MIN(273.15,XTS_WATER)

XUNIF_T_ICE

REAL

MIN(273.15,XTS_WATER)

XUNIF_T_WML

REAL

MIN(273.15,XTS_WATER)

XUNIF_T_BOT

REAL

TS_WATER or 277.15 if TS_WATER \(\leq\) 273.15

XUNIF_T_B1

REAL

TS_WATER-0.1 or 277.05 if TS_WATER \(\leq\) 273.15

XUNIF_CT

REAL

0.5

XUNIF_H_SNOW

REAL

0.0

XUNIF_H_ICE

REAL

0.0 or 0.01 if XTS_WATER XTS_WATER \(\leq\) 273.15

XUNIF_H_ML

REAL

XWATER_DEPTH or XWATER_DEPTH/2 if TS_WATER \(\leq\) 273.15

XUNIF_H_B1

REAL

0.0

CFILE_FLAKE

CHARACTER(LEN=28)

CFILE in NAM_PREP_SURF_ATM

CTYPE

CHARACTER(LEN=6)

CFILETYPE in NAM_PREP_SURF_ATM

CFILEPGD_FLAKE

CHARACTER(LEN=28)

CFILEPGD in NAM_PREP_SURF_ATM

CTYPEPGD

CHARACTER(LEN=6)

CFILEPGDTYPE in NAM_PREP_SURF_ATM

LCLIM_LAKE

LOGICAL

.FALSE.

NYEAR

INTEGER

NMONTH

INTEGER

NDAY

INTEGER

XTIME

REAL

LWAT_SBL

LOGICAL

.FALSE.

  • XTS_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_FLAKE data.

  • XUNIF_T_SNOW : surface temperature of snow (K)

  • XUNIF_T_ICE : surface temperature at the ice-atmosphere or at the ice-snow interface (K)

  • XUNIF_T_WML: mixed-layer temperature (K)

  • XUNIF_T_BOT: water temperature at the bottom of the lake (K)

  • XUNIF_T_B1: temperature at the bottom of the upper layer of sediments (K)

  • XUNIF_CT: shape factor (thermocline)

  • XUNIF_H_SNOW: snow layer thickness (m)

  • XUNIF_H_ICE: ice layer thickness (m)

  • XUNIF_H_ML: thickness of the mixed-layer (m)

  • XUNIF_H_B1: thickness of the upper level of the active sediments (m)

  • CFILE_FLAKE / CFILEPGD_FLAKE : name of the PREP and 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_FLAKE file.

  • CTYPE / CTYPEPGD : type of the CFILE_FLAKE / CFILEPGD_FLAKE 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”: Surfex PREP / PGD ASCII file

    • “LFI “: Surfex PREP / PGD LFI file

  • LCLIM_LAKE : to use the climatological lake database to initialise FLAKE pronostic variables. Needs to link with LAKE_LTA_NEW.nc.

  • 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).

  • LWAT_SBL : activates surface boundary multi layer scheme over inland water.

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_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).

Free-format part

Vertical grid

If you want to define your own vertical grid (CZGRID_TYPE = ‘MANUAL’ in NAM_VER_GRID), you must give the heights of the vertical velocity levels. You must NKMAX + 1 values from the surface to the top of the domain.

Note

Example of free part of PRE_REAL1.nam :

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

ZHAT
0.
1050.
2100.
3250.
4300.
5200.
6100.
7000.
8000.
9000.
10000.

Chemical species

This part is only used if you have previously run extractarpege with MOCAGE outputs. This part must begin by the keyword MOC2MESONH. Then, the list of the Meso-NH chemical species, and their corresponding GRIB code in the GRIB file, is specified as follows:

Note

Example of free part of PRE_REAL1.nam :

MOC2MESONH
transfer mocage/RACM variables (default values)
2 # NUMBER OF OPTIONAL GRIB VARIABLES
(A4,1X,I5)
O3   180
NO2  183

Read extra AROME variables

AROME GRIB files obtained with Extract Meteo-France data contain fields which aren’t read by default by PREP_REAL_CASE program. To use these fields, you have to use LDUMMY_REAL=T in NAM_REAL_CONF and fill a free formatted part at the end of PRE_REAL1.nam. This part must begin by the keyword DUMMY_2D (for 2D fields). Then, the list of the Meso-NH dummy fields, and their corresponding GRIB code in the GRIB file (parameter code, and if ambiguity with another variables type and value of level), is specified as follows :

DUMMY_2D
diagnostics
10
ACCPLUIE 62
ACCNEIGE 79
ACCGRAUPEL 78
INSPLUIE 169
INSNEIGE 64
INSGRAUPEL 63
UM05 33 105 10
VM05 34 105 10
CLSHUMI.RELA 52 105 2
CLSTEMPE 11 105 2