Single Column Model (SCM) output:

For the GABLS4, we will use the NetCdf format with the names, the units used in the DICE intercomparison case (http://appconv.metoffice.com/dice/dice.html), so the description below is almost just a copy of the DICE one. Nevertheless due to the specific focus on turbulent parameters and the observations available at the tower , we ask for additional variables (in red) and two NetCdf files :

1. Time series variables or fluxes (every time step)

2. Profiles of state variables, fluxes, forcing (10mn average or time step if greater than 10mn)

   The additional variables required for GABLS4 are in red

   For the profile, the data should be given on the model vertical grid

Filename

The filename should be composed as:

gabls4_profile_scm_<institute>_<model>_<exp>_<stageX>_<ver>.nc

gabls4_time_scm_<institute>_<model>_<exp>_<stageX>_<ver>.nc

where

Note: Spaces and underscores are not allowed in these strings.

example: gabls4_profile_MF_AROME_OPER_stage1_v01.nc and gabls4_time_MF_AROME_OPER_stage1_v01.nc

Global attributes

All relevant meta information should be added to the file as global attributes. This should include:

General:

reference to the model
contact person.
type of model where the SCM is derived from (climate model, mesoscale weather prediction model, regional model) ?
time step

Turbulence scheme:

Turbulence scheme  (e.g., K profile, TKE-l, EDMF ...)
Formulation of eddy diffusivity K.
Formulation of massflux, if used.
For E-l and Louis-type scheme: give formulation for length scale.
For K-profile: how is this  profile determined ? (e.g., based on Richardson, Brunt-Vaisala frequency (N^2),  Parcel method, other)

Surface/snow scheme:
Number of layer, prognostic variables such as temperature, snow albedo, snow density ...

and more:
Any other model specific aspects you find relevent for this intercomparison.
Any deviation from the prescribed set up that you had to make because of the specific structure of your model.

Initial profiles
Include the initial profiles in the mean state section as the first time step at 0 seconds.
Here we assume that the model is in eta coordinates and has a specific structure in the use of full levels and half level. If your model differs in this respect please adapt the output to fit with your model.

Sign convention
Surface energy fluxes (shf, lhf, g) are positive when directed away from the surface. Surface radiation fluxes (qdw, qup, ldw, lup) are all positive.

Vertical interpolation:
for the time series output at a given height between two full level please use linear interpolation and between surface and first model level please
specify the method in the general comment: (log, linear etc ...)

NetCdf dimensions and variables

Variables and dimensions should have the names as specified below between curled brackets (all lower case) . Exclude the curled brackets from the name.
For the fluxes like u'w', it is the total flux so in case of EDMF it is the diffusivity part + mass flux.

Each variable should have an attribute "units" with the unit prescribed as below between brackets. Exclude the brackets from the unit.

Each variable should have an attribute "long_name" which explains the meaning of the variable. The exact formulation is free, but could be taken from the description below. If a variable is not available for your model, use the attribute _FillValue to prescribe the numerical value that defines not available.

All physical variables should be of type float.

For the time series : gabls4_time_scm_<institute>_<model>_<stageX>_<ver>.nc
Dimensions:

{time} output times

Variables:

Time series output {time}

{time} in seconds since the beginning 00 UTC on December 11, 2009 [s]
{lwdw} long wave downward radiation at surface [W/m2]
{lwup} long wave upward radiation at surface [W/m2]
{swdw} short wave downward radiation at surface [W/m2]
{swup} short wave upward radiation at surface [W/m2]
{shf} sensible heat flux [W/m2]
{lhf} latent (Liq+sol) heat flux [W/m2]
{evap} evaporation+sublimation flux [mm/day]
{ustar} friction velocity [m/s]
{rain} precipitation (liq+sol) rate [mm/day]
{psurf} surface pressure [Pa]
{hpbl} boundary layer height [m]
{tsurf} surface temperature [K]
{trad} radiative temperature if different from tsurf [K]
{alb} surface albedo [0-1]
{z0m} momentum roughness length [m]
{z0h} heat roughness length [m]
{emis} surface emissivity [0-1]
{t2m} 2m temperature [K]
{q2m} 2m specific humidity [kg/kg]
{rh2m} 2m relative humidity [0-100] (computed with es(T) or ei(T))
{u10m} 10m u-component wind [m/s]
{v10m} 10m v-component wind [m/s]
{t3m} temperature at 3.30 meter above the surface [K]
{q3m} specific humidity at 3.30 meter above the surface [kg/kg]
{rh3m} relative humidity [0-100] at 3.30 meter above the surface (computed with the saturated vapor tension vs ice for T<0 or water if not)
{u3m} u-component wind at 3.30 meter above the surface [m/s]
{v3m} v-component wind at 3.30 meter above the surface [m/s]
{t9m} temperature at 8.80 meter above the surface [K]
{q9m} specific humidity at 8.80 meter above the surface [kg/kg]
{rh9m} relative humidity [0-100] at 8.80 meter above the surface (computed with the saturated vapor tension vs ice for T<0 or water if not)
{u9m} u-component wind at 8.80 meter above the surface [m/s]
{v9m} v-component wind at 8.80 meter above the surface [m/s]
{t18m} temperature at 17.90 meter above the surface [K]
{q18m} specific humidity at 17.90 meter above the surface [kg/kg]
{rh18m} relative humidity [0-100] at 17.90 meter above the surface (computed with the saturated vapor tension vs ice for T<0 or water if not)
{u18m} u-component wind at 17.90 meter above the surface [m/s]
{v18m} v-component wind at 17.90 meter above the surface [m/s]
{t25m} temperature at 25.30 meter above the surface [K]
{q25m} specific humidity at 25.30 meter above the surface [kg/kg]
{rh25m} relative humidity [0-100] at 25.30 meter above the surface (computed with the saturated vapor tension vs ice for T<0 or water if not)
{u25m} u-component wind at 25.30 meter above the surface [m/s]
{v25m} v-component wind at 25.30 meter above the surface [m/s]
{t33m} temperature at 32.70 meter above the surface [K]
{q33m} specific humidity at 32.70 meter above the surface [kg/kg]
{rh33m} relative humidity [0-100] at 32.70 meter above the surface (computed with the saturated vapor tension vs ice for T<0 or water if not)
{u33m} u-component wind at 32.70 meter above the surface [m/s]
{v33m} v-component wind at 32.70 meter above the surface [m/s]
{t42m} temperature at 41.90 meter above the surface [K]
{q42m} specific humidity at 41.90 meter above the surface [kg/kg]
{rh42m} relative humidity [0-100] at 41.90 meter above the surface (computed with the saturated vapor tension vs ice for T<0 or water if not)
{u42m} u-component wind at 41.90 meter above the surface [m/s]
{v42m} v-component wind at 41.90 meter above the surface [m/s]
{cc} total cloudcover fraction [0 1]
{alb} surface albedo [0 1]
{uw_3m} vertical flux u-component momentum at 3.30 meter above the surface[m2/s2]
{vw_3m} vertical flux v-component momentum at 3.30 meter above the surface[m2/s2]
{wt_3m} vertical temperature flux at 3.30 meter above the surface [Km/s]
{TKE_3m} turbulent kinetic energy at 3.30 meter above the surface [m^2/s^2]
{uw_7m} vertical flux u-component momentum at 7,03 meter above the surface[m2/s2]
{vw_7m} vertical flux v-component momentum at 7,03 meter above the surface[m2/s2]
{wt_7m} vertical temperature flux at 7,03 meter above the surface [Km/s]
{TKE_7m} turbulent kinetic energy at 7,03 meter above the surface [m^2/s^2]
{uw_15m} vertical flux u-component momentum at 15,43 meter above the surface[m2/s2]
{vw_15m} vertical flux v-component momentum at 15,43 meter above the surface[m2/s2]
{wt_15m} vertical temperature flux at 15,43 meter above the surface [Km/s]
{TKE_15m} turbulent kinetic energy at 15,43 meter above the surface [m^2/s^2]
{uw_23m} vertical flux u-component momentum at 22,79 meter above the surface[m2/s2]
{vw_23m} vertical flux v-component momentum at 22,79 meter above the surface[m2/s2]
{wt_23m} vertical temperature flux at 22,79 meter above the surface [Km/s]
{TKE_23m} turbulent kinetic energy at 22,79 meter above the surface [m^2/s^2]
{uw_30m} vertical flux u-component momentum at 30,15 meter above the surface[m2/s2]
{vw_30m} vertical flux v-component momentum at 30,15 meter above the surface[m2/s2]
{wt_30m} vertical temperature flux at 30,15 meter above the surface [Km/s]
{TKE_30m} turbulent kinetic energy at 30,15 meter above the surface [m^2/s^2]
{uw_38m} vertical flux u-component momentum at 37,51 meter above the surface[m2/s2]
{vw_38m} vertical flux v-component momentum at 37,51 meter above the surface[m2/s2]
{wt_38m} vertical temperature flux at 37,51 meter above the surface [Km/s]
{TKE_38m} turbulent kinetic energy at 37,51 meter above the surface [m^2/s^2]

For the profile : gabls4_profile_scm_<institute>_<model>_<stageX>_<ver>.nc

the profile at the time 300 is the average on the period ]0;600], at time 900 is ]600;1200] etc ...

Dimensions:

{time} output times
{levf} full levels
{levh} half levels
{levs} soil levels

4. Mean state {time} {levf}
{zf} height of full level [m]
{pf} pressure at full level [Pa]
{t} temperature [K]
{th} potential temperature [K]
{q} specific humidity [kg/kg]
{qc} cloud water and ice [kg/kg]
{u} zonal component wind [m/s]
{v} meridional component wind [m/s]

5. Prescribed forcings {time} ({levf} or {levh})
{ugeo} u-component geostrophic wind [m/s]
{vgeo} v-component geostrophic wind [m/s]
{dudt_ls} u-component momentum advection [m/s/s] (not used in GABLS4)
{dvdt_ls} v-component momentum advection [m/s/s] (not used in GABLS4)
{dtdt_ls} temperature advection [K/s]
{dqdt_ls} moisture advection [kg/kg/s]
{w} vertical movement [m/s] (not used in GABLS4)

6. Fluxes and increments {time} ({levf} or {levh})
{zh} height of half level [m]
{ph} pressure at half level [Pa]
{wt} vertical temperature flux [Km/s]
{wq} vertical moisture flux [kg/kg m/s]
{uw} vertical flux u-component momentum [m2/s2]
{vw} vertical flux v-component momentum [m2/s2]
{uu} u-variance [m2/s2]
{vv} v-variance [m2/s2]
{ww} w-variance [m2/s2]
{θ^2} potential temperature variance [K^2]
{Km} eddy diffusivity momentum [m2/s]
{Kh} eddy diffusivity heat [m2/s]
{mf} massflux [kg/m2/s]
{dT_dt_rad} temperature tendency from radiation [K/s]
{dT_dt_swrad} temperature tendency from short wave radiation [K/s]
{dT_dt_lwrad} temperature tendency from long wave radiation [K/s]
{TKE} turbulent kinetic energy [m^2/s^2]
{TTE} total turbulent energy [m^2/s^2]
{shear} shear production [m2/s3]
{buoy} buoyancy production [m2/s3]
{trans} total transport [m2/s3]
{dissi} dissipation [m2/s3]

7. Soil/snow variables {time} {levs}
{tsn} snow temperature [K]