Regional historical changes for REF (#3926)

Co-authored-by: Axel Lauer <axel.lauer@dlr.de>
Co-authored-by: Diego Cammarano <diegocamm@gmail.com>
Co-authored-by: Lisa Bock <lisa.bock@dlr.de>
Co-authored-by: Manuel Schlund <32543114+schlunma@users.noreply.github.com>
Co-authored-by: Birgit Hassler <birgit.hassler@dlr.de>
Co-authored-by: Alistair Sellar <16133375+alistairsellar@users.noreply.github.com>
Co-authored-by: Bettina Gier <gier@uni-bremen.de>
Co-authored-by: Manuel Schlund <manuel.schlund@dlr.de>
Co-authored-by: Bouwe Andela <b.andela@esciencecenter.nl>

Author: 9 people

doc/sphinx/source/recipes/figures/ref/timeseries_tas_ambiguous_dataset_Amon_historical_r1i1p1f1.png

@@ -13,7 +13,7 @@ Here ESMValTool recipes are collected which will be used in the CMIP
 Available recipes
 -----------------
 
-Recipes are stored in recipes
+Recipes are stored in `recipes/`
 
 * :ref:`recipe_ecs.yml <recipes_ecs>`:
   Calculate equilibrium climate sensitivity (ECS)
@@ -23,6 +23,12 @@ Recipes are stored in recipes
   Calculate transient climate response to cumulative CO2 emissions (TCRE)
 * ref/recipe_ref_cre.yml:
   Maps and zonal means of longwave and shortwave cloud radiative effect
+* ref/recipe_ref_timeseries_region.yml:
+  Time series in comparison with reference data for a selected IPCC region defined through a shape file, based on :ref:`recipe_ref_timeseries.yml <recipe_benchmarking>`
+* ref/recipe_ref_annual_cycle_region.yml:
+  Annual cycle in comparison with reference data for a selected IPCC region defined through a shape file, based on :ref:`recipe_ref_annual_cycle.yml <recipe_benchmarking>`
+* ref/recipe_ref_trend_regions.yml:
+  Linear Trends for all IPCC land regions compared with reference data, based on :ref:`recipe_ref_trend_regions.yml <recipes_seaborn_diag>`
 * ref/recipe_ref_scatterplot.yml:
   2D histograms with focus on clouds
 * ref/recipe_ref_sea_ice_area_basic.yml:
@@ -82,3 +88,11 @@ Example plots:
 
    Zonal mean vertically resolved ozone climatology from CNRM-ESM2-1 compared with ESACCI-OZONE
    for the years 1990 to 2000. Created with recipe_ref_ozone.yml.
+
+.. _fig_ref_7:
+.. figure::  /recipes/figures/ref/timeseries_tas_ambiguous_dataset_Amon_historical_r1i1p1f1.png
+   :align:   center
+   :width:   8cm
+
+   Time series of near-surface air temperature anomalies from MIROC6 compared with HadCRUT5
+   for N.Europe for the years 1980 to 2014 (reference period 1980 to 2009). Created with recipe_ref_timeseries_region.yml.

doc/sphinx/source/recipes/recipe_ref.rst

@@ -0,0 +1,149 @@
+# ESMValTool
+---
+documentation:
+  title: Benchmarking of a single model for a region through extract_shape.
+  description: >
+    Benchmarking: annual cycle.
+  authors:
+    - lauer_axel
+    - bock_lisa
+    - hassler_birgit
+    - lindenlaub_lukas
+    - schlund_manuel
+    - weigel_katja
+  maintainer:
+    - lauer_axel
+  references:
+    - lauer25gmd
+  projects:
+    - dlrmabak
+
+
+datasets:
+  # Dataset to be benchmarked
+  - {dataset: MIROC6, grid: gn, benchmark_dataset: true, plot_label: 'MIROC6'}
+
+preprocessors:
+  pp_abs:
+    regrid:
+      target_grid: 1x1
+      scheme: linear
+    extract_shape:
+      shapefile: ar6
+      crop: true
+      decomposed: true
+      ids:
+        Name:
+          - N.Europe
+    climate_statistics:
+      period: monthly
+    area_statistics:
+      operator: mean
+
+  pp_abs_200:
+    extract_levels:
+      levels: 20000
+      scheme: linear
+      coordinate: air_pressure
+    regrid:
+      target_grid: 1x1
+      scheme: linear
+    extract_shape:
+      shapefile: ar6
+      crop: true
+      decomposed: true
+      ids:
+        Name:
+          - N.Europe
+    climate_statistics:
+      period: monthly
+    area_statistics:
+      operator: mean
+
+
+diagnostics:
+  anncyc:
+    description: Plot global mean mean annual cycle for the given period including a reference dataset.
+    variables:
+      tas: &var_default_abs
+        start_year: 1980
+        end_year: 2009
+        preprocessor: pp_abs
+        project: CMIP6
+        mip: Amon
+        exp: historical
+        ensemble: r1i1p1f1
+        grid: gr
+        reference_dataset: HadCRUT5
+        additional_datasets:
+          - {dataset: HadCRUT5, project: OBS, type: ground, version: 5.0.1.0-analysis,
+             tier: 2, reference_for_metric: true'}
+      pr:
+        <<: *var_default_abs
+        reference_dataset: GPCP-V2.3
+        additional_datasets:
+          - {dataset: GPCP-V2.3, project: obs4MIPs, tier: 1,
+             reference_for_metric: true'}
+      psl:
+        <<: *var_default_abs
+        reference_dataset: ERA-5
+        additional_datasets:
+          - {dataset: ERA-5, project: obs4MIPs, tier: 1,
+             reference_for_metric: true'}
+      hus_200:
+        <<: *var_default_abs
+        short_name: hus
+        reference_dataset: ERA5
+        preprocessor: pp_abs_200
+        additional_datasets:
+          - {dataset: ERA5, project: native6, type: reanaly, version: v1,
+             tier: 3, reference_for_metric: true'}
+      ua_200:
+        <<: *var_default_abs
+        short_name: ua
+        reference_dataset: ERA-5
+        preprocessor: pp_abs_200
+        additional_datasets:
+          - {dataset: ERA-5, project: obs4MIPs, tier: 1,
+             reference_for_metric: true'}
+    scripts:
+      allplots:
+        script: monitor/multi_datasets.py
+        plot_folder: '{plot_dir}'
+        group_variables_by: variable_group
+        plots:
+          annual_cycle:
+            annual_mean_kwargs: False
+            plot_kwargs:
+              default:
+                color: red
+                label: '{dataset}'
+                linestyle: '-'
+                linewidth: 2
+                zorder: 4
+              HadCRUT5:
+                color: black
+                label: '{dataset}'
+                linestyle: '-'
+                linewidth: 2
+                zorder: 3
+              ERA5:
+                color: black
+                label: '{dataset}'
+                linestyle: '-'
+                linewidth: 2
+                zorder: 3
+              ERA-5:
+                color: black
+                label: '{dataset}'
+                linestyle: '-'
+                linewidth: 2
+                zorder: 3
+              GPCP-V2.3:
+                color: black
+                label: '{dataset}'
+                linestyle: '-'
+                linewidth: 2
+                zorder: 3
+            pyplot_kwargs:
+              title: '{long_name} N.Europe'