from pyatmos import coesa76
coesa76_geom = coesa76([0,20,40,60,80]) # geometric altitudes by default
print(coesa76_geom.rho) # [kg/m^3]
print(coesa76_geom.T) # [K]
print(coesa76_geom.P) # [Pa]
# coesa76_geop = coesa76([0,20,40,60,80],'geopotential') # geopotential altitudes