SciTools
diff --git a/‎docs/src/userguide/cube_maths.rst‎
Lines changed: 2 additions & 2 deletions b/‎docs/src/userguide/cube_maths.rst‎
Lines changed: 2 additions & 2 deletions
diff --git a/‎docs/src/whatsnew/3.0.1.rst‎
Lines changed: 1 addition & 2 deletions b/‎docs/src/whatsnew/3.0.1.rst‎
Lines changed: 1 addition & 2 deletions
diff --git a/‎docs/src/whatsnew/3.0.2.rst‎
Lines changed: 571 additions & 0 deletions b/‎docs/src/whatsnew/3.0.2.rst‎
Lines changed: 571 additions & 0 deletions
diff --git a/‎docs/src/whatsnew/3.0.3.rst‎
Lines changed: 586 additions & 0 deletions b/‎docs/src/whatsnew/3.0.3.rst‎
Lines changed: 586 additions & 0 deletions
diff --git a/‎docs/src/whatsnew/3.0.4.rst‎
Lines changed: 620 additions & 0 deletions b/‎docs/src/whatsnew/3.0.4.rst‎
Lines changed: 620 additions & 0 deletions
diff --git a/‎docs/src/whatsnew/index.rst‎
Lines changed: 3 additions & 0 deletions b/‎docs/src/whatsnew/index.rst‎
Lines changed: 3 additions & 0 deletions
diff --git a/‎lib/iris/_lazy_data.py‎
Lines changed: 3 additions & 1 deletion b/‎lib/iris/_lazy_data.py‎
Lines changed: 3 additions & 1 deletion
diff --git a/‎lib/iris/coords.py‎
Lines changed: 14 additions & 23 deletions b/‎lib/iris/coords.py‎
Lines changed: 14 additions & 23 deletions
diff --git a/‎lib/iris/cube.py‎
Lines changed: 38 additions & 24 deletions b/‎lib/iris/cube.py‎
Lines changed: 38 additions & 24 deletions
diff --git a/‎lib/iris/fileformats/_nc_load_rules/actions.py‎
Lines changed: 17 additions & 3 deletions b/‎lib/iris/fileformats/_nc_load_rules/actions.py‎
Lines changed: 17 additions & 3 deletions
@@ -157,7 +157,7 @@ Let's extend this example slightly, by taking a slice from the middle
     air_temperature / (K)               (longitude: 49; time: 240)
 
 Compared to our original time-series, the *air_temp_T_slice* cube has one
-less dimension *and* it's shape if different. However, this doesn't prevent
+less dimension *and* its shape is different. However, this doesn't prevent
 us from performing cube arithmetic with it, thanks to the extended cube
 broadcasting behaviour::
 
@@ -237,7 +237,7 @@ by a cube with unit ``'1'`` will preserve units, so the cube ``temperature``
 will be given the same units as are in ``pot_temperature``. It should be
 noted that some combinations of units, particularly those involving power
 operations, will not result in a valid unit and will cause the calculation
-to fail. For example, a cube ``a`` had units ``'m'`` then ``a ** 0.5``
+to fail. For example, if a cube ``a`` had units ``'m'`` then ``a ** 0.5``
 would result in an error since the square root of a meter has no meaningful
 unit (if ``a`` had units ``'m2'`` then ``a ** 0.5`` would result in a cube
 with units ``'m'``).
 
@@ -24,12 +24,11 @@ This document explains the changes made to Iris for this release
       factory, and the associated derived coordinate will be missing. (:pull:`3965`)
 
 
-.. dropdown:: :opticon:`report` Release Highlights
+.. dropdown:: :opticon:`report` v3.0.0 Release Highlights
    :container: + shadow
    :title: text-primary text-center font-weight-bold
    :body: bg-light
    :animate: fade-in
-   :open:
 
    The highlights for this major release of Iris include:
 
 
@@ -11,6 +11,9 @@ Iris versions.
    :maxdepth: 1
 
    latest.rst
+   3.0.4.rst
+   3.0.3.rst
+   3.0.2.rst
    3.0.1.rst
    3.0.rst
    2.4.rst
 
@@ -191,7 +191,9 @@ def as_lazy_data(data, chunks=None, asarray=False):
     if isinstance(data, ma.core.MaskedConstant):
         data = ma.masked_array(data.data, mask=data.mask)
     if not is_lazy_data(data):
-        data = da.from_array(data, chunks=chunks, asarray=asarray)
+        data = da.from_array(
+            data, chunks=chunks, asarray=asarray, meta=np.ndarray
+        )
     return data
 
 
 
@@ -403,28 +403,15 @@ def __binary_operator__(self, other, mode_constant):
         # Note: this method includes bounds handling code, but it only runs
         # within Coord type instances, as only these allow bounds to be set.
 
-        if isinstance(other, _DimensionalMetadata) or not isinstance(
-            other, (int, float, np.number)
-        ):
-
-            def typename(obj):
-                if isinstance(obj, Coord):
-                    result = "Coord"
-                else:
-                    # We don't really expect this, but do something anyway.
-                    result = self.__class__.__name__
-                return result
-
-            emsg = "{selftype} {operator} {othertype}".format(
-                selftype=typename(self),
-                operator=self._MODE_SYMBOL[mode_constant],
-                othertype=typename(other),
+        if isinstance(other, _DimensionalMetadata):
+            emsg = (
+                f"{self.__class__.__name__} "
+                f"{self._MODE_SYMBOL[mode_constant]} "
+                f"{other.__class__.__name__}"
             )
             raise iris.exceptions.NotYetImplementedError(emsg)
 
-        else:
-            # 'Other' is an array type : adjust points, and bounds if any.
-            result = NotImplemented
+        if isinstance(other, (int, float, np.number)):
 
             def op(values):
                 if mode_constant == self._MODE_ADD:
@@ -441,8 +428,14 @@ def op(values):
 
             new_values = op(self._values_dm.core_data())
             result = self.copy(new_values)
+
             if self.has_bounds():
                 result.bounds = op(self._bounds_dm.core_data())
+        else:
+            # must return NotImplemented to ensure invocation of any
+            # associated reflected operator on the "other" operand
+            # see https://docs.python.org/3/reference/datamodel.html#emulating-numeric-types
+            result = NotImplemented
 
         return result
 
@@ -461,8 +454,7 @@ def __div__(self, other):
     def __truediv__(self, other):
         return self.__binary_operator__(other, self._MODE_DIV)
 
-    def __radd__(self, other):
-        return self + other
+    __radd__ = __add__
 
     def __rsub__(self, other):
         return (-self) + other
@@ -473,8 +465,7 @@ def __rdiv__(self, other):
     def __rtruediv__(self, other):
         return self.__binary_operator__(other, self._MODE_RDIV)
 
-    def __rmul__(self, other):
-        return self * other
+    __rmul__ = __mul__
 
     def __neg__(self):
         values = -self._core_values()
 
@@ -3193,24 +3193,26 @@ def _intersect_modulus(
             # and call the new bounds = the new points + the difference.
             pre_wrap_delta = np.diff(coord.bounds[inside_indices])
             post_wrap_delta = np.diff(bounds[inside_indices])
-            close_enough = np.allclose(pre_wrap_delta, post_wrap_delta)
-            if not close_enough:
-                split_cell_indices, _ = np.where(
-                    pre_wrap_delta != post_wrap_delta
-                )
-
-                # Recalculate the extended minimum.
+            split_cell_indices, _ = np.where(
+                ~np.isclose(pre_wrap_delta, post_wrap_delta)
+            )
+            if split_cell_indices.size:
                 indices = inside_indices[split_cell_indices]
                 cells = bounds[indices]
-                cells_delta = np.diff(coord.bounds[indices])
-
-                # Watch out for ascending/descending bounds
-                if cells_delta[0, 0] > 0:
-                    cells[:, 0] = cells[:, 1] - cells_delta[:, 0]
-                    minimum = np.min(cells[:, 0])
-                else:
-                    cells[:, 1] = cells[:, 0] + cells_delta[:, 0]
-                    minimum = np.min(cells[:, 1])
+                if maximum % modulus not in cells:
+                    # Recalculate the extended minimum only if the output bounds
+                    # do not span the requested (minimum, maximum) range.  If
+                    # they do span that range, this adjustment would give unexpected
+                    # results (see #3391).
+                    cells_delta = np.diff(coord.bounds[indices])
+
+                    # Watch out for ascending/descending bounds.
+                    if cells_delta[0, 0] > 0:
+                        cells[:, 0] = cells[:, 1] - cells_delta[:, 0]
+                        minimum = np.min(cells[:, 0])
+                    else:
+                        cells[:, 1] = cells[:, 0] + cells_delta[:, 0]
+                        minimum = np.min(cells[:, 1])
 
             points = wrap_lons(coord.points, minimum, modulus)
 
@@ -3764,37 +3766,49 @@ def __ne__(self, other):
     def __hash__(self):
         return hash(id(self))
 
-    def __add__(self, other):
-        return iris.analysis.maths.add(self, other)
+    __add__ = iris.analysis.maths.add
 
     def __iadd__(self, other):
         return iris.analysis.maths.add(self, other, in_place=True)
 
     __radd__ = __add__
 
-    def __sub__(self, other):
-        return iris.analysis.maths.subtract(self, other)
+    __sub__ = iris.analysis.maths.subtract
 
     def __isub__(self, other):
         return iris.analysis.maths.subtract(self, other, in_place=True)
 
+    def __rsub__(self, other):
+        return (-self) + other
+
     __mul__ = iris.analysis.maths.multiply
-    __rmul__ = iris.analysis.maths.multiply
 
     def __imul__(self, other):
         return iris.analysis.maths.multiply(self, other, in_place=True)
 
+    __rmul__ = __mul__
+
     __div__ = iris.analysis.maths.divide
 
     def __idiv__(self, other):
         return iris.analysis.maths.divide(self, other, in_place=True)
 
-    __truediv__ = iris.analysis.maths.divide
+    def __rdiv__(self, other):
+        data = 1 / self.core_data()
+        reciprocal = self.copy(data=data)
+        return iris.analysis.maths.multiply(reciprocal, other)
 
-    def __itruediv__(self, other):
-        return iris.analysis.maths.divide(self, other, in_place=True)
+    __truediv__ = __div__
+
+    __itruediv__ = __idiv__
+
+    __rtruediv__ = __rdiv__
 
     __pow__ = iris.analysis.maths.exponentiate
+
+    def __neg__(self):
+        return self.copy(data=-self.core_data())
+
     # END OPERATOR OVERLOADS
 
     def collapsed(self, coords, aggregator, **kwargs):
 
@@ -43,11 +43,15 @@
 from functools import wraps
 import warnings
 
+from iris.config import get_logger
 import iris.fileformats.cf
 import iris.fileformats.pp as pp
 
 from . import helpers as hh
 
+# Configure the logger.
+logger = get_logger(__name__, fmt="[%(funcName)s]")
+
 
 def _default_rulenamesfunc(func_name):
     # A simple default function to deduce the rules-name from an action-name.
@@ -385,13 +389,23 @@ def action_ukmo_stash(engine):
     attr_name = "ukmo__um_stash_source"
     attr_value = getattr(var, attr_name, None)
     if attr_value is None:
-        attr_altname = "um_stash_source"  # legacy form
-        attr_value = getattr(var, attr_altname, None)
+        attr_name = "um_stash_source"  # legacy form
+        attr_value = getattr(var, attr_name, None)
     if attr_value is None:
         rule_name += "(NOT-TRIGGERED)"
     else:
         # No helper routine : just do it
-        engine.cube.attributes["STASH"] = pp.STASH.from_msi(attr_value)
+        try:
+            stash_code = pp.STASH.from_msi(attr_value)
+        except (TypeError, ValueError):
+            engine.cube.attributes[attr_name] = attr_value
+            msg = (
+                "Unable to set attribute STASH as not a valid MSI "
+                f'string "mXXsXXiXXX", got "{attr_value}"'
+            )
+            logger.debug(msg)
+        else:
+            engine.cube.attributes["STASH"] = stash_code
 
     return rule_name