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

Trouble with ungribbing ERA5 data - Mismatch between namelist and input file dimensions

Pechudin

Member
Hello,

After a year-long hiatus, I have started running WRF simulations again. Last time I used WRF, I downloaded data from the Copernicus CDS, and ungribbed it using the Vtable.ECMWF. I used the cdsapi python script to download the data. All worked fine. The WRF version was 4.3.3, and I ran it on a supercomputer.

I once again started working on the supercomputer, again using WRF 4.3.3 and the same Vtable. I downloaded the pressure-level and single-level data (same data as previously, although the python API seems to have changed somewhat in the meantime), one file for 3D data (pressure level) and one for the 2D data (single-level data). I ran ungrib and metgrid, producing met_em files.

However, running real.exe produces the following:
Input data is acceptable to use: met_em.d01.2023-09-22_00:00:00.nc
metgrid input_wrf.F first_date_input = 2023-09-22_00:00:00
metgrid input_wrf.F first_date_nml = 2023-09-21_00:00:00
input_wrf.F: SIZE MISMATCH: namelist num_metgrid_levels = 38
input_wrf.F: SIZE MISMATCH: input file BOTTOM-TOP_GRID_DIMENSION = 1
---- ERROR: Mismatch between namelist and input file dimensions
NOTE: 1 namelist vs input data inconsistencies found.
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE: <stdin> LINE: 1301
NOTE: Please check and reset these options
-------------------------------------------
MPICH ERROR [Rank 0] [job id ] [Tue Jun 24 09:06:35 2025] [x8000c0s0b0n0] - Abort(1) (rank 0 in comm 0): application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0

What seems to have happened is that some of the met_em files have 1/3 of the size of others. More precisely, I took data from 21.09.2023. to 26.09.2023., 3-hourly spacing, and for the data for 22.09.2023., metgrid put the UU, VV, GHT, RH and TT fields into the 2D category, meaning only one level was read??

It seems to me that the issue is a change in how ERA5 data is interacting with WPS, as using data downloaded previously (in 2023) with WRF yields the correct results.

So, my question is - what changed and what is the correct way to use ERA5 data with WRF nowadays?

I will post the cdsapi requests to this thread shortly.
 

Attachments

  • ERA5_request_2D.txt
    1.2 KB · Views: 0
  • ERA5_request_3D.txt
    1.1 KB · Views: 0
I figured out what the issue was. I had to ungrib the data twice; once for the pressure-level data, and once for the surface-level data, and then use metgrid. This thread can thus be closed, if that is necessary.
 
Top