Budget_closure.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 54,
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "import xarray as xr\n",
    "import matplotlib.pyplot as plt"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [],
   "source": [
    "yr2sec = 24*3600*365.25"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [],
   "source": [
    "surface_hfds = xr.open_mfdataset('~/UNSW_work/ACCESS_CM2/hfds_Omon_ACCESS-CM2_historical_r1i1p1f1_gn_185001-201412.nc')\n",
    "surface_area = xr.open_mfdataset('~/UNSW_work/ACCESS_CM2/area*')\n",
    "OHC = xr.open_mfdataset('~/UNSW_work/ACCESS_CM2/OHC_ACCESS_CM2_CMIP.nc')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 52,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div><svg style=\"position: absolute; width: 0; height: 0; overflow: hidden\">\n",
       "<defs>\n",
       "<symbol id=\"icon-database\" viewBox=\"0 0 32 32\">\n",
       "<path d=\"M16 0c-8.837 0-16 2.239-16 5v4c0 2.761 7.163 5 16 5s16-2.239 16-5v-4c0-2.761-7.163-5-16-5z\"></path>\n",
       "<path d=\"M16 17c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n",
       "<path d=\"M16 26c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n",
       "</symbol>\n",
       "<symbol id=\"icon-file-text2\" viewBox=\"0 0 32 32\">\n",
       "<path d=\"M28.681 7.159c-0.694-0.947-1.662-2.053-2.724-3.116s-2.169-2.030-3.116-2.724c-1.612-1.182-2.393-1.319-2.841-1.319h-15.5c-1.378 0-2.5 1.121-2.5 2.5v27c0 1.378 1.122 2.5 2.5 2.5h23c1.378 0 2.5-1.122 2.5-2.5v-19.5c0-0.448-0.137-1.23-1.319-2.841zM24.543 5.457c0.959 0.959 1.712 1.825 2.268 2.543h-4.811v-4.811c0.718 0.556 1.584 1.309 2.543 2.268zM28 29.5c0 0.271-0.229 0.5-0.5 0.5h-23c-0.271 0-0.5-0.229-0.5-0.5v-27c0-0.271 0.229-0.5 0.5-0.5 0 0 15.499-0 15.5 0v7c0 0.552 0.448 1 1 1h7v19.5z\"></path>\n",
       "<path d=\"M23 26h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "<path d=\"M23 22h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "<path d=\"M23 18h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "</symbol>\n",
       "</defs>\n",
       "</svg>\n",
       "<style>/* CSS stylesheet for displaying xarray objects in jupyterlab.\n",
       " *\n",
       " */\n",
       "\n",
       ":root {\n",
       "  --xr-font-color0: var(--jp-content-font-color0, rgba(0, 0, 0, 1));\n",
       "  --xr-font-color2: var(--jp-content-font-color2, rgba(0, 0, 0, 0.54));\n",
       "  --xr-font-color3: var(--jp-content-font-color3, rgba(0, 0, 0, 0.38));\n",
       "  --xr-border-color: var(--jp-border-color2, #e0e0e0);\n",
       "  --xr-disabled-color: var(--jp-layout-color3, #bdbdbd);\n",
       "  --xr-background-color: var(--jp-layout-color0, white);\n",
       "  --xr-background-color-row-even: var(--jp-layout-color1, white);\n",
       "  --xr-background-color-row-odd: var(--jp-layout-color2, #eeeeee);\n",
       "}\n",
       "\n",
       "html[theme=dark],\n",
       "body.vscode-dark {\n",
       "  --xr-font-color0: rgba(255, 255, 255, 1);\n",
       "  --xr-font-color2: rgba(255, 255, 255, 0.54);\n",
       "  --xr-font-color3: rgba(255, 255, 255, 0.38);\n",
       "  --xr-border-color: #1F1F1F;\n",
       "  --xr-disabled-color: #515151;\n",
       "  --xr-background-color: #111111;\n",
       "  --xr-background-color-row-even: #111111;\n",
       "  --xr-background-color-row-odd: #313131;\n",
       "}\n",
       "\n",
       ".xr-wrap {\n",
       "  display: block;\n",
       "  min-width: 300px;\n",
       "  max-width: 700px;\n",
       "}\n",
       "\n",
       ".xr-text-repr-fallback {\n",
       "  /* fallback to plain text repr when CSS is not injected (untrusted notebook) */\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-header {\n",
       "  padding-top: 6px;\n",
       "  padding-bottom: 6px;\n",
       "  margin-bottom: 4px;\n",
       "  border-bottom: solid 1px var(--xr-border-color);\n",
       "}\n",
       "\n",
       ".xr-header > div,\n",
       ".xr-header > ul {\n",
       "  display: inline;\n",
       "  margin-top: 0;\n",
       "  margin-bottom: 0;\n",
       "}\n",
       "\n",
       ".xr-obj-type,\n",
       ".xr-array-name {\n",
       "  margin-left: 2px;\n",
       "  margin-right: 10px;\n",
       "}\n",
       "\n",
       ".xr-obj-type {\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-sections {\n",
       "  padding-left: 0 !important;\n",
       "  display: grid;\n",
       "  grid-template-columns: 150px auto auto 1fr 20px 20px;\n",
       "}\n",
       "\n",
       ".xr-section-item {\n",
       "  display: contents;\n",
       "}\n",
       "\n",
       ".xr-section-item input {\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-section-item input + label {\n",
       "  color: var(--xr-disabled-color);\n",
       "}\n",
       "\n",
       ".xr-section-item input:enabled + label {\n",
       "  cursor: pointer;\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-section-item input:enabled + label:hover {\n",
       "  color: var(--xr-font-color0);\n",
       "}\n",
       "\n",
       ".xr-section-summary {\n",
       "  grid-column: 1;\n",
       "  color: var(--xr-font-color2);\n",
       "  font-weight: 500;\n",
       "}\n",
       "\n",
       ".xr-section-summary > span {\n",
       "  display: inline-block;\n",
       "  padding-left: 0.5em;\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:disabled + label {\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-section-summary-in + label:before {\n",
       "  display: inline-block;\n",
       "  content: '►';\n",
       "  font-size: 11px;\n",
       "  width: 15px;\n",
       "  text-align: center;\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:disabled + label:before {\n",
       "  color: var(--xr-disabled-color);\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:checked + label:before {\n",
       "  content: '▼';\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:checked + label > span {\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-section-summary,\n",
       ".xr-section-inline-details {\n",
       "  padding-top: 4px;\n",
       "  padding-bottom: 4px;\n",
       "}\n",
       "\n",
       ".xr-section-inline-details {\n",
       "  grid-column: 2 / -1;\n",
       "}\n",
       "\n",
       ".xr-section-details {\n",
       "  display: none;\n",
       "  grid-column: 1 / -1;\n",
       "  margin-bottom: 5px;\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:checked ~ .xr-section-details {\n",
       "  display: contents;\n",
       "}\n",
       "\n",
       ".xr-array-wrap {\n",
       "  grid-column: 1 / -1;\n",
       "  display: grid;\n",
       "  grid-template-columns: 20px auto;\n",
       "}\n",
       "\n",
       ".xr-array-wrap > label {\n",
       "  grid-column: 1;\n",
       "  vertical-align: top;\n",
       "}\n",
       "\n",
       ".xr-preview {\n",
       "  color: var(--xr-font-color3);\n",
       "}\n",
       "\n",
       ".xr-array-preview,\n",
       ".xr-array-data {\n",
       "  padding: 0 5px !important;\n",
       "  grid-column: 2;\n",
       "}\n",
       "\n",
       ".xr-array-data,\n",
       ".xr-array-in:checked ~ .xr-array-preview {\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-array-in:checked ~ .xr-array-data,\n",
       ".xr-array-preview {\n",
       "  display: inline-block;\n",
       "}\n",
       "\n",
       ".xr-dim-list {\n",
       "  display: inline-block !important;\n",
       "  list-style: none;\n",
       "  padding: 0 !important;\n",
       "  margin: 0;\n",
       "}\n",
       "\n",
       ".xr-dim-list li {\n",
       "  display: inline-block;\n",
       "  padding: 0;\n",
       "  margin: 0;\n",
       "}\n",
       "\n",
       ".xr-dim-list:before {\n",
       "  content: '(';\n",
       "}\n",
       "\n",
       ".xr-dim-list:after {\n",
       "  content: ')';\n",
       "}\n",
       "\n",
       ".xr-dim-list li:not(:last-child):after {\n",
       "  content: ',';\n",
       "  padding-right: 5px;\n",
       "}\n",
       "\n",
       ".xr-has-index {\n",
       "  font-weight: bold;\n",
       "}\n",
       "\n",
       ".xr-var-list,\n",
       ".xr-var-item {\n",
       "  display: contents;\n",
       "}\n",
       "\n",
       ".xr-var-item > div,\n",
       ".xr-var-item label,\n",
       ".xr-var-item > .xr-var-name span {\n",
       "  background-color: var(--xr-background-color-row-even);\n",
       "  margin-bottom: 0;\n",
       "}\n",
       "\n",
       ".xr-var-item > .xr-var-name:hover span {\n",
       "  padding-right: 5px;\n",
       "}\n",
       "\n",
       ".xr-var-list > li:nth-child(odd) > div,\n",
       ".xr-var-list > li:nth-child(odd) > label,\n",
       ".xr-var-list > li:nth-child(odd) > .xr-var-name span {\n",
       "  background-color: var(--xr-background-color-row-odd);\n",
       "}\n",
       "\n",
       ".xr-var-name {\n",
       "  grid-column: 1;\n",
       "}\n",
       "\n",
       ".xr-var-dims {\n",
       "  grid-column: 2;\n",
       "}\n",
       "\n",
       ".xr-var-dtype {\n",
       "  grid-column: 3;\n",
       "  text-align: right;\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-var-preview {\n",
       "  grid-column: 4;\n",
       "}\n",
       "\n",
       ".xr-var-name,\n",
       ".xr-var-dims,\n",
       ".xr-var-dtype,\n",
       ".xr-preview,\n",
       ".xr-attrs dt {\n",
       "  white-space: nowrap;\n",
       "  overflow: hidden;\n",
       "  text-overflow: ellipsis;\n",
       "  padding-right: 10px;\n",
       "}\n",
       "\n",
       ".xr-var-name:hover,\n",
       ".xr-var-dims:hover,\n",
       ".xr-var-dtype:hover,\n",
       ".xr-attrs dt:hover {\n",
       "  overflow: visible;\n",
       "  width: auto;\n",
       "  z-index: 1;\n",
       "}\n",
       "\n",
       ".xr-var-attrs,\n",
       ".xr-var-data {\n",
       "  display: none;\n",
       "  background-color: var(--xr-background-color) !important;\n",
       "  padding-bottom: 5px !important;\n",
       "}\n",
       "\n",
       ".xr-var-attrs-in:checked ~ .xr-var-attrs,\n",
       ".xr-var-data-in:checked ~ .xr-var-data {\n",
       "  display: block;\n",
       "}\n",
       "\n",
       ".xr-var-data > table {\n",
       "  float: right;\n",
       "}\n",
       "\n",
       ".xr-var-name span,\n",
       ".xr-var-data,\n",
       ".xr-attrs {\n",
       "  padding-left: 25px !important;\n",
       "}\n",
       "\n",
       ".xr-attrs,\n",
       ".xr-var-attrs,\n",
       ".xr-var-data {\n",
       "  grid-column: 1 / -1;\n",
       "}\n",
       "\n",
       "dl.xr-attrs {\n",
       "  padding: 0;\n",
       "  margin: 0;\n",
       "  display: grid;\n",
       "  grid-template-columns: 125px auto;\n",
       "}\n",
       "\n",
       ".xr-attrs dt,\n",
       ".xr-attrs dd {\n",
       "  padding: 0;\n",
       "  margin: 0;\n",
       "  float: left;\n",
       "  padding-right: 10px;\n",
       "  width: auto;\n",
       "}\n",
       "\n",
       ".xr-attrs dt {\n",
       "  font-weight: normal;\n",
       "  grid-column: 1;\n",
       "}\n",
       "\n",
       ".xr-attrs dt:hover span {\n",
       "  display: inline-block;\n",
       "  background: var(--xr-background-color);\n",
       "  padding-right: 10px;\n",
       "}\n",
       "\n",
       ".xr-attrs dd {\n",
       "  grid-column: 2;\n",
       "  white-space: pre-wrap;\n",
       "  word-break: break-all;\n",
       "}\n",
       "\n",
       ".xr-icon-database,\n",
       ".xr-icon-file-text2 {\n",
       "  display: inline-block;\n",
       "  vertical-align: middle;\n",
       "  width: 1em;\n",
       "  height: 1.5em !important;\n",
       "  stroke-width: 0;\n",
       "  stroke: currentColor;\n",
       "  fill: currentColor;\n",
       "}\n",
       "</style><pre class='xr-text-repr-fallback'>&lt;xarray.DataArray (time: 1980)&gt;\n",
       "array([1.4319747e+22, 2.6224199e+22, 3.0627833e+22, ..., 3.2706899e+23,\n",
       "       3.4030244e+23, 3.5309187e+23], dtype=float32)\n",
       "Coordinates:\n",
       "  * time     (time) datetime64[ns] 1850-01-16T12:00:00 ... 2014-12-16T12:00:00</pre><div class='xr-wrap' hidden><div class='xr-header'><div class='xr-obj-type'>xarray.DataArray</div><div class='xr-array-name'></div><ul class='xr-dim-list'><li><span class='xr-has-index'>time</span>: 1980</li></ul></div><ul class='xr-sections'><li class='xr-section-item'><div class='xr-array-wrap'><input id='section-792f8248-4e5a-4d14-9af3-5b34a2a9913e' class='xr-array-in' type='checkbox' checked><label for='section-792f8248-4e5a-4d14-9af3-5b34a2a9913e' title='Show/hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-array-preview xr-preview'><span>1.432e+22 2.622e+22 3.063e+22 ... 3.271e+23 3.403e+23 3.531e+23</span></div><div class='xr-array-data'><pre>array([1.4319747e+22, 2.6224199e+22, 3.0627833e+22, ..., 3.2706899e+23,\n",
       "       3.4030244e+23, 3.5309187e+23], dtype=float32)</pre></div></div></li><li class='xr-section-item'><input id='section-f10e10c8-81ea-48a4-ab6f-a7cf4a716eab' class='xr-section-summary-in' type='checkbox'  checked><label for='section-f10e10c8-81ea-48a4-ab6f-a7cf4a716eab' class='xr-section-summary' >Coordinates: <span>(1)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><ul class='xr-var-list'><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>time</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>datetime64[ns]</div><div class='xr-var-preview xr-preview'>1850-01-16T12:00:00 ... 2014-12-...</div><input id='attrs-562df21d-303c-4fb9-8410-a28a6ef5fc7e' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-562df21d-303c-4fb9-8410-a28a6ef5fc7e' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-dfeace95-29df-4fd6-8511-22db847061a1' class='xr-var-data-in' type='checkbox'><label for='data-dfeace95-29df-4fd6-8511-22db847061a1' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>bounds :</span></dt><dd>time_bnds</dd><dt><span>axis :</span></dt><dd>T</dd><dt><span>long_name :</span></dt><dd>time</dd><dt><span>standard_name :</span></dt><dd>time</dd></dl></div><div class='xr-var-data'><pre>array([&#x27;1850-01-16T12:00:00.000000000&#x27;, &#x27;1850-02-15T00:00:00.000000000&#x27;,\n",
       "       &#x27;1850-03-16T12:00:00.000000000&#x27;, ..., &#x27;2014-10-16T12:00:00.000000000&#x27;,\n",
       "       &#x27;2014-11-16T00:00:00.000000000&#x27;, &#x27;2014-12-16T12:00:00.000000000&#x27;],\n",
       "      dtype=&#x27;datetime64[ns]&#x27;)</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-1f43acee-1276-4648-b61a-576a28cd4c63' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-1f43acee-1276-4648-b61a-576a28cd4c63' class='xr-section-summary'  title='Expand/collapse section'>Attributes: <span>(0)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><dl class='xr-attrs'></dl></div></li></ul></div></div>"
      ],
      "text/plain": [
       "<xarray.DataArray (time: 1980)>\n",
       "array([1.4319747e+22, 2.6224199e+22, 3.0627833e+22, ..., 3.2706899e+23,\n",
       "       3.4030244e+23, 3.5309187e+23], dtype=float32)\n",
       "Coordinates:\n",
       "  * time     (time) datetime64[ns] 1850-01-16T12:00:00 ... 2014-12-16T12:00:00"
      ]
     },
     "execution_count": 52,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "hfds_integrated = ((surface_hfds.hfds*surface_area.areacello)*(yr2sec/12)).cumsum('time').sum('j').sum('i')\n",
    "hfds_integrated.load()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 39,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div><svg style=\"position: absolute; width: 0; height: 0; overflow: hidden\">\n",
       "<defs>\n",
       "<symbol id=\"icon-database\" viewBox=\"0 0 32 32\">\n",
       "<path d=\"M16 0c-8.837 0-16 2.239-16 5v4c0 2.761 7.163 5 16 5s16-2.239 16-5v-4c0-2.761-7.163-5-16-5z\"></path>\n",
       "<path d=\"M16 17c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n",
       "<path d=\"M16 26c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n",
       "</symbol>\n",
       "<symbol id=\"icon-file-text2\" viewBox=\"0 0 32 32\">\n",
       "<path d=\"M28.681 7.159c-0.694-0.947-1.662-2.053-2.724-3.116s-2.169-2.030-3.116-2.724c-1.612-1.182-2.393-1.319-2.841-1.319h-15.5c-1.378 0-2.5 1.121-2.5 2.5v27c0 1.378 1.122 2.5 2.5 2.5h23c1.378 0 2.5-1.122 2.5-2.5v-19.5c0-0.448-0.137-1.23-1.319-2.841zM24.543 5.457c0.959 0.959 1.712 1.825 2.268 2.543h-4.811v-4.811c0.718 0.556 1.584 1.309 2.543 2.268zM28 29.5c0 0.271-0.229 0.5-0.5 0.5h-23c-0.271 0-0.5-0.229-0.5-0.5v-27c0-0.271 0.229-0.5 0.5-0.5 0 0 15.499-0 15.5 0v7c0 0.552 0.448 1 1 1h7v19.5z\"></path>\n",
       "<path d=\"M23 26h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "<path d=\"M23 22h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "<path d=\"M23 18h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "</symbol>\n",
       "</defs>\n",
       "</svg>\n",
       "<style>/* CSS stylesheet for displaying xarray objects in jupyterlab.\n",
       " *\n",
       " */\n",
       "\n",
       ":root {\n",
       "  --xr-font-color0: var(--jp-content-font-color0, rgba(0, 0, 0, 1));\n",
       "  --xr-font-color2: var(--jp-content-font-color2, rgba(0, 0, 0, 0.54));\n",
       "  --xr-font-color3: var(--jp-content-font-color3, rgba(0, 0, 0, 0.38));\n",
       "  --xr-border-color: var(--jp-border-color2, #e0e0e0);\n",
       "  --xr-disabled-color: var(--jp-layout-color3, #bdbdbd);\n",
       "  --xr-background-color: var(--jp-layout-color0, white);\n",
       "  --xr-background-color-row-even: var(--jp-layout-color1, white);\n",
       "  --xr-background-color-row-odd: var(--jp-layout-color2, #eeeeee);\n",
       "}\n",
       "\n",
       "html[theme=dark],\n",
       "body.vscode-dark {\n",
       "  --xr-font-color0: rgba(255, 255, 255, 1);\n",
       "  --xr-font-color2: rgba(255, 255, 255, 0.54);\n",
       "  --xr-font-color3: rgba(255, 255, 255, 0.38);\n",
       "  --xr-border-color: #1F1F1F;\n",
       "  --xr-disabled-color: #515151;\n",
       "  --xr-background-color: #111111;\n",
       "  --xr-background-color-row-even: #111111;\n",
       "  --xr-background-color-row-odd: #313131;\n",
       "}\n",
       "\n",
       ".xr-wrap {\n",
       "  display: block;\n",
       "  min-width: 300px;\n",
       "  max-width: 700px;\n",
       "}\n",
       "\n",
       ".xr-text-repr-fallback {\n",
       "  /* fallback to plain text repr when CSS is not injected (untrusted notebook) */\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-header {\n",
       "  padding-top: 6px;\n",
       "  padding-bottom: 6px;\n",
       "  margin-bottom: 4px;\n",
       "  border-bottom: solid 1px var(--xr-border-color);\n",
       "}\n",
       "\n",
       ".xr-header > div,\n",
       ".xr-header > ul {\n",
       "  display: inline;\n",
       "  margin-top: 0;\n",
       "  margin-bottom: 0;\n",
       "}\n",
       "\n",
       ".xr-obj-type,\n",
       ".xr-array-name {\n",
       "  margin-left: 2px;\n",
       "  margin-right: 10px;\n",
       "}\n",
       "\n",
       ".xr-obj-type {\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-sections {\n",
       "  padding-left: 0 !important;\n",
       "  display: grid;\n",
       "  grid-template-columns: 150px auto auto 1fr 20px 20px;\n",
       "}\n",
       "\n",
       ".xr-section-item {\n",
       "  display: contents;\n",
       "}\n",
       "\n",
       ".xr-section-item input {\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-section-item input + label {\n",
       "  color: var(--xr-disabled-color);\n",
       "}\n",
       "\n",
       ".xr-section-item input:enabled + label {\n",
       "  cursor: pointer;\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-section-item input:enabled + label:hover {\n",
       "  color: var(--xr-font-color0);\n",
       "}\n",
       "\n",
       ".xr-section-summary {\n",
       "  grid-column: 1;\n",
       "  color: var(--xr-font-color2);\n",
       "  font-weight: 500;\n",
       "}\n",
       "\n",
       ".xr-section-summary > span {\n",
       "  display: inline-block;\n",
       "  padding-left: 0.5em;\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:disabled + label {\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-section-summary-in + label:before {\n",
       "  display: inline-block;\n",
       "  content: '►';\n",
       "  font-size: 11px;\n",
       "  width: 15px;\n",
       "  text-align: center;\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:disabled + label:before {\n",
       "  color: var(--xr-disabled-color);\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:checked + label:before {\n",
       "  content: '▼';\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:checked + label > span {\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-section-summary,\n",
       ".xr-section-inline-details {\n",
       "  padding-top: 4px;\n",
       "  padding-bottom: 4px;\n",
       "}\n",
       "\n",
       ".xr-section-inline-details {\n",
       "  grid-column: 2 / -1;\n",
       "}\n",
       "\n",
       ".xr-section-details {\n",
       "  display: none;\n",
       "  grid-column: 1 / -1;\n",
       "  margin-bottom: 5px;\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:checked ~ .xr-section-details {\n",
       "  display: contents;\n",
       "}\n",
       "\n",
       ".xr-array-wrap {\n",
       "  grid-column: 1 / -1;\n",
       "  display: grid;\n",
       "  grid-template-columns: 20px auto;\n",
       "}\n",
       "\n",
       ".xr-array-wrap > label {\n",
       "  grid-column: 1;\n",
       "  vertical-align: top;\n",
       "}\n",
       "\n",
       ".xr-preview {\n",
       "  color: var(--xr-font-color3);\n",
       "}\n",
       "\n",
       ".xr-array-preview,\n",
       ".xr-array-data {\n",
       "  padding: 0 5px !important;\n",
       "  grid-column: 2;\n",
       "}\n",
       "\n",
       ".xr-array-data,\n",
       ".xr-array-in:checked ~ .xr-array-preview {\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-array-in:checked ~ .xr-array-data,\n",
       ".xr-array-preview {\n",
       "  display: inline-block;\n",
       "}\n",
       "\n",
       ".xr-dim-list {\n",
       "  display: inline-block !important;\n",
       "  list-style: none;\n",
       "  padding: 0 !important;\n",
       "  margin: 0;\n",
       "}\n",
       "\n",
       ".xr-dim-list li {\n",
       "  display: inline-block;\n",
       "  padding: 0;\n",
       "  margin: 0;\n",
       "}\n",
       "\n",
       ".xr-dim-list:before {\n",
       "  content: '(';\n",
       "}\n",
       "\n",
       ".xr-dim-list:after {\n",
       "  content: ')';\n",
       "}\n",
       "\n",
       ".xr-dim-list li:not(:last-child):after {\n",
       "  content: ',';\n",
       "  padding-right: 5px;\n",
       "}\n",
       "\n",
       ".xr-has-index {\n",
       "  font-weight: bold;\n",
       "}\n",
       "\n",
       ".xr-var-list,\n",
       ".xr-var-item {\n",
       "  display: contents;\n",
       "}\n",
       "\n",
       ".xr-var-item > div,\n",
       ".xr-var-item label,\n",
       ".xr-var-item > .xr-var-name span {\n",
       "  background-color: var(--xr-background-color-row-even);\n",
       "  margin-bottom: 0;\n",
       "}\n",
       "\n",
       ".xr-var-item > .xr-var-name:hover span {\n",
       "  padding-right: 5px;\n",
       "}\n",
       "\n",
       ".xr-var-list > li:nth-child(odd) > div,\n",
       ".xr-var-list > li:nth-child(odd) > label,\n",
       ".xr-var-list > li:nth-child(odd) > .xr-var-name span {\n",
       "  background-color: var(--xr-background-color-row-odd);\n",
       "}\n",
       "\n",
       ".xr-var-name {\n",
       "  grid-column: 1;\n",
       "}\n",
       "\n",
       ".xr-var-dims {\n",
       "  grid-column: 2;\n",
       "}\n",
       "\n",
       ".xr-var-dtype {\n",
       "  grid-column: 3;\n",
       "  text-align: right;\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-var-preview {\n",
       "  grid-column: 4;\n",
       "}\n",
       "\n",
       ".xr-var-name,\n",
       ".xr-var-dims,\n",
       ".xr-var-dtype,\n",
       ".xr-preview,\n",
       ".xr-attrs dt {\n",
       "  white-space: nowrap;\n",
       "  overflow: hidden;\n",
       "  text-overflow: ellipsis;\n",
       "  padding-right: 10px;\n",
       "}\n",
       "\n",
       ".xr-var-name:hover,\n",
       ".xr-var-dims:hover,\n",
       ".xr-var-dtype:hover,\n",
       ".xr-attrs dt:hover {\n",
       "  overflow: visible;\n",
       "  width: auto;\n",
       "  z-index: 1;\n",
       "}\n",
       "\n",
       ".xr-var-attrs,\n",
       ".xr-var-data {\n",
       "  display: none;\n",
       "  background-color: var(--xr-background-color) !important;\n",
       "  padding-bottom: 5px !important;\n",
       "}\n",
       "\n",
       ".xr-var-attrs-in:checked ~ .xr-var-attrs,\n",
       ".xr-var-data-in:checked ~ .xr-var-data {\n",
       "  display: block;\n",
       "}\n",
       "\n",
       ".xr-var-data > table {\n",
       "  float: right;\n",
       "}\n",
       "\n",
       ".xr-var-name span,\n",
       ".xr-var-data,\n",
       ".xr-attrs {\n",
       "  padding-left: 25px !important;\n",
       "}\n",
       "\n",
       ".xr-attrs,\n",
       ".xr-var-attrs,\n",
       ".xr-var-data {\n",
       "  grid-column: 1 / -1;\n",
       "}\n",
       "\n",
       "dl.xr-attrs {\n",
       "  padding: 0;\n",
       "  margin: 0;\n",
       "  display: grid;\n",
       "  grid-template-columns: 125px auto;\n",
       "}\n",
       "\n",
       ".xr-attrs dt,\n",
       ".xr-attrs dd {\n",
       "  padding: 0;\n",
       "  margin: 0;\n",
       "  float: left;\n",
       "  padding-right: 10px;\n",
       "  width: auto;\n",
       "}\n",
       "\n",
       ".xr-attrs dt {\n",
       "  font-weight: normal;\n",
       "  grid-column: 1;\n",
       "}\n",
       "\n",
       ".xr-attrs dt:hover span {\n",
       "  display: inline-block;\n",
       "  background: var(--xr-background-color);\n",
       "  padding-right: 10px;\n",
       "}\n",
       "\n",
       ".xr-attrs dd {\n",
       "  grid-column: 2;\n",
       "  white-space: pre-wrap;\n",
       "  word-break: break-all;\n",
       "}\n",
       "\n",
       ".xr-icon-database,\n",
       ".xr-icon-file-text2 {\n",
       "  display: inline-block;\n",
       "  vertical-align: middle;\n",
       "  width: 1em;\n",
       "  height: 1.5em !important;\n",
       "  stroke-width: 0;\n",
       "  stroke: currentColor;\n",
       "  fill: currentColor;\n",
       "}\n",
       "</style><pre class='xr-text-repr-fallback'>&lt;xarray.DataArray &#x27;__xarray_dataarray_variable__&#x27; (time: 1980)&gt;\n",
       "array([2.59977719e+25, 2.60113345e+25, 2.60189363e+25, ...,\n",
       "       2.66573444e+25, 2.66691887e+25, 2.66822825e+25])\n",
       "Coordinates:\n",
       "  * time             (time) datetime64[ns] 1850-01-16T12:00:00 ... 2014-12-16...\n",
       "    lev              float64 5.831e+03\n",
       "    ensemble_member  &lt;U8 &#x27;r1i1p1f1&#x27;</pre><div class='xr-wrap' hidden><div class='xr-header'><div class='xr-obj-type'>xarray.DataArray</div><div class='xr-array-name'>'__xarray_dataarray_variable__'</div><ul class='xr-dim-list'><li><span class='xr-has-index'>time</span>: 1980</li></ul></div><ul class='xr-sections'><li class='xr-section-item'><div class='xr-array-wrap'><input id='section-c8a2abda-ee97-4d92-8738-9e2664ccb1d2' class='xr-array-in' type='checkbox' checked><label for='section-c8a2abda-ee97-4d92-8738-9e2664ccb1d2' title='Show/hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-array-preview xr-preview'><span>2.6e+25 2.601e+25 2.602e+25 ... 2.666e+25 2.667e+25 2.668e+25</span></div><div class='xr-array-data'><pre>array([2.59977719e+25, 2.60113345e+25, 2.60189363e+25, ...,\n",
       "       2.66573444e+25, 2.66691887e+25, 2.66822825e+25])</pre></div></div></li><li class='xr-section-item'><input id='section-8dc03a84-0789-456f-93e1-327ab4c040ca' class='xr-section-summary-in' type='checkbox'  checked><label for='section-8dc03a84-0789-456f-93e1-327ab4c040ca' class='xr-section-summary' >Coordinates: <span>(3)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><ul class='xr-var-list'><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>time</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>datetime64[ns]</div><div class='xr-var-preview xr-preview'>1850-01-16T12:00:00 ... 2014-12-...</div><input id='attrs-c26922a9-c0a8-42d0-a062-74b17545547d' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-c26922a9-c0a8-42d0-a062-74b17545547d' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-7015b603-352a-4bd0-914b-9016cda254c2' class='xr-var-data-in' type='checkbox'><label for='data-7015b603-352a-4bd0-914b-9016cda254c2' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>bounds :</span></dt><dd>time_bnds</dd><dt><span>axis :</span></dt><dd>T</dd><dt><span>long_name :</span></dt><dd>time</dd><dt><span>standard_name :</span></dt><dd>time</dd></dl></div><div class='xr-var-data'><pre>array([&#x27;1850-01-16T12:00:00.000000000&#x27;, &#x27;1850-02-15T00:00:00.000000000&#x27;,\n",
       "       &#x27;1850-03-16T12:00:00.000000000&#x27;, ..., &#x27;2014-10-16T12:00:00.000000000&#x27;,\n",
       "       &#x27;2014-11-16T00:00:00.000000000&#x27;, &#x27;2014-12-16T12:00:00.000000000&#x27;],\n",
       "      dtype=&#x27;datetime64[ns]&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>lev</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>5.831e+03</div><input id='attrs-1258e578-57ec-49ba-a1ea-f6234fe66875' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-1258e578-57ec-49ba-a1ea-f6234fe66875' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-c810cc58-59d5-4b33-9b8d-1b56661a47e1' class='xr-var-data-in' type='checkbox'><label for='data-c810cc58-59d5-4b33-9b8d-1b56661a47e1' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>bounds :</span></dt><dd>lev_bnds</dd><dt><span>units :</span></dt><dd>m</dd><dt><span>axis :</span></dt><dd>Z</dd><dt><span>positive :</span></dt><dd>down</dd><dt><span>long_name :</span></dt><dd>ocean depth coordinate</dd><dt><span>standard_name :</span></dt><dd>depth</dd></dl></div><div class='xr-var-data'><pre>array(5831.29443359)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>ensemble_member</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>&lt;U8</div><div class='xr-var-preview xr-preview'>&#x27;r1i1p1f1&#x27;</div><input id='attrs-ff17f09e-3a09-446d-b805-27bf6ea8dab0' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-ff17f09e-3a09-446d-b805-27bf6ea8dab0' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-47e18339-bfc7-4594-9a61-bc96fa08d502' class='xr-var-data-in' type='checkbox'><label for='data-47e18339-bfc7-4594-9a61-bc96fa08d502' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array(&#x27;r1i1p1f1&#x27;, dtype=&#x27;&lt;U8&#x27;)</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-eb004d82-afeb-436a-8d32-6941b8092cbd' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-eb004d82-afeb-436a-8d32-6941b8092cbd' class='xr-section-summary'  title='Expand/collapse section'>Attributes: <span>(0)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><dl class='xr-attrs'></dl></div></li></ul></div></div>"
      ],
      "text/plain": [
       "<xarray.DataArray '__xarray_dataarray_variable__' (time: 1980)>\n",
       "array([2.59977719e+25, 2.60113345e+25, 2.60189363e+25, ...,\n",
       "       2.66573444e+25, 2.66691887e+25, 2.66822825e+25])\n",
       "Coordinates:\n",
       "  * time             (time) datetime64[ns] 1850-01-16T12:00:00 ... 2014-12-16...\n",
       "    lev              float64 5.831e+03\n",
       "    ensemble_member  <U8 'r1i1p1f1'"
      ]
     },
     "execution_count": 39,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "OHC_integrated = OHC.__xarray_dataarray_variable__.isel(ensemble_member=0).isel(lev=-1)\n",
    "OHC_integrated.load()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 63,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAWoAAAEwCAYAAACezTtyAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAABbr0lEQVR4nO2dd3gVxfrHP5NOEggQIPTepIbeIQEpioJdsWIBOyrX7s9ernotXLvYrw0v9iuoiBI6Il16Db0GEtLr+/tjNslJcpJsktMS5vM85zlnd2dnvlvOu7PvzLyjRASDwWAw+C5+3hZgMBgMhrIxhtpgMBh8HGOoDQaDwccxhtpgMBh8HGOoDQaDwccxhtpgMBh8HGOoDYZyUEq1VkqJUipFKTW1lDRPKqVSrXQBntZYU1BKxSmlbvJwmfFKqbOt3w8rpd73ZPl28FlDbV2wU0qpYCfbrlRKrbL+OIeVUj8rpYY6bO+olJqtlDqhlEpSSm1QSk1XSvkX+9M5fi639m2ulPrGYd+/lVKTHfK+USm1VSmVrJQ6qpSao5SqXcZx1FFKzVBK7bPK2WktN7C2xyulsvKXHfZbZ+lsbS3fp5TaaJW7Ryl1Xxll5h9jQLH1Hyulninv3JeHUipGKXXARrr+Sqm5SqlEpdRJpdRKpdT1Lih/slJqSVXzccjP7nmpKyIzrX2KnGMReRzoWsFyn7burxyl1BNOtt9pXevT1v3ueI9/bN03jvewv8P2aKXUaqVUmvUdXY6WmUqpbUqpPMf73dqmlFLPKKUOWv+JOKVU12JprlBKbbEeVruUUsMcto2y/jNpSqkFSqlWFTlPnkREnhMRjz4o7OCThtoyTsMAASYU2zYdmAE8B0QBLYG3gInW9nbAn8B+oLuIRACXAn0BR4NaV0TCHT5fWes/tfZtBUQC1wJHrbxHWOVOEpHawFnAf8s4jiDgd/QfeBxQBxgMJAD9HZLuASY57NcdqFU8O0tLPSuvO5RSV5RWtrdRSg0C/gAWAu3R5/JW4Bxv6vIxdgL3A3OKb1BKDQCeBy4BIoAPgO8cjTHwYrF7ONfaNwj4AfgMfb98AvxgrS+N9cBtwBon2y4FbkD/J+sDy9H/k3yto4EXgOvR/7HhwG5rWwPgW+BRa99VwFcYKoaI+NwHeAxYCrwC/OSwPgJIAS4tY9/PgDllbG+NfgAElLI9BYguZdu9wPcVOI6b0EY+vIw08cD/AX85rHsJeMTS2bqU/V4DXq/IMQIfA884LA8ElgGJ6D9qjMO264EtQDL6T3eztT4MSAfyrHOVAjR1omEJ8GY552cK2lidBH50zMfSfwuwAzgFvIl+WJ0FZAC5VtmJVvpg67zts875O0Ata1sMcAD4B3AMOAxcb22bCmQDWVZ+/7NzPu2us3mffAY8UWzd5cBKh+UwK+8mzq5lsX3HAAcB5bBuHzDOhpYlwORi6x4A/uuw3BXIcFheBtxYSn5TgWXFjiMd6FxK+jjgn8BKIAn9wKnvsH02cMTatgjo6rDtXGCzdc8eBO512HYesA59ry8DehT7D55t/X4C+KzY9bzOOn8ngEcc9vMDHgR2oStf/3XU6sqPT9ao0TXHz63PWKVUlLV+EBACfFfGvmcDX1eh7BXAm9arXMti2/609DyplBqinLhlnGj5RURSbJRZRyl1llVjuhz953WKUkqhazebysm3VJRSzdA1uWfQNZ17gW+UUg2tJMfQN3cdtNF+VSnVW0RS0bXiQ1JYkztULO9Q9LUq9ToopUai/5CXAU2AvcCsYsnOA/oBPa10Y0VkC9qAL7fKrmulfQHoCESja/DN0A/8fBqjH/TNgBvR17ieaFfG5xTWTs8v88RVAKXUW0qptyq5+8+Av1JqgHVP3IA2NEcc0txmuZRWK6UudljfFdggljWx2EAFXTMOzALaK+1SDEQbrl8ALG19gYaWW++AUuoNpVT+G2FXdCUAAOv+2VWOlmvRx9sUyEFXSvL5GegANELX/j932PYBukJRG+iGfqNDKdUb+BC4Gf1m9y7wo43/bz5DgU7AKOAxpdRZ1vppwAXACEtrfoXC5ficobb8cK3QT/DV6It6pbU5EjghIjllZBGJrjGVxwnLd5r/yT/5lwKL0a9qeyxfcT8AEVkMXAT0Rhu5BKXUK8VeRyujBfSr5LXAaGArukZQGk+gr91H5eRZ5BgpPI8AVwNzRWSuiOSJyG/o19JzAURkjojsEs1CYB764WCHepa+so79KuBDEVkjIpnAQ8Agy+2Vz/Mikigi+4AFaCNcAuvBNQW4R0ROikgy2kXl6BrKBp4SkWwRmYuuPXeyeTyVQkRuE5HbKrl7MvANuoabCTwOTHUwvq9RaLAeBT5WSg2xtoWja5yOJFHU9VcRDqP/E9vQteFLgXusbVFAINpFMwx9jXqh3xIrq+VTEdloGfVHgcvy/2Mi8qGIJFv3zBNAT6VUhLVfNtBFKVVHRE6JSL4bZwrwroj8KSK5IvIJ+pwOtHn8T4pIuoisRz90elrrb0bXsA846LlEuaEx2ecMNfppPU9ETljLX1jrQL9eNCjnRCSga2jl0UBE6jp8tgBYF/hBEemKvgnXAd9bxgAR+dmqddVH+8UnAzcppVo6NuxUUAtoQ32lld9/SkuklLoDbdDHWzeH7WNEn8t8WgGXFjPkQ/P1KqXOUUqtsGpsiWgD3gB7nEK7Rso69qboWjQA1ltHArrGm49j7TEN/ad3RkMgFFjtcCy/WOvzSSj2gC8rP1/gJnStsisQhH6w/qSUagpgPeASRCTHevB8jq5EgH4I1SmWXx208adYA2Txt0ZnPI5+s2mBfqN9EvjDenNKt9K8LiKHrf/tK1gP/PK0lMJ+h9970Q+CBkp3Bnjeaqw8jXZZQOF9ebFV7l6l1EKrnQT0vf6PYvd6C/Q9aIfS7sNW6HaD/Dy3oF1yUbgYnzLU1uvSZcAIpdQRpdQR9JO7p1KqJ7oRIwP9ulEa89EXrMpYN91L6Atav9i2PBH5Hf161U1E9jm4AvIv5Hy0qyTMRll70Y2K56IbX0qglLoB7RMbJSLl9rooh/3omovjwypMRJ63Xgm/QR97lGXk56J9xKD9dmUdSxr6WpV1HQ6hb3QArHMUSdlvEgVFFFs+gTYYXR2OJcLhOlQ0P1+gJ9pfvt26135B12wHl5JeKLw+m4Ae+ZULix7WeqRoA+Q+m1q+smqOOSLyMfqtqYuInEL7/0s7h5sorIHmX+d2lO22a+HwuyW6pnwCXZGZiHYpRqB9yGAdt4j8JSIT0W8Z31PY0L8feLbYvR4qIl/aOPay2A+cUyzfEBGxcw9XCJ8y1GgDnAt0Qb9CRaMbjxYD14pIEtrv+KZS6gKlVKhSKtCq/b1o5fE4MFgp9S+lVGMApVR7pdRnSqm65QlQSr2glOqmlApQutvdrcBOEUlQSk20fNf1lKY/2j+1opTs8nuQfKOU6qyU8lNKRSrdV/NcJ+lvBEZar3zFdV2Ffp0fLSK7yzsOG3wGnK+UGmvVVEKU7nbXHF2DCwaOAzlKqXPQDVT5HAUiHV45nXE/MFnpboWR1jH0VErl+6G/AK5XuhtZsHVsf4pIvA3tR4HmyurFICJ5wHtoP3ojq6xmSqmxts6Ezq+tzbQuw7p3Q9D/wwDrGuS70f4Cxiul2lr32mi0D36jte8lSqlw654ag65x/2jtG4f+H01TSgVbb2Fg+WxL0RJkaVFAoKUl3z78hX77irLKuwZdy91pbf8IuFMp1UgpVQ+4G/jJ2vYd0E0pdbGV/2No//nWMk7N1UqpLlaN/Snga9E9WmqjXRYJ6Deo54rpv0opFSEi2cBp6xyAvjduUdrfr5RSYUqp8aqMbrU2eQd4VlndDZVSDZVSE6uYp3PEDS2Ulf2gX1dfdrL+MvTrR4C1fBXan5pqrZ8DDHZI3wndOpyA9oetR988/hS25KYU+0y39n0d3dMgBW2ofgLOsrYNR3e3O4F+ddsO3F/OMUWguxPut/LchX41jJRiLc7F9gvAodcHuradXUzzO6WUmX+M5fX6GIDuPnfSOtY5QEtr2+1oA5aIfuDMKrbvh9b5TcRJrw8rTX9040+SVcaf6Adu/vZbrPNx0jrPzR22CdDemXb0g2SOtd8Ja10I+o+7G/0n3QJMs7bFAAeKaSs472hf7zrrWL63cz4d1vmXlg79R3Z6jRyOSYp9JlvbFNpI7UPfa1uAaxz2XWyd19Po+/uKYnn3Alaj3zTWAL3KuU/jnGiJcTi3b6Jr9Ket/MY57BuI7iKbiP4/vgaEOGw/G93ukm6V07ocHfm9Pk4D/0O78EC7HH6wzsdetAtQ0I3HQWj7ccra7y9gqEO+46x1idZxzAZqO7kXnqBkr4+AYvpusn77AdPRvvtk9L38nDtso7IKNBgMpWDVmLah3W73ich7SqkewCKxep0opR5H/2mDgTCx+jQbDK7AGGqDoYJYLoGngIEicra39RhqPiYmgcFQcfahX/GrPBzeYLCDqVEbDAaDj+NrvT4MBoPBUAxjqA0Gg8HHMYbaYDAYfBxjqA0Gg8HHMYbaYDAYfBy3GWql1IdKqWNKqY020k5XSm1WeiaW3x2GZLZSOoTjOqXUJqXULe7SazAYDL6K27rnKaWGo4c5/0dEupWTNhYd5yFNKXUreujq5VYsByUimUqpcHScg8FSLP6xwWAw1GTcVqMWkUXoWAwFKKXaKaV+sWrJi5VSna20C0RHXAMd4Ki5tT5LCkN5BrtTr8FgMPgqnjZ8M4E7RaQPekYRZ7Nf3IgO5AOAUqqFUmoDOqjRC6Y2bTAYzjQ8NoTccl0MBmY7hMkNLpbmavS0PiPy14nIfnRs3aboAP5fi8hRz6g2GAwG7+PJWB9+6IlIo51tVEqdjZ7QdYQ4mblERA4ppTahp/upypyIBoPBUK3wmOtDRE6j5yC8FPQ8d0rP2oJSqhd6wskJInIsfx+lVHNlTZJpBSQfgg43aTAYDGcM7uz18SU6YHsDdAD6x9EzTLyNnksvEJglIk8ppeYD3SmcDHWfiEywZrV4mcJpht4QPWu0wWAwnDGY6HkGg8Hg45jubgaDweDjGENtMBgMPo5ben00aNBAWrdu7Y6sbZOamkpYWJhXNRTHaLKH0WQPo8kevqapND2rV68+ISINne7kjhlz+/TpI95mwYIF3pZQAqPJHkaTPYwme/iaptL0AKukFJtqXB8Gg8Hg4xhDbTAYDD6OMdQGg8Hg43hsCHl2djYHDhwgIyPDI+VFRESwZcsWj5RlF6PJHkaTPUrTFBISQvPmzQkMDPSCKoM78JihPnDgALVr16Z169Y4BGVyG8nJydSuXdvt5VQEo8keRpM9nGkSERISEjhw4ABt2rTxkjKDqynX9aGU6mTNsJL/Oa2UuruiBWVkZBAZGekRI20wnKkopYiMjPTYm6vBM5RboxaRbUA0gFLKHzgIfFeZwoyRNhjcj/mf1Twq2pg4CtglInvdIcZdJCQkEB0dTXR0NI0bN6ZZs2ZER0cTHh7Obbfd5pYyBw8eXG6aGTNmkJaWVm66qhIfH88XX3xR4f0mT57M11+XjCgrIjzzzDN06NCBjh07Ehsby6ZNm1whtcI40xgeHl7p/J577rlSt6WkpHDzzTfTrl07unbtyvDhw/nzzz8rVc7HH3/MoUOVnwNj3bp1/Prrr5Xe3+Babv98De8v3u22/Cvqo74C+NLZBqXUVGAqQFRUFHFxcUW2R0REkJycXAmJlSM3N7egvKCgIBYvXgzoP2J4eDjTpk0rSOsOXb/++muJfB01Abz66qtccMEFREZGVrm8nJwcAgKcX87Nmzfzn//8h/PPP7/EtuKaHMnOziY9Pb3E9nfffZdFixaxZMkSQkND+f333znvvPNYuXIlISEhVT6WsjTZ1VjZa/rcc89x5513OtV044030rp1a9asWYOfnx979uxh27ZtdOnSpcLlfPDBB7Rp06bSfu8VK1awevVqxo4d63R7RkZGif+gJ0hJSfFKuWXhbk15Isz5O405fx+mfe4+9+gpbSRM8Q8QBJwAospL62xk4ubNm10wpsc+p0+fdrr+8ccfl3/9618iokcIjR8/vmD9tddeK6NHj5ZWrVrJN998I/fdd59069ZNxo4dK1lZWSIismrVKhk+fLj07t1bxowZI4cOHXJaTlhYWEEZI0aMkIsvvlg6dOggV155peTl5cm///1vCQwMlG7duklMTIyIiPz6668ycOBA6dWrl1xyySWSnJwsIiJz5syRTp06yZAhQ+TOO+8sonnKlCkyevRomTRpkuzZs0eGDh0qvXr1kl69esnSpUtFRGTAgAFSp04d6dmzp7zyyiuSk5Mj9957r/Tt21e6du0q77zzjoiI5OXlye233y5nnXWWnHvuuXLOOefI7NmzSxxb8+bNZefOnUXWXX311fL++++LiMjPP/8svXr1kh49esjIkSNFRCQlJUWuv/566du3r0RHR8v3338vIuJU8+nTp4uct06dOhWct+Jcd911JTTmn3sRkRdffFH69u0r3bt3l8cee6xg/cSJE6V3797SpUsXeffdd0VE5IEHHhA/Pz/p2bOnXHnllUXyXLdunbRu3VpycnKcXu+XX35ZunbtKl27dpVXX3214Ng6d+4sN910k3Tp0kVGjx4taWlpMnv2bAkLC5OOHTtKz549JS0trdT7asSIEXL//fdLv379pEOHDrJo0SLJzMyUFi1aSGRkpPTs2VNmzZpVQo+n/2/5+NooQBH3a0rNzJZWD/wkrR74qUp6KGNkYkUM9URgnp205RnqJ37cKJe9s8ylnyd+3FikvMoY6iFDhkhWVpasW7dOatWqJXPnzhURkQsuuEC+++47ycrKkkGDBsmxY8dERGTWrFly/fXXOy3H0VDXqVNH9u/fL4mJiTJw4EBZvHixiIi0atVKjh8/LiIix48fl2HDhklKSoqIiDz//PPy5JNPSnp6ujRv3lx2794tIiJXXHFFEc29e/eWtLQ0ERFJTU2V9PR0ERHZvn275F8Hx+MUEXn33Xfl6aefLii3T58+snv3bvnmm2/k7LPPlpycHDl48KBERESUMIJJSUlSr169Esc7Y8YMueeee+TYsWNF9CYkJIiIyEMPPSSffvqpiIicOnVKOnToICkpKU415xvq/POWm5tb5Lw5ct1110nr1q2lZ8+eBZ/8c//rr7/KlClTJC8vT3Jzc2X8+PGycOHCIrrS0tKka9eucuLEiSLXrTizZs2SCy64wOm2VatWSbdu3SQlJUWSk5OlS5cusmbNGtmzZ4/4+/vL2rVrRUTk0ksvLTgHI0aMkL/++ktEpMz7asSIETJ9+nQR0Q/sUaNGiYjIRx99JFOmTHGqR8QYakfcrel4cobbDXVFXB+TKMXtUVM455xzCAwMpHv37uTm5jJu3DgAunfvTnx8PNu2bWPjxo2MHj0a0K/DTZo0KTff/v3707x5c5KTk4mOjiY+Pp6hQ4cWSbNixQo2b97MkCFDAMjKymLQoEFs3bqVtm3bFnS1mjRpEjNnFs6dMGHCBGrVqgVoN8Add9zBunXr8Pf3Z/v27U71zJs3jw0bNvD111+Tl5dHcnIyO3bsYNGiRUyaNAl/f3+aNm3KyJEjbZ87EUEpxYoVKxg+fHiB3vr16xeU+eOPP/LSSy8B+tV83759NG3atFTN+ecNKPW8AfzrX//ikksuKVjO91HPmzePefPm0atXL0C/cu7YsYPhw4fz2muv8d13uk18//797Nixo9IuqCVLlnDhhRcWBNq56KKLWLx4MRMmTKBNmzZER0cD0KdPH+Lj40vsX959ddFFF5W5v8G7pGflur0MW4ZaKRUKjAZudkWhj5/f1RXZuJzgYD3Xrp+fH4GBgQWt535+fuTk5CAidO3aleXLlxfZb//+/QX+31tuuYVbbrnFab4A/v7+5OTklChbRBg9ejRffln0Wbh27doyNTtG4Xr11VeJiopi/fr15OXlleovFhFef/11xo4dW6Qv7ty5c8vtMVCnTh3CwsLYvXs3bdu2LVi/Zs0aRowYUWCwnZX5zTff0KlTpyLrn3jiiVI12zlvZSEiPPTQQ9x8c9HbNi4ujvnz57N8+XJCQ0OJiYkptztb586dCzT6+RVtg5cyJt8ofgzp6elOdTq7r4rnUZlzYHAvmTm5nEgpnOL1pk/+4v3r+rm8HFu9PkQkTUQiRSTJ5QqqEZ06deL48eMFf6js7Gw2bdpEixYtWLduHevWrSthpMuidu3aBY1eAwcOZOnSpezcuROAtLQ0tm/fTufOndm9e3dBTeqrr74qNb+kpCSaNGmCn58fn376Kbm5uSXKARg7dixvv/022dnZAGzfvp3U1FSGDx/OrFmzyM3N5fDhwyxYsMBpOffddx/Tpk0rMDrz589nyZIlXHnllQwaNIiFCxeyZ88eAE6ePFlQ5uuvv15g1PIfQKVpdgVjx47lww8/JCUlBYCDBw9y7NgxkpKSqFevHqGhoWzdupUVK1YU7BMYGFhwXhxp27Ytffv25fHHHy84hh07dvDDDz8wfPhwvv/+e9LS0khNTeW7775j2LBhZWpzvCal3Vfl7Z9/XAbvMerlhVz41rKC5c2HTrulHBProwIEBQXx9ddf88ADD9CzZ0+io6NZtmxZ+TuWwtSpUznnnHOIjY2lYcOGfPzxx0yaNIkePXowcOBAtm7dSq1atXjrrbcYN24cQ4cOJSoqioiICKf53XbbbXzyyScMHDiQ7du3F9S2e/ToQUBAAD179uTVV1/lpptuokuXLvTu3ZsBAwZw8803k5OTw4UXXkiHDh3o3r07t956KyNGjHBazp133km/fv3o3r07nTp14umnn+aHH36gVq1aNGzYkJkzZ3LRRRfRs2dPLr/8cgAeffRRsrOz6dGjB926dePRRx8tU7MrGDNmTMHDo3v37lxyySUkJyczbtw4cnJy6NGjB48++igDBw4sck169OjBVVddVSK/999/nyNHjtC+fXu6d+/OlClTaNq0Kb1792by5Mn079+fAQMGcNNNNxW4W0pj8uTJ3HLLLURHR5Obm1vh+yo2NpatW7cSHR1d5sPb4F4OnCr6hhQS5O+WctwyZ2Lfvn1l1apVRdZt2bKFs846y+VllUZ1GfJrh5SUFMLDwxERbr/9djp06MA999zjVU3uxGiyR1maPP1/yycuLo6YmBiPl1sW7tTU+sE5RZbbNwpn/nTnFZzy9CilVotIX2f7mBp1NeC9994jOjqarl27kpSUVMLnajAYPMviHceZvWp/ifXuGhPqsaBMhspzzz33uKwGbTAYqs41H6x0ut7fzz2m2hhqg8FgsMmhxHTW7U8sdXuAvzHUBoPB4FVenredb9YcKHV7gJ97vMnGR20wGAw2WRmfUOb2QDfVqI2hNhgMBpsE+pdtMnu3queWcs8IQ23CnLo2zGlMTAzFu1/mk98P/NVXXy2hoVu3bhXWYDD4Alk5eazeexLK6M38+U0DuD22vVvKPyN81JGRkaxbtw7QQ5bDw8O599573VqmnYEwM2bM4OqrryY0NLTK5ZUV5jTfUF955ZVVLqcsjhw5wrJly9i7t1qFKzcYyuX7dQe5/+sNZaYZ0r6B28o/I2rUpREXF8d5550HaAN+3XXXMWbMGFq3bs23337L/fffT/fu3Rk3blzBsOLVq1czYsQI+vTpw9ixYzl8+LDTvPMDA+V3br/kkkvo06cPV111FSLCa6+9xqFDh4iNjSU2NhbQQYQGDRpE7969ufTSSwuGCM+dO5fOnTszdOhQpk2bVkTz1KlTGTNmDNdeey3x8fEMGzaM3r1707t374KHxYMPPsjixYuJjo7m1VdfJTc3l/vuu49+/foxaNAg3n33XUDHnLjjjjvo0qUL48eP59ixY6Weu9mzZ9O/f386duxYEOt7zJgxHDt2jOjoaBYvXszq1avp2bMngwYN4s033yzYd9OmTfTv35/o6Gh69OjBjh07KncBDQYPsemgd6NneKdG/fODcORv1+bZuDuc83yVsti1axcLFixg8+bNDBo0iG+++YYXX3yRCy+8kDlz5jB+/HjuvPNOfvjhBxo2bMhXX33FI488wocfflhmvmvXrmXTpk3Url2bcePGsXTpUqZNm8Yrr7zCggULaNCgASdOnOCZZ55h/vz5hIWF8cILL/DKK69w//33c/PNN7No0SLatGnDpEmTiuS9evVqlixZQq1atUhLS+O3334jJCSEHTt2MGnSJFatWsXzzz/PSy+9xE8//QTAzJkziYiI4K+//uLEiROMGzeOMWPGsHbtWrZt28bff//N0aNH6dKlCzfccIPTY8rJyWHlypXMnTuXJ598kvnz5/Pjjz9y3nnnFby99OjRg9dff50RI0Zw3333Fez7zjvvcNddd3HVVVeRlZXl0vgeBoM7KM837W7OCNeHXUyYU/thTssLvZmUlERiYmJBvJBrrrmGn3/+GYBBgwbx7LPPcuDAAS666CI6dOhQ7jk0GLxBbp7w3Nwt7Dru3QBY3jHUVaz5ugsT5rT8MKfFj6ms4yktryuvvJIBAwYwZ84cxo4dy/vvv1+h2NcGg6fYdCiJD5bs8baMM9tHXVFMmFP71K1bl4iICJYsWQLA559/XrAtP5b1tGnTmDBhAhs2lN1IYzB4i5TM8uN/j+0axeuTyo6WWFWM66MC5Ic5nTZtGklJSeTk5HD33XfTtWvlJkLID3PapEkTFixYUBDmNDNTByJ/5pln6NixY0GY0wYNGtC/f/9S87vtttu4+OKLmT17NrGxsU7DnE6ePJm77rqL+Ph4evfuTW5uLlFRUXz//fdceOGF/PHHH3Tv3p2OHTuWGubULh999BE33HADoaGhRSZh/eqrr/jss88IDAykcePGPPbYY1Uqx2BwF2mZ5befvHuN04B3LsWEOfUgJsypPYwme5gwp/aorKadx5JZsy+x3G558c+Pd4messKc2p2Kqy7wPtAN3eX7BhFxPm+QweW89957fPLJJ2RlZdGrVy8T5tRgcDNpWTmc/coib8sowK7r49/ALyJyiVIqCKj6CA2DbUyYU4PBs5xO9625Kcs11EqpOsBwYDKAiGQBWe6VZTAYDN4jLat8Q/3R9f1oGlHLA2rs1ajbAseBj5RSPYHVwF0iklrRwsrqsmUwGFyDO9qdziQysnM5lVZ6XfScbo3p3LgOsZ0aeUxTuY2JSqm+wApgiIj8qZT6N3BaRB4tlm4qMBUgKiqqz6xZs4rkEx4eXjAxqyeMdW5uLv7+7plosrIYTfYwmuzhTJOIkJSUxNGjR70yS3l+w7cvUVFNt85PpSzPxyMDQuhQr/L3Qml6YmNjS21MtGOoGwMrRKS1tTwMeFBESm3qdNbrIzs7mwMHDpCRkVHecbiEjIyMUgd8eAujyR5Gkz1K0xQSEkLz5s0JDAz0uKaa0Ouj+IS1xZl3z3A6RlW+B5Bben2IyBGl1H6lVCcR2QaMAjZXVFxgYGDBMGhPEBcXR69e7u2EXlGMJnsYTfbwRU1nAnVCPP8AtNvr407gc6vHx27gevdJMhgMBs/zdtwu9p0svemtY1Q4tYICqB8W5EFVGluGWkTWAe4ffmMwGAweJjkjm/0n03nhl61lprt2UGuuHtjKQ6qKYoaQGwyGM5q7Z63j962lx17PJyTQe43JJiiTwWA4o7FjpAHCg71nqE2N2mAwnJGczshmw/7yZ24Z3rEh7RuGM6Kj5/pNF8cYaoPBcEby2vwdvG8j1nSj2sE8dn4XDygqHeP6MBgMZyQbD9mbBzHQ3/ujqY2hNhgMZyR2R9qHBXnf8eB9BQaDweBBkjOy+WRZPMkZ5QdeumZgK+4Y2d4DqsrGGGqDwXBG8dvmo7w0z/nEz44EBfjx9AXdPKCofIzrw2AwnFEkpmXbShfo533fdD7GUBsMhjOGtKwc8mw6p/19yFAb14fBYDgjyMjOpctjvxLkb69+OqxjQzcrso8x1AaD4YzgYGI6AFm5eeWm/e2e4TSO8J2wtsZQGwyGM4K0zFzbaTtUId60OzA+aoPBUONJzsgmOdNeI6IvYmrUBoOhxtP9iXkEBdirl/rCSMTiGENtMBhqNDl5updHVk75vukXL+5Bk7q+45vOxxhqg8FQo7Hjmr49th11QgK5rF8L9wuqBMZQGwyGGsufuxPYcap8S33DkDZEhgd7QFHlsGWolVLxQDKQC+SUNlOuwWAw+BKXz1xhK13dUM/Pg1gRKlKjjhWRE25TYjAYDB5m2sj2RLes61OjEJ1hXB8Gg6HGEX8ilfUHEstNN7BdJIPbNXC/oCqixMa4d6XUHuAUIMC7IjLTSZqpwFSAqKioPrNmzXKx1IqRkpJCeHi4VzUUx2iyh9FkD6OpdG6al4qNTh48N7QWTcM9O5yktHMUGxu7ujS3sl1D3VREDimlGgG/AXeKyKLS0vft21dWrVplX7kbiIuLIyYmxqsaimM02cNosofRVDqtH5xT5vZJ/VsyoE19JkY3RSnPuj1KO0dKqVINta1HiYgcsr6PAd8B/Ssv02AwGLxLnZAALujVzONGurKUa6iVUmFKqdr5v4ExwEZ3CzMYDIaKIiIcT84sN11woL8H1LgOO42JUcB31pMnAPhCRH5xqyqDwWCoBK/9vpNX55c/e0vnxr4VdKk8yjXUIrIb6OkBLQaDwVAlflh/sMztLerX4rkLuzOwbaSHFLkG0z3PYDDUCNKzcgn0K9ubG+jnx7AOvjMhgF2MoTYYDNWe37cc5cZPyu9pFmhzdhdfo3qqNhgMBgcWbT9uK9053Ru7WYl7MDVqg8FQ7bETa/rb2wbTpUkdD6hxPcZQGwyGaktCSiYT3lhKq8jQctP2blnPA4rcgzHUBoOh2rJq7ykOJqYXTFxbUzE+aoPB4JMs2HaMQf/8vcwBLEHVtHGwopgatcFg8ClEhMNJGXzx5z4OJ2Ww5fBpGtZ23qUuI7v8SQGu6hzE+cN7u1qmRzGG2mAw+BQrdp9k0nuFAf9znQSOy80Tzn5lIcE2GhFHtw6sdgNcinNmvDcYDIZqw85jyUWW9yWksW5/YpF1iWlZ7DmRytYjRdMWZ1TnRq6W5xVMjdpgMPgU/sVGFz7+4yYAdj13Lv5+ih1Hk9mbkFZuPisfHkVkeDCLFy10i05PYgy1wWDwCfLyhKfnbCYhJcvp9nmbjpAncPsXa2zlVy8syOen2LKLMdQGg8EnOHAqnY+Wxpe6/dbP7RnofKrrcHFnGENtMBh8gpTMHJfkUzc0kFaRYS7Jy1cwhtpgMHidD5bs4eAp1wxaeevK3gxqV717eRTHGGqDweB1nv5ps8vyqlMrsNpMsWWXmuPEMRgM1ZIsO9OFV4CGtYNdmp8vYLtGrZTyB1YBB0XkPPdJMhgMZwq5eeIy3/SdI9tTJySQRmeyoQbuArYA1TNOoMFg8DrZuXkkZ+SwbNcJ4k+k8sfWYySkOu+OV1G6NKnDOd2buCQvX8OWoVZKNQfGA88C092qyGAw1BhEhDX7EukYFc6yXQks23mCT5bvdUtZNak7XnHs1qhnAPcD1WvqXoPB4FU2Hz7NxW8vIyTQj4xs1/qi82lYO5iuTevQu1X1jTddHkqcBDwpkkCp84BzReQ2pVQMcK8zH7VSaiowFSAqKqrPrFmzXK+2AqSkpBAeHu5VDcUxmuxhNNnDVzWFhYWhlGLRgWyycuGzLa5xbZRGdEN/7u4TUqYmXzpPpemJjY1dLSJ9ne1jp0Y9BJiglDoXCAHqKKU+E5GrHROJyExgJkDfvn0lJiamgvJdS1xcHN7WUByjyR5Gkz18VdNrmwMJ8Pdj5Z5Uj5TZtHFDYmL6lKnJl85TZfSUa6hF5CHgIQCHGvXVZe1jMBjOPE6mZnE8LY81+xI9Wu5ZjWt+/wYz4MVgMFSJvDwhIyeXS95Zxv4E106JNTG6Kav3nuKAk1GL/n6Kf17UnbFdqufM4hWhQoZaROKAOLcoMRgM1ZLHftzIrJX7yckru72rIlzRrwWN6oRw16gO+PspWj84p0SaWoH+XNa3hcvK9GVMjdpgMFSJz1bsc3meUXVCuGd0xzLT9GpZ1+Xl+irGUBsMhkqx9chpxs1Y7Ja8I8ODytz+6uU9q/30WhXBGGqDwVBhcvOE79YcdHm+bRqEMXV4Wy7s1azMdBf2au7ysqvE2s+h5UCIbOeW7I2hNhiqIXtOpLLjVC7NjiYTHOBPy8hQj5SbkJLJydQsRr+6iK5NXd/bolagP5P6tyyx/qVLe7LvZBojOzciPNjf5eVWibw8+OE28A+CR4+7pQhjqA2GasiE15eQnJnDs38uAiD++fFuLS8vT0hIzaLfs/NpEqEHl2w6dNrl5XRvFuF0/SV9fKwG7Ui2NX9jrvsG9hhDbTBUA7Jy8li99xTdm0cQfyKV5GIR5z5dHk9EaBATejZ1abl5eUKeCO0f+Zl+rfUQ7cNJGS4tI59nLujGyOo4a3hWituLMIbaYKgG/G/9If4xez0t64ey72TJGbgf/UHP1L14+3GmDm9Lh6iqheXZeDCJg4npfP7nPpLSswH4K/5UlfIsj6sHtnJr/m4jy2EE5jdT4OL3XF6EMdQGgw+TkJLJhDeWEhSgI8M5M9KOzF59gGW7Elj64MhS09zz1Tpa1KvF9DGdyMsTvl17kHO7NyY0KICfNhzixV+2kZiWxekM18SJrtF8fxukOzzA9i51SzHGUBsMPkx8QioHEys22i8713mUutw8QUT4bq3urbHuQBJdmtThnYW7mLfpCAPaRvLd2gPlPgxczZ0j29OzeV2Pluky1n1edDnIPZPqGkNtMPgIS3ee4O+DSVzapzkhgf7847/raVG/VoXz8fdzPl/gP/67juW7EwqWF20/zqLtupfCvM1Hmbf5aOWEV5F/jOnklXLdgjHUBkPNJCM7l/X7E7nq/T8BeP7nrUTUCizwDVcUP6X4K/4kfVrWw8/BaH+/7pBL9FaVsCB/UrNyvS2japzYCZlOer0EGkNtMNQoluw4wY2f/MVZTeqwbn9ikW2VNdIABxPTufSd5Tx6XhduHNqGbUeS+XjZniqqrTwjOzdiVfxJujWLYNmuBBbcGwMKdhxNIa+cePg+yxulhFUNdk/ca2OoDQY3ciw5g5xcoWndQhfGiZRM7v96A9uOJJOZk1fCSLuK/yyP571FuzmZluXymb7t8uHkvgxsG0loUAC5ecL+k2k0qqP7YTeqXXqwf59l47ew6sPStwe6Z+CRMdQGgxs577Ul5OYJq/7vbJRSLNh6jENJ6fyx9Zjby96b4NlGweLcfXYHRnaOKlj291O0buAe14DH+Pr6srf7B7qlWGOoDQY3ciw5E4AJbyylbcMwfvARP7G72f3cuUX842cMfsZQGww+T16esPnwaboVGwr998Ek/j6Y5CVVnuOiDoEM6NG55hnplGPw1/vlp/NzTxwSY6gNBhfy7qLdvPDLViYPbo2fqmHGqgzuH9eJ3Fyhu/9BYvqVDKpU7Zl7L2z+ofx04p62AGOoDQYXkJmTy89/H2H1Xj1K7eNl8d4V5GGmDmtLgL8fcXGuD33qE6Qn2kvXc5Jbii/XUCulQoBFQLCV/msRedwtagwGHyYhJZPsXGHxjuO0qB9aELh+w4FE/tx9kmfnbvGyQs/z/rV9iW5ZlwB/P29L8T5PuM+1ZadGnQmMFJEUpVQgsEQp9bOIrHCbKoPBB7nji7VFRvZtfXocoBsKz1Qa1QmmQXiwt2W4j7ST8NnFkHrCqzLKNdQiIkB+HL9A61NNe6kbDJXH0UgD9HrqN0aeVQ3DcrqASf1bEOjvR+fGrp88wKfY9jMcWuNtFSixMTJIKeUPrAbaA2+KyANO0kwFpgJERUX1mTVrloulVoyUlBTCw90zSqiyGE328DVN8/dm88POTJKzfatx8OxmQlhIED/sqvwoxopSJwiGNQtkfNtAQgNLng9fu3ZQNU0Njy2m6+aXbKWNi7HR2FiGntjY2NUi0tfZPrYaE0UkF4hWStUFvlNKdRORjcXSzARmAvTt21diYmJsiXYXcXFxeFtDcYwme/iKptTMHD5csodvdu4i3ceMNMBV3cKIjY3lxgOJHEnKYOqnq91a3uguUdw4tE2Zk8r6yrVzpNKaNn0HwZGw2U5iZbuMyuipUK8PEUlUSsUB44CN5SQ3GKo1L83bxkdL470to4DHzuvC6YxsrhnYiqOnMzm2Xb+S92hel45RudQOCaBH8wiW7kwoJ6eyaRUZSkiAP9uOJgMwbWR7ejSvy9ldosrZswaRmgCzJ9tLGxQO57zoVjl2en00BLItI10LOBt4wa2qDAYvk5qZw3FrVKG3ie3UkPE9mnJBdNOC3hWR4cEc216YJiTQn7+fGMt3aw9U2VC/dkUvujeL4KK3lxEW7M/0mhSG1C7ZqeWnyedh93dJtFOjbgJ8Yvmp/YD/ishP7pVlMHiP/SfTGPbiAm/LKKBh7WDbk7v6+9nrJvfwuZ1Jzsjh9T92EhLoR0Z2HkH+fmTl5tGoTjB+forvbhuMOoMG7RQh0/3zIFYEO70+NgC9PKDFYPA6SWnZbDjgO0O9ezaPqNBcgkPaRdIqMpSk9GwS05w3MrZvFM7U4e3IzMklOSOH64e0JijAj9DAABbuOE6TCB3p74w10rsW6CHjPoQZmWgwoKepSs3KoedT82hcx/PhN/9v/FmkZOYwqX9LsnLyGPbiAno0j+CHO4ZWKJ/I8GAW3hfLZe8sZ2X8yRLblz80kkDLfRIc4M8TE7oW2e7qWcyrJZ9eYD9t835uk+GIMdSGM5rdx1M4kZLF3bPWFsRJPnI6w+M6rh3UumACW4DPbxpA24aVDwmanVc05sSzF3bDT6mC2rLBCTlZ9tJFtoc6zeCi99w2UUBxjKE2nJEs2HaMbUeSef7nrQXrDiV53kD/9+ZBhAT6FTHSAEPaN6hSvtNGdeD6j/4qWL6wVzNCg8zf3Sl5eZCXAzO62wuqdNVsqN/W/bocMFfOcMaRk5tXxIi5m3vHdCQnT5gxfwcAQf5+zL5lEGlZufRvU98tZcZ2akT88+NZu+8U/1t/mFqB7gm/WSP4erK9yHj5BHl+QI8x1IYzio+W7uHJ/9kaweAy7hjZgeSMbGbM38G0UR2YPrqjx8ru1bIevVrW81h51ZKKGGmAkLpukVEWxlAbzig+W7HX7WWM79GE3Fxh0Y7j1AsNAqB2SCAbnhhDmHE/+A4px2DHPPvpe1yhJwYICHKfplIwd43hjEBEWLDtGNm57o8nduuIdnRrFkFqZk6R6GV1QtwzTZOhksy6Eg5UwAXW6ypoM9x9esrAGGrDGcGP6w9x16x1Hikryuo9EhZs/l4+TcLOiqX3gm86H3MnGWo82bl5HEp0f4+Oly/tSZemdWhYuwbHZ64JpJ2E35+E3Bx76Zv1Bb8AiOpaflo3YQy1oUZz7HQG/Z/7nXZV6JPsyPKHRuKvFP2f+52J0U35fcsxZlweTVJ6NhMdYnEYfJi4f8Lqj+2nbzUIxjzjNjl2MIbaUOP4Yd1BImoFUj8sqCCw0q7jFQiy44TLOgbSrl27ggEje/55LkopROTMHWpdXclKq1h65f2HrzHUhhqHO3zRI1sGMm5Eu4LlfONsjHQ1Y+sckNyK7RPRwj1aKoAx1AZDGTx/UXeiIkJQhz3b99rgBo5u0j097BIeBUOnQ5/r3afJJsZQG2oMz87ZzHoXR75rXi+UoR0aEGcMdfUn12Ysj3wCQmDgLe7RUkGMoTbUGN5bvMdleY3v0YTWkaEMaOueId4GL5CdXrH0PuCbzscYakO1Z+WekyzZcdyleUbVDuG+sZ1dmqfBi/w0HQIrGDmwgeeG+peHMdSGasumQ0k889MWVuxJQFw84LBXy7quzdDgNVReDqz6oGI7nfsS9LjMPYIqgZ05E1sA/wEaA3nATBH5t7uFGQzl8cmyeJbvrtr8gMVpEhHCl1MG0rJ+qEvzNXiPgJxKdM3sP8X1QqqAHSdMDvAPETkLGAjcrpTq4l5ZBkPZpGflkmsjdLBdhlrxn+8Y2Z7WDcLw8zPd7moEG7+h4fFl3lZRZezMmXgYOGz9TlZKbQGaAaYZ3OAVth1JZuyMRYQEVq2x5/VJvfBTihGdGhLk74dSFExTZaghfH0DFfY097vJHUqqRIV81Eqp1uiJbv90ixqDoRQS07II9Pfj2g9XUi9UR6HLyK5alfq8Hk3MgJWayoHVkHbCfvpWQ6BuK5j4JvjgPaHEZiuMUiocWAg8KyLfOtk+FZgKEBUV1WfWrFmu1FlhUlJSCA/3XrQrZxhN9nCmafIvqdQPUZzMqHqr4ZWdg1DA6Nb2w45Wl/PkbXxFU0zcxAqlj4up4OQBVaC0cxQbG7taRPo628dWjVopFQh8A3zuzEgDiMhMYCZA3759JSYmxqZs9xAXF4e3NRTHaLKHo6Z3F+6y1m51iZEGuOHcgbRvVLvSmnwFo6kM4iqW3JOaK3OO7PT6UMAHwBYReaVy0gyGyvFPh8lnq8qUYW3IysmjdaRrIukZagC3LoN6bbytolzs1KiHANcAfyul1lnrHhaRuW5TZThjScvKYU9SLv47jvPd2oMuzbtT4zpc0qe5S/M0+BgHV8P2X+2nr9cGgny/K6adXh9LAN/zrhtqJO8t2sOryzNg+UqX5dm4TghtG4YxvEMDl+Vp8FE+Pg+ybYQxHX6fnhCgGhhpMCMTDT5ASmYOuXnCp8vjmfP3IZfn3yEqnE9vHODyfA0+iB0jDVCnKXQa514tLsQYaoNXyZ+BJTIsiITUCkY3s0nDcDM1Vo1nx3z4/GL76f2r1z1hDLXBqxw5recydLWRXvnwKFDwwNcbuG5wa5fmbfAQOVmQuBcadCg/7ZqPK5Z32xGVkuQtjKE2eJyktGyCAvyY8/dhDidWMPRkOQT6K9o1DKeRNRP4R9f3d2n+Bg/y8316bsP7dkNYpPM0IrDlf+Bns0982xi41nN9pl2FMdQGl5GTm8cFby3lyv6tqBsayLAODagdUvgH+mjpHnYcS+GLP/fRqHYwx6z5DF3J0gdGEhFqfyCLwQdJPwXzn4RN3+vllCOlG+oN/4XvptrPW/lXWZ43MIba4BJGvRxH68gwNh48zcPf/Q3Aud0bk5MrPHNBN4ID/Xnyf4XhYdxhpIGCmrShGrP0NVj9UeFyRhmz9iTsrFjeHUZXTpOXMYbaUCWe+t9mElIz2XU8tcRM33P/PgLAH1uPkZPn4oDRxfj3FdFk5rgwnJ7BixS7V0rryZG4r2IT1d6yFKK6Vl6WFzGG2kB6Vi5r952iZWQo244k46cUqVk5nNejaan7rNufyO7jKXy4tPzpr9xppN+8sjd9W9cjytSkqz+ZyfB6XwhvWHT9sa3aZdEutmjaGd0rln/jblXX6CWMoTYwbdZaftt8tMT6vq3q0yA8iACH0J9r9p3i8xX7mPP3oSpHr3MF9UIDjZGuKSTs1P7olCNF1897RH8/YblAslIh7WTF8rbb2OijGEN9BrNyz0lu+S2VjFznM2AM/OfvDOvQgHev6UNKZg4z5u9g0fbjHDjl2p4aVaFeWJC3JXiHnEz8cjO8rcJ1HF5fvvFd8Q74B8Ccf0B4lO2sDzQ7n+bn3F01fV7GGOozhISUTGoF+RMaFICIMPfvI3z+514yynHxLd5xgi6P/UqryFD2Jtgc9eVmGtcJ4d1r+iDAWU3qeFuOZ8jLg4Qd0LCTXn5/FMOP/A21noXQ+hB9pXf1VZV3h5ef5pcHCn+nlHwDLI2dHW6ieXOn0UOrDcZQnyGMnbGYuqGB5InQo1kE36+r2FBtXzHS/dvU5783D/K2DM+SekJ3Q/v1IThvBiQdgCO6Z02BW6A6G+oc9/QAAiDSxmCZaoAx1DWcEymZ5IlwIiWTEyn6D7H7eCUm+/QBzj6rEZMH+35IyiL88hCERsLwe0tuy0yB+CXQcSzk5cL8xyHzNEx4vTDN4fVFa5s/3e28nFe76QlZh9zlUvlVJuUYfDIBLvkQTh+E1kPh9CF9nKkJeoRgWd3vqsL0rRAcDstXuyd/D2IMdQ0lPSuXlMwcpn66irX7Er0tp8pc1LsZr1wW7V0RuTmAgH8FGqZWvKW/HQ11ZjLEPQ/Ht8LO+dBmBOxZWLi9dhNo0BG+uREadrZXTtJ++O0x3zDU676EyPaQegyObITjW+Bt6y0oOAIyHQxzSAQEVWwSB9vUaeKefL2AMdReZtfxFESgfSPXTl907+z1zPn7sEvz9BZbnx5HkC9MOvvWAEg+oo1hoy5w1nkl0yz6F/zxDPzfsaLrX2wLfSbD4pehTjNdu8zH0UgDLHyh8PfxCk6c8PWNOoRnI5sG3pWIgOTB97eUniazWO05I8l9NeoahDHUXmbiG0tJycwh/vnxVconN09YvfcUWTl5PPrDRvacqJ7uDUd+u2c4Gdl5hAR6edhvdgaBWYmFo+AWPKu/87uL7VkEterB7Ot1gx/Am/3hVHxhHmkJ2khDUSPtajZ+rQ3f1V+7r4ziHFwDRzfCyvfon1SBCWXdxVVfQ73W3lbhUoyh9hJz/z7MTxsOkZKZU+k8UjNzCA7wI8Dfj982H+GWz9a4UKH3uKt3MP16RdMhyk2vxHbJSoV9K2DlewzZ/nPJ7U9EwFkTYMuPJbc5GmlPk+PhbnvvFQ5E8Ykw/E16Qngjb6twKcZQe4nbPq+8Uf1u7QHqhQZx2+drGHVWFK9P6uVTfZurSrNwP4b6wmwsP98Paz8rO40zI+1tlIcmZDq6CZa/5ZmyKkJoKQGcqjF2Jrf9EDgPOCYi1XcMpo+QnpXLx8viS6zPyc0rMgKwLO75an3B7/+tP8Tp9GwCfcGHW0XuH9eJOiGBNEzf7V0hudnw1/twooIBf3yF0iLEpSdqH/JXV0PbWBhxX+XyP7UXPr8UslLc68apKBNehwadwK96RsgrCzs16o+BN4D/uFdKzeZ4ciYzN2Ty84lNfLVqf4ntGTl5hJdjbPefTOOQk/jNC7cfd5lOb3Bx7+YMaR/JhJ5NCfD3Iy6u/PghbmX1x/DLg97VUBUOr9NumVuXQ1SXwvUvttGGGmDvUmg9BOq3g/lPwNlPQO0oyE7X3efqNNUPLL8AUH66gbNZH/jqKqjdFJJdP2ValWl4FrTo520VbsHO5LaLlFKtPaClRiIivBW3iyNJGSw7lAOHShppgFd/286Ijg0Z3rFhiW0bDybx4/pDzNt0hHgfGXjiKi7q1YybhrXx/AjDzGRtlAKCtSE6sV0PS/7lQfd1F/MU6af0967fixpqKRab5aNzHBYExj0P306FHb/q7oEntutNteoV5gm+Z6Tz9UW287YSt6FEyo9sZhnqn8pyfSilpgJTAaKiovrMmjXLVRorRUpKCuHhru3yVlFEhOPpwv2L7PuPPx4XVmLd8yvT2XrS+wGQ3IGz4/XEtRuwYiq1Mo4iKEChqHnnd3/ziWSENCQoK4lmB/9HQE2KDeLA7jZXs7/FBUgpgZd8wRY4Upqe2NjY1SLidKy7yxoTRWQmMBOgb9++EhMT46qsK0VcXBze1tD50Z/p3bIeYN9QO2reeSyZp37awrZTvt3Vrl/renRuXIdPV+ylb6t6TB/TkQ+XxDN/S8l4DK9e3pOQAH8e/WEj1w9pQ0xM+xJp3HrtDq6G9V9BhtamEErEP64htDhQ/aacKmDg7dBqsHbBJB+GWaUPkW/bshltY0ufEMAXbIEjldFjen24gSd+3ISIkJGdx7JdCZXO56u/9rOoGvifZ98yGIDuzSIY2DaSlpGhNKodzImUTO46uwO//H2kwC9/Ya/mAJzT3Uujxj4YA3mV7xJp8BD12zgfUOSMzlUbg1AdMIbaDTjr1WGXHk/8ytUDW5GWlctxN01X5YyJ0U3JzM6jdkgAoUH+fLJ8b6lpgwP8nM6mclm/FgW/2zeqzfe3DwEgtlMjbo9tz9FkL756pyfq4dzGSPs2YQ0htAF0vbD8tJHt4c7qH8fDDna6530JxAANlFIHgMdF5AN3C6uOnErNIjE9u0p5nM7I4a24XS5SVD5BfpCVBy9c3KPICMAf1x+iTYMw1jiJE/LNrYNpUS+UlfEnOW3zeFtGhtIy0ovDIRa/DMte8175BntEdYNrv7eXNqxmDWopCzu9PiZ5Qkh1Izkjm/DgAJRSJGdkc+BUOi/+spUF23zfVeHIlB7BTJ0wosQw7bWPjSExLYuL315GnVqBrN2XyLndG7P/ZDptG4YRGhTA6C72g7d7lP9MhOb9YeQjcPowvN5bGwCD7xNis/fPgFth0O3u1eJDGNeHTWbM307HqNrc9vkabhrahveX7OH+cZ1oVT+Mv+JPVsnd4U0ahSoiQp23ltcNDeL3f8SQlpXDOwt3c/PwtoQF+/gtIwK74/Rn8UvQtJeeHPXASm8rM5RHhzEw+inn2y58F4LrwKxJ0H8qnPO8Z7V5GR//1/kOM+bvKPj9/hI9IOPFX7Z5S06lGN+jCWef1YjZqw7w2PldCAsKYNeG8g1YaFAA00d39IDCKvLB2MK+v6D7DR+sQT7MNiNg7zIdBjVpH6t7v0ifAcMKQ4hWZwJqwVWzS9/e8wr9/XiiR+T4GsZQl0P8iVTybPQ1rw68eWVvoLDnBYDnvOEeYP8KbyuoPH6BkFeKvz+qm37oXPEFBIVB4j5IP0Xy9kQ9oGXU4zqu85zpHpXsUtqPspfOU3FMfAxjqEshM0f3urh85nKOnnZf74sJfsvIxY85eQPdVsYXNw2gbmgNnARWRE9L9df7sHSGt9WUztDpOpqbs2HpkR10aNTpm/WEBGs/g1r14YfbCtO0i4UxzxQu12ulP9vj9PKw6XoCAl8loiUk7StcbjkYDq2FHGt8wc2LdQ8OQ6kYQ10K987ewP/Wu2+obHN1nCYk8FrQGwDMyXCfoR7c3gci0bmD35+EJa96W0VJxr8MgWGFAfTPfhy2z9O/67bS0d26XwILX4Qpf+gug6H19fbBd8JJK9ZJULgOfNT9svLLbGBNets2FnYvcO3xVJbe10Gf66BxTz3byy8Pwebv4Zpv9RvCc011uiY9vCqzOmAMdSm4z0gLCmFB0HQCVTlTgFeQy/o2JydX+HnjEdKzc1l0XyxZua4twycQgR2/wZafPFZkRnBDQupEwvhXtPth5ojSE/e6RscQkVxtmEHPFdhlIox+WteGofReC/Vaw7gXtDEPs/mQrdtCT2Sw9N/eN9QNOsIF72gDnD9tWZ2mcPEHMPFNCKyl1132nxoX4N9dGEONDjHq76c4cjqDurWCWH8g0eVlBJLDxf6LiCCVhwK/LLH92Qu78ch3GyuUZ2iQP2lZhYa4Q6PaTBnelnvHprPp0Gn7/ZYzUyDlqG8HtUk7CYGh8PtT2gguecWjxa8Y+C4xI2IKQ2he8YWOLjf7Ohh0h54oYPi9WmdAsE7T6+rCDIJCtWGyg1IwsIzprMoiN6ty+1WV6KshuDYr87rR/+yJEOwksJV/APg7xLjoMtFz+qo5Z7yhTkrLpt9z8xnavgF/bD1Gx6hwth9NcUneYaQTTDYrg29jZV5nBvtvLjVtk4iQCuW96L5YkjOzGf/aEvq3qc+JlMyCfs1N69aiad1aRXeY/yS0HKhnvM4nN4fap3fAly9B/GKIeUjHVeh7A4TULaz5eZMDq7VxnDkCGneHI397tvxJs6B2Y9ieVDTOcefxumZ/y1Jo7EN9tJv08mx5LQfrqb8ueBOAtLg450baUCVqvKHecvg0nRvXRjlpLd6bkMq+k2lk5eTxx1Y9GakrjPR5fstJoA5fBj1LfF4UASqvTCMNEOxfsdbs/Nryt7cNpnPj2oQGlXIpM07rV3DHGuh5MyBxLyTuo8/GbwrXx/1Tf6/+WH/nzwlYERJ26dd9/0reWjmZhKQfgZTjkJEI748s3OZJI935PKhVV/ft9fMvbLhzRCnfMtIAHc7WjXPfTqn4xLiVYeIbvv0mVkPwbUN9aB3UbVnY0GKD+BOp1LUGcCzbeYIr3/+T2E4NOZ2Rw5Rhbdh/Mp0BbesTGhTA2a8sLCe3Qn4MeoRfc/sxKzeWU9TGjzxySjl9bwS9XvC7tV/JCHLOqJ1xhBAyySC4zHQt64fSvVmEXkhPpHez2s6N4o7f4PB6+OPpktt+utuWJlITdLevsozuobU6uPzWuboG/t3NukFs8DR789blZsO8R2HINH29N33HwL//C+siSs5Y7UnajYR+N3qv/KrQpEfJ2NPOqN8OTlayg2b70bDzt0IfvMGt+KahPrFTNzjMHAGNuugGnKa9ILCke2DX8RTiT6TSsn4oEbUCiXkpjj6t6pGTlkHPpCMABcO6V+89VWL/8qhLMn4IPfz20MNvD/cF/peFuT0Y4b+Ba7IeZHFeYYt1FxVPhKpcSNIe3wxjVlA7LshyYlgd+OMfI/Bb9AKs2QQ/3qFb1ie8pgOnL35FP9TmP1EpDSX4V1sY9RgM+0fhutwc3YIfFKZr6zNjSu637HX9mfQVtBmu/bOlsfN3+PNt/XHEm0YadIjN6kx2KQGwGnWF7hfr0X0b/uu87/W5L2nX17c3Fa6r3xZiH9FvGEkHdP9tkTO2X7On8S1DnZ2hO/2/0adw3bHN8NE4ZOSjqDYjSky1c9nLP6KAE+haZnN1nPV7c8ghgK3HdzDGbz3z8io+PU8UJ2msTvJD8GMlto3w3wDA+X7LOSXhnO2/hrsDvq1wGcWJ9iu9dvPuNX1IO7qLgBNbYKHD8Nm1n2lDvfBFHR3O1Sz4J+z7U9cul7+p/Y9bf9LDeTNPl73vNzfq7mXTt+hW/4/O1a6RG+fprmd+/t79ow++U9cIF78MrYbAxq+1f37cCxBQzfud974WFjj0vR79lO6f3fuawnXRV0Jagn57WPE2HNsCxzZBn8m6t8aPd+jeKt0ugbYxUMcKTRtizRpjjLTH8C1DPTMGjm9xukn98TTwNGkXf0Zo9/PZsu8oa7bvZXXIrQDkimJhXk9G+q/j45wxPJEzmfsDvuLGgJ+5JPMxVknnCkn5NfgB6pZTO74sYCGXBdh3n1SGpyZ2pXFmPGPS5sCi6bCoWALJhS8u17Vbd5CXradm2vFr0fXlGWnQRhp0/9n0k3qePoB/99A1u2Ob9PRX3iJ/EEn/Kbp22OkcOOv86m+kQU9cO2w6vD9Ku6eG3FUyTWAtGHG//n3JB7pdIPlQYZe6hw/rb7/qP3Fydce3DHUpRtqRP/9cSpeghuT9926uyi1M76+Ekf7rALjMfyGd1AFqKz2/YKz/OvrIDlbmdaa+Os3veYU19kiSWB1yK3dl3YYC5uYNoIlKKNdIu4vYqHQWHQ0iF3+6qd1cuOIJaidtL3un7b94Rlxl2fx9yXXHNunvFHs+fJcy5Q/9au+IUrrfck3Czx+u/xlybMYBD2+oPwX7GwPtK/iWobZB7IG34cu3KaseFqoyGeTQy+L2gB+LbO+f8SZZBDDGf5U1Zx78O0i7DZ6VDwhTngvYX5yPkm7k84BRzG19PzftfaF8I22wT3iUfjA07X3mvLYH1iocYGKotlQ7Q+0KVoaUHsfWm0Y6n6sCfueyWsEEtq4FzictN1SGC9+BVkPPHCNtqDGckYa6OhC4Y663JdQcajcB/yBo3q9m+J8NZxy2nFBKqXFKqW1KqZ1KKSchwFxADQklavARRjwAPa2Zq9uMgLs3mBFzhmpLuYZaKeUPvAmcA3QBJimlurhciVLcmXWHy7M1nIH4BULswxDzoI5iN9jcV4bqjR3XR39gp4jsBlBKzQImAmWPia4E6eWMyjMYyuXaH/TgDNAjJR9xX6hag8FT2DHUzSjapHUAGOAOMWnGUBsqS5OekJOpB2YYDDUMO4baWRN5CYeyUmoqMBUgKiqKuLi4CovJENPQY6gYf3d7hPRaTUkLbQbkQSXuOzukpKRU6p52J0aTPXxNU2X02DHUB4AWDsvNgRLvkyIyE5gJ0LdvX4mJiamQEIAjyelQg+YiLZcpf0BwhB6YsOlbHWvZUCG6x17skehtcXFxVOaedidGkz18TVNl9Ngx1H8BHZRSbYCDwBXAlRVWZ4Pe7ZqeGYZ62D90L5dmDjFNvDmUujoRUld3t4uepGczqdvS24oMBrdTrqEWkRyl1B3Ar4A/8KGIbHKHmMBa4eUnqo4Un2F66HQILnasARWbOMDlTJoFKPjycj090ql47+opjf5TYeQj+rez+BUGQw3E1oAXEZkLuH0ERnBoDTPUTXrqmNBXfAEdx/DXTx/Tr97pkkYaPDPM96wJuub+13sFqzZ2fYhu/YbqKGkADx3UAftf7ep+PRWh1zXQ93rPz2BiMPgAPjUyMcTXa9RD7tLzC676oPy0d6zWNVOHUJ6p4a1hSIzz9O6uUUe2h8s/haw07dPtMAa2/cyJzK6FRhr0QyQoDGIe1qP4XBXbugIcbzCIhuEBWle91jqqXX5YVIPhDMSnwmPVCnFurA6IzZmY3c3opwonLK3dpPR0tepBg/Z6ZhS7cSUqWqMe9Tic92rRdVHdIchh9J3yAxTc8Cvc+JteFxQKA2/VxnrwHc71KQUxD0C9Nva0XPoJ3LZCl93ner0utPLXbGvnaTD5Jz1gpecVepYZY6QNZzA+VaMODih8bjyffQW/5PUjXprQVcXzQdC/qM9pglRuGTm4S1gdGGQFcmrWW88lmHQAtvwEvzyg15/9BGSl6iDroZEVL6OiNeph07WG8MZQvw3sWw5jnoZ2sfDtVD0V123Ldb616lZcD2hXycQ3tftm5czC9fft1t9vDdQTAnS9QC8/fEB/nz8DPhoPe5fYK6dpb90oOOIBOLGd3OM2Z083GM4QfMpQ+/np2t1pqcU7uRPw91PMv2cYDWuP4bMVYxmyahrRqUuL7lS3JSTuK1xuMQDqtiR113LC0g7AwNt1APs1n1Rc0J1rIDC0cGYLRyKaw8BbdED89bNg6D0Vz9+Rihjqxt0LNdy7Tc+Ms39F4WCPi2a6ZpokPz/9BtH1Qmg3Sj+kTu6GMOtBNH0zzrvZA7EP6anBjm6ClCOllxHVDS5+v7CLXVQXt/WFNhiqKz5lqAFGZL5CooTz9AXd6BRVm/aN9Kv87bHtyez3OX9tWEe/eRcAkPPgIQICAmDhC9D+bDixXRuVkAj2fvUUXba8rF/hQyLKNtT12ujZpLf8r+h6O/1zYx/Wn6pSu7H+rt9WG8PijH5K+5f7T9EPD0cCQ0qOyHNlKM+gMOg0Tv92nLA2fyYQZ7S2GiiTDkDKMXgvVq+/eZGeDPftwTqa3U3zXafTYKih+Jyh3ivaYF0zsOTsxsHh9eg3OJbMDX0JPrKKgJAwvWGUNa+hw4Skx6KG0+Vyh/kOr/4WstPgq6tLFtqkB1z2H/j2ZtgwC+7ZbH9WDFcRWh8eTYClM4rOHH7xB9Cgo9ZYHYlorj8XvK1DjTbpqddf8A50HOtdbQZDNcHnDPVLl/YkL6/skKfBU36BvJyKZdx+lP6+4VcIawiv94ZBd8DyN7R7BGDC69rP61hr9CT+AYW1+Ka99Fx3rQZrP3B1J7rYGKnoSd7RYTBUQ3zOUF/Sp3n5ifwDy37tLouWA/X3E0n6e+yzhdsCgrxnpPPpcgFMWwsRLfS0UTXBSBsMhirhc4b6jEepwjCdETYeWgaDocbjU/2oDQaDwVASY6gNBoPBxzGG2mAwGHwcY6gNBoPBxzGG2mAwGHwcY6gNBoPBxzGG2mAwGHwcJVL2KMBKZarUcWCvyzOuGA2AE17WUByjyR5Gkz2MJnv4mqbS9LQSkYbOdnCLofYFlFKrRKSvt3U4YjTZw2iyh9FkD1/TVBk9xvVhMBgMPo4x1AaDweDj1GRDPbP8JB7HaLKH0WQPo8kevqapwnpqrI/aYDAYago1uUZtMBgMNQJjqA0GQ7VGKVfOO+caXK2p2hvqM+EiuQKjyR5Gkz18TFMlZxFxKy7VVC0NtVKqq1IqBkB8xMluNNnDaLKH0VQ+SqlBSqnZwEtKqS5KKf+aqqlaNSYqpfyAN4CRwD7gT+AHEVmllPITkTyjyWgyms4ITY2Any1dLYBmwCoReU8ppbzxIHGnpupWo64H1AbOAq4CEoB/KKXCvXGz+LCmukC4D2rytfPki9fOaLJHT2C7iHwEvAx8C0xUSnUUEfGSa8ZtmnzeUCulLlFK3WYt1gEGAaEichz4BjgJ3G6l9cjFUUr1Vkp1tBYjgME+oKmtUirUWqwPDPEBTa2UUiEOmnzh2g1VSrW3Fuv6iCZzj5evZ5JS6kml1ARr1Vqgr1KqnYikAn8Bq4CbwTOuGaXUBUqph5VS461V69ylyWcNtVIqXCn1DXAvcEopFSAie4ClwN1WssPop1YvpVRTd18cpVQbpdQc4E3gU6XUaBHZDSzzoqYmSqlFwGfAD0qp7iKyE4gDpntJUxel1PfAx8D/rBrFTmAFXjpPlq5oYBEwSSlVR0R2Acu9pcnc47b0KKXULcD9QDzwL6XUTUAK8B/gLitpIjAfCFVKNXGXHktTQ+v+no5+YH2klLrE4SF2p6s1+ayhRvt4jorIQBH5Esi11n8MDFFKtRWRHOAokAHU8oCme4F1IjII+AG4wVr/oaWpjSc0FavBXA78JSKDgd+BB5VSvdHnaaCnzlO+JqVUZ+BtYIGIxALr0T478O55Au0z/A1934/wAU0+cY8X0+QT93g+1kNgEPC85VK4HYgBRgFzgfZKqbMtF0wC+honuVMT0A5YKiLDReQd4B/APda2L4HOSqlRrtTkU4ZaKdVYFbaS9gSaW+tvAx5XSg0FNqOf7i8BiMhGoBWQ6UZN+V1t0oFs63cdYIv1Gr0U/ZrjEU1ASLHfgVa5zwPH0DfxUXSjz4se1pQEPCgi/7aWn0bXKBoCK4E1XtCUzylgB5AH9FNK1UK/eXjr2vXAB+7xfE2WwU7Fy/e4UupapdQIpVR9a9UWoJn1xjEf2AgMBI4DXwAzLI2jAAUEuUlTjOVeXI2uzWPZq83WB+BvYBbwb1dq8glDrZQapZRajH7dettavQY4pJT6EP1ETQQeASYCrwCNlFJvKKU2omNfJ7nSV1ZM02vW6sXoJ/haYBzgj75RRgD/BJq4WdNopdRv6Ne/K6zVe4AEpVRLa/kroDvar/hPoLlS6nUPabpMRA6LyHKHMroDWSJyXERSgKc8qOlFh/OUr2UNOtZCCPAwcAn6PDX1oKZJ1uo1wGEv3uPFr50AS4AOnr7HLRdHE6XUAuA6dKPl60qpOsB+oBGQ37YwC+gGRIrIZ8DnwIPAFcD9IpJYVT2laLoSeA+IEJGjSil/EclFN7JGAIhInoh8jDbkrtMkIl79AB3RNb9L0BfjZ3RDWEt0zWs1EGilvQaYaf1uhG7gmOABTb8Aw6xt7YBvHdI+Crxh/Y5yo6b2lqaJQC/0n+cOoAnwEXAehd0tPwYes3439qCmz4CHrW3512xM/vlx2M+d186Zpv+ztp0PXA20RbtkTgPXeuE8fY5+XQ5A9w7wxj3u7H6619rWyZP3OOBvfXcEPrN+BwBvAZ+g3xg/tM5NhMM9/qxDHkEuPj+laXo9/9w4pPkPcFn+feQOTQF4AaX7ZSLahxMNrBSRr62nZyoQLyIHlVJzgb7ApegbaT1wsdJ9N4+hX/M9oSkF2KOUCrJ+71dKnSUiW4AFwN1KKSUiR9EuB3doGgCsFpEfrG2/o//kn6BdCsMsbXHAT+iHHSJyBDjiIU1/AK8opd63rg/oV78/re2PAh+JyAHcd+1K0/Q2+qE2BXgcmAPMQ7tl/L1wnl4GPkD7gXvh+Xu8uKb56PP0KbqBzO33uFIqAP2G5W/91+tg+elFJEcpdQe64bIL+txcgHYV/RPtvlqWn5eIZHlI0zT0m/4IEVlo7ZZvH54CLlJKjRORA67SBF5wfSilrgcOoH2XABuAPkqp99D+nUboV7H3RSQO+De6z+YD6FeeJVY+rnwFtKPpBfQT/gQQCUxTSt0FvItu2XUpTjT9je6t0NpaDkC3gr9gaTgAvKyUehCYgTbYntYUCOzC8mNa16gPuhFqIdAZ7Sf2tKY9aOP8NfpaDRKRu9G+zmS0D9GTmgIsTS+KyCL09fL0Pe7sPO22tieju1O67R5XSo1Av0nUA3Za5WYDsUqp/lDwQHkKeEG0b3omMFQp9ae1X5wXNIml6QlrH390g+vXaKMea1VEXIsrXxdsvE6EA9+ju9SsATpZ6xsC9wG3WsshaIM4yFoegO6LOMjLmhLQvrGWVvpPgIEe0NTZWj8D3aq8FP1K3x3d8t3Q2n4O8Bgw1Mua5qBbupWV9g+gl5c1/Qw0KLZ/oA9cu8bW9n4evMfLO09haL/rnW68x4cB1zgsvwXcCkxG1/RBVyQbA7OB1ta6ukAzV+uphKb/ohtT21nnsbc7NBVocWfmpZyMltb388CXDgf/AZYf2Fr3BnCej2l6EzjXC5q+sn77o2s6Q63lFmhfXYgPavJHP9zcewPb1/QREJx/bX1Eky9eu09wsb+3FD2hQDCFft6rgH9av9cBd1q/++b/J31M0yxPaMr/eNz1ISL7rJ8zgHZKqXNEv+LsBGYqpToppR4GhgKbfEzTEHRXIU9raqOUGiu6hTlJRJZY224B0ijsTuVLmpSIZIjIGh/RlA7kWPu4ddhzNb92qRT253annjQRybQ0AIxGd7cDuB44Syn1E7rG79Z7qJKaVoMHowh68qng5Al2M7DYYfkldIv4l0ALo6mIpoUOy/3RjVAFr85Gk9FUHTWha/V+aJdLe2tde7SLYyhucnNUN01ei55ntWrnKaW+Rrdsp6H9Pn+LSLrR5FTTYfQAg/nADtFDoI0mo6naarJqpEHA+8B36Ia5BLSb4bSn9fiqJq8NeLFullB0j4rLgH0istJbBrGaaJpkafrFW390o8locrEeQXdPvAodO+M7EbnOWwbRVzV5pR+1A7eh/U+jRcRdw2MritFkD6PJHkZT+RxAj8h8xUf0gI9p8urEAcpLQcfLwmiyh9FkD6PJ4Aqq1QwvBoPBcCbiE0GZDAaDwVA6xlAbDAaDj2MMtcFgMPg4xlAbDAaDj2MMtaHao5Sqq6zJYZVSTa3BHAZDjcH0+jBUe6xQnT+JSDdvazEY3IG3B7wYDK7geXQwrXXoORHPEpFuSqnJ6GDz/ujwtC+jhwZfgx46fa6InFRKtUNHRmyIDhswRUS2evogDIbSMK4PQ03gQWCXiESjY4g70g09111/4FkgTUR6AcuBa600M9FxHPqgZ+F+yxOiDQa7mBq1oaazQESSgWSlVBLwP2v930APpVQ4eg7A2Q4RK4M9L9NgKB1jqA01Hcc4DXkOy3no+98PSLRq4waDT2JcH4aaQDJQuzI7WhHR9iilLgUd4lIp1dOV4gyGqmIMtaHaIyIJwFKl1EbgX5XI4irgRqXUevSsQhNdqc9gqCqme57BYDD4OKZGbTAYDD6OMdQGg8Hg4xhDbTAYDD6OMdQGg8Hg4xhDbTAYDD6OMdQGg8Hg4xhDbTAYDD6OMdQGg8Hg4/w/uQhD5QjJa8IAAAAASUVORK5CYII=",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "(OHC_integrated-OHC_integrated.isel(time=slice(0,120)).mean('time')).plot(label = 'Time-integrated Ocean Heat Content')\n",
    "(hfds_integrated-hfds_integrated.isel(time=slice(0,120)).mean('time')).plot(label = 'Time-integrated hfds')\n",
    "plt.title('ACCESS-CM2 Heat Content [J]: 1850-1860 baseline', pad=20)\n",
    "plt.legend()\n",
    "plt.savefig('../EN4_Data/Figures/hfds_nonclosure.png')\n",
    "plt.grid()\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3.7.6 ('base')",
   "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 (default, Jan  8 2020, 13:42:34) \n[Clang 4.0.1 (tags/RELEASE_401/final)]"
  },
  "orig_nbformat": 4,
  "vscode": {
   "interpreter": {
    "hash": "00ec6b5807eb4cae4bdd23e758d78199b5749ea4aa251c45b926838befc3bdc7"
   }
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}