Scheduled Downtime
On Friday 21 April 2023 @ 5pm MT, this website will be down for maintenance and expected to return online the morning of 24 April 2023 at the latest

ECMWF open data 0p25 grib2 Vtable?

M

meteoadriatic

Member
Hello,

Do we have Vtable suitable for ungribbing this new dataset?

Code: 
https://data.ecmwf.int/forecasts/<YYYYMMDD>/00z/ifs/0p25/oper/

Thank you.
 
I'm not sure if the existing Vtables (perhaps Vtable.ECMWF?) are suitable. If you'd like to test them out and let us know, that would be great! If they don't work, you can try contacting the input data support group to see if they are able to help. You can also use the utilities, such as g2print.exe and g1print.exe to look at the codes given in the new data type, which can help you to modify an existing Vtable. If you are able to find one that works, feel free to share it so that it may help other users who are trying to use these data.
 
Hello Kelly,

Provided Vtable.ECMWF unfortunatelly doesn't work as it expects grib1 format. I tried to convert grib2 to grib1 with cnvgrib but it fails too.

I hoped that someone already has new Vtable, but ok, it looks I will need to create new one using tools you mentioned. Thank you anyway. If I succeed I will post here working solution later.

Ivan
 
Hello again,

Is there anything we can do about this problem?


Code: 
./g2print.exe 20240402000000-0h-oper-fc.grib2
 ungrib - grib edition num           2
 reading from grib file =
 20240402000000-0h-oper-fc.grib2                                               
                                          
 getdrstemplate: DRS Template           42  not defined.
      ECMWF                           
---------------------------------------------------------------------------------------
 rec Prod Cat Param  Lvl    Lvl      Lvl     Prod    Name            Time          Fcst
 num Disc     num    code   one      two     Templ                                 hour
---------------------------------------------------------------------------------------
 getdrstemplate: DRS Template           42  not defined.
  ERROR extracting field gf_getfld =           12
 getdrstemplate: DRS Template           42  not defined.
  ERROR extracting field gf_getfld =           12
 getdrstemplate: DRS Template           42  not defined.
  ERROR extracting field gf_getfld =           12
 getdrstemplate: DRS Template           42  not defined.
  ERROR extracting field gf_getfld =           12
 getdrstemplate: DRS Template           42  not defined.
  ERROR extracting field gf_getfld =           12
 getdrstemplate: DRS Template           42  not defined.

These error messages repeat until the end of the output, so I don't get anything useful with g2print.exe tool.

Thank you.
 
I think this could be related to CCSDS compression. You can change the packing type by using eccodes for instance:

grib_set -r -w packingType=grid_ccsds -s packingType=grid_simple input.grb output.grb

Then try again with the uncompressed gribs.
 
use wgrib2 IN.grb -set_grib_type c3 -grib_out OUT.grb to change the compression format
Regarding the rest, I am also in the process of trying to use ifs data. It's almost working doing:
- download
- change soil variables with -set_var : wgrib2 test.grb2 -set_var TMP -grib out.grb2 and wgrib2 test.grb2 -set_var SOILW -grib out2.grb2
- change the level, once downloaded the level that are supposed to be cm are put in m (0-7 cm appears 0-7 m)
-change the compression format
- run ungrib with Vtable:
GRIB1| Level| From | To | metgrid | metgrid | metgrid |GRIB2|GRIB2|GRIB2|GRIB2|
Param| Type |Level1|Level2| Name | Units | Description |Discp|Catgy|Param|Level|
-----+------+------+------+----------+---------+-----------------------------------------+-----------------------+
129 | 100 | * | | HGT | K | Height | 0 | 3 | 5 | 100 |
130 | 100 | * | | TT | K | Temperature | 0 | 0 | 0 | 100 |
157 | 100 | * | | RH | % | Relative Humidity | 0 | 1 | 1 | 100 |
131 | 100 | * | | UU | m s-1 | U | 0 | 2 | 2 | 100 |
132 | 100 | * | | VV | m s-1 | V | 0 | 2 | 3 | 100 |
167 | 105 | 2 | | TT | K | Temperature at 2 m | 0 | 0 | 0 | 103 |
151 | 102 | 0 | | MSL | Pa | Mean Sea Level Pressure | 0 | 3 | 0 | 101 |
134 | 1 | 0 | | PSFC | Pa | Surface Pressure | 0 | 3 | 0 | 1 |
165 | 105 | 10 | | UU | m s-1 | U at 10 m | 0 | 2 | 2 | 103 |
166 | 105 | 10 | | VV | m s-1 | V at 10 m | 0 | 2 | 3 | 103 |
172 | 1 | 0 | | LANDSEA | proprtn | Land Sea Mask 0/1 Flag | 2 | 0 | 0 | 1 |
139 | 112 | 0 | 7 | ST000007 | K | T 0-07 cm below ground layer | 0 | 0 | 0 | 106 |
170 | 112 | 7 | 28 | ST007028 | K | T 07-28 cm below ground layer | 0 | 0 | 0 | 106 |
183 | 112 | 28 | 100 | ST028100 | K | T 28-100 cm below ground layer | 0 | 0 | 0 | 106 |
236 | 112 | 100 | 255 | ST100255 | K | T100-255 cm below ground layer | 0 | 0 | 0 | 106 |
39 | 112 | 0 | 7 | SM000007 | fraction| Soil Moist 0-07 cm below grn layer | 2 | 0 | 9 | 106 |
40 | 112 | 7 | 28 | SM007028 | fraction| Soil Moist 7-28 cm below grn layer | 2 | 0 | 9 | 106 |
41 | 112 | 28 | 100 | SM028100 | fraction| Soil Moist 028-100 cm below grn layer | 2 | 0 | 9 | 106 |
42 | 112 | 100 | 255 | SM100255 | fraction| Soil Moist 100-255 cm below grn layer | 2 | 0 | 9 | 106 |
-----+------+------+------+----------+---------+-----------------------------------------+-----------------------+
 
kwerner said:
I'm not sure if the existing Vtables (perhaps Vtable.ECMWF?) are suitable. If you'd like to test them out and let us know, that would be great! If they don't work, you can try contacting the input data support group to see if they are able to help. You can also use the utilities, such as g2print.exe and g1print.exe to look at the codes given in the new data type, which can help you to modify an existing Vtable. If you are able to find one that works, feel free to share it so that it may help other users who are trying to use these data.
Click to expand...
Could you please check the data from Data Sample link and tell us the solution?
 
Avinash said:
Have you got Vtable for the IFS data?
Click to expand...
The one that is posted by user laetis should be it. But I'm not able to proceed further because metgrid crashes, apparently data inside intermediate files is still missing some mandatory fields. I gave up until someone comes up with complete solution.
 
forum.mmm.ucar.edu

What fields are mandatory for running WRF?

The following are mandatory input data fields: 3D Data (for data on constant pressure levels) Temperature U and V components of wind Geopotential height Relative humidity (or specific humidity) 3D Data (for data on native model levels) Temperature U and V components of wind Geopotential...
forum.mmm.ucar.edu forum.mmm.ucar.edu

Looking at the data provided by ECMWF they are missing the mandatory fields.


Code: 
3D Data (for data on constant pressure levels)
Temperature
U and V components of wind
Geopotential height
Relative humidity (or specific humidity)


3D Data (for data on native model levels)
Temperature
U and V components of wind
Geopotential height
Relative humidity (or specific humidity)
Pressure


2D Data
Surface pressure
Mean sea-level pressure
Skin temperature/SST
2 meter temperature
2 meter relative humidity
10 meter U and V components of wind
Soil temperature
Soil moisture
Soil height (or terrain height)
*Note: the 2m temperature, RH, and wind fields may be optional if one sets the namelist option use_surface = F before running real.exe)
 
laetis said:
use wgrib2 IN.grb -set_grib_type c3 -grib_out OUT.grb to change the compression format
Regarding the rest, I am also in the process of trying to use ifs data. It's almost working doing:
- download
- change soil variables with -set_var : wgrib2 test.grb2 -set_var TMP -grib out.grb2 and wgrib2 test.grb2 -set_var SOILW -grib out2.grb2
- change the level, once downloaded the level that are supposed to be cm are put in m (0-7 cm appears 0-7 m)
-change the compression format
- run ungrib with Vtable:
GRIB1| Level| From | To | metgrid | metgrid | metgrid |GRIB2|GRIB2|GRIB2|GRIB2|
Param| Type |Level1|Level2| Name | Units | Description |Discp|Catgy|Param|Level|
-----+------+------+------+----------+---------+-----------------------------------------+-----------------------+
129 | 100 | * | | HGT | K | Height | 0 | 3 | 5 | 100 |
130 | 100 | * | | TT | K | Temperature | 0 | 0 | 0 | 100 |
157 | 100 | * | | RH | % | Relative Humidity | 0 | 1 | 1 | 100 |
131 | 100 | * | | UU | m s-1 | U | 0 | 2 | 2 | 100 |
132 | 100 | * | | VV | m s-1 | V | 0 | 2 | 3 | 100 |
167 | 105 | 2 | | TT | K | Temperature at 2 m | 0 | 0 | 0 | 103 |
151 | 102 | 0 | | MSL | Pa | Mean Sea Level Pressure | 0 | 3 | 0 | 101 |
134 | 1 | 0 | | PSFC | Pa | Surface Pressure | 0 | 3 | 0 | 1 |
165 | 105 | 10 | | UU | m s-1 | U at 10 m | 0 | 2 | 2 | 103 |
166 | 105 | 10 | | VV | m s-1 | V at 10 m | 0 | 2 | 3 | 103 |
172 | 1 | 0 | | LANDSEA | proprtn | Land Sea Mask 0/1 Flag | 2 | 0 | 0 | 1 |
139 | 112 | 0 | 7 | ST000007 | K | T 0-07 cm below ground layer | 0 | 0 | 0 | 106 |
170 | 112 | 7 | 28 | ST007028 | K | T 07-28 cm below ground layer | 0 | 0 | 0 | 106 |
183 | 112 | 28 | 100 | ST028100 | K | T 28-100 cm below ground layer | 0 | 0 | 0 | 106 |
236 | 112 | 100 | 255 | ST100255 | K | T100-255 cm below ground layer | 0 | 0 | 0 | 106 |
39 | 112 | 0 | 7 | SM000007 | fraction| Soil Moist 0-07 cm below grn layer | 2 | 0 | 9 | 106 |
40 | 112 | 7 | 28 | SM007028 | fraction| Soil Moist 7-28 cm below grn layer | 2 | 0 | 9 | 106 |
41 | 112 | 28 | 100 | SM028100 | fraction| Soil Moist 028-100 cm below grn layer | 2 | 0 | 9 | 106 |
42 | 112 | 100 | 255 | SM100255 | fraction| Soil Moist 100-255 cm below grn layer | 2 | 0 | 9 | 106 |
-----+------+------+------+----------+---------+-----------------------------------------+-----------------------+
Click to expand...
Hi Laetis,

Thanks for the info! Any updates on this? Are you now able to run ungrib and metgrid with IFS data?
 
laetis said:
Hello,
Yes it works fine.
Click to expand...
Laetis, thanks for the quick reply. I'm trying to follow your strategy:

» Ran wgrib2 IN.grb -set_grib_type c3 -grib_out OUT.grb to change the compression format. After this, I ran ungrib with the Vtable you supplied above and it went through without errors, but seems like the soil vars (temp and humidity at 4 levels) are not being read by ungrib.

» So, I tried to to as you explain, to change soil variables with: wgrib2 test.grb2 -set_var TMP -grib out.grb2
wgrib2 test.grb2 -set_var SOILW -grib out2.grb2.

But I'm with some doubts about this step, because seems like from the commands you suggest I'll get 2 output files: out.grib2 and out2.grib2. Can you please explain how can I change the soil vars, and also how to change the levels units from meters to cm?

Thanks in advance for your help
 
dmoreno said:
set_var
Click to expand...
Hi,

Here an update with more explanation on what I posted before:
- download
- change soil variables with -set_var and -set_level (I noticed when running wgrib2 -s on the data downloaded that the soil depths appear in meters instead of cm! no idea why it's doing that but obviously ungrib cannot find the data with the wrong depth! set_level allows to change it.
wgrib2 test.grb2 -set_var TMP -grib temp1.grb2
wgrib2 temp1.grb2 -if ":0-7 m below ground" -set_lev "0-0.07 m below ground" -elseif ":7-28 m below ground" -set_lev "0.07-0.28 m below ground" -elseif ":28-100 m below ground" -set_lev "0.28-1 m below ground" -elseif ":100-" -set_lev "1-2.89 m below ground" -endif -grib res.grb2
wgrib2 test2.grb2 -set_var SOILW -grib temp1.grb2
wgrib2 temp1.grb2 -if ":0-7 m below ground" -set_lev "0-0.07 m below ground" -elseif ":7-28 m below ground" -set_lev "0.07-0.28 m below ground" -elseif ":28-100 m below ground" -set_lev "0.28-1 m below ground" -elseif ":100-" -set_lev "1-2.89 m below ground" -endif -grib res2.grb2
-change the compression format:
wgrib2 ifs_file.grb2 -set_grib_type c3 -grib_out newfilename
- run ungrib with Vtable (sightly different from previous post after reading more carefully the documentation I had some depths wrong):
GRIB1| Level| From | To | metgrid | metgrid | metgrid |GRIB2|GRIB2|GRIB2|GRIB2|
Param| Type |Level1|Level2| Name | Units | Description |Discp|Catgy|Param|Level|
-----+------+------+------+----------+---------+-----------------------------------------+-----------------------+
129 | 100 | * | | HGT | K | Height | 0 | 3 | 5 | 100 |
130 | 100 | * | | TT | K | Temperature | 0 | 0 | 0 | 100 |
157 | 100 | * | | RH | % | Relative Humidity | 0 | 1 | 1 | 100 |
131 | 100 | * | | UU | m s-1 | U | 0 | 2 | 2 | 100 |
132 | 100 | * | | VV | m s-1 | V | 0 | 2 | 3 | 100 |
167 | 105 | 2 | | TT | K | Temperature at 2 m | 0 | 0 | 0 | 103 |
151 | 102 | 0 | | PMSL | Pa | Mean Sea Level Pressure | 0 | 3 | 0 | 101 |
134 | 1 | 0 | | PSFC | Pa | Surface Pressure | 0 | 3 | 0 | 1 |
165 | 105 | 10 | | UU | m s-1 | U at 10 m | 0 | 2 | 2 | 103 |
166 | 105 | 10 | | VV | m s-1 | V at 10 m | 0 | 2 | 3 | 103 |
172 | 1 | 0 | | LANDSEA | proprtn | Land Sea Mask 0/1 Flag | 2 | 0 | 0 | 1 |
139 | 112 | 0 | 7 | ST000007 | K | T 0-07 cm below ground layer | 0 | 0 | 0 | 106 |
170 | 112 | 7 | 28 | ST007028 | K | T 07-28 cm below ground layer | 0 | 0 | 0 | 106 |
183 | 112 | 28 | 100 | ST028100 | K | T 28-100 cm below ground layer | 0 | 0 | 0 | 106 |
236 | 112 | 100 | 289 | ST100289 | K | T100-289 cm below ground layer | 0 | 0 | 0 | 106 |
39 | 112 | 0 | 7 | SM000007 | fraction| Soil Moist 0-07 cm below grn layer | 2 | 0 | 9 | 106 |
40 | 112 | 7 | 28 | SM007028 | fraction| Soil Moist 7-28 cm below grn layer | 2 | 0 | 9 | 106 |
41 | 112 | 28 | 100 | SM028100 | fraction| Soil Moist 028-100 cm below grn layer | 2 | 0 | 9 | 106 |
42 | 112 | 100 | 289 | SM100289 | fraction| Soil Moist 100-289 cm below grn layer | 2 | 0 | 9 | 106 |
-----+------+------+------+----------+---------+-----------------------------------------+-----------------------+
Then it should be fine.

About your remark on the files names, I have the ifs data in several files (test.grb2 and test2.grb2 in my example) and after my modification I cat all the results files (res.grb2 and res2.grb2 in my example). If you have only 1 file apply all change making the previous output file the next input file:
wgrib2 test.grb2 -set_var TMP -grib temp1.grb2
wgrib2 temp1.grb2 -if ":0-7 m below ground" -set_lev "0-0.07 m below ground" -elseif ":7-28 m below ground" -set_lev "0.07-0.28 m below ground" -elseif ":28-100 m below ground" -set_lev "0.28-1 m below ground" -elseif ":100-" -set_lev "1-2.89 m below ground" -endif -grib res.grb2
wgrib2 res.grb2 -set_var SOILW -grib temp1.grb2
wgrib2 temp1.grb2 -if ":0-7 m below ground" -set_lev "0-0.07 m below ground" -elseif ":7-28 m below ground" -set_lev "0.07-0.28 m below ground" -elseif ":28-100 m below ground" -set_lev "0.28-1 m below ground" -elseif ":100-" -set_lev "1-2.89 m below ground" -endif -grib res.grb2

To check everything is all right, run wgrib2 -s on your result file, the output should be something like:
1:0:d=2024041100:TMP:0-0.07 m below ground:120 hour fcst:
2:696788:d=2024041100:TMP:0.28-1 m below ground:120 hour fcst:
3:1382077:d=2024041100:TMP:0.07-0.28 m below ground:120 hour fcst:
4:2076965:d=2024041100:TMP:1-2.89 m below ground:120 hour fcst:
5:2759813:d=2024041100:SOILW:0-0.07 m below ground:120 hour fcst:
6:3230018:d=2024041100:SOILW:0.07-0.28 m below ground:120 hour fcst:
7:3697958:d=2024041100:SOILW:0.28-1 m below ground:120 hour fcst:
8:4167993:d=2024041100:SOILW:1-2.89 m below ground:120 hour fcst:
9:4566592:d=2024041100:HGT:1000 mb:120 hour fcst:
10:5053854:d=2024041100:HGT:700 mb:120 hour fcst:
11:5434160:d=2024041100:HGT:250 mb:120 hour fcst:
12:5786053:d=2024041100:HGT:100 mb:120 hour fcst:
13:6135174:d=2024041100:HGT:400 mb:120 hour fcst:
14:6467120:d=2024041100:HGT:300 mb:120 hour fcst:
15:6810162:d=2024041100:HGT:850 mb:120 hour fcst:
16:7219145:d=2024041100:HGT:200 mb:120 hour fcst:
17:7574065:d=2024041100:HGT:500 mb:120 hour fcst:
18:7895701:d=2024041100:HGT:600 mb:120 hour fcst:
19:8275293:d=2024041100:HGT:925 mb:120 hour fcst:
20:8717764:d=2024041100:HGT:150 mb:120 hour fcst:
21:9073792:d=2024041100:HGT:50 mb:120 hour fcst:
22:9395573:d=2024041100:TMP:500 mb:120 hour fcst:
23:9964163:d=2024041100:TMP:300 mb:120 hour fcst:
24:10508422:d=2024041100:TMP:250 mb:120 hour fcst:
25:11070179:d=2024041100:TMP:400 mb:120 hour fcst:
26:11627248:d=2024041100:TMP:100 mb:120 hour fcst:
27:12239054:d=2024041100:TMP:200 mb:120 hour fcst:
28:12923766:d=2024041100:TMP:850 mb:120 hour fcst:
29:13529006:d=2024041100:TMP:1000 mb:120 hour fcst:
30:14166396:d=2024041100:TMP:700 mb:120 hour fcst:
31:14701799:d=2024041100:TMP:925 mb:120 hour fcst:
32:15313958:d=2024041100:TMP:150 mb:120 hour fcst:
33:16008784:d=2024041100:TMP:600 mb:120 hour fcst:
34:16507143:d=2024041100:TMP:50 mb:120 hour fcst:
35:17213155:d=2024041100:RH:300 mb:120 hour fcst:
36:17762905:d=2024041100:RH:250 mb:120 hour fcst:
37:18204670:d=2024041100:RH:1000 mb:120 hour fcst:
38:18737781:d=2024041100:RH:500 mb:120 hour fcst:
39:19276476:d=2024041100:RH:400 mb:120 hour fcst:
40:19830525:d=2024041100:RH:925 mb:120 hour fcst:
41:20397885:d=2024041100:RH:850 mb:120 hour fcst:
42:20992796:d=2024041100:RH:200 mb:120 hour fcst:
43:21456572:d=2024041100:RH:700 mb:120 hour fcst:
44:21920707:d=2024041100:RH:600 mb:120 hour fcst:
45:22467483:d=2024041100:RH:100 mb:120 hour fcst:
46:22699927:d=2024041100:RH:150 mb:120 hour fcst:
47:23059934:d=2024041100:RH:50 mb:120 hour fcst:
48:23341488:d=2024041100:UGRD:600 mb:120 hour fcst:
49:24009649:d=2024041100:UGRD:50 mb:120 hour fcst:
50:24601595:d=2024041100:UGRD:150 mb:120 hour fcst:
51:25246639:d=2024041100:UGRD:100 mb:120 hour fcst:
52:25988233:d=2024041100:UGRD:400 mb:120 hour fcst:
53:26642712:d=2024041100:UGRD:200 mb:120 hour fcst:
54:27291926:d=2024041100:UGRD:250 mb:120 hour fcst:
55:27949428:d=2024041100:UGRD:300 mb:120 hour fcst:
56:28605153:d=2024041100:UGRD:925 mb:120 hour fcst:
57:29451390:d=2024041100:UGRD:850 mb:120 hour fcst:
58:30297006:d=2024041100:UGRD:1000 mb:120 hour fcst:
59:31145535:d=2024041100:UGRD:500 mb:120 hour fcst:
60:31800522:d=2024041100:UGRD:700 mb:120 hour fcst:
61:32486900:d=2024041100:VGRD:600 mb:120 hour fcst:
62:33136879:d=2024041100:VGRD:50 mb:120 hour fcst:
63:33759256:d=2024041100:VGRD:150 mb:120 hour fcst:
64:34350511:d=2024041100:VGRD:100 mb:120 hour fcst:
65:34916746:d=2024041100:VGRD:400 mb:120 hour fcst:
66:35565981:d=2024041100:VGRD:200 mb:120 hour fcst:
67:36182403:d=2024041100:VGRD:250 mb:120 hour fcst:
68:36719462:d=2024041100:VGRD:300 mb:120 hour fcst:
69:37265183:d=2024041100:VGRD:925 mb:120 hour fcst:
70:38107825:d=2024041100:VGRD:850 mb:120 hour fcst:
71:38836090:d=2024041100:VGRD:1000 mb:120 hour fcst:
72:39656988:d=2024041100:VGRD:500 mb:120 hour fcst:
73:40297303:d=2024041100:VGRD:700 mb:120 hour fcst:
74:41080214:d=2024041100: PRES:surface:120 hour fcst:
75:41688879:d=2024041100:TMP:2 m above ground:120 hour fcst:
76:42438061:d=2024041100: PRES:mean sea level:120 hour fcst:
77:42942388:d=2024041100:VGRD:10 m above ground:120 hour fcst:
78:43866285:d=2024041100:UGRD:10 m above ground:120 hour fcst:
79:44816209:d=2024041100:LAND:surface:120 hour fcst:


Best
 
Much appreciated for all your help Laetis, and I'm sorry for all my questions, but there are still some issues.

I've downloaded the file 20240603000000-0h-oper-fc.grib2, from https://data.ecmwf.int/forecasts/20240603/00z/ifs/0p25/oper/, as a test file. Let's call it test.grb2 . So:

1 » When I run the wgrib2 commands to change the TMP levels from m to cm, wgrib2 test.grb2 -set_var TMP -grib temp1.grb2, the output files comes with ALL variables called TMP:

» wgrib2 -s temp1.grb2
1:0:d=2024060300:TMP:surface:anl:
2:639113:d=2024060300:TMP:atmos col:anl:
3:1349091:d=2024060300:TMP:200 mb:anl:
4:1781578:d=2024060300:TMP:10 m above ground:anl:
5:2649879:d=2024060300:TMP:700 mb:anl:
6:3362358:d=2024060300:TMP:10 m above ground:anl:
7:4236641:d=2024060300:TMP:925 mb:anl:
8:4843790:d=2024060300:TMP:700 mb:anl:
9:5582167:d=2024060300:TMP:50 mb:anl:
10:6170924:d=2024060300:TMP:50 mb:anl:
(...)


The same happens for the SOILW... So, the resulting res.grib2 only has SOILW vars, and ungrib only sees those...
Anyway, I came up with a workaround:

1 » Change compression format:

wgrib2 test.grb2 -set_grib_type c3 -grib_out vars_except_soil.grib2

2 » Fix soil temp
wgrib2 vars_except_soil.grib2 -set_var TMP -grib temp1.grb2

wgrib2 temp1.grb2 -if ":0-7 m below ground" -set_lev "0-0.07 m below ground" -elseif ":7-28 m below ground" -set_lev "0.07-0.28 m below ground" -elseif ":28-100 m below ground" -set_lev "0.28-1 m below ground" -elseif ":100-" -set_lev "1-2.89 m below ground" -endif -grib soil_temp.grib2

2 » Fix soil wv
wgrib2 temp1.grb2 -set_var SOILW -grib temp2.grb2

wgrib2 temp2.grb2 -if ":0-7 m below ground" -set_lev "0-0.07 m below ground" -elseif ":7-28 m below ground" -set_lev "0.07-0.28 m below ground" -elseif ":28-100 m below ground" -set_lev "0.28-1 m below ground" -elseif ":100-" -set_lev "1-2.89 m below ground" -endif -grib soil_wv.grib2


3 » Run ungrib.exe separately for vars_except_soil.grib2, soil_temp.grib2 and soil_wv.grib2: runs fine
4 » Run metgrid.exe with the 3 ungribbed files: runs fine

BUT, the resulting met_em files has wrong values for many variables... (attached) . Any idea why? You get "normal" met_em fields?
 

Attachments

  • met_em.d01.2024-06-03_00:00:00.nc.zip
    1.5 MB · Views: 8
I'm having the same problem now, please, did anyone succeed? Regarding this Vtable file, can you send me a copy, the data I want to process is also 0p25 ifs data.
 
You must log in or register to reply here.
Top