From 64e09be8ba47c38db8dc4c34700ae0f16e2b4de1 Mon Sep 17 00:00:00 2001
From: Lionel GUEZ <guez@lmd.ens.fr>
Date: Tue, 21 Mar 2023 17:01:39 +0100
Subject: [PATCH] Bug fix: plot geodesics
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit
Otherwise we can have a line across the globe between two points with
longitudes close to -180° and 180°.
---
Trajectories/plot_traj.py | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/Trajectories/plot_traj.py b/Trajectories/plot_traj.py
index 0df6f206..e0dc3c13 100755
--- a/Trajectories/plot_traj.py
+++ b/Trajectories/plot_traj.py
@@ -61,7 +61,7 @@ if __name__ == "__main__":
e_overestim = int(f.read())
SHPC = util_eddies.SHPC_class(args.SHPC, def_orient = args.orientation)
- src_crs = ccrs.PlateCarree()
+ src_crs = ccrs.Geodetic()
projection = ccrs.PlateCarree(central_longitude = 110)
fig, ax = plt.subplots(subplot_kw = {"projection": projection})
random.seed(0)
--
GitLab