Fix for projections using a round earth (#336)

mpvginde · web-flow · commit 8dca14afbf74 · 2023-08-16T20:14:55.000+02:00
diff --git a/pysteps/visualization/utils.py b/pysteps/visualization/utils.py
@@ -170,6 +170,9 @@ def proj4_to_cartopy(proj4str):
             proj_dict["lat_1"],
             proj_dict["lat_2"],
         )
+    if "R" in proj_dict.keys():
+        globe_kwargs["semimajor_axis"] = proj_dict["R"]
+        globe_kwargs["semiminor_axis"] = proj_dict["R"]
 
     if globe_kwargs:
         globe = ccrs.Globe(**globe_kwargs)

Original file line number	Diff line number	Diff line change
`@@ -170,6 +170,9 @@ def proj4_to_cartopy(proj4str):`
`170`	`170`	`proj_dict["lat_1"],`
`171`	`171`	`proj_dict["lat_2"],`
`172`	`172`	`)`
	`173`	`+ if "R" in proj_dict.keys():`
	`174`	`+ globe_kwargs["semimajor_axis"] = proj_dict["R"]`
	`175`	`+ globe_kwargs["semiminor_axis"] = proj_dict["R"]`
`173`	`176`
`174`	`177`	`if globe_kwargs:`
`175`	`178`	`globe = ccrs.Globe(**globe_kwargs)`