namelists: GEM 5.0

switch:.T.:MONO(CLIPPING) after Bermejo-Conde (MASS-CONSERVATION for Chemical Tracers)

switch:.T.:

East boundary in GY for an embedded LAM (MASS-CONSERVATION for Chemical Tracers)

North boundary in GY for an embedded LAM (MASS-CONSERVATION for Chemical Tracers)

South boundary in GY for an embedded LAM (MASS-CONSERVATION for Chemical Tracers)

West boundary in GY for an embedded LAM (MASS-CONSERVATION for Chemical Tracers)

switch:.T.:MONO(CLIPPING) of RHS

Scaling for mass of tracer

• adv_scaling=0 (none)
• adv_scaling=1 (CO2)
• adv_scaling=2 (O3)

Type of rebuild in SLICE

• SLICE_rebuild=1 (LP)
• SLICE_rebuild=2 (CW)

(MASS-CONSERVATION for Chemical Tracers)

Use surface layer winds for advection of lowest thermodynamic level

Activate conservation diagnostics if /=0 (MASS-CONSERVATION for Chemical Tracers)

Main selector for dynamical kernel

• "DYNAMICS_FISL_P" : Fully implicit SL in pressure
• "DYNAMICS_FISL_H" : Fully implicit SL in height
• "DYNAMICS_EXPO_H" : Exponential integrators in height

Switch to activate generation of Markov chains, use of SKEB and use of PTP

Switch to activate the first initialisation of Matkov chaines in PTP & SKEB

Switch to activate the in PTP & SKEB

Seed of the random number generator usually we put DAY and member number (3D MARKOV CHAINES)

CAPE value in Kain-Fritsch scheme to stop perturbing the physical tendencies (2D MARKOV CHAINES)

switch to activate PTP (perturb tendencies of physics) (2D MARKOV CHAINES)

vertical velocity value (m/s) above which we stop perturbing the physical tendencies (2D MARKOV CHAINES)

bottom value of transition zone of vertical envelope in sigma for PTP (below that no perturbation) (2D MARKOV CHAINES)

upper value of transition zone of vertical envelope in sigma for PTP (above that full perturbation) (2D MARKOV CHAINES)

factor of reduction of the perturbation the physical tendencies (in PTP) when convection occurs (2D MARKOV CHAINES)

(ignored) Ens_ptp_l = Ens_ptp_trnh-Ens_ptp_trnl+1 (2D MARKOV CHAINES)

(ignored) Ens_ptp_lmax = maxval(Ens_ptp_l) (2D MARKOV CHAINES)

(ignored) Ens_ptp_m = Ens_ptp_trnh+1 (2D MARKOV CHAINES)

maximum value of the 2D Markov chains (2D MARKOV CHAINES)

minimum value of the 2D Markov chain (2D MARKOV CHAINES)

(ignored) Ens_ptp_mmax = maxval(Ens_ptp_m) (2D MARKOV CHAINES)

number of 2d Markov chains (2D MARKOV CHAINES)

no. of latitudes for 2D Markov chains (2D MARKOV CHAINES)

no. of longitudes for 2D Markov chains (2D MARKOV CHAINES)

standard deviation value for 2D Markov chains (2D MARKOV CHAINES)

value of stretch for Markov chains (2D MARKOV CHAINES)

decorrelation time (seconds) for 2D Markov chains (2D MARKOV CHAINES)

TLC value (convective precipitation) in Kuo (OLDKUO) scheme to stop perturbing the physical tendencies (2D MARKOV CHAINES)

high wave number horizontal truncation limit for 2D Markov chains (2D MARKOV CHAINES)

low wave number horizontal truncation limit for 2D Markov chains (2D MARKOV CHAINES)

coefficient Alpha for momentum in SKEB (3D MARKOV CHAINES)

coefficient Alpha for temperature in SKEB (3D MARKOV CHAINES)

coefficient for Gaussian filter used in SKEB (3D MARKOV CHAINES)

Switch to activate SKEB (3D MARKOV CHAINES)

switch to do SKEB calculation based on diffusion (3D MARKOV CHAINES)

switch to do the calculation of the divergence due to SKEB forcing (3D MARKOV CHAINES)

switch to do SKEB calculation based on gravity wave drag (3D MARKOV CHAINES)

wavelength for Gaussian filter in SKEB (3D MARKOV CHAINES)

maximum value of the 3D Markov chain (used by SKEB) (3D MARKOV CHAINES)

minimum value of the 3D Markov chain (used by SKEB) (3D MARKOV CHAINES)

(3D MARKOV CHAINES)

number of latitudes of the gaussian grid used for the 3D Markov chains (used in the SKEB calculation) (3D MARKOV CHAINES)

number of longitudes of the gaussian grid used for the 3D Markov chains (in the SKEB calculation) (3D MARKOV CHAINES)

std. dev. value for the 3D Markov chain (used by SKEB) (3D MARKOV CHAINES)

value of stretch for 3D Markov chain (used by SKEB) (3D MARKOV CHAINES)

decorrelation time (seconds) for 3D Markov chain (used by SKEB) (3D MARKOV CHAINES)

high wave number truncation limit used in 3D Markov chain (used by SKEB) (3D MARKOV CHAINES)

low wave number truncation limit used in 3D Markov chain (used by SKEB) (3D MARKOV CHAINES)

switch to print global stat related to Markov chains, SKEB and PTP (3D MARKOV CHAINES)

SL off-centering parameter for hydrostatic

SL off-centering parameter for the momentum equations

SL off-centering parameter for nonhydrostatic

Fraction of adjustment to be given to the ocean

another reference pressure

Parameter controlling modified epsilon (Ver_epsi_8) [nonhydrostatic part]

T* basic state temperature (K)

Latitudinal ramping of equatorial sponge

Coefficients that multiply KM to simulate sponge layer near the top of the model. Warning! if this parameter is used, the EPONGE in the physics namelist should be removed.

number of points for the halo on X

number of points for the halo on Y

Heap memory will be painted to NaN using an array wrk01(G_ni,G_nj,Heap_nk)

array of model levels , 0.0 < HYB < 1.0

pair of coefficients (min,max) to control the flattenning of the vertical coordinate

Fraction of the maximum divergence damping - range(0.0-1.0)

Background 2 delta-x removal ratio - range(0.0-1.0)

Theta 2 delta-x removal ratio - range(0.0-1.0).

Tracers 2 delta-x removal ratio - range(0.0-1.0)

Order of the background diffusion operator 2, 4, 6, 8

Order of the background diffusion operator on theta 2, 4, 6, 8

Order of the background diffusion operator on tracers

Frictional heating is considered when Hzd_smago_fric_heat>0.

The levels (bot,top) in the hybrid coordinate where the background diffusion coefficient varies between the value defined by Hzd_smago_lnr(1) and Hzd_smago_lnr(2).

Coefficient of background diffusion added to the coefficient computed using the Smagorinsky approach. The first element of the array determines the constant value of background diffusion coeff. below Hzd_smago_lev(1). The second element represents the value at Hzd_smago_lev(2). The third element determines the maximum coefficient at the model top. Two ramps of COS^2-type are used between Hzd_smago_lev(1) and Hzd_smago_lev (2), and between Hzd_smago_lev(2) and the model lid.

Main Smagorinsky control parameter (usual range 0.1-0.3)

Apply Smago diffusion on theta using Hzd_smago_param/Hzd_smago_prandtl parameter

Apply Smago diffusion on HU using Hzd_smago_param/Hzd_smago_prandtl_hu parameter

If TRUE then background diffusion is applied to THETA and HU.

Filter cutoff period for Iau_weight_S='sin' in hours

IAU Input TYPE

• 'OLD ' : GEM 4.8 physic's input, MPI blocking based
• 'IO ' : new input system, RPN_COMM_IO/RPN_COMM_ezshuf_dist based
• 'BLOC' : new input system, RPN_COMM_bloc based

The number of seconds between increment fields

IAU Input PE blocking along npex

IAU Input PE blocking along npey

The number of seconds over which IAU will be will be run (typically the length of the assimilation window). Default < 0 means that no IAUs are applied.

IAU Input Stats

An optional list of tracers to be incremented.

The type of weighting function to be applied to the analysis increments:

• 'constant' (default) uniform increments
• 'sin' DF-style weights (Fillion et al. 1995)

true -> Digital filter initialization is performed

number of points for digital filter (equals the number of timesteps +1)

period limit of digital filter units D,H,M,S

• true -> passive tracers digitally filtered
• false-> passive tracers set to result obtained at mid-period during initialization period (no filtering)

true -> Windowing is applied

List of variables to NOT process during input

Number of PEs to use for input

Type of vertical interpolation scheme

True-> for blending to zero the physics tendency in blending area

Number of points for horizontal blending

Number of levels for top blending

True-> to blend the model topography with the pilot topography

True-> to force constant (fixed) boundary conditions

Number of levels for top piloting

Type of horizontal interpolation to model grid

• 'CUB_LAG'
• 'LINEAR'
• 'NEAREST'

True-> The plane of the top temperature layer is completely overwritten from the 2D pilot data

True->to check for time left in job

True->to print more information to std output

precision in print glbstats

• 'LCL_4'
• 'GLB_8'

True-> to clip humidity variables on output

Interval of output file name change

'etiket' used for output fields

Default value for IP3 is 0, -1 for IP3 to contain step number, >0 for given IP3

Precision criteria for the Liebman procedure

List of levels for underground extrapolation

Maximum number of iterations for the Liebman procedure

number of iterations to exchange halo for the Liebman procedure

Number of layers close to the bottom of the model within which a linear interpolation of GZ will be performed

Packing factor used for all variables except for those defined in Out_xnbits_s

Minimum of digits used to represent output units

Total number of PEs for output using MFV collector

Total number of PEs along npex for output using MID collector

Total number of PEs along npey for output using MID collector

Sortie jobs lauched every Out3_postproc_fact*Out3_close_interval_S

Vertical interpolation scheme for output

Latitude at which the multiplication factor becomes P_lmvd_weigh_high_lat

latitude at which the multiplication factor becomes P_lmvd_weigh_low_lat

Multiplication factor of P_pbl_spng at latitude P_lmvd_high_lat

Multiplication factor of P_pbl_spng at latitude P_lmvd_low_lat

Number of bits to perturb on initial conditions

Stride for perturbation on initial conditions

• True-> cubic interpolation
• False-> linear interpolation

• 0 -> NO advection
• 1 -> traditional advection
• 2 -> consistent advection with respect to off-centering
• 3 -> reversed advection with respect to off-centering

True-> auto barotropic option

True-> averaging B and C in SLEVE scheme

True-> Modify slightly code behaviour to ensure bitpattern reproduction in restart mode using FST file

True-> variable cappa in thermodynamic equation

Use cubic interpolation in trajectory computation

True->

True-> Eulerian treatment of mountains in the continuity equation

• True-> hydrostatic
• False-> non-hydrostatic

True-> horizontal diffusion of momentum at each CN iteration

Number of iterations for Crank-Nicholson

Number of iterations to solve non-linear Helmholtz problem

Number of iterations to compute trajectories

• -1: no blending between Yin and Yang
• 0: blending at init only
• > 0: blending at every nblendyy timestep

Physics coupling strategy

• 0 -> No conservation of surface pressure
• 1 -> Conservation of Total air mass Pressure
• 2 -> Conservation of Dry air mass Pressure

Confirmation to use psadjust with a LAM configuration

True-> print dry/wet air masses

True-> to use topography

Use trapezoidal average for advection winds

Apply water loading in the calculations

2D preconditioner for iterative solver

Krylov method for 3d iterative solver (FGMRES or FBICGSTAB)

3D preconditioner for iterative solver

• True-> use FFT solver if possible
• False-> use MXMA solver (slower,less precise)

Epsilon convergence criteria for none Yin-Yang iterative solver

maximum number of iterations allowed for none Yin-Yang iterative solver

size of Krylov subspace in iterative solver - should not exceed 100

Type of solver

• 'ITERATIF'
• 'DIRECT'

Epsilon convergence criteria for the Yin-Yang iterative solver

maximum number of iterations allowed for the Yin-Yang iterative solver

The filter will be set zero for smaller scales (in km)

The filter will be set 1.0 for larger scales (in km) between Spn_cutoff_scale_small and Spn_cutoff_scale_large, the filter will have a COS2 transition.

Spectral nudging list of variables (eg. 'UVT' or 'UV')

Nudging relaxation timescale (eg. 10 hours )

Nudging profile lower end in hyb level (eg. 1.0 or 0.8) If use 0.8, the profile will be set zero when hyb > 0.8

Nudging interval in seconds (eg. 1800, means nudging is performed every every 30 minutes)

Nudging profile transition shape('COS2' or 'LINEAR') Set the shape between Spn_start_lev and Spn_up_const_lev

Nudging profile upper end in hyb level (eg. 0.0 or 0.2) If use 0.2, the profile wll be set 1.0 when hyb < 0.2

Nudging weight in temporal space (.true. or .false.). If the driving fields are available every 6 hours and Spn_step is set to 30 minutes then nudging will have more weight every six hours when the driving fields are available

The weight factor when Spn_weight_L=.true. (The weigh factor is COS2**(Spn_wt_pwr), Spn_wt_pwr could be set as 0, 2, 4, 6. If Spn_wt_pwr = 2, weight factor is COS2)

list of variables to do blocstat. Any gmm variable name, or predefine lists :

• 'ALL'
• 'ALL_DYN_T0'
• 'ALL_DYN_T1'
• 'ALL_TR_T0'
• 'ALL_TR_T1'

True-> tracers validity time does not have to match analysis

Override for default tracers attributes

list of tracers to be read from analyse

Top coefficient for del-2 diffusion (m2/s)

Number of levels from the top of the model

True-> Riley diffusion on vertical motion on Vspng_nk levels

On which length of time to evolve topography

Time at which to start evolving topography toward target

True-> divergence high level modulation in initial computation of Zdot

SL off-centering parameter for hydrostatic

SL off-centering parameter for the momentum equations

SL off-centering parameter for nonhydrostatic

Fraction of adjustment to be given to the ocean

another reference pressure

Parameter controlling modified epsilon (Ver_epsi_8) [nonhydrostatic part]

T* basic state temperature (K)

Latitudinal ramping of equatorial sponge

Coefficients that multiply KM to simulate sponge layer near the top of the model. Warning! if this parameter is used, the EPONGE in the physics namelist should be removed.

number of points for the halo on X

number of points for the halo on Y

Heap memory will be painted to NaN using an array wrk01(G_ni,G_nj,Heap_nk)

array of model levels , 0.0 < HYB < 1.0

pair of coefficients (min,max) to control the flattenning of the vertical coordinate

Fraction of the maximum divergence damping - range(0.0-1.0)

Background 2 delta-x removal ratio - range(0.0-1.0)

Theta 2 delta-x removal ratio - range(0.0-1.0).

Tracers 2 delta-x removal ratio - range(0.0-1.0)

Order of the background diffusion operator 2, 4, 6, 8

Order of the background diffusion operator on theta 2, 4, 6, 8

Order of the background diffusion operator on tracers

Frictional heating is considered when Hzd_smago_fric_heat>0.

The levels (bot,top) in the hybrid coordinate where the background diffusion coefficient varies between the value defined by Hzd_smago_lnr(1) and Hzd_smago_lnr(2).

Coefficient of background diffusion added to the coefficient computed using the Smagorinsky approach. The first element of the array determines the constant value of background diffusion coeff. below Hzd_smago_lev(1). The second element represents the value at Hzd_smago_lev(2). The third element determines the maximum coefficient at the model top. Two ramps of COS^2-type are used between Hzd_smago_lev(1) and Hzd_smago_lev (2), and between Hzd_smago_lev(2) and the model lid.

Main Smagorinsky control parameter (usual range 0.1-0.3)

Apply Smago diffusion on theta using Hzd_smago_param/Hzd_smago_prandtl parameter

Apply Smago diffusion on HU using Hzd_smago_param/Hzd_smago_prandtl_hu parameter

If TRUE then background diffusion is applied to THETA and HU.

Filter cutoff period for Iau_weight_S='sin' in hours

IAU Input TYPE

• 'OLD ' : GEM 4.8 physic's input, MPI blocking based
• 'IO ' : new input system, RPN_COMM_IO/RPN_COMM_ezshuf_dist based
• 'BLOC' : new input system, RPN_COMM_bloc based

The number of seconds between increment fields

IAU Input PE blocking along npex

IAU Input PE blocking along npey

The number of seconds over which IAU will be will be run (typically the length of the assimilation window). Default < 0 means that no IAUs are applied.

IAU Input Stats

An optional list of tracers to be incremented.

The type of weighting function to be applied to the analysis increments:

• 'constant' (default) uniform increments
• 'sin' DF-style weights (Fillion et al. 1995)

true -> Digital filter initialization is performed

number of points for digital filter (equals the number of timesteps +1)

period limit of digital filter units D,H,M,S

• true -> passive tracers digitally filtered
• false-> passive tracers set to result obtained at mid-period during initialization period (no filtering)

true -> Windowing is applied

List of variables to NOT process during input

Number of PEs to use for input

Type of vertical interpolation scheme

True-> for blending to zero the physics tendency in blending area

Number of points for horizontal blending

Number of levels for top blending

True-> to blend the model topography with the pilot topography

True-> to force constant (fixed) boundary conditions

Number of levels for top piloting

Type of horizontal interpolation to model grid

• 'CUB_LAG'
• 'LINEAR'
• 'NEAREST'

True-> The plane of the top temperature layer is completely overwritten from the 2D pilot data

True->to check for time left in job

True->to print more information to std output

precision in print glbstats

• 'LCL_4'
• 'GLB_8'

True-> to clip humidity variables on output

Interval of output file name change

'etiket' used for output fields

Default value for IP3 is 0, -1 for IP3 to contain step number, >0 for given IP3

Precision criteria for the Liebman procedure

List of levels for underground extrapolation

Maximum number of iterations for the Liebman procedure

number of iterations to exchange halo for the Liebman procedure

Number of layers close to the bottom of the model within which a linear interpolation of GZ will be performed

Packing factor used for all variables except for those defined in Out_xnbits_s

Minimum of digits used to represent output units

Total number of PEs for output using MFV collector

Total number of PEs along npex for output using MID collector

Total number of PEs along npey for output using MID collector

Sortie jobs lauched every Out3_postproc_fact*Out3_close_interval_S

Vertical interpolation scheme for output

Latitude at which the multiplication factor becomes P_lmvd_weigh_high_lat

latitude at which the multiplication factor becomes P_lmvd_weigh_low_lat

Multiplication factor of P_pbl_spng at latitude P_lmvd_high_lat

Multiplication factor of P_pbl_spng at latitude P_lmvd_low_lat

Number of bits to perturb on initial conditions

Stride for perturbation on initial conditions

• True-> cubic interpolation
• False-> linear interpolation

• 0 -> NO advection
• 1 -> traditional advection
• 2 -> consistent advection with respect to off-centering
• 3 -> reversed advection with respect to off-centering

True-> auto barotropic option

True-> averaging B and C in SLEVE scheme

True-> Modify slightly code behaviour to ensure bitpattern reproduction in restart mode using FST file

True-> variable cappa in thermodynamic equation

Use cubic interpolation in trajectory computation

True->

True-> Eulerian treatment of mountains in the continuity equation

• True-> hydrostatic
• False-> non-hydrostatic

True-> horizontal diffusion of momentum at each CN iteration

Number of iterations for Crank-Nicholson

Number of iterations to solve non-linear Helmholtz problem

Number of iterations to compute trajectories

• -1: no blending between Yin and Yang
• 0: blending at init only
• > 0: blending at every nblendyy timestep

Physics coupling strategy

• 0 -> No conservation of surface pressure
• 1 -> Conservation of Total air mass Pressure
• 2 -> Conservation of Dry air mass Pressure

Confirmation to use psadjust with a LAM configuration

True-> print dry/wet air masses

True-> to use topography

Use trapezoidal average for advection winds

Apply water loading in the calculations

2D preconditioner for iterative solver

Krylov method for 3d iterative solver (FGMRES or FBICGSTAB)

3D preconditioner for iterative solver

• True-> use FFT solver if possible
• False-> use MXMA solver (slower,less precise)

Epsilon convergence criteria for none Yin-Yang iterative solver

maximum number of iterations allowed for none Yin-Yang iterative solver

size of Krylov subspace in iterative solver - should not exceed 100

Type of solver

• 'ITERATIF'
• 'DIRECT'

Epsilon convergence criteria for the Yin-Yang iterative solver

maximum number of iterations allowed for the Yin-Yang iterative solver

The filter will be set zero for smaller scales (in km)

The filter will be set 1.0 for larger scales (in km) between Spn_cutoff_scale_small and Spn_cutoff_scale_large, the filter will have a COS2 transition.

Spectral nudging list of variables (eg. 'UVT' or 'UV')

Nudging relaxation timescale (eg. 10 hours )

Nudging profile lower end in hyb level (eg. 1.0 or 0.8) If use 0.8, the profile will be set zero when hyb > 0.8

Nudging interval in seconds (eg. 1800, means nudging is performed every every 30 minutes)

Nudging profile transition shape('COS2' or 'LINEAR') Set the shape between Spn_start_lev and Spn_up_const_lev

Nudging profile upper end in hyb level (eg. 0.0 or 0.2) If use 0.2, the profile wll be set 1.0 when hyb < 0.2

Nudging weight in temporal space (.true. or .false.). If the driving fields are available every 6 hours and Spn_step is set to 30 minutes then nudging will have more weight every six hours when the driving fields are available

The weight factor when Spn_weight_L=.true. (The weigh factor is COS2**(Spn_wt_pwr), Spn_wt_pwr could be set as 0, 2, 4, 6. If Spn_wt_pwr = 2, weight factor is COS2)

list of variables to do blocstat. Any gmm variable name, or predefine lists :

• 'ALL'
• 'ALL_DYN_T0'
• 'ALL_DYN_T1'
• 'ALL_TR_T0'
• 'ALL_TR_T1'

True-> tracers validity time does not have to match analysis

Override for default tracers attributes

list of tracers to be read from analyse

Top coefficient for del-2 diffusion (m2/s)

Number of levels from the top of the model

True-> Riley diffusion on vertical motion on Vspng_nk levels

On which length of time to evolve topography

Time at which to start evolving topography toward target

True-> divergence high level modulation in initial computation of Zdot

x horizontal resolution of target cascade grid (degrees)

y horizontal resolution of target cascade grid (degrees)

Time string (units D, H, M or S) from the start of the run to stop producing the cascade files

Number of points for the blending zone (Hblen_x)

TRUE to dump out permanent bus for cascade mode

Reference Point I on rotated cascade grid, 1 < Grdc_iref < Grdc_ni

Reference Point J on rotated cascade grid, 1 < Grdc_jref < Grdc_nj

Latitude on rotated grid of ref point, Grdc_iref,Grdc_jref (degrees)

Longitude on rotated grid of ref point, Grdc_iref,Grdc_jref (degrees)

Max Supported Courrant number; Pilot area=Grdc_maxcfl +Grdc_bsc_base+Grdc_bsc_ext1

Number of bits for the packing factor

Nesting interval specified with digits ending with one character for the units:

• S : seconds
• D : days
• M : minutes
• H : hours

Number of points along X

Number of points along Y

Time string (units D, H, M or S) from the start of the run to start producing the cascade files

List of tracers to be written from piloting run

(LU only) Mesh length (resolution) in x-direction (degrees)

(LU only) Mesh length (resolution) in y-direction (degrees)

Reference Point I on rotated grid, 1 < Grd_iref < Grd_ni

Reference Point J on rotated grid, 1 < Grd_jref < Grd_nj

Latitude on rotated grid of reference point, Grd_iref,Grd_jref (degrees)

Longitude on rotated grid of reference point, Grd_iref,Grd_jref (degrees)

Max Supported Courrant number; Pilot area=Grd_maxcfl +Grd_bsc_base+Grd_bsc_ext1

Number of points along NI

Number of points along NJ

(GY only) Overlap extent along latitude axis for GY grid (degrees)

Type of grid described using 2 characters:

• "GY" : Global Yin-Yang
• "LU" : LAM Uniform

Geographic longitude of the center of the computational domain (degrees)

Geographic latitude of a point on the equator of the computational domain east of Grd_xlon1,Grd_xlat1 (degrees)

Geographic latitude of the center of the computational domain (degrees)

Geographic longitude of a point on the equator of the computational domain east of Grd_xlon1,Grd_xlat1 (degrees)

Save a restart file + continue at that time

Save a restart file + continue every Fcst_bkup_S

End date for model run slice (yyyymmdd.hhmmss)

Output global stat (glbstat) every Fcst_gstat_S

Read nesting data every Fcst_nesdt_S

Save a restart file + stop every Fcst_rstrt_S

Starting date for model run slice (yyyymmdd.hhmmss)

Setting for Fortran alarm time

Length of model timestep (sec)

Account for leap years

Starting date for model run (yyyymmdd.hhmmss)

Dcmip case selector

• Dcmip_case=0 (none)
• Dcmip_case=11 (3D deformational flow)
• Dcmip_case=12 (3D Hadley-like meridional circulation)
• Dcmip_case=13 (2D solid-body rotation of thin cloud-like tracer in the presence of orography)
• Dcmip_case=20 (Steady-state at rest in presence of oro.)
• Dcmip_case=21 (Mountain waves over a Schaer-type mountain)
• Dcmip_case=22 (As 21 but with wind shear)
• Dcmip_case=31 (Gravity wave along the equator)
• Dcmip_case=41 (Dry Baroclinic Instability Small Planet)
• Dcmip_case=43 (Moist Baroclinic Instability Simple physics)
• Dcmip_case=161 (Baroclinic wave with Toy Terminal Chemistry)
• Dcmip_case=162 (Tropical cyclone)
• Dcmip_case=163 (Supercell -Small Planet-)

Set lower value of Tracer in Terminator

• Dcmip_lower_value=0 (free)
• Dcmip_lower_value=1 (0)
• Dcmip_lower_value=2 (1.0e-15)

Account for moisture

• Dcmip_moist=0 (dry)
• Dcmip_moist=1 (moist)

Vertical Diffusion Theta (if <0,we remove REF)

Vertical Diffusion Tracers (if <0,we remove REF)

Vertical Diffusion Winds (if <0,we remove REF)

Type of planetary boundary layer

• Dcmip_pbl_type=-1 (none)
• Dcmip_pbl_type=0 (Reed-Jablonowski Boundary layer)
• Dcmip_pbl_type=1 (Georges Bryan Planetary Boundary Layer)

Type of precipitation/microphysics

• Dcmip_prec_type=-1 (none)
• Dcmip_prec_type=0 (Large-scale precipitation -Kessler-)
• Dcmip_prec_type=1 (Large-scale precipitation -Reed-Jablonowski-)

Do Rayleigh friction if T

Do Terminator chemistry if T

Earth's radius reduction factor

Rotation angle in DEGREE - W_NAIR=0)

Williamson case selector

• Williamson_case=0 (none)
• Williamson_case=1 (Advection 2D - Williamson_NAIR/Terminator_L)
• Williamson_case=2 (Steady-state nonlinear zonal geostrophic flow - Williamson et al.,1992,JCP,102,211-224)
• Williamson_case=5 (Zonal flow over an isolated mountain - Williamson et al.,1992,JCP,102,211-224)
• Williamson_case=6 (Rossby-Haurwitz wave - Williamson et al.,1992,JCP,102,211-224)
• Williamson_case=7 (The 21 December 1978 case - Williamson et al.,1992,JCP,102,211-224)
• Williamson_case=8 (Galewsky's barotropic wave - Galewsky et al.,2004,Tellus,56A,429-440)

LAT cosine Bell in DEGREE - W_NAIR=0)

LON cosine Bell in DEGREE - W_NAIR=0)

Set lower value of Tracer in Terminator

• Williamson_lower_value=0 (free)
• Williamson_lower_value=1 (0)
• Williamson_lower_value=2 (1.0e-15)

Used when Williamson_case=1

• Williamson_NAIR=0 (Solid body rotation of a cosine bell - Williamson et al.,1992,JCP,102,211-224)
• Williamson_NAIR=1 (Deformational Non-divergent winds - Lauritzen et al.,2012,GMD,5,887-901)
• Williamson_NAIR=2 (Deformational divergent winds - Lauritzen et al.,2012,GMD,5,887-901)
• Williamson_NAIR=3 (Deformational Flow for Circular vortex - Nair and Machenhauer,2002,MWR,130,649-667)

rotation period in SECS - W_NAIR=0)

LAT rot.POLE in DEGREE - W_NAIR=3)

LON rot.POLE in DEGREE - W_NAIR=0/3)

Do Terminator chemistry if T

Select closures for shallow convection

• 'CAPE'
• 'EQUILIBRIUM'

Select formulation of fractional detrainment rate for shallow convection

• 'BECHTOLD01'
• 'CUIJPERS95'
• 'DEROOY10'

Select formulation of fractional entrainment rate for shallow convection

• 'BECHTOLD01'
• 'BECHTOLD08'
• 'DEROOY11'
• 'SIEBESMA03'

Evaporate detrained condensate in shallow convection

Number of species for convective transport (never tested)

Number of additional ensemble members (max 3) for deep bkf convection

Take ice phase into account in deep bkf (yes=1)

Limit vertical computation by ktdia-1 levels

Activate convective transport of species for deep and shallow bkf

Allow downdrafts in deep bkf

Activate shallow convective momentum transport

Cloud radii at LCL for bkf_shallow from bkf_rads(1) to bkf_rads(2) with increment bkf_rads(3)

Temperature perturbations at LCL for triggering bkf_shallow An ensemble of shall. cumuli will be generated using perturbations starting from bkf_tperts(1) to bkf_tperts(2) with increment bkf_tperts(3)

Deep convection scheme name

• 'NIL ' :
• 'SEC ' :
• 'OLDKUO ' :
• 'KUOSTD ' :
• 'KFC ' :
• 'KFC2 ' :
• 'KFC3 ' :
• 'BECHTOLD' :

Treat convective clouds as cloud objects

Decay timescale for convective cloud objects (seconds)

Conservation corrections for deep convective scheme

• 'NIL ' : No conservation correction applied
• 'TEND ' : Temperature and moisture tendencies corrected
• 'PRECIP' : Surface precipitation rate corrected

Minimum depth of conv. updraft for KFC trigger (m)

Maximum depth of the downdraft detrainment layer (Pa) for 'kfc2' and 'kfc3'

generate wind tendencies in KFC or deep BKF if .true.

Compute production terms for Kain-Fritsch scheme

Initial convective updraft radius in KFC scheme(m)

Convective updraft radius over water in KFC scheme(m)

Varies convective timescale as a function of CAPE for Kain-Fritsch scheme KFCTAUCAPE = time1, time2, cmean, dcape

• time1 (s): max kfctimec
• time2 (s): min kfctimec
• cmean (J/Kg): cape value at which kfctimec will be mean of time1 and time2
• dcape (J/Kg): decrease in kfctimec from time1 to time2 will occur over range cmean-dcape to cmean+dcape

Trigger parameter of Kain-Fritsch convection scheme (WKLCL). Trigger parameter will increase from kfctrig4(3) to kfctrig4(4) [m/s] between timestep kfctrig4(1) and timestep kfctrig4(2)

Nominal resolution for which KFCTRIG4 is set. This is inactive if value <= 0.

Over land and lakes we keep the value set by the "ramp" above over sea water:

• for :lat: >= TRIGLAT(2) we keep value set by the "ramp" KFCTRIG4
• for :lat: <= TRIGLAT(1) we use the new value KFCTRIGL [m/s]
• and linear interpolation in between TRIGLAT(1) and TRIGLAT(2)

Logical key for variation of the trigger function depending on latitude and land-sea-lake mask

Relaxation timescale for trigger velocity

Trigger parameter of Kain-Fritsch convection scheme (WKLCL). Trigger parameter will increase from kfctrigw(3) to kfctrigw(4) [m/s] between wstar values kfctrigw(1) and kfctrigw(2)

Switch for mid-level convection

• 'NIL' : No mid-level convective scheme
• 'KF ' : Kain-Fritsch-based mid-level convective scheme

Conservation corrections for mid-level convective scheme

• 'NIL ' : No conservation correction applied
• 'TEND ' : Temperature and moisture tendencies corrected
• 'PRECIP' : Surface precipitation rate corrected

Minimum cloud depth for mid-level convection (m)

Downdraft detrainment depth for mid-level convection (Pa)

Fraction of environmental mass flux that enters updrafts

Modulation of the minimum environmental mass flux for mid-level convection

Maximum deep CAPE (J/kg/m2) for mid-level convective triggering

Minimum parcel departure level for mid-level convection (m)

Minimum environmental mass flux for mid-level convection (kg/s)

Switch for shallow convection

• 'NIL'
• 'KTRSNT'
• 'KTRSNT_MG'
• 'BECHTOLD'

Conservation corrections for shallow convective scheme

• 'NIL ' : No conservation correction applied
• 'TEND' : Temperature and moisture tendencies corrected

Over land and lakes we keep the value set by the "ramp" above over sea water:

• for :lat: >= TRIGLAT(2) we keep value set by the "ramp" KFCTRIG4
• for :lat: <= TRIGLAT(1) we use the new value KFCTRIGL
• and linear interpolation in between TRIGLAT(1) and TRIGLAT(2)

Time length (hours) for special time accumulated physics variables

Turbulent kinetic energy advect. is active if .true.

Surface heat flux from oceans is active if .true.

Clip tracers negative values

Conservation corrections for gridscale condensation

• 'NIL ' : No conservation correction applied
• 'TEND' : Temperature and moisture tendencies corrected

Evaporation parameter for Sunqvist gridscale condensation

Minimum cloud mixing ratio (kg/kg) for autoconversion in Sunqvist gridscale condensation

Max allowed values of modified hu00 (threshold relative humidity for stratiform condensation, Sunqvist gridscale condensation)

Min allowed values of modified hu00 (threshold relative humidity for stratiform condensation, Sunqvist gridscale condensation)

Standard deviation length scale (gridpoints) of Gaussian smoother applied to temperature and humidity inputs for Sunqvist gridscale condensation)

Activate Debug memory mode

Print a trace of the phy functions (MSG verbosity = debug)

Obtain estimate of surface wind gusts if .true.

Diffuse vertical motion if .true.

Surface friction is active if .true. Uses Schuman-Newell lapse rate if .false.

Minimal value for TKE in stable case (for 'CLEF')

Surface evaporation is active if .true.

Boundary layer processes

• 'NIL ': no vertical diffusion
• 'CLEF ': non-cloudy boundary layer formulation
• 'MOISTKE': cloudy boundary layer formulation
• 'SURFACE': TODO
• 'SIMPLE ': a very simple mixing scheme for neutral PBLs

(MOISTKE only) Apply factor fnn_reduc

• .false.: everywhere
• .true.: over water only

(MOISTKE only) Reduction factor (between 0. and 1.) to be applied to the parameter FNN (turbulent flux enhancement due to boundary layer clouds)

(CLEF+CONRES only) Non-dimensional parameter (must be >= 1.) that controls the value of the flux enhancement factor in CONRES

Use Fomichev radiation code if .true.

Gravity wave drag formulation

• 'NIL ': no Gravity wave drag
• 'GWD86': gravity wave drag + low-level blocking
• 'SGO16': new formulation (2016) of GWD86

Number of times the 3-point filter will be applied to smooth the GW flux profiles

Consider heating from non-orog. drag if = 1

Comma-separated list of diagnostic level inputs to read. Default: indiag_list_s(1) = 'DEFAULT LIST', expanded to: UU, VV, TT, HU + all dynamic Tracers

Initialize water content and cloud fraction seen by radiation for time 0 if .true.

Update ozone climatology during the run

Time between full radiation calculation (units D,H,M,S,P)

Compute ice fraction in KTRSNT_MG if .true.

Compute lightning diagnostics if .true. (currently for Milbrandt-Yau microphysics only)

Add methane oxydation as source of humidity in the stratosphere if .true.

Mixing length calc. scheme

• 'BLAC62 ': mixing length calc. using Blackadar
• 'BOUJO ': mixing length calc. using Bougeault
• 'TURBOUJO': mixing length calc. using Bougeault in turbulent regimes (otherwise Blackadar)
• 'LH ': mixing length calc. using Lenderink and Holtslag

Time length (hours) for special time averaged physics variables

Switch for aerosol activation scheme (1 = default, 2 = ARG + Aerosol climatology)

Switch for airmass type (1 = maritime, 2 = continental)

Double-moment for cloud (for 'mp_my' only)

Double-moment for graupel (for 'mp_my' only)

Double-moment for hail (for 'mp_my' only)

Double-moment for ice (for 'mp_my' only)

Double-moment for rain (for 'mp_my' only)

Double-moment for snow (for 'mp_my' only)

Compute MY Diagnostic fields if .true.

Ice-phase switched on if .true.

Initialize the number concentration for each category

Autoconversion (cloud to rain) switched on

Sedimentation switched on

Snow initiation switched on

Parameters for three-component freezing term

Warm-phase switched on

Physic input blocking along X

Physic input blocking along Y

Hines non-orographic GWD scheme is active if .true.

Pressure (in Pa) that defines the bottom emission level for gravity waves

Number of timesteps for which surface fluxes "FC" and "FV" are gradually set from 0 to their full value in a "slow start fashion" at the beginning of a time integration

switch for real-time debugging in microphysics (P3)

calibration factor for ice deposition in microphysics (P3)

Maximum time step (s) to be taken by the microphysics (P3) scheme, with time-splitting used to reduce step to below this value if necessary

Number of ice-phase hydrometeor categories to use in the P3 microphysics scheme (currently limited to <5)

precipitation fraction factor used by SCPF in microphysics (P3)

model resolution factor used by SCPF in microphysics (P3)

switch for subgrid cloud/precipitation fraction scheme (SCPF) in microphysics (P3)

calibration factor for ice sublimation in microphysics (P3)

Vectoc lenght physics memory space folding for openMP

Time-averaging of transfer coefficient for momentum to reduce 2-dt oscillations in fluxes

Conservation corrections for PBL scheme

• 'NIL ' : No conservation correction applied
• 'TEND' : Temperature and moisture tendencies corrected

Include the turbulent effects of trade wind cumulus clouds

Diffuse condensate fields

Run with a modified closure for the dissipation length scale

• 'NIL ' : No modified closure for the dissipation length scale
• 'LIM50' : A maximum value of 50m is imposed on dissipation length

Dissipative heating tendencies are computed for the PBL scheme such that total energy (kinetic + internal) is conserved

• 'NIL ' : No dissipative heating is computed
• 'LOCAL_K ' : Local total energy conservation based on diffusion coefficients
• 'LOCAL_TEND' : Local total energy conservation based on wind tendencies

Call surface scheme immediately before the PBL

Class of stability functions (stable case) to use in the PBL

• 'DELAGE97 ' : Use functions described by Delage (1997; BLM)
• 'LOUIS79 ' : Use functions described by Louis (1979; BLM)
• 'DERBY97 ' : Use functions described by Derbyshire (1997; Cardignton Tech Note)
• 'BELJAARS99' : Use functions described by Beljaars and Viterbo (1999; Clear and Cloudy Boundary Layers)
• 'BELJAARS91' : Use functions described by Beljaars and Holtslag (1991; JAM)
• 'LOCK07 ' : Use functions described by Lock (2007; Tech Report) employed at UKMO

Class of stability functions (unstable case) to use in the PBL

• 'DELAGE92' : Use functions described by Delage and Girard (1992; BLM)
• 'DYER74 ' : Use functions described by Dyer (1974; BLM)

Choose form of asymptotic mixing length for Blacadar-type estimates

• 'BLAC62' : Asymptotic 200 m proposed by Blackadar (1962; JGR) with clipping
• 'LOCK07' : Diagnosed asymptotic scale of Lock (2007; Tech Report) used at UKMO

Apply "turboujo" turbulence conditions to dissipation length scale

Run with legacy moistke clouds (no limits on cloud effects)

Use the non-local PBL cloud formulation

• 'NIL ' : no non-local PBL cloud formulation
• 'LOCK06' : Non-local cloud scheme of Lock and Mailhot (2006)

Use the mixing length to average the Richardson number profile of (potentially) many layers to derive a "background" Ri estimate

Richardson num. critical values for hysteresis

PBL representation of boundary layer clouds

• 'NIL ': No Shallow convection
• 'CONRES ': Bulk Richardson number-based turbulent enhancement
• 'SHALOW ': Deprecated (see 1998 RPN physics doc)
• 'SHALODQC': Deprecated (see 1998 RPN physics doc)
• 'GELEYN ': Deprecated (see 1998 RPN physics doc)

Layer over which to adjust from SL to PBL stability functions [(bot,top) in m]

Adjustment to coefficient for TKE diffusion

Control of time scale for TKE diffusion

Depth (Pa) of the always-turbulent near-surface layer in the PBL

Use true (motionless) surface boundary conditions for TKE diffusion

Relaxation timescale (s) for mixing length smoothing

Scheme to determine precipitation type

• 'NIL ': no call to bourge
• 'BOURGE ': use Bourgouin algorithm (bourge1) to determine precip. types.
• 'BOURGE3D':

Physic statistics output for 3d varables:

• .false. : mean, var, min and max for the whole 3d fiels
• .true.  : mean, var, min and max are done for each levels independently

Use double presision for physic statistics output

Physic statistics output Frequency

Print stats for phy_input read var

Physic statistics output: bus variable list that should be included in physics "block" stats. Possible values:

• Long varnames
• Short varnames
• 'ALLVARS=EDPV': all variables from E, D, P, V buses (any combination of the 4 letters);

CFC11 bckgrnd atmospheric concentration (PPMV)

CFC12 bckgrnd atmospheric concentration (PPMV)

CH4 bckgrnd atmospheric concentration (PPMV)

CO2 bckgrnd atmospheric concentration (PPMV)

N2O bckgrnd atmospheric concentration (PPMV)

Atmospheric path length for solar radiation

• 'RODGERS67' : Formulation used by Li and Barker (2005)
• 'LI06' : Estimate of Li and Shibata (2006)

Conservation corrections for radiation scheme

• 'NIL ' : No conservation correction applied
• 'TEND' : Temperature and moisture tendencies corrected

Use emissivity computed by the surface schemes

format of radiation files to be read

• 'STD': RPN standard file
• 'UNF': unformatted

Radiation fixes near the model top(for newrad only) if .true.

Vertical smoothing on radiative fluxes(for newrad only) if .true.

Use climatological values of GHG in radiation (CCCMARAD2 only)

Radiation scheme

• 'NIL ': no radiation scheme
• 'NEWRAD ': complete radiation scheme
• 'CCCMARAD ': most advanced radiation scheme
• 'CCCMARAD2': most advanced radiation scheme v2

List of levels on which IR and VIS radiation calculations are performed (to save on CPU time) (for newrad only)

Key for activation of the radiation along slopes

Additional output for low level refraction

Launching level value of GW RMS wind (m/s) from non-orographic origin

water/ice phase for saturation calc. if .true.; water phase only for saturation calc. if .false.

Tuning factor for blocking height

Sets the minimum value of the drag coefficient in the orographic blocking scheme.

Turns on/off the non-linear amplification factor (depending on wind direction) of the drag coefficient in the orographic blocking scheme

Critical phase for blocking height

Turns on/off the amplification factor (due to stability) of the drag coefficient in the orographic blocking scheme

Standard deviation length scale (gridpoints) of Gaussian smoother applied to wind GWD tendencies

Description of threshold for mean wind speed for blocking

Run ISCCP cloud simulator (cccmarad only) if .true.

Condensation scheme name

• 'NIL ' : No explicit condensation scheme used
• 'CONSUN ' : Sunqvist type condensation scheme
• 'NEWSUND ' : Sunqvist type condensation scheme
• 'MP_MY2_OLD' : Milbrandtl and Yau microphysics scheme (old formulation)
• 'MP_MY2 ' : Milbrandtl and Yau microphysics scheme
• 'MP_P3 ' : P3 microphysics scheme

Special treatment of stratosphere; if .true. ignore convection/condensation tendencies where pressure is lower than topc or specific humidity is lower than minq as specified in nocld.cdk

Factor used in the gwd formulation = 1/(LENGTH SCALE)

Run the physics in test harness mode

Print runtime timings

Select a turbulent orographic form drag scheme

• 'NIL'  : No turbulent orographic form drag scheme
• 'BELJAARS04' : Form drag scheme described by Beljaars et al. (2006; QJRMS)

(newrad only) Use TT(12000) instead of skin temp in downward IR flux calculation if .true.

Times series package stops at this timestep

List of time series for profile variables

List of time series for surface variables

Number of timesteps between time-series writeout

Stations chosen in lat,lon for time-series Format: "STN1_NAME",lat1,lon1, "STN2_NAME",lat2,lon2, ...

Adjust surface temperature over snow after reading (coherency check)

Prandtl number for neutral stability (initialized by SL module)

Diurnal SST scheme

• 'NIL ' : No Diurnal SST scheme
• 'FAIRALL' : #TODO: define

Diurnal SST scheme active coolskin if .true.

Diurnal SST scheme active coolskin over freshwater lakes if .true.

Diurnal SST scheme active warmlayer if .true.

Diurnal SST scheme active warmlayer over freshwater lakes if .true.

Depth of soil layers in [METERS] in SVS land surface scheme (schmsol=SVS)

Set water temperature of ice-covered lakes to 0C for points north of ice line if .true. needs an initialization file otherwise the model stops

Sea ice melting

Implicit surface fluxes if .true.; explicit fluxes if .false.

If .true. apply temporary fix to ISBA

• timestep dependent KCOEF
• No PSN factor for meting and freezing

Use the vegetation-only roughness length to compute vegetation snow fraction

If .true., freeze precipitation reaching the ground in sub-zero conditions

Deepest active (permeable) soil layer in SVS land surface scheme (schmsol=SVS)

Vegetation field update frequency (units D,H,M,S,P)

Minimum fraction of leads in sea ice.  Multiply ice fraction by (1.-leadfrac)

Limit snow depth to 10 cm for calculation of heat conductivity of snow over sea-ice and glacier if .true.

(coupling) fluxes over ocean are taken from ocean model if .true.

read-in land surface emissivity if .true.

Takes into account effect of ocean salinity on saturation specific humidity at ocean surface (boundary condition for LH flux calculation)

Land surface processes

• 'NIL ' : No Land surface processes
• 'ISBA' : Interaction Soil Biosphere Atmosphere (ISBA) land sfc scheme
• 'SVS ' : Soil, Vegetation, and Snow (SVS) (Multibudget) land sfc scheme

Urban surface processes

• 'NIL' : No Urban surface processes
• 'TEB' : Town Energy Balance (TEB) urban scheme

Class of stability functions (stable case) to use in the surface layer

• 'DELAGE97 ' : Use functions described by Delage (1997; BLM)
• 'BELJAARS91' : Use functions described by Beljaars and Holtslag (1991; JAM)
• 'LOCK07 ' : Use functions described by Lock (2007; Tech Report) employed at UKMO

Class of stability functions (unstable case) to use in the surface layer

• 'DELAGE92' : Use functions described by Delage and Girard (1992; BLM)
• 'DYER74 ' : Use functions described by Dyer (1974; BLM)

Mimimum Obukhov length (L) for soil surfaces

Define bulk Ri values for near-neutral regime in the surface layer

Use a reference roughness for surface layer calculations

Use snow albedo "I6" directly if .true.; Use snow age "XA" to calculate snow albedo if .false.

Soil texture database/calculations for SVS land surface scheme

• 'GSDE '  : 8 layers of sand & clay info from Global Soil Dataset for ESMs (GSDE)
• 'SLC '  : 5 layers of sand & clay info from Soil Landscape of Canada (SLC)
• 'SOILGRIDS' : 7 layers of sand & clay info from ISRIC — World Soil Information

Limit temperature inversions to 8K/40m in surface layer if .true.

OPTION TO USE PHOTOSYNTHESIS CODE FOR STOMATAL RESISTANCE in SVS

Factor multiplying stomatal resistance in ISBA

Use directional roughness length if .true.

Constant value of thermal roughness length (m) applied over water within latitudinal band defined by z0tlat

Minimum value of momentum roughness length (m)

Momentum roughness length formulation over water

• 'CHARNOCK' : #TODO: define
• 'BELJAARS' : #TODO: define

Thermal roughness length formulation over vegetation

• 'FIXED' : Uses z0h = z0m
• 'ZILI95': evolves with u*

Latitude (2 elements, in degrees) used to specify Z0T over water

• If :lat: <= Z0TLAT(1) constant Z0T.
• If :lat: >= Z0TLAT(2) Charnock's relation.
• In between, linear interpolation is used.

Thermal roughness length formulation over water

• 'MOMENTUM' : Uses z0h = z0m (replaces key z0trdps300=.false.)
• 'DEACU12'  : #TODO: define (replaces key z0trdps300=.true.)
• 'ECMWF'  : #TODO: define (New formulation used by ECMWF)

Height at which to compute screen-level temperature (m)

Height at which to compute anemomenter-level winds (m)