diff --git a/Trajectories/plot_traj.py b/Trajectories/plot_traj.py index f9e93f4a779c1bea32625606f3ccea922afcb1b8..0df6f2065a90f2730fbccdd0281974e4d0de194b 100755 --- a/Trajectories/plot_traj.py +++ b/Trajectories/plot_traj.py @@ -20,8 +20,9 @@ def get_extr_coord(traj, e_overestim, SHPC): def plot_single_traj(traj, e_overestim, SHPC, ax, src_crs, annotate_flag): x, y = get_extr_coord(traj, e_overestim, SHPC) - ax.plot(x, y, color = "red", transform = src_crs) - ax.plot(x[0], y[0], marker = "s", color = "black", transform = src_crs) + ax.plot(x, y, color = "red", linewidth = 0.5, transform = src_crs) + ax.plot(x[0], y[0], marker = "s", markersize = 2, color = "black", + transform = src_crs) if annotate_flag: ax.annotate(str(traj[0]), @@ -35,6 +36,7 @@ if __name__ == "__main__": import json import cartopy.crs as ccrs import argparse + import cartopy.feature as cfeature parser = argparse.ArgumentParser() parser.add_argument("expanded_traj", help = "JSon file") @@ -60,7 +62,7 @@ if __name__ == "__main__": SHPC = util_eddies.SHPC_class(args.SHPC, def_orient = args.orientation) src_crs = ccrs.PlateCarree() - projection = ccrs.PlateCarree(central_longitude = 90) + projection = ccrs.PlateCarree(central_longitude = 110) fig, ax = plt.subplots(subplot_kw = {"projection": projection}) random.seed(0) n_long_traj = 0 @@ -76,7 +78,7 @@ if __name__ == "__main__": plot_single_traj(traj, e_overestim, SHPC, ax, src_crs, args.annotate) - ax.coastlines() + ax.add_feature(cfeature.LAND, edgecolor = "black") ax.gridlines(draw_labels = True) if args.min_duration > 1: