Skip to content

Instantly share code, notes, and snippets.

@ocefpaf
Created July 12, 2022 20:57
Show Gist options
  • Save ocefpaf/07b13c841297e9415b07a4f3b82721e0 to your computer and use it in GitHub Desktop.
Save ocefpaf/07b13c841297e9415b07a4f3b82721e0 to your computer and use it in GitHub Desktop.
pyugrid_uxarray_comparison
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"metadata": {},
"id": "ee65a365",
"cell_type": "markdown",
"source": "# pyugrid"
},
{
"metadata": {
"trusted": true
},
"id": "fb13a11b",
"cell_type": "code",
"source": "import pyugrid\n\n\nurl = \"http://www.smast.umassd.edu:8080/thredds/dodsC/FVCOM/NECOFS/Forecasts/NECOFS_GOM3_FORECAST.nc\"\n\nugrid = pyugrid.UGrid.from_ncfile(url)\n\nlon = ugrid.nodes[:, 0]\nlat = ugrid.nodes[:, 1]\ntriangles = ugrid.faces[:]",
"execution_count": 1,
"outputs": []
},
{
"metadata": {
"trusted": true
},
"id": "7e54f590",
"cell_type": "code",
"source": "print(ugrid.edges)\nprint(ugrid.edge_coordinates)",
"execution_count": 2,
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": "None\nNone\n"
}
]
},
{
"metadata": {
"trusted": true
},
"id": "8f48d25d",
"cell_type": "code",
"source": "ugrid.build_edges()",
"execution_count": 3,
"outputs": []
},
{
"metadata": {
"trusted": true
},
"id": "6db3f874",
"cell_type": "code",
"source": "print(ugrid.edges)\nprint(ugrid.edge_coordinates)",
"execution_count": 4,
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": "[[ 3014 3015]\n [16623 16624]\n [ 8282 9039]\n ...\n [26778 27095]\n [39666 40190]\n [25310 25600]]\nNone\n"
}
]
},
{
"metadata": {
"trusted": true
},
"id": "d9971562",
"cell_type": "code",
"source": "import matplotlib.pyplot as plt\nimport cartopy.crs as ccrs\nfrom cartopy.mpl.gridliner import LATITUDE_FORMATTER, LONGITUDE_FORMATTER\n\n\ndef make_map(projection=ccrs.PlateCarree()):\n fig, ax = plt.subplots(figsize=(8, 6), subplot_kw=dict(projection=projection))\n ax.coastlines(resolution=\"50m\")\n gl = ax.gridlines(draw_labels=True)\n gl.top_labels = gl.right_labels = False\n gl.xformatter = LONGITUDE_FORMATTER\n gl.yformatter = LATITUDE_FORMATTER\n return fig, ax\n\n\nfig, ax = make_map()\n\nkw = {\"marker\": \".\", \"linestyle\": \"-\", \"alpha\": 0.25, \"color\": \"darkgray\"}\nax.triplot(lon, lat, triangles, **kw)\nax.coastlines();",
"execution_count": 5,
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAfQAAAEfCAYAAAC3VTTLAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAAsTAAALEwEAmpwYAAEAAElEQVR4nOy9eZicdZU2fD+171Vd1Uu6sxEIQRZBEBgXUFBUEFBfdRj0c/R1GddBvoFRR0XF+dAXHcZteB0dnXFcxlGHGQbHBUEgRECIBJCEnexJJ73UvnWtv++Pzn361NPVSSfpTjrJc64rV7qrn3rq2ep3zrnPfe5jGWPgmGOOOeaYY44d2eY63AfgmGOOOeaYY44dvDkO3THHHHPMMceOAnMcumOOOeaYY44dBeY4dMccc8wxxxw7Csxx6I455phjjjl2FJjj0B1zzDHHHHPsKDDPbDe0LMsN4GEAO40xl+157SoAfwmgCeCXxpiPd3nfFgBFAC0ATWPM2XteHwLwoz1/+3+MMSXLsq4H8HEAxxljRvdsVzLGRLod08UXX2zGx8dnewrzavV6HT6f73AfxmG18fFxbN26dZ/bBQIB9PT0IJVKwe/3H4IjO3R2rD8H833+27dvRzabRaPRAAAMDAxgyZIlHdvUajVs2LABS5cuhdfrhdvtht/vl2dtYmICfr8fO3bswNjYGHw+H0466SR4vd45OcYj5RnIZrPYtGkThoaGMDg4OKv37NixAyMjI1i0aBEWL17cdZsj5fzn0+b7Gqxbt+43xpiLp/3BGDOrfwCuAfBjAL/Y8/uFAH4LwL/n9/4Z3rcFQG+X128EcCqAywF8cM9r1wPYBuBLarvSTMf04he/2CwUu+eeew73IRxya7fb5uGHHzbXXHON8fv95qabbjIATDKZNKeccorp6+szbrfbDA4Omj/90z81t99+u6nVaof7sOfVjsXnQNuhOv9Wq2VWrVpl3vzmN0/7W7vdNmeddZYBYHw+n/F6vQaAWbRokXnVq15lAJidO3caY4z59Kc/bQAYAGbdunVzcmxHwjNwySWXGADmrLPOMq1Wa9bv+7u/+zsDwPz85z+fcZsj4fzn2+b7GgB42HTxibOC3C3LWgLgUgDfVS9/CMCNxpjansBgdDb7UuYG0N7zz1Kv/wuAP7MsK7mf+3PsENkf//hHnHbaaXC5XDj77LPxla98BbVaDT6fD4888gjS6TSeeOIJjI6OotlsYnh4GD/72c/wute97piP3B2bG3O5XHC5XMhkMtP+ZlkW1q1bh3a7jVqthnq9js2bN+PSSy/Fgw8+CAD47W9/CwD4+Mc/jqVLlwIANm3ahGaziQ996ENYuXIlzj//fHzsYx/DunXrDt2JHSLbuHEjAOC+++6DyzX7yusPfvADRKNRHH/88fN1aI4dhM32Tn4Nk1B4W722CsD5lmU9ZFnWvZZlnTPDew2AOyzLWmdZ1vvV6zcD+DaAD2ISeqeVMOnUr57lsTl2iO1Nb3oTnnjiCQDAkiVL8Fd/9Vf4zW9+gxe+8IU488wzD/PROXas2Nvf/nasWbMGpVKp698taypPOO644/Dd734Xxx9/PHp7e3HRRRcBAH79619j+/btAIBms4kvfvGL+Na3voXly5ej0Wjg+9//Ps4++2z09vZi+fLluPjii/Gd73wH7Xa762ceCfab3/wGzz//PL72ta8hGAzO+n2//OUv8cQTT+C3v/0tTj311Hk8QscO1CyzD+lXy7IuA/B6Y8yHLcu6AMBfG2MusyxrA4C7Mel4zwHwUwDHG9sOLcsaMsYMW5bVD+BOAFcZY9bM8FnXY9KhfxfAYwBOBzBsZqihn3766eYb3/jGLE91fq1UKiES6XqYR41NTExgx44dyOfzAIATTzwRsVhM/n4sXIN92bF+DQ7l+Y+OjmL79u17refabf369YjFYli+fDkAdGTflmXB5/OhVqvhRS96EdxuN4BObkgoFEKlUkE4HMayZcsQCoWmfYb9GtTrdTQaDezevRvNZhPNZhPGGLRaLXsZEgA6fna73YhEIojH40ilUvuVTXezVquFDRs2wOfz4eSTT96v9z7xxBPweDw46aST9rrdsf4dAOb/Glx44YXrzB4+mrbZkOJeDuANlmW9HkAAQMyyrB8B2AHgv/Y48LWWZbUB9AIY0282xgzv+X/UsqxbAZwLoKtDV+/JWZb1YwAf3tt2Pp8PF1xwwSxOYf5t9erVC+ZY5tLy+Ty+/vWv4+abb8bY2NStve666/CGN7yhY9uj9Rrsjx3r1+BQnn86ncYHP/hB/PznP8cPf/hDXHHFFXvdftu2bbjwwgtx0kkn4UMf+hCGhobwwAMPYGBgAD6fD9/73vcwMTGBer2O5557DitXrpT3PvLII8hms3j1q1+Nu+66C5dddhlcLhdqtRpSqRSSySQsy0K9XsfVV1+Nz3zmM6hUKkLes9vFF1+MM888E6lUColEAuFwGK1WC61WC8cddxxWrlyJp59+Gv/6r/+KH/7whwCA3//+93jJS15ywNfLGIOXv/zleOyxx3Dbbbft1336/e9/j6uvvhpr1qzBy172sr1ue6x/B4DDdw326dCNMZ8E8EkAUBn6OyzL+iCAVwFYbVnWKgA+AB2Uc8uywgBcxpjinp9fC+BvZ3lsXwHwh9kco2Nzb8YYfOlLX8InP/lJeS0cDuNlL3sZrrnmGlx88XSCpWOOHUpLpVK46aabcMcdd+DOO+9ET08PXvOa18y4fW9vLy699FI89dRTuO6661Cv15FIJDAxMSGOHACuueaaDmcOAGeddZb8/OpXvxqbN2/GZZddhre85S3I5XLYtm0bACAYDCIej+OTn/wkTjrpJDz44IP4+Mc/jmQyiYmJCQwPD2Pr1q04//zz4fFMLW1PPfUUKpUKPvOZz+BHP/oRHnnkEXziE5/Ao48+ikWLFuGnP/3pQTlzAHjXu96Fhx9+GHfddRfOP//8Wb/PGIP3vOc9WLVqFV760pce1DE4Ns/WjSk30z8AF2CK5e7DZO17A4BHALxqz+tDAH615+fjAfxxz78nAHx6H/u/HpMBA3//yuQhOiz3Q2WtVsv87Gc/M9FoVNi/V1xxhbCC92ZHyzU4GDvWr8GhPP9nn31WnlH+u+uuuw5oX+12W/axePHiA97H6tWrzU9+8pP9Yo4PDw93nINlWcblcplzzjnH3HfffQd0LHb71a9+ZSzLMr/61a/2+72/+MUvjNfrNVdccYX5zGc+Y0ZGRsw//dM/mW9961tdtz/WvwPGHD6W+35lv8aY1QBW7/m5DuAdXbYZBvD6PT9vAnDGfuz/etvv12CyXc6xebK1a9filltuwdq1azEyMoKdO3eiWCwCAE455RSsW7cOgUDgMB+lY45Nt0AggAsuuACLFy/Gtm3b8Lvf/U6e3f01y7Jw/fXX44YbboDH40Eul0Mikdjn+9atW4cVK1bgJz/5CT7ykY8AAG666aZZfeZ//ud/4je/+Q2+853vAADOPfdcPPvss3jve9+LT33qU0gm56bR549//CPe+MY34sorr8Qll1yyX+/dvXs3LrvsMgCT5/qzn/0MX/rSl9But9FsNuH1evGe97xnTo7TsYM3B84+Bm3r1q349Kc/jdtuuw3lchm9vb044YQTcPLJJ+Ppp5+Gy+XCV7/6VXz0ox893IfqmGMz2tKlS3HPPfcAAL7//e/jd7/7HT72sY/huuuuQzgcxnHHHYdmswmXy4VTTjkF11xzTQeJ026f+9zn8N73vhdnn322dG+88pWvxPHHH498Pj+tg2Pz5s04++xOXtKf/dmf4cwzz9wree0///M/8ZOf/ARPP/00NmzYAGAqOLntttuwaNGiA70k0+xf//Vf8f73vx/nnHMO/u3f/m2/33/hhRfKz2x1u/TSS3Hqqadi7dq1eN/73odbbrkFzzzzjPztzW9+89wcvGP7b93S9iPlnwO577+1Wi0DwPT19ZlPf/rTplAomHq9bv75n//ZpFIpEwqFzNjY2AHt+0i5BvNpx/o1OFzn/9xzz5lLL73UvPOd7zRvfetbzXnnnWdWrVplTj31VHPiiSeaaDRqLMsy5557rvngBz9o3vnOd5oPf/jD5jvf+Y6pVqsd+yqVSubd7363icfjxrIs43a7jd/vN6eccor527/9W7NlyxbzgQ98YBrc7/P5TE9Pj7npppvMZz/72a7HWavVpr3nne98p2k2m3N6PbZs2WJe8IIXGJfLZd797nebiYmJ/d7HxMREx7H29fV1/P7Sl75Ufl6+fLn8fNNNN5lly5aZb33rW2b16tXm+eefN6VSydTr9Tk9x4VshwtyP+xO+WD+OQ59/+173/uecblc5m//9m/N//7f/9uceuqpxuVyGZ/PZy666CKTy+UOeN9HyjWYTzvWr8FCPv9LL71UnM6KFSvM4OCg8fl8xuVymT/5kz8xX/va10ypVOp4T6VSMWvXrjVf/vKXzcte9jITi8UMAON2u2Vfg4ODZmxszPziF78wZ555pigmDg0NmUql0rG/QqHQ4RSvuuqqOT/Pu+++2wQCAbNixQqzefPmA9rHyMiI+djHPjYtaAFgvv71r5v7779ffj/hhBPMk08+abZt22YeeOABc8sttxi/329CoZBZvHixecUrXmEef/xx8+yzz5pyuTy3J7tAzXHojkOfN5uYmDD/8R//YV73utfJl9DtdpulS5eaSy65xPz85z837Xb7oD9nIV+DQ2XH+jU4VOf/xz/+0fzHf/zHAT23OlttNBrmL/7iL+R7cfXVV+/z/eVy2dx8883muuuuM9u2bev4265du8zNN98sGW0+n++6j3g8Lp85107uVa96lUkkEvsVnG/fvt3s2rXL3HLLLeb00083Ho9nmiMPBALm/vvvN8ZMIn3j4+Nm3bp10yRz77nnHpPNZs2nPvUpc+mll5oPf/jDZs2aNebee+81GzZsOCac+hFBinPsyLNWqzWN1PahD30IN954417riY45tpDtl7/8JT71qU/h3HPPxUMPPTTr9xlj8P73vx/r1q1DLpfD2NgYLMvCC1/4Qrz97W/HBz7wgX3uIxQKCQFO2+9+9zu84hWvEFLc2Wef3fU7lslk4PF4cMYZZyCdTs/ZUJhKpYLPf/7zWL16Nb75zW8iHo/v8z2tVgvveMc78JOf/KTj9SuvvLLjtfe+97244YYbpL7vcrmQSqWQSqW67vexxx7DH/7wByxZsgQveclLUCgU0Gw2ZWjJcccdN2fn7diUOQ79KLdvfvOb8vPHPvYxfPaznz3mVZwcO/Ltk5/8JJ5//nn8y7/8C8rlMsLh8Kzet337dvzgBz+Q36+44gp84AMfwIUXXtghFXsgdt555+HXv/41yuUygElZ2SeffBKnnHJKx3b33Xcf0uk00uk0AMyZY/vxj3+M7373u/D5fHj/+9+/122NMbjlllvw7ne/W443GAyiWq1i2bJlsCwLbrcb8Xgct99+O8455xxks1l873vfw5ve9Cb09PTMuO92u423vvWtuPLKK7FixQo0m02sX78ezWYT4XAYiUQC5XJ5Vl0Eju2fOfPQj2LbuHGjMNU3bdqEL3/5y44zd+yoMT7bkUgEIyMjs3rPsmXLsHbtWlx++eU488wz8V//9V947Wtfi4985COTNciDMMuycPHFFyOVSmHbtm0IhUI47bTTEAqFEIlEEAqF4Pf7O1jgX/jCFw7qM7W95S1vAQAMDQ3NGJw0Gg28/e1vh8vlwhVXXIFyuSx67tdcM9khPDo6in//93/H29/+dvyf//N/8OUvfxnHHXcckskk3vOe9yCXy+31OFwuF77yla+gp6cHLpcL1WoVqVQKixcvxsqVK1GpVERCeiYlPccOzJwM/Si1QqEgqk5//ud/jhUrVhzmI3LMsbm1M844A4888gjOOussnHzyyfjhD3+I448/HkNDQwiFQsjlcvB6vQiHwx1Z8DnnnIOf//znACYHstxwww34/Oc/j//6r//CrbfeOidqaEuXLkW5XMbGjRtx2223YXh4GKeddhrC4TD6+vowODiIVatW7Rcq0Gg00Gw2Zxyo0tPTg4GBAcm4u9mjjz6Kf//3f5fffT4fGo0G+vr68KUvfQnA5MwGALj99tvx05/+VBT0gsEgHnjgga5rSaPRQKlUQiAQQDabxb333ouBgQHRrU8kEojH47AsC36/H5lMBvV6Hb29vRgcHHTg9zkyx6EfhbZ27VpccsklMlryhhtuOMxH5Jhj82MvetGLcO655+K1r30tdu3ahVwuh40bN8Lv92PNmjUwxuC8887DS1/60q41ZY/Hg+uvv15Gpr7sZS9DIpFANBpFo9FANptFq9VCs9kEMKkHz3Grs7ETTjhBMt8DtYmJCbz73e+Wmnar1Zqxz/1zn/sc3va2t6HRaHR1kueeey6+853v4C/+4i8AQJz1Jz/5Sbzvfe+D3+/H9ddfj1tvvRUvfOEL8eCDD2J4eBgf//jH8YUvfKFrAJJOp/GP//iP2LFjB0ZHR3HppZciFArhhBNOQCqVgs/nQ6vVgmVZqFarACYh/56eHlSrVeTzecTjccepz4E5Dv0oszvvvBOve93rYIyBZVn4+c9/jmXLlh3uw3LMsXmxfD6Pj3zkI1LztSwLhUIBxhjJgHfs2IHHH38cL37xi7tORwOAgYEBFAoFuFwu5HI55HI5vP71r8c73/lOfP3rX8fvf/97AMD555+PT33qU2g2m8hkMrj22mvh9Xpx8803Y8WKFVixYsVeR6u22228613vwo9+9CO8+tWvlrnsM9m///u/4wMf+ACMMTj++OPRbrf3KlrT29sLj8eDxx57DOec032itdvtRjQalVnxwOQo2mg0CmCyDPCqV70Kb3nLW+ByufDAAw/g3HPP7bqvSqWCO++8E9u3b0ez2cSZZ56Jnp4evOENb4DH40G1WoXH44FlWTLNzuVywbIsVCoVFAoFVKtVVCoVJ1OfA3Mc+lFkX/3qVyUbSCQSeOCBB/Z7RKJjjh0p1mg0cOeddwKYzLTdbrdkgpZlweVywRgDv9+PUqmEfD4/o0MHJmvgxkyONX344Ydx7rnnwrIsPProo1i7di1arRa2bt3awYT//ve/j4GBAdx///0AJgcYff7zn8db3vIW9Pf34/jjj8fy5ctxzz33IBQK4cc//jH+8Ic/AJgc8jKT3XbbbfjIRz6CXbt24YorrsAPfvAD/PCHP8QHP/hBtFotGe1qt1e+8pWIxWL4wQ9+0NWhF4tF/OQnP0G5XMYLXvACvO1tb8PKlSsxMjKC/v5+tNttXHTRRXjwwQdx8skn4+677+5KXms0GqhUKnjuuefw7LPPYsmSJTDG4KSTToLX60W9XhckwePxoNVqodFowO12o16vw+/3o9lswu12w+12I51OIxaLOUS5gzTHoR+BNjo6iv/+7//Go48+iueffx47d+7Ejh07RMf6hBNOwNq1a+dMC9oxxxairV27FtlsFpFIRJy5NpLcdHY4G3O73fiTP/kT+f3GG2/EjTfeCAC45JJLcPvtt8vfnn/+eTz//PM477zzMDAwgPe///0ol8u49tprsWHDBjzzzDNYvXq11KXf+973Yv369V0z0VKphO985zv4yle+gp07d+KSSy7BunXrMDAwAGByWtqHPvQh3Hzzzbj66qu7HrvH48GZZ56JBx98sOvfv/zlL+OOO+7AWWedhRe84AW44YYbUK/XYYyBy+VCMBhEuVzGTTfdhGuuuWYaxM5a+datW/Hss8/ivvvuQ6PRwEtf+lL09PTA4/Gg0WggEAgIj6Fer0uWTuSQ29XrdRQKBQQCAezcuROWZSEUCjmZ+gGa49CPIDPG4NOf/rQMR6DpL0pPTw/uuusux5k7dtQZIdpWq4XNmzdjw4YNiEaj05yOdhrAJJlr6dKlc6K78Otf/xqZTAa33norLrvsMpTLZSxZsqQjWFi9ejWuu+46+f2b3/wmrr/+eoyNjeH3v/89XvziF+NP//RP8Za3vAW1Wq1jNKvP58PFF1+Mhx56CENDQx2f7Xa7sWrVKlxzzTW477778M///M9yTtVqFR/5yEfwm9/8BqVSCaVSCffcc49osdN5jo6OAgAef/xx7Ny5E1/4whfw3ve+Fx6PB2eddRZ27NgBYHJeu2VZ8j5jDMrlMnbt2oUHH3wQTzzxBJ588kmcddZZuPDCC+H1euWas9+80Wig3W7D7XZLwOB2uyVjtywL0WgUhUIBoVAI2WwWtVrNIcodhDkO/Qixer2Oc845B48//jg8Hg9e85rX4De/+Q2AyalozzzzDJrNJr72ta9h+fLlh/loHXNs7qxSqWDr1q1YvXo1Go0GIpEIPB4PYrFYx0xxmq6nDw4OYuXKlUilUnPmIJLJJN773vfOatsdO3ZME6EZHh7GDTfcgM9+9rMdrz/xxBM4+eST98p8X7t2Lf7t3/4N/+//+/8ikUggHA6j1WqhWq3C7/fjyiuvxPj4OH75y1/ioosuwre//W38r//1v7B161YYY/BXf/VX+Ou//mtks1m0223s3r0b3/72t/HUU09h8eLFGBwcxFVXXYXly5cjn8/j+eefx9q1a7F161YMDw8jFAqhWq3irLPOwmte8xr4/X60Wi15HcC0kgfvRywWQz6fFw5AKBRCpVJBq9VCqVSCy+VCIBBwiHIHYY5DP4y2L4ILrVQqYdmyZchms4hGo/j4xz+Oz3/+8wCAv/7rv8by5ctx1VVXIRgM4sorr5zvw3bMsXm3RqOBcrmMarWKe+65B/l8HuFwuAO6pTMn1K5fi0ajOOWUUzA0NHRYncLQ0BBuv/12/OVf/iVuuukmVCoVvO1tb8PDDz+Ml770pTjxxBPx61//etZBeDAYxPve9z68613vwq9+9Sv87Gc/w+9//3ts3rwZmzdvxuDgIADgv//7v/H2t78d3/ve97Bz507E43E88cQTWLVqFZrNJtLpNOr1OkKhEEZHRwUiP+6443DHHXdg3bp1eP755xGPx5HNZnH22WcjkUhgcHAQAwMD8Hq9qNVq8Hg8Qqzzer0dMqQulwvNZlOy9kqlAmOMoAoejwculws9PT3yt2q1iomJCYcod4DmOPTDaCS2fOhDH8LAwACSySTe9a53yQOfyWTw93//9/jqV78qNbiJiQl85jOfAQB84hOfwI033ohEIgHLsvCtb33LefgdO+S2ceNGjI2Noa+v76D31Wg0MD4+jt///vfYuHEjarUaBgcHEYlE4Pf7O7ZtNpuSzfK75PF4MDQ0hJNPPnlW0qfzbS6XC6973evw3HPPdbx+9tlnH5SoitfrxRvf+Ea88Y1v7Pr3N73pTfja176GH//4x+LQTz/9dCGrJRIJGGNQr9eRSqUQDofxute9DsViEevXr0ckEsEpp5yCM888E+FwGBMTE4jH4wK/AxDYnA6b7HWXy9Xxj73ufr8flmWh2WxKP321WkWz2USj0UA4HEaz2UQgEEAwGMTExISTqe+nOQ79MNozzzyDk046Cf/4j/8or330ox9FNBpFMpnEyMiIsENZf4rH43jFK16BL37xizjppJNQr9eRz+fhdrtxxRVXHMazcexYs61bt3boqG/btg1erxc///nP8YY3vGHW+6nX6ygWi3jooYdw1113AZhUdIvH4wiFQtPIbNqRAxB4PZFI4JRTTkF/f/8x7wAajQb+7M/+DK1WC5lMBsFgUORa8/k8arUa2u02ms0m8vk8zj33XOkUeNWrXoV6vY5gMIhms4lQKASXyyUsdYrDMDP3+XwIBAIIBALIZDJSNw+FQoKaNBoNyejb7bZsXyqV0Gq1YIwRuJ0Q/MTEBCqVipOp74c5Dv0w2qpVq5DNZvHggw/im9/8JtauXYuRkREUi0VhrPf19WFsbAztdhterxfj4+O49dZb8etf/xq9vb0YGRmB1+tFq9WalsE45th82YYNG/Dkk092vEZ4dWJiYlZOndn4mjVr8PjjjyMcDmPp0qXo6emZUQ2tG7weCoUwMDAgCmQceEInxNY1vp/vDYfDUgOmEIsmdfl8vo4avcfjOSIY2JVKRdCNpUuXSg2blkql0Gw2sXHjRrTbbZx11lmSObtcLrRarY5Ws3w+D5/PJyWPdDot6AIzaZLnmHhwfwAkEeH1q9fraDabssaxtl4ulzsCCWbqtVrNydRnaY5DP8yWSCRw8cUX4+KLLwYwucj99re/xc0334xf/epXGBsbg9vtxvLly/HRj34Uw8PD+NrXvoZms4kdO3YI/HXZZZcd9HAJxxybjaXT6WnOnPArjSWibtZoNDA6Oor7778fTz75JFKpFE444YR9zhmgg9BOtlarYXh4GLt27RKHDEAIWbRutV32RFMVrtFoSI93u92W95OEFwqF0Nvbi1WrVi0IOL+bNRoNPPbYY1i3bh1GR0cRDAaRSqUQDAblmvAacGhMMBiUhKDdbqPdbkvvPoMBYwxqtRqq1aoEP81mE+VyGa1WCz09PSgWi3C5XGi32wgGg9Lz7/F4xBHzudCEuXq9Ls7b6/UiGAxK5s5Mndm6k6nv3RyHvsDM6/XikksuwUte8hKcd955ePLJJ9FqtbBp0yZcd911mJiYwHve8x5cfPHFWLlyJY4//niUSiX09/cf7kN37BgxwuIApvV+780ajQYymQwefPBBrFu3DqlUCkuWLJlVO5n+HCIBtVoNlmWhv79fggk6CQB7DXCZsbrdboHwZ+pTb7VaaLVamJiYwLZt27B+/XpceeWVexWpOVxWr9fxuc99DmeddRZ6e3sFmahUKlKv9ng8CAaD4ngZvNRqtWnIBp2/y+WC1+uF1+vtKAGGQiEUi0WUSiU0m00kk0lUKhVYloWJiQnU63XE43EYYwSKp+BMIBCQTJ5CM5SHnZiY6MjUKR9br9cdh74Xcxz6ArL7778fl112WddpRh6PB5dffjn+6Z/+aVomM9vRkY45Npdmd+aEWmmEYWu1GtavX4+77roLExMTGBoawpIlS/ZLFYyOl8YMG+jM2LuJy3RrbbOfBx3N3iauWZaFcDgsTmshOnSfz4e/+qu/wvr164WBblmWTFWbmJhAJBKRfvRyuYyXv/zl8Pv94tAp1+rz+WS4CqF4Zu+8Ztw2Eokgm82iXC7LNgy6CK2zXNFut1EqleSaRyIRIc7pa8rjbbfb8hzFYjH4fD7Hqc9gjkNfALZ27Vpce+21uO+++wBM1tZ3796NQqEAADj11FOxdu3aBbmAOLawLZ1OY3R0FE8//bTUPV/0ohdh1apVB7S/SqUiP9vV2Zj90X77299i9erV6OvrQ7PZxMDAAMLh8AE/x90ct10CVWfc+zJdj9fBAoBppDuem8vlklGoC9G8Xi8uuOACnHjiicIFePbZZ1GpVJDP55HP56WEMDExIa1mzJ6BKYU9Y4xMqnO73QiHw2g0Gh0a7VyjqtUqvF4vAoEAisWiBHNerxehUAjtdhuVSgW1Wk361gOBgIjPcJAUMDk1TuvM83gLhQKGh4eRz+exdOlSx6l3McehH0b74Q9/iI9//OMYGRnBwMAAVq5cic2bN+PZZ58FMJl5X3fddfjEJz7h1Mcd228bHx/Hvffe2zEtjIM7KEQUCoXwghe8YMY+aLKM6cDuvvvuGZ0f/2db0tjYGJYvX45wOAyfzzej/viBGPelAwqdXc/k2PkaoXktfGI/J7fbLTVh/pxIJHDqqacuSIdOfXWOK3388cfxL//yL0ilUsjn8+jt7QUAPPXUU1LuSKVS6O3tldo5iW3MzMmE5zMUCoWkzWxiYkLU3phB12o1GGMQCoWE2T6TEydrPhqNCgTfbreFr2CMkRo6YX6fz4dSqYRcLodEIuE4dZs5Dv0QW7vdxmc/+1l861vfQjqdltdHR0fhdrvxqU99CpdeeilOOukkmR/smGMHYo8//rg4PD5HhEoJg1arVTz00EN49tlncfLJJyOZTCIUCiGfz2Pjxo3Yvn07arUafD4f+vv7pT5KoxOlQ6RKmMfjmZO+9G5mDyj053bL2GmauGeH4YPBIPx+P4aGhhAMBuHxeBAIBIT1zgy0p6dnwTjzRqOBH//4x/ibv/kbXHXVVTj++OOlnp3P55HL5XDcccfhjDPOQDQaRTabRT6fRzAYRDgcltnpxWJR6tTMxEulEhKJhPAU2IvOtrN4PI5CoSBEOdbBvV4v8vm8EOrC4TACgQAmJiZk8I39HpGAx/LA2NiY1PHpxNvtNmq1GgAIrF+pVA67cNBCszl16JZluQE8DGCnMeYyy7KuB/AXAMb2bPIpY8yv9mz7dwAuBHCtMeZey7KOA7AZwEeNMf+wZ5ubATxsjPnXuTzOw2WtVguXX3457rzzzg5pyi996Ut4+9vfPivVOMccm60RsrSbPXOlctgDDzwAYNLZEVrXoy63bNkCoBOO7laf3t9MXLOv92Z2R27/bNZyXS4XhoeHO+rudgEaAIjH44jFYujv70c4HEYsFlswznpvVq/XMTw8jPHxcezatQtvfetbkc1msXHjRsmABwYGsGjRIiSTSWGvh8Nh6UVnoEIN9mAwKEptlmUJEZCZOp2ty+USFEAPdCmVSgKTu1wuRKNRec5IZmMGHo1GUS6XpctgfHxckJBoNIqJiQkhwfGYyAHQjHi2uTkOfcrmOkO/GsBTADRt9avGmJv0RpZlvWDPj68A8K8A7t3z+yiAqy3L+rYxpvtqdITa97//fXzsYx/D2NiYvPbZz34Wn/vc55ws3LF5sSVLlkh7ma5360yWWSuzJ9aT6cj1SNKZ4OmDtX05c6qT0aFoRx6NRhEMBkVpMRAIwOv14uGHH8bGjRtFnYz7bzQaCAaDOO200zAwMHBE9JWzHt1sNrFlyxb8/d//fcexr1y5En6/H/39/R1Sq8FgUNjqzHZbrZb0hFPgJZPJoFgsyqSzQCAgjpP180KhIOx4BorBYFCmqbndbkQiEQkIqOvOunqr1RInDUxm2XzmIpFIh7Nmrd3v9yMcDktmTtZ7vV7vIPsBcIhye2zOHLplWUsAXArgCwCu2cfmbgBtAAaA9mZjAO4H8C4A35mrYzuc9stf/hKXX355x4L1spe9DF/96ldx7rnnHsYjc+xot9NOOw3AZM3UGDONxKZJYfY6uHbafN++2OIHajM583q9jna7LT3SbJvy+XyIRCJYvHgxent7uzrlM844AwBkghcA5HI55PN5AMCb3/zmeTmXuTLWw8vlMm6//XY89NBD0lP+kpe8RNTZdKsdZ0NQ7x5ARwDk8XjEubfbbViWhVKpBK/Xi0QigWazKdebLPdarYZSqQRjDCYmJgSC1zVwt9sNr9eLUqkEYPIZojY8IXbqs7tcLoTDYYTDYUFQ2u02yuUy/H6/PAtEBjiSldl/u91GJBIRmdjR0VEYYxCPx50edcxthv41AB8HELW9/peWZb0Tk1D8tcaYrDHmCcuyQgDuA/Ax2/Y3Avi1ZVn/MofHdlis2WzisssuAwC84AUvwJ//+Z/jqquuQjRqv0SOOTY/dtppp2Hp0qV48MEHxZkBnWSy2TjquczKdcbczZiRN5tNYVkHg0H09vair68PkUgEsVhsr4t3PB7H2WefjUwmI6SsL37xi7jjjjvwta99bc7OZa6MDnxiYgJjY2PYtGkTduzYIQS0E044AbFYDMFgcJqKnkZPtCCOnjpH4/UkR4DtY8CkSA+zejphstrD4TAqlQqazSZisVjHz5qbwL50wuk8FjrwVqsln0/IXjtxokPM6r1eL8rlMvr7+1GtVpHL5SSwIHHO5/M5uu97zNpX3WpWO7GsywC83hjzYcuyLgDw13tq6AMAxjGZif9/AAaNMe+ZYR/HAfiFMeY0y7J+AOBOAH+CvdTQTz/9dPONb3zjoI9/LqxUKnVVutLZz9FuM12DY8kW2jWg2pZe6O02l8/nvpz13kxD7xQz0Rmg3Tntr5GVPd8SyfvzDDBjpvNst9sAOp20zra7Gbelc+R2drEdojR8nZ/Fa8338lkh2302xvfxPYFAoIOroc9FHx9hd/vn8z36HFky0Cpz3FZ3IpAXcDhtvteBCy+8cJ0x5mz763OVob8cwBssy3o9gACAmGVZPzLGvIMbWJb1HQC/mOX+vgjgFgBr9raRz+fDBRdccGBHPMe2evXqBXMsh8uca7DwrsGOHTvwwAMP7FPRba4ycCqJzdYot0pRE9Zwo9Eoent7MTQ0tM9sfKHZ3p6BRqOBQqGATCaDNWvWiIwqVdOAToEenW3TdNDDNi+S06rVKsLhMKrVKnw+n0Dcfr9f7gt/Zs28Xq/LsBWyy+koOWWtUqnIvoj0RCIRIcgFAgHkcjkEg0FkMhmsXLlSiG+ZTEZKJszu6byj0SgKhYK0yEWjUXg8HumTL5fL8Pl8MMYgEokgl8vBGINCoYBwOCwBGpXl2I1xuMW2Dtc6MCcO3RjzSQCfBACVob/DsqxBY8yuPZv9LwAbZrm/py3LehLAZQDWzsUxOubYobZnn31WNAVOOOEEnHzyyYf8GFje2RtSNNckN20zZeyEVwmxshbe19cHn8+HwcHBI86Rz2TUrt+1axe2b9+OTCYDj8cDv98v8PlMZENmtfZMnTVwXSPnNDM7bD0xMSGsdRLQ2F9OmVXul0x0OnbWwYmWkJxHp+r3+1GtVqXGnUqlkMlkhCFP/kI4HJZjIEzO3nWOcs1kMiiXy7AsCwMDAwAmdRAI3TcaDcRiManzk0jH3njuLx6PH7MkufnuQ/+yZVkvwiTkvgXAB/bjvV8A8Og8HJNjjs27rVmzBrt375bf169fDwCH3KnH43EsXrwYO3fu3GuWzoV2rlnsdmdOkRIuypFIBP39/YjH41i2bJkQvY6kxXh8fBz33HMPLr30UnFgzMSfe+45PPbYYwgGgwIJB4PBaUp0M5m9v56OldeQbHXWpJmta/g7GAwiEomIY6dYCx0+28CYwbP32+12S52dMDbLFaxbs2cfgCjEARC9AorM6OfL2jOQpVAoCDoTDoflc0l+M8agt7cXlmVh9+7d0ttOlj6Z7hSi8fv9GB8fx86dOxEIBHDccccdUc/RXNicO3RjzGoAq/f8/Of78b4tAE5Tv/8RgNOY7dgRZ0899VSHM6eNjIzMu0NnNrhhwwbk83m4XC4MDQ1J3VLXNe1qa3NZd9S1UB4Xa+Ttdlsyw/Xr1+MDH/gAli1btuB6wLVKXrdje/rpp+Hz+XDffffhl7/8JX7+85+j0WjgrW99K7797W9LCSEWi01z4HS4+zJqn5O5zjo44XUy0fWEODpfnbGz1cvuxNk+RicKoON3Et74eXTKbEMj3O1yuRCLxRCNRqWvXAvVAJMOv1qtitZ7OByWLJtQPK/L7t27pXWS89jp7NlrzylxmUwGxhjkcjnZjsQ5x6E75phjB2VPPPFE19cJI861NRoNPPnkk6L/r7XJ3W43duzYAQCy8PN1ux0MgdPeEqeduRa4obPYvn07BgYG8KlPfeqQS3jm83ns3r0b2WwWlUoFuVwO559/foey3Wte8xq8+tWvximnnAKPx4MVK1ZgYGBApEefe+45/P73v0dvby98Ph/e+MY3yntdLheSySSA7m2A3GYm08QxOlJqrvf09AjJkW1fnFSms1sS3ZLJpNTLtRNvNpuinx4KhRAOh1EsFkXRLRgMCuueI1IJzVMBjvrqdPC9vb0dSAT13vVIW/aix2IxCS5dLpcIxzCDLxaLHQ68UqlIlj88PCwCN5rpzjIBM3fqzLN+fyyY49Adc2yOTTszWiwWm5PsfOfOnXjyySdlMVu1ahW2b9+Obdu2yTa63qqzLo/HM82x6F70g7Fu58xJWbovmp91+eWX4/TTTz/kWXk+n8fq1atFkIR277334sILL0QqlQIwORCpt7dXeqsfeeQRjIyMoFQqYcWKFXC73ViyZAmA6cgGYeD9MTpBrcAGTF5Xn8+HWCwmvdz8vHg8jnK5LHA6RVjYklYoFMQZ0ymWy2VhgutecQYBPp8P1WpV6tZUkNMBQb1el+0I6VNchigCM/Vms4lCoYB4PI5AIIBarYZ2u42JiQkkEgl5bzqdls/hdadj1jrw9Xq9w9kDU2UcTQAEgEwmg3Q6fUz1qDsO3THHDoHRMcxk+4J3gUmn8vzzz8vvtVoNDz30kEzN2pvtqzaunVI3vW272TNyHRDU63WpjQYCgY73BQIB9PT0IBwOH3QQMRvL5/PYuXMn0um0kLpIvNLWbrexefNmcejvfve78cQTT3QI7ixatEi2nynzno3RWWudcjpyr9crmSsDAz10JRgMirobyxa6vczlcnVMKgMmr3kkEpEBLJRXJRTPz2EpRI9aZWbO33m8dMCVSkXU+FqtlpRWAIgIEIMD9qIXi0WUy2VhzpNkx23pwBmg5PN5+axkMilCNABE4Y5BSbce9WKxiFAodEz0qDsO3THH5sn0iM92u43/+Z//weWXXw5gCvblYrp+/XohVPn9fvT29mL58uUYGBhApVLBhg0bsHPnzmk67My8Z4LL+fd92VwQ4er1eodgCFnTdC7sJaekqFY5mw8bHx/HL37xixn7ku1Byfj4uARVqVRK6s/aqR/IdeL15+exbYztWMzO6fjo8LVyGmePs2TBLByATDvjaFcGj8zcjdJo93q9qNfrqFarwnYHIMN6uE9C50R1mKnzntK58/zC4bCUe6g0x88io57wu66Hh8Nhybq9Xi8mJiakBm+H1MkDoDErZxkBgGzHljrOY280GiiVSli8ePFR7dQdh+6YY3No3aBnZrzUt87n81izZo38DnSO9WS7z86dOxGNRmWx7TYOVPcrH4iz2d/3aDEQbZQE5TaEcBOJhPQxsw2qp6cHy5Ytm9eFtdFo4JFHHoHP55t2zegk7bPcR0ZGcOutt+K0004T9rbmIszG7BwFBi6EzsvlshDdSFDTsroUZaED1vA9W7VIVGs0Gmg2m9Lv3Ww2kc1m5fpns1m4XC6psbP0QUi6myOkIAwDBL/fj1arJWiAng3P7J6v0fFSJpbtZDojj8Vi8Pl88iwzCNCTJdl7rp8nBjPUpOc1YoADQORoXS6XfD5RiWAwiGKxeNQPc3EcumOOzZFt3LgRhUJhrxO+NmzYgG3btqFQKEgm1k1LnUFANpuV1/X+dPY/WzWvubBuMLmexMVj9Xg8iMfjWLRoEVauXClQKluM5ntRzWQyGB8f78roZ1Di8Xg6GP48pvXr13dMIpuN0QHRMZLExoyXGTg/k//rrJm65pp8ptu82HNOaF4T0OyZOX9mMEUYWw+6IeNci8m022309PRIZkunqTP3er0u/fN01twP98lhLnaGusfjQSwWk/Y2wuksaViWJb3nLA8kEgmMjo5K1g5Ajisej8tx8BonEokOsRpeP8uykE6nJZg6Gh2749Adc2wObGxsDA899FCH4hdN16Q5/UxvM5PT6NZOtjeW+uEwZke6TY3OZmhoCCtXrtwrL2C+LJfLdR0oYx80Y5+XDkzVtoHZcQ8YyAQCgQ7naFmWBFscYUpHH41Gpa+aGTwRARLb6PQikYgECdxfq9VCOp2W/bKWrn+mgyUSpAMcOm9m7YFAAOVyWTJp7YgJ29MxahRKE+k4DY0wOJ8NZvM8tkgkIlrxzOpZA282m+jt7UWlUpFgge2XmvFPBr0+Vu6TQQuDCD57RASo+X40EuUch+6YY3NgjPxnMjp1TXjSpv9mr+3uy+ZrChqtm9obJVt5vFo0JBqNYvHixeLMD4fZAysdVM3koGcTJDGzZo80IeZQKCQOhH9jHZntVnQ4brdbiHlacpWO0ev1CiGMx2RZFnp6esTJ8zjYbsZgQQcQXq8XLpdLjo1ZM4liup+dZDbNXiecXSqVRGZVs+V1Vk7ontecDpxCQayTM2AYHBxEKBRCOp2W4GXXrl3Tes95zYgoEI6v1WpIJpNSsqjVakKyo/FalUolCbQYGBytRDnHoTvm2BwY+4Bnsr0xurloM0Oj89FTrA5nRm536LpNiPVQADLe9Mwzz0QqlTpsC+Xdd9+N//t//y8uuuiijtrsgRrf7/P5pFeabWF6qhidMTAVUGjyGwBxosxkCbMTTtYEtEAgIHVjy7I6Wu3Y/81ea00+pMOjDCodsGVZHSNQAXTA5NVqdVq/uc/n6xCA4XvsWTmAjvY1BhrpdHqaOEyj0UA+nxc1uWq1Ksx7HVBwOwDStsbzq1arQqRjlwCVBwm18zukuwk4eY/DcBhcHC1O3XHojjl2kJbP5/Hss88KWaibsbWmm+mae7fM/XCbHXkgnMwAhjVd9gDr9q7DYfF4HD09PfucqqazZ83O11PImPlyDjedJVu6qH1O8lg0GkWlUkEwGJQZ4fF4XPqvWed1uVzyut/vh9vt7tgW6KyFk8BGpjizc9aUuwUGDLD4TMbjcRSLRRl7yrpyIBCQNkISzuicOXfczsAnGQ+YKjtokpydzc7rWy6XxTGzDY9BAD+HAQWvPe8jkSASBvX0Na/XKy1uvF4kHVKznu14JB4SnTiaxq46Dt0xxw7CGo0GHn30UYEFad1qs8DM9fJumfhCcOZ20+pe/J3ZZzQa3a9Ja/Nh6XQaTzzxBF75yld2TC2zX0s6UR2U0Lnzd01WZKsZB8rwfZRHpcMuFoswxqBUKkk2ykwegGTrtVoNpVKpo8xSr9eRzWbFAWvGOrN01p1brZZk3T09PSID6/F4UCwWZeLa2NiYOCoKxuh9h8NhcaC635xku0ajIcSzTCYj8qwM4OiQOWAFAIaHhwFAau48H/690WhIX3y3fnLC9LVaraMcwICFBD8GPFpwh0FYtVpFoVCQHnxN6mPHhX7v0ZKpOw7dMccwqbOeyWSQTCb3S6I1k8lgdHRUfp9tr3K3OvlCdOB2mwmB8Pv9SCQSSKfTh/iIpmzz5s34wx/+IM5opuvp8/k6sm22FJKcpZnwusbdarVkPCxRFRLLotGoIBeVSgUejwfBYBB+v1/q2JVKRdjz9sAnHA7LawxEmMXrrJr74NhRAMLo1j83Gg3JlBOJhJyDvZ5uz8yZ1RIFYNCgIXO2zzGjJsGOJYBWq4VCoYC+vj5x2HoCHBn9NA2Lk5uRyWQk2IlEIkI0pJaBMUZY+AyuQqGQXMNEIoF8Pi/lC83sJ4KiM/VWq3VUtLQ5Dt2xY95GRkZw7733ApiEWc8///yuTp060Xoa2Pbt2zuyLKAzI7c7lW7bzDepbb6Mjo911qGhIdEHP5R2yy234Lvf/S7e8pa3dJCpuhnhW2a+0WhUsjpC7XSE+r4SSgcg7Wxszctmsx2SrO12W7alE83n8x2DTlj3ZebP54oBAX8uFAoCwY+Pj8s+yfxOJBLToHgKvzBLZxDA7JROmoQ3ZuZaGpbnScjc7/dLOxmvXTKZxMTEBIrFIgqFgpDOyB3gsTO79ng80m6mofJyuSzHAkwGXD6fD6FQSPbBoIsIAQVneE+JPmhSKcsCdPK8ZywxUJqYMP/RMHb1yFxJHHNsDo3OHJhcjO+//374fD4kEgkcd9xxQtj5wx/+gGaziYGBAZx99tnwer0YGRnp2Fc38puehLWvbY9U83q9As0eSqtUKti+fTve9ra3iaPt5szpuKPRKKrVKsLhcAfkS2iZ2ui8Ny6XC9lsFrVaraOmTlnbarUKv9/fodZG8Ri+n5PJCIdzKhilUt1ut9S8G40GvF6v9FJ30xtg7b/RaIh4DDPlRqMhpLdGoyFCMVToI7rAGrrOzBk4xOPxDnjesiYnt6VSKfnMUqmEcrk8LXNnSYIDYNjCyNY58kyo3sbrSifOLJ1StxMTExKEkQ1PSJ2Ev2w2K10EiUQCjUYDfr9fZHwLhYK8x647H4lE5DnauXPnEd/O5jh0x45ZazQauPXWWwFMLTT8udlsolKpSD1Qb7N9+3Zs374dfr9feolphGrtAjFanYwL+NFkdAaH0n7yk59gw4YNOOmkk2SB7nZdNfmNIiu5XA5ut1ucXqPRQLlcxsTExLRhMqw1ky9ASD8Wi8l5s42L2Sj7qumMtEIca892iVZmsVQ6Ayah+Gw2K3VjypkyC7e3q9GBh0IhydJ5zuPj4wKv6xq6ztoBiCqclobt1p+udRKYeUciEZnyptn/ut0sHo+jVCrJ94LnzyCCDpkBBXkDAARpYGBBY2kik8l0qMyx/FAqldBoNBCLxQBMtf2R7a6Z+0cySc6ZN+7YMWl2Z64XcLZl8X+9jT0LITNdOxIyenUGPpts/GDbqw618ZxZz2RNd76tUqng6aefRqlUwimnnNKhOW431lz593a7LfKjzM64nd/vl5pqIBCA3++XmnOpVJJWKTp+itcAkCwcmBrkovXHqebGWns2m5UaNMe4kp3OZ8cY06EoCEwNdqlWqx01aSrOVatV1Go1ZLNZgb3JiGeWnkgkhBAYCASkXz6VSonDJixNI5mP0sSlUkkEcGKxmNTQiQ7QmcfjcYTDYWzatAmf/vSn5TwBSMCQzWY7RGgAdHQS8JyIsmjYn6+xFx+AlAd4vPF4HJFIBF6vF7lcTsh27Lvn97Ver2N8fBw7duzArl27DjvB80DMydAdOyaN40Z1rzfQ6dz5v93hz6a1zF4X31ednKMzu4m4LHRjRsoa73zazp07cffdd4uYy95IiOwF9/v94nxIztLENt7fZrM5LUPXDpZQOwBRfIvH4x1tVnrsaSgUgtvtht/vRz6fl1o5HUgsFhMoOhgMYnR0VLLxRCIh5xEMBiUA4DGxTY1Oh+fEnwF0THKjqhyV6Lr1pOvMnEx3BrGE9EulkrD4I5EICoUCQqEQ7r77buzevRunnnoqenp6sGPHDkxMTODJJ5/Ev/3bv2H58uVS1ohEIqhWq4JQUGs9n8/Ls9TT0yP3isgGAx5e23w+j1KpJCQ9v9+PdDot8ssk/FHrHoCw4wndM/CIRCIS+BzJJDnHoTt2TFmlUsHatWs7mOk0u3O321wy0+2Om87gSHPm2vY1IvZgLZ1O45577tlrOxoAyZQpi6pJWXwPoVgu7pFIREaIsr7OgSoABJalcyeZTbd42aV6iQ6QqBYMBjvUzwj7M7v1+/0ytSyTyYgjq9VqUmsnyYz1Zo/HI+1bPEcKtpB4FovFUCgUEIlEBF3SjPlu9XS73nyj0RAuCa8BA4tarYb169fDsiysXr0a559/Pn72s5+hWCwiGo1i+fLliEajuPbaa+WcKPwSDAaRTqeFy0AiH0tZFObRhLZMJiMlDr/fL6NgyT8IBoPI5XLSCkdRGrYWMsPX9XPK5RJdqNVqguQcSU7dceiOHROWz+exfv16qYnrmrmuec9kc62hzmzRDsvPt0NniWBvMrUHamwjeuSRR+B2u7FkyRIhJh2MtVotXH311UgmkzjxxBP36szZLkZdAC7cdJrA1EhZzbSmEpsupfCc6NRZm2ZrGwMDjjn1+XyoVCoolUoChZPZ7fF4RJ6U1z8YDEoLGYMOQsdshaOKG9EBlg7oWEksY2DA+6qfIzLs+TmWZXUw5jmulEx+O9OdpDY6TdajGZAAwFVXXYWNGzfi9ttvh9frxVlnnQVjDBYvXoxTTjkFAwMD4lAJowOQLJgZOj8HmOr95/AaIhwul0s4DQxKSNijQh75Abzv1Apg/Z3XkAESx7j6/X74/X4Ui0UMDw+jVCphaGjoiHHqjkN37Ki3DRs24Mknn5ymyNZt+Ek3m6vaNrNyOhAumjwOLsIzwe76OPTf6VjsTtoudKMXN62INle2Y8cOBINBPP300wCAZ555Bi9+8YtxwgknzHof1AKgffe730WxWEQwGMSJJ56414BKi71Q+KRYLMoCH4vFxJnwGWBWymtDaVdmv1R/65bZhsNhySRZAyaEy0EozAypAMd2rFqthnK5LA6IzofnoRnsfD7sA1+08ybLnnVlapUTNeFx6KCGwYxmz2tVOZLUtBwsyXk9PT2yb5ILV6xYgfe+973I5XK45JJLBP2oVqtSMuB1bzQa0vbGdjQeE+vkRADYqsYuBJIGdTbO+0AUhB0M5AdoYRky7L1eb8f0Nl16IEmuUCggHA4fMSQ5x6E7dlTb2rVrsWXLlq6zxPdlc8lE13V5u+qVnrsNTC222sHzddZfuT0XL7vz1sEKHRezs705c/1Zsz0v3VfMzwMmM9p169bB4/Fg+fLle93P/fffD5fLha9//eu46KKL8OpXvxoejwdbtmxBT08PzjjjjL2iJJxfTnibJDYAwphmdhwOh+V68FrxmpONzRq1zmzdbre0SJEpDkw6PzqiUCjU0bbGPmut2KahfRLJSF5jNkydd83QJ4Od56QzTRL0WCfXmTSzWR4jM1XC5QwgeIwcf6rf6/P5Oq5VrVaTIGbXrl3w+/3iPBk8sj2P8DqAjpIG28kajYYIwBBiJ8eBI1jZJtctGwcgfeksgfA1TqrTwjJ8LohysFtBozucDMcafblcPiIydcehO3bU2sMPP9zhzLuJuthtPtvJuKAyy+LCx9e1k6XT4HFzQWItVzsk+5Q2Ou1uimSsKWvIma/rISLA1GjUfZ3TTNeTjv3RRx9FIpGYsU/92muvhd/vx6mnnorLL78cwKQ2QLPZxIknnigktb0pv/l8PrmGvC4UJmH2x/Gm5XJZWM2RSER+ZxaqR5cys2WfOO8LiWrAdBlVjgElhO5yuQRCZ5uUz+cTx88+9mKxCAAy3CSXy8lxMGgKBoOSafNzWJdnvVhntyT68fkhxEw5VWBqGhlZ8OVyWcoKfAbZLseggY5dq7MxG2fLHrXg9XXV4jCsi7OtjufJ7DkSiSCfz0v5JBqNSg88rzcDGgZr9no5iX58HskPoFANyyB8PvTIWpbAmNUfCSQ5x6E7dlTa1q1bsWnTpmnOHJjKYO0Ofj6duf0zmdnSCWknwr/p/nUAHSQl7ofiJsAUJKtRAC5K/AwNzXMf/GzC8SR5MVPVPfV7Ozd7bZtOoV6v47nnnsPZZ5/d9f2rVq2SrJnm9/s7+olnujesKWtEg9KrbE3i4k1HQOfXaDSQyWSEra4zXF5b1myZ5RNRoZOgjjq3d7lcoiHOexsKhaS/mY6YfdyEuHXbFrNIlmfsmTfvje4T1+1x2qEzk2fGTyfNgIPPj3bkdG6lUkkQCTL2mVVT2IZkNu3smIVruVZm4uz7Bzo5IwxMksmktD9SVIbBDvvTE4mE8AX0azwnktrs8q68LloPPxwOS+DLOe3FYlECM5ZT2u02xsfHFzxJbtZ96JZluS3LetSyrF/s+f3/syzrccuyHrMs6w7LsoZmeN8Wy7LW79nuYfX6kGVZd1uWdZtlWZE9r11vWVbFsqx+td38UmcdO+rsj3/8Ix566KG9wux2Bz/fQi+6DqpFZuzHQGfDLJswLp00HRIdGBd0HQRQnERntkQFWB/W/1M2lH3YXNjYcsVj4X67mT52uzWbTWzZsqWrznulUkEymezYrw409nZveK7US6dDc7vdSCaTHUxwzVnQSAjrv8zSWIdnGYSBELNdHg/3USgURJs8k8nI/eXwFULD3A/JV1SfY5Clp77RWetnRwdarPkTlibDXfeVk7zHrBqYmlHO17SmAoMVt3tySEwul5MAoN1uI5PJSOudzvB5HYgE8PdMJoNWqyUzygln6/ZNBpBerxfJZFLEctxuN3p7e2VIC4CO51mXm2i8zrrsEg6HBXUiCkP0Ix6PIxqNwrIs6anXcr18H68p7/UzzzyzoPvT9ydDvxrAUwBie37/O2PMZwDAsqyPAvgsgA/O8N4LjTH2JtWPArgKwPEA3gHgW3teHwdwLYBP7MexOeYYgEln/swzzwCYymgPhe2NyKbbfJiV0TnrjInOiOx3ncHoOnij0ZB90Ikyi2S9ko5OM+n1iFB+Hp2IruEDEFUy9jfz//1h4evMvd1u47HHHkNfXx9GRkbQ09OD9evXo1gsIpVKdUVR9nW9iWDoe8xMmux0niOJURz+wevCgCYUCol2OZ0mJ6axLcyu4tbT0yPSpYSW6YCohsZWK9ZiORQlFAohl8shl8tJpszpa8zEAUgdnZBzOp2WaxoKhToQKAYOhPVZ+2YWbecM6D52BiK8LsymmQHzGWDLF8sQRHY02Y1BHRGOQCCAYrHYIdbkdrtFVY4sfDptEuHq9bpA+oFAAF6vV8RmNM+Ez6bWaNfoAR0/n3V+BssyvN925Tu+X5Pk6vW6DH1ZiNn6rBy6ZVlLAFwK4AsArgEAY4yWhQoD2N9+GzeA9p5/Oqz/FwD/27KsLxljDv2kB8eOWNu4caM480PlyGkztb7ZmenaweuFUDs/bkMonNAlsxnWGJk96bq4z+eTvmP26brdbsk+aHZEQBOaarWaOBH2LvMz9lYz35s1m02MjIxg9+7dsCwLY2Nj8Pv98jn7IixqUh8dOQVAotGoyKGyVYwOjw6U2aE9MGIdularyXlqKVSS33QgxGyb7WkABM5lFqzfQwiZ98UYg1wuBwDSjsa+agYNzHS1WAyzegYNdFx0QHTYGtUgvEzHz2CQx6ufSTpqOm8+W3Tkmp/Afel9a5Inn1G2AWqhHZY3iAzw2ff7/QKN8zrzGOlc2THA7wXLCrwmLF3oDhKy4fm3WCwmrXoMTPjdCwQCiMfjGB0dFc5CqVSS42y329i5c6fU+Bea7vtsM/SvAfg4gKh+0bKsLwB4J4A8gAtneK8BcIdlWQbAt40x/7Tn9ZsB/HDPe9+uti9h0qlfDeBzszw+x45x27hxI9atWwdgurLbfMPpNJ29aqgb6BSt0Zmhhi6ZPQDoyHw0+Y1ZOM+L2TgJW8YYISbRwRE6Z/ZFx0bYnY6NTpFQvyZUUXWMx8Fj3JfpxZIOzv732RqvJ4MZIhx0xmwpoxIZmd/MvMkKz+fzErzE43Fx+HTm1PQOBAKIxWIdmuuaQEWGNVnrJL3xfgDoIM/RsUejUeTzebmmJMyl02nZJx0lnanuJKDD5THourmeCMf7xPdqx6Ohct4jEuB0FwUw1YrH58eO5ujPAKY6L7QWAGVoedxED1qtllxXkuoYrJXLZeRyOXG03G8ulxPkg98ZBnB03OxWYDsa90vjM8FShcvlEtW/Xbt2ybPPmjv713mPAUgWr1X9DrdZ+4LQLMu6DMDrjTEftizrAgB/bYy5zLbNJwEEjDHTHLBlWUPGmOE9dfE7AVxljFkzw2ddj0mH/l0AjwE4HcCwMSbSbfvTTz/dfOMb39jr8R8qoxTisWyH8xpwBCMw3bHOl3WD2fVrM2XsM33ndPbYbb+6PU23q9Hha1KcFj7Rr+lj4+d1EyPh6xrW1DwA+/azOb+Ztp+N6e00v0C3//Ea6fPm8eqAgtvo99rbBPW11GbvRqAT1H+jw9XnrFsJ9ed0Q3b4mj5ufV35v0Yt9Lb2c9D77HYP9tbGeCCmj9l+3/is6iBVI1L62bcTOfX7dZcHgA64XAfKDKx1D3y3e8H36ueh23PD66o1A4gyaZvvtfDCCy9cZ4yZxjKdTYb+cgBvsCzr9QACAGKWZf3IGPMOtc2PAfwSXTJqY8zwnv9HLcu6FcC5ALo6dPWenGVZPwbw4b1t5/P5cMEFF8ziFObfVq9evWCO5XDZ4boGjUYDv/3tbwWC45frYDNzLgDaSQBTC2A4HBZ4mpmRrs9pJ6EXTb0Y6TY01hCZETNDoIIVxz+S/Vur1aRNSU+KYy8zsxdmb2TykjxEchSdEuuV/F2zgpnts5WIcL89W+vWKkc7kG4C3RWgs0xmZSwzsI5dLBaFXKXb2KhExgEmzOTIhCYMTHSj0WgglUp1DFZhvzRbpJrNJqrVqjwbZOoPDw+jr68PbrdbxGFoFKOJRCKiCqcdL0seOoiy8y2ICDBb1tm43o/mCWin1O053x9uxL5MPwO6rU0/Z9Rlt6xJwR/eN/7d6/Uin893BDzxeHxaJ0dvby+MMYJ48BllX7/bPTnLnhl7o9FAKBSSZ4glAb43EAjIsVFEh6hMuVwW9j6Ps9Vq4fjjj+8YZAMcvrVwnyx3Y8wnjTFLjDHHAbgSwN3GmHdYlnWi2uwNAJ62v9eyrLBlWVH+DOC1ADbM8ti+AuADcFrrHNuLNRoN/M///E/H8ApgbmB2Olm9Py6MOiLX0T+3ZQZIB8QaH1/ngqwDBpdrUn4TmMoqo9FoB4lO11gZFGjtcpKqqAtOwQ06IEposmYJTI3cJPQei8Wkvk1omk6TxCvd+zxb259uAjounm8wGITX65XpZ5FIBD6fDz09PQK16lYqqq/x2rA2TKdI58ESBfvV/X4/YrEYXC6XTFYjK50ZI/ueGRjwmLSASqVSQbFY7GgdIxO+VqvJIBJCu7r1kMesM0yNItiJffYaNlECu5OeKQufS2fezXhsfNYZpPC8eO6hUEjEdkgUjEQiSCaTHa2VWsqXjtyyJnvltVYAAyxN+my1WnJfSZxk8FooFKRWzy4Dwusk3PE55LG321MiPAvBDsZZ3mhZ1kmYJLVtxR6GuzXZvvZdY8zrAQwAuHXPl94D4MfGmNtns3NjzPiejP6vDuIYHTvK7amnnpJaHW0uMnOdMfM1ZkaacKOJVlrwg194Zti6B1jDomxjYmamFxQtuUnGPrPESCQik7l4zoQtuZBRq5oEJdYdaVyE6LCZnegWH56z1rzW9VH2GGs4/2BNt+6xzk+Cn3ZqPE6eMzN41v/HxsbEOWpBEU2y4j0mO54KaAyISHxjlkYnwAyYGTrvA1nYbPerVCodvdzsi6aYis7QdUeCFpRhEMjnmgpqhKK5nTb9+3w77L2ZrqsTDdEES7fbLcQzZsZ0pgwigalRq5ZliVgPURZeS2OMBHJ00iQoWtakVj9r4NRYIJOezzD3qUmo/J/PBr+3wGTgNj4+jomJCSxZsuSwE+T2y6EbY1YDWL3n57fMsM0wgNfv+XkTgDP2Y//X236/BntY9Y45ZrdKpYLnnntOfp8LR96tds2MmIs+Xye0rmumDCx0y5zOErRwBzCVibtcro7Inw5UZ57cpzFGYHCOfSRUTMY238uAgn/jZKqenh4UCgWBRynIQccDTC3GkUgEgUAAuVyuw2HxWEjeOtgsRaMhzEyZTTOY4DHZmdo8JvZkUzCFThWYIhqS8U6nQlJZvV5HMpkUJ0lRHAYQzIA1JBsOh2X+NrNrXicSCpmdMmvmvdDHoGve9tnjNB6Dzso1MrQ3bsPhMl1v5j2i4hsDXJaHOCue2TUhcABSRiHRjciKrnmzLESo3A6pN5vNDiW9iYkJCQ7sIjLFYlGCWxph9lKpJB0GbGcrlUrI5/OHXfPdgbMdO2Isn8/LSMZ4PC6LAhf/gzWtxEZHR2elYULWLAmnMnu3k520cAfZ6XTuzH49Ho/AuswGPB6PZOBkMbOuXqvVhGxDWUv24BJN4EhJZtGECwuFgkySoqQmJ2wxg2fm22w2RZHLsixhX3OyGIMATTQ6UGJVNzIXF17WjnXQwHonnR+FYch217K2vNYMUji9rFQqybAWljU4lpNZYqvVEng/nU7LOE4GBLwGhMx5nRhM0emyBY7HS4dERIcoAxEEljWYJWqnCHT2VfM5W6jGzgKaMUaet3a7jWQyKVyCarUqcDodKb8zDKD5/WAfPzDFwuczkc/n5T2xWEwcNZ8lIgD22jpRALfbjXK5LMem0QAeD59BBvhU2atUKhgcHDy0F1mZ49AdOyJszZo12L17t/y+aNGiDinMgzGdDWpnQEdqZ4/TgXOhAKb6fQEIlEjnrxm4zBSAKTEPPZCDjp0LDUlYGtZlFkJxC04UI7Q/MTEh2SYDES78dCbMdukoGMAQvtakJQYtdFzMcAg7a4e7P9mhZkPzGGjkH8RiMVSrVYFfNbzOyWi8VwwGiFrwnBkI2EsG9Xpd2szY2sf3ulwuCYy42MdiMcn4eD0oZMJrRedCSVmNwAAQh89/PHc6LD5jLHXwWtA072IhZeJ7Mx4/2wNdLpc4ZJIXOXCFz16rNaUlzyCIv/N7wNYzIiW81nxeeI9jsZhk+HqEbG9vL/L5vHznqWnP/bCljtC8vh8MpunoKS9LVOdwmePQHVvw9otf/GKaMMqOHTsATDnjgzF7iw1hXH65uXgygyJ0yIXG7/d3wHq61Uv3fNPBs0bLLz+DCZLaSKKi06Bz5kJBVTE9zzkajUr2yiyP0qOxWKxDdUvX5rlwEtYnm75er6NUKonT4DbM5Mnu1zyBboGVDkS06azevg2dNzMeMvd5vZgNc0IXWefM0ljvZs2dwiJ05DxXLtI6e2dgxn51fd/4HDCL1DPCNbTP68VnS19jZnj8XyMPfI8uzXRz2nPNSp9v08Ex7yVlYsmP0N8tBl5UiyO6xYCg1Wp1sNE1/K2duuZYEL0i0lWtVjE8PCyQv66tkymvOQvGGPnOAegYx8pz0aWtw3V/HIfu2IK2hx9+eJoz14v//rCstdkhYt3+Q8dHqUlmRHSYdMpc5LlgcDFimxNbZ7xeL6LRqCzoXCjoFAnfcmGi4wgEApIl0jkzaKATJ5GLGYw2MrKZWXNRpWIXg5JSqSRZOGd/+3w+aY3jdSKkzACEx8Hz1vAogBkdkl1bne+heAdryCxL0AlbliUa7XSmZCkzU2YbGzMxXUYAIGM+6cwbjUbH9C6S58iWZmbH7NC+yDM4YADC0osmtfF/noP+H5hCauyvdbOFDrHTNPLALDuRSEhgw7/xGSY/ZWJiQpTlKO5TrVaRSCTEwRPNYHuby+XqaD1jLR6AZNB89nmd+T3V7Yt8TtjlAUCQFl0vZylG69f39vbC5XJhfHwcw8PDomR3qOvpjkN3bEHbtm3bur6us+kDMTsRCZgi8BBiIxlLs73pgHXNm19wLuKs6dJJMLJ3uVwCIWsRHPaRs6bLbI1ZpIbBOZ9aO3G20rAmywydBDjWlumkdH8yFboAiHY7rweDDDoaPYmMzl3Xd7sxq4lO6M4AcgE0qYttZJycxUwcgHAKeA+oi04ZVMLcAKT0wEWa14+lC+7X6/UiHA4jm81KlqfZ+3TkLFGQEKWncPG6MMDRcDCPVWul6+vS7VodSVm3Nu28eQ5EihgUuVwuCZwsy0JPT49o7jNwZYDGwIkOl5oPwWBQCHW6Rz2ZTCKfzwOAIF0AZJhMu93GyMiIOHWPx4OhoSGMjY2hXq8LSsB7VS6XhfnOjJ5BK5+HaDSKvr4+7Ny5Ey7XpB49deH5rB+OcauOQ3dsQZuuOdvtYKB2knU0oxrorOsym9YtU8CULCfbvfiPvfBcOFiXa7fb4ly0tCWdKRc9zscmUY3tPIQF6UCZZTOrICGMkCFJfYSSiSIAkMWTGbquT2vCEK+BHsZBwhlfsyui6QxdQ/XayBPg/olycH9er1fIZyw58PP1ORFeJcuczHZjjAQ3OsslwpFMJsWx2CFysqQZWNDBxGIx5PN5jI+Py3lFo1G57/l8XkR6mAXqFkZt+vk6Uh040CmgxN91cEykhGgHAEF9KpWKQNd8P4fnMDvn4BYGtuVyWaRueY2J1IyMjMjz09/fLxwGPr92OF2PW6U2fLlc7mDKA5BnpdlsymS2XC7X4cB1HR9AB4+CLZOH0hyH7tiCtpNPPlk02oEpEZe5aFEDpli0rOfpKJ7ZJOFWZriMxOnQmClzodYiGPwsQr7cLxcmO2TP/Wn9bNbPNSueGQVHnjLw0Bk8nYxWjOPCx3+61svaNJ0N69NcWIlc6ADLTtjS/2vTJC7CnbyfbNdj4EQSGjOnWCzW0QvPLJkLPY+ZQzR47CRKsdxQLpelJ7ynp0dgWl275vPFLJ/3gP3T/MdrxXuha6c6SLHbkZ6NMwOlghoD7p6eHgk+6/W6kFb1PdXjVYkqMdAl54H7qFQqQkLkfeao3UgkIhwPfgf5mi4Fcd/RaBTValWCZMuykE6n5VmoVCpIpVLi8BmUshOCwQcwRcTT+hRE8wKBAHp6elCtVrF7924JDA5llu44dMcWtJGcMtd1Q82AZiZFtiyzK35RmamFQiFkMpkOUpbX60WhUEAwGEQgEMDY2FgHJE8hEV03Zpag9cSJGHBBoNhJOp0Wp8EAgJk0a/IkuXExYrZJx6J7enl+RAW08hvLAqxFE6HQQhuzvbYzmb3lyu12yzWms9S9wAA6hG64D4/H05FZsVNAK9ox0LITHjm4xRiDRCKBYrEoLH/W0DUngQ6d5RUGlHT2MzHOj1SnrU3zHBj06nvHNkaSO1OpFHbv3t3R183rFgqF5PqGQiFUq1Upq3D/zOx1qUTXx3kveBx8XumoidIUi0VBEXSQx/ZLOn3u1+12Y8mSJRgZGUE4HEYsFpMgnCRLButE2Xp7e7FixQqMjIzI9SK7vtVqYdOmTYjH44d0Ipvj0B1bkNZoNLB69WqZVEWbi+xcGx0Vs1LWikOhEIrFomTrhLy5GAGQTJpZLBdwkumAqYWDTpOLkl2BjouKng2t6+eEiem0dc98JBJBuVwWJbhSqdRRSwYgTpoQPzNdDY3bHZDOvg82oNKZKRdiLrjMyogS8LxbrRYSiURHCUTLdXJxJ9tfHzMXawYvRAI0+jAxMYHR0VEA6CBF6XPWOgC8j5r0djQZnao2zpAnOY0ZMANUlpMSiYT0YjNI5j/dCsjniL386XRaoOzBwUEUCgXhOACQa83yGD+XQTcDfnIiyBuJRqOCALAMps/F7/ejUCjAsixBafL5vIjJjI2NIZlMyjHoIIPXIZ/Pw+12i8NnCUcHC8ChraU7Dt2xBWfPP/88HnnkkWmvzxe7l8Qm/mPGS4INSVNanpULBDA1EIWCGAw6SAaiTCQdEwDJmllXLJVKCAaDMr+brWEkE2mSFxcg1sDJZCcj3OVydWT63Jcm7jGjBSBsdp4nz4kBzkwchgMxTZZjpqNbwjS0X6/XZdFlm5/L5UIikZCARZPRGDQxeyZ5LplMyudqzXodxJDrwGBRs9SPBFb5XBizcV0qYgbOZ1qT/dijDUCeJV1TZ0DFzgWqrBHyppobnaUe+gNAtO+JuPT19QEARkZGRHaXNXoSSslKtywL/f390lam9db5HMViMUEBLMuSzggeD/kT5XIZExMTQszU9Xg9XpjIGHkf3MehrKU7Dt2xBWUjIyN49NFHu/6NzpYOeK5M9zBrEQxmAgA69NDJgLVLi9KREGZkVs+/caFkJsOskXC5nTHPRYe92ACkXshj1uUIZg+6bUuTdIDOVjKytAmPco44oXi32y1qXHNt2kl269ulg9DXgc6AtU9mb8y26dw1M73RaCCdTgvXgPdMw+k6az+SBFsOxpiN6yFCdEz8neUEt9stjo1sf3JEPB4Penp6pCc8Go1idHRUJvox0CUMrlvFiHyUy2VhpjNQ1lC5rrPr+jg7UtxuN8bGxmQQkYbnQ6GQMOE5ra1WqyGbzU4TmeEEPB4bWfh9fX3IZDKSZROBIMrE7ysRML/fLyJFPK9DZY5Dd2xBGSPsQ2macUzHSYdGaJCROJnXACS71Jk9gI7paISXdV1bQ3PsM9f96LovmWQgZtlED1gDZ+2Xn01EgO+nE9PkLmbkhBt12xqlUUlgOhQiGTOR6Oh4mb1pJTHdJ0/nRNSBzoYlCQAygc3eamf//GPFmbMzgLB5KBQSFIc1aTplY4woqtEhs0OB5EoGU3w2+Sxrhjm/NyQ96jIVuSZ8vjmAiN9DrUyoCadUU6zX64jH4x1lMcuykMlk5Puh9RLsNfRAIIBEIgGPx4NcLoe+vj4J1vmPAc2iRYvEmdOR63WDKEa5XJahM8uXLz8ksLvj0B1bUJZMJmdFgrP//WAydr2oELLlolEqlYRoRagXgBDWSJYCpovc6AWTGQbhODohOh3Czj6fT1TfAEjG6fF4RAOb9WMukGzZ4nHojJwLIAl8pVJJmPB0+CTd6dqo/docDmMQojsStENmfZzb8m+8Phq65bU/lo1Okc6XjoeOOhKJIBaLSdkHmMpGJyYmJNvkaF22VWqBH35XyDqnNDG7LThDnKRO3SrWbDY7RF2okKhllguFgpBIE4mEfE8ZwFHFj0Iv2WxWkoR8Po/e3l5pNysUCvKdZx2cDHaOxmU5TDv/SqUiKAIAuVb8HApK8T3Uk3AcumPHnA0MDOC8887Dhg0bUCgUEA6HkUgksGXLFgCQKL1bS9CBGrNo/tPZHWFxAAL5UnCk3W6Lshrrt0DnBDVdF6Yj1YQ41iP5Gj+HNXFdw2WwQdEa1ijJiKdzI5Od0ClrffxnV34jAkD4UgunHG7jOfF8umXfGjnpRlY7Vmrg3UxzCviMcfAN73l/f39HFwWzcPaM644QAJK9M7NlwKjFmFjKYd+2MUZa3Rh8DQ0NCaExGo2KoBAH5fC+JpNJ4aNouJ0jf5vNptTzKThTKpXQ09MzLRPn+QOQ/XI2QKVSQT6fn/bZrMMTDerr65NrqYcCETXTmha8Joeqju44dMcWnA0MDGBgYAAA8Oyzz+Kxxx6TvzEq1gSzuaqnaxiWNXKyX8l8J3RIZ8/WKhLfyuUyIpGILDya0MX9c1sy24EpwRWSa+hUNZlNt0mROcw6JrNVZiU0zoKmUArb0BhwAFNz3bXAi5ZHPVymmdK67YgICh0AiVfHAlw+W+M9Zw3c5XIhHo+jUCgAgLRIcrIeg1YGTHTmdOwcYMO6+NDQEEZGRiRYZXBKEiOnAmoCGdDJGOe0P8uyhNTWDQ7n8QUCAWQyGVQqFSQSCfT19aHVaklboyaokURKIh5JpMaYjnq4LmHpGj330Wq1EIvFEA6HRR9eo21afKhYLEo2blkWjj/+eGQyGUEEDoU5Dt2xBW2bN28GgK6w9ly3r+mMWNfByW5nTU1n43TEbF1rtyelSfXMZnvrDbMUO8RN58VMU9fk9fAJsn7plEulksCmhPd53PY+aULVun7J9/Ca6qDpcBoXWZZgGFwQch8fHz+sx7dQjdkjA0kSChlAUomQWTdRGY/HI9rq6XRa0CrOMOc2zLr5HPMzOY6UkDtr4Pp51KJBzWZTiGnDw8Po6+uT9xUKBbjdbvmd3wkGbqVSCel0Wr4/DDoKhUKH3rvb7UZfXx+azabA8GyN05l4X1+fqC+ydZJiMeS8tFotFAoFGdqTTCYRjUalRBCNRlGr1eRYCN+XSiVpo5v3ez/vn+CYY/tpW7duFXEKajRrm0tHro11Lzo1fo6udeseabaEsf+bRmepe6J53HTs3B/Z2VokxhgjsDn7flnr4+czm9bXgmx4XVfm8XNxZtahiXeEsBnMMBA5HKYJijQ78/xwBxoLyTRKxftsjBGnAkAcXDgcRrFYRLVale1DoZB0N2hlN6rB1et1GVqTy+WEBc+smJms1+tFKpXC8PCwoAI+nw+pVErY42xBKxaLSKVSADozdsuypLZOqJsSsBSh0TX3VqslCm98tpPJpCg5ks3eaDRE78Cu6qczca/XC7/fj6GhIWSzWam3M0inGJFGDnhcOjghBD86OipBgzEGQ0ND8+7UHYfu2IKyrVu3Yu3atR2Ltq6BzpczBzodIAARqtDHQuiSNTfWdfmzNjpRGh0TM3bC2rqmycWQ8LnX6+2QZOXiRCOZjsS2bg4RmFrU7UGGbmPTrx8u63b8h/uYFpIxmCQJksEg0RoK5JCBTjg4l8sJubLdbiOVSoljo5gOUSdOMKOTC4fDHcxtih2xA4EEu1KphEQiIUIv1F/XDpQOmy1dWtyl2WwinU53QN+EuAEIgVNn/MViEcFgUEbp8tio559OpzvIaosWLUI0GkU8HkexWBRCqdvtFsIdgyF+FhnyJKTy81l+YHkinU7LUBoiBAwAgEMjMOM4dMcWlO3evfuwOHNt/HzNiuZCytcpE0sny0wagCyczPZprA/SsbK2SadK9m4wGOwQOGEGzc8n5MzMpdux2+1w1sL3x46U4zwUpoMbPid0nmxzZFAYjUY7FPcIt7Mtkg45Ho8L6ZKMdGCq5ERlNQ7pGRgYQKs1OWSIZoxBKpXC4OAgRkZGEAqFkEgkRC/dsiwhpxUKBfT29kqNnPB1OByWYTnRaFQQrrGxsQ4HnEwmEYvFRD45kUhgdHS0A1lgayeDXfIEumXirO37fD709/ej0WgIya1er6NcLguaQARrYGBA2Pd03qlUSrJzBgVavZEEOU74Y5lhvs1x6I4tKOvp6cHWrVsBHB5nro0Zjs4odDZCx64Jesyk2f9NeJ2kIkKRWolMC2Swzsnf+Rl26UzAcX5Hm2mlNjuPg68zY9YTzDjlj0Q3tol5vV7JKuPxOHK5nIi7UCCGCohEePjMcbstW7YIW1zD3dlsFj09PUISZVBLdMkOZ7vdbvT29kq7ZbFYlHozgw6SzOywNp01Vdc4bIWvEcafKROnkJNuV2P5i/wR1vpZ3tLHoDsnGPRwdoMmyrFMRnVJ3TqXSqWcGrpjx56tWLGig9VOoyM7FI5dt2wxm2b2TWizXq+LfjVVtejctWwqANGP1vA2RWu4uOq2MdbgmLHzfcdy+9XRZgz0mFFrNj9LMF6vtyNg5DPJ540kTaqSkc0dDoc7auDUItC1cTpZzVNh7ZnkNj1uNJVKCUQdiUQQCAQ6WrwymQySyaQ4M92yRYetyzxE4rgd69eEzqnZn0qlYIyR7Jms/Ewmg7GxMSGn2YMICs2wts8yAztCMpmMHHsul5PaeyQSkYCBx084Xn839chh3SPPGnq73UY2m5XzHB0dxZIlSxzI3bFjy5hVaPKYtkPh2Lk4cEHVsDjrgIzM2T5VKpVkoeUXn5A6MDVUhjOcAcgCyW0Ir7NmyGNx7Mg2OmoGblp3gMp31EDg86NhdTplsqm5L0quat0D/cyxZMTfybwGJtEdjpfVNXk6QULFrAVTDIYjQgmf27Ppnp4euN1uLF++XDJTKhdSSrVbFk4Y3Ov1YmhoCJlMRkh7LGvZe761iI3f75cWTQY02WxWSG3cB/vJux0D9RzcbjcWLVqEfD6P/v5+CUi0SA3LAayZ877oaXSaOEi+gePQHTvmbPHixR1CMnSqzFDmQ8+9m5HExqieiyuzJsKPzBoIkTJyZwZESUlgSrmNUDu31W1ZjhNf2NaNuMfXtK6AdrRan59OgyUWLvhaXIhOvVariWgLtQSCwSAqlYqotFFbndkwoWQeA2FiEjCByeBxdHRUnCbljXlcrG0Tzt60aROMmZooRtU3ws2EzrXDZmmKQYMmt+mAgbKzDED0uW/fvh1+v18QCupB0MFT+rXZbCKZTCIcDiMUCkkZoVAoiBPmBLVwOIyenh60220ZncvAiogHj4GvzcSMZ6sbxXDIVeA1YiBxRNbQLctyA3gYwE5jzGWWZf0dgMsB1AFsBPBuY0xuz7Z/B+BCANcaY+61LOs4AJsBfNQY8w97trkZwMPGmH+dy+N0bGHbueeeCwAdTl3bfA1psZse4MFaG0UuqtWqSEJyEdWkNb3g6xGeFKXQk58I6zs18YVlM3UM2AM2Oi+yyzXhkXVvMrnp4OlwyQhnlwMA6d/WQR6le+3POxEjXZcm54NiKCRPEk7nz1rHnHAzWx9ZJyZ8rB0ZAKRSKTl2ZrIMcKmNYFdci0aj4vx6enqE0Z5MJgUJ4AjbRqMx7f3xeFx6u5ctW4ZcLidEOD0dTbeO2Z0wa/DkA4yNjQks32q1hOhm//yhoSHpZyd7n6OKOYSJaAJLaZwO5/P5BPmYb3PN8f6uBvCU+v1OAKcZY04H8CyATwKAZVkv2PP3VwD4iNp+FMDVlmUdunlzji1IO/fcc3HFFVfgla98pSiraePCNpu6smaIH4hpcRYqS3GxJrSmgwsu6Fy0NQNeE294XFp33LGFY2wbZHbI7JVtXcBUiYgqgoFAQBZ1Lvrs406lUuKYdT3c7XaLeAv5GczuY7GY/D0YDAojnD3almUJCS4SiSAej4uICYlaVC5st9sYGRnpyNgpCMOgMpPJIJ/PY2RkBKOjo9KmxlIQzzuXy3V0gWiHrHkkhKBJhiuVStImx9p3NpuV7wVr7TqIYABULpeRzWYxPj6OcrkswUE2m0U+nxeGezqdlnJWLBbr0G+gaA5Jf7FYDMCUVgPLHPaWM7YELl68GPF4HH19fejr65P5DplMRmr8rLuT5Md7fihszjJ0y7KWALgUwBcAXAMAxpg71CYPAnjrnp/dANoADAAtGD0G4H4A7wLwnbk6NseOXBsYGMCll14KYBImXL16tUibagbwTKYd+YHW35mlEWYDphTiuBDS0euygG5PIzxq73V3bOEZnZbOsklC4/3lpDduT84FyYwAOsZqMutzuSan4fn9fuRyOQkOSOSyLKtDYY2tUISFmUHqeQIkoTHg5OxxqsSxXKSh4mQyiVwuJ3VjZqc6m/V4PJKJ00H19/eL3KomlhWLRfT09AjKkEgkJNh1uVzivJmBk1XfTbM9Go0KC71UKiEQCGBgYEBIovxHNKEbEsBJhAAkGFi0aBGAqamC1HdoNBrYvXu3XEdqTbCcQFU8AKJcxzo5EQJeNx5zo9HA1q1b5ftOfsB821xC7l8D8HEA0Rn+/h4APwUAY8wTlmWFANwH4GO27W4E8GvLsv5lDo/NsaPAvF4vzjrrLNx1110d4ir7C70fKFTPRUMzztm6wuhfk+CYWbCu6MDqC8d4H+kwuZhriWE6c9a96SxJ2OTrzO78fr/A2az1Es0hHEz+hdZMB6bkWqlKRmTH7XYLQkUH7vF4kEwmZbtQKIRMJiP7sg9BGRgYEAdIxGBiYkK2KRQKQl6jdCrLCuyhBiadGTNYy7KmqcVxWprbPTk/PRgMolgsynhgHh/72tmrzvcbY4QcBwBLlizB2NiYoCLAZJBFAl+lUkE2m512DLqbBIDU2Hk9KQzVTQJ22bJlglhYloVcLofe3l7hD2QyGQAQlTheF5LtYrEYent78cwzz3QgFIdq2po1FwQcy7IuA/B6Y8yHLcu6AMBfG2MuU3//NICzAbzZzPCBe2rovzDGnGZZ1g8wCdf/CfZSQz/99NPNN77xjYM+/rmwUqmESCRyuA/jsNqhuAaM8O2OcV+ZOjCdbDab93Qz/b6Z6qz7+uyj1WZ7PQ6l7Y3EptsTuyE++nX9v66Da9OT7DRUT8SG+6fTYd+2fp4ZRPBz9LHruj1ft/etayIpAwA6aAYKLPVwG7a0MdDQn8P9M+jRx0oHSgY/ABk9yoBIo1p2Z6vfz/3zH7fXn81rrK+znnbIa2C/BzxXe8cK0TV9rXSngH4fgGkSzQzWWFZgsMOuF12aOdD1pptdeOGF64wxZ9tfn6sM/eUA3mBZ1usBBADELMv6kTHmHZZlvQvAZQBePZMz72JfBHALgDV728jn8+GCCy44iMOeO1u9evWCOZbDZft7DZ599lmMjo6iv78fq1atktc5sagbTJXP57FmzZoO+Js2m6xb9/UeSJauFxoAwlDmgq0HiejF+1gxthweDtPKexqK5n0hqqMdhJ77zfukeQ9kZBOJIdSue7VJhON7i8Ui+vr6ZFQt/065VDKoScgi05ukMmOmZrdTgYw1b0qrkgSnHQnPf9myZdi2bZvU9zm21w5NE9ZvNBpYsWKF9FsnEglxUtlsFsPDwx3vY8mArHu2mA0MDMAYgzVr1uCss86Sc7F/NkVfKPeq1eY4hY3kOxJQycbnPSSRj+WMsbExJBIJxGIxGUBEwZtu5x4KhSSI2bhxo+yT7+fvVMiLxWJCvpuYmEA+nxdUoL+/Hz09PRgfH5fn6Mknn8R5550nPetHVA3dGPNJTBHeLsBkhv4Oy7IuBvAJAK80xlT2Y39PW5b1JCYDgbVzcYyOLSzTY1GHh4cBAKtWrUKj0cAvf/lLgfouvvjiji9DPB7Hueeei4ceeqiDlLM/pDdd5z5Q6J1GBjtfZ2Zj386xuTWtUcBFlK/xnjJbsncd0Ikyc+XfNGRNlT/C8uwLj0QiAtMTiiYxSrOq+f5UKiVOgXVzBn9kP2sYnrOztSqb2+0WeJfOmcdOp5JIJMTBUt+crPJ2uz1tsIjL5UJvby+AyZoyW9ZIDKNDZCeHhrSZtZLBn0wmkc1mMTEx0cH0tiyrq/qbDra8Xi+i0WiHEh0h9Wq12tF2RgGYVCol94rXVc8hZ8mD6ID98wOBABKJhLy/v78f9XodfX19woOoVqsd6nOZTAa9vb2IRqPweDwifANAJv9xW12O0cd5KGy++9BvBuAHcOeeRfRBY8wHZ/neLwB4dL4OzLHDa3Y1uMceewyrVq3CnXfeKdF5tVrFXXfdhYsvvrhj24GBAVx00UXYsWMHNm3aJApNs7FuYjX7Y3Y4z6mJHx7rBqsSDtdwKslqJKdpONrr9XZAysDUfGzqjFN6lGqBrJ9zKEcmkxFyGaVJK5UKQqGQ1E7pUDKZDCzLQiqV6mBz66Ci2WwiHo9L+xkdPEd2cswoh4voIIBGx1KpVJBMJiUQYX2c2zJIcLlcEpB4PB6k02k5rmg0KvVqksWIIJBUxoCb0rLNZhO5XA6WZUnNnq/TYWYyGekh5yCUUqkkyECtVptRACafz0sw4/V6BSVjTZ+1fzrXaDTawQ2gSEw4HJZ7ps+Dz4/9s8PhMMLhsDh/vh8Adu3aJdvyHA6HzblDN8asBrB6z88r9+N9WwCcpn7/I+a+rc6xBWBkjJI5zMj6Zz/72bTXOWhB28jICNasWdPxpdGL8kxZ91xJp3Kh0HVGxw69cRFmFq0n1fHvfC7oXKnslUwmOyB0wr3AlD56uVwW1jn3R3U1AEJKoxwq688ABK7X4iNspSK5TgcRGlZnO1W9XkcsFpPhKtqpBQIBDA4OCnubk8cYKOiRpKlUSsoI1DmnYyV0340pTsU1BgO1Wk0mrxFR0Brs3AePj6NPgckgnI6QTrVarUq7WSaTkRnsPT098Hq9AquTuBcIBNDf3y/3kvC3HVLv7e2VsgAdLq+hHhLDLgOeB0sAvI88b3YnVKtVKZvw2Cg8Q8SExLu5rJfvjzlKcY4dcqND19C33Tlrp/zggw9KrX3FihW49957O1juOvPi/vZm+9PDPpPxvZoo5dj8mP36anITzU4SY51ak780K133bjM75LAPZpKRSKSjHZHZOiV86YQ5K5vqaXT8Pp9PskbtiI0xwjSnY6cwERnZRJ3YKkYInDV16q+73W6Mj4/Lz8FgULJROt50Oi290FSEIwelWxZMuJvsdDL0CSvrLJ1BkL3vnNrnhMJ53iyD0OEWCoVpQYjL5ZIgBEBHJs6pbBzD2u397AzQhEU96pTvr1aryOVyHXKufX190vbG+eq8dpVKBWNjY1JmY/nA7/ejp6cHfX19EjRu2rRpvr4OezXHoS9wYwQYCoVEtvFIt/vvvx9AZx/53iLabdu2yf/btm2TDN5ue3PU/BtrgHMRQTMq1/3ljs2f0VHTMWuGdCAQkMl0rMUCUzO0ufAyq6djYy00HA5Ll4bOPJlhkyBFsZJGo9EhcxoMBkWSlMgAW5z4fmDKETGL4/sJbQcCgWnKbIThAaCvrw+BQACFQqEjM/V6vVi2bJlolwOTwi3A1BqiM2GdyRKBYM94f38/AEgWTA0Fftbo6KjomPPcjJnsRQemWsvYb64Z9Jx/TuVFXhMNibfbbaTTaalBU91NlynYSw5AFO7YGtZoNJBOpwUZsXci0Bnbe8jJ2+E6W6vV4PP5kEgkJKhkb3mhUEB/f78QCGmHivw2kzkOfQFbuVzGHXfcIVrAr371q48Kp86sYCan2g0yt2fxM5nO3O2v67ahgzHd0nI4IHetE36oPpufeaiNECYhdc6e1lk3M812uy3ZHO8Rp4xx4SYkb1mW9BsDEAiZvdmcKMZatl788/m8TDjjIB4AkjWTkU1BGNbdqVXAc0gmkyIqw4y82Wxi0aJFoj5IBz86OgpjDHbs2CF1cbugCdvFqCxHdj3/bndeACTzHRwcFMiehD1eN8LT9iw+Go12kL76+vrknDkqlZk+pVQZTPT19clQFZIOi8UiBgYGhExHVELfyx07dkhfOp0vy1+UoHW73ahUKiiXy3JNqG+vFRp1OYG/VyoVuYaUcyXJkfeIqIquoS8Ucxz6ArZf/vKX8nO1WsUDDzyAiy666DAe0eExLfYBzJyJ05nvLXvX7ztQHXgypJlpzTbYmCszprMP+mAdrb2nutv+DldNUBv7hO29xXSgrGUS0vV6vchkMpIxcoBHs9mU7KtWq8nfEonEtNYxMrcbjYbIu1YqFUSjUdRqNclK6QBJANNOBoAw3jlxrFwuo1KpIJVKAZia1MXAs7e3F41GAz09PdPEU0KhEOLxOMbHx4VIxlatXC4nmTVhfAYysVhsGjGOHBUKv0xMTCCTyUhAwhZEMv8JaVPDnVA2B8OQwc9zotmDCaIsbBvU6IYmljHYsWurr1ixokMSmrrzJDvG43HRrNBjXi3LwsDAAGKxmAgBJZNJEc1h8MKggrVzoiAcMcuyy0Izx6EvULvtttumvUbJ04VkTz31FEZGRjAwMDCn++XCpmF5oLsTtmff8z2FjbVAfYyH0vTnaf34g9mfFu3Q15ufpcU3DpXxONgipbkSur+cx8zsXGfR1OD2er0yFU87cZK3OBCEBDVd19aCKnTQzN6YRTPjZTBAjXSdpbOuTYY24XuKkPCZ0qTRbDYr56CvCTkDg4OD4gDr9fo05zU4OIhYLNbB4C6VSujv75cAgv3o9veyTYt18EgkgtHRUSQSCYHA2d5FmVRmwOl0uqM9LxAICCuf10Mz3Xl8evocA6pujPN0Oi0taJyOBkwF62S629vWGPyxBEJeAjXdfT4fRkZGOurqg4OD0jff09Mzo0bGQjDHoS9QI4ynbaEp0T311FNYv349AGB0dFS+FAdrGo7v5qTtmfmBiMQcjNNnhqwXmvkyvYDr4IH/61GNhKQPJMigYA4wlbHzPHU72KE0Hg+vM49RC8IwK9TXRcPqGn6PRCJdnTihWTthzev1olAoIJfLwefzCURO9rrWOGd7FNvLCoWCZIEMxN1uN/r7+8UZMmiwLAuLFy+WoME+qWxoaEh0yHVdG4DMGu82K7xerwtqwVp0OByWrFtfz25O0/5s0zEzwOFENmb1PGZqrcfjcQmqAGDRokVC8rMz3clFyWQyUuYAIMI4iURC9PCj0aiw1YvFIkqlkpQHxsbGhInO95NzQMGniYkJ6c9vNpuizc7j0deB58xsnGI9C9Uch74ALZ1Od32dY0UPhz388MPCND/77EnFwaeffrpjGy30cjCmFdf25si17Y9IzMFm8Fz87NnsXJsdoWC9mAuxHd7sBnXOxuisGRzoIIF/m2+Hru83s2RC7Aws6MRTqZQQo+i0e3p6pOZKp0JoWPdUd3Pi7OWORCKyoLO9zO/3S6BKmJ3Qrp5fDky1yfFnYAoKZhBAB0JyHMlqOgjTpDj2V7OcxGdXE9XIJQiFQohGo6hUKmg0GjJCtdFooFarybWiOpvd8VEgh2UGzl1nex8RCQYJhPm189dQea1W64CmGXzpsgihfDvbvLe3F6lUqsN5Dg4OIpfLoa+vTxALcgd0EFQsFjE0NCQiMNlsFv39/SLKw9IEnTjnnvN542AZ3sOFCK3PZI5DX4CWy+W6vv7www/jda973Zx9zsaNGzE2Noa+vj6ccMIJM2533333iZobM4NkMjlN4tMYg1tuuaVjpOTSpUtRr9exc+dO5HI5iYrt1s3J7m972cHWxmdjWk2LtbX5Mg1x6/5qYOpcdSY918fCGvWhaM3TrWZsv2Lmq5nr9XpdBE3Y51wqlVAqlWQfiURC4HLdVkXFNZ/PJy1jmiFdrVZltKdlWejp6ZH7rXUOeB20uhnhXZLZcrkcotGo1PYJvYfD4Q5BIu3EBgYGxJnUajVUKhVEIhEpZzHjtvddDw4OoqenRzLhwcHBjhnlmvDHASP8zKGhIcRisY7gJJfLyXs1sU07zZGRESSTSeENsLzA3nxKpbKdj0QzKtmxfVDLs2p0oF6vCzLCDJnDbmgkQ9qRCV5rXTqhKAyvo1Z2I6yuGevUBFio0PpM5jj0BWhallAb63cHYo1GA7/73e9kelAgEMCWLVsATLWFdXPqjUYDw8PDHWIvmzZtmrHPklF7rVZDqVTC6OhoVwe+vzD5TApvXFxnmykfrNNnpkonM5+Zqya/cfHWMLT9ePbHZlsPn28ZW40+8HedpbMljbVpZtrGTI7m1PBvKBRCPp+X7JzZKt/P+jpbF7kdM3US5TTRkbPFOV9cO3CPxyM1YdbGObGMziORSCCRSGB4eBjhcFgY193GbnJSG8eFjo6OSlZMERjeE/0+3ZbHTJqkOwZihKvtjrNcLkvwzWyave7AJErDoMA+XpWyr7yHQ0NDMkedgYmus3cTkGFAQA31YrEos+P1WqJbFcfHxyXbBiA1dLbEsQQRCoUkoPP7/bOG1YGFD63PZI5DX4CWSqUksp8LazQauPXWW+X33bt3T9tm3bp1eP755yWD5pfULgIDdLLJdfbGiFibndRG219nPpPty5FrGLfb6/trhKY1834+nB1hVC7ahI01AfBggolDXQ+fybTSGzB5PdkKRecIQDJvws9sOdKjN9mixsCHkC9bxQh/U0CGzG6SzsrlMorFogiHcLAKmdEMMLg/Etaq1SpSqZSIlQCTGZ7P55PaLduoxsbGkEqlEIlEEIvFJKOPRqNIJBIyu5vkOP0ewuAApGYej8c7nB9r88w+uzlOOndC8pwNzhY5toxpJTZC5LpfnMQ2ksz8fr9ouPNeMjPvpmLn8Xg6kIXFixfLoBYGRJFIRNATrWiXz+cxNDSEeDwOt9stUDzr4YTxCbFTSOhogNX3Zo5DX6BGR3qwGWU+n8cDDzww7XXttLm4z1S7pxF21E50X8Iw9jrjwZh2zPtzPeZK8hWYulZ0svNl9sxbB05zAbEfbodOGJv3lLVu1niBKQZyJBIRJTdgquwBoEPIBZhqKdQtaWw/qlQq6Onpke3piDgXu9lsdkw9A9DRd+3xeEThbGJiQqataefk8/mwc+dOHHfccfB6vSgWi9PmfnM/LpdL1MU4GITnaHd+xhgEAgEJMAYHB5HNZtHb2yvaFKFQCMVisaNful6vw+/3d2TSfX19qFQqGBgYgNfrlYCjVCohnU4LDJ1MJuX6BwIBeebq9Tp6e3uFx0CWO9GEdruNarUqkDjXDHtA0Gw2uwrI8P7zs7rB8kRG2OWgtdZ53zSsns/njxpYfW/mOPQjzFbvx4jSdDqN1atXdzg0QondnPJMti/nqfc/m20PtsY9l1D9gZgmkB3ssJe9fYa9dq1LF4eifWy+2tToyNlGxLq5hvf7+/s7iIcakqdD1aNJ2R7FhZ9cD7tOutvtxsDAgMDXXq8XqVQKu3fvFsY6IWg6XwAd4i/lchl9fX2S6RIVyGQyco+oKmeMmSbHSiUzksQotUpYmvyAUqkkDsflcgkbnKIxHB3KMoJ2/NpxMkMNh8PCU9CtYjqw0oFAu92WkgONg0zoFFmqYP92JpNBo9HAxo0bZUypHlWqAwIAgiow6CJCODIyIoEryxAaVrcsS2B1BkGE1Wu12lENq+/NHIe+QI1faBod4ejoKNLptAhSdLORkRE8+eSTGBsb69rPfSBtXnsz7ms28PdcfdZ8v2df1q1Hfq6Nmcd8M8wPlWl4k4xrZnRclFn3Zp1bs9BZx6bamGVZ0gJGMRfLmhpGomvn1WpVnCAAUYSjU2BWTz12Zo1UJNPOQRPyGo0G+vr6YFkWRkdHkclk0Gq1sGvXLmFaE3Vwu93o6ekBgGlDUfL5fIeOOFn0JKj5/X7JVrsNVCGMz+sETD4/PT09Up9nuxYZ6IT7WWJglszAgwxxEtJ47vpZZKskMCUeo9Ejvo/GejyRBhITa7VaVznbJUuWCJEwk8kIW533kERbng+DQLfbjXg83iGYdLTA6nszx6EvUDvjjDNw//33d80An3nmGbzsZS/reI3ktSeeeEJaWoDpkPiBOjf7Mcy3eMtCNGaHuu95vpTimK3MZybOxbnbMzZf2Tmfm0AgIAs6YWFmqszWdb84HTYz3lAoJOM6dT2dkpyaTEUnpmu+OhMFpqB7yqYyG3e5XFi8eHFHBsnjZBDC+qzdoTFYYL+6loqdidlNYRYGJHTuOismpK3fGwgEEI1GO64z69qsbfM60mFms1kJBDjGNBwOY3x8XLLriYmJrpPJSqWSdBHwGYrFYtL+RdRCs9xJXiPCRdY7gw2gU1GOARp5C81mU2a8A5MZth1WX7RoEeLxuAQfRFqOJlh9b+Y49AVqixcvnsZwpnG6Uj6fx/DwMHbs2IFsNtsByWpHbh9OsL82U2bdDWrfW437SA8CmDXZ69ms6dI5HIhp2HM+WtC00UEww+x2f+famfP8eF4MggiZ659TqZQ4SGaRdNhaWx2YXk/Xn6GvITNwAOJ4+DwyKGg0Gl2lWMmI7+vrQygUkmBBE88GBweRSqWExMfjZdbJf5oEaO/7ZiCix6U2m02RWCWaQaa/fbQnWfb6GeW50FHagwhdzwemhFgoyEMiXrde8W6oADDZNbN8+XJ4vV6BvykAQzW8sbGxjntFNESjMo1GQ/7OVsV9wep6ohrP51hw5DTHoS9gu+CCC3D33Xd3hXWfeOIJPPHEE12d+Fw5zn0Judh/1wS4+YajD7XZxVy4OLKmyfrsgZyzZq1rxzdf2Tn3Ox+1/27Gc6IGtr6OHB4CTLZ4FYtFVCoVcSh0Ygx2SHSjtroWg9FDc5jh8rMty5IggLKtrOuS3NZqtZBMJqWVy+PxoFarYefOnRIEJZNJcUoagqcueTKZxMjICJYsWQK32410Oj0NGk+lUnLOLpcLmUxGyG2EjyuVCgqFwrR+c2afZHL7/X5Eo1HpM6/VaoJIMIjQgUAkEkGxWOyYTEbUgpmvNt061k1Jzp4w8JryPlNNjhyEarXaEQgtXrwY0WhUdO85NY5lAiIR+Xxe7nW1Wu2A1RkEaxTjWDXHoS9gI3vV7hwnJibEmevXD8aRz+XifqRn4t2MDpfOkBAtMzudFdq7B2ZrZP4eTKa/0EwHKESctGIYfwfQkXWT6Mb6OsVmAEgvN5nmZMlTMKbVanVA0ZyM5vf7kc/nZdiJZr3TcfD+DQ0NyX7z+bzsz+ebnA2eyWTk+LXACQMV1rrt0DjHrgYCAckmSfgCIPvrlk0zeASmAmhO/CLpyw6rW5aFRYsWiQwtrymZ9QBkHnylUpFebzLXeSx0nmxpo8obxWoI6xNdKpVKsn5xbjoJbzwnj8cjTHW2l7GLgP30Y2NjHSNLWQpgW+CxCKvvzRyHvoCNk5zspuvjGuqeqc49185a14w1KjDXZLuFYkQedMbI8+YCpclC+9tWRkdyOMaTzqfpFjRg6nkMhUKiKEZ5VT1ExefzTRuWQj1yPufMyrjfdrst88fdbjcGBwelLqvZ3nYmd6vVEqcPQJjmJKHaiWI6I1y6dKlk/lqBjVrrdHpUUAuHw6KwSOEXXqN0Oi1wPDCZ6cbj8Y4+fMLOVGUDpgR5GODYM2kGQxxGAkC6Cuj8EomEHLM9GBgaGpLyAS2fz2NgYKCD96A/u9VqYfv27dIDz3uqAzmWVPhMlMtl6VsHIP3wrNFrMqJ9ZsSxBqvvzRyHvoDt/vvvn5aF23/fH1gcmF32bP8MGrMG2tHqwLWxr7mbMpzuCbdzB44257w/puVSdRmCOgacuc0Rl6yfM1iiY2ZtloIuAGQxZ32dzpaDi/QULY/H01ValUIl4XAYPT09MMaIgEo3SdVkMolCoSABAYl7VCuzZ54k7vFYFy1ahLGxMfT09MhxMgvWgYAdgq7Vakgmk6Iwx3POZDIChbMTht9ZZusk9RF9ACAjQu2BJrXeNbFPt67pUpqusfNes45OIRf7dUgkEtJm5na7sWzZMpHJpS48yyBUxCPznqp1zOa7lQUcmzLHoS9QS6fT07KbbhKqQPd69Wwcrb0FhTaT4hsAIa0cjabr2LwGlmXJIsV7oBng88UGP5KNWvckf7E0QeW7RCKBXC4n+usUHWFPNwAZZsIWJULAlmWhr69P+rI1/wCA1IrtdW62fFFWlVO5+B4A0zJbDefzeNhrHggEhIhqF02h8wIm0QjtBDWhq5tYCp8zHg/HsPLa+f1+aV1loNLf398h/lKv15FMJhEMBmUAS7lcRjabhd/vR61Wk+tEJ04eAIMBthWyTZC1cF1O0up0LJdoNT12KfD+VCoVmfgGQEoglIfmMVmWJSUJkgcdSH125jj0BWp33XVXR6Zsb0Oj7U+GvGjRIixfvlzYusAkfLZlyxbs2LED9957L4wxOP3006ftl1F6twj/SDftEDTJRw8m0bXNQ8FEPxKNTptOqV6vy5hKPcZUT+tKJpMoFotIp9Md+uH27FDL4Gqdco7prFQqwmdg6xLV4ciS51hUQuQ6QCCUSwcVjUaRSqUEQmePs874tSMFJsl2fr8f4+PjiMfjwgwnN2J0dFREYpjhs0WOpSxm/5FIRJ4zOl7W5en4dRatx5QykGJGDkw6TzvRjn3d0WhU3q+DAaIJ1WpVauUk2pGHwHsdDofh9Xqlva2vr0+6F7gPziLnMTJY5vF4PB4sXrxYevxpDqQ+e3Mc+gK0DRs2THttb448HA4jHo/LsAoON0gmkyiVSnC7J+cwh0IhPP/887jrrrtQrVZx5513ymhBQlrxeBxjY2MywMVuR6MT01mHztI1vN6t9c/JzqdMO3Mt5EHInA6GEHEkEpFBKmzRouqX1+uVvu9CoYAlS5aIY9YOIRgMSoDF31kXp5OjY+3r65M6tJ2oRiTG6/Vi8eLFGB8fR09PT8cAkG5wNDNJPgfst2bGHgwGkc/nO5jdHo8HS5cuFbGUUqmEVCqFWCwmrPNyudzRl55OpwUx46hYBhgUaiFTnUx53hMee7f6Os9BBwMAhGxHJ0qoe2xsbNqkNl1fZ2bP9SaXy6FQKCCRSHR0LpRKpWkQP1EI3Ybo2P6b49AXmG3btg0PPPCAQFW0mTLxQCCACy+8UBY6bY8++ij+53/+B4lEAj/4wQ+QSCSwfft2nHbaaVJPXLZsGTweD/r6+mSMol0ohc6NC9dMUP1cWjdnuS8H2s0xd9umGxPdfn7dBtAQEuTvOps/VC1gC83oyAFIHZyOiXVgYBLW5gxq1tLpfHw+n7Rc2Z0jyVrs/+b0L2aMHElKBxEOh+XvOhvdunWrDCmJxWIoFouiQOfxeGQEKwOQbDYreuoAhEVOERmeI+FoPWqULG/K09phdY288Ri0BrnH4+lodysUClJbt8+EpypeqVRCuVyWY6YgDMVfeH/YGmg/BzLH7WRXto91CwZ0bz+DBta4Xa7JEarZbFYgdcre8jzD4fC0yZLHetvZwdqsHbplWW4ADwPYaYy5zLKsJICfAjgOwBYAVxhjsl3etwVAEUALQNMYc/ae14cA/GjP3/4fY0zJsqzrAXwcwHHGmNE925WMMZEDPcEjye644w5ceeWV+OIXv9jxur2m5vF4EIvF0NfXh4GBARhjBDbfunUr1qxZg3vvvRe5XA6vfOUrEQ6H8fKXvxyJRALValUWhsHBQQkcCI9pmAzoRAbmcpSmrn3O5HjtNW2+3s152mupfI2/8/2cuOV2u6X1jPVWknMIcWrhC94DfbyERGdy5gfavnYkmK6L0+hICZkTVo9EIvLsUKpVs9k1P4Q91RMTE6Jg1m63heTFEae6/MFsmgNH3G631IztDohDStgeyKEoJKnpbJqONBaLyf0na5tEL8L/OnigjCmDC6/Xi0Kh0DEIhVmwzrJncp6E2DU5jSgAs+h4PI5arTYtEOjt7ZWyAN9PHXVei3K5LOfBYIAdNgw+9Hnw3lerVYHHNWeAxmBmJtIfbXBw0KmRz5HtT4Z+NYCnALBn4G8A3GWMudGyrL/Z8/snZnjvhcaYcdtrHwVwFYDjAbwDwLf2vD4O4Nq97OuotNtuuw1vetObcOaZZ4pzAqbIcK1WC3fccQfWr1+PG2+8EbfffruQRTjGkQuV3+/Ha1/7Wsl+uOCRpMJsoFurlP33bs6IdXQuqMDs5pHrTF87SjvZj8xoLkBaLU8LfjCbo8KYdsTaIfNn/TlEIrgtnQ9beggdUomLGacOfBhccFv7eTBQOhqN87+NMUJk40JPGFdDyKyzkmHOei2JccYYGZVZLpclG7Yrt9XrdamXU6udrV2cp60h+G7TvZhFmj1dG6yba56EnhfOHms6VRLdAAjKoKejsbWMUqwMNKgZT+EYDl9hXzgDRD0Fjtd0YmIC+Xxe9qlLG6zz25XT2PJF0Rsa2/F0MNBoNDA+Pj6NbKfr65zups8rn88jFosJ4ZD3gWUCHg/7xtl3rs2pkc+dzcqhW5a1BMClAL4A4Jo9L78RwAV7fv4+gNXYPyfsBtDe8097g38B8L8ty/qSMSazH/s7Yq3VauEv//IvsXjxYnzoQx/qaM2gkxgZGUEikcAb3/hG3HXXXQAmW3N8Ph9OPPFEgbIoIcnJUZRDpAPWGaXOau3MbQAdv3Px1j3Aejvt0PU+gc72Lio8aceoj4XkJ80w52JEp8mWGDJ2Q6GQzHPmwkJyEVnplMPkcTOLZNbk9/tRrVY79L0ZRHTrP9fXi47BbrMJco5U09wCY0xHDZgzqVmjTSQSaLVaHaMxfT4fisWiQOatVguxWEyezWAwKBA67yH/T6VSQtCiEAozeNZ2iUKxJh+Px2HMpLpaPp/v0DOnWAnr4QwKKEZTKpVQrValdYyOUzPENbvb/uw3m01hnAOQYx4eHpZnbGBgoEP8pVqtIplMSnbNwIicg4mJCfj9fpTLZSGcuVyTkrAcNUulPPuAFW26TcxOttPqcZRijcVikoVz+tnIyIiw1BkYMTDhdac5kPr82mwz9K9hEgqPqtcGjDG7AMAYs8uyrP4Z3msA3GFZlgHwbWPMP+15/WYAPwSQB/B2tX0Jk079agCfm+XxHdF200034fTTT8cJJ5wgC4TO7trtNmKxGE4//XTE4/GO93KB5IKi68daGYpfKmah9r5qOlq9HU07LpLnmLF2c1qEsLXTs5NddHarWeVseeLvfC+dLI+Tx6Nhe32ezOYty5KFmMfGth3+0+/h73qf2lnva2AKr8eh4BkcDuP90rVluzOjWIsWNpqYmJDnhb3mGuUIBoMCt9dqNektZ8C1ZMkSAJNdGcwK+bl2R8Qecz7r3A/Rlm4KbBoGJiqQTCZFj7xYLHZomdMJM4PN5/Po7+9HIBDAyMiIMOgZSPI82LJWLpe7ZtJ8fljjjkajsCxLWtYymQx27NghwUs35TStPc/AXg9YIcph5whQyEar3ZG8x64Cu0CNHVIfGhqS2fA0B1I/dLZPh25Z1mUARo0x6yzLuuAAPuPlxpjhPQ7/TsuynjbGrDHGbAXwihne8w0Aj1mW9fcH8HlHlDUaDTz//PNYuXIlTjjhhI5WE922oyNb7TQ4U1o7EO2UuYjSueoaMB20HXrn4sfPIaytnZwOHuh4uYBxgIIOELQT4HuZaes6N8+di4UOSijtSXiSkDnPhWIU1FXncbJmC3SS15ipNxqNDuduh9V5zHQG+3LWvGcz6QYcicayAuF0TkijE2Hpgs8kSWK8H+wr5vsAyDaEXBuNhrSLjY2NwbIs5PN5gbyr1WqHYAqHg5ANb1mWfH/y+XzHEJdyuSz3kxPc6GQIEXOKm0aRqtVqV9lSt9stjpTytNSW53eAUw/twa1mqhMdYisYAxWqqGkH3y2LJvlMB8tut7vDefI47UpwxWKxo6ZdLpeRTCYRjUalll4ul4V3EgqF5Ptk76PXKn92cyD1Q2fWvrIIy7L+D4A/B9AEEMBkDf2/AJwD4II92fkggNXGmJP2sa/rAZSMMTft6++WZX0RQAHAdWYGUtzpp59uvvGNb+z1+A+VlUol6dvcH2u328jn8yK6wBqktr1Bt/asvNvftWmnare9/a3bPmf63L2x4bu9V38ef9eOlK/bA5aZjl+/T7dS2Vnq3Oe+PkMfy75q4t3ux9GQrdsdhkZKeL8Iu2qluG4lHjp/YPp1phNjps9rqfkQvJd0rIS1ifDwWLQsr+ZGAJASCrN8HdDZP1cHxDpI4/eV7+G2XAsoXct92evHemQrj4nXV0u76ucS6Jw7rq8z75G+FrzGOkjVsrlaH4DBqlY71CihDtB5L4mQ0LxeL8rl8gGthUeTHag/mK1deOGF68wegrm2fWboxphPAvgkAOzJ0P/aGPMOy7L+DsC7ANy45//b7O+1LCsMwGWMKe75+bUA/naWx/wVAH/Y2zH6fD5ccMEFs9zd/Nrq1asP6FgajQY2b94sY1CZ0TDb0TrTdrN/wewsdN0CxIzZDnXrDJXvtct1Ap31dBqhRNZKuS2HNpD4pD9XZx9U3wLQQQri53IB5KLBbfmzPkZ7Nsw6POuPrK/rQQ6so7NmCXQy17mw+v3+aaQ4Lqh2hrtGN450J64hY54vlbvcbrfcG7KjqVm+ZMmSjg4CktwIxbIGTNPTzvr7+6UNi0NRAKC/vx9+vx+bN2+WZyQejyMWiwl5LZlMdsDaeqhKPB5HPB4XQluxWEStVkNvby+MMVKjLhaL8j5jDAYHBxGLxWBZFtLpNIDJZ4vPSywWE1Z+KpVCo9HAfffdh1NOOUWIbGSqL1myRNq0Go2GZOSEsVut1jQ5WI1EUDktk8mgWq2iv79frh+/Z+FwWEiyRED4TAeDQcnQgcnv7/HHH49gMIiJiQkUCgXE43FZb9rtNnbv3j1NDjcajUo5i+eiIfUDXQuPJjtc1+Bg+tBvBPAzy7LeC2AbgD8FpB3tu8aY1wMYAHDrnhvvAfBjY8zts9m5MWbcsqxbAfzVQRzjgjev14slS5bAsianIvHLAwDDw8Mi19jNmAXzH6N1XQPmdjqS1wQxQqSs4zGK105dO37CzqwJ2tu7gKlMl5kaoX0AHf/bIUy2kxFK1D3e2qFrRjuAacEHj4EBhWVZ0g9NoiCZxXTUJONVKhWBdxkE6AxNX1/9mXbnfTQ4cz3iUwu6aGfOmd3AZCBXLpdRKBTkPV6vVyBzksg0bN7f3w/LslAsFuXZoygM74PP5xMGPPXQA4GAtIfxO5PP55FKpaYNRWF9vl6vo1arSV28Xq8jl8vJc8wAmsx4Bp0MGEmUI/E0n893zDsfHR2V593OVGfgw/PSnA5eV7vUMIOcWCwmhDreG82018NVMpnMtDo/yXb8Xi1duhQAhOzH55xQPTkPdtY8J8PZYXUHUl84tl8O3RizGpNsdhhj0gBe3WWbYQCv3/PzJgBn7Mf+r7f9fg2mWPVHrYVCIRx33HGo1+syC7jZbGJ8fBzlcnma+AKADidDB6cJZQA6iGR2MpkdymONXEPGzE65CJJhTsdGR87AgBmZPdvnMbBGThhSZxHcnsEFW8Z0Rs9tuD8eJ49HO11d/9aQKLN1CoDwvJh9MzPj++ztaRqGPBrFZFyuSbUutqUFAgEJfniPqehWLBblGaDDYWBE5joncVFOlc+ddhy6/j46OopFixYhmUwimUwil8sJglStVhEKhcQ52+valJrlUJTBwUGMjo4ikUggFosJC1uT2wqFgowX1cIvdICcja4FUWh8jrgvBuck4hE1YAdAtVqVaW7ManXbHK+BJvsBEKerv0t2B8rvk74eXq9XAgvux+VyifaEz+dDPp9HOp1GpVJBKBSSITQ812Qy2fE5Dkt9YZujFLdAzB7l5nI5jI+Pz1jP5sLLRYCOmE7L7/eLE+diDEx98XVNGYAsZtYe0gtHHXLMIhd5CnTwmHWfNhfreDze0dMbiUQkwyFxjY6UC5CeQkXnCnQy7DXrndKg+u/awepsmteIf+dkLWbuzAh10MNzYkcA93GkZ94zGc+TjpytZ8BUnZXywiynUPaTWSSNIi18ppPJpEDSvb29Ivnq9/sRDodRLpdFdYwMbH4us34OKwoEAkin0/LsRiKRjl5zwv9s0aIjNMaI45xJwhWYmjPe09ODZrOJfD6PXC4nQQ1Z9GSI633Z2zcZYPDZppCORhVyuZww1antns/nRc+dAjiZTEbkZBlot9tt6WknrK47DjTKwv2QtQ5A+tPZPlcoFKbB+9zOYakfGeY49AVqlUoFxWKx67hALbxiZ5zz77oGrIkvwHSiGBc8Om86SdY6mS0wENDOn1mEDhiKxSICgQBCoZCwjZmVcOEjKYnOOZvNdmTlzBII1xszOUqRbGU7A90O+fN/Tb7S9Xotv6mvITM/WjduwtFi+vrx+SB5jA7O5XKJk6HTYImC3AI+i9xXMplErVZDKpWSNqhgMAiXyyX1W3ImIpGIwOnMRHWdmapp9ow6Ho+Ls2IZhXPCWT+nDkE6nUaxWOzgADSbTWktq9VqItqi0R62co2Pj3c4YML6gUAA+XxeZrtz6AuADuY7jc+fHVLXM75Zz6eD5/uKxSJKpZIEL7rnm8qRrK8vWbJEnnE+s5ylrsVfOHRFl72i0eg0MpcDqR855jj0BWiNRgPbt29HrVYT4gsdjiajsH5Jp8N2MQ2xs4bJWjEdvXZcXKj5GaxPc2Hj+3T7Gr/8ds1zOn0urPwsHgM/hxmhdty6DUZn21SFY5uP/jy+d3+vL9+nWc/2UoR25keT8doRhWEGRmdORMbtdovGOq8FW7OYsbGsUygUhDzWarXg8/mQSqU6Bq+wn5tQPksrsVgMkUgEu3fvRjQaRSwWEwGXbupnmlEOQJ5xck1I/mJNmfdy+fLlonLm9/sxMjIicDiHnZCrwT5x7RjJaqfTY9lnaGhIAuFqtSoOmCRMZvTdRozag1F+lnagbrcb8XgcmUwGo6OjHQS1RCIhgQ0H37B3PZvNolgsIpFICOoGQI6NQTIJi1SQc+zINcehL0Arl8vYvn27yJzS7Kx1nTGRFc7Fl/909k3yDAAhgAFTsDYzGk38IvzNBZh1OsKHjOwJ+/H4uDDaRWO04ho/h5/P9+g+YP137cDnKlvmfjThzv63o8E0F4HIjl2+lBA1gzb2WLfbbalbsyOA10YP6OCcgWXLlgnDmwEf7znnhhM5YvbN97LGTKdFWJ7qbUR4MpkM/H6/kN84qIV1eZZJdCCQy+WEA+D3+yUTDYfD8Pl8HdPEtIKcdsC615rBqNZUJxS+fft2Qas4FIZZPbNlBkYaIaDT5/7IoCe0rs+HOgC8vtVqFX6/X8pG0WgUuVxumpIbANEFoIqeA6kfHeY49AVo4+PjePbZZ0XIgYsbM1tmA6z3AZOLsm4PI1TOxciyLFGK0tkona6d8Mb6uCaZMdOqVqvitLmN3o/WtaYsJRnszNRZy2TLjEYBCF/ayWlHk4M9FEaEhdeNCzadOWuolF9loETHzVnjhL45UCSfzwtrm7wKHSDwedHZ8dDQEBKJhDwLgUBASHLlcrkjm+ZQFM40ByDtZalUCm63W+Dicrks7xsbG5NZ3uFwWGrqlKNlOyWhfZfLJYNINDrEIJVEuXq9jqVLl3YQ2LSwixaC4j64H8LnmthKoR2iBRyOQk5JvV5HsViU/bEWHovFpEOAx86ghNwUXaLjvbYruelhM9zOceRHhzkO/SBtdHQU1113Hf7sz/5sTvbXarXwiU98Ai960YtktrKGgumgqbrFGrTOspndJpPJjglt+XxeZlLrujEzbMKZZHqT4cwsnwsYFw/dCqezetYlmXGzTECInQSjdrstMDqDAl4DuznOfP+M98ReotHKg9VqVfgVsVhMoGqqoFEXX3MegKmsnLVxZqJ0SpQJ1U6tVCp1KOjxdQZ29myao0yDwWBHGYCOh8dln23u8/lkEhowyXTP5XLo6+uTdjE6daIDY2NjcLlc6Ovrkzp4MBgUrXTC7Aw+eJ7lcllGtzKLBqZEVzSkrks4zWazo7bODF/X6vP5fIe0rMvlQi6Xk+lp/C6Rnc4gm0E5/8Z2VLadMVhx7Og0x6EfoBlj8OY3vxn//d//DQC48sor52S/t9xyi2Te/CLb27n4PwAhw5BAQxidETz7WulE8/k8PB4P4vG4ZCp68tjExESHHCwzF0KxGtrv1srGDF2Lq+ja5NFYk15oxgCM94XOj1miLnVQjEUHbLrzgQQ4lly4r1KpJLA126EI/1J+lV0anDxGtjqz8GQyKZkqGfQkiXEcKoeiECUgqsTSEVuteOxkhHNGOREg3SGhUSJNTvP5fFi2bBnGxsbk+0VmON9PB05meD6fRzKZlCFIDEAooUoUK5fLCTeA1xaAcAXsI1PJ8NfDUer1OqLRqHzXc7kcyuWy7CccDnfMYbdrquvPdezoNMehH6D99Kc/FWe+detWbNq0aU72++ijj6K/v18WGGbAbCeKRqOoVCqyaHPqFKFttuswWwgGgzKIIRKJyMJRLpeFKMOFjZ+j22VCoZDUNen09WxxLjQ6A+Fna9PykI4dmOl2Q81x4P8A5Lmh8202m6LsxfGlRDuI5LBOnM/nOwbakEm+ZcsWQVoACIJTq9UwMjKC/v5+ITUCEHieNWjWue0DSVwuFxKJRAesbp9Rbp81nk6n0dfXJ+IyRKSKxSIWLVok34lSqSRCLhSz0WJIVE5jPZ/vM8ZIK93u3bvh8XhEbZBCRPbzCAQC6Onp6QiUwuGwQOrVahWlUknms1O0hdPTAMg4U2b2vM40MvI1pB6LxZDJZJDL5QSuZyLAYS0AJBhzauRHvzkOfZaWTqfx5S9/Gc888wxuv/12cVqf+cxnsGzZsjl16Geeeaa0b3HBJoTIoQitVktq7GSmklxDGBSYiux1fyzVntrttsCsHo9H5CxzuZwQ5LjYk9XLrJ9sXraicUylA43PjxEZIZqiSZFERkh2I+rC38lziMfjgqhoqdtqtSrjRjUSRGfG4R8c2lEoFDq25VAVsrNbrRYymYxknIVCQTgYzKJ5bOyv1kxr1n75nM/UN84ghg6Xzp2sfQq5zKScxrIQx6TyWmixGN0Kyv/JI9BiO9xOO2gae7o16Y7zxhmYMEBavHhxRz2f+7QjBURUiEjw+ixbtmxaq6tTIz92zHHoe7FGo4F8Po8LLrgATzzxRMffEokErr32Wlx33XVz9nnGGGzbtg3Lly8X0Q4SlwgTkmRDWJs921zA2Y7EDI4kMwYgjP45tYw1vXg8LrV1LkrtdltqqoTq2X9M+JNQqGMHb7qLgSQzZrKa4Mb/SSTk9nqELtuRyKZmqYOZL7enUAnJb3TORHPY+2xMZw82YW5+DtnwLLFox6Xbq3w+n/SSAxD9dLbH6bZFniPheGaZ3D+1DjQ6QdSBGasOBqgSyF5sclLC4bDsIxgMYteuXdPKXGxvY8mL5QIAoprHawZMZdf8ftgFbRgYsJZO6Jy1eracJZNJ4STwXjHD10gNSydOSevYNsehz2CVSgW//e1v8b3vfQ/Lly9HpVLB+973Pvzpn/4pVq5cOaOC28HYP/zDP2BsbEzgbf0ZXOD1RCSSnPQCyIjd3krERZ6ZHBnKXKDK5bLUQvXoVLKg6cA5TpLBhJORz50xC2eHAbM8ZuWEj9lyyNfZumiM6RicEg6HBW7mfeRwFL7OgIHiJNVqFb29vUKUtGfHPp9P2rAymYw4eELWAKY5UT6rwGSGq+Fgn29yJK7WIC8UChgaGpKgFoC0pg0MDAhCxH+6/Y3ZPdEC+4zvUqkkgTKfXT7rJOPxNWbgXq+3Yw4429dYLiC8r1XudF0cgOig85rr9jfC8LpNNRKJIJ/PI5PJoFwuS+lLBwwkzbKMwNcdO3bNcehdrFwu449//COefPJJvOQlL0EymcTf//3f44QTTphX6Oof/uEfsHTpUoyOjuI1r3mN1NWYOfFLW6lUpC7GQSk8LrbqMHPjgqJnfgNTzoMLLnuR7Q5aO3G7wAZriY4dnNEBa/U1zUgmRM7Fm/dOIzFUSaNzI1TNfbCOHQgE0NfXh2w2KySuWq3W0f61detWcdqUVuVIzlqthnQ6LchRsViU7ggSI6PRqMDsbG8DIHKyGp7uRgir1+sdLVnAlIgSUYhoNIpQKIRSqdRx7C6XSyaksbWS2TSHvbAspeHrfD4vrZUkkvE5Z92cZDQSQvl+jkrduXOnIBT8DM1ST6fTGBgYkOvI7y6vnx7AwvcQLaAsazQalfIEbXBw0KmROwbAcejTrF6v48knn8Rjjz2GVCoFl8uF/v5+IYnN1xeG9byhoSEMDAygr69PamSEzrmQcwEnqYb1a8Lr7AenyASHWlAzW6vL0aEXCgUAUy1jeqwqxSgI1+fzeYGBHds/0xr6vAfsM6bTZB2Z/AQ9slRP42O7mW5NAiDs6ImJCclQGTRwKhh7wbkvu0Nl7zQhXr/fj1gsJtmkXYWNsDqRAU2IYysY6+o8BgAd4jEMROlsqZtgrx0DEBjbHmiSFEiHx8CWDj4QCMiUMX0shN6paMdOEAYthMT5GaFQSIiAvB48Bn4/SYoDpkSTotEogsEgcrmctOaRhc8AjfucjSwr4NTIHZsyx6HbLJ1O4/vf/z5WrlyJQCAgXyAuuvNh7XYbF198MVatWoW+vj6ce+65Ug+kw2X2QTicte14PC5CH/aaKmFZndmwVgpMOW0uKhTeCAaDKJVKkjnMdMyO7b8xE+f9JIkMgDgVBldkoOs2wEAgIPeQ95qZJGe9swTjdrtRKpUEok6lUgCAsbGxjqyPmSThYMuyZEIaYWCSI/lM2AMAPRGP9flIJCLHZ++zLhQKWLJkiRDuAKCvrw+1Wg19fX1C0qxUKh0sbr6f149dHfqa6fnrRKh4bFSK08Q3l8slc8EBdJQhiAYAk8gd1e9YtiJiYifK2b8fJJOyjBKNRpHJZJBOp1Gr1SSg0VPXIpGIiNjQ4Tvm2N7smHPozWYTH/3oR/H6178eF1xwAXK5HJYsWQJgMop+4IEHsGjRImHRer1exGIxDAwMzEsUvHnzZlx55ZXIZDI48cQTkUqlsGzZso4aOhdnZt9kKwNTrHcy0pkhRCIRZLNZGbDCunehUOgYkqKJTMxa2Oc7kzN3bP+NARfZ43SKdDQMqpjpaTi52Wyiv79f1PMo9qMnw+nJe3ryGeV5WYOnKJDOxBkgEJrmdszsm82mDM/hecTjcZFjJUM8k8lI2yJJYMw2GYDYxWNYWyb6xGsBQCRa7YHA4OCgtOIxuE2lUnLe1WoV2WwWsVhMvjsAhPnO49N1fn6fyENhyyjvXTQanaZmx7Gr5Dbw+hEBqdVqknUT5bD3w2/btq2jL59DWbRuvAOnOzZbO2Yc+g033ICbb74Z2WxWoOPLL78cK1euxHPPPQdgsjbd09Mji2a1WsWiRYtkutNc2k9/+lP86Ec/wv33348/+ZM/QSqVwgtf+EJcdNFFAolzkWEmRwZzOp2W+icj92w2KzVJXV8jU5i63MCUUhfbh7SmO0k+jjM/OCOqQrSEv9OYfbEuzoCNbHY6ACr71Wo1CdaAqSyYioGsvxIetiwL/f39sCxLfgcgZDFC/NVqFeFwWPZHR8kSje4Bd7vdWLp0qaBWPT09Qg4j4azVanXVMY9EItPEY/x+vzDcSeAjXwOA/GwPQlh75jlR+IjXgFPQmE3zH68tZ4QXCgVpOeN1IBfATlizD4nhd4pkOcL0LDt4PB4UCgUpW7D7hMc7MTEh3QMMHuwwPeDA6Y7tnx0TDn3jxo3427/9W5x//vnYtWsXnnrqKfzjP/4jAOCrX/2qbMdFiyIRlUoFqVRqXrLzH/3oR7jrrrtw3nnnCVHmJS95CaLRKIrFomzHFiMNuWqJVDpjnW3odicuZlyQ6BTYR873MlN3ZCFnZ2yP0lr7wFTwRbJZIBCQcgYzNrKS2XbFrJz1dbvcqc7A6cQpecr/GcRRsY3EMmaIxWIRIyMjsKzJqWNEe8rlMnK5XMe4VNa87WQ13X1B2JmqcMYY0T/Qztfv9wvUD0wXj9GEPC17ymCBjpqBAJ//fD4v093o0GkulwvhcBjpdFr26fF4sGTJEkQikY6snpPJCKfze+7xeKTzg8EOIXhm9DpwJr+GGXg8Hkc2m5XZ7YThGXgDkPvIbhN9Do45diB21Dr0QqGAyy67DB/+8IexbNkyIXg99dRTHdtt2bJFfg6Hw4hEIiiXy4hGo3C73di1axeOO+44cZRzYay1EWI7+eST8cIXvhAAkM1mZYFkdsfFUi+oesElM1qLznDB0b25XIyAqfng3D+hVIpYHK1GB2z/eW+m+5w1GZBte7onHJhSdPP7/UgkElLf1cpsvBfkSFB8hMQx1nD1SFA9ypRBA/kTHGjCZyAWi8lccu6PTjabzUpNnHVxMtk10S0UCsHr9XZA75VKRa6B5l0wQOTftO5/sViUfnE6Pa06B6CDVGZXX6MoDtGAZrOJYrGIbDYrzHmy+1n714NStOIi73mtVpNRsLw3hUJBAl/Kv3IgkcfjQTqdxtDQkHSSMGDh900T1thZsGnTJrmezNKJTrC/3YHUHZsrO6oceqvVwiOPPILbbrsN7XYbv/vd73D22WfjgQcegGVZOPXUU7F69Wp4PB7ccMMN6Ovrw+mnny7v93q9oiTF2uCGDRtQrVZx6qmn4rjjjutoLTlQe+1rX4vnnnsOL33pS9Hf34++vr4Op6DHh1L1ivVW3XPOcyYEydYkOm3CuQAkIGH9nP2seqgLh3DM1tEdStvfY+omwsJryCAmEAjIsBjun06aUCydBfdh/xuDLgZcDL5IvGKfdKvVQjKZlJo4M28O2mC2zW4FlniYkZL1TBEXQrNasa1YLCIajQqxqlQqiQ47HZDP50M8HpfM1k42KxQKWLp0qQic6Po0kR8Nw7P2ywlnHo8HQ0ND0s/OujWdPwNGHiMzYDpdksp04ETNAzp1Ihf5fB7pdFpQLGbHnI2u75mukzO41yQzolladW3p0qUi/ELGOUsSLFWRnwB09oBXq1XkcrmunQPaHEjdsbm0o8qh/+Y3v8Gll17a8dpXv/pVaX351re+BWASxv6bv/kbAJO17LPPPlu2j0QiWLlyJTZv3oxms4m+vj5x6tlsFmeeeeYBO/VarYY///M/x86dO/HSl74Ug4ODOOuss7B06VLJZuiEuXjqmqnb7e7IAlgDZxbD9jIypUlKIhTPbJIZG4AO2J0L3kI0neXYza6iRgfNAIkZM6FrbgNMZdO8Nsz0NMTN7QB09PKTmEayGAeHsMWRnRE+nw/pdFqEe+jw8/m83B8Sprxer7xOZ8ShHgwWmMFynCbJU6zL0/lUKhUMDg4iHo8jkUhgdHRUNNrpTLsx1mu1mgQM9vo074HenkRNTfgiuY664qVSSYRSjDHIZDKCUJGASidLx14qlaTHG4DcC34+98Nng/V9iiVZloUlS5bI0CGS4Vi+oqgNBXT0dbCLJjHYovOlLn42m+3gSHDuAeWUOaaVAY9jjs2nHRUOvV6vY82aNfjP//zPjtd9Ph/e/e534z/+4z9Qq9Vw2WWXYcuWLfKlDoVCOOWUUzreEw6HMTQ0hEajgccffxyhUAhLly5FPp/Hfffdh7GxMbzyla+USH02tnv3bvziF7/AV77yFfT19eG0007DqlWrcP7554sOtZ6pzOyJWQsJN8wmmbHYnRLVvUqlUocOO7flgtVtgdRZ/+HOzjXEbX8NQIcAjnbyOkPTOtssSdi15rWimFZm03O9ee30UBRuT/U1ZtqVSgWVSkXgcAZZvN7UO2cNXYv5EDFpNBpS+uG8b8uyRAKUUqSEv4kmUQVQM8k5EEW3t3FSGACB9DmxjNeKrPZYLCZBZqlUkhKBHVZvt9siNsPAsFaryT1gYGqvx0ciEWGmUyPBrr7Ga0g4nWND+bren2b+k3AYiUSECJfL5aSkwfKGDpT5M507vx/s+NAjTynFyqlrOhB3uVzCAejp6XEgdccOmR3RDj2bzeJd73oX/uIv/gKvec1r5PW1a9fi29/+Nt7//vfj0ksvRSKRwEMPPYTe3t597tPr9WLFihVwuVxIJpNYs2YNvF4vFi1ahOHhYTz++OP41a9+hfe973046aSThGS2Nyd4/PHH44wzzsDSpUuxZMkSXHTRRUgkEh0OglE8dZu16Ihe8HXWqBcywuWaoMVsktmeXSKSAYHWgD6Uzlyfp/11OlFm3fo687xYHqBTYT2T50bHTiiUjphO1a59T2dGYhg/h/VuwuHcnmxnkqrowEgqAyC1Y60vTkEXwuNUgtPCQEQSmAXHYjH5ube3F81mE7lcTpjS5XIZS5YsmcYk7+/vl3GndkiZgzzYshYMBqUmXSwWO4aJWJaFvr4+USqLx+MYHx+XVi3C6s1mU9TNstmsSM0yW9W97gzAgCm0iQ6ZkH+3rJ7DXOz1fTLV+X1hdkyjE+f5sHtFIw1Lly6VXnaSGKnSpzXgOYrYTgDs6enp4Ag4kLpjh9KOaIe+adMmbNq0CT/4wQ9kUe/t7cU555yDc845BxdccAFKpRIeeeT/b+/N4+Sqqu3xdaq65rGreu6k00kgkIEhQh6gMgR5DKKEScQBcVZ8gD6/Av7EJyqgPIGHqCiT8FAmBZMHMvMRQhiDhEnmMUCmHmuuruqqrvP7o2rtnKp0BiAM3Tnr88kn6ep7b917q3L32WuvvfbjWxTMCb/fj+222w6Dg4OYM2cOVqxYgXg8jpkzZ0oN8e9//ztWrlwplGkkEkFfXx9isdgG/4Gvv/56XHTRRXA6nTjqqKPqZkvzQcCHPYMbAwkzAT5wGIhIrzMYNw5V4R8+IJld8WFGdyqTln6/zWI2thDiObH2bPZbm3VrYH2PPjNX3i9eE1mORvqd2/E+ezwecVUzKXzeb2bwLEu43W7EYjH53EjjMpPjLGz2ZTOo06yInz8FZVww8NwpNjPtU7mAMb0HGqnyYDAoBi0tLS11pZjGDHl0dFQCOt3UeJ/cbnfdmNBSqSSlAgAy5ITiPQZ6k44nc8CsO5FIoKWlRdiHkZERobfJTPB7wetrtBsOh8N1i47Ozk7RFXC6GulvXhvpdtN5jddOVob/H6iCp8qdzBj1KabozryXkUikziTHwuKDwIQO6LvtthvK5TKeeuopecAvXrwYAHD33XfjkUceEX/0tws+XCORCFavXo3Zs2fLyvuTn/wknnjiCdx3332IxWL417/+hSOPPBKnnHIK9t9/fxxzzDHy4Fu9ejUuv/xyfOQjHxEWgVkEgysDObNRsx/Z4/FIxsVgQjqQFDuDH4Mjsx3W4blPI6XObHK8LPn9RuOD0KxLU7xGoRkXNHxgM1hyPwZ3ZtDMEDkRiw9qAKJ0ZoBhFs5julwu6Xjgv2m7W6lURLDIgFQul4VpASD1cXMWOZ3MGvURJm0+OjqKaDSK5uZmEbg1jgGlzzqvtVKp1NmrDg8P1wVc3gPqBMyAxQUB6XKKBs2MmtP4OFyFQjK6EpJdMjNe2rZyQRkOh8WaNplMIpVKSY2bve+mWxo1IlwAmQG4UqnIFDcAovbn/49UKiXHcbvdojlo/M6Txmf3Ab8XxWJRXOrMFlFqAzwejwj/LKVu8WHAVgvoSikvgGUAPLXj3qi1PkMptQuAiwEEAawE8AWtdbq2z7kAFgL4f1rr+5RSvQBeB3Cy1vq3tW1+B+AxrfX/jve+jz32GD72sY/h0Ucfxec+9znsvffeAICrrroKoVAI0WgUH/3oR9HU1IQ99tgD55577hZfEz2p2VLE8aHlchmzZs0SynL33XdHuVxGNBrFsmXLsMsuu2CHHXbAE088gfPPPx9Tp07FoYceKg90tjEppRCPx0XYZParGvdVjECYRQCoC+7MoEghc1FAsQ8XDY0UtxnkP0jwIWoKy0wbUdLswHrrVA4vYUmCr9OmlAGKgY6Bl21XbBdqFEnxQQ5A6q65XE6GoXChxGBDFoDjPLPZLGKxGDKZjGSObCVsDCKk98nUMDBlMhlpR+N7NWaqXDDGYjFxKwOwAU3Ozgxm+e3t7fJ9LhQKdZQ2x6MGg0HphWcHBYe/cLgJF04M5LwOl8uFGTNmSP2fvfAM9qStnU6nCPX6+vrGVcwzC582bZq0tfH7yoUd7xkZqcZFjzkcxRxwwhbAQCAgDAatX0mpl0qlOlOd7u5uEd0RllK3+DBha/JDRQD7a613AbArgIOVUnsCuBzAD7XWOwFYAuAUAFBK7Vjbbx8A/2Ecpx/Ad5VSW2ScvmLFCixatAj77LMPrrvuOgwNDQEAzj77bAwODuJLX/oSHn74Ydx///3Yfvvt39YFuVwudHV1YfXq1WL2QtrR7/ejpaUFLpcLnZ2d8Hg8WLBggdRBX331VfzpT39Cc3MzDj74YBEX0UmK9GwqlQIAoZH5ADc9voH1AjEGOz5E+HDhQ4f0MilVZkemcOfDIn4jGql31scZvEwfa7PsQP0BVdJut1tU0VwENDc3S821kZ4dHR0VBzPuRw973n8aubAWHAqFZH426eyxsTFEo1FRerPXmJqFYDCI5uZmEVKajmZkUdh3bS4SeF9oBUpGhZk0RWFOp1O6MMzAr2utcOY9ZbA03eLM7b1er2TzLEXQnjUQCEhALBaLGB4eRiqVwuDgYJ3mAQDC4TCi0WjddvzMCoWCLFr4WTHD5b1iXZ6Ctfb2dng8HqTTaanVk7lh+5rZJqa1FiaA3xlOcON7sFZO6h+oivOGhobqzo37mDV+C4sPI7Zahq6r/5OztR9dtT8awA6oZu4AcDeAOwH8FwAngEptG3Ns1wCABwEcD+CyTb3n2NgY9txzz7rXXnvtNWlPASB/n3POOfjmN7/5tq/L4/Ggra0NTz/9NGKxWF2fMlXKfGiuXr0aXV1dyOfzeOmll/Diiy/i0EMPlQwJgGSOZrZMmpcPQwZdBqzR0dG6dik+yIH1s7KZKfAhxIycPwPrg+GHrc+cCxFgPf0MoI4VYSAklcoHNa+FD3/2+1IwRV2BuQBigOR9Nh/cfD0QCMjsbGZz1BxQsMUgEAgEJNiMjY0hnU7LdK9MJiNe+q2traJU5+ff0tKCdDottXNeLxdjVEvn83l0dXWJUt0cGZpIJCSzpmKd51cqlWQ6HssU2WxWKGuTVnc4HOI8R1tTfmcAiPCMAdWsR5N2HxsbQyKRkHtdLlenmTHbZmZtMk+k1FnmAiCLFjq5KaXEfY1lBdMy1+v1ypAWfsdNS1xTpMjvg8/nkzYzGuwA1TIJr91ky96r4UwWFlsLW7WGrpRyAlgBYDsAF2mtlyulngFwGICbAHwGwFQA0Fo/q5TyA3gAtazdwDkAbldKXbGp93M4HDj44INxxx13yGvXXXcdbrjhBqxatapu21NPPfUdXZPWGl/4whdw5ZVX4sknn8Ree+0lphK0iqQwbdWqVfjYxz6GNWvW4Mknn8QhhxyCXXfdFQDq6q6cAkWRFltbGJzN4SiNw1TMCWk0tmAmRaU0szGgXlxHmv7DFMwBSJnA7B/nNZCaNdXjzJQY/IrFoqitWdZggOY9U6o6tITHYODnQoIZPkV3Ho9HFOi8p6Sdea4jIyOSuTLAMejTg52BlQKxWCwmrWJkBMz2tKlTp6KlpUWobFNfwcyZgi0GVADSh57P5xGLxUSkRjX80NBQHa0dCASkFs7vAzsvKC6jcp/fVVNAyECZy+WkI4QlCbbSNVrHknXhWFKgunhtb2+XAEyKP5/P19mh8ngmDe5yuTBlyhRZMCeTSRHdkTKnxqJQKAiTA0B0D+bEte7ubtFLEHbWuMVEgnovHu5KqSiq9PpJAMoAfgMgDuBmVOvj8Y3s1wvgFq31PKXUn1DN6PfARmroO+20k/7Od74jq+9p06bB5XLh1VdfFcoOAGbOnLmBQ9OWgvU0UprNzc2SUZqBiGIis/2LPbDM6BksGg1c+Dvzb/P9+TcDc+Nr/Hfj66b4zTzWe4V3kvmb19p4/ea9M9vTzEEnVGXzswDqr9Pcr1GnYDq/mcwAYXYc8FhmnZ/7su2QNDjNZBrbAXl883PkwoXnxsUKz8l8b9PlTmtddw40ZOFCzizFUCTH96BYkGjcx9zPfA8uohq/c2amTDMYnrt5nmbXAffnvmY5iP82TYBMAaFpMctFHY/FEg2ADXQj4+kyxjvWuwXvwbaKbf36gff+HixcuHCF1nr3xtffE5W71jqplFoK4GCt9XkADgQApdQsAIdual8DvwBwI9bT9RvA4/HgiSeewGWXVZn5G264Abfeeqv83NHRgXXr1gGAKJzfCfL5PN5880389re/hcvlwgEHHCAUKVftg4ODeOCBB1AsFvHxj38c3d3dANYHF1PUZj58TVOSQCCARCIhv2PGSSqZ2TkXBWYwNLOoRsW6qQJ/L4M6Kf8tAQOjKYhjsDX7sJnJNTU1oaWlBalUSmhTPqRpicp7yaDKrI5BqK2tDUNDQ0Jt+/1+JBIJodibm5vrerUrlYq0j/Fntj+mUimhm9mypJTCP//5T8yaNQsAJBtlnX10dHSDPmi32401a9aI0ZHpT9DorU4WgudEQWShUBBrVA5aYZbLDH9wcFCCKmvUHE+aTCZF38GZAMViUQRwvM5oNCrfXZZ4+P78//XYY49hp512wtjYmGgJKpUKkskkSqUSWlpahFWgy1osFpPv/MDAgJwrr5k1fLfbjcHBwbqZ5DNnzhR3vdHRUcTjcQnKhUKhzkjGtPBlyaa/v1++k52dnVslC1+6dCn222+/d32ciYpt/fqBD+4ebE2VeyuAUi2Y+wAcAOC/lVJtWut+pZQDwI9RVbxvFlrrF5RSzwH4FIBHN7YdvdiPOuooLFq0SOrkV155JZqamnDccccBqGbpfX197+ja/H4/Ojo6cPjhh+PPf/4zli1bhgMOOADxeFx6XH0+Hw477DAxHWFwAtZPTDP9vmkIo5SSudLM0JhxcFuzBt7omMZMhjVAM8gDG2ZdHxQaGQNqEMxBMOZ10sOchjsMNAwqVFbTQMR0+iJ9S+qbHuAjIyMSkBhoCdaFTac1HsOkzXnOZp07FAoJS0Na28z8uShprDtT8NbR0SHztSlUY+3d3J4Bnxk3a+blchmpVAoDAwNCJTPTJ/XOenY8HpcWLNbazdY7njfvKQApHbS1tQGA+LHncjnxph8eHsbQ0JAcq7m5WQK30+lES0sLhoeHhfIfG6vOHGDfN/v1eV9MfQJ97XO5HJxOJ3p6emTRYbJjfr9fgjlFcOaCqKurS3QShKXULSYTtmaG3gngqlod3QHgr1rrW5RS31VKUcW+GMCVb+OYZwN4YlMb7LXXXgDW0+1HHnkkXnnlFXzlK1+RbQ488EB89KMffRtvuyECgQB6e3tx1FFH4cYbb8Qdd9yBnXbaCR6PRyYwsec2EomI2CYQCEhGAVQzTXqtU3jFeiCz6EaTDbY+UZhF2rRYLIpojpQo6VMGdlMp/35gY/S+Ka5ixmzS6czMuR+zTwCiUzBb2kxRGkVxlUoFwWAQ06dPR19fnyyUxsbGEAqFkEgk6gIWz8tkQihSMxdKXAxQTe73+zFz5kyk02n5DLig4gLM7MWORqNoampCMpmUPnCK8YaHh4XVYIA0hXxcnPDzBSBtdVT+8/6Zwau7u1u0BwCkM4L3jGNdKTLjQqSlpUXMkZRS6OzslMUF39Oc69248GCJijoFZvJsFUwkEvKdZ6ZvHsPsO+eCxfyM2Vro8/nkWDxXltm4MDDPi34AjcI223ZmMZmwNVXuTwOYP87rFwK4cAuPsRLAPOPnp7CZ1rrrr78eSimceOKJAIDLL79clLKtra0YGBjAXXfdhUcffRRnnHHGFl7NhnC5XEKjn3DCCbj00kuxePFipFIp/PSnP5Vs2XzIKaWkn7itrU2cv5h9U11s1oHb2trErYyZPd3k6DhGipVZPrMStqgZ9+8dX+87RWNt26TWWbM2a+NOpxOxWEyyRbpuFYtFhMNhWbxQlcwHPWlstpBRyc0eabP9y8w+zcw7Go0K7csMfGhoSFoTzSyeg0w4hxuoLkz6+/uFso7FYsIQ0POcizEGnba2Nvj9fhHbmTO70+k02traRDVPxXkwGBQK3OVyySKChjPs6uD1lUolabOkot+ESas3sgajo6NIJpPSWcCyQaVSkUUqsF6kprWWQMqFRyAQENU7wQWFUkpsaJuamkQRb6rUOa3NbBOjKp8BXqmqT/vg4KAY5HDBxk4GDpghrErdYrJjQjvFAdU2tWg0iunTpwOADF8BqiYbADB37lw88sgj7/q9/H4/ent70dXVhQsvvBCPPfYYrrvuOqmx0omK9e3m5mYJ4rSyJJXMoMWsjs5hzOoYpAOBgFCLTqdTHmgMWlQFs3WIme/7GcxNOp30sjm1DFgvgCLMYEvjFrIZVDNns1kZuEEdAYM+gxwV5qbJiJm50uIzEonI8WjhmkwmJXhTlc76LT83s4bPoJdKpSTDNRXcVEiz7suFB7NueoKbveSNAZXBiCYudKarVCpiAsPrpSGO3+9HOByWIB4KhRCJRFAqlZBMJuFyuaSkMB6trmstZW63W2h1BneyAY31fLbJsV7f3NyMVCqFdevWCa1NOpzv26hSJ5NAy1meP30BaOpDjYnJMPD70UipNw5NspS6xbaECR/QjzzySNx6663y8yuvvAIA+NznPoehoSHcddddePbZZxGJRCSDeOKJJ7Dbbru9o/fjQ2HZsmW46KKLMGXKFBkaQVW91+sVlyzTAMPMMk2KmEHZVL+TiuTigMGO9D2PxWzd7E0H3pni/O1ivLo4sL5uz3vVqMZn1m4GZKWUiNmi0agsdJhl0e2LWTPr1MViER0dHZKtcsHT1NSEGTNmyP3K5XKIRqN107zYH83gbdbPKepi6xcXB16vF62trWhqasLAwIC8Tn0DaWO2Znk8HmSz2boZ4szUaViSyWSk7EK7US7ifD6fBMLGYGjOLWfd2uPxyOxy0tScMkYWoL29XdgEAOjq6pLJYzTzYUmCLWzj1f9NJb7L5RJxWy6Xk+NQ1W46G/I4SilR3NOkh/V+CtfoW8/PL5/P101f47FYCmgM2pZSt9iWMOED+n777YdisYhVq1ZhypQp2HfffbFkyRLk83ncfvvtsl2lUsHHP/5xtLe346abbkIikXjHrWz5fB533HEHMplMnbEMsL7Wx8yJWQXNSUzjE2Y/fIBRzGQ+QGl+wfcgDUn1Nunr9yOYMxCTPTDdwUzFOkVsDNrmFDRm6dyPdDQXNrwXDPQsMfB4uVxObEEZXCiWY5ApFouIRCLSL/3mm2/WLSTMBZCZGcfjcZTLZVkYmOYyM2fOxPDwsDjTpdNpNDU1Yfr06SIyo7GKmd2Ol4lTb+H3+1EoFNDc3Czug9lsFqtWrZLz5QIBwAbnnM/nRbjG7w4XEszwxzOBGRkZwdDQkHw33W43mpubMTg4iGQyKZT5yMiILGa48DCzXZNSZ+BUSsl3lUI8HsPr9UKp9X7r/D4DkPKCOaI0EAggl8vJiFIK6sgsUCxHWErdYlvHhA/o3d3dcLlcWLFiBaZMmQKXy4XDDz8cgUCg7oEDQGh3ery/UyxYsAB77rknpk+fDq/Xi7Vr16Knp0cyBI7t5AON9fR4PC4mMWyxyWar5nperxeZTEayNWbkfDiShjXFWsx82GPc6Bf+btC4KGCgNXurzX5pM7CbPuU8BoVRpH9JYZNK5yJo6tSpUjLh4sXlcknWNzo6Kl7cIyMjIpwaGRmR0gfNTsiSNIrUOLGL4jczAHZ0dGBgYEDeQ2uNaDQKl8uFeDyO4eFhrF27FkCVGqaVLOv45XJ1lngymawbScrFCdXdXNjQ65wUOevSPF+XyyWuc4ODg9L2yOw6kUggFArVtTGaIjOPxyPKdi5U4vG4lHcoxGv0QWfgNL3Q/X4/+vv7EQwGEYvFpJTEP2Q6GpkIHoPXy/Yw2uy63e66yWsA5P+NqXrnwoIaBsBS6hYWJiZ0QOeQE7fbjVdffbXud2ZQ++1vf4uddtoJf/zjH5HL5XDxxVvUObdRHH/88XjhhRcwb948vPTSS3XtTgBkZCMDHjMeBmNmV2Z7Gx9KzIRI+QYCAWQyGXFII03J1iG+r2k4sjXQmPEzoNOAg0GMdVwzO6cmgD3ipGSDwaAwE1yosB5ObUGhUJB7w15tLs6odDf7q0lHN2ahfPibdqDmokNrjVgsJpkxaXmTJWHGyyBMX3KzhJDL5eQ+NdLJrH9zAdPT0wMAMoJzbGwMmUxGesHNITEM5lxMZDIZKFUdtgJA3iOdTstCgrVqLir9fr8EyO7ubvT394uXO8FzYAmD780FCLUPzLDpz85FHM+DbIup8q9UKuJIx3tFrYmZzXNRSIqdjAwA6QjhooLHIiylbmGxHhM6oL/++uvI5/OYP38+fvazn+HLX/4yYrEYfvSjH0mW95e//AVHH300HA4H9t13363yvqeeeirWrl2LbDYLn8+H6667Dh/96EeFHmWtlQ8oBmNm3oVCAbFYTAazUPRjWp9SIW2yDFS6U+3MLGw8M5mtBbNFjAGLYjxmtuZrtDQFICUDZlSsh1JAxnYx03mtsQ98bGwM4XBYsl2XyyVjQc36duOoT2bunAvOKVmsKTMA8hoogBsZGdkgaG+M8qaQj4s0s9ZOK1kAYmNKi9JAIIDBwUGsXr1a3qelpQXhcBherxfxeFyMUiiiGxkZqQvGzI5ZPjDPubW1Vc6D10VtQjqdFnEgF1CVSkXYJfP+sb8bgCzk2EmxatWqOkc4fgamD4Dpm0AK38ykWWIYHh5GU1OTLCp4r3g8TsazWbiFxaYxoQN6uVzGgQceiD//+c9YsGAB9thjD7z00ktSh/vKV76CY445Zqu/LwVAd9xxB5YsWSIPWU7molkH55ibYji3241MJlOXVSaTSQCQTJEudADESSuRSNR5tTMwbq1gPl7dnbVyBnNzAhxV+T6fr+4caLRi1v4ZmP1+P7xer2SiFAz6fD5p7yuVSnIs1t4HBweFgmXJgkEjHA5L339HR4cIEtlmlslkEIlEpK7fGABZlyW9b2bZAIQyV0phaGhIFh90pTNFjnRemzJlCkZGRiT7JAtBRoGLssZ2MXM4SCwWw9DQkBi3mCIxshIM4I3sBAN1Y8mEorqhoSGh1ZlB8ztNcRzFg2QnyAiZC550Oo3Ozk6Ew+E6S9e2tra6ITDsCDDbziqVikw2M1XqjcYv/P9mA7mFxeYxoQM6ADz44IOYMWMGDjnkENxzzz047rjjcNVVV+Gyyy7DjTfeiN122w3f+c535GEDVG07I5EIjjjiCBx66KH4+te//rbes1gsStBasGCBZHAMMqR1gfVCHWZL7PGlOI7BDEBdcGDwZ20XgLQBkeJmRtU4FeqdgveI585pclRRkyEg9cnWKV4PAzyDL4M7FwSEaYBjiqMYFDhKlBnc2rVr5aHf2dkp2ezg4KBQwADEwMXr9aK/vx/r1q2TliyOJuV7MgByJCoHerAbgcGfYzO54Ojt7UWlUpH+b9qOcigJ2xWpiWB5AIDoJbgtzVcCgYCMW+U2pLsZ0JuamjA8PCy/83q94npXKpVk8aNU1QWPmgsyEfxOmYHf9LOnyp6Lx0KhAL/fLwp0LrTMe0d7YiriWf8n0+B0OpHNZoWlIs3OGr+5qNmY8YuFhcWWY0IH9J122glutxsrVqzA7bffDqUUlixZgquvvhp33nknZs+ejRNPPBG33XYbbr75ZpTLZRxyyCG499578elPfxp///vfNxi/uiXweDy44IILcNZZZ4mjmVLrp1CxzYYPPQrCAMiENACyHWnO4eHhOgU8H6J8qLJ2zge3SQu/G5jtZqQ8zYWD6T3PiVrM2tjHTDtWKsl5PFKuFFVRbGaK3ADU9dKHQiHJABszWdasfT6fZIHMgMvlMvL5/AY1dXPKF1XfPB4XAJxs5/V6xTOepQOKt7iI47Ulk0lRu/t8Pimb0KzG7N1Op9MyzYsLCwDSRkbRGA1beF0mlT5lyhRZYBDsM582bZp8nzgFjp8nmR2T2WBLXaNSnd/NZDIpdX6+Pz97ahLMoSn8bvP/Bz+jYrEoc9rN3nlzYhthg7mFxbvDhA7obrcbjz32GL761a/iyiuvFOVuoVDAjjvuiG9/+9u4+eabsXTpUrS3t+NjH/sY7r33XgDA3//+dwDA9773vbf9vqVSCTNnzhSF8+mnn47BwUEJFpx7TaEQs1hmQqaVJx9wzc3NdQr4aDSKoaEhABCq2VSTA+uFawyib1cU19hHzgBAyp/2snRnY40ZqGa5VIvTvGVgYEBasXgdfGgzKysWi6L2JwsRjUZFKMaxnZz/3dLSIrX3cDiMaDQqmWxjwKRlKUeMMnixxYnBLRQKoaWlBQMDA6K051CcSCQiLWx9fX1Yu3atBH1m+6wtj1drJ1PC7NhcjGSzWTkWhWter1cCNMWU1H+M5+JmTgPja8FgUPqwvV4vEolE3dwC1vr9fj9isRj6+/vh8XjQ2tpap1Tngsm8rw6HA1OmTJF+d44o5ThgZvKk5s32RbIDPJbT6ZRFDa/DqtQtLLYeJnRAZx3xiiuuwA477IAzzjijbkziRRddhAsvvBCrVq3C5z73Odx88811+zudTixfvhz77LPPFr9nOp3GySefjNWrV2Pq1KlYvXq11B5JmefzeZm17Xa7JRslNUqwvsh9mKUy06Mgjepgs1ZtZk5mq9p441mB9cHbpOdNy1gAImgzfddDoZAIwCiWolEO2+hYD6VNrcPhkGszh5mYGXl7ezuGh4eFrqWxjrldpVIRkdjo6KiIvbjAGa+3m259LpcLM2bMQCKREDFXIpEQwRWnjiWTSaxevVqsYM2gatZ3Y7GYsBambzs/V9L2rMWTvWDHQjAYFEEkRY2sdbNbw8xefT4fcrkc0um0ZK5siWMZg+151DEwKI9nQMMecNbnzZYxZt58D/O+8nf8TpCpoKVsIpGQz54LS+oHuPBh14N5LMLWxy0sth4mdEB/+umn0d3djdNOOw1z5szBqaeeirPPPhsvvvgi5s2bJ8FhxowZuOeee7DffvvhC1/4AsrlMq6++mo89dRT2HfffaVeuiX49re/jTfeeAPJZBLRaBSXXXYZ4vE4Xn31VQmW9CVnHy4ftEopmRJmBmZg/YjTpqYm8TYHIK5m/DeDsjkSlRQqM3hS3Y3znpklNbZwAetd3fjQZbbudrsl6DKYmlPOisWiWJzyurhIYWYdi8Wkzs72IxqxMCh1dnYiEonIvePiZnh4WBTkyWRSGAW32y0iMXOQB+vgzEjZv7169WqpXVOgtbEsG1i/WOTigi1b1BMUCgW4XC4J3qYArlAoIBAIiAFOoVBAW1ubDOApFApIpVKyYEin02htbRXlPz8TKtYpFCNDYKr0XS6XUNrcB0DdYoN2uhw0w8UZFws8Fj9/LqxMkR6AOidEfmcofiSlXi6XpeXPZEYIS6tbWLx3mNABferUqfD7/fj5z38ugqE99tgDn/jEJ+QBTvzqV7/CY489ht133x29vb146qmnsOuuu+Lggw/Gn/70J3z1q1/d7PutXbsWDz30EL71rW/hlltuQaVSwezZs2W8JOuWZn2YrmUUR7GGTjqaTmIMRsxMKTgqlUpiW0sRFrBeIMes0czWWQog7QxA6sSkRslihMNhGYTCc+QoSrO3nkGWdDZV17yOxrY70/ucgzISiQRaWlrgcDjEdtXUHzATZmDTWmNwcHADWp0Wpw6HA9OnT8fAwIBkfhSokW3I5XJb1IZWKpXESIV96KThTdMenofL5UJvby/GxsaQy+XkeqjcZ6mB953vyc+tccHAIGpmr43HAqotjuZAF7a80cOcTEcqlZKWP85wZzmK2bNJq3NBxc+/p6enrtTDhQlLFLyWRiMZpapz1c3AbWl1C4v3BxM6oL/11lv47Gc/i7/85S/yms/nQ39/P5577jmZlQ4Ahx12GNLpNBYtWoQ777wTAPDSSy/hySefBFA13jjooIM2+l4PPfQQzjzzTMyaNQsvvfQS3njjDZx99tnSO0uKPBQKiaqXPcrsJ6d5iWnQYip+gfWta2znIt1OMHNm9kcamZPYSCnzNfYN87xI65KmNSlzp9MpQzYYmLlfJpOR8+7q6hIxn0nxs4/c6/WKgpslA2b+g4ODdUI0Cq8Y0GhAw6EmJv1Lmp2OenSRCwaDyOVysojzer0YGhqSbNjcn4pyh8OBgYEBocQBSG2bpYP29nZR+o+MjCCRSMiCiFazJq3MLDqbzUoXArN2Bj0ueniNxWJR2sYoAmTPPO8lr4OteI29+vxjeui3t7cjn8+jvb1dgigd2RjQTVqd7nNcnNJnPpFIIJPJSNsduxhoANOoVCetb8LS6hYW7w8mdECPxWK48cYbMXfuXDz99NM499xz8ZOf/AQAMGXKlLptd9ttNxnIsscee+Dkk0/Geeedh//5n/8BADz22GMbDehaaxx55JGYP38+Dj/8cFxzzTUol8tYuHAhUqmUPMjp6c2HJR/Opi87AMlYR0ZG4HA4EIvFJEAFg0FRJZPGZobPuqX5IDUzebPflwIv0uIM9m63G0NDQ5Jp0UqUgY9jK6nM5+ukfdlP3tzcjOHhYfHbJtXvcrnEY9vtdostKSeFme1nzc3NUvs1Pc/7+/slIMbjcRnd6Xa70d7eLq1aLGWQmjfr3RxBynNhkGUrFtXyPT09UErVqcPJNjQ3N0sgKpVK4iluZvqs3dMIhmWejWXCnC7GbJwqd+6Ty+WEQqcbHs/V5XLJWFd6HLAPv1wuS7mCCxPWs+my10irUzRIdTsV9tR/ABBmhVPeTDaIbnQmtkbXhYWFxTvDhA7o06dPx49+9CMcddRR2G233XDbbbdh/vz5OOigg5DNZhGLxeq2/8UvfoHTTz8dzzzzDObOnSu16V/+8pf44Q9/uNH3Of/885HP53H66afjtNNOQ6VSweLFizeoNbJ1iQ9VZrButxuDg4OSlUejUemxZtbm8XikFstaL49Ly1hmf3wwsy2IbWOBQECOzwDChzQDOzNvs5WO6nq/349MJiPnzvcxHb/orGbWwGnxCUB6tMvlsvRCNw4IYW89M7nGLFap6ljNaDQqGR8NT0x/dxq0NKrBed7crrOzs87fnKUPj8dTN80tkUhItwKV/qSd+VnxPchmsKxBCn1wcBAtLS1S8zfPzbyfvJcUi/HzdLvd6O/v38APneY3RFtbW10Gzo4KBnDeT5rIcLobvwfFYhG9vb1IpVJSSjFbAonGzyaTydSp3gHIPXG73Xjttdc2/x/XwsLiPcGEDugAcOihh+LCCy/Ed7/7XXR3d+PQQw/FCy+8IL7ZJk4//XQAwKWXXopjjz0Wl156KYLBIG644QYsWbIEs2bNwtFHH41FixbJPrfccgvOPfdcfPazn8VXv/pVDA8P46STTkI4HJYsiQ/qdevWobm5WWqJdCUzRXLshWa2x1p0JpORB7Pf7xcqH0Cd4IgK6IGBgTrHMk7CMm1MqdoGqg9pM5BTUMZFCQDpSyYNbjrcdXR0YGRkREZkUtBlDj3h4JOBgYG6mrfP50MsFpN2M6/Xi/b2drFM5fQt83h00yNbEAwGkUqlxI+c2avT6RS/dR6bynxzZC3bAMkEAJBJZWz3q1QqMg89nU6LLzoHpNBalgslit9YZ6eSnGYqADb4LNmrz/59c1FFdX1jfZ0iOp/PJ6I50+OAMNkD836abWtaa7EipiqfVD8HxZCZYYmGx2IpxyrVLSw+nJjwAd3j8eDkk0/Gt771LZx33nm45JJLMGfOHHzmM5/Br3/9a3R0dGywz0EHHYSZM2cCgGTy9957Lx599FFcffXVEkhKpRIuvvhizJ07F9tvvz2uueYatLe34ytf+YpQ7aTNqRonZcn+YmaFDIwUEpkPbAZms63LrPua6nStq+5hXDCYTnK0naX6PZvNiuiNbAEp1PFsac12N4/HIz3zrMPTrYyZJQMcA/brr79ed+48JoV0nZ2dIl40W8+UUnKe7Gsmu8JA19iKFY/HEYvFpG7r8/mwbt066VU3687pdFoYDiqy/X6/ZLJc9DQ6qdGxjkGSn1UsFhMVP8+X9D+zf4rZqPJnrziz6OHhYbl+LoBMwxXqLYLBIJqbm1EuV8e6mi2QIyMjwhxwYApFkeNR6gCk3ZD++l6vt87f3tQfUI9hTjcDrFLdwuLDigkf0AmPx4PTTz8dp59+Oq6//np8/etfx/3334/Vq1fLNrfddhuWL1+OQw45BEopLFu2DDNnzkQikcC8efMAAHfccYdQm8lkEolEAqlUCjfddBNisRhuvfVWAJCRnWZ/OGu9bMtqb2+Hz+erGxrCzJjCKwZfBnkGj8YRm6TL4/E4AMiCg4rudDotFqimxzizZgYoM7NyOp3weDxyPlxYMNNsampCc3OzCNvMgNrb24sZM2aI/SprqaahC93hTNEfgDoquKurS9rBOJ6T/ukAJEA2OsaRBqeZCksGjc50Ho9HDFAajWji8bj0ypN2JpvCLNX0FOC9oYiMn62pD6BrHe8H28UoxOP5mrQ6XeQikYh8Ph0dHRgcHJTFAO8FBYqN18FFJKn9QqGAqVOnitiR94WDYkjzkzlo9GePRCJ1/uxWqW5h8eHHhA7oG/Mv/7//+z+4XC6sWbMGiUQCzc3NAIBDDjkEhxxyiGy39957AwC6urrw9a9/HZdffrnMui6VSvja176GpqYm7LDDDhgcHMRZZ50lFqHM3DjWEsAGLmkcqsEsZ2RkBH19ffIQpYqdzmAMSgwqWmupVTJDZ+8zX+PMaABimWqq1tkXz7YjbkM1N1kELkIaLViplDeV1cyW2eedTCZFjR+NRsWFjaYwrO+OR61TUEdqPRwOI5vNIpFISIAkRW3S6szCWcun2QuzbtZ+ze9I43szAwcg+oSOjg7xqadIjd79ZEuoW+AiwTwuAMnQTbaE5j9cnIxXYuAixHSRMxdgXDQ27ktTG4I98bSqTSaT8tmWy2X5DrC8ws+OSnWzrk9YWt3C4sOPCR3Qn3zySfzrX//CTjvtVPf6zJkzkUwmcdxxxyEajW7yGHfccYcE+dbWVixatAilUgnnnXceEokETjnlFMnezf5ys5ZoUpAMlmxvYq05l8tJDbOx7mz6ipPmZlsXgLqecrYVaV2d500hE4MmAz5h0vumvSczbQZlc9gGgDoL1qlTp8roTa2rwzeam5vrMkVTN+ByuaSOzQBI8ZzJWJCWNsfENorjurq6xI41FAqhr68Pfr9faGAK9Bqz1tbWVgnM/JNMJqWbgApwLhbIUHCwCLB+XG06ncbAwIAcv7+/H83NzXWsAGl9oErNDw0NiQCtUChIy5spYOTigy1irPvTG540O6l5MiihUKiO9eH3hJ83RZD8brE8QnGlKfJzuVyyqCEspW5hMTExoQO61hp77LHHBiK4t956C16vF3/60582e4wrrrgCAHDZZZfh8MMPR3NzMx566CE88MAD2H333TF9+nS43W7Jjig6M5XUrF+ytg1UJ381Krs5mIQPYk7kKhQKdTQ/szj2irOGTRV9Y7DneTAjZEBg5s6AyWBPcdPY2Jhsl8vlxN2MWbmZPfLBz2yVDIR5feY0Od6jt956S2rukUhEgtng4CAikYi0lY2NjSGVSm0wczyXy8kgEQYnUzjG+9nYm00BoZlpVioVGbHKbJWlBGbaXPiYnzWZAvNa2T7GgEpHt7feegutra3jUuNdXV0Ih8PyWTcq1c26P1vFuLBpb28Xw51yuYxYLCbXxkE+LpdL7h8XJQCkS8BceLEsw3NhZ4Ol1C0sJi4mdEB3Op0YGRnB7NmzsWDBAlx88cVob2/H4sWL5aG/qYfTm2++iRtuuAELFy6UEaovvvgifvGLX6Crqwuf/vSnJcsm9cmeb7qIUcxF5bjpdZ1MJiUAUEnNiWPmoBNzlCj70kmPm21knC/Nui6DPWn+QCBQJ8yixzmw3vSF1D+FeubipFAoyKAY0sIUcyWTSQkKuVwOra2tMlAmn8+LYtocC9qYabPn3OfzIRQKoVAoyGQ3Bmj24Y+OVsfLcoFB5T0ZD5ZSeN1UoSulpH0vlUrVTfViGcGc4e5yuZDP5zeoSVPtbdrgUhtgDrXh58f3IbU+3nAWcyQp++w9Ho98R1mnp70s93U6nTJqF1hfl2dmPTw8LGUg+sRzwVIoFMQylwvTRooesJS6hcVkwGYDulLKC2AZAE9t+xu11mcopf4CYIfaZlEASa31ruPsvxJABsAYgLLWevfa610Arq797gta66xS6qcATgXQq7Xur22X1VoHG48LANtttx1ee+015PN53HfffZg9e7aIzo444ohNPqBee+017LBD9fRvueUWAFVr18985jOYNWsWZs+eLWIv9mlTDd04QMSkuJVS4sZFZTdru7QHHRkZEeUyhVZ0AAMgQYTztNkzzXnVZp2eNW9S7n6/Xyh5Zp1mXzxd4zg8hIsR/p4BkRPjaOTCmeRmySEQCMiQlZaWFrzxxhtCDzdmtOZccPZtDw8Pi4kMhVik/4eHhxGPx8URr1gs1tnFMpuMx+OSJdMSllkw7x2vm4p7Lvb4WTa60QWDQelPJ0qlEtra2sSil+fv9XqF3aCTWyqVEkMgOtHRFpjUvM/nE6U6FysUJ/L3pNXZbkjwc+diLRwOI5lMYmBgQHzsk8lknT+CWeJprI9bWFhMDmxJhl4EsH8t4LoAPKCUul1r/VluoJQ6H0BqE8dYqLUebHjtZAAnAZgB4IsALq69Pgjg/wE4bXMnFgwGsddee2HZsmXyGrPeT3ziExvdr1AoYObMmRJY+PBftGgRDj74YCxfvhwLFy4U4RndyPgAZXtQpVKRhzewflgFa82N3uFm0PD5fDIgg73YzPDcbrcYfTCAm33QPCc6mgFV1zwzSJstcQxkvF6e38jIiNC+DJBkH0g7m4M3zDYoZuBmfZoByufzibKebnlcGLAuPJ4ZDGu87LdmJkxWwjwHXgvrxfRLp6c5AMmGh4eH6zLw9vZ2cd1jTZkqd9LwJh1tMgvMbDnnmwudVColrXbhcFi86dnWSBMamuGYjAf79ZmFu91uZDIZUeGPjIzIZ8rFkplh8/tCb4N0Oi3+7lxI8POwlLqFxeTFZgO6rnKL2dqPrtof4RtVNWodA2D/t/neTgCV2h/TL/IKAF9WSv231np4cwe599578dxzz6FSqeCSSy7B73//ewDAiSeeiBNOOKFOJWzuAwBnnnmmPLRPO+00zJgxA3fffbfMySZlSXqazmgAJHhQKEfTGAZ5BirSy3SKM7NUzt0GUFfPpsLZHKNKNTgNPgqFgojiuHAwp6wB62vnDCytra1C/Zs979FoFLlcTnr2zTYstlRxqIrWWs6Zojlua1LPDocDU6dORX9/PwKBgCwCmIWafeykgukhz3s/ODgogZE2sD6fD9lsFi0tLWLaY6rOOeKTbmikv81r5vs7nU5p4+rp6ZFFWD6fRyKRQDQahcvlqjN6IcxShnkvnU5nXf3aVK4DEJvgTCZTt28wGBShH4WRrNOzPEPLVgCSZbO8YjrZUUBoDnWxlLqFxeTHFtXQlVJOACsAbAfgIq31cuPXewPo01q/vJHdNYC7lFIawCVa60trr/8OwJ9Rzew/b2yfRTWofxfAGZs7N4fDISr0E044QQI6ADz44IPSmmbivPPOg9PplAlr1113Ha677jrEYjHE43F86lOfkp5uBjH+25wgxiBrUuWpVErUw+wbZrY7ffp0ABCjF2ZbACTgMdCz3k4RmKmw93g8UpMNh8OSlZOi5/FcLtcGynxSzKyvMyhqrfH666/L3O/GQEW6loIyc6gKt2WmSu941ogLhQJWr14twZqZo8vlQjabFVEYA3E2m92g593MjJnJMns1R5FmMhnxYOfgEaWULBaYzScSCbjdbqTT6Q08yd1uN1KpFJLJpGS2XFSYDmpaa2mXY4Cnb31jTzwXh9QvmAsafuYE2+0YtKmNMO1vOaSG+7Otj46EVqluYbHtQZnins1urFQUwBIAJ2mtn6m99gcAr2itz9/IPl1a6zVKqTYAd9f2XbaRbX+KakC/HMCTAHYGsGZjNfSdd95Z/+Y3v5Gfmdm89dZbSKVSCAQC2HHHHTfYb8WKFWhqasIuu+yCgYEBvPXWWxLgvF6vCLjMHmKqkBvvF1u1KEJjEGILGrNDvsZs2swcGx/oADb4HRXXDKyNx+V7mUp3ADIqVWstv2OWzm3M+eiszXJBA6zP/Llw4TVQMc3gxkyTnwX/mOc6nlCPIjFd66GnnsC8FmaX5vx4HsPUClDoZvqecxveOzIlXIzw8zWvy3wvXgfPm9fNe8ayjcPhQCaTqfNC577ch/vxvvG9ANRR/OZnyNeoleB3whTRjfdeHwSy2ewGgrttDdv6PdjWrx947+/BwoULV1CPZuJtqdy11kml1FIABwN4RinVBOBIALttYp81tb/7lVJLAPwbqiK7zb3PtQC+s7lz2m+//eTfSikcccQRWLx4MY499lgsXrxYMkzigQcewA9+8APcc8892HPPPeHz+TBv3jxcf/31+M53voOjjz5alN3MyLXWQvGOlzFHIhG88MIL6OzsrHsNQB39Tkcvh8OBZDJZ97u2tjYZlZnJZOp+FwgERFxFOrXx/YPBIFauXCn0P3uVGbRIk/O4pVIJU6ZMwejoKAYGBkQT0NvbC5/Ph76+PhlewuzfzJp7enoQCoWEoh4cHMTLL7+M/fffv85AhRk03zMcDsPn88k41WQyKcyB1lomwr355puiMqeCPhgMCptAW1Vms/RgdzqdiEajknXz88pkMojFYtLqlkql6u4xFfONaOyLnzZtGoLBYN3CbHBwUNzhli5din322QfpdLquHOHxeNDT0yN0PJkRGhQlEgnx32cHBVXsPIe+vj453+7u7jpb2g8Lli5dWvd/clvEtn4PtvXrBz64e7AlKvdWAKVakPUBOADAf9d+fQCAF7TWqzaybwCAQ2udqf37QAA/38Jz+x8A/9zUOT777LPo6+tDe3s7brrpJgDAkiVLEAwGsWDBgnEf0H/4wx8AAHPmzMHcuXMBVKn6RYsWYebMmdh7771RqVQQiUTE1AMABgYGxE0tFovVtVYxsFGw5fV60dbWJi5tpu0r+9NzuZz0hTPj8nq9SCQSAOopeDrKUYnObJnDNABgzZo10q/O4M3a/9jYmCwKzBIClecdHR0ymIVmK9QImBmwScOnUilxruNYT601hoaGUCqVhO1g7z1b5mKxmJifmNQ6a+68l7FYTCaaqZojn9mbPTg4iLa2NhGgeTwecXrjgoLmOCxv0EeeCx6q751OJ9rb29HU1CQMxnhzw0dHR5FMJgGsn53OgMp7yrKDOdSErIRpAMPPlqwNB+wkEgnx4Dfpfd7rUqmESCQi98nCwsKC2JIMvRPAVbU6ugPAX7XWt9R+dyyA68yNa+1ol2utPwmgHcCS2oOnCcC1Wus7tuTEtNaDtYz+Pze2DTPKJ554Qqj1YDCIuXPn4v7779/AQQ4Arr32WnzsYx/DhRdeiClTpmDHHXfEN77xDTQ3NyMejyObzdYppfk+Zv28XC6jtbVVxq9Sybz99ttLz7rb7UYul0M8HhcjFirIGTR7e3vFUS6fz+Pll18W2pvtbpz7Teqc5i08L7Z1MXhThc5WLh5/1apVQs8ymDE45vN5acli9soASKU/W7DYGkfBlenwVi6XZWExffp0qUu3t7cjk8mgo6ND2th430zl/9jYmDAFbIvjBDqPx1MnJOP8d2a8rJczQNLlbWhoqE7g19LSglgsJhPcTFEdaf9sNivOffF4HJlMRr5bFCSOjIzIIqxcLiORSMjChb3equbsx8/ABMWADMpstzSV6jSP4djTlpYWq1S3sLDYKLZE5f40gPkb+d2Xx3ltDYBP1v79GoBdtvRktNY/bfj5+wC+v6l9CoUCZs+ejYMOOkiytMWLF+O2227DgQceWLctA/D8+fNx1VVXIZ/P48c//jFGR0exyy7V08zlcqImZpbJVrhcLifq4WKxiHw+LwGRaGlpQTKZRDablZ5gWosODAxg9erVdXVOU0zVmAEzWNHbvVgsIhgM1vWhO51OCUz8ma1pAMbdB4DU1Rlgc7mcBGWTYu7s7JQBLQ6HA4lEoq4/nHR24/kPDw9jbGwMbrdbglUul5P30TUTGVqS0oqVNDvLG+wlb7SNdTgcyGazUErVCdpYhx8vwy6Xy5Jls25OkRwXSeVyWYI9WwtbW1sxNjaGlpYW0VLQ7MbsCOA9IrWfyWTQ2dkp4kT6CQDVxZjJIFEQaH4mZHQY9B0Ohw3kFhYWG8WEdoqbOXMmurq6cP/99+POO++U11tbW8X5zQQHqcRiMfT396NcLmP33XeH2+1GMBiUfmxON2NwZLYLQHqCGwOF2XYFQHqtAWBoaAhjY2N1D2wKnFhTN8ddmn3LHK1J6jwWi2Ht2rV1ympOVONYUB7bpNuprOe2TU1NSCQS0vcdjUY3uK5SqSQe7FSsNzU1IZvNikiL/d+ck05bVJ/PJ5nseIsE2qaybECzm3Q6LYp1c556OBwWypmfl2nMQre5RCIhtDlb6DjUJRAIoK2tDQDE196k41mDnzlzZp35Chd4DKYMrI32txSx8XMEIP3uXq9XdBNcqFE1bw6UoWCO/u6WVrewsNhSTOiA/vrrr+OVV17BtddeizPOOAP5fB5r1qzBtddei+OPP36D7Ts7O9HR0YErrrhCgu1+++2Hxx9/HPPmzRMvcVLXHFzCmrfpUMY6OoVeTU1NdbS0SfN2dHSIT7qZYTJLZ2/7jBkzZMFAAxIAIgBjbZiUNxch6XRaxGL5fF58xUlpM4uNRqPiSMbXOLedznaxWEyyz1AohHg8LsG70V+8paVFBGjMeMPhsLShMUs1rUfZq26Of2XQow2ruegxh4dQmc+WNa/XC5/PJ+NYeb/j8bhk0wDQ3d2Nvr4+BAKBusErnCfeyC4MDQ0hEolIrzr9+pn9m1oDljcY3BmAzVY3AGKLSzEe++zJSIRCIXg8HsTjcUurW1hYvCNM6IBeqVSwfPlyfP7zn8fnP19tZd9tt93w5S9/GXfddReuueaaDfb50Y9+hJNOOkke3ttvvz3mzJkDl8uFSCQiIyVZh2V7FB/abHliFs2s+ZVXXpHMtTHLdblckimyLur3+5HJZPD6669LgOQ862KxiHXr1snrdAxjjzMDN7PjxoBULBYlaFKoRnEb1fpm/Twej4ttqlJV3/WhoSGh1k1zG/N9yGYAEOEX68Kc6MWM2rx/XBxRSMaAmcvlhBEx3eCSyaT4s5ugW17jQBfuQ1aBve98D2bzZB5I5VOASIqc5wlUqXiOviWr4ff7kUqlZAGTz+fR19cn52aK5vgdMhd7LpcLU6ZMqWtvsbS6hYXFO8WEDugOhwP77rsvTjvtNPz0pz+FUgoPP/wwPB4Prr32Wpx++umYM2dO3T7f+MY3sGzZMtxwww3Ye++9cd9998lD1xSjFQoFsfVkUKVAjMI0s+ZsCrnYj0w6Pp/Pw+fzSYYZCoXkOGYgouVoY4AmNc8s3DSYCYfDkukxGLLtjrVh7sOARvGXGQCpBne73UIfswTB86H/OCe3UVzHYMjFzJo1a8STnQJDurPRaKbx/lHpH4vFEAwGkUwmxdmOzm2kvlOplLAaLCFUKhVks1nxdyerQEEb34fe8e3t7YhGo7KwiMVicLvdiEQidYNNyMKYQXjatGnyOdGGlaUTTkxjDZz9+fl8vk75HgwGpcPBwsLCYmtgQgf0UCiE2bNn45xzzsE//vEP3HrrrYhEIvjDH/6A73znO3j88cc3COgejwd//etfxZ+8EaZFZmNQZVsTUN+fzPYzUtqVSgVdXV1iw0nlOEVkVD6bx3c4HGhpaYHX6xWqGoCYsdCxzJznTXq4VCqhu7tbps/x2lhW4LZOpxNtbW0oFosyvINZKR3PWI82qfWOjg6xQY1GoxgYGBAqnMYqmUwGlUql7p5Mnz5dXO3YR06RXCKRkGErSilMnToV4XBYFgZkK6gkZ2nBPC9zn0AggIGBAaHwCQoMG1kT06yG1Hjj5DPC3Nfj8ch3gOwE6+3MtDmOlaUMiv6CwaD4ExDW0c3CwmJrYUIH9FQqhUceeQRnnHEGLrroIvT29mKnnXbC8uXL4fP5cMABB2x03/GCuQlma8ViUfqFSXuPJx7jQ55173A4LHS9y+VCX1/fuO1TsVhMsvNMJiM09dSpUyWLHxkZwbp161AulxGPx4VK52hXLgBaW1tFIf7mm29K4Ovu7hZlNq1Z4/E4wuGwZLrBYFCYBTOrZ3ZJK1dm3dlsVsxu+Iftc9yPLAXvDTNYKtcb2Qlm8ASZBGoKGpkLKvnZh87RrY295CbDwcWT6SlPQx8OgCmVShKsOS6XiyK+PwAZmGJavDocjg0GwZjDWgAIC2Tr5BYWFlsTEzqg77zzzohGo/j5z3+O888/H6eccgrWrFmDf/u3f8PAwAA+85nP4LDDDsMXv/hFdHR0SOtUa2vrZo9t0u8MqAyQsVgMkUhEsnGq3IeHh6XtjdT0eONEWX9OpaoD6hjo6HbGjNEcg2pmpqzDNx6T5iONVD5nqgOQrJiZL1vK+vr6xImO/eh0qQuHw0J9U5E9PDwsA2c6OjpkAUNGAYC0iJGZYI2c/fP8mb7ndG6jaQ7b4UwlPd+fPeGcuMZgWS6XsWrVKpkCR70CaXQ6sHGhlMvl6lrPBgYGEI/HhUbnufT09Ig+gcwHFf0A5PM2PxMa7pjDWvjdsoHcwsJia2NCB3SXy4X77rsPn/jEJ/Dzn/8cp556KtasWYOrrrpKaOqf/exnOO200yR70lrjqquuEhHd5o4/Hv3OLJ1B1efzYc2aNWIKQuEcPcX9fr8EE9a5Wb9ldtuYsZZKJanHN/6Ok7nYjgVARHykwUmpOxwOqf83Zt7ZbFaugQHN6/WitbUVPp8PQ0NDUtMGIK15pgjNVJ0DQFtbm6jWOQTFVKC3t7fL8BQAMnOd9HQ6nZaFTTqdRltbGyKRiGTuHR0dchyPxyNMRTabrRP70abVnGvOY3DBQXc8nh/r6aFQSK6Zve78OZVKyWvFYlHKB/xcac3LlkJLqVtYWLxfmNABHagGhHvuuQcAcMkll0jGtGDBAjzwwAMAgJdeegn//Oc/obXGlVdeiVNOOWWLAjpB+p2CskAgILao5tCMRsEae7NJ5U6bNg0AJDNlGxUV0bQHpeKdLEA2m62b101av1AoiG0rM8e1a9eiVCqhvb1dzGiYqZrT4Zjp+3w+rFq1Ss69UCiIeYyZWbJOXi6XEQgEROzH7JmsAFvqeG8atQLmIollDPPnxu35h6wFWwrN9jBmxua+7EmnFS+1CSMjI5J5s75tqvBZWgDWD0ThbHagqtswyyepVAotLS0Ih8PS5tbd3W0pdQsLi/cdEz6ge71e/O///i8WLFggArizzjoLa9askW1mzZqFWbNmAQA+8pGPYN68eejv7xeTkc3BpN8ZOBmkmbFz7jWzc2bZbPdi7ZnnnEwmMTw8DJfLhZGREbjdbsycOVMMaJg106Ckt7dXqP1sNis18nw+L7PKzUx7dHQUU6dOBbA+y6xUKpgyZYqcExXc5XJZ2rsAiJVpuVyWujOvkdky6+8jIyOi4ue9Ms1SzD5tsg3AegEZzWjorMfhLAyIpkEN2+LK5TL6+/vlfE3BGQWHoVBIRIqNtHoikZDFTCAQkHtD4xqa3ZCJ4MKGiwCTVmdGbrrVWUrdwsLig8CED+gOh2MDE5kf//jHG91+9uzZCAaD+M1vfoOzzjpri9/HHEKSSCQkq2Wm63A40NbWVtfrPTg4KGpn9kEzc9Va16m8Ozs7xe+bQdH8/aZsYmmCY74+MjKC/v5+yWi5qKANLSl4agPYZhYOhwFgA/o6FosJzc/gyto0A3ulUkEikRCzFIrBnE4nRkdH0dLSIrVmKsG5+KE6na51/f39iEQiiEQikimzQ8A8rylTpojrHFD1jSeNz6A6Ojoq18l7RqMdAMKAtLS0IJ/Pi01roVAQ8WQ6nRahHAfdMHO3tLqFhcWHARM+oL9dUPV93nnnva2Abu7P2q1ZU6fICoBkvWbgAeqV9Y1mKKyTE41Bm3VcCs84GtWk51kXZ8sXW8bM8zDnbDcq9SnyI8vQKLgzPetN8FrHxsYwODiI3t5eyXwBiBELe9y9Xm8dbc0WNJ/PJ/75ZDt4z4Fqhjye0p1ZNNXu5XIZTqdTmBQGZtOwhosjALKgYDmlXC5LnzzLHGzvY4sf2wgtrW5hYfFhwTYX0DOZDF544QUAwPXXX49jjz32bR+DmSxr6pFIRFTYrNM2CtDMUaCsq6fTaRFT+f1+qbkDwODgoAwuYQAya9pdXV0SsEulEtauXYtyuYxwOCy2ssCGCwMzSzf9471eLzo6OuByuWQ6GNX3FPEFg8E6hoGUPxcnFB2mUilR0bP1jAHUdGszFxOpVEq0B1wwZbNZaRdjcG1ubhY6PxQKCWvAlj9e8+rVq6UeTm95ngNb+9h3Tg98AHJdjYsNU+gGWFrdwsLiw4dtLqAzmAPAiSee+I4CusvlQldXl7S0FYtF6WMeHR2VyWAMxlprUbNTnDU6Oiqe61TNp9NphMNhoa+nTp0qgjMuEoBqr7vZjz4yMlIXgNra2sTCtKmpSZzc+Mfj8WB4eBhutxvbb7+9DDThn2QyCafTienTp0NrjVwuh0KhIO1lDNwcM+rz+WS0KWl4ziL3+/11rXt01BvPM576AZNWHx4eRmtrq9TueX25XE4WIABksEuxWKxT7Pf09NRNz6NnejKZRDqdlkBulg7MrgJ63ZvB3MLCwuLDiG0qoP/617/G66+/LlnZ0NAQ+vr60N7e/raPxQyN7V2jo6N44403xAs9EokgGo0KncuMHIAEeA42Aapq/Ww2K5m2x+MR61CO2mTtmS1dnOvNY5psQCgUkiDKbJ4LB74H929tbUUymcTQ0JAI2eLxeF17VyKRQDqdRjQaRaVSEXaBqnaHw4E1a9agp6dHJpHl83kR/43n7tboGQ9ArGMb6X5mzlyQmE5/HGXbSMmThaAvAD83ivs4/Yz3iTqCeDxe91nbGrmFhcVEwDYV0H/84x+LMxpxySWX4Cc/+ck7Piap8Mb+btLZrOXSgcy0BaW6m9l7Y92dLWGEGaxoakJl9/DwsFDMzFIp4qLXOWvjjYsC1onN92arGzNyzpqna5zT6ZRWLQDS+mW2oI03OIbiPI/HI+1krO1zihmDLin4pqYmodT5fmbdHoAsYqhWZymhWCwKS0CxHrB+brppVdvT0yMZu1misNS6hYXFRMA2FdCPO+44XHXVVSJeA6p19HcT0AFI/zGDOGld0xa1ce62ORDGHPpiBj+v1yu1WzrJsT0NgAjkKpUKWlpaZBBKOp2uM7kBILSxGWCp0CbL0Ghsw2MD2ODcOzs74XQ661T4AKQEQAvWYDAoffekrxmcmWkPDg7KtDav14tYLIZAIIBkMom2tjaxX+Xsc3NYisPhkBo3afWuri5R1rMffWRkBF6vV8bSMthrrWWSmylKtDVyCwuLiYZtKqAfddRR+OMf/4j58+fj4IMPxi9/+Us8//zzdYHvnYAP/56enjprUrM3ebzWKZ/PJ7O+AUhAYlZs/pzNZtHU1ISZM2eKuQx7orPZrPSBZzKZDZzlPB4PotEoSqWSlADYuhUKhWQWOjNZ2qOa08I4OtTMtKkRYJtdpVLBunXrpJc8EAjIgsV0d2OvdyaTqVskTJ8+Xah3l8tVZzpDS1tzYhlV81zYFItFAOt96QOBABKJBFavXl3Xfsdzi0QidUHcUusWFhYTGdtUQD/ggAOwzz774B//+AeeeOIJ7Lfffli6dCnuv/9+7LPPPu/q2MzM/X6/1Gaj0aiYuOTzeXg8HgSDQaG+OdmM1HSlUkF3d3edLWoymURTUxNSqZT0dXNS2NDQkPSqp1IpxGIx6a+mYI0LA7avuVwuzJw5E4ODg6hUKhgaGgIAxGIxmUuezWbrZruTevZ6vSgUCvB6vWhvbxfb1Xw+j8HBQYyNjWFoaAhtbW0IBoNCb9MNzxSWjdeCxuDu8/nqMn5TGe90OqXFzJxpDqwflsL3Gc+pzuv1SrcBYAelWFhYTB5sUwEdAG666SZxKTv66KOxdOlSXHDBBe86oBNOp1MEV/39/RKwg8EggsGg0O8dHR3IZrNIJpNobm6WMgBnpfNYw8PDMped2fd4feSmWYrWGs3NzTKfPZ1O49VXX5VaMfdnax172UnDOxwOycg5fKW1tRXxeBzr1q2TeetANVhzohnV7zSlYR2eC5ZUKlXHPEQiEXGrI6NhDoApl8sYGBiQmn0oFILb7cbQ0BCmTJkis8cLhYJoDmjryna6UqkkDnCRSAShUKiOjbHUuoWFxWTBpA7ot9xyCz796U9j8eLFOPzww6GUwq9//WtMmzYN69atw9VXXw0AuPPOO7fq+5K+Nie00UOdwQlY7zg2MDAg2SOd35h9NvZrs4/c4XCgUChIrzpHgtLClQGQYrxGurzRe54qcQbkxqyW5+33++sybfq7s3fc5XKhra1NxrUyszap9e7ubkQiEVkApNNpdHR0iK1rozK+p6dHAjEXPmQAcrkccrmcCBE5AhaAZOOtra02C7ewsJj0mNQBfe3atQCAI488Er29vbjmmmswb948vPHGG+jt7cXy5cvhdDoxMjKCV155Bdttt91We29mrBRnUVxm/huoF5vR/Y21YQYfBm1gfW0ZgIyCpYPcyMiIBFhS+PQ/p+o7HA6jvb1dLFrpN9/W1ibTxmhqQ1c1KtO5UKB4L5vNYmRkBOFwGK2trXjzzTfr+sV9Ph+y2SwSicQGgjueG9Xu5uz3Rg0Aa/UAxO+dLAMXRbyHmUwGXV1ddfPoed8sLCwsJjMmTUC/55570NvbixkzZqBUKuHMM8/ESSedhGOOOQatra1YuXIlPv7xj2PGjBkAgJUrV4rI7MUXX8Txxx+PBx98cKudDzNp1oo5QITiLr4+MDCwUY9xDmLp6uqC2+1GsVhEsVhEJpOBUgqZTAbRaFRoZtLf42XDwWAQw8PDiMfjosbv7e2tM8dxu92SvdPUBqja1CaTSaltU4jHyWUsYVD1TmMaLl78fj8ymYwwDFzEcIFATwAG+WAwKIsYuro1ivmA9f3nnHbHdjge18LCwmJbwlYL6EopL4BlADy1496otT6j9ruTAJwIoAzgVq31qbXXzwWwEMD/01rfp5TqBfA6gJO11r+tbfM7AI9prf93U+///e9/H0899RSOOOIIXHLJJTjzzDNx5pln4oQTThD1cqlUwrp162SfcrmMl19+GQDw0EMPYd26dejo6Ngq98Oc0MYpY8lkEtFoVOxVWQ/3er1SWy4Wi0LNs62K7WiBQAC5XA7pdFqCdjablXo5gA1EYLlcDh6PB36/X1TxpJ8DgYAo0dPptNDunLdOMNPm/jwG3ejY7272ddPRjXPflVIoFApoa2uTcgGzce7DEbMulwuxWExEhJyW5nQ6USgU4HK5pJ+ftXgTVq1uYWGxLWJrZuhFAPtrrbNKKReAB5RStwPwAVgEYGetdVEp1QYASqkda/vtA+B/AdxX+7kfwHeVUpdorUe39M3vu+8+HHTQQViyZAmWLFkiQqg//OEPss1PfvITnHHGGXj++edx22234Yc//CEqlQoWLlyIe++9F1/60pdw1113vdv7IDAFV263G8lkEgMDA3VUeGtrq9R4SW8nEglEo1Ghl7kAYC97oxiOJi+sl1N8FggE6gxhaONaLpfh8/nEsrZxgAtNXVhnJxpr4bFYTBzYgA0XE1Tj8/pN1rKDiQAAFSFJREFURTrb10x6nZQ/92WXQFNTEwYGBkSol06n0draKsY29MK3dXILC4ttGVstoOvqkz9b+9FV+6MBnADgHK11sbZdf20bJ4BKbRuzCXwAwIMAjgdw2abe03R9c7vdeOaZZ+TnadOmYWhoCD09Pbj88ssxd+5c6dWeO3cu5s6di3A4jG9/+9t47rnnAAB33303Hn74Yey1117v5BZsElScU7GulEI8HpcBJOyl9vl8yOfzopB3uVxIJBISOJ1Op9SouQigZSkz31gshlAoJH3cFM9tzCWuUQRHS1fT/10phcHBwbrttNaIRqPy+8bZ51S7syTA3nn2ofPcG1vQcrkcmpqahMKnUJDv7XQ6EQwG5fMErFrdwsLCYqvW0JVSTgArAGwH4CKt9XKl1CwAeyulzgZQAPADrfU/tdbPKqX8AB4AcErDoc4BcLtS6opNvd8LL7wgw1Vuv/125HI5PProo5g/f36dZerG8K1vfQu77rorQqEQdt11V5RKJRx22GHo7+9/V0YzG8PY2Bg8Ho9kxXRPI4XNbLgxE25ra5Nsm8rusbExhMNhlEolcU7zer1IpVIiNKMNLRcDjT7nNJYh/U+BHLNi0xgnk8mI29zo6KgE71wuJ4p1pRRaWlrgcDikfS6fzyORSNQtJOLxOILBILxeLxKJRF0LGnvgSbVzUI3X60U+n5eAbml1CwsLi3ps1YCutR4DsKtSKgpgiVJqXu09mgHsCWABgL8qpWboKk7ayHFeV0o9CuDzm3q/qVOn4oEHHkCpVEI0GsWxxx6L1atXo6enBxdccAFuuukmrFy5UkRo06ZNw+23345Zs2YBAJYsWYJKpYKjjjoKO+64I/71r39hcHAQv/rVr3DaaadtvRtTA0VbVIkzs6Y7Gs+z0VGOlLI5EpX1ZQBib7py5UpZBJiUdyAQgNvtFr93zj03XekYlKPRqARWLkA4GzwcDiMcDmN4eFgCdyaTQSKRkDq2w+GA3++XgNvU1CR96ryeYDCIUCgkDAt774eHh8WylgyC3++X4zU3N1ta3cLCwmIjUGaNdKseWKkzAOQAHIAq5b609vqrAPbUWg+Ms08vgFu01vNqNfYbURXaPTqeKC4ajeqf/OQn2Ng1BAIBtLS0oFKpIJVKIZ1OIxQKYerUqXj99delpzkWi0kfM1vK5s+f/7ZGZjKr3BzMqWEcEOJwOOS9GGw5TIXtbwy+3AaA/Ky1lto6fd/N+jW3Yb86ALlOZtYU4nHRwMUF6+Ks+fP9uS2PC1THunLuOM+Dixceh/7x7Ms3r6NcLssQFfa0TzS1+pZ+DyYrtvXrB+w92NavH3jv78HChQtXaK13b3x9a6rcWwGUtNZJpZQP1UD+36jW1fcHsLRGv7sBDG7ueFrrF5RSzwH4FIBHx9smFouhpaUFixYtQiQSQaVSwdVXX41oNIqzzz4bv/vd77BgwQIAwKWXXoof/OAH6O7uliwQAPbff3/ce++9mD17Np577jnJYA855BDcdtttW3z9S5cuxX777bfF2xMUrJXLZQmUtGFlNjo2NiatahxlygEwHIKitZYWOJfLJYKxYDAIt9uNTCaDQqGAeDxeZ+hSLBZlMRCLxaQ3nYK1VColgjX2qjMjp1vc2NgYUqkUHn/8ccyZM0fc8ajgZ4udz+eTgM3WtGAwKP3nZBoo6Ovs7Jxwmfg7/R5MFmzr1w/Ye7CtXz/wwd2DrUm5dwK4qlZHdwD4q9b6FqWUG8AVSqlnAIwCOF5vOS1wNoAnNvbLWCyGL33pS/Kzw+GQnw877LC6be+//34AQHt7O7bbbjucc845uOyyy3DBBRdg7733xtDQEJqbm5FIJABUa/K0DX0vEQwGUS6XxZOdKnPWqDkshdQ16XhS0j6fTyxRTfU5e9jZE57L5eD1ekUtT2ZgPKEce7k595zq+o19bKTZtdZyLLfbjWnTpsHj8WBoaEgG0ZTLZampm3PIraubhYWFxbvD1lS5Pw1g/jivjwL44hYeYyWAecbPT6G6OHhXePLJJ6XfPJvNYvvtt8eee+6JPffcEwDwyU9+EhdccAEuvvhiLF68GNtttx0uuOACXHTRRfje9773bt9+s3C5XHU1ZlLmpOdJVze2q7GuDEDodJrX8LhDQ0OirKfSnYJBjkc1rWVDoRBCoRCAKk3f0tKCQCBQNxK1XC6LAQ3V6zyO2YI2NjYmU9h4TC5SzMDf09OzgWLdwsLCwuLtYdI4xW0KSiksX74cQHUgSWtra93vGUCOO+44fPnLX4bWGq+++ir+8z//E4cffjh6e3vf0/Nzu92IRCLiHsdzMlXmzI6ZMVNtnkqlJMvnfgykNG9p7FmnEY3D4ZAsnvVs+r4z8JIBoMVqNpvdoBfd7/cLo8CFAal/quObmprqzGm01mJZa1rhWlhYWFi8M7zr7Hci4Nlnn5Xe8v7+/g0mqz300EMoFot4/vnnAVQXAMzMf/rTn77n50dXuba2Nmy33Xbw+/3IZrNIp9OS9abTaVQqFfT09GDGjBnYfvvtEYlExBRmaGgIhUIBqVQKQ0NDyOVycDgcaG5uFpOdpqYmEbpRHKe1Rjwex7Rp0zBlyhSUy2XJ5CuVSp0QjhS9mYUHg0FEo9E6O9nOzk6Ew2GkUinkcjnJyDlCtrW1FZFIRARvtgXNwsLC4t1jm8jQb7zxRjz88MPy89FHH41KpYJ8Pg+fz4dly5Zh+vTpdcNZ9t13XwDAHXfc8b6cY6OrXD6flz5sCsmYCXOSGZXwptNbOp2W4SSk7c2WtHw+LwsFn88nFq5sncvn88hms+Io19TUJJawAMQIhu1nPp9PLGs5bMXtdqOvr69uBKt5ToCdQ25hYWGxtbHNBPTnnnsOjz/+OObNmwelFE488UT8/ve/x4IFC9DS0oLtt99eppgBkAy2r6/vfT9fpRQCgQC01nW95UB1Whu3oQEMPd9N/3aXyyWqcvbpezwe8W3PZDIYGRlBLpcTExhm2Y1Cufb2doRCIfGOnzp1KjKZTN3xKOYrl8sYGRmpc3Zj/71p9mOd3SwsLCy2LiZ0QB8dHcXdd9+NefPmiVnJeHA4HJg3bx7mzavq7ZLJJG6//XYAwD//+U84HA4sXbp0g/1aW1sxMLBBu/z7BgZt1svL5TKam5vh8XjqeuTZf87sOBAIoFAo1I0VZc89+9NpL6trI0ebm5sl4Jq0Omv5Xq8Xw8PDaGpqQjgclla4XC6HQqEgqvjh4WEopaSuz351S6tbWFhYvLeY0AH9X//6Fw488ED5ec2aNejs7NzkPo899hj23ntvuFwu3HnnnSgWi9hll13Q09OzwbZDQ0Nb/ZzfDkwXObrKcQxpIBCoG1E6ZcoUEZ8BGFcQl0wmhYUwg/bY2BhcLhdCoRBKpZIEYrrNsQQwNjYm3u0c1ZpKpermr/v9fvh8PjidToRCIUurW1hYWLxPmNABnfD5fDjvvPPwwAMPYMGCBRtVpedyOey1116IRCJ45ZVXZO74xvBBq6/NEaz0MM/n8xvMLweqrnilUkl+pl0sVfGhUAhtbW0YGxuT4SeshUciEUSjUWlBo8Nef3+/OL9R0MZ+92KxKMNXtNbweDwyNMU8fxvILSwsLN4fTIqAPjIygueeew4XXXQRdt55Z1x11VXYdddd67Z54IEHsPfeewMAVq5cuVlbPtaqG1vc3m80BkUONUkmk3jrrbfq6uter1d8z10uF4rForSVMcsul8twu93yOmvkrJfTZU4phUgkgkKhID3nPp+vbmAKFfixWEx65S0sLCwsPhhM+La1k06qzne56KKLAABPP/005s+fj0WLFuH73/8+nn/+eZx//vnYf//9AVRV61visfvII48AQB2l/2FBU1MTfD6fUOYUtHF2OnvCi8Wi+KFns1lUKhWEw2FEo1E0NTXJVDYASKVSaGpqQiQSkQUCe8o5nIVtaNyO7WdUwL8XE+osLCwsLLYMEzpDdzqdyGQyuOuuuzA4OIh7770Xl112GXbaaSfcfPPNAIBf//rXUErhmGOOwX/9139hzpw5W3Tsv/71rwCAT33qU+/Z+b8b0IyG/2b9m+1jHG3K1jW2ppliOira+/r64Ha7pZ+cveompT5e6xlQ33722muvve/3wcLCwsKiigkd0OPxOBYvXoy//OUvWLRoEa699lqcc845aG5uliCUz+fhcDjqrEU3h5dffhmXX345AODwww9/j87+3aGxvk77WM5Ap90r3dxGRkZQKBSkvUwphZGRkbqATctZYL1FK1kAj8ezQesZz8PWyS0sLCw+eEzogD516lSEw2E89NBDuPXWW6WeC6ynf6kS3xxGR0fxpS99CTfeeKNQyD09PXW96R82NAZTTmMrl8syiIXtbB6PR2re/GMq3QHA4/HIbHL2lNNSFrCObhYWFhYfZkzogP7cc89Jf/Xf/va3tzW/HABWrFiBk08+GW+++SZWr14toq5p06Zht912w6WXXrrVz/n9AM1ciHA4PG4W3bjYadxuPBbAwsLCwuLDiQkd0BnMH3nkEeyxxx5btE+pVMLPfvYznH322fJaOBzGqaeein//93/HwoUL3/bC4MOGLQ3EW7KdpdQtLCwsJgYmdEAnZsyYAQB1LVyN4KjPGTNmYNWqVQCAr33ta/jZz36G7u7u9+1c3y9saSC2AdvCwsJicmBip6IAnnnmGcTjcey33374xje+AQDj9kN/8YtfhMvlwqpVq3DppZeiUqng8ssvn5TB3MLCwsJi28OEDuhz5szBCy+8gP/4j//Afffdh2XLluG2226Dw+HA4sWLZbs99tgD1157Lb75zW8ik8ngG9/4hu2ZtrCwsLCYVJjQlPtzzz2Ho48+Wn5++eWXJUs/6qijcPvtt6O3txdTpkzBueeeu8EcdAsLCwsLi8mCCR3QOzs7sXbtWpx11ll48MEH8eSTT2LNmjUAgB122AGHHHKIDBGxsLCwsLCYzJjQAb2rq0sCuNYas2bNAlD1NH/xxRcBVLN2CwsLCwuLyY4JXUMfHBxEsVjE6OgoPv3pT+OVV14BUB0d2t3djVKphClTpnzAZ2lhYWFhYfHeY0Jn6G+88QZ23HFHNDU14Y033pDXjzjiCPztb3+zwjcLCwsLi20GEzpDB6qjUPv7+xGPxwFUvdcXL15sg7mFhYWFxTYFNZFnWCulHvugz8HCwsLCwuJ9xqDW+uDGFyd0QLewsLCwsLCoYsJT7hYWFhYWFhY2oFtYWFhYWEwK2IBuYWFhYWExCWADuoWFhYWFxSSADejjQCm1g1LqSeNPWin1PaXUT5VSq43XP2nsc65S6jGl1L61n5copQ43fv+iUurHxs9/U0od+b5e2NvAJu7BuUqpF5RST9euMWrsM9nuQVQpdWPtep9XSu2llNpFKfWwUupfSqm/K6XCxvaT6vqB8e9B7fWTatfzrFLqV8b228Q9UEqdWfs/8KRS6i6lVJex/aS6Bxu5/phS6m6l1Mu1v5uN7SfV9QOAUmpl7f/8k+yuUkr9xXg+rlRKPWls/4HcAxvQx4HW+kWt9a5a610B7AYgD2BJ7dcX8Hda69sAQCm1Y+13+wD4j9q/HwLw0drv4wCyAPYy3mav2jYfSmziHtwNYJ7WemcALwH4/4DJeQ8AXAjgDq31jgB2AfA8gMsB/FBrvROq9+MUYNJePzDOPVBKLQSwCMDOWuu5AM4Dtq17AOBcrfXOtf8ftwD4CTBp78F41/9DAP/QWm8P4B+1nyfr9RMLa8/E3QFAa/1Z4xn5NwCLgQ/2HtiAvnl8AsCrWus3NrGNE0AFgAZAR5sHUfsAa3/fAqBVVTEdwIjWet17dM5bG3IPtNZ3aa3LtdcfAUBv3Ul1D2qZ9z4A/ggAWutRrXUSwA4AltU2uxvAUbV/T6rrBzZ5D04AcI7Wulh7vb+2yzZzD7TWaWOzAKrXDEyye7CJ78AiAFfVNrsKwOG1f0+q698SKKUUgGMAXFd76QO7Bzagbx7HYv0HBQAn1qi2K0gzaa2fBeAH8ACAP9S2WwFgnlLKjeoH+DCAFwHMrv384Pt0/lsDjfeA+CqA24FJeQ9mABgAcKVS6gml1OVKqQCAZwAcVtvmMwCmApPy+oGN34NZAPZWSi1XSt2nlFoAbHP3AEqps5VSbwH4AmoZ+iS8Bxu7/nat9VoAqP3dVvv3ZLt+QgO4Sym1Qin1zYbf7Q2gT2v9MvDB3gMb0DeB2s0/DMANtZf+AGAmgF0BrAVwPrfVWp+ktd5Na31P7ecigGcBfATAngCWo/ohfrT2ZyJQTOPdA75+OoAygGv42iS7B02onvcftNbzAeRQpRW/CuA/lFIrAIQAjHKHSXb9wMbvQROAZlSv6RQAf61lKdvSPYDW+nSt9VRU/w+cyB0m2T3Y6PVvDJPs+omPaa0/AuAQVP//72P87nNoSHg+qHtgA/qmcQiAx7XWfQCgte7TWo9prSsALgPwb5vZ/yFU6aqQ1jqBKkXND3AirEqBhnsAAEqp4wF8CsAX9OatBifqPVgFYJXWennt5xsBfERr/YLW+kCt9W6o/id+dTPHmajXD2zkHtReX6yreBRVerFlE8eZjPfAxLVYX3rZGCbqPdjY9fcppToBoPZ3/0b2Jybq9QMAtNZran/3o6qd+TcAUEo1ATgSwF+24DDv+T2wAX3TqFt58QtcwxGo0q+bwoMAvgXgqdrPT6O6OutBdbU2EdB4Dw4GcBqAw7TW+S3Yf0Leg1o96y2l1A61lz4B4DmlVBsAKKUcAH4M4OLNHGpCXj+w8XsA4P8A7A8ASqlZANwABjdxqEl3D5RS2xubHQbghc0cakLeg018B24GcHztteMB3LSZQ03I6wcApVRAKRXivwEciPXP/gMAvKC1XrUFh3rP78GEHp/6XkIp5Qfw76h+AMSvlFK7olpPWdnwu/HwEKo1qF8CgNa6rJTqB/BWLcv/UGMj9+B3ADwA7q6xrI9orb+9icNM5HtwEoBramWH1wB8BcCXlFJUri4GcOVmjjGRrx8Y/x7kAFyhlHoG1ZLD8ZthaibjPbi8FuQqAN4AsKn/A8DEvgfjXb8D1VLL1wC8iaqeZFOYyNffDmBJ7XnXBOBarfUdtd9tTF80Ht7ze2CHs1hYWFhYWEwCWMrdwsLCwsJiEsAGdAsLCwsLi0kAG9AtLCwsLCwmAWxAt7CwsLCwmASwAd3CwsLCwmISwAZ0CwsLCwuLSQAb0C0sLCwsLCYB/n/HpKWKnNdqdQAAAABJRU5ErkJggg==\n",
"text/plain": "<Figure size 576x432 with 1 Axes>"
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
]
},
{
"metadata": {},
"id": "63b192dc",
"cell_type": "markdown",
"source": "# xugrid\n\nhttps://deltares.github.io/xugrid/"
},
{
"metadata": {
"trusted": true
},
"id": "60289131",
"cell_type": "code",
"source": "import xugrid\n\n\nds = xugrid.open_dataset(url, drop_variables=\"siglay\")",
"execution_count": 6,
"outputs": [
{
"ename": "ValueError",
"evalue": "face_dimension from connectivity does not match dimension lonc: three versus nele",
"output_type": "error",
"traceback": [
"\u001b[0;31m---------------------------------------------------------------------------\u001b[0m",
"\u001b[0;31mValueError\u001b[0m Traceback (most recent call last)",
"Input \u001b[0;32mIn [6]\u001b[0m, in \u001b[0;36m<cell line: 4>\u001b[0;34m()\u001b[0m\n\u001b[1;32m 1\u001b[0m \u001b[38;5;28;01mimport\u001b[39;00m \u001b[38;5;21;01mxugrid\u001b[39;00m\n\u001b[0;32m----> 4\u001b[0m ds \u001b[38;5;241m=\u001b[39m \u001b[43mxugrid\u001b[49m\u001b[38;5;241;43m.\u001b[39;49m\u001b[43mopen_dataset\u001b[49m\u001b[43m(\u001b[49m\u001b[43murl\u001b[49m\u001b[43m,\u001b[49m\u001b[43m \u001b[49m\u001b[43mdrop_variables\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[38;5;124;43m\"\u001b[39;49m\u001b[38;5;124;43msiglay\u001b[39;49m\u001b[38;5;124;43m\"\u001b[39;49m\u001b[43m)\u001b[49m\n",
"File \u001b[0;32m~/micromamba/envs/UXARRAY/lib/python3.10/site-packages/xugrid/ugrid_dataset.py:722\u001b[0m, in \u001b[0;36mopen_dataset\u001b[0;34m(*args, **kwargs)\u001b[0m\n\u001b[1;32m 720\u001b[0m \u001b[38;5;28;01mdef\u001b[39;00m \u001b[38;5;21mopen_dataset\u001b[39m(\u001b[38;5;241m*\u001b[39margs, \u001b[38;5;241m*\u001b[39m\u001b[38;5;241m*\u001b[39mkwargs):\n\u001b[1;32m 721\u001b[0m ds \u001b[38;5;241m=\u001b[39m xr\u001b[38;5;241m.\u001b[39mopen_dataset(\u001b[38;5;241m*\u001b[39margs, \u001b[38;5;241m*\u001b[39m\u001b[38;5;241m*\u001b[39mkwargs)\n\u001b[0;32m--> 722\u001b[0m \u001b[38;5;28;01mreturn\u001b[39;00m \u001b[43mUgridDataset\u001b[49m\u001b[43m(\u001b[49m\u001b[43mds\u001b[49m\u001b[43m)\u001b[49m\n",
"File \u001b[0;32m~/micromamba/envs/UXARRAY/lib/python3.10/site-packages/xugrid/ugrid_dataset.py:165\u001b[0m, in \u001b[0;36mUgridDataset.__init__\u001b[0;34m(self, obj, grid)\u001b[0m\n\u001b[1;32m 163\u001b[0m \u001b[38;5;28;01mraise\u001b[39;00m \u001b[38;5;167;01mValueError\u001b[39;00m(\u001b[38;5;124m\"\u001b[39m\u001b[38;5;124mAt least either obj or grid is required\u001b[39m\u001b[38;5;124m\"\u001b[39m)\n\u001b[1;32m 164\u001b[0m \u001b[38;5;66;03m# TODO: find out whether it's 1D or 2D topology\u001b[39;00m\n\u001b[0;32m--> 165\u001b[0m \u001b[38;5;28mself\u001b[39m\u001b[38;5;241m.\u001b[39mgrid \u001b[38;5;241m=\u001b[39m \u001b[43mUgrid2d\u001b[49m\u001b[38;5;241;43m.\u001b[39;49m\u001b[43mfrom_dataset\u001b[49m\u001b[43m(\u001b[49m\u001b[43mobj\u001b[49m\u001b[43m)\u001b[49m\n\u001b[1;32m 166\u001b[0m \u001b[38;5;28;01melse\u001b[39;00m:\n\u001b[1;32m 167\u001b[0m \u001b[38;5;28mself\u001b[39m\u001b[38;5;241m.\u001b[39mgrid \u001b[38;5;241m=\u001b[39m grid\n",
"File \u001b[0;32m~/micromamba/envs/UXARRAY/lib/python3.10/site-packages/xugrid/ugrid/ugrid2d.py:177\u001b[0m, in \u001b[0;36mUgrid2d.from_dataset\u001b[0;34m(dataset)\u001b[0m\n\u001b[1;32m 175\u001b[0m \u001b[38;5;66;03m# Collect names\u001b[39;00m\n\u001b[1;32m 176\u001b[0m mesh_topology \u001b[38;5;241m=\u001b[39m ugrid_io\u001b[38;5;241m.\u001b[39mget_topology_variable(ds)\n\u001b[0;32m--> 177\u001b[0m ugrid_roles \u001b[38;5;241m=\u001b[39m \u001b[43mugrid_io\u001b[49m\u001b[38;5;241;43m.\u001b[39;49m\u001b[43mget_ugrid2d_variables\u001b[49m\u001b[43m(\u001b[49m\u001b[43mdataset\u001b[49m\u001b[43m,\u001b[49m\u001b[43m \u001b[49m\u001b[43mmesh_topology\u001b[49m\u001b[43m)\u001b[49m\n\u001b[1;32m 178\u001b[0m \u001b[38;5;66;03m# Rename the arrays to their standard UGRID roles\u001b[39;00m\n\u001b[1;32m 179\u001b[0m topology_ds \u001b[38;5;241m=\u001b[39m ds[\u001b[38;5;28mset\u001b[39m(ugrid_roles\u001b[38;5;241m.\u001b[39mkeys())]\n",
"File \u001b[0;32m~/micromamba/envs/UXARRAY/lib/python3.10/site-packages/xugrid/ugrid/ugrid_io.py:194\u001b[0m, in \u001b[0;36mget_ugrid2d_variables\u001b[0;34m(dataset, mesh_topology)\u001b[0m\n\u001b[1;32m 193\u001b[0m \u001b[38;5;28;01mdef\u001b[39;00m \u001b[38;5;21mget_ugrid2d_variables\u001b[39m(dataset, mesh_topology) \u001b[38;5;241m-\u001b[39m\u001b[38;5;241m>\u001b[39m Dict[\u001b[38;5;28mstr\u001b[39m, \u001b[38;5;28mstr\u001b[39m]:\n\u001b[0;32m--> 194\u001b[0m \u001b[38;5;28;01mreturn\u001b[39;00m \u001b[43m_extract_topology_variables\u001b[49m\u001b[43m(\u001b[49m\n\u001b[1;32m 195\u001b[0m \u001b[43m \u001b[49m\u001b[43mdataset\u001b[49m\u001b[43m,\u001b[49m\u001b[43m \u001b[49m\u001b[43mmesh_topology\u001b[49m\u001b[43m,\u001b[49m\u001b[43m \u001b[49m\u001b[43mUGRID2D_TOPOLOGY_VARIABLES\u001b[49m\n\u001b[1;32m 196\u001b[0m \u001b[43m \u001b[49m\u001b[43m)\u001b[49m\n",
"File \u001b[0;32m~/micromamba/envs/UXARRAY/lib/python3.10/site-packages/xugrid/ugrid/ugrid_io.py:176\u001b[0m, in \u001b[0;36m_extract_topology_variables\u001b[0;34m(dataset, mesh_topology, ugrid_attrs)\u001b[0m\n\u001b[1;32m 174\u001b[0m \u001b[38;5;28;01mif\u001b[39;00m dimrole \u001b[38;5;129;01min\u001b[39;00m attrs:\n\u001b[1;32m 175\u001b[0m \u001b[38;5;28;01mif\u001b[39;00m attrs[dimrole] \u001b[38;5;241m!=\u001b[39m dimname_x:\n\u001b[0;32m--> 176\u001b[0m \u001b[38;5;28;01mraise\u001b[39;00m \u001b[38;5;167;01mValueError\u001b[39;00m(\n\u001b[1;32m 177\u001b[0m \u001b[38;5;124mf\u001b[39m\u001b[38;5;124m\"\u001b[39m\u001b[38;5;132;01m{\u001b[39;00mdimrole\u001b[38;5;132;01m}\u001b[39;00m\u001b[38;5;124m from connectivity does not match dimension \u001b[39m\u001b[38;5;124m\"\u001b[39m\n\u001b[1;32m 178\u001b[0m \u001b[38;5;124mf\u001b[39m\u001b[38;5;124m\"\u001b[39m\u001b[38;5;132;01m{\u001b[39;00mname_x\u001b[38;5;132;01m}\u001b[39;00m\u001b[38;5;124m: \u001b[39m\u001b[38;5;132;01m{\u001b[39;00mattrs[dimrole]\u001b[38;5;132;01m}\u001b[39;00m\u001b[38;5;124m versus \u001b[39m\u001b[38;5;132;01m{\u001b[39;00mdimname_x\u001b[38;5;132;01m}\u001b[39;00m\u001b[38;5;124m\"\u001b[39m\n\u001b[1;32m 179\u001b[0m )\n\u001b[1;32m 180\u001b[0m \u001b[38;5;28;01melse\u001b[39;00m:\n\u001b[1;32m 181\u001b[0m variables[dimname_x] \u001b[38;5;241m=\u001b[39m dimrole\n",
"\u001b[0;31mValueError\u001b[0m: face_dimension from connectivity does not match dimension lonc: three versus nele"
]
}
]
}
],
"metadata": {
"gist": {
"id": "",
"data": {
"description": "pyugrid_uxarray_comparison",
"public": true
}
},
"kernelspec": {
"name": "python3",
"display_name": "Python 3 (ipykernel)",
"language": "python"
},
"language_info": {
"name": "python",
"version": "3.10.5",
"mimetype": "text/x-python",
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"pygments_lexer": "ipython3",
"nbconvert_exporter": "python",
"file_extension": ".py"
}
},
"nbformat": 4,
"nbformat_minor": 5
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment