Icon Database Main
Icon Database Main
The CC license “BY-NC-ND” allows others only to download the publication and share it with others
as long as they credit the publication, but they can’t change it in any way or use it commercially.
Publisher Editors
Deutscher Wetterdienst Daniel Reinert, FE13,
Business Area “Research and Development” Tel. +49 (69) 8062-2060, [email protected]
Frankfurter Straße 135 Helmut Frank, FE13,
63067 Offenbach Tel. +49 (69) 8062-2742, [email protected]
www.dwd.de Florian Prill, FE13,
Tel. +49 (69) 8062-2727, [email protected]
0.1.0 10.01.13 DR, FP Generated preliminary list of available GRIB2 output fields
0.2.0 12.07.13 DR, FP Added a short section describing the horizontal ICON grid.
AUMFL_S, AVMFL_S added to the list of available output fields
0.2.1 15.07.13 DR Provide newly available output fields in tabulated form. Change
levelType of 3D atmospheric fields from 105 (Hybrid) to 150 (
Generalized vertical height coordinate)
0.2.2 16.07.13 FP Short description of ICON’s vertical grid.
0.2.3 25.09.13 DR Added description of available First Guess and analysis fields
0.2.4 17.12.13 DR Added description of external paramater fields
0.3.0 24.01.14 DR Added information about horizontal output grids
0.3.1 24.01.14 DR Added information about newly available output field OMEGA
0.4.0 22.05.14 HF Added SKY-database documentation
0.4.1 15.07.14 DR Some documentation on statistical processing and minor up-
dates. New output fields ASWDIR_S, ASWDIFD_S, ASWDIFU_S,
DTKE_CON
0.4.2 10.09.14 DR New output fields CLCT_MOD, CLDEPTH
0.5.0 01.10.14 DR Description of IAU initialization method
0.5.1 15.10.14 DR Updated description of necessary input fields
0.5.2 31.10.14 DR Add full table with model half level heights
0.6.0 05.12.14 DR Add short introduction and fix some minor bugs
0.6.1 10.12.14 DR New output field APAB_S
0.7.0 16.12.14 DR Revised documentation of time invariant fields and a couple of
bug fixes
0.7.2 09.01.15 DR General GRIB2 description
0.8.0 15.01.15 FP, DR Couple of bug fixes regarding the available fields on triangular
and regular grids
0.8.1 16.01.15 FP, DR List of pressure-level variables available on triangular grids
0.8.2 16.01.15 FP List of height-level variables available on regular grids
0.8.3 16.01.15 DR List of variables exclusively available for V V = 0
0.8.4 06.02.15 FP, DR Details of internal interpolation onto lon-lat grids. Details re-
garding output frequency.
0.8.5 18.02.15 FP Additional pressure levels for regular grid output.
0.8.6 23.02.15 FP Formula for computing non-zero topography level height.
1.0.0 23.02.15 FP Additional table of model full levels.
1.0.1 24.02.15 DR Update on available forecast runs and time span.
iii
Revision History
1.0.2 27.02.15 FP Added tables for grid point with maximum topo height.
1.0.3 13.03.15 DR, FP Section on statistically processed fields.
1.1.0 15.04.15 FP, DR Section on ICON EU nest (preliminary).
1.1.1 07.07.15 HF Added SMA list, list of half levels for EU nest, modified output
lists to automatically write model level variables in the namelist
templates.
1.1.1 17.07.15 HF Preliminary add T_S because it is already written in operations.
Some other minor modifications.
1.1.2 14.08.15 FP Added note on ICON’s earth radius and a table summarizing
regular grids.
1.1.3 04.12.15 FP Added WW code table 6.1.
1.1.4 11.01.16 HF Updated examples how to retrieve ICON data from SKY.
1.1.5 22.01.16 AR Description of En-Var.
1.1.6 28.01.16 DR Extend tables by field specific lat-lon interpolation method.
1.1.7 11.04.16 DR, FP Add timeline of model changes.
1.1.8 06.07.16 HF Add DTKE_HSH and other minor corrections.
1.1.9 27.09.16 DR Update intro and timeline.
1.1.10 03.02.17 DR update lat-lon interpolation methods and timeline.
1.1.11 08.05.17 DR update version history.
1.1.12 13.07.17 DR Update description for output variable SOILTYP.
1.1.13 25.10.17 FP, DR Remove references to COSMO-EU
1.1.14 10.01.18 DR Bug fix regarding availability of CLCT_MOD on global domain
1.2.0 26.01.18 MD Documentation for ICON-EPS products added
1.2.1 31.01.18 FP Updated height tables (appendix)
1.2.2 12.03.18 DR Added new output fields EVAP_PL, and SMI. Further adaptions
to the list of available fields and updated timeline; info on down-
load of grids.
1.2.3 07.07.18 FP Added output field CAPE_ML, EVAP_PL, SMI for global domain,
native grid. Added output field ALB_SEAICE.
1.2.4 31.08.18 MD Updated table of probability products.
1.2.5 09.10.18 DR Updated history of model changes; updated description of out-
put products; added output fields SNOAG and HSNOW_MAX
1.2.6 01.11.18 HF Updated section on Sky database examples.
1.2.7 27.02.19 DR Updated history of model changes.
1.2.8 01.04.19 HF Updated list of pressure levels of EPS output.
1.2.9 30.04.19 HF Updated list of pressure levels of deterministic global ICON on
triangular grid; hourly output of 5 fields of ICON-EU on regular
grid until the end of the forecast.
1.2.10 08.08.19 DR Updated history of model changes and description for HZEROCL
1.2.11 30.09.19 HF, DR Updated output list of ICON-EPS. Updated description of T_G
and T_SO
1.2.12 20.04.20 DR Updated history of model changes and added new diagnostic
output field CEILING
1.2.13 10.06.20 DR Updated history of model changes
iv
Revision History
2.0.0 21.01.20 MB, FP, DR, MD, Common database reference for all models ICON global/-EU/-
CS, CG D2; partly restructured chapters about output variables
2.0.1 17.02.20 MD, CG, FP, DR, Ensemble model output is now contained in chapters 10-12
MB
2.1.0 15.05.20 MD, CM, MB, FP New tables for all EPS products
2.1.1 05.06.20 MB, DR New native grid for ICON-D2
2.1.2 19.01.21 DR, CG, MB Update a few product descriptions; hints about the rotated lat-
lon grid output
2.1.3 16.02.21 DR Add description for DBZ_850, DBZ_CMAX
2.1.4 30.06.21 FP,DR Updated forecast range for ICON-D2(-EPS)
2.1.5 09.09.21 HF New output field LPI_CON_CI_MAX for ICON-EU(-EPS)
2.1.6 20.10.21 HF Output of ICON-EU runs to 51 h
2.1.7 17.11.21 HF Output of LPI_CON_CI_MAX for ICON-EU also on regular grid
2.1.8 15.08.22 DR Fix units of AUMFL_S and AVMFL_S and remove non-existing DOI
2.2.0 23.11.22 HF Adapted for ICON with 120 levels. Output of LPI_CON_MAX
instead of LPI_CON_CI_MAX for ICON-EU on regular grid
2.2.1 12.01.23 HF Added LPI_CON_MAX to output of global ICON. Corrected sev-
eral time ranges, resolution of MODIS data.
2.2.2 03.04.23 HF Added ASOB_S_CS to output of global ICON and ICON-EU.
2.3.0 27.03.24 LS, JF, DR Added ICON-ART mixed ensemble and related variables plus
ICON-EU-NA2 domain.
2.3.1 23.04.24 HF Added VIS to output of ICON-EU and ICON-D2.
v
Simulations are believed by no one except those who
conducted them.
2. Grid geometry 3
2.1. Horizontal grid . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
2.2. Vertical grid . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
2.3. Refined subregion over Europe (“local nest”) . . . . . . . . . . . . . . . . . . . . . . . . . 6
2.3.1. ICON-EU . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
2.3.2. ICON-EU-NA2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
2.4. ICON-D2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
4. Analysis fields 17
4.1. Ensemble Data Assimilation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
4.2. Incremental analysis update . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
4.3. Initial state for the convection-resolving ICON-D2 . . . . . . . . . . . . . . . . . . . . . 22
4.3.1. Interpolated initial conditions and initialisation . . . . . . . . . . . . . . . . . . . 22
4.3.2. Data assimilation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
vii
Contents
11.Ensemble forecasts with ICON (global, nested, and limited area mode) 129
11.1. Initial Perturbations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 129
11.2. Ensemble Physics Perturbations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 130
11.3. Lateral boundary perturbations (limited area EPS) . . . . . . . . . . . . . . . . . . . . . 132
11.4. ICON Ensemble output fields in DWD databases . . . . . . . . . . . . . . . . . . . . . . 133
11.4.1. Model Output . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 133
11.4.2. Ensemble Products for the ICON-EPS (global) and ICON-EU-EPS . . . . . . . . 133
11.4.3. Ensemble Products for the ICON-D2-EPS . . . . . . . . . . . . . . . . . . . . . . 137
viii
Contents
Bibliography 159
Glossary 161
ix
1. Introduction
The ICOsahedral Nonhydrostatic model ICON (Zängl et al., 2015, Zängl et al., 2022) is the global and
regional numerical weather prediction model at DWD. It became operational at 2015-01-20, replacing
the former operational global model GME. In June 2015 a refined sub-region (nest) over Europe
was activated (ICON-EU), in order to replace the regional model COSMO-EU. On 2021-02-10 the
convection-permitting model setup ICON-D2 (i.e. using the limited-area mode of ICON) has replaced
COSMO-D2. Since that date, the entire NWP system at DWD is based on ICON.
The ICON modelling system as a whole is developed jointly by DWD, the Max-Planck Institute for
Meteorology (MPI-M), the German Climate Computing Center (DKRZ), MeteoSwiss, and the Karlsruhe
Institute for Technology (KIT). While ICON is the new working horse for short and medium range
weather forecast at DWD and MeteoSwiss, it embodies the core of a new climate modelling system at
MPI-M.
ICON analysis and forecast fields serve as initial and boundary data for a couple of different limited
area models: Since 2015-01-20, analysis and forecast fields of the deterministic forecast run at 13 km
horizontal resolution serve as initial and boundary data for
• RLMs (Relocatable Local Model) of the German armed forces,
• DWD’s wave models.
ICON-D2 (-EPS) is driven by the deterministic (ensemble) forecasts of the ICON-EU nest.
This document provides an overview of all ICON analysis and forecast fields that are stored in the
database SKY at DWD. A subset of these data is publicly available from DWD’s Open Data Server
https://opendata.dwd.de
For additional information on the Open Data Server, we refer to this webpage. The document at
hand also provides some selected information on ICON’s grid structure and the data assimilation sys-
tem. For more detailed information, in particular regarding ICON’s numerical algorithms and physical
parameterizations, the reader is referred to the ICON Tutorial (Prill et al., 2023).
http://www.dwd.de/DE/fachnutzer/forschung_lehre/numerische_wettervorhersage/nwv_
aenderungen/nwv_aenderungen_node.html
You can receive regular information about changes to the forecasting environment by subscribing to
this mailinglist.
If you encounter bugs or inconsistencies, or if you have suggestions for improving this document, please
contact one of the following colleagues:
1
Chapter 1. Introduction
2
2. Grid geometry
The horizontal ICON grid consists of a set of spherical triangles that seamlessly span the entire sphere.
The grid is constructed from an icosahedron (see Figure 2.1a) which is projected onto a sphere. The
spherical icosahedron (Figure 2.1b) consists of 20 equilateral spherical triangles. The edges of each
triangle are bisected into equal halves or more generally into n equal sections. Connecting the new edge
points by great circle arcs yields 4 or more generally n2 spherical triangles within the original triangle
(Figure 2.2a, 2.2b).
(a) (b)
Figure 2.1.: Icosahedron before (a) and after (b) projection onto a sphere
(a) (b)
Figure 2.2.: (a) Bisection of the original triangle edges (b) More general division into n equal sections
ICON grids are constructed by an initial root division into n sections (Rn) followed by k bisection steps
(Bk), resulting in a RnBk grid. Figures 2.3a and 2.3b show R2B00 and R2B02 ICON grids. Such
grids avoid polar singularities of latitude-longitude grids (Figure 2.3c) and allow a high uniformity in
resolution over the whole sphere.
Throughout this document, the grid is referred to as the “RnBk grid” or “RnBk resolution”. For a
3
Chapter 2. Grid geometry
Figure 2.3.: (a) R2B00 grid. (b) R2B02 grid. (c) traditional regular latitude-longitude grid with polar
singularities
given resolution RnBk, the total number of cells, edges, and vertices can be computed from
nc = 20 n2 4k
ne = 30 n2 4k
nv = 10 n2 4k + 2
4π re2
∆A = ,
nc
with the earth radius re , and nc the total number of cells. ICON uses an earth radius of
re = 6.371229 · 106 m.
Based on ∆A one can derive an estimate of the average grid resolution ∆x:
5050 · 103
r
p π re
∆x = ∆A = k
≈ [m]
5 n2 n 2k
Visually speaking, ∆x is the edge length of a square which has the same area as our triangular cell.
In Table 2.1, some characteristics of frequently used global ICON grids are given. The table contains
information about the total number of triangles (nc ), the average resolution ∆x, and the maximum/min-
imum cell area. The latter may be interpreted as the area for which the prognosed meteorological quan-
tities (like temperature, pressure, . . . ) are representative. Some additional information about ICON’s
horizontal grid can be found in Wan et al. (2013).
4
2.2. Vertical grid
Table 2.1.: Characteristics of frequently used global ICON grids. ∆Amax and ∆Amin refer to the maximum
and minimum area of the grid cells, respectively.
Grid number of cells avg. resolution [km] ∆Amax [km2 ] ∆Amin [km2 ]
(nc )
The operational deterministic version of ICON is based on the R3B07 grid (∆x ≈
13 km), while the ensemble version (ICON-EPS) is based on the R3B06 grid (∆x ≈
26 km) Until 2022-11-22 the ensemble version was based on the R2B06 grid (∆x ≈ 40 km)
The vertical grid consists of a set of vertical layers with height-based vertical coordinates. Each of these
layers carries the horizontal 2D grid structure, thus forming the 3D structure of the grid. Close to
the surface the vertical layers are terrain following, while with increasing distance from the surface the
terrain signal is smoothed out and they gradually evolve into layers of constant height. The ICON grid
employs a Lorenz-type staggering with the vertical velocity defined at the boundaries of layers (half
levels) and the other prognostic variables in the center of the layer (full levels) (see Figure 2.4).
To improve simulations of flow past complex topography, the ICON model employs a smooth level
vertical (SLEVE) coordinate (Leuenberger et al., 2010). It allows for a faster transition to smooth
levels in the upper troposphere and lower stratosphere, as compared to the classical height-based Gal-
Chen coordinate. In the operational setup (deterministic and ensemble), the transition from terrain
following levels in the lower atmosphere to constant height levels is completed at z = 16 km. Model
levels above are flat. The required smooth large-scale contribution of the model topography is generated
by digital filtering with a ∇2 -diffusion operator.
Figure 2.5 shows the (half) level heights and layer thicknesses of the operational ICON setup with 120
vertical levels. This figure applies to the deterministic and the ensemble system, as both share the same
vertical grid.The table to the right shows the height above ground of selected half levels (for zero height
topography) and the corresponding pressure, assuming the US standard atmosphere. Standard heights
for all 121 half levels are given in Table A.1.
Please note that for grid cells with non-zero topography these values only represent rough
estimates of the true level height. Actual heights and layer thicknesses may vary consid-
erably from location to location, due to grid level stretching/compression over non-zero
topography (see e.g. the layer compression which is visible in Figure 2.4).
5
Chapter 2. Grid geometry
half level 1
full level 1
half level 2
w
full level 2
vn , ρ, θv
half level 3
full level 3 w
half level 4
Figure 2.4.: Illustration of ICON’s vertical levels. With num_lev layers, there are num_lev + 1 so-called
half levels. The half levels k − 1/2, k + 1/2 enclose layer k at the centers of which are the
corresponding full levels k, for k = 1, . . . , num_lev. Layer 1 is at the top of the atmosphere
and layer n at the bottom of the atmosphere. Half level num_lev + 1 coincides with the Earth’s
surface.
2.3.1. ICON-EU
The operational ICON has one refined subregion over Europe (ICON-EU). Key figures like edge coor-
dinates and mesh size of the ICON-EU nest are given in Table 2.2. The geographical location of the
6
2.4. ICON-D2
level [m] [P a]
70,000 1 75 000 2.1
1,400 5 68 806 5.6
10 61 485 16.5
60,000
15 54 656 41.8
z [m] ∆z [m] 1,200
20 48 323 93.8
25 42 491 198.7
30 37 159 411.1
50,000
1,000
35 32 328 826.6
40 27 998 1 586.7
45 24 170 2 854.4
40,000 50 20 795 4 830.8
800
55 17 908 7 614.0
60 15 637 10 893.4
30,000 65 13 933 14 251.2
600
70 12 433 18 054.1
75 10 933 22 871.6
20,000 80 9 433 28 812.4
400
85 7 933 35 945.4
90 6 433 44 445.2
95 4 933 54 504.1
10,000 200
100 3 449 66 200.5
105 2 164 77 893.2
110 1 142 88 332.9
0 0
20 40 60 80 100 120 115 409 96 507.3
level 120 20 101 085.0
Figure 2.5.: Vertical half levels (blue) and layer thickness (red) of the ICON operational setup (deterministic
and ensemble). The table of selected pressure values (for zero height) is based on the 1976 US
standard atmosphere.
Model simulations including the nested region over Europe are running regularly,
starting from
Main forecasts starting at 00, 06, 12, 18 UTC reach out to 120 h, while additional short-
range forecasts starting at 03, 09, 15, 21 UTC provide data until +51 h.
2.3.2. ICON-EU-NA2
The deterministic forecast and the ensemble members that include prognostic mineral dust have a
refined subregion over Europe, Northern Africa and the North Atlantic (ICON-EU-NA2 ) instead of
the ICON-EU domain. This refined subregion is considerably larger than the ICON-EU domain. Key
figures like edge coordinates and mesh size of the ICON-EU-NA2 nest are given in Table 2.4. The
geographical location of the nest is visualized in Fig. 2.7 (bottom).
Model simulations including prognostic mineral dust are running regularly, starting
from
Main forecasts starting at 00 and 12 UTC reach out to 180 h on the global and to 120 h
in the ICON-EU-NA2 domain. Main forecasts starting at 06 and 18 UTC provide data
until +120 h, and the additional short-range forecasts starting at 03, 09, 15, 21 UTC until
+51 h.
Chapter 2. Grid geometry
Figure 2.6.: ICON grid refinement (zoom view). Blue and red dots indicate the cell circumcenters for the
global (“parent”) and the refined (“child”) domain, respectively.
Table 2.4.: Key figures of the ICON-NA2 and ICON-NA2 -EPS nest.
8
2.4. ICON-D2
Figure 2.7.: Horizontal extent of the R3B8/R3B7 (∆x ≈ 6.5/13 km) ICON-EU/ICON-EU-EPS nest
(top; greenish blue area) and the R3B7 (∆x ≈ 13 km) ICON-EU-NA2 /ICON-EU-NA2 -EPS
nest (bottom) in a cylindrical equidistant projection.
9
Chapter 2. Grid geometry
Figure 2.8.: ICON-D2 domain in cylindrical equidistant projection. It is comprised of a regional R19B07
grid with 542040 cells and a horizontal resolution of ∆x ≈ 2 km. The highlighted frame depicts
the lateral boundary zone where the model is forced by externally specified data sets. It has a
total width of 14 cell rows and consists of two sub-zones named interpolation zone (rows 1-4)
and nudging zone (rows 5-14). The prognostic region starts at cell row 5. In the nudging zone
the prognostic solution is nudged towards the data of the driving model.
22,000
2,400
20,000
2,200
18,000
2,000
z [m] ∆z [m]
16,000
1,800
14,000
1,600
12,000 1,400
10,000 1,200
1,000
8,000
800
6,000
600
4,000
400
2,000
200
0 0
10 20 30 40 50 60
level
Figure 2.9.: Vertical half levels (blue) and layer thickness (red) of the ICON-D2 operational setup (deter-
ministic and ensemble).
10
3. Mandatory input fields
Several input files are needed to perform runs of the ICON model. These can be divided into three
classes: Grid files, external parameters, and initialization (analysis) files. The latter will be described
in Chapter 4.
11
Chapter 3. Mandatory input fields
Figure 3.1.: Screenshots of the ICON download server hosted by the Max Planck Institute for Meteorology
in Hamburg.
For fixed domain sizes and resolutions a list of grid files has been pre-built for the ICON model together
with the corresponding reduced radiation grids and the external parameters.
The contents of the primary storage directory are regularly mirrored to a public web site for download,
see Figure 3.1 for a screenshot of the ICON grid file server. The download server can be accessed via
http://icon-downloads.mpimet.mpg.de
The pre-defined grids are identified by a centre number, a subcentre number and a numberOfGridUsed,
the latter being simply an integer number, increased by one with every new grid that is registered in the
download list. Also contained in the download list is a tree-like illustration which provides information
on parent-child relationships between global and local grids, and global and radiation grids, respectively.
Note that the grid information of some of the older grids (no. 23 – 40) is split over two files: The users
need to download the main grid file itself and a grid connectivity file (suffix -grfinfo.nc).
12
3.2. External Parameters
set of raw datafields which are listed in Table 3.1. For a more detailed overview of ExtPar, the reader
is referred to the User and Implementation Guide of Extpar.
Table 3.1.: Raw datasets from which the ICON external parameter fields are derived.
MERIT and REMA are used for ICON, and ASTER is used for ICON-D2. GlobCover 2009 is a land
cover database covering the whole globe, except for Antarctica. Therefore, we make use of GlobCover
2009 for latitudes 90◦ > ϕ > −56◦ and switch to the coarser, however globally available dataset GLCC
for −56◦ ≥ ϕ > −90◦ .
The products generated by the ExtPar software package are listed in Table 3.2 together with the
underlying raw dataset. These are mandatory input fields for assimilation- and forecast runs.
13
Chapter 3. Mandatory input fields
Table 3.2.: External parameter fields for ICON, produced by the ExtPar software package (in alphabetical
order)
14
3.2. External Parameters
Note that fields marked with (*) are not required in operational model runs. I.e. the surface roughness
Z0 is only required, if the additional contribution from sub-grid scale orography shall be taken into
account (i.e. for itype_z0=1). In operational runs this is not the case. Instead, land-cover class specific
roughness lengths are taken from a GlobCover-based lookup table. FOR_D, FOR_E, LAI_MX, PLCOV_MX,
RSMIN, and ROOTDP became obsolete with the activation of the surface tile approach (2015-03-04). The
latter 4 fields are replaced by land-cover class specific values taken from lookup tables.
Remarks on post-processing
Some of the external parameter fields are further modified by ICON. The following fields are affected:
Hence, for post-processing tasks the modified external parameter fields should be used rather than the
original fields, for consistency. See Section 5.1.1 for more details.
15
4. Analysis fields
Numerical weather prediction (NWP) is an initial value problem. The ability to make a skillful forecast
relies heavily on an accurate estimate of the present atmospheric state, known as the analysis. In
general, an analysis is generated by optimally combining all available observations with a short-range
model forecast, known as first guess (FG) or background. Currently an atmospheric analysis is created
every 3 h. The 3-hourly first guess output provided by ICON comprises the following fields:
Table 4.1.: Available 3h first guess output fields from the forecast database
CAT_NAME=$model_ass_fc_$suite
Atmospheric analysis fields are computed every 3 hours (00, 03, 06,. . . 21 UTC) by the 3DVar data
assimilation system, which has recently been upgraded to an En-Var system (see Section 4.1). Sea
surface temperature T_SEA and sea ice cover FR_ICE are provided once per day (00 UTC) by the SST-
Analysis. A snow analysis is conducted every 3 hours, providing updated information on the snow height
H_SNOW, and fresh snow factor FRESHSNW. In addition a soil moisture analysis (SMA) is conducted once
per day (00 UTC). It basically modifies the soil moisture content W_SO, in order to improve the 2 m
temperature forecast.
For the 3-hourly assimilation cycle and forecast runs, ICON must be provided with 2 input files: One
containing the First Guess (FG) and the other containing analysis (AN) fields, only. Variables for which
no analysis is available are always read from the first guess file (e.g. TKE). Other variables may be read
either from the first guess or the analysis file, depending on the starting time. E.g. for T_SEA the first
guess is read at 03, 06, 09, 12, 15, 18, 21 UTC, however, the analysis is read at 00 UTC when a new
SST analysis is available. In Table 4.2 the available and employed first guess and analysis fields are
listed as a function of starting time.
17
Chapter 4. Analysis fields
Table 4.2.: The leftmost column shows variables that are mandatory for the assimilation cycle and forecast
runs. Column 2 indicates, whether or not an analysis is performed for these variables. Columns
3 to 10 show the origin of these variables (analysis or first guess), depending on the starting
time.
ShortName Analysis 00 03 06 09 12 15 18 21
Atmosphere
VN – FG FG FG FG FG FG FG FG
THETA_V – FG FG FG FG FG FG FG FG
DEN – FG FG FG FG FG FG FG FG
W – FG FG FG FG FG FG FG FG
TKE – FG FG FG FG FG FG FG FG
QC, QI, QR, QS – FG FG FG FG FG FG FG FG
QV 3DVar AN AN AN AN AN AN AN AN
T 3DVar AN AN AN AN AN AN AN AN
P 3DVar AN AN AN AN AN AN AN AN
U, V 3DVar AN AN AN AN AN AN AN AN
Surface
Z0 – FG FG FG FG FG FG FG FG
T_G – FG FG FG FG FG FG FG FG
QV_S – FG FG FG FG FG FG FG FG
T_SO(0:nlevsoil) – FG FG FG FG FG FG FG FG
W_SO_ICE – FG FG FG FG FG FG FG FG
W_SO SMA AN FG FG FG FG FG FG FG
W_I – FG FG FG FG FG FG FG FG
T_SNOW – FG FG FG FG FG FG FG FG
RHO_SNOW1 Ana_SNOW AN AN AN AN AN AN AN AN
H_SNOW Ana_SNOW AN AN AN AN AN AN AN AN
FRESHSNW Ana_SNOW AN AN AN AN AN AN AN AN
SNOWC – FG FG FG FG FG FG FG FG
HSNOW_MAX – FG FG FG FG FG FG FG FG
W_I – FG FG FG FG FG FG FG FG
EVAP_PL – FG FG FG FG FG FG FG FG
Sea ice/Lake
T_ICE – FG FG FG FG FG FG FG FG
H_ICE – FG FG FG FG FG FG FG FG
Continued on next page
18
4.1. Ensemble Data Assimilation
FR_ICE Ana_SST AN FG FG FG FG FG FG FG
ALB_SEAICE – FG FG FG FG FG FG FG FG
T_SEA Ana_SST AN FG FG FG FG FG FG FG
T_MNW_LK – FG FG FG FG FG FG FG FG
T_WML_LK – FG FG FG FG FG FG FG FG
H_ML_LK – FG FG FG FG FG FG FG FG
T_BOT_LK – FG FG FG FG FG FG FG FG
C_T_LK – FG FG FG FG FG FG FG FG
1 Note that RHO_SNOW is read from the analysis, however it does not contain any new/independent information compared
to the model first guess, except for an initialization of newly generated snow points and a limitation over glacier points.
W_SNOW is re-diagnosed within the ICON-code based on the analysed snow height H_SNOW and the former mentioned
snow density RHO_SNOW.
19
Chapter 4. Analysis fields
analysis is run once a day in the deterministic forecast system. By adjusting only the ensemble
mean the ensemble spread is preserved.
Snow
For each ensemble member the mean ensemble snow cover is adjusted to its deterministic value.
The data assimilation system also provides a couple of fields, which are not modified with respect to
their guess values, so that a full set of nominal analysis fields is available.
Table 4.3.: Fields provided by the ensemble analysis system. The column Increment indicates if an analysis
increment is provided. Analysis indicates if the field is analysed by the LETKF (letkf), taken from
the first guess (fg), interpolated (det) from, or (mean) adjusted to the respective deterministic
quantity, or additionally perturbed (per).
20
4.2. Incremental analysis update
et al., 2004) is applied. The basic idea is quite simple: Rather than adding the analysis increments
∆xA = xA − xF G ( i.e. the difference between the analysis xA and the model first guess xF G ) in one
go, they are incorporated into the model in small drips over many timesteps (see Figure 4.1).
Figure 4.1.: Incremental Analysis Update. Analysis increments are added to the background state (FG) in
small drips over some time interval rather than in one go. Currently, increments for U, V, P, T,
QV are treated in this way.
Mathematically speaking, during forward integration the model is forced with appropriately weighted
analysis increments:
Z
dx
= Ax + g(t)∆x A
, with g(t) dt = 1 (4.1)
dt
x is the discrete model state, A is a matrix representing the (non)-linear dynamics of the system and
g(t) is a weighting function, which is non-zero over some time-interval ∆t.
This drip by drip incorporation acts as a low pass filter in frequency domain on the analysis increments
such that small scale unbalanced modes are effectively filtered (see Bloom et al. (1996)). The filter
characteristic depends on the weighting function g(t). It should be noted that IAU only filters the
increments and not the background state, such that regions where analysis increments are zero remain
unaffected. This method is currently applied to the prognostic atmospheric fields π, ρ, vn , qv , based on
analysis increments provided for u, v, p, t and qv . π denotes the Exner pressure.
The method sounds incredibly simple, however there are a few technical aspects to be taken care of when
implementing this into an operational system: Figure 4.2 shows how the IAU-method is implemented
in ICON for a 3 h assimilation run starting at midnight. Analysis increments are applied over a 3 h time
window, centered at the actual model start time. As indicated by the blue line, constant weights are
used:
∆t
g(t) = , for − T /2 < t < T /2 (4.2)
T
T is the window width and ∆t is the fast physics time step. The key point in terms of technical
implementation is that the model must be started 90 minutes prior to the actual starting time of the
assimilation run. The model is started from the 22:30 UTC first guess. The analysis increments for
U, V, P, T, QV, whose validity time is 00:00 UTC are added over 3 hours until at 1:30 the free forecast
starts. Then, two first guess data sets are written into the database. One at 1:30 UTC, which will be
used for starting the next 3 h assimilation run, and a second one at 3:00 UTC, which serves as input
21
Chapter 4. Analysis fields
Figure 4.2.: Time line for an ICON assimilation run starting at 00:00 UTC.
for the assimilation system itself. Thus in general, using the IAU method requires some care in terms
of reading and writing the right fields at the right times.
This method is not restricted to atmospheric fields, but also applicated to assimilated soil and surface
fields, specifically soil moisture W_SO, and snow quantities H_SNOW and FRESHSNW.
For experimental predictions and simulations with the ICON-D2 the initial state can be determined by
interpolation from the analysis of a driving model (normally ICON / ICON-EU, or also the IFS). With
interpolated initial conditions one should generally note that the calculated initial state is not very
well defined due to the difference in the horizontal and vertical resolution. Therefore, a settling period
occurs (spin-up, approx. 3–6 hours), during which the flow adjusts to the high-resolution topography1 .
KENDA-LETKF
The initial conditions will be generated with the data assimilation system KENDA (’Kilometer-scale
Ensemble Data Assimilation’, Schraff et al. (2016)) which is based on the method of the “Local Ensemble
Transform Kalman Filter” (LETKF, Hunt et al. (2007)). With this method, one can simultaneously and
consistently provide initial conditions both for the deterministic ICON-D2 in the form of a deterministic
analysis as well as for the ICON-D2-EPS (COSMO-D2-EPS before 2021-02-10) using a whole ensemble
of suitably disturbed analyses.
For the calculation of the analysis, the information is combined from the current observations and the
previous short-term forecast, in the case of the current configuration of KENDA this is a 1-hour forecast.
1 The digital filter initialization (DFI) by Lynch (1997) used earlier, is difficult to adapt for the 2-time-level based
integration used in ICON-D2 according to previous experience.
22
4.3. Initial state for the convection-resolving ICON-D2
The weighting of these components is based on the estimation of the respective uncertainties, where the
errors of the predictions in particular on the convective scale depends to a high degree on the situation
and weather. In KENDA, these forecasting uncertainties can be estimated with the help of an ensemble
from (currently) 40 appropriately slightly different 1-hour forecasts. The analysis procedure allows
not only to estimate the most likely actual state of the atmosphere, but also the analysis error. This
estimate is used in the generation of an entire ensemble of different analyses (with the same number
as the incoming prediction Ensemble) in such a way that the analysis ensemble mean corresponds to
the most likely current state and that the spread of the analysis ensemble corresponds to the estimated
analysis error. As a result, the estimation of the analysis error influences directly the spread of the
subsequent ensemble forecasts (in the data assimilation cycle or as an actual short-term forecast ICON-
D2-EPS), which in turn serves as a measure of the forecasting uncertainty. However, it should be
mentioned that the ensemble spread generally only describes random errors, but not the systematic
analysis errors or forecast errors.
To sufficiently take into account the uncertainty in the heat flux from the surface of the earth into
the atmosphere, additional explicit random errors in the sea surface temperature and soil moisture
are applied, so that ensemble members have a spread of 1 K or approx. 15 % relative soil moisture
(between wilting point and field capacity). Without these disturbances, the ensemble spread and thus
the estimate of the uncertainty of analysis and prediction in the planetary boundary layer would become
underestimated.
Because the analysis ensemble mean in the atmosphere is not a very well balanced model state, and,
as used as an initial state for a deterministic forecast, would lead to a slightly increased spin-up in
the first forecast hours, an additional undistorted model run (’Control Run’ or ’deterministic run’) is
determined. Based on this 1 hour forecast (’deterministic first guess’), the deviations of the observations
from this run, and the estimation of the forecast errors from the LETKF (in the form of the ’Kalman
Gain’ for the ensemble mean) a ’deterministic’ analysis is calculated. This serves as the initial condition
for the subsequent ’Control Run’ or for the actual deterministic short-term prediction.
The so-called ’control variables’ of the LETKF, i.e. the variables which are changed (’analysed’) by
LETKF are currently: 3-D wind components, temperature, specific humidity, cloud water, and cloud
ice on all model levels, as well as pressure at the bottom model surface. In the areas above it, the pressure
is adjusted in such a way that the entire analysis corrections (analysis increments) are hydrostatically
balanced. For all other model variables, the analysis is just the 1 hour forecast (’First Guess’).
Only conventional observations are currently assimilated in the LETKF (data from radio sondes, aircraft,
wind profilers, ground stations). To assimilate precipitation rates from radar data, the LETKF is
combined with ’Latent Heat Nudging’ (LHN). This is done by application of LHN between the LETKF
analysis steps and at the beginning of the predictions in the deterministic model run.
In order to have a sufficiently good forecast quality of precipitation, especially in the short range, the
use of radar-based precipitation information is essential for the determination of a reasonable initial
condition. Currently, quality proven products of near-ground precipitation rates are used in a temporal
resolution of 5 minutes and a horizontal resolution of 1 km × 1 km from the DWD radar network and
foreign radar stations. These data are aggregated to the ICON-D2 model grid and are brought in GRIB
format as precipitation analyses into the database. With the help of the “latent heat nudging” method
these radar precipitation data are assimilated during the forward integration of the (ICON) model into
the model state (Stephan et al., 2008). To do this, one determines temperature increments from the
ratio between observed and modelled precipitation as well as from model based latent heating rates. The
temperature changes take place while maintaining the relative humidity, whereby the specific humidity
is adjusted accordingly. The increments introduced influence the dynamics of the model in that the
model precipitation adjusts to the observation.
23
Chapter 4. Analysis fields
The soil moisture is adapted by relaxing the soil moisture index (SMI) of ICON-D2 towards the SMI
of ICON-EU (which uses a soil moisture analysis, see above).
Once a day there is an analysis of the sea surface temperature carried out. Based on the previous
analysis as ’first guess’, the new analysis is produced by using all observations from ships and buoys
of the previous 2 days with the aid of a correction procedure. In low data areas this is complemented
through the global analysis, based on the analysis by NCEP, which is also based on satellite data.
Furthermore, a snow depth analysis is carried out every 6 hours. It is based on a simple weighted
averaging of SYNOP snow depth observations. The weighting depends on the horizontal and vertical
distances the target grid points. In areas with low data density, an attempt is made to derive the snow
depth increments from SYNOP precipitation and temperatures.
24
5. Output fields of the ICON model:
General description
ICON output fields are exclusively available in the General Regularly-distributed Information in Binary
Form, 2nd edition (GRIB2), with the exception of meteogram data (NetCDF). GRIB is a bit-oriented
data storage format which was developed by World Meteorological Organization (WMO) to facilitate
the exchange of large volumes of gridded data between weather prediction centres.
augmented by a large number of additional metadata in order to uniquely describe the nature of the
data. Noteworthy examples of additional metadata are
• typeOfStatisticalProcessing, former known as stepType (instant, accum, avg, max, min, diff, rms,
sd, cov, . . . ), describing the statistical process used to calculate the field
https://community.wmo.int/activity-areas/wis/latest-version
For decoding and encoding GRIB2 messages, the DWD in general and ICON in particular makes use
of the ecCodes package developed by ECMWF. ecCodes includes both programming interfaces for
reading and writing GRIB2- (and the older GRIB1-) data in Fortran-, C-, and Python-programs and
command line tools for analysing and further processing of GRIB-fields. Examples for the latter are
for extensive information about the single grib fields. To see the meta data in ’pure’ form (i.e. only
so-called ’coded keys’ are displayed) and ordered by GRIB sections then
should be used.
https://confluence.ecmwf.int/display/ECC/ecCodes+Home
25
Chapter 5. Output fields of the ICON model: General description
Table 5.1 provides an overview of the available time-constant fields. As mentioned in Section 3.2, there
are two types of such variables. The one type is delivered as an external (invariant) field; such fields
are available from the database category CAT_NAME=$model_const_an_$suite. In the later tables
8.1, 9.2, and 10.2 they are denoted by ’invar’. The other type of variables (in particular DEPTH_LK,
HSURF, FR_LAND, FR_LAKE and Z0) is modified by ICON. Thus, the latter should not be taken from
the const_an database category, unless you definitely know what you are doing. For convenience, the
modified invariant fields (and some more) are stored in the forecast database categories for step s[h] = 0
(CAT_NAME=$model_$run_fc_$suite) (such variables are denoted by ’t=0’ in the above mentioned
tables).
See Section 12.1 for more details on the database categories and Section 12 for sample retrievals.
Table 5.1.: Time-constant fields or variables exclusively available for V V = 0 from the forecast databases
LL IntpType
Lev-Typ 1/2
ShortName
Discipline
stepType
Category
Number
Description
Unit
26
5.1. Available output fields
27
Chapter 5. Output fields of the ICON model: General description
Anisotropy of sub-gridscale
SSO_GAMMA 0/3/24 1/– inst – 1
orography
28
5.1. Available output fields
Table 5.2.: Hybrid multi-level forecast (V V > 0) and initialised analysis (V V = 0) products
LL IntpType
Lev-Typ 1/2
ShortName
Discipline
stepType
Category
Number
Description
Unit
CLC Cloud cover 0/6/22 150/150 inst %
DEN Density of moist air 0/3/10 150/150 inst – kg m−3
Buoyancy-production of TKE due
DTKE_CON 0/19/219 150/– inst – m2 s−3
to sub grid scale convection
3 for the time being, erroneously encoded as mixing ratios instead of specific quantities
29
Chapter 5. Output fields of the ICON model: General description
Table 5.3.: Hybrid multi-level forecast (V V > 0) and initialised analysis (V V = 0) products for forecasts
employing prognostic mineral dust
LL IntpType
Lev-Typ 1/2
ShortName
Discipline
stepType
Category
Number
Description
Unit
DUST_TOTAL Diagnostic total mass concentration
0/20/0 150/150 inst – kg m−3
_MC of mineral dust aerosol
30
5.1. Available output fields
The variables SAT_BSC_DUST and CEIL_BSC_DUST share the same grib key triplet for discipline,
parameterCategory and parameterNumber, but are available for different wavelengths, 532 and 1064 nm.
They can be distinguished by the grib key scaledValueOfFirstWavelength = 532/1064. The vari-
ables DUSTA, DUSTB and DUSTC, and likewise DUSTA0, DUSTB0 and DUSTC0 also share the same grib
key triplet, but represent different modes. They can be distinguished by the additional grib keys
numberOfModeOfDistribution, that indicates the number of modes used, and modeNumber = 1,2,3
that indicates which mode number (A: 1=fine, B: 2=medium, C: 3=coarse) is encoded.
Table 5.4.: Multi-level forecast (V V > 0) and initialised analysis (V V = 0) products interpolated to
pressure levels
LL IntpType
Lev-Typ 1/2
ShortName
Discipline
stepType
Category
Number
Description
Unit
CLC Cloud cover 0/6/22 100/– inst %
FI Geopotential 0/3/4 100/– inst m2 s−2
Vertical velocity in pressure
OMEGA 0/2/8 100/– inst Pa s−1
coordinates (ω = dp/dt)
4 for the time being, erroneously encoded as mixing ratios instead of specific quantities
31
Chapter 5. Output fields of the ICON model: General description
Table 5.5.: Multi-level forecast (V V > 0) and initialised analysis (V V = 0) products from forecasts
employing prognostic mineral dust
LL IntpType
Lev-Typ 1/2
ShortName
Discipline
stepType
Category
Number
Description
Unit
DUST_MAX Vertical maximum total mass
_TOTAL concentration of mineral dust 0/20/61 100/100 inst kg m−3
_MC_LAYER aerosol in a layer
Currently output for seven different layers is generated. The first layer starts at the surface (SFC).
Further above the bottom and top of each layer is a certain flight level (FL). Since the FLs are defined
as the heights of pressure levels in the ICAO standard atmosphere, these are encoded as the values
of scaledValueOfFirstFixedSurface and scaledValueOfSecondFixedSurface. The following levels
are used to define the layers: SFC: 101325, FL050: 84307, FL100: 69682, FL140: 59524, FL180: 50600,
FL250: 37601, FL350: 23842, FL450: 14748 Pa.
LL IntpType
Lev-Typ 1/2
ShortName
Discipline
stepType
Category
Number
Description
Unit
32
5.1. Available output fields
33
Chapter 5. Output fields of the ICON model: General description
34
5.1. Available output fields
SNOWLMT Height of snowfall limit above MSL 0/1/204 4/101 inst NNB m
Net short-wave radiation flux at
SOBS_RAD 0/4/9 1/– inst W m−2
surface (instantaneous)
35
Chapter 5. Output fields of the ICON model: General description
Maximum temperature at 2m
TMAX_2M 0/0/0 103/– max K
above ground
TQR Column integrated rain (grid scale) 0/1/45 1/– inst kg m−2
Column integrated snow (grid
TQS 0/1/46 1/– inst kg m−2
scale)
36
5.1. Available output fields
5 Note that the unit which is displayed, when inspecting the GRIB2 message with grib_dump is kg m−2 s−1 rather
than kg m−2 . Mathematically this is wrong, however, it is in accordance with the GRIB2 standard. To get the
mathematically correct unit for accumulated fields (typeOfStatisticalProcessing=1), the unit displayed by grib_dump
must be multiplied by s.
6 T_S is identical to T_SO at level 0. It will no longer be available in the future. Use T_SO(0) instead of T_S.
37
Chapter 5. Output fields of the ICON model: General description
Table 5.7.: Single-level forecast (V V > 0) and initialised analysis (V V = 0) products of the lake model
model
LL IntpType
Lev-Typ 1/2
ShortName
Discipline
stepType
Category
Number
Description
Unit
Shape factor with respect to the
C_T_LK temperature profile in the 1/2/10 162/166 inst – 1
thermocline
38
5.1. Available output fields
Table 5.8.: Single-level forecast (V V > 0) and initialised analysis (V V = 0) products of the ensemble
members with prognostic dust
LL IntpType
Lev-Typ 1/2
ShortName
Discipline
stepType
Category
Number
Description
Unit
DUST_TOTAL Column integrated mineral dust
0/20/1 10/- inst – kg m−2
_MC_VI aerosol
39
Chapter 5. Output fields of the ICON model: General description
Table 5.9.: Multi-level forecast (V V > 0) and initialised analysis (V V = 0) products of the soil model
LL IntpType
Lev-Typ 1/2
ShortName
Discipline
stepType
Category
Number
Description
Unit
T_SO Soil temperature 2/3/18 106/– inst K
Soil moisture integrated over
W_SO 2/3/20 106/106 inst kg m−2
individual soil layers (ice + liquid)
Soil temperature is defined at the soil depths given in Table 5.10 (column 2). Levels 1 to 8 define the
full levels of the soil model. A zero gradient condition is assumed between levels 0 and 1, meaning that
temperatures at the surface-atmosphere interface are set equal to the temperature at the first full level
depth (0.5 cm). Temperatures are prognosed for layers 1 to 7. At the lowermost layer (mid-level height
1458 cm) the temperature is fixed to the climatological average 2 m-temperature.
Soil moisture W_SO is prognosed for layers 1 to 6. In the two lowermost layers W_SO is filled with W_SO(6)
(zero gradient condition).
7 The grib key modeNumber is used as additional identifier to distinguish the different modes, with A: 1, B: 2, C: 3.
40
5.1. Available output fields
0 0.0
1 0.5 1 0.0 — 1.0
2 2.0 2 1.0 — 3.0
3 6.0 3 3.0 — 9.0
4 18.0 4 9.0 — 27.0
5 54.0 5 27.0 — 81.0
6 162.0 6 81.0 — 243.0
7 486.0 7 243.0 — 729.0
8 1458.0 8 729.0 — 2187.0
The soil moisture analysis (SMA) requires the following fields from the main run at 00 UTC. They are
written only by this run and from forecast hour 2 to 24. As a soil moisture analysis is made for the
global and the nest domain, these fields are available for both domains, but only on the native grid.
Table 5.11.: Fields for SMA from 00 UTC run for forecast hours 2 to 24.
LL IntpType
Lev-Typ 1/2
ShortName
Discipline
stepType
Category
Number
Description
Unit
ALHFL_BS Latent heat flux from bare soil 2/0/193 1/– avg – W m−2
ALHFL_PL Latent heat flux from plants 2/0/194 106/106 avg – W m−2
RSTOM Stomatal resistance 2/0/195 1/– inst – s m−1
The latent heat flux from plants is defined at the same soil layers as the soil moisture W_SO.
41
6. Extended description of available output
fields
In order to facilitate the selection and interpretation of fields and to guard against possible misinterpre-
tation or misusage, the following section provides a more thorough description of the available output
fields.
CEILING Ceiling is that height above MSL (in m), where the large scale cloud coverage (more
precise: scale and sub-scale, but without the convective contribution) first exceeds
50% when starting from ground.
CLCT_MOD Modified total cloud cover (0 ≤ CLCT_MOD ≤ 1). Used for visualisation purpose (i.e.
gray-scale figures) in the media. It is derived from CLC, neglecting cirrus clouds
if there are only high clouds present at a given grid point. The reason for this
treatment is that the general public does not regard transparent cirrus clouds as
‘real’ clouds.
CLDEPTH Modified cloud depth (0 ≤ CLDEPTH ≤ 1). Used for visualisation purpose (i.e. gray-
scale figures) in the media. A cloud reaching a vertical extent of 700 hPa or more,
has CLDEPTH= 1.
HBAS_CON Height of the convective cloud base in m above MSL. HBAS_CON is initialised with
−500 m at points where no convection is diagnosed.
HBAS_SC Height of the convective cloud base in m above MSL, but only the shallow convection
part is active.
HTOP_DC Height of the top of dry convection above MSL. It is the upper limit of dry thermals
rising from near the surface. At grid points without dry convection the value is zero,
or the surface height at points below MSL.
43
Chapter 6. Extended description of available output fields
HZEROCL Height of the 0◦ C isotherm above MSL. In case of multiple 0◦ C isotherms, HZEROCL
contains the uppermost one. If the temperature is below 0◦ C throughout the entire
atmospheric column, HZEROCL is set equal to the topography height (fill value).
Note that prior to 2019-07-30, HZEROCL contains the height of the lowermost 0◦ C
isotherm. At grid points where no 0◦ C isotherm could be diagnosed, a fill value of
0 is used.
SNOWLMT Height of snow fall limit above MSL. It is defined as the height where the wet
bulb temperature Tw first exceeds 1.3◦ C (scanning mode from top to bottom). If
this threshold is never reached within the entire atmospheric column, SNOWLMT is
undefined (GRIB2 bitmap).
es (Td ) = e
TMIN_2M Minimum temperature at 2 m above ground. Minima are collected over 6-hourly
intervals on all domains. (Prior to 2015-07-07 minima were collected over 3-hourly
intervals on the global grid.) Especially in situations with partial snow cover the
minimum temperature TMIN_2M of a grid point and time interval can be much lower
than any instantaneous 2 m temperature T_2M during that time interval. The reason
is that T_2M is defined as the average over all tiles of a grid point, while TMIN_2M is
based on the minimum temperature of all tiles.
TMAX_2M Same, but for maximum 2 m temperature.
VIS Near surface visibility in m.
VMAX_10M Maximum wind gust at 10 m above ground. It is diagnosed from the turbulence
state in the atmospheric boundary layer, including a potential enhancement by the
SSO parameterization over mountainous terrain. In the presence of deep convection,
it contains an additional contribution due to convective gusts.
Maxima are collected over hourly intervals on all domains. (Prior to 2015-07-07
maxima were collected over 3-hourly intervals on the global grid.)
44
6.4. Surface products
↑ dif
ASWDIFU_S =
b F sw,s Upward solar diffuse radiation flux at the surface, averaged over forecast time.
↓ dir
ASWDIR_S =
b F sw,s Downward solar direct radiation flux at the surface, averaged over forecast
time. This quantity is not directly provided by the radiation scheme. It is a-
posteriori diagnosed from the definition of the surface net shortwave radiation
flux Fsw,s
net
45
Chapter 6. Extended description of available output fields
RHO_SNOW Snow density in kg/m3 . It can vary between 50 kg/m3 for fresh snow and 400 kg/m3
for compacted old snow. At snow-free points over land and over water RHO_SNOW is
set to 0 kg/m3 . Note that prior to 2019-07-30 RHO_SNOW was set to 50 kg/m3 over
snow-free land points.
SDI2 The supercell detection index detects the mesocyclone of a supercell. It is based
on the product of a correlation between vertical velocity and vorticity and the local
vorticity Wicker et al. (2005).
T_ICE Ice temperature over sea-ice and frozen lake points. Melting ice has a temperature
of 273.15 K. Ice-free points over land, sea, and lakes are set to T_SO(0).
T_G Temperature at the atmosphere-surface interface, i.e. the temperature of those parts
of the ground which are in direct contact with the atmosphere. E.g. at snow-free
land points it is the temperature of the soil surface, whereas at snow covered land
points it is the temperature of the snow surface.
At snow-free land points T_G is equal to T_SO(0). Likewise, at open water points
T_G is equal to T_SO(0), and represents the sea-surface temperature SST (for more
details on SST see description of T_SO(0) in Section 6.5). At other grid points one
has
• T_G = T_SNOW+(1−f_snow)∗(T_SO(0)−T_SNOW) over (partially) snow covered
grid points. f_snow is the grid point fraction that is snow covered.
• T_G = T_ICE over frozen sea and fresh water lakes
TOT_PREC Total precipitation accumulated since model start. In global simulations (with and
without nests) it is
TOT_PREC = RAIN_GSP + SNOW_GSP + RAIN_CON + SNOW_CON,
whereas for ICON-D2 it is
TOT_PREC = RAIN_GSP + SNOW_GSP + GRAU_GSP + RAIN_CON + SNOW_CON.
T_SNOW Temperature of snow surface. At snow-free points (H_SNOW = 0), T_SNOW contains
the temperature of the soil surface T_SO(0).
WW Significant weather of the last hour. The predicted weather will be diagnosed hourly
at each model grid point and coded as a key number. The latter is called ww-code
and represents weather phenomena within the last hour. The interpretation of such
weather phenomena from raw model output relies on an independent post-processing
method. This technique applies a number of thresholding processes based on WMO
criteria. Therefore, a couple of ww-codes may differ from the direct model output
(e.g. ww-category snow vs. SNOW_GSP/SNOW_CON). Due to limitations in temporal
and spatial resolution, not all ww-codes as defined by the WMO criteria can be
determined. However, the simulated ww-code is able to take the following values:
no significant weather/ cloud cover (0, 1, 2, 3), fog (45, 48), drizzle (51, 53, 55, 56,
57), rain (61, 63, 65, 66, 67), solid precip not in showers (71, 73, 75, 77), showery
precip (liquid & solid) (80, 81, 82, 85, 86), thunderstorm (95, 96, 99 (only ICON-
D2)) (see also Table 6.1).
W_I Water content of interception layer, i.e. the amount of precipitation intercepted by
vegetation canopies. Over water points, W_I is set to 0.
W_SNOW Snow depth water equivalent in kg/m2 . Set to 0 above water surfaces and snow-free
land points.
W_CTMAX updraft velocity; delivered is the maximum value between ground and 10 km above
ground and during the last hour.
46
6.5. Soil products
Z0 Surface roughness length. Constant over land, where it depends only on the type of
land cover. I.e. it does not contain any contribution from subgrid-scale orography.
Over water, the roughness length usually varies with time. It is computed by the so
called Charnock-formula, which parameterizes the impact of waves on the roughness
length. Note that this field differs significantly from the external parameter field Z0
(see Table 3.2 or 8.1).
Table 6.1.: Weather interpretation (WW) code table for the ICON model. This table is a subset of the
WMO code table FM 94 BUFR/FM 95 CREX code table 0 20 003 – present weather. In the
case that none of the values provided in Table 6.1 is returned, the WW output contains the
total cloud cover, encoded in the following form: 0: clear sky 1: mainly clear 2: partly/generally
cloudy 3: cloudy/overcast.
RUNOFF_S Surface water runoff from interception and snow reservoir and from limited infiltra-
tion rate. Sum over forecast.
SOILTYP Characterizes the dominant soiltype in a grid cell. The soiltype is assumed to be
the same for all soil levels. Currently 9 soiltypes are distinguished and encoded by
1-digit integers 1-9. The mapping between these integer numbers and soiltype short
names is given in Table 6.2, together with some soil-dependent hydraulic parameters.
For the full list of hydraulic and thermal parameters, the reader is referred to Doms
et al. (2011).
47
Chapter 6. Extended description of available output fields
Table 6.2.: Mapping between the the soiltype index stored in the field SOILTYP and
soiltype short names. The hydraulic parameters porosity and field capacity,
currently used by ICON, are given in terms of volume fractions.
1 ice – –
2 rock – –
3 sand 0.364 0.196
4 sandyloam 0.445 0.260
5 loam 0.455 0.340
6 clayloam 0.475 0.370
7 clay 0.507 0.463
8 peat 0.863 0.763
9 sea water – –
T_SO Temperature of the soil and sea water. At land points T_SO(1:7) provides the prog-
nostic temperature of the soil. The full level depths at which the soil temperature
is defined are given in Table 5.10. The temperature at the uppermost level T_SO(0)
is not prognostic. It is rather set equal to the temperature at the first prognos-
tic level T_SO(1). The temperature at the lowermost level T_SO(8) is set to the
climatological 2 m temperature T_2M_CL.
At sea-points T_SO(0:7) provides the sea-surface temperature SST (same value at
all levels). So far, the SST in ICON is not prognostic. It is read from the analysis at
model start and is updated incrementally each day at 00 UTC based on its annual
climatological cycle.
Note that T_SO(0) does not necessarily represent the temperature at the interface
soil-atmosphere. I.e. over snow/ice covered surfaces, T_SO(0) represents the tem-
perature below snow/ice.
DUST_TOTAL_MC_VI Vertical integral of the mineral dust aerosol mass. Only calculated in fore-
casts with prognostic mineral dust.
TCOND_MAX The column integrated condensate (i.e. C, I, R, S, G), delivered is the maximum over
the last hour.
TCOND10_MAX As TCOND_MAX, but the vertical integration is restricted to heights above z(T =
−10◦ C).
TQx Column integrated water species x, derived from the 3D grid-scale prognostic quan-
tities Qx, with x ∈ {V, C, I, R, S, G}. TQx is based on the assumption that there would
be no sub-grid-scale variability. That assumption is particularly problematic for
precipitation generation, moist turbulence and radiation.
48
6.6. Vertical Integrals
TQx_DIA Total column integrated water species x, with x ∈ {C, I, V}. Takes into account
the sub-grid-scale variability that includes simple treatments of turbulent motion
and convective detrainment. These cloud variables attempt to represent all model
included physical processes. They are also consistent with the cloud cover variables
CLC, CLCT, CLCH, CLCM and CLCL.
TWATER This is just the sum over all TQx (but can be independently calculated).
UH_MAX Updraft helicity (i.e. the product of vertical velocity and vorticity) that is vertically
averaged between 2000 m and 8000 m above ground. Delivered is its maximum value
(either positive or negative) over the last hour.
UH_MAX_LOW same as UH_MAX but vertically averaged over the interval [0 m, 3000 m]
UH_MAX_MED same as UH_MAX but vertically averaged over the interval [2000 m, 5000 m]
VORW_CTMAX Vorticity, vertically averaged between the surface and 1500 m above ground. Deliv-
ered is its maximum value (either positive or negative) over the last hour.
49
7. Remarks on statistical processing and
horizontal interpolation
constitute time averages over the respective forecast time. The averaging process is performed from
forecast start (t0 = 0 s) till forecast end. Thus, time averaged fields which are written to the database
at t = ti contain averages for the elapsed time interval [t0 , ti ].
Let Ψ denote the instantaneous value of one of the above fields. The time average Ψ at time t stored
in the database is given as
1 t
Z
Ψ(t) = Ψ dt , for t > 0.
t 0
For t = 0, the average Ψ is equal to 0. If time averages are required for other time intervals [t1 , t2 ], with
t1 > 0, these can be computed as follows:
Z t2
1
Ψ(t2 − t1 ) = Ψ dt
t2 − t1 t1
Z t2 Z t1
1
= Ψ dt − Ψ dt
t2 − t1 0 0
1
= t2 Ψ(t2 ) − t1 Ψ(t1 )
t2 − t1
For this equation to work, it is of course necessary that the fields Ψ(t1 ) and Ψ(t2 ) are available from
the database.
The averaging process is fully reflected by the field’s GRIB2 metainfo. In order to check whether a
field contains the desired time average, it is advisable to check the content of the GRIB2 keys listed in
Table 7.1. I.e. productDefinitionTemplateNumber=8 indicates that the field in question is statistically
processed. The statistical process itself is specified by the key typeOfStatisticalProcessing. The
averaging interval (relative to the start of the forecast) is given by
51
Chapter 7. Remarks on statistical processing and horizontal interpolation
[forecastTime, forecastTime+lengthOfTimeRange].
Table 7.1.: List of GRIB2 keys which provide information about the averaging process
Octet(s) Key Value Meaning
as well as
ACCEMISS_DUST[ABC] ACCDRYDEPO_DUST[ABC]
ACCWETDEPO_GSP_DUST[ABC] ACCWETDEPO_CON_DUST[ABC]
ACCSEDIM_DUST[ABC]
for the forecasts including prognostic mineral dust
are accumulated over the respective forecast time. The accumulation process is performed from forecast
start (t0 = 0 s) till forecast end. Thus, fields which are written to the database at t = ti are accumulated
for the elapsed time interval [t0 , ti ].
Let Ψ denote the instantaneous value of one of the above fields. The accumulation Ψ̂ at time t stored
in the database is given as
Z t
Ψ̂(t) = Ψ dt , for t > 0.
0
For t = 0, the accumulation Ψ̂ is equal to 0. If accumulations are required for other time intervals
[t1 , t2 ], with t1 > 0, these can be computed as follows:
Z t2
Ψ̂(t2 − t1 ) = Ψ dt
t1
Zt2 Z t1
= Ψ dt − Ψ dt
0 0
= Ψ̂(t2 ) − Ψ̂(t1 )
52
7.1. Statistically processed output fields
For this equation to work, it is of course necessary that the fields Ψ̂(t1 ) and Ψ̂(t2 ) are available from
the database.
The accumulation process is fully reflected by the field’s GRIB2 metainfo. In order to check whether a
field contains the desired accumulation, it is advisable to check the content of the GRIB2 keys listed in
Table 7.2. I.e. productDefinitionTemplateNumber=8 indicates that the field in question is statistically
processed. The statistical process itself is specified by the key typeOfStatisticalProcessing. The
accumulation interval (relative to the start of the forecast) is given by
[forecastTime, forecastTime+lengthOfTimeRange].
Table 7.2.: List of GRIB2 keys which provide information about the accumulation process
Octet(s) Key Value Meaning
represent extreme values, which are collected over certain time intervals χ, starting from the beginning
of the forecast. The interval χ is variable dependent:
• χ = 6 h for TMAX_2M, TMIN_2M
• χ = 1 h for LPI_MAX, TCOND_MAX, TCOND10_MAX, UH_MAX, UH_MAX_LOW,
UH_MAX_MED, VMAX_10M, VORW_CTMAX, W_CTMAX, DBZ_CTMAX
• χ = 1, 3, or 6 h, depending on the forecast hour for LPI_CON_MAX
After χ hours of forecast the fields are re-initialized with 0 and the next χ-hourly collection phase is
started. This procedure is repeated till the end of the forecast.
Let Ψ denote the instantaneous value of one of the above fields. The maximum value Ψmax at time t
stored in the database is given as
Ψmax (t) = max(Ψ(t), Ψmax (t)) , for ti < t < ti + χ
Here, ti indicates the time when Ψmax was (re)-initialized the last time. For t = 0, the extreme value
Ψmax is equal to the instantaneous value Ψ.
53
Chapter 7. Remarks on statistical processing and horizontal interpolation
Please note: Even though a 6 hour time window is used for temperatures, the database contains
hourly, 2-hourly, etc. extreme temperatures. This is because the extreme temperatures are written to
the database hourly, irrespective of the start/end of the 6-hourly time windows. Example: Extreme
temperatures which are written into the database after a forecast time of 8 hours, contain extreme values
collected over the last 2 hours. On the other hand, extreme temperatures written into the database
after 12 hours contain values collected over the last 6 hours. Thus, when dealing with those fields it is
very important to check the GRIB2 keys listed in Table 7.3.
productDefinitionTemplateNumber=8 indicates that the field in question is statistically processed.
The statistical process itself is specified by the key typeOfStatisticalProcessing. The time interval
(relative to the start of the forecast) over which the extreme value collection was performed is given by
[forecastTime, forecastTime+lengthOfTimeRange]. Since the collection process is restarted every χ
hours, the key forecastTime can differ from 0.
Table 7.3.: List of GRIB2 keys which provide information about the extreme value process
Octet(s) Key Value Meaning
54
7.2. Technical Details of the Horizontal Interpolation
A small number of output fields is treated differently, with a nearest-neighbor interpolation (NNB). The
nearest neighbor algorithm selects the value of the nearest point and does not consider the values of
neighboring points at all, yielding a piecewise-constant interpolant.
55
8. Global output fields
ICON forecasts are performed multiple times a day with varying forecast periods. An overview of the
forecast runs, including its forecast period and output intervals is provided in Figure 8.1.
Figure 8.1.: Time span covered by the various global ICON forecasts which are launched every three hours.
Output on the native (triangular) grid ( ) and the regular grid ( ) is generally available until
forecast end, as indicated by the lenght of the two bars shown for each forecast run. Output
fields are available hourly up to V V = 78 h and 3-hourly for larger forecast times (for exceptions
see the following tables).
Main forecasts are performed 4 times a day at 0, 6, 12, 18 UTC, covering a forecast time span of 180 h
for the 0 und 12 UTC runs and 120 h for the 6 und 18 UTC runs. Prior to 2015-02-25 the 6 and 18
UTC runs were restricted to 78 h. Additional short-range forecasts are performed at 3, 9, 15 and 21
UTC. The forecast time covered by these runs is limited to 51 h since one main purpose of these runs
is to provide boundary data for the high resolution ICON-D2 runs from the ICON nest. See Chapter 9
for more details on the ICON nest and the available output fields.
57
Chapter 8. Global output fields
In general, all time-dependent output fields are available hourly up to V V = 78 h and 3-hourly for larger
forecast times2 . Please note that for ICON fields the time unit is minutes rather than hours, and thus
differs from the previously used global model GME (hours).
Output is available on two distinct horizontal grids:
• The native triangular grid with an average resolution of 13 km that covers the earth with 2949120
triangles
• a regular latitude-longitude grid with a resolution of ∆λ = ∆Φ = 0.25◦ (see table 8.1)
On the native grid most output fields are defined on triangle cell (circum-)centers, except for VN, which
is defined on cell edges. On the lat-lon grid, all fields are defined on cell centers.
global lat-lon
Table 8.1.: Summary of the latitude-longitude grid for ICON global output.
For details regarding the available fields, please see the tables below. A few remarks about the column
’Time range’: listed is the output time range in hours, followed by the output intervall (also in hours).
The time range is given for the longest runs (i.e. the 00 and 12 UTC runs); of course, for the shorter
runs at 03, 06, 09, 15, 18, 21 UTC, output is only available until the end of the forecast range.
Time range
EPS
Det.
ALB_SEAICE t=0 ✓ ✓
CLAT t=0 ✓ ✓
t=0 ✓ ✓
CLON t=0 ✓ ✓
t=0 ✓ ✓
Continued on next page
2 An exception here are the output fields VMAX_10M, U_10M and V_10M, which are available hourly throughout the forecast.
For the latter two this is because U_10M and V_10M are needed as input by the wave models.
58
8.1. Time-constant (external parameter) fields
C_T_LK t=0 ✓ ✓
DEPTH_LK t=0 ✓ ✓ ✓
ELAT t=0 ✓ ✓
t=0 ✓ ✓
ELON t=0 ✓ ✓
t=0 ✓ ✓
EVAP_PL t=0 ✓ ✓
FRESHSNW t=0 ✓ ✓
FR_ICE t=0 ✓ ✓
FR_LAKE t=0 ✓ ✓ ✓
FR_LAND t=0 ✓ ✓ ✓
t=0 ✓ ✓
HHL t=0 ✓ ✓ ✓ m
t=0 ✓ ✓ m
HSNOW_MAX t=0 ✓ ✓
HSURF t=0 ✓ ✓ ✓
t=0 ✓ ✓
H_ICE t=0 ✓ ✓
H_ML_LK t=0 ✓ ✓
H_SNOW t=0 ✓ ✓
LAI t=0 ✓ ✓ ✓
NDVIRATIO t=0 ✓ ✓
P t=0 ✓ ✓ m
PLCOV t=0 ✓ ✓ ✓
QC t=0 ✓ ✓ m
QI t=0 ✓ ✓ m
QR t=0 ✓ ✓ m
QS t=0 ✓ ✓ m
QV t=0 ✓ ✓ m
QV_S t=0 ✓ ✓
RHO_SNOW t=0 ✓ ✓
RLAT t=0 ✓ ✓
RLON t=0 ✓ ✓
ROOTDP t=0 ✓ ✓ ✓
Continued on next page
59
Chapter 8. Global output fields
Time range
EPS
Det.
AER_BC12 invar ✓ ✓ ✓
AER_DIF12 invar ✓ ✓ ✓
AER_DUST12 invar ✓ ✓ ✓
AER_MRAT invar ✓ ✓ ✓
AER_NI12 invar ✓ ✓ ✓
AER_ORG12 invar ✓ ✓ ✓
AER_SO412 invar ✓ ✓ ✓
AER_SS12 invar ✓ ✓ ✓
AER_UV12 invar ✓ ✓ ✓
Continued on next page
60
8.2. Multi-level fields on native hybrid vertical levels
CLAT invar ✓ ✓
CLON invar ✓ ✓
DEPTH_LK invar ✓ ✓ ✓
EMIS_RAD invar ✓ ✓ ✓
FOR_D invar ✓ ✓ ✓
FOR_E invar ✓ ✓ ✓
FR_LAKE invar ✓ ✓ ✓
FR_LAND invar ✓ ✓ ✓
FR_LUC invar ✓ ✓ ✓
HSURF invar ✓ ✓ ✓
LAI_MX invar ✓ ✓ ✓
NDVI_MAX invar ✓ ✓ ✓
PLCOV_MX invar ✓ ✓ ✓
ROOTDP invar ✓ ✓ ✓
RSMIN invar ✓ ✓ ✓
SOILTYP invar ✓ ✓ ✓
SSO_GAMMA invar ✓ ✓ ✓
SSO_SIGMA invar ✓ ✓ ✓
SSO_STDH invar ✓ ✓ ✓
SSO_THETA invar ✓ ✓ ✓
T_2M_CL invar ✓ ✓ ✓
Z0 invar ✓ ✓ ✓
61
Chapter 8. Global output fields
native
latlon
Time range
EPS
Det.
ShortName level type
CLC 0–78, 1 h ✓ ✓ ✓ m
81–180, 3 h ✓ ✓ ✓ m
DEN 0–78, 1 h ✓ ✓ m
81–180, 3 h ✓ ✓ m
DTKE_CON 0–180, 1 h ✓ ✓ m61-ke1
0–36, 1 h ✓ ✓ m_3
DTKE_HSH 0–180, 1 h ✓ ✓ m61-ke1
0–36, 1 h ✓ ✓ m_3
P 0–78, 1 h ✓ ✓ m
81–180, 3 h ✓ ✓ m
0–78, 1 h ✓ ✓ m39-ke1
81–180, 3 h ✓ ✓ m39-ke1
0–180, 6 h ✓ ✓ m
QC 0–78, 1 h ✓ ✓ m
81–180, 3 h ✓ ✓ m
0–78, 1 h ✓ ✓ m39-ke1
81–180, 3 h ✓ ✓ m39-ke1
0–180, 6 h ✓ ✓ m
QI 0–78, 1 h ✓ ✓ m
81–180, 3 h ✓ ✓ m
0–78, 1 h ✓ ✓ m39-ke1
81–180, 3 h ✓ ✓ m39-ke1
0–180, 6 h ✓ ✓ m
QR 0–78, 1 h ✓ ✓ m
81–180, 3 h ✓ ✓ m
QS 0–78, 1 h ✓ ✓ m
81–180, 3 h ✓ ✓ m
QV 0–78, 1 h ✓ ✓ m
81–180, 3 h ✓ ✓ m
0–78, 1 h ✓ ✓ m39-ke1
81–180, 3 h ✓ ✓ m39-ke1
Continued on next page
62
8.2. Multi-level fields on native hybrid vertical levels
0–180, 6 h ✓ ✓ m
T 0–78, 1 h ✓ ✓ m
81–180, 3 h ✓ ✓ m
0–78, 1 h ✓ ✓ m39-ke1
81–180, 3 h ✓ ✓ m39-ke1
0–180, 6 h ✓ ✓ m
TKE 0–180, 1 h ✓ ✓ m61-ke1
0–36, 1 h ✓ ✓ m_3
U 0–78, 1 h ✓ ✓ m
81–180, 3 h ✓ ✓ m
0–78, 1 h ✓ ✓ m39-ke1
81–180, 3 h ✓ ✓ m39-ke1
0–180, 6 h ✓ ✓ m
V 0–78, 1 h ✓ ✓ m
81–180, 3 h ✓ ✓ m
0–78, 1 h ✓ ✓ m39-ke1
81–180, 3 h ✓ ✓ m39-ke1
0–180, 6 h ✓ ✓ m
W 0–78, 1 h ✓ ✓ m
81–180, 3 h ✓ ✓ m
0–78, 1 h ✓ ✓ m39-ke1
81–180, 3 h ✓ ✓ m39-ke1
Forecasts employing prognostic mineral dust also provide the “standard” meteorological variables.
However, they are only made available on the native grid, there is no interpolation to lat-lon.
Time range
EPS
Det.
AOD_DUST 0–48, 6 h ✓ ✓ m
60–180, 12 h ✓ ✓ m
Continued on next page
63
Chapter 8. Global output fields
0–48, 6 h ✓ ✓ m
60–180, 12 h ✓ ✓ m
CEIL_BSC_DUST
(1064 nm) 0–48, 6 h ✓ ✓ m
60–180, 12 h ✓ ✓ m
0–48, 6 h ✓ ✓ m
60–180, 12 h ✓ ✓ m
CEIL_BSC_DUST
(532 nm) 0–48, 6 h ✓ ✓ m
60–180, 12 h ✓ ✓ m
0–48, 6 h ✓ ✓ m
60–180, 12 h ✓ ✓ m
CLC 0–78, 1 h ✓ ✓ m39-ke1
81–180, 3 h ✓ ✓ m39-ke1
DEN 0–78, 1 h ✓ ✓ m39-ke1
81–180, 3 h ✓ ✓ m39-ke1
DTKE_CON 0–48, 1 h ✓ ✓ m61-ke1
0–36, 1 h ✓ ✓ m_3
DTKE_HSH 0–48, 1 h ✓ ✓ m61-ke1
0–36, 1 h ✓ ✓ m_3
DUSTA 0–48, 6 h ✓ ✓ m
60–180, 12 h ✓ ✓ m
0–48, 6 h ✓ ✓ m
60–180, 12 h ✓ ✓ m
DUSTA0 0–48, 6 h ✓ ✓ m
60–180, 12 h ✓ ✓ m
0–48, 6 h ✓ ✓ m
60–180, 12 h ✓ ✓ m
DUSTB 0–48, 6 h ✓ ✓ m
60–180, 12 h ✓ ✓ m
0–48, 6 h ✓ ✓ m
60–180, 12 h ✓ ✓ m
DUSTB0 0–48, 6 h ✓ ✓ m
60–180, 12 h ✓ ✓ m
Continued on next page
64
8.2. Multi-level fields on native hybrid vertical levels
0–48, 6 h ✓ ✓ m
60–180, 12 h ✓ ✓ m
DUSTC 0–48, 6 h ✓ ✓ m
60–180, 12 h ✓ ✓ m
0–48, 6 h ✓ ✓ m
60–180, 12 h ✓ ✓ m
DUSTC0 0–48, 6 h ✓ ✓ m
60–180, 12 h ✓ ✓ m
0–48, 6 h ✓ ✓ m
60–180, 12 h ✓ ✓ m
DUST_TOTAL_MC 0–180, 3 h ✓ ✓ m
0–180, 3 h ✓ ✓ m
P 0–78, 1 h ✓ ✓ m
81–180, 3 h ✓ ✓ m
0–180, 6 h ✓ ✓ m
QC 0–78, 1 h ✓ ✓ m
81–180, 3 h ✓ ✓ m
0–180, 6 h ✓ ✓ m
QI 0–78, 1 h ✓ ✓ m
81–180, 3 h ✓ ✓ m
0–180, 6 h ✓ ✓ m
QR 0–78, 1 h ✓ ✓ m
81–180, 3 h ✓ ✓ m
QS 0–78, 1 h ✓ ✓ m
81–180, 3 h ✓ ✓ m
QV 0–78, 1 h ✓ ✓ m
81–180, 3 h ✓ ✓ m
0–180, 6 h ✓ ✓ m
SAT_BSC_DUST
(1064 nm) 0–48, 6 h ✓ ✓ m
60–180, 12 h ✓ ✓ m
0–48, 6 h ✓ ✓ m
60–180, 12 h ✓ ✓ m
SAT_BSC_DUST
Continued on next page
65
Chapter 8. Global output fields
There are several ’level types’ for output on pressure levels. In the following table 8.6 they are denoted
as
• p2: output on pressure levels 30, 50, 70, 100, 150, 200, 250, 300, 400, 500, 600, 700, 800, 850, 900,
925, 950, 1000 hPa
• p3: output on pressure levels 5, 10, 30, 50, 70, 100, 125, 150, 175, 200, 225, 250, 275, 300, 350,
400, 450, 500, 550, 600, 650, 700, 750, 775, 800, 825, 850, 875, 900, 925, 950, 975, 1000 hPa
• pe1: output on pressure levels 1, 2, 5, 10, 30, 50, 70, 100, 200 hPa
• pe3: output on pressure levels 300, 400, 500, 700, 850, 900, 925, 950, 1000 hPa
• pd7: output of the maximum in a layer defined by certain flight levels (FL), which correspond to
pressure levels in the ICAO standard atmosphere
SFC: 101325, FL050: 84307, FL100: 69682, FL140: 59524, FL180: 50600, FL250: 37601, FL350:
23842, FL450: 14748 Pa
66
8.3. Multi-level fields interpolated to pressure levels
native
latlon
Time range
EPS
Det.
ShortName level type
CLC 0–78, 1 h ✓ ✓ p3
81–180, 3 h ✓ ✓ p3
FI 0–78, 1 h ✓ ✓ p2
81–180, 3 h ✓ ✓ p2
0–78, 1 h ✓ ✓ p3, p4
81–180, 3 h ✓ ✓ p3, p4
0–75, 1 h ✓ ✓ pe3
75–180, 3 h ✓ ✓ pe3
0–180, 6 h ✓ ✓ pe2
0–72, 6 h ✓ ✓ pe1
72–180, 12 h ✓ ✓ pe1
OMEGA 0–78, 1 h ✓ ✓ p3
81–180, 3 h ✓ ✓ p3
0–75, 1 h ✓ ✓ pe4
75–180, 3 h ✓ ✓ pe4
RELHUM 0–78, 1 h ✓ ✓ p2
81–180, 3 h ✓ ✓ p2
0–78, 1 h ✓ ✓ p3
81–180, 3 h ✓ ✓ p3
0–75, 1 h ✓ ✓ pe3
75–180, 3 h ✓ ✓ pe3
0–180, 6 h ✓ ✓ pe2
0–72, 6 h ✓ ✓ pe1
72–180, 12 h ✓ ✓ pe1
T 0–78, 1 h ✓ ✓ p2
81–180, 3 h ✓ ✓ p2
0–78, 1 h ✓ ✓ p3, p4
81–180, 3 h ✓ ✓ p3, p4
0–75, 1 h ✓ ✓ pe3
75–180, 3 h ✓ ✓ pe3
Continued on next page
67
Chapter 8. Global output fields
0–180, 6 h ✓ ✓ pe2
0–72, 6 h ✓ ✓ pe1
72–180, 12 h ✓ ✓ pe1
U 0–78, 1 h ✓ ✓ p2
81–180, 3 h ✓ ✓ p2
0–78, 1 h ✓ ✓ p3, p4
81–180, 3 h ✓ ✓ p3, p4
0–75, 1 h ✓ ✓ pe3
75–180, 3 h ✓ ✓ pe3
0–180, 6 h ✓ ✓ pe2
0–72, 6 h ✓ ✓ pe1
72–180, 12 h ✓ ✓ pe1
V 0–78, 1 h ✓ ✓ p2
81–180, 3 h ✓ ✓ p2
0–78, 1 h ✓ ✓ p3, p4
81–180, 3 h ✓ ✓ p3, p4
0–75, 1 h ✓ ✓ pe3
75–180, 3 h ✓ ✓ pe3
0–180, 6 h ✓ ✓ pe2
0–72, 6 h ✓ ✓ pe1
72–180, 12 h ✓ ✓ pe1
Time range
EPS
Det.
DUST_MAX_TOTAL
_MC_LAYER 0–180, 3 h ✓ ✓ pd7
FI 0–78, 1 h ✓ ✓ p2
81–180, 3 h ✓ ✓ p2
0–75, 1 h ✓ ✓ pe3
Continued on next page
68
8.3. Multi-level fields interpolated to pressure levels
75–180, 3 h ✓ ✓ pe3
0–180, 6 h ✓ ✓ pe2
0–72, 6 h ✓ ✓ pe1
72–180, 12 h ✓ ✓ pe1
OMEGA 0–75, 1 h ✓ ✓ pe4
75–180, 3 h ✓ ✓ pe4
RELHUM 0–78, 1 h ✓ ✓ p2
81–180, 3 h ✓ ✓ p2
0–75, 1 h ✓ ✓ pe3
75–180, 3 h ✓ ✓ pe3
0–180, 6 h ✓ ✓ pe2
0–72, 6 h ✓ ✓ pe1
72–180, 12 h ✓ ✓ pe1
T 0–78, 1 h ✓ ✓ p2
81–180, 3 h ✓ ✓ p2
0–75, 1 h ✓ ✓ pe3
75–180, 3 h ✓ ✓ pe3
0–180, 6 h ✓ ✓ pe2
0–72, 6 h ✓ ✓ pe1
72–180, 12 h ✓ ✓ pe1
U 0–78, 1 h ✓ ✓ p2
81–180, 3 h ✓ ✓ p2
0–75, 1 h ✓ ✓ pe3
75–180, 3 h ✓ ✓ pe3
0–180, 6 h ✓ ✓ pe2
0–72, 6 h ✓ ✓ pe1
72–180, 12 h ✓ ✓ pe1
V 0–78, 1 h ✓ ✓ p2
81–180, 3 h ✓ ✓ p2
0–75, 1 h ✓ ✓ pe3
75–180, 3 h ✓ ✓ pe3
0–180, 6 h ✓ ✓ pe2
0–72, 6 h ✓ ✓ pe1
72–180, 12 h ✓ ✓ pe1
Continued on next page
69
Chapter 8. Global output fields
native
latlon
Time range
EPS
Det.
ShortName
AER_BC12 invar, 0 h ✓ ✓ ✓
AER_DIF12 invar, 0 h ✓ ✓ ✓
AER_DUST12 invar, 0 h ✓ ✓ ✓
AER_MRAT invar, 0 h ✓ ✓ ✓
AER_NI12 invar, 0 h ✓ ✓ ✓
AER_ORG12 invar, 0 h ✓ ✓ ✓
AER_SO412 invar, 0 h ✓ ✓ ✓
AER_SS12 invar, 0 h ✓ ✓ ✓
AER_UV12 invar, 0 h ✓ ✓ ✓
ALB_RAD 0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
ALHFL_BS 2–24, 1 h ✓ ✓
ALHFL_PL 2–24, 1 h ✓ ✓
ALHFL_S 0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
6–180, 6 h ✓ ✓
APAB_S 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
ASHFL_S 0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
6–180, 6 h ✓ ✓
ASOB_S 0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
Continued on next page
70
8.4. Single-level fields
ASOB_S_CS 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
ASOB_T 0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
ASWDIFD_S 0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
3–180, 3 h ✓ ✓
ASWDIFU_S 0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
6–180, 6 h ✓ ✓
ASWDIR_S 0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
3–180, 3 h ✓ ✓
ATHB_S 0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
ATHB_T 0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
AUMFL_S 0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
AVMFL_S 0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
CAPE_CON 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
CAPE_ML 0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
CIN_ML 0–78, 1 h ✓ ✓
Continued on next page
71
Chapter 8. Global output fields
81–180, 3 h ✓ ✓
CLAT invar, 0 h ✓ ✓
CLCH 0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
CLCL 0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
CLCM 0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
CLCT 0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
CLCT_MOD 0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
0–180, 6 h ✓ ✓
CLDEPTH 0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
CLON invar, 0 h ✓ ✓
C_T_LK 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
DEPTH_LK invar, 0 h ✓ ✓ ✓
EMIS_RAD invar, 0 h ✓ ✓ ✓
FOR_D invar, 0 h ✓ ✓ ✓
FOR_E invar, 0 h ✓ ✓ ✓
FRESHSNW 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
FR_ICE 0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
Continued on next page
72
8.4. Single-level fields
0–72, 6 h ✓ ✓
72–180, 12 h ✓ ✓
FR_LAKE invar, 0 h ✓ ✓ ✓
FR_LAND invar, 0 h ✓ ✓ ✓
FR_LUC invar, 0 h ✓ ✓ ✓
HBAS_CON 0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
HSURF invar, 0 h ✓ ✓ ✓
HTOP_CON 0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
HTOP_DC 0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
HZEROCL 0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
H_ICE 0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
0–72, 6 h ✓ ✓
72–180, 12 h ✓ ✓
H_ML_LK 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
H_SNOW 0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
0–72, 6 h ✓ ✓
72–180, 12 h ✓ ✓
LAI_MX invar, 0 h ✓ ✓ ✓
LPI_CON_MAX 0–48, 1 h ✓ ✓
0–48, 1 h ✓ ✓
NDVI_MAX invar, 0 h ✓ ✓ ✓
PLCOV_MX invar, 0 h ✓ ✓ ✓
PMSL 0–78, 1 h ✓ ✓ ✓
Continued on next page
73
Chapter 8. Global output fields
81–180, 3 h ✓ ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
PS 0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
0–180, 6 h ✓ ✓
QV_2M 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
QV_S 0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
RAIN_CON 0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
RAIN_GSP 0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
RELHUM_2M 0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
0–180, 6 h ✓ ✓
RHO_SNOW 0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
ROOTDP invar, 0 h ✓ ✓ ✓
RSMIN invar, 0 h ✓ ✓ ✓
RSTOM 2–24, 1 h ✓ ✓
RUNOFF_G 0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
RUNOFF_S 0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
SNOW_CON 0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
Continued on next page
74
8.4. Single-level fields
SNOW_GSP 0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
SOBS_RAD 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
SOILTYP invar, 0 h ✓ ✓ ✓
SSO_GAMMA invar, 0 h ✓ ✓ ✓
SSO_SIGMA invar, 0 h ✓ ✓ ✓
SSO_STDH invar, 0 h ✓ ✓ ✓
SSO_THETA invar, 0 h ✓ ✓ ✓
TCH 0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
TCM 0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
TD_2M 0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
THBS_RAD 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
TMAX_2M 0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
6–180, 6 h ✓ ✓
TMIN_2M 0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
6–180, 6 h ✓ ✓
TOT_PREC 0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
0–75, 1 h ✓ ✓
Continued on next page
75
Chapter 8. Global output fields
75–180, 3 h ✓ ✓
TQC 0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
TQC_DIA 0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
TQI 0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
TQI_DIA 0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
TQR 0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
TQS 0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
TQV 0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
T_2M 0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
T_2M_CL invar, 0 h ✓ ✓ ✓
T_BOT_LK 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
T_G 0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
T_ICE 0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
Continued on next page
76
8.4. Single-level fields
T_MNW_LK 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
T_S 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
T_SNOW 0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
0–72, 6 h ✓ ✓
72–180, 12 h ✓ ✓
T_SO 0–180, 6 h ✓ ✓
T_WML_LK 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
U_10M 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
0–180, 1 h ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
VMAX_10M 1–180, 1 h ✓ ✓
0–180, 1 h ✓ ✓
1–180, 1 h ✓ ✓
V_10M 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
0–180, 1 h ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
WW 0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
W_I 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
W_SNOW 0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
W_SO 0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
Continued on next page
77
Chapter 8. Global output fields
W_SO_ICE 81–180, 3 h ✓ ✓ ✓
Z0 invar, 0 h ✓ ✓ ✓
0–78, 1 h ✓ ✓ ✓
81–180, 3 h ✓ ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
Table 8.9.: Single-level fields from forecasts employing prognostic mineral dust
native
latlon
Time range
EPS
Det.
ShortName
ACCDRYDEPO_DUSTA 0–180, 3 h ✓ ✓
0–180, 3 h ✓ ✓
ACCDRYDEPO_DUSTB 0–180, 3 h ✓ ✓
0–180, 3 h ✓ ✓
ACCDRYDEPO_DUSTC 0–180, 3 h ✓ ✓
0–180, 3 h ✓ ✓
ACCEMISS_DUSTA 0–180, 3 h ✓ ✓
0–180, 3 h ✓ ✓
ACCEMISS_DUSTB 0–180, 3 h ✓ ✓
0–180, 3 h ✓ ✓
ACCEMISS_DUSTC 0–180, 3 h ✓ ✓
0–180, 3 h ✓ ✓
ACCSEDIM_DUSTA 0–180, 3 h ✓ ✓
0–180, 3 h ✓ ✓
ACCSEDIM_DUSTB 0–180, 3 h ✓ ✓
0–180, 3 h ✓ ✓
ACCSEDIM_DUSTC 0–180, 3 h ✓ ✓
0–180, 3 h ✓ ✓
ACCWETDEPO
_CON_DUSTA 0–180, 3 h ✓ ✓
Continued on next page
78
8.4. Single-level fields
0–180, 3 h ✓ ✓
ACCWETDEPO
_CON_DUSTB 0–180, 3 h ✓ ✓
0–180, 3 h ✓ ✓
ACCWETDEPO
_CON_DUSTC 0–180, 3 h ✓ ✓
0–180, 3 h ✓ ✓
ACCWETDEPO
_GSP_DUSTA 0–180, 3 h ✓ ✓
0–180, 3 h ✓ ✓
ACCWETDEPO
_GSP_DUSTB 0–180, 3 h ✓ ✓
0–180, 3 h ✓ ✓
ACCWETDEPO
_GSP_DUSTC 0–180, 3 h ✓ ✓
0–180, 3 h ✓ ✓
AER_BC12 invar, 0 h ✓ ✓
AER_DIF12 invar, 0 h ✓ ✓
AER_DUST12 invar, 0 h ✓ ✓
AER_NI12 invar, 0 h ✓ ✓
AER_ORG12 invar, 0 h ✓ ✓
AER_SO412 invar, 0 h ✓ ✓
AER_SS12 invar, 0 h ✓ ✓
AER_UV12 invar, 0 h ✓ ✓
ALB_RAD 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
ALHFL_BS 2–24, 3 h ✓ ✓
ALHFL_PL 2–24, 3 h ✓ ✓
ALHFL_S 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
6–180, 6 h ✓ ✓
ASHFL_S 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
6–180, 6 h ✓ ✓
Continued on next page
79
Chapter 8. Global output fields
ASOB_S 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
ASOB_S_CS 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
ASOB_T 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
ASWDIFD_S 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
3–180, 3 h ✓ ✓
ASWDIFU_S 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
6–180, 6 h ✓ ✓
ASWDIR_S 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
3–180, 3 h ✓ ✓
ATHB_S 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
ATHB_T 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
AUMFL_S 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
AVMFL_S 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
CAPE_CON 0–78, 1 h ✓ ✓
Continued on next page
80
8.4. Single-level fields
81–180, 3 h ✓ ✓
CAPE_ML 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
CLAT invar, 0 h ✓ ✓
CLCH 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
CLCL 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
CLCM 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
CLCT 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
CLCT_MOD 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
0–180, 6 h ✓ ✓
CLDEPTH 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
CLON invar, 0 h ✓ ✓
C_T_LK 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
DEPTH_LK invar, 0 h ✓ ✓
DUST_TOTAL_MC_VI 0–180, 3 h ✓ ✓
0–180, 3 h ✓ ✓
EMIS_RAD invar, 0 h ✓ ✓
Continued on next page
81
Chapter 8. Global output fields
FOR_D invar, 0 h ✓ ✓
FOR_E invar, 0 h ✓ ✓
FRESHSNW 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
FR_ICE 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
0–72, 6 h ✓ ✓
72–180, 12 h ✓ ✓
FR_LAKE invar, 0 h ✓ ✓
FR_LAND invar, 0 h ✓ ✓
FR_LUC invar, 0 h ✓ ✓
HBAS_CON 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
HSURF invar, 0 h ✓ ✓
HTOP_CON 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
HTOP_DC 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
HZEROCL 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
H_ICE 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
0–72, 6 h ✓ ✓
72–180, 12 h ✓ ✓
H_ML_LK 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
H_SNOW 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
0–72, 6 h ✓ ✓
72–180, 12 h ✓ ✓
Continued on next page
82
8.4. Single-level fields
LAI_MX invar, 0 h ✓ ✓
LPI_CON_MAX 0–48, 1 h ✓ ✓
0–48, 1 h ✓ ✓
NDVI_MAX invar, 0 h ✓ ✓
NDVI_MRAT invar, 0 h ✓ ✓
PLCOV_MX invar, 0 h ✓ ✓
PMSL 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
PS 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
0–180, 6 h ✓ ✓
QV_S 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
RAIN_CON 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
RAIN_GSP 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
RELHUM_2M 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
0–180, 6 h ✓ ✓
RHO_SNOW 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
ROOTDP invar, 0 h ✓ ✓
RSMIN invar, 0 h ✓ ✓
RSTOM 2–24, 3 h ✓ ✓
RUNOFF_G 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
RUNOFF_S 0–78, 1 h ✓ ✓
Continued on next page
83
Chapter 8. Global output fields
81–180, 3 h ✓ ✓
SKD invar, 0 h ✓ ✓
SNOW_CON 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
SNOW_GSP 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
SOBS_RAD 0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
SOILTYP invar, 0 h ✓ ✓
SSO_GAMMA invar, 0 h ✓ ✓
SSO_OROMAX invar, 0 h ✓ ✓
SSO_OROMIN invar, 0 h ✓ ✓
SSO_SIGMA invar, 0 h ✓ ✓
SSO_STDH invar, 0 h ✓ ✓
SSO_THETA invar, 0 h ✓ ✓
TAOD_DUST 0–180, 3 h ✓ ✓
0–180, 3 h ✓ ✓
TCH 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
TCM 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
TD_2M 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
THBS_RAD 0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
TMAX_2M 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
6–180, 6 h ✓ ✓
Continued on next page
84
8.4. Single-level fields
TMIN_2M 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
6–180, 6 h ✓ ✓
TOT_PREC 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
TQC 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
TQC_DIA 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
TQI 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
TQI_DIA 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
TQR 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
TQS 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
TQV 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
T_2M 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
T_2M_CL invar, 0 h ✓ ✓
T_BOT_LK 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
Continued on next page
85
Chapter 8. Global output fields
T_G 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
T_ICE 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
T_MNW_LK 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
T_S 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
T_SEA invar, 0 h ✓ ✓
T_SNOW 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
0–72, 6 h ✓ ✓
72–180, 12 h ✓ ✓
T_SO 0–180, 6 h ✓ ✓
T_WML_LK 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
USTAR 0–180, 3 h ✓ ✓
USTAR_THRES 0–180, 3 h ✓ ✓
U_10M 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
VMAX_10M 1–180, 1 h ✓ ✓
1–180, 1 h ✓ ✓
V_10M 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
WW 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
W_I 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
Continued on next page
86
8.5. Soil-specific multi-level fields
W_SNOW 0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
W_SO 0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
W_SO_ICE 81–180, 3 h ✓ ✓
Z0 invar, 0 h ✓ ✓
0–78, 1 h ✓ ✓
81–180, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–180, 3 h ✓ ✓
fr_cloa invar, 0 h ✓ ✓
fr_hcla invar, 0 h ✓ ✓
fr_lcla invar, 0 h ✓ ✓
fr_loam invar, 0 h ✓ ✓
fr_lsan invar, 0 h ✓ ✓
fr_sand invar, 0 h ✓ ✓
fr_scla invar, 0 h ✓ ✓
fr_sclo invar, 0 h ✓ ✓
fr_sicl invar, 0 h ✓ ✓
fr_silc invar, 0 h ✓ ✓
fr_silo invar, 0 h ✓ ✓
fr_silt invar, 0 h ✓ ✓
fr_sloa invar, 0 h ✓ ✓
fr_udef invar, 0 h ✓ ✓
87
Chapter 8. Global output fields
native
latlon
Time range
EPS
Det.
ShortName level type
Table 8.11.: Soil-specific multi-level fields from forecasts employing prognostic mineral dust
native
latlon
Time range
EPS
Det.
88
9. EU Nest output fields
This section contains a list of output fields that are available with the launch of the ICON-EU nest. See
Fig. 2.7 (top) on page 9 for details regarding the nest location and extent. In the forecasts employing
prognostic mineral dust a bigger nest domain ICON-EU-NA2 is used. For details see Fig. 2.7 (bottom)
on page 9.
Forecasts on the EU-nest or respectivly the EU-NA2 -nest are performed multiple times a day with
varying forecast periods. Forecasts reaching out to 120 h are performed at 00, 06, 12, and 18 UTC.
Additional short-range forecasts reaching out to 51 h are performed at 03, 09, 15 and 21 UTC. Its
main purpose is to provide boundary data for the high resolution ICON-D2 (formerly: COSMO-D2
or COSMO-DE) runs. A schematic overview of the various forecasts, including its forecast period and
output intervals is provided in Figure 9.1.
Output is available on two distinct horizontal grids:
• a native triangular grid with an average resolution of 6.5 km, and
• a regular latitude-longitude grid with a resolution of ∆λ = ∆Φ = 0.0625◦ .
See Table 9.1 for a summary.
Output on the native (triangular) grid is hourly to 51 h, and every 6 hours for verification from
forecast time ≥ 54 h until the forecast end at 120 h. Output on the regular grid is hourly to 78 h,
and every 3 hours until forecast end. See also Figure 9.1. Output of the 10m wind U_10M, V_10M,
VMAX_10M and the solar radiation ASWDIR_S and ASWDIFD_S on the regular grid is hourly until
the end of the forecast.
Output by the forecasts employing prognostic mineral dust is available on the native (triangular) grid
only.
In the subsequent tables the availability of specific fields on the native grid, on the lat-lon grid, or on
both grids is denoted. A few remarks about the column ’Time range’: listed is the output time range
in hours, followed by the output intervall (also in hours). The time range is given for the longest runs
(i.e. the 00 and 12 UTC runs); of course, for the shorter runs at 03, 06, 09, 15, 18, 21 UTC, output is
only available until the end of the forecast range.
EU nest lat-lon
Table 9.1.: Summary of the latitude-longitude grid for the ICON-EU nest output.
89
Chapter 9. EU Nest output fields
Figure 9.1.: Time span covered by the various EU-nest (or EU-NA2 -nest) forecasts which are launched every
three hours. Output on the native (triangular) grid ( ) and the regular grid ( ) is generally
available until forecast end, as indicated by the lenght of the two bars shown for each forecast
run. Output on the native grid is available hourly to 51 h, and every 6 hours for later forecast
times (forecast time ≥ 54 h). Output on the regular grid is available hourly to 78 h, and every
3 hours for later forecast times.
90
9.1. Time-constant (external parameter) fields
native
latlon
Time range
EPS
Det.
ShortName level type
ALB_SEAICE t=0 ✓ ✓
CLAT t=0 ✓ ✓
t=0 ✓ ✓
CLON t=0 ✓ ✓
t=0 ✓ ✓
DEPTH_LK t=0 ✓ ✓ ✓
ELAT t=0 ✓ ✓
t=0 ✓ ✓
ELON t=0 ✓ ✓
t=0 ✓ ✓
EVAP_PL t=0 ✓ ✓
FR_LAKE t=0 ✓ ✓ ✓
FR_LAND t=0 ✓ ✓ ✓
t=0 ✓ ✓
HHL t=0 ✓ ✓ ✓ m
t=0 ✓ ✓ m_5
HSNOW_MAX t=0 ✓ ✓
HSURF t=0 ✓ ✓ ✓
t=0 ✓ ✓
H_SNOW t=0 ✓ ✓
LAI t=0 ✓ ✓ ✓
PLCOV t=0 ✓ ✓ ✓
RLAT t=0 ✓ ✓
RLON t=0 ✓ ✓
ROOTDP t=0 ✓ ✓ ✓
SMI t=0 ✓ ✓ soil
SNOAG t=0 ✓ ✓
SOILTYP t=0 ✓ ✓ ✓
Continued on next page
91
Chapter 9. EU Nest output fields
Z0 t=0 ✓ ✓
Table 9.3.: Time-constant (external parameter) fields in forecasts employing prognostic mineral dust
native
latlon
Time range
EPS
Det.
ShortName level type
ALB_SEAICE t=0 ✓ ✓
CLAT t=0 ✓ ✓
t=0 ✓ ✓
CLON t=0 ✓ ✓
t=0 ✓ ✓
DEPTH_LK t=0 ✓ ✓
ELAT t=0 ✓ ✓
t=0 ✓ ✓
ELON t=0 ✓ ✓
t=0 ✓ ✓
EVAP_PL t=0 ✓ ✓
FR_LAKE t=0 ✓ ✓
FR_LAND t=0 ✓ ✓
t=0 ✓ ✓
HHL t=0 ✓ ✓ m
t=0 ✓ ✓ m
HSNOW_MAX t=0 ✓ ✓
HSURF t=0 ✓ ✓
t=0 ✓ ✓
H_SNOW t=0 ✓ ✓
t=0 ✓ ✓
LAI t=0 ✓ ✓
PLCOV t=0 ✓ ✓
ROOTDP t=0 ✓ ✓
SMI t=0 ✓ ✓ soil
Continued on next page
92
9.2. Multi-level fields on native hybrid vertical levels
SNOAG t=0 ✓ ✓
SOILTYP t=0 ✓ ✓
SSO_STDH t=0 ✓ ✓
t=0 ✓ ✓
Z0 t=0 ✓ ✓
t=0 ✓ ✓
Time range
EPS
Det.
CLC 0–78, 1 h ✓ ✓ m
81–120, 3 h ✓ ✓ m
DTKE_CON 0–48, 1 h ✓ ✓ m15-ke1
0–36, 1 h ✓ ✓ m_6
DTKE_HSH 0–48, 1 h ✓ ✓ m15-ke1
0–36, 1 h ✓ ✓ m_6
P 0–51, 1 h ✓ ✓ m
54–120, 6 h ✓ ✓ m
0–78, 1 h ✓ ✓ m
81–120, 3 h ✓ ✓ m
0–51, 1 h ✓ ✓ m_4
54–78, 3 h ✓ ✓ m_4
Continued on next page
93
Chapter 9. EU Nest output fields
78–120, 6 h ✓ ✓ m_4
QC 0–51, 1 h ✓ ✓ m
54–120, 6 h ✓ ✓ m
0–78, 1 h ✓ ✓ m
81–120, 3 h ✓ ✓ m
0–51, 1 h ✓ ✓ m_4
54–78, 3 h ✓ ✓ m_4
78–120, 6 h ✓ ✓ m_4
QI 0–51, 1 h ✓ ✓ m
54–120, 6 h ✓ ✓ m
0–78, 1 h ✓ ✓ m
81–120, 3 h ✓ ✓ m
0–51, 1 h ✓ ✓ m_4
51–120, 6 h ✓ ✓ m_4
QR 0–51, 1 h ✓ ✓ m
54–120, 6 h ✓ ✓ m
0–51, 1 h ✓ ✓ m_4
QS 0–51, 1 h ✓ ✓ m
54–120, 6 h ✓ ✓ m
0–51, 1 h ✓ ✓ m_4
QV 0–51, 1 h ✓ ✓ m
54–120, 6 h ✓ ✓ m
0–78, 1 h ✓ ✓ m
81–120, 3 h ✓ ✓ m
0–51, 1 h ✓ ✓ m_4
54–78, 3 h ✓ ✓ m_4
78–120, 6 h ✓ ✓ m_4
T 0–51, 1 h ✓ ✓ m
54–120, 6 h ✓ ✓ m
0–78, 1 h ✓ ✓ m
81–120, 3 h ✓ ✓ m
0–51, 1 h ✓ ✓ m_4
54–78, 3 h ✓ ✓ m_4
78–120, 6 h ✓ ✓ m_4
Continued on next page
94
9.2. Multi-level fields on native hybrid vertical levels
Table 9.5.: Multi-level fields on native hybrid vertical levels in forecasts employing prognostic mineral dust
native
latlon
Time range
EPS
Det.
AOD_DUST 0–48, 6 h ✓ ✓ m
60–120, 12 h ✓ ✓ m
0–48, 6 h ✓ ✓ m
60–120, 12 h ✓ ✓ m
CEIL_BSC_DUST
Continued on next page
95
Chapter 9. EU Nest output fields
96
9.2. Multi-level fields on native hybrid vertical levels
60–120, 12 h ✓ ✓ m
0–48, 6 h ✓ ✓ m
60–120, 12 h ✓ ✓ m
DUST_TOTAL_MC 0–120, 3 h ✓ ✓ m
0–120, 3 h ✓ ✓ m
P 0–51, 1 h ✓ ✓ m
54–120, 6 h ✓ ✓ m
0–51, 1 h ✓ ✓ m_4
54–78, 3 h ✓ ✓ m_4
78–120, 6 h ✓ ✓ m_4
QC 0–51, 1 h ✓ ✓ m
54–120, 6 h ✓ ✓ m
0–51, 1 h ✓ ✓ m_4
54–78, 3 h ✓ ✓ m_4
78–120, 6 h ✓ ✓ m_4
QI 0–51, 1 h ✓ ✓ m
54–120, 6 h ✓ ✓ m
0–51, 1 h ✓ ✓ m_4
51–120, 6 h ✓ ✓ m_4
QR 0–51, 1 h ✓ ✓ m
54–120, 6 h ✓ ✓ m
0–51, 1 h ✓ ✓ m_4
QS 0–51, 1 h ✓ ✓ m
54–120, 6 h ✓ ✓ m
0–51, 1 h ✓ ✓ m_4
QV 0–51, 1 h ✓ ✓ m
54–120, 6 h ✓ ✓ m
0–51, 1 h ✓ ✓ m_4
54–78, 3 h ✓ ✓ m_4
78–120, 6 h ✓ ✓ m_4
SAT_BSC_DUST
(1064 nm) 0–48, 6 h ✓ ✓ m
60–120, 12 h ✓ ✓ m
0–48, 6 h ✓ ✓ m
Continued on next page
97
Chapter 9. EU Nest output fields
60–120, 12 h ✓ ✓ m
SAT_BSC_DUST
(532 nm) 0–48, 6 h ✓ ✓ m
60–120, 12 h ✓ ✓ m
0–48, 6 h ✓ ✓ m
60–120, 12 h ✓ ✓ m
T 0–51, 1 h ✓ ✓ m
54–120, 6 h ✓ ✓ m
0–51, 1 h ✓ ✓ m_4
54–78, 3 h ✓ ✓ m_4
78–120, 6 h ✓ ✓ m_4
TKE 0–48, 1 h ✓ ✓ m15-ke1
0–36, 1 h ✓ ✓ m_6
U 0–51, 1 h ✓ ✓ m
54–120, 6 h ✓ ✓ m
0–51, 1 h ✓ ✓ m_4
54–120, 6 h ✓ ✓ m_4
V 0–51, 1 h ✓ ✓ m
54–120, 6 h ✓ ✓ m
0–51, 1 h ✓ ✓ m_4
54–120, 6 h ✓ ✓ m_4
W 0–51, 1 h ✓ ✓ m
54–120, 6 h ✓ ✓ m
0–3, 3 h ✓ ✓ m_5
98
9.3. Multi-level fields interpolated to pressure levels
native
latlon
Time range
EPS
Det.
ShortName level type
CLC 0–78, 1 h ✓ ✓ p5
81–120, 3 h ✓ ✓ p5
FI 0–78, 1 h ✓ ✓ p5
81–120, 3 h ✓ ✓ p5
0–75, 1 h ✓ ✓ pe5
75–120, 3 h ✓ ✓ pe5
OMEGA 0–78, 1 h ✓ ✓ p5
81–120, 3 h ✓ ✓ p5
0–75, 1 h ✓ ✓ pe4
75–120, 3 h ✓ ✓ pe4
RELHUM 0–78, 1 h ✓ ✓ p5
81–120, 3 h ✓ ✓ p5
0–75, 1 h ✓ ✓ pe5
75–120, 3 h ✓ ✓ pe5
T 0–78, 1 h ✓ ✓ p5
81–120, 3 h ✓ ✓ p5
0–75, 1 h ✓ ✓ pe5
75–120, 3 h ✓ ✓ pe5
U 0–78, 1 h ✓ ✓ p5
81–120, 3 h ✓ ✓ p5
0–75, 1 h ✓ ✓ pe5
75–120, 3 h ✓ ✓ pe5
V 0–78, 1 h ✓ ✓ p5
81–120, 3 h ✓ ✓ p5
0–75, 1 h ✓ ✓ pe5
75–120, 3 h ✓ ✓ pe5
99
Chapter 9. EU Nest output fields
native
latlon
Time range
EPS
Det.
ShortName level type
DUST_MAX_TOTAL
_MC_LAYER 0–120, 3 h ✓ ✓ pd7
FI 0–75, 1 h ✓ ✓ pe5
75–120, 3 h ✓ ✓ pe5
OMEGA 0–75, 1 h ✓ ✓ pe4
75–120, 3 h ✓ ✓ pe4
RELHUM 0–75, 1 h ✓ ✓ pe5
75–120, 3 h ✓ ✓ pe5
T 0–75, 1 h ✓ ✓ pe5
75–120, 3 h ✓ ✓ pe5
U 0–75, 1 h ✓ ✓ pe5
75–120, 3 h ✓ ✓ pe5
V 0–75, 1 h ✓ ✓ pe5
75–120, 3 h ✓ ✓ pe5
Time range
EPS
Det.
ShortName
ALB_RAD 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
ALHFL_BS 2–24, 1 h ✓ ✓
ALHFL_PL 2–24, 1 h ✓ ✓
Continued on next page
100
9.4. Single-level fields
ALHFL_S 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
APAB_S 0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
ASHFL_S 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
ASOB_S 0–51, 1 h ✓ ✓
54–120, 3 h ✓ ✓
0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
ASOB_S_CS 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
ASOB_T 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
ASWDIFD_S 0–51, 1 h ✓ ✓
54–120, 3 h ✓ ✓
0–120, 1 h ✓ ✓
0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
ASWDIFU_S 0–51, 1 h ✓ ✓
Continued on next page
101
Chapter 9. EU Nest output fields
54–120, 6 h ✓ ✓
0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
6–120, 6 h ✓ ✓
ASWDIR_S 0–51, 1 h ✓ ✓
54–120, 3 h ✓ ✓
0–120, 1 h ✓ ✓
0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
ATHB_S 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
ATHB_T 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
AUMFL_S 0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
AVMFL_S 0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
CAPE_CON 0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
CAPE_ML 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
CEILING 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–78, 1 h ✓ ✓
Continued on next page
102
9.4. Single-level fields
81–120, 3 h ✓ ✓
CIN_ML 0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
CLCH 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
CLCL 0–51, 1 h ✓ ✓
54–120, 3 h ✓ ✓
0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
CLCM 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
CLCT 0–51, 1 h ✓ ✓
54–120, 3 h ✓ ✓
0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
CLCT_MOD 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
CLDEPTH 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–78, 1 h ✓ ✓
Continued on next page
103
Chapter 9. EU Nest output fields
81–120, 3 h ✓ ✓
C_T_LK 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
FRESHSNW 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–3, 3 h ✓ ✓
FR_ICE 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–120, 6 h ✓ ✓
HBAS_CON 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
HTOP_CON 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
HTOP_DC 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
HZEROCL 0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
HZERO_CL 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
H_ICE 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
0–6, 3 h ✓ ✓
Continued on next page
104
9.4. Single-level fields
6–120, 6 h ✓ ✓
H_ML_LK 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
H_SNOW 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
0–120, 6 h ✓ ✓
LPI_CON_MAX 0–51, 1 h ✓ ✓
54–120, 3 h ✓ ✓
0–48, 1 h ✓ ✓
51–72, 3 h ✓ ✓
78–120, 6 h ✓ ✓
0–48, 1 h ✓ ✓
51–72, 3 h ✓ ✓
78–120, 6 h ✓ ✓
PMSL 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
PS 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
QV_2M 0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
QV_S 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
Continued on next page
105
Chapter 9. EU Nest output fields
0–51, 1 h ✓ ✓
RAIN_CON 0–51, 1 h ✓ ✓
54–120, 3 h ✓ ✓
0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
RAIN_GSP 0–51, 1 h ✓ ✓
54–120, 3 h ✓ ✓
0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
RELHUM_2M 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
0–120, 6 h ✓ ✓
RHO_SNOW 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
RSTOM 2–24, 1 h ✓ ✓
RUNOFF_G 0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
RUNOFF_S 0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
SNOWLMT 0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
SNOW_CON 0–51, 1 h ✓ ✓
54–120, 3 h ✓ ✓
0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
0–75, 1 h ✓ ✓
Continued on next page
106
9.4. Single-level fields
75–120, 3 h ✓ ✓
SNOW_GSP 0–51, 1 h ✓ ✓
54–120, 3 h ✓ ✓
0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
SOBS_RAD 0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
SYNMSG_BT_CL_IR10.8 0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
SYNMSG_BT_CL_WV6.2 0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
TCH 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
TCM 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
TD_2M 0–51, 1 h ✓ ✓
54–120, 3 h ✓ ✓
0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
THBS_RAD 0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
TMAX_2M 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
6–120, 6 h ✓ ✓
Continued on next page
107
Chapter 9. EU Nest output fields
TMIN_2M 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
6–120, 6 h ✓ ✓
TOT_PREC 0–51, 1 h ✓ ✓
54–120, 3 h ✓ ✓
0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
TQC 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
TQI 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
TQR 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
TQS 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
TQV 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–78, 1 h ✓ ✓
Continued on next page
108
9.4. Single-level fields
81–120, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
T_2M 0–51, 1 h ✓ ✓
54–120, 3 h ✓ ✓
0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
T_BOT_LK 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
T_G 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
T_ICE 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
0–3, 3 h ✓ ✓
T_MNW_LK 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
T_SNOW 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
0–51, 1 h ✓ ✓
T_SO 0–51, 1 h ✓ ✓
T_WML_LK 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
U_10M 0–51, 1 h ✓ ✓
54–120, 3 h ✓ ✓
Continued on next page
109
Chapter 9. EU Nest output fields
0–120, 1 h ✓ ✓
0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
VIS 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
VMAX_10M 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–120, 1 h ✓ ✓
1–120, 1 h ✓ ✓
V_10M 0–51, 1 h ✓ ✓
54–120, 3 h ✓ ✓
0–120, 1 h ✓ ✓
0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
WW 0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
W_I 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–3, 3 h ✓ ✓
W_SNOW 0–51, 1 h ✓ ✓
54–120, 3 h ✓ ✓
0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
W_SO 0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
W_SO_ICE 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
Z0 0–78, 1 h ✓ ✓
81–120, 3 h ✓ ✓
0–75, 1 h ✓ ✓
Continued on next page
110
9.4. Single-level fields
75–120, 3 h ✓ ✓
native
latlon
Time range
EPS
Det.
ShortName
ACCDRYDEPO_DUSTA 0–120, 3 h ✓ ✓
0–120, 3 h ✓ ✓
ACCDRYDEPO_DUSTB 0–120, 3 h ✓ ✓
0–120, 3 h ✓ ✓
ACCDRYDEPO_DUSTC 0–120, 3 h ✓ ✓
0–120, 3 h ✓ ✓
ACCEMISS_DUSTA 0–120, 3 h ✓ ✓
0–120, 3 h ✓ ✓
ACCEMISS_DUSTB 0–120, 3 h ✓ ✓
0–120, 3 h ✓ ✓
ACCEMISS_DUSTC 0–120, 3 h ✓ ✓
0–120, 3 h ✓ ✓
ACCSEDIM_DUSTA 0–120, 3 h ✓ ✓
0–120, 3 h ✓ ✓
ACCSEDIM_DUSTB 0–120, 3 h ✓ ✓
0–120, 3 h ✓ ✓
ACCSEDIM_DUSTC 0–120, 3 h ✓ ✓
0–120, 3 h ✓ ✓
ACCWETDEPO
_CON_DUSTA 0–120, 3 h ✓ ✓
0–120, 3 h ✓ ✓
ACCWETDEPO
_CON_DUSTB 0–120, 3 h ✓ ✓
0–120, 3 h ✓ ✓
ACCWETDEPO
Continued on next page
111
Chapter 9. EU Nest output fields
_CON_DUSTC 0–120, 3 h ✓ ✓
0–120, 3 h ✓ ✓
ACCWETDEPO
_GSP_DUSTA 0–120, 3 h ✓ ✓
0–120, 3 h ✓ ✓
ACCWETDEPO
_GSP_DUSTB 0–120, 3 h ✓ ✓
0–120, 3 h ✓ ✓
ACCWETDEPO
_GSP_DUSTC 0–120, 3 h ✓ ✓
0–120, 3 h ✓ ✓
ALB_RAD 0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
ALHFL_BS 2–24, 1 h ✓ ✓
ALHFL_PL 2–24, 1 h ✓ ✓
ALHFL_S 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
ASHFL_S 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
ASOB_S 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
ASOB_S_CS 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
ASOB_T 0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
ASWDIFD_S 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
ASWDIFU_S 6–120, 6 h ✓ ✓
Continued on next page
112
9.4. Single-level fields
ASWDIR_S 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
ATHB_S 0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
ATHB_T 0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
CEILING 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
CLCH 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
CLCL 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
CLCM 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
CLCT 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
C_T_LK 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
DUST_TOTAL_MC_VI 0–120, 3 h ✓ ✓
0–120, 3 h ✓ ✓
FRESHSNW 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–3, 3 h ✓ ✓
FR_ICE 0–51, 1 h ✓ ✓
Continued on next page
113
Chapter 9. EU Nest output fields
54–120, 6 h ✓ ✓
0–120, 6 h ✓ ✓
HBAS_CON 0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
HTOP_CON 0–48, 1 h ✓ ✓
0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
H_ICE 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–6, 3 h ✓ ✓
6–120, 6 h ✓ ✓
H_ML_LK 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
H_SNOW 0–120, 6 h ✓ ✓
LPI_CON_MA 0–48, 1 h ✓ ✓
LPI_CON_MAX 0–48, 1 h ✓ ✓
51–72, 3 h ✓ ✓
78–120, 6 h ✓ ✓
51–72, 3 h ✓ ✓
78–120, 6 h ✓ ✓
PMSL 0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
PS 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
QV_S 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–51, 1 h ✓ ✓
RAIN_CON 0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
RAIN_GSP 0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
RELHUM_2M 0–120, 6 h ✓ ✓
Continued on next page
114
9.4. Single-level fields
RHO_SNOW 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
RSTOM 2–24, 1 h ✓ ✓
SNOW_CON 0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
SNOW_GSP 0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
SOBS_RAD 0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
TAOD_DUST 0–120, 3 h ✓ ✓
0–120, 3 h ✓ ✓
TCH 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
TCM 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
TD_2M 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
THBS_RAD 0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
TMAX_2M 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
6–120, 6 h ✓ ✓
TMIN_2M 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
6–120, 6 h ✓ ✓
TOT_PREC 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
TQC 0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
TQI 0–75, 1 h ✓ ✓
Continued on next page
115
Chapter 9. EU Nest output fields
75–120, 3 h ✓ ✓
TQV 0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
T_2M 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
T_BOT_LK 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
T_G 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
T_ICE 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–3, 3 h ✓ ✓
T_MNW_LK 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
T_SNOW 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–51, 1 h ✓ ✓
T_SO 0–51, 1 h ✓ ✓
T_WML_LK 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
USTAR 0–120, 3 h ✓ ✓
0–120, 3 h ✓ ✓
USTAR_THRES 0–120, 3 h ✓ ✓
0–120, 3 h ✓ ✓
U_10M 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
VMAX_10M 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
Continued on next page
116
9.5. Soil-specific multi-level fields
1–120, 1 h ✓ ✓
V_10M 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
W_I 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–3, 3 h ✓ ✓
W_SNOW 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
W_SO 0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
W_SO_ICE 0–51, 1 h ✓ ✓
54–120, 6 h ✓ ✓
Z0 0–75, 1 h ✓ ✓
75–120, 3 h ✓ ✓
Time range
EPS
Det.
117
Chapter 9. EU Nest output fields
54–120, 6 h ✓ ✓ soil
0–78, 1 h ✓ ✓ soil
81–120, 3 h ✓ ✓ soil
W_SO_ICE 0–78, 1 h ✓ ✓ soil
81–120, 3 h ✓ ✓ soil
118
10. ICON-D2 output fields
This section contains a list of output fields that are available with the launch of ICON-D2. See Fig. 2.8
for details regarding its location and extent. Forecasts of ICON-D2 are performed 8 times a day for the
forecast times 00, 03, 06, 09, 12, 15, 18, 21 UTC, with a forecast range of 48h. Prior to 2021-06-23, the
forecast range was limited to 27h with the exception of the 03 UTC run which reached 45h. During
the pre-operational phase (i.e. prior to 2021-02-10) forecasts of ICON-D2 have been performed 2 times
a day for the forecast times 00 and 12 UTC, with a forecast range of 27h.
Output is available on two distinct horizontal grids:
• a native triangular grid with an average resolution of about 2.1 km, and
• a rotated (!) latitude-longitude grid with a resolution of ∆λ = ∆Φ = 0.02◦ (remark: this horizontal
lat-lon grid is exactly the same as for the former COSMO-D2!). See Table 10.1 for a summary.
The geographical coordinates of every rotated grid point can be found in the fields RLON and RLAT.
This information should be sufficient for the most users (otherwise some more details can be found
in appendix B).
Note that there are a few differences to some of the former COSMO-D2 fields:
– although the velocity components u and v are given on the rotated lat-lon grid points,
too, their components now are the purely (i.e. unrotated) zonal and meridional components,
respectively (in the Grib-Metadata ResolutionAndComponentFlags the 5th bit is 0, whereas
in the former COSMO it was 1).
– Now every variable is interpolated to the same cell center point (whereas in the former
COSMO again the velocity components have been staggered by the half grid mesh size) (in
the Grib-Metadata scanningMode the last four bits are all zero, i.e. no staggering).
– The vertical model levels for 3D fields are slightly different to COSMO. In any case the height
values are given by the vertical averaging of the two neighbouring HHL-values (HHL-fields are
delivered on the rotated lat-lon grid, too).
The output of the most variables takes place hourly. A few variables, which are of particular interest
in the cases of deep convection, are delivered every 15 min.
The model area of ICON-D2 (Fig. 2.8) completely contains the areas of Germany, Switzerland and
Austria and also parts of the neighbouring coutries.
The rotated latitude-longitude output grid contains 651 × 716 = 466116 grid points with a grid mesh
size of 0.02◦ (∼ 2,2 km).
Table 10.1.: Rotated coordinates (λ, φ) and geographical coordinates (λg , φg ) of the four corner points of
the lat-lon grid.
In the subsequent tables the availability of specific fields on the native grid, on the lat-lon grid, or on
both grids is denoted.
119
Chapter 10. ICON-D2 output fields
native
latlon
Time range
EPS
Det.
ShortName level type
CLAT t=0 ✓ ✓
t=0 ✓ ✓
CLON t=0 ✓ ✓
t=0 ✓ ✓
DEPTH_LK t=0 ✓ ✓ ✓
t=0 ✓ ✓
ELAT t=0 ✓ ✓
t=0 ✓ ✓
ELON t=0 ✓ ✓
t=0 ✓ ✓
FR_ICE t=0 ✓ ✓ ✓
t=0 ✓ ✓
FR_LAKE t=0 ✓ ✓ ✓
t=0 ✓ ✓
FR_LAND t=0 ✓ ✓ ✓
t=0 ✓ ✓
HHL t=0 ✓ ✓ ✓ m
t=0 ✓ ✓ m
HSURF t=0 ✓ ✓ ✓
t=0 ✓ ✓
LAI t=0 ✓ ✓ ✓
t=0 ✓ ✓
PLCOV t=0 ✓ ✓ ✓
t=0 ✓ ✓
RLAT t=0 ✓ ✓
RLON t=0 ✓ ✓
ROOTDP t=0 ✓ ✓ ✓
t=0 ✓ ✓
SOILTYP t=0 ✓ ✓ ✓
t=0 ✓ ✓
Continued on next page
120
10.2. Multi-level fields on native hybrid vertical levels
native
latlon
Time range
EPS
Det.
ShortName level type
CLC 0–48, 1 h ✓ ✓ ✓ m
0–48, 1 h ✓ ✓ m
P 0–48, 1 h ✓ ✓ ✓ m
0–48, 1 h ✓ ✓ m
QC 0–48, 1 h ✓ ✓ ✓ m
0–48, 1 h ✓ ✓ m
QG 0–48, 1 h ✓ ✓ ✓ m
0–48, 1 h ✓ ✓ m
QI 0–48, 1 h ✓ ✓ ✓ m
0–48, 1 h ✓ ✓ m
QR 0–48, 1 h ✓ ✓ ✓ m
QS 0–48, 1 h ✓ ✓ ✓ m
QV 0–48, 1 h ✓ ✓ ✓ m
0–48, 1 h ✓ ✓ m
Q_SEDIM 0–48, 1 h ✓ ✓ ✓ m
0–48, 1 h ✓ ✓ m
T 0–48, 1 h ✓ ✓ ✓ m
0–48, 1 h ✓ ✓ m
TKE 0–48, 1 h ✓ ✓ ✓ m
U 0–48, 1 h ✓ ✓ ✓ m
0–48, 1 h ✓ ✓ m
V 0–48, 1 h ✓ ✓ ✓ m
0–48, 1 h ✓ ✓ m
W 0–48, 1 h ✓ ✓ ✓ m
0–48, 1 h ✓ ✓ m
Continued on next page
121
Chapter 10. ICON-D2 output fields
native
latlon
Time range
EPS
Det.
ShortName level type
FI 0–48, 1 h ✓ ✓ ✓ p1
0–48, 1 h ✓ ✓ p1
OMEGA 0–48, 1 h ✓ ✓ ✓ p1
0–48, 1 h ✓ ✓ p1
RELHUM 0–48, 1 h ✓ ✓ ✓ p1
0–48, 1 h ✓ ✓ p1
T 0–48, 1 h ✓ ✓ ✓ p1
0–48, 1 h ✓ ✓ p1
U 0–48, 1 h ✓ ✓ ✓ p1
0–48, 1 h ✓ ✓ p1
V 0–48, 1 h ✓ ✓ ✓ p1
0–48, 1 h ✓ ✓ p1
Time range
EPS
Det.
ShortName
ALB_RAD 0–48, 1 h ✓ ✓ ✓
0–48, 1 h ✓ ✓
ALHFL_S 0–48, 1 h ✓ ✓ ✓
0–48, 1 h ✓ ✓
Continued on next page
122
10.4. Single-level fields
APAB_S 0–48, 1 h ✓ ✓ ✓
ASHFL_S 0–48, 1 h ✓ ✓ ✓
0–48, 1 h ✓ ✓
ASOB_S 0–48, 1 h ✓ ✓ ✓
0–48, 1 h ✓ ✓
ASOB_T 0–48, 1 h ✓ ✓ ✓
0–48, 1 h ✓ ✓
ASWDIFD_S 0–48, 0.25 h ✓ ✓ ✓
0–48, 1 h ✓ ✓
ASWDIFU_S 0–48, 1 h ✓ ✓ ✓
0–48, 1 h ✓ ✓
ASWDIR_S 0–48, 0.25 h ✓ ✓ ✓
0–48, 1 h ✓ ✓
ATHB_S 0–48, 1 h ✓ ✓ ✓
0–48, 1 h ✓ ✓
ATHB_T 0–48, 1 h ✓ ✓ ✓
0–48, 1 h ✓ ✓
AUMFL_S 0–48, 1 h ✓ ✓ ✓
AVMFL_S 0–48, 1 h ✓ ✓ ✓
CAPE_CON 0–48, 1 h ✓ ✓
CAPE_ML 0–48, 0.25 h ✓ ✓ ✓
0–48, 1 h ✓ ✓
CEILING 0–48, 1 h ✓ ✓ ✓
0–48, 1 h ✓ ✓
CIN_ML 0–48, 0.25 h ✓ ✓ ✓
0–48, 1 h ✓ ✓
CLCH 0–48, 1 h ✓ ✓ ✓
0–48, 1 h ✓ ✓
CLCL 0–48, 1 h ✓ ✓ ✓
0–48, 1 h ✓ ✓
CLCM 0–48, 1 h ✓ ✓ ✓
0–48, 1 h ✓ ✓
CLCT 0–48, 1 h ✓ ✓ ✓
0–48, 1 h ✓ ✓
Continued on next page
123
Chapter 10. ICON-D2 output fields
CLCT_MOD 0–48, 1 h ✓ ✓ ✓
0–48, 1 h ✓ ✓
CLDEPTH 0–48, 1 h ✓ ✓ ✓
C_T_LK 0–48, 1 h ✓ ✓
DBZ_850 0–48, 0.25 h ✓ ✓
0–48, 0.25 h ✓ ✓
DBZ_CMAX 0–48, 0.25 h ✓ ✓
0–48, 0.25 h ✓ ✓
0–48, 1 h ✓ ✓
DBZ_CTMAX 0–48, 1 h ✓ ✓
0–48, 1 h ✓ ✓
ECHOTOP 0–48, 0.25 h ✓ ✓
0–48, 0.25 h ✓ ✓
FRESHSNW 0–48, 1 h ✓ ✓
0–48, 1 h ✓ ✓
GRAU_GSP 0–48, 0.25 h ✓ ✓ ✓
0–48, 1 h ✓ ✓
HBAS_SC 0–48, 0.25 h ✓ ✓ ✓
0–48, 1 h ✓ ✓
HTOP_DC 0–48, 1 h ✓ ✓ ✓
HTOP_SC 0–48, 0.25 h ✓ ✓ ✓
0–48, 1 h ✓ ✓
HZEROCL 0–48, 0.25 h ✓ ✓ ✓
H_ICE 0–48, 1 h ✓ ✓ ✓
H_ML_LK 0–48, 1 h ✓ ✓
H_SNOW 0–48, 1 h ✓ ✓ ✓
0–48, 1 h ✓ ✓
LPI 0–48, 0.25 h ✓ ✓ ✓
0–48, 0.25 h ✓ ✓
0–48, 1 h ✓ ✓
LPI_MAX 0–48, 1 h ✓ ✓ ✓
0–48, 1 h ✓ ✓
PMSL 0–48, 1 h ✓ ✓ ✓
0–48, 1 h ✓ ✓
Continued on next page
124
10.4. Single-level fields
125
Chapter 10. ICON-D2 output fields
TCOND_MAX 0–48, 1 h ✓ ✓ ✓
0–48, 1 h ✓ ✓
TD_2M 0–48, 1 h ✓ ✓ ✓
0–48, 1 h ✓ ✓
TMAX_2M 0–48, 1 h ✓ ✓ ✓
0–48, 1 h ✓ ✓
TMIN_2M 0–48, 1 h ✓ ✓ ✓
0–48, 1 h ✓ ✓
TOT_PREC 0–48, 0.25 h ✓ ✓ ✓
0–48, 0.25 h ✓ ✓ ✓
0–48, 1 h ✓ ✓
TQC 0–48, 0.25 h ✓ ✓ ✓
0–48, 1 h ✓ ✓
TQC_DIA 0–48, 1 h ✓ ✓
0–48, 1 h ✓ ✓
TQG 0–48, 0.25 h ✓ ✓ ✓
0–48, 1 h ✓ ✓
TQI 0–48, 0.25 h ✓ ✓ ✓
0–48, 1 h ✓ ✓
TQI_DIA 0–48, 1 h ✓ ✓
0–48, 1 h ✓ ✓
TQR 0–48, 0.25 h ✓ ✓ ✓
TQS 0–48, 0.25 h ✓ ✓ ✓
TQV 0–48, 0.25 h ✓ ✓ ✓
0–48, 1 h ✓ ✓
TQV_DIA 0–48, 1 h ✓ ✓
0–48, 1 h ✓ ✓
TWATER 0–48, 1 h ✓ ✓ ✓
0–48, 1 h ✓ ✓
T_2M 0–48, 1 h ✓ ✓ ✓
0–48, 1 h ✓ ✓
T_BOT_LK 0–48, 1 h ✓ ✓
T_G 0–48, 1 h ✓ ✓ ✓
0–48, 1 h ✓ ✓
Continued on next page
126
10.5. Soil-specific multi-level fields
T_ICE 0–48, 1 h ✓ ✓ ✓
T_MNW_LK 0–48, 1 h ✓ ✓
T_SNOW 0–48, 1 h ✓ ✓ ✓
0–48, 1 h ✓ ✓
T_WML_LK 0–48, 1 h ✓ ✓
UH_MAX 0–48, 1 h ✓ ✓ ✓
0–48, 1 h ✓ ✓
UH_MAX_LOW 0–48, 1 h ✓ ✓ ✓
UH_MAX_MED 0–48, 1 h ✓ ✓ ✓
U_10M 0–48, 1 h ✓ ✓ ✓
0–48, 1 h ✓ ✓
VIS 0–48, 1 h ✓ ✓ ✓
0–48, 1 h ✓ ✓
VMAX_10M 0–48, 1 h ✓ ✓ ✓
0–48, 1 h ✓ ✓ ✓
VORW_CTMAX 0–48, 1 h ✓ ✓ ✓
0–48, 1 h ✓ ✓
V_10M 0–48, 1 h ✓ ✓ ✓
0–48, 1 h ✓ ✓
WW 0–48, 1 h ✓ ✓ ✓
W_CTMAX 0–48, 1 h ✓ ✓ ✓
0–48, 1 h ✓ ✓
W_I 0–48, 1 h ✓ ✓ ✓
0–48, 1 h ✓ ✓
W_SNOW 0–48, 1 h ✓ ✓ ✓
0–48, 1 h ✓ ✓
Z0 0–48, 1 h ✓ ✓ ✓
127
Chapter 10. ICON-D2 output fields
native
latlon
Time range
EPS
Det.
ShortName level type
128
11. Ensemble forecasts with ICON (global,
nested, and limited area mode)
Two ensemble systems based on ICON are running under operational conditions at DWD. The global
ICON ensemble suite started operational in January 2018 now providing short to medium range forecasts
at approx. 26, 5 km (R3B06) horizontal resolution on the global scale with a 13, 2 km (R3B07) nesting
area over Europe. The number of vertical model levels in the actual version is 120. Before 2022-11-23,
the ensemble used a grid resolution of 40 km (R2B06) for the global domain and 20 km (R2B07) for the
nest, both on 90 levels.
The ICON-EPS with its EU-nest runs 8 times a day providing boundary conditions for the ICON-
D2-EPS. At 00/06/12/18 UTC, the whole system including the EU-nest is integrated up to 120h. In
addition, at 00/12 UTC, the global system (without nest) is further integrated to 180 h lead time. At
03/09/15/21 UTC the forecast lead time is limited to +51 h for both, the global domain and the EU
nest.
With the operational start of ICON-D2 (2021-02-10), a convection-permitting ensemble ICON-D2-EPS
(analogous to the former COSMO-D2-EPS) with 20 members and the same resolution and grid as the
deterministic ICON-D2 (≈ 2 km, R19B07) is performed. It runs 8 times a day with 48 hours of forecasts
for the 00, 06, 09, 12, 15, 18, 21 UTC runs. Prior to 2021-06-23 the forecast range was limited to 27
hours, with the exception of the 03 UTC run (45 hours).
The main purpose of an ensemble system is to estimate forecast uncertainty by running a number
of possible physically consistent scenarios of future development. The different scenarios arise from
uncertainties in initial conditions and model error. For limited area ensembles, an additional source of
forecast uncertainty is the uncertainty in the boundary conditions. In the following sections we explain
the techniques used in the ICON-EPS to simulate the effects of those error sources on the forecast and
describe its output data.
In the context of ensemble forecasting it is important to note that the LETKF establishes a square
root filter with multiple variance inflation techniques (see Anlauf et al., 2017, Freitag and Potthast,
2013, Schraff et al., 2016) The ”Kalman gain” from adding observations reduces the uncertainty in the
analysis and thus the variance in the analysis ensemble. By the time this would lead to underestimation
129
Chapter 11. Ensemble forecasts with ICON (global, nested, and limited area mode)
of the true background error compared to the observation error and the analysis ensemble must be
re-inflated in each analysis step to stabilise the ensemble variance. The analysis increments as well as
the partly random variance inflation techniques introduce imbalances in the initial states of the forecast
ensemble, which are damped using an incremental analysis update scheme (IAU; Section 4.2). All the
modifications from the analysis cycle lead to a new analysis ensemble. The new properties and relative
arrangements of the analysis members determine the spread growth and thus the quality of the forecast
uncertainty estimation.
Table 11.1.: List of parameters which are perturbed in the global/EU ICON-EPS (glo) and/or ICON-D2-
EPS (D2). The perturbation mode is either additive (a) or multiplicative (m), possibly limited
to positive (+) deviations from the default value.
130
11.2. Ensemble Physics Perturbations
131
Chapter 11. Ensemble forecasts with ICON (global, nested, and limited area mode)
For models running in limited area mode, the provision of lateral boundaries is an additional source of
forecast uncertainty.
This is accounted for in ICON-D2-EPS by using as hourly lateral boundary conditions forecast mem-
bers of the ICON-EU-EPS which have started 3 hours before the start of the respective ICON-D2-EPS.
Only the members 1–20 of the 40 members of ICON-EU-EPS are used for driving the 20 members of
ICON-D2-EPS. As for the reduction to KENDA members 1–20 for the perturbation of initial conditions
(Section 11.1), this is statistically equivalent to a random selection of 20 members of ICON-EU-EPS
over a large number of ICON-D2-EPS runs. However, there are promising approaches in the research
community to select a suitable subset of 20 out of 40 members at each forecast following certain op-
timisation criteria. The application and suitability of such methods has not yet been tested for the
combination of ICON-EU-EPS and ICON-D2-EPS.
132
11.4. ICON Ensemble output fields in DWD databases
1 https://www.dwd.de/DE/wetter/warnungen_aktuell/kriterien/warnkriterien.html
2 http://www.wmo.int/pages/prog/www/DPS/Publications/WMO_485_Vol_I.pdf
133
Chapter 11. Ensemble forecasts with ICON (global, nested, and limited area mode)
products are generated every 6 hours up to 120 h lead time on the EU domain and for the global fields
up to 180 h twice a day (00/12 UTC). This is done for different accumulation periods depending on
the forecast parameter. In the following tables, the accumulation time range is given in hours. The
meaning of the level types is:
• (no key): 2D field
• pe4: pressure level 500 hPa
• pe7: pressure level 850 hPa
• pe8: pressure levels 500, 850 hPa
• pe9: pressure levels 250, 500, 850 hPa
All products are delivered only on the lat/lon grid (i.e. not on the ICON native grid).
Table 11.2.: EPS products from ICON global. See table 11.4 for a description of the various product types.
level type
Time range
ShortName Accum. Product
timerange type
134
11.4. ICON Ensemble output fields in DWD databases
135
Chapter 11. Ensemble forecasts with ICON (global, nested, and limited area mode)
Table 11.3.: EPS products from ICON-EU. See table 11.4 for a description of the various product types.
level type
Time range
ShortName Accum. Product
timerange type
136
11.4. ICON Ensemble output fields in DWD databases
Most of the parameters are available on both domains, but there are exceptions: SP250, FI500,T_SO
and the temperature anomaly are available on the global domain only. The latter is calculated for thresh-
olds of ±1, ±1.5 and ±2 standard deviations with respect to the reanalysis climatology ERA_INTERIM3 .
The global products are available via the WMO WIS/WMS system or directly as grib files and charts
on the opendata server of DWD4 in /weather/wmc/icon-eps. A graphical user interface for direct access
to the charts is available on the DWD website5 . The dissemination of the EU-Nest ensemble product
grib files via the opendata server of DWD is planned to start in October 2018.
3 https://www.ecmwf.int/en/forecasts/datasets/archive-datasets/reanalysis-datasets/era-interim
4 https://www.dwd.de/opendata
5 https://www.dwd.de/EN/ourservices/wmc/wmc.html
137
Chapter 11. Ensemble forecasts with ICON (global, nested, and limited area mode)
Table 11.4.: Meaning of the ’product type’ key for EPS-products from ICON (global) and ICON-EU EPS.
perct_1 perct; mean; spread; min; max; 10; 25; 50; 75; 90 (ensemble distribution)
probt_1 1.0; 5.0; 10.0 (precipitation thresholds (in mm))
probt_2 1.0; 5.0; 10.0; 15.0; 20.0; 25.0; 30.0; 50.0 (precipitation thresholds (in mm))
probt_3 1.0; 5.0; 10.0; 20.0; 25.0; 30.0; 50.0; 80.0; 100.0 (precipitation thresholds (in mm))
probt_4 20.0; 35.0; 60.0 (precipitation thresholds (in mm))
probt_5 25.0; 40.0; 70.0 (precipitation thresholds (in mm))
probt_6 40.0; 60.0; 90.0; 150.0 (precipitation thresholds (in mm))
probt_7 50.0; 100.0; 150.0; 250.0 (precipitation thresholds (in mm))
probt_8 11.0; 12.5; 14.0; 18.0; 21.0; 25.0; 29.0; 33.0; 39.0 (wind speed thresholds (in m/s))
probt_9 273.15 (temperature thresholds (in K))
probt_10 253.15; 263.15; 273.15; 293.15 (temperature thresholds (in K))
probt_11 263.15; 273.15; 298.15; 303.15; 308.15 (temperature thresholds (in K))
probt_12 750.0; 1000.0; 1500.0; 2000.0; 3000.0; 4000.0 (CAPE thresholds)
probt_13 50.0; 87.5 (cloud cover thresholds (in %))
probt_14 -1.0; -1.5; -2.0; 1.0; 1.5; 2.0 (temperature anomalie thresholds)
2. Percentiles,
i.e. physical values of a forecast parameter (e.g. T_2M,. . . ), which define the
perc=10,25,30,50,75,90 [%] parts of the ensemble distribution.
3. Exceedance Probabilities
Probability of event above lower limit (probt=3)
Probability of event below upper limit (probt=4)
4. Grid specification
Products on the rotated latlon COSMO-D2 grid (localTypeOfEnsembleProductGeneration=1)
Products on the upscaled (10x10 grid point) grid (localTypeOfEnsembleProductGeneration=101)
The thresholds for the exceedance probabilities are given by the DWD alert thresholds which are used
138
11.4. ICON Ensemble output fields in DWD databases
Table 11.5.: EPS products from ICON-D2. See table 11.6 for a description of the various product types.
level type
Time range
ShortName Accum. Product
timerange type
6 https://www.dwd.de/DE/wetter/warnungen_aktuell/kriterien/warnkriterien.html
7 http://www.wmo.int/pages/prog/www/DPS/Publications/WMO_485_Vol_I.pdf
139
Chapter 11. Ensemble forecasts with ICON (global, nested, and limited area mode)
1–48, 1 h 60 perct_1d
1–48, 1 h 60 probt_1d
12–48, 1 h 720 perct_1d
12–48, 1 h 720 probt_9d
12–48, 1 h 720 probt_9d_ups
1–48, 1 h 60 probt_1d_ups
6–48, 1 h 360 probt_3d_ups
TOT_SNOW 6–48, 1 h 360 perct_1d
6–48, 1 h 360 probt_5d
1–48, 1 h 60 perct_1d
1–48, 1 h 60 probt_4d
12–48, 1 h 720 perct_1d
12–48, 1 h 720 probt_2d
T_2M 1–48, 1 h perct_1d
1–48, 1 h probt_11d
T_G 1–48, 1 h perct_1d
1–48, 1 h probt_10d
VMAX_10M 6–24, 6 h 360 perct_1d
6–24, 6 h 360 probt_8d
1–48, 1 h 60 perct_1d
1–48, 1 h 60 probt_8d
12–24, 12 h 720 perct_1d
12–24, 12 h 720 probt_8d
140
11.4. ICON Ensemble output fields in DWD databases
Table 11.6.: Meaning of the ’product type’ key for EPS-products from ICON-D2-EPS.
perct_1d perct; mean; spread; min; max; 10; 25; 30; 50; 75; 90 (ensemble distribution)
probt_1d_ups as probt_1d, but upscaled to 10x10 gridpoints
probt_1d > 0.1; > 1.0; > 2.0; > 5.0; > 10.0; > 15.0; > 25.0; > 40.0 (probability thresholds)
probt_2d > 0.1; > 1.0; > 5.0; > 10.0; > 15.0; > 20.0; > 25.0; > 30.0; > 50.0 (probability
thresholds)
probt_3d > 0.1; > 1.0; > 2.0; > 5.0; > 10.0; > 20.0; > 35.0; > 60.0 (probability thresholds)
probt_3d_ups as probt_3d, but upscaled to 10x10 gridpoints
probt_4d > 0.1; > 1.0; > 2.0; > 5.0 (probability thresholds)
probt_5d > 0.1; > 5.0; > 10.0 (probability thresholds)
probt_6d < 50.0; > 87.5 (probability thresholds)
probt_7d > 750.0; > 1000.0; > 1500.0; > 2000.0 (probability thresholds)
probt_8d > 11.0; > 12.5; > 14.0; > 18.0; > 21.0; > 25.0; > 29.0; > 33.0; > 39.0
(probability thresholds)
probt_9d > 25.0; > 40.0; > 70.0 (probability thresholds)
probt_9d_ups as probt_9d, but upscaled to 10x10 gridpoints
probt_10d < 273.15 (probability thresholds)
probt_11d < 253.15; < 263.15; < 273.15; >= 293.15; >= 298.15; >= 303.15; >= 308.15
(probability thresholds)
probt_12d < 263.15; < 273.15; >= 298.15; >= 303.15; >= 308.15 (probability thresholds)
probt_13d > 28.0; > 37.0; > 46.0; > 54.0 (probability thresholds)
probt_14d > 0.1; > 1.0; > 5.0; > 10.0; > 20.0 (probability thresholds)
probt_14d_ups as probt_14d, but upscaled to 10x10 gridpoints
141
12. ICON data in the
SKY data bases of DWD
GRIB data of the numerical weather prediction models are stored in the data base SKY at DWD.
Documentation on the SKY system is available in the intranet of DWD at IT/Messnetz/Technik
→ Preprocessing → Management der DWD Fachdaten - Dokumentation → SKY. Below, some re-
marks are given on the SKY categories for ICON data, and some examples are given how to retrieve
data from the data base.
In SKY the data is stored in different categories and data base subsystems. These are identified by the
cat=CAT_NAME parameter. The name of a category is made up of 4 parts:
$model_$run_$type_$suite
run, type, and suite are general for all forecast models of DWD. They can have the following values:
• run: main for main forecast runs, ass for assimilation runs, pre for pre-assimilation runs, const
for invariant data.
• type: an for analysis data, fc for forecast data, fcprod for EPS products.
For forecasts employing prognostic mineral dust, the suffix aero is added, i.e. we have anaero,
fcaero and fcprodaero respectively.
• suite:
– rout for operational data in db=roma,
– para1 for pre-operational data in db=parma, or vera for pre-operational data in db=vera.
– exp or exp1 for data from experiments in db=numex. The category extension exp1 is used
for experiments of the NUMEX wizard, a special NUMEX user.
Data from experiments is additionally identified by the parameter exp=N U M where N U M
is the experiment number.
The model part of the sky-categories for ICON itself is constructed by several substrings. It starts
• with the string ico for ICON data on the native ICON grid, or
• with icr for data on a regular lat-lon or a rotated lat-lon grid.
Next follows a two-letter string to identify the domain of ICON;
• gl for the global domain,
• eu for the nest over Europe,
• la for the limited-area model ICON-D2.
143
Chapter 12. ICON data in the SKY data bases of DWD
Until 2022-11-22 the category names for the deterministic runs of the global domain included the
resolution and the number of levels, i.e. icogl130l90 or icrgl130l90.
For ensemble forecasts or ensemble analyses the first part of the category is extended by an e (for instance
icogle). Except for 5 fields for ICON-D2 there is no output of ensemble forecasts on the regular grid.
Ensemble members or ranges of ensemble members are specified by the parameter enum=N U M or
enum=N U M 1/to/N U M 2 where N U M is the member id. enum must be given to get output from the
ensemble categories. To get all members use enum=1/to/. Ensemble products are available only on a
regular grid.
Hence, the full category name for data from a global operational deterministic forecast run of ICON on
a regular grid is icrgl_main_fc_rout. The initial analysis for this run is in category
icogl_main_an_rout.
Since 2014-08-12 12 UTC ICON is running pre-operationally at DWD. Hence, forecast data
was available in the sky database db=parma in categories icogl130l90_main_fc_para and
icrgl_main_fc_para. Data of the present pre-operational ICON runs is in db=parma in cat-
egories icogl_main_fc_para1 and icrgl_main_fc_para1.
Since 2015-01-20 06 UTC the global ICON model is running operationally at DWD. Forecast
data was available in the sky database db=roma in categories icogl130l90_main_fc_rout
and icrgl130l90_main_fc_rout. Analysis data is available in icogl130l90_ass_an_rout.
Present analysis data is in icogl_ass_an_rout.
Since 2016-01-20 06 UTC an ensemble data assimilation for ICON is running operationally
at DWD. Analysis data is available in the sky database db=roma in the ensemble cat-
egories icogle_main_an_rout and icoeue_main_an_rout. First guess data is in
icogle_pre_fc_rout and icogle_pre_fc_rout. The ensemble runs of ICON write data only
on the native ICON grid. Data on regular grids must be interpolated from the native grid.
Since 2018-01-17 06 UTC the global ICON-EPS is running operationally at DWD. Forecast
data is available in the sky database db=roma in categories icogle_main_fc_rout and
icoeue_main_fc_rout for the EU domain.
Since 2023-11-27 00 UTC the global ICON(-ART) model employing prognostic mineral dust is run-
ning operationally at DWD. This includes an ensemble data assimilation and the deterministic as
well as ensemble forecast runs for a ten member ensemble. Data of the operational routine is saved
in the sky database db=roma.
Global deterministic forecast data is available in the category icogl_main_fcaero_rout. Anal-
ysis data is available in icogl_ass_anaero_rout.
Data for the ICON-EU-NA2 nest is available in the categories icoeu_main_fcaero_rout and
icoeu_ass_anaero_rout respectively.
Forecast data of the ten member ensemble employing prognostic mineral dust is available in the
categories icogle_main_fcaero_rout for the global domain and icoeue_main_fcaero_rout
for the nest domain. Analysis data is avalable in the global icogle_main_anaero_rout
and nest icoeue_main_anaero_rout ensemble categories. Likewise first guess data is in
144
12.2. Retrieving ICON data from SKY
icogle_pre_fcaero_rout and icogle_pre_fcaero_rout. All output for the runs with prog-
nostic mineral dust write data only on the native ICON grid. Data on regular grids must be
interpolated from the native grid.
Data from the operational forecast runs of ICON on the native ICON grid, cat =
ico*_main_fc_rout, is kept in the database roma only for 15 months! Analysis and first
guess data is kept forever.
ICON global
• Retrieve the 2 m temperature and dew point temperature for forecast hours 3 to 78 every 3 hours
of today’s run at 00 UTC on the global domain from an ICON run on a R3B07 grid with 120
levels to file icon2mdat
read db = roma cat = i c o g l _ m a i n _ f c _ r o u t d = t00 s [ h ]=3/ to /78/ by /3 p = t_2m , td_2m bin f =
icon2mdat
• Retrieve the analysis of T on the native grid and model levels (lv=genv) from yesterday 18 UTC:
read db = roma cat = i c o g l _ m a i n _ a n _ r o u t d = t18 -1 d p = T lv = genv gptype =0 bin f = t_icon_ana
• Get the 6, 12, 18, and 24 hour forecast of the 2 m temperature from a forecast in experiment 11503
on 2022-08-29 at 00 UTC from an ICON run on a R3B07 grid with 120 levels. Retrieve data on
the regular lat/lon grid:
read db = numex cat = i c r g l _ m a i n _ f c _ e x p exp =11503 d =2022082900 s [ h ]=6 ,12 ,18 ,24 p = t_2m
bin f = t_2m_fc . grb
• Retrieve the 12 hour forecast on 2024-03-13 at 00 UTC of the column integrated mineral dust
optical depth (for a wavelength of 550 nm) and the total mineral dust mass concentration on
model levels to file icon-art_dust:
read db = roma cat = i c o g l _ m a i n _ f c a e r o _ r o u t d =2024031300 s [ h ]=12 p = taod_dust ,
dust_total_mc bin f = icon - art_dust
ICON-EU
• Retrieve accumulated precipitation of the ICON-EU nest on the regular grid every 6 hours to 72
hours from yesterday’s operational run at 12 UTC:
read db = roma cat = i c r e u _ m a i n _ f c _ r o u t d = t12 -1 d s [ h ]=6/ to /72/ by /6 p = tot_prec bin f =
tot_prec_ieu
145
Chapter 12. ICON data in the SKY data bases of DWD
• List the data on pressure levels of the 18 hours forecast from 06 UTC of ICON-EU nest on the
regular grid. Write reference date (d), forecast step (s), level type (lv), value of first level (lv1),
decoding date (dedat), and store date (stdat) in information file icr.info.
read db = roma cat = i c r e u _ m a i n _ f c _ r o u t d =06 step [ h ]=18 lv = P info = metaData metaArray =d ,
s ,p , lv , lv1 , dedat , stdat sort =p , lv1 infof = icr . info
• Retrieve the 12 hour forecast on 2024-03-13 at 00 UTC for the ICON-EU-NA2 domain of the
attenuated backscatter for mineral dust seen from the ground for a wavelength of 1064 nm on
model levels to file icon-euna2-art_ceil_bsc_dust:
read db = roma cat = i c o e u _ m a i n _ f c a e r o _ r o u t d =2024031300 s [ h ]=12 p = ceil_bsc_dust wvl1
=1.064 E -6 bin f = icon - euna2 - a r t _ c e i l _ b s c _ d u s t
ICON-D2
• Retrieve accumulated precipitation of the operational ICON-D2 on the rotated lat-lon grid every
3 hours to 27 hours from the first day of the pre-operational start at 12 UTC:
read db = roma cat = i c r l a _ m a i n _ f c _ r o m a d =2019112512 s [ h ]=3/ to /27/ by /3 p = tot_prec bin f
= tot_prec_id2
• Retrieve the surface net short-wave radiation flux (averaged since model start) of the operational
ICON-D2 on the native grid every hour to 12 hours for today’s run at 00 UTC:
read db = roma cat = i c o l a _ m a i n _ f c _ r o u t d = t00 s [ h ]=0/ to /12/ by /1 p = asob_s bin f =
asob_s_id2
ICON global
• ICON-EPS: Retrieve the 2 m temperature and dew point temperature for forecast hours 3 to
78 every 3 hours of today’s run at 00 UTC on the global domain from an ICON-EPS run on a
R3B06 grid to file iconEPS2mdat (use cat=icoeue_main_fc_rout for corresponding forecasts of
ICON-EU-EPS on R3B07)
read db = roma cat = i c o g l e _ m a i n _ f c _ r o u t d = t00 s [ h ]=3/ to /78/ by /3 p = t_2m , td_2m bin f =
iconEPS2mdat
• Retrieve temperature in 850 hPa from the first guess of the 40 ensemble members of the EDA on
the 26 km grid in the parallel suite yesterday at 21 UTC. Sort the data by ensemble member.
read db = parma cat = i c o g l e _ a s s _ f c _ p a r a 1 enum =1/ to / d = t21 -1 d s =3 p = T lv = P lv1 =85000
bin f = T850_eps sort = enum
ICON-EU
• ICON-EPS: Retrieve 90% percentile (on regular lat/lon grid) of accumulated precipitation (avail-
able accumulation periods) at forecast hour 72 of today’s run at 06 UTC on the EU domain.
read db = roma cat = i c r e u e _ m a i n _ f c p r o d _ r o u t d = t06 s [ h ]=72 perc =90 p = TOT_PREC bin f =
i co nE P S_ RR 7 2_ 9 0
146
12.2. Retrieving ICON data from SKY
• Retrieve ensemble spread (deriv=4) of CAPE for forecast hours 6 to 120 every 6 hours of today’s
run at 00 UTC on the EU domain (on regular lat/lon grid).
read db = roma cat = i c r e u e _ m a i n _ f c p r o d _ r o u t d = t00 s [ h ]=6/ to /120/ by /6 deriv =4 p = CAPE_ML
bin f = i c o n E P S _ C A P E _ s p r e a d
• Retrieve probabilities of TMIN_2M of the last 12 h and 24 h for any available threshold, where the
probability of event is above lower limit (probt=3), for all available forecast hours of today’s run
at 00 UTC on the EU domain.
read db = roma cat = i c r e u e _ m a i n _ f c p r o d _ r o u t d = t00 probt =3 p = TMIN_2M bin f =
iconEPS_TMIN_2M_probt3
ICON-D2
• Retrieve ensemble spread (deriv=4) of T_2M for all forecast hours (from 0 to 27 every hour) of
a run starting at 00 UTC.
read db = roma cat = i c r l a e _ m a i n _ f c p r o d _ r o m a d =2020050100 s [ h ]=0/ to /27/ by /1 deriv =4 p =
T_2M bin f = i c o l a e _ T 2 M _ s p r e a d
• Retrieve probabilities of TOT_PREC for any available threshold, where the probability of event
is above lower limit (probt=3), for the forecast hours between 0 and 24 with 6-hourly step, all
available accumulations (1h, 6h, 12h), of one run at 00 UTC.
read db = roma cat = i c r l a e _ m a i n _ f c p r o d _ r o m a d =2020050100 s [ h ]=0/ to /24/ by /6 probt =3 p =
TOT_PREC bin f = icolae _TP_pr ob
• Retrieve 90% percentile of CAPE_ML (6h accumulation period) at forecast hour 12 of today’s
run at 00 UTC.
read db = roma cat = i c r l a e _ m a i n _ f c p r o d _ r o m a d = t00 s [ h ]=12 perc =90 p = CAPE_ML bin f =
icolae_CAPE_perc
• The corresponding request to retrieve ICON-D2-EPS products from the parallel suite is
read db = vera cat = i c r l a e _ m a i n _ f c p r o d _ v e r a d = t00 s [ h ]=12 perc =90 p = CAPE_ML bin f =
icolae_CAPE_perc
147
149
zih0 + zi+1
h0
zif 0 = .
2
See Table A.2 for a list of all full level heights of the operational setup for ICON global/ICON-EU and
Table A.4 for ICON-D2.
As an example, Tables A.5 and A.6 show these model heights for a special grid point over India with a
quite high surface elevation.
150
Table A.1.: Standard heights zih0 (i.e. for zero topography height) for all 121 vertical half levels of the
global 13 km domain and the 75 vertical half levels for the 6.5 km EU nest.
Table A.2.: Standard heights zif 0 (i.e. for zero topography height) for all 120 vertical full levels of the global
13 km domain and the 74 full levels of the 6.5 km EU nest.
Table A.3.: Standard heights zih0 (i.e. for zero topography height) for all 66 vertical half levels of ICON-D2.
level idx. height [m] level idx. height [m] level idx. height [m] level idx. height [m]
1 22 000.000 19 8 256.329 37 3 333.549 55 702.132
2 19 401.852 20 7 890.952 38 3 138.402 56 606.827
3 18 013.409 21 7 539.748 39 2 949.656 57 516.885
4 16 906.264 22 7 201.825 40 2 767.143 58 432.419
5 15 958.169 23 6 876.388 41 2 590.708 59 353.586
6 15 118.009 24 6 562.725 42 2 420.213 60 280.598
7 14 358.139 25 6 260.200 43 2 255.527 61 213.746
8 13 661.439 26 5 968.239 44 2 096.537 62 153.438
9 13 016.363 27 5 686.321 45 1 943.136 63 100.277
10 12 414.654 28 5 413.976 46 1 795.234 64 55.212
11 11 850.143 29 5 150.773 47 1 652.748 65 20.000
12 11 318.068 30 4 896.323 48 1 515.610 66 0.000
13 10 814.653 31 4 650.265 49 1 383.761
14 10 336.841 32 4 412.272 50 1 257.155
15 9 882.112 33 4 182.043 51 1 135.760
16 9 448.359 34 3 959.301 52 1 019.556
17 9 033.796 35 3 743.791 53 908.539
18 8 636.893 36 3 535.279 54 802.721
153
Table A.4.: Standard heights zif 0 (i.e. for zero topography height) for all 65 vertical full levels of ICON-D2.
level idx. height [m] level idx. height [m] level idx. height [m] level idx. height [m]
1 20 700.926 18 8 446.611 35 3 639.535 52 964.048
2 18 707.630 19 8 073.640 36 3 434.414 53 855.630
3 17 459.836 20 7 715.350 37 3 235.976 54 752.427
4 16 432.216 21 7 370.787 38 3 044.029 55 654.479
5 15 538.089 22 7 039.106 39 2 858.399 56 561.856
6 14 738.074 23 6 719.557 40 2 678.926 57 474.652
7 14 009.789 24 6 411.462 41 2 505.461 58 393.002
8 13 338.901 25 6 114.219 42 2 337.870 59 317.092
9 12 715.508 26 5 827.280 43 2 176.032 60 247.172
10 12 132.398 27 5 550.148 44 2 019.836 61 183.592
11 11 584.105 28 5 282.374 45 1 869.185 62 126.857
12 11 066.360 29 5 023.548 46 1 723.991 63 77.745
13 10 575.747 30 4 773.294 47 1 584.179 64 37.606
14 10 109.477 31 4 531.269 48 1 449.686 65 10.000
15 9 665.235 32 4 297.157 49 1 320.458
16 9 241.077 33 4 070.672 50 1 196.457
17 8 835.344 34 3 851.546 51 1 077.658
154
Table A.5.: Height above ground zih (x) (half levels) for the grid point with maximum topography height in
the operational setup R03B07, 13 km spatial resolution.
level idx. height [m] level idx. height [m] level idx. height [m] level idx. height [m]
1 68 779.785 32 28 946.129 63 8 390.156 94 2 225.703
2 67 200.387 33 27 979.969 64 8 078.022 95 2 088.177
3 65 649.395 34 27 033.848 65 7 776.084 96 1 950.706
4 64 107.981 35 26 107.772 66 7 483.826 97 1 813.210
5 62 585.699 36 25 201.733 67 7 192.479 98 1 675.732
6 61 082.684 37 24 315.734 68 6 902.146 99 1 539.776
7 59 599.020 38 23 449.777 69 6 612.939 100 1 409.701
8 58 134.805 39 22 603.867 70 6 324.984 101 1 285.344
9 56 690.113 40 21 778.006 71 6 038.416 102 1 166.802
10 55 265.024 41 20 972.199 72 5 753.389 103 1 054.064
11 53 859.598 42 20 186.453 73 5 470.070 104 947.087
12 52 473.891 43 19 420.775 74 5 188.639 105 845.981
13 51 107.957 44 18 675.178 75 4 941.922 106 750.672
14 49 761.840 45 17 949.674 76 4 735.653 107 661.273
15 48 435.574 46 17 243.703 77 4 563.203 108 577.713
16 47 129.195 47 16 550.117 78 4 425.707 109 499.860
17 45 842.738 48 15 876.354 79 4 288.199 110 427.890
18 44 576.219 49 15 215.272 80 4 150.699 111 361.787
19 43 329.668 50 14 574.893 81 4 013.218 112 301.589
20 42 103.098 51 13 955.242 82 3 875.690 113 247.134
21 40 896.524 52 13 356.359 83 3 738.201 114 198.440
22 39 709.957 53 12 778.283 84 3 600.717 115 155.705
23 38 543.410 54 12 221.057 85 3 463.198 116 118.833
24 37 396.891 55 11 688.190 86 3 325.699 117 87.462
25 36 270.406 56 11 180.248 87 3 188.214 118 61.043
26 35 163.957 57 10 699.899 88 3 050.701 119 40.010
27 34 077.547 58 10 246.899 89 2 913.215 120 19.982
28 33 011.180 59 9 819.694 90 2 775.681 121 0.000
29 31 964.852 60 9 449.262 91 2 638.213
30 30 938.570 61 9 074.562 92 2 500.675
31 29 932.328 62 8 721.942 93 2 363.176
155
Table A.6.: Height above ground zif (x) (full levels) for the grid point with maximum topography height in
the operational setup R03B07, 13 km spatial resolution.
level idx. height [m] level idx. height [m] level idx. height [m] level idx. height [m]
1 67 990.086 31 29 439.229 61 8 898.252 91 2 569.444
2 66 424.891 32 28 463.049 62 8 556.049 92 2 431.926
3 64 878.688 33 27 506.909 63 8 234.089 93 2 294.439
4 63 346.840 34 26 570.810 64 7 927.053 94 2 156.940
5 61 834.192 35 25 654.753 65 7 629.955 95 2 019.441
6 60 340.852 36 24 758.734 66 7 338.153 96 1 881.958
7 58 866.912 37 23 882.755 67 7 047.312 97 1 744.471
8 57 412.459 38 23 026.822 68 6 757.542 98 1 607.754
9 55 977.568 39 22 190.936 69 6 468.962 99 1 474.738
10 54 562.311 40 21 375.103 70 6 181.700 100 1 347.523
11 53 166.745 41 20 579.326 71 5 895.903 101 1 226.073
12 51 790.924 42 19 803.614 72 5 611.729 102 1 110.433
13 50 434.898 43 19 047.977 73 5 329.354 103 1 000.576
14 49 098.707 44 18 312.426 74 5 065.280 104 896.534
15 47 782.385 45 17 596.689 75 4 838.788 105 798.327
16 46 485.966 46 16 896.910 76 4 649.428 106 705.973
17 45 209.478 47 16 213.235 77 4 494.455 107 619.493
18 43 952.943 48 15 545.813 78 4 356.953 108 538.786
19 42 716.383 49 14 895.083 79 4 219.449 109 463.875
20 41 499.811 50 14 265.068 80 4 081.958 110 394.838
21 40 303.240 51 13 655.801 81 3 944.454 111 331.688
22 39 126.683 52 13 067.321 82 3 806.945 112 274.361
23 37 970.151 53 12 499.670 83 3 669.459 113 222.787
24 36 833.649 54 11 954.624 84 3 531.957 114 177.072
25 35 717.182 55 11 434.219 85 3 394.448 115 137.269
26 34 620.752 56 10 940.073 86 3 256.957 116 103.148
27 33 544.363 57 10 473.399 87 3 119.457 117 74.252
28 32 488.016 58 10 033.297 88 2 981.958 118 50.526
29 31 451.711 59 9 634.478 89 2 844.448 119 29.996
30 30 435.449 60 9 261.912 90 2 706.947 120 9.991
157
Note, that all angles are given in arcs (not in degrees). To get the angle in degrees, one has to multiply
by 180/π ≈ 57, 2957795. Take care that the arctan is correctly evaluated in all 4 quadrants.2 .
In the dwdlib (in particular in the library libmisc.a) the four Fortran functions RLSTORL, PHSTOPH,
RLTORLS and PHTOPHS are contained, which calculate the transformtions. These programs give and
expect angles in degrees.
1 In former models that used a spherical coordinate system, not only for output but as the base for their numerical grid
(e.g. the COSMO model), the use of a rotated grid was necessary to avoid too narrow grid cells near the poles (the so
called ’pole problem’). Narrow grid cells induce strong time step restrictions and therefore would result in inefficient
code.
2 Most programming languages have an extension of the standard arctan-function, e.g. in Fortran the function
• RLSTORL calculates geographic longitude (RL) from rotated longitude and latitude.
• PHSTOPH calculates geographic latitude (PH) from rotated longitude and latitude.
• RLTORLS calculates rotated longitude (RLS) from geographic longitude and latitude.
• PHTOPHS calculates rotated latitude (PHS) from geographic longitude and latitude.
For the transformation of many points or even whole fields dwdlib also contains the better optimized rou-
tines PLSTOPL, PLTOPLS, APLSTPL, and APLTPLS. An online description can be get via disdwd PLSTOPL
or man libmisc.
158
Bibliography
Anlauf, H., A. Rhodin, A. Cress, R. Faulwetter, M. Lange, G. Paul, C. Kopken-Watts, O. O. Stiller,
and R. Potthast, 2017: Deterministic global data assimilation for the ICON model - system suite
and evaluation. Mon. Wea. Rev., in preparation.
Bloom, S. C., L. L. Takacs, A. M. D. Silva, and D. Ledvina, 1996: Data assimilation using incremental
analysis updates. Mon. Wea. Rev., 124, 1256–1270.
Doms, G., J. Förstner, E. Heise, H.-J. Herzog, D. Mironov, M. Raschendorfer, T. Reinhardt, B. Ritter,
R. Schrodin, J.-P. Schulz, and G. Vogel, 2011: A description of the nonhydrostatic regional COSMO
model. Part II: Physical parameterization. Technical report, Deutscher Wetterdienst. Available at
http://www.cosmo-model.org/.
Freitag, M. and R. Potthast, 2013: Large Scale Inverse Problems Computational Methods and Applica-
tions In the Earth Sciences. De Gruyter.
Hunt, B. R., E. J. Kostelich, and I. Szunyogh, 2007: Efficient data assimilation for spatiotemporal
chaos: A local ensemble transform Kalman filter. Physica D: Nonlinear Phenomena, 230, 112–126.
Leuenberger, D., M. Koller, and C. Schär, 2010: A generalization of the sleve vertical coordinate. Mon.
Wea. Rev., 138, 3683–3689.
Lopez, P., 2016: A Lightning Parameterization for the ECMWF Integrated Forecasting System. Mon.
Wea. Rev., 144(9), 3057–3075.
Lynch, P., 1997: The Dolph-Chebyshev window: A simple optimal filter. Mon. Wea. Rev., 125, 655–660.
Lynn, B. and Y. Yair, 2010: Prediction of lightning flash density with the WRF model. Adv. Geosci.,
23, 11–16.
Polavarapu, S., S. Ren, A. M. Clayton, D. Sankey, and Y. Rochon, 2004: On the relationship between
incremental analysis updating and incremental digital filtering. Mon. Wea. Rev., 132, 2495–2502.
Prill, F., D. Reinert, D. Rieger, and G. Zängl, 2023: ICON Tutorial: Working with the ICON model.
Technical report, Deutscher Wetterdienst.
Ruppert, T., 2007: Diplomarbeit: Vector field reconstruction by radial basis functions. Master’s thesis,
Technical University Darmstadt, Department of Mathematics.
Schraff, C., H. Reich, A. Rhodin, A. Schomburg, K. Stephan, A. Perianez, and R. Potthast, 2016:
Kilometre-scale ensemble data assimilation for the COSMO model (KENDA). Q.J.R. Meteorol. Soc.,
142, 1453–1472.
Stephan, K., S. Klink, and C. Schraff, 2008: Assimilation of radar-derived rain rates into the
convective-scale model COSMO-DE at DWD. Quart. J. Royal Met. Soc., 134(634), 1315–1326.
DOI: 10.1002/qj.269.
159
BIBLIOGRAPHY
Wicker, L., J. Kain, S. Weiss, and D. Bright, 2005: A Brief Description of the Supercell Detection
Index. http://www.spc.noaa.gov/exper/Spring_2005/SDI-docs.pdf.
Zängl, G., D. Reinert, P. Ripodas, and M. Baldauf, 2015: The ICON (ICOsahedral Non-hydrostatic)
modelling framework of DWD and MPI-M: Description of the non-hydrostatic dynamical core. Q.J.R.
Meteorol. Soc., 141, 563–579.
Zängl, G., D. Reinert, and F. Prill, 2022: Grid refinement in icon v2.6.4. Geosci. Model Dev., 15(18),
7153–7176.
160
Glossary
ASTER Advanced Spaceborne Thermal Emission and Reflection Radiometer Global Digital Elevation
Model. 13–15
CRU-CL Climate Research Unit - Gridded climatology of 1961-1990 monthly means. 13, 15
CRU-UEA Climate Research Unit - University of east Anglia. 13
161
Glossary
162
Glossary
163
Deutscher Wetterdienst
Business Area “Research and Development”
Frankfurter Straße 135
63067 Offenbach
Germany