From d3919797a21987e0e481c62695634a5d637e0ce3 Mon Sep 17 00:00:00 2001
From: Paul RASCLE <paul.rascle@openfields.fr>
Date: Mon, 31 Aug 2020 11:21:27 +0200
Subject: [PATCH] fix fitShapesPointsToMeshEdges with double points, no split

---
 src/HYDROTools/shapesGroups.py | 2 ++
 1 file changed, 2 insertions(+)

diff --git a/src/HYDROTools/shapesGroups.py b/src/HYDROTools/shapesGroups.py
index 8bcdd4d0..f4673d80 100644
--- a/src/HYDROTools/shapesGroups.py
+++ b/src/HYDROTools/shapesGroups.py
@@ -340,6 +340,8 @@ def fitShapePointsToMesh(freeBorderShapefile, shapefileToAdjust, outputDirectory
                 chaincoords.append(sfta.points[i])
         for i in range(i1):
             chaincoords.append(sfta.points[i])
+        if discountLastSftaPoint:
+            chaincoords.append(fbs.points[ifb1]) # close shape when first point if superposed with last
         w.line([chaincoords])
         w.record(chainName)
 
-- 
2.39.2