{ "cells": [ { "cell_type": "markdown", "id": "21de9f01", "metadata": {}, "source": [ "# Lumped capsule\n", "\n", "By default the packed-bed model resolves the temperature across each capsule. With\n", "`options={\"storage element\": \"lumped\"}` each capsule has a uniform temperature instead: one\n", "enthalpy per capsule, heated through its surface. This is much cheaper, and accurate when\n", "conduction inside a capsule is fast compared with heat transfer at its surface, that is when\n", "the capsule's Biot number, `Bi = h R / k`, is small.\n", "\n", "`\"lumped with Jeffreson\"` puts part of the neglected internal resistance back, through the\n", "heat transfer coefficient: `h_eff = h / (1 + Bi / 5)`, whichever coefficient the \"heat\n", "transfer coefficient\" option gives." ] }, { "cell_type": "code", "execution_count": 1, "id": "863f9410", "metadata": { "execution": { "iopub.execute_input": "2026-11-01T11:16:48.683482Z", "iopub.status.busy": "2026-11-01T11:16:48.683315Z", "iopub.status.idle": "2026-11-01T11:16:51.996816Z", "shell.execute_reply": "2026-11-01T11:16:51.995670Z" } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Biot number, solid: 6.9\n", "Biot number, liquid: 18.3\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "resolved: solve time 1.713 s\n", "lumped: solve time 15.434 ms\n", "lumped with Jeffreson: solve time 6.557 ms\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAhQAAAFzCAYAAABmY5CaAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjExLjIsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvgI3uAAAAAAlwSFlzAAAPYQAAD2EBqD+naQAAdKlJREFUeJzt3XdcVfX/B/DX3ZfLuAgoQ5YKqLhw77ScWY5yj4bZULM0m9r3Z7Y0K0tbVmqmmZpZ4kjN1JzhQMG9QfZeF7jc/fn9cbgHrgzhcuEy3s+H98E9n7Pe54KcN5/zGQLGGAMhhBBCSA0I7R0AIYQQQho+SigIIYQQUmOUUBBCCCGkxiihIIQQQkiNUUJBCCGEkBqjhIIQQgghNUYJBSGEEEJqTGzvAOqCyWRCcnIynJ2dIRAI7B0OIYQQ0mAwxpCfnw8fHx8IhRXXQzSJhCI5ORl+fn72DoMQQghpsBISEuDr61vh+iaRUDg7OwPgPgwXFxc7R0MIIYQ0HCqVCn5+fvy9tCJNIqEwP+ZwcXGhhIIQQgixwoOaDFCjTEIIIYTUGCUUhBBCCKkxSigIIYQQUmOUUBBCCCGkxiihIIQQQkiNUUJBCCGEkBqza7fRiIgIGI3GMuVeXl4ICgqyKEtISEBaWhpCQkKo6ychhBBSz9g1oXjvvfegVqv5ZZ1Oh3PnzmHRokVYtmwZAECj0WD69OnYv38/AgICEBcXhxUrVuCVV16xV9iEEEIIuY9dE4qDBw9aLG/duhXTpk3DM888w5e9//77OHv2LO7evQtvb2+Eh4fjiSeeQK9evdC7d++6DpkQQggh5ahXbSjWr1+Phx56CG3btuXLNmzYgOeffx7e3t4AgHHjxqFjx47YsGGDvcIkhBBCyH3qzdDb9+7dw5EjR7Bp0ya+LDk5GWlpaejevbvFtr169UJUVFSFx9JqtdBqtfyySqWyfcCEEFKfMAbo1UBRLqBVAfoiwKAFjFrAoCv+qi0pMxm5fZip1Mt43zIr3q54Gazic5e/ovJ4q7NPhdtbu08j1nowEDyszk9bbxKKn376Ca6urpgwYQJflp2dDQBwd3e32Nbd3Z1fV57ly5fj/fffr51ACSGkrjEG5KcC6VeBzNuAKglQpQCqZKAgFdDkcS+Twd6RkvpALGu6CYXJZMLPP/+MGTNmQC6X8+USiQQA1zCztKKiIkil0gqPt2jRIixcuJBfNs+URgghDYJeAySeBWKPA/GngbQrQFFO1fYVigGZCyBRAGIpIJYDIil3kxHLAJGMWxaKuJdAeN9LBAgElmVCEYDisopUOHFUJRNKPWCyqeptX81jNWb+/exy2nqRUBw8eBAJCQl44YUXLMr9/PwgFAqRlJRkUZ6UlAR/f/8KjyeTySCTyWolVkIIqRVGAxBzFLj0G3DjL0BfaLleIATc2gAt2gFKf8DFh3s5ewEOzQC5kntJFNW/URNiA/UioVi3bh169+6NTp06WZQrFAr069cPu3fvxowZMwAAhYWFOHToEJYuXWqHSAkhxMYMWiB6C3DyCyA3vqTcyQto9RAQOADw7gI0bwtIHOwXJyEPYPeEIjMzE3v27MF3331X7vqPPvoIw4YNw6JFi9C3b198/fXXaNGiBV588cU6jpQQQmyIMeDqTuDg/wGqRK5M4Q50HA90mgT49qCaBtKg2D2h+O+//9C3b19Mnjy53PWDBg3Cv//+i2+//RZnz55Fp06d8Msvv8DJyamOIyWEEBtRpQD73gBu7OWWnb2B/guA7s9QLQRpsASMNf5+NSqVCkqlEnl5eTRsNyHEfgw64OwPwLFPua6dQjEw8A1gwGuARP7g/Qmxg6reQ+1eQ0EIIU3CncPAvjeB7Lvcsk83YOw3gGcH+8ZFiI1QQkEIIbVJrwEOLQXOrOGWnTyBIUuALtMAYb0arJg0AowxGJkRYmHd394poSCEkNqSdg34YxaQfo1b7vUil0zInO0bF6lXGGPQGxm0BiO0BhP30lu+1xlN0OqLlw1GaPRGZBalIVUTjwxNPDJ18cjRJ0BlSMLj/s/ig4dn1/l1UEJBCCG2xhhw9keuB4dRCzg2B8atscvohaTqGGPQGU3Q6LmbuEZvgsZgRJGOu4FrDCbuq94IbfE6TfF2RcVl5ScFxV+L19+fHGgNpkpGCTdBIMmBUJYGkSwdQmk6hMVfBSJduXtcybhZa59RZSihIIQQW8pPBXbNA+78wy0HDwfGfgs4tbBvXA2U+SZfpDNCrTOiSG8s815b6sauuS8R4NaVJAIaPhHgbvhF95Xbr5uCEUJpFoSyNEjkmRDL0yGQpoFJMgCBvtw9BBBBIfCEs8gXrmJfuEn90FwWgKFB9mmXQwkFIYTYytVwYO9rQFE2N8T18A+5xxyNeDwJxhi0huIbfvENnrvhG/gbfpG+OAGweG/g32uKv5b/3gCTHW7yQgEgl4i4l1gIuUQEmUQEuUQIuVgEB2nJe75cIoJMLIRMXPxVUuq9WAiZRAShUI9sXRLSNQlILbqHpMJ7SCyIQ2JBPIzMci4W82VLhVK0UrZCa2VrtHZtjTaubdBa2Rr+zv6QiCR1/+FUgBIKQgipqaJcrgfH5e3cslcn4IkfAc9Qu4Z1P/PNv1BrQKHWiEKdAYVaAwq0Bqh1Ru6r1oBCnbF4m1LvLcoMfOJQpDfW2Q1fIhLAQcLdzBVSMeQSERyKb+TcexFk5mVxyU1eXrpMWpIglFkn4W76DhIRJCIBBDVIBNV6NWLzYnE37y7u5N5FTFYMYnJjkFiQCBMzlbuPg9gBrZUlCYP5fUunlhAJRVbHUlcooSCEkJq4+y+w62VuBlCBEBiwEBj0Njcxlw0YTQwFGgNUGj0KtAbkawzIL35fqOX+gr8/ISgoLi+TLOiMMNbi3V8qEsJByt2QFVIR/55LAMzvxWXWW74Xw0EqhINEXFIuNd/k61+vmDxtHpc45N5FTF4M7ubdRUxuDFIKUyrcx1nqjDbKNiWJg2trtFG2gaejJ4SVTcBWz1FCQQgh1tCpue6gZ3/glt1aA0/8APj14jfRG03I1xj4hMCcDFh81ZYuK/u+UGeslfAVxX/lO8nMX8VQyERwlInhKDV/FXNfZaLi91y5Qiou3r8kYXCQiCCuhzd8W2CMIVuTzSUMxYlDTC6XPGQWZVa4n5vcrUxtQxvXNnCXu9eo9qO+ooSCEEIqwBiDWmdEbpEeuWod8tR65BbpIUyOQu/od9CsKA4AcFw5BpscZiFjtx75mqN8MqDRl1+1bQ2ZWAhnuQQucjGc5FwCYHHz5xMA801fxCUJpZIF8z4OEhFEwsZ3Q6spo8mIlMIUxObFci9VLGJyYxCTF4NcbW6F+3kqPMvUNrRWtoar3LXOYq8PKKEghDR6JhNDvtZQnBDokFucGOSpS97nqvXIKyq7rDeWPCIQw4B54nDME4VDLDAhjbniLf1LOJbWBUgrrPD8CqkIznIxnOUSOMnEcJaL4SKXFJdx5ZZfxXCWWa6XihvnX//2oNarEavikoZ7eff45CEuLw46U/ldMQUQoKVTy5KEofhrK2UrOElpbimAEgpCSAOjN5qQU6hDZoEO2YU6ZBVqkV2oK04AuJqEkoSguGahSF+jhoMSkQBdHdLwoekbtDXdAQBEuzyMw63fRg9nDzziUDYhMCcMTjJxo30UUJ8xxpCmTiupbShOGu7l3UOaOq3C/SRCCQJcAtBK2QqBLoF8zUOgMhAOYpq4rTKUUBBC7IoxhrwiPTILtMgs0CGrQIfswuL3hVpkFZdlFWqRVZw4WMtBIoKrQgKlgwSuCglcHaTccqn3rg73LUv0cIj4EoL/vgJMBkCuBB77AmGdJiDMdh8DsZLGoEGcKq5MjcM91T0UGYoq3M9N7sYnDa2UrfiXj6NPg+hRUR9RQkEIqRWFWgMy8rXILNAiI1+LDPPX0mX5XOKgM1avrYFQALg5SuHuKIOboxRujsU3/+JEQFmcGLgqSpIEFwcJ5JJq3ihuHQT2vQ7kxnPLIY8Cj60ElC2rdxxSI3qTHkn5SYjPj0ecKg5xqjjEq+IRnx+P5IJkMJRf/SQWiOHn4meRNJjfK2XKOr6Kxo8SCkJItWj0RqTmaZCSp0FKXhFS8jRIzdOUJArFyYK6mr0TXORieDjL4F6cKLg7SeHuVLzsxJV5OJmTB2ntNirMSwIOvANc310cnC8w6lOg3WO1d84mzmAyIKUwBfGqeNxT3UO8Kh5x+VzikFyQDCOr+OfJRepiUctgThp8nX0hEdafgZ8aO0ooCCE8vdFkkSwk52qQmleEZHPykKtBVmH5jdbKo5CK0NxZhuZOMjR3lsGj+Ku5zMPZXC6FTFwPqpmNBm4Ojn8/BnQFgEAE9J0LDHoHkFHDu5oymAxILUxFQn6CRcIQp4pDYkEiDCZDhfs6iB3g7+wPfxd/BLgEwN+Z+xqoDEQzWbNG2Q2zoaGEgpAmgjGGjHwtEnKK+OQgufiruaYho0BbpbkM5BIhfJQO8HaVw1vpAG+l3CJxMCcPjrIG9Csm/jSw7w0g9TK37NsLePxLwKujfeNqQBhjUOlUSMxPREJBAhLzE5FUkITE/EQk5icipTCl0poGqVAKfxd/PlkonTy0ULSgpKGea0D/2wkhD2IyMaTla3AvU424rELcy7L8WpXHEFKxEN5KefGLSxa8XR3gU7zs4yqH0kHSeH65Z8dwA1Rd28Uty12BYe8DXZ8GhNQ7435qvRophSlIKUxBUn4SEgu4ZCGxIBFJ+UnI1+dXur9UKEVL55YIcC6VMLj4I8A5oMGPFNnUUUJBSANjNDGk5BUhLkuN2MxCi4QhLksNraHiBo5CAfikwLu4hsFH6QAvpZyvcXB3lDaeZKEyRTnA8c+BMz8AJj03bHbXp4BH/g9wam7v6OzCYDIgsyiTSxgKUvjEIa0wjX+v0qkeeJzmDs3h6+wLXydf+Dr7oqVTS365uaI5JQ2NFCUUhNRDBqMJSblFJTUMfI1DIRKyiyrtFSEWCuDnpkCAuwKB7o4WX32bKWiAJKMeiPwJOLqcSyoAoM0jwPCPAE/7TPtc20zMhBxNDjKLMpFRlIHMokzuvTqDX04tTEW6Or3SRxJmzlJneDl6cYlCcdLg5+wHXydf+Dj5QC6W18FVkfqGEgpC7IgxhqTcItxMzcfNtHzua2o+7mYUWIzQeD+pSAg/N4fiRMERgR5c0hDo7ggfVzkNpFQexoCb+4F//g/I4ganQvP2XCIRPNS+sVWT3qRHnjYPKq0Kudpc5GpzkafNs3jPJw/qTGRpsqqUKABcV0tPR094O3rD29EbXo5e8HYqfq/wgpejF40MScpFCQUhdURrMOJ2WgGuJatwNTkP11JUuJGSj3xt+S3bZWJhSQ2Dh2VNg7fSgeZiqCrGgJijwL/LgMSzXJljc+DhxVw7CVHd/RpkjKHIUIRCfSEK9AUo0BWgQF+AQn0h8nX55Zfr81GoKy4vtc4abnI3eDh4oLlDc3g4eHDvFdx7c/LgLnengZ2IVSihIKQWqDR6XEtWFScPXAJxJ70AhnLGfxYLBWjT3AltvZzR1ssZ7Yq/+igdIKSkoWbuneQSibhT3LJYDvSZCwx4DZC7VPkwjDFojVr+Jm++qd9/87+/vLwEoao1BQ8igADOUme4ylyhlCmhlCnhKnOFq8wVLjKXMomDu4M7jclAahUlFITUAGMMKXkaLnlIKU4gUvKQkF3+kL9KBwlCvV3QwccFHVq6oL23C1p7OFG7BluLP82NJRF7HDoABRI5CjpPREGXySiUylGQfp6/8VdaO1CqvLIxEqpLKBDCUeIIJ4kTnKROcJI4wVHiCGeJMxylxeXF68orV8qUcJG6UE0CqVcooSCkGrQGIy4l5uFsbDbOxmbjYmJuhXNL+CjlCPVRooOPC0J9uCSipatD0+hBUQOlHwsU6gtRaChEoa7kvVqv5m/4ar2av+Gr9WoUFqahMC8ehXo1CoUCFAb6wWD+vLOPAf8eq1FsAgjgKHEsPxmQOpdbbn5futxBTD8HpPGpFwlFamoqTp48CYVCgYcffhgODiUzup06dQrXr1+32N7NzQ1PPvlkXYdJmiCN3ogLcTk4HZOFM7HZiErIhe6+bpkioQDBLZwQ6s0lDu29XRDq7YJmjlI7RV33jCZjyc1eV8AlAfoH3PxLJwultlcb1DCx6s3tUYak7K82hVjB3eSlxX/xl7rZ33/jt9imVCKgkCioyyMhFbB7QrFixQp88MEHGDx4MBQKBd59913s2rUL/v7+AIBffvkFf//9N4YMGcLv07JlS0ooSK3QGoyIjs9FREwWIu5mlZtAuDtK0auVG3oGuqFHYDOEeDpXf9IpO9Ob9CgyFPE38CI9VyOgNqj5MrVeXe4N3yIhKH5pjBqbx2iuDVBIFPxNXyFRwFHM1RA4ihVwzEuEY+J5KFTJcDIxODJAETgIjt2ehpN7iMU+9HiAkNpl14Ti999/x7vvvovDhw9j0KBBAICYmBgYjZaNlrp3745169bZI0TSyOkMJlxKzEXE3SxExGThfFxOmYGhWjjL0LeNO/q0dkevVm5o7eFYZ9XVjDHoTXqLm7zaoC5z4y+dHNyfGBTpi8psrzNVfT6O6hALxfxf9RY3/1Kv8hIEJ2nZ7St8LKBTA1GbgYivS2YBlSiAbs9w8264+tfKtRFCKmfXhGLFihUYP348n0wAQOvWrctsl5GRgS1btkCpVKJHjx7w9PSsyzBJI8IYw43UfBy7lYFTdzIReS8HRXrLBNbDSYo+rd3Rt407+rZ2R6sqJhAmZoLGoOFu4oaikpt8BUlApTf94sSgSF8EA7NdY8D7iQViKCQK7iUufhW/d5A4FNcEOMJR6lhucnB/giAV1eJjnoJ0bkCqMz8ARdlcmcId6D0b6Pk8oHCrvXMTQh7IbgmFWq1GVFQUZs+ejWvXriEqKgo+Pj7o378/pFLLX0p3795FeHg4kpKSEBUVhc8//xxz586t8NharRZarZZfVqkePFQsabxUGj1O3s7E0ZvpOHorBekF+RAI9YBQB4FIB1cnE9r5yBDkKUOAhxguigJojMlQGdQIjy9C0d0iPkG4/6UxaEre10K1f2kykYy/4TuIHcpNAkp/NW/jKHYsd71CrIBEVM+7ETIGJJwFzq0FroZzQ2QDQLNAoN8rQJdpgFRhzwgJIcUEjFVlbkHbS0pKgq+vL0aNGoU7d+6ge/fuOH/+PEwmE/7++2++piIiIgI9evSARML94lu7di3mzJmDs2fPolu3buUee+nSpXj//ffLlOfl5cHFpep9z0ndM/+VX94Nu/SrdC2AxUvP3djVejUyCvORpS5AkaEIEOi4BEJQw8Z+VeQgduBu6PfdxM1V+RUlAffXEJReFgvt3uSp7ujUwJUdwNm1QOqlknLfXkCf2UD7sXU6IBUhTZlKpYJSqXzgPdRuCUVWVhY8PDwQFhaGM2fOQCqVwmAwoH///vDx8cHOnTsr3LdFixZYuHAh3nnnnXLXl1dD4efnRwlFLWGMIaMoA8kFyUgtTEW+Pp9PBPivxpLl+/+qL71tbf+VbyYWiPmbvoPEoeT9fS+5WG6xrBAryt9W4gC5SM7vQz0BrJQdA5xbz7WR0ORyZWI50GkC0PMFwCfMntER0iRVNaGwW4rv7u4ONzc3DBs2jH/EIRaLMXz4cGzevLnSfR0cHJCdnV3heplMBplMZtN4iSXGGC5mXMQft//AicQTyNJk2fwc/E1dJK/05s9MUtzL0ONaUhHS8xiYSQowKVxkjngoyAcjQv3R3tPDYn8aMbAe0RcB1/cCUZuA2OMl5a4BQM9Z3Ayg1D6CkHrPrnWGY8eOxaVLlyzKLl26hFatWgEAjEYjkpOT4efnx6+PiIhAQkIC+vbtW6exkhJ3c+9ixdkViEiJ4MuEAiE8FdyEQkqZEnKxHAqxAnKx3OIv99J/+ZdOFMzrzOUP+itfazDiyPV07DifiKO3MmAsHtJaKhJieKgnxndviYeCm9MkWfUVY0BKNFcTcfl3QJNXvELAzfzZ60UgeBhAXT0JaTDsmlB88MEH6N27NyZPnoz+/fvjzJkzOHToEA4fPgyA+yv40UcfRb9+/dChQwfEx8fjxx9/xKRJkzB27Fh7ht4kMcbw+63fseLsCuhMOkiEEoxqNQqj24xG1xZda7eFf/H5Lybm4Y/zidh9MRl5RSUjVIb5uWJ8d1+M7uwNV0XTGVCqwclPBa78AURvAdKulJQr/YGu04EuU4FmAfaLjxBiNbu1oTDLzMzEhg0bEBcXB39/f0ybNg2+vr78eq1Wi99++w1RUVFo1qwZBg4ciIcffrha56jq8x9SMRMzYWXkSmy6tgkAMKDlACzuvRh+zn4P2LPmsgq0+ONCIrZHJuJOesksi14ucjzRrSXGd/NFUAuaTrne0uQB13ZzNRH3TgDmUTBFMqD9aKDrDKDVIEBItUmE1Ef1vlFmXaKEomYYY/jk7CfYcmMLAGB+t/mY1XFWrQ/uFJtZiC//uYUDV1KhM3I3IZlYiJEdvTC+my/6B3nQFN71lV4D3P6bSyJuHQSMJY2k4dcb6DSRa2jp0Mx+MRJCqqTeN8okDceGqxv4ZOLD/h9iXNC4Wj1fZoEWP52MxbqTsfyw151aKjG1lz8e7+INFzk1qKyX9EXA3SPA9T3Ajb8AbanxX5q3BzpPBDqO58aQIIQ0OpRQkEpFJEdg1flVAIC3e75dq8mEycSw5Ww8lu+7jkIdN3rlwGAPvDWiHTr5KmvtvKQGNCrg9kHg+m7g9j+AXl2yzsWXq4XoNBHw7ADQ7JqENGqUUJAKZRVl4Z0T74CBYXzweMwInVFr57qXWYh3/ryE0zFcd+BOLZWY90gQhod60jTP9U1BOnDrAFcTEXMUMJaaF0Tpz7WLaD+ae7RB7SIIaTIooSAV+uTsJ8jWZCO4WTDe6VX+IGI1pdEb8f2xu/ju6F3oDCY4SER4e2RbPN03EEJqH1E/mIxAchRXE3H7IPe+NI8QoP0YLonw7kI1EYQ0UZRQkHIdSziGA/cOQCQQ4aP+H0Eultv+HLcysGTXFcRlcdXkA4M98PG4TvB3p7kZ7E6dzbWHuH0QuHMIUN83cJl3GND+cS6RaN7WLiESQuoXSihIGXqjHp9FfgYAeDr0aYS6h9r0+LGZhVi+7zoOXksDAHi6yLDk8Q4Y1cmLHm/Yi0EHJF8AYo4Bdw8DiedKuncCgMyFG3AqeDgQNBRwphl/CSGWKKEgZWy/tR1xqji4yd3wUpeXbHbcnEIdVh++jc2n42AwMYiEAjzbLxCvDQuBk4x+FOuUychNuhV7nHvFRQD6QsttWnTgRqsMHg749QLq+8ykhBC7ot/ixIJar8aPl34EAMzrOg+OEscaH7NAa8CmiHv4/uhdqDQGAMAj7Vpg0aPtEOzpXOPjkyowmYDMm0DsCSD2GHDvZMnkW2YKdyBwINB6EBA0DHCt/UHLCCGNByUUxMLvt35HtiYbvk6+eCLoiRodq0BrwMb/7mHtiRjkqrlhstt5OeN/j4ViQLCHLcIlFdEXAUkXgITTQPwZIOFM2QRC5gIE9AdaPcS9WoRSrwxCiNUooSA8jUGDDVc2AABe7PwixELrfjxS8zTYciYOm07H8YlEaw9HzHskCGPDWtLolrWhIKM4eTjNJQ/J0YBJb7mN2IF7dNF6EDfUtXcYIKJfAYQQ26DfJoS3L3YfsjRZ8Hb0xuNtHq/WvgajCSfuZGJHZCL+vpoKQ/Hsn609HPHqkGCM7uJDiYStaFTcTJ1JF7iGlElRQF582e2cvAD/3oBfH+6rV2dqB0EIqTWUUBAA3Hwdv1z7BQAwvf10SIRVu/HcTM3HHxcSsTMqCRn5JfM19Grlhmf7BWJEBy9KJGpCp+Zm5Uy6wI3/kHwByLwNoJwpeFqEcoNJ+ffhXq4BNCYEIaTOUEJBAABnUs/gTu4dOIgd8ERw5W0n0vM12B2djD8vJOFaSsl8DW6OUozp4oPJPf3Q3psmYasWxoD8FCD1MvdKuwKkXgGy71p23zRT+gE+XYGW3QCfboBPGCCn4ckJIfZT5YTi+eefr9aB161bV+1giP38fvN3AMCYNmPgIi0/GTh3Lxvf/XsHx25loPiJBiQiAR5p1wLju/licNsWkIqpUd8DaVRA5i0g4waQdg1Iu8wlD0XZ5W/v2JxLHny6lSQQTs3rNmZCCHmAKk9fLhAIMHny5Cod9LfffkN9mhWdpi+vXLYmG0N+HwKDyYAdo3egrZvlyIdR8Tn47O+b+O9uyWiJ3fxd8WQ3Xzze2RuuCmldh9wwFGZxXTUzbgAZN0te+cnlby8QAR7BgGdHwKsj4NmJ++rkSY8uCCF2UyvTl2/btq1K2/3222/VOSyxsz1398BgMqCDeweLZKJIZ8TKgzex/lQsGAPEQgEm9vDFCwNbo3VzJztGXE8wxg1JnR0LZMeUvHKKl+8frro0Jy+geQg3rbdXR8CrE/deYvshzgkhpC5UOaHYs2dPlQ9anW2J/YXfCQcAPBn8JF92JSkPr2yNQmwmN3riE11b4vXhIfBt1kTm2WAM0KqAwkxAlQSokoG8RO6rKgnISwJy47htKuPqD3i05ea7aN4WaN6Om0zLwbVOLoMQQupKlR95NGT0yKNit3Nu48ndT0IsFOPopKNQypQ4cCUVC36LgkZvgp+zCF/3K0SY8C6Ql8DN+WAylLzKU+ZHilWwrqJya/a5b//q7sMYYNAARTnFr1yAGVElLr6AW6viV2vu1awV4N4GkNZ8pFFCCLGnWnnkodFo8Ndff2H8+PHlrv/jjz/w2GOPQS6natuG4sC9AwCAAT4DoJQpER6VhIXboyFkBnzs9R+mardDeKyCxoJNgcQRcPEBlC0BF/PLh/vq6g80C6THFIQQgmomFOvWrcOdO3cqTCiOHz+O1NRUvPzyyzYJjtQuxhgOxHIJxchWI/HPtTQs3B4NF5aPP9zWoE1uNLehkxfQejD3l7fEARCKi19CQFC6V8d9DQctGhLaYF2Zhok2XieWAQ5ugEOzkhclC4QQUiXVSig2btyIH374ocL1zzzzDGbPnk0JRQNxLesa4vPjIRfJ0ULYDU9tvQAFU2O/cgW81TGA1BkY8TEQNp2GaCaEEFKpat0lbt26hZCQkArXh4SE4NatWzUOitSN/bH7AQD9vAfitd+uQ6s3YJvrWnhrYriuik/vAlq0t3OUhBBCGoJqjUJkNBphMFTQEA+ATqeD0VjFhmzErkzMxLefSEluh8ScIsx1PokwzRlALAembqVkghBCSJVVK6Ho0KED9u3bV+H6AwcOoEOHDjUOitS+6PRopKnTIBMqcPZaC3gJ87BQsIVbOXQp0LK7XeMjhBDSsFQroXjuuefw2muv4ciRI2XWHTlyBK+99lq1h+gm9nEo/hAAQKdqDzAJvvc/DJFOxQ3x3OtFO0dHCCGkoalWG4oXX3wRR48exZAhQxAWFoa2bduCMYZbt24hOjoa06dPx6xZs6odxH///YcjR45AoVBg8uTJaNmypcX6jIwMbNu2DWlpaejUqRMmTJgAkUhU7fMQDmMMRxOOAgDUue3Rt7kOXTKKByMb9iEgpM+WEEJI9VSrhkIgEGDLli3Ytm0b/P39cfHiRVy+fBn+/v747bffsHnzZgiqMeeAyWTCzJkzMWbMGOTk5CAnJwejRo3CjRs3+G3u3r2LTp06ITw8HEajEYsWLcKoUaOorUYNxOTFICE/AcwkgqEgGF/4HofAqAX8+wKBA+wdHiGEkAbIriNlfvPNN3j77bcRHR2N4OBgANyIXEVFRfD09AQATJgwAampqTh+/DiEQiHu3buHkJAQbNiwAdOnT6/SeWikTEvrLq/D6gurYShoixHK17EqYTKgzQOm7wCCh9k7PEIIIfVIVe+hdp1r+ptvvsG0adP4ZAIAXFxc+GTCYDDgr7/+wvTp0yEUcqEGBgZi0KBBCA8Pt0fIjcKe21z7CVNBKP7X+g6XTLj6A22G2DkyQgghDVW1EopDhw4hKCgIQUFBOHToUI1OnJ+fj5s3b2LgwIEIDw/H//3f/2HNmjVIS0vjt4mLi4NGo0GbNm0s9m3Tpk2l411otVqoVCqLF+FkqjMRo7oGABgTMgQet4pnhg2bwY18SQghhFihWneQOXPm4Pfff8fvv/+O2bNn1+jEeXl5AICVK1dizZo1kEgk2Lt3L4KDg3H27FkAgFqtBoAyVSxKpRKFhYUVHnv58uVQKpX8y8/Pr0axNiYbo/cDAgaTpiXe6uUL3DsBQACETbN3aIQQQhqwavXyyMzMREBAABhjyMzMrNGJnZ2dAQBubm74+++/+fJHH30UixYtwuHDh+Hk5AQAyM3Ntdg3NzeX3788ixYtwsKFC/lllUpFSUWxnbcOAgKgvbIPmicW1zL59wFc6fMhhBBivWrVUMybNw9jxozBmDFjMHfu3BqdWKlUwsfHBz179rQo79mzJ2JiYgAA/v7+cHR0xM2bNy22uXHjBtq3r3gUR5lMBhcXF4sXAa4kZyKXXQUAvNxrLHBjL7ei3eN2jIoQQkhjUK0aig8//BBnz56FyWRCnz59anzyKVOm4NixYzCZTBAKhdz4CEePonPnzgAAkUiEJ598Ehs3bsTs2bMhk8lw5coVnDx5Ejt27Kjx+Zuab/77BwKhHlK4YrCnP3DvFLei3Sj7BkYIIaTBs2u30by8PAwZMgRGoxF9+/bF2bNnkZGRgSNHjvANMVNSUjBw4EA4OTmhW7du2Lt3L4YPH47NmzdX+TzUbRQo0BrQ5/vXIHA9ir4tRuJHnx7AzheBFh2Auf/ZOzxCCCH1lM27jX7++edVPnlVt1UqlYiIiMAHH3yA0NBQLFmyBDdv3rTo1eHt7Y1Lly5h8eLF6Nq1K7Zv316tZIJw9l5MhknOPToa23YwEPMvtyJkuP2CIoQQ0mhUuYZCIBCgqpUZ1dm2LlANBTD+xwO4JXsTAHBs0lG4rRkIqJKAp8KBNg/bNzhCCCH1VlXvodVqQxEWFlbTuIgdpOQV4VJWJOQ+QJCyHdzUuVwyIZICfr3tHR4hhJBGoMoJxWeffVblg86YMcOqYEjt2HMxGSIFNxDYIL/+QMxRboVvL0CqsF9ghBBCGo0qJxRvvPFGbcZBatFfl1MgcroNAOjn0w848S23ovUgO0ZFCCGkMaGxlhu5zAItrmTcgFBcALlIjrAWYUACNxIp/PvaNTZCCCGNByUUjdzRmxkQKrjaiZ5ePSEtyADykwGBEGjZzc7REUIIaSwooWjk/r2RDrHjHQBAX5++QFIkt8KzAyB1tGNkhBBCGhNKKBoxg9GE47dSIXK4BwDo5dULSDzHrfTtWfGOhBBCSDXVKKHIysrCuXPnbBULsbGrySoUCuIhEOngLHVGcLNgILG4hqJlD/sGRwghpFGxKqHIzc3FuHHj4OHhgV69evHl48aNQ0REhM2CIzVzJjYLIkUsAKB7i+4QmoxAcjS3kmooCCGE2JBVCcVbb70FrVZbZhbQuXPn4sMPP7RJYKTmzsZmQ1ycUPTw6gFk3gYMRYDUGXAPsnN0hBBCGpNqjZRptnfvXpw+fRr+/v4W5T179sSxY8dsEhipGaOJ4UxsJkT+xTUUnt2BpMvcSq+OgJCazxBCCLEdq+4qOTk5aNasGQBu3g6zgoICCOlGVS/cTM1HIUuEQKSBQqxAO7d2QJo5oehk3+AIIYQ0Olbd/bt3746dO3cCsEwoVqxYgX79+tkmMlIj5+Nz+PYTXVt0hVgoBlKLEwrPjnaMjBBCSGNk1SOP5cuXY9SoUTh+/DgYY/i///s//P3337h8+TI98qgnLifm8glFD68eAGNA6hVuJdVQEEIIsTGraigGDhyIEydOoKioCMHBwdi2bRtatWqFiIgIi14fxH4ulkoount2B/JTAXUmIBABLdrbOTpCCCGNjVU1FN9//z1mz56NX3/91dbxEBso0hlxNzcGDs0KIRXK0NG9I3D3X26lRzAgcbBvgIQQQhodq2ooXn31VRiNRlvHQmzkWooKAvk9AECX5p0hEUmA9Gvcyhah9guMEEJIo2VVQhEaGorz58/bOhZiI5cTcyF0iAcAbnZRAMi8xX1t3s4+QRFCCGnUrHrk8dxzz2HKlClYsmQJQkNDIZVKLdaHhYXZIjZipctJKoiKE4rOzTtzhRnFg5A1D7FTVIQQQhozqxKK+fPnAwBmzpxZ7nrGmPURkRq7npYKkUs6AKCTRyeuh4e5hsKjrR0jI4QQ0lhZlVBkZGTYOg5iIyYTQ2z+DYhdAE8HH7g7uAOqZECr4np4uLexd4iEEEIaIasSCg8PD1vHQWwkMacIRmkcxAC6enbhCs2PO9xaAWKZ3WIjhBDSeFmVUBw6dKjS9UOHDrUqGFJzt9PzIZInAADCWhQnFPS4gxBCSC2zKqEYNmxYpeupDYX93EzN53t4dPagBpmEEELqhlXdRhljFi+j0YgbN25g0KBBNNiVnV1KuwuhWA0RJNyEYACQdZv76kEJBSGEkNphVQ3F/YRCIdq2bYt169bhiSeewLRp06q033fffYd9+/ZZlAUEBODbb7+t1jakxM3cK4AD4OsYzA1oBQDZ3BDccKMGmYQQQmqHTRIKs+bNmyM2NrbK21+6dAm5ubl45513+DIXF5dqb0M4jDGka29B6AB0bl48AZheA+Qlcu+phwchhJBaYlVCkZiYWKYsJycHn332GUJDqze0s5eXFx5//PEab0OAjHwtmCwOANDPtxtXmBsHgAEyF0Dhbr/gCCGENGpWJRR+fn7lloeEhGDr1q3VOtbFixcxadIkKJVKDBw4EDNmzIBQKKz2NgS4k5kDoTwFQKkuo9kx3Fe3VoBAYKfICCGENHZWJRTXr18vU9asWTN4enpW6zgSiQSPPPIIBg0ahKSkJLz11lvYtm0b/vrrLwiKb35V2eZ+Wq0WWq2WX1apVNWKq6GKTL4OgcAEEXOCj6MPV8gnFK3tFxghhJBGz6qEYsqUKYiOji53XVhYWIXr7vfxxx9btIcYOXIkwsLCsHPnTjz55JNV3uZ+y5cvx/vvv1+1i2lErmZyiV4zcWBJskUJBSGEkDpg1XODixcvlltuNBpx9erVKh/n/saVnTp1QkBAgMVMplXZ5n6LFi1CXl4e/0pISKhyTA1ZXAE3gJWfY3BJISUUhBBC6kC1aihOnz5d7nsAMJlMiIiIgL+/v9XBMMaQm5sLiURSo21kMhlksqY3xHSWLhYQAe3c2pcUUkJBCCGkDlQroejbt2+57828vLywevXqKh1Lr9dj06ZNeO655/jq+WXLliEvLw9jx46t8jaEYzAZoBFwvW96enfkCo16ILe4dqZZKztFRgghpCmoVkKRn58PgJscLDMz02KdRCKpVq2ASCRCZGQk3nvvPQQHByMhIQH5+fn45Zdf0LVr1ypvQzjXM+8CQj2YUYoevsUjYuYlAswIiOWAs5d9AySEENKoVSuhcHJyAgBoNJoan1goFGLNmjX46KOPcOXKFTRr1gwhISGQy+XV2oZwTideBgAI9D5opihO7HK5OT3g6k9dRgkhhNQqm46UaQ13d3cMGjSoxts0ddcz7wAAnIS+JYWlEwpCCCGkFlnVy0Oj0eC9995DaGgonJycIJfLLV6k7sXlc40vPaSlBh2jhIIQQkgdsSqhWLJkCf78808sXrwYhYWF2Lx5M1599VUIhUK8+eabto6RVEFaEZc8+CgCSwopoSCEEFJHrHrk8dtvv2HPnj3o3LkznnrqKYwfPx4TJkxA9+7d8c0339g6RvIAepMeKmMqAKC1a6nuoeaEQln+UOmEEEKIrVhVQ5GQkIAOHToAABwdHZGXlwcAeOyxxxAZGWm76EiVJKgSwGAEM0rRulnLkhV8DUWAfQIjhBDSZFiVUDDGIBKJAABBQUH4559/AADnz5+nqcXt4G7eXQCASdcC3koHrtCgA/KTuff0yIMQQkgtsyqhUCqV/PvXXnsNTz/9NHr06IFHH30UL774os2CI1UTk8s1yDRpW8BLWdxlVJUEMBM3BoVTCztGRwghpCmwqg1Fbm4u//6ZZ55BmzZtcPr0abRr1w6PP/64rWIjVXQnNxYAV0Ph6VLcy6Z0+wkag4IQQkgtsyqh+P777zF79mx+ecCAARgwYIDNgiLVE5+XBACQmNzhJCv+luYVD7ntSg0yCSGE1D6rHnm8+uqrMBqNto6FWClNnQIAcJN5lkxbXpxkwKVlBXsRQgghtmNVQhEaGlrp9OGk7hhNRuRoMwAAnopS83WouInCoPQtZy9CCCHEtqx65PHcc89hypQpWLJkCUJDQyGVSi3Wh4WF2SI2UgUZRRkwwQTGhPB18SxZQTUUhBBC6pBVCcX8+fMBADNnzix3PWPM+ohItaQWcgNaMb0SXp6KkhWq4oRCSQkFIYSQ2mdVQpGRkWHrOIiVkgu4sSZMelc0dy41fTxfQ0GPPAghhNQ+qxIKDw8PW8dBrJRSyDXIZHpXuDsWP3rS5AG6fO491VAQQgipA1Y1ygSAzMxMrF+/Hu+++y5fdvr0aer9UcfMCYXJ4Ipm5oTCXDshdwWkjvYJjBBCSJNiVUIRHR2N9u3bY+XKlVi2bBlf/vPPP2Pjxo02C448WEkbCle4KYoTCr79BD3uIIQQUjesSihef/11LFy4ENeuXbMonz17Nr788kubBEaqhq+h0LuimaOEK8wr7jLq4mOnqAghhDQ1VrWhiIyMRHh4OACUDKQEIDg4GLdu3bJJYKRqUgpK2lC4Od5XQ0FdRgkhhNQRq2ooRCIR8vPzy5TfuHEDbm5uNQ6KVI3epEe+nvs+SOACBwk3AyxUxbOMUoNMQgghdcSqhGLUqFFYunQpTCYTX0MRFxeHOXPmYPTo0TYNkFQsT5sHAGBMgGYOypLaInNC4UyPPAghhNQNqxKKlStXIiIiAt7e3jCZTOjUqROCgoKg0WiwfPlyW8dIKmBOKGCSw00hL1lRkMZ9dfYquxMhhBBSC6xqQ+Hp6Ynz589j165diIyMhMlkwuLFizF+/Pgyw3CT2pOrzQUAMIMj3J1Kfe75XLsKOHvXfVCEEEKaJKsSCgCQSqWYOHEiJk6caMt4SDXwCYVRgWbOxQmFQQsU5XDvqYaCEEJIHbE6oYiJicG3336L69evA+BmIJ03bx4CAwNtFRt5AL4NhVFR0sMjnxuXAiIp4NDMTpERQghpaqxqQ7Fnzx60a9cOx44dQ0BAAAICAnD06FG0bdsW+/bts3WMpAI5Gq4mghkVaGYe1MrcfsLJCyjVpZcQQgipTVbVULz++ut4//33sWjRIovy5cuX47XXXsOoUaOqdJwlS5Zgy5YtFmUhISFlkpJdu3bhq6++QlpaGjp16oQPP/wQQUFB1oTeqFjWUBQPasW3n6DHHYQQQuqOVTUUSUlJmDt3bpnyOXPmICkpqcrHSU9PR3BwMA4cOMC/fvjhB4tt9u7diwkTJmD06NFYv349GGMYOHAgsrOzrQm9UbFoQ8E/8jD38PC0T1CEEEKaJKsSitDQUFy8eLFM+cWLFxEaGlqtYzk6OiIoKIh/+fn5WaxfunQpZsyYgQULFqB3797YtGkTdDod1qxZY03ojUpJQuFYMo8H9fAghBBiB1YlFNOnT8fEiRPxxRdf4MSJEzh+/Di++OILTJo0CTNmzEB0dDT/epCTJ0+ic+fOGDhwIBYvXmwxAmd+fj4uXLiAESNG8GVSqRRDhw7F0aNHrQm9USn9yIOvoeDbUFANBSGEkLpjVRuK1157DQDXluJ+8+fPt1hmjFV4HDc3N7z99tsYNGgQkpKS8L///Q979+5FZGQkpFIpEhMTwRiDl5dlewAvLy9cunSpwuNqtVpotVp+WaVSVem6GprSjzyUDve3oaAaCkIIIXXHqoQiIyPDJif/8MMPIRJx809069YNXbt2RatWrbBt2zY8/fTTMJlMAACJRGKxn1QqhdForPC4y5cvx/vvv2+TGOuz0gmFs7z4W0ltKAghhNiBVQmFh4eHTU5uTibMfH19ERgYiKtXr1qcJysry2K7zMzMSmNYtGgRFi5cyC+rVKoybTMaOsYYcks98nCUmhMKqqEghBBS96we2AoACgsLkZOTU6bc19fXquNptVqkpqbC1dUVADfEt5+fHyIiIjBmzBh+u1OnTmHkyJEVHkcmk0Emk1kVQ0ORr8+HiXG1NI5iFwiFguJRMot7vzhRt1FCCCF1x6pGmdHR0ejatSucnJzg5+dX5lUVOp0Ob775Jl/7oFarMXfuXBiNRkyePJnfbvbs2Vi3bh1u3LgBAPjxxx8RExODF154wZrQG408TXHthEkKZ5kDV2hukCmUAAqaRp4QQkjdsaqG4rnnnkPHjh2xZs0avjahuiQSCXx9fdG5c2eYTCbk5OSgS5cuOHz4MFq3bs1v9/bbbyMxMRFdunSBo6MjRCIRNm/ejE6dOll13saidPsJJ9n97SdolExCCCF1y6qE4tq1azh27BicnZ2tPrFAIMD8+fMxf/58pKWlQalUQi6Xl9lOJBLhu+++w2effYbs7Gx4e3tDLK7Rk5pGQaXjeq4wo0OpBpk0SiYhhBD7sOqRR3BwMOLj420WhKenZ7nJRGmOjo7w8/OjZKKY2qAGADCTDE7y4l4wNAYFIYQQO7Hq7rxixQo8++yzeO+999CmTRsI7qteb9eunU2CIxUrMhRxb0xSOMuohwchhBD7siqhMBqNuH79OkaPHl3u+soGsyK2odaX1FDQGBSEEELszaqEYsGCBZg+fTpeeeUVqxtlkpoxP/KASVKqUSbVUBBCCLEPqxKK1NRUrFy5Ek5OTraOh1RR6RoKJ3MNRUGpXh6EEEJIHbKqUWa7du1w584dW8dCqqGkUaYUzvL75vGgQa0IIYTUMatqKMaPH4+pU6figw8+QFBQUJlGmWFhYbaIjVTCXEPBN8o06AB18RDl9MiDEEJIHbMqoXj33XcBAJMmTSp3PTXKrH3mXh58o8zCdG4FjZJJCCHEDuw62yixHt8ok0m5NhT5SdwyjZJJCCHEDuw62yixXpHeXEMh5Xp55KdyK5xa2DEqQgghTZVVjTIBbgrx9evX848/AOD06dMwGo02CYxUrkyjzMLiWiNHSigIIYTUPatnG23fvj1WrlyJZcuW8eU///wzNm7caLPgSMVKGmWa21BkcstOze0XFCGEkCbLqoTi9ddfx8KFC3Ht2jWL8tmzZ+PLL7+0SWCkcoX6QgClHnnwNRSUUBBCCKl7VrWhiIyMRHh4OABYdBkNDg7GrVu3bBIYqZy6uJeHgMmgkIoooSCEEGJXVtVQiEQi5Ofnlym/ceMG3Nyoy2JdMHcbVYgduKSOEgpCCCF2ZFVCMWrUKCxduhQmk4mvoYiLi8OcOXMqnDCM2I7epIfepAMAOEsduUJzGwpH6oFDCCGk7lmVUKxcuRIRERHw9vaGyWRCp06dEBQUBI1Gg+XLl9s6RnIffupyAE58QkE1FIQQQuzHqjYUnp6eOH/+PHbt2oXIyEiYTCYsXrwY48ePh1QqtXWM5D78xGBMBEeZDDAZS4bdpoSCEEKIHViVUISFhSE6OhoTJ07ExIkTy11Hak/J1OVSOEhEgDobAAMgAByoDQshhJC6Z9Ujj4sXL5ZbbjQacfXq1RoFRB6sZB4PqWUPD4UbILIqRySEEEJqpFp3n9OnT5f7HgBMJhMiIiLg7+9vm8hIhfhHHiYp5BLqMkoIIcT+qpVQ9O3bt9z3Zl5eXli9enXNoyKV4htlmmgMCkIIIfVDtRIK89gTHh4eyMzMtFgnkUggk8lsFxmpUOkaCgeJiLqMEkIIsbtqJRROTk4AAI1GUyvBkKop3ShTTjUUhBBC6gGrZxsl9lO6hkIhoXk8CCGE2F+96hJgMBggFAohFJbkOSaTCSaTyWI7gUAAkUhU1+HVG3wvDyaFg1QIZNAjD0IIIfZVb2ooXn/9dUgkEixcuNCifO7cuZBKpZDL5fyre/fudoqyfih55CGDg5RqKAghhNifVQlFWFiYVesqsn//fhw4cAAhISHlrn/yySdhMBj4V1MfOKtso0xKKAghhNiX3Qe2SklJwQsvvIDNmzfDwcHBmnCaHI2xuFEsP7CV+ZEHJRSEEELsw64DW5lMJsyYMQPz589H165dK9xu//79kMlkcHFxwcCBA/HZZ5+hTZs21Qm9UdEatQAAxsRwFOgBXfFU8tSGghBCiJ3YdWCrZcuWwWQy4fXXX69wm5CQEOzYsQODBg1CUlIS5s+fj0GDBuHKlStwdXUtdx+tVgutVssvq1SqKsfUEGgNxdfGxHA25XLvRVJA5mK3mAghhDRtdhvY6uzZs1i1ahU/W6m5JwdjDAaDAWIxF1rpRprBwcHYunUrvLy88Ntvv+Gll14q99jLly/H+++/X51La1D4GgqTBE6GHK7QsTkgENgxKkIIIU2Z3Qa2ioyMRG5uLoKCgvgyo9GIS5cu4dtvv4VWqy23a6hSqYSfnx/u3LlT4bEXLVpkkYioVCr4+fnVOOb6QmfUcW+YGI58QkGPOwghhNiP1d1GMzMzsX79erz77rt82enTp2E0Gqu0/9y5cy16bhgMBnTp0gWvvvoqDAZDheNM5OTkID4+Hi1btqzw2Ob2FqVfjUlJo0wJHHTZ3HtqkEkIIcSOrEoooqOj0b59e6xcuRLLli3jy3/++Wds3LjRZsFptVqMHj0a//33H/Ly8nD58mVMmjQJSqUS06dPt9l5GprSjTJl2iyukBIKQgghdmRVQvH6669j4cKFuHbtmkX57Nmz8eWXX1odjFgstqiZkMlkmDdvHv73v/8hMDAQY8eOhY+PD86dO4fmzZvuDVRTqlGmVGuuoaBHHoQQQuzHqqG3IyMjER4eDoAbBtssODgYt27dsjqYyMjIMmUjRozAiBEjrD5mY2ROKJhJAnERjUFBCCHE/qyqoRCJRHyPj9Ju3LgBNze3GgdFKmd+5CEVSSFQU0JBCCHE/qxKKEaNGoWlS5fCZDLxNRRxcXGYM2cORo8ebdMASVnmXh5ykYyG3SaEEFIvWJVQrFy5EhEREfD29obJZEKnTp0QFBQEjUaD5cuX2zpGch+diauhcBDLSw27TW0oCCGE2I9VbSg8PT1x/vx57Nq1ix+YavHixRg/fjykUqmtYySlMMagN5lrKKQlCYWCEgpCCCH2Y1VCAQBSqRQTJ07ExIkTbRkPeQBdcTIBAM0kAsCk5xYU1HaFEEKI/VQ5oTh58mSVDzpgwACrgiEPpjGUjFLaQlScXIikgERhp4gIIYSQaiQUAwcOrPJBGWNWBUMezNwgkzEBPETFyYWDG83jQQghxK6qnFAUFRXVZhykisxdRsHEcBeoufcOzewXECGEEIJq9PKQy+X8q0+fPhbL968jtaf0TKNuwuKxQKj9BCGEEDuzqtvoxYsXyy03Go24evVqjQIilStdQ6FEIfeeaigIIYTYWbV6eZw+fbrc9wBgMpkQEREBf39/20RGylUydbkESlZcQ0EJBSGEEDurVkLRt2/fct+beXl5YfXq1TWPilTIPHU5M4nhDEooCCGE1A/VSijM83d4eHggMzPTYp1EIoFMJrNdZKRcpWsoHJmKe09tKAghhNhZtRIKJycnAIBGo3nAlqS28I0ymRgKY3FCQTUUhBBC7MyqkTIPHTpU6fqhQ4daFQx5MH5gK5MECkMW996BaigIIYTYl1UJxbBhwypdTwNb1Z6SRx5iyAxUQ0EIIaR+sKrbKGPM4mU0GnHjxg0MGjQIv/76q61jJKWUPPKQQKbP5QqpDQUhhBA7syqhKHMQoRBt27bFunXraPryWsaPQ2ESQaLL495TDQUhhBA7s0lCYda8eXPExsba8pDkPuaEQsIEEDIjV0gJBSGEEDuzqg1FYmJimbKcnBx89tlnCA0NrXFQpGLmhEJqbqcidgAkDnaMiBBCCLEyofDz8yu3PCQkBFu3bq1RQKRy5oRCzkxcAdVOEEIIqQesSiiuX79epqxZs2bw9PSscUCkcuZeHg6UUBBCCKlHrEoo2rVrZ+s4SBWZx6FwMLefkCvtGA0hhBDCqXajTMYYtm7discffxxt27ZFSEgIRo8eja1bt9L4E3XA/MiDr6GghIIQQkg9UK2Ewmg0YuLEiZg2bRoyMzPRp08f9O3bFxkZGZg2bRomTZoEo9FYW7ESAEUGLqFwZHquQO5ix2gIIYQQTrUeeaxevRpnzpxBREQE+vTpY7EuIiICkyZNwldffYXXXnut2oGo1WrExMTAw8MDXl5eZdbn5uYiMzMT/v7+kEql1T5+Y1Gk5x55ODIDVyCjhIIQQoj9VauGYsOGDfjhhx/KJBMAN535999/jw0bNlgVyMyZM9GpUyd88sknFuUGgwGzZs2Cp6cnHnroIbRo0QKbNm2y6hyNgaa4hsIFxQkFPfIghBBSD1Qrobh9+zYeeeSRCtc/8sgjuHXrVrWD+PHHH5GamopOnTqVWbds2TLs3bsX165dQ3JyMr766ivMnDkT0dHR1T5PY6ApbkOhBD3yIIQQUn9U65GHTCZDfn4+5HJ5uesrW1eRq1evYunSpTh9+jTGjBlTZv2PP/6IF154AW3atAEAPP3001i+fDnWrVuHb775plrnagy0BnNCUTxJGD3yIKRSjDEYDAZq30VIBUQiEcRiMQQCQY2OU62EolevXli7di0WL15c7vq1a9eiV69eVT5eUVERJk+ejM8//xz+/v5l1qempiIpKQm9e/e2KO/bty/Onz9fndAbDXMvD1cUz+lBNRSEVEin0yElJQVqtdreoRBSrykUCnh7e9eojWK1EopFixZhxIgRSE9Px4IFCxAQEAAAiIuLw6pVq/Dtt9/i4MGDVT7eggULEBYWhmnTppW7PisrCwDg7u5uUe7u7o7MzMwKj6vVaqHVavlllUpV5ZjqO52Jq5lQMnNCQW0oCCmPyWRCbGwsRCIRfHx8IJVKa/wXGCGNDWMMOp0OGRkZiI2NRXBwMIRC66b5qlZC8cgjj2DdunV4+eWXsXr1av7xhkajgaOjI9avX4+HH364Ssf6+++/sX37duzduxdXrlzhj5OVlYUrV66gY8eOEIu58HQ6ncW+Wq0WEomkwmMvX74c77//fnUurcHQmWsoGNfbAzJKKAgpj06ng8lkgp+fHxQKhb3DIaTecnBwgEQiQVxcHHQ6XbWbLphVe6TMZ555Bo8//jh2797NN8AMCQnBmDFjytQkVCY7OxstW7bESy+9xJfFxcUhMzMTUVFRuHjxIlq2bAmBQICUlBSLfVNSUiqcTwTgalIWLlzIL6tUqkq3b0h0puKEwlRchUuPPAiplLV/bRHSlNji/4lVQ2+7u7tj5syZNTrx1KlTMXXqVIuysLAwDB48GKtWrQIAODk5oWfPnti3bx+/rVarxaFDh/D2229XeGyZTAaZTFaj+Oojxhj0xQlFMz6hoBoKQggh9lfvU/f3338fW7duxYoVK3D8+HFMmTIFjo6OmD17tr1Dq3MGkwEM3JDbCvNcHtTLgxBiI4sXL67xH4tV8fHHH1fYdo40XPUqoQgKCoK3t7dF2ciRI7F3714cO3YMr732GpRKJU6ePAlXV1f7BGlHGqOGfy83MUAoBiQOdoyIENKYqFQq5OTk1Pp58vPzkZ2dXevnIXXLqkcetWXHjh3llo8cORIjR46s42jqH/NMo2CABOAed1CrdUIIIfVAvaqhIJUz11AImBgCgB53ENIIHThwAG3atMHu3bvRp08ftGzZEjdv3gQA7N+/H0OGDEFgYCAGDhxYZhqC69evY8KECQgODkafPn3w1VdfwWQy8et/+ukn9OrVCwEBARg6dCiOHTtWYRx///03AgMDUVhYaFH+ySefWIyY/KCYAODLL79Ehw4d0LlzZ8yfP7/MMUnjUK9qKEjlzKNkikwiroB6eBBSLYwxFOnrfsRMB4moymNgaDQaxMTE4MMPP8Tq1avRunVreHh44M8//8RLL72Er776Cr1798b169fxwgsvQKfT4fnnnwcAjB49GsOGDcNff/2FgoIC/Prrr9i/fz8ee+wxbNy4Ea+++iq+++479O7dG1u2bMGwYcMQFRWFDh06lIlj8ODByMvLQ3h4OKZPn86Xr1+/nm9nUZWY1q5di6VLl+LHH39Ely5d8P333+Pbb7/F8OHDa/qxknpGwBhj9g6itqlUKiiVSuTl5cHFpeHehK9kXsHUv6ZCppcjMvEW0Ooh4Jk99g6LkHpJo9EgNjYWrVq14vvVq3UGhC75u85jufbBCCikVfv7LTw8HE888QQuXbpkMb9R27ZtsWDBAsyZM4cvW7NmDb777jtcvnwZGo0GDg4OOHv2LHr27MlvwxiDQCBAmzZtMH36dHzwwQf8uiFDhsDb2xubN28GAMybNw+JiYkIDw8HADz//PNITk7Gvn37AABnzpxB3759ERsbi4CAgAfGBACtW7fGrFmz8O677/LxdOzYEX5+fjhw4EB1PkZSi8r7/2JW1Xso1VA0IOY2FBJW/KSKHnkQ0ihJJBJ07NiRX87OzsatW7fw3nvv4eOPPwZjDIwxaDQaFBUVAQDkcjnGjh2LSZMm4bnnnsPDDz+Mvn37QiQSobCwEDExMRg4cKDFeQYNGlRh2zUAmD59OoYPH4709HS0aNECv/76KwYMGICAgIAqxVRQUIDY2Fj079+fP6ZAIED//v0RHx9vy4+M1AOUUDQg5nk8JKbiqlMag4KQanGQiHDtgxF2OW913D9MuHkqgW+++QYPPfSQxbalt/vzzz+xd+9e/P3333j22WdhMBgQHh7OT5Nw/zwNUqm0zEjEpQ0ePBheXl7Yvn075syZg99++w0ffvhhlWPS67lZke8f2bgm80WQ+osSigbE3ChTan5KRQkFIdUiEAiq/OihPvHy8oKHhwcuXbqESZMmVbidUCjEmDFjMGbMGJhMJjz66KP49NNPsWXLFjRv3hxRUVEYNGgQv/2FCxcQEhJS4fEEAgGmTp2KX3/9FUFBQcjNzcXEiROrHFOzZs3g4eGBixcvWtRSREdHw8nJqbofA6nnqJdHA2J+5CE3JxT0yIOQJkEgEOCtt97Cl19+iZ07d8JoNEKn0+Gff/7hawySkpKwYMEC/lFCXl4eMjIy4OXlBQB49dVX8fnnn/NzJ+3atQs7d+7E/PnzKz33jBkzcPr0aXz00Ud47LHH0KxZsyrHBABz5szBp59+ips3b4Ixhg0bNuDUqVM2/4yI/TW8VL0JMz/y4BMK6uVBSJPxxhtvQCwWY968eZg2bRrEYjEGDBiAZcuWAQA8PT3RqlUrPPTQQ8jJyYFer8fo0aOxdOlSAMDbb7+N9PR09OrVi6upUSiwatUqDBkypNLzdu7cGZ06dcKpU6fwxx9/VCsmgBt9MyYmBh06dICjoyO6dOmCyZMnIzc316afD7E/6uXRgGy5vgXLzy5HWL4Uv2TeAcZ+C3SdYe+wCKmXKmu1Xp9ptVrk5eWhRYsWFW6Tm5sLZ2dniETlt83Iz8+Ho6NjuRM+GQwGqFQqNGvWrExX1vz8fBiNxjIjEefl5aGoqAjNmzev8JwPikmtVkMsFkMqlaKgoAB6vZ6v7SD2R708mhhzGwpHZuAK6JEHIY2OTCarNJkA8MCpB5ydnStcJxaL4ebmVq39lEollMrK22w9KKbSU8hT+4nGidpQNCDmga0czROD0SMPQggh9QQlFA2IuYbCiXFdsaiXByGEkPqCEooGxNzLw8WcUNAjD0IIIfUEJRQNiDmhcKYaCkIIIfUMJRQNiLo4oZDROBSEEELqGUooGpBCHTc+vgNjYGIHQEzD1xJCCKkfKKFoQNT6khoKAfXwIIQQUo9QQtGAFJV+5EGPOwghhNQjlFA0IEWG4kceJhM1yCSkkbp27RqeffZZe4dRY3fv3sWzzz7LzzhKGj9KKBoQTfHAVjLGaFArQhqp5ORkbNy40d5h1FhGRgY2btwIo9Fo71BIHaGEogHRlZ4cjB55EEIIqUdoLo8GRGsqNX05PfIgpEk4fvw4du/ejc8//5wvu337Nj7++GOsXbsWEokEFy5cwI8//oi3334bv//+O2JjY9G1a1fMmjULcXFx2LhxI7KysvDwww9j/Pjx/HHM+/3vf//D5s2bERcXhx49euDZZ5+1mOTLYDBg+/btOHXqFJycnDBs2DAMHTrUIs7s7GysWbMGKSkp6NixI9q3b1/7Hw6pV6iGogHRm+iRByFNza1bt7Bjxw6LsrS0NIvHCfHx8Vi/fj2GDRsGnU6HwMBAvP3223jiiScwcuRIyGQyeHt749lnn8X69ev545j369evH/Lz8xEYGIj//e9/eP755/lt9Ho9RowYga+++grBwcFwdXXFM888gw8++IDfRq1Wo0+fPti7dy+Cg4Nx8uRJTJkypZY/GVLfUA1FA2Iw6QAAchMDZFRDQUi1MQbo1XV/XokCuG+qcFszGAzYsGEDBg4cCADQ6XRYsmQJzp49i549ewLgpqHeuHEjZs2aZbHfBx98wDcEHTJkCHr16oX58+cjLCwM33//PVJTUxEVFQWplBv7ZujQoejXrx/mzJmD5s2bY82aNVCr1bh06RI/9fWUKVPw22+/1eo1k/qFEooGwmgywgiutTQ98iDESno1sMyn7s+7OBmQOtbqKSQSCfr3788vt2nTBo6OjnwyYS67v7YDgMVjkB49eiAwMBAnTpxAWFgY/vrrL5hMJsydOxeMMTDGYDKZYDAYcPXqVQwePBjHjh3DY489xicTADBp0iRKKJoYuycUGo0GUVFRKCgoQPv27eHr62uxPioqCnfv3rUoUyqVGDZsWF2GaXfa4gaZAD3yIISUJZVKIRSWPMUWCoUWN3hz2f29LiQSCZydnS3K3NzckJGRAQDIysqCv78/BgwYYLHN4MGDERQUBIDr0dGhQ4cyxyBNi10Tip9//hnvv/8+/P39IRaLERERgRkzZuD777/n/2P88MMP2Lt3L/r06cPv5+/v36QTCurlQYiVJAqutsAe57WSTCaDTqezKMvJyalpRDy9Xo+UlBR4e3sDABhjSEhIgJ+fHwDA19cXarW60rEx/P39ER8fb1EWFxdnsxhJw2DXhMLR0RHR0dFQKrnq++joaHTt2hVjxozB448/zm/Xp0+fcqvpmhJzQiExMa4lLT3yIKT6BIJaf/Rga8HBwUhNTcWdO3cQFBQEk8mEDRs22PQcX375JT799FMAwC+//AKVSoXHHnsMAPD0009j/Pjx2L9/Px599FEAXNLx22+/8Q0vn3zyScyaNQuxsbFo1aoVNBoNvvvuO5vGSOo/uyYUEydOtFgOCQmBWCxGXl6eRXleXh727dsHpVKJTp06wcWl6f11bh4lk59plB55ENIk9OnTB6NHj0b//v0xaNAgXL16FS1btrTZ8WUyGU6dOoUBAwbA2dkZR44cwWeffQYfH66tyRNPPIGlS5fiiSeeQO/eveHu7o5Lly6hV69efEJhbi/RvXt3DBo0CFeuXIG/v7/NYiQNg93bUKSlpeHEiRNQqVTYtGkThg8fjgkTJlhsc/78eaxatQpJSUlITk7GV199haeeeqrCY2q1Wmi1JY8IVCpVrcVfV7SlB7UCAIWHHaMhhNSWDh06lKmBCA8Px8mTJ5Geno4uXbrAxcUF+/btg0QiAQB069YNP/zwg8U+vXv3xjfffGNRNmjQILi6ulqUicViHDt2DOfOnUNcXBxWr16NkJAQi22WLFmCF198ERERERAIBOjatSsCAgL49QKBAH/++SdOnjyJ1NRUdOjQAR4eHvjrr7/4GEnjJ2DMfIeyj+joaHz00UfIzMzEtWvXsHDhQrz99tsQFHexOnLkCPr06QOFgnsG+fnnn2Px4sW4cOECOnbsWO4xly5divfff79MeV5eXoOt3YhOj8ZT+5+Cr16P/YkpwP8yaPpyQiqh0Wj4Kvj7GycSTnh4OGbMmIGCggJ7h0LsrLL/LyqVCkql8oH3ULsPbBUWFoYdO3bg6NGjOHLkCD788EP89NNP/PpHHnmETyYA4PXXX4eLiwv2799f4TEXLVqEvLw8/pWQkFCr11AXcrW5AAClyQSD1IWSCUIIIfWK3ROK0jp27IgePXrg2LFjFW4jEAjg7OyM9PT0CreRyWRwcXGxeDV0ORquVXczowkmetxBCLGB8h6VEGItuyUURqMR2dnZFmWFhYW4ffs2PxaFyWRCZmamxTZRUVH8ePNNSbaG+6zcjEZA0dzO0RBCGgN/f39Mnz7d3mGQRsJujTL1ej369++PsWPHIjQ0FLm5ufj5558hl8vx6quvAuCSjgEDBmDMmDHo0KED4uPjsXr1aowYMaJMw83GLrOISyiamUwQOFENBSGEkPrFbjUUcrkcZ8+ehY+PD/7991/cunULL7/8Mq5duwYvLy8A3Ahu586dg6+vL44ePYqcnBysXbsW+/fvt5gJrynIVHMJhavRBJET1VAQQgipX+zabdTZ2ZmvjajJNk1BVlHJIw+hUws7R0MIIYRYqleNMknFcjS5ALhHHnCkGgpCCCH1CyUUDUSuztzLwwg4uts5GkIIIcQSJRQNhEqXC4DrNko1FIQ0Xmlpadi2bZu9w6ixrKwsbNu2DSaTqcbHqspnkpSUhD/++KPG56qqCxcuYMeOHThz5gwAboTmI0eO4LfffkNiYmKdxVGfUELRAGiNWuhM3FwezUxGSigIacQuX76MqVOn2juMGrt9+zamTp1aZqZUa9z/maSmpmL79u0W25w7dw7PPPNMtY6bk5ODbdu2lZk/6kEWLFiA0aNHY/v27bhw4QK0Wi26d++ON998Ezt37kRSUlK1jtdY2H0uD/Jg5kGtxIzB2cQooSCENCleXl6YPHkyvxwdHY1p06Zh0qRJNTru3bt3MXXqVFy+fJmf9boq1q5dix07dvCzrx45cgR3795FXl4epNKmO4oxJRQNQMkomUYAAsChmX0DIoTUmZiYGFy9ehWjR4/myzIzM3Ho0CFMmjQJQqEQSUlJOH36NJ544glcuXIFsbGx6Ny5M1q1agWTyYQzZ84gKysLPXr04LvlA7DY79KlS4iLi0NYWJjFxF9mqampOHfuHJycnNCtW7dyb8Dnz59HSkoKOnTo8MDr+uOPP9CnTx9+5tTTp08jOTkZTz75JABuivTt27djyJAhaN68OcaNGweAq1U4duwYGGP8Y5DSk5kxxnDlyhXExcUhNDQUrVu3rsKnXLVrTU9Px/79+6FWq3H+/Hm+ZiMiIgISiQR//vknAGDKlCn8Z/vkk0/izJkzSEhIwPDhw/ljXb58GXfv3oWfnx/CwsLKDIUQFxeHK1euwN3dHd26dbNIVIxGI86cOYO0tDQEBweXmdfqzJkzYIyhc+fOiIqKQkFBAXr37l1mYjhbo4SiASg97LZB7gaJsGmNwUFIU3bkyBF89NFHFgnFjRs3MHXqVIwbNw5yuRznzp3D9OnT0bVrV4hEIkgkEpw4cQJffPEFtmzZAqlUCsYYLl68iIMHD6JPnz4AwO/Xr18/fuKn//77D9988w1eeOEF/nzLly/HihUr0KdPHxQVFeHatWvYvHkzRowYwW8zY8YM7NmzB/369cP169cRGhpa6XV9/vnnGD16NBYvXgwAmDNnDqKjo5GUlAQfHx9cunQJU6dORUZGBv/IY8qUKcjJyUFERAQYYwgPDwcAjBgxAs2aNYPBYMDIkSOhUqng4uKCo0ePlrmWB6nsWjMzM7Fnzx4AwLFjx3DlyhUAQHx8PLRaLR/PlClT+M924MCByMvLQ+vWrdG3b18IhUJMmjQJV65cQdeuXXHr1i04Oztjz549fLL30Ucf4dNPP8WAAQNQUFAAlUqFP/74A23atEFaWhpGjhyJrKwsdOjQAadPn8awYcOwdetWPilZs2YNLl68CI1Gg4CAAKSmpiIlJQXHjh1Du3btqvxZVBtrAvLy8hgAlpeXZ+9QrLLr9m7W8eeObNaaIKb7doC9wyGkQSgqKmLXrl1jRUVFfJnJZGKFusI6f5lMpirH/c8//7DSv5rXrl3LAgICLLY5ceIEA8Bf286dOxkA9vXXX/PbvPTSSwwAW7t2LV82ffp0Nnr0aH7ZvN/8+fMtzqdQKFhqaipjjLF9+/YxDw8PFhsby2/z888/sxYtWrDCwkLGGGN//PEHc3BwYDdv3mSMMVZYWMi6detmEeP93nnnHTZs2DDGGGM5OTlMIpGwrl27sl9//ZUxxtiqVatY586dy/1M9u/fz0QikcXxzNeyatUqvmzVqlXMzc2t3PMzxti5c+cYAHb58uUqX2t+fj4DwCIiIvhttm7dytzd3cuN58MPP7Qof/HFF9mjjz7KtFotY4wxo9HInnzySTZ9+nTGGGM6nY5JpVJ28OBBfp9bt27xMT777LOse/furKCggDHGWExMDHNxcWE//PADv/0zzzzDHBwc2LVr1xhj3M/9I488wp577rkKP4vy/r+YVfUeSjUUDUBcLjcRWjOjESL/pjWHCSG2VGQoQu8tvev8vGemnYFConjwhjUgFArx/PPP88t9+/bFTz/9hJkzZ1qUff3112X2ffPNN/n3zz33HBYvXoy9e/di1qxZ2LBhA0JDQxEZGYlz586BMQaRSIT09HRcvXoVPXv2xPbt2zF27Fj+0YNCocC8efPw3HPPVRjv4MGD8dVXX0Gn0+HYsWPo3LkzRo8ejX///RfTpk3D0aNHMXjw4Gp/Bi+99JLFORYsWIDMzEx4eDx4yoKqXGt1lR6Y0WAwYPPmzZg7dy52794NxhgYY/D19cWOHTv4a5DJZLh8+TKGDBkCoVCI4OBgACWPgb7//ns4OjoCAFq1aoVp06Zh27ZtePHFFy2uvX379gC4STUfeugh/PPPP9WOvzoooWgAElQZALhHHkJfSigIIWU5ODhALpfzy+ZZl0s/m5fJZNBoNBb7SSQSvh0DwN3Q/P39ERcXBwC4d+8eCgoK+Bue2eTJkyEWc7eQ+Ph4DBo0yGJ9q1atKo13wIAB0Ov1OHv2LJ88DB48GLNmzQJjDCdOnMDatWur8QmU/xkAKHPNFanKtVaHQqGwmO06PT0darUa0dHRSEhIsNjWnDyJRCL88ssvWLhwIZYvX45BgwZhypQpmDBhAjIyMqBWq8u0C2nTpg0OHjxoUebm5maxXN733tYooWgArmZGAQBa6/UAJRSEWM1B7IAz087Y5bzWEgqFZcZysOWNQa/Xo6CgAE5OTnxZTk4O/xe9i4sLgoKCsGXLlgqP4e7ujpycHIuy+5fv5+joiO7du+Po0aM4evQoPvroI/Tp0wdJSUnYt28fsrOz8dBDD9XgyqqvKtdaHQKBwGLZ2dkZAoEAL7zwQqU9VMaOHYuxY8fi9u3b2L9/P55//nnEx8fjlVdegUgkKjNTd3Z2dpVqYGobjUNRz2VrspFQdB0A0EcjBNyD7RwRIQ2XQCCAQqKo89f9N5bqaNmyJdLS0pCfn8+X/fvvv7b4OHjmxoQA1yUzNjYW/fv3BwCMHDkSe/bsKTO2QunlAQMGYN++fRZjTph7PFRm8ODB2LlzJ65evYqBAwdCJpOhT58+WLp0KTp16gR39/JHBXZycoLRaIRer6/OZT5QVa61JpydndGvXz/88MMPYIyVe46ioiI+GQsODsarr76KsWPH4vTp05BIJOjVq5fFZ2s0GrFz504MGDDAJjHWBNVQ1HNH4o6BgaGdVgcHZSggpByQkKZk0KBB8Pb2xoQJEzB58mRcvHixTJV8TYhEIrzzzju4d+8enJ2d8fnnn2PSpEno3r07AGDevHnYs2cPevfujZdffhnu7u6IiorC4cOHcevWLQDAyy+/jDVr1mDYsGGYNm0aIiMjceDAgQeee/Dgwfjkk0/Qs2dP/tHA4MGDsXTp0konhQwNDYVCocDbb7+NXr16WXQbrYmqXGtNfffddxgyZAgeeeQRTJw4ERqNBocPH0br1q3x9ddfIz8/H/369cOYMWPQqVMnJCYm4o8//sBPP/0EAFi5ciUeeeQRiEQi9OrVC9u3b0d+fj4WLVpkk/hqgu5O9ZjBZMDmaG40uMHqIrh3GGLniAghte3+QZzkcjkiIiLQo0cPREREICQkBPv378fkyZP59hG+vr6YMGGCxXECAgL4MR3M2rRpY9H91Hz8EydOQK/XIyoqCm+//TZ++eUXi/WHDx/GihUrcO/ePVy4cAFdu3bFxYsX+W2cnJxw5swZDBw4EOfOnUNoaCgOHz5sEWN5BgwYgClTpmDOnDl82dixYzF58mRMnDixws/Ezc0Nhw4dgsFgwK5du3Dx4sVyPwMXFxdMnjwZCkX5DWLNNRzmdhdVuVaJRILJkydbPGIo77MuLx4A6Ny5M65evYpHH30UZ8+eRXJyMl599VW+sWyLFi0QGRkJf39/nDx5Enl5eTh48CD/iKRv3744f/483NzccPLkSQwePBhRUVEW8fTu3Rt9+/a1OG9oaKhFN9/aIGD317s0QiqVCkqlku9nXVMno/Zg5/lvAAAM5o+v9MfIiv9xZYLi7Up/BQOYgAHsvv2KvxrAcFdYhASJDgLG8FO2GD1ePgnInEAIeTCNRoPY2Fi0atXKoqEeKREeHo4ZM2agoKDA3qHYxa+//oqnn34ahYWFTf5npLL/L1W9h9IjDyvcSY7GQUlynZxLYTLhzWwdOk3fQckEIYTYQHJyMnbv3o1PP/0U06ZNa/LJhK1QQmGFoJZdMSL1v+IlAf9VAAEYAK79laDU+lINsgQl25n34rcTCMCKSwQQwkvshIdaD0L3PpMACf3AE0Jsq6Jq+cYuLS0NkZGReOONN6o1iiapHD3yIIQ0SvTIg5Cqs8UjD2qUSQghhJAao4SCEEIIITVGCQUhpFFrAk91CakxW/w/oYSCENIoSSQSAIBarbZzJITUf+b/J+b/N9agXh6EkEZJJBLB1dUV6encbL0KRc2GwCakMWKMQa1WIz09Ha6urpUORPYglFAQQhotLy8vAOCTCkJI+VxdXfn/L9aihIIQ0mgJBAJ4e3ujRYsWNp9IipDGQiKR1KhmwowSCkJIoycSiWzyC5MQUjFqlEkIIYSQGqOEghBCCCE11iQeeZj716pUKjtHQgghhDQs5nvng8aqaBIJRX5+PgDAz8/PzpEQQgghDVN+fj6USmWF65vE5GAmkwnJyclwdna2WT90lUoFPz8/JCQkNJoJxxrbNTW26wHomhoKuqaGobFdU21dD2MM+fn58PHxgVBYcUuJJlFDIRQK4evrWyvHdnFxaRQ/iKU1tmtqbNcD0DU1FHRNDUNju6bauJ7KaibMqFEmIYQQQmqMEgpCCCGE1BglFFaSyWR47733IJPJ7B2KzTS2a2ps1wPQNTUUdE0NQ2O7JntfT5NolEkIIYSQ2kU1FIQQQgipMUooCCGEEFJjlFAQQgghpMYoobBCfn4+IiMjce/ePXuHYhW1Wo3o6GgkJSVVuA1jDFevXsXFixdhNBrrMLqauXbtGk6ePAmdTldmnV6vR1RUFK5fv26HyKyTlJSE6Ojocq8HAAoLC3H+/HnExsbWcWTWSUlJwfnz5xEXF1fhNllZWTh37hxSU1PrMLKqYYzhwoULiI6OrnAbvV6P6OhoXL9+vcKhiquyTV0xGAw4ffo0bt68WeE26enpiI6ORl5eXoXbmH8WY2JiaiPMaiksLMSpU6cQHx9f6XaMMfz3338Vfj9zcnJw7tw5JCcn10KU1ZOZmYmTJ08iMzOzwm0YY7h+/Xql38vU1FScO3cOWVlZtg+SkWpZu3YtUygUrG3btszR0ZE9+uijrKCgwN5hVUlKSgp75plnmFKpZGFhYaxZs2asf//+LDY21mK769evs5CQEObp6cl8fX2Zr68vO3PmjH2CroYLFy4wBwcHBoAlJCRYrDty5Ajz9PRkgYGBzN3dnXXu3Jndu3fPTpE+WFJSEhsyZAhTKpWsR48erFWrVmzfvn0W22zatIk5OTmxkJAQ5uTkxIYOHcry8vLsFHHlsrOz2ZAhQ5izszPr3r07c3NzY926dSvzs7dkyRImk8lYaGgok8lkbNasWcxoNNon6Pt88cUXLCQkhLm6urIuXbqUu82xY8eYl5cXCwgIYB4eHqxjx44sJiam2tvUhYKCAvZ///d/zM/Pjzk7O7PJkyeX2ebkyZOsX79+rEWLFiwsLIw5ODiwl156iRkMBovtfv31V+bs7MxCQkKYs7Mze/jhh1lubm5dXQovOTmZvfLKK8zb25vJZDL23nvvVbr9ihUrmFAoLPf7+fHHH/M/i3K5nD311FNMr9fXTuCVuHbtGnvqqaeYt7c3A8C2bt1a7naHDx9mgYGBzM/Pj4WFhbE+ffqwuLg4fr3BYGCzZs1icrmc//+1ZMkSm8ZKCUU1XLp0iQmFQrZlyxbGGGPp6eksMDCQvfLKK3aOrGpOnz7NNm7cyP+nyM/PZw899BAbOHAgv43JZGKdO3dm48aN43+Rz5o1i/n5+TGNRmOXuKsiPz+ftW3blr355ptlEoq8vDzm7u7O3nrrLcYYYzqdjj388MMW112f6HQ61rlzZzZ8+HCWn5/PGGMsKyvL4hfJjRs3mFgsZj/99BNjjLthBwcHsxdeeMEuMT/IwoULmZ+fH8vKymKMMaZWq1m3bt3Y+PHj+W3Cw8OZRCJhp06dYoxx16hUKtnq1avtEnNpBoOBLViwgN24cYO9/vrr5d6A8vPzWfPmzdnChQsZY4zp9Xo2dOhQ1rdv32ptU1diY2PZ0qVLWVJSEnvsscfKTSg2bNjA/vvvP3752rVrzNXVla1YsYIvu337NpNIJOzHH39kjDGWk5PD2rZty2bOnFn7F3GfU6dOsdWrV7OcnBzWpk2bShOKM2fOsICAADZjxowy388DBw4wkUjE/v33X8YYY3fv3mVubm4W111X/vjjD7Zx40aWn59fYUJx5coVJpPJ2LJly/iy8+fPs4iICH551apVzNXVld28eZMxxtiJEyeYWCxmu3btslmslFBUw8KFC1lQUJBF2SeffMKUSmWZjL2hWLduHROLxXzycPbsWQaARUZG8tvcu3ePAWB79uyxV5gP9NRTT7EFCxawf/75p0xCsWnTJiaRSFhOTg5fduDAAQaA3b592w7RVu7XX39lAoGg0hqUxYsXM19fX4uyVatWMYVCUS8Tv6eeeooNHTrUomzWrFnsoYce4pfHjBnDRo4cabHN888/X2FtgL1UlFBs2bKFiUQilpmZyZcdOnSIAWDXr1+v8jb2UFFCUZ5x48ax0aNH88tLlixh3t7ezGQy8WXffPMNk8vlTK1W2zzWqqosocjNzWVt2rRhBw8eZPPnzy/z/Zw0aRIbPHiwRdm8efNY27ZtaynaB9Pr9RUmFFOnTn3g/5POnTuz2bNnW5QNHTqUjR071mYxUhuKaoiKikL37t0tynr16oW8vLx68dzQGufOnUOrVq34CV+ioqIgFArRtWtXfpuAgAC0aNECUVFR9gqzUps3b0Z0dDQ++eSTctdHRUWhdevWcHV15ct69erFr6tvDh8+jM6dO8Pf3x9Xr17FjRs3oNfrLbap6GdRrVbj1q1bdRlulSxcuBBXr17FsmXLcPjwYXz99dfYu3cvlixZwm9T0TVduXKlzPXXR1FRUQgMDIS7uztfdv/PWVW2qc/0ej0uXryIoKAgviwqKgrdunWzmHixV69e0Gg0uHHjhj3CfKAXX3wRY8aMwbBhw8pdX9HP4q1bt6BWq+sixGo5fPgwHn/8cRQVFeH8+fNISEiwWK/X63H16tVyr8mWP3dNYnIwW8nOzkb79u0tysy/GLKzs+0RUo0cPnwY69atw4YNG/iy7OxsuLq6lplRzt3dvV5e4+3bt/Haa6/h33//rXB0uOzsbItf4AD4a6yP15ScnAxHR0f069cPKpUKarUaWq0W69atw6hRowBw19ShQweL/erzz2JoaCieeuopfPrpp2jTpg1iY2Mxbtw4/mYKlP99cnd3h9FohEqlKrOuvikvfmdnZ0gkEv57UpVt6rNFixYhKysLr776Kl+WnZ2NNm3aWGxXn38W165dixs3bmDTpk0VblPRzyJjDDk5OVAoFLUdZpUZjUakp6cjISEBbdu2hbu7O+Li4tCuXTts27YN/v7+yMvLg9FoLPeabPk9ohqKapBIJNBoNBZlRUVFAACpVGqPkKx27tw5PPnkk3jjjTfw1FNP8eXlXSPAXWd9vMZZs2Zh1KhRyM3NxcmTJ3HlyhUA3PWZez6Ud006nQ4mk6leXpNEIsF///2HV199FVevXkVsbCyefvppTJ06FTk5Ofw2Delncf78+QgPD8edO3f4v6Du3r2LqVOn8ts0tGu6X3nxGwwGGAwGPv6qbFNfrVy5Et9++y127NiBwMBAvrwhfd9SU1OxYMECzJkzB+fOncPJkyeRnJyMwsJCnDx5ku/F0pCuSSQSQSgU4q+//sLx48cRFRWF+Ph4mEwmvPTSSwC46wFQ7jXZ8noooaiGgICAMl0tzcv+/v72CMkqkZGRGD58OF544YUyjwkCAgKgVquRm5vLlxkMBqSnp9fLa/T09MTdu3fxzjvv4J133sH69esBAMuXL8eePXsANLzvW2BgIFxcXCxuti+99BJUKhUuXrwIoOFd0969ezFlyhR4eHgAABwdHfHss89i//79/OOMiq7J1dUVzs7OdR5zdQUEBCA5OdmiG6h52fw9qco29dGqVavwv//9Dzt37izzmKAh/SwWFRWha9eu2Lx5M/8749SpU0hJScE777yD27dvA6j4mhQKRb2sKQsICMCIESP4RM/JyQnTp0/HiRMnAHBTj7u6upZ7Tbb8HlFCUQ3Dhg3D8ePHLfpi79q1C127dq2XP2TlOX/+PIYNG4bnnnsOn3/+eZn1gwcPhkQiwe7du/myf/75B2q1GkOHDq3LUKvk999/x8mTJ/nXl19+CQD4888/+WrZYcOGIS0tDWfPnuX327VrF5ycnNC3b1+7xF2ZESNGlEnqEhMTAQDNmzcHwF3Tf//9Z9EnfdeuXWjfvj1atmxZp/FWRfPmzflrMEtISICrqyv/19OwYcOwb98+i3FPdu3aVeFz7vpm2LBhyMzMREREBF+2a9cuKBQK9O/fv8rb1DdfffUVFi1ahD///BMjR44ss37YsGE4c+YM0tPT+bJdu3YhODgYAQEBdRnqA7Vq1cri98XJkycxceJEBAUF4eTJk+jRowcA7poOHDhg0XZn165dGDJkSJnHwfXBiBEjyiQLiYmJ/O8LABg6dCj/RxbA/aH4119/2fb/l82adzYBRUVFLDQ0lD300EMsPDycLVmyhIlEIrZ//357h1Yl169fZ25ubuzhhx9mJ06csHjpdDp+u0WLFrFmzZqx9evXs82bN7OWLVuy5557zo6RV115vTwYY+yJJ55gQUFBbNu2bey7775jjo6O7NNPP7VTlJUzmUxsyJAhbNCgQWz37t1s69atLCQkhD3++OP8NlqtlnXp0oX169eP7dy5k33wwQdMJBKx8PBwO0ZesQ0bNjCRSMSWLl3KDh48yL744gvm5ORk0Qo/OTmZtWjRgk2YMIHt3r2bvfTSS0yhULDLly/bL/BSoqKi2IkTJ9iUKVNYUFAQ/3+ndA+viRMnstatW7OtW7ey77//njk5OVl05avqNnXl1KlT7MSJE6xfv35syJAh7MSJExZjzqxbt44BYO+++67F74vo6Gh+G71ez7p168b69OnD/vzzT/bxxx8zkUjEduzYUefXo9Fo+BhbtmzJZs6cyU6cOMEuXbpU4T7l9fLIyMhgPj4+bOzYsWz37t1s3rx5TC6Xs/Pnz9fyFZSVk5PDTpw4wY4ePcoAsKVLl7ITJ06wO3fu8NvExcUxDw8PtmDBAvb333+zzz//nDk4OLA1a9bw21y6dIkpFAo2e/Zstnv3bjZ+/HjWokULlpKSYrNYabbRasrMzMSKFSsQFRUFd3d3zJ49Gw8//LC9w6qS/fv34+OPPy533e7du+Hm5gaAG21t/fr1CA8Ph8FgwMiRIzFv3jyIxfW/DW9kZCQWLFiAnTt3WmTnWq0Wq1atwpEjRyCTyTB58mRMnz7djpFWrqioCF9++SWOHTsGR0dHDB48GLNnz7Z43pmTk4MVK1YgMjISbm5ueOGFF+r1X/P//PMPfv31VyQlJaFFixYYN24cJkyYYNE7IDY2Fp9++ilu3rwJf39/vPbaa+jSpYsdoy4xc+ZMvkq8tAMHDsDJyQkA1zZn9erVOHToEGQyGSZOnGjRRqmq29SVIUOGQKvVWpS5u7tj165dAID/+7//w7///ltmv3bt2mHdunX8cm5uLlasWIFz586hWbNmeP755zFixIjaDb4cqampmDBhQpnyrl274uuvvy53n6+++gpRUVEWjdMBID4+HitWrMCNGzfQsmVLzJ8/v0wvibpw5swZvP7662XKx40bhzfeeINfjomJweeff46bN2/Cx8cHM2bMKPM9uHjxIr788kvEx8ejbdu2eOutt9CqVSubxUoJBSGEEEJqrP49DCKEEEJIg0MJBSGEEEJqjBIKQgghhNQYJRSEEEIIqTFKKAghhBBSY5RQEEIIIaTGKKEghBBCSI3V/5GKCCF1KjU1FcePH8ekSZPq/NwxMTH8EOkdOnRAp06drD6WNdeRmZmJQ4cOAQD8/Pzq7XDYhNRHlFAQ0kRkZWXhn3/+qXSb7t274+7du5g2bZpdEoojR45g/vz5GD16NIRCYY0Siujo6GpfR1ZWFsLDwxEZGYmwsDBKKAipBkooCGkicnJyEB4ezi9HRkYiOzsbw4cP58vc3Nzg7e2NyZMn2yFCTvPmzbFt27YaH8ea62jbti22bduG2bNnW0y8Rgh5MEooCGkigoKCLG7Us2fPRmRkZJmbd2pqKsaOHcsvJyYm4uzZsxg3bhwuX76MuLg4dO7cGYGBgTAajThz5gyys7PRs2dPeHp6ljlvSkoKIiMj4ezsjG7dusHFxaVacVt7fk9Pz3Kv44knnsDly5cRHx+PDh062HQuA0KaMkooCCEW7n9UcPr0aTzzzDPo2LEjpFIphEIhTp06hVWrVmHTpk1wcHCA0WjE5cuXcejQIfTs2ZM/1ocffoiVK1eiT58+UKvVuHHjBrZs2YKhQ4dWOR5rz1/RdfTr1w+FhYVwcnLCsWPH8P3332PmzJk2/AQJaZqolwch5IHUajVmzpyJEydO4NixY3j22WfxyiuvYO7cuTh27BhOnjyJRx99FMuXL+f32bNnD7755htcunQJBw4cwPHjx7F8+XI89dRT0Gg0tX7+io4zbtw4/Pfffzh48CCWLVuGd955x6rPhBBiiRIKQsgDiUQizJo1i1/u27cv5HI5nn76aYuyW7du8csbNmxAhw4dcPbsWfz+++/Yvn07xGIxUlNTcf369Vo/f0XHefHFF/nlwYMHIz09Hbm5udWKhxBSFj3yIIQ8kJOTEyQSCb8sk8mgVCohFAotykrXPNy7dw9arRY7duywONbkyZMt9qut81f1OACqXWNCCCmLEgpCSK1wcXGBv78/Nm3aZO9QCCF1gB55EEJqxciRI7Fr1y6kpKRYlCclJdkpIkJIbaIaCkJIrZg/fz727t2L3r17Y+7cuXBzc8OFCxdw/PhxXLt2zd7hEUJsjGooCGmievbsiREjRpQpv39AKD8/P4wfP95im8DAQIwbN86iLCgoCI8//ji/7ODggH///RfLli1DbGwsLly4gB49euDChQuVxlVYWIht27bh8uXLNTp/Va7D1dUVkydPhoODAwBu6O1t27bh7t27lcZICClLwBhj9g6CEEIA4N9//8UPP/wAAHjyySfrfPjvmzdv4r333gMA9OnTBwsWLKjT8xPSkFFCQQghhJAao0cehBBCCKkxSigIIYQQUmOUUBBCCCGkxiihIIQQQkiNUUJBCCGEkBqjhIIQQgghNUYJBSGEEEJqjBIKQgghhNQYJRSEEEIIqTFKKAghhBBSY5RQEEIIIaTG/h9uQFS8KR5a2AAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "import matplotlib.pyplot as plt\n", "import numpy as np\n", "import pybamm\n", "\n", "import pytes\n", "\n", "variants = {\n", " \"resolved\": {},\n", " \"lumped\": {\"storage element\": \"lumped\"},\n", " \"lumped with Jeffreson\": {\"storage element\": \"lumped with Jeffreson\"},\n", "}\n", "\n", "\n", "def charge(options, conductivity=None):\n", " \"\"\"The default charge (Nallusamy et al., 2007), optionally with other capsules\"\"\"\n", " model = pytes.PackedBedModel(options=options)\n", " parameter_values = model.default_parameter_values\n", " if conductivity is not None:\n", " for phase in (\"solid\", \"liquid\"):\n", " parameter_values[f\"PCM {phase} phase conductivity [W.m-1.K-1]\"] = (\n", " conductivity\n", " )\n", " t_end = 160 * 60\n", " return pybamm.Simulation(model, parameter_values=parameter_values).solve(\n", " [0, t_end], t_interp=np.linspace(0, t_end, 321)\n", " )\n", "\n", "\n", "def compare(conductivity=None):\n", " fig, ax = plt.subplots(figsize=(6, 4))\n", " for label, options in variants.items():\n", " solution = charge(options, conductivity)\n", " print(f\"{label}: solve time {solution.solve_time}\")\n", " ax.plot(\n", " solution[\"Time [min]\"].entries,\n", " solution[\"Outlet temperature [degC]\"].entries,\n", " label=label,\n", " )\n", " ax.set_xlabel(\"Time [min]\")\n", " ax.set_ylabel(\"Outlet temperature [°C]\")\n", " ax.legend()\n", " plt.show()\n", "\n", "\n", "values = pytes.PackedBedModel().default_parameter_values\n", "h = values[\"Heat transfer coefficient [W.m-2.K-1]\"]\n", "R = values[\"Capsule radius [m]\"]\n", "for phase in (\"solid\", \"liquid\"):\n", " k = values[f\"PCM {phase} phase conductivity [W.m-1.K-1]\"]\n", " print(f\"Biot number, {phase}: {h * R / k:.1f}\")\n", "compare()" ] }, { "cell_type": "markdown", "id": "47a82a75", "metadata": {}, "source": [ "The paraffin capsules of the default set have Biot numbers of 7 to 18: far too large to\n", "lump. The lumped capsule takes heat much faster than the resolved one, as nothing slows the\n", "heat on its way into the capsule: its outlet stays colder while the PCM melts, and the bed\n", "is fully charged after about 105 minutes instead of more than 160. Jeffreson's correction\n", "follows the resolved outlet for the first half hour, but then charges too slowly: one\n", "effective coefficient cannot follow a melting front moving into the capsule.\n", "\n", "With capsules a hundred times more conductive (`k` = 40 W.m-1.K-1, `Bi` about 0.07), as\n", "with a conductive filler, the three agree closely, and the lumped capsules solve hundreds\n", "of times faster." ] }, { "cell_type": "code", "execution_count": 2, "id": "2b60f4e2", "metadata": { "execution": { "iopub.execute_input": "2026-11-01T11:16:51.998432Z", "iopub.status.busy": "2026-11-01T11:16:51.998165Z", "iopub.status.idle": "2026-11-01T11:16:59.135377Z", "shell.execute_reply": "2026-11-01T11:16:59.134472Z" } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "resolved: solve time 6.849 s\n", "lumped: solve time 11.158 ms\n", "lumped with Jeffreson: solve time 10.678 ms\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAhQAAAFzCAYAAABmY5CaAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjExLjIsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvgI3uAAAAAAlwSFlzAAAPYQAAD2EBqD+naQAAZEdJREFUeJzt3Xd8U1UbB/Bfkqbp3lAKdABtgbLK3kumKHsPByAyREBwgb4IDhAVGQ5UQEQREFHKEBDZBcootOwllAItpXTvkeS8f4RGQlto0rS34/f1E+k99+Te56RN8/Tcc8+RCSEEiIiIiIpBLnUAREREVP4xoSAiIqJiY0JBRERExcaEgoiIiIqNCQUREREVGxMKIiIiKjYmFERERFRsFlIHUBq0Wi2io6Nhb28PmUwmdThERETlhhACqampqF69OuTywvshKkVCER0dDU9PT6nDICIiKrfu3LmDmjVrFrq/UiQU9vb2AHQvhoODg8TREBERlR8pKSnw9PTUf5YWplIkFHmXORwcHJhQEBERmeBpQwY4KJOIiIiKjQkFERERFRsTCiIiIio2JhRERERUbEwoiIiIqNiYUBAREVGxSXrbaEhICDQaTb7yatWqwdfX16Dszp07uH//Pvz9/XnrJxERURkjaULxwQcfICMjQ7+dk5ODU6dOYfbs2ViwYAEAICsrC6NHj8auXbvg7e2NyMhILFq0CK+//rpUYRMREdFjJE0o9uzZY7C9YcMGjBo1Ci+99JK+bP78+Th58iRu3LgBDw8PBAUFYeDAgWjVqhVat25d2iETERFRAcrUGIrVq1ejU6dOqFu3rr5szZo1eOWVV+Dh4QEAGDBgABo2bIg1a9ZIFSYRERE9psxMvX3r1i3s378fP//8s74sOjoa9+/fR/PmzQ3qtmrVCmFhYYUeKzs7G9nZ2frtlJQU8wdMRFSWCAHkZiA3PQ4pqfcQm/YASWnJSM1KQ3ZuBrJyM5Gdm4FsdRZy1FnI1uRCrdVAK7TQCA2EEBDQQCsEtEIDAQEhBLRCC63QQsi0um0IiPwnN/inkL1PLDEsFfnKZQU8raBIinKeiq5J9c54sc+cUj9vmUkofvzxRzg5OWHIkCH6soSEBACAq6urQV1XV1f9voIsXLgQ8+fPL5lAiYhKmxBAagwy7p3BtaiTuBh3HXfS7iM+NxkJyEK8XItUOZAqlyHzCctLG0328EHliiLmmCTnLRMJhVarxU8//YQxY8bAyspKX65UKgHoBmY+KjMzE5aWloUeb/bs2Zg5c6Z+O2+lNCKiciE3C9o7x/HvtR04Fn0cZ3LicM1ChijlI7+yFQ8fkOV9YcBGK6DSAkohgwV0/yqEHBaQw0LIoYQMCsghy/tPJnt0CzI83Jbpv3pYBsgKyjIKLJL917tQIFlhTy3ksPlr5sXy5PNULo2qd5DkvGUiodizZw/u3LmDCRMmGJR7enpCLpcjKirKoDwqKgpeXl6FHk+lUkGlUpVIrEREJUKjhrhxAFfD1+Cv2BPYbW2JGAsLXa5grdRXs1fL4ZxjDQeZExyUVeFkXR1V7WvC3akmqjjXQFUnD3jYO8PJxgpWyvyJBlFJKRMJxapVq9C6dWs0atTIoNzGxgbt2rXDtm3bMGbMGABAeno69u7di3nz5kkQKRGRmamzIcJ+RfCJJVilzEKYlRVgbwMAUGplcMp0gYMsALWrdUCHWs3Q0ssLNZysIZfz73EqWyRPKOLi4rB9+3Z8++23Be7/+OOP0aNHD8yePRtt27bFV199hapVq+LVV18t5UiJiMxICODiFoTv/x8WWKlx2d4SgBUUWhms0rxhr+yJ4Q16on+gD6rYs8eVyj7JE4pjx46hbdu2GD58eIH7O3fujAMHDuCbb77ByZMn0ahRI/zyyy+ws7Mr5UiJiMwk5R4S/pqBJfEnEeRoB8ASCq0cOYlt0dx1CGYMaIFATyepoyQyikwIUeHvq0lJSYGjoyOSk5M5bTcRSUedA82JFdgcugzL7a2QotCNcdAkNUOg/Yt4v3dLBFTn7ygqW4r6GSp5DwURUaXw7z5c/HsWPlJm4KKTLQDAIssNruqXMa/Pc+jg5yZxgETFw4SCiKgk5WYh+5//Yem/m/Crgz2ETAULjQUs0vpidsdXMKCJJwdYUoXAhIKIqKTcv4Sbf47FW4okXHN82FWc3AjP1ZqK98a0hrUlb+ukioMJBRGRuQkBceJ7bD62EJ852yFLbgkLtQq26WOxfMAoNPNyljpCIrNjQkFEZE6pMUgOmoz5qefwj6uuV0Ke7oOBPnPwdo9WnGyKKiwmFERE5nIxCKG7Z2G2kyVibG0gEzLYZ/bHin4z0bgmeyWoYmNCQURUXJlJyP5rFr6L+gc/ujpAK5PBIscRz9V4D3N79oKlhRkX7CIqo5hQEBEVx40DuLhtCv5np8V1J0cAgH1OK3zz7MdoWtND4uCISg8TCiIiU+RkIOef/2H19U34wdkRapkFFGor9PCYjoW9RsFCwV4JqlyYUBARGSvqNC78+QrmqTJw1dkJAOCgboxlPReihWfhKyETVWRMKIiIikqTi4wDC7HmwmqscrKHWmYJC40KPT1exyc9X2CvBFVqTCiIiIoi9goO/P4yvrRMwi1n3e2gLppGWNb7UwRWZ68EERMKIqInycnA7b/n4dvbf+AvOxsASqg0KgysPQuzOw2HXM5eCSKACQURUaFyL+/Eb3vexAp7IMXOBjIB1LPsiGUDP4aHg4vU4RGVKUwoiIgelxyFU39MwfLcywh3UgEAnNUumNXuU/Sv31bi4IjKJiYURER5NGpE7VuC9Ve+w3oHa6itVFBq5XjG/UUs6DkNlhZKqSMkKrOYUBARAUi9Goy/dk/DGtssRDvaAAB8ZAH4/PkvUK+Kp8TREZV9TCiIqFLLfXADe/58HevxL845qQBYwF5rg6nNP8Coxn2kDo+o3GBCQUSVkshIxPGgd7ExaR/221oDUMFSK0PvGsPwXteZsFHaSB0iUbnChIKIKhdNLs7v/hxBN9fgT3sV1LbWkAmgpW0bfPLsJ6hmV1XqCInKJSYURFQ5CIFrRzdg9+lPsNEeSHWwAgDUkflgXq/PEFitvsQBEpVvTCiIqGITAhGh27At5CME2WYgzkkBAHDXOmJGu/l4vm43iQMkqhiYUBBRhXUjdCf+CpmHrTapiHW0AKCAs1aFFxtOxdjmL0AhV0gdIlGFwYSCiCqcq6G7sfvYPGy3TcZ9BwsAFnDSKDGi7ni82uZVKBWcT4LI3JhQEFGFcTF0D/45Ng87bRNxz/FhIqFVYpjvy5jYbhIsFZZSh0hUYZWJhCImJgZHjhyBjY0NunbtCmtra/2+o0eP4vLlywb1XVxcMGjQoNIOk4jKIKHV4vjBjdh3cSn22qQh3lEBwAKOWgsMq/MCJrZ/DSqFSuowiSo8yROKRYsW4cMPP0SXLl1gY2OD9957D1u3boWXl2454F9++QV///03unX7b+BUjRo1mFAQVXI5OTnYu/MrHLmzFnvtNMh0kANQwEVjgcG+o/Bqu9dhZWEldZhElYakCcXvv/+O9957D/v27UPnzp0BADdv3oRGozGo17x5c6xatUqKEImojElJTcbubR/jWOJfOGQnh9pBBkCOGlobjA54BSNavgylnGMkiEqbpAnFokWLMHjwYH0yAQC1a9fOV+/BgwdYv349HB0d0aJFC7i7u5dmmERUBly/eR279n6Is5pTOGmjBOx1d2jUE654tfWb6F7vOchkMomjJKq8JEsoMjIyEBYWhkmTJuHSpUsICwtD9erV0b59e1haGg6cunHjBoKCghAVFYWwsDB88cUXmDJlSqHHzs7ORnZ2tn47JSWlxNpBRCVHrdZg/6E/EHz5G5y0ikW0ygKAEnIh0MLCB1O7zkXTGq2kDpOIAMiEEEKKE0dFRaFmzZro06cP/v33XzRv3hynT5+GVqvF33//re+pCAkJQYsWLaBU6rowV65cicmTJ+PkyZNo1qxZgceeN28e5s+fn688OTkZDg4OJdcoIjKLBwmJ2L7zU5xJ3oUQGy1y5LqeB3st0NWhBSZ1mwtPp1oSR0lUOaSkpMDR0fGpn6GSJRTx8fFwc3NDYGAgTpw4AUtLS6jVarRv3x7Vq1fHli1bCn1u1apVMXPmTLz77rsF7i+oh8LT05MJBVEZptUKHDt9DIdCP0O4/BquWP3XgVpLY40h/sMxrM1rHGhJVMqKmlBIdsnD1dUVLi4u6NGjh/4Sh4WFBXr27Il169Y98bnW1tZISEgodL9KpYJKxdvEiMqDO/fjsO2fxbiQ/DdO22Qj00YOwAJKIdDGwgfjO7yNZt4dOT6CqIyTdFBm//79ce7cOYOyc+fOoVYtXVemRqNBdHQ0PD099ftDQkJw584dtG3btlRjJSLzycpRY/fhLQi5+gPOWUbhrqUCsAMAOTw0CvSo2gXjur4HV9sqUodKREUkaULx4YcfonXr1hg+fDjat2+PEydOYO/evdi3bx8AQAiBZ599Fu3atUODBg1w+/Zt/PDDDxg2bBj69+8vZehEZCStVuDE+XDsP7kMV9WncdZaQGsnA6CAlVagrYUXRrd6Da38+7A3gqgckmwMRZ64uDisWbMGkZGR8PLywqhRo1CzZk39/uzsbPz2228ICwuDs7MzOnbsiK5duxp1jqJe/yEi8xJC4EJEBHYHL8HVtKM4Z5ONTLlcv7++xgbP+zyPIe3fgI3KTsJIiagwZX5QZmliQkFUuq5H38eug8txMXEvzlunIVXxXxJRTS1DR/vGeLHzO/BxbyRhlERUFGV+UCYRVRxCCFy8HY2/j36H64n7cdkqCQkWcv24CBcN0E7liyGtJqKZby9e0iCqgJhQEJFJhBAIvf4v/jn+Df5NC8EVq4c9EQ+TCAeNQCtFTQwMfAntGw6FQsFfN0QVGd/hRFRkmTkaHDx3EiHnfsTN7DO4YpWFbLkcsAUAOZw0Ai0U1dGn/iB0aTYWSgvevk1UWTChIKInuh2XjF3H1uFi9E7clt3GTZWAUMgAGwCQo5oaaGlZC30bj0Qr9kQQVVp85xORgVyNFscuhuNQ+BrcTAvFv1YpSFbIAeu8GjLUyZWhpW099G8+Fg38enNMBBExoSCq7IQQOHvnLvaHbsD1B4cQhTu4pdJCyGT6Sxm2Wi0aax3Q1r0t+rQaB/eqDaUOm4jKGCYURJXQrbhE/H1iAy5F7cVtEYFbqlyoZTJAv0yGDLVygUBLH3Sv3x9tA0dBqbSRMmQiKuOYUBBVcEIIXImJw76wINyI2osozQ1EqLKQJZcB+jGTMrirtWggnNCsagv0bPESPKoXvJovEVFBmFAQVTBqjRYnb1zHkfO/42b8cUSLKNxV5SJXJgMs82rJ4KzRooHGDk2cG6NHk6GoXbsbZAqFlKETUTnGhIKonItJTsf+i8G4fHMn7qadR5QiHvcsH06A+0gPhKtag7paGwTY10fXBv3RsH4/yC0sCzssEZFRmFAQlSOZOWocvBKKM1f/QmRSGO5r7yHKMhvZ8od3WVj/V9cnRwN/mTMaOjdGp4bPo3adHpAxgSCiEsKEgqiMSsnKwZHrZ3Hu+t+4mxCGGM0dRCkzkKZ4mDw8cvnCVquFb64cdZQeaOLeGl2aDoVLtUYAb+ckolLChIJIYkII3ElMQvDlA7h6OxhRaVfwQMTinjJbN3ASAJQPH5BBKQR8cwR8ZM6o51QPrX27on69PpBbO0vYCiKq7IqcULzyyitGHXjVqlVGB0NU0d1LSUXIv6G4EnkU0UkXEZ9zF3HyFMQqNdDK8vc8WGoFauVq4QkH1LGrjabeHdCyUX9YOtSQqglERAUq8vLlMpkMw4cPL9JBf/vtN5SlVdG5fDmVJo1W4Mr9+wi/dRq3ok/iXsplxOVG44E8FXEWjyQOj3HRaFArVwFPuSvqOPoh0KcdGtbtCQuH6rx0QUSSKepnqFEJRVGTBGPqlgYmFGRuyZnZuHw/CteizuN2bDjup/yLhOxoJCEZ8RbZSH/C3Zf2Gi281QLVhS1qWFVHHbcGaO7XDTVrtQeUVoU/kYhIAkX9DC3yJY/t27cX+eTG1CV6VF4iKiDyleUrf+Rrwy8FhBAQAtAILQQArRDQagW0D8vVGi1yNAK5Gi1yNZqH/wpkqXOQmJWM+6nxeJAShdTUKKRnRCEtJwap6jikilSkyrOQbKGF5vFOg8cW1qyqVsNDLYcH7B8mDgEIrN0BNb1bQ2bD8Q5EVLEUuYeiPGMPhelO3j2OTSdX41rKFSSINKihhQaAFgJamcFHup4o5GuDOhWgB18uBKpqNPBQy1AFtnBXuqKmfS34uQegfq1WsHNvCFjaSh0mEVGxmL2HAgCysrLw119/YfDgwQXu/+OPP/Dcc8/ByordtuXd7cRbeGf7BFwQMboC2cNHJWCl1cJRq4WjRgsHLWADSzjJbOBi4YSqth6o6eSFWu71UdOjMSxc6/AyBRERjEwoVq1ahX///bfQhOLw4cOIiYnBa6+9ZpbgSBohNw7irUPTkKwQsBACz6SpUVvmjaq2tWBlaQtLuSUUCiUUcgvIZHJdriEDIJNDJtPlHXKZ7OGS1rKHZbr/5HlfP3ySHP89Ty6T6XIWmUxfR/7I15DJ9HV0++QPxyrK9OfT7dN9LX/4NWT/xSCXATKZHAoZoJDLIZPJ/2u4hQoKG1eobKsC1s66B5MFIqIiMSqhWLt2Lb7//vtC97/00kuYNGkSE4py7Er0Wbxx+HWkKwC/7Fy87DgCfUbPgYWSMywSEVHhjEoorl27Bn9//0L3+/v749q1a8UOiqSRnp2K13e9hHQLICBLjfda/4DGgZ2kDouIiMoB+dOr/Eej0UCtVhe6PycnBxqNpthBkTTm/TEBMRYauKs1mN7kUyYTRERUZEYlFA0aNMDOnTsL3b979240aNCg2EFR6Tt5/SD25FwAAAyx7Ix2bfpLHBEREZUnRiUU48aNwxtvvIH9+/fn27d//3688cYbRk/RTWXDV4ffh1YmQ9t0OcaPWC51OEREVM4YNYbi1VdfxcGDB9GtWzcEBgaibt26EELg2rVrCA8Px+jRozF+/Hijgzh27Bj2798PGxsbDB8+HDVqGK5T8ODBA2zcuBH3799Ho0aNMGTIECgUT5iKkIxy4toBhFskQyYEBtadAaVSKXVIRERUzhjVQyGTybB+/Xps3LgRXl5eOHv2LM6fPw8vLy/89ttvWLdu3cNbBYtGq9Vi7Nix6NevHxITE5GYmIg+ffrgypUr+jo3btxAo0aNEBQUBI1Gg9mzZ6NPnz4cq2FGPwR/CABom6FE764vSRwNERGVR5LOlPn111/jnXfeQXh4OPz8/ADoZuTKzMyEu7s7AGDIkCGIiYnB4cOHIZfLcevWLfj7+2PNmjUYPXp0kc7DmTILF5N4B89ufRZqmQzzq7yCQX2mSx0SERGVIUX9DDWqh8Lcvv76a4waNUqfTACAg4ODPplQq9X466+/MHr0aMjlulB9fHzQuXNnBAUFSRFyhbNm7wKoZTL4ZWvRv9dUqcMhIqJyyqiEYu/evfD19YWvry/27t1brBOnpqbi6tWr6NixI4KCgvC///0PK1aswP379/V1IiMjkZWVhTp16hg8t06dOk+c7yI7OxspKSkGD8pPCIHDyccAAC0sGnJcChERmcyohGLy5Mn4/fff8fvvv2PSpEnFOnFycjIAYPHixVixYgWUSiV27NgBPz8/nDx5EgCQkZEBAPm6WBwdHZGenl7osRcuXAhHR0f9w9PTs1ixVlQnru7HXaUWllqBwZ3elTocIiIqx4y6yyMuLg7e3t4QQiAuLq5YJ7a3twcAuLi44O+//9aXP/vss5g9ezb27dsHOzs7AEBSUpLBc5OSkvTPL8js2bMxc+ZM/XZKSgqTigLsOL0aANA4W4m6vk0ljoaIiMozo3oopk6din79+qFfv36YMmVKsU7s6OiI6tWro2XLlgblLVu2xM2bNwEAXl5esLW1xdWrVw3qXLlyBfXr1y/02CqVCg4ODgYPyu9sxiUAQIBNoLSBEBFRuWdUD8VHH32EkydPQqvVok2bNsU++YgRI3Do0CFotVrI5XIIIXDw4EE0btwYAKBQKDBo0CCsXbsWkyZNgkqlwoULF3DkyBFs3ry52OevzK7dPYdblhrIhUDvVhOkDoeIiMo5SW8bTU5ORrdu3aDRaNC2bVucPHkSDx48wP79+/UDMe/du4eOHTvCzs4OzZo1w44dO9CzZ0+sW7euyOfhbaP5Lfp9CtZlBKNBFrBx4nmpwyEiojLK7LeNfvHFF0U+eVHrOjo6IiQkBB9++CECAgIwd+5cXL161eCuDg8PD5w7dw5z5sxB06ZNsWnTJqOSCSrYpaQwAICfopbEkRARUUVQ5B4KmUyGonZmGFO3NLCHwpBanYuOPzdFmkKG+dUnYlAPzj9BREQFK+pnqFFjKAIDA4sbF5UBh8/vRJpCBlutFt3bjJI6HCIiqgCKnFB8/vnnRT7omDFjTAqGSkfw5T8BAPWyLeFg7yJxNEREVBEUOaF48803SzIOKkVX0y8DlkAdld/TKxMRERWBpGt5UOnTCi0iFLoZSJt6dZc4GiIiqiiYUFQyp68FI00hg5VWi04tB0sdDhERVRBMKCqZY5e2AwBq58jh4OAqcTRERFRRMKGoZK4lnAUAeMqqShwJERFVJMVKKOLj43Hq1ClzxUKlIFIbCwDwdW4scSRERFSRmJRQJCUlYcCAAXBzc0OrVq305QMGDEBISIjZgiPzSslIwh2lBgDQul5fiaMhIqKKxKSE4u2330Z2dna+VUCnTJmCjz76yCyBkfkdPb8LWpkMrmoNAht0ljocIiKqQIyaKTPPjh07cPz4cXh5eRmUt2zZEocOHTJLYGR+5yOPAgA81ZaQyRUSR0NERBWJST0UiYmJcHZ2BqBbtyNPWloa5HKO8yyrIlN0PUrV5ByQSURE5mXSp3/z5s2xZcsWAIYJxaJFi9CuXTvzREZmFy3iAADeDnUljoSIiCoaky55LFy4EH369MHhw4chhMD//vc//P333zh//jwveZRRanUuoixyAcjQ0Luj1OEQEVEFY1IPRceOHREcHIzMzEz4+flh48aNqFWrFkJCQgzu+qCyI+zmCWTKdTNktmrUS+pwiIiogjGph+K7777DpEmT8Ouvv5o7HiohZ67+AwDwypXBxtZR4miIiKiiMamHYtq0adBoNOaOhUrQrYSLAIBqwkHiSIiIqCIyKaEICAjA6dOnzR0LlaCY7HsAAHfLGhJHQkREFZFJlzzGjRuHESNGYO7cuQgICIClpaXB/sDAQHPERmZ0X5YKAPByqSdxJEREVBGZlFBMnz4dADB27NgC9wshTI+IzC4nNwcxSi0AGRr6tJc6HCIiqoBMSigePHhg7jioBJ29eQK5Mt0dHo3r8pZRIiIyP5MSCjc3N3PHQSXo3M1gAECNXBksrWwljoaIiCoikxKKvXv3PnF/9+7dTQqGSkZknO4OD3fBZIKIiEqGSQlFjx49nrifYyjKlpisu4Al4G7hLnUoRERUQZl026gQwuCh0Whw5coVdO7cmZNdlUEPkAwAqOHoK3EkRERUUZnUQ/E4uVyOunXrYtWqVRg4cCBGjRpVpOd9++232Llzp0GZt7c3vvnmG6Pq0JPFKtQAZKhdrYnUoRARUQVlloQiT5UqVRAREVHk+ufOnUNSUhLeffddfZmDg4PRdahw9xKjkaLQrQjbyI93eBARUckwKaG4e/duvrLExER8/vnnCAgIMOpY1apVw/PPP1/sOlSw8OtHAQBuag2qVa0lcTRERFRRmZRQeHp6Flju7++PDRs2GHWss2fPYtiwYXB0dETHjh0xZswYyOVyo+tQwW5GhwMAqqotAJlM2mCIiKjCMimhuHz5cr4yZ2dnuLsbdxeBUqnEM888g86dOyMqKgpvv/02Nm7ciL/++guyhx9+RanzuOzsbGRnZ+u3U1JSjIqrIolK+RcA4Ao7iSMhIqKKzKSEYsSIEQgPDy9wX2BgYKH7HvfJJ58YjIfo3bs3AgMDsWXLFgwaNKjIdR63cOFCzJ8/v2iNqeAeZMcAlkAVZVWpQyEiogrMpOsGZ8+eLbBco9Hg4sWLRT7O44MrGzVqBG9vb4OVTItS53GzZ89GcnKy/nHnzp0ix1TRxAtd70x1e46fICKikmNUD8Xx48cL/BoAtFotQkJC4OXlZXIwQggkJSVBqVQWq45KpYJKpTI5jook1iIXgAy13BtJHQoREVVgRiUUbdu2LfDrPNWqVcOyZcuKdKzc3Fz8/PPPGDdunH4sxIIFC5CcnIz+/fsXuQ4VLi4lHsl5t4zWaSdxNEREVJEZlVCkpqYC0C0OFhcXZ7BPqVQa1SugUCgQGhqKDz74AH5+frhz5w5SU1Pxyy+/oGnTpkWuQ4U7f/MEAMBRo4VH9boSR0NERBWZTEi88EZ8fDwuXLgAZ2dn+Pv7w8rKyqQ6T5KSkgJHR0ckJydXqkmxfty1CEti18E3W2DLqxekDoeIiMqhon6GmnWmTFO4urqic+fOxa5D+UUn6m4ZdRHGJWBERETGMukuj6ysLHzwwQcICAiAnZ0drKysDB5UNsRlRgMAnOVO0gZCREQVnkkJxdy5c/Hnn39izpw5SE9Px7p16zBt2jTI5XK89dZb5o6RTJSgSQAAuKm4bDkREZUsky55/Pbbb9i+fTsaN26MF154AYMHD8aQIUPQvHlzfP311+aOkUyUKMsAAFRz9JE2ECIiqvBM6qG4c+cOGjRoAACwtbVFcnIyAOC5555DaGio+aKjYolXaAAAtdwbShwJERFVdCYlFEIIKBQKAICvry/++ecfAMDp06cr1V0UZVlMcixSH85BUd+nhcTREBFRRWdSQuHo6Kj/+o033sCLL76IFi1a4Nlnn8Wrr75qtuDIdBdungQAOGk0qOruK3E0RERU0Zk0hiIpKUn/9UsvvYQ6derg+PHjqFevHp5//nlzxUbFEBF9HgDgppZz2XIiIipxJiUU3333HSZNmqTf7tChAzp06GC2oKj47qfcBMA5KIiIqHSYdMlj2rRp0Gg05o6FzCghMwYA4CizlzgSIiKqDExKKAICAp64fDhJL0mTCABwtnSTOBIiIqoMTLrkMW7cOIwYMQJz585FQEAALC0tDfYHBgaaIzYqhiSkAwDcbKpLHAkREVUGJiUU06dPBwCMHTu2wP0SrzdGABLluQBkqO5SR+pQiIioEjApoXjw4IG54yAzUmvUSNRNEwLvapzUioiISp5JCYWbG6/Ll2UR969DI5NBIQT8fJpKHQ4REVUCJg3KBIC4uDisXr0a7733nr7s+PHjvPujDLh+NwwA4KbRwtbOWeJoiIioMjApoQgPD0f9+vWxePFiLFiwQF/+008/Ye3atWYLjkxzJ/YaAMBFrZA4EiIiqixMSihmzZqFmTNn4tKlSwblkyZNwpIlS8wSGJkuNjUSAODISa2IiKiUmDSGIjQ0FEFBQQAA2SPTOvv5+eHatWtmCYxMF595H1AAjnIu1EZERKXDpB4KhUKB1NTUfOVXrlyBi4tLsYOi4knS6JaTd1Zy8CwREZUOkxKKPn36YN68edBqtfoeisjISEyePBl9+/Y1a4BkvBRkAACcbTmpFRERlQ6TEorFixcjJCQEHh4e0Gq1aNSoEXx9fZGVlYWFCxeaO0YyUrJcDQBwd/KRNhAiIqo0TBpD4e7ujtOnT2Pr1q0IDQ2FVqvFnDlzMHjw4HzTcFPpEkIgUSEAyOBd1V/qcIiIqJIwKaEAAEtLSwwdOhRDhw41ZzxUTNGJ95D78DJU7ZqNJI6GiIgqC5MTips3b+Kbb77B5cuXAehWIJ06dSp8fHzMFRuZ4Pqd8wAAR40GLq41JI6GiIgqC5PGUGzfvh316tXDoUOH4O3tDW9vbxw8eBB169bFzp07zR0jGeFO7FUAgItGBjxySy8REVFJMqmHYtasWZg/fz5mz55tUL5w4UK88cYb6NOnT5GOM3fuXKxfv96gzN/fP19SsnXrVixfvhz3799Ho0aN8NFHH8HX19eU0Cu8+8kPJ7XSKiWOhIiIKhOTeiiioqIwZcqUfOWTJ09GVFRUkY8TGxsLPz8/7N69W//4/vvvDers2LEDQ4YMQd++fbF69WoIIdCxY0ckJCSYEnqFl5BxDwDgAGuJIyEiosrEpIQiICAAZ8+ezVd+9uxZBAQEGHUsW1tb+Pr66h+enp4G++fNm4cxY8ZgxowZaN26NX7++Wfk5ORgxYoVpoRe4SXlxAMAHBWcJZOIiEqPSQnF6NGjMXToUHz55ZcIDg7G4cOH8eWXX2LYsGEYM2YMwsPD9Y+nOXLkCBo3boyOHTtizpw5BjNwpqam4syZM+jVq5e+zNLSEt27d8fBgwdNCb3CS9XqXj9nFWfJJCKi0mPSGIo33ngDgG4sxeOmT59usC2EKPQ4Li4ueOedd9C5c2dERUXh/fffx44dOxAaGgpLS0vcvXsXQghUq1bN4HnVqlXDuXPnCj1udnY2srOz9dspKSlFaldFkIpMAEAVOw+JIyEiosrEpITiwYMHZjn5Rx99BIVCt8R2s2bN0LRpU9SqVQsbN27Eiy++CK1WCwBQKg0HGFpaWkKj0RR63IULF2L+/PlmibG8SVboZsms7uIjbSBERFSpmJRQuLmZpzs9L5nIU7NmTfj4+ODixYsG54mPjzeoFxcX98QYZs+ejZkzZ+q3U1JS8o3NqIiyczVIypsl072e1OEQEVElYvLEVgCQnp6OxMTEfOU1a9Y06XjZ2dmIiYmBk5MTAN0U356enggJCUG/fv309Y4ePYrevXsXehyVSgWVSmVSDOXZzQd3oH4494RP9QYSR0NERJWJSYMyw8PD0bRpU9jZ2cHT0zPfoyhycnLw1ltv6XsfMjIyMGXKFGg0GgwfPlxfb9KkSVi1ahWuXLkCAPjhhx9w8+ZNTJgwwZTQK7Rb0bqeHSeNBir7ak+pTUREZD4m9VCMGzcODRs2xIoVK/S9CcZSKpWoWbMmGjduDK1Wi8TERDRp0gT79u1D7dq19fXeeecd3L17F02aNIGtrS0UCgXWrVuHRo24TsXj7sXfAMBZMomIqPTJxJNuwyiElZUVHjx4AHt7e7MEcf/+fTg6OsLKyqrQOunp6UhISICHhwcsLIzLg1JSUuDo6Ijk5GQ4OFTc+Rk++30WfsnYg6ZZCvw8MVzqcIiIqAIo6meoSZc8/Pz8cPv2bZODe5y7u/sTkwlANwGWp6en0clEZZKYyVkyiYhIGiZ9Oi9atAgvv/wyPvjgA9SpUweyx7rX69XjHQZSSMmJB5ScJZOIiEqfSQmFRqPB5cuX0bdv3wL3m3AVhcwg5eEsmU6WrhJHQkRElY1JCcWMGTMwevRovP766yYPyiTzS5XpZsl0seUsmUREVLpMSihiYmKwePFi2NnZmTseKoYUuW6WTHcnb4kjISKiysakQZn16tXDv//+a+5YqBg0Wg0SFbpLTTWr1JU4GiIiqmxM6qEYPHgwRo4ciQ8//BC+vr75BmUGBgaaIzYyQnRyzH+zZNbgLJlERFS6TEoo3nvvPQDAsGHDCtzPQZml7+Zd3SyZzhoNnFxqSBwNERFVNpKuNkrmc/fBNQCAM2fJJCIiCUi62iiZz4PkOwAABy0n/iIiotJn0qBMQLeE+OrVq/WXPwDg+PHj0Gg0ZgmMjJOYEQMAsBOVb5VVIiKSnsmrjdavXx+LFy/GggUL9OU//fQT1q5da7bgqOhScnTLyNvJeSsvERGVPpMSilmzZmHmzJm4dOmSQfmkSZOwZMkSswRGxklVpwAA7C0cJY6EiIgqI5MuuIeGhiIoKAgADG4Z9fPzw7Vr18wSGBknXWQAABxVnHabiIhKn0k9FAqFAqmpqfnKr1y5AhcXl2IHRcZLleUAAFzsqkkcCRERVUYmJRR9+vTBvHnzoNVq9T0UkZGRmDx5cqELhlHJSns47XZVx5oSR0JERJWRSQnF4sWLERISAg8PD2i1WjRq1Ai+vr7IysrCwoULzR0jPYVGo0WyQvd1jSp1pA2GiIgqJZPGULi7u+P06dPYunUrQkNDodVqMWfOHAwePBiWlpbmjpGe4k7iA/20297VuI4HERGVPpMSisDAQISHh2Po0KEYOnRogfuo9ETe1w2EtdNqYe9UXeJoiIioMjLpksfZs2cLLNdoNLh48WKxAiLjRcfdBAA4aQSg4EyZRERU+oz69Dl+/HiBXwOAVqtFSEgIvLy8zBMZFdmDpEgAgINWIXEkRERUWRmVULRt27bAr/NUq1YNy5YtK35UZJTE9Lxpt5USR0JERJWVUQlF3twTbm5uiIuLM9inVCqhUnEdCSmkZMUBCsBeZi11KEREVEkZlVDY2enWicjKyiqRYMg0aZokQAE4WthLHQoREVVSJq82SmVHhjYdAOBgyVlKiYhIGmXqlgC1Wg25XA65/L88R6vVQqvVGtSTyWRQKDgAMU+6TNdj5GJdReJIiIiosiozPRSzZs2CUqnEzJkzDcqnTJkCS0tLWFlZ6R/NmzeXKMqyKU2eCwBwteccFEREJA2TEorAwECT9hVm165d2L17N/z9/QvcP2jQIKjVav2DE2f9RwiBVLmuB6eaM2/ZJSIiaUg+sdW9e/cwYcIErFu3DtbWvEvBWIkZmUhT6Kbd9nT3lTgaIiKqrCSd2Eqr1WLMmDGYPn06mjZtWmi9Xbt2QaVSwcHBAR07dsTnn3+OOnW4CBYARMbdBgAohIB7FSYUREQkDUkntlqwYAG0Wi1mzZpVaB1/f39s3rwZnTt3RlRUFKZPn47OnTvjwoULcHJyKvA52dnZyM7O1m+npKQUOaby5m7sdQCAk0YLuZWTtMEQEVGlJdnEVidPnsTSpUv1q5Xm3ckhhIBarYaFhS60Rwdp+vn5YcOGDahWrRp+++03TJw4scBjL1y4EPPnzzemaeXWg0TdtNuOWhnwcMVRIiKi0ibZxFahoaFISkqCr+9/3fQajQbnzp3DN998g+zs7AJvDXV0dISnpyf+/fffQo89e/Zsg0QkJSUFnp6exY65LIpPvQuA63gQEZG0TL5tNC4uDqtXr8Z7772nLzt+/Dg0Gk2Rnj9lyhSDOzfUajWaNGmCadOmQa1WFzrPRGJiIm7fvo0aNWoUeuy88RaPPiqq5MxYAIAdOO05ERFJx6SEIjw8HPXr18fixYuxYMECfflPP/2EtWvXmi247Oxs9O3bF8eOHUNycjLOnz+PYcOGwdHREaNHjzbbecqz1JwEAICdzFbiSIiIqDIzKaGYNWsWZs6ciUuXLhmUT5o0CUuWLDE5GAsLC4OeCZVKhalTp+L999+Hj48P+vfvj+rVq+PUqVOoUoWzQgJAukY34NTeouL2whARUdln0tTboaGhCAoKAqCbBjuPn58frl27ZnIwoaGh+cp69eqFXr16mXzMii5d6NbxcFJxHQ8iIpKOST0UCoVCf8fHo65cuQIXF36wlaZ0me72WCcbd4kjISKiysykhKJPnz6YN28etFqtvociMjISkydPRt++fc0aID1Zqlw3CLaKY02JIyEiosrMpIRi8eLFCAkJgYeHB7RaLRo1agRfX19kZWVh4cKF5o6RniBFLgAA7i4+0gZCRESVmkljKNzd3XH69Gls3bpVPzHVnDlzMHjwYFhaWpo7RipEUkYacuS6HqLqVWpLHA0REVVmJiUUAGBpaYmhQ4di6NCh5oyHjHA37g4AwEIIVHXlSqNERCSdIicUR44cKfJBO3ToYFIwZJx78bcAAA5aLeQqO2mDISKiSq3ICUXHjh2LfFAhhEnBkHHikqMAAHYacB0PIiKSVJETiszMzJKMg0yQmBoDALDjOh5ERCSxIt/lYWVlpX+0adPGYPvxfVQ68tbxsDF9KAwREZFZmHTb6NmzZwss12g0uHjxYrECoqJLy04EANhwYTAiIpKYUX/aHj9+vMCvAUCr1SIkJAReXrzboLSk5yYDCsBGbiN1KEREVMkZlVC0bdu2wK/zVKtWDcuWLSt+VFQkGZo0QAHYKbjSKBERScuohCJv/Q43NzfExcUZ7FMqlVCp2PVemjKEbqCsndJJ2kCIiKjSMyqhsLPTzXWQlZVVIsGQcTKRAwBwsOKCbEREJC2Tbg/Yu3fvE/d3797dpGDIOBmyXACAk21ViSMhIqLKzqSEokePHk/cz4mtSke6XAtABhf7alKHQkRElZxJt40KIQweGo0GV65cQefOnfHrr7+aO0YqgFarRerD+ayqOHtKGwwREVV6JiUU+Q4il6Nu3bpYtWoVly8vJQkZqVA/nG67hpuPtMEQEVGlZ5aEIk+VKlUQERFhzkNSIe7GRQIALLUCzs41JI6GiIgqO5PGUNy9ezdfWWJiIj7//HMEBAQUOyh6uph4XUKhW2mU81AQEZG0TEooPD0Lvmbv7++PDRs2FCsgKpr4ZF1SZ6flKqNERCQ9kxKKy5cv5ytzdnaGu7t7sQOioklJfwAAsBFmvWpFRERkEpMSinr16pk7DjJSWpZuYTArwZVGiYhIekb/eSuEwIYNG/D888+jbt268Pf3R9++fbFhwwbOP1GK0nOSAQBWUEocCRERkZEJhUajwdChQzFq1CjExcWhTZs2aNu2LR48eIBRo0Zh2LBh0Gg0JRUrPSIjNwUAYCXj+ilERCQ9o/rLly1bhhMnTiAkJARt2rQx2BcSEoJhw4Zh+fLleOONN4wOJCMjAzdv3oSbmxuqVcs/82NSUhLi4uLg5eUFS0tLo49f0WRp0gELwFpuLXUoRERExvVQrFmzBt9//32+ZALQLWf+3XffYc2aNSYFMnbsWDRq1AiffvqpQblarcb48ePh7u6OTp06oWrVqvj5559NOkdFkqXVLdBmzaXLiYioDDAqobh+/TqeeeaZQvc/88wzuHbtmtFB/PDDD4iJiUGjRo3y7VuwYAF27NiBS5cuITo6GsuXL8fYsWMRHh5u9HkqkiyhSyhsLOwkjoSIiMjISx4qlQqpqamwsrIqcP+T9hXm4sWLmDdvHo4fP45+/frl2//DDz9gwoQJqFOnDgDgxRdfxMKFC7Fq1Sp8/fXXRp2rIsl6uHS5rcpR4kiIyjYhBNRqNcd3ERVCoVDAwsICMlnx5jUyKqFo1aoVVq5ciTlz5hS4f+XKlWjVqlWRj5eZmYnhw4fjiy++gJeXV779MTExiIqKQuvWrQ3K27Zti9OnTxsTeoWTJdP9crSzcpI2EKIyLCcnB/fu3UNGRobUoRCVaTY2NvDw8CjWGEWjEorZs2ejV69eiI2NxYwZM+Dt7Q0AiIyMxNKlS/HNN99gz549RT7ejBkzEBgYiFGjRhW4Pz4+HgDg6upqUO7q6oq4uLhCj5udnY3s7Gz9dkpKSpFjKi8yHyYUjjZuEkdCVDZptVpERERAoVCgevXqsLS0LPZfYEQVjRACOTk5ePDgASIiIuDn5we53LQJE41KKJ555hmsWrUKr732GpYtW6a/vJGVlQVbW1usXr0aXbt2LdKx/v77b2zatAk7duzAhQsX9MeJj4/HhQsX0LBhQ1hY6MLLyckxeG52djaUysLnX1i4cCHmz59vTNPKnUy5ACCDk11VqUMhKpNycnKg1Wrh6ekJGxsbqcMhKrOsra2hVCoRGRmJnJwco4cu5DF6msWXXnoJzz//PLZt26YfgOnv749+/frl60l4koSEBNSoUQMTJ07Ul0VGRiIuLg5hYWE4e/YsatSoAZlMhnv37hk89969e4WuJwLoelJmzpyp305JSXli/fJGq9Ui/WEC6eqY/xZbIvqPqX9tEVUm5nifmDRvs6urK8aOHVusE48cORIjR440KAsMDESXLl2wdOlSAICdnR1atmyJnTt36utmZ2dj7969eOeddwo9tkqlgkpVcSd8SsnKgPph120VJy5dTkRE0ivzqfv8+fOxYcMGLFq0CIcPH8aIESNga2uLSZMmSR2aZO4nRgMAZELA1bm6xNEQUUUxZ86cYv+xWBSffPJJoWPnqPwqUwmFr68vPDw8DMp69+6NHTt24NChQ3jjjTfg6OiII0eOwMnJSZogy4DYpLylywUsVJyHgojMIyUlBYmJiSV+ntTUVCQkJJT4eah0lamlKjdv3lxgee/evdG7d+9SjqbsSkq+DwCw1QLgqHUiIioDylQPBRVNUnosAMBaMJkgqmh2796NOnXqYNu2bWjTpg1q1KiBq1evAgB27dqFbt26wcfHBx07dsy3DMHly5cxZMgQ+Pn5oU2bNli+fDm0Wq1+/48//ohWrVrB29sb3bt3x6FDhwqN4++//4aPjw/S09MNyj/99FODGZOfFhMALFmyBA0aNEDjxo0xffr0fMekiqFM9VBQ0aRk6LoKrbXMB4mMIYRAZm7pz5hprVQUeQ6MrKws3Lx5Ex999BGWLVuG2rVrw83NDX/++ScmTpyI5cuXo3Xr1rh8+TImTJiAnJwcvPLKKwCAvn37okePHvjrr7+QlpaGX3/9Fbt27cJzzz2HtWvXYtq0afj222/RunVrrF+/Hj169EBYWBgaNGiQL44uXbogOTkZQUFBGD16tL589erV+nEWRYlp5cqVmDdvHn744Qc0adIE3333Hb755hv07NmzuC8rlTEyIYSQOoiSlpKSAkdHRyQnJ8PBwUHqcIrts99exy9ZB9EsS4m1E89IHQ5RmZSVlYWIiAjUqlVLf199Ro4aAXP/LvVYLn3YCzaWRfv7LSgoCAMHDsS5c+cM1jeqW7cuZsyYgcmTJ+vLVqxYgW+//Rbnz59HVlYWrK2tcfLkSbRs2VJfRwgBmUyGOnXqYPTo0fjwww/1+7p16wYPDw+sW7cOADB16lTcvXsXQUFBAIBXXnkF0dHR2LlzJwDgxIkTaNu2LSIiIuDt7f3UmACgdu3aGD9+PN577z19PA0bNoSnpyd2795tzMtIJaig90ueon6GsoeiHMrI0c38aQUu405UESmVSjRs2FC/nZCQgGvXruGDDz7AJ598AiEEhBDIyspCZmYmAMDKygr9+/fHsGHDMG7cOHTt2hVt27aFQqFAeno6bt68iY4dOxqcp3PnzoWOXQOA0aNHo2fPnoiNjUXVqlXx66+/okOHDvD29i5STGlpaYiIiED79u31x5TJZGjfvj1u375tzpeMygAmFOVQhjoVkAPWsoo71wZRSbBWKnDpw16SnNcYj08TnreUwNdff41OnToZ1H203p9//okdO3bg77//xssvvwy1Wo2goCD9MgmPr9NgaWmZbybiR3Xp0gXVqlXDpk2bMHnyZPz222/46KOPihxTbm4uAOSb2bg460VQ2cWEohzK1KTrEgqFtdShEJUrMpmsyJceypJq1arBzc0N586dw7BhwwqtJ5fL0a9fP/Tr1w9arRbPPvssPvvsM6xfvx5VqlRBWFgYOnfurK9/5swZ+Pv7F3o8mUyGkSNH4tdff4Wvry+SkpIwdOjQIsfk7OwMNzc3nD171qCXIjw8HHZ2vOW9ouGovnIoS+i6E60VthJHQkSlQSaT4e2338aSJUuwZcsWaDQa5OTk4J9//tH3GERFRWHGjBn6SwnJycl48OABqlXTTc8/bdo0fPHFF/q1k7Zu3YotW7Zg+vTpTzz3mDFjcPz4cXz88cd47rnn4OzsXOSYAGDy5Mn47LPPcPXqVQghsGbNGhw9etTsrxFJr/yl6oQsoeuitFHaSxwJEZWWN998ExYWFpg6dSpGjRoFCwsLdOjQAQsWLAAAuLu7o1atWujUqRMSExORm5uLvn37Yt68eQCAd955B7GxsWjVqpWup8bGBkuXLkW3bt2eeN7GjRujUaNGOHr0KP744w+jYgJ0s2/evHkTDRo0gK2tLZo0aYLhw4cjKSnJrK8PSY93eZRDQ39ohiuqXEyx7YbJQ5ZKHQ5RmfSkUetlWXZ2NpKTk1G1auErCSclJcHe3h4KRcFjM1JTU2Fra1vggk9qtRopKSlwdnbOdytramoqNBpNvpmIk5OTkZmZiSpVqhR6zqfFlJGRAQsLC1haWiItLQ25ubn63g6SHu/yqKSyZLr76O2t+GYkqmhUKtUTkwkAT116wN6+8N5LCwsLuLi4GPU8R0dHODo6FiumR5eQ5/iJioljKMqhzIcJhaONm8SREBER6TChKIcyHn7XnO2f/FcMERFRaWFCUc5otVp9QuHiyKXLiYiobGBCUc4kZqRA83AgVRVnJhRERFQ2MKEoZ2Li7wIAFELA1clD4miIiIh0mFCUM3FJUQAAW62AXFl+boUjIqKKjQlFOZOYch8AYKuVOBAiIqJHMKEoZ5IzHgAArIXsKTWJiIhKDxOKciY1Ix4AYC2MW72QiMqHS5cu4eWXX5Y6jGK7ceMGXn75Zf2Ko1TxMaEoZ9KykwAAVoKTnBJVRNHR0Vi7dq3UYRTbgwcPsHbtWmg0GqlDoVLChKKcychNAQBYySwljoSIiOg//DO3nMlUpwJywBq8w4OoMjh8+DC2bduGL774Ql92/fp1fPLJJ1i5ciWUSiXOnDmDH374Ae+88w5+//13REREoGnTphg/fjwiIyOxdu1axMfHo2vXrhg8eLD+OHnPe//997Fu3TpERkaiRYsWePnllw0W+VKr1di0aROOHj0KOzs79OjRA927dzeIMyEhAStWrMC9e/fQsGFD1K9fv+RfHCpT2ENRzmRqMgAAVgqbp9Qkoorg2rVr2Lx5s0HZ/fv3DS4n3L59G6tXr0aPHj2Qk5MDHx8fvPPOOxg4cCB69+4NlUoFDw8PvPzyy1i9erX+OHnPa9euHVJTU+Hj44P3338fr7zyir5Obm4uevXqheXLl8PPzw9OTk546aWX8OGHH+rrZGRkoE2bNtixYwf8/Pxw5MgRjBgxooRfGSpr2ENRzmRpswAANha2EkdCVA4JAeRmlP55lTaArGTvzFKr1VizZg06duwIAMjJycHcuXNx8uRJtGzZEoBuGeq1a9di/PjxBs/78MMP9QNBu3XrhlatWmH69OkIDAzEd999h5iYGISFhcHSUneptXv37mjXrh0mT56MKlWqYMWKFcjIyMC5c+f0S1+PGDECv/32W4m2mcoWJhTlTDayAQA2ysLXpCeiQuRmAAskmLJ+TjRgWbJ/BCiVSrRv316/XadOHdja2uqTibyyx3s7ABhcBmnRogV8fHwQHByMwMBA/PXXX9BqtZgyZQqEEBBCQKvVQq1W4+LFi+jSpQsOHTqE5557Tp9MAMCwYcOYUFQykicUWVlZCAsLQ1paGurXr4+aNWsa7A8LC8ONGzcMyhwdHdGjR4/SDLPMyILuFiw7laPEkRBRWWJpaQm5/L+r2HK53OADPq/s8bsulEol7O3tDcpcXFzw4IFuzpv4+Hh4eXmhQ4cOBnW6dOkCX19fALo7Oho0aJDvGFS5SJpQ/PTTT5g/fz68vLxgYWGBkJAQjBkzBt99953+jfH9999jx44daNOmjf55Xl5elTahyJTpfhnYW/PNSmQ0pY2ut0CK85pIpVIhJyfHoCwxMbG4Eenl5ubi3r178PDQrQ0khMCdO3fg6ekJAKhZsyYyMjKeODeGl5cXbt++bVAWGRlpthipfJA0obC1tUV4eDgcHXV/bYeHh6Np06bo168fnn/+eX29Nm3aFNhNVxllyrUAZHCwrSJ1KETlj0xW4pcezM3Pzw8xMTH4999/4evrC61WizVr1pj1HEuWLMFnn30GAPjll1+QkpKC5557DgDw4osvYvDgwdi1axeeffZZALqk47ffftMPvBw0aBDGjx+PiIgI1KpVC1lZWfj222/NGiOVfZImFEOHDjXY9vf3h4WFBZKTkw3Kk5OTsXPnTjg6OqJRo0ZwcKi84wfSH/Zoutm7SxsIEZWKNm3aoG/fvmjfvj06d+6MixcvokaNGmY7vkqlwtGjR9GhQwfY29tj//79+Pzzz1G9um6sycCBAzFv3jwMHDgQrVu3hqurK86dO4dWrVrpE4q88RLNmzdH586dceHCBXh5eZktRiofJB9Dcf/+fQQHByMlJQU///wzevbsiSFDhhjUOX36NJYuXYqoqChER0dj+fLleOGFFwo9ZnZ2NrKzs/XbKSkpJRZ/aUrJTEe6XDdSvKa7r8TREFFJaNCgQb4eiKCgIBw5cgSxsbFo0qQJHBwcsHPnTiiVSgBAs2bN8P333xs8p3Xr1vj6668Nyjp37gwnJyeDMgsLCxw6dAinTp1CZGQkli1bBn9/f4M6c+fOxauvvoqQkBDIZDI0bdoU3t7e+v0ymQx//vknjhw5gpiYGDRo0ABubm7466+/9DFSJSAkFhYWJgYPHiw6d+4sqlSpIhYuXCi0Wq1+/759+0R6erp++/PPPxdKpVKcP3++0GN+8MEHAkC+R3Jycom2paSdvRkmGv7UUASuaSA0OZlSh0NUpmVmZopLly6JzEy+VwqzZcsWYWtrK3UYVAY86f2SnJxcpM9QySe2CgwMxObNm3Hw4EHs378fH330EX788Uf9/meeeQY2Nv8NaJo1axYcHBywa9euQo85e/ZsJCcn6x937twp0TaUluhY3d0ujhoBuZIzZRIRUdkheULxqIYNG6JFixY4dOhQoXVkMhns7e0RGxtbaB2VSgUHBweDR0VwP0k3atpBy6XLiaj4CrpUQmQqyRIKjUaDhIQEg7L09HRcv35dPxeFVqtFXFycQZ2wsDD9fPOVTWKq7nY3O63kQ1+IqALw8vLC6NGjpQ6DKgjJPplyc3PRvn179O/fHwEBAUhKSsJPP/0EKysrTJs2DYAu6ejQoQP69euHBg0a4Pbt21i2bBl69eqVb+BmZZCcqZtoxg4qiSMhIiIyJFkPhZWVFU6ePInq1avjwIEDuHbtGl577TVcunQJ1apVA6Cbwe3UqVOoWbMmDh48iMTERKxcuRK7du0yWAmvskjN1U1mYysvX/fRExFRxSdp37m9vb2+N6I4dSqLNHUKoADsLTjtNhERlS1lalAmPVk6dKskOqo47TYREZUtTCjKkXTo5vN3tOEsmUREVLYwoShH0uRqAICbg/mm3SWisuX+/fvYuHGj1GEUW3x8PDZu3AitVlvsYxXlNYmKisIff/xR7HMV1ZkzZ7B582acOHECgG6G5v379+O3337D3bt3Sy2OsoQJRTmS/HAcqodLLWkDIaISc/78eYwcOVLqMIrt+vXrGDlyZL6VUk3x+GsSExODTZs2GdQ5deoUXnrpJaOOm5iYiI0bN+ZbP+ppZsyYgb59+2LTpk04c+YMsrOz0bx5c7z11lvYsmULoqKijDpeRcEJDcqJlIxUZOSt41GN63gQUeVRrVo1DB8+XL8dHh6OUaNGYdiwYcU67o0bNzBy5EicP39ev+p1UaxcuRKbN2/Wr766f/9+3LhxA8nJybC0tCxWTOUZE4py4ta9awAACyHg4cYeCqLK4ubNm7h48SL69u2rL4uLi8PevXsxbNgwyOVyREVF4fjx4xg4cCAuXLiAiIgING7cGLVq1YJWq8WJEycQHx+PFi1a6G/LB2DwvHPnziEyMhKBgYEGC3/liYmJwalTp2BnZ4dmzZoV+AF8+vRp3Lt3Dw0aNHhqu/744w+0adNGv3Lq8ePHER0djUGDBgHQLZG+adMmdOvWDVWqVMGAAQMA6HoVDh06BCGE/jLIo4uZCSFw4cIFREZGIiAgALVr1y7Cq1y0tsbGxmLXrl3IyMjA6dOn9T0bISEhUCqV+PPPPwEAI0aM0L+2gwYNwokTJ3Dnzh307NlTf6zz58/jxo0b8PT0RGBgYL6pECIjI3HhwgW4urqiWbNmBomKRqPBiRMncP/+ffj5+aFhw4YGzz1x4gSEEGjcuDHCwsKQlpaG1q1b51sYztyYUJQTUbHXATxcx8OCq/cRVRb79+/Hxx9/bJBQXLlyBSNHjsSAAQNgZWWFU6dOYfTo0WjatCkUCgWUSiWCg4Px5ZdfYv369bC0tIQQAmfPnsWePXvQpk0bANA/r127dkhOToaDgwOOHTuGr7/+GhMmTNCfb+HChVi0aBHatGmDzMxMXLp0CevWrUOvXr30dcaMGYPt27ejXbt2uHz5MgICAp7Yri+++AJ9+/bFnDlzAACTJ09GeHg4oqKiUL16dZw7dw4jR47EgwcP9Jc8RowYgcTERISEhEAIgaCgIABAr1694OzsDLVajd69eyMlJQUODg44ePBgvrY8zZPaGhcXh+3btwMADh06hAsXLgAAbt++jezsbH08I0aM0L+2HTt2RHJyMmrXro22bdtCLpdj2LBhuHDhApo2bYpr167B3t4e27dv1yd7H3/8MT777DN06NABaWlpSElJwR9//IE6derg/v376N27N+Lj49GgQQMcP34cPXr0wIYNG/RJyYoVK3D27FlkZWXB29sbMTExuHfvHg4dOoR69eoV+bUwmvnXLCt7irpSWlm2eNMs0fCnhmLI942kDoWoXCho9UStVivSc9JL/fHoCspP888//4hHfzWvXLlSeHt7G9QJDg4WAPRt27JliwAgvvrqK32diRMnCgBi5cqV+rLRo0eLvn376rfznjd9+nSD89nY2IiYmBghhBA7d+4Ubm5uIiIiQl/np59+ElWrVtWvBP3HH38Ia2trcfXqVSGEEOnp6aJZs2YGMT7u3XffFT169BBCCJGYmCiUSqVo2rSp+PXXX4UQQixdulQ0bty4wNdk165dQqFQGBwvry1Lly7Vly1dulS4uLgUeH4hhDh16pQAoF+9uihtTU1NFQBESEiIvs6GDRuEq6trgfF89NFHBuWvvvqqePbZZ0V2drYQQgiNRiMGDRokRo8eLYQQIicnR1haWoo9e/bon3Pt2jV9jC+//LJo3ry5SEtLE0IIcfPmTeHg4CC+//57ff2XXnpJWFtbi0uXLgkhdD/3zzzzjBg3blyhr4U5VhtlD0U5cTPpHGAJeMurSB0KUbmVqc5E6/WtS/28J0adgI3S5ukVi0Eul+OVV17Rb7dt2xY//vgjxo4da1D21Vdf5XvuW2+9pf963LhxmDNnDnbs2IHx48djzZo1CAgIQGhoKE6dOgUhBBQKBWJjY3Hx4kW0bNkSmzZtQv/+/fWXHmxsbDB16lSMGzeu0Hi7dOmC5cuXIycnB4cOHULjxo3Rt29fHDhwAKNGjcLBgwfRpUsXo1+DiRMnGpxjxowZiIuLg5ub21OfX5S2GuvRiRnVajXWrVuHKVOmYNu2bRBCQAiBmjVrYvPmzfo2qFQqnD9/Ht26dYNcLoefnx+A/y4Dfffdd7C11c2YXKtWLYwaNQobN27Eq6++atD2+vXrA9AtqtmpUyf8888/RsdvDCYU5cRdcR8A4OvUWOJIiKgssra2hpWVlX47b9XlR6/Nq1QqZGVlGTxPqVTqxzEAug80Ly8vREbqVje+desW0tLS9B94eYYPHw4LC91HyO3bt9G5c2eD/bVqPXmsV4cOHZCbm4uTJ0/qk4cuXbpg/PjxEEIgODgYK1euNOIVKPg1AJCvzYUpSluNYWNjY7DadWxsLDIyMhAeHo47d+4Y1M1LnhQKBX755RfMnDkTCxcuROfOnTFixAgMGTIEDx48QEZGRr5xIXXq1MGePXsMylxcDCdALOh7b25MKMqBlPQU3LbUAJChdf1npQ6HqNyytrDGiVEnJDmvqeRyeb65HMz5wZCbm4u0tDTY2dnpyxITE/V/0Ts4OMDX1xfr168v9Biurq5ITEw0KHt8+3G2trZo3rw5Dh48iIMHD+Ljjz9GmzZtEBUVhZ07dyIhIQGdOnUqRsuMV5S2GkMmkxls29vbQyaTYcKECU+8Q6V///7o378/rl+/jl27duGVV17B7du38frrr0OhUORbqTshIaFIPTAljfNQlAOHwrciVyaDo0aLJvWekToconJLJpPBRmlT6o/HP1iMUaNGDdy/fx+pqan6sgMHDpjj5dDLG0wI6G7JjIiIQPv27QEAvXv3xvbt2/PNrfDodocOHbBz506DOSfy7nh4ki5dumDLli24ePEiOnbsCJVKhTZt2mDevHlo1KgRXF1dC3yenZ0dNBoNcnNzjWnmUxWlrcVhb2+Pdu3a4fvvv4cQosBzZGZm6pMxPz8/TJs2Df3798fx48ehVCrRqlUrg9dWo9Fgy5Yt6NChg1liLA72UJQD4RH7AAC1clWQK/gtI6pMOnfuDA8PDwwZMgTDhw/H2bNn83XJF4dCocC7776LW7duwd7eHl988QWGDRuG5s2bAwCmTp2K7du3o3Xr1njttdfg6uqKsLAw7Nu3D9eu6W5nf+2117BixQr06NEDo0aNQmhoKHbv3v3Uc3fp0gWffvopWrZsqb800KVLF8ybN++Ji0IGBATAxsYG77zzDlq1amVw22hxFKWtxfXtt9+iW7dueOaZZzB06FBkZWVh3759qF27Nr766iukpqaiXbt26NevHxo1aoS7d+/ijz/+wI8//ggAWLx4MZ555hkoFAq0atUKmzZtQmpqKmbPnm2W+IqDPRRl3IWbZ7An5xQAwNeqjsTREFFJe3wSJysrK4SEhKBFixYICQmBv78/du3aheHDh+vHR9SsWRNDhgwxOI63t7d+Toc8derUMbj9NO/4wcHByM3NRVhYGN555x388ssvBvv37duHRYsW4datWzhz5gyaNm2Ks2fP6uvY2dnhxIkT6NixI06dOoWAgADs27fPIMaCdOjQASNGjMDkyZP1Zf3798fw4cMxdOjQQl8TFxcX7N27F2q1Glu3bsXZs2cLfA0cHBwwfPhw2NgUPCA2r4cjb9xFUdqqVCoxfPhwg0sMBb3WBcUDAI0bN8bFixfx7LPP4uTJk4iOjsa0adP0g2WrVq2K0NBQeHl54ciRI0hOTsaePXv0l0jatm2L06dPw8XFBUeOHEGXLl0QFhZmEE/r1q3Rtm1bg/MGBAQY3OZbEmTi8X6XCiglJQWOjo76+6yLK/jMNgSd+QYCgIDu5RMQgHh0y2DPY/sKrmN4LEALLa5bpCLRQo5aOcAvI/bB0b5qseMnqgyysrIQERGBWrVqGQzUo/8EBQVhzJgxSEtLkzoUSfz666948cUXkZ6eXul/Rp70finqZyj7z01w495Z7FFGl9LZ5HBTa/F+my+ZTBARmUF0dDS2bduGzz77DKNGjar0yYS5MKEwgW+N5uh17xggA2SQIe//kMnw33//bePx/8vy/n10j65cBpnBcbzdAjCw02TY2TpJ1FoiqqgK65av6O7fv4/Q0FC8+eabRs2iSU/GSx5EVCHxkgdR0ZnjkgcHZRIREVGxMaEgIiKiYmNCQUQVWiW4qktUbOZ4nzChIKIKSalUAgAyMjIkjoSo7Mt7n+S9b0zBuzyIqEJSKBRwcnJCbGwsAN1CTcWZApuoIhJCICMjA7GxsXBycnriRGRPw4SCiCqsatWqAYA+qSCigjk5OenfL6ZiQkFEFZZMJoOHhweqVq1q9oWkiCoKpVJZrJ6JPEwoiKjCUygUZvmFSUSF46BMIiIiKjYmFERERFRsleKSR979tSkpKRJHQkREVL7kfXY+ba6KSpFQpKamAgA8PT0ljoSIiKh8Sk1NhaOjY6H7K8XiYFqtFtHR0bC3tzfbfegpKSnw9PTEnTt3KsyCYxWtTRWtPQDbVF6wTeVDRWtTSbVHCIHU1FRUr14dcnnhIyUqRQ+FXC5HzZo1S+TYDg4OFeIH8VEVrU0VrT0A21ResE3lQ0VrU0m050k9E3k4KJOIiIiKjQkFERERFRsTChOpVCp88MEHUKlUUodiNhWtTRWtPQDbVF6wTeVDRWuT1O2pFIMyiYiIqGSxh4KIiIiKjQkFERERFRsTCiIiIio2JhQmSE1NRWhoKG7duiV1KCbJyMhAeHg4oqKiCq0jhMDFixdx9uxZaDSaUoyueC5duoQjR44gJycn377c3FyEhYXh8uXLEkRmmqioKISHhxfYHgBIT0/H6dOnERERUcqRmebevXs4ffo0IiMjC60THx+PU6dOISYmphQjKxohBM6cOYPw8PBC6+Tm5iI8PByXL18udKriotQpLWq1GsePH8fVq1cLrRMbG4vw8HAkJycXWifvZ/HmzZslEaZR0tPTcfToUdy+ffuJ9YQQOHbsWKHfz8TERJw6dQrR0dElEKVx4uLicOTIEcTFxRVaRwiBy5cvP/F7GRMTg1OnTiE+Pt78QQoyysqVK4WNjY2oW7eusLW1Fc8++6xIS0uTOqwiuXfvnnjppZeEo6OjCAwMFM7OzqJ9+/YiIiLCoN7ly5eFv7+/cHd3FzVr1hQ1a9YUJ06ckCZoI5w5c0ZYW1sLAOLOnTsG+/bv3y/c3d2Fj4+PcHV1FY0bNxa3bt2SKNKni4qKEt26dROOjo6iRYsWolatWmLnzp0GdX7++WdhZ2cn/P39hZ2dnejevbtITk6WKOInS0hIEN26dRP29vaiefPmwsXFRTRr1izfz97cuXOFSqUSAQEBQqVSifHjxwuNRiNN0I/58ssvhb+/v3BychJNmjQpsM6hQ4dEtWrVhLe3t3BzcxMNGzYUN2/eNLpOaUhLSxP/+9//hKenp7C3txfDhw/PV+fIkSOiXbt2omrVqiIwMFBYW1uLiRMnCrVabVDv119/Ffb29sLf31/Y29uLrl27iqSkpNJqil50dLR4/fXXhYeHh1CpVOKDDz54Yv1FixYJuVxe4Pfzk08+0f8sWllZiRdeeEHk5uaWTOBPcOnSJfHCCy8IDw8PAUBs2LChwHr79u0TPj4+wtPTUwQGBoo2bdqIyMhI/X61Wi3Gjx8vrKys9O+vuXPnmjVWJhRGOHfunJDL5WL9+vVCCCFiY2OFj4+PeP311yWOrGiOHz8u1q5dq39TpKamik6dOomOHTvq62i1WtG4cWMxYMAA/S/y8ePHC09PT5GVlSVJ3EWRmpoq6tatK9566618CUVycrJwdXUVb7/9thBCiJycHNG1a1eDdpclOTk5onHjxqJnz54iNTVVCCFEfHy8wS+SK1euCAsLC/Hjjz8KIXQf2H5+fmLChAmSxPw0M2fOFJ6eniI+Pl4IIURGRoZo1qyZGDx4sL5OUFCQUCqV4ujRo0IIXRsdHR3FsmXLJIn5UWq1WsyYMUNcuXJFzJo1q8APoNTUVFGlShUxc+ZMIYQQubm5onv37qJt27ZG1SktERERYt68eSIqKko899xzBSYUa9asEceOHdNvX7p0STg5OYlFixbpy65fvy6USqX44YcfhBBCJCYmirp164qxY8eWfCMec/ToUbFs2TKRmJgo6tSp88SE4sSJE8Lb21uMGTMm3/dz9+7dQqFQiAMHDgghhLhx44ZwcXExaHdp+eOPP8TatWtFampqoQnFhQsXhEqlEgsWLNCXnT59WoSEhOi3ly5dKpycnMTVq1eFEEIEBwcLCwsLsXXrVrPFyoTCCDNnzhS+vr4GZZ9++qlwdHTMl7GXF6tWrRIWFhb65OHkyZMCgAgNDdXXuXXrlgAgtm/fLlWYT/XCCy+IGTNmiH/++SdfQvHzzz8LpVIpEhMT9WW7d+8WAMT169cliPbJfv31VyGTyZ7YgzJnzhxRs2ZNg7KlS5cKGxubMpn4vfDCC6J79+4GZePHjxedOnXSb/fr10/07t3boM4rr7xSaG+AVApLKNavXy8UCoWIi4vTl+3du1cAEJcvXy5yHSkUllAUZMCAAaJv37767blz5woPDw+h1Wr1ZV9//bWwsrISGRkZZo+1qJ6UUCQlJYk6deqIPXv2iOnTp+f7fg4bNkx06dLFoGzq1Kmibt26JRTt0+Xm5haaUIwcOfKp75PGjRuLSZMmGZR1795d9O/f32wxcgyFEcLCwtC8eXODslatWiE5OblMXDc0xalTp1CrVi39gi9hYWGQy+Vo2rSpvo63tzeqVq2KsLAwqcJ8onXr1iE8PByffvppgfvDwsJQu3ZtODk56ctatWql31fW7Nu3D40bN4aXlxcuXryIK1euIDc316BOYT+LGRkZuHbtWmmGWyQzZ87ExYsXsWDBAuzbtw9fffUVduzYgblz5+rrFNamCxcu5Gt/WRQWFgYfHx+4urrqyx7/OStKnbIsNzcXZ8+eha+vr74sLCwMzZo1M1h4sVWrVsjKysKVK1ekCPOpXn31VfTr1w89evQocH9hP4vXrl1DRkZGaYRolH379uH5559HZmYmTp8+jTt37hjsz83NxcWLFwtskzl/7irF4mDmkpCQgPr16xuU5f1iSEhIkCKkYtm3bx9WrVqFNWvW6MsSEhLg5OSUb0U5V1fXMtnG69ev44033sCBAwcKnR0uISHB4Bc4AH0by2KboqOjYWtri3bt2iElJQUZGRnIzs7GqlWr0KdPHwC6NjVo0MDgeWX5ZzEgIAAvvPACPvvsM9SpUwcREREYMGCA/sMUKPj75OrqCo1Gg5SUlHz7ypqC4re3t4dSqdR/T4pSpyybPXs24uPjMW3aNH1ZQkIC6tSpY1CvLP8srly5EleuXMHPP/9caJ3CfhaFEEhMTISNjU1Jh1lkGo0GsbGxuHPnDurWrQtXV1dERkaiXr162LhxI7y8vJCcnAyNRlNgm8z5PWIPhRGUSiWysrIMyjIzMwEAlpaWUoRkslOnTmHQoEF488038cILL+jLC2ojoGtnWWzj+PHj0adPHyQlJeHIkSO4cOECAF378u58KKhNOTk50Gq1ZbJNSqUSx44dw7Rp03Dx4kVERETgxRdfxMiRI5GYmKivU55+FqdPn46goCD8+++/+r+gbty4gZEjR+rrlLc2Pa6g+NVqNdRqtT7+otQpqxYvXoxvvvkGmzdvho+Pj768PH3fYmJiMGPGDEyePBmnTp3CkSNHEB0djfT0dBw5ckR/F0t5apNCoYBcLsdff/2Fw4cPIywsDLdv34ZWq8XEiRMB6NoDoMA2mbM9TCiM4O3tne9Wy7xtLy8vKUIySWhoKHr27IkJEybku0zg7e2NjIwMJCUl6cvUajViY2PLZBvd3d1x48YNvPvuu3j33XexevVqAMDChQuxfft2AOXv++bj4wMHBweDD9uJEyciJSUFZ8+eBVD+2rRjxw6MGDECbm5uAABbW1u8/PLL2LVrl/5yRmFtcnJygr29fanHbCxvb29ER0cb3Aaat533PSlKnbJo6dKleP/997Fly5Z8lwnK089iZmYmmjZtinXr1ul/Zxw9ehT37t3Du+++i+vXrwMovE02NjZlsqfM29sbvXr10id6dnZ2GD16NIKDgwHolh53cnIqsE3m/B4xoTBCjx49cPjwYYN7sbdu3YqmTZuWyR+ygpw+fRo9evTAuHHj8MUXX+Tb36VLFyiVSmzbtk1f9s8//yAjIwPdu3cvzVCL5Pfff8eRI0f0jyVLlgAA/vzzT323bI8ePXD//n2cPHlS/7ytW7fCzs4Obdu2lSTuJ+nVq1e+pO7u3bsAgCpVqgDQtenYsWMG96Rv3boV9evXR40aNUo13qKoUqWKvg157ty5AycnJ/1fTz169MDOnTsN5j3ZunVrode5y5oePXogLi4OISEh+rKtW7fCxsYG7du3L3Kdsmb58uWYPXs2/vzzT/Tu3Tvf/h49euDEiROIjY3Vl23duhV+fn7w9vYuzVCfqlatWga/L44cOYKhQ4fC19cXR44cQYsWLQDo2rR7926DsTtbt25Ft27d8l0OLgt69eqVL1m4e/eu/vcFAHTv3l3/Rxag+0Pxr7/+Mu/7y2zDOyuBzMxMERAQIDp16iSCgoLE3LlzhUKhELt27ZI6tCK5fPmycHFxEV27dhXBwcEGj5ycHH292bNnC2dnZ7F69Wqxbt06UaNGDTFu3DgJIy+6gu7yEEKIgQMHCl9fX7Fx40bx7bffCltbW/HZZ59JFOWTabVa0a1bN9G5c2exbds2sWHDBuHv7y+ef/55fZ3s7GzRpEkT0a5dO7Flyxbx4YcfCoVCIYKCgiSMvHBr1qwRCoVCzJs3T+zZs0d8+eWXws7OzmAUfnR0tKhataoYMmSI2LZtm5g4caKwsbER58+fly7wR4SFhYng4GAxYsQI4evrq3/vPHqH19ChQ0Xt2rXFhg0bxHfffSfs7OwMbuUrap3ScvToUREcHCzatWsnunXrJoKDgw3mnFm1apUAIN577z2D3xfh4eH6Orm5uaJZs2aiTZs24s8//xSffPKJUCgUYvPmzaXenqysLH2MNWrUEGPHjhXBwcHi3LlzhT6noLs8Hjx4IKpXry769+8vtm3bJqZOnSqsrKzE6dOnS7gF+SUmJorg4GBx8OBBAUDMmzdPBAcHi3///VdfJzIyUri5uYkZM2aIv//+W3zxxRfC2tparFixQl/n3LlzwsbGRkyaNEls27ZNDB48WFStWlXcu3fPbLFytVEjxcXFYdGiRQgLC4OrqysmTZqErl27Sh1WkezatQuffPJJgfu2bdsGFxcXALrZ1lavXo2goCCo1Wr07t0bU6dOhYVF2R/DGxoaihkzZmDLli0G2Xl2djaWLl2K/fv3Q6VSYfjw4Rg9erSEkT5ZZmYmlixZgkOHDsHW1hZdunTBpEmTDK53JiYmYtGiRQgNDYWLiwsmTJhQpv+a/+eff/Drr78iKioKVatWxYABAzBkyBCDuwMiIiLw2Wef4erVq/Dy8sIbb7yBJk2aSBj1f8aOHavvEn/U7t27YWdnB0A3NmfZsmXYu3cvVCoVhg4dajBGqah1Sku3bt2QnZ1tUObq6oqtW7cCAP73v//hwIED+Z5Xr149rFq1Sr+dlJSERYsW4dSpU3B2dsYrr7yCXr16lWzwBYiJicGQIUPylTdt2hRfffVVgc9Zvnw5wsLCDAanA8Dt27exaNEiXLlyBTVq1MD06dPz3SVRGk6cOIFZs2blKx8wYADefPNN/fbNmzfxxRdf4OrVq6hevTrGjBmT73tw9uxZLFmyBLdv30bdunXx9ttvo1atWmaLlQkFERERFVvZuxhERERE5Q4TCiIiIio2JhRERERUbEwoiIiIqNiYUBAREVGxMaEgIiKiYmNCQURERMVW9mcqIqJSFRMTg8OHD2PYsGGlfu6bN2/qp0hv0KABGjVqZPKxTGlHXFwc9u7dCwDw9PQss9NhE5VFTCiIKon4+Hj8888/T6zTvHlz3LhxA6NGjZIkodi/fz+mT5+Ovn37Qi6XFyuhCA8PN7od8fHxCAoKQmhoKAIDA5lQEBmBCQVRJZGYmIigoCD9dmhoKBISEtCzZ099mYuLCzw8PDB8+HAJItSpUqUKNm7cWOzjmNKOunXrYuPGjZg0aZLBwmtE9HRMKIgqCV9fX4MP6kmTJiE0NDTfh3dMTAz69++v37579y5OnjyJAQMG4Pz584iMjETjxo3h4+MDjUaDEydOICEhAS1btoS7u3u+8967dw+hoaGwt7dHs2bN4ODgYFTcpp7f3d29wHYMHDgQ58+fx+3bt9GgQQOzrmVAVJkxoSAiA49fKjh+/DheeuklNGzYEJaWlpDL5Th69CiWLl2Kn3/+GdbW1tBoNDh//jz27t2Lli1b6o/10UcfYfHixWjTpg0yMjJw5coVrF+/Ht27dy9yPKaev7B2tGvXDunp6bCzs8OhQ4fw3XffYezYsWZ8BYkqJ97lQURPlZGRgbFjxyI4OBiHDh3Cyy+/jNdffx1TpkzBoUOHcOTIETz77LNYuHCh/jnbt2/H119/jXPnzmH37t04fPgwFi5ciBdeeAFZWVklfv7CjjNgwAAcO3YMe/bswYIFC/Duu++a9JoQkSEmFET0VAqFAuPHj9dvt23bFlZWVnjxxRcNyq5du6bfXrNmDRo0aICTJ0/i999/x6ZNm2BhYYGYmBhcvny5xM9f2HFeffVV/XaXLl0QGxuLpKQko+Ihovx4yYOInsrOzg5KpVK/rVKp4OjoCLlcblD2aM/DrVu3kJ2djc2bNxsca/jw4QbPK6nzF/U4AIzuMSGi/JhQEFGJcHBwgJeXF37++WepQyGiUsBLHkRUInr37o2tW7fi3r17BuVRUVESRUREJYk9FERUIqZPn44dO3agdevWmDJlClxcXHDmzBkcPnwYly5dkjo8IjIz9lAQVVItW7ZEr1698pU/PiGUp6cnBg8ebFDHx8cHAwYMMCjz9fXF888/r9+2trbGgQMHsGDBAkRERODMmTNo0aIFzpw588S40tPTsXHjRpw/f75Y5y9KO5ycnDB8+HBYW1sD0E29vXHjRty4ceOJMRJRfjIhhJA6CCIiADhw4AC+//57AMCgQYNKffrvq1ev4oMPPgAAtGnTBjNmzCjV8xOVZ0woiIiIqNh4yYOIiIiKjQkFERERFRsTCiIiIio2JhRERERUbEwoiIiIqNiYUBAREVGxMaEgIiKiYmNCQURERMXGhIKIiIiKjQkFERERFRsTCiIiIiq2/wMLH0dheufzOAAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "compare(conductivity=40)" ] }, { "cell_type": "markdown", "id": "14e89708", "metadata": {}, "source": [ "Energy is conserved exactly whichever variant is used, and the lumped option combines with\n", "`\"storage along flow path\": \"averaged\"` (one well-mixed capsule for the whole bed) and with\n", "any heat transfer coefficient. On Schumann's benchmark, with sensible capsules, the lumped\n", "bed is Schumann's model itself (see the [validation notebook](validation.ipynb))." ] } ], "metadata": { "kernelspec": { "display_name": "pytes (3.14.4)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.14.4" } }, "nbformat": 4, "nbformat_minor": 5 }