real(kind=kind_phys), parameter physcons::con_amch4 =16.043_kind_phys