Skip to content

Instantly share code, notes, and snippets.

@peterm790
Created May 24, 2023 16:25
Show Gist options
  • Save peterm790/d4ed17800d1175dda273850f409394c5 to your computer and use it in GitHub Desktop.
Save peterm790/d4ed17800d1175dda273850f409394c5 to your computer and use it in GitHub Desktop.
identify_exceeding_26
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"id": "880f73ec-cf80-4ca5-ab87-420c231064d6",
"metadata": {},
"outputs": [],
"source": [
"import intake\n",
"from pathlib import Path\n",
"import xarray as xr\n",
"import dask\n",
"import numpy as np\n",
"import pandas as pd\n",
"import matplotlib.pyplot as plt"
]
},
{
"cell_type": "code",
"execution_count": 2,
"id": "6937e239-e10a-40d3-88e1-58800fe3b26d",
"metadata": {},
"outputs": [],
"source": [
"catalog = intake.open_catalog(Path(Path.home(),'heat_center/data/climate/reanalysis/intake/reanalysis.yaml'))"
]
},
{
"cell_type": "code",
"execution_count": 3,
"id": "1da6580a-8337-4a0e-b499-a284daefc54d",
"metadata": {},
"outputs": [],
"source": [
"lat_, lon_ = -26.198513, 28.029066 #JHB\n",
"#lat_, lon_ = 5.383200, -4.038460 #ABJ"
]
},
{
"cell_type": "code",
"execution_count": 4,
"id": "9b04b0a6-6e5a-4a9e-9e49-2fbf435b490e",
"metadata": {},
"outputs": [],
"source": [
"variables = ['hus','va','ua']"
]
},
{
"cell_type": "code",
"execution_count": 5,
"id": "2284d184-7a73-4e43-8bd9-74d767cb53f6",
"metadata": {},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"ERROR 1: PROJ: proj_create_from_database: Open of /share/apps/anaconda3-2022.05/envs/pangeo/share/proj failed\n"
]
}
],
"source": [
"datasets_to_merge = []\n",
"for var in variables:\n",
" ds = catalog['ERA5']['day'][var].to_dask()\n",
" ds = ds.sel(latitude = lat_,longitude = lon_, method = 'nearest')\n",
" #ds = ds.sel(time = slice('2000-01-01', '2001-12-31')) \n",
" ds = ds[[var]] #just to get rid of time_bnds if present\n",
" datasets_to_merge.append(ds)"
]
},
{
"cell_type": "code",
"execution_count": 6,
"id": "ac6a4a2a-61d6-4b80-b513-1e8ddf279304",
"metadata": {},
"outputs": [],
"source": [
"ds = xr.merge(datasets_to_merge)\n",
"#ds = ds.load() #load into memory to speed things up later"
]
},
{
"cell_type": "code",
"execution_count": 7,
"id": "c39f18f3-a15b-446d-bd7d-9e36c6b81ca8",
"metadata": {},
"outputs": [],
"source": [
"ds_HEAT = catalog['ERA5-HEAT']['day']['UTCImax'].to_dask()\n",
"ds_HEAT = ds_HEAT.sel(lat = lat_,lon = lon_, method = 'nearest')\n",
"#ds_HEAT = ds_HEAT.sel(time = slice('2000-01-01', '2001-12-31'))\n",
"ds_HEAT = ds_HEAT.UTCImax - 273.15\n",
"ds_HEAT = ds_HEAT.load() #load into memory to speed things up later"
]
},
{
"cell_type": "code",
"execution_count": 8,
"id": "feff923b-f5a3-4a3e-817b-07b415dda06f",
"metadata": {},
"outputs": [],
"source": [
"ds_HEAT = ds_HEAT.where(ds_HEAT >= 26).dropna(dim = 'time')"
]
},
{
"cell_type": "code",
"execution_count": 9,
"id": "5370eb46-fa4d-4733-a644-146fc796ddb3",
"metadata": {},
"outputs": [],
"source": [
"ds_above_26 = ds.sel(time = [y for x,y in zip(ds.time.dt.date.values, ds.time.values) if x in ds_HEAT.time.dt.date.values])"
]
},
{
"cell_type": "code",
"execution_count": 10,
"id": "2d6089c7-f706-4730-bcef-2634f49ebfc0",
"metadata": {},
"outputs": [],
"source": [
"import seaborn as sns"
]
},
{
"cell_type": "code",
"execution_count": 11,
"id": "0ecf94f9-e49c-45a3-ae0b-18e9cf6cad80",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.legend.Legend at 0x7efdaf7af220>"
]
},
"execution_count": 11,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAjcAAAGdCAYAAADuR1K7AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/P9b71AAAACXBIWXMAAA9hAAAPYQGoP6dpAABn5klEQVR4nO3dd3yV9d3/8dfZ2XtDgLDDxgAKiIoiioqzStUqKlq5cVPtD+pdV63UqtxoVYqKIq0ite4WRRwMxcWUvSEhJITsnZNzzvX745BATMIISU5yeD8fj+sRcp3rnPM5J8B55ztNhmEYiIiIiPgJs68LEBEREWlOCjciIiLiVxRuRERExK8o3IiIiIhfUbgRERERv6JwIyIiIn5F4UZERET8isKNiIiI+BWrrwtobR6PhwMHDhAaGorJZPJ1OSIiInICDMOgpKSEpKQkzOZjt82cduHmwIEDJCcn+7oMERERaYKMjAw6dux4zGtOu3ATGhoKeN+csLAwH1cjIiIiJ6K4uJjk5OTaz/FjOe3CTU1XVFhYmMKNiIhIO3MiQ0o0oFhERET8isKNiIiI+BWFGxEREfErp92YGxER8V9ut5vq6mpflyFNZLPZsFgsp/w4CjciIuIXSktL2b9/P4Zh+LoUaSKTyUTHjh0JCQk5pcdRuBERkXbP7Xazf/9+goKCiI2N1SKt7ZBhGBw6dIj9+/fTo0ePU2rBUbgREZF2r7q6GsMwiI2NJTAw0NflSBPFxsayd+9eqqurTyncaECxiIj4DbXYtG/N9fNTuBERERG/onAjIiIifkVjbkRExG9dduU15OTmt9rzxcVE8Z8P32u2x9u7dy8pKSmsXbuWQYMGsXTpUkaPHk1BQQERERHN9jz+RuFGRET8Vk5uPtc/OqfVnm/B43e22nNJ49QtJSIiIn5F4UZERMSHPvvsM84++2wiIiKIjo7msssuY9euXb4uq11Tt9Rp5lj9z83dVywiIsdXVlbG1KlT6d+/P2VlZTzyyCNcddVVrFu3zteltVsKN6eZY/U/q69YRKT1XXPNNXW+nzt3LnFxcWzevPmUtyE4XalbSkRExId27drFDTfcQNeuXQkLCyMlJQWA9PR0H1fWfqnlRkRExIfGjx9PcnIyr776KklJSXg8Hvr164fT6fR1ae2Wwo2IiIiP5OXlsWXLFubMmcOoUaMA+Oabb3xcVfuncCMiIuIjkZGRREdH88orr5CYmEh6ejrTpk3zdVntnsKNiIj4rbiYqFadLBEXE3VS15vNZt555x3uvfde+vXrR69evXjhhRc477zzWqbA04TCjYiI+K32sLzFmDFj2Lx5c51zhmE0+OfzzjuvzvfSMM2WEhEREb+icCMiIiJ+ReFGRERE/IrCjYiIiPgVhRsRERHxKwo3IiIi4lcUbkRERMSvKNyIiIiIX9EifiIi4rfcbnerLnpnMpmwWCzN9nh79+4lJSWFtWvXMmjQoGZ7XH+nlhsREfFLbreblM6dsNlsrXakdO6E2+329UtvFa+++iqjRo0iMjKSyMhIxowZw48//ljvuszMTH7zm98QHR1NUFAQgwYNYvXq1S1am1puRETELxmGQUbmASo+exyrpeV/l3e5PQRe/Ohpsz3C0qVLuf766xkxYgQBAQH89a9/ZezYsWzatIkOHToAUFBQwMiRIxk9ejSffvopcXFx7Nq1i4iIiBatTS03IiLi16wWM1aLpRWOk/9I/eyzzzj77LOJiIggOjqayy67jF27dtW7buvWrbUhom/fvixdurTO7cuWLWPYsGE4HA4SExOZNm0aLpcLgDlz5tChQwc8Hk+d+1x++eVMnDix9vtPPvmEtLQ0AgIC6Nq1K48//njtYzTkrbfeYsqUKQwaNIjevXvz6quv4vF4+PLLL2uvefrpp0lOTuaNN95g2LBhdOnShQsuuIBu3bqd9Ht1MhRuREREfKSsrIypU6fy008/8eWXX2I2m7nqqqvqBZGHHnqI3/3ud6xdu5YRI0Zw+eWXk5eXB3i7fS655BKGDh3K+vXrmT17NnPnzuXJJ58E4NprryU3N5evv/669vEKCgpYvHgxN954IwCLFy/mN7/5Dffeey+bN29mzpw5zJs3jz//+c8n/FrKy8uprq4mKurIzugff/wxQ4YM4dprryUuLo7Bgwfz6quvNvn9OlEKNyIiIj5yzTXXcPXVV9OjRw8GDRrE3Llz2bBhQ71dwu+++26uueYaUlNTmT17NuHh4cydOxeAl19+meTkZF588UV69+7NlVdeyeOPP85zzz2Hx+MhKiqKiy++mLfffrv28d59912ioqK44IILAPjzn//MtGnTmDhxIl27duXCCy/kT3/6E3PmzDnh1zJt2jQ6dOjAmDFjas/t3r2b2bNn06NHDxYvXszkyZO59957mT9//qm8bcelcCMiIuIju3bt4oYbbqBr166EhYWRkpICQHp6ep3rhg8fXvtnq9XKkCFD2LJlCwBbtmxh+PDhmEym2mtGjhxJaWkp+/fvB+DGG2/kvffeo6qqCvB2Kf3617+undm1evVqnnjiCUJCQmqPO+64g6ysLMrLy4/7Ov7617+yYMEC3n//fQICAmrPezwezjjjDJ566ikGDx7MnXfeyR133MHs2bOb8nadMA0oFhER8ZHx48eTnJzMq6++SlJSEh6Ph379+uF0Oo9735owYxhGnWBTc+7oa8aPH4/H4+G///0vQ4cOZcWKFcycObP2eo/Hw+OPP87VV19d73mODisNefbZZ3nqqaf44osvGDBgQJ3bEhMT6dOnT51zqampvPfee8d9fadC4UZERMQH8vLy2LJlC3PmzGHUqFEAfPPNNw1e+/3333POOecA4HK5WL16NXfffTcAffr04b333qsTclauXEloaGjtrKXAwECuvvpq3nrrLXbu3EnPnj1JS0urffwzzjiDbdu20b1795N6Dc888wxPPvkkixcvZsiQIfVuHzlyJNu2batzbvv27XTu3PmknudkKdyIiIj4QGRkJNHR0bzyyiskJiaSnp7OtGnTGrz2pZdeokePHqSmpvJ///d/FBQUcNtttwEwZcoUZs2axT333MPdd9/Ntm3bePTRR5k6dSpm85HRJzfeeCPjx49n06ZN/OY3v6nz+I888giXXXYZycnJXHvttZjNZn7++Wc2bNhQOzD5l/7617/yxz/+kbfffpsuXbqQnZ0NUNutBfDAAw8wYsQInnrqKa677jp+/PFHXnnlFV555ZVTfv+ORWNuRETEr7ncHlxudyscnuMXcxSz2cw777zD6tWr6devHw888ADPPPNMg9f+5S9/4emnn2bgwIGsWLGCjz76iJiYGAA6dOjAokWL+PHHHxk4cCCTJ09m0qRJ/O///m+dxzj//POJiopi27Zt3HDDDXVuu+iii/jPf/7DkiVLGDp0KGeddRYzZ848ZgvLyy+/jNPp5Fe/+hWJiYm1x7PPPlt7zdChQ/nggw9YsGAB/fr1409/+hOzZs2qnaXVUkzG6bLa0GHFxcWEh4dTVFREWFiYr8tpdcPOHs31jzY8+n3B43fy4zdfN3ibiEhbVllZyZ49e0hJSakdI1KzQnFG5oFWqyO5QxJ79qU36xYMp5OGfo41TubzW91SIiLilywWC3v2pbfrvaWkaRRuRETEbylonJ58OuZm+fLljB8/nqSkJEwmEx9++OExr3///fe58MILiY2NJSwsjOHDh7N48eLWKVZERETaBZ+Gm7KyMgYOHMiLL754QtcvX76cCy+8kEWLFrF69WpGjx7N+PHjWbt2bQtXKiIiIu2FT7ulxo0bx7hx4074+lmzZtX5/qmnnuKjjz7ik08+YfDgwc1cnYiIiLRH7XrMjcfjoaSkpM4mXb9UVVVVu9w0eEdbi4iIfzrNJgD7neb6+bXrdW6ee+45ysrKuO666xq9ZsaMGYSHh9ceycnJrVihiIi0hpqBwyeybYG0XTU/v1MdCN5uW24WLFjAY489xkcffURcXFyj102fPp2pU6fWfl9cXKyAIyLiZ6xWK0FBQRw6dAibzVZnZV5pHzweD4cOHSIoKAir9dTiSbsMNwsXLmTSpEm8++67dbZWb4jD4cDhcLRSZSIi4gsmk4nExET27NnDvn37fF2ONJHZbKZTp071NgI9We0u3CxYsIDbbruNBQsWcOmll/q6HBERaSPsdjs9evRQ11Q7Zrfbm6XVzafhprS0lJ07d9Z+v2fPHtatW0dUVBSdOnVi+vTpZGZmMn/+fMAbbG6++Waef/55zjrrrNpNugIDAwkPD/fJaxARkbbDbDbXW7ZfTj8+7ZRctWoVgwcPrp3GPXXqVAYPHswjjzwCQFZWFunp6bXXz5kzB5fLxV133VVnk6777rvPJ/WLiIhI2+PTlpvzzjvvmNO+5s2bV+f7pUuXtmxBIiIi0u5pOLmIiIj4FYUbERER8SsKNyIiIuJXFG5ERETEryjciIiIiF9RuBERERG/onAjIiIifkXhRkRERPyKwo2IiIj4FYUbERER8SsKNyIiIuJXFG5ERETEryjciIiIiF9RuBERERG/onAjIiIifkXhRkRERPyKwo2IiIj4FYUbERER8SsKNyIiIuJXFG5ERETEryjciIiIiF9RuBERERG/onAjIiIifkXhRkRERPyKwo2IiIj4FYUbERER8SsKNyIiIuJXFG5ERETEryjciIiIiF9RuBERERG/onAjIiIifkXhRkRERPyKwo2IiIj4FYUbERER8SsKNyIiIuJXFG5ERETEryjciIiIiF9RuBERERG/onAjIiIifkXhRkRERPyKwo2IiIj4FYUbERER8Ss+DTfLly9n/PjxJCUlYTKZ+PDDD497n2XLlpGWlkZAQABdu3bl73//e8sXKiIiIu2GT8NNWVkZAwcO5MUXXzyh6/fs2cMll1zCqFGjWLt2LX/4wx+49957ee+991q4UhEREWkvrL588nHjxjFu3LgTvv7vf/87nTp1YtasWQCkpqayatUqnn32Wa655poWqlJERETak3Y15ua7775j7Nixdc5ddNFFrFq1iurq6gbvU1VVRXFxcZ1DRERE/Fe7CjfZ2dnEx8fXORcfH4/L5SI3N7fB+8yYMYPw8PDaIzk5uTVKFRERER/xabdUU5hMpjrfG4bR4Pka06dPZ+rUqbXfFxcXK+A0YseO7Qw7e3SDt8XFRPGfDzW2SURE2r52FW4SEhLIzs6ucy4nJwer1Up0dHSD93E4HDgcjtYor91zeQyuf3ROg7ctePzOVq5GRESkadpVt9Tw4cNZsmRJnXOff/45Q4YMwWaz+agqERERaUt8Gm5KS0tZt24d69atA7xTvdetW0d6ejrg7VK6+eaba6+fPHky+/btY+rUqWzZsoXXX3+duXPn8uCDD/qifBEREWmDfNottWrVKkaPPjLGo2ZszMSJE5k3bx5ZWVm1QQcgJSWFRYsW8cADD/DSSy+RlJTECy+8oGngIiIiUsun4ea8886rHRDckHnz5tU7d+6557JmzZoWrEpERETas3Y15kZERETkeBRuRERExK8o3IiIiIhfUbgRERERv6JwIyIiIn5F4UZERET8isKNiIiI+BWFGxEREfErCjciIiLiVxRuRERExK8o3IiIiIhfUbgRERERv6JwIyIiIn5F4UZERET8isKNiIiI+BWFGxEREfErCjciIiLiVxRuRERExK8o3IiIiIhfUbgRERERv6JwIyIiIn5F4UZERET8isKNiIiI+BWFGxEREfErCjdyhMmMYRi+rkJEROSUKNwIAE6Xh/Ar/pfXv91LZkGFr8sRERFpMoUbAWDVvnys0cmUVrl4f+1+1mcU+rokERGRJlG4EQrLnazZVwhAYngAHgOWbj/Evrwy3xYmIiLSBAo3wvIdubgNA2fmJq5N60hqYigAW7NLfFyZiIjIyVO4Oc1lF1WyJ7cMswnKv3sHk8lE/w7hAOw6VIrL7fFxhSIiIidH4eY0l1FQDkBKTDDuomwAEsICCA2wUu022JOrrikREWlfFG5Oc1lFlQAkRQTWnjOZTPSM93ZNbT9Y6pO6REREmkrh5jRmGAZZhd5p30nhgXVu63U43OzJK6PK5W712kRERJpK4eY0VlBeTaXLg8VsIjbUUee2mBA7kUE23B6D3YfUNSUiIu2Hws1p7ECRt9UmISwAi9lU5zaTyUSPOG/rzb688lavTUREpKkUbk5jWYXe8TaJ4QEN3p4U4T2fVaQVi0VEpP1QuDmN1YSWxIiGw03C4dBTXOnCYwtqtbpEREROhcLNaarC6aagvBqoP5i4hsNqITrYDkB1aFKr1SYiInIqFG5OUzWtNlFBdgJslkavq2m9qQ5TuBERkfZB4eY0lV18eLxNI11SNWrG41SHJrZ4TSIiIs1B4eY0VdMlVdPt1JjEw11W1SEJVGsrBhERaQcUbk5TheVOACKCjh1uIoNsOKxmsNjYklXcGqWJiIicEp+Hm5dffpmUlBQCAgJIS0tjxYoVx7z+rbfeYuDAgQQFBZGYmMitt95KXl5eK1XrHwzDoKjC23ITEWg75rUmk6l23M2afQUtXpuIiMip8mm4WbhwIffffz8PP/wwa9euZdSoUYwbN4709PQGr//mm2+4+eabmTRpEps2beLdd9/lp59+4vbbb2/lytu3MqebareBCQg7TrgBSAw7HG7SC1u2MBERkWbg03Azc+ZMJk2axO23305qaiqzZs0iOTmZ2bNnN3j9999/T5cuXbj33ntJSUnh7LPP5s4772TVqlWtXHn7VtMlFRZoq7cycUNqWm5+3l/YkmWJiIg0C5+FG6fTyerVqxk7dmyd82PHjmXlypUN3mfEiBHs37+fRYsWYRgGBw8e5N///jeXXnppo89TVVVFcXFxneN0V1h+Yl1SNWr2ndqXX05ZlavF6hIREWkOPgs3ubm5uN1u4uPj65yPj48nOzu7wfuMGDGCt956iwkTJmC320lISCAiIoK//e1vjT7PjBkzCA8Prz2Sk5Ob9XW0R4U1422CTizcBNmtmJ2lGAZsO1jSkqWJiIicMp8PKDaZ6naLGIZR71yNzZs3c++99/LII4+wevVqPvvsM/bs2cPkyZMbffzp06dTVFRUe2RkZDRr/e3Ric6UOpq17BAAW7MUbkREpG2z+uqJY2JisFgs9VppcnJy6rXm1JgxYwYjR47koYceAmDAgAEEBwczatQonnzySRIT6y8053A4cDgczf8C2rGT7ZYCb7hxRqZoOriIiLR5Pmu5sdvtpKWlsWTJkjrnlyxZwogRIxq8T3l5OWZz3ZItFu/WAYZhtEyhfsbg5Lul4EjLjcKNiIi0dU0KN3v27GmWJ586dSqvvfYar7/+Olu2bOGBBx4gPT29tptp+vTp3HzzzbXXjx8/nvfff5/Zs2eze/duvv32W+69916GDRtGUpL2PjoRHnsobo+B2QRhAScfbrZmlyhIiohIm9akbqnu3btzzjnnMGnSJH71q18REHDs/YkaM2HCBPLy8njiiSfIysqiX79+LFq0iM6dOwOQlZVVZ82bW265hZKSEl588UV+97vfERERwfnnn8/TTz/dpOc/HbkDIwBvsDGfwDTwGpaKfGwWE6VVLvYXVJAcFdRCFYqIiJyaJrXcrF+/nsGDB/O73/2OhIQE7rzzTn788ccmFTBlyhT27t1LVVUVq1ev5pxzzqm9bd68eSxdurTO9ffccw+bNm2ivLycAwcO8M9//pMOHTo06blPR+7ASODkuqQATIaH7nGhgLqmRESkbWtSuOnXrx8zZ84kMzOTN954g+zsbM4++2z69u3LzJkzOXToUHPXKc3EFVATbk58plSN1MSacKMZUyIi0nad0oBiq9XKVVddxb/+9S+efvppdu3axYMPPkjHjh25+eabycrKaq46pZnUttycxEypGqkJYQBszVbLjYiItF2nFG5WrVrFlClTSExMZObMmTz44IPs2rWLr776iszMTK644ormqlOaiTsgHIDwpoSbRG+4UbeUiIi0ZU0aUDxz5kzeeOMNtm3bxiWXXML8+fO55JJLaqdpp6SkMGfOHHr37t2sxcqp89hDAAgJOPkffe/D3VI12zAEO3y2TJKIiEijmvTpNHv2bG677TZuvfVWEhISGrymU6dOzJ0795SKk+ZV7nRh2AIBCG1CuIkJcRATYie31MnOnFIGJkc0c4UiIiKnrknhZsmSJXTq1KnegnqGYZCRkUGnTp2w2+1MnDixWYqU5pFdVAmAzWLCbmlaj2SPuFByS/PYfrBE4UZERNqkJn3CdevWjdzc3Hrn8/PzSUlJOeWipGXUhJsQh7XR/buOp2e8t1truzbQFBGRNqpJ4aaxFWpLS0ubvKCftLwDNeGmCV1SNXomeMfdbD9Y2iw1iYiINLeT+pSbOnUq4N3J+5FHHiEo6MgqtW63mx9++IFBgwY1a4HSfLKLKgBvy01T9Yz3hpsdarkREZE26qQ+5dauXQt4W242bNiA3X5kITi73c7AgQN58MEHm7dCaTZZh1tuQh0nPw28Rs/DqxQfKKqkpLKa0JPYn0pERKQ1nFS4+frrrwG49dZbef755wkLC2uRoqRlZDdDt1R4kI34MAcHi6vYfrCUtM6RzVWeiIhIs2jSmJs33nhDwaYdyjpqQPGpUNeUiIi0ZSf8KXf11Vczb948wsLCuPrqq4957fvvv3/KhUnzyy5unnDTIy6UFTtyNahYRETapBP+lAsPD6+dPhweHt5iBUnLqKx2k1/mBJq2gN/RNB1cRETashP+lHvjjTca/LO0DwcPt9rgrsZhPaUtxY6aDq5wIyIibU+TPuUqKiooLy+v/X7fvn3MmjWLzz//vNkKk+Z1oNAbbixVJU1ewK9Gjzhvy01OSRVF5dWnXJuIiEhzalK4ueKKK5g/fz4AhYWFDBs2jOeee44rrriC2bNnN2uB0jyyi71r3Jidp97aEhpgIyncu1jj9hy13oiISNvSpHCzZs0aRo0aBcC///1vEhIS2LdvH/Pnz+eFF15o1gKledTMlLJUNc8g4JquqW3ZCjciItK2NCnclJeXExrq/XD7/PPPufrqqzGbzZx11lns27evWQuU5lGzxo25qrhZHk/TwUVEpK1qUrjp3r07H374IRkZGSxevJixY8cCkJOTo/Vv2qjalhtn87Tc1Iy70XRwERFpa5oUbh555BEefPBBunTpwplnnsnw4cMBbyvO4MGDm7VAaR5HWm6ap6WltuVGY25ERKSNadKCJ7/61a84++yzycrKYuDAgbXnL7jgAq666qpmK06aT03LTXMMKAbocXitm9xSJ3mlVUSHOJrlcUVERE5Vk1dzS0hIICEhoc65YcOGnXJB0vycLg+5pVVA8w0oDrJbSY4KJCO/gu0HSxmucCMiIm1Ek8JNWVkZf/nLX/jyyy/JycnB4/HUuX337t3NUpw0j5pgY7OYMLkqmu1xe8aFkpFfwY6cEoZ3i262xxURETkVTQo3t99+O8uWLeOmm24iMTHxlBeFk5Z1qMQbbmJDHLia8XF7JoTy5dYcrVQsIiJtSpPCzaeffsp///tfRo4c2dz1SAuoDTehDrKa8XFr95jK1owpERFpO5oUbiIjI4mKimruWqSFHCo99XCzY8d2hp09us656uA4GHwzP+3I5NIrr+G/H753ipWKiIicuiaFmz/96U888sgjvPnmmwQFBTV3TdLMjm65aSqXx+D6R+fUPef28PLSXRi2wNqp5iIiIr7WpHDz3HPPsWvXLuLj4+nSpQs2m63O7WvWrGmW4qR5HD3mpjlZLWbCA20UVlTjDopp1scWERFpqiaFmyuvvLKZy5CW1BwtN42JDrFTWFGNK1jhRkRE2oYmhZtHH320ueuQFnT0mJvmFhVsZ9ehMlxBmgouIiJtQ5O2XwAoLCzktddeY/r06eTn5wPe7qjMzMxmK06aR06JdzxMi7TcBHsfU91SIiLSVjSp5ebnn39mzJgxhIeHs3fvXu644w6ioqL44IMP2LdvH/Pnz2/uOqWJDMM4asxNQLM/fnSIHQBXUAyGYWjNIxER8bkmtdxMnTqVW265hR07dhAQcOQDc9y4cSxfvrzZipNTV1rlorLau4J0TKi92R8/IsiGyQSG1UF2sWZMiYiI7zUp3Pz000/ceeed9c536NCB7OzsUy5Kmk9Nq02Iw0qQvclbiTXKajYTEeidLbf9oBbzExER32tSuAkICKC4uLje+W3bthEbG3vKRUnzqQk3cS0w3qZGzY7gO7QNg4iItAFNCjdXXHEFTzzxBNXV1QCYTCbS09OZNm0a11xzTbMWKKemZqZUTEuGm2Bvd9e2bIUbERHxvSaFm2effZZDhw4RFxdHRUUF5557Lt27dyc0NJQ///nPzV2jnIKWXOOmRk242Z6jbikREfG9Jg3CCAsL45tvvuHrr79m9erVeDwezjjjDMaMGdPc9ckpaqnViY8WdTjc7DxYohlTIiLicycdbjweD/PmzeP9999n7969mEwmUlJSSEhI0AdbG9QaLTcRQXbwuClzQmZhBR0jtd+YiIj4zkl1SxmGweWXX87tt99OZmYm/fv3p2/fvuzbt49bbrmFq666qqXqlCbKaYVwYzGbsFR4F3LcrkHFIiLiYyfVcjNv3jyWL1/Ol19+yejRo+vc9tVXX3HllVcyf/58br755mYtUpquNVpuAKzlebiDY9l+sJTze8e36HOJiIgcy0m13CxYsIA//OEP9YINwPnnn8+0adN46623TqqAl19+mZSUFAICAkhLS2PFihXHvL6qqoqHH36Yzp0743A46NatG6+//vpJPefppHZfqRYccwNgLc8F1HIjIiK+d1Lh5ueff+biiy9u9PZx48axfv36E368hQsXcv/99/Pwww+zdu1aRo0axbhx40hPT2/0Ptdddx1ffvklc+fOZdu2bSxYsIDevXufzMs4bbg9BnmlLb/ODRwJN5oOLiIivnZS3VL5+fnExzfe5RAfH09BQcEJP97MmTOZNGkSt99+OwCzZs1i8eLFzJ49mxkzZtS7/rPPPmPZsmXs3r2bqKgoALp06XIyL+G0kl/mxGOAyXRkRlNLsZYeAmDHwVKq3R5slibvySoiInJKTuoTyO12Y7U2nocsFgsul+uEHsvpdLJ69WrGjh1b5/zYsWNZuXJlg/f5+OOPGTJkCH/961/p0KEDPXv25MEHH6SioqLR56mqqqK4uLjOcbqoGW8THezA2sJhw1xVRKjDitPtYafWuxERER86qZYbwzC45ZZbcDga7uKoqqo64cfKzc3F7XbXawmKj49vdH+q3bt388033xAQEMAHH3xAbm4uU6ZMIT8/v9FxNzNmzODxxx8/4br8Se3qxCEt22oDYAJSk8L4cU8+mw4Uk5oY1uLPKSIi0pCT+nV+4sSJxMXFER4e3uARFxd30jOlfrkuzrHWyvF4PJhMJt566y2GDRvGJZdcwsyZM5k3b16jrTfTp0+nqKio9sjIyDip+tqz3FaaKVWjb5I30Gw6UNQqzyciItKQk2q5eeONN5rtiWNiYrBYLPVaaXJychod15OYmEiHDh0IDw+vPZeamophGOzfv58ePXrUu4/D4Wi0pcnf5ZXVtNy0Vrjx/lw2Hzh9uv5ERKTt8dmoT7vdTlpaGkuWLKlzfsmSJYwYMaLB+4wcOZIDBw5QWnpkTMf27dsxm8107NixRettj3JLncCRvZ9aWp/DXVGbs4oxDKNVnlNEROSXfDqlZerUqbz22mu8/vrrbNmyhQceeID09HQmT54MeLuUju7muuGGG4iOjubWW29l8+bNLF++nIceeojbbruNwMBAX72MNiv38Jib6OO03HjcbtxuF4Zh4Ha7GjwMw8Djdh/zcXrEh2C3mCmpdJGR3/ggbxERkZbUpI0zm8uECRPIy8vjiSeeICsri379+rFo0SI6d+4MQFZWVp01b0JCQliyZAn33HMPQ4YMITo6muuuu44nn3zSVy+hTatpuTnWgGKP281TN51Hfm4OAA+N69votU/ddB5/+MdSzBZLg7fbLGZ6JoSwMbOYTQeK6BStPaZERKT1+TTcAEyZMoUpU6Y0eNu8efPqnevdu3e9riyp67IrryEnN5/8QTdBSDxPPv4ozxTsAWDnrp11rjUwyM/N4evZ97Hiw39w3tUTG3zMr957k8c+z8Hg2N1NfRPD2ZhZzOasYsb1T2yeFyQiInISfB5upPnl5OZz/aNzeO2b3biq3Iy75X7iwwIAeOSG8xq8j8VsxmI2NboejvkEN3vvUztjSoOKRUTEN7SMrJ8yDIMKp3eMTJC94W6klqDp4CIi4msKN36qyuXBc7gHKbAVw01qYhgmExwsriKnpLLVnldERKSGwo2fKj/camO3mrGaW+/HHOyw0j02BID1GWq9ERGR1qdw46d80SVVY1ByBADrMk58E1UREZHmonDjp8qd3g1Mg2w+CDedIgBYl1HY6s8tIiKicOOnymtbblp/QlxNy83PGUV4PFqpWEREWpfCjZ+qCTetOZi4Rq/4UAJtFkqqXOzOLT3+HURERJqRwo2fKq8+3C3lg3BjtZjp38G7ieba9MJWf34RETm9Kdz4KV8OKAaNuxEREd9RuPFTvhxzA0fPmCr0yfOLiMjpS+HGT/lyzA0cCTdbs0tqW5FERERag8KNn/J1t1RieABxoQ7cHoON2opBRERakcKNHzLMVpxuD+C7cGMymTijUyQAP+7J90kNIiJyelK48UMeWxAAFpMJeyO7fLeGM7tGAfCDwo2IiLQihRs/VBNuAu0WTCaTz+o4MyUagNV783EdbkkSERFpaQo3fqgm3PiqS6pG74RQwgNtlDndbDxQ7NNaRETk9KFw44c89rYRbsxmE0O7HO6a2p3n01pEROT0oXDjh47ulvK1szTuRkREWplvVniTFnWkW6r1frw7dmxn2Nmj652vDo6HwTfx05583B4Di9l3Y4BEROT0oHDjhzy2YKB1u6VcHoPrH51TvxbD4MXPN1ECbMkqpt/hPadERERairql/FBbGXMDYDaZsBVnAvC9xt2IiEgrULjxQ7Vjbmy+DzcA9qJ0AL7ZmevjSkRE5HSgcOOHfDHm5ljsBXsB+G5XHpXV2mdKRERalsKNn3F7DAxbINA2uqUALOW5JIYHUOXyqGtKRERanMKNnykod4LJ+2NtK91SJuC8XrEALN12yLfFiIiI31O48TO5pVWAN9iY29C063N7xgGwdFuOjysRERF/p3DjZ/JKnUDbWMDvaCO7R2OzmNibV87e3DJflyMiIn5M4cbP1LTctJXxNjVCA2wM6exdrVitNyIi0pIUbvxM7uGWm6BmHm/jsBicn2IhuWg18SWbMBknP+upZtzNVxp3IyIiLahtzBWWZnOk5aaZfrSGh/iyrdw9yM3vhwbDlrsByAvswg/Jk9gecyGG6cSC1AWpccz4dCvf7cqlqKKa8EBb89QoIiJyFLXc+Jm8mgHFzdAtZTJc9D30X7oWfEuQDfYXe8gN6EKVJZjoir1csv2PXL1xCrbKPAzDwO124Xa78LgbbtXpHhdKj7gQqt0GSzYfPOX6REREGqJw42dqBhQ3x5iblIKVhFVl4zLZWLTHRLfnS4md/jOxT2bx8FeVlFQZdCpew9jPLyTJyOGhcX15aFxfnrrpvEYDzqUDEgFYtCHrlOsTERFpiLql/ExzDSiOLd1GfNk2DGB7zIX8+O1inB746qV7sFq9f222V+fTN38JPaNLWTs5hC1J11BqCWP0/zyPgdHg417aP5FZX+xgxY5D6poSEZEWoXDjZ2oHFJ/CmBuHq5iuhd8CkBE+hKKADrW3WcxmrBZvg5/TEsPG+CvofegzQsijf/5nrI++tN7j7dixnWFnjz7yGINvoTo4hrNvvJ/Onmz+8+F7Ta5VRETklxRu/IhhGEcW8TuFlpsOxeswG26KHIlkhg465rXVliC2xI4jZdc/iQkso3/eZ8QE1V080OUxuP7RObXff787jx/25BMz4lpyPny8yXWKiIg0RGNu/EiZ002VywM0vVvK7ioltmwHAOnhQ8F0/FWOXZZA3t5qpdISQqC7mA8nBGLxVDV6fY+4EAD25ZXhsQY0qU4REZHGKNz4kdwSb6AwuZ3YLE370XYoWY8ZD0WOJEod8Sd8v9JqE1tjL6baZGdkJysX7/ozGJ4Gr40OcRATYsdjQGVsapPqFBERaYzCjR/JKzscbpzlTbp/iM0grnQbAPvDBp/0/StskWyJuoBqt0Fq3hKGp89p9Nq+SeEAVMb3xzAaHnwsIiLSFAo3fuRQiXcwsbm6aeFmUKwHM26K7fEUOxKb9BhFjiTu+KQSgLP2v06fg580eF3vhFAsZhOukDg2ZhY36blEREQaonDjR2pabpoUbgyDvjHebqSDIaknNNamMW+ur+a7DrcAMGbXU4zqUH/NmwCbhW6xwQAsXJXe5OcSERH5JYUbP5Jb23Jz8rtuh1TnEh0AbpOF/MDOp1zLtx3vYGvMWCyGi7curiSyfG+9a2q6pj5ad4DK6pPfq0pERKQhPg83L7/8MikpKQQEBJCWlsaKFStO6H7ffvstVquVQYMGtWyB7Uhty00TxtzEVewCoCCwCx6z/dSLMZn5vMcjHAgdQGQAXLnlfgKrC+pckhwZiLmykJJKFx+vP3DqzykiIoKPw83ChQu5//77efjhh1m7di2jRo1i3LhxpKcfu5uiqKiIm2++mQsuuKCVKm0fata4OdluKYsJYit2A3AoqFuz1eM2O/g49Vl2F5mIqMzkis0PYHMdaVUymUwEZq0HYO6KPRpYLCIizcKn4WbmzJlMmjSJ22+/ndTUVGbNmkVycjKzZ88+5v3uvPNObrjhBoYPH95KlbYPNasTn2y4OT/Fgt1TQXk1FAUkN2tNFbZIfvWfACqs4SSWbuLyrQ/WWQMnMPtnguwWth0sYcWO3GZ9bhEROT35LNw4nU5Wr17N2LFj65wfO3YsK1eubPR+b7zxBrt27eLRRx89oeepqqqiuLi4zuGvata5MTtPbszNhH7e/Z225JsxTM3/V2JHoZkP+rxAlSWYTkWruGzrtNqAY3ZXcd0Qb6B6dcXuZn9uERE5/fgs3OTm5uJ2u4mPr7tQXHx8PNnZ2Q3eZ8eOHUybNo233nqrdvPG45kxYwbh4eG1R3Jy87ZMtCWHasLNyQwoNjxc0t37Xm4raPoMqeM5GNqHj1JnUm120LXgGy7f8iBWt3fK+G0jUzCbYMWOXLZll7RYDSIicnrw+YBi0y+mHBuGUe8cgNvt5oYbbuDxxx+nZ8+eJ/z406dPp6ioqPbIyMg45Zrbogqnm5IqF3ByLTex5TtJDDXjNlnJKGm5cAOQGX4GH6X+H05zIF0Kv+fKzfcRbHHTKTqIi/omADB76c4WrUFERPyfz8JNTEwMFoulXitNTk5OvdYcgJKSElatWsXdd9+N1WrFarXyxBNPsH79eqxWK1999VWDz+NwOAgLC6tz+KOawcQOqxmT23nC90sp/A6AQnsibqNlww1ARsRQPuj7N6oswSQXr+G1odugMIMp53UH4KP1B9iZo9YbERFpOp+FG7vdTlpaGkuWLKlzfsmSJYwYMaLe9WFhYWzYsIF169bVHpMnT6ZXr16sW7eOM888s7VKb5NyDndJxYY6OJmIklL4PQAFAR2btR6321V7GIZR5/uM4L4s7PMSpbYYuoVUwmsX0J8dXNgnHsOA//tiR7PWIiIip5cTG7jSQqZOncpNN93EkCFDGD58OK+88grp6elMnjwZ8HYpZWZmMn/+fMxmM/369atz/7i4OAICAuqdPx0dOircZJ7gfeyuUjqUbAAg35EMrDrlOtweDxYTTLtsYJ3zD43rW+/ah8NMfHpjMP04CG+M44kRj7Fkc2f++3MWd48uJjXRP1vZRESkZfk03EyYMIG8vDyeeOIJsrKy6NevH4sWLaJzZ+8KuVlZWcdd80a8Dh3ulooNOfFw06nwR8y42ZrrpioptFnqMAxwG/DVS/fUDvpe+v6bnHf1xHrXuj0eRtz3PIVv3oB5639IXPEH3os5j0m51/Pc59t5beKQZqlJREROLz4fUDxlyhT27t1LVVUVq1ev5pxzzqm9bd68eSxdurTR+z722GOsW7eu5YtsB45uuTlRXQq8U+4/2+lq9nosZjNWi/ewmE21fz76sJjNlDjBc808GPM4mCyklS5lieP3GNsW8e1OrXsjIiInz+fhRppHU8JN56IfgJYJNyfFZIKz74fbl0BML2JNRcy1P0fZwt/iKis47t1FRESOpnDjJ0423IRWZhFWlY3bZGFFum83rXS5XN4jfiCu27+iPO1/8GBibPWXOJ8fgnvNP3G7qn1ao4iItB8KN37i6DE3J6Jj8VoADgb3ptxHuaFm8HFgYCA2m817BIYSfPnTXPRNf/Z44gly5mL5+C423BuHe2/jK1eLiIjU8OmAYmk+uSfZctOhaA0A+0MHAd+1UFXHVjP4+MXru2Ex183ZHg4xsfx+xjl+5n7bBwyKq4R546DfNXDedIjp4ZOaRUSk7VO48QOGYZx0t1RNy83+sMEtVteJ6n/OZdhs9f8q5q58j1djH+A99zlMyZjGrT2KMW18Dza+D6nj4ewHoMMZPqhYRETaMnVL+YHiChdOtweAmBPolgp25hJZmY6BiczQAS1dXpMlVu/nji6HyCWcRyKepOCGz/D0HAcYsOVjeHU0nnnjcW/9DJezCpfLhdvt2/FDIiLie2q58QOHSr0bUIYFWAmwWY57fYfDrTaHgntQZW2e9W1aygPdD/LZwTD2EU2vB98j95OF9Ik18/sRdm4cYMO6dznsXU5GkYd//FzNN7kRfLJ6P5YT2FjV7XZjGMZxrzOZTFgsx39fRUSkbVC48QM5DXRJedxuDOp/cBuGQVLhagAyQgfhdvt4GvhxBFgM/pK6h1+v6klwn/N4+douXJV4eHp4ZRGeA6sx5WwiObySP4xyABUYfxsM3UZDyjnebquILvCLMT1ut5uUzp3IyDxw3BqSOySxZ1+6Ao6ISDuhcOMHfjnexjAMnrrpPPJzcxq83vLzvyDewlMvzOP9La957+PxtE6xJyEvL49X35iH0+WhaG8aEefcxMMbEtj91dskO8q4YcIE6HGhN8jk7sCTs5nK7O0EFaXDmje9B4A9BOL7eo+YnhCWhBGcgKsoi4pPH61dSbkhLreHwIsfPaEWHhERaRsUbvzAkXATcPiMQX5uDl/Pvq/eLKTvP5rHgHjvuJSH7rmdKR4bY+7+W2uWe8I8BvQ7dzyVThdF7/yNMy+8jG1VkSxKvocJ+548cqHZCnGpeKJ7Evu7P1K8/r9Y9n0D+76FnC3gLIWMH7zHYVbgwO9CMb6fhckRCo5QCAg/6oiAwEiwBrf66xYRkVOjcOMHGlvjpmYLhKN1PLwXZaU1DMMejNnZtrulahke/id2A49kncU+ZyiLwq/jfwzv4sZHK68Go8dYSL3Ee8LtgrwdkL0RDm6Egr1QfACjaD+eokwsZqCq2HsU19+Vy2INZPFvgjB/+Zi3i6vLKAiJbelXKyIip0Dhxg+czDTwDsHe7pUSe/v7gI6yVnFfwib+nDmQjUFpvJGeyW2d8+pd53L9IrBF9fAefa6qc01IUCCVn/weq6sMqkqgsuiooxAqCjC5KhjbzQrfHWndMhIGYnQdjafb+dBxGFhsgAYei4i0FQo3fuBkwk1SiDfclNrjWrSmltInqJAbY3bxj9wePLk1iU6BTsbElQDgdh9Z8fhEeezBEBTRyI0unIXZ3PPUqwxKtDC8o4VBCRZM2esxZa/HvHIWeeUe/r3FxdsbqtnrimP3vgwFHBERH1O48QMnHG4Mg6Tg9h1uAMZF7GfDvjzWBZ/FPT93ZuHQXQwIr8DAwG1A2aJHsTewKODRKp3VhF76BA1MKDvCbMUTEscra6opW/QH7DYrLmcppsJ9mAr2YCrcR3RQOXem2bkzzc7+4hJMS/4Ig2/0Dl4WERGfULjxA7knuK+Uw11CkA08mCmzR7dGaS3CZIKLi/5NWKe+LM8L5bY1XXh32C4SbU4ArBYz1uO0nlgtJ7fYX+1jBoZD4ABIHACGBwrTIWczxqFtdAyrgu9f8h6Jg2Dwb7zbRQRFNfWliohIE2iF4nbO5faQV+b9UI8JtR/z2pAq79TwMns0hql9d51Y8PDSoH30Ca0g12njxlVdyaw89utvdiYzRHaBXpfgPnMKV75Tjqf3ZWC2QdY6WPQgPNcL/jURdizxDm4WEZEWp5abdu5QaRWGAVaziejgY7fchDoPAe27S+pooVYP89N2M+GnbuwqC+CWtT2xhPqoRcps5aNtLpxXzsXqLMK08d+Y1y/AdHADbP4QNn+IEZqA0f86PAOuxxTbS2NzRERaiMJNO5dd5N16IS7UgcVsOua1IU5vy42/hBuAGIebt4bs5rofu5FeEUDCDU+TXn6I7uGtu+heY4OZB8abuXWwjRv724ghG9PKFzCvfIHVB80Mvu1ZzP1/5V1XR0REmo3CTTt3sNgbbuLDA455nclwE+z0TpsubYfTwI8lIcDFgqG7ueGnFPaRwA1rovjHkD30Cq1qtRqOO5jZ48advwvTwY2YCnaTFu+B/06FxQ97dzgfdAOknFtvmwgRETl5CjftXNbhlpuEsGOHmyBnPmbclFd7F/Br72q2ZjjaJZ5gnndMIie2C1d824nXhmZxdnRpq9bV6GBmiwXiUyE+FVdFEdOe/D+euWEgptxtsOFf3iM8Gfpc4d1SotNwsB5/ar+IiNSncNPOZR9uuUk4TstN6OEuqQNlpvrL+rZDNVszHK3S6eLg1GmM+f1stlVFcsvqFJ5IzeSG5HwfVdkIewjPfefkL8u+xXrwZ1j3Fmz8NxRlwHcveg9bMHQ9F0/X0Xg6DoPYVDA3PkZHCwiKiByhcNPOHTzBlpua8TZZZSa//qF7Kkv5fcJqZm6JZFNQGn/Y3JH1RUE8nppJgKWNbX5pMkHHNO9x0VOw/VPY/jns/ALKcmDbIszbFmEGiqsMftjvZuV+Nz9lulmd5Sa79Mjr0c7lIiJH+PPn3GnhRFtuasJNZqmJzi1elW/ZTAZn73qB2M7XsSx0HAszo1i2p4wrC/5JtPsQwUGBXH3VVcd/oBZWZ5sIkxV6jfcehgcObsTY/jlfv/EnLuwZRJjDyYXdrFzY7cg/WcMejBGcgCc4jquf+RKjKBMik/2iZU5E5FQo3LRzNbOl4o/RcmN1VxLoKga8LTf+Hm4ADAPuTAtieNnP/C27D9n2ZN5I+D2/id1J/Lo5Pq3tZLeJqFr8EPaqIu/GnsX7oeQglOdhcpZhcu7CXLCLj68Pghf6YwTHYSQOhISBGIneg9CkOoFHXVgi4u8UbtoxwzCOtNwcI9yEHF7fpsIaRoWrolVqaysGBBfwdOefmJ2dysaKKObm9KJz9BTOLqv0WU0nuk1E7RYRmCEkznskDfbe6HZCaQ6UZOMuPsDmjRvoE2vGUpaDaecS2Lmk9nEOlnpYk+VhdZabHzLd7HLGsmGn9sASEf+lcNOOFVe4qKz2AMfulqq7vs2+1iitTYmyOpneYT2fFXZkYV5X9jm6c8VPHsLPvpEqtwm7zTd1HW+biGNuEWGxQ3hHCO9ItdPJgLu/p+yTP2B35mMqPXj4yIbyPOJDzIzrYWZcD+8/92p3KeZ5l0L386Hb+ZB0hqagi4hfUbhpx2pabSKCbATYGv+QDKldmTiW0zHcAJhNcEnkfoaE5PLC1nB2BfQhYuT1XPZjJY+nHuDcmNJ2P1TFandgDe4EkZ2OnHRXewcnlxyE0iyMwgxslYWQ8Z33+PrP3pWTe4/H0/tySD6zdlaWuq9EpL1SuGnHTqRLCsPwy5WJmyrOVsmE/HeJu/Au7v4+gn1Ec8uaroyMKmF6r2z6hflZt53FBmEdvAfgdFbT7/rHOL+rlbFdrYzpaiWcbEw/vYr5p1fJKvHwr83VzF1TTaEtQTOwRKRdUrhpx7KLvB/ExxpMHOAqxuapavc7gTcnE3BxXAEHXpvOozP/xlv7Y/k2P5TLvgvlisQCHuyeTXJQta/LbBEGBjsLDNb/6Q/e8T4eF+7CfZhyt2HK20liaBX3nengvjMd/JhZBKvnwYBrIaD9L/woIqcPhZt2LLvIu71A4jHH23i7pMrsMe1+J/DmkpeXx9w3/4nhLCdu5QzusMewNHQcm4LS+CgrkkXZ4UzomM+UlENEWZy+LrdF1I73sVggtqf38LihYA9kb8DI28GwDsCiqRhL/hej37V4hkyC+L71HkvdVyLS1mgUYTtW0y11rJabI11S/rWf1KnwGNBn1KUA9B11KWefPZL/HVjMU8k/kVK1jWrDzD8zYjh3RS8e39YJS2iMjytuJWYLRHeHvlfhTLuThz6vZMshN6bqcsxr38T66jl8d0c4Nw4KJMhhw2bzHimdO+F2H2Pws4hIK1PLTTt28AQW8DsymFjjbY4nJaCUG/IW0O+Ku3h+ZzzfF4Sw4EAcHe58lce2FdBx8zxsZdkN3tdm9699oAxbIM9+5+Tx/52GqyIbc9ZaTLnbGdXZyqjOVgxbMEbCAJxxAwi8/BkMo42t/iwipzWFm3bseJtm+vNO4C0lLy+PjR+9xAVAd3s3loWMJSOgB29nxmEOfYALOh3kish0om11dxxf8+VHvim4hVmtFqxRXSCqC1SVQNY6yFrnXUAw4zsc+39g4a8CYd+30PUcrY4sIm2Cwk07dvA43VLB1QWYceMyO/xiJ/DWcPSGnP2AC5x7GTfjdUbf9hBbKqNYUtSRr4uTOC8siysi9xHzi5Dj1xyh0GUUdBoBudvhwGpMRfu5rq8N/nE5xPWBobfDgAngCPF1tSJyGtOYm3aqyuUmv8w72LWxAcWh1UeNt9Fv1E1WlbGB6QmruWzXDPoEFuAyzHxR1IH7957FKwd7cbD62Pt6+R2zBeJSYdBvcA2eyJzVTgxbEORshv9OhZmp8On/g9wdvq5URE5TCjftVE6xt8XAbjUTEdTwEruhzlxAXVLNJalsK3/suI4/dlhLn8AC3Jj5ujiJqXvPZFHUDVgjk3xdYusLjmPyfypx37cRLv4LRHWDqmL44e/w4hCYfwVs+Q+4Xcd/LBGRZqJuqXbq6PE2pkZaZUKrvYOJSzSYuFn1CSqkT9A6tlWE835+Z34uj2ZTyDCSbk9j6qYi7u1+iJ4hp1F3FUBAOJz1PzDsTtj9Nfz0Gmz/DHYv9R7hyTDkVjhjIgSfJrPPRMRnFG7aqQOF3gX8GuuSCnNAoKsQUMtNS+kVWMT0Dj+zszKUedsD2BXYj/8cjOK/ByMZF1/E3V1z6BPmuw06W5PLdVTLTJdzvUdhOuY1b2Ba+09MRRnw5RMYS/+C0edKPGm3QYchjXaXau0cETkVCjftVEZ+OQDJUUEN3p6WaMEEVFpCcFkCW7Gy00/3gBKuPvRPHvjSwy2/f5zPD0Wy6GAEiw5GMCa2iHu75TAg3M+2dTjM7fZgMUFgYON/xxwWuK6vjbuH2RnWwYlpw78wb/gX67LdzF7l5K2fqyn7xYLQyR2StPWDiDSZwk07tb/A+2GZHNlwuDmzo/dDQevbtJ7qnN282H83uytDeXF3HP/NDueLQ97jvJhi7umWQ1pEua/LbFYGBm4DyhY96t3OoRGVzmpCL32Cin/dif3QBky5WxmUAHMuC+TvV4RjxPfH02EIOMJwuT0EXvyo1s4RkSZTuGmnMgq8H5IdIxv+jXlY0uFw41CXVGvrHVrJiwPTub+bg5f3xPFRVgRLc8NYmhtG56qdnFH8A2Bi7pv/wGoxExwUyA0TJvi67FNSu51Do7d7VzA2hydhjukM1RfAwQ1wYC2migJMB1ZjzlrrnU6eNLS1yhYRP6Vw007VhJvGuqWGdVDLja91D6liZv8M7ut2kHs/K2Rj8DD2ObqzL7Y7SXdcSlb3IkZHHmLnig99XWrrswVCx2HQYah3P6uM76EwHQ5uxHpwIx9OCISDm6DDQF9XKiLtkM+ngr/88sukpKQQEBBAWloaK1asaPTa999/nwsvvJDY2FjCwsIYPnw4ixcvbsVq2waX20NWoXeganJU/ZabWIeTDmFmDEyU2bQTuK91DnJyWdG/eL7L94yP3EeQuRpbVBLz81O5a89wPgu/mo3Fp9laOTVMJojqCgNvgMETIaYXBnBFbxuWV8+FDyZ7Q4+IyEnwabhZuHAh999/Pw8//DBr165l1KhRjBs3jvT0hv8zW758ORdeeCGLFi1i9erVjB49mvHjx7N27dpWrty3sosrcXkMbBYTcaH1PxT7hXtbdcqskXjMDa+BI60v2lbFDTG7ea7DCvKX/J14aznlHhurg8/msu96cunKHsxPj6bAeZoOog1LhL5X4T5jEgs3VmPCgPULMP6WhufTabhKDuFyueoc2rBTRBri026pmTNnMmnSJG6//XYAZs2axeLFi5k9ezYzZsyod/2sWbPqfP/UU0/x0Ucf8cknnzB48ODWKLlNqBlM3CEiEIu5/lTavmHecFOiKeCtJi/Pu4dXzTiaX8rPz6/9c4DZTcma//D07d3YVh3LR7sMdgYNZFNJII9s6cDjW5M4M7KU86PzsYSefi1vbkcEN75fwTMrq3h6TAAXdAXTD7PJ//olpn1Rxetrq6kZaqxZVSLSEJ+FG6fTyerVq5k2bVqd82PHjmXlypUn9Bgej4eSkhKioqIavaaqqoqqqiMLqhUXFzet4DbkeNPA+xxuuSmxKdy0Fs/hT9u+oy7F1sCsoa/fm1fvnNkEA4ILMBd8wq8usfNhVgTvZkaxuSSQlfmhrMwPpeOUN7n6pzLOiSljRHQpaRFlBFr8exZRzQys5XP/iN1qwV24F/OepcSQy2uXB/LqjV1xdxuDKzBOs6pEpEE+Cze5ubm43W7i4+PrnI+Pjyc7O/uEHuO5556jrKyM6667rtFrZsyYweOPP35KtbY1GYdbbhqcKeVxk6qWm3Yn0u7m1s553No5j/RyO4tzwliUFcaaoiA2lgSzsSSYl/fEYTd5GBxRzplRZZwZWcbg8DLfD5xrIVaLGavVCjHdISoFDqyBvSswlWRhXfcPzAkDiThNhyqJyLH5fLbUL7cOMAyj0e0EjrZgwQIee+wxPvroI+LiGp8RNH36dKZOnVr7fXFxMcnJyU0vuA3YXzsNvIGWm9wdhFg9lDoNyq0Rvv8By3Hl5eXx6hvz6p0f4/Lw8ce53Hn9JWQE9mCvowcllgh+KAjhhwLvrttmw028cz8Ro29j2sI1dHHtI9Ao94vp5XWYLdBxKMSmerd3yNmEOXs9W+4KwbTxPRh4nTaHFZFaPvvsi4mJwWKx1GulycnJqdea80sLFy5k0qRJvPvuu4wZM+aY1zocDhwOxynX25bszz9Gy03mKgBWHXBj6eKvv9P7F48B/c4dX+98pdOFZ8Hz/DotFputFMNYS3Z1IJsqItlaEc66PBtl9miyHJ0JH9aZ9w/fr6O9lLiC9QRnRTAssozEgOp6j91uOUIgdTwkDsDYvpgE8uHD38KGd+DS57wzr0TktOezTz+73U5aWhpLliypc37JkiWMGDGi0fstWLCAW265hbfffptLL720pctsk/Yfa42bzNUA/JipWST+xmSCRHsFY8IPcHfCFm7cOpUXuqzkjugNlKz7lERbKQD7nSGsCR7JfT93YviyVM5b0YvHtiSxPDcEp8dPWjciOuMePJE/fl2JYXHArq/g5eGw/FlwOX1dnYj4mE9/tZ86dSqvvfYar7/+Olu2bOGBBx4gPT2dyZMnA94upZtvvrn2+gULFnDzzTfz3HPPcdZZZ5GdnU12djZFRUW+egmtzunykFV8eI2bhrqlFG5OK7G2KkaEZJO/+CWe7vAdc1K+YWriBoaVLqV/WDlmDPaWO5iXHsPNq7ty1jeDiL3qYRZmxpBd2c47Lc1WnlzuxP3b5ZByDrgq4as/wZxRsO87X1cnIj7k0//dJkyYQF5eHk888QRZWVn069ePRYsW0blzZwCysrLqrHkzZ84cXC4Xd911F3fddVft+YkTJzJv3rzWLt8nDhRWYBgQYDMTE2Kve2NVKWRvBOAHhZvTUpi1mqEhuZTvmUd08cdcZHKw19GDnY4+7AxIpZRwgnoO54/b4I/boF9YOWPjihkbV0yvkMr2OWwlujvc/DH8/C9YPB0ObYU3LoYzboYxj0NQ47MpRcQ/+fxXtylTpjBlypQGb/tlYFm6dGnLF9TG7a+dKRVUf+B15mow3GRV2NhfrOmxp7Ojx/GkAVCBYaxhe1kg9y/cwXmXXcX64mA2FgexsTiImTsT6BRY5Q068UX0DWpnXTsmEwycAD0uhC8ehTXzvcfWRXDRUzBAA45FTic+Dzdyco65YWb69wCsLwwG8lqxKmkPTCbobC+h6LuFXJSymvNsYexw9GF7YD92O3qRXuHgtX2xvLYvliB3CdHj7mX6wjV0q96JDe+g5LY4C8vlch35xh4Gl/wf9LsOy6LfYcrdBh/8FmPdW5gu+z+I7ua7QkWk1SjctDPpNQv4NTTeJqMm3IS0ZknSDtUsNngmAIeo9OTzc3kkq0pj+aEglHJrKCEDxvIe4DC5GRiUx5CQXCrWveXbwo/idnuwmCAwsIGgD9jM8OAIO388x0HgnmUYLw/HM/IBjBH3gLXxBXJMJpNWPBZp5xRu2pldOd4ZMV1jg+ve4HZBxo9ATcuNyIkLMLsZFpLLsJBcen37DyLH3MUfPz1Al2EXku8O4MeyOH4si8OU8AQ//OQdp3NhXDFRFt91X9WsZFy26FHsDawKXcNZnMvnH7zK2G5gWf4Xdn74FPd9VsmiHa4Gr9eWDiLtn8JNO7PzkDfcdI/7RetMziZwloIjjF2lWrZVms6Mmz6B+RR8MYcPrg5gvyeCVWWxrCqNIcMZwrf5oXybH8qjWzvQLaiCyDGTWZwTwdmxlUTaW38gu9VixnqMIOIKjuSif5ZTOX8C9ozldI8q4783BOGJ6oYnZTQERh651u3Rlg4ifkDhph1xujzsy/N2S9ULN+k/eL92HIpH422kmZhM0DWglK4BpVwXvYfl33xH1Dm38PnBcFYVBrGrPJCwtMu4ZyOYMOgTWsnwqFJGRJcyNLKMUKvH1y+hlim+D6aE3rDvW8hchTl/F+aCvZB8JnQaDhabr0sUkWaicNOO7Msrw+0xCHFYSQj7RetM+uF1PToNB/7T6rXJ6SHKncsdXbxHgdPCikMObnlzK2eccwE7ywLZVOI9XtsXi8VkMCCsnJHRpZSu+QjM1gZ3TW/VQcpWB3Q7HxIGwM4lULgP0lfCwY3Q/QKI0IBjEX+gcNOO7Dg83qZbXEjdaeCGUTtTik5nonAjLeWX+2A5XR4KvtjJNXFfkFMOBXFp7LV3Z5+jBwXWGNYWBbO2KBgifkvyfTfxWXgpg0IKGBKcS4ytCoCNyz5p/RcSHAMDfg2527yrG1cVw6YPMEd0oWe0ti0Rae8UbtqRnYfDTffYX3RJFaZDyQEwW6FDmg8qk9PFL/fBqnS6YMHz9B11KQUf/5MJwzoClcAGDlU72FQeyYbyKNYUBFJpD2N9RSDrK2J581BPUhwlDAk+RIQ1EcPwwTI0JhPE9vbuR5X+PWT8gLlwLxv+Jxjzl4/Bef8PHKGtXJSINAf9itKO1IabX4632bPM+7VDGtg1U0rahlhbFeeFZ3NP4mYmbLibA2/cw4TI7fQOKMSEwZ6qUN7N78qrcQ9x0cqe/H1PrG+2hLDYvds3DL0dT2RX7BYT5u/+Bi8OhQ3/9raMiki7opabdqTRcLP7cLhJObeVKxI5MSYMqnP2cGn4Pq6MyaTIZWNNWQyrymJYXxLG9tIA/rI9kb9si6dz5XaCUs9lzptvE2AxWm9MTmAknr7XcPmUx/lkSh9MhXvhvUmw6g245BmI79PyNYhIs1DLTTvh9hjsOjwNvMfR4cYwYM9y759TzvFBZSInL9xazejwLB5K2sBvNt/H7XFb6RVQCCYz+wJ7E3v5Q7za+Sk29b2XrEpHq9b23x0uqm5fhvvc6RjWQNj3Dcbfz8az6Pe4SvNwuVy4XC7cbu3fJtJWqeWmncgsqKDK5cFuNZMcddTqxIe2QlkOWAMheZjvChRpIoennNHhWVwQnsXB6gC+Lojjvf3hFIfG8GFBF0zxD7NjbSk3d8pleFRZi47NqV31ONS79k2ncBPPjQ3gV33A9OMcDn01m//3RRXz11fTUYv9ibRZCjftxM5DJQB0jQnGYj7qf/eaLqlOZ3mnuYq0Y/G2Sq6M2M0LD/+NZ5/5E1+WdmJzRSSf5YTzWU443YMrmdgpj6uTCmiJSNHYqsfugr2Yd39JPPnMuzKQ13/TlSHP7dBifyJtlLql2onjDibuqvE24kcMD0ODc/hjx3X8NudpbkrOJdjiZmdZAH/c0oGzlqXyl50dsYbHt8jT16x6XHNYYrphGjIJuo4Gix1zaRarfhuMedHvoDy/RWoQkaZTuGknGgw3bhfs/cb7Zw0mFj8V6zrIn/oc4PvztvBY70xSgqoocVmYl5FA0m9fYcrPXfkuP7jlJzWZLd7VjIfegSc2FbPJhHnNPPjbGbDqdfBoDI5IW6FuqXZi28EGwk3WOu/iYwHhkDjQN4WJtLBfLhx4PSZ2OXrzQ9Ao9gb25ovcSL7IjaR3SAW3ds7lisRCAiwtmHQcoXh6XcboZ1ax9PdDMeVshv88AKvfhEueheShLffcInJCFG7agSqXmy0HigHo3yH8yA07Pvd+TTnH+1uliB/65cKBAP2Bi50ZXPTH57nk179hU9AQtpYG8v82JfPYz5EMKv+R4cY67r12TIvVtXyfm6pblmBbPx/z0hmYstbB3DF4+lyN55yHIKYnACaTSYOORVqZwk07sCWrBKfbQ1SwnU5Hz5Tatsj7tdclvilMxMeq8zKY2q+cKvOPfF2cyOeFHcklhO9Cz+c7zufbH0u5JiEHk9XerM9bO6sq2LuCcWyQiRljHEwabMe8+X2Mje/x1oZqnlhWhTMoUbOqRFqZwk07sC69AICBHcOP7ClVtB+yN4DJDD3G+rA6Ed8LsbgYH5nBJRH7WVMWzVdFSawri+SHghB+KAih413z+cOWSq5IKmZ4VCnWUxxt2NisKlfpQczpK7Hk7+TmgXZuGujgzXX5GAc3QdKAU3yVInKiFG7agXUZhQAMSo48cnLbp96vyWd6NwEUESwmg6EhuQwNyWXlihVYR/6WhfsjOUAI/84K4d9ZMcTYq7k0oYjLEwoZHFGO+RTWzamZVVUrPAn6/wpKsmDvN5jyd3HLIDu8Mso76P+sKd5fRsyayyHSkhRu2oHacNMp4sjJmnDTa1yr1yPSHlQf2kXYir9yvcvgoW8CuPjyy9gWOJBcZwhvpsfwZnoMHQKcXJZQyEXxxQwKL2++Jw9NhP7X4irM4P2353FtPwemPcu8SzdEdYO0id5dyUNbZiq7yOlO4aaNKyhzsjfP+5/uoI4R3pOVxUe2XNB4G5EG1QxErnS6qFrwPFP7nY/JuooN5ZF8utPJ7tAzyKy0M2dvHHP2xhHvqOaCmAICOg3A5YFmGaUTmsSEf1dw9exNWNfMhdXzIX8XLHkEvnjc24oz+EbocRE087ggkdOZwk0bt25/IeBdmTg8yOY9uesr8FR7fwOM6eG74kTaGavJYHBwPsW75jEuOoGdAalsDRjAzoA+HKwK4O3MOOKvf4oR37i4ML6Yi+KKODu69NSnlkckw9gn4dz/591pfN1bsP8n2P6p9wiIgNTx0O9q6HIOWPRfs8ip0L+gNm5deiEAg5Ijjpzc/KH3q7qkRJrEY8Dgcy9mMABFVHu+Z2NFJN8Vx7A0J4TCoHDezYzi3cwo7J5KulVtoXfFBgaY93DrhKtO+vlcLpf3D5ZAGHST98jdhnn9Akwb/oWp9CCs/Yf3CIqBPpdD36uh8wgt8yDSBAo3bVy98TYVhbD18BTwAdf5oiQRv2Mze1t0Um05/PORv/H35x5hTWUCP5XFkO8KYEvgYLYEDuZjo5pv1lRwfnQe5sCw4z5u7ZTxwMBGrzGb4JzOFib0tXFtXzvR5HpXPF71OkZIPEbv8Xj6XOXdGNfkHYistXNEjk3hpg3zeIyjZkpFeE9u/hDcVRDXBxI0tVSk2RkeUgMLGBBWwkRjB7urQvmxNJYfS2PJrg7iy0M2vjwURse7/8HNa8oYl1DC2PhiEgOq6z9UI1PGG1LldBE1/nHO7WJhQj8bV/e2EclBTKtew7zqNfYXe3hvSzXvbnKR7o5lx+69xw04CkFyulK4acM2ZBZRVFFNsN1C74TDvyWuf8f7deCvwXQKc1hF5LhMJugWUEK3gBJ+Hb2bDxctITthFFsD+pNj78j3hWF8XxjGo1s7kOhMp7dnF1POSWZweHmdtXTqTRlvgMvixumBD1/8X28Q8rhxF+7FdGgrpvyddAxzct+ZDu4708GBklJeuSqCdzdV822GG08jQ4KSOyRpAUE5LSnctGFfbjkIwDk9Y7FbzZC/G9K/8zZN91eXlEhrMpkgojKTq9ICqHRuYez0R3lw+u9ZUxHP9spwsuydyKITX/8IoVY3I6JKSQsvwhbXFfdJjEeuDUIWC8T29B4eF+TvgdytGLk7SAp1cs8wO/cMs2PYgjFiemBE98II71jbdeVyewi8+FGMFt9RVKTtUbhpw77YkgPABamH18JYv9D7tetoCEv0UVUiAuAqOsi48HQujzlAgcvOz+VRLN9TwYHwARRWW1mcE87inHCSbn2BocvdpEWWkRZRTv+wCvqFVRDrcJ34k5mt3pmRMT2oqiznmvv+wsf3DcVSsBNTdZl3X6usdWALgpheENsLQju01EsXafMUbtqoA4UVbM4qxmSC0b1iweWENW96bxx4vW+LE5E6Iq1Ozg3LJrrwE2678hY2FAXyTV4o3+cHsjzLRqkjiGW5YSzLPTIIOd5RTf+wCvoeDjs9AopO7MnMVhbtcOHucTEWixkK9kLuVsjdAdXlkLUWstZisQUx+9IATLuXQrfzNL1cTiv6295GfbnV22pzRqdIokMcsPYt75LuoYneaaIi0ubk5eXx+rx5ADiA4S4Pb72zm0dvHsUuOlAQ1otsWwdyrXEcrLJx8JCNLw4dCTwd73mLW9aa6B9eSZ+wSvqGVpASXIWlseF1ZgtEd/MePdxQuA8ObYXc7Ziqy5k8xA5vXwNB0dD7MuhzBXQ5G6yOFn8vRHxJ4aaN+mKzd7zNmNR48Hhg5QveG876H/3HJNJG1ayKXKPS6YIFz3P+WYOwf/xPRo9KBIqp9GxnX1UIe6pC2FMZysZDLgocSViCwllZACsLjgQem+GkX0R1bdjpEVjc8C7nZgtEdfUePS7Cnb+XufPf4o5RiZjK87wtv2veBGugd/2crudBt9EQ11d7XYnfUbhpg8qqXHy3Kw+AMalxsHOJ97cxeyik3eLb4kTklAWY3fQKLKJXoLcr6utv55F26c1c/uQn/GHqb8lwhbOvKoT0qhCqsLO2yM7aouDa+yc/8C7jf6iib3gl/Q53a6WGVhJq9XgvMFswIrtw538que29zVj3fw+bPoRti6D0IOz60nsswbs6csch0HGo92uHNAiMrFezSHuicNMGfbYxG6fbQ6eoILrHBsN//897w5BbISDct8WJSIuwmz04s3cyOjQTm83bcusx4ONFi6mK7ctBaweybd6jwhLCtrIgtpUF8f6BI4+RElRF37AK+oZWkBpS5l1o0Gz1ttJ0PQ8u+z/vL0q7vobdX8Peb6GyEHZ+4T1qxPTyBp2kwd4jvi/YGl+IUKStUbhpYwzDYO43ewCYMDQZ09ZPvNO/LQ5vl5SInDbMJgivzGL0sIsOnzlARVU6Y6f9k2cfe4BvNuyhKKwH2fYOlFgi2VPuYE+5g/9kRwCQfO/bnP3XpXSPC6F7bDBdY0PoFhtN99SJRA25A5PHBTmbMGWugv2rMGeuxlSwG3K3eY91bwFgmCwQl4qRMBCSBmHucAbE9wNbgG/eGJHjULhpY77bncfmrGICbRZuHBwNr0/33nD2/RCW5NPaRMT3TCZwl+RxRtAhyg9+yOizbwEOUeyysbcqhD1Vod6vlaEcdAWRXVxFdnEV3+zMq/M47opiqvP2U52XgSt/P0ZxNpV5hUQ4nQyLc3FmRwtpiWaGJFmICwYObsR0cCOsPxx4zFaI7Y2RMBAjcSBG4iBvC4814HCdWh1ZfEfhpo2Zu8LbavOrtI5ErHoeijMhohOc/YCPKxORtizMWs0AawEDggsA72DmC+5/lf93w0gKbfHk2uLJs8aTZ4unyBKJJTAMS8c+BHTsU++x9nhKsUaaOBRYxeaASnrZD9HD2EtCVTrfr1zBGYlm4oJd9QJPtdtg0yEPq7Pc7CoL5U/zFmNJ7K9p6NLq9DeuDdl1qJQvt+ZgMsGdPYrgvRe9N1z8tPq7ReSkGc5yxp3ZB1vtvlYHgYNUecx8/OkXJJ3zaw44g9hfGciKncVEdexCmcdOmTmEdUWwrnYQcxIwEABP/G/pGuxmQEAOZ1h209u0l87ufUQ7M7FTzqAEC4MSLEAVvHaed3ZW0mDomAYdhnjH8oR10PYx0qIUbtoIwzCYsWgLAFf1sNFx8R3gqfauTdFrnI+rExF/4jB7iKlMZ2Sodz2tSqeLdx99noWz78NpDmDlDz8yZMyvSK+ws6/cTsbhrwcq7WBzsLcS9lZ24WO6HPWoBknk0d+8m6HWXfSo2sTQoGyCXOWQvtJ71FwZEo/RIQ2jwxBISoMOg7EEarKENB+FmzbiX6sy+GJLDsEWD09WP+PtjoruAVe+rN9wRKTVBFtc2LN/JvPTTCxA18MHQIULfreoiBcf/R8KCSLPFUBetYNcVwB5LgcHK4JZzDAWO4eBCUwVHrqashhs3skg004GmXfR25SOtfQgpm2LvFPTAbdhIje4G8XRA/DE9cWW2JewzgOJiEnEpP//pAkUbtqA9LxynvhkMyGU82ncawRl/QiOMLh+gaZ+i0ir++VihDUqnS5cC54nNbAAm62k3u1fvTePM6+8nQMVVia/sJjf/fY6CowgDrh6ssE1gL+WGnhsQaSyl0HmnYePXXQw5RFTvpOY8p2Q8X7t4x0ywtlj7kx2QFeKg7vgDuuEObIzAbFdiIkMIzbEQWyog6hgOzaLFiKUIxRufCyrqILb5/9EmDOHt0NmkVywC2zBcN1870Z5IiLthAlvy0+yvZKK3au4MHzkUeN9YNn784iMjqHUHEaxJYKvLOF8aDqXFbvKGJls4Yyocnrac+lpO0hn8yFiTUXEGj9Dxc9QAeQeea5sI5IMI5aVRiyHjAjKbVFUB0RhBMViDo3FHhZPYGQCUWGhxIY6iDkqCFnMag3ydz4PNy+//DLPPPMMWVlZ9O3bl1mzZjFq1KhGr1+2bBlTp05l06ZNJCUl8fvf/57Jkye3YsXNZ2t2MZPmfsfF5R8xNeA9gl2VEBwLN/wLOpzh6/JERJqVx4AB515W51xVtYs3Xp3HNgNeP+p8sMNK/07hDOgQSP94K90iDFLCPHQJqiSAKhJMBSSYChjKdu8dDLwBqAI4atZ7sRFEnhFKIaFsMEIoJIRKazguRziegCjMQVHYQqMJCI0hODKOsKg4IiOiiA0LICLQhllBqF3yabhZuHAh999/Py+//DIjR45kzpw5jBs3js2bN9OpU6d61+/Zs4dLLrmEO+64g3/+8598++23TJkyhdjYWK655hofvIKmySmu5B+ff49p/dssNH9JR9vhX0eSz4Sr5kBUim8LFBFpJYYBbgMeGxvJ+dfc0uA1S99/E/s5N9Lv7r/hNiAmyESXCBMpEWY6R1iID3eQEGYnLtRCXBDEBnqItTuxmw3CTOWEmcrxzhSreVKg8vBRWP/5nIaFIkLYaYRQbgmj3BpGlTUcpz0ClyMCT0AEBEZjDorEGhKNLTSGgNAYHEHBBNktBNktBNotBNmtBNosainyAZ+Gm5kzZzJp0iRuv/12AGbNmsXixYuZPXs2M2bMqHf93//+dzp16sSsWbMASE1NZdWqVTz77LNtM9x4PFSWF1FwMIOczL0U799E5f6NdClby+9MmXB4fStPYDTmCx+DQb/RBnYicloym8DayLgZi9mE2WTGbcBXL92D1Xrko2vp+29y5pUTa7/Pq3Yx+O6/8dWLd7Nq0VuMvvRKbO5KrJ4qLJ4q3G4Xhqua/AN7yC/3EBUI0QEeIm0uIq1VBJpd2E1uYiki1lQERiZU4z0qjv0aKgw7pQRSYdgpwEEWdiqxU4WDarODanMAbosDtzkAtzUQwxrgXebDGojZ5sBsC8BitWO1B3gPRwB2eyAWuwOLLQCr3YHFHoDN7sBmD8BmD8LmCMBut2O3mrFbzNgsJg3Cxofhxul0snr1aqZNm1bn/NixY1m5cmWD9/nuu+8YO3ZsnXMXXXQRc+fOpbq6GpvNVu8+VVVVVFVV1X5fVOTdqK64uPhUX0JdJQfhnRuguhyPs4KqimICcQIQDPyyLaYQE6XRAwg76xZIvcz7F7y0tFlKcbtcABSXVWD5RVgqc3ooKj3yL9RZ7T58bWW9245WUW3UXme11l91tOa+Rz9ezXWNPW5D157KY57oaylzeigpq2zwuVv6/WnodbfW+3Mij1lz39Z8f471mG3t/TnWYx7vtbT238kTfS0NXdta78/Rj1lebZzQ+1NW4cRqddfeVukyKKs48n98zeOVVVZTWAW51YHA4XXCzIcPG3y9YQ1/XVrIh0/fQYnNwt7D9zcbbmyeKtYt+4R+wy/E5Xbhqnbyw7qtRAaYiAyAqAAPETYPkbZqImxVRFiqsJoMoAoHVTiAiAZfSctwGybKsFKAlWpsVGPBwIyBCQ8mDJMJMB3+3nvewOw9Z/J+j8lUe59fztA18I6nAngr+Ca221Lh8LmaS2turwlWEUE2XryheYdX1HxuG4Zx/IsNH8nMzDQA49tvv61z/s9//rPRs2fPBu/To0cP489//nOdc99++60BGAcOHGjwPo8++qiB92ejQ4cOHTp06GjnR0ZGxnEzhs8HFP+y+cwwjGM2qTV0fUPna0yfPp2pU6fWfu/xeMjPzyc6Otrvmu6Ki4tJTk4mIyODsLAwX5fTpui9OTa9P43Te9M4vTeN03tzbE15fwzDoKSkhKSk4++z6LNwExMTg8ViITs7u875nJwc4uPjG7xPQkJCg9dbrVaio6MbvI/D4cDhcNQ5FxER0fTC24GwsDD9Y2qE3ptj0/vTOL03jdN70zi9N8d2su9PeHj4CV3ns9GrdrudtLQ0lixZUuf8kiVLGDFiRIP3GT58eL3rP//8c4YMGdLgeBsRERE5/fh0as7UqVN57bXXeP3119myZQsPPPAA6enptevWTJ8+nZtvvrn2+smTJ7Nv3z6mTp3Kli1beP3115k7dy4PPvigr16CiIiItDE+HXMzYcIE8vLyeOKJJ8jKyqJfv34sWrSIzp07A5CVlUV6enrt9SkpKSxatIgHHniAl156iaSkJF544YW2OQ3cBxwOB48++mi9bjjRe3M8en8ap/emcXpvGqf35tha+v0xGcaJzKkSERERaR+0YpyIiIj4FYUbERER8SsKNyIiIuJXFG5ERETEryjc+ImXX36ZlJQUAgICSEtLY8WKFb4uqU1Yvnw548ePJykpCZPJxIcffujrktqMGTNmMHToUEJDQ4mLi+PKK69k27Ztvi6rTZg9ezYDBgyoXWBs+PDhfPrpp74uq02aMWMGJpOJ+++/39eltAmPPfYYJpOpzpGQkODrstqMzMxMfvOb3xAdHU1QUBCDBg1i9erVzf48Cjd+YOHChdx///08/PDDrF27llGjRjFu3Lg60+hPV2VlZQwcOJAXX3zR16W0OcuWLeOuu+7i+++/Z8mSJbhcLsaOHUtZWZmvS/O5jh078pe//IVVq1axatUqzj//fK644go2bdrk69LalJ9++olXXnmFAQMG+LqUNqVv375kZWXVHhs2bPB1SW1CQUEBI0eOxGaz8emnn7J582aee+65Ftk1QFPB/cCZZ57JGWecwezZs2vPpaamcuWVVzJjxgwfVta2mEwmPvjgA6688kpfl9ImHTp0iLi4OJYtW8Y555zj63LanKioKJ555hkmTZrk61LahNLSUs444wxefvllnnzySQYNGsSsWbN8XZbPPfbYY3z44YesW7fO16W0OdOmTePbb79tlZ4Ftdy0c06nk9WrVzN27Ng658eOHcvKlSt9VJW0R0VFRYD3Q1yOcLvdvPPOO5SVlTF8+HBfl9Nm3HXXXVx66aWMGTPG16W0OTt27CApKYmUlBR+/etfs3v3bl+X1CZ8/PHHDBkyhGuvvZa4uDgGDx7Mq6++2iLPpXDTzuXm5uJ2u+ttNhofH19vk1GRxhiGwdSpUzn77LPp16+fr8tpEzZs2EBISAgOh4PJkyfzwQcf0KdPH1+X1Sa88847rFmzRi3DDTjzzDOZP38+ixcv5tVXXyU7O5sRI0aQl5fn69J8bvfu3cyePZsePXqwePFiJk+ezL333sv8+fOb/bl8uv2CNB+TyVTne8Mw6p0Taczdd9/Nzz//zDfffOPrUtqMXr16sW7dOgoLC3nvvfeYOHEiy5YtO+0DTkZGBvfddx+ff/45AQEBvi6nzRk3blztn/v378/w4cPp1q0bb775JlOnTvVhZb7n8XgYMmQITz31FACDBw9m06ZNzJ49u84+ks1BLTftXExMDBaLpV4rTU5OTr3WHJGG3HPPPXz88cd8/fXXdOzY0dfltBl2u53u3bszZMgQZsyYwcCBA3n++ed9XZbPrV69mpycHNLS0rBarVitVpYtW8YLL7yA1WrF7Xb7usQ2JTg4mP79+7Njxw5fl+JziYmJ9X45SE1NbZHJLwo37ZzdbictLY0lS5bUOb9kyRJGjBjho6qkPTAMg7vvvpv333+fr776ipSUFF+X1KYZhkFVVZWvy/C5Cy64gA0bNrBu3braY8iQIdx4442sW7cOi8Xi6xLblKqqKrZs2UJiYqKvS/G5kSNH1ltuYvv27bWbZTcndUv5galTp3LTTTcxZMgQhg8fziuvvEJ6ejqTJ0/2dWk+V1pays6dO2u/37NnD+vWrSMqKopOnTr5sDLfu+uuu3j77bf56KOPCA0NrW39Cw8PJzAw0MfV+dYf/vAHxo0bR3JyMiUlJbzzzjssXbqUzz77zNel+VxoaGi9cVnBwcFER0drvBbw4IMPMn78eDp16kROTg5PPvkkxcXFTJw40del+dwDDzzAiBEjeOqpp7juuuv48ccfeeWVV3jllVea/8kM8QsvvfSS0blzZ8NutxtnnHGGsWzZMl+X1CZ8/fXXBlDvmDhxoq9L87mG3hfAeOONN3xdms/ddttttf+eYmNjjQsuuMD4/PPPfV1Wm3Xuueca9913n6/LaBMmTJhgJCYmGjabzUhKSjKuvvpqY9OmTb4uq8345JNPjH79+hkOh8Po3bu38corr7TI82idGxEREfErGnMjIiIifkXhRkRERPyKwo2IiIj4FYUbERER8SsKNyIiIuJXFG5ERETEryjciIiIiF9RuBERERG/onAjIiIifkXhRkRERPyKwo2IiIj4FYUbERER8Sv/H8Fl68IDd5UwAAAAAElFTkSuQmCC\n",
"text/plain": [
"<Figure size 640x480 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"sns.histplot(\n",
" ds.hus.sel(level = 500), kde=True,\n",
" stat=\"density\", kde_kws=dict(cut=3), label = 'all'\n",
")\n",
"sns.histplot(\n",
" ds_above_26.hus.sel(level = 500), kde=True,\n",
" stat=\"density\", kde_kws=dict(cut=3), label = 'above 26'\n",
")\n",
"\n",
"plt.legend()"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "3d84ade2-6f0b-4db6-ab00-51d74f6d67da",
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Pangeo",
"language": "python",
"name": "pangeo"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.10.9"
}
},
"nbformat": 4,
"nbformat_minor": 5
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment