Adds diagnostics to identify properties of the MAUL

docs/source/reference/operators.rst

@@ -122,6 +122,15 @@ CSET.operators.mesoscale
 .. automodule:: CSET.operators.mesoscale
    :members:
 
+Precipitation Operators
+~~~~~~~~~~~~~~~~~~~~~~~
+
+CSET.operators.precipitation
+----------------------------
+
+.. automodule:: CSET.operators.precipitation
+   :members:
+
 Pressure Operators
 ~~~~~~~~~~~~~~~~~~
 

src/CSET/cset_workflow/meta/diagnostics/rose-meta.conf

@@ -1807,7 +1807,7 @@ title= Presence of moist absolutely unstable layers (MAULs) in a column
 description=PROCESS-BASED DIAGNOSTIC.
             Determines if a moist absolutely unstable layer (MAUL) is present
             in a grid column regardless of depth.
-            Requires "air_temperature", "air_pressure", and "specfic_humidity" on model levels.
+            Requires "air_temperature", "air_pressure", and "specific_humidity" on model levels.
 help=This diagnostic identifies only whether a MAUL is present within the column.
      It does not provide any information about the properties of the MAUL. A MAUL
      is defined as an unstable layer (equivalent potential temperature gradient
@@ -1817,6 +1817,55 @@ type=python_boolean
 compulsory=true
 sort-key=xp-maulp
 
+[template variables=MAUL_NUMBER]
+ns = Diagnostics/Derived/XPPN
+title= Number of moist absolutely unstable layers (MAULs) in a column
+description=PROCESS-BASED DIAGNOSTIC.
+            Determines the number of moist absolutely unstable layers (MAUL) present
+            in a grid column regardless of depth.
+            Requires "air_temperature", "air_pressure", and "specific_humidity" on model levels.
+help=This diagnostic identifies the numbers of MAUL that are present within the column.
+     It does not provide any information about the properties of the MAUL. A MAUL
+     is defined as an unstable layer (equivalent potential temperature gradient
+     with height less than zero) that has a relatively humidity above '90%' following
+     Takemi and Unuma (2020).
+
+type=python_boolean
+compulsory=true
+sort-key=xp-maulp1
+
+[template variables=MAUL_DEPTH]
+ns = Diagnostics/Derived/XPPN
+title= Depth of moist absolutely unstable layer (MAUL) in a column
+description=PROCESS-BASED DIAGNOSTIC.
+            Determines the depth of the deepest moist absolutely unstable layer (MAUL)
+            in a grid column.
+            Requires "air_temperature", "air_pressure", and "specific_humidity" on model levels.
+help=This diagnostic identifies only the deepest MAUL within the column. A MAUL
+     is defined as an unstable layer (equivalent potential temperature gradient
+     with height less than zero) that has a relatively humidity above '90%' following
+     Takemi and Unuma (2020).
+
+type=python_boolean
+compulsory=true
+sort-key=xp-maulp2
+
+[template variables=MAUL_BASE]
+ns = Diagnostics/Derived/XPPN
+title= Base (AGL) of moist absolutely unstable layer (MAUL) in a column
+description=PROCESS-BASED DIAGNOSTIC.
+            Determines the base height (AGL) of the deepest moist absolutely unstable layer (MAUL)
+            in a grid column.
+            Requires "air_temperature", "air_pressure", and "specific_humidity" on model levels.
+help=This diagnostic identifies the base of the deepest MAUL within the column. A MAUL
+     is defined as an unstable layer (equivalent potential temperature gradient
+     with height less than zero) that has a relatively humidity above '90%' following
+     Takemi and Unuma (2020).
+
+type=python_boolean
+compulsory=true
+sort-key=xp-maulp3
+
 ###################################
 # Ensembles
 [Diagnostics/Ensembles]

src/CSET/cset_workflow/rose-suite.conf.example

@@ -74,6 +74,9 @@ LIGHT_RAIN_PRESENCE_SPATIAL_DIFFERENCE=False
 LIGHT_RAIN_PRESENCE_SPATIAL_PLOT=False
 LOGLEVEL="DEBUG"
 !!LONGITUDE_POINT=0
+MAUL_BASE=False
+MAUL_DEPTH=False
+MAUL_NUMBER=False
 MAUL_PRESENCE=False
 MEAN_STRUCTURAL_SIMILARITY_MLEVEL=False
 MEAN_STRUCTURAL_SIMILARITY_PLEVEL_FIELD=False

src/CSET/loaders/spatial_field.py

@@ -572,6 +572,54 @@ def load(conf: Config):
                 aggregation=False,
             )
 
+    # Number of Moist Absolutely Unstable Layer in a Column
+    if conf.MAUL_NUMBER:
+        for model in models:
+            yield RawRecipe(
+                recipe="number_of_moist_absolutely_unstable_layers_present_in_a_spatial_column.yaml",
+                variables={
+                    "MODEL_NAME": model["name"],
+                    "SUBAREA_TYPE": conf.SUBAREA_TYPE if conf.SELECT_SUBAREA else None,
+                    "SUBAREA_EXTENT": conf.SUBAREA_EXTENT
+                    if conf.SELECT_SUBAREA
+                    else None,
+                },
+                model_ids=model["id"],
+                aggregation=False,
+            )
+
+    # Moist Absolutely Unstable Layer depth (of deepest)
+    if conf.MAUL_DEPTH:
+        for model in models:
+            yield RawRecipe(
+                recipe="depth_of_moist_absolutely_unstable_layers_spatial_plot.yaml",
+                variables={
+                    "MODEL_NAME": model["name"],
+                    "SUBAREA_TYPE": conf.SUBAREA_TYPE if conf.SELECT_SUBAREA else None,
+                    "SUBAREA_EXTENT": conf.SUBAREA_EXTENT
+                    if conf.SELECT_SUBAREA
+                    else None,
+                },
+                model_ids=model["id"],
+                aggregation=False,
+            )
+
+    # Moist Absolutely Unstable Layer base height (of deepest)
+    if conf.MAUL_BASE:
+        for model in models:
+            yield RawRecipe(
+                recipe="base_of_moist_absolutely_unstable_layers_spatial_plot.yaml",
+                variables={
+                    "MODEL_NAME": model["name"],
+                    "SUBAREA_TYPE": conf.SUBAREA_TYPE if conf.SELECT_SUBAREA else None,
+                    "SUBAREA_EXTENT": conf.SUBAREA_EXTENT
+                    if conf.SELECT_SUBAREA
+                    else None,
+                },
+                model_ids=model["id"],
+                aggregation=False,
+            )
+
     # Screen-level temperature probabilities
     for model, condition, threshold in itertools.product(
         models,

src/CSET/operators/__init__.py

@@ -39,6 +39,7 @@
     mesoscale,
     misc,
     plot,
+    precipitation,
     pressure,
     read,
     regrid,
@@ -65,6 +66,7 @@
     "mesoscale",
     "misc",
     "plot",
+    "precipitation",
     "pressure",
     "read",
     "regrid",