real(kind=kind_phys), parameter physcons::rhoair = 1.28_kind_phys