real(kind=kind_phys), parameter physcons::qamin = 1.e-16_kind_phys