def METHOD_NAME(x):
    y = (y0 - D[ii]) - (x - x0) * wake_slope
    if isinstance(y, (float, np.float64, np.float32)):
        if x < (x0 + 0.01):
            y = -np.Inf
    else:
        y[x < x0 + 0.01] = -np.Inf
    return y