WRF driven by ERA-5 model level

[MKG]

Hi,

I am trying to run WRF simulation by using ERA-5 model level to simulate the upper stratosphere and lower mesosphere.
I download the ERA-5 model data from CDS, and I successfully run ungrib, metgrid, and real.
However, my simulation created nan values after a few steps and crashed. Note that the model top is set to 1 hPa.

Could you give me some suggestions?
I am using WRF Model Version 4.7.1 and WPS 4.6.0.
I downloaded values as follows:
3D data-> geopotential, T, u, V, q, lnsp,
Surface data-> mean_sea_level_pressure, skin_temperature, 2m_temperature, 10m_u_component_of_wind, 10m_v_component_of_wind, 2m_dewpoint_temperature, geopotential, land_sea_mask, sea_surface_temperature, sea_ice_cover, snow_depth, soil_temperature_level_1, soil_temperature_level_2, soil_temperature_level_3, soil_temperature_level_4, volumetric_soil_water_layer_1, volumetric_soil_water_layer_2, volumetric_soil_water_layer_3, volumetric_soil_water_layer_4, snow_density, surface_pressure

My vtable is as follows:
GRIB | Level| Level| Level| metgrid | metgrid | metgrid |GRIB2|GRIB2|GRIB2|GRIB2|
Code | Code | 1 | 2 | Name | Units | Description |Discp|Catgy|Param|Level|
-----+------+------+------+----------+----------+------------------------------------------+-----------------------+
130 | 109 | * | | TT | K | Temperature | 0 | 0 | 0 | 105 |
131 | 109 | * | | UU | m s-1 | U | 0 | 2 | 2 | 105 |
132 | 109 | * | | VV | m s-1 | V | 0 | 2 | 3 | 105 |
133 | 109 | * | | SPECHUMD | kg kg-1 | Specific humidity | 0 | 1 | 0 | 105 |
157 | 109 | * | | RHUM | % | Relative humidity | 0 | 1 | 1 | 105 |
129 | 1 | 0 | | SOILGEO | m2 s-2 | | 0 | 0 | | 103 |
| 1 | 0 | | SOILHGT | m | Terrain field of source analysis | 0 | 0 | | 106 |
165 | 1 | 0 | | UU | m s-1 | U At 10 m | 0 | 2 | 2 | 103 |
166 | 1 | 0 | | VV | m s-1 | V At 10 m | 0 | 2 | 3 | 103 |
167 | 1 | 0 | | TT | K | Temperature At 2 m | 0 | 0 | 0 | 103 |
168 | 1 | 0 | | DEWPT | K | At 2 m | 0 | 0 | 6 | 103 |
| 1 | 0 | | RH | % | Relative Humidity at 2 m | 0 | 0 | | 103 |
# 172 | 1 | 0 | | LANDSEA | 0/1 Flag | Land/Sea flag | 2 | 0 | 0 | 1 |
134 | 1 | 0 | | PSFC | Pa | Surface Pressure | 0 | 3 | 0 | 1 |
134 | 109 | 1 | | PSFCH | Pa | | 0 | 3 | 0 | 1 |
151 | 1 | 0 | | PMSL | Pa | Sea-level Pressure | 0 | 3 | 0 | 101 |
235 | 1 | 0 | | SKINTEMP | K | Sea-Surface Temperature | 0 | 3 | | 101 |
31 | 1 | 0 | | SEAICE | fraction | Sea-Ice-Fraction | 10 | 2 | 0 | 1 |
34 | 1 | 0 | | SST | K | Sea-Surface Temperature | 0 | 3 | | 101 |
141 | 1 | 0 | | SNOW_EC | m | | 0 | 3 | | 101 |
| 1 | 0 | | SNOW | kg m-2 |Water Equivalent of Accumulated Snow Depth| 0 | 3 | | 101 |
139 | 112 | 0 | 7 | ST000007 | K | T of 0-7 cm ground layer | 2 | 0 | 2 | 106 |
170 | 112 | 7 | 28 | ST007028 | K | T of 7-28 cm ground layer | 192 | 128 | 170 | 106 |
183 | 112 | 28 | 100 | ST028100 | K | T of 28-100 cm ground layer | 192 | 128 | 183 | 106 |
236 | 112 | 100 | 255 | ST100289 | K | T of 100-289 cm ground layer | 192 | 128 | 236 | 106 |
39 | 112 | 0 | 7 | SM000007 | m3 m-3 | Soil moisture of 0-7 cm ground layer | 192 | 128 | 39 | 106 |
40 | 112 | 7 | 28 | SM007028 | m3 m-3 | Soil moisture of 7-28 cm ground layer | 192 | 128 | 40 | 106 |
41 | 112 | 28 | 100 | SM028100 | m3 m-3 | Soil moisture of 28-100 cm ground layer | 192 | 128 | 41 | 106 |
42 | 112 | 100 | 255 | SM100289 | fraction | Soil moisture of 100-289 cm ground layer | 192 | 128 | 42 | 106 |
152 | 109 | * | | LOGSFP | Pa | Log surface pressure | 0 | 3 | 25 | 105 |
129 | 105 | * | | GEOPT | m2 s-2 | | 0 | 3 | 4 | 105 |
156 | 105 | * | | HGT | m | Height | 0 | 3 | 5 | 105 |
33 | 1 | 0 | | SNOW_DEN | kg m-3 | | 0 | 1 | 61 | 1 |
| 1 | 0 | | SNOWH | m | Physical Snow Depth | 0 | 1 | 11 | 1 |
-----+------+------+------+----------+----------+------------------------------------------+
#
# For use with ERA-interim model-level output.
#
# Grib codes are from Table 128
# http://www.ecmwf.int/services/archive/d/parameters/order=grib_parameter/table=128/
#
# snow depth is converted to the proper units in rrpr.F
#
# For ERA-interim data at NCAR, use the ml (sc and uv) and sfc sc files.


Also, namelists are attached.

Thank you in advance

 

[William.Hatheway]

Using the ECMWF Climate Data Store (CDS) to Download ERA5 Data in a Conda Environment

Using the ECMWF Climate Data Store (CDS) to Download ERA5 Data in a Conda Environment This guide shows you how to: Create a Conda environment called CDS with pip and cdsapi Activate the environment Set up your CDS API token (.cdsapirc) Use the Python script Download_CDS_ERA5_Data_Hourly.py to...

forum.mmm.ucar.edu

This might help

 

[Ming Chen]

Please download the utility we provide , then follow the instruction to process ERA5 model level data. This utility is specifically developed for ERA5 data processing. It works as expected and we have well tested it.

Let me know if you have any issues during the process.