real(kind=kind_phys), parameter physcons::con_epsqs =1.0E-10_kind_phys