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

ERA5 Downscaling

camtono

New member
Hello.

I'm an undergrad taking meteorology. For my research, I need to downscale my ERA5 data to probably around 3km. I was informed by an alumna who tried to downscale their ERA5 data, however, they mentioned that WRF would not accept the input file of ERA5 because it lacked mandatory meteorological variables. May I ask what the variables are?

If there is also an ERA5 downscaling guide in WRF, I'd highly appreciate it if you could leave a link to it.

Thank you in advance!
 
Good day.

I ran WPS today and was able to successfully ungrib my ERA5 from CDS data. geogrid was also a success. However, I encountered problems with metgrid.exe. When I ran metgrid, it gave this error:

Processing domain 1 of 1
Processing 1951-01-01_00
FILE
WARNING: Field PRES has missing values at level 200100 at (i,j)=(1,1)
WARNING: Field PMSL has missing values at level 200100 at (i,j)=(1,1)
WARNING: Field PSFC has missing values at level 200100 at (i,j)=(1,1)
WARNING: Field SOILHGT has missing values at level 200100 at (i,j)=(1,1)
ERROR: Missing values encountered in interpolated fields. Stopping.
--------------------------------------------------------------------------
MPI_ABORT was invoked on rank 0 in communicator MPI_COMM_WORLD
with errorcode 22032.

NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
You may or may not see output from other processes, depending on
exactly when Open MPI kills them.
--------------------------------------------------------------------------


The ERA5 data that I downloaded from CDS does not have Soil Height. However, I am sure that I chose to download mean sea level pressure (PMSL) and surface pressure (PSFC). I am sure that PRES means pressure, but I am unsure of what pressure. Anyway, I will have here attached several files that could hopefully help in solving this error.

Thanks in advance.
 

Attachments

  • metgrid.log
    15.9 KB · Views: 0
  • ungrib.log
    668.7 KB · Views: 0
  • namelist.wps
    614 bytes · Views: 0
@camtono

These are the variables to download for the ERA5 to work from CDS.

Pressure level:
'divergence', 'fraction_of_cloud_cover', 'geopotential',
'ozone_mass_mixing_ratio', 'potential_vorticity', 'relative_humidity',
'specific_cloud_ice_water_content', 'specific_cloud_liquid_water_content',
'specific_humidity', 'specific_rain_water_content', 'specific_snow_water_content',
'temperature', 'u_component_of_wind', 'v_component_of_wind', 'vertical_velocity',
'vorticity',

'divergence', 'fraction_of_cloud_cover', 'geopotential',
'ozone_mass_mixing_ratio', 'potential_vorticity', 'relative_humidity',
'specific_cloud_ice_water_content', 'specific_cloud_liquid_water_content',
'specific_humidity', 'specific_rain_water_content', 'specific_snow_water_content',
'temperature', 'u_component_of_wind', 'v_component_of_wind', 'vertical_velocity',
'vorticity',


Surface:

'10m_u_component_of_wind', '10m_v_component_of_wind', '2m_dewpoint_temperature',
'2m_temperature', 'land_sea_mask', 'mean_sea_level_pressure',
'sea_ice_cover', 'sea_surface_temperature', 'skin_temperature',
'snow_density', 'snow_depth', 'soil_temperature_level_1',
'soil_temperature_level_2', 'soil_temperature_level_3', 'soil_temperature_level_4',
'surface_pressure', 'volumetric_soil_water_layer_1', 'volumetric_soil_water_layer_2',
'volumetric_soil_water_layer_3', 'volumetric_soil_water_layer_4',
 
Hello,

I downloaded ERA5 data with the mandatory variables. I ran WPS, but I encountered the same problem as I did last time which was:

-----------------------------------------------------

WARNING: Field PRES has missing values at level 200100 at (i,j)=(1,1)
WARNING: Field PMSL has missing values at level 200100 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.
-----------------------------------------------------

I thought that maybe the Vtable I was using was wrong, so I changed it from Vtable.ECMWF to Vtable.ERA5-Interim.pl. However, the same problem occurs.
I also followed what @kwerner suggested where I run metgrid alongside the invariant file they provided here. But, when I try running metgrid with the invariant file, it says the following:

Processing domain 1 of 1
ERA5_INVARIANT
Processing 1951-01-01_00
FILE
ERROR: Cannot combine time-independent data with time-dependent data for field LANDSEA.mask

I am currently thinking that maybe the reason why it is not working is because I only used the 'Surface' ERA5 Data and not the 'Pressure level' one. However, I find it redundant using the data for pressure levels since I only need surface data. Hoping someone can help with this.
 

Attachments

  • metgrid.log
    5.4 KB · Views: 0
  • namelist.wps
    649 bytes · Views: 0
  • ungrib.log
    668.7 KB · Views: 0
@camtono,
Are you downloading the full global dataset, or are you only grabbing a subset of data? Often we see this error with a subset of data because the projection can make it tough to tell exactly where the appropriate edges of the domain should be. It's recommended to just download the full global set of data.
 
@camtono,
Are you downloading the full global dataset, or are you only grabbing a subset of data? Often we see this error with a subset of data because the projection can make it tough to tell exactly where the appropriate edges of the domain should be. It's recommended to just download the full global set of data.
I see. I do indeed only have a subset of the data. I will do this later after class and update here again.

Thanks again.
 
Hello again.

I was finally able to run metgrid thanks to @kwerner, but I am currently facing a problem with running real.exe. It runs into an error where the p_top_requested < grid%p_top possible from data. I assume that this means my p_top_requested in my namelist.input is less than the possible pressure top that is available in the met_em* files.

Where do I find this information? I tried running ncdump -h on one of the met_em* files, but it did not show any information on pressure top whatsoever. I also tried putting in 100000 Pa since this ERA5 data was for the surface only but it still ran into the error. I will be attaching files here that would hopefully help in solving this matter.
 

Attachments

  • met_em.d01.1951-12-31_00:00:00.nc.txt
    1.2 MB · Views: 0
  • namelist.input
    3.5 KB · Views: 0
  • rsl.error.0000
    525.7 KB · Views: 1
  • rsl.out.0000
    544.6 KB · Views: 0
I'm glad you were able to get past the WPS errors. In an effort to keep each thread to a specific issue, we ask that you repost your new question as a new thread - this time in the WRF section of the forum. We will address it there. This method makes searching and reading through threads easier for other users in the future. Thanks!
 
Top