namelists: GEM 4.8.12 dynamics

Type of grid described using 2 characters:

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

Number of points along NI

Number of points along NJ

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

I reference point (within Grd_ni)

J reference point (within Grd_nj)

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

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

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

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

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

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

Geographic longitude 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)

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

x horizontal resolution of target cascade grid (degrees)

y horizontal resolution of target cascade grid (degrees)

Number of points for the blending zone (Hblen_x)

TRUE to dump out permanent bus for cascade mode

I reference point (within Grdc_ni)

J reference point (within 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

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

List of tracers to be written from piloting run

Starting date for model run (yyyymmdd.hhmmss)

Starting date for model run slice (yyyymmdd.hhmmss)

End date for model run slice (yyyymmdd.hhmmss)

Read nesting data every Fcst_nesdt_S

Output global stat (glbstat) every Fcst_gstat_S

Save a restart file + stop every Fcst_rstrt_S

Save a restart file + continue every Fcst_bkup_S

Save a restart file + continue at that time

Stop model if timestep not done within Fcst_alarm_S

Account for leap years

Length of model timestep (sec)

switch: .true. : prints stats of traj and interp-var

switch: .true. : check if upstream points from north and south pe's for which an interpolation is requested are inside own advection source grid

switch: .true. : print diagnostics on the number of upstream positions exchanged among processors

switch: .true. : advection of tracers is monotonic

switch: .false.: positivity-preserving advection

Switch:.T.:MONO(CLIPPING) after Bermejo-Conde (Mass-conservation for Chemical Tracers)

Number of neighborhood zones in ILMC (Mass-conservation for Chemical Tracers)

Switch:.T.:MONO(CLIPPING) after ILMC (Mass-conservation for Chemical Tracers)

Activate conservation diagnostics > 0 (Mass-conservation for Chemical Tracers)

South 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)

West boundary in GY for an embedded LAM (Mass-conservation for Chemical Tracers)

East boundary in GY for an embedded LAM (Mass-conservation for Chemical Tracers)

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

Number of neighborhood zones in ILMC (MASS-CONSERVATION for Chemical Tracers)

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

Type of rebuild in SLICE

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

(MASS-CONSERVATION for Chemical Tracers)

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

South 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)

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

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

value for time off-centering of Semi-lag (hydrostatic part >=.5)

T* (basic state temperature) (K)

value for time off-centering parameter for nonhydrostatic part >=.5

Parameter controlling modified epsilon (Ver_epsi_8) [nonhydrostatic part]

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

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

array of model levels , 0.0 < HYB < 1.0

Logical control for latitudinal modulation of the horizontal diffusion coefficient

U and V 2 delta-x removal ratio when using HO_IMP range(0.0-1.0)

Theta 2 delta-x removal ratio when using HO_IMP range(0.0-1.0). If Hzd_lnr_theta <= 0 If Hzd_lnr_theta <= 0 no diffusion on theta

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

Hor. diffusion order on U and V when using HO_IMP del-n available: 2, 4, 6, 8

Hor. diffusion order on theta when using HO_IMP del-n available: 2, 4, 6, 8

Tracers Hor. diffusion order when using HO del-n (n=2,4,6,8)

Hor. diffusion type: HO_IMP, HO_EXP5P, HO_EXP9P or NIL

Filter cutoff period for Iau_weight_S='sin' in hours

The number of seconds between increment fields

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.

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 -> initialization is performed

true -> windowing is applied

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)

number of levels for top piloting

number of points for horizontal blending

number of levels for top blending

true for Zero physics tendency in blending area

true-> to force constant (fixed) boundary conditions

Types of horizontal interpolation from 3DF/BCS files onto model grid in GEMDM:

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

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

true to blend the topography with the pilot

.true. to print more information to std output

precistion in print glbstats

• 'LCL_4'
• 'LCL_8'

.true. to clip humidity variables on output

interval of output file name change

        case('P') !pas
           units_S = 'STEPS'
        case('S')
           units_S = 'SECONDS'
        case('M')
           units_S = 'MINUTES'
        case('H')
           units_S = 'HOURS'
        case('D')
           units_S = 'DAYS'

'etiket' used in output fields

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

number of layers close to the bottom of the model within which a linear interpolation of GZ

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

minimum of digits used to represent output units

list of levels for underground extrapolation and VT is done even if flcubic is .true.

multiplication factor of P_pbl_spng at latitude P_lmvd_high_lat

multiplication factor of P_pbl_spng at latitude P_lmvd_low_lat

latitude at which the multiplication factor becomes P_lmvd_weigh_high_lat

latitude at which the multiplication factor becomes P_lmvd_weigh_low_lat

Number of bits to perturb on initial conditions

Stride for perturbation on initial conditions

.true. -> NO-LOG

• .true. => cubic interpolation
• .false.=> linear interpolation

.true. => auto barotropic option

Use cubic interpolation for Coriolis winds

• true -> hydrostatic run
• false-> non-hydrostatic run

true to trigger 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
• 0: blend at init only
• > 0: blend at every nblendyy timestep

0: Not activated (DEFAULT)

1: Wet psadj

2: Dry psadj (source_ps required)

3: Dry psadj 

• =1 Last thermo level at momentum level
• =0 Last thermo level is staggered

true to have topography

.true. -> activating SETTLS scheme

.true. -> variable cappa in thermod. eqn.

• 0 -> NO advection
• 1 -> traditional advection
• 2 -> consistant advection with respect to off-centering

Use super winds for advection

true. => Modify slightly code behaviour to ensure bitpattern reproduction in restart mode using FST file

Apply water loading in the calculations

Use cubic interpolation in trajectory comp.

Use trapezoidal average for advection winds

Epsilon value for non Yin-Yang iterative solver

Epsilon value for the Yin-Yang iterative solver

• .true. to use FFT solver if possible
• .false. to use MXMA solver (slower,less precise)

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

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

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

type of solver:'ITERATIF','DIRECT'

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

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 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 profile transition shape('COS2' or 'LINEAR') Set the shape between Spn_start_lev and Spn_up_const_lev

Nudging relaxation timescale (eg. 10 hours )

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.

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

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'

list of tracers to be read from analyse

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

(non-LAM) Number of points along j on which to apply modulation to 0. of the diff. coefficient

number of levels from the top of the model (thickness of the vertical sponge layer)

.true. will subtract the zonal mean of U component before diffusion and then added back

Number of timesteps for topography to evolve

Timestep to start topography which varies in time

.true. to apply vertical sponge (mtn_cfg only?)

thickness (in meters) of vertical sponge (mtn_cfg only?)

f .T. and Zblen_L=.T., apply sponge also to temperature field (mtn_cfg only?)

true for divergence high level modulation in initial computation of Zdot

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

Switch to activate SKEB (3D MARKOV CHAINES)

switch to print global stat related to Markov chains, SKEB and PTP (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 diffusion (3D MARKOV CHAINES)

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

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

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

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

(3D MARKOV CHAINES)

low wave number truncation limit used in 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)

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)

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

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

value of stretch for 3D Markov chain (used by SKEB) (3D 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)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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