WRF output in MetPY

[Armeteo]

I want to plot different variables of a .nc WRF model running output file. But when I write:

import cartopy.crs as ccrs
import cartopy.feature as cfeature
import matplotlib.pyplot as plt
from netCDF4 import Dataset, num2date

from metpy.cbook import get_test_data
from metpy.interpolate import log_interpolate_1d
from metpy.plots import add_metpy_logo, add_timestamp
from metpy.units import units

data = Dataset('/home/arnaitz/Build_WRF/WRF-4.1.5/run/wrfout_d02_2021-04-10_06:00:00')
lat = data.variables['lat'][:]
lon = data.variables['lon'][:]
time = data.variables['time']
vtimes = num2date(time[:], time.units)
temperature = units.Quantity(data.variables['temperature'][:], 'degC')
pres = units.Quantity(data.variables['pressure'][:], 'Pa')
hgt = units.Quantity(data.variables['height'][:], 'meter')

AN ERROR OCCUR:

KeyError Traceback (most recent call last)
in
1 data = Dataset('/home/arnaitz/Build_WRF/WRF-4.1.5/run/wrfout_d02_2021-04-10_06:00:00')
----> 2 lat = data.variables['lat'][:]
3 lon = data.variables['lon'][:]
4 time = data.variables['time']
5 vtimes = num2date(time[:], time.units)

KeyError: 'lat'

I think is because in the output file is not defined lat like in MetPy. Do you know how to fix this?

[stavrosdafis]

I suppose you have not changed the name of the default output variables, so the latitude should be 'XLAT' and longitude 'XLONG' and are 2D variables. Also, you don't have a 'time' variable but 'Times' and temperature, pressure, hgt all have different names from your script.

You can check the names of the WRF output variables with ncdump, e.g.:

ncdump -h wrfout_d02_2021-04-10_06:00:00