Add velocity() function

AlexeyPechnikov · Apr 26, 2024 · 2aef4e6 · 2aef4e6
1 parent 3d93878
commit 2aef4e6
Showing 1 changed file with 6 additions and 0 deletions.
diff --git a/pygmtsar/pygmtsar/Stack_stl.py b/pygmtsar/pygmtsar/Stack_stl.py
@@ -12,6 +12,12 @@
 
 class Stack_stl(Stack_tidal):
 
+    def velocity(self, data):
+        years = ((data.date.max() - data.date.min()).dt.days/365.25).item()
+        #print ('years', np.round(years, 3))
+        velocity = data.mean('date')/years
+        return velocity
+
     def trend(self, data, deg=1):
         import xarray as xr