3
$\begingroup$

I have calculated the $∂θ_\varepsilon/∂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 dthdz. 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!

Improve this question
$\endgroup$
4
  • $\begingroup$ For programming related questions unless you provide the dataset very few here can help you $\endgroup$
    – gansub
    Sep 14, 2021 at 7:44
  • $\begingroup$ Why not use MetPy's thickness API - unidata.github.io/MetPy/latest/api/generated/… ? $\endgroup$
    – gansub
    Sep 14, 2021 at 7:45
  • $\begingroup$ thanks @gansub! this uses conventional wrfout data. i tried the thickness_hydrostatic but it outputs in 1D. i was hoping to make a horizontal plot of the depth $\endgroup$
    – Kas
    Sep 14, 2021 at 8:36
  • $\begingroup$ Suggest you take this over to MetPy gitter and then ask the question over there. $\endgroup$
    – gansub
    Sep 14, 2021 at 8:49

0

Reset to default

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Browse other questions tagged or ask your own question.