Merge pull request NCAR#174 from NCAR/mercator_fix

michaelavs · web-flow · commit 74dcbb2a3b91 · 2022-05-03T10:22:03.000-06:00
Mercator fix
diff --git a/src/wrf/projection.py b/src/wrf/projection.py
@@ -18,6 +18,7 @@
     from Ngl import Resources
 
 
+
 if cartopy_enabled():
     class MercatorWithLatTS(crs.Mercator):
         """A :class:`cartopy.crs.Mercator` subclass that adds support for
@@ -60,6 +61,10 @@ def __init__(self, central_longitude=0.0,
                             ("units", "m")]
             super(crs.Mercator, self).__init__(proj4_params, globe=globe)
 
+            # Need to have x/y limits defined for the initial hash which
+            # gets used within transform_points for caching
+            self._x_limits = self._y_limits = None
+
             # Calculate limits.
             limits = self.transform_points(
                 crs.Geodetic(),
@@ -80,7 +85,8 @@ def __init__(self, central_longitude=0.0,
             self._x_limits = tuple(xlimits)
             self._y_limits =  tuple(limits[..., 1])
 
-            self._threshold = np.diff(self.x_limits)[0] / 720
+            self._threshold = min(np.diff(self.x_limits)[0] / 720,
+                                  np.diff(self.y_limits)[0] / 360)
 
 
 def _ismissing(val, islat=True):
