real, parameter radcons::prsmin = 1.0e-6