scratch_plot_linear_trend.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "459530f8",
   "metadata": {},
   "outputs": [],
   "source": [
    "import datetime as dt\n",
    "import glob\n",
    "import numpy as np\n",
    "import os\n",
    "\n",
    "import pandas as pd\n",
    "import xarray as xr\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "7aee0111",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div><svg style=\"position: absolute; width: 0; height: 0; overflow: hidden\">\n",
       "<defs>\n",
       "<symbol id=\"icon-database\" viewBox=\"0 0 32 32\">\n",
       "<path d=\"M16 0c-8.837 0-16 2.239-16 5v4c0 2.761 7.163 5 16 5s16-2.239 16-5v-4c0-2.761-7.163-5-16-5z\"></path>\n",
       "<path d=\"M16 17c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n",
       "<path d=\"M16 26c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n",
       "</symbol>\n",
       "<symbol id=\"icon-file-text2\" viewBox=\"0 0 32 32\">\n",
       "<path d=\"M28.681 7.159c-0.694-0.947-1.662-2.053-2.724-3.116s-2.169-2.030-3.116-2.724c-1.612-1.182-2.393-1.319-2.841-1.319h-15.5c-1.378 0-2.5 1.121-2.5 2.5v27c0 1.378 1.122 2.5 2.5 2.5h23c1.378 0 2.5-1.122 2.5-2.5v-19.5c0-0.448-0.137-1.23-1.319-2.841zM24.543 5.457c0.959 0.959 1.712 1.825 2.268 2.543h-4.811v-4.811c0.718 0.556 1.584 1.309 2.543 2.268zM28 29.5c0 0.271-0.229 0.5-0.5 0.5h-23c-0.271 0-0.5-0.229-0.5-0.5v-27c0-0.271 0.229-0.5 0.5-0.5 0 0 15.499-0 15.5 0v7c0 0.552 0.448 1 1 1h7v19.5z\"></path>\n",
       "<path d=\"M23 26h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "<path d=\"M23 22h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "<path d=\"M23 18h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "</symbol>\n",
       "</defs>\n",
       "</svg>\n",
       "<style>/* CSS stylesheet for displaying xarray objects in jupyterlab.\n",
       " *\n",
       " */\n",
       "\n",
       ":root {\n",
       "  --xr-font-color0: var(--jp-content-font-color0, rgba(0, 0, 0, 1));\n",
       "  --xr-font-color2: var(--jp-content-font-color2, rgba(0, 0, 0, 0.54));\n",
       "  --xr-font-color3: var(--jp-content-font-color3, rgba(0, 0, 0, 0.38));\n",
       "  --xr-border-color: var(--jp-border-color2, #e0e0e0);\n",
       "  --xr-disabled-color: var(--jp-layout-color3, #bdbdbd);\n",
       "  --xr-background-color: var(--jp-layout-color0, white);\n",
       "  --xr-background-color-row-even: var(--jp-layout-color1, white);\n",
       "  --xr-background-color-row-odd: var(--jp-layout-color2, #eeeeee);\n",
       "}\n",
       "\n",
       "html[theme=dark],\n",
       "body.vscode-dark {\n",
       "  --xr-font-color0: rgba(255, 255, 255, 1);\n",
       "  --xr-font-color2: rgba(255, 255, 255, 0.54);\n",
       "  --xr-font-color3: rgba(255, 255, 255, 0.38);\n",
       "  --xr-border-color: #1F1F1F;\n",
       "  --xr-disabled-color: #515151;\n",
       "  --xr-background-color: #111111;\n",
       "  --xr-background-color-row-even: #111111;\n",
       "  --xr-background-color-row-odd: #313131;\n",
       "}\n",
       "\n",
       ".xr-wrap {\n",
       "  display: block;\n",
       "  min-width: 300px;\n",
       "  max-width: 700px;\n",
       "}\n",
       "\n",
       ".xr-text-repr-fallback {\n",
       "  /* fallback to plain text repr when CSS is not injected (untrusted notebook) */\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-header {\n",
       "  padding-top: 6px;\n",
       "  padding-bottom: 6px;\n",
       "  margin-bottom: 4px;\n",
       "  border-bottom: solid 1px var(--xr-border-color);\n",
       "}\n",
       "\n",
       ".xr-header > div,\n",
       ".xr-header > ul {\n",
       "  display: inline;\n",
       "  margin-top: 0;\n",
       "  margin-bottom: 0;\n",
       "}\n",
       "\n",
       ".xr-obj-type,\n",
       ".xr-array-name {\n",
       "  margin-left: 2px;\n",
       "  margin-right: 10px;\n",
       "}\n",
       "\n",
       ".xr-obj-type {\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-sections {\n",
       "  padding-left: 0 !important;\n",
       "  display: grid;\n",
       "  grid-template-columns: 150px auto auto 1fr 20px 20px;\n",
       "}\n",
       "\n",
       ".xr-section-item {\n",
       "  display: contents;\n",
       "}\n",
       "\n",
       ".xr-section-item input {\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-section-item input + label {\n",
       "  color: var(--xr-disabled-color);\n",
       "}\n",
       "\n",
       ".xr-section-item input:enabled + label {\n",
       "  cursor: pointer;\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-section-item input:enabled + label:hover {\n",
       "  color: var(--xr-font-color0);\n",
       "}\n",
       "\n",
       ".xr-section-summary {\n",
       "  grid-column: 1;\n",
       "  color: var(--xr-font-color2);\n",
       "  font-weight: 500;\n",
       "}\n",
       "\n",
       ".xr-section-summary > span {\n",
       "  display: inline-block;\n",
       "  padding-left: 0.5em;\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:disabled + label {\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-section-summary-in + label:before {\n",
       "  display: inline-block;\n",
       "  content: '►';\n",
       "  font-size: 11px;\n",
       "  width: 15px;\n",
       "  text-align: center;\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:disabled + label:before {\n",
       "  color: var(--xr-disabled-color);\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:checked + label:before {\n",
       "  content: '▼';\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:checked + label > span {\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-section-summary,\n",
       ".xr-section-inline-details {\n",
       "  padding-top: 4px;\n",
       "  padding-bottom: 4px;\n",
       "}\n",
       "\n",
       ".xr-section-inline-details {\n",
       "  grid-column: 2 / -1;\n",
       "}\n",
       "\n",
       ".xr-section-details {\n",
       "  display: none;\n",
       "  grid-column: 1 / -1;\n",
       "  margin-bottom: 5px;\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:checked ~ .xr-section-details {\n",
       "  display: contents;\n",
       "}\n",
       "\n",
       ".xr-array-wrap {\n",
       "  grid-column: 1 / -1;\n",
       "  display: grid;\n",
       "  grid-template-columns: 20px auto;\n",
       "}\n",
       "\n",
       ".xr-array-wrap > label {\n",
       "  grid-column: 1;\n",
       "  vertical-align: top;\n",
       "}\n",
       "\n",
       ".xr-preview {\n",
       "  color: var(--xr-font-color3);\n",
       "}\n",
       "\n",
       ".xr-array-preview,\n",
       ".xr-array-data {\n",
       "  padding: 0 5px !important;\n",
       "  grid-column: 2;\n",
       "}\n",
       "\n",
       ".xr-array-data,\n",
       ".xr-array-in:checked ~ .xr-array-preview {\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-array-in:checked ~ .xr-array-data,\n",
       ".xr-array-preview {\n",
       "  display: inline-block;\n",
       "}\n",
       "\n",
       ".xr-dim-list {\n",
       "  display: inline-block !important;\n",
       "  list-style: none;\n",
       "  padding: 0 !important;\n",
       "  margin: 0;\n",
       "}\n",
       "\n",
       ".xr-dim-list li {\n",
       "  display: inline-block;\n",
       "  padding: 0;\n",
       "  margin: 0;\n",
       "}\n",
       "\n",
       ".xr-dim-list:before {\n",
       "  content: '(';\n",
       "}\n",
       "\n",
       ".xr-dim-list:after {\n",
       "  content: ')';\n",
       "}\n",
       "\n",
       ".xr-dim-list li:not(:last-child):after {\n",
       "  content: ',';\n",
       "  padding-right: 5px;\n",
       "}\n",
       "\n",
       ".xr-has-index {\n",
       "  font-weight: bold;\n",
       "}\n",
       "\n",
       ".xr-var-list,\n",
       ".xr-var-item {\n",
       "  display: contents;\n",
       "}\n",
       "\n",
       ".xr-var-item > div,\n",
       ".xr-var-item label,\n",
       ".xr-var-item > .xr-var-name span {\n",
       "  background-color: var(--xr-background-color-row-even);\n",
       "  margin-bottom: 0;\n",
       "}\n",
       "\n",
       ".xr-var-item > .xr-var-name:hover span {\n",
       "  padding-right: 5px;\n",
       "}\n",
       "\n",
       ".xr-var-list > li:nth-child(odd) > div,\n",
       ".xr-var-list > li:nth-child(odd) > label,\n",
       ".xr-var-list > li:nth-child(odd) > .xr-var-name span {\n",
       "  background-color: var(--xr-background-color-row-odd);\n",
       "}\n",
       "\n",
       ".xr-var-name {\n",
       "  grid-column: 1;\n",
       "}\n",
       "\n",
       ".xr-var-dims {\n",
       "  grid-column: 2;\n",
       "}\n",
       "\n",
       ".xr-var-dtype {\n",
       "  grid-column: 3;\n",
       "  text-align: right;\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-var-preview {\n",
       "  grid-column: 4;\n",
       "}\n",
       "\n",
       ".xr-var-name,\n",
       ".xr-var-dims,\n",
       ".xr-var-dtype,\n",
       ".xr-preview,\n",
       ".xr-attrs dt {\n",
       "  white-space: nowrap;\n",
       "  overflow: hidden;\n",
       "  text-overflow: ellipsis;\n",
       "  padding-right: 10px;\n",
       "}\n",
       "\n",
       ".xr-var-name:hover,\n",
       ".xr-var-dims:hover,\n",
       ".xr-var-dtype:hover,\n",
       ".xr-attrs dt:hover {\n",
       "  overflow: visible;\n",
       "  width: auto;\n",
       "  z-index: 1;\n",
       "}\n",
       "\n",
       ".xr-var-attrs,\n",
       ".xr-var-data {\n",
       "  display: none;\n",
       "  background-color: var(--xr-background-color) !important;\n",
       "  padding-bottom: 5px !important;\n",
       "}\n",
       "\n",
       ".xr-var-attrs-in:checked ~ .xr-var-attrs,\n",
       ".xr-var-data-in:checked ~ .xr-var-data {\n",
       "  display: block;\n",
       "}\n",
       "\n",
       ".xr-var-data > table {\n",
       "  float: right;\n",
       "}\n",
       "\n",
       ".xr-var-name span,\n",
       ".xr-var-data,\n",
       ".xr-attrs {\n",
       "  padding-left: 25px !important;\n",
       "}\n",
       "\n",
       ".xr-attrs,\n",
       ".xr-var-attrs,\n",
       ".xr-var-data {\n",
       "  grid-column: 1 / -1;\n",
       "}\n",
       "\n",
       "dl.xr-attrs {\n",
       "  padding: 0;\n",
       "  margin: 0;\n",
       "  display: grid;\n",
       "  grid-template-columns: 125px auto;\n",
       "}\n",
       "\n",
       ".xr-attrs dt,\n",
       ".xr-attrs dd {\n",
       "  padding: 0;\n",
       "  margin: 0;\n",
       "  float: left;\n",
       "  padding-right: 10px;\n",
       "  width: auto;\n",
       "}\n",
       "\n",
       ".xr-attrs dt {\n",
       "  font-weight: normal;\n",
       "  grid-column: 1;\n",
       "}\n",
       "\n",
       ".xr-attrs dt:hover span {\n",
       "  display: inline-block;\n",
       "  background: var(--xr-background-color);\n",
       "  padding-right: 10px;\n",
       "}\n",
       "\n",
       ".xr-attrs dd {\n",
       "  grid-column: 2;\n",
       "  white-space: pre-wrap;\n",
       "  word-break: break-all;\n",
       "}\n",
       "\n",
       ".xr-icon-database,\n",
       ".xr-icon-file-text2 {\n",
       "  display: inline-block;\n",
       "  vertical-align: middle;\n",
       "  width: 1em;\n",
       "  height: 1.5em !important;\n",
       "  stroke-width: 0;\n",
       "  stroke: currentColor;\n",
       "  fill: currentColor;\n",
       "}\n",
       "</style><pre class='xr-text-repr-fallback'>&lt;xarray.DataArray &#x27;siconca&#x27; (time: 31, yc: 432, xc: 432)&gt;\n",
       "[5785344 values with dtype=float32]\n",
       "Coordinates:\n",
       "  * time     (time) datetime64[ns] 2010-01-05 2010-01-06 ... 2010-02-04\n",
       "    lat      (yc, xc) float32 ...\n",
       "    lon      (yc, xc) float32 ...\n",
       "  * xc       (xc) float64 -5.388e+03 -5.362e+03 ... 5.362e+03 5.388e+03\n",
       "  * yc       (yc) float64 5.388e+03 5.362e+03 ... -5.362e+03 -5.388e+03</pre><div class='xr-wrap' hidden><div class='xr-header'><div class='xr-obj-type'>xarray.DataArray</div><div class='xr-array-name'>'siconca'</div><ul class='xr-dim-list'><li><span class='xr-has-index'>time</span>: 31</li><li><span class='xr-has-index'>yc</span>: 432</li><li><span class='xr-has-index'>xc</span>: 432</li></ul></div><ul class='xr-sections'><li class='xr-section-item'><div class='xr-array-wrap'><input id='section-8f42f5df-5832-46d1-8240-e79f30077357' class='xr-array-in' type='checkbox' checked><label for='section-8f42f5df-5832-46d1-8240-e79f30077357' 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>...</span></div><div class='xr-array-data'><pre>[5785344 values with dtype=float32]</pre></div></div></li><li class='xr-section-item'><input id='section-d2ae2bd9-5997-4aa2-86d7-32f21266dd43' class='xr-section-summary-in' type='checkbox'  checked><label for='section-d2ae2bd9-5997-4aa2-86d7-32f21266dd43' 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'>2010-01-05 ... 2010-02-04</div><input id='attrs-b520f809-5df9-4cb9-bf21-e870b3c28f04' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-b520f809-5df9-4cb9-bf21-e870b3c28f04' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-657329a0-6daf-46c4-9003-b0be68374632' class='xr-var-data-in' type='checkbox'><label for='data-657329a0-6daf-46c4-9003-b0be68374632' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([&#x27;2010-01-05T00:00:00.000000000&#x27;, &#x27;2010-01-06T00:00:00.000000000&#x27;,\n",
       "       &#x27;2010-01-07T00:00:00.000000000&#x27;, &#x27;2010-01-08T00:00:00.000000000&#x27;,\n",
       "       &#x27;2010-01-09T00:00:00.000000000&#x27;, &#x27;2010-01-10T00:00:00.000000000&#x27;,\n",
       "       &#x27;2010-01-11T00:00:00.000000000&#x27;, &#x27;2010-01-12T00:00:00.000000000&#x27;,\n",
       "       &#x27;2010-01-13T00:00:00.000000000&#x27;, &#x27;2010-01-14T00:00:00.000000000&#x27;,\n",
       "       &#x27;2010-01-15T00:00:00.000000000&#x27;, &#x27;2010-01-16T00:00:00.000000000&#x27;,\n",
       "       &#x27;2010-01-17T00:00:00.000000000&#x27;, &#x27;2010-01-18T00:00:00.000000000&#x27;,\n",
       "       &#x27;2010-01-19T00:00:00.000000000&#x27;, &#x27;2010-01-20T00:00:00.000000000&#x27;,\n",
       "       &#x27;2010-01-21T00:00:00.000000000&#x27;, &#x27;2010-01-22T00:00:00.000000000&#x27;,\n",
       "       &#x27;2010-01-23T00:00:00.000000000&#x27;, &#x27;2010-01-24T00:00:00.000000000&#x27;,\n",
       "       &#x27;2010-01-25T00:00:00.000000000&#x27;, &#x27;2010-01-26T00:00:00.000000000&#x27;,\n",
       "       &#x27;2010-01-27T00:00:00.000000000&#x27;, &#x27;2010-01-28T00:00:00.000000000&#x27;,\n",
       "       &#x27;2010-01-29T00:00:00.000000000&#x27;, &#x27;2010-01-30T00:00:00.000000000&#x27;,\n",
       "       &#x27;2010-01-31T00:00:00.000000000&#x27;, &#x27;2010-02-01T00:00:00.000000000&#x27;,\n",
       "       &#x27;2010-02-02T00:00:00.000000000&#x27;, &#x27;2010-02-03T00:00:00.000000000&#x27;,\n",
       "       &#x27;2010-02-04T00:00:00.000000000&#x27;], dtype=&#x27;datetime64[ns]&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>lat</span></div><div class='xr-var-dims'>(yc, xc)</div><div class='xr-var-dtype'>float32</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-e5ba8db1-43ee-4253-8080-ea9da4365fc0' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-e5ba8db1-43ee-4253-8080-ea9da4365fc0' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-4bee12c1-e84c-470f-81e3-057308aae5ae' class='xr-var-data-in' type='checkbox'><label for='data-4bee12c1-e84c-470f-81e3-057308aae5ae' 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>degrees_north</dd><dt><span>long_name :</span></dt><dd>latitude coordinate</dd><dt><span>standard_name :</span></dt><dd>latitude</dd></dl></div><div class='xr-var-data'><pre>[186624 values with dtype=float32]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>lon</span></div><div class='xr-var-dims'>(yc, xc)</div><div class='xr-var-dtype'>float32</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-1bccb8b8-4009-4158-81f8-2a2c5b115638' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-1bccb8b8-4009-4158-81f8-2a2c5b115638' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-81f3d03a-88f1-4a14-ad3f-c7ab9b5a17d8' class='xr-var-data-in' type='checkbox'><label for='data-81f3d03a-88f1-4a14-ad3f-c7ab9b5a17d8' 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>degrees_east</dd><dt><span>long_name :</span></dt><dd>longitude coordinate</dd><dt><span>standard_name :</span></dt><dd>longitude</dd></dl></div><div class='xr-var-data'><pre>[186624 values with dtype=float32]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>xc</span></div><div class='xr-var-dims'>(xc)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>-5.388e+03 -5.362e+03 ... 5.388e+03</div><input id='attrs-08e851cb-d955-469f-b72d-e718b9b267fb' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-08e851cb-d955-469f-b72d-e718b9b267fb' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-882ec5b9-1239-448a-a241-f5247f444715' class='xr-var-data-in' type='checkbox'><label for='data-882ec5b9-1239-448a-a241-f5247f444715' 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>km</dd><dt><span>long_name :</span></dt><dd>x coordinate of projection (eastings)</dd><dt><span>standard_name :</span></dt><dd>projection_x_coordinate</dd></dl></div><div class='xr-var-data'><pre>array([-5387.5, -5362.5, -5337.5, ...,  5337.5,  5362.5,  5387.5])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>yc</span></div><div class='xr-var-dims'>(yc)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>5.388e+03 5.362e+03 ... -5.388e+03</div><input id='attrs-a27d2c03-f882-4bf1-bb76-a119a8ac8c9c' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-a27d2c03-f882-4bf1-bb76-a119a8ac8c9c' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-7b5794af-dc14-4c65-b91f-6bdced323f9c' class='xr-var-data-in' type='checkbox'><label for='data-7b5794af-dc14-4c65-b91f-6bdced323f9c' 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>km</dd><dt><span>long_name :</span></dt><dd>y coordinate of projection (northings)</dd><dt><span>standard_name :</span></dt><dd>projection_y_coordinate</dd></dl></div><div class='xr-var-data'><pre>array([ 5387.5,  5362.5,  5337.5, ..., -5337.5, -5362.5, -5387.5])</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-09ddcc54-57cd-4e1d-b003-dde7352ab2cc' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-09ddcc54-57cd-4e1d-b003-dde7352ab2cc' class='xr-section-summary'  title='Expand/collapse section'>Attributes: <span>(0)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><dl class='xr-attrs'></dl></div></li></ul></div></div>"
      ],
      "text/plain": [
       "<xarray.DataArray 'siconca' (time: 31, yc: 432, xc: 432)>\n",
       "[5785344 values with dtype=float32]\n",
       "Coordinates:\n",
       "  * time     (time) datetime64[ns] 2010-01-05 2010-01-06 ... 2010-02-04\n",
       "    lat      (yc, xc) float32 ...\n",
       "    lon      (yc, xc) float32 ...\n",
       "  * xc       (xc) float64 -5.388e+03 -5.362e+03 ... 5.362e+03 5.388e+03\n",
       "  * yc       (yc) float64 5.388e+03 5.362e+03 ... -5.362e+03 -5.388e+03"
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "file = \"../lt_arr.nc\"\n",
    "arr = xr.open_dataset(file)\n",
    "arr = arr.siconca\n",
    "arr"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "daf6c607",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.contour.QuadContourSet at 0x7f4207ad2df0>"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZ0AAAEiCAYAAAAiQw8CAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAABFyklEQVR4nO2de5gcVZn/P18CwoiMinhhgHCZgCa4RgXJ+HOTqIgDeMGoG0AUGNZNXOMl6y2wrqvrLrtGXY2XCGSVAK6KAYmii0ZAScbFDhcxXIJoBjDGQVgUCWBEAu/vj3Oqp7qnurt6prune+b9PE8903XqVNWp7uR86z3nPe8rM8NxHMdxWsEuE90Ax3EcZ+rgouM4juO0DBcdx3Ecp2W46DiO4zgtw0XHcRzHaRkuOo7jOE7LcNFxWoak6ZIeljRtotviOM7E4KLjNA1Jd0t6VbJvZlvN7Clm9vhEtqsSknaX9BVJv5b0kKSbJB1XVudoSb+Q9CdJP5Z0YOrYK2LZg5Luzrj+QfH4n+I1XlVeJ299Sa+R9BNJf5T0O0n/JWmvGtd7S3y2RyR9W9LeqWOflvSr+Ny/kHRqjq/McerGRcdxRtgV+A0wH3gq8BFgjaSDACTtA1wWy/cGbgC+mTr/EeB84IMVrv8N4CbgGcCHgUslPbNKe6rVfyrwb0APMBPYH/hUpQtJOhw4D3gb8GzgT8CXytr+unjd04DPSfp/VdrmOGPDzHzzreEb8FXgCWAH8DDwIeAgwIBdY51rCB3ntbHOdwkd7NeA7cD1wEGpaz4PuBL4A3AHsLAFz3Ez8Kb4eRFwberYnvH5nld2zquAu8vKDgMeBfZKlQ0C76hw33rrvxG4pcpz/Dvw9dR+L/CX9PXL6l8OvH+i/x35Nvk2t3ScpmBmbwO2Aq+zMKT2yQpVTyK8fe9H6Ah/CqwmWBK3Ax8FkLQnQXC+DjwLOBn4UnyDH4WkL8Whp6zt5jzPIOnZhM7/tlh0OLAp9YyPAEOxvBaHA3ea2UOpsk1Vzq23/rxUOytdL932IYLoHFZeUVIX8JIa13OcMeGi40w0q81syMweBL4PDJnZVWa2E7gEeFGs91qC9bDazHaa2c+AbwFvzrqomb3TzJ5WYXtBrUZJ2o1gcV1oZr+IxU8BHiyr+iBQdS5ljOfmri/pGMKQ2D836P7nEgRqXZXrOc6YcNFxJpp7U593ZOw/JX4+EJiTtliAU4DnNLpBknYhDA/+BXhX6tDDQHdZ9W7gIWpT9VxJt0XPvoclzc17L0l9BOvvzWb2y1g2N3WtxFrJe71PAc8nDF16NGCn4ew60Q1wJjWN7LR+A6w3s2PyVJZ0LvDWCod/bWaVhuUEfIUw2X68mT2WOnwbwaJI6u5JGBLMMwx1G3CIpL1SQ2azCYJBeXskHVatfqzzIsLcyxlmdnVSbmaDjIh1+v6zU+ceAuwO/DJV9i/AccB8M9ue45kcp27c0nGayb3AIQ261veAwyS9TdJucXuJpJlZlc3sHXEuKWurNgdzDsEb7HVmtqPs2Frg+ZLeJGkPwnDWzcnwm6RdYvluYVd7SHpSbM8vgZ8DH43lC4AXEIYIs9pftb6k5wM/AN5tZt+t/tUBYajwddEK2hP4OHBZImiSzgLeAhxjZr/PcT3HGRMuOk4z+Q/gn+Jw2AfGc6HYOb6a4HgwDPwOWE54W28Icc3NYuCFwO9SQ1SnxDb8H/Am4GzgAWBObE/CPMKQ4BXA9Pj5h6njJwFHxnM/QRgS+78qTapW//3AM4GvZAyljcLMbgPeQRCf+whzOe9MVfn32OZfpa73j1Xa5jhjQj5s6ziO47QKt3Qcx3GcluGi4ziOM4mRdL6k+yTdWuG4JH1e0hZJN0t6cTPb46LjOI4zubkAOLbK8eOAQ+O2iOBM0zRcdBzHcSYxZraBEDqqEicAF1mgADxN0r7Nao+LjuM4ztRmP8I6uIRtsawp+OLQnDxJu9se7DnRzXAcpwN4iAfuN7NqEcSr0v+KPe33f8iXAeTGmx+9DfhzqmiVma2q43bKKGuaW7OLTk72YE/m6OiJbobjOB3AVXbpr8dz/u//8DjXrZueq+60fX/1ZzM7chy32wYckNrfn7AWrin48JrjOM7U5nLg1OjF1gc8aGb3NOtmbuk4juNMYiR9A3g5sI+kbYR0IbsBmNm5hAgaxwNbCMn9BprZHhcdx3GcSYyZnVzjuAFLWtQcH15zHMdxWoeLjuM4jtMyXHQcx3GcluGi4ziO47QMFx3HcRynZbjoOI7jOC3DRcdxHMdpGS46juM4Tstw0XEcx3FahouO4ziO0zJcdBzHcZyW4aLjOI7jtAwXHcdxHKdluOg4juM4LcNFx3Ecx2kZLjqO4zhOy3DRcRzHcVqGi47jOI7TMlx0HMdxnJbRcaIjaZqkmyR9L+7vLelKSb+Kf5+eqnuWpC2S7pDUnyo/QtIt8djnJWkinsVxJis7Fsypq2499Z3OpuNEB3gvcHtq/0zgajM7FLg67iNpFnAScDhwLPAlSdPiOecAi4BD43Zsa5ruOJOfHQvmMDxfuYVkeL7oXXZ77YrOpGDXiW5APUjaH3gNcDbwvlh8AvDy+PlC4BpgWSy/2MweBe6StAU4StLdQLeZ/TRe8yLgDcD3W/IQjjPJ6Vq7kd61+ev3rDcGmQUrwn7v0kJT2uW0Bx0lOoR/lh8C9kqVPdvM7gEws3skPSuW7wek//Vui2WPxc/l5Y7jNIAdC+bQtXYjAEMr+uhZb8X9coZW9AEuNFOJjhEdSa8F7jOzGyW9PM8pGWVWpTzrnosIw3DswZPzNdRxpjDJ0FoPYWitZ72F4bM6LB9nctMxogO8DHi9pOOBPYBuSf8N3Ctp32jl7AvcF+tvAw5Inb8/MBzL988oH4WZrQJWAXRr70xhchwnkAgOhHmaxMJJBCgLt3CmHh3jSGBmZ5nZ/mZ2EMFB4Edm9lbgcuC0WO004Dvx8+XASZJ2l3QwwWHgujgU95Ckvui1dmrqHMdxcpJ4ne1YMIehFX1FwSmna+3GMMxW6G5xC512pJMsnUp8Algj6W+BrcDfAJjZbZLWAJuBncASM3s8nvP3wAVAF8GBwJ0InElBej6l2fdJi8zcvs0lx1dPH2Te+kXFuuXHnalLR4qOmV1D8FLDzH4PHF2h3tkET7fy8huA5zevhY4zMTRz/iRxgU7EJi0kq6cPZp4ztKKP3qUFhn1Ox4l0pOg4jpNNo+dIsoSml8oik9DfM5suNsL8voa2x+l8XHQcp4NIRKBVQ2jJos08QpOFOwo45XSMI4HjZOHhU5pHYt2snj5Yl+DMWLO4WU1yJgEuOk7Hk0d40p5WnSpUyULKvFZOUr8ZDGydW/IXwpBaf8/sMVs3iRdc+e/Tqb+Xk40PrzkdTbIeJIvy+YiEuX2bGSyba0iu0Yphq7FSb2fezKGtxPIZ7ttOP7OB0SJX7/0rhc9J/yat8s5zmoeLjtPRZHVsaXfeuX2b2RA7yOStfPX0QSgrI/aXQ7Rfp5asbxnu217XOfXUz7xGjIk2EPfLh9hmrFlMb4w01VPopofNJfX6l84e1/2zaLffxqkfFx2nY6n11ju3b3NJR5k1L1FeNrCs/YSnXvHYsWAOq6evKlogYyWJJjDIrFHrbAa2zi0K/rrhTcWyZpMeamun38jJj8/pOB1JtRXw42H19EF6l92eObfQaoZW9LFueFNxy0NPoZveZbc3bTJ/YOvcEnHJEpwZaxbT39N4KwdGC02136gdfsN2QdKxMa/YFklnZhx/qqTvStok6TZJA1nXaQRu6TgdR7LCfbAwq6RTqbZCvh5WTx9kABhkFj0ttHrKLbctC8/LfW4zQswkw5TlFiOUWl9pwRkszGq6m3T6O6r22yTtSJ5jqrpvxzxiK4FjCLEnr5d0uZml/5MsATab2eskPRO4Q9LXzOwvjW6Pi47T0ZQvWhwslA4FDWydO6b1Jcm8z0DfXAbnVw/PXy/lIWQSJ4YNK1eFriGDeUsWhcWWkaEVfRWFdbAwq6JzRb1UE5yeQneJ4Az3bS/O8bQT9eb3mYQcBWwxszsBJF1MyDeW/gdkwF4xHuVTgD8Qwoc1HBcdp2PJ6nR7lxag7K0/EZ4SR4KcjFd8yj3oklD/JW3vC0KRxYw1i+ldWigKTiI2PYx+9kRsemiMQA7PV9EJA8L3OLR8JsMrVLTC5i1Z5HMrTeDuvzyljjmyX+0j6YZUwaoYIT9hP+A3qf1tMCr09xcJQZKHCfnKTjSzJ+psdi5cdJyOZbhve01ngqHlM+lddvuYLZ6EtPjkcTQo96AbjoKQvHUPpVb7D/dtL2bNTJNM1qetmkpi02iygnSunj4IK0tFqFMFp6fQ3ZKhwBZxv5kdWeV4nhxi/cDPgVcSAlBcKWnQzMbnApmBi47TcXSt3VgMIJnV6aXFqItQd2hFHwOMLZRLmtXTBxlYBizrLhnGSuePSUSiN56T1bmlnwGChZZ2MU7WvMwtBKtmdZnFkaZkCGxhON7oIcEsxuuS7bSMSrnF0gwAnzAzA7ZIugt4HnBdoxvjouNMSso72571VlyLM16KHXy0fIAwDLUQ5oVEs0CwsgB619b3Nt1T6Gbd9OpOBNXEM3GEoI+i+CTUY6H1Vq3Z2QwtnwnzJ7oVLeN64NCYV+y3hHxkbymrs5UQrX9Q0rOB5wJ3NqMxNUVH0huB5cCzCGaaADMzz8jkdAxFyyIzR+zYWT19MLgnRxFIhs+GoWTiPy9DK/oyBSc9PFhNcPp7ZhcDdZaITzIElxKhrrUbSzzUEieMySw2UxEz2ynpXcA6YBpwfsw39o54/FzgX4ELJN1C6OOXmdn9zWiPgjVVpYK0BXidmd3ejAZ0Ct3a2+YoM22P00HkXe9SL41YlzK0oq84QZ9M2gMljgd5UgokrBveNGooLhGfRHgScRovA1vn+nBbiqvs0htrzLNUZZ+Z+9hrLjwhV92L5pw/rnu1mjzDa/dOdcFxnFaQnrgf7ttetJTmLqs8qFCts5+3ZFHRWSERlmTF3yBBfIYLs4pW2nhJErY5TjXyiM4Nkr4JfBt4NCk0s8ua1SjHaSTJEFI9iy3r9XZrRIc73LedeQvCnFB6aC4RlR0L5hS9x0ZclStbF11rNzLEnLD+J1IyH5WKFN0IBwumDzYl3lozSRw2XCxbR57htdUZxWZmZzSnSe2JD691Jukhq2aTrKlpd5oxxNiJw2vpJHVDy2dWdLIYS+I8H16rTE1Lx8yaFoPHcZrJeAVn3pJgddSa90jinG1ZeB7z1rf/Ysn+ntkNEePy+aJGRLZ2Jj95vNf2B74AvIywoOgnwHvNbFuT2+Y4E0oiHsNrYaCQHc2gv2d2MfRL/9LZY/JYmwiSdUH1Wj1pB4csOuX5YWT4MbF2qtVzGkeeOZ3VwNeBv4n7b41lxzSrUY7TCHqXFmDhyH5ikZSH5M+iPNZZsuB0HjOLlk95nU4k8XbLyjKaPVS4veOfOU3iSl/pmZL5wGYvtJ1K5BGdZ5pZel7nAklLm9Qex2koaaEpD0bZ35P9pt/fM2KxpCMUpy2fGSsW173os53phLmoVuNpEZpDnnw690t6q6RpcXsr8PtmN8xxGkHv0kJJhzq0oq/4Vl/J0lk3vKnY4WxYuSpz7sM76anBqOCszrjJY+mcQYhA+lnCnM61scxxOoaeQnecj4lCs7Bq9VFpBnwNytRjvNElnGzyeK9tBV7fgrY4TlMYEZyxs2XhecygdE7IcZz6qSg6kj5kZp+U9AVGh8HGzN7T1JY5ToMYWj6TGfNnjdtFuHj+wsaEvXGcqUg1SyfxI7yhSh3HaXuSYZL0avnxLpBcN7zJE5g5zhio6EhgZt+NH/9kZhemN+BPrWme4zSH/p7Z9PfMLnq3jYUNK1fVFK+eQjfrhjfRU/Cg7I4D+bzXzspZ5jgdxY4FcxoSImfd8KYSj7csGhHJ2XEmA9XmdI4Djgf2k/T51KFuYGezG+Y4zSYdCLNh10t5vM1bsojV00fukXZo8KE5Z6pSzdIZJszn/Bm4MbVdTsin7ThOFcpFLW3tJENzOxbM8UWIzpSioqVjZpuATZLWAo+Y2eMAkqYBu7eofUUkHQBcBDwHeAJYZWafk7Q38E3gIOBuYKGZPRDPOQv4W+Bx4D1mti6WHwFcAHQBVxBiyVUPt+1MKpqVzK1eisK00q0fZ2qQZ07nh4TOOaELuKo5zanKTuD9ZjaTkO1+iaRZwJnA1WZ2KHB13CceOwk4HDgW+FIUTIBzgEXAoXE7tpUP4kwsjRCcJAK14zj1kUd09jCzh5Od+PnJzWtSNmZ2j5n9LH5+iODSvR9wAnBhrHYh8Ib4+QTgYjN71MzuArYAR0naF+g2s59G6+ai1DnOJKdRXmSNng8Cj2bsTA3yiM4jkl6c7MShqR3Na1JtJB0EvAjYCDzbzO6BIEzAs2K1/YDfpE7bFsv2i5/Ly7Pus0jSDZJueGwkaarToSSCU54Hpl3IivTsOJONPLHXlgKXSBqO+/sCJzatRTWQ9BTgW8BSM9suqWLVjDKrUj660GwVsApC5tD6W+tMBMnwWToN84w1i+lhc1u7Lm9ZeF4xJlySt8atH2eykSf22vWSngc8l9Bh/8LMHiuvJ6lWykAB95jZYWNqabjHbgTB+ZqZXRaL75W0r5ndE4fO7ovl24ADUqfvT/DI2xY/l5c7HcyOBXOqeou1KmV1o1g9fRBWDsJKD7njTC7yZA59MvA+4EAz+ztJh0p6rpl9r6zqkJm9qMa1bhprQxVMmq8At5vZZ1KHLgdOAz4R/34nVf51SZ8BeggOA9eZ2eOSHpLURxieO5WQGdXpMEoDebaHN5rjONXJmzn0RuClcX8bcAlQLjpvynGtPHUq8TLgbcAtkn4ey/6RIDZrJP0tsJWY4dTMbpO0BthM8Hxbkrh9A3/PiMv09+PmdACNiBjdafQUuhnuqzWQ4DidQR7R6TWzEyWdDGBmO5QxkWJmd6b3JXWnr29mfyivUw9m9hOy52MAjq5wztnA2RnlNwDPH2tbnNaQnpuZygwWZo3Keuo4nUoe77W/SOoiTrZL6oXKrlySFku6F7iZkSgGHqnaaQhTzcqBMB/lgUOdyUIeS+ejwA+AAyR9jTDMdXqV+h8ADjez+8ffPMcJTEWxSTPVn9+ZPNS0dMzsSuCNBKH5BnCkmV1T5ZQhPPWBMw6GVvSVRA3wDrcybv04nUa1KNPPM7NfpBaG3hP/To9x0P5gZr/OOPUs4FpJG0kNw3mmUccZP6unDzJvQYjRlgjz6umDMOyu1U5nUG147X2E+GT/WeH4MyRtMrO3lZWfB/wIuIUQmNNxnAZSnkIhYWhFH71L3eHAaW+qRZleFP++olIdST/MKN5pZu9rQNucKcbIAk9fc5PFjDWLqy5yTSIaeLRqp53Jszh0N8K6lnmx6BrgPDN7zMxenXHKjyUtAr5L6fDaH8bfXMeZuuSNqpC2hGasWezWj9NW5PFeOwfYDfhS3H9bLHt7hfpviX/TKa0NOGQsDXQmPyNOAyMWTjpumuM4k4c8ovMSM0vPUP5IUrXxj0PKE6JJ2mNMrXMmPZVy27jgBPp7Zo8r/49bOU67kWdx6ONxQSgAkg4hZOKsxFfSO5L2BP5nbM1znKnNuuFN9PfMHpNnmnuzOQmSjpV0h6Qtks6sUOflkn4u6TZJ65vVljyWzgcI8zR3EsLQHAgMVKn/W0nnmNnfS3o6QXD+a/xNdSYTWVGhnWyGVvTVFSXbxcZJEzMmrwSOIcTOvF7S5Wa2OVXnaYQplGPNbKukZ2VerAFUFZ3Y2NmECM3p1AYVw+CY2UckLZd0LnAE8Akz+1YD2+x0OGG4yD3U8tJpaRmctuMoYEsS+1LSxYTMyptTdd4CXGZmWwHM7L5RV2kQVUUnpgF4vZl9lhBLrSKS3pjavQ74SPxrkt6Yyn/jTFHyzE0MbJ1bElF5PPMZU5UdC+a4y7STJiuL8pyyOocBu0m6BtgL+JyZXdSMxuQZXrtW0heBbwKPJIVm9rOyeq8r27+J4PX2OoL3movOFCMRjHrWjayePkg/I8ND451In4qkXaZ9qK0zeeiRLgYLs/JW30dSOqjyqpj1OCFPtuRdCSNTRxNSvvxUUsHMfpm3EXnJIzr/L/79eKrMgFemK5lZtXkeZ4oxtKKPZAitfAV9vULiwjN2PErBlOB+MzuyyvFKWZTL69xvZo8Aj0jaQJhaabjo5An4+YqM7ZXl9eKC0KrkqeNMDirNQ/T3zI6CNLKfMGPNYiBYSDsWzBl1nlM/Ph/kANcDh0o6WNKTgJMImZXTfAeYK2nXmC16DnB7MxqTJyLBMwjpDf6aYOH8BPi4mf2+rOqZkqqlMxDwXsBdlqYY5cNrSUcYkrNtLwpKLwX6l5ZaNTsWzKF32e2+bsdxxoiZ7ZT0LmAdMA04P2ZWfkc8fq6Z3S7pB4S5+yeAL5vZrc1oT57htYuBDYykmj6FML/zqrJ66xk9r1POlXW1zulIykPtV5rPyZuC2QVn7MxYs9izjjqY2RXAFWVl55btfwr4VLPbkkd09jazf03t/5ukN5RX8jmdqU29czSVrjGwdS4DW+em5oGyr5uksHZBqs6WhefRv9SHJp32IU9Egh9LOknSLnFbiEcYcCowb8n4pu1WTx/MJSSDhVn1ePdMabLmyBxnPEh6pqRPS7pC0o+SLc+5eURnMfB1QsToRwnDbe+T9JCkfOMjzqQnCdXStXZj0QqphzyWUrmguVdWfnzdjtNgvkZwNDgY+BfgboLDQk3yeK/tZWa7mNlucdsllu1lZp4n18kkLTyV3rKHVvTRU+iumm45EZpE0BIPty0Lz3M36jroKXS7teM0kmeY2VeAx8xsvZmdAfTVOgnyzenUhaQjgblAD7ADuBW4yvPpTA3KhaC/ZzZdZL9lz+3bnGsobd6SRfQWEg82F5pKlEdz6Cl0F7/f1dMHYeUgA8vm5nbgcJwqPBb/3iPpNYR1P/vnObFhoiPpdOA9wF3AjcAdwB4EV+tlkm4FPpLE9nEmH+WCk+U5lVg1w33bGe7bXhJ9IMty8aCg+SkXk+T7TQdXLY/44Dhj5N8kPRV4P/AFoBv4hzwnVhQdSQeb2V11NGJP4GVmtqPC9V5ICBzqojNFSKdPHp6v0oWK5euhyU7HPGPNYnrWG73Lwjo191arn661G+lfG4Rmx4I5FS1Px8mLmX0vfnwQeEU951ab07kUQNLVORuxspLgxOM/N7Nc13ImFxtWrsq1Mj6rzpaF5xUXh7rgVCfPHJc7FDiNQNKFMR1Csv90SefnObfa8Noukj4KHCbpfeUHzewzFRrzSeDfCPM5PyDE71lqZv+dp0FOZ9LMSX0Xm3zMW7LIrRinVbzAzP6Y7JjZA5JelOfEapbOScCfCcK0V8ZWiVeb2XbgtYQgcocBH8zTGGfykMdtet6SRUVvNGfszFuyqK5I3o7TAHaJSToBkLQ3OX0EKlYyszuA5ZJuNrPv19GY3eLf44FvmNkfpKzI2k4rSbslJ4sqG7XOJR1RGpKIA9sz523SdbrYSO9amMFiepcWKlpLHmV6NPOWLCo6ByR/B5bNZbAwy9cvOa3gPwlpby6N+38DnJ3nxLz5dD4DzIv76wkBPx+sUP+7kn5BGF57p6RnEiwmZwIpGaKaPhgskUI3g4VZ9Ky3cb0ll8/F5BGIdcObUo4Dm5jBYsrdoRMngmpDRlnOB5OdGWsWM3fZZga2zi35XVdPH4Tpg0XnDbd8nGZhZhdJupHgRCDgjen019WQWXkun7IK0rcIa20ujEVvA2ab2RurnPN0YHvMPLonsJeZ/S5Pg9qVbu1tc3T0RDejLpIUAtU65fQw2FjXbySLDrPcm5POz6NFN4a8MefK1+w4reUqu/TGGjluqrL79ANsvw8szVX3rvd+YFz3GiuSpgHPJmW85FkSkycMTq+ZfdTM7ozbvwCHVGnIGwnqd0L83A8cLulZOe7VMiQdK+kOSVsknTnR7Wkk64Y3sW54E1sWnlfTCkh7haWH4HYsmFOS96YaXWs30rV2YzF6QFrIkrftrrUbGe7bXnJsxprF9PfMLs5JpBlvDLfJyLwlixju2+4x55wJR9K7gXsJmQO+R4jH+b2qJyXn5rB0fgp80Mx+EvdfBnzazF5aof7/AC8FfhyLXg4UCA4FHzezr+ZpWDOJCv1L4BiCs8P1wMnVzMNOsnTGOv+RCMLQ8pnFdTXjSZ6WtKN8GChN+vrjnbdJh8iZTKS/o6EVfbmfz62diWOyWzqStgBzMvKq1STPnM47gIvi6lOAB4DTqtR/AphpZvfGxj0bOIeQiW4DMOGiAxwFbDGzOwEkXQycAOQak5ysJMIwsAx6GZsLbtaQXrWhoERogmUzIlJDy2cWh+vSk+blpAUtmUAP80OTR3xKxTi/MHv0AaeJ/IawMLRuaoqOmW0CZkvqjvu1Xp0OSgQnch9wWPRie6zSSS1mP8KXlrCNIIolxPTaiwD24Mmtadk4CfMr47MYiiKxcpAZ8xcXO/MknMqMNYszPaTqeQsvJy0qg4VZ9K4tMI8gNhtWrsq0lpJ2pDvWnkI366ZPDrFpBEMr+tybzWkGdwLXxJGtR5PCSus30+SOvZZDbBIGJX0PuCTuvwnYEB0K/pj3fk0my4d71Dijma0iptfu1t7VxyHbhCRcTKOY27e56Plc7dqJhVPLvTktWDsWzBkdHid1z/T9sqylLQvPYwbBw638Oo0ctutkku/IhcdpMFvj9qS45abmnE69KCzKeRPwMkLn/hPgW9boG40DSS8FPmZm/XH/LAAz+49K53TCnM54LI1aJPMl5e7V6UjG1Ugslf6e2ZkOClmCkRaLpCx9v3KBS7tPl89FTWXhKWc883ROPib7nE6CpL0AM7OH857T8NQGUVwujVu7cj1wqKSDgd8Soi+8ZWKbND5Cp9q8jrUoCgspDnvNW7KIoeXAynzuu/3MHnM71w1vCoK1fCYz5sfFrTGCdXL9XgoM9AVxS+onJGI3WeZ5HGcikfR8wvz83nH/fuBUM7ut1rk1RUfSkwnhq6eb2d9JOhR4birKaHn9NwLLgWcRLB0RtKhtEr6Z2U5J7wLWAdOA8/N8WVONga1zM91zt6wc6biTCMaJ9ZLOkZNYHkPLZxYdEpIFn8PzVRzyKbeWEsHIGqrrWhuiGCSUL4Ic7qscCaF3aQEW1vEFTCF87sepk1XA+8zsxwCSXg78F/D/ap2Yx9JZTciPk7hIbyPM11Tyyf4k8Doza+zkQoMxsyuAKya6HY2gWUNH5SFVysUh3dkn9YahOLG/ZXi0VZHUSwvHcN92BgrZbtXlDgTp7JfJ+qCkXYllU234yEPqlFLibbjQh96c3OyZCA6AmV0T5+1rkkd0es3sREknx4vvUPVgave2u+BMJhrhrVaJLQvPY976EUuiXBR6Ct0l60ASQUgiECTtqhSOJenwst6ws4QhiMuIl9vAsuBanZT523o+slywkzm78qR7jlOBOyV9hJElMG8lJPCsSR7R+YukLqJ3l6ReUi5yGdwg6ZvAtyl1pbssT4Oc9mLDylX0r822DlZPH4RhSmKoAbAS0kK4Y8GcUcKzY8Gc4vzKDBbTk2OJVHEdUSqqQdrVOnlbT+oMLZ9ZvG9WGzqdvCFxatHfM9vFxqmXM4B/AZJ+fQMwkOfEPKLzMUJenAMkfY3glVbt4t3An4BXp8os1Ting5i3ZBEbhqunjK41Od+77HaG11atkrvjLE+ZUC3aQe+y22FZd7z+qiiGnU3y/Om5tgHGLjzJAuDJKMpO8zCzB4D3jOXcmrHXzOyHwBuB04FvAEemx/Iy6g9kbGeMpXHOxFMpEkA9rJ4+WBLXDaIDQs/suqJEp9/sk21o+czMnDzJsU6NU1YpH1FaXHrWGz3rbVzPuGHlKoZW9LFh5apRv5HjVELSlRmZQ9flOTeP99rVZnY0IaBbeVm63ofM7JOSvkD2QssxqaKTTTJ/Mjy/M3IVpT3YEup1n856m0+82fqXzi5xdJixZjG9awv0ri11boD8a4valdXTBxkABplVFJ7EVbxeEtHPijKRRKBISK+1cqY8+2RkDs0V1Lmi6EjaA3gysE9MVZD0bt1AT8YpifPADXlu7NTPWGNwTTRJwrY0Y32rTtbkJMKRHhZK1gJB9Qnxam7V7UItASkXnsHCrJBLp0561lvFiAXF5HDRdb4oTj7/48ATkqYnqQwkHUiGsZFFNUtnMbCUIDA3MiI628kYHTez78a/F5Yfc8ZH+RvnRFFt/gRG1syko0sPFmZldlJjtTSSIJaJwNQbkDQhcZ1u1GT8RDI8X/Sst7rWICW/ZfHfVTy3fN1TeuFt/1J3OHCKfBj4iaT1cX8eMU5lLaqlq/4c8DlJ7zazL+RtiaTDgA8AB1Ga3OeVea/hjFDuJjxRjI4qMJoNK0sn64eWz6R3bXYnVUvAqtEo1+i0h1t5NOaJWMtTz3dSaeFuHsrTiY/yBlxZ5hmXquepEhwAM/uBpBcDfQSD5B/M7P485+aJMv2FGPJgFrBHqvyiCqdcApwLfBl4PE8jnGxaNfeQ9ojqXVqo2OH2FLoZWj6TvEN7lbyhxvtcjQpiOdy3vaKl1Oo02JUcB6rVTSI7QPyuc3rnJeurRuZmRoRkxorglDG3b7QLeyJAPqczdZH0PDP7RRQcGBmonh6H235W6xp5HAk+SkjENouwgv84QhDPSqKz08zOqXVdpzFUyzWTVTcRgmROJby5hk6n2tDJ6umDDGydG9btpFb1V3o777SOKT2EmYTqaWXInLxBUxPrpmf9yPB58jmvpZSIR/lQGmRHloBSy698UbAzpXgfYRjtP1Nl6bmcmiNaedbpvBmYDdxkZgMxKduXyytJ2jt+/K6kdwJrKV0c+occ93JSVJrwTg9zdLGRgWX5OptRMcrqJH2P9Fty0jnl6YzaNQRN19rwPULpItN2IRkKTIYrk7QQaeq1HjesXMW8OAxfa41O+kXDk8NNXcwsmbc5B/iBmW2PkQleDPxrnmvUXKcD7DCzJ4CdMZHbfcAhGfVuJHiunQZ8ELg2liXlzhgY2DqXeUsWMbB1LgNb59LfM3tUxz7ctz1m3sxmxprFDbM8qt0nS3B6Ct2sG97EjgVzGio4zQh3M9y3vWhJ1DPc1QoGC7MqCkOSd2gsbU6S5K0b3jRqK09BkbUeypmy/FMUnL8GjgEuIAhRTfJYOjfERUD/RRCQh4HryiuZ2cEQXK3N7M/pY9H92hkDybxDrRX9SbRnCB19OlhnvR5HWUExi55oFRwDKpG8fTfS+66ZQ3e9S4OX1rrhxs2ljddDrpaYJB6CjXYF37LwPPqXjnzXyTxPaI8Pr01xkvn61wDnmtl3JH0sz4l5HAneGT+eK+kHQLeZ3VzllGsJplatMqdJJC6uY6Wi4LRBMM1WdXhpYRuL40OlcDUJta5XPn9TLvZdazfSw+jQNel5nfF4CGYJ+2BhFgOMbWjWmVgkHQt8jpDK5ctm9okK9V4CFIATzaxaTrTfSjoPeBWwXNLu5Bs5y5fETdJ+wIFJfUnzzGxDWZ3nAPsBXZJeROli0ifnuY8z8WQJDlQfzkqfkzU53Uha2eGNZX1UWmx61hs9jHiYDRZmleQbqnaNau7mCVnfc5ImIn2teoSnmhXZu7TQ7mtqnQwkTSP4Nh5DSE1zvaTLzWxzRr3lhDxjtVgIHAt82sz+KGlfwrRKTfJ4ry0HTgQ2M2JSGSGqaJp+Qny2/QmeDenFpP+YpzFOe5F+2+6lMCpvDTTXMSA9LNXqNSL1huhJD4ElgpHO/dO7NO6PztQ96jpjGcZMM7R8JhtWriq2KY/wdJq3oVMXRwFbzOxOAEkXAyfAqNDu7wa+Bbyk1gXN7E+kgjib2T3APXkak8fSeQMhU2i1dAaY2YWSvgqcbGZfy3Nzp/2Z27e5KDxJ519peGtg69zi2/d4xSgtOFlxwZpJPW0vHwbrWruxZO3P8HyFuZC4yDfvZP9Yoz4nC2dnzF+cy6oCF5x2ZLeHSt3iq3FXCFWWdtZaZWZpE30/4Dep/W3AnNR+Mpq1gODyXFN0xkMe0bkT2I3qOXQAMLMnJC0GXHQ6iKEVfZkLIROR6V0arJzyDrNaArmxLq5M32OwMIv+cc5P1Uvw2MonOrWGwZJ5l7R1U00EGuExV3QeWVpgaMEcWFlbdBKRdfHpWO43syOrHM+KClyuaCuAZWb2ePUcneMnj+j8Cfi5pKspXXdTKWr0lZI+AHwTeCRV39fptBFDK/qK3kjrplcXh6EVfaybfl7JEFd5eJ5GWCOJ1ZC0a27f5pbOIdQ7pJZnGKxrbfA8rBbgtNxa6sHGPS+WWEqJRyPUdojwRZ+Tlm3AAan9/Rnt53gkcHEUnH2A4yXtNLNvN7oxeUTn8rjlJcmdsyRVZmSv7XEmiN6lBQYrWDhpVk8fLIleXKnjGr2gsnrnXW7RJCSCMzKU157kaVt5Ou6suZVErJthzSULSLNC2pRTLTir0/FcDxwq6WDgt8BJwFvSFZIlLwCSLgC+1wzBgXwu03VFjU433mlv6h3+alQcuKSDS+ZAeimUWF5JnVZ7qjVqWC2hfEx+sDCLectnFoWgmW7oxXBH1HZegPG72Tvti5ntlPQuglfaNOB8M7tN0jvi8XNb2Z5q+XTWmNlCSbeQnZTtBRXO2w34e0Koa4BrgPPM7LHxN9fpZMoDi6bpXVqAuKg1EcNWhlqpJ2BmEpk61zUjw33b6VkQ/hv1rDeG1s+kJ1/6kXHRtXYjQ8yhd9ntFV8afC5n8mNmVxBiZ6bLMsXGzE5vZluqWTrvjX9fW+c1zyE4Hnwp7r8tlr29zus4TaaeYKH1kB56Snd0iXWTZSHsWDCHwcLIUNC8JYvGnCunXur1VhsraRFKu1M3kx0LqguO47Saavl07ol/f13nNV9iZulXpx9Jas8oj1Oc8olmGOmkIN+q+XISYYHgLjwjNV+TDKclDK3oKwnPnwhOVqbRdmGwMKshFkozF9CmGZ4vNrjgOG1EteG1h6iSftTMKrnjPC6p18yG4nUOwfPqdAxZlk/5pH8iKlkdZ3peoDcjXlx6/cncvs3QB8NRmCbKcyor1lwWA1vnjhLOdqeaE4EPqzkTQTVLZy8ASR8Hfgd8leDvfQqwV5VrfhD4saQ7Y/0DKQ075aRI1si0Q1bGrI43PekPjGulfGJFlQcvneiOPK/g5J3L6QQ8aKczUeRxme43s/QA9DmSNgKfzKpsZldLOhR4LkF0flErmsFUJolq3FMoXcvRygCblRaHlnuZNZp0TqB2plMFJ8wb3T7RzXCcEvJEBX1c0imSpknaRdIp1B4uOwJ4PiH524mSTh1vQyc7w33bSyyduX2b6Sl001Pobtmkc5pmCE7X2o1tYc3tWDCnmOenlpUzY83ijhWcZG4ui4n+HZypSx7ReQth2d+9cfsbyhYWpYnx1z4N/DUhhs9LCKtdnRyUiw+EJF3NEp6eQneJlTOwdW6xo+1dWui4zjaLRGASNqxcVbPTTb6HiR76Gw+rpw9WXGRbnqDNcVpF1eG1GOp6iZmdUMc1jwRmmVnzFyFMYtLhZiB4IWVNzI+VSpEFGm3dFOdx2uDNOi081SycZg8rtoLyVNZ5nSUcp9lUFZ0Y/O2IOq95K/Accoa5dqoz3Lc9DAe1YCFho72zyqMMTBTDfdvpJ7+HWrskrBsPPeuNgb5UxtLhkTA87RxeyJn85HEkuEnS5cAllAbwvKxC/X2AzZKuozRA6OvH2khJnwJeB/wFGAIGzOyP8dhZwN8S5pneY2brYvkRhLzdXYSVuO81M4sZ7i4izDv9npAh7+6xtq2ZJGHqm/G2nRUDrNFv9p2W9GuyCE4lkt88Twgfx2kWeURnb0Ln/MpUmZFK4FPGx8bZpiyuBM6KMYSWA2cByyTNIgSvOxzoAa6SdJiZPU6IgrCIkHr1CkKWu+8TBOoBM5sh6SRCprwTm9DmcVMeZbl8XmesApGVHbSTh5KySIfrXze8iXlLFtHfU31YrRMEp948Ox6JwGk3ajoSmNlAxnZGeT3FmNhmtj5rS9epFzP7oZntjLsFQmhuCNnvLjazR83sLmALcFRMndptZj+Nc0sXEZLRJeckQUwvBY4ea7uaReJdBSMTvknEYCcfycLHvPMYiaU3WUg8BectWVQsG9g6N0R7mEQvF07nUVN0JO0vaa2k+yTdK+lbkvbPqPpjSe+WNL3s/CdJeqWkC4HTGtDmMwgWC2RnxNsvbtsyykvOiUL2IPCMBrSrLnYsmFPRgyhxdR1aPrNk9X/68/B85fZoS1yvs0jncpnMbFi5quLi13lLFnWMp169bexau7Ek0KrjTDR5XKZXE/Lp9BA67O/GsnKOJcyrfEPSsKTNMSrBr4CTgc+a2QWVbiLpKkm3ZmwnpOp8GNjJSGbSShnxqmXKy5NFL7nfIkk3SLrhsdqJU8fM0Iq+ooCUi0O6k+lau3EkZH0db+WJ6+y64U0lGUA73UurWmI0CNZO+k2/nE5dg1MvPsTmtBN55nSeaWZpkblA0tLySmb2Z0Jk6S/F9Ab7ADuSCf9amNmrqh2XdBoh4vXRKXfsShnxtjEyBJcuT5+zTdKuwFOBzKymMc/4KoBu7d3QsZeutRtLXKCT8DBJBzFvyaKanWHezjI91NSu3kuJp1vSvuG+7awb3jQqPtio6AnxV00/03Df9pRLePbw2mScx0qTCHJacFqdidVxssgjOvdLeivwjbh/MsGxoCIxd07DXKYlHQssA+ab2Z9Shy4Hvi7pMwRL7FDguujq/ZCkPmAjcCrwhdQ5pwE/Bd4M/Gii1xQluWRqvbmnGUtnme6Y20V4knU8PQTX6qRdyXexbnhTsayaiJS8zefoWfOkmu4UsgQmi8R1Os8LjeM0izyicwbwReCzcf9/GUlJ3Sq+COwOXBnn/Atm9o6Y/W4NsJkw7LYkeq5BSCR3AcFl+vuMzAN9BfiqpC0EC+eklj1FFYqLJ8f5KlrNuyndKZVYBvNFD/V5RdVqw/B8VfQES6/fGSyI3rLjg4VZJet7sjrTrNTPWdRKvzAZyDt8lnwXLjjORJInXfVWYMxrbBqBmc2ocuxs4OyM8hsI8d/Ky/9MCOXTlvT3zA6d8rLNo6Ix5yEREBjduSTXLg/u2egOOEnHPJQhgOX3HyA4TJS0B2OQWcwoE5+EpJMttYBGUx4hO2GyWDgJtYYgq2VsdZxWU1N0oqfaF4CXESbcf0JYaLmt6onOuEksgnomu3uXFopOCcnfrrUbi55bA1vDUEy1LJ7jJUnHnNXm9Er5ak4MyXzXMKXrk9KJ4RIxmTF/tDglnlpTsZPtXVoI0RIjRZFm3Ia044wb1ZrOkHQl8HVCPh2AtwKnmNkxTW5bW9GtvW2Ojm7Z/RJ36rFMdpd30lCazKsVHmv1LmIcz7WzXMd9CGl0fD2fy2kdV9mlN5rZmAMdP+XpB9gLX/neXHX/97IPjuteraZh3mtOY0ne0HcsmJPZyVbqPNLH0p1xMk/SKhfpRs8TpSm/pnek2SQx5xLaPW+RMzVoivea0zjq7VDTnX3aLXvHgjkMrZ9JD61xE56Kw1rNpBihYgqsK3ImN/V6rxlwLa33XnNSVOt0KnX23lF1JjsWzGHDylUjBSsH6V87u/IJjtPmdIT3mlMfzZxPcVpLieA4ziQgT+y1CyU9LbX/dEnnN7VVTm58En1q0S6Leh1nrOSJvfaCdCgbM3sAeFHTWuTkJlmE2axU1s7E4yLjTDbyzOnsIunpUWyQtHfO85wm07V2Y3EhaKNy7TjtxdDymbByxO25HdJ+O854yCMe/wlcK+lSgiPBQjIiADgTQ9fajcVFpFN5QeRkJUlz4TiThTyOBBdJuoGQOVTAG81s4hPfO0WStNA9CyZv1OSpSN4EdI7TSeSZ08HMNpvZF83sCy447YsLzuQha56unvmdJKJFpUSBjjNR+NyM47QpM9YsLuYYysuIdbQpxl/bRP9SX9fjtA8uOo7ThiTRJIaBeQtS2U8XZFu0WdHDYXQEaseZaPJEmX4X8LXEe81xnNaSiEziMFKe8mLHgjmZqQzc081pR/JYOs8Brpf0M+B8YN1EZ9p0nMlOOi1Fsl8cakvlJ0iG4GBEbJL8RB7g02lHajoSmNk/EdJAfwU4HfiVpH+XVJ7w0XGcBtG1dmMI3lroLgpQ1txOWnAGC7MY7tteDPbqOAmSjpV0h6Qtks7MOH6KpJvjdq2kpo3L5vVeM+B3cdsJPB24VNInm9Uwx5nq9Kw3hpbPZHi+QvK6NYuBIDADW+cW9xOx8fVZThaSpgErgeOAWcDJkmaVVbsLmG9mLwD+FWha0L88czrvAU4D7ge+DHzQzB6TtAvwK+BDzWqc40xlEmuld+3IcNu89SmngvnR+pk+WMwU6o4DTgZHAVvM7E4ASRcDJwDF5S9mdm2qfgHYv1mNyTOnsw9hQeiv04Vm9oSk1zanWY7jpMkcLpvva3CcXOwH/Ca1vw2oFrDxb4HvN6sxeSIS/HOVYx6jw3GaSLU0Fb1LC0ULJ6Gn0O1ea5OAXf74SD3zcvvEqDEJq8wsPTymjHMyncEkvYIgOn+d9+b14ut0HKeN6V12O4Pz+yrO1yTzOolDwWBhFr343M4U434zO7LK8W3AAan9/SnxgQxIegFhCuU4M2tadmgXHcdpY4b7ttOzwIppLMrFp2d9iLeX9CAuOE4G1wOHSjoY+C1wEvCWdAVJ04HLgLeZ2S+b2RgXHcfpAIbnZ42QOE5tzGxnXOS/DpgGnG9mt0l6Rzx+LvDPwDOAL0kC2FnDehozLjqO08bsWDCH3mW30+XzNM44MLMrgCvKys5NfX478PZWtCXXOh3HcSaGrrUbi3mSKh13nE7CRcdx2pxqiz53LJhDT6GbnkJ3C1vkOGPHh9ccp8MpZox1JwKnA3DRcZwOJUnQ5uFvnE7Ch9ccp0PIyibqOJ2GWzqO0yGUOw24heN0Im7pOI7jOC2jo0RH0gckmaR9UmVnxRwRd0jqT5UfIemWeOzziiueJO0u6ZuxfKOkgybgURzHcaYkHSM6kg4AjgG2pspmEUI6HA4cS1hNOy0ePgdYREhAd2g8DiGY3QNmNgP4LLC8JQ/gOI7jdI7oEATiQ5RGRz0BuNjMHjWzu4AtwFGS9gW6zeynMQHdRcAbUudcGD9fChydWEGO4zhOc+kI0ZH0euC3Zrap7FBWnoj94rYto7zkHDPbCTxIiDmUdd9Fkm6QdMNjPDru53Acx5nqtI33mqSrgOdkHPow8I/Aq7NOyyizKuXVzhldGHJSrALo1t6ZdRzHcZz8tI3omNmrssol/RVwMLApjoLtD/xM0lFUzhOxjdJ0q+n8Eck52yTtCjwV+EPjnsRxHMepRNsPr5nZLWb2LDM7yMwOIojGi83sd8DlwEnRI+1ggsPAdWZ2D/CQpL44X3Mq8J14ycuB0+LnNwM/ivM+juM4TpNpG0tnLMScEGuAzcBOYImZPR4P/z1wAdBFyPed5Pz+CvBVSVsIFs5JLW204zjOFKbjRCdaO+n9s4GzM+rdADw/o/zPwN80q32O4zhOZdp+eM1xHMeZPLjoOI7jOC3DRcdxHMdpGS46juM4Tstw0XEcx3FahouO4ziO0zJcdBzHcZyW4aLjOI7jtAwXHcdxHKdluOg4juM4LcNFx3Ecx2kZLjqO4zhOy3DRcRzHcVqGi47jOI7TMlx0HMdxnJbhouM4juO0DBcdx3Ecp2W46DiO4zgtw0XHcRzHaRkuOo7jOJMcScdKukPSFklnZhyXpM/H4zdLenGz2uKi4ziOM4mRNA1YCRwHzAJOljSrrNpxwKFxWwSc06z2uOg4juNMbo4CtpjZnWb2F+Bi4ISyOicAF1mgADxN0r7NaMyuzbjoZOQhHrj/Krv01xPdDmAf4P6JbkQD8OdoHybDM0B7PceB4zn5IR5Yd5Vduk/O6ntIuiG1v8rMVqX29wN+k9rfBswpu0ZWnf2Ae3K2ITcuOjkxs2dOdBsAJN1gZkdOdDvGiz9H+zAZngEmz3MAmNmxDbycsm4xhjoNwYfXHMdxJjfbgANS+/sDw2Oo0xBcdBzHcSY31wOHSjpY0pOAk4DLy+pcDpwavdj6gAfNrOFDa+DDa53IqtpVOgJ/jvZhMjwDTJ7naChmtlPSu4B1wDTgfDO7TdI74vFzgSuA44EtwJ+AgWa1R2ZNGbZzHMdxnFH48JrjOI7TMlx0HMdxnJbhotOGSPqAJJO0T6rsrBii4g5J/anyIyTdEo99XpJi+e6SvhnLN0o6qIXt/5SkX8RwGmslPa0Tn6MatcKKTCSSDpD0Y0m3S7pN0ntj+d6SrpT0q/j36alz6vpdWvw80yTdJOl7nfwcTsTMfGujjeC2uA74NbBPLJsFbAJ2Bw4GhoBp8dh1wEsJfvbfB46L5e8Ezo2fTwK+2cJneDWwa/y8HFjeic9R5fmmxbYfAjwpPtOsiW5Xqn37Ai+On/cCfhm/+08CZ8byM8fzu7T4ed4HfB34XtzvyOfwLWxu6bQfnwU+ROnCrBOAi83sUTO7i+BhclQMU9FtZj+18D/rIuANqXMujJ8vBY5u1dudmf3QzHbG3QLB57/jnqMKecKKTBhmdo+Z/Sx+fgi4nbC6PP1dXkjpd1zv79ISJO0PvAb4cqq4457DGcFFp42Q9Hrgt2a2qexQpRAV+8XP5eUl50QBeBB4RhOaXYszCG+WJW2KdNJzpKn0HG1HHI58EbAReLbFtRfx77NitbH8Lq1iBeEl7IlUWSc+hxPxdTotRtJVwHMyDn0Y+EfC0NSo0zLKrEp5tXMaQrXnMLPvxDofBnYCX6vRpgl7jjHSjm0ahaSnAN8ClprZ9ioG4lh+l6Yj6bXAfWZ2o6SX5zklo2zCn8MpxUWnxZjZq7LKJf0VYRx6U+wc9gd+JukoKoeo2MbI0FW6nNQ52yTtCjwV+EOznyNB0mnAa4Gj45BGuk3l7Z2w5xgjLQsZMlYk7UYQnK+Z2WWx+F5J+5rZPXHI6b5YPpbfpRW8DHi9pOOBPYBuSf9N5z2Hk2aiJ5V8y96AuxlxJDic0gnSOxmZIL0e6GNkgvT4WL6E0gn4NS1s+7HAZuCZZeUd9RxVnm/X2PaDGXEkOHyi25VqnwjzFivKyj9F6QT8J8f6u0zAM72cEUeCjn0O38xFp123tOjE/Q8TvHHuIOV5AxwJ3BqPfZGRKBN7AJcQJlOvAw5pYdu3EMbWfx63czvxOWo84/EEr7AhwpDihLcp1ba/Jgwf3Zz6DY4nzIVdDfwq/t17rL/LBDxTWnQ69jl8Mw+D4ziO47QO915zHMdxWoaLjuM4jtMyXHQcx3GcluGi4ziO47QMFx3HcRynZbjoOJMGSQ/Hvz2SLh3HdZZKenKD2vQ8ST+PUZJ7G3HN1LW/LGnWGM57YVxwmey/vt0iZTuTF3eZdjoSSbvaSFDRpOxhM3tKA659N3Ckmd3fgGudCXSZ2UfHeP6o52xAm04nPN+7Gnldx8mDWzrOuJD0kpg3Zw9Je8b8Lc/PqHdqrLdJ0ldj2YGSro7lV0uaXqP8AkmfkfRjYLmkgyX9VNL1kv41da+DJN0aP58u6TJJP4j5Vz6ZqneOpBtim/8llr0H6AF+HO+DpFfH+/xM0iUxpln5871QUkEjOYSeHq2JpcDbk2uVnfOwpP+M171a0jNj+TWS/l3SeuC9ko6OltItks6XtHuq3pHV2hh/n2vj936dpKcCHwdOjBbYifE7+mKO7/7z8Vp3SnpzHf9MHGeEiV6d6lvnb8C/AZ8GVgJnZRw/nLBCPAnrs3f8+13gtPj5DODbNcovAL7HSGiTy4FT4+clwMPx80HArfHz6YRwKE8lRDf4NXBAWTumAdcAL4j7d6faug+wAdgz7i8D/jnjGW8G5sfPHyeGoAE+BnygwvdmwCnx8z8DX4yfrwG+FD/vQYjucFjcv4gQwDOpd2SlNhJC9NwJvCSWdxNC+Jye3Cv1HSX3rvbdX0J4UZ1FSO0w4f/2fOu8zS0dpxF8HDiG0AF+MuP4K4FLLQ5XmVkSsPOlhORcAF8lhG+pVg5wiZk9Hj+/DPhGql4lrjazB83sz4SYcAfG8oWSfgbcRBDGrPmRvlj+v5J+DpyWOh+AaD08zczWx6ILgXlV2pPwBPDN+Pm/KX3OpPy5wF1m9ssq167UxucC95jZ9QBmtt1qD9VV++6/bWZPmNlm4Nm1H89xRuNRpp1GsDfwFGA3wpv5I2XHRb5Q8pXqpMvLr53nuo+mPj8O7CrpYOADBCvgAUkXENpejoArzezkHPcZL1nPmSdhXWYbJb2A8YfwT5+f/h4nOpGe06G4peM0glXARwh5c5ZnHL+aYFU8A0KO+1h+LSFyNMApwE9qlJfzv2X16qGb0LE/KOnZwHGpYw8R0jxDyHz6MkkzYtufLOmw9IXM7EHgAUlzY9HbgPXUZhcgmRt5C9nP+QvgoOT+Fa5dqY2/AHokvSSW76WQHiL9fOXk/e4dZ0y4peOMC0mnAjvN7OuSpgHXSnqlmf0oqWNmt0k6G1gv6XHCcNbpwHuA8yV9EPg/YCCeUqm8nPcCX5f0XkLumNyY2SZJNwG3EeY9/jd1eBXwfUn3mNkrFLy9vpFM4AP/RIgwneY04FwFV+s7q7Q5zSPA4ZJuJGREPTGjnX+WNABcEgXjeuDc0ir2f1ltNLNfSjoR+IKkLmAH8Crgx8CZcSjuP8pumfe7d5wx4S7TjjNBaJwu3pJuAV5vZnc1sFmO01R8eM1xOhBJVwK3uOA4nYZbOo7jOE7LcEvH6QgkdUlaH+eNGnXN0yX1pPbHFFamxj2ukvT0Cscejws0e+L+w+O4zz9I2pos8nScdsVFx+kUzgAuS63RaQSnE6IPAGBmb49rUBrJV4F3Vji2w8xeaGbD472JmX2WsCDUcdoaFx2nUzgF+E6yI+mDCuFvbk5C2MTyb0u6USG0zaJYNi2Gcbk1hpL5hxjG5Ujga9Ha6CoLK/OwpLNj+JhCdKtGUm/cv17SxzUSZHRfSRvitW5NuU9fDtS1xkfSPjGkzWskvTxaeGsk/VLSJySdohDS5hY1OIio4zQbFx2n7ZH0JOAQM7s77r8aOBQ4CnghcISkZJX+GWZ2BEFQ3hPXBr0Q2M/Mnm9mfwWsNrNLgRsIYWheaGY7ym67J1Aws9mEEDN/F8s/B3zOzF4CpC2UtwDrzOyFwGzg5wBm9gCwe7JGKcezPhv4H0Konf+JxbMJ7uF/RVinc5iZHQV8GXh3nus6TrvgouN0AvsAf0ztvzpuNwE/A55HECEIQrOJsGDygFh+J3CIpC9IOhbYnuOefyHEeQO4kRDPDUKYmEvi56+n6l8PDEj6GPBXZvZQ6th9pIbxqrAbYSHth8zsyvS1zeweM3sUGAJ+GMtvSbXLcToCFx2nE9hBaYgaAf8RLZQXmtkMM/uKpJcTFj++NFooNwF7RGtjNiFA5hKChVCLx2zEtfNxaiykNrMNhJhovwW+GhfNJuwRn6EWOwkC119Wng4/80Rq/4la7XKcdsNFx2l7omhMk5QIzzrgDI2E799P0rMIkaQfMLM/SXoeIRAmkvYBdjGzbxHC9bw4XqdaOJhKFIA3xc9JuBgkHQjcZ2b/BXwluYckAc8hRK6u+agEh4nnyZOqOZMUf0tyOoUfEiIeX2VmP5Q0E/hp6NN5GHgr8APgHZJuJqRSKMRz9wNWS0pess6Kfy8ghK7ZQRg2y8NS4L8lvZ8w9/JgLH858EFJj8X2JJbOEYS5oVyJ2MzscUknAd+VtJ0QFdtxJg2+ONTpCCS9CHifmb1tgtvxZIKrs0VxONnMTqhS/3PA5WZ2dcaxhmQ6TV3vdDwjqNPm+PCa0xGY2U2EbJ4NWxw6Ro4Afh6tqXcC769R/9YswYlsTy8OHQ+S/oFgweVxknCcCcMtHcdxHKdluKXjOI7jtAwXHcdxHKdluOg4juM4LcNFx3Ecx2kZLjqO4zhOy3DRcRzHcVrG/wdVSvu9d9sGaAAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 2 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "parr = arr.isel(time=-3, yc=slice(1))\n",
    "parr.plot.contourf()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "id": "c982fb53",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div><svg style=\"position: absolute; width: 0; height: 0; overflow: hidden\">\n",
       "<defs>\n",
       "<symbol id=\"icon-database\" viewBox=\"0 0 32 32\">\n",
       "<path d=\"M16 0c-8.837 0-16 2.239-16 5v4c0 2.761 7.163 5 16 5s16-2.239 16-5v-4c0-2.761-7.163-5-16-5z\"></path>\n",
       "<path d=\"M16 17c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n",
       "<path d=\"M16 26c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n",
       "</symbol>\n",
       "<symbol id=\"icon-file-text2\" viewBox=\"0 0 32 32\">\n",
       "<path d=\"M28.681 7.159c-0.694-0.947-1.662-2.053-2.724-3.116s-2.169-2.030-3.116-2.724c-1.612-1.182-2.393-1.319-2.841-1.319h-15.5c-1.378 0-2.5 1.121-2.5 2.5v27c0 1.378 1.122 2.5 2.5 2.5h23c1.378 0 2.5-1.122 2.5-2.5v-19.5c0-0.448-0.137-1.23-1.319-2.841zM24.543 5.457c0.959 0.959 1.712 1.825 2.268 2.543h-4.811v-4.811c0.718 0.556 1.584 1.309 2.543 2.268zM28 29.5c0 0.271-0.229 0.5-0.5 0.5h-23c-0.271 0-0.5-0.229-0.5-0.5v-27c0-0.271 0.229-0.5 0.5-0.5 0 0 15.499-0 15.5 0v7c0 0.552 0.448 1 1 1h7v19.5z\"></path>\n",
       "<path d=\"M23 26h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "<path d=\"M23 22h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "<path d=\"M23 18h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "</symbol>\n",
       "</defs>\n",
       "</svg>\n",
       "<style>/* CSS stylesheet for displaying xarray objects in jupyterlab.\n",
       " *\n",
       " */\n",
       "\n",
       ":root {\n",
       "  --xr-font-color0: var(--jp-content-font-color0, rgba(0, 0, 0, 1));\n",
       "  --xr-font-color2: var(--jp-content-font-color2, rgba(0, 0, 0, 0.54));\n",
       "  --xr-font-color3: var(--jp-content-font-color3, rgba(0, 0, 0, 0.38));\n",
       "  --xr-border-color: var(--jp-border-color2, #e0e0e0);\n",
       "  --xr-disabled-color: var(--jp-layout-color3, #bdbdbd);\n",
       "  --xr-background-color: var(--jp-layout-color0, white);\n",
       "  --xr-background-color-row-even: var(--jp-layout-color1, white);\n",
       "  --xr-background-color-row-odd: var(--jp-layout-color2, #eeeeee);\n",
       "}\n",
       "\n",
       "html[theme=dark],\n",
       "body.vscode-dark {\n",
       "  --xr-font-color0: rgba(255, 255, 255, 1);\n",
       "  --xr-font-color2: rgba(255, 255, 255, 0.54);\n",
       "  --xr-font-color3: rgba(255, 255, 255, 0.38);\n",
       "  --xr-border-color: #1F1F1F;\n",
       "  --xr-disabled-color: #515151;\n",
       "  --xr-background-color: #111111;\n",
       "  --xr-background-color-row-even: #111111;\n",
       "  --xr-background-color-row-odd: #313131;\n",
       "}\n",
       "\n",
       ".xr-wrap {\n",
       "  display: block;\n",
       "  min-width: 300px;\n",
       "  max-width: 700px;\n",
       "}\n",
       "\n",
       ".xr-text-repr-fallback {\n",
       "  /* fallback to plain text repr when CSS is not injected (untrusted notebook) */\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-header {\n",
       "  padding-top: 6px;\n",
       "  padding-bottom: 6px;\n",
       "  margin-bottom: 4px;\n",
       "  border-bottom: solid 1px var(--xr-border-color);\n",
       "}\n",
       "\n",
       ".xr-header > div,\n",
       ".xr-header > ul {\n",
       "  display: inline;\n",
       "  margin-top: 0;\n",
       "  margin-bottom: 0;\n",
       "}\n",
       "\n",
       ".xr-obj-type,\n",
       ".xr-array-name {\n",
       "  margin-left: 2px;\n",
       "  margin-right: 10px;\n",
       "}\n",
       "\n",
       ".xr-obj-type {\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-sections {\n",
       "  padding-left: 0 !important;\n",
       "  display: grid;\n",
       "  grid-template-columns: 150px auto auto 1fr 20px 20px;\n",
       "}\n",
       "\n",
       ".xr-section-item {\n",
       "  display: contents;\n",
       "}\n",
       "\n",
       ".xr-section-item input {\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-section-item input + label {\n",
       "  color: var(--xr-disabled-color);\n",
       "}\n",
       "\n",
       ".xr-section-item input:enabled + label {\n",
       "  cursor: pointer;\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-section-item input:enabled + label:hover {\n",
       "  color: var(--xr-font-color0);\n",
       "}\n",
       "\n",
       ".xr-section-summary {\n",
       "  grid-column: 1;\n",
       "  color: var(--xr-font-color2);\n",
       "  font-weight: 500;\n",
       "}\n",
       "\n",
       ".xr-section-summary > span {\n",
       "  display: inline-block;\n",
       "  padding-left: 0.5em;\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:disabled + label {\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-section-summary-in + label:before {\n",
       "  display: inline-block;\n",
       "  content: '►';\n",
       "  font-size: 11px;\n",
       "  width: 15px;\n",
       "  text-align: center;\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:disabled + label:before {\n",
       "  color: var(--xr-disabled-color);\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:checked + label:before {\n",
       "  content: '▼';\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:checked + label > span {\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-section-summary,\n",
       ".xr-section-inline-details {\n",
       "  padding-top: 4px;\n",
       "  padding-bottom: 4px;\n",
       "}\n",
       "\n",
       ".xr-section-inline-details {\n",
       "  grid-column: 2 / -1;\n",
       "}\n",
       "\n",
       ".xr-section-details {\n",
       "  display: none;\n",
       "  grid-column: 1 / -1;\n",
       "  margin-bottom: 5px;\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:checked ~ .xr-section-details {\n",
       "  display: contents;\n",
       "}\n",
       "\n",
       ".xr-array-wrap {\n",
       "  grid-column: 1 / -1;\n",
       "  display: grid;\n",
       "  grid-template-columns: 20px auto;\n",
       "}\n",
       "\n",
       ".xr-array-wrap > label {\n",
       "  grid-column: 1;\n",
       "  vertical-align: top;\n",
       "}\n",
       "\n",
       ".xr-preview {\n",
       "  color: var(--xr-font-color3);\n",
       "}\n",
       "\n",
       ".xr-array-preview,\n",
       ".xr-array-data {\n",
       "  padding: 0 5px !important;\n",
       "  grid-column: 2;\n",
       "}\n",
       "\n",
       ".xr-array-data,\n",
       ".xr-array-in:checked ~ .xr-array-preview {\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-array-in:checked ~ .xr-array-data,\n",
       ".xr-array-preview {\n",
       "  display: inline-block;\n",
       "}\n",
       "\n",
       ".xr-dim-list {\n",
       "  display: inline-block !important;\n",
       "  list-style: none;\n",
       "  padding: 0 !important;\n",
       "  margin: 0;\n",
       "}\n",
       "\n",
       ".xr-dim-list li {\n",
       "  display: inline-block;\n",
       "  padding: 0;\n",
       "  margin: 0;\n",
       "}\n",
       "\n",
       ".xr-dim-list:before {\n",
       "  content: '(';\n",
       "}\n",
       "\n",
       ".xr-dim-list:after {\n",
       "  content: ')';\n",
       "}\n",
       "\n",
       ".xr-dim-list li:not(:last-child):after {\n",
       "  content: ',';\n",
       "  padding-right: 5px;\n",
       "}\n",
       "\n",
       ".xr-has-index {\n",
       "  font-weight: bold;\n",
       "}\n",
       "\n",
       ".xr-var-list,\n",
       ".xr-var-item {\n",
       "  display: contents;\n",
       "}\n",
       "\n",
       ".xr-var-item > div,\n",
       ".xr-var-item label,\n",
       ".xr-var-item > .xr-var-name span {\n",
       "  background-color: var(--xr-background-color-row-even);\n",
       "  margin-bottom: 0;\n",
       "}\n",
       "\n",
       ".xr-var-item > .xr-var-name:hover span {\n",
       "  padding-right: 5px;\n",
       "}\n",
       "\n",
       ".xr-var-list > li:nth-child(odd) > div,\n",
       ".xr-var-list > li:nth-child(odd) > label,\n",
       ".xr-var-list > li:nth-child(odd) > .xr-var-name span {\n",
       "  background-color: var(--xr-background-color-row-odd);\n",
       "}\n",
       "\n",
       ".xr-var-name {\n",
       "  grid-column: 1;\n",
       "}\n",
       "\n",
       ".xr-var-dims {\n",
       "  grid-column: 2;\n",
       "}\n",
       "\n",
       ".xr-var-dtype {\n",
       "  grid-column: 3;\n",
       "  text-align: right;\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-var-preview {\n",
       "  grid-column: 4;\n",
       "}\n",
       "\n",
       ".xr-var-name,\n",
       ".xr-var-dims,\n",
       ".xr-var-dtype,\n",
       ".xr-preview,\n",
       ".xr-attrs dt {\n",
       "  white-space: nowrap;\n",
       "  overflow: hidden;\n",
       "  text-overflow: ellipsis;\n",
       "  padding-right: 10px;\n",
       "}\n",
       "\n",
       ".xr-var-name:hover,\n",
       ".xr-var-dims:hover,\n",
       ".xr-var-dtype:hover,\n",
       ".xr-attrs dt:hover {\n",
       "  overflow: visible;\n",
       "  width: auto;\n",
       "  z-index: 1;\n",
       "}\n",
       "\n",
       ".xr-var-attrs,\n",
       ".xr-var-data {\n",
       "  display: none;\n",
       "  background-color: var(--xr-background-color) !important;\n",
       "  padding-bottom: 5px !important;\n",
       "}\n",
       "\n",
       ".xr-var-attrs-in:checked ~ .xr-var-attrs,\n",
       ".xr-var-data-in:checked ~ .xr-var-data {\n",
       "  display: block;\n",
       "}\n",
       "\n",
       ".xr-var-data > table {\n",
       "  float: right;\n",
       "}\n",
       "\n",
       ".xr-var-name span,\n",
       ".xr-var-data,\n",
       ".xr-attrs {\n",
       "  padding-left: 25px !important;\n",
       "}\n",
       "\n",
       ".xr-attrs,\n",
       ".xr-var-attrs,\n",
       ".xr-var-data {\n",
       "  grid-column: 1 / -1;\n",
       "}\n",
       "\n",
       "dl.xr-attrs {\n",
       "  padding: 0;\n",
       "  margin: 0;\n",
       "  display: grid;\n",
       "  grid-template-columns: 125px auto;\n",
       "}\n",
       "\n",
       ".xr-attrs dt,\n",
       ".xr-attrs dd {\n",
       "  padding: 0;\n",
       "  margin: 0;\n",
       "  float: left;\n",
       "  padding-right: 10px;\n",
       "  width: auto;\n",
       "}\n",
       "\n",
       ".xr-attrs dt {\n",
       "  font-weight: normal;\n",
       "  grid-column: 1;\n",
       "}\n",
       "\n",
       ".xr-attrs dt:hover span {\n",
       "  display: inline-block;\n",
       "  background: var(--xr-background-color);\n",
       "  padding-right: 10px;\n",
       "}\n",
       "\n",
       ".xr-attrs dd {\n",
       "  grid-column: 2;\n",
       "  white-space: pre-wrap;\n",
       "  word-break: break-all;\n",
       "}\n",
       "\n",
       ".xr-icon-database,\n",
       ".xr-icon-file-text2 {\n",
       "  display: inline-block;\n",
       "  vertical-align: middle;\n",
       "  width: 1em;\n",
       "  height: 1.5em !important;\n",
       "  stroke-width: 0;\n",
       "  stroke: currentColor;\n",
       "  fill: currentColor;\n",
       "}\n",
       "</style><pre class='xr-text-repr-fallback'>&lt;xarray.DataArray &#x27;siconca&#x27; (time: 31, yc: 432, xc: 432)&gt;\n",
       "[5785344 values with dtype=float32]\n",
       "Coordinates:\n",
       "  * time     (time) datetime64[ns] 2010-01-05 2010-01-06 ... 2010-02-04\n",
       "    lat      (yc, xc) float32 ...\n",
       "    lon      (yc, xc) float32 ...\n",
       "  * xc       (xc) float64 -5.388e+03 -5.362e+03 ... 5.362e+03 5.388e+03\n",
       "  * yc       (yc) float64 5.388e+03 5.362e+03 ... -5.362e+03 -5.388e+03</pre><div class='xr-wrap' hidden><div class='xr-header'><div class='xr-obj-type'>xarray.DataArray</div><div class='xr-array-name'>'siconca'</div><ul class='xr-dim-list'><li><span class='xr-has-index'>time</span>: 31</li><li><span class='xr-has-index'>yc</span>: 432</li><li><span class='xr-has-index'>xc</span>: 432</li></ul></div><ul class='xr-sections'><li class='xr-section-item'><div class='xr-array-wrap'><input id='section-c3b42b79-34c7-4e39-9a28-81115f3b2dc5' class='xr-array-in' type='checkbox' checked><label for='section-c3b42b79-34c7-4e39-9a28-81115f3b2dc5' 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>...</span></div><div class='xr-array-data'><pre>[5785344 values with dtype=float32]</pre></div></div></li><li class='xr-section-item'><input id='section-53f0561c-7f23-4332-b7ad-cc09d1602c74' class='xr-section-summary-in' type='checkbox'  checked><label for='section-53f0561c-7f23-4332-b7ad-cc09d1602c74' 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'>2010-01-05 ... 2010-02-04</div><input id='attrs-687a54b2-7c23-4c21-976f-1b6a9fb9f3e2' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-687a54b2-7c23-4c21-976f-1b6a9fb9f3e2' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-6b09ac5e-028f-48a5-be2b-77d4446ca300' class='xr-var-data-in' type='checkbox'><label for='data-6b09ac5e-028f-48a5-be2b-77d4446ca300' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([&#x27;2010-01-05T00:00:00.000000000&#x27;, &#x27;2010-01-06T00:00:00.000000000&#x27;,\n",
       "       &#x27;2010-01-07T00:00:00.000000000&#x27;, &#x27;2010-01-08T00:00:00.000000000&#x27;,\n",
       "       &#x27;2010-01-09T00:00:00.000000000&#x27;, &#x27;2010-01-10T00:00:00.000000000&#x27;,\n",
       "       &#x27;2010-01-11T00:00:00.000000000&#x27;, &#x27;2010-01-12T00:00:00.000000000&#x27;,\n",
       "       &#x27;2010-01-13T00:00:00.000000000&#x27;, &#x27;2010-01-14T00:00:00.000000000&#x27;,\n",
       "       &#x27;2010-01-15T00:00:00.000000000&#x27;, &#x27;2010-01-16T00:00:00.000000000&#x27;,\n",
       "       &#x27;2010-01-17T00:00:00.000000000&#x27;, &#x27;2010-01-18T00:00:00.000000000&#x27;,\n",
       "       &#x27;2010-01-19T00:00:00.000000000&#x27;, &#x27;2010-01-20T00:00:00.000000000&#x27;,\n",
       "       &#x27;2010-01-21T00:00:00.000000000&#x27;, &#x27;2010-01-22T00:00:00.000000000&#x27;,\n",
       "       &#x27;2010-01-23T00:00:00.000000000&#x27;, &#x27;2010-01-24T00:00:00.000000000&#x27;,\n",
       "       &#x27;2010-01-25T00:00:00.000000000&#x27;, &#x27;2010-01-26T00:00:00.000000000&#x27;,\n",
       "       &#x27;2010-01-27T00:00:00.000000000&#x27;, &#x27;2010-01-28T00:00:00.000000000&#x27;,\n",
       "       &#x27;2010-01-29T00:00:00.000000000&#x27;, &#x27;2010-01-30T00:00:00.000000000&#x27;,\n",
       "       &#x27;2010-01-31T00:00:00.000000000&#x27;, &#x27;2010-02-01T00:00:00.000000000&#x27;,\n",
       "       &#x27;2010-02-02T00:00:00.000000000&#x27;, &#x27;2010-02-03T00:00:00.000000000&#x27;,\n",
       "       &#x27;2010-02-04T00:00:00.000000000&#x27;], dtype=&#x27;datetime64[ns]&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>lat</span></div><div class='xr-var-dims'>(yc, xc)</div><div class='xr-var-dtype'>float32</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-e4f913c9-1308-41a1-be7d-efe70c122727' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-e4f913c9-1308-41a1-be7d-efe70c122727' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-aaf6a98e-f224-4d41-9766-f11b9cd69c5f' class='xr-var-data-in' type='checkbox'><label for='data-aaf6a98e-f224-4d41-9766-f11b9cd69c5f' 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>degrees_north</dd><dt><span>long_name :</span></dt><dd>latitude coordinate</dd><dt><span>standard_name :</span></dt><dd>latitude</dd></dl></div><div class='xr-var-data'><pre>[186624 values with dtype=float32]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>lon</span></div><div class='xr-var-dims'>(yc, xc)</div><div class='xr-var-dtype'>float32</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-bcbe8b41-0a43-4bba-ae3b-8a2bed5a5afd' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-bcbe8b41-0a43-4bba-ae3b-8a2bed5a5afd' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-52825347-638a-4272-8177-a36088318aa4' class='xr-var-data-in' type='checkbox'><label for='data-52825347-638a-4272-8177-a36088318aa4' 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>degrees_east</dd><dt><span>long_name :</span></dt><dd>longitude coordinate</dd><dt><span>standard_name :</span></dt><dd>longitude</dd></dl></div><div class='xr-var-data'><pre>[186624 values with dtype=float32]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>xc</span></div><div class='xr-var-dims'>(xc)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>-5.388e+03 -5.362e+03 ... 5.388e+03</div><input id='attrs-c1ce5538-7fb4-41b2-90c8-bc830224cfb1' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-c1ce5538-7fb4-41b2-90c8-bc830224cfb1' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-86e93fae-31c1-4ac8-a808-f2e3e9dcbb28' class='xr-var-data-in' type='checkbox'><label for='data-86e93fae-31c1-4ac8-a808-f2e3e9dcbb28' 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>km</dd><dt><span>long_name :</span></dt><dd>x coordinate of projection (eastings)</dd><dt><span>standard_name :</span></dt><dd>projection_x_coordinate</dd></dl></div><div class='xr-var-data'><pre>array([-5387.5, -5362.5, -5337.5, ...,  5337.5,  5362.5,  5387.5])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>yc</span></div><div class='xr-var-dims'>(yc)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>5.388e+03 5.362e+03 ... -5.388e+03</div><input id='attrs-e20705f7-c8ba-44ae-81d1-fe74581ac1c6' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-e20705f7-c8ba-44ae-81d1-fe74581ac1c6' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-ad0bfd20-a934-4fbd-9762-e2ac1568ee2d' class='xr-var-data-in' type='checkbox'><label for='data-ad0bfd20-a934-4fbd-9762-e2ac1568ee2d' 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>km</dd><dt><span>long_name :</span></dt><dd>y coordinate of projection (northings)</dd><dt><span>standard_name :</span></dt><dd>projection_y_coordinate</dd></dl></div><div class='xr-var-data'><pre>array([ 5387.5,  5362.5,  5337.5, ..., -5337.5, -5362.5, -5387.5])</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-5eb3e06f-462a-4321-b888-6f270ca5a36a' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-5eb3e06f-462a-4321-b888-6f270ca5a36a' class='xr-section-summary'  title='Expand/collapse section'>Attributes: <span>(0)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><dl class='xr-attrs'></dl></div></li></ul></div></div>"
      ],
      "text/plain": [
       "<xarray.DataArray 'siconca' (time: 31, yc: 432, xc: 432)>\n",
       "[5785344 values with dtype=float32]\n",
       "Coordinates:\n",
       "  * time     (time) datetime64[ns] 2010-01-05 2010-01-06 ... 2010-02-04\n",
       "    lat      (yc, xc) float32 ...\n",
       "    lon      (yc, xc) float32 ...\n",
       "  * xc       (xc) float64 -5.388e+03 -5.362e+03 ... 5.362e+03 5.388e+03\n",
       "  * yc       (yc) float64 5.388e+03 5.362e+03 ... -5.362e+03 -5.388e+03"
      ]
     },
     "execution_count": 23,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "arr.sel(yc=)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 29,
   "id": "074fbcc8",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "['../processed/laptop/osisaf/nh/siconca/2010_01_05.npy',\n",
       " '../processed/laptop/osisaf/nh/siconca/2010_01_06.npy',\n",
       " '../processed/laptop/osisaf/nh/siconca/2010_01_07.npy',\n",
       " '../processed/laptop/osisaf/nh/siconca/2010_01_08.npy',\n",
       " '../processed/laptop/osisaf/nh/siconca/2010_01_09.npy',\n",
       " '../processed/laptop/osisaf/nh/siconca/2010_01_10.npy',\n",
       " '../processed/laptop/osisaf/nh/siconca/2010_01_11.npy',\n",
       " '../processed/laptop/osisaf/nh/siconca/2010_01_12.npy',\n",
       " '../processed/laptop/osisaf/nh/siconca/2010_01_13.npy',\n",
       " '../processed/laptop/osisaf/nh/siconca/2010_01_14.npy',\n",
       " '../processed/laptop/osisaf/nh/siconca/2010_01_15.npy',\n",
       " '../processed/laptop/osisaf/nh/siconca/2010_01_16.npy',\n",
       " '../processed/laptop/osisaf/nh/siconca/2010_01_17.npy',\n",
       " '../processed/laptop/osisaf/nh/siconca/2010_01_18.npy',\n",
       " '../processed/laptop/osisaf/nh/siconca/2010_01_19.npy',\n",
       " '../processed/laptop/osisaf/nh/siconca/2010_01_20.npy',\n",
       " '../processed/laptop/osisaf/nh/siconca/2010_01_21.npy',\n",
       " '../processed/laptop/osisaf/nh/siconca/2010_01_22.npy',\n",
       " '../processed/laptop/osisaf/nh/siconca/2010_01_23.npy',\n",
       " '../processed/laptop/osisaf/nh/siconca/2010_01_24.npy',\n",
       " '../processed/laptop/osisaf/nh/siconca/2010_01_25.npy',\n",
       " '../processed/laptop/osisaf/nh/siconca/2010_01_26.npy',\n",
       " '../processed/laptop/osisaf/nh/siconca/2010_01_27.npy',\n",
       " '../processed/laptop/osisaf/nh/siconca/2010_01_28.npy',\n",
       " '../processed/laptop/osisaf/nh/siconca/2010_01_29.npy',\n",
       " '../processed/laptop/osisaf/nh/siconca/2010_01_30.npy',\n",
       " '../processed/laptop/osisaf/nh/siconca/2010_01_31.npy',\n",
       " '../processed/laptop/osisaf/nh/siconca/2010_02_01.npy',\n",
       " '../processed/laptop/osisaf/nh/siconca/2010_02_02.npy',\n",
       " '../processed/laptop/osisaf/nh/siconca/2010_02_03.npy',\n",
       " '../processed/laptop/osisaf/nh/siconca/2010_02_04.npy']"
      ]
     },
     "execution_count": 29,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "dfs = glob.glob(\"../processed/laptop/osisaf/nh/siconca/*.npy\")\n",
    "dfs"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 32,
   "id": "0a82f791",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div><svg style=\"position: absolute; width: 0; height: 0; overflow: hidden\">\n",
       "<defs>\n",
       "<symbol id=\"icon-database\" viewBox=\"0 0 32 32\">\n",
       "<path d=\"M16 0c-8.837 0-16 2.239-16 5v4c0 2.761 7.163 5 16 5s16-2.239 16-5v-4c0-2.761-7.163-5-16-5z\"></path>\n",
       "<path d=\"M16 17c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n",
       "<path d=\"M16 26c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n",
       "</symbol>\n",
       "<symbol id=\"icon-file-text2\" viewBox=\"0 0 32 32\">\n",
       "<path d=\"M28.681 7.159c-0.694-0.947-1.662-2.053-2.724-3.116s-2.169-2.030-3.116-2.724c-1.612-1.182-2.393-1.319-2.841-1.319h-15.5c-1.378 0-2.5 1.121-2.5 2.5v27c0 1.378 1.122 2.5 2.5 2.5h23c1.378 0 2.5-1.122 2.5-2.5v-19.5c0-0.448-0.137-1.23-1.319-2.841zM24.543 5.457c0.959 0.959 1.712 1.825 2.268 2.543h-4.811v-4.811c0.718 0.556 1.584 1.309 2.543 2.268zM28 29.5c0 0.271-0.229 0.5-0.5 0.5h-23c-0.271 0-0.5-0.229-0.5-0.5v-27c0-0.271 0.229-0.5 0.5-0.5 0 0 15.499-0 15.5 0v7c0 0.552 0.448 1 1 1h7v19.5z\"></path>\n",
       "<path d=\"M23 26h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "<path d=\"M23 22h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "<path d=\"M23 18h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "</symbol>\n",
       "</defs>\n",
       "</svg>\n",
       "<style>/* CSS stylesheet for displaying xarray objects in jupyterlab.\n",
       " *\n",
       " */\n",
       "\n",
       ":root {\n",
       "  --xr-font-color0: var(--jp-content-font-color0, rgba(0, 0, 0, 1));\n",
       "  --xr-font-color2: var(--jp-content-font-color2, rgba(0, 0, 0, 0.54));\n",
       "  --xr-font-color3: var(--jp-content-font-color3, rgba(0, 0, 0, 0.38));\n",
       "  --xr-border-color: var(--jp-border-color2, #e0e0e0);\n",
       "  --xr-disabled-color: var(--jp-layout-color3, #bdbdbd);\n",
       "  --xr-background-color: var(--jp-layout-color0, white);\n",
       "  --xr-background-color-row-even: var(--jp-layout-color1, white);\n",
       "  --xr-background-color-row-odd: var(--jp-layout-color2, #eeeeee);\n",
       "}\n",
       "\n",
       "html[theme=dark],\n",
       "body.vscode-dark {\n",
       "  --xr-font-color0: rgba(255, 255, 255, 1);\n",
       "  --xr-font-color2: rgba(255, 255, 255, 0.54);\n",
       "  --xr-font-color3: rgba(255, 255, 255, 0.38);\n",
       "  --xr-border-color: #1F1F1F;\n",
       "  --xr-disabled-color: #515151;\n",
       "  --xr-background-color: #111111;\n",
       "  --xr-background-color-row-even: #111111;\n",
       "  --xr-background-color-row-odd: #313131;\n",
       "}\n",
       "\n",
       ".xr-wrap {\n",
       "  display: block;\n",
       "  min-width: 300px;\n",
       "  max-width: 700px;\n",
       "}\n",
       "\n",
       ".xr-text-repr-fallback {\n",
       "  /* fallback to plain text repr when CSS is not injected (untrusted notebook) */\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-header {\n",
       "  padding-top: 6px;\n",
       "  padding-bottom: 6px;\n",
       "  margin-bottom: 4px;\n",
       "  border-bottom: solid 1px var(--xr-border-color);\n",
       "}\n",
       "\n",
       ".xr-header > div,\n",
       ".xr-header > ul {\n",
       "  display: inline;\n",
       "  margin-top: 0;\n",
       "  margin-bottom: 0;\n",
       "}\n",
       "\n",
       ".xr-obj-type,\n",
       ".xr-array-name {\n",
       "  margin-left: 2px;\n",
       "  margin-right: 10px;\n",
       "}\n",
       "\n",
       ".xr-obj-type {\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-sections {\n",
       "  padding-left: 0 !important;\n",
       "  display: grid;\n",
       "  grid-template-columns: 150px auto auto 1fr 20px 20px;\n",
       "}\n",
       "\n",
       ".xr-section-item {\n",
       "  display: contents;\n",
       "}\n",
       "\n",
       ".xr-section-item input {\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-section-item input + label {\n",
       "  color: var(--xr-disabled-color);\n",
       "}\n",
       "\n",
       ".xr-section-item input:enabled + label {\n",
       "  cursor: pointer;\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-section-item input:enabled + label:hover {\n",
       "  color: var(--xr-font-color0);\n",
       "}\n",
       "\n",
       ".xr-section-summary {\n",
       "  grid-column: 1;\n",
       "  color: var(--xr-font-color2);\n",
       "  font-weight: 500;\n",
       "}\n",
       "\n",
       ".xr-section-summary > span {\n",
       "  display: inline-block;\n",
       "  padding-left: 0.5em;\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:disabled + label {\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-section-summary-in + label:before {\n",
       "  display: inline-block;\n",
       "  content: '►';\n",
       "  font-size: 11px;\n",
       "  width: 15px;\n",
       "  text-align: center;\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:disabled + label:before {\n",
       "  color: var(--xr-disabled-color);\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:checked + label:before {\n",
       "  content: '▼';\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:checked + label > span {\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-section-summary,\n",
       ".xr-section-inline-details {\n",
       "  padding-top: 4px;\n",
       "  padding-bottom: 4px;\n",
       "}\n",
       "\n",
       ".xr-section-inline-details {\n",
       "  grid-column: 2 / -1;\n",
       "}\n",
       "\n",
       ".xr-section-details {\n",
       "  display: none;\n",
       "  grid-column: 1 / -1;\n",
       "  margin-bottom: 5px;\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:checked ~ .xr-section-details {\n",
       "  display: contents;\n",
       "}\n",
       "\n",
       ".xr-array-wrap {\n",
       "  grid-column: 1 / -1;\n",
       "  display: grid;\n",
       "  grid-template-columns: 20px auto;\n",
       "}\n",
       "\n",
       ".xr-array-wrap > label {\n",
       "  grid-column: 1;\n",
       "  vertical-align: top;\n",
       "}\n",
       "\n",
       ".xr-preview {\n",
       "  color: var(--xr-font-color3);\n",
       "}\n",
       "\n",
       ".xr-array-preview,\n",
       ".xr-array-data {\n",
       "  padding: 0 5px !important;\n",
       "  grid-column: 2;\n",
       "}\n",
       "\n",
       ".xr-array-data,\n",
       ".xr-array-in:checked ~ .xr-array-preview {\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-array-in:checked ~ .xr-array-data,\n",
       ".xr-array-preview {\n",
       "  display: inline-block;\n",
       "}\n",
       "\n",
       ".xr-dim-list {\n",
       "  display: inline-block !important;\n",
       "  list-style: none;\n",
       "  padding: 0 !important;\n",
       "  margin: 0;\n",
       "}\n",
       "\n",
       ".xr-dim-list li {\n",
       "  display: inline-block;\n",
       "  padding: 0;\n",
       "  margin: 0;\n",
       "}\n",
       "\n",
       ".xr-dim-list:before {\n",
       "  content: '(';\n",
       "}\n",
       "\n",
       ".xr-dim-list:after {\n",
       "  content: ')';\n",
       "}\n",
       "\n",
       ".xr-dim-list li:not(:last-child):after {\n",
       "  content: ',';\n",
       "  padding-right: 5px;\n",
       "}\n",
       "\n",
       ".xr-has-index {\n",
       "  font-weight: bold;\n",
       "}\n",
       "\n",
       ".xr-var-list,\n",
       ".xr-var-item {\n",
       "  display: contents;\n",
       "}\n",
       "\n",
       ".xr-var-item > div,\n",
       ".xr-var-item label,\n",
       ".xr-var-item > .xr-var-name span {\n",
       "  background-color: var(--xr-background-color-row-even);\n",
       "  margin-bottom: 0;\n",
       "}\n",
       "\n",
       ".xr-var-item > .xr-var-name:hover span {\n",
       "  padding-right: 5px;\n",
       "}\n",
       "\n",
       ".xr-var-list > li:nth-child(odd) > div,\n",
       ".xr-var-list > li:nth-child(odd) > label,\n",
       ".xr-var-list > li:nth-child(odd) > .xr-var-name span {\n",
       "  background-color: var(--xr-background-color-row-odd);\n",
       "}\n",
       "\n",
       ".xr-var-name {\n",
       "  grid-column: 1;\n",
       "}\n",
       "\n",
       ".xr-var-dims {\n",
       "  grid-column: 2;\n",
       "}\n",
       "\n",
       ".xr-var-dtype {\n",
       "  grid-column: 3;\n",
       "  text-align: right;\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-var-preview {\n",
       "  grid-column: 4;\n",
       "}\n",
       "\n",
       ".xr-var-name,\n",
       ".xr-var-dims,\n",
       ".xr-var-dtype,\n",
       ".xr-preview,\n",
       ".xr-attrs dt {\n",
       "  white-space: nowrap;\n",
       "  overflow: hidden;\n",
       "  text-overflow: ellipsis;\n",
       "  padding-right: 10px;\n",
       "}\n",
       "\n",
       ".xr-var-name:hover,\n",
       ".xr-var-dims:hover,\n",
       ".xr-var-dtype:hover,\n",
       ".xr-attrs dt:hover {\n",
       "  overflow: visible;\n",
       "  width: auto;\n",
       "  z-index: 1;\n",
       "}\n",
       "\n",
       ".xr-var-attrs,\n",
       ".xr-var-data {\n",
       "  display: none;\n",
       "  background-color: var(--xr-background-color) !important;\n",
       "  padding-bottom: 5px !important;\n",
       "}\n",
       "\n",
       ".xr-var-attrs-in:checked ~ .xr-var-attrs,\n",
       ".xr-var-data-in:checked ~ .xr-var-data {\n",
       "  display: block;\n",
       "}\n",
       "\n",
       ".xr-var-data > table {\n",
       "  float: right;\n",
       "}\n",
       "\n",
       ".xr-var-name span,\n",
       ".xr-var-data,\n",
       ".xr-attrs {\n",
       "  padding-left: 25px !important;\n",
       "}\n",
       "\n",
       ".xr-attrs,\n",
       ".xr-var-attrs,\n",
       ".xr-var-data {\n",
       "  grid-column: 1 / -1;\n",
       "}\n",
       "\n",
       "dl.xr-attrs {\n",
       "  padding: 0;\n",
       "  margin: 0;\n",
       "  display: grid;\n",
       "  grid-template-columns: 125px auto;\n",
       "}\n",
       "\n",
       ".xr-attrs dt,\n",
       ".xr-attrs dd {\n",
       "  padding: 0;\n",
       "  margin: 0;\n",
       "  float: left;\n",
       "  padding-right: 10px;\n",
       "  width: auto;\n",
       "}\n",
       "\n",
       ".xr-attrs dt {\n",
       "  font-weight: normal;\n",
       "  grid-column: 1;\n",
       "}\n",
       "\n",
       ".xr-attrs dt:hover span {\n",
       "  display: inline-block;\n",
       "  background: var(--xr-background-color);\n",
       "  padding-right: 10px;\n",
       "}\n",
       "\n",
       ".xr-attrs dd {\n",
       "  grid-column: 2;\n",
       "  white-space: pre-wrap;\n",
       "  word-break: break-all;\n",
       "}\n",
       "\n",
       ".xr-icon-database,\n",
       ".xr-icon-file-text2 {\n",
       "  display: inline-block;\n",
       "  vertical-align: middle;\n",
       "  width: 1em;\n",
       "  height: 1.5em !important;\n",
       "  stroke-width: 0;\n",
       "  stroke: currentColor;\n",
       "  fill: currentColor;\n",
       "}\n",
       "</style><pre class='xr-text-repr-fallback'>&lt;xarray.DataArray (time: 31, yc: 432, xc: 432)&gt;\n",
       "array([[[0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        ...,\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.]],\n",
       "\n",
       "       [[0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        ...,\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.]],\n",
       "\n",
       "       [[0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        ...,\n",
       "...\n",
       "        ...,\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.]],\n",
       "\n",
       "       [[0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        ...,\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.]],\n",
       "\n",
       "       [[0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        ...,\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.]]], dtype=float32)\n",
       "Dimensions without coordinates: time, yc, xc</pre><div class='xr-wrap' hidden><div class='xr-header'><div class='xr-obj-type'>xarray.DataArray</div><div class='xr-array-name'></div><ul class='xr-dim-list'><li><span>time</span>: 31</li><li><span>yc</span>: 432</li><li><span>xc</span>: 432</li></ul></div><ul class='xr-sections'><li class='xr-section-item'><div class='xr-array-wrap'><input id='section-286100d0-45bd-431b-85c6-54fefb078813' class='xr-array-in' type='checkbox' checked><label for='section-286100d0-45bd-431b-85c6-54fefb078813' 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>0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0</span></div><div class='xr-array-data'><pre>array([[[0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        ...,\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.]],\n",
       "\n",
       "       [[0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        ...,\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.]],\n",
       "\n",
       "       [[0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        ...,\n",
       "...\n",
       "        ...,\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.]],\n",
       "\n",
       "       [[0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        ...,\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.]],\n",
       "\n",
       "       [[0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        ...,\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.]]], dtype=float32)</pre></div></div></li><li class='xr-section-item'><input id='section-90b576bb-2310-4593-ae01-853fbd867d02' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-90b576bb-2310-4593-ae01-853fbd867d02' class='xr-section-summary'  title='Expand/collapse section'>Coordinates: <span>(0)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><ul class='xr-var-list'></ul></div></li><li class='xr-section-item'><input id='section-2e8eef31-d8dc-47b5-9c05-568152771e2e' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-2e8eef31-d8dc-47b5-9c05-568152771e2e' class='xr-section-summary'  title='Expand/collapse section'>Attributes: <span>(0)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><dl class='xr-attrs'></dl></div></li></ul></div></div>"
      ],
      "text/plain": [
       "<xarray.DataArray (time: 31, yc: 432, xc: 432)>\n",
       "array([[[0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        ...,\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.]],\n",
       "\n",
       "       [[0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        ...,\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.]],\n",
       "\n",
       "       [[0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        ...,\n",
       "...\n",
       "        ...,\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.]],\n",
       "\n",
       "       [[0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        ...,\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.]],\n",
       "\n",
       "       [[0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        ...,\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.],\n",
       "        [0., 0., 0., ..., 0., 0., 0.]]], dtype=float32)\n",
       "Dimensions without coordinates: time, yc, xc"
      ]
     },
     "execution_count": 32,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "narr = xr.DataArray([np.load(file) for file in dfs], dims=[\"time\", \"yc\", \"xc\"])\n",
    "narr"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 42,
   "id": "19374d50",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.collections.QuadMesh at 0x7f74fcb205b0>"
      ]
     },
     "execution_count": 42,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYMAAAEGCAYAAACHGfl5AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAABi1UlEQVR4nO29eZwkV3Xn+z03InKpvbp639RSSy20IoQszGaQMbuNbGbskZnB+A1+jOfB2Mx4g/Hz2GNG7/G84vfGy8iAjTFClm3AMotBgDG7hBBCaEFSq6Vu9d5dXV1rLhFxz/vjRmRmVWV1Zalr7/v9fPJTmZERkTeqKu+Je5bfEVXF4/F4PBc2ZqUH4PF4PJ6VxxsDj8fj8Xhj4PF4PB5vDDwej8eDNwYej8fjAcKVHsD5sHHjRr1o9+6VHobH41kD3P+d75xW1U3nc45dUtYqtqN9T1P/rKq+5nw+bzlZ08bgot27+drXvrbSw/B4PGuAclfXwfM9Rw3LT8m2jvb9Ez248Xw/bzlZ08bA4/F4lhMBApHOdl5jJVzeGHg8Hs8CCDq0BWsNbww8Ho+nQxa0MlhjeGPg8Xg8HSICBeONgcfj8VzQuJXBSo9iafDGwOPxeDpGvJvI4/F4LnSE9Vup642Bx+PxLAC/MvB4PJ4LHBEfM/B4PJ4LHsFnE3k8Hs8Fz3quM1jyWIiIBCLyHRH5ZPZ6g4jcLSJPZD8HW/Z9t4jsF5HHROTVSz02j8fjWSiBdPZYayxHYPwXgUdbXr8L+IKqXgZ8IXuNiFwJ3AJcBbwG+BMRCZZhfB6Px9MRLmYgHT3WGktqDERkJ/B64P0tm28GPpQ9/xDw4y3b71DVmqo+BewHblzK8Xk8Hs9CWa8rg6WOGbwP+FWgt2XbFlU9BqCqx0Rkc7Z9B/DNlv0OZ9umISJvA94GsGvXriUYssfj8bTHIOs2gLxkKwMR+VHgpKp+u9ND2mybJQKrqrep6g2qesOmjWtKLtzj8awD/Mpg4bwYeIOIvA4oAX0i8tfACRHZlq0KtgEns/0PA623+juBo0s4Po/H41kQecxgPbJkKwNVfbeq7lTVPbjA8BdV9d8BdwFvyXZ7C/AP2fO7gFtEpCgiFwOXAfcu1fg8Ho9noeRCdX5lsDi8F7hTRN4KHAJ+EkBVHxaRO4FHgAR4u6qmKzA+j8fjmZP1ujJYFmOgql8CvpQ9HwZeMcd+twK3LseYPB6PZ6F4CWuPx+PxIAKRWZ+6pd4YeDweT8cIsk6XBuvTxHk8Hs9SIGAC6egx76lEXpNJ7+wXkXe1ef9XROSB7PGQiKQisiF772kR+V723n2LcWl+ZeDxeDwdIoAE538PnUnt/DHwSlxa/bdE5C5VfSTfR1V/F/jdbP8fA/6zqp5pOc1Nqnr6vAeT4Y2Bx+PxdIqwWG6iG4H9qnoAQETuwEnyPDLH/j8NfHQxPnguvJvI4/F4OkU6cxFlbqKNInJfy+NtLWfaATzT8rqt/I77SOnCiXf+fctmBT4nIt+ecd5njV8ZeDweT4eIQBB1LKZ8WlVvmOtUbbbNkt/J+DHgazNcRC9W1aOZttvdIvJ9Vf1ypwNrh18ZeDwezwKQQDp6zMNC5HduYYaLSFWPZj9PAh9nERSevTHweDyeThFBAtPRYx6+BVwmIheLSAE34d81++OkH3gZTdkeRKRbRHrz58CrgIfO99K8m8jj8Xg6RKCjtNH5UNVERN4BfBYIgA9mkjw/n73/Z9muPwF8TlUnWw7fAnxcnCxGCNyuqv90vmPyxsDj8Xg6RUAWqZ+Bqn4a+PSMbX824/VfAn85Y9sB4LmLMogWvDHweDyeThEhKKzPbrzeGHg854mKIDpXIohnPSGLV2ew6vDGwOM5T2YaAjM5jO0eWqHReJYaswgVyKsRbww8nkXGG4J1jKxfoTpvDDwej6dDBDCLFEBebSzZekdESiJyr4h8V0QeFpH/nm3/LRE50qLG97qWY96dKfg9JiKvXqqxeTwez7NCWKw6g1XHUq4MasAPq+qEiETAV0XkM9l7f6iqv9e6s4hciSu8uArYDnxeRPb51pcej2fVIEJQWHsTfScs2VWpYyJ7GWWPc6Vc3Azcoao1VX0K2M8ilFh7PEtF4ej3iJ66Z6WH4VlGZB2vDJZ0xCISiMgDwEngblXNvznvEJEHReSDIjKYbetIxU9E3parAJ46vWhS3h7Pgok3XUp88QtmbTeVkRUYjWe5WKzmNquNJTUGqpqq6nU4EaYbReRq4E+BvcB1wDHg97PdO1LxU9XbVPUGVb1h08aNSzJuj2c+arffikbltu/Z8uC018FDdy/HkDzLQVaB3MljrbEsaxlVPQt8CXiNqp7IjIQF/pymK2ghKn4ez4pSfNOvz/leeOZp7Cf/Px57yxv59itfxcR9X1/GkXmWEkEwgenosdZYymyiTSIykD0vAz8CfF9EtrXs9hM01fbuAm4RkaKIXAxcBty7VOPzeJaMkwcxP/qf6BrqZseLL6H8s7+50iPyLBayaBLWq46lzCbaBnwo6/VpgDtV9ZMi8mERuQ7nAnoa+A8AmWLfnbi2bwnwdp9J5FkrFI49jJqQeMvlJM95GQC7/uDDKzwqz6IjgonWZ3nWkl2Vqj4IPK/N9jef45hbgVuXakwez2ISHX+UtH8r4dmjPPHe/wsbJ+z4oetIq3WCUoHCLe9e6SF6FhmR9StHsT6vyuNZQqITjxGePUw6uBNbHuS+/+NX2XjtXnb9yA18/X/8IxIYirsuhn/+y5UeqmfRWbTmNquO9bne8XiWiGD8BPGWyxuvC8ce5rm//Gb0xbcA8PI3/grgfKDJx36X5MPvoeua52Ove12703nWGlmdwXrEGwOPZwGkvVumva5vuwq2XdV23/CNv0IIJJ/+Yzj6xwT9Qw2j4VmrCGLWpzFYn1fl8SwD0cH7MJURzOTwOfczr3s7ptwNJiB60qeZrmVEBFOIOnqsNfzKwONZAIVn7qd+4GGS4eMcf/BJ4rHJzrKGbvpZFIiXfISeJUXALNLKQEReA/wRrgfy+1X1vTPefznwD8BT2aaPqepvd3Lss8EbA49nAdR3Xc9Xf+KdlAdLVEaqvOS+L6/0kDzLzGLEDLKU+z8GXokruP2WiNylqo/M2PUrqvqjz/LYBeHdRMvEcNWu9BA8i8RL7vsyW67bxrX/24sIRn2R/AWFLFo20Y3AflU9oKp14A6cWGcnnM+xc+KNwRIylSiTiRIrbH7c69OsJ3b+7ocYes0bMNXxlR6KZxkRQIzp6AFszEU1s8fbWk7VkTAn8MKsJ8xnRCTPVOj02AXh3URLSFfYLElPr301pj6JLXSv4Ig8i8nRO26nONBD39vP213rWStkK4MOOa2qN8x1pjbbZgpz3g9clPWEeR3wCZxMT0eingvFrwyWEW8I1hfV4VFGnzxC4ch3V3oonuVCICiEHT3mYV5hTlUdy3vCqOqngUhENnZy7LPBrwyWGYtgzt+Ie1YBF/27n8KOnyXZsHulh+JZJkQWrc7gW8BlmSjnEVyXxzfN+KytwAlVVRG5EXfzPgycne/YZ4M3BsuMNwTw9SMTvGhHT+P1eKz0RmtP5dFc/FzEhKjxX6MLicXIJlLVRETeAXwWlx76wUys8+ez9/8M+NfAfxSRBKgAt6iqAm2PPd8x+f9iz7LTagiANWkIAJLB3ZjKCBJXodi70sPxLAcLixmck8z18+kZ2/6s5fn/BP5np8eeL94YLCIHx2Mu6m1feTiZKL21MxAWsC0TR3R6PxoUSAa9q2EtMrOrmWf9s17lKLwxWCBPj8Xs6YtcJequ66e9N5chAOgOBRsOTdtWv+P/Ji0VSKaqlPZeQXjx1cSb9y3JuHOC8ROkvVsoHHuY+rarODyRsLNncf8NxuqWkVrKV54e4eholYuHutkzUKa/FLJvsLCon+XxLCciggmClR7GkuCNwQLZ0+cm/JmGYD6CiVPw9ANoEjcULHO9+6d+7l+zM7WUXvDGxR1sG3KhtXomrrazJyRWODIeUwyEbd0L/5f41rEprCpPn63QVwwJjPDUyBTHz1YJjHCmUufQyBTlQsC+G7Yv6vWsZVQEUR9DWlMImPkzhdYkS3ZVIlICvgwUs8/5O1X9TRHZAPwNsAfX6eynVHUkO+bdwFuBFPgFVf3sUo1vuUl7NmGqU2g6u3nbpe//O7fPJ36faPc+7PU/tqxjiwT+9OsH+Z3+7zB6/7fpu/oqCpdeS33Hc9vuP5m4CSy16lYB1Zg4tdx08SAnJmJGawn7hrrZ0lME4IU7+wkEJmLLM+MJ9VRJVS/IVUI4cghb6MZ2D1GJLYlV+grr0+2wPlm/qqVLaeJqwA9nBRMR8FUR+QzwRuALqvpeEXkX8C7g10TkSlyK1FXAduDzIrJvNbW+jJ78OvHeFz27g79yO/al587+Ova1Byl9/yCbKpOY3kGCTTsbd/BLzf/zmktRLqXvhT9J9NQ91B//DlF1EsICUihT33YV0cnHqW2+nB6pI9VxCAt0dfdST7swAluKSk9UpJYWqCaWwAipVbYWUxBDGARM1C3j9YQjY1WsKlt7IjZOHr5gYibJA18kuupF2O6hrChxbQbPL1TE9zNYOFkK1ET2MsoeitPQeHm2/UPAl4Bfy7bfoao14CkR2Y/T4PjGUo3xXJyYStnSNd03aCfPQ3pgHkPwletfwsbLh4gnq5hv3kv/3h2Utl787D/vPIgvfgFc/ALs+AnkmYdJjh8iMiHJxospnn6CtHcLac9GTH2KUGCwFFBPldFEiK3SHQmDxRDUEhOQimA0ZdDEaBRR6C2wuTuiYIRSICQDOzGTw9iuQZD180WTtI4GBYIHP8vE/d9k8vgwQanAhlI3bLqUZyYsw5WYrii4IFdJa5JFzCZabSyp8ytT1/s2cCnwx6p6j4hsUdVjAKp6TEQ2Z7vvAL7Zcvii6G08W3JDEIweJe13fu702lcv2ee99P6v8uBPvJ70xCRR9zBdmwcJljiYPB9p7xYm920mEOiqn0WDAvHGSwEQVTQqAzBYdF+OYOw4k91bqCRKMTSYuEbRJiAGW+hGwyKhKsVAMAKlQLAIoGAMk6lgdW3WHLRDAzfBp9e+mvK1r6YMFI5+D1vspi4BO3oDNneHFNfn3LJuWa9uoiW9KlVNVfU6XLn0jSJy9Tl270hvQ0Telgs/nTp9epFGOje5IQCofvg98+4fjJ941p917cc/RWWkigSGqLfrWZ9nMekO3d17JymUad9WSoGwsX6qYSxssXeaDEdvJHRl5wSIJk5iauNosZfuQKmn2ohJrEfq268hGboEgxJWznhDsNYQQcJCR4+1xrL8K6rqWZw76DXACRHZBpD9PJnt1pHehqrepqo3qOoNmzZuXJLxho98kX85NM7Xj0xM2147Oz6nDo2pOReSRqWOu1kFE6emvb735a8gLIdEXWXUWqKTjz+L0benmjr11OVgZmvIuZC0jha6qBV6qUtIOHKIraNP0Fc5Of/B6wBb6s9WRp61g4AxnT3WGEs2YhHZJCID2fMy8CPA94G7gLdku70F18mHbPstIlLMNDcuA+5dqvHNRazw2PaX8rLdvbMqZfv/0+80MmxiBTM5jMQVotP7Cc46u3XW9HYcZE57Nk17PXW6QlJJOPP4UU4/+CTxo/cswhU5SoGwmrwvktQAJ94XictoSvu3k/ZsROqThI98cYVHuPQEk8MESXWlh+FZCAISBB091hpLGTPYBnwoixsY4E5V/aSIfAO4U0TeChwCfhIg0+W4E3gESIC3r0QmUSSwd+Dc/Uujk4/D5n3YbldEFj96L4VLruJkJSU9j1TB0mCJ6tkq8WRMZaTKpusuY+39S3WGhkUkrU/fZkK3xE5iavsfhKcfpXjFDS6gvYp5tuKDKgbMev0Lr1dk3f7NljKb6EHgeW22DwOvmOOYW4Fbl2pMi0W86dLpG176JurA5rZ7d84PfuWLfONFL6c25ibJ4QefZMfz7iO+aC5J9LVNHmAFt9KKxKBhkWTDbooveC3pkcdJh48T9e5vBK5XIwsxBMHoUWzPRkxlFGyKjJ+EIGKyvLERR/GsYgRvDDwtLGH6Y1gKSeuW0ZOT7P/Uw5SG/o6BFx5b9kK0ZUUtkRhQi9QnkaSORkXMJc/DlvupFftcfvc6qNb9nYdT/s21hotsiqT1Rn1FaYXH5ekMWcdFZ+vzqpaRmUHg8+XG9/wc8WSdQIRDR8Z54hP3c+Tj/0jh2Hkr1K5eMuNq4grB+ClMZRTbPYTtGkSjMmF9gujMwRUe5Plzz9FJdg6USVXBBA1D8MTZOk+ejTk2mazwCD3zIgJhobPHGsMbg2eJ3PMxgofuJu0emn/nBaAvvoWbHv0mrzrwLbb0Fvj+I6d46nOPMfb5jxGdeGzRPufgeLxo51osbKGbePM+4s37kPoU4OIIpHVXozB8gOipxQuqLzf3Hj7LB7/0JM+M1qYlD1zWH7K3PyAw3k20FlhAD+Q1xdob8QoTPXUPevefo9VJ0qtfuaQuo5se/SYTifLIwVEe/quvMP65O+ErtxOePXze5z6XwuqKoLbxVJIaqEJcxdQn0WIvac8mtNgLXQNrNtNoqLvA5bsGuOfQCIcnEk5XUmKF8OxhzNQIm8vr0xe9rpAsgNzJY43hYwYLQS3xxS9AljG7JRAYjS1fuv8448c+w4bLBrn2v0C07/nEWy5ftnEsOS1GVcMiad9Wl22k6n7vEmLCArZrABMERCcfX3K578UmMoYvfPFJejeUeeFFg7xoexfB+HFMdRxs0shO86xm1m82kV8ZzENeTAasiG7Ozx1/gECgai2PnZpi/7ePc+RTd1O75zPrO44AoIqGTvk0wrpq5q5BbNRF2j2EqYys8AAXxvXb+th28QYmzlZ5ZrSKqY1jy/3YUi+YkJOVVaPJ6JkLX2dw4WJXQTvDd55+kDt2Xs9EYqmkzp0yfvAoowf+gvJQP4X+XorXvBDt20QydMm0Y4erlqHSwo2YxJWG9tBKoWERUxlBC90gxqWfBgW01ItUx0n7t3P/iSlGqwk3XdS3omPthL0DETu2dCMGpuKUA/V+9vRFJIO7CcZP0B35e7PVjyxadbGIvAb4I1wf4/er6ntnvP9vcSKe4EQ//6Oqfjd772lgHCf3n6jqeeef+/++NcIth+/nedds5rk/uINCbzdpPWHsqWPs/8Q9jB44jBRKaFhC0joSV0AtRycThisJ3zo2xV9/7yS3fOS7DYmNYHSW0kdjmyQ1NCoTnXyc8MzTy3mZs7DlQZePn9Ybwg1qQsQmBKNH2dTlsjb+5dB5KMouI5dt6aU6FXPPgTM8PjzFgdGYsbrFdg/hywzWAIukTZQV4/4x8FrgSuCnMxn/Vp4CXqaq1wLvAW6b8f5NqnrdYhgC8CuDNcXz7/4cx37jbfRdvI3ijt1sevVW0lNHkEKJZMNutNjLaGLYMHkcgA19u5ioW0qhoZpart7Rj1VXMStdg4TDB5A0afjec1E+DYs8eTamp2cvIudfTPdsySuUc1dRiCVRQ1DoQgrdSBpzcf0QF20epF5aG72ID5yaYGq0xjfveYb+ckR86RCvuHgANWGj1mCkZhtKsJ5VyOKsDG4E9qvqAQARuQMn4/9IvoOqtoqcfROn17ZkeGOwSokO3odERRCDTo0y8sVPc+Srj2DrKVt+4EqiK17A8MClFK90/5iJm+XpioT64EVYVWyq7OiNSK1SDDdQ3CtsLIeMVFMCU2RzdZz6gYeRpx5CZ7TcbJXkqKZKKRCCseOkfVuX7XfQWqHsNlhCAVSxXYOYuAImQOoVwmfZmP7geMxIJaErCtjdF9F76F7shl3T1GoXi+GqpRAG9G3s4olvfIt7Njk118uGuukvBmzpCjhbt9i1X1u3fhFBFieAvAN4puX1YeBcmSlvBT7T8lqBz4mIAv9LVWeuGhaMNwarkMIz92Mnx5DBLagJ0VqVqLvM5uddQhCFhDv3AtBXMEwmSpzNHgUjhAKxKnGqFAL3GoSeQpFqYqlbpZZaQhXq268hLPcz9bmP0r35HugaaJuhVNY6SpFq71ZWIiHVpZpaJ4mNkKoSSIAt9BAEEZLUnrU+UCkwDJRCTkzUKYbCfdFVXBKWWIpOzUMlw3cePM6ZIyfZctlVbNrczcHhST792EleuHuQnkIXfYWASmLnP5lnhVhQNtFGEbmv5fVtLZN2R5L9ACJyE84YvKRl84tV9WjWD+ZuEfm+qn6504G1wxuDVUh91/WuV64YzNQIFEsUd+ymvM+1g9DKJOnhxzC9m+guD9L6f2UmhykUuojCImdjRRUKgSAoRqAcGgrG9RSQuEKyYQ/dN72R5MCDYI4QnDqE2XpxQwsoOngfdnAnWigTFXtZqTCThkVUBNSJCVogVcUEBTQoEI0fR4s9qJh5A9+tQfWdJ+4j3X4le/pcosDx8TrffGaUx05O8KbnbccghAa2dZ/fV2W4avnL+49wz2/e1Nj2mSfPcmS8ypEzFQ6NVrlioxt3MRBXd7GOur6tG4SFuIlOn8Of35Fkv4hcC7wfeG2m6waAqh7Nfp4UkY/j3E7eGKxH4g0XUUss3SaE7iHXMzR2VcP1k4cxSQzlQU5WUvqLQaNJiu0ewtQniTFMxQnb20xipnIGG7gWkxYh3nI5kVr0zDHqTz9K/TvfpDI8Sv8lO4n7h+DQ45y59z5OfucAV37gw8ueYSVp3d2NWcAEWHWrACMgNkFNSLVnCwVNkLjC8TYtS1tpza7KFVHzhjo/sL2H8XpKuRDwl986zKa+IpcOdVMb6mJP37NbF339yATvvv0BTh8d45de1Pz+v3bvAACf2j9CLbFUEqWnAEYEUx3Hlvo5OB6vvgLBCxgRQaJFkZr4FnBZJtd/BNf/fVpvXBHZDXwMeLOqPt6yvRswqjqePX8V8NvnOyBvDFYpos5Pn/ZtxUwOQ/9mtNCDRkXCi6+fNiEnVim2SBnYQjdRWmdrd5F2K09Jao2MoVSVaqp0b74MNu0lGD9L7cmDPH33Qxx/5EsMbOvh6p95EWGpwHP++hOshAMjv1ZJ64hal02U1rFhEVGbGYsCdQmRQi+FeuejHI+ViXpKXzGgO1BiVV52UR/XbunmqwfPEgaGvYNdbCgFPDOePKtVwkfuO0ytGvPwH76u7fuvv9TFO0xlBBmvo1HJifbZxBuCVcfiFJ2paiIi7wA+i0st/WAm4//z2ft/Bvw3YAj4ExGBZgrpFuDj2bYQuF1V/+l8xyS6hpUgn3/99fq1r31tpYexbASjR0n7tnJwIiVJm0Hek5V0tpRB5mYYj2f3FA4mTpH2bCJRZyoEiDQBtZiJ05gzz1A/8DCP/9UnOf7ACZJKwo4bt9G9uRebWi55688saT/ouZC44i4tLDZWNbmaaaLujjq2iqoyVrdM1C2JdbGTdnf1pjKC1CuITagM7CYQIaxPYAvdJAqjNctEnNITBQyWAkZrKZ978gw3XTy4IOmIK37xUzz6R6+ftu0Tj5/hmbMVhroLBCJcMtjFZUMlBp5lLwzP/JS7ur59vmmYN1z9HL33Yx/oaN/g8pec9+ctJ35lsEYwk8PIiScJkyql8kWcrjUVLttOTNlk2e7GMu3ZhCQ1wiBq7Iu1YFOnFFrqJdi6j72jwwzufYozjx/l7MExjtx7jEtfeykHP3w7u358BH3xLUt0te1pxALU4rw6zpoFRhqRjKKB8MxBeusV4q1X8M8Hx9jYVYCW0LekdTQouBqGLAup+8h3STZdSi3qJsKttgID9VQ5MFGlKwqIAmGwvPA79VZD8PDpGo+dnmD/6Ukq9ZTRqZgrt/YSW8tjpyvZWOdvsORZIWRxVgarEW8M1ghBZQTp3UBS6qeaWq7aWJz3GKMppj5JcPborCwhSetuUjTuX0DDIuSpnGkdbEjh9T/P1lfX2WoM9t5PUjtyiMKmzWhlEk1iwpFDDRnmZUUMJguaBzZGqhV3HdmX1Jb7CWxCePYwL9+zC7EpwfgJJKm58YqBf/kwvOzNjVOmIycJojKyaR/V1K2WDdBfDIiMcGKijhFhR2+p41VBnAW7W3ns9ARTcUo9sfR3RfSXI8LAMFFPMSLsGTAUQ+F0JWWjF65bnazTwP6SGQMR2QX8FbAVF/q7TVX/SER+C/jfgbwRwH9V1U9nx7wbl0KVAr+gqp9dqvGtNeKNl7qK4LHj7NnSz4l2rqGZiHF6Pm3SRc3EKeTscae7Xu4jGdzluoxhwBQwgXOZYIyTg/ihf0vXxCnMmWdIjh8iHT5GfO+nMT0Dy75CANddzNU/hGBCNCqjIq4BTgFsveKMRiaFHUycJn76EXTiswRbdqO9A1Ruv5WRx59h+xteT7BlNyQ1SmcOUNlwCZOxRRUGSwFbCimXJu7f1RZ7oeaM6OOT7vd/2UD7gOJkbLnt3sP86kuafQuKoZtIrtzay9beIr2FkFJkqCdKbC3FUIiMsEEnkbFJECHt3bLUv05Px4g3Bs+CBPglVb1fRHqBb4vI3dl7f6iqv9e6c1aKfQtwFbAd+LyI7FuJPsirlXjzPgpHv0d4+gBbu4c4NNHPzp5n9ydMhi6Bh76GhBFarxLc8FpseQACt+IQcYFoSSOwKcHUCNgUyn1Euy8n6B0gOXHIHfvgZ1ckhlAKxN19R2VMbby5OlB1qqdZ5zTEkAztwZw6TP3wAaRQYuLxx4gnXTN6rVWpb78GxCBxhWJtlEKhm7E0YDK2EBUonzxIfPBRNE0pXfNC1Fo2bbuep8/WeOBkhes2N9NZzeQwtnuIgYJpGIK8cc1gOaKWWAbLEVu6IwqBm1hS1ax/dkBYn8BMnMJMnUWLPc6tZ5PZRXgZq0FH6oJBaKym1xtL2QP5GHAsez4uIo/iqu7m4mbgDlWtAU+JyH5c7uw3lmqMy00wevS8K1vr268hPPM0ktTY2BO0dUV0zMvejOJqCZg6i5qQsCtyVc8tu4lN0EIXElddM58ghKE9hNsudV3Jxs4QHVz+Xs1Jy7W3ZldJUgMx1DFExV6XPaUWe9VNFC9/EWZqhMGde0mOH+LARz7OyS/+C1uv/iGqPVuY1CJBUMJY6I4MI9WUeqqUt1xMMDXG+He+RXz2MxS2bKPvoh9gsBwSSPMPEH7/Xzjzhc8w8vgzXPRHH2mO1boajx19RQxCT8FQDIQokEaVSIglOvH9rB3mTrTQTdq7hclE6Q4LDSMzE43Kc77nWWzE3SmtQ5bFxInIHuB5wD3Ai4F3iMjPAPfhVg8jOEPxzZbDDtPGeIjI24C3AezatWvm28vHAoqCzOSwc7eU+zGVEcIzh6gfeBhe+qb5D26h8Mz9pCOnoNSFGdgCPXB6KiEycl7+5fiiGygce5hw7DhxqRcNCogqYhMn/VDozp5bpDYJWW6/c0N1IxuKaKHrWX9+p+SB3/x3H7Z+J1v/HiKITYjEuHEHEZrGSOrqNJKBncQDuynsfh57i2UmH3oAnnqA0kXXUAoibNcgEzYgVaWnYEiscrp7J9G1u+jfcw0cP4AZ2MTBSSdlcdmGZgdj0zsAgI2bAf4nztYpBYbuyFAIhEIgBDK7YjocOUx67EnsddNTULuzC51rsg/GjrtUW28Mloc12MWsE5bcGIhID/D3wDtVdUxE/hSnwKfZz98H/j0dlmdn5dy3gUstXapxz4eVAKtKoTZKcOz72PGz7V0lap0GvwQYTRExJBsvwQzsJHjqnkbRUydo4iYzrVddv2AjpOqSlM9FJ4He+rarADexaKGMLfZm6qCpm1DBrRgKZbAWLXRhjQHpdXflYjD1SWyhu+Pr6RRTGQETNu/+xZAo041Ba3OcTNVU0hhJqi6mEBadvVCLpDFBWHJB8+e+it5LnouMnoDqOMmG3Uh9kt40cddtU/f7KPWTIEz0X0Rpwx4A+hIlMkJohOjEY24FFZXp2bWV0w8eaIznbCVhd3+RIKsAN5oiaeLGrNatZJI6nHhqliHo6PdTHUPXYM/dtYgC6mMGC0dEIpwh+IiqfgxAVU+0vP/nwCezlx2VZ68WxuopkRGiqIz2bSIolJj88HuIJyuUNw8SvvFXGnnxmACCgKoVClEX0cRJNAihdwi+cvu0FUJ45mmSbLKZiZR6MCZAezaQFroRYFNXSIlkulsnrrhJ0FowBts1SDVVAiPzupRMdQyrFvKJN8vQUTFoIWxU/FoJCPLJLJ80l+hLYmeI0OU1BaCNz8YE7X25JgSbIIkzEnnBWphlUklax/ZsIohrzoAkVWypH1FLTUKK8WRjZRQaQ4LLEgJ3x94dCIXjD1N75B7OfPtBNr/u9ZjeAfoubgr6RdlqoBiaxmpA0hjSOpImmInTiFrqz3nZgn83wdhx0t5Ns35HniVCfAB5wYgrj/sA8Kiq/kHL9m1ZPAHgJ4CHsud3AbeLyB/gAsiXAfcu1fjOl4GCYaxuEa1RH3JSz11v3EkwOQyVMdLqKBoUCM8cchPWxksoRWVAscUeMAGpCQn6h5BHvojuvgYALXQTjB7F1CZI+7cRjJ9o6ARpVIQgwvZuwRa6Sa3Lq9eWP2N45ml3V1/oQqMyNQtGoCROfTSaRzQ/3rzPpWFmAUs1IRJXsJFzA0kYIDbFaIqNSm7loBaxtWWTqagkSnckWIQgjRGbYMPss3PjlBVTqgmRVqMFkNTd3Xih7IyFWuLNl7l+yyZETUBsDXGqRNlKR02IRSiGQi2x1FOFyCAYIhOy/yOf4qkvHeLlu7dQuvRKBp9zUWO8A6WQcigYcYYkkAAJIoKpEUxltLEqOxfhyCFnuKIuCFxQH7VOnsMbgmVE3P/MOmQpr+rFwJuB74nIA9m2/4pr4nAdbsX1NPAfALJS7Dtxet4J8PbVnknUVzBYejEosQXt3kjaszGfhzCaNiYZyO7YbeJcFkHBySlc+2pXRIVxmjTZ3a7t2Qg2Je3d4ozDmWeg3Oe+/IVuUEtRYJZw3OlnCPo3o4UyyeBuWmXxSx12T5H6JGJCyOIHGpWb/m11rhYASZyrQ8W4ZvXLQDVVIgNW1Y0o6moZWxY3EOOcjirNzdC4E89XEq2oCbDFXiSuYJIaYVgkEEHiqYYERpjWIalT7NlIYLQROI43X8bWGy+nMlLl5Le/z47BQZKpasN9Z7PCOMEZskCUrrCALXaTDOycHX9Si6lPIllabNq7ZZqbLzx72B3nWRG8m2iBqOpXaR8H+PQ5jrkVuHWpxrSUWAXBqYQ6VU3BipsOJJv8AVRtY7u1znUjMyIfqXHVp8akSK2OFruxG3YhSRVJE+fXj1wKaH5XaGrjrkq5WEYLZecyaRfknifwHY4cQo89idm2FxsWsurkFA2iuY9bpi9HrBBbKLUEC6wqiDg/fObCUjHOYLVKrahFknrjORI07/DEIDZ1P5MqWujGJDVnrOtT7rzlATQoYOKqK2LL/oYGBTH0XnU1W0+e5eR3D9H/1CE0teRlgXnAGMCIE9izCCYLhksaI/EUtjzYcHvZYm/TVTcDDXx18orijYHnXDTuwNViapNUwx4C4yaBNMvdt6qAaUxSRqQRoC2kceNLLmERVVyQ+MQTLmWz1O+yknD+ZlObxEZlgolTSG0cSWIX6BzYjpbcJGLFZcQIUJgaRtJ43uY0yeBuyHryQlMPCLUNTaDc1567VCbqdpr+UTB+4vwLpTKjlaiThkjV3Vnki5vEKjaT584n5DxwnAdmEXHb1Lrx5+9nqwKJq9hSn8s+ytJR837LqMXEVbTcj82MhqmNU+3ehFUoijqBPFXnLnvpm9ix7/l03fVh4skqtu5cUvefmOKSwVJj9VIO2wW66y6+A5ipEdKeTef81fgitBVEfGqpp1Oyqt8S7k7WqmJEGm6N1CoiLqAIOHmdoNAwGPk2EUiLPejOa50vO4gaLqjYKkXjJqe00I3JUygBVJHqOBoWCfUsQVQmDbLUT5tQTZWRako5MucURUt7txA+8kVM/xDxluc070a1qQiqJmAqdv7z4VQZKhkKR75LfcdzZ53PVEY68m1LXGm0ucx+FUDTCASZOmu+osoDyfnvPk89zb+uotYt6zO3F3EFqU267UHUWE00P1DBxs4dF5VcfKBuKUeuZmGsaumOhMaiN/t9HByPSYqXsOktv83QI/+E1qpMJEpXFDTcc5bm/wGADYqEWEQrjXTFtGcTNQu+6+Uqxq8MPAslEjJBteYkFhpB1c05BkVNQGoz18E0t4ZLndSseQuASWPUhM3eBYVurARIPsnFVXc3HJWbvnObEmjVyVUjlLKpqBJbIiONHPZ2JFf+MMHEKTc553GCNAa1pFEXT43Uia2lmli6ooDhCuzbfk2jTWaORUhKgwzP6DMQKxQ0oS4hVkFVkaBMMXOj5B3c8lVBa/BbgMh00N0sjZu+SsUZhSCCpOaMQFxFolLWPMf9YiWpOdeYwljdklplMnYGoZ5a+oqZnpMJqRIxMhHz2OkpjoxV+c6hs/yXl72STV2uvWhPwRkDFcFapZI4VdX8fyIy7veqxd5GumxRLGBcVXWLi9GzOvAxA8+zIl8V0LiLba4WnCuihgkirIlmVVXkLhLIpBckJGmdaMVkBiWEoOB84mrRLPOHXOsfSMMSE/WUU1MpxVDoK5g5DYGpjKBhyRmQ3GWRfQE0iBC1nKkkdEVCqgGRMZyeqnNsosanH5uiUk959w9dxGTi5KNtJrVQDIXh7M66FLg0V5XQ6YlKltZa6s8yoJyBzH3tgTR/fyLOmKZ5vABtG4R1F+PkKZwrJslqDgqN/H4Ni06YL4gaMQQNi0hSI0hr9EQR43VLbJVKbNnYFWbjFpCAapxyeKzOSCXGiFCOAj7xyAl2D5S5bKibjV0uG0twoX4TdRFbbWQkjWlAPxCMHEK6Bt3qKbuO5W4i5OkAEVeBvw5Zn1e1imi927ctsQJwPn0pdCHZXWhuOFrdStXUBaVLgSswE5GGqmbD/ZCdJ+8R3Hq3rFEZsQkGpa9g6CsYphKl6xwrAi32UtGAUpv3gslhpD7Jpv7tmMkzYALGejYy1BWwo6/I3sEujAgnKynl0FBP1RkdI3S1zSdo+V0Ve1ERUmtJcO0fDUpRbKOIzLl8AtKg2LaKF8h8/llqaWa8lELz08U4gyAGiasuNdemiCSN2EKevlmODFEgxGnTOFmAzBhFRtjcE7G5px9VuHJTD7G1lELDlp6InshgI6dHBEDqfnZFhjRbDZFYlxgwfsoFkb3W0CrG1xl4zoOZK4JWA2EQ0iyo3LpfTnfk8t2nEqWWKvVUZ7V0NChTsZvg202OM4uxurWKMvdkoyacZggkqTWyidxKYZPLeS/3o1GZbmAygZ094TThvMlE55n+p2PFuczKkWmsiNwdfZrFAyLUhC4Nd+bBM76gjbhDvkLIgslAFkR22Va55pLm+6o6oxNEEJhm5lDgXHwmrjZ+F4GtE0ZFBm1AOXKfv6M3cimp2YWn2Sog1Wa8yIirGo+MYDPXVbx5X3PsWWLAfIFkzwrhjYHnfJlpCGD6amHme3loNAqESmwZPEdU8Vx3+jNZ6F2nqY7NmpjqpsBIPWVzFldu53I6VzyiHYGNCQC1hiDLhKpKwRURZ5lDiXX5/Wn+u5prdZCPPavx0KDg7vbBuYvyjKMs8C5qUTWNLKk8K0kzV5wRGkZCTeC0m+qTBNVxurqGCPOGoGIQ62I7ibrmOI3rM+JWO63CdkLbySVvQNQaTPesDtZrzGB9XtUqwrRMXPnKoHUyyJlpCGYSGuHEVMrRyWTOfWrWtcCccyyVkcbz5Nwf1yB68uuuJ+8MTkwm09Q6O6GRpjrX+7Vxt484/Z5I3N1zhKWkdUxSo5DWMJoSyLmlNZwkh+vzrGGxKVvRusoICpjapAsYZxXKrauovF4hz9TKU2mnYovNK1GNaf49MwOTmoiabRryunUxk0Dctvz/IYinZmVPteINwSokl6Po5LHGWHsjXmMYlNC03Mm2oDN+zkVqXfzAqhJbdSmrNGMHY1kD+KKZ3gLT1Mabz7OgsKmOzvmZk4nrHVxNtXHueO+L2mr+9BQMQ6W5/30m21ibc61InJ+8yymmxhVMXEHiCkFaa8hLOKVUp0hqsuL0fFVgWx1SajNBurzITLPgcdDMEc+NQqs7Ka27L3FrdXJWo9AYp6oTm8OpudpSf/O97POsuoyhQFydScEIQeYSSm3TEBBEzdWKZ+2Q1xrM95j3NPIaEXlMRPaLyLvavC8i8v9m7z8oItd3euyzYV5jICL/l4gMtLweFJH/sRgffqEgqkQt2TCQFZzN2E/bPBrnAHoKARtKoQtm0swq6pujXsAWe4kO3kd0en8jS0UL3VikMSm10p1lGZUCmZYa2u4OdYMdn7WtlbjN+dshNpnW6D4N3SpEs4k6r8GwUal5t9WQ17ZznNS4ibrY29Twyc+ZidWR1WZoyypBw+YKKN8n/9lKwy2V3QFOe41z/biAv7ogcyZUFwVO4TQQIY26GmNZr81S1icy7e92rsc5zyISAH8MvBa4EifTc+WM3V6L02i7DCfb/6cLOHbBdLIyeK2qns1fZL0HFq6z6yGIp9r6twMRlLzCVhqP1Cpx6iptXSaLyz4Zrc12BZ2tt58Y44tuaArdhUVI48ZqJaeatp+4g7HjmPokT4wmfOvY1LT3pDbJ02NOUhu1qLhc+mqqJMo5C9qmYVM36UdlYoV6al3/gaiMmsClmaIu40okm7SLbSfptuRunrzat8WIuMwk06xSbqTPFtBcrdUEDXkQcO61eMYjX5XoDPdAJNMfDS0jq/PXR3hWL4vjJroR2K+qB1S1DtyBa/DVys3AX6njm8CAiGzr8NgF08ktSSAixawDGSJSBrwz81lgqmNoUkW7Nky7M7fazLrJ3UlpJrcQBc4oVBLnIlKFbd2z/2zzTb6mOtrQ59H6FFGxp1FMVgwDYuvGEGaFcvVUOckQewoRfUk6q89v2r+dPdnzGEOkOqvYrB2SOsVQW+yd1p3LIsRZ2mWFgEBdxk1J6xBnwnJimqmzc92tt3xOTuMuTdVJXYfZ5E/QOCqXhdCg0OilrCLM9O7VslTZRJ3onBEoBGA6MExWXTbYQoL9ntVFftPTIRtF5L6W17dl/VjANe56puW9w8DM5ibt9tnR4bELphNj8GHgCyLyFzjPxb8HPnS+H3whkvZuwVRGENVGY5b8rjI3AkJTd0dxncwmYsvZSsIPbHv23cRafdumMoKpnMWW+pxCarG3EYytWRivpWwsB+zpc3fExVBm9fnNSdRVM6eBWxWIiKugnWNy1KBAMH7SSUJnhiBRUJxoX16dGxqn2wRNN5UkNbeKyipyZ9ZUzP4wbfQ5yOWpsSBim1W9WZBYoaFh1Ho4TDc0xbzvsjBNj6mT9qNGBPUrgrWNMusG4RycVtW5esG2+2+Zeea59unk2AXTiTGo4vxTz88G8R5V/ez5fvCFSDBxymW51CsNwTjJg5U6/S+swOlKwvHxOrFVXrKzZ9HG0aoRlP8Hmdo4ttjL4bGYM5WYjeWm4RkoGIbKTVdJa9FbKFCODIEIpcAFt2varANonSBNZQSCAmnvZmJTwDSuOCsQU2cQCraOVFxWUWvQOQ2d4Fui+TfC3Zm3o7U9JtDy0xkTizSC0GpCl8ZrUxegzlYF1iozb+KN5MdNN3ad9KEO0hrdviPZGkfnzfzrkE6aec21T6GDYxdMJ8ZgC/CLwP3AB4HPne+HXqjkufqFI991RVvFXggLWavIZj+DqVQ5U005PFZbNCOQuz4afYRnkEsf7B2I2DswO8NloNTMsCkFTkrDZhNiRJrNzoZcw02zmIdFCGtjSHUcUUvaXXIpnTNubsIshVRs3EgvzZvr5Ksnq0rqTt+iN9TM759FI+BsG5lIrnZAGuN1KaeJSy+FhutMyFNGp3/xDUowdtztW+qdVzIiH7vR1GsMrQNcavCinOpbwGUicjFwBLgFmNkU/S5cv/g7cG6gUVU9JiKnOjh2wcxrDFT1/xSR3wBeBfxvwP/MmtB8QFWfPN8BXIi0qnqayWHEhERRmZFaymgt5amRCl1RsGiGwN0Fq5vAwyIyx52NRaa5q4wIgXXCdOWw6Uu3CGZmFk2bzJ7I1jGVUczUCLbYQ9qzqaF+2mz/WCcQg6pz5ZDUISo1eymrMzlqApKsdef0u/UOgsh5umjmBhKbIriKZzNT9tq6PtXK3HEALXYjNTfWoH4CLfZAWseWB6cpjlZTzWpLoGTwMhPrBF2ElYGqJiLyDuCzuPDYB7MGXz+fvf9nuN4vrwP2A1O4+XfOY893TB3ltKmqishx4DiuC9kg8Hcicreq/mq7Y0RkF/BXwFbA4oInfyQiG4C/AfbgOp39VJahhIi8G3grkAK/cEG4o4zLPDg8kXBkrIYR4aaL+hbl1HlAt3XSFnVulsDINJ947n+f7nbRRvC1UBvPMnzC2X76NoZAVJGkhtTGSXs3ubTNlv6xeWUw4O7M1TpN/yBsGoIZ5ysZ7Sh/exYt/QyaJ5zeVwAzPQW12RNhtkGwpX4o9btVlk0JJk41OpEVjXMHopZqcRORcemkDTVZz5rG1fss0rlUP82MZl+ZEcifK/D2To89X+Y1BiLyC8BbgNPA+4FfUdVYRAzwBNDWGOCMxi+p6v0i0gt8W0TuBn4W+IKqvjcrlngX8GtZnuwtwFW4HsifF5F9q7315fmS++939jBN12cxyN05Mye0MPeSiJtkUXvObJhw+AC2ZxOk8WwffAu5lIN7HmOLzo0SjB13TWOyfXKZh7ziV7LJV8Nw9oSZZxBp2tQMmqlO2kEmT6shsBI07tjzc0zrHjbHNebZRbkx1KCAUEfD4rR+DRqEaKGbgcDXdK5H1msKQCezz0bgjap6sHWjqloR+dG5Dsqa3h/Lno+LyKO4lKibgZdnu30I+BLwa9n2O7IU1qdEZD8un/YbC7kgzwyyCS1fDeTB0fxn6z4z5Rhyt4Yt9bsm8iZxFbNqp/m/8wYyrUFbDSIKxx4m7R7Cdm+Yft6WCVxskvVsmLutZhC7God2d9aN5jVzMaMZT2LVXTtOC4nUNbJpBJVbvuoNf38mkZ11qJ5+ehNiS32YyigSVpxBLHT7+MB6RRdvZbDa6CRm8N/O8d6jnXyIiOwBngfcA2zJDAVZMGRzttsO4Jsth+U5tTPP9TZcNR67du2a+banDc7H33ydt9qktXArm1RzldBc+lpq4xA2J3lJas6N01qkld3ZWwQkcJ4cCbDFbpc5lAeapUUfaVqBztyywJLWG+mh7Ziz0rM1sJw9z5VQQ6FRfJanrbbKWeRKGkZm1y/M/AxJYwgibM/GRtW0zX4Hc8VmPGubxYgZrEaWfB0rIj3A3wPvVNWxc+3aZtus37qq3qaqN6jqDZs2blysYa57Wl1CksYNbR6gWaULjWBrntGDWpcxkyZO6llM8z2abp/GZ4jLwxZVkg17kKRGmOnyWARb6icpb3B3zvmEfI676Eant4UIf2XaRDlWAqwEBJL57/PYQB48pikU2GoU8m2NfWb+i+aSGS0SF7k8heSSGZ51RZ5N1MljrbGkoigiEuEMwUdU9WPZ5hMisi1bFWwDTmbbO8m79TwLGnnxuQtnptbQDJcNAMEMSxwW0DzrpkXIbbq7qJmJBK5ITJIqNVNqSGqYTJFUkipYixaXJqja6jrKA+SpKgHSyBZqve4Imr0SWvoRNLSNshXIQmQkvObQ+mS9uomWbGUgIgJ8AHhUVf+g5a27cAFpsp//0LL9FhEpZvmzlwH3LtX4Ligy98+5xN1aey3POjyXng4Kc2fzZPpErZjKWSeMlxkCaNHxKfZiuwaXZsJsow8jQkOPKVa3WgCmiYpJUiMUGm06G5fWGnxGZq8Q5mA8nn/WSJQ5/yae1YfrX64dPdYaS3nr8mLgzcD3ROSBbNt/Bd4L3CkibwUOAT8JkOXY3gk8gstEevt6zyRaTjQoTG/0Mg8NA2ACF6TNVxNZJtDsA0xjZZAbhWrPlka7SMmK0IJAXH7/Atw+88pOdED+3bSZtPRMpdE0KELgVlFhizpqMy7iVGaztstYFYK0NqdgntiEPmNRbf9+npmUqhIan3W0llivpnvJjIGqfpX2cQCAV8xxzK3ArUs1pgsd25LeWTeFc0ooNDJ38kBx7mI61518HmTO/OWxbQrBGYH5upK1HXPLv1Br6qqVYEHncsqnCcaEQDOjauZntAZ/E6VRFQ5Ng5Kf71zBazUhktQwcQUzOexWIIUyktSdRlVSQ4MC52he51mlrMGb/o7wTs0LEDVhR1o6wLRsnDzg2lojMC0I3ZrmmcaUIpe7v5DMmnarAKuKGOfrz8c0nyGY6coxaDNbCqZnUrXsY2kaiRCL2Gxiz5YEmgmVNc7eqqQ6091jAifL0bul8Ro7SjB+AokrJBv2dPIr8awiXNHZ+rQG3hhcAOR6RA2J5qyHb6dM85kHUSYG2iJBka0IbMt5c7eSmyTbn3daHn9GY1LNK5VRDLa5Nu/AvZTo9LTQ/HWuP4QJGhIUea8Cm032IjiZ7+ya1YSN6xQFySb/XJIjl8uYSSM+Y538tiQ1iJMsEG+aBsKz5liLmUKd4I3BBUAeI9Cg0NThWUiqZgvS0nshF32TLG20VXnVZg3tLcxS/mwln4SBhoroXFpHnQabQwFJa6RBkSCtEUFjUm40srEp2BgyY5C7g1oNQWs3tFYaxqthJKYXvrX2TBYgb9kpSR2JJ+cOwnvWBOt0YeCNwQVHGiNB1D4IPI+0w0zF00RdOprQnCBjzXoR4O6k2yl/5hic+0c1NwrNOEBjgj3HeM4ZWLYpYTrhWlqa0GVC5ZeZnVeS2jRdprzgruESMuHs31M7I9VujJkYnrsYQVKnv6RRcZqEuGdtoSh2nQpSeGNwoZGpd7aVsj6HOBtAbArT/mHcHb82g7lqicQ0Kn3zu+28onfmCqGxmhAwCLaDlcS0Y1t/6vS2MWGhqymWZxPSoDjtjl5sgoZFxCYE0FI0ps2K6PkM5gzDMG2FkNdi2BRTm8j0irqcWynrHeFZg6hfGXjWCXnsAJopoK3xAzUB9VTbZrnMNUmLANYFkFMJSKylEJjG5D7r2GxCNeqqoBVpGASXtknTfXOOVQW0xB1EGoG9xLo6gcAIBEXXzlPTWRXBuaFouI6yi2m2yWyV3Mh+V7Q3BG6fbN9WNxN1t2+u2+RZ86zXojNvDC5QGhkyM0itUgicfEQ960nspKPP4T7KpC3SbIWQF3eZbAJ1PRFoNpbPCr4aPzM3kWnJvc/7KcxXY9AaJM4JMvkLwRmGoFWYKUt7Pae4Hcy63lzKI98+Z2prm5WDRiXncgqLflWwxnFyFOvTGnhjcAHSruhsMptNA3HZEoHQnETnCHhO6/urliDLrw/FJf+IOOMS4O6a202erQ11Zrp68m2NAPM8vtrWdyPjVhmpZndyLeO0UdlpFdkYNUFj3d9YIeQTfqshmpE5lccaZo2rTSCZzBAA3hCsA9apLfDGwOPoDqUR/M01+891Vx4rBCJu0jdBdowjb3kZkTWBl3P/m+UhZoWme6eFhrFgehFYK6HQ6I0MEGdr+chIVmMQNFJL88KyvDDM1QO0KLBmtJ3kM6PQGqswQvtYQpbB5Fk/+DoDzwVBlOktBJPDaLEHiUqN2XWsbikEQilwRiO1zdk/IBd2y4TfMsmHRJVQzv3FMbi+xsB0PSCaK4XW7YFhRje2Jm5FIg3lVKvaCGwb20wDDdJa8yBp0wBontf5Z4gqcZ4x1VK5PDN9t13nNs8aRSFdp6EfXwzvmYXtHnJun+oYktaJFboiQ2CEqUSJ21TdNLR/soIskayjmkiz6nfm5yANzZ9W8sl/5qfkK4ZzicXl8tGtktQxWXFZ3gfZpjRKiVWzdNIZ3/As46o16GsRqtZpEuXxlqKZfqeoeSaSbdZMzHX9nrVHfpPSyeN8EJENInK3iDyR/ZyVjywiu0Tkn0XkURF5WER+seW93xKRIyLyQPZ43Xyf6Y2Bpy2S1p1/O5vUaqlST5XEukd+sx6aPAOo+c9vVZuVyOpE4GJtTuKtk3l+mOAm+1kuopZVQf688YVr1RRCGp8BzhC1ZhhVE0saFJuFd1l6aWOintEDofF7yPz+scJEPSXV2TGX1jFLLnmdx1laYg3eKKwHlFQ7e5wn78K1B74M+EL2eiZ5a+ErgB8E3p61D875Q1W9LnvM2y/Zu4k8bbGlfpehEzrBujOVlLPVhK4ooLdo6DbSmJwLmmS6QS6w2ipXnU/GgUgzlbRFJE5Fpi0BUquI5P74NuPKpCUCaa4SnGRF/rnTD8ozmgiaiqVAVoXsur01s4ucKWnEDRqpogFYpacwPYNIkhoaRA2ZCkndeRqS2Hnaqk0agWlTHXVtRD1rElXaroyXgJtp3x64ZSxzthZ+5Nl8oDcGnjlprQ0YKgdY3OqgEEgj2ygQISZsTLr55G+kmSXUmLhzv3qrW0YMIc3+AhhnEFq/b3mqqDuvO19+boHZ2UZqnWvI2kaaa8MY0Vpg1tKox6YgQUOEr6FVDZg0JgwiN+6sNmFak6CsWG9mumprT2hJau55mmAmh50rzrPmWMYA8lztgdsyo7VwzjtE5GeA+3AriJFzncMbA09HdIXCrt4Cp6YSDG6CzjOPoH1aaG4EWhHbMiHPCMzmd/dixDWszyZ9waW6tmaIBi3GJ1UXx+jKrZcYAgSMG6fY1N21BwV3p571ZGgor+Lu7KcVhal11aaZlLdkvQ0kqSKxm9i1e8hFtMG1EgU0iBpjyOUtMEEziDyjyZxn7bEAF9BGEbmv5fVtqnpb/kJEPg9sbXPcry9kPHO0Fv5T4D24r817gN8H/v25zuONgadjDMqxiTqD5ZByZEitKxTLVwWxdZO5u4vXxt14LtNgJSDGNJvL5HfhuHaYjdVBXl3cUnyWf76bsJspnBY3DmV63YPJdC5Sq81/8rwnQ/5ypvDdrEIzC3EFUYstdLuGP3lHuCwu0GgVmqUwSRo7g6DWvVaLDZamtadn+XErg453P62qN8x5LtUfmes9EZmrPfDM/dq1FkZVT7Ts8+fAJ+cb7FK2vfygiJwUkYdats0Z4RaRd4vIfhF5TERevVTj8pwfQ+WIJHV+05mZFaqaGQjn6skNAWnspCpaFU8zWpvJN33/zv2UGxqDEtg466lcm9a+M/9eBtmJE3WPmQ1r0qCIlcBtb6M7lLt4xCYuu8gmkNaR+lRDT0iSqns+OUxw9mijbqEhDR4Wnd5RGmPiCmITn1a63sjSqjt5nCdztQducI7WwmQGJOcngIeYh6XMJvpL4DVtts+KcGcR8FuAq7Jj/kREfLXOKmRnX8SmroCTUwn1dLpryCqN4K+IZNW9Fo3KxBI2MoZgdvOZVhpGIXs0JuiWzB/XdS1oxBds9iXNi9fy84Qyu3K5kf5qk0Zwd1omUUNsLgRjkPqUMwrZexoW0WJ3s+fBjOwiDYuoCb0hWIconaWVLkJc4b3AK0XkCeCV2WtEZLuI5JlBeWvhH25zg/07IvI9EXkQuAn4z/N94FK2vfxyFtTohJuBO1S1BjwlIvuBG4FvLNX4PM8OUaWvYDg+mXBqKqG/GFAMmhN7nkmUrw40yJ3kzVRTq0pIM4B8TtQ2XDKtk2uioOp0kGyau5UgIHcV6eziMaYHuBsqrdNiBe64XJpCwxIalprXn7m1bKkXgsitemyl2SY0P03ogwPrEaVZ3b6kn6M6TJv2wKp6FHhd9nzO1sKq+uaFfuZK1Bm8Q0QezNxIeSHFDuCZln0OZ9tmISJvE5H7ROS+U6dPL/VYPe1QSzk0jNdSRmspxyYSjkzE1FIXyI1TbcQT8pxrI9LITgppTtKS1M5ZRJb3XXZ+eRdXSLJ7/XwlEAXiOptJS11Cm5qB3AUFWSOdTGBv+k7NBamKmV17kGcaZQZDo3Lb+gTPOmX53ETLznIbgz8F9gLX4fJjfz/b3m4maPvbVNXbVPUGVb1h08aNSzJIzzyIwWYT62Tdct/RUR4+OcHpSoLSXB3UW+r206xILLAu66YxWYsQxFME8dS0WIALDrtCsgRDoi5AnLZkLRUCQ5QFrPNHYCQr8Ipdhs+MCmLIitEkcI1v8rqA3CjkekKZUWhtsCNqXbeytD6tQY13B104LFcF8kqwrNlE54hwHwZ2tey6Ezi6jEPzLJCLeiPIBDiv2uhcImN1N/GmCsVAGnIUjU5muB7KJo3bt7BMYzddm6Dhhw+z4jZlerGPEbCS1x1AXSFqvbVp06AmF9/L5bIbu2bZTtJwGZnm8XkjnCwdFRNCCKYy4juWXaCs1x7Iy7oyOEeE+y7gFhEpisjFwGXAvcs5Ns/501cwdIdCdyhESQVTGcHEVYymhNiGpHWuFipxBchy+U0wrZAr708gLRN37gpqBIjFGYvASKO7WkOpuk3LymlSFa39lmmRmxbTKDbLu6GlxR7isOzGZgwaFrwhuEDxK4NngYh8FFdOvVFEDgO/CbxcRK7D/U6fBv4DgKo+LCJ34sqoE+Dtqvm31bMW0ajcDKo2GsdrI26QT/rVVKnWLRASGSGwUAiKRGktk64wiAkIs3RTq9oIEktaRwkblc2QfVmzxjOtchiKTiuCUxNgklqjEC2n0Q85SxtFDKJKlGSGKyzNH/T2rFtU2ws1rgeWMpvop9ts/sA59r8VuHWpxuNZQfK777iZdZOou0MvBcL3h6v0FkJ6CobQCPU0pRw5N1GUyVcgLm7QuPMXQcOia3GZ2sYStzX45AwCjWBzFLQcmwnoAQ3padSdpxG7SOtIEk6LZfjmNJ61eNffCb4C2bNstKZftuoejVRiRiox/aWQ/mLEQCkgSF2BWipQCNwsHolOUzttVbsoR2aaRAVZsVquXBoaaQrSSdCQFGhIaog4wbm8FwFpllZaRL2GhCdD8W0vPZ4l46aL+gC4/8QUqSqlUBoVxZKlpwZZExoDDUmLvIlOEEgjHjEVW1KF/qIzIK5wLc8uqkPoKpFzl1Me3G64kFSReBLEkJb62vaJ9lzAKNg1mDbaCd4YeFYNV24sU0+bkhUi0hCom0lgY9dXOcv4MbVJpNhLX8EwVreM1VP6CkFjP3AunkRpNFHIC4iMkxUixRBlWUS20E1itdnj2eMhXxms9CiWBh8J86waSoHQV3AFZYER179YhEhw3cUy6YhAZFpBmJoAkmbjmL6CoZ4qY/W0aQiiErE2v8z597mWNOsQjIhrflOfIhz1mc2e9vhsIo9nmagmroKzEAiFIG8wEzpFUNMsAFMxLoagYLsGMZURJy8NbC12YwvNYG/epQ2gkPmg6ln3tsAo3UGzR3LauwWAaHku17OGUNVpxZTrCW8MPKuOQiDZJJ0F67KG8wFOuTTIKoenkRWEaWDRQpfrIYA0vripuvO2un2iUOgOc/mJLDDt00Y95yDPTFuPeGPgWXX0RoIWDLXEOomJZhCh0Y0MMnkIE7jMolzRNHQpqRoUMCilYG6nf15T4PF0iqo3Bh7PsiLqJvLJxPVJEBEkqyvIm+MoYOKqKz7LxOyaqaHz4wT0lvY6POsPbww8nmXEVEbQqIvQFBsTdv4ljK1SUrcSMNUxt1IohLNkrufDZwp5FoqyNhVJO8EbA8+qIhw5hJkaIRnagwYRUdY6U03gehnnzeiz1pO21NcQlWsrfufxLCKqUE98ANnjWTIKz9yPnRzHdPdiuzc4l09Sw4RFTFyBtJk6qmEpUxhV30TGs6z4mIHHs4RIWifefnVDtjoYPepWAFnvAImrAGhYAGtdG0q12GLPSg7bc4HijYHHs0TM6iFcKCO1SRBxBgBAdbpsdJuWlp4LhyfPxuwdWP5KkPUcM/DfJs+qw5YH0ajo9ITiKmnvFtK+rdN38obgguVT+0e4nBPz77gEqEJitaPH+SAiG0TkbhF5IvvZtoGGiDydNb5/QETuW+jxrfhvlGdFCEaPNprbtEPiGlKfcplCC+gxfK5zetYHr790ED386Ip9/jL1QH4X8AVVvQz4QvZ6Lm5S1etU9YZneTzg3USeFcLUJjC1CSSNseV+koGdzffqk05uOggx9Um3MSpxljJ9hdn3LzULReOKyOqmRGm5LmKJqd1+KxOHT6HW0r11iMJAD8GP/9JKD2vZGalZBovT/+7p1a9ckbGoslxyFDfjmoMBfAj4EvBrS3n8UnY6+yDwo8BJVb0627YB+BtgD67T2U+p6kj23ruBtwIp8Auq+tmlGptn5dEzx9Ckjh0dJnjOjYQjh0h7NrmeB3EVqU+5/TJNIokrbGCEpLBz2nkkrVNOY4irVEsbMi2jte/TffSnbyatp3z5K88QiFAwwu7eArvv+jrbXngFPfsuJ9y0AzO4haRvK7WuIU5MJq6QDifnUU0sPUXj+lWvIfJe2p8/MMIbnzM0yxDkhGcPT7uJWA4WGDPY2Oq6AW5T1ds6PHaLqh4DUNVjIrJ5ziHB50REgf/Vcv5Oj2+wlCuDvwT+J/BXLdvypct7ReRd2etfE5ErgVuAq4DtwOdFZJ9vfbl+SZ7zMgCig/dBXEOLoZOXqIw0VgNa6EKDiGDsOOnIScym3Y3jY4XJ2GII6SsUMHGVwGRtMdMawcQp1ISzYw2rnG+/8lX0X9THie+d4qnhCqnC01N1+iPDo+M13vjECEnle2ydrFIeeprSzl2YH7yZNJPiriUQGEWENWkIwKnOTiXKqak6Z+uWgTarQQAZPQ7LbQwWllp6eobrZhoi8nmg3T/ory9gSC9W1aPZZH+3iHxfVb+8gOMbLGXbyy+LyJ4Zm+dautwM3KGqNeApEdkP3Ah8Y6nG51kdxBfdQHTycUx13PURiMqueCwqYcv9YBOnSBqVSfu3OaG6sWMwuNv1K0iqSD2BqERB805mTuJ6rRmCh/7Vj7Lxyk0UB3t5zr/eSPjJRzh4bJyj1aaG/t/ef4yeB09Q+KcnAfi1D7yFAk7yG1z93Z6+iKfHYgxrt8S6KxT+w/XbADCTw9juoVn7xBfNOc8uKYuVTaSqPzLXeyJyQkS2ZXf124CTc5zjaPbzpIh8HDdvfhno6PhWljtmMNfSZQfwzZb9DmfbZiEibwPeBrBr164lHKpnuYg37yN66h7skf2E2/aQbLwE4mqz37AJscVewpFDSL2CjhwnqrtAcbLxkmblsU0RjTHVMdL+7Ys6Rklqcxa4SVwhOn0AO3q6seKZD737z6mfOslXf+sutlyzif6LBrnyra+lcMnV2KkxTn7ubqbOVDhTtxSMMBgFbC+FnKm7xXI5MJQD4W/+0+28+R+fQ9f2S5ks7WW8lvCkQm/RtI2vrEXaGYKVwqmWLkvM4C7gLcB7s5//MHMHEekGjKqOZ89fBfx2p8fPZLUEkNvdwrQ1v5lP7DaA519//dp3DnsaSKFEOnwcNl6C1MYJ44pTI508i6lVML0bwCZQLKHFbmypFw0i516KK6SlPmqJpSde/LRDiacwUyOzVhuSZFlPJkQGtxIdvA+dHEOTmPTaV7c9l7n/H3ngTz5B11CZS197KeWhfsqbBwg372wYk6Hh41z00sPszoKVmirVsRpJJWHk6VFO1RIqqXLZZRvQyiSilk1dIRf3F6b1l16rJMrqvA5dtjqD9wJ3ishbgUPATwKIyHbg/ar6OmAL8HFxLVtD4HZV/adzHX8ultsYzLV0OQy03ubvBHyrqQuJrgGMTUlHThGMn8BGXUhSRSpjaLEbKfVii92uoU3XIKYyCkAwfhJEMFNn0Y1FymJINuxZ9OHZ8iBh9RBmcphw7Dj1bVcBoGERMzmMZmMLkippvYpOjhE9+XUA4r0vmn6u63+Ma//LOPUjBwm7ygRDWyEsEG6/hFx0o3jVD7Lv34HpG0Krk9jJMSYPHWH0ySP07eqDrzzDoamYK255Aea5r6Des4n+rC/0emBVGgLA6vTueEuFqg4Dr2iz/Sjwuuz5AeC5Czn+XCy3MZhr6XIXcLuI/AEugHwZcO8yj82zgtiuAWTkOHZyDH3yQaKdexvbtdANNnWFZmoJJoeRkSPYyTFOfvpTjDx+mIFLtzP0/GsIhrYiJuALgy+hKwo4NFpld3+JF2zvXM10LpJBF8COiz2Y2jjB8cdIh4+RTo4TXPtyZ6RqE+jkGOOPPMSxD/wt+//5IK996j4e/embGbh0C7v//VupPXof5nVvpzx+grR3S2MJXG/5rPr2a5Dt1zTeM8DQ0e/R/cCXOfbl+5hIDhKIMPCil7nfyfgJbKnPZWOtYSYTpcekq1Z00De3eRaIyEdxweKNInIY+E3mWLqo6sMicifwCJAAb/eZRBcOFsGEBezwMezoMEyOEW7egQ5sw5Z6mSoMUDJKeOZp183MhNS++1WGH3ySI/cc4F/uOcob/pVhYN8uTN8QwaVXsTMo0V8MiALh2HiN931zhLc8bztDdhSNus5L4E7DIhoWCZI6Wpnk6Ge+gPnClyn0dlHo6yKerFIa6GXrjZez7YVXULv9Vp7zs6+lsO95SLEL87q3AzTaa3ZKffs1BM88QdRVZjS2TCSW9NpXuwK+tI4UutZ8Um2PxKgpzL/jSuGF6haOqv70HG+1Xbqo6q3ArUs1Hs/qxWiKxDWnQxRGkMTYoYtI+7ZSTZVCYAhGj8CpZ7BTY2gcM/rkEcYOneSL3zzKdUNldr/6Byg//yYo9RIPXcJl2bkLQYlNXQX2DHQBTvFUv3Yn3/mdj/C8T/3T3IPqhJ4hTO84ai2jT51k6vQUQSFAU2Xg4kH2/NSPodZiunuRQgkp9TTcS8+WoHeAqK+L1tiw03JyAn4Ue+c+eA0g1XF0FQWMZ7KetYlW51rMsy4Zj5XI0GhFaWrjLmMoSwMNL7ue4FKLxDWmercSJ0pv5ST65LepP/0oT33SZRpXR6p895tH2FQMefMvvoRNL7gWefXPE7f5zIGCcXnqve5fXSnDy97M87v7mPrweyi9+TdmHXP8t36erb/1Z4DrryBnDhPvfRHytTvQJIaXvRmAeMvlsOVydlz7akbe98t8/877qI5U6dvZB0A6Oky0ex+y8wqSgZ0sxlLXDGyi/5Kd7N7SDScm+dju5/PGQ992ek7rIGawmjKH5sIbA4+nA/TuPycdP0v4xl9h/8/9a8JyxJ7/76MARAaOTiSUQ0GBQLo4drLCzt4C/aWAY8WdbOwKKU4N03XiUahNcvxv/poH/+KbhOWQU6crVFJLfxTwQ//majY9b1/bybwT7A03U7rhZqKTjxNv3tfYXrv9Vra8+PmE3/8Xkue8jPo3/pGH/vzTnHlihE1XbqR7Szd7MmPQyuA7f48XvhPsXe9j9MBhxp46xoGP/TMbrniSoZeME+6zjZjD+ZB2D1HY9zw2X/Nl4CT1iaYJXOuGYC2gColvbuPxzM/TH/88W2+8kuTD72HrC55DPFnl2G+8jUJvF2m1zr63/zJKyWkPmZCdlaPY8i4YS7lo7DjJY0c4+42vcfrBA9Qn69z5yf2kqvSFAdcNldn7nCEued11DN38b8/b5QJMMwQAxTf9OuHwAZKhS9zrG1/NDdf9EPXt1zT2sXe9Dyl3I6/832edz7zhnQwCuUTkN170cnYeOcWOlx+ieNWNszKLFkrav52oNsHm512Crac8c89R7n35K7jxS1+Ytp+9632EO/dir/+x8/o8z3RUwfqVgcfTnlb3RKGvi64rrnHZQGLQ2hRSKIEY4mee4InfeDfHHzjB4WfG2Lqpi2JfkY2XDxF1F7Gp5YlPPsFIPWE0dndf20shPaFhS2+BK/7NdWy+/jkUXvwG6tlkvRQkLeeON146+3ptyoG/+Sz7Nu3AXve6c57rhV//EuC0hgYv/z67/o+NswzQQok372PoJS+lNNDL0W8fpzRY4tGfvpkrPtqsK6qdGcX0DlM49vCiGM2lIDrxmHO1rSkUXacrMG8MPOdNq3ti13//f7HF3mlpkjkFMfTv3cHkiUnSQ2M8dnQCjk7Q/+QIPaHhTD2lkiqpQt0q5UDojwzX/vAedrzkSja87iep73guyfJdWluCH/8lLv1xWIiz4IqP/gPDv/PO8zYEjTHsvoJyZZK+nb0Mf3+Yvp19SFxppJaWf/Y3GXnfL3PiA5/g8je/hmDL7nkN1/my0Ml97RkCh67TlcH6qFf3rAo+tX8EW+x1VcNt0CCie9sQ/RcNsnNrDzvKITvKIRvKEYUooGCEciBsKBi3vRBw6bWbufxn39AwBGuZoV9936Kdy548SHLqCGE55PREnZP7z/D5q14+bZ/Bd/4ePdv6efJvP8+jf/h+Hv3pmxft89uxVif3BZG5iTp5rDX8ysCzaLz+Uucpn9nGMifevI/im36dS154Dzte/g3OPHKAyWPDpNWY+mSdTZMxhe6I7s29lDcP0rNjE6U9e+Glb2q70lhpwpFDixIUXijRU/fwzN/9PaNPnWTi2ARn6pYzdUt/NHutsvN3PwRfud0V7WWMx0pvtLglvkFWQLfeURbUa2lN4Y2BZ0mIDt43p6pkfPELMBe/gI1vgG1Pfh2tVV3KJoAxBBu2TgvYrlZWxBCcfJz48JOc/t5hTj1ymsfPVAEoGKEnNNyx83puOXz/9INe+iaC0aMN8b7FNgSw8AK6NYtCujzNbZYdbww8i46KdCwv3C67Zq2Vni9nkxVJak4CI7bYVNlRDsnKNghE2NLrVmXBQ3czds9X6H6rE7FcbBXXCxf1MQOPp1MutHz3ZGAn3ztVXZbP0rBItHsfGy4dYseN2ygXAiqpUkm1IW998Bf/LSc++Y88/OGv8be7rud9G69dlrFdCDg3kXb0WGt4Y+DxLALXbFqezsvx5n2Yi65i7y/9Mlf8x1u47udeyJU/tIuLLh9iQyEA4LF/fJxCbxddQ2WOVBLqVvntgauXZXzrHgWr2tFjreHdRB7PGqNRB7HjuQy9+BY23fcPpMPHGN//NCOPP0P1bIUH/+LrxJN1dpRDJhIlXYOT02plLd71d4I3Bh7PCtAa0D1f7A03E509zMarh+nb/yDpyEkG9z7J5PFhBg6OMnliCrWW47/18+z4z/+NtH87wfgJtNDV7Cbn6RhvDDwez6Kx2AHdZGAnDOwk7BokOnuMTd19DE2OMXTkFJPHzyCBYWDvDtIHvwQvfVMj+8fUJ7GF8+/1cKGgqj6byOPxrH6Swd0wuJtocAdSm6RsEwaO7Cc+vD/bYXrFhjcEC2e91hn4ALLHsw5JBnYSb7mc+rarCHZfQeGSqwn6hwi3XQyA3POxtscVjj28nMNcc+gyVSCLyAYRuVtEnsh+DrbZ53IReaDlMSYi78ze+y0ROdLy3rxaJCuyMhCRp4FxXEp5oqo3iMgG4G+APcDTwE+p6shKjM/jWU/Em/dBpomko661uL7gjW33zUXtLIJBCSZOkfZsWp6BrhGWKWbwLuALqvpeEXlX9vrXpo1D9THgOgARCYAjwMdbdvlDVf29Tj9wJVcGN6nqdaqaVyflF38Z8IXstcfjWUQ6jVWYrIFm2rMJsSstDbiK0GWrM7gZ+FD2/EPAj8+z/yuAJ1X14LP9wNXkJlroxXs8ng6ZSzywlWrafgJbrc3pVwLFBZA7eeD6v9/X8njbAj5qi6oeA8h+bp5n/1uAj87Y9g4ReVBEPtjOzTSTlTIGCnxORL7d8gvq6OJF5G35L/fU6dPLNFyPZ22jQQGLEIyfaGw7MZViaeoU5e1IPedgYSuD06p6Q8vjttZTicjnReShNo8FycuKSAF4A/C3LZv/FNiLcyMdA35/vvOslMl/saoeFZHNwN0i8v1OD8x+obcBPP/669dnwq/HswQYdJqg3JauAPBfoYWyWPLUqvojc70nIidEZJuqHhORbcDJc5zqtcD9qtqw9K3PReTPgU/ON54VWRmo6tHs50lcwONG4ER20XRw8R6Px7MiqGpHj/PkLuAt2fO3AP9wjn1/mhkuonwuzfgJ4KH5PnDZjYGIdItIb/4ceBVuoAu5eI/H41l2VDtzES1CAPm9wCtF5AngldlrRGS7iHw630lEurL3Z+YK/46IfE9EHgRuAv7zfB+4Em6iLcDHRST//NtV9Z9E5FvAnSLyVuAQ8JMrMDaPx+M5J8vRxUxVh3EZQjO3HwVe1/J6Chhqs9+bF/qZy24MVPUAMKt/4VwX7/F4PKsGVWyyGvvunT8+Z8zj8Xg6RFHUrrX2S53hjYHH4/F0ioKm3hh4PB7PBY5fGXg8Ho9HvTHweDweD3hj4PF4PBc66rOJPB6PxwOK9SsDj8fjucDxMQOPx+PxKD5m4PF4PB5VX2fg8Xg8Fzw+gOzxeDweX3Tm8Xg8nixmYFd6GEuCNwYej8fTKT6byOPxeDzgs4k8Ho/Ho+u36GxFeiCfCxF5jYg8JiL7ReRdKz0ej8fjyVFVbFzv6LHWWFUrAxEJgD/G9fQ8DHxLRO5S1UdWdmQej8cDPpto+bgR2J+1xkRE7gBuBrwx8Hg8qwJvDJaHHcAzLa8PAy9YobF4PB7PdHw20bIhbbbptB1E3ga8LXs5Ue7qGgZOL/XAViEb8dd9IeGv+/y56HxPoJXhz8YP/MXGDndfU3+v1WYMDgO7Wl7vBI627qCqtwG35a9F5D5VvWF5hrd68Nd9YeGve3Wgqq9Z6TEsFastm+hbwGUicrGIFIBbgLtWeEwej8ez7llVKwNVTUTkHcBngQD4oKo+vMLD8ng8nnXPqjIGAKr6aeDTCzjktvl3WZf4676w8NftWVJEVeffy+PxeDzrmtUWM/B4PB7PCuCNgcfj8XjWrjFYzxpGIvJBETkpIg+1bNsgIneLyBPZz8GW996d/R4eE5FXr8yozx8R2SUi/ywij4rIwyLyi9n2dX3tIlISkXtF5LvZdf/3bPu6vu4cEQlE5Dsi8sns9QVx3asOVV1zD1ym0ZPAJUAB+C5w5UqPaxGv74eA64GHWrb9DvCu7Pm7gP8ne35ldv1F4OLs9xKs9DU8y+veBlyfPe8FHs+ub11fO67Ysid7HgH3AD+43q+75fr/C3A78Mns9QVx3avtsVZXBg0NI1WtA7mG0bpAVb8MnJmx+WbgQ9nzDwE/3rL9DlWtqepTwH7c72fNoarHVPX+7Pk48ChOomRdX7s6JrKXUfZQ1vl1A4jITuD1wPtbNq/7616NrFVj0E7DaMcKjWW52KKqx8BNmsDmbPu6/F2IyB7gebi75HV/7Zmr5AHgJHC3ql4Q1w28D/hVoLWX5IVw3auOtWoM5tUwuoBYd78LEekB/h54p6qOnWvXNtvW5LWraqqq1+EkWG4UkavPsfu6uG4R+VHgpKp+u9ND2mxbc9e9WlmrxmBeDaN1yAkR2QaQ/TyZbV9XvwsRiXCG4COq+rFs8wVx7QCqehb4EvAa1v91vxh4g4g8jXP1/rCI/DXr/7pXJWvVGFyIGkZ3AW/Jnr8F+IeW7beISFFELgYuA+5dgfGdNyIiwAeAR1X1D1reWtfXLiKbRGQge14GfgT4Puv8ulX13aq6U1X34L7DX1TVf8c6v+7VyqqTo+gEXecaRiLyUeDlwEYROQz8JvBe4E4ReStwCPhJAFV9WETuxDUASoC3q+paFVx/MfBm4HuZ/xzgv7L+r30b8KGs058B7lTVT4rIN1jf1z0X6/3vvSrxchQej8fjWbNuIo/H4/EsIt4YeDwej8cbA4/H4/F4Y+DxeDwevDHweDweD94YeDwejwdvDDwej8eDNwaeNY6I/ICIPJj1BOjO+gE8V0R+T0S+l733n1Z6nB7PascXnXnWPCLyP4ASUMbp15zFSTr8m6xafYOqzpQE93g8LXhj4FnzZPpU3wKqwIuAO4E/U9W7V3RgHs8aYk1qE3k8M9gA9OCawpRwUsf+LsfjWQB+ZeBZ84jIXTgJ5Itxom8P4dxEt3g3kcfTGT6A7FnTiMjPAImq3o5Tu/wB4ABO7fJBEfku8KYVHKLHsybwKwOPx+Px+JWBx+PxeLwx8Hg8Hg/eGHg8Ho8Hbww8Ho/HgzcGHo/H48EbA4/H4/HgjYHH4/F4gP8fCJPbn7kbwPkAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 2 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "(narr.isel(time=-1) - narr.isel(time=0)).plot()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "0585b436",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "0efadc8c",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.8.0"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}