Hi guys,

I used that nice script Compute geopotential on model levels to compute the geopotential height on model levels of some ERA5 data. In doing so, I came across some strange behavior: At those locations on Earth where the orography has negative geopotential height values, the geopotential of the model levels decreases (gets more negative) when going upward in the atmosphere (from model level 137 to 0). Did anybody else came across that? Does anybody know how to fix that?

Thanks,

Amelie