From 3482d91de5840a51ddab06fb3cd4b04a539b6d72 Mon Sep 17 00:00:00 2001 From: Mark Bakker Date: Wed, 7 Feb 2024 16:59:58 +0100 Subject: [PATCH] fixed large diameter well --- notebooks/large_diameter_experimental.ipynb | 130 ++++++++++++++++++-- timml/well.py | 49 ++++++++ 2 files changed, 167 insertions(+), 12 deletions(-) diff --git a/notebooks/large_diameter_experimental.ipynb b/notebooks/large_diameter_experimental.ipynb index 337be76c..02e88320 100644 --- a/notebooks/large_diameter_experimental.ipynb +++ b/notebooks/large_diameter_experimental.ipynb @@ -11,7 +11,7 @@ }, { "cell_type": "code", - "execution_count": 1, + "execution_count": 2, "id": "ed2bfa2a-5ffc-48bc-8224-66edf4410140", "metadata": {}, "outputs": [], @@ -23,7 +23,7 @@ }, { "cell_type": "code", - "execution_count": 2, + "execution_count": 29, "id": "0c7a38da-e288-4fb4-9d73-5cfea442d38c", "metadata": {}, "outputs": [ @@ -39,21 +39,21 @@ ], "source": [ "rw = 5\n", - "ml = tml.Model3D(kaq=10, z=np.arange(20, -1, -2))\n", - "w = tml.LargeDiameterWell(ml, Qw=100, layers=[0, 1, 2], rw=rw)\n", + "ml = tml.Model3D(kaq=10, z=np.arange(20, -1, -2), kzoverkh=50)\n", + "w = tml.LargeDiameterWell(ml, Qw=100, layers=[0, 1, 2, 3, 4], rw=rw)\n", "rf = tml.Constant(ml, 100, 0, 20)\n", "ml.solve()" ] }, { "cell_type": "code", - "execution_count": 7, + "execution_count": 30, "id": "dc1a5d30-8a1a-46d6-bfd4-67e9a360dc4c", "metadata": {}, "outputs": [ { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAj4AAAGdCAYAAAASUnlxAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAACSmElEQVR4nOzdd3zV1f348dfdN3vvPYBAAmGDgAwZAoriVlRQcduiorVa9efCWm2LVlvU+tXauusCBwiIIrKH7DASkpC9k5vcvT6/Py4EUjaGkeT9fDzysLnnfM4973AJ756pUhRFQQghhBCiC1Cf6w4IIYQQQpwtkvgIIYQQosuQxEcIIYQQXYYkPkIIIYToMiTxEUIIIUSXIYmPEEIIIboMSXyEEEII0WVI4iOEEEKILkN7rjtwNnm9XioqKggKCkKlUp3r7gghhBDiJCiKQktLC/Hx8ajVv27MpkslPhUVFSQlJZ3rbgghhBDiNJSWlpKYmPir2uhSiU9QUBAARUVFhIeHn+PenB6Xy8WSJUuYMGECOp3uXHfntHWGODpDDCBxnE86QwzQOeLoDDFA54mjoaGBtLS01n/Hf40ulfgcnN4KCgoiODj4HPfm9LhcLvz9/QkODu7QH+LOEEdniAEkjvNJZ4gBOkccnSEG6FxxAO2yTEUWNwshhBCiy5DERwghhBBdhiQ+QgghhOgyJPERQgghRJdxyonPihUrmDJlCvHx8ahUKubPn9+mvLq6mltuuYX4+Hj8/f2ZOHEi+fn5x23T5XLx7LPPkpGRgdFoJDc3l++++65NnaeffhqVStXmKzY29lS7L4QQQogu7JQTH4vFQm5uLn//+9+PKFMUhalTp1JYWMiCBQvYvHkzKSkpjBs3DovFcsw2n3jiCd58801ee+018vLyuPvuu7niiivYvHlzm3rZ2dlUVla2fm3fvv1Uuy+EEEKILuyUt7NPmjSJSZMmHbUsPz+ftWvXsmPHDrKzswGYN28e0dHRfPTRR9x+++1Hfe69997j8ccfZ/LkyQDcc889LF68mL/+9a+8//77hzqr1coojxBCCCFOW7ue4+NwOAAwGo2tr2k0GvR6PStXrjxm4uNwONo8A+Dn58fKlSvbvJafn098fDwGg4EhQ4bwxz/+kfT09OP252CfAJqbmwHf1NrBMwE6moP97qj9P6gzxNEZYgCJ43zSGWKAzhFHZ4gBOl8c7UGlKIpy2g+rVHz55ZdMnToV8HWsW7duDB48mDfffJOAgADmzp3LY489xoQJE1i8ePFR25k2bRpbt25l/vz5ZGRksGzZMi6//HI8Hk9r4rJo0SKsVivdu3enurqaOXPmsHv3bnbu3ElERMRR23366ad55plnjnj9ww8/xN/f/3TDFkIIIcRZZLVamTZtGiaT6VcfQNyuiQ/Apk2bmDlzJlu3bkWj0TBu3LjWC8UWLlx41HZqa2u54447+Prrr1GpVGRkZDBu3Dj+9a9/YbVaj/qMxWIhIyODRx55hNmzZx+1ztFGfJKSkqisrDxmsnS+c7lcLF26lPHjx3f4Uzg7ehydIQaQOM4nnSEG6BxxdIYYoPPEUV9fT1xcXLskPu1+ZcWAAQPYsmULJpMJp9NJVFQUQ4YMYeDAgcd8Jioqivnz52O326mvryc+Pp5HH32UtLS0Yz4TEBBA7969j7tjzGAwYDAYjnhdp9N16A8AdI4YoHPE0RliAInjfNIZYoDOEUdniAE6fhzt2fczdo5PSEgIUVFR5Ofns3HjRi6//PITPmM0GklISMDtdvP5558f9xmHw8GuXbuIi4trz24LIYQQohM75REfs9lMQUFB6/dFRUVs2bKF8PBwkpOT+fTTT4mKiiI5OZnt27dz//33M3XqVCZMmND6zPTp00lISOCFF14AYN26dZSXl9O3b1/Ky8t5+umn8Xq9PPLII63PPPzww0yZMoXk5GRqamqYM2cOzc3NzJgx49fEL4QQQojzUElZMQuWzadQ7aJE034TVKfc0saNGxkzZkzr9wfX18yYMYN3332XyspKZs+eTXV1NXFxcUyfPp0nn3yyTRslJSWt634A7HY7TzzxBIWFhQQGBjJ58mTee+89QkNDW+uUlZVxww03UFdXR1RUFEOHDmXt2rWkpKScaghCCCGEOE84XU6+XvQZm+tKqQj2pyIglAp9NLWqKJTk0QB4LeZ2e79TTnxGjx7N8dZDz5o1i1mzZh23jeXLl7f5ftSoUeTl5R33mY8//vik+yiEEEKI88/e/J18tfI79uuhMjiYcmMk5Zo4nEG9IKjXEfWDlSYSXFVENlfwWTv1od0XNwshhBCia3O4nHz13Rdsri+lIsiP8oAwKnQx1KsjIXXsEfV1ipN4byXx9lrim5uJabSS5R/K5KtuJjA4hPr6ej7joXbpmyQ+QgghhDhtpeUlfPnDlxSo3VQEBVPud2AUJzALArOOqB/hrSPBVUW8pYnYxhaSbF4uGj6Bnv1PvAmqPUjiI4QQQogT8nq9rNywgh93baDUX095UBjl+hhq1DGQOOqI+nrFQbynkoQDozixjRayQ6KZeM0MjEa/cxCBjyQ+QgghhGjD7nKwYPEXbKovozw4gHL/CMp08ZhVoZAy/oj6YUoDCc5K4i1NxDU2k2TxMGbYWLIHTT3rfT8RSXyEEEKILqy2qZbPF3/BLreZ8uBgyoxRlGvicQX0hICebeqqFQ9x3koSHDXEt5iIrzeTqQ3gkmunExJ20TmK4NRI4iOEEEJ0EYXlxXz541cUaDyUB4dSaoihSh2DEj3kiLpGxUaip5wEax0Jpmbi623kpvRgzJSr0eo6bvrQcXsuhBBCiGPKK9jBl2uWUmRQ88bSjyg1xFKnjoaEkUfUDVaaSHJVkGhuIL6xmYRmFxcOH0PukKvPQc/PLEl8hBBCiA5u8+4tfLPxR4oMWsqCwinVx9GoCofEMUfUjfTWkuSsJLGlkbj6FlIcKsZfdjXJ6aPPfsfPAUl8hBBCiA5CURS27NnKN5uWU2jQUBoYQak+HpMqFBLaJjkqxUusUkWio5pEUxPxdc2kq4xcfO10IqOPXKDcVUjiI4QQQpynthXu4ut1SynQqSgNjKBEH0+zKhTiR7epp1I8xHsrSXTUkNTUSFydmUx9AGOvuolVa9Yx+eqbOvTt7O1JEh8hhBDiPJBfXsQXq75jr8pNaVA4Jfp4mlRhENt2Tc7BnVXJ9ioSTCYSalrIMAZx8Q23EBIa0aauy+U6myF0CJL4CCGEEGdZRX0VX/y8kO3OFkoCwykxxPmuc4i6oE29gyM5yfYqEk0m4mtayND5MX7arYSFR5+j3ndskvgIIYQQZ1CzzcKClYtY11DG/sBQSoyxVKtjIaR/m3oH1+Qk2atINjUSV9NCqqJj3A03ExN76TnqfecjiY8QQgjRTlweNz9sXcWyvVsp9Pdjv38M5eoEvNpMiM5sUzfSW0OKo4Kk5kYSqltIsrkZdcXVpPWYfI563zVI4iOEEEKcBkVR2FG2l683LCdP5WF/QBQl2iQcqhCIabsuJ1hpItVZRnJLPYm1zSTU2xhw0Wj6jZmOSq0+RxF0TZL4CCGEECeh3tLE/DVLWd9URXFAGMWGRN828rC2px4bFRsp7hJSLHUk1ZmIr2yhW89MRl12EwaD4dx0XrSSxEcIIYT4Hy6Pm/3OOp5c8BZ7DXr2+8VRoY5D0XSDiG6t9dSKh0RvGSm2KpIbTCRUmEg0GBl1063ERCeewwjEsUjiI4QQossraahiwcYf2GRuoiggkv26JOyRA46oF6HUkuooJ9nUQGJFMwkmG/0uvZg+F02RKasOQhIfIYQQXYrT4+bHXWv5Pn8He3QGiv3iqVHHgK4XhB2qd3DKKtVcS1JtM/GlTSRlxDHq6tsIDgo5dwGIX0USHyGEEJ1auamW+b8sZ4OpnkK/CIp1SThVgRA6tE29OG8FqfYKUhsaSSwz4V9XxwU33kDu+GtRqVTnqPeivUniI4QQotPweL2sKdrGd7s2sRMVhX7xvjNz6AYhh9bm+ClW0tz7SW2pJaWqmdjSBmK7xTHq6pmEhoTjcrlYuHAh2b2HSdLTyUjiI4QQosMy2S0s3PYzKyr3k28IplCfjFUVAAFt1+ccHM1Ja2gkscREmLmZrMkTGDL2OtQaWZvTlUjiI4QQosPY11DBV1tXstHcTIExilJNIl5VLATHttbRK3bS3PtJM9eQUt1CXFE9IZF6LrhuBqkpWeew9+J8IImPEEKI85LH62VdyU4W7dnMNo/CPkM8deoooDsEHqoXrtSR7iglramB5JJmwirriBmYwcjLbiUkOOyY7YuuSRIfIYQQ5wWb28n3ezewrHgPO9V+7Ds4baXv01pHpXhI8paRbqsgtbaZxEIT/rYG0ieMYNhVv0Wv15/DCERHIImPEEKIc6LBZmbh7jX8VFHGbl0QxdpkXKoA8Dt0eadesZPuLibdXENqZQtxBfVoA23kTLmSAZPHyvocccok8RFCCHFWVJgb+HrnalbV17JbF0GpJhFFFQX+Ua11ghQTmc5i0k31pJa2EFFYiy5RRf+pN9HrisGyw0r8apL4CCGEOCOKm2tYuGcDa5qa2KOPoVITBySC36GrHKKUGtLtJWQ0NJJSZCGwrBJ9dz+GXHYr3a7tI4mOaHeS+AghhGgXexvL+SpvHWubzewNjKVmWy2QCn6H6iR4S8mwlpNeZyJlnwVddTnGnGCGTrmDzJt6nauuiy5EEh8hhBCnTFEU9jSW89WuDaxtsbBHF0e9OgJIB6OvjkrxkHwg0cmoaSZ5jxlMFRh6hTB0yh10u6WXjOiIs04SHyGEECd0MNFZsHsja1vM7NbG06gOB9LA4KujUdykePeTaa4kvaaFpF0WPC1lGHr6M2jyTHrO7C+JjjjnJPERQghxVPlNlSzYs4E1TS3s0sbSoI4AUuHAjnGN4iLNU0yGpYrMSjOJuyy4zCWou2lRYnK46YXZGIyGcxmCEEeQxEcIIQQA+1vqWbB3PSvrG8jTxFCnjgSS/yfRKaKbuYqMKjOJu2w4m0tRJTvpMf4aht4+Bb1O23rPlWw1F+cjSXyEEKKLqrY28/W+DSyvqma7OoJqdQwQB7o4wDd1leotJtNcQbeqFhJ3OXE3luOKbSRp5DhG3vRbgoOCzm0QQpwiSXyEEKKLMDntLCz6hWXlJWzxBlGmSQAiQBsB+BYjp3hL6GYto1tVCym73Lhrq3CGVxAysDcjf38XifGJx38TIc5zkvgIIUQnZfd4WF62g4XFe9jkMlCkTsSr8gdVFmh8dRK9JXS3lpJZ20zaHjfq8kZs/oVoe4bT97Y76J09FLVaFiSLzkMSHyGE6CS8isLG2mK+2reVtRYPe9WJOFUGoHtrohOtVNHdVky3uiYyC9wYiq1Y1AV4kx3Ejb2CESMeJtDf77jvI0RHJomPEEJ0YIUt9SzYt4mf6k3sVMXSogoCUlsTnWCliR7OAro31NOt2EnIXg9m536ckZXo+vVh+My7SUtNP5chCHFWSeIjhBAdSKPTwaL9m1lSXsZmTwjV6iggGtTRABgUG93d+XQ31dCjzEbUbjU2UzWOgEI8GXribpvO0KGPY9Trzm0gQpwjkvgIIcR5zO1VWF1dwFfFeay1qtmnSkBRGYFMUINa8ZDmLaSHuZwelWaS9qjxVpuxqotwxlbhGTGMi0bNJjU59VyHIsR5QRIfIYQ4zxSZG/m6ZBvL65vZQSw2lT+QAgeOxYlVyulhKyartomMAgVtqYLNXo4zqJDGDCNZl93MBcOewl8ODxTiCJL4CCHEOWZ2uVhcvpNvSwvZGBBEzZYKIBJUkQAEKi30cO0lq7GWrGInQUUGnI31OLT7aY4tw2/UIIaNmEVWt55yJYQQJyCJjxBCnGWKorC1qYr5RdtY0eRgD/F4VFogHTS+gwPTvfvo0VJGr0or8Xu1uOtduGyV2IMLqE72kDz5KiZd+AeiwkPPdThCdCiS+AghxFnQ4HSwsHQ7iyrK2OgKw6QKAWLgwABNlFJNlj2fnnXNdC/woqkw4DLV41KXURe5FwamkjvgNgYNHodRL7+6hThdp3yRyooVK5gyZQrx8fGoVCrmz5/fpry6uppbbrmF+Ph4/P39mThxIvn5+cdt0+Vy8eyzz5KRkYHRaCQ3N5fvvvvuiHrz5s0jLS0No9HIgAED+Pnnn0+1+0IIcVZ4FYUN9RU8vmkxF/64kOyVO3m4RMsydyomVQgGxU5v91aurl/M01u+5skvtnDpx01kzq/BXLiG6rCPcU2yMOSpR7n3tXX85v7/cuGIiZL0CPErnfLfIIvFQm5uLrfeeitXXXVVmzJFUZg6dSo6nY4FCxYQHBzM3LlzGTduHHl5eQQEBBy1zSeeeIL333+ft956i6ysLBYvXswVV1zB6tWr6devHwCffPIJDzzwAPPmzWP48OG8+eabTJo0iby8PJKTk08jdCGEaF/1TgeLynbybUUZG51hB87UifEVqiDBW0qWbR+9as1k7lXjrdNgb7GgcVZTG7oDe7oXdcgALrtyHhkpKec0FiE6q1NOfCZNmsSkSZOOWpafn8/atWvZsWMH2dnZgG+UJjo6mo8++ojbb7/9qM+99957PP7440yePBmAe+65h8WLF/PXv/6V999/H4C5c+cyc+bM1jZeeeUVFi9ezOuvv84LL7xwqmEIIcSvpigKm5tqWbB/Gz82OslXYlFUaiAZVGBUrPR076KHqZreJR7CinW0mD2oms1YVaW0ROXh6htBZs9rmHjhi4QFBbBw4UKS4+PPdWhCdFrtOmbqcDgAMBqNra9pNBr0ej0rV648ZuLjcDjaPAPg5+fHypUrAXA6nWzatIlHH320TZ0JEyawevXq4/bnYJ8AmpubAd/UmsvlOoXIzh8H+91R+39QZ4ijM8QAEsepana5+b46n4WVJayzB9OkCgZ8u69QQaKynyzrPnrVWuhRoMNbq9Bk8+LX0kydYS+WuHyUnHR69pzGlGFvEBJw6Hef/FmcPzpDDND54mgP7Zr4ZGVlkZKSwmOPPcabb75JQEAAc+fOpaqqisrKymM+d/HFFzN37lxGjhxJRkYGy5YtY8GCBXg8HgDq6urweDzExMS0eS4mJoaqqqpjtvvCCy/wzDPPHPH6jz/+iL+//2lGeX5YunTpue5Cu+gMcXSGGEDiOBYFqFbDDp2VbRo/CjSJB3ZgJYLKd1JyT88uepgq6FMGEYV6zBYnNpuCp6WS+qA8bPFlNPTOIDhwNAkxt2LQqnC6YNVPP5yVGM6VzhBHZ4gBOn4cVqu13dpq18RHp9Px+eefM3PmTMLDw9FoNIwbN+6YU2MH/e1vf+OOO+4gKysLlUpFRkYGt956K//617/a1Pvf8ykURTnumRWPPfYYs2fPbv2+ubmZpKQkxowZQ0RExGlEeO65XC6WLl3K+PHj0ek67pHznSGOzhADSBxH4/R6WVVfwdfle/nZrKOK8DblsUoFPe176FXbTM8iP7TVXhrsoFi9uCyFNIbuwJFSjy0hl5zshxgwaCT+hhP3Sf4szh+dIQboPHHU19e3W1vtvj1gwIABbNmyBZPJhNPpJCoqiiFDhjBw4MBjPhMVFcX8+fOx2+3U19cTHx/Po48+SlpaGgCRkZFoNJojRndqamqOGAU6nMFgwGA48uRSnU7XoT8A0DligM4RR2eIASSOWoeLbyt2sbCinA2OMGwqIwcXJmsUFz28e+jesp8+FR6SioJxmlqodWowW+0EWotpjtyBq1sLtth+5PZ5hkFDRuFnOL1fsV39z+J80hligI4fR3v2/YztiwwJCQF8C543btzIc889d8JnjEYjCQkJuFwuPv/8c6699loA9Ho9AwYMYOnSpVxxxRWt9ZcuXcrll19+ZgIQQnRqiqKQ19LMV2V5LK5tYY8n8sDC5DhQQYjSSC9XHlkNNeTu9yOo0kCz2UWDW4fJ3IzRno8leifNvWxYI/vSt8/zDB46+qRGdoQQ584pJz5ms5mCgoLW74uKitiyZQvh4eEkJyfz6aefEhUVRXJyMtu3b+f+++9n6tSpTJgwofWZ6dOnk5CQ0Loba926dZSXl9O3b1/Ky8t5+umn8Xq9PPLII63PzJ49m5tvvpmBAwdywQUX8M9//pOSkhLuvvvuXxO/EKILcXkVfq6r5KvyPfzQpKKGUMDP96WCZKWInta9ZFeb6bE/El29kzqrh1qPgmI2Y7DtxRKXR0tvO0Hhvenb6zkGDx9PkJ/+HEcmhDhZp5z4bNy4kTFjxrR+f3ANzYwZM3j33XeprKxk9uzZVFdXExcXx/Tp03nyySfbtFFSUoJafejsRLvdzhNPPEFhYSGBgYFMnjyZ9957j9DQ0NY61113HfX19Tz77LNUVlaSk5PDwoULSZGzLoQQx2FyuVlUWcjXlSWssQZgxQ8IA0CnOOnpzaN7czF9KiCpNBaPyUO13UOJ10l0i51gWwGVcXkU97ETGNyd3j3+wJBRlxEeaDz+GwshzkunnPiMHj0aRVGOWT5r1ixmzZp13DaWL1/e5vtRo0aRl5d3wve+9957uffee0+qn0KIrqvEamd++W4W1dSxzRmGBw0Ht5wHK030dm2nW1MlfUoCiayOwm5SUem0UYCV6GYv0ZZ97I7bRX6ujaDAJHqkzmLg6GtIjAw5t4EJIX41OftcCNHhKYrC1uZmvja6+dPqlRR5ww6U+JKdBKWEbPsuutU2klMSQ1BDIM0tGqqcTdSr/Ilu0ZFs2sWu2EXs6m0lJCiC1OibyBl9M1nJsXLjuRCdiCQ+QogOyeVVWFVfx/zyPXzfpFCnBIE+A7ygUjz0UHaTZcmnR6WV7uUZ+JlCaTQ3Ue4sR9EkE9USREZjIQWRS9jW00p4qIHowIkMHz6TgTk90Kgl2RGiM5LERwjRYVjcHpbWVPJVZSE/NeuxYAQCATAodnK82+hmKqZnuUJqdQ8MLXHUmPdR6NqDok0jwhJDj9otlIX8wKbuNqKjPURoBjO6/+0MGzqCgNPcfi6E6Djkb7kQ4rxW73TzbVUJX1eWss7qjxMdEAxAsGIi17OFjIZSskr9SKrvgdacTJW5gD2uzaDLJNSWSlbNFhr1/+GXHhb8BzmIUdIZlnYzgy66gsSIoHMboBDirJLERwhx3imzO1lQUcQ31VVstQfhRQ34FhZHK1XkOreSUVdFZlkYCaZuaFoMVFr2stO1FkWXSZCrF91qtoLjY9b2aGR7PycZqiD6hFxDj9G30r97iqzbEaKLksRHCHFe2GO2saBiH9/WNrLHeXAUxpfspCiF9LbvIL22ltSyBBLNGajNYZSb97DNuRx03TDSn7SGPIKavmVzajkbR3jJCnCT7h1K5IDbGTlsOEFGOVxQiK5OEh8hxDmhKArbWqx8WV7AoroW9rsDD5QEoVK8dGcXOdZdpFY1kFiRQZK1F2pLE2WWPWyx70XRZaLTDiLRUkRUzSqKw/eyqI9CdKyNFFsiFyZMY8C46+geH37cfgghuhZJfIQQZ41XUdhgsjC/vIBF9VaqPP4HSgLRKC6y2U4v8x5SKk3EVmaR4uiH2tJMmWU3m6xbUHTpaHT9iHA3EF++Hrv3E1b0dtM40sUgL+Tah5E44D7GDu2PUac5p7EKIc5PkvgIIc4ot1dhTVMLX5YXsLjBSb334InH/hgUO72VrfRs2UtihYXo6mxSnUPRWs2+ZMeyHo82CbW+J0Ga3sRVrSek8TU2ZppYPEEhO8BGjiUNVcxNDBh9Bbs3/syEoX3RSdIjhDgGSXyEEO3O6fXyc0ML8yv2saTRhclrANSAET/FQl/lF7KaC4grdxBZk0Oa+0KMVjtllt38YlmNUx2FRt8Tvf8gEmu3Elf5CWWhxXyTq8ad4mSkDcapxhM5+HYuGpSLUafB5XKx+1wHLoQ470niI4RoFw6vl5/qm/miopDvG92YFT2+ZMdAoNJMf2UTPZr2EV3hIqyuN+mu0QTZFcosu9lm/gyLYkCj74km4AZiTKXElazCYNnGihw3749UGKazMqYlEWvgjfS94jqyEyPOdchCiA5IEh8hxGmze7z82GDii/IiljV5sCo6fMmOnhClkf7ejfQwFRJR7iGkLod090VEOvWUW/aw2/wNjS6zL9nxv4wQl4q4yrXEVn1JflwjHw5T4Ux2MdniJMk+Al2fOxl/4QWE+MnOLCHE6ZPERwhxSmweL8vqTcyvKGZZkwebosWX7KgJU+oZ6F1P96YigsshuL4X6e6xxLmDqbYUUGRexnp7OSpdBhrDMIzGBKLrtxOf/xEa+x6W58K8i2GwzsoVpiAq3TcSO/kWLuyZjFqukBBCtANJfIQQJ2Q9kOx8WbGfH5vcB5IdFaAlQqltTXaCKhT86rLJ8IwlyR1Bo62UYvNaNlsL8Koj0Riy0YdcQoC9ifj9q4mtfoOCWDP/GammKs3LVRYLt7f0oD7zVrKuu4zJ0XKqshCifUniI4Q4KmvryE4pPzQ5DyQ7AFoilRoGeteR1VSEX4UKfUMW6a6LSPNEY3c0UGzeyWLzLpyAWtcTTdA0DKoQomo2k1A5D71tHz/nqPjbRDWpAS6uNLmorh2NdsidjLlwsExnCSHOGEl8hBCtrB4vP9Q382VlGT80OrApB7eF+5KdQcpashoL8atUQ303MlxjSPfGoHV7KGrZznLzYlpcDah1aWj8xqHXpRFgrSW+aDlxVeupDrHw2RA1G3upmOi08pjJyA6uxTvhdqblpqLVqM9p/EKIzk8SHyG6OJvHy2atH1/t3MMPTYcnOxoilRoGK2vIMhXgV6nBVZdBmmsUGZ5YQrxGSi272dSylhp7CSpNCGpdNgb/bNQqP6Jqt5BQ8RrBpnx+yVTx7jUqGuPhxhYTk6sSyYu5B7+bbuCetMhzGr8QomuRxEeILsju8bK8oYUvqypZWm/B6h8NjW4OT3Z6NufjX6nGXJ9KsmMkmZ4YYpRQ6hzlFDSvoNSyBw8e1LpMtIFXodYmYXQ2kVD8A3FVa3Cpm1mWq2JJfw0JRiczTGYcpX0o7XE7w2+eyNiowBP2Uwgh2pskPkJ0EQ6vl58aWphfVc3i+hYs3kMjOxFKLUNYQ6/m3fhXqTDVpRLtGEamJ5ZEbwQubOxr3MJG8w4s7iY0+ghUxgvQ63uhUvsR1ribxPJ/ElG/g4pwL/8Zo2ZltpoLXTb+ZLKzt2Ekxf3u5Mqxw7g8yHjcfgohxJkkiY8QnZjT6+XnRjMLqmtZWNuE+bBkJ1ypYwhryDHnEVClUF+bQpBjCJmeWNI80WhVasqte1nRtIwa+35Uai1afTf0xgmotAnocBFbvpqEihUEWKvZmqrinxep2Juq4nKLhXeqPaxSJpE37C6uGtFHbkYXQpwXJPERopNxexVWNZmZX13HwtpGTJ6DC4Y1hCn1DGYtfazbCazy0FCbitYxgCRXDKM9sQRipEXVyLa6ZZRYduHyOtD7RaD1G4XmwOiOv7OBhIJPiatah4KdldkqvhmkoSXCy7TmZp4o1bJEM5Xto+9kxtAeclmoEOK8IomPEJ2AR1FY02Tmq+pGvqmpo6E12VETojQymLX0tW8mpMpFfW0qbkcuEc5ohnjiiFACceGi1LablfUbMLnqUKnUGIN7gCcblTYJrUpFREs+iUVLCG/YhdUIXw2BRQM0+Bk93GJqZOB+f74y3MyOiTO5e0Aaeq3s0BJCnH8k8RGig/IqChtMFhbUNPJ1dR217oMnG6sJUkwMYi0DnJsIq7ZTX5OK1ZpDtBLNSGcMCd5wVKhp1jWwpvwryi35ePGg9w/Bz38EXnqhqAPRa7zEVa0hYf/3BFirqQtR8e+xKn7IVRGrcvM7UxNZ1cF84Xc7AZfdyoN9k2VLuhDivCaJjxAdiKIobG62sqCmiQXVdVS5lAMlKgKUFgaynkHuDUTWtNBQk0Jzc0+CdJH0t0aT5olGjxaPwUuRdSd5FSuxupsB0AYmoNcMQFGloag0+KkdJBR+TXzZCnRuKyXRKt4eq2ZNTxUZbhdzmkwkmSP4zP8elKkzeKRvEhq5UkII0QFI4iPEeU5RFLabba3JTpnD21rmp1gYyHoGe9cRU9tEY00STU3d0ehD6OaIIdMRTZDdD9RgDmphY8mPlDXtRkFBqzcQljAES0sP1BrfWTqh6mYS8r4kunoTasXD3kQ1n1+gZnOGihynk1dqm4m3hPOh331Yp87gcUl4hBAdjCQ+QpyHFEVhl8V+INmpp9jubi0zKDb6s5GhyhoS6qtpqkmmoSGTcrU/Gdp4RtgjiLaFoEKFEqiizLuPTXsWYXebAQgIi8EY1B9zczo2qwG1BkKcJaTtmk9Y4x5UwNYMDZ9doGFPkoredgfzqk3EWcJ4z3AP2ZdO5/EBKehkSksI0QFJ4iPEeWSPxc6Cmka+qm6gwOZqfV2vOOjLJoYqq0k3ldNYE099XQr7PD1ICYqjjzeSZGs4WjSgVuGMdrOjYgX529e3thGVmoNX6Y3ZFIulRYVaA4nKfuLWv0egpRKAjT20fHqBQlHcgYSnykSKJYi3tbeTcfFt/GFIKgat7NISQnRckvgIcY7ts/pGdr6qaWS3xdH6uk5xkstmhrKKHtYSmiqjqa1NZbczk4igMAYZ40mtCyPAZgBAHaKjPqCG9dsX0LSv2teGwUhsN990VovJd1KyTq8ixbWHmFX/xuhoQlHBmp5aPhuuUBoFPR0u/lHVRHerkbe4iYiRM/n9yO746+XXhRCi45PfZEKcA8U2B18dSHZ2mO2tr2sUF33YylBW0cdVQGN5JLW1aey0pWM0GOgZmkpaXTgRtX6oUIEa1ClGCs3b2fjLV3jcvimxoIgoYjJGUFuRQm25b0rK6Kcm1b6dqB//hc5tQ1GrWN/byEdDXZRHQqbTxcvVJgZaVPzTeyVrBt/Bby7KISxAf05+RkIIcSZI4iPEWVJic/B1rYkF1Y1sM9taX1crHnLYxlBW0d+7i5aqCKqrUthqvgiVSk1abDLdAmKJK/dDa/JNM6mD9biSFbYULqXgh3WtbcVm9CAscTjlBeGU5fsWQQcGa0izbiF88dtovC4UtYpNuQG8N8RORYSbRJebF2pMXGR28R/PRH7odRf3ThpAYpj/2f0BCSHEWSCJjxBnUKndyTc1TXxV08TmFmvr6yrFQzY7GMJqBrEVV0MUZWXxbGkaC6gJDwtnRHo6qRUhGIsOtafPCKExuI51Gz6hdmuxry2VmvQBQwkIG0LRDi1N2zyAl5BwHRm2Xwj+9v9Qe3wjQdtzQ3h7iJmKCAcRHi+P15m4otnCAs9Ifhc/k7svG83diSFn7wckhBBnmSQ+QrSzsoPJTm0TvzS3TXZ6kncg2fkFnTWa4qIotjaMR1E06PV6+vboQTdHDMEFCqrKA8/5aTHmhlNsz2Pj8rcwNzYAoDP60XPEWLTGfuzdaMe9zwN4CI8xkOncSsBXb6By+xZIF+RG8vrARkqjLQR6YVZDEzc2t7DJnc19ATNJi43k5WkT0etlWksI0blJ4iNEOzg4svNNbROb2iQ7XrJak50NBHmiKS6OYWfVWDweX5KRmpJKr7B04kv9YKsd8B1KqEsIRNsniJ37f2br/IU4bb7pscCwcHpfdAluT092rWnA7bQAEBnvT3dlJ8Yv/wF2X93K7BheG1xPQWwTWkXFNFMLdzWZqHfH8bD6XgZOuI5XByWwdPF3qFRyHo8QovOTxEeI01Rid7JEH8zrWwrZctiaHRVespQ8BrOGwawlXBVJRWUSu4suwuXyAyAkJIQ+WTl0c8Sg3WbBu8cN2EGrwr9PFJ4MDb9s+Ja8N3/Ee2CaKiIxmb4TLsfSksr2n6twO+oAiEoOJMtYhOGTl/GaTAA0d4tj3lATvyTWAzDeYuX+hiaCXUZeds9ANfAW5kzoRXiAHpfLhRBCdBWS+AhxCoqsDr6pbeLr2ia2tdjAGAZmW5tkZxBridaG0tjUnd27LsJmCwBAq9XSu3dPsqMziSjU4ljRAIoJL6AJNRAwNA5zhJmfF39GwRdrQfGN/CRkZdP34qk01kSxfmE5Lkc54Et4sqPr0H/0R9zl5XgBZ3Is/xnpYUliDahU5DicPFLfQG+7m/c841mRcDuPTB1Kz7jgc/QTFEKIc0sSHyGOQ1EU9ljtfFtjYmFdEzsP23quwktPZSdDWMNA1hGt88fh6MPu3RPYUWdorRcfH0/f3n1Jd0fh3tCAe0MzB0/rMWSGEnBBHLWuUlZ+/TolO7a1Ppc5aCh9J06ltiSAnz8txWkrASAyKZDc7m70H72IIy8PN0BkOIvGhfJuyn4UtYpoj4oH6uu4xGJltSebmwy3c8MVF/NObrxMaQkhujRJfIT4Hwfvxvq21sS3tU0UWA8dKqjGQ7ayncGsYSDrCfQqGAzD2V98Gcv2eQFfUuHv709ubi590nrit9eNZXEVdnsZACq9Gv/+MQQMjWV/6XaWfPA8Vfvyfe1rNPQcMZr+k6+gokDNsn/vx9ZSA0B4fAD9hwTgt+ANLO8twwGoAgPYMbEbL6btxK5txoia2xobmWFqockbxiz3ncQPu47/G9edQIP8dRdCCPlNKATgURQ2mCwsrPWN7JTZD6170eIiR9nKYNbSnw2EabX4+Q2nsuJ61my14fUCKKhUarp160bfvn1J0UZjW1ONffl+zAcuUNdGGAm4IB6/fpHkb17D+r+8TF3pfl+Z3kDvsRPoP2kqFQUeFr1ZhLnRl3CFRPkxcGwMIas+pumhj7C43aDR0DBxIM/3KqBUuwOA8VY7D9fXE+WC//Ncyqq4W3niqkFkxcq0lhBCHCSJj+iy7B4vPze28F2dicV1zdS5DrsIFAd9lM0MYg39+IUgjZqQkFE0Nc1ky2YXjY0tgG9Bc2RkJP369aN3zxzU+VbM31XQWLPrUFvdwwgcFo8uPYhdq35k/R8+panKt1dd7+dP34svof+ky6gqcrHw9UIaq3y7wgLDDAyclExMyc80/P5Bmg4sXGbYQN640MoP6k0ApHvUPFZTyVC7gzWeXtynvYPrLx3HfwYkoZab04UQog1JfESX0uRys6y+me/qmvmhoRmLx9taFoCVfsp6BrGO3mzBT60mInwMDsdv2bVLS2FhCeA7Q0ev15OdnY3FYuGKcVNwbKqj5bU8FJsveVLpNfgPiCZwWDyqUB07ly9l/d8/o7nWN21lDApmwKTL6DvxUurLXSx6cx/VRc0AGAK0DJyUSpqxgvqXfkPt3r0A6DIzWHFlBn/TLserePFHw331ddzQ3EKLEshs160ofa7n/y7tRUSgASGEEEeSxEd0eiU2B4vrmllcZ2KtyYxbOVQWQRP9lDUMYj1Z5KFXa4iIGI1W+xiF+wJYu2YvDkdFa/3U1FT69etHz549cZWYKfxiM/WvbIED+ZMm3EjgsHgCBsbgVXvZtmwxG776DHODb1u5f0goA6dcSe74SZgbvCz79z6Kt/m2pWv1avqOSyYnR0fj3/5KxeLFAKhDgqm/aQJPRa6iyvkDKDDB4eWR6nJiPB7+6x7FuwG38ci0YYzuEX1WfqZCCNFRSeIjOh2PovBLs5WldSaW1Dez22JvU56kqqK/dxUDWEcahWjVOsLDRxIcfB1lpRGs+GkPdXWFrfVDQkLo27cvffv2JTQoBOv2Ohr/uRNXmZkwfIcQGjJCCByegDErHLfbyZZlC1m/4DMsB05ZDgyPYNBlV9F77MU4bbDqsyJ2rapAUUClVpE9Ip4BExKwf/4hpU+9gWKzgVqN4erLmDewgcWN88EJSSoDj1eWMtxmp9gbw/XuO+h5wWQ+ndCDAFm8LIQQJyS/KUWnYHK5Wd7Ywvf1zfxQ30L9Yet11ChkqfbRz7uSAawnRqlGrdYTHj6SyMh7aKhPZNu2vRQU7EE5cHaOVqulZ8+e9OvXj9TUVBSrG8u6KirX7sHbcmDhs1ZFXbiN7tcMxj8pBJfTwebvvvIlPE2NgO+W9MFTryFnzHgURc2WpSX8sqQEt8MDQHrfKIZOTUdXsIXqmx7AWVwMgN+AAWy+eSAvNnyMpdGCVqXm1hYHd9WVoPWqeNNzCfNDZzDnmsEMSAk7ez9oIYTo4CTxER2SoijssthZVt/MsvpmNjRb8Bw2hRWoctFHtZW+npXksplAxYxabSQiYjTRURNxurLYvm0vixZux27f2vpcUlISffv2JTs7G6PRiKvKQtMXBVi31HBwjkwdrCdwaByG/pFs+mkpPcK1/LLoa9Yv+LR1hCc4KpohU68le/RY1Gote9dXsWZ+IZYm306tmLRghl+VSWSQk+oXnqDlu+8A0ERGovrNDJ4I+pGt1W8DkKsO4KmSArq5XOzyJvGo+04GDx/HlxN6YNRpzsaPWwghOo1TTnxWrFjBn//8ZzZt2kRlZSVffvklU6dObS2vrq7m97//PUuWLKGpqYmRI0fy2muv0a1bt+O2+8orr/D6669TUlJCZGQkV199NS+88AJGoxGAp59+mmeeeabNMzExMVRVVZ1qCKKDanK5WdFo5seGZn6sb6HK2faqhRSNid7etfTxrqaHsgstHjSaACIjxhAVfTF63QB27sxnxU9bqKvb1PpcUFAQubm59O3bl8jISBSvgn1vI7Ur83EUNLXW0yUGEjQ8Ab/ekai0amxWK6a9efx70RdYGn1reIIioxh65XVkjxqLRqujoqCJlf/Np7akxVceYeSCKzLI6BtB08efUPjyy3gtFtBoCL3xBhZfFMZr+a/jrHMSoDFwf5OZ6+p24VXU/M1zBd+E3MgL1wxgYGr4mf+BCyFEJ3TKiY/FYiE3N5dbb72Vq666qk2ZoihMnToVnU7HggULCA4OZu7cuYwbN468vDwCAgKO2uYHH3zAo48+yjvvvMOwYcPYu3cvt9xyCwAvv/xya73s7Gy+//771u81Gvl/u52Zy6uwqdnCTw0t/NTYwpZmK97Dyo0qL320pfRy/Uiuso5ot2/HlE4XRmTklURHXUxg4ED27i1i8XdbKCx8vfXZg1NZffv2JS0tDbVajdfhwbymAvOqCtx1B+7eUoFfTiSBIxLQJwehUqnwuN3sXLaENZ9/jLnetzA5KCKKIVdcS86YcWi0Oloa7Kz+Yg8FGw/0yahh4KRU+lyUiDt/D/uv/w32Hb7zd4x9+uB+aCaza//Fjj2+10Zow3iqcAexHg/53gRmu+6h/9AxfDWpJ356+dwLIcTpOuXEZ9KkSUyaNOmoZfn5+axdu5YdO3aQnZ0NwLx584iOjuajjz7i9ttvP+pza9asYfjw4UybNg3w7Zy54YYbWL9+fdvOarXExsaeapdFB+E9MH31c0MLPzeaWWsyt9luDpCqs5LLNno6l9BDyUN/YNTHaEwkKupWoiLHExTUj5KSMn7+eSu7dr2G0+lsfT45OZm+ffvSq1ev1tFEd5OD5tUVWNZXodgPbEc3aggYFEvgsHi0Yb56Xq+HvBXLWfP5R5iqfSONGj9/Lrz2JnLHT0Kr0+Fyetj0dSGbl5TgdnlBBb1GxDNkSjpGnYfaV+bS8K93wetFHRhI5IMP8GW2hXnbH8XldRGk9eMRk43La7aiKCre8FzKe8Ybef7GgbJjSwgh2kG7rvFxOHzrFw7+gwK+URm9Xs/KlSuPmfiMGDGC999/n/Xr1zN48GAKCwtZuHAhM2bMaFMvPz+f+Ph4DAYDQ4YM4Y9//CPp6enH7c/BPgE0N/vOSXG5XB32RuqD/e6o/T/I5XLhBXaYzGy0OFhrsrK22UK9y9OmXpjGS19dGT1dK8hy/0SEs6G1LDAwh4iIi4gIH0NAQBa1tbWsX7+DHTteo6WlpbVeaGgovXv3pnfv3oSFHVoIbClswLq6CkdefZvt6P4XxGLsF4XaoEEBnA4HBRvWsPbzj2is8F0Q6hccQr9LplLlUdNz9Di8isKe9RWs+bIIy4ETl+MyQ7jgqnQiEwOxrltHxTPP4iot9fX94otx3ncT9+e/wtatvjVGI41x/L/8TcS43ZR6o3jQdQ/hPUfyxWVn9hb1zvSZOvy/HVFniAE6RxydIQbofHG0B5VycBvL6TysUrVZ4+NyuejWrRuDBw/mzTffJCAggLlz5/LYY48xYcIEFh84l+RoXnvtNR566CEURcHtdnPPPfcwb9681vJFixZhtVrp3r071dXVzJkzh927d7Nz504iIiKO2ubR1gUBfPjhh/j7+59u2OI0eYBStZ4CrYF8jZF8jQGLuu20jUHx0N1bRS+20Vv1EymqAtT4PqKKosXj6YbbnY3H3QtFCcXpdNLY2EhDQwN2+6Ft6xqNhtDQUMLDwwkICDh0MacCYfV6oiuNBJoP5f0twS6q4+yYwlwHr9tCURSsFaXUb9uI88AaHrXeQFivXEK690Kt1QHgalHTtMuAo97XnsboJSTLgV+sG43dRtS33xKyYaOvbkgI1VMvZ3mGme9s3+HEiQEds01ebmgoQgV87hnB854ZXJxmZEiUgtwpKoTo6qxWK9OmTcNkMhEc/Ouu4WnXxAdg06ZNzJw5k61bt6LRaBg3bhxqtRqAhQsXHrWd5cuXc/311zNnzhyGDBlCQUEB999/P3fccQdPPvnkUZ+xWCxkZGTwyCOPMHv27KPWOdqIT1JSEpWVlcdMls53LpeLpUuXMn78eHQ63bnuznE1uNxsbrHxS4uVTc2+/9q8bT9ufmroY7STrcon3fkDSc7VaDns6ghDHOFhowiPGE1oyBA0Gj+sVmtr0ltSUtJaV6PRkJmZSU5ODpmZmWi1hxIbr82NbVMN1rVVeE0Hpr40Kox9IvC/IA5dXNv1Z+W7d7L6k/epzN8NgM7oR//Jl9N34hQMB5JmS4uNBW+txrLfgOJV0GhV5I5Pou+4RLR6DZYVK6h55lk8Nb51PiHXXYfqrpt5ZsefWVW5CoBBgak8V7CFBFszzYo/j7tuIz9qAq9cl0tG1NHXxLW3jvSZOp7OEEdniAE6RxydIQboPHHU19cTFxfXLolPu29nHzBgAFu2bMFkMuF0OomKimLIkCEMHDjwmM88+eST3Hzzza1TYb1798ZisXDnnXfy+OOPtyZOhwsICKB3797k5+cfs12DwYDBcOTR/TqdrkN/AOD8i8Hi8bCzxcbWFhtbWqz80myhyOY8ol6oVk1ffxc9KCS6YSGZmrVoLYfqqVR6QkOHERExkojwUQQEdEOlUuFwONizZw87duygoKAAr/fQ2p/k5GRyc3Pp1asXfn5+bd7PVWfDvKoc66ZqFKfvGXWAjoAhsQReEI8mSN+mfnXRPlZ+/B+Kt/h2fWl1evpOvJTBl1+NX5DvL5uiKBRsqmHlp/lYTXpAIS03kuFXdyMkyg+PyUT1c3/CNH8+APrUVOKen8PGGAtP/nwbDfYG9Go9DxhTuHH7MtTAem8PHnDex/gLBjB/cs9zsk39fPtMna7OEEdniAE6RxydIQbo+HG0Z9/P2Dk+ISEhgG9dzsaNG3nuueeOWddqtR6R3Gg0GhRF4VgDUg6Hg127dnHhhRe2X6fFSal3uskz29hptrHjwNdeix3vUepm+unp7Wejm6qQVMcKQszLwHRgFO7Av+3+/pmEhw8jPHwEYaFD0Wp9Ix1Op5O8vDx27NhBfn4+bvehkaDY2Fh69+5NdnY2oaGhbd5TURQc+5owr6zAvqeBAzNl6GL9CRyegH/faFS6tp+3hopyVv33ffau+RkAtUZD74suZuiV1xEYfmh0sKnayoqP91C6y3dAocbfy/jpvcnoGwNAy/LlVP2/p3DX1IBKRfgttxD8m7t4decbvL/sfQC6B6XwYnUNmfuW4UXFq+7L+Zf2Ol68uT8TsmXxvhBCnEmnnPiYzWYKCgpavy8qKmLLli2Eh4eTnJzMp59+SlRUFMnJyWzfvp3777+fqVOnMmHChNZnpk+fTkJCAi+88AIAU6ZMYe7cufTr1691quvJJ5/ksssua92y/vDDDzNlyhSSk5Opqalhzpw5NDc3H7EAWrSfBpebfVYHeyx29lrs7Dnw9b/n5xwUq9fRO1BHd1096d5dJNh/wtuyHsXaduTHoI8hJGQwJSUBjBp1N4GBSa1lTqeTnTt3kpeXx969e9ssaAsPDycnJ4fevXsTFRV1xPsrLg/WzbW0rCrHXW1tfd2YFU7giHgMGaGH1vocYG6oZ81nH7H9xyUoXi+oVGQNG8nwa28iNDautZ7b5WHTov38smQ/XreCRqum7/hEqjw7Sc4Ox2O2UPPin2j69DPgwCjPH/9IVXoI9y67jb2NvotGb4oawv2bv8XoslKrhPCA616a40bw1Y39SQqXdWdCCHGmnXLis3HjRsaMGdP6/cH1NTNmzODdd9+lsrKS2bNnU11dTVxcHNOnTz9inU5JSUmbEZ4nnngClUrFE088QXl5OVFRUUyZMoXnn3++tU5ZWRk33HADdXV1REVFMXToUNauXUtKSsopBy18XF6FSoeTUruTMruLUruTYpuDQpuDQquDJrfnmM+m+unJDjCSabCSqhST4NqIzrwaa/2+1joHn9browgNHURY2AWEh12An18qbrebwsKFGAyx2O128vPz2bVrF/n5+W2SnZCQEHJycsjJySE2NvaIxAXAY3JgXluJZV0lXuvB29HV+A+IIXB4ArpIvyOesVvMbPjqc35Z+BVup28EKr3/IIZfdzPRqW13CpbtbmD5h3sw1fjO9knuFc6F13cnIEzHwoU7sf3yCzWPP4GrrMw3yjN9OlEPPsA35UuY8+1d2Nw2wg1hPGdIZ+T6TwH42ZPDg677uHhIb568tJecwCyEEGfJKSc+o0ePPub0E8CsWbOYNWvWcdtYvnx5205otTz11FM89dRTx3zm448/PqV+djVur4LV68Xi8WDxeGlxe2l0uX1fbg8NLje1TjfVDhdVThc1Djc1TtdRp6cOF2/Q0SPASDc/HanaZuLZT4xzCx7zZlrqd+H1+nZSuQ58Afj5pRIaMoDQ0EGEhg7Czy/liITFYrFQX1/Pxx9/THFxMR7PoSQrNDSUXr16kZ2dTXx8/FGTHUVRcJa0YF5dgW17HRxYNK0JNfhuRx8Ui9rvyI+32+lky+JvWPflf7FbzADEdc9i5LRbSOyZ06auzexk1WcF7FnrO7PHP0TPhdd2J6N/FCqVCqfFQuSiRZT/tAIUBV18PHEvvIB6QG+eWvc88wvmAzA0qh8vlJcQuftLAP7mvoJ/qq5lznV9uKJf4gn+BIQQQrQnuaurg7lyWxG7AxP5/drdeBQFt+K7jdx5mpvzDGoVCQY9SUY9iUYdKUYdCdoWYpQKotz5uG27MZt3Y2nYh6L4UpuGw57XaAIJDsohJKQfISH9CQ7ui15/5HUKiqJQV1fHnj172LNnD6UHzrM5KCIigqysLLKzs4mLiztqsgOguL1Yt9ViXl2Bq8zc+ro+LZig4QkYe0ag0hz5rNfrYdfPy1n13/dpqav1vWdiMiNumEHGgMFt3k9RFPasq2LVpwXYLb7t7TkjExg6NQPDgWTKUVhE+UOzCd/l2/UVcsUVxDz+B4rd1Tz87TQKmgpQq9TcmzyZ29d/gsZaj0kJ4AHXvewLHc7n0weQFfvrdiYIIYQ4dZL4dDCNLjctag0cYxpKo4IAjZpAjYYwnYYwrZYwnZYwnYZovY4Yg5ZIrZcwmgj1VuPv2o/DXoLVth+bqQRbVTFerxMHUPY/bWu1QQQG9iQoKJvgoN4EBfXG3z8VlerIXXfg20a5f/9+8vPzyc/Pp6GhoU25n58fgwcPJjs7m6ioqGMmOwCe5gPTWeur8JoP3Y7unxtN4LB49AmBR31OURSKt/7Czx/8i9qSYgACIyIZds00skeNRf0/5wg119lY/uEeSvN8fY1ICGT0jT2ITQ9pba/ps8+o/uMLKDYbHn9/Ev74PGETJ7KoaBFPrX4Km9tGpF8kL0WOYOBP/0CleNjhTeVu1wOkZfbiqxv6EerfdjeZEEKIs0MSnw7mje4xrFi1ghHDBqPXqNHgQY0bPXYMigW1x4LHa8HtbsblbMDpavD911aPs6kGu6MSt7sZC2A5xnuo1Qb8/TMICMggwD+TwMAsAgN7YjQefdrpIEVRqK+vp7CwkPz8fIqKitrsxNJoNKSlpdGjRw/S09NZuXIlI0aMOOY2RUVRcBY3Y15TgW1H/aHprGA9ARfEETAoFk3gsROI6sICVnzwL0p2+E5GNvgHMHjqNfSbNAWdvu0xB16vwvYfy1j7VSFuhweNVs2gS1PpOz4ZjcaX2Hmamqj8f0/RsmQJAH5DhrBj3FjSx4zkxfUv8v4u366tobGDecGiJnLFq4DvQMI/uG7nlpFZ/O7iHmg1R08UhRBCnHmS+HQwpt3X0CNgP7Vbf107Gk0gRmMcfn4p+Pkl4+eXhJ9fMgH+6RiNCahUJ7fY1mw2U1RURGFhIfv27Wu9FuSgoKAgunXrRmZmJhkZGa3nKh3v+HGvw4N1Sw2WtZW4Kg+lZ/rUYAIviMcvJwLVcZIHU00VKz9+j92rfvLFqtXS9+JLGXLFta1n8RyuodLCD//ZRXWRr+/x3UIZc1MWoTGHdllZf/mF8ocexl1ZCVot0Q8+QNBNN7Fy4afcvexuNtduBuCOHjdw345laMo24kHNH13TeF91CS9dn8vlfRNO9OMUQghxhkni08G0TUhUqFQ6VCoNGo0fWk0gGm0AGk0AWm0gel04On2477+6CAyGaAzGOIyGWLTaoNN6/+bmZvbv309xcTH79++nrq6uTblGoyEpKYmMjAy6detGTEzMcUeJDueqsWJZW4llUzWKwzeVp9Kp8e8bTcAFcejjjz6ddZDN3MK6Lz5hy+Jv8BwYaeo5YjTDr7uZkOiYI+p7PV42Ly1hwzfFeNxe9EYNF1yZSfaIeFRqX58Vr5f6t9+m9pW/gceDPiWF+L/8Bb/eOWys2Mi8lnm0NLcQoAvg+V63M3bZX6ClEhOB3Of8Lbv9B/LJjIH0TQo9qZ+BEEKIM0sSnw6mf7/PWbx4KZMmXYpef+Sp1O3J7XZTXV1NaWkpZWVllJaWYjKZjqgXExNDeno66enppKSkoNef/PoVxe3FmleLZV0ljsJDbWsjjAQMjSNgQAxq/+Of2Ol2udjy3ddtdmol5+Qy8sZbiUnPPOoz9eVmfvjPLmr2+y4zTcmJYPSNPQgMO3TBrruhgYpHH8WywneoYfCllxL79NNoAgOYXzCfZ9c8i0txkR6SzispV5D27WPgslKgJDDT+RB+Md1YcMsgEkKP3E4vhBDi3JDEp4PRaPwA7TEXFJ8uh8NBbW0tlZWVrV/V1dVtroYA3/1ssbGxpKSkkJqaSnJy8mld+OpusJOw34/av/yCYjmwDkgFxp4RBA6Nw5AZ2jrqciyKorBn9QpWfvwfTDXVAEQmpTDyxltJ7TvgqCNNXo+XX5aUsOGbIrweBYO/lhHXdqPHkLZnBFk3baJ89kO4q6tRGQzEPPE4oVdfjUfx8NKGl3gv7z0Aeml78n/hQwhc8ACgsMLTm/tc9zM4K5W/3dCPQIP8FRNCiPOJ/FbuQrxeL83Nza23mdfX11NbW0tNTc1RR3LAt/MqMTGRxMREkpKSSEhIOOr9ZydDcXmx7azDsr4KR6GJWPxQcKMO1hMwKJaAQbFoQ0+u7bK8Hfz0/ttU7fPd1RYQFs7wa28ie/SRO7UOaqyy8P27u6gp9q3lSe0TyehpPQg47D0VRaHxvfeofunP4HajT0sj4ZWXMfboQbOzmUd+eoRVFb4LRu/Mvp3Ja38kKP9ZAD5wj+Up9wxuGpbJk5f2QnOCxE0IIcTZJ4lPB+P1evF6vTgcDpxOJ16vF7fbjdPpxG63t35ZrVZaWlpav5qbmzGZTG0OCvxfgYGBxMTEEB8fT1xcHHFxcYSGHnnNw6lQFAVXhQXLxiqsW2pRbIdGd0whTpIn5xCYHX3Us3eOpqGijJ8/fJeCDWsB363pgy67koGXXIHOaDzqM4pXYduPZayZvw+Py4veT8uF1x05yuO1WKh88kmaFy4CIHjyZOKeexZ1QAClzaXcu+xeipuLMWqMzBn6JOPX/Bt13Y94UfG8axpveybz+ORe3DEy/aj9EEIIce5J4tPBvPHGGzQ2NrJ16+lt61Kr1YSGhhIeHk54eDhRUVFER0cTFRV1WlNWx+JpcWLdUoN1Uw2uqkM7szShBgIGxqDPjWDT6mV07xl+UkmPtdnEms8+ZOvSRSheLyq1mj5jL+aCq6cREBp2zOea62388O9dlO9tAiCpVzgX3ZzVZi0PgKOwkLJZs3AW7AOtlphHHiHs5ptQqVRsqdnCrB9m0ehoJDYglleHPkvPRY9D+Sbs6Pmt8zcsVw3ib7JzSwghznuS+HQwBy9tPZxWq8VgMGA0Glu//Pz8CAoKavMVGhpKcHBwm3vS2pPX6cG+qx7r5lrsextovQ9Do8KvVwQBg2Jb1+4cbzv74VxOB798u4D1Cz7FafPdlZU+YDAjp91KRGLSMZ87ePryzx/vxWn3oDVoGH5VJtkXHnkWUcsPP1Dxu0fwWixoo6NJeOVl/Pv3B+C7ou94fOXjOL1OsiOy+fuA3xP56Uxo2EezKogZ9ofZq8/i3zcPZFhm5Kn/0IQQQpxVkvh0MNOnT+f7779n0qRJGAyGM5bEnCzF48We34RtSw22vHoU56HF0PrkIPz7x+DfJ/KEO7OOaNfrZdfK5az8+D1a6n1XTESnZTDqppkk5/Q57rN2s4vlH+5m3y++52LTQxh3a09CotqOaCleL3VvvEHdq68B4D9wIAkvz0UbFYWiKLy9423+9svfABiTNIY/dbsR/w+uA3M1VapobrT/jhptPB/NHEyf5COv6RBCCHH+kcSng/Hz80Oj0aDVas9Z0qO4vdgLmrBtr/MlO7bDTmeOMOKfG4V/32h00ac3dbZ/2xZ++uAdaosLAQiKjGLE9dPpOXwUqhPEXJJXz7J/78JqcqJWqxh0aRr9L05G/T8HHnotFioe+0PrKcxhN95IzKO/R6XT4fF6mLNuDp/t/QyAm3rexMPRw9G8dxU4W8hXpXCj7RF0oXHcn2qmZ9zpnYkkhBDi7JPER5wUr92NfW8j9rx6bLsaWg8YBFAH6vDvE4Vf3yj0SUGnvRi6rqSYFR/8i6ItmwDQ+/kzeOo19J982RFXTPwvj8vLmvn72LrMd/lpWKw/427tRXTKkSc1O8vKKLv3Phx794JOR+z/e5Kwa64BwOFx8PsVv2dZyTJUqPj94N9zoy4G3r8a3DY20otbbbOJjoriXzMGsHnVD6cVqxBCiHNDEh9xVIqi4K6zYd/TiH13g+9wQe+hG+DVwXr8cyLxy4lEnxp8wjN3jqeloY7V//2Qncu/R1G8qDUacidMZuiV1+MfHHLC5xsqLSx5eyf1B25rzxmZwLCrM9Hpj1wPZf3lF8ru+w2exkY0kZEkvvoq/v37AdDsbGbWD7PYVL0JnVrHiyNfZLzVDh9eD14XK+jHHfb7yYyP5D+3DSbYoGbzaUcthBDiXJDER7TymJ04Cpqw5zfhKGjCY3K0KddG+WHsGY5frwj0yb8u2QHwOJ2s/u/7bPnuG9xO33t1HzKcEdNmEBYbf8LnFUVh588VrPo0H7fLizFQx0U3Z5GWG3XU+qavvqLy8SdQXC6MvXqROO8f6GJjAaix1nD393eT35hPoC6QVy96lUHVhTD/HlA8LFaG8hvHveSmRPHOrYMINupOeoG2EEKI84ckPl2U4vWN6DiLm3Hsb8a5vxl3na1tJY0KQ2owxqwI/HqGo41sn6sXPG4XWxZ/w/6vP6HIYQcgIasXI2+8lfjuPU+qDbvFxY/v7aZwi28Bc1LPMMbe0ouAkCOnxBSvl9rXXqP+9TcACBo/nvgX/4T6wPb90uZS7lh6B+XmciL9Inlj3Bv0KFoDXz8AKHypjOZhx0wGpEbxr1sHESCnMQshRIclv8G7AK/djbvOhqvSgrPCjKvSgqvS0madzkG6uAAM3UIxZoahTw1GfZTpotOleL3sXr2CVZ+813rFRFhcAiNvvJWMgUNOem1QZUETS97ZibnBgVqjYujUDPqOTTrqCJTX4aDyscdaDyWMuOMOoh58oHWRdGFTIXcsuYMaWw3JQcm8Mf4NkvZ8D988AMAHykSecNzE4LRI3rlFkh4hhOjo5Ld4B6O4vajd4DW7cKs8KC4vilvBa3HhNTvxmF14zS48JgfuehvuejteyzGmZLRq9EmBGFJD0KcEY0gOOuVt5yfVZ0Vh/9ZfWPHRv1t3avmHhBHQPZvrfvMAhmOcuPy/vF6FX74rZv03xSheheAoPy6+PfuoC5gBPE1NlN73G2ybNoFOR9zTTxN61ZWt5bsbdnPX0rtosDeQGZrJWxPeIjLvm9ak59/KJTzlmMbQ9AjeuWUQ/nr56yKEEB2d/CbvYOpf3Uq/xnBqN2w6pefUgTp0sQHo4gLQxQeijwtAG+WHSnNmt8RX7N3Fyo/fo3TnNuDATq3Lr6b3uEks/eEH1Ec5kPFoLCYHS9/Jo3xPIwDdB8cw6oYe6P2O/hF2lpVTeuedOAsLUQcFkfj3vxMwZHBr+bbabdz9/d20OFvoFdGLN8e9SejO+fD1/QC8p0zmKcc0hmVE8vaMQfi148iXEEKIc0cSn45Gq27zv1VaFSqtGrW/Dk2gDnWgDk2gHnWwHm2EEW2EH9oII+qzPEVTU1zIqk/eo/CXDQBotFr6XnwpQ664Fr+g4FNaGFy2u4El7+Rha3ai1asZdUMPegyNPebUmD0vj5K77sJTW4c2Npbkt/6JoVu31vKNVRu5b9l9WN1W+kb1Zd64eQRt/6I16fmAyTzpuJHBaRGS9AghRCcjiU8HE35XDouXLmbiJZPQ6/XnujtHqC8vZc1nH7Fn9QoAVGo12aPGccHV1xMcGX1KbXm9Chu/LWLDwmJQIDw+gIvvyCE8LuCYz5hXrqJ81iy8ViuGHj1I+ueb6GJiWss3VW/i3mX3YnPbGBI3hFfHvIp/3lfw9SwAPlJfwuPWaeQmhfHOLZL0CCFEZyOJTwejNmhQ1PyqG9PPhIaKctZ+/hG7V61AUXzXVvQYNpJh19xIePypX9z5v1NbPYfHceF13Y96Ns9BzQsXUv77R8Hlwv+CoSS++iqaoEOnKm+u2cw939+DzW1jWPww/jbmbxjzl8L8ewH4TDOZxyzTyIoN5t+3DiJQFjILIUSnI7/Zxa/SWFnO2i8+YdfPy1sTnoyBQxl2zTSiU9NPq83yvY0s+b+dWJudaA0aRk/rQY8hscfvx0cfUfXsc6AoBE+eRPyf/oTqsBGxLTVbuHvp3djcNobGDfUlPftXw2e3geJhkfYifmeeRnpkIO/NHEKo//k3miaEEOLXk8RHnJbakmLWfflf9q5Z2ZrwpA8YzLCrpxGTnnlabSpehc1LS1g7fx/KgamtiXfmEBZ77KktRVGof+MNav/2KgChN1xP7BNPoDps0fTW2q3c/f3dWN1WhsQO4dWLXsVYsRU+vhE8Tlbph/Ob5luJDw3g/duHEBV0/OsxhBBCdFyS+IhTUpm/h7VffkLhpvWtr6X3H8QFV91AbGb3027XbnGx7N08irfXA9BjSCyjpvVAZzj21Jbi9VLz4ks0/PvfAETeew+Rv/1tm2nAvPo87l56NxaXhUGxg3ht7Gv41RXAB9eAy8o240BubbqTkAA/3ps5mPjQ9jmkUQghxPlJEh9xQl6vh32b1rPpmy8p353ne1GlovvQEQyZes1pT2kdVFvSwqI3t9NSb0ejVXPhdd3oNSL+uOuYFK+XqqeepunTTwGI+cNjhE+f3qZOoamQu5fejdllpn90f/5+0d/xa6mB968Ch4lCv95c23gfGp2Rd24ZRHpU4K+KQwghxPlPEh9xTC67nR0/fc8vCxfQVFUJgFqjpeeI0QyeejXh8Ym/+j12r6li1X/34XF7CY40MvHO3kQlBx33GcXtpvLxxzEt+ArUauKef57QK6a2qVNhruDOJXfS6GikV0Qv/jH2H/i77PDB1WCupsYvk8sb78elNvJ/N/anb1Lor45FCCHE+U8SH3GEupJitn6/iLwVP+K0WQEwBgTSZ/wk+l18KYHhEb/6PdxODw3bDZSV5QOQ2ieScbf0xHCCk6MVl4vy3z1Cy3ffgUZDwp9fInjy5Lb9t9Vxx5I7qLZWkx6Szhvj3iBQpYGPboC6vZiNsUxpfIAW/Pnzlb0Zk3Vq2+yFEEJ0XJL4CMA3upO/YQ1bly6iYk9e6+uhsXH0n3w5OaPGoTvJqyVOxFRrY9Gb27CW6VGpYMjl6fSfkHLC2969TiflD87GvGwZ6HQkzP0rwePHt6nT7Gzm7qV3U9JSQnxAPG+Of5MwfTB8OgNK1+LSBXNV82yqCefhCd25ZmBSu8QkhBCiY5DEpwvzejyUbN9C3srlFKxfg+vATekqtZrMgUPJHT+Z5Jw+rRd6tof9O+pZ+s5OHFY3ar2XSXf2ITXnxCMuXqeTst/+FstPK1Dp9SS+9iqBo0a1qePwOPjtst+yp3EPEcYI3prwFrH+MbDo97Dra7xqPXc4Z7PHm8i1AxO5b8zp7T4TQgjRcUni08W4nA5Kd26jcNN68tevwWpqai0LiYkle+RYel80oV2msw6neBU2LCxmw7dFoEBUShCatEoSeoSd8Fmv00n5rPt9SY/RSNK8fxAwbFjbOoqXx1c+zi81vxCkC+LN8W+SHJwM696E9W8C8JRmFsut3RmSFs6cqb3Pu0MghRBCnHmS+HRyitdLfVkJZbvzKNqykZLtW3E7Ha3lxqBgsoZdSM8Ro4nrlnVGkgG7xcX3/8pj/w7fVvWckQkMmZrK4qUVJ+6/00n5Aw9iXr4clcFA0uvzCLjggiPqzd04l8XFi9Gqtbwy5hV6hPeAgmXw3aMA/DvgNt6r709qhD9v3DQAvfbMXs4qhBDi/CSJTyficjowVVXSWF1JbXERFXt3UVWwF4fV0qZeUEQU6f0HkTFgMMm9+6LRnrmPQV1ZC4ve2E5znR2NTs3oaT3IuiDupC4pVVwuyh96CPMPP/imt+b946hJzwe7PuDfeb6zfOYMn8PguMFQuxc+vRUUL2tDJvFU9ViCjVrevmUQYQFyKrMQQnRVkvh0ML8sXEDtpo0sqyzG7XDgsttwWK2Yaqsx19cd9RmdwUhsZndSevclvf8gIpNTz8o0z551VSx/fzdul5egCCOT7jrxVvWDFLeb8od/R8vS731Jzz/+QeDw4UfU+37/97y4/kUA7u9/P5ekXwLWBvjoOnCYqAzpy/Tq69Gq1bx+0wAy5KweIYTo0iTx6WB2/LAEU1UFpj07jlpu8A8gNDae8IRE4rtlEd+jJ5FJKag1Z++WcY/Hy6rPCtj+YxkAydnhjL8tG2PA8beqH6R4vVT84Q+0LF6MSqcj8e+vEXjhiCPqba/dzqM/P4qCwnU9rmNmzkzwuOC/06GhEHtAIlNq7saJjuem9GJ4ZmS7ximEEKLjkcSng+k58iJ2bd9GVq9sjAGB6P380Pv5ERQRRWhsHH5Bwed00a7F5GDxWzuoLDABMHByKoMuTUN9gq3qBymKQtVzz9H81de+c3r+9gqBI0ceUa/aUs39P96Pw+NgVOIoHh38qC/u7x6F4p/x6gK4yfogdUowVw9I5KahKe0apxBCiI5JEp8OZtBlV1Gr9WPw5MnodCc3gnK2VOQ3sfitHVibneiNGsbdlk1an1MbZamd+zJNH30MKhXxL75I0EUXHVHH5rZx/4/3U2urJTM0kxdHvohWrYUtH8GG/0NBxXOG2WxsiaNPYghzpubIDi4hhBCAJD6iHSiKwrYfylj9eQFer0J4fACT7upNaIz/KbVT9+Y/qX/rLQBin36akEsvOep7PbXqKXbW7yTUEMprF71GgC4AKrfBNw8AsDjqFv5V2pPwAD2v3zQAo+7sTfMJIYQ4v0niI34Vl8PDj+/vJn9DNQDdBsUw5qas496qfjQNH35I7csvAxD9u98Rdt21R6331va3WFS8CK1Ky8ujXyYxKBFsjfDfm8FtpyzyQu4pHYtGreLv0/qRILetCyGEOIwkPuK0NVZZ+O6fO2iosKBWqxh2VSZ9Lko85Wml5u++o/q5OQBE3HM3ETNvO2q9ZfuX8drm1wB4YugTDIwdCF4vfHEXNBbjDEriisoZKKh5dGIWwzJkMbMQQoi2JPERp6VgUw0//GcXLocH/2A9F9+RQ3y30FNux7phA5W/ewQUhdAbridq1qyj1isyFfH4qscBuLHnjVzV/Spfwc9/gfzFKFoj97oepNbjz4ReMdx+YdrphiaEEKITk8RHnBKP28vqLwrY9oNvq3p8t1Am3J5NQIjhlNvSV1RQ9dwcFJeLoPHjiX3iiaOOFlldVmYvn43FZWFgzEAeHviwr6BwOfz4RwA+iLyf74tjSQj1489X58piZiGEEEcliY84aS0Ndpb8306qCn1b1ftfnMyQy9JRa079+gdXeTmJ7/wLr9mM/8CBxP/lz6iOctaQoig8s+YZCpoKiPKL4s+j/uzbwWWuhS/uBBT2JV3FE/m5aNQqXr2hHyH+59duNyGEEOePU/4Xa8WKFUyZMoX4+HhUKhXz589vU15dXc0tt9xCfHw8/v7+TJw4kfz8/BO2+8orr9CjRw/8/PxISkriwQcfxG63t6kzb9480tLSMBqNDBgwgJ9//vlUuy9OU/G2Oj55fj1VhSb0flom39ObC67IPK2kx93YSMXd96BtaUGfmUnivH+gNhx9xOjjPR+zsGghGpWGv4z6C5F+kb51PfPvAXM1jrAeXFl0GQC/u7gHA1JOfOmpEEKIruuU/9WyWCzk5uby97///YgyRVGYOnUqhYWFLFiwgM2bN5OSksK4ceOwWCxHac3ngw8+4NFHH+Wpp55i165dvP3223zyySc89thjrXU++eQTHnjgAR5//HE2b97MhRdeyKRJkygpKTnVEMQp8Li9rPosn2/nbcNhcROVHMS1fxhIWm7UabXndTop++1vcRUX4woNJf6N19EEBx+17tbarby04SUAZg+YTf+Y/r6CtfOgYCmK1sh9zt9gcusY1T2KOy9MP60+CSGE6DpOeapr0qRJTJo06ahl+fn5rF27lh07dpCdnQ34Rmmio6P56KOPuP3224/63Jo1axg+fDjTpk0DIDU1lRtuuIH169e31pk7dy4zZ85sbeOVV15h8eLFvP7667zwwgunGoY4Cc11Npa8vZPqomYA+lyUyLArMtHoTu9mc0VRqHz8CWwbN6EOCqL8tlvpGRNz1LqN9kYeWv4Qbq+b8SnjubnXzb6Cis3w/dMAfBP7G74viCA6yMDca3NP+nRoIYQQXVe7rvFxOBwAGI3G1tc0Gg16vZ6VK1ceM/EZMWIE77//PuvXr2fw4MEUFhaycOFCZsyYAYDT6WTTpk08+uijbZ6bMGECq1evPm5/DvYJoLnZ9w+4y+U6qdvBz0cH+30m+68oCgUba1n53wJcdg96Py2jb+xGam4kXjx4XZ7Tard+3jyav/4atFoiX3qR3c3NR41DURSeXPkk1dZqUoJSeHLwk7jdbnC0oP30VlReFzUJ4/ltQT8AXrwyh2CD+qz/mZ6NP4uzQeI4f3SGGKBzxNEZYoDOF0d7aNfEJysri5SUFB577DHefPNNAgICmDt3LlVVVVRWVh7zueuvv57a2lpGjBiBoii43W7uueee1kSnrq4Oj8dDzP+MDsTExFBVVXXMdl944QWeeeaZI17/8ccf8fc/tVOFzzdLly49I+16XdC404it0rdAWB/qITzXTF75evLKT7/doE2biPvvpwBUXX45ew8koUeLY51jHT/ZfkKDhkuVS1mxdAUA/fa/RXJjERZdOFeWXA2oGBnrpXnvOhbuPf2+/Vpn6s/ibJM4zh+dIQboHHF0hhig48dhtVrbra12TXx0Oh2ff/45M2fOJDw8HI1Gw7hx4445NXbQ8uXLef7555k3bx5DhgyhoKCA+++/n7i4OJ588snWev+7RVlRlONuW37ssceYPXt26/fNzc0kJSUxZswYIiIiTjPKc8vlcrF06VLGjx/f7nd1Vext4sf39mBrcqJSw4BJKfQdn4Ra8+umkKwbNlDxxZcAhN0+k8z77z9mHAVNBTy3+DkAHuj/ADdm3QiAas9CtJt/RkHF61H/j7LCINIjA/j7HUPx05+bKynO5J/F2SRxnD86QwzQOeLoDDFA54mjvr6+3dpq9+3sAwYMYMuWLZhMJpxOJ1FRUQwZMoSBAwce85knn3ySm2++uXUqrHfv3lgsFu68804ef/xxIiMj0Wg0R4zu1NTUHDEKdDiDwYDhKLuFdDpdh/4AQPvG4LS5WTN/Hzt+8g3phET5Mf62bGLSjr7o+JTaLi2l+sHZ4HYTNGkiMbNno1IfWiN0eBx2t50/rP4DDo+DEQkjmJEzw5fYWuph0UMA5Gfext93RKNVq3jl+r4EBxiP+r5nU2f4PIHEcT7pDDFA54ijM8QAHT+O9uz76a1SPQkhISFERUWRn5/Pxo0bufzyy49Z12q1ola37YpGo0FRFBRFQa/XM2DAgCOG6pYuXcqwYcPOSP+7iuJtdXz07LrWpKfXiHiufXxQuyQ9HrOZsnvvxWMyYezdm/gXXmiT9Pyvv278KwVNBUQYI5gzfI4v6VEU+PZBsNTiiujB9QW+29pnje1Gn8TQX91HIYQQXcspj/iYzWYKCgpavy8qKmLLli2Eh4eTnJzMp59+SlRUFMnJyWzfvp3777+fqVOnMmHChNZnpk+fTkJCQuturClTpjB37lz69evXOtX15JNPctlll6E5cKjd7Nmzufnmmxk4cCAXXHAB//znPykpKeHuu+/+tT+DLsna7GTlf/eSv7EGgOBII6NvyiIpK7xd2le8Xip+9wiO/AK0UVEk/v3vqI3HHp1ZXrqcj/d8DMDzI54nwu/AVOSOzyFvAYpayzPaWTTYVfRNCuXe0Rnt0k8hhBBdyyknPhs3bmTMmDGt3x9cQzNjxgzeffddKisrmT17NtXV1cTFxTF9+vQ263QASkpK2ozwPHHgqoInnniC8vJyoqKimDJlCs8//3xrneuuu476+nqeffZZKisrycnJYeHChaSkpJxy0F2Z2+lh6w+lbPpuPy67B5UKcsclM3hKGrp2XCtT+8rfMP/4Iyq9nsR//B1dTPQx6zbaG3lq9VMATO81neEJw30FzZXwrW+Ka2fmXby/LQyDVs3ca3PRnsbBiUIIIcQpJz6jR49GUZRjls+aNYtZx7ho8qDly5e37YRWy1NPPcVTTz113Ofuvfde7r333pPuqzhE8Srs3VDN2gX7MDf4tvhHJQcx+sYeRKf8+mmtw5m++Zb6f/4TgLjn5+DXp89x6/9x3R9psDeQGZrJ/f3vP9BhBb6eBfYmXDG5TN/rS4YeHN+d9KjAdu2vEEKIrkPu6urkvB4vhVvq2LxkPzX7WwAIDDMwdGoG3QfFoGrnQ/9sO3dS+bjvFvWIO24nZMqU49ZfVrKM74q/Q6PSMGf4HPQava9g2yeQvwQ0Bl7Qz6LBDjkJwdw+Qm5dF0IIcfok8emkHDY3eSsr2P5jGS0NvjvPdAYN/Sem0HdsEtozsAXc3dhI+W9noTgcBI4aRdQDDxy3vsVr4a8b/grAbTm3kR2ZfaCgHr7zXVeyt+e9vLPRD61axUtXyRSXEEKIX0cSn07EYXVRuquR/dvr2Le5FpfDd8KyMVBHzsgEeo9OxD9Yf0beW/F4qHjoYVwVFehSkon/80tHvW39cF/bvqbR1UhmaCZ35x62SH3J42BrwBPVi+m7hwIe7hqVTq/49p2SE0II0fVI4tMBeV3QVGXFbvFgbXLQXG+nbHcjlftMKN5D66/C4wPIvSiJ7oNjzsgIz+Fq//YqltWrUfn5kfjaa8e8ePSgpSVL2eHagUal4fkRzx+a4tr3I2z9CFAxL2gWVaUe0qMC+O1F3c5o/4UQQnQNkvh0MB8/s4HmuiD++/2mo5aHxfqTkhNBap9I4ruFHvdk6/bSvHTpocXMc57D2L37ces32Bt4YYPvKIPbsm+jV0QvX4HLBt88CEBlj+n8dWswKhW8dFUfjLpzczqzEEKIzkUSnw5GZ/QlAHo/DQEhBgJCfV8xqcGk5EQQHOl3VvvjKCyk8lHfepzwGTMIueSSEz7z141/pcnRRIw6htuzD7u49qeXoLEIJTieuyomAV6mDU5mYGr7nC0khBBCSOLTwUy6O4cffvqeSy+bfM6PH/darZTNmoXXYsF/0CCiH37ohM9sqNrAV/u+QoWKqf5T0WkOxFC1A1a/CsCSlN+xbYOX8AA9j1ycdSZDEEII0cXIFpkOxj9Ej/o8SFcVRaHqmWdwFuxDGxVFwstzUZ0gEXN5XDy31ncB6VWZV5GkTTrYGHw7G7xu7JmXMntrPACPTsoixL/j3i0jhBDi/COJjzgtTZ99hmnBV6DRkDD3r2gjI0/4zLs736XIVES4MZzf9P3NoYJt/4XSdaAL4DnPdCxOD/2TQ7m6f+IZjEAIIURXJImPOGX2Xbuofm4OAFEP3I//oEEnfKa0pZQ3t70JwO8G/Y5g/YFdX44WWPr/ACjOuY8PdrlRq+C5qTmo2/lwRSGEEEISH3FKPGYzZQ88gOJ0EjhqFBEzZ57wGUVReH7d8zg8DobEDeGStEMLoNWrXgZzFUpYOnfnDwZg+gWpZMeHnLEYhBBCdF2S+IiTpigKlU88iWt/Cdr4OOL+9AIq9Yk/Qkv3L2VV+Sp0ah1PDHmidYt9gL0K9brXAViYeD+765xEBup5cPzxt8MLIYQQp0sSH3HSGj/6iJbvvgOdjsSXX0YbFnbCZ2xuGy9teAmA23vfTmpIamtZTvkHqLwuHGnjeGRrLACPTupJiJ8saBZCCHFmSOIjTop9925q/vQiANEPzcYvN/eknvv3zn9Tba0mLiCO23Jua31dlb+E2OatKGodfzfchsXpoU9iCFf2Szgj/RdCCCFAEh9xErxWK+UPzvat6xk9mvAZM07quRprDe/seAeABwc8iFFr9BW4nWi+fwKA+pzb+MdW38tPXNJLFjQLIYQ4oyTxESdUNed5nEVFaGNiiHvhjyd9Dcarv7yKzW0jNyqXiakTDxVsfBtVQyF2bQh/qJ+IV4GJ2bEMTpMTmoUQQpxZkviI4zJ9/TWmL74AtZr4P790Uut6APLq8/hq31cAPDLokUPJkt3ku5oC+CH0apbss6HTqHh0kpzQLIQQ4syTxEcck7O4mKqnngYg8t57CRg8+KSeUxSFlza8hILC5LTJ9Inqc6hw1atga8Ab0Y0/1o0C4OahqaRGBrR394UQQogjSOIjjkpxOil/6GG8Viv+gwYRec/dJ/3sspJlbKrehEFj4IH+DxwqaKmCNf8AYEXSPZTZtYT4aZk1NrOdey+EEEIcnSQ+4qhqXvkb9p070YSEEP/nl1BpNCf1nNPj5K8b/wrAjOwZxAXGHSpc/gK4bXgSBvPwdt89XfeNziDUX9/u/RdCCCGORhIfcQTzylU0vOPbjRX3x+fRxcae9LOf7f2MMnMZkX6RzMw57FTn2r3wy3sAfBFxJ3UWF5EGhRsHJ7Vr34UQQojjkcRHtOGur6fi0UcBCJt2A0Fjx570sza3jbe2vwXAPbn34K/zP1S47BlQPLgyJ/LsVt89XZOTvei18hEUQghx9si/OqKVoihU/OEPeOrqMHTrRvQjj5zS85/s/oQ6Wx0JgQlckXnFoYLS9bD7G1Cp+Y//DFocbnrEBNIvQmnnCIQQQojjk8RHtGp87z0sP61AZTAQ/9e/oDYaT/pZi8vC2zveBuDu3LvRaQ67dmLZswDYsm/gL5t9H7kHxmYiZxUKIYQ42yTxEcCBKyn+/BcAon//CMbup3ZR6Pt579PkaCI1OJVL0y89VFD0MxT/DBo9b6iuxubykJsUytisqPbsvhBCCHFSJPERvispZj+E4nIReNFFhN1wwyk9b3KY+PfOfwNwb9970aq1hwp/8t3vZcm+gdd/cQDw8ITuJ336sxBCCNGeJPERVL/wAs7CQrTR0cQ9P+eUk5J/7/w3La4WMkMzuTj14kMFxat8oz1qHa85p+D0eBmSFs6IzMh2jkAIIYQ4OZL4dHHNixbR9OlnoFIR/9LJX0lxUIO9gfd3vQ/Ab/r9BrXqsI/UT3/yvUevG3hrmwuA313cQ0Z7hBBCnDOS+HRhzrJyKv/fUwBE3HknAUOHnHIb72x/B5vbRq+IXlyUdNGhgv2roWgFqHW8Yr8Uj1dhdI8oBqbKRaRCCCHOHUl8uijF7abid7/D29KCX24uUb+575TbaLI38d+9/wXgvr73tR3JWe4b7WnpdT3v7nQDMHv8qS2YFkIIIdqbJD5dVO0//oFt82bUgYHE//UvqHS6Ez/0Pz7a8xE2t42s8CwuTLjwUEHJWij6CdQ65rkvx6vAqO5R9EkMbb8AhBBCiNMgiU8XZF65ivo33gQg9umn0ScmnnIbNreNj3Z9BMBtObcddbTHmn09b2/3jfb85iK5iFQIIcS5J4lPF+OqrqbikUdAUQi95hpCLr3ktNr5Mv9LGh2NJAQmMD5l/KGC0g1Q+COotbylXIHT42VwWjiDZG2PEEKI84AkPl2I4nZT/tBDeBoaMGRlEfP4H06rHbfXzX/y/gPALdm3tD23Z/XfALD3uobXtzoB+M0YGe0RQghxfpDEpwup/dur2DZuQh0QQOIrL5/SlRSHW1y8mHJzOeHGcKZmTj1UUL8Pdn0DwIeay7C7vPRJDOHCbnJujxBCiPODJD5dRMvy5dS/5bs5PW7Oc+hTU0+rHUVR+NeOfwEwLWsaRu1hydOafwAKrozxvLxFA8B9YzLl3B4hhBDnDUl8ugBnWTmVv38UgLBp0wieNOm021pVsYo9jXvw0/pxfdb1hwosdbDlAwC+DriaFoeb7jGBjO8Z86v6LoQQQrQnSXw6OY/ZTNk9d+MxmTBmZxP96O9/VXvv7HgHgKu7X02IIeRQwfq3wG3HE9eXOTt8C5nvHZ2JWq5gF0IIcR6RxKcTU9xuyh+cjSO/AG1UFIn/+Dtqvf6029tZt5MNVRvQqrRM7zX9UIHTCht802g/R02jweoiOdyfS/vE/doQhBBCiHYliU8nVv3Cn7D8/DMqo5HEefPQxcb+qvY+3P0hABPTJhIbcFhbWz8Eaz1KaArPFWQAcPuFaWg18vESQghxfpF/mTqphvc/oPED35qb+JdexK93zq9qr8nexHdF3wG0Xdvj9cDqvwOwO/Vm9jU4CPHTcfWAUz8UUQghhDjTJPHphMwrVlD9xz8CEDV7NsETJvzqNucXzMfpddIzvCd9IvscKtj9LTQWgTGUFyoHAnDD4GT89dpjtCSEEEKcO5L4dDLN331H2X2/Aa+XkCuvJOKO2391m17Fyyd7PgF8oz1ttqev+QcANVk3s2K/Fa1axYxhKb/6PYUQQogzQRKfTqThvfcpf3A2istF0PhxxD39VLucobOqfBVl5jKCdEFMSjtsK3zlNihdC2ot/zCPAeCSPnHEhfj96vcUQgghzoRTTnxWrFjBlClTiI+PR6VSMX/+/Dbl1dXV3HLLLcTHx+Pv78/EiRPJz88/bpujR49GpVId8XXJJYfukXr66aePKI/9lYt1OwtFUaj561yqn38eFIWwaTeQ8MorqH7FDq7DHRztuTzzcvy0hyU1G/4PAFvmJXyQ5wBg5oi0dnlPIYQQ4kw45YUYFouF3Nxcbr31Vq666qo2ZYqiMHXqVHQ6HQsWLCA4OJi5c+cybtw48vLyCAgIOGqbX3zxBU6ns/X7+vp6cnNzueaaa9rUy87O5vvvv2/9XqPRnGr3Ox13bS01f/kLpgVfARD1wP1E3HVXu52WXG4uZ0XZCgCu63HdoQJbE2z/FIAvtZNwexUGp4XTJzG0Xd5XCCGEOBNOOfGZNGkSk45x8m9+fj5r165lx44dZGdnAzBv3jyio6P56KOPuP32o683CQ9ve3P3xx9/jL+//xGJj1arlVEeQOV2Y16yBPNXX2NeuRI8HtBoiHv2GUL/Jxn9tT7d8ykKChfEXUBqSOqhgq0fgcuKN6onL+aFA24Z7RFCCHHea9c1Pg6Hb7rDeNjllxqNBr1ez8qVK0+6nbfffpvrr7/+iBGi/Px84uPjSUtL4/rrr6ewsLB9Ot5BeO12ap//I+lznqfqoYcx//QTeDz49e1L0j/fbPekx+lx8kX+FwBcl3XYaI/X2zrNtSHqSkx2NykR/oyT6ymEEEKc59p1z3FWVhYpKSk89thjvPnmmwQEBDB37lyqqqqorKw8qTbWr1/Pjh07ePvtt9u8PmTIEP7zn//QvXt3qqurmTNnDsOGDWPnzp1EREQctS2Hw9GajAE0NzcD4HK5cLlcpxnluaOo1VjXr0djs6GJiSH4sikETZmCPs030tLeMS0sWkijo5EY/xiGxQxrbV9V9BPa+gIUfSDP7s8BFKYPTcbrceP1nFzbB9vqiH8OB3WGGEDiOJ90hhigc8TRGWKAzhdHe1ApiqKc9sMqFV9++SVTp05tfW3Tpk3MnDmTrVu3otFoGDduHGq1b2Bp4cKFJ2zzrrvuYvXq1Wzfvv249SwWCxkZGTzyyCPMnj37qHWefvppnnnmmSNe//DDD/H39z9hX85HAbt2o2g0WDMzQH1mN+X9s+WflHhKGGscyxjjmNbXBxf+jTjTJjYHj+OKmtswqBWeHejBKEuuhBBCnAFWq5Vp06ZhMpkIDg7+VW21+ylzAwYMYMuWLZhMJpxOJ1FRUQwZMoSBAwee8Fmr1crHH3/Ms88+e8K6AQEB9O7d+7g7xh577LE2SVFzczNJSUmMGTPmmKNE5zvX+PEsXbqU8ePHo9Ppztj7lLaUUvJ1CWqVmkcmPkKUf5SvoLkc7ZbNAHwdegPUwJUDkrhySq9Tat/lcp2VOM6kzhADSBznk84QA3SOODpDDNB54qivr2+3ts7Y8bohIb6bu/Pz89m4cSPPPffcCZ/573//i8Ph4KabbjphXYfDwa5du7jwwguPWcdgMGAwGI54XafTdegPAJz5GBaXLAZgaNxQ4kPiDxVseQ8UL86kYfynwB9QuPmCtNPui/xZnD8kjvNHZ4gBOkccnSEG6PhxtGffTznxMZvNFBQUtH5fVFTEli1bCA8PJzk5mU8//ZSoqCiSk5PZvn07999/P1OnTmXCYdcmTJ8+nYSEBF544YU2bb/99ttMnTr1qKMxDz/8MFOmTCE5OZmamhrmzJlDc3MzM2bMONUQxAkoisLXhV8DMCVjyqECtxN++TcAywIvw+1V6J8cSq/4XzfsKIQQQpwtp5z4bNy4kTFjDq33ODiVNGPGDN59910qKyuZPXs21dXVxMXFMX36dJ588sk2bZSUlLSu+zlo7969rFy5kiVLlhz1fcvKyrjhhhuoq6sjKiqKoUOHsnbtWlJS5HqE9ra1diulLaX4af24KOmiQwV7vgVLLUpgLC8UZgAubhoqP38hhBAdxyknPqNHj+Z466FnzZrFrFmzjtvG8uXLj3ite/fux233448/Puk+il/nq32+wxDHp4zHX3fYIvDN7wNQnDSVks0uQv11TO4ddy66KIQQQpwWuatLtOH0OPmu+Dvgf6a5mitg3w8AvGkaAsA1AxIx6mQrlxBCiI5DEh/Rxk9lP9HibCHGP4ZBMYMOFWz9GBQvjvjBfFLkWzA+bYhMcwkhhOhYJPERbXy9z7eo+ZL0S9CoD4zmKAps+QCAH4zjURS4sFskaZFHv3tNCCGEOF9J4iNaNdob+bnsZwCmpB82zVW6HuoLUHT+vLA/C4AbZbRHCCFEBySJj2j1XfF3uBU3PcN7khmWeahgi29Rc1ncBEosGmKCDYztGX2OeimEEEKcPkl8RKuD01xtFjU7rbDjSwA+cIwA4JoBSeg08tERQgjR8ci/XgKAYlMx2+u2o1FpmJQ26VDBrq/B2YI7JIW3SmIBuGpA4jnqpRBCCPHrSOIjAFiy33dw5ND4oUT6RR4qODDNtTlsEh5FzcCUMFnULIQQosOSxEcAsKxkGQDjkscderFxPxStQEHF3NoBAFwzUEZ7hBBCdFyS+AgqzZXk1eehQsWYpEPXkbD1IwDMccNYUx+AUaeWk5qFEEJ0aJL4CH4o9Z3I3C+6HxF+By6IVRTfoYXAQt1YACbnxBFk7Li3+wohhBCS+IjWaa6xyWMPvVixGRqLUHT+/LXEt7X9alnULIQQooOTxKeLa7Q3sql6EwBjUw5LfHZ+AUB59Chq7FoSQv0Ymh5xLroohBBCtBtJfLq45aXL8SpessKzSAhM8L3o9bae3fOZfTDg28KuVqvOUS+FEEKI9iGJTxf3Q4lvfc9FyRcderFsPTSX4dUH8XpFGgBX95dpLiGEEB2fJD5dmNVlZXXFauB/1vfs+ByAvWEjcSh6hqSFkxzhfy66KIQQQrQrSXy6sJXlK3F6nSQFJdEttJvvRa8Hds4H4N1m39k9sqhZCCFEZyGJTxd2+G4ulerA+p3ilWCpwW0I5bPGTAxaNZPk7B4hhBCdhCQ+XZTL42JF2Qrg6NNc24NH4UbL2J7RBBq056KLQgghRLuTxKeLWle1DrPLTKRfJH2i+vhe9Lhg11cAvNPUH4ApfeLPVReFEEKIdieJTxd1cJprTNIY1KoDH4PC5WBrxGWM5NuWDAINWsZkRZ+7TgohhBDtTBKfLkhRlNZprjbb2A9Mc/0SOAovaiZkx2DUac5FF4UQQogzQhKfLqjIVESNtQa9Ws/AmIG+F1122PUNAP9s7AfAlFyZ5hJCCNG5SOLTBa2pXANA/5j+GLVG34v7fgBnCw7/WH6wpBLmr2NEZuQ57KUQQgjR/iTx6YLWVPgSnwviLzj04u5vAVhvHI6Cmsm949Bp5OMhhBCic5F/2boYl9fFhqoNAFwQdyDx8Xpg7yIA/tWQDcg0lxBCiM5JEp8uZlvtNqxuK+HGcHqE9/C9WLoOrPW4dCH8ZO9GbLCRwanh57ajQgghxBkgiU8Xc3Caa0jskEPb2A9Mc202DsaDhkv7xMlN7EIIITolSXy6mIMLm1vX9ygK7FkIwAdNOQBc1lemuYQQQnROkvh0Ic3OZnbU7QAOS3xq90BDIR61ju9dOaRE+NM7IeQc9lIIIYQ4cyTx6UI2VG7Aq3hJDU4lNiDW9+Ju39k9u/z6Y8GPSTlxhy4sFUIIIToZSXy6kCOmuaB1muuzFt99XRNzYs96v4QQQoizRRKfLqT1/J6D29ibK6F8EwD/v717D46qvv8//txsNgkJSeSSK4Q0QMItIX5NBINaQSE1tQ4ZW8dbS6zaGQULNPZbpY4VHUecdibTdixWp7S2f1SsVcHfiJX4E4KURrlFA6iEEghoLhIgCYkkm+zn+8eSLUtuBFN3z9nXYyYz5uzu8f3OO2Ffc87n7Hmz63KS46KYrdNcIiJiYwo+IeKzM59R11aH0+EkP/ncbSrOfXbP0VEz+YIxfGtWkq7mEhERW1PwCRG9R3tyxucQGxHr3XjuMvaNZ7335vrWLJ3mEhERe1PwCRF9blPR2Qa13ju0bzx7OZdFu5iToQ8tFBERe1PwCQE9nh7eb3gfOC/4HHoHerpojkzj3yaVhTOSCNe9uURExOb0ThcCPjn1CS2dLcS4Ysgen31uo/dqrn90XwE4dJpLRERCgoJPCNjd4L1yKz8pH1eYy3tT0kPlAGzoyCU6wsm1meMDWaKIiMjXQsEnBHz4xYcAXJ54uXfD51Xw5SnOOkezx2Qyf1oCUS5nwOoTERH5uij4hIDe4JObkOvdcOgdAD5w5NCDU6e5REQkZCj42FxDewONHY04HU5mjZvl3fjv/w/Api9n4XI6WDA9MYAVioiIfH0UfGyu6osqALLGZBHtioYvT8HxnQBs65nNvCnjiYtyBbBCERGRr4+Cj8192HTBaa7DFWA8HHNO5HPG695cIiISUhR8bM63vifxXPA5d5prc2cOADfoNJeIiISQYQefbdu2cfPNN5OamorD4WDDhg1+jzc2NnL33XeTmppKdHQ0N954IzU1NYPuc/78+Tgcjj5fN910k9/z1q5dS0ZGBlFRUeTl5fHee+8Nt/yQ0tnTyccnPwbg8oTLwRg49C4A2zyzyZkQT2JcVAArFBER+XoNO/i0t7eTm5vLs88+2+cxYwzFxcUcPnyYjRs3snfvXtLT01m4cCHt7e0D7vO1116jvr7e97Vv3z6cTie33nqr7zkvv/wyK1eu5NFHH2Xv3r1ce+21FBUVUVdXN9wWQsaB5gN0e7oZFzWOCaMnwBefQutxuhwRVHpmaFGziIiEnPDhvqCoqIiioqJ+H6upqaGyspJ9+/Yxa5b3CqK1a9eSmJjISy+9xH333dfv68aO9b9H1Pr164mOjvYLPmVlZdx7772+ffz617/m7bff5rnnnmPNmjXDbSMknL++x+Fw+E5z7fTMoJMIFkxLCGR5IiIiX7thB5/BdHZ2AhAV9Z/TJ06nk4iICLZv3z5g8LnQunXruP3224mJiQGgq6uL3bt388gjj/g9r7CwkB07dgxaT29NAK2trQC43W7cbvfFNRVkeuu+mPr3NO4BIHtcNm63G2fNO4QB73ZnMzbGxcykmID9HIbTR7CyQw+gPoKJHXoAe/Rhhx7Afn2MhBENPtOnTyc9PZ1Vq1bx/PPPExMTQ1lZGQ0NDdTX11/UPj744AP27dvHunXrfNtOnDhBT08PSUlJfs9NSkqioaFhwH2tWbOGJ554os/2LVu2EB0dfZFdBafy8vJBHzfGsLPVe9l6x6EO/nF4A9+u9a6JqvDkMiW6k3/8463/ep1DGaoPK7BDD6A+gokdegB79GGHHsD6fXR0dIzYvkY0+LhcLl599VXuvfdexo4di9PpZOHChQOeGuvPunXryM7OZs6cOX0eczgcft8bY/psO9+qVasoLS31fd/a2kpaWhoLFixg3LhxF11TMHG73ZSXl7No0SJcroE/f+fzM59z5o0zhDvCuefb9zDq6A6cH7ppcozjkJnAbxbk8u2cwF3KfrF9BDM79ADqI5jYoQewRx926AHs00dzc/OI7WtEgw9AXl4eVVVVtLS00NXVRUJCAnPnziU/P3/I13Z0dLB+/XqefPJJv+3jx4/H6XT2ObrT1NTU5yjQ+SIjI4mMjOyz3eVyWfoXAIbuYd+pfQBMHzud2FGxcKQCgHfdOTjDwpg/IzkofgahMAurUB/Bww49gD36sEMPYP0+RrL2/9rn+MTHx5OQkEBNTQ27du1i8eLFQ77mb3/7G52dnXz/+9/32x4REUFeXl6fQ3Xl5eXMmzdvROu2iz43Jj23sLnCk0te+hjiR1n3D0BERORSDfuIz5kzZzh06JDv+9raWqqqqhg7diyTJk3ilVdeISEhgUmTJlFdXc2KFSsoLi6msLDQ95olS5YwYcKEPldjrVu3juLi4n5PQ5WWlvKDH/yA/Px8CgoKeOGFF6irq+P+++8fbgshwe/GpC3H4YtP8BDGPz2zWKrL2EVEJEQNO/js2rWLBQsW+L7vXUNTUlLCiy++SH19PaWlpTQ2NpKSksKSJUt47LHH/PZRV1dHWJj/waaDBw+yfft2Nm/e3O//97bbbqO5uZknn3yS+vp6srOz2bRpE+np6cNtwfY63B18evJT4FzwqfEe7akyU2hlNNcr+IiISIgadvCZP38+xpgBH1++fDnLly8fdB9bt27tsy0rK2vQ/QIsXbqUpUuXXlSdoWx/8356TA+J0YkkxyTDke0A/LNnFhMuG0Vm4ugAVygiIhIYuleXDZ1/mssBvuDj/bTmhEGvhBMREbEzBR8b8rsj+6kj0HIMN+Hs8WTqNJeIiIQ0BR8bOtB8AIDZCbN9R3uqPJPxhEdTMHl8IEsTEREJKAUfmzl19hRNXzYBkDUm67zTXDO5avI4RkU4A1meiIhIQCn42Mynp7xXc6XFphETHu0LPv/yzOTaTB3tERGR0KbgYzMHTx4Ezh3tOXUEWo/TZZzs8WRybabuxi4iIqFNwcdmeo/4TBsz7T/re8xUYmPjyErSZewiIhLaFHxspuZUDdC7vsd7N/ZKzwyumTpel7GLiEjIU/CxEbfHzaHT3tuJXLiw+eqpWt8jIiKi4GMjR1uO4va4iXHFMMHdCa2f+db3XKPgIyIiouBjJ73rezIvyyTs6A7Au75nYuI4kuOjAlmaiIhIUFDwsZGDp7xXdE0bO83vNhU62iMiIuKl4GMjvUd8si7L9Fvfo+AjIiLipeBjIzUnz13RFR4LrZ/RacL5kCyumjIuwJWJiIgEBwUfmzj/VhWZJz8DvOt7ZqQlMToyPJCliYiIBA0FH5voXd+TFptGzLEPgHPre3SbChERER8FH5v49OS59T1jsjDnruh6XwubRURE/Cj42ITviq7oFBytx+kxDmpc08lNuyywhYmIiAQRBR+b6A0+We5uAD426eROTsXl1IhFRER66V3RBro93f+5VcXpRgB269OaRURE+lDwsYEjLUd8t6pI+bwagN2eLN2fS0RE5AIKPjbQe5orM34KYQ3e4FMbNZOpiaMDWZaIiEjQUfCxgd5PbJ4WMYYwj5smcxkTM6bjcDgCXJmIiEhwUfCxAd/C5i434F3fM2eyPq1ZRETkQgo+NnDw5LlTXacaAO/6nrkZCj4iIiIXUvCxuPNvVZHx2T4ADrpmMD05NpBliYiIBCUFH4vz3aoiOpnLvjxBpwknOj2PsDCt7xEREbmQgo/F1bbUAjAl3HuEZ5/JIH9qciBLEhERCVoKPhZ3tPUoAJO6OgHY48lkTsbYQJYkIiIStBR8LK6urQ6A1HMLm/c7pzMzJS6QJYmIiAQtBR+Lq2v1Bp/Jpz/zbph4JeG6P5eIiEi/9A5pYd2ebo6fOQ7AN9xdHPMkMC1zWoCrEhERCV4KPhZW315Pt6ebCMJI6ulhj8lk7mSt7xERERmIgo+F9Z7mmmCchAEfOaaRMyE+sEWJiIgEMQUfC+td2Dzpy3YAOpPzcWl9j4iIyID0LmlhvUd8MrrO0mEiSc7KC3BFIiIiwU3Bx8J8n+HT7eYjM5k5U5ICXJGIiEhwU/CxsGNtxwCY5O6mmqnMnqj1PSIiIoNR8LGobk83x9u8l7Knu7s5MzaHKJczwFWJiIgENwUfi6pvr6fbdBPpMST29BD9jSsDXZKIiEjQU/CxqN7TXGndblrMaKZmzQxwRSIiIsFPwceifJeyu7up9mTwP+n64EIREZGhKPhYVO8Rn3R3N3VR0xgbExHgikRERIKfgo9F9R7xSet24076nwBXIyIiYg3DDj7btm3j5ptvJjU1FYfDwYYNG/web2xs5O677yY1NZXo6GhuvPFGampqhtzv6dOnWbZsGSkpKURFRTFjxgw2bdrke3z16tU4HA6/r+Tk5OGWbxvH2ryf4ZPu7uayqXMCXI2IiIg1hA/3Be3t7eTm5vLDH/6Q7373u36PGWMoLi7G5XKxceNG4uLiKCsrY+HChRw4cICYmJh+99nV1cWiRYtITEzk73//OxMnTuTYsWPExsb6PW/WrFm88847vu+dztC8fLvH9PBZ++cAjO4axYxpuiO7iIjIxRh28CkqKqKoqKjfx2pqaqisrGTfvn3MmjULgLVr15KYmMhLL73Efffd1+/r/vjHP3Ly5El27NiBy+UCID09vW+x4eEhfZSn12nPaXqMh0iPh8aeDL6ZFBfokkRERCxh2MFnMJ2dnQBERUX5tjmdTiIiIti+ffuAweeNN96goKCAZcuWsXHjRhISErjzzjt5+OGH/Y7q1NTUkJqaSmRkJHPnzuXpp59m8uTJg9bTWxNAa2srAG63G7fb/ZV6DRS3202zpxmAtO5uTsZl4+npxtMT4MKGqffnb9U5gD16APURTOzQA9ijDzv0APbrYyQ4jDHmkl/scPD6669TXFwMeAvLzMxkzpw5PP/888TExFBWVsaqVasoLCzk7bff7nc/06dP58iRI9x1110sXbqUmpoali1bxooVK/jFL34BwFtvvUVHRwdZWVk0Njby1FNP8cknn7B//37GjRvX735Xr17NE0880Wf7X//6V6Kjoy+17YD7V+e/ePPLN7mhvYOpZ0tIm5Ib6JJERET+azo6OrjzzjtpaWkhLu6rneUY0eADsHv3bu69914+/PBDnE4nCxcuJCzMu4b6/MXK58vKyuLs2bPU1tb6jvCUlZXxq1/9ivr6+n5f097ezpQpU/jZz35GaWlpv8/p74hPWloa9fX1A4alYOd2u/nJhmXscO/ih6dbmXPN/2Pu7BmBLmvY3G435eXlLFq0yHd602rs0AOoj2Bihx7AHn3YoQewTx/Nzc2kpKSMSPAZ0VNdAHl5eVRVVdHS0kJXVxcJCQnMnTuX/Pz8AV+TkpKCy+XyO601Y8YMGhoa6OrqIiKi72fUxMTEkJOTM+gVY5GRkURGRvbZ7nK5LP0LcNLtXdgc544gZ+YMS/di9VmAPXoA9RFM7NAD2KMPO/QA1u9jJGv/r32OT3x8PAkJCdTU1LBr1y4WL1484HOvvvpqDh06hMfj8W07ePAgKSkp/YYe8B7N+fjjj0lJSRnx2oPdiZ6TAISZCcSPsu4vsoiIyNdt2MHnzJkzVFVVUVVVBUBtbS1VVVXU1Xk/UO+VV15h69atHD58mI0bN7Jo0SKKi4spLCz07WPJkiWsWrXK9/0DDzxAc3MzK1as4ODBg7z55ps8/fTTLFu2zPecn/70p1RUVFBbW8v777/P9773PVpbWykpKbnU3i3J7XFzMuwsADGx2QGuRkRExFqGfapr165dLFiwwPd97/qakpISXnzxRerr6yktLaWxsZGUlBSWLFnCY4895rePuro637ofgLS0NDZv3sxPfvITZs+ezYQJE1ixYgUPP/yw7znHjx/njjvu4MSJEyQkJHDVVVdRWVnZ72XvdlZ/ph6PA6I8HhInXR3ockRERCxl2MFn/vz5DLYeevny5SxfvnzQfWzdurXPtoKCAiorKwd8zfr16y+6RjurPXEA8F7Knp6t4CMiIjIculeXxez79z8BSHA7+cbEtABXIyIiYi0KPhZT+4X3iE8cYwkLcwS4GhEREWtR8LGYL7q8n2s0blRGgCsRERGxHgUfiznh6ADgG6lXBLgSERER61HwsRBPj4fro6/i2rZ45ub0f6NYERERGZiCj4WEOcN46Nbf8620/2ViwpRAlyMiImI5Cj4iIiISMhR8REREJGQo+IiIiEjIUPARERGRkKHgIyIiIiFDwUdERERChoKPiIiIhAwFHxEREQkZCj4iIiISMhR8REREJGQo+IiIiEjIUPARERGRkKHgIyIiIiEjPNAFfJ2MMQC0tbXhcrkCXM2lcbvddHR00NraatkewB592KEHUB/BxA49gD36sEMPYJ8+2tragP+8j38VIRV8mpubAcjIyAhwJSIiIjJczc3NxMfHf6V9hFTwGTt2LAB1dXVf+QcXKK2traSlpXHs2DHi4uICXc4ls0MfdugB1EcwsUMPYI8+7NAD2KePlpYWJk2a5Hsf/ypCKviEhXmXNMXHx1v6FwAgLi7O8j2APfqwQw+gPoKJHXoAe/Rhhx7APn30vo9/pX2MQB0iIiIilqDgIyIiIiEjpIJPZGQkjz/+OJGRkYEu5ZLZoQewRx926AHURzCxQw9gjz7s0AOoj/44zEhcGyYiIiJiASF1xEdERERCm4KPiIiIhAwFHxEREQkZCj4iIiISMkIm+Kxdu5aMjAyioqLIy8vjvffeC3RJw7J69WocDoffV3JycqDLGtS2bdu4+eabSU1NxeFwsGHDBr/HjTGsXr2a1NRURo0axfz589m/f39gih3EUH3cfffdfWZz1VVXBabYAaxZs4Yrr7yS2NhYEhMTKS4u5tNPP/V7jhXmcTF9BPs8nnvuOWbPnu37QLmCggLeeust3+NWmAMM3Uewz6E/a9asweFwsHLlSt82q8zjfP31YYV5DPU+N1KzCIng8/LLL7Ny5UoeffRR9u7dy7XXXktRURF1dXWBLm1YZs2aRX19ve+ruro60CUNqr29ndzcXJ599tl+H//lL39JWVkZzz77LDt37iQ5OZlFixb5bkYXLIbqA+DGG2/0m82mTZu+xgqHVlFRwbJly6isrKS8vJzu7m4KCwtpb2/3PccK87iYPiC45zFx4kSeeeYZdu3axa5du7j++utZvHix7x9wK8wBhu4DgnsOF9q5cycvvPACs2fP9ttulXn0GqgPsMY8BnufG7FZmBAwZ84cc//99/ttmz59unnkkUcCVNHwPf744yY3NzfQZVwywLz++uu+7z0ej0lOTjbPPPOMb9vZs2dNfHy8+f3vfx+ACi/OhX0YY0xJSYlZvHhxQOq5VE1NTQYwFRUVxhjrzuPCPoyx5jzGjBlj/vCHP1h2Dr16+zDGWnNoa2szmZmZpry83Fx33XVmxYoVxhjr/V0M1Icx1pjHYO9zIzkL2x/x6erqYvfu3RQWFvptLywsZMeOHQGq6tLU1NSQmppKRkYGt99+O4cPHw50SZestraWhoYGv7lERkZy3XXXWW4uAFu3biUxMZGsrCx+9KMf0dTUFOiSBtXS0gL858a9Vp3HhX30sso8enp6WL9+Pe3t7RQUFFh2Dhf20csqc1i2bBk33XQTCxcu9NtutXkM1EcvK8xjoPe5kZyF7W9SeuLECXp6ekhKSvLbnpSURENDQ4CqGr65c+fyl7/8haysLBobG3nqqaeYN28e+/fvZ9y4cYEub9h6f/b9zeXo0aOBKOmSFRUVceutt5Kenk5tbS2PPfYY119/Pbt37w7KT0s1xlBaWso111xDdnY2YM159NcHWGMe1dXVFBQUcPbsWUaPHs3rr7/OzJkzff+AW2UOA/UB1pgDwPr169mzZw87d+7s85iV/i4G6wOsMY/B3udGcha2Dz69HA6H3/fGmD7bgllRUZHvv3NycigoKGDKlCn8+c9/prS0NICVfTVWnwvAbbfd5vvv7Oxs8vPzSU9P58033+SWW24JYGX9e/DBB/noo4/Yvn17n8esNI+B+rDCPKZNm0ZVVRWnT5/m1VdfpaSkhIqKCt/jVpnDQH3MnDnTEnM4duwYK1asYPPmzURFRQ34vGCfx8X0YYV5DPY+17sQeyRmYftTXePHj8fpdPY5utPU1NQnOVpJTEwMOTk51NTUBLqUS9K7Ut9ucwFISUkhPT09KGfz4x//mDfeeIMtW7YwceJE33arzWOgPvoTjPOIiIhg6tSp5Ofns2bNGnJzc/nNb35juTkM1Ed/gnEOu3fvpqmpiby8PMLDwwkPD6eiooLf/va3hIeH+37mwT6Pofro6enp85pgnMeFzn+fG8m/DdsHn4iICPLy8igvL/fbXl5ezrx58wJU1VfX2dnJxx9/TEpKSqBLuSQZGRkkJyf7zaWrq4uKigpLzwWgubmZY8eOBdVsjDE8+OCDvPbaa7z77rtkZGT4PW6VeQzVR3+CcR4XMsbQ2dlpmTkMpLeP/gTjHG644Qaqq6upqqryfeXn53PXXXdRVVXF5MmTLTGPofpwOp19XhOM87jQ+e9zI/q3MdxV11a0fv1643K5zLp168yBAwfMypUrTUxMjDly5EigS7toDz30kNm6das5fPiwqaysNN/5zndMbGxsUPfQ1tZm9u7da/bu3WsAU1ZWZvbu3WuOHj1qjDHmmWeeMfHx8ea1114z1dXV5o477jApKSmmtbU1wJX7G6yPtrY289BDD5kdO3aY2tpas2XLFlNQUGAmTJgQVH088MADJj4+3mzdutXU19f7vjo6OnzPscI8hurDCvNYtWqV2bZtm6mtrTUfffSR+fnPf27CwsLM5s2bjTHWmIMxg/dhhTkM5MKroawyjwud34dV5jHU+9xIzSIkgo8xxvzud78z6enpJiIiwlxxxRV+l79awW233WZSUlKMy+Uyqamp5pZbbjH79+8PdFmD2rJliwH6fJWUlBhjvJcnPv744yY5OdlERkaab37zm6a6ujqwRfdjsD46OjpMYWGhSUhIMC6Xy0yaNMmUlJSYurq6QJftp7/6AfOnP/3J9xwrzGOoPqwwj3vuucf3b1FCQoK54YYbfKHHGGvMwZjB+7DCHAZyYfCxyjwudH4fVpnHUO9zIzULhzHGXOJRKBERERFLsf0aHxEREZFeCj4iIiISMhR8REREJGQo+IiIiEjIUPARERGRkKHgIyIiIiFDwUdERERChoKPiIiIhAwFHxEREQkZCj4iIiISMhR8REREJGQo+IiIiEjI+D/ip73/yNFLGgAAAABJRU5ErkJggg==", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkcAAAGdCAYAAAAYDtcjAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAABll0lEQVR4nO3de1zW9f3/8cclZxAPKEcFVMQsNS1N1CzRkMRGktUqV2pp5aHUqDVN/UkttW9bzDZnWXOWlYdapVtzqR3QzEwxMURTFBJFQEXlKHABn98f5JUkqBh6HXjebzdu6/oc3rxevCme+xxNhmEYiIiIiAgAzaxdgIiIiIgtUTgSEREROYfCkYiIiMg5FI5EREREzqFwJCIiInIOhSMRERGRcygciYiIiJxD4UhERETkHM7WLsCWVFdXc/ToUby9vTGZTNYuR0RERC6BYRgUFRURFBREs2a//riPwtE5jh49SnBwsLXLEBERkctw+PBh2rdv/6vHUTg6h7e3NwCZmZn4+PhYuZrLZzabWb9+PdHR0bi4uFi7nMvmCH04Qg/gGH04Qg+gPmyJI/QAjtHHyZMn6dixo+Xv+K+lcHSOs6fSvL29adGihZWruXxmsxlPT09atGhht7/o4Bh9OEIP4Bh9OEIPoD5siSP0AI7Rh9lsBmi0S2J0QbaIiIjIORocjjZt2kRsbCxBQUGYTCZWr15da31eXh5jx44lKCgIT09Phg0bRnp6+gXHNJvNvPDCC4SFheHu7k7Pnj359NNPa22TkJCAyWSq9RUQEFBrG8MwSEhIICgoCA8PDyIjI0lLS2toiyIiItKENTgclZSU0LNnTxYuXHjeOsMwiIuLIyMjgzVr1rBz505CQ0OJioqipKSk3jFnzZrF4sWL+dvf/saePXuYMGECd911Fzt37qy1Xbdu3cjJybF8paam1lr/8ssvk5iYyMKFC9m+fTsBAQEMHTqUoqKihrYpIiIiTVSDrzmKiYkhJiamznXp6els3bqV3bt3061bNwAWLVqEn58fK1asYPz48XXu98477zBz5kyGDx8OwMSJE1m3bh2vvPIK77777s/FOjufd7ToLMMwWLBgATNnzmTkyJEAvP322/j7+7N8+XIef/zxhrYqIiIiTVCjXpBdXl4OgLu7u2WZk5MTrq6ubN68ud5wVF5eXmsfAA8PDzZv3lxrWXp6OkFBQbi5uREREcG8efPo1KkTUHOHWW5uLtHR0Zbt3dzcGDRoEFu2bKkzHJWXl1tqBigsLARqTvOdvbjLHp2t3Z57AMfowxF6AMfowxF6APVhSxyhB3CMPhq7dpNhGMZl72wy8fHHHxMXFwfUFBceHk7fvn1ZvHgxXl5eJCYmMmPGDKKjo1m3bl2d44waNYpdu3axevVqwsLC+PzzzxkxYgRVVVWW8PK///2P0tJSunTpQl5eHi+++CI//PADaWlptGnThi1btnDzzTeTnZ1NUFCQZezHHnuMQ4cO1fm9ExISeP75589bvnz5cjw9PS/3xyIiIiJXUWlpKaNGjaKgoKBR7jZv1CNHLi4ufPjhh4wbNw4fHx+cnJyIioqq9zTcWa+++iqPPvooXbt2xWQyERYWxsMPP8zSpUst25w7Ro8ePejfvz9hYWG8/fbbxMfHW9b98jY+wzDqvbVvxowZtfYtLCwkODiYwYMH06ZNmwb1bkvMZjMbNmxg6NChdntbJjhGH47QAzhGH47QA6gPW+IIPYBj9JGfn9+o4zX6c4569+5NSkoKBQUFVFRU4OvrS0REBH369Kl3H19fX1avXk1ZWRn5+fkEBQUxffp0OnbsWO8+Xl5e9OjRw3In3NlrkXJzcwkMDLRsd+zYMfz9/escw83NDTc3t/OWu7i42O0vyLnUh+1whB7AMfpwhB5AfdgSR+gB7LuPxq77ij3nqGXLlvj6+pKenk5ycjIjRoy46D7u7u60a9eOyspKPvzwwwvuU15ezt69ey1BqGPHjgQEBLBhwwbLNhUVFWzcuJEBAwb8+oZERESkSWjwkaPi4mIOHDhg+ZyZmUlKSgo+Pj6EhITwwQcf4OvrS0hICKmpqUydOpW4uLhaF0qPHj2adu3aMX/+fAC+/fZbsrOz6dWrF9nZ2SQkJFBdXc2zzz5r2eeZZ54hNjaWkJAQjh07xosvvkhhYSFjxowBak6nTZs2jXnz5hEeHk54eDjz5s3D09OTUaNGXfYPSERERJqWBoej5ORkBg8ebPl89pqdMWPG8NZbb5GTk0N8fDx5eXkEBgYyevRoZs+eXWuMrKysWm/NLSsrY9asWWRkZNC8eXOGDx/OO++8Q6tWrSzbHDlyhAceeIATJ07g6+tLv3792Lp1K6GhoZZtnn32Wc6cOcOkSZM4deoUERERrF+/vtHetSIiIiKOr8HhKDIykgvd4DZlyhSmTJlywTGSkpJqfR40aBB79uy54D4rV668aG0mk4mEhAQSEhIuuq2IiIhIXfTiWREREbFbm5PW8dHuLY06psKRiIiI2I1tW75g4/db2e/jyd5W7clw7khlh6HAC432PRSORERExGa98vp8vvNvRY5nK3Kc/TjVrA1cM7zWNgHVuRxvxO+pcCQiIiI2JSf7MK+ue49NwV3IuOb8B0m3Mk4SUnGUrsdziTB5MyzmHto24vdXOBIRERGrKS4sYO3af3HgdC4nPF041rw5O1p15VTHYQA4G2Z6n0mlU34+gWcqGXDNjQwcfHutMWz+CdkiIiIiF7Ns2d/52MeNZM/rMfv3hl+8zKKFcZqB+d9zt6sfd4x45KrWpnAkIiIiV83St//KmjZebGvfj2qTEwBORiVtjHzaVp6kTXkhHfNPM+nWEXQYEmmVGhWORERE5Io4fCiDT9Z/RFYzM3ktvTjS3IfU4JsxfgpF3cv3EH3oEA/HjcHXv/53sF5tCkciIiLSKI7nHWXp6rf5wceLfa3a8aNTKFWdo87brkd5GjE/ZhE/YYYVqrw4hSMRERG5bDu2buL9PVv4PjCIH9w6c6ZL7bvLPIxSAqtyCSg7iX9xIdcWVjLl8WfrGc02KByJiIhIg+xO+ZZlO5PYFdCOPW5dMP90ZxmAl1HENWUZdDlxjN54cc89Y/Hw9LRitQ2ncCQiIiKXZOXKN/mvawVft+pJaYefb6f3r86l1+kD9DxezNh7x+PT5hYrVvnrKRyJiIhIvU7mn2DxB2/wVUg7UvxutNxh5lt9jD6n9tHnZBmPPDgZD89hFxnJfigciYiISC3btnzBx3u/ZU+AP2ke4RSf87qOruZ93PJjBk/d/Sg+baKtWOWVo3AkIiLSxB3L+ZG5S18mo20rDnoH1Nxl1unnC6u9jGJ6lezltpwCJj32jBUrvToUjkRERJqgbVu+4L2DO9gW0Jkfuwy0PHvoLP/qXLoXZtA97zTj7xqLr/9AK1V69SkciYiINBFnSkv567IFbG3vzw7PblS0H2pZ51edR9iZw3TKP0lPkxejR0+2YqXWpXAkIiLiwNL3p/Fe0hrSAtqy2yuMU+dcPxRUlU2/4/u47tgZHp/0LC4uLlas1HYoHImIiDiQM6WlrHj/TXY6m9nXNoAfXDtTEf5zIHIzyuhdspv+2cd4YvQ0nF2iWbt2rRUrtj0KRyIiInZsd8q3rN/6BYc9nTnUujX7PDuQHzqo1jatq/PpXnKAbrn5/C5yBOFdxlvWmc3mq12yzVM4EhERsSNnSktJXLaAb4MDOOQexDGTH8Y1tV/Z4WyYCavM5JpTR+lRUMH4B5/Aw/M2K1VsfxSORERE7MAX69fwQf5BtvhdR9451w0BtDBO086cS1DpSbocO83dEbfRvddvrVSp/VM4EhERsVFnSkv5y7K/sDU4iJ0e3TAHhALgbpzhpuJUrs85wY0BnbhjxP1WrtSxKByJiIjYmNUfvsPayny2tr2WY9fcYVkeUpXFgKP7eKhLBL2HPGbFCh2bwpGIiIiVbfxiLRsO7uKAb2vSm7cn26eHZZ27cYY+xbsZcDSfiaOn4OF5pxUrbRoUjkRERK6yk/knWPrBm6S19WJvy2B+dArF6BxUa5uwyoP0O3qQh7oPpNeQR61UadOkcCQiInKVvPL6fL5q70eaZzhFv7jDLKA6h/CSLDofP0lkQDi333GPlaoUhSMREZErKCf7MK+ue49NwV3IOCcQeRildC0/wLXHc7nFw4+77h1rvSKlFoUjERGRRnKmtJRP//s+qSezyfP2INfbmz1eYZzqOAyoef5Qn9LviTiSx9jYhwhsN8DKFUtdFI5ERER+pdUfvsNq4zRbfHpQ2LYXtO1Va30L4zQD87/nblc/7hgxzio1yqVTOBIREblMaT/uZMn6k+xo3YMqU82f1GZGFf5GHgEVJwgoKSDodDGPDr6LDkMirVusXDKFIxERkQY4nneUxE/eZlNIFw5eH2dZ3tl8gEFZ6TVBKGx4/QOIzWvW0B02bdpEbGwsQUFBmEwmVq9eXWt9Xl4eY8eOJSgoCE9PT4YNG0Z6evoFxzSbzbzwwguEhYXh7u5Oz549+fTTT2ttM3/+fG666Sa8vb3x8/MjLi6Offv21dpm7NixmEymWl/9+vVraIsiIiLnSfr8Ex5f+RcG7d7L0k4xHHQOw8moJKL0O6bvX8fm6HuYO34GHcK6WrtU+ZUafOSopKSEnj178vDDD3P33XfXWmcYBnFxcbi4uLBmzRpatGhBYmIiUVFR7NmzBy8vrzrHnDVrFu+++y5vvvkmXbt2Zd26ddx1111s2bKFG264AYCNGzcyefJkbrrpJiorK5k5cybR0dHnjTts2DCWLl1q+ezq6trQFkVERDhTWsq7Kxez08Ngr08Q6S5hVPq3B8DbKGDgye/pc7SYxyY+g4uLi5WrlcbU4HAUExNDTExMnevS09PZunUru3fvplu3bgAsWrQIPz8/VqxYwfjx4+vc75133mHmzJkMH15zGHLixImsW7eOV155hXfffRfgvCNJS5cuxc/Pjx07dnDrrbdalru5uREQENDQtkRERABY8s9X+bKVCyktwznRcXCtdcFVh7k1ey8T+t9Bh1smsHbtWitVKVdSo15zVF5eDoC7u7tlmZOTE66urmzevLnecFReXl5rHwAPDw82b95c7/cqKCgAwMfHp9bypKQk/Pz8aNWqFYMGDWLu3Ln4+fldVj8iItI07E39jiU7NvBtUBjpHQdZljsbZsLNB7n25FFuOGPiwfsfx8MzFqi5JEQcU6OGo65duxIaGsqMGTNYvHgxXl5eJCYmkpubS05OTr373X777SQmJnLrrbcSFhbG559/zpo1a6iqqqpze8MwiI+PZ+DAgXTv3t2yPCYmhnvvvZfQ0FAyMzOZPXs2Q4YMYceOHbi5uZ03Tnl5uSXQARQWFgI1v/D2/Et/tnZ77gEcow9H6AEcow9H6AHUR2NK/+F7lm9dz/dBAXzvcQ1lobcDYDKq6FGxl5uOZHFPz1vp3qv2JSS/rF1zYX2NXbvJMAzjsnc2mfj444+Ji4uzLNuxYwfjxo1j165dODk5ERUVRbNmNdd913f48fjx4zz66KP85z//wWQyERYWRlRUFEuXLqW0tPS87SdPnsx///tfNm/eTPv27eutLycnh9DQUFauXMnIkSPPW5+QkMDzzz9/3vLly5fj6el5sfZFRMTOlBQXsKcoi13B7Uj16EqZycOyzrf6GP1O7KFnbgnBHbpfYBSxNaWlpYwaNYqCggJatGjxq8dr9Fv5e/fuTUpKCgUFBVRUVODr60tERAR9+vSpdx9fX19Wr15NWVkZ+fn5BAUFMX36dDp27Hjetk8++ST//ve/2bRp0wWDEUBgYCChoaH13i03Y8YM4uPjLZ8LCwsJDg5m8ODBtGnT5hI7tj1ms5kNGzYwdOhQu75I0BH6cIQewDH6cIQeQH1crpUr3+Bzj2q+DezGqaDrLct9qvPpVbSfG3JPMf7ex2jeYvAFRqlNc2E78vPzG3W8K/aco5YtWwI1F2knJyfzxz/+8aL7uLu7065dO8xmMx9++CG//e1vLesMw+DJJ5/k448/Jikpqc7g9Ev5+fkcPnyYwMDAOte7ubnVebrNxcXFbn9BzqU+bIcj9ACO0Ycj9ADq42LOlJayfOWbfOdeRVrbdvzQ/mbLuuZGEX0Ld3Nj7mkmjJpE8xa3/arvpbmwvsauu8HhqLi4mAMHDlg+Z2ZmkpKSgo+PDyEhIXzwwQf4+voSEhJCamoqU6dOJS4ujujoaMs+o0ePpl27dsyfPx+Ab7/9luzsbHr16kV2djYJCQlUV1fz7LPPWvaZPHkyy5cvZ82aNXh7e5ObmwvUhDAPDw+Ki4tJSEjg7rvvJjAwkB9//JHnnnuOtm3bctddd132D0hEROxDcWEBC1b8ndSANqQ178iJcy6sNhnVXGveR8SRH3k88i46hE20YqVi6xocjpKTkxk8+OfDjmdPS40ZM4a33nqLnJwc4uPjycvLIzAwkNGjRzN79uxaY2RlZVmuQwIoKytj1qxZZGRk0Lx5c4YPH84777xDq1atLNu89tprAERGRtYaa+nSpYwdOxYnJydSU1NZtmwZp0+fJjAwkMGDB7Nq1Sq8vb0b2qaIiNiJZcv+zpee1WxrfS35XX5+MvXZO826njzKwCovfvfg41asUuxJg8NRZGQkF7qGe8qUKUyZMuWCYyQlJdX6PGjQIPbs2XPBfS523biHhwfr1q274DYiIuIYfr71vhPpwT+fMvM0SrihZA/dck4wsvcgevW+34pVir3Su9VERMRuvLr4ZTYHtiS5eQ/OWG69rzll1jf7EGP7DaPrdY9auUqxdwpHIiJi01J2bOGttK/ZEtSFrC4/X7/atvoY/fL3MsxowT33PWzFCsXRKByJiIhNWvTGn0kK8GZb8x6UBQ8Faq4j6lWWRv+so0wdNZnmLaIvMopIwykciYiITUl8fT5rO4SyOzzKssyvOo+bj+3h7tYdiRo+1nrFSZOgcCQiIjbhz6/P538dO5B2Tc3LzU1GFb3K0xh46MhPR4lut3KF0lQoHImIiNV8/K9lbDQKSPHtwMGfQlEzo4q+pSmMyC/h4TEXvvtZ5EpQOBIRkatq61ef887hXez07UBGi96W5c2MKiJKd3LniTM8PPZJK1YoTZ3CkYiIXBWbk9ax9Pg+vmxzI6WBQ4Ca2/DDKjPoefwQUa5tuOve8VauUkThSERErrAv1q/h3cJDJPn0prTtrQAEVh/l5tw9DHZty933jrVugSK/oHAkIiKN7r9rVrK+8DDf+wWz36UzVW1CAQiqyiY6K5Wnho9ha7GZ4cOHX2QkkatP4UhERBrFti1fsPJAMskBHUn37oLRoqtlXUhVFrcdSmPmvZNoHnUHZrPZipWKXJjCkYiIXLbiwgIWrPg729oHkuJxHRXBPz+UsWNlJtef/JF+ZU4/3XV2p/UKFWkAhSMREWmwHVs38c+D29kU0IPjXX4+NRZYfZS+x/dxGy347f3jrFihyOVTOBIRkUu2cuWb/MfNzDcte1IadBsAnkYxvYvT6H/0JBNHT8XDU9cRiX1TOBIRkQs6U1rKK+8s4OuQ9uzyu5FqkxMA/tW5ROakMqHnbVw75HErVynSeBSORESkTpuT1rHiaBqb/a8j75xTZ9eY9zPo0AGevW8yzVsMs2KFIleGwpGIiNTy3zUrWVWVz6ZWN1D208Ma3Y1SIopSGXy8mAnjn7ZyhSJXlsKRiIgA8P7KJXzkYeZr7xsxm2puw29XdYSBOXsZHR5Bb506kyZC4UhEpIn7y+v/x5fBvuzw60mVqebPQsfKTKIyf+C5B6fi4fkbK1cocnUpHImINEHH846S+MnbbArpwsFrbrcs72Lez+0ZB4gfPQ2PoXdZsUIR61E4EhFpQrZ+9TlvH/meTX7dye8UA4CTUUnvM98z+PBxnprwBytXKGJ9CkciIk3AO8sWsbaFiW9a9KIsYDAA3kYhA0/u4q5mrblz5CNWrlDEdigciYg4sKVv/Y3Vvl5sax+B8dPziYKqsok8msbEfncQPuRJK1coYnsUjkREHNCK997ggxbN2BoywPLQxusq9hL5Ywa/f+gpPDzvsHKFIrZL4UhExEGczD/Bax+8wfb2/mwPvNFy51n38j385tBhpj2u64lELoXCkYiIHTtTWspry/5KckBLdrS4loJrfn6S9bUVPzA8M4PfT3jOihWK2B+FIxERO7T9m02sPLSTLYHXcOSan1/h4WUUcWPxXvoePaFQJHKZFI5EROzIm2+9yhf+Ldle6U95+ygAXIwKepbt4YYjOTx2+30Eh95i5SpF7JvCkYiIHfj4g7dY5lHJ1rAhGKZmAARWH2VA3l5+G3Atg4aPtW6BIg5E4UhExIal7NjCq5k7+LzNTVSY3AG4oex7bjl0hKdGT8PDc/hFRhCRhlI4EhGxQUvf/iubvExs8elBQZua02RhlQe56+B+wgOvY/i4B3BxcbFylSKOSeFIRMRGbPxiLR8dTSPZP4yDIbdalvtV5xGb+R2zRk3FOfJO1q5da8UqRRyfwpGIiJX97Y2X+TLQh2SvHlS0GwqAyaiiW8U++hzNYvJtvyX4tpqXw5rNZmuWKtIkNGvoDps2bSI2NpagoCBMJhOrV6+utT4vL4+xY8cSFBSEp6cnw4YNIz09/YJjms1mXnjhBcLCwnB3d6dnz558+umn5223aNEiOnbsiLu7O7179+arr76qtd4wDBISEggKCsLDw4PIyEjS0tIa2qKIyBV3Mv8E0//5Ejdv+JC54dFsad6HCpMb7aqOcFfuF7xRuI/Pho3ipUemExzaydrlijQpDQ5HJSUl9OzZk4ULF563zjAM4uLiyMjIYM2aNezcuZPQ0FCioqIoKSmpd8xZs2axePFi/va3v7Fnzx4mTJjAXXfdxc6dOy3brFq1imnTpjFz5kx27tzJLbfcQkxMDFlZWZZtXn75ZRITE1m4cCHbt28nICCAoUOHUlRU1NA2RUSumIWL/0TMd5t5q+MwDjqH4WRU0udMCk/tW8vmAUN47YF4YuNGWbtMkSarwafVYmJiiImJqXNdeno6W7duZffu3XTr1g2oOdrj5+fHihUrGD9+fJ37vfPOO8ycOZPhw2vuupg4cSLr1q3jlVde4d133wUgMTGRcePGWcZYsGAB69at47XXXmP+/PkYhsGCBQuYOXMmI0eOBODtt9/G39+f5cuX8/jjjze0VRGRRvXDnhTm/fA1n4UPodrkhKdRQmT+DuKcfLhz5FhrlyciP2nUa47Ky8sBcHd3tyxzcnLC1dWVzZs31xuOysvLa+0D4OHhwebNmwGoqKhgx44dTJ8+vdY20dHRbNmyBYDMzExyc3OJjo62rHdzc2PQoEFs2bKlznBUXl5uqRmgsLAQqDnNZ8/n9c/Wbs89gGP04Qg9gGP0Yc0ezpSWsmDlQv7VuQ/HWt8MwI1luxhfYiI2bmKD6nKEuQDH6MMRegDH6KOxa2/UcNS1a1dCQ0OZMWMGixcvxsvLi8TERHJzc8nJyal3v9tvv53ExERuvfVWwsLC+Pzzz1mzZg1VVVUAnDhxgqqqKvz9/Wvt5+/vT25uLoDlf+va5tChQ3V+3/nz5/P888+ft/zLL7/E09Pz0hu3URs2bLB2CY3CEfpwhB7AMfq4mj2cOVNKasFBvujUjcwudwDQ2jjJyPQt9A3oAS5c9p1njjAX4Bh9OEIPYN99lJaWNup4jRqOXFxc+PDDDxk3bhw+Pj44OTkRFRVV72m4s1599VUeffRRunbtislkIiwsjIcffpilS5fW2s5kMtX6bBjGecsuZZuzZsyYQXx8vOVzYWEhwcHBDB48mDZt2ly0X1tlNpvZsGEDQ4cOtevnoDhCH47QAzhGH1ezh+N5R1mwYQWfBffgqP9vgJpXfAws3MGTPl246ZE/XPbYjjAX4Bh9OEIP4Bh95OfnN+p4jX4rf+/evUlJSaGgoICKigp8fX2JiIigT58+9e7j6+vL6tWrKSsrIz8/n6CgIKZPn07Hjh0BaNu2LU5OTpajQ2cdO3bMcqQoICAAqDmCFBgYWOc2v+Tm5oabm9t5y11cXOz2F+Rc6sN2OEIP4Bh9XMkecrIP8+fPVrChfU+Odah5GayHUcotp3dyn1Mb7vjpFFpjcIS5AMfowxF6APvuo7HrbvDdapeqZcuW+Pr6kp6eTnJyMiNGjLjoPu7u7rRr147Kyko+/PBDyz6urq707t37vEN+GzZsYMCAAQB07NiRgICAWttUVFSwceNGyzYiIlfC4UMZPLXsTwz9YS/vhURzrJk/3kYhsceSWOFyhmUjJ3PHiPutXaaIXKIGHzkqLi7mwIEDls+ZmZmkpKTg4+NDSEgIH3zwAb6+voSEhJCamsrUqVOJi4urdaH06NGjadeuHfPnzwfg22+/JTs7m169epGdnU1CQgLV1dU8++yzln3i4+N56KGH6NOnD/379+eNN94gKyuLCRMmADWn06ZNm8a8efMIDw8nPDycefPm4enpyahRuiVWRBrfti1fsPRQChv9enIyuObhjS2MAoYc+44nug6g+5Bp1i1QRC5Lg8NRcnIygwcPtnw+e83OmDFjeOutt8jJySE+Pp68vDwCAwMZPXo0s2fPrjVGVlYWzZr9fNCqrKyMWbNmkZGRQfPmzRk+fDjvvPMOrVq1smxz3333kZ+fzwsvvEBOTg7du3dn7dq1hIaGWrZ59tlnOXPmDJMmTeLUqVNERESwfv16vL29G9qmiEi93l+5hDVu5Wxp2YszAUMAaGGc5ra875jSPZJrhzxl5QpF5NdocDiKjIzEMIx610+ZMoUpU6ZccIykpKRanwcNGsSePXsu+r0nTZrEpEmT6l1vMplISEggISHhomOJiDTUjq2beCU3jY1+N1FlqvnPZ2D1USKz05jUbzjhQ+IvMoKI2AO9W01E5CLOlJby4vJX+ahTBKda9gegi3k/g388wO/vn0zz24ZbuUIRaUwKRyIiF/DOskW87deK3WE1jyTxrT7GPQeTmfPYc1auTESuFIUjEZE6vPnPv7DWrwXb2velyuSMk1FJ9Mlv+MN1t9D1NgUjEUemcCQico4Fi/+PdSGBpHQYhGGquXHk2oofGJWdw6OP6EJrkaZA4UhEhJrTZ8t9m7Ozy+2WZT3LdjP00GGemTDDipWJyNWmcCQiTdq2LV+wMHcvX7TvS6XJBZNRxU1ndhFz9CQTH33G2uWJiBUoHIlIk5STfZj/+3wla9v3pbD1zUDN6bPfZh1RKBJp4hSORKRJOZl/glf+t4xPQ3px7KenWvtX5zIi4ztm/G4aHp6eVq5QRKxN4UhEmoQzpaVsPZbG87udOdqh5roib6OA6NwdTL85juDbhlm5QhGxFQpHIuLwXnvzzywPDSW9c83DGj2NYm47kczksL700lOtReQXFI5ExGGl7NjCn7JSSAqLpMrkjKtRxpCT2xnftgsDfzvN2uWJiI1SOBIRh3OmtJQXVrzK6g59OdVqAAC9ylIZvv8QEyf9ARcXFytXKCK2TOFIRBzK39/4E+93CGVfp5rXfbStPsbdB5OZ9fDvWVux1srViYg9UDgSEYfw2aerWVyRx+bOQzBMTrgYFUSd/Nbyug+z2WztEkXETigciYhdO5l/goT/vcUnQf0odesAwI1ndvG7k2X87sEnrVuciNglhSMRsVuJr7/EqrBrOdQuCoDgqsOMPJDKjAl6MayIXD6FIxGxOxu/WMvC4iNs7hKNYWqGp1HCHUe/4fmYsfhExVq7PBGxcwpHImI30ven8Zcd61jvfxPF3n0BiCj9jscqPbnjQb3yQ0Qah8KRiNi8k/knmPvJP/m0/Y3kBwwBIKgqm3sP7NIpNBFpdApHImKzzpSWMu+9Bfy3Y0+OhkQD0No4ScyRZGbc/iC+UXdYuUIRcUQKRyJik15ePJ81nbpw8KdXfngZxQw9tp1p1w+m65BnrVydiDgyhSMRsSlv/vMvfBjYlpQuNQ9xdDXKiTydzGM+nRl4/1NWrk5EmgKFIxGxCWdKS3lm9eus7nArVSZnTEYVA0q+48EyF+66d7K1yxORJkThSESsbt1//8WfnCrYHVhzsXX38j3ceySXx8fHW7kyEWmKFI5ExKrmvj6PZV0GUGBqhYtRwV1HN/F/Iyfh4elp7dJEpIlSOBIRq9j4xVoWFR1mU5dhGKZmBFYf5dGDe5j0mJ5XJCLWpXAkIlfVyfwTvPDfpXzSvh/FLSIAGFi8nT926M21tw23cnUiIgpHInIVvbx4Pv8K60ZW8FCg5kGOdx/YxUw9yFFEbIjCkYhccT8e/IGZ33/BF+G3W96FFpOzlVmRvyVQD3IUERujcCQiV9Tf3/gT/+x0LdmtBgDQvySZCaZW3P67p61cmYhI3RSOROSKKC4s4Jm1/+STzpFUmlzwNgp44MAWXnhshrVLExG5IIUjEWl0K1e+ySKfluz3HwzA9eW7mVbpynAFIxGxA80ausOmTZuIjY0lKCgIk8nE6tWra63Py8tj7NixBAUF4enpybBhw0hPT7/ouAsWLOCaa67Bw8OD4OBgnnrqKcrKyizrO3TogMlkOu9r8uSfn5w7duzY89b369evoS2KyGU6U1rKs0v/j+f8rmO/SxfcjDLuP7yBNbeOZPhvfmvt8kRELkmDjxyVlJTQs2dPHn74Ye6+++5a6wzDIC4uDhcXF9asWUOLFi1ITEwkKiqKPXv24OXlVeeY7733HtOnT+ef//wnAwYMYP/+/YwdOxaAv/zlLwBs376dqqoqyz67d+9m6NCh3HvvvbXGGjZsGEuXLrV8dnV1bWiLInIZkj7/hD+V57Ojw+0AdKrM4NHsbB4e+3srVyYi0jANDkcxMTHExMTUuS49PZ2tW7eye/duunXrBsCiRYvw8/NjxYoVjB8/vs79vvnmG26++WZGjRoF1BwleuCBB9i2bZtlG19f31r7vPTSS4SFhTFo0KBay93c3AgICGhoWyLyKzy/eB4rw/txyqM9zYwqhud/zfxb78F3aJC1SxMRabBGveaovLwcAHd3d8syJycnXF1d2bx5c73haODAgbz77rts27aNvn37kpGRwdq1axkzZkyd21dUVPDuu+8SHx+PyWSqtS4pKQk/Pz9atWrFoEGDmDt3Ln5+fvXWe7ZmgMLCQgDMZjNms/nSG7cxZ2u35x7AMfpwhB6g/j6+S/6aP51IZ3OXmoc3+lfnMiZ9J0+Oe6bO7a3J0efC3jhCH47QAzhGH41du8kwDOOydzaZ+Pjjj4mLiwNqigsPD6dv374sXrwYLy8vEhMTmTFjBtHR0axbt67esf72t7/x9NNPYxgGlZWVTJw4kUWLFtW57fvvv8+oUaPIysoiKOjn/2e6atUqmjdvTmhoKJmZmcyePZvKykp27NiBm5vbeeMkJCTw/PPPn7d8+fLleOq9TiIXlHJ4Fx9c258TzWqO6kYWbGVYvonWbXXkVkSurtLSUkaNGkVBQQEtWrT41eM1ajgC2LFjB+PGjWPXrl04OTkRFRVFs2Y1132vXbu2znGSkpK4//77efHFF4mIiODAgQNMnTqVRx99lNmzZ5+3/e23346rqyv/+c9/LlhfTk4OoaGhrFy5kpEjR563vq4jR8HBweTk5NCmTZtL+RHYJLPZzIYNGxg6dCguLi7WLueyOUIfjtAD1O5j57avWFD4I19517z6o231MUalJ/P7cbZ9bZEjzoX6sC5H6AEco4/8/HwCAwMbLRw1+q38vXv3JiUlhYKCAioqKvD19SUiIoI+ffrUu8/s2bN56KGHLKfdevToQUlJCY899hgzZ860hCuAQ4cO8dlnn/HRRx9dtJbAwEBCQ0PrvVvOzc2tziNKLi4udvsLci71YTscoYeKsgpeXP4qH3Xqx+mfgtEtRduY2a47vezo9R+OMBegPmyJI/QA9t1HY9d9xZ5z1LJlS6DmIu3k5GT++Mc/1rttaWlprQAENdcqGYbBLw9sLV26FD8/P+644+KvHMjPz+fw4cMEBgZeRgcictbH/1rGP1q78b1fzc0Y/tW5/Db9O70TTUQcUoPDUXFxMQcOHLB8zszMJCUlBR8fH0JCQvjggw/w9fUlJCSE1NRUpk6dSlxcHNHR0ZZ9Ro8eTbt27Zg/fz4AsbGxJCYmcsMNN1hOq82ePZs777wTJycny37V1dUsXbqUMWPG4Oxcu/Ti4mISEhK4++67CQwM5Mcff+S5556jbdu23HXXXQ3+wYhIzXOLZv5rIR+1H0iZyRMno5Lok9/w3PVDCL9tmLXLExG5IhocjpKTkxk8eLDlc3x8PABjxozhrbfeIicnh/j4ePLy8ggMDGT06NHnXTeUlZVV60jRrFmzMJlMzJo1i+zsbHx9fYmNjWXu3Lm19vvss8/IysrikUceOa8uJycnUlNTWbZsGadPnyYwMJDBgwezatUqvL29G9qmSJO39pP3WeBcwffBNf/HJqwyg4cOZTJhvN6JJiKOrcHhKDIy8rxTXeeaMmUKU6ZMueAYSUlJtYtwdmbOnDnMmTPngvtFR0fX+709PDwueDeciFy62W/OY2XYzRSZWuJsmLnj2GYG48/dYy7877aIiCNo8OtDRMRx7U39jrv/8yZvdh5Okakl7aqOMONAEgvvfgIPj7qfcC8i4mj04lkRAeC1N//Mmx2v5WjzmzAZ1USd/oY/9ryNDlG/seuHw4mINJTCkUgTd6a0lBn/+jsfhg3CbHKlpXGah9K3MOtx3YkmIk2TwpFIE7Ztyxc8f/owO4KHAtCtYi9Pm50YrmAkIk2YwpFIE/X6P17htY7dyPPoiZNRSeyxr/jz8Edo3qKltUsTEbEqhSORJujZpS+xslMkFSZ3WlfnMy59G89MmGHtskREbILCkUgTcvhQBs+kbGBjh5oHOF5b8QPTK525XcFIRMRC4UikiXjtzT+ztGMXslrUvBctJv8rXhn8AD5t2lq5MhER26JwJOLgTuaf4LkN7/BJ2CAqTS54GUU8eOArnn9MF12LiNRF4UjEgS1b9nfeCPDngH/NK396lKfxZKnBnQpGIiL1UjgScVAz/zGP9zoNpszkgbtxhpFHvmLuPU/g4elp7dJERGyawpGIgzlTWsqT//0Hn4QNByDcnM6EYyf53ehnrVyZiIh9UDgScSC7U77l2Zwf+K7trQBEFmzlb/3vxNc/yMqViYjYD4UjEQex4r03SPQP5LB7zUMd7zvyBYk6WiQi0mAKRyIOYM4b83i38y2UmLxpYRQwbv/X/GGCLroWEbkcCkcidux43lGe/vpj1ofXXF8UWvkjT+Uf534FIxGRy6ZwJGKnVn/4Dq94u5Pe+mYAIgu38lL3SDoMjbNuYSIidk7hSMQOvbx4Pv8IH0ChqSXuRikPZG5k/ji9AkREpDEoHInYmWeX/h/Lw6OoNLnQvuowk45k8IiCkYhIo1E4ErETxYUFPLHhHT7tcDsAvcpS+XNgF7pHxVq5MhERx6JwJGIH9qZ+x9NHUvnOZyAA0ae+ZlHUgzRv0dLKlYmIOB6FIxEbl/T5JzxXXUHGT88vuv/wF7wyRs8vEhG5UppZuwARqd+/P3qXpw3IcO6El1HMk/vXKxiJiFxhOnIkYqPeWbaIl9t15ngzP1oZp5h6cAcT9fwiEZErTuFIxAa99uafeTWsN6dNrfGtPsaz2Qd46NFnrF2WiEiToHAkYmOef2MeyzrfSompOe2qjjCn8DR3jp5k7bJERJoMhSMRG3Ey/wTPfLmc/3UehmFqRqfKDOY1cyVy5IPWLk1EpElROBKxAev++y/+z9nMnja3AnBL0Tbmdb2Z8C7drFyZiEjTo3AkYmWvvD6ff4T35VSzNrga5dx3KIk/PfwHa5clItJkKRyJWNHMf8xjWZcozCZX/KtzmZCxm4mPKhiJiFiTwpGIlTzx3it82Knm+qIe5Wm81CqY3rojTUTE6hSORK6y4sICHv/8PT4Pug2oub5oUcRv8PUPsnJlIiICCkciV9WPB39g8r5v2dFqAAB3Hkvi1djH8PD0tHJlIiJyVoNfH7Jp0yZiY2MJCgrCZDKxevXqWuvz8vIYO3YsQUFBeHp6MmzYMNLT0y867oIFC7jmmmvw8PAgODiYp556irKyMsv6hIQETCZTra+AgIBaYxiGQUJCAkFBQXh4eBAZGUlaWlpDWxS5ItL3pzE+/Tt2eNS8I230j5/yxn3TFIxERGxMg8NRSUkJPXv2ZOHCheetMwyDuLg4MjIyWLNmDTt37iQ0NJSoqChKSkrqHfO9995j+vTpzJkzh71797JkyRJWrVrFjBkzam3XrVs3cnJyLF+pqam11r/88sskJiaycOFCtm/fTkBAAEOHDqWoqKihbYo0qh927+SxzFR2u12Hq1HOxPT1vPzwdGuXJSIidWjwabWYmBhiYmLqXJeens7WrVvZvXs33brVPJ9l0aJF+Pn5sWLFCsaPH1/nft988w0333wzo0aNAqBDhw488MADbNu2rXaxzs7nHS06yzAMFixYwMyZMxk5ciQAb7/9Nv7+/ixfvpzHH3+8oa2KNIpTeUeY7FHOfteuuBllTNj/BTP0jjQREZvVqNcclZeXA+Du7m5Z5uTkhKurK5s3b643HA0cOJB3332Xbdu20bdvXzIyMli7di1jxoyptV16ejpBQUG4ubkRERHBvHnz6NSpEwCZmZnk5uYSHR1t2d7NzY1BgwaxZcuWOsNReXm5pWaAwsJCAMxmM2az+TJ/CtZ3tnZ77gEco49vv/6Sf3Tw54BLZ9yNUibuSyJ+/O/tridHmAtH6AHUhy1xhB7AMfpo7NobNRx17dqV0NBQZsyYweLFi/Hy8iIxMZHc3FxycnLq3e/+++/n+PHjDBw4EMMwqKysZOLEiUyf/vNph4iICJYtW0aXLl3Iy8vjxRdfZMCAAaSlpdGmTRtyc3MB8Pf3rzW2v78/hw4dqvP7zp8/n+eff/685V9++SWeDnAdyIYNG6xdQqOw1z7yc7NY0qkdB10642mUMC51A1073MDatWutXdpls9e5OJcj9ADqw5Y4Qg9g332UlpY26niNGo5cXFz48MMPGTduHD4+Pjg5OREVFVXvabizkpKSmDt3LosWLSIiIoIDBw4wdepUAgMDmT17NkCtMXr06EH//v0JCwvj7bffJj4+3rLOZDLVGtswjPOWnTVjxoxa+xYWFhIcHMzgwYNp06ZNg/u3FWazmQ0bNjB06FBcXFysXc5ls+c+dqd8y589zBx07oSnUcyT+zfzxKSZ1i7rstnzXJzlCD2A+rAljtADOEYf+fn5jTpeo9/K37t3b1JSUigoKKCiogJfX18iIiLo06dPvfvMnj2bhx56yHLarUePHpSUlPDYY48xc+ZMmjU7/7pxLy8vevToYbkT7uy1SLm5uQQGBlq2O3bs2HlHk85yc3PDzc3tvOUuLi52+wtyLvVhHXtTv2PaqaOku3TBwyhlfOpnPDFppl31UB97m4u6OEIPoD5siSP0APbdR2PX3eC71S5Vy5Yt8fX1JT09neTkZEaMGFHvtqWlpecFICcnJwzDwDCMOvcpLy9n7969liDUsWNHAgICah0WrKioYOPGjQwYMKAROhK5uPT9aUw6up99Ll1wN0qZvH8j3TrcYO2yRESkARp85Ki4uJgDBw5YPmdmZpKSkoKPjw8hISF88MEH+Pr6EhISQmpqKlOnTiUuLq7WhdKjR4+mXbt2zJ8/H4DY2FgSExO54YYbLKfVZs+ezZ133omTkxMAzzzzDLGxsYSEhHDs2DFefPFFCgsLLRdtm0wmpk2bxrx58wgPDyc8PJx58+bh6elpuQtO5EpK35/GxIxd7HW7DjejjIn7NzJ13DN2fY2RiEhT1OBwlJyczODBgy2fz16zM2bMGN566y1ycnKIj48nLy+PwMBARo8ebblu6KysrKxaR4pmzZqFyWRi1qxZZGdn4+vrS2xsLHPnzrVsc+TIER544AFOnDiBr68v/fr1Y+vWrYSGhlq2efbZZzlz5gyTJk3i1KlTREREsH79ery9vRvapkiDJH3+CbOqyjjw03OMJuz/gj9MeM6u7/4QEWmqGhyOIiMj6z3VBTBlyhSmTJlywTGSkpJqF+HszJw5c5gzZ069+6xcufKitZlMJhISEkhISLjotiKNZcV7b/ByQHtyfror7dH9G/UcIxERO6Z3q4n8Cn9742X+3vkmTpta09o4yZMHvmOSgpGIiF1TOBK5THNfn8ebXSIpM3nSruoIM/KPc89jz1i7LBER+ZUUjkQuwx8Xz+ONLlGYTa50Me9nnmtLBt73sLXLEhGRRqBwJNJALy6exxvhNcHo+vLdvBHehw5hXa1dloiINBKFI5EGmPv6PBZ3+TkYLbl2AMGhnaxdloiINCKFI5FL9OLiebzR5TbMJld6lKfxRngfBSMREQd0xZ6QLeJIXnp9Hm+E30aFyY0e5Wm8Gd5bp9JERByUwpHIRSx44/94vUskFSY3upfvUTASEXFwCkciF7D4H4ks7DyAMpMnXczpLAy5TsFIRMTBKRyJ1OO9d17nL516UWzypkPlj/y5eRu6XtfL2mWJiMgVpnAkUofVH77DS0GdOG3yIagqm3mGib4Dhli7LBERuQoUjkR+4bNPV/NCSx+ON/PDt/oYCYWnGBI9wtpliYjIVaJwJHKOb7/6guecTRx1akcr4yTTj2Zw58gHrV2WiIhcRQpHIj9JS/2Op8tOkeUUirdRSHzmLn730ARrlyUiIleZwpEI8OPBH3jy6D4OOIfhaZTwZPpWHhv3lLXLEhERK1A4kibveN5RJuzfzh7Xa3EzypiQvokpjz9r7bJERMRKFI6kSTtTWsrj29aS4t4DZ8PMY/u/4NnHZ1i7LBERsSKFI2nSHv90KVua98FkVDH24AZmTnjO2iWJiIiVKRxJkzV5eSLrW98MwANHPufFRxWMRERE4UiaqGeXvsSHgTUPdYw9lkTiaF1jJCIiNRSOpMl5YfE83guNAiCyYCt/jX3MyhWJiIgtUTiSJuWvb/yJf4YPpsrkzE2lO/nHkAfw8PS0dlkiImJDFI6kyVi16h8s6tyHMpMHXSv2sfC6m2neoqW1yxIRERujcCRNwjdffcb/tQnktKk17aqOMNetNaEdu1i7LBERsUEKR+Lwso9kMeNMvuV9aX84lsPNkdHWLktERGyUwpE4tDOlpUza9Rk/uFyDu3GGyfu389tRj1q7LBERsWEKR+LQnvzkTb71vJFmRhWPpH/JkxP+YO2SRETExikcicN6bsl8PvEdBMBvj3zB/3tcD3kUEZGLUzgSh/TqGy/zbseahzwOKfiGBaN/b+WKRETEXigcicP590fv8FrnPlSY3OhensZf+4+wdkkiImJHFI7EofywZxdzvb05bfIhqCqbFzz8aOsfZO2yRETEjigcicM4U1rKMz/u5JBzB5obRUzJzmTAoKHWLktEROxMg8PRpk2biI2NJSgoCJPJxOrVq2utz8vLY+zYsQQFBeHp6cmwYcNIT0+/6LgLFizgmmuuwcPDg+DgYJ566inKysos6+fPn89NN92Et7c3fn5+xMXFsW/fvlpjjB07FpPJVOurX79+DW1R7NRT/15MskcvnIxKxu3bxNgxT1i7JBERsUMNDkclJSX07NmThQsXnrfOMAzi4uLIyMhgzZo17Ny5k9DQUKKioigpKal3zPfee4/p06czZ84c9u7dy5IlS1i1ahUzZsywbLNx40YmT57M1q1b2bBhA5WVlURHR5837rBhw8jJybF8rV27tqEtih168Y15rPG7FYCRR5OYMXGmlSsSERF75dzQHWJiYoiJialzXXp6Olu3bmX37t1069YNgEWLFuHn58eKFSsYP358nft988033HzzzYwaNQqADh068MADD7Bt2zbLNp9++mmtfZYuXYqfnx87duzg1ltvtSx3c3MjICCgoW2JHXv//bd4q/MtGCYn+pck8/LISdYuSURE7FiDw9GFlJeXA+Du7m5Z5uTkhKurK5s3b643HA0cOJB3332Xbdu20bdvXzIyMli7di1jxoyp93sVFBQA4OPjU2t5UlISfn5+tGrVikGDBjF37lz8/PzqrfdszQCFhYUAmM1mzGbzJXRsm87Wbs89wKX1kZGxn1dat6bY5E1I1SGeb9cdZxcXm+m9Kc2FrXOEHkB92BJH6AEco4/Grt1kGIZx2TubTHz88cfExcUBNcWFh4fTt29fFi9ejJeXF4mJicyYMYPo6GjWrVtX71h/+9vfePrppzEMg8rKSiZOnMiiRYvq3NYwDEaMGMGpU6f46quvLMtXrVpF8+bNCQ0NJTMzk9mzZ1NZWcmOHTtwc3M7b5yEhASef/7585YvX74cT0/PBv405Gqrrq5mmctJvvHqTXOjiCd2JtGx8/XWLktERK6y0tJSRo0aRUFBAS1atPjV4zVqOALYsWMH48aNY9euXTg5OREVFUWzZjWXNtV3/U9SUhL3338/L774IhERERw4cICpU6fy6KOPMnv27PO2nzx5Mv/973/ZvHkz7du3r7e+nJwcQkNDWblyJSNHjjxvfV1HjoKDg8nJyaFNmzaX+mOwOWazmQ0bNjB06FBcXFysXc5lu1gfz6x8lVXtozAZVUzc9ykzxj9rhSovrKnMhT1whB5AfdgSR+gBHKOP/Px8AgMDGy0cNeppNYDevXuTkpJCQUEBFRUV+Pr6EhERQZ8+ferdZ/bs2Tz00EOW0249evSgpKSExx57jJkzZ1rCFcCTTz7Jv//9bzZt2nTBYAQQGBhIaGhovXfLubm51XlEycXFxW5/Qc7lyH0kvvky/wqreQL2ncc28f9s/AJsR54Le+MIPYD6sCWO0APYdx+NXfcVe85Ry5Yt8fX1JT09neTkZEaMqP8pxaWlpbUCENRcq2QYBmcPbBmGwRNPPMFHH33EF198QceOHS9aQ35+PocPHyYwMPDXNSM25cvPPuHNTr2pMjlzQ9n3JA4fZ+2SRETEgTT4yFFxcTEHDhywfM7MzCQlJQUfHx9CQkL44IMP8PX1JSQkhNTUVKZOnUpcXBzR0dGWfUaPHk27du2YP38+ALGxsSQmJnLDDTdYTqvNnj2bO++8EycnJ6DmVNry5ctZs2YN3t7e5ObmAjUhzMPDg+LiYhISErj77rsJDAzkxx9/5LnnnqNt27bcddddv+qHJLbj9OmTPF9dwimX9gRU5zDL3RevRjiEKiIiclaDw1FycjKDBw+2fI6PjwdgzJgxvPXWW+Tk5BAfH09eXh6BgYGMHj36vOuGsrKyah0pmjVrFiaTiVmzZpGdnY2vry+xsbHMnTvXss1rr70GQGRkZK2xli5dytixY3FyciI1NZVly5Zx+vRpAgMDGTx4MKtWrcLb27uhbYqNmvLFCn5ofTNuRhnj07/n5gl/sHZJIiLiYBocjiIjI7nQNdxTpkxhypQpFxwjKSmpdhHOzsyZM4c5c+bUu8/Frhv38PC44N1wYv9mLXmJ9Z2GAfBAZhJPTJhu5YpERMQR6d1qYhfe/+At3utY87DPIQXfMP8RHTESEZErQ+FIbN6J/OP8pVVLzpg8Cas8yLwegzGZTNYuS0REHJTCkdi8329ZTaZzRzyNYh4/fIQOYV2tXZKIiDgwhSOxaTtydvNZ6wEA3H9wE6MfmWrlikRExNEpHInN+vzz/7Ciyy0ADCr8lhd1nZGIiFwFCkdik0pKipnfrJIiUwuCqw7zXLseNPvpmVciIiJXksKR2KSnPvkH+1y64GaUMfaHVHr2GWDtkkREpIlQOBKbs/jdv/Nf34EA3JO1kcce/72VKxIRkaZE4UhsyvETx3gjIJQqkzO9ylPp6x5i7ZJERKSJUTgSm/LUV/8i26k9LYwCHsktwt3D09oliYhIE6NwJDbjlbcT+bxlPwB+u38zd/3uUStXJCIiTZHCkdiE9Ix9LG3fA8PUjAElyfy/h5+xdkkiItJEKRyJTfhD2lecaOZL2+rjPFrujqubm7VLEhGRJkrhSKxuztuvsKV5H0xGNQ+kfUPM3Q9auyQREWnCFI7EqrZ99w3Lg28CIKrgG6ZPes7KFYmISFOncCRWU1VVxaxjGRSZWtC+6jCTvENwcna2dlkiItLEKRyJ1fx+5QK+d+uGs2FmVEoy/aNirV2SiIiIwpFYx/82ruXDwJsBiD22maeeTrBuQSIiIj9ROJKrrsxczosVZZSb3Oli3s+kTjdhMpmsXZaIiAigcCRW8NRHr3PQuRMeRin3p6TRo+9Aa5ckIiJioXAkV9VnWz7nE9/+AIw8vJlJzz5v5YpERERqUziSq6aqqooXi05gNrnS1fwD47vfYu2SREREzqNwJFfNzA8W8oPrNbgYFYxM/p5re/e3dkkiIiLnUTiSq+KHQ+l84HcjAMNPfM2k3+t0moiI2CaFI7kqfr/7a0pM3oRUHWJEuQfOLi7WLklERKROCkdyxf1l9T/Z7tkLk1HFPd/tYPhDE6xdkoiISL0UjuSKOlF0in94hwBwW+FWxox+1MoViYiIXJjCkVxRUz9fRX6ztrStPk5UWjb+AcHWLklEROSCFI7kivnP9s/5osVNANyz71vGPDfXyhWJiIhcnMKRXBEVVZW8eLoYw+REn7KdDO3cVa8IERERu6BwJFfE9E/+ySHnUDyNYmK+SuPm6HusXZKIiMglUTiSRrc7ex8fel8PwJ1Hv+Hep6ZZtyAREZEGaHA42rRpE7GxsQQFBWEymVi9enWt9Xl5eYwdO5agoCA8PT0ZNmwY6enpFx13wYIFXHPNNXh4eBAcHMxTTz1FWVlZrW0WLVpEx44dcXd3p3fv3nz11Ve11huGQUJCAkFBQXh4eBAZGUlaWlpDW5RfwTAMnv4+mXKTO50rD3DDwVz8fNtbuywREZFL1uBwVFJSQs+ePVm4cOF56wzDIC4ujoyMDNasWcPOnTsJDQ0lKiqKkpKSesd87733mD59OnPmzGHv3r0sWbKEVatWMWPGDMs2q1atYtq0acycOZOdO3dyyy23EBMTQ1ZWlmWbl19+mcTERBYuXMj27dsJCAhg6NChFBUVNbRNuUwLkz5il3s3nAwzI7el8tCsP1m7JBERkQZxbugOMTExxMTE1LkuPT2drVu3snv3brp16wbUHO3x8/NjxYoVjB8/vs79vvnmG26++WZGjRoFQIcOHXjggQfYtm2bZZvExETGjRtnGWPBggWsW7eO1157jfnz52MYBgsWLGDmzJmMHDkSgLfffht/f3+WL1/O448/3tBWpYFOlxXzWnUraAZRBVu5sVc4zZx05lZEROxLg8PRhZSXlwPg7u5uWebk5ISrqyubN2+uNxwNHDiQd999l23bttG3b18yMjJYu3YtY8aMAaCiooIdO3Ywffr0WvtFR0ezZcsWADIzM8nNzSU6Otqy3s3NjUGDBrFly5Y6w1F5ebmlZoDCwkIAzGYzZrP5cn4ENuFs7Ve7h/jP3+ek5420MY4zIGk/N/858VfVYK0+GpMj9ACO0Ycj9ADqw5Y4Qg/gGH00du2NGo66du1KaGgoM2bMYPHixXh5eZGYmEhubi45OTn17nf//fdz/PhxBg4ciGEYVFZWMnHiREsYOnHiBFVVVfj7+9faz9/fn9zcXADL/9a1zaFDh+r8vvPnz+f5589/AeqXX36Jp6fnpTduozZs2HDVvld2dQHrWtZchH33gWQ8evdg7dq1jTL21ezjSnGEHsAx+nCEHkB92BJH6AHsu4/S0tJGHa9Rw5GLiwsffvgh48aNw8fHBycnJ6Kiouo9DXdWUlISc+fOZdGiRURERHDgwAGmTp1KYGAgs2fPtmz3y+fkGIZx3rJL2easGTNmEB8fb/lcWFhIcHAwgwcPpk2bNpfUsy0ym81s2LCBoUOH4nIVXvBqGAZDk/5NlcmZ6ytSCdm7nwdefuNXj3u1+7gSHKEHcIw+HKEHUB+2xBF6AMfoIz8/v1HHa9RwBNC7d29SUlIoKCigoqICX19fIiIi6NOnT737zJ49m4ceeshy2q1Hjx6UlJTw2GOPMXPmTNq2bYuTk5Pl6NBZx44dsxwpCggIAGqOIAUGBta5zS+5ubnh5uZ23nIXFxe7/QU519Xq49Uta9jn0hlXo5w7t2Rwd8Irjfp9HWE+HKEHcIw+HKEHUB+2xBF6APvuo7HrvmJXy7Zs2RJfX1/S09NJTk5mxIgR9W5bWlpKs2a1S3FycsIwDAzDwNXVld69e593yG/Dhg0MGDAAgI4dOxIQEFBrm4qKCjZu3GjZRhpf/pkiFpW1BOD2U9/g7XGc1i19rFyViIjI5WvwkaPi4mIOHDhg+ZyZmUlKSgo+Pj6EhITwwQcf4OvrS0hICKmpqUydOpW4uLhaF0qPHj2adu3aMX/+fABiY2NJTEzkhhtusJxWmz17NnfeeSdOTk4AxMfH89BDD9GnTx/69+/PG2+8QVZWFhMmTABqTqdNmzaNefPmER4eTnh4OPPmzcPT09NyF5w0vqe/WkOBS3f8q3Pov/4ID/wt0doliYiI/CoNDkfJyckMHjzY8vnsNTtjxozhrbfeIicnh/j4ePLy8ggMDGT06NG1rhsCyMrKqnWkaNasWZhMJmbNmkV2dja+vr7ExsYyd+7PLyq97777yM/P54UXXiAnJ4fu3buzdu1aQkNDLds8++yznDlzhkmTJnHq1CkiIiJYv3493t7eDW1TLsG32ftY73wtAPcc2Emrm9vg6upq5apERER+nQaHo8jISAzDqHf9lClTmDJlygXHSEpKql2EszNz5sxhzpw5F9xv0qRJTJo0qd71JpOJhIQEEhISLjiO/HqGYTAjLZVql87cUJFC0NZ07vrnP61dloiIyK+mJ/TJZXln90b2uHTGxaggdvsROo8ZXu9dgSIiIvZE4UgarKyqkj8fqwBgSNE3OGXuZlDkvVauSkREpHEoHEmD/XHrfzjWzI9WxkmGbCzk2rEPWLskERGRRqNwJA1ypPgk75XXPFPqztxtnCn6jltvqf8xDSIiIvZG4Uga5PfbP6fM5EHHqoNc/3kZ140eY+2SREREGpXCkVyyzTn7+ZJwAO7ev5cil2RuvSXOukWJiIg0MoUjuSTVhsEf9uwHoG/5NgK+Kefa+x+2clUiIiKNT+FILslre77iYLP2uBtniNtxnILm2xg0aKS1yxIREWl0CkdyUSfLy1iQV/PPw05vwnl3MV3vG6/nGomIiENSOJKLmvHdZxSZWhBUfYQhX1dR2Go7QwbdY+2yRERErgiFI7mgnSeP8u8zQQCMPPwdpYcLCYt7UEeNRETEYSkcSb0MwyD++10Ypmb0NifT9VtXiltv5fao0dYuTURE5IpROJJ6LctIYa8RiJtRxsjdeZQcP4XvkKE0c9KvjYiIOC5naxcgtqnQXMn8Q0VgasHQ4i/x2muQ3+ob7ov7zNqliYiIXFE6BCB1+n+pmzhtaoG/cZTh212pOFmMZ79euLu5Wrs0ERGRK0pHjuQ8uwtO8v7plmCCETlfU3TEhZIWX3P/PZ9YuzQREZErTkeOpBbDMHjq+++oNjlxQ2Uy/Xb641RQATeE4NOqpbXLExERueJ05EhqWXn4AKmVbXE1yrljfyZ5+S6Uem/iN3HLrV2aiIjIVaEjR2JRVFnFCwdrHoUdfWY916S3w7mwGnOPFoSGhFi5OhERkatDR47E4sW9yZyiOf5GDoO/dyHr1GEMlxQG3T7X2qWJiIhcNTpyJADsKSrinRM1Wfk3J9YTcNQbyltS0OU4N/S40crViYiIXD06ciSYqw0mfr+LappzY9V2+qZdw8GCdJqfOUrrG/Q0bBERaVp05Eh4+cAP7KtoTnOjiKEZO2lxwkx1dTtyQ3ZxW/Sj1i5PRETkqlI4auJ2FBTz9+wzAPy2+H16/DiQfYX7CMhPx7lLD1p6eVi5QhERkatL4agJK62qZmJqGtU0o1/1ZkJT29Ds9DEqm3XhR/+t3DJourVLFBERueoUjpqwhH37yTK70drIJzIjmd6FN/BD4R7an8ihLNyJXj16WbtEERGRq07hqIlKyi9kWV4ZAL8rfo+AQwM4VbCfSudrOem6kfDwUZhMJitXKSIicvUpHDVBmaXlPJG2H4DbqtfhuacN11cE80NhKgGF5eRdm8+g23UhtoiINE0KR01MdlkF9+zcy4kqZ4KNQwz4cRc9im7kUMFOql2uxbXka1oE9aJtS29rlyoiImIVCkdNyPEKM/fu3Ed2BQQYR5l4+nWqs24gtKI1+wt20rwqgL0d0+kT8bS1SxUREbEahaMm4rS5kvtS0skoq6KNcZyppYlkp93AgIpr2Xf6WwzXawjK244p1Iubet9k7XJFRESsRuGoCThaVsHvdh1kT0kFLY1TPGN+idzUXlxX3hlvczMOFqXg7Ho9xzy3E9r+Tpo104XYIiLSdOn1IQ6qGth0qph3806z/kQBVUBzo4g/VM/jROp1tChvR5/KMHaeXo/h3JGA/Ex+CC/ggUHjrF26iIiIVTX4yNGmTZuIjY0lKCgIk8nE6tWra63Py8tj7NixBAUF4enpybBhw0hPT7/gmJGRkZhMpvO+7rjjDss2HTp0qHObyZMnW7YZO3bseev79evX0BbtWnFlFW9kn2COVxC/SzvE/34KRtcau5lhPE9RWgfMZwIYXN6NUvNpMopScXbvQ+sTm3EJ9CUspL21WxAREbGqBh85KikpoWfPnjz88MPcfffdtdYZhkFcXBwuLi6sWbOGFi1akJiYSFRUFHv27MHLy6vOMT/66CMqKiosn/Pz8+nZsyf33nuvZdn27dupqqqyfN69ezdDhw6ttQ3AsGHDWLp0qeWzq6trQ1u0axWGwcuHjlHu5IInZxhofMltrCfYlMe+fb05daodQyq70sLwYFP+J+AcSPNy2B18kPB2OmokIiLS4HAUExNDTExMnevS09PZunUru3fvplu3bgAsWrQIPz8/VqxYwfjx4+vcz8fHp9bnlStX4unpWSv4+Pr61trmpZdeIiwsjEGDBtVa7ubmRkBAQEPbchgelTnEsZrmRg4D+Irmzm60bXMPGzZUcvp0Nd28O9HpuB955YfIOZOBi9cIgrI28+OASu6KHG3t8kVERKyuUa85Ki8vB8Dd3d2yzMnJCVdXVzZv3lxvOPqlJUuWcP/999d7pKmiooJ3332X+Pj4857inJSUhJ+fH61atWLQoEHMnTsXPz+/eus9WzNAYWEhAGazGbPZfEm12honJ3/ud/uG4uLjhIZM5vCRDnz88XdUVlbTpqUPN+WFYGCw4/gGTM1a4eQUDOWLadGqHYFtWtlU32drsaWaGsoRegDH6MMRegD1YUscoQdwjD4au3aTYRjGZe9sMvHxxx8TFxcH1BQXHh5O3759Wbx4MV5eXiQmJjJjxgyio6NZt27dRcfctm0bERERfPvtt/Tt27fObd5//31GjRpFVlYWQUFBluWrVq2iefPmhIaGkpmZyezZs6msrGTHjh24ubmdN05CQgLPP//8ecuXL1+Op6fnJf4UbNEJTp82kZ2dazld2bx5c24rvg7fMm8OlKWyI2ctzh6RtCuoYo/v27S9diTtOukWfhERsT+lpaWMGjWKgoICWrRo8avHa9RwBLBjxw7GjRvHrl27cHJyIioqimbNaq77Xrt27UXHfPzxx9myZQupqan1bnP77bfj6urKf/7znwuOlZOTQ2hoKCtXrmTkyJHnra/ryFFwcDA5OTm0adPmorXaoqKiIv7zn/+QmZkJgLe3N1FRUYQWtqb4v4cwXOHj9AWYjWrcWj7GDSmL+DTuBx595BsCfVpaufrazGYzGzZsYOjQobi4uFi7nMviCD2AY/ThCD2A+rAljtADOEYf+fn5BAYGNlo4avRb+Xv37k1KSgoFBQVUVFTg6+tLREQEffr0uei+paWlrFy5khdeeKHebQ4dOsRnn33GRx99dNHxAgMDCQ0NrfduOTc3tzqPKLm4uNjtL4iXlxfHjx/HZDLRv39/Bg0ahEtVM3L/nAxAlks65upynNx64l5ewuHWB/F170yIf1srV14/e56PsxyhB3CMPhyhB1AftsQRegD77qOx675iD4Fs2bIlvr6+pKenk5yczIgRIy66z/vvv095eTkPPvhgvdssXboUPz+/Wrf51yc/P5/Dhw8TGBjYoNrtmZubGyNGjKBr164MHjwY58pmHF+aRnVpJc3auLItbQ0ATm434HsihR3XGYSH/c7KVYuIiNiOBoej4uJiUlJSSElJASAzM5OUlBSysrIA+OCDD0hKSiIjI4M1a9YwdOhQ4uLiiI6OtowxevRoZsyYcd7YS5YsIS4urt5TWtXV1SxdupQxY8bg7Fz7oFdxcTHPPPMM33zzDT/++CNJSUnExsbStm1b7rrrroa2adc6dOiAu7s7VafLOf76LsyHi2jm6cyhVgeoNqpwcg2lmZMPPvkpeLarovegOGuXLCIiYjMafFotOTmZwYMHWz7Hx8cDMGbMGN566y1ycnKIj48nLy+PwMBARo8ezezZs2uNkZWVZbkO6az9+/ezefNm1q9fX+/3/uyzz8jKyuKRRx45b52TkxOpqaksW7aM06dPExgYyODBg1m1ahXe3k3vDfPuJU6cfGM31UVmnFq60eqhcD7+f68C0MylF64VhRxulUmg03X4tWxu5WpFRERsR4PDUWRkJBe6hnvKlClMmTLlgmMkJSWdt6xLly4XHBcgOjq63m08PDwu6W64pqDix0KuSfOmusqMs78nbR/pzp7kLykrKcbV0weTS0d8j37NzrBqBlyjZxuJiIicSy+edTDmY6WcensvzlXNcAnxxu/x63Fq4crOTz8BoJlLT0ymZvie2IkpqJqbbhlu5YpFRERsi8KRg3H29cCjly+nW1fQeuy1NPN04XBaKieyfsTJxRWDa3E2F1PslE4LzxB8mrtffFAREZEmROHIwZhMJrx/05GD1xRjcmmGYRhsXvE2AD7t+2Bq5o7viVR2hRl0DNBRIxERkV9SOHJAJicT/PRWlR++3kjOgX24uLlTab4RAN/jOykIruTaiLsvMIqIiEjTpHDkwMzl5Xy1vOao0XWDYjlT4opT5Rlcz+zD3ac513QItnKFIiIitkfhyIHtXLuGovzjeLf1xa15zRPK2+bvJrVDFe28Is57aa+IiIgoHDmsytISkj+pecXKraMeJvP7AqDmlNqBTgbXdL/PmuWJiIjYLIUjB5W/azuV5eUEdbkWn/a9KMovo1lVOa1P7cE1wKBX75utXaKIiIhNUjhyQHkZByjKrHnZ7uAxj5K15yQAbU7u5WCgGR/Xzni52+fLBUVERK40hSMHYxgGm95dAkDXgZEEdO7CicPFALQozCQlzETn4FhrligiImLTFI4czOm8HE4eOYzJyZn+9/4OgPwjhQA0L86mvL2ZHv1HWrNEERERm6Zw5GBaBwQx+pVFBNwShXebtpgrqig4XgZAuSkb7+atCQ30tXKVIiIitkvhyAF5eLfAK6jmGUanckowDHCpKCItpIgQ74FWrk5ERMS2KRw5uBNHaq43al6STX5ANV1veMDKFYmIiNg2hSMHl3+4CKi53si9lYnre95o5YpERERsm8KRgzt+4DgALmVHcW7TGTdnJytXJCIiYtsUjhyYYRjk59ZcjH3SM5vO7YZYuSIRERHbp3DkwEoLKqiobIbJqOJk6xx63Khb+EVERC5G4ciB5WeXAOBZeoyKttWEhYZYuSIRERHbp3DkwPIPnQbAqzgbI6AdJpPJugWJiIjYAYUjB3Zsbw4AzSqzad9psJWrERERsQ8KRw7seHbNa0MK3bO5qfc9Vq5GRETEPigcOSijCkorPAE40+IonTp0tnJFIiIi9kHhyEGZS5qBqRnO5hKq/Jx0vZGIiMglUjhyUJXHygHwKsmmxXUDrFyNiIiI/VA4clDlR0/X/EPVUfoPGG3VWkREROyJwpGDqiisOY1W4ZpNeKfrrFyNiIiI/VA4clDVRlsAzC1O6XojERGRBlA4ckAlp8uocvYGoxqPcF9rlyMiImJXFI4c0IGvdwDgceY4N94+xsrViIiI2BeFIwe0/+vtAJiqsul2XX8rVyMiImJfFI4c0Jn8mmuMql3ydL2RiIhIAzU4HG3atInY2FiCgoIwmUysXr261vq8vDzGjh1LUFAQnp6eDBs2jPT09AuOGRkZiclkOu/rjjvusGyTkJBw3vqAgIBa4xiGQUJCAkFBQXh4eBAZGUlaWlpDW7R7RrUfAM6+1VauRERExP40OByVlJTQs2dPFi5ceN46wzCIi4sjIyODNWvWsHPnTkJDQ4mKiqKkpKTeMT/66CNycnIsX7t378bJyYl777231nbdunWrtV1qamqt9S+//DKJiYksXLiQ7du3ExAQwNChQykqKmpom3brRF42Fa41obHjzX2tXI2IiIj9cW7oDjExMcTExNS5Lj09na1bt7J79266desGwKJFi/Dz82PFihWMHz++zv18fHxqfV65ciWenp7nhSNnZ+fzjhadZRgGCxYsYObMmYwcORKAt99+G39/f5YvX87jjz/eoD7tVWV5GR6uG6kucuem6KetXY6IiIjdaXA4upDy8ppXVri7u1uWOTk54erqyubNm+sNR7+0ZMkS7r//fry8vGotT09PJygoCDc3NyIiIpg3bx6dOnUCIDMzk9zcXKKjoy3bu7m5MWjQILZs2VJnOCovL7fUDFBYWPMWe7PZjNlsvsSubUubwBDue3kmGzZswAC77QN+rl09WJ8j9OEIPYD6sCWO0AM4Rh+NXXujhqOuXbsSGhrKjBkzWLx4MV5eXiQmJpKbm0tOTs4ljbFt2zZ2797NkiVLai2PiIhg2bJldOnShby8PF588UUGDBhAWloabdq0ITc3FwB/f/9a+/n7+3Po0KE6v9f8+fN5/vnnz1v+5Zdf4unpeUn12rINGzZYu4RG4Qh9OEIP4Bh9OEIPoD5siSP0APbdR2lpaaOO16jhyMXFhQ8//JBx48bh4+ODk5MTUVFR9Z6Gq8uSJUvo3r07ffvWvl7m3DF69OhB//79CQsL4+233yY+Pt6y7pd3ZxmGUe8dWzNmzKi1b2FhIcHBwQwePJg2bdpccs22xmw2s2HDBoYOHYqLi4u1y7lsjtCHI/QAjtGHI/QA6sOWOEIP4Bh95OfnN+p4jRqOAHr37k1KSgoFBQVUVFTg6+tLREQEffr0uei+paWlrFy5khdeeOGi23p5edGjRw/LnXBnr0XKzc0lMDDQst2xY8fOO5p0lpubG25ubuctd3FxsdtfkHOpD9vhCD2AY/ThCD2A+rAljtAD2HcfjV33FXvOUcuWLfH19SU9PZ3k5GRGjBhx0X3ef/99ysvLefDBBy+6bXl5OXv37rUEoY4dOxIQEFDrsGBFRQUbN25kwIABl9+IiIiINCkNPnJUXFzMgQMHLJ8zMzNJSUnBx8eHkJAQPvjgA3x9fQkJCSE1NZWpU6cSFxdX60Lp0aNH065dO+bPn19r7CVLlhAXF1fnKa1nnnmG2NhYQkJCOHbsGC+++CKFhYWMGVPzegyTycS0adOYN28e4eHhhIeHM2/ePDw9PRk1alRD2xQREZEmqsHhKDk5mcGDB1s+n71mZ8yYMbz11lvk5OQQHx9PXl4egYGBjB49mtmzZ9caIysri2bNah+02r9/P5s3b2b9+vV1ft8jR47wwAMPcOLECXx9fenXrx9bt24lNDTUss2zzz7LmTNnmDRpEqdOnSIiIoL169fj7e3d0DZFRESkiWpwOIqMjMQwjHrXT5kyhSlTplxwjKSkpPOWdenS5YLjrly58qK1mUwmEhISSEhIuOi2IiIiInXRu9VEREREzqFwJCIiInIOhSMRERGRcygciYiIiJxD4UhERETkHApHIiIiIudo9NeH2LOzjxIoKiqy20eoQ817ckpLSyksLFQfVuYIPYBj9OEIPYD6sCWO0AM4Rh9FRUUAF3wkUEMoHJ3j7IvrOnbsaOVKREREpKHy8/Np2bLlrx5H4egcPj4+QM0TvBvjh2sthYWFBAcHc/jwYVq0aGHtci6bI/ThCD2AY/ThCD2A+rAljtADOEYfBQUFhISEWP6O/1oKR+c4+0qTli1b2u0vyLlatGihPmyEI/QAjtGHI/QA6sOWOEIP4Bh9/PLVZJc9TqOMIiIiIuIgFI5EREREzqFwdA43NzfmzJmDm5ubtUv5VdSH7XCEHsAx+nCEHkB92BJH6AEco4/G7sFkNNZ9byIiIiIOQEeORERERM6hcCQiIiJyDoUjERERkXMoHImIiIicQ+HoHIsWLaJjx464u7vTu3dvvvrqK2uX1CAJCQmYTKZaXwEBAdYu64I2bdpEbGwsQUFBmEwmVq9eXWu9YRgkJCQQFBSEh4cHkZGRpKWlWafYC7hYH2PHjj1vbvr162edYusxf/58brrpJry9vfHz8yMuLo59+/bV2sYe5uNS+rD1+Xjttde4/vrrLQ/l69+/P//73/8s6+1hHuDifdj6PNRl/vz5mEwmpk2bZllmL/Nxrrr6sIf5uNjfucaaC4Wjn6xatYpp06Yxc+ZMdu7cyS233EJMTAxZWVnWLq1BunXrRk5OjuUrNTXV2iVdUElJCT179mThwoV1rn/55ZdJTExk4cKFbN++nYCAAIYOHWp5yaCtuFgfAMOGDas1N2vXrr2KFV7cxo0bmTx5Mlu3bmXDhg1UVlYSHR1NSUmJZRt7mI9L6QNsez7at2/PSy+9RHJyMsnJyQwZMoQRI0ZY/iNvD/MAF+8DbHsefmn79u288cYbXH/99bWW28t8nFVfH2Af83Ghv3ONNheGGIZhGH379jUmTJhQa1nXrl2N6dOnW6mihpszZ47Rs2dPa5dx2QDj448/tnyurq42AgICjJdeesmyrKyszGjZsqXx+uuvW6HCS/PLPgzDMMaMGWOMGDHCKvVcrmPHjhmAsXHjRsMw7Hc+ftmHYdjnfLRu3dr4xz/+YbfzcNbZPgzDvuahqKjICA8PNzZs2GAMGjTImDp1qmEY9vfvRX19GIZ9zMeF/s415lzoyBFQUVHBjh07iI6OrrU8OjqaLVu2WKmqy5Oenk5QUBAdO3bk/vvvJyMjw9olXbbMzExyc3NrzYubmxuDBg2yu3kBSEpKws/Pjy5duvDoo49y7Ngxa5d0QQUFBcDPL2S21/n4ZR9n2ct8VFVVsXLlSkpKSujfv7/dzsMv+zjLXuZh8uTJ3HHHHURFRdVabm/zUV8fZ9nDfNT3d64x50IvngVOnDhBVVUV/v7+tZb7+/uTm5trpaoaLiIigmXLltGlSxfy8vJ48cUXGTBgAGlpabRp08ba5TXY2Z99XfNy6NAha5R02WJiYrj33nsJDQ0lMzOT2bNnM2TIEHbs2GGTT6U1DIP4+HgGDhxI9+7dAfucj7r6APuYj9TUVPr3709ZWRnNmzfn448/5rrrrrP8R95e5qG+PsA+5gFg5cqVfPfdd2zfvv28dfb078WF+gD7mI8L/Z1rzLlQODqHyWSq9dkwjPOW2bKYmBjLP/fo0YP+/fsTFhbG22+/TXx8vBUr+3XsfV4A7rvvPss/d+/enT59+hAaGsp///tfRo4cacXK6vbEE0/w/fffs3nz5vPW2dN81NeHPczHNddcQ0pKCqdPn+bDDz9kzJgxbNy40bLeXuahvj6uu+46u5iHw4cPM3XqVNavX4+7u3u929n6fFxKH/YwHxf6O3f24vHGmAudVgPatm2Lk5PTeUeJjh07dl4CtSdeXl706NGD9PR0a5dyWc7egeBo8wIQGBhIaGioTc7Nk08+yb///W++/PJL2rdvb1lub/NRXx91scX5cHV1pXPnzvTp04f58+fTs2dPXn31Vbubh/r6qIstzsOOHTs4duwYvXv3xtnZGWdnZzZu3Mhf//pXnJ2dLT9zW5+Pi/VRVVV13j62OB+/dO7fucb8d0PhiJp/eXv37s2GDRtqLd+wYQMDBgywUlW/Xnl5OXv37iUwMNDapVyWjh07EhAQUGteKioq2Lhxo13PC0B+fj6HDx+2qbkxDIMnnniCjz76iC+++IKOHTvWWm8v83GxPupii/PxS4ZhUF5ebjfzUJ+zfdTFFufhtttuIzU1lZSUFMtXnz59+N3vfkdKSgqdOnWyi/m4WB9OTk7n7WOL8/FL5/6da9R/Nxp6pbijWrlypeHi4mIsWbLE2LNnjzFt2jTDy8vL+PHHH61d2iV7+umnjaSkJCMjI8PYunWr8Zvf/Mbw9va26R6KioqMnTt3Gjt37jQAIzEx0di5c6dx6NAhwzAM46WXXjJatmxpfPTRR0ZqaqrxwAMPGIGBgUZhYaGVK6/tQn0UFRUZTz/9tLFlyxYjMzPT+PLLL43+/fsb7dq1s6k+Jk6caLRs2dJISkoycnJyLF+lpaWWbexhPi7Whz3Mx4wZM4xNmzYZmZmZxvfff28899xzRrNmzYz169cbhmEf82AYF+7DHuahPr+8y8te5uOXzu3DXubjYn/nGmsuFI7O8fe//90IDQ01XF1djRtvvLHWrb/24L777jMCAwMNFxcXIygoyBg5cqSRlpZm7bIu6MsvvzSA877GjBljGEbNrZlz5swxAgICDDc3N+PWW281UlNTrVt0HS7UR2lpqREdHW34+voaLi4uRkhIiDFmzBgjKyvL2mXXUlf9gLF06VLLNvYwHxfrwx7m45FHHrH8t8jX19e47bbbLMHIMOxjHgzjwn3YwzzU55fhyF7m45fO7cNe5uNif+caay5MhmEYl3k0S0RERMTh6JojERERkXMoHImIiIicQ+FIRERE5BwKRyIiIiLnUDgSEREROYfCkYiIiMg5FI5EREREzqFwJCIiInIOhSMRERGRcygciYiIiJxD4UhERETkHApHIiIiIuf4///SmNEUvFOsAAAAAElFTkSuQmCC", "text/plain": [ "
" ] @@ -74,7 +74,7 @@ }, { "cell_type": "code", - "execution_count": 11, + "execution_count": 31, "id": "ef53616a-88a3-4722-84e2-92727e5399a1", "metadata": {}, "outputs": [ @@ -82,10 +82,10 @@ "name": "stdout", "output_type": "stream", "text": [ - "[-2.63160819e+01 -2.90289292e+01 -4.46549889e+01 -3.95001107e-13\n", - " 2.51126532e-13 -3.89769305e-13 -5.37131749e-13 -2.90365052e-13\n", - " -1.76573343e-13 1.18914517e-13]\n", - "-100.00000000000131\n" + "[-1.50511966e+01 -1.54865781e+01 -1.66133444e+01 -1.94936336e+01\n", + " -3.32803095e+01 -6.18665346e-03 -2.60712180e-02 2.81040181e-02\n", + " 7.22067838e-03 1.38265742e-02]\n", + "-99.90816871499192\n" ] } ], @@ -97,9 +97,115 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 32, "id": "8ced6e89-e8fd-4a7b-8b77-ca5e610581e9", "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "array([19.7515328 , 19.75153236, 19.75153185, 19.75153275, 19.75153281,\n", + " 19.76498845, 19.77043627, 19.77362338, 19.77547235, 19.7763392 ])" + ] + }, + "execution_count": 32, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "ml.head(rw, 0)" + ] + }, + { + "cell_type": "code", + "execution_count": 37, + "id": "801170e6-538e-4400-81fe-817734bc04e7", + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "array([2.06829668e-03, 6.74900159e-06, 3.32024764e-08, 2.57897091e-10,\n", + " 3.46416156e-12, 8.82738725e-14, 4.63575933e-15, 5.37054606e-16,\n", + " 1.44345016e-16])" + ] + }, + "execution_count": 37, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "from scipy.special import k0\n", + "k0(5 / ml.aq.lab[1:])" + ] + }, + { + "cell_type": "code", + "execution_count": 38, + "id": "8a558620-5871-44ee-9d3a-9a21f4d4b6e0", + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "array([[ 1.00000000e-01, 4.41707654e-01, 4.25325404e-01,\n", + " -3.98470231e-01, 3.61803399e-01, -3.16227766e-01,\n", + " 2.62865556e-01, -2.03030724e-01, -1.38196601e-01,\n", + " -6.99596196e-02],\n", + " [ 1.00000000e-01, 3.98470231e-01, 2.62865556e-01,\n", + " -6.99596196e-02, -1.38196601e-01, 3.16227766e-01,\n", + " -4.25325404e-01, 4.41707654e-01, 3.61803399e-01,\n", + " 2.03030724e-01],\n", + " [ 1.00000000e-01, 3.16227766e-01, -3.10776816e-17,\n", + " 3.16227766e-01, -4.47213595e-01, 3.16227766e-01,\n", + " -2.38891966e-16, -3.16227766e-01, -4.47213595e-01,\n", + " -3.16227766e-01],\n", + " [ 1.00000000e-01, 2.03030724e-01, -2.62865556e-01,\n", + " 4.41707654e-01, -1.38196601e-01, -3.16227766e-01,\n", + " 4.25325404e-01, -6.99596196e-02, 3.61803399e-01,\n", + " 3.98470231e-01],\n", + " [ 1.00000000e-01, 6.99596196e-02, -4.25325404e-01,\n", + " 2.03030724e-01, 3.61803399e-01, -3.16227766e-01,\n", + " -2.62865556e-01, 3.98470231e-01, -1.38196601e-01,\n", + " -4.41707654e-01],\n", + " [ 1.00000000e-01, -6.99596196e-02, -4.25325404e-01,\n", + " -2.03030724e-01, 3.61803399e-01, 3.16227766e-01,\n", + " -2.62865556e-01, -3.98470231e-01, -1.38196601e-01,\n", + " 4.41707654e-01],\n", + " [ 1.00000000e-01, -2.03030724e-01, -2.62865556e-01,\n", + " -4.41707654e-01, -1.38196601e-01, 3.16227766e-01,\n", + " 4.25325404e-01, 6.99596196e-02, 3.61803399e-01,\n", + " -3.98470231e-01],\n", + " [ 1.00000000e-01, -3.16227766e-01, -2.90992345e-16,\n", + " -3.16227766e-01, -4.47213595e-01, -3.16227766e-01,\n", + " -5.07992125e-17, 3.16227766e-01, -4.47213595e-01,\n", + " 3.16227766e-01],\n", + " [ 1.00000000e-01, -3.98470231e-01, 2.62865556e-01,\n", + " 6.99596196e-02, -1.38196601e-01, -3.16227766e-01,\n", + " -4.25325404e-01, -4.41707654e-01, 3.61803399e-01,\n", + " -2.03030724e-01],\n", + " [ 1.00000000e-01, -4.41707654e-01, 4.25325404e-01,\n", + " 3.98470231e-01, 3.61803399e-01, 3.16227766e-01,\n", + " 2.62865556e-01, 2.03030724e-01, -1.38196601e-01,\n", + " 6.99596196e-02]])" + ] + }, + "execution_count": 38, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "ml.aq.eigvec" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "575756e5-99b3-40cf-9b70-690268c257bf", + "metadata": {}, "outputs": [], "source": [] } diff --git a/timml/well.py b/timml/well.py index 5c23d9c2..f429c5a9 100644 --- a/timml/well.py +++ b/timml/well.py @@ -541,3 +541,52 @@ def initialize(self): def setparams(self, sol): self.parameters[:, 0] = sol + + def potinf(self, x, y, aq=None): + if aq is None: + aq = self.model.aq.find_aquifer_data(x, y) + rv = np.zeros((self.nparam, aq.naq)) + if aq == self.aq: + pot = np.zeros(aq.naq) + r = np.sqrt((x - self.xw) ** 2 + (y - self.yw) ** 2) + if r < self.rw: + r = self.rw # If at well, set to at radius + if aq.ilap: + pot[0] = np.log(r / self.rw) / (2 * np.pi) + pot[1:] = -k0(r / aq.lab[1:]) / (2 * np.pi) / k0(self.rw / aq.lab[1:]) + else: + pot[:] = -k0(r / aq.lab) / (2 * np.pi) / k0(self.rw / aq.lab) + rv[:] = self.aq.coef[self.layers] * pot + #rv[:] = pot + return rv + + def disvecinf(self, x, y, aq=None): + if aq is None: + aq = self.model.aq.find_aquifer_data(x, y) + rv = np.zeros((2, self.nparam, aq.naq)) + if aq == self.aq: + qx = np.zeros(aq.naq) + qy = np.zeros(aq.naq) + rsq = (x - self.xw) ** 2 + (y - self.yw) ** 2 + r = np.sqrt(rsq) + xminxw = x - self.xw + yminyw = y - self.yw + if r < self.rw: + r = self.rw + rsq = r**2 + xminxw = self.rw + yminyw = 0.0 + if aq.ilap: + qx[0] = -1 / (2 * np.pi) * xminxw / rsq + qy[0] = -1 / (2 * np.pi) * yminyw / rsq + kone = k1(r / aq.lab[1:]) / k0(self.rw / aq.lab[1:]) + qx[1:] = -kone * xminxw / (r * aq.lab[1:]) / (2 * np.pi) + qy[1:] = -kone * yminyw / (r * aq.lab[1:]) / (2 * np.pi) + else: + kone = k1(r / aq.lab) / k0(self.rw / aq.lab) + qx[:] = -kone * xminxw / (r * aq.lab) / (2 * np.pi) + qy[:] = -kone * yminyw / (r * aq.lab) / (2 * np.pi) + rv[0] = self.aq.coef[self.layers] * qx + rv[1] = self.aq.coef[self.layers] * qy + #rv[0], rv[1] = qx, qy + return rv