Calculation of derived coords points and bounds is always lazy.

Author: pp-mo

lib/iris/aux_factory.py

@@ -188,7 +188,7 @@ def _nd_bounds(coord, dims, ndim):
         # Transpose to be consistent with the Cube.
         sorted_pairs = sorted(enumerate(dims), key=lambda pair: pair[1])
         transpose_order = [pair[0] for pair in sorted_pairs] + [len(dims)]
-        bounds = coord.core_bounds()
+        bounds = coord.lazy_bounds()
         if dims:
             bounds = bounds.transpose(transpose_order)
 
@@ -216,7 +216,7 @@ def _nd_points(coord, dims, ndim):
         # Transpose to be consistent with the Cube.
         sorted_pairs = sorted(enumerate(dims), key=lambda pair: pair[1])
         transpose_order = [pair[0] for pair in sorted_pairs]
-        points = coord.core_points()
+        points = coord.lazy_points()
         if dims and transpose_order != list(range(len(dims))):
             points = points.transpose(transpose_order)