I have calculated the ∂θe/∂z with 'Height', 'Latitude', 'Longitude' dimensions.
import numpy as np import xarray as xr from wrf import (getvar) eth = np.array(getvar(ncfile, "eth", timeidx=time)) z = np.array(getvar(ncfile, "z", timeidx=time, units="km")) deth = np.gradient(eth, axis=0) dz = np.gradient(z, axis=0) deth_dz = deth/dz dthdz = xr.DataArray(data=deth_dz, dims=['Height', 'Latitude', 'Longitude'])
I would like to plot the horizontal depth or thickness (in km) of the
It is a bit similar to this question. But when I do
dthdz.sum(dim='Height'), I get very low depth values. I just want to know where I'm doing it wrong. Thank you!