def METHOD_NAME(keep_original=False, **kwargs):
    """Perform velocity dealiasing on reference data."""
    radar = pyart.testing.make_velocity_aliased_radar()
    # speckling that will not be dealiased.
    radar.fields["velocity"]["data"][13, -4:] = [-7.5, 8.5, 0, 0]
    height = np.linspace(150, 250, 10).astype("float32")
    speed = np.ones((10), dtype="float32") * 0.5
    direction = np.ones((10), dtype="float32") * 5.0
    profile = pyart.core.HorizontalWindProfile(height, speed, direction)
    dealias_vel = pyart.correct.dealias_fourdd(
        radar, sonde_profile=profile, keep_original=keep_original, **kwargs
    )
    return radar, dealias_vel