Missing values encountered in interpolated fields when running metgrid with UCAR ERA5 Data

[fcarotenuto]

Dear all,

I'm (desperately) try to run WRF with the ERA5 data available from UCAR through the era5_to_int Python script.
Through extensive research and trial and error I think I've understood the variables needed to run WRF with the --isobaric option of era5_to_int which are the following:

- Dataset "d633000". On pressure levels (mbar): Z/Q/T/U/V
- Dataset "d633000". Surface variables: SOILGEO/MSL/VAR_2T/VAR_2D/VAR_10U/VAR_10V/RSN/SD/SSTK/SKT/SWVL1/SWVL2/SWVL3/SWVL4/STL1/STL2/STL3/STL4/CI
- Dataset "d633006". On model levels: SP
- Invariants: LSM (d633000, 1979/01/01 00:00 - 1979/01/01 00:00), Z (d633006, 2016/01/01 00:00 - 2016/01/01 00:00)

This set of variables allows era5_to_int to run as era5_to_int.py -i -p /home/wrf/Software/Data/ERA5 2023-07-01_00 2023-07-15_23 1 &> era5_to_int.out and produce a complete set of "ERA5:" data (see era5_to_int.out.txt attached).

Still, If I try to run metgrid.exe it complaints about encountering missing values in interpolated fields (PRES, SOILHGT, PSFC, PMSL, GHT) and does not run (see attached metgrid.out.txt and metgrid.log).

The things that I've tried to do up until now were:

1) Increase the domain over which ERA5 data were subsetted, to be sure that my WRF domain was more than fully included. Error is the same.

2) Use model levels data instead than isobaric ones. No matter how I do the request it downloads only the "Q" data and even if the request is from 00:00 of 2023/07/01 the first file is 2023/07/01 06. Thus era5_to_int complains on missing 2023/07/01 00.

Any help would be super appreciated

Thanks!

 

[Ming Chen]

I just run a test case using era5_to_int.py to create intermediate files. All work fine and my case can run through WPS to WRF.

Would you please upload your namelist.wps for me to take a look?

Also, please tell me what files you download from NCAR RDA and your command usage.

Thanks.

 

[fcarotenuto]

Dear Ming Chen,

Thank you for your reply!

In attachment I've put:
1. my namelist.wps
2. the list of files I've downloaded from NCAR RDA
3. the ipynb I use to obtain said files

The command I use is:
era5_to_int.py -i -p /path/to/files 2023-07-01_00 2023-07-15_23 1 &> era5_to_int.out

When I run it for model levels files I remove "-i" from the command.

Thank you again for your support!

 

[Ming Chen]

I just noticed that you didn't download global ERA5 data? Please let me know if I am wrong.

Note that some of the data are on Gaussian grids, and our script can only handle global data.

Please download the entire ERA5 data on global grids and rerun the case. Let me know if you still have issues.