def METHOD_NAME(x, width, centre, steepness):
    """
    BOUT++ TanhHat function
    """
    return 0.5 * (
        np.tanh(steepness * (x - (centre - width / 2.0)))
        + np.tanh(-steepness * (x - (centre + width / 2.0)))
    )