Add Hovmoeller Z vs. time plot to monitoring diagnostic (for Changelo…

…g v2.10: authors: @cubeme and @helgehr) (#3345) Co-authored-by: Helge Heuer <[email protected]> Co-authored-by: helgehr <[email protected]> Co-authored-by: Manuel Schlund <[email protected]> Co-authored-by: Manuel Schlund <[email protected]>
ESMValGroup · Sep 7, 2023 · 41acdb3 · 41acdb3
1 parent f9c7e4a
commit 41acdb3
Show file tree

Hide file tree

Showing 7 changed files with 497 additions and 27 deletions.
diff --git a/.zenodo.json b/.zenodo.json
@@ -181,6 +181,11 @@
             "name": "Hassler, Birgit",
             "orcid": "0000-0003-2724-709X"
         },
+        {
+            "affiliation": "DLR, Germany",
+            "name": "Heuer, Helge",
+            "orcid": "0000-0003-2411-7150"
+        },
         {
             "affiliation": "BSC, Spain",
             "name": "Hunter, Alasdair",
@@ -199,6 +204,10 @@
             "affiliation": "MPI for Biogeochemistry, Germany",
             "name": "Koirala, Sujan"
         },
+        {
+            "affiliation": "DLR, Germany",
+            "name": "Kuehbacher, Birgit"
+        },
         {
             "affiliation": "BSC, Spain",
             "name": "Lledó, Llorenç"

diff --git a/CITATION.cff b/CITATION.cff
@@ -181,6 +181,11 @@ authors:
     family-names: Hassler
     given-names: Birgit
     orcid: "https://orcid.org/0000-0003-2724-709X"
+  -
+    affiliation: "DLR, Germany"
+    family-names: Heuer
+    given-names: Helge
+    orcid: "https://orcid.org/0000-0003-2411-7150"
   -
     affiliation: "BSC, Spain"
     family-names: Hunter
@@ -199,6 +204,10 @@ authors:
     affiliation: "MPI for Biogeochemistry, Germany"
     family-names: Koirala
     given-names: Sujan
+  -
+    affiliation: "DLR, Germany"
+    family-names: Kuehbacher
+    given-names: Birgit
   -
     affiliation: "BSC, Spain"
     family-names: Lledó

diff --git a/doc/sphinx/source/recipes/figures/monitor/hovmoeller_z_vs_time_with_ref.png b/doc/sphinx/source/recipes/figures/monitor/hovmoeller_z_vs_time_with_ref.png
diff --git a/doc/sphinx/source/recipes/recipe_monitor.rst b/doc/sphinx/source/recipes/recipe_monitor.rst
@@ -216,4 +216,11 @@ Zonal mean profile of ta including a reference dataset.
    :align:   center
    :width:   14cm
 
-1D profile of pr over latitude.
+Zonal mean pr including a reference dataset.
+
+.. _fig_hovmoeller_z_vs_time_with_ref:
+.. figure::  /recipes/figures/monitor/hovmoeller_z_vs_time_with_ref.png
+   :align:   center
+   :width:   14cm
+
+Hovmoeller plot (pressure vs time) of ta including a reference dataset.
diff --git a/esmvaltool/config-references.yml b/esmvaltool/config-references.yml
@@ -257,6 +257,11 @@ authors:
     name: Hempelmann, Nils
     institute: IPSL, France
     orcid:
+  heuer_helge:
+    name: Heuer, Helge
+    institute: DLR, Germany
+    email: [email protected]
+    orcid: https://orcid.org/0000-0003-2411-7150
   hogan_emma:
     name: Hogan, Emma
     institute: MetOffice, UK
@@ -297,6 +302,11 @@ authors:
     name: Krasting, John
     institute: NOAA, USA
     orcid: https://orcid.org/0000-0002-4650-9844
+  kuehbacher_birgit:
+    name: Kuehbacher, Birgit
+    institute: DLR, Germany
+    email: [email protected]
+    orcid:
   lejeune_quentin:
     name: Lejeune, Quentin
     institute: Climate Analytics, Germany

diff --git a/esmvaltool/diag_scripts/monitor/multi_datasets.py b/esmvaltool/diag_scripts/monitor/multi_datasets.py
diff --git a/esmvaltool/recipes/monitor/recipe_monitor_with_refs.yml b/esmvaltool/recipes/monitor/recipe_monitor_with_refs.yml
@@ -8,14 +8,17 @@ documentation:
   authors:
     - schlund_manuel
     - winterstein_franziska
+    - sarauer_ellen
+    - kuehbacher_birgit
+    - heuer_helge
   maintainer:
     - schlund_manuel
 
 
 datasets:
   # Note: plot_label currently only used by diagnostic plot_multiple_annual_cycles
-  - {project: CMIP6, dataset: EC-Earth3, exp: historical, ensemble: r1i1p1f1, grid: gr, plot_label: 'EC-Earth3 historical'}
-  - {project: CMIP6, dataset: CanESM5,   exp: historical, ensemble: r1i1p1f1, grid: gn, plot_label: 'Reference (CanESM5 historical)', reference_for_monitor_diags: true}
+  - {project: CMIP6, dataset: MPI-ESM1-2-HR, exp: historical, ensemble: r1i1p1f1, grid: gn, plot_label: 'MPI-ESM1-2-HR historical'}
+  - {project: CMIP6, dataset: MPI-ESM1-2-LR, exp: historical, ensemble: r1i1p1f1, grid: gn, plot_label: 'Reference (MPI-ESM1-2-LR historical)', reference_for_monitor_diags: true}
 
 preprocessors:
 
@@ -71,10 +74,23 @@ preprocessors:
       operator: mean
     regrid:
       target_grid: 2x2
-      scheme:
-        reference: esmf_regrid.schemes:ESMFAreaWeighted
+      scheme: linear
     zonal_statistics:
       operator: mean
+    convert_units:
+      units: mm day-1
+
+  global_mean_extract_levels:
+    custom_order: true
+    extract_levels:
+      levels: {cmor_table: CMIP6, coordinate: alt16}
+      scheme: linear
+      coordinate: altitude
+    regrid:
+      target_grid: 2x2
+      scheme: linear
+    area_statistics:
+      operator: mean
 
 
 diagnostics:
@@ -96,9 +112,9 @@ diagnostics:
             annual_mean_kwargs:
               linestyle: '--'
             plot_kwargs:
-              EC-Earth3:  # = dataset since 'facet_used_for_labels' is 'dataset' by default
+              MPI-ESM1-2-HR:  # = dataset since 'facet_used_for_labels' is 'dataset' by default
                 color: C0
-              CanESM5:
+              MPI-ESM1-2-LR:
                 color: black
 
   plot_multiple_annual_cycles:
@@ -117,9 +133,9 @@ diagnostics:
             legend_kwargs:
               loc: upper right
             plot_kwargs:
-              'EC-Earth3 historical':  # = plot_label since 'facet_used_for_labels: plot_label'
+              'MPI-ESM1-2-HR historical':  # = plot_label since 'facet_used_for_labels: plot_label'
                 color: C0
-              'Reference (CanESM5 historical)':
+              'Reference (MPI-ESM1-2-LR historical)':
                 color: black
             pyplot_kwargs:
               title: Near-Surface Air Temperature on Northern Hemisphere
@@ -173,9 +189,9 @@ diagnostics:
         plots:
           1d_profile:
             plot_kwargs:
-              EC-Earth3:  # = dataset since 'facet_used_for_labels' is 'dataset' by default
+              MPI-ESM1-2-HR:  # = dataset since 'facet_used_for_labels' is 'dataset' by default
                 color: C0
-              CanESM5:
+              MPI-ESM1-2-LR:
                 color: black
 
   plot_variable_vs_latitude:
@@ -187,6 +203,27 @@ diagnostics:
         timerange: '20000101/20030101'
     scripts:
       plot:
+        <<: *plot_multi_dataset_default
         script: monitor/multi_datasets.py
         plots:
           variable_vs_lat:
+
+  plot_hovmoeller_z_vs_time:
+    description: Plot Hovmoeller Z vs. time including reference datasets.
+    variables:
+      ta:
+        preprocessor: global_mean_extract_levels
+        mip: Amon
+        timerange: '2000/2004'
+    scripts:
+      plot:
+        <<: *plot_multi_dataset_default
+        script: monitor/multi_datasets.py
+        plots:
+          hovmoeller_z_vs_time:
+            plot_func: contourf
+            common_cbar: true
+            time_format: '%Y'
+            log_y: false
+            pyplot_kwargs:
+              ylim: [0, 20000]