Split up 2d and 3d flexible geometry

[ecomodeller]

No description provided.

[ecomodeller]

A 3d geometry is still fairly complex,, but that is expected, but a 2d geometry is now easier to understand.

[ecomodeller]

Current state of FM geometries

classDiagram

  class _Geometry{
+projection
}

class _GeometryFM{
+node_coordinates
+element_table
}

class GeometryFM2D{
+interp2d()
+get_element_area()
}

class _GeometryFMLayered{
- _n_layers
- _n_sigma
+ to_2d_geometry()
}

class GeometryFM3D{
+ plot()
+ boundary_polylines()
+ find_index()
}

class GeometryFMVerticalProfile{
+plot()
+ find_index()
}
  _Geometry <|-- GeometryFMPointSpectrum
  _Geometry <|-- _GeometryFM
  _GeometryFM <|-- GeometryFM2D
  GeometryFM2D   <|--  _GeometryFMSpectrum
  _GeometryFMSpectrum <|-- GeometryFMAreaSpectrum
  _GeometryFMSpectrum <|-- GeometryFMLineSpectrum
  _GeometryFM <|-- _GeometryFMLayered
  _GeometryFMLayered <|-- GeometryFM3D
  _GeometryFMLayered <|-- GeometryFMVerticalProfile
  GeometryFM3D <|-- GeometryFMVerticalColumn

Loading

Not entirely sure what the GeometryFMVerticalColumn does 🤔

[ecomodeller]

It is possible to subset with a polygon, but the definition if an element is part of the area or not is not yet possible.

Currently an element is selected if the element center is inside the polygon. I would like to be able to select elements which are either entirely within the polygon or has any part that intersects with the polygon.

dss = ds.sel(area=polygon) # no option on how