Weirs: 5000-point limit on polyline features

[MartijnBregman]

Hello everyone!

I'm considering developing a large-scale SFINCS model in which I want to represent elevated features (such as levees and roads) as weirs. Some of the individual polylines already contain several hundred points, making it easy to reach the 5,000-point limit (as indicated in the link below). Even after simplifying the polylines, staying below that limit would still be challenging. Is the 5,000-point limit still applicable in the latest model version? If so, is the limit intended to simplify model logistics, or is it necessary to prevent adding too much computational burden from the additional weir calculations?

Looking forward to hearing more about this!
Thanks, Martijn

https://github.com/Deltares/SFINCS/blob/4590505d43c146fd1e95c0deb61217ec8db14d22/docs/input_structures.rst#L87

[Leynse]

Hi Martijn, thanks for asking. We don't limit to 5000 points anymore. The model checks how many points there are, and allocates that amount. While running SFINCS the number of structure uv points found is displayed, e.g.:
Info : reading weir file
Info : 7932 structure u/v points found

So you can check whether this number matches (roughly) what you expect.
Note that SFINCS displays the points found after snapping to the grid (max 2 per grid cell), not how many were specified in the input.
I'm changing it in the docs now, thanks!

[Leynse]

To add - the snapped coordinates are available in sfincs_his.nc as structure_x, structure_y & structure_height from SFINCS v2.0.2 onwards.
So you can use that as well after running SFINCS to check the found/snapped weirs

[Leynse]

Update in docs: #161