...
So Mark wrote his own script and kindly provided it to us. The script takes temperature and specific humidity (t and q) on model levels as inputs, along with geopotential and the pressure (z and lnsp) on the surface, and it creates as output calculates the geopotential in m^2/s^2 for each model level. It then also calculates the height in meters by dividing the geopotential by the gravity of Earth (9.80665 m/s^2).
...