HandyPythonScripts.md

@@ -1959,6 +1959,14 @@ def linregress_3D(y_array):
 
 
 def mmm_skirving(cube):
+    try:
+        iris.coord_categorisation.add_year(cube, 'time', name='year')
+    except:
+        pass
+    try:
+        iris.coord_categorisation.add_month(cube, 'time', name='month')
+    except:
+        pass
     cube = cube.aggregated_by(['year','month'], iris.analysis.MEAN)
     print('calculating NOAA Skirving MMM for month:')
 #     missing_data_value_greater_than = -32768.0