Hi everyone,
I'm trying to run WRF 4.7.1 with WPS 4.6.0 using ERA5 data for a small domain over the UK (Kelmarsh area), but metgrid keeps crashing with missing values at grid point (1,1).
The error I get is:
WARNING: Field PRES has missing values at level 200100 at (i,j)=(1,1)
WARNING: Field GEOPT has missing values at level 100000 at (i,j)=(1,1)
WARNING: Field PSFC has missing values at level 200100 at (i,j)=(1,1)
ERROR: Missing values encountered in interpolated fields. Stopping.
What's puzzling me is that when I check the intermediate files with rd_intermediate.exe, the data looks perfectly fine at (1,1). For example, PSFC shows 99653.69 Pa which seems reasonable for surface pressure. All the other fields also have sensible values.
I've tried different Vtables (ECMWF, ERA-interim.pl, even created a custom one) but keep hitting the same wall. geogrid runs fine, ungrib completes successfully and generates the ERA5:* files, but metgrid just won't cooperate.
My domain setup is pretty standard - 41x41 grid points, 9km resolution, Lambert projection centered at 52.399°N, -0.942°W. I downloaded the ERA5 data from CDS with both surface fields (2m temp, 10m winds, surface pressure, etc.) and pressure levels from 1000-100 hPa.
Has anyone else run into this with ERA5 data? I'm wondering if it's something specific about the grid boundaries or if there's some trick with ERA5 that I'm missing. The data itself seems solid, so I'm thinking it might be an interpolation issue during the metgrid step.
Any ideas would be really appreciated. Thanks in advance!
