def Weir(Q, C, b, p):
    yA = (Q / (C * b  * 2 / 3 * (2 * g / 3) ** 0.5 )) ** (2/3) + p
    return yA





Yupstream2 = EulerMethod(100, LC, b, m, k, S0, Qupstream / 2, Weir(1.5, CA, bA, pA))

def crest(Q, C, b, y):
    h = y - (Q / (C * b * 2 / 3 * (2 * g / 3) ** 0.5)) ** (2/3)
    return h

print('The iterations of the height of the crest at C is', crest(1.5, CC, bC, Yupstream2))