1750 cardington near surface diagnostics

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

@@ -93,20 +93,11 @@ ns=Observations/Cardington
 sort-key=sec-d4
 title=Cardington UK observations
 
-# Air temperature timeseries
-[template variables=CARDINGTON_AIR_TEMPERATURE_SINGLE_POINT_TIME_SERIES]
+# Generic timeseries for variables at various heights
+[template variables=CARDINGTON_SINGLE_POINT_TIME_SERIES]
 ns=Observations/Cardington
-title=Cardington air temperature timeseries
+title=Cardington multi-variable timeseries
 description=Create Cardington diagnostic variable timeseries
 type=python_boolean
 compulsory=true
 sort-key=0card1
-
-# Relative humidity timeseries
-[template variables=CARDINGTON_RELATIVE_HUMIDITY_SINGLE_POINT_TIME_SERIES]
-ns=Observations/Cardington
-title=Cardington relative humidity timeseries
-description=Create Cardington diagnostic variable timeseries
-type=python_boolean
-compulsory=true
-sort-key=0card2

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

@@ -12,8 +12,10 @@ AVIATION_COLOUR_STATE_VISIBILITY=False
 AVIATION_FOG_PRESENCE_DOMAIN_MEAN_TIMESERIES=False
 AVIATION_FOG_PRESENCE_SPATIAL_DIFFERENCE=False
 AVIATION_FOG_PRESENCE_SPATIAL_PLOT=False
-CARDINGTON_AIR_TEMPERATURE_SINGLE_POINT_TIME_SERIES=False
-CARDINGTON_RELATIVE_HUMIDITY_SINGLE_POINT_TIME_SERIES=False
+BASIC_QQ_PLOT=False
+CARDINGTON_AIR_SINGLE_POINT_TIME_SERIES=False
+!!CARDINGTON_AIR_TEMPERATURE_SINGLE_POINT_TIME_SERIES=False
+!!CARDINGTON_RELATIVE_HUMIDITY_SINGLE_POINT_TIME_SERIES=False
 CLOUD_BASE_HEIGHT_LESS_THAN_50_M_DOMAIN_MEAN_TIMESERIES=False
 CLOUD_BASE_HEIGHT_LESS_THAN_50_M_SPATIAL_DIFFERENCE=False
 CLOUD_BASE_HEIGHT_LESS_THAN_50_M_SPATIAL_PLOT=False

src/CSET/loaders/__init__.py

@@ -21,6 +21,7 @@
 
 from CSET.loaders import (
     aoa,
+    cardington,
     histograms,
     observations,
     power_spectrum,
@@ -35,6 +36,7 @@
 
 __all__ = [
     "aoa",
+    "cardington",
     "histograms",
     "observations",
     "power_spectrum",

src/CSET/loaders/cardington.py

@@ -0,0 +1,259 @@
+# © Crown copyright, Met Office (2022-2025) and CSET contributors.
+#
+# Licensed under the Apache License, Version 2.0 (the "License");
+# you may not use this file except in compliance with the License.
+# You may obtain a copy of the License at
+#
+#     http://www.apache.org/licenses/LICENSE-2.0
+#
+# Unless required by applicable law or agreed to in writing, software
+# distributed under the License is distributed on an "AS IS" BASIS,
+# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+# See the License for the specific language governing permissions and
+# limitations under the License.
+
+"""Load cardington recipes."""
+# long_name='VIS AT 1.5M  (incl precip)         M'
+#'SURFACE SENSIBLE  HEAT FLUX     W/M2'
+#'SURFACE LATENT HEAT FLUX        W/M2'
+#'DEWPOINT AT 1.5M (K)'
+
+from CSET.recipes import Config, RawRecipe, get_models
+
+
+def height_token(ht: float) -> str:
+    """Convert height to the token used in Cardington variable names."""
+    # Cardington uses "1p2m" for 1.2 m
+    if abs(ht - 1.2) < 1e-9:
+        return "1p2m"
+    # Integers become "2m", "10m", etc.
+    if float(ht).is_integer():
+        return f"{int(ht)}m"
+    # Fallback for any other non-integer heights (e.g. 0.5 -> "0p5m")
+    s = str(ht).replace(".", "p")
+    return f"{s}m"
+
+
+def load(conf: Config):
+    """Yield recipes from the given workflow configuration."""
+    # Load a list of model detail dictionaries.
+    models = get_models(conf.asdict())
+    T30 = "30min"
+    T05 = "05min"
+    if conf.CARDINGTON_SINGLE_POINT_TIME_SERIES:
+        cardington_models = [m for m in models if "Cardington" in m["name"]]
+        um_model = next(m for m in models if "Cardington" not in m["name"])
+
+        fields = [("visibility", "m01s03i281", "visibility_belfort", [2])]
+        for plot_field, model_field, card_field, heights in fields:
+            for ht in heights:
+                card_label_30 = f"Cardington {T30} ({ht} m)"
+                card_label_05 = f"Cardington {T05} ({ht} m)"
+                varname = None
+                if isinstance(card_field, str):
+                    varname = f"{card_field}_{height_token(ht)}"
+                    yield RawRecipe(
+                        recipe="cardington_visibility_single_point_time_series.yaml",
+                        variables={
+                            "MODEL_NAME": um_model["name"],
+                            "VARNAME": varname,
+                            "CARDINGTON_30_LABEL": card_label_30,
+                            "CARDINGTON_05_LABEL": card_label_05,
+                            "UM_VARNAME": model_field,
+                            "HEIGHT": ht,
+                            "PLOTNAME": plot_field,
+                        },
+                        model_ids=[
+                            cardington_models[0]["id"],
+                            cardington_models[1]["id"],
+                            um_model["id"],
+                        ],
+                        aggregation=False,
+                    )
+
+        field = "surface_temperature"
+        card_label_30 = f"Cardington {T30}"
+        card_label_05 = f"Cardington {T05}"
+        yield RawRecipe(
+            recipe="cardington_surface_temperature_time_series.yaml",
+            variables={
+                "MODEL_NAME": um_model["name"],
+                "VARNAME": field,
+                "CARDINGTON_30_LABEL": card_label_30,
+                "CARDINGTON_05_LABEL": card_label_05,
+                "UM_VARNAME": "m01s00i024",
+                "PLOTNAME": field,
+            },
+            model_ids=[
+                cardington_models[0]["id"],
+                cardington_models[1]["id"],
+                um_model["id"],
+            ],
+            aggregation=False,
+        )
+
+        fields = [("surface_air_pressure", "air_pressure")]
+        for model_field, card_field in fields:
+            card_label_30 = f"Cardington {T30}"
+            card_label_05 = f"Cardington {T05}"
+            yield RawRecipe(
+                recipe="cardington_surface_pressure_time_series.yaml",
+                variables={
+                    "MODEL_NAME": um_model["name"],
+                    "VARNAME": card_field,
+                    "CARDINGTON_30_LABEL": card_label_30,
+                    "CARDINGTON_05_LABEL": card_label_05,
+                    "UM_VARNAME": model_field,
+                    "PLOTNAME": card_field,
+                },
+                model_ids=[
+                    cardington_models[0]["id"],
+                    cardington_models[1]["id"],
+                    um_model["id"],
+                ],
+                aggregation=False,
+            )
+
+        fields = [
+            ("wind_speed", "wind_speed_vector", [2, 10, 25, 50]),
+            ("wind_direction", "wind_direction", [2, 10, 25, 50]),
+        ]
+        for plot_field, card_field, heights in fields:
+            for ht in heights:
+                card_label_30 = f"Cardington {T30} ({ht} m)"
+                card_label_05 = f"Cardington {T05} ({ht} m)"
+                yield RawRecipe(
+                    recipe="cardington_wind_single_point_time_series.yaml",
+                    variables={
+                        "MODEL_NAME": um_model["name"],
+                        "VARNAME": f"{card_field}_{height_token(ht)}",
+                        "CARDINGTON_30_LABEL": card_label_30,
+                        "CARDINGTON_05_LABEL": card_label_05,
+                        "UM_U_VARNAME": "x_wind",
+                        "UM_V_VARNAME": "y_wind",
+                        "PLOT_VARNAME": plot_field,
+                        "HEIGHT": ht,
+                    },
+                    model_ids=[
+                        cardington_models[0]["id"],
+                        cardington_models[1]["id"],
+                        um_model["id"],
+                    ],
+                    aggregation=False,
+                )
+
+        fields = [
+            (
+                "air temperature",
+                "air_temperature",
+                "air_temperature_rtd",
+                [1.2, 10, 25, 50],
+            ),
+            (
+                "relative humidity",
+                "relative_humidity",
+                "relative_humidity",
+                [1.2, 10, 25, 50],
+            ),
+            ("dew point", "m01s03i250", "dewpoint_temperature", [1.2, 10, 25, 50]),
+            ("wind gust", "m01s03i463", "wind_speed_max", [2, 10, 25, 50]),
+            ("latent heat flux", "m01s03i234", "wq_covariance", [10]),
+            (
+                "sensible heat flux",
+                "m01s03i217",
+                {
+                    "wt_covariance": "HEIGHT",
+                    "air_temperature_rtd": "HEIGHT_EXCEPT_2M",
+                    "pressure_barometric": "FIXED",
+                },
+                [2, 10, 25, 50],
+            ),
+        ]
+        for plot_field, model_field, card_field, heights in fields:
+            for ht in heights:
+                card_label_30 = f"Cardington {T30} ({ht} m)"
+                card_label_05 = f"Cardington {T05} ({ht} m)"
+                varname = None
+                if isinstance(card_field, str):
+                    varname = f"{card_field}_{height_token(ht)}"
+
+                if "latent" in plot_field:
+                    yield RawRecipe(
+                        recipe="cardington_latent_heat_single_point_time_series.yaml",
+                        variables={
+                            "MODEL_NAME": um_model["name"],
+                            "VARNAME": varname,
+                            "CARDINGTON_30_LABEL": card_label_30,
+                            "UM_VARNAME": model_field,
+                            "HEIGHT": ht,
+                            "PLOTNAME": plot_field,
+                        },
+                        model_ids=[cardington_models[1]["id"], um_model["id"]],
+                        aggregation=False,
+                    )
+                elif "sensible" in plot_field and isinstance(card_field, dict):
+                    card_varnames = []
+                    for var, rule in card_field.items():
+                        if rule == "HEIGHT":
+                            card_varnames.append(f"{var}_{height_token(ht)}")
+                        elif rule == "HEIGHT_EXCEPT_2M":
+                            if ht == 2:
+                                card_varnames.append(f"{var}_1p2m")
+                            else:
+                                card_varnames.append(f"{var}_{height_token(ht)}")
+                        elif rule == "FIXED":
+                            card_varnames.append(var)
+                        else:
+                            card_varnames.append(f"{var}_{rule}")
+
+                    yield RawRecipe(
+                        recipe="cardington_sensible_heat_single_point_time_series.yaml",
+                        variables={
+                            "MODEL_NAME": um_model["name"],
+                            "CARDINGTON_30_LABEL": card_label_30,
+                            "UM_VARNAME": model_field,  # m01s03i217 (unchanged)
+                            "HEIGHT": ht,
+                            "PLOTNAME": plot_field,
+                            "CARDINGTON_VARNAMES": ",".join(card_varnames),
+                        },
+                        model_ids=[cardington_models[1]["id"], um_model["id"]],
+                        aggregation=False,
+                    )
+                elif "temperature" in model_field or "temperature" in card_field:
+                    yield RawRecipe(
+                        recipe="cardington_temperature_single_point_time_series.yaml",
+                        variables={
+                            "MODEL_NAME": um_model["name"],
+                            "VARNAME": varname,
+                            "CARDINGTON_30_LABEL": card_label_30,
+                            "CARDINGTON_05_LABEL": card_label_05,
+                            "UM_VARNAME": model_field,
+                            "HEIGHT": ht,
+                            "PLOTNAME": plot_field,
+                        },
+                        model_ids=[
+                            cardington_models[0]["id"],
+                            cardington_models[1]["id"],
+                            um_model["id"],
+                        ],
+                        aggregation=False,
+                    )
+                elif "gust" in plot_field:
+                    yield RawRecipe(
+                        recipe="cardington_wind_gust_single_point_time_series.yaml",
+                        variables={
+                            "MODEL_NAME": um_model["name"],
+                            "VARNAME": varname,
+                            "CARDINGTON_30_LABEL": card_label_30,
+                            "CARDINGTON_05_LABEL": card_label_05,
+                            "UM_VARNAME": model_field,
+                            "HEIGHT": ht,
+                            "PLOTNAME": plot_field,
+                        },
+                        model_ids=[
+                            cardington_models[0]["id"],
+                            cardington_models[1]["id"],
+                            um_model["id"],
+                        ],
+                        aggregation=False,
+                    )

src/CSET/loaders/timeseries.py

@@ -368,31 +368,6 @@ def load(conf: Config):
             aggregation=False,
         )
 
-    # Cardington air temperature time series.
-    if conf.CARDINGTON_AIR_TEMPERATURE_SINGLE_POINT_TIME_SERIES:
-        base_model = models[0]
-        for model in models[1:]:
-            yield RawRecipe(
-                recipe="cardington_air_temperature_single_point_time_series.yaml",
-                variables={
-                    "MODEL_NAME": model["name"],
-                },
-                model_ids=[base_model["id"], model["id"]],
-                aggregation=False,
-            )
-
-    # Cardington relative humidity time series.
-    if conf.CARDINGTON_RELATIVE_HUMIDITY_SINGLE_POINT_TIME_SERIES:
-        base_model = models[0]
-        for model in models[1:]:
-            yield RawRecipe(
-                recipe="cardington_relative_humidity_single_point_time_series.yaml",
-                variables={
-                    "MODEL_NAME": model["name"],
-                },
-                model_ids=[base_model["id"], model["id"]],
-                aggregation=False,
-            )
     # Surface wind gusts on Beaufort Scale
     if conf.SFC_WIND_GUSTS_BEAUFORT_SCALE_DOMAIN_MEAN_TIMESERIES:
         yield RawRecipe(

src/CSET/operators/_colorbar_definition.json

@@ -273,8 +273,8 @@
     "cmap": "bwr",
     "max": 250.0,
     "min": -250.0,
-    "ymax": 250.0,
-    "ymin": -50.0
+    "ymax": "auto",
+    "ymin": "auto"
   },
   "grid_surface_upward_sensible_heat_flux_difference": {
     "cmap": "PuOr",
@@ -1030,7 +1030,7 @@
   },
   "visibility_in_air": {
     "cmap": "viridis",
-    "max": 30000.0,
+    "max": "auto",
     "min": 0.0
   },
   "visibility_in_air_difference": {