real(kind=kind_phys), parameter physcons::con_omega =7.2921e-5_kind_phys