Last active
August 29, 2015 14:08
-
-
Save rabernat/dae4ebccaabbd72c424c to your computer and use it in GitHub Desktop.
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
{ | |
"cells": [ | |
{ | |
"metadata": {}, | |
"cell_type": "markdown", | |
"source": "# Surface Water Mass Transformation#\n" | |
}, | |
{ | |
"metadata": { | |
"collapsed": false, | |
"trusted": true | |
}, | |
"cell_type": "code", | |
"source": "import numpy as np\nfrom matplotlib import pyplot as plt\nfrom data_helper import load_data\nimport brewer2mpl\n%matplotlib inline\n\nSv = 1e6\nbmap = brewer2mpl.get_map('Set1', 'Qualitative', 4)\nplt.rcParams['axes.color_cycle'] = bmap.mpl_colors\nplt.rcParams['font.size'] = 11\nplt.rcParams['legend.fontsize'] = 10", | |
"execution_count": 1, | |
"outputs": [] | |
}, | |
{ | |
"metadata": {}, | |
"cell_type": "markdown", | |
"source": "## Load pre-computed data for plotting ##" | |
}, | |
{ | |
"metadata": { | |
"collapsed": false, | |
"trusted": true | |
}, | |
"cell_type": "code", | |
"source": "runs = ['control', 'perturb']\n\n# average results from each ensenble member\nwmt = dict()\nrholevs = dict()\nfor r in runs:\n wmt[r] = dict()\n rholevs[r] = dict()\n res = load_data('WMT_stu_%s_sigma0_hconst050' % r, year0=55)\n for reg in res.keys():\n # keep only heat flux and salt flux terms\n wmt[r][reg] = res[reg]\n rholevs[r][reg] = res[reg]['rholevs']\n \n\nyear = res[reg]['year']", | |
"execution_count": 5, | |
"outputs": [] | |
}, | |
{ | |
"metadata": {}, | |
"cell_type": "markdown", | |
"source": "## Southern Ocean Transformation Rates ##" | |
}, | |
{ | |
"metadata": { | |
"trusted": true, | |
"collapsed": false | |
}, | |
"cell_type": "code", | |
"source": "reg = 'so'\n#Alevs = np.arange(-40,40,5)+2.5\nAlevs = np.arange(-5,5,0.5)+0.25\n\nfig = plt.figure(figsize=(12,5))\n\nfor n,expt in enumerate(runs):\n \n wmt_mean = -wmt[expt][reg]['annmean']/Sv\n\n ls = '-' if n==0 else '--'\n\n rho = rholevs[expt][reg]\n\n ax = fig.add_subplot(111)\n ax.plot(rho, wmt_mean[0], ls, rho, wmt_mean[1], ls, linewidth=2)\n ax.plot(rho, wmt_mean[:2].sum(axis=0), ls, color='k', linewidth=2)\n ax.set_color_cycle(plt.rcParams['axes.color_cycle'])\n\nax.set_xlim([24,28.2])\nax.grid()\nax.set_xlabel(r'$\\sigma_0$')\n\nax.legend([r'$\\Omega_{HF}$ (control)', r'$\\Omega_{FWF}$ (control)',\n r'$\\Omega$ (cotrol)', r'$\\Omega_{HF}$ (perturb)',\n r'$\\Omega_{FWF}$ (perturb)', r'$\\Omega$ (perturb)',], loc='lower left')\nax.set_ylabel('Sv')\n\nplt.title('Southern Ocean (South of 30S) - Annual Mean Surface Transformation (Years 56-65)')", | |
"execution_count": 29, | |
"outputs": [ | |
{ | |
"output_type": "execute_result", | |
"data": { | |
"text/plain": "<matplotlib.text.Text at 0x10dd63f50>" | |
}, | |
"metadata": {}, | |
"execution_count": 29 | |
}, | |
{ | |
"output_type": "display_data", | |
"data": { | |
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAtQAAAFUCAYAAAAEdau1AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzsnWeYHMXRgN/dvaisUx41KCGByFnkDBbZgOHDgEnGgRxM\nMAYLIRsnwIBNMMEgcjA5ypggMkeOkpCEBKI1QjlLp0v7/agZdnZv925vw+3uqd7nmedudmZ6emp6\nuqurq6tBURRFURRFURRFURRFURRFURRFURRFURRFURRFURRFURRFURRFURRFURRFURRFURRFURRF\nURRFURRFURSlJBkPvFroTHQSzgImFjoTSB7uynGaPYBbgfnARzlOuz1MBg4q4P0LwWTgikJnokjZ\nAngcWAP8vcB5KQZ+DLwBNAGHFTgv6dAM7NEB9/kA2LMD7tOZmcz6V/e2IFzoDChp0R+4DHgHaRy+\nBu4FynN4j8kUf8M8GLgfef4lSJ5LoWEIA+cgSqfPMOAq4BNgNTANuCHH9/0GODHht6i35ZKDgH28\nbd8U59wGzEKe9Q3gUqAi4ZzdgVpgGdKR2zjh+K7Av4A5wFLgfeDngeO3A+dl+hBJ6AmsBT7MYZq5\npq33ORlRTH6a8LtBFKvm/GSrXRwH/Af5pucDLwG75CDdM5HnGw2My0F67WGyd+9UW+J32RH8BXgB\nqUcnFeD+qdiL5OVwINLm5ZPdge7Aa8AYYB2wfcI52wN1wA55zks6fEPLsrRlwjn9gQcAC6xAni3S\nRrrdkLr1a2AV8B4if5/EezYhhhSfXNe9JYkq1KXBOOBQ4M/AZsAvgZW0/ZG0h1wrWalIVKLSZSTw\nKdAFaYx2AZ4BHgEuyE3W8sYhSEUdbBz+CWwFXIxY0vLxDFEglId0E9kTebYvEUU3Ge8BxyPl93rg\ndODUwPGhwIvAXOBwpGy/BlR7x6uAJxGF/CfAzsA/iK/UHwO2BTbN8nl8jkEamE1ymGZHE0VkekLC\n78d7v3fUd5+K7YB/E1OiD0SUvpos0vTrmD2B54FvESWhIzkCUUgGAkd5vw0MbI8Ezs20TmwPA4FR\nwEPA90B9BmlU5jRHbbMAaMjzPc5HlEGQzvyNwB1AmfdbmXf8n0gHPpdk8t6jwLnEl6UvA8d7Ihb3\nHsAZwDZIPdnadx4BXkfq5nFIe3Ql0mYFOTJwz0GIsu6T67pXUfJCNaI879XGeYcj1s4VwNtIb9vn\nZGB2wvnjiblzTCS+9zkr4ZxTEAvqx7S0rPRAKpzvEIX310DXwPHJwATko56LdArSSTeRp4GpSX6/\nDGksncBvOyFW0OWA693TZxDSmH0PvIUoFkFL/6XefVZ4541NuN9kpLK5znueZxDLaWs8isjAZwjS\nwx/axnW/9PKyDFE4RgWOjaelO85EYu4ck4l/p68knPNbpEy8jij8rdETuBuR2RzgWmIyS7xPupbA\nK4A3A/tXIeUgyFKkjICUj8QynIz7iH/f2fAG4qrzAGLdC3Kyl58Dke9uunduWeCcb5Dv4W7E8no/\n0lj5jKf1dwhtl8dXaV3mryIdmJVAv8Dvn3rPlGgZPAhRHBYjiu52gWP9kbI8D/ne/4SMtPgM9dLb\nGylv84C/Jtw3kTtp3RXKT3PDwG97JeR7PPKcvwSmIPVfokVtD6Rj9AKwEHlflwIDEu63CWK9Xeyd\nd0/gWFt1XSoS8+vv74x0ROsQZfXnXrorvXweTbzRayKtf7sh4A+I29UapNwcHrhf0LqIl/blSKdx\nMfAE8VbJ8cTkOhWpL/f00tgFKScuUh/ipTULeA7YOpBOa3IfSst35bcFiS4fuwDvIt/CxwnP7qfT\nnrLXE2gEhgd+qwZmAr/z9i/x8lxF2+8/3ff3e0Tm95P6naViNnBSK8fHe2m3h5O9e7fWYWqmbbeY\nXNa9ipI3piAKXKrKeyjSk/8Hokjfj1Q6vtJzMq0r1D2QyvJvSKPZJ3DOUuApYDdkCLWJeMXuIaQi\n3hcZ8n8NacB9Jntp3AxsjjSMV6SRbpBq5IO+Nsmxrbxjv/T2eyMN1L1Ig7UNcFrg/HcRC8TOSMX1\nKfFDVecjFfdIREFaB2yU8DyLgWsQZeM2RLFujS+JrwSrkOHtC0ldie2JWGMnIDKahFTsPuNpqYzd\nhSgoIHKYA5yNvNNe3u8TvfxPBHZEKvN1iOU/FXcBXyA+mKciDYrf4PRGysCD3n3SUTA2RZTQywO/\nvYgonkHeQIYhQcp1PW1X6lci5SpbhiHlqAY4AHnmoLX/ZOT9vOQdPwHpwB0QOGc2UjbGIcPGTyDK\nns94Wn+H0HZ5TEehvgKxIJ3t/bYV0jnal3hFbzRSLs9CvtULEJeXvt7xId6xzZD38SDxnaKhXnqf\nIOVkb6Tc/KmV/F2MKBEjUhz302xLoV6GlKE9kXpkAPHlvxz5Xk/2ju/jnX9fIJ0IIpdngf285zwn\ncLytui4Vifn19/1Rm02RsnWal+5wRH4LvXz4TKT1b/cwRBncF5HbWOQdlCPlr9mTQX/v/JOARUgZ\n2x+x0L4YuN94WsrVz/trSOfrFOS7fAhxadsB6YB+EEinNbmHEWt+s5ev/sTqxKBCXY18b/ciBpMb\nvfsa7/hQ2l/2dkTazcSR+r0RBfMgxFjj1zltvf903t9SpGO8A1LmU72zVMxGOjHfITpBorvHR0j7\n/wDSubmW+PoiGY8j9djN3jW3IXVEkGbvnl8j5S7Z95qruldR8sphxHybkvV8z0YUvSAWsZ5B2wo1\nJG+YxyM9+ODw62fIUBJIJTcX6en77IB8eD6TESt0e9JNZAvkg07ln/sdcLX3/ymk9nndEKm8gq4y\nRyOdiVQ8ggyx+UxGKm2fPkhnINVQVxhRShIn15yGWFGWIhXajxKOX0u85a4L0nhu4u2PJ7kyFrRu\nzqal5X8i0oj6ymHIy0OqCSVhxLq6V+C3XyCNbzDNdCY6/g5pFJuBmxKOfYZ0MILcjIxM+NyENHTf\nIo33Zknu4Y96ZIuvhILIwBLfOJ6MPEfQEvc00in1+QZxU/HZFmnAu3n742n7HSaSWB7TVaiPIPbO\nrkYaY9/a6DOOmLXR51nkfSejBqmT/O94qJdesIN6AaI4pqI/Uqc1Id/VeOKtin6abSnUzcRbVyF5\n+Q+yGfHfwt6IEpSMdOq6VCTm198/to3r/oa8J5+JJP92/Xr+SkSRSsZQWsrxMUR2PiOROqa7tz+e\nlnL18x60Dn+A1BE+myP1e3eSkyh3P81Eggr1obRsw95CXMcgs7J3fJI0fW5DyqTfoc/k/Sd7fyuI\nb39ae2fJOAcxBu2FWMubiFdulyLtzc3I932391tr860+9q55xEv370g9FXzvlyCdsgORUarVxOox\nn1zVvSWL+lCXBk8jvcxjgBnIxz6DmGXiUGJD+j4vk5sJex8hViufTxDFAKQ3PQBRNlZ62yuI+8U2\n3jnRJHlrK932ErQcHo34TSbjKKRBWRbI712Itc1XCnZHKtHp3vHDiVekosRbcRYjFe02JKcvYnH5\nPuH3O7y8/AKp8J5ArOc+hxAvtzXe8daGA9PlFWI+dVHgc1LLfidkBCNoWX0ZqVwTh8vb4hbvPqcj\nFrG/pnFN0PfvTGADxOq0ITK68LeE8+d556RiFbF3/9tWzvsZseH+ZsTq87OEcxYR37lKLMOJZeVT\npDFNtP60RlvlMR2iyDD8RkiH7FjEQpjoX38UcBEx+axE3lNQeToLsUItRDo21bS0qCU+c2vf9QKk\nkR6NKBZHIpblI9J6shjTaPmNJVKNWMRfROqed5Hv3v92jyb1ZL106rr2klgvbokoYJ956Z9Ny3ed\n7Nv13XIeRBTQLxFFzSE1Fci7DeZhBiLDAwO/pZJr8LoviTfoTEF0C996mkru7anvD6Vl5/MlWrZx\n7Sl7g0ldZv6GfB9+/ZLO+0/n/b1JzOUG2vfOQEah30Hk/QsvveD8iDIvj2cgFvSTEeXYz8eXgfzf\nHLimwUtnMtIR+Yz4b/CvSMfpBWQOywJicwN82qp7Oz2qUJcOUURRPAexJAxCCrZ/LLFxDO430fJd\nJ/YuU7E6Yb85kNYGyIe1VcI2knh/52SVVrJ0U02ynIk848FJjm2NVEK+O0RrE/E2QCrZYF63RIYi\nlyMVywOIdeVspKJ8kJaySpzg1NRK3hchVp9E6xlIJfYo4ns3BhmCTLRkB2nrnaayCAWJ0lL2ydJK\nh/ZOeFwKfIVYly8HfkXM53ge8f7FePtuwm+LvesPRKxRpxPvrjIIUchSsSWxd39rinN2RYZtH0Xe\nUQPwG0TZC94rnTIcLCt+QxoJ7Lf2DtMtj+lQjzzP7UhjmmwUxyCNevD7GI34ioIoFVchCvUhXn5W\n0tLNJ/GZ0ylb0xHlZUtE+fddLYL+vj7Jnn9+kt8SOQlRNO5BRrv8oXw/vbbqjnTquvawIGH/NqSe\n+B3SYb2W+Gdt69udhrgq/RGpF6cjEVTaQ6IMUsl1TcI1wX3f2uyX81RyT8c9LFW+SLIP7St7luT1\nMsg3F/ybzvtP5/0ltoXZvrMXiJ9XMZf4Edco0oHxZX5gIO/jAtd8RPxE1bdofd7WJFrO52ir7u30\nqEJdmixDFAvfF+5ZxHcryL7Ehss/RhrLXoHjBxBv/VtFLKJCawSv+RypQHogvlfBrS6NtFKlm8ha\npJFN5pZwCNLIPOftP5biPJBe99aI8pyY3ybE0tYHUdT+iyjyu7SRt7ZoRtx1hrVxnkVkFnynQReX\nrojSHXynmweOlyFlIJN32hrvIkpTcOLlfoi1wm8cMpFPOaJg+wrT+8T7A4aQoeFaUjMTsf4HI30M\nQ5T2VATfeaqIJCcCD9Oy8fwOUapzRVvvMNfl8T7kPd6f4vjnyHBy4rcxzzv+M2RY+E7kvfQnXva5\nYiYxtw+LdEqDVv0DWlyRHicgHYr7ENlvknD8McT1Kpmi9hm5q+uSMRLpUF+IfPtfefvtfdcNSKfr\ncGTo/pcpzqsH/kd8HTMKMU680M57tpXHn9G63H0luLWoF8/SciRkP+JdwtrLDKRdTEcPauv9p/v+\nkpWtdN9ZMvZH6gafd5FRxeD9diA2ojAnkO9FgWu2Jd4tZCdkJDLd+0LbdW+nRxXq4mcAUnB/glRE\nuyEhfHohbgMgkSb8MGI7I1at7sQ+iC+Q3vWfifWcBxP/cb+FWEdHIBPNUhEKXPcqMkR0F2Jl3Qhp\n7G5JcX5rtHXOhUgj+ySiFIxGrIaXIz1t35L5JKKI3YfIYltiPqAPIMrt48iw90hkyNB3PfgI+SbO\nQ+R0G6LQBPOW7vMEmUb8xJDtkbi7h3h52N+710IvbyDv9CfIpMTdvd/nELOGvIa84/GIon0nLa04\nbyHvY0Niik+yvLf2TM3IyMiNxCYl/p74hqwteQxGfPC2RhrtU5F3djOxBud2JO70Y4hlZDJi2X/I\nO3460tjsicjyaCRKxX+It/qMIHOLIYiCfgyidE4JbF8iZas1n9z2kuod+uSiPAaPv4G4IKVytRmP\nKJT/9PKzNWIZ94d230HK095IWbgG6exmw7+Qb3tH5Jv+BTJy4cdk913Gzkca+fOJd0loD+8izzIG\nsZqem3D8NaS8P4N8k5sHzplM23VdNsxGrMEXIHL/K9LBTHzXiQR/Ox+pM0YideRYkrvb+TyDuFGd\njzzLvYgMVrRyTTLa+v7foXW5+7H4j0Xq+GSK9cveMX+y+Y3IN/FMO/MaZDqS96FpnDuZ1t9/pu+v\nPe/MD6+6JfLsf0fa7GAkmn8io2tXI4r0zUgd8lIrz3abl7fbke9wgnevR73jeyPf5aZIW3QXopf8\nJyGdbOvekkcV6uJnOeKm8HtkgoXf6B6BWKpBJj8dg3xc/0V6ij8ifgjn/5Bhq2eQj+dG4hvvhxGF\n7X1iw8FRWvawE387DAnfdJmXzxuI941OJ41UvwWZjlQk/kzvtxC/uqOJn/ixBLG6DEP86Z5BohOA\nyGN3RF5XIwrLlcSU8cXI0ORxiE/zSlrKqa18JuNe4v3RvkWUwKsRi821iPXzJ8SGSycjlefRiMUj\nQrwv60rv2D6I1XAKouQG83YbYvWcQmxyXCayPxcpE7ci8nqUeN/ltq6vR8rmi0gZPgR55psD53yL\nlFnj5bUJUax9618tMqR4D9IAn4ko35cE0qjy0ri3lby0xWGIzP6X5NhTSOMyyNtvrxwTSfUOfXJR\nHhOPLyE2jE3CscmIW9W2yLt6EVFev/WO34t8Tw8iHaI/0NJtIVleWsvfy4jy+jTyTR+E1HXBd/g7\nxPL1NFImLiGzb/IGpI56FlHaz0y4rgl5F5WIsvAK8ZNO26rrWiNZWQnS6OVpB6QOH4BMJm3rORNH\npH6P1Cn/QpTjG1KcCzJh7QbEL/5BpB48MeH8dN5nW/JvS+6NSCi9cxGlNNlkzTXItz0Skc+uiJJu\nW8lXqt98liPlOdEXONW1rb3/TN9fW+8sSD1SX7yOfAtlSL06K3DOh4ile2/v2dYgsko26dPH9c7f\nCDGeDEQ6Lf6zNSAGjVpi4VN3Id79KBd1r5IDnkQayE+QRmW09/uGiAV0DlKxmaRXK0rxE0aGF3Ox\n+puSmhOIn5SkKIrSGruznrsp5Aite4uEoA/eKcQiNDyE9ATLkB5ce0LLKEqxcRbphZZTMmcyqf3n\nFUVRkvE+bce3V1pnMlr3FhURZIjbj/36HTFluyfr+exRRVEURVEURWmNO5AZp9OQ+JTDifnt+cwh\n5gurKIqiKIqiKEVBsUxKPA2ZMHcn4kifTZgyRVEURVEURekw2hv+K990QWa21yCTuDZDZuL2QmKk\nJluF51vil1NVFEVRFEVRlFzzKfHRf4qGrsSU5AgSIsqPffgQEmO4ApmU+GCKNNSanTnjC52BEmd8\noTNQ4owvdAZKmPGFzkCJM77QGShxxhc6AyXO+EJnoIQZX+D7p9Q5y1Id6CC6IqHyuiKLOLxNLK7s\nxUjMw18iK2edVIgMdmJGIYtEKJmh8ssOlV/mqOyyQ+WXHSq/7FD5ZU5Ry67QCvUCJIB4MubQcqlR\nJXdMR9xrlMxQ+WWHyi9zVHbZofLLDpVfdqj8MqeoZVcskxKVwjC50BkocSYXOgMlzuRCZ6CEmVzo\nDJQ4kwudgRJncqEzUOJMLnQGSpjJhc5AZ0Z9qBVFURRFUZR8k1LnVAv1+ssob1MyQ+WXHSq/zFHZ\nZYfKLztUftmh8ssclV2eUQu1oiiKoiiKkm/UQq0oiqIoiqIo+UAVakVRFEVRFEXJAlWo11/UFyk7\nVH7ZofLLHJVddqj8skPllx0qv8xR2eUZ9aFWFEVRFEVR8o36UCuKoiiKoihKPlCFWlEURVEURVGy\nQBXq9Rf1RcoOlV92qPwyR2WXHSq/7FD5ZYfKL3NUdnlGfagVRVEURVGUfKM+1IqiKIqiKIqSD1Sh\nVhRFURRFUZQsUIV6/UV9kbJD5ZcdKr/MUdllh8ovO1R+2aHyyxyVXZ5RH2pFURRFURQl36gPtaIo\niqIoiqLkA1WoFUVRFEVRFCULVKFef1FfpOxQ+WWHyi9zVHbZofLLDpVfdqj8Mkdll2fUh1pRFEVR\nFEXJN+pDrSiKoiiKoij5QBVqRVEURVEURckCVajXX9QXKTtUftmh8ssclV12qPyyQ+WXHSq/zClq\n2YUKnYEcEKVzPIeiKMp6g3VMGDga+NC4dmah86MoipIGKXVOtVAriqIoheA3wEPA59Yxv7GOiRQ6\nQ4qiKJnSGSy7aqFWFEUpIaxjBgAzgO6Bn98FTjGunVaYXCmKorRJ0Vqoa4DngWnA+8B1xCrYDYFX\ngTnAK4ApRAY7MUXti1QCqPyyQ+WXOZ1BdhOQuv454GDABXYCPrGOuTDP1urOIL9CovLLDpVf5hS1\n7Apt2e0NbAG8DkSA24EvgWvxhgKBvwKXAqOB45KkoRZqRVGUEsE6ZivgI6AZ2MK4dpp1TC/g78Ap\n3mlqrVYUpRgpGZ3zaGCS9/93QA/v/56IpToZurCLoihKCWAdE7KOedk6Jmodc32S4wdZx1jveJ11\nzEXqW60oShFREgu7hIGzgCeA4Yj1YoV3bLn3d0gB8qUoiqLkhkOBfYCliNtHHMa1zwObA3cClcDf\ngLesY3bqyEwqiqK0l2JSqG8EFgO3ZplOon+N7iffH0VyfyTdT29f5Zfdvsov8/2g3IohP2ntW8dU\nfNXQ8A9v9wrj2iXJzjeuXWZc+3PgoOkNDfOBMcA71jHPXNy9x2E5yE9Jyq+I9lV+2e2r/DLfPyDh\nt0LnJ46y1g52INcgVulDvf1ZiE91T8Q63QvxWfk2xfXjA//PBKbnJZedC19GrRYQJSUqv+xQ+WVO\nqdZvZ5WHGIJMQv9XWycb176wc0XlQf/p2+8o4FzgkEOqqw85p3uP/yAK+dQM81Gq8isWVH7ZofLL\nnG/oePnt5W0AfVo7sRgcq69CZncfDNQFfn8I+AIZ8rsE2BT4aZLrS8ZBXFEUZX3EOqYfEiavJ3CQ\nce0L7by+P/Bb4AzEFaQZuA+40rh2Vo6zqyiKkoqUOmehFdHNkEgeXxFTpmcBRyFh8+4GRiBW55OQ\niYqJqEKtKIqSR7xVDfsCTca1izO4/mbgdGCSce2BWeRjMHA5cBoywtoI/Bv4o3GtzTRdRVGUNCla\nhToXqEKdGf5Quw4/ZYbKLztUfpmTd9lZxxwG/AIY5G0DEDe8PxnXXpbk/MOBw4EPve1T49q13rHN\ngU+RunpL49op1jETkMmJg4G/Gdfe0s78DQOuAH6GzAVaBuySphuIlr3sUPllh8ovc4pBdqpQK4qi\nKOlhHXM6cHPCz4uBfxnXXp7kfN8C7dMETAVORdz69gduMq49yzv/CeDH3rnrgM2Ma7/OIJ+bAP8E\n9kNcSnY0rl3W3nQURVHSRBVqRVGUzo6nYPYxrn0ry3SGAVsB87xtvnHtulbO3xLYA9ge2A6Z8xJG\nojedhViQRxrXLvLOHwN0AX4JHIusP3CQcW271xWwjukKvA1siay8e5hxbVN701EURUkDVagVRVE6\nK56P8wXAnxHldbSvvBYoP12AbYB7gWHA+ca1yRZy6Y/MoekJbGdc+3GG9xsGfADUkMItRVEUJQek\n1DmLKQ610rEkiwGspI/KLztUfpkTJzvrmAGIZfZqoKwuGj2rkMo0gHHtGsRaPQxxxUh0H/HPW4C4\nhYzJVJn20pkNHIO4mvzOOuboVk7XspcdKr/sUPllTlHLrjNYdtVCrSjKeol1zAHAPcCAl+rWrjhr\n6ZIVq6JRkMVQ3HakU2Zc25hlXkLA1sAh3rYDUjcfZlz7TDZptyMP5wHXAWuAnY1rP+uI+yqKst6Q\nUucsloVdFEVRlPazw9zGxgHnLluy6N36+r5AD+/3vyO+yW3iuWe8bx1zL3CtcW1Dujf3/Jf3RRTo\ngwEncHgdcAvwbLrp5YAbEFeTE4EnrWN2yCTMn6IoSntRhVpRFKV0+cuMxsa179bXXwWsQhTpi4Bd\nEH/iJWmkcSWxhbOubetk65hewE+Q9QL2RhZa8ZkLPIco0a8Y165O/1Gyx7g2ah3za+R5tgcesY75\nUbbWd0VRlLboDK4S6vKRGcUQz7GUUfllh8pP8Bcn+QHPdaLs7XV14V0qq0JI/RYGmr3Yzslk92Og\nFonIsTfwPrDKm6x4MfCMce2XiTe3jtnWOxdgJ+Pa9xNOMcD3XjoHAicAhxJToqPAe4gC/SwSf7rd\nkTqS5KvGuDadzkCq6w0SD7s/cJ1x7QWBw1r2skPllx0qv8wpBtlplA9FUZQiIAxsgbhJ7APsiYSn\n+2H5bOuY3YA3frRgPmUh2Lq8gq0rKti4rPzjgxYt2JEEBbw1rGPOR6zW04AdjGtXBY6VIcrwNrRU\nOrGOCW06b+5nDTByv8qq0IHV1RV7VVbRIxyOAq8CDwDPGtfOz0gSyfNbiYTaOxKJVLIgi7R2A14B\nyoGTjGvvyU0uFUVZj+nUOmfW1hBFUZQO4M/AIqTOCm4/DZ5kHbPL9IFOY6TleVHgLdpRmVvHdLWO\n+cI6Jmod84Bn/faPXej9PtvzhfZ/38Q65sqZgwZ/PbKsLO7+YWiqIvQG4k6Sc6xjQtYx//XydXcO\n0vuVl1addcy+wedXFEXJgJQ6Z2eoXDp1b0FRlE7DtcAF3UOhFftUVZXtVVlV+XlDw3Z3rl71aYrz\neyKLpOwA7Oj9HeD9TTt6hbfYy/tAN+AMf5lv65gRwK1IuL2vmqPR/wuHQscikTp83HfXrXv+98uX\nLZ/a2LADsBviVrIBeTJmWMdsBHyBuJXsZVz7Wpbp/Qv4lbc7A3jE2z5P1z3FOqYG8UsPI+4zashR\nlPUTdflQWlAMvkiljMovO9Yr+VnHbDOtoWF8RYgDh0XKykOhEIgbxolJ/JZbIwxs5P2fruwq/9Gr\n5swju3S5FqgHdjWu/cA6ZhBwNHDsU2vX7Pz4mjXc0LuGXuHwcuBx4D7gtYRVB2u8+7/Xjjy3G+uY\nK4DxyPLlWxvX1meRVgWy/PmJQP9ZjRLEZHhZ+TREsX7YuHZK4PwQMBLYNbBtEkjyBOPa+zPNTydg\nvfp284DKL3OKQXaqUCuKohQCb0Lfu4hluRl4CrgJiYKRb0vnhsBLQNVnAwZNqolEtkSU5R8BewHh\nKQ31HLZwIXVEObiq+h+31vS5uLVlxlthBDAbecassI6pAj5HlPcLjWvbjD6SRpoRxGf9GCRCSd/A\n4S+R5c83QizR/RIur0Ma8S2BhYh/t4bjU5T1D1WoFUVROpBqYK2/Yx3TAzgLuNe49rsOzEekHD5s\ngK1+2qXLnKt71RhiK+TWz29qemmPBd9vtzoaHYAsEHMymbly9EfcND4DTgGyfkbrmP0Rxfd32UT8\nSJF2GRIN5f+QCZC9E06Zj/irvwW8DXwENCCTHPcC7jKuPTWXeVIUpSRQhVpRFKWD+BEwETiMWEi6\nDsU6ZjCiKP7kk/r63Q9dtCAUASb1G9C4SXn5/4D/zGhoeGrvhfPvB8YiCuNuBDoB7WQX4EnEsrsc\nOBOJAlL0vsbWMeXAfsjzf4Uo0bOSjR5Yx4xCOg2VwL7Gta90ZF4VRSk4nVrnLPoKu0gZRcwfSWk/\nKr/s6Ky04H3iAAAgAElEQVTy2wtRSqPAX/J0j6Sys45xrGPOtY5504ts4W/rDq6qng1Ey8Ta6jcG\nv/TyuQgYkoN8DQCeJhYV5BFaWn5zgnXMJOuY26xj9vRcatpDVmXPOuZyT64zrGOqM02nhOms325H\nofLLnGKQnUb5UBRFyTM7A/8DugK3I5El8trht47pg7hFHIso835dWAe8ADyKxIqOINbXfohF+r9I\nfOarkQVZXspRlkLAz4HrEDkMIQfuH0GsY4Yivto+3yHuKjca136fy3uluH8F8DGyGuNVxrWX5/ue\nilJorGNCGt0GUJcPRVGUjKlAFlNpbbLddoh/bY+ty8s/eqpv/68iodCpxrV1uc6M5499OKJEH4Cs\ntggSweN54CHgueAiLh5HInXl4+R/ZG84cCqQF2XTOmYz4HjgOGLW9dnARsa1WU+KTOP+uyCuIY3A\nNsa1X+T7nopSKLzv7UVgonHtZYXOT4FRhVpRFKWd7A2cgyy1HUZ8g5cgludEi+4JwN1DI5Hayf0H\n7lwWCkWBA4xrc2X5xTpmDHARcDBQ5f3c5OXlQeBJ49rlubpfHtkd+CtwN/AwsCzThDx3j92B84C3\njWuvzkkO07v3zcDpSASXXTtCkVeUjsY6pjsyx2It8Cvj2ncKnKVC06l1Th2CyIxi8EUqZVR+2VEK\n8nsUqV+aiV+tcL9kJ/+0S5fTZw8avNbzr70oV5mwjhlsHXOv7xP9Wv8Bzf/t1/8965hfW8ckhncr\nBW4lJss6pDOSNalWQbSO6ZdwLCdlzzqmp3WM672XM7JNr4QohW+3mCkZ+Xkrlz7ilfHPgyuqFohi\nkJ36UCuKorSTXYB9EX/oRUAvZJKdC6wOnuj59dYi4eNuRyw50cDx7kDX9vj4ehPefgNcCnQB1iG+\nyTcZ19qMn6rwdEXcT04G9kHiOm+IKNc5xYs9/SViBX8ImOJtc3PhD2odcyTwGLASiU09N9s0FaVY\nsI45D6lzVgLbG9fqYjTq8qEoipKUg4FtgD8mO+hZNjdClECQ+maGce2ahPPuRlbiewk4yLi2IXD9\nr4BrvDQ+QSJhPAV8nCI0WwiZaHgNMf/gx4GLjGtnZfykxUcIGUreGpnIeGeub2AdszHi69wn4dBc\nYINslWrvXT2B+LQ/aVx7RDbpKUqxYB2zG/AqMkfjKOPaxwucpWJBFWpFUZQEfgHchkwsG4JYnn/A\nW1jkCmTp6SDbGdd+lHBuF2AC8Efj2mXebxsAd3xYv+6AG1auJArNF3bvEd6qosK/bC6iXD8NvGpc\nu846ZmvgemRFP5DVAs81rn01J09cfJwA/BTp0OTFN9Mbpv4psD0wGonO8bVx7U45St8gS6R3A440\nrn0iF+kqSiGxjhmBjL68aFx7caHzU0SoQq20wPdD0iGczFD5ZUeh5Xck8B9ksuF44FogLiqGdcyr\nSCi6ZcC3xOqZY41rp6ZK2LNa/uyz+vqbrl25otvL6+I9GXasqJh5Ta/e3YeXlQ8I/LwK+BCZYBcG\nFiMRMu4wrm1MuEWhZZcrRhKbnPhUB9531Oiy8qqpjQ2fJR6wjtkEmN7eCYbWMWcB/0Q6ZZuWyOTQ\nTOks5a9QlIz8PLezBr8O8iYBjwK+Na7NdBGobCgG2XVqnVMnJSqK0h72QfyRo8C4VCdZx+xmHfNb\nz/85LaxjBjzdt99LB1dV/zCJMST+1n8C/ob4CUeB+iGRyANv9B9wtXXMJ4FFWBqsY66zjkm2IEoZ\nMALYAHFh6EJsGfFSY0OkkxJF5DOysNkB65g7vHdwSAbXRqxj3vWuvzUf+VOUQmMd84pXxvcudF4K\nSKfWOTv1wymKklPKkXjFUeCf53frHrKOGZ2LhO/o3efMw6ur68KeIh2WuNB/RyYq+myIWGT9yCHL\ngd890qfvxtYxx3lLW/t0B/YHrkR8s1cRH23E39YgkybnANMQ14kHEDeKU5Hwf0OASC6eMwcMQBaZ\niRJbVfIdYvG0C4J1zIWesvC/DK/fwusQRa1jLsx1/hSl0FjH/NMr378rdF4KSFHrnNcAs5AGZtPA\n7xsiDvFzkAUTTIrri/rhFEUpOkYD1703YGBvb5nuNdYxA9q8KjlddquoPGz/yqqZZZ6CG4HmXqHw\nPcDgVq7bCphETCmeiywFfiziOvAxEmM6UXn+ztsWE1NG093qgRnIAg3XIYuijKRjhy97IRMzo97f\nYYD19gvaSFvH9LaOWe0pDJu2fUXSNE4IjDb8PNd5VJR8YR1Tkyr0ZOCc472y/XRH5asIKWqdc1dE\nWZ5NvEL9EHAZYrX4PWJxSUZRP1wRUwzxHEsZlV92FFR+1jFl1jGTvMbhe+uYPdK4rAuyvPhZXUOh\nh2rC4bnhgMIahujm5eXvlIuSmC77Ir7TyRTgBiQU37WIz7ev9AdlF0Gih/RDrNCjkQmNJyOW7XuR\nKBfzUtwjCixFlkz/E3AEUh/nQ8nu6uUliliofcv9/sSedz+gZx7uHSRl2bOOudkrE//KNHHrmLO9\nNJqsY36ScS6LF637sqPo5Gcd09WLM/1Iay5u1jHDvbK9oC3lO08Ug+xKIg71bCSE1RRv/ztgM2AF\nUsF+jlitE+nUDuKKouQe65jrgXMRV4kdjGu/SXHqjkg0kDFIhz/ObSICjCwrjw4pi8weFIn8euLq\n1Zm4C4QRy/RZiAvIm4ji+R7izpErugBDkTCAWyNRL3YABiY5dz5SF3+VsH2LWM7bSyXwDKI8zwF2\nQ+p4nxuBM5H6fDwSMaXD8SYlTkWs/8a4dkmG6YxDOjQNwCHGtS/mLpeKkj3WMbsj3+EwpB7YGvnG\ndzSuXZHimhBSN/QDRnSyMJ7pUhJRPoIK9XDE3WNI4PgcZAb8twnXqUKtKEoqwog72Q9Yx/jh8hqA\nfYxr30y4JoRYji9FJjACPyjPbFVezublFU1DyiLvbFFecWe/SOSxVA1QCRBCXFN2IKZgb48sYJOM\nemAm0vB+jsTH/ozWRwrLgEcQ6/d8pB6fkXBOFy9Ng0xSrPHu1eFYx9yLtDN/z0KhDiH+8+chnaL9\ndMlmpZiwjrkBOCfw0wpgF+PaL9u47hGkzjg7MXzoekJJ6JxBl4/htFSc5xCvYPskVuSJwwG6r/u6\nv37uVwCTgQuCx61jTvMmj52ScH4YsZK+j+cOUQVNp3frFn2qb7/ojEGDG1/o2/9165hTAlE4iul5\nc7UfQqxWBwJXAbcALxPzdU7cpiDRUpINx26MTML0XUsObeX+OwbS/GOS4yW1bx0Tfrpvv8e8IfKl\n1jFbFlP+dL9z71vHhB6o6XuNZ0BIdvzgR/v0u8M65gzrmAOtY/qmSK8nMQVy1OBIpLIj8l/E+ymN\nBwWdVd0KsxCDUE9kCLQX8kITlWyf8YH/Z1IC8R2LgMSCorQPlV92pCu/4Yjr1whv2wIYhLgOJBtu\nfBTxI65HKr6u3vX3IO4dGNfeYR3zqnHt114+KoDjgUsQBZAeoVD96d26V+xRWRXeqqJiGjKJ7/ED\nFy2oofD1iy+7fOUjihg4ZgNfJ9xna6TjsTHiq30MIu8rve1LRIF+GHHpuBxZvGU1oqC3ZvF9D5mA\nvg9wERJmMB/kW34AGNc214TDl3820AkDPwZe/EXXbifcvnpVoctPtnSI/DoxeZefdUwPYKIpixwB\nrLOOeda4dl7wHOPa55CRomA+apIk9zASpvNk4Ddzm5qOQOqBVPpYPhnq/e3IsreXt0HLFVeLltlI\nw+fzEFIZVyCTEh9McZ1OSlSUzkstya2iB6Q4/9OE8xqALxAXsucRhfse4FZESb4GGf2KAtFeofCK\nCT16Nc0Y6EStY5ZZx5xrHVOer4frBJQDY4GJiPEjKPup3t91iAtNOnRD3lkUeC7HeS0I1jFV1jEv\ne5bq2dYxrUV/UZSssI4ZbR0zzStvyzKJqx5gG+RbXIUok/d7+1flIKulSlHrnP9ALBn1yEz0z73f\nE8PmbZDi+qJ+OEVRsuJviGL1D2QS4SGIRbQixfm+T/AbtCOkXO9Q2L2mV++VswcN9qMz3GId0y9/\nj9UpqUL8pB9G/IajyLLuh7czneu866IZXFuUWMd0t46p9ZScL61jSsbSpZQO1jE/so5Z6ZWzL6xj\nsl0w6UHkO/y7t7+btz+f1HVwZ6ckonxkSkk4iCuK0ipdEbeAbNkWeBLpgM8DTjuoqjr6fN3aOmTi\nW3UIuvQNh/v2DUcG1YTDg4/v2nXnQ6qqh4ZDIRCf6/OMaz/NQV7WZ7ohLh7zkKgl7eU8RLFeCGwO\nLMhd1tLHm1w4yrj2qxyk1Qd4HZkr9JhxbWcMqacUEOuYjRHXqUnAz41rV2WR3AjEtaIZcb37DtG1\nPkO+yWORzvP6RqfWOdVCnRnJJhAp6aPyyw5ffuWI7+1cJBRTNhxHYOW9y3r03NE65j7rmHXWMZO9\nbYZ1zNrA4hv+Nts65qgCxVZtL+tD2QsjI5NR4Aly24ClJT/rmErrmA+98pPpwj+JaW7oLSQUtY7Z\nNRdpFoD1ofzlk7zKzzpmoxzVYzcj399dCb+fCUT7hcOfW8ccn4P7tIdiKHtqoVYUpejYBFl4ZHvk\nOz4Z8W9uLxHgL8CFAPtXVj19U++a97uEw5chbgjJWIlErZiLLOv9D+PatRncW8kfGyIugD2QsnF3\nR2fAOuZJxO1kvHHtlTlKcwIyL+gdYFfjWjUKKW1iHRNG3Nn8CdoVxrW35PGWvwJ+CxyEzIfw6RGG\n78ZWVfe4sXfN3OHz5qZaxbqzUhJxqDNFFWpFKS388HR/QxTeOcBJiLtFUjyLi0GW7N4K2BI40bi2\nK+LndwDid3vud4MG/zEUCgXjKL8F3I4ozxaYa1y7MsfPpOSHk4hNeNwMeYcdhnXM3oilfD4wxLh2\nXQ7S7I5Eo+oPHG1c+2i2aSqdF+uYQch8smHE+y0vNK7tn/yqnNEijj/AkEik+1sDBn2HRGLbwLjW\n5jkfxURKnTPcwRlRFEXZErgBUabv9vYnpzrZOuY/wGJE8X4GiVF8zMTVqw5F/AUPQHxt97OOiQSU\n6elIKLy9jGvvNq59ybh2mirTJcU9iDLRE+kUdbTxZDLiMzoACRGYNV75u8Lb/Yt1zPo6uUtJj8XA\nSESZng+8jYzs3exZrfNJC2Ua4NumppVIBCaAnfKcB6UD0eGyzCgGX6RSRuWXHdciVuo2sY550Tom\n+t2gwYs/GzDog4f69H31nG7d36wJh1ch3/9HFbChdcyfAj7Rl5SIP3QmrE9l7ybkHftRQ07NQZrt\nkp91zM+9MvVBrsqUdUyZdcxUL91zc5FmB7I+lb980G75eaHwuuUpPxlhHTPeK7/XdOBti6HsdWqd\ns1M/nKKsp/QA9gYu3Ka84hknEpmJWEsSQ949cGBVdQ/rmLu8yr3ROuakwmVbyTFHEP/elyO+1R2G\ndUy1dczb1jFnW8dEcpjuoV6ZXWwd0ytX6Sqlh3VMH+uYm6xjflXovKSLdcxYr/y+Vei8dDCdWufs\n1A+nKOsDXoSNBwaGwzWIb3UdLZXneuBDZOj/18C2nwwY1NU65jmvYl9tHXNggR5ByR9nEl8OXqQT\nzJuxjgl5kWei1jF/LXR+lI7HOiZiHXO616mKWsfMt45JNZG6IxhBmitoW8f0to650zrmZDrB99gO\nNMqHoijFh9d4/B04HeDXSxaverZurT+0+YG3fYQo0l8iq+751/YFngXGIH6GBxvX1qJ0Rn4PTAjs\n/wq4rUB5yRnWMdsD7yPlemPj2kIs59ypCUxoHu1tDjALmAJMMa5dXKB87QT8C5lkDTL59Rzj2i8L\nkR8kWtJUJJTpAciy5G2xMfAHpP49PX9ZKyo0yofSAt8PaXpBc1G6qPzSI4RMQHwImUzjM+qMbt2H\nXdq9x19DodBW9dEoE5YvY+Ka1SCT0C5ClOgWeJO49vbS3Rj4FvhRLhbfKBHWx7IXQjpexyHRMVYB\nWwDfZJBWUcnPOuZ+5LnuN649odD5SYOikp+P544zHFk4x1eeN0XCc7bmf7wAT7kObN8gMe3XIqNl\n9TkMbziqOhQKzRg0+EFkae9vgQuAJwocQvEo4FHk2UciUZPaYiRSDtYiIf2W5itzHsVQ9lShVhSl\nIFyATEBcDAxFFCGsY7ZrikbfjIRCVbMbGzh96RK+aGj4ArgYWeUrrmHxLNkHIJX+YYDvc/oZcKBx\nrdsRD6MUlDDQG7Hq/QSx6O1PikgE+cY6JpQLBcg6ZijwFRLFYXvj2qQdyc6MdcwoJN73QiQ0ouv9\nXZ5Mxp7P+Zbe5ofR3BxZDTUZixDr61QvXV/x3hRZpbUtoohivTbw92PgPuAF49r6dJ4z4Rk2QeKr\nTzCuXdPe63NMCImYtD1wFjIZOF3+B+wHnAv8I/dZKzpUoVYUpcPZE3gZGUo8ElnxDmCbCvjTf/r2\nGzunsYlLly+dtyoavQwJkdbkX2wd0xVZVOAo4GDirUxTgP8A1xnXLs//oyhFRD/E/acf4l99c0fe\n3LOEXgV0M649K0dpXo0sTPQqsO/6tNiLdcxAxLVrcJLDaxDl2lewuyPKc6qJqRYpG1ODm3HtohT3\nDiPuIJsGts0Qt5BqJLRnNa37FS8BHkGU67dL9N3tg9TVCxHDR3sUfN+yPRWRXSk+f3tQhVpRlA5l\nMOL73B/4K+Jndwzi+zoGoGsotHJNNPrnqLhu/FCBW8cMQ1Y+PBRpzHw+Ah4DHjeundYRD6EULX4j\nvgZx/ZjVUTe2jtkMKYsVwK+Ma7P25baO6Y0s9lIDHGJc+1y2aZYCnvvWy8BuwKfAF4gyO9jbUlmP\n67xzP0VGqT4FPjeuXZKnfJYhyrWvYPdAOvsnIAq+z2zgfsR9Z5p3bRioLPKVWF8C9gUuRzqL7aEc\ncVsZBOwFvJbTnBUfqlArLSgGX6RSRuWXmjDwOrArEvz/I+B4pBECCX32OHArscUB/Ibn10iUD78h\nfRdRnB43rp3dEZkvAbTsCQ8CxyKTp3ZCLIXpkLX8vNCMExE/032Na1/PNK1AmucB1yGWvi2Na9Px\nYS0EOSt/1jE3AWcg1uftjWu/Tzjeg3gFuw5RoGcY1zZRBFjHbInUb8cTb2X/ELG874osSnWkt9Jm\nMX6/+yCGjT1opx+0dcwxlyxbeuWDa1YPbRZ3kX/nI4MexSA7VagVRekwugB/An5BvE/jO0hkhkdI\nGFL0/EjvRCYb4p1zoXHtd/nOrFKy9EEmUHVDLIObIb6tHYJ1zLXIHIFFiDKYVYQOz1o7FfHvzYnl\nu5ixjjkVUb7qgT1KPUKP5wq0B/AzxMe/e8IpTcA0xJrub58Y187P8tblSMcuW1eLpMuMt4V1zJ+W\nNzdfWheNXrvd/HkXZpmHUkAVakVRssZBfKH7IRbkrsgksa6Ij6H/20bErNErKuC+Sf0GNI0qL78s\ncdlvzyr9K+Bq79pFwOnGtY92wPMopc/PgTu8/59AXEE6xIfTcwN4Dpks+5Rx7Y9zkOYxwMPIEtMb\nd9b5AdYxY5BRrArgVOPauwqcpZxiHXM6cCOipC4EViAdpWS6yjeIscHfPjWubWjjFhsikziPQJT4\n1cR8xqcE/v+GwLyUfGAdcxjwFPCqce0++bxXkaAKtaIoGRFGZnD/Gomuke5KcbWIS8cj1jGXA78F\nJgP7+JN2PKv0v5HhRpBJhmca1y7MVeaV9YLXgd29/68gPl51XvF8n/8OXJyLcuvFTH4bcWH5DDjU\nuHZOtukWE94kxA+RDvpNuZrYWSx4CuaTiF7yJ+AK49pG65guSCQSPyrJVkjYvMSQfmsRV5EflGzj\n2gXIhMkjgB8D26WZnTokgswUxCL+sbflrI61jhkAfI8o9b2K2FUpV6hCrbSgGHyRSpnOLr9+wCmI\n9Xi491szomAHqUMUmhuQCnUNUll/A2AdczjSuDQhyvTr1jGh/9Wt/f2wSNnFG5WX+1bpM4xr/5Pf\nR+o0dPay114OBJ4P7B+N+N2noqjlZx0zHHgByef3wGHGte8XNldxZCw/z63lFcSv+A3E/7wta2xJ\n4fl9Pwc8bVx7dZJTfpCf5yayGbBzYPOPM62hgUfXrOa5urUN3zU1lQfSWI2UkSe9e1UQi7s9OrAl\ni5wCMA+Z2/IJMSV7NhmO7ljHzAKGAVsb136aSRppUgzfrirUiqK0SQgZPvw1MnTuV+BzEN/nh4DL\ngLeQiWDTkaHppJWwdcxGiKWlJ3CRce011jF7AH8kZlFUq7SSLWGkPA5HOnJ7IVEHShbP8v0YMqdg\nLXCCce3jhc1V9ljH3ILULxbxO8/Wf7gosY4pz7SjYB3TB9hp4upVx41fvuynjZ5+UxMOs39VFQdW\nVbNrZeXq6lD4PWJW7NoUdWhPZGGbzYGtEYv4ViRf6GYFolh/5G0fI9btNi3OgcWJfm1ceysykvl7\n4F7g6/SfviRQhVpRlFbZBAn3tK23H0UsH/9CFlpplx+eN7z5DjK0+ThwPTCemHuHWqWVXHIh4of/\nIvCjAuclJ3jW3JsRP3EQt6m/lWicY6xjTgNuR5ZZ373IrO7Fxq5IWe4SgoePqu7y4p979aquDoXH\nALsAI5JcMxtZnMXfPkqxYMxdyIIys5HJ377ryYAk564l5iryEfAsMmoSh3XMdoiS/oFx7Wrgd0j4\nvSmI1X1Feo9dEqhCrShKUgYhqxOeScwiDXAO8M9ME/VWNvwHMBaxUOzlHVqOhAa7wbh2WabpK0oC\nfRCrZxUyKbZgVjFvsuKfgVuMa7OKj+35VF+ExHIHiYRzeiYr8xUS65idkPjEFcApxrUTC5ujomYb\nZIGfnojyexoJ0TesY/ojfva7IArr9rRcJbIJ+BxRrj8AZoxdOH+TLxoabkFc9bZFJi76DPLuva23\nbYMs8hJkORIa735adw/pgYQ8HY0o4T8mz5MjOxBVqJUWFIMvUinTGeRXiVgeRgd+W4ZU4jeSxWIZ\n1jG7Ihbp/byfViCK9PWeIt0Z5FcoVHbJ8S1v1yIW61TkVX7WMVcgZX86sItx7eIcpHkkshJfNTK5\n96h8LWKSBu2Sn3XMIGQS4iDgRuPas/OVsY7GOmYscCJwcjs6Oa3JbxNkTko/xOXnWNJzuShD6vEd\nkYWzdkTcPH6YRD6vqYn9F8xnWbSZy3r0XHB6t+4fInX8LKQD+kmS0I81iKvItshcBX+E8XHEdac1\nV72NEGW+N9Ih/G1bz5EGxVD3qUKtKOsjoUjZ8Ehll0O7DBi+ZvTP/vIlXmi7pdNrh85+5vrzGlYv\nHQJEqwcM+3rgmB9/VrPJbgvLqrrWI9aEZFsjsoDGfGTob/6eM95edsbrEzcoizYPQvxYDwP297Kw\nAnH3uN64tl0LBihKO9kBacCXIpOx1iJtw77ISnAdgjcp7XVkKP1tYL9crJJnHbM98AwwEFEoDjGu\nnZFtuvnEOqYbIottvL/7dZZJiNYx+yHW10rgF8a1d7RxSVsMAd5ElkL/LxIWb10W+euKyH1HYJsT\nFy8a+8q6ur57V1ZxT00fQqGkatPLSHSmp5J0EELIRPXrkRjbC5C1Bp4GGDNuUgiZFPl17YSxdd41\n+yCuKxFkpdzO4OKnCrWidDbGjJvUD2m0N0V86gY0rl05YMm0t2tWzP5kwMrvvuxTv3xBGUCvjXZg\n4+P+AEDD6uWsW76ASHkFkapuNDfWU9V7UHo3jUYZvGwem30/HbPUZYNl89hgqUs42sy5R02AEFQ0\n1lPR2NDcUFY+f2l1z5lNkbKnkErVBZbUThgb9fK/K2JxWJWwfVs7YWyhrG9KafMeolifAtztbT/z\n9id2VCasYwYjQ94GsTT+Xy5W9rOO2QBR4rZErIPbFeviR57V9GnEsjkTsdZ3isnH1jF7IlE2qpF5\nJmdk6ds+EIl6shGiVP+IhMWvcsAA4Op9KquuuKdP366I8WOE93cjYE/keUCU5buA241rE92nhiLf\n0p7e/kTgvDHjJp2EF+0pFG2eFA2FnwCer50w9jgk8s5PkLkzpY4q1IpS7Hg9/B7ASGBjpKIbgqyu\ntRqpyIJb/8Q0Vnz7OVPvvuiH/UhVN3oM3bKp54jtFvff5sDmUDg8MMmtX0J89lYjVugIMYvCmMST\nT3vrPg6cOjnut7qyCk762Q00RsoTT09kHRKyySIV/Mgk55xdO2HsjYk/jhk36Qok+oiveK9ElIoH\naieMzXrpZ6VTcDKiCHyAKNa/QhSeBsT9qMPKiXXM5ohy1BOJU50shFom6XZHFtLYG5n4u1ex+VR7\nvt+3IhbMRcDOxrUzC5ur7PGe6yrgfMRf/07EOt3uFQYD1CBuPFsgLnh7I77KHYp1TC/gBOSb2Txw\n6CViVmt/dCGMzLP5MyKHOZW9B/1667PvOhlpN3wao9Hmx977w0EnkIbrSomgCrXSgmLwRWo3ntIZ\nqp0wtkUFNmbcpKHIc1V6PzV729e1E8a2eM4x4yb54YSqE7a3kyloY8ZN8lemiqya+1UvINxt8Mar\ngMdrJ4x9OCGPFYjy9+NAulXeNhuY39zU0Hf5zA9GRJubevYcsX0kUlE1iJiF4AcaVi9j7cJvWbvY\nUrfIUrfY0rB2JZv//PompPyvRZTLVc2N9Ss/v+3MmnBZxctVvZ0XBu169Etf3nHOCmCDbhts+mR5\nl57bNtXXNVf06HPTsEPOuzUcKVsFLKqdMHZ18J7WMZGFXWuOCBE9M0por4ZIGY2RMlZWdvty0Ir5\nlTVrli/Em/DSGIq8ftYxVy1Y2L3vPojFvA/SSPhbM6JUON5f1i62AFT3MYmPC+ImMhMJgTY7sJ3m\nyT+R02onjP13kvd1PTJs+l3C9mLthLFFPVTeBiX57XYQ1UhnrQYZ6n4fGaI+FylX+xFTVvIuP+uY\nvYHfAMcZ1+Ys0oF1TF9E+TKIO9X5uUo7Ddosf9YxlyKLmtQh8eff6YiMdQTWMXcjftO3IiE/2zvy\nEJRfd+B/iOFiGhK2tKBWfK/TsPPcngPOtb2cH3/dd0jFzH5DmdV3SHNDpHz1urKKr6Oh8FzAXfb1\nh3UtnnkAACAASURBVA2znrzmYM91kLKqbv9+ecjIHadtuOUW7wzdrn7qwJERQqF7aieMPTXJrUK0\nP+51MdR9qlArmTFm3KQIUFY7Yey6hN/LkBnGuyGTDnoi76ERsQ69gShSTYG/23jXdEEUS1/JfBex\n5JQhkSYqvO0wb6sIbOWI+8BTQH3CdhRwfJLHeMLbyhPSOoCYr2+QtxHLT1kgT+XIxIytkpw/D2ms\nu3mbvxR3UtYs+Ab3zYdYNvMDmupWARAur2ToQWfTd/O9IBSKhsKReqCuce3KVR9efczgZPVOefc+\nNQ0rF7fllzwKiRJwovfM84D/Q95PC6xj+iGzuM8Benk/1wEPADcb137Yxv1aZcy4SV0Rxdog1vfE\nbUPio40kYwkyJLkYsXh/hCj33we2RchQe7LloI+rnTD2wSR5uxgZSl/kpb3I296snTB2XnueUyko\nVyOTEicirh4RJIb6T5BJt3siqxCWNF7kjDeQuuZo49rWFrPpMKxjjiMWBeJo49rHCpylnGIdsyEw\nIAdh/6qQ0KT7IHHTd0M6g+3CMwwNQpTzVFsVYnRZ4W3Lk/zfBxnV8Tcnnfs3NzXivvkw7hsPEG1u\nonsoxJ19+0d3Li8PLars2nBnj0F//8enk69EDD8/UDN6twfXzJ99QGWvAW8RCt2+/OsPn6c0IoGU\nrEK9IeIDNwKxWJ1IywKnCnUKPGW4ErHW9EM+mL7IcPkC71glouD2QkKb7Ykohb7SW4YoKKsQRbGL\nt7U5tr+eU48ooo31KxfXrZr7VXjN9zO71C2d133N/G9CaxfMBiBcUW2j0ejaaEPdyECHvQlpKJ/y\ntvuQg18l2VJVQDsClyAWXT/hRxFleUHiydYxNYjF5QhiM8NXI0s539VREQXGjJsURvwJhyBuLcMS\n/g4hvbLXhFh6FiBK8Upk2L8MuLp2wti3k9z7RZJ3sA6tnTD22STn34o0PEsRBXw2Mlv++doJY9vd\nMCo5w28v6pCO22KkzDwMbIB0pDvFBFnrmHMQv9WVyEIpBR218HyLX0Q67+cb115fyPxkinVMGNjJ\nuLZFPZEjwkgM6KOQ9nV3pMymzZhxk7YG/kL+4q4vR4xj74eize9f8PKtVRst+mbciqruGy/p0otv\n+pgpr47cdfKCHv0qAWflnC9GzHvnsWGr7LTy0OqlXN+7hsOqu7AuGuXsZct4pUvvJRU9+31RVtXt\nDULhj+qXzb9t1dxpffybRSq7NlX3GzKzut+Gv1348aQn8/RMuaBkFeqHEMvTX4FLkbAwxyWcU/IK\nteci0B2x9NYE/vYkZvXs5v3en1ivsxsxF4FmRDmuCmz5VHqbkQYLYhbiKFJR1HlbxNv3/1YhFa0f\nLcL/uwr5eBu9bZ2XHsi7rUN6t/4WJd7SnGjBrkCUp6D1OnHfv1eqaBZNgXMaErbgb6u9/Pt/V835\n3+0N8955bDjiE3c5EgrJpwmZnBfx/u+CdHRAGsZNER+6oIX7M8T3c0pgSzYsGEIq10uIxX2uB+4B\nrkEU8BZYx5zgnRP8jj4FLjOufS7ZNYXC6yQOQhTsIYjynWzrkyoNjyXI6nozENeSuYj/uu8uVEWs\nA3pF7YSxXybJyxuIVSmR/WonjH05yfmHIWV3KvBN7YSxncWnsBh5AYmBfhFS9kHqhWoK4J+aL7zh\n+YeRSV+fI0pgriezpZuXTZFVVHshcebPK0Q+ssU6ZjckBv9WwDZ5Wkp7HHAlYhneDXl3aTFm3KRh\nwB+IjcauREbpVga3xrpVa6fec8nBvUft9MbgPY7/MBQOd0PquB6IbpH4/1o8BdrbZiS6VnoTTU8H\nJiDvuRFZb2CCce3yMeMmhZob64cv//rDg9fOnT72AvvpHoevW9H1+MULeac+PmBJRa8BdOk/DCC6\ndsE3oXXLZL2Y/tsd/L9oc9P90eamj9fMnz1lzfczGwH6bX3AgQ1rVny2bPq7c9OVVZ4oWYX6O2Sd\n+xXIS/8csVoHyYtC7VnKqolZZP3Nd1PwlbfywP89kYlWXYkf/l+HFPJgOl2R0E79vOvz/S58f+JG\nYNma+bMWRqPR+q4DRyxHPqRliIJXhlj15iOK3zxE+ViJKI1rvL8NfrSG9ZSgL9e+iAK7GaIQb4Qo\nzH7ZXIY08M8gqw7+i/iJG2uQhuiPiOLcC5kZ/2Pvb/ck919EvIK9FrE++y4pK4FbtimvuOWZfv0H\nIZX2MOPaM/wEvMrxcmSJ2DCifD8HXJmnRiRIXn3hxoybVIF0PocgEx8Tt65tJNGAlP25yKjYbMSC\nNBOxQlvE+tkP6egOQJT84cDvayeMbVHpjxk36XNik33qEWV+KnBR7YSx37Tj8YrBj7DYOQT53mYh\n7zuoGBRUftYxlUgdcLVx7ZQcpNcDUYBGISO6p+R5NcUW8vNiTb+DfG9PIK4epTB8H4d1zNaIy181\n8u3/3Lj2vzm+zRnATUj7cAjwfDoXeVGdLkcU2nKkDrkJ+FPthLHJomechLg9fYiMpOWsTHiugX/4\n1Bn9y2VdeoZ2+/q9BZFo86XAxOAETeuY0IrKbpseudtJPVbO+fIn65bO26N+xaJhDWuWd69fsajC\nd3kkFKHroBGUd+nFsq8/gKiXRChMeddejWVdeq5pqlvZo37FIiLV3aNVvQetjVR1Wx4pr1xc3rP/\nGY0rF7+/ZOqbdS1zmhdKUqEejkQeGBL4bQ4yNBIMPh4dM27SU8T7uq5Bhvl8l4Yq729vpAGMJGwN\niNJbRsyKUZWn50qH+chHvYz4cGIjEAVrLfKMa7xz3gYeJGYZ9rcwQCAmpJIevRE/zG2JlaFKROG9\nNcn5tyGz2X2CH9zvEAtZMPbqKESh833YlpJ6BnQlMlHFD4/nb8mUbID55XDDZwOdft3DYd8XrjJw\nvLdx7TLrmKGIn+MuXn6vB35bbNEC8oE3IjSQmHK9IdK5DW41bSSzjngleyaiHE8Fvk/W2RwzbtLf\nkEUSNkFcD3wG1U4Y22I53zHjJvnhxmau553XTIggHZ8hwEFIh7YosI7xrZOLgP2Naz/JQZpbALVI\n23WacW2LSbr5wos1/RpSX9YikxALYiXPBk9JfB8pM/cj0Tuyjh+ewOZIx6MbstDJX1s/HcaMm9QN\niSpyEVLvR4F7gXG1E8YmLsQSpBZx/TsViXyTU8aMmxSqbKj7cl151ejTX5/IftPfBJko+wnxI7r1\ntBzhbW4mNHZxeVXV95HyHgtDkW7LotEuK5ubq29oDofW1q2uqlu9tLxhtbeYbihM9/5DWL1kHs0N\nCepMpByaGgiXV1Fe1bW5vKK6sbyiuqFu7Yqu0WiUSKScSFk5kUhZNBIuo1eX7t9EQuGmSIimcCjc\nTJRoKER0XVPTqBBEQ4QIhYiGQyFCEC0LhacFY3Z/PPWtrenkCvUPO2sX27jIAVnur1u72FYmOT4f\nGWJpWOVOr+zmjFqCFJSBaxfO2aG634brEIV23Zr5s6NdBgz7FPErW7Pws5d79dty368RZbjL0um1\nu/QeNeZzZBhytX3tvr5mzxM+CDSwo4i3pOh+/vcvJF5B9rkUWSEq8fyRyPKvM4Ffev/XIyG8Pgye\nPzgS2ax2wKCu3jm7TG9o2H1UeXk9Mmx7t3FtY6r8WcdUnr90yQ5X9+q98qP6+k3cpqYt/7xiWfXc\n5uZuyASSp5HO1RfWMTOBEbMaGxheVj4FmfT5ziEL53/ybL8BmyAdgx7TGxoWjCov/6lx7St5lGfJ\n7Y8ZN6kacL5//+kxA3c4DLxYrau//3rzrgNHDEYU8lT1x3JEsZ6ybOYHC3tttP0bxNw8NvLS7wZs\nvHjKG3v22XT36wIKs3//CmCFV/8sQeIrv7vwkxfn9tv6gH8nnl9oeRXp/iWIj+mzwKEpzh+NWPxu\nRqIs5D1/1jHfIRNmD5zZ0LBio/Ly/QIT3DJO3zrmxFmNDXcPLytfh7h+fJLv5+kRCm0yZdDga4CD\nga9PW7Lo+El1dbX5ul8+961jnp/V2HDg8LLy94A9jWvrcny/GkTZ3ABxZz0OaTtSXj/y6N9fXDN6\n1wuQETBWz5vxWtdBI8+pnTD2s2TnB/a3RzoHy5G2YU0b57d7f8y4SQesXWz/W10zeMEdD/zm4t5r\nV1w1q7Fh8PCymKep1/5kvP9pFNZ1rWFGpJzjV85nWmMzi3v154NIFS4hvlu3hne+n03D2lUxqzZA\n+P/ZO+8wSaqqD789YWc2R9alKdgls2QW2CIv2QIEJCnhQxBFJAiSg1piSVRQAUFBBIkKkmOBBHFB\nKERyhl1S0Ute2Dg7qb8/ftVOb2/PTMeZ6pn7Pk8/3dVddev27eqqU+ee8zsN0FnViLq8tnO3SgQx\nYDbyMoxGB8UY9CWWuSN778ErPoJ0mnS6s2ns8guHjreeQsb4koWpt8YNHW+9gzxKVvvi+d9AnoEl\nQGvrvM+GDR1vBUgHsu39B69Yac0DvJeQQTxx3uzndhw63ppNFHv76XMPJFfa8XvPB57zcbT7og/I\n2Xf8Onv5vdzPP3zs+mW8VYY+5XQUO/sbdPy9iY6XL8if+XwvMnjuAjZDswbfBB7brXnoZpePG78n\n8gRvPqu9bTpZHuOGrr/lKcApDy438aW1G4f8EfhbdiJgmLQeBHb60ciRNCQSTG9SEwvSnYec9tWX\n12b1JTMdezqw+OdfffXxo0ta/hO1MdIbPeYauqTn7jhn3lfnPrSk5emiR2iAE3jOYuThrH/v/svy\nGdwjgFUWpt7aeuh4a3j0/lTSneugY2YzYLOmsZNAhh3AV+sefsnzw5df3UdJp8+8fcvZ89/OPxU7\nFniAzs4t0LHoAM6IFdb6PPCcZSqy2a7fDOz4/sNXta+0w2HvBp4z4GcaCuAqFGu6G0pofTfPOt8D\nfoQUQLbti05ZqXBxmLT2Am6qS7An8HCYtHaxUuETZbZ77R0TltttFYWT3RImrU2sVHVzY/88bsKP\n0fh+Duzqt9T0ZOjPW9JpC9g7MqYrSQNdSbHPoeOux1kn2/X3Xfx5mPFgPw2c+vKffpSisFClo6Pn\n26h8gZgMpwOQSPx2vVmvXhMmrVtfbmv73ioNjfOJwmDfa+9IrtLQ+GVm+eOOzomrNDAf2XKbt6TT\nQ5FtlQASS7SckaxNtLYtGT190Rfzpmt5S9IdTVsu+Kx1y8igfb2tbchaY0b9s3XMmMQ7TcOGP9rB\npLrhoxNz0nXDp6Zb13y3tbV+aPPwxCLqEosTicQnrYvrPm9o/riDdF17Ol2/sGVRQ3NTc0c6nUgM\nq68bv7B1caJ5SDMdaf04i1oXJVpaW+YtWDR/CEB7R1ti7rzPuo1eiLOHGnQX9zLwK3RRWhs4IGed\nmk9K7CdMHGZ5ZI/famhKeTU0i7JL5CH+E1KmyeUNlNndiUI+JrFsGEcrMtCvQXHXf0fxdp/lPM7r\nTcouTFr1SL3i9yhsaDGaQryiyrGWPTEgj78onGQi8nyuHT1nXueWo2xB07Iz0UXkycBzFuRpbyWk\nU2sDQwLPuST6+M2s9dalK7FpEbrJuwUpjizV5iDjOlSs4nw0xQ5LH3vD0Fhti3JGtiCP06YahEmr\nMerft4HbrFS4TwXabEazVBsBdyDjsNL/8TWifa2LPO2dqKT4oxXeT58TJq1Elc6Jv0Hn3E/QDc8c\nejj3RVJ4/0GhIWcA5xUR9jUE3TwujzzgFS+oY7v+Zih0ZR6wUuA5fZXoG4frRk3GUMOysnmHoETF\nbIxBbSiHOpZOWCqWzYC7125onLBdc/MHR40Y+ejouro10JRbA0r0/Bc6+fwbeMpKhZ+HSetqFBKS\nTRq4EoU57UTXcf0xujjeDTxaaHximLTWQv+Zg1FcMEi948BKJEMZisN2/RVRyFrmsU7OKh0oFvVv\nwK3FlF+3XX99FPe/OkqOzPBk4DlblNPvGidz4f8MeQfzeR6Ho8SwbVCy4JF91bnoZvcE4NJKxR2H\nSWtVFGpW0QqNOfvINvhOslLhhZXexwAikxzYhjSnH+9p5Wj262l0I34TcEAJORTD0E1iQQmPxWK7\n/h+AHwLnBp5zRjX2EWNq1qAuhPTN4yc03NeyuP7Tjs72e1sWl2McGQYPCXRC+DbSpV1qijxMWj8C\nNqaruEvmcbOVCm+MVtsLFTxpvmzsuPf2GDosO94f9Mf7iZUKz83deZi0dkblxeej6dK3gdlWKlwS\nfW4hz9qh0XoZWlA4033AfVYqnJ3T7tjoOx3K0mXDZwN/Bn5ThelMQwnYrj8B2BIZ19ugpK6MBng7\n8AAyru8MPGd+Ee1OBvZGYQx3B55zXp51EoMk0TGBjMuN0I3l9d2stxHKi/kI3XzW9HUkTFp7opvw\nNLB7JeUvo9LnT6Pk2puB/ftxpivu2OgmuQmV9L6ix5U1I/VXdA5/DZgexxmmSL50H+CfgecsU9dg\ngDOwDerMi/UbG7lvua/dYqXC/XLWWWdkInH9lIaGiSMTde81JxKPf9TRcder7W1Pk2NIGQYFyyGd\n0W9Hy/uh6fH/ESYtD8nJ5XKWlQovAs5GiYsJ4Mq3l1/h/uZE4gfIK9CEPMG7WamwLM3MSGd2ExTu\nsWv0Ops3kHH9LEq82pOuGO35KFTkL8Dj5qIXb2zXH4Nu0vZHUowZ47oFJdb9DYVwFKw80J3hbLv+\nr5ECwN3AXYHnDKjQmxy+j8KvXkE3yUvyrJNARuLLwI8ZAFrVWWoiC4AtrFRYsNZxD20m0LlybzSe\nm1mpMHYGX29EoTE7VFlnP4l0nZcH/oDk8nrEdv3jkOLSAmDTwHNer2L/DKUx8A3qOmDTIUO4dcLE\nm6xUuH/OOjNQ0uFSrN3QuPjBiV/7hZUKe5WuGYDEIRapr9kFSditHS0vQEbx33JXjDSaD0PxqB1A\ne0s63Xn03M83eKCl5RiUNNYOuMePGHneiaNGH4+m3OuQEXtoNaSjwqT1NZSgtisq4jI6Z5U08BAK\nlbo9xvJVg/H4Kxjb9SciD/P+yIOdYf78D169f8GHb9y6/GZ7/b1UL7Pt+q+iKeUMb6KY/d/l09Cu\ncYai+PJVUSz1VdH7ucfe/0qV9jdh0hoKbGWlwn+U0UYCzaDtj+LCp1upsCxvYpi0Tp3d3nZeZ5oF\nqzU2TrNS4VvltNcfROPyZ1SW/lQrFf6qCrtpRjaHjUL+dqLLeZf33Ge7/pbRNg3AfoHnxKKUfMyI\nw3VjwBvU9UDnfRMmNqw/ZEhDnintMTOamr6+buOQLed2dm4yp6NjzXfb28dt1tTEr8eMvcBKhSfn\nrL93M1wytq6+pT7BJ2n4YGFn56x56fSjnZK/y2U0ksRpzXrMp8anDQtgEvKk7YBUKvKJyzeTP26x\nWjQhiaH383y2DZp+a0EnuWPppnpgN9teAqwfLf8DOC5MWrOQ9+Gw6P0zUdWoql+YI6N/c2Rcb4S+\n2/VWKszNMzDUMFHs9X7IMNo066PXkXF4bZbqUKFtjkY3ZLsjlYax0UcrBZ4zEI+fzVHsagL9l3uM\nY+1Pov/1vcCOwHesVHhDGW0NRSFiNsrh2KHUkK8wae2IwpDqgD2tVHhXqf3qT6JwvotRcvaWVip8\nrgq7ORP4OboObYqSEbvFdv1JaJZxeeA3geecWIU+GSrDgDeoi/4eYdJqWpzu3HBoou6LPHfZPwZ+\nm7vNQcOGt58/ZuyJViq8OOejo5GCQjYdCbg0Dcfl2f3K0eMpqidpUy22Q2EFO7J0UtW3kGc2l5nI\nI/wWuqvMPD+A5OXKYQgKgZgWPTZCwvlvsmzCF8i43xCduFoBwqS1GdBkpcLHutnHCsjznFGXeQ9l\na98RJq0JKMt9a3RyPsRKhfnGwGCoCFH2/3dRotPXorc7UEjIVcD9gee0dbN5d202IHWL6YHnXJDn\n8wQwIfCcfOXua4mzkWLCbHQeKDguva/JCtfopHyjehIKZ1kRqYocUuwNf5i0JqNY9PHAL61U6Jba\nn/4kTFozgIeRE25/KxXeVIXdLIeOsRFodvxfPa0c/f8eitadCexQ7H844gh0HbqZvnViDTYGtkG9\nQ1Pz+Z92dtzzYlvFYqKHXTxm7KEJEjvPS3euOK+zc9L8dHr8JkOGNO3cPPR7Viq8Kmf9746rq7uo\nEUa2AUvSaRam0xw5YkTHT0aNOdpKhbnV9U5EoQeZBLN7kEeiT+SayuTPdHljF6GTxcPIsHwnz/pv\no6nWXNZGSRe57Il0dzNJgPXR8+9Z9rddjmXv/NPIczeNHk4qUXb9aeii9SmwvpUKsw2GJmQ4/xSp\nALSgIhG/AhaHSevrSBFgCipRu2dv8nUGQ6WwXb8RhTAdhuLrM/HWH6GY+d9XKnTDdv3N0f/8TpRU\n9VDgObU4+zYEyRRuiGKqf9C/3emZChvVGyKv/HDgjHyJ0j1s2xxtuzGSB929RsuKD0Mx8isD51up\n8LReNimV3yKn3H1o9qdHogqqJ6P/7rTAc+aUsM9mIEQ3PFsgZZuKYrv+yehYvKKYJOkByMA2qDMv\n6oEJdXULxtbVhUMSiVcWptNPvtve/kSHDKxyvaGESWsU0GGlwoV5PjscxUmNAcYtSadX7EinJw6r\nq/tWHq/lgZPr6y96v6NjQo6b4CSgr+SHuotFGoemgj8nf3jLjmjK9CHkYe/tBia3zPPq0b4PIH9y\n0OfkL/s8Gmle5uIj/djnkOf5BRQbnZdoOnUj5HWeEb19IVLjWIJOtociL2CmPPRt6Cbo3Ugu6sLZ\n7W27AqzS0PgfZEyXchIczMQhFq5WWWrsounig5FxvVb0WRtStLgg8JyyJBJt1/8RMhIyRvs7yCC9\nOl/J9JizLjpPNKIbkWompZVNjlG9s5UKHy6jrT2B29E5eR8rFd5WwDbZ8cbvAJtYqXBC9HFN/XfD\npLU2MnK/Qt+jFC9wb6yEZmGHoOtMvrLy//v/2q6/N3JGdQDbBZ4zs8T9Hgxci66DG1PhXADb9VdH\nMxQjATvwnP4qBhaH68bANqi3a2qe+05725j3OjoS3R1FIxKJhSMTde8Dz3zS2fFQhw6O19GBXBWi\nO+LOfDFrYdJ68NOOjp0eXdLCQy0tPLqkpcMdNfrSg4eP+EV2hbyIDZAWcTUuXvXAeigOdzek21qH\nDNVdetowSpDLFJzYEDjWSoWzyuhLE7pQt6GEv47ouR15issOjwmTVsYbABrT71ip8HGUtX4YCmnJ\n8BoK2flHmLTGoZi4o5DHfD6qwnZxRurOYOhPotCMzZF3bB/0PwbNgP0KeLyMJMYkMqoORzrpAD8N\nPOfssjrdP5yEbqg/RgZ2bu7HQchAOY7CcyyqRmRU26hIS1nnmjBpnYISMxcBW1up8Nlu1mtA4XTf\nQjN1i4HNrVT4Qjn772/CpDUcmFTmdaonrkRVEG9COQ/dYrv+mkjLeyRwQuA5y4SZFsFT6Bg5POpD\nxbBdfwZyKo1DKkO9et0HOAPboAYSYdIaMqu9bf3nWlt3TnV0bP5ZZ+fan3R0WO+0tw+Z1dFOS3rZ\n60gjtI2sq5vVnk4/MS+dfhx5Ll5FBlxViZJFdkHSbd9oTaeH1QP1icTKVip8N2f115Dn6RN0x/tC\n9LiLwuMAG8j/vdaP2srQhqZ3bwcu7abvVyGB+lzd5QlWKvw8z/qPIm/0LCS19CLwspUKy541yNlP\nAsU8rwe058uQj0r+/jqdTv/ju198ftdDS1q+iU58o6JVWpAs1NXAP6PQkCOQl2gc8hRdCbhWKiwq\nEcxg6Cts118NFQz5LpoOBl10f4Vk8kpyJET6szuhCqA/DDwn36xR3KkHHkEzbbeihM/sC8S1yKA+\nHYV59Tth0mqwUmHZ16XoHHk1isH/ECl/pMKkVYccN9tHj61ZunLrwVYq7E7D2yDWRPZDGinodKuA\nYrt+E4prX59Iy7sMXfiNkTzfV0iqr2J5WbbrH4rCvBrRjfkBcdTF7mMGvkGd74OMgdWaTm/0Ulvr\nNrPb2zd7r719vTfa20e/1NZK2LHsNaUe2ocmEq8vSqdnduog/S/6k1RjeijTz+Fo+nGalQpPzfm4\nrg4erYfpbV0XxgwTUfxvLs8i43khCslYBVVOWiHPus3IYH8ITYc9HCat+dE2n1ipcBmDPUxa/0In\n3AVojAI0TrfkJruESWsI8m7U5bYDjM01qqPfrAEZxL2eYMKktQbyXk9FNx0joo9mWqlwm9z1G2Gt\nNnnvDmJp2bAAJXXdRKRBGyYtB5WMzaz3KHB8rXtpDIOHSH7vmOiRUfJ4EzgXuK5Uw7qH/dUBUwLP\nmd3ryv3LFCSlN4JlC758EzkU/oO0ugcUYdJqQipFWyPnxiyknz82Z9U30TnvNisV5gv/MyzNzejm\nrNf4fNv1L0Q3vG+juOlyYpIzbf0OzSZUhOjm+TFUfOp3wEmVPl/UKIPToO6OMGktD9jvtbfPeKmt\nddt32tvXfqWtbchLba28l8fIrpNR+kKnTrD/RYkNIZoyrPoBFiatrdLp9MwPOjp4qa110ZOtS+a8\n0tY2//YJE4+zUmFuBnE9+T3RnShhIWPA/i8WKTLoN0XhHptHj+XoJs4uTFrTP+7oaJvxyUdDFqTT\n09FU08boLnYxukNeDCyug0UbDxkyZK2Gxsa1GhuHrd7QMG6l+oYJzYnEkA0/nnMZGr9MaEfHtMYh\ndXctN/EcgHQ63Z5IJNrRzcwHVipcRrkjTFpT0Q1Phs/QDcJjVir8GTo2NkQhHfuwtBH9CfJGXZ1p\nIwrT2RMVg9g+Wm8WiqG+K8vIj0MsVy1jxq90ih67qJzxYejCm5lZehE4OfCcihlLtut/E3l9bwTO\njmlhisz4bYXig79Cs1oZucBhyFExDI1VPgnOfidMWk2lhoCESWs55ETILlP/HvLcPwI82kNRKvPf\nXZZpyDZYAqyG7IO8rPyN4w4ZOXm9vwwdb3UAW1QgHrkehWy+TH5hgJKxXX85YNfAc66pZLtlEIdj\nzxjUPRFN608F7I86Oma83Na63Ztt7dZLba282I2RHdEBzEF/ng+znj9Ef6z6Hh6gJLu5yMjN9BTu\n5AAAIABJREFUfl7qJBkmrWnIoNsCeVYy3Gmlwm/m9Knu4GHDnVUbGjZ5s739q7fb2+ZvMqRp3mHD\nR7ROqq+fZaXCV/J8/ytQ7FU2nyHR+6vQ+FrI4M48NqarIl/FWKG+nicmTqIhsfRP+llHx6INP57z\nR+BdZAh/Bcyz6usX3jh+ws4JEs9PrKt7bs2PUp8hb/hmyIjem6UvGnORWsGtSLqvLZru3ApNY3+L\nrqnOeShO+hITJ20YCEQSXQei43ql6O0HkWH9YgXaPxU4C80ypdGMzy9ialgnUHnuPZBS0c501Q64\nBd2AH4c0i2NF5GV+HBUCOaOUBLswaa2GZi5eAh6xUmFFjbE4ECatg9HMacHVRUvkflRw6zfoWp0X\n2/XHoxvZJPCzwHPOKmIfCRTz/zIxKUA0SDEGdbFE2p1bAVt/3tGx7avtbeu+1NpW91JbK+90tDOn\no4PPO6umHNWCDL+P0J/nxczj7eVXqGtOJDZHxvWz+aSUwqR1DCpCkst5Vir8CVLMGI9igsddOGbs\nnts1NX/jjfa2OcGSJV8+tKRlySttbU1prZN5jMjT3usoLvMpFA+2AHl1hkaP7NeZR0YKr76b12OA\nyXUwpR4m1ycSTY3IQv4qTxx8Dp3ICE6z9PTlx2gK91Y0hdUWjdOqyIg+mKWN7qdRpcG/5UkQNRhq\nHtv1h6LCRmeg/IE0ktv7Wblye7brT0bxx4ehWatO4OuB5zxUTrtVYiI6xy7H0sbzgcANKG509/7p\nWvdERVZ8dM78N9JUHogFeUomype5DXmOp1upsFoX7EzBsEyoZL4CZ5mk4UzZ9ieAGQWEUIxD3ued\no8fXUGhjvyfLDmKMQV0uYdIajYzYrVHRkJWXpNMrf9zRMeKjjg7mdHYwp6ODj6JHG2nqSdAA1CWg\nngT1Wa8BvujsXPxpR0fL552dHV+lOxMLOzsbW2F4usuDnY9PyTKwkQE5Ehm8I4ERhw4fvu7uzcPW\nH19XN2psXd3weel0+uOOjvRdixd1Xrto4TDyxzP3xpd0Gc8ZA3pu1viMQEbxEmBJhU5edegEMhl5\n5jOP8eimIPsxChnwGd5DJ9PbkCZnR9TPKegC+W0UG/a/r4CKHlxrpcI4etMMhopju/4E4Gd0qdcs\nRl6288vVmrVdfyVksG8PrBt4TiVqBFSDTMx0C6r6+m90PtkEJWhXPUm9FMKktRXwN5Qb8zlKHLy/\nf3sVD6Lr9avIE3yMlQrzJthXgAQqxrIlqo74i+5WtF3/MBRiNB/YIPCc3mYEJqPr7MSs91JI1rXk\nkvTd9M1B6iQHBJ4Ty+M9RhiDuhpECXTj6ap8mHlMQQZuE9KjbMp6ZJbzGrbpdJrF6TRfpjtbZrW3\nfzBzyZI5TyxpaX2zvX3M4nR6DbrUKMrhK/SnnotCVr5AJ+SeHvOAdJi0GlE28/oo7jDzvOLSu6Ad\nxZ4v6ebR0s3ruSxdUfHDIozzRjQ+Q1HYTToK59kMJX3uztIVFBchr/U1wD+LLFYQh1iuWsaMX+lU\nfOwiVZBzgX2jtz5FlQUvDzynrKprtus3BZ4Tp5CpfON3NTJU0kjV4Ax0Xow1USz0dah2AKiU9r+r\nvNvY/3fDpPV7VME4QGNSrVyn3dAsxmfIO533JjT6fz0PDP/itcdPeevvZ91J7+OXSTZ8AeX6PECX\nikjFsF1/A+QkawYOCzzn6kq2X2HicOwZgzpuRIaehf6Eq6BqgqtkPcbnbpNOpz99s739uQdaFs/x\nWxa3vtrWNra9SxN5QeY5AQsaYVFzIrFoRKKuZcWG+sXrNA5ptYcMad2qqbljdF1dM6qYlXnkhmLk\ne6yAjOepyHDNJWMQZ24cKsVili1b/ilZiY85j0VR/3ZGRvSuqPJihvloqvQu4A4rFQ52CSCD4X9E\nVREvQLNxoBvTs4CrKu1htl3/EOQFPivwnP6WoGxG3sUT0LnwM+BUFAYT66qQUQ7IaUj2bv9iy4oP\nNMKktTkKqegANrZSYdm5Ad1QhxS1NkDHTV4d6aiq6UyUvH8T8gIX8hvVR+1ehZxaFcd2/TFIqWtV\ndKwfVoZ832DBGNS1Rpi0xiDPbyZeegsU55dNO/LoNiAjMvPcU8hIJZhNV8jJS9HzrIwXIEv6rqmH\nR3POc+b1RJauqJg93VVOf++OHjOtVBjXqWeDod+JYj13BzxkLICSgX8JXFuJKeFIkuttNJv3JSq2\nclUMLuZrA5fRVUX138jTma/iXawIk1ZisBvTsFTxrvOsVHh6FXe1P/BXFC64OpppXQbb9X8BuEhF\nZoPAc+bmW6+viSQub0dJuc8BWwaeU+3kzYGAMahrnchIXYUu43oLlPHbXTx0RoquDXltFyIv9sKs\nR2Z5MV3ydd09PkPG88t96dWNbiwyxnXmeSxLJzrmJj82oCmsu9F03GvmQmMwFEd0wd0bGdYZucm3\nkSf3r+Vq0tquvx6qWJgJV3gMOCLwnP5OuEqgpMQLUQ5HJ/B7ZBR91Y/9MhRImLS+ATxcRXWPRhR+\nsRo9VCe0XX9LFIefALYPPOefVepP0diufwTwR3RDu3ENaMfHBWNQD0SiRMDhdBnOmeeOAuKO4xCL\nVMuY8SsPM36l06djF3mT9wd+jm5qQQo/LnBLOV7lyBu+P3ARmoF7HhW6qOYNcKHjNxrdTByDHBcf\nISWQm6vXtcoSXSOGWqkwXwGwUjH/XRVuuRyNwTrkSVy1XX8Uin+egpJ8T4s+isX4Rdr0l6Ob43v6\nsy9FEIexMwa1wWAwGEon0rD+P2RIZ2Qm/4G8ymVpGNuuPw7Fbl8deM7MsjpaWSaixG3omg3cEoWC\nxJowac1AIQn/slLh/v3dn2oQzdzWV6IsexEMR7J1KyC1qLw3WLbrX4NkWZ8FNi8gB+FY4O90HW+G\neGIMaoPBYDCUj+36Q4DvAucgndxFSHrv4gEqufU0qiR7H0pyvgHdWMSaMGmthGYShgJbW6nw8X7u\nUkUJk9a2wPmoXkDehMAqkVEQeRYdF8vMBtuuvztKfF+MZlx6k2LN6J6/j1S0ylLWMVSVbm3OUvSI\nDQaDwTBICTynNfCcy1Fc9V9RDsOFwJORBFdFsV1/mO3621a63SK4LXpuQRfT/Vg2QTx2WKnwfRSj\nDvC7SA2k5gmT1gZh0roPeBSYDnyvD7/bTsiYbkO6zd2FVp4RPf+0AGN6KpJqBN2kGmO6RhkQfzBD\nSaxBVzySoXjM+JWHGb/SicXYBZ7zSeA5ByJ5yg+QBN4ztuufE1VirBQXAI/arn+B7fqVkOQsdvwy\nBvX2qMT0EKRZXQv8CkkfbozCDypBvxx/YdIaFiat65AixS5IAvVnwGZVrIKYzRikVw4q4pJX+cV2\n/c1Q7YO5KEY5l+zxG47CPIYDN9JlWFcc2/U3t11/WO9rxppYnPu6wxjUg5c3GdxJJeVixq88zPiV\nTqzGLvCce1Fi1iVIsvN04AXb9Wf0uGEBREmLIVIaOhF5wdcss9lix+9N4BVkUD0VvXcENXD9tFLh\nQqRRDXBumLSaK9Bsfx1/i1H1wHaUxLqqlQrP6kPVqYtR3PRT6EalO46Pni8PPGdhns8z45dAEo3r\noNCcI6hw0ZYMtutPBR4EnrJdf2w19tFHxOrcl8tAiD02MdQGg8EQA6LCMFciPWeQh+6kwHPKMnoi\nr9+NKBlyEVLb+HMf6lYfhxRO/gjcC6yEikdVtAR0NYjCIa4ErrdS4SP93Z9yCJPWVKDFSoVlJcGW\nwN6oqu5iYEO6Meps15+M6h50AlMCz/mwhzYTwE+BU1C9iZcr2eGsPo1AeQCZEK2DYqD3XsuYGGqD\nwWAwVJfAc54EpqEp8TbkdXsh0uMtp92nkCFzHYrZ3r68nhbNRUg+72W6puWP7OM+lISVCjutVHhY\nrRvTAFYqfK0fjOmv0RW6cQo9e0h/hOyqm3sxpkGG2S9RfYmqGNMRv0PG9KvAD4wxXT0K8ewORXdl\nccV4qEsjDnqOtYwZv/Iw41c6NTF2tuuvC1yPqi12IkWGM8stYW67/n7Ag4HnlFpkpdzxWx6pMSRQ\nCEJvhtNAoyaOvwqRAO5A1QQfQkWI8sZr264/EoUnjQI2DTznmW7a7LPxs11/a1RYphXYKPCcV6u9\nzyoTh2OvLA/1+2i6aKtK9gjJDr2IvBhH53zWhOrXv4Pu3Lar8L4NMY9FqgHM+JWHGb/SqYmxCzzn\nZcAGzoveOh0IIkO7nHb/XoYxDeWP3xxkZNUjpYfBRk0cfxXiEGRMfwUcRveqHkSfjwJm9mBMQ9+O\n317R83kDwJiGmB97hRjUa6FqP78D3kJyMFYF9v0cEkW/kWUD8b+LyktvgBJRrsZgMBgMNUXgOUsC\nzzkd2AbFlm6IlEBOiEqbVwzb9Ve3XX/l3tesCH+Mng8HGvpon4OOMGntHiattXtfsypMRqE+oKIr\nH3S3YlRN9LhosS81sXvjRKTCc25/d2QwUMgJ7XOUvb0JCsxfHXmOy+UV4DV0x5frPt8KeajnAQ8g\n+ZlVK7BPg8FgMPQxgec8gYzpK9EM5IXAw1ESV9nYrt+M5Meet13/gEq02QuLkafMAnbrg/1VjDBp\njQmT1oVh0qqUjF5ViMqmXwW8HCatiuub90IdcuSNAm5Hsfs9sQdKmJ2NCrr0xIpl965AAs9JB55z\nb+A5Rtu6DyjUQ1AH7A78HJ08/lKtDkVsDbyUtfxS9F4h5GoUmuX8y2uQX9PRLBe2bMavvGUzfqUv\nZ49bHPpT0HLgOfMDz/k1Mj4+AbZd/NkHL9uuv38h2/ey3ATMXvx5OAq40Xb9v0QxrdUYv8uBJ1Cl\nPIAfVqD/fbnsACcAv5o+ZMhGJWzfV8ffUbPb2yYghYoX+2B/2cvHoFDTL1hazq679U8A+PKtp28M\nPKejh/anIYm8F1BybRyOh1pa3jnnvf7uz1IUYlD/BgXaH4vu1Kagaa7eeBb4NM/jE0pLIuwpM/XM\nrMf0EtoejMQ6FqkGMONXHmb8Sqemxy7wnLuBdYE7SCRGAH+1Xf93tus3ltHmV8A+LZ+FLvIeHwI8\nu+IO38sXr13u+M2MnrdEyV5fR0oNtcJNSEv5ayePHH1ECdtX/fgLk9Zw4KRo8RdWKuxLZYq1UAIt\nSNbu055Wtl1/EzSrPu/DmX+9rad1UUzzMOAzJP9oKI536ftz37Z02Zc/6mnFQgzbnwLXoEzmGcB/\nUShGpbga+A8SOM9wA/rTZ6ZOngP2BWbl2d6ofBgMBkONERVtORLl5zQiQ/Vbged8VGa7ayO93fWR\nTNifyu1rDvXIw7gO8p5OR4mXp1d4P1UjTFrTgQDdEEy1UuHsfu7SUoRJ6yRUNv1pVAmxrwzqEUgV\nYyNk9xza2wa2698AHAhcGHjOST2sWofCXNcA9qGrAmfFsF1/NLBc4DlvV7ptw/8oSeXjBhTzdhbK\nan4dHWCzkHFbKRJ5OjcTJSaOQnf/Y8lvTBsMBoOhBoniOy9DjpoUCut71nb9Lcps91WkLvI9FLNd\naTqQowlgtej5eyjspCawUuHTKC54CPCXMGnV93OX/kdUiOaYaLEvvdONwM3ImJ5NV5Jht9iubwHf\nQrlgl/Sy+g7ImA7pPc66VM4CXrZdP9bx8QOVngzqaXTVqv82Sk6cjE5+lfixDkBZs/sicfMP0FQL\nKEb7S+QFuBDJ0RgqS774VUPhmPErDzN+pTOgxi4qBrMx8gwuDzxmu/4xkQe71DZbAs+5qpsiFpUY\nvzvRzOo45HBaji6JslrhROAjZDwWczNQ1ePPSoWdqHLgycD91dpPDgngD8AuKBzDQVJ5vXEMUnm5\nJfCc9wpYF5Q8W/EQoSj05OioP8/3snqtEutzX08nrOeRhxpUcvMR4NJo+Tl0FxcHTMiHwWAw1DhR\nDPV5RAleqCjMEYHnVDTW1Hb9cYHnfFGBprYDtgAWIqm0x1C8Zc0QJq2JVir8pL/7EQN+jmJkF6Pf\nNehtA9v1hyNH4Fhg86iaZ3ck0HF9OJKQrOiY267fgPo8Dbgg8JyTK9m+YSlKCvl4B5WrHImmKu7M\n+mxoxbpmMBgMhkFP4DltgeecCOyPjNT/A560Xb9ikqm2668DvB3pYJfriHkUOBtJuy1Es7dTy2yz\nTzHGNKAZ8DNR2Mb+FGBMRxyCjOmnejGmQUbYhej4qMaYH42M6Q+AX1ShfUMB9GRQ/wG4F2VV3oTi\nfkBTMb1NbRgMBoPBUDSB59yEYqDfQomFz9iuv02Fmp+BjKALgT+UoyySxTzkTYelJfQM8WcX4Iro\n9dEUGNscFSX6cbT4myL2V/F48EiD/Yxo8ZjAcxZUeh+GwujJoH4QxTRPQzqMGd6jMNk8Q7yJdSxS\nDWDGrzzM+JXOgB+7wHNeATZFM6NjAN92/R0r0O5ln7/82I8Xf/ZBK7qu3Wu7/phy20Xa1CCv5bAK\ntBdnBsrxtzGKZ65HlQT/2PPqS7EbKnL3HpITLoaKjl9UtGVz4CeB51Qr2TEuxPrYGwixxyaG2mAw\nGAYgUUnnPyHVpyXAXoHnlJ2oZrv+ZshYn4iS3zcJPKe9zGafQp71w5AcbM0RJq31Ua7UQVYqfL8P\n99uMxAgus1Lhv/pglysDTwJfQ2onh1CE99h2/UdQrPVJgedcWJUedu2rCfV3tejxVuA591Zzn4Ye\nKSmG2mAwGAyGfiOqOvd95D1sAu6wXX/PCrSbMX5fAX5XAWM6gVQ/QGoOterk+SkqUnJ9H0vpHYXU\nxC4Kk1a1x24C4CNj+iF0fBVjTO+AjOkF9C7LOAZ5wIvGdv01bNefiRIlXwPuRsmv/1dKe4bqYwxq\ng8FgMMSWwHM6kcF1EdJNvsV2/f0q0O67wMaB5/yl3LaQDNpR0etp1J6EXoajkZTe1sCpfbHDMGl9\nG/hVtFht3emhKE56DVTOfB9U3KYgbNcfhZJQAc6LKnT2xMWost+M4rvKVHRzk0bShg+i3La7S2jL\n0AfU6l10NibkozQycUg1W8K4nzHjVx5m/EpnUI5dpMpxHnAKUmT4TuA5N5TQVLXG72q6Kuu9icqr\nt1V4H1UnTFo7IeOtHdgyKgKTTcXGL0xaB6CEzjqkmPKzKhvUf0CJox8Am6GCQgVju/7lwA9QxejN\nepnZWA6JOTSiUI1MNcqCxy8q0HJnAYb7YCEO575ubc6BYIgag9pgMBgGAZFRfSbgonP/9wPPuarH\njUrbzwzg+SINmSnoQp9RDjkKGXA1R5i0LkS6ybOA9a1UWFEt8Ggfk5A8bzPgAWdW2ZiejuLcO1BC\n4ovFbGy7/s7AA8ijvXHgOS/3ssnpwDlILe0bRffWEFeMQW0wGAyGgYHt+mcgjybAkYHnFKPQ0Fvb\nGwFPAO8DewSeU4w37GLgR9HrT5Bncn6l+tZXhEmrCcUXXwlcWy1DN0xaeyKD/ZfVaD+LeuBpFI5z\nPnBaMRvbrj8aeAlYETgj8JxzC9jfbGAlJM3n99B2Atgo8Jxni+mTod8wSYkGg8FgGBgEnnMOcFK0\n+Afb9Y+tYPNfIs/smsDTtut/vYhtz0ZFXpYgBZETK9ivPsNKhUuAbaxUeE01vcZWKryzD4xpUJhH\npvBJKfu7ABnT/wF+XcD6uyFj+m0UPpMX2/UnoJjuwHb96SX0yxAjjEE9eIm1nmMNYMavPMz4lY4Z\nOyCSKzsmWrzIdv2DCty0x/ELPOcdpOt7OzAauK+IyoofA9sjryTI6J9UYL9iRQ+GdK0df1+jazbj\nx+iGp2Bs13eQEkgrcGiBijAdqOLi5SjeP5s1gDUi6cYXUDjIAmB8Mf0apMT62BsIoRL53O9foGpY\nhtpgLjCuvzthMBhqD9v1jwN+hwyeHQPPmVmhdutQrPbPkXTZOpGxXSh3AHsiyb8jK9GnOBAmrTFW\nKvyyhO0mWanwo2r0qReuAb4D3I88x8VI5I0BXgZWAE4LPOf8IvedyLe/SC3kdWB54N/AgYHnmArU\ntcGgC/kYi76wedTGw9z8GAyGkgg85yLg90hS73bb9VerULudgeecCewLHFKkMQ1KSutElYXXqkSf\n+pswae0AvBcmrX2K3O4wYFakINKXbIOM6SUotr3Y8JXfIGM6QOXqi6W7/f0CGdMBMMMY0wODgWpQ\nGwwGg2HwcDxwH5o2v9d2/YrNeAWec2vgOX8vYdPXUFJfprT1QGAPYBRwS5i0zgyTVrc2RJi0VgyT\n1rFh0noM+DMqyb5+H/UTpLZyWfT6PBQXXzC26+9KV4XOQkM9CuVBFF99ZIXbNfQjAzXkwyh/9E4c\n9Bwz1OLvFafxq0XM+JWOGbs82K4/EpgJbAA8BuwceE6+oh0VGb8opnpM4Dlze1hteWQ4DUNFOp4o\nZ5/9TZi0Ev9saTlnxYb6U1ZtaKwDbgUOsVLhwpz1piGt5gytwGlWKvxtH3b3JJRAOAtYD4XtFITt\n+mNRqEcSODnwnAsq2K81on3MiiqBGgonDue+QSebV4sG2mDG/F4Gg6FsbNe30DR6ErgWeRarolJh\nu/7hwFlIC7u76nV7o4Ivo4AngS0pPuwgdoRJaxfgb+h7PQTsnJ3EGJUtfxN4DrgNuNdKhX1ZnMRC\nMcrDgV1R/HTB2K7/F+AQ9JttbQxfQxbGoDbEGvN7GQyGimC7/jTkqR4GuIHnVFyWLfJO3wnsHr31\nJ+CEwHMW5Kz6bWR4dqIQy32QgVnzhElrKvJQH22lwkfzfF5vpcL+MkRvBvZDSi17F7Oh7frfQOW9\nW4ANA895o4jNf41i+X8DmLjogYkxqPuJDYDvoZP7SsAzaCqyVOqj7bemcl6O0cCBlFbRawEwAlW6\n+kfUr1KIy+9lMBgGALbr74FUNhJIQeGvVdhHHZJhOxcZUbOBgwPP+XfWanUo9GHDaLlmS5LnI0xa\nDVYqjFsMcKai4SJgKirQUxC26yeBZ5HU3omB5/ymiP02oWI+o5CGuQnJGpgYg7ofWB54HpiM7nTr\nUNnTw5H2ZCnshUqm/rQC/cvEIrWiu/H18qyTGcPujPf5wMjo9aXALcAynooCiMPvVSxxiOWqZcz4\nlY4ZuwLIkdPbPvCcTPxyRcfPdv11geuRA+UZYHpOmMkuKGEy46Wu2ZLkEXE+/ppRRcPVUDXEgmXu\nbNdvBB5Bse6PoBj8YjzsGc/2i+hYyG57K8ABzg48Z8Xo7TiOX9yJw7E36GTz4sBu6E63JVruRHqT\nRckN5XA4cGPW8tfRtOPzKF4ww77ohPAIMsIBpgCvIsP3VeAE4B2U/bwqinU7H90AvAZcgU4MK3bT\nXi43Rv0bLLyJOSGWgxm/0jFjVxgXo/PdEOAO2/VXjd6v6PgFnvMyYCNPdb6YbR/4F13X2zPpckTU\nInE+/k5GxvRrQLEJkOchYzoFHFBC3PS+0fNSijCRof4H4CeoEFGcxy/uxHrsas0rmI+CPdRh0ior\nTMJKhcWM13eRxFC2AXop8Ck6oZZCiAzeDhQf+CLyfryFtJznRs9PAtuhE/gjwKaocMpsYCfkRb4P\nOAfFed1Dl4d6SrTeXshYz9feJsg7ne2hHo8y2EvRW61FD7XBYIg5tus3oNLOmfPkloHnfNoPXdkC\nnf+bUBiCS2klsA3dswrwCvJSb08Rs6W26++LDOF2pAv97142yWUIqpI5Bv2+r2e1nVEbmQ2sG3hO\nwWojhlhiPNT9wD3A6jnvbUDO3WsRjEKGdOaueTeUXf1WtJyRbtoFaVzOAT6M1smUwf0QeBh5yx9D\n5XXz8Tkyprtrb9duthmP4rwNBoOh34k0fr+NZvFWRxrVI/qyD7brD7VdPwCm0VUq/YdAQ1/2Y4CT\nQDcszWi2tBhjek2kxAJwUgnGNMhplamqmG1Mr0iXA+1oY0wPbAbVH7pID3O5fIqM57XQH2w00iN9\nJfp8HeACFCaxLzJGz0bhF7+KHtOQ4fss+e+K8n2f3PUS0XtpILtc7CgkKZSP7PKw3bXXHTUvCVUg\ncYjlqmXM+JWOGbsiCDxnvu36u6CQu00XznnrnsYR445qW/DFq9Xed6QGciXwNdv1Dwo851H0u62B\nHBN3VbsPVSCOx99+KEb5SxTOWBC26w9HSiUjkDLIxSXu/wkk1bdCzvsXoevsLYHn+NF7cRy/WiHW\nY2c81NVjB1QOd6toeRtkMG8XLY9AXuAQGApcjuKth0fvv4numjPlbucj72/mJuheYEe6DrBMZbD7\no/cnIS3W7aP3co3vz5BX+WNkXHdHd+3lMj5qr7OHtgYSsY7lqgHM+JWOGbsiCTznI2RwfTZ8+dVn\nTDvhxlMiY7faJNH5cwfgedv1Z6D8FIAf9MH+q0Hcjr8xyHAFJSJ+XMhG0e9/BXJuvY70xMtxCH0I\nPJ3VfjOysRagSp4Z4jZ+tUSsx84Y1NVjEgr7mBwtW8BfkEoHyNBegJJZ1kOGdfb7B6G77ewqXC8i\nOR6QJNCRKPHiBeTtJtrmZ8BfUeb56cgYh2W9x2lkxN+EvODn0+XNzlBoe1OjNgwGgyF2BJ7zJgqV\nW4SKdpzTB/v8ENgIzTROAh6eduJN40gk2lA43Yo9bW8oiLPR2D6J9MAL5UgkGbsQ2CfwnPm9rF8U\ngee0BJ7zTWC9wHPCXjcw1DwDIREsrrJ5vXEVcBjyOl+CpJRAJ4TD0Qnio5xt9kIJgT/poz4Ww2Uo\nxGWwyOYZDIYaJAr/uBude48LPKfUaf5i9tmAYml/AvD+I9c8P+fxv26IKi3+rNr7H8BMR3K0HShE\n8qVCNrJd30b1IRqpkk65YcASy6TES5G0zfPAdcDaWZ81IYPzHRTkv90yW9c2G6BYq3HAoXR5ejdA\nHu0VWdaYBoWCbEtljM816AoXKZdmYH1KM6ZrlUqO32DEjF/pmLErg8BzZn3x2hNnRIu/s13/W32w\nz/bAc36Kwk5mznnib5mPjqX2Ernjcvw1oJCNBKpMWKgxPQE5fxqBS/rBmI7L+NUisR67/vQK7ob0\nOTuQp3Yvusq4/hBpLB+ClCguR3Ju+ahVD7WhC/N7GQyGPsV2/dOQdnQr4EQJg32x30TNuITyAAAg\nAElEQVTgOTPockAcjhIXDcVxIgp1fBdVn1zY2wa269cjydidkWd7RuA5rWX0YSqqqjiTLgUuw8Am\nlh7qe+k6AO9GSRsZtkIe6nmohOhcpH5hMBgMBkMlOB+F22UKv2zQy/oVIUp8+yddOSdnAwnb9adE\nBp+hd1YCvOj10RRgTEecjozpz4D9yjSmQTMMjxKF8tiuP8J2/ck9b2IYqMQlKfE4unSPAbZm6emb\nl6L3CiF3OsAsm2WzbJbNslleajnwnPR/zt3zUiSXNgrwJ22294592J/vRM8T65tH/Ah4ZNHHs5+1\nXX/tbtY3y1pOIAWtYWiW+75Ctrddf8PFn4U/jxYPihIFy+lPPZLrg0j+8KtZz54JvGW7/vF51jfL\nA295KaqpQ/0s+TOY02iKJBM3fADSYd4qz7q523XHmVmv3y6wf4OdHg8MQ6+Y8SsPM36lkxm72MpH\nxZz/jV9n25I0MmwnAttO3Mi5efLOP/hm4Dn/6oN+vIJyhNataxhyKNCYaBiyPvCc7foecH7gOX3Q\njaLp7+Pvmyg8dD7y7veK7fpDgGtI0ABcGnjOgxXox9ZILnYW8ILt+iMWf/r+YSg2++Uetuvv8atl\npkTPfTl220YP0O/dLf0dt/pNFAO1PfB+1vs3ICm3jOj9c8jonpWnDRNDXfuY38tgMPQbtuuPQtcc\nB2gDjgw85899sOudUVjjexsed92GTaOX+xWKqQZ4HClQfNAH/agVRiIxgxWAHyFPda9ENyg/Q+W/\nNwg8Z0EF+vJ7FG5yHnC67fqnRq+fRCXuB0uRs8FGt/ZKfyclXoxE79/J+cwkJQ4uzO9lMBj6lUja\n7tfAj6O3LgRODTynmslmdWhWdWVkzD9gu/6OwDWoKMy+gefcWsX91xq/Rb/PM8BmFJAIaLv+xkCA\nQjRmVGj2oQ4VcpkEbGK7/hvIjpmAElwfqMA+DPEklkmJV6NkkNuQB/o5uqr9/QUVFHkBndQO64f+\nGQwGg2GQEEnbHY8qGLYjFYk7I+91teikS+HjB1E/HgI2RF5yY0x3MQ0lAXYCR1CYMd2Ebk7qgYsq\nGMrTiJIib0ThrUchY/opoBLhJIYaZCB4BY2HujTiFMdVi79XnMavFjHjVzpm7Mqj1/GzXX9b4Fbk\n5HkZ2D3wnHer1J/lgQ/QeXAlYE6V9lMp+uP4a0ShFBsjL/UJhWxku/45SNnjbRTqsaganbNdf1Xg\nDOCmAuKzzf+3dOIwdrEM+agUxqCufczvZTAYYoXt+qshSde1kMzaXoHnPF6l3d2GajGcgZQr5pC/\nuBe2608FXh9kMbq/Ak5GuVbrooTEHrFdfzoywhPA1oHnPFHVHhoGC7EM+RgMbIDixPdD04czqrSf\nepTAUimjdDRwZBnbz0fVLmdWpjsGg8HQtwSe8zbK4XkATec/Yrv+oVXa3RXR8/EoPvjCfCtF8cAv\nANfZrj+ySn2JG7sgY7oDOJDCjOlmFOpRB1xojGlDX2AM6uqxPIqlOgWVOf0tSnipRvGAPVChgEp5\nLMaimLB8JCjMcF8CvMjAKxtvMBgGCYHnfAl8AzlGGoGrbdf3bNev9IzaP5D3dTkUv30gUr/KZTJS\nITkIyettXuF+xI0VgGuj1z8DCjWMf4lmFl4H3Cr0y2BYBmNQV4/dULJCS7TcCfwb2KcK+zocJUeA\n1FHuBJ6n60QEkh18JHrshWKRtgNeBS6Nnv+ILhrnocqUz6FqYpORVNEVyEheMU97+biRLgmogcYa\nGC3lcjDjVzpm7MqjqPGLkhWPQ7N2nciwuybSNq4UHXQlJ2biQy8ip1ZE4Dm3oTjiF9A5+nHb9X9p\nu35jBfvSG311/NUjCd0J6Ibj/EI2sl1/CzQj3AkcGnjO4gr2aQs0g1sO5v9bOrEeu4EQtxrXGOrv\nIs9xtrF5KfApSxeiqQQhMnqbkMG7C/AW8jTPjZ6fRAZ0HTKCN0UJN7OBnVD51PuAc4D3gHuA9aL2\np0Tr7YWM9e7amxetPx/phY5HHoW1eul/HH4vg8Fg6BHb9XdDetXD0Xlv78BzvqpQ8ysgL3UHiqFe\nCc0U/iFPP5qQysTJ6Hy77gDUqz4T+DmKJd8Q+Li3DWzXH4acSasD5wWec3oF+zMa1cJII6m8dYEw\n8JwXKrgPQ/wxSYkAtuuXFRIReE4x47UcMlLXzXrvcST380o5/chhFCrNPhnFau+AdLyzORBpdh4b\nLV8K/Atpc86kq6Ll6eiu/m8sa1D/B32nntq7KVrOGNSgG4hJ9CxxZAxqg8FQE0RxzPeiir8vA7tW\n0Ji9EzlirgMORtrGq9PN+dN2/W2A8YHn3F6h/ceF7YGHotc7AQ8XspHt+hmd6leAjQPPWVLBPp0L\nnAbM3Oj4G74+ZOT42ejaZgee83QF92OINyYpsR/4FMVOZ7yzo5F0T8aYXgedlGcApwKLkNfjBlQ9\nawyqxDQqZ93jkR5nhtwfN98PnW+ddPT4Muv9VqC5m++TnXHeXXvdMZiy0Q0GwwAm8Jz/IofCa8hh\n8pTt+htWqPk/Rc+bo+vCFvTgjAg8518D0JieiK6DCeAsCjemZwDHofE6tMLG9Ep0Ffw5acjI8T9A\nxvSzyNlkMAwI8hlrcTDgdkAhD9+PlndHnoZMkt50JOS/PDKe70AxYzcij/HadBm3ueuOzdnXhyjW\nbjiaksrEGGUK5YxBJ/9JqPrW68BGUV9eymrnRDTF1ozCPjJMyVkvX3vZGeeZLOzx0We9EYffq1hi\nHctVA5jxKx0zduVRkfGzXX+s7fr/tF0/bbv+fNv1v16BvjXQpUm9bxl9S9iuX62E8Goef3VINjAN\nPEZODHl32K4/xnb996Lf4pdV6Ne1UZ/+Zrv+SNv1U9G+9iihLfP/LZ04jF239orxUFePSShsYnK0\nbKEKkBtHy1sBC1C29pfIoz0DhYksBwyjK6Exd925Oft6EVgTWIgSZ36LklYuiD7/EiXS/BW4HoV2\nPIemE3MPjnS035vQ3ff5dHmzM+Rrb35OGwBTozYGIm9ihPnLwYxf6ZixK4+KjF/gOXNREvjfgBHA\nvbbrl1vVtx1pLgNcRe/5J91xFJL5u8d2/a0rrEpSzePvZDSmn6PQwvYCt/s98iI/g2LLK8nKwP8B\nbYn6hjOAPyPn1jNIp7xYzP+3dGI9dgMhbjWuSYm98SekgDEJhVNcjwzhKciTnTFk862by17AJsBP\nqtrj4rkMhb082st6tfB7GQwGwzLYrl8HnI3ia0GJ3W7gOb2Wxu6GBDLSvwW8AdhAUYmPtut/H6mE\nDIveCpBs6x1l9KvabIFyceqRVOG9hWxku/7+yLmzCNgo8JxqGFxbAhvYrn8XCvEYDmxSpX0Z4o2J\noY4ZGyDP9Yp0GcjvA/9Fd+YBXcZ0vnVzuRPYlngZpc3A+vRuTBsMBkPNEnhOZ6QmkZHVOwPwbdef\nWGKTaeAwFGa3JkpQrEPn94JCIALPuRJ5Vn8JfIGM8luIb12Accgorkczq4Ua0yvSpYJyYhUN3CeA\nywLPCVEO0x7GmDbkEicDrFRq1UPd32TikOJwUqjF3ytO41eLmPErHTN25VG18bNdf0eUB7MckAL2\nDzyn1IqxqyJv6FjgcqS6dBvdVFHsoU/DkYzrzsCeFShZXunxSwA3o5jxp4GtUYJ8j0QzA/9AiiD3\nArvXSDl28/8tnTiMnZHNM8Qa83sZDIYBge36K6CQja2Q4sQZwAWB53SW0NzXUX2AzGzyPGRU9KrJ\nXCiROsbRdMWnvgW8GXjO55XaRy8cgG5CFqBZzXcK2ch2/eOB36D8o/UCz6nYmBgMPWBCPgwGg8Fg\nqDaB53yIQivORyEM5wN32K4/rscN8/MAMshBCXqjkJRcJdkU1TD4CXANquj7me36v6/wfvKxPKpj\nAHAChRvT66GKvgDfN8a0IQ4MBK+g8VDXPub3MhgMAw7b9XdHRupYJEW6X+A5xeoWZycpgs6XGyOl\npkr0cRWUELg6XbJkawBnBZ5TULnvEkmg/J/d0Y3DLhQgoWq7fjMKDVkPuCLwnCOq0LfTgRXW+f7F\nL45IrjGZ8pJMDQOLbu2VghIcDAOSOMQi1TJm/MrDjF/pmLErjz4bv8Bz7rZdfxqKEd4UeNx2/RNQ\ngluh8b6ZJMWpyIhMIAWPGVRAwz/wnNnA7Oz3Ipm9xnzrr7b3acfOe/eFFz559v7Hytz1wciY/grV\nayj0u5yFxuFtVDuh0kxC3vrhna2LW1CC/X9QrYhKYP6/pRPrsRsIXkHjoa59zO9lMBgGLLbrNyH1\nimOit24DflBknPKqSPt4DCp+8g16qKJYDWzXTyLjuw7JuZ4bKV8Ui4XKto8GDkVe/EL2nylJ3gls\nGXhOUMK+e+Ny4AejV9144VoHnT0cFXU5tEYSHg3VxyQlGmKN+b0MBsOAx3b9bwFXosqyc5Ch9mAR\nTWQnKX4ThUz0GZFM3fnA/uicvQS4Ajgv8JxUgc0kgPvRd7kH2IPCQj3GoiJmFnBm4Dm/KPoL9M40\n4BkSCdY/8orE0AkrvgBsEXjOoirsy1CbmKREg8FgMBj6k8BzbkZKFo+jhLwHbNe/yHb9oQU2kZ2k\neG3URp8ReM4HgecciEIubgaagB/RlSBYCN9DxvRc4AcUHupxKTKmA1RIp9LUoYqLiUnTv5kYOmHF\nL4F9jDFtKBRjUA9eMoknhtIw41ceZvxKx4xdefTr+AWe8y4qxHUGUu44FnjGdv0NC2ziV0h3eRQK\nI+lr1gg8py3wnG+jwmO3ouqQhTAZSd2Bwl/mFLJRVA3xAGAh8H+B5xRakrwYmoHXgDkrzDjofeDg\nwHNmVWE/5v9bOrEeu4EwzW5CPmof83sZDIZBh+36GwM3oIqIbSgZ7sICNKtXBl5FRuB2wD+r2M2S\nsV1/dOA5mbLp2YVYbkOFXAoJ9UiieOuxwBGB51xRob4lgLo86h2jbddvNlJ8hm4wMdT9xAZoemsm\nsBJKKMmXGV0fvb81xWVtjwYOpKv0arHMByagJI+tS2yjEsTl9zIYDIY+xXb9YcCvgaOit/4JfCfw\nnA962fSnqLR4Cp3DD6UCqh+Vwnb9qeiadyVwfuA530RhG58B6wCfFNBGAsWMOyjuerfAc9K2649A\n37cxejREz/QWW227fgMy5k8Args855JSvp9h0GIM6n5geeB5NMXVgu7OnwIOB17IWXcvpCv60yL3\nMQW4G8Wz5ZL5/j2dYOej5JhLgVuAR4vcf6WIw+9lMBgM/Ybt+rsBVwETUUXEPyOd5de72aQZealX\njpYPBq6vdj8LxXb944DfAaTT6dZ57zxfPz98tb5t/henfPLfe3/dy7YJZHQfhyT1vgDWDTxnTvR5\nEvgwz6YfBZ6zTFy57frjgVmowuRIumLPnw48xy7pCxoGK0aHuh/YDXgWGdMgmZ9/A/uwrEF9OHBS\n1vLXkbdicrTuIeiOOuPBuAS4HSWCrIoE/h8ELkNySjOBzaM+TM+zHSgOKXNQ3IhKz/aXQV2LxFoP\nswYw41c6ZuzKI5bjF3jOvVEFwD8h5YvjgeNt1/8XknK7LfCclqxNWtC5/f5o+VI009mbZ7tcChq/\nwHMusl3/0XQ67SYSiX1Gr7IRo1fZCGD+J/+9d5n1bdefBOwI7BQ9sg3jH2aM6Yj56Pu2o1CZzHN3\nYRoT0Yzu6Gj5DeC3wHU9fYcqEcvjr0aI9dgNBK9gMR7qcqfDihmv76KT4l5Z710KfAqcmbNuiIzn\nDmAYkgbaBXgLaY4mgCdRrFwd8AgqEjAWyQ5lPNRTkEboXkhOaWw3282L1s94qMcDTwBrFfH9Konx\nUBsMBkNEFFt9BArpGx69/QXSa/5T4DmvZa1+K7B39PphYGfkwIkDdcAljSPGHTVypXW/Wvkbx17X\n0Dzi94HnvJG7ou36tyMpwAxtKIzj9sBz9s5dvxhs169D19KJyLP/YhSnvhNSG3mmnPYNg4pYhnz8\nBJVSTaNYqktQ+AJIiucPyBBciGR5uvOextWgXg71ed2s9x5HJ8lXst4bBbyEDGqA/YAdgB9mrXMg\nsBnKBgcZ5v9C8kHZIR9TUEWn5XrZ7qZoOWNQgwz9SfRxoYAIY1AbDAZDDrbrj0TqFkcgjeQMM5En\n+47Ac8YidYph0WfHoutpf9OIQlj+D2hFVRG71dy2Xf814F2UuLgi8GMUH75u4Dlzq9C/EchTvTyq\nOjmzCvswDDxiGfJxCV1akmsg4/MhYDHy7o5FSX2bA1cjY7Fc+tJo+xT4O/L6vo6mmt5maWMa8v84\nvd0gJKL38t0gfFTAdt0Rm4QWg8FgGOwEnjMfFU65IvJa/wA5SraOHott17/rtWtPvWXeuy98BxVa\nebzfOtxFM3Lc7IGcYnsi73m3BJ4zFcB2/Q2QYwjgu1UypgF+BiSjfT1RpX0YBhH9qUM9L+v1uOi5\nNXreCt3ZzkNC9nNRrHAtsQMSid8qWt4G3UBsl7PefKTykbm5uRfFkWVihcaiGLkdkQc5iWSH7kee\n/VE99KG77WDpGOrxwOfEZ5qwFoi1HmYNYMavdMzYlUdNjl/gOf8NPOcIdC4/AnlUhwLfXvOgs77T\nPG6FTqBpxAprHWW7fn0Vu9Lb+I1Eyhx7oGv3DvRiTGeISrRfh7zblxVZRbIY1kIx6mmUP9SX176a\nPP5iQqzHrr+TEn8InIak23agK9xga5ZWvHgpeq8aIuvVYhKKb34oWrbQd72dZcNXXkQ6pK8Ai4Aj\nUcKEBfwXOAzdTf8VnQBOR4Y4yAvwLJomu4ylvcxf9rDdm3SdRKZGbRgKJ5ZJETWEGb/SMWNXHjU9\nfjle65WA/evqGw+csusxG7x+/eks+vid77fM/Wh32/VvQNUM/1OArnUx9DR+45ExPR3Nlu6Mrt+F\n4qEQxreBU0rtYC8kgIuR0X4lXd7wvqKmj79+ZtCO3bMo7CH38QnLhjRsh4y/TAbueywd4nEtUrrI\nR26Ywho57+XezcRxeS+WDn/p6/1fRpfnvD/2X2u/l1k2y2bZLMdq2Xb9dYYuN+UtID12rS3Stuun\n1z/6yrTt+qHt+pfYrr9dw7DRU6vYnyQyeNLAO2hWuZj+b7X+0Vd22q7fYbv+ZlXoX2Z5MkrwnAvk\nSubF5vc0y7Fd7jY0tpoe6mm9r/I/HkV6mpsij+7jwPooQQF0x9qTWPuZWa/fLmK/ceFOJJvXH4l5\nQ9BYH9XbigaDwWCIJ4HnvILyj+6f+/q/R3z8zD13jZy8ng2sgMp8H7P2Ib+eO3S5lW5HlQofDjyn\nUrtfBc2SroISJHdCOtG5BkleouTLa9E18JzAc56qVMfy8F7Ur/WQUW0w9MS20QM0A9Mt/amssDYy\nokHG9x3oz9iOQkG+jrzSmyMNzindtBPXwi5xJ3Oii8MUSi3+XnEav1rEjF/pmLErj4E+ficCFwDv\nA/dP3Hi3/668249WQTUQVs9abwFyYC1BsdiZR3POcgMKx2wH2hd98m4C6Bg2ccpioH1B6s2G12/4\nyaodi+c31zePeHO5DXbc9aPgjh7DMyMZu6kox2hrpLJhoWJoduA5rT1sXusM9OOvmsRh7GIpm3cT\nqoTUhmKH/0xXbHEzks3bFmUIH4s0lPNhDOrax/xeBoPBUBkaUbGvdaLlu4E9ouqDayPN6r2BDcvd\nUcvcObxy5XG0L57HyMnrscb+Z9LQNBw0U/w0ik/+D7rGr40M6K2ALekSI8gQAk7kaTcY4kosDepK\nYQzq2sf8XgaDwVA5ZgD/zFqeiuRb/4ft+qsgw7YDydV292hHSlQZNaoGoH7Bh2+MfO2aU27obF+y\n+pDRE19a74jL/tXQPGIamnFuKqCPIVIqeTx6vFzh5EmDoRoYg9oQa8zvZTAYDJXlcqRbDZoB/n4F\n204g9ahvo+n36cBXALbrN6KCZptG72+KvNNv0GU8zwTeDzynr2ofbA38m/4pXGYYWAxoeyXfH9IU\nKOmdNSgwYaQPqMXfK07jV4uY8SsdM3blMVjGbxSqNJhGXuZJFWp3DeC8qN35yFiOM5ug7/8ECofp\nbwbL8VcN4jB2tWivFIwxqGsf83sZDAZD5fk6XVV1f1yhNndEnt40knyNM03Ay6ivv+nnvhgGBgPa\nXjEGde1jfi+DwWCoDn9H59gnURx0OayMquqmgV+W2VZfcA7q65vAsH7ui2FgMKDtFWNQ1z7m9zIY\nDIbqMAZpQqeBE8poZxhSD0kD91K+cV5tNkWe9E6UfGkwVIIBba8Yg7o04hCLlKEWf684jV8tYsav\ndMzYlcdgHL/d0Hl2MaV99wRwfdTGu8DGFetZ9bgC9ffC/u5IDoPx+KsUcRi7WrRXCibOBvUGwMXA\nfkhsf0aZ7dWjDOlKZZiOBo4sYbspwEt53m9G2dvFEpffy2AwGAYq16Bz7UygrshtfxxtuwApeNQC\ndcBhqDiNwVApBrS9EleDenngY2Rkgv7cTyMju1T2As4qs18ZEqgyZT7DuJztLgW2K7LNOPxeBoPB\nMJAZB8xB59tji9huO6SSkQb2rUK/DIZaYkDbK3E1qL8P3J/z3u8Ar4w276NLomgKKt3+Z+A14Od0\nienviBJRngTOyNp+SrTuFcggfgRYhOLizgcms7ShfFLU7uSc7bahq7pl7r63BG4s8nvF4fcyGAyG\ngc5J6HzbAqxawPorAZ9G25xbxX4ZDLVCt/ZKsdM+NU+YtNL5Hr2tX8KuOtBJK5tGlCBRKusjcfwM\nawH3oBKy6wPfQIkjpwIHI+N2XcDO2mZNlFCyDzLuZwEbRdvkK5CTztluPeB9VHkrd9+galzTyviO\ntUIcYrlqGTN+pWPGrjwG8/iNjp6bkEMknw0wBTgane/fACYADwI/jT4fzONXCcz4lU6sx66hvzsw\ngLkHxU1nswFwRIntjUJGenalp6+A26PXfwUcZACvjbzToJCT7YAgWv4cuDN63VrAfjNGdvZ23e37\n1mi98SjeeyBXpXqzvztQ45jxKx0zduUxmMfv98DJyKCeAfwQ+BOwFbArSl6cmrPNo8ABdJ3P4zh+\nI5FT6Cl07Y0zcRy/WiHWYzfoDGorFRaV0Ffs+ll8isIu1kJe29HA2yhUAmAd4P/bu/soO+rygOPf\nTYJEEAKhaiDRBFJeBYKh5aUksBDEeETkRQ5a1HPAQkGLpy2cQ0EKyJsCnqIIUd5TsVCKQKFAgJ72\n0CJSrCQQEEJMeFNpIqUIYgUxpH88c7t3787c3b2zd2fmzvdzzp69M3fmzm+f+e3uc3/3md98DTiO\nqEvbAjif+BjuouRrLvBvwFJGfrvLicRowjEZz69ps+8bDJRukLTpzRHsl8YyDkkql7XAd4j/OxAz\nYHyFGLBpeI34H3I3UbY42r/942kSUV75ZeA9RMK1hN4ezJG6pqw11AuIhPRPkuWPAdsycMHensTo\nAMBpwIzk8R7A8cRFjZsBmze95s8ZeBM0iygfOZRIgm8GDidGpFcxMMowlaiDa+zTXCP9TuD5lnav\nAqYn+/0YOJOhtdVZxyb5mVcwOmU4X5JUBzsQf3MbdztcT1yPczHQTzluzz2cPqLM8EkGfoYHgb2K\nbJRqwRrqAkwjPnqamSzPABYzMH/nPGIKoj2JuuSftaw/Gvgl8ErTay4napkbVgCHAI8St1e9ixhl\nPo64i9VyYqRhWtM+jc6wHfA+4CZiBPzCZP0ZxMjE7cD9Kfs1HqcdGyKRXzokGr2n1LVcFWD8Omfs\n8ql7/FYAdxAXpJ9HzNq0E1EKcj/w1jD7lyF+k4BvEP9vVhOf8s4jSj7Krgzxqypj12VlHaEezrXJ\n94nAoqb1VyXfpzHUYURZCGTPBV20RThtniSV2UziE9AqO4SYH/sdRTdEtdLT+UoVE+o5wL1EWcXn\niLmbG+vvI0aO00wgPtbqIxLq5V1t5ehNJm48M1plP1+SJEk9na9UMaHWYJ4vSZJUdtZQawhrkfIx\nfvkYv84Zu3yMXz7jHb8+orTjveN4zG6y/3XO2HWZI9TV5/mSpOJMJKZyLaODiP8RqxnZ1LFSN/V0\nvmJCXX2eL0kqxmbEjb9eIaZrLZu7if8RpxXdEIkez1dMqKvP8yVJxehjIGm9peC2tGrMmf0b4h4H\nUtF6Ol8xoe5MmWqRqni+yhS/KjJ+nTN2+Ri/od4P/Ir4W3zEMNuOZ/wWEW26YpyONx7sf50rQ+yq\nmK+MmAl19Xm+JKlYXyD+Fv8Xg+/QW5R3Ab8m2rRTwW2RGno6XzGhrj7PlyQVawJxH4Ey1SvvAvxl\n0Y2QmvR0vmJCXX2eL0kq3vbA53FKXSlLT+crZU6o5wCXAkcCJwP7dek4E4mRhdFMKZRVizQFOLGD\nNswi+1bok4EH2uxblvM1GmWo5aoy49c5Y5eP8cvH+OVj/DpXhtiVOl/pB9YR9VsNGwLXAs8CTwD7\nt9m/rAn1lsBaIpmEeMf/QyLJHmuHAeeNwev0AduQnRjn2e9yss9jGc6XJElSO6W9U+ImwFeBu1rW\nH0NcFDGHGNm9bpzbNRY+CiwF3kiW3wZ+wPBXUHfiOOAGYpT4SeAa4CngLOLNScOBwM3AQ8DpybpZ\nybZXAsuBq4HZwDLgQmAmgxPlU5LXndmy3wyiP2Ud+4aknZIkST2l6IT6b4CLgJdb1s8jRqhfA+4l\nJpyfPb5Ny20dA8l0wwZEYj3WdgWeTh7vANwJ7JasPzhZvxFwKvAZYB9gZ2DP5LntiTc1uxBvZlYD\nH0y2by0jWc/AO7Tm/V4Adsw4NsAKYG6un1KSNN62ZPDgSLcdS/wvkSplUoHH/ggxQn0r8DEGD6PP\nB85oWn48Wbc670H3PPOe1OH6h89ZmFp/3Ng+6/k27iRG15vNAf50lK8znE2J5H1dsvwqcFvy+EZg\nITFZ/0eIqYceatrvSOAy4g3N7cn6kfycjW2a92t37Ma2WxD13uuovkYd18pCW1Fdxq9zxi4f4zdy\nRxCfOn4H+GKyrpvxm8HAnNNbAS914RhFs/91rtSx6+YI9VLil6H16xfAVKLU42+esPYAAAz4SURB\nVKRk2z6GT+RGWmfbWrA+2gL24bYf7vUbyy8R5RULk+UpwCrgrWT5A8To7tHESPD/AhsDdwAHEbeD\nvYwYKW5sux9wAYNHemczOHat53TTpvX3AUclrzmbKN/YGljTtP3WLftvxeDRiW2bHq9h5PFo+P1R\nbl/W5ZXJV1naU7Vl49f5ciN2ZWlP1ZaN38iXIT7dPIlIrrsdvzOIgb5bif+hRf/89r9yLcPgZLoM\n7fl/3Ryhbvfx/jxgGnGRHsDvEeUBmxMX132fKBl4Lnl+F+DLbV7v7KbHq9o16uFzFm7P0BOS+W4n\nZfuRWkAkxCcA9wD7AucDeyWvtzExuvsckSzvTZSI9BE1yFsRCe/7ieT69mS/14mLNRt+TYz6Ns7l\nJsChwBIieb4/WX8n8BXgpuR1phIT57f6bwaS8MbyBGA6cfvXBcD1GT/zFKJO+/nk2Dc2PbcFMUrt\nBYiSVA2PE/+HvkGUYR5O90YHNyL+bwB8vUvHkEarP/mCyGNK7zpi7suGE4jSgU2BDzOQWKcp6ywf\nRxPlFecmyycCDxJ/nCAmq//jpuWrgAOIC/fmAn/Q9Fqt27ZaQoxizyIuSryW9AsD9we+R1xE+CNg\nj2Sf5S2vdxHxCcOFyfIngceIqe8uB84kLkps3m/mMMeeR1yYmKYM50uSNFQfUbq3HniE7tVTH58c\n42FGNwWsNJ5Kn6+0JtSTk3WNafMOaLNvWRPq4VyVfJ+WfP8usDvxsdpfMfgPSuu2rQ4jRr9bZ+Ro\nZzuG+fhijC2it6bNG+/49Rrj1zljl4/xG73NgGeA3wKfpjvxe5T4X/CpLrx2mdj/OleG2FUxXxmx\nKibUc4h65vc1rbsg+d7P4MQzbdtWE4jR720YOtpcBpOJMp4sZT9fklR3c4lPNbtlBvDXxGxYUln1\ndL5SxYRag3m+JElS2ZX2xi6SJElSpZlQ11cZapGqzPjlY/w6Z+zyMX75jEX8NiZKPOrI/tc5Y9dl\nlnxUn+dLkqrnCGIqvdHYjrh4/jFiqjypSno6XzGhrj7PlyRVy3zib/drDL7hVzuHJ9uvB1YQ07ZK\nVdLT+YoJdfV5viSpWvqAfyD+fq8jEuSbGXwPhYZJwMXJtuuT7TYZn2ZKY6qn8xUT6s6UqRapiuer\nTPGrIuPXOWOXj/HLpzl+U4hpXX/HQLI8L2WfI5Pn3gL+nHrfuMX+17kyxK6K+cqImVBXn+dLkqpr\nMrAb8BnSR577gEtIT7alKunpfMWEuvo8X5Ikqeych1qSJEnqBhPq7poDXErUj50M7Jex3UTi1txj\nUVc2BThxBNul1SLNIqYzajUZeCBfs3pOGWq5qsz4dc7Y5WP88jF++Ri/zhm7LitryceWwFoiGYV4\n8/JDIsludRhw3hgcsw/YhvSkOO++lwP7d9iu4ZThfEmSJLVjyUcBPgosBd5Ilt8GfkBMhN/qOOCG\n5PEs4EngGuAp4Cxgw+S5A4nphh4CTm/a/ingSmA5cDUwG1gGXAjMZHCSfErymmn7ziD6RNqxb0ja\nKUmSpCZ1TKjXZ3wNt/1orWMgmW7YgEisW+0KPN20vANwJ3HV9K7AwcQdpU4lrqLeB9gZ2DPZfnvg\nLmAX4BhgNfDBZPvWMpLWn6V53xeAHVOODTHH6NzsH1eSJKme6phQj5c7GXr3qDnECHOzTYnke13T\nuleB24A3gRuBhcnXTsTo9CNEwtwowXgZuD15PNI67O2ArVv2zTp24xhbEPXespYrL+PXOWOXj/HL\nx/jlY/w6V+rYTSq6AQUY7YV/nV4o+BKRPO9AjO5OAVYBP27Zbv0IjzGRmED/mJb1s4A1bfZ7g4Gy\nDYik+E1gJfDbYfZNY71zWFl0AyrO+HXO2OVj/PIxfvkYv86VOnaOUHfPAuAyBiay3xc4n6EX9v2K\nSJab39xMAQ4lEuGjgCXAPwHziZIMgKnA+1OOu5YY9W5YQ5zn6ck+H6d9Upx2bIhE/GXSS1YkSZJq\ny4S6e6YRZR8zk+UZwGJg95RtlxO1zA0rgEOAR4EniBrnN4iLAs9Ntr8vOQYMTpB/A9xEXBB5YbLu\nDOBuorTj/pZjr295nHZsiER+adYPK0mSpOoq67R5o3EYMXoN2XNBj7XR1iItwmnzmpW6lqsCjF/n\njF0+xi8f45eP8etcGWJXxXxlxHohoZ4APEjUUs8iRqDLZDJx45luqdr5kiRJ9dPT+UovJNR15/mS\nJEll541dJEmSpG4woa6vMtQiVZnxy8f4dc7Y5WP88jF++Ri/zhm7LrPko/o8X5Ikqews+ZAkSZK6\nwYRakiRJyqHIhHox8FNgWfJ1WtNzGwLXAs8SNxcZ7fzHrxDD8n5V4+uV9NNYatZy5WP8Omfs8jF+\n+Ri/fIxf54xdhuuAz2c8dwJwG3EL7Q8Dz7V5HetvO9dfdAMqrr/oBlRcf9ENqLD+ohtQcf1FN6Di\n+otuQMX1F92ACusv+PilraHuy1g/jxihfg24lxjBnD1ejaqR/qIbUHH9RTeg4vqLbkCF9RfdgIrr\nL7oBFddfdAMqrr/oBlRYf9ENyFJ0Qn0y8AhwAbB10/r5DL799uPJOkmSJKlUuplQLwVeSvn6RXLc\nLwHbAPsSI9GLhnm9kZZ2tNbXuJy+vB2wRYnaU7Vl45dv2fh1vtyIXVnaU7Vl45dv2fjlWzZ+nS8f\nxEDsytCeQbJKLsbbO4D/AaYDrwJ/B9wE3JE8vwz4BLA6Zd9VWA4iSZKk7noM2K3oRrSannyfAHwO\n+Nem50ZzUaIkSZJUS/8MLCdGn68Edm56bjIxC0hj2rwDxr11kiRJkiRJkqT6mQrcDawA/hO4BNik\nZZuzgLeBnTJeow+4iKi5fhr4ZFdaWk5jEb+zgbUM3HTnm91oaEm1i99zwFMMxOVDGa9R1/43FrE7\nG/teWvw2AC4GVhKf7F2R8Rp17XswNvE7m3r2v6zYzWIgFsuI3+OXM17DvpcvfmdTz74H7X93DwTu\nIWJyD9k3+qtz/8u0OTHzB8BEYj7qk5uen0sE/hmyE8KFwIPANOAPgReBjbrR2BIai/idRXTMOmoX\nv2fJjlmzuva/sYidfS+0xu9LwPVEnwJ4T8Zr1LXvwdjEr679b7j/Gw2XAJdmvIZ9L3Qav7r2PWgf\nv58A+ySP9yWS5TSF97+i56FO8wrw78njdcASBkazNgQuA06k/Qwl84G/B9YQ73YeBvbuRmNLaCzi\nxwie71Xt4gcji0td+99YxG402/WadvH7OHAu0acgph9NU9e+B2MTP6hn/xvudxdiNq6jiWQnjX0v\ndBo/qGffg/bxexV4NxGbdwO/zHiNwvtfGRPqZhOAPyNm/AA4hxhleH6Y/ebhjWGg8/gBfAp4FPg6\nsEtXWld+rfEDuBF4APgL4F0Z+9n/Oo8d2PdgcPy2IubsP4742PMKYE7Gfva90Gn8wP6X9rsLcAjw\nMyI2aex7odP4gX0PhsbvE8C3iaT7m8CRGfsV3v/KnlBfRtQbXUG809gd+FbT8+3ezbU+N9Ibw/SS\nTuP3bWAmsAfwI+JdXx01xw/iF3ZX4hd6b+ALbfate//rNHb2vdAcvw2JGsMXiY8y7yLm6s9S974H\nncfP/jf0d7fhWNqProJ9DzqPn30vNMdvAnA/8an6VOCkZDkrd7H/ZfgaUYC+QbJ8KvBzohbzWeAt\n4t3egSn7ng98sWn5H4EFXWtpOeWJX6tV1O/dcmv8Ws0n7gaapu79L0/sWtn3wnMMrgd8gfQbWtW9\n70G++LWqW//L+t2dDrxO1Lpmse/li1+ruvU9GBq/HYnf1WY/Jf2Ohfa/DOcD/0LMR52l3UVOhRen\nFyxv/KY3Pf4QcVFAnaTFbyNgSvL4ncDVwJkZ+9e5/+WNnX0v/Xf3b4HPJo/3ImZMSVPnvgf541fn\n/tfu/8bpDD9iat/LF7869z1Ij18fMWvHXsnyPEp8UWIZfYCY0q15iq1bUrZrTQjvImawgDgJFxMn\nYiVRl1QXYxG/xUT90aPAd4lOXBdZ8duaGFV9DPgP4KsMHm2w/41N7BZj30v73d0KuBd4Mlm3T9N+\n9r0wFvFbTD3733D/N54GDkrZz74XxiJ+i6ln34P28VsA3ErE5Ragv2k/+58kSZIkSZIkSZIkSZIk\nSZIkSZIkSZIkSZIkSZIkSZIkSZIkSZIkSZJUpD8i7hT4DHG3wFuBTQptkSRJklQRWwJLgInAe4Hr\nk/WzgU8DpwIbFNM0SeoNE4pugCSpq44FvgWsAzYFfpesPx74HrAMOLiYpklSbzChlqTetjlR6gHw\nWeCa5PFWwBvAWmDbAtolST1jUtENkCR11RXAUcCewBPA95P1ryffJxCj15KkDplQS1Jv+wlwXsr6\nF4na6enJNpIkSZJGYTpwAHAKccGiJEmSJEmSJEmSJEmSJEmSJEmSJEmSJEmSJEmSJEmSJEmSJEmS\nJEmSJEmSJKlw/wevEK5B1vbWOAAAAABJRU5ErkJggg==\n", | |
"text/plain": "<matplotlib.figure.Figure at 0x10dc08e90>" | |
}, | |
"metadata": {} | |
} | |
] | |
}, | |
{ | |
"metadata": {}, | |
"cell_type": "markdown", | |
"source": "## Discussion ##\n\nThe transformations rates $\\Omega_{HF}(\\sigma_0)$ and $\\Omega_{FWF}(\\sigma_0)$ quantify the rate at which surface water is made denser (positive values) or lighter due to surface forcing (heat flux and fresh-water flux respectively). Their sum $\\Omega(\\sigma_0)$ is therefore directly related to the mean overturning circulation in isopycnal coordinates (modulo the uncalculated interior mixing terms). The value of computing transformation rates is that it allows you to understanding the connection between surface thermodynamic forcing and the interior circulation. The mean structure of $\\Omega(\\sigma_0)$ reflects the well-known sturcutre of the Southern Ocean overturning cells in density coordinates.\n\nThe wind perturbation seems to have the strongest effect on the heat-flux-related transformation $\\Omega_{HF}$. In the ACC densities (~$\\sigma_0 = 27.0$), this transformation goes from near zero in the control to significantly negative for the perturbation. A negative $\\Omega_{HF}$ corresponds with __stronger surface heating in the ACC__. This is constistent with what we expect when Ekman pumping is increased, more cold water being exposed to the surface. The net transformation rate $\\Omega$ becomes more negative, corresponding with a stronger upper MOC cell (approx. 8 Sv). It is possible that these changes are canceled out by increased interior mixing and do not actually correspond to a change in the isopycnal MOC. I haven't calculated the interior mixing.\n\nThere are also some interesting changes for the densest water. Increased negative heat fluxes (presumably in the Antarctic zone) appear to be driving more porduction of very dense water. This could correspond with a slightly stronger AABW cell." | |
}, | |
{ | |
"metadata": { | |
"collapsed": false, | |
"trusted": false | |
}, | |
"cell_type": "code", | |
"source": "", | |
"execution_count": null, | |
"outputs": [] | |
} | |
], | |
"metadata": { | |
"gist_id": "dae4ebccaabbd72c424c", | |
"kernelspec": { | |
"name": "python2", | |
"display_name": "Python 2", | |
"language": "python" | |
}, | |
"language_info": { | |
"mimetype": "text/x-python", | |
"nbconvert_exporter": "python", | |
"name": "python", | |
"pygments_lexer": "ipython2", | |
"version": "2.7.3", | |
"file_extension": ".py", | |
"codemirror_mode": { | |
"version": 2, | |
"name": "ipython" | |
} | |
} | |
}, | |
"nbformat": 3, | |
"nbformat_minor": 0 | |
} |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment