support for multiple topologies

[veenstrajelmer]

Todo:

• merge_partitions fails for datasets with multiple topologies xugrid#134
• merge_partitions fails for grid with long_culverts xugrid#86
• update minimal xugrid version to 0.9.0 in pyproject.toml (includes above fixes)
• prevent calls to uds.grid in dfmt #363
• 1D2D grid generation overview #782
• the contacts are not read as topology by xugrid

Usecase is GTSM + rivers:

Additional 1D issue (network is 2D saved as 1D), but add_network_cellinfo was added:

file_nc = r'p:\1230882-emodnet_hrsm\GTSMv5.0\grid\gtsm_200s_5000m_eu5000m_ca5000m_v4_net.nc'
uds = xu.open_dataset(file_nc)
# uds_sel = uds.ugrid.to_nonperiodic(xmax=180) # TODO: raises "AttributeError: 'Ugrid1d' object has no attribute 'to_nonperiodic'"
uds_sel = uds.ugrid.sel(x=slice(-75, 2.5), y=slice(10, 60)) #TODO: does not subset

Plot gtsm partitioned mapoutput (requires xugrid 0.9.0):

import matplotlib.pyplot as plt
plt.close("all")
import dfm_tools as dfmt

file_nc = r'p:\1230882-emodnet_hrsm\GTSMv5.0\SO_NHrivGTSM\computations\BD014_fix_mapformat4\output\gtsm_model_0*_map.nc'
uds = dfmt.open_partitioned_dataset(file_nc, remove_edges=False, remove_ghost=False)

fig,ax = plt.subplots()
uds.grids[0].plot(ax=ax)
uds.grids[1].plot(ax=ax, edgecolor='crimson')

# re-adjust axis limits
x_min = min([x.bounds[0] for x in uds.grids])
y_min = min([x.bounds[1] for x in uds.grids])
x_max = max([x.bounds[2] for x in uds.grids])
y_max = max([x.bounds[3] for x in uds.grids])
ax.set_ylim(y_min, y_max)
ax.set_xlim(x_min, x_max)

# plot improvements
dfmt.plot_coastlines(zorder=2, min_area=50)
dfmt.plot_borders()