def AeroReist(um, zm, z0, d, zmp=zm):
    k = 0.4
    r_a = 1.0/(k**2*um)*log((zm - d)/z0)*log((d-zmp)/(z0/10.0))
    return r_a
dataset/ETHPy150Open LEAF-BoiseState/SPEED/Module05/PenmanMonteithEx.py/AeroReist/swapped operands of `-`