python/cleanup_ETHZ_ClimWIP_atlas_netcdf.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Cleaning up Atlas data - ETHZ ClimWIP\n",
    "**Function**      : Preprocess netCDF files and restructure the dataset<br>\n",
    "**Description**   : In this notebook serves to clean up Atlas data which is given in netcdf format and aggregate the data into a single file.<br>\n",
    "**Return Values   : .nc files**<br>\n",
    "**Note**          : All the data is saved to netCDF4 format. Note that data from different models may vary concerning the resolution and coordinates.<br>"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [],
   "source": [
    "import os\n",
    "from pathlib import Path\n",
    "\n",
    "import numpy as np\n",
    "import pandas as pd\n",
    "import xarray as xr"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Path\n",
    "Specify the path to the dataset and the place to save the outputs. <br>"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [],
   "source": [
    "# please specify data path\n",
    "datapath = Path(\"./AtlasData/raw\")\n",
    "\n",
    "# please specify output path\n",
    "output_path = Path(\"./AtlasData/preprocess\")\n",
    "os.makedirs(output_path, exist_ok=True)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Components used to create the output file names. Here, only `institution_id` and `cmor_var` is based on on CMIP DRS conventions."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "output_file_name = {\n",
    "    \"prefix\": \"atlas\",\n",
    "    \"activity\": \"EUCP\",  # project name e.g. EUCP\n",
    "    \"institution_id\": \"ETHZ\",  # ETHZ\n",
    "    \"source\": \"CMIP6\",  # e.g. CMIP6 or CMIP5\n",
    "    \"method\": \"ClimWIP\",  # e.g. ClimWIP\n",
    "    \"sub_method\": \"cons\",  # e.g. cons or uncons\n",
    "    \"cmor_var\": \"tas\",  # e.g. tas or pr\n",
    "}"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Load and process raw data"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Make some functions to combining multiple dimensions with a preprocessor and load data"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [],
   "source": [
    "INSTITUTION_ID = output_file_name[\"institution_id\"]\n",
    "METHOD = output_file_name[\"method\"]\n",
    "\n",
    "\n",
    "def add_percentile_tas(ds):\n",
    "    filename = ds.encoding[\"source\"]\n",
    "    percentile = int(filename.split(\"_\")[-1][:2])\n",
    "    _, variable, future, season, dataset, _ = filename.split(\"_\")[1:]\n",
    "\n",
    "    return ds.assign_coords(percentile=percentile).expand_dims(\"percentile\")\n",
    "\n",
    "\n",
    "def add_percentile_pr(ds):\n",
    "    filename = ds.encoding[\"source\"]\n",
    "    percentile = int(filename.split(\"_\")[-2][:2])\n",
    "    _, variable, future, season, dataset, _ = filename.split(\"_\")[1:-1]\n",
    "\n",
    "    return ds.assign_coords(percentile=percentile).expand_dims(\"percentile\")\n",
    "\n",
    "\n",
    "# data loader and batch processing\n",
    "def load_data(project, season, variable):\n",
    "    # open multiple files with xarray\n",
    "    if variable == \"tas\":\n",
    "        ds = xr.open_mfdataset(\n",
    "            str(\n",
    "                Path(\n",
    "                    datapath,\n",
    "                    f\"{INSTITUTION_ID}_{METHOD}\",\n",
    "                    f\"eur_{variable}_41-60_{season}_{project}_*perc.nc\",\n",
    "                )\n",
    "            ),\n",
    "            preprocess=add_percentile_tas,\n",
    "        )\n",
    "    elif variable == \"pr\":\n",
    "        ds = xr.open_mfdataset(\n",
    "            str(\n",
    "                Path(\n",
    "                    datapath,\n",
    "                    f\"{INSTITUTION_ID}_{METHOD}\",\n",
    "                    f\"eur_{variable}_41-60_{season}_{project}_*perc_rel.nc\",\n",
    "                )\n",
    "            ),\n",
    "            preprocess=add_percentile_pr,\n",
    "        )\n",
    "    else:\n",
    "        raise ValueError(\"Given variable is not valid.\")\n",
    "\n",
    "    weighted = (\n",
    "        ds[f\"{variable}_mean_weighted\"]\n",
    "        .rename(variable)\n",
    "        .assign_coords(constrained=1)\n",
    "        .expand_dims(\"constrained\")\n",
    "    )\n",
    "    unweighted = (\n",
    "        ds[f\"{variable}_mean\"]\n",
    "        .rename(variable)\n",
    "        .assign_coords(constrained=0)\n",
    "        .expand_dims(\"constrained\")\n",
    "    )\n",
    "    return xr.concat([weighted, unweighted], dim=\"constrained\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Call functions"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div><svg style=\"position: absolute; width: 0; height: 0; overflow: hidden\">\n",
       "<defs>\n",
       "<symbol id=\"icon-database\" viewBox=\"0 0 32 32\">\n",
       "<path d=\"M16 0c-8.837 0-16 2.239-16 5v4c0 2.761 7.163 5 16 5s16-2.239 16-5v-4c0-2.761-7.163-5-16-5z\"></path>\n",
       "<path d=\"M16 17c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n",
       "<path d=\"M16 26c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n",
       "</symbol>\n",
       "<symbol id=\"icon-file-text2\" viewBox=\"0 0 32 32\">\n",
       "<path d=\"M28.681 7.159c-0.694-0.947-1.662-2.053-2.724-3.116s-2.169-2.030-3.116-2.724c-1.612-1.182-2.393-1.319-2.841-1.319h-15.5c-1.378 0-2.5 1.121-2.5 2.5v27c0 1.378 1.122 2.5 2.5 2.5h23c1.378 0 2.5-1.122 2.5-2.5v-19.5c0-0.448-0.137-1.23-1.319-2.841zM24.543 5.457c0.959 0.959 1.712 1.825 2.268 2.543h-4.811v-4.811c0.718 0.556 1.584 1.309 2.543 2.268zM28 29.5c0 0.271-0.229 0.5-0.5 0.5h-23c-0.271 0-0.5-0.229-0.5-0.5v-27c0-0.271 0.229-0.5 0.5-0.5 0 0 15.499-0 15.5 0v7c0 0.552 0.448 1 1 1h7v19.5z\"></path>\n",
       "<path d=\"M23 26h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "<path d=\"M23 22h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "<path d=\"M23 18h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "</symbol>\n",
       "</defs>\n",
       "</svg>\n",
       "<style>/* CSS stylesheet for displaying xarray objects in jupyterlab.\n",
       " *\n",
       " */\n",
       "\n",
       ":root {\n",
       "  --xr-font-color0: var(--jp-content-font-color0, rgba(0, 0, 0, 1));\n",
       "  --xr-font-color2: var(--jp-content-font-color2, rgba(0, 0, 0, 0.54));\n",
       "  --xr-font-color3: var(--jp-content-font-color3, rgba(0, 0, 0, 0.38));\n",
       "  --xr-border-color: var(--jp-border-color2, #e0e0e0);\n",
       "  --xr-disabled-color: var(--jp-layout-color3, #bdbdbd);\n",
       "  --xr-background-color: var(--jp-layout-color0, white);\n",
       "  --xr-background-color-row-even: var(--jp-layout-color1, white);\n",
       "  --xr-background-color-row-odd: var(--jp-layout-color2, #eeeeee);\n",
       "}\n",
       "\n",
       "html[theme=dark],\n",
       "body.vscode-dark {\n",
       "  --xr-font-color0: rgba(255, 255, 255, 1);\n",
       "  --xr-font-color2: rgba(255, 255, 255, 0.54);\n",
       "  --xr-font-color3: rgba(255, 255, 255, 0.38);\n",
       "  --xr-border-color: #1F1F1F;\n",
       "  --xr-disabled-color: #515151;\n",
       "  --xr-background-color: #111111;\n",
       "  --xr-background-color-row-even: #111111;\n",
       "  --xr-background-color-row-odd: #313131;\n",
       "}\n",
       "\n",
       ".xr-wrap {\n",
       "  display: block;\n",
       "  min-width: 300px;\n",
       "  max-width: 700px;\n",
       "}\n",
       "\n",
       ".xr-text-repr-fallback {\n",
       "  /* fallback to plain text repr when CSS is not injected (untrusted notebook) */\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-header {\n",
       "  padding-top: 6px;\n",
       "  padding-bottom: 6px;\n",
       "  margin-bottom: 4px;\n",
       "  border-bottom: solid 1px var(--xr-border-color);\n",
       "}\n",
       "\n",
       ".xr-header > div,\n",
       ".xr-header > ul {\n",
       "  display: inline;\n",
       "  margin-top: 0;\n",
       "  margin-bottom: 0;\n",
       "}\n",
       "\n",
       ".xr-obj-type,\n",
       ".xr-array-name {\n",
       "  margin-left: 2px;\n",
       "  margin-right: 10px;\n",
       "}\n",
       "\n",
       ".xr-obj-type {\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-sections {\n",
       "  padding-left: 0 !important;\n",
       "  display: grid;\n",
       "  grid-template-columns: 150px auto auto 1fr 20px 20px;\n",
       "}\n",
       "\n",
       ".xr-section-item {\n",
       "  display: contents;\n",
       "}\n",
       "\n",
       ".xr-section-item input {\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-section-item input + label {\n",
       "  color: var(--xr-disabled-color);\n",
       "}\n",
       "\n",
       ".xr-section-item input:enabled + label {\n",
       "  cursor: pointer;\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-section-item input:enabled + label:hover {\n",
       "  color: var(--xr-font-color0);\n",
       "}\n",
       "\n",
       ".xr-section-summary {\n",
       "  grid-column: 1;\n",
       "  color: var(--xr-font-color2);\n",
       "  font-weight: 500;\n",
       "}\n",
       "\n",
       ".xr-section-summary > span {\n",
       "  display: inline-block;\n",
       "  padding-left: 0.5em;\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:disabled + label {\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-section-summary-in + label:before {\n",
       "  display: inline-block;\n",
       "  content: '►';\n",
       "  font-size: 11px;\n",
       "  width: 15px;\n",
       "  text-align: center;\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:disabled + label:before {\n",
       "  color: var(--xr-disabled-color);\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:checked + label:before {\n",
       "  content: '▼';\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:checked + label > span {\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-section-summary,\n",
       ".xr-section-inline-details {\n",
       "  padding-top: 4px;\n",
       "  padding-bottom: 4px;\n",
       "}\n",
       "\n",
       ".xr-section-inline-details {\n",
       "  grid-column: 2 / -1;\n",
       "}\n",
       "\n",
       ".xr-section-details {\n",
       "  display: none;\n",
       "  grid-column: 1 / -1;\n",
       "  margin-bottom: 5px;\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:checked ~ .xr-section-details {\n",
       "  display: contents;\n",
       "}\n",
       "\n",
       ".xr-array-wrap {\n",
       "  grid-column: 1 / -1;\n",
       "  display: grid;\n",
       "  grid-template-columns: 20px auto;\n",
       "}\n",
       "\n",
       ".xr-array-wrap > label {\n",
       "  grid-column: 1;\n",
       "  vertical-align: top;\n",
       "}\n",
       "\n",
       ".xr-preview {\n",
       "  color: var(--xr-font-color3);\n",
       "}\n",
       "\n",
       ".xr-array-preview,\n",
       ".xr-array-data {\n",
       "  padding: 0 5px !important;\n",
       "  grid-column: 2;\n",
       "}\n",
       "\n",
       ".xr-array-data,\n",
       ".xr-array-in:checked ~ .xr-array-preview {\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-array-in:checked ~ .xr-array-data,\n",
       ".xr-array-preview {\n",
       "  display: inline-block;\n",
       "}\n",
       "\n",
       ".xr-dim-list {\n",
       "  display: inline-block !important;\n",
       "  list-style: none;\n",
       "  padding: 0 !important;\n",
       "  margin: 0;\n",
       "}\n",
       "\n",
       ".xr-dim-list li {\n",
       "  display: inline-block;\n",
       "  padding: 0;\n",
       "  margin: 0;\n",
       "}\n",
       "\n",
       ".xr-dim-list:before {\n",
       "  content: '(';\n",
       "}\n",
       "\n",
       ".xr-dim-list:after {\n",
       "  content: ')';\n",
       "}\n",
       "\n",
       ".xr-dim-list li:not(:last-child):after {\n",
       "  content: ',';\n",
       "  padding-right: 5px;\n",
       "}\n",
       "\n",
       ".xr-has-index {\n",
       "  font-weight: bold;\n",
       "}\n",
       "\n",
       ".xr-var-list,\n",
       ".xr-var-item {\n",
       "  display: contents;\n",
       "}\n",
       "\n",
       ".xr-var-item > div,\n",
       ".xr-var-item label,\n",
       ".xr-var-item > .xr-var-name span {\n",
       "  background-color: var(--xr-background-color-row-even);\n",
       "  margin-bottom: 0;\n",
       "}\n",
       "\n",
       ".xr-var-item > .xr-var-name:hover span {\n",
       "  padding-right: 5px;\n",
       "}\n",
       "\n",
       ".xr-var-list > li:nth-child(odd) > div,\n",
       ".xr-var-list > li:nth-child(odd) > label,\n",
       ".xr-var-list > li:nth-child(odd) > .xr-var-name span {\n",
       "  background-color: var(--xr-background-color-row-odd);\n",
       "}\n",
       "\n",
       ".xr-var-name {\n",
       "  grid-column: 1;\n",
       "}\n",
       "\n",
       ".xr-var-dims {\n",
       "  grid-column: 2;\n",
       "}\n",
       "\n",
       ".xr-var-dtype {\n",
       "  grid-column: 3;\n",
       "  text-align: right;\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-var-preview {\n",
       "  grid-column: 4;\n",
       "}\n",
       "\n",
       ".xr-var-name,\n",
       ".xr-var-dims,\n",
       ".xr-var-dtype,\n",
       ".xr-preview,\n",
       ".xr-attrs dt {\n",
       "  white-space: nowrap;\n",
       "  overflow: hidden;\n",
       "  text-overflow: ellipsis;\n",
       "  padding-right: 10px;\n",
       "}\n",
       "\n",
       ".xr-var-name:hover,\n",
       ".xr-var-dims:hover,\n",
       ".xr-var-dtype:hover,\n",
       ".xr-attrs dt:hover {\n",
       "  overflow: visible;\n",
       "  width: auto;\n",
       "  z-index: 1;\n",
       "}\n",
       "\n",
       ".xr-var-attrs,\n",
       ".xr-var-data {\n",
       "  display: none;\n",
       "  background-color: var(--xr-background-color) !important;\n",
       "  padding-bottom: 5px !important;\n",
       "}\n",
       "\n",
       ".xr-var-attrs-in:checked ~ .xr-var-attrs,\n",
       ".xr-var-data-in:checked ~ .xr-var-data {\n",
       "  display: block;\n",
       "}\n",
       "\n",
       ".xr-var-data > table {\n",
       "  float: right;\n",
       "}\n",
       "\n",
       ".xr-var-name span,\n",
       ".xr-var-data,\n",
       ".xr-attrs {\n",
       "  padding-left: 25px !important;\n",
       "}\n",
       "\n",
       ".xr-attrs,\n",
       ".xr-var-attrs,\n",
       ".xr-var-data {\n",
       "  grid-column: 1 / -1;\n",
       "}\n",
       "\n",
       "dl.xr-attrs {\n",
       "  padding: 0;\n",
       "  margin: 0;\n",
       "  display: grid;\n",
       "  grid-template-columns: 125px auto;\n",
       "}\n",
       "\n",
       ".xr-attrs dt,\n",
       ".xr-attrs dd {\n",
       "  padding: 0;\n",
       "  margin: 0;\n",
       "  float: left;\n",
       "  padding-right: 10px;\n",
       "  width: auto;\n",
       "}\n",
       "\n",
       ".xr-attrs dt {\n",
       "  font-weight: normal;\n",
       "  grid-column: 1;\n",
       "}\n",
       "\n",
       ".xr-attrs dt:hover span {\n",
       "  display: inline-block;\n",
       "  background: var(--xr-background-color);\n",
       "  padding-right: 10px;\n",
       "}\n",
       "\n",
       ".xr-attrs dd {\n",
       "  grid-column: 2;\n",
       "  white-space: pre-wrap;\n",
       "  word-break: break-all;\n",
       "}\n",
       "\n",
       ".xr-icon-database,\n",
       ".xr-icon-file-text2 {\n",
       "  display: inline-block;\n",
       "  vertical-align: middle;\n",
       "  width: 1em;\n",
       "  height: 1.5em !important;\n",
       "  stroke-width: 0;\n",
       "  stroke: currentColor;\n",
       "  fill: currentColor;\n",
       "}\n",
       "</style><pre class='xr-text-repr-fallback'>&lt;xarray.Dataset&gt;\n",
       "Dimensions:      (lon: 20, lat: 19, percentile: 5, constrained: 2, season: 2)\n",
       "Coordinates:\n",
       "  * lon          (lon) float64 -8.75 -6.25 -3.75 -1.25 ... 33.75 36.25 38.75\n",
       "  * lat          (lat) float64 31.25 33.75 36.25 38.75 ... 71.25 73.75 76.25\n",
       "  * percentile   (percentile) int64 10 25 50 75 90\n",
       "  * constrained  (constrained) int64 1 0\n",
       "  * season       (season) &lt;U3 &#x27;DJF&#x27; &#x27;JJA&#x27;\n",
       "Data variables:\n",
       "    tas          (season, constrained, percentile, lat, lon) float64 dask.array&lt;chunksize=(1, 1, 1, 19, 20), meta=np.ndarray&gt;\n",
       "    pr           (season, constrained, percentile, lat, lon) float64 dask.array&lt;chunksize=(1, 1, 1, 19, 20), meta=np.ndarray&gt;</pre><div class='xr-wrap' hidden><div class='xr-header'><div class='xr-obj-type'>xarray.Dataset</div></div><ul class='xr-sections'><li class='xr-section-item'><input id='section-be483d0b-369d-4337-9c24-3b2810524a2e' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-be483d0b-369d-4337-9c24-3b2810524a2e' class='xr-section-summary'  title='Expand/collapse section'>Dimensions:</label><div class='xr-section-inline-details'><ul class='xr-dim-list'><li><span class='xr-has-index'>lon</span>: 20</li><li><span class='xr-has-index'>lat</span>: 19</li><li><span class='xr-has-index'>percentile</span>: 5</li><li><span class='xr-has-index'>constrained</span>: 2</li><li><span class='xr-has-index'>season</span>: 2</li></ul></div><div class='xr-section-details'></div></li><li class='xr-section-item'><input id='section-b7f229b3-85fe-4f12-be39-4123838ff0d6' class='xr-section-summary-in' type='checkbox'  checked><label for='section-b7f229b3-85fe-4f12-be39-4123838ff0d6' class='xr-section-summary' >Coordinates: <span>(5)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><ul class='xr-var-list'><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>lon</span></div><div class='xr-var-dims'>(lon)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>-8.75 -6.25 -3.75 ... 36.25 38.75</div><input id='attrs-400d1eee-7c44-442c-9b29-e5c153b952e1' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-400d1eee-7c44-442c-9b29-e5c153b952e1' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-b6ceb29a-aa28-48e2-9ff4-d671c279b5f8' class='xr-var-data-in' type='checkbox'><label for='data-b6ceb29a-aa28-48e2-9ff4-d671c279b5f8' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>standard_name :</span></dt><dd>longitude</dd><dt><span>long_name :</span></dt><dd>longitude</dd><dt><span>units :</span></dt><dd>degrees_east</dd><dt><span>axis :</span></dt><dd>X</dd></dl></div><div class='xr-var-data'><pre>array([-8.75, -6.25, -3.75, -1.25,  1.25,  3.75,  6.25,  8.75, 11.25, 13.75,\n",
       "       16.25, 18.75, 21.25, 23.75, 26.25, 28.75, 31.25, 33.75, 36.25, 38.75])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>lat</span></div><div class='xr-var-dims'>(lat)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>31.25 33.75 36.25 ... 73.75 76.25</div><input id='attrs-db96d92d-a990-4016-9db3-ed6c232340e8' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-db96d92d-a990-4016-9db3-ed6c232340e8' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-ca209575-ac06-429f-90c0-50626313875d' class='xr-var-data-in' type='checkbox'><label for='data-ca209575-ac06-429f-90c0-50626313875d' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>standard_name :</span></dt><dd>latitude</dd><dt><span>long_name :</span></dt><dd>latitude</dd><dt><span>units :</span></dt><dd>degrees_north</dd><dt><span>axis :</span></dt><dd>Y</dd></dl></div><div class='xr-var-data'><pre>array([31.25, 33.75, 36.25, 38.75, 41.25, 43.75, 46.25, 48.75, 51.25, 53.75,\n",
       "       56.25, 58.75, 61.25, 63.75, 66.25, 68.75, 71.25, 73.75, 76.25])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>percentile</span></div><div class='xr-var-dims'>(percentile)</div><div class='xr-var-dtype'>int64</div><div class='xr-var-preview xr-preview'>10 25 50 75 90</div><input id='attrs-0a2daf2d-10e8-4b56-b91d-b1c20379bce6' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-0a2daf2d-10e8-4b56-b91d-b1c20379bce6' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-1124445f-63d8-4dc5-8f93-3efc85460ae0' class='xr-var-data-in' type='checkbox'><label for='data-1124445f-63d8-4dc5-8f93-3efc85460ae0' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([10, 25, 50, 75, 90])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>constrained</span></div><div class='xr-var-dims'>(constrained)</div><div class='xr-var-dtype'>int64</div><div class='xr-var-preview xr-preview'>1 0</div><input id='attrs-0db1ef0d-9c5c-4aeb-a8d5-8a145c3cbab4' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-0db1ef0d-9c5c-4aeb-a8d5-8a145c3cbab4' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-b9d932a5-398b-4ecc-bbe4-87ad619eb95d' class='xr-var-data-in' type='checkbox'><label for='data-b9d932a5-398b-4ecc-bbe4-87ad619eb95d' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([1, 0])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>season</span></div><div class='xr-var-dims'>(season)</div><div class='xr-var-dtype'>&lt;U3</div><div class='xr-var-preview xr-preview'>&#x27;DJF&#x27; &#x27;JJA&#x27;</div><input id='attrs-a3a9f38d-0fad-41b3-8128-447e473d9dc6' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-a3a9f38d-0fad-41b3-8128-447e473d9dc6' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-3111a16e-5e9e-458b-9369-8cd3079c6b72' class='xr-var-data-in' type='checkbox'><label for='data-3111a16e-5e9e-458b-9369-8cd3079c6b72' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([&#x27;DJF&#x27;, &#x27;JJA&#x27;], dtype=&#x27;&lt;U3&#x27;)</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-15288e02-0f40-4c8c-9999-d3268557eb48' class='xr-section-summary-in' type='checkbox'  checked><label for='section-15288e02-0f40-4c8c-9999-d3268557eb48' class='xr-section-summary' >Data variables: <span>(2)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><ul class='xr-var-list'><li class='xr-var-item'><div class='xr-var-name'><span>tas</span></div><div class='xr-var-dims'>(season, constrained, percentile, lat, lon)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>dask.array&lt;chunksize=(1, 1, 1, 19, 20), meta=np.ndarray&gt;</div><input id='attrs-f1959ad9-8b42-4c0e-a255-99f3b973f742' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-f1959ad9-8b42-4c0e-a255-99f3b973f742' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-1366e00f-a811-48dd-bd5a-f8e55af5ebe9' class='xr-var-data-in' type='checkbox'><label for='data-1366e00f-a811-48dd-bd5a-f8e55af5ebe9' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><table>\n",
       "    <tr>\n",
       "        <td>\n",
       "            <table>\n",
       "                <thead>\n",
       "                    <tr>\n",
       "                        <td> </td>\n",
       "                        <th> Array </th>\n",
       "                        <th> Chunk </th>\n",
       "                    </tr>\n",
       "                </thead>\n",
       "                <tbody>\n",
       "                    \n",
       "                    <tr>\n",
       "                        <th> Bytes </th>\n",
       "                        <td> 59.38 kiB </td>\n",
       "                        <td> 2.97 kiB </td>\n",
       "                    </tr>\n",
       "                    \n",
       "                    <tr>\n",
       "                        <th> Shape </th>\n",
       "                        <td> (2, 2, 5, 19, 20) </td>\n",
       "                        <td> (1, 1, 1, 19, 20) </td>\n",
       "                    </tr>\n",
       "                    <tr>\n",
       "                        <th> Count </th>\n",
       "                        <td> 160 Tasks </td>\n",
       "                        <td> 20 Chunks </td>\n",
       "                    </tr>\n",
       "                    <tr>\n",
       "                    <th> Type </th>\n",
       "                    <td> float64 </td>\n",
       "                    <td> numpy.ndarray </td>\n",
       "                    </tr>\n",
       "                </tbody>\n",
       "            </table>\n",
       "        </td>\n",
       "        <td>\n",
       "        <svg width=\"411\" height=\"189\" style=\"stroke:rgb(0,0,0);stroke-width:1\" >\n",
       "\n",
       "  <!-- Horizontal lines -->\n",
       "  <line x1=\"0\" y1=\"0\" x2=\"38\" y2=\"0\" style=\"stroke-width:2\" />\n",
       "  <line x1=\"0\" y1=\"19\" x2=\"38\" y2=\"19\" />\n",
       "  <line x1=\"0\" y1=\"38\" x2=\"38\" y2=\"38\" style=\"stroke-width:2\" />\n",
       "\n",
       "  <!-- Vertical lines -->\n",
       "  <line x1=\"0\" y1=\"0\" x2=\"0\" y2=\"38\" style=\"stroke-width:2\" />\n",
       "  <line x1=\"19\" y1=\"0\" x2=\"19\" y2=\"38\" />\n",
       "  <line x1=\"38\" y1=\"0\" x2=\"38\" y2=\"38\" style=\"stroke-width:2\" />\n",
       "\n",
       "  <!-- Colored Rectangle -->\n",
       "  <polygon points=\"0.0,0.0 38.596863036086,0.0 38.596863036086,38.596863036086 0.0,38.596863036086\" style=\"fill:#ECB172A0;stroke-width:0\"/>\n",
       "\n",
       "  <!-- Text -->\n",
       "  <text x=\"19.298432\" y=\"58.596863\" font-size=\"1.0rem\" font-weight=\"100\" text-anchor=\"middle\" >2</text>\n",
       "  <text x=\"58.596863\" y=\"19.298432\" font-size=\"1.0rem\" font-weight=\"100\" text-anchor=\"middle\" transform=\"rotate(0,58.596863,19.298432)\">2</text>\n",
       "\n",
       "\n",
       "  <!-- Horizontal lines -->\n",
       "  <line x1=\"108\" y1=\"0\" x2=\"133\" y2=\"25\" style=\"stroke-width:2\" />\n",
       "  <line x1=\"108\" y1=\"114\" x2=\"133\" y2=\"139\" style=\"stroke-width:2\" />\n",
       "\n",
       "  <!-- Vertical lines -->\n",
       "  <line x1=\"108\" y1=\"0\" x2=\"108\" y2=\"114\" style=\"stroke-width:2\" />\n",
       "  <line x1=\"113\" y1=\"5\" x2=\"113\" y2=\"119\" />\n",
       "  <line x1=\"118\" y1=\"10\" x2=\"118\" y2=\"124\" />\n",
       "  <line x1=\"123\" y1=\"15\" x2=\"123\" y2=\"129\" />\n",
       "  <line x1=\"128\" y1=\"20\" x2=\"128\" y2=\"134\" />\n",
       "  <line x1=\"133\" y1=\"25\" x2=\"133\" y2=\"139\" style=\"stroke-width:2\" />\n",
       "\n",
       "  <!-- Colored Rectangle -->\n",
       "  <polygon points=\"108.0,0.0 133.23458560109532,25.23458560109531 133.23458560109532,139.23458560109532 108.0,114.0\" style=\"fill:#ECB172A0;stroke-width:0\"/>\n",
       "\n",
       "  <!-- Horizontal lines -->\n",
       "  <line x1=\"108\" y1=\"0\" x2=\"228\" y2=\"0\" style=\"stroke-width:2\" />\n",
       "  <line x1=\"113\" y1=\"5\" x2=\"233\" y2=\"5\" />\n",
       "  <line x1=\"118\" y1=\"10\" x2=\"238\" y2=\"10\" />\n",
       "  <line x1=\"123\" y1=\"15\" x2=\"243\" y2=\"15\" />\n",
       "  <line x1=\"128\" y1=\"20\" x2=\"248\" y2=\"20\" />\n",
       "  <line x1=\"133\" y1=\"25\" x2=\"253\" y2=\"25\" style=\"stroke-width:2\" />\n",
       "\n",
       "  <!-- Vertical lines -->\n",
       "  <line x1=\"108\" y1=\"0\" x2=\"133\" y2=\"25\" style=\"stroke-width:2\" />\n",
       "  <line x1=\"228\" y1=\"0\" x2=\"253\" y2=\"25\" style=\"stroke-width:2\" />\n",
       "\n",
       "  <!-- Colored Rectangle -->\n",
       "  <polygon points=\"108.0,0.0 228.0,0.0 253.23458560109532,25.23458560109531 133.23458560109532,25.23458560109531\" style=\"fill:#ECB172A0;stroke-width:0\"/>\n",
       "\n",
       "  <!-- Horizontal lines -->\n",
       "  <line x1=\"133\" y1=\"25\" x2=\"253\" y2=\"25\" style=\"stroke-width:2\" />\n",
       "  <line x1=\"133\" y1=\"139\" x2=\"253\" y2=\"139\" style=\"stroke-width:2\" />\n",
       "\n",
       "  <!-- Vertical lines -->\n",
       "  <line x1=\"133\" y1=\"25\" x2=\"133\" y2=\"139\" style=\"stroke-width:2\" />\n",
       "  <line x1=\"253\" y1=\"25\" x2=\"253\" y2=\"139\" style=\"stroke-width:2\" />\n",
       "\n",
       "  <!-- Colored Rectangle -->\n",
       "  <polygon points=\"133.23458560109532,25.23458560109531 253.23458560109532,25.23458560109531 253.23458560109532,139.23458560109532 133.23458560109532,139.23458560109532\" style=\"fill:#ECB172A0;stroke-width:0\"/>\n",
       "\n",
       "  <!-- Text -->\n",
       "  <text x=\"193.234586\" y=\"159.234586\" font-size=\"1.0rem\" font-weight=\"100\" text-anchor=\"middle\" >20</text>\n",
       "  <text x=\"273.234586\" y=\"82.234586\" font-size=\"1.0rem\" font-weight=\"100\" text-anchor=\"middle\" transform=\"rotate(0,273.234586,82.234586)\">19</text>\n",
       "  <text x=\"110.617293\" y=\"146.617293\" font-size=\"1.0rem\" font-weight=\"100\" text-anchor=\"middle\" transform=\"rotate(45,110.617293,146.617293)\">5</text>\n",
       "</svg>\n",
       "        </td>\n",
       "    </tr>\n",
       "</table></div></li><li class='xr-var-item'><div class='xr-var-name'><span>pr</span></div><div class='xr-var-dims'>(season, constrained, percentile, lat, lon)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>dask.array&lt;chunksize=(1, 1, 1, 19, 20), meta=np.ndarray&gt;</div><input id='attrs-e3f6052c-be11-45fe-961e-fd844c333f66' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-e3f6052c-be11-45fe-961e-fd844c333f66' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-5ef34a42-70ce-4687-a0ee-427e9bdaa614' class='xr-var-data-in' type='checkbox'><label for='data-5ef34a42-70ce-4687-a0ee-427e9bdaa614' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><table>\n",
       "    <tr>\n",
       "        <td>\n",
       "            <table>\n",
       "                <thead>\n",
       "                    <tr>\n",
       "                        <td> </td>\n",
       "                        <th> Array </th>\n",
       "                        <th> Chunk </th>\n",
       "                    </tr>\n",
       "                </thead>\n",
       "                <tbody>\n",
       "                    \n",
       "                    <tr>\n",
       "                        <th> Bytes </th>\n",
       "                        <td> 59.38 kiB </td>\n",
       "                        <td> 2.97 kiB </td>\n",
       "                    </tr>\n",
       "                    \n",
       "                    <tr>\n",
       "                        <th> Shape </th>\n",
       "                        <td> (2, 2, 5, 19, 20) </td>\n",
       "                        <td> (1, 1, 1, 19, 20) </td>\n",
       "                    </tr>\n",
       "                    <tr>\n",
       "                        <th> Count </th>\n",
       "                        <td> 160 Tasks </td>\n",
       "                        <td> 20 Chunks </td>\n",
       "                    </tr>\n",
       "                    <tr>\n",
       "                    <th> Type </th>\n",
       "                    <td> float64 </td>\n",
       "                    <td> numpy.ndarray </td>\n",
       "                    </tr>\n",
       "                </tbody>\n",
       "            </table>\n",
       "        </td>\n",
       "        <td>\n",
       "        <svg width=\"411\" height=\"189\" style=\"stroke:rgb(0,0,0);stroke-width:1\" >\n",
       "\n",
       "  <!-- Horizontal lines -->\n",
       "  <line x1=\"0\" y1=\"0\" x2=\"38\" y2=\"0\" style=\"stroke-width:2\" />\n",
       "  <line x1=\"0\" y1=\"19\" x2=\"38\" y2=\"19\" />\n",
       "  <line x1=\"0\" y1=\"38\" x2=\"38\" y2=\"38\" style=\"stroke-width:2\" />\n",
       "\n",
       "  <!-- Vertical lines -->\n",
       "  <line x1=\"0\" y1=\"0\" x2=\"0\" y2=\"38\" style=\"stroke-width:2\" />\n",
       "  <line x1=\"19\" y1=\"0\" x2=\"19\" y2=\"38\" />\n",
       "  <line x1=\"38\" y1=\"0\" x2=\"38\" y2=\"38\" style=\"stroke-width:2\" />\n",
       "\n",
       "  <!-- Colored Rectangle -->\n",
       "  <polygon points=\"0.0,0.0 38.596863036086,0.0 38.596863036086,38.596863036086 0.0,38.596863036086\" style=\"fill:#ECB172A0;stroke-width:0\"/>\n",
       "\n",
       "  <!-- Text -->\n",
       "  <text x=\"19.298432\" y=\"58.596863\" font-size=\"1.0rem\" font-weight=\"100\" text-anchor=\"middle\" >2</text>\n",
       "  <text x=\"58.596863\" y=\"19.298432\" font-size=\"1.0rem\" font-weight=\"100\" text-anchor=\"middle\" transform=\"rotate(0,58.596863,19.298432)\">2</text>\n",
       "\n",
       "\n",
       "  <!-- Horizontal lines -->\n",
       "  <line x1=\"108\" y1=\"0\" x2=\"133\" y2=\"25\" style=\"stroke-width:2\" />\n",
       "  <line x1=\"108\" y1=\"114\" x2=\"133\" y2=\"139\" style=\"stroke-width:2\" />\n",
       "\n",
       "  <!-- Vertical lines -->\n",
       "  <line x1=\"108\" y1=\"0\" x2=\"108\" y2=\"114\" style=\"stroke-width:2\" />\n",
       "  <line x1=\"113\" y1=\"5\" x2=\"113\" y2=\"119\" />\n",
       "  <line x1=\"118\" y1=\"10\" x2=\"118\" y2=\"124\" />\n",
       "  <line x1=\"123\" y1=\"15\" x2=\"123\" y2=\"129\" />\n",
       "  <line x1=\"128\" y1=\"20\" x2=\"128\" y2=\"134\" />\n",
       "  <line x1=\"133\" y1=\"25\" x2=\"133\" y2=\"139\" style=\"stroke-width:2\" />\n",
       "\n",
       "  <!-- Colored Rectangle -->\n",
       "  <polygon points=\"108.0,0.0 133.23458560109532,25.23458560109531 133.23458560109532,139.23458560109532 108.0,114.0\" style=\"fill:#ECB172A0;stroke-width:0\"/>\n",
       "\n",
       "  <!-- Horizontal lines -->\n",
       "  <line x1=\"108\" y1=\"0\" x2=\"228\" y2=\"0\" style=\"stroke-width:2\" />\n",
       "  <line x1=\"113\" y1=\"5\" x2=\"233\" y2=\"5\" />\n",
       "  <line x1=\"118\" y1=\"10\" x2=\"238\" y2=\"10\" />\n",
       "  <line x1=\"123\" y1=\"15\" x2=\"243\" y2=\"15\" />\n",
       "  <line x1=\"128\" y1=\"20\" x2=\"248\" y2=\"20\" />\n",
       "  <line x1=\"133\" y1=\"25\" x2=\"253\" y2=\"25\" style=\"stroke-width:2\" />\n",
       "\n",
       "  <!-- Vertical lines -->\n",
       "  <line x1=\"108\" y1=\"0\" x2=\"133\" y2=\"25\" style=\"stroke-width:2\" />\n",
       "  <line x1=\"228\" y1=\"0\" x2=\"253\" y2=\"25\" style=\"stroke-width:2\" />\n",
       "\n",
       "  <!-- Colored Rectangle -->\n",
       "  <polygon points=\"108.0,0.0 228.0,0.0 253.23458560109532,25.23458560109531 133.23458560109532,25.23458560109531\" style=\"fill:#ECB172A0;stroke-width:0\"/>\n",
       "\n",
       "  <!-- Horizontal lines -->\n",
       "  <line x1=\"133\" y1=\"25\" x2=\"253\" y2=\"25\" style=\"stroke-width:2\" />\n",
       "  <line x1=\"133\" y1=\"139\" x2=\"253\" y2=\"139\" style=\"stroke-width:2\" />\n",
       "\n",
       "  <!-- Vertical lines -->\n",
       "  <line x1=\"133\" y1=\"25\" x2=\"133\" y2=\"139\" style=\"stroke-width:2\" />\n",
       "  <line x1=\"253\" y1=\"25\" x2=\"253\" y2=\"139\" style=\"stroke-width:2\" />\n",
       "\n",
       "  <!-- Colored Rectangle -->\n",
       "  <polygon points=\"133.23458560109532,25.23458560109531 253.23458560109532,25.23458560109531 253.23458560109532,139.23458560109532 133.23458560109532,139.23458560109532\" style=\"fill:#ECB172A0;stroke-width:0\"/>\n",
       "\n",
       "  <!-- Text -->\n",
       "  <text x=\"193.234586\" y=\"159.234586\" font-size=\"1.0rem\" font-weight=\"100\" text-anchor=\"middle\" >20</text>\n",
       "  <text x=\"273.234586\" y=\"82.234586\" font-size=\"1.0rem\" font-weight=\"100\" text-anchor=\"middle\" transform=\"rotate(0,273.234586,82.234586)\">19</text>\n",
       "  <text x=\"110.617293\" y=\"146.617293\" font-size=\"1.0rem\" font-weight=\"100\" text-anchor=\"middle\" transform=\"rotate(45,110.617293,146.617293)\">5</text>\n",
       "</svg>\n",
       "        </td>\n",
       "    </tr>\n",
       "</table></div></li></ul></div></li><li class='xr-section-item'><input id='section-9e177d91-1613-4409-8d55-4309f9ab8bff' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-9e177d91-1613-4409-8d55-4309f9ab8bff' class='xr-section-summary'  title='Expand/collapse section'>Attributes: <span>(0)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><dl class='xr-attrs'></dl></div></li></ul></div></div>"
      ],
      "text/plain": [
       "<xarray.Dataset>\n",
       "Dimensions:      (lon: 20, lat: 19, percentile: 5, constrained: 2, season: 2)\n",
       "Coordinates:\n",
       "  * lon          (lon) float64 -8.75 -6.25 -3.75 -1.25 ... 33.75 36.25 38.75\n",
       "  * lat          (lat) float64 31.25 33.75 36.25 38.75 ... 71.25 73.75 76.25\n",
       "  * percentile   (percentile) int64 10 25 50 75 90\n",
       "  * constrained  (constrained) int64 1 0\n",
       "  * season       (season) <U3 'DJF' 'JJA'\n",
       "Data variables:\n",
       "    tas          (season, constrained, percentile, lat, lon) float64 dask.array<chunksize=(1, 1, 1, 19, 20), meta=np.ndarray>\n",
       "    pr           (season, constrained, percentile, lat, lon) float64 dask.array<chunksize=(1, 1, 1, 19, 20), meta=np.ndarray>"
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "project = output_file_name[\"source\"].lower()\n",
    "seasons = []\n",
    "for season in [\"djf\", \"jja\"]:\n",
    "    tas = load_data(project, season, \"tas\")\n",
    "    pr = load_data(project, season, \"pr\")\n",
    "    ds = xr.merge([tas, pr]).assign_coords(season=season.upper())\n",
    "    seasons.append(ds)\n",
    "ethz_climwip_ds = xr.concat(seasons, dim=\"season\")\n",
    "ethz_climwip_ds"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Make some metadata. Here, we follow CF-conventions as much as possible."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [],
   "source": [
    "attrs = {\n",
    "    \"tas\": {\n",
    "        \"description\": \"Change in Air Temperature\",\n",
    "        \"standard_name\": \"Change in Air Temperature\",\n",
    "        \"long_name\": \"Change in Near-Surface Air Temperature\",\n",
    "        \"units\": \"K\",  # in line with raw data\n",
    "        \"cell_methods\": \"time: mean changes over 20 years 2041-2060 vs 1995-2014\",\n",
    "    },\n",
    "    \"pr\": {\n",
    "        \"description\": \"Relative precipitation\",\n",
    "        \"standard_name\": \"Relative precipitation\",\n",
    "        \"long_name\": \"Relative precipitation\",\n",
    "        \"units\": \"%\",  # in line with raw data\n",
    "        \"cell_methods\": \"time: mean changes over 20 years 2041-2060 vs 1995-2014\",\n",
    "    },\n",
    "    \"latitude\": {\"units\": \"degrees_north\", \"long_name\": \"latitude\", \"axis\": \"Y\"},\n",
    "    \"longitude\": {\"units\": \"degrees_east\", \"long_name\": \"longitude\", \"axis\": \"X\"},\n",
    "    \"time\": {\n",
    "        \"climatology\": \"climatology_bounds\",\n",
    "        \"long_name\": \"time\",\n",
    "        \"axis\": \"T\",\n",
    "        \"climatology_bounds\": [\"2050-6-1\", \"2050-9-1\", \"2050-12-1\", \"2051-3-1\"],\n",
    "        \"description\": \"mean changes over 20 years 2041-2060 vs 1995-2014. The mid point 2050 is chosen as the representative time.\",\n",
    "    },\n",
    "    \"percentile\": {\"units\": \"%\", \"long_name\": \"percentile\", \"axis\": \"Z\"},\n",
    "}"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Assemble data and save to netcdf"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Make a function to assemble and save data "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [],
   "source": [
    "TIMES = {\n",
    "    \"JJA\": \"2050-7-16\",\n",
    "    \"DJF\": \"2051-1-16\",\n",
    "}  # \"0000-4-16\", \"0000-7-16\", \"0000-10-16\", \"0000-1-16\" MAM JJA SON DJF\n",
    "PERCENTILES = [10, 25, 50, 75, 90]\n",
    "\n",
    "LAT = ethz_climwip_ds.coords[\"lat\"]\n",
    "LON = ethz_climwip_ds.coords[\"lon\"]\n",
    "\n",
    "\n",
    "def assembly(ds_original, var, cons):\n",
    "    \"\"\"\n",
    "    Select data from original nc files and save the target fields.\n",
    "    \"\"\"\n",
    "    ds_target = xr.Dataset(\n",
    "        {\n",
    "            var: (\n",
    "                (\"time\", \"latitude\", \"longitude\", \"percentile\"),\n",
    "                np.full([len(TIMES), len(LAT), len(LON), len(PERCENTILES)], np.nan),\n",
    "            ),\n",
    "            \"climatology_bounds\": (\n",
    "                pd.to_datetime([\"2050-6-1\", \"2050-9-1\", \"2050-12-1\", \"2051-3-1\"])\n",
    "            ),\n",
    "        },\n",
    "        coords={\n",
    "            \"time\": pd.to_datetime(list(TIMES.values())),\n",
    "            \"latitude\": LAT.values,\n",
    "            \"longitude\": LON.values,\n",
    "            \"percentile\": PERCENTILES,\n",
    "        },\n",
    "        attrs={\n",
    "            \"description\": f\"Contains modified {INSTITUTION_ID} {METHOD} data used for Atlas in EUCP project.\",\n",
    "            \"history\": f\"original {INSTITUTION_ID} {METHOD} data files eur_pr_41-60_djf_cmip6_*perc_rel.nc, eur_tas_41-60_djf_cmip6_*perc.nc\",\n",
    "        },\n",
    "    )\n",
    "    for season in [\"JJA\", \"DJF\"]:\n",
    "        for j, p in enumerate(PERCENTILES):\n",
    "            ds_target[var].values[list(TIMES).index(season), :, :, j] = (\n",
    "                ds_original[var]\n",
    "                .sel(percentile=p, season=season, constrained=cons)\n",
    "                .values\n",
    "            )\n",
    "    return ds_target"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Call the function"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "one dataset is saved to atlas_EUCP_ETHZ_CMIP6_ClimWIP_uncons_tas.nc\n",
      "one dataset is saved to atlas_EUCP_ETHZ_CMIP6_ClimWIP_cons_tas.nc\n",
      "one dataset is saved to atlas_EUCP_ETHZ_CMIP6_ClimWIP_uncons_pr.nc\n",
      "one dataset is saved to atlas_EUCP_ETHZ_CMIP6_ClimWIP_cons_pr.nc\n"
     ]
    }
   ],
   "source": [
    "for VAR_NAME in [\"tas\", \"pr\"]:\n",
    "    output_file_name[\"cmor_var\"] = VAR_NAME\n",
    "    for i, sub_method in enumerate([\"uncons\", \"cons\"]):\n",
    "        output_file_name[\"sub_method\"] = sub_method\n",
    "        new_ds = assembly(ethz_climwip_ds, VAR_NAME, i)\n",
    "\n",
    "        # Fix attributes\n",
    "        for key in new_ds.keys():\n",
    "            new_ds[key].attrs = attrs[key]\n",
    "\n",
    "        file_name = f\"{'_'.join(output_file_name.values())}.nc\"\n",
    "        print(f\"one dataset is saved to {file_name}\")\n",
    "        new_ds.to_netcdf(output_path / file_name)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Check output"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Load one of the saved data."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div><svg style=\"position: absolute; width: 0; height: 0; overflow: hidden\">\n",
       "<defs>\n",
       "<symbol id=\"icon-database\" viewBox=\"0 0 32 32\">\n",
       "<path d=\"M16 0c-8.837 0-16 2.239-16 5v4c0 2.761 7.163 5 16 5s16-2.239 16-5v-4c0-2.761-7.163-5-16-5z\"></path>\n",
       "<path d=\"M16 17c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n",
       "<path d=\"M16 26c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n",
       "</symbol>\n",
       "<symbol id=\"icon-file-text2\" viewBox=\"0 0 32 32\">\n",
       "<path d=\"M28.681 7.159c-0.694-0.947-1.662-2.053-2.724-3.116s-2.169-2.030-3.116-2.724c-1.612-1.182-2.393-1.319-2.841-1.319h-15.5c-1.378 0-2.5 1.121-2.5 2.5v27c0 1.378 1.122 2.5 2.5 2.5h23c1.378 0 2.5-1.122 2.5-2.5v-19.5c0-0.448-0.137-1.23-1.319-2.841zM24.543 5.457c0.959 0.959 1.712 1.825 2.268 2.543h-4.811v-4.811c0.718 0.556 1.584 1.309 2.543 2.268zM28 29.5c0 0.271-0.229 0.5-0.5 0.5h-23c-0.271 0-0.5-0.229-0.5-0.5v-27c0-0.271 0.229-0.5 0.5-0.5 0 0 15.499-0 15.5 0v7c0 0.552 0.448 1 1 1h7v19.5z\"></path>\n",
       "<path d=\"M23 26h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "<path d=\"M23 22h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "<path d=\"M23 18h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "</symbol>\n",
       "</defs>\n",
       "</svg>\n",
       "<style>/* CSS stylesheet for displaying xarray objects in jupyterlab.\n",
       " *\n",
       " */\n",
       "\n",
       ":root {\n",
       "  --xr-font-color0: var(--jp-content-font-color0, rgba(0, 0, 0, 1));\n",
       "  --xr-font-color2: var(--jp-content-font-color2, rgba(0, 0, 0, 0.54));\n",
       "  --xr-font-color3: var(--jp-content-font-color3, rgba(0, 0, 0, 0.38));\n",
       "  --xr-border-color: var(--jp-border-color2, #e0e0e0);\n",
       "  --xr-disabled-color: var(--jp-layout-color3, #bdbdbd);\n",
       "  --xr-background-color: var(--jp-layout-color0, white);\n",
       "  --xr-background-color-row-even: var(--jp-layout-color1, white);\n",
       "  --xr-background-color-row-odd: var(--jp-layout-color2, #eeeeee);\n",
       "}\n",
       "\n",
       "html[theme=dark],\n",
       "body.vscode-dark {\n",
       "  --xr-font-color0: rgba(255, 255, 255, 1);\n",
       "  --xr-font-color2: rgba(255, 255, 255, 0.54);\n",
       "  --xr-font-color3: rgba(255, 255, 255, 0.38);\n",
       "  --xr-border-color: #1F1F1F;\n",
       "  --xr-disabled-color: #515151;\n",
       "  --xr-background-color: #111111;\n",
       "  --xr-background-color-row-even: #111111;\n",
       "  --xr-background-color-row-odd: #313131;\n",
       "}\n",
       "\n",
       ".xr-wrap {\n",
       "  display: block;\n",
       "  min-width: 300px;\n",
       "  max-width: 700px;\n",
       "}\n",
       "\n",
       ".xr-text-repr-fallback {\n",
       "  /* fallback to plain text repr when CSS is not injected (untrusted notebook) */\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-header {\n",
       "  padding-top: 6px;\n",
       "  padding-bottom: 6px;\n",
       "  margin-bottom: 4px;\n",
       "  border-bottom: solid 1px var(--xr-border-color);\n",
       "}\n",
       "\n",
       ".xr-header > div,\n",
       ".xr-header > ul {\n",
       "  display: inline;\n",
       "  margin-top: 0;\n",
       "  margin-bottom: 0;\n",
       "}\n",
       "\n",
       ".xr-obj-type,\n",
       ".xr-array-name {\n",
       "  margin-left: 2px;\n",
       "  margin-right: 10px;\n",
       "}\n",
       "\n",
       ".xr-obj-type {\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-sections {\n",
       "  padding-left: 0 !important;\n",
       "  display: grid;\n",
       "  grid-template-columns: 150px auto auto 1fr 20px 20px;\n",
       "}\n",
       "\n",
       ".xr-section-item {\n",
       "  display: contents;\n",
       "}\n",
       "\n",
       ".xr-section-item input {\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-section-item input + label {\n",
       "  color: var(--xr-disabled-color);\n",
       "}\n",
       "\n",
       ".xr-section-item input:enabled + label {\n",
       "  cursor: pointer;\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-section-item input:enabled + label:hover {\n",
       "  color: var(--xr-font-color0);\n",
       "}\n",
       "\n",
       ".xr-section-summary {\n",
       "  grid-column: 1;\n",
       "  color: var(--xr-font-color2);\n",
       "  font-weight: 500;\n",
       "}\n",
       "\n",
       ".xr-section-summary > span {\n",
       "  display: inline-block;\n",
       "  padding-left: 0.5em;\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:disabled + label {\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-section-summary-in + label:before {\n",
       "  display: inline-block;\n",
       "  content: '►';\n",
       "  font-size: 11px;\n",
       "  width: 15px;\n",
       "  text-align: center;\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:disabled + label:before {\n",
       "  color: var(--xr-disabled-color);\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:checked + label:before {\n",
       "  content: '▼';\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:checked + label > span {\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-section-summary,\n",
       ".xr-section-inline-details {\n",
       "  padding-top: 4px;\n",
       "  padding-bottom: 4px;\n",
       "}\n",
       "\n",
       ".xr-section-inline-details {\n",
       "  grid-column: 2 / -1;\n",
       "}\n",
       "\n",
       ".xr-section-details {\n",
       "  display: none;\n",
       "  grid-column: 1 / -1;\n",
       "  margin-bottom: 5px;\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:checked ~ .xr-section-details {\n",
       "  display: contents;\n",
       "}\n",
       "\n",
       ".xr-array-wrap {\n",
       "  grid-column: 1 / -1;\n",
       "  display: grid;\n",
       "  grid-template-columns: 20px auto;\n",
       "}\n",
       "\n",
       ".xr-array-wrap > label {\n",
       "  grid-column: 1;\n",
       "  vertical-align: top;\n",
       "}\n",
       "\n",
       ".xr-preview {\n",
       "  color: var(--xr-font-color3);\n",
       "}\n",
       "\n",
       ".xr-array-preview,\n",
       ".xr-array-data {\n",
       "  padding: 0 5px !important;\n",
       "  grid-column: 2;\n",
       "}\n",
       "\n",
       ".xr-array-data,\n",
       ".xr-array-in:checked ~ .xr-array-preview {\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-array-in:checked ~ .xr-array-data,\n",
       ".xr-array-preview {\n",
       "  display: inline-block;\n",
       "}\n",
       "\n",
       ".xr-dim-list {\n",
       "  display: inline-block !important;\n",
       "  list-style: none;\n",
       "  padding: 0 !important;\n",
       "  margin: 0;\n",
       "}\n",
       "\n",
       ".xr-dim-list li {\n",
       "  display: inline-block;\n",
       "  padding: 0;\n",
       "  margin: 0;\n",
       "}\n",
       "\n",
       ".xr-dim-list:before {\n",
       "  content: '(';\n",
       "}\n",
       "\n",
       ".xr-dim-list:after {\n",
       "  content: ')';\n",
       "}\n",
       "\n",
       ".xr-dim-list li:not(:last-child):after {\n",
       "  content: ',';\n",
       "  padding-right: 5px;\n",
       "}\n",
       "\n",
       ".xr-has-index {\n",
       "  font-weight: bold;\n",
       "}\n",
       "\n",
       ".xr-var-list,\n",
       ".xr-var-item {\n",
       "  display: contents;\n",
       "}\n",
       "\n",
       ".xr-var-item > div,\n",
       ".xr-var-item label,\n",
       ".xr-var-item > .xr-var-name span {\n",
       "  background-color: var(--xr-background-color-row-even);\n",
       "  margin-bottom: 0;\n",
       "}\n",
       "\n",
       ".xr-var-item > .xr-var-name:hover span {\n",
       "  padding-right: 5px;\n",
       "}\n",
       "\n",
       ".xr-var-list > li:nth-child(odd) > div,\n",
       ".xr-var-list > li:nth-child(odd) > label,\n",
       ".xr-var-list > li:nth-child(odd) > .xr-var-name span {\n",
       "  background-color: var(--xr-background-color-row-odd);\n",
       "}\n",
       "\n",
       ".xr-var-name {\n",
       "  grid-column: 1;\n",
       "}\n",
       "\n",
       ".xr-var-dims {\n",
       "  grid-column: 2;\n",
       "}\n",
       "\n",
       ".xr-var-dtype {\n",
       "  grid-column: 3;\n",
       "  text-align: right;\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-var-preview {\n",
       "  grid-column: 4;\n",
       "}\n",
       "\n",
       ".xr-var-name,\n",
       ".xr-var-dims,\n",
       ".xr-var-dtype,\n",
       ".xr-preview,\n",
       ".xr-attrs dt {\n",
       "  white-space: nowrap;\n",
       "  overflow: hidden;\n",
       "  text-overflow: ellipsis;\n",
       "  padding-right: 10px;\n",
       "}\n",
       "\n",
       ".xr-var-name:hover,\n",
       ".xr-var-dims:hover,\n",
       ".xr-var-dtype:hover,\n",
       ".xr-attrs dt:hover {\n",
       "  overflow: visible;\n",
       "  width: auto;\n",
       "  z-index: 1;\n",
       "}\n",
       "\n",
       ".xr-var-attrs,\n",
       ".xr-var-data {\n",
       "  display: none;\n",
       "  background-color: var(--xr-background-color) !important;\n",
       "  padding-bottom: 5px !important;\n",
       "}\n",
       "\n",
       ".xr-var-attrs-in:checked ~ .xr-var-attrs,\n",
       ".xr-var-data-in:checked ~ .xr-var-data {\n",
       "  display: block;\n",
       "}\n",
       "\n",
       ".xr-var-data > table {\n",
       "  float: right;\n",
       "}\n",
       "\n",
       ".xr-var-name span,\n",
       ".xr-var-data,\n",
       ".xr-attrs {\n",
       "  padding-left: 25px !important;\n",
       "}\n",
       "\n",
       ".xr-attrs,\n",
       ".xr-var-attrs,\n",
       ".xr-var-data {\n",
       "  grid-column: 1 / -1;\n",
       "}\n",
       "\n",
       "dl.xr-attrs {\n",
       "  padding: 0;\n",
       "  margin: 0;\n",
       "  display: grid;\n",
       "  grid-template-columns: 125px auto;\n",
       "}\n",
       "\n",
       ".xr-attrs dt,\n",
       ".xr-attrs dd {\n",
       "  padding: 0;\n",
       "  margin: 0;\n",
       "  float: left;\n",
       "  padding-right: 10px;\n",
       "  width: auto;\n",
       "}\n",
       "\n",
       ".xr-attrs dt {\n",
       "  font-weight: normal;\n",
       "  grid-column: 1;\n",
       "}\n",
       "\n",
       ".xr-attrs dt:hover span {\n",
       "  display: inline-block;\n",
       "  background: var(--xr-background-color);\n",
       "  padding-right: 10px;\n",
       "}\n",
       "\n",
       ".xr-attrs dd {\n",
       "  grid-column: 2;\n",
       "  white-space: pre-wrap;\n",
       "  word-break: break-all;\n",
       "}\n",
       "\n",
       ".xr-icon-database,\n",
       ".xr-icon-file-text2 {\n",
       "  display: inline-block;\n",
       "  vertical-align: middle;\n",
       "  width: 1em;\n",
       "  height: 1.5em !important;\n",
       "  stroke-width: 0;\n",
       "  stroke: currentColor;\n",
       "  fill: currentColor;\n",
       "}\n",
       "</style><pre class='xr-text-repr-fallback'>&lt;xarray.Dataset&gt;\n",
       "Dimensions:             (time: 2, latitude: 19, longitude: 20, percentile: 5, climatology_bounds: 4)\n",
       "Coordinates:\n",
       "  * climatology_bounds  (climatology_bounds) datetime64[ns] 2050-06-01 ... 20...\n",
       "  * time                (time) datetime64[ns] 2050-07-16 2051-01-16\n",
       "  * latitude            (latitude) float64 31.25 33.75 36.25 ... 73.75 76.25\n",
       "  * longitude           (longitude) float64 -8.75 -6.25 -3.75 ... 36.25 38.75\n",
       "  * percentile          (percentile) int64 10 25 50 75 90\n",
       "Data variables:\n",
       "    tas                 (time, latitude, longitude, percentile) float64 1.586...\n",
       "Attributes:\n",
       "    description:  Contains modified ETHZ ClimWIP data used for Atlas in EUCP ...\n",
       "    history:      original ETHZ ClimWIP data files eur_pr_41-60_djf_cmip6_*pe...</pre><div class='xr-wrap' hidden><div class='xr-header'><div class='xr-obj-type'>xarray.Dataset</div></div><ul class='xr-sections'><li class='xr-section-item'><input id='section-f987b0a6-b2da-43bc-81e4-7d73e3fe8bff' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-f987b0a6-b2da-43bc-81e4-7d73e3fe8bff' class='xr-section-summary'  title='Expand/collapse section'>Dimensions:</label><div class='xr-section-inline-details'><ul class='xr-dim-list'><li><span class='xr-has-index'>time</span>: 2</li><li><span class='xr-has-index'>latitude</span>: 19</li><li><span class='xr-has-index'>longitude</span>: 20</li><li><span class='xr-has-index'>percentile</span>: 5</li><li><span class='xr-has-index'>climatology_bounds</span>: 4</li></ul></div><div class='xr-section-details'></div></li><li class='xr-section-item'><input id='section-bc90c9c3-4253-4551-96a7-4221f6650eed' class='xr-section-summary-in' type='checkbox'  checked><label for='section-bc90c9c3-4253-4551-96a7-4221f6650eed' class='xr-section-summary' >Coordinates: <span>(5)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><ul class='xr-var-list'><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>climatology_bounds</span></div><div class='xr-var-dims'>(climatology_bounds)</div><div class='xr-var-dtype'>datetime64[ns]</div><div class='xr-var-preview xr-preview'>2050-06-01 ... 2051-03-01</div><input id='attrs-561ea5c1-751a-40f5-afcb-6c74e44ff90d' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-561ea5c1-751a-40f5-afcb-6c74e44ff90d' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-0f268075-45cc-4460-a061-69c47ddb329c' class='xr-var-data-in' type='checkbox'><label for='data-0f268075-45cc-4460-a061-69c47ddb329c' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([&#x27;2050-06-01T00:00:00.000000000&#x27;, &#x27;2050-09-01T00:00:00.000000000&#x27;,\n",
       "       &#x27;2050-12-01T00:00:00.000000000&#x27;, &#x27;2051-03-01T00:00:00.000000000&#x27;],\n",
       "      dtype=&#x27;datetime64[ns]&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>time</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>datetime64[ns]</div><div class='xr-var-preview xr-preview'>2050-07-16 2051-01-16</div><input id='attrs-a54f440c-2b91-4341-9d28-93c2065bc06f' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-a54f440c-2b91-4341-9d28-93c2065bc06f' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-b380cd35-d6cc-469b-a2e7-5fc4764e5504' class='xr-var-data-in' type='checkbox'><label for='data-b380cd35-d6cc-469b-a2e7-5fc4764e5504' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([&#x27;2050-07-16T00:00:00.000000000&#x27;, &#x27;2051-01-16T00:00:00.000000000&#x27;],\n",
       "      dtype=&#x27;datetime64[ns]&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>latitude</span></div><div class='xr-var-dims'>(latitude)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>31.25 33.75 36.25 ... 73.75 76.25</div><input id='attrs-ffdc8391-260a-44c2-ae98-1ebb1a924fbf' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-ffdc8391-260a-44c2-ae98-1ebb1a924fbf' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-10d64347-e4d1-4d4b-b512-38d8f00ccdbd' class='xr-var-data-in' type='checkbox'><label for='data-10d64347-e4d1-4d4b-b512-38d8f00ccdbd' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([31.25, 33.75, 36.25, 38.75, 41.25, 43.75, 46.25, 48.75, 51.25, 53.75,\n",
       "       56.25, 58.75, 61.25, 63.75, 66.25, 68.75, 71.25, 73.75, 76.25])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>longitude</span></div><div class='xr-var-dims'>(longitude)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>-8.75 -6.25 -3.75 ... 36.25 38.75</div><input id='attrs-edbf5c8e-a74f-4a86-b8ab-37dba030680d' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-edbf5c8e-a74f-4a86-b8ab-37dba030680d' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-2f5b446c-e3f0-4c57-a9ba-43eb86200f7e' class='xr-var-data-in' type='checkbox'><label for='data-2f5b446c-e3f0-4c57-a9ba-43eb86200f7e' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([-8.75, -6.25, -3.75, -1.25,  1.25,  3.75,  6.25,  8.75, 11.25, 13.75,\n",
       "       16.25, 18.75, 21.25, 23.75, 26.25, 28.75, 31.25, 33.75, 36.25, 38.75])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>percentile</span></div><div class='xr-var-dims'>(percentile)</div><div class='xr-var-dtype'>int64</div><div class='xr-var-preview xr-preview'>10 25 50 75 90</div><input id='attrs-272eae0e-9faa-487c-830a-e895140e6b4c' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-272eae0e-9faa-487c-830a-e895140e6b4c' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-315940da-0323-43a2-8ef8-8ac3de85c510' class='xr-var-data-in' type='checkbox'><label for='data-315940da-0323-43a2-8ef8-8ac3de85c510' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([10, 25, 50, 75, 90])</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-9c843791-4398-4ab4-8484-090bde728c03' class='xr-section-summary-in' type='checkbox'  checked><label for='section-9c843791-4398-4ab4-8484-090bde728c03' class='xr-section-summary' >Data variables: <span>(1)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><ul class='xr-var-list'><li class='xr-var-item'><div class='xr-var-name'><span>tas</span></div><div class='xr-var-dims'>(time, latitude, longitude, percentile)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-017ebf6e-23dd-4ce2-bbe4-86209350fdbc' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-017ebf6e-23dd-4ce2-bbe4-86209350fdbc' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-905d2530-05f5-43ec-bdd2-87ed606bf05f' class='xr-var-data-in' type='checkbox'><label for='data-905d2530-05f5-43ec-bdd2-87ed606bf05f' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>description :</span></dt><dd>Change in Air Temperature</dd><dt><span>standard_name :</span></dt><dd>Change in Air Temperature</dd><dt><span>long_name :</span></dt><dd>Change in Near-Surface Air Temperature</dd><dt><span>units :</span></dt><dd>K</dd><dt><span>cell_methods :</span></dt><dd>time: mean changes over 20 years 2041-2060 vs 1995-2014</dd></dl></div><div class='xr-var-data'><pre>array([[[[1.586   , ..., 2.587486],\n",
       "         ...,\n",
       "         [2.176319, ..., 3.586963]],\n",
       "\n",
       "        ...,\n",
       "\n",
       "        [[     nan, ...,      nan],\n",
       "         ...,\n",
       "         [     nan, ...,      nan]]],\n",
       "\n",
       "\n",
       "       [[[1.518471, ..., 2.417952],\n",
       "         ...,\n",
       "         [1.553173, ..., 2.943531]],\n",
       "\n",
       "        ...,\n",
       "\n",
       "        [[     nan, ...,      nan],\n",
       "         ...,\n",
       "         [     nan, ...,      nan]]]])</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-6c4e48fa-22ef-4dda-ad86-81c9d1de78f3' class='xr-section-summary-in' type='checkbox'  checked><label for='section-6c4e48fa-22ef-4dda-ad86-81c9d1de78f3' class='xr-section-summary' >Attributes: <span>(2)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><dl class='xr-attrs'><dt><span>description :</span></dt><dd>Contains modified ETHZ ClimWIP data used for Atlas in EUCP project.</dd><dt><span>history :</span></dt><dd>original ETHZ ClimWIP data files eur_pr_41-60_djf_cmip6_*perc_rel.nc, eur_tas_41-60_djf_cmip6_*perc.nc</dd></dl></div></li></ul></div></div>"
      ],
      "text/plain": [
       "<xarray.Dataset>\n",
       "Dimensions:             (time: 2, latitude: 19, longitude: 20, percentile: 5, climatology_bounds: 4)\n",
       "Coordinates:\n",
       "  * climatology_bounds  (climatology_bounds) datetime64[ns] 2050-06-01 ... 20...\n",
       "  * time                (time) datetime64[ns] 2050-07-16 2051-01-16\n",
       "  * latitude            (latitude) float64 31.25 33.75 36.25 ... 73.75 76.25\n",
       "  * longitude           (longitude) float64 -8.75 -6.25 -3.75 ... 36.25 38.75\n",
       "  * percentile          (percentile) int64 10 25 50 75 90\n",
       "Data variables:\n",
       "    tas                 (time, latitude, longitude, percentile) float64 ...\n",
       "Attributes:\n",
       "    description:  Contains modified ETHZ ClimWIP data used for Atlas in EUCP ...\n",
       "    history:      original ETHZ ClimWIP data files eur_pr_41-60_djf_cmip6_*pe..."
      ]
     },
     "execution_count": 16,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "ds = xr.open_dataset(output_path / \"atlas_EUCP_ETHZ_CMIP6_ClimWIP_cons_tas.nc\")\n",
    "ds"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "interpreter": {
   "hash": "90e02b2587fbb2ca467fc85381f5522fddb4a9e5fbb8605712260c849ecf752b"
  },
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.9.6"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}