Open RADOLAN RW with xarray backend, georeference, plot and select

radolan-rw.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Imports"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "import xarray as xr\n",
    "import wradlib as wrl"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Get file and open with RADOLAN xarray backend"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "fname = \"/home/kai/projects/data/wradlib-data/radolan/showcase/raa01-rw_10000-1408102050-dwd---bin.gz\""
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/kai/mambaforge/envs/wradlib/lib/python3.10/site-packages/xarray/conventions.py:521: SerializationWarning: variable 'RW' has multiple fill values {2490, 2500, 65535}, decoding all values to NaN.\n",
      "  new_vars[k] = decode_cf_variable(\n"
     ]
    }
   ],
   "source": [
    "ds = xr.open_dataset(fname, engine=\"radolan\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Inspect RADOLAN xarray Dataset"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "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[data-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 !important;\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:  (y: 900, x: 900, time: 1)\n",
       "Coordinates:\n",
       "  * time     (time) datetime64[ns] 2014-08-10T20:50:00\n",
       "  * y        (y) float64 -4.658e+06 -4.657e+06 ... -3.76e+06 -3.759e+06\n",
       "  * x        (x) float64 -5.23e+05 -5.22e+05 -5.21e+05 ... 3.75e+05 3.76e+05\n",
       "Data variables:\n",
       "    RW       (y, x) float32 ...\n",
       "Attributes:\n",
       "    radarid:         10000\n",
       "    formatversion:   3\n",
       "    radolanversion:  2.13.1\n",
       "    radarlocations:  [&#x27;boo&#x27;, &#x27;ros&#x27;, &#x27;emd&#x27;, &#x27;hnr&#x27;, &#x27;umd&#x27;, &#x27;pro&#x27;, &#x27;ess&#x27;, &#x27;asd&#x27;,...</pre><div class='xr-wrap' style='display:none'><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-f04882d1-dc6c-423d-b767-054b4846c603' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-f04882d1-dc6c-423d-b767-054b4846c603' 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'>y</span>: 900</li><li><span class='xr-has-index'>x</span>: 900</li><li><span class='xr-has-index'>time</span>: 1</li></ul></div><div class='xr-section-details'></div></li><li class='xr-section-item'><input id='section-d0cff812-09f2-455a-a59f-4a6913852e79' class='xr-section-summary-in' type='checkbox'  checked><label for='section-d0cff812-09f2-455a-a59f-4a6913852e79' class='xr-section-summary' >Coordinates: <span>(3)</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'>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'>2014-08-10T20:50:00</div><input id='attrs-a525e666-7c82-4431-8e94-a91ba9f9054f' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-a525e666-7c82-4431-8e94-a91ba9f9054f' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-e942042b-fed4-4fe4-9112-bb3c79c5713b' class='xr-var-data-in' type='checkbox'><label for='data-e942042b-fed4-4fe4-9112-bb3c79c5713b' 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>time</dd></dl></div><div class='xr-var-data'><pre>array([&#x27;2014-08-10T20:50:00.000000000&#x27;], dtype=&#x27;datetime64[ns]&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>y</span></div><div class='xr-var-dims'>(y)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>-4.658e+06 ... -3.759e+06</div><input id='attrs-5a67e1d9-690b-448e-be86-a9cfcada5fcf' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-5a67e1d9-690b-448e-be86-a9cfcada5fcf' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-27a6a2a1-66ee-44cf-8b67-654799bfde76' class='xr-var-data-in' type='checkbox'><label for='data-27a6a2a1-66ee-44cf-8b67-654799bfde76' 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>units :</span></dt><dd>m</dd><dt><span>long_name :</span></dt><dd>y coordinate of projection</dd><dt><span>standard_name :</span></dt><dd>projection_y_coordinate</dd></dl></div><div class='xr-var-data'><pre>array([-4658144.724266, -4657144.724266, -4656144.724266, ..., -3761144.724266,\n",
       "       -3760144.724266, -3759144.724266])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>x</span></div><div class='xr-var-dims'>(x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>-5.23e+05 -5.22e+05 ... 3.76e+05</div><input id='attrs-b540d87b-bb3b-4b9e-80ad-51757b3efde3' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-b540d87b-bb3b-4b9e-80ad-51757b3efde3' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-4795d687-529f-4b97-aedd-6ed656a98ed6' class='xr-var-data-in' type='checkbox'><label for='data-4795d687-529f-4b97-aedd-6ed656a98ed6' 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>units :</span></dt><dd>m</dd><dt><span>long_name :</span></dt><dd>x coordinate of projection</dd><dt><span>standard_name :</span></dt><dd>projection_x_coordinate</dd></dl></div><div class='xr-var-data'><pre>array([-522962.166922, -521962.166922, -520962.166922, ...,  374037.833078,\n",
       "        375037.833078,  376037.833078])</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-ef2fa023-ffc5-4a90-a344-1eabb0c69d53' class='xr-section-summary-in' type='checkbox'  checked><label for='section-ef2fa023-ffc5-4a90-a344-1eabb0c69d53' 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>RW</span></div><div class='xr-var-dims'>(y, x)</div><div class='xr-var-dtype'>float32</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-cb84ec86-4831-4f0f-94b3-3127431bbeec' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-cb84ec86-4831-4f0f-94b3-3127431bbeec' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-7812919b-938b-4931-8831-d077e3671b7a' class='xr-var-data-in' type='checkbox'><label for='data-7812919b-938b-4931-8831-d077e3671b7a' 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>valid_min :</span></dt><dd>0</dd><dt><span>valid_max :</span></dt><dd>4095</dd><dt><span>standard_name :</span></dt><dd>rainfall_rate</dd><dt><span>long_name :</span></dt><dd>RW</dd><dt><span>unit :</span></dt><dd>mm h-1</dd></dl></div><div class='xr-var-data'><pre>[810000 values with dtype=float32]</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-b29700d0-8a79-49f1-bfaa-4bee79b64b43' class='xr-section-summary-in' type='checkbox'  checked><label for='section-b29700d0-8a79-49f1-bfaa-4bee79b64b43' class='xr-section-summary' >Attributes: <span>(4)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><dl class='xr-attrs'><dt><span>radarid :</span></dt><dd>10000</dd><dt><span>formatversion :</span></dt><dd>3</dd><dt><span>radolanversion :</span></dt><dd>2.13.1</dd><dt><span>radarlocations :</span></dt><dd>[&#x27;boo&#x27;, &#x27;ros&#x27;, &#x27;emd&#x27;, &#x27;hnr&#x27;, &#x27;umd&#x27;, &#x27;pro&#x27;, &#x27;ess&#x27;, &#x27;asd&#x27;, &#x27;neu&#x27;, &#x27;nhb&#x27;, &#x27;oft&#x27;, &#x27;tur&#x27;, &#x27;isn&#x27;, &#x27;fbg&#x27;, &#x27;mem&#x27;]</dd></dl></div></li></ul></div></div>"
      ],
      "text/plain": [
       "<xarray.Dataset>\n",
       "Dimensions:  (y: 900, x: 900, time: 1)\n",
       "Coordinates:\n",
       "  * time     (time) datetime64[ns] 2014-08-10T20:50:00\n",
       "  * y        (y) float64 -4.658e+06 -4.657e+06 ... -3.76e+06 -3.759e+06\n",
       "  * x        (x) float64 -5.23e+05 -5.22e+05 -5.21e+05 ... 3.75e+05 3.76e+05\n",
       "Data variables:\n",
       "    RW       (y, x) float32 ...\n",
       "Attributes:\n",
       "    radarid:         10000\n",
       "    formatversion:   3\n",
       "    radolanversion:  2.13.1\n",
       "    radarlocations:  ['boo', 'ros', 'emd', 'hnr', 'umd', 'pro', 'ess', 'asd',..."
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "display(ds)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Create RADOLAN grid in wgs84 and assign"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "radolan_grid = wrl.georef.get_radolan_grid(nrows=ds.dims[\"y\"], ncols=ds.dims[\"x\"], wgs84=True, mode=\"radolan\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "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[data-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 !important;\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:  (y: 900, x: 900, time: 1)\n",
       "Coordinates:\n",
       "  * time     (time) datetime64[ns] 2014-08-10T20:50:00\n",
       "  * y        (y) float64 -4.658e+06 -4.657e+06 ... -3.76e+06 -3.759e+06\n",
       "  * x        (x) float64 -5.23e+05 -5.22e+05 -5.21e+05 ... 3.75e+05 3.76e+05\n",
       "    lon      (y, x) float64 3.589 3.601 3.613 3.625 ... 15.66 15.67 15.69 15.7\n",
       "    lat      (y, x) float64 46.95 46.95 46.95 46.96 ... 54.74 54.73 54.73 54.73\n",
       "Data variables:\n",
       "    RW       (y, x) float32 ...\n",
       "Attributes:\n",
       "    radarid:         10000\n",
       "    formatversion:   3\n",
       "    radolanversion:  2.13.1\n",
       "    radarlocations:  [&#x27;boo&#x27;, &#x27;ros&#x27;, &#x27;emd&#x27;, &#x27;hnr&#x27;, &#x27;umd&#x27;, &#x27;pro&#x27;, &#x27;ess&#x27;, &#x27;asd&#x27;,...</pre><div class='xr-wrap' style='display:none'><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-1198588a-551c-4eff-92f2-0d241fe7fa84' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-1198588a-551c-4eff-92f2-0d241fe7fa84' 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'>y</span>: 900</li><li><span class='xr-has-index'>x</span>: 900</li><li><span class='xr-has-index'>time</span>: 1</li></ul></div><div class='xr-section-details'></div></li><li class='xr-section-item'><input id='section-bb6b075d-1cb3-4f27-bed5-f6a995eebf57' class='xr-section-summary-in' type='checkbox'  checked><label for='section-bb6b075d-1cb3-4f27-bed5-f6a995eebf57' 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'>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'>2014-08-10T20:50:00</div><input id='attrs-6274c8eb-ac92-43d4-b578-391319a80081' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-6274c8eb-ac92-43d4-b578-391319a80081' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-dd7583c5-dbaf-4c3b-ae1f-81d12884a288' class='xr-var-data-in' type='checkbox'><label for='data-dd7583c5-dbaf-4c3b-ae1f-81d12884a288' 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>time</dd></dl></div><div class='xr-var-data'><pre>array([&#x27;2014-08-10T20:50:00.000000000&#x27;], dtype=&#x27;datetime64[ns]&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>y</span></div><div class='xr-var-dims'>(y)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>-4.658e+06 ... -3.759e+06</div><input id='attrs-ac9bf6ab-d45b-4fc2-95fd-082ec7de1841' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-ac9bf6ab-d45b-4fc2-95fd-082ec7de1841' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-67243869-54d2-4fd1-b2bf-fa108c43c624' class='xr-var-data-in' type='checkbox'><label for='data-67243869-54d2-4fd1-b2bf-fa108c43c624' 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>units :</span></dt><dd>m</dd><dt><span>long_name :</span></dt><dd>y coordinate of projection</dd><dt><span>standard_name :</span></dt><dd>projection_y_coordinate</dd></dl></div><div class='xr-var-data'><pre>array([-4658144.724266, -4657144.724266, -4656144.724266, ..., -3761144.724266,\n",
       "       -3760144.724266, -3759144.724266])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>x</span></div><div class='xr-var-dims'>(x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>-5.23e+05 -5.22e+05 ... 3.76e+05</div><input id='attrs-c93b46c2-ce93-41c6-8453-c8f241fce467' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-c93b46c2-ce93-41c6-8453-c8f241fce467' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-93b36881-9382-47ab-8edb-d37946e1668c' class='xr-var-data-in' type='checkbox'><label for='data-93b36881-9382-47ab-8edb-d37946e1668c' 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>units :</span></dt><dd>m</dd><dt><span>long_name :</span></dt><dd>x coordinate of projection</dd><dt><span>standard_name :</span></dt><dd>projection_x_coordinate</dd></dl></div><div class='xr-var-data'><pre>array([-522962.166922, -521962.166922, -520962.166922, ...,  374037.833078,\n",
       "        375037.833078,  376037.833078])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>lon</span></div><div class='xr-var-dims'>(y, x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>3.589 3.601 3.613 ... 15.69 15.7</div><input id='attrs-0882cc2e-238b-4565-8256-234041b73c0e' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-0882cc2e-238b-4565-8256-234041b73c0e' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-319d4e4a-4349-4b2c-af15-2bf9179c5a73' class='xr-var-data-in' type='checkbox'><label for='data-319d4e4a-4349-4b2c-af15-2bf9179c5a73' 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([[ 3.58892995,  3.60107571,  3.61322204, ..., 14.58426285,\n",
       "        14.59648289, 14.6087025 ],\n",
       "       [ 3.58756495,  3.59971325,  3.61186212, ..., 14.5852429 ,\n",
       "        14.59746553, 14.60968773],\n",
       "       [ 3.58619938,  3.59835022,  3.61050164, ..., 14.58622337,\n",
       "        14.59844858, 14.61067338],\n",
       "       ...,\n",
       "       [ 2.07772252,  2.09266529,  2.10760913, ..., 15.67097866,\n",
       "        15.68606112, 15.70114278],\n",
       "       [ 2.07564265,  2.09058924,  2.10553691, ..., 15.67247681,\n",
       "        15.6875632 , 15.70264879],\n",
       "       [ 2.0735617 ,  2.08851211,  2.10346361, ..., 15.67397575,\n",
       "        15.68906607, 15.70415559]])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>lat</span></div><div class='xr-var-dims'>(y, x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>46.95 46.95 46.95 ... 54.73 54.73</div><input id='attrs-14f244b5-9e27-4b58-86fe-94991fd6725b' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-14f244b5-9e27-4b58-86fe-94991fd6725b' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-79e381d1-016a-4bea-9df6-dac26f337bc0' class='xr-var-data-in' type='checkbox'><label for='data-79e381d1-016a-4bea-9df6-dac26f337bc0' 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([[46.95258041, 46.95351109, 46.95444002, ..., 47.07247581,\n",
       "        47.07180758, 47.07113758],\n",
       "       [46.96087117, 46.9618021 , 46.96273128, ..., 47.08079892,\n",
       "        47.08013051, 47.07946033],\n",
       "       [46.96916238, 46.97009356, 46.97102299, ..., 47.08912249,\n",
       "        47.0884539 , 47.08778353],\n",
       "       ...,\n",
       "       [54.56171718, 54.5629217 , 54.56412397, ..., 54.71701774,\n",
       "        54.71615147, 54.71528292],\n",
       "       [54.57038125, 54.57158615, 54.57278879, ..., 54.72573003,\n",
       "        54.72486349, 54.72399466],\n",
       "       [54.5790457 , 54.58025097, 54.58145398, ..., 54.7344427 ,\n",
       "        54.7335759 , 54.7327068 ]])</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-b1b35797-51ad-4a13-8d4f-5d3a1edf845f' class='xr-section-summary-in' type='checkbox'  checked><label for='section-b1b35797-51ad-4a13-8d4f-5d3a1edf845f' 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>RW</span></div><div class='xr-var-dims'>(y, x)</div><div class='xr-var-dtype'>float32</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-ece1e501-277e-48f9-bced-281e46f852a0' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-ece1e501-277e-48f9-bced-281e46f852a0' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-8892bf59-5417-4eb7-8536-6b6c4c3aef6d' class='xr-var-data-in' type='checkbox'><label for='data-8892bf59-5417-4eb7-8536-6b6c4c3aef6d' 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>valid_min :</span></dt><dd>0</dd><dt><span>valid_max :</span></dt><dd>4095</dd><dt><span>standard_name :</span></dt><dd>rainfall_rate</dd><dt><span>long_name :</span></dt><dd>RW</dd><dt><span>unit :</span></dt><dd>mm h-1</dd></dl></div><div class='xr-var-data'><pre>[810000 values with dtype=float32]</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-a13a96ed-505e-4e8f-be17-8f37a8b7d336' class='xr-section-summary-in' type='checkbox'  checked><label for='section-a13a96ed-505e-4e8f-be17-8f37a8b7d336' class='xr-section-summary' >Attributes: <span>(4)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><dl class='xr-attrs'><dt><span>radarid :</span></dt><dd>10000</dd><dt><span>formatversion :</span></dt><dd>3</dd><dt><span>radolanversion :</span></dt><dd>2.13.1</dd><dt><span>radarlocations :</span></dt><dd>[&#x27;boo&#x27;, &#x27;ros&#x27;, &#x27;emd&#x27;, &#x27;hnr&#x27;, &#x27;umd&#x27;, &#x27;pro&#x27;, &#x27;ess&#x27;, &#x27;asd&#x27;, &#x27;neu&#x27;, &#x27;nhb&#x27;, &#x27;oft&#x27;, &#x27;tur&#x27;, &#x27;isn&#x27;, &#x27;fbg&#x27;, &#x27;mem&#x27;]</dd></dl></div></li></ul></div></div>"
      ],
      "text/plain": [
       "<xarray.Dataset>\n",
       "Dimensions:  (y: 900, x: 900, time: 1)\n",
       "Coordinates:\n",
       "  * time     (time) datetime64[ns] 2014-08-10T20:50:00\n",
       "  * y        (y) float64 -4.658e+06 -4.657e+06 ... -3.76e+06 -3.759e+06\n",
       "  * x        (x) float64 -5.23e+05 -5.22e+05 -5.21e+05 ... 3.75e+05 3.76e+05\n",
       "    lon      (y, x) float64 3.589 3.601 3.613 3.625 ... 15.66 15.67 15.69 15.7\n",
       "    lat      (y, x) float64 46.95 46.95 46.95 46.96 ... 54.74 54.73 54.73 54.73\n",
       "Data variables:\n",
       "    RW       (y, x) float32 ...\n",
       "Attributes:\n",
       "    radarid:         10000\n",
       "    formatversion:   3\n",
       "    radolanversion:  2.13.1\n",
       "    radarlocations:  ['boo', 'ros', 'emd', 'hnr', 'umd', 'pro', 'ess', 'asd',..."
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "dsg = ds.assign_coords({\"lon\": ([\"y\", \"x\"], radolan_grid[..., 0]),\n",
    "                        \"lat\": ([\"y\", \"x\"], radolan_grid[..., 1])})\n",
    "display(dsg)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Simple plot in latlon-space"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.collections.QuadMesh at 0x7f475c13eb30>"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXkAAAEGCAYAAACAd+UpAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8/fFQqAAAACXBIWXMAAAsTAAALEwEAmpwYAAB9LklEQVR4nO39e5h1+1XXiX7Gb17Wta7vbd+ys3dC5NJoAmxpJWpzMTEognoaEcVOt3bHVjitHjkI2q22tv3gEUWf0/2oW6CNShBRc0DahmxRpFEEEwghkEBI2DvZe7/3ute6zTl/4/wxfnOuVfXW9X1rVdVa7/w+z3qq1m2u31q16vsbc4zv+A5RVWrUqFGjxnzCXfQCatSoUaPG9FCTfI0aNWrMMWqSr1GjRo05Rk3yNWrUqDHHqEm+Ro0aNeYY8UUv4CS4evWqPvfccxe9jBqXFL/yoU890vNFBBopDEccqzYTscfvh3MgAgLEEeQFeAVV3vK2Nz7S+mqcDh/60Ifuqeq1RznG7/iyjt5fK45/rY8Mf1RV3/UorzVtzATJP/fcc3zwgx+86GXUmCLe4b72+AfJxImn+urX/1yef7gXFQfqia9cgeVF/Gdew49GRz8lTkA90miAKn4wJOq0katXYDiEPEdHGdKO0aKAokB/MQOwY0+s+6Lwkv+Bi17CVCEirzzqMe6tFfz0jz5z7OOSJz959VFfa9qYCZKvMd94gOAD+T4A9UgUGXmeBcJr+J1dHOCz/Pin5EbYUnj8cIBrt43w+31QBa9It0Nx+y6SxGiWI07s2JeA4OGEG+oRmPdNwqAUl+Tv9aioSb7GheJAwtn/zzVB+uqn0LwXOfL79x/66X57G2m3xxvQ6iKRKtrrW8oGECfoGe1NNaYPBTzz0Shak3yNC8OJI0r1e1M1cHi0f1qIw/d6p36aH41AHDocEi0vo+rxO7t2yNdu43f7oL7alKayOdWYKjzzEcnX6poaswH1iJOzT3kctIGc8HkuTYmeuGHR+yir8vk6GBrBFwWaZ3b/AeuWOBm/9sOs4YLwqOmeWYCiZOqPvcwC6ki+xoE4yT/yo+RmH4YotCjGZHiW/2CHHeuYswU/HMC9NYgcOsrsseoBZ1F+cXgRV+IE12mjo9E4b39M0bfG+UGBok7X1Jh1HEu0R5GcuCOfP0vFOYmTqqC6B8dtJOKQVhPt9apisM/yB4u4+44jUYREEXhvhdooQpoNCCkgiaKD11PjXFHn5GvMHE4cPZfkfhTJHUOA+19rkvQf+nT/rPLw+/CwhOrSFFT3qn3UHx+Ri0M6LQiKG2k00NGoIv8zUw9NEe9wXztTG/lpoUAxJw69NcnPOU5N7DAVIj2TPO5ly4E6we/2Tk3KkiZQeHQ4RL3i0hQtCpNiFh71x0s59xwvToCH36xqHIxL9m17aNQkP6c4NaleNgKdBXg9NbGW0XqxtW2F26VOdXtZqD3t3+KiyH2eo3lF65x8jcuHEzcV1XhkWKrGn/ozVm/pHZeEfz0n1ZmA63bI+wO7vf7bXShUIZsPjq9Jfh7wyKmQmlBOBzH1zKltCsQRdS1y98Mh0ojRwZBie5v4yhX85vZDbRwXhXmN4g1CwQEeRTOImuRnGEeS+6T+uiSMGSGPS4+yKH1KgndJjAY9vUtT3EIXv7FF1O3i+32TZM4w5kVtBaHjtY7ka1wkThK9i5O9nZazSPCXcGOSODmRz80ehO5XzY3IpdXEb20b4Ucmm5x87Cxg3pui6ki+xoVgKgXVy0gq02h6OiMc1sF6kudJnFSE7vt9u56m+O3tQ5/nGs2Zj/JnDdYMNR8kPzu91DVOR/Chrf5MMUOt99NCKVd8mOdVtgzqIc9xaWpnW8PhkX+rkxL8rKVELjMUyNQde5kFzMYqa5zvqfFBZmBw/lH1ZYziH0LiaE+c2HTFoVk+Njk7AzuDMyH4ehOvoAgF7tjLcRCRpoj8jIj8vIj8ooj8z+H2vyQir4nIh8Pld07rvdR/1RnAmRG8uJNFovtJ7LzJdn8ePqhZTvS8SwaJIrNNmCB4mIjOz+iMq/yOPBLZX8JN9SLhVY69nABD4MtV9a3A24B3ichvCvd9p6q+LVz+1ZTeRp2Tv+w4E4KfiMRn0tM8DAtBknE+fPLs4rLn72Vc/JYkxg+G9nsUVQXc+MknyG/eeqTXmvdC6HnirHLyavMkd8LVJFzOVbdz+UKfGhXO7J92Uu53CYnwARw0NCQQuUtTc3BM4orgp2JBfJaYkLD6ft/Wy94C7qMSfI2zhlCoO/ZyoiOJRCLyYeAO8JKq/nS465tE5CMi8j0isjKlN1KT/GXFVKKyg9IZlzDFcRC0KCpydM2GRcBh7bNg6DXZt7A/dVPj8sEmQ7ljL8BVEfngxOU9DxxLtVDVtwHPAF8sIp8P/B3gzVgK5ybwN6b1XqaarhGRl4FtoAByVX1h4r5vBv46cE1V701zHbOGczvtvoQa9AMxsUZpNdHRCNdqVgZfhw34vjQ4bH2Xca01AFAVRhod/0C4N8lrRx9TN0Tkx4F3qep3lLeLyN8HfvihFnoCnEdO/sv2k7iIvAF4B/Dpc3j9mcLUCP6gYuasIShRJE3tahTh+4OxPa961F/CjWuGrApqjOHPICcvIteALBB8C/jtwF8TkSdV9WZ42O8FPvrIL3YILqrw+p3AtwA/eEGv//hhH8G4JJ6tSUSBJCUK6qCiQNIU5xx+t/fgBnbZCPWyrafGkbDC65kEQk8C7xWRCEuP/1NV/WER+Uci8rbwUi8Df+wsXuwgTJvkFfiAiCjw91T1RRH5auA1Vf15kcN3ypDbeg/As88+O+VlXg6cV5rGNZqzRTrltKQwSxWXm/e6EyRNkEEELoagZJmJHP1jhNls0pITF1aPgqp+BPiCA27/w4988BNi2iT/dlV9XUSuAy+JyMeBPw+887gnquqLwIsAL7zwwpxYBR2O85S/SZrg+4PLGfHug0tTiGMoPLgYzXIoQMpBG+UIvtGosvEVJ7MpFZ1TTH63Z4Xwy8LrPGCqJK+qr4efd0Tk/cB/ATwPlFH8M8DPisgXq2qtITsPBHXHw/qvnDek20EHQ4vaowgk/IwcEt6LpIkZgKvaBtbrXfSyaxyCWSL84mTNTpceUyN5EekATlW3w+/vBP6yql6feMzLwAuPu7rmPKP4aHHBSHMGCD5aWDAixyJ3fFhzFFUED0DhkUaKZhk6ymbiDOVhsd+szDWap/e1vyQ4ag7wRUMRMp2PXtFpvosbwPtDxB4D71PVH5ni69U4CSI3EwXXygqgKHCtlk1Q6g+QZsPG7gVSkzSB4OooSSjKjrB0zQwS31FwjebeUX8yG3/Lk+JMrBnOCGdYeL1wTI3kVfVTwFuPecxz03r9WcF5RvEuTfFbOzNBfiWZuW4HRNBe33LzgDRSyAvUeUvTjLLqdg0dpRKlc2fP67od8vv3qw3QpSmaZ2dbaD7MImKy+D3l789lIHtF6nRNjdlDNUBkRtIZWhQ2Es8JiEBuqhrNMshyiCLU56aqiSIYjWwjUDWJZUlKsw5xREuL1gTWbuP7A6JuZ9wMdlYvE7zuD90cg4fQeZ0lXTTZ14XXGo+Mc4niA6FHqyuQ5RQ7uzNB8EC1Ts1yXLMBcYx02mi/D87+ASWJUSdG+qpBhSOhwDwfqYzSc15aTTtrCc6VZ5qqEXf0WUFlbleMvYTO6Xt0EWSvyplIKC8D5uNd1DgULomJb1xHd3YvN8HLvhF4AX44MAsDb+StO7tIHCON1Ag+yChxDj8YVnl8c388+JizBrNuEHQwrIjd9882FVV9TpPzgMvPb3JecHiMhI30PHGeqU0rvEbHXmYBdSQ/rxBH1GlDHFPcu394hHZJUjcSRUHfPl5n2bSlIfWCV6TTtvQMmJzSOYvwM4eMRlaUDXljF+wP1OuleI8PC5emSBRRlCMCz/pvFsha4qSaRStRhCSxpclGGa67iPYHJmWN4/FazhnnGdXPS+F1Pt5FjT2QOLGcbVFQbG4dnZe+DOQXSMtn+Z6hJpImRjhpinQ7lq7Ic2i1LAcfx/bTe/tdBGk1kSiqFDl7bJZnEeIgcnYWVuKMCb5096yuR5F1EsexWUd0O3a92cAtLiCLC/b5XiCmHdUrxw8MOeHQkAtHTfIXhGl9SV2jiVvsWs62358pgnNJXE1Scq2W5aHLNEIeCqztNmzv2PXIQbuFXl22makTZKRZbimOsMHNbNpGvaVmpvB3LM+exr5Adl1C45nmNtBEFrtW6/DeNtjSJO6CTe6mTfRnMf7vMmA2VlnjeIgzSV1RUKytz1bXZ0gRqFfrZE0T86aJY/s9dLqSJNDvQ6mNbzWtQuY9LHTt/vJ65CrrBpemZ6pCOXdMieDtFze+7sQ+80ajuk37A8gLaKTwxqfAF+j1VdtgLwGmRfQKeHXHXmYBdU5+DlCmOGa9McY1G6gPEWUZSUYRpImRzSiz5qckhsEQ30zxnZRos2fknufgFZyzxqkoMr39SebaPg4IqZlJD/4qkg8yVXHOPvOiMMIXu11XF5H+CF3qIps7cP0qUZbje735kKk+ADmT8X+XATXJzyomTpVnxYfmKLjEzMfKDlYNmncdDo3ARZB2q5JO6hNXKRYbRLsjuy2JodmAUYaGs5iSfOaThE6JcEZj3vY67iGIxxQgaWpk79XqH4W3s6VuG+kNoddHO03rLO71x8+Lk72duHMAhZlRzxyHmuRnFFVENjnIekZhhmPjfygdDo3U09R8axoNSw80G+hiBxnlyMh08b4RE28P0GaC9AcWwcex1SOqA872BvjIKKWQkTOZexKi9sgZsRcFOGedxM4ZsXsFJxTXFpFBDg50pUu0YQVgXd+0M62gr583qMrMpGOOQ03ys4ZQJCt14MDZy+nOixQnpXshZ65FYb9HMvataTbQ3R6SJKhz+KUWeMg7CfFuZgSfFZBY/l5V96hGHudIvkx9WaFUbDMtvX7SBM1zJEohjiBJ0IUO4r19zotNek+2aN4bkrdjopEnujNCN7dtM7gEdY5pSinnpRmqJvlZQknA0yD3EudE8FXePRCFtJr4nV1r3W9a4a9MQ8lohCwv4Vc6yCBHvCdbbSNecZn516gIEhQ3EqVVj4A0Uvz2zmwVos8KZWdqHO/1/XEOyuKqL6Dw6Ooi2WqbaDezDRYYrjRo3R7g0wiNhHhzAIOhbQ6jzFI7BZem1+IsYX7ydU6+xnmj/Eea9X+oknyiCNe13G+xuQVY449mOZLESJkvTqzYKqMcbSX4RowbFvjUkbcT0q0BsrkdZJWRPbbZABfhb95+vCP5xLTu0kih24Ess4i93cA37AxIioJ8qYkbFRTdhLwZEWUeCYF63o1RAU0iaLfRrS3zD5poPJs/nM1kqMuAmuRrnC/E4YL0URoN/NbOWJddFGie2e9Z0GgnsaUSigLJBXyGOEfRSYgGBa6fQ+yMvESg0yZ7YpHk1XV0YzPkoYu5jDaPhDiibsc2uzS1PHtsaRq/3MEnERoLDvDthKybEGWeInGIV3xsP7fe3EY8LLzcJ1rbhXBGJFFEMSX9/mWASSjrSL5GjdNB3LjhKU3xO7vjAjKMid7npgQRa6Enn4jEWyniPW5U4HoZ4j1yb6OK9ovVDvFG36L5pQXY2LJ0xBwWB4+CS+Jx3lzE6hY7A7TbtPuHGZob6Q+uNWisZ2TdmGjgwYEUnqwb4zJl8RPbuO0BbG2blFXVNuE5JXgYe9fMA2qSr3FuKMldR2GCE+YrUxJ95VFjdyBp09IvQatNGPyuSUS0PUCGwXmy3UKbCdpIiO5sQeTQ9Q0j98Fw7nzlj4W4caF1oYu2GvZZZRnSd2gSoamlvxCheW+EFJ68lZJujCiaRgtZx9G8n+N6I1jbsDm65d8tn+8JXHA2VsMi0gR+AmhgfPvPVPUvisgq8P3Ac8DLwO9X1fVHfsEDMB9JpxqXH2KdrFoUNgjEiaVmkhhfpmYmVB/SaplmHoKlgZGJbzVMQpkVaDMBEUs/dJtGaj5E/VlupPQ45uPVV66dZBmiCrs9S105ByK4fobrjcgXGqiD4WqDxmYOCj61M4BopDRfXoP1TXQ0wvcHeyyOZ9Yq4gQwq2E59nICDIEvV9W3Am8D3iUivwn4VuDHVPUtwI+F61NBTfI1zgVl3l3i2BqdRCrpZKmmQX1VKKza5t1YHaKdJlIUSFbgO0GB00zoPd3GJw632UOXurC1bTn/OH7s0jQwQb6F2TuwvWusFccwHOJ6I2SUk692iPoZRTMm60ZIYdr4/rUYlysLv3TfCD7LjOAn/ebF4ZoNcwqdU5yFQZkadsLVJFwU+BrgveH29wK/ZwpvAZhyuiYM6t7GhFa5qr4gIn8Fe4MeuAP816r6+jTXUeMCMeFqCCDNhkklMzMU09HImp/Kx5aKmiyHZiCQyEEcmboj92gjQWNH0U4YXG3QvDfEjQro9ZEt03BLp43f2LqAN3zxMJtgO8shy0Nu3sPOrhm8DUbghGh3xOCJDgCt20Nc5inasRVZN3uwZdykIYIvew/K5jX1F+vwOU2NvLlQnigGvioiH5y4/qKqvjj5ABGJgA8BnwX876r60yJyQ1VvAqjqTRG5flZr34/zyMl/marem7j+11X1fwIQkf8B+AvAf38O67hwnOfQg8uCysa28BBF5PfWjCjixOwHxFX+5a7ZsLy8qlkY5HnIxSfmfDiI0UbCaLVJ0XDE/YJ0O7dCq0hQ0YgRWjA7e+wQ8vE6GuEWuraBloXtbsfSWS5CGymo0nxtC99qoImjf6NJYyMj+fQ9GNjGq6NR5d9f/d2KAs2G1W3zCLM1ONH3556qvnDksVQL4G0isgy8X0Q+/9FXeHKce+FVVSfDqw72ec49HkuCD4QgkeXZy2lGEkV7zNSqxqgs2AdXDonlBlFAo0W+0sYnEYMrMe2bI0ZLMZ1PbSI7A9jdHTtVAuT5Y9kA5ZLYTN5K9VIUWepqccHcO4chn54XwRYigUjoP9Ek2SlI7u1a9F94O8vK8z0Onnt08epRP68b6dnbGqjqhoj8OPAu4LaIPBmi+CexrMZUMO2/kAIfEJEPich7yhtF5K+KyGeAP4RF8g9ARN4jIh8UkQ/evXt3ysucLh5HgrcZq0UY4NGk2Nm1/Lh6/GhUmWW5ZG+cIWlqKZvCW5phcYHi6WtsvfU6RTNmcMVI3A0LOp/cQNMYGsk4d99uQZZR3Lt/Ee/6wuGz3Da43V6Qjubj4eZ5blr5NLWh56PMahyDnHi3oHFrx4zIcitaV8omgtd/mu6tccy9ukaOvRwHEbkWInhEpAX8duDjwA8B7w4Pezfwg9N5F9Mn+ber6hcCXwl8o4j8NgBV/fOq+gbge4FvOuiJqvqiqr6gqi9cu3ZtysucHk5E8DNuMLYf5Ri5cgDFQbNl/XCAazSrVI3mmUXxeW7E1G1DmjB6aonhtSbNtQwEup/p0X2lz/BqAxkVuNvr6K27lqJJYugP8Jvbs+0f/ygIm2iJyic+DAAhSeyzUjW3ybUt3MY2zZs7yCCzz38wHNsN57kdI46rprIqmp9jgj9Ddc2TwL8VkY8A/wl4SVV/GPh24B0i8gngHeH6VDDVdE1ZUFXVOyLyfuCLMc1oifcB/yfwF6e5jksP9dOzaz3vaKt0JQzWBXsIZ+I9VhtBGgZ9iKuid2m3LApd6hBvDHBZyuabWyx8ZohPIyTztD+5YXa3o5GRVpqi2ztQBJKbYwI6Eco+A7AzozQ1gs+LcXfwbhB9tFsWwQ+Gdl9In/nhcOx5kyRz3eF6EM4iXaOqHwG+4IDb7wNf8cgvcAJMLYQUkY6ILJS/A+8EPioib5l42Fdjpy5zi5MoAKpIdho473/KYKDmkviBodySxLh2G4kTi/LbbZNLeg0+8uHrWBQWfQL5SpOdZ5ssvjwg6ue4fk58fwdtJZBlptSpCq5Ksbn5WBHRoSjTZWk6Ts/0B0bi3oc6h00SY33Tbo+jcT9Cv28F28KbE6jIkZ9rOTR9XjBPM16nGcnfwCrJ5eu8T1V/RET+uYh8NiahfIXHRFlzFOatYUeiyHLDgRQkTsyDRgQdZbhmwwqBqmivPxHNj61wdaGNjAqGywndV/r4NGJ4tUnz5i69N6/S/rlX8JvbSKeFdDuMnr9O9JM/f8Hv/BIicuPPtijGaZjwU9IEHVrHK2XzWEjtaFFY81pz1XoPDsH+Qvo8QIG8Nig7Gqr6KeCtB9z+/5jWa15WvOR/4PEpvu4fIiEO12mD9yGyTEzKB2PJYxKPpxF5b41Sa1votRW6n9wEEUbLXdzQo0lE+2O3wDnc8iKjz32GeHNA8uFPUszZZvkoECdmThZFRtol0eelciYLufYwTGQwnGh0ss04Wlq02blZRrG9fWDqr1RQzSPmZWjIfLyLGcCBaRtxuEYT126f/4KmhQkSkCgyJ8QQqUuaWgS/2A0j+xLk2pXKwqA0vyILefv7m8ggQ2NH69Vtmp/ZILq9CeIonr7G9pc8T/qJW/DLL1uapoYh2DiT5TDK7LOOIvupwSOoEdIreWHDWMoxi0VhG0CaQKcFvsDfW7PHHpCuKRVT541pNkIBcIJUTZ2uqfEAyi/m7+i++4H7XLs9d7pu124jC110a9s84uMwh3UwtIi93bLiaeTQfl4No9DRyOoUowxpNXG9Ntrvw/ISutBidL3D1htTrv70fYrbdy0inXM536lRFFaMLgrER9UZUmUCpxOGcHmBdNoU99eRUEtxrSYkCcUrrz6W1hD10JAaNY6BC0VV3dq29vfIWTomD9ElwOaW+ZwPR8HTJqtSC6UZliM0SzUa+IU2wyc6SKFc+zevU7z6ujlYztPgijPYrMSJHSfISk0SaZYRUsp1Q4FV1dvm6803CDAvoaLAf+a1g9Mzjwnpz0qkfhxqkp8y3pn+wT3Xq4HKB2BeonmJE6TZQIdD/GBoJN1pVTr2Kh2QpujOrjXlXL8KG5tIu4WWjVNgjpWNFH9thWIhJd0YIh/+lT3597nJCZ9Rv4R6HcegweETEZvnWm6mQb0kElI4owy3vGSNVFs7YwfQ/cd+TAi+HhpS40TYT/BghCRg/2Ty4Jdopom+HMwdCqk6GuEWu1ZUXVk2gyzGLom6uYV6j1tahM3tSrJX2h+4RtM2xPA5xVtD9OOfDMqQ8ZCRuYni4fj3MjHnt4zYHyDectD7yAzhtLQbFrHZri43ck/Cd7AoJZUN/OaWPX84NI//x7SpTBHyObFsqEl+SjiI4EtoUSBuDj/60NSFiBVRxVmUvrpsBO/s+qTPu0RRdRtRZNFmFKEFlqu/sgKRIJ/6jCk8CGcKUWRkNE+R5Uk2q3KQu3qi60/g1zYO/AzKz1d3e/acMEREwyzdscLJPGr85zxP9Nrd8BI6Pjuapw30lKhz8jUOxVEEX8GraZvmKZoPkbyOMjTPLC9//aq911YL3dpr/WuF1agifSF41zgBFUvV3F/Hb+9Yc075vBC9P45RpkQR6pVoedk2xcPm16pH83HEbx7+48jUb2xW06Pk+lWiz9w2a+Zo3LFsZwqPTw5+D3R+0jXzcT4yg9CiqCxgD8LUJWJTQtlAE3W7po8fDGE4hF4PabVsYlNR4PsDoqVFpNkwBc3VVSOdbjsYYzm0P6C4v7aH4O1F/HykaYL1w4kfE7z549UVZHkJ3dkdD1w5DOGzKqPzcli6z/Lg+hmhaYLf3hk3Qk2kgqoi7mOGMidfSyhrPBJMwrb3H+hHd957yKNnAKWlQbOBLHQr2R7DEXRa6PrmOMpvNc36dpRBM7TeNxrozdvV5lemZ+YGE2RZ+rAfWjQu6xvloHNxuFYTcQ4aKfT7yGIXXITL8uPn2B4Q6YPVT/wrr9rfpSiAaOwMGvL9pmA69budecwKiR+HmuQvCFXKIZwKz2rkfhBKSR7djrXDNxsVwaMet7BgOeGGEXtxY4losw+b2/h+f37UMhMotfxaFLhyru1hZyLlNK3gaVT+7jptMwzrtNFmguwO0PvrDzeoPGwifnPrgRm7iDO9fJZbjSVySOHPbCB6VTB/SJzH/4oiFHNSeJ2Pd3HJ8IHR+w6/syT2ff/g82B7IHFC1Gmbi2SzYRG8c2jwPdGisO7Wq6vkb7iK77bwq12iu1vorbsUd+/PJcFPplxcq3k0wTMxRKUokDixM6M0QZ+8RvHsdbQRo87Bzq75xj8k/GhUOXaWIxpdq4l79mlkccHm7XZaVZH70LTNSdM5Zc1mRv7GZ+EnfxlQk/yUsJ/oy1zoUV/w/UQ/S9G9xIkVSkvv8SSxdExeICtLFDu7RCtL+Lc8izZiop0hW5+zaMO3797Hb2/PdYFPnIQmo70+7FX0XF6Pkwfuk0YDnn8GTSJ8GuE7Ddy9Dfz2Dg/gIfLn5d/OdTtmO+E92uvb3zHUUI4scp+0NjJDNRTV+cnJ1yQ/RXxg9L6K3GfpC/4wkNRkjRLH+PUNs7XFNjd/+x7xtSvQaiJFwebnrzJ8osvCJ3fR++sU29szE92dGmX06hU/GD7g1jientUyO+ZyUDZAZTLmkP4Id38bnzqie9uVJ3xl6XBacg/fR9doVmcKrC5bDWV7Z4+PjQ6Hj8V3eD9U5djLLKDOyZ8x3hF93Zkeb2YcLPPc8umjzPK5uz1kNEKHQ9yNq/Q/5wkat3u4zR5LH82RO2v4rR2KObOoPQgueMgcuJGFwdvAuMtUQ9QsglxdRVsp0htSXF8k6hfWcxBsCQhnCJrliIsPtfyNFhasA7lMFQXPf2k1kWtX0GaCOofrmQOlDoY2ZeqY1NKxmFlPodmJ1I9DTfJniEMJ/pRf9JLUZ4bgsfxu2UKvRUHUaFCsrRM99yxrv/kJ0x03usS9Fo2PfoZibd109Hlmyo1Sz+11RklhH8ocfDIm3nGH6t73p5mZs0maWodvmdZptyiW2kR3N+x5zhH/2q2xbXMYWi5hepM/wA7YNZq4Kyv4++t7CF6iyIa2XLtCfrVL/Np9JMstau/1THJ5FumzGf5bzkqkfhxqkj8jHBvBn5DoJ0l9Vgi+wsT70/6A6Mkn2HrbDZZ/aQsZZsidNbiyQnHvPtGN69aNKQ7XblSuk3pGCo4LxWTT0UQTkmu1rPFr8rugHolSI/hOy2oZcYSmMVufvUzndbNfzt6wQvKp2+GgYS7uaIQsLaLbO0i79UDUHS0sIKvL6Ob2nr4C12riOm30xhUGTy6QbI1gt28NeF4fPXqfA6hC4WuSr3FazOyp6+kgcYJb7FI8tcriR++hr98eD/O+t4Y4Ib95m3h1hWhlCR0M0VF2ZhK9S4NSix7y7FpKFZ0AEa7dptjZxV27gr++DMMM321SNGN23tCgc3OEjApGz10jfeWeMU+eW1fqKLN8vSqaZTB0xE9cJ795CwguoFdW0Y1NK2qHzljXauK6HUaf+wY0Fpq/fAtd28B7szeYmzOpM8BZqGdE5A3APwSewKbhvaiqf1tE/hLw3wF3w0P/nKr+q0d+wQNQk/wj4ET2BfvxGBB99NQNRm+6TvSTv0B+oHkWxNevmW5+VKB5Pnfj4/ac1XhF4tDc1G7jQmHV9/vET1zHX1kCEYZPLuITYeu5hKVPZRSJg6WUxi/fMgOxskmp20GGQ3ARur6BW+iabcSGDU5x7Tbu6l6Clygiun4Nkhh/dckGpN9ZQ7d3KnKvCX4M5czSNTnwZ1T1Z8PM6w+JyEvhvu9U1e84ixc5CjXJPwQOIvcqn3wCzFKu/WGg65tE//fNAwuN0dKSFRi7HfTWnXGEOy/kcpiPTBh9SFFAsHtwS4vo1WWKhRSfRGgk3PnChKf+/YCNNze58vPbuN1hZU2tC21kmNuglWYDNraQbgdd7pqscmPTUjRt6y6uOobFmQwziSmurxC9Yn8bP7BpUJUh2bz8Dc4EZ1N4VdWbwM3w+7aIfAx4+pEPfApMVUIpIi+LyC+IyIdF5IPhtr8uIh8XkY+IyPtFZHmaazhrHBW9l63qR+Gl4vvnmuARd6gkMlpeRlaWyT/radjt4Ra6Y4npvOAQoiz9eiinMg2HIIJvJLjMk3ci+tdiGhuQtyIWPjPC9YfoK6+CKsUTq2gjwS+2oJGiTSu66mIXjWMkDw1N166gu72xY2fZWDUc4u/ex3381yjWN/FbO3YmNWMEf569I6rHX4CrIvLBict7DjueiDwHfAHw0+Gmbwo8+D0isjKt93EekfyXqeq9iesvAd+mqrmI/DXg24A/ew7rmC5CBFdNKjoEc03w8CBZBCVHdOMaJAl+tUt8c4NysPRMum0+JMpiabGxAeKIFxeIbt6HZoNoMaV9y9O8H7zz+4UNK3/yBn7JZgBLVpBdaZNuhs+s2bAIfpihOPTJq8jt+xahM+HbX34fRyO81z1nnXO1wZ4xTpiuuaeqLxz3IBHpAv8c+FOquiUifwf4K1hm6K8AfwP4Iwc87/edYA2Do/L5556uUdUPTFz9j8B/ed5rmAoCuT2O9rdHwSUx7qkn0YUWeI+7vW6qmjyn2DmgY3PeUBqNRRGu1USHQ0uplLYPzQa+26JxaweNHUW3gY8dyd0dNI3pf9ZVomFBtJMxutohvbNtvvBxjH/6Cq6f4VsJrjdChrn5+DP2ypEoQoLevfKYLz1ZJjfkx6BWdBqYuuZsEh0ikmAE/72q+i/s+Hp74v6/D/zwIU//+8APwpFV4N8GXBjJK/ABEVHg76nqi/vu/yPA9x/0xHDa8x6AZ599dqqLPFPU/ygVouVlc0psJMite+Y0mdkc1+IRPFdmBkH7X1k+gPn6lM1OrQa61EFU0djh+tbFGhVKvtKmf6NJup0T9TIkK5DCI4MM4ohstYVGQmPb0j6ytQuDIZrlJo8cjZBGwwaC9PcVtQ/6jj7M91YO2CzmCHoG8ZqICPDdwMdU9W9O3P5kyNcD/F7go4cc4v9S1Qci/H2v8Y+Pun/aJP92VX1dRK4DL4nIx1X1J8LC/jxWef7eg54YNoQXAV544YU6PJ4xRMvLSCM14lnftMJjSCPoRSlpLiBalTixgmunhS4vIKMcdsIG12qicYzb6dtn026BCKPrLZKtEY21kVkbDAuKhZRkYwDOUSy36d9IWP7wffKVNtH2wD7nUYa7dgVGI6SRosOR5d1hOu97Tsm9xBmpa94O/GHgF0Tkw+G2Pwd8vYi8DQuEXwb+2MFr0G84fp1HP2aqJK+qr4efd0Tk/cAXAz8hIu8Gvgr4CtWz2C/PDx8Yve/g4mt9uluhHM/nN7YqaaRLUxtQMcqml9I66m9wEYMv1NtQlKurkMTIVg+S2KL3zV1LX61vgbPh7qOnlxguJ0RDz+4zbZr3R8TrA4iE5OYm+Y0lIu/ZebZN625GsdgMr6MQRUgnsWNvqA1DHwwfrqh6ku/ynH/flbPxplHVn+TgVMsja+JF5B2q+tJxj5vaN19EOkEXioh0gHcCHxWRd2GF1q9W1bk6Zz8Jec2Ss+RpIXFSEXyxuWXNTWqui340wvd60zO6mhhwPTlJqbyUKpPzhGs0zXKg17foXT14jxRqE7DyHO20zRAMGC4npNs54kEF4s0BUhT4RszgTVfIOzHZ1Q4+huRe3xQ1mUd2w1lAkiB5eI04NmOzEyi+HsC+v8+kM+Zhj5lH6AkuF4zvPsmDphnJ3wDebykpYuB9qvojIvKrQANL3wD8R1X976e4jjNHaSM8GdGXtx2lnikJfi518oFkj5JEViqPs4wCZUxkOlEoK2ehVna9TqrhGNOGxAnRlRUbnOKcaePBfooDZ0RMmiKjDBa6DJ+/St4SXBER9wpa9wo0inA7QyIZopHgBgXDa006N0doEiFZQXTzLiwt4Nup2Ta/ftvOFnb7488ndo+konksFTgKeglsDUTkhw67C7hykmNMjeRV9VPAWw+4/bOm9ZrnjYOGg5RE/g73tUdG7ZOPmwuETtaj7j/w94fBhGKleu1wFlVZCHhFogkZ4TkRlWs0cavL1p3atxy67uza2pYW0U4DKRS/0MT1RvjFFoPrLTbenLD0ck7j/pDh1QbJphVbsxsLRL2MqJeRd1OKhiNdH5EvpTQ+vWGbyNoGbs38gnSUVQoboOp2nff0yjRwSQzKfivwDcB+KZpg6e9jUXe8TgnznJa5MOwz9tJi79CNcjZpRfhholHVsl8eIzz/rNcmUYS0mhaxO2dSx+HQGqAA+gNLozQt6vbLHYZXTSXTvuuJhh6fOJo3+2jiKBZS3Khg941dmvdGFK2YxkaGb0Skt7ZhtweZkbq0W/aaJSY2vnmJxM/7f+qSVAv/I9BT1X+3/w4R+eWTHKAm+QvC3ETwF4XJQqo4JLKfVeRapo4Oakw7S4IP5A6hu7TTNploWfQMt5MmlqKJIzR2SOHIF1JUIOvae0m2RqgTfDOcfThh89d1aWwUlsdXJdrJcP0hsh0IfhA6Z3d2kSQ5dvpYjZPhDL1rHm0dql95xH2/7STHqCdDXQBqgn9I7JszOplvLzuNJYn3kHiZr99TgHxUpc1EQRfGKZFyWIr2+mOZqKp5z0SRRdqbW8jWLvm1BUZLCUXToU7ovDqgaESIV4pmjA+/JzueZCuziVCDAlE1rTzgd3tjQlel2NysCf6soFj1+7jLBUBEvuo0j68j+XNGTfCPiD1ReIh4y7TEPoKz9n2HS8KwjtFovBnEyZl5tpRDt1Fvg1OCoZhANfWJyMFoBI2GSSa9Eu/aeqOhKW6ifk7etRSNJo7+9QYaCZo44p0MtzNA+iNT5WztWNdslptN80X0Hsx5nv+SpGsOwl/m8A7ZB1BH8kegJuRLCnFjWd++tMwk0fuJSVWTBUiXpg8vp5wgtsr8KxB85eZYBLOwZgPKM4hRZoOxGyl+qYPknng3IxrkpJuZEXvscMMCn0ZI5on7nua9zAi+nyGjAkYZur5px81yfH9wcdH7HBM8COqPv1zY4k6BOpKfwEGkPjmKb94RX7lCfv/+RS9jL/ZHi3JAcXEijTOZlpksvEIx1so7CTa7p28SKs8OJtdgP8JruPJUXm2sXxRZKFUUsNCt3CPdMIMhVpxNInwjRh24UYFPBFHHaDGiNQhrLAo0jZDNkYWYWQ6quGajMiSrcca4vJH8gd2xh6GO5AOOi9ofh6j+UhD8ZCMTPJCHrzChHpmUVB4UoZeySomss/ShCL562bFUU4vi4AY4MWUNGiJ7r+bAudS22UDlsZyzx0RSEXzRTkjXrQCb9Dw4iDZ6yCBDdkz7Lp02EoZ362g0fi8X0dU7r1ArvB53OU+IyJeIyB8EPkdE/isR+a9O8rz6W3EKPA5Ef2EQh2s0D7z9KBIrI3eXxOOmqMnUzUTHK5EzSeMjpBkkio4eECOuKrSWA7ZpWF7ebfXBgSbmKe+GGb6RULQT086noZCsaicDTkhvbdtmsL1jEfxCF/IibCJ+vNGEzuIaZ4hL1PIqIv8I+A7gtwC/MVyOtTiGOl0DnC95z2W36yPCNZpImhgB78cBxCVOxrJJ9USdNup9ZcZV5sknUzdaFHuj3odEedzjJoHpaGRDuSOL1hGBvEC8x3nQZmwdraMcHTryToLLrQCbLaYMriYs//x9k05uhT6YRgPWNuDaKtxdgzi+GLO3/RtvebY1d5vMxUsoJ/AC8HkP4/VVR/Kcf779ccjvnxTxlStImuD7A1OITDQ77cGE/8ykLt2121VqpIyyy/sQV3nlPJKSZtL/JqR+EGdnDwcVcPNgnyCCtEKTUuSsCzZN7L/Og3hPttKyIutOhhsWZl8wLFj52XsmlSxdJLsda35KEyP4cHZwLima/a+xn+Anb5sSLuR/xp/gcn74KDYM/NSoI/lT4ji7ghrHIBBGtLhQdYQ+MDxkX1S4n0jLNIwksenPi5C+QMPtYvnq8LhHUZ+UFgn2wnvrAw8ct7J2iJAy4Moym8daBL+ayFEstXG75j8T9TPyhRSXFUS7I+J7O4yeXLTofxiKrHF4/yLWTDUcIiLj1NE0xTUHRejlbYfdPg8odfIXDBH5l9hqFoBfEpGfwUr2AKjqVx93jJrkHwKT6ZaHJfyZS9ucwT+wxAmu00bimGJz6/B2+/1qmv2RZFWclaB9d4hzaEiL6Ch7cPTdw6x3kuDD+rUoTAoZIuwqHz6xLkkTJElsEwIYDI3oe310dZFoq4+msXW+ZgXitVLZaCMmvbNj6hmws4L2AmxuQZKgm9vIyjJshxmt4oApsfxhyqaDJkrNYdrmkujkv+NRD1CTfMDDku5jE9k/4j9uadzlNzYpNjdP/Lw9WnTGssjKwiBNjUyTBF3fqNI2OEGHj+DZIu4BF0vARvhVY/Qe9MORKLI1pYmdZYhYmmYwhG4HjWOTTLZi4o0+6hzJWg/JAlH3+mYb7JxtJmkCO7t2nCyz43qrL0izgT/FZ3lqHDUe8CSEP+u4BCR/kGfNaVGT/ATmzhnyEkDihGhpER2NyG/dOdU//9iDZnxdkhgNXaWSpuYVEzl0oW0FpuEojLwbHHXo46G+8sav1tNsmGImNDztOUsI5Cah0UpHGdJJxsXX5QbaSPCthHwhwWXehnK/tmlpnGEooDZSVASaCTIU65KdsC7AOWuGmqhPVCqiaRLs/mMed30ecAnSNWeBuvB6AF7yP3Cq6LzeFA5GtLBAdOMafnubYnv75E8MKZrKPTJEyFoUY4JvNiyC9x68Ihs70GnbSL3RGZt0RSH/H9JDpQpofyHYNRtjZU+7ZakWX0Ce45uppWdCDkBFkNyjcWSRfiOFOMZ3WxQrbTSK7L2NslB38JXLpCwuIEmM7/fH6aywKbk0PVEx9pE6fo/DnBC+6PGXWUAdydd4dOyLHuMrYZbBtVX01ZtjX5WT/vOXUfFEmsaPRmYlPNnYFMdGfN6DCP7V15E0PTNPGtQj6YR2P9gjlEVdV/rSqEcmNP6SJlYsbTWt8JokEAnZUpPhakKyU0AMWgh+sYmMciQMAcmutmi8sh66WjN7fyOb9UojNY38aBQKzYzfa/U3CNF9fvT7f+hN8KjPdZ4Ge6vAJRgachaoI/kzwFzm5E8jzdufl3VinaWvvGq544eIGstUzeRz1at1j0YRpKkRnfeQZfj1DSRNKba2z84/fcJlsvoZ0jTiZPw64qzYGlu0TxGMyvIc2i18p0neTdFYSLdyxCvxbk68M4JCbbi39xRLLRsE0h+gLRuCjve2SSwu2DGHo6q4DBCtLJmPz6PKRM8C89aQdQbNUCLyBhH5tyLyMRH5RRH5k+H2VRF5SUQ+EX6uHHOcrxKRnxORNRHZEpFtEdk6yduoSf4ITJu8L/Xm8LD/rOrJ794jv3PXcuOj0emjxuARU5qQlQXQsmuVohhb+vb6NglJFb+zc3YkEyJjH6Yt+f4AHQ5tbm35VstRh2XTVTk0pBMKp60m2dUORKFoOyyIdzKifo4ECwMA32qAc7idgZH6YAiv3h6vJU2N3Ksxgt42lTRFB0PECa7VMvVPaZJ2VjitDn+erBXOpuM1B/6Mqn4u8JuAbxSRzwO+FfgxVX0L8GPh+lH4W8C7gSuquqiqC6q6eJIFTDVdIyIvA9uYxitX1RdE5GuBvwR8LvDFqvrBaa5h2rjURH3BeFiyMTKPxlr3cKyqoNmKjNgjZ7a7Xs9++lFV8D14Lm3lIR90+ipi+fbIbISJI8hyktvbaBpBNyXZGFSaft9pEK/3rGGqH9JZg6ENAikK2yzAulwjBwSlTRgpSORs0lSaooXVJcp1nXk0fZpibq2u2XsI1ZvAzfD7toh8DHga+BrgS8PD3gv8OPBnjzjUZ4CPPkzH64lIXkT+mqr+2eNuOwRfpqr3Jq5/FPh9wN87+TIvDgcpbmpinx5KLfqkZLKy8s0z3MICOMGlKb7Xmw7BH4T9M2rLgmccV5tSFWkH2wVyc4qUvpJ4kF7f8vTbO0TbvfFzWk0YhManKEKi8FoukGWSQjEygvfB9XK3Z171wY9Hs9w+h7OOpC8BWV/I/9vJm6GuishkoPqiqr540ANF5DngC4CfBm6EDQBVvSki1495nW8B/pWI/Dv2NkP9zeMWeNJI/h08uMt85QG3HQtV/RiAyGwVNab1RZu5pqizwiE665KwJ4d0V7p4r+hohGu18GE49h754DliHMV7tLze7Rix5wVQmNoHYDhE+kMj8btrls7JCyP5OIbtoIPvtEEVvZpanj4Ps2JHWfgMJmyUGyl+e6daR2U3fAlI+VKs4QxwQvXMPVU91ihMRLrAPwf+lKpuPQT//VVsmHcTSE/zxCNJXkT+OPAngDeJyEcm7loA/v0Jjq/AB0REgb932A53yGu/B3gPwLPPPnvSp9WYBZSSvygyb/b9/ieB0F0SV8qc0jNGRxm+OGJ+63mtHyqlT+mdQ1FY/jxJoN+H4dCmQBUedvuoBknkcIgsLVqTU+lUGQaMZNcWiHZHaJySLzRI1vtIz9vxW01rlirllYDm+dikbJqb3bykYE6DM5JIikiCEfz3quq/CDffFpEnQxT/JHDnmMOsquo7H+b1jzu3ex/wu4EfCj/Lyxep6jec4PhvV9UvxKL+bxSR33bShanqi6r6gqq+cO3atZM+bebwWEbxhwz52A9fKnMmJkG5ZmNirun5mGPtR7VmEUgTpNU0ZU2SWLeqL+z3LIfhEN3t2eNDjUEajdDk5McyySTGd5vEW0PyhQZFtwGAb6fj7tcssxz9aGSOluXZzrSli5ONVid9/BzgLHTyYiH7dwMf25da+SGskEr4+YPHHOpfi8jZk7yqbqrqy6r69ar6CtDH9reuiBwbXqvq6+HnHeD9wBc/zCLnFY8lwZeYGPpRjfKDvRFjIAuXpuOUTenpckEEX66r8qUJM1ZJEyuc7uwaeReF5cr7g9CBaxG4a7dDiklCJ+wSdNvmWxOJdbs6waf23n0a4Rdb1jCV5Wingea5NZjt7Fb2ya7ZGG+G08jLn7b4Og84m0Hebwf+MPDlIvLhcPmdwLcD7xCRT2Dp8G8/5jjfCPyIiPRPK6E8aeH1dwN/E3gKO614I/Ax4D874jkdwIWKcgd4JzaAtkaNPUS0J6++3+xq4ndpNfE7u+NI/lFTCOU4v1OmfSSJxwNBCMXhnV2k1TTVD0wUYfcSru/3barTKDNLhuEQf22ZvGN2B43bfdzQNgHJC9zusDIzY3URXruz5zOohpJnOS6Jj/S4fyTMC3GfFGc0FERVfxIONab/ilMcZ+Fh13DSwuv/gmk8/7WqfoGIfBnw9cc85wbw/lBgiIH3qeqPiMjvBf6/wDXg/xSRD6vq73i45c8uHtso/oEZrQekAgKhmNlXYrrwZsOKrftVLg+7jChCGoFsY3cqhY5muVkpBH94zXMj+FFmks5JM7WiQJIEHQyrDUAHQ1ynHWyIu2jsyDsx0aCwaVFYh6wb5FWxVda3raibxPa6ZfNTNXdW906JumyYxZz+JbMtEJHfADzHBG9P5PgPxUlJPlPV+yLiRMSp6r8Vkb921BNU9VPAWw+4/f1Y6uaxxWNL8PsxORg7eK9MtulLGMqhgN/tPWgI9oivDePOWj0pCYU1V8sYDKvInHBWII2GpWm8+ezg/T7iD12zcQxxRHR/h9bOAG0kqHO4UY70hvbcYYb0LKdfNX1NzKitJlRNpL8uJWaN4AG5REsWke8BfgPwi4zHlShwZiS/ESRAPwF8r4jcwTq5apwSjy3B78sT7/d7l9j81yvjLbACZVFULfxn0mgz4YkjReigzU/Xjq9eEefRoHAhy6smKLe8BKMM6XbsfRS+irqrQmnkkPLzyAsrurYaRDsDfCNB7m+aT83mlkXyIviNrXGH7eRaimJvBXAGyfTS4nLtmb9JVT/vYZ54UpL/GmAA/GngDwFL1Pn1GqdBGXmW0ToTUXTp2+4VSROihQXTi/cHFuGXefNHKShO1gBC05CfmI/q0tS8b/ZPqTrgfVTTnwpvzUhFYd2ucWwplTQB73GtVjVb1sUx4pw1O3U7YUpUkFdmfVzsbPTfYGT3bW5bKigUbu0tyMFnHDWxnzkuocvkT4nI56nqL532iScieVXdnbj63tO+SA3DYxvFT2LCdKzMIYsT84n33tr08xwdDE1CGbpHT930NFmwDWPyfJaHjUUqHb6GjaXU4J/kNSbXL2BnA3GMNFLTs3sfGpkEIXTDgtkPq1qUniamme8PoNtBNnbQlQWY2GS0bzYIduZh78Ul8QMb1IWR/FGf1Rnl4C+0u/xy+cm/FyP6W1jHqwCqqr/huCce1wy1zcEnLeULnMggp0YNmOhiDSgjammYJtx12mN1ihNcqxlIU/d6shxQpH0Ak8XbMGikmvKUpkTt0I0aOYumwUj+OGI6QDNuG4eig5BHz3JwDokt9eS6HXuNvBiT+/aOed8Ph9X7lztraJZZrj5y9k/mpDquODdOEV0kDpvxOol5OLu4XJH892BSzF/glCPEjyT5R5Ht1NiLxz6Kn5RDwp5Cq4Yctg6GlqIYDIPh1oTp1sRxJouNR8E1mpbjD06NxLENGyntkFtNay4aDsk//erx66dcugJmN0xk6hqzQHaIj6oIX4sC1+1UNQUpO1WdQ6+vIhs7SNtklEb+OpZmlpOnSj+qYJ9AYR4+c0GilxyXLF3zaVX9oYd5Yj00pMa5o1TUmFVBOk6XgBFw+L2MsCeeWKVvxI3HAh7wAmFwRlZ5v0ujYRtJfwCthOHnPE3j46/Bbr/ywTnwUGEE4H4bhSqNNMrsTCHPoHCQYrn36oGCNMwSWNVD5hFV5O46tNsW0XdalrYpRwUOhrZxAHhffTYl8V+4iuZx2GD0cqlrgI+LyPuAf8leg7IzU9fUeAQ89lH8BCROkFBQ9FmOS9MgPwzEGGa3MhqhIVVTqWFCbv1IoptUz5Rj+ZYW7ewAkNVltNOk8UufQXt984nPM1yjuccrvlyr63YsFVNaEwRYTt82j2Jnt8qVGxGbrh9VuxR+HKE30jBMJKSL4sg6ZSd/D344VeoqpGr86ORa/hpngMsVybcwcp+0NjhTCWWNGo8G9UjSMDlgf7DXZTKOx4Zbk3rzOK7mtbpmw5qQxJn+/KBoUpxZEA8H1bQkt7RoaZqVJn5lEXd3Hdncpti0jvBSkjhJ8BInNq91oVupW0wbHz3wer7XMxfIkPOvZsFmmalrRpl1v5aulSJG8KMRLHYtcu+ZZFTXNuy4SVz55Zd5eA3DV8TJnt6CGlPEJSJ5Vf1vHva58+EkVOPyQ9zYDre6KfiuT85BKArTm7da9pgktgKsGLn54NVyGMo2/+jKiqVoFrqwtGDR9Cuvo9s7FJsHa84Rh2u3iVaXzSUyrKecJ+uajWpYt6RJdUaxRxKKFZDL9ZPnVU5e1UO3bcSexOYhXxqURTbQW5oNK9yCpavCOMFSDVQ6dNaYPi7DIO/gxvtIj6kj+RpTw6SaRqLIRgGWhJdnFm07a/bRPLeURkhxqPoq+i3tDHx2cP+dWRYrUbdjuWwR6HZwiwvQH+A//Rq+TH0cskG4NDUiX1oMrpCJ5cvLCFy9kW2WIWlSpXmq1y/VP4V150oY8q3e41aXbU2tJgyzsaKn/HhUYTSyM5fSIqEozDe/15voAk7QYnTg+i8Us2hZMDv4VhG5d8T9AvxJ4FAb95rkp4zHNh9f2gakqaUavFZkUDU/gUXucTD8mjT+KjxKUXW/6kRufs9rBJWOazVtZF5/YHr0nV2016fY3t7z2AOXWkbq3Y7JG6NwghtHSGw6d8mC3DNNKNbW90TTWhS2SYShJwJG1o2GkX3Z+BQmQ/mlNu5+bqRfeMh7Nqh7MESSBLIMSVPytfXKQmFPHeKyEeplW89Z4XKcMP07zN79KLx01J01ydc4W+x3lwx59LJz1R4SmkzKZiXnjJhbLStKDofjLs/QCFQ2M+3XykcLC9UYPBa7pnRZ30SLAt/rVYqc8kzioPVa+iUdyxbz3LpZOy10sYPkHgYjdGsL3RmEM4fx+yrXUub1yzMUHWUQhVROp21nLN0WbrOHLnbQJML1RjYRygdveRH85naVfz9Tv54aJ8clUdc8Si6+RE3yNc4MpS/8ZK67+j2QafXYUKCUJLE8dBQZyWVZNZQaF8y9ghSyem5wkPT9AbK4YEXMOIZ765VfjB+Nqujaj0bjNMe+iN41G3b2ECSSZbqEhS660ELTGOntWOoGbG0TEs2K4EMqqdT+i5ou33Xa9pzBEBppNbRb0xjJwnuKnPnQt5r4u/cf+AwfV1z4LOXLEck/MurCa40zwUFacrvDFC8uTceDNqBKZUgjNYJPU3Rr2wi63zfTrt0exfb2g81QaYrvD6wgC1X0rUVhx0wTXBLbRjAc7I2AJ84ESrWLhoHbWhS4dts6UNUjvSFuY9c07I3UFC9ldB3sGSY3LtdojhuVosgaocLGQZrgr63YsUSQwQgZ5Whqw0Bot9C1DRtveNAZR41zhXA5Cq9ngTqSr/FoCJHsQYM3JIosss0zS29gBU73xDUYZWiQDqKKrm+AqkkS4xjd7e85C6iOF4i51NcXt25bqiVNkaUFdGcXHRUWWR9CltVIwdIgDFCCRXCeI9dWTfmiCrs96HbQO/cqe4X9Bmsl/HBQrVHC+nQ0wl2/Zvn9UU7xxCq+GRGv9/CNhHypQXp3Hb+5ZQqgNIERaH4JcgWPO2aExI9DTfJTxGNRdD3KoAos8p0oUkoSQ26j8awDtW/5b7H0Rtm0VA3ECA1QYCkRcYIOs7EPTZyYiqbbgY0tfJBHHrYmiSLzyMlswAdgUffqiv3eSKHXR9c3reDaaaP31tDhcK+6R/d13E6eITSsH0CiCJ5/Bt9IcJs9xHui7T7RekF+fZH43g7pnTW0Pwia+pGZs9WpmovHJYnUReQ+8B+B/wD8e+BnVLV39LP2oib5KeIl/wOPB9EfhDD1SUtv+Dip/GNIEvTGCnJ/C2m10N0efme3SueURVeXxLiVZfCeYn3TcuyANMy217Wa0OlAluE/89rRaQ6xx0vTBnrg1XxjwKLnbtvy5HfW0OGoKub69Y0H7RXC8fZscJM1hyRGFhfJnlkh3gzSTRG0kSCbO/gri8R3LHInNEuV8tJTu23WmB4ux5/geWwq35cAfw74IhH5FIH0VfWfHneAOic/RTy2BF9igqi0KMaRsC9wt9bQ7R10a9tSGt1ONZRa4gRpNHBXVu3xIlbAjBxEDre6Yk1Oy0v4O3fJX3t9L8Hv85137baliZYWIcuRNMV1O2NSXuiSr3aQ7T6o2sCPXq+axFSuvzyr2J+mGctFE6KlReTKKvlTK6gTim6D0WqLYrljfjVxjLuzYQXm3V00DBSpOnvDRvZI3vk1zgRnlZMXke8RkTsi8tGJ2/6SiLy2b7j3A1DVLVX9gKr+JVV9J/AsZjv8u4DvO8nrTzWSF5GXgW2gAHJVfUFEVoHvx2YVvgz8flVdn+Y6LgoXHcmX+fBLAfVGYEVRebVIbBbAAH5n14qeTUt1AGHakrdi5M4Ort1Gh0MYjSjur6E3bx36WvYCbkzu7Ra+20J2dq1bNcuRpQVklOG7LeJfedWcIINvjsRxlaIpc/h2aHOgrBBuL18ne/MTFM0YHwuuUMQrycaQ6N6WHT9YDLPTs9RMeaaTxHamM+HSWeOCcXbpmn8A/G/AP9x3+3eq6ncc9UQReQqL4r8E+I3h5g8B/yPwUyd58fNI13yZqk52bH0r8GOq+u0i8q3h+p89h3U8drg0BA9UVsOJ+dFII7WRdnmGa7chpTLmqtIjZct/aCzSUYZ6PdI1snq5sgDaaUG3Q36lCwLuLc/ik4jo9qbVC1ZbuFdvA6DbO8jiAn5n12oGwYHS7vSHRteu1cR1OxTPXkcKZXAlpn1rhMZCcnOT/Poi2mkg2z3bRLBmKS2ln4V505RzbTmks7fGOWLsdP3oh1L9CRF57iGf/irws8B3At+qqqeWXl3EOeHXMJ4u9V7g91zAGmqcAyZJ0oUctwSvGs2yynFS89zseYPFAU7MV344qpwYzdXRBlX4Xu9Y/xYNee5ibR29v8bG57TZeWOL3hu61tzUTM0/5rU7tmnEtk5/304qSxklMOGxM55qVUXwrSZuZZneFz3H7jNtdp9pkex68lZEereHX2wRbfUtVZPldmaytlF55uvIbBKk0ahGINa4HDhhuuaqiHxw4nKs18wEvklEPhLSOSuHPObtwPuA34tNhvrnIvLNIvJ2EWmc5EWmTfIKfEBEPjTx5m+o6k2A8PP6QU8UkfeUH9zdu3envMwaZ4rJASHB9AvMYVKHQ1OeSIjQQ2OTOVGqTU8SsWHY7VZ1v/YHe3X4p0hniDjadwv+3P/yXlQw7fvWDv7XPl354uhu8IgpfWWKsht3Qh0UIu1q42o0TR650CHZzhisOIqGIIUS93Ky5SaSe2Rjp0o9+fUN89GZmFxlTpW+8qw5U8xgbv8l/wMX3wgF42j+qAvcU9UXJi6Hesjsw98B3gy8DbgJ/I0Dl6D6U6r6N1X1v1TVLwL+DGY5/F5g8yQvNO10zdtV9XURuQ68JCIfP+kTw4f1IsALL7xQhzezBPW4RtMi8tCoBBAtLZpMcDRCWq2xt0vkzH53dzekL4zoSh29DoeVh83DoHj+KV75SuGbv/e/4flP3MO/erPqiBXnq6lU5Rg+KwCHsYNeLYVSEn346ZoN3I3rZM+s0L/RZONNjvYdpble0Lg7INoZmNbeh43DexsAEkYCEpvZmEviSlGjZRrqLPPxM5bbvxTkHjBNWwNVvV29jsjfB3740HWIfA7jvPzbgRUsH/93T/JaUyV5VX09/LwjIu8Hvhi4LSJPqupNEXkSuDPNNdS4AJSDtwtTnJSe6L4f1CtxDJ2WFUB3dpErq5Y6KX3lCw+RyS91OHyk6NZ9wefRf7LF0rOb3PhHKf6Xf23P2YANKMn2TnMKBF+lZ0SoRPFxjEvNvmH07Cq7zzRZ/2wh2YL27ZzmzV3cvQ1Y6BjBF94GoAyG4Mx/3rXbtvGpVGcnEh3cUPa44DKRO3CmOfmDUHJguPp7gY8e8rh7WKT/H4D/G/h2Vf3V07zW1EheRDqAU9Xt8Ps7gb8M/BDwbuDbw88fnNYaalwQ1AMTMsOglql8auIYf2Wx8m+R7YGZkvUH5iET5IR7hnc/JNytNZKlJk/98S3y12+NTdGiyHTyzvxzyIvKHsGeKCDB+rhM3QQPeYKhWv96gzu/a0D751pc/WhG65P30SQ29UyvD7v98HhBmg381k5lQbx/+MlZE/ylUlYdgktH7BOQcDmTY4l8H/ClWP7+VeAvAl8qIm/DtpKXgT92yNPfrKoPpGVEJAL+gKp+73GvP81I/gbwfrF/8Bh4n6r+iIj8J+CfisgfBT4NPOZi8vlEad4lui8ccg5d7iKDHLe+awqaRmqySidoPzt2MMip1tHrE/+Hj5JP6OgrU7KJx5jro51xSByPxweqGhmHrlu3tAijjM3f8hwbn+V44v/XYOk/vGwpmSxHiqKyDS498n2/b121i12TjBYFRe9UTYunf98nIPiLlPheZoKvcHbqmq8/4ObvPunTReTbgKexAPkl4JuAbwY+DFwcyavqp4C3HnD7feArpvW6lw0XrZW/SEgUVTp4qVIxBTLILNINUbPuWuPRfrfJs0CxOQ6CyrF+ZRHUzhqM3M2/niptI0HTb6ZniUX0UcTdr3yee79txOq/d7zxX9yBoJSRlSVoNND7azAcVSP/dDi0ouqzT8Pmtg0x6fcvvKu1JNny57S+ozNB5ofgMtgaAP8IWMdy8P8t8P/GBMdfo6ofPskBaluDGlNDFQ0XBSpixNkMqq+iQAdDS81MgdwfQGk7EEVhduxo3M3qvaVUqmHiprFXTPqpeQ6RY+Ndn8POG4TP/Yt3YKcXCqqKPHEdvb9mZJ8kVliOzUVTogh346oRfJrYvNaQ758nzDKZH4rLQfJvUtVfDyAi3wXcA55V1e2THqAm+XnCZfM8UY9EadC3K7LYhf4A3drBb2+fG9m5NA0qnmAtXKp3imJcL4hja7wqjcWcIM5UNdLtoKuLLH1sk5UPjfDLC7idkG5ZWULv3g8WCN6OAdWAb3d1bM2g99eh8Lgk3jsAZYYwl2R+EC7J0BCgyrupaiEiv3Yagoea5M8F55ayuYyk4QRJW9bhevc+iNjEpvN6+VZrPDvW7/XSqUbzldp32VdqK5+XZcjuABkMKd5wDbczgGurFItNolduj713hsMqPeXabRs+kiTQMwM23x+MfW8u49+qxl5cjkj+rSKyFX4XoBWuC6CqunjcAWqSrzFVaJbbZWt7bA1wHmccwbemKp5mlnIpJZySprYBNcyLHifgIrTbNPuBoLgxGWRhmvkrS/aeGgnZQoP09U1LPz11Fe5sUKyt27CSdtuOu7xoUszgLx+Fxic/vHyql5MEIo9NFB9wGXLyqhod/6ijMXvtcDUuD07QTal5Np6WBOP5rlOGK22LwyYDjD1x1IaCS5IYSTtBFzvkTy4juccvL5jNQTPMfb22gi62KToNKJThtTbp7W3Y7VE8sYLbskYtl6Y2lDyJYaFrBD8aVWcQRvCDc3n/D4OjSPyg++ZeUHCyjtdLjzqSPyfMjcqmJPbTpBwmh3ufUx7ej0bVzNnKb8bJOIqPXJWDJ00pFltEO0N8O0WDlFKdg8YyrjekWG0S9UZkqy0at3Zhcxu9umxzWje2zIcniZGlxWDLECSbWT7O019igi9xXLS+/zs8eX3eIv3LEMmfBepIvsbJMUnwp31OeJ5L4nHKZsrY30xV2SiAFWGTGG23GLxxBVTxrQR1DvEe30gQ73E7ffIrXVzPZJHJ/Z51tC4tILfv4z5zOzhKZsjTT1QpIZNRhqEgIjb9KQwWv2g/mWmR8VwEMSUUGxpy3GUGUJP8OWKWI52HKRjuj6T3zFY9wrr3zBAGkFjULpXuncgieF3skN3okmwMcaMCGeRIVqBJhOsPkawwr/lbG+TLLYvQvbcUz/qm+e3Esalv3vQ0DDNwUUgBORv+/cQVdDAwy+NgeDarRdeTkPi8EP08DfKuSb7GiVCNpTsJxFXj/lyaIklsTUjB0+ZhnCQfCmWBd6LYKk9cg4UuxVNX0Sgi3h6ZtUJvhKji24l5zZeODHkOjYT43i6aRPDanXFuf3fXGrq6bWSUWw4/vF7+xhvmiXN33czYosh8fKKHqKOd8WY4L0Q8dcxJTr4m+XPGzEbzpyEa9VVbvRbFWGFTFHuKsOcB9aHIeu0KXFu1iUxpghtmEDvcfZMcS17g05h4bdcKsNu7lrIZFWjskK0d5JWbyPKSzZTd3oFOh+L5J/Ht1BQ8SQStJtlTy8S3NtDl7nghe7T6p/y3m9HIf9YhqsdeZgF14fUCMHNF2IfJxcPRg7XPARJFNsu11UQ3NpGVZWg20FaC643MXqGRIHmOtlKitW38QhvXsw3K7Q7QJEbWtkxL32rB9o7p4peXKK4uos5RLMRAg2RjQLHUIvnkLVjsIuvb6LUVpD9CNzbR3b4VhKPITC0vkLznuWB6JpihSP041JF8jeMxBTLaM1pvitDRCO31kcVF88tRxd1eh+1dyHK0kSDDHNncwa900WaM2x1Y3r4/QO6vj+fShk1LnrrB9hc/S76UgoAmjqiXI31L/finrppHzY0VpFCzUc5ypNUMBH+5LIXf4b72REHHSTaDedow5iUnX0fyF4SZi+bPCBInZvhV+GDRPl3CkzRFOm1raspz2NqGRgO6bcgLi7IbMYM3X6H5+jbya7fh+lUrrEbO0ivOheJtRHF9kcH1FgBu5CmaMZJ5S/+oIqMc2dqFdssUOVluTpuq5mkTRZfWAvgd7msfiaT3P/eg7/csbQKXxNbgkVGT/AVipoi+8oZ/uG5ViSJcq7XHO2aq6YrS/z3MiC0JVtqtYG08gmZKvtzGZQWtT9w1uWOSoK/dQrqdsbd8tw1bO9Bt4NOIZDsnGuT4xJGs9XHbPUgT873Z2oXdvm0kW9umsNnaMf385hZaXG7mOI7oJ+8rv7unaZQ66jmXDjMSqR+HmuQvGJea6CfmtAKHk/IxxO8aTRsWUhSVf4vm2bHDuB8aYb2+HAJeFLhW07xkuh00TSCyoeGun+F2+rb5bG6Zxj2ObRRfcKzk3jq9L34TPhHivifZGkGhJBs74NXy+Otb0G7ZJvFEG9np45+4inv1tkXxm1s2LGVnZyYGepwEhxH1pf0+nwYzlI45DnVO/hLg0kU1R/nL7L9tQikicbJXOSIO17Dh1DoY4gfDKg9/FlOfDsWEnl+c4DptS7kUHnp9xHvTvYOlWXZ20dt3x+qXIPM0L/gCfeoazdd3GS5FZJ2IbDFFVNHYUay0rTmq0bCzgzhCCkUXO/CJl9E8x4ch4drrWUPUYT0Cl2To9kE5+vK2k+bvT/Ialx61hLLGWeKRif4wgngU4pgk+0lCLm8Px5YoqhqfJq+X/jHl2DwbvuEfbO+fYheoxAmaZXYZjaoianJ7G9ncMXuCLKvcMTXLbKJTMBUrnSg3fv0S228UGusZLvO4rT6axkT3d9CVBXCCXzW/Gt9I0F95GbeyTLG1jSQ2HUrDYPBDN7hLJpU8itRngqQfAXUzVI2p4JGI/qAIuyTO05Lo5LGOIR4z+jJTLokiom4Ht7RItLSIW1hAmg2b7ZrE+F7vYFnllKx3XWlhUNhwcB1l+J1d6A/sAjb6b2JUXznEW/McaaTolSV6zy3hcuWJnx4xWoqJ+uZo6e5ukl+3WbXqHL4R2+237uGeuIZf3wCwM5iQorlsyppHwdwTvddjL7OAmuQvGY4l+uPIepLQH4Y8D4ve9z8siqpuVjD/dHfjKnL9Kjx9AzotSBObwtTv28i7c4YfDqpLGT3bpCqP9nro+sR8ZBHL4Zf/uF6h2yFbaYHAwq/u0ri1Q/fnb0KhFEtWwHUjs0Egdqavv78OrZZJKHs9O5spm8Mu21CXC8SlS1Hux0lSNSfkeBH5HhG5IyIfnbhtVUReEpFPhJ8rZ/wOKkyd5EUkEpGfE5EfDtffKiI/JSK/ICL/UkSONb1/3HDkP8ARxc/KK6b0RzmIsI/bJI4j95CGkTCY2i0u4K6swvIi+dNXGD4XJiElCX5zi/z+/QtviqoQNOq+30fE5JGahQJoZUNs7186LfxCm/SVe3R+ZY3o3hZ85hYkCdHadjAs84xWmkRru2ZWdn8TOh10fYNibQPXaO5973NI8LMukzwK4o+/nBD/AHjXvtu+FfgxVX0L8GPh+lRwHpH8nwQ+NnH9u4BvDXML348Npq2xDy/5Hxj/sxyXbpmI3EsTsEMffxDRnCSVE4zFxAnSauGWFqDVZPs3P48udcivL4JC4xc+jbx+F//6rXOdAHUUynVX19PUJI95buSe5+ORgF7NrXJ5CTcYWQPVYAjDIfrGp+z3PEd6Q0himq+sg/f4nV10YxN/83ZVaJ0Fa+Fp4TCi3/O9vuw4o0heVX8CWNt389cA7w2/vxf4PY+63MMwVQmliDwD/C7grwL/r3DzZwM/EX5/CfhR4H+a5jpmGUdKLCeVLc4GZetodHobgqMeV5J76STZbFS5a11o0369h+wOiddsQpmOrMCpmRUvK5nkRUax4lCvuEb4uhcFurNruXInRuzhZ7SyhN/tWdrFe2RlyYZ2N1J8O8EN2+ide8HDJodej/zOXTuLSpPHjthPO2hklnDCwupVEfngxPUXVfXFEzzvhqreBFDVmyJy/SGWeCJMWyf/t4BvARYmbvso8NXADwJfC7zhoCeKyHuA9wA8++yzU13kZUf5z7KH7CcI3iWxRaKBUM+suCfO3CNhXJQkNBa1msjukGi7Z6+bJlUxs/JsB1wUUezsns16Hgalva8483xvNdHhWM45qdWXJMZv7di81tHIJJQ7u0i3g19ZJL65gXYaSLttKpqbt/HDQWVVoL0z1L5P5u/rXP75Q6m+78fgnqq+MOXVPBKmlq4Rka8C7qjqh/bd9UeAbxSRD2Hkf2DCVlVfVNUXVPWFa9euTWuZM4XqVDcQV6U5n5DmnUmDUcjvu+Cc6BYXTD8eujopPWcy6xCtBmWE4dWSJmYMFmaanitCqqqUcVbpKyf44YBic6s6y5BQFJUkrvTxOMFvb9u6xaGDIdmzV3F31sAXyMYOjEb4O3dNFhon42lX+3oEHmr5cYJrt/ce45IS/KxH6sfhDHPyB+G2iDwJEH7eOYs1H4Rp5uTfDny1iLwM/BPgy0XkH6vqx1X1nar6RcD3AZ+c4hrmEi8V388Hsn9iRD+Ra9aiGBddH5YYxBF12kbSnZbJCIcjy197RRphOHaSGLEXYYiG10D6Rpg6ysZ68HMiqaoo7ARpNKrb/XBg6p59RWhJU9yS1f01y/FZjh8MLX3TaeGefRo+903E6z1rptrto7s9/Oa25eUn8/xl0TtO9vQMHL1gNy6Ylx78YfB4Wdy+zJhnCeU56OR/CHh3+P3dWGZjKpgayavqt6nqM6r6HPAHgH+jqt9Q5p5ExAH/I/B3p7WGecePDr6XD2T/ZA+5P0qqxjWaREuLSLtlAy7KFI0q0mmHyDc084ysu5PhMExEMpsA+oNxVJzn01fWTBClpIEsGw0kzFXdM+wkRNvlYzXLgoa+qD5D12oSPfMUcmWVYrWDTyN49Tb+3po1NIX3M3nckoyjsGGcZFOb3AjK15Y03FYU1edWnYnUOF+E7/2xlxNARL4P+Cngs0XkVRH5o8C3A+8QkU8A7wjXp4KL8K75ehH5xvD7vwD+jwtYw1zhpeL7q99PHV0FAnHNhpFjp10Zc2nPtO2SplWTEJEY0Y8yuz2JoW8WvoiAc5UHzFQJvkxZNRoPpIV0OLSNpiTbfaSrRWGbEFYoliiy999q2ii/TsO6WT99B+0P0NEIPxpVtQ8pZ7WWqR6v4IRicyukuJoU29uHL30iQndpWn1OlddOeb96K6ZfQI/BSfGozpWXGWfV0aqqX3/IXV9xNq9wNM6F5FX1x4EfD7//beBvn8frPo44id1rhZIo09TIst0yglfF7+xaNBzbV0TLCH00ssi+3bJIPk1trmmJ4RBpNcnv3p/G29vbCyAOca4yPau07uLwWX5wF7B6xNkwcc2zagas67TtvcQx3LyLiyL81o699zybkKYGy+ByHWmK39kZrwkOr0NMyFDLtfjRyM7AvBsrmJzgut1K03/ZMbdEPxsNrceidqGccxz1z/eO6Ov2pAJ0MDBZoKpZEXQ76NY2qlqlYKTdgpUlfMNG6GkSQTPBt1PUCfEnXye/c3cq78U1mhY5i6DDodkEhE5SLQokiasxg4d5w1SKGs0qVYxrNa3G0G7hX7uJRJHp3kMhe1znkNAn0EUHQyPp3d44bRNefzJXX2JPfj2sd/J6VQiWUFCPItxCF7e8RH7z1ll/lDVOgFnxpjkONck/pnhH9HVGWCVplhF6o4EsLUBeVDlHuXYFsozR89dIP73G6HqXbDEh3RhVRclos4//9GvkU2iAcu22RceRQzPzjSlVRDbD1R7nA/EeZdamRYFEZdNTbAW2KyswGOI//dqBKSb1Ok7VNBpmYlYUoX/Aegd0NKqKu5r7ByWQpZQTILLfNctxaVTdRuFBC7TwaDFA+rOjuZ+7cYIKFPPB8jXJP4Z4R/R1pn8PUaOWyplmwwi9KKDwFE9eIdrusPuWVaKhR3KluLZItpjQ+fg9m57UbMBgaNLEh/VIP0gHPtmEVXgrlJZpkFAoLVH5wlQ3HO25Y8fG0lEi+M+8NpZBTsAlsW0G4bXL3L8ksWnph0P7ORja4xtNgCoFA/vSSyH1JUnw9FEPWW65+PAeygKsay3Y2cRwePrP84IxU4NBjkAdydeYWUSd9rhIOhiCE9xi124rO2abDaKdAdpIaL+yhdzbQHd2ketXaX/iM9ZMlKbo2rqlSB6G4A/x1ansE7odKCwlooe5V4bnVWP1jmgcMvvgCVWL92O1UKnKKXyV+pE4rgaduMTy+K5l0lGJIjNj6/UhiXGry/i798dWDmFNqDc5Z+ERZ6+hRYHmebU5VMVb9bjlJRtssrFpMtQZcTqcS5xQPXPZUZP8Y4Z3Lf0RI/MkhhCVWgoiR66uwvomxWc9DYXiPvFp8GbRq1iEKv1BlXYoHjU1MxntTvjuuG7H7h6NzOP9MGuEEPmWeXNJU2sKyw/256kKm2k6Phvwaq/vJMydncjnRxGiims1QyRvyhy3tIiuLiJrW7C6DLfuWvdruRGVqZmyBjAcWoNT+F29x7VaSDe2syiRajMp1jbQu/etTjDjmPWCbB3J15g5/I7uu3ELXXQ4MhfGdiv4rI+C5FBheRH3i59CvTcCgqqzs8xL+yzk6gMpP0qapkxPlE1A0miY57v6I8m9fP0yNw4YKVOgB0TzJm1MrHAMEMdIllUe8xIl4ECIIEoqO2Jpt/FBGumWFmB1meGTizQ+dcdqGbfv2ZnGvs+gPDtQnxOvrli+vigsmi/Xvdi1M4et7er9SKtpm+3WNjQb5J9+9eE+2xqPhhma/HQcapJ/jOBarbG1buTQ4Wh8Pc/RV2+OtdnD4Vhe2GgY4Reho7MsEkauav8/bVerRJGlMYITpCwtVpp0YDxB6aDnlgXMkOrRAiRtIs6hgHhFy6eGaNrduGpSzzy34dq37yKtJr7UsxcmZyyjdtds2EbXH+BaTdz1awyfv4rLCtJ7u/hbd/fk3vegPEOJonC8UegrSOyMqdNCVpbtbGowgrK71TlkdRk2t2yj29g0u+LHzPTsMkAAqQuvNWYJX/nEnzBVB66akqRZNibzCbLSLK9sAcrCrMTxWDM/HNrvxcHNRsdi0vu+0bC8eH9gl2MGfI8lkPaa0YJ532lR2JASH4eIOQ7HCVH+bt90/Y0Uvb8OYEQaziZKwva9nm2GoTbhuh30iSv03rBENPTEv/QKuts/kngr35xQsJ0825BWE1kMnbF9k6zKtSuwtYNu7VB86uU9n9Nl9a05CWY5VQMgdU6+xkwiiqoIXoeheBoMu8BIqOpeDT9L1bfmedU4VezsjodSnxQTjpDl3FfzuRlBfzA2WFP/IMHJhE0B7PF5IU2QUAzV0aDagCQJVguFkTciyCAoW4rCJIzNhqVSyk2j260siOWznmX4xALx9oj2z71iHj7lhnDI+4ufuI7u7FoaqBwxWK45csEMzaGdJvlig+TeLmxu4ze3Htw4ZpjgZx51uqbGLOErb/zxyg5Y+wPzSU/TkGopm29cVZCk9H3JMlyrBUlcSQgpCohjom4HH4j5RChJW5JxagLGFgr7558epnefkFaWShkRh4rpzquu3cKjWZA9lptVOa4wjq0zt8yNT5iLaUjZyPWrDJ5aJN7NiV67H3xuCorNg22T4zc/z+7nXiP+yOtmZZzlYblSDTR3nTY0Gvbe728Sf7oHraZtCjMolTwKsx7Fw8m9aS47apJ/HCCCjoZjtQdmU1BqvaXUzCexbQRh8EdZlJU0tU2gYRuD9noPdm0eh9IcLDUJYtkxWuXeDymuApUUUb2O58oGhZCURdhWE1pNGGWB5IM3TVFYqimJzQc+jmF7x5qqnKla3PVroTegwG9uI6vLjJ5eJr3fx91eh8HQUkmhm3X/+46vXQWg/e8+ht+fnokiUwt1TV3DcFR11Uq3Y7WQ0bj7dh4w+wRvqNU1NWYCX3njj1v0niZIkqC9ftXQI1GELC6MbYPzHPIC7Q8snQJBq15Ab4iGpiRZWqRY3zx1qsal6bj5Cg73v5+QO1aEGSU4Z3NZpRGGkqgaaU9aBkSRpWO8IiEdJFEE3Ta6WRZZi8o1060sjwl+YzP45qe4XHGbPcgy/PZ2FZkfpOnXUUbxqVeIV1fQPA8FbuvMlcTOEnTDJmdJHOOurqJrG3acRgNpmZUEl2UWbg1DHcnXmAl4rYjObwcjrVbTiK9RukgOxhrxPK+ibUlTKKP6sjFot2cqlpPKJitTMKkKm2ZFsHcq057IPahnxjNlrU4gofBb/fMlia2/fH6Wja0YWk1LOwU1DaHxiDy3XHy3g6qldMhyNDc/+ejXPQ9bu8Sf2MRvbgXdfVatB0msUBsM0DT3FJubxNevmSpnZRm/sWmbYVDq6Ghkm2XZnxDHyJM3oNe3syLMGoHN0/5xa0wNWqtraswIymabaqxdHvLWIrDQgbvmk152Y/qQntHRCBdFsLxo7f87uxbxPoyfSpixqsW4Weig1ESVZw8FXvU+6NcFEPN3KaP2snM1itBO0wZre2eeO82G9QD0+kBs0slmw2oAUWRyzd1e1YCkeUZ8/ZpF9evb+HtroSs1q3xnXEhZ6XCIHz24dr+9Y01lg4Glf3wBYcQgV1dhc8s+87icqpXj1zeM3J2Q37lXfTb2h6uLrheO+eD4muTnGe9a+W/HrfPNhhGdT2A4MsJ69aa19hfeot5GinNmdyCdtkXGhYcsQ3t9/M7OqfPGB+WwgQcsCfZ0o5YujKUpWDzxNY0iO6soC7POIdsTdryLXXPSzHL81SXc2o49Ni9gcYH8Spf4ldvBOtkRrS7b8fKcoiTasL5JjfqxWnUfhqu0O9BMYXvX1tpswK4pe7TXR8O8W9/rWRds4Sm2tueK1Ge907VELaGscemho5G14Ydh3AyGlSmW5nll/uWevGEkk6YWqfcHlpsfjdCNzUca/nEowU94u5SDsyuyjyKLeJ35xZOHHHpSRvXY7zJh6dtIYAgaCYJJKMucOmW6aTgkfn1o9YlOC7l2BU0iZGvX8vHYdKdic8si+VNsaNJqmptlrw9318JwkaQqAPvN7Wrtvt+v5rj63d4Def55IPy5MCmbE5Kv54rNM6rct0cnCb6MjiUYk/nCbr97H4ZD/M4uxeYWvtebynQniaJKJVPKFiVJxsPCozCspNTqp4kpZ5wDF6Gthj0GwHuLnIcZ2mogo8KINsvHefjVZXtuq0n+zFX8sjVQsb6J7AwgSXCdNtGVFRvifVqbBnHmGHn7rq0nssHhOhxVtgmot2JsURA//ZSdWeX5hIY+1B8u+VzX02Jm58Aq4E9wmQHUkfyc4p3JHxg3H+EsPTMY4pYWkVEGSwuWYtjewd9bMz18nlfa94f2ozkGpR1vqcsvJZrSalZRL0GRgnPj39O0ymeLKr6dIqMcJEZyI3oZZUbwjYaN8Gsamcsoxz9zlc23dFj8ZA+3bgXV4v4aURTBKChoHmVDK90mnRvbMxRhPCAE+ahtTP7uPfZPrxqfNcyHjHISs5i+EbRO15wUIhIBHwReU9WvEpG3YcO7m0AO/AlV/Zlpr+NxQzlMQ73CKDOtdlGMI81hILQ4MrItCmRpEXZ7Vd74zFDOkS0JTz3gzH4gzJYliY3Uq0g4GitnyiHheYG20ipNM7q+QHprGxVBhqNxw9dim9GVNkXTMVyJad7L6F9L6L46JF7v4a8s4uIY3dwiv3Xn0dIjoShbjeoLE6X2oBwLOBwebro255hFosefzd9IRF4GtrEdPFfVF87kwCfEeaRr/iTwsYnr/x/gf1bVtwF/IVyvcdZQb9G4+srlUW5cG6tPSjL1HlleskKkCMXGxtlH8WW0GhqfpN0OSpoUaQVXSBeN1TOheYnIoXGExqEDtxFXBF90Uho3txDv7VvsFV1ooe0W+WKDu1/QYOu5hK1nBZ86Fn+tT9GIKJZauLVt8Eq0tHgmBE/kxtbAlWlaUf3ULJ/uUPMZwMwR/Nmna75MVd923gQPUyZ5EXkG+F3Ad03crEBwaGIJeH2aa3isIQ63sIBbXYFrq1ZQTRNIEvTqsqVCFrqwvYNbWsBPafh2SYbR0iLR6rIVJJ0bR+hJbAXSsikLrDbQTu0bKmIzZOMY6Q3RSIg3B2hkt4+eWMCvdtl58xK7b1nm9d/aJu6BCjTXIN0YkS0kNO72iG6tQ5qiG5umIHqI91KN8ysln8HwDVXbyIL/zmTxVpxU9gmPWxQ/qxDVYy+zgGmna/4W8C3AwsRtfwr4URH5Duxf+EsOeqKIvAd4D8Czzz471UXOI1yrhcSxEVlRwNqGKWpWV9BWivRH1mJ/567pu4fDqUSbZQ5eJnTtZkcQiqpRUM94b1r8JLGCaZri7m+RP7lK1Bsho5xioYnkeRVBaSNheLVF/1rC9m9skOxCuq2s/Ipn50kHDtxAGa6muEwp2gmSddBffcWmWW1unexN7DNLc2lqc12DH73NbA3xUlFgzpdjKwZgzxhAnb+0+3zi7EhcgQ+IiAJ/T1VfPKsDnwRTI3kR+Srgjqp+SES+dOKuPw78aVX95yLy+4HvBn77/ueHD+JFgBdeeGE2tsxLBNftmBlWq2kEmmVIexn6A2R9E33iCrK1DeIogvrjLFG6RLpuZ9y4FAqiMrAu03I6Fa2mEXsc2z+Wqrk9Li8Sv34f7bTJbnRJ7u7gWw0kK+i9cZGdp2NQaN33LH9SyVvCcFHYfqMgBbgMhstCNIpItz3JJ2+hq4vjgd/YJlRtbof550xIHEvtPllW+e4D1pQ1GIbn7C1ci4sPf43HADOXqgFOYVB2VUQ+OHH9xQNI/O2q+rqIXAdeEpGPq+pPnNlSj8E0I/m3A18tIr8TK7Iuisg/Bn43lqcH+AH2pnJqnBVaTXRrB3o9y3tnmRUc765Z5HnLUjOHDr54WIgjWlwwM7PSUsA5M+jKcku55DYXlXbLmoaGI1POZLmRvhMrFt/fQpe6oEpyv4fvmr1BdrXFvc+PEYXlTxRkXUc0UtTBcAWWP6EMVoTuawWNjYyoH7xyri7jW0m1zvjaFXSUES8tUNxfOzjCnnDDdGmK5plF5YUNFKly76MJF819n+fjnI+fTYLHYu+T2RrcOy7Prqqvh593ROT9wBcDs0/yqvptwLcBhEj+m1X1G0TkY8B/Afw48OXAJ6a1hscZur5ZRZdlUdX/3F0I5lkHOj+WqGyBT9mcI46o064MxKTRgMUuxVKLbDGl+ZlN3E4fDSP4RNXSM6V/zHAYtP2hiShykHvyK22i3RFFK2bz+QZbzwudm5DsKsmup2g6BsuOlV8e0r7t0EjovJ4zXEnImxFFI0IFUhGi3SHy69403lyKAu3lhw8qmVAGlV2vlZa9mjE7kc457mxgzrGf1Cd18rNG+GeRcxeRDuBUdTv8/k7gLz/ygU+Bi9DJ/3fA3xaRGBgQ8u41zhZaFAiEgqCvokzfn7AAOCAdMTm67jTzWyVOzId9ccHOIhpxZfAUre2iYtE5YV3EEeockllOm6ArL4d8kCRBG+9QJwyvd8hbjv51oXUHooHSvp3jU6F/VWhsqM1o9VA0hLwV4TIlHhSMFmOS7QJRJVtpkX78VYrnnyR6+Rbay/ekb/a8p9Cc5Ecji8ZLJ80kto0yfDZHyiIfE3Lfj4OaoGauC/ZscvI3gPeLqcJi4H2q+iNnceCT4lxIXlV/HIvcUdWfBL7oPF73cYYP7oZHYh+529SmCZfF4xpzwvNcmBYlC11Y6OAbiTUgNRKIHajiRrkpWhYXzRkyy4zsB0PbkDptK17u+8ey40Q0bu1SvLHL0qc8rbsj8o7l74s0YvkTGfFujiaOqJ8zXGkQZYLkSvqxV+FznsEnzpQ1T6xQPP8krjfCb21bM9ZBxdADTNVcEiOdlm1WsGcMYY0xjutynQnNvGJ9F496GNVPAW995AM9AuqO18cRB80ODbJAk/1xMvJSbxH80qIVVRc64ECTCN9JzbAr80RbfUu/THav7uyiO7vW3j/M0XCGUZmRZRn0TYGT3m9A7IiGSut2j41f12HhFUudpJsjVIS8m5AtWPTefaVPvN6D9U1IEuKNAflyk9vvegPXf2od2dyheO1W9b6PNFCrPh4bW6i7ts4zr2XUuGSYn8lQtXfN44JS3w0VOUuc7HGBPLBwuH8EX3k9DON23Q50OuhSB2KH7zQYXWkSbRsJu1FuqZg8N01+HMFwWI0YxCt+MER3+zblqcp329Qm3e0hn3oNt7ELqvSebLHyi1u4zFO0YvpPtIh6I5KNIe3X+zTWRkS//Gm4v25qnW6bX/2GJdyo4Pp/WEPuruNv3a0kncfWJgK0KPDDQXWpCf5BXPro/LQolV5HXWYAdSQ/zwgEvidSLbs0J1HquY/IKUucVMVYLcC1muZ348wjXQrFNyLc9pAkcqhzRDtDk2+2muYp0x+YyiZJINYwO3ZUNRT5fh9XTn0KtgNlh6u/eZuWCP03X+HeFy7RulfQ/eQm979oBSkW6PzCTVO8hHmp0u1As8Gnv/oav+677th9t+/ih0OIomrI9oHvuTzTOeq+Gg9gZs3IDoJitaE5QE3yc4pS/10NvggE7ZLYRvtNSP6Os9SNVleQlWX09l0b5EFSjbVjoWvk3Wzgtj2+3UAKjxtmlrvuD8ZRu6rJJEVsMIg4GzFY6s0LH6YqpbYpiE2F0v4At7TI9q+/RjRSln5tSHpnFxnlXPnp+5Yi2ty2jtN+H7ewgH/mOrd/8zLP/uAd2LUNRovCjMFKd8rDyPooEi9VRzXRzznmp95Sp2vmFD/a/0d7vqQSJxbBixt3aSbx2Gv+AEicEL/lzcjqCoTRgTDupiWOLA0TOYqrC+CN3CX3aBoHGWQ0tgUWG0hCFhQ7C11koVt5yaN+rD/PcnN07JmmvnjDNdKtHB9BvD1C7m/aJtAfIjv9Stbo3vgG/DM3uPdFyyx9KoMtW7cPEf6Jm5Im01v7b6/xeKBO19SYBUgUmWd70J2Lc0agBDnlITl3V86B7fXRrW3bGMDUJU4sgt/tVc6P0WbfHCIxRYwMc7TbQnaH5lef52ECU2GXNDWy77SRwcDy8XFcaeUBk1veuEb21LI1RK31abyWoXfu2VBxVds4ssxM15av4OOYwTNdpFBar9hG4O+vjYeGH0XukxH6w0T5NU6Emcjdn5G65jKgJvk5hmu3cZ12NepOnBvnuCdy0nuae4JjJZFF6X5t3RqUOiF6V7Uc+84utNtj22Kv0Gog/aEN8Sij91IHXw7VXugY6ceRFWR7oeAK1iXbH9iUpTSxTaHXJ15LyK52ceu7+Dt3w9mDNXpJmiCry2i3RXalQ9TLibczrn56C7m/joaGpwdqE5N4zBqWzgszQeZHYUYi9eNQk/w8Q9U6XjNLqWho6DHCD52aAZKmUBS469fQ+2tGtsGSWJqpkbiIEXx/YASfZTD01WtJ5IzgMRmlDALBD4Zj+WSvb7bCpW/9c08j9zasWBo3cFdXx5OhBsMww7VHurmLbm5VnuwuNClFaYJ2W2gjwY08bjAiureF3l/Dl41gpd69js5rnAZzQvJ1gnGO8aO7/7CK2CVJKk24RfRjPxZJU9xCF3ftKv7m7bEfS5ZZgTQoXnQ4gp3QZNXrobvh92CARuHNLmBgrpH0B3ZbqX0vC7DbOzDKkMUF5NY9e43FRfJnr0G7hS600VaKri6OzxRGIyP8KMI1G6bEaTaQJ67ju03c7pD4ldvIzgB/83b1vksP+5rIzxdzEcUXxfGXGUAdyc85XLNhw0C2tiuiK3Ysv+7SFJwgcYzftNmmGubCSmre7u7aVRgMbYxdbMVUzfPxC5SpHe8tDZMkpoMvi6uDIXRaRtRRZA1QpbXx7bu40usmy/CJY/jsMlknonnXFDBxVthZRFHYmURQx7h2Gz73TYwWUpK7O7C2AV4p7t63jSyKqvdTE/z5YuYJvsScRPI1yc85ip2xKmbShmDSjwUd+9lIFJnqZqELzuHXNyxtM/LmAV/YMHBpNqDdCDr58DUqCTUoaHSUmV49FHq11zd/+8hR3Lpjxd0nr1tnarPB1puaLP1qn9gJw6sN2q9s2yYRuZAOivBFQbS0iKwu07/aIt7NkLUtivXNvRG7ynl9xDUmMDcEDzXJ15gNvOR/YNykEjTePsv3KkmCfp4oskg9y/Fb25ajX1iohm1XCpvFBUrPd1pNuLc+Ht0Hls8vUzy9nrlRqiJLC7DTQ4dD4nJa1eY2dNr0P+sqq7+wTb6Q2kSnm31kJ6h/khh0CFdXLAe/uohvprQ+8hm01yefOEuBCXO1Ok1zLpgrYq+gc6OuqXPyjwH2/BMeIhH0o5HZC4R0jTUitVD1SKtpxdqrq0irib+xYiQ/GKK371n6pvDW/DTKrIGpN9au027B0iKMMjNOi2MbJHJ3Db2yxM7nX0e8ki03iLdHNF/bQrKCrS94guwNq9DrM3ruGr7bZPh5TyO7Q/Tnf5n85i2Kzc0H3ot6rebb1pgu5pPgCdY1/tjLLKCO5B83lHLFYFNQdoFKFBnBp8GPpgiDvhe6+G4LTSKiO+vQaOBur+Pvr+EWFiAP3vTe0jTSsO5VWVwwp8kozG3dMM26W1zAP3UVdY6dN3ZAhHi3wA0LomGBDDNkd8idd17lxr96BZyQv+Ea6c0t2Nwi/vAa+VFSyMPsCGpMBTPhKPmwqG0NaswiXJpOyAqLPc1QLk2N7IcjZGW50rq7u+voyoJF7v2BKVfUk9+/T9TtWtTf7RjR56FQOhpZoXSUmVQyipBWk97b3gDA6781oftp6NwqEMUKwIMc307Jnlzg6nf/dEXm0c4u+fqmLfI4y4Ea545Joj/Iv2YmNwFV+z7PAWqSf0zwkv8B3hF9HX442DPGTouiMiyTNLHUzELX9OxFYdFMs4H+6qftFHUi8netFsQxrtWsPNZJE9POe4U0DsdNYaHL4C3X2Hpjgo+hdQdcDhoJw5WI9us5OIhub8DP3trTuFSsrR/8puompkuDo8zJZjbarwuvNWYNLxXfz+9o/eFKVVMRqRMb+rHYRdc3kcJbsbXVNOuDXi/4zNvjXWpROar2s5zXWnjY3jFNfjnWr/DQSBk9d4XN51P6NyDZBJ8AXYhGQrrtwQlufRfd2EKLYjZJ4QwwV06OE5hFotc6kq8xi6iGXYizCLywqFuSGH9vzbxtRiN8f4AURZWakSjipeL7L3r5c4/zIsJ53UzODrNjQHYcpk7yIhIBHwReU9WvEpHvBz473L0MbKjq26a9jhqGmqhrwNltJnO7WdQGZafCnwQ+BiwCqOrXlXeIyN8ANs9hDTVq1JgCZi0Fc1Iox89ZmBVMVScvIs8Avwv4rgPuE+D3A983zTXUqFGjxqmhOpbjHnU5AUTkXSLyyyLyqyLyrVNe+QOYdjPU3wK+BTjo0/itwG1V/cRBTxSR94jIB0Xkg3fv3p3iEmvUqFHjQajXYy/HIaSr/3fgK4HPA75eRD5vykvfg6mRvIh8FXBHVT90yEO+niOieFV9UVVfUNUXrl27NpU11qhRo8ahOJtI/ouBX1XVT6nqCPgnwNdMdd37MM2c/NuBrxaR3wk0gUUR+ceq+g0iEgO/D/iikxzoQx/60D0ReWWKa31YXAXuXfQiDkG9tofHZV5fvbbj8cZHPcA26z/6r/WfXT3BQ5si8sGJ6y+q6osT158GPjNx/VXgP3/U9Z0GUyN5Vf024NsARORLgW9W1W8Id/924OOq+uoJj3UpQ3kR+aCqvnDR6zgI9doeHpd5ffXazgeq+q4zOtRBdqjnKtu5KIOyP0BdcK1Ro8b841XgDRPXnwFeP88FnEszlKr+OPDjE9f/6/N43Ro1atS4YPwn4C0i8jzwGhbg/sHzXEDd8fpoePH4h1wY6rU9PC7z+uq1zRBUNReRbwJ+FIiA71HVXzzPNYjOSetujRo1atR4EPXQkBo1atSYY9QkX6NGjRpzjJrkHxIiEonIz4nID1/0WvZDRJZF5J+JyMdF5GMi8psvek0lRORPi8gvishHReT7RKR5gWv5HhG5IyIfnbhtVUReEpFPhJ8rl2x9fz38XT8iIu8XkeXLsraJ+75ZRFRETqIzrzFl1CT/8CiN1y4j/jbwI6r6OcBbuSTrFJGngf8BeEFVPx8rRP2BC1zSPwD266G/FfgxVX0L8GPh+kXhH/Dg+l4CPl9VfwPwK4RelAvAP+DBtSEibwDeAXz6vBdU42DUJP8QOMp47aIhIovAbwO+G0BVR6q6caGL2osYaIWu5zbnrBmehKr+BLC27+avAd4bfn8v8HvOc02TOGh9qvoBVc3D1f+I6a7PHYd8dgDfiflV1YqOS4Ka5B8Of4vDjdcuGm8C7gL/R0gnfZeIdC56UQCq+hrwHViUdxPYVNUPXOyqHsANVb0JEH5ev+D1HIU/AvxfF72IEiLy1djciJ+/6LXUGKMm+VPiBMZrF40Y+ELg76jqFwC7XGzKoULIb38N8DzwFNARkW84+lk1DoKI/HkgB773otcCICJt4M8Df+Gi11JjL2qSPz1K47WXMUe5LxeRf3yxS9qDV4FXVfWnw/V/hpH+ZcBvB35NVe+qagb8C+BLLnhN+3FbRJ4ECD/vXPB6HoCIvBv4KuAP6eVpdHkztnn/fPjfeAb4WRF54kJXVaMm+dNCVb9NVZ9R1eewouG/mTBeu3Co6i3gMyJSjlj8CuCXLnBJk/g08JtEpB2GxnwFl6QoPIEfAt4dfn838IMXuJYHICLvAv4s8NWq2rvo9ZRQ1V9Q1euq+lz433gV+MLwfaxxgahJfj7x/wS+V0Q+ArwN+F8vdjmGcHbxz4CfBX4B+/5dWCu8iHwf8FPAZ4vIqyLyR4FvB94hIp/AVCLffsnW978BC8BLIvJhEfm7l2htNS4haluDGjVq1Jhj1JF8jRo1aswxapKvUaNGjTlGTfI1atSoMceoSb5GjRo15hg1ydeoUaPGHKMm+RqXGiKyc9FrqFFjllGTfI0aNWrMMWqSrzETEMNfDz70vyAiXxdu/1IR+fEJ//zvDd20NWrUoB7kXWN28Puw7t23AleB/yQiPxHu+wLgP8Nsi/895i/0kxewxho1Lh3qSL7GrOC3AN+nqoWq3gb+HfAbw30/o6qvqqoHPgw8dzFLrFHj8qEm+RqzgqNSMMOJ3wvqM9QaNSrUJF9jVvATwNeF2brXsOlXP3PBa6pR49KjjnhqzAreD/xm4Oex0XLfoqq3RORzLnZZNWpcbtQulDVq1Kgxx6jTNTVq1Kgxx6hJvkaNGjXmGDXJ16hRo8Ycoyb5GjVq1Jhj1CRfo0aNGnOMmuRr1KhRY45Rk3yNGjVqzDH+/xWR6FvRJUZDAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 2 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "import matplotlib.pyplot as plt\n",
    "dsg.RW.plot(x=\"lon\", y=\"lat\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Get Point Coordinate\n",
    "\n",
    "Adapted from [StackOverflow](https://stackoverflow.com/a/58774123).\n",
    "\n",
    "There are other solutions, but they involve further packages, see [xarray-discussions](https://github.com/pydata/xarray/discussions/5009).\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "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[data-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 !important;\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.DataArray &#x27;RW&#x27; ()&gt;\n",
       "array(1.1, dtype=float32)\n",
       "Coordinates:\n",
       "    y        float64 -4.326e+06\n",
       "    x        float64 37.83\n",
       "    lon      float64 9.994\n",
       "    lat      float64 50.0\n",
       "Attributes:\n",
       "    valid_min:      0\n",
       "    valid_max:      4095\n",
       "    standard_name:  rainfall_rate\n",
       "    long_name:      RW\n",
       "    unit:           mm h-1</pre><div class='xr-wrap' style='display:none'><div class='xr-header'><div class='xr-obj-type'>xarray.DataArray</div><div class='xr-array-name'>'RW'</div></div><ul class='xr-sections'><li class='xr-section-item'><div class='xr-array-wrap'><input id='section-c2b19131-3ed5-43cc-87d2-0668faee6332' class='xr-array-in' type='checkbox' checked><label for='section-c2b19131-3ed5-43cc-87d2-0668faee6332' title='Show/hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-array-preview xr-preview'><span>1.1</span></div><div class='xr-array-data'><pre>array(1.1, dtype=float32)</pre></div></div></li><li class='xr-section-item'><input id='section-318c8dd7-ab5e-4a2f-8087-e0efd659d0a2' class='xr-section-summary-in' type='checkbox'  checked><label for='section-318c8dd7-ab5e-4a2f-8087-e0efd659d0a2' class='xr-section-summary' >Coordinates: <span>(4)</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>y</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>-4.326e+06</div><input id='attrs-d0a4ca16-5a6b-4b86-ad68-4d21dc8b7b29' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-d0a4ca16-5a6b-4b86-ad68-4d21dc8b7b29' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-77c29c0d-111d-4c14-9af2-7eb8721605e0' class='xr-var-data-in' type='checkbox'><label for='data-77c29c0d-111d-4c14-9af2-7eb8721605e0' 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>units :</span></dt><dd>m</dd><dt><span>long_name :</span></dt><dd>y coordinate of projection</dd><dt><span>standard_name :</span></dt><dd>projection_y_coordinate</dd></dl></div><div class='xr-var-data'><pre>array(-4326144.72426557)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>x</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>37.83</div><input id='attrs-02f10e2a-a18e-46ec-9508-df459dd228c1' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-02f10e2a-a18e-46ec-9508-df459dd228c1' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-0bcf18fe-5a5d-495d-bec3-840bc1c2d309' class='xr-var-data-in' type='checkbox'><label for='data-0bcf18fe-5a5d-495d-bec3-840bc1c2d309' 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>units :</span></dt><dd>m</dd><dt><span>long_name :</span></dt><dd>x coordinate of projection</dd><dt><span>standard_name :</span></dt><dd>projection_x_coordinate</dd></dl></div><div class='xr-var-data'><pre>array(37.83307814)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>lon</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>9.994</div><input id='attrs-8c35a67e-28df-4183-9471-0772d207996b' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-8c35a67e-28df-4183-9471-0772d207996b' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-630f06db-028e-45b0-b9d7-31ac61c279d7' class='xr-var-data-in' type='checkbox'><label for='data-630f06db-028e-45b0-b9d7-31ac61c279d7' 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(9.99387973)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>lat</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>50.0</div><input id='attrs-2ffcc16d-470c-4d65-a09f-940d6d0b910c' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-2ffcc16d-470c-4d65-a09f-940d6d0b910c' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-f1b6e837-3208-4ae8-87c3-40120e88e714' class='xr-var-data-in' type='checkbox'><label for='data-f1b6e837-3208-4ae8-87c3-40120e88e714' 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(49.99782409)</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-7ad2572b-10b1-4715-8ca8-60d3e68e20e9' class='xr-section-summary-in' type='checkbox'  checked><label for='section-7ad2572b-10b1-4715-8ca8-60d3e68e20e9' class='xr-section-summary' >Attributes: <span>(5)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><dl class='xr-attrs'><dt><span>valid_min :</span></dt><dd>0</dd><dt><span>valid_max :</span></dt><dd>4095</dd><dt><span>standard_name :</span></dt><dd>rainfall_rate</dd><dt><span>long_name :</span></dt><dd>RW</dd><dt><span>unit :</span></dt><dd>mm h-1</dd></dl></div></li></ul></div></div>"
      ],
      "text/plain": [
       "<xarray.DataArray 'RW' ()>\n",
       "array(1.1, dtype=float32)\n",
       "Coordinates:\n",
       "    y        float64 -4.326e+06\n",
       "    x        float64 37.83\n",
       "    lon      float64 9.994\n",
       "    lat      float64 50.0\n",
       "Attributes:\n",
       "    valid_min:      0\n",
       "    valid_max:      4095\n",
       "    standard_name:  rainfall_rate\n",
       "    long_name:      RW\n",
       "    unit:           mm h-1"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "lat = 50.0\n",
    "lon = 10.0\n",
    "\n",
    "# find nearest xy-grid point to a specific latlon-coordinate   \n",
    "abslat = np.abs(dsg.lat-lat)\n",
    "abslon = np.abs(dsg.lon-lon)\n",
    "c = np.maximum(abslon, abslat)\n",
    "\n",
    "# Attention: y/lat is first dim, get\n",
    "([yidx], [xidx]) = np.where(c == np.min(c))\n",
    "    \n",
    "# Select index location at the x/y dimension\n",
    "# use isel as we select with index \n",
    "point_ds = dsg.isel(x=xidx, y=yidx)\n",
    "\n",
    "display(point_ds.RW)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "{'source': '/home/kai/projects/data/wradlib-data/radolan/showcase/raa01-rw_10000-1408102050-dwd---bin.gz',\n",
       " '_FillValue': array([ 2490,  2500, 65535], dtype=int32),\n",
       " 'scale_factor': 0.1,\n",
       " 'add_offset': 0.0,\n",
       " 'dtype': dtype('uint16'),\n",
       " 'coordinates': 'time y x'}"
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "point_ds.RW.encoding"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "wradlib",
   "language": "python",
   "name": "wradlib"
  },
  "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.10.0"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}