bugfix: missing computation of origin uncertainty

floenavi/__init__.py

@@ -876,7 +876,7 @@ class FloeNaviGridSolver(object):
 
         # NOTE: Currently, this is a simple point-wise mean of the parameters
         # TODO: Add quality check (e.g. the min-max distance)
-        parameter_names = ["origin_longitude", "origin_latitude", "axes_heading"]
+        parameter_names = ["origin_longitude", "origin_latitude", "axes_heading", "origin_uncertainty_m"]
         parameter_units = dict(origin_longitude="deg", origin_latitude="deg", origin_uncertainty_m="m",
                                axes_heading="deg")
         n_records = self.solutions[0].n_records
@@ -908,6 +908,10 @@ class FloeNaviGridSolver(object):
             if parameter_name == "axes_heading":
                 self.heading_uncertainty_deg = np.nanstd(var_array, axis=1)
 
+            # Compute origin uncertainty (point-wise mean of different solutions)
+            if parameter_name == "origin_uncertainty_m":
+                self.origin_uncertainty_m = np.nanmean(var_array, axis=1)
+
     @staticmethod
     def basestations_grid_distance(bs1, bs2):
         return np.sqrt((bs2.xc_grid-bs1.xc_grid)**2.0 + (bs2.yc_grid-bs1.yc_grid)**2.0)