{ "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-10-23T17:48:41.417925Z", "iopub.status.busy": "2026-10-23T17:48:41.417752Z", "iopub.status.idle": "2026-10-23T17:48:44.640058Z", "shell.execute_reply": "2026-10-23T17:48:44.639247Z" } }, "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.668 s\n", "lumped: solve time 9.730 ms\n", "lumped with Jeffreson: solve time 5.921 ms\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAhQAAAFzCAYAAABmY5CaAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjExLjIsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvgI3uAAAAAAlwSFlzAAAPYQAAD2EBqD+naQAAdLZJREFUeJzt3Xd8U/X6wPFPkibpbmkZLXSw994gCAoIDobKBicO8Kog3nsV7/0pTi6OK3rdorjFzRIcgCAgm5a9KQVKC3SmM834/v5IGxraQpuOtOV5v8grOd9zcvKcDs7T79QopRRCCCGEEBWg9XQAQgghhKj9JKEQQgghRIVJQiGEEEKICpOEQgghhBAVJgmFEEIIISpMEgohhBBCVJgkFEIIIYSoMC9PB1Ad7HY7Z8+eJSAgAI1G4+lwhBBCiFpDKUVmZiaNGzdGqy29HuKqSCjOnj1LZGSkp8MQQgghaq3Tp08TERFR6v6rIqEICAgAHF+MwMBAD0cjhBBC1B4mk4nIyEjnvbQ0V0VCUdjMERgYKAmFEEII4YYrdRmQTplCCCGEqDBJKIQQQghRYZJQCCGEEKLCJKEQQgghRIVJQiGEEEKICpOEQgghhBAV5tFho5s3b8ZmsxUrDwsLo2XLli5lp0+f5ty5c7Ru3VqGfgohhBA1jEcTimeeeYacnBzndn5+Ptu3b2fOnDm89NJLAOTl5TFlyhRWrVpFdHQ08fHxzJ8/n0ceecRTYQshhBDiEh5NKH777TeX7a+//prJkydz1113OcueffZZtm3bxvHjxwkPD2fJkiXceuut9O7dmz59+lR3yEIIIYQoQY3qQ/HRRx9x7bXX0qZNG2fZokWLuO+++wgPDwdgzJgxdOzYkUWLFnkqTCGEEEJcosZMvX3y5EnWrl3LZ5995iw7e/Ys586do0ePHi7H9u7dm5iYmFLPZTabMZvNzm2TyVT5AQshRE2iFFhyIC8D8kxgzQWr+eLDdslru83xHmUv8rBdsq0KjivYRpX8uSUHVHqclXL8Zd5ztWs2CFoNq/aPrTEJxccff0xwcDBjx451lqWmpgIQGhrqcmxoaKhzX0nmzZvHs88+WzWBCiFEdVMKss7Buf2QfBRMCWA663hkJRUkERlgt3o6UlET6AxXb0Jht9v55JNPmDp1Kt7e3s5yvV4PODpmFpWbm4vBYCj1fHPmzGH27NnO7cKV0oQQolaw5MGZ7RD3J5za7Egkckv/I8qF1guMgaD3BS8DeHk7bjBeRsdDZ3Rsa3WOh0Z7yUMHGo1rmVYHFJSVpNRFo0opv8IiU+U7vpznuhpE9ffIx9aIhOK3337j9OnT3H///S7lkZGRaLVaEhISXMoTEhKIiooq9XxGoxGj0VglsQohRJWwWSFuHez5Fg6uAEu2636NFkJaQMO2EBQFgY0dj4Aw8KkH3kGOh963/DdsISpBjUgoFi5cSJ8+fejUqZNLua+vL/3792fZsmVMnToVgOzsbFavXs3cuXM9EKkQQlQyaz7EfgkbX4f0+Ivl/o0cbeFNB0B4F2jQBvQ+notTiCvweEKRnJzM8uXLeeedd0rc/8ILLzBs2DDmzJlDv379+N///kfDhg154IEHqjlSIYSoRErB/p/gt/8D0xlHmU8IdLwdOo+HiF5S0yBqFY8nFH/99Rf9+vVjwoQJJe4fNGgQf/zxB2+//Tbbtm2jU6dOfP755/j7+1dzpEIIUUlMifDz43D4Z8e2fxgMmAXd7wKDr0dDE8JdGqUuOyanTjCZTAQFBZGRkSHTdgshPMeaD9veh/Uvg9nk6EA58O8w4DHQe1/5/UJ4QFnvoR6voRBCiKvC0dWw6h+QesKx3bg7jH4LGnXwbFxCVBJJKIQQoirlmeDXpyDmc8e2X0MY8jR0nQLaGjVZsagDlFLYlA0vbfXf3iWhEEKIqhL/F/z4AGScBjTQZzpc9xR4S9OrqBiL1cbpzCQOpxzleMZx4jJOEJ8Zx+msOO5q+yAP9bi72mOShEIIISqb3QZ/vgrr/+OYsjo4Gsa8C02v8XRkopJZbHbyLDbyLIXPBa+tjte5+TbyrI59Zqsdc+Gz1Y7ZasNsKfrseJ1XWGa1k2exkmtPJpdELLqzWHXnUPoktIbzaHTmEmNac2IPD/UocVeVkoRCCCEqk+mso1bi5AbHdpdJcNMrYAzwbFxXEYvNTk6+42aek28lt+Amb7bYyLPayM0vuPlbL00EHNu5RRODwtcFx+bm25w3/TyLDau9ssY12NAYUtEZzqM1nkdrPIfWcAGt33k0WovzqKKNZEppsefXh/xG6KxheNnCMapwru/Zs5JiKh9JKIQQorIcWgnLHoacFDD4w83/hS4lD4m/mtnsilxLwc0+31bwujABcJTnWQpfXyzPtdjIzbcWeX2xPKegPM9iw2LzzOBFb70Wb70OH70Ob70Oo9fFbaNei7eXDr2XDavuPGZNIjnqLFkqAZM1gQzrWeyUvBaLl0ZPuF8UUf5NaRrYnObBzWldrxUtgqPxMxjRamvGfCWSUAghREXlZzs6Xu78xLEd3gXGLoLQFh4Nq6KUUpitdrLMVnLMNsdzvrXg2bGdXeR1jtlKlvnizf3iTd+ROOQUbJut9mqJX6fV4KvX4W3Q4a3XOm/03l6OG7xzW1/0xn/JsQWJgLfB8b5Lk4bCbaOXFk2RichyLDnEmeI4kX6CExknOJ5+nLiMOE5nnsaWbysxXh8vH5oFNaN5UHNaBLdwPjfxb+KRTpblVfMjFEKImixhF/x4P6QcAzTQ/xG4/t+Ohbg8JN9qJzPPQmae1fEwF3ntLLc4b/7ZZivZZhvZl7zOybdhq7Qq/eI0GvDR6/A1OG7OvgYdPgYvfAvLDDrnax+Dl+NZr8PH4CgrLPdxHlNQrvfC26DFoHO9yVcFU76JQ2kniMuI43j68YIOknEkZCWU+p4AfQDNg5u7JA7Ng5sT7heOtrQF2GoBSSiEEMId1nzY8CpseM2xbHhAY7jtfWh2rdunVEqRZylIBszFE4CLSUGRbXPx8qqoAXDcwL3wNxY+e+Fr1OFn9MLPUPjs5Xg26gpu8l74GLT46L2KJAAX93nrq/6GX1lS81KdtQzOxCE9jvO550t9T4h3iCNZKEgYCpOHBj4Nas11l4ckFEIIUV4JO2Hpw3D+gGO7/Ri45XWUTz0y8yxk5FhIz7GQnptPes6lCYHj2ZRnJauEmoPK6+TnSAICvL0I8Na7Phu9CPD2wt+ox68wKSiWGFws99Hr0NWQdvqqZLPbOJt9lriMuGKPNHNaqe9r6NvQtbahIIEI8Q6pxug9TxIKIYS4hM2uyMwrTAospOfkk5FrITMzk7aH3qZ7whdosWPSBvO+/wxWxfch/bWdZORaKqWJQKMBf6MXgc5EwJEM+Bu9XBKDwFLL9fh7e10VSYA7Cvs3XJo0nDKdIt+eX+r7mvg3KdZM0TyoOQEGGcEDklAIIa4CNrsiLSeflKx8UrLNpGTlk5qdT0qWmZTsfOd2craZ1GxH8nDpKkc9NId5Wf8BLbSJAPxku4bn8u4gLScQyHE51luvJdjHQLCvnkAfPYHe+oKbv2ttgX9BWeAl5b56XY3puV9bKaU4l3POJWE4aTpJXEYc53LOlfo+g9ZAdFA0zQKb0Szo4qNpYFN89bJw2+VIQiGEqHXsdkV6roXUbDPJpSUHBdup2fmk5eQXSxDKws+gI8zHzqN8xUjzCrQoMrzqs7rFk6Q0uZ45PgaCfPUE++gJ9nUkEEE+erz1usq/aFEis81MvCmekxmOZKGw5uFkxklyrDmlvi/EO+RiwlCQPDQNakpjv8botPL9c4ckFEKIGsFis3M+00xSRi4XMktPFFIKahHcaVmo56snxM9AqL+RUD8Dof4GQvyM1Pc3OMr9jIT6G6jnayDIR4/h1AZY9gikxztO0O0Ogm54gdt9giv12sXlWewWEjITiDfFE2+K51TmKU6aTnLKdIqk7CQUJf8w6DQ6IgMii9U0NAtqRpAxqJqvou6ThEIIUeWyzVYSM/I4Z8pzPidl5LmUpWSby12LEOjtRX1/Y0FiUCRR8DMQ4m+kvp+BEH9HolDPV4+XroxD8jLPwfK5sPsrx3ZQJIx8A1oOKV+AosxsdhuJ2YklJg1ns85iUyXP3QCOYZjNgl1rGpoFNSPSPxK9Tl+NV3F1k4RCCOE2u12RmpNP0iXJgstzRh6Z5pJnALyUXqehUaA3DQOMZapFMHhV8ph9az5sex/WzYf8TEdZr/tg6FyZOrsSWOwWkrKSOJ11mjOZZzhlOuVIIDLjOZN5BovdUup7fbx8iAqIIjowmujAaKICo2ga2JSowCjqGevVyWGYtY0kFEKIUmXkWIhPzeZ0ai6JGbkkZeSRZLqYLJw3mcm3lW3OgwCjF42CvAkP8qZRoDdhgd6EBRV5DvImxNfguc6Ix9bAL09C8hHHduPujjU4IjyzLkJtpJQiw5zBmawznMk8c/G54HVidiJ2VfrPi0FrIDIg0iVpKHxdV+duqEskoRDiKqaUIiU7n/iUbOJTcjiZkkN8SrbzOT2n9L8YC2k0EOpndCYK4QXJQeHrRgUJg7+xhv53kxoHv/0bDq1wbPvWd9RIdJ0C2to7a2FVybHkkJidSGJ2IgmZCcWShyxL1mXfb9QZifCPoElAE5cah+jAaBr5NpIOkbVYDf0NF0JUFrtdcT7TzMmUbGeycColp2A7h6wrNEc0DDASFeJLeLAPYYHGgkTBh7AgI2FBPjQMMKIva9+EmiQ/GzYugE1vgM0MGh30eRAGPQFXaadLq91Kcm6yI2HISnQmDueyzzlfm/JNVzxPQ5+GRAREOB7+ES6vQ31Ca/X00qJ0klAIUQfY7Iqz6bkFtQzZLrUMp1JzyLOUXs2s0UDjIB+iQ32JDvWjacGzY9sXX0Md+2/CZnEs4rX+ZcgumDa52SC48WVo2NajoVUVpRTp5nQu5F4gOSeZC7kXHK9zk7mQ43idlJ3E+Zzzl+38WChAH0Ajv0YuyUJkQCQR/hE09m+Mt5d3NVyVqGnq2P8UQtR9adn5HEw0cSDRxMHETA4mmjh2PuuyfRl0Wg2R9XyKJAoXE4eIej5Xx7wJdjvs/xHWvgBpcY6yek1h2HPQbpQjs6ol7MqOyWwizZxGujmd9Lx00s3pLttp5jRSc1OdyYPVXraOsV4aLxr5NSLcL5wwvzCX58LX/gb/Kr5CURtJQiFEDWWzK+KSszhQkDQUPs6ZzCUeb9BpiQr1dSYKTUN9iSp4bhzsUzubJSqDUnB8Dax+FpL2OMr8GjiaNrrfBV6GagxFkWvNJcuSRVZ+luO54HW2JZvM/EyyLdnFyy2ZZOcXPFuyybZku/X5wcZg6vvUp6FvQ+r71KeBTwMa+Dagvk99Z9IQ6h0q/RiEWyShEKIGyMi1cMiZNGRyMMnE4aTMUleNjA71pV1YIO3CA2kXHkC78EAaB/vI2g2XOrMDVs+Fkxsc24YAuGYm9J0BxrL/la2Uwmwzu9zkS00KCm7+zvIix2Vbsi87yqG8AvQBBHsHE2x0POp513O+DvYOJsQ7xJE0+DiSBpmTQVQlSSiEqGYpWWb2nMlg95l09iU4koiE9NwSj/U16GgTFlCQOATSPjyANmGBNXfERE1xehv8+Srq6K/kayBbbyS7yySyuk4iy0tH1rntzpt9tiW7WGJQWF60xqCsTQZlodVo8df7Ox6Gi89+ej/XsiLlAfoA/AyOZ1+9L0GGIEkQRI0i/yuJKmO12fntwDkOJ2Vi1Gu5tlUDOjQOvKrGkmebrexNyGDPmXR2n8lg9+l0zqSVnDw0CfZx1jYUPqJDfK+6RaIKmwUKq/ZdHtZssvMLni3Z5FhynDf/HEs2WVmJZJvOkGPJJVurIbtpJNbCn7eUtbBmbYVi06DBT+/nuMEbAkpMAPwMfsWThYLyAL3jPT5ePlfV74G4OtSIhCIpKYmNGzfi6+vLddddh4+Pj3Pfpk2bOHjwoMvxISEh3HbbbdUdpiiHg4kmZn+7m4OJF4eYvfzLYXo3DeGfI9rQs2mIB6OrOllmK9vjUvnreDJ/HU/hYKKpxDUnWjTwo0tEMB2bBNG+cSDtwgIJ8q29f23a7DayrY4bfOHNP8uS5dwu+vrSBME1KXAcU9raDGWmL/5fm4+Xj+PGXiQZKFNSUGSfr95XhjwKUQqNUu6swVd55s+fz3PPPcfgwYPx9fXl2LFjLF26lKioKACmT5/Or7/+ypAhF+fQb9KkCc8++2yZP8NkMhEUFERGRgaBgYGVfg3C1b6EDKYs3EpGroUgHz03dgwjNTufdUcukF/QJ+C2bk148sa2NAys3cPL8iw2dp1KY/PxFDYdS2b3mQxsl2QQ4UHedI4IoktkMF0jgukYEUSgt2eTB4vdQo4lh1xrLjmWHHKsOc7nbEu2y/aVkoIcq+M8la2wNsBX7+tMBHz1vvh5+eFv8MfXywc/0zn8T2/DN/0Mfkrhhw6/FkPw6zIZv5AWztoEXy9f6WgohJvKeg/1aA3Fd999x7/+9S/WrFnDoEGDADhx4gQ2m+s46B49erBw4UJPhCjKKSE9lzs/3kZGroXuUcF8eGdPQv2NACRl5LFg9RG+2XGaH2MS+O3AOWYNbcVd/ZvWmhEIVpudPQkZ/HXMUQOxIz7NmSQVigzxoX/z+vRvGUrf5qE0qmDSZLFZit/kL7nhFzYRlLQ/15J7MTEoqBG43JoJFeGl9XLcxL388DMUPJeUFBTUCrgkCAWvC5OAUpsF8rMh9ivY+h6kHHOU6X2h573Q72EIDK+SaxNCXJ5Hayh69uxJixYt+Oabb0o9Zvr06Rw8eJAHH3yQoKAgevbsSaNGjcr1OVJDUT3yrXbGvb+Z3afT6dA4kMUP9CWghL/Ed59O5+ll+9l9Oh2AVg39mTuqA9e0rF/NEV+Z3a44lJTpbMLYFpdabGbJhgFG+rcIpV/zELo29SU0QOO8iedac4vd1C+94Wdbsi+7vzI7A15Kr9Xjq/fF16vgUfDaR++Dr5ev86/7KyYFBa8NuiocgpmRANs+cExKlZfuKDMGQa9p0O9v4Ffzfn6EqAvKeg/1WEKRk5NDQEAA77//Pv379ycmJobGjRtzzTXXYDBc/E9p+vTprFixgv79+5OQkEBMTAyvvvoqDz30UKnnNpvNmM0Xx+qbTCYiIyMloahir/x6iLf/OE6Qj54VjwwgMsS31GPtdsV3O08z/5fDpGbnAzC4TQOeGNGWduFV9z2y2CzOG31pj5TsLPYnXeDw+VTiU9Mw2/LQaPNBm49Gm4+XlxV/bztGgxWt1oJF5ZFrzcVsK3l+iMpi0BqcN3MfL58SEwFf/SWvS9pf5LlWjBI4sxO2vA37l0DhLI71mjmGfnadUq7hn0KI8qvxCUVCQgIRERHcdNNNHDt2jB49erBz507sdju//vorzZs3B2Dz5s307NkTvd7xH9+HH37IjBkz2LZtG927dy/x3HPnzi2xj4UkFFVn9+l0bn1nE3YF707pzo2dylbtnJFj4fXVR/hiSzxWu0KjsTO8Y33u7B9OdAM9uZaSb/olJQV51rySk4Qi57CqqvtrvygfLx/no/DmXfjXfmkJQdH9RWsJCo/Xa2vBzb+yWHLhwDLY8RGc3nqxPHoA9HsIWo8A6RMhRLWo8QlFSkoK9evXp2vXrmzduhWDwYDVauWaa66hcePG/PTTT6W+t2HDhsyePZsnn3yyxP1SQ1G9cvOt3PT2z8SbEujVQsPY3vWdN/eiN/k8W16x8jyb4zk7P5dsSy528qslZp3GC73WiLIbyLd4YbV6oewGsBtQdgOBRl8igoNoFlKPqHrB+BkKkgAvX5dkwfnQX3ztrfOWIYHuunDY0aQR+9XFZg2tHjre7kgkwrt4Mjohrko1vlNmaGgoISEhDBs2zNnE4eXlxQ033MAXX3xx2ff6+PiQmppa6n6j0YjRaKzUeIUrpRS7L+zmh6M/8MvxP8gLycAvBA7Y4LnNlfQZ9oKbvDJg1HkTYPAjwOBLsI8fob4BBBTc5Atv5r5evnjrvF1u7j5ePhg03qTnaIg7n8+BhDz2nM7lUGKOy3BOH72Oa1rW57q2DbiuTUMaB/uUHpioXJY8OLDUkUic+utieVAkdL/T8QgI81h4Qoiy8egoj9GjR7Nnzx6Xsj179tCsWTMAbDYbZ8+eJTIy0rl/8+bNnD59mn79+lVrrOKi4+nHmb9tPpsTL2YOSmkIMTagWXAEQcYgvL28HTd4L2/HTd7LB2+vS5513q7HFSm/YLLz3Y4Elu0+y6nUHLKBS1PI+v5G6vsbaBBgpL6/kUBvL2fNgCnXwoUsM+dMWZxMPl/iwllRIb5c37Yh17VtSJ9mIVfHAlk1hVKQtNdRE7H764u1ERqdozmj5z3Q4npp1hCiFvHoKI8zZ87Qp08fBgwYwDXXXMPWrVtZsmQJa9asoW/fvlitVrp27Ur//v3p0KEDp06d4oMPPuDmm2/mq6++Qqst21BDGeVROZRSfHfkO+Zvm0++PR+9Vo82pztp5zozsk0/FkzoVSWfeyYth83HU9hzJoNj57M4diGLC5nl6wDprdfSJiyQHlH16BHteIQF1e45MGolUyLs/RZ2fwPn918sL6yN6DYVAht7Lj4hRDE1vg9FoeTkZBYtWkR8fDxRUVFMnjyZiIgI536z2cw333xDTEwM9erVY+DAgVx33XXl+gxJKCrOruy8tuM1PjvwGQADmgwgIGs8327JoVGgkd8eG0SQT/V1GszIsXA6LYfkLDMpWfkkZ5nJMltRChSKAG+9swajRQN/mgT7XHVTWNcY+dlw6GdHTcSJdVC4OJbOCG1udIzUaDlEaiOEqKFqTUJRHSShqBilFP/Z9h++OvQVADO7z6Sj7xgmLdyKUrDonl5c16ahh6MUNYolF47+Bvt/giO/giXn4r6oftB5AnQYAz71PBaiEKJsanynTFF7LNq/yJlMPH/N8wyLvIUb39iAUjChZ6QkE8LBaoZjqx1JxOFVkJ91cV+9ZtBlInQeDyHNPRejEKLKSEIhLmtL4hYW7FwAwBO9nmBMyzH835J9nErNoXGQN/++pZ1nAxSelWeCY7/DoZWOGgnzxcXgCIpy1EJ0vA3Cu4IMpRWiTpOEQpQqJTeFJ/58AoXi9la3M7X9VFYfOMfnW+IBmD+2c4lTa4s6zpQIh1c6+kXE/QlF1wUJaAwdbnU8InpKEiHEVUQSClGq+dvmk5qXSqt6rXiy95OcM+Xxj+93AzBtQDMGtmrg4QhFtVAKkvbA0d8diUTCTtf9oa2g7c2OR5OeUMbRV0KIukUSClGi9afXs+rkKnQaHS9c8wJ6rZHHvtlKWo6FDo0D+eeINp4OUVSlnFQ4vhaOrYHjayDrXJGdGojoBW1vgjY3Q4PWHgtTCFFzSEIhirHYLLyy4xUA7mx/J+1D2/POumP8dTwFH72ONyd1w+glQ/zqFJsVEmMdnSqPrXbUQqgik4Hp/aDZtdBmBLS+EQLKt+KvEKLuk4RCFPPtkW+JN8UT4h3Cg10eJOZUGq/9dgSAZ0d1oEUDWd2x1rPb4dw+Rx+Ikxvg5CbIz3Q9pmEHx/wQLYdCVF/wkunshRClk4RCuMix5PDBng8AeLjbw9hsBh5dvBWbXXFL53DG9Yy4whlEjWS3Q/JhOLkR4tY7nnPTXI/xDoJmg6DVMGgxBIKaeCZWIUStJAmFcPHdke9IzUslwj+CMS3G8Pfv9nE6NZeIej68eGsnWUWztrDkQsIuOL0FTm11LAFeuF5GIYM/RPeHpgMdzRlhnWS2SiGE2yShEE551jwW7VsEwAOdH+DX/RdYGnsWnVbDGxO7VevU2qKcsi4UJA9bHMnD2VjX4ZwAXj4Q2ctRC9HsWmjcDXTyPRVCVA5JKITTyriVpOSlEO4XTv9GN3DTG46lpP82uAU9omWK5Bojz+ToQJmwC87ugoQYyDhV/Dj/RhDZx9H/IbIvhHeWBEIIUWUkoRCAY72Ozw98DsCUdlN44ecjpGTn06ZRAA9f38rD0V3FLLmQtK8gcShIIJKPAiUswdOwfZEEog/UayoTSwkhqo0kFAKArUlbOZZ+DB8vH4Ks17B892F0Wg2vjOuMwUsmKqpySkFmoiN5OLe34HkfpBxzHb5ZKCjS0WTRpDs07g6Nuzo6VQohhIeUOaG47777ynXihQsXljsY4TnfHf4OgJuajeQ/Pzum1n7g2uZ0jgj2YFR1VJ4Jko/AhUNw7sDFBCI3teTj/Ro4kofG3S8mEP4yS6kQomYpc0Lx0UcfMWHChDId+80330hCUYuk5qWy9vRaAKxpvUky5RFRz4eZQ6Spo0KyUxxDNS8cgguHLz4yz5Z8vEYH9VtBo44Q1hEadXI8+zeSpgshRI1XriaPxYsXl+m4b775xq1ghGcsP74cq91Kq6B2fLfZUb3+r5va4a2XIYSXpRTkpEBqHKSeuPhIK9jOSSn9vf5hjimrG7RzJA1hnRyv9d7VF78QQlSiMicUy5cvL/NJy3Os8Lwlx5YAoDJ7Y7ba6ds8hBEdwzwblKcp5ViKOzsZTAlgOgsZZxzPpgTISID0eNfluksSHAX120CDwkdbqN8afIKr5TKEEKK6lDmhuOWWW8p80vIcKzzraNpRjqUfQ6fxIuZgNFoNPH1Lh4sTWFnNjlkVE3ZCxmmw5oPdevFREnXpCARVyr7Syt15zyXvL+97lAJrnmP2yNw0yE0HZaNMAiMgpFnBo7njUa8ZhLYAg1/ZziGEELVcuZo88vLy+Pnnn7n99ttL3P/DDz9w88034+0t1ba1xS8nfwHAx9qedLsvE3pH0b5xINgssO0D+PPV0jsLXg30fhDY2DENdWDho7HjOTjKMTRTmimEEKJ8CcXChQs5duxYqQnFn3/+SVJSEn/7298qJThRtZRS/BLnSCguJLXH4KV1dMTMSYVv7oD4jY4D/cOg+WDHX956H9B6FTy0oCk6pPSSjoMuHQkrYV+xjomVvM/LCD4h4FPv4kOSBSGEKJNyJRSffvop77//fqn777rrLqZPny4JRS1xIOUApzJPoVF6rJntuKd/FGHGfPj4Fji/HwwBMPxF6DoFdDJliRBCiNKV6y5x5MgRWrduXer+1q1bc+TIkQoHJarHqrhVAORntsPby4cZg5rBD3c7kgn/RnDnUmjYzrNBCiGEqBXKNQWizWbDai2lIx6Qn5+PzVbGjmzCo+zK7uw/YTV15s5+TWl4ZDEc/RW8vGHS15JMCCGEKLNyJRQdOnRg5cqVpe7/5Zdf6NChQ4WDElUv9nws53LOoWxGDPntmd7DD36f69g5dC406eHJ8IQQQtQy5Uoo7r33Xh577DHWrl1bbN/atWt57LHHyj1Ft/CM1adWA2DNbM+Ens0J2fEGmDMcUzz3fsDD0QkhhKhtytWH4oEHHmDdunUMGTKErl270qZNG5RSHDlyhNjYWKZMmcK0adPKHcRff/3F2rVr8fX1ZcKECTRp0sRl/4ULF1i8eDHnzp2jU6dOjB07Fp1OZnF0l1KK3+IcSaE9uyP3dfGBTx0rjTLsedDK11YIIUT5lKuGQqPR8NVXX7F48WKioqLYvXs3e/fuJSoqim+++YYvvvji4oRIZWC327nnnnsYNWoUaWlppKWlcdNNN3Ho0CHnMcePH6dTp04sWbIEm83GnDlzuOmmm6SvRgWcyDjBudwElF3HkKYDaXJgIdjMENUPmg7wdHhCCCFqIY1SxaYorDZvvfUWTzzxBLGxsbRq5ViIymQykZubS6NGjQAYO3YsSUlJ/Pnnn2i1Wk6ePEnr1q1ZtGgRU6ZMKdPnmEwmgoKCyMjIIDAwsMqup7Z4fdt7fHzwbaxZbfju5rfp+HUvR3PHlO+h1TBPhyeEEKIGKes9tFw1FJXtrbfeYvLkyc5kAiAwMNCZTFitVn7++WemTJmCVusItWnTpgwaNIglS5Z4IuQ6YenR3wGIMvako2m9I5kIjoIWQzwcmRBCiNqqXAnF6tWradmyJS1btmT16tUV+uDMzEwOHz7MwIEDWbJkCf/3f//Hu+++y7lz55zHxMfHk5eXR4sWLVze26JFi8vOd2E2mzGZTC4P4ZBgOk+K5SgAD/a6BXZ95tjRdapj5kshhBDCDeW6g8yYMYPvvvuO7777junTp1fogzMyMgB47bXXePfdd9Hr9axYsYJWrVqxbds2AHJycgCKVbEEBQWRnZ1d6rnnzZtHUFCQ8xEZGVmhWOuSd7YtB41Ca4lgdFQQnNwAaKDrZE+HJoQQohYr1yiP5ORkoqOjUUqRnJxcoQ8OCAgAICQkhF9//dVZfuONNzJnzhzWrFmDv78/AOnp6S7vTU9Pd76/JHPmzGH27NnObZPJJElFgTXxf4AXdK8/AN2RgjlFovpCsHx9hBBCuK9cNRQPP/wwo0aNYtSoUTz00EMV+uCgoCAaN25Mr169XMp79erFiRMnAIiKisLPz4/Dhw+7HHPo0CHatSt9Fkej0UhgYKDLQ8CehAtkaQ8AcF+PW+DQCseOtrLcvBBCiIopVw3F888/z7Zt27Db7fTt27fCHz5x4kTWr1+P3W5Hq9WilGLdunV07twZAJ1Ox2233cann37K9OnTMRqN7Nu3j40bN/L9999X+POvNu9sXo1Ga8FAMP1DIuDkJseOtjd5NjAhhBC1nkeHjWZkZDBkyBBsNhv9+vVj27ZtXLhwgbVr1zo7YiYmJjJw4ED8/f3p3r07K1as4IYbbuCLL74o8+fIsFHIybfS593HIHgd/RqM4IMmPeGnB6BhB3joL0+HJ4QQooaq9GGjr776apk/vKzHBgUFsXnzZp577jnat2/P008/zeHDh11GdYSHh7Nnzx6eeuopunXrxrfffluuZEI4rNqbhM3b0XQ0us1gOPGHY0frGzwXlBBCiDqjzDUUGo2GslZmlOfY6iA1FDDp49/Zp3N0VF0/fh0h7w4EUwLcsQRaXOfZ4IQQQtRYZb2HlqsPRdeuXSsal/CAlCwzO89vxRgOLQLbEJKT7kgmdAaI7OPp8IQQQtQBZU4oXnnllTKfdOrUqW4FI6rGyr2JaH0cE4ENjhoAJ9Y5dkT0BoOv5wITQghRZ5Q5ofj73/9elXGIKvTz3kR0/o7ZMfs37g8b3nbsaD7Ig1EJIYSoS2Su5TouI9fCzsQDaL2yMOq86dqwK5x2zERKVD+PxiaEEKLukISijttw9AIUNHf0DuuFIesCZJ4FjRaadPdwdEIIIeoKSSjquLUHz+PldwyAfo37QcIOx45GHcDg58HIhBBC1CWSUNRhdrvijyNJ6HxOAtA7rDec2e7YGdGr9DcKIYQQ5VShhCIlJYXt27dXViyikh1KyiTDdhKNLp8AQwCt6rWCMwU1FE16ejY4IYQQdYpbCUV6ejpjxoyhfv369O7d21k+ZswYNm/eXGnBiYrZciIFnW8cAD0a9kBrt8HZGMdOqaEQQghRidxKKP75z39iNpuLrQL60EMP8fzzz1dKYKLitsal4FWQUPQM6wnJR8GaB4YACG3p4eiEEELUJeWaKbPQihUr2LJlC1FRUS7lvXr1Yv369ZUSmKgYu12xNS4ZXURBDUWjHpCw17EzrCNopfuMEEKIyuPWXSUtLY169eoBjnU7CmVlZaGVG1WNcPR8FibbaTS6PHy9fGkb0haS9jh2hnXybHBCCCHqHLfu/j169OCnn34CXBOK+fPn079//8qJTFTIjvhUZ/+Jbg274aX1gnP7HDsbdfRgZEIIIeoit5o85s2bx0033cSff/6JUor/+7//49dff2Xv3r3S5FFD7D2T4Uwoeob1BKUgqbDJQ2oohBBCVC63aigGDhzIhg0byM3NpVWrVixevJhmzZqxefNml1EfwnN2n0m/OMKjUQ/ITIKcFNDooGE7D0cnhBCirnGrhuK9995j+vTpfPnll5Udj6gEeRYbR9OO4xOcjUFrpGNoRzj+h2Nn/Vag9/FsgEIIIeoct2ooHn30UWw2W2XHIirJgUQTGu+TAHRp0Bm9Tg/nDzh2NmzvucCEEELUWW4lFO3bt2fnzp2VHYuoJHvPZKD1OQXgWF0UINmxQBgN2nomKCGEEHWaW00e9957LxMnTuTpp5+mffv2GAwGl/1du3atjNiEm/acyUBXkFB0btDZUXihYBKyBq09FJUQQoi6zK2EYubMmQDcc889Je5XSrkfkaiwg+fOoQs8D0Cn+p0cIzwKayjqt/FgZEIIIeoqtxKKCxcuVHYcopLY7YqTmYfQBUIjn8aE+oSC6SyYTY4RHqEtPB2iEEKIOsithKJ+/fqVHYeoJAnpuVgNJ9EB3Rp1cRQWNneENAMvo8diE0IIUXe5lVCsXr36svuHDh3qVjCi4o6cy3T2n+jasCChkOYOIYQQVcythGLYsGGX3S99KDzncFImWu/TAHSuLx0yhRBCVA+3ho0qpVweNpuNQ4cOMWjQIJnsysP2nDuO1isHHXrHgmAAKUcdz/UloRBCCFE13KqhuJRWq6VNmzYsXLiQW2+9lcmTJ5fpfe+88w4rV650KYuOjubtt98u1zHiosPp+8AHIvxaOSa0Akh1TMFNiHTIFEIIUTUqJaEo1KBBA+Li4sp8/J49e0hPT+fJJ590lgUGBpb7GOGglOK8+QhaH+jcoGABMEseZJxxvJYRHkIIIaqIWwnFmTNnipWlpaXxyiuv0L59+aZ2DgsL45ZbbqnwMQIuZJpRxngA+kd0dxSmxwMKjIHgG+q54IQQQtRpbiUUkZGRJZa3bt2ar7/+ulzn2r17N+PHjycoKIiBAwcydepUtFptuY8RcCw5Da13IlBkyGjqCcdzSDPQaDwUmRBCiLrOrYTi4MGDxcrq1atHo0aNynUevV7P9ddfz6BBg0hISOCf//wnixcv5ueff0ZTcPMryzGXMpvNmM1m57bJZCpXXLXVjrMH0Wjs6JQ/jf0aOwqdCUVzzwUmhBCiznMroZg4cSKxsbEl7uvatWup+y714osvuvSHGDFiBF27duWnn37itttuK/Mxl5o3bx7PPvts2S6mDtmf7FhRtJ5X04vJliQUQgghqoFb7Qa7d+8usdxms7F///4yn+fSzpWdOnUiOjraZSXTshxzqTlz5pCRkeF8nD59uswx1WbxWY7hoZF+rS4WSkIhhBCiGpSrhmLLli0lvgaw2+1s3ryZqKgot4NRSpGeno5er6/QMUajEaPx6ptiOiU/DnTQNqTdxUJJKIQQQlSDciUU/fr1K/F1obCwMN54440ynctisfDZZ59x7733OqvnX3rpJTIyMhg9enSZjxEOVruVPI1j9E2v8I6OQpsF0gtqZyShEEIIUYXKlVBkZmYCjsXBkpOTXfbp9fpy1QrodDp27NjBM888Q6tWrTh9+jSZmZl8/vnndOvWrczHCIeDycdBa0HZDPSMKJgRM+MMKBt4eYN/+TrMCiGEEOVRroTC398fgLy8vAp/sFar5d133+WFF15g37591KtXj9atW+Pt7V2uY4TDljN7AdBYGlPPtyCxS3csEkZwlAwZFUIIUaUqdaZMd4SGhjJo0KAKH3O1O5h8DAB/bcTFwqIJhRBCCFGF3BrlkZeXxzPPPEP79u3x9/fH29vb5SGqX3ymo/NlfUORScckoRBCCFFN3Eoonn76aX788UeeeuopsrOz+eKLL3j00UfRarX84x//qOwYRRmcy3UkD439ml4slIRCCCFENXGryeObb75h+fLldO7cmTvuuIPbb7+dsWPH0qNHD956663KjlFcgcVuwWRLAqB5UJHRHJJQCCGEqCZu1VCcPn2aDh06AODn50dGRgYAN998Mzt27Ki86ESZnDadRmFD2Qw0r9fk4o7ChCJIEgohhBBVy62EQimFTqcDoGXLlvz+++8A7Ny5U5YW94DjGccBsOc3IjzIx1FozYfMs47XUkMhhBCiirmVUAQFBTlfP/bYY9x555307NmTG2+8kQceeKDSghNlcyLd0SHTbm5AWFDBkFFTAih7wRwUDT0YnRBCiKuBW30o0tPTna/vuusuWrRowZYtW2jbti233HJLZcUmyuhYehwA9vyGNAosGGXjbO6IlDkohBBCVDm3Eor33nuP6dOnO7cHDBjAgAEDKi0oUT6nMhIA0NtDCfAuWOMko2DK7eDIUt4lhBBCVB63mjweffRRbDZbZcci3HQuJxGAEGOR6bULkgwCm5TwDiGEEKJyuZVQtG/f/rLLh4vqY7PbSDNfAKCRb9jFHSbHQmEERZTwLiGEEKJyudXkce+99zJx4kSefvpp2rdvj8FgcNnftWvXyohNlMGF3AvYsaOUlohAqaEQQgjhGW4lFDNnzgTgnnvuKXG/Usr9iES5JGU7JrRSlmDCGvle3GEqSCiCJKEQQghR9dxKKC5cuFDZcQg3nc1yzDVhtwTRIKDI8vHOGgpp8hBCCFH13Eoo6tevX9lxCDclZjs6ZCpLMKF+BU1PeRmQn+l4LTUUQgghqoFbnTIBkpOT+eijj/jXv/7lLNuyZYuM/qhmhQmF3RpMvcKEorB2wjsYDH6eCUwIIcRVxa2EIjY2lnbt2vHaa6/x0ksvOcs/+eQTPv3000oLTlxZ0T4UIb4FCYWz/4Q0dwghhKgebiUUjz/+OLNnz+bAgQMu5dOnT+f111+vlMBE2ThrKCzB1PMrnNSqYMiojPAQQghRTdzqQ7Fjxw6WLFkCgKbItM6tWrXiyJEjlRKYKJvErKJ9KIqs4wEQ2NhDUQkhhLjauFVDodPpyMzMLFZ+6NAhQkJCKhyUKBuL3UKmxfF9MGgC8TE4VoDFVLDKqHTIFEIIUU3cSihuuukm5s6di91ud9ZQxMfHM2PGDEaOHFmpAYrSZZgzAFBKQz3viyvAOhOKAKmhEEIIUT3cSihee+01Nm/eTHh4OHa7nU6dOtGyZUvy8vKYN29eZccoSlGYUGD3JsTP++KOrHOO54Cw4m8SQgghqoBbfSgaNWrEzp07Wbp0KTt27MBut/PUU09x++23F5uGW1SddHM6AMrqR4hfka97pqNfBQHh1R+UEEKIq5JbCQWAwWBg3LhxjBs3rjLjEeXgTChsvtQLKEgorGbITXO8lhoKIYQQ1cTthOLEiRO8/fbbHDx4EHCsQPrwww/TtGnTyopNXIGzD4XN92INRaZjXgp0BvCp56HIhBBCXG3c6kOxfPly2rZty/r164mOjiY6Opp169bRpk0bVq5cWdkxilKk5TlqIpTNl3qFk1oV9p/wD4MiQ3qFEEKIquRWDcXjjz/Os88+y5w5c1zK582bx2OPPcZNN91UpvM8/fTTfPXVVy5lrVu3LpaULF26lDfffJNz587RqVMnnn/+eVq2bOlO6HWKaw1FwaRWzv4T0twhhBCi+rhVQ5GQkMBDDz1UrHzGjBkkJCSU+Tznz5+nVatW/PLLL87H+++/73LMihUrGDt2LCNHjuSjjz5CKcXAgQNJTU11J/Q6xaUPhbPJo3CERyPPBCWEEOKq5FZC0b59e3bv3l2sfPfu3bRv375c5/Lz86Nly5bOR2RkpMv+uXPnMnXqVGbNmkWfPn347LPPyM/P591333Un9DrlYkLhd3EdDxnhIYQQwgPcSiimTJnCuHHj+O9//8uGDRv4888/+e9//8v48eOZOnUqsbGxzseVbNy4kc6dOzNw4ECeeuoplxk4MzMz2bVrF8OHD3eWGQwGhg4dyrp169wJvU4p2uThrKFw9qGQGgohhBDVx60+FI899hjg6EtxqZkzZ7psK6VKPU9ISAhPPPEEgwYNIiEhgX//+9+sWLGCHTt2YDAYOHPmDEopwsJc+wOEhYWxZ8+eUs9rNpsxm83ObZPJVKbrqm2KNnkE+Vzah0JqKIQQQlQftxKKCxcuVMqHP//88+h0jvUnunfvTrdu3WjWrBmLFy/mzjvvxG63A6DX613eZzAYsNlspZ533rx5PPvss5USY01WNKEI8C74VkofCiGEEB7gVkJRv379SvnwwmSiUEREBE2bNmX//v0un5OSkuJyXHJy8mVjmDNnDrNnz3Zum0ymYn0zajulFOnOqbd98TMUJhRSQyGEEKL6uT2xFUB2djZpaWnFyiMiItw6n9lsJikpieDgYMAxxXdkZCSbN29m1KhRzuM2bdrEiBEjSj2P0WjEaDS6FUNtkWnJxK4ctTR+ukC0Wk3BLJkFo1/8ZdioEEKI6uNWp8zY2Fi6deuGv78/kZGRxR5lkZ+fzz/+8Q9n7UNOTg4PPfQQNpuNCRMmOI+bPn06Cxcu5NChQwB88MEHnDhxgvvvv9+d0OuMjLyCDpl2A/5GH0dhYYdMrR58ZRl5IYQQ1cetGop7772Xjh078u677zprE8pLr9cTERFB586dsdvtpKWl0aVLF9asWUPz5s2dxz3xxBOcOXOGLl264Ofnh06n44svvqBTp05ufW5dUbT/hL/x0v4TMkumEEKI6uVWQnHgwAHWr19PQECA2x+s0WiYOXMmM2fO5Ny5cwQFBeHt7V3sOJ1OxzvvvMMrr7xCamoq4eHheHlVqKWmTjDlO0auKJtPkQ6ZMkumEEIIz3CryaNVq1acOnWq0oJo1KhRiclEUX5+fkRGRkoyUSDHmgOAshvx9y4YBZNVpIZCCCGEqEZu3Z3nz5/P3XffzTPPPEOLFi3QXFK93rZt20oJTpQu15rreGE3FK+hkA6ZQgghqplbCYXNZuPgwYOMHDmyxP2Xm8xKVI4cy8UaioCS+lAIIYQQ1cithGLWrFlMmTKFRx55xO1OmaJiCps8sOuLdMqUPhRCCCE8w62EIikpiddeew1/f//KjkeUkUsNhfShEEII4WFudcps27Ytx44dq+xYRDlc7JRpwF/6UAghhPAwt2oobr/9diZNmsRzzz1Hy5Yti3XK7Nq1a2XEJi6jsIYCu8HRh8KaDzkFU5TLtNtCCCGqmVsJxb/+9S8Axo8fX+J+6ZRZ9QpHeTiaPLwg+7xjh8ySKYQQwgM8utqocJ+zU6YqaPLITHBsyyyZQgghPMCjq40K9+VaCmsoDI5RHplJjh3+DT0YlRBCiKuVW50ywbGE+EcffeRs/gDYsmULNputUgITl1e0U2aAtx6yC2qN/CShEEIIUf3cXm20Xbt2vPbaa7z00kvO8k8++YRPP/200oITpbvYKbOwD0WyY9u/geeCEkIIcdVyK6F4/PHHmT17NgcOHHApnz59Oq+//nqlBCYuL9uSDRTWUHgVqaGQhEIIIUT1c6sPxY4dO1iyZAmAy5DRVq1aceTIkUoJTFxeTsEoD40y4KPXSUIhhBDCo9yqodDpdGRmZhYrP3ToECEhMmSxOhQOG/X18nUkdZJQCCGE8CC3EoqbbrqJuXPnYrfbnTUU8fHxzJgxo9QFw0TlsdgtWOz5AAQY/ByFhX0o/GQEjhBCiOrnVkLx2muvsXnzZsLDw7Hb7XTq1ImWLVuSl5fHvHnzKjtGcQnn0uWAvzOhkBoKIYQQnuNWH4pGjRqxc+dOli5dyo4dO7Db7Tz11FPcfvvtGAyGyo5RXMK5MJjS4Wc0gt12cdptSSiEEEJ4gFsJRdeuXYmNjWXcuHGMGzeuxH2i6lxcurygQ2ZOKqAADfhIHxYhhBDVz60mj927d5dYbrPZ2L9/f4UCEld2cR0PA76GIiM8fENA51aOKIQQQlRIue4+W7ZsKfE1gN1uZ/PmzURFRVVOZKJUziYPuwFvGTIqhBCiBihXQtGvX78SXxcKCwvjjTfeqHhU4rKcnTLtRtcaCkkohBBCeEi5EorCuSfq169PcnKyyz69Xo/RaKy8yESpitZQOCa1kiGjQgghPKtcCYW/vz8AeXl5VRKMKBuXTpkGmXZbCCGE57m92qjwnOI1FJJQCCGE8KwaNSTAarWi1WrRai/mOXa7Hbvd7nKcRqNBp9NVd3g1hnOUhzLgY9DCBWnyEEII4Vk1pobi8ccfR6/XM3v2bJfyhx56CIPBgLe3t/PRo0cPD0VZM1xs8jBKk4cQQogawa2EomvXrm7tK82qVav45ZdfaN26dYn7b7vtNqxWq/NxtU+cJU0eQgghahqPT2yVmJjI/fffzxdffIGPj4874Vx18mwFnWKdE1sVNnlIQiGEEMIzPDqxld1uZ+rUqcycOZNu3bqVetyqVaswGo0EBgYycOBAXnnlFVq0aFGe0OsUs80MgFJe+GkskF+wlLz0oRBCCOEhHp3Y6qWXXsJut/P444+Xekzr1q35/vvvGTRoEAkJCcycOZNBgwaxb98+goODS3yP2WzGbDY7t00mU5ljqg3ybY6ly1FeBNjTHa91BjAGeiwmIYQQVzePTWy1bds2FixY4FyttHAkh1IKq9WKl5cjtKKdNFu1asXXX39NWFgY33zzDQ8++GCJ5543bx7PPvtseS6tVsmzOpo8lF2PvzXdUejXADQazwUlhBDiquaxia127NhBeno6LVu2dJbZbDb27NnD22+/jdlsLnFoaFBQEJGRkRw7dqzUc8+ZM8clETGZTERGRlY45pqiaA2FnzXV8VqaO4QQQniQ28NGk5OT+eijj/jXv/7lLNuyZQs2m61M73/ooYdcRm5YrVa6dOnCo48+itVqLXWeibS0NE6dOkWTJk1KPXdhf4uij7rkYqdMPT75hQmFdMgUQgjhOW4lFLGxsbRr147XXnuNl156yVn+ySef8Omnn1ZacGazmZEjR/LXX3+RkZHB3r17GT9+PEFBQUyZMqXSPqe2Kdop02hOcRRKQiGEEMKD3EooHn/8cWbPns2BAwdcyqdPn87rr7/udjBeXl4uNRNGo5GHH36Yf//73zRt2pTRo0fTuHFjtm/fToMGV+8NNM9a0OFUeWEwS5OHEEIIz3Nr6u0dO3awZMkSwDENdqFWrVpx5MgRt4PZsWNHsbLhw4czfPhwt89ZFxUmFMquxytX5qAQQgjheW7VUOh0OueIj6IOHTpESEhIhYMSl1fY5GHQGdDkSEIhhBDC89xKKG666Sbmzp2L3W531lDEx8czY8YMRo4cWakBiuIKR3l464wy7bYQQogawa2E4rXXXmPz5s2Eh4djt9vp1KkTLVu2JC8vj3nz5lV2jOIS+XZHDYWPl3eRabelD4UQQgjPcasPRaNGjdi5cydLly51Tkz11FNPcfvtt2MwGCo7RlGEUgqLvbCGwnAxofCVhEIIIYTnuJVQABgMBsaNG8e4ceMqMx5xBfkFyQRAPb0G7BbHhq/0XRFCCOE5ZU4oNm7cWOaTDhgwwK1gxJUVTrsN0EBXkEzoDKD39VBEQgghRDkSioEDB5b5pEopt4IRV1bYIVMpDQ10uY5CnxBZx0MIIYRHlTmhyM3Nrco4RBkVDhlFeRGqyXG89qnnuYCEEEIIyjHKw9vb2/no27evy/al+0TVcU67bdcToi2YC0T6TwghhPAwt4aN7t69u8Rym83G/v37KxSQuLyiNRRBZDteSw2FEEIIDyvXKI8tW7aU+BrAbrezefNmoqKiKicyUaKLS5frCVIFNRQ+wR6LRwghhIByJhT9+vUr8XWhsLAw3njjjYpHJUpVuHS5snsRQGFCIU0eQgghPKtcCUXh+h3169cnOTnZZZ9er8doNFZeZKJERWso/JTJ8VqaPIQQQnhYuRIKf39/APLy8q5wpKgqzk6Zygtfm3TKFEIIUTO4NVPm6tWrL7t/6NChbgUjrsw5sZVdj681xfFaaiiEEEJ4mFsJxbBhwy67Xya2qjoXmzy8MFoLmzykhkIIIYRnuTVsVCnl8rDZbBw6dIhBgwbx5ZdfVnaMooiLTR56jJZ0R6HUUAghhPAwtxKKYifRamnTpg0LFy6U5curmHMeCrsOfX6G47X0oRBCCOFhlZJQFGrQoAFxcXGVeUpxicKEQq80aJXNUSg1FEIIITzMrT4UZ86cKVaWlpbGK6+8Qvv27SsclChdYUJhKOyn4uUDeh8PRiSEEEK4mVBERkaWWN66dWu+/vrrCgUkLq8wofBWdkeB1E4IIYSoAdxKKA4ePFisrF69ejRq1KjCAYnLKxzl4SMJhRBCiBrErYSibdu2lR2HKKPCeSh8CvtPeAd5MBohhBDCodydMpVSfP3119xyyy20adOG1q1bM3LkSL7++muZf6IaFDZ5OGsoJKEQQghRA5QrobDZbIwbN47JkyeTnJxM37596devHxcuXGDy5MmMHz8em81WVbEKINfqSCj8lMVR4B3owWiEEEIIh3I1ebzxxhts3bqVzZs307dvX5d9mzdvZvz48bz55ps89thj5Q4kJyeHEydOUL9+fcLCwortT09PJzk5maioKAwGQ7nPX1fkWhxNHn7K6igwSkIhhBDC88pVQ7Fo0SLef//9YskEOJYzf++991i0aJFbgdxzzz106tSJ//znPy7lVquVadOm0ahRI6699loaNmzIZ5995tZn1AV5BTUUgRQkFNLkIYQQogYoV0Jx9OhRrr/++lL3X3/99Rw5cqTcQXzwwQckJSXRqVOnYvteeuklVqxYwYEDBzh79ixvvvkm99xzD7GxseX+nLogr6APRRDS5CGEEKLmKFeTh9FoJDMzE29v7xL3X25fafbv38/cuXPZsmULo0aNKrb/gw8+4P7776dFixYA3HnnncybN4+FCxfy1ltvleuz6oLCGoogChYJkyYPIS5LKYXVapX+XUKUQqfT4eXlhUajqdB5ypVQ9O7dmw8//JCnnnqqxP0ffvghvXv3LvP5cnNzmTBhAq+++ipRUVHF9iclJZGQkECfPn1cyvv168fOnTvLE3qdkV9QQxFMwZoeUkMhRKny8/NJTEwkJyfH06EIUaP5+voSHh5eoT6K5Uoo5syZw/Dhwzl//jyzZs0iOjoagPj4eBYsWMDbb7/Nb7/9VubzzZo1i65duzJ58uQS96ekpAAQGhrqUh4aGkpycnKp5zWbzZjNZue2yWQqc0w1Xb7dUTMRpAoTCulDIURJ7HY7cXFx6HQ6GjdujMFgqPBfYELUNUop8vPzuXDhAnFxcbRq1Qqt1r1lvsqVUFx//fUsXLiQv/3tb7zxxhvO5o28vDz8/Pz46KOPuO6668p0rl9//ZVvv/2WFStWsG/fPud5UlJS2LdvHx07dsTLyxFefn6+y3vNZjN6vb7Uc8+bN49nn322PJdWazhrKJRjtAdGSSiEKEl+fj52u53IyEh8fX09HY4QNZaPjw96vZ74+Hjy8/PL3XWhULlnyrzrrru45ZZbWLZsmbMDZuvWrRk1alSxmoTLSU1NpUmTJjz44IPOsvj4eJKTk4mJiWH37t00adIEjUZDYmKiy3sTExNLXU8EHDUps2fPdm6bTKbLHl+b5NsLEgp7QRWuNHkIcVnu/rUlxNWkMn5P3Jp6OzQ0lHvuuadCHzxp0iQmTZrkUta1a1cGDx7MggULAPD396dXr16sXLnSeazZbGb16tU88cQTpZ7baDRiNBorFF9NpJTCUpBQ1HMmFFJDIYQQwvNqfOr+7LPP8vXXXzN//nz+/PNPJk6ciJ+fH9OnT/d0aNXOareicEy57Vu4loeM8hBCVJKnnnqqwn8slsWLL75Yat85UXvVqISiZcuWhIeHu5SNGDGCFStWsH79eh577DGCgoLYuHEjwcHBngnSg/Jsec7X3nYFWi/Q+3gwIiFEXWIymUhLS6vyz8nMzCQ1NbXKP0dUL7eaPKrK999/X2L5iBEjGDFiRDVHU/MUrjSKAj04mjuk17oQQogaoEbVUIjLK6yh0CgvNCDNHULUQb/88gstWrRg2bJl9O3blyZNmnD48GEAVq1axZAhQ2jatCkDBw4stgzBwYMHGTt2LK1ataJv3768+eab2O125/6PP/6Y3r17Ex0dzdChQ1m/fn2pcfz66680bdqU7Oxsl/L//Oc/LjMmXykmgNdff50OHTrQuXNnZs6cWeycom6oUTUU4vLMBbNk6uw6R4GM8BCiXJRS5Fqqf8ZMH72uzHNg5OXlceLECZ5//nneeOMNmjdvTv369fnxxx958MEHefPNN+nTpw8HDx7k/vvvJz8/n/vuuw+AkSNHMmzYMH7++WeysrL48ssvWbVqFTfffDOffvopjz76KO+88w59+vThq6++YtiwYcTExNChQ4dicQwePJiMjAyWLFnClClTnOUfffSRs59FWWL68MMPmTt3Lh988AFdunThvffe4+233+aGG26o6JdV1DAapZTydBBVzWQyERQUREZGBoGBtfcmvC95H5N+noTR4s2OM0eg2bVw13JPhyVEjZSXl0dcXBzNmjVzjqvPybfS/ulfqz2WA88Nx9dQtr/flixZwq233sqePXtc1jdq06YNs2bNYsaMGc6yd999l3feeYe9e/eSl5eHj48P27Zto1evXs5jlFJoNBpatGjBlClTeO6555z7hgwZQnh4OF988QUADz/8MGfOnGHJkiUA3HfffZw9e5aVK1cCsHXrVvr160dcXBzR0dFXjAmgefPmTJs2jX/961/OeDp27EhkZCS//PJLeb6MogqV9PtSqKz3UKmhqEUK+1DoVUFLlTR5CFEn6fV6Onbs6NxOTU3lyJEjPPPMM7z44osopVBKkZeXR25uLgDe3t6MHj2a8ePHc++993LdddfRr18/dDod2dnZnDhxgoEDB7p8zqBBg0rtuwYwZcoUbrjhBs6fP0/Dhg358ssvGTBgANHR0WWKKSsri7i4OK655hrnOTUaDddccw2nTp2qzC+ZqAEkoahFzAWzZOrtBVWnMgeFEOXio9dx4LnhHvnc8rh0mvDCpQTeeustrr32Wpdjix73448/smLFCn799VfuvvturFYrS5YscS6TcOk6DQaDodhMxEUNHjyYsLAwvv32W2bMmME333zD888/X+aYLBbHqsiXzmxckfUiRM0lCUUtUtgp01DYSiUJhRDlotFoytz0UJOEhYVRv3599uzZw/jx40s9TqvVMmrUKEaNGoXdbufGG2/k5Zdf5quvvqJBgwbExMQwaNAg5/G7du2idevWpZ5Po9EwadIkvvzyS1q2bEl6ejrjxo0rc0z16tWjfv367N6926WWIjY2Fn9///J+GUQNJ6M8apHCJg/vwoRCmjyEuCpoNBr++c9/8vrrr/PTTz9hs9nIz8/n999/d9YYJCQkMGvWLGdTQkZGBhcuXCAsLAyARx99lFdffdW5dtLSpUv56aefmDlz5mU/e+rUqWzZsoUXXniBm2++mXr16pU5JoAZM2bw8ssvc/jwYZRSLFq0iE2bNlX610h4Xu1L1a9ihU0ezoRCRnkIcdX4+9//jpeXFw8//DCTJ0/Gy8uLAQMG8NJLLwHQqFEjmjVrxrXXXktaWhoWi4WRI0cyd+5cAJ544gnOnz9P7969HTU1vr4sWLCAIUOGXPZzO3fuTKdOndi0aRM//PBDuWICx+ybJ06coEOHDvj5+dGlSxcmTJhAenp6pX59hOfJKI9a5KuDXzFv2zy6Zhr4PPkYjH4buk31dFhC1EiX67Vek5nNZjIyMmjYsGGpx6SnpxMQEIBOV3LfjMzMTPz8/Epc8MlqtWIymahXr16xoayZmZnYbLZiMxFnZGSQm5tLgwYNSv3MK8WUk5ODl5cXBoOBrKwsLBaLs7ZDeJ6M8rjKFPah8FNWR4E0eQhR5xiNxssmE8AVlx4ICAgodZ+XlxchISHlel9QUBBBQZfvs3WlmIouIS/9J+om6UNRixRObOVXuDCYNHkIIYSoISShqEUKayj8lWMoltRQCCGEqCkkoahFCkd5BBYmFDJsVAghRA0hCUUtUphQBEhCIYQQooaRhKIWySlIKIwyD4UQQogaRhKKWiQ73zE/vo9SKC8f8JLpa4UQQtQMklDUIjmWizUUGhnhIYQQogaRhKIWyS3a5CHNHUIIIWoQSShqkVxrQZOH3S4dMoWoow4cOMDdd9/t6TAq7Pjx49x9993OFUdF3ScJRS2SVzCxlVEpmdRKiDrq7NmzfPrpp54Oo8IuXLjAp59+is1m83QooppIQlGL5BddHEyaPIQQQtQgspZHLWK2F1m+XJo8hLgq/PnnnyxbtoxXX33VWXb06FFefPFFPvzwQ/R6Pbt27eKDDz7giSee4LvvviMuLo5u3boxbdo04uPj+fTTT0lJSeG6667j9ttvd56n8H3//ve/+eKLL4iPj6dnz57cfffdLot8Wa1Wvv32WzZt2oS/vz/Dhg1j6NChLnGmpqby7rvvkpiYSMeOHWnXrl3Vf3FEjSI1FLWIxS5NHkJcbY4cOcL333/vUnbu3DmX5oRTp07x0UcfMWzYMPLz82natClPPPEEt956KyNGjMBoNBIeHs7dd9/NRx995DxP4fv69+9PZmYmTZs25d///jf33Xef8xiLxcLw4cN58803adWqFcHBwdx1110899xzzmNycnLo27cvK1asoFWrVmzcuJGJEydW8VdG1DRSQ1GLWO35AHjbFRilhkKIclMKLDnV/7l6X7hkqfDKZrVaWbRoEQMHDgQgPz+fp59+mm3bttGrVy/AsQz1p59+yrRp01ze99xzzzk7gg4ZMoTevXszc+ZMunbtynvvvUdSUhIxMTEYDI65b4YOHUr//v2ZMWMGDRo04N133yUnJ4c9e/Y4l76eOHEi33zzTZVes6hZJKGoJWx2GzYcvaWlyUMIN1ly4KXG1f+5T50Fg1+VfoRer+eaa65xbrdo0QI/Pz9nMlFYdmltB+DSDNKzZ0+aNm3Khg0b6Nq1Kz///DN2u52HHnoIpRRKKex2O1arlf379zN48GDWr1/PzTff7EwmAMaPHy8JxVXG4wlFXl4eMTExZGVl0a5dOyIiIlz2x8TEcPz4cZeyoKAghg0bVp1hepy5oEMmSJOHEKI4g8GAVnuxFVur1brc4AvLLh11odfrCQgIcCkLCQnhwoULAKSkpBAVFcWAAQNcjhk8eDAtW7YEHCM6OnToUOwc4uri0YTik08+4dlnnyUqKgovLy82b97M1KlTee+995y/GO+//z4rVqygb9++zvdFRUVd1QmFjPIQwk16X0dtgSc+101Go5H8/HyXsrS0tIpG5GSxWEhMTCQ8PBwApRSnT58mMjISgIiICHJyci47N0ZUVBSnTp1yKYuPj6+0GEXt4NGEws/Pj9jYWIKCHNX3sbGxdOvWjVGjRnHLLbc4j+vbt2+J1XRXk8KEQm9Xjp600uQhRPlpNFXe9FDZWrVqRVJSEseOHaNly5bY7XYWLVpUqZ/x+uuv8/LLLwPw+eefYzKZuPnmmwG48847uf3221m1ahU33ngj4Eg6vvnmG2fHy9tuu41p06YRFxdHs2bNyMvL45133qnUGEXN59GEYty4cS7brVu3xsvLi4yMDJfyjIwMVq5cSVBQEJ06dSIw8Or767xwlkxjwUKj0uQhxNWhb9++jBw5kmuuuYZBgwaxf/9+mjRpUmnnNxqNbNq0iQEDBhAQEMDatWt55ZVXaNzY0dfk1ltvZe7cudx666306dOH0NBQ9uzZQ+/evZ0JRWF/iR49ejBo0CD27dtHVFRUpcUoageP96E4d+4cGzZswGQy8dlnn3HDDTcwduxYl2N27tzJggULSEhI4OzZs7z55pvccccdpZ7TbDZjNl9sIjCZTFUWf3UxOye1sjsKfEM9GI0Qoqp06NChWA3EkiVL2LhxI+fPn6dLly4EBgaycuVK9Ho9AN27d+f99993eU+fPn146623XMoGDRpEcHCwS5mXlxfr169n+/btxMfH88Ybb9C6dWuXY55++mkeeOABNm/ejEajoVu3bkRHRzv3azQafvzxRzZu3EhSUhIdOnSgfv36/Pzzz84YRd2nUUqpKx9WdWJjY3nhhRdITk7mwIEDzJ49myeeeAJNwRCrtWvX0rdvX3x9HW2Qr776Kk899RS7du2iY8eOJZ5z7ty5PPvss8XKMzIyam3tRuz5WO5YdQcRFgurziTCvy/I8uVCXEZeXp6zCv7SzonCYcmSJUydOpWsrCxPhyI87HK/LyaTiaCgoCveQz0+sVXXrl35/vvvWbduHWvXruX555/n448/du6//vrrnckEwOOPP05gYCCrVq0q9Zxz5swhIyPD+Th9+nSVXkN1SDenAxBkt2M1BEoyIYQQokbxeEJRVMeOHenZsyfr168v9RiNRkNAQADnz58v9Rij0UhgYKDLo7ZLy3P06q5ns2P3re/haIQQdUFJTSVCuMtjCYXNZiM1NdWlLDs7m6NHjzrnorDb7SQnJ7scExMT45xv/mqSmuf4WoXYbCAJhRCiEkRFRTFlyhRPhyHqCI91yrRYLFxzzTWMHj2a9u3bk56ezieffIK3tzePPvoo4Eg6BgwYwKhRo+jQoQOnTp3ijTfeYPjw4cU6btZ1ybmOhKKe3Y7Gv4GHoxFCCCFceayGwtvbm23bttG4cWP++OMPjhw5wt/+9jcOHDhAWFgY4JjBbfv27URERLBu3TrS0tL48MMPWbVqlctKeFeD5BxHQhFss6OThEIIIUQN49FhowEBAc7aiIocczVIyb3Y5KH1b+jhaIQQQghXNapTpihdWl464GjywE9qKIQQQtQsklDUEun5haM8bOAnk1oJIYSoWSShqCVM+emAY9io1FAIUXedO3eOxYsXezqMCktJSWHx4sXY7fYKn6ssX5OEhAR++OGHCn9WWe3atYvvv/+erVu3Ao4ZmteuXcs333zDmTNnqi2OmkQSilrAbDOTb3es5VHPbpOEQog6bO/evUyaNMnTYVTY0aNHmTRpUrGVUt1x6dckKSmJb7/91uWY7du3c9ddd5XrvGlpaSxevLjY+lFXMmvWLEaOHMm3337Lrl27MJvN9OjRg3/84x/89NNPJCQklOt8dYXH1/IQV1Y4qZWXUgTYlSQUQoirSlhYGBMmTHBux8bGMnnyZMaPH1+h8x4/fpxJkyaxd+9e56rXZfHhhx/y/fffO1dfXbt2LcePHycjIwOD4eqdxVgSilrg4iyZNkADPvU8G5AQotqcOHGC/fv3M3LkSGdZcnIyq1evZvz48Wi1WhISEtiyZQu33nor+/btIy4ujs6dO9OsWTPsdjtbt24lJSWFnj17OoflAy7v27NnD/Hx8XTt2tVl4a9CSUlJbN++HX9/f7p3717iDXjnzp0kJibSoUOHK17XDz/8QN++fZ0rp27ZsoWzZ89y2223AY4l0r/99luGDBlCgwYNGDNmDOCoVVi/fj1KKWczSNHFzJRS7Nu3j/j4eNq3b0/z5s3L8FUu27WeP3+eVatWkZOTw86dO501G5s3b0av1/Pjjz8CMHHiROfX9rbbbmPr1q2cPn2aG264wXmuvXv3cvz4cSIjI+natWuxqRDi4+PZt28foaGhdO/e3SVRsdlsbN26lXPnztGqVati61pt3boVpRSdO3cmJiaGrKws+vTpU2xhuMomCUUtUHTabat3CHrt1TUHhxBXs7Vr1/LCCy+4JBSHDh1i0qRJjBkzBm9vb7Zv386UKVPo1q0bOp0OvV7Phg0b+O9//8tXX32FwWBAKcXu3bv57bff6Nu3L4Dzff3793cu/PTXX3/x1ltvcf/99zs/b968ecyfP5++ffuSm5vLgQMH+OKLLxg+fLjzmKlTp7J8+XL69+/PwYMHad++/WWv69VXX2XkyJE89dRTAMyYMYPY2FgSEhJo3Lgxe/bsYdKkSVy4cMHZ5DFx4kTS0tLYvHkzSimWLFkCwPDhw6lXrx5Wq5URI0ZgMpkIDAxk3bp1xa7lSi53rcnJySxfvhyA9evXs2/fPgBOnTqF2Wx2xjNx4kTn13bgwIFkZGTQvHlz+vXrh1arZfz48ezbt49u3bpx5MgRAgICWL58uTPZe+GFF3j55ZcZMGAAWVlZmEwmfvjhB1q0aMG5c+cYMWIEKSkpdOjQgS1btjBs2DC+/vprZ1Ly7rvvsnv3bvLy8oiOjiYpKYnExETWr19P27Zty/y1KDd1FcjIyFCAysjI8HQobll6dLnq+ElHNe3dlsr89gBPhyNErZCbm6sOHDigcnNznWV2u11l52dX+8Nut5c57t9//10V/a/5ww8/VNHR0S7HbNiwQQHOa/vpp58UoP73v/85j3nwwQcVoD788ENn2ZQpU9TIkSOd24Xvmzlzpsvn+fr6qqSkJKWUUitXrlT169dXcXFxzmM++eQT1bBhQ5Wdna2UUuqHH35QPj4+6vDhw0oppbKzs1X37t1dYrzUk08+qYYNG6aUUiotLU3p9XrVrVs39eWXXyqllFqwYIHq3LlziV+TVatWKZ1O53K+wmtZsGCBs2zBggUqJCSkxM9XSqnt27crQO3du7fM15qZmakAtXnzZucxX3/9tQoNDS0xnueff96l/IEHHlA33nijMpvNSimlbDabuu2229SUKVOUUkrl5+crg8GgfvvtN+d7jhw54ozx7rvvVj169FBZWVlKKaVOnDihAgMD1fvvv+88/q677lI+Pj7qwIEDSinHz/3111+v7r333lK/FiX9vhQq6z1Uaihqgfj0c4CjycMr6upaw0SIypRrzaXPV32q/XO3Tt6Kr973ygdWgFar5b777nNu9+vXj48//ph77rnHpex///tfsff+4x//cL6+9957eeqpp1ixYgXTpk1j0aJFtG/fnh07drB9+3aUUuh0Os6fP8/+/fvp1asX3377LaNHj3Y2Pfj6+vLwww9z7733lhrv4MGDefPNN8nPz2f9+vV07tyZkSNH8scffzB58mTWrVvH4MGDy/01ePDBB10+Y9asWSQnJ1O//pXXQCrLtZZX0YkZrVYrX3zxBQ899BDLli1DKYVSioiICL7//nvnNRiNRvbu3cuQIUPQarW0atUKuNgM9N577+Hn5wdAs2bNmDx5MosXL+aBBx5wufZ27doBjkU1r732Wn7//fdyx18eklDUAqdNFwBHk4c2QhIKIURxPj4+eHt7O7cLV10u2jZvNBrJy8tzeZ9er3f2YwDHDS0qKor4+HgATp48SVZWlvOGV2jChAl4eTluIadOnWLQoEEu+5s1a3bZeAcMGIDFYmHbtm3O5GHw4MFMmzYNpRQbNmzgww8/LMdXoOSvAVDsmktTlmstD19fX5fVrs+fP09OTg6xsbGcPn3a5djC5Emn0/H5558ze/Zs5s2bx6BBg5g4cSJjx47lwoUL5OTkFOsX0qJFC3777TeXspCQEJftkr73lU0Silpgf3IMAM0tFpCEQgi3+Xj5sHXyVo98rru0Wm2xuRwq88ZgsVjIysrC39/fWZaWlub8iz4wMJCWLVvy1VdflXqO0NBQ0tLSXMou3b6Un58fPXr0YN26daxbt44XXniBvn37kpCQwMqVK0lNTeXaa6+twJWVX1mutTw0Go3LdkBAABqNhvvvv/+yI1RGjx7N6NGjOXr0KKtWreK+++7j1KlTPPLII+h0umIrdaemppapBqaqyTwUNVxqXiqncw8C0DdPC6GtPByRELWXRqPBV+9b7Y9Lbyzl0aRJE86dO0dmZqaz7I8//qiML4dTYWdCcAzJjIuL45prrgFgxIgRLF++vNjcCkW3BwwYwMqVK13mnCgc8XA5gwcP5qeffmL//v0MHDgQo9FI3759mTt3Lp06dSI0tORZgf39/bHZbFgslvJc5hWV5VorIiAggP79+/P++++jlCrxM3Jzc53JWKtWrXj00UcZPXo0W7ZsQa/X07t3b5evrc1m46effmLAgAGVEmNFSA1FDffHqfUoFG3N+fgEtQet5IBCXE0GDRpEeHg4Y8eOZcKECezevbtYlXxF6HQ6nnzySU6ePElAQACvvvoq48ePp0ePHgA8/PDDLF++nD59+vC3v/2N0NBQYmJiWLNmDUeOHAHgb3/7G++++y7Dhg1j8uTJ7Nixg19++eWKnz148GD+85//0KtXL2fTwODBg5k7d+5lF4Vs3749vr6+PPHEE/Tu3dtl2GhFlOVaK+qdd95hyJAhXH/99YwbN468vDzWrFlD8+bN+d///kdmZib9+/dn1KhRdOrUiTNnzvDDDz/w8ccfA/Daa69x/fXXo9Pp6N27N99++y2ZmZnMmTOnUuKrCLk71WA2u43PYxyzwQ3OySW0wxAPRySEqGqXTuLk7e3N5s2b6dmzJ5s3b6Z169asWrWKCRMmOPtHREREMHbsWJfzREdHO+d0KNSiRQuX4aeF59+wYQMWi4WYmBieeOIJPv/8c5f9a9asYf78+Zw8eZJdu3bRrVs3du/e7TzG39+frVu3MnDgQLZv30779u1Zs2aNS4wlGTBgABMnTmTGjBnOstGjRzNhwgTGjRtX6tckJCSE1atXY7VaWbp0Kbt37y7xaxAYGMiECRPw9S25Q2xhDUdhv4uyXKter2fChAkuTQwlfa1Ligegc+fO7N+/nxtvvJFt27Zx9uxZHn30UWdn2YYNG7Jjxw6ioqLYuHEjGRkZ/Pbbb84mkn79+rFz505CQkLYuHEjgwcPJiYmxiWePn360K9fP5fPbd++vcsw36qgUZfWu9RBJpOJoKAg5zjritoYs5yfdr4FgKLwy1f0y6gK/jnKNAXHFX1GgdIoUJe8r+DZiuK4NpfT+nw0SvFxqp6ef9sARn+EEFeWl5dHXFwczZo1c+moJy5asmQJU6dOJSsry9OheMSXX37JnXfeSXZ29lX/M3K535ey3kOlycMNx87G8pv+bLV8lq/dzt9T8+k05XtJJoQQohKcPXuWZcuW8fLLLzN58uSrPpmoLJJQuKFlk24MT/qrYEvjfNagQQGO/leaIvuLdMjSXDyu8F3O4zQaVEGJBi1hXv5c23wQPfqOB738wAshKldp1fJ13blz59ixYwd///vfyzWLprg8afIQQtRJ0uQhRNlVRpOHdMoUQgghRIVJQiGEEEKICpOEQghRp10FrbpCVFhl/J5IQiGEqJP0ej0AOTk5Ho5EiJqv8Pek8PfGHTLKQwhRJ+l0OoKDgzl//jzgWKipIlNgC1EXKaXIycnh/PnzBAcHX3YisiuRhEIIUWeFhYUBOJMKIUTJgoODnb8v7pKEQghRZ2k0GsLDw2nYsGGlLyQlRF2h1+srVDNRSBIKIUSdp9PpKuU/TCFE6aRTphBCCCEqTBIKIYQQQlTYVdHkUTi+1mQyeTgSIYQQonYpvHdeaa6KqyKhyMzMBCAyMtLDkQghhBC1U2ZmJkFBQaXuvyoWB7Pb7Zw9e5aAgIBKG4duMpmIjIzk9OnTdWbBsbp2TXXtekCuqbaQa6od6to1VdX1KKXIzMykcePGaLWl95S4KmootFotERERVXLuwMDAOvGDWFRdu6a6dj0g11RbyDXVDnXtmqriei5XM1FIOmUKIYQQosIkoRBCCCFEhUlC4Saj0cgzzzyD0Wj0dCiVpq5dU127HpBrqi3kmmqHunZNnr6eq6JTphBCCCGqltRQCCGEEKLCJKEQQgghRIVJQiGEEEKICpOEwg2ZmZns2LGDkydPejoUt+Tk5BAbG0tCQkKpxyil2L9/P7t378Zms1VjdBVz4MABNm7cSH5+frF9FouFmJgYDh486IHI3JOQkEBsbGyJ1wOQnZ3Nzp07iYuLq+bI3JOYmMjOnTuJj48v9ZiUlBS2b99OUlJSNUZWNkopdu3aRWxsbKnHWCwWYmNjOXjwYKlTFZflmOpitVrZsmULhw8fLvWY8+fPExsbS0ZGRqnHFP4snjhxoirCLJfs7Gw2bdrEqVOnLnucUoq//vqr1O9nWloa27dv5+zZs1UQZfkkJyezceNGkpOTSz1GKcXBgwcv+71MSkpi+/btpKSkVH6QSpTLhx9+qHx9fVWbNm2Un5+fuvHGG1VWVpanwyqTxMREddddd6mgoCDVtWtXVa9ePXXNNdeouLg4l+MOHjyoWrdurRo1aqQiIiJURESE2rp1q2eCLoddu3YpHx8fBajTp0+77Fu7dq1q1KiRatq0qQoNDVWdO3dWJ0+e9FCkV5aQkKCGDBmigoKCVM+ePVWzZs3UypUrXY757LPPlL+/v2rdurXy9/dXQ4cOVRkZGR6K+PJSU1PVkCFDVEBAgOrRo4cKCQlR3bt3L/az9/TTTyuj0ajat2+vjEajmjZtmrLZbJ4J+hL//e9/VevWrVVwcLDq0qVLicesX79ehYWFqejoaFW/fn3VsWNHdeLEiXIfUx2ysrLU//3f/6nIyEgVEBCgJkyYUOyYjRs3qv79+6uGDRuqrl27Kh8fH/Xggw8qq9XqctyXX36pAgICVOvWrVVAQIC67rrrVHp6enVditPZs2fVI488osLDw5XRaFTPPPPMZY+fP3++0mq1JX4/X3zxRefPore3t7rjjjuUxWKpmsAv48CBA+qOO+5Q4eHhClBff/11icetWbNGNW3aVEVGRqquXbuqvn37qvj4eOd+q9Wqpk2bpry9vZ2/X08//XSlxioJRTns2bNHabVa9dVXXymllDp//rxq2rSpeuSRRzwcWdls2bJFffrpp85fiszMTHXttdeqgQMHOo+x2+2qc+fOasyYMc7/yKdNm6YiIyNVXl6eR+Iui8zMTNWmTRv1j3/8o1hCkZGRoUJDQ9U///lPpZRS+fn56rrrrnO57pokPz9fde7cWd1www0qMzNTKaVUSkqKy38khw4dUl5eXurjjz9WSjlu2K1atVL333+/R2K+ktmzZ6vIyEiVkpKilFIqJydHde/eXd1+++3OY5YsWaL0er3atGmTUspxjUFBQeqNN97wSMxFWa1WNWvWLHXo0CH1+OOPl3gDyszMVA0aNFCzZ89WSillsVjU0KFDVb9+/cp1THWJi4tTc+fOVQkJCermm28uMaFYtGiR+uuvv5zbBw4cUMHBwWr+/PnOsqNHjyq9Xq8++OADpZRSaWlpqk2bNuqee+6p+ou4xKZNm9Qbb7yh0tLSVIsWLS6bUGzdulVFR0erqVOnFvt+/vLLL0qn06k//vhDKaXU8ePHVUhIiMt1V5cffvhBffrppyozM7PUhGLfvn3KaDSql156yVm2c+dOtXnzZuf2ggULVHBwsDp8+LBSSqkNGzYoLy8vtXTp0kqLVRKKcpg9e7Zq2bKlS9l//vMfFRQUVCxjry0WLlyovLy8nMnDtm3bFKB27NjhPObkyZMKUMuXL/dUmFd0xx13qFmzZqnff/+9WELx2WefKb1er9LS0pxlv/zyiwLU0aNHPRDt5X355ZdKo9FctgblqaeeUhERES5lCxYsUL6+vjUy8bvjjjvU0KFDXcqmTZumrr32Wuf2qFGj1IgRI1yOue+++0qtDfCU0hKKr776Sul0OpWcnOwsW716tQLUwYMHy3yMJ5SWUJRkzJgxauTIkc7tp59+WoWHhyu73e4se+utt5S3t7fKycmp9FjL6nIJRXp6umrRooX67bff1MyZM4t9P8ePH68GDx7sUvbwww+rNm3aVFG0V2axWEpNKCZNmnTF35POnTur6dOnu5QNHTpUjR49utJilD4U5RATE0OPHj1cynr37k1GRkaNaDd0x/bt22nWrJlzwZeYmBi0Wi3dunVzHhMdHU3Dhg2JiYnxVJiX9cUXXxAbG8t//vOfEvfHxMTQvHlzgoODnWW9e/d27qtp1qxZQ+fOnYmKimL//v0cOnQIi8XickxpP4s5OTkcOXKkOsMtk9mzZ7N//35eeukl1qxZw//+9z9WrFjB008/7TymtGvat29fseuviWJiYmjatCmhoaHOskt/zspyTE1msVjYvXs3LVu2dJbFxMTQvXt3l4UXe/fuTV5eHocOHfJEmFf0wAMPMGrUKIYNG1bi/tJ+Fo8cOUJOTk51hFgua9as4ZZbbiE3N5edO3dy+vRpl/0Wi4X9+/eXeE2V+XN3VSwOVllSU1Np166dS1nhfwypqameCKlC1qxZw8KFC1m0aJGzLDU1leDg4GIryoWGhtbIazx69CiPPfYYf/zxR6mzw6Wmprr8Bw44r7EmXtPZs2fx8/Ojf//+mEwmcnJyMJvNLFy4kJtuuglwXFOHDh1c3leTfxbbt2/PHXfcwcsvv0yLFi2Ii4tjzJgxzpsplPx9Cg0NxWazYTKZiu2raUqKPyAgAL1e7/yelOWYmmzOnDmkpKTw6KOPOstSU1Np0aKFy3E1+Wfxww8/5NChQ3z22WelHlPaz6JSirS0NHx9fas6zDKz2WycP3+e06dP06ZNG0JDQ4mPj6dt27YsXryYqKgoMjIysNlsJV5TZX6PpIaiHPR6PXl5eS5lubm5ABgMBk+E5Lbt27dz22238fe//5077rjDWV7SNYLjOmviNU6bNo2bbrqJ9PR0Nm7cyL59+wDH9RWOfCjpmvLz87Hb7TXymvR6PX/99RePPvoo+/fvJy4ujjvvvJNJkyaRlpbmPKY2/SzOnDmTJUuWcOzYMedfUMePH2fSpEnOY2rbNV2qpPitVitWq9UZf1mOqalee+013n77bb7//nuaNm3qLK9N37ekpCRmzZrFjBkz2L59Oxs3buTs2bNkZ2ezceNG5yiW2nRNOp0OrVbLzz//zJ9//klMTAynTp3Cbrfz4IMPAo7rAUq8psq8HkkoyiE6OrrYUMvC7aioKE+E5JYdO3Zwww03cP/99xdrJoiOjiYnJ4f09HRnmdVq5fz58zXyGhs1asTx48d58sknefLJJ/noo48AmDdvHsuXLwdq3/etadOmBAYGutxsH3zwQUwmE7t37wZq3zWtWLGCiRMnUr9+fQD8/Py4++67WbVqlbM5o7RrCg4OJiAgoNpjLq/o6GjOnj3rMgy0cLvwe1KWY2qiBQsW8O9//5uffvqpWDNBbfpZzM3NpVu3bnzxxRfO/zM2bdpEYmIiTz75JEePHgVKvyZfX98aWVMWHR3N8OHDnYmev78/U6ZMYcOGDYBj6fHg4OASr6kyv0eSUJTDsGHD+PPPP13GYi9dupRu3brVyB+ykuzcuZNhw4Zx77338uqrrxbbP3jwYPR6PcuWLXOW/f777+Tk5DB06NDqDLVMvvvuOzZu3Oh8vP766wD8+OOPzmrZYcOGce7cObZt2+Z839KlS/H396dfv34eiftyhg8fXiypO3PmDAANGjQAHNf0119/uYxJX7p0Ke3ataNJkybVGm9ZNGjQwHkNhU6fPk1wcLDzr6dhw4axcuVKl3lPli5dWmo7d00zbNgwkpOT2bx5s7Ns6dKl+Pr6cs0115T5mJrmzTffZM6cOfz444+MGDGi2P5hw4axdetWzp8/7yxbunQprVq1Ijo6ujpDvaJmzZq5/H+xceNGxo0bR8uWLdm4cSM9e/YEHNf0yy+/uPTdWbp0KUOGDCnWHFwTDB8+vFiycObMGef/FwBDhw51/pEFjj8Uf/7558r9/aq07p1XgdzcXNW+fXt17bXXqiVLlqinn35a6XQ6tWrVKk+HViYHDx5UISEh6rrrrlMbNmxweeTn5zuPmzNnjqpXr5766KOP1BdffKGaNGmi7r33Xg9GXnYljfJQSqlbb71VtWzZUi1evFi98847ys/PT7388sseivLy7Ha7GjJkiBo0aJBatmyZ+vrrr1Xr1q3VLbfc4jzGbDarLl26qP79+6uffvpJPffcc0qn06klS5Z4MPLSLVq0SOl0OjV37lz122+/qf/+97/K39/fpRf+2bNnVcOGDdXYsWPVsmXL1IMPPqh8fX3V3r17PRd4ETExMWrDhg1q4sSJqmXLls7fnaIjvMaNG6eaN2+uvv76a/Xee+8pf39/l6F8ZT2mumzatElt2LBB9e/fXw0ZMkRt2LDBZc6ZhQsXKkD961//cvn/IjY21nmMxWJR3bt3V3379lU//vijevHFF5VOp1Pff/99tV9PXl6eM8YmTZqoe+65R23YsEHt2bOn1PeUNMrjwoULqnHjxmr06NFq2bJl6uGHH1be3t5q586dVXwFxaWlpakNGzaodevWKUDNnTtXbdiwQR07dsx5THx8vKpfv76aNWuW+vXXX9Wrr76qfHx81Lvvvus8Zs+ePcrX11dNnz5dLVu2TN1+++2qYcOGKjExsdJildVGyyk5OZn58+cTExNDaGgo06dP57rrrvN0WGWyatUqXnzxxRL3LVu2jJCQEMAx29pHH33EkiVLsFqtjBgxgocffhgvr5rfh3fHjh3MmjWLn376ySU7N5vNLFiwgLVr12I0GpkwYQJTpkzxYKSXl5uby+uvv8769evx8/Nj8ODBTJ8+3aW9My0tjfnz57Njxw5CQkK4//77a/Rf87///jtffvklCQkJNGzYkDFjxjB27FiX0QFxcXG8/PLLHD58mKioKB577DG6dOniwagvuueee5xV4kX98ssv+Pv7A46+OW+88QarV6/GaDQybtw4lz5KZT2mugwZMgSz2exSFhoaytKlSwH4v//7P/74449i72vbti0LFy50bqenpzN//ny2b99OvXr1uO+++xg+fHjVBl+CpKQkxo4dW6y8W7du/O9//yvxPW+++SYxMTEundMBTp06xfz58zl06BBNmjRh5syZxUZJVIetW7fy+OOPFysfM2YMf//7353bJ06c4NVXX+Xw4cM0btyYqVOnFvse7N69m9dff51Tp07Rpk0b/vnPf9KsWbNKi1USCiGEEEJUWM1rDBJCCCFErSMJhRBCCCEqTBIKIYQQQlSYJBRCCCGEqDBJKIQQQghRYZJQCCGEEKLCJKEQQgghRIXV/JmKhBDVKikpiT///JPx48dX+2efOHHCOUV6hw4d6NSpk9vncuc6kpOTWb16NQCRkZE1djpsIWoiSSiEuEqkpKTw+++/X/aYHj16cPz4cSZPnuyRhGLt2rXMnDmTkSNHotVqK5RQxMbGlvs6UlJSWLJkCTt27KBr166SUAhRDpJQCHGVSEtLY8mSJc7tHTt2kJqayg033OAsCwkJITw8nAkTJnggQocGDRqwePHiCp/Hneto06YNixcvZvr06S4LrwkhrkwSCiGuEi1btnS5UU+fPp0dO3YUu3knJSUxevRo5/aZM2fYtm0bY8aMYe/evcTHx9O5c2eaNm2KzWZj69atpKam0qtXLxo1alTscxMTE9mxYwcBAQF0796dwMDAcsXt7uc3atSoxOu49dZb2bt3L6dOnaJDhw6VupaBEFczSSiEEC4ubSrYsmULd911Fx07dsRgMKDVatm0aRMLFizgs88+w8fHB5vNxt69e1m9ejW9evVynuv555/ntddeo2/fvuTk5HDo0CG++uorhg4dWuZ43P380q6jf//+ZGdn4+/vz/r163nvvfe45557KvErKMTVSUZ5CCGuKCcnh3vuuYcNGzawfv167r77bh555BEeeugh1q9fz8aNG7nxxhuZN2+e8z3Lly/nrbfeYs+ePfzyyy/8+eefzJs3jzvuuIO8vLwq//zSzjNmzBj++usvfvvtN1566SWefPJJt74mQghXklAIIa5Ip9Mxbdo053a/fv3w9vbmzjvvdCk7cuSIc3vRokV06NCBbdu28d133/Htt9/i5eVFUlISBw8erPLPL+08DzzwgHN78ODBnD9/nvT09HLFI4QoTpo8hBBX5O/vj16vd24bjUaCgoLQarUuZUVrHk6ePInZbOb77793OdeECRNc3ldVn1/W8wDlrjERQhQnCYUQokoEBgYSFRXFZ5995ulQhBDVQJo8hBBVYsSIESxdupTExESX8oSEBA9FJISoSlJDIYSoEjNnzmTFihX06dOHhx56iJCQEHbt2sWff/7JgQMHPB2eEKKSSQ2FEFepXr16MXz48GLll04IFRkZye233+5yTNOmTRkzZoxLWcuWLbnllluc2z4+Pvzxxx+89NJLxMXFsWvXLnr27MmuXbsuG1d2djaLFy9m7969Ffr8slxHcHAwEyZMwMfHB3BMvb148WKOHz9+2RiFEMVplFLK00EIIQTAH3/8wfvvvw/AbbfdVu3Tfx8+fJhnnnkGgL59+zJr1qxq/XwhajNJKIQQQghRYdLkIYQQQogKk4RCCCGEEBUmCYUQQgghKkwSCiGEEEJUmCQUQgghhKgwSSiEEEIIUWGSUAghhBCiwiShEEIIIUSFSUIhhBBCiAqThEIIIYQQFSYJhRBCCCEq7P8Bsl5B07ENY5kAAAAASUVORK5CYII=", "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-10-23T17:48:44.641815Z", "iopub.status.busy": "2026-10-23T17:48:44.641495Z", "iopub.status.idle": "2026-10-23T17:48:51.814645Z", "shell.execute_reply": "2026-10-23T17:48:51.813655Z" } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "resolved: solve time 6.883 s\n", "lumped: solve time 9.710 ms\n", "lumped with Jeffreson: solve time 10.630 ms\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAhQAAAFzCAYAAABmY5CaAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjExLjIsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvgI3uAAAAAAlwSFlzAAAPYQAAD2EBqD+naQAAZERJREFUeJzt3XlcVFUbB/DfzDAMOwgooLKogIqCuOOuuVsuuW8tauZSqVm9pfWa1ptmZaktVmpmuVWWuJu5iyKKgvsuooKI7PsyM+f9A50cAWWGgcvy+/qZmnvumXufMyzzcO6558iEEAJEREREpSCXOgAiIiKq/JhQEBERUakxoSAiIqJSY0JBREREpcaEgoiIiEqNCQURERGVGhMKIiIiKjUzqQMoD1qtFrGxsbC1tYVMJpM6HCIiokpDCIH09HTUrl0bcnnx/RDVIqGIjY2Fu7u71GEQERFVWrdv30bdunWL3V8tEgpbW1sABW+GnZ2dxNEQERFVHmlpaXB3d9d9lhanWiQUDy9z2NnZMaEgIiIywtOGDHBQJhEREZUaEwoiIiIqNSYUREREVGpMKIiIiKjUmFAQERFRqTGhICIiolKT9LbR0NBQaDSaQuWurq7w9vbWK7t9+zbu3bsHX19f3vpJRERUwUiaUHz44YfIysrSbefl5eHEiROYNWsW5s+fDwDIycnBmDFjsHPnTnh6eiI6OhoLFy7EG2+8IVXYRERE9BhJE4rdu3frba9fvx6jR4/GSy+9pCubN28ejh8/juvXr8PNzQ3BwcF4/vnn0aZNG7Rt27a8QyYiIqIiVKgxFCtXrkTnzp3RsGFDXdmqVavwyiuvwM3NDQAwaNAgNG3aFKtWrZIqTCIiInpMhZl6++bNm9i3bx9++eUXXVlsbCzu3buHli1b6tVt06YNIiIiij1Wbm4ucnNzddtpaWmmD5iIqCIRAsjPQn5WItLS7yI+4z5SMlKQnpOJ3LxM5ORnIzc/G7nqbOSpc5CryYdaq4FWaKERGgghIKCBVghohQYCAkIIaIUWWqGFkGkLtiEgHjmtDPi3RBQR1mPbevWLrFuCg+iKi9lRzTWr3Rkv9ptd7uetMAnFTz/9BAcHBwwdOlRXlpSUBABwcnLSq+vk5KTbV5QFCxZg3rx5ZRMoEVF5EwLIuIes2FO4EnMc5xOu4HbGPSTmpyIJOUiUa5EuB9LlMmQ/YXlpg8kePKhSUcQdleS8FSKh0Gq1+PnnnzF27FhYWFjoypVKJYCCgZmPys7Ohrm5ebHHmzVrFmbOnKnbfrhSGhFRpZCfA+3tMFy7sg2hscdwMu8+rpoBdx78TgQAKB48IHv4RI+VVkClBZRCBjMU/F8h5DCDHGZCDiVkUEAO2cN/MtmjW5DhwbZM9+xBGSCDrESJhqzYSrInbBUuE088EjOex/nX7ijJeStEQrF7927cvn0bEydO1Ct3d3eHXC5HTEyMXnlMTAw8PDyKPZ5KpYJKpSqTWImIyoRGDXFjPy5H/ozt945hl6U54szMCnIFy38TCVu1HDXyLGEnc4CdshYcLGujlm1duDjURc0adVDLwQ1utjXgYGUBC2XhRIOorFSIhGLFihVo27Yt/P399cqtrKzQvn17bNmyBWPHjgUAZGZmYs+ePZg7d64EkRIRmZg6DyJiDQ6HfYUVyhxEWFgAtlYAAKVWBodsR9jJ/FDftSM61muB1h4eqONgCbmcf5lTxSJ5QpGQkICtW7fiu+++K3L///73P/Ts2ROzZs1Cu3bt8PXXX6NWrVp49dVXyzlSIiITEgI4vwmR+/6L+RZqXLQ1B2ABhVYGywxP2Cp7YUTT3hjQzBM1bdnjShWf5AnF0aNH0a5dO4wYMaLI/V26dMH+/fvx7bff4vjx4/D398evv/4KGxubco6UiMhE0u4iadt0fJV0AsH2NgDModDKkZcchJZOwzBjUCsEujtIHSWRQWRCiCp/301aWhrs7e2RmprKabuJSDrqPGjClmFj+BIstbVAmqJgjIMmpQUCbV/EB31aw682f0dRxVLSz1DJeyiIiKqFq3twfvfb+FiZhfMO1gAAsxxnOKlfxtx+z6Kjj7PEARKVDhMKIqKylJOG3F3vYfGt7VhrZwshU8FMYwazjP6Y1ekVDGrmzgGWVCUwoSAiKivRR3Ej+FW8Y5mPK/YPuopT/fFc/Tcwe2wbWJrztk6qOphQEBGZmlYDcfBzbDz1DT5ztEeO3BwKtQVss8Zh6aBRaO5RQ+oIiUyOCQURkSmlxSL1rwmYl3UF/zg/SBwy62Ow1yy826s1VGbslaCqiQkFEZGpXNqB8B3TMMteiThrK8iEDHY5g/DdgDcRUJe9ElS1MaEgIiqtvEzk7nwPP9zcjJWOdtDKZFDk2aN/3dn4b88+MDcz4YJdRBUUEwoiotKIOYXzmybgv6osXHWwBwDY5rXGt30/QfO6bhIHR1R+mFAQERlDnYe8g59ixdmVWO5gC7XMHAqNCj1dp2NB7zEwU7BXgqoXJhRERIaKOYmzwZMxT5mGyzUKbge1VzfFkt4L0bJu8SshE1VlTCiIiEoqPxtZez/GqktrseJBr4SZxhy93N7AJ71eZK8EVWtMKIiISuLWMRwKnoTPLXNx80GvhKPGH0v6fIrA2uyVIGJCQUT0JHmZuL3tXXx7dwe221sDUMJcY4HB9d/ErM4jIZezV4IIYEJBRFSs/Kv78NuON7DMViDNxhoQgK95N3z9/FzUtnOUOjyiCoUJBRHR49LvIWzTdHyTfQqRDioAgIPaCW+3X4CBjdtJHBxRxcSEgojoIXUeYvZ+ibWXl2O9nSXUFiootTI84/IC5veaAXMzpdQRElVYTCiIiABknPsbO/a8jZU2asTaWwEA6qEhPnvuSzSqyUGXRE/DhIKIqrX8hBvYs3EK1siicKaGCoAZbLVWeK3lHIwJeFbq8IgqDSYURFQtidwMHN30Pn5P2o591pYAVDDXytG77gh80HUGrJRWUodIVKkwoSCi6kWTjzN/f4ng6z9hk60SamtLyATQyro95vf9GK42taSOkKhSYkJBRNWDVotrB1dh55lFWG8rQ7qdOQDAR+aBOb0/R6Crn8QBElVuTCiIqGoTAjdPbMK2Yx/hT+s8JDgoAACuWntMbz8XzzXsIXGARFUDEwoiqrKuRezB9pD3scUqDfH2ZgAUcNBa4KWmb2BcyzFQyBVSh0hUZTChIKIq5/LJ3dh+9CPssErEPTszAGZw0JhjZKOJeLXtBCgVnE+CyNSYUBBR1SAELoZtxe7j87HdOhV3HyQSNTRmGOb9MiZ1mAJzhbnUURJVWRUioYiLi0NISAisrKzQrVs3WFpa6vYdOXIEFy9e1Kvv6OiIwYMHl3eYRFQBCa0Wx/avxf7zi7HbOhuJ9goAZnDQmmF4/dF4teM0qBQqqcMkqvIkTygWLlyIjz76CF27doWVlRXef/99bN68GR4eBTPT/frrr/j777/RvXt33Wvq1KnDhIKomsvLV2P3zu9xJHo59tqokW0vB6BADa0Sw7zHYmK7qbAws5A6TKJqQ9KE4o8//sD777+PvXv3okuXLgCAGzduQKPR6NVr2bIlVqxYIUWIRFTBpKWlYOfWBTiSvA2HbWRQ28kAyFFXY4Ux/pMwouULUMo5RoKovEmaUCxcuBBDhgzRJRMAUL9+/UL17t+/j3Xr1sHe3h6tWrWCi4tLeYZJRBXAtWsX8Pfej3BKG4HjVuaArRwA0FhbAxOD3kGPRs9BJpNJHCVR9SVZQpGVlYWIiAhMnjwZFy5cQEREBGrXro0OHTrA3Fx/4NT169cRHByMmJgYRERE4IsvvsDUqVOLPXZubi5yc3N122lpaWXWDiIqO2qNFvsOBiPkwhKEWcQj1sIMgDnkQqC1mRde7/ZfBNZpK3WYRARAJoQQUpw4JiYGdevWRb9+/XDt2jW0bNkSJ0+ehFarxd9//63rqQgNDUWrVq2gVBZ0YS5fvhxTpkzB8ePH0aJFiyKPPXfuXMybN69QeWpqKuzs7MquUURkEveTUrB55+c4lbIdYVZq5MkLeh5stDJ0s2+DKc98AHcHL2mDJKom0tLSYG9v/9TPUMkSisTERDg7OyMwMBBhYWEwNzeHWq1Ghw4dULt2bWzatKnY19aqVQszZ87Ee++9V+T+onoo3N3dmVAQVWBarUDoiYM4eOpLRMiv4pLFvx2o9TUWGOI7CsODONCSqLyVNKGQ7JKHk5MTHB0d0bNnT90lDjMzM/Tq1Qtr1qx54mstLS2RlJRU7H6VSgWVireJEVUGt+4lYts/n+N8yt84YZ2HbCs5ADMohUCQwh0TOryNFvWe4fgIogpO0kGZAwcOxJkzZ/TKzpw5g3r16gEANBoNYmNj4e7urtsfGhqK27dvo127duUaKxGZTk6eGjsOBSP0ynKcM7+NO0oFYAsAcrhqzNDLpRvGd50FJ+uaUodKRCUkaULx0UcfoW3bthgxYgQ6dOiAsLAw7NmzB3v37gUACCHQt29ftG/fHk2aNMGtW7fw448/Yvjw4Rg4cKCUoRORgbRagbCzEdh/fDEuqSNw2lJAay0DoICFVqCdmSdGB72Btt692RtBVAlJNobioYSEBKxatQrR0dHw8PDA6NGjUbduXd3+3Nxc/Pbbb4iIiECNGjXQqVMndOvWzaBzlPT6DxGZlhACZ2/exq5DX+FyxiGctcpBtlyu2++ntsRzXv0wpONbsFLZShgpERWnwg/KLE9MKIjK19XYOOw68C3OJ+/BWcs0pCn+TSJc1HJ0smuGlzq/Ay8XfwmjJKKSqPCDMomo6hBC4Nztu/g75DtcTd6HixapSDaTAzYAIEcNjUB78wYY1mYSWvj05SUNoiqICQURGUUIgZNXr2F32DJcST+KyxbpyFD8m0TYaAVay2vj+YAx6BwwBgoFf90QVWX8CSeiEsvO02Bf5BGEnf0Z0bmROG+Zi1y5HLAGADkcNQKtFXXQr8lwdGr+ApRcLpyo2mBCQURPdCsxHduPrsO5mO24Lb+JKJUo+M1hBgByuKoF2ph7oX/TEWjtPxIKMy7MRVQdMaEgIj35Gi2OnIvE4dOrcSPjBK5ZpCJFIQcs/63jky9Hc+tGGNDiJQT4ckwEETGhIKr2hBA4fScWe8LX41r8PtyR3UG06sHNXw8uZVhrtQjU2CCoVhCebTcJNWv6SRkyEVVATCiIqqHoxFTsCvsNF2L+xi1tFG6q8qCWyYBHlsnwygeamXuiZ8Pn0L75C1CaW0sXMBFVeEwoiKo4IQQuxiVgb8QmXIvZi7ua64hW5SBLLgN0YyZlcFVr0FTUQEuXIPRu/TJqunKOCCIqOSYURFWMWqPFiRvXEHLmD1xPPIZYcQe3VfkFPRCPJBD2Gi2aaqwR6OCPngFDUd+7J2S8tZOIjMTfHkSVXFxqJvacO4hLUTsRm3EesYoExJg/GAOhW3RXhppqDRpprNDYvjE6NeqPgCYDITfjqrxEZBpMKIgqkew8NQ5dOomTl7cjOiUCcdpYxJjnIlf+4C6LR+7E8MzXwlc4IMDRH12a9IeXd0/IzDgvBBGVDSYURBVUanYuDl2NxNlr/+BOUgTua24jVpmFNMWD5OGRyxc2Wi188uWop6wNf5dW6BI4DDXdmgG8nZOIygkTCiKJCSFwOzkFhy/ux6VbhxGbcQnxIh53lY/0PCgfPCCDmRDwzhOoJ6uBRg6NEOT9DBo16gu5ZQ0JW0FE1V2JE4pXXnnFoAOvWLHC4GCIqrq7aekIvRaOS9FHEJtyAYl5t5EgT0O8UgOtrHDPg4VWi/r5Au6wg5dNfQR4tEcb/0GwsK8rVROIiIpU4uXLZTIZRowYUaKD/vbbb6hIq6Jz+XIqTxqtwKV79xB58yRuxh7H3bSLSMiPxX15OhLMHkkcHuOo0aB+vgLucic0sPdFi3od4NewBxS2tXnpgogkU9LPUIMSipImCYbULQ9MKMjUUrNzcfFeDK7EnMWt+EjcS7uGpNxYpCAViWa5yFQU/1pbjRaeaoHawhp1LGqjgXMTtPTpjrr1OgBKi+JfSEQkgZJ+hpb4ksfWrVtLfHJD6hI96mEiKiAKlRUqf+S5/lMBIQSEADRCCwFAKwS0WgHtg3K1Ros8jUC+Rot8jebB/wVy1HlIzknFvfRE3E+LQXp6DDKzYpCRF4d0dQLSRTrS5TlINdNC83inwWN3YNZSq+GmlsMNtg8SBz8E1u+Iup5tIbPieAciqlpK3ENRmbGHwnhht4/hjxMrcSXtEpJEBtTQQgNACwGtTO8jXUcU81yvThXowZcLgVoaDdzUMtSENVyUTqhrWw8+Ln5oXK8NbFyaApyumogqOZP3UABATk4Otm/fjiFDhhS5/88//8Szzz4LCwt221Z2t5Jv4t2tE3FOxBUUyB48qgFLrRZ2Wi0cNFrYaQErmMNBZgVHMwfUsnZDHQcP1HNpjLpuATBzasDLFEREMDChWLFiBa5du1ZsQnHo0CHExcXhtddeM0lwJI3Q6wfwzsFpSFUImAmBZzLUqC/zRC3rerAwt4a53BwKhRIKuRlkMnlBriEDIJNDJivIO+Qy2YMlrWUPygr+yR8+f/AiOf59nVwmK8hZZDJdHfkjzyGT6eoU7JM/GKso052vYF/Bc/mD55D9G4NcBshkcihkgEIuh0wm/7fhZioorJygsq4FWNYoeDBZICIqEYMSitWrV+OHH34odv9LL72EyZMnM6GoxC7GRuLNQ28gUwH45ObjZYeR6DdmNsyUnGGRiIiKZ1BCceXKFfj6+ha739fXF1euXCl1UCSNzNx0TNv5MjLNAL8cNT4I+hH+zTpLHRYREVUC8qdX+ZdGo4FarS52f15eHjQaTamDImnM/XMi4sw0cFFrML3Zp0wmiIioxAxKKJo0aYIdO3YUu3/Xrl1o0qRJqYOi8nf86gHszjsHABhq3gXtgwZKHBEREVUmBiUU48ePx5tvvol9+/YV2rdv3z68+eabBk/RTRXD14c+gFYmQ7tMOSaMXCp1OEREVMkYNIbi1VdfxYEDB9C9e3cEBgaiYcOGEELgypUriIyMxJgxYzBhwgSDgzh69Cj27dsHKysrjBgxAnXq1NHbf//+fWzYsAH37t2Dv78/hg4dCoXiCVMRkkHCruxHpFkqZELg+YYzoFQqpQ6JiIgqGYN6KGQyGdatW4cNGzbAw8MDp0+fxtmzZ+Hh4YHffvsNa9aseXCrYMlotVqMGzcOAwYMQHJyMpKTk9GvXz9cunRJV+f69evw9/dHcHAwNBoNZs2ahX79+nGshgn9ePgjAEC7LCX6dHtJ4miIiKgyknSmzG+++QbvvvsuIiMj4ePjA6BgRq7s7Gy4uLgAAIYOHYq4uDgcOnQIcrkcN2/ehK+vL1atWoUxY8aU6DycKbN4ccm30XdzX6hlMsyr+QoG95sudUhERFSBlPQz1KAeClP75ptvMHr0aF0yAQB2dna6ZEKtVmP79u0YM2YM5PKCUL28vNClSxcEBwdLEXKVs2rPfKhlMvjkajGw9+tSh0NERJWUQQnFnj174O3tDW9vb+zZs6dUJ05PT8fly5fRqVMnBAcH47///S+WLVuGe/fu6epER0cjJycHDRo00HttgwYNnjjfRW5uLtLS0vQeVJgQAodSjwIAWpk15bgUIiIymkEJxZQpU/DHH3/gjz/+wOTJk0t14tTUVADAokWLsGzZMiiVSmzbtg0+Pj44fvw4ACArKwsACnWx2NvbIzMzs9hjL1iwAPb29rqHu7t7qWKtqsIu78MdpRbmWoEhnd+TOhwiIqrEDLrLIyEhAZ6enhBCICEhoVQntrW1BQA4Ojri77//1pX37dsXs2bNwt69e2FjYwMASElJ0XttSkqK7vVFmTVrFmbOnKnbTktLY1JRhG0nVwIAAnKVaOjdXOJoiIioMjOoh+L111/HgAEDMGDAAEydOrVUJ7a3t0ft2rXRunVrvfLWrVvjxo0bAAAPDw9YW1vj8uXLenUuXbqExo0bF3tslUoFOzs7vQcVdjrrAgDAzypQ2kCIiKjSM6iH4uOPP8bx48eh1WoRFBRU6pOPHDkSBw8ehFarhVwuhxACBw4cQEBAAABAoVBg8ODBWL16NSZPngyVSoVz584hJCQEGzduLPX5q7Mrd87gprkGciHQt+2rUodDRESVnKS3jaampqJ79+7QaDRo164djh8/jvv372Pfvn26gZh3795Fp06dYGNjgxYtWmDbtm3o1asX1qxZU+Lz8LbRwhb+MRVrsg6jSQ6wYdJZqcMhIqIKyuS3jX7xxRclPnlJ69rb2yM0NBQfffQR/Pz8MGfOHFy+fFnvrg43NzecOXMGs2fPRvPmzfH7778blExQ0S6kRAAAfBT1JI6EiIiqghL3UMhkMpS0M8OQuuWBPRT61Op8dPqlOTIUMnxUexKe78n5J4iIqGgl/Qw1aAxFYGBgaeOiCuDQ2R3IUMhgrdWie9BoqcMhIqIqoMQJxeeff17ig44dO9aoYKh8HL74FwCgUa457GwdJY6GiIiqghInFG+//XZZxkHl6HLmRcAcaKDyeXplIiKiEpB0LQ8qf1qhRZSiYAbS5h49JI6GiIiqCiYU1czJK4eRoZDBQqtF59ZDpA6HiIiqCCYU1czRC1sBAPXy5LCzc5I4GiIiqiqYUFQzV5JOAwA8ZLUkjoSIiKqSUiUUiYmJOHHihKlioXJwSxsPAPCuESBxJEREVJUYlVCkpKRg0KBBcHZ2Rps2bXTlgwYNQmhoqMmCI9NKy0rFLaUGANC2UX+JoyEioqrEqITiP//5D3JzcwutAjp16lR8/PHHJgmMTO/I2Z3QymRwUmsQ2KSL1OEQEVEVYtBMmQ9t27YNx44dg4eHh15569atcfDgQZMERqZ3NjoEAFBXbQ6ZXCFxNEREVJUY1UORnJyMGjVqAChYt+OhjIwMyOUc51lR3Uor6FFyk3NAJhERmZZRn/4tW7bEpk2bAOgnFAsXLkT79u1NExmZXIxIAAB42jWUOBIiIqpqjLrksWDBAvTr1w+HDh2CEAL//e9/8ffff+Ps2bO85FFBqdX5iDHLByBDU89OUodDRERVjFE9FJ06dcLhw4eRnZ0NHx8fbNiwAfXq1UNoaKjeXR9UcUTcCEO2vGCGzDb+vaUOh4iIqhijeii+//57TJ48GWvXrjV1PFRGTl3+BwDgkS+DlbW9xNEQEVFVY1QPxbRp06DRaEwdC5Whm0nnAQCuwk7iSIiIqCoyKqHw8/PDyZMnTR0LlaG43LsAABfzOhJHQkREVZFRlzzGjx+PkSNHYs6cOfDz84O5ubne/sDAQFPERiZ0T5YOAPBwbCRxJEREVBUZlVBMnz4dADBu3Lgi9wshjI+ITC4vPw9xSi0AGZp6dZA6HCIiqoKMSiju379v6jioDJ2+EYZ8WcEdHgENecsoERGZnlEJhbOzs6njoDJ05sZhAECdfBnMLawljoaIiKoioxKKPXv2PHF/jx49jAqGykZ0QsEdHi6CyQQREZUNoxKKnj17PnE/x1BULHE5dwBzwMXMRepQiIioijLqtlEhhN5Do9Hg0qVL6NKlCye7qoDuIxUAUMfeW+JIiIioqjKqh+JxcrkcDRs2xIoVK/D8889j9OjRJXrdd999hx07duiVeXp64ttvvzWoDj1ZvEINQIb6rs2kDoWIiKookyQUD9WsWRNRUVElrn/mzBmkpKTgvffe05XZ2dkZXIeKF5scizRFwYqw/j68w4OIiMqGUQnFnTt3CpUlJyfj888/h5+fn0HHcnV1xXPPPVfqOlS001ePAACc1Rq41qoncTRERFRVGZVQuLu7F1nu6+uL9evXG3Ss06dPY/jw4bC3t0enTp0wduxYyOVyg+tQ0W7ERgIAaqnNAJlM2mCIiKjKMiqhuHjxYqGyGjVqwMXFsLsIlEolnnnmGXTp0gUxMTH4z3/+gw0bNmD79u2QPfjwK0mdx+Xm5iI3N1e3nZaWZlBcVUlM2jUAgBNsJI6EiIiqMqMSipEjRyIyMrLIfYGBgcXue9wnn3yiNx6iT58+CAwMxKZNmzB48OAS13ncggULMG/evJI1poq7nxsHmAO1lLWkDoWIiKowo64bnD59ushyjUaD8+fPl/g4jw+u9Pf3h6enp95KpiWp87hZs2YhNTVV97h9+3aJY6pqEkVB74ybLcdPEBFR2TGoh+LYsWNFPgcArVaL0NBQeHh4GB2MEAIpKSlQKpWlqqNSqaBSqYyOoyqJN8sHIEM9F3+pQyEioirMoISiXbt2RT5/yNXVFUuWLCnRsfLz8/HLL79g/PjxurEQ8+fPR2pqKgYOHFjiOlS8+2kJSH14y2gDrjJKRERlx6CEIj09HUDB4mAJCQl6+5RKpUG9AgqFAuHh4fjwww/h4+OD27dvIz09Hb/++iuaN29e4jpUvHM3jgMA7DVauNX2lTgaIiKqymRC4oU3EhMTce7cOdSoUQO+vr6wsLAwqs6TpKWlwd7eHqmpqdVqUqyfdi7EV/Fr4J0rsOnVc1KHQ0RElVBJP0NNOlOmMZycnNClS5dS16HCYpMLbhl1FIYlYERERIYy6i6PnJwcfPjhh/Dz84ONjQ0sLCz0HlQxJGTHAgBqyB2kDYSIiKo8oxKKOXPm4K+//sLs2bORmZmJNWvWYNq0aZDL5XjnnXdMHSMZKUmTBABwVnHZciIiKltGXfL47bffsHXrVgQEBOCFF17AkCFDMHToULRs2RLffPONqWMkIyXLsgAArvacg4KIiMqWUT0Ut2/fRpMmTQAA1tbWSE1NBQA8++yzCA8PN110VCqJCg0AoJ5LE4kjISKiqs6ohEIIAYVCAQDw9vbGP//8AwA4efJktbqLoiKLS41H+oM5KBp7tZI4GiIiquqMSijs7e11z9988028+OKLaNWqFfr27YtXX33VZMGR8R7OQeGg0aCWi7fE0RARUVVn1BiKlJQU3fOXXnoJDRo0wLFjx9CoUSM899xzpoqNSiEq9iwAwFkt57LlRERU5oxKKL7//ntMnjxZt92xY0d07NjRZEFR6cWnRQHgHBRERFQ+jLrkMW3aNGg0GlPHQiaUlB0HALCX2UgcCRERVQdGJRR+fn5PXD6cpJemTgYAOJg7SxwJERFVB0Zd8hg/fjxGjhyJOXPmwM/PD+bm5nr7AwMDTREblUIqMgEATla1JY6EiIiqA6MSiunTpwMAxo0bV+R+idcbIwAp8jwAMrjWqC91KEREVA0YlVDcv3/f1HGQCWm1WiQWTBMCT1c/aYMhIqJqwaiEwtmZ1+Ursjv3byJPXnCrqI9nc4mjISKi6sCoQZkAkJCQgJUrV+L999/XlR07dox3f1QA125HAAAcNFrY29WUOBoiIqoOjEooIiMj0bhxYyxatAjz58/Xlf/8889YvXq1yYIj49yOvwwAqKExOl8kIiIyiFGfOG+99RZmzpyJCxcu6JVPnjwZX331lUkCI+PFp94EADhoVdIGQkRE1YZRYyjCw8MRHBwMAJA9Mq2zj48Prly5YpLAyHhJ2XGAArCT20odChERVRNG9VAoFAqkp6cXKr906RIcHR1LHRSVTsrDSa3M+LUgIqLyYVRC0a9fP8ydOxdarVbXQxEdHY0pU6agf//+Jg2QDJcqOKkVERGVL6MSikWLFiE0NBRubm7QarXw9/eHt7c3cnJysGDBAlPHSAZKkecDAFwcPCSOhIiIqgujxlC4uLjg5MmT2Lx5M8LDw6HVajF79mwMGTKk0DTcVP6SzLQA5PCs2VDqUIiIqJowKqEAAHNzcwwbNgzDhg0zZTxUSvdS4pEpL+h48nYPkDgaIiKqLoxOKG7cuIFvv/0WFy9eBFCwAunrr78OLy8vU8VGRrhy6zQAwEajhYuzp8TREBFRdWHUGIqtW7eiUaNGOHjwIDw9PeHp6YkDBw6gYcOG2LFjh6ljJANE37sEAHDWAHjkll4iIqKyZFQPxVtvvYV58+Zh1qxZeuULFizAm2++iX79+pXoOHPmzMG6dev0ynx9fQslJZs3b8bSpUtx7949+Pv74+OPP4a3t7cxoVd5cbpJrZTSBkJERNWKUT0UMTExmDp1aqHyKVOmICYmpsTHiY+Ph4+PD3bt2qV7/PDDD3p1tm3bhqFDh6J///5YuXIlhBDo1KkTkpKSjAm9ykvMigUAOMBS4kiIiKg6MSqh8PPzw+nTpwuVnz59Gn5+hi2XbW1tDW9vb93D3d1db//cuXMxduxYzJgxA23btsUvv/yCvLw8LFu2zJjQq7zkvEQAgIOZvcSREBFRdWJUQjFmzBgMGzYMX375JQ4fPoxDhw7hyy+/xPDhwzF27FhERkbqHk8TEhKCgIAAdOrUCbNnz9abgTM9PR2nTp1C7969dWXm5ubo0aMHDhw4YEzoVV6qtuD9c1RxlVEiIio/Ro2hePPNNwEUjKV43PTp0/W2hRDFHsfR0RHvvvsuunTpgpiYGHzwwQfYtm0bwsPDYW5ujjt37kAIAVdXV73Xubq64syZM8UeNzc3F7m5ubrttLS0ErWrKkiVZQMAXKzrSBwJERFVJ0YlFPfv3zfJyT/++GMoFAoAQIsWLdC8eXPUq1cPGzZswIsvvgitVgsAUCr1Bxiam5tDo9EUe9wFCxZg3rx5JomxsklRFLwvdZ3qSxwJERFVJ0YlFM7OziY5+cNk4qG6devCy8sL58+f1ztPYmKiXr2EhIQnxjBr1izMnDlTt52WllZobEZVlJGbg/QHb6mXW2NpgyEiomrF6ImtACAzMxPJycmFyuvWrWvU8XJzcxEXFwcHBwcABVN8u7u7IzQ0FAMGDNDVO3LkCPr06VPscVQqFVQqlVExVGbX4qMBAGZCoLYrp90mIqLyY9SgzMjISDRv3hw2NjZwd3cv9CiJvLw8vPPOO7reh6ysLEydOhUajQYjRozQ1Zs8eTJWrFiBS5cKJmz68ccfcePGDUycONGY0Ku0W3cLZi110mhhZsNBmUREVH6M6qEYP348mjZtimXLlul6EwylVCpRt25dBAQEQKvVIjk5Gc2aNcPevXtRv/6/1//fffdd3LlzB82aNYO1tTUUCgXWrFkDf39/o85blcUl3QAAOGhknCWTiIjKlUw86TaMYlhYWOD+/fuwtbU1SRD37t2Dvb09LCwsiq2TmZmJpKQkuLm5wczMsDwoLS0N9vb2SE1NhZ2dXWnDrbAW/PEm1mXtQascBVZNipQ6HCIiqgJK+hlq1CUPHx8f3Lp1y+jgHufi4vLEZAIomADL3d3d4GSiOknNjgMA2HKWTCIiKmdGfTovXLgQL7/8Mj788EM0aNAAsse61xs1amSS4MgwaXlJgBKwVVTdXhgiIqqYjEooNBoNLl68iP79+xe534irKGQCGQ9mybRXOUkcCRERVTdGJRQzZszAmDFj8MYbbxg9KJNMLwMFs2Q6Wrk+pSYREZFpGZVQxMXFYdGiRbCxsTF1PFQKaXI1AKCWvYfEkRARUXVj1KDMRo0a4dq1a6aOhUpBq9UiVVFwqal2TR+JoyEiourGqB6KIUOGYNSoUfjoo4/g7e1daFBmYGCgKWIjA9xPT0SOvODr4Fm7icTREBFRdWNUQvH+++8DAIYPH17kfg7KLH83YgrWP7HSalHT2VPiaIiIqLqRdLVRMp2Y+CsAAAcNOEsmERGVO0lXGyXTuZ9SsDCYvZYTfxERUfkzalAmULCE+MqVK3WXPwDg2LFj0Gg0JgmMDJOYeRcAYCvMJY6EiIiqI6NXG23cuDEWLVqE+fPn68p//vlnrF692mTBUcml5CUBAGzl1hJHQkRE1ZFRCcVbb72FmTNn4sKFC3rlkydPxldffWWSwMgw6epUAIC9wl7iSIiIqDoy6oJ7eHg4goODAUDvllEfHx9cuXLFJIGRYdJFJgDAwYLjW4iIqPwZ1UOhUCiQnp5eqPzSpUtwdHQsdVBkuHRZHgDA2dpN4kiIiKg6Miqh6NevH+bOnQutVqvroYiOjsaUKVOKXTCMytbDabddHNwljoSIiKojoxKKRYsWITQ0FG5ubtBqtfD394e3tzdycnKwYMECU8dIT6HWaJCiKHhet2YDaYMhIqJqyagxFC4uLjh58iQ2b96M8PBwaLVazJ49G0OGDIG5OW9bLG/RSXHQPugp8nRtJHE0RERUHRmVUAQGBiIyMhLDhg3DsGHDitxH5Sfq3mUAQA2NBtb2XLqciIjKn1GXPE6fPl1kuUajwfnz50sVEBku5v51AEANDQAFZ8okIqLyZ9Cnz7Fjx4p8DhQsnx0aGgoPDw/TREYldj/1FgDAXquQOBIiIqquDEoo2rVrV+Tzh1xdXbFkyZLSR0UGScqMAwDYcdptIiKSiEEJxcO5J5ydnZGQkKC3T6lUQqVSmS4yKrHUvERADtjJrKQOhYiIqimDEgobGxsAQE5OTpkEQ8bJUKcC5oCtmZ3UoRARUTVl9GqjVHFkomDabXsVZyklIiJpVKhbAtRqNeRyOeTyf/McrVYLrVarV08mk0Gh4ADEh7KQCwCoYeUicSRERFRdVZgeirfeegtKpRIzZ87UK586dSrMzc1hYWGhe7Rs2VKiKCumjAfTbjvZ1ZY4EiIiqq6MSigCAwON2lecnTt3YteuXfD19S1y/+DBg6FWq3UPTpz1L61Wi3SFAAC4OXpJGwwREVVbkk9sdffuXUycOBFr1qyBpaWlMeFUa/cyUqB+MO22u4uPxNEQEVF1JenEVlqtFmPHjsX06dPRvHnzYuvt3LkTKpUKdnZ26NSpEz7//HM0aMBFsAAgOuE2AMBSq4WDEycVIyIiaUg6sdX8+fOh1Wrx1ltvFVvH19cXGzduRJcuXRATE4Pp06ejS5cuOHfuHBwcHIp8TW5uLnJzc3XbaWlpJY6psvl32m0toOJto0REJA3JJrY6fvw4Fi9erFut9OGdHEIIqNVqmJkVhPboIE0fHx+sX78erq6u+O233zBp0qQij71gwQLMmzfPkKZVWgkpBdNu22plwINLH0REROVNsomtwsPDkZKSAm9vb12ZRqPBmTNn8O233yI3N7fIW0Pt7e3h7u6Oa9euFXvsWbNm6SUiaWlpcHd3L3XMFVFyegwAwFYoJY6EiIiqM6NvG01ISMDKlSvx/vvv68qOHTsGjUZTotdPnTpV784NtVqNZs2aYdq0aVCr1cXOM5GcnIxbt26hTp06xR774XiLRx9VVWrOfQCAteC050REJB2jEorIyEg0btwYixYtwvz583XlP//8M1avXm2y4HJzc9G/f38cPXoUqampOHv2LIYPHw57e3uMGTPGZOepzNLzkgAANnJriSMhIqLqzKiE4q233sLMmTNx4cIFvfLJkyfjq6++MjoYMzMzvZ4JlUqF119/HR988AG8vLwwcOBA1K5dGydOnEDNmjWNPk9VkqUpGNfCdTyIiEhKRk29HR4ejuDgYAAF02A/5OPjgytXrhgdTHh4eKGy3r17o3fv3kYfs6rLFFkAADuu40FERBIyqodCoVDo7vh41KVLl+DoyA+28pQpe7COh3UtiSMhIqLqzKiEol+/fpg7dy60Wq2uhyI6OhpTpkxB//79TRogPVmGvGAQrJNdXYkjISKi6syohGLRokUIDQ2Fm5sbtFot/P394e3tjZycHCxYsMDUMdITpMsL1vFw4ToeREQkIaPGULi4uODkyZPYvHmzbmKq2bNnY8iQITA3Nzd1jFSMlKwM5MkLeohq16wvcTRERFSdGZVQAIC5uTmGDRuGYcOGmTIeMsCdB+t4mAmBWlzHg4iIJFTihCIkJKTEB+3YsaNRwZBh7iZGAwDstFrIVTYSR0NERNVZiROKTp06lfigQgijgiHDJKTeAQDYaMB1PIiISFIlTiiys7PLMg4yQnJ6HADARlv0NOVERETlpcR3eVhYWOgeQUFBetuP76PykZodDwCwMn4oDBERkUkYddvo6dOniyzXaDQ4f/58qQKiksvITQYAWIELgxERkbQM+tP22LFjRT4HAK1Wi9DQUHh48G6D8pKZnwooAGu5pdShEBFRNWdQQtGuXbsinz/k6uqKJUuWlD4qKpEsTUZBQqHgHR5ERCQtgxKKh+t3ODs7IyEhQW+fUqmESsWu9/KUJQoGytooHaQNhIiIqj2DEgobm4K/hHNycsokGDJMNvIAAHYWXJCNiIikZdTtAXv27Hni/h49ehgVDBkmS5YPAHCwrilxJEREVN0ZlVD07Nnzifs5sVX5yJRrAcjgaOsmdShERFTNGXXbqBBC76HRaHDp0iV06dIFa9euNXWMVAStVov0B/NZ1XTg0uVERCQtoxKKQgeRy9GwYUOsWLGCy5eXk6SsdKgfTLddp2Y9iaMhIqLqziQJxUM1a9ZEVFSUKQ9JxYhJKFgYTCkEatSoI3E0RERU3Rk1huLOnTuFypKTk/H555/Dz8+v1EHR0z1cadReo4VcZS1xNEREVN0ZlVC4u7sXWe7r64v169eXKiAqmcSHK41qucooERFJz6iE4uLFi4XKatSoARcXl1IHRCWTlnkfAGAlTHrVioiIyChGJRSNGjUydRxkoIycgoXBLARXGiUiIukZ/OetEALr16/Hc889h4YNG8LX1xf9+/fH+vXrOf9EOcrMSwUAWEApcSREREQGJhQajQbDhg3D6NGjkZCQgKCgILRr1w7379/H6NGjMXz4cGg0mrKKlR6RnZ8GALCQcf0UIiKSnkH95UuWLEFYWBhCQ0MRFBSkty80NBTDhw/H0qVL8eabbxocSFZWFm7cuAFnZ2e4uroW2p+SkoKEhAR4eHjA3Nzc4ONXNdmaTMAMsOTS5UREVAEY1EOxatUq/PDDD4WSCaBgOfPvv/8eq1atMiqQcePGwd/fH59++qleuVqtxoQJE+Di4oLOnTujVq1a+OWXX4w6R1WSoy1YoM1SwVtGiYhIegYlFFevXsUzzzxT7P5nnnkGV65cMTiIH3/8EXFxcfD39y+0b/78+di2bRsuXLiA2NhYLF26FOPGjUNkZKTB56lKckRBQmFlZiNxJERERAZe8lCpVEhPT4eFhUWR+5+0rzjnz5/H3LlzcezYMQwYMKDQ/h9//BETJ05EgwYNAAAvvvgiFixYgBUrVuCbb74x6FxVSc6DpcutVfYSR0JUsQkhoFarOb6LqBgKhQJmZmaQyUo3r5FBCUWbNm2wfPlyzJ49u8j9y5cvR5s2bUp8vOzsbIwYMQJffPEFPDw8Cu2Pi4tDTEwM2rZtq1ferl07nDx50pDQq5wcWcEvRxsLB2kDIarA8vLycPfuXWRlZUkdClGFZmVlBTc3t1KNUTQooZg1axZ69+6N+Ph4zJgxA56engCA6OhoLF68GN9++y12795d4uPNmDEDgYGBGD16dJH7ExMTAQBOTk565U5OTkhISCj2uLm5ucjNzdVtp6WllTimyiL7QUJhb+UscSREFZNWq0VUVBQUCgVq164Nc3PzUv8FRlTVCCGQl5eH+/fvIyoqCj4+PpDLjZsw0aCE4plnnsGKFSvw2muvYcmSJbrLGzk5ObC2tsbKlSvRrVu3Eh3r77//xu+//45t27bh3LlzuuMkJibi3LlzaNq0KczMCsLLy8vTe21ubi6UyuLnX1iwYAHmzZtnSNMqnWy5ACCDg00tqUMhqpDy8vKg1Wrh7u4OKysrqcMhqrAsLS2hVCoRHR2NvLw8g4cuPGTwNIsvvfQSnnvuOWzZskU3ANPX1xcDBgwo1JPwJElJSahTpw4mTZqkK4uOjkZCQgIiIiJw+vRp1KlTBzKZDHfv3tV77d27d4tdTwQo6EmZOXOmbjstLe2J9SsbrVaLzAcJpJN94Vtsiehfxv61RVSdmOLnxKh5m52cnDBu3LhSnXjUqFEYNWqUXllgYCC6du2KxYsXAwBsbGzQunVr7NixQ1c3NzcXe/bswbvvvlvssVUqFVSqqjvhU1pOFtQPum5rOnDpciIikl6FT93nzZuH9evXY+HChTh06BBGjhwJa2trTJ48WerQJHMvORYAIBMCTjVqSxwNEVUVs2fPLvUfiyXxySefFDt2jiqvCpVQeHt7w83NTa+sT58+2LZtGw4ePIg333wT9vb2CAkJgYODgzRBVgD3U2IAADZaATMV56EgItNIS0tDcnJymZ8nPT0dSUlJZX4eKl8VaqnKjRs3Flnep08f9OnTp5yjqbiSU+MAANZaABy1TkREFUCF6qGgkknJjAcAWAomE0RVza5du9CgQQNs2bIFQUFBqFOnDi5fvgwA2LlzJ7p37w4vLy906tSp0DIEFy9exNChQ+Hj44OgoCAsXboUWq1Wt/+nn35CmzZt4OnpiR49euDgwYPFxvH333/Dy8sLmZmZeuWffvqp3ozJT4sJAL766is0adIEAQEBmD59eqFjUtVQoXooqGTSsgq6Ci21zAeJDCGEQHZ++c+YaalUlHgOjJycHNy4cQMff/wxlixZgvr168PZ2Rl//fUXJk2ahKVLl6Jt27a4ePEiJk6ciLy8PLzyyisAgP79+6Nnz57Yvn07MjIysHbtWuzcuRPPPvssVq9ejWnTpuG7775D27ZtsW7dOvTs2RMRERFo0qRJoTi6du2K1NRUBAcHY8yYMbrylStX6sZZlCSm5cuXY+7cufjxxx/RrFkzfP/99/j222/Rq1ev0r6tVMHIhBBC6iDKWlpaGuzt7ZGamgo7Ozupwym1z357A7/mHECLHCVWTzoldThEFVJOTg6ioqJQr1493X31WXlq+M35u9xjufBRb1iZl+zvt+DgYDz//PM4c+aM3vpGDRs2xIwZMzBlyhRd2bJly/Ddd9/h7NmzyMnJgaWlJY4fP47WrVvr6gghIJPJ0KBBA4wZMwYfffSRbl/37t3h5uaGNWvWAABef/113LlzB8HBwQCAV155BbGxsdixYwcAICwsDO3atUNUVBQ8PT2fGhMA1K9fHxMmTMD777+vi6dp06Zwd3fHrl27DHkbqQwV9fPyUEk/Q9lDUQll5RXM/GkBLuNOVBUplUo0bdpUt52UlIQrV67gww8/xCeffAIhBIQQyMnJQXZ2NgDAwsICAwcOxPDhwzF+/Hh069YN7dq1g0KhQGZmJm7cuIFOnTrpnadLly7Fjl0DgDFjxqBXr16Ij49HrVq1sHbtWnTs2BGenp4liikjIwNRUVHo0KGD7pgymQwdOnTArVu3TPmWUQXAhKISylKnA3LAUlZ159ogKguWSgUufNRbkvMa4vFpwh8uJfDNN9+gc+fOenUfrffXX39h27Zt+Pvvv/Hyyy9DrVYjODhYt0zC4+s0mJubF5qJ+FFdu3aFq6srfv/9d0yZMgW//fYbPv744xLHlJ+fDwCFZjYuzXoRVHExoaiEsjVZBQmFwlLqUIgqFZlMVuJLDxWJq6srnJ2dcebMGQwfPrzYenK5HAMGDMCAAQOg1WrRt29ffPbZZ1i3bh1q1qyJiIgIdOnSRVf/1KlT8PX1LfZ4MpkMo0aNwtq1a+Ht7Y2UlBQMGzasxDHVqFEDzs7OOH36tF4vRWRkJGxseMt7VcNRfZVQjihYOdFSYS1xJERUHmQyGf7zn//gq6++wqZNm6DRaJCXl4d//vlH12MQExODGTNm6C4lpKam4v79+3B1LZief9q0afjiiy90aydt3rwZmzZtwvTp05947rFjx+LYsWP43//+h2effRY1atQocUwAMGXKFHz22We4fPkyhBBYtWoVjhw5YvL3iKRX+VJ1Qo4o6KK0UtpKHAkRlZe3334bZmZmeP311zF69GiYmZmhY8eOmD9/PgDAxcUF9erVQ+fOnZGcnIz8/Hz0798fc+fOBQC8++67iI+PR5s2bQp6aqyssHjxYnTv3v2J5w0ICIC/vz+OHDmCP//806CYgILZN2/cuIEmTZrA2toazZo1w4gRI5CSkmLS94ekx7s8KqFhP7bAJVU+plp3x5Shi6UOh6hCetKo9YosNzcXqampqFWr+JWEU1JSYGtrC4Wi6LEZ6enpsLa2LnLBJ7VajbS0NNSoUaPQrazp6enQaDSFZiJOTU1FdnY2atasWew5nxZTVlYWzMzMYG5ujoyMDOTn5+t6O0h6vMujmsqRFdxHb2vBH0aiqkalUj0xmQDw1KUHbG2L7700MzODo6OjQa+zt7eHvb19qWJ6dAl5jp+omjiGohLKfpBQ2Fs5SxwJERFRASYUlVDWg69aDZsn/xVDRERUXphQVDJarVaXUDg6cOlyIiKqGJhQVDLJWWnQPBhIVbMGEwoiIqoYmFBUMnGJdwAACiHg5OAmcTREREQFmFBUMgkpMQAAa62AXFl5boUjIqKqjQlFJZOcdg8AYK2VOBAiIqJHMKGoZFKz7gMALIXsKTWJiIjKDxOKSiY9KxEAYCkMW72QiCqHCxcu4OWXX5Y6jFK7fv06Xn75Zd2Ko1T1MaGoZDJyUwAAFoKTnBJVRbGxsVi9erXUYZTa/fv3sXr1amg0GqlDoXLChKKSycpPAwBYwlziSIiIiP7FP3MrmWx1OiAHLGS8w4OoOjh06BC2bNmCL774Qld29epVfPLJJ1i+fDmUSiVOnTqFH3/8Ee+++y7++OMPREVFoXnz5pgwYQKio6OxevVqJCYmolu3bhgyZIjuOA9f98EHH2DNmjWIjo5Gq1at8PLLL+st8qVWq/H777/jyJEjsLGxQc+ePdGjRw+9OJOSkrBs2TLcvXsXTZs2RePGjcv+zaEKhT0UlUy2JgsAYKGwekpNIqoKrly5go0bN+qV3bt3T+9ywq1bt7By5Ur07NkTeXl58PLywrvvvovnn38effr0gUqlgpubG15++WWsXLlSd5yHr2vfvj3S09Ph5eWFDz74AK+88oquTn5+Pnr37o2lS5fCx8cHDg4OeOmll/DRRx/p6mRlZSEoKAjbtm2Dj48PQkJCMHLkyDJ+Z6iiYQ9FJZOjzQEAWJlZSxwJUSUkBJCfVf7nVVoBsrK9M0utVmPVqlXo1KkTACAvLw9z5szB8ePH0bp1awAFy1CvXr0aEyZM0HvdRx99pBsI2r17d7Rp0wbTp09HYGAgvv/+e8TFxSEiIgLm5gWXWnv06IH27dtjypQpqFmzJpYtW4asrCycOXNGt/T1yJEj8dtvv5Vpm6liYUJRyeQiFwBgpSx+TXoiKkZ+FjBfginrZ8cC5mX7R4BSqUSHDh102w0aNIC1tbUumXhY9nhvBwC9yyCtWrWCl5cXDh8+jMDAQGzfvh1arRZTp06FEAJCCGi1WqjVapw/fx5du3bFwYMH8eyzz+qSCQAYPnw4E4pqRvKEIicnBxEREcjIyEDjxo1Rt25dvf0RERG4fv26Xpm9vT169uxZnmFWGNkouAXLRmUvcSREVJGYm5tDLv/3KrZcLtf7gH9Y9vhdF0qlEra2tnpljo6OuH+/YM6bxMREeHh4oGPHjnp1unbtCm9vbwAFd3Q0adKk0DGoepE0ofj5558xb948eHh4wMzMDKGhoRg7diy+//573Q/GDz/8gG3btiEoKEj3Og8Pj2qbUOTICn4Z2Fryh5XIYEqrgt4CKc5rJJVKhby8PL2y5OTk0kakk5+fj7t378LNrWBtICEEbt++DXd3dwBA3bp1kZWV9cS5MTw8PHDr1i29sujoaJPFSJWDpAmFtbU1IiMjYW9f8Nd2ZGQkmjdvjgEDBuC5557T1QsKCiqym646ypZrAchgb11T6lCIKh+ZrMwvPZiaj48P4uLicO3aNXh7e0Or1WLVqlUmPcdXX32Fzz77DADw66+/Ii0tDc8++ywA4MUXX8SQIUOwc+dO9O3bF0BB0vHbb7/pBl4OHjwYEyZMQFRUFOrVq4ecnBx89913Jo2RKj5JE4phw4bpbfv6+sLMzAypqal65ampqdixYwfs7e3h7+8PO7vqO34g40GPprOdq7SBEFG5CAoKQv/+/dGhQwd06dIF58+fR506dUx2fJVKhSNHjqBjx46wtbXFvn378Pnnn6N27YKxJs8//zzmzp2L559/Hm3btoWTkxPOnDmDNm3a6BKKh+MlWrZsiS5duuDcuXPw8PAwWYxUOUg+huLevXs4fPgw0tLS8Msvv6BXr14YOnSoXp2TJ09i8eLFiImJQWxsLJYuXYoXXnih2GPm5uYiNzdXt52WllZm8ZentOxMZMkLRorXqdVA4miIqCw0adKkUA9EcHAwQkJCEB8fj2bNmsHOzg47duyAUqkEALRo0QI//PCD3mvatm2Lb775Rq+sS5cucHBw0CszMzPDwYMHceLECURHR2PJkiXw9fXVqzNnzhy8+uqrCA0NhUwmQ/PmzeHp6anbL5PJ8NdffyEkJARxcXFo0qQJnJ2dsX37dl2MVA0IiUVERIghQ4aILl26iJo1a4oFCxYIrVar2793716RmZmp2/7888+FUqkUZ8+eLfaYH374oQBQ6JGamlqmbSlrp29EiKY/NxWBq5oITV621OEQVWjZ2dniwoULIjubPyvF2bRpk7C2tpY6DKoAnvTzkpqaWqLPUMkntgoMDMTGjRtx4MAB7Nu3Dx9//DF++ukn3f5nnnkGVlb/Dmh66623YGdnh507dxZ7zFmzZiE1NVX3uH37dpm2obzExl8DANhrBORKzpRJREQVh+QJxaOaNm2KVq1a4eDBg8XWkclksLW1RXx8fLF1VCoV7Ozs9B5Vwb2UglHUdlouXU5EpVfUpRIiY0mWUGg0GiQlJemVZWZm4urVq7q5KLRaLRISEvTqRERE6Oabr26S0wtud7PRSj70hYiqAA8PD4wZM0bqMKiKkOyTKT8/Hx06dMDAgQPh5+eHlJQU/Pzzz7CwsMC0adMAFCQdHTt2xIABA9CkSRPcunULS5YsQe/evQsN3KwOUrILJpqxgUriSIiIiPRJ1kNhYWGB48ePo3bt2ti/fz+uXLmC1157DRcuXICra8EtkUqlEidOnEDdunVx4MABJCcnY/ny5di5c6feSnjVRUZ+wWQ21vLKdR89ERFVfZL2ndva2up6I0pTp7rIUKcBCsDWjNNuExFRxVKhBmXSk2WiYJVEexWn3SYiooqFCUUlkomC+fztrVwkjoSIiEgfE4pKJEOuBgDUtDPdtLtEVLHcu3cPGzZskDqMUktMTMSGDRug1WpLfaySvCcxMTH4888/S32ukjp16hQ2btyIsLAwAAUzNO/btw+//fYb7ty5U25xVCRMKCqR1AfjUF0d60kbCBGVmbNnz2LUqFFSh1FqV69exahRowqtlGqMx9+TuLg4/P7773p1Tpw4gZdeesmg4yYnJ2PDhg2F1o96mhkzZqB///74/fffcerUKeTm5qJly5Z45513sGnTJsTExBh0vKqCExpUEmlZ6bp1POq6ekscDRFR+XF1dcWIESN025GRkRg9ejSGDx9equNev34do0aNwtmzZ3WrXpfE8uXLsXHjRt3qq/v27cP169eRmpoKc3PzUsVUmTGhqCRu3r0CADATAm7O7KEgqi5u3LiB8+fPo3///rqyhIQE7NmzB8OHD4dcLkdMTAyOHTuG559/HufOnUNUVBQCAgJQr149aLVahIWFITExEa1atdLdlg9A73VnzpxBdHQ0AgMD9Rb+eiguLg4nTpyAjY0NWrRoUeQH8MmTJ3H37l00adLkqe36888/ERQUpFs59dixY4iNjcXgwYMBFCyR/vvvv6N79+6oWbMmBg0aBKCgV+HgwYMQQugugzy6mJkQAufOnUN0dDT8/PxQv379ErzLJWtrfHw8du7ciaysLJw8eVLXsxEaGgqlUom//voLADBy5Ejdezt48GCEhYXh9u3b6NWrl+5YZ8+exfXr1+Hu7o7AwMBCUyFER0fj3LlzcHJyQosWLfQSFY1Gg7CwMNy7dw8+Pj5o2rSp3mvDwsIghEBAQAAiIiKQkZGBtm3bFloYztSYUFQSMfFXATxYx8OMq/cRVRf79u3D//73P72E4tKlSxg1ahQGDRoECwsLnDhxAmPGjEHz5s2hUCigVCpx+PBhfPnll1i3bh3Mzc0hhMDp06exe/duBAUFAYDude3bt0dqairs7Oxw9OhRfPPNN5g4caLufAsWLMDChQsRFBSE7OxsXLhwAWvWrEHv3r11dcaOHYutW7eiffv2uHjxIvz8/J7Yri+++AL9+/fH7NmzAQBTpkxBZGQkYmJiULt2bZw5cwajRo3C/fv3dZc8Ro4cieTkZISGhkIIgeDgYABA7969UaNGDajVavTp0wdpaWmws7PDgQMHCrXlaZ7U1oSEBGzduhUAcPDgQZw7dw4AcOvWLeTm5uriGTlypO697dSpE1JTU1G/fn20a9cOcrkcw4cPx7lz59C8eXNcuXIFtra22Lp1qy7Z+9///ofPPvsMHTt2REZGBtLS0vDnn3+iQYMGuHfvHvr06YPExEQ0adIEx44dQ8+ePbF+/XpdUrJs2TKcPn0aOTk58PT0RFxcHO7evYuDBw+iUaNGJX4vDGb6NcsqnpKulFaRLfr9LdH056Zi2A/+UodCVCkUtXqiVqsVmXmZ5f54dAXlp/nnn3/Eo7+aly9fLjw9PfXqHD58WADQtW3Tpk0CgPj66691dSZNmiQAiOXLl+vKxowZI/r376/bfvi66dOn653PyspKxMXFCSGE2LFjh3B2dhZRUVG6Oj///LOoVauWbiXoP//8U1haWorLly8LIYTIzMwULVq00Ivxce+9957o2bOnEEKI5ORkoVQqRfPmzcXatWuFEEIsXrxYBAQEFPme7Ny5UygUCr3jPWzL4sWLdWWLFy8Wjo6ORZ5fCCFOnDghAOhWry5JW9PT0wUAERoaqquzfv164eTkVGQ8H3/8sV75q6++Kvr27Styc3OFEEJoNBoxePBgMWbMGCGEEHl5ecLc3Fzs3r1b95orV67oYnz55ZdFy5YtRUZGhhBCiBs3bgg7Ozvxww8/6Oq/9NJLwtLSUly4cEEIUfB9/8wzz4jx48cX+16YYrVR9lBUEjdSzgDmgIe8ptShEFVa2epstF3XttzPGzY6DFZKq6dXLAW5XI5XXnlFt92uXTv89NNPGDdunF7Z119/Xei177zzju75+PHjMXv2bGzbtg0TJkzAqlWr4Ofnh/DwcJw4cQJCCCgUCsTHx+P8+fNo3bo1fv/9dwwcOFB36cHKygqvv/46xo8fX2y8Xbt2xdKlS5GXl4eDBw8iICAA/fv3x/79+zF69GgcOHAAXbt2Nfg9mDRpkt45ZsyYgYSEBDg7Oz/19SVpq6EenZhRrVZjzZo1mDp1KrZs2QIhBIQQqFu3LjZu3Khrg0qlwtmzZ9G9e3fI5XL4+PgA+Pcy0Pfffw9r64IZk+vVq4fRo0djw4YNePXVV/Xa3rhxYwAFi2p27twZ//zzj8HxG4IJRSVxR9wDAHg7BEgcCRFVRJaWlrCwsNBtP1x1+dFr8yqVCjk5OXqvUyqVunEMQMEHmoeHB6KjowEAN2/eREZGhu4D76ERI0bAzKzgI+TWrVvo0qWL3v569Z481qtjx47Iz8/H8ePHdclD165dMWHCBAghcPjwYSxfvtyAd6Do9wBAoTYXpyRtNYSVlZXeatfx8fHIyspCZGQkbt++rVf3YfKkUCjw66+/YubMmViwYAG6dOmCkSNHYujQobh//z6ysrIKjQtp0KABdu/erVfm6Kg/AWJRX3tTY0JRCaRlpuGWuQaADG0b95M6HKJKy9LMEmGjwyQ5r7HkcnmhuRxM+cGQn5+PjIwM2NjY6MqSk5N1f9Hb2dnB29sb69atK/YYTk5OSE5O1it7fPtx1tbWaNmyJQ4cOIADBw7gf//7H4KCghATE4MdO3YgKSkJnTt3LkXLDFeSthpCJpPpbdva2kImk2HixIlPvENl4MCBGDhwIK5evYqdO3filVdewa1bt/DGG29AoVAUWqk7KSmpRD0wZY3zUFQCByM3I18mg71Gi2aNukkdDlGlJZPJYKW0KvfH4x8shqhTpw7u3buH9PR0Xdn+/ftN8XboPBxMCBTckhkVFYUOHToAAPr06YOtW7cWmlvh0e2OHTtix44denNOPLzj4Um6du2KTZs24fz58+jUqRNUKhWCgoIwd+5c+Pv7w8nJqcjX2djYQKPRID8/35BmPlVJ2loatra2aN++PX744QcIIYo8R3Z2ti4Z8/HxwbRp0zBw4EAcO3YMSqUSbdq00XtvNRoNNm3ahI4dO5okxtJgD0UlEBm1FwDgla+CXMEvGVF10qVLF7i5uWHo0KEYMWIETp8+XahLvjQUCgXee+893Lx5E7a2tvjiiy8wfPhwtGzZEgDw+uuvY+vWrWjbti1ee+01ODk5ISIiAnv37sWVKwW3s7/22mtYtmwZevbsidGjRyM8PBy7du166rm7du2KTz/9FK1bt9ZdGujatSvmzp37xEUh/fz8YGVlhXfffRdt2rTRu220NErS1tL67rvv0L17dzzzzDMYNmwYcnJysHfvXtSvXx9ff/010tPT0b59ewwYMAD+/v64c+cO/vzzT/z0008AgEWLFuGZZ56BQqFAmzZt8PvvvyM9PR2zZs0ySXylwR6KCu7cjVPYnXcCAOBj0UDiaIiorD0+iZOFhQVCQ0PRqlUrhIaGwtfXFzt37sSIESN04yPq1q2LoUOH6h3H09NTN6fDQw0aNNC7/fTh8Q8fPoz8/HxERETg3Xffxa+//qq3f+/evVi4cCFu3ryJU6dOoXnz5jh9+rSujo2NDcLCwtCpUyecOHECfn5+2Lt3r16MRenYsSNGjhyJKVOm6MoGDhyIESNGYNiwYcW+J46OjtizZw/UajU2b96M06dPF/ke2NnZYcSIEbCyKnpA7MMejofjLkrSVqVSiREjRuhdYijqvS4qHgAICAjA+fPn0bdvXxw/fhyxsbGYNm2abrBsrVq1EB4eDg8PD4SEhCA1NRW7d+/WXSJp164dTp48CUdHR4SEhKBr166IiIjQi6dt27Zo166d3nn9/Pz0bvMtCzLxeL9LFZSWlgZ7e3vdfdaldfjUFgSf+hYCgEDB2ycgAPHolt6ex/YVXUf/WIAWWlw1S0eymRxeecCakXthb1ur1PETVQc5OTmIiopCvXr19Abq0b+Cg4MxduxYZGRkSB2KJNauXYsXX3wRmZmZ1f575Ek/LyX9DGX/uRGu3z2N3crYcjqbHE5qLT4IWsRkgojIBGJjY7FlyxZ89tlnGD16dLVPJkyFCYURvOu0RO+7RwEZIIMMD/8LmQz//vt3G4//V/bw/4/uKSiXQaZ3HE9nPzzfeQpsrB0kai0RVVXFdctXdffu3UN4eDjefvttg2bRpCfjJQ8iqpJ4yYOo5ExxyYODMomIiKjUmFAQERFRqTGhIKIqrRpc1SUqNVP8nDChIKIqSalUAgCysrIkjoSo4nv4c/Lw58YYvMuDiKokhUIBBwcHxMfHAyhYqKk0U2ATVUVCCGRlZSE+Ph4ODg5PnIjsaZhQEFGV5erqCgC6pIKIiubg4KD7eTEWEwoiqrJkMhnc3NxQq1Ytky8kRVRVKJXKUvVMPMSEgoiqPIVCYZJfmERUPA7KJCIiolJjQkFERESlVi0ueTy8vzYtLU3iSIiIiCqXh5+dT5urolokFOnp6QAAd3d3iSMhIiKqnNLT02Fvb1/s/mqxOJhWq0VsbCxsbW1Ndh96Wloa3N3dcfv27Sqz4FhVa1NVaw/ANlUWbFPlUNXaVFbtEUIgPT0dtWvXhlxe/EiJatFDIZfLUbdu3TI5tp2dXZX4RnxUVWtTVWsPwDZVFmxT5VDV2lQW7XlSz8RDHJRJREREpcaEgoiIiEqNCYWRVCoVPvzwQ6hUKqlDMZmq1qaq1h6Abaos2KbKoaq1Ser2VItBmURERFS22ENBREREpcaEgoiIiEqNCQURERGVGhMKI6SnpyM8PBw3b96UOhSjZGVlITIyEjExMcXWEULg/PnzOH36NDQaTTlGVzoXLlxASEgI8vLyCu3Lz89HREQELl68KEFkxomJiUFkZGSR7QGAzMxMnDx5ElFRUeUcmXHu3r2LkydPIjo6utg6iYmJOHHiBOLi4soxspIRQuDUqVOIjIwstk5+fj4iIyNx8eLFYqcqLkmd8qJWq3Hs2DFcvny52Drx8fGIjIxEampqsXUefi/euHGjLMI0SGZmJo4cOYJbt249sZ4QAkePHi3265mcnIwTJ04gNja2DKI0TEJCAkJCQpCQkFBsHSEELl68+MSvZVxcHE6cOIHExETTBynIIMuXLxdWVlaiYcOGwtraWvTt21dkZGRIHVaJ3L17V7z00kvC3t5eBAYGiho1aogOHTqIqKgovXoXL14Uvr6+wsXFRdStW1fUrVtXhIWFSRO0AU6dOiUsLS0FAHH79m29ffv27RMuLi7Cy8tLODk5iYCAAHHz5k2JIn26mJgY0b17d2Fvby9atWol6tWrJ3bs2KFX55dffhE2NjbC19dX2NjYiB49eojU1FSJIn6ypKQk0b17d2FraytatmwpHB0dRYsWLQp9782ZM0eoVCrh5+cnVCqVmDBhgtBoNNIE/Zgvv/xS+Pr6CgcHB9GsWbMi6xw8eFC4uroKT09P4ezsLJo2bSpu3LhhcJ3ykJGRIf773/8Kd3d3YWtrK0aMGFGoTkhIiGjfvr2oVauWCAwMFJaWlmLSpElCrVbr1Vu7dq2wtbUVvr6+wtbWVnTr1k2kpKSUV1N0YmNjxRtvvCHc3NyESqUSH3744RPrL1y4UMjl8iK/np988onue9HCwkK88MILIj8/v2wCf4ILFy6IF154Qbi5uQkAYv369UXW27t3r/Dy8hLu7u4iMDBQBAUFiejoaN1+tVotJkyYICwsLHQ/X3PmzDFprEwoDHDmzBkhl8vFunXrhBBCxMfHCy8vL/HGG29IHFnJHDt2TKxevVr3Q5Geni46d+4sOnXqpKuj1WpFQECAGDRokO4X+YQJE4S7u7vIycmRJO6SSE9PFw0bNhTvvPNOoYQiNTVVODk5if/85z9CCCHy8vJEt27d9NpdkeTl5YmAgADRq1cvkZ6eLoQQIjExUe8XyaVLl4SZmZn46aefhBAFH9g+Pj5i4sSJksT8NDNnzhTu7u4iMTFRCCFEVlaWaNGihRgyZIiuTnBwsFAqleLIkSNCiII22tvbiyVLlkgS86PUarWYMWOGuHTpknjrrbeK/ABKT08XNWvWFDNnzhRCCJGfny969Ogh2rVrZ1Cd8hIVFSXmzp0rYmJixLPPPltkQrFq1Spx9OhR3faFCxeEg4ODWLhwoa7s6tWrQqlUih9//FEIIURycrJo2LChGDduXNk34jFHjhwRS5YsEcnJyaJBgwZPTCjCwsKEp6enGDt2bKGv565du4RCoRD79+8XQghx/fp14ejoqNfu8vLnn3+K1atXi/T09GITinPnzgmVSiXmz5+vKzt58qQIDQ3VbS9evFg4ODiIy5cvCyGEOHz4sDAzMxObN282WaxMKAwwc+ZM4e3trVf26aefCnt7+0IZe2WxYsUKYWZmpksejh8/LgCI8PBwXZ2bN28KAGLr1q1ShflUL7zwgpgxY4b4559/CiUUv/zyi1AqlSI5OVlXtmvXLgFAXL16VYJon2zt2rVCJpM9sQdl9uzZom7dunplixcvFlZWVhUy8XvhhRdEjx499MomTJggOnfurNseMGCA6NOnj16dV155pdjeAKkUl1CsW7dOKBQKkZCQoCvbs2ePACAuXrxY4jpSKC6hKMqgQYNE//79ddtz5swRbm5uQqvV6sq++eYbYWFhIbKyskwea0k9KaFISUkRDRo0ELt37xbTp08v9PUcPny46Nq1q17Z66+/Lho2bFhG0T5dfn5+sQnFqFGjnvpzEhAQICZPnqxX1qNHDzFw4ECTxcgxFAaIiIhAy5Yt9cratGmD1NTUCnHd0BgnTpxAvXr1dAu+REREQC6Xo3nz5ro6np6eqFWrFiIiIqQK84nWrFmDyMhIfPrpp0Xuj4iIQP369eHg4KAra9OmjW5fRbN3714EBATAw8MD58+fx6VLl5Cfn69Xp7jvxaysLFy5cqU8wy2RmTNn4vz585g/fz727t2Lr7/+Gtu2bcOcOXN0dYpr07lz5wq1vyKKiIiAl5cXnJycdGWPf5+VpE5Flp+fj9OnT8Pb21tXFhERgRYtWugtvNimTRvk5OTg0qVLUoT5VK+++ioGDBiAnj17Frm/uO/FK1euICsrqzxCNMjevXvx3HPPITs7GydPnsTt27f19ufn5+P8+fNFtsmU33fVYnEwU0lKSkLjxo31yh7+YkhKSpIipFLZu3cvVqxYgVWrVunKkpKS4ODgUGhFOScnpwrZxqtXr+LNN9/E/v37i50dLikpSe8XOABdGytim2JjY2FtbY327dsjLS0NWVlZyM3NxYoVK9CvXz8ABW1q0qSJ3usq8vein58fXnjhBXz22Wdo0KABoqKiMGjQIN2HKVD018nJyQkajQZpaWmF9lU0RcVva2sLpVKp+5qUpE5FNmvWLCQmJmLatGm6sqSkJDRo0ECvXkX+Xly+fDkuXbqEX375pdg6xX0vCiGQnJwMKyursg6zxDQaDeLj43H79m00bNgQTk5OiI6ORqNGjbBhwwZ4eHggNTUVGo2myDaZ8mvEHgoDKJVK5OTk6JVlZ2cDAMzNzaUIyWgnTpzA4MGD8fbbb+OFF17QlRfVRqCgnRWxjRMmTEC/fv2QkpKCkJAQnDt3DkBB+x7e+VBUm/Ly8qDVaitkm5RKJY4ePYpp06bh/PnziIqKwosvvohRo0YhOTlZV6cyfS9Onz4dwcHBuHbtmu4vqOvXr2PUqFG6OpWtTY8rKn61Wg21Wq2LvyR1KqpFixbh22+/xcaNG+Hl5aUrr0xft7i4OMyYMQNTpkzBiRMnEBISgtjYWGRmZiIkJER3F0tlapNCoYBcLsf27dtx6NAhRERE4NatW9BqtZg0aRKAgvYAKLJNpmwPEwoDeHp6FrrV8uG2h4eHFCEZJTw8HL169cLEiRMLXSbw9PREVlYWUlJSdGVqtRrx8fEVso0uLi64fv063nvvPbz33ntYuXIlAGDBggXYunUrgMr3dfPy8oKdnZ3eh+2kSZOQlpaG06dPA6h8bdq2bRtGjhwJZ2dnAIC1tTVefvll7Ny5U3c5o7g2OTg4wNbWttxjNpSnpydiY2P1bgN9uP3wa1KSOhXR4sWL8cEHH2DTpk2FLhNUpu/F7OxsNG/eHGvWrNH9zjhy5Aju3r2L9957D1evXgVQfJusrKwqZE+Zp6cnevfurUv0bGxsMGbMGBw+fBhAwdLjDg4ORbbJlF8jJhQG6NmzJw4dOqR3L/bmzZvRvHnzCvlNVpSTJ0+iZ8+eGD9+PL744otC+7t27QqlUoktW7boyv755x9kZWWhR48e5Rlqifzxxx8ICQnRPb766isAwF9//aXrlu3Zsyfu3buH48eP6163efNm2NjYoF27dpLE/SS9e/culNTduXMHAFCzZk0ABW06evSo3j3pmzdvRuPGjVGnTp1yjbckatasqWvDQ7dv34aDg4Pur6eePXtix44devOebN68udjr3BVNz549kZCQgNDQUF3Z5s2bYWVlhQ4dOpS4TkWzdOlSzJo1C3/99Rf69OlTaH/Pnj0RFhaG+Ph4XdnmzZvh4+MDT0/P8gz1qerVq6f3+yIkJATDhg2Dt7c3QkJC0KpVKwAFbdq1a5fe2J3Nmzeje/fuhS4HVwS9e/culCzcuXNH9/sCAHr06KH7Iwso+ENx+/btpv35MtnwzmogOztb+Pn5ic6dO4vg4GAxZ84coVAoxM6dO6UOrUQuXrwoHB0dRbdu3cThw4f1Hnl5ebp6s2bNEjVq1BArV64Ua9asEXXq1BHjx4+XMPKSK+ouDyGEeP7554W3t7fYsGGD+O6774S1tbX47LPPJIryybRarejevbvo0qWL2LJli1i/fr3w9fUVzz33nK5Obm6uaNasmWjfvr3YtGmT+Oijj4RCoRDBwcESRl68VatWCYVCIebOnSt2794tvvzyS2FjY6M3Cj82NlbUqlVLDB06VGzZskVMmjRJWFlZibNnz0oX+CMiIiLE4cOHxciRI4W3t7fuZ+fRO7yGDRsm6tevL9avXy++//57YWNjo3crX0nrlJcjR46Iw4cPi/bt24vu3buLw4cP6805s2LFCgFAvP/++3q/LyIjI3V18vPzRYsWLURQUJD466+/xCeffCIUCoXYuHFjubcnJydHF2OdOnXEuHHjxOHDh8WZM2eKfU1Rd3ncv39f1K5dWwwcOFBs2bJFvP7668LCwkKcPHmyjFtQWHJysjh8+LA4cOCAACDmzp0rDh8+LK5du6arEx0dLZydncWMGTPE33//Lb744gthaWkpli1bpqtz5swZYWVlJSZPniy2bNkihgwZImrVqiXu3r1rsli52qiBEhISsHDhQkRERMDJyQmTJ09Gt27dpA6rRHbu3IlPPvmkyH1btmyBo6MjgILZ1lauXIng4GCo1Wr06dMHr7/+OszMKv4Y3vDwcMyYMQObNm3Sy85zc3OxePFi7Nu3DyqVCiNGjMCYMWMkjPTJsrOz8dVXX+HgwYOwtrZG165dMXnyZL3rncnJyVi4cCHCw8Ph6OiIiRMnVui/5v/55x+sXbsWMTExqFWrFgYNGoShQ4fq3R0QFRWFzz77DJcvX4aHhwfefPNNNGvWTMKo/zVu3Dhdl/ijdu3aBRsbGwAFY3OWLFmCPXv2QKVSYdiwYXpjlEpap7x0794dubm5emVOTk7YvHkzAOC///0v9u/fX+h1jRo1wooVK3TbKSkpWLhwIU6cOIEaNWrglVdeQe/evcs2+CLExcVh6NChhcqbN2+Or7/+usjXLF26FBEREXqD0wHg1q1bWLhwIS5duoQ6depg+vTphe6SKA9hYWF46623CpUPGjQIb7/9tm77xo0b+OKLL3D58mXUrl0bY8eOLfQ1OH36NL766ivcunULDRs2xH/+8x/Uq1fPZLEyoSAiIqJSq3gXg4iIiKjSYUJBREREpcaEgoiIiEqNCQURERGVGhMKIiIiKjUmFERERFRqTCiIiIio1Cr+TEVEVK7i4uJw6NAhDB8+vNzPfePGDd0U6U2aNIG/v7/RxzKmHQkJCdizZw8AwN3dvcJOh01UETGhIKomEhMT8c8//zyxTsuWLXH9+nWMHj1akoRi3759mD59Ovr37w+5XF6qhCIyMtLgdiQmJiI4OBjh4eEIDAxkQkFkACYURNVEcnIygoODddvh4eFISkpCr169dGWOjo5wc3PDiBEjJIiwQM2aNbFhw4ZSH8eYdjRs2BAbNmzA5MmT9RZeI6KnY0JBVE14e3vrfVBPnjwZ4eHhhT684+LiMHDgQN32nTt3cPz4cQwaNAhnz55FdHQ0AgIC4OXlBY1Gg7CwMCQlJaF169ZwcXEpdN67d+8iPDwctra2aNGiBezs7AyK29jzu7i4FNmO559/HmfPnsWtW7fQpEkTk65lQFSdMaEgIj2PXyo4duwYXnrpJTRt2hTm5uaQy+U4cuQIFi9ejF9++QWWlpbQaDQ4e/Ys9uzZg9atW+uO9fHHH2PRokUICgpCVlYWLl26hHXr1qFHjx4ljsfY8xfXjvbt2yMzMxM2NjY4ePAgvv/+e4wbN86E7yBR9cS7PIjoqbKysjBu3DgcPnwYBw8exMsvv4w33ngDU6dOxcGDBxESEoK+fftiwYIFutds3boV33zzDc6cOYNdu3bh0KFDWLBgAV544QXk5OSU+fmLO86gQYNw9OhR7N69G/Pnz8d7771n1HtCRPqYUBDRUykUCkyYMEG33a5dO1hYWODFF1/UK7ty5Ypue9WqVWjSpAmOHz+OP/74A7///jvMzMwQFxeHixcvlvn5izvOq6++qtvu2rUr4uPjkZKSYlA8RFQYL3kQ0VPZ2NhAqVTqtlUqFezt7SGXy/XKHu15uHnzJnJzc7Fx40a9Y40YMULvdWV1/pIeB4DBPSZEVBgTCiIqE3Z2dvDw8MAvv/widShEVA54yYOIykSfPn2wefNm3L17V688JiZGooiIqCyxh4KIysT06dOxbds2tG3bFlOnToWjoyNOnTqFQ4cO4cKFC1KHR0Qmxh4KomqqdevW6N27d6HyxyeEcnd3x5AhQ/TqeHl5YdCgQXpl3t7eeO6553TblpaW2L9/P+bPn4+oqCicOnUKrVq1wqlTp54YV2ZmJjZs2ICzZ8+W6vwlaYeDgwNGjBgBS0tLAAVTb2/YsAHXr19/YoxEVJhMCCGkDoKICAD279+PH374AQAwePDgcp/++/Lly/jwww8BAEFBQZgxY0a5np+oMmNCQURERKXGSx5ERERUakwoiIiIqNSYUBAREVGpMaEgIiKiUmNCQURERKXGhIKIiIhKjQkFERERlRoTCiIiIio1JhRERERUakwoiIiIqNSYUBAREVGp/R/WLUVfEiG3ewAAAABJRU5ErkJggg==", "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 }