{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {
    "toc": true
   },
   "source": [
    "<h1>Table of Contents<span class=\"tocSkip\"></span></h1>\n",
    "<div class=\"toc\"><ul class=\"toc-item\"></ul></div>"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Crashing on first time step on the western boundary at 2 465 14 in fortran."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Going to try smoothing the restart file."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import matplotlib.pyplot as plt\n",
    "import numpy as np\n",
    "import xarray as xr\n",
    "\n",
    "%matplotlib inline"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "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",
       "<title>Show/Hide data repr</title>\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",
       "<title>Show/Hide attributes</title>\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",
       ".xr-wrap {\n",
       "  min-width: 300px;\n",
       "  max-width: 700px;\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, 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><div class='xr-wrap'><div class='xr-header'><div class='xr-obj-type'>xarray.Dataset</div></div><ul class='xr-sections'><li class='xr-section-item'><input id='section-0f9b7713-3fa4-4b6b-8cc8-67b4b53c1581' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-0f9b7713-3fa4-4b6b-8cc8-67b4b53c1581' class='xr-section-summary'  title='Expand/collapse section'>Dimensions:</label><div class='xr-section-inline-details'><ul class='xr-dim-list'><li><span>t</span>: 1</li><li><span>x</span>: 398</li><li><span>y</span>: 898</li><li><span>z</span>: 40</li></ul></div><div class='xr-section-details'></div></li><li class='xr-section-item'><input id='section-a5401a99-d730-4719-a691-73c6805936f2' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-a5401a99-d730-4719-a691-73c6805936f2' 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-a89a5016-0660-475e-89ac-1b2b811898d5' class='xr-section-summary-in' type='checkbox'  ><label for='section-a89a5016-0660-475e-89ac-1b2b811898d5' class='xr-section-summary' >Data variables: <span>(47)</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>nav_lon</span></div><div class='xr-var-dims'>(y, x)</div><div class='xr-var-dtype'>float32</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-c434187b-4fae-4b49-af7a-ed64d90c8d27' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-c434187b-4fae-4b49-af7a-ed64d90c8d27' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-4376a441-e71a-472c-bb88-7e351a9463f5' class='xr-var-data-in' type='checkbox'><label for='data-4376a441-e71a-472c-bb88-7e351a9463f5' 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><pre class='xr-var-data'>[357404 values with dtype=float32]</pre></li><li class='xr-var-item'><div class='xr-var-name'><span>nav_lat</span></div><div class='xr-var-dims'>(y, x)</div><div class='xr-var-dtype'>float32</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-5db456cf-47cd-4864-bda6-04872bca691f' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-5db456cf-47cd-4864-bda6-04872bca691f' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-847f5794-a70b-4d39-b41b-4ac7d15a9b05' class='xr-var-data-in' type='checkbox'><label for='data-847f5794-a70b-4d39-b41b-4ac7d15a9b05' 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><pre class='xr-var-data'>[357404 values with dtype=float32]</pre></li><li class='xr-var-item'><div class='xr-var-name'><span>nav_lev</span></div><div class='xr-var-dims'>(z)</div><div class='xr-var-dtype'>float32</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-70ab56f5-32cd-40a4-8434-3a98ac943400' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-70ab56f5-32cd-40a4-8434-3a98ac943400' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-b9e89efb-c436-4cd0-8e0d-daff9110a176' class='xr-var-data-in' type='checkbox'><label for='data-b9e89efb-c436-4cd0-8e0d-daff9110a176' 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><pre class='xr-var-data'>array([  0.5     ,   1.500003,   2.500011,   3.500031,   4.500071,   5.500151,\n",
       "         6.50031 ,   7.500623,   8.501236,   9.502433,  10.504766,  11.509312,\n",
       "        12.518167,  13.535412,  14.568982,  15.634288,  16.761173,  18.007135,\n",
       "        19.481785,  21.389978,  24.100256,  28.229916,  34.685757,  44.517723,\n",
       "        58.484333,  76.58559 ,  98.06296 , 121.866516, 147.08946 , 173.11449 ,\n",
       "       199.57304 , 226.2603  , 253.06664 , 279.93454 , 306.8342  , 333.75018 ,\n",
       "       360.67453 , 387.6032  , 414.5341  , 441.4661  ], dtype=float32)</pre></li><li class='xr-var-item'><div class='xr-var-name'><span>time_counter</span></div><div class='xr-var-dims'>(t)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-e96eb2fb-1100-48f4-af31-af39d89d718b' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-e96eb2fb-1100-48f4-af31-af39d89d718b' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-3eae18fa-387c-41e6-a719-9552143cee4b' class='xr-var-data-in' type='checkbox'><label for='data-3eae18fa-387c-41e6-a719-9552143cee4b' 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><pre class='xr-var-data'>array([4914000.])</pre></li><li class='xr-var-item'><div class='xr-var-name'><span>kt</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-4a3d499f-f08a-4f7a-85ca-c120d8f3d8a5' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-4a3d499f-f08a-4f7a-85ca-c120d8f3d8a5' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-65b77b25-9473-4f8c-9bc2-30944eac8802' class='xr-var-data-in' type='checkbox'><label for='data-65b77b25-9473-4f8c-9bc2-30944eac8802' 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><pre class='xr-var-data'>array(4914000.)</pre></li><li class='xr-var-item'><div class='xr-var-name'><span>ndastp</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-413191e1-280b-477e-be48-afdda238b2b1' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-413191e1-280b-477e-be48-afdda238b2b1' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-9a002992-50d4-422b-bede-f4071415388c' class='xr-var-data-in' type='checkbox'><label for='data-9a002992-50d4-422b-bede-f4071415388c' 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><pre class='xr-var-data'>array(20201113.)</pre></li><li class='xr-var-item'><div class='xr-var-name'><span>adatrj</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-52354f19-e479-4cb0-acbe-a60b2f4af7a4' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-52354f19-e479-4cb0-acbe-a60b2f4af7a4' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-2e6d87e8-17b8-4dfd-b9f0-7ff54266d480' class='xr-var-data-in' type='checkbox'><label for='data-2e6d87e8-17b8-4dfd-b9f0-7ff54266d480' 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><pre class='xr-var-data'>array(2275.)</pre></li><li class='xr-var-item'><div class='xr-var-name'><span>ssh_ibb</span></div><div class='xr-var-dims'>(t, y, x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-abf70191-26e5-453b-ac8b-8e3df5114760' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-abf70191-26e5-453b-ac8b-8e3df5114760' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-c6c9f121-138a-4103-af4d-88a3a0cc9143' class='xr-var-data-in' type='checkbox'><label for='data-c6c9f121-138a-4103-af4d-88a3a0cc9143' 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><pre class='xr-var-data'>[357404 values with dtype=float64]</pre></li><li class='xr-var-item'><div class='xr-var-name'><span>rnf_b</span></div><div class='xr-var-dims'>(t, y, x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-77078af8-6243-44fa-a126-b88d010d99d3' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-77078af8-6243-44fa-a126-b88d010d99d3' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-61ebeacd-65ef-46bc-96d7-e261065e9abd' class='xr-var-data-in' type='checkbox'><label for='data-61ebeacd-65ef-46bc-96d7-e261065e9abd' 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><pre class='xr-var-data'>[357404 values with dtype=float64]</pre></li><li class='xr-var-item'><div class='xr-var-name'><span>rnf_hc_b</span></div><div class='xr-var-dims'>(t, y, x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-3c6be837-d64c-4931-914d-6eea905bae8c' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-3c6be837-d64c-4931-914d-6eea905bae8c' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-433d6c8f-051a-40b6-b4e1-d26dd5541a92' class='xr-var-data-in' type='checkbox'><label for='data-433d6c8f-051a-40b6-b4e1-d26dd5541a92' 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><pre class='xr-var-data'>[357404 values with dtype=float64]</pre></li><li class='xr-var-item'><div class='xr-var-name'><span>rnf_sc_b</span></div><div class='xr-var-dims'>(t, y, x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-9b62bd20-25d5-438d-9157-c47928fe93aa' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-9b62bd20-25d5-438d-9157-c47928fe93aa' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-fd52fe12-4e22-4c11-9353-34d635f7226a' class='xr-var-data-in' type='checkbox'><label for='data-fd52fe12-4e22-4c11-9353-34d635f7226a' 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><pre class='xr-var-data'>[357404 values with dtype=float64]</pre></li><li class='xr-var-item'><div class='xr-var-name'><span>utau_b</span></div><div class='xr-var-dims'>(t, y, x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-91d9f1df-2028-4752-ad38-414a72240cbd' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-91d9f1df-2028-4752-ad38-414a72240cbd' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-fe0f7861-0caf-4d11-97a5-32e0ab3b43f0' class='xr-var-data-in' type='checkbox'><label for='data-fe0f7861-0caf-4d11-97a5-32e0ab3b43f0' 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><pre class='xr-var-data'>[357404 values with dtype=float64]</pre></li><li class='xr-var-item'><div class='xr-var-name'><span>vtau_b</span></div><div class='xr-var-dims'>(t, y, x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-22360a25-1075-4e82-aea6-059986585b17' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-22360a25-1075-4e82-aea6-059986585b17' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-251e7a8d-90f4-4313-8701-86649231fb4a' class='xr-var-data-in' type='checkbox'><label for='data-251e7a8d-90f4-4313-8701-86649231fb4a' 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><pre class='xr-var-data'>[357404 values with dtype=float64]</pre></li><li class='xr-var-item'><div class='xr-var-name'><span>qns_b</span></div><div class='xr-var-dims'>(t, y, x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-9cb2ec68-e0b3-4d25-adfe-dfd2ed0e852b' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-9cb2ec68-e0b3-4d25-adfe-dfd2ed0e852b' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-119e1f4b-e7de-4162-92c0-dd62e996d10b' class='xr-var-data-in' type='checkbox'><label for='data-119e1f4b-e7de-4162-92c0-dd62e996d10b' 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><pre class='xr-var-data'>[357404 values with dtype=float64]</pre></li><li class='xr-var-item'><div class='xr-var-name'><span>emp_b</span></div><div class='xr-var-dims'>(t, y, x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-b187f363-5bbb-470d-ba7b-b1321fee7d13' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-b187f363-5bbb-470d-ba7b-b1321fee7d13' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-de6b25f5-e709-4b19-b8b3-647ce9d0e3c9' class='xr-var-data-in' type='checkbox'><label for='data-de6b25f5-e709-4b19-b8b3-647ce9d0e3c9' 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><pre class='xr-var-data'>[357404 values with dtype=float64]</pre></li><li class='xr-var-item'><div class='xr-var-name'><span>sfx_b</span></div><div class='xr-var-dims'>(t, y, x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-12026a6f-dcc6-472c-be11-ca6527b2a0e0' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-12026a6f-dcc6-472c-be11-ca6527b2a0e0' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-45cddfa1-f0c8-444f-87a6-a22a69715763' class='xr-var-data-in' type='checkbox'><label for='data-45cddfa1-f0c8-444f-87a6-a22a69715763' 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><pre class='xr-var-data'>[357404 values with dtype=float64]</pre></li><li class='xr-var-item'><div class='xr-var-name'><span>en</span></div><div class='xr-var-dims'>(t, z, y, x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-89839907-b3a6-4bed-9fc9-a711be321a45' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-89839907-b3a6-4bed-9fc9-a711be321a45' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-0338a392-7a87-4606-9f23-f8eb8b6146e8' class='xr-var-data-in' type='checkbox'><label for='data-0338a392-7a87-4606-9f23-f8eb8b6146e8' 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><pre class='xr-var-data'>[14296160 values with dtype=float64]</pre></li><li class='xr-var-item'><div class='xr-var-name'><span>avt</span></div><div class='xr-var-dims'>(t, z, y, x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-3762b66a-92ff-4502-84f9-043e65026d4b' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-3762b66a-92ff-4502-84f9-043e65026d4b' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-ff16ada7-2bb5-4c0d-947e-97acb96a410e' class='xr-var-data-in' type='checkbox'><label for='data-ff16ada7-2bb5-4c0d-947e-97acb96a410e' 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><pre class='xr-var-data'>[14296160 values with dtype=float64]</pre></li><li class='xr-var-item'><div class='xr-var-name'><span>avm</span></div><div class='xr-var-dims'>(t, z, y, x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-aed9e6fa-3d63-44ff-933e-0c54bb960f61' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-aed9e6fa-3d63-44ff-933e-0c54bb960f61' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-5a661199-59a2-4b49-b821-fafda121e4ae' class='xr-var-data-in' type='checkbox'><label for='data-5a661199-59a2-4b49-b821-fafda121e4ae' 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><pre class='xr-var-data'>[14296160 values with dtype=float64]</pre></li><li class='xr-var-item'><div class='xr-var-name'><span>avmu</span></div><div class='xr-var-dims'>(t, z, y, x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-d5edf48b-c533-4893-809b-538bd226d3b9' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-d5edf48b-c533-4893-809b-538bd226d3b9' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-bd8be8d3-5ce2-4dc8-bfed-3aa680a8ae31' class='xr-var-data-in' type='checkbox'><label for='data-bd8be8d3-5ce2-4dc8-bfed-3aa680a8ae31' 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><pre class='xr-var-data'>[14296160 values with dtype=float64]</pre></li><li class='xr-var-item'><div class='xr-var-name'><span>avmv</span></div><div class='xr-var-dims'>(t, z, y, x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-2444f801-f8c4-459b-8a10-189543c46ac4' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-2444f801-f8c4-459b-8a10-189543c46ac4' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-4e265b94-6495-438e-9739-d697d305e6e1' class='xr-var-data-in' type='checkbox'><label for='data-4e265b94-6495-438e-9739-d697d305e6e1' 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><pre class='xr-var-data'>[14296160 values with dtype=float64]</pre></li><li class='xr-var-item'><div class='xr-var-name'><span>mxln</span></div><div class='xr-var-dims'>(t, z, y, x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-426d2239-9535-4939-9076-81501bb6f412' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-426d2239-9535-4939-9076-81501bb6f412' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-8b35dcd1-6026-41e1-93c8-c9938787ae6a' class='xr-var-data-in' type='checkbox'><label for='data-8b35dcd1-6026-41e1-93c8-c9938787ae6a' 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><pre class='xr-var-data'>[14296160 values with dtype=float64]</pre></li><li class='xr-var-item'><div class='xr-var-name'><span>ub2_b</span></div><div class='xr-var-dims'>(t, y, x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-55bfca38-884b-459b-b5ac-f65c94914f42' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-55bfca38-884b-459b-b5ac-f65c94914f42' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-348fa61a-40d6-4e5e-9667-22f02a42db9c' class='xr-var-data-in' type='checkbox'><label for='data-348fa61a-40d6-4e5e-9667-22f02a42db9c' 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><pre class='xr-var-data'>[357404 values with dtype=float64]</pre></li><li class='xr-var-item'><div class='xr-var-name'><span>vb2_b</span></div><div class='xr-var-dims'>(t, y, x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-ee0a2a48-3e6b-4500-bde1-af192ebcdb5b' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-ee0a2a48-3e6b-4500-bde1-af192ebcdb5b' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-f78fe208-0280-4b88-947a-a7b093ad87f7' class='xr-var-data-in' type='checkbox'><label for='data-f78fe208-0280-4b88-947a-a7b093ad87f7' 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><pre class='xr-var-data'>[357404 values with dtype=float64]</pre></li><li class='xr-var-item'><div class='xr-var-name'><span>sbc_hc_b</span></div><div class='xr-var-dims'>(t, y, x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-3b115916-df01-4daa-bc28-5e0716b0324d' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-3b115916-df01-4daa-bc28-5e0716b0324d' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-43f59da7-3c64-466d-84c6-73628e6e88a9' class='xr-var-data-in' type='checkbox'><label for='data-43f59da7-3c64-466d-84c6-73628e6e88a9' 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><pre class='xr-var-data'>[357404 values with dtype=float64]</pre></li><li class='xr-var-item'><div class='xr-var-name'><span>sbc_sc_b</span></div><div class='xr-var-dims'>(t, y, x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-a4d50068-40f9-42df-bbe0-fc9d5d98da23' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-a4d50068-40f9-42df-bbe0-fc9d5d98da23' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-ef354aaf-1001-47d0-9b29-a40958044e9f' class='xr-var-data-in' type='checkbox'><label for='data-ef354aaf-1001-47d0-9b29-a40958044e9f' 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><pre class='xr-var-data'>[357404 values with dtype=float64]</pre></li><li class='xr-var-item'><div class='xr-var-name'><span>qsr_hc_b</span></div><div class='xr-var-dims'>(t, z, y, x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-d02b3a35-24a0-4e82-b11b-c9eb8448b21f' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-d02b3a35-24a0-4e82-b11b-c9eb8448b21f' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-e3132e27-17f8-4e33-ad5f-f0f722b4e479' class='xr-var-data-in' type='checkbox'><label for='data-e3132e27-17f8-4e33-ad5f-f0f722b4e479' 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><pre class='xr-var-data'>[14296160 values with dtype=float64]</pre></li><li class='xr-var-item'><div class='xr-var-name'><span>fraqsr_1lev</span></div><div class='xr-var-dims'>(t, y, x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-49d90f72-b484-4083-b825-f65137f56fdd' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-49d90f72-b484-4083-b825-f65137f56fdd' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-b5134e7a-b5b8-4c35-9a2c-891e5c406d4d' class='xr-var-data-in' type='checkbox'><label for='data-b5134e7a-b5b8-4c35-9a2c-891e5c406d4d' 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><pre class='xr-var-data'>[357404 values with dtype=float64]</pre></li><li class='xr-var-item'><div class='xr-var-name'><span>fse3t_b</span></div><div class='xr-var-dims'>(t, z, y, x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-3d332100-ee82-403e-9de4-d6ff90f147ea' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-3d332100-ee82-403e-9de4-d6ff90f147ea' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-b6496c67-7052-47bd-98e1-aa433ac3e0da' class='xr-var-data-in' type='checkbox'><label for='data-b6496c67-7052-47bd-98e1-aa433ac3e0da' 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><pre class='xr-var-data'>[14296160 values with dtype=float64]</pre></li><li class='xr-var-item'><div class='xr-var-name'><span>fse3t_n</span></div><div class='xr-var-dims'>(t, z, y, x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-daa9e2af-1edb-4b8b-aa57-81ce6b60e3a5' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-daa9e2af-1edb-4b8b-aa57-81ce6b60e3a5' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-958c0967-95d6-4f28-bc20-33107b7ee54c' class='xr-var-data-in' type='checkbox'><label for='data-958c0967-95d6-4f28-bc20-33107b7ee54c' 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><pre class='xr-var-data'>[14296160 values with dtype=float64]</pre></li><li class='xr-var-item'><div class='xr-var-name'><span>rdt</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-4614921e-1492-443a-83a4-7ef7b03ea2ae' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-4614921e-1492-443a-83a4-7ef7b03ea2ae' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-da47a7da-c23f-49a7-8f04-19e3a0ac948f' class='xr-var-data-in' type='checkbox'><label for='data-da47a7da-c23f-49a7-8f04-19e3a0ac948f' 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><pre class='xr-var-data'>array(40.)</pre></li><li class='xr-var-item'><div class='xr-var-name'><span>rdttra1</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-79fc410b-46bc-49af-b49e-4ef399af74d4' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-79fc410b-46bc-49af-b49e-4ef399af74d4' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-0124b9d6-f988-4e33-b621-0224460ff4b8' class='xr-var-data-in' type='checkbox'><label for='data-0124b9d6-f988-4e33-b621-0224460ff4b8' 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><pre class='xr-var-data'>array(40.)</pre></li><li class='xr-var-item'><div class='xr-var-name'><span>ub</span></div><div class='xr-var-dims'>(t, z, y, x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-00ed2d75-e107-4e9f-8e0f-c82bb4736f60' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-00ed2d75-e107-4e9f-8e0f-c82bb4736f60' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-268d30d5-43d9-40b0-8c83-4dde05934872' class='xr-var-data-in' type='checkbox'><label for='data-268d30d5-43d9-40b0-8c83-4dde05934872' 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><pre class='xr-var-data'>[14296160 values with dtype=float64]</pre></li><li class='xr-var-item'><div class='xr-var-name'><span>vb</span></div><div class='xr-var-dims'>(t, z, y, x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-6b974821-a9b0-4cde-bf67-0b233b7164cb' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-6b974821-a9b0-4cde-bf67-0b233b7164cb' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-51a1c82c-29c1-4cd3-8df9-6f07e9a7277c' class='xr-var-data-in' type='checkbox'><label for='data-51a1c82c-29c1-4cd3-8df9-6f07e9a7277c' 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><pre class='xr-var-data'>[14296160 values with dtype=float64]</pre></li><li class='xr-var-item'><div class='xr-var-name'><span>tb</span></div><div class='xr-var-dims'>(t, z, y, x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-c67d56f1-3688-40a7-8d5c-6e6aba3b4f60' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-c67d56f1-3688-40a7-8d5c-6e6aba3b4f60' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-736f97df-a0db-4287-b09b-bca9c7ed522d' class='xr-var-data-in' type='checkbox'><label for='data-736f97df-a0db-4287-b09b-bca9c7ed522d' 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><pre class='xr-var-data'>[14296160 values with dtype=float64]</pre></li><li class='xr-var-item'><div class='xr-var-name'><span>sb</span></div><div class='xr-var-dims'>(t, z, y, x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-37fa1654-2a59-46ac-a6a6-cf3161b9a7d5' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-37fa1654-2a59-46ac-a6a6-cf3161b9a7d5' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-3809028c-5871-4295-abba-f0397926274e' class='xr-var-data-in' type='checkbox'><label for='data-3809028c-5871-4295-abba-f0397926274e' 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><pre class='xr-var-data'>[14296160 values with dtype=float64]</pre></li><li class='xr-var-item'><div class='xr-var-name'><span>rotb</span></div><div class='xr-var-dims'>(t, z, y, x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-65c32baf-e173-4b31-98b6-c90a352fa32f' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-65c32baf-e173-4b31-98b6-c90a352fa32f' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-86c43ac0-d396-4fdd-9c65-fb66eac9a70a' class='xr-var-data-in' type='checkbox'><label for='data-86c43ac0-d396-4fdd-9c65-fb66eac9a70a' 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><pre class='xr-var-data'>[14296160 values with dtype=float64]</pre></li><li class='xr-var-item'><div class='xr-var-name'><span>hdivb</span></div><div class='xr-var-dims'>(t, z, y, x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-18f39ad1-2b1a-49d9-84f8-1cc95ce5a2b5' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-18f39ad1-2b1a-49d9-84f8-1cc95ce5a2b5' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-f89a89a9-f778-443c-a566-034f4d4da10b' class='xr-var-data-in' type='checkbox'><label for='data-f89a89a9-f778-443c-a566-034f4d4da10b' 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><pre class='xr-var-data'>[14296160 values with dtype=float64]</pre></li><li class='xr-var-item'><div class='xr-var-name'><span>sshb</span></div><div class='xr-var-dims'>(t, y, x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-80028829-40aa-4674-a7dd-a20e3d7f1372' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-80028829-40aa-4674-a7dd-a20e3d7f1372' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-7b86097d-f925-4bb9-b797-914b2af7a03a' class='xr-var-data-in' type='checkbox'><label for='data-7b86097d-f925-4bb9-b797-914b2af7a03a' 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><pre class='xr-var-data'>[357404 values with dtype=float64]</pre></li><li class='xr-var-item'><div class='xr-var-name'><span>un</span></div><div class='xr-var-dims'>(t, z, y, x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-bdfa95fa-893f-4f22-b26e-f65ee01e6b7e' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-bdfa95fa-893f-4f22-b26e-f65ee01e6b7e' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-2e009fbb-aa0b-4ccc-82fc-9794142eda61' class='xr-var-data-in' type='checkbox'><label for='data-2e009fbb-aa0b-4ccc-82fc-9794142eda61' 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><pre class='xr-var-data'>[14296160 values with dtype=float64]</pre></li><li class='xr-var-item'><div class='xr-var-name'><span>vn</span></div><div class='xr-var-dims'>(t, z, y, x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-4489c3c7-fdb2-4cf1-8a69-fc9c76953201' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-4489c3c7-fdb2-4cf1-8a69-fc9c76953201' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-63aec75e-21c7-4d35-abe5-056aa46c3b7c' class='xr-var-data-in' type='checkbox'><label for='data-63aec75e-21c7-4d35-abe5-056aa46c3b7c' 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><pre class='xr-var-data'>[14296160 values with dtype=float64]</pre></li><li class='xr-var-item'><div class='xr-var-name'><span>tn</span></div><div class='xr-var-dims'>(t, z, y, x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-cd83abf4-3643-4395-9977-71268bbf6351' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-cd83abf4-3643-4395-9977-71268bbf6351' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-16e10f63-5ef7-4149-b9a2-448779dfe3ef' class='xr-var-data-in' type='checkbox'><label for='data-16e10f63-5ef7-4149-b9a2-448779dfe3ef' 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><pre class='xr-var-data'>[14296160 values with dtype=float64]</pre></li><li class='xr-var-item'><div class='xr-var-name'><span>sn</span></div><div class='xr-var-dims'>(t, z, y, x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-1bd2b6b1-6bd2-43bc-848b-ac5da19c4167' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-1bd2b6b1-6bd2-43bc-848b-ac5da19c4167' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-68514823-2afc-4a12-b04d-6c0fbb07bba6' class='xr-var-data-in' type='checkbox'><label for='data-68514823-2afc-4a12-b04d-6c0fbb07bba6' 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><pre class='xr-var-data'>[14296160 values with dtype=float64]</pre></li><li class='xr-var-item'><div class='xr-var-name'><span>rotn</span></div><div class='xr-var-dims'>(t, z, y, x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-6d44481f-a712-42ef-aa0a-16071a74a903' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-6d44481f-a712-42ef-aa0a-16071a74a903' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-cdb411e1-ec14-477c-a531-c6ef59416181' class='xr-var-data-in' type='checkbox'><label for='data-cdb411e1-ec14-477c-a531-c6ef59416181' 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><pre class='xr-var-data'>[14296160 values with dtype=float64]</pre></li><li class='xr-var-item'><div class='xr-var-name'><span>hdivn</span></div><div class='xr-var-dims'>(t, z, y, x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-5a3f5a26-864d-4f30-b8e7-dec5d63dc971' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-5a3f5a26-864d-4f30-b8e7-dec5d63dc971' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-b61dc266-6a56-4d47-82b2-b724230af7ba' class='xr-var-data-in' type='checkbox'><label for='data-b61dc266-6a56-4d47-82b2-b724230af7ba' 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><pre class='xr-var-data'>[14296160 values with dtype=float64]</pre></li><li class='xr-var-item'><div class='xr-var-name'><span>sshn</span></div><div class='xr-var-dims'>(t, y, x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-b5004a74-b52b-401c-bf09-f88457a0c08c' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-b5004a74-b52b-401c-bf09-f88457a0c08c' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-ecd0cd83-fc65-4b3e-b273-ae0cb7ec053f' class='xr-var-data-in' type='checkbox'><label for='data-ecd0cd83-fc65-4b3e-b273-ae0cb7ec053f' 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><pre class='xr-var-data'>[357404 values with dtype=float64]</pre></li><li class='xr-var-item'><div class='xr-var-name'><span>rhop</span></div><div class='xr-var-dims'>(t, z, y, x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-78128fcc-e854-46f3-9026-b1dd2ac139f1' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-78128fcc-e854-46f3-9026-b1dd2ac139f1' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-f3834387-efd9-47c1-90eb-4ee2f7de640f' class='xr-var-data-in' type='checkbox'><label for='data-f3834387-efd9-47c1-90eb-4ee2f7de640f' 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><pre class='xr-var-data'>[14296160 values with dtype=float64]</pre></li></ul></div></li><li class='xr-section-item'><input id='section-765717ac-46b4-4b05-8168-33e0c93fdd33' class='xr-section-summary-in' type='checkbox'  checked><label for='section-765717ac-46b4-4b05-8168-33e0c93fdd33' class='xr-section-summary' >Attributes: <span>(2)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><dl class='xr-attrs'><dt><span>file_name :</span></dt><dd>SalishSea_04914000_restart.nc</dd><dt><span>TimeStamp :</span></dt><dd>13/11/2020 10:56:10 -0800</dd></dl></div></li></ul></div></div>"
      ],
      "text/plain": [
       "<xarray.Dataset>\n",
       "Dimensions:       (t: 1, x: 398, y: 898, z: 40)\n",
       "Dimensions without coordinates: t, x, y, z\n",
       "Data variables:\n",
       "    nav_lon       (y, x) float32 ...\n",
       "    nav_lat       (y, x) float32 ...\n",
       "    nav_lev       (z) float32 ...\n",
       "    time_counter  (t) float64 ...\n",
       "    kt            float64 ...\n",
       "    ndastp        float64 ...\n",
       "    adatrj        float64 ...\n",
       "    ssh_ibb       (t, y, x) float64 ...\n",
       "    rnf_b         (t, y, x) float64 ...\n",
       "    rnf_hc_b      (t, y, x) float64 ...\n",
       "    rnf_sc_b      (t, y, x) float64 ...\n",
       "    utau_b        (t, y, x) float64 ...\n",
       "    vtau_b        (t, y, x) float64 ...\n",
       "    qns_b         (t, y, x) float64 ...\n",
       "    emp_b         (t, y, x) float64 ...\n",
       "    sfx_b         (t, y, x) float64 ...\n",
       "    en            (t, z, y, x) float64 ...\n",
       "    avt           (t, z, y, x) float64 ...\n",
       "    avm           (t, z, y, x) float64 ...\n",
       "    avmu          (t, z, y, x) float64 ...\n",
       "    avmv          (t, z, y, x) float64 ...\n",
       "    mxln          (t, z, y, x) float64 ...\n",
       "    ub2_b         (t, y, x) float64 ...\n",
       "    vb2_b         (t, y, x) float64 ...\n",
       "    sbc_hc_b      (t, y, x) float64 ...\n",
       "    sbc_sc_b      (t, y, x) float64 ...\n",
       "    qsr_hc_b      (t, z, y, x) float64 ...\n",
       "    fraqsr_1lev   (t, y, x) float64 ...\n",
       "    fse3t_b       (t, z, y, x) float64 ...\n",
       "    fse3t_n       (t, z, y, x) float64 ...\n",
       "    rdt           float64 ...\n",
       "    rdttra1       float64 ...\n",
       "    ub            (t, z, y, x) float64 ...\n",
       "    vb            (t, z, y, x) float64 ...\n",
       "    tb            (t, z, y, x) float64 ...\n",
       "    sb            (t, z, y, x) float64 ...\n",
       "    rotb          (t, z, y, x) float64 ...\n",
       "    hdivb         (t, z, y, x) float64 ...\n",
       "    sshb          (t, y, x) float64 ...\n",
       "    un            (t, z, y, x) float64 ...\n",
       "    vn            (t, z, y, x) float64 ...\n",
       "    tn            (t, z, y, x) float64 ...\n",
       "    sn            (t, z, y, x) float64 ...\n",
       "    rotn          (t, z, y, x) float64 ...\n",
       "    hdivn         (t, z, y, x) float64 ...\n",
       "    sshn          (t, y, x) float64 ...\n",
       "    rhop          (t, z, y, x) float64 ...\n",
       "Attributes:\n",
       "    file_name:  SalishSea_04914000_restart.nc\n",
       "    TimeStamp:  13/11/2020 10:56:10 -0800"
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "restart_file = xr.open_dataset('/results/SalishSea/nowcast-agrif.201702/13nov20/SalishSea_04914000_restart.nc')\n",
    "restart_file"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [],
   "source": [
    "prob = 465-1\n",
    "pz = 14-1"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Look at velocity"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYoAAAEGCAYAAAB7DNKzAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nOzdeXxV1bn4/8/aZ848DxCSACqDEMKoCKLggCPWobbK76LiFelVr7a2Xr29vcXOA35Vai3S6i1qW62K1glFREEUZJ7nIRMJZJ7PvNfvjx1CkBARkpwAz/v1Oq+cYe2z1wlkP2dNz1Jaa4QQQojjMSJdASGEED2bBAohhBAdkkAhhBCiQxIohBBCdEgChRBCiA7ZI12BrpCSkqJzc3MjXQ0hhDhtrF27tlJrndrea2dkoMjNzWXNmjWRroYQQpw2lFKFx3tNup6EEEJ0SAKFEEKIDkmgEEII0aEzcoyiPcFgkJKSEnw+X6Srctpwu91kZWXhcDgiXRUhRASdNYGipKSE2NhYcnNzUUpFujo9ntaaqqoqSkpK6Nu3b6SrI4SIoIh1PSml+iilPlFKbVdKbVVKPdhOGaWUmqOU2qOU2qSUGnGy5/P5fCQnJ0uQOEFKKZKTk6UFJoSIaIsiBDystV6nlIoF1iqlPtJab2tT5mrg3JbbBcCfWn6eFAkS34z8voQQEMEWhda6TGu9ruV+A7Ad6P2VYjcAL2rLSiBBKZXZzVUVQoge4aGHHuKhhx7q9vP2iFlPSqlcYDjw5Vde6g0Ut3lcwrHB5PB7zFBKrVFKramoqOiKap4ym81Gfn4+w4YNY8SIEXzxxRdfe8ycOXMYNGgQU6dO7dS6VFVVMXHiRGJiYrj//vs79b2FEF1jw4YNbNiwodvPG/HBbKVUDPAG8JDWuv6rL7dzSLs7LWmt5wHzAEaNGtUjd2PyeDyt/8gffvghjz32GEuXLu3wmGeffZaFCxee8IByKBTCbv/6f1a3283Pf/5ztmzZwpYtW07ovYUQZ6eIBgqllAMrSPxNa72gnSIlQJ82j7OA0lM5Z+6j753K4Sek4DfXfm2Z+vp6EhMTWx///ve/55///Cd+v58bb7yRxx9/nJkzZ7Jv3z6mTJnC9OnTueOOO5g+fTr79u0jKiqKefPmkZeXx6xZsygtLaWgoICUlBSefvppZs6cSVFREQBPPfUU48aNO+r80dHRjB8/nj179nTuhxdCnHEiFiiUNVL6PLBda/3/jlPsbeB+pdQrWIPYdVrrsu6qY2fzer3k5+fj8/koKytjyZIlACxatIjdu3ezatUqtNZMmTKFZcuWMXfuXD744AM++eQTUlJSeOCBBxg+fDhvvfUWS5YsYdq0aa0tlLVr17J8+XI8Hg+333473//+9xk/fjxFRUVMnjyZ7du3R/KjCyFOY5FsUYwD/g3YrJQ63On230A2gNZ6LvA+cA2wB2gG7opAPTtN266nFStWMG3aNLZs2cKiRYtYtGgRw4cPB6CxsZHdu3czYcKEo45fvnw5b7zxBgCTJk2iqqqKuro6AKZMmYLH4wFg8eLFbNt2ZPJYfX09DQ0NxMbGdvlnFEKceSIWKLTWy2l/DKJtGQ3c1z016l5jx46lsrKSiooKtNY89thj3HvvvR0eY/06jnZ4Cmt0dHTrc6ZpsmLFitbAIYQQp6JHzHo6G+3YsYNwOExycjKTJ0/mhRdeoLGxEYADBw5QXl5+zDETJkzgb3/7GwCffvopKSkpxMXFHVPuyiuv5Jlnnml9HIlZEkKIM0fEZz11txMZaO4qh8cowGodzJ8/H5vNxpVXXsn27dsZO3YsADExMbz88sukpaUddfysWbO46667yMvLIyoqivnz57d7njlz5nDfffeRl5dHKBRiwoQJzJ0795hyubm51NfXEwgEeOutt1i0aBGDBw/u5E8thDjdqfa6M053o0aN0l/duGj79u0MGjQoQjU6fcnvTYie49JLLwWsHoXOppRaq7Ue1d5r0vUkhBCiQxIohBBCdEgChRBCiA5JoBBCCNEhCRRCCCE6JIFCCCFEhyRQdCOlFA8//HDr49mzZzNr1qyTeq+//vWvlJa2nx/xtdde4/zzz8cwDL46TVgIIb4pCRTdyOVysWDBAiorK0/5vToKFEOGDGHBggXH5IoSQoiTIYGiG9ntdmbMmMGTTz55zGsVFRXcfPPNjB49mtGjR/P5558DcMMNN/Diiy8C8NxzzzF16lRef/111qxZw9SpU8nPz8fr9R71XoMGDWLAgAFd/4GEEGeFsy6FBwALH4WDmzv3PTOGwtW/+dpih1NrPPLII0c9/+CDD7abGnzevHmMGzeOvn378sQTT7By5UqSkpJ45plnmD17NqNGtbuQUgghOs3ZGSgiKC4ujmnTpjFnzpyjsrseLzV4eno6P/vZz5g4cSJvvvkmSUlJkai2EOIsdnYGihP45t+VHnroIUaMGMFddx3ZXqOj1OCbN28mOTn5uGMSQgjRlWSMIgKSkpK49dZbef7551ufO15q8FWrVrFw4ULWr1/P7Nmz2b9/PwCxsbE0NDR0b8WFEGclCRQR8vDDDx81+2nOnDmsWbOGvLw8Bg8ezNy5c/H7/dxzzz288MIL9OrViyeeeILp06ejtebOO+9k5syZ7Q5mv/nmm2RlZbFixQquvfZaJk+e3N0fTwhxBpE046JD8nsToueQNONCCCF6JAkUQgghOiSBQgghRIckUAghhOiQBAohhBAdkkAhhBCiQxIoupHNZiM/P59hw4YxYsQIvvjii689Zs6cOQwaNIipU6d2al1WrVpFfn5+a33efPPNTn1/IcSZ4+xM4REhHo+ndcX1hx9+yGOPPcbSpUs7PObZZ59l4cKF9O3b94TOEQqFsNu//p91yJAhrFmzBrvdTllZGcOGDeP6668/oWOFEGeXs+6qMHT+0C4/x+Y7vj4zbX19PYmJia2Pf//73/PPf/4Tv9/PjTfeyOOPP87MmTPZt28fU6ZMYfr06dxxxx1Mnz6dffv2ERUVxbx588jLy2PWrFmUlpZSUFBASkoKTz/9NDNnzqSoqAiAp556inHjxh11/qioqNb7Pp8PpVQnfXohxJkmooFCKfUCcB1QrrUe0s7rlwL/Ava3PLVAa/2z7qth5/J6veTn5+Pz+SgrK2PJkiUALFq0iN27d7Nq1Sq01kyZMoVly5Yxd+5cPvjgAz755BNSUlJ44IEHGD58OG+99RZLlixh2rRprS2UtWvXsnz5cjweD7fffnu7Kcu/6ssvv2T69OkUFhby0ksvSWtCCNGuSF8Z/go8A7zYQZnPtNbXdU91ulbbrqcVK1Ywbdo0tmzZwqJFi1i0aBHDhw8HoLGxkd27dx+zQ93y5ct54403AJg0aRJVVVXU1dUBMGXKlNbMs8dLWR4bG3vU+11wwQVs3bqV7du3c8cdd3D11Vfjdru75sMLIU5bEQ0UWutlSqncSNYhUsaOHUtlZSUVFRVorXnssce49957Ozymvbxch7uMoqOjW5/rKGV5ewYNGkR0dDRbtmyRjZCEEMc4HWY9jVVKbVRKLVRKnX+8QkqpGUqpNUqpNRUVFd1Zv5OyY8cOwuEwycnJTJ48mRdeeIHGxkYADhw4QHl5+THHTJgwgb/97W+AlRQsJSWFuLi4Y8odL2V5W/v37ycUCgFQWFjIzp07yc3N7YyPJoQ4w0S66+nrrANytNaNSqlrgLeAc9srqLWeB8wDK3vs8d7wRAaau8rhMQqwWgfz58/HZrNx5ZVXsn37dsaOHQtATEwML7/8MmlpaUcdP2vWLO666y7y8vKIiopi/vz57Z5nzpw5rVuuhkIhJkyYwNy5c48qs3z5cn7zm9/gcDgwDINnn32WlJSULvjUQojTXcTTjLd0Pb3b3mB2O2ULgFFa68qOykma8c4jvzcheg5JM94OpVSGaumEV0qNwapvVWRrJYQQZ5dIT4/9B3ApkKKUKgF+CjgAtNZzgVuA7ymlQoAX+K6OdBNICCHOMpGe9XTb17z+DNb0WSGEEBHSo7uehBBCRJ4ECiGEEB2SQCGEEKJDEii6kVKKhx9+uPXx7NmzmTVr1km911//+ldKS0vbfe0nP/kJeXl55Ofnc+WVVx63nBBCnAgJFN3I5XKxYMECKis7XAZyQjoKFD/60Y/YtGkTGzZs4LrrruNnPztt8ygKIXoACRTdyG63M2PGDJ588sljXquoqODmm29m9OjRjB49ms8//xyAG264gRdftHImPvfcc0ydOpXXX3+dNWvWMHXqVPLz8/F6vUe9V9u0Hk1NTZJCXAhxSnp6Co8ucfBXv8K/fUenvqdr0EAy/vu/v7bc4dQajzzyyFHPP/jgg+2mBp83bx7jxo2jb9++PPHEE6xcuZKkpCSeeeYZZs+efdwkfj/+8Y958cUXiY+P55NPPumUzyiEODtJi6KbxcXFMW3aNObMmXPU84sXL+b+++8nPz+fKVOmtKYGT09P52c/+xkTJ07kiSeeICkp6YTO88tf/pLi4mKmTp16VIJAIYT4ps7KFsWJfPPvSg899BAjRozgrrvuan2uo9TgmzdvJjk5+aQGpW+//XauvfZaHn/88VOqsxDi7CUtighISkri1ltv5fnnn2997nipwVetWsXChQtZv349s2fPZv9+a7O/2NhYGhoa2n3/3bt3t95/++23GThwYFd8DCHEWUICRYQ8/PDDR81+mjNnDmvWrCEvL4/Bgwczd+5c/H4/99xzDy+88AK9evXiiSeeYPr06WitufPOO5k5c2a7g9mPPvooQ4YMIS8vj0WLFvH0009398cTQpxBIp5mvCtImvHOI783IXoOSTMuhBCiR5JAIYQQokNnVaA4E7vZupL8voQQcBYFCrfbTVVVlVz8TpDWmqqqKtxud6SrIoSIsLNmHUVWVhYlJSVUVFREuiqnDbfbTVZWVqSrIYSIsLMmUDgcDvr27RvpagghxGnnrOl6EkIIcXIkUAghhOiQBAohhBAdkkAhhBCiQxIohBBCdEgChRBCiA5JoBBCCNEhCRRCCCE6JIFCCCFEhyIaKJRSLyilypVSW47zulJKzVFK7VFKbVJKjejuOgohxNku0i2KvwJXdfD61cC5LbcZwJ+6oU5CCCHaiGig0FovA6o7KHID8KK2rAQSlFKZ3VM7IYQQEPkWxdfpDRS3eVzS8twxlFIzlFJrlFJrJEOsEEJ0np4eKFQ7z7W7oYTWep7WepTWelRqamoXV0sIIc4ePT1QlAB92jzOAkojVBchhDgr9fRA8TYwrWX204VAnda6LNKVEkKIs0lENy5SSv0DuBRIUUqVAD8FHABa67nA+8A1wB6gGbgrMjUVQoizV0QDhdb6tq95XQP3dVN1hBBCtKOndz0JIYSIMAkUQgghOiSBQgghRIckUAghhOiQBAohhBAdkkAhhBCngbCp8QbDhM12k1N0KQkUQgjRwxVXN3PZE5+ysbiWtYU1fLDlYLeeXwKFEEL0cD99eyvFVQ1kqQrsOsj/vLWZUNjstvNLoBBCiB6svN7Hkh3lXGxsIktVkqgaqGwMsL+yqdvqIIFCCCF6sPc3W+ntrretBMBOGIA6b7Db6iCBQggherB3N5XhIsCVxhrgSKCo90mgEEKIs15prZc1hTVcYmwkVnkBsGONTXRniyKiSQGFEEIc33ubrG6n62wrqdKxeAliO9z11CwtCiGEOOu9s6kUDz4uN9axMDyGEPY2YxShbquHBAohhOiBCqua2FRSxyRjA1HKz7vmWELawK6sricZoxBCiLPcu63dTis4pBNYZQ4khO1I15PMehJCiLPbOxtLiaGZScYG3g9fgIlBGJtMjxVCCAF7yhvYcbCBy411uFSQd8JjAQhhYMPEwOx5s56UUuOAWUBOyzEKa6fSfl1XNSGEODs9v7wAgCm2Lzigk1mvzwEghA2AOJqo98Z3W31OdHrs88D3gbXQ0u4RQgjR6dYUVPOPVUVcaGxjkm0DT4VuQrd0/oRbAkW8aqK+p7UogDqt9cIurYkQQpzlAiGTxxZsxkWAX9qfp9BM40+hKa2vH2lRNLOvB45RfKKU+r1SaqxSasThW5fWTAghzjLzlu1ld3kj/2F/m/5GGT8O3Y0fJz+5bjBKQUhbl+x41URTIEywmzLInmiL4oKWnyNbfipAA5M6vUZCCHEW2l/ZxJwle+ivDvA92794MzyO5eZQhmXFc+dFuTy9eBd1h7uesDLH1nuDJMe4urxuJxooPm3nue7fZkkIIc5AWmv+563NBEMhfuV8niY8/CL4/2EzFL+6aSg2QxEf5aCgzRgFWFNkuyNQnGjXU2ObWwi4CsjtojoJIcRZZVtZPZ/vqeJm22dcYOzgV6HbqSKeu8f35fxe1uymeI+DUMslO+5wi8LXPWk8TqhFobV+ou1jpdRs4O0uqZEQQpxlNhTXAnCpsZESncJr4UvoneDhocvPbS0T73FgYqBRR7UousPJLriLAmQNhRBCdIKapgAACTRwUCcBiuvyMolyHvkuH+9xANaiu8NjFD0qUCilNiulNrXctgI7gadP9eRKqauUUjuVUnuUUo+28/qlSqk6pdSGltv/nuo5hRCip6lusi74SaqRGh1r3Y92HlUmzn04UNiI6+YWxYkOZl/X5n4IOKS1PqXOMaWUDfgjcAVQAqxWSr2ttd72laKfaa2vO+YNhBDiDFHb3NKiUA1sMXMBSIw6OlAcblGEsRFHM0C3Lbo70TGKwi449xhgj9Z6H4BS6hXgBuCrgUIIIc5o1S2BIokGqrFaFInRTuoD9XjsHhyGg7jDXU/a1jpG0aMCRRfpDRS3eVzCkfUabY1VSm0ESoEfaq23tvdmSqkZwAyA7OzsTq6qEEJ0nZrmIG78uFWQWh0DmLxe9Dt+uGYxCsX43uNJNi8GIjNGEclAodp57qtrM9YBOVrrRqXUNcBbwLnHHgZa63nAPIBRo0bJGg8hxGmjpilAEg0AVBOLI+lzvji0iNE7TZwh+CL4KUH7UmyeIqrDiugeOkbRFUqAPm0eZ2G1Glpprevb3H9fKfWsUipFa13ZTXUUQoguV9MUIFs1Wvd1NM6kJcQ1aR76l4kjDPUeWDpU8ZtwiDo7/CrNDcW6Z8166iKrgXOVUn2VUk7gu3xlbYZSKkMppVruj8Gqb1W311QIIbpIMGzS4A+RqKwWRXF0LYajlkkbNY4w/Okag605iqvXaM4p1UT74eNoDzG2ijO/RaG1Diml7gc+BGzAC1rrrUqpmS2vzwVuAb6nlAoBXuC7WmvpVhJCnDFq2gxkA1Qk7EGZmivWm2zOUXwyzOCTYZDYoOERiPKBVopo5wHqfTndUsdIdj2htX4feP8rz81tc/8Z4JnurpcQQnSX2marVZCgGii02wnEFDFqlya1HuZfbg3lJroSqaEGhw0cLTsCuZ0V1DSf+V1PQghx1qtuWZWdpBp4NS4GgMlrNZWxsOZcRXZsNlf3vRqAoA0cLSvY7I4aGvwhTLPrO1kkUAghRAQdTt8Rrep5MyaGzCrNsALN4uEGpqG4dcCt9I7pDUDQrnCErcCgHLVoDQ3dkBhQAoUQQkTQ4e6j4thyGm0GV64zCRnwcb7CbXPzrXO+RWZMJgBBe0uLQmvCDmtMozsGtCVQCCFEBFmD2Zo18fW4AppLN2tWDlTURSuu7ns18a54MqNbAoUNDA2xXgjYW9J4+CRQCCHEGa2mKYDhPkCx22T8Vmv664cjrUvzdwd+F+BIoGiZfpRSD16HD+ietRQSKIQQIoKqmwPYPAUAjNirKUuEnb0hLzWPwcmDAUhyJ+GyuQhaG9yRUqcJGiYYXgkUQghxpqttDmJzlwGQU67Zl6FAKSb1mdRaRilFZnTmUS0KAMNRK4FCCCHOdNVNAWyuUjw+TVodFKZZaycGJA04qlxGdAYhA0wFKfVHZj5JoBBCiDNcdbMXw1VOdksGu6JU6+d5iecdVa7tOEXbFkV3pBqXQCGEEBFUEygFI0x2udVKKExTxDsTSPWkHlWudYqsTZFSZ5WVrichhDjDhcImXlUCWOMTTS6oioMBSefRkg+1VdsWRWpLi0K6noQQ4gxX6w1iuKyB7OxyTWEaoNQx3U7AUWspEprAEdIYdgkUQghxRqtpCmBzHURpTXYFFKVarYj2AkWv6F7AkbUUSfVgOGpkjEIIIc5k1U0BDHcZKXUQFTj+jCeA9Oh0oO2iO41yNFDn83Z5PSVQCCFEhByor8Jw1JHTMpBdlKZAK/on9D+mrNPmxGE4WhfdHR6nqAt0/V5uEiiEECJCdtXsBiCn3HpclArRRiYum6vd8s6W1dkmkFJnPdcYrqCr93OTQCGEEBFS2LAHgOwKzcEE8DsVyY7j71rnsjnRCmpjjiy60/YamgLhLq2nBAohhIiQ0ua9wOEZT9b4REZUv+OWd9qcAFTGdW8aDwkUQggRIZWBQpxBTWYNFKVZz+XGnHPc8q2BIv7IojvVDTOfJFAIIUQEhM0wjWYJWZXWHhOHWxTntjM19rDDYxetLQqtMex10qIQQogzUXFDMaYKtM54KkxV6LCbvgm9j3uM0zjc9aRwhiGuuXtWZ0ugEEKICNhZsxOwxid8DihPhLAvk6Ro53GPOdz1VBFvPU6pt8YoapsDXVpXCRRCCBEBu2p2AZBTYU2L1Uph+jNJ7CBQ2A07CoPKOKubKrVOo4wgh5q6di2FBAohhIiAndW7QGuyy7W10A4w/RkkeBwdHmdTDirjrPuHZz6VNZV1ZVUlUAghRCTsqN5JYiPEea3xCQC3mYXd1vFl2W44aXKDz3FkLUVF88EurasECiGE6GZ1/joONZe17kFhpe6AREefrz3WYThBKSrij7Qoqv2HurK6kQ0USqmrlFI7lVJ7lFKPtvO6UkrNaXl9k1JqRCTqKYQQnWntobUAnHdAYyooSAMzkExSVOzXHus0rK6pyrgjaynqQxVdV1nA3qXv3gGllA34I3AFUAKsVkq9rbXe1qbY1cC5LbcLgD+1/OxUE+bfQbW56RsfZzMTmJh5E09d80BnV0kIcQZbWbYSgKEFmj2Z4HUrQjX9SUo4/kD2Yc42ayn6tfQ4HWQxQ/5veGuZTPsFfPRvczutvhELFMAYYI/Weh+AUuoV4AagbaC4AXhRWxmvViqlEpRSmVrrTh25MXUQZYToX6pR3yC31qHECj6umMefV5/LPaOv6swqCSHOYCvLVuLxa84thbfGWuMT4aZzSOjV8UA2gNvuAkJUxinimzUDijVhQwMmYRvsz1CEdahT6xvJQNEbKG7zuIRjWwvtlekNdMkQ/0//Hsb9Ddat7OoF/3OHnff2LpZAIYQ4IYeaDrG/bj8jizQ2DVtyWwJFc3+Sor6+ReGxu/ARoizJevzzl48kBKyNhhn/2fmX9UgGCtXOc1/9Pn8iZayCSs0AZgBkZ2efVIVm32RgnGCL4vINmiGFGrSm2t+1U9OEEGeOLw9+CcCQAo3fDrt6K7J8NraHoztcQ3FYtNNNKBTLqgFlPH6bgaNN4tiQrWvqHMlAUQK0HeLPAkpPogwAWut5wDyAUaNGnVRy9k39TnxsP6vSZMwuTbQPmlxdO5AkhDhzrCw9Mj6xvY8iaFfkNESxHUg8gRYFwBOXzuZHSx9jS04J6pv0l5+kSAaK1cC5Sqm+wAHgu8DtXynzNnB/y/jFBUBdZ49PACy87XnCpnlCZTeVlXLfZzdRnmA9Tq2DgrQqtNYo1V4DSAghLFprVpatJKFRk10JS4da14z0ZuuCkhj19WMUAJP6D2Nt//dp8HvbvXY5jM69tEcsUGitQ0qp+4EPARvwgtZ6q1JqZsvrc4H3gWuAPUAzcFdX1CXW5TnhssMyc9Cmg4p4K7dKap2mICNAta+aZE9yV1RPCHGG2Fe3jwpvBeMLrFbAllyFQ2s8zVaO8RPpemrrm1y7TkUkWxRord/HCgZtn5vb5r4G7uvuenUkzu2AYBLl8da8tLSW7Qh3VxWTnHVmBopQ2OTzvVVUN/k5Ny2WgRmxX7t6VAhxrLbTYhvcUJAOo3x+GrXVougoIWAkRTRQnI6UUjhIpsl9kGan1aIA2FlVyIVZ+RGuXecKm5q3Nx5gzsd72F/Z1Pq822EwtHc8o3OT+O7obLKToyJYSyFOHytLV4LWDC3QbMlVaKW40Otjo7YW2iWcYNdTd5NAcRJibKnUtiyhT6u1nttbUxTZSnWiqkY/H+8oZ+7SveyrsAJEvtrDYKOQrWYO24M5rC6oYXVBDX/+bB+3j8nmgcvOJSWm/Q3hhRAQMkOsPrSazGpIaYA3WqbFXujz8SkxwIkPZnc3CRQnIdGZQa0JFfGqtUVR0nAgwrU6NXvKG3h7YxlLd5az6UAdumUiRZ7ay/ftrzPRtrG1rF/b2a5zWBQexf+FJzN/RSGvrS3hzotyuWxQOnlZ8Tika0qIo2yp3EJTsIlxLeMTm3MVtrCdwf4ANTqWWJe9x/7dSKA4CWmeTPY3QXkCDC4CtOaQ9/RcS7GvopGnFu/mnU2lrcEhjkYuNrZwo+0zLretp0bH8Jvgd1lojmGQKiLf2MtIYyePOF7lTvuHzAndyCuBiTz76V6e/XQvHoeNkTmJ5KZEkRLjIjnGRazLTiBk4guF8QdNnHaDjHg3mfFuMuM9JEc7MQyZNSbOXG3HJ8rjoTxRcU6jG7TBQZ1IWlzPbZFLoDgJWbG9+LLJalFEBay1FDW2rs3e2Nn2lDfw3NJ9vLGuBFNDLM3cbvuYy21rGa72YFcm1TqG3wW/zfzoTHSvL3E4t/B5IJFPfSmY/on0DYznkdByfuH4P+62vc9b4fF8qQexPngOy/eEWb7nxOvjtBlkJrjpFe+hV4KH/mnRDMqIY1BmHOlxLpl6LE5rWms+Kf4EZVoLdVcOtP4/3+Qr4dXwRGqI49uD0iNcy+OTQHES+iVkw0Eob9mOMK0OCj0VPX4tRZ03yDsbS3l9bQkbiq3BFQ8+7rQt4l77OySoJjabuTwbnsKn4Xw2ux04+ryPzWOtJDUB3PU43IWAtfLxIcAVGMxIfzOTfB/ye9+/SAvCFt2XMp1MlY6lSsfThBsXQVwqgJsAPlyU6SQO6iTKdBKl4RQKq0wKqwpfiMUAACAASURBVJqPqXdClINzUmM4Jy2G/qkx9E+Lpl9KDFmJHpl9JU4LS4qWsK1qGxfs0kT7YX1/6zox3BvirtDNRDttzJjQL8K1PD4JFCchJyEVbTqpiPcD1syn/T18LcU/1xQz6+2tNAes9f6p1HKzbRl3298nVdWzODycJ0PfZpuRii16N/Epa3G7rHGJxAbNd5aZDCjRHEpUlCZBWZKiME2xPwP8zka+cMIXsVbymbigg3xviPN9RQwMNTE51EhaONya096rnbgIYnxlRWm99lCqUzigU9its9hmZrFTZ7HPm8GaQh9rCquwsrpYf2QOmyI7KYqBmXFc2DeJsf2T6Z8a06ODtTj7hMwQT617CsPU3LbUpCgFVp+rGOgPsNg3mQoSePiS/j16MogEipOQEe/BDCRScXgtRcvMpwONB3pkoPjXhgM88vom7ISYbKzj27alXGpsxK5MPgsP4d/D32JbfDWexPeJcRUDmgDgCGquX6X51goTexg29lUkN2jOLwRXyLrIhwwoTIPdvRTbshVbcxT1UUGWOWBZnB2It25aoU0HWjtA20ErbMrERhiDMEqFQYVBaUyqMFUlYWUFKnfL7TDDNFCmk7Dppsx0U1oTy+LyBPTyBGLsaYxIP5+rBw7jkvPSe/Qfnzg7vLnnTQrqC7hsk6ZXNfz2FgNtKKZVB3g0dB1psS7uvrhvpKvZIQkUJyEt1oUOJtIUc5Bm15G1FKWNpeSl5kW4dkfbVFLLI69vwsDkOceTXGZbzyGdwJ/D1/JP80IOJBbiSn4Dt60JDdjCmvMOQP4+k4u3alLq4csBipcnGhxKtL6pK61Jroe+hzTnlFqpki/ZorlqnfV72JcOW3MUB5IVhxLhYKKiOlaDLYAi0Fo3DbSXDFlpjS0MLlNjM8EwQSswlfUzaA8TtvsAH1YOtCPbQAaBL4Owcn0U4S9ySHUM4pI+47g1bwxDesdLa0O08gXDVDcFqGkOUNscpM4bpN4bpN4XpN4bIsZtZ1z/FIb0jkMpRcgM4Q15aQ42EwgHiHfHE+uI7fD/VHOwmT9t+BOOoObbn5ns7A1rz1GM8vr4suF6mvDwkyvOI8rZsy/FPbt2PVRClAMVTgKlKI8/sjp7f11xxwd2s/J6HzNeXIs/ZPKY/R9cZlvPz4NTma8uxEhciyfpZVyGtU7i3BLN9atMhu3XeAJWS2FHH8Uz1ym25Ricl3ge/zviQRSK/XX72Vu3l82Vm3mlZg8ajWFq+pdZMzqGFphctVbjCB/dtRSwQdBu3UxlBQBDWzebCbYw2Fue+zpBm7VnsNdlpVauiFdUxkN5vGJPL0VBWhMqdjs1bOetigUsWJiAKzCUMekX8x8XXkneGbqK/kyjtaaouplV+6tZXVDN9rIGDAWpsa6WGXXO1imlqqVLUqNbZ/DplvcwtSZkag7W+SisaqaoupnqpkB7Z0Q5qrG5SzBc5Ty1sQJ3VBXKWUVI+44p7bF7yIjOICMqgyEpQxiTOYb81Hzcdjdaa17a9hIV3gquX6tJaoSnb7CBUny3yuC+8GWckxbDLSOzuui313kkUJwEpRRRRip+rAtUem3LN+manhMofMEwM15ay8F6HzcZy5hhf49Zzot4JaUST8zsli4eGFSkuflzk7wCTb0HPjtfsbGfYkuOwutSJLmT+N/h93PTOTdhM6zv7xdnXdx6nsZAI5srN7OufB0r0lfwr6wtLBgXRplWaySjRpNea41zOMNgD4EzZI0yHG4lmArCBoRtVoAKG6r1vmlYr8ORoGIPgSdgBTSPHxIboX+ZZswuWoNTs9NK37w9W7Guv6IwrYag8zM+b/qMzz78A+nqImaMuJXvDBstrYwuFDY1Ow82UOcNEu9xEB/lINZtp7LBT1F1M8XVzRyo9VHTFKC6OUBNU4AGXwhvMIw/FKbZH6bBf6TdqTDRbcapTo4mgUaGqgqSVB2Gu4Tm6BLqPDVUuZvw2cJHlTZNTUoNxHjBFQJXQOMIQ120ojKumYLgPvbX7WdF2Qr+vPnPOAwHmdG9KG8uxxf2EuXT3LjCZF1/6//jFU3NvN40jRB2Hr1q4GkxIUMCxUlKcKRzCGvm05BCQOses+iuqtHPQ69uYENxLcPVbh50v8B3UvuyPaoEB+AIaS7coblincnAA9Y38hcnGXw0XOF3KnLjcrm593jG9R7H6IzRuGzH7+ePccYwttdYxvYay33599EQaGDVwVVsq9pGcX0xRQ1FfFlfREOw4aQ+i03ZUEq1fCs00e1vRwJYXVYpddY+xAOLNYNKNLct1dy2FCriYN05ipUDFNty6qlQH/LLjR/yu7XZ3JM3nXtH3Yihev4f7OlAa822snreXHeAtzeWUt7g/8bvYSeEhwAe/Aw1DjDG2MEYtYPhhjXnuoo4KnU8VTqOEDY0ChOF5vB0B+v/iQKMlldsmKSrGtJVORuiYUlUFJ9Huam2HdnEIbZZM6jUah1nV2h6VWsyqsEZPraOh4UMOJQA27IVW3IVW3ICFIULiPZB3zq4cp1JjA/+cYmBoTXFh25ltTmKywamcdmgtG/8u4kECRQnKT0qk0PaalF4ApoYLxxqjvyiu1X7q3ngH+s4VO+njzrIrYl/4Dsp6TTZwiTXa65ZbXLpJk2sD0oT4YUrDD4epgg7bVzb91ruHno3/RP6n/T5Y52xXJZ9GZdlX3bU88FwEH/Y33oztYmhjNaLs92wY1d27IYdm2FrffzVb/taa/xhP43BRpqDzTQEGzjYdJCyxjKK6ktYf3AHXyZs4fPzrYtTfKNmxF7NqN2aSzZrJq/TVMbBsiGKpUMMypKLeHbbLF7b+388cuH9XJlzpQSMk1BW52XV/mpW7a9mxb6q1tQvAOOMzQxXe3CqIE5C1jTpNlOlPQSIV00k0EiiaiAWLw519JU5rBXbdA6vhCcSxE6yqiOVOtJVDbaWQKDQGC1fJXTLMwAmBiFgpxtejXGwNiYdv81Kze3xaS7YbTJ8r7W+4XA3sgkcTITSZMWGflCapKiLBr8D/HZFyAYJTVarOaVek1UJ47ZrrthgBSivEzxteraW5CkK0xW2muGs9o/nurxMfntz3mnTmpVAcZKyYnqzqQEqWvalSKuDwuhDEVtLYZqaPy3dy//7aBdhUzPJ8RlxGW/wq5go0JpLNpnc9ZGJMwSrBig+yrdmKNkNB9f1v457ht5DdtzJ7Qx4Ihw2Bw6bg5iWnDYnSymF2+7GbXdDS4bl85PPP6pM0Ayyo2oHi/d/wUcFS/kkeiufDDNxBI8EjG+t0Nz0RZiNuYrXLzbYmVXIj5b+iLnxc/npRT9leNrwds4u2qpo8PPm+hJeX1vCrkONx7weRxM/dbzIzbbPWp/zawd+7Phx4seBTzvx4qROR1NGDjVmDHVE49UufDjx4qJUJ7PWPI8GvmnySRObpwh73GbssVswHHVAgCif5rLNmnHbrJan3YQml5VS48ORij2Ziv3p4HN95e9Yg13bcJng0iYl6Zpgm+8UhqnpVwZDCzVxTZqKeEVFgjVuVpQGZiiaUO11/O7mPL49Kuu0CRIggeKk9Y5LRtc6KW+zlmJfZoAqXxUpnpRurUu9L8gPXt3I4u2HcBBgStKfWJVSQr3NTVyT5t6FJqN3a7b1gT9eZ6MiQZEelc79532bm869idSo1G6tb1dzGA6Gpg5laOpQvj/mXur8dXxUsJRXtrzHikFfsmJwmIRGK2Bcu8rk5y+F2ZSreH28wQ72csfCO5g2eBr3D7/fCkiiVZM/xLJdFby5/gBLdpQTMjUGJueqUrw4qdUxNOJhgrGJ3zr+TCq1/CH0Lf4YugEfTk5mbMFhU+TlJDCmbxJjcpOIaRnjqGj0U90UwLRGrDF1mOpQEeWBHZQHdnAwsBWvac1dt4c0w3ZrJmzRjNxtjZeVJMM7FyjW9zfY1RvMlhQybpubQcmDGJw8mBRXHw6Ux7C9yM2egwaVDW0HwDXK1oyy1+L0FBMVtYOC9EL29PYeKWHa0aF4wvVZZKnree6+a+ifempfliJBAsVJSo93YwbbrKVoabIeaDzQrYFi96EG7n1pLfsqm0hzbaNXr5dY7NagDcZvMbnjYxOPH+ZfZvD+aEVqVDpPjPkvJmVPwt7Ju2D1VPGueG4ZMIVbBkyh2lvDkytf5b39b/OvscUsHKW4Yr3mhpUmP3s5zGeDFf93hcH8bfNZdmAZvxj3ix435bm71XmDfLjlIIu2HWTZ7koCIavbJkcd5Bb7Mm62LaOXqm4tH9A2nCrMbrM39wa/z3bjHAb3iccbCLVMQQ0R57GTnRRFn8QospKiSIt1kRTtJDHKSbzHgcdpw+0wcNttxLiPJMvTWlPhrcDwlOJzF1HjKqK4roCC+gKK6osImEcu5I6gZtR+azxu1B5NlB/qPbB4uGLZEIN9GYBSuG1uxmWMZnzv8YxMH0n/hP7H/duobgqw42A9+yubqPeGqPcFafAFqWocyYFaLyWlzTQGDqIMP2lRqQxIy2Rg7zhG5SYxaWAattM0n9nZcaXoAmmxbnQwkebYQzS6IbX2yFqKYanDuqUOCzeX8fBrG2kOBOiX+k+qkjewVyl6V2ru/tBkSJFmdyY8e52NAymKybmT+cmFPyHeFd8t9euJkjyJ/HziTH48/t+597VXWdX8Cu+NKWDRcMUNKzU3fWEypDDMvKsN1p67n2kLp/HDUT9k6qCpp1VXQdjUbCutp9EfYkROAi677esPakNrzeqCGl5ZXcT7m8vwBa3g0JsKJtvWcJVtFWOMnYS14jMzjyfNWwBIoJFYo55DhoOtmVcwdkACM/o4cTu9bd7djqEM7EYIm2rCZvhwGA6chhOHzdqPoTbkxdvopSnURHFDMXtr97Kndg/76/bjDXm/Wl3QmuQGyDlkTWIYUKzpfxAcYWh0w8oBihUDrcHmsE3hsXu4Pvtyrup7FWMyxpxwyzEp2slF/VO4qP/xvwx6A2E0usevjfgmzpxP0s1SY12YQStlRUW8tXc2WC2K7vD+5jLu+/s6sFeQmzuPCk8DjiDctjzMtas0Pic8d5XBknxFjCuO31zwY67pe81pdbHrSm6Hnb/edjvPLBnNU59/QEzaO7x+8UFWn6e4790w//W6ydIhmr9M1vx29W/ZWrWV/x37v3js3bP15MnQWrOxpI63N5Ty3uZSDtW3dIvGuvjBFefx7ZFZ7U7FNE3N/qomtpXWs+tQA7sONbC1tJ6SGuuCnEotd9g+41rbSvKM/QDsMPvwq9AtLLAPoC66AcN1kLi4WgxnFV6zsmVIeS17i4FOnDUe16TJroGMak1GjSajBnpVazKrwR20yoQM2JsJ74+2pnlvybGCg13ZuTDzAq7tdy2XZV9GlKNrNtzyOL9ZUD4dSKA4SelxVtcTWDOfMqqPtCi62vqiGr7/6gbscauJzVhAlaHpe1DzwNthsqrg42GKv19q0BCluL7f9fxg1A+6fdzkdKCU4oHLzqN/Wizff7Uf9viPKUlbzGN32rj5c5ObvtD0PRTm9zfZeJd32Vu7lycnPknvmN6Rrvoxluw4xC/e294628hJkMuMTXgIsLBhDI8t2Mzzy/cz4+J+2AxFgy9IZWOAjSW1bCyupd4XAhVEGT6w+XAYzYyM3s54+wYG2fYSNDRvG+nMMi6gwO7B52jCcK1BqVWt6VV8gAprkuqtjXmifFZ3T7TP+mavtHUzdMv6GaNl/YwBoZZ1MyG7tUjOFQR3wLr4W7OKNH0qIK5NY8JU1pe00iTFtmzrZ3GqYm8GBB3WFyK7YWdcr4u4IucKJvaZeFa3pk+FBIqTlBztRIWOtCjy9gNaU9zFaymKq5u558U1mLGf4c54h5CpufELKz1AfRT84jsGm/oZ9Ivvx9MX/g+jM0Z3aX3OBNcMzaRXwnj+fb6bmv2Dyek1j39O8LG9j8lDb5n85q9h5kwxWH/Odr7z7nf41fhfMSFrQqSrDVgzj3727jbe2VgKaC4ytjLF+IKrbauIV1Ym3j1mL34X+g6LykfxyBubWo9V9jpsUfuITtxOctReAs6mo957V8sNDq9iDwPWFHCbtqaG5h7S5B6CnApNZrX1Dd/ZXl6WU9DshOJUWH2eojhFcSAZDiVaO0yGbUe3kOOccYxJHcqItBEMTxvOkJQhPboVeLqQQHGSDEMR70jFhzX9zR3UxHqhuL6ky85Z5w0y/a+rqTVW4sl4B2dQ88jr1qrqFQMV864y8EbZuGfI3Xxv2Pda+3vF18vvk8Bb913Ev893smP/TxiaMZfNfYt59C7FDxeE+a/XTF67GBZcVMt9H9/H3UPu5v7h90dsQoDWmtfXlvCL97ZT5w3SV5Xxa8dfuNDYToP2sMgcyTvhsbgI8Yj9FeY6n+R1xzl84MykzlNLubuJWod1RTcBM6TpWwYxviMr3j0B66LvDIIzZP3/Tq6H5AYr11dUy7ixCZQlWWsONva1vtlXxkOjW9HsgmYXBKw8kK03pa0ULraWm920VtwfTuHic4DPaa1bCNiBNl2mUfYosuOyGRrbh6zYLHLjcq1bfC6JrkTpXu0CEihOQaq7F8VYLQqwssgWRh9sXUzWmcKm5v6/r2Nf8yqisl7DFtZ8/02TIQWaP11j8EmeIiu2D3Mv/jX5afmdeu6zRVZiFK/NHMsD/1jPpzvvY6DvHSrSl/OTf7Nxzwcm3/nMZGCx4g9TDJ7f8jwbKjbw6/G/JjMms1vrWVrr5dEFm1m2qwI7If7D9i4P2t/Ej4P/Dt7NG+GL8ePEbijCRjWfxl1FfMLneB0+wNpLxB7SnHsAzi+0FpoNLNEdtgSCNmtQuCrOCgRbcqA4VVGQpihOBb/z6ItzqieVeFc86Y4YYpwxR63uP7zKPqRDhM0wIR0iZIYIhoMEzSAmJol2D1H2KDx2D0nuJPrF96N/Qn/6J/QnPSpdgkE3k0BxCjJiEykKuyhPsJKFpdZp9vYKUOWt6vS1CX9fVcQXB74kqs/fUJj8x7smI/dq5l1l8Mkwg6tyr2LWRbOIdkR36nnPNrFuB3+ZNorfLNzBX5ZfT6Y/G1fvv/PH6wy299FM/8jkdy+EefoGG2tZyzULrmNyn28zM/8echKSuvQCprXmldXF/PK97TT6Q5yv9jPb8RyDjCLeC49hVvAOKkhkRE4CowZWsMf/HmsOfYlG4zc1w/dYQeG8A5q+B4+kpShIg4+GK3ZkKWqjFV6XtbLY57S+zQfsoI8zrTM9Kp2LkgczKGkQ5ySeQ3ZsNtlx2dLdc4aRQHEK0uM8mBWJVMRZ/bZtZz51ZqCoaPDzu4+X4cl6EVSIuxaZXLxN8/dLDBYPN7g8+3J+ffGvz5p1EV3NbjP4n+sGc15GLP/zpoFjf28G9p7DkvwAe3opfvBmmJ/+Pcy7YxQLLvLzXtHfeGffmzgaJzEk/hLG9OnP0Kx4kqKcRDlteJw2kqKdJzVdMhg2WVtYw8fbD/Hx9nL2VTbhJMgP7G/yH7a3qSKOGYHvs8gcTXyUnZkTvGxuepZ/FFljEZlVmks3m1yy2cpeGrDD3gz4YJRq3UOkIcoKAoYySHAlEO2IJsURg8fuwWP34La7cdlcxDnjyIzJJCMqg8yYTLJjs3vk/iui88mV5RSkxbowSxPxxh6kwW1lSgUoqC/o1O6fWe+uxkx9AcPm5+bl1r4P74xRvDVWcWHmhfx2wm8lSHSBW0f1oX9qNPe+tJa1RbPol/IaxanrefROG3cuNrnuS83EjWEWjDP4cEQToYR32ci7rNuTSWj9YEx/JmYwAR1MBDOKUTlJXD0kk6uHZpAZf+w37lDYpNYbZG95I2sKa1hdUM3agpo22VM1I9Uuful4gYFGMW+EL+bx4L/RYIORQ7eiYtfyt8JtgJU2/pbPTYbv04SVtfXmX4Yp1vdXrQPAhjIYkDiAMRljGJM5hpHpI6VFKtolV5dTkBbnwvSnQex2ilOhT4UVKHbV7Oq0cyzffYgl1U9hj6lk9C6rn3zZ+YqXJhnkpQ7j6YlP47Q5O+184mgjc5J4/8GL+cPHe3h19W2Em0dh6/Uyc6/1snCUZuon1ur3q9dYiQbXnmOwL7MUm/srCSJNG9tMN9u2eZi9xY3L5kFpJ0o7QTsJ+KPx+WLRwTjMUDw6mIAOWxftgaqY62wruNZYSV/jEMUkcKtxF+s9UUQlvE6cexe7QibUwIASzS2fmQxrSRv/jwkGn+YpamIVhjK4IOMCxmSOIS8lj/NTzpfAIE6IBIpTkBbrxvRbA5mFaYpLNmuU1p0WKPyhMD9Y/FvsMTvJqtDc/47Jnkx47hqDrNgcnr382S5bNCSOSIt18/NvDeF7l/bnT5/24dW1WZDwCQWpy/nVd0MM3W9y0xfWuotbPg9THWMlmDuYqKiMg8o4qI3ReJ2NNLsa8bsh9NWxDM/R270COExFeihMmhlgBYqFeKhWg2h0ekF9jBtrxlFajWbcds24bSbZFVAXBS9NNFg0wkob3ye2D/92zo1c3/96MqIzuum3Js4kEQkUSqkk4FUgFygAbtVa17RTrgBowJrAHdJaj+q+Wn69tFgXps8KFEVpiqiAtR/CLveuTski+8P3X8Qb9RHRXs2P3gjjd8Dsm2wYzmj+ePkfZPFQN+uV4OHn3xrC/ZPO4aNtI9ly8Nt8Wfc3NueuYHNfg9hmTf5eK6/QsP2aS7a0v3eGCfid4LdDwGFNBa2NUVTHQlUsVMdaawQq4hUVcYoSt7v1SGU20qsa+pdqzimzBqb7t+wEuyMLnr/S4NOhVoAYkDiA7w37HhOzJ0rqdHFKItWieBT4WGv9G6XUoy2P/+s4ZSdqrSu7r2onLi3OhRlIQZkGhanWRSGnXLMmoYZKb+UpDWj/ZcV6Pqn8I4bSPPgvk9Q6eHyqjepYgzmX/JZ+8f0662OIbyg9zs3/d2EOkANMYF/tPt7d9z4f7v+Yz6L28NlQq5wjdHiVsia+idb1CVF+jStI680dsPY2GFIIiQ3WYra2AvaWbWNNaHu59zlgX4bVelgxSFEZb30xGZQ0iHuH3cvEPhIgROeIVKC4Abi05f584FOOHyh6rJQYF0rZcAUSKE6twgRyymHNedY4xckGinc2lvLE2t9ij/Vzwxea/P2a564y2JmluO3ce5iYfWmnfg5xavol9OM/R9zPf464n9LGUpaWLGVn9U4ONh2ktKmUvY1l+MLH7rfcHmVqEpqsGXSpdZrUOmsRXNg4kvaiIl6xN1NRknxk2mqv6F7c3fdqru57NeclnifrDESnilSgSNdalwForcuUUsfbD1ADi5RSGnhOaz3veG+olJoBzADIzu66DXjactgMkqOd4E/DH1/NoURr+0SwAsW43uO+8Xt+tO0QD7/3Mq7eW8mo1tyy3OTLAYqPhxtkuy7ksbH3d/bHEJ2oV0wvbht421HPaa3xhX00BhppCDRQXFdNjbcRX9iLL+TDG2qiKVxNbaCCCm85ZU1lFCceZFd7WVJbpEWlMSllKENThjIyfSR5qXnSejgL5OdHZjFtlwUKpdRioL2Rsx9/g7cZp7UubQkkHymldmitl7VXsCWIzAMYNWrU8TdW7mSpsW5KfH0gfgdFaYqcQ9apd9bs/Mbv9cGWMv7z1RU4ct4CrbnnA5OQ3dqu1EEcL075nXxTPA0ppVrXJKRGpdIv4eu7DbXW1PnrKG0qpSnYhMNwWNvDGnaS3cln3GZT4sQ89dRTETlvlwUKrfXlx3tNKXVIKZXZ0prIBMqP8x6lLT/LlVJvAmOAdgNFpKTFuthVl00UUJiqGL1T4wp8s5lPWlvbmP7ug524Mt7HcDRwySbN0ELNnycb1MQqfj3uxyRHJXbdBxE9ilKKBHcCCe6ESFdFCCLVVn0buKPl/h3Av75aQCkVrZSK/f/bu/MgKcozjuPfZ4Zd7mOXZRFYwOXaCMgVVBCUiKCAKTFBLI1nNIiKQS015REPSkViohVIJdHEKzGGlCheiAQCRqNGEEXighwrx3LJckoICsvOkz+6N6yKEzCz0435faq6+mCn+E3vzD7T7/tOv9XbwGlAadYSHqJebZvVGCIbnNCirbBq5yr2Ve1L/2CCIbDXT1vMfbOWU6/+SnLz5tN4j3PR3BTLioLZuAa0HsgZHYfX8jMRETm4qArFJGComa0Ehob7mFlrM5sZ/kxL4HUzWwwsAF5y91mRpE3j/BPaUT/ZFKtsSHlh0CzUvsKp8ipWf7w67WM/+vhTzv/tfKa/u4Hmto12rR8B4OK/pKi/Dx4anqReTgPu6H+7mpxEJDKRFAp33+bup7p753C9PTy+0d1HhNur3L1nuHRz93uiyPrfFDapx7hTOlG5tw0VzYIhi+0rgn6K55cu/NLHzVu2meGTX2Ph2h2UWDnfaTmRTbnQ88MUJy9xnutvbCgwxvcen/W7k4qI1KRhEhlw2cBiGlCEm1He4sDIp2dKF7C/KvWZn923P8XdM5Zy6eML2bGnklGJ15jU6C6mNatDvb3O5bNSrG8O009M0L159y+MoBERyTYVigyol5PkrG7BTHJrC432FYA7u6rWMWXuSlKpoHCs2rKbsx98k4dfX01DPuH+nF8xKfdBbm3RgpQZ572aovkueHBEEnJzmDBgAsnE12/+XRE5suheTxlyTo/jeGptcCuPoe85+f+EbQ02MWVeGQvWbGfIMS25f/YKPqmsoput4Rc5U2hvm/lB4/6sq7uBknXO8HecmX2NFUXGlceOoUtel6ifloiICkWmdGjagaTVYW1hJRD0U2zvtJuSOkt5a1VX3lq1nQQpLk++xA11nmIHjRll4yjLn0nufufKmVVUNIWpgxJ0ataJMceOifgZiYgE1PSUITnJHDo160h5+D2odluC9Y8bPcAPk9Npa5t5Mmcit+RMZV6qD6dX3snKNm9iif2c/XqK1tvhoeEJKusmuWvAXZrvWkRiQ1cUGdQlrwvLdyxnS5MDI59eyOnIT3OeY25htwAAB2hJREFU5jp/hj3U5cbKy5lWdRL5Hf5IIncrPVelGPl3Z14P4/3iBJd0vYjuBd0jfiYiIgeoUGRQSX4JL656kfJCo11YKF6s04Hd+wZzevJtHtg/mnJvSe8e8ymrXELhjuDOsOWF8NjQBMVNi7mq11URPwsRkc9S01MGdc7rDATf0G6zDersd5rnb2MO/bi28mr2NGzH1WfspazyWeruc26YHsxuf/93k9Rp2IjJp0zWpPQiEju6osig6lFKa1sYSXeKtsKao8oZe9Zyhhx1CSv2zOG+t38C7ox9OUW7Crj3nASb84wpA++luGlxxM9AROSLVCgyqKB+AYX1CylrvZmUwffnVDFpdJLff/Aob256jbKdZeDO6L+lGLjUmTooweKOif/MQiYiEkdqesqwUV1GsaWZMfnMBJ03wh1Tq2i8xynbWUaDT50bn0kx+g3n1e7Gs/2NQUWDuKLnFVHHFhH5UrqiyLCxPcay9ZOtTGMan+bC9c+muPPJKp4YnODS2SkKdsFjQxK83NfoVtCdiSdN1IQzIhJr+guVYclEktv63caYY8ewqFOCieckKNgFtzyVIqcK7jg/ycvHJRjV5WweH/Y4TXKbRB1ZRCQtXVHUAjNjfJ/x5NXL4z7uY8L3jJNKU0wfkGBf43rc3e82RnYaGXVMEZFDokJRiy7seiEt6rfg3gX38rtW2+nWvBsTTpxASX5J1NFERA6ZCkUtG1Y8jMHtBrO7cjd5dfM0AZGIHHFUKLIgN5lLfjI/6hgiIl+JOrNFRCQtFQoREUlLhUJERNJSoRARkbRUKEREJC0VChERScvcPeoMGWdmW4C1X/HhBcDWDMbJhDhmAuU6HHHMBPHMFcdMEM9cmczU3t1bHOwfvpaF4n9hZgvdvW/UOWqKYyZQrsMRx0wQz1xxzATxzJWtTGp6EhGRtFQoREQkLRWKL/pN1AEOIo6ZQLkORxwzQTxzxTETxDNXVjKpj0JERNLSFYWIiKSlQiEiImmpUITMbJiZLTezMjO7KcIcj5pZhZmV1jiWb2ZzzGxluM7Lcqa2ZvaKmX1gZkvM7JqY5KpnZgvMbHGYa0IccoUZkma2yMxmxCjTGjN738zeM7OFMcrVzMyeNrNl4Wusf5S5zKwkPEfVyy4zuzYm5+q68LVeamZTw/dAredSoSB4UwO/BIYDXYHzzKxrRHEeB4Z97thNwFx37wzMDfezaT9wvbsfA/QDxoXnJ+pce4HB7t4T6AUMM7N+McgFcA3wQY39OGQCOMXde9UYex+HXJOBWe7+DaAnwXmLLJe7Lw/PUS/gm8Ae4NkoMwGYWRtgPNDX3bsDSeDcrORy9//7BegP/LnG/s3AzRHmORoorbG/HGgVbrcClkd8vp4HhsYpF9AAeBc4IepcQFH4hh0MzIjL7xBYAxR87ljU56oJsJpwYE1cctXIcRrwRhwyAW2AdUA+waRzM8J8tZ5LVxSB6l9AtfXhsbho6e6bAMJ1YVRBzOxooDcwPw65wiae94AKYI67xyHXz4EfAakax6LOBODAbDN7x8wuj0muDsAW4LGwqe5hM2sYg1zVzgWmhtuRZnL3DcDPgHJgE/Cxu8/ORi4VisDBJrLWuOHPMbNGwDPAte6+K+o8AO5e5UETQRFwvJl1jzKPmX0bqHD3d6LM8SUGuHsfgibWcWZ2ctSBCD4Z9wF+7e69gX8RXbPcZ5hZLnAmMC3qLABh38NIoBhoDTQ0swuy8X+rUATWA21r7BcBGyPKcjCbzawVQLiuyHYAM8shKBJPuvv0uOSq5u47gb8S9O9EmWsAcKaZrQH+BAw2sz9EnAkAd98YrisI2tyPj0Gu9cD68EoQ4GmCwhF1LggK6rvuvjncjzrTEGC1u29x90pgOnBiNnKpUATeBjqbWXH4KeJc4IWIM9X0AnBxuH0xQR9B1piZAY8AH7j7AzHK1cLMmoXb9QneSMuizOXuN7t7kbsfTfA6mufuF0SZCcDMGppZ4+ptgrbt0qhzuftHwDozKwkPnQosjTpX6DwONDtB9JnKgX5m1iB8T55K0PFf+7mi6CCK4wKMAFYAHwK3RphjKkH7YyXBp63LgOYEnaMrw3V+ljMNJGiK+wfwXriMiEGuHsCiMFcpcHt4PNJcNfJ9iwOd2VGfqw7A4nBZUv0ajzpXmKEXsDD8PT4H5EWdi2BwxDagaY1jcThXEwg+DJUCTwB1s5FLt/AQEZG01PQkIiJpqVCIiEhaKhQiIpKWCoWIiKSlQiEiImmpUIiISFoqFCIikpYKhUgtM7O7qufwCPfvMbPxUWYSORz6wp1ILQvvuDvd3fuYWYLgG7THu/u2SIOJHKI6UQcQ+bpz9zVmts3MegMtgUUqEnIkUaEQyY6HgUuAo4BHo40icnjU9CSSBeFdid8HcoDO7l4VcSSRQ6YrCpEscPd9ZvYKsFNFQo40KhQiWRB2YvcDRkedReRwaXisSC0zs65AGTDX3VdGnUfkcKmPQkRE0tIVhYiIpKVCISIiaalQiIhIWioUIiKSlgqFiIik9W9t4UL9391TSgAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "ib, ie = 385, 467\n",
    "uveln19 = restart_file.un[0, pz]\n",
    "uvelb19 = restart_file.ub[0, pz]\n",
    "uvelb19[ib:ie, 1].plot(label=\"Before 1\", linewidth=5)\n",
    "uveln19[ib:ie, 1].plot(label=\"Next 1\")\n",
    "uvelb19[ib:ie, 3].plot(label=\"Before 3\", linewidth=5)\n",
    "uveln19[ib:ie, 3].plot(label=\"Next 3\")\n",
    "plt.plot([prob-ib, prob-ib], [-0.5, 2], color='k');\n",
    "plt.legend();"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYoAAAEGCAYAAAB7DNKzAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nOzdd3hUVcLH8e+ZmfRGKgESIISSUEISQKoISFGBoKKIsNRdEQuKy+qKbbGgroIIKiKvsohgpSMiSBepAQIEEghJIAkhvfeZuef9YzCChBAwDXI+z8MDM/fOvedOwvzmniqklCiKoijKtejqugCKoihK/aaCQlEURamUCgpFURSlUiooFEVRlEqpoFAURVEqZajrAtQEDw8P2bJly7ouhqIoyi3j8OHDGVJKz4q23ZZB0bJlS8LDw+u6GIqiKLcMIcT5a21TVU+KoihKpVRQKIqiKJVSQaEoiqJU6rZso6iI0WgkKSmJkpKSui7KLcPW1hYfHx+srKzquiiKotShBhMUSUlJODk50bJlS4QQdV2cek9KSWZmJklJSfj5+dV1cRRFqUN1VvUkhPAVQuwQQkQJIU4KIZ6tYB8hhFgghDgrhDguhAi92fOVlJTg7u6uQqKKhBC4u7urOzBFUer0jsIEzJBSHhFCOAGHhRC/SClPXbbPvUCbS3+6A59e+vumqJC4Mer9UhQF6vCOQkp5UUp55NK/84EooNmfdhsBLJMW+4FGQogmtVxURVGUemH69OlMnz691s9bL3o9CSFaAiHAgT9tagYkXvY4iavD5PdjTBFChAshwtPT02uimH+ZXq8nODiYzp07Exoayt69e6/7mgULFhAYGMjYsWOrtSyZmZn0798fR0dHnn766Wo9tqIoNSMiIoKIiIhaP2+dN2YLIRyBVcB0KWXenzdX8JIKV1qSUi4GFgN07dq1Xq7GZGdnV/5D3rx5MzNnzmTXrl2VvmbhwoVs2rSpyg3KJpMJg+H6P1ZbW1vefPNNIiMjiYyMrNKxFUVpmOo0KIQQVlhCYoWUcnUFuyQBvpc99gGS/8o5W7648a+8vErOvTv0uvvk5eXh6upa/vj999/n+++/p7S0lAceeIDXX3+dqVOnEhcXR1hYGJMnT2bChAlMnjyZuLg47O3tWbx4MUFBQcyaNYvk5GTOnTuHh4cH8+fPZ+rUqSQkJADw4Ycf0rt37yvO7+DgQJ8+fTh79mz1XryiKLedOgsKYWkp/QKIklJ+cI3d1gNPCyG+xdKInSulvFhbZaxuxcXFBAcHU1JSwsWLF9m+fTsAW7ZsISYmhoMHDyKlJCwsjN27d7No0SJ+/vlnduzYgYeHB9OmTSMkJIS1a9eyfft2xo8fX36HcvjwYfbs2YOdnR1jxozhueeeo0+fPiQkJDBkyBCioqLq8tIVRbmF1eUdRW9gHHBCCPF7pdtLQHMAKeUi4CfgPuAsUARMqoNyVpvLq5727dvH+PHjiYyMZMuWLWzZsoWQkBAACgoKiImJoW/fvle8fs+ePaxatQqAAQMGkJmZSW5uLgBhYWHY2dkBsHXrVk6d+qPzWF5eHvn5+Tg5OdX4NSqKcvups6CQUu6h4jaIy/eRwFO1U6La1bNnTzIyMkhPT0dKycyZM3n88ccrfY3l7bjS711YHRwcyp/TNI19+/aVB4eiKMpfUS96PTVE0dHRmM1m3N3dGTJkCEuWLKGgoACACxcukJaWdtVr+vbty4oVKwDYuXMnHh4eODs7X7Xf4MGD+fjjj8sf10UvCUVRbh913uuptlWlobmm/N5GAZa7gy+//BK9Xs/gwYOJioqiZ8+eADg6OrJ8+XK8vLyueP2sWbOYNGkSQUFB2Nvb8+WXX1Z4ngULFvDUU08RFBSEyWSib9++LFq06Kr9WrZsSV5eHmVlZaxdu5YtW7bQvn37ar5qRVFudaKi6oxbXdeuXeWfFy6KiooiMDCwjkp061Lvm6LUH/369QMsNQrVTQhxWErZtaJtqupJURRFqZQKCkVRFKVSKigURVGUSqmgUBRFUSqlgkJRFEWplAoKRVEUpVIqKGqREIIZM2aUP54zZw6zZs26qWMtXbqU5OSK50f84Ycf6NChAzqdjj93E1YURblRKihqkY2NDatXryYjI+MvH6uyoOjYsSOrV6++aq4oRVGUm6GCohYZDAamTJnCvHnzrtqWnp7OyJEj6datG926deO3334DYMSIESxbtgyAzz77jLFjx7Jy5UrCw8MZO3YswcHBFBcXX3GswMBA2rVrV/MXpChKg9DgpvAAYNOLkHKieo/p3Qnuffe6u/0+tcYLL7xwxfPPPvtshVODL168mN69e+Pn58fcuXPZv38/bm5ufPzxx8yZM4euXSscSKkoilJtGmZQ1CFnZ2fGjx/PggULrpjd9VpTgzdu3Jg33niD/v37s2bNGtzc3Oqi2IqiNGANMyiq8M2/Jk2fPp3Q0FAmTfpjeY3KpgY/ceIE7u7u12yTUBRFqUmqjaIOuLm5MWrUKL744ovy5641NfjBgwfZtGkTR48eZc6cOcTHxwPg5OREfn5+7RZcUZQGSQVFHZkxY8YVvZ8WLFhAeHg4QUFBtG/fnkWLFlFaWspjjz3GkiVLaNq0KXPnzmXy5MlIKZk4cSJTp06tsDF7zZo1+Pj4sG/fPoYOHcqQIUNq+/IURbmNqGnGlUqp901R6g81zbiiKIpSL6mgUBRFUSqlgkJRFEWplAoKRVEUpVIqKBRFUZRKqaBQFEVRKqWCohbp9XqCg4Pp3LkzoaGh7N2797qvWbBgAYGBgYwdO7Zay3Lw4EGCg4PLy7NmzZpqPb6iKLePhjmFRx2xs7MrH3G9efNmZs6cya5duyp9zcKFC9m0aRN+fn5VOofJZMJguP6PtWPHjoSHh2MwGLh48SKdO3dm+PDhVXqtoigNS4P7VOj0ZacaP8eJCdefmTYvLw9XV9fyx++//z7ff/89paWlPPDAA7z++utMnTqVuLg4wsLCmDx5MhMmTGDy5MnExcVhb2/P4sWLCQoKYtasWSQnJ3Pu3Dk8PDyYP38+U6dOJSEhAYAPP/yQ3r17X3F+e3v78n+XlJQghKimq1cU5XbT4IKiLhUXFxMcHExJSQkXL15k+/btAGzZsoWYmBgOHjyIlJKwsDB2797NokWL+Pnnn9mxYwceHh5MmzaNkJAQ1q5dy/bt2xk/fnz5Hcrhw4fZs2cPdnZ2jBkzpsIpy//swIEDTJ48mfPnz/PVV1+puwlFUSqkPhlq0eVVT/v27WP8+PFERkayZcsWtmzZQkhICAAFBQXExMRctULdnj17WLVqFQADBgwgMzOT3NxcAMLCwspnnr3WlOVOTk5XHK979+6cPHmSqKgoJkyYwL333outrW3NXLyiKLcsFRR1pGfPnmRkZJCeno6UkpkzZ/L4449X+pqK5uX6vcrIwcGh/LnKpiyvSGBgIA4ODkRGRqqFkBRFuYrq9VRHoqOjMZvNuLu7M2TIEJYsWUJBQQEAFy5cIC0t7arX9O3blxUrVgCWScE8PDxwdna+ar9rTVl+ufj4eEwmEwDnz5/n9OnTtGzZsjouTVGU20yd3lEIIZYAw4A0KWXHCrb3A9YB8ZeeWi2lfOOvnLMqDc015fc2CrDcHXz55Zfo9XoGDx5MVFQUPXv2BMDR0ZHly5fj5eV1xetnzZrFpEmTCAoKwt7eni+//LLC8yxYsKB8yVWTyUTfvn1ZtGjRFfvs2bOHd999FysrK3Q6HQsXLsTDw6MGrlpRlFtdnU4zLoToCxQAyyoJin9JKYfdyHHVNOPVR71vilJ/NMhpxqWUu4GsuiyDoiiKUrlboY2ipxDimBBikxCiw7V2EkJMEUKECyHC09PTa7N8iqIot7X6HhRHgBZSys7AR8Daa+0opVwspewqpezq6elZawVUFEW53dXroJBS5kkpCy79+yfASgihWlwVRVFqUb0OCiGEt7g0UEAIcQeW8mbWbakURVEalrruHvsN0A/wEEIkAf8BrACklIuAh4AnhBAmoBgYLeuym5aiKEoDVNe9nh6VUjaRUlpJKX2klF9IKRddCgmklB9LKTtIKTtLKXtIKa8/L3c9JoRgxowZ5Y/nzJnDrFmzbupYS5cuJTk5ucJtr776KkFBQQQHBzN48OBr7qcoilIV9brq6XZjY2PD6tWrycjI+MvHqiwonn/+eY4fP05ERATDhg3jjTf+0hhFRVEaOBUUtchgMDBlyhTmzZt31bb09HRGjhxJt27d6NatG7/99hsAI0aMYNmyZQB89tlnjB07lpUrVxIeHs7YsWMJDg6muLj4imNdPq1HYWGhmkJcUZS/pEFOCpjy9tuURkVX6zFtAgPwfuml6+73+9QaL7zwwhXPP/vssxVODb548WJ69+6Nn58fc+fOZf/+/bi5ufHxxx8zZ86ca07i9/LLL7Ns2TJcXFzYsWNHtVyjoigNk7qjqGXOzs6MHz+eBQsWXPH81q1befrppwkODiYsLKx8avDGjRvzxhtv0L9/f+bOnYubm1uVzjN79mwSExMZO3bsFRMEKoqi3KgGeUdRlW/+NWn69OmEhoYyadKk8ucqmxr8xIkTuLu731Sj9JgxYxg6dCivv/76XyqzoigNl7qjqANubm6MGjWKL774ovy5a00NfvDgQTZt2sTRo0eZM2cO8fGWiXSdnJzIz8+v8PgxMTHl/16/fj0BAQE1cRmKojQQKijqyIwZM67o/bRgwQLCw8MJCgqiffv2LFq0iNLSUh577DGWLFlC06ZNmTt3LpMnT0ZKycSJE5k6dWqFjdkvvvgiHTt2JCgoiC1btjB//vzavjxFUW4jdTrNeE1R04xXH/W+KUr90SCnGVcURVHqPxUUiqIoSqUaVFDcjtVsNUm9X4qiQAMKCltbWzIzM9WHXxVJKcnMzMTW1raui6IoSh1rMOMofHx8SEpKQq1+V3W2trb4+PjUdTEURaljDSYorKys8PPzq+tiKIqi3HIaTNWToiiKcnNUUCiKoiiVUkGhKIqiVEoFhaIoilIpFRSKoihKpVRQKIqiKJVSQaEoiqJUSgWFoiiKUikVFIqiKEqlVFAoiqIolVJBoSiKolRKBYWiKIpSKRUUiqIoSqVUUCiKoiiVUkGhKIqiVEoFhaIoilKpOg0KIcQSIUSaECLyGtuFEGKBEOKsEOK4ECK0tsuoKIrS0NX1HcVS4J5Ktt8LtLn0ZwrwaS2USVEURblMnQaFlHI3kFXJLiOAZdJiP9BICNGkdkqnKIqiQN3fUVxPMyDxssdJl567ihBiihAiXAgRnp6eXiuFUxRFaQjqe1CICp6TFe0opVwspewqpezq6elZw8VSFEVpOOp7UCQBvpc99gGS66gsiqIoDVJ9D4r1wPhLvZ96ALlSyot1XShFUZSGxFCXJxdCfAP0AzyEEEnAfwArACnlIuAn4D7gLFAETKqbkiqKojRcdRoUUspHr7NdAk/VUnEURVGUCtT3qidFURSljqmgUBRFUSqlgkJRFEWplAoKRVEUpVIqKBRFUZRKqaBQFEVRKqWCQlEURamUCgpFURSlUiooFEVRlEqpoFAURVEqpYJCURRFqVSV5noSQvQCWl6+v5RyWQ2VSVEURalHrhsUQoivAH8gAjBfeloCKigURVEagKrcUXQF2l+ayVVRFEVpYKrSRhEJeNd0QRRFUZT6qSp3FB7AKSHEQaD09yellGE1VipFURSl3qhKUGwAvgASAFGzxVEURVHqm6oEhRPwdyAL+BZYKaVMrdFSKYqiKPXGddsopJSvSyk7YFmStCmwSwixtcZLpiiKotQLNzLgLg1IATIBr5opjqIoilLfXDcohBBPCCF2AtuwNGw/JqUMqumCKYqiKPVDVdooWgDTpZQRNV0YRVEUpf65blBIKV+sjYIoiqIo9ZOaFFBRFEWplAoKRVEUpVIqKBRFUZRKVWmacUVR6kZGcQZ5pXk4WjviZO2Erd4WIdQECUrtUkGhKPXIudxzHLh4gKPpRzmaepTkwuQrtht0Bnwcfejr05dBLQYR5BmETqiKAaVmqaBQlDqWU5LDpnOb2BC7gRMZJwDQmyX+F2FEosQ9T1JiDSXWghKrMhI941iRHc+yU8vwsvPiHr97mNxxMu527nV8JcrtSgWFotSRk5knWRq5lK0JWzFpJhyKJYOiJD2jJW0uSGxMlv3ybcHGBNamP5aEKbCFo/6Cg21TWFHwJWvPruXZ0Gd5qO1D6g5DqXYqKBSlFkkp2Zu8l/9F/o8DKQcQUhIUJ+l/QtLtjMTKDAkesDVEcMpXEO0ryLe3tEnoNIldKQQkSe44I+kSI7nzpCTJHb4YksObZW+y9uxaXunxCu3d29fxlSq3kzoNCiHEPcB8QA98LqV890/b+wHrgPhLT62WUr5Rq4VUbkkms8be2ExOXMilsNREUZmZEqOZZo3sGNLRm7aNnWq1PIXGQjbEbuC7099xNucserPkrpOSsAMavhmQZwe/hAh2BOk43/iPxmq90OPj0IQiUxF5ZXkU6kwcbiM43MYSHKFnJRO2afzna4097SVfDTjO6IzRPNjmQZ4OeRoPO49avU7l9lRnQSGE0AOfAIOAJOCQEGK9lPLUn3b9VUo5rNYLqNyy9sRk8NbGU0Sn5Fe4fd7WMzzSrTkzBrfFw9GmRssSnRXNqjOr2BC3gUJjIQ7FkuHHJPeFa7jnwzkvWDBcx75AgVlvCYgAtwAG+A4gpHEIQR5B2FvZA5a7kaySLHYl7eKX87+w/+J+wtuaOOYnuH+fxoj9ktCzZtb0kqw3reTncz/zj07/YFz7cdjoa/Y6ldubqKulsIUQPYFZUsohlx7PBJBSvnPZPv2Af91oUHTt2lWGh4dXY2mVW0FcegFv/xTF1qg0bCjj7/pNDNXvx4ES7EUptpQSJ5vyjnEMB2QgjjYGnh7Qmsm9/bA2VF+9fk5JDhvjN7Lu7DqisqIAaHVRMuSwRu8oibUJIpsL1vUQHGslQAjcbN0Y1moYYf5htHNrV6XzZBRnMO/wPNbHrgegcZbl7qLrWUlqI1jeX8eBdoJWjfxZNHARTRybVNs1KnWjX79+AOzcubPajy2EOCyl7FrhtjoMioeAe6SU/7j0eBzQXUr59GX79ANWYbnjSMYSGievcbwpwBSA5s2bdzl//nzNXoBSbxjNGgt3xPLxjhiMZo2hugPMtPoaH5HBAS2AFOlGkbShBGsG68NpJjL52dyNt01jSJCN6dTMhYVjQ/F1s7+p80spicmJYXfSbn5N+pVj6ccwSzMOxZI+pyT9jmv4p0CJFezuKNgSqiPBy3L30NK5JZM6TmJYq2FY661v6vzhKeG8tf8tYnNjAegUrzFhm0bzdDjWUjD3QR1ubs34fPDn+Dr73tQ5lPqhIQbFw8CQPwXFHVLKaZft4wxoUsoCIcR9wHwpZZvrHVvdUTQckRdyeX7lcaIu5tFGJDHb6gvu0J0mSmvOG6Zx7NM6XLG/DWX8Q/8TTxrWYYWJ+aaRLDSH4Wxnw7xHOjMgoHGl55NSklqUyqnMU0RmRHIy8yQnM0+SW5oLgF2pJChe0iPa0jhtbYb4xrA9SMfujoJiW0tAhHqFMr7DePr79q+WXkpGs5Gvo7/ms2OfkW/MR6dJBh2RTNyqEdsE3n5Ej6NrYz4f8jl+Ln5/+XxK3WiIQXHdqqcKXnMO6CqlzKjs2Coobm8ms8ahc9lsPJHMtwcTQTPyuP5HnjGspgA73jON5gfzXWjoCPZtxN0BXthZ68krMfG/3+LJLzHhSTavWX3FcP1+tplDeM74BHk48nT/1jw3qC16neUDvcxcxtG0o+xL3sfJzJNEZ0WTU5pTXhZro8QvFdomSULiJAGJEoNm6dK6p4Olcfqct+VYdgY7hrcaziMBj9DWtW2NvDdZJVksjFjID2d+QJMa3U5rPLdW47wXvDVaj62rB+/1fY9u3t3UCO9bUEMMCgNwBrgbuAAcAsZcXrUkhPAGUqWUUghxB7ASaCGvU2gVFLcHKSWlJo3UvBLi0guJTS8g6mI+O06nkVVYBkAHcY7/Wi2mo+4cP5p78JpxIlk408TFlhfvDWB4UFN0uj8+EDMLSpm39QxfH0hAk5K/6bfymmEZKdKNJ4zPcVK2pKufPfd0T+FY5l4OpByg2FQMWEKhRRq0TJW0SpH4X5T4poP+0m9jgicc8Rccaa3jTDPQLp23k0cn7m99P/f63YuTde30tjqVeYontj5BVkkWoTEaM9ZoXHC3hEWeg8DH0Yfh/sO51+9eWjq3VKFxi2hwQQFwqTrpQyzdY5dIKWcLIaYCSCkXCSGeBp4ATEAx8E8p5d7rHVcFxa0js6CUExdyib0UBLFpBSRlF5NfYqSozIxJ+/PvpyRQJHCP/hCDdYcI1CWSLp15xTiZzdod6HWCx/u2YtqANthZ66953mOJOTy54ggXcooJETF8Yj2fUkMRzzl25axLGkJfgtAkbZIhJFajc5ykVSroLhUn3xZimwhim8DZpoKzTQS5jn982DZzbMbA5gO5v/X9tHZtXQPv3PXF5cbx2ObHSCtOIyhO44VVGgCH2gh2dRIc9xNoOoGXvRfdvbvTzbsbd/rcqbrU1mMNMihqigqK+is2vYBdp9M5mphDRGI2iVnFV2x3oYAWIhVnUYQDxThSgrvIpZW4iL8uGX+RjJsoQJOCQ7IdW8xdWWnuSy6OBHg78f5Dnenk41KlsmQXljH9+yPsSdyPg+secDqNFJbeQ/eFa9x5UuJYApqAM03hZAtBXBNBXGNBpjNw2bdwW70tnT07c6fPndzZ7E78XPzqxbf0xLxE/r7l71wsvIhvmmTwEY1eURKnEshxgKOtBJEtBCdbCLKcBQadgaF+Q5nYYWKdBZxybSooqpEKitpTWGoiOiWfpOwiXOysCPJphJvDH713couMHEvKYdeZdLZHpxGfUQiADg1/kUyw7izBIpbWugv4i2Q8RF6F50mXzsTJpsRqTTkm/dlqDiUTSyDYGHQ8fpc/T/dvXeVurlklWaw/u54fzvxAQn4CSEmH85KhhyyD2DQd7A8QHGpj+eZdaPfHh76N3oY2jdoQ6B5IR4+OdHDvgH8jfwy6+jnRQXJBMk9te4qzOWcByzxSobGSPict1+x8KasvuMGWUB3bOwtKrQV9mvVhauepdPbsXIelVy6ngqIa1UVQlBjNWOt1V9SH346klPwak8H34YmcuJDL+cyiq/Zp7maPv6cDsemFJGT9sb0Z6fTXR3C37ghddWdwEpZPqDxpT5RsTpzWhFjZlHPSm1zpQAF2FGBLrnQkD4crzmFrpeOutp4M6eDN3QGNcbG3um7ZTZqJvcl7WROzhp1JO8vnV7rrhGTQUY1mWZBrD1tCBL+E6si5VJUkzXaIknY82mkgD3fshZ+LX70NhWsxmo3lEw8euHgAieX/vZCS5mnQ4bykx2mNgCTLPFKbQwWbuurId9Axvv14poVOU4P26gEVFNWoNoOiuMzMy2tPsPrIBRys9bwxoiMju/jUyrlr24mkXN7ZFMXe2EwMmGgtkgkU5wnQJdBaJJMqXTkqW3NM8ydeNqGVSKaTLp6OIp7uuigCdYkAxGuN2aN14qjWmgjZmnjpjbzGGlpWeoGLnRUt3B3w93Sglacj7Ro70aOVe6VtEL+TUhKVFcXGuI1sit9EenE6SElgIgw4ptEz2jIA7kxT+CVEx972AqNBIDUD3fKtebYwjtSijjxvfIJ84cAzA9rw7N1tbukvBCmFKfwU/xM7E3dyIv0EJmkq39YmyTKtSLczlhlr/zdIx65OgtaubXi7z9sEugfWYckVFRTVqDaD4qt953h13UkG6w4RK5sSK5ux6G9duKejd62cvzZkF5bx2vqTbDiWjB4zI/W7+adhJd4iG4BSaUWc9KapyMRFXH2HUSBtOa61YrsWwnYthDjZtHybTkDbxk50aOpCay/HS2HggJuDDQ42emwM1w+DiqQVpbH27Fp+jPuR+FzLVGHueZK+JyT9j2t450CRtaUL6y8hf8yv5GXnRTuHe9i6348yoy0T9Zt52bCCC9KDJ4zTiZItGBjYmA8e6Yyz7fXvYuq7ImMREWkR7E3ey9rYteXjQZpmSqZsMtM+EQ62FXx2r45iRyueC32O8R3G13GpGy4VFNWoNoNi/JKD/HomlVM2k8nGkWGlb5OFMzv+1Q8/D4frH6CeKyozMfLTfURdzGWA7igvGr6hre4CR7TWLDUN4aRsyTnpjRk9Ao2WIpVgcRY/3UVitaZESj/iZJPyOwYhIMS3EQMCvOjSwo0gHxccbKqnGsesmdmbvJcfzvzA7qTdmKUZ+xJJ99OSOyMl7RMspYhsLtgRJDgQICizsgREr6a9GNV2FH19+2Kls+J4Ug5PLLf0iuoiTvOJ9QJcKOQF4xQ2aL1o5enA4nFdae3lWC1lrw+KjEWsPbuWZaeWcaHgAkJKhh6UPLpLo9AGPhmu41grHY91eoxpIdPqRWN9Q6OCohrVZlD0fW8HJVkXOGj7FAC/mjsywfgibb1dWPNk7ypVj9Rn87fGMG/rGWYYvmeaYS2xWhPeMz3CZq0bYPmgaO5mT2svR5Kyi4hJK+DyXykhwN/TkU7NXOjT2oN+7Txxr+aJ+DSp8VP8TyyMWEhifiJISygMjJDccdoyOjrZDX7toGNPB0Gqq6XcHnYehPmH8VCbhyqc2iK7sIxnv4tg95l0PMjlE+v5dNdF85lpKO+ZRmNnY8OwoCZ0a+nGHX5u+Lja3RYfnibNxLJTy/jo6EeYNBO+6ZJp6834ZMBHw3Xsa6/jb4F/44VuL9wW13srqauguLVa5OoZTZOk5JbQSaQB8KO5O8P0B5guV/JByiheWRvJnIeDbtn/TGn5JXy2OxYfkcYU/Y+sN/fkn8YnMGGgpbs9Mwa3o3+AF46X3REUlJo4kZRLWn4JTVzs6NDUudruGP5MSsmvF35l/pH5nMk+g02ZZGiEZOClhukCW9gWLNjdUUdsE0BYun8O9OnHA20eoFfTXpU2Srs6WPO/id2Yv/UMC7af5W9lL/GK4SseN2ykvTjPq2WT+PaQkW8PWdpeDDpLe4qLnRVuDtaEtnBlSIfGhPi63lJtGu3WLqkAACAASURBVAadgckdJ9O7aW9e2vMSZzjDf/6m598/mHl2vYaNCZaznBJzCa/2eFUtlNQAqKD4C9LySykza/jq0gGYZ3qIQmnHM4a1HNHasOoIBHg78VjfVnVc0psz75cYisrMvGv1HRo6ZhvHgs6KN4a359E7mmOlv/oDwtHGQE//ml2S06SZ2JawjWUnl3E84zh6s2RwhOSh3zQaFcLpZvDJUMvU3b9XLYV6hTK01VAGtxhMI9tGVT6XXif45+B2dGjmwj+/i+A/ZZOIlH68ZVjCTpsZpMpGHNICCNfaclG6kVfsQE6RI/EZjQg/n83i3XF4ONrQu7U7VnodmpRomqTYaKaw1ExBqYkSoxlXe2t8XO3wcbWnubsd3f3cadrIrqbewipp59aOb4Z+w3uH3uO709/x9iN6/rVK48mNGrZlsJKVlJnLeKPXG+h1t/ads1I5FRR/we9dP30v3VEkSU9eM02koy6eD60WMrJsFrN/gsTsIl4b1h5DBR+s9VVMaj7fHUogSMQSpt/HAtP9pOLGxB4tGN+zZZ2UKb8snzUxa1gRtYLkwmR0mmV21lG/WhqnT/nCnAf1nPG51DBt78XINiMZ0XoEzRyb/aVzD+ngzdqnejPlq8P8kNGP/Vogd+mO0013mm66aIbp91/1mhNaS7aYu/JLYVfWRZTwe1UdgDVGHCnGSRThQAmx0oX9uFzRltOzlTsjQ30Y0tEbKSX5JSayi8qISS3gZHIuJ5PzyCosw9/TkddHdKiRtTWs9da83P1lHKwcWBK5hPce0jF9ncbkXzSsTLCe9UgpebP3myosbmMqKP6CxPKgSCdNNqIUy0CzJ4zTWWX9H9ZY/4enjM+wbB/EZxTy8aOhVervXx+8uykaTUpetl5BunTmM9NwnGwNPHP3dSfvrXanMk/x/env+Sn+J4pNxdiVSIYdk9wbruGZZ1n85+1ROiIure3Qs0lPxgSOoU+zPtU63qFNYyfWPtWbNzacYm2EYLm5McvNgwCJFzl4iFxcRCEuFOInUhioP8xzhlXMECsplZafu0BDh8QgtKuOXyoNXJTuxEtvftK681Nsd/bGZjLjh2MVlseNPJqLNH5O8SM6JY+Nz9yJrVX1f1gLIZgeOh1bgy0LIxYy734dT2/QGLfDcg0b2ICGxuzes1VY3KZUUPwFidl/BEWi9Cx/PkE25v7SN/ncei5Lrf7LG6bxfBkzmJA3t9DY2ZYmLrb4eTgytkdzQpu7VulcWYVl7I/LJPxcNsVGE34eDnRo6kKHps40sr+5dQyklOSVmLiQXUxyTjEXcix/n88sYlt0GoN0h+mui+Zl42QKsePF/q2vGHVdk4pNxWw+t5nvT3/PiYwTADTLsAyM639cYldm6b20ZLDgSGuBFILOnp15NvRZunl3q7FyudhZMXdUZ14b1p7w81kcjM/i4Lksoi7qSTO6wmUN+Z+aw/Akh7v1R2gpUpAINAQSQbG0IR878qU9RdjgKXJpJjJoJjLoJOJ432oxbxiW8rPWjUNaAI4U0UgU0ujSFCdtdYl4XhrFvsB0Px+kj+KLPfE81b9mpt0QQvBE5yew1dvyweEP+ChMB1wKCwEb2IiUktl9Zt9ygxGV61M/0b/g93mKfHVphGuWaaNfHdaes2n5fHMQRpbN4kOrT3jd6kuCdLHsNXckOd+d5Dx3jid4sDbiAgtGhzA0qMllxyxie3Qa6fml5BYbyS02EpNWQNTF36e2+P2T6I9qjDZejjzRz58HQppdt+G8sNTEztPpbD6Zwu6YdHKKjFfto8dMU7J50fANZ7WmfGvuT7NGdkzs1fKm36uqis+N5/vT37Mudh35ZfnozZKeZyxrK3RMkJh0sC9Q8GM3HfFN/lg69MnOT9LPt1+tdRxwsbfi7sDG3B34x/oVJUYzucVGcoqMRCRms/lkKntidHxrHnCDR5eEihhG6n9lmH4fD+h/A6BM6snFgQvSg+3mUM5IH3roTvG4fiPfm/vzyQ49D3fxwcvZthqv9EqTOk7CoDPw3qH3/giL7b/fWfxEmbmM//b9700vwqTUTyoo/oLE7CIMmGhCJonSCwA/D3sm925JGy8n3tp4iseNz/Ev+T1T9RsYqd9T/tpYrQljy17imW+PohPQP8CLeb+cYfGvcfy5x7IjRQzRneQu3THu0h/DjXyiZHNOaH5ESj92pIXwz+8L2BSZwjsPdrqqrtpo1tgencYP4UnsjkmnzKQh0AgWsbTTJ9JMZOAj0mkmMmgqMvEmq7xqZHLZvzCj5/kh7WqkWgMs3Vv3XNjDiqgV7E22TA7snit5JELj7mOSRoWQ5gIr+unYESTIc7CEQT+ffozvMJ6ujbvWi55ltlZ6bK30NHa2pZ23E490a05BqYn9sZmkF5SiFwKdTqDXgY1Bj4ONAUcbPdZ6PWn5JSRmFZGQVcz26FSOZLbliKktr5vG40EuOThShA2Xf0EA2GjuznabfzHTsIKnyqbz3ubTzHm4ZudmGtd+HDqh492D7/JRmA5xKSysTLBa/sKT2/KZ338+Dla3/jgixUIFxV+QlFVEE5GJXsjyqidfV3uEEEzu40crTwemfX2U90tHs8D0IN4ii6YiEz+RwkzD13xtPZvRZa8y7Zuj+LrZE59RSKA4z3NWK2kiMi81dhbTiAIMQiNf2vGb1pGL0o32uvM8qN/DBPELOQYHXjFO5sdTPTlyPpvpA9vg7miDQSc4lpTD9+FJpOeXYsBET90phhgOMVh/GC9hWYDHJHVclO5cwIMDWiAXpDvJ0oNorTkRsjVBPi6EdW5a2VtxU0pMJayOWc3yqOXl4x86nZPcc1jS5awECUdbC7aECCJaCaRO4GjlyBj/MEYHjL4lVmpztDEwsH3lq+ZZ/DHj7avDAjmSkM3Kwxf45VQKyQVW2Fnp8bI14GRrwNvFlvZNnEnIKmLzSfjUNJwZVitZZopi5WEY16MFnX2r3rPrZowNHAvAuwffZUGYDqNeY/RuDftSwfL++/nH5n+wcOBCXG2rVrWq1G9qwN1NKjWZCXj1Z3qKSL62fptHy15mn9aBqDfuuWKQXVx6AZ/siGXn6TQyLy22AxAqzrDM+l1SpSujy14hByeeMqzlKf06cnHguNaKfOwpkHZk4MJecweOyDYYL8t2gUZ7kcBsqy8I1sWyztyLV40TyePK0cK+IpXR+h2M0u/CU+RSJG3YqXVms7kr4Vo7UnDDTMV3CwHeTnw+oSs+rje3nnRFioxF/HDmB5aeXEpGcQY2ZZaJ+e45rOGTaZmYb1uwYGuwjgwXyzfoQLdAHmn3CPf63Yu9VfWV5VZg1mT5inuXyy4so9+cnZQUF7LN5l/kSQeGlc0mpIU7K6f2rJW7rG+jv2X2gdkIKZm0ReOeI5KtwYL/G6LDz9WfT+7+BB+n23Pus7qgBtzdYpJzSpCS8jEUidITD0ebq0Zit/J0ZO4oS1VAidFMck4xr607yZ6zMLHsBb60/i/fWM/GiIFAXQJrzL153TieHK5eCa21lyO9/N1p4mLH6ZQ8TibncTKtJSPLZvGkfh3PGNbQ3SaKCK01JnSY0ONFDr30pzBLwXYtlB/MfdmldS7voWWt1+HbyJZmrnY0cbGjWSM7mrna4dPIDl83+2odbZxfls93p7/jq1NfkVWShUOx5KFwy9oPjiUQ6w0fDbOMfzAZBDZ6G+73u5dRbUfR0aNjvaheqgsVhQRYBgT+c1Bb/rP+JO8Yx/CJ9QJG6Xfy7fkBLNwZy5P9/Gv8PRsdMBoHKwde/e1Vlgw2UWirMXKvxL5U45NhsYzZOIZ5/efRpXGXGi2HUrNUUNykxMvGUPxeddPZrfIBUrZWelp5OvJ/47vy9y8PsTc2gEllL7DU+j3yseMfZTPYqnXBSi+Y1KMFzd3scbGzwtXemg5NnStspPw58iIvrYnko8IH2aEFM9PwDS1EClaY0WOmDCvmGh/iB/NdpGAZCOfhaM2D7b0Z0qExvfw9qryGw83KLM5kedRyvo3+lgJjAY5FkkcOadx7WGJfapl0bn13y/KhCIGvky+j241mROsRuNhUbRGihmps9+Ys33+ejWndGa8F8C/D9xzS2vH+ZohIzOH9h4JuuldcVQ33H46ztTMzds3gu7tKKbTVGL9dwyPXzPsPZfGPLf/gtR6v8UCbB2q0HErNUVVPN2nFgfO8vCaS+VYfEypiuLNsPiOCmzJ/dEiVXl9cZmbS0oPsj8vCm0zysacQOwKbOPPBqM4ENnGuclnS8kuYueoE26LTrrmPQScY1L4xo+9oTp/WHtf8llqdzmSf4euor/kx7kdKzaW45kuGHdQYdFRibYQDAYLVvf6YubV3097l4x/UtBBVt/tMOuOXHKSdSOBr69nYUsYrxsms0e6kqYstH40JoUsLtxovR3hKONO2T6PAWED3aI2nN2jk2cO7D+tJ9BKMaz+O57o8h5Xu1hhLVB+pqqdbTHnXWJF2RUN2VdlZ61kysRuvrI1k9REsg9l6+93QKm2/83Ky5fMJXdkdk8HWU6lkF5VhMktMmoZBpyOkeSMeDPXB06nmF57JKM7g4MWDrIpZxcGUgwB4Z0mGH9Dod0Ki0+C39oK1PXUkeQoEgntaDmFK0BTauNb+YL7bQd+2njwQ0ow1R+G+0ndYYP0x86w/pafpFK/lTmT04v18MCqY4TXQIeFyXb27svSepUzbPo0DARdJcxH8e6WZt74ysyBMx1d8xfH048y5aw7eDrfPNPwNgQqKm3T5YLvtmuUuwvc6VU9/Zm9t4INRwbzzYCeAm157ASwDou5q68ldbT2vv3M1Mmtm9l/cz47EHRxKOURcbpxlg5QExUvuC7csu2nUw44gSxVTmqtAJ3QM87uPx4Ieo5XLrTkXVn3y35FBuNhZsXQvjCl7mWcNq3hav45OujgmlL3IM98eJb/ExJjuzWu0HO3c2vH10K+ZvmM6xzjGzIl6Xlhp5t8rNTbcIfm6XwQPbXiIt/u8TV+fvjVaFqX6qKC4SYlZRdhSiqfIvak7isv9lYCoKwl5Caw9u5b1setJLUoFwLZUEpwk6ZAg6RIj8cmEHHv4oY9gS4iOXEeBlc6KUa0fYGLHifg6XT21t3JzrA06ZoV1oEcrd15YeYwPSkYRrrVjodV8frB+nb8ZZ/LSGsgtNvJEP/8aLYuHnQdfDPmC1/e+zoa4Dbw2Ts+4bRrDD0oCE83Muz+Hp7Y9xaQOk5gWOk1VRd0CVFDcpMSsInzEHz2eAHzd/ggKTWq3ZT17SmEK7x16j1/O/wKAXank7ihJ3xMabS+AXoJJB2eawcc9dey91IPJycqJSW0fYlz7cXja1+5dT0NyT0dvOjR1Zto3R9md2JmxZS+x1Po9Vlq/zriymfz3Zzh8PptH7/Dlrrae5JeY2HE6je3RaVzIKcbdwYbAJk6083YisIkzrTwcbqrnlI3ehtl9ZtPGtQ3zj8xnyRAzJ1toTP1J470lZhbdp+N//I8jaUd4v+/7NHFscv2DKnVGBcVNKCg1kV1kJLi8a6wXep2giYstmcWZzPx1JodSD9GmURve7P0m7dza1XGJ/zqTZuLrqK/5JOITikxF+Cdbxj30iJbYmCDRA9b3EES2EJxpJii1tny4+Ln4MTZgLMP9hze48Q91xdfNnm+n9OCpFUfYFg2jyl7jK+t3+N76Df5tnMKOqGC2RqXSyN6K/BITZk0CEnfyiMCJrVGp5cfy83Dg/uBmPBDSjObuN/bzE0IwqeMkOnt25vndz3MgII04b8H0tWZmrNHYHCpZdrelKuqt3m/Rv3n/an4nlOqiej3dhKiLedw7/1fG6zfzhtWXdCtZiK1bE3b+6y4m/jyREylHufOk5FAbgXUjN1bct6LCFdRuFVGZUbz626uczj6NY5Fk7E7L1BpF1paG6R1BOs42xTI3NtDWtS3dvLtxl89ddG/S/ba8s7oVGM0az/9wjLURyfiINJZZvUsrXQp50o6tWhe2mUNpLLK5QxdNN1007iKfAmnLGenDac2X47IVm8x3lI/p6drClX8Obksvf48bLkt2STYz98zktwu/oTdLxuy0VEWd84J59+u56C6Y2GEiz4Q+o6qiKqGWQq1GNR0UW06mMOWrw7xsWM7f9FsJLP0fvfw9GNYnjncOvsOIfRpjd2pEN4M3H9XTzN2P5fctvyXHBPx87mde3fMqJaZi7johGbddw6EENt4hWNlbR4mNQC/03NnsTob6D6WHd48bWhhIqVmaJpm14STL9p3HGiO9dJHcqzvEYH04rqIAgPOaFwe1AKJlc3xEOu1EIu10ibiLfIxSz06tM+vMvflF60KZsGbagDY8e3ebG+5irUmNJZFL+Pjox5ilmZCzGk/9qGFtgk/vsyyxGuoVypy75qjqyWtQQVGNajoovtgTz5s/nmKR1Tz8RTKDyt4nLNSOA8aZOKYX8sHnZi66gm+6ZazAh/fr6NrkDj4b+BlW+lvj25KUkkXHF7EwYiF2pZJp6zW6npVEN4PP79GT4CVws3VjQocJhPmH4WF3498yldohpWT5gQTmb40ho6AUAAMmOol4Lkq38oGYf3oVgSKB+/V7GKHfi7fIJlZrwpPGZzktLWNxPhwdfFOLJR1JPcILu18gtSgVtzzJ9LVmAi7Axq6C5QN0NHLw4P273q/R6eJvVSooqlFNB8Ws9SdZuvccG61nkipdmWx8ng4hP5BQfJiXvtNoe0Hy3BQ9vU9Jxm/XWN9dsHyAnjD/MN7q/Va9n4oiuySbdw68w6Zzm/DKlryw0kyzTFh2t46fuwoQOka1G8Uzoc/gbF31gYFK3TKZNXbHpLPycBK/nErFaLb83/dxtWNQ+8b09vcgo6CU6JR8Ii/kcjghGylBh0Y/XQTvWn2OE0W8YpzMKq0vjZ1tWDA6hO6tbnzp25ySHF7+7WV2J+1Gb5b8bYfG0EOWLyLzHtCT42T5HXs29FmcrK+ezqahUgPubiGXr2wXrrXF4HyMhJIj9IqSBMdLlgzSke0k+PEOaJwjCDsgSXPRWM96AtwCGNd+XB1fQcWS8pNYdmoZa2LWUGIuITBBMmO1GZ2E2aN1RLbU0ca1Da/3fJ1Onp3qurjKDTLodQwIaMyAgMbkFJVx4kIu3s62tPZyrPDLS3JOMeuPJbPqcBLb00IZWvo2860+Zq71IrqaTjMrbwKP/t9+pg1ow7QBrW9oqd9Gto34aMBHLD25lPlH5vPlQMGZZhpPbNR4/wszXw6UfCe/ZVvCNv59x78Z0mJIvf+CdTtTdxQ3YfC8XaSkpnDcdgoz5UjWNz+Fk7GAef9nJsMJXp6gx8XOlYKyAjSzkedXaYTESmZO0JPia8+aEWv+8hrO1el01mmWRC5h87nNmKUZzxzLSOqBEZLURvDfh/WkuAkGNh/I7D6zVe+lBsasST7aHsP8bTHopJnnDCt52rCOaM2XacZpxEgfurZw5cPRwTc1y/ChlEM8v+t5MksyaZYheWKjmbbJcLyl4PMhOlLcLNWcXRp3IdQrlJ5Ne+LfqGbHgtRXquqpGtVkUEgpaf/aZlqZzvKm8xv83as5Jn0Jj2+yLNE5c6KeeG/Bh/0/JKckh1n7ZmFfIpm32EyWE7w0QU8f3758cvcndfoNSUrJkbQjfHHiC3698CsAzdMkI/Zp9IqSSAG7Ogm+GqCjyFbwWKfHeDrkadWDqQH7NSadZ7+NIKuwjL66Y8y1+hQninnTNI4V5rtxtrVi3iPBV6z6V1XpRek8v/t5DqceRkjJwKOW3nUGE2zvLDjZQhDl+8eiVb2b9WZa8DQ6eHSo7sus1xpkUAgh7gHmA3rgcynlu3/aLi5tvw8oAiZKKY9c77g3GhQxGRcpKCup0r55xWVM/F84nZ1+5nzjo2CGKT9b5jFa112wYoCeQS0G8UG/DwB4a/9bfHf6O3qe0nhuncb/BurY1E3HnLvmMKTlkCqXsToYNSOHUw+zK3EXu5J2WRYLAtolSh7YpxEaKym2hl9CBBu7WarP7Ax2vNrjVYb7D6/Vsir108XcYp755iiHzmXjQS5zrT7lLv1xNpu78rzxcfJwYPrANjwzoA26G+wVZdJMfBP9DR8d/YhiUzGu+ZZedneckVibLPtccINdnXRs6moZq9Pftz9PBT9Vp2OVpJQUlJqwtzbcUE+w40k5fHcokfOZRfT0d2dczxY42167s4vJrNG7bx/MUuP/Vn5T6bFdbO1p5XZjgV0vg0IIoQfOAIOAJOAQ8KiU8tRl+9wHTMMSFN2B+VLK7tc79o0GRZ+lY8gVJ27sAgDnQsm/VpsJSLJMU7Gyj45Gtm6sHrG6vBdQQVkBI9aNIK0w9YqGbr2XJ+vuX1crjcFSSlbFrGLBkQVkl2YDYFMmCY6TDD2kEZAEeXbwUzcdm0MFhXYCe4M9D7d9mHHtx9HY4ca/ISq3L5NZ46PtZ/loewxSakzWb+Lfhm+5ID143PhPzkhfBgZ68cEjwZV+8F1LSmEKbx94mx2JOwDQmyWtLkL7RMv8YZ3OS3LsYU0vHVtDBEaDoEvjLjzS7hEGNh9YKz0LS4xmdp9JZ/PJVLZHp5JdZEQIcLGzws3eGk8nGzr7NiK0eSNCm7vi6WRDiVGjsMxE+LksvtgTz6Fz2YDEgBkTBlztrZg2oA1jezQnMauYg/FZHEnIJiGziMS8FHKdv+D8h5b3pNXMyudH8xTd2T7+8xu6pvoaFD2BWVLKIZcezwSQUr5z2T6fATullN9cenwa6CelvFjZsW82KHqf1NBrVXuNXoORv2k0KoRPhunYF6jDy96Ljwd8TKB74BX7bju/jek7p+OVLfngczNH/AUfPKjnkXaP8EqPV6pczptRZi7j7QNvsypmFQ7Fkl5RlnmYOp6XWJsh3Rk2dNexvbOgzMpSFzyu/ThGtRulejQplToYn8X0b4+SnFtCF3GaT63n40Axzxsf5yetB34eDnwyJpT2TW/u92hn4k6Wn1rO0bSjlGl/rA7ZNknyyG6NTuclGU6wuYvl9zff3vL7O6zVMIa0HEInj07VWr0rpeRIQjbL9p1ny8lUio1mQNJWJNFJxJOHPZnSmSycuCjdyxcHA9AJ0C77qLXCxH26/Uw2/EwHcY6vzXfzoWkkWThjpRflPdLA0uvMpclyTI1OkfafOAxmuP+RlhWW0WiAfYG62yooHgLukVL+49LjcUB3KeXTl+3zI/CulHLPpcfbgH9LKa9KASHEFGAKQPPmzbucP3++ymX5PSiWzTFha6z6NWQ5wnsP6YlrIujs2ZkP+39Y4XgCKSXP7HiGnYk7uX+vxphdGu88rCOitZ75/efX2NQFGcUZPLfjOSLSI2iRKnl+lRmvXEhpBIfbCMJbC6J9BWa9wMfRh0kdJxHmH4at4eoFkhSlIjlFZTz3XQQ7TqfjRTafWn9IF10Mi0zDed80Cr3BilnDO/DoHb43/aFdZi4jMiOSw6mHWR2zmqSCJAA6nLOsptfxvGV24n0Bgl9CdZy+tABWU4emDGoxiLtb3E2QRxB63c1NvllcZmbD8WSW7TtH5IU8QNJZxHKv/hCDdYdopUu56jUmqSNKNidCa81RrTV5OOBAMQ6iFG+RySj9rvKxKcekP2G6vRRhw0LTCDZpdxAizl4aMX+aFiKFAS28sc3TkfeSZXbmL5u3qLCsOQ4w5RnDbRUUDwND/hQUd0gpp122z0bgnT8FxQtSysOVHftm7yg8cyTiBt6ObEcwWgn6NL6H+YNmY62/9kpiKYUphK0No6y0iP8uMeNUDK+M15PrZsPCgQvp3uS6NWpVUmou5XDqYfZe2MtP8T+RXpxOjyiNJzdqFNrCh/fry/8j6YSOYM9gRgeMZlCLQRh0qre0cuM0Tf5/e3ceX1V97nv886y9M4dMEMYwT0GGBGWeFUVKnbDqkaJFPA7nqtXWaiv19N6rnfvyWrVarUexOAA9WgfUyqwiCjJjGEIIEBJIGAKEkDnZ67l/rM0kEIgnyd62z/v12q89sPbaX/aQZ63fb63fj6cW5/DM0lwiqeX/+F9lqn8Jq9ze3F9zH/toydUZ7fnN5H60aEBT1Lr8I7y8fBcoTB7YgfF9WlOndbyX+x4vbHzhxKjFHYqVCetcxm7yZkw8mAAr+ggr0h12tAPE29MYmzaWcR3HMbz9cGL89U8JoKqs2X2Et9bs4cOsIsqq64immmt9XzDNt5CLnN3Uqo8V7kUscAez0u1DLNW0lFJSKKWrs4+Bsp0MZyctpPKM9S8L9Gdm4Dt86g5Acegue5nhn83lvvUnlinVWFa7vVnjtOX1bpsZmu3y1dO7yE8V+v6469k/CweKE+WfqlCETdPTla/fQ1F1w/soIqUFk7rcwGOXTr+graXXtrzGH1b/gbSDyuOvByiLhl/c6qMmKZaXJrzEgNQBDc6gquQcyWFl0UpWFK5g7f61VAW8jvmIOuX6z72truwO8P+u93E0XhiXNo4JXSYwqsMokqOTG/yaxpzNgs37+Ml/b6Ssuo5rnM/5TcTL1ODnwdr/xSfuQBJjIrh1WGemjehy3km0Psoq4odz1uO6AbwhCx0yOibx8ITejOzRkupANXO3vs2bOXPJL9sFeP1uQ7YpI7KVjJ2K34UDibC6p7C6l7f37DpCpBPF0HZDuazTpYxNG3tiuJDy6jq+3HWIZTnFLM0+QH7wfKneks+Nvk+50fcpiVLBVrcTrwau4MPAUEqJB7wZJOvcM/+WOrh0l0KiqKGCaMo0hjJiqMDba09v24IereP5MKsIVRgqW+np7GGt24tt2hEXB1/cNmI7vcKNywK8PyeP7I5Clxm9kED8Od+/LrGDmTflyQZ9fuFaKPx4ndnjgb14ndnfV9XNpyzzXeA+TnZmP6OqQ8637uaYCvWbqHPruGfxPawoWkHPvcov5gTYlwz/d6oPf0IiM6+ceUFHb6gqWw5v4b3c91iYt5BDVYeO/wMdDkHGjHmc1AAAFE5JREFUTiVjl3JRvjey65IM4eUJDkREMGPoDG7qfVMT/0/Nv6odB8u49411ZO87Rlcp4rmIZ7jI2c3sust4IXA1+dqGSL/DNRntGd2zFYO7pNA+6fSt+4+yirhvznpauKW8FvlbWkopf6qbzJuBsdThp1V8FKWVtdQEXEDxxewiIvlL/AmbEAkAEFepDN7uFY4Bu7z+uGPRsK6HsKansKHbyRGOpSYNf3Vfyo70pKa8PeDQiqN817eSG3yf0t/Jo0Z9LHQHM6tuAqu1NyAkxkRweZ82TOzXltE9vemFSypqOVxew9aiUtbuPsK6/CNsLSrFVW/OkLhIH/HRfgZ0SOL7QzsxontLRIStRaX8fn42n2zzRqSO8nszUw7p2pJi30e8X/ASD/09wJOf5pHbTnjw5Qf55chfNupnF5aFAk4c1fQU3uGxM1X11yLyHwCq+kLw8NhngYl4h8dOP1v/xNeFa6EAqKit4M6Fd/JV8VcM2OnyyJsu29vDr2/24Y+J46eDf8r1Pa8/sYeiqpRUl1BYXkhhWSG7ju5ift58th/ZDkBiuTJgpzIgT+mfp6R447yxNwU2dhPW9hCyugjJ0Sk8Oe5JBrU96/fAmEZTVRvgsfc3M2dVAVHU8Ih/DlN9i/HhMt8dzIt1V7FRe5xYvkNSDBd3TiYjLZEIn8PjH2whwT3KG5G/pbsUkq0dyXB2stttzVN13yNH02glpbTkKElSTgVRlGoshx0/xfGFlCfsoix2H8fbkaNqvA2nwTnKJblKfBXU+GBTF+/3sbaHcDjB+73F1vkZXl7HDZUHGFJZxTa3C28FxjIvMJwjeJ3ymR2TmDaiM5P6t7ugScfqAi4KRFzAmesFhys4VlVH99ZxJ9b9s2U/4x+7/sEzz9fxH7n57GkFL7z9QqOP8BC2haKphHOhADhafZTpC6az/ch2hm91eeBdlx3t4JlrfOxPEcakjWFcx3Gs3reaNfvWcLDy4GnPj6j1tpbGZHlDhjjqHd66qYuQ1Vn4qqtwMOlkU9ioDqP4xbBf0D6+aedMNuZU767fy8/fyaKiJkBrjnCbfwG3+BaTIBXku6ms1nRWuemsdnuzU9sB3nc2mVJmR/6GrlLEnbU/4TO3P5c6G3jI/9/0dS7sIJUjjsPS2Bg+iY1hRUw01Y73R9pxlfQC7/czaLvSpsRbflcbWNNDWN3LIa8NIIIEIqguzaS2ZCix6hWGqUM7k9GxeUdHnvzeZPbs386rTwb4Tmk+BxNg6cdLG61f8zgrFGGouLKYaR9NI/9YPkOzXe7+yMUfgL9e7h3qh8iJpqRuRUpqKbQqVVqWQq+9Slw1FCfAsn7Cyt4Ou9uABvdCkqOSGd5+OCM7jGRE+xE2sqsJmYLDFTy7NJd31u+lJuASRyWTfcsZ5Ww6MQcGQInG8ZXbjQ3anSuctXSVfdxR+xDL3f60iPJzrLoOwWWUs4loajikCRSTyFGNI4ZqEqSCBMqJkRqqNYIqIqnGuy5HqIrbS6BFLr74bMTv9T0c/30N2q5cnOvSey846vVrrOolfNHn5DwrvZPTuan3jc1+VGBNoIahbwyla0Etv34twKWV+RyLgT0b9jR6/6IVijBVWFbIHQvvoOBYASmlyr0feMeGr+sm1EZAeoGSWHFy+ZJYOJQA+a2FZX2FLZ0FFSHSiWRgm4EMbzec4e2Hk56SbkNtmLByoLSKV77I4/WVuzlWFTzNGqW7FDLIySFDcsl0dtJb8qkhgjtqf8Lnbn9G9mjJf/1gEB9nH+TpJTnk7C87sc7oCIfk2EiqagMcq6o7ozM5widE+33ERHqXaL+PyrpajgZ2UBWZhS9uK77ok7P5tajw9jKGbPOaciMCkNcaFg10WN5XqIzyDrn90SU/YmKXic0yBM+2w9u44f0buHy9y13zXYYH8pHISI5sPtLor2WFIoyV1ZTxu1W/470d7yGqTFqt3PypS0k8bO3ojW+T00E4kOgdinuqzNRMrutxHRO6TLChmM23Qnl1HZ9tL2Z13mHW5B1mU2FpcCpWTwxV+HApI5ZRPVrxXz8YREzkyX6Ag8eqqaipo1V8FLGRvtP68qpqXapqA0RFOET5ffUOp+G6SllNHQVH9/JF0XJWFH3GhuLV1LreiVQx1cqILcoV61267YeqCG9om7+P9MY+G5A6gIcHPUxm68wmeqc883bM49Hlj/LvCwKM3qyMqc4nMSqRwo2Fjf5aVii+BZbkL+HxFY9zuOow4ip6ypfc7/jp2KIj7ePb0yGuA50SOjEmbQxdE89+LLUx3xbl1XVk7T3KxoISNu4p4as9R6muc7l6QHsevrL3aUWiybPUlrN873KW5i/l0z2fUl5bDqp0L4KJa11Gb1LKYuBvYxyWZHqH2t7e73buH3j/Nz6Z73yeWP0Es7bM4rHXvL2wW/fk0y6uHblrcxv9taxQfEscqjzEsxue5cuiL0mNSWVQ20EMbjuYjNSM854gZIxpPBW1FSzIW8BbOW/xVfFXAHTZp0xbEqBvPuS3gmev9pHXVhjRfgR/GPOHJpnq+K6Fd7Gi8Av++scAn/UVHl23m25J3cha2fDzvs7HCoUxxnxD2YezeWbdM95w/KoMyVFuW+TSohJemOTweV+HtPg0nrr0qUYfxXbc38bh7CvmuecDvDjR4cX5efRr1Y8vl3/ZqK8D9RcK6/E0xph6pKek8+fL/8xfLv8LPZJ7sqq3w4zpPna2hQfmuUz9OMDe0gJu/ehWVu9b3WivW1xZzKGqQ3Q66G3M704VRCQkY7FZoTDGmAswosMI3rz6Te4feD+lcQ6Pf9/HgoHCtSuVR950qaus4IGlD7CjZEejvN7xk2o7H/DuF6RCjD8GofknPLNCYYwxF8jv+LlzwJ08N/45YmMSeHmijxcnOmTsVO77wKWsppR7Ft/DwYqD51/ZeeQcyQGg8wFlXxJURXnzxISCFQpjjGmg0WmjmfvdufRI6sHigQ5vXOYwPFuZ8olLYXkh9y65l4raivOvqB7HC0WnA0p+a28vIlQHtVihMMaYb6BTQideufIVOid05v0hwsKBwnUrlfEbXLYe3srDyx4m4Aa+8fpzjuQQWau0OwK7W3uP2R6FMcZ8yyRFJ/H8+OdJjk5h5gSH9d2EO+a7ZOx0WbZnGXOy65/b+lxq3Vp2lOwgrdgbVmT38T2KCNujMMaYb52OCR350/g/ERERzR+vcyhIhR/Oc4mrVP684c8cqjzU4HXuPrqbWreWzgdOHvGUEp1ChNP084GfjRUKY4z5H8pIzeD3o39PdZTDc1f5iK+Cmz5zOVZ7jKfXPd3g9Z3aP1EVAQeSoVdyr8aOfcGsUBhjTCMY33k81/e8nt1thEUDhQnrlI4HlHdy3yHrYMPOpM4q9pbvfBDyU72Roa1QGGPMP4H7L76fFpEt+Ntoh4oomL7YBVV+u+q3uOpe0Doq6yqZt2Me8RVKr73KjnZe/0Rjn/XdEFYojDGmkaREp3Bv5r2UxQpzxzr0260My1ayirOYt2PeBa3jg50fUFpTyviNSmQdLM50iHAiGNVhVBOnPzcrFMYY04j+rfe/eedXZAp5reEHS10ia5U/rv0jZTVl9T5XVZm9dTa+gHLlWpeszkJBa2FS10mkRKc00//gTFYojDGmEfkdPzOGzEAdYeYEH61K4dqVLoerDvPqllfrfe7KopXkluQyOEdpdQz+MdhrdrrloluaI/o5WaEwxphGNqTdECZ0nkB2R2FFunD1l0pCuTJr8ywOVx0+5/Pe2PoGAJPWuOxPgnXdhUvaXEJ6SnpzRT8rKxTGGNMEfnzJj/E7fuaOcYiog8lfuFTUVfBy1stnXT6/NJ9le5bRtUhJ3wPzL3FQR7ilT2j3JsAKhTHGNIm0Fmnc0PMGiloKnwwQJqxXUkuUudlz2Ve+74zlZ2fPRlEmrXGpjISlA7w5ui/teGkI0p/OCoUxxjSRuzPuJsYfw5ujHBS4cblLjVvDCxtfOG25spoy3s19l8QyZeQW5ZP+QmW0MCV9SpNNs9oQViiMMaaJtIppxdQ+UzmcICy4RBizSUk7qLyb+y55R/NOLPfa1tcory1nwjoXv+s1O8X4Y5jcc3Lowp/CCoUxxjSh2/reRovIFrwz3KEqEm5e5hLQAM+sfwaAkqoSZm2eRVyl8t01yuqeQlFL4epuVzfJPNzfhBUKY4xpQolRidze73bKYoV5Qx2G5CjpBcqi3YuYt2MeMzfNpLy2nOtWukRXw9yxDn7Hz/R+00Md/QQrFMYY08Sm9plKq5hWfDhYOJAI970fILZK+dXKXzE7ezbJx5SJa5TlfYWCVOF7Pb9HWou0UMc+wQqFMcY0sRh/DA9e8iDVkcLT1/poWQp3f+RSWVtBdaCa6z938bnw5miHaF80dw+4O9SRT2OFwhhjmsFV3a7iqm5Xsb2DNw7U8Gzlso1K6yPK+I3Kkgxhf7Iwpc8UUmNTQx33NP5QBzDGmH8FIsJ/DvtPb4DAYXn0zxOmL3LJ6SAEHPj7KIf4iHhu73t7qKOewfYojDGmmcRFxPHE2Cfw+yJ59mrvKKj+u5X5g4SSeOG2vreRFJ0U6phnCEmhEJEUEVkkItuD18nnWC5PRLJEZIOIrGnunMYY09jSU9J5aNBDlMQLT13rsK678O4wh5TolJAP/ncuodqjeARYoqo9gSXB++dyqapmquqg5olmjDFNa0r6FK7pfg2buzj87iYf5THCz4f+nLiIuFBHO6tQ9VFcC4wL3p4FfAL8LERZjDGmWYkIvxr5K0anjWZnyU5GtB9BZuvMUMc6p1AVijaqWgSgqkUi0vocyymwUEQU+IuqvniuFYrIXcBdAJ06dWrsvMYY06hEhIldJjboOZmZoSkmoqpNs2KRxUDbs/zTo8AsVU06ZdkjqnpGP4WItFfVwmAhWQT8UFWXne+1Bw0apGvWWJeGMcZcKBFZe64m/ibbo1DVy+sJtF9E2gX3JtoBB86xjsLg9QEReQcYApy3UBhjjGk8oerMngdMC96eBrz39QVEJE5EWhy/DUwANjVbQmOMMUDoCsXvgCtEZDtwRfA+ItJeRP4RXKYNsFxENgKrgA9VdX5I0hpjzL+wkHRmq+ohYPxZHi8EJgVv7wQymjmaMcaYr7Ezs40xxtTLCoUxxph6WaEwxhhTLysUxhhj6tVkJ9yFkogcBHZ/w6e3AoobMU5jCMdMYLkaIhwzQXjmCsdMEJ65GjNTZ1U960QY/5SF4n9CRNaE2wCE4ZgJLFdDhGMmCM9c4ZgJwjNXc2WypidjjDH1skJhjDGmXlYoznTOEWpDKBwzgeVqiHDMBOGZKxwzQXjmapZM1kdhjDGmXrZHYYwxpl5WKIwxxtTLCkWQiEwUkW0ikisi9c3h3dQ5ZorIARHZdMpjKSKySES2B6/PmOSpiTN1FJGPRWSriGwWkQfCJFe0iKwSkY3BXI+FQ65gBp+IrBeRD8IoU56IZInIBhFZE0a5kkTkLRHJDn7Hhocyl4j0Dr5Hxy+lIvKjMHmvfhz8rm8SkTnB30CT57JCgfejBp4DvgNcBEwRkYtCFOevwNfnR3wEWKKqPYElwfvNqQ74iar2AYYB9wbfn1DnqgYuU9UMIBOYKCLDwiAXwAPA1lPuh0MmgEtVNfOUY+/DIdfTwHxVTccbMXprKHOp6rbge5QJXAJUAO+EMhOAiHQA7gcGqWo/wAfc3Cy5VPVf/gIMBxaccn8GMCOEeboAm065vw1oF7zdDtgW4vfrPbx5RMImFxALrAOGhjoXkBb8wV4GfBAunyGQB7T62mOhfq8SgF0ED6wJl1yn5JgAfB4OmYAOQAGQgjdFxAfBfE2ey/YoPMc/gOP2BB8LF21UtQggeN06VEFEpAswEPgyHHIFm3g24E2nu0hVwyHXU8BPAfeUx0KdCUCBhSKyVkTuCpNc3YCDwCvBprqXgjNahjrXcTcDc4K3Q5pJVfcCTwD5QBFwVFUXNkcuKxQeOctjdtzw14hIPPB34EeqWhrqPACqGlCviSANGCIi/UKZR0SuAg6o6tpQ5jiHkap6MV4T670iMibUgfC2jC8GnlfVgUA5oWuWO42IRALXAG+GOgtAsO/hWqAr0B6IE5FbmuO1rVB49gAdT7mfBhSGKMvZ7BeRdgDB6wPNHUBEIvCKxBuq+na45DpOVUuAT/D6d0KZayRwjYjkAXOBy0Tk9RBnAk7MIImqHsBrcx8SBrn2AHuCe4IAb+EVjlDnAq+grlPV/cH7oc50ObBLVQ+qai3wNjCiOXJZofCsBnqKSNfgVsTNwLwQZzrVPGBa8PY0vD6CZiMiArwMbFXVJ8MoV6qIJAVvx+D9kLJDmUtVZ6hqmqp2wfseLVXVW0KZCUBE4kSkxfHbeG3bm0KdS1X3AQUi0jv40HhgS6hzBU3hZLMThD5TPjBMRGKDv8nxeB3/TZ8rFB1E4XjBm6s7B9gBPBrCHHPw2h9r8ba2/h1oidc5uj14ndLMmUbhNcV9BWwIXiaFQa4BwPpgrk3A/w4+HtJcp+Qbx8nO7FC/V92AjcHL5uPf8VDnCmbIBNYEP8d3geRQ58I7OOIQkHjKY+HwXj2GtzG0CXgNiGqOXDaEhzHGmHpZ05Mxxph6WaEwxhhTLysUxhhj6mWFwhhjTL2sUBhjjKmXFQpjjDH1skJhjDGmXlYojGliIvLL43N4BO//WkTuD2UmYxrCTrgzpokFR9x9W1UvFhEH7wzaIap6KKTBjLlA/lAHMOafnarmicghERkItAHWW5Ew3yZWKIxpHi8BtwFtgZmhjWJMw1jTkzHNIDgqcRYQAfRU1UCIIxlzwWyPwphmoKo1IvIxUGJFwnzbWKEwphkEO7GHATeGOosxDWWHxxrTxETkIiAXWKKq20Odx5iGsj4KY4wx9bI9CmOMMfWyQmGMMaZeViiMMcbUywqFMcaYelmhMMYYU6//D6BEO18ahHmiAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "ib, ie = 385, 467\n",
    "vveln19 = restart_file.vn[0, pz]\n",
    "vvelb19 = restart_file.vb[0, pz]\n",
    "vvelb19[ib:ie, 1].plot(label=\"Before 1\", linewidth=5)\n",
    "vveln19[ib:ie, 1].plot(label=\"Next 1\")\n",
    "vvelb19[ib:ie, 3].plot(label=\"Before 3\", linewidth=5)\n",
    "vveln19[ib:ie, 3].plot(label=\"Next 3\")\n",
    "plt.plot([prob-ib, prob-ib], [-0.5, 2], color='k');\n",
    "plt.legend();"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [],
   "source": [
    "restart_file.un[0, :, :, 1] = restart_file.un[0, :, :, 2]\n",
    "restart_file.ub[0, :, :, 1] = restart_file.ub[0, :, :, 2]\n",
    "restart_file.vn[0, :, :, 1] = restart_file.vn[0, :, :, 2]\n",
    "restart_file.vb[0, :, :, 1] = restart_file.vb[0, :, :, 2]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "{'file_name': 'SalishSea_04914000_restart.nc',\n",
       " 'TimeStamp': '13/11/2020 10:56:10 -0800',\n",
       " 'Modified': 'copied velocity one away from west boundary to boundary'}"
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "restart_file.attrs['Modified'] = 'copied velocity one away from west boundary to boundary'\n",
    "restart_file.attrs"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [],
   "source": [
    "restart_file.to_netcdf('SalishSea_04914000_restart.nc')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "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.7.6"
  },
  "toc": {
   "base_numbering": 1,
   "nav_menu": {},
   "number_sections": true,
   "sideBar": true,
   "skip_h1_title": false,
   "title_cell": "Table of Contents",
   "title_sidebar": "Contents",
   "toc_cell": true,
   "toc_position": {},
   "toc_section_display": true,
   "toc_window_display": false
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}