Data accompanying the paper "Stable Numerical Implementation of a Turbulence Scheme with Two Prognostic Turbulence Energies".

Data accompanying the paper "Stable Numerical Implementation of a Turbulence Scheme with Two Prognostic Turbulence Energies", submitted to Monthly Weather Review.

This repository contains results from the idealized simulations of GABLS1 case,
obtained by the OpenIFS single column model version with the TOUCANS two-energy
turbulence scheme.

Each simulation has its own folder, containing single NetCDF file 'progvar.nc'.
Experimental settings, explained in the paper, are following:

-----------------------------------------------------------------
folder timestep delta beta beta_tau
[s] [1] [1] [1]
-----------------------------------------------------------------
GABLS1_2TE-R_dt001s_delta00_bt15/ 1 0.00 1.0 1.5
GABLS1_2TE-R_dt045s_delta00_bt15/ 45 0.00 1.0 1.5
GABLS1_2TE-R_dt090s_delta00_bt15/ 90 0.00 1.0 1.5
-----------------------------------------------------------------
GABLS1_2TE-R_dt090s_delta25_bt15/ 90 0.25 1.0 1.5
GABLS1_2TE-R_dt090s_delta25_bt10/ 90 0.25 1.0 1.0
GABLS1_2TE-R_dt180s_delta25_bt10/ 180 0.25 1.0 1.0
-----------------------------------------------------------------

Names of the NetCDF fields examined in the paper are:

tke - turbulence kinetic energy [J/kg]
tte - turbulence total energy [J/kg]
efb3 - turbulent heat flux [W/m^2]
t - thermodynamic temperature [K]

They are all 2D fields, indexed by the timestep number and the model level.
Turbulent heat flux is given on the model half levels, remaining quantities
on the model full levels. Heights of the model levels are stored in the
NetCDF fields:

height_f - full level heights [m]
height_h - half level heights [m]

Simulation time in seconds is stored in the NetCDF field 'time'.

CAUTION:
Since the simulations were run with the prescribed forcings and only with
the turbulence scheme activated, some fields in NetCDF files are meaningless
(e.g. skin temperature).