Add more interp test for mac os issue with tide amp and phase

oceanmodeling · Aug 1, 2024 · 77a8027 · 77a8027
1 parent 81141d6
commit 77a8027
Showing 1 changed file with 19 additions and 0 deletions.
diff --git a/tests/test_tidal_run.py b/tests/test_tidal_run.py
@@ -70,6 +70,25 @@ def test_tidal_run(shinnecock_mesh_directory):
             df_adjy = pd.DataFrame(yi[:, 0])
             df_adjy.to_csv(output_directory / f'debug_adjy_{constituent}_{index}.csv')
             ###############
+            dataset = ds._get_dataset('elevation', constituent)
+            zi = dataset['PHAS'][yidx, xidx]
+            mask = ~zi.mask
+            if mask.size == 1:
+                mask = np.array(zi * mask, dtype=bool)
+            xi2 = xi[mask].flatten()
+            yi2 = yi[mask].flatten()
+            zi = zi[mask].flatten()
+            pd.DataFrame({'xi': xi2, 'yi': yi2, 'zi': zi}).to_csv(output_directory / f'debug_adjxyz_phase_{constituent}_{index}.csv')
+
+            zi = dataset['AMPL'][yidx, xidx]
+            mask = ~zi.mask
+            if mask.size == 1:
+                mask = np.array(zi * mask, dtype=bool)
+            xi2 = xi[mask].flatten()
+            yi2 = yi[mask].flatten()
+            zi = zi[mask].flatten()
+            pd.DataFrame({'xi': xi2, 'yi': yi2, 'zi': zi}).to_csv(output_directory / f'debug_adjxyz_ampl_{constituent}_{index}.csv')
+            ###############
             amp, phase = mesh.forcings.tides.tidal_dataset(constituent, vertices)
             df_tide = pd.DataFrame({'amp': amp, 'phase': phase})
             df_tide.to_csv(output_directory / f'debug_tide_{constituent}_{index}.csv')