real(kind=kind_phys), parameter physcons::con_rgas =8.314472_kind_phys