DIAG

After running the model, useful quantities can be diagnosed from prognostic variables contained in the synchronous backup files. It is done by the program DIAG which computes diagnostic variables.

DIAG program and its corresponding namelist and function

Executable	Namelist	Function	Output
DIAG	DIAG1.nam	Compute diagnostics after simulation	DIAGFILE.{des,nc}

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

• NAM_CONFIO

• NAM_DIAG - Default

• NAM_DIAG - General

• NAM_DIAG - Deep convection

• NAM_DIAG - Shallow convection

• NAM_DIAG - Microphysics

• NAM_DIAG - Turbulence

• NAM_DIAG - Radiation

• NAM_DIAG - Aerosols

• NAM_DIAG - Chemistry

• NAM_DIAG - Lightnings

• NAM_DIAG - Lagrangian tracers

• NAM_DIAG - GPS simulator

• NAM_DIAG - Satellite simulator

• NAM_DIAG - RADAR simulator

• NAM_DIAG - LIDAR simulator

• NAM_DIAG - Interpolation levels

• NAM_DIAG - Clustering

• NAM_DIAG - Coarse graining

• NAM_DIAG_BLANK

• NAM_DIAG_FILE

• NAM_STO_FILE

• NAM_CONF_DIAG

• NAM_PARAM_KAFRn

• NAM_PARAM_RADn

• NAM_DIAG_SURF_ATMn

• NAM_WRITE_DIAG_SURFn

• NAM_DIAG_SURFn

• NAM_DIAG_ISBAn

• NAM_DIAG_TEBn

• NAM_DIAG_FLAKEn

• NAM_DIAG_OCEANn

Note

For additional SURFEX namelists, please go to SURFEX documentation https://www.umr-cnrm.fr/surfex/IMG/pdf/surfex_tecdoc.pdf.

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_DIAG - Default

By default, the following variables are written in the output file created by the DIAG program:

Name	Meaning [Unit]	Dimension
ZS	orography [m]	2D
ZSMT	smoothed orography for SLEVE vertical coordinate [m]	2D
RHODREF	dry density for reference state with orography [\(kg/m^3\)]	3D
THVREF	\(\theta_v\) for reference state with orography [K]	3D
RHOREFZ	\(\rho_d(z)\) for reference state without orography [\(kg/m^3\)]	1D
THVREFZ	\(\theta_v(z)\) for reference state with orography [K]	1D
EXNTOP	exner function at model top [-]
SVPPn	passive pollutant n concentration only if LPASPOL=T in YINIFILE.des [\(g/m^3\)]	3D

Surface variables by default only if CSURF='EXTE' in YINIFILE.des:

Name	Meaning [Unit]	Dimension
UM10, VM10	zonal and meridional wind at 10 m height AGL [m/s]	2D
FF10MAX	wind gusts at 10 m height AGL, with TKE component taken at first mass level \(4\sqrt{TKE_{IKB}}\) (only if CTURB=’TKEL’) [m/s]	2D
FF10MAX2	wind gusts at 10 m height AGL, with TKE component taken at 10m \(4\sqrt{TKE_{10m}}\) (only if CTURB=’TKEL’) [m/s]	2D
FF10MAX_AROME	wind gusts at 10 m height AGL, with TKE component taken at 20m \(3.8\sqrt{TKE_{20m}}\) as in AROME-France (only if CTURB=’TKEL’)	2D
SFCO2	\(CO_2\) flux if present in YINIFILE [mg/m²/s]	2D

NAM_DIAG - General

NAM_DIAG - default option

Fortran name	Fortran type	Default value
XDTSTEP	REAL	XTSTEP
CISO	CHARACTER(LEN=NFILENAMELGTMAX)	PREVTK
LVAR_RS	LOGICAL	TRUE

• XDTSTEP : time step of the DIAG program (one time step is performed). By default time step of the simulation is used (XTSTEP).

• add CISO="PREVTK" in NAM_DIAG to store following variables:

Name	Meaning [Unit]	Dimension
PABST	absolute pression [Pa]	3D
THT	potential temperature [K]	3D
POVOT	potential vorticity [PVU]	3D

Note

Other options for CISO are :

• “PR” to store PABST

• “TK” to store THT

• “EV” to store POVOT

• “PRTK” to store PABST + THT

• “PREV” to store PABST + POVOT

• “TKEV” to store THT + POVOT

• add LVAR_RS=T in NAM_DIAG to store following variables:

Name	Meaning [Unit]	Dimension
UT	u-wind speed [m/s]	3D
VT	v-wind speed [m/s]	3D
WT	w-wind speed [m/s]	3D
RVT	water vapor mixing ratio [kg/kg]	3D

Tip

Add LWIND_ZM=T (with LVAR_RS=T) in NAM_DIAG to store following variables :

Name	Meaning [Unit]	Dimension
UT_ZM	zonal wind speed [m/s]	3D
VT_ZM	meridian wind speed [m/s]	3D

• add LVAR_LS=T in NAM_DIAG to store following variables :

Name	Meaning [Unit]	Dimension
LSUM	large scale u-wind speed [m/s]	3D
LSVM	large scale v-wind speed [m/s]	3D
LSWM	large scale w-wind speed [m/s]	3D
LSTHM	large scale potential temperature [K]	3D
LSRVM	large scale water vapor mixing ratio [kg/kg]	3D

• add LVAR_LS=T and LWIND_ZM=T in NAM_DIAG to store following variables :

Name	Meaning [Unit]	Dimension
LSUM_ZM	large scale zonal wind speed [m/s]	3D
LSVM_ZM	large scale meridian wind speed [m/s]	3D

• add LVAR_FRC=T in NAM_DIAG (LFORCING has to be T in YINIFILE.des) to store following variables :

Name	Meaning [Unit]	Dimension
UFRCn	zonal component of horizontal forcing wind [m/s]	1D
VFRCn	meridian component of horizontal forcing wind [m/s]	1D
WFRCn	vertical forcing wind [m/s]	1D
THFRCn	\(\theta_{\rm frc}\) forcing potential temperature [K]	1D
RVFRCn	\(r_{v,\rm frc}\) forcing vapor mixing ratio [kg/kg]	1D
TENDTHFRCn	\((\partial \theta/\partial t)_{\rm frc}\) (K/s)	1D
TENDRVFRCn	\((\partial r_v/\partial t)_{\rm frc}\) ((kg/kg)/s)	1D
GXTHFRCn	\((\partial \theta/\partial x)_{\rm frc}\) (K/m)	1D
GYRVFRCn	\((\partial \theta/\partial y)_{\rm frc}\) (K/m)	1D
PGROUNDFRCn	forcing ground pressure (Pa)	0D

• add LTPZH=T in NAM_DIAG to store following variables :

Name	Meaning [Unit]	Dimension
TEMP	temperature [C]	3D
PRES	pressure [hPa]	3D
ALT	height of model levels (geopotentiel in pressure level) [m]	3D
REHU	relative Humidity [%] (if LUSERV=T)	3D
VPRES	vapor Pressure [hPa] (if LUSERV=T)	3D

• add LCOREF=T in NAM_DIAG to store following variables :

Name	Meaning [Unit]	Dimension
COREF	refraction coindex (if LUSERV=T)	3D
MCOREF	modified refraction coindex (if LUSERV=T)	3D

• add LMOIST_V=T in NAM_DIAG to store following variables :

Name	Meaning [Unit]	Dimension
THETAV	virtual potential Temperature [K]	3D
POVOV	virtual Potential Vorticity [PVU]	3D

• add LMOIST_V=T and LMEAN_POVO=T in NAM_DIAG to store following variables :

Name	Meaning [Unit]	Dimension
MEAN_POVOV	mean virtual potential vorticity (PVU)	2D

• add LMOIST_E=T in NAM_DIAG to store following variables :

Name	Meaning [Unit]	Dimension
THETAE	equivalent potential Temperature [K]	3D
POVOE	equivalent Potential Vorticity [PVU]	3D

• add LMOIST_E=T and LMEAN_POVO=T in NAM_DIAG to store following variables :

Name	Meaning [Unit]	Dimension
MEAN_POVOE	mean equivalent potential vorticity (PVU)	2D

• add LMOIST_ES=T in NAM_DIAG to store following variables :

Name	Meaning [Unit]	Dimension
THETAES	equivalent saturated potential temperature [K]	3D
POVOES	equivalent saturated potential vorticity [PVU]	3D

• add LMOIST_S1=T in NAM_DIAG to store following variables :

Name	Meaning [Unit]	Dimension
THETAS1	Moist air Entropy (1st order) potential temperature [K]	3D

• add LMOIST_S2=T in NAM_DIAG to store following variables :

Name	Meaning [Unit]	Dimension
THETAS2	Moist air Entropy (2nd order) potential temperature [K]	3D

• add LMOIST_L=T in NAM_DIAG to store following variables :

Name	Meaning [Unit]	Dimension
THETAL	Liquid water potential temperature [K]	3D

• add LMEAN_POVO=T in NAM_DIAG to store following variables :

Name	Meaning [Unit]	Dimension
MEAN_POVO	Mean Potential Vorticity (PVU)	2D

• add LMEAN_POVO=T and LMOIST_V=T in NAM_DIAG to store following variables :

Name	Meaning [Unit]	Dimension
MEAN_POVOV	Mean virtual Potential Vorticity (PVU)	2D

• add LMEAN_POVO=T and LMOIST_E=T NAM_DIAG to store following variables :

Name	Meaning [Unit]	Dimension
MEAN_POVOE	Mean equivalent Potential Vorticity (PVU)	2D

Note

Add XMEAN_POVO(1:2) in NAM_DIAG to chose averaged between two isobaric levels in Pa (XMEAN_POVO(1),XMEAN_POVO(2)) (by default (15000,50000))

• add LVORT=T in NAM_DIAG to store following variables :

Name	Meaning [Unit]	Dimension
ABVOR	Mean equivalent Potential Vorticity (PVU)	2D
UM1	u-relative vorticity components (/s)	3D
VM1	v-relative vorticity components (/s)	3D
WM1	w-relative vorticity components (/s)	3D

• add LVORT=T and LWIND_ZM=T in NAM_DIAG to store following variables :

Name	Meaning [Unit]	Dimension
UM1_ZM	zonal relative vorticity components (m/s)	3D
VM1_ZM	meridian relative vorticity components (m/s)	3D

• add LDIV=T in NAM_DIAG to store following variables :

Name	Meaning [Unit]	Dimension
HDIV	Horizontal divergence (/s)	3D
HMDIV	Horizontal Moisture divergence (kg/m3/s)	3D

• add LGEO=T in NAM_DIAG to store following variables :

Name	Meaning [Unit]	Dimension
UM88	Geostrophic u-wind components (m/s)	3D
VM88	Geostrophic v-wind components (m/s)	3D
WM88	Geostrophic w-wind components (m/s)	3D

• add LGEO=T and LWIND_ZM=T in NAM_DIAG to store following variables :

Name	Meaning [Unit]	Dimension
UM88_ZM	Geostrophic zonal wind components (m/s)	3D
VM88_ZM	Geostrophic meridian wind components (m/s)	3D

• add LAGEO=T in NAM_DIAG to store following variables :

Name	Meaning [Unit]	Dimension
UM89	Ageostrophic u-wind components (m/s)	3D
VM89	Ageostrophic v-wind components (m/s)	3D
WM89	Ageostrophic w-wind components (m/s)	3D

• add LAGEO=T and LWIND_ZM=T in NAM_DIAG to store following variables :

Name	Meaning [Unit]	Dimension
UM89_ZM	Ageostrophic zonal wind components (m/s)	3D
VM89_ZM	Ageostrophic meridian wind components (m/s)	3D

• add LMSLP=T in NAM_DIAG to store following variables :

Name	Meaning [Unit]	Dimension
MSLP	Mean Sea Level Pressure (hPa)	2D

• add LBV_FR=T in NAM_DIAG to store following variables :

Name	Meaning [Unit]	Dimension
BV	Brunt-Vaissala frequency (/s)	3D
BVE	Equivalent Brunt-Vaissala frequency (/s)	3D

• add LVAR_MRSV=T in NAM_DIAG to store following variables :

Name	Meaning [Unit]	Dimension
MRSVnnn	Mixing Ratio for User Scalar Variable n (g/kg)	3D

• define CBLTOP="RICHA" in NAM_DIAG to store following variable :

Name	Meaning [Unit]	Dimension
HBLTOP	Height of boundary layer top (m) computed with bulk Richardson number method	2D

• define CBLTOP="THETA" in NAM_DIAG to store following variable :

Name	Meaning [Unit]	Dimension
HBLTOP	Height of boundary layer top (m) computed with parcel method	2D

• add LHU_FLX=T in NAM_DIAG to store following variables :

Name	Meaning [Unit]	Dimension
UM90	u-wind components of moisture ground flux (kg/s/m2)	3D
VM90	v-wind components of moisture ground flux (kg/s/m2)	3D
UM91	u-wind components of moisture ground flux integrated on 3000 meters (kg/s/m)	2D
VM91	v-wind components of moisture ground flux integrated on 3000 meters (kg/s/m)	2D
HMCONV	Horizontal CONVergence of moisture flux (kg/s/m2)	2D
HMCONV3000	Horizontal CONVergence of moisture flux integrated on 3000 meters (kg/s/m2)	2D
UM92	u-wind components of hydrometeores ground flux (if CCLOUD=ICE3 or ICE4) (kg/s/m2)	2D
VM92	v-wind components of hydrometeores ground flux (if CCLOUD=ICE3 or ICE4) (kg/s/m2)	2D
UM93	u-wind components of hydrometeor ground flux (if CCLOUD=ICE3 or ICE4) (kg/s/m)	2D
VM93	v-wind components of hydrometeor ground flux (if CCLOUD=ICE3 or ICE4) (kg/s/m)	2D
HMCONV_TT	Horizontal CONVergence of hydrometeor flux (kg/s/m2)	2D
HMCONV3000_TT	Horizontal CONVergence of hydrometeor flux integrated on 3000 meters (kg/s/m2)	2D

• define NCAPE=0 in NAM_DIAG to store following variables :

Name	Meaning [Unit]	Dimension
CAPEMAX	maximum of CAPE3D (J/kg)	2D
CINMAX	value of CIN3D corresponding to CAPEMAX (J/kg)	2D

• define NCAPE=1 in NAM_DIAG to store following variables :

Name	Meaning [Unit]	Dimension
CAPEMAX	maximum of CAPE3D (J/kg)	2D
CINMAX	value of CIN3D corresponding to CAPEMAX (J/kg)	2D
CAPE3D	Convective Available Potential Energy (J/kg)	3D
CIN3D	Convective INhibition energy (J/kg)	3D
DCAPE3D	Downdraft cape (J/kg)	3D

• define NCAPE=2 in NAM_DIAG to store following variables :

Name	Meaning [Unit]	Dimension
CAPEMAX	maximum of CAPE3D (J/kg)	2D
CINMAX	value of CIN3D corresponding to CAPEMAX (J/kg)	2D
CAPE3D	Convective Available Potential Energy (J/kg)	3D
CIN3D	Convective INhibition energy (J/kg)	3D
DCAPE3D	Downdraft cape (J/kg)	3D
VKE	Vertical Kinetic Energy (from explicit vertical motion) (J/kg)	3D

NAM_DIAG - Deep convection

• define NCONV_KF=0 in NAM_DIAG to store following variables :

Name	Meaning [Unit]	Dimension
CAPE	Convective Available Potentiel Energy (J/kg)	2D
CLTOPCONV	top of convective clouds(km)	2D
CLBASCONV	base of convective clouds(km)	2D

• define NCONV_KF=1 in NAM_DIAG to store following variables :

Name	Meaning [Unit]	Dimension
CAPE	Convective Available Potentiel Energy (J/kg)	2D
CLTOPCONV	top of convective clouds(km)	2D
CLBASCONV	base of convective clouds(km)	2D
DTHCONV	Convective tendency for potential temperature (K/s)	3D
DRVCONV	Convective tendency for vapor (/s)	3D
DRCCONV	Convective tendency for cloud (/s)	3D
DRICONV	Convective tendency for ice (/s)	3D
DSVCONVnn	Convective tendency for scalar variables (/s)	3D

• define NCONV_KF=2 in NAM_DIAG to store following variables :

Name	Meaning [Unit]	Dimension
CAPE	Convective Available Potentiel Energy (J/kg)	2D
CLTOPCONV	top of convective clouds(km)	2D
CLBASCONV	base of convective clouds(km)	2D
DTHCONV	Convective tendency for potential temperature (K/s)	3D
DRVCONV	Convective tendency for vapor (/s)	3D
DRCCONV	Convective tendency for cloud (/s)	3D
DRICONV	Convective tendency for ice (/s)	3D
DSVCONVnnn*	Convective tendency for scalar variables (/s)	3D
UMFCONV	Updraft Convective Mass Flux (m2 kg/s)	3D
DMFCONV	Downdraft Convective Mass Flux (m2 kg/s)	3D
PRLFLXCONV	Liquid PRecipitation Convective FLuX (m/s)	3D
PRSFLXCONV	Solid PRecipitation Convective FLuX (m/s)	3D

NAM_DIAG - Shallow convection

• add LMFFLX=T in NAM_DIAG to store following variables :

Name	Meaning [Unit]	Dimension
MF_THW_FLX	Conservative potential temperature vertical flux (K*m/s)	3D
MF_RCONSW_FLX	Conservative mixing ratio vertical flux (kg/kg*m/s)	3D
MF_THVW_FLX	Theta_v vertical flux (K*m/s)	3D
MF_UW_VFLX	U momentum vertical flux (m2/s2)	3D
MF_VW_VFLX	V momentum vertical flux (m2/s2)	3D

NAM_DIAG - Microphysics

• add LTHW=T in NAM_DIAG to store following variables :

Name	Meaning [Unit]	Dimension
THVW	Thickness of Vapor Water (mm) (if LUSERV=T)	2D
THCW	Thickness of Cloud Water (mm) (if LUSERC=T)	2D
THRW	Thickness of Rain Water (mm) (if LUSERR=T)	2D
THIC	Thickness of Ice (mm) (if LUSERI=T)	2D
THSN	Thickness of Snow (mm) (if LUSERS=T)	2D
THGR	Thickness of Graupel (mm) (if LUSERG=T)	2D
THHA	Thickness of Hail (mm) (if LUSERH=T)	2D

• add LVAR_MRW=T in NAM_DIAG to store following variables :

Name	Meaning [Unit]	Dimension
MRV	Mixing Ratio for Vapor (g/kg)(if LUSERV=T)	3D
MRC	Mixing Ratio for Cloud (g/kg) (if LUSERC=T)	3D
MRR	Mixing Ratio for Rain (g/kg) (if LUSERR=T)	3D
MRI	Mixing Ratio for Ice (g/kg) (if LUSERI=T)	3D
CIT	Ice concentration (m-3 (if LUSECI=T)	3D
MRS	Mixing Ratio for Snow (g/kg) (if LUSERS=T)	3D
MRG	Mixing Ratio for Graupel (g/kg) (if LUSERG=T)	3D
MRH	Mixing Ratio for Hail (g/kg) (if LUSERH=T)	3D
CCCN	if CCLOUD=’C2R2’	3D
CCLOUD	if CCLOUD=’C2R2’	3D
CRAIN	if CCLOUD=’C2R2’	3D
SUPSAT	if CCLOUD=’C2R2’ and LSUPSAT=T	3D
CICE	if CCLOUD=’C1R3’	3D
CIN	if CCLOUD=’C1R3’	3D

• add LVAR_PR=T in NAM_DIAG to store following variables :

Name	Meaning [Unit]	Dimension
ACPRR	Accumulated explicit Precipitation Rate for Rain (mm) (accumulated from the beginning of the simulation) (if LUSERR=T)	2D
INPRR	Instantaneous explicit Precipitation Rate (mm/h) (if LUSERR=T)	2D
INPRR3D	Instantaneous explicit 3D Rain Precipitation flux (m/s) (if LUSERR=T)	3D
EVAP3D	Instantaneous 3D Rain Evaporation flux (kg/kg/s) (if LUSERR=T)	3D
ACPRC	Accumulated Cloud Precipitation Rate (mm) (if LUSERC=T)	2D
INPRC	Instantaneous Cloud Precipitation Rate (mm/h) (if LUSERC=T)	2D
ACPRS	Accumulated explicit Precipitation Rate for Snow (mm) (if LUSERS=T)	2D
INPRS	Instantaneous explicit Precipitation Rate for Snow (mm/h) (if LUSERS=T)	2D
ACPRG	Accumulated explicit Precipitation Rate for Graupel (mm) (if LUSERG=T)	2D
INPRG	Instantaneous explicit Precipitation Rate for Graupel (mm/h) (if LUSERG=T)	2D
ACPRH	Accumulated explicit Precipitation Rate for Hail (mm) (if LUSERH=T)	2D
INPRH	Instantaneous explicit Precipitation Rate for Hail (mm/h) (if LUSERH=T)	2D
ACPRT	[2D] Total Accumulated explicit Precipitation Rate (mm) (if LUSERR=T)	2D
INPRT	Total Instantaneous explicit Precipitation Rate (mm/h) (if CCLOUD not NONE)	2D
PACCONV	Convective Accumulated Precipitation Rate (mm) (if CDCONV not NONE)	2D
PRCONV	Convective Instantaneous Precipitation Rate (mm/h) (if CDCONV not NONE)	2D
PRSCONV	Convective instantaneous Precipitation Rate for Snow (mm/h) (if CDCONV not NONE)	2D
PRECIP_WAT	Precipitable water (kg/m2)	2D

• add LCHAQDIAG=T in NAM_DIAG to store following variables :

Name	Meaning [Unit]	Dimension
WC_O3	Chemical scalar variables in aqueous phase (cloud and rain) as defined in BASIC.f90 (M)	3D

• add LTOTAL_PR=T in NAM_DIAG to store following variables :

Name	Meaning [Unit]	Dimension
ACTOPR	Accumulated Total Precipitation (mm)	2D
INTOPR	Instantaneous Total Precipitation (mm/h)	2D

• add LTOTAL_PR=T and LMEAN_PR=T in NAM_DIAG to store following variables :

Name	Meaning [Unit]	Dimension
LS_ACTOPR	Large scale accumulated Total Precipitation (mm)	2D
LS_INTOPR	Large scale instantaneous Total Precipitation (mm/h)	2D

• XMEAN_PR (1,1) nb of grid points of the small-scale model inside the LS grid mesh along x, y for LMEAN_PR

• add LCLD_COV=T in NAM_DIAG to store following variables :

Name	Meaning [Unit]	Dimension
HECL	Height of Explicit CLoud top (km)	2D
HCL	Height of maximum CLoud top (km)	2D
TCL	Temperature of maximum Cloud top)	2D
CLDFR	Cloud Fraction (_)	3D
VISI_HOR	Visibility (m)	3D

• add LLIMA_DIAG=T in NAM_DIAG to store following variables :

Name	Meaning [Unit]	Dimension
MRV	Mixing Ratio for Vapor (g/kg)(if LUSERV=T)	3D
MRC	Mixing Ratio for Cloud (g/kg)	3D
MRR	Mixing Ratio for Rain (g/kg)	3D
MRI	Mixing Ratio for Ice (g/kg)	3D
MRS	Mixing Ratio for Snow (g/kg)	3D
MRG	Mixing Ratio for Graupel (g/kg)	3D
MRH	Mixing Ratio for Hail (g/kg) (if LUSERH=T)	3D
NCT	Cloud concentration (\(\rm cm^{-3}\))	3D
NRT	Rain concentration (\(\rm cm^{-3}\))	3D
NFREE	Free CCN concentration (\(\rm cm^{-3}\))	3D
NCCN	CCN concentration (\(\rm cm^{-3}\))	3D
MASSAP	Scavenging (\(\rm kg.cm^{-3}\))	3D
CICE	Ice concentration (\(\rm cm^{-3}\))	3D
CIFNFREE	Free IFN concentration (\(\rm cm^{-3}\))	3D
CIFNNUCL	Nucleated IFN concentration (\(\rm cm^{-3}\))	3D
CCNINIMM	Nucleated IMM concentration (\(\rm cm^{-3}\))	3D
CCCNNUCL	Homogeneous Freezing of CCN (\(\rm cm^{-3}\))	3D
LWC	Liquid Water content (\(\rm g.m^{-3}\)) (if LUSERC=T)	3D
IWC	Ice Water content (\(\rm g.m^{-3}\)) (if LUSERC=T)	3D

• add LVISI=T in NAM_DIAG to store following variables :

Name	Meaning [Unit]	Dimension
VISIKUN	Visibility from Kunkel (m) (if CCLOUD/=REVE or NONE)	3D
VISIGUL	Visibility from Gultepe (m) (if CCLOUD=C2R2 or KHKO)	3D
VISIZHA	Visibility from Zhang (m) (if CCLOUD=C2R2 or KHKO)	3D

NAM_DIAG - Turbulence

• define LVAR_TURB=T in NAM_DIAG to store following variables :

Name	Meaning [Unit]	Dimension
TKET	Turbulent Kinetic Energy (m2/s2)	3D
SIGS	Sigma_s from turbulence scheme (kg/kg2)	3D
SRCM	Normalized 2nd_order moment (kg/kg2)	3D
BL_DEPTH	Boundary Layer Depth if CTOM=’TM06’ (m)	3D

• define LTURBDIAG=T in NAM_DIAG to store following variables :

Name	Meaning [Unit]	Dimension
AMOIST	See Scientific documentation part III, chap 7 equation 7.29	3D
ATHETA	See Scientific documentation part III, chap 7 equation 7.30	3D
RED_TH1, RED_R1, RED2_TH3, RED2_R3, RED2_THR3	Redelsperger numbers	3D
TKE_DP	dynamical production of TKE (m2/s3)	3D
TKE_TP	thermal production of TKE (m2/s3)	3D
TKE_TR	transport of TKE (m2/s3)	3D
TKE_DISS	dissipation of TKE (m2/s3)	3D
LM_CLEAR_SKY	mixing length in clear sky (m)	3D
COEF_AMPL	amplification of the mixing length (-)	3D
LM_CLOUD	mixing length in the clouds (m)	3D
LM	mixing length (m)	3D
THLM	conservative potential temperature (K)	3D
RNPM	conservative mixing ratio (kg/kg)	3D
RVCI	rv + rc + ri (kg/kg)	3D
GX_RVCI,GY_RVCI	x and y gradient of RVCI (kg/kg/m)	3D
GNORM_RVCI	Horizontal norm of the gradient of RVCI (kg/kg/m)	3D
QX_RVCI	x gradient of the advection of RVCI (kg/kg/m)	3D
QY_RVCI	y gradient of the advection of RVCI (kg/kg/m)	3D
QNORM_RVCI	Horizontal norm of the gradient of the advection of RVCI (kg/kg/m)	3D
CEI	Cloud entrainment instability index (kg/kg/m/s)	3D

• define LTURBFLX=T in NAM_DIAG to store following variables :

Name	Meaning [Unit]	Dimension
PHI3	Turbulent Prandtl number (-)	3D
PSI3	Turbulent Schmidt number (-)	3D
PSI_SV_n	Turbulent Schmidt number for the scalar variables (-)	3D

• define LTURBFLX=T in NAM_DIAG and if CTURBDIM=’1D’ in YINIFILE.des following variables will stored :

Name	Meaning [Unit]	Dimension
THW_FLX	theta vertical flux (K*m/s)	3D
RCONSW_FLX	rv vertical flux (kg*m/s/kg)	3D
RCW_FLX	liquid water mixing ratio vertical flux (kg*m/s/kg)	3D
THL_VVAR	< T Hl, T Hl > (K2)	3D
THLRCONS_VCOR	< T Hl, Rnp >( K*kg/kg)	3D
RTOT_VVAR	< Rnp, Rnp > ((kg/kg)2)	3D
UW_VFLX, VW_VFLX	wind component vertical flux (m2/s2)	3D
WSV_FLX_n	< W, SV th >(SVunit m/s)	3D

• define LTURBFLX=T in NAM_DIAG and if CTURBDIM=’3D’ in YINIFILE.des following variables will stored :

Name	Meaning [Unit]	Dimension
THW_FLX	theta vertical flux (K*m/s)	3D
RCONSW_FLX	rv vertical flux (kg*m/s/kg)	3D
RCW_FLX	liquid water mixing ratio vertical flux (kg*m/s/kg)	3D
THL_VVAR	< T Hl, T Hl > (K2)	3D
THLRCONS_VCOR	< T Hl, Rnp >( K*kg/kg)	3D
RTOT_VVAR	< Rnp, Rnp > ((kg/kg)2)	3D
UW_VFLX, VW_VFLX	wind component vertical flux (m2/s2)	3D
WSV_FLX_n	< W, SV th >(SVunit m/s)	3D
U_VAR	U variance ((m/s)2)	3D
V_VAR	V variance ((m/s)2)	3D
W_VAR	W variance ((m/s)2)	3D
UV_FLX	< U, V >((m/s)2,)	3D
UW_HFLX	< U, W > ((m/s)2)	3D
VW_HFLX	< V, W > ((m/s)2)	3D
USV_FLX_n	< U, SV th > ( SVunit m/s)	3D
VSV_FLX_n	< V, SV th > ( SVunit m/s)	3D
THL_HVAR	< T Hl, T Hl > (K2)	3D
THLR_HCOR	< T Hl, Rnp > (K*kg/kg)	3D
R_HVAR	< Rnp, Rnp > ( (kg/kg)2)	3D
UTHL_FLX	horizontal < U, T Hl > (K*m/s)	3D
VTHL_FLX	horizontal < V, T Hl > (K*m/s)	3D
UR_FLX	horizontal < U, Rnp > (kg/kg*m/s)	3D
VR_FLX	horizontal < V, Rnp > (kg/kg*m/s)	3D

NAM_DIAG - Radiation

Only available if CRAD not NONE in YINIFILE.des.

• add NRAD_3D=0 in NAM_DIAG to store following variables :

Name	Meaning [Unit]	Dimension
DTHRAD	Radiative heating/cooling rate (K/s)	3D
FLALWD	Downward LW on FLAT surface (W/m2)	2D
DIRFLASWD	Direct Downward SW on FLAT surface (W/m2)	2D
SCAFLASWD	Scattered Downward SW on FLAT surface (W/m2)	2D
DIRSRFSWD	Direct Downward SW (W/m2)	2D
CLEARCOL_TM1	trace of cloud (-)	2D
EMIS	emmissivity (-)	2D
ZENITH	solar zenithal angle (RAD)	2D
AZIM	azimuthal angle (RAD)	2D
DIR_ALB	direct albedo(-)	2D
SCA_ALB	scattered albedo (-)	2D
TSRAD	radiative surface temperature (K)	2D

• add NRAD_3D=1 in NAM_DIAG to store following variables (in addition to NRAD_3D=0 variables):

Name	Meaning [Unit]	Dimension
SWU	Upward SW radiative fluxes (W/m2)	2D
SWD	Downward SW radiative fluxes (W/m2)	2D
LWU	Upward LW radiative fluxes (W/m2)	2D
LWD	Downward LW radiative fluxes (W/m2)	2D
SWF_DOWN	Downward SW radiative fluxes (W/m2)	3D
SWF_UP	Upward SW radiative fluxes (W/m2)	3D
LWF_DOWN	Downward LW radiative fluxes (W/m2)	3D
LWF_UP	Upward LW radiative fluxes (W/m2)	3D
LWF_NET	Total LW net radiative fluxes (W/m2)	3D
SWF_NET	Total SW radiative fluxes (W/m2)	3D
DTRAD_LW	LW radiative tendency for T (K/day)	3D
DTRAD_SW	SW radiative tendency for T (K/day)	3D
RADSWD_VIS	surface radiative flux in visible (W/m2)	2D
RADSWD_NIR	surface radiative flux in near-infrared (W/m2)	2D
RADLWD	LW surface radiative flux (W/m2)	2D

• add NRAD_3D=2 in NAM_DIAG to store following variables (in addition to NRAD_3D=1,0 variables):

Name	Meaning [Unit]	Dimension
SWF_DOWN_CS	Clear-Sky Downward SW radiative fluxes (W/m2)	3D
SWF_UP_CS	Clear-Sky Upward SW radiative fluxes (W/m2)	3D
LWF_DOWN_CS	Clear-Sky Downward LW radiative fluxes (W/m2)	3D
LWF_UP_CS	Clear-Sky Upward LW radiative fluxes (W/m2)	3D
LWF_NET_CS	Clear-Sky Total LW net radiative fluxes (W/m2)	3D
SWF_NET_CS	Clear-Sky Total SW radiative fluxes (W/m2)	3D
DTRAD_LW_CS	Clear-Sky LW radiative tendency for T (K/day)	3D
DTRAD_SW_CS	Clear-Sky SW radiative tendency for T (K/day)	3D
RADSWD_NIR_CS	Clear-Sky surface radiative flux in near-infrared (W/m2)	2D
RADLWD_CS	Clear-Sky LW surface radiative flux (W/m2)	2D

• add NRAD_3D=3 in NAM_DIAG to store following variables (in addition to NRAD_3D=2,1,0 variables):

Name	Meaning [Unit]	Dimension
PLAN_ALB_VIS	Planetary albedo in visible (-)	[2D]
PLAN_ALB_NIR	Planetary albedo in near-infrared (-)	[2D]
PLAN_TRA_VIS	Planetary transmission in visible (-)	[2D]
PLAN_TRA_NIR	Planetary transmission in near-infrared (-)	[2D]
PLAN_ABS_VIS	Planetary absorption in visible (-)	[2D]
PLAN_ABS_NIR	Planetary absorption in near-infrared (-)	[2D]

• add NRAD_3D=4 in NAM_DIAG to store following variables (in addition to NRAD_3D=3,2,1,0 variables):

Name	Meaning [Unit]	Dimension
EFNEB_UP	Upward equivalent emissivity (Morcrette scheme)(-)	3D
EFNEB_DOWN	Downward equivalent emissivity (-)	3D
FLWP	Liquid water path (g/m2)	3D
FIWP	Ice water path (g/m2)	3D
EFRADL	Cloud liquid water effective radius (\(mu\) m)	3D
EFRADI	Cloud ice effective radius (\(mu\) m)	3D
SW_NEB	Effective cloud fraction (-)	3D
RRTM_LW_NEB	Effective cloud fraction (-)	3D
OTH_VIS	Cloud optical thickness (-)	3D
OTH_NI1	Cloud optical thickness (-)	3D
OTH_NI2	Cloud optical thickness (-)	3D
OTH_NI3	Cloud optical thickness (-)	3D
SSA_VIS	Cloud single scattering albedo (-)	3D
SSA_NI1	Cloud single scattering albedo (-)	3D
SSA_NI2	Cloud single scattering albedo (-)	3D
SSA_NI3	Cloud single scattering albedo (-)	3D
ASF_VIS	Cloud asymetry factor (-)	3D
ASF_NIR1	Cloud asymetry factor (-)	3D
ASF_NIR2	Cloud asymetry factor (-)	3D
ASF_NIR3	Cloud asymetry factor (-)	3D
ODAER_VIS		3D
ODAER_NIR1		3D
ODAER_NIR2		3D
ODAER_NIR3		3D
SSAAER_VIS		3D
SSAAER_NIR1		3D
SSAAER_NIR2		3D
SSAAER_NIR3		3D
GAER_VIS		3D
GAER_NIR1		3D
GAER_NIR2		3D
GAER_NIR3		3D

• add NRAD_3D=5 in NAM_DIAG to store following variables (in addition to NRAD_3D=4,3,2,1,0 variables):

Name	Meaning [Unit]	Dimension
O3CLIM	Climatological ozone content (Pa/Pa)	3D
CUM_AER_LAND	Cumulated optical thickness of the different aerosols from the top of the domain	3D
CUM_AER_SEA	Cumulated optical thickness of the different aerosols from the top of the domain	3D
CUM_AER_DES	Cumulated optical thickness of the different aerosols from the top of the domain	3D
CUM_AER_URB	cumulated optical thickness of the different aerosols from the top of the domain	3D
CUM_AER_VOL	cumulated optical thickness of the different aerosols from the top of the domain	3D
CUM_AER_STRB	cumulated optical thickness of the different aerosols from the top of the domain	3D

NAM_DIAG - Aerosols

• ORILAM (Only available if LUSECHEM=T and LORILAM=T in YINIFILE.des.), add LCHEMDIAG=T in NAM_DIAG to store following variables :

Name	Meaning [Unit]	Dimension
SOAI…	Aerosol scalar variable as defined in ch_aer_init_soa.f90 (ppb)	3D
RGAn	Aerosol number mean Radius of the lognormal mode n (\(mu\) m)	3D
RGAMn	Aerosol Mass mean Radius of the lognormal mode n (\(mu\) m)	3D
N0An	Aerosol Number of the lognormal mode n (/cc)	3D
SIGAn	Aerosol Standard deviation of the lognormal mode n (-)	3D
MSO4n	Mass SO4 aerosol mode n (\(mu\) m/m3)	3D
MNO3n	Mass NO3 aerosol mode n (\(mu\) m/m3)	3D
MNH3n	Mass NH3 aerosol mode n (\(mu\) m/m3)	3D
MH2On	Mass H2O aerosol mode n (\(mu\) m/m3)	3D
MSOA1n	Mass SOA1 aerosol mode n (\(mu\) m/m3)	3D
MSOA2n	Mass SOA2 aerosol mode n (\(mu\) m/m3)	3D
MSOA3n	Mass SOA3 aerosol mode n (\(mu\) m/m3)	3D
MSOA4n	Mass SOA4 aerosol mode n (\(mu\) m/m3)	3D
MSOA5n	Mass SOA5 aerosol mode n (\(mu\) m/m3)	3D
MSOA6n	Mass SOA6 aerosol mode n (\(mu\) m/m3)	3D
MSOA7n	Mass SOA7 aerosol mode n (\(mu\) m/m3)	3D
MSOA8n	Mass SOA8 aerosol mode n (\(mu\) m/m3)	3D
MSOA9n	Mass SOA9 aerosol mode n (\(mu\) m/m3)	3D
MSOA10n:	Mass SOA10 aerosol mode n (\(mu\) m/m3)	3D
MOCn	Mass OC aerosol mode n (\(mu\) m/m3)	3D
MBCn	Mass BC aerosol mode n (\(mu\) m/m3	3D

• Dust (Only available if LDUST=T in YINIFILE.des.)

Name	Meaning [Unit]	Dimension
DSTM0n	Dust 0-order moment of the lognormal mode n (ppb)	3D
DSTM3n	Dust 3rd-order moment of the lognormal mode n (ppb)	3D
DSTM6n	Dust 6rd-order moment of mode n (if LVARSIG) (ppb)	3D
DSTRGAn	Dust number mean Radius of the lognormal mode n (\(mu\) m)	3D
DSTRGAMn	Dust Mass mean Radius of the lognormal mode n (\(mu\) m)	3D
DSTN0An	Dust Number of the lognormal mode n (/m3)	3D
DSTSIGAn	Dust Standard deviation of the lognormal mode n (-)	3D
DSTMSSn	Dust Mass concentration of the lognormal mode n (\(mu\) g/m3)	3D
DSTBRDNn	Dust Burden of the lognormal mode n (g/m2)	2D
DEDSTM3nC	Dust Mass of mode n in cloud water only if LDEPOS_DST (ppb)	3D
DEDSTM3nR	Dust Mass of mode n in rain only if LDEPOS_DST=T (ppb)	3D
DEDSTN0An	Number of dust particles in cloud water (for n=1,2,3) or in rain (for n=4,5,6) only if LDEPOS_DST=T (/m3)	3D
DEDSTMSSn	Dust mass in cloud water (for n=1,2,3) or in rain (for n=4,5,6) only if LDEPOS_DST=T(\(mu\) g/m3)	3D
DSTAOD2D	Dust Optical Depth (-) if NRAD_3D \(geq\) 1 in NAM_DIAG - Radiation	2D
DSTAOD3D	Dust Optical Depth between two vertical levels (-) if NRAD_3D \(geq\) 1 in NAM_DIAG - Radiation	3D
DSTEXT	Dust Extinction (1/km) if NRAD_3D \(geq\) 1 in NAM_DIAG - Radiation	3D

• Salt (Only available if LSALT=T in YINIFILE.des.)

Name	Meaning [Unit]	Dimension
SLTM0n	Salt 0-order moment of the lognormal mode n (ppb)	3D
SLTM3n	Salt 3rd-order moment of the lognormal mode n (ppb)	3D
SLTM6n	Salt 6rd-order moment of mode n (if LVARSIG_SLT) (ppb)	3D
SLTRGAn	Salt number mean Radius of the lognormal mode n (\(mu\) m)	3D
SLTRGAMn	Salt Mass mean Radius of the lognormal mode n (\(mu\) m)	3D
SLTN0An	Salt Number of the lognormal mode n (/m3)	3D
SLTSIGAn	Salt Standard deviation of the lognormal mode n (-)	3D
SLTMSSn	Salt Mass concentration of the lognormal mode n (\(mu\) g/m3)	3D
SLTBRDNn	Salt Burden of the lognormal mode n (g/m2)	2D
DESLTM3nC	Salt Mass of mode n in cloud water only if LDEPOS_SLT=T (ppb)	3D
DESLTM3nR	Salt Mass of mode n in rain only if LDEPOS_SLT=T (ppb)	3D
DESLTN0An	Number of salt particles in cloud water (for n=1,2,3) or in rain (for n=4,5,6) only if LDEPOS_SLT=T (/m3)	3D
DESLTMSSn	Salt mass in cloud water (for n=1,2,3) or in rain (for n=4,5,6) only if LDEPOS_SLT=T (\(mu\) g/m3)	3D
SLTAOD2D	Salt Optical Depth (-) if NRAD_3D \(geq\) 1 in NAM_DIAG - Radiation	2D
SLTAOD3D	Salt Optical Depth between two vertical levels (_) if NRAD_3D \(geq\) 1 in NAM_DIAG - Radiation	3D
SLTEXT	Salt Extinction (1/km) if NRAD_3D \(geq\) 1 in NAM_DIAG - Radiation	3D

NAM_DIAG - Chemistry

Only available if LUSECHEM=T in YINIFILE.des

• add LCHEMDIAG=T in NAM_DIAG to store following variables :

Name	Meaning [Unit]	Dimension
O3 …	Chemical scalar variables as defined in BASIC.f90 (ppb)	3D

Note

Following options can be used in NAM_DIAG when LCHEMDIAG is activated :

Fortran name	Fortran type	Default value
XCHEMLAT	REAL	XUNDEF
XCHEMLON	REAL	XUNDEF
CSPEC_DIAG	CHARACTER(LEN=1024)
CSPEC_BU_DIAG	CHARACTER(LEN=1024)

• write chemicals species on vertical profile defined by (XCHEMLAT,XCHEMLON)

• CSPEC_DIAG : list of the chemical species for production/loss terms computation. Each species is separated by a comma. Ex : CSPEC_DIAG=’O3,CO,BIO’

• CSPEC_BU_DIAG : list of the chemical species for production/loss terms computation in all the reactions where the species is involved. Each species is separated by a comma. Ex : CSPEC_BU_DIAG=’O3,CO,BIO’

NAM_DIAG - Lightnings

If LCH_CONV_LINOX=T and LUSECHEM=F in YINIFILE.des with LCHEMDIAG=F in DIAG1.nam, following variables will be stored :

Name	Meaning [Unit]	Dimension
LINOX	Linox scalar variables (ppb)	3D
IC_RATE	IntraCloud lightning Rate (/s)	2D
CG_RATE	CloudGround lightning Rate (/s)	2D
IC_TOTAL_NB	IntraCloud lightning Number (-)	2D
CG_TOTAL_NB	CloudGround lightning Number (-)	2D

NAM_DIAG - Lagrangian tracers

Only available if LLG=T in YINIFILE.des.

• add LTRAJ=T in NAM_DIAG to store following variables :

Name	Meaning [Unit]	Dimension
X	X coordinates (km)	3D
Y	Y coordinates (km)	3D
LGX	X Lagrangian tracers coordinates (m)	3D
LGY	Y Lagrangian tracers coordinates (m)	3D
LGZ	Z Lagrangian tracers coordinates (m)	3D
X_TRAJ	X Lagrangian tracers coordinates at time origin n	4D
Y_TRAJ	Y Lagrangian tracers coordinates at time origin n	4D
Z_TRAJ	Z Lagrangian tracers coordinates at time origin n	4D
THT_TRAJ	corresponding Theta (K)	4D
MRV_TRAJ	corresponding Vapor mixing Ratio (g/kg)	4D

NAM_DIAG - GPS simulator

• put NGPS=0 in NAM_DIAG to store following variable :

Name	Meaning [Unit]	Dimension
ZTD	Zenithal Total Delay (m)	2D

• put NGPS=1 in NAM_DIAG to store following variables :

Name	Meaning [Unit]	Dimension
ZTD	Zenithal Total Delay (m)	2D
ZHD	Zenithal Hydrostatic Delays (m)	2D
ZWD	Zenithal Wet Delays (m)	2D

Warning

Following options has to be set in NAM_DIAG when NGPS > -1 :

Fortran name	Fortran type	Default value
CNAM_GPS	ARRAY(CHARACTER)	50*
XLAT_GPS	ARRAY(REAL)	50* XUNDEF
XLON_GPS	ARRAY(REAL)	50* XUNDEF
XZS_GPS	ARRAY(REAL)	50* -999.0
XDIFFORO	REAL	150.0

• CNAM_GPS : name of the GPS stations

• XLAT_GPS : latitude of the GPS stations

• XLON_GPS : longitude of the GPS stations

• XZS_GPS : height of the GPS stations (m)

• XDIFFORO : maximum difference between model orography and station height accepted when computing interpolated delays value (m)

For stations where latitude, longitude and height are different from default values, the interpolated values of GPS delays are written in ASCII files YINIFILEYSUFFIXGPS[.P00n] (where n is the number of processor).

NAM_DIAG - Satellite simulator

Since RTTOV requires a license agreement, no RTTOV package is included in the open source version of Meso-NH. However, a subroutine calling RTTOV14 is included in Meso-NH version 60. To compile and use Meso-NH with RTTOV, you must follow the instructions in Radiative computation (RTTOV).

Fortran name	Fortran type	Default value
NRTTOVinfo	INTEGER(1:4,10)	999

For each instrument nb you want (1 \(\leq\) nb \(\leq\) 10) :

• NRTTOVinfo(1,nb) : Plt = Plateforme

• NRTTOVinfo(2,nb) : Sat = Satellite

• NRTTOVinfo(3,nb) : Sen = Sensor

• NRTTOVinfo(4,nb) : 0

Following variable will be store :

Name	Meaning [Unit]	Dimension
BT	Brightness temperature (K)	2D

To simulate an instrument, use the code given in the following tables reproduced from the RTTOV website. For example, to simulate both all the SEVIRI channels of MSG-2 and all the AMSU-B channels of NOAA-16, you can use the following namelist in DIAG1.nam :

&NAM_DIAG NRTTOVinfo(:,1)= 12 2 21 0,
          NRTTOVinfo(:,2)= 1 16 4 0 /

For this specific choice, you need to have four files: rtcoef_msg_2_seviri.dat, sccldcoef_msg_2_seviri.dat, rtcoef_noaa_16_amsub.dat and hydrotable_noaa_amsub.dat. The AMSU-B coefficient files can be found in directories rttov13pred54L and hydrotable (or to be downloaded from the RTTOV website).

Tip

To download all RTTOV files (~ 30 Go), you can use the script RTTOV_14.1/rtcoef_rttov14/rttov_coef_download.sh :

./rttov_coef_download.sh

The MSG-2 coefficient files are aliased following:

• for IR only calculation, do:

ln -sf $(SRC_MESONH)/src/LIB/RTTOV_14.1/rtcoef_rttov14/rttov13pred54L/rtcoef_msg_2_seviri_7gas_ironly.dat rtcoef_msg_2_seviri.dat
ln -sf $(SRC_MESONH)/src/LIB/RTTOV_14.1/cldaer_ir/sccldcoef_msg_2_seviri_ironly.dat sccldcoef_msg_2_seviri.dat

• for visible and IR calculation, do:

ln -sf $(SRC_MESONH)/src/LIB/RTTOV_14.1/rttov13pred54L/rtcoef_msg_2_seviri_o3co2.dat rtcoef_msg_2_seviri.dat
ln -sf $(SRC_MESONH)/src/LIB/RTTOV_14.1/cldaer_visir/sccldcoef_msg_2_seviri.dat .

In addition, download the brdf_data file and set NRTTOVinfo(4,1) to 1.

List of Plt and Sat variables

Platform	Plt	Sat range
NOAA	1	1 to 18
DMSP	2	8 to 16
Meteosat	3	5 to 7
GOES	4	8 to 12
GMS	5	5
FY-2	6	2 to 3
TRMM	7	1
ERS	8	1 to 2
EOS	9	1 to 2
METOP	10	1 to 3
ENVISAT	11	1
MSG	12	1 to 2
FY-1	13	3
ADEOS	14	1 to 2
MTSAT	15	1
CORIOLIS	16	1

List of Sen variables

Sensor	RTTOVid (Sen)	Sensor Channel	RTTOV-8 Channel
HIRS	0	1 to 19	1 to 19
MSU	1	1 to 4	1 to 4
SSU	2	1 to 3	1 to 3
AMSU-A	3	1 to 15	1 to 15
AMSU-B	4	1 to 5	1 to 5
AVHRR	5	3b to 5	1 to 3
SSMI	6	1 to 7	1 to 4
VTPR1	7	1 to 8	1 to 8
VTPR2	8	1 to 8	1 to 8
TMI	9	1 to 9	1 to 9
SSMIS	10	1 to 24	1 to 21
AIRS	11	1 to 2378	1 to 2378
HSB	12	1 to 4	1 to 4
MODIS	13	1 to 17	1 to 17
ATSR	14	1 to 3	1 to 3
MHS	15	1 to 5	1 to 5
IASI	16	1 to 8461	1 to 8461
AMSR	17	1 to 14	1 to 7
MVIRI	20	1 to 2	1 to 2
SEVIRI	21	4 to 11	1 to 8
GOES-Imager	22	1 to 4	1 to 4
GOES-Sounder	23	1 to 18	1 to 18
GMS/MTSAT imager	24	1 to 4	1 to 4
FY2-VISSR	25	1 to 2	1 to 2
FY1-MVISR	26	1 to 3	1 to 3
CriS	27	TBD	TBD
CMISS	28	TBD	TBD
VIIRS	29	TBD	TBD
WINDSAT	30	1 to 10	1 to 5

NAM_DIAG - RADAR simulator

Note

A radar simulator already existed in Meso-NH [Richard et al., 2003] that computes reflectivities in the Rayleigh approximation on each grid points of the model (NVERSION_RAD=1). However, with the view to code an observation operator for radar reflectivities, this simulator was not sufficient. That is why a new simulator was built, while the original version is still available. This new simulator (NVERSION_RAD=2) simulates Plan Position Indicators (PPI), which are cones usually projected on a horizontal plane obtained by scanning the atmosphere at constant elevation. New features are:

• possibility to choose among several scattering models,

• beam bending taken into account,

• possibility to take attenuation into account,

• antenna’s radiation pattern (beam broadening) modeled,

• ouptut on operational (Cartesian) grids of the Aramis French radar network.

• add LRADAR=T and NVERSION_RAD=1 in NAM_DIAG to store following variables :

Name	Meaning [Unit]	Dimension
RARE	Radar Reflectivity (dBz)	3D
VDOP	Radar Doppler fall speed (m/s)	3D
ZDR	Radar Differential Reflectivity (dBZ)	3D
KDP	Radar Differential Phase shift (degree/km)	3D

• add LRADAR=T and NVERSION_RAD=2 in NAM_DIAG to store following variables :

Note

As output fields are not on the model grid, they have to be written in specific files. Therefore the following files are written in the following format: AAABBBCC.CDDDX for cartesian coordinates and PAAABBBCC.CX for polar coordinates, where AAA is the descriptor of the field (3 characters, see below for further explanations), BBB is the name of the radar (3 characters), CC.C is the elevation (in degrees), DDD is half the number of pixels on each row or column (3 characters), and X is the name of the input file. Example of file name for cartesian coordinates : ZHHBOL00.4300BOG12.2.SEG04.004RD

Name	Meaning [Unit]	Dimension
ZHH	overall reflectivity (dBZ)
ZER	reflectivity due to rain (dBZ)
ZEI	reflectivity due to pristine ice (dBZ)
ZES	reflectivity due to snow (dBZ)
ZEG	reflectivity due to graupel (dBZ)
KDP	specific differential phase (km-1)
ZDR	differential reflectivity (dB)
VRU	Doppler velocity (m s-1)
HAS	height of middle of beam MSL (m)
M_R	rainwater contents in the middle of the beam (kg kg-1)
M_I	primary ice contents in the middle of the beam (kg kg-1)
M_S	snow contents in the middle of the beam (kg kg-1)
M_G	graupel contents in the middle of the beam (kg kg-1)
CIT	pristine ice concentration in the middle of the beam (kg m-3)
AET	overall two-way specific attenuation (dB km-1) (if LATT=T)
AER	two-way specific attenuation due to rain (dB km-1) (if LATT=T)
AEI	two-way specific attenuation due to pristine ice (dB km-1) (if LATT=T)
AES	two-way specific attenuation due to snow (dB km-1) (if LATT=T)
AEG	two-way specific attenuation due to graupel (dB km-1) (if LATT=T)
ATT	overall two-way path-integrated attenuation (PIA) (dB) (if LATT=T)
ATR	two-way PIA due to rain (dB) (if LATT=T)
ATI	two-way PIA due to pristine ice (dB) (if LATT=T)
ATS	two-way PIA due to snow (dB) (if LATT=T)
ATG	two-way PIA due to graupel (dB) (if LATT=T)
RFR	refractivity in the middle of the beam (if LREFR=T)
DNZ	vertical gradient of refractivity in the middle of the beam (km-1) (if LDNDZ=T)
CSR	index characterizing the pixel
PDP	differential phase beetween horizontal and vertical polarizations (◦)
KDR,KDS,KDG	specific differential phase due to rain, snow or graupel(km-1)
ZDA,ZDS,ZDG	differential reflectivity due to rain, snow or graupel (dB)
RHV,RHR,RHS,RHG	copolar correlation coefficient due to all hydrometeors, rain, snow or graupel (/)
TEM	model temperature (C)
DHV	backscattering differentiel phase ()

Note

Following options can be used in NAM_DIAG when LRADAR=T and NVERSION_RAD=2 is activated :

Fortran name	Fortran type	Default value
XLAT_RAD	array of reals	XUNDEF
XLON_RAD	array of reals	XUNDEF
XALT_RAD	array of reals	XUNDEF
CNAME_RAD	array of strings	XUNDEF
XLAM_RAD	array of reals	XUNDEF
XDT_RAD	array of reals	XUNDEF
XELEV	2-dim array of reals	XUNDEF
NBSTEPMAX	integer	-1
XSTEP_RAD	real	XUNDEF
LATT	logical	.FALSE.
LQUAD	logical	.FALSE.
NPTS_H	integer	1
NPTS_V	integer	1
CARF	character(5)	PB70
LREFR	logical	.FALSE.
LDNDZ	logical	.FALSE.
NCURV_INTERPOL	integer	0
LCART_RAD	integer	.TRUE.
NBAZIM	logical	.720
NDIFF	integer	0
NPTS_GAULAG	integer	7
XGRID	real	2000.0
LFALL	logical	.FALSE.
LWREFL	logical	.FALSE.
LWBSCS	logical	.FALSE.
XREFLMIN	real	-30.
XREFLVDOPMIN	real	-990.
LSNRT	logical	.TRUE.
XSNRMIN	real	0

• XLAT_RAD : latitude of each radar

• XLON_RAD : longitude of each radar

• XALT_RAD : altitudes of radars (m)

• CNAME_RAD : names of radars

• XLAM_RAD : radar wavelengths

• XDT_RAD : beam width to the -3 dB level for one-way transmission (\(\Delta\theta\))

• XELEV : radar elevations (\(\theta\)). First dimension: radar; second: site number

• NBSTEPMAX : number of gates

• XSTEP_RAD : gate length (m)

• LATT : attenuation is taken into account if true

• LQUAD : if true Gauss-Legendre quadrature if false Gauss-Hermite quadrature

• NPTS_H : number of angles for the quadrature in horizontal

• NPTS_V : number of angles for the quadrature in vertical

• CARF :

  • “PB70” : Pruppacher and Beard [1970].

  • “AND99” : axis ratio of raindrops : Andsager et al. [1999].

  • “BR02” : axis ratio of raindrops : Brandes et al. [2002].

  • “SPHE” : axis ratio for spheres (r=1)

• LREFR : if true writes out refractivity (\(N\equiv(n-1)\times10^6\))

• LDNDZ : if true writes out vertical gradient of refractivity (\(\partial N/\partial z\))

• NCURV_INTERPOL :

  • 0 : use an average beam bending equivalent to 4/3 of the Earth’s radius

  • 1 : compute the beam bending at each gate by using model variables

• LCART_RAD : if true interpolation of reflectivity on a cartesian grid ; false if polar

• NBAZIM : Number of azimuths in polar coordinates (used only if LCART_RAD=.FALSE)

• NDIFF :

  • 0 : Rayleigh scattering

  • 1 : Mie scattering

  • 3 : Rayleigh for spheroids scattering

  • 4 : Rayleigh with 6th order for attenuation calculations

  • 7 : T-matrix scattering (from lookup tables reading)

• NPTS_GAULAG : number of points of the quadrature

• XGRID : size of the Cartesian grid (m)

• LFALL : if true takes into account hydrometeor fall speeds

• LWREFL : if true takes into account the weighting by reflectivities

• LWBSCS : if true takes into account the weighting by hydrometeor concentrations

• XREFLMIN : minimum detectable reflectivity (in dBZ)

• XREFLVDOPMIN : minimum detectable reflectivity to compute Doppler velocities (in dBZ; useless when LWREFL=.FALSE.)

• LSNRT : if true ZHH ZER ZEI ZES ZEG and doppler velocity are thresholded when NR < XSNRMIN

• XSNRMIN : minimum SNR (used only if LSNRT=T)

NAM_DIAG - LIDAR simulator

• add LLIDAR=T in NAM_DIAG to store following variables :

Name	Meaning [Unit]	Dimension
LIDAR	Total backscatter coefficient (1/m/sr)	3D
LIPAR	Particle backscatter coefficient (1/m/sr)	3D

Note

Following options can be used in NAM_DIAG when LLIDAR is activated :

Fortran name	Fortran type	Default value
CVIEW_LIDAR	CHARACTER(LEN=5)	‘NADIR’
XALT_LIDAR	REAL	0
XWVL_LIDAR	REAL	0.532E-6

• CVIEW_LIDAR : gives the lidar point of view : either ‘NADIR’ or ‘ZENIT’

• XALT_LIDAR : gives the altitude of the lidar source (in meters) (by default, the altitude of the ground will be used for zenithal view, and the altitude of the model top will be used for nadir view)

• XWVL_LIDAR : gives the wavelength of the lidar source (in meters)

NAM_DIAG - Interpolation levels

Interpolation on altitude, isobaric and isentropic levels

Fortran name	Fortran type	Default value
LISOAL	LOGICAL	.FALSE.
XISOAL	ARRAY(REAL)	99*-1
LISOPR	LOGICAL	.FALSE.
XISOPR	ARRAY(REAL)	10*0
LISOTH	LOGICAL	.FALSE.
XISOTH	ARRAY(REAL)	10*0

• LISOAL : flag to interpolate on altitude levels the following variables: potential vorticity, wind, cloud (liquid water and ice) and precipitation (rain, snow and graupel) mixing ratio, dust extinction (if available). The outputs are 3D fields named: ALT_CLOUD, ALT_PRECIP, ALT_PRESSURE, ALT_PV, ALT_U, ALT_V and ALT_DSTEXT (if available).

• XISOAL : altitude of the isobaric levels

• LISOPR : flag to interpolate on pressure levels the following variables: potential temperature, wind, water vapour mixing ratio, geopotential (in meters). The outputs are 2D fields named with suffix ‘xxxxHPA’ where ‘xxxx’ stands for the pressure value

• XISOPR : altitude of the isobaric levels

• LISOTH : flag to interpolate on isentropic levels the following variables: pressure, potential vorticity, wind, water. The outputs are 2D fields named with suffix ‘xxxK’ where ‘xxx’ stands for the temperature value

• XISOTH : altitude of the isentropic levels

NAM_DIAG - Clustering

Fortran name	Fortran type	Default value
LCLSTR	LOGICAL	.FALSE.
LBOTUP	LOGICAL	.TRUE.
CFIELD	CHARACTER(LEN=8)	‘CLOUD’
XTHRES	REAL	0.00001

• LCLSTR : flag for 3D clustering

• LBOTUP : to propagate clustering from bottom to top (when TRUE); otherwise from top to bottom

• CFIELD : field on which clustering is applied, could be ‘W’ or ‘CLOUD’

• XTHRES : threshold value to detect the 3D structures

Following variables will be stored :

Name	Meaning [Unit]	Dimension
CLUSTERID	identity number	3D
CLUSTERLV	level where the object has been identified for the first time (at its bottom if LBOTUP is true, at its top otherwise)	3D
CLDSIZE	horizontal section of the object at the current level	3D

Note

Together, CLUSTERID and CLUSTERLV refers univoqually to a unique 3D object. Their value is homogeneous inside each identified object. CLOUDSIZE is homogeneous at each level inside each object.

NAM_DIAG - Coarse graining

Fortran name	Fortran type	Default value
LCOARSE	LOGICAL	.FALSE.
NDXCOARSE	INTEGER	1

• LCOARSE : flag to compute TKE (summation of the gridscale and the subgridscale parts) using coarse graining by both block and moving average

• NDXCOARSE : number of gridpoints over which the averaging is done

Following variables will be stored :

Name	Meaning [Unit]	Dimension
TKE_BLOCKAVGxx	TKE averaged block by block	3D
TKE_MOVINGAVGxx	TKE averaged over a moving block	3D

Note

The suffix xx stands for the number NDXCOARSE.

NAM_DIAG_BLANK

NAM_DIAG_BLANK content

Fortran name	Fortran type	Default value
XDUMMY_DIAG	ARRAY(REAL)	20* 0.0
NDUMMY_DIAG	ARRAY(INTEGER)	20* 0
LDUMMY_DIAG	ARRAY(LOGICAL)	20* .TRUE.
CDUMMY_DIAG	ARRAY(CHARACTER(LEN=80))	20* ‘’

Note

Similar use than NAM_BLANKn. Add USE MODD_DIAG_BLANK in a DIAG subroutine to use any of these variables.

NAM_DIAG_FILE

NAM_DIAG_FILE content

Fortran name	Fortran type	Default value
YINIFILE	ARRAY(CHARACTER(LEN=128))
YINIFILEPGD	ARRAY(CHARACTER(LEN=128))
YSUFFIX	ARRAY(CHARACTER(LEN=5))	_DIAG

• YINIFILE : name of the input synchronous backup files.

• YINIFILEPGD : name of the PGD file associated to YINIFILE.

• YSUFFIX : suffix appended to input file name to form output file name.

NAM_STO_FILE

Controls trajectories computation, only read if LTRAJ=.TRUE. in NAM_DIAG.

NAM_STO_FILE content

Fortran name	Fortran type	Default value
CFILES	ARRAY(CHARACTER(LEN=128)
NSTART_SUPP	ARRAY(INTEGER)	100*NUNDEF

• CFILES : name of all the input synchronous backup files used to compute trajectories. They must be in inverse chronological order, and correspond to a reinitialisation of Lagrangian tracers (see MESONH).

• NSTART_SUPP : extra origins for trajectory computations. In the second example below, the output files will contain the set of variables (X_TRAJ, Y_TRAJ, Z_TRAJ, THT_TRAJ, MRV_TRAJ).

Note

• Example 1 :

 &NAM_DIAG LVAR_LS=T, NCONV_KF=2, NRAD_3D=1,
           LVAR_MRW=T, LVAR_MRSV=T, LMOIST_V=T, LMOIST_E=F,
           LTPZH=T, LVORT=F, LMSLP=F, LGEO=T, LAGEO=T, LWIND_ZM=F,
           LTHW=T, LCLD_COV=T,
           LVAR_PR=F, LTOTAL_PR=F, LMEAN_PR=F, XMEAN_PR(1,2)=4. ,
           NCAPE=1, LRADAR=T, LTRAJ=F /

&NAM_DIAG_FILE YINIFILE(1) = "F9801.1.06A12.002" ,
               YINIFILEPGD(1) = "PGD_F9801" ,
               YSUFFIX = "diag" /

&NAM_DIAG_ISBAn N2M=2, LSURF_BUDGET=T /

• Example 2 : Namelist file for 6 files using trajectories computation

&NAM_DIAG LVAR_PR=T, LTOTAL_PR=T, LTPZH=T, LVAR_MRSV=T, LTRAJ=T /

&NAM_DIAG_FILE YSUFFIX='d18-6',
               YINIFILE(1) = "NAPE2.1.APE05.001" ,
               YINIFILEPGD(1) = "PGD_NAPE2" /

&NAM_STO_FILE CFILES(1) = "NAPE2.1.APE05.001" ,
              CFILES(2) = "NAPE2.1.APE04.001" ,
              CFILES(3) = "NAPE2.1.APE03.001" ,
              CFILES(4) = "NAPE2.1.APE02.001" ,
              CFILES(5) = "NAPE2.1.APE01.001" ,
              CFILES(6) = "APE10_ARP19990919.18" ,
              NSTART_SUPP(1)= 4 ,
              NSTART_SUPP(2)= 2 /

• Example 3 : Namelist file for simulator of radar To simulate the radar of Nancy, with T-matrix scattering, for 1 elevation (1.3)

&NAM_DIAG LRADAR=T,NVERSION_RAD=2,NCURV_INTERPOL=0,LCART_RAD=T,
          LQUAD=F,LWBSCS=T,LDNDZ=F, LFALL=F,LWREFL=F,LREFR=F,
          NPTS_GAULAG=7,NPTS_H=1,NPTS_V=1,CARF="AND99",
          NDIFF=0,NBSTEPMAX=400,XSTEP_RAD=700.,XGRID=2000.,LATT=F,
          XELEV(1,1)=01.3, XLAT_RAD(1)=48.7167,XLON_RAD(1)=6.5825,XALT_RAD(1)=297.55,
          CNAME_RAD(1)="NANCY",XLAM_RAD(1)=0.0535,XDT_RAD(1)=1.3 /

&NAM_DIAG_FILE YSUFFIX = "RD",
               YINIFILE(1) = "ALD00.2.SOG12.004",
               YINIFILEPGD(1) = "PGD_ALD00.2" /

NAM_CONF_DIAG

NAM_CONF_DIAG content

Fortran name	Fortran type	Default value
NHALO	INTEGER	1
JPHEXT	INTEGER	1

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

NAM_PARAM_KAFRn

It contains the options for deep and shallow convection parameterizations used by the model (CSCONV = “KAFR” or CDCONV = “KAFR” in NAM_PARAMn).

NAM_PARAM_KAFRN content

Fortran name	Fortran type	Default value
XDTCONV	REAL	MAX(300.0,XTSTEP)
NICE	INTEGER	1
LREFRESH_ALL	LOGICAL	.TRUE.
LCHTRANS	LOGICAL	.FALSE.
LDOWN	LOGICAL	.TRUE.
LSETTADJ	LOGICAL	.FALSE.
XTADJD	REAL	3600
XTADJS	REAL	10800
LDIAGCONV	LOGICAL	.FALSE.
NENSM	INTEGER	0

• XDTCONV : timestep for the call of the convective scheme. Maximum value is 300s.

• NICE : flag to include ice proceses in convection scheme ( 1 = yes, 0 = no ice ).

• LREFRESH_ALL : flag to refresh convective columns at every call of the convection scheme.

• LCHTRANS : flag to take into account the convective transport for scalar variables (chemical variables, passive pollutants, …). Can only be used with the options CDCONV=’KAFR’.

• LDOWN : flag to use downdrafts in deep convection.

• LSETTADJ : flag to allow user to define adjustment time.

• XTADJD : deep convective adjustment time (if LSETTADJ=TRUE).

• XTADJS : shallow convective adjustment time (if LSETTADJ=TRUE).

• LDIAGCONV : flag to store diagnostic variables in module MODD_DEEP_CONVECTIONn (CAPE, deep and shallow convective cloud top and base levels, up-and downdraft mass fluxes).

• NENSM : number of additional convective ensemble members for deep convection (for the moment limited to 3).

NAM_PARAM_RADn

It contains some options retained for the radiative scheme used by the model n CRAD = “ECMWF”; and some options common for both CRAD = “ECMWF” or “ECRA” in NAM_PARAMn)

NAM_PARAM_RADN content

Fortran name	Fortran type	Default value
XDTRAD	REAL	60.0
XDTRAD_CLONLY	REAL	60.0
NRAD_AGG	INTEGER	1
LCLEAR_SKY	LOGICAL	FALSE
NRAD_COLNBR	INTEGER	1000
NRAD_DIAG	INTEGER	0
LAERO_FT	LOGICAL	FALSE
LFIX_DAT	LOGICAL	FALSE
CLW	CHARACTER(LEN=4)	‘RRTM’
CAER	CHARACTER(LEN=4)	‘SURF’
CEFRADL	CHARACTER(LEN=4)	‘MART’
CEFRADI	CHARACTER(LEN=4)	‘LIOU’
COPWLW	CHARACTER(LEN=4)	‘SMSH’
COPILW	CHARACTER(LEN=4)	‘EBCU’
COPWSW	CHARACTER(LEN=4)	‘FOUQ’
COPISW	CHARACTER(LEN=4)	‘EBCU’
CAOP	CHARACTER(LEN=4)	‘CLIM’
XFUDG	REAL	1.0

Note

The following options are common for CRAD = “ECMWF” or “ECRA”:

• XDTRAD : Interval of time (in seconds) between two full radiation computations. ( the radiative tendency is computed for all verticals levels). This is done to save CPU time because the radiation scheme is very expensive and the radiative tendency is not evolving too much, in some cases, during periods greater than the model timestep XTSTEP. In this case, the “radiation timestep” is increased to XDTRAD

• XDTRAD_CLONLY : Interval of time (in seconds) between two radiation computations for the cloudy columns only. This is based on the same principle as the intermittent full radiation call: the cloudy column radiative tendency may, in some cases, evolve faster than the dry ones but still slower than the timestep XTSTEP. In this case, the “cloudy radiation timestep” is increased from XDTRAD to XDTRAD_CLONLY. Of course, when all and part of the radiative tendencies must be refreshed at the same MESONH timestep, only the full radiation call is performed.

• LCLEAR_SKY : When this flag is set to .TRUE., the radiative computations are made for a mean clear-sky and for the whole cloudy columns. This is still the way to spare some CPU time, by postulating that the clear sky columns do not lead to very different radiative tendencies. This hypothesis is only valid in academical cases.

• NRAD_COLNBR : Maximal number of air columns called by a single call of the radiative subroutine. This is performed in order to save memory, because the radiation subroutine allocate for every column of size NKMAX , NKMAX working arrays . This leads to a quadratic dependency of the memory with the number of vertical levels of the model. A way to limit the necessary memory is to split the number of columns passed to the radiation subroutine in several sets of NRAD_COLNBR column. Finally, all the desired columns (depending on the preceding parameters ) will be treated by sequentially calling the radiation subroutine for every set of column.

• NRAD_DIAG : number of diagnostic fields related to the radiative scheme stored in every output synchronous file (same fields as NRAD_3D in DIAG program). WARNING, a lot of variables are written only if NAM_DIAG_SURFn N2M=2.

• LAERO_FT : for a temporal interpolation of aerosol and ozone distribution. By default, they consist of monhtly averages kept constant for each month If true, the climatology of O3 and aerosols (only in TEGE case) are interpolated at each call of phys_paramnn. It is not usefull if your simulation lasts less than a month or does not contain any restard. It is necessary for long-term simulation with several segments to avoid too strong a perturbation at the beginning of each month.

• LFIX_DAT : flag to fix the date to a constant perpetual day. It is set by the initial SOUNDING date (RSOU). Note that the diurnal cycle is still considered.

• NRAD_AGG : side of a square of aggregated columns on which the radiation code will be called. This allows cheaper numerical cost of the radiation code and reduce its cost by \({NRAD\_AGG}^{2}\). If NRAD_AGG = 1, the radiation code is called on every columns (historical version). May be useful for very high resolution LES on which calling radiation on every columns is not necessary

Note

The following options are only used by CRAD = “ECMWF”

• CLW : choice of long wave radiative code

  • ‘RRTM’: RAPID RADIATIVE TRANSFER MODEL

  • ‘MORC’: MORCRETTE model

• CAER : type of aerosol distribution

  • ‘SURF’: deduced from cover data

  • ‘TEGE’: computed from Tegen et al. (1997) mensual climatology (horizontal resolution is 4 degrees of latitude by 5 degrees fo longitude

  • ‘TANR’: computed from ECMWF T5 climatology

  • ‘NONE’: no aerosol

• CEFRADL : liquid effective radius calculation

  • ‘MART’ : based on Martin et al. (1994, JAS)

  • ‘2MOM’ : based on the prediction of the number concentrations. Recommended with the 2-moment microphysical schemes.

  • ‘PRES’ : very old parametrization as f(pressure)

  • ‘OCLN’ : simple distinction between land (10) and ocean (13)

• CEFRADI : ice water effective radius calculation

  • ‘LIOU’ : ice particle effective radius =f(T) from Liou and Ou (1994)

  • ‘SURI’ : ice particle effective radius =f(T,IWC) from Sun and Rikus (1999)

  • ‘2MOM’ : based on the prediction of the number concentrations. Recommended with the 2-moment microphysical schemes (not yet available for mixed clouds).

  • ‘FX40’ : fixed 40 micron effective radius

• COPWLW : cloud water LW optical properties

  • ‘SMSH’: Smith-Shi formulation

  • ‘SAVI’: Savijarvi formulation (recommended only with 1-moment microphysical schemes with small precipitation)

  • ‘MALA’: Malavelle formulation (recommended only with 2-moment microphysical schemes with small precipitation)

• COPILW : ice water LW optical properties

  • ‘EBCU’: Ebert-Curry formulation

  • ‘SMSH’: Smith-Shi formulation, only with CLW=’RRTM’

  • ‘FULI’: Fu-Liou formulation, only with CLW=’MORC’

• COPWSW : cloud water short wave optical properties

  • ‘FOUQ’: Fouquart, 1991 formulation

  • ‘SLIN’: Slingo, 1989 formulation

  • ‘MALA’: Only for 2-moment microphysical schemes. According to Malavelle.

• COPISW : ice water short wave optical properties

  • ‘EBCU’: Ebert-Curry formulation

  • ‘FULI’: Fu-Liou formulation

• CAOP : type of aerosol optical properties calculation

  • ‘CLIM’: climatological aerosols

  • ‘EXPL’: explicit aerosols (if LORILAM=.T. in NAM_CH_ORILAM or LDUST=.T. in NAM_DUST)

• XFUDG : subgrid cloud inhomogeneity factor.

Note

The cloud overlap assumption for CRAD=’ECMW’ is defined in the routine ini_radconf.f90. The different assumptions are :

• NOVLP=5 : Random overlap for Clear Sky fraction and Effective Zenithal Angle. It is the best choice without subgrid condensation.

• NOVLP=6 : Maximum Random Overlap for Clear Sky fraction, and Random Overlap for Effective Zenithal Angle (DEFAULT VALUE). This option is well adapted to multi-layer clouds.

• NOVLP=7 : Maximum overlap for Clear Sky fraction and Random Overlap for Effective Zenithal Angle. This option is well adapted in the absence of multi-layer clouds.

• NOVLP=8 : Maximum Random overlap for Clear Sky fraction and Effective Zenithal Angle.

NAM_DIAG_SURF_ATMn

Warning

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

Diagnostics relative to each grid cell.

NAM_DIAG_SURF_ATMn content

Fortran name	Fortran type	Default value
LFRAC	LOGICAL	.FALSE.
LDIAG_GRID	LOGICAL	.FALSE.
LT2MMW	LOGICAL	.FALSE.
LDIAG_MIP	LOGICAL	.FALSE.

• LFRAC : flag to save in the output file the sea, inland water, natural covers and town fractions.

• LDIAG_GRID : flag for mean grid diagnostics

• LT2MMW : alternative weighting of grid average T2M giving more weight to the land tile.

• LDIAG_MIP : flag to perform intercomparison of land surface model diagnostics as required by several MIP (such as CMIP, SnowMIP, GCP, GSWP, etc.). These diag are only implemented for general surf atm diags, seaflux, Flake and ISBA.

NAM_WRITE_DIAG_SURFn

Warning

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

Diagnostics for to each grid cell and each tile.

NAM_WRITE_DIAG_SURFn content

Fortran name	Fortran type	Default value
LSELECT	LOGICAL	.FALSE.
CSELECT	ARRAY(CHARACTER)
LPROVAR_TO_DIAG	LOGICAL	.FALSE.
LSNOWDIMNC	LOGICAL	.FALSE.
LRESETCUMUL	LOGICAL	.FALSE.

• LSELECT : if T it indicates that a selection will be used as output.

• CSELECT : array containing the list of output fields.

• LPROVAR_TO_DIAG : used to write out prognostic variables like diagnostic one, on average over all patches.

• LSNOWDIMNC : in case of OFFLIN output files, to write the output snow fields in 2D (number of points / number of snow layers).

• LRESETCUMUL : in OFFLINE mode, for the ISBA scheme, replaces the instantaneous fields by the averaged cumulated fields on the output writing time step. Then the cumulated fields are cumulated during the writing time steps and reset at the end of each of them.

NAM_DIAG_SURFn

Warning

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

Diagnostics for each grid cell and each tile.

NAM_DIAG_SURFN content

Fortran name	Fortran type	Default value
N2M	INTEGER	2
LSURF_BUDGET	LOGICAL	.FALSE.
LSURF_BUDGETC	LOGICAL	.FALSE.
LRESET_BUDGETC	LOGICAL	.FALSE.
LRAD_BUDGET	LOGICAL	.FALSE.
LCOEF	LOGICAL	.FALSE.
LSURF_VARS	LOGICAL	.FALSE.
L2M_MIN_ZS	LOGICAL	.FALSE.

• N2M : flag to compute surface boundary layer characteristics:

  • N2M=2 : computes temperature at 2 m, specific humidity at 2 m, relative humidity, zonal and meridian wind at 10 m, and Richardson number. 2m and 10m quantities are calculated interpolating between atmospheric forcing variables and surface temperature and humidity.

• LSURF_BUDGET : flag to save in the output file the terms of the surface energy balance (net radiation, sensible heat flux, latent heat flux, ground flux), for each scheme (on the four separate tiles), on each patch of the vegetation scheme if existing, and aggregated for the whole surface. The diagnosed fields are (* stands for the scheme considered (* = nothing: gfield aggregated on the whole surface; * = name of a scheme : field for this scheme):

  • RN_*: net radiation

  • H_*: turbulent sensible heat flux

  • LE_*: turbulent latent heat flux

  • GFLUX_*: round or storage heat flux

  • FMU_*: zonal wind stress

  • FMV_*: meridian wind stress

Note

If both LSURF_BUDGET and LRAD_BUDGET are T then downward and upward short-wave radiation per spectral band will be written into output file (they are computed even if LRAD_BUDGET is False). The following output fields are then available:

• SWD_*: downward short wave radiation

• SWU_*: upward short wave radiation

• SWBD_*: downward short wave radiation for each spectral band

• SWBU_*: upward short wave radiation for each spectral band

• LWD_*: downward long wave radiation

• LWU_*: upward long wave radiation

• LSURF_BUDGETC : flag to save in the output file the time integrated values of all budget terms that have been activated

• LRESET_BUDGETC : flag to reset cumulatives variables at the beginning of a run

• LCOEF : flag to save in the output file the transfer coefficients used in the computation of the surface energy fluxes, for each scheme (on the four separate tiles) and aggregated for the whole surface. The diagnosed fields are (* stands for the scheme considered * = nothing: field aggregated on the whole surface; * = name of a scheme : field for this scheme):

  • CD_*: gdrag coefficient for momentum

  • CH_*: gdrag coefficient for heat

  • CE_*: gdrag coefficient for evaporation (differs from CH only over sea)

  • Z0_*: groughness length

  • Z0H_*: gthermal roughness length

• LSURF_VARS : flag to save in the output file the surface specific humidity for each scheme (on the four separate tiles), on each patch of the vegetation scheme if existing. The diagnosed fields are (* stands for the scheme considered (* = nothing: gfield aggregated on the whole surface; * = name of a scheme :< field for this scheme):

  • QS_*: specific humidity

• L2M_MIN_ZS : flag for 2 meters quantities evaluated on the minimum orography of the grid

NAM_DIAG_ISBAn

Warning

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

ISBA diagnostics.

NAM_DIAG_ISBAn content

Fortran name	Fortran type	Default value
LPGD	LOGICAL	.FALSE.
LSURF_EVAP_BUDGET	LOGICAL	.FALSE.
LSURF_MISC_BUDGET	LOGICAL	.FALSE.
LSURF_DIAG_ALBEDO	LOGICAL	.FALSE.
LPATCH_BUDGET	LOGICAL	.TRUE.
LSURF_MISC_DIF	LOGICAL	.FALSE.
LWATER_BUDGET	LOGICAL	.FALSE.
LLUTILES_BUDGET	LOGICAL	.FALSE.
LPROSNOW	LOGICAL	.FALSE.
LPROBANDS	LOGICAL	.FALSE.
LVOLUMETRIC_SNOWLIQ	LOGICAL	.FALSE.
LUTCI	LOGICAL	.FALSE.

• LPGD : flag to save in the output file the physiographic fields of ISBA scheme that are computed from ecoclimap data from the ecosystem fractions.

• LSURF_EVAP_BUDGET : flag to save in the output file the detailed terms of the water vapor fluxes, on each patch of the vegetation scheme if existing, and aggregated for the natural surface. The diagnosed fields are:

  • GPP: Gross primary production

• LSURF_MISC_BUDGET : flag to save in the output file miscelleaneous fields. The diagnosed fields are:

  • HV: Halstead coefficient

  • SNG: snow fraction over bare ground

  • SNV: snow fraction over vegetation

  • SN: total snow fraction

  • SWI: soil wetness index for each ground layer (wg - wwilt)/(wfc - wwilt) where wg is the volumic water content, wfc is the porosity and wwilt corresponds to the plant wilting point.

• LSURF_DIAG_ALBEDO : to write ALB…_ISBA et ALB…_S.

• LPATCH_BUDGET : flag to save in the output file the diagnostics for each patch (default is .T.)

• LSURF_MISC_DIF : to calculate and write specific DIF diagnostics

• LWATER_BUDGET : to calculate and write the water budget

• LLUTILES_BUDGET : flag to bring together diag from the ISBA patches into 4 surface covers type required for land-use-land-cover purpose (not implemented for ECOCLIMAP-SG) :

  • 1. Primary and secondary natural land (Forest, grassland, bare ground, etc.)

  • 2. Cropland (Agriculture)

  • 3. Pastureland (not yet implemented in ISBA)

  • 4. Urban settlement (not yet implemented; should be implemented if TEB is used).

• LPROSNOW : add new diagnostic fields for the CROCUS snow scheme, reproject the snow mantel and other diagnostic fields on the vertical, according to the subgrid slope, and merges ISBA_PROGNOSTIC.OUT.nc and ISBA_DIAGNOSTICS.OUT.nc in ISBA_PROGNOSTIC.OUT.nc, in case of CTIMESERIES_FILETYPE = ‘OFFLIN’.

• LPROBANDS : enable the spectral resolution of Crocus diagnostics, necessary if you want to get spectral albedo and spectral direct/diffuse ratio diagnostics

• LVOLUMETRIC_SNOWLIQ: convert the SNOWLIQ diagnostic field in kg / m3 (instead of m).

• LUTCI : flag to compute UTCI (human thermal comfort indicator) quantities in rural areas

NAM_DIAG_TEBn

Warning

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

TEB diagnostics.

NAM_DIAG_TEBn content

Fortran name	Fortran type	Default value
LPGD	LOGICAL	.FALSE.
LSURF_MISC_BUDGET	LOGICAL	.FALSE.
LSURF_DIAG_ALBEDO	LOGICAL	.FALSE.
LUTCI	LOGICAL	.FALSE.

• LPGD : flag to save PGD fields if TEB garden is activated

• LSURF_MISC_BUDGET : flag to save in the output file miscelleaneous fields. The diagnosed fields are:

  • Z0_TOWN: roughness length for town

  • QF_BLD: domestic heating

  • QF_BLDWFR: domestic heating

  • FLX_BLD: heat flux from bld

  • TI_BLD_EQ: internal temperature without heating

  • TI_BLDWFR: internal temperature without heating

  • QF_TOWN: total anthropogenic heat

  • DQS_TOWN: storage inside building

  • H_WALL: wall sensible heat flux

  • H_ROOF: roof sensible heat flux

  • H_ROAD: road sensible heat flux

  • RN_WALL: net radiation at wall

  • RN_ROOF: net radiation at roof

  • RN_ROAD: net radiation at road

  • GFLUX_WALL: net wall conduction flux

  • GFLUX_ROOF: net roof conduction flux

  • GFLUX_ROAD: net road conduction flux

  • LE_ROOF: roof latent heat flux

  • LE_ROAD: road latent heat flux

• LSURF_DIAG_ALBEDO : flag to save in the output file albedo diagnostics

• LUTCI : to calculate and write UTCI diagnostics

NAM_DIAG_FLAKEn

NAM_DIAG_FLAKEn content

Fortran name	Fortran type	Default value
LWATER_PROFILE	LOGICAL	F
XZWAT_PROFILE	REAL
LSEDIM_PROFILE	LOGICAL	F
XZSED_PROFILE	REAL
LFLKFLUX	LOGICAL	F
LFLKWATER	LOGICAL	F

• LWATER_PROFILE : flag to save in the output file miscelleaneous fields. The diagnostic is temperature at the depths defined by:

• XZWAT _PROFILE : depth of output levels (m) in namelist

• LSEDIM_PROFILE : flag for sediment diagnostics

• XZSED_PROFILE : depth of output levels (m) in namelist

• LFLKFLUX : flag for heat and radiative diagnostics

• LFLKWATER : flag for water budget P-E diagnostics

NAM_DIAG_OCEANn

Warning

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

NAM_DIAG_OCEANn content

Fortran name	Fortran type	Default value
LDIAG_OCEAN	LOGICAL	.FALSE.

• LDIAG_OCEAN : flag for ocean variables