[Wrf-users] Screen level variables not ungribbing in ungrib.exe from custom built Vtable
Elisabeth Donnell
edonnell at weatherzone.com.au
Wed Sep 8 19:45:32 MDT 2010
Hello,
I am trying to ultimately drive WRF from the Australian Bureau of Meteorology new NWP model ACCESS-G which is based on the UKMO unified model.
I have created a custom Vtable with help from utilities g1print and wgrib.
The problem is that the variables at screen height are not ungribbing with ungrib.exe
Is it possible my troubles are caused by the level type(LEVEL TYPE =125) for these variables not being handled by ungrib?. (This is height level above ground in centimeters.) If I put the level1=150cm (as specified by g1print) then ungrib crashes with the following error :
Could not find a level for this Vtable entry
iparm = 167 ktype = 125 lvl1 = 150 lvl2 =
-99
Fix the Vtable or modify rd_grib1.F
I have tried changing the LEVEL TYPE = 105 (and lvl1=2), and while ungrib.exe runs successfully the screen height variables are not ungribbed, similarly for ktype =125 and lvl1=200.
For the cases where ungrib.exe runs without error (even though the screen level variables are not ungribbed) I get an error running metgrid because I think it requires these fields (??)
WRF_DEBUG: Warning DIM 4 , NAME num_metgrid_levels REDIFINED by var TT 6 5 in wrf_io.F90
line 2424
ERROR: Error in ext_pkg_write_field
Any suggestions/help would be most welcome!
Below is the Vtable I use (this version runs but as I said I have tried different Level1 values and Types):
GRIB1| Level| From | To | metgrid | metgrid | metgrid |
Param| Type |Level1|Level2| Name | Units | Description |
-----+------+------+------+----------+---------+------------------------------------------+
130 | 100 | * | | TT | K | Temperature |
131 | 100 | * | | UU | m s-1 | U |
132 | 100 | * | | VV | m s-1 | V |
17 | 100 | * | | DEWPT | K | |
156 | 100 | * | | HGT | m | Geopotential Height
167 | 125 | 200 | | TT | K | Temperature at 2 m |
168 | 125 | 200 | | DEWPT | K | Dew point temperature at 2 m |
101 | 125 | 200 | | RH | % | Relative Humidity at 2 m |
165 | 105 | 10 | | UU | m s-1 | U at 10 m |
166 | 105 | 10 | | VV | m s-1 | V at 10 m |
134 | 1 | 0 | | PSFC | Pa | Surface Pressure |
151 | 102 | 0 | | PMSL | Pa | Sea-level Pressure |
125 | 1 | 0 | | SKINTEMP | K | Skin Temperature |
159 | 1 | 0 | | SOILHGT | m | Topography Elevation |
172 | 1 | 0 | | LANDSEA | proprtn | Land/Sea flag (1=land,0=sea) |
140 | 112 | 0 | 10 | SM000010 | kg m-3 | Soil Moist 0-10 cm below grn layer (Up) |
-----+------+------+------+----------+---------+------------------------------------------+
This is the output from g1print.exe from an ACCESS-G grib file.
----------------------------------------------------
rec GRIB GRIB Lvl Lvl Lvl Time Fcst
Num Code name Code one two hour
----------------------------------------------------
1 134 1 0 0 2010-09-07_12:00 +00
2 167 125 150 0 2010-09-07_12:00 +00
3 168 125 150 0 2010-09-07_12:00 +00
4 101 125 150 0 2010-09-07_12:00 +00
5 165 105 10 0 2010-09-07_12:00 +00
6 166 105 10 0 2010-09-07_12:00 +00
7 201 125 150 0 2010-09-07_12:00 +00
8 151 102 0 0 2010-09-07_12:00 +00
9 202 125 150 0 2010-09-07_12:00 +00
10 125 1 0 0 2010-09-07_12:00 +00
11 164 200 0 0 2010-09-07_12:00 +00
12 61 1 0 0 2010-09-07_12:00 +00
13 137 1 0 0 2010-09-07_12:00 +00
14 159 1 0 0 2010-09-07_12:00 +00
15 188 200 0 0 2010-09-07_12:00 +00
16 186 200 0 0 2010-09-07_12:00 +00
17 187 200 0 0 2010-09-07_12:00 +00
18 140 112 0 10 2010-09-07_12:00 +00
19 58 1 0 0 2010-09-07_12:00 +00
20 57 1 0 0 2010-09-07_12:00 +00
21 172 1 0 0 2010-09-07_12:00 +00
22 131 100 1000 0 2010-09-07_12:00 +00
23 132 100 1000 0 2010-09-07_12:00 +00
24 130 100 1000 0 2010-09-07_12:00 +00
25 134 100 1000 0 2010-09-07_12:00 +00
26 156 100 1000 0 2010-09-07_12:00 +00
27 17 100 1000 0 2010-09-07_12:00 +00
28 131 100 300 0 2010-09-07_12:00 +00
29 132 100 300 0 2010-09-07_12:00 +00
30 130 100 300 0 2010-09-07_12:00 +00
31 134 100 300 0 2010-09-07_12:00 +00
32 156 100 300 0 2010-09-07_12:00 +00
33 17 100 300 0 2010-09-07_12:00 +00
34 131 100 500 0 2010-09-07_12:00 +00
35 132 100 500 0 2010-09-07_12:00 +00
36 130 100 500 0 2010-09-07_12:00 +00
37 134 100 500 0 2010-09-07_12:00 +00
38 156 100 500 0 2010-09-07_12:00 +00
39 17 100 500 0 2010-09-07_12:00 +00
40 131 100 700 0 2010-09-07_12:00 +00
41 132 100 700 0 2010-09-07_12:00 +00
42 130 100 700 0 2010-09-07_12:00 +00
43 134 100 700 0 2010-09-07_12:00 +00
44 156 100 700 0 2010-09-07_12:00 +00
45 17 100 700 0 2010-09-07_12:00 +00
46 131 100 850 0 2010-09-07_12:00 +00
47 132 100 850 0 2010-09-07_12:00 +00
48 130 100 850 0 2010-09-07_12:00 +00
49 134 100 850 0 2010-09-07_12:00 +00
50 156 100 850 0 2010-09-07_12:00 +00
51 17 100 850 0 2010-09-07_12:00 +00
and an extract from ungrib.log
###############################################################################
Inventory for date = 2010-09-11 06:00:00
PRES TT UU VV DEWPT HGT RH PSFC PMSL SKINTEMP SOILHGT LANDSEA SM000010
-------------------------------------------------------------------------------
2013.0 O O O O O O O X O O O O
2001.0 O X X O O O X O X X O X
1000.0 X X X X X
850.0 X X X X X
700.0 X X X X X
500.0 X X X X X
300.0 X X X X X
-------------------------------------------------------------------------------
Regards,
Dr. Elisabeth Donnell | NWP Scientist
weatherzone°
M 0425322032
weatherzone.com.au
More information about the Wrf-users
mailing list