Skip to content

Instantly share code, notes, and snippets.

Show Gist options
  • Save gustavo-marques/84b402e9fe48c857976d9a6851f8ff54 to your computer and use it in GitHub Desktop.
Save gustavo-marques/84b402e9fe48c857976d9a6851f8ff54 to your computer and use it in GitHub Desktop.
barotropic_streamfunction_aqua_ridge.ipynb
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Basemap module not found. Some regional plots may not function properly\n"
]
}
],
"source": [
"from mom6_tools.MOM6grid import MOM6grid\n",
"from mom6_tools.m6plot import xyplot\n",
"import xarray as xr\n",
"import matplotlib\n",
"import numpy as np\n",
"import matplotlib.pyplot as plt\n",
"import warnings\n",
"warnings.filterwarnings(\"ignore\")\n",
"%matplotlib inline"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"### AQUA"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"MOM6 grid successfully loaded... \n",
"\n"
]
}
],
"source": [
"grd_aqua = MOM6grid('/glade/scratch/xiaoning/BMOM_aqua_production_branch_3h/run/BMOM_aqua_production_branch_3h.mom6.static.nc')"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {},
"outputs": [],
"source": [
"ds1 = xr.open_mfdataset('/glade/scratch/xiaoning/BMOM_aqua_production_branch_3h/run/BMOM_aqua_production_branch_3h.mom6.hm_04??.nc')"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {},
"outputs": [],
"source": [
"psi_aqua = -ds1.uhbt.cumsum(dim='yh')"
]
},
{
"cell_type": "code",
"execution_count": 16,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"array(8.880302e+08, dtype=float32)"
]
},
"execution_count": 16,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"psi_aqua.mean('time').mean().values"
]
},
{
"cell_type": "code",
"execution_count": 18,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAvgAAAF8CAYAAAC668+wAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nOzdeZhkVX3/8fdngBGRXSBEBBFREY1BHdQYN8QFUCFRECQmLlGiiSgajbhgBFwhigsxBHHj54IrRFFEkc2FqCggCig74j5hAEeWweH7++PeZmpqqvfq7qqa9+t5+um+555769TtMz3fe+p7zk1VIUmSJGk0LFroBkiSJEnqHwN8SZIkaYQY4EuSJEkjxABfkiRJGiEG+JIkSdIIMcCXJEmSRogB/gJKsm2Ss5JcmuSnSV7Ztf/gJD9r9x3VUf7QJOe15RcnWb8tf267/eMkX0uyRY/XfG2SC9uvnyRZmWTzuX+3kiRpPiXZJMmXk1zUxgwv7Nj3rjYO+EmS/Xsc+4Eky8c57+IkH21jjouSPLFj39uS/GK8Y6fR9u2T/GQa9Xu+bpJXJ7mkjY2+meQ+Hfu2S/L1Ng67JMn2XceOew0GnQH+wvoT8K9V9SDg0cC/JNkZIMluwD7AQ6vqwcB/tOXrAp8AXtqWPxG4oy1/H7BbVT0U+DHw8u4XrKqjq2qXqtoFeD1wTlXdMMfvU5Ikzb9/AS6pqr+kiRfe3QbnTwceDuwCPAp4bZKNxw5KsgTYdILzvgSgqv4CeEp73rGY8svAI/v9RqZgvNe9AFjSxkafB47q2HcicHQbhz0S+N3Yjilcg4FmgN9q7xQvS3JCezf7ySRPTvKdJJcneWT79d0kF7TfH9ge++okH2l//ov2+A0me82q+nVV/aj9+Q/ApcA27e6XAe+sqtvb/WOd7qnAj6vqorb8/6pqJZD26x5JAmwM/GqSJjwX+PSUL5IkSRrXFGOJeyT5SJIftPHEPh3HfivJj9qvx7TlT0xydpLPt+f+ZPv//FQUsFFbf0PgBprBxZ1pBvj+VFV/BC4C9mhfbx3gaODfJjjvzsA34a745EZgSbv9v1X160mu0zpJPtZeo4uTvKotf0T7icB5NDcnUzbe61bVWVV1S7v5v8C929faGVi3qr7R1ls+Vm+K12CgGeCvbkeaUfCHAjsBBwKPBV4DvAG4DHh8VT0MeDPw9va49wI7Jvlb4KPAP1XVLUl260iH6fz6bvcLtx8LPQz4Xlv0AOBxSb6X5Jwku3aUV5LT2z8A/wZQVXfQ3BRcTBPY7wx8eLw32t6A7AF8YfqXSZIkjWOyWOKNwJlVtSuwG3B0knvQjB4/paoeDuwPvL/jnA8DDqH5v30H4K8BkhwzTpxxaHvcscCDaOKCi4FXVtWdNAH9nkk2SJPOuxuwbXvMy4EvTRKkXwTsk2TdJPcFHtFx/FTsAmxTVQ9pPwX4aFv+UeAVVfVXnZWTPHCc93lhkumMsv8jcFr78wOAG5N8sb3ROroN7GFq12CgrbvQDRgwV1fVxQBJfgp8s6oqycXA9sAmwMeT3J/mrng9gKq6M8kLaNJi/ruqvtOWn0XTiSeUZEOaQPuQqrq5LV4X2IwmdWdX4LNJdmjLH9uW3QJ8M8kPgXNpAvyHAVcBH6BJwXnrOC/7TOA7pudIktRXk8US9wb2TvKatv76wHY0QfixSXYBVtIEoGO+X1XXt+e8sD3Pt6vqVZO05WnAhcCTgPsB30jyrar6ejtw+F3g98B5wJ+S3AvYjyadZyIfoblxOB+4tj3PnyY5ptNVwA5JPgB8Bfh6kk2ATavqnLbO/wP2BKiqnzGFeGoiSZ5H8ynDE9qidYHH0cRN1wGfAV6Q5DSmdg0GmiP4q7u94+c7O7bvpOkIRwJnVdVDaALk9Tvq3x9YDtxrrGAqI/hJ1qMJ7j9ZVV/sON/1wBer8f22DVu05edU1dL2o6SvsiqPjqq6sqoK+CzwmAne6wGYniP1VZK3JKl0TdQaZO3H5DXFuk9s398L5rhZ0jCbLJYI8Oyx+XBVtV1VXQq8Cvgt8Jc0gejicc65sj3PVEbwX8iqWOIK4GqaTxWoqre1r/+Utk2X0wS7OwJXJLkG2CDJFd1vsE3teVV7/D40ueqXT/UCVdWy9n2eTZOKc0Lbhp5/i2Y7gp/kyTSfnOw9lvpME09dUFVXVdWfgFNo4qkpXYNB5wj+9GwC/LL9+QVjhe1d5/uAx9Pcfe9bVZ+fbAS/zYn7MHBpVb2na/cpNHfcZyd5AM0/9KXA6cC/tSk2K2juRI9p27Vzki2r6vc0k14uHed1N2mPe9403rskTUt7I7BpVb13odsiDZDTgYOTHNyO7D+sqi6giTGub7MCng+sM/FpYAoj+NcBuwPfSvJnwAOBq9pUlE2r6v+SPJQmnejrbaC79djBSZZX1Y7dJ21jkFTVH5M8BfhTVV0ypXffHL8FsKKqvpDkSuBjVXVjkpuSPLaqvg38Xcf7nPEIfpKHAf8N7NExnxHgB8BmHXHTk4Dzq+orTOEaDDpH8KfnKOAdSb7D6v/wjgE+WFU/p8nvemeSraZwvr8G/h54Used6F7tvo/QfHz1E+Ak4PntHfgy4D00HfNC4EdV9ZWq+hVwOHBukh/T/EN4O0CSlyZ5acfr/i3NP+Q/zugqSBrPW4G703xkrWYg5JCFboQ0YI6kSfH9cft//JFt+QeB5yf5X5r0nH78H30k8Jg2PeibwOuqamn7+t9KcglwPPC8NrgfV5K9kxzRbm4F/CjJpcDraGKZsXpHJbmeZuT7+iRv6XG6bWgGMC8EPkaTUgzNJw7/mWaS7a3TeaMTvO7RNBOMP9fGWV8CaBcoeQ1NqvPFNJ8gfGg6rznI0mRzSJLWRkk+RjOAMOmqHGnWuj4LeGFVfWwK9c8Gtq+q7WfVSEnStDiCL2mtNFG+fJJr2uC0s+zpaVa0Wprk1iTXpVl94QEdddY4Z0fZA5O8vR1Zuj3NUnB70SXNqhbvSfLrJLck+d8ku081Vz7NUnvVa9Rskve8SZL/SvK7JLelWdbvURO8zsFJft7W/XmSg7v2X0OTCnif9jXHvp442XuQJM2OOfiSNIkkTwC+RLPM3Dto1ny+F/BkmslYP5/CaT4O3EHz0LrFNKkrpyR5QFVd01Hvc8BeNPNwzgDuC5xMMzluLp1Os5rGEcA9gVcDX02yffucjk4H0+So/jfwB5pnarw/yeZVdXhb5xCaa7UFzeTBMT3nBkmS+scAX5Imtw/NJ55P7ZqkdeQ49XtZCjyzXeWKJGcB3wf+iTb/tB3R3ws4oapeMnZgkjNplpKbSz+qqn/ueM1LaFbjOpAmkO/0AOBBHcv2/SfwbeBNST5cVddX1SlJDgHuXlWfmOO2S5I6mKIzDe0Ek0Mnrznt874rzdPcfpJk/47ylye5ov1Ye4uO8s2SnJzkx0m+n+Qh45z3SWkehvWTJB9P4g2dNDM3td+fPYt/R++rjklPVfUDmtHv+3fUeWb7fbVVtarqq8z9yPcxXdtntt/v312RZlnf68c2qmpFe/y6rHoPkhZYkv2S/DTJnUmWdJRv36Yaji3wcVyPY7/UTgIe296lTRm8MMn5SR7Z45jdsvoSlrcl+Zu5e4cajwH+NFTVl6rqnf08Z5Kns2od+0cBr02ycbv7OzQpAN0rcrwBuLCqHgr8A80Snd3nXUSTEnBAu27/tcDz+9l2aS1yLHABzSoXNyT5apJXJNlyGue4qkfZDTTpMGPuS7NWdq81l382jdeaidXaV1X/1/54zx51e91sjC2Rt0M/GyWNsjTmMhb7CfAsmodhdruyYy3+zpX2SPIsmmf7dDoKOLyqdgHe3G6vpqrOGjsnzbKTtwBf78P70DQZ4Lfau9nLkpzQjnh/MsmT24lmlyd5ZJIXJDm2rf+xJO9P8t0kVyXZd4YvvTPNg6v+1C5beRGwB0BVXdCVm9t5zDfbOpcB26dZ37bTPYHb26U7Ab4BPHuGbZRG0UQTVlcbpW+D3bHHyn8A2IhmxPrnSf5qzcN7WjlOeXr8PJvlzab8vlY7qFkyrpdeq+v0eo1JV+GRdFe8cWmSDwI/Ala2n+T/MMkZbbxxdhtb7N0e8+D2E/sL20/ve32ytoaqurRdQ3467duQZg7OW7tPB4wNQG5C8+TdiewLnNY+lFPzzAB/dTvSjIY/lOZJbwcCj6VZJ/UNPer/ebv/GcA7AZJslPGftrZzj3NcBOyZZuWMLWgCiG0naedFNHfktB+R3Yfm0dedlgLrdXwkt+8UziutTW5ov2/eWZhkfZp/26upqpVVdXZVvbGqxh5vviHwpj626Wqav8u9/vN+4BTP0fN9tfo1ut7rb9mD2u+dnwS4DrPU2wOBE6vqYe322VX1CJq0vbfSPKzyb2kmvQO8lCbNbxeap9yOzX/51jjxxpOn0Ib7Jrkgzepgj+soPxJ4N83oe6dDgKOT/IJmsYDXM7EDgE9PoR2aA+Zkr+7qqroYIMlPgW+2T5m7GNi+R/1TqupO4JKxEfR2tYkpP22tqr6eZFfguzQrWJwHTPiwCZqbifeleUDExTSpA6sd07b7AOCYJHej+YhssvNKa5OxT7eeTDOKNuZVdA1+JNmifThMp8toHsTSK5CeqS/TTLp9FXBQx+vvxaoAekJV9Yckv6F5gF46JvXuAPQrF/bvkry1Y5Lt4rbNK4FTO+otp3lS5F3tkATAtVX1v+3PK4CvtT9fTPPp+x1dscd5wBuT3Bv4YlVdDtAONszEr4Ht2ifZPoJmRa8H0wwC7FhVr8qay+m+DHhV+/TZ5wAfpvn7uYYkfw78Bc3qXFoABviru73j5zs7tu+k97XqrB9oRvCBb41z/gNpPtofW5HizW1e/9uAt7XHfwq4fKJGVtXNNE97I0loRv3WWEKvqs4DHtfWeyrNyheSGmfQBOlHJLknzb+hxwKPpvkErNOH2v9Yv04zn+XuwP40/55P7GObvkrzH+JL2k/0xpbJPAj4Mc2ni1NxLM0o4GlJTqFZ0vOlNPm4u/ahnT8HvtdOzPsDzd+2XYEjq+oXHfX+l+YTzmOTfJfmBuDMrpWIpLVR51Nq7+i4Ab4r9qiqO9NO6q+qTyX5HvB04PQkL66qM5N8i+bvULfXVNUZ4714Vd3e8To/THIlTYywK/CINM+xWBfYKsnZVfVEmnl8r2xP8TnghAne33OAk6vqjgnqaA4Z4PfZFEfw79qfZB1g0/Yu+qE0/4FPOCElyabALe3KFS8Gzm2D/u56W1XV79oR/NfR3kRIalJukuwDvJ9mXfcVNP/2nkAzwb3T/wNeQPMf3JbAzTSTSvetqi/0sU2V5Nk0/1afC+xJE9j/LfDP9E7d6eVdNDmyfw88sW3rPwKPoD8B/gdocnEPBrYDrgMOqaruCf/vpRkR3JfmBmMRTRqiAb40De0ncFdV1fvbnx9Kc7M8oxH8doGAG9q/gzvQ/G25qqrOB/6rrbM9cGob3EOTc/8E4GyaCbQTDUY+l8lTeDSH4qemjY6O/JB2+2Pt9ufH9tHknC2pqpd37m/rL6+qDWfwuuuzKj3gZuClVXVhu+8VwL/RPFDmd8BXq+rF7aS+E2lGwy4B/rGqlrXHfBV4cVX9KsnRNKNni4D/qqr3Trd9kgZD+3H9elW100K3RdLM9Yg37oof0jyBenlV/UfnviSvB55H87C83wAHVtUNPU7f/Vp/S3NDviXNA/ourKqntQMJR9Ck7q4E/r2qvjxJOx9LM09xXeA24J/b0f8lNLHLizuO+w6wbZvGrAVggC9JAyTJ3avq1q6yp9MMMhxbVQcvTMskScPCAF+SBkiSd9Cs0HMWzQO2dgFeRPMJ3y6dD5iSJKkXA3xJGiDtijmH0ixFuQnNspdnAodVVa8HYEmStBoDfEmSJGmEuIoOsMnGqa23WuhWSJIkjaafX8nSqtpyOsdsl9RtM3y938PpVbXHDA8fegb4wNZbwX+9Z6FbIUmSNJp234drp3vMbcB+68zs9T64ki1mduRoWDR5FUmSJEnDwgBfkiRJGiEG+JIkSdIIMcCXJEmSRogBviRJkjRCDPAlSZKkEWKAL0mSJI0QA3xJkiRphPigK0mSJA2cxYth+61nePB1fW3K0HEEX5IkSRohBviSJEnSCDHAlyRJkkaIOfjAnSvgtklytdbfbn7aovk12e99Pti3Rs8g9Cuwb42iQehb9qvRNAh9S/3jCL4kSZI0QhzBn6Kp3tk6sjEYhmkkwr41PIapX8H47bUvDZ5h6lv+zRouw9S31D+O4EuSJEkjxBH8PpvpnfKojHQ4UjB37FsL3YLR0Y9rab/SZPybtdAt0NrMEXxJkiRphDiCPyDm8k5/rkdDHKUYbP3+/czX6Jr9arANa78C+9ags29pzLqLYYuZ/v7W8t+FAf5awD846qfZ9qfu/2ztn4L+p3PYrzRmJn1hopsC+5aGgSk6kiRJ0ghxBF/SvHL0S/1kf9JcsF9p2DmCL0mSJI0QA3xJkiStVZLsl+RLSX6ZZHmSHyZ57gT1D0lSST7fY982SU5uz7M0ybFJNuhR7yVJLk9yW/t6u/f7fY0xwJckSdLa5tXAcuBVwN7AWcCnkhzcXTHJVsCbgd/32LcucDpwH2B/4JXAfsDxXfUOAI4DTgT2BH4KnJrkIf17S6uYgy9JkqS1zTOramnH9plJ7kUT+H+gq+47gK8A2/Y4z37Ag4Adq+pqgCR3ACclObyqLm/rHQ58vKqObOucAzwMOBR4Xp/e010cwZckSdJapSu4H3MBsFVnQZJdgefQBOK97An8YCy4b50CrAD2aM+xA/AA4LMdr38n8Ln2+L4zwJckSZLgMcAlYxtJAhwLHFVVvxznmJ2AyzoLqmoFcGW7j47vq9UDLgU2T7LlLNu9BlN0JEmSNHAWLZ7Vk4i3SHJ+x/bxVXX8eJXbCa/7AC/qKH4hsDXwHxO8zmbAjT3Kl7X76PjeXW9Zx/418vtnYyhG8JOsm+TQdubx7UmuT3JMV50keUOSXyS5Ncm5SXZZqDZLkiRpwSytqiUdXxMF99sDnwL+p6o+1pZtArwdeG1V3TrJa1Wv0/Yo797OBMfPyrCM4H8U2J1mgsJlNJMcdu6qcyhwGPDats6rgTOSPKSqfjOPbZUkSdIQSLI5cBpwHatPdn0D8Avg60k2bcvWBdZrt/9QVStpRuE3ZU2bsmrEfllH2U1ddaD3JwCzMvABfpI9gAOAv6yqS8apsz5NgP+Oqjq2LTsPuAZ4OfCm+Wnt8Ll9vIwy4G7bzF87NNo6+5n9SrPh3yzNFfvW2qddq/5UYDHw9Kr6Y8fuBwJLWBWcd1oGPA74Ns2g8k6dO5MsBnagWRYTVuXe7wRc21F1J+CGqupreg4MR4rOi4AzxwvuW48BNmb12cl/BL7MHM1OliRJ0nBq16//HHB/YM+q+l1XlTcBu3V9XQSc2/58cVvvNGDXJPfpOHZv4G7A1wCq6irg5zRLao69/qJ2+7S+vrHWwI/gA48CvpTkWOAfaNr8NeDlVfWrts5OwErg8q5jL6V56MBAmmi0YBDMpH2Ociy8UexXYN8aBKPYt+xXg2EU+1Yn+9lA+iCwF82DqTZP8uiOfRdU1U+6D0hyI01u/9kdxZ8H3gh8MclhwCbAMcCnOtbAB3gL8Ikk1wDfAZ5Pc3NxYL/eUKdhCPC3Bl5Ac9d0ALARcBRwcpJHV1XRzD5e3uZCdVoGbJBkcbtk0V2SHAQcBLDlZkiSJGnt8dT2+/t67LsvTZr3pKrqjjad/FiaTJLbgZNo5oR21vt0kg2B19HMGf0p8IxeNxL9MAwBftqvfarq/wCS/Bo4B3gS8M223ngzmHvua2dTHw9w/+3Sc/byoI8oDKJhu2YLNaoybNdpEAzTNbNfDY9hu2b2reE0bNdvbfjEoaq2n8ExTxyn/Hrgb6Zw/IeAD033dWdiGHLwlwEXjwX3rW/TPCFs5446GyVZp+vYTYFbquqOuW+mJEmStPCGYQT/UpqJCt0C3Nn+fBmwDrAj8LOOOms8XayXumP47q7VH/7eNRfsV5or9i3NB/vZ8BuGAP9U4PAkW1TV0rbs8cB6NHn5AN8FbqaZjfxWuGvpo2fSpuFIkiRpeMzySbZrtWEI8I8HXgF8OcnbaSbZvgs4o6q+DVBVtyV5J3BYkmWsetDVIuADC9NsSZIkaf4NfIBfVTcneRLwfppZySuA/wFe1VX1nTQB/euBewLnA0+pqt/OY3MlSZKkBTXwAT5AVV1Bs1bpRHUKeFv7JUmSJK2VhmEVHUmSJElTZIAvSZIkjZChSNHRKrddt9At6I9BnBU/k2s7iO9jpuxbc2dt7lv2q7kz02s7iO9lJuxbc2dt/ps1KhzBlyRJkkaII/jAnSv6NxIw0R3sqIw29MNUr8VMRwTm61rP5+/UvjU1c9m35vM6D0Lfsl+tMtej5aPYt/ybNTX+f6i54Ai+JEmSNEIcwe8z72D7y+u5iteiv7yeq3gt5s7afG3X5vc+F9bG65n14G7bLHQrhpMj+JIkSdIIMcCXJEmSRogBviRJkjRCDPAlSZKkEWKAL0mSJI0QV9HRguvHygA+QU+92LfUL3Oxgol9S93W9qcTq38cwZckSZJGiCP49PdJtvNpUO7YB+HaDUIbRslC9a1B/D0OYpuGmX1rlUFs07Ba2/8/tC+pmyP4kiRJ0ghxBH+IeceuuWLf0lyxb2ku2K9G06LFg/PpzLBxBF+SJEkaIQb4kiRJ0ggxwJckSZJGiAG+JEmSNEIM8CVJkrTWSbJjkv9OclGSlUnOHqfeXyQ5NclNSf6Q5PtJHtFVZ+ck30xyS5JfJTkiyTpddZLkDUl+keTWJOcm2WUu3psBviRJktZGDwb2An7efq2hDcC/C9wI7A/sB3wZuHtHnc2AM4AC9gGOAP4VOLzrdIcChwHvAp4JLAfOSLJ1395Ry2UyJUmStDb6clX9D0CSzwNb9KhzXFvveR1lX+uq81KagP9ZVXUz8I0kGwNvSXJUVd2cZH2aAP8dVXVs+5rnAdcALwfe1Mf35Qi+JEmS1j5VdedE+5PsDDwK+MAkp9oTOL0N7secRBP0P6HdfgywMfDZjtf/I82nAXtOr+WTG7oAP8k2SZYnqSQbdpTPW16TJEmSRt6j2u+btXn6f0pyZZJ/7Kq3E3BZZ0FVXQfc0u4bq7MSuLzr2Es76vTNMKboHE2Ts3SPrvKxvKbX0lzkV9PkNT2kqn4zv03UfFja5ycXbuHT8tSyb2ku9KNf2ZfUy6j+zcrsnmS7RZLzO7aPr6rjp3mOsdz4E4GjgB8A+wInJPl1VX213b8ZTY5+t2XtvrE6y6tqZY86GyRZXFUrptm+cQ1VgJ/kccAewNtpAv2x8lnlNf1pRf//cWj42Ac0V+xb6hf7kubDiPSzpVW1ZJbnGMt0OaGqjmp/PivJg4DXA1/tqFs9jk9X+Xh1xts3Y0OTotMuNfQBmpnJS7t2z2tekyRJkkbeDe33s7rKzwR27theBmza4/hNWDWyvwzYqHvpzPa4W6rqjlm2dTVDE+DTzFBeH/jPHvvmNa9JkiRJI+/SccoDdE7QvYyueDPJtjTp5Jd11FkH2LHrXGvk7/fDUAT4Se4JHAm8epw7nEnzmnqc86Ak5yc5/4/dR0mSJGlt912aWHL3rvLdgYs6tk8DnpZko46y/YFbgXM6znUzzTr6ACTZgGY9/NP62+zhycF/G/C9jskMvUwrr6mdaHE8wL3vlr7mPUmSJGmwtQH2Xu3mNsDGSfZtt79aVbckOQI4KsmNNJNsnw08nlXLX0KzVv4rgC8meRewA/AW4D1jS2dW1W1J3gkclmQZqxaEWcTky3BO28AH+EkeDLwIeHySsfymDdrvmyRZSUdeU9co/pzkNUmSJGnobQV8rqtsbPu+wDVV9d4ki4CDaYL2nwH7VtW3xg6oqmVJdgeOpZn/eSNwTFu/0ztpAvrXA/cEzgeeUlW/7eN7AoYgwAfuD6wHnNdj3/XAh4FPsSqv6Wcd++ckr0mSJEnDraquYVW2x0T13gO8Z5I6lwBPmqRO0WSlvG3qrZyZYQjwvw3s1lW2B/A6mo9VrgKuZVVe01thtbym6a55KkmSJA2tgQ/wq2opcHZnWZLt2x+/VVXL27J5y2uSJEmSBtXAB/jTMG95TZIkSZpbi2b3JNu12lAsk9mtqj5WVRkbvW/LqqreVlX3rqq7V9XjquqChWynJEmSNN+GMsCXJEmS1JsBviRJkjRCDPAlSZKkEWKAL0mSJI0QA3xJkiRphBjgS5IkSSPEAF+SJEkaIQb4kiRJ0ggZpSfZSpIkaUQsWndd1t98ixke/Zu+tmXYGOBr4GzR8VjqpdctXDskSZKGkSk6kiRJ0ghxBB9Yd/Hqo8YaHP5eJEkaAX4iP68cwZckSZJGiAG+JEmSNEIM8CVJkqQRYoAvSZIkjRADfEmSJGmEGOBLkiRJI8RlMiVJkjRwwobcjUfP8OhT+tqWYeMIviRJkjRCHMEHFi2G9X2gkiRJkkaAI/iSJEnSCDHAlyRJ0lonyQFJfpRkeZJfJjkxyb266iTJG5L8IsmtSc5NskuPc+2c5JtJbknyqyRHJFln/t7N6gzwJUmStFZJsjfwaeC7wD7A64DHA6cm6YyPDwUOA94FPBNYDpyRZOuOc20GnAFUe64jgH8FDp/7d9LbwAf4SfZL8qX2zmp5kh8meW6Pei9JcnmS29o6uy9EeyVJkjTwDgR+VFUvr6pvVtUngFcADwMeCJBkfZoA/x1VdWxVnQHsRxPIv7zjXC8F7g48q6q+UVXH0QT3r06y8fy9pVUGPsAHXk1zt/QqYG/gLOBTSQ4eq5DkAOA44ERgT+CnNHdgD5n/5kqSJGnArQfc1FV2Y/s97ffHABsDnx2rUFV/BL5ME2+O2RM4vapu7ig7iSbof0If2zxlw7CKzjOramnH9pltftSrgQ+0ZYcDH6+qIwGSnENzB3Yo8Lz5bKwkSZIG3keAU5L8A82i+VsDbwXOqqpL2jo7ASuBy7uOvRTYv2N7J+DMzgpVdV2SW9p9X+5/8yc28AF+V3A/5gKaHCeS7AA8AHhlxzF3JvlcZ9lEXCZTkiRppGyR5PyO7cB15LEAACAASURBVOOr6vixjar6SpIXAB8GPt4Wf5cmW2TMZsDyqlrZde5lwAZJFlfVirbejaxpWbtv3g18gD+OxwCdd1cAl3XVuRTYPMmWVfX7eWuZJEmSZm0RG7E+T5zh0acsraol4+1NshtNevf7gNOAPwPeApyc5MkdQX31OrzHvvHq9Sqfc0MX4LeTZ/cBXtQWjd0Zdd85LevYv0aAn+Qg4CCALRdk+oMkSZIWyLuBL1XV68YKklxIM2C8D/BFmlhyoyTrdI3ibwrcUlV3tNvL2rJum9B7ZH/ODcMk27sk2R74FPA/VfWxrt3dd0i97q5WVa46vqqWVNWSjTfoZyslSZI04HYCLuwsqKqfAbcC92uLLgPWAXbscWxn5shlrMooASDJtsA9WDPDZF4MTYCfZHOaj1CuY/WJs2Mj9d13TmPbC3LnJEmSpIF1LfDwzoIkD6JZ+eaatui7wM00S2OO1dmAZj380zoOPQ14WpKNOsr2p7lZOKffDZ+KSVN0kpw7xXPdVlVPnWV7xmvDBsCpwGLg6e0SRWPG7ox2ovll0bF9g/n3kiRJ6nIccEySX7EqB//NNMH9VwGq6rYk7wQOS7KMJuZ8Nc0A+Qe6zvUK4ItJ3gXsQJPP/56upTPnzVRy8HelWcB/IqGZpNB3SdYFPgfcH/jrqvpd5/6quirJz2nurk5vj1nUbp/GFGQ9uNs2fW22ZuH2X46/z9+TZqq7X9mXJGmt9n5gBfAymjj3RuDbwOu7BpLfSRPQvx64J3A+8JSq+u1Yhapa1s4RPZZmScwbgWNogvwFMZUA/7tV9fHJKiU5sA/t6eWDwF40S15unuTRHfsuqKrbaS7gJ5JcA3wHeD7NDcFctUlzyMBLc8F+JUkaU1UF/Ff7NVm9t7VfE9W7BHhS3xo4S5MG+FW1+1RONFfpOcDYeXt9QnBf4Jqq+nSSDYHXAYfRPMn2GVX1kzlqkyRJkjSQBn6ZzKrafor1PgR8aG5bI0mSJA22KQX4Sc5i4oX6a6oj/ZIkSZLmzlRH8D8xTvk2NLOGXUlekiRJfRPuzvr85UI3YyhNKcCvqg93bie5J81s4pcAnwGO6H/T5s+ixbD+dgvdCkmSJGn2pvWgqyQbJzkSuIJmvdCHV9VBVXX9nLROkiRJ0rRMKcBPcvckrweuAh4EPLaq/r6qrpzT1kmSJEmalqnm4F8NrAMcRbPA/58l+bPOClV1Zp/bJkmSJGmaphrg30azis7LxtlfNI/llSRJkrSApjrJdvs5bockSZKkPpjWJFtJkiRJg22qk2wfnuT0JGckedRcN0qSJEnSzEw1B/9EYG9gMXAyzUo6IyOugy9JkqQRMdUAfx2aibR3TuMYSZIkaUYWsZj12WahmzGUphqsPx84nial50Vz1xxJkiRJszHVVXS+DzxljtsiSZIkaZYmnWSb5MipnCjJ4bNvjiRJkqTZmMoI/iFJPgJkknqvAP599k2SJEmSNFNTCfDvAVzB5AH+bbNvzsJY5Co6kiRJGhGTBvhV5cOwJEmSpCFh8C5JkiSNEAN8SZIkaYQY4EuSJEkjxKfSSpIkaeAsYj2fZDtD0x7BT/KUJB9O8uV2e0mSJ/W/aZIkSZKma1oj+EkOBl4JnADs2xbfCrwfeEx/mzZ/Fq27LutvvsVCN0OSJGlE/WahG7BWme4I/iHAk6vqncCdbdllwAP72ipJkiRJMzLdAH8j4Bftz9V+Xw9Y0bcWzVCSnZN8M8ktSX6V5Igk6yx0uyRJkjR4Rjl2nG6Afy5waFfZK4Cz+tOcmUmyGXAGzU3HPsARwL8Chy9kuyRJkjR4Rj12nO4qOgcDX07yEmCjJD8Dbgae2feWTc9LgbsDz6qqm4FvJNkYeEuSo9oySZIkCUY8dpzWCH5V/RrYFXgOcCDwfOBRVbXQMyf2BE7v+mWcRPOLe8LCNEmSJEkDaqRjx2mvg19VBXy//RoUOwFndhZU1XVJbmn3fXmig8OG3I1Hz2HzJEmS1manLHQDus0qdhx0kwb4SY6Yyomq6s2zb86MbQbc2KN8WbtvDUkOAg4CuNd2d5+7lkmSJGm+bZHk/I7t46vq+I7taceOw2QqI/jbdvy8PvBs4AfAtcB2wCOBL/S/adNWPcoyTjntL/l4gL9YslnPOpIkSVoYd7KIFWww08OXVtWSSepMK3YcJpMG+FX1wrGfk5wEPLeqvtBR9ixgv7lp3pQtAzbtUb4Jve/OJEmStPYa6dhxustk7smaSVT/A+zVn+bM2GU0+VJ3SbItcI92nyRJkjRmpGPH6U6yvQL4F+D9HWX/DFzZtxbNzGnAa5NsVFV/aMv2B24Fzpns4EVsxPo8cQ6bJ0mStDYbuEm2s4odB910R/BfDLw6yfVJvpfkepqHAry4/02bluOA24EvJnlyO4H2LcB7hn0dU0mSJPXdSMeO0xrBr6oLktwfeDRwL+DXwHlVdcdcNG4a7VqWZHfgWJpljW4EjqH5RUmSJEl3GfXYcSbr4N8BfGsO2jIrVXUJ8KSFbockSZIG3yjHjtMK8CdaE3+B18GXJEmSxPRH8Lft2t6a5nG+J/enOQsj3J31+cuFboYkSZI0a9PNwX9hd1mSPYDn9q1FkiRJkmZs2jn4PXwd+EwfziNJkiQBsBK4YaEbMaSmm4O/Q1fRBsCBwC/61iJJkiRJMzaTB10VkHb7FuAC4B/62ShJkiRJMzPdB129rqrWqapF7deGVfU48DGwkiRJ0iCYboB/2Djlb5ptQyRJkiTN3pRSdJKMPQRg3SS7sSpFB2AH4A/9bth8WsRi1mebhW6GJEmSNGtTzcH/cPv9bsBHOsoL+A1wcD8bJUmSJGlmphTgV9V9AZKcWFVOqJUkSZIG1LRy8A3uJUmSpME26Qh+ksdX1bntz08ar15VndnPhkmSJEmavqmk6HwQeEj784fHqVM0k20lSZKkWfsTPsl2piYN8KvqIR2bH6yqo7vrJPnXvrZKkiRJ0oxM90m2hwFrBPjAG4F3z745C2MR67lMpiRJkkaC6+BLkiRJI8R18CVJkqQR4jr4kiRJ0ghxHXxJkiRphEx3ki1J/gx4JLAFHbn4VfWRcQ+SJEmSNC+mFeAn+RvgE8DlwIOBn9Kskf9tVs/NHyp3sogVbDDnr7OYW+b8NcbMx/tZCLO9hl6X2RvFa9iP6+d1mZ1RvH5g3xqPfWv2ZnINR/VaaE3TStEB3gq8sKoeBvyx/X4Q8MO+t0ySJEnStE03RWe7qvpcV9nHaVbSeU1/mjS6vHOePa9hb16X2fH69eZ1mT2vYW9el9lbG67hihVw3XUL3YrhNN0R/N+1OfgA1yT5K+B+wDr9bZYkSZKkmZhugP8h4LHtz8cAZwEXAR/sZ6PGJNk4yeFJvp/kpiS/SXJykgf0qLtJko8mWdbW/WSSe85FuyRJkqRBNd1lMt9VVV9ofz4ReADwaJoHXs2F7YCXAKcD+wL/BPw58L0k23bV/QzwRODFwAuAXYFT5qhdkiRJ0kCa7gj+aqrqOuBnwBv705w1XA3cr6oOq6pvVNX/AHsB6wEvGqvUpgo9DXh+VX2hqk4Gngc8NsmT56htkiRJWgu0WSXvTXJNkluSXJrkkCTpqrdNm22yPMnSJMcmWWPCRJKXJLk8yW1Jfphk9362d9rr4I8jk1eZvqr6Y4+yG5JcC2zVUbwn8NuqOrej3veTXN3uO2Oi11kJ3NCfJkuSJGn0fAx4PPAG4ApgN+A9NDHwMQBJ1qXJOlkB7A9s2tbZlGbgmbbeAcBxwFtolpp/IXBqkl2r6if9aGy/Avy5StFZQ5ItgR1ZPe9/J+CyHtUvbfdJkiRJ09aOwO8DHFJVx7fFZyZ5MHAAbYAP7Ac8CNixqq5uj70DOCnJ4VV1eVvvcODjVXVkW+cc4GHAoXTcCMzGlAL8JE+aYPfifjRkGt4NLAdO6ijbDLixR91lwA69TpLkIJo1/Nlmu+363ERJkiSNiHVp0tpv6iq/Ebh3x/aewA/GgvvWKTQj+nsAlyfZgWYO6yvHKlTVnUk+11nWjwZPxYcn2T/lVUqTbEIzUXZCVbXGiHySl9Hc2Ty7qv6v+5BeLzdOOe0d2PEAf7lkybx9AiFJkqThUVU3J/ks8G9JfgxcCTwBeA7wLx1VdwIu6Tp2RZIrWZVRMva9O869FNg8yZZV9fvZtnlKAX5V3Xe2L9RhP5rlNifTPWlhb+ADwOvaSbSdlgFb9jjHpvQe2ZckSdLo2iLJ+R3bx3ek18zEPwCfBC5otwt4fVV9vKPORBklm3XUoUe9ZR375yfA76eqOgE4YTrHJHkMTUrOcVV1dI8qlwGP61G+Ey6VKUmSNHRuvxWuuHjGhy+tqiXj7ZxBRskxwKNoJsReRfNcqLckWVpVnZkuU80o6d7OOOUzMu8B/nS1ExhOBb4GvGKcaqcBhyV5bFV9uz1uCU3+/Wnz0lBJkiQNiylnlCR5KPAy4KlV9Y22/NwkGwH/keSjVXUnzSj8pj3O0ZlRsqyj7KauOtCnzJOBDvCTbEUT2C8H3g88smO50Zur6hKAqjovyenAiUleA9wJvAv4dlVNuEQmwJ9wmUxJkqS1xTQzSsby5i/sKr+AJjC/J01azWV0rd6YZDHNgPNxbdHYJwI7Add2vcYN/ci/h1k+6Goe7EwzO3lb4CzgvI6vD3bVPQA4B/gIcCLwQ+Bv562lkiRJGkVjgfjDu8ofAfwRWNpunwbsmuQ+HXX2Bu5GM2BNVV0F/JzmEwQAkixqt/uWdTLQI/hVdTZTfIhWVd1Ikxf1wrlskyRJktYq57dfH0nyZuBqmhz8Q4D3VdVY3vzngTcCX0xyGLAJTe7+pzrWwIfmAVefSHIN8B3g+cD9gQP71eCBDvAlSZKkhVRVK5M8E3gr8GaalRuvpQnU391R744kewDHAp8FbqdZJOa1Xef7dJINgdcBhwE/BZ7Rr6fYggG+JEmSNKGq+g3w4inUux74mynU+xBTm+Q7I4Oegy9JkiRpGhzBB1asgOum/CxeDYPttpvd8fYH9TLTfmV/0mRm0rfsV+pltv//dbOfDScDfI0k/yBpLtivNFfsW+oX+5LAAF+SJEkD6Pbb4JrLJq+nNZmDL0mSJI0QA3xJkiRphBjgS5IkSSPEAF+SJEkaIQb4kiRJ0ggxwJckSZJGiMtkArffCldcvNCtkCRJkmbPEXxJkiRphBjgS5IkSSPEFB1JkiQNnNvvgCt+u9CtGE6O4EuSJEkjxABfkiRJGiEG+JIkSdIIMcCXJEmSRogBviRJkjRCDPAlSZKkEWKAL0mSJI0QA3xJkiRphBjgS5IkSSNkqJ5km+QQ4BjgC1W1b9e+bYBjgacAtwEnAf9WVbfMe0MlSZI0K7evhCtuWuhWDKehCfCTbAW8Gfh9j33rAqcDK4D9gU2B97TfnzePzZQkSZIW1NAE+MA7gK8A2/bYtx/wIGDHqroaIMkdwElJDq+qy+evmZIkSdLCGYoc/CS7As8BDh2nyp7AD8aC+9YpNCP6e8xx8yRJkqSBMfABfpLQ5NYfVVW/HKfaTsBlnQVVtQK4st0nSZIkzUiS/ZN8Mcmvk1SSF0xQ96AkP0lyW5LfJvlMjzovSXJ5W+eHSXbvUWebJCcnWZ5kaZJjk2wwlfYOfIAPvBDYGviPCepsBtzYo3xZu28N7cU/P8n5y29aI61fkiRJGrMvsD1w6kSVkrwVOAr4GPA04GBgaVedA4DjgBNpslB+Cpya5CEddcbml96HZn7pK2lS0o+fSmPnPQc/ySbAn09Wr6oua+u+HXhFVd062SG9Xm6ccqrqeNqLtN39l/SsI0mSJAH7V9WdSTYEXtyrQpIHA68H9qiqb3Ts+mxX1cOBj1fVke1x5wAPo0lFH1scZlbzSxdiku1+wIemUC/AG4BfAF9Psmlbvi6wXrv9h6paSTNSv2mPc2xK75F9SZIkaUqq6s4pVHs+cGVXcL+aJDsAD6AZkb/r3Ek+11nG5PNLJwzw5z1Fp6pOqKpM9tVWfyCwhCaAH/v6a2Dv9ue/autdRleufZLFwA505eZLkiRJc+BRwMVJ3tLmzN+e5IwkD+qoMxavdsenlwKbJ9myo96M55cOeg7+m4Ddur4uAs5tf764rXcasGuS+3QcuzdwN+Br89ZaSZIkDYItxuZatl8HzcNrbg08lWblx5fQ5O1vCXwtyfptnbG5od0ZJsu69k97fmmngV4Hv6p+0l2W5EZgaVWd3VH8eeCNwBeTHAZsQvPE20+5Br4kSdLwuW12T7JdWlVLxts5nTmh03jNRcA9gGdX1aXt6/wU+Dnwd8CHO0/d3aQe5dOaX9ppoAP8qaqqO5LsQbOc5meB24GTgNcuaMMkSZI0iKYzJ3SqbgA2HAvuAarqqiTXADu3RWMj9ZsCnbcvY3NJb+yoN+P5pUMX4FfVE8cpvx74m/ltjSRJkoZNVZ0AnNDn014KbNejPMDYJN2xTwR2Aq7tqLMTcENV/b6j3njzS4+brCGDnoMvSZIkDYNTga2TjI3Wk+R+NGvZXwTNiD5Nys5+HXUWtdundZxrVvNLh24Ef67t+BcL3QL1wxUXT15nPtmvRsOg9Suwb40K+5bmyiD2rWHUBu07A2OTZZckWQ78vqrOactOBn5EMyf0TcBK4AiagL7zabZvAT7Rpu58h2Z5zfsDB3bUmdX8UkfwJUmSpIk9B/gc8P/a7X9ptw8fq9A+m2kv4AKaCbUnAlfQPPjqjo56nwZeCryAZjT+ocAzOheXaevvQfM8qM/SzDP9AjCl1YBS5UNc7//QJfXeU89f6GZIkiSNpGfcJz+caFWbXrL5klr3qTOLz/70mem/3ihxBF+SJEkaIQb4kiRJ0ggxwJckSZJGiKvoSJIkaeCst/IOtrrplzM6dmZHjQ5H8CVJkqQRYoAvSZIkjRADfEmSJGmEmIO/Fthuu/6e77rr+ns+DS/7luaKfUtzwX6ltYUBPrB4cf//0Y+yYfsD6e92eNi35sdUr8uwvr9ehqlvjdJ179R5zUblPQ5Tv4LRue6anCk6kiRJ0ghxBF8LzhEFzRX7Vm9el9nzGk6f12xyXiP1iyP4kiRJ0ghxBF+SJEkDZ72VK7jXzb+a0bE+6EqSJEnSyDDAlyRJkkaIAb4kSZI0QszBn6HNp1jvhjltxdptot/BsF53+9VgGO/3MMzX3b618EbxbxbYtwZN9+/D6752cgRfkiRJGiGO4NNchKmOQExXP867tt199+OazdXvc1DMxfsb9X5mv5qafr/HUe9XYN+aKvvW9Nm3NFOO4EuSJEkjxBH8IeDdt+aD/UxzwX6luWLfksY3FCP4Se6T5NNJbkhyS5KLkuzRVWebJCcnWZ5kaZJjk2ywUG2WJEmSFsLAj+An2RY4D7gIeCHwR2AX4O4dddYFTgdWAPsDmwLvab8/b56bLEmSpFlavHIF29103YyO/UGf2zJsBj7AB44GrgSeXlV3tmVndNXZD3gQsGNVXQ2Q5A7gpCSHV9Xl89ZaSZIkaQENdIpOkk2AZwEf7Ajue9kT+MFYcN86hWZEf4/eh0iSJEmjZ6ADfODhwHpAJflOkjuSXJ/k9UnSUW8n4LLOA6tqBc3I/07z11xJkiRpYQ16gL91+/2/gW8BTwU+ArwVeFlHvc2AG3scv6zdt4YkByU5P8n5//f73/evxZIkSdICmvcAP8kmSXaa7KurfadV1aFVdVZVvRn4OPD6rlNXr5cbp5yqOr6qllTVkntuuWVf3pskSZJGS5KNkxye5PtJbkrym3blxgdMcMyiJD9MUkme0WP/PkkuTnJbkkuS7N+jziZJPppkWfu6n0xyz6m0eSFG8PcDLp3CF6x6UN1ZXec4E7h3ko3b7WU0K+Z025TeI/uSJEnSVGwHvIRmxcZ9gX8C/hz4XrvaYy8vBrbptSPJY4Ev0MS3ewJfAT6d5KldVT8DPLE91wuAXWnmmE5q3lfRqaoTgBOmWP3SccrH8u/HJt5eRleufZLFwA7AcdNtoyRJktS6GrhfVd06VpDkW8B1wIuAwzsrJ9kMeBtwKL1j3sOAc6vqFe32WUkeDLwZ+Hp7jr8CngY8oarObct+SXNT8eSq6l5RcjUDvUxmVV2T5KfA7jR5+GN2B66squXt9mnAgUnuU1XXtmV7A3cDvjZvDe6h30/au2HyKgNrPp86OMzXaSrm4loO6zWzX83eXF7DYb5m9q3ZmevrN8zXzL41XKrqjz3KbkhyLbBVj0OOBL4DfLN7R5K7AbsBr+jadRLw0SSbVNVNNCP7vx0L7tvX/H6Sq9t9Ewb4gz7JFpq7nH2THJ3kKUneBvw98O8ddT5PM4r/xSR7JXkucCzwKdfAlyRJUj8l2RLYEbikq/yhNA9mfc04h96PZoXIy7rKL6WJy8fy+tdYIbKj3qQrRA70CD5AVZ2c5B+ANwKvpPk45F+q6pMdde5IsgdNUP9Z4HaaO6HXzua15/MOe6pm0qZ+370P4nXpNtU2LtTIxiBeQ/vW5OxX0zfTNtm3erNvreLfrKkZ9L41kdk8yRbYIsn5HdvHV9XxfWjWmHcDy2nizU4fAP6zqq5Isn2P48ZWd+yeI7qsa/9EK0TuMFnjBj7AB6iqTwCfmKTO9cDfzOT86zAc/0hnapTf22x5bWbH69eb12X2vIa9eV1mx+s3vhG8Nkurasl4O9uHqf75ZCepqjVG0ZO8DHge8Oyq+r+O8gOABwLPnEL7uld5TI/yaa0Q2WkoAnxJkiSpj/YDPjSFelltI9mbZpT+dVV1ckf5esDRwLuARUk2BcZWe7xHko2q6g+sGqnvXv1xbHts1H4Z0Gsd9ymtEDkMOfiSJElS31TVCVWVyb46j0nyGJqUnOOq6uiuU94DuDfwHprgfBlwUbvvJOCC9ucrgTtYM49+J5rVIX/ebq+xQmRHvV65+asxwJckSZIm0C5jeSrN6ozdK+BAk4+/W9fXc9t9bwD+DqCqbqdZ/36/ruP3B85rV9CBZoXIrds188fasIQm//60ydprio4kSZI0jiRb0QT2y4H3A49M7hrcv7mqLqmqPwFndx23ffvjxVX1vY5dRwJnJ3kvzYOr9mq/9hirUFXnJTkdODHJa2hG998FfHuyNfDBAF+SJEmayM406TfQjL53OofmabNTVlXfTrIv8FbgZTQP0jqwqr7eVfUA4BjgIzRZN6fS+9ODNRjgS5IkSeOoqrPpmmw7xeOuGe+4qjqFZvR+ouNvpFlT/4XTfW1z8CVJkqQRYoAvSZIkjRBTdCRJkjR4Vq7gTzddvdCtGEqO4EuSJEkjxABfkiRJGiEG+JIkSdIIMcCXJEmSRogBviRJkjRCDPAlSZKkEWKAL0mSJI0QA3xJkiRphBjgS5IkSSPEJ9lKkiRp4Nx+5+1c84drFroZQ8kRfEmSJGmEGOBLkiRJI8QAX5IkSRohBviSJEnSCDHAlyRJkkaIq+jM0GJu6Vm+gg3muSUaduP1pYnYzzQV9i1JWjs5gi9JkiSNEEfwgUXcOaORrl76dR5pIvYzzRX7liQNP0fwJUmSpBGSqlroNiy4JH8AfrbQ7RghWwBLF7oRI8Jr2V9ez/7yevaP17K/vJ791Y/r+f/bu/tgu6ryjuPfX94IkZcIQjGUFykIIoPoACW2lSgMEoYXFfCCaM1My2hnbLSWQi0gEbFCO1WKlGaqI9Y3goRABzRE0xiJyhBEKBYMVUwCSGB4uUnMCwHC0z/WOrLu5tx7z7333Jy332dmz7137Wevs/aTvXfWWWftfQ6IiL1GsoGkO/Jrj8YzEXHyKLfteO7gA5J+FhFHt7od3cL5bB7nsrmcz+ZyPpvHuWwu57O5nM/O4yk6ZmZmZmZdxB18MzMzM7Mu4g5+8h+tbkCXcT6bx7lsLuezuZzP5nEum8v5bC7ns8N4Dr6ZmZmZWRfxCL6ZmZmZWRdxB9/MzMzMrIv0TAdf0nxJqyRtktQv6U5JJ9aJ21fSLTnuGUnXSppWJ+58Sb+S9LykeyWdsGP2pPUk7SbpM5JWStog6cmcszdW4g6UFHWWBXXqPEPSL3I+H5LUt+P2qHUazWWO3V3S9fn43SDpW5L2rBPXk7mskdQnaZGkdfl4m1MnZtYgx+aVdWJ79lyHxvKZ43ztHAVJywc5FqdW4hrKby+TdLik/5a0RdITki6XNLHV7Wp3kuYMcgx+tIiRpH+Q9JikrbkPdVQr221D65kOPrAzcC3wXuCDpC9sWCzpuFqApEnAEuAAoA/4OHA2lZtLJJ0DzAe+DswGHgRul3TE+O9GW9gfOJ+Uq7OAjwCvB+6WtF+d+AuAmcVySblS0p8CNwM/JOXzu8ANkk4arx1oIyPJ5Y3ALOAvgTnAMcCtZUCP57LmLOBA4PYGYs9j4LH5b+VKn+tAA/n0tXPMfsjA43AmsK22stH89jJJrwWWAgGcAVwO/C3wmVa2q8O8i4HH4KJi3d8DlwJXAacBm4ClkvbZ0Y20BkVETy7AROBR4Jqi7FxgO/CGouz9wMvAIUXZw8BXi78nAL8Avtnq/dpBuXsNsHOlbA/SCX9ZUXYg6WJ76jD1LQGWVcq+B/y41fvaRrmcmXP5jqLs2Fx2onM5YH8n5J+75PzMqRMzK687Ypi6evpcH0E+fe0cfX6XAwuHiWkov728AJ8C+oHdirILgS1lmZe6uZuTz+1dBlk/FdgAfLooew3wNHBFq9vvpf7SSyP4A0TEdmA9MKUong3cExGri7JbgReAkwEkHQS8EfhOUdfLwE15+64XEZsjYmul7DlgLbD3SOqStBPwTop8ZguAmZJ2H0tb290IcjkbeCoi7iziVgKr87qez2VNPh/HzOd60mA+fe0cX8Pm15gNLImIjUXZAtKn98e3pkld4+3Abgw8dzcDt+Fzt231VAc/zyGbJGlPSX8DHAJ8tQg5DFhVbhMRLwCP5HUUPwfEAb8E9pC0V/Nb3v7yfh8MPFRngNhMcAAACWFJREFU9fWStuc5vF+QtHOx7o+AydTP5wRSh6CnDJLLVx2b2S955Zh0LkduWT4210i6pDJf1+d643ztHJuT8rzxLZKWSDqysr6R/Pa6ejl6lDSC7xw15hFJL0l6WNJHivLDSJ8g/aoSX/7/Y21mUqsbsIP1ATfk3zcDfXkUtOa1pFH9qv68juJnNa6/WP/02Jvacf6FNK2kvIF2G2lO8/eBjaRpEReROqJn5JhG8tlr6uVyqGPzoCKGOnG9nMvBbACuBFaQRkFPJc3V3Ys0vxl8ro+Er52j9yPgP4Ffk+bYXwyskPSWiFiTYxrJb69zjkZvHWl+/UrS9OVzgfmSpkXEF0n525RnPpT6gWmSpuQ3nNZGOraDn6cbvH64uIgo39EvId2Y+DrSzXULJJ0SEcvLTeq9XJ3y6t8aYvu2N8p81rb9K9KNy2dGxLNF7DrgY0XocklPAddJOioi7i+rrlY7SHnbG49c1jap93J1yrsmlzC2fA4Rex9wX1G0VNI24JOSPhsRz5Th1SYNUt4RxiOftU3qvVyd8q7KZ9VI8xsRlxXFKyQtJY1EfyIvv9+k3ssNUt6rnKNRiIglpP5RzeI85fMSSf9aC6uzaVedu92mYzv4pCcIfLmBuNoBSET0Az/Lf94haQbpTvt35LJ+YHqdOqbzyshAf1G2oRID9UcQOsGI8wkg6XTgS8BFEXFLA9svBK4D3gbcz8B8ljo5n+ORy37S6HLVYMdmNQY6M5cwynyOwkLSTXlHAsvwuT6SfPbytbNqTPmNiCcl/YR0jaxpJL+9brAc7Y5zNBoLSTdyH0jK7a6SJlZG8acDWyLixRa0z4bRsXPwI+IrEaHhlmGquY9XpjdAGjUZMJ9M0pQcs6qIoRqX/34uIjryI+bR5FPS20nTSOZHxD83+lKVn48AL1I/ny8D/ze6PWqdccrlq47NrJx32nW5hKad6yN6yfzT53rjevbaWdXE/Jajoo3kt9fVy9F+pKe9OEejF6T8TSTdG1Ya7N4wawMd28EfK0kiPXqwfCrBYuAYSQcUZacDOwF3AETEb0gdpbOLuibkvxePc7PbhqQ3k56LfQcwdwSbnpV/3gsQEdtIz4A+uxLXB9wVERvocg3mcjGwT37OfW27o0n/wS8G57IJzgReAh4An+sj5Gtnk0j6A+BPyNfIbNj8GouBd0vatSjrA7aS7nOwkTmT9H1Ba4Gfku6jK8/daaTn4fvcbVOdPEWnYZL+jPSFF4tIz77fE/gwcBzpAK1ZSLrBaZGkS0kf7X0R+HZElHePzwO+KWkN8JNc1yHAB8Z1R9qEpL1J/6lsAq4Bjk3vlwDYGBEP5bh5wK6kHG0kTYX6O2BRRDxQVPlZ0vz8q0mPfjslL13/+LdGcxkRd0laAnxd0gWkEfmrSM+3X1pU2bO5rJF0OHA46dnNAEdL2gQ8HRE/yjH/Trqh8x7STbankO4Xubpy78M8evhch8byia+do5KflvN50qNC15K++O5TpPP76iK00fz2svmkAZJFkq4iDX7MA75QeXSmVUi6mXSD7QOkkfq+vMzNj7J9Xulbvi+V1E8atf8kaZD4S61ptQ0r2uBh/OO9kOaQLQQeJz3Z5XHSiOnMOrF/SOoYbQKeJT0FZlqduPNJTz3YBvwcOKHV+7kD8zmL9LFdvWV5EXcO6Z6HDaRO1K9J9zzsVKfO9wD/m/O5Cjin1fvZTrnMsdOB60nzSTcC3wZe51y+av/nNXBsziX9Z/a7nKcHSTc0TqhTX8+e643mM8f52jny3O5L+iK6dfka+Szpm6gPqxPbUH57eSG9EV1GGrVfRxrwmNjqdrX7Avwj6UvotuTc3Qt8qBIj0pvMx3PMCuCtrW67l8EX5X84MzMzMzPrAj07B9/MzMzMrBu5g29mZmZm1kXcwTczMzMz6yLu4JuZmZmZdRF38M3MzMzMuog7+GZmZmZmXcQdfDMzMzOzLuIOvpmZmZlZF3EH38yshSQ9KGlWk+r6mqQrRrltSNos6XNjbMMySc9L+vFY6jEzs9FzB9/MLJO0RtKJO/I1I+LNEbG8Va9f8ZaIuLj2h6Q/lrRK0kOSji7K10jaKmlTsVwLEBHvAj7agrabmVk2qdUNMDOztvU54NT8+3XAScW60yJi6Y5vkpmZDccj+GZmw5D0JknLJa3PU2pOr6xfI+kCSQ9I2iDpRklT87q3SbpP0u8k3ZTXXVHZ9kRJ3wD2B27LI+IX5vUh6eAi/vfTcCS9VdLPc903AlMr7Zoh6WZJT0taLWnuCHd9Yl4m4QEhM7OO4Q6+mdkQJE0GbgO+D+wN/DXwLUmHVkLfD5wMvAE4EpgjaQpwC/A1YA/gBuC99V4nIj4EPEoaGd8lIv5pmHZNAW4FvpHrvgk4s1g/Ibf7f4B9gROAT0h6d6P7DlwGfBf4L+DiYWLNzKxNuINvZja044BdgCsj4oWIWAbcDpxbibsmIp6IiOdIHeuj8raT8roXI2IRsLKJ7ZoMXJ3rXgjcU6w/BtgrIi7P7f4N8GXgnEZfICLujIiDI+LQiLirsvrW/IlGbTl/rDtkZmbN4Y9czcyGNgN4LCJeLsrWkkbFS08Wv2/J280AfhsRUax7rIntqta9tvj9AGCGpPVF2URgRZNe/z2eg29m1p48gm9mNrQngP3ylJea/YHfNrDtOmBfSSrK9hsiPuqUbQGmFX/vM0Td+xe/PwasjojpxbJrRJzSQLvNzKyDuYNvZjbQZElTawtwN7AZuFDS5PzM+tOABQ3UdRewHfiYpEmSzgCOHSL+KeCgStn9wAckTZR0MnB8UfdLwNxc9/sqda8ENkq6SNLOefsjJB3TQLvNzKyDuYNvZjbQ94CtxfJp4HRgNvAM6XGRfx4Rq4arKCJeAN4H/AWwHvggaf7+tkE2+TxwSZ7TfkEu+zjpDcV64DzSjbVl3XOAfqAPWFS89va83VHA6tz2rwC7D5+ChtSe9lNbbmlSvWZmNkYaOH3TzMzGk6S7gfkRcX2r21KS9Dzpjcc1EXHpGOr5AekG4JURcUKz2mdmZo1zB9/MbBxJOh54mDSCfh4wHzgoIta1tGFmZta1/BQdM7PxdSjwHdKjNh8BznLn3szMxpNH8M3MzMzMuohvsjUzMzMz6yLu4JuZmZmZdRF38M3MzMzMuog7+GZmZmZmXcQdfDMzMzOzLuIOvpmZmZlZF3EH38zMzMysi/w/309tyba2I54AAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 864x396 with 2 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"plt.rcParams.update({'font.size': 15})\n",
"fig, ax = plt.subplots(nrows=1, ncols=1, figsize=(12,5.5))\n",
"dummy = np.ma.masked_invalid(((psi_aqua.mean('time'))/1.0e6).values)\n",
"xyplot(dummy, grd_aqua.geolon, grd_aqua.geolat, grd_aqua.area_t,axis=ax, clim=(-3000, 3000),title='using uhbt')"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"### Ridge"
]
},
{
"cell_type": "code",
"execution_count": 20,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"MOM6 grid successfully loaded... \n",
"\n"
]
}
],
"source": [
"grd_ridge = MOM6grid('/glade/scratch/xiaoning/BMOM_brdg_production_branch_3h/run/BMOM_brdg_production_branch_3h.mom6.static.nc')"
]
},
{
"cell_type": "code",
"execution_count": 21,
"metadata": {},
"outputs": [],
"source": [
"ds2 = xr.open_mfdataset('/glade/scratch/xiaoning/BMOM_brdg_production_branch_3h/run/BMOM_brdg_production_branch_3h.mom6.hm_04??.nc')"
]
},
{
"cell_type": "code",
"execution_count": 22,
"metadata": {},
"outputs": [],
"source": [
"psi_ridge1 = ds2.vhbt.cumsum(dim='xh')"
]
},
{
"cell_type": "code",
"execution_count": 23,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAu8AAAF8CAYAAAB/mVUhAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nOzde9zcZJn/8c9VoFQUWkoRFKmIIFVRkC2uuqig4A9wAQ8grrKeZXVVPKwnRBAFFETFA8tiQQVWEBUFjwgiFEVBRQ6CUKUIFDzSpRQrlGJ7/f64M33ypDPzZDLJzJ3k+369+upMJpO5k5ln5so3d+6YuyMiIiIiIvGbNu4GiIiIiIhIPireRURERERqQsW7iIiIiEhNqHgXEREREakJFe8iIiIiIjWh4l1EREREpCZUvI+ZmX3RzP5qZjdmph9jZr82s+vM7GIze3QyfXczW55Mv87Mjuqx3DPM7LbUfDuPYn1ERERkdCz4rJktTuqGXXrMd5yZ3WlmKzLTn2Nm15jZP8zswMxjq1N1xLdT05+XPOdGMzvTzNZPpr8yacOvzexnZrbTEOu1TbY26jPvRmb2PTNbZGa/MbPjU4+9y8xuStr0IzN7bI9l/MDMrk+ef6qZrZdM/2pqG9xuZtcl0zczs8vMbIWZnVx0PYtQ8T5+ZwB7d5l+ors/1d13Br4LpIv0n7j7zsm/j/RZ9ntS811XYptFREQkDvsA2yf/DgX+p8d83wGe3mX6EuA1wDldHnsgVUfsD2Bm04AzgZe7+47AHcCrk/lvA57r7k8FjgEWFFqjYj7h7vOApwH/Ymb7JNOvBeYnbToP+HiP57/M3XcCdgQ2Bw4CcPeDO9sA+AbwzWT+lcCRwLsrWZs+VLwnkj28RWZ2erInebaZ7WlmPzWzW8zs6cm/n5nZtcn/OyTPfZeZfTG5/ZTk+RvleV13/zFwT5fp96XuPhzQ1bRERETGLGe98PDkyPovk5rhgNRzf5Kk1teY2bOS6bub2UIzOy9Z9tlmZjmbdABwlgdXAbPM7FHZmdz9Knf/U5fpt7v7r4E1OV9vM+BBd/9dcv+HwEuTZf3M3Zcl068CHtNtAWa2XtJD4EYzu8HM3plM/6ck/b4SeEvO9uDu97v7ZcntVcA1ndd298vc/f6p2pSqu9YHppOpu5L342XAV5L5/+7uVxCK+JFaf9QvGLntCHtahwK/BF4B7AbsD3wAeBXwHHf/h5ntCXyU8IH9NLDQzF4MHAH8h7vfb2Z7ACd1eZ373f1ZUzXGzI5LXnM5sEfqoWea2fXAH4F3u/tveiziuKRbzY+A97v7g1O9poiIiExpqnrhJuBSd3+dmc0CfmFmlwB/BfZy95Vmtj2hEJyfLPNpwJMJv+0/Bf4FuMLMTmJyDdBxrrsfD2wF3JmaflcybZ1CvYAZZnY18A/geHe/AFgKbGBm8939auBAYOsuz309cGGP5e4MbJUk9yTbCOBLwNvc/XIzO7EzcxKWfrXHsnZ393tT884C9gM+M2CbMLOLCEcnLiSk9GnPBv7i7rf0ev6oqHif7DZ3vwHAzH4D/Mjd3cxuALYBZgJnJn9wDmwA4O5rzOw1wK+Bz7v7T5PplxE+oIW4+xHAEWZ2OPBW4EOEvcnHuvsKM9sXuIBwqCzrcODPhL3HBcD7gH5dbERERCSfqeqFxwD7m1mnS8UMYC6hMD/Zwnloq4EnpJb5C3e/K1nmdclyrnD3d07Rlm4JfVlH6+e6+x/NbFvgUjO7wd1vNbOXAyeZ2YbAxYTifqJBIbx8PWGHppvfA9ua2eeA7wEXm9lMYJa7X57M87+ELkG4+2/JUU8lfe+/AnzW3X+feewQwo7Sc3s9393/n5nNAM4Gnkc4qtDxb8myx07dZiZLJ9NrUvfXEHZ0jgEuS/YU9yP8MXZsD6wAHt2ZYGZ7pE5ySP/72YDtOoeJQ1L3ufuK5Pb3CXu/c7JPcPc/JYfQHiTsyXbr5yYiJTOzo83MzWybcbclDzN7TdLe3XPO72Z2RrWtEoneVPWCAS9N9Ref6+43A+8E/gLsRCgkp/dY5upkOZjZST1qifcn897F5OT7MYSdhKG5+x+T/38PLCQcHcDdr3T3Z7v704EfA2vTaDN7KnA6cIC7/1+P5S4jbIOFhO4xpxO2WdedDjPbocc2uC6V2kMIK29x909nnr8noWfE/lP1QnD3lcC3Cd2ROs9fH3gJvdP/kVLyPpiZwB+S26/pTEz2Fj8DPIewR32gu583TPJuZtunDs3sDyxKpm9JOGzjZvZ0wg7YOn8cZvYod/9T0kfrRUCuM7ZFRMqS7BDsDnw6fVhbpAUuAt5mZm9Lfq+f5u7XEuqIu5Ij9q8G1ptqQTmS928DbzWzc4F/BpZ369s+KDPblNDN98EkJPwXkpM9zeyR7v7XJHl/H3BcMn0u4YTOf0/1ie+27DnAKnf/hpndCpzh7vdaGE1vt6Qv+Ss78+dJ3s3sWML2fUNm+tOAzwN7u/tfezz3EcDGSd20PrAv8JPULHsCizpHRsZNyftgPg58zMx+yuQ/uJOAU5IP6uuB483skXkWaGZfAa4EdjCzu8zs9clDxycncvwaeAHw9mT6gcCNSZ/3zxLO9vZkWd+3ZEhJ4Ozk8N0NwBzg2ILrLCKDORZ4GGEEhrbbndDdb9YU84k0zTGErrW/tjDc4THJ9FOAV5vZVYQuM38v4bW+T+iGshg4DfjPzgNJ95vO7Y+b2V3ARkm9cXQyfddk+kHA55NuQABPBK5O6o3LCH3eb0oee4+Z3UzoLvwdd780mX4U4YTWU5JU/Ooebd6KcK7gdYRR9w5Ppr8W+G8LJ6w+kHcDmNljCMn6k4BrktfuFPEnAo8Avm7rDnnZ2T4PB76d1FzXE85NODX1Ei+nS5cZM7sd+BTwmmSbPilvm4dhSd0nIiItlJyv8yVgD3dfmGN+B85099fkmPdoQvH+OHe/fZh2iohIoORdRFqrX/90CxfjWJiZ9kIzu9zMlprZA2a2xMy+aWZPSM2zzjJT03Yws48mCc2DFoZE27fLa29kZp8ysz+Z2f1mdpWZPd/C0Gp9ExcLQ7D9wcyu6fH4fyRteVHmoWlm9m4zuzVp2++Sw/q9XmfPpF33m9mfzewzZvbw1ONnEAp3gNuS1/RO2iciIsWoz7uISA5m9lxC39IbgI8B9xJOUN+TMGxcz/6dKWcCDwGfIJyo9g7gAjN7QiaZ/jqhz+UFwCXA44DzCRdA6cvdV5vZ2YTD2ju6e/Z8l1cRhnr7Xmb6RwndfT5POHHuzcAZZra4M4JWyi6ELnynAWcRhrE7DNjRzPZy9zXJcjYBXkw4SW9p8txfT7UOIiLSm4p3EZF8DiAcrXxB5qSnY3rM381SYL/UeSqXAb8A/oOkz2eSxO8LnO7ub+w80cwuZd2Cu5czgfcQCvX3ppbxeOBZwOfc/aHMczYEdk0ucIKZnUfoR/tWwpjTaU8BXpyM+Qyhf+tnCAX8ywjjT1+Z9B99MXCBus2IiJRD3WYGYGb728TQTGUu9wdmdq+ZfbfH458zsxWp+xua2VfNbLGZ/bzbIf9kvrcnJ73+xszeUXa7RVpmefL/S5PRCIr4jKdONHL3XwJ/Y/K1GvZL/v9U+onJ0LA353mR5MJtvwJeaeFS5h2vSv4/s8vTTukU7sky/kA4mtDtOhK/TRXuHccn/784TxtFpHpmdqKFK7b+2szOt2RYRTPby8x+ZeHqpr8ys+f1eP4xyXOvM7OLLRkUw8zeYxNDNd5oZqvNbLaZbW1ml5nZzUnt8fZuy5XhqHgfgLt/O7maWdlOBP692wNmNp91R2p4PbDM3bcjjHRzQpfn7Qi8kTC++07Av1q4uJSIFHMycC1htIh7LIzudJiZbT7AMn7fZdo9hNEZOh5HGCt6cZd5fzvAa53FRLeejkOA37j7r3K27f8ybetYZyciGZruXmDbAdoo0moWVFmL/RDY0d2fStgZ74zq0jkK+BTg1YQLInVzors/1d13Br5LGE0Gdz+xM4Z9sszL3f0ewsWa/svdnwg8A3jLqEZgaRMV7wkz2ybZOz092Ys8Ozkh66dmdouZPd3CxUxOTuY/w8w+a2Y/M7Pfm9mBRV/b3X9ESN+ybVqPUNi/N/PQAUwkZ+cBzzez7BXWnghc5e73u/s/gMtRIiaS1e/kz0npenKxkV0J/bs/B2xM2Hn+nZk9M+frre4x3brcHnYosHMI/etfBWBmzyYU1mcN0baOXm3rNq+IpCT1xs1mdgrhqumrzeyEJAG/JKk3Fia1xf7Jc55sZr9Iku5f5w3j3P3ipAYAuIpwASfc/drOBZiA3wAzLIzZnn3+fam7D6f73/7aK48mF4i8Jrn9N8KO/lZ52ir5qXifbDvCxZaeCswDXkG4tO+7gQ90mf9RyeP/SnLI2Mw2tt5XARt07/OtwLe7XGxhK+BOgOSPcjnrpmM3As8xs83MbCNCH9qtEZG0e5L/Z6cnWrg89qOyM7v7andf6O5HuPuzCVcbfATwwRLbdBvhu7nbj/MOeRfi7ksJ4z+/2MIFSF5FSPS/XEIb1/kuM7NHES6Qkk7wNRaxSHc7AGe5+9OS+wvd/Z8IQd6xwF6EwO0jyeNvInS725lwZda7AMzsJz3qjT1Z1+uAC7tMfylwba8rj5rZcWZ2J+GiSUdlHtsI2Bv4RpfnbUP4jvx5j20gBemE1cluc/cbACxcpOBHyZXRbgC26TL/BcmoCjeZ2Rawdk+z0FVV05J+ZQcRLnKyzsNdpk36kXT3m83sBMIhsxWEiw78o8vzRNqsM0LMnoQErOOdZMINM5uTFMRpiwgXEplNeb5DOIH1ncChqdffl3BEbRBnEo7UHUL4PvlhKm0bxg5m9qJMv/f3Jf+np3XO1ZkN3F7C64o0xR3uflVyexXwg+T2DcCD7v5Qpva4EjjCwsWIvtm5AnsSIkzJzI4g1ABnZ6Y/mdD19gW9nuvuRySvfTghVPxQ6uH9gJ8mXWbSy30EoaB/Rya9lxKoeJ8svde5JnV/Dd23VXp+g5C8M/mSummvIBxq/3xy/yh3/3aPeZ9GOBKwOOkRs5GFIdu2I+xxbw3clZw4N5OJBHEtd/8C8IWkXR9NniciEy4hFOAfMbPNCKn3boS+mtlC/bTkh/NiwtVTHwYcTPib7tUVpYjvEy6t/kYLlxDvDBV5KGGYxacOsKzvEfqtn0AYtrHbiapF3AB82cxOA24hdCU6kNA976up+TrFyQkWhq9cCdzYZfhKkbZJX1n1odSJ7GtrD3df0zk53t3PMbOfAy8ELjKzN7j7pWb2E8J3UNa73f0SAAvXa/hX4PnpE+aT77PzgVe5+6052nwO4TslXbyvc+VRM9uAULif7e7fzLFcGZCK95LlTN6nTObd/XvAlp37ZrYiKdwhjDX9asKe+IHApek/yNRzHunufzWzucBLgLz9ckVaIRkT/QDgs8DbCAnYxcBzWXd4xP8FXkP429scuA+4CTjQ3dc5ZDxEm9zMXgocR+hLug+haH8x4bLnuU88d/dVZvYVQlp2H5NT8WFcA7wraeObkmWfDHwgORrZef2fmtn7knlOI/zmfJjQrU9EcjKzbYHfu/tnk9tPJfz2903ezWxvwlGx57r7/anpswiF+OFdruOQfv72nZQf2J8QdnQem0n4rjwkNc0IoeHN7j5pxCwpj3Wp+Vop6Zv1XXffMbl/RnL/vM5jhAurzHf3t6YfT+Zf4e6PKPjaPyH0sX8EISV7vbtflJln7fKT/rj/S0jn7wFe7u6/T7ranO7u+6aWuxnhpLV3JSfGikhNJYfRN3D3eeNui4gU06XeSP++Hw2scPdPpB9LuqwcQvg9/zPwimxXlR6vtZhwDYf/SyZd5e5vMrMPEkaJuSU1+wuSwO904FR3v9rMvkHon7+GcMTxTckwspjZa4C93f3lqdfbjdD74IbkORB26r8/yDaS/lS8i4hExswe5u4PZKa9kBAinOzubxtPy0REZNxUvIuIRMbMPkY4snYZYTSpnQkjRdwH7OzuOn9FRKSlVLyLiEQmGVnm/YQhGTsnpF8KHOnu3S7eJCIiLaHiXURERESkJjTaDDBnzhzfZpu5426GiIiISCP96lfXLnX3zQd5zlwzX1nw9e6Gi9x974JPj5qKd2Cbbeby86uvGHczRERERBppfXv4HYM+ZyVw0HrFXu+U1cwp9sz4TZt6FhERERERiYGKdxERERGRmlC3mYJWsdG4m8B07p96JmmkcX7+qvzcFVmvKa9SUqLZJS9vXH/DMXx/wWjXf9h1ruJzVvbnKa1O23aUytguZa9v1d9hVX7OZDyUvIuIiIiI1ISS95xGmSwoUY/PSv4wltedwVZjed20GD+Po0zb+72uEq3Bxfh56qXqz1l6+XX+LKXf01hS+Ng+Z+P6zsq+dp0/ZzJBybuIiIiISE0oec8ptr14mWxcyXheZbQvhhQ+Fun0qM593rMp5ai+Z/q9ThnJaYzfl0XS4X7vdxmfuyamoFV/tvK+VpVi/CwVfW2pJyXvIiIiIiI1oeS9BWJPpasWw/qv5PpSl9cvhS873YqlD2svTUqVimzrstPHGFPzshVZx+x7E/vnLu9naZTvdxM/W3nXqd/7EftnaZymT4dttiz45CWlNiUqSt5FRERERGpCyXsJYkh2q1bndSw79S7WhoWVLj+dxMc4jnERdUjpYthO/dow7DYc5O++6edkNPXzmH5O1etY5HekKZ+rOnx+eonhe04mU/IuIiIiIlITSt5bKLYUPYZkfCpVJucPclXPx2awe9fpRd/DGFOsOidSw45vXXSUibx9ZMeVqsb4OUuL7Tswq+j2G/bzWOVRnKL6vVexf87yin00sux7ryR+/JS8i4iIiIjUhJL3EqT3eGNPdCCO9saetpedtPdL14sYJGXJO++40pRRpsNFlLFdyh7HucgVE4ts56Ykm1nj+g6senvGnoiWsd2b+JmMfZ1i/1y1kZJ3EREREZGaUPJesuwedOxJfL89/irbPoOdhl5Glel9tq/5sEn8hjxj0v28SXx6O6Xfq0GSmipTkyIJcD+xJDyjvGrrOMTYt3lcbYol9Yzlsz8qsWz3NmjbZ6sNlLyLiIiIiNSEkveKxZAuxNi3sJyz6/Ol92Uk9Okkvoz+8Okkvv9oMxPvwTS2X3v7z0O3oHxVj5zS9DQ8FjGmdOO6WmiM22JcYrxiaxGxHHXSZyuf9afDnLkFn9zgK6yqeG+BcXWN6SfvjkHs3Y7KkO1Sk5beTv0K9rIL21FerrvKojyWgr/K7Vnn975sKojGr87vQZ3bLu2ibjMiIiIiIjWh5D1jlHvesR9erLPYh6LMK33kYXaq20w2be2XlhZJZos8Z1yJbSzpelqd0+sytmed119EJHZK3kVEREREakLJO7CGaWPp6zaqi9PE3m+8jPZVnbSP66JN9/L9tbdnse/a21v2OWeg32e5Df3Lq6REWUSkOcxsO+A9wDOAHYGfuPvuqccfBbwLeAHweGAZcClwuLv/MbOsrYCTgb2AlcC5wHvdvfQCT8W7iIiIiLTRk4F9gauA6V0e/yfgxcDpwM+BLYCjgZ+Z2Y7uvgLAzNYHLgJWAQcDs4BPJf8fUnajVbw31LjS9qpfd1R92ceVtPeT3rbZ7dz/Ak4TSXw6OW5DUl5UXRP2cQ67Wddt1lQa9lAkl++4+7cAzOw8YE7m8SuAee7+j84EM7sG+C3wUuDMZPJBwBOB7dz9tmS+h4BzzezD7n5LmY1Wn3cRERERaR13XzPF4/emC/dk2u+A+4FHpibvA/yyU7gnLiAk8XuX1Ny1lLw3VIxju9dJ+qJMMHwSnx3LfdgkPvv+pu9nUzCl7c0zbMpd9uhEMh6xjFjWqx1K5KWJzOyphEPaN6Umz8vcx91XmdmtyWOlUvEuIiIiItGZNh1mFL3CKswxs6tT9xe4+4Jh2mNm04DPALcAF6ce2hS4t8tTliWPlaoWxXtyIsC7gdcDc4G7ga+7+ztT8xhwOPBmQp+lXwKHuft1Uy1/GmuGSjCKpgvjSk2qvLppetlVJPwz2Gng59R5zPf0CDP90vVVfZZxT4/bWUsKXEp6bvEv1aj12k7q1y0xyfvbE0tCLzJiS919fsnL/BjwTOC57v5Q5jHvMr/1mD6UWhTvwJeA5wMfBhYBWwNPyszzfuBIwpA/iwhD+1ySnA3c78ryIiIiIiI9mdl/EmrMf3P3n2ceXkYYWSZrFt0T+aFEX7yb2d7Ay4Gd3P2mHvPMIBTvH3P3k5NpVwK3A28FPlhlG5uaagybomcT/nH1te+X1udN5dN94MsYiSbdB75f//f0NvxzwSM8vVLkIkl72w1yZdu6auI6yXj0+21Uf3ipEzN7KfA5wrjtX+0yyyIyfdvNbDqwLXBq2e2pw2gzrwMu7VW4J54FbAJ8rTPB3f8OfIdwBrCIiIiIyEDMbHfgbOBkd/9Ej9kuBHY1s8empu0PbAj8oOw2RZ+8A/8MfNvMTgZeRWjzD4C3pq5uNQ9YTTiBIO1mwmD5UsCwSXlTR7UpO4Xvp0g6lXe0kNj7qysBlmGMa5zzph6JLZu20wQdhRgfM9sI1p5cthWwiZkdmNz/PvBYwpCPi4Cvmll66Li73f3W5PZ5wBHAN83sSGAmcBJwTtljvEM9ivctgdcA1xO6z2wMfBw438ye4e5OOJN3hbuvzjx3GbCRmU1390nn9JnZocChAHPnbl3tGoiIiIhIbB4JfD0zrXP/cYQAeSawE/DTzHxnEupT3P2hpJv3yYReIA8C5xL6yJeuDsW7Jf8OcPf/AzCzPwGXA88DfpTM1+ss366PJcMFLQCYP3+X0s8Erosy0vE6J+zp/vBFRqUpYzz47Bjww8om1rGN261EvT6qSEeHTRljTGxjbJPUR9mfHyX5+bn77UzUit2ckfzLs6y7gBcN3agc6tDnfRlwQ6dwT1xBGB3vSal5Njaz9TLPnQXc32U4HxERERGR2qlD8n4zocN/lgGdy9ouAtYDtgN+m5pnXvJYraWT7bxjtPd6vlSr7P7wa1Knccxm+7W3B0nTeyXdVSfyZSfssSf2sR3hKKrqFFkptUi1NMpP89WheP8u8GEzm+PuS5NpzwE2gLX9HH4G3AccBBwLa09C2I+ka4yIiIiI1MeQV1htrDoU7wuAw4DvmNlHCSesngBc4u5XALj7SjM7HjjSzJYxcZGmaYRxOaOXNx3Pm8KPMm2v8oqtozRs//d1l7d7z8eGTeXL6NceY5IdY5vyKqPt6fexztuiTsr4XipyRFRkHPIe+VJCH7foi3d3v8/Mngd8lnDm7irgW8A7M7MeTyjWDwc2A64G9nL3v4ywuSIiIiIilYm+eAdw98VMjMPZax4Hjkv+DWQN09buZY6rP2Ybkpuy17HKJL/fVVknt6FYQt8vlZ8838Q2S491Wuf+1UqUe4th22jM8u7a8B0t0qGEPm51GG1GRERERESoSfI+Snn3ImNPpNIpUdkJdRkJ1Lj6mZa/LfIl9P30S+/Tn8d02t4veR9lehtDUtxPGUcoYl/Htin7u0OJusSkiiS7ynolvWyl8KOj4r2gOhX5Mf44jasLzSi3Rf429d4ByFuw93oODF98Fn3+uA67lt2lqIyTSIt8D8TwQ5htQ5XfZ7Gf0N52ZX8eY/htHKVx/j33eu22vQdNom4zIiIiIiI1oeS9Yjr5azRiHK6yjJT/jyW0o4giCXNTLwxS9rYoexl13rbVnnQ+uqNs/d4Dfe9314ZtFvvf5iiPrEm5lLyLiIiIiNSEkvcRGuVeeJV70FWvR5Vtj7H/fz/p1Lfqk1SHXUYsKdMoTzCNIamKoQ1F1e3vsZc6vwd5lXEkqA3bKZbvwSKK9BSoen1tA9iwGV8TpVLyLiIiIiJSE0reG6pte/9NVSRF1tCGE9r++REZtbr9zdX5t7LOR9hlOEreRURERERqQsk7MI01pe3Bam+1XG3fnunP5ewWbIumjCPelL7c0ltTRsppqhh+O3QUQqqi5F1EREREpCaUvJdAe6vNV/ZVO/PaMnU77zjqdRtvfdh0qg1X5hzXOioBnlD1e6BtPbgYv8/qRNuvvpS8i4iIiIjUhJL3gmLfYx1XUiyD6/dezeqR9mVTurzpdd36YKbVKWHPtrVXqlqndWo7JeNxiP23Ny1vW0f5vVyn7Se9KXkXEREREakJJe/AGqY1Ym9UaXt5YtyW/VLadCoY46gqbU+Y67r+eY8giFSlCb/NUyk7oW/SNps2HWbMHXcr4qPkXURERESkJpS8SyPEmJSnLVlS7HmPnpsvsS0j2S2yjLzPUWLbXdH3bVzbM91evacyDuM6byf2NDv29km5lLyLiIiIiNSEkvcGmZ26HXsSnVfs61E0Uc/rXj6z9vYMdq/2xXpYyfWFnjeDnVLLqGef71jFcMQj9ve0jUcGer0ndd4WsYyQVWU7ykjN0+3rt7zY10PyUfIuIiIiIq1jZi83s2vMbIWZ/cHMzjKzR2fmMTP7gJndaWYPmNmPzWzncbUZlLw3Vp1S+NjbB9Un7L0sv+eCtbcfnH1V6pH3j74xAyqa2A8rnfiXocr1KLutWU1MYvOK/cjAKOUdqUrGo+w0PJajFbEzs/2BrwD/DbwHeBRwLPBdM5vv7muSWd8PHJnMswh4F3CJme3o7n8efctVvIuIiIhI+7wCuMbd39qZYGb3Ad8CdgBuNrMZhOL9Y+5+cjLPlcDtwFuBD4660aDiXWStcaXri2/o/djDVk/cnrXbxA7+zNmT51vJwnIbNSZl9OsfV+JfRL+2VpnK1y2VVjpcHY0gJC22AbA8M+3e5H9L/n8WsAnwtc4M7v53M/sOsA9jKt5r1+fdzLZK+ia5mT0iNT26PkkiIiIiEqUvAs82s1eZ2SZm9gRCt5nL3P2mZJ55wGrglsxzb04eG4s6Ju8nAiuAh2emR9cnSSaLsW/7qNL2ful6Pyt7tC+btD/IVd1nHKENecbQy4j9CMIoR/wZVyofoyJHCpQiD65u1x2Q5rPhrrA6x8yuTt1f4O4LOnfc/Xtm9hrgC8CZyeSfAfunnrMpsMLdU8fBAVgGbGRm0919VeEWFlSr4t3Mng3sDXyUUMR3pkfZJ6lqMZGZqw0AACAASURBVBbDMbYpxu4wRaQL+RmzexfrK+/Jt686Y/aWwzZpkqp3IMrYORhW3p2Lqov8Op9EOyo6SXN0Yu+Gpfe7OpEPD7nU3ef3etDM9gBOBT4DXAhsARwNnG9me6YKdu/29D6PVa423WbMbD3gc8BHgKWZh7v2SQI6fZJERERERDo+CXzb3d/n7gvd/avAi4DdgQOSeZYBGyc1aNos4H53f2hkrU2pU/L+JmAGYUifV2Ye69cn6eDqm1atGNPstBja15R0vai8afuwz4HyE/u8iiT740rr+yX047rYVl55U/06J/RlJMVKc+tDR2Gkh3mEoSLXcvffmtkDwOOTSYuA9YDtgN9mnrtoFI3sphbJu5ltBhwDvKvHXs6UfZK6LPNQM7vazK5eevfd5TdaRERERGJ1B7BLeoKZPRF4GKHbNYQ+8PcBB6Xm2QjYj9DVZizqkrwfB/zc3b/fZ56B+iQlJy0sANhl/vyx9FnqJYYku5/Y29ckvU7UKZqal6Hs164yyc+b1o8yoS9yUm6MaX0bEvp+dBJtM+R9H/XeNdKpwElm9kcm+rwfRSjcvw/g7ivN7HjgSDNbxsSAKNMIXbnHIvri3cyeDLwOeI6ZzUomd86QmGlmq0n1Scqk72PtkyQiIiIiUfossAp4M6Fr9r3AFcDhyXmTHccTivXDgc2Aq4G93P0vo23uhOiLd2B7wkD6V3Z57C7CED/nEGGfpKIy199R0h2x7Z4y+X7ZfeCXpz71/YbL6jWk5CCGGI5rKOMaHSetX0Jfp1Fu8oplyMt+mpjYayjG+ur13um9qS93d+B/kn9TzXdc8i8KdSjerwD2yEzbG3gfsC/we0K/pU6fpGNhUp+kBYiIiIiINED0xbu7L4XJsZOZbZPc/Im7r0imRdcnqSzpJD6GFD7GIwNzM6nxuEafSSfxZaTw6US9jHQ972v1E1tCX/XoNzH2mx9WLOPV91P2WPZ1TvJjH0e9lzak0hrJRsYh+uJ9ANH1SRIRERGRYqYNd4XVxqpl8e7uZwBnZKZF1yepLWI7MgCTk/g6p/B3XTFxe+YzJ25nv8z6peZlf/FNutJrBF+qMfSZh/Yl9DGOgNPPsEl+nZP7cWl7Kp1e/6au73TuX3s78qutNkotxnkXEREREZGaJu8SL6Xw3WVHpUmL5SqtTZdN6GO/UmzsCX2TUvk82j6ufdmyqXxTk+mOtq2vVEvJu4iIiIhITSh5l8poVJp8+qXyN6Vu5+1fPsp+6LH1fx9Ev77y40rl09IJfewpfFY6lW9iCt+PxrUXkaopeRcRERERqQkl7zIysfeHT4shkQeYE2Ha3kt2xJsY2lRUOpWPLYWHeiXxdRhTXqRq6uMuZVLyLiIiIiJSE0reaybbjzwtljQ7j9jXo1cinxVLQh+j2PrDx5CgS29tG72m7dqQRLdhHWU8VLyLiIiISHSmrb8+M2bPKfjsfBfwqyMV7w3SK82OIckeRL9UPi2G9ao6oU8n1jGk13XT1LS9ziPRtI1GkZnQhiS6Deso46c+7yIiIiIiNaHkvQXqlGQPok7rlTehz7qx3GaULsajAU1N29OamLbXuV9729P1NqfNbV53GR8l7yIiIiIiNaHkXdaKfQSYouqU0GdVmWzHmJrn1dR0vSmJep1T9LzakLa3OVVu87pL/FS8Sy55C+C0GIvhfmJfxzoX22kqvOtDRXh9tb34bPv6S7Op24yIiIiISE0oeZfK1Lm7Sl5Vr+OGIwyPmpqIFxF7it6GRDytqel4Xm1Okdu87iK9KHkXEREREakJJe8ydkX6mkO9Evui63h3+iJNLUjGY0+882pbMt6PUvPmJ8dtWMc2W8VGY3tt4xFD/C5cUGpbYqLkXURERESkJpS8S23FPjpMGdIjzKTThwe5qtDympJsl6HN6bjS8GYmxU1dLynPOFN0KY+SdxERERGRmlDyLq1StwtRzdp5op/7KJPiNqfSg2h7gp2H0uDetG2qo4R59KZz/7ibMBQz2wr4LfBwYGN3X5FMN+Bw4M3AHOCXwGHuft242qrkXURERETa7kRgRZfp7weOBE4A9kvmucTMxjaKRPTJu5kdBPw78E/ATMJe0Sfc/SuZ+d4IvBfYGvgN8F53/1Ge15jGmsr2GLX3Xx8xjkuf7qNet5S3bu2tK6W3g9M2Gw/9HjZD3RP2bszs2cDewEcJRXxn+gxC8f4xdz85mXYlcDvwVuCDI28s9Uje30XYy3knsD9wGXCOmb2tM4OZvRw4FTgL2IdQvH/XzHYcfXNFREREpA7MbD3gc8BHgKWZh58FbAJ8rTPB3f8OfIdQb45F9Mk7sJ+7pzfmpWb2aEJR/7lk2oeBM939GAAzuxx4GmFv6ZBRNjarjD1UpRVxKXuUm/798Hdfe7ttaWHb1jcW2u7NoN8NkdzeBMwA/ht4ZeaxecBq4JbM9JuBg6tvWnfRF++Zwr3jWuAAADPbFngC8PbUc9aY2dfT00TGqehFmmaxb9fpbS+w2r7+MhwVtiL51bibzBwzuzp1f4G7L0jPYGabAccAh7j7Q+Hc1Ek2BVa4++rM9GXARmY23d1Xld3wqURfvPfwLOCm5Pa85P9FmXluBmab2ebufvfIWiYiIiIiQ5vGxkOMfnbBUnefP8VMxwE/d/fv95nHu0yzPo9VrnbFu5k9n5C6vy6ZtGny/72ZWZelHl+neDezQ4FDAbaeO3dtElNkD1MpjlRFCXPz6PtCRGJQ40S9FGb2ZEIt+Rwzm5VM7nxBzzSz1YRacmMzWy+Tvs8C7nf3h0bX4gl1OGF1LTPbBjgH+Ja7n5F5OLv303evyN0XuPt8d58/Z/PNy2ymiIiIiMRte2AD4EpCkb6M0O8d4C7CeZWLgPWA7TLPnce6PT5GpjbJu5nNBi4EljD5JNROwj4LWJ6a3tmLyibyfSkVk5jo8ygiIoNoe6I+gCuAPTLT9gbeB+wL/B64A7gPOAg4FsDMNiKM976AMZmyeDezH+dc1kp3f8GQ7enVho2A7wLTgRcmw/R0dPZ85hE2Mqn796i/u4iIiIikJQOiLExPS3p4APwkdYXV44EjzWwZoeZ8F6HnyucYkzzJ+66EYXT6MeAzwzeny4LN1ge+Tji88S/u/tf04+7+ezP7HWGv6KLkOdOS+xdW0SaRUVmcul10xBoRqTf97beHUvMoHU8o1g8HNgOuBvZy97+Mq0F5ivefufuZU81kZq8ooT3dnEI4fPF2wugxz0g9dq27PwgcDXzZzG4Hfgq8mlDsV9UmkdL0+7JesiTVbWbuCBojkqKiMQ55r+qs96s6KqrbITmf8ozMNCeMSnPcGJrU1ZTFu7s/P8+CquoyA3SW2y3Zfxxwu7t/xcweQeindCThCqv/6u43VtQmEREREZGRi/6EVXffJud8pwGnVdsakf7KTmcW3zBxe24qeV+ypNSXkYy5OsoB5E98q6ZEOZ/0+6VtNkGpuTRNruLdzC6j/0D0njehFxERERGRYvIm71/uMX0r4DDQeHZSfzGmM5dcPnF7u6dM3E4n8qOWbkdT1enIRhuOEhQ5AqDkuZli/J6W6hgPYwY7jbsZ0clVvLv7F9L3zWwzwlm3bwS+Cnyk/KaJlKPOX/aLl3effnvm0hDbzKu+LR3j3HEYRlN3OorsaLS94G9qYV/n9arz97TIqA10hVUz28TMjiGMYLcFsIu7H+rud1XSOhERERERWStvn/eHAe8A/oswoP1u7v6bCtsl0lWb05l+iXc2iR/WKJP8Uan6iEGdkv1+aX0bUvmmUNIu0k55+7zfBqwHfJwwOP0WZrZFegZ3v7TktklL6Eu8t17dZqpWZGegiQX/INI7B3Uq5NuizoVunduepu96kXLkLd5XEkabeXOPxx3YtpQWiYiIiIhIV3lPWN2m4nZIjSlNGc5K/tDn0e1H1o5hjfMk2hjUKW1vUteYpqTSaU1cJxEpT/QXaZLyqdgevf4Fen21oUBXUT4eTSxgm7hOg1iVGlVav0MixeUabcbMdjGzi8zsEjP756obJSIiIiIi68qbvJ8F7A9MB84HnlhZi1pEyUP9jTJRTyfAZY8uA81P0euUoEO9U/Q2JMxtWMcqrSp4bUf9borkL97XI5yUumaA59TGNNboC0H6ir3bSxML77oV22l1K7zbVoi2bX2bpGjR32Z1rm+mMZ0ZbDXuZkQnbyH+amABoZvN66prjsj4xF6gx67OxXZWnYrvphSiTVkPkdgMu8NT5+K/qfKONvMLYK+K2yIiIiIiIn1MWbyb2THufmSO+T7s7h8qp1ki1ahbur7dzInb6TS4im4yTUnOY0/Nm5IwN2U96izGRFTdWpqn33sa42ewDfIk7+8wsy8CNsV8hwEq3iUKdSvSh1Xnwjv2YrufOhWwdWprLFSYDG6U20w7CuOn4T/HI0/x/nBgMVMX7yuHb46I9BJj8aXCezTq1NayqBCQqVT5GdGOgcRsyuLd3XONBS8i5Ut3m5k0vYSkvc6Fd1rdCtu6tTcvFdvSJE36PGtHpHkaN+xjEWuYtvbD3aQ/2DZLDy3Vhi40KsTHr85t1/eeSHPp77t5VLxn5N1D1R9DffQbI7YNhf241LmYTavzeuh7SkTGRYl/dVS8F1TGh1I/rOMXe0K/3Rbdp1edtNe5YM2rKeuo75Hxi/G7owhdDEekHlS8j1GRHQD9UFenzgl9UwrRMjRxW+jvfnix/w3HoA7bSDsY7TKNDfSedzFw8W5mewEvBx7p7vuZ2XxgE3e/tPTWyTo0LNN4jCuh7zWeexML1Km0YZ31Nz21OhSYUp0q338ViVIXAxXvZvY24O3A6cCByeQHgM8Czyq3aTKVbHKvH/7RGFdCny5e76nsVcrThmK7CP2ddqeiXMYt72dQRb6M26DJ+zuA57v77Wb2vmTaImCHcpslUk9lJ/RlX3xJBXV1VJTnoyJd6q7fZ1iFvYzCoMX7xsCdyW1P/t8AWFVaiwoysycBnwOeCdxLODrwYXdfPdaGVUjFQtyq7GqjInx09Hc2HBXrIhKrutaOgxbvPwbeDxyXmnYYcFlpLSrAzDYFLgFuAg4AHg98EpgGfHCMTSudCol6yqYxeQuapozfHiP9LYmItFeda8dBi/e3Ad8xszcCG5vZb4H7gP1Kb9lg3gQ8DHiJu98H/NDMNgGONrOPJ9NEoqFDq+VRER63Oo/iJNKh7+xGqm3tOFDx7u5/MrNdgV2BxxK60PzC3ddU0bgB7ANclNnQ5wInAM8FvjOWVpVExUl7taF7jD7f7TXKgkg7CnFTcSxjUNvaceChIt3dgV8k/2IxD5g0VKW7LzGz+5PHon0DRPrpVbyr4BUZjIpDEcmobe04ZfFuZh/JsyB3P2r45hS2KeFEg6xlyWPrMLNDgUMBto68Y7HGdhfQey8iIjKAOWZ2der+AndfkLo/cO0YizzJ+9ap2zOAlwK/BO4A5gJPB75RftMG5l2mWY/pJG/gAoD583fxqgqjIldRHeXyJG4q2Cfosy8iEq/s71UZv19rmDbMd/9Sd58/xTwD1Y6xmLJ4d/fXdm6b2bnAv7n7N1LTXgIcVE3zclsGzOoyfSbd96pGJu+HV4WJNI0+0yIi7VHDC0dGWztOZdA+7/sAr8xM+xbwpXKaU9giQv+ktcxsa+DhyWMitZQ+ya5ufXZ7fXGrqBcRkQjUtnYctHhfDLwF+Gxq2n8Ct5bWomIuBN5jZhu7+9+SaQcDDwCXj69Z+aULHRU30mQ6GiV1cc+4GxCRNox8Ja1T29px0OL9DcD5ZvZe4A/AVsA/gJeU3bABnUq4WNQ3zewEYFvgaOBTMY/T2Uu2uFER0151Tt6HVcYhV/3t1JcK57i07f3Qzkor1LZ2HHSc92vNbHvgGcCjgT8BV7r7Q1U0boB2LTOz5wMnE4b2uRc4ifAm1J5S+for+sM3o9RWtI9S/ri1rSCU+ij62VTRXx91rh2LjPP+EPCTCtoyFHe/CXjeuNtRNRXycSu7GFnJ9al7Ty956dKhIn90VLBLk6U/3yrk41fX2nGg4r3fmO9jHuddZGRGWXysZGHq3htT0+O/WmQTu/moK8/w8hY0KvKlDlSgyzgMmrxvnbm/JeESsueX0xyROKmQGFzZOxhN2RlQyp9Pr6JIf4syCirKJWaD9nl/bXaame0N/FtpLZK+2v6DPiqxFAgPctXa2+liOFsYN6Ww7afIzkCdt0u/Ir/N3wNK7iUPFd/SZAP3ee/iYuCrJSxHEm3+YR6nOv/Yl5Fy17nQ7SXvdqnbumsM/amVXbzV+fthXFRAy7BWo7+9bgbt875tZtJGwCuAO0trUcPoxzRu+lKYMOwOQN0K4LR+616n9VJaXx0VoiISiyIXaXLAkvv3A9cCryqzUbHQj13z1K1YX3nPnyfupKqHyaPQwAx2GlGLemtq+t+U9L7Iybb6DhSRovT9UZ1Bi/f3ufuJ2Ylm9i7gmnKaNHprmKYPWYPVrWAvIlvM9xJDkd9Pnfu117ntvRQdXUffpyL11m9Yag1ZPX7TBpz/yB7TPzhsQ0REREREpL9cybuZdQawX9/M9mCi2wyEy8n+reyGiQisXJK6M6nbzMJJ881g93zLy5nQ9xNbel+0u04MqXdTuuRklTEefprSPRm3sj/TddLmdY9V3m4zX0j+3xD4Ymq6A38G3lZmo0RkXbFcmKmMHYC0ce0M1KmPfp3aWoUYiwftUIxfjJ8LkVHIVby7++MAzOwsd2/kyanSXP1GiWhKf/hsEp9H3rS+anXeGahTP3ddNKtcdS4ch93xqPO6izTBoBdpUuEujZIu7GMs5NPdZlbuPFHkpi/eBLAhzxh82QUK/n6atDNQ5Q5AU1L0qo8ExbCOTaXiWwaho0zxmbJ4N7PnuPuPk9vP6zWfu19aZsNERERERGSyPMn7KcCOye0v9JjHCSeuitSWLrs+nCYl+UXS+9i76/QTY8o9qnM8Ylx3karULUX/B/rN7WbK4t3dd0zdPaXHOO//VWqrRmwaa7oeRkx/yLOP1+0PQMoTYx/6bDeaPIp0tRmlMnYGRrkDEHt3nX7asDPQS5vXve70OyxtNehFmo4E1inegSOATw7fnLj06xdY5z6D+sKrTq/CPsbkoEjBD/EX/Wl5dwDq3F8/tqE7oTn9+osYZN3ruo5V0O+SSH4a511EREREpCY0znsLVXnUQOlJd2V0tSm7T3lRvRL7OiXyWUW3bQyJfd60PsaEvp+mXsAqrQ3r2It+K0SK0zjvUqoiOwZt/xLPe6KsxCdd9MdQyPfTr8ivW2Gf1oYCOL2OdV6Ptn/Xi5Rl2iAzq3AXERERERmfQU9Yxcy2AJ4OzCHV993dv9jzSSJ9FO3G04YUZ8bccbdgXXXuHlO22NN2mdCU9LrO65H+rm/D97dIVQYq3s3sRcCXgVuAJwO/IYwBfwWT+8KLVC5v0d+GHwkV1INrc+Fd524yMiHbZahOxbyGXxYpbqBuM8CxwGvd/WnA35P/DwV+VXrLRERERERkkkG7zcx1969npp1JGHHm3eU0SaRcdT6JNt1tps5JcZ3bHjul6N3VKYUWHUmV7latgiVLxt2K+AyavP816fMOcLuZPRN4PLBeuc0SEREREZGsQZP304DdgG8AJwGXAWuo6OqqZrYJ8F/APsAOwAPAlcD73P13mXlnAp8GXkTYKfkucJi7/18VbZNmK2Ms/DISohmzt5y4rYR17PQejI6Sc+mmzkdSRcoyUPHu7iekbp9lZguBWcCBJberYy7wRsJFoo4ANgIOB35uZk919ztT836VUOC/gbBDcQJwAfDsitom0lcZh4Fn8v6u01VEyiiogI6b3p98YgljRMoyaLeZSdx9CfBbQmFdhduAx7v7ke7+Q3f/FrAvsAHwus5MSfed/we82t2/4e7nA4cAu5nZnhW1TURERERaxszeYWZuZud1eWwrMzvfzFaY2VIzO9nMSt37G3ic9x5s6lkG5+5/7zLtHjO7A3hkavI+wF/c/cep+X5hZrclj11SRftEytAvFUon7ErZRNpBf+vxKSO9L5uOBoyHmT0SOAq4u8tj6wMXAauAgwm9Uz6V/H9IWW0oq3j3kpYzJTPbHNgOOCU1eR6wqMvsNyePiYiIiIgM62PA94Ctuzx2EPBEYDt3vw3AzB4CzjWzD7v7LWU0IFfxbmbP6/Pw9DIaMoBPAiuAc1PTNgXu7TLvMmDbbgsxs0MJY9Qzd2637S8yfrN4eq75shdrSVOK1179PhfD0udKZLyGPRqg5H5wZrYr8DJCMHx2l1n2AX7ZKdwTFxCS+L0JFzkdWt7k/QtTPJ57FM5kVJhHTTWfu6+TpJvZmwmHHV7aZRSZbum/9ZiOuy8AFgDMn7/LyI4ciIiIiEi9mJkBJwMfd/c/hLvrmAfclJ7g7qvM7FZK7AmSq3h398eV9YKEQwqn5Zhv0lYxs/2BzxGGiTw/M+8yYPMuy5hF90RepBbSyUg6ZckmJv1S0LqmKzH2Ma0bpeMi0ksLLow1x8yuTt1fkAS3Rb0W2BL4RJ95+vUE2XSI156krD7vubn76cDpgzzHzJ5F6CZzqruf2GWWRXQfEnIe4XBFX2uYFvWHU0VMe92Tuj079Rm9JzPf7D6f3+y8MZuduh3z36TET9+bIuXo97dU9ff0gw/A4hsKP32pu8/v9eAgPUGSeT9KuH7QA1M9pdvL9ZheyMiL90GZ2ZMJF1z6AXBYj9kuBI40s93c/YrkefMJ/d0vHElDRURERKQuBukJ8gHgTuBiM5uVTF8f2CC5/zd3X01I2Gd1WUapPUGiLt6T4Xh+QDhB9bPA01N9jO5z95sA3P1KM7sIOMvM3s3ERZqucPfaDxPZq+uENN91qbNJdp47cTubptc1Xc+q03rUWb/3oCmadORG3/sSq37dOWM3YE+QHYD5hOI8axmh98cVhJ4gk/q2m9l0Qph8auHGZkRdvANPAh6T3L4s89jlwO6p+y8HTgK+SLj41HfpndSLiIiIiOTxQeDTmWmfBpYDHwI6nXsuBF5hZo919zuSafsDGxLC6FJEXby7+0JyXgDK3e8lnEzw2irbJJJW5VB8AItv2H7t7bmp5H1J7vGdypduRxFK1wdXdlLexvegzkcbYks0dSRA2sbdb8xOM7N7Cf3qF6YmnwccAXzTzI4EZhKC5XPKGuMdQkItIiIiIiJDcPeHCOO53wl8jTC05DdIritUlqiT9zbp129MKcfoVZ2o53V76moHS54ycXuIs+9z2e4pvR8bZ+o/bsMedSiqzkl5LIn3uLZhLOtfpnEeCdDvYbyy701sR4zK5u6795h+F/CiKl9bybuIiIiISE0oeQemsSaqvfmY2tJ0sSTsvfzgdxO3t0ml7bevc/3hfLbJeX23qpP9GPQ7utBLnY86tP2owbgS8LLXv4lJ/iCanuY26fe/SesSGxXv0nixF+j9LF4+cTtdsC/+S+/nbLdF78eKFv0xyLvjkVeVOyhFdgyqNsodj3HtKPQzqp2IqovrGHaG2r4DUSV1SZI8VLyLiIiISHQeXFnv0KkqKt6ltuqcqOdt+xbLJ1KYxX/ZKtdz+qXyRfRL8kdpXF/gRRL/cZ5QHIOyU/4Yk/xeqkjGY0u6R5n+x7buTTZs6q/kfnR0wqqIiIiISE0oeZexq3OC3k8Z6/Xo+/6YujeRvKf7wk9lu5nDtaHsJD+vOif+ZffPz4rhhOJRpv/qrz+cOqfXZaT8dV5/kW6UvIuIiIiI1ISSdylVU1P0tFGu49zlE5Hjb5bvWmgZg6T03Qyb3BdVNPGPIbEvo39+1en9sPKm/7H3z88qkvLHmNantT29bvv6S/MoeRcRERERqQkl75JLUxP1pq5XmYZN7mG06X06sY8hhS8qb3rflIS+n9jT+3RaH3sKX1Te9LqpCXWv9W/q+krclLyLiIiIiNSEkveWa2ryHPt6reT6nHM+Zu2tdAI+rn7oRfVL76tcl6JXoq2TdEIfewpfVL/0PrZUvujIOE1J7Psl9E1Mqdu2vhIHFe8iIiIiEp0HHxrfcMUxU/HeULEnz0XFvl75E/V80qPN9NOUVH5cfeP7qVNCn+0n39Qkvon6JfZNTOXbkEq3bX1ldNTnXURERESkJpS811jsKXRedViPshP1KmX7l8eexMfevqaMXtMUsfVxH4UmjmbTtr7i2fVt4jrK6Ch5FxERERGpCSXvNVOHlDqPGNcjxnT9H8tvW3t7i+WpbTZzq9zLqHN/+Nhk+8kriS9PGxN1aS/1h5dhKHkXEREREakJJe9SGaXr+axkYc/Hbv/b6lJfa1wpfFMS/9iT9jqMLqOEvbum9GWXdpnO/eNuQispeRcRERERqQkl71Iqpe292rCw0PO22XibtbevKaEdStu7U6JenbYn7UrU81G/b5H8alW8m9k7gJOAb7j7gZnHtgJOBvYCVgLnAu9199of04mxIE6LsX11LtjT1p/5uKGXMaoiOvZiPZYCvc6FeFobinIV3uVSgT5B2yKfB1evO/yx1Kh4N7NHAkcBd3d5bH3gImAVcDAwC/hU8v8hI2ymiIiIiEhlalO8Ax8Dvgds3eWxg4AnAtu5+20AZvYQcK6ZfdjdbxldM4cXY5KdFmP7YkjaoZy0Pe2Kubutvf2X1PCQG2fma3N3GCXq1alzuq7UPB8lwNXRtpWq1OKEVTPbFXgZ8P4es+wD/LJTuCcuICTxe1fcPBERERGRkYg+eTczI/Rl/7i7/yHcXcc84Kb0BHdfZWa3Jo9JA8WQtpedtGddu/Wua2+nE+9s+l1Gn0Al6t01JVFXil4fSmzjpvenncxsE+AjwIuARwJ3AJ8HPuPunpqv8nMwoy/egdcCWwKf6DPPpsC9XaYvSx5bh5kdChwKMHdut544IiIiIiIAnAE8B/gAsBjYg3B+pREGUxnZOZgjL97NbCbwqKnmc/dFybwfBQ5z9wemekq3l+sxHXdfACwAMKtwnwAAHWJJREFUmD9/l67zSNxmsNPa2zGk8FXYYnnq/IJUn/dYUul+6tDGjljS9Tqn4720LTXvR4ltfem9azcz2wg4AHhHUj8CXGpmTwZeTlK8M6JzMMeRvB8EnJZjPiPs3dwJXGxms5Lp6wMbJPf/5u6rCQn7rC7LmEX3RF5EREREJI/1CeeJZjup3gs8JnV/qnMw61m8u/vpwOk5Z98BmE8ozrOWAc8GrgAWkenbbmbTgW2BUws3dkxmsFXPx2IY6SXG9qVT+NCO0STxM9g987oLS13+o+/749rbf+uz3ctWp9Q8K5YUvZcmpuvQ7oRdqaxIs7n7fWb2NeC9ZvZr4FbguYTBVN6SmnUk52DG3uf9g8CnM9M+Tdjz+RBwQzLtQuAVZvZYd78jmbY/sCHwg1E0VERERESiMcfMrk7dX5Dq8lLEq4CzgWuT+w4c7u5npuYZ+BzMIqIu3t39xuw0M7sXWOruC1OTzwOOAL5pZkcCMwn9j86p2xjvU+mXeqeNLwGPI5XPJvETbag2kc8m8ROvu3DoZadHg4k9XR5Ek9alI5Z0vc1peFFK0aWbe1K39RmZsIqN1t6eTvkXtF853BVWl7r7/F4PDnIOZnLzJOCfCQOp/B7YDTjazJa6+xfST+n2cj2mFxJ18Z6Xuz9kZnsThub5GvAgYWie94y1YSIiIiISo9znYJrZU4E3Ay9w9x8m039sZhsDnzCzL7n7GkZ0Dmbtind3373H9LsIY28KkxPwGPrJQxxt6pXIQ7WpfNG+8elx3vtdbayJ6XVesaTco6REfWpKR8djlNv9nqlnKUyfn+Yb8BzMzq/sdZnp1xIK882AuxnROZi1uMKqiIiIiMiYdM6n3CUz/Z+AvwNLk/sXArua2WNT85R+DmbtkneRqoxy3PhefeOz9k4lrOkRYEaZNi++Yep5umljIl6muqXrvZLKMtJRpaDjF+N7EGObpLGuTv590cyOAm4j9Hl/B5OvsDqSczBVvIuIiIiI9ODuq81sP+BY4Chgc0IafzTwydR8IzkHU8W7jEUM/d/7GVff+Kw9n9t9etWp7JIlE7fLTtDrlig3xbhSyqamo3mPNBRZ/yr7cg+iqe+dSBHu/mfgDTnmq/wcTPV5FxERERGpCSXvDRJjgp1H3rHrs2IcsaaXomn9817YfXoViVg67VM63gxKTrurcruUsWy9bzKMKsZbl7goeRcRERERqQkl71JbRRP7tFGl90XSeuidwA2SzOXtP1sk7cv7nFj68Ep5lA5LHfVLpdNXC62z7HrUOYnfYPVDPHJ5sd/pevZFyEfJu4iIiIhITSh5B9Ywbco97jrvuUpvZaT3w6o6/S87IS3ytzC7hERL6f34KW3vruhnU9uzPGX8RhdZRh3S+nQbVcs0g5J3EREREZGaUPKeUx36kMWQIscuxhF5+r1vvTKdGD9//eRtb78Uq8oreFZt2IS1DuvYdFW8B+llKoWXUVAK3wwq3gvSH0A91W8Hp12frfTfUt7D0UWLnjoVxCrsxi/7HpT9+VEhP5x+3xf6je5O26y+1G1GRERERKQmlLwD01ijvUzJbZQnKK3hlpG9VloMRyiq/pss4yRaaS+l4/VR5Ds7+/1ThxNTy6TeBXFT8i4iIiIiUhNK3luojBRiXNqWfqQVPdm2SIqefq0YUvgqlPGZjv3zWPbfbezrW6UqvgPbvD1jV8Z7E8vvZl691jnvtqhifTdYvYpH3/fHQs+Nb3iK8ih5FxERERGpCSXvY1Sn1KVOba3CuEYmmZHKDlZyfZ/5dur5WK/EvimJ+iiH/8xuszKGwMwjlgRvlJeWj2H0laq3e6/lx/h9q6M4g1O/camKkncRERERkZpQ8l6xNqQLMYh9zO6i7Xt0qa0oJptsx5DYj+tiW0W3xbhSt7zbqYz3tMgY/f2k0/Yy/r6LpPejHAc7xt8KpcXlifFCjzF+5iQfJe8iIiIiIjWh5L0Ese+9Vp1K5020YkzHY2jTkiW9H9t2bnUpd9EUeZRprgTjOtLQzyi/92IcUz327/1x0XbJp8pR3/QeNJ+SdxERERGRmlDyntMo92T7pcExptzjSq9jSM2z+qXoRUxj+7W3Z/SZr4xkNsZ0N3ZFxsOPYTsXPSpS5fdg269oGSO9B3GJ8f2Ioe9+G9UieTezx5rZV8zsHjO738yuN7O9M/NsZWbnm9kKM1tqZiebWXyfdBERERGRgqJP3s1sa+BK4HrgtcDfgZ2Bh6XmWR+4CFgFHAzMAj6V/H/IVK+xhmlj36PNmyLHmDbnVae2l52gZy2+Id98O8+duL1lKi0dpL96LCOzdDS1L3wMiXod5E3q6jQGet2UPTKQNNs40/Xpq1cxd3mxH+RfltyWmERfvAMnArcCL3T3Ncm0SzLzHAQ8EdjO3W8DMLOHgHPN7MPufsvIWisiIiIiUpGoi3czmwm8BHh1qnDvZh/gl53CPXEBIYnfG4i+eO/Xl33YKw1WnXjXKVHvp+y0PW+63k+R9zv2BDjGcePbYFxjuZed2o3yKq8xqnJ7tmH7tZ36qDdD7H3edwE2ANzMfmpmD5nZXWZ2uJlZar55wKL0E919FSGxnze65oqIiIiIVCfq5B3YMvn/88D/AB8E9gCOBZYDpySPbwrc2+X5y5LH1mFmhwKHAmw9d263WaIx7BjHeVP9MpZftxQ+xrQ9LZ2SpBPrqtP1lVw/9DJmsFPO16rPuPFVb/cY1jGvGBO8uqXysW3DNoz4U7fPSF6xfZakWiNP3s1sppnNm+pfpn0Xuvv73f0ydz8KOBM4PLNo7/ZyPabj7gvcfb67z5+z+ealrJuIiIiINI+ZHWxm3zSzP5mZm9lrusyzp5l91czuSEZHvNHM3mpm63WZ9wAzu8HMVprZTWZ2cN62jCN5Pwg4Lcd8xkSQe1nmsUuB15rZJu5+HyFhn9VlGbPonshLouxUvs4pfF5lp+u3L+r92MoXlpu2l5GoD/taeRP5dZcXd1/+2BXZfnU6EjAVJZODK3ubVXlV0TLoMyJTOBDYBvgu8IYe8xwKbEToKXInsBvwSeBxwH91ZjKz3YBvEHqQHAbsC3zFzJa5+8VTNWTkxbu7nw6cnnP2m3tM7/R375zEuohM33Yzmw5sC5w6aBtFRERERFIOdvc1ZvYIehfv/+nuS1P3FybXHHqnmX3A3R9Mph8J/NjdD0vuX2ZmTwaOAuIr3gfh7reb2W+A5xP6vXc8H7jV3Vck9y8EXmFmj3X3O5Jp+wMbAj8YWYMbZtgUPZvqNzWJr1KR5DP2hLpf+l80lZdqxPJZatIRgDZTsi11NsWoh515lnaZfC3hIumbAHeb2YaE8zcPy8x3LvAlM5vp7sv7vU7so81A2Ds50MxONLO9zOw44N+BD6XmOY+Qvn/TzPY1s38DTgbO0RjvIiIiIjImzwKWuvvdyf3HE0ZSzHaavZlQlz9hqgVGnbwDuPv5ZvYq4Ajg7cAS4C3ufnZqnofMbG9Cwf414EHCHsx7Rt3eQZKFOp/ZnkcsSXsZI8qMsp972p9Tn5HZbL/2dhkJ1ij7v+fVtlReiXI+6SMA2mYi7THMFVaBOWZ2der+AndfUEKzcjOzJwFvAk5ITe6Mgpg9J3NZ5vGeoi/eAdz9y8CXp5jnLuBFo2nRZEULqRgujlF2gV11wV720I69jKtYz7outb7pEU1nZz4v6S5K2eImlq4PvTSxKM9SwdmdtouIVGipu8/v9WByIdBHTbUQdy/0C25mmxJOSv018NFui84+pcf0ddSieBcRERERKdEgox8OxMxmAN8inHu5f3Lh0I5Owp4dJbFzf8pRElW8A9NY04gTaWLoptKGoSKbIpt4lz20YxsoOR6ctpmIxGDA0Q9zS8Z0Pwd4MvAsd/9LZpZbgYcIoyRenpo+jzCK4u+meo06nLAqIiIiIlIHpwB7A/u5+2+zDybDRV5GSP7TDgaunGqkGVDy3ii9Lrg0rgS83wWg0gZpX7rfd1rZfeG3e8rk+8P2gd9m3uT7efvAT+7n3v02lHMCqxL2welEysE1ZZvlPZekzusoIhOSk0+fRBj2EWC+ma0A7nb3y5N5PkC4UNPHgDVm9ozUIm5KLiwKcAxhDPhPAxcQLtK0L6Hon5KKdxERERGR/l7G5GHK35L8uxzYPZn2guT/w5N/aXsACwHc/QozOxA4FngzcBvwijxXVwUV7wCsYVrX0V6q7gc/qhFmqkjAy9SvfXnblE3kq0ziyxiJJpvE95LeNlumPo/Z1G9ln2WMa7SZKhPHGEfQ6dcmpa/tVeSzqs+LSHzc/Wjg6Cnm2X2A5V1ASN0Hpj7vIiIiIiI1oeS9j6ZfRAmaOyJMOomPPYXv1/999gsnbqcTvFEmzzGmgHnbFEtC35R+3v00db16Sa9v2Z8zHcURkX5UvIuIiIhIfFav4h/Lbxt3K6Kj4r0FYk/XY29fVtlXX81Lidvg+m2zKlP5WN6rWNrRdNrOIjJK6vMuIiIiIlITSt5LVrcUuZ/Y16XsvuyjlHec98Wp29ulzsGYxvY9n1P1KElNGd+6yj7L/7+9ew+ypCzvOP797S63lZvcghAuEhCCFKIFBEwiKBQCxUUFXBBNtiqhNFUGjSEQAwgiRkgqSpCQrWiJ8QbIspACXdZs1hU0FIsIwYBLFFkuslBchl2WheWyT/54e9ye5syZnjN9prvP+X2qumbm7bf7vP1Md5/3PP12n7pUHfN+3/fTy77atnuRBuHbuyfStv/JIBqG/cw2cObdzMzMzKwlnHmvQNMz1N3U1famZ82rHtdeNtNeNCZO43zb6uuNzYKVfc5/2czNdGXUpzOb1+8rGU0bX9+UTGlT2tFPVW9jEzKsw/B/a5t+/k+asM/ZWM68m5mZmZm1hDPvVot+Poe9CvlnucPUM/HdvlG1l6x88YpJ2ex6N/nMTV2ZliZm9Mq2qVvMxsuOT+dY+15i2+uVuSr2xyaYziuTZWPWhOO0V1XEc1D2rTYpnjvatt8NImfezczMzMxawpn3IddLFqPqbFQ+C98PVWT2i5n4yeqWue+WlR9P8f9WRTaqCdmUbm1oYlZ+qjHr9/0DU43ZdGY5e41lP/cLZ3k3qCLOTYhn284xdWnC+4GNz513MzMzM2ucdevXseL5FXU3o3Hcea9AE7IJveoli9627d0ml9mv6+k6xasLZa8GHFDyCTPDkCUZhm2cqjZnDpvSdu9nGzTlf1KlbuO3/b+3tvCYdzMzMzOzlnDmfci1LYte1ngZ9sZs7xTH+deZIWrCc6ub+I2oTf9206ZkFevK5k51+/u9z03nNxIP6r6VV3Ybe4lF07/7oZumf/O1lePMu5mZmZlZSzjz3kWbswnDrp8Z9ul8VnFjrhS0VNnsVlOyUcNwfOfPq03c3royovnXbfpVnDZo2n7WlKuF/dzPbPq4895FtwO+ro592ddtwsmqihg1YTuK3KGuX/5Np9ubYr/fnPzmNzVN62BB+X2rn69r1ari8ZBVPxbW+5ZNhYfNmJmZmZm1hDPvwAzWN3KIzFQMyvYMynYMkqb9TwYpq9S02E6nJm77oOxbTYxtU9QVm0HZt6wezrybmZmZmbWEIqLuNtRO0vPAA3W3Y4BsBzxddyMGhGNZLcezWo5ndRzLajme1aoinrtFxPaTWUDSLdlr9+LpiDi6x2UbzZ13QNJPI+LAutsxKBzP6jiW1XI8q+V4VsexrJbjWS3Hs1k8bMbMzMzMrCXceTczMzMzawl33pN/q7sBA8bxrI5jWS3Hs1qOZ3Ucy2o5ntVyPBvEY97NzMzMzFrCmXczMzMzs5Zw593MzMzMrCWGpvMuaZ6k5ZLWSBqRdKukIzvU21nSDVm9pyVdIWl2h3pnSPqlpJck3SXpiOnZkvpJ2lLSZyUtk7RK0hNZzN5SqLe7pOgwXdNhnSdK+nkWz/slzZm+LapP2VhmdbeSdFW2/66S9G1J23aoN5SxHCVpjqQFklZm+9vcDnUOH2ffvKRD3aE91qFcPLN6Pnf2QNLScfbFTQv1SsV3mEnaV9J/SVor6XFJF0maWXe7mk7S3HH2wY/l6kjS30l6VNKLWR/qgDrbPcyGpvMObAZcAbwf+DDpywYWSjpktIKkWcAiYDdgDvAJ4BQKN2pIOhWYB3wDOAa4D7hZ0n7934xG2BU4gxSrk4GPAm8C7pC0S4f6ZwGH5qbz8jMl/RFwPfBDUjy/B1wt6ah+bUCDTCaW1wKHA38OzAUOAm7MVxjyWI46GdgduLlE3dMZu2/+S36mj3WgRDx97pyyHzJ2PzwUWDc6s2x8h5mkNwKLgQBOBC4C/hr4bJ3tapn3MHYfXJCb97fA+cClwPHAGmCxpB2nu5EGRMRQTsBM4BHg8lzZacBrwJtzZR8E1gN75coeAL6W+3sG8HPgW3Vv1zTF7g3AZoWybUgH8wW5st1JJ9LjJljfImBJoez7wI/r3tYGxfLQLJbvypUdnJUd6ViO2d4Z2c/Ns/jM7VDn8GzefhOsa6iP9UnE0+fO3uO7FJg/QZ1S8R3mCfg0MAJsmSs7G1ibL/PUMXZzs2N783HmbwqsAj6TK3sD8BRwcd3tH8ZpmDLvY0TEa8BzwMa54mOAOyPioVzZjcDLwNEAkvYA3gJ8N7eu9cB12fIDLyJeiIgXC2XPAg8DO0xmXZI2Ad5NLp6Za4BDJW01lbY23SRieQzwZETcmqu3DHgomzf0sRyVHY9T5mM9KRlPnzv7a8L4GscAiyJida7sGtJV98PqadLAeCewJWOP3ReAm/CxW4uh6rxnY7ZmSdpW0l8BewFfy1XZB1ieXyYiXgYezOaR+zmmHvALYBtJ21ff8ubLtntP4P4Os6+S9Fo2ZvaLkjbLzfs9YCM6x3MG6c1+qIwTy9ftm5lfsGGfdCwnb0m2b66QdF5hfKyP9fJ87pyao7Jx2mslLZK0f2F+mfgOu04xeoSUeXeMynlQ0quSHpD00Vz5PqQrP78s1M+//9g0mlV3A6bZHODq7PcXgDlZ9nLUG0nZ+KKRbB65n8V6I7n5T029qa3zT6ShHvmbUdeRxhD/AFhNGqpwDqmTeWJWp0w8h02nWHbbN/fI1aFDvWGO5XhWAZcAt5Gyl8eRxsZuTxpPDD7WJ8Pnzt79CPh34FekMe3nArdJeltErMjqlInvsHOMereSNJ59GWlI8WnAPEmzI+JLpPityUYs5I0AsyVtnH2YtGnS2s57NgTgTRPVi4j8J/FFpJv8tiPdqHaNpGMjYml+kU4v16G8+Le6LN94PcZzdNm/IN0EfFJEPJOruxL4eK7qUklPAldKOiAi7smvurjaccobrx+xHF2k08t1KB+YWMLU4tml7t3A3bmixZLWAZ+S9LmIeDpfvdikccpboR/xHF2k08t1KB+oeBZNNr4RcUGu+DZJi0kZ5E9m028X6fRy45QPK8eoBxGxiNQ/GrUwG4Z5nqR/Hq3WYdGBOnbbpLWdd9Kd9l8pUW905yIiRoCfZn/eImkn0h3p78rKRoCtO6xjazZ8oh/Jla0q1IHOn/zbYNLxBJB0AvBl4JyIuKHE8vOBK4F3APcwNp55bY5nP2I5QsoKF423bxbrQDtjCT3GswfzSTe47Q8swcf6ZOI5zOfOoinFNyKekPQT0jlyVJn4DrvxYrQVjlEv5pNuit6dFNstJM0sZN+3BtZGxCs1tG+otXbMe0R8NSI00TTBau5mw5ADSNmOMeO3JG2c1Vmeq0OxXvb3sxHRysu+vcRT0jtJQzvmRcQ/ln2pws8HgVfoHM/1wP/1tkX16VMsX7dvZvLjPAcullDZsT6pl8x++lgvb2jPnUUVxjefzSwT32HXKUa7kJ6K4hj1Lkjxm0m6FytvvHuxrM9a23mfKkkiPX4vf/f+QuAgSbvlyk4ANgFuAYiIX5M6Qafk1jUj+3thn5vdGJLeSnru8y3AmZNY9OTs510AEbGO9IzjUwr15gC3R8QqBlzJWC4Edsye4z663IGkN++F4FhW4CTgVeBe8LE+ST53VkTS7wB/SHaOzEwYX2Mh8F5JW+TK5gAvku4rsMk5ifR9OA8D/026by1/7M4mPe/dx24N2jxspjRJf0z6soYFpGe7bwv8KXAIaecbNZ90s9ACSeeTLrd9CfhOROTvsr4Q+JakFcBPsnXtBXyorxvSEJJ2IL1hrAEuBw5On4UAWB0R92f1LgS2IMVoNWl40t8ACyLi3twqP0caD38Z6fFnx2bTwD8CrWwsI+J2SYuAb0g6i5RJv5T0/PbFuVUObSxHSdoX2Jf0bGKAAyWtAZ6KiB9ldf6VdHPknaQbVo8l3Z9xWeFegwsZ4mMdysUTnzt7kj1V5gukx2U+TPrStk+Tju/LclXLxneYzSMlPxZIupSU2LgQ+GLh8ZFWIOl60s2q95Iy7HOy6czsca4vKX379PmSRkjZ9k+REsBfrqfVQy4a8LD5fk+kMVvzgcdIT0B5jJTpPLRD3d8ldXrWAM+QnpYyu0O9M0hPB1gH/Aw4ou7tnMZ4Hk66lNZpWpqrdyrpHoNVpA7Sr0j3GGzSYZ3vA/43i+dy4NS6t7NJsczqbg1cRRq/uRr4DrCdY/m67b+wxL55JumN6vksTveRbg6c0WF9Q3usl41nVs/nzsnHdmfSl6itzM6Rz5C+IXmfDnVLxXeYJ9KHzCWkbPtKUjJjZt3tavoE/D3pC9TWZrG7C/hIoY5IHyAfy+rcBry97rYP66Tsn2JmZmZmZg03tGPezczMzMzaxp13MzMzM7OWcOfdzMzMzKwl3Hk3MzMzM2sJd97NzMzMzFrCnXczMzMzs5Zw593MzMzMrCXceTczMzMzawl33s3MaiTpPkmHV7Sur0u6uMdlQ9ILkj4/xTYskfSSpB9PZT1mZtaZO+9mZhlJKyQdOZ2vGRFvjYildb1+wdsi4tzRPyT9gaTlku6XdGCufIWkFyWtyU1XAETEe4CP1dB2M7OhMKvuBpiZWWN9Hjgu+/1K4KjcvOMjYvH0N8nMbLg5825mNgFJvy9pqaTnsmEuJxTmr5B0lqR7Ja2SdK2kTbN575B0t6TnJV2Xzbu4sOyRkr4J7ArclGWyz87mh6Q9c/V/OzRG0tsl/Sxb97XApoV27STpeklPSXpI0pmT3PSZ2TQLJ3vMzBrBnXczsy4kbQTcBPwA2AH4S+DbkvYuVP0gcDTwZmB/YK6kjYEbgK8D2wBXA+/v9DoR8RHgEVJGe/OI+IcJ2rUxcCPwzWzd1wEn5ebPyNr9P8DOwBHAJyW9t+y2AxcA3wP+Azh3grpmZjYN3Hk3M+vuEGBz4JKIeDkilgA3A6cV6l0eEY9HxLOkTvMB2bKzsnmvRMQCYFmF7doIuCxb93zgztz8g4DtI+KirN2/Br4CnFr2BSLi1ojYMyL2jojbC7NvzK5EjE5nTHWDzMxsYr4MambW3U7AoxGxPlf2MCmbnfdE7ve12XI7Ab+JiMjNe7TCdhXX/XDu992AnSQ9lyubCdxW0eu/z2PezcymnzPvZmbdPQ7skg1DGbUr8JsSy64EdpakXNkuXepHh7K1wOzc3zt2Wfeuud8fBR6KiK1z0xYRcWyJdpuZWUO5825mNtZGkjYdnYA7gBeAsyVtlD2T/XjgmhLruh14Dfi4pFmSTgQO7lL/SWCPQtk9wIckzZR0NHBYbt2vAmdm6/5AYd3LgNWSzpG0Wbb8fpIOKtFuMzNrKHfezczG+j7wYm76DHACcAzwNOmRiX8SEcsnWlFEvAx8APgz4Dngw6Tx8uvGWeQLwHnZGPKzsrJPkD4sPAecTrpJNb/uucAIMAdYkHvt17LlDgAeytr+VWCriUNQyuhTcUanGypar5mZdaGxwyXNzKyfJN0BzIuIq+puS56kl0gfKi6PiPOnsJ7/JN1MuywijqiqfWZmlrjzbmbWR5IOAx4gZb5PB+YBe0TEylobZmZmreSnzZiZ9dfewHdJj5t8EDjZHXczM+uVM+9mZmZmZi3hG1bNzMzMzFrCnXczMzMzs5Zw593MzMzMrCXceTczMzMzawl33s3MzMzMWsKddzMzMzOzlnDn3czMzMysJf4fv7V6GPp9i2AAAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 864x396 with 2 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"plt.rcParams.update({'font.size': 15})\n",
"fig, ax = plt.subplots(nrows=1, ncols=1, figsize=(12,5.5))\n",
"dummy = np.ma.masked_invalid((psi_ridge1.mean('time')/1.0e6).values)\n",
"xyplot(dummy, grd_ridge.geolon, grd_ridge.geolat, grd_ridge.area_t,axis=ax, clim=(-150, 150),title='using vhbt')"
]
},
{
"cell_type": "code",
"execution_count": 25,
"metadata": {},
"outputs": [],
"source": [
"psi_ridge2 = -ds2.uhbt.cumsum(dim='yh')"
]
},
{
"cell_type": "code",
"execution_count": 26,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAu0AAAF8CAYAAAB7bIUcAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nOzde9zcZJ338c+vQClVaCknBakVi1bwgAoqLioq7gI+gCcWl3VVXOXRVVlwPYCKC+IBjyiy6iK64OMBRQVPIC4CKgoIKohAhSrlKAhSwFJKS/t7/rgyvXOnM3NncieTK8n3/Xr11ZlMJnMlM/fML99cuWLujoiIiIiIxGtG3Q0QEREREZHhVLSLiIiIiERORbuIiIiISORUtIuIiIiIRE5Fu4iIiIhI5FS0i4iIiIhETkV7zczsS2b2FzP7fWb6N8zsiuTfUjO7IvXYk83sYjO72syuMrNZfZY7z8z+18yuT/7fPJm+uZmdaWa/M7NfmdkTq19LERERmS4ze4yZXZr8tn/DzGYOmO8oM1tiZn8ws39ITV+a1A1XmNnlqemDagYzsxOTZf3OzJ6Wes5Hkzrk2mQeS6bPNLOTzew6M1tsZi8vuK4LsrVRzud9L/08MzswaedaM9t1iuduYGa/NbMfpKb13eZm9mgz+0myXS40s0eN2tZRqWiv36nA3tmJ7n6Qu+/i7rsA3wa+A2BmGwJfAd7o7jsDewKr+yz3SOAn7r4j8JPkPsC7gSvc/cnAq4FPl7o2IiIiUpWPACckv+3LgH/NzmBmOwGvBHYm1BefNbMNUrM8P6kv0gXsoJphH2DH5N+hwOeS13g28HfAk4EnArsBz0ue8x7gL+7+OGAn4KfTXem8zOxlwPLM5N8DLwN+lmMR/w5cm5k2aJt/HPhyUk+9H/hw0XbnpaI9kezRLTazU8zs92b2VTPby8x+kexdPSP598tkL+yXZvb45LlvM7MvJbeflDx/dp7XdfefAXcPaZcB/wh8PZn098Dv3P3K5Pl/dfc1fZ56AHBacvs04CXJ7Z0If5C4+2JggZltk6etIiIibZazFnhYcpT8sqQeOCD13J+b2W+Sf89Opu+ZJLHfSpb91V4qPWLbDHgB8K1kUvq3Pe0A4HR3f9DdbwCWAM+YYvGDaoYDCIWpu/slwFwzeyTgwCxgJrAxsBFwR/Kc15EUsO6+1t3v6rMuG5jZqck2vsrMjkimP93MrjSzi4E3T9Hm7DIfDrwN+EB6urtf6+5/yPH8RwEvBk5JTRu2zdfVU8AFhG1VKRXtky0kJM9PBhYBBwN7AG8nJNSLgee6+1OB9wEfSp73KWChmb0U+B/g/7r7CjN7vk10cUn/++UIbXoOcIe7X5/cfxzgZnZu8qXwzgHP28bd/wyQ/L91Mv1Kwh4nZvYM4NFA5Yd0REREGmKqWuA9wPnuvhvwfOBjZvYw4C/Ai9z9acBBwImpZT4VOJxQ6O1ASKkxsxMG1AlHsr4tgHvc/aHk/i3Adn3m2w64OXU/PZ8DPzazX5vZoal5BtUMfZfl7hcTCtU/J//OdfdrzWxuMt9xSY1yxoBgcJdkOU909ycRaieS/w9z993TM5vZ4wdspyvSrwl8AljR5/Xy+BTwTmBtatqwbX4l0Ov681JgUzPbouBr57JhlQtvoBvc/SoAM7uacKjIzewqYAEwBzjNzHYkfPA3grAnaWavBX4H/Le7/yKZfgHhgzkd/8REyg7hPduDcChqBfATM/u1u/+k35P7OB74tIU+8lcBvwUeGv4UERGRzpiqFngUsL+ZvT2ZfxYwH7gNOMnMdgHWEEK2nl+5+y3JMq9IlnORux8xQrv6pfM+4nx/5+63mdnWwP+a2eLkiP9Ir2lmC4EnMBH6/a+ZPRe4Jpn2C3d/m5m9jdCN5F8yy/gTsIOZfQb4IWFHYg4w19173Wn+H6F7DklSPrCeSrb5Qnc/wswWDFmfQc//P4QuPb82sz3TD/WZvbct3054v19L6HpzKxXXU0raJ3swdXtt6v5aQrF8HHCBuz8R2I/wh9qzI6Ef1ba9CdNN2i30X38Z8I3U5FuAn7r7Xe6+AjgbeFqfp9+RHMIi+f8vAO5+n7sfkvSVfzWwFXBDnvaIyHBmdoyZeZEfjbokh6j7/fD3m3fPZP1eW3GzROo0VS1gwMt75525+3x3vxY4gtBF5CnAroSuI/2WuSZZzpRJe3JU/QozOwW4i9A9pRe4Poqwo5B1C7B96v66+dy99/9fgDOZ6DbTt2YYsqyXApe4+3J3Xw6cAzwL+CshUDwzmf8M+tQo7r4s2U4XErrBnJJs177fRTmS9t2Bp5vZUuAi4HFmdmG/ZQ3wd4QdsaXA6cALzOwrDNnm7n6bu78s6X3xnmTavSO85shUtI9mDmFPCuC1vYnJ3uGngecCW5jZKyAk7ak/6vS/Z+d8vb2Axb2988S5wJPNbHbyIXoeYc8263vAa5LbrwG+m7R1rk2cbf564Gfufl/O9oiI5GZmrzWzw+tuh0jJzgXemvR3xsyemkyfA/zZ3dcSkuUNBjx/HXc/YkCdcHzy+D8k91/v7k7okvKK5Onrftszvge80sw2NrPHEELFX1noi79p0uaHEc6R+33qOevVDMn0V1vwLODepPvMTcDzzGxDM9uIUItcm7Tx+4RBMgBeSJ8axcy2BGa4+7eBo4Gnufs9wL1mtkcy2z+nttMfBmynXdz9Hnf/nLtv6+4LCL0RrnP3PcnJ3Y9y90clz38lofvTq4ZtczPb0sx6dfRRwJfyvl5RKtpH81Hgw2b2Cyb/MZ4AfNbdryOcVXx8cuhpSmb2deBi4PFmdouZpc8EfyWTu8b09k4/CVwGXAH8xt1/mCzrFJsYzuh44EVmdj3wouQ+hMNZV5vZYsJhp3/Pt+oiksMHgE2AG+tuSCReS+jHK9ImxxG6x/7OwtCCxyXTPwu8xswuIXSNub+C134X8DYzW0Lob/1FADPb38zeD+DuVwPfJBTLPwLenAxYsQ1wkZldCfwK+KG7/yhZ7qCa4WxCV5YlwBeAf0umfwv4I6Gb7ZXAle7+/VQbjzGz3xF2Xv6jz3psB1yYdBU6lVD0AhwC/JeFE1EfKLSFMszspWZ2CyGN/6GZnZtM39bMzs6xiL7bnLBj8gczu46wbT9YRnuHsbATISIiXWRmpwKvcfcpR7NI+npeABzi7qfmmP9CYEGSXomIyDQoaReRThrW/9zCBUguzEx7sZn91MzuMrMHzOwmM/uOmT0uNc96y0xNe7yZfSg5ovaghWHN9u3z2rPN7JNm9mczW2Fml5jZC/P2Pbcw7Jyb2TEjrvMcM/uchYu9rbQwxN0zh7zOWy1cPGVl8v9bM48vJRwyf3Tymr1/e061DiIisj6NHiMiMgUzex6hb+dVhPGH7yGcdL4XYXi463Is5jTChdA+TjhB7XDgLDN7nLsvTc13BrAvcBZwHvAYwkldVZ8wfi5wJ+EiIVsQxjs+28wWuPvfMvO+FXgE8N/A3wijXJ1oZvPc/dhknsMJ22pLwgl6PdkLl4iISA4q2kVEpnYA4cjk3yejLvQcN2D+fu4C9ktObMLMLiD0K/2/JP05k+R9X+AUd39D74lmdj5hWLQq/cbde/1VMbNrCP1iDyYU52mPA56QGsLuvwgjNrzXzL7o7re4+1nJSaibuPtXKm67iEjrqXvMCJITPfpd8GC6y/2Rmd1jZj/ITH+hhYsTXGFmF1kYFxUze24y/SFLRqrps8zZZvZDC1dfu9rMjk89NuXzRWSS3jBeL7eJob9G9WlPnUTk7pcRUuodU/Psl/z/yfQT3f1sqk+oT8jcPz/5f8fsjMBX06Naufuq5PkbMrEOIhIJM3t70j1ty+T+IjO7OOmq9/Ycz/+MmS1P3d/YzL5hZkvM7NJelzszm2lm/2PhKqdXqjtcuVS0j8Ddv9cbhqlkH2P9Cw8AfA7452RM9a8B702m30QYleFrUyz34+6+iHAltr8zs31GfL6IBCcRLkT2WeBuMzvbzA4zs61GWMaf+ky7m9AVpecxhLGgl/SZd8rLcE/TpPa5+1+Tm/2u8NdvB6I3rNsOZTZKpM0sqLQWM7PtCSPC3JSafDdwGKG73lTP3xWYm5n8r8Ayd19I2GH/SDL9DQDJVU5fBHyi6vXrEm3IRHLy1mILwyb+3sy+amZ7JSdjXW9mz7Aw5vBJyfynmtmJZvZLM/vTdBLr5Gqm2T6jEC4ysFlyew4TA/ovdfffMflSu9llrkiuyNpLwX5DcuWyPM8X6YBhJ3VOStOTArZ3yfLPAJsSfqiuM7Pd1396X2sGTLc+t6czrFfu9Zr0pDAkXD/TuQqjiGQk9ca1ZvZZwm/zGjP7iJn92szOS+qNC5PaYv/kOTub2a+SI++/s3Bl9rxOAN5J6u/W3f+SHO1bPUVbNyAEi+/MPHQA4TwdCMM/vtDMDNgJ+EnvNQjn/+yKlEJF+2QLCRdJejKwiNCXcw/CpWrf3Wf+RyaP/x+SMU3NbFMbfNWunUZsz+sJJ4LdQkjiC6X8Fq4Wth/JH5KIACFpApiXnmhmswh/25O4+xp3v9Dd3+PuzyEcwXo4E0fAynAD4Xu53w/y43Muo+96JcpKwft9lz0h+T+d2GtMYZH+Hg98ObmaJsCF7v50QoD3AUJK/VLCieEAbyR0sduFUAT3zif5+YB6Y6/k8f2BW939yoLtfAvwveSCSmnbATcDuPtDhC6EWxDGbD/AwkWXHgM8nclXVJVp0Imok93g7lcBmNnVwE/c3c3sKmBBn/nPSq58do2ZbQOQjLKwS0ntOQLY190vNbN3EPq5vn6UBST9b78OnOju/Q7Pi3RVb8SXvQhpV88RZAINM9vS3e/KPH8x4eIf/Yrjor5PODH1CODQ1Ovvy0RRPJS7/83MbidchttSJ77uALykpHb+s5l9IHUi6sykzWuA9Lk5y4HN0+0QEQBudPdLkturCBdBgjBC1YPuvjpTe1wMvMfMHgV8x92vB0gChL7MbDbwHsKVT0dmZtsCBzJxddNJD/eZ5oSrgj4BuJxwkblfAg8VeX1Zn4r2yR5M3V6bur+W/tsqPX/vcsabAj8fsPyDCYfVeyMxvM/dv9dvxqSv7FPc/dJk0jeY+KMexcnA9e7+qQLPFWmz8wiF9/vNbAtCyr0H8CzCSC9pX0h+LH9M+CHaBDiI8Pf85RLbdDZh6MU3JCeM9YZ8PBT4HeEoYB4nEdK6c8zsLMLwlG8kXLJ8txLaeR1wqZl9npAMHpws9zh3vzk13yWEI5EnmdkvCUX9+ZkReES6KH211NWpndp1tYe7r+2d+O7uXzOzS4EXA+ea2evd/Xwz+znheyjr7cAdhO+PK0PPFR4F/MbMnuHut+do41MJPRCWJM+fbWZLkn7stxAS9FuSNs4B7k7WY90Qr8nf/fU5XktyUNFespxJe54kfhkwx8IYztcRDpWNNHqEmX2A8Ic0Ujov0gXuvsbMDgBOJIw7vopQlD8P+EVm9v9HOHn7NcBWwH2EEy9f4e7fLrFNbmYvJ1wO+5+AfQjF+ksJlw/P24/1I4S//X8hpGTXEE4cezrlFO2fIZxv81ZgPuEEt8Pd/dOZ+T5F6JLzCsJOwwzCeQEq2kVGkBwp+5O7n5jcfjJhB3hg0p7YOrWMpcCufY4a9uXuPyRcj6H3/OVJwQ7huhWvIRwBeEXSFk/SfXP3+83sRcBD7n5NdtlSjOmIZZAMV/QDd39icv/U5P63eo8RzrLe1d3fkn48mX+5uz+84Gv/nNCH/uHAX4F/dfdzzazXn20toYh/nbv/ycx2I1xsZXNgJXC7u++cLOsKd98lSQVvJiSJvSMCJ7n7KcOeLyJxSg6Vb5SMCCUiDdWn3lhXP1i4kvFyd/94+jEzOwp4FeHE0duBg9397j6LH/a6S0mKdjN7BKELy2aEGmM5sJO732dmZwOvd/fbMs9Pt3MWIcx4KuE8mlcm9ckCwtHCtcCthHrmxlHaKYOpaBcRiYiZbeLuD2SmvZgQHJzk7m+tp2UiIlInFe0iIhExsw8T0qsLCCMy7AK8jtAlZ5f0RY1ERKQ7VLSLiEQkGSnmSMKwinMIh57PB452934XXRIRkQ5Q0S4iIiIijZGMQ3+Zu99bd1vGSUU7sOWWW/r8BQvqboa0zAxdcFYaYq2usyciFfvtr399l7tvNd3lmNl2G8AtWwO3uXfqSswa8hGYv2ABv7j88rqbIS0zkxV1N0Ekl1XMrrsJItJys81KGUVmK7jl4YQhdMxsrrvfU8Zym0BFu4iIiIhEz8y2m024qtMaYHUYDrszabuOiYqIiIhI9LaCWx5BqNK3AO4hpO31tmp8lLRXLO+VD+ZV2gqpQ94uB8O60eRdxkhX2Mghhs9j0XWKfVvEsG2zyt5mw8Sw/mWsb4zf7TFs27Rxfq6yxrUtxvlZyqutn7l0yg4hdd6KbqXtStpFREREJGrplL2na2m7knZCv6g6EwEieP2miy1hqkLez0gZCWCRz2PZyyu6jJtuKuHF0uZPfxHpbdP1v/XYPltFn5/7c1bC5yet6r+ztqhrWxR53Ri/s9Ji+MxlU/aerqXtStpFREREJFr9UvaeLqXtStqBh4groehCaly2Ye9fG7dnnf0oB23Puvqgl55SFVz+/CHpVkzfL3Ur8vcYQ7/hop+zMj4/aVV+lmJIVKtQxm9A2etf5fdWkz5zeQxK2Xu6lLYraRcRERGRKA1L2Xu6krYraQdWreq/Z5p3L7Rsde/Vjir2JLvJKfzdA25nVZ02T3qt8b1Uo4zzPSiiru+zrHH1+c27vrG8bzG0I9uEWD4z03VFBNs2RjF85qYyVcre05W0XUm7iIiIiEQnT8re04W0XUk78OADsOSq9aenpy180vRfp4zUosrEqagmjzE7qO2xJ/BlJCT9PvNSvaWLp7+MBYumv4wY6DM4unFus0GvVcbvYdXGtZ3GuS2GrdOwdjT17yxvyt7ThbRdSbuIiIiIxObgrRit+u6l7W2lpB14cGX/9CudZuXdUx22t1tX/7Fhrxtjn8V0+l1X4t3kfvBZRVKWMtLgIpqWIBfZTkvumP7rppexcJt8z0m3teh2blKaV7Q9daXIZWy/uv5uiyi6ndvyOStb2e2I5LM0c9EMeOyImfkta6ppTAyUtIuIiIiIRE5J+xBFkqmifc7qEnsKH2PiHcORgGHyJi5Fk5RBSXHexHeYcfb3rjJJKiNNL+O1Br0nw9Z92PaLIc2r+mhMLMnpIDEkoNk2FHlPYt/OXRDDZ0lGo6RdRERERCRyStqHKDvRqStZKJrwxzhSTVqMKXzZ0ut4/g/zPafs9CRvapydr4zkvYi86z/ONLxJ8m6/uo5olJHyVqmL6WUZ50tUqezPdJW6+PmR/JS0i4iIiIhETkl7Rgx72mXLJvxV9q2PJZ2vst/5OFP8dLqeTmCakBI3oY2DLLk333wL51TbjlGN8+hGLIngoHYU/S6PZb2aqq7rEJTxumUfMdBnScqmoh3YeFY7i/VBYrtgRuwnw8YidzeVnAVnVmwFaNmKbpciy6xrWxbdUaqrK1OVVDA1VwzvXQxtEMlS9xgRERERkcgpaWfwxZWq1qV0f9zaciJq+ujHsM/opIvtZFLevAlzkSQ6xnS+ikQ9jxi3RV5ld2VqY3IvIuO39TxYMHu058z8czVtiYGSdhERERGRyClpr1HZJ1AVEcMFn9RvfbD0ttnrwInb2fMS0p+ZoYl8ySl0Xal2LJqcrldp0pEfpe4iIqVQ0i4iIiIiEjkl7R2XTmzHmborXc9nUt/8IdusjAt3pVNzJchSllguuiUi0nRK2kVEREREIqekPUJ1XaZ7WFpbRgo/znT97gHTyxhVpsoLN43TsP7obRlJRuKj/u4iIsUoaRcRERERiVwjknYz2xB4O/CvhJ69dwJnuPsRqXkMOAp4E7AlcBlwmLtfMf4WT0+M47eX0fc9feXTuvq0ZxP4Ikn5ONP1K1LbLP0e5B2zfZyGpfNtTeFjuyJqEyhdFxEpphFFO/A/wAuBY4HFwPbATpl5jgSOBt6RzPM24Dwze6K73z7GtoqIiIjING22JWy55WjP2fCv1bQlBtEX7Wa2N/BK4Cnufs2AeWYRivYPu/tJybSLgaXAW4D3jqe1o4kxUe+aIv3Tx9mn/bwzJm4XTdDHNZZ619Pl2NdfCbf0U/bvUB1XFxfpiib0aX8dcP6ggj3xbGAz4Ju9Ce5+P/B9YJ9qmyciIiIiUq3ok3bgmcD3zOwk4NWENv8IeIu735bMswhYA1yfee61wEFlNKKLqXjZ47ZrbPbxyJusR58MR96+WChBH13e7/O2pMbj/P2K8beyLe+jSBOS9kcArwV2IXSTOQR4OnBmcvIpwObAcndfk3nuMmC2mc3MLtTMDjWzy83s8gfuv7OyxouIiIhIPMxsoZn9t5ldaWZrzOzCzOOPNLOPJY8vN7Obzew0M9u2z7K2M7Mzk/nuMrOTzGx2Fe1uQtJuyb8D3P2vAGb2Z+CnwAuAnyTz+YDn9n3M3U8GTgaYv+OuHmM6MA5NG3+9bLGPs57+XObt055NqMfVp30YpeblGle6XvX3YpEEtOupsYyua++jjizksjOwL3AJsF6wSwiHXwqcAlwKbAMcA/wyGeBkOawb3fBcYBWhZ8dc4JPJ/68qu9FNKNqXAX/qFeyJiwgbaCdC0b4M2NTMNsik7XOBFe6+emytFREREZGYfd/dvwtgZt8iDBWedhGwyN0f6k0ws98AfwBeDpyWTD4QeAKw0N1vSOZbDZxuZse6e7bb9rQ0oWi/Fti4z3QD1ia3FwMbAAsJG7RnUfJY9MruP162JqfpabEn61npz0XRsdmblHKrf/Z4xJI8xtCOsr97h11Zuq3q+v3q4rbOo+oRgWL4u50ud187xeP39Jl2nZmtALZOTd4HuKxXsCfOIgTLe7P+uZbT0oQ+7T8Anmxm6b2g5wIbAVcm938J3EfY4wEg6U+0H3DOmNopIiIiIi1kZk8GZgPp0QzXC4fdfRXwx+SxUjUhaT8ZOAz4vpl9CNgU+AhwnrtfBODuK83seOBoM1vGxMWVZgCfqafZQewJeloT0vSmJeXTlX5P0ulGNvnIm1DXdbXUMhL0KtOdvH1Ah7Uhxn6kZW+zMr7PppuOxv6dGnv72kTbekKVRx3qTNY32gpmbT/ac+zqMNgIcGhq8snJuYyFmNkM4NOE5PzHqYc2B9ZL5Qndtjcv+nqDRF+0u/t9ZvYC4ETgdMIhh+8CR2RmPZ5QpB8FbAFcDrzI3WsqU0RERERk3NKDjZTkw8DuwPP6nCc5aCCUftOnJfqiHcDdlxDO8h02jwMfTP6NZONN6t9jb0LK3WV1Jfw33VTu8sbZZ7xJ/R7LaGuT1jdrnN9/Xbv+Q9l/w22V933U9uyv7hqmzczs34B3AP/k7pdmHl5GGPQkay79E/hpaUKfdhERERGRsTKzlxO6Wb/T3b/RZ5bFZPquJ9cG2oEKBkJpRNLeJLGnPlnDEuS7I2hDk19rnOpKeWNPd4r284x9vcYlxu+zMv6Gx/XdBnFuwyZr+/as4kjCuLZZ246CmNmewFeBk9z94wNmOwc42Mwe7e43JtP2J4x6+KOy26SiHZg5s54vgtgLyNjbN0yT256WLjrTJzoWLdKLFKOx/0gO+6GIsStG+rM5zuKxalX+zZVdqJfd1qLvY+zfUzF8PmPfRsMU2X6xf98O06S2JyMM9rpdbwdsZmavSO6fDTyaMHTjYuAbZvas1NPvdPc/Jre/BbwH+I6ZHQ3MAU4Avlb2GO2gol1EREREumVr4IzMtN79xwDPJBTgTwF+kZnvNOC1AO6+2sz2Bk4Cvgk8SBg05R1VNFpFO2EjNHlvvqu68J7lTYqrTNBj3M7pBKvOdKfItpk34HZWDClnUTF+ZvK2qch2r3p9x7U9s+se4/s4yDjbmvczUuVnbhRNeh/Hxd2XEkZ4GeTU5F+eZd0CvGTajcpBJ6KKiIiIiEROSTuwAf33RMve+82+RoyJjsQlnSKnb2f7cRdJm8v4LFX5eRz291HX34FOmh7dTFaUurxVzC51eVkxbvdxbcMY1z1GRbZTLN9nTXuPZ24Js0b8fZsxs5q2xEBJu4iIiIhI5JS0D9GWfooyWNkJVtkWpm5PGv1ihOShqZ+zstsd+3st+eh9nL4ub8Oqj9QM0tTvYYmLknYRERERkcgpaQdmsLYxyUNdKUFTtk/bpLf7I2psR2xWcmuB50jMZrFdqcsr8hmR+JT9uWjrb1n68172NpN4KGkXEREREYmcknZgLavHksqUsffb1pRA+ruHX/WdXnWSEmPqqeS0OYp8fsb5/jb5s9S1FLXoezWu7RTjZynGNkk5lLSLiIiIiEROSfsYVb3327UEpgtWcmVNrxtfUlPXtpAJs3hKrvli/PyMyyif07zbc/Lyu7ttIf/vXOzbqa72qU5oNhXtIiIiIhKdDefBxiPuZ1iLL66koh1Yy6pa9npj7DfcVG1ND1ZyYd1NqFwX1rFJZrHnwMd0tKNcXd6eRY4yQLd/58qg7dds6tMuIiIiIhI5Je2A80DfxKNoEpCX+rSVJ++2bNq6P8glUS9PmmNjnpVrPh356G/YEYhYlP3eVbnOXT7KIFKUknYRERERkcgpaR+i7CSg6uQ+ryoT/tiT7GHrHmPbb//e7etuz5pP6vZZNbRGYjFr3ujXx636KEveJD9vO/Iub1y6eARinOtcdqpfpO1NOJoi3aakXUREREQkckrax6ityX1ak/uWp9seS/vuvXji9sqbJm7P3WP8bWmK9BGJsqXfgzIUbevKu2+feqZhr1sgqZ9K28+/qDr5r2t9YzmiEcORjBja0GQ6UlE9Je0iIiIiIpFT0t5gRZP7GBL6pvUtr8s9F03cTqfrZSe+TTMsoW7Stsnb1rKPHkw3qY9JFUcN+okt+S9LkfWKJZ2XuFRxpGLmVqN//81o8cWVGpe0m9l2ZrbczNzMHp6abmb2bjO72cweMLOfmdkudbZVRERERKQMTUzaPwYsBx6WmX4kcDTwDmAx8DbgPDN7orsPjZXW8rcp9xDb1FdrUEIfQwIvk12e6tO+RyptmLP75PmalC4PkzdRafL6FusqPgsAACAASURBVEnNq17fKs8DqFqVRw3GleI3TVuPOpRNRySkbI0q2s3sOcDewIcIxXtv+ixC0f5hdz8pmXYxsBR4C/De6b523sM+bSruJV7ZIu7Bmq5MvXHJPZnqKsbHWbTW1SVmmC5s9yJi7EakHYnmKHvnRjsB0pjuMWa2AfAZ4P3AXZmHnw1sBnyzN8Hd7we+D+wzrjaKiIiISPzM7JVm9puky/WtZvZlM9s2M09UXa+blLS/EZgF/Bfwz5nHFgFrgOsz068FDqq+aROyiXzsybu6xMRtwYAkO5usl5GUFkk9x5nwl53qpxXdfk0aXnKYuhLvWIbQbJIY0/80HQmozrDkXin8aMxsf+DrhJryHcAjgQ8APzCzXd19bTJr4a7XVWhE0W5mWwDHAa9y99Vmlp1lc2C5u6/JTF8GzDazme6+KrPMQ4FDAbadv0k1DRcRERGR2BwM/Mbd39KbYGb3Ad8FHg9cO46u16NqRNEOfBC41N3PHjKP95lmgx5z95OBkwGetOvm/Z5bSOzJeow0xGMciqSe40w2i6T6VabzEGf/9CK6th7DxL6OsSv7SICSe6nIRsC9mWn3JP/3ase+Xa/NrNf1euxFe/R92s1sZ+B1wLFmNtfM5gKzk4fnmNkmhER906Tfe9pcYIW7rx5fi0VEREQkYl8CnmNmrzazzczscYTuMRe4+zXJPMO6Xi8aX1MnNCFp35GwR3Rxn8duAb4IfA3YAFgI/CH1+CJCH6ROi7HfeozpeoxtWppKl7eMcJjD2BPadDpfdeo+zKDt1LRUd9j73bR1GST2z3TXZJN7Je/dstGWxS6ulO4CnTg56WEBgLv/0MxeS6ghT0sm/xLYP/WckbteV60JRftFwPMz0/YG3gXsC/wJuBG4DziQsKeEmc0G9iPpAiMiIiIi7ZfuAt2PmT0f+DzwaeAcYBvgGOBMM9srVaiP1PW6atEX7e5+F0weksXMFiQ3f+7uy5NpxwNHm9kyJs7wnUEYJrJSMfZjjzFdl9Fdmtq/X5BKAB+VSR6UDk5tWJ/4ulL4Nr1v6XVpQnunqy1HT5omnbx3PXXXiDHT8gnge+7+rt4EM7uCUD8eAHyHVNfrTNpeW9fr6Iv2ERxPKNKPArYALgde5O531NoqEREREYnJIsKQj+u4+x/M7AHgscmkxUTW9bqRRbu7nwqcmpnmhFFmPljW68SYoKc1LU2Psc94jG1qCyWO3dW11D0tm8B3bf2r1rV0XWl6ZW4EnpaeYGZPADYhDOkIoY97VF2vG1m0i4iIiIgU9HngBDO7jYk+7e8jFOxnA7j7yjq7Xvejoh2YwaZRp+pNS9R7YkyxY2xT1caV9ClRFJEyKE1vpqnrqLPG0Yy8TgRWAW8C3kgYo/0i4Ch3vz81X1Rdr1W0i4iIiEhnJF2qP5f8m2q+UrteT4eK9kg0NU2HbqbXddhySJJdV8rd5HS9znHb82jatm1ae6ui7TBYGxP0tqTkw8TcE6FrVLSLiIiISHRmFri4ks2spi0xUNE+Rk1O09Oalqw3rb2xiy1JjD0xHya2bZkVe/uq0MV1zqONKXlWF1LzNCXozTOj7gaIiIiIiMhwStoLaktqPowS6rikE8BsGpj3ypqDltcETUrUY9+2sbeval1f/zQl6M2hZFyUtIuIiIiIRE5Je0YXEvS0tqbpbVmvZ24wcXtYOtjk5LCuBL3J22yQNq7TKNqy/l1Iv7PakoYPopRcyqCkXUREREQkckraAWOTVifsbUmdJQ6x9y1vS9qa1qZ1asu6dDEN72l7Kt6PkvIJba6XYqekXUREREQkckraW6TLiXpb131BarViT7iLak3yqvWIWteS8S6k4V1Lv7uYcG+06VxmzZs10nNmbHBXRa2pn4r2GrW10ByXLmy/LSssoNpanKU1bR2b1t4iVDw3X9eK5WG6WEhLfdQ9RkREREQkckrac+pCqtsEbXkfVjF73e2ZrBg437ALKo1LF9LfvJq8LbqQcLcx1c7qesqtZFu6TEm7iIiIiEjklLQDM5jZmgS3bbr+vsSQtJetLeuR1YUke5gupNyDtDX9Vqo9Hl3/nZP8lLSLiIiIiEROSbvUTinDYHN2n7jd1oQ6ryYn2V1OoYdpa0Kdl5Ls/vSbINKfknYRERERkcgpaZdSKSHpLz1azCgm9WlvcNJcNiXXE7qeVqd1PbnW969Iu6loFxEREZHobMTj2JhtR3qOcUFFralf9EW7mR0I/AvwdGAO8Afg4+7+9cx8bwDeCWwPXA28091/Mubm1kopS1yKputpc3eZSNe7ni4rUR5d15PnIvQ9Wj+9B9VZya11N0GmoQl92t8GLAeOAPYHLgC+ZmZv7c1gZq8EPg98GdiHULT/wMyeOP7mioiIiIiUK/qkHdjP3e9K3T/fzLYlFPOfSaYdC5zm7scBmNlPgacCRwKvmuoFZrCR9uylFGWk62npdH0u/77u9kquLPV16qQ0eDz0HRc3vT8yDvqcrc/MtiP04ngYsKm7L0+mG3AU8CZgS+Ay4DB3v6KutkZftGcK9p7fAgcAmNkOwONgoqJx97VmdkZ6mkhZyi7Mh0kX6sO+bPN+ERc5NKov+fbTeywiHfYxQo+Oh2WmHwkcDbwDWEwIi88zsye6++3jbWLQhO4x/TwbuCa5vSj5f3FmnmuBeWa21dhaJSIiIiKNYGbPAfYGPp6ZPotQtH/Y3U9y9/OAAwEH3jL2hiaiT9qzzOyFhJT9dcmkzZP/78nMuiz1+J19lnMocCjA/Pnbl99QicY4k/GypRPQQbeLLk+kbZr8ty71mMmKupsgNTGzDQjdrN/P+jXks4HNgG/2Jrj7/Wb2fcK5k+8dVzvTGpW0m9kC4GvAd9391MzDnp19wPQw0f1kd9/V3Xfdaqsty2ymiIiIiMTtjcAs4L/6PLYIWANcn5l+LRM9PMauMUm7mc0DzgFuYvLJpb1EfS5wb2r63OT/7N6TTJPSrPGZwY7rbq+qsR0iIm00rt8zJfrjle5NkTjZ3U9OPb4FcBzwKndfHc45nWRzYLm7r8lMXwbMNrOZ7j72n+Upi3Yz+1nOZa1097+fZnsGtWE28ANgJvBid78/9XCvL/si4MbU9EXA3e6+XtcYEREREYnbRuw88ghjM7gc93tOBk4eMtsHgUvd/ewh8/TrqTG0F0fV8iTtuxEOIQxjwKen35w+CzbbEDgD2BH4O3f/S/pxd/+TmV1HOEHg3OQ5M5L75+R5jbXMUHosUVqSuj2vtlaISCz0PdBMddYYSvknM7OdCedFPtfMer0yem/QHDNbQ0jUNzWzDTJp+1xghbuvHl+LJ+Qp2n/p7qdNNZOZHVxCe/r5LLAvYfjGeWaWvizkb939QeAY4CtmthT4BfAaQpFfVZtExuKmm1J35tfWDJHSqOicnrtrel29b82lUHI9OwIbARf3eewW4IuE8yc3ABYSxnDvWcT6oxWOzZRFu7u/MM+CquoaA/SW2y/Jfwyw1N2/bmYPB95FGFPzauD/uPvvK2qTiIiIiDTPRcDzM9P2JtSQ+wJ/InS3vo/Qa+MDsK6r9n4M73ZTqehPRHX3BTnn+wLwhWpbIzJe550xcXuvA+trh8RrfsOOwNSVFJeta8nzsPeta9tCmi25aOeF6WnJ6IQAP09dEfV44GgzW8bExZVmEIaJrEWuot3MLmB4p3vPm8iLiIiIiETueEKRfhSwBXA58CJ3v6OuBuVN2r8yYPp2wGGgDlMiVViS+mrYKz39qrE3pbMWPqnuFgw36byHCDXtSEBeZRwxaEtCnd4WbVkn6Zbk2j+nZqY5YZSZD9bQpL5yFe3u/sX0/WR8y6OANwDfIFxNSkRKtiR15YF0ob60ttNgpm9BbZelKCb2HaSu7VS0aSdgUOHf5MJXBbxIdUa6IqqZbWZmxxFGotsGeJq7H+rut1TSOhERERERyd2nfRPgcOA/CJ3393D3qytsV2O15SSrMihlmb500p5O15dketQt3GY87SlDE44SNOloQJEjAbGn88PkTe6bnMi35aTP7Ho0qe0Sh5ksGPniSsYmFbWmfnn7tN9AGK/yo4SO+NuY2aQywd3PL7ltY7MGFdtVyLtN9UU+2Db33pq6t93A+bJFfB2atOMwlenuWMRe9A8r9Jtc0Kdli/smF/Ft0ZadEelPdVT18hbtKwmjx7xpwOMO7FBKi0REREREZJK8J6IuqLgd0mE6hDq6dLcZgIVz6mlHWt60v02J/CDZpD725D0tncK3JXVvk7ae6FllStum7TSIUu5uiP7iSiLSLk3uj59Xk4p0aGdxru4w0qOCVtoi1+gxZvY0MzvXzM4zs2dW3SgREREREZmQN2n/MrA/MBM4E3hCZS0SkXW2ve+21L3BJ6LGri1petMS9LQ2pumgRF1EuiNv0b4B4WTTtSM8R0RKNKlgvK62ZrSmAE9rWjHelgJcBffoutA/W0T6y1uAvwY4mdCd5nXVNUe6SD9Cg9222bZ9p5dx4mkbi29oVgHe5OJbBff06btPREaRd/SYXwEvqrgtIiIiIiIAbMQ8Zo3YNXQGMytqTf2mLNrN7Dh3PzrHfMe6+3+W0yxpOyVM5WpLat6klByam5QrJc9H31MiEpM8SfvhZvYlwKaY7zBARbusox+8cqULxPN+Wl870ppWZA8Se/HdxSJb3x8iIpPlKdofBixh6qJ95fSbI02nH9py3TGn/2HBMpJ1Fdz1aEsBrr91EZHxmrJod/dcY7lLd+nHu35NKsCbVmSnNbng1t+piEizafhGGZl+/Mdn0Cgx4yzSm1xkp6ngbj9tp+bQVUpFRqeiXXLRj2H90kXnkqsmP9aWwnqYJhfdaV37W+ra+ko+Tf5caIdD6qKiXSRig5L2phXpbSm4s5pceAzSxnUSKVNdfyPaWRAV7TKQfry7pa2FdR5t/ay3db1Euqjo37OK/fYYuWg3sxcBrwS2dvf9zGxXYDN3P7/01ol0XN5RYrpccGd1oVDtwjq2xUxW1N2EUqxidt1NkILyfl/EWNzPZIsCF1faqKLW1G+kot3M3gr8O3AK8Ipk8gPAicCzy22ajJsKgeZoa5Guz+AEbYvqtKWQHqfYt5l2KqYv/Z0TYwEvoyfthwMvdPelZvauZNpi4PHlNkvGRYVB3NKjxNT1harPSLm0PcsVezEp41HG50CF/wQV8HEatWjfFLg5ue3J/xsBq0prUUFmthPwGWB34B7C0YBj3X1NrQ2LkIqG5tN7GAe9D9VRMS7jNuwz1+WCPvs915Yivol146hF+8+AI4EPpqYdBlxQWosKMLPNgfOAa4ADgMcCnwBmAO+tsWm1UTHRDk0bJaap9PdSPxXpIs3QhhS+qXXjqEX7W4Hvm9kbgE3N7A/AfcB+pbdsNG8ENgFe5u73Af9rZpsBx5jZR5NpraDiQrpKn/32UaEu0mzDvpcjL+gbWTeOVLS7+5/NbDdgN+DRhK4yv3L3tVU0bgT7AOdmNvLpwEeA5wHfH/bkDVBBIHEadMKpPq/SVCrUpQm63B2mLJH/Tk2rbqzLyEM+ursDv0r+xWIRMGnISXe/ycxWJI9FufFFpjJvwG2RmKkwlyZSod4pjawbpyzazez9eRbk7u+bfnMK25xwEkHWsuSx9ZjZocChANu3dfw8aTwV6jIdKp5FpIvSNV7iZHc/OXV/5LoxBnmS9u1Tt2cBLwcuA24E5gPPAL5dftNG5n2m2YDpJG/eyQC77vo014+bxEifS5H2KDvJ1fdDubQ947OajUf+u1nLjEk13hAj1Y0xmLJod/dDerfN7HTgn9z926lpLwMOrKZ5uS0D5vaZPof+e1IijbCSW9fdHvWqcCJtpC4ME7QtRqfCXBKNrBtH7dO+D/DPmWnfBf6nnOYUtpjQB2kdM9seeFjymIiINJSKUynLsM+SCvpOaWTdOGrRvgR4M3Biatq/AX8srUXFnAO8w8w2dfe/JdMOAh4Aflpfs0SmJ520j5NSfambCnUZt/RnTgV86zWybhy1aH89cKaZvRO4FdgOeAh4WdkNG9HnCRd5+o6ZfQTYATgG+OR0xtocZ8GkIkn6WcmV626P8zNS185C2fR3JSJFKJFvvUrqxqqNOk77b81sR+BZwLbAn4GL3X11FY0boV3LzOyFwEmEYXruAU4gvAFTWsvq2ouUIq+vgkRkuLL/rvU3JyLSX9111CimWzfWpcg47auBn1fQlmlx92uAF9TdjnHK+weiQqO5VnJh6t6+qemT33u9x+NR9EdJ78/o0mmmusqIxKlJhXpWE+vGkYr2YWO21zxOuwwx7I9KxUTcHuSSdbeHvY9N/uLMq8mf1em+P01e9zKU0R1Bhb/0o64uU+vC70tTjJq0b5+5/wjC5V7PLKc5IpJHuq87wCyeUlNLxqfLha+Oqk1f2cWZdgKqo0K6firU4zRqn/ZDstPMbG/gn0prkYyVxgGP28q7b5+4M+TyqNkivoi2F/6j/Ag19W9Bxf34qLCUpmtCYb4cuHvE56ypoiGRGLlPex8/Br5RwnJEpEZlFP6DNG2HoO0nhqu4F2mvJhTjUsyofdp3yEyaDRwM3Fxai6Q2OrkxPitvSt1JJe2TT1CFWew5juYUVvYOQYw7AW0shMv48W/S+oqMiwprKaLIxZUcsOT+CuC3wKvLbJSIjCZbxE9XG3cCYin0257iZ1VdnDR520h5VARLF4xatL/L3T+WnWhmbwN+U06TJBbq716/dNK+cpeJQjU9qgzAxjyr3Ndt4U5Ak/v9K/EerEnFWozvQZO2n0jXzRhx/qMHTH/vdBsiIiIiIiL95Urazaw3+PyGZvZ8JrrHQLj069/KbpjERal73LLJe1rZKXwRbUnui6b1MXTN0dVh66dUW0SmI2/3mC8m/28MfCk13YHbgbeW2SiJm05YHZ9JJ6Kmp6eHggRmzXvEwGUMK+iLaMtOwDgL/+l2zYmh6M9SX3URkfHKVbS7+2MAzOzL7q6TTmUSpfDjUXZaXVRbdgKaVPg3uT9+UXWl0voOE5FYjXpxJRXsMtSwH1r9GFYnm7wPMiyRr8t0dwLqTP6LFP5N69ozSNN2AvIa586CvhNFhitycaWHqmhIJKYs2s3sue7+s+T2CwbN5+7nl9kwEREREREJ8iTtnwWemNz+4oB5nHBCqshA6kZTv7yJfFqM6Xxa0aS+Sd1yYhgyM6uLXXbK1rUx+0VkeqYs2t39iam7nx0wTvt/lNoqaT0V8KNLF6fZE1Rnza/udYsU+tCeYr/JJ97GWOynNfkiWXXRd6dId416caWjgfWKduA9wCem3xwRaYtBxX7sxXxW7MNpDjOo2I+9mB8mb6HfheJeI3mJdIvGaRcRERERiZzGaZdaKBGaniq7w0h+6RQ+9tQ9LZvANzl5HySbyLcxee/692jX139cdFGweGicdmAGG035x68P7fTpC3Z0TSrOm9btRaSJ2vI92pb16AK9V/GYMcrMbSzYRURERERiN+qJqJjZNsAzgC1J9W139y8NfFIL1LmnGVvKr73u+KRT7uwJoErAp9akri1FtbELTFYbu8BAc79zm9puicddd8JNN009X9qqVdW0JQYjFe1m9hLgK8D1wM7A1YQx3C9icl93KZG++Lor3T1mUmE5b/CIJk0r0rtQMJetCwX4ICrM49aW9RCJ0UjdY4APAIe4+1OB+5P/DwV+XXrLREREREQEGL17zHx3PyMz7TTCCDJvL6dJItKT90RUpdXj0eWEexRtTcOL6ELy3IV1FInBqEn7X5I+7QBLzWx34LHABuU2S0REREQkLmZ2uJm5mX2rz2PbmdmZZrbczO4ys5PMbHZZrz1q0v4FYA/g28AJwAXAWiq6GqqZbQb8B7AP8HjgAeBi4F3ufl1m3jnAp4CXEHZGfgAc5u5/raJtIuOQ7p8+LOVVAtwtSrLHQwnyBG0LETCzrYH3AXf2eWxD4FxgFXAQMBf4ZPL/q8p4/ZGKdnf/SOr2l83swqQxryijMX3MB95AuLjTe4DZwFHApWb2ZHe/OTXvNwiF/esJOxIfAc4CnlNR20QqN4cj190eVqipiJOYqMBrJr1vIlP6MPBDYPs+jx0IPAFY6O43AJjZauB0MzvW3a+f7ouPPORjmrvfZGZ3EArq9023MX3cADzW3R/oTTCznwM3Aa8Djk2m7Q78A/A8d/9ZMu1WQnG/l7ufV0HbRERERKQDzGw34B+BRcBX+8yyD3BZr2BPnEVI3vcmjLw4LdMq2lNs6llG5+7395l2t5ndCGydmrwPcEevYE/m+5WZ3ZA8NrRoX8sMVlFOl6OZrChlOSIwOUFPp2BKxETaRX/TMpWy6hQZnZkZcBLwUXe/NdxdzyLgmvQEd19lZn9MHpu2sop2L2k5UzKzrYCFwGdTkxcBi/vMfi0lbSgRERERGZ8/3wRLrhrtOQ8+AGZ2KGFI8p6T3f3kaTTlEOARwMeHzLM5cE+f6cuSx6YtV9FuZi8Y8vDMMhoygk8Ay4HTU9OGbagd+i0k/YZuPz/nuHo5dHFPWEcXqjOXZ9TdBBERmaYu1gZ1Sgr0gUV6MnjJI3MsZ3Ey74cIg5s8MNVT+r3cgOkjy5u0f3GKx3NfZHaUDdXnuW8inIH78j6jwoy0odJv6NN23XVsRwpEREREpFYHEkZEnIoB7wZuBn5sZnOT6RsCGyX3/+buawhB8dw+y5hL/2B5ZLmKdnd/TBkvlhhlQ03cMdsf+AxhuMczM/MuA7bqs4zSNpQMFnuC0OQjAbFvW5GYNPlvXeKi7952c/dTgFNyzv54YFdCrZm1jDBK4UWEbtqTumSb2UxCj4/PF25sSll92nMbcUMBYGbPJnSH+by7f6zPLIvpP7TjIsKZu9JhZX/5jrMwWJK6PW9srypt1IXPjwotEanAewnXAUr7FHAv8J9Ar9f9OcDBZvZod78xmbY/sDHwozIaMvaifVRmtjPhQkk/Ag4bMNs5wNFmtoe7X5Q8b1fC3s05Y2moiIiIiLSKu/8+O83M7gHucvcLU5O/RRgC/TtmdjQwh3Ah0q+VMUY7RF60J1ee+hHhxNMTgWekhtm5z92vAXD3i83sXODLZvZ2Ji6udJHGaJeyDUvzyk7hb0qfLTJ/wHSgyLnU2WVMd3lSv2Fp+t1ja4WkdeEIh4iAu682s70JQ0N+E3iQ0EvkHWW9RtRFO7AT8Kjk9gWZx34Kk67d/krCHs2XgBmEdH5QMi8iIiIiMjJ333PA9FuAl1T1ulEX7clhh1wXbnL3ewjjaB4y6uusoVsplJKf6qRT+DJS9/POmLi914ETt7Pj1o46ju1Uiixv4ZPKbYPkkz4q0qXvsekY53dg194T/b60Q9c+t00RddEuIiIiIt10562wtN+lM4dYtbKatsRARXsHjWsPuuuJSxmp+4+uS905Y+BsURj1i7WfBR2/fnGRoxXDzk1oi7LPsWhaitik79KmbVuRJplRdwNERERERGQ4Je1SmaKJS5NSpbyyI87kTd6X3Ju6fUeZLZKylXGUoOxzE2JQxrkOTT6aUMZRgjam1238nhepmop2ic6wH6iufdE/9ebL1t1eMme3GlvS38I55S6vaTsmC7eZuF1G96DYFdkxacuOSNGdj7p2OGIftrWNOyJZXfu9kuqpe4yIiIiISOSUtEujpNOZLqQY8++diOm+e289SfuwND3dfaethq5/TUcG0gn/ONV1NCGGE5TrPGLQtROUYz9KkFcXjiakdeE3uW5K2kVEREREIqekHXiIuPaItbeaT5NT97zDQaaT9m3uvbXSNt0xZ7u+05uWppfez77k9S+jfU3q+1/GUYEmnS9QxVGBulL+ui6YFvtRgrYcCShbTHVUW6loFxEREZHo/PkeYMSQ4sGHKmlKFFS0R6iMvdWmJc/Tld1mbVn/h+69Idd8295327rbt222baHXqjrJH2RQwl/UOI8MFEnNYzhyUfbRiGHU7398yk75Yxz5p670P22cRwKU6kua+rSLiIiIiEROSXtLaazzdlj6t6XrbqfT9GHyzleGoql+WtkJf9nJ/TDjSs1j76c/inGl/FUn/HUl+cNUmfLHMIIPqH9/HZT2x0NJu4iIiIhI5JS0d9CgFL5NCXyTR5ZJ+zW/Xnd799RIMjfNmRx9zB/yWJXKTvXHmdyPM5GfrqLJ+Dj7rudVZcqvvvrVqbqvfixJ/iDTTfhj6ItfVAxpvwRK2kVEREREIqekHVi1avQ9yTb28VI/+HZIp+5VqDLJz5vcx9iXPq/Y+9zHmM7n1YX1rTLhr7Of/rhG3akr0a+iL36T03spRkm7iIiIiEjklLQXVKSPV5PTeaXw9as6Qc9rUDti7EtfRiJfttj73Lep/3weede3qes3irwpfowj5+RVNNGPsc99F/rZ374C7h/xO2nlmmraEgMl7SIiIiIikVPSPkZ50/mmJfKxj9QSe/vKMCyFH1cCHkMbssoY3aautD72RD6r7Yl11WPbN2m7lNGvvmlpfZGEPsZ0Pi1vUt+ERL4rlLSLiIiIiEROSXuEhiXysafwsafa2b75MbZxkIfuvWHizgjJdQx90JssndY3rY98LCl8WjqxblK6XLVhSX4bt5P6z+cTQ1pf11VoY2JmmwHvB14CbA3cCPw38Gl399R82wEnAS8CVgKnA+909xVltENFu4iIiIjIYKcCzwXeDSwBng98EjDgBAAz2xA4F1gFHATMTeaZC7yqjEaoaAcefGDqPclY+nQ1KYVvcqoduzJGkskuo8rkPcb+7l2QTuFjTN1FpjIskW9yCp9X3rQ+hkS+rcxsNnAAcLi7n5xMPt/MdgZeSVK0AwcCTwAWuvsNyXNXA6eb2bHufv1029Koot3MDidsnG+7+ysyj1V6SKIJJ2ykC/rYCniIs+tMjG3qmnRBH2MBH2OXmCJiLODVVWZC19e/iK4X9GnDinsV9NO2IeEc0GwntnuAR6Xu7wNc1ivYE2cRkve9Y/hJsAAAHvFJREFUge4U7Wa2NfA+4M4+j1V+SEJEREREusXd7zOzbwLvNLPfAX8Engf8I/Dm1KyLgGsyz11lZn9MHpu2xhTtwIeBHwLb93ms8kMSeWUT+bqS99i70cSYcMfYpqfz9HW3N5zzmHW3s4l0ke4yMabasWhquh5Lgp5X19Llrq1vndIpfNdS9zb58wq4Y8ThVh9aA2Z2KHBoavLJqa4tRbwa+Crw2+S+A0e5+2mpeTYnpO9Zy5LHpq0RRbuZ7UbYo1lE2GhZlR+SEBEREZH4JQX6wCLdzOYAj8yxnF7HoxOAZwKHAH8C9gCOMbO73P2L6af0e7kB00cWfdFuZkboq/5Rd7813F1P5Yck2iS2vu86YXWwBZsuWHd7WDLepNQ8lrbGmKY3LSkfRInyBG2L8VCaPtig/u4d7+t+IPCFHPOZmT0ZeBPw9+7+v8n0n5nZpsDHzex/3H0tIVGf22cZc+mfwI+sCRdXOgR4BPDxIfOMfEjCzA41s8vN7PLl967XTV5EREREWsjdT3F3m+pfMntv9+aKzGJ+SyjIt0juLyYTFJvZTGCH5LFpG3vSPsohiWTeDwGHufsDUz2l38sNmD7p0Mn8HXct5bBFP+k+7rEMG5kWW+ou+ZTRp73oa8VOCXp1upYad219Y6QEvTrZBL7jyfswNyb/P40w6EnP04H7gbuS++cAB5vZo92995z9gY2BH5XRkDq6x+Q+JEEYxP5m4Mdm1jvksCGwUXL/b+6+hjEckhARERGRzrk8+fclM3sfcAOhT/vhTL4i6reA9wDfMbOjgTmEvvBfK2tAlLEX7e5+CnBKztkfD+xKKMqzlgHPAS5i+CGJzxdubAliTNelOZb+bem62zN53rrbMabJRTVpXWJPzJUMT9C2KJcS7/ZRsp6Pu68xs/2ADxCGHt+KkL4fA3wiNd9qM9ubcB7mN4EHCdcMekdZbYn9RNT3Ap/KTPsUYYD7/wR6nU8qPyQhIiIiIt3j7rcDr88x3y3AS6pqR9RFu7v/PjvNzO4B7nL3C1OTKz8kkaUEXWLSpLS6bLGn36Pocjrc5XWvgpLx9lNS3j1RF+15jeOQhIiIiIiMz9z7/8om99460nP+smZ1Ra2pX+OKdnffc8D0Ug5JdC1Bj2HEGI3LPlh6nPZLUqO4tCldjkGTU94mtz0GSqRlOpR2yzg1YZx2EREREZFOa1zSXoWNN2lfwh5Dgj6M0vV8LnnSv6y7ne63Xme6uuTe/tOV+JYrxu1ZJJVeckc9ryvlUqIso2pbXRUDJe0iIiIiIpFT0l6j2NPwIpqcoMfY9t9uv9u62+nktc4UNsYEWIIYE+m62qRkWGKi1FnKoKRdRERERCRyStpzamMqnhVj0lyGJq/XW3brP71oirh0cfG2SHsohZZxULosUi4l7SIiIiIikVPSDsycWX+S3uQ0uC5d2GaDEtGiCVaMydeSq+pugQwS4+dFuqXu32ap1+YPLGPefbeN9JxlLb64kpJ2EREREZHIKWkvKPaUN/b2ST6Dks5Y0qe8n7O7hzw2aF1uumnU1sg4xfIZrNI4P4NN2p7Z7RJ728v+PRz2fdYk2i7No6RdRERERCRyStoJG6GOZFppuEwlnWCV/XkZ5+cv/Vp505h5mfSuyhSnaclhk8T4PZf3sxT756CMbVvk76qM7RLj5yKvJre9Stou1VPRDmxAMz9sM1lR6fJXMbvS5cvUFg6YPuy9H/a+Vf2ZyWNewc9VkcI/97IjL866YJyH6pv4fV+VJm+LJre9rq4kTd5mou4xIiIiIiLRU9IOzGBtFAlkbLRN6reW6/tOX1lweUWeN4vtCr5af2V8rh5RQjt0JKndlCjKqIoewSyiSZ9P1QLxUNIuIiIiIhI5Je0FreTWXPOVnVJKt6Q/Zyu5MtdzZvGUQq816LOa97PeBfp7Llfe9FJJX39F019tz9GNa5tV8X3b5O+trVbcyXYbjjb26vVrVlXUmvopaRcRERERiZySdmAtqytLE5VSjkeTk4Rh0un6Si7M+ZzJ881iz/IaxPi2dYx/O9k2xf65K7INy16nYcuLIfGN8XM2THp7lr39mrwtYhHjNoyxTVKMknYRERERkcgpaQfWsmose6IxpgJtMez9a/J2T6fmD3JJrudszLMGLqPJ8vbp7zJto3LVlVAO+86qsk11JrJFvqeblCA3oa1tfw+mw8wOAg4CdicMYHaIu5+amWcv4A3As4CtgD8Bnwc+5+5rMvMeAHwA2DGZ71h3/8ZU7VDSLiIiIiIy2CuABcAPhsxzKPAw4L3AvsDpwCeAj6ZnMrM9gG8DFwD7AD8Evm5mfz9VI5S0A84DpSVUw0buiDG1KUMM/WaH6dpIP3kT+aLyJvdl96Uv0oa2mu76D3tvyh6lqCtJXJli3GZVj14V4zp3TdnvQcuO/B3k7mvN7OHA6wfM82/uflfq/oVmNhs4wsze7e4PJtOPBn7m7ocl9y8ws52B9wE/HtYIJe0iIiIiIgO4+9oc89zVZ/JvgVnAZgBmtjHwfOCbmflOB3Y3sznDXkNJe8lG2bMsmkiMKsYEI8b0O8Z+8YMS0SrS9LKXWbQP/rgMa1/eNtWV+E9325ZxtKTJKVqR794mr+84lbGdhh+xrud9KPv3Wp+nzng2cJe735ncfyywEbA4M9+1hCD9ccBlgxbWiKTdzB5tZl83s7vNbIWZXWlme2fm2c7MzjSz5WZ2l5mdlByWEBEREZGOMLNDzezy1L9Da2jDTsAbgf9KTd48+f+ezOzLMo/3FX3SbmbbAxcDVwKHAPcDuwCbpObZEDgXWEU4u3cu8Mnk/1dN9Rpr+Vtpidko/Xhj2NMeV9pfVCzjYg9K4WPpB7/y7tvrbkJx86rtg19E1ecFjMu4EvlY5P3+jeG7t6gY35Oyz1+J8f2JsU1FxPj5GWbLFXcx30e7IurMNatw95OBkwfNk3RDeeRUy3L3bCKei5ltTjjZ9HfAh/otOvuUAdMnib5oBz4G/BF4capP0XmZeQ4EngAsdPcbAMxsNXC6mR3r7tePrbUiIiIiErMDgS/kmM+mniXzBLNZwHeBjYH93X1V6uFeoj4387Te/WwCP0nURXuyJ/Qy4DVTnASwD3BZr2BPnEVI3vcGxla0j7IXW+XoGnkNSw9iTOFjT7yrlE5NR0nWV44WUlRu1vzJ9xt9lCAys+Y9otDzBiXy4zzfoMnnGHTdOEeU0nssZXD3U4BTyl6umW0AfA3YGXi2u9+RmeWPwGpgEfDT1PRFwFrgumHLj71P+9MIHfbdzH5hZqvN7BYzO8rM0ns/i8h06k/2bP6YPCYiIiIiUqXPEsLi/dz9D9kHk2EfLyAk/WkHARe7+73DFh510k646hTAfwOfIwxY/3zCVaTuJWwcCB33+x1SWMaATv3JSQmHAmw7f5N+s1SursSg7H6fMSTyMY5GU4Z7OX7d7XQiHVt6Poq8bc8m8jK17FGLosl7Tyx9+2MfiahsdW33qrefUvJ6jOvz1Ja/v36Sk0p3IgzfCLCrmS0H7nT3nybzvJtQV34YWGtm6Q1yjbvfl9w+jjCG+6cIvUL2Tf5NGmCln7EX7SN2/u8dCTjH3Y9Mbl9gZo8CjmKiaIf+nfdtwPRJJyk8adfNh3b8FxEREZHO+kfgP1P335z8+ymsS0J7VzQ9KvmX9nwIe63ufpGZvYIQQL8JuAE42N2HXlgJ6knaR+n8f3dy+4LMY+cDh5jZZsmeyzLW79RPMm1op/4uKrv/YZMT+diT97LT9RgS+rwJegxtHbeyjy5M93yB6Sb14xbLkYGmKmP7VZ226j2OV5vfG3c/Bjhminn2HGF5ZxFS9pGMvU+7u5/i7jbVv2T2awcspvd47+TUxWT6rpvZTGAH1h/AXkRERESkUaLu0+7uS83sauCFhH7tPS8E/ujuy5P75wAHm9mj3f3GZNr+hOF2fjTl67B8WnuIbe3HNSiRLzoCQDqRjyF1L2qc6fzkpDNfahp7Qq0+7c0R48g+TUv/u6bNaatI3aIu2hNHA982s48BPyb0HfoX4NWpeb4FvAf4jpkdDcwBTgC+pjHaRURERJpnzf138NCqG6aeMcXXPFhRa+oXfdHu7mea2asJRfm/AzcBb3b3r6bmWW1mewMnAd8EHgROB94xjjaWnSy0NbmPQex92LPmcOTEnXnHD55xkm4l8kUpyW+mGNP/NB0JEJGqRF+0A7j7V4CvTDHPLcBLxtOi7injohh1dYlpWqGeNnmbHTlwvkldmeZN3omMvcgZFxXpMg5lD7spItIT+8WVREREREQ6rxFJu4xHGWn6MOM8EbXJ6XraoPXIO8xmE3QtAe/a+nZdOnlX6j7YoCOC2mYiE5S0i4iIiIj8//buPViSsrzj+Pe3N24LLNcgyEUCQoAAWgsBkygRCoHiogIuiCZblVCaKgNECcQAgogRTEUJErIVLTHeuC3LpkCWNQRX0FDcBDEgRJE7CwXs2QPL7rLL8uSPtw/0mZ05p2dOz5me6d+nquvM6X675+3nvN3nnaff7qk4Z9prrtvZ9VbGyhT38+Mgy5b/Mqh8zBofx5m/GbqMMezOBpfL8ZyYbmZbJ/OeD493b82xMBufM+1mZmZmZhXnTDsgZjZ9zGK3vySiCo92zGdse5V1bzQoX8JUhnwshnn7kY/tZAdbPTqx29nfbm6/6o+rtNaqllHttD5lZOg9jtsG2fp9nIU9qccgcafdzMzMzCrnxdUvMm11e13VNazpUm16z532MVQhE95tVcmuly0/FrzsJ8l0c9vrv9eSt1/nsnKTmWmu4njsonWqSkY+X48qxrMMdcsO5/e37HHxHvtug6DboxXqyGPazczMzMwqzpn2Gqp6dr3u49hbGZ1t87ecFjFWVrtuVyucre2eMsbF++9jZuNxpt3MzMzMrOKcaQeCFW+NvarCOPaimfDGZ3V3so2yOUterlmc/tbr5bn5G2w5ulzZYwfLHqNbhSzieuOEK5D97rYqxH0yTeb5u4xjzn+f8pR9Duy3tmT14Ey7mZmZmVnFOdPeIP+Jt4xP2t3MeHc7m171rHm3n9xSBfl9zGfdG79RNt8WxroCkzdWdqeMJ2NULYvY7fpM5jdrdlO3M4z9nBGdrPft58xrVf4encSwKnWfLP3czurKmXYzMzMzq5onhxlua4VVrGLKAHdtFRG9rkPP7TN7elx779ZAZ594q/40lrFUPZveqGrZ9e4/p/3ZFq9/2ax4i20sKVSualmXKtxfMp5OYlaF+wW6cRWxaDtrZaxY9kNb6IVeHbOTeRV6ou3KJq6sdra3nr8vImYXLS9p+jSmrdmarQt3xJeznA3YgGWxTB1XtMIG9+OImZmZmfWliFi7CZuwghWFyq9jHWtZyxBDA9u3daYd+MPZW8TCew9Zb34ZGfR+y2RXTVUy672qRxmZ9tHbWzLRKvXVWNGq6GZGtN+uDvYqc9rPV0Srnm2uQmyrEqOyY1F0v/rhysVuWthWph3ay7YPepYdnGk3MzMzswoqmm2vQ5YdnGkHYL/ZO8XN934OcGa8iKpkv8tQ9X1plWlfv1xnmfc6aSfDVIXMYTf1w3nObbpaqtJm3C761zt1RtuZdiiWba9Dlh2caTczMzOzihov216XLDv4Oe0AiI2aZhGqnoW1Yvr579iq7o1Z98nMgk00nmNdMSi+jfazbYP6vQb5v0cZsZ1MY7Wlso/bTmLTq6xuFdpSo6q0rU5iU0ZbKnv/69a2JmqY4RnTmLZmJjPXy7a/yqvMZCbLYtnADx0Z+E8lZmZmZta/WmXb65RlB2faAZjCjL7Oxlp/Z9M7UUZGrFcxK/q+Y+1HPlvU7YxVGZmpyYr1oB4H3WzvRdvZ+ut1crWnf9pSVd63Kso4b43eXrnnsLIz6FX8ezfLttcpyw7utAMwhemVbKBmeUU7LmW05X46HhrrWvZl7LJjUcXYTubQhyrEs+zjpahuv1fV2lY/t6ui1h+q2H49Oh3iVrcPcBGxdpZmsYIVbMZmb2XZV7KyFll28PAYMzMzM+sDwwzPWM1q3uTNt7LsUaPHIDrTbjYAqpZhK0On+zSIsei2usWsblejeqUOMarClaM6Gcm2DzPMG7xRqyw7ONNuZmZmZn1imOEZa1hTuyw7+MuVAJD0KvBor+sxQLYGXup1JQaEY1kux7Ncjmd5HMtyOZ7lKiOeO0fENmVUpq48PCZ5tJNv6bLmJN3reJbDsSyX41kux7M8jmW5HM9yOZ7V4OExZmZmZmYV5067mZmZmVnFudOe/HuvKzBgHM/yOJblcjzL5XiWx7Esl+NZLsezAnwjqpmZmZlZxTnTbmZmZmZWce60m5mZmZlVXG067ZLmSXpE0gpJQ5Jul3RYk3I7SLohK/eSpMslbdyk3KmSfiNptaT7JB06OXvSe5I2k/RFSXdLGpb0fBazdzeU20VSNJmubrLN4yT9Kovnw5LmTN4e9U7RWGZlN5d0ZdZ+hyX9QNJWTcrVMpYjJM2RtEDS0qy9zW1S5pAWbfPiJmVre6xDsXhm5Xzu7ICkJS3a4oYN5QrFt84k7SXpvyWtlPScpAslTe11vapO0twWbfDTuTKS9A+Snpa0KutD7d/LetdRbTrtwEbA5cBHgE+QviRgkaSDRgpImgYsBnYG5gCnAyfScAOGpJOAecB3gSOBh4CbJO3T/d2ohJ2AU0mxOgH4FPAO4C5JOzYpfyZwcG46N79Q0p8A1wM/IcXzR8BVkg7v1g5USDuxvAY4BPgrYC5wALAwX6DmsRxxArALcFOBsqcwum3+a36hj3WgQDx97pywnzC6HR4MvD6ysGh860zSFsCtQADHARcCnwO+2Mt69ZkPMroNLsgt+3vgPOAS4BhgBXCrpO0mu5K1FhG1nICpwFPAZbl5JwPrgHfl5n0MeBPYPTfvUeDbud+nAL8Cvt/r/Zqk2G0CbNQwb0vSQXx+bt4upBPo0eNsbzFwW8O8m4Gf9XpfKxTLg7NYvj8378Bs3mGO5aj9nZL9nJnFZ26TModky/YZZ1u1PtbbiKfPnZ3Hdwkwf5wyheJb5wn4PDAEbJabdxawMj/PU9PYzc2O7Zktlm8IDANfyM3bBHgRuKjX9a/TVKdM+ygRsQ5YDszIzT4SuCciHs/NWwisAY4AkLQr8G7g2ty23gSuy9YfeBHxWkSsapi3DHgS2LadbUnaAPgzcvHMXA0cLGnzidS16tqI5ZHACxFxe67c3cDj2bLax3JEdjxOmI/1pGA8fe7srnHjaxwJLI6IV3LzriZdZf9Ab6o0MN4HbMboY/c14EZ87E6qWnXaszFZ0yRtJelvgd2Bb+eK7Ak8kl8nItYAj2XLyP0cVQ74NbClpG3Kr3n1Zfu9G/Bwk8VXSlqXjYn9mqSNcst+H5hO83hOIf2Tr5UWsVyvbWZ+zdtt0rFs321Z23xC0rkN4199rBfnc+fEHJ6Nw14pabGkfRuWF4lv3TWL0VOkTLtjVMxjkt6Q9KikT+Xm70m60vObhvL5/z82Cab1ugKTbA5wVfb6NWBOlq0csQUp+95oKFtG7mdjuaHc8hcnXtW+88+kIR35m0xfJ40R/jHwCmlIwtmkzuVxWZki8aybZrEcq23umitDk3J1jmUrw8DFwB2kbOXRpLGv25DGC4OP9Xb43Nm5nwL/AfyWNGb9HOAOSftFxBNZmSLxrTvHqHNLSePV7yYNHT4ZmCdp44j4Oil+K7IRCnlDwMaSZmQfIq3L+rbTnl3qf8d45SIi/8l7Menmva1JN6BdLemoiFiSX6XZ2zWZ3/i7xli/8jqM58i6f026uff4iHg5V3Yp8Jlc0SWSXgCukLR/RDyQ33TjZlvMr7xuxHJklWZv12T+wMQSJhbPMcreD9yfm3WrpNeBz0r6UkS8lC/eWKUW8/tCN+I5skqzt2syf6Di2ajd+EbE+bnZd0i6lZQxPiOb3lql2du1mF9XjlEHImIxqX80YlE23PJcSf8yUqzJqgN17PaDvu20k+6c/2aBciONiogYAu7Nfr1F0vakO8zfn80bAmY12cYs3v4EP5SbN9xQBpp/0u8HbccTQNKxwDeAsyPihgLrzweuAN4LPMDoeOb1czy7EcshUha4Uau22VgG+jOW0GE8OzCfdOPavsBt+FhvJ551Pnc2mlB8I+J5ST8nnSNHFIlv3bWK0eY4Rp2YT7rZeRdSbDeVNLUh2z4LWBkRa3tQv1rq2zHtEfGtiNB40zibuZ+3hxZAym6MGp8laUZW5pFcGRrLZb8vi4i+vLzbSTwlvY80hGNeRPxT0bdq+PkYsJbm8XwT+L/O9qh3uhTL9dpmJj+Oc+BiCaUd6229ZfbTx3pxtT13NioxvvnsZZH41l2zGO1IesqJY9S5IMVvKuleq7xW91pZl/Rtp32iJIn0GL383fiLgAMk7ZybdyywAXALQET8jtT5OTG3rSnZ74u6XO3KkLQ36bnNtwCntbHqCdnP+wAi4nXSM4pPbCg3B7gzIoYZcAVjuQjYLnsO+8h6s0n/tBeBY1mC44E3gAfBx3qbfO4siaTfA/6Y7ByZGTe+xiLgQ5I2zc2bA6wi3Tdg7Tme9H02TwL/Q7ovLX/sbkx6XruP3UnUz8NjCpP0p6QvWVhAejb7VsBfAAeRGt2I+aSbgBZIOo90We3rwA8jIn/X9AXA9yU9Afw829buwMe7uiMVIWlb0j+KFcBlwIHpMxAAr0TEw1m5C4BNSTF6hTQM6e+ABRHxYG6TXyKNd7+U9Bizo7Jp4B9lVjSWEXGnpMXAdyWdScqcX0J6/vqtuU3WNpYjJO0F7EV6tjDAbEkrgBcj4qdZmX8j3fR4D+lG1KNI919c2nAvwQXU+FiHYvHE586OZE+J+QrpsZdPkr5s7fOk4/vSXNGi8a2zeaSkxwJJl5ASGhcAX2t4DKQ1kHQ96SbUB0kZ9TnZdFr2WNbVSt8WfZ6kIVJ2/bOkxO83elPrmooKPCy+2xNpTNZ84BnSE02eIWU2D25S9p2kzs4K4GXS0082blLuVNLd/q8DvwAO7fV+TmI8DyFdMms2LcmVO4l0D8EwqWP0W9I9BBs02eaHgf/N4vkIcFKv97NKsczKzgKuJI3PfAX4IbC1Y7ne/l9QoG2eRvoH9WoWp4dIN/1NabK92h7rReOZlfO5s/3Y7kD68rOl2TnyZdI3Gu/ZpGyh+NZ5In24vI2UXV9KSmJM7XW9qj4B/0j64rOVWezuAz7ZUEakD47PZGXuAN7T67rXbVL2xzAzMzMzs4qq7Zh2MzMzM7N+4U67mZmZmVnFudNuZmZmZlZx7rSbmZmZmVWcO+1mZmZmZhXnTruZmZmZWcW5025mZmZmVnHutJuZmZmZVZw77WZmPSTpIUmHlLSt70i6qMN1Q9Jrkr48wTrcJmm1pJ9NZDtmZjaaO+1mZhlJT0g6bDLfMyL2joglvXr/BvtFxDkjv0j6I0mPSHpY0uzc/CckrZK0IjddDhARHwQ+3YO6m5kNtGm9roCZmVXWl4Gjs9dXAIfnlh0TEbdOfpXMzOrJmXYzs3FI+gNJSyQtz4azHNuw/AlJZ0p6UNKwpGskbZgte6+k+yW9Kum6bNlFDeseJul7wE7AjVnm+qxseUjaLVf+rSEwkt4j6RfZtq8BNmyo1/aSrpf0oqTHJZ3W5q5PzaZpOMljZtZT7rSbmY1B0nTgRuDHwLbA3wA/kLRHQ9GPAUcA7wL2BeZKmgHcAHwH2BK4CvhIs/eJiE8CT5Ey2DMj4qvj1GsGsBD4Xrbt64Djc8unZPX+JbADcChwhqQPFd134HzgR8B/AueMU9bMzLrInXYzs7EdBMwELo6INRFxG3ATcHJDucsi4rmIWEbqLO+frTstW7Y2IhYAd5dYr+nApdm25wP35JYfAGwTERdm9f4d8E3gpKJvEBG3R8RuEbFHRNzZsHhhduVhZDp1ojtkZmat+XKnmdnYtgeejog3c/OeJGWv857PvV6Zrbc98GxERG7Z0yXWq3HbT+Ze7wxsL2l5bt5U4I6S3v/DHtNuZjZ5nGk3Mxvbc8CO2XCTETsBzxZYdymwgyTl5u04RvloMm8lsHHu9+3G2PZOuddPA49HxKzctGlEHFWg3mZmVjHutJuZjTZd0oYjE3AX8BpwlqTp2TPVjwGuLrCtO4F1wGckTZN0HHDgGOVfAHZtmPcA8HFJUyUdAXwgt+03gNOybX+0Ydt3A69IOlvSRtn6+0g6oEC9zcysYtxpNzMb7WZgVW76AnAscCTwEunRh38eEY+Mt6GIWAN8FPhLYDnwCdJ4+NdbrPIV4NxsjPiZ2bzTSR8SlgOnkG4+zW97LjAEzAEW5N57Xbbe/sDjWd2/BWw+fggKGXnKzch0Q0nbNTOzJjR6OKSZmXWTpLuAeRFxZa/rkidpNenDxGURcd4EtvNfpJtk746IQ8uqn5lZ3bnTbmbWRZI+ADxKynSfAswDdo2IpT2tmJmZ9RU/PcbMrLv2AK4lPTbyMeAEd9jNzKxdzrSbmZmZmVWcb0Q1MzMzM6s4d9rNzMzMzCrOnXYzMzMzs4pzp93MzMzMrOLcaTczMzMzqzh32s3MzMzMKs6ddjMzMzOzivt/eREXZFsC7E8AAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 864x396 with 2 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"plt.rcParams.update({'font.size': 15})\n",
"fig, ax = plt.subplots(nrows=1, ncols=1, figsize=(12,5.5))\n",
"dummy = np.ma.masked_invalid((psi_ridge2.mean('time')/1.0e6).values)\n",
"xyplot(dummy, grd_ridge.geolon, grd_ridge.geolat, grd_ridge.area_t,axis=ax, clim=(-150, 150),title='using uhbt')"
]
}
],
"metadata": {
"kernelspec": {
"display_name": "analysis",
"language": "python",
"name": "analysis"
},
"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.3"
}
},
"nbformat": 4,
"nbformat_minor": 4
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment