Parameters Reference Guide

This page defines the configurable parameters available in FastEddy®. Parameters are organized into logical sets by model feature or model configuration component.

Each table provides the:

  • Parameter name: Description

  • Minimum value (if applicable)

  • Maximum value (if applicable)

  • Default value (if applicable)

  • Req. (Requirement)

    • M (Mandatory)

    • O (Optional)

    • C-M (Conditionally Mandatory, with the condition listed in the description)

    • C-O (Conditionally Optional, with the condition listed in the description)

These definitions serve as a reference to ensure correct configuration and valid inputs for FastEddy simulations.

FEMPI

Parameter Name: Description

Min

Max

Default

Req.

numProcsX:
Number of cores to be used for horizontal domain decomposition in X

1

INT_MAX

1

M

numProcsY:
Number of cores to be used for horizontal domain decomposition in Y

1

INT_MAX

1

M

FECUDA

tBx:
Number of threads in x-dimension

1

INT_MAX

1

M

tBy:
Number of threads in y-dimension

1

INT_MAX

8

M

tBz:
Number of threads in z-dimension

1

INT_MAX

32

M

IO

inPath:
Path where initial/restart file is read in from

N/A

N/A

O

inFile:
Name of the input file for coordinate system and initial or restart conditions

N/A

N/A

O

outPath:
Path where output files are to be written

N/A

N/A

M

outFileBase:
Base name of the output file series as in (outFileBase).element-in-series

N/A

N/A

M

frqOutput:
Frequency (in timesteps) at which to produce output
Should be an even multiple of NtBatch

0

INT_MAX

M

ioOutputMode:
0: N-to-1 gather and write to a netcdf file
1: N-to-N writes of FastEddy binary files

0

1

0

O

GRID

Nx:
Number of discretised domain elements in the
x (zonal) direction

1

INT_MAX

128

M

Ny:
Number of discretised domain elements in the
y (meridional) direction

1

INT_MAX

122

M

Nz:
Number of discretised domain elements in the
z (vertical) direction

1

INT_MAX

122

M

Nh:
Number of halo cells to be used
(dependent on largest stencil extent)

0

INT_MAX

3

M

d_xi:
Computational domain fixed resolution in the ‘i’ direction

FLT_MIN

FLT_MAX

10.0

M

d_eta:
Computational domain fixed resolution in the ‘j’ direction

FLT_MIN

FLT_MAX

10.0

M

d_zeta:
Computational domain fixed resolution in the ‘k’ direction

FLT_MIN

FLT_MAX

10.0

M

coordHorizHalos:
Switch to setup coordinate halos as periodic=1 or gradient-following=0

0

1

1

M

topoFile:
A file containing topography
(surface elevation in meters ASL)

N/A

N/A

O

verticalDeformSwitch:
Switch to use vertical coordinate deformation 0=off; 1=on

0

1

M

verticalDeformFactor:
Deformation factor
0.0=max compression
1.0=no compression

0.0

1.0

M

verticalDeformQuadCoeff:
Quadratic term coefficient in the deformation scheme

-2.0

2.0

0.0

M

TIME INTEGRATION

timeMethod:
Selector for time integration method.
[0=RK3-WS2002 (default)]

0

0

0

O

Nt:
Number of timesteps to perform

1

INT_MAX

1000

M

dt:
Timestep resolution in seconds

FLT_MIN

FLT_MAX

1.0

M

NtBatch:
Number of timesteps to compute in batch launch
Must have NtBatch <= Nt
Should be equal to or a factor of frqOutput

1

Nt

1

M

HYDRO CORE

BOUNDARY CONDITIONS

hydroBCs:
Selector for hydro BC set.
1=Dirichlet lateral ceiling and surface boundary conditions (Limited Area Domain–LAD)
2=periodicHorizVerticalAbl

1

2

2

M

ceilingAdvectionBC:
Selector for normal advection at the domain ceiling
0=zero-valued w-BC
1=non-zero w-BC

0

1

0

O

hydroBndysFileBase:
Base file name LAD BC set
C-M if hydroBCs==1

N/A

N/A

C-M

hydroBndysFileStart:
Start counter value for LAD BC set files
C-M if hydroBCs==1

0

500000

0

C-M

hydroBndysFileEnd:
End counter value for LAD BC set files
C-M if hydroBCs==1

0

500000

0

C-M

dtBdyPlaneBCs:
Delta in time (seconds) between BdyPlane sets
C-M if hydroBCs==1

0.0

6e5

0.0

C-M

HYDRO_IO/LOGGING

hydroForcingWrite:
Switch for writing output of hydroFldsFrhs for prognostic fields
0=off; 1=on

0

1

0

O

hydroSubGridWrite:
Switch for writing output of SGS fields
0=off; 1=on

0

1

0

M

hydroForcingLog:
Switch for logging Frhs min/max etc
0=off; 1=on

0

1

0

M

ADVECTION

advectionSelector:
Advection scheme selector:
0=1st-order upwind
1=3rd-order QUICK
2=hybrid 3rd-4th order
3=hybrid 5th-6th order
4=3rd-order WENO
5=5th-order WENO
6=2nd-order centered

0

6

3

O

b_hyb:
Hybrid advection scheme parameter:
0.0=lower-order upwind
1.0=higher-order centered
0.0 < b_hyb < 1.0 = hybrid
C-M if advectionSelector==2 or 3

0.0

1.0

0.8

O; C-M

MOISTURE

moistureSelector:
Moisture selector:
0=off; 1=on

0

1

0

M

moistureNvars:
Number of moisture species
C-M if moistureSelector==1

0

2

0

C-M

moistureAdvSelectorQv:
Water vapor advection scheme selector
C-M if moistureSelector==1

0

6

0

C-M

moistureAdvSelectorQv_b:
Hybrid advection scheme parameter for water vapor
C-M if moistureSelector==1

0.0

1.0

0.0

C-M

moistureAdvSelectorQi:
Moisture advection scheme selector for non-qv fields (non-oscillatory schemes)
C-M if moistureSelector==1

0

2

0

C-M

moistureSGSturb:
Selector to apply sub-grid scale diffusion to moisture fields
C-M if moistureSelector==1

0

1

0

C-M

moistureCond:
Selector to apply condensation to moisture fields
C-M if moistureSelector==1

1

4

1

C-M

moistureCondTscale:
Relaxation time in seconds
C-M if moistureSelector==1

1e-4

1000.0

1.0

C-M

moistureCondBasePres:
Selector to use base pressure for microphysics
C-M if moistureSelector==1

0

1

0

C-M

moistureMPcallTscale:
Time scale for microphysics to be called
C-M if moistureSelector==1

1e-4

1000.0

1.0

C-M

CORIOLIS

coriolisSelector:
Coriolis Force selector:
0=None
1=Horizontal-only
2=Horizontal & Vertical

0

2

0

M

coriolisLatitude:
Characteristic latitude in degrees from equator of the LES domain

-90.0

90.0

54.0

M

TURBULENCE

turbulenceSelector:
Turbulence scheme selector:
0=None
1=Lilly/Smagorinsky

0

1

0

M

TKESelector:
Prognostic TKE selector:
0=None
1=Prognostic
2=requires canopySelector=1

0

2

1

O

TKEAdvSelector:
SGSTKE advection selector

0

6

0

O

TKEAdvSelector_b_hyb:
Hybrid advection scheme parameter

0.0

1.0

0.0

O

c_s:
Smagorinsky turbulence model constant used for turbulenceSelector=1 with TKESelector=0

1e-6

1e6

0.18

O

c_k:
Lilly turbulence model constant used for turbulenceSelector=1 with TKESelector>0

1e-6

1e6

0.10

O

CANOPY

canopySelector:
Canopy selector:
0=off; 1=on

0

1

0

O

canopySkinOpt:
Canopy selector to use additional skin friction effect on drag coefficient:
0=off; 1=on
C-M if canopySelector==1

0

1

0

C-M

canopy_cd:
Non-dimensional canopy drag coefficient
C-M if canopySelector==1

0.0

1e+2

0.15

C-M

canopy_lf:
Representative canopy element length scale
C-M if canopySelector==1

0.0

1e+2

0.1

C-M

DIFFUSION

diffusionSelector:
Diffusion term-type selector:
0=none
1=constant
2=scalar turbulent-diffusivity

0

1

0

O

nu_0:
Constant diffusivity used when diffusionSelector=1

0

FLT_MAX

1.0

O

AUXILIARY SCALARS AND SOURCES

NhydroAuxScalars:
Number of prognostic auxiliary scalar fields

0

256

0

O

AuxScAdvSelector:
Advection scheme for auxiliary scalar fields
C-M if NhydroAuxScalars > 0

0

6

0

C-M

AuxScAdvSelector_b_hyb:
Hybrid advection scheme parameter
C-M if NhydroAuxScalars > 0

0.0

1.0

0.0

C-M

AuxScSGSturb:
Selector to apply sub-grid scale diffusion to auxiliary scalar fields
C-M if NhydroAuxScalars > 0

0

1

0

C-M

srcAuxScFile:
The path and filename to an Auxilliary Scalar Sources specification file
C-O if NhydroAuxScalars > 0

N/A

N/A

C-O

srcAuxScTemporalType:
Temporal characterization of source (0=instantaneous; 1=continuous)
C-M if srcAuxScFile == NULL and
NhydroAuxScalars == 1

0

1

0

C-M

srcAuxScStartSeconds:
Source start time in seconds
C-M if srcAuxScFile == NULL and
NhydroAuxScalars == 1

0

FLT_MAX

0

C-M

srcAuxScDurationSeconds:
Source duration in seconds
C-M if srcAuxScFile == NULL and
NhydroAuxScalars == 1

0

FLT_MAX

30

C-M

srcAuxScGeometryType:
0=point (single cell volume)
1=line (line of surface cells)
C-M if srcAuxScFile == NULL and
NhydroAuxScalars == 1

0

0

0

C-M

srcAuxScMassSpecType:
Mass specification type
0=strict mass in kg
1=mass source rate in kg/s
C-M if srcAuxScFile == NULL and
NhydroAuxScalars == 1

0

0

0

C-M

srcAuxScMassSpecValue:
Mass specification value in kg or kg/s given by srcAuxScMassSpecType 0 or 1
C-M if srcAuxScFile == NULL and
NhydroAuxScalars == 1

0.0

FLT_MAX

1.0

C-M

srcAuxScLocation_X:
Source geometry centroid position in x (west-east)
C-M if srcAuxScFile == NULL and
NhydroAuxScalars == 1

-FLT_MAX

FLT_MAX

0.0

C-M

srcAuxScLocation_Y:
Source geometry centroid position in y (south-north)
C-M if srcAuxScFile == NULL and
NhydroAuxScalars == 1

-FLT_MAX

FLT_MAX

0.0

C-M

srcAuxScLocation_Z:
Source geometry centroid position in z (vertical above the surface)
C-M if srcAuxScFile == NULL and
NhydroAuxScalars == 1

-FLT_MAX

FLT_MAX

0.0

C-M

EXPLICIT FILTERS

filterSelector:
Explicit filter selector:
0=off; 1=on

0

1

0

O

filter_6thdiff_vert:
Vertical 6th-order filter on w selector:
0=off; 1=on
C-O if filterSelector == 1

0

1

0

C-O

filter_6thdiff_hori:
Horizontal 6th-order filter on rho; theta; qv selector:
0=off; 1=on
C-O if filterSelector == 1

0

1

0

C-O

filter_6thdiff_vert_coeff:
Vertical 6th-order filter factor:
0.0=off; 1.0=full
C-M if filterSelector == 1 and
filter_6thdiff_vert == 1

0.0

1.0

0.03

C-M

filter_6thdiff_hori_coeff:
Horizontal 6th-order filter factor:
0.0=off; 1.0=full
C-M if filterSelector == 1 and
filter_6thdiff_hori == 1

0.0

1.0

0.03

C-M

filter_divdamp:
Divergence damping selector:
0=off; 1=on
C-O if filterSelector == 1

0

1

0

C-O

RAYLEIGH DAMPING LAYER

dampingLayerSelector:
Rayleigh Damping Layer selector

0

1

0

O

dampingLayerDepth:
Rayleigh Damping Layer Depth
C-O if dampingLayerSelector==1

0.0

FLT_MAX

100.0

C-O

SURFACE LAYER

surflayerSelector:
Monin-Obukhov surface layer selector:
0=off
1=surface kinematic heat flux (surflayer_wth)
2=skin temperature rate (surflayer_tr)
3=skin temperature and moisture from boundary condition files (hydroBCs=1)

0

3

0

M

surflayer_z0:
Roughness length (momentum) when surflayerSelector > 0

1e-12

1e+0

0.1

M

surflayer_z0t:
Roughness length (temperature) when surflayerSelector > 0

1e-6

1e+1

0.1

M

surflayer_z0tdyn:
Dynamic z0t calculation following Zilitinkevich (1995) approach:
0=off
1=constant Zilitinkevich coeff
2=variable Zilitinkevich coeff

0

2

1

O

surflayer_wth:
Kinematic sensible heat flux at the surface
C-M if surflayerSelector==1

-5e+0

5e+0

0.0

C-M

surflayer_tr:
Temperature rate at the surface
(>0 for warming; <0 for cooling)
C-M if surflayerSelector==2

-1e+1

1e+1

0.0

C-M

surflayer_wq:
Sensible heat flux at the surface (kg/kg m s-1) when surflayerSelector = 1
C-M if moistureSelector==1

-5e+0

5e+0

0.0

C-M

surflayer_qr:
Water vapor rate (kg/kg h-1) when surflayerSelector = 2
C-M if moistureSelector==1

-1e+1

1e+1

0.0

C-M

surflayer_idealsine:
Selector for idealized sinusoidal surface heat flux or skin temperature forcing

0

1

0

O

surflayer_ideal_ts:
Start time in seconds for the idealized sinusoidal surface forcing
C-M if surflayer_idealsine > 0

0

1e+5

0.0

C-M

surflayer_ideal_te:
End time in seconds for the idealized sinusoidal surface forcing
C-M if surflayer_idealsine > 0

0

1e+5

0.0

C-M

surflayer_ideal_amp:
Maximum amplitude of the idealized sinusoidal surface forcing
C-M if surflayer_idealsine > 0

0

1e+3

0.1

C-M

surflayer_ideal_qts:
Start time in seconds for the idealized sinusoidal surface forcing (qv)
C-M if moistureSelector==1 and
surflayer_idealsine > 0

0

1e+5

0.0

C-M

surflayer_ideal_qte:
End time in seconds for the idealized sinusoidal surface forcing (qv)
C-M if moistureSelector==1 and
surflayer_idealsine > 0

0

1e+5

0.0

C-M

surflayer_ideal_qamp:
Maximum amplitude of the idealized sinusoidal surface forcing (qv)
C-M if moistureSelector==1 and
surflayer_idealsine > 0

0

1e+3

0.1

C-M

surflayer_qskin_input:
Selector for file input (restart) value for qskin under surflayerSelector == 2
C-O if moistureSelector==1

0

1

0

C-O

surflayer_offshore:
Offshore selector:
0=off;1=on

0

1

1

O

surflayer_offshore_visc:
Viscous term on z0m:
0=off; 1=on

0

1

1

O

surflayer_offshore_opt:
Offshore roughness parameterization:
==0 (Charnock)
==1 (Charnock with variable alpha)
==2 (Taylor & Yelland)
==3 (Donelan)
==4 (Drennan)
==5 (Porchetta)
C-M if surflayer_offshore > 0

0

5

4

O

surflayer_offshore_dyn:
Selector to use parameterized ocean parameters:
0=off; 1=on
C-O if surflayer_offshore > 0

0

1

1

C-O

surflayer_offshore_hs:
Significant wave height
C-M if surflayer_offshore > 0 and
surflayer_offshore_dyn == 0 and
surflayer_offshore_opt == 2 or 3 or 4 or 5

0

1e+2

0.0

C-M

surflayer_offshore_lp:
Peak wavelength
C-M if surflayer_offshore > 0 and
surflayer_offshore_dyn == 0 and
surflayer_offshore_opt == 2

0.1

1e+3

0.1

C-M

surflayer_offshore_cp:
Wave phase speed
C-M if surflayer_offshore > 0 and
surflayer_offshore_dyn == 0 and
surflayer_offshore_opt == 3 or 4 or 5

0.1

1e+2

0.1

C-M

surflayer_offshore_theta:
Wave/wind angle
C-M if surflayer_offshore > 0 and
surflayer_offshore_dyn == 0 and
surflayer_offshore_opt == 5 or
surflayer_offshore > 0 and
surflayer_offshore_dyn==1 and
surflayer_offshore_opt == 5

0.0

180.0

0.0

C-M

CELL PERTURBATION METHOD

cellpertSelector:
CP method selector:
0=off; 1=on

0

1

0

O

cellpert_nts:
Number of time steps after which perturbations are seeded

0

1e+6

500

O

cellpert_sw2b:
Only south and west bnds perturbed (for idealized cases):
0=off
1=south & west
2=south
C-O if cellpertSelector > 0

0

3

0

C-O

cellpert_amp:
Maximum amplitude for the potential temperature perturbations
C-O if cellpertSelector > 0

0.0

20.0

0.5

C-O

cellpert_gppc:
Number of grid points conforming the cell
C-O if cellpertSelector > 0

0

50

8

C-O

cellpert_ndbc:
Number of cells normal to domain lateral boundaries
C-O if cellpertSelector > 0

0

10

3

C-O

cellpert_kbottom:
Z-grid point where the perturbations start
C-O if cellpertSelector > 0

1

10

1

C-O

cellpert_ktop:
Z-grid point where the perturbations end
C-O if cellpertSelector > 0

0

200

20

C-O

cellpert_tvcp:
Time-varying CP method selector:
0=off; 1=on (hydroBCs==1)
C-O if cellpertSelector > 0

0

1

0

C-O

cellpert_eckert:
Eckert number for the potential temperature perturbations (hydroBCs==1)
C-O if cellpertSelector > 0

0.0

10.0

0.2

C-O

cellpert_tsfact:
Factor on the refreshing perturbation time scale (hydroBCs==1)
C-O if cellpertSelector > 0

0.0

10.0

1.0

C-O

BASE-STATE

stabilityScheme:
Scheme used to set hydrostatic stability-dependent Base-State EOS fields
0=none
1=profile
2=linear in theta

2

2

2

M

temp_grnd:
Air Temperature (K) at the ground used to set hydrostatic Base-State EOS fields

FLT_MIN

FLT_MAX

300.0

M

pres_grnd:
Pressure (Pa) at the ground used to set hydrostatic Base-State EOS fields

FLT_MIN

FLT_MAX

1.0e5

M

zStableBottom:
Height (m) of the first stable upper-layer when stabilityScheme = 1 or 2

0

FLT_MAX

1000.0

M

stableGradient:
Vertical gradient (K/m) of the first stable upper-layer when stabilityScheme = 1 or 2

FLT_MIN

FLT_MAX

0.1

M

zStableBottom2:
Height (m) of the second stable upper-layer when stabilityScheme = 2

0

FLT_MAX

1100.0

M

stableGradient2:
Vertical gradient (K/m) of the second stable upper-layer when stabilityScheme = 2

FLT_MIN

FLT_MAX

0.03

M

zStableBottom3:
Height (m) of the third stable upper-layer when stabilityScheme = 2

0

FLT_MAX

1500.0

M

stableGradient3:
Vertical gradient (K/m) of the third stable upper-layer when stabilityScheme = 2

FLT_MIN

FLT_MAX

0.03

M

thetaPerturbationSwitch:
Switch to include initial theta perturbations:
0=off; 1=on

0

1

0

M

thetaHeight:
Height below which to include initial theta perturbations:
(meters)

0.0

FLT_MAX

0.0

M

thetaAmplitude:
Maximum amplitude for theta perturbations:
thetaAmplitude*[-1;+1] K

0.0

2.0

0.0

M

U_g:
Zonal (West-East) component of the geostrophic wind (m/s)

-FLT_MAX

FLT_MAX

0.0

M

V_g:
Meridional (South-North) component of the geostrophic wind (m/s)

-FLT_MAX

FLT_MAX

0.0

M

z_Ug:
Height (m) above ground for linear geostrophic wind gradient (zonal component)

0.0

FLT_MAX

10000.0

M

z_Vg:
Height (m) above ground for linear geostrophic wind gradient (meridional component)

0.0

FLT_MAX

10000.0

M

Ug_grad:
Geostrophic wind gradient for U component
U_g gradient above z_Ug (ms-1/m)

-1e2

1e2

0.0

M

Vg_grad:
Geostrophic wind gradient for V component
V_g gradient above z_Vg (ms-1/m)

-1e2

1e2

0.0

M

physics_oneRKonly:
Selector to apply physics RHS forcing only at the latest RK stage:
0=off; 1=on

0

1

1

O

LARGE SCALE FORCINGS

lsfSelector:
Large-scale forcings selector:
0=off; 1=on

0

1

0

M

lsf_horMnSubTerms:
Large-scale subsidence terms switch:
0=off; 1=on
C-M if lsfSelector > 0

0

1

0

C-M

lsf_freq:
Large-scale forcing frequency (seconds)
C-M if lsfSelector > 0

1e-3

1e+3

1.0

C-M

lsf_w_surf:
Large-scale forcing to w at the surface
C-M if lsfSelector > 0

-1e+4

1e+4

0.0

C-M

lsf_w_lev1:
Large-scale forcing to w at the first specified level
C-M if lsfSelector > 0

-1e+4

1e+4

0.0

C-M

lsf_w_lev2:
Large-scale forcing to w at the second specified level
C-M if lsfSelector > 0

-1e+4

1e+4

0.0

C-M

lsf_w_zlev1:
Large-scale forcing to w height 1
C-M if lsfSelector > 0

0.0

1e+4

100.0

C-M

lsf_w_zlev2:
Large-scale forcing to w height 2
C-M if lsfSelector > 0

0.0

1e+4

200.0

C-M

lsf_th_surf:
Large-scale forcing to theta at the surface
C-M if lsfSelector > 0

-1e+4

1e+4

0.0

C-M

lsf_th_lev1:
Large-scale forcing to theta at the first specified level
C-M if lsfSelector > 0

-1e+4

1e+4

0.0

C-M

lsf_th_lev2:
Large-scale forcing to theta at the second specified level
C-M if lsfSelector > 0

-1e+4

1e+4

0.0

C-M

lsf_th_zlev1:
Large-scale forcing to theta height 1
C-M if lsfSelector > 0

0.0

1e+4

100.0

C-M

lsf_th_zlev2:
Large-scale forcing to theta height 2
C-M if lsfSelector > 0

0.0

1e+4

200.0

C-M

lsf_qv_surf:
Large-scale forcing to qv at the first specified level
C-M if lsfSelector > 0

-1e+4

1e+4

0.0

C-M

lsf_qv_lev1:
Large-scale forcing to qv at the first specified level
C-M if lsfSelector > 0

-1e+4

1e+4

0.0

C-M

lsf_qv_lev2:
Large-scale forcing to qv at the second specified level
C-M if lsfSelector > 0

-1e+4

1e+4

0.0

C-M

lsf_qv_zlev1:
Large-scale forcing to qv height 1
C-M if lsfSelector > 0

0.0

1e+4

100.0

C-M

lsf_qv_zlev2:
Large-scale forcing to qv height 2
C-M if lsfSelector > 0

0.0

1e+4

200.0

C-M

EXTENSIONS

GAD

GADSelector:
Generalized Actuator Disk Selector:
0=off; 1=on

0

1

0

O

turbineSpecsFile:
The path+filename to a turbine specifications file
C-O if GADSelector > 0

N/A

N/A

C-O

GADoutputForces:
Flag to include GAD forces in the output:
0=off; 1=on
C-O if GADSelector > 0

0

1

0

C-O

GADofflineForces:
Flag to compute GAD forces in an offline mode:
0=off; 1=on
C-O if GADSelector > 0

0

1

0

C-O

GADaxialInduction:
Flag to compute axial induction factor:
0=off (uses prescribed GADaxialIndVal)
1=on
C-O if GADSelector > 0

0

1

1

C-O

GADaxialIndVal:
Prescribed constant axial induction factor when GADaxialInduction==0
C-O if GADSelector > 0

0.0

1.0

0.02

C-O

GADrefSwitch:
Switch to use reference windspeed:
0=off; 1=on
C-O if GADSelector > 0

0

1

0

C-O

GADrefU:
Prescribed constant reference hub-height windspeed
C-M if GADSelector > 0 and
GADrefSwitch == 1

0.0

50.0

0.0

C-M

GADrefSampleWindow:
Sample duration (in seconds) over which to average per-timestep values (filtering out highest frequencies)
C-O if GADSelector > 0

1.0

300.0

30.0

C-O

GADrefSeriesLength:
Number of sampling windows over which to average again for reference velocity magnitude and direction
C-O if GADSelector > 0

1

360

10

C-O

GADForcingSwitch:
Switch to use the GADrefU-based or local windspeed in computing GAD forces:
0=local
1=ref
C-O if GADSelector > 0

0

1

0

C-O

URBAN

urbanSelector:
Urban selector:
0=off
1=on with no thermal contribution
2=on with thermal relaxation towards base state

0

2

0

O

cd_build:
c_d coefficient (m-1) used by the drag-based building formulation:
-c_d|u_i|u_i
C-O if urbanSelector > 0

0.0

1e+8

100.0

C-O

ct_build:
c_t coefficient (s-1) used by the drag-based building formulation:
-c_t(rho*theta-rho_b*theta_b) & -c_t(rho-rho_b)
C-O if urbanSelector > 0

0.0

1e+8

10.0

C-O

urban_heatRedis:
Selector to activate surface heat redistribution
C-O if urbanSelector > 0

0

1

0

C-O