from pyatmos import expo
expo_geom = expo([0,20,40,60,80]) # geometric altitudes by default
print(expo_geom.rho) # [kg/m^3]
# expo_geop = expo([0,20,40,60,80],'geopotential') # geopotential altitudes