def melt_bond_length(p, coeffs):
    """
    Return the bond length (in m) given
    a polynomial expansion in pressure.
    """

    r = 0.0000000001 * (coeffs[0] + p * coeffs[1] + p**2 * coeffs[2])
    return r