CSS Stylesheet |
---|
.cell_nopadding {
padding: 3px;
border: 0px;
}
.table_noborder {
border: 0px;
}
.vertical {
height: 100%;
}
.vertical td:first-child {
position: relative;
width: 35px;
padding: 0;
border: 1px solid #ddd;
}
.vertical_text {
display: block;
white-space: nowrap;
-moz-transform: rotate(-90deg);
-moz-transform-origin: center center;
-webkit-transform: rotate(-90deg);
-webkit-transform-origin: center center;
-ms-transform: rotate(-90deg);
-ms-transform-origin: center center;
position: absolute;
width:25px;
color: #205081;
}
.table_line {
padding: 0;
height: 2px;
border-left: 0px;
border-right: 0px;
}
.table_line td:first-child {
padding: 0;
height: 2px;
border-left: 0px;
border-right: 0px;
}
.font_background_cyan {
background-color: cyan;
}
.font_background_yellow {
padding: 2px;
background-color: #ffffb3;
}
.font_background_grey {
padding: 2px;
background-color: #dddddd;
}
.font_code_black {
color: #000000;
font-family: Monospace;
font-family: Menlo,Monaco,Consolas,"Courier New",monospace;
font-size: 16px;
}
.li_withspace {
margin-bottom: 10px;
} |
Introduction
Meteorological evaluation of the OpenIFS outputs is demonstrated here on case studies. Two weather events have been chosen in order to show the model's capabilities on the one hand, and on the other hand, to provide reference cases for the users (to compare the results of their installed model version with the ones run at the ECMWF or to inter-compare the results of different model versions). Important aspects at selection of the cases were their geographical scope and their physical background, i.e., those storm events were preferred which had severe impact over Europe and were governed by large-scale dynamical forcings. Cases both with weak and good operational forecast skill were considered.
Panel | ||||||||
---|---|---|---|---|---|---|---|---|
| ||||||||
|
Case studies
Storm Xaver
A large and violent cyclonic storm hit the North Sea region and several adjacent countries on 5 December 2013. Problems were caused both by the high wind speeds and the related storm surge. The surge reached 6 metres in Hamburg for example and was the highest along the England east-coast for 60 years. In the aftermaths of the cyclone a blizzard hit Sweden. The storm system was named Xaver by Berlin's Free University; other names assigned elsewhere include Bodil, Sven and St. Nicholas (Hewson et al., 2014). The cyclone developed around 00 UTC on 4 December northeast of Newfoundland and it was situated between converging northerly and southerly airstreams. Due to the westerly wind jet accelerated by the convergence, the cyclone moved to northeast and east, deepening explosively. It had an intense meso-vortex hanging back to west, which enhanced the strong wind (see Figure 1). The cyclone was presented in the operational forecasts 8-9 days before the event and the forecasts indicated the very strong wind gust 3-4 days in advance. (Although some strength overestimation over Germany as well as timing error in surface pressure were concluded.)
Figure 1: 24-hour maximum wind gust (m/s) on 5 December based on ECMWF operational IFS forecasts at 00 UTC on 3 and 5 December (left and middle, respectively; with white contours for the mean sea level pressure in hPa on 12 UTC on 5 December) and from the observations (right).
Storm Desmond
Storm Desmond caused severe flooding, travel disruption and a power outage across northern England, parts of Scotland and Ireland on 5 December 2015. Cumbria in northwestern part of England is one of the worst affected regions with more than 200 mm of rain in 24 hours recorded in that area. Storm Desmond broke the United Kingdom's 24-hour rainfall record, with 341.4 mm of rain falling in Honister Pass, Cumbria. On Saturday, 5 December, UK Met office issued a red warning of heavy rain for Cumbria. The cyclone also led to flooding in southern Norway.
Orographical enhancement of precipitation played a major role in the event and the operational model of the ECMWF picked well the highest rainfall amounts over the orographical barriers. However, the forecast underestimated the peak values of about 100 mm in 24 hours in Cumbria and overestimated the precipitation amount in lee of the hills (Figure 2).
Figure 2: 24-hour precipitation amount (mm) between 6 UTC on 5 December and 6 UTC on 6 December, based on ECMWF operational IFS forecast at 00 UTC on 5 December (left; with cyan contours for the mean sea level pressure in hPa at 12 UTC on 5 December) and observations (right).
Model experiments
Several experiments have been conducted with OpenIFS for both cases with the aim and to details of experiments are summarized in Table 1.1 Dec, 2 Dec, 3 Dec, 4 Dec, 5 Dec 2015
3 Dec, 4 Dec 2015
1 Dec, 2 Dec, 3 Dec, 4 Dec, 5 Dec 2015
3 Dec, 4 Dec 2015
3 Dec, 4 Dec 2015
!! The input data and the namelists needed to run the experiments can be downloaded from the ECMWF FTP server: download.ecmwf.int/test-data. !!
The files are packed in .tar.gz files and structured into directories named after the case studies (i.e., Xaver_201312, Desmond_201512) and subdirectories indicating the main experiment characteristics (e.g., T255L91_ERA5). The archive files were prepared by starting dates (e.g., gs0c_2015120300.tar.gz). Typical content of a .tar.gz file is as follows:
borderColor | #cccccc |
---|---|
borderWidth | 1 |
borderStyle | solid |
title | Sample directory list of an archive file |
File name
Description Size
: directory for a given starting date
0 MB2015120300/ecmwf/ : subdirectory containing the namelists and some outputs
0 MB2015120300/ecmwf/namelistfc : namelist with detailed experiment setup
4.6 KB 2015120300/ecmwf/NODE.001_01.model.1 : text output (log) file including all the important information about the model run 3.8 MB
2015120300/ecmwf/ifs.stat.model.1 : information about model steps (useful for debugging)
26 KB 2015120300/ecmwf/wam_namelist : namelist of the coupled wave
3.1 KB 2015120300/ICMCLgs0cINIT : input file containing surface and soil information (albedo, soil temperature etc.) 9.3 MB
2015120300/ICMGGgs0cINIT : input file containing gridpoint surface initial data
7 MB 2015120300/ICMSHgs0cINIT : input file containing initial data for the prognostic variables in spectral representation 35 MB
2015120300/ICMGGgs0cINIUA : input file containing initial data for the prognostic variables in gridpoint representation 101 MB
2015120300/wam_grid_tables : model grid and tables for the wave model
52 MB 2015120300/wam_subgrid_0 : information for model advection, including sub-grid parametrisation for the wave model
12 MB
2015120300/wam_subgrid_1 : information for model advection, including sub-grid parametrisation for the wave model
25 MB
2015120300/wam_subgrid_2 : information for model advection, including sub-grid parametrisation for the wave model 25 MB
2015120300/cdwavein : initial value of drag coefficient for the wave model
63 KB
2015120300/specwavein : initial wave spectra for the wave model
7.6 MB 2015120300/uwavein : initial value of wind speed for the wave model
63 KB
2015120300/sfcwindin : initial value of 10-meter horizontal wind components and sea ice fraction for the wave model 2.2 MB
The namelist file highlighted in green in the box above is needed to run OpenIFS. It controls the necessary settings (e.g., time step, experiment ID) as well as the post-processing. The most important namelist elements are listed below with their explanation:
Panel | ||||||||
---|---|---|---|---|---|---|---|---|
| ||||||||
|
In the evaluation, we would like to investigate the following variables:
- 2-meter temperature: its GRIB code number is 167;
- precipitation: it is composed from large-scale and convective precipitation with code numbers 142 and 143, respectively;
- mean sea level pressure: its code number is 151;
- 10-meter wind gust: its code number is 49;
- temperature at 850 hPa level: its code number is 130;
- relative humidity at 700 hPa level: its code number is 157;
- geopotential at 500 hPa level: its code number is 129;
u and v wind components at 250 and 100 hPa with code numbers 131 and 132, respectively.
HTML |
---|
<p>The listed variables have to be included in the namelist among the post-processing variables (see the code numbers and levels highlighted with <b>bold</b> characters in the box above). More information about the namelist settings and GRIB field codes can be found in the <b><a href="https://software.ecmwf.int/wiki/display/OIFS/How-to+articles" target="_blank">how-to articles</a></b>: <a href="https://software.ecmwf.int/wiki/display/OIFS/How+to+control+OpenIFS+output" target="_blank">How to control OpenIFS output</a>.</p>
<p><b><font color=red>!!</font></b> To run the model, the paths of the input data and the namelist have to be set in the job. ..... <a href="https://software.ecmwf.int/wiki/display/OIFS/5.6+Acceptance+testing+OpenIFS+after+installation" target="_blank">Acceptance testing OpenIFS after installation</a>.</p> |
Preparation of data for visualization
After the job run is completed, the outputs are split into 2 files (located in the working directory): ICMGG
for the gridpoint fields and ICMSH
for the spectral ones. The two files should include all the necessary variables set in the namelist. However, before visualisation of the results, some steps are still needed.
Post-processing of model outputs
The Metview macros prepared for visualisation of experiment results desire the meteorological variables in GRIB files separated by variables and day with the following names:
t2_${day}.grib
for 2-meter temperature,p_${day}.grib
for precipitation,mslp_${day}.grib
for mean sea level pressure,gust_${day}.grib
for 10-meter wind gust,t850_${day}.grib
for temperature at 850 hPa,q700_${day}.grib
for relative humidity at 700 hPa,z500_${day}.grib
for geopotential at 500 hPa,u250_${day}.grib
for horizontal wind components at 250 hPa,u100_${day}.grib
for horizontal wind components at 100 hPa,
where day is the given date in format of yyyymmdd (e.g., 20151203). All the files should contain data for 8 timesteps per day (i.e., in every 3 hours).
Note |
---|
Please note the lowercase letters in the filenames. |
The 2-meter temperature, the precipitation, the mean sea level pressure and the wind gust are expected in gridpoint representation, so will be taken from the ICMGG files. To prepare the needed input files for the Metview macros, the next operations are needed on the raw ICMGG outputs:
Panel |
---|
|
where expID is the 4-digit experiment ID and step is the post-processing step (every 3 hours). For precipitation, both the convective and large-scale precipitation components have to be in the same file:
Panel |
---|
|
The pressure level data are required in spectral representation, so they will be taken from ICMSH files. To prepare the needed input files for the Metview macros, the next operations are needed on the raw ICMSH outputs:
Panel |
---|
|
For wind, both the u and v components have to be in the same file:
Panel |
---|
|
The size of the resulted files varies by the spatial resolution and the representation of the data. For instance, the file size at T255L91 resolution is 10 MB and 8 MB per variables for gridpoint and spectral fields, respectively, whereas these values increases to 35 MB and 26 MB at T639L137, to 233 MB and 179 MB at T1279L137.
Preparation of reference data
As reference data, we will apply the ECMWF re-analyses. Both ERA-Interim and ERA5 datasets are available for the users and can be downloaded from the ECMWF MARS (Meteorological Archival and Retrieval System) system. Re-analyses are created by optimal combination of available observational information and short-range numerical weather predictions using data assimilation techniques. They provide the most comprehensive description of the past and current states of the 3-dimensional atmosphere or the Earth system.
ERA-Interim dataset (Berrisford et al., 2011) was prepared on 79 km horizontal resolution with 60 vertical levels starting from 1979. Analysis fields were constructed in every 6 hours using variety of observations (conventional measurements, remote sensing data, air craft measurements etc.), the 4D-Var data assimilation technique and the IFS model version which was operational in 2009 (cycle 31r2). The forecasts initialized from the analysis produced 3 hourly outputs up to 24 hours.
HTML |
---|
<p><b>ERA5</b> (<a href="#hersbach">Hersbach and Dee</a>, 2016) is being constructed on higher, 32 km horizontal resolution with 137 vertical levels from 1950. Analysis fields are is being prepared hourly with inclusion of newly reprocessed observational data, using the 4D-Var data assimilation technique and the IFS cycle 42r1 model version. ERA5 forecasts initialized from the hourly analyses produce hourly outputs up to 18 hours and give an estimation of forecast uncertainty. There is an important difference between ERA-Interim and ERA5 in handling of the accumulated parameters: in ERA5 the accumulation is calculated from the previous post-processing step (i.e., along one hour), while in ERA-Interim it is from the beginning of the forecast – this feature will have relevance in evaluation of the precipitation amount and wind gust. (More information about the characteristics of ERA-Interim and ERA5 can be found in the <b><a href="https://software.ecmwf.int/wiki/display/CKB/Copernicus+Knowledge+Base" target="_blank">Copernicus Knowledge Base</a></b>: <a href="https://software.ecmwf.int/wiki/pages/viewpage.action?pageId=74764925" target="_blank">What are the changes from ERA-Interim to ERA5?</a>)</p> |
The Metview macros prepared for visualisation of re-analysis data desire the meteorological variables in GRIB files separated by variables for a time period with the following names:
${dataset}
_t2_${period}.grib${dataset}
_${period}
.grib${dataset}
_${period}
.grib${dataset}
_${period}
.grib${dataset}
_${period}
.grib${dataset}
_${period}
.grib${dataset}
_${period}
.grib${dataset}
_${period}
.grib${dataset}
_${period}
.grib
where dataset is a 2-digit identifier of the re-analysis data, ei for ERA-Interim and eafor ERA5; period is the investigated time period in format of yyyymmdd-yyyymmdd (e.g., 20151201-20151206 for Desmond).
Note |
---|
Please note again the lowercase letters in the filenames. Furthermore, the re-analysis data should not be split by day, because data for the whole period will be handled together by the Metview macros. |
The 2-meter temperature, the precipitation, the mean sea level pressure and the wind gust are expected in gridpoint representation, while the pressure level data are required in spectral representation. Total precipitation and wind gust as accumulated parameters derive from forecasts, all the other variables are real analyses. Consequently, the daily quantities for precipitation and wind gust are composed of 8 and 24 timesteps from ERA-Interim and ERA5, respectively, the other variables have 4 and 8 timesteps (recall that output frequency of the forecast experiment is 3 hours). Besides the two (large-scale and convective) precipitation components, total precipitation is also available for direct retrieve both in ERA-Interim and ERA5, with GRIB code 228.
HTML |
---|
<p>
To download the necessary data from MARS, the following steps have to be accomplished:
<ol>
<li>To have access to the ECMWF public datasets, you will need to have an account on ECMWF web site: <a href="https://apps.ecmwf.int/registration/" target="_blank">https://apps.ecmwf.int/registration/</a>.
<li>To retrieve data from MARS, you will need to download and install an ECMWF key. This page shows that step by step: <a href="https://software.ecmwf.int/wiki/display/WEBAPI/Access+ECMWF+Public+Datasets" target="_blank">Access ECMWF Public Datasets</a>.
<li>Having your key, you can run the <a href="https://software.ecmwf.int/stash/projects/OIFS/repos/oifs40r1/browse/casestudies/mars/scr_download_re-analysis?at=evaluation" target="_blank"><font style="font-family:monospace;">scr_download_re-analysis</font></a> shell script available in the <a href="https://software.ecmwf.int/stash/projects/OIFS/repos/oifs40r1/browse" target="_blank">OpenIFS repository (for cycle 40r1)</a> and using <font style="font-family:monospace;">getmars</font> for retrieve.
</ol>
</p> |
The re-analysis source (i.e., ERA-Interim or ERA5), the surface and pressure level variables to be retrieved, the period of the data and the output directory have to/can be specified for the scr_download_re-analysis script. Calling it with -h option, it gives a detailed help with some examples at the end (or calling it without any option, it also gives short instructions to its configuration):
Code Block | ||
---|---|---|
| ||
% ./scr_download_re-analysis -h
----------------------------------------------------------------------
This script downloads surface and pressure level ERA-Interim or ERA5 re-analysis data for a given time range from MARS.
Usage:
-c[class],-s[surface_variables],-p[plevel_variables],-f[firstdate],-l[lastdate],-o[output_directory]
-h/-help
Examples:
./scr_download_re-analysis -cei -s"t2 p mslp gust" -p"t850 q700 z500 u250 u100" -f20151203 -l20151205 -o"../reference"
./scr_download_re-analysis -cea -sall -p"t850 q700 z500 u250 u100" -f20151203 -l20151205 -o"/home/rd/digs/metview/paper_OIFS/input/reference"
./scr_download_re-analysis -ce5 -s" " -p"t850" -f20151203 -o"/home/rd/digs/metview/paper_OIFS/input/reference"
---------------------------------------------------------------------- |
Note |
---|
Please take into account that the script is able to retrieve only the variables listed above. For further parameters the program has to be modified manually. Furthermore, please note the all option which can be used with -s and -p switches. That makes possible to take automatically all the surface and/or pressure level variables discussed above (instead of listing them in the command line). |
The retrieved ERA-Interim fields occupy approximately 54 MB, while ERA5 fields take 700 MB for the case Desmond (i.e., for 1–6 December 2015).
Visualization with Metview
The visualization package contains the following folder and file:
borderColor | #cccccc |
---|---|
borderWidth | 1 |
borderStyle | solid |
title | Content of the visualization package |
Description macros : Metview visualization macros
definitions :
place of some external functions and macros used by the macros
input
: input directory for the input files (downloaded and prepared based on the instruction above)
figs : output directory for the figures
README :
description about the content of the visualization package
Using the Metview macros
borderColor | #cccccc |
---|---|
borderWidth | 1 |
borderStyle | solid |
title | The available Metview macros |
Description plot_raw_IC.mv : Metview macro to visualize the raw initial conditions for 2 experiments and the difference between them
plot_forecastrun.mv : Metview macro to visualize the results of forecast experiments
plot_ERAI_ERA5.mv : Metview macro to visualize the ERA-Interim or ERA5 data
raw_IC_dialog : Metview dialogue box (in macro language) to
plot_raw_IC.mv
forecastrun_dialog : Metview dialogue box (in macro language)
to plot_forecastrun.mv
ERAI_ERA5_dialog : Metview dialogue box (in macro language)
to plot_ERAI_ERA5.mv
help_plot_raw_IC : help file
to
plot_raw_IC.mv
help_plot_forecastrun : help file
to
plot_forecastrun.mv
help_plot_ERAI_ERA5 : help file
to
plot_ERAI_ERA5.mv
scr_run_macros : shell script to execute the Metview macros
The visualisation is done by the Metview macros ( *.mv files) interactively using a dialogue box or in batch mode.
Interactive dialogue box
With right click on the macro and then selecting the Execute option from the menu, the settings can be seen in a dialogue box:
Expand | ||
---|---|---|
| ||
|
title | Settings of surface and upper level parameters |
---|
surface parameters
in raw_IC_dialog : 3 variables are available for visualisation, soil level temperature 2, surface pressure, (model) orography;
model level parameters
in raw_IC_dialog : 4 variables are available for visualisation, temperature, wind, specific humidity, cloud cover;
model levels
for model level variables: from 1 (uppest) to 91 (lowest);
surface parameters
in forecastrun_dialog and ERAI_ERA5_dialog : 4 variables are available, 2-meter temperature, mean sea level pressure, precipitation, 10-meter wind gust;
pressure level parameters
in forecastrun_dialog and ERAI_ERA5_dialog : 5 variables are available, temperature at 850 hPa, geopotential at 500 hPa, wind at 250 and 100 hPa, relative humidity at 700 hPa;
pressure levels
for pressure level variables: 850, 700, 500, 250, 100 hPa;
Multiple variables (both surface and model/pressure level ones) can be selected at the same time. In case of choosing any model or pressure level parameters, selecting also (at least) one model or pressure level should not be forgotten (multiply options are possible also here).
Note |
---|
Please note that to visualise different atmospheric variables on different levels (e.g., temperature at 850 hPa and humidity at 700 hPa), the macro has to be run separately with the two settings. |
Expand | ||
---|---|---|
| ||
|
Expand | ||
---|---|---|
| ||
|
Expand | ||
---|---|---|
| ||
|
Batch mode
In batch mode the macro can be executed following the next syntax:
Panel |
---|
% metview -b macro option1 option2 option3 ... |
where macro is the macro to be run ( plot_raw_IC.mv, plot_forecastrun.mv or plot_ERAI_ERA5.mv). A detailed help together with some useful examples is provided with simple execution of the macro:
Panel |
---|
% metview -b macro |
The shell script scr_run_macros executes the macros from the Unix/Linux shell and it can be tailored for the own needs.
External functions and macros
borderColor | #cccccc |
---|---|
borderWidth | 1 |
borderStyle | solid |
title | The external functions located in definitions directory |
Description build_layout_2plus1 : layout definition with 2 left panels and 1 right panel
build_layout_single : layout definition a single panel
titlemain : title style for the main plot
titlemain_2L : 2-line title style for the main plot
titlepanels : title style for the individual panels
legend_main : legend definition for a single page
legend_shade : legend definition for left panels
legend_diff : legend definition for right panel (for the difference field)
base_visdef : colour definitions for the different variables
diff_range : dynamic colour definitions for the difference fields
Note |
---|
Please note that there are 2 include statements in the plot_raw_IC.mv and 1 in the plot_forecastrun.mv and the plot_ERAI_ERA5.mv, taking the two colour definitions (base_visdef and diff_range) from this directory. The path of the definitions directory has to be set in the downloaded macros according to the local working tree (it is necessary because using dynamic path with include is not possible in the macro language). |
Input data
The input data are requested with the following content, format and name convention:
- Macro plot_raw_IC.mv expects the raw ICM* files as input: ICMCL${expID}INIT, ICMGG${expID}INIT, ICMGG${expID}INIUA, ICMSH${expID}INIT, where expID is the 4-digit experiment ID.
- Macros plot_forecastrun.mv and plot_ERAI_ERA5.mv expect grib files as input with the following file names: ${variable}_${date}.grib, where variable can be t2, mslp, p, gust, t850, q700, z500, u250, u100; date is day in format yyyymmdd.
Output figures
All the macros produce figures in single-page .ps files. The file name can be as follows:
- plot_raw_IC.mv: ${variable}_${level}_ERAI-ERA5${date}+${timestep}.ps,
where variable can be stl2, lnsp, z, t, cc, u, q; level can be 0 (in case of surface variables) or from 1 to 91; date is day in format yyyymmdd; timestep is forecast lead time in hours, e.g., 0, 3, 6 etc. - plot_forecastrun.mv: ${variable}_${level}_${expID}_${date}+${timestep}.ps,
where variable can be t2, mslp, p, gust, t, q, z, u, q; level can be 0 (in case of surface variables) or 850, 700, 500, 250, 100; expID is the 4-digit experiment ID; date is day in format yyyymmdd; timestep is forecast lead time in hours, e.g., 0, 3, 6 etc. - plot_ERAI_ERA5.mv: ${variable}_${level}_${reference}_${date}.ps,
where variable can be t2, mslp, p, gust, t, q, z, u; level can be 0 (in case of surface variables) or 850, 700, 500, 250, 100; reference can be ERAI or ERA5; date is day in format yyyymmdd.
title | To have an overview on the large amount of figures, a catalogue can be prepared with the steps below: |
---|
- Convert the .ps files into .png files (images with 120 DPI are sharp enough with limited file size): convert -density 120 psfile pngfile
- Choose a variable (e.g., T850) and open one of the .docx files in the figs/docs directory (initial_condition_t850_table.docx). Delete all figures from the table, but keep the table itself as it is.
- Click on the
View macros
menu item in Macros menu point on View tab (in MS Office 2013) and edit the macro aainsertpicsOIFS_initialconditions. The macro source code can be seen in the opening window. The directory paths and file names can be replaced here to paths and file names set in the local environment. This replacement can be done automatically (withCTRL+H
). Afterwards the macro has to be saved. - Then return to the the main document, stand in the cell of the first picture and run the macro (with clicking on the
View macros
menu item in Macros menu point on View tab and running macro aainsertpicsOIFS_initialconditions). - Save the document at the end of the macro run. The same can be repeated for all the variables.
Info |
---|
To prepare similar document for comparison of impact of resolution and startdate, the corresponding *resolution.docx, *startdate.docx files and aainsertpicsOIFS_resolutions, aainsertpicsOIFS_startdates macros have to be used. {{holmap}} |
Info |
---|
There are separate macros for precipitation figures as 24-hour amount is processed instead of precipitation between 2 time steps (see macros aainsertprecpics*). |
References
HTML |
---|
<p><a name="hersbach"></a>Hersbach, H., Dee, D.P., 2016: <a href="https://www.ecmwf.int/en/newsletter/147/news/era5-reanalysis-production" target="_blank">ERA5 reanalysis is in production.</a> <em>ECMWF Newsletter 147,</em> p. 7.</p> |