From 88e47d88861eba33d33099eb4bea7417db9cee2b Mon Sep 17 00:00:00 2001 From: svangennip Date: Thu, 3 Nov 2022 16:50:26 +0100 Subject: [PATCH 1/2] modify grid.py for longitudinal quadrant --- src/py_eddy_tracker/dataset/grid.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/py_eddy_tracker/dataset/grid.py b/src/py_eddy_tracker/dataset/grid.py index 9345bf45..f40731ca 100644 --- a/src/py_eddy_tracker/dataset/grid.py +++ b/src/py_eddy_tracker/dataset/grid.py @@ -177,7 +177,7 @@ def _fit_circle_path(vertice): return 0, -90, nan, nan centlon, centlat, eddy_radius, err = fit_circle(c_x, c_y) centlon, centlat = local_to_coordinates(centlon, centlat, lon0, lat0) - centlon = (centlon - lon0 + 180) % 360 + lon0 - 180 + #centlon = (centlon - lon0 + 180) % 360 + lon0 - 180 return centlon, centlat, eddy_radius, err From dfbced0bcc47376d7e40f9b1b9f0273b57e0334a Mon Sep 17 00:00:00 2001 From: svangennip Date: Mon, 7 Nov 2022 16:24:49 +0100 Subject: [PATCH 2/2] Change of lat-lon completed --- src/py_eddy_tracker/dataset/grid.py | 10 +++++----- 1 file changed, 5 insertions(+), 5 deletions(-) diff --git a/src/py_eddy_tracker/dataset/grid.py b/src/py_eddy_tracker/dataset/grid.py index f40731ca..3041913a 100644 --- a/src/py_eddy_tracker/dataset/grid.py +++ b/src/py_eddy_tracker/dataset/grid.py @@ -950,11 +950,11 @@ def eddy_identification( eddies.time[:] = (date - datetime(1950, 1, 1)).total_seconds() / 86400.0 # normalization longitude between 0 - 360, because storage have an offset on 180 - eddies.lon_max[:] %= 360 - eddies.lon[:] %= 360 - ref = eddies.lon - 180 - eddies.contour_lon_e[:] = ((eddies.contour_lon_e.T - ref) % 360 + ref).T - eddies.contour_lon_s[:] = ((eddies.contour_lon_s.T - ref) % 360 + ref).T + #eddies.lon_max[:] %= 360 + #eddies.lon[:] %= 360 + #ref = eddies.lon - 180 + #eddies.contour_lon_e[:] = ((eddies.contour_lon_e.T - ref) % 360 + ref).T + #eddies.contour_lon_s[:] = ((eddies.contour_lon_s.T - ref) % 360 + ref).T a_and_c.append(eddies) if in_h_unit is not None: