def METHOD_NAME(par):
    """Create small dataset with several linear events and check output
    contains the events at correct time and correct amplitude
    """
    # Data creation
    v = 1
    t0 = (50, 130)
    theta = (0.0, 0.0)
    phi = (0.0, 0.0)
    amp = (0.6, 1)
    # Create axes
    t, _, x, y = makeaxis(par)
    # Create data
    d, dwav = linear3d(x, y, t, v, t0, theta, phi, amp, wav)
    # Assert shape
    assert d.shape[0] == par["ny"]
    assert d.shape[1] == par["nx"]
    assert d.shape[2] == par["nt"]
    assert dwav.shape[0] == par["ny"]
    assert dwav.shape[1] == par["nx"]
    assert dwav.shape[2] == par["nt"]
    # Assert correct position of event
    assert_array_equal(d[:, :, t0[0]], amp[0] * np.ones((par["ny"], par["nx"])))
    assert_array_equal(d[:, :, t0[1]], amp[1] * np.ones((par["ny"], par["nx"])))