Hello everyone,
I'm currently working with WRF v4.1.1 and WPS v4.1, and I've encountered the following fatal error when running real.exe:
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE: <stdin> LINE: 5857
troubles, could not find trapping x locations
-------------------------------------------
My data sources:
Geogrid data: https://www2.mmm.ucar.edu/wrf/src/wps_files/geog_high_res_mandatory.tar.gz
Meteorological data: ECMWF Operational Model Analysis (NCAR RDA Dataset d113001), 6 hourly atmospheric "surface" and "isobaric" analysis, (subsets of Latitudes: -28/-38 & Longitudes: -77/-65)
"25 isobaric pressure levels 1.0, 0.95, 0.925, 0.9, 0.85, 0.8, 0.7, 0.6, 0.5, 0.4, 0.3, 0.25, 0.2, 0.15, 0.1, 0.07, 0.05, 0.03, 0.02, 0.01, 0.007, 0.005, 0.003, 0.002, and 0.001"
I converted the ECMWF GRIB1 files (Gaussian grid) to a regular lat-lon grid using cdo remapnn. Then, ungrib.exe and metgrid.exe both completed successfully. I think that all essential parameters have been accounted for. However, when running real.exe, I see multiple warnings such as [ t(i,j,k) was 0 at 1 1 1 , setting Qv to 0] before encountering the fatal error mentioned above. Based on various discussions and suggestions I found in the forum, I've tried the followings modifications (regardless if they were matching my conditions precisely or not, I just wanted to try):
=============================
&domains
# e_vert = different values > 30 up to 50
auto_levels_opt = 2,
dzbot = 50
dzstretch_s = 1.3,
dzstretch_u = 1.1,
max_dz = 10000,
=============================
&domains
# p_top_requested = 5000 or 10000
=============================
&domains
e_vert = 26
eta_levels = 1.0, 0.95, 0.925, 0.9, 0.85,
0.8, 0.7, 0.6, 0.5, 0.4,
0.3, 0.25, 0.2, 0.15, 0.1,
0.07, 0.05, 0.03, 0.02, 0.01,
0.007, 0.005, 0.003, 0.002, 0.001,
0.0,
=============================
&domains
sfcp_to_sfcp = .true.
=============================
&domains
lowest_lev_from_sfc = .false.
force_sfc_in_vinterp = 0
use_surface = .false.
=============================
&physics
Changed physics_suite to CONUS in order to isolate potential issues related to physics configuration.
=============================
In the METGRID.TBL
modifying the intrepolation options (removing the sixteen terms)
=============================
Despite all these modifications, the error still persists. If anyone has experienced a similar issue or has insights into what else I might try, I would be truly grateful for your help.
Best regards,
Mohsen
Sample intermediate file by ungrib.exe: FILEïº2024-04-30_00
Sample met_em.nc file by metgrid.exe: met_em.d01.2024-04-30_00ïº00ïº00.nc
I'm currently working with WRF v4.1.1 and WPS v4.1, and I've encountered the following fatal error when running real.exe:
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE: <stdin> LINE: 5857
troubles, could not find trapping x locations
-------------------------------------------
My data sources:
Geogrid data: https://www2.mmm.ucar.edu/wrf/src/wps_files/geog_high_res_mandatory.tar.gz
Meteorological data: ECMWF Operational Model Analysis (NCAR RDA Dataset d113001), 6 hourly atmospheric "surface" and "isobaric" analysis, (subsets of Latitudes: -28/-38 & Longitudes: -77/-65)
"25 isobaric pressure levels 1.0, 0.95, 0.925, 0.9, 0.85, 0.8, 0.7, 0.6, 0.5, 0.4, 0.3, 0.25, 0.2, 0.15, 0.1, 0.07, 0.05, 0.03, 0.02, 0.01, 0.007, 0.005, 0.003, 0.002, and 0.001"
I converted the ECMWF GRIB1 files (Gaussian grid) to a regular lat-lon grid using cdo remapnn. Then, ungrib.exe and metgrid.exe both completed successfully. I think that all essential parameters have been accounted for. However, when running real.exe, I see multiple warnings such as [ t(i,j,k) was 0 at 1 1 1 , setting Qv to 0] before encountering the fatal error mentioned above. Based on various discussions and suggestions I found in the forum, I've tried the followings modifications (regardless if they were matching my conditions precisely or not, I just wanted to try):
=============================
&domains
# e_vert = different values > 30 up to 50
auto_levels_opt = 2,
dzbot = 50
dzstretch_s = 1.3,
dzstretch_u = 1.1,
max_dz = 10000,
=============================
&domains
# p_top_requested = 5000 or 10000
=============================
&domains
e_vert = 26
eta_levels = 1.0, 0.95, 0.925, 0.9, 0.85,
0.8, 0.7, 0.6, 0.5, 0.4,
0.3, 0.25, 0.2, 0.15, 0.1,
0.07, 0.05, 0.03, 0.02, 0.01,
0.007, 0.005, 0.003, 0.002, 0.001,
0.0,
=============================
&domains
sfcp_to_sfcp = .true.
=============================
&domains
lowest_lev_from_sfc = .false.
force_sfc_in_vinterp = 0
use_surface = .false.
=============================
&physics
Changed physics_suite to CONUS in order to isolate potential issues related to physics configuration.
=============================
In the METGRID.TBL
modifying the intrepolation options (removing the sixteen terms)
=============================
Despite all these modifications, the error still persists. If anyone has experienced a similar issue or has insights into what else I might try, I would be truly grateful for your help.
Best regards,
Mohsen
Sample intermediate file by ungrib.exe: FILEïº2024-04-30_00
Sample met_em.nc file by metgrid.exe: met_em.d01.2024-04-30_00ïº00ïº00.nc