Add total travel time in trace meta data

README.md

@@ -54,5 +54,5 @@ Some of the papers that you may want to cite when using `timflow` are:
   * M. Bakker. 2013. Semi-analytic modeling of transient multi-layer flow with TTim.
   Hydrogeology Journal, 21: 935-943. [https://doi.org/10.1007/s10040-013-0975-2](https://doi.org/10.1007/s10040-013-0975-2)
   * M .Bakker. 2013. Analytic modeling of transient multi-layer flow. In: Advances in
-  Hydrogeology, edited by P Mishra and K Kuhlman, Springer, Heidelberg, 95-114. Available [here](https://github.com/mbakker7/timflow/blob/master/papers/2013_bakker_ttim_theory.pdf).
+  Hydrogeology, edited by P Mishra and K Kuhlman, Springer, Heidelberg, 95-114. Available [here](https://github.com/timflow-org/timflow/blob/main/papers/2013_bakker_ttim_theory.pdf).
 

timflow/steady/trace.py

@@ -57,10 +57,11 @@ def timtraceline(
         if True, return layers numbers
     metadata: boolean
         if False, return xyzt array or xyzt array plus layer array
-        if True, return list of result dictionaries with three entries:
+        if True, return list of result dictionaries with four entries:
         - "trace": np.array(xyzt)
         - "message": termination message
         - "complete": True if terminated correctly
+        - "total_travel_time": final time value in trace
     """
     verbose = False  # used for debugging
     if win is None:
@@ -272,7 +273,12 @@ def timtraceline(
     if not silent:
         print(message)
     if metadata:
-        result = {"trace": np.array(xyzt), "message": message, "complete": terminate}
+        result = {
+            "trace": np.array(xyzt),
+            "message": message,
+            "complete": terminate,
+            "total_travel_time": xyzt[-1][-1],
+        }
         if returnlayers:
             result["layers"] = layerlist
     else: