1D2D grid generation overview

[veenstrajelmer]

Compute contact points between mesh1d and mesh2d:

import matplotlib.pyplot as plt
import xugrid as xu
from meshkernel import MeshKernel, ProjectionType

file_nc = r'p:\1230882-emodnet_hrsm\GTSMv5.0\SO_NHrivGTSM\UGgtsm_snipbd_v2_withriv_cellinfo_net.nc'
uds = xu.open_dataset(file_nc)
fig,ax = plt.subplots()
uds.grids[0].plot(ax=ax, edgecolor='crimson')
uds.grids[1].plot(ax=ax)
ax.set_ylim(None,uds.grids[0].bounds[-1])

mk = MeshKernel(projection=ProjectionType.SPHERICAL)

mk.mesh1d_set(uds.grids[0].meshkernel.mesh1d_get())
mk.mesh2d_set(uds.grids[1].meshkernel.mesh2d_get())

mk.contacts_compute_boundary(node_mask=None, polygons=None, search_radius=0.2)
# TODO: mk.contacts_compute_boundary() requires 'node_mask', 'polygons', and 'search_radius', alternative?

Request to make the arguments optional: Deltares/MeshKernelPy#81

1D2D grid generation:

• 1D: mesh.py#L92
• 1D2D: mesh.py#L503
• connection: mesh.py#L470
• mesh1d_set: meshkernel.py#L919
  
• contacts_compute_single: meshkernel.py#L1006
  
• contacts_compute_boundary: meshkernel.py#L1096 (can connect multiple 2d faces to the closest 1d node)
• maybe easiest way is to do this with hydromt_delft3dfm