ICON-EU grib data

Dataset specific topics involving WPS.

ICON-EU grib data

Postby beni » Sat Feb 03, 2018 9:24 am

Hello,
I've been trying to use the german ICON-EU model's grib files for my WRF runs. The free output from the model is updated every 3 hours (+120 hours at 00z,06z,12z,18z and +30 hours az 03z,09z,15z,21z) with 6.5 km resolution over Europe (data can be accessed here: https://opendata.dwd.de/weather/icon/eu_nest/grib/). Each parameter has its separate grib files, which has to be concatenated before ungrib can process it.

My only problem is with the soil data - ICON has it for the following levels:
Temperature: 0 cm, 6 cm and 1458 cm
Moisture is column-integrated: 0-1 cm, 1-3 cm, 3-9 cm, 9-27 cm, 27-81 cm, 81-243 cm, 243-729 cm, 729-2187 cm.

Is there a proper way to use it in WRF? Currently I bind these values to the nearest metgrid levels, which is not correct but shown to be enough for RUC land physics (Noah and Noah-MP crashes). Maybe an official Vtable would be fine for the next WRF release as ICON-EU seems to be the best free data source for realtime european WRF runs, with the highest resolution and most frequent refresh rate.

Code: Select all
GRIB1| Level| From |  To  | metgrid  | metgrid | metgrid                                 |GRIB2|GRIB2|GRIB2|GRIB2|
Param| Type |Level1|Level2| Name     | Units   | Description                             |Discp|Catgy|Param|Level|
-----+------+------+------+----------+---------+-----------------------------------------+-----------------------+
  11 | 100  |   *  |      | TT       | K       | Temperature                             |  0  |  0  |  0  | 100 |
  33 | 100  |   *  |      | UU       | m s-1   | U                                       |  0  |  2  |  2  | 100 |
  34 | 100  |   *  |      | VV       | m s-1   | V                                       |  0  |  2  |  3  | 100 |
  52 | 100  |   *  |      | RH       | %       | Relative Humidity                       |  0  |  1  |  1  | 100 |
  11 | 105  |   2  |      | TT       | K       | Temperature       at 2 m                |  0  |  0  |  0  | 103 |
  33 | 105  |  10  |      | UU       | m s-1   | U                 at 10 m               |  0  |  2  |  2  | 103 |
  34 | 105  |  10  |      | VV       | m s-1   | V                 at 10 m               |  0  |  2  |  3  | 103 |
  52 | 105  |   2  |      | RH       | %       | Relative Humidity at 2 m                |  0  |  1  |  1  | 103 |
   1 |   1  |   0  |      | PSFC     | Pa      | Surface Pressure                        |  0  |  3  |  0  |   1 |
   2 | 102  |   0  |      | PMSL     | Pa      | Sea-level Pressure                      |  0  |  3  |  1  | 101 |
 144 | 111  |   0  |      | SOILM000 | kg m-2  | Soil Moist 0 cm                         |  2  |  3  | 20  | 106 |
 144 | 111  |   1  |      | SOILM001 | kg m-2  | Soil Moist 1 cm                         |  2  |  3  | 20  | 106 |
 144 | 111  |   3  |      | SOILM004 | kg m-2  | Soil Moist 3 cm                         |  2  |  3  | 20  | 106 |
 144 | 111  |   9  |      | SOILM010 | kg m-2  | Soil Moist 9 cm                         |  2  |  3  | 20  | 106 |
 144 | 111  |  27  |      | SOILM030 | kg m-2  | Soil Moist 27 cm                        |  2  |  3  | 20  | 106 |
 144 | 111  |  81  |      | SOILM100 | kg m-2  | Soil Moist 81 cm                        |  2  |  3  | 20  | 106 |
 144 | 111  | 243  |      | SOILM160 | kg m-2  | Soil Moist 243 cm                       |  2  |  3  | 20  | 106 |
 144 | 111  | 729  |      | SOILM300 | kg m-2  | Soil Moist 729 cm                       |  2  |  3  | 20  | 106 |
  11 | 111  |   0  |      | SOILT000 | K       | Soil Temperature 0 cm                   |  2  |  3  | 18  | 106 |
  11 | 111  |   6  |      | SOILT005 | K       | Soil Temperature 6 cm                   |  2  |  3  | 18  | 106 |
  11 | 111  |1458  |      | SOILT300 | K       | Soil Temperature 1458 cm                |  2  |  3  | 18  | 106 |
  81 |   1  |   0  |      | LANDSEA  | proprtn | Land/Sea flag (1=land, 0 or 2=sea)      |  2  |  0  |  0  |   1 |
   7 |   1  |   0  |      | SOILHGT  | m       | Terrain field of source analysis        |  0  |  3  |  6  |   1 |
  11 |   1  |   0  |      | SKINTEMP | K       | Skin temperature                        |  0  |  0  |  0  |   1 |
  65 |   1  |   0  |      | SNOW     | kg m-2  | Snow depth water equivalent             |  0  |  1  | 60  |   1 |
     |   1  |   0  |      | SNOWH    | m       | Physical Snow Depth                     |  0  |  1  | 11  |   1 |
-----+------+------+------+----------+---------+-----------------------------------------+-----------------------+
beni
 
Posts: 22
Joined: Tue Apr 22, 2014 10:52 am

Re: ICON-EU grib data

Postby Gippox » Sat Feb 17, 2018 6:28 pm

Hello Beni,
i tried to initialize the wrf-arw with the ICON-EU grib2 (RUC land surface model) but real.exe gives me this error:

Using sfcprs to compute psfc
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE: <stdin> LINE: 4900
troubles, could not find trapping x locations
-------------------------------------------
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0

I used the VTable you posted. How did you solve it?

Thanks.
Gippox
 
Posts: 9
Joined: Fri May 10, 2013 8:30 am

Re: ICON-EU grib data

Postby beni » Sat Feb 24, 2018 1:45 pm

Hello, how much is your "p_top_requested" and "e_vert" in the namelist?
Are you using the pressure level or the model level GRIB's?
beni
 
Posts: 22
Joined: Tue Apr 22, 2014 10:52 am

Re: ICON-EU grib data

Postby beni » Sun Mar 11, 2018 7:38 am

I got a reply from wrfhelp which I share with all of you -

Convert the soil temperature and moisture data into the intermediate format on their native vertical levels, giving the fields names like 'ICON_SM1', 'ICON_SM2', etc. for the levels, from top to bottom. The exact names of the 2-d layers of data are not important, since the 2-d layers will be combined into a 3-d array.

Then, add lines like this in the METGRID.TBL entry for 'SM':

fill_lev = 1 : ICON_SM1(200100)
fill_lev = 2 : ICON_SM2(200100)
...
fill_lev = 8 : ICON_SM10(200100)

it will also be necessary to add entries for 'ICON_SM1', etc., in the METGRID.TBL file to describe how these fields will be horizontally interpolated following what is done for, e.g., SM000010.

Running metgrid should create a 3-d 'SM' field, as well as a 3-d 'SOIL_LAYERS' field, that contains the values 1, 2, ..., 8 for each of the 8 layers. These are of course not the real bottom depths of the ICON soil layers; so, it will then be necessary to edit the met_em netCDF files to fill in the actual bottom depths in the 'SOIL_LAYER' array with the actual depths in units of centimeters. Since SOIL_LAYERS is a real-valued array, the values can be, e.g., 243 cm for the layer whose bottom is at 243 cm.

After editing the SOIL_LAYERS field, I think the real program should work as usual.
beni
 
Posts: 22
Joined: Tue Apr 22, 2014 10:52 am


Return to Working with Various Datasets

Who is online

Users browsing this forum: No registered users and 1 guest