function rho = rhodryair(p,T)
% input: T [^o C], p [mbar]
% output: kg/m^3
constants

rho = 100*p / Rd ./ (T+273.15);
