Issue plotting contour levels

[tomadams-ssf]

I have encountered an issue with plotting contour levels. This occurs in both the version of the package available on pip and the development version on GitHub.

When using the “contourf” option with levels, the colour shading is incorrect. If I set e.g. “levels = [0,1,250, 2000]”, all values less than 1 are shown with the shading for the 1-250 category. I have verified this by checking values in the original file using MIKE Data Manager, and also by generating a DataArray class with zeros everywhere (as below) - in addition to other values less than 1, or even negative.

Code to create this behaviour is given below (uses DataArray class from development version). Note that the data (sedMean2) used to construct the DataArray below is 0 everywhere, but the map shows 1 everywhere. The same thing occurs when using the old “dfs.plot()” method.

I'm interested to understand what I’m missing – or whether there is a bug. I can't attach the source DFSU file here but presume the same behaviour would be seen using another domain, as I'm only using it to set the geometry.

Thanks
Tom

Example code:

dfs = Dfsu('meanDep.dfsu')
ds = dfs.read()
sedMean = ds[1].to_numpy() # Values not used for plot in present DataArray created
sedMean2 = np.zeros(np.shape(sedMean))

da = mikeio.DataArray(data=sedMean2,
geometry=ds.geometry,
item=ItemInfo('Mean sedimented mass',EUMType.Mass_per_Unit_Area,EUMUnit.gram_per_meter_pow_2))
ax = da.plot(plot_type='contourf', label="Solids (g m$^{-2}$)", figsize=(10,10), show_mesh=False, levels=[0,1,250,2000])
ax.set_xlim([490000,520000])
ax.set_ylim([6520000,6580000])
ax.set_xlabel("Easting (m)")
ax.set_ylabel("Northing (m)")
plt.savefig('meanDep.png',dpi=300,bbox_inches="tight")

[tomadams-ssf]

Have just installed and tested latest updates (version 1.0b0) and confirm this has now been resolved - thanks!