Some problems in real.exe with the ECMWF

Hello,

I am running WRF 4.1.4 with the ECMWF boundary conditions. It has been successful until the real.exe. I obtain a 174 exit code and checking, for instance, the rsl.error.0000 I see lots of messages like this: Warning: vapor pressure exceeds total pressure, setting Qv to 1.E-6. It does not seem a good signal.

I have visualized the generated wrfinput_d01, and there are things looking weird. For instance, the TSK has huge oscillations, with values in the shores above 330 K and below 170 K. There seems to be also high values of T in the vertical, with 70 vertical levels, the level 70 has T above 450 K and the 0 level is below 100 K, which it is physically impossible.

I am running a namelist very similar to the namelist for 4 Km (I run with 2.5 Km in dx and dy). Do you have any idea of what is going wrong and what could I do? Thx.

What ECMWF data are you using? Please also send me your namelist.wps to take a look.

Hi, thanks for your answer. I am using the ECMWF data in grib 2, but I have included grib 2 information in the Vtable (ungrib.exe has a succesful execution). This is the namelist.wps you requested:

&share
wrf_core = 'ARW',
max_dom = 1,
start_date = '2020-06-29_00:00:00',
end_date = '2020-06-30_00:00:00',
interval_seconds = 10800,
io_form_geogrid = 2,
/

&geogrid
parent_id = 1,
parent_grid_ratio = 1,
i_parent_start = 1,
j_parent_start = 1,
e_we = 566,
e_sn = 470,
!
!!!!!!!!!!!!!!!!!!!!!!!!!!!! IMPORTANT NOTE !!!!!!!!!!!!!!!!!!!!!!!!!!!!
! The default datasets used to produce the MAXSNOALB and ALBEDO12M
! fields have changed in WPS v4.0. These fields are now interpolated
! from MODIS-based datasets.
!
! To match the output given by the default namelist.wps in WPS v3.9.1,
! the following setting for geog_data_res may be used:
!
! geog_data_res = 'maxsnowalb_ncep+albedo_ncep+default', 'maxsnowalb_ncep+albedo_ncep+default',
!
!!!!!!!!!!!!!!!!!!!!!!!!!!!! IMPORTANT NOTE !!!!!!!!!!!!!!!!!!!!!!!!!!!!
!
geog_data_res = 'modis_landuse_20class_30s_with_lakes',
dx = 2500,
dy = 2500,
map_proj = 'lambert',
ref_lat = 29.00,
ref_lon = -17.50,
truelat1 = 29.00,
truelat2 = 29.00,
stand_lon = -17.50,
geog_data_path = '/scratch/WRF4.0/WPS_GEOG',
/

&ungrib
out_format = 'WPS',
prefix = 'FILE',
/

&metgrid
fg_name = 'FILE'
io_form_metgrid = 2,
/

Let me add the Vtable I have been using, there are minor differences with the standard Vtable.ECMWF_sigma provided (losing format when pasting, try to fix it, sorry):

GRIB | Level| Level| Level| metgrid | metgrid | metgrid |GRIB2|GRIB2|GRIB2|GRIB2|
Code | Code | 1 | 2 | Name | Units | Description |Discp|Catgy|Param|Level|
-----+------+------+------+----------+----------+------------------------------------------+-----------------------+
130 | 109 | * | | TT | K | Temperature | 0 | 0 | 0 | 105 |
131 | 109 | * | | UU | m s-1 | U | 0 | 2 | 2 | 105 |
132 | 109 | * | | VV | m s-1 | V | 0 | 2 | 3 | 105 |
133 | 109 | * | | SPECHUMD | kg kg-1 | Specific humidity | 0 | 1 | 0 | 105 |
152 | 109 | * | | LOGSFP | Pa | Log surface pressure | 0 | 3 | 25 | 105 |
| 109 | 1 | | SOILGEO | m2 s-2 | Surface geopotential (Orography) | 0 | 3 | | 105 |
129 | 1 | 0 | | SOILHGT | m | Terrain field of source analysis | 0 | 3 | 5 | 1 |
165 | 1 | 0 | | UU | m s-1 | U at 10 m | 0 | 2 | 2 | 103 |
166 | 1 | 0 | | VV | m s-1 | V at 10 m | 0 | 2 | 3 | 103 |
167 | 1 | 0 | | TT | K | Temperature at 2 m | 0 | 0 | 0 | 103 |
168 | 1 | 0 | | DEWPT | K | Dew point temperature at 2 m | 0 | 0 | 6 | 103 |
| 1 | 0 | | RH | % | Relative Humidity at 2 m | 0 | 1 | | 103 |
172 | 1 | 0 | | LANDSEA | 0/1 Flag | Land/Sea flag | 2 | 0 | 0 | 1 |
134 | 1 | 0 | | PSFC | Pa | Surface Pressure | 0 | 3 | | 1 |
| 109 | 1 | | PSFCH | Pa | | 0 | 3 | | 1 |
151 | 1 | 0 | | PMSL | Pa | Sea-level Pressure | 0 | 3 | 0 | 101 |
235 | 1 | 0 | | SKINTEMP | K | Sea-Surface Temperature | 0 | 0 | 17 | 1 |
31 | 1 | 0 | | SEAICE | 0/1 Flag | Sea-Ice-Flag | 10 | 2 | 0 | 1 |
34 | 1 | 0 | | SST | K | Sea-Surface Temperature | 10 | 3 | 0 | 1 |
141 | 1 | 0 | | SNOWH | m | Snow depth in m of Water Equivalent | 0 | 1 | 11 | 1 |
65 | 1 | 0 | | SNOW | kg m-2 |Water Equivalent of Accumulated Snow Depth | 0 | 1 | 13 | 1 |
139 | 112 | 0 | 7 | ST000007 | K | T of 0-7 cm ground layer | 192 | 128 | 139 | 106 |
170 | 112 | 7 | 28 | ST007028 | K | T of 7-28 cm ground layer | 192 | 128 | 170 | 106 |
183 | 112 | 28 | 100 | ST028100 | K | T of 28-100 cm ground layer | 192 | 128 | 183 | 106 |
236 | 112 | 100 | 255 | ST100289 | K | T of 100-289 cm ground layer | 192 | 128 | 236 | 106 |
39 | 112 | 0 | 7 | SM000007 | m3 m-3 | Soil moisture of 0-7 cm ground layer | 192 | 128 | 39 | 106 |
40 | 112 | 7 | 28 | SM007028 | m3 m-3 | Soil moisture of 7-28 cm ground layer | 192 | 128 | 40 | 106 |
41 | 112 | 28 | 100 | SM028100 | m3 m-3 | Soil moisture of 28-100 cm ground layer | 192 | 128 | 41 | 106 |
42 | 112 | 100 | 255 | SM100289 | m3 m-3 | Soil moisture of 100-289 cm ground layer | 192 | 128 | 42 | 106 |
-----+------+------+------+----------+----------+------------------------------------------+-----------------------+

Looking with more detail the meesages of the rsl.error.0000, I have found this notification:

d01 2020-07-07_00:00:00 Fatal error BAD VARIABLE DIMENSION in ext_ncd_read_field RH2020-07-07_00:00:00
d01 2020-07-07_00:00:00 StoredDim 3 .NE. NDim+1 4
input_wrf.F reading 3d real rh_gc Status = -1004

It could be important. The first time I only noticed the message

d01 2020-07-07_00:00:00 Warning: vapor pressure exceeds total pressure, setting Qv to 1.E-6

that now I am thinking it could be related to that fatal error. Thx.

Hello, I just wanted to say that this problem was solved running the software "calc_ecmwf_p.exe" after the "ungrib.exe". "calc_ecmwf_p.exe" can be found under /WPS/util/