[NOT FOR CHANGELOG; PR NOT IN MAIN] Add global mass checks

esmvaltool/recipes/monitor/recipe_create_ranges_cmip6.yml

@@ -192,4 +192,3 @@ diagnostics:
         ensemble: r1i1p1f1
         grid: gr
     scripts: null
-

esmvaltool/recipes/monitor/recipe_create_ranges_obs.yml

@@ -70,9 +70,53 @@ preprocessors:
     climate_statistics:
       operator: max
 
+  pp_min_sum:
+    custom_order: true
+    area_statistics:
+      operator: sum
+    multi_model_statistics:
+      span: full
+      statistics:
+        - operator: min
+    climate_statistics:
+      operator: min
+
+  pp_max_sum:
+    custom_order: true
+    area_statistics:
+      operator: sum
+    multi_model_statistics:
+      span: full
+      statistics:
+        - operator: max
+    climate_statistics:
+      operator: max
 
 diagnostics:
 
+  prw_global:
+    description: Calculate range of reasonable monthly values for global water vapour mass
+    variables:
+      prw_min:
+        short_name: prw
+        preprocessor: pp_min_sum
+        mip: Amon
+      prw_max:
+        short_name: prw
+        preprocessor: pp_max_sum
+        mip: Amon
+    additional_datasets:
+      - {dataset: ERA5, project: native6, type: reanaly, version: v1,
+          tier: 3, start_year: 1995, end_year: 2014}
+      - {dataset: JRA-55, project: OBS6, type: reanaly, version: 1,
+          tier: 2, start_year: 1958, end_year: 2022}
+      - {dataset: MERRA2, project: OBS6, type: reanaly, version: 5.12.4,
+          tier: 3, start_year: 1980, end_year: 2021}
+      - {dataset: ESACCI-WATERVAPOUR, project: OBS6, type: sat, version: CDR2-L3-COMBI-05deg-fv3.1,
+          tier: 3, start_year: 2002, end_year: 2017}
+    scripts: null
+
+
 #  clt:
 #    description: Calculate range of reasonable monthly values for total cloud cover.
 #    variables:

esmvaltool/recipes/monitor/recipe_monitor_sanity_checks.yml

@@ -5,6 +5,7 @@ documentation:
   description: >
     Calculate basic sanity checks e.g. for monitoring new model runs.
   authors:
+    - bock_lisa
     - lauer_axel
   maintainer:
     - lauer_axel
@@ -42,6 +43,27 @@ preprocessors:
     area_statistics:
       operator: max
 
+  global_sum:
+    area_statistics:
+      operator: sum
+
+  global_sum_air_anom:
+    custom_order: true
+    convert_units:
+      units: kg m-2
+    area_statistics:
+      operator: sum
+    anomalies:
+      period: full
+      reference:
+        start_year: 1980
+        start_month: 1
+        start_day: 1
+        end_year: 1980
+        end_month: 12
+        end_day: 31
+
+
 
 diagnostics:
 
@@ -245,3 +267,106 @@ diagnostics:
                 linestyle: '--'
                 linewidth: 1
                 zorder: 3
+
+
+  # *** global sum ***
+
+  prw_mass:
+    description: time series of global water vapour mass
+    variables:
+      prw_sum:
+        <<: *time_period
+        short_name: prw
+        preprocessor: global_sum
+        legend_label: 'global sum'
+    scripts:
+      timeseries:
+        script: monitor/multi_datasets.py
+        plot_folder: '{plot_dir}'
+        group_variables_by: short_name
+        facet_used_for_labels: legend_label
+        plots:
+          timeseries:
+            annual_mean_kwargs: false
+            pyplot_kwargs:
+              title: "Global sum of water vapour"
+            hlines:
+              - y: 14206879936806912.
+                color: red
+                linewidth: 2
+              - y: 10902071180001280.
+                color: red
+                linewidth: 2
+            plot_kwargs:
+              'global sum':
+                color: blue
+                label: null
+                linestyle: '-'
+                linewidth: 2
+                label: '{dataset}'
+                zorder: 3
+
+
+  air_mass:
+    description: time series of global air mass anomaly
+    variables:
+      ps:
+        <<: *time_period
+        preprocessor: global_sum_air_anom
+        legend_label: 'global sum'
+    scripts:
+      timeseries:
+        script: monitor/multi_datasets.py
+        plot_folder: '{plot_dir}'
+        group_variables_by: short_name
+        facet_used_for_labels: legend_label
+        plots:
+          timeseries:
+            annual_mean_kwargs: false
+            pyplot_kwargs:
+              title: "Anomaly of global air mass"
+            hlines:
+              - y: 0.
+                color: black
+                linewidth: 2
+            plot_kwargs:
+              'global sum':
+                color: blue
+                label: null
+                linestyle: '-'
+                linewidth: 2
+                label: '{dataset}'
+                zorder: 3
+
+
+  moisture_flux:
+    description: time series of moisture flux into the atmosphere
+    variables:
+      qep:
+        <<: *time_period
+        derive: true
+        preprocessor: global_sum
+        legend_label: 'global sum'
+    scripts:
+      timeseries:
+        script: monitor/multi_datasets.py
+        plot_folder: '{plot_dir}'
+        group_variables_by: short_name
+        facet_used_for_labels: legend_label
+        plots:
+          timeseries:
+            annual_mean_kwargs: false
+            pyplot_kwargs:
+              title: "Anomaly of global moisture flux into the atmosphere"
+            hlines:
+              - y: 0.
+                color: black
+                linewidth: 2
+            plot_kwargs:
+              'global sum':
+                color: blue
+                label: null
+                linestyle: '-'
+                linewidth: 2
+                label: '{dataset}'
+                zorder: 3