Deal with schism bug in station output

Manually set the time index of staout files
Fix #174

pyposeidon/schism.py

@@ -1903,14 +1903,22 @@ def get_station_sim_data(self, **kwargs):
 
         skiprows = kwargs.get("station_skiprows", 0)
 
+        try:
+            dt = self.parameters["dt"]
+        except:
+            dt = self.params["core"]["dt"]
+
         dfs = []
         for idx in vals.index:
             obs = np.loadtxt(sfiles[idx], skiprows=skiprows)
             df = pd.DataFrame(obs)
             df = df.set_index(0)
             df.index.name = "time"
             df.columns.name = vals.loc[idx, "variable"]
-            df.index = pd.to_datetime(dstamp) + pd.to_timedelta(df.index, unit="S")
+            # deal with schism bug
+            ns = np.arange(1, df.shape[0] + 1)
+            df.index = pd.to_datetime(dstamp) + pd.to_timedelta(ns * dt, unit="S")
+            #            df.index = pd.to_datetime(dstamp) + pd.to_timedelta(df.index, unit="S")
             pindex = pd.MultiIndex.from_product([df.T.columns, df.T.index])
 
             r = pd.DataFrame(df.values.flatten(), index=pindex, columns=[vals.loc[idx, "variable"]])