namelists: GEM 4.8.12 physics

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.

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':
• 'PHYSIMP':

(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

• 'GWD86': gravity wave drag + low-level blocking
• 'GWD95': gravity wave drag (1995 version)

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

Time length (hours) for special time averaged physics variables

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

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

Vectoc lenght physics memory space folding for openMP

Run with a coupled high vertical resolution boundary layer (PBL).

Run with a modified closure for the dissipation length scale

• 'NIL'
• 'LIM50'

Limit the vertical extent of the PBL scheme to below the model level listed here

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

• 'NIL ':
• 'CONRES ':
• 'SHALOW ':
• 'SHALODQC':
• 'GELEYN ':

Time splitting (number of steps per model step) for the coupled PBL scheme

Relaxation timescale (s) for mixing length smoothing

Number of layers to split model layers into for the high vertical resolution coupled PBL

Scheme to determine precipitation type

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

Use double presision for physic statistics output

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

Physic statistics output Frequency

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)

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.

Radiation scheme

• 'NIL ': no radiation scheme
• 'OLDRAD ': simple 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.

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

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

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

Condensation scheme name

• 'NIL ' : No explicit condensation scheme used
• 'CONDS ' :
• 'CONSUN ' :
• 'NEWSUND ' :
• 'MP_MY2_OLD' :
• 'MP_MY2 ' :
• 'MP_P3 ' :

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

Use correct vertical levels (thermo) and offset (none) for Bougeault-Lacarrere mixing length calculation

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

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

Force re-calculation of deep bkf at every timestep

Set convective timescales for deep and shallow

Activate shallow convective momentum transport

Deep bkf timescale (s) if bkf_lsettadj= .false.

Shallow bkf timescale (s)if bkf_lsettadj= .false.

Deep convection scheme name

• 'NIL ' :
• 'SEC ' :
• 'OLDKUO ' :
• 'KFC ' :
• 'BECHTOLD' :

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

Total forced detrainment in KFC scheme

Init. level of forced detrainment in KFC scheme

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)

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

Time interval for refresh of tendencies in Kain-Fritsch scheme (s)

Convective time scale in Kain-Fritsch(s)

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

Switch for shallow convection

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

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)

Adjust surface temperature over snow after reading (coherency check)

Parameter for stability function

Diurnal SST scheme

• 'NIL ' :
• 'FAIRALL' :

Diurnal SST scheme active coolskin if .true.

Diurnal SST scheme active warmlayer if .true.

Uses dry adiabat if .true.

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. make sure there is soil water where MG > critmask (0.1%); WSoil = 0.3 if < critwater (0.001)

If .true. apply temporary fix to ISBA

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

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.

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

Land surface processes

• 'NIL ' :
• 'ISBA' :

Urban surface processes

• 'NIL' :
• 'TEB' :

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

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

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 roughness length (m) over water

Momentum roughness length formulation over water

• 'CHARNOCK' :
• 'BELJAARS' :

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 described by eq. 3 of

Height (m) of T and Q input for sfc fluxes calc.

Height (m) of wind input for sfc fluxes calc.

Format: "STN_NAME",lat,lon,
         ^         ^    ^ 
         string    2 reals 

type station_latlon
   character(STN_NAME_LENGTH) :: name
   real :: lat,lon
end type station_latlon