def METHOD_NAME():
    filename = os.path.join("tests", "testdata", "HD2D.dfsu")
    dfs = mikeio.open(filename)
    da = dfs.read(items="Surface elevation", time=0)[0]
    elem40 = np.arange(40)
    ax = da.plot(elements=elem40) # this is ok
    assert ax is not None
    wl_40 = da.values[elem40]
    with pytest.warns(FutureWarning):
        dfs.plot(wl_40, elements=elem40, plot_type="shaded", levels=5)