def extrapolate_wind_direction(hgt, vr, ref_wind_direction, ref_hgt):
    """Extrapolate wind direction

    Extrapolate wind direction using the linear law and veer

    Parameters
    ----------
    hgt : float
        Target height
    vr : float
        Vertical wind veer
    ref_wind_direction : float
        Reference wind direction
    ref_hgt : float
        Reference height

    Returns
    -------
    float
        Wind direction at target height
    """
    return ref_wind_direction + vr * (hgt - ref_hgt)
