def ConcentrationOxonium (c, Ka=1.0e-3):
    return (-Ka+sqrt(Ka**2+4*c*Ka))/2.

def p(x):
    return -log10(x)