option added to make cross-shore coords of figure start at offsetted …

…value

xbTools/xbeachpost.py

@@ -21,7 +21,7 @@
 import warnings
 import re
 # toolbox specific import
-from .general.geometry import rotate_grid
+from .general.geometry import rotate_grid, path_distance
 
 class XBeachModelAnalysis():
     '''
@@ -60,7 +60,8 @@ def __init__(self, fname, model_path):
                           'u_max', 'v_max', 'ue_max', 've_max', 'Subg_max', 'Svbg_max', 'Susg_max', 'Svsg_max',
                           'u_min', 'v_min', 'ue_min', 've_min', 'Subg_min', 'Svbg_min', 'Susg_min', 'Svsg_min']
 
-
+        self._cross_offset = 0
+
     def get_metadata(self):
         '''
         function to print metadata from XBlog.txt
@@ -186,9 +187,17 @@ def get_params(self):
         if params['npointvar'] > 0:
             i0 = [i for i, var in enumerate(dat) if 'npoints' in var][0]
             points = dat[i0 + 1:i0 + int(params['npoints'] + 1)]
-            x = [float( re.split('\t| ',t)[0] ) for t in points]
-            y = [float(re.split('\t| ',t)[1] ) for t in points]
-            name = [t[2].strip() for t in points]
+
+            # get points 
+            p_spl = [ re.split('\t| ',t) for t in points]
+            # remove empty entries
+            p_spl = [[p for p in p_spli if p.split()] for p_spli in p_spl]
+            # xy
+            x = [float(p[0] ) for p in p_spl]
+            y = [float(p[1]) for p in p_spl]
+            # if names are there list of names else index
+            name = [p[2].strip() if len(p)>2 else str(ip) for ip, p in enumerate(p_spl)]
+
             params['points'] = dict(zip(name, zip(x, y)))
         else:
             params['points'] = {}
@@ -387,31 +396,9 @@ def load_output_coordinates(self):
             self.var['station_x'] = ds.variables['pointx'][:]
             self.var['station_y'] = ds.variables['pointy'][:]
 
-        def path_distance(polx, poly):
-            '''
-            computes the distance along a polyline
-            ----------
-            polx : TYPE array
-                X COORDINATES.
-            poly : TYPE array
-                Y COORDINATES.
-
-            Returns: TYPE array
-                PATHDISTANCE
-            -------
-            python alternative to matlab function.
-
-            '''
-            dx = np.diff(polx)
-            dy = np.diff(poly)
-
-            dr = np.sqrt(dx ** 2 + dy ** 2)
+        cross = path_distance(x[0, :], y[0, :])
+        self.var['cross'] = cross + self._cross_offset
 
-            pathdistance = np.insert(np.cumsum(dr), 0, 0, axis=0)
-
-            return pathdistance
-
-        self.var['cross'] = path_distance(x[0, :], y[0, :])
         self.var['along'] = path_distance(x[:, 0], y[:, 0])
 
         # self.var['localy'], self.var['localx'] = np.meshgrid(self.var['cross'], self.var['along'])
@@ -435,6 +422,13 @@ def path_distance(polx, poly):
 
         return
 
+    def add_cross_offset(self, offset):
+        """
+        adds an offset value to the computed cross-shore local coordinate system (such that it does not start at zero at offshore boundary)
+        """
+        self._cross_offset = offset
+        self.load_output_coordinates()
+
     def load_modeloutput(self, var):
         """_summary_