Design simple FM/FM nesting workflow

[veenstrajelmer]

• Nesthd1:
  • input: rooster grof (evt ook epsg code, of uit rooster lezen)
  • input: pli-line met bnd langs rooster-fijn (incl epsg code, default is WGS84 maar check dan: coords.min()>=-180 en coords.max()<=180)
  • action: kdtree van cellcenters grove rooster
  • action: find 3 nearest neighbors (cartesian) or find cell numbers and connected cells via meshgrid oid
  • action: generate list of cellscenters and connected cellcenters
  • action: drop duplicate coordinates, drop coordinates with distance larger than np.sqrt(cellarea.mean() of nestcells)
  • returns: obsfile met coordinaten van deze cellcenter values
• Run grof model met obsfile
• Nesthd2:
  • input: pli-line
  • input: hisfile grove model
  • (input: obsfile, to check if all obspoints are present in ncfile (geen adminfile nodig))
  • action: make kdtree van obspoints, query per pli punt (incl distance, cartesian), invdist avg of neighbors, results in timeseries
    • alternative: hisfile.interp() met alle pli punten (zoals interpolate_nc_to_bc() function) (lijkt niet zomaar mogelijk, maar kan ongetwijfeld)
  • action: put timeseries in forcingmodel and write to bc file
  • eerste opzet in "workinprogress_interpolate_his.py" example script
• er is een idee om interpolatie in te bouwen in de FM kernel, dan zou het nog veel makkelijker worden (plifile omzetten naar xyn, grof model draaien met xyn, history inlezen en omzetten naar nc/bc

First setup: https://github.com/Deltares/dfm_tools/blob/main/tests/examples_workinprogress/workinprogress_nestingFMtoFM.py

Also request from RS/PK to add Neuman boundary, where inflowing direction of grid should be known (possible?). Or alternatively, combination of waterlevel/ux/uy boundary, which is simpler to implement.

Related, nesting SFINCS in FM is done in GC with:

• Generated FM obsfile with: Add function to modelbuilder to generate observations along coastline #578
• Nest SFINCS from FM hisfile: https://github.com/Deltares/gc-japan-cyclone/blob/issue_22/src/sfincs/sfincs_update_forcing.py
• Uses functions from: https://github.com/Deltares/gc-japan-cyclone/blob/issue_22/src/sfincs/sfincs_utils.py

Also request from WO to add Riemann boundary, so combination of waterlevel and discharge (latter from cross-sections). For morphology we also require nesting bed level changes (form obsfiles) and sediment transport (form cross-sections).

Also align with nesting in ocean models like CMEMS, maybe also consider nesting in other existing models (shyfem, schism, delft3d4, etc).