calculation of height from geopotential

[met-sree]

Dear all,
I am trying to convert geopotential data (era5 data with units: m2/s2) into height (m) by using metpy. But, it showing errors like the above attached. How to rectify this error. please suggest me.

Thank you,
Srinivas.

[dcamron]

Hello! Can you edit your post or reply here with some of your original code and the full error you received? A screenshot or copy-and-paste will do. Thanks!

[met-sree]

Thank you for your reply... I added an error screenshot.

[dcamron]

I see you have a commented out line constructing a Quantity. You will need to create one with associated units for your geop so that geopotential_to_height() can understand the units of your values, e.g.

geop = gpdata['z'][:]
geop_w_units = units.Quantity(geop, 'm ** 2 / s ** 2')
# geop_w_units = geop * units('m ** 2 / s ** 2') also can work most of the time
hgt = metpy.calc.geopotential_to_height(geop_w_units)

We have a more helpful error to catch this and are not using it here! So thanks for catching that 😄