From 77a8027cbc37fd25d12678f2f229fba1fc14c043 Mon Sep 17 00:00:00 2001 From: SorooshMani-NOAA Date: Thu, 1 Aug 2024 14:20:03 -0400 Subject: [PATCH] Add more interp test for mac os issue with tide amp and phase --- tests/test_tidal_run.py | 19 +++++++++++++++++++ 1 file changed, 19 insertions(+) diff --git a/tests/test_tidal_run.py b/tests/test_tidal_run.py index 0804426..fdc9e79 100644 --- 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')