{ "cells": [ { "cell_type": "markdown", "source": [ "# Graphene band structure" ], "metadata": {} }, { "cell_type": "markdown", "source": [ "This example plots the band structure of graphene, a 2D material. 2D band\n", "structures are not supported natively (yet), so we manually build a custom\n", "path in reciprocal space." ], "metadata": {} }, { "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "n Energy log10(ΔE) log10(Δρ) Diag\n", "--- --------------- --------- --------- ----\n", " 1 -11.15663652643 -0.60 5.9\n", " 2 -11.16017635789 -2.45 -1.30 1.0\n", " 3 -11.16040527607 -3.64 -2.29 3.1\n", " 4 -11.16041638521 -4.95 -2.94 3.3\n", " 5 -11.16041688914 -6.30 -3.24 3.1\n", "Computing bands along kpath:\n", " Γ -> M -> K -> Γ\n", "\rDiagonalising Hamiltonian kblocks: 5%|▊ | ETA: 0:00:07\u001b[K\rDiagonalising Hamiltonian kblocks: 10%|█▌ | ETA: 0:00:05\u001b[K\rDiagonalising Hamiltonian kblocks: 15%|██▍ | ETA: 0:00:04\u001b[K\rDiagonalising Hamiltonian kblocks: 20%|███▏ | ETA: 0:00:04\u001b[K\rDiagonalising Hamiltonian kblocks: 24%|███▉ | ETA: 0:00:03\u001b[K\rDiagonalising Hamiltonian kblocks: 29%|████▋ | ETA: 0:00:03\u001b[K\rDiagonalising Hamiltonian kblocks: 34%|█████▌ | ETA: 0:00:03\u001b[K\rDiagonalising Hamiltonian kblocks: 39%|██████▎ | ETA: 0:00:03\u001b[K\rDiagonalising Hamiltonian kblocks: 44%|███████ | ETA: 0:00:02\u001b[K\rDiagonalising Hamiltonian kblocks: 49%|███████▊ | ETA: 0:00:02\u001b[K\rDiagonalising Hamiltonian kblocks: 54%|████████▋ | ETA: 0:00:02\u001b[K\rDiagonalising Hamiltonian kblocks: 59%|█████████▍ | ETA: 0:00:02\u001b[K\rDiagonalising Hamiltonian kblocks: 63%|██████████▏ | ETA: 0:00:01\u001b[K\rDiagonalising Hamiltonian kblocks: 68%|██████████▉ | ETA: 0:00:01\u001b[K\rDiagonalising Hamiltonian kblocks: 73%|███████████▊ | ETA: 0:00:01\u001b[K\rDiagonalising Hamiltonian kblocks: 78%|████████████▌ | ETA: 0:00:01\u001b[K\rDiagonalising Hamiltonian kblocks: 83%|█████████████▎ | ETA: 0:00:01\u001b[K\rDiagonalising Hamiltonian kblocks: 88%|██████████████ | ETA: 0:00:00\u001b[K\rDiagonalising Hamiltonian kblocks: 93%|██████████████▉ | ETA: 0:00:00\u001b[K\rDiagonalising Hamiltonian kblocks: 95%|███████████████▎| ETA: 0:00:00\u001b[K\rDiagonalising Hamiltonian kblocks: 100%|████████████████| Time: 0:00:03\u001b[K\n" ] }, { "output_type": "execute_result", "data": { "text/plain": "Plot{Plots.GRBackend() n=25}", "image/png": "iVBORw0KGgoAAAANSUhEUgAAAlgAAAGQCAIAAAD9V4nPAAAABmJLR0QA/wD/AP+gvaeTAAAgAElEQVR4nOydd1hTWROHfwm9tyQURVCxooI0EQRFFFGx97rqKtbVxV37WtaKrqtgx7Z2V8QGNop97Qp2URRUFAWkE0JIme8PXD9XKUkIRb3vk4cn3MyZmZty554558xhEREYGBgYGBi+V9jV7QADAwMDA0N1wgRCBgYGBobvGiYQMjAwMDB81zCBkIGBgYHhu4YJhAwMDAwM3zVMIGRgYGBg+K5hAiEDAwMDw3cNEwgZGBgYGL5rmEDIwMDAwPBdwwRCBgYGBobvmq87EN69e3fHjh1VaZGI0tLSiuvSSaW4cwf37/9H4NUrPH4svXUr7uzZ28nJyQCEQowdi6goSKVV6WmNRijE1Kng88uXlEpx+zYEAkWsiES4fRsiUakCRFTDSwxevRp/6xbt24cTJzLfvUurbnfKQiKRVKU5Pp8fFfX01i307y/4558npYmJRKLIyHu3b+OPP3IWLHimmK337zFjBqr2/CCRYPx4PHxYpUYV4MSJ+15eIldX6a5dd6SfXOOWLsWqVXj58v+SDx8+vH79urTKr4My/cbpa2b37t2DBw+uMnMvXrywtnYyNe1iZeUUGPiiUSNq1YpOnPi/wP79xOPRyZP05MmTqKio/Px8IpJKaft2cnQkGxtasYLS06vM3xqKWEx9+lC/fiSRlC88dSq1bUtCoXwm4uLioqKiJk4U+vmRVFqqWGFhYVFRkXyqqxYrKxdAAhAQNnfusup2pyxyc3Or0pxIJGrbtheP56Opae/ufqQMyXHjZnC5dhYWTnXqnJ88mcRiuW0VFVHHjuTvr7i3irF+PZma0oQJdPjwWXt7HxcXvxs3bla1E2WSnp5naOgGEEAczrA7d+58fOnKFRo9mjgccnam5cupf/9fOZyeJiajXF27iBX4DBRFIpHI8htnAqEcDBw4kcWKBgiIrlVr0vnz/39JKqXAQLK0pNjYUpvfukX+/mRsTP36UXQ0Xb163d6+Y+PGnocPR1SB8zUEqZTGjKH27amwsHzhkBBq0oQyM+UzMWnSHBOT3trav6qru799W1CGZM0PhEOHzmWxJgCHOZxWd+/erW53yqKKA2Exqamp2dmFTk60alX5wrm51LUr+fhQdrbchvh8at2afvtNAR8rREYGjRuXxWY7Aq+BBHNze5FIVDWmpVLpuXPnoqKiPotbr19TaCgFBJCrK2lrC1RVHYrv1TQ0ut++Hf+ZErGYYmJozBgBm92qOF7q6485e/Za1ZwCMYFQXqRSaXJyMp/P/+x4Tg6dPk1z51L79qSq+iNwAyDgRvfuP36UKSykoUPJ1ZXevSvfUFYWBQVR48YSNTV74AWQweW6vpOl5TfBjBnUqhXl5ZUvefw4mZlRQoJ8+qVSKY/nAEgB0tGZffLkyTKEa3ggzMujxo0lEycemDlzUQ2PglRNgbCYly/J3JzOnStfUiymSZOoWTNKSpLbSno6NW5Mf/whiIiIuHjxotztFSUuLs7YePS/UaT7y5dvqsauj89AI6OxhoZTnJ19r16VBAXRgAFkaUk8HnXvTsuW0YULxOfTsmVruVxHDsetceNfGzQouScgEok4HDtABJCaWl9t7Xv29jRpEu3dS69efS785MmT69evS2TJF8kAEwjlgM/n29m15/G68XgtIyJOvXlDoaE0eTK5u5OBAbm70+TJFBpK587Fmpq2NDGZbGraMi4urrhtSgq5uNCAAVRQVt/jc96/zzAy8in+chsb/3TlypWKn0XNZ8UKsrWl9+/Ll7x1i7hcunpVESumpvZAHkDGxj/8888/ZUjW8EA4bFg1pOMUphoDIRHFxJCFBSUnyyQcEkJmZqRALHv2TKCm5qapOdPY+IchQybK3V4hBAJBnTqOLNYhFZU92tpuDRtK9+yRaVihIqSlpfF4HYsvUCzWkEaNHo4fT7t2lXxjyufzc3JyiOjAAeLxKDi4hPGI4OAtHI4dl+s6cmSASEQPHlBICA0bRnXqkLk5+flRYCBdukQ//TSPw+liYjLK2dlXKX1fJhB+oKCg4MCBAxEREaUlpgUCWrRou5raCoCATFVVNzMz6t2bVq2ia9fos/cwKyvrn3/+ycrKKv733j2ytqb588saiCqN5s3bqaltZLH2amg4ZmZ+3hP99ti5k+rUoZcvy5dMSiILCzp8WBErYjHZ2R3T0WnB5ToPGza5bOGaHAi3bqVmzeiLDEXNpXoDIREtW0atWsmUciei06eJx6Pdu+UzcebMGX39X4rDA4/Xqsq+PMnJyT//PHfmzEVpaWlnz5K7O9na0sGDilx2ZCQ3l6+l5QAUARIOp+3r169lbPj8OTk7U48elJHx+UtCoTCvpFzQ06e0YweNGUNNmohYLMd/0zkT1669mJJSshWpVHr8+PFdu3Zll5fm/l4CoZmZ4+HD4aUJCIXCxo3dtbQW6utP9fbul5lJN2/S33/TkiU0ahS1bUu1a5OGBpmabmazgwEC8mvXdpXR+smTxOPR3r0KOp+dnb106arZs5d07fpmwIBKv8WrLvz9p/N49lyus4nJxfjPRxBKICeHmjen4GAFzQUEUMeOJBJJhTJMsKmxgfDBA+Jy6dGj6vZDHqo9EEql1K8fjR0rq/z9+2RtTTNmyBFO7ty5w+H0BqRAnpGRnbTyAlF5REeTiws1a0ahocoPhzk51L072dhs53Cac7l2v/8uwwDsJ4hENH8+1akjd59bIpFwOHZAEUAaGkNbtLjB45GODtnZUZ8+NGMGbdlC585RcjL16TNaX3+SuvoKKyvHMr54iYmJHToMvHTpQbmmWVSzp4+XzZ49e8aOddTQ8N+5c4+WlhWfDz4f+fnIyUFBAfh8PHt28+TJ7ULhRgBstreu7mEbG4P69VGvHurX//CwtEReXo6LS+esrEbA46CgXwcP7luaxZSUlLi4uJYtWx48aLFyJQ4fhrNzRc+isBA+PnBxwcqVFVVV07h48WLPntuysnYCaebm3VJSrpchLBaLiVQ7d4adHf78UxFz+/dj9mzcugUTE5nkhUIhm81WU1NTxFilUViIVq0wdSp++KG6XZGHvLw8PT296vUhPx+urggIwI8/yiT//j169YKFBXbsgJaWTE2mTJm3b98xIpZYvOjQoW7e3hXxt6LExGDmTIhE+O039OsHADk5OSwWS19fX2GdT5+iZ094eGDdOqiqEhGx2YqssouJwYgRGD0ac+dCRUXWVhs37pg/PxjQ8/JqeuDAJgC5uXj27D+PhARpamoropsAVFXne3u7NG7cVVsbBgbQ04O2NnR0YGQEbW0MH+71/PmSixddPDxUy7b71QfCESP6Ey2qX9/dysq3+C3Q04O+PnR0oK0NieTlsmUj8/KigEIOx/3t29uqqiW/IyKR6NGjRxYWFlwutzRzly9f6dVrSmFhF5HohKXlupgY1zp1lHMiWVlo0wZjxuDnn5WjsIYQFhY2dOhToXA2AFNTx3fvbpcolpub6+3d/9WrHIGA5eS0KzraRvZfzkfu3YO3NyIj4eAga5OaGQh//BEFBdi/v7r9kJOaEAgBPH0KT09ERMh6hyoU4scf8fw5/v5bcO9ejImJiZubmywNL19G797Yvx/t21fI4QpChOPHMX8+NDTA4cy9fv0MkXTUqG7Ll89RQNvJkxg1CsuWYeRIJfiWmorhwyESYc8eWFjI2kooFAoEAkNDwzJkatVyTEk5DZjo6fWZMmUeh9OSz0duLvLywOejoABZWeDzcfOms1h8c9Qo6bZt5cVy+fquNYzdu3drax80N3d4X/oEjMDAdWZm9rVrO4aFlZpBlRE/v5FAHEBAXOfOIyuo7TNev6Y6dWjXLuVqrWZOnMhSUXHU0FhqZDR03LgZpYn99tsyVdWtAAHXvb0Vmf2UmUn169POnfK1qoGp0b//pgYNqLqzjIpQ7anRjxw7RlZWlJYmq7xUSrNn56uquurozDY2HjJq1FQZG166RDwenT2roJ9KRCKhdeveqKi0Lx6/5HLblHFJLBGplBYtIktLun5dmY6JxbRgAVlY0L59b0aNmtq7t/+naw0V5syZ89bWLmZm9nPmBJYh5uXVV109+OLFnHIVfvWBsHlz9ydPnlSNOT+/n4EIgFisiB9++Fnp+h88IB6PoqOVrrh6iIwkDodOnMg5cODApUuXypAcO3YGcBogIMXBobO8hiQS6tKFfvlFbg9rWiBMSCAul/6dj/yVUXMCIRHNmkXt25Pssw5PnjyppTX731kwTrLP3b94kXg8mVZuVDZJSUlcbrfiU1BR8Zk/P+XL6SqlkZ9PfftS69ZU2uSUCnL+PKmqurNY4cBlMzPHDNk9qxj5+fmLFq18+DCxXMlyMqc1n+bNrRo2bFgFhvbuxbVr02vXHlBUtMzERCUwMFTpJmxtERqK/v0RGQl7e6Wrr1IiIzF0KMLC0LatPtC/bOF69X5QURmsqXlVSyt67txZ8tqaOxcCAQIDFfW1ZiAUon9/LF781X/0NYHFi+HnhzlzsHy5TPImJiba2okCAYAcsVgs+6iYhwfCwtC3L0JD0batwv4qAWtrazc3kytX+gOSFi3qJSaa29igRw+MHw8Xlw8yT58+zcvL+6xhSgp++QVNm+LPP5GSgpQU5fvGZhfo6RVmZVkAyMtrEhYW5ujoqIAeFotlb28v+6ejo6Mze3aALMX/vvoxwlOnTu3du7dSrYjF+O03HDmCw4dha4uioiJ1dfXKMxcail9+wT//wMqq8oxULqdPY8QIHD0KV9fyhcPDMW4cDh5Mz8y81rx5c2tra7lsHTuGKVNw8yZKH9stlRo1RjhhAt6/R6jy76+qiBoyRviRzEy4uGDpUvQv5zbsA+PGzTx06BTAlkiWhYT4Fs89kZHoaAwejEOH4OmpmLNK4+HDhyoqKo0bNwaQnY3QUAQHQ1MTY8di8GBYWhpZW1urfDICL5GgoADq6tDQqFzH8vP5UqkmwAIE6uramposBZTEx8efO3fOWZ4JilKpVCKRlP8br/zuaSVSBSXW0tLIy4u6dqV/lw5WBevXU9OmJSzE+So4eZJMTemabEWUzp4lHo9u3VLQVnw88Xh044aCzWtOavTgQbKxUaT0V82hRqVGi7lzh3g8elD+5Pn/8PAh1a5N27fL1yoyUo6vfVUilVJ0NPXrR8bGpKWlVwM/JtlxdXW9KmeVDRnXEX7du09UNrGxcHGBiwvCw1HmJCYlM2ECOndGly4oKKg6o0rhxAmMHInwcLRqVb7wjRsYOBAHD0KhNAny8tC7N5YuVcIKluolMRETJ+LAARgYVLcr3xZ2dli1Cj4+p1xdew8YMD5Ftqxf06aIicH8+VizRg5bPj7YtQs9euDGDQW9rSRYLHTogNBQ3L1b3a7UYJhAWCq7d6NzZ6xahcBAKLSQpkL88QcaNcKgQVW9+UtFOHQIo0YhIuL/YxJl8OABevTAX38pmE0iwqhRaNdO1hVjNRaRCEOG4Lff5Fj1wSA7zs5PMzOXXb++8uDBbn5+sq4JaNQIly5h7Vr51rP6+GDHDnTvXuNiYTG1a6OUtWMMYN6YEvg4KHj+PJo0qR4fWCxs3YquXTFq1Dtb279NTAyHDRtcqWOTFSQsDBMn4sQJODmVL/z8OTp3xurV6NJFQXNLluDNG1Ty6HBVMH06OBxMmlTdfnyj3L17l8XyA+oR1UtJmS97QysrXLqEDh2Qni7HPCxfX/z1F7p3R0QEOJwkHR0dHo+niN8MVQsTCP/Pu3fvjI2N379X79sXJia4caOaU1Vqati+PadevS5i8SRNzeR9+4acOXOwOh0qndBQBAQgOhotWpQvnJICHx/MnYuBA+U2lJ+f/+bNm6Sk+ps2qV6/jhp8Y1AOeXl5iYmJz57ZHD6sExsLliJTBxjKx9nZWVf3T4GgHYv1WF/fXK62ZmY4exY+PpBKsXy5rJ9R587Yvp3atBmmo5Onqpo1bly3hQunKeI6QxXCBEIAEAgEHh49k5NVRKIUFZWVP/7YYenSakiHfsnDh9c0NTvn5Y0SCPDggZtIJKohUxwBSCSSqKgoIsrI8Jk5UzUqCra25bfKyICPD8aOhb+/3BbPn784cGCAWNw4O/vZkSPHa9WSf55ozeDu3bu+viNFopZZWbG7d+8zMammtMN3gLW19bFjQX/8sVVPz/zUqZ3Xrsk0k/kjPB7On0fnzpgwAevXy3pNMDWN1dBgZ2UdA6SbNrn89tuUmpzLYQATCIvZvz/0wQMvoXAmkGlh0T0wsEN1e/SBunXramj8kZdXCKQCRTUnCgLw9Oz58GFDkYgtEq29ffukrW35N8x8Prp3R8+emD5dEYtTpy5NTT0G1Gax9l6/vrVbN7lXHNYQ5swJevduM+AE/HPs2NrBgzdUt0ffMq1bux4+7Arg1Cn07YurV2FpKUdzQ0NER6NHDwwfjh07ZB1m09BA8YI9kegrXp/2/VADej3VTVERoqOFRUW6AABNdXVxNTv0CQ0bNvz99yFWVp5WViOl0k1JSdXt0L8kJyc/e6aSk/NnQcEf+vr6mprPym1SVITevWFnh8WLFTT6cckriwXgK76+qKqqAkIALFahmhpzM1pFdO6MgAB07w4+X76GurqIiEB6OoYMgUhUvryDg4Ozs5TH625o2FYkGhQeznQHZYKI0j8hMzOzykx/1z/CoiLs2IHFi1G/fv9atXwLCx+wWPcXLKhZda8nTBg5YcJIAOvXw88PV67UiEn2enr6+fkpgAhgsdnJRkZGpUnm5uZGRkaampqvWdOGy8W6dYobtbGZ/ehRDwODJnp6CVOmHFdcUbVCBDW1XzU0BhsYWGtqvly8OKy6PfqO+OUXJCRg2DCEhck39qGtjYgIDBwIe/uQ3NytOjpaW7cGtmlTcnluFot18uSepKQkHR2d1FRe167IysKYMco5hW+Y7OxsHo/H4/GKa8eYmZnFxcVVjenvNBAKhdi5EwsXokkThIXBxcWwqOji3bt3LS0tzczMqtu7kpk4EfHx6N8fJ05U8zRooRATJxqYmo7Pz3disTBr1iQOh1OiZHZ2tr19h7S03hJJuLn58YSEQIVHXrdswcOHnomJF/LyXtvY2JS2i0jNZ9o0vH7d6N27azk5KbVq1fp6T+QrZe1a+PpizhwsWyZfQ3V1LFr0xNHxkFB4GXg/aFC35OSSd1Mppm7dugB4PJw9i06dkJmJGTMq4rgSkEqr2YEvyc/PV1VV1dTU/HgkLi7OQva9KpTEd/cj5POxdSv++AMtW+Lo0f/P9VdXV5erck+1EBSEnj0xfjy2bKk2HzIy0KsXTE3x8OFILa1yFmadOXMmLa23QDAboMJCJzU1BeuBnj6NefNw6RJq1dIFGiumpCawYAGio3H+PAwN1QwNv9oael8zamo4eBCurqhfH6NHy9c2JeWVlpajUKgOWAiFGjJWW7SxwT//wNcXaWlYubLaZgjHx0Ncg4Z9cOvWrYEDBxYVFUmlUg8Pj/0ybDwmECA1tVKc+fYD4cGDx+bM+UNLS2v16sX377dasQIODggP/3z9cnZ29oULF+rVq9e8efMytB05Eh4ZecXX171nz26V63dJqKhg3z60aYM1azB5ctXbx7Nn6NoVvr5YvVqmzJKxMUcqjQAAvNfUlH+PQQDAw4cYORKHD8PGRjEFNYX167F3Ly5dQumJZIaqwNgYERHw9ISNDdq1k6Ohq6urnt7MvDwrFiulsNAyJ0ddxgq35ua4cAHdu2PECGzbVg0ZnfR0+Pl9vkFuekH62utriyRFVeCAtaH1OKdxnx4JDg4eOnToggULAHxaCvz333/X0dEB0KFDhy7/XWjMYmHECHTrhg4dLsyfP72wsGjUqAFLlsysuHvfeCB89+7dxInL09NPAbkdO3YfNOj2mTPsxl/0KN6+fevs3CUnp6eGxqYZM/ymTZtYora//tobEHAsJ2fs339vCgrKGzFicIlieXl5Cxb8+eDB84kTB3bv3lW5Z6Snh/BwtG6NunXRrWpj8ZUr6NMH8+dj3LjyhQFIJNi1q62R0XGx2F5bW+2vvxTZdf7tW3TpglWr0Lq1Aq1rEHv2YPlyXLyImpp6/75o1Ah79mDgQPzzjxw3WHp6erGxkfv2hZqYNHn5cq6nJ6KiZJ2DamiIyEgMGICePREaCm1thX2Xm8JC9OqFQYMQHPyf44/SHwX+EyiSyjD/p8JY6FmMdhityv5/xOnZs2dQUJC1tXWDBg3c3d0/Hre0tNTX1wfw5d68mpr4+29cuYKRI6dKpacBkw0bevbtG9eyZcuK+qdY8dMawu7du52dnfPy8j47npJCx4/TwoXk6XkFmFC8R5emptvHfbDy8ykxka5epfBw2r6d/Pw2sFibAQKKtLVdOnemDh0+PFxcyNHxw0NLawgQDxAQb2w8uG9f8venmTNpxQraupUOH6Zz5+juXfL2Hq6mthG4b2zc4YbCBaHL5OZNMjWlu3crQ3fJHDhAXC6dPCmrvEhEgwdThw7E5ytulM8nFxdavlxxDWVTZUW3jx0jc3N69KgKTFUbX2M15y1bqEkTxWudr11L1tb09KkcTUQiGjGC3N0pM1NBo/IildKgQTRoEEmlpKdXg4puh4eHOzk5TZw4cceOHcVHiqeJvnnzprQmjo6uW7Zcff5coqrasPiqDsy1tj7m708hIXTrFgmFnzfZtWtXQkJCuc589T3C27c9OJwm9+49evJE7/Zt3LiB27dRUABra1hZwcDAmMWKJNoN5AiFrzt1Uk9PR3o6iMDlwswMXC64XLBYRmpq94uKALwzMtL66Sd8XLCnr///fMK2bbbbth0pKpqqrn64U6dmPXsiMxOZmUhLQ3w8MjORlYXMTDx69Egi2QkgM/OHH3640rGjc926KH7Uqwdd3Q/adu0KnTFjqVRKEyb8MH/+VLnO2skJa9eiRw9cuwZTU+W8k2UQHIw//0R0NOzsZJIvrp+Zn4+ICHwyCi4fUikGD0bz5gouOqw5nD2LMWNw8mS1letjKI3Ro3HnDgYMwPHjiqQrJ02CgQG8vHDypEw1lQCoqmL7dsyYAU9PREaiCiaFzJuHV68QE1PjqhfNnTt31apV7du3l71JURGCg7FkCcRiCTAJsAL+NjBolZSEBw+wciWSk2FhAQcHuLvDyQkTJ/rev08rVjj8+ms5mqstEGZlZQEobdo9Eb1//15FRcXY2LhsPVLpSqEQTZpcVVf3Ka5PbWwMa2sYG0NNDSoqpKHRoLAwBdDS1TX+80/V2rVhagodnf8oEYv79uhx/OZNJy0tVljYxtImzTRrNlUonH3xYrt27dzWrVtS2g5eXbo0jo7eLha3MTDYM2LEAnV1JCXhzBkkJSEpCTo6qFsXVlaiEyeWFhRcBTTWru00fHiv4jlmstOvHx49gp8fLlyoxDRLURHGjEF8PG7elDXiFhVhwABIpThypEKbnAUEICfnK96fr5gbNzBokOI7bDBUNsHB6NED06Zh9WpFmg8bBnV1dOyIo0dlzd6zWFixAjwe2rTBqVOUm3tLS0urWbNmipgvj127sH8/rl5V/H608rC3t1+3bp2mpqaqqmpcXNzYsWPLbaKjg9Wr4erKbtjQLCGhIaCur88bObIhi4V375CSgrdvkZiIqCgcPgw2G2JxGnA7IgLlBsJqSI0KhcL+/fubmJhwOJx+/foJv+jNXr58mcPhmJubGxsbN2nS5Fbpu9Xt3r2bxZKwWJ2XLt376hXl55cgM3JkAJfrzuU6Ll++vmzHpFKp/GdTAtnZ2RMmzPTyGnjw4JEvX333jq5coS1bsnV12xX37lksfxOTGx06UEAAbd1K16/Tp7negweP2tp6eXr2efz48RcO09Ch1LcvKcnx/xMbGxsREfH6db6XF/XqJUd6UyCgrl2pTx+qYMYxJISaNq30PSArOzV6/z6ZmdHx45VnoQZRc3Ju8pKTQ82a0fr10idPnqSmpiqg4cQJ4nAoKkq+Vn/9JVFX725oOILD6TV8+BQF7JbN+fNkakqfXjZqVGo0Pz9/xYoVffv27dmz54IFC4iooKDg559/zsnJKa3Jx/0Inzx50q5dP1tbrwMHDpcoKRZTSgqx2ZZA1pQpknKdqYZAuHHjRgcHB4FAUFhY6OzsvH795/EpIyMjLS2NiCQSyYwZM5ycnEpTtXv3bja7ua1tu7ItZmdnFxQUVNxz5eLrO9jQcIKe3uzmzdslJYlOn6YVK+iHH8jRkbS1qV496t6dxo9P1NNrA6QDsfXrt/5SiUBAbm40b54yHZs1a6mJSXdd3Tlqak5jx2ZJyv8WfaCggHx8aOBAEokq5MDJk2RmRs+eVUiJLFRqIHz+nGrXpv37K0l9jaPmXGEVICFBpK7e3dCwL5frHhi4TgENFy4Ql0uHS74sl0xcXJy+/vDiu2Eer41y38D4eDIzo7Nn/3OwRgVCBZB3Y97585ey2bVWrowpV7IaAqG7u/umTZuKn2/durV16xKu7x8JDw+vX79+aa/u3r174MCBSvavqpBKpdHR0ceOHfuyTywW09OndOgQDRt2Sk3tt+KfioZG6zVrir4cmU9Npbp1ad26lH379sXFxVXcMR7PEZAApKq6Yu/efTK2ys+n9u1p2DASiytk/f59MjUlObehVpDKC4SpqdSoEf37Nf8u+KqvsGfOnNHV/QkgQMTj2SuWHLp1i8zM6N+ZH+UTHx/P5XYDCBDr6joIBAIFjJZIejrZ2NBff31+/HsLhEQkkUhk6QVVwxhhUlJSo0aNip83bNgwqaQCmkVFRTt27EhJSQkPD//jjz9KUyWVSgsKChITEz8esbKyUlFRcL1aFcNisTp0KLm6t4oKGjRAgwZo29Y5JmbO27et1dReWVoaPXig9uefEArh4YEOHeDnBwsL8HhYty6hW7f+bPZgHZ0dixf3nzRJ7s1qxWL88w+OH0d4ODIzVYEMgKuj89LYuKxVlR/h89GtG/tq6TkAACAASURBVOrVw+bNCu7akZOTc+fOHQODRj16mAUFybdFQM3hzZs3586ds7Rs+vPPDiNHQoZRD4YaAYvFUlP7UL1WKlWwjK2j44ciMjk5Mi3zbdSoUc+ejY8ebS2RCLW1x0yapLlhQ4V2FsvOzt62bTegcuzY0AED9EeMUFzVt4QsxZuqIRDm5uZqaWkVP9fR0cnJyflShogSExMzMzP5fD6/9BK5iYmJZ86c8fb2Lv5XVVU1JCTERZb90b8SNDQ0Tp7ctn79HjMz4wkTturp5f/5J168YJ87p3L6tMr06Srm5tS5s/jJk41S6W9SaZ+cnEmzZnlaWg62saHataVf3hKcP3/+6tU7Pj4ejo6OAHJyWNHRKidPqsbEqFhbU+fO4u3bxVlZ88eO7SgSqXh42Lu7u+fn55ftZG4uq1cvzaZNpatWCQsKFDnNxMTETp2GFxa2z8+/PGTIAj8/j/JsKgehUMhms5W1p8fTp099fX/MyhoO/Obm1mvixEFVcxY1BD6fz6ppExNlxsHBwdY26MmTngUFqSoq49+8KTA0VCQcWlri9Gl29+6aCQlxKSkhRkZ6s2ZNMC19mtmqVfMWLsxXU1OTSDT9/cWenqx9+wp5PEVMi8ViJ6dOL1+OlErFOjqdIyJO5+d/rR9HaRCRQCAo94r0KVKpVEOWOXuKdVErQr169U6fPl38PCYmxtraugzhK1euaGlp8UuZrbF79+7Bgwcr38WvhKIiuniR5s0jQ8M+wEKAgLsqKm28valOHdLUpKZNqWdPmj6dtmyhCxcoMHC7oWEv4IChYduxY894eZGeHvn5UUgIfbl0R1xeivPw4eP16rlaWbWqX//o5MkVmrAzYcIs4DhAQLKTU1fFFcmJclOj06b9DhwFCOA3bOihLLVfC191zq2Y4pvvgADy9KSK5CkfP85QVXUELrHZh5o08ZSxlVRK8+aRtTXduaOY0ccczpB/Rxz7JSYmfinzfaZGZfmNV0OPsEWLFjdv3uzUqROAmzdvtihzAY6VlVXxtBrtqqzE8JWgpgYPD3h4wNjYa+rUvVLpUUCVy9WJjISKCgQCJCTg6VMkJODKFfz1F27cOCwWbwHMsrOtoqN3BwW19/YudelF2RlmPp8/duzc9PRzAEtb23vhwnYsluKbYmRlabNYmUQA3uvp6ZQrXzPJzDRlsR4R9QAem5p+rZsGf88UL2H680+MGIEBA3DokIK10FJT7+vqtsvObiOVIiNjbX5+vu7H5cOlw2Lh99/RrBl8fLBxI3r3lsNiRgY2bDDPzIwHsgEJ8KyMbihDCSgamxUnMjKSy+WeO3fu/Pnzpqamp06dKj7epk2ba9euEdGBAwcOHjx47969s2fP+vj4dOzYsTRV33mP8COFhYVeXn24XAdr69Zubo/c3Eqecjls2GQVlX0AaWouXbJklcLm7t+n6dMT1dT6/nv7OeTLpR0yIhDQjBlkZpZZp46Hqam3pWXLhw8fKuyYvCirR5iTQ/7+VLeusFWrwaamjk2btntWBXNeaxhfdVfjM4qKqFMn8vdXsHlqaiqP5ww8Bi5qabUqfTlAydy5Q9bWNGMGyTJhm8+nwEDicsnfn7ZvP25j49awYZuTJ0teycH0CEujGnqEPj4+K1eunDt3LoDly5f7+voWH7e2ti4eOzQ2Nt60adOLFy/09fU9PT0DAgKq3smvCw0NjbNnP2xrR4QtW+Dqit9+w5Qp/xELDl6Qmjrx0aM/27Rx+fXXoLJ1CoXCT3PrEgmuXsXx4zh6FIWF6NTJ2tw89d27VQDb2PhFgwYNFHD7wQMMHYoGDfDggZGJycWcnByDmrDXopxcvYrhw9G6Ne7dU9fV3Vvd7jAogeIdKry88PvvmD9f7uY8Hu/gwVVz5841NNTX09vr4YGTJ1GrlqzN7exw8yb69sXAgdixo9ScjUiEv/7CwoVwc8Ply2jQAEDXkSPLKW585MiRj1M0vjoyMjIqSTOL6Cve6XvPnj2nTp3au5e5+nzO48cYOhQ8HrZvh7m5fG0jI8+MGDFVItFyc2uyZcu28+fZERE4fhz16sHPD926wcEBLBb4fP6uXfsADBs2SJbMz6dIpVi7FkuWYPFi+PvL554SqeBkGZEIS5Zg82aEhFR1AfQaSF5enp6eXnV7oUzev0ebNpg4ET/9VCE9wcFYswbHj8tXY08oxPjxiI3FsWOw+u+GXUQIC8OcObCyQmCgHEWLJk2alJaWJocTNQwWixUUFGQuzxVNKpVKJJLyf+MK9lFrBkxqtAyEQpo1iywsKCJCvobW1i5ABkAqKpN1dE727k3bt1OJBTdycnKy5a9YnJREnp7k7U2vXsnbVMlUJDX64AG1bEm9e1N6unKd+lr5qnNupfH8OdWqRQcPVlRPaCiZm9O5c3I3DAkhMzM6ejR13rzApUtXZWVlRUeToyM5O1NM+cvEGWRNjTKB8Bvn0iWqV4/8/UuuP/cp2dkUFkajR5OKigMgBkhNLXDnzlIX1E+btojLdeVyW0+ZIkdhm507iculwECZxj8qG8UCoVRKQUHE5VJISGU49bXyTQZCIrp3j3g8io6uqJ6zZ4nHoz175G4YEVGoouLEYu1SVQ3R0fFo1IhCQ5VfVfFbpeaOETJUJW3a4M4d/PormjfHrl3Q1Lx1+3Zs27aejf/dlfHhQxw/jpgY3LwJZ2d06IBRowYePtxdLG5uYnK2R4/oEtVmZGTs2BGZnn4VwM6d3j/88Nbe3rzEVWQZGRnHjx83NTW1t+80Zgzr7VtcuPAV78Pw6hVGjEBhIa5eRf361e0NQ+XTvDkOHMCAAYiMhL294nq8vHDuHLp2RUICFiyQo6GlZbyBgV1m5jCxGNraJ6KiUurUqfxNK74zmED47aOnh5AQHDkCP78DRUXbCwv7GBuPnD79z+fP3U6ehLo6OnSAvz+OHPm4RdS0KVO6vH792tNzfonj6llZWL++KDPzwyA+n6/j6yvMzf2wz9SnDyOjTFfXjqmpQ9XUrhCd/uWXoMOHoaT161XH/v1hf/yxrXZtsy5dFs6fbzl+PObO/Xyzb4ZvmHbtEBKC7t1x8SKsrRXX07QprlxBt25ITkZIiKxrMywtLVVU7gLvgUItrZfm5szKnEqgCjqnlQeTGpWLli27Aa8BAq7UqjVpwwYqadFtWdy4QSNHkpERDR1KXl5juNw+XG6/Hj1GElFBAT14QOHhFBREkyeTnx81aUJqaqEs1rLiVRbGxo5laC4oKIiLi1NgxLEiyJIavXfvnomJN5ACxGhqdqzKzZC/Lr7V1OhHNm4kGxt6966ievLzyc+PfHxI9mUVERGnGjf2bN7c+9KlyxU1/53BpEYZPqdBg1p378ZKpbXU1W+NGlV7/HhZGxYUYP9+bNyI7GyMHYsVK8DhANj86NEjIrK1tQWgpQVbW9ja/qfhP/9YdOt2JDubgBQDg1K/bK9evXJ37ykUOrLZsaGhwZ6ebUqTTE5OTkhIcHR0rIKFFkVFuHYNwcH3s7K6AuaAuYHBTBn3X2X49hg3Dq9fo02bcDW1IC7XeNOmRU0USvHr6ODoUfz0Ezw8cOIEMjPvaWtr29jYlNHEz8/Xz89XUccZZKDyQ3IlwvQI5SI1NdXV1c/c3LFTp8FlVGR//Pixg0OnOnWcFi9e/eQJzZhBXC516EChoYrsLBEQMN/MzL5+fbcrV66VJjNp0m8sVnFxssR69bofOkTnztG9e/T6NX3q5sGDRzmc1kZG0yws7JOSkkrTlpmZGRAwb9Cgibdv3y7bt0ePHj1//vzTIyIRXbtGy5ZRx46kp0cuLjRuXJKBgQtwTl19ffv2/WQ86++Qb75HSESJiUmamsXbosXVrduqgtpWr5ZqafU1MhrM5Xb195+uDAcZPkfGHiGzjpDhc2xtvR49Wgc0VFMbpK8/dfx4N39/WFoq31B6OiIjceIEjh79XShsTDQAeGxq+pub26GMDGRmIjMTGRlQUYGxMYyN8fy5N58fCpgAR9zd7/fuPU9XF0ZG0NX98DAwgIEBfH273LkzUCyuy+P9Eht7pFYpK5l79hx15cp7IkHfvk7jxi07exZnz+LSJVhZoX17tG8PT08Udztv374dHLzb2tp8+vSJ8q6Y/H749tYRfklkZOSAAZdychYD4PHck5PPqVdgt4i7d+96eKzMy9sNgMNxT0qKZL5dSkfGdYRMavT7Ii0tLT4+vkWLFoaGhp8eLyhAXBxu3cKtW3j6NB+wBSCVeq5Y8WzUKLeKWHz58uXu3QctLc2GDBmoqqpKhNhYnDyJEycQHw9vb3Tpglmzxvfu3S0//wCbnXTixPaWLf+jgc//EBGHDNF99CgdMGGz0zQ0dF+/Rn4+srORl4f8fOTnIzcX2dmUmfkeGA4gPb2Ho+MNLreXhga0tfHp38LCp1FROQJBOIDNm9ucOZPToYPB8OHYvh3cL+YiODo67tol86Jlhm8XJycnbe3ZOTlubPZrsdiAza7AnkmAhoaGhkZ+Xh4AaUaGsEcPlQED0KtXCd9AhsqGCYTfEefOXRgw4FeJxENdfUp09H4Wq/Ht2yh+3LkDKys4OsLJCSkpTtevzygosONwdnTqFFERi+/evXN17ZWa+ouW1uP16/2bN99+8iQMDdG1K5YsgYfHx93XeE+eXHnz5o2ZmdmXt9g6OtDRgaUl9u9f1LXr4KIivVq11I8ePVxK94PVrJnOw4fnAGsjo9MHDgwyMYFQiIKC//xNTmafOCEGAJChoeTWLZa+fkVOlOG7wMTE5MKFv1ev3s7hGN+48feQIdi7V8HC3AAaN27cpYvVqVOeRIKffhresqXWwYOYMQO2tujXD/37y10TCoBQKGSxWBXpp36nVH6SthJhxghlp6CA7Oz6AAkAATGqqlMdHMjfnzZvpthY+jSLLhaLd+3au3DhihJ3cpGLv//+W1NzRfGsUQ0Np7Vr6b9DcoqQl5dXtsDLly979Bjl6trj+PFTZYgNGTKJy/XmcNpMm7aooj4xfB9jhJ9SWEhdutAPP1S0NERmZuanX2mBgMLDadgwMjIid3cKCqI3bygkZFetWg516jiFhYWXoWrWrGUcjj2HY/f774qX1P/GYMYIvy+WLFm9b194ixZNNmxYYmRkBCA3F3fuIC4OsbGIjcXz51BTG5abOxlwZrHCxo+/s3794sr26vLlW15ei0Sig0Biw4Zjnjy5VNkW5eLly5fq6upylS5kKI3vYYzwMwQC+PmhTh1s2wY2W8nKCwsRGYmwMBw/nsbn9xSJzgIiAwOP/fuvq6lpiET4bHva1NR3s2YNyM8/DxCX2+bJkxPF14HvHGaM8Dvi0KEjK1bcy80Nj4+PuH//12bNtsXGIiUFzZvDwQEeHpgyBc2aISnpNx+fwUKhiaEhf8GCo5XtFZ+P+fOdmjXzychoY2JivGvXpsq2KC9mZmZspV/AGL4btLQQEQE/P4wZgy1blBwLNTXRowd69MDt26ne3o1ycjQBTYGgVmBglrq6mZoaPptYIxYXSCTGAAtg8flGmZkFTCCUHSYQfgtcvnwvN7c3oCeVDnj9es3MmZg7F40bf179pFGjRi9f3s7MzDQ2Nq5sl3Jy0KULbG0RFTWRzZ5Y2eYYGKoFbW0cP46uXeHvjy1bUGKVwQpiZ9fE3DyhsHC5ioqgWTPphQtmJYoR1fXxUY+LG0kk0dU18fGptWEDOnVSvj/fJEwg/Oq5cAFHj3ZQU5snEhno6JwYPrzD0KFlyVdBFMzKQufOcHDA+vWVcmlgYKg5FMfCLl0wdixCQpT/hVdVVb19Oyos7JC6unrv3rNLE2OxWFFRf9+4cUNFRcXJyen8eYwdCzs7rFsHHk/JLn17MGOEXzHv3mH6dJw7h1WrYGJybufOcBcX27FjR6gqPI9NGaSlwccHnTph+fJq9EImKrgfIcOnfIdjhJ+Sm4tOneDiguDg6nblXwQCLF+OjRuxaBHGjPlOb0llHCNkBki+SqRS7NoFe3sYGeHxY/Trh/btvXbuXD1x4ujqjYLv3sHbG76+X0EUZGBQIvr6OH0a164hIED5yqOjz9rZdXR27nrjxk3ZW2lpYcECREdj+3a0a4f4eOU79s3ABMKvj9hYuLlh2zacOYPg4M/HzKuRV6/g4YH+/REYWN2uMDBUOQYGiIzEpUuYPl2ZarOzs4cOnX7v3s5bt9b07OkvFovlat6iBa5cwZAhaNsWCxZAKFSmb98MTCD8msjOxpQp6NYNEybg/PnPK1xXLy9eoH17TJyIuXOr2xUGhmrC0BBRUYiJwfjxzydOnP3bb8uysrIqqPPFixciUUvAAqhfVFQnLS1NXg1sNvz9EReHBw/QogWWL49u1MijYcM2x4+frqBv3wxMIKzpzJq1zMyspZWV8+zZ0U2aICsL9+5h+HBFMv7Z2dldugyrU8dp6NCfioqKlOjkkyfw9ERAAH7+WYlaGRi+PoyNcfRozrZt/TdscFu+3NTLq18FFd6+3Tg3N47FCmOzd2dnvw8NNVdsXoeFBcLCMG9ewezZ058+PZqQEP7jj7/l5uZW0L1vA2bWaI3m/v37W7dee//+NpCzcmX7q1fjHCtQ83LKlAVRUb4Sye5Dh5Y0bbph9uyKRq13794RUXq6eefOWL4cZc9WZWD4TkhKuqOj45Wd7ScW482b3fn5+YpV087JwZgxePZM8/z5Y0ePbtHUVO/W7dikSaxz57B1q4IlSd3cUjmcBmlpJgAKCpokJaXY2THVBZkeYc0mLi41N7cxwAaMjIzUW7QQVUTbnTvPJRJvAIWFHe7cSaigb+PGzWzRYoit7TBX11+CgpgoyMDwgYYNG6qpXQbeAfcyM/ODgnQFArmV3LwJBweYmuLqVbRpU2vlygWLF89u1Ypz5QqaNkXLloiKUsQ3KysrI6M3amqbVFU3A8+6dWuweTPkHHb8BmECYQ3l9WtMmYKAADdt7fOamksNDSd4eLRQeKK/QIAFC5CY2EtTcyJwSE1t5r17vZ8/V9y99PT0w4evpaefycqKUVd/4Ob2RnFdDAzfFubm5tu3z7ezG+7puTA6esfTp2jYELt2QcaUJhGCg9G1K5Yuxdq10ND4z6tqali2DAcPYvx4jB2LggL5fGOz2devn1yypGjJEsGrVyfDwlQOHEDz5jh4UFb3vk0qu+ZppfJNFt1OT6cZM8jEhGbMoMxMysvL279/f2RkpFQqVUzhsWNUty7160fJyRQZGTV37rKrV68FBxOXS3v2KOhkbOw7DY32xdW0udzOL1++VFBRtVJYWChLQV4GWfjeim7LxbVr1Lo1OTvT5cvlSGZnU9++5Oxcfnn67GwaMoRsbenu3Yq6Fx1Njo7k4kIxMRVVVdOQseg2EwhrEBkZNH8+cTjk708pKUpQ+OwZ+flRgwZ0qqRtGB48oObNqV8/ysqSQ2dREQUFEY9H9vYBHE4nLrfzkCGTlOBrdcAEQiXCBMKykUopNJSsrKhfP3rxomSZ69epXj2aPJmEQlnV7txJXC4FBZGi98n/ca9BA+rQgW7dqpCqGoWMgZBJjVYbRHTp0qWoqCixWMznY/lyNG6MxETcuIGQEEW2IvuU4lyoiwscHXH/Pnx9S5CxtcWNGzA3R8uWuHJFJrUxMbCzQ0wMrl1DXNyq27c337y5cc+etRXylYHhO4DFQr9+ePgQTZvC0REzZ37YPuL9+/e5ubnF6dCePbFmDYKDIft+gsOH48YNHDwIX1/Ex2etWBG8atU6BeaCfnSvXz/06IH+/bFuXYSVlYulpdO2bd9B6a7KD8mVyFfdI+zRY4Sx8Y8GBgGWlt7m5qJ+/ejpU8W1FRYW3r17NzMzk4jCw6luXfLzIxkTlocPE49H8+eXtbNaQgL5+VHDhnTypOJO1jSYHqESYXqEspOcTMOGUe3a5OY2ictta2Li0qDBKk9Pev1aQYUiEc2dK1JRaaWiEqKmtr5xY3eFR1KIiM+nhQsL2Ww7IAco4HJbp6WlKaytemFSozWanJwcDsezeIxNXX10aGiFkhFpaWl16zpzuSNMTBxdXCIbNqTTp+XT8OoVeXpS+/b05s3nL/H5NH8+mZjQ/PlUWFgRN2scTCBUIkwglJfdux+pqvYGCJBoazvl5RVURNvjx4+NjIb+O2zfr4K7ar99+5bL7VKsTVNzxJ49FR6HrCZkDITMOsJq4OFDrFmjmZGRAwgBdQODZAcHw4oo3LRp58uXE6XSH4B0VdUByck+8k4vtbTE2bNYtAhOTli3jn///obnz9/89NOQxETnadPQrh3u369otpaBgeFTWrQQGRpqvX8PgK2jowZIKqLNwsJCTe0xkAlIMzISLl0yq1tXcW1mZmZWVhI+f6FYrGVomLB4cdPVqzF5MgYNwjdZo54ZI6w6pFLExKBbN3ToACMj9eXLp/F4zlxuy9Gj29WvX78imp8+VSEqrhQj4vHYin1TVVSwYAEOHMCQIWMWLdLYvbunu/tPCxc+378fu3YxUZCBQck0b97c3l7E5Q7kcjv37dtOsUX3H9HX19+6dVGDBt0bN+6zZs3KtWu1vLzw5IniCv/559jmzTabNvGePYt8+FA1MBAHD8LKCgsWID29Ip7WSCq/b1qJfC2p0exsCgoia2tydKSQEBII/v+SWCyuiOaEBOrcmRo1yrK0dDM17W1m1vL8+YsV9NbU1Kk4JcJmr922bUcFtdVkmNSoEmFSowoQHx/fo8ewYcPGv337VrmaRSIKCiIuV8kjGk+e0OTJZGREw4b9f9nG/fv3b968WZFRycqDmTVabTx//tzZuYulpdMvvyyIj8eUKahbF5cvY+9e3LoFf39oav5fWOWzXeRlRiDAvHlo3RpeXrh/3zAp6eL166uSkq60betRQf8tLbksVjSQZmQU4eBgV0FtDAwMJZKXl9e+/eBjx4bv3evdtm0f5SpXVcWUKYiLw/37aNECZ88qR23DhggORlISHB3RrRvatEG7dlO8vOZ37ry+bdteUqlUOWaqniqIyZVHzewROjl1Bm4BEhWVHzmc40uWkNKnXMk7L1QuHj58WLeuq75+4xkz5ilfe02C6REqEaZHKC9Xr141Np5UnH0xNfVNT0+vJEPh4VSnDg0bRsq1UFREf/2Vq6bmXnwKxsbD4+LilGlAGTA9wmpAKEREBB49SgccALZU6jlt2vPZsxUsj1siiYno1g2//IKNGxERgTp1lKb5I+PHz331anxu7t/btp2NZ3bzZGCoHAwMGuTmXgOeAreys9OlUpNKMtStGx49goUFmjbF5s0QicTh4eHHjh0TiSpUu1hNDQMHqhkYFBRP88nOzvrhB+1Vq5CcrCS/qxAmECqBwkJERGD4cJibY/lyODh01NGZAOzhcoN79+5SQeWvX7+eNy9w/frNmZmCBQvg7PxhjXynTkrxvQQSEl5JJMMBu6ysAVevXqssMwwM3zHx8eja1WTEiNXu7jN9fIJGjNjdqhXr3r3KMqejg8BAREZi61YyMek+ePCVoUOvt27dlSpWYFRTU3PatNEcjjOX23r06Kbr1zd8+RJOTrC1xfLlSElRlvuVT+X3TSuRKk6NSiSSoKBNvr7Dg4NDpFKpQEDh4TRsGBkZkbs7BQV9WIQnlUoPHTq8bNmqhISEClrMzMy0sGjJYu1TU1uuqdnLz6/U4kzKIjOTjIx6sFihwBMVlbYXLjyoXHvVCpMaVSJMalR2Ll4kMzPates/B/fvJy6XjhypXNPPn7/U0+tZnMzkcAY8rUgVj38RCAR5eXkf/xWL6dIlmjyZTE0/XBhTUujKlWvNmrWvX7/1rl0HKm5Rdmr0gvo3b97MnDlzzJgxYWFhX74qEAjCwsICAgLGjx+/Y8cOkUhUmp4qDoSrVm3U1R0PPNDSGufktMnQ8P8fcyVx+vRpXd05xd9aQ8PWlX3Vjo+nxo1p+PD0H34IaNeu/4gRJxs1KmGJ/TcDEwiVCBMIZSQsjLjcksv/xsWRlRXNmFHRwqFlkJOTw+U6AYWAkM12Xbjw/aeT2JWLUEgRETR0KBkZSTU07IFXQB6H0+bVq1eVZfILam4g5PP5VlZWEyZM2LFjh5WVVUhIyGcCJ06c8PDw+OOPP0JCQmxtbYcMGVKaqioLhHl5FBVF9eoNBp4ABDxp0mRwZVcdeviQfH0T2GwvIB9IsLJyrlRzx48Tj0d//fWfgytWUMOG32wsZAKhEmECoSwEBVHt2lTGnJKUFHJxof79ic+vLB9CQnaZm9ubmdnPnbu9Xz+ytKSQECq9u6EE0tPzDAzaFt/Qq6hM9fO7sGtXpSe3iqm5gXD79u3Ozh+u6eHh4fXr1/9sAcqnXcA7d+6oqKgUFJRcfEiJgfDMmbMDB05cuHDlR1upqRQeTjNmkLs7aWuToyO5uQVpaPwCJOvoTF22LFgpdkvkxQvy9ycejwIDacOGPfXrt3Zw6BQbW1kzsqRSCgwkS0u6fr2EV//4gxo2VLwKYk2GCYRKhAmEZSMW06RJ1KwZldsdKiyk4cPJ3r5S5oR/yY0b5O1NjRrRzp1lVRuuIM7Ovhoaq1RUdpqbO65ZkzdsGFlZkbk59etHQUF069aHTvCLFy969x7Ttm2/c+cuKMVuzQ2EP/744/Tp04uf5+XlAShjMenZs2f19fUlpXw+ygqEsbGxJiaewBV19RWurv6TJ5OjI+npUYcONH8+RUd/WAIvFosXLlzZunXPRYv+rOBC+NJITqbJkz9sRpidXRkWPicvj3r3Jnd3KmNF78qV1KABJSdXhT9VCRMIlQgTCMtAIKB+/cjbW44fdXHfscR708ogOpqcnKh5cwoNrRT9fD5/zZqNv/++4tOr/Zs3FBpK/v7UtCnxeOTnRxyOO4t1Bojn8Vq9UUYmSsZAyKIq35bYz8+vbdu206ZNK/5XW1v78uXLLVu2/FJSIBC4uroOGTJk+vTpJaoKDAxcu3atq6vrxyMz8pDsKwAAIABJREFUZ860tbX9VOb06Zh584K0tbWDg3+zs2vx8Xh+Puv5c9azZ6zERHZ4+LrYWC4wFICWlvOiRZfbtJE2bSpVdLG7ImRksIKCVHfsUB04UDJ9uojLrYrP5flz1oABGq1bS//8s6jsnV/WrFHdtEnt1KlCK6tvZx9roVDIZrPVvsniiVVOfn5+BYuEfatkZrL699eoV0+6fn2RXN+1yEiVsWPVly4VDR4srjTv/sO5cyozZ6rp6eH334vs7fP27j0gEomHDh1gYGBQ2aZTUljnz4snT/YVCC4DYLPnODu3cXXtUL8+2dhQ/frSWrWIxfogHB8fP27c3Kys7GnTxgwd2r8MtZmZmbq6uoaG5RRzroai2xoaGkVFxYUxQUQikUjz01Ir/1JUVNSnT5/mzZv/+uuvpakyNjbOymqWkJAxZ854FosFoH79+p9qS0tLGzt2yfv3x4AcP79BP/98qzj4PXuGvDzY2MDGhmxs4O3t8PTpqvz8Dmz2pdat6wYEVMXbEhNzftKkeUVFoqlTJ2dkDN6wgdWrF927R+bmbECjChyIimINH86aM4d++okNlPARfMr06dDQQJcuWjEx0ooU861RsFgsJhAqi9J+yN8tz54909bWFggs/PzYPXvS0qVsFku+96dHD9SrR716qR85cvrp0+Vqampr1y7w8vKsJIcBdO6MTp1w+DBr7FjN9PTOQmE3Is2QkF4PH16o7J9JvXqoV08jOFjj0aMwqbSWkVGUv//ktDTVu3dx6BDr6VNkZ6NBA9jYUIMG2L59Ynr6VqDOrFk93dxaNW3aqESdGzbsWLhw8/79+7y9a14grFWr1uvXr4ufv337ViKRWFhYfCYjEon69++vra29Y8cONrvUxY7a2tpSaVhCwq/79tXW0nLn87FpE3JzkZcHPh98PrKzk1isVoAZYMbnczMzsz08jEeOhI0NatUq1lF8j9G2QYOXISE/Nm5cLzh4YxkWlYVEIhk+PCA19QygHRDQaeDAdrGxtSwtWf/6U1k8e/bs/v37zs4ue/fWWrsWR47A3V1WowEBYLPRqRP77NlKWchf9bD/pbod+RZg3smPEFGnToPu3i0qKsoUi33XrJk5cqSCP207O5w+ndms2WKRKAoQDhnSKTn5ZqXGJDYb/fvDxeVtixa6ItFUAJmZt+Pj4+3sqqLaYkzM/t9/X52eHjN79uaWLU0/fSk/HwkJSEhgPX0qzcuTArYAsrLa2tsnqKo20dGBoSF0daGri4/PDxxYV1BwLT5exdu7HLvVEAh79eo1bNiwFStW6Onp7du3r3379sX97qtXrxoaGjZp0kQikQwfPryoqOjAgQOqquV6SGw2q1EjcnCAjg50dGBg8OG90NGBmlrzJk0mvXlzREUl28ZGuHq1cWlaxowZPmbMcKWeaFncvp2bl8cDOAAMDOx++inZ0rJWZRs9duzE6NHLCgo6SySLGzTYev16y1py2pwyBWw2vLxw9iysrCrHSwaGr5zY2Ni4OM337/8GyMDAeciQqYDMW85/gUCQbGjYIj3dEIBYXCc9Pf3LnoPSsbAw0dJKzsvLANRzcx+MHl37l1/Qp0+l78Fkamq6YUNgiS/p6qJlS7RsCYAdHV3r+vVVRUV1zMyO3b07Xk8PfD5ycj50gfLzPzw/dIgAlkSW7a0qPhopL1KpdMCAAQ0bNuzevTuPx7v+73Cwr6/v3LlziSg0NBRAs2bNHP+ltHUnu3fv1tfv7uXVp7TZNET04sWLn3+eO29eYPHu7dWLRELR0eTnR6amZG7uq6GxVFV1Xb16LoLKW8vzCa1b9wReAARcHDRoisJ61q2junUpKUl5nlUTzGQZJcJMlvnI7du3DQyK98iVcLkthUJhRbQJhcK6dZ1VVLapqq5TVfVcs6YSVxl+yqlT0TY2bvXque7bd7j4qmVuTvPnK7lgqWLw+fyVK9dOm/b78+fPyxALCtrM5baOiSl/j+JqmCxTHH3j4uJSU1NdXV2NjIyKD75+/VpTU5PD4eTn56elpX0qb2lpWWI2YM+ePYcOHTpy5EhVOF0xcnPx998ICoK2Nvz9MWwYWKzCvXv/LigQDB068OObUKm0bPnjnTvDgHaqqjsDAl6uWDFPYVUbN2Lu3J0s1moVFcyZ89NPP/2oRD+rDGayjBLJy8vT09Orbi9qBAkJ1Lz5CC2t92pq2RMn9pk/f2oFFWZkZISE7NTUVO/Y8YcxY/R4PGzfDg5HKc7KwdOnWL8eu3fDzw+//ooWLcpvUu2kpqZqamqWO9mn5ED4/v37Cxcu3Lt37/3796qqqjwer1WrVm5ubtra2pXjrYLs2bPn1KlTe/furW5HPicyMvLZs6Ru3f7H3n2HNXX1cQD/JkzZGwRBBEURJ26GExUUtU4QFEcVt6j1de+Jq4qtG60I1oETF05AQEDFDThQhou9R4Ak5/0j1loNQwi5Ac7n4ekTkpN7v6lJftx7zxhkZGT09d1jZwd3d9jZMZAnOxvu7oiNTeZwxuXnFxsba928ebImPcHy8vIMDPoVFIQDLG3tnnFxlzU1a2vK4NpDC6EI0UIokJkJKyv89hvs7ZMVFRVF/rkoK8PKlfj7b/j4oE8f0W67SnJzcfQofv8dhobw8ICeXuTs2Ss4nJLNmxeNGDGEgUAV4vP5PB6v8s/4d0eIN27cGDp06NdF8hQVFeXk5L7enjJlyvPnz6t/QCtqkrkM04IFq1VVJ7LZB9XUOvTtm6CnRxYvZnIQXlAQMTIiCxaIcn3OpKQkXd0vMxZqa4tmxkLxo6dGRYieGiWEcDjE1pb8M0y6Ft28SQwMyPLltTsjTAVKS8mJE6R7d56UVHsgGcjS1u6ekpLCTJry/fQyTImJiQMGDHBwcOBwOHv27Hn69CmHwykoKOBwOHl5eWFhYatWrYqOju7QocPUqVMLCwtrp37XB8ePX8nNPcLnT83Jmaevfzk5GZ6eaNKEgSRcLtasgYsL9u3Djh2QE92gDCMjo2bNSlVUFsnKLikryzExaS6yTVNU3UQIpk6Fri42b671fdnZ4fFjPH0Ka2u8fVvru/uRjAycnXHlSo6qqi5gCKjn5nY6eDChjlaGfwthVFSUqalpQkLC9evXp02b1q5du6/HgsrKytbW1osWLXr06FFUVFRKSsrHjx8ZCiy5njzB8uVo2RI5ORpADECUle+PHGnI1Im3hAT06oXISDx6hEE1XQxKiNDQi8eP9zx1yrpdu8sbN9buqA+KknyrVuHlS/j4QDyjSLS1ERAAFxf06IETJ8Sxxx9paGg0aVImI7OPzfZTUYmKiGinr48hQ3DsGOpYRRTDwWntEf+p0aKiopiYmG/nPn3xgqxeTVq1Ik2bkrlzSWgoiYt72bZtv8aNLadNW8QXTwevH/j4EC0tsmuXODqYpaQQQ0Ny9myt70jk6KlREWrgp0b/+os0a0YYOTUYHU3MzMj48SQnp+zKlSuBgYG1NAGkULm5uZs371y5crNg7rSsLOLjQxwdiYoKGT2aBAQQQZ9ZPp9//fr1s2fPiqeH/FeSO8WaCIm5s0xcXJydnSuX21Za+vmBA38/fNjq5EmUlGDoUIweDWtrsJg7LiotLX358qWRkRGLpTZzJp4+xd9/i69bV1QUHB0REoLWrcW0R5GgnWVEqCF3lgkJgbMzgoPRUvgMJ7WuuBhLlpCDBwdJS7eRkuK2afM2LCyAmSj/yMzElSvw90dkJBwc8PKl25s3KmVlWoaGt548uS0nwus0FapiZ5nKB9SHh4fv27cvNTVV8OvNmzdFkK5uWrJk16dPfwDWQPjYsTvnzj1w8iQ6dGA6FpCent6t2+CCAouysmeysptdXQdER4vyimClunWDpydGjUJUFBrqlyHVQMXGYswYnD7NWBUE0KgR5sx5e/SoSl7eNgBv3vySnJxsxOj8T5qacHODmxs+fICfX9HJkwllZaEAPnzIunz54ciR1gxm+1ElhZDL5U6dOtXT09PgZ+cgqYPS0tLU1dW/+9shIwNhYbh7F6GhePKEBfABADx7e9bGjYzEFOLAAZ+kpFl8/gQgvWlTp99/HyD+DL/+iocP4eaGc+eYPDKmKHFKTYWjI3bsQK9eDCdRU1OTk3sPcAGSlfUpOVlNQuZBbNIEixbJ79yZl5ZWACiUlb2eMUN7wQL07ImePWFri1at/tM+OTm5oKDA3NycJcbvkUoKYXFxcefOnYcOHSqeNEzhcDi2tr+8f88CPh87tqNt235hYQgLQ3g43rxB166wtsbmzdDVXeDg4FxW1kpW9uWmTaeYTv2vhAQpQsoAAFwVFcbme/TyQt++2LABK1cyFYGixKeoCEOGYOJEjBvHdBRAS0tryZIJW7ZYslisAQPmjxql4uaGNWsgCWO/2Wz2n3+umzevJ5fLnzFj/Jo1Zp8+ITwct25h1y5kZKBrV9jYwM4OZ854Hj58E9Bu2bI4OPiclNgWAKr0KuLcuXMTxbOW8M+rSmcZHo8XFhb26NGj8huQrVuPyspuAgiQISNjratLRo0iXl7k8WPy3VXnkpKSV69eSU4Pi8+fyfjxxNg4q0kTK13dkbq6HYOD7zKYJyWFGBmRS5cYjPATaGcZEWponWV4PDJiBHF2FtOEZz8rM5O4u5NmzUhgINNRKvP+PfHzI9OmkVatylisjgAPIEpK0w8fDq3gPfX48eObN29WOn2dCBbmNTExUVdXF1wAV/9HpVsUJ19fX1PT7vfv3y+vAZfL7dZtkKbmNC0tp7FjZ/J4JDGR3LpF9u0jv/1Ghg0jrVsTeXmipnaAxfoDIEChgUE3cb6EauPziY8P0dMjc+eS/HzC5XITEhK+7c7KlIgIoqVFYmOZzlEFtBCKUMMphHl5eSUlJfPnk549RTlJRW0ICiJmZmT0aJKaynSUKuDxeFpa7YEygMjKTmzRIlJBgTRuTHr1IlOmkC1byLlz5MULwuGQmTOXamqOUFNbYG5uU0E31PT0dHf3/z16VPl0HxWdGn3LyEDNn5SSsmHw4Fm3bgUoKurl5KCg4MsCTNnZKCrCmzdPnz3TLS7eD8Df3/b8+XwtLeXmzQXrWsHKCoLbHM7oLl3sc3Ji2exnGzZ4MP2aKvfsGaZPh5QUbt3CP+sQSxkbGzMa6ovu3bFjB4YPR1QUan85T4oSKyen6UFBLzicAkXF32Jjx4uzS1o19O6NJ0+wZQvatsX69Zg6VaKv37PZ7KVLZ3h6WrFYKt26Nbl4sSuADx8Eqy/hzRuEheH1ayQl8cvKrvN4DwFWUdHSyZND2rcfqKb2ZcUhZWWoqn65PWLEhKdP3dzcKl+s4yeGTwiOMRs1alSj1ypSfn5+EyaM5fPXqKtbqas7qKt/ef1KShD8fykpeX3s2G+FhZcArrp6l3fvItXUhL9zS0pKnjx5YmRk1LhxYzG/ip9SWIjVq+Hnh02bMGmS5L6tp09HejrOnJHchKDDJ0SqIQyfePDggYPD7sxMX6BUU7NrWtqjurIE4+PHmDoVamrYvx/Nm6OwsJAQoqSkxHQuIfLy8jgcjo6OTnkNuFzo63dMTw8DFBUU3FxdZ2ho9MjO/nIIVFCAnJwvt9+/70LIg7Fj+X//Xdk/07eHh926dfvzzz8Ft/l8vouLS0RExNdHBcP1qnfMW0t8fX0VFcO0tbt+/PixvDazZi3T0rLU1m6/e7e3OLOJCpfLjYiIePnyJSHk4kXStCmZMIGkpTEdqzIlJcTKimzYwHSOCtFToyLUEE6NBgcHKyvPAQjA19buVLfePFwu2bGDaGmR3r03aWl10dbuumDBWqZDVdPJk+d1ddtpa3dxc5tbQTNLy4Es1oWIiMqH8P+nsBkZGW3atElwm8fjAThx4sTXRyWzEBobd7t7N7TiZnX3+660tLRjRztNzakaGoNMTVe3akWCgpjOVGWfPpEmTciVK0znKF/dfWNIoIZQCJ8/L5WR6a+sPEVLa7CHxyqm41THo0cZsrLWAB8g2tq9BdPB1EV8Pp9T2RXaz58/jx8/5+HDV5VujYEV6kXLysrU1tam4jZim8VA5MLDwxMSWufkeAGkpKRLauoKRcU680/WuDFOn0a/fksVFAKVlWV9fHb07FnJvxRFSayPHzFkiMzBg4Hm5g/U1NRaMjh+vgZ0dUvU1JTT0lgA8vJUnj0r1tNjOlO1sFisSr/Y9fT0jh7dxavCEvV15lu1ASorw/XrCnl52QCAUkVFroKCuEbViAiXG8pifczMfJSZmebmNjQxMYrpRBRVHTk5cHDAjBmYOJENdGM6TvXp6+tbW+uHhjoDbB0dVTe3Zg4OWLsWEjL6nil140pvQ8Pn4/hxmJvj8eOu/fopaGvbamt327BhoTinWhCJ1NRUQswBFqDL4XCZjkNR1cHhYNgw9OuHRYuYjiIK584dvn172a1bi2Nijr1+DUNDWFrit9+Qmcl0MubQI0KJc+sWFi2CnBwOHkTfvgD25+XlycvLy8rKMh3tp9nZ2Wlrb05NZfN4sTo6/ZiOQ1E/jcfDuHHQ18eOHUxHEZ12/8zHr6KCdeswdy62b0fLlpgyBcuWQUWF2XQM+P6IcPny5VJSUlJSUoI+5S4uLlL/GCcJ8wjVa+Hh6NUL8+Zh6VJERAiqIACoqKjUxSoIQE1N7dmzO0ePNg0I+JWQrXv3Mh2Ion7SvHnIzsbRo2JaZZARWlrw9MSjR8jOhpkZtmxBSQkAZGdnZzaM48T/HBGOGTMmIyODqSgNTXZ2tpfXoYKCIg+PX1NSDNetw4sXWLoUv/4KsU2wJwaqqqrOzs4AzM1hZQUjIzg6Mp2Joqpm9WpERuLOHbGu5cIUIyMcOIC5c7F2LVq2hLn55ujoS4D00KFdvb23M52udv2nEG7bto2pHA2Qre3wV6/G83gme/cO19IKWbVK8fx5SNffc9XGxjhzBsOH486dr7PhUJTk2r8fJ04gLKxhrSxmYYHTp3H1at6wYRe53AiAFRAwOCEhoVmzZkxHq0X192hfsmVlZX3+LMfl/krIKDbbytf32ZQp9bkKClhZYdcuDBmCf1a3pCgJde4cNmxAYCDKn+GkPuvRg6euLguwAGRlyf/1V1lBAdOZatO/hTA+Pj49Pb0qz4mLi8vOzq61SPUch4MjR9Cvn1peXhoQD6QpKt43NzdlOpeYjB2LCRPwyy8oLmY6CkWVIyQEM2bg8mWYmDAdhSHq6uqOjp21tR21tUd376706pWZsTEWLkRiItPJase/hfDp06cmJib/+9//YmJihDYlhISEhLi4uHTo0CE3N1dcCeuPz5+xZg2aNsWJE1i3jn3vnneXLh5t2ow9cmRtBRPr1T+rVqFFC7i5ocrT3FKU+Lx4AScn+PmhQwemozDqyJHf793bFRKyLizM59QpPH8OJSV07YohQ3DrFtPhRO0/k24HBgYuWrTo+fPnrVq16t69u5mZmYaGBpfLzcrKevr0aURExKdPn+zt7bdt29amTRsGQ3/l5+d37do1wdxvkqa4uPjrBOXR0fDywpUrGDUKHh5o3ZrZaMwrLcXAgbC1xbp1TMagk26LUF2fdDskJNTT85CSkkZk5JLff9cbPZrpQBKJw8Hp09i6FbKymD4dbm5gs0v9/c8UFRU7OY1WkbyBF3w+n8fjVfoZ/371CUJIcHCwj49PUFBQcnLyl0YsloWFhZ2d3ZQpUywkqZ+DZBbC0NBwJ6dZXK5i69aNp08/4eUlk5GBKVPg7g51dabDSYzMTPTogRUr4ObGWAZaCEWoThfC5OTkzp2d0tP3A0n6+rs/fqx3hzwiRQhu34aXFx4+BJvtmJ3dhcdTMTI6/eJFsKTNZ1nFQvh99wwWi9WnT58+ffoAyM3NTU1NlZOT09bWVlBQqK2k9c7UqUs/fw4E9O7eXZudfdrT03XgwPo8CKl6NDUREIDevWFsjJ49mU5DNWxPnjwpKhoEtAfac7kb+Xx+XVlfiREsFuzsYGeHyMjsvn1LiotXA0hLe/Pw4VNr665Mp6uOiv6xVVVVzczMmjZtSqtgFcXHY/NmJCSUAOoA2GzdqVPzHRxoFRSuVSscPw4nJ8THMx2FathatuxYWnoFeMJiXWzcWIVWwSrq3FlZWTkNyAY4xcUvhgzRnzwZ16+DW9emU6T/3iKQlAQvL9jYoEcPxMVh4sSpWlqD1dQWGxh4OzuPYjqdROvXDxs3wt4edCIHiimlpZg3z9DObru9/c6JE0OuX/djOlGdIS0tfeTItmbNHAwNbf/4Y8rz503at8f69dDTg5sbLl2qMxXxJ1aol0BivkZYUlIye/by4OCI3r17/PnnxpQUuQsX4O+PV6/g4IDRo+Hg8GUs4Js3b5KSknr06KGoqCiebHXaggWco0cny8m9NjDQCgg4oq+vL5790muEIlRHrxGWlWHECMjL4+TJejWjE7Pev8e5c/9+N7q5IT39/G+/reVyibu7y4YNi8WWpJqdZeoWMRfC5cs379jBKilZIC29Q0uLxWItGTkSo0fDxoae/KyRjRt3rlpF+PwFLNatIUP+vnjxiHj2SwuhCNXFQsjlwskJAE6dqv/TWTAiPh6nTuHkSV5cnCWPdw9opKnpeO/eLjMzM/EEqGIhpN/flSMEMTHYtw9HjsSUlAwHZLncEQYGzz98wB9/oGdPWgVrKjHxM5/fAQAhHZKSPjIdh2oQ+HxMnAgOBydO0CpYW5o3x/LliIoqUVVVBhQBdlaWqbNz+po1uHMHRUVM5/sH/Qr/4vBh386dB48fP1cw2zqPh5gYHDyIMWOgo4P+/REUhIEDHZSVVwJ3VFVXzJ49mNY/UZk2zUlLa7GU1B5paVdT0wlMx6HqP0IwcyY+f8aZM6ibK7vUJQoKCjY2zdXVpysprWzV6tmmTV0AbN4MXV107gwPD/j7IysLAB49emxpObB5cysfn5PiTEhPjQJAaGjosGHbs7P3sVihzZqdNzE5ef8+TExga/vlR0/vS8vz5wOuXg0bNMhm+PChNU1PfePt27d37941NGw/a5bl9OmYP18cO6WnRkWobp0anTsXjx7h+nXQi/jiQQgJCgrKz8+3t7f/OtawuBj37yMkBHfv4sEDNG+ON2+65OefAnS0tBzv3/+r5jN9N5RrhH///felS5ekyr/MffCg78aNf8jKyuzZs37AgC9L/JWUID4eb958+e+dO17x8ZrAOADKyl1OnHhgYwNVVTG9Cupb79+jVy8sWQJ391rfFy2EIlSHCuGSJbh9G7du0c+4BCkrQ0RE6eDBfQoKwgEAy0xNe1taDmjRAmZmMDNDixbQ0vrS+O7d8MmT/1dczJk2bdyqVQsq2GxcXJyurq6GhkbFe2esEGZlZRUVFTVp0qQmG/Hz85s5M8jUNOn+/Ws/fp3l5SE+/kP//mOzsm4AhcrKdi4u0W/fSr15g9RUNG2KFi2+/ACPly2bn5OzRU4ubMiQ1/7+B2qSiqqht2/Rpw82bUJtLwVNC6EI1ZVCuG4dzp5FUBAq+26kGGBtPTQ6uldpqb6e3ra//76TkqL2+jVev8abN3jzBsCXb+zLlzvn5l4GtNTURvj5rbeyaq+oKOQU94oVWw4cCD558kC/fkYV75eBa8SEkDlz5pw4cUJZWVlXV/fy5cva2trfNRgxYsTDhw8/fPgQFhZmbW1dwdZKSnbHxi63tQ2Xl+9dUIDcXBQWorAQeXlQVoas7MecnHZAI6ARj6dtapo3YoR68+Zo2vS7rtIdzcxW7d9/tE0bk0WLdtXKy6aqzNQU16+jXz8oKGDECKbTUPXI9u04cQIhIbQKSqibN08eOuSTkZE4bVpAkyZq3z2akYHXr/HqFf/CBTagByAvr+OUKe9LStoLVolSUoK6OhQVoagIJSWEhZ3icB4WF1d+sPdvIXz16tWzZ89G1/5cszdu3AgICHj9+rWmpubYsWPXr1+/e/fubxuwWKyhQ4euW7euV69elW6NzebKyGQ6OSl07AhlZaiqQkEBiopfTnpwOO0tLB6+f/+HtHRuu3Zy//tfuXN99uvXt1+/vjV+cZRomJvj6lXY20NBAfb2TKeh6rj09HRlZeXDh+X370dISANdYrBOUFBQ8PCYUd6jWlrQ0oKVFfviRYs7dxZwOM0aN77y/Plvgrm+S0tRUICcHBQUoLAQBQWIjgaHU1yl6U/JP44cOdK6dWvBbV1d3fDwcFI73NzcFi1aJLgdGhqqoaFRXksNDY2wsLAKNuXr66us3N3FZVYFbbKzs//8c//Ro8c4HE71AlNMiYwkOjokKKi2ts/hcEpLS2tr6w1MXl4e0xGE4PF4/fs76ejYKSt30NQ89e4d04EoUeDxeOfPnz940DszM7OCZidPntPWbnfr1qtKN/jvEaG6unpWVpYYZptNTEy0sbER3G7evHlWVla1ry7weLwePVQWLJgUHR0NQFpaunXr1t9d8lFTU5s1a1rNY1Pi160bTp2CkxMCAtCtG9NpqDro9u3bDx7o5OScBIp1dKybNRvDdCJKBNhs9i+//FJpMyen4cOHD+bxeJW2/LcQdu7cOScnZ8CAAaampnl5eVu3btXV1f3xCQcO1LQjSX5+/teF+gQzkOXm5lavECYmJj548GDKlClf79mxY0eXLl1qmJCSHJ06Yd8+6aFD5c+dK2rXji/ajdPOMiJUILhKI2FycnI4HCUAgAxA8vPzGQ5EiRefz6/KqdF/C2GTJk38/f23bt166dKlkpKSO3fuSAubbqHmhVBXVzcnJ0dwOysri8ViVXt9dlNTUwcHB0lbj5ASrREjwGZj9GjFW7dEvKaxrKwsLYQiJIG9RktLh/J4+9TUpsjIvJs3b5IEJqRqlWAcYaXN/lPqHB0dHR0dAejp6Z07d87Kyqo2krVv3z4qKmr27NkAoqKiWrduLUundqAq9MsvKCiAvT1//vzTiYkxo0fb29hU1JeYogCcOoVFi+QePrxZWvpYR0fHyKgrLKqDAAAgAElEQVSSPvRUgyX8cuChQ4datmxZwdPi4+P3799fvV1OnTr14sWLfn5+ERERK1euFFREAKNGjTp58su0OtevX/f39y8tLQ0KCvL395fMsy6UOI0bh3bttixcGLx7t9WwYcvDwsKZTkRJtFOnsGABbtxAu3ZSnTt3plWQqoDwcYRDhgyp+GmvXr3atWvX9OnTq7FLU1PTixcv7tixIy8vb+bMmdOmfenJ0qpVq68DCq9fv/7hwwcHB4dnz549e/bM1tZWSUmpGvui6pM3b67y+TcB+awsqZMnr9KDQqo8vr5YuhS3b6NVK6ajUHUBM5Ou9+nTp0+fPt/duWHDhq+3f//9d/EmouqANm1avX17gccbxWJdKi7uznQcSkL5+GDFCty6RasgVVV0AQWqzjh0aMuQITdNTGwnT9Z++NDFzQ2FhUxnoiTMX39h2TJcv06rIPUT6DJcVJ2hoaFx/vxhwW0OBx4e6NwZp0+jbVtmc1GS4q+/sGIFbt4Uce9iqt6jR4RUnSQvjwMHsHQp+vaFtzfTaSgJcOTIlzOitApSP4sWQqoOc3NDaCi8vODmJkGrXVPid+QI1q5FSAjMzZmOQtVBtBBSdVurVrh/H7Ky6NwZMTFMp6HEiBCSkJCQm5t7+DDWrsXt22jenOlMVN1UzWuEbDZb6LwzFCV+jRrB2xvHjqFfP+zcibFjmQ5E1b7S0lJb22FJSQpFRe9lZedFRbmYmjKdiaqzfuKIsKio6OtMfQ4ODi9evKidSBRVHW5uCAzE6tUYNSq6des+TZt22bZtL9OhqNpy9erV2FjL1NSz+fkhMjK/0ypI1YTwQujs7Oz9Qw+EM2fOmJqaEoZWtKeoSnXogIcPERg4My7uaHJy+JYtV54/f850KKpW8Pn80lLB4tpsNpt+KVE1IqQQlpWVnT9/3sDAAEBycvK1a9cE9/fv3z89PT0xMVGc+SjqpygrE0VFAjQFZLOzu69Z8y4sDFWYdJeqYx4+HMxmR2hqjtXW7r1qlQfTcai6TUghzMjIKC0tNTQ0BBAaGurh8eVNpqOjw2KxMjIyxBqQon4Gi8Xq1MlMQWEFm+2tpXW+TRvbJUugq4sxY3DsGLKz/9O4sLCQw+EwlJSqJkKwYAGuX5dLTLweGroqLu7ytGluTIei6jYhhVCwFoTgcmBmZmZWVpbg/rS0NEKIpqamOPNR1M8KCDi6f7/5tm3FT59eXbtWIywMsbFwdMTlyzAxQefOWLMG0dGYM2dFy5b2Jia9Nm7cXfEGJ09eoKPTwcioc1DQ3Ypbnj8fMHLkNE9Pr9LSUtG9IOpfPB6mTkV4OG7ehK4u29zcnH4jUTUnpOenpqamvr7+9u3b16xZ4+3tLS0tffr06TFjxuzatUtTU7NJkybiT0lRVSctLT1+vOu39+jowM0Nbm4oLkZQEC5fxtChiampD3m8+wDf07O7isqvcnKKgsaqqmB/8/dhTEzQ6dO5hYVPgDQXl6EnTkQKlrRTVoa0NGRkIJgNXlkZISF3fv31QHb2umvXApKTV+/du7mCkIGBN27fjnRwsO3b9/tJd6nylJbC1RX5+QgKgoIC02moekT4EIgNGzZMmTLl3Llz/fv3nzt3rpOT0/Tp07Ozs7dt20bXDqTqrkaNMGgQBg3CnDnF1tZq2dkA2DyeYnR0yddCmJsLPv/fp3z4kMHhtAAAaGdlla1dC0HX6fx8cLkoK4NgibC8PPB4wcBsoFNxcfsjR2wTE6Gl9eVHR+ff21pauHTp+IIFZ3NyJh0+vP3AgdzRo38R6/+FuqmoCKNGoVEjXLyIKiw5TlE/QXghnDRpUo8ePZKTk/v06SMjI6Ompvb48WMrK6vBgweLOR9F1YZWrVp17Mh79mwCUGRnZ370qEZ5LXNy+nfosCU9nScvH+Pq6rC7/NOo5851mjTpeF5eB2npi716Wc6ahcxMZGQgPR1v3yIjAxkZyMxEejqyswMI2QaYZmebLFy47d27X5o2RdOmMDZG48b/2WZqauqkSQtfv453dR2+du0i0f0PqGPy8jBkCIyNcfgw6ABmSuRYdXo4hJ+f37Vr144fP850EKruIYRERkbKy8t37Nix4pZ5eXk3btzQ19e3srKquKWX1wE/v0udO7fZtm1FBStozpmz8sABvbKyWTIy2/v2LWvXbmlSEpKSkJiInBwIiqLg5+RJ59hYN0LsVFSmHzv2y7BhQ6vzUsUiPz9fWXDWWNQyMuDggO7dsXs3WKza2ANVb/H5fB6PJyMjU3EzWgiphqukpITNZlf6IRG5wsLCyZN/i4h42LNnD2/vbfLy8l8f4nCQmAhBXUxKgpdX18LCSIANnG3b9vWMGUttbGBh8Z+rmBKilgrhp08YMADDhmHjRpFvm6r/qlgI6VkGihI3RUXFU6f2C31IXh6tWv27ll5Ojp2Pz8LCQjt19Z2jR//x5An+/BMfPqBrV1hbw8YG1tZo1AgA8vLytm7d8+lTxvz5E9vWl4WpkpJgZ4dx47B6NdNRqHqNFkKKkly7d6/v1On406cR48bt7NLlyynctDRERSE8HGvW4MkTtGoFa2sEBo57924Il9vl6tXJjx8HNP7uYmOdsnLlFh8ffxUV7awsr8WLzTzocHmqltFCSFGSS0pKavLk70eL6+hgyBAMGQIARUV48AB375LExBQudyqA7OyowMCoSZPqak/U4ODgP/98mpMTCcQZGs708LjFdCKq/pO8Sw0URVWZggJ69cLKlSxT00ZAOJDCZt+YP9/CwQEnTqC4mOl8P+/t23d5ebaANNCWx8thOg7VIFRSCAsLC729vTMzM8WThqKo6rly5ejAgfs6dJh4/PhvKSktZs7E+fMwNISbG27dQh3qEpeebgccYrF8lZTm9+9vw3QcqkGo5NRoZmbm1KlTO3XqROcxoihJ1qxZs8BAv6+/Cs6dfv6M06exeDHS0uDqismTYWaGvLw8b+9jhJBffx2vpqbGYObvcLmYPx/BwUZBQWfCwy+YmFiNHj2S6VBUg0CvEVJUvdW4MTw84OGBmBj4+qJXLxgY8N+/d8jKcmWx2AcPDoqNDZWSkmI6JgDk58PZGTwewsKgqmrSs+cCphNRDQi9RkhR9Z+FBTw98f49Zs1Kyssz4HJnlpVNz8w0iY+PZzoaALx9i65dYWaGK1egqsp0GqrhqaQQSktLN27cWPwjjimKEjlpaTg56aqovAEygazs7LidO/VTUxlOdfcubGzg4YGdOyEZR6dUg1NJIdTX1//06VObNm3Ek4aiqFqloKBw9OhWM7NfWrQY6uu7UU1N2cICS5Z8v1Kj2Jw8iVGj4O2N6dOZCUBRoKdGKaqhcXDo/+pV6OvXYS4u9p6eePIE2dlo2RJr1nxZWEM8CMGaNVi+HMHBoJP5U8yihZCiGrQmTXDgAO7fx+fPaNECW7aAw6n1nXI4cHXFzZuIiEDr1rW+O4qqGC2EFEXB2BgHDuD2bURHo2VLHDwIHg+EkBcvXrx//15Ue1m6dJORUedOnYZ06/aOENy+DR0dUW2boqqPDp+gKOoLCwucPo2wMCxbBi8vLpc7LDtbDfg8Y8agtWsX1nDj169f37fvVW5u1Pv3jw0MZp84cVUkmSmq5ugRIUVR/2Fjg7t34eYW+vZt8/T04+npt/bvP8Hn82u42Vev3ubl9QKkgM4sVrpIolKUSFTpiPDDhw/Pnj1r3769gYFBbQeiKEoSWFvLqKoWZ2UB4BUW8ng8Vk0WQSwrw6NHA9nsMTxeIyWliEGD+ogsKEXVmPC3tpOT05o1awS3g4ODzczMBg8ebGpqeuHCBfFFoyiKOdbW1t26FenoDNDUtDIxmWttzXr5spqbSk+HvT3S000jI0+sW/fh8OGe+/Z5ijQsRdWIkBXquVyukpJSQEDAgAEDAPTo0aOoqGjPnj1//fXXrVu33r17JyFzMoGuUE/VDFMr1NchKSkpioqKysrKx45h4UL89hv+9z8IPTQsb4X6hw8xahScnbFpk/AnUlTtqeIK9ULemFlZWSUlJaampgDS0tLu37+/ePFiGxubDRs2JCcnJycn10peiqIkj56enqC8ubkhKgqXL8PeHh8/VvXpBw/CwQG7dsHTk1ZBSnIJeW8KimdJSQmAa9euEUL69esHQENDAwBdkomiGqZmzRAUhH79YGkJP79KGnM4mDIFe/YgMhK/1NVFgqmGQkghVFdXNzAw8Pb2zs3N9fb27tixo66uLoCkpCQA2traNd/rhQsXJkyYMGfOnJflXHaIi4ubM2fOhAkTAgICar47iqJEQloaixfjyhVs3IgxY5CVJbzZ+/fo2RNFRbh3D6am4o1IUT9P+NmKDRs2eHl5qamp3bt3b/ny5YI7L1++rK2tbWRkVMNdnj59evr06QMGDNDS0rKxsUn9YdLflJQUGxsbbW3tAQMGuLu7nzlzpoZ7pChKhDp3xuPHMDGBpSWCgr5/NDgY3btj5Ej8/TcUFZnIR1E/SfjwiYkTJ3bq1Onx48cdOnRo166d4E59ff3du3ezWKwa7nLbtm2bN292dXUF8OjRo8OHDy9btuzbBt7e3ra2tqtWrQJQUlKybdu2UaNG1XCnFEWJkLw8PD1hZ4cJE+DgABOT3Xv2/G1h0bR79yP79yseP46+fZmOSFFVVu44wrZt27Zt2/bbe5ydnWu+v7KyskePHvXq1Uvwa+/evYN++JMyMjLSzs5OcLtXr15Tp07lcrnS0nQSHIqSLHZ2ePIEgwYFHTx4Btjx/v2doKDZb9/+RccbU3VLudXlw4cP3t7eMTExxcXFly9fBnD58mVlZeWvNax60tLS+Hy+pqam4Fdtbe3Pnz9/1yYlJeVrAy0tLT6fn5qaKnQs//v3789onvFf6f/1Hnl5eckZ3UFRDUGBXSF6ywJDAJQgt40vHW1MSZCbrjc7G3auuI3wQvjgwYMBAwawWCxDQ8Ov3URfvHjh4+MTFxdXk0yNGjUCUFpaKviVw+EoKCj82OZrA0Hn1R/bCGhpaenH6Lc3ay/4lcViGRsby8rK1iQh1XDweDwWi8Wm/fprgM+Hv//9pKRkwBFIkpK6N2nuJDoyk5IQhBAdpcpndhdeCKdPn96xY8fz588/fvx43LhxgjsdHR2XLl2alpamU4MZ49XV1RUUFJKTkwW9T5OTk3881DMwMPg6WjE5OVlRUVFNTU3o1ho1amRFrI7PpgPqqeqgA+prKDYWbm4w10Zb+cm3bl3R0dHs3j386oIWZ8/CwoLpcBT1z4D6SpsJ+Vs4Ozv70aNH69evV1VV/bZrTNOmTQF8rPpgWmFYLNaIESN8fX0BFBcX+/v7jxw5EkBRUdGJEyeKi4sBjBw50t/fX3Db19d35MiRNe+hQ1GUCBGCgwfRrx/c3XH1Ki5dOpKW9jgpKerUqRbLlqFvX/j7V74RipIQQo4IBWcjf5wtKTs7G0DNO62sXr26b9++z549+/Tpk5mZ2bBhwwBkZma6uLh8+PDBwMBg+PDhPj4+lpaWjRs3fvv27Z07d2q4R4qiRCgpCZMmoaQEoaFo3vz7R93c0LYtRo3CrVv480/Q421K8gmpajo6Otra2levXm3Xrt23h2InT55UVFQ0MzOr4S6bN2/+6tWrhw8fKisrd+jQQXCnvr5+YmKinp4eAGlp6UuXLj19+jQ/P79Lly7y8vI13CNFUaLi74+5czFvHhYuRHn90jp2xIMHGDcO/frh1Ck0bizeiBT1k4QUQjabPX/+/DVr1vB4PAMDAz6fHxsbe+rUKU9Pz/nz58vJydV8r40aNbK1tf32HikpKcGpVwEWi/W1RlIUJQnS0uDujnfvEBiI9u0raayhgStXsHUrLC1BhxVSEk74ec7FixdnZGSsWbOGy+UCsLCwYLFYEyZMWLdunXjjURTFpKKiIhkZGRkZmbNnMXs2JkzA6dOoYr9sFguLF6NLF4wfj7lzsXhxLWelqOoSXgjZbPaOHTs8PDyCgoJSUlLU1NR69uxpbm4u5nAURTGFEOLsPCMo6CkhRW3aLE1Odvb3h43NT2+nb19ERmL0aERGYvXqt9evB7RoYTx8+C+0BxwlOSrq+WJkZDRhwgSxRaEoSnJERkbevs3JzIwAOPfvd09Lc1JUrGbpMjREcDAmTHjXpctoHm+ektLNmzej6Nq8lOQQXgiTkpLKG3thYmJSm3koipIIb98W5eUJJniSU1KSkpfnA9Wfs0leHjY21/39ZxPilp/vdvlyl337RJWUompKeCHs1q3bj4tCCPy4oj1FUfVJURE2b8b+/bZaWptKSqaz2R9dXYfWfObCFi2aqaicys2dCDwFNEWRlKJEQ3ghPHToEIfD+fprfn5+SEjIxYsXN23aJK5gFEUx4NIlzJ2LNm0QHS2rr3/9/v376urqIukfYG9v/+uvUadPd1VX1+bx/nR3x969oHPpU5KAVfUjPE9Pz8DAwODg4NrM83P8/PyuXbt2/DidYo2qDjrF2rdevoSHB9LTsXt3dTrF5Ofn/zgLR3kKCjB2LEpLcfo0VFV/el8UVUWCKdYq/Yz/xHTDrq6uISEh7969q1kwiqIkS3Y2PDzQsycGDcKDB9Wpgj9LSQkXLsDMDDY2SEqq9d1RVMV+ohDm5uYCKCoqqrUwFEWJFZ+PY8dgbo7sbMTGwsOj3MliRE5KCn/8gSlTYG2NBw/EtFOKEqpKvUbLysrevn27evVqNTW1mk+xRlEUgw4fPr5q1Q5C4ObmERQ0oVEj3LiBdu2YCePhgZYt4egILy+IYuVviqqOn+g1qq+v7+vrS1f7o6i6Kz8/f8mSnRkZ9wDW9u09vb2HTpyozmwke3tcu4Zhw/DxI377jdksVANVpV6j0tLSTZo0adu2LZ3/mqLqtJcvswsLDQF5ABoaxtbWGQDDhRCApSUiIjBkCN68wZ9/0q6klLgJf8cNGTJEzDkoiqpV6enYsQPe3oYqKsVs9jKAbWiY1fzHVZQY0qQJwsLg4gIHB/z1V/6jR0GGhoYdO3ZkOhfVIPxEZxmKouqijAysWYPWrZGdjefPWe/fX/L17XLsmGVExBWJmvBTURHnzqF160wTk94uLg/791+zatVWpkNRDcK/R4TXrl3z9Kx89r+QkJDazENRlMhkZmLHDhw8iLFj8fQp9PUFd8sMHz6c2WDlkZKCpeUVPn98YeG8wkLekSPd161bxHQoqv77txDKyMgoKSkxGIWiKFHJycHvv2PvXowejcePYWjIdKAq09BQV1B4lJ8PIJ0Q2jWPEod/C6GdnZ2dnR2DUSiKqrbQ0PBJk/5XVFQ8fvxYDY1FO3bA3h5RUTA1ZTrZTxo8ePCAAQF371qyWKySkj+OHMHkyUxnouo72j2LouqDcePmJScHADrbt48aPPhxRETHOlcCBdhs9pkzhwghLBYrNhYODsjJwYIFTMei6rWKCmFaWtq7d+8KCgq+vZMeNVKUpHnxgqSnA2gMQF6+w9SpH0xN63Z/S0EvntatER4Oe3t8/oytWyFJPXuoekV4IUxLS3Nxcbl9+/aPD9FlmChKQvD5uHIFu3fjxQuWkVGHjx89SkpMNDUDevacx3Q0kWnSBKGhGDYMbm44cgR0gnSqNggfPuHu7h4TE3P8+HFHR8fJkycHBgbOnTtXTU3N19dXzPkoivpRWhq2bIGJCdauxfjxSE5GbOyBv/+227dP5cmTm6r1a0EHdXXcvAkOBw4OyM9nOg1VHwk5IuTz+Tdu3PD29nZxcbl586aent7AgQMHDhxoYmLi6enp6uoqUWOPKKp+43K5ubm5mppfVrKNjsbBg/D3h6MjLl5E+/ZfG7Lr8TwYcnI4eRIzZ6JfP1y+DB0dpgNR9YuQI8L09PTi4uJOnToBkJOTy8vLE9w/fvz4mJiY+Ph4sQakqAYsNDTc0LCLhYVLx479/fyKrK3h5AQTE7x5g2PHvq2C9Z+UFA4cwMiRsLIC/RKiREvIEaGamhqLxRIsumRgYHD37l3B/cXFxaDLMFGUGM2YsTIl5SrQOC1t15Ytvlu3Ths4EOwGPB/U4sVQUkKfPrhyhbEVM6j6R8hHSk5OzsLC4sGDBwAcHByCgoK2bt1669atKVOmqKqqtmjRQuwhKarBSU7Gtm2Ijy8FlACwWKrjxhU7ODToKigwaxZ+/x39++PIkWczZy7dvv0Pwd/oFFVtwnuNLlmyJCMjA0Dnzp1nzZq1ZMkSQoiiouLhw4cVFBTEm5CiGpDsbFy6BF9fPHqEwYMxe/Y8H5+BfH5nZeV7kyZdZTqdpBg9GjxeoqvrFD5/nbz88zt33K9epf34qOpjVWU4REZGRkJCgpmZmaT1RvPz87t27drx48eZDkLVSSUlJWw2W4aJLvlZWVlxcXEWFhZqamoACgpw4QJOnMC9exg0CGPHYsAACJb+/PTpU0JCQqdOnSR8EbT8/HxlZWWx7c7X13fy5Cwu1wOAnl6Xz5/pIveUEHw+n8fjVfoZF35E+PHjRwMDg6+/amlpaWlpiTIdRTVg9+8/GDp0elmZrYzMnOXLjz540C4gAF26YPx4nDyJ76qJvr6+/j+zZVNfWVhYqKisyMqaBDzn89WYjkPVbcIvOHTq1Klz584HDx7Mp8N2KErUliz5MzX1cFbWrtTUPevX77Wxwbt3uHkTbm7fV0GqPJaWlps3O7ds6diz5x4dHe/Zs8HnM52JqrOEF8J169axWKxp06Y1btx4woQJQUFBfPouo6gaeP0aR45g4kSYmuLePXkgFwCLlTtoUCN3d2hoMJ2vDnJ3d3v58m5IyN/h4U1jYzF2LEpLmc5E1U3lzizz4MGDuLi4hQsXhoSE9O3bt2nTpkuWLKGDCCmqPFwu18lpur5+p86dHZKTkwG8e4eDB+HmhqZN0bs3AgPRsSNOn0ZMzGJj44W6uoOaNVu9YcNvTAev81RUEBgIAPb2+GfYM0X9hMo7y/B4vBs3bhw7duzChQslJSUSdWhIO8tQNVHFzjKEkPj4eC0tLXV19Qqa7dt3eMGCJA5nHYsVpqn5B493Sk8PtrawtUXPnjAy+n6b6enpOvVoihQxd5b5EY+HWbMQHY2rV6GtzWAQSoLUqLPMt6SkpAwNDQ0MDFRUVNLS0kQUj6IYFh0dvXjxdnl52R07lrVs2bK8ZqWlpTY2w5KSlAlJWr9+1rRpbgBSUpCY+J+fhAS8ffuBx+sCgJBOMjIfY2IqmgmMxWLVpyooCaSksH8/Vq2CjQ1u3EDTpkwHouqOigphRkbG33//7ePj8+jRIxUVldGjR0+cOFFcwSiqFuXn5zs6uqekHAGKnjxxDQ29Twg7Nxd8PnJyQMiX/2Zn4+nTwGfPOpSUbAZKPDysvLzcEhKgogJj4y8/7dtj2DAYG6OwcISDw6TMzDQ1tWsLFrjQMseIdeugowNbW1y9ijZtmE5D1RHCC+Hly5cPHz589epVLpfbp08fX1/fESNG0KH0VD1QXIyHD+HvH5+V1QloD+DzZwNb2xQ5OX0VFUhJQVUVbDbU1MBiQV0dSUn4evVAUZH4+6NZM5TzUWgbEXEyMPBm+/Yetra2YntF1Hdmz4aWFuzscPYsrK2ZTkPVBcIL4ZQpU+Tk5ObPn+/u7m5iYlIbO05JScnLy2vRogVdy4ISrYSEhFmzVqWlZS1fPn348CEA0tIQFYXoaISHIyICTZuie3czBYXo0tK7LFaRkdHnd+8al/c2LC21t7Xdl5g4mpDk9evnWVhUtOvmzZvPnt28Fl4T9XOcnaGjg+HDsX8/RoxgOg0l8YQXwitXrlhaWtZSieLz+e7u7pcuXdLQ0JCTkwsMDNTT0/u2ASFk0KBB0dHR6enpYWFh1vSPOupnDB48MS5uI9DMzW3swIGtnz83TUuDlRWsrLBiBbp0ERzPKS5c6LdqlVejRrLr1/tX8FaXlZWNjLxalc4ylETp2xdXrmDYMDx4EPLkyWFDQ92NGxdp0140lDDCC6FgDaZaEhgYePv27VevXqmpqbm5ua1bt27v3r3fNmCxWG5ubrt37+7WrVvtxaDqJUJIamoxYAOAwxmorf387FnT1q2FzFVtbm7u5+dVlV6jLBaLzjVfF3XpAh+fd/b2S/n8fVJSr588mXj//hWmQ1GSqNzOMleuXLlw4cLHjx/Lysq+vf/mzZs13OXJkyednZ0FUyxOnz7d0dHxu0IIYOzYsQDoWVPqp+TmYvVqVn6+Lpt9mM831tS8sHLlJDo9WUOWkxMtK/sLh9Oex2ufnLyV6TiUhBI+oH7BggWOjo7Xrl2rjdUHExMTTU1NBbdNTU2zs7PzqjsIlsvlxsbGRkVFRf+jpKREdEmpOoMQHDsGc3NkZSEmxnfx4pSJE6/cuHGYztLZwFlaWiorBwAvWKxzXK72f/+qp+q/27dvf/78udJmQo4IeTze/v37Z86cuXv3bikpqWrs+8WLF7t37/7x/s2bN2tqahYUFHydR79Ro0YA8vPzVVRUqrGj5OTkJ0/aWVn90ratLpvNlpaW3rZt24/ndePi4hQUFJrSgUX11LNn7N9+k+NwWH5+JV278gDpZcs8BA8VFBRU8EQGV5+ofwoLCyXwFI6uru5ffy3bsWND48Y6mZneAwfyTpzgKCpWvuQOJeEyMzPz8vKaNWtWQRt7e5d799K2bt23cKFRBc0gtBBmZGQUFxdPnjy5elUQgJqaWvfu3X+8X05ODoCurm52drbgnqysLDabXe2RxSYmJizWSj5/lo5Oj8mTxxsaQksLZWX42qeBEDJggHN0NJ/FynF27rJnz6bq7YiSTNnZWLEC585h/XpMngw2u9FPPV1GRoYWQlEhhCgpKTGdQojBgwcNHjwIAJeLKVMwfLji5ct0cte6bdu2vVu3+hKi26YN+/Zt/6+lqqQE6elIT8enT4iLw717rwiJq8pZQiGFUEtLS1dX9+PHj9XuMtOkSZPJkyeX96ilpWVERMScOXMAREZGtt9BwqMAACAASURBVGnTpsbfRFkRETp37kBaGlJS4HLB40FNDTo6UFJ6Fh3N4vP9ARw+3GHz5uUqKoo12xfFpNLS0vT0dH19fUJYR45gxQqMGoXYWNDunFSlpKXx119YsgS2trh+HU2aMB2IqhZCyLp1fxQUPAekQ0Nn2NqGENL30yekp6O0FLKyIASlpVBQACFsgCMjI1fpNoUUQikpqZ07d65cubJdu3bGxsYifxlTp05t3769t7e3iYnJ8uXLly1bJrh/6NChTk5Orq6uAPz9/bOzs0tKSgICAmJiYpydncs7d8pidTc0JElJAwlBUhJevkRsLGJj8fw5Xr1CaWkqny84OcYvKWEbGrKsrWFlBVvbr93ovygtLX3x4oWBgYGurq7IXzIlEpGRUcOHTyfEUF4+U1PzcqNG6oGB6NCB6VhU3cFiYcsW6OrCxgbXrsHcnOlAVJV9+IDgYISE4O5dFBR8qW18Pj82NpHDgYEB+vdH27awsECrVjA3h7w8Zs8eum9fK2lpP6BnxRsXPun2mDFjQkJCcnJyzM3Nvxt5U/NeowDCw8N37tyZn58/cuRId3d3wZ3r16+3tbXt3bs3gDVr1nx7hXP9+vVCT5/6+flduHDhzJkz5e0oNDSmXz+XsjINoJjNztmy5UmzZvIREbh3D0+fol07WFvD1hYWFjkODg45OW1YrJjt2z3c3Jxq/hopkevSZfDDh38CzYCjw4ennT27qIbXpOg1QhFifNLtn+Lnh//9D+fPQ9g1HIox167dnD59WWkp193dde3ahUlJCAn5Uv8KCtCzJ3r1grQ0mTOnJZerDehJScXv3btt0qQB5X2ICwoKeDyeqqpqxfsttxB+vYz3HZEUQlGpyuoTTk7Tb9+OZbE4Q4Y4paf/9vAhpk3DvHmQk8P9+wgLQ3g4goMPcTgcQuYARcbG/RISIsT2EqiqMzMb8ObNMUCPxTq1ePGbzZtX1HCDtBCKUN0qhACuXcPEiTh6FA4OTEeh/qGv3+Hz5yBARU5umKbmVkJaC4pfr15o1QpXrmD9enA4aNFiZ3DwGUCzRQteWNhFaelyxwFWcfUJkLrM19fXxcWl0mYpKSlZWVmC20+fkvHjiYYGmTuXfPjwpcGhQ3/JynoCBMiUkrKaOZOEhRE+v/aCUz+HzydeXkRN7aaqamctLXcjo06fPn2q+WY5HE5paWnNt0MRQvLy8piO8NMiIoieHjl2jOkcFCF8PgkMLJWT6w4QgCgo/O+vv24JHuJwiI8PadGCWFuTgIAv38zv379/8eIFv7KvaR6PV5XPeOXrEUqyaq9HmJiInTvh6wtHRyxbhqZNi62shiQlSUlJff799x1ZWf19fJCRgeHDMXEiOnasjexUVeXlYepUvH2LU6egrp6VkJBgYWHxdQROTdAjQhGqc0eEAnFxcHDAzJlYtIjpKA1Vdjb8/bF7N/h8ABM/fFApLdVr0uTis2dBRUUKf/6JvXvRrRuWLUOPHj+98SoeEQofUA8gMjLS2dm5bdu2Fv9MM+zl5XXkyJGfDiKRjI3h5YXXr2Figp498csvRSkpOXy+ISGqRUWfPDzw6BGuXYO6OkaOhIUF1qxBYiL4fP727XsGDBi3ffseiVqguB57+BAdO0JNDffuwdQUGhoanTp1EkkVpCgA5uaIiMDx4+jf/7iZmU23bo6xsbFMh2oooqMxbRpMTXHrFnbtQmwsZsywZbOvysn5mZq2Xrq0kZkZ3r3D3bu4dKk6VfAnCD1ODAgIkJaWbtu27dChQw0MDAR37tu3z9DQsNJDUXGq4qnRihUUkFGj9rNYewECFDdt2u3bR3k8EhpK5s4lOjrE2PhPObnZQJyS0qwdO/bWcL9Uxfh8smsX0dYmp0/X1i7oqVERqounRr+KioqTkekH5AMvTE17MB2nHnr16tXAga5WVr/cuROclUUOHCAWFsTcnHh6kszML214PJ6OTnugFCAsltukSfdTUmq63yqeGhV+RDhv3jxnZ+cnT57Mnz//6519+/Z9//79p0+farMuM0BREY6OCvLymQCAvOJimW/PFrPZsLGBlxeSkqCicq+kZA7QqqBg7tmzYQzlbRAyMzF0KE6dwsOHGD2a6TRUfZea+kZOriegBFjk5fHo+R6RGzhw3PXrc+7d2+bgsKhZs9TISHh7IzYWixf/O7kBn88vLpYSDOpTU1OaNIkjtrFsQgphWlrau3fv5s+fz2azv50zycDAAEBKSoqYoomRs/MYS8v7Ojq9tbUHaGtvtLfHj69SXh6urj0UFfcAb2Vl9zx7ZjViBMLDmYhb34WFwdISrVsjJARGlUyNRFEi0L17dyWlCyzWGSmpHcXFTT5/LveaEVUNBQUFOTmKQDeguZSU7fHjcUeOfD9wJScHbm7ScnLD1NUdtLQmmZkl9qjdk6H/IaTXqaD4kR860QiOBevlOvVycnJhYQH5+fmKiop8PnvDBlhawtsbgwb9p9nChTP5/D8CA5c5OFjPnTvD3x9TpkBBAR4ecHFB+T14qcr5+58/c+aWjU17LvfXbdukDh3C4MFMZ6IaDG1t7fDwM/v3+zVurMXj+XXvjoAA2ktONF68wLp1SkVFfDb7FJ+vqaoabGW1/Ls2t29j0iQMH44PH1Z9+PA2Jyen9hbEFU7oCdNmzZp5eHgQQkJCQr5eI/Tw8NDW1uZyuTU9ays6IrlGKNSdO8TQkMydSziciprxeCQggFhbk2bNyK5dpKCgNrLUf+fPX1JTGwZESknNMzTc9HVYS22j1whFqE5fI/zOmTNEW5tcucJ0jjouMZG4uxNdXeLpST58SJ0zZ/mECfPi4uK+bcPhkMWLib4+uXatVjJU8Rqh8EJ47NgxAK6urmvWrNHR0Tlz5szo0aMBeHl5iTpnjdReISSEZGcTJyfSpg15/rzyxuHhZMQIoq1Nli8nERHxnTrZGxp2njdvVS1lq2cmTJgP3AUIkNOmTT+x7ZcWQhGqT4WQEBIaSvT0yOHDTOeomz5/JjNmEG1tsmkTKSwst1lMDOnQgYwc+W9/GZGrUSEkhBw6dOjbydWUlJQ8PT0lqssoqeVCKODjQ7S0yK5dVWr85g2ZOZNISdkD0QBPRWXKpUuXazVe/TBunC+LNQ1Ik5X9fdq0RWLbLy2EIlTPCiEh5NUrYmpKVqygc2v8hNxcsmIF0dQkv/1GMjLKbcbnkwMHiK4uOXCgdvPUqNcogClTpnz48CEiIuLcuXPBwcGfP39evHixBK43Vtvc3BAaCh8fjByJrCwAyM/PL69x8+bYswe6uhlAR4Cdn28bFvZOfFnrIA4Hkybh+XPXefPMOnZ0mzIl7fffVzMdiqIAwMwMERG4dQtubigtZTqNpHrx4kWfPmO6dnW8fj1k506YmeHjRzx+jO3boakp/CmpqRgyBEeOICwM/0w1zbTaLce1TAxHhAKCE9lNmrw3MOimo9Pb1LTbh/IvZM2evUxZeRqLdUxJyVJHJ370aPLqlRgy1j3v35OuXcnYsYxdW6VHhCJU/44IBTgc4uxMrKxIejrTUSQPn883MOgIPAUSpKS69OqV8eSJ8Ja3bwe3aGFtZNRl2rS/dHXJ4sVEPJ+8Kh4RCu/pePfu3VJhfwKpqKgYGxtXex3duktODp6euH9/a1DQWmBgenrg0qVbjx3zEtp49+4NfftefP06YdSo0zo6pl5esLHB4MFYuRImJmIOLrlu3cL48ViyBB4eTEehqPLJyeH4cSxciF69cOUKamFhujosKyuLw9EB2gFQUurm6fmmfXshh4GEEFdXj5SUG4CKt7fDmTO9f/nFWNxZKyS8EI4ZMyY1NbW851hbW/v4+JiamtZaKgmlpFQIaAEgRDszs9wTpCwWa/jwX77+umIF5s3Dnj3o3h329li1Cs2biyOtxCIEu3djyxYcP46+fZlOQ1GVYbPx++84dAi2tjh0KOXZs5Oammrjx7vIysoyHY1JfD7OntXMySlgsU4BykpKEW3abBbasqioiMtVA3QAqKh00NR8DxiLNWtlhF8j3Ldvn6am5syZM69du/bgwYOAgAAXFxcDA4OLFy/u3bs3ISFhyJAhPB5PzFkZt2LFdB0ddzW1RXJy7unpM/PyqvpEJSUsXoy3b2FhASsruLnh3TtwudwLFy6cPHmyqKioNlNLloICODnB1xcREbQKUnXJ1Knw8sp1dBy8ZInKnDnvBwwYy3QiJsXHo3dv+PkhPPych0eMu3tISMgpJSUloY1fvFDMy1OVkVknLb1HVTW8U6dOYk5bOaEnTHv06LF169bv7pwxY4bggtyDBw8AhIaG1vj8bU2J7RrhV2lpabdv3/78OW36dNKpE0lN/ektZGSQpUuJpiYxMBihqLioUaP1LVtalZSU1EJYifPyJWndmkyfTiTk5dJrhCJUX68Rfuv69etKSksF6wTp6DSUj+13uFzi6Um0tcm+fVXqT3v6NNHRIWfPlhw9euyPP/Z9XRFPPKo/fCIzMxPAy5cvv7v/2rVrqqqqgtuNGzc+JgGreIm/EH7L05OYmJA3b6rz3Pj4HAWFvoJPlIbG9KioKFGnkxTx8fFLl27Yv//wmTMlOjqSNTCLFkIRagiF8PXr11pa/YBiIFlJqVMDfO/ExJCuXYmdHUlIqLwxn088PYmREXn0qNaDlaf6nWUIIQDi4+Nbtmz57f3x8fHkn3nXZGVl6VI4ixdDWxt9+uDSJXTo8HPPbdpUUVExo6goH5D/f3t3HldT/v8B/HWX9v22K6XEhAit9qXQQtmyq8FYxm4wWcZgGGQZ2xBZs8s2GKJIqEghWqxRRJRS3bZb3fv5/XF9m36Ilns73fo8H/NH99xzP+d1G7f3Ped8lpycJ9u26RsaomlT6QRlTnp6epcuXhkZi7jcRC534o0bgba2TGeiqJpq0aLF8uWj167trqysqqvr37cvTp6sdJBAA1NWhvXr8ddfWLECkybhuyPpBAL89BOSk3HnDups7uwa+8o9Qm1tbQcHh2nTpt24cUO8hRBy7ty5JUuWuLm5AcjMzExLSzMzM6vTpPXS+PHYuhV9+yI0tHov5HK5/v4rmzbtZWho//vvwwwMTNu3x5QpSE2VTlCG3LhxIzd3FCFDS0uXqqkl0SpIybqpU8elpNxJSgq7ds3O0REODmjAyxcKBIKioiIASUno0gVXryImBpMnf78Kpqeje3eUluLqVRmogkAl9wifPHlibm4OQE1NzdzcXElJCUDHjh3fv39PCImMjPz555/rw6SjzF4aLXftGtHTq+2yeZmZZNEioq1NJkwgyckSSsa0w4fvstnuQBEQ/8MP3ZiO8zl6aVSCGsOl0S8dPEj09MjFi0znkIKVKzfr6LTX0bHt0mWpri7ZubOqM+w8eECaNSNLl9aLGXlqO8VaQUFBYGDg3LlzR48evWDBghMnTpSWlko0oQTUk0JICImPJyYmxN+/tu1kZZElS4iODvnxx//uPvL5fFm8Lb9nD9HVJePG7TQzc7C3d/9sst36gBZCCWqchZAQcusWMTIiGzYwnUOi+Hy+jo4tIASIvHyfqKjUKr7wwgWip0eOHJFqumqo1YB6AMrKyt7e3nV3ZirjrKwQHo5+/ZCSgjVrat4Oj4c//sAvv2DzZnTuDBcXfPgwJTY2HihYtWr+Tz+NllxkKRIIMGMGIiJw4wYsLScB9WQaJYqSPEdH3LoFDw8kJsLfHw1jbKFAIBAI1MT3zjQ0eJqaBVV51ebNWL8eFy5A5m6C0PUnJcbMDDdvIjQU06aRM2fO+vltTE5OrllTmppYuhTPn0NF5c6lS4WZmZGZmdGLF68jXywSWQ+9eoVu3ZCTg+hoWFoynYaipK9pU0REICcHzs7IzGQ6Ta3x+Zg5U5vNttDSGqGj49O2rdCy8k/yiRNnjY1tDA072Nn579+PqCjZq4KoWAiPHz/epEmTv//+G4C1tXWTSjAXVQbo6+P6dfz772/Dhl1cuFC7c+dhNa6FANTVMXx4oYaGFgBAPjeXEx9f3ycxuHIFDg4YPhzHj0NNjek0FFVXVFRw8iR69ICjIxISmE5TC3fvomNHqKkhPT0gLGxBcPCMK1eCKltuobS0dNq039+8uf7uXUx8/Iljx17JaNf3/y6NmpmZeXl5iSu/h4dHXtXnTaEqUFWFSBRSVnYHYGVm4uzZi7/8MqPGrXXp0sXM7M9XryYBb9u2HeDqyu3YEX/8UR/XziYEa9di40YcPgwnJ6bTUFSdY7GwYgXatoWTExYsuH/jxnYtLbWVK+fJyvmDeO7DP//Eli0YMQIA2n9vWFhBQQEhOoAqADW1lnl57wGTOogqcf8VQnt7e3t7e/HPK1asYChPQ6Cvr5OWFgdYczg3zc09a9OUnJxcTMylmJgYTU1NS0vLggLs3o3+/dGxI5YvR8eOkopcW3w+xo1Daiqio2FqynQaimLOsGHQ08tydp4oFP7NYmVGRQ1//Pgm06G+LysL48YhPR23bqHq00jn5WkWFGgpKMzjcNT19JK+WzjrLXqPUPKOHdtqY7PI2NjBxKTJqVP9RaJatcbhcBwdHcVn6ioqmDULT5+iZ0+4u2PoUDx8KJnMNfDx48fVqzeuWvVXZGS2rS0MDBAZSasgRYHFSlBV7Qk4EjIgO5ubn5/PdKLvuHEDHTrA0hJRUdWogi9folcvrFx5/MSJHvv3t7p7N0ROTk6aMaWpsu6k//zzT9euXbW0tIyMjMRb/Pz8Nm7cKKlerRJRf4ZPfFVhIenZk0ydKpXGi4vJzp3EyIj0709iY8nmzTttbd0nTpxXN73YhUJhy5aduNwANns3h+Owa1e9G1pTFXT4hAQ12uETX8rIyNDTswXigatcbqeoKKYDVU4kIps2EUNDcuFC9V746BFp2pTs3i2dWJJTqxXqAwMDBw4cqKio6On535U9AwOD1atXN8JFJ2pMSQlnz+LOHfz+u+QbV1DApEl4+hROTujb999588JjY3ft29diwoT5kj/YF1JTU3NyjMvKJopEEzQ0WnTt+rwODkpRMkFXV/fMmS1OTqu8vIK2bDkycCA2bUL96fGdl5e3efP2TZu2PXmS6+SEc+dw9y7c3KrRQlyc+FwQEyZILWUd+7I2ikSiJk2azJo1ixASHh5efkb4/PlzAKmpVR1ZWQfq+RmhWGYmad2afLGYhyTNmrUMOA8QoMzQ0EGKR/qfBw8Kudz2QAbwQU+vo4yeDdAzQgmS0X8DdeDlS+LgQAYMIB8+MB2FEKFQaGnZRU5uK5e7jcvttHx5mVBYvRZiYoihITl5Ujr5JK3mZ4Tv379/+/btuHHjPttuYGAAICMjQ/rVuUHR0UFICHbswK5d0jrEwIE9tbS2AbFyciuKirp36oQzZ1DLe5PfcOAAnJ2Vpk7dYGk51NJy8P79a9ToUAmKqkSzZrh5E5aW6NgRN5nuN5OampqdbVhaOr2sbKqaWvNhw56xq9NR5OZNuLlhxw4MGSK1iEz4yswy4mWXxXOtVpSSkgJAQ0ND+qkaGiMjhISgZ0/x0EDJt9+zZ4/9+/MPHNhtZ9d61qxFFy/Czw8LF2LePIwdCwUFiR0oPx/TpiE2FmFhsLLqDdCldSnq++TksHYtevXC8OGYOhWLFqFa5UeCCNH/+PE5kA2wudxH1RraER6OESNw5AicnaUXkBlf+b/B4/Fat269fft2Qkj5OEpCiJ+fn7GxsYWFRd0mbCCaN0dwMGbPRnCwVNr38HA/eXKHr+9MRUX5wYNx+zZ278bZszA1xbJl+PhRAodITISjIwjBnTuwspJAgxTVqLi6Ii4OERFwcsLbtwwEuHED3bope3qusbAYYGHhvmvXH+rq6lV87YULGDECp083wCoIVNJr9Ny5c2w2u3fv3tOnT+fxeFu2bOnWrRuAgwcPSvYCbi3JxD3CiqKjiZ4euXGj7o748CEZO5bweGTmTPL4MX/kyKnNmztOm7awulOoBwYSfX1y+LCUYjKD3iOUIHqPsIo+66hZWloqkv4yDeUHvXSpJi8/fpzo6ZHbtyUdS/pqu/rE+fPnK84vZ2xsXN+qIJHBQkgIuXqV6OmRu3fr9KApKWTWLKKgMJ/N9gdKlJUXbtiwrYqvzcoinp7Ezo48fy7VjAyghVCCaCGslvBw0rQpsbdfoaPTXk+v3c6dB6R3rOxs4uFBOnUir19X41WFhYV79uzz9w/YsyevSRPy4IHU8klTbQuh2OvXr2NjY58+fVoH31lqQBYLISHk7FliaEjqflWi9u3dgXcAAaKbNPl5xQpy6RLJzv7KnsXFxRkZGYSQmzeJiQmZO5fI4DJQ30cLoQTRQlhdMTHP5eX7ASJAoKvbXkpLrd29S8zNyaxZpLr/0tu3d1JQWM3hbORyO8XFyerHpLbLMJWfCBobG0v8euzz5883bNiQlZXVr1+/8ePHfzajK5/PDwoKioyMFAgE9vb2kyZNEq8M3GB4eCAvD7163Tcw+F0oLFq9er67e786OO7o0W4vXvyalzdGQ2OTt/fU/Hz4+eHuXejrw94ednawt0eHDjh//uy0aUtZLD1lZVWBIGj3bm61xhhRFFUVbHauurrRhw8sQF4g0CgqKpKX9BpOu3Zh8WJs2wYvr+q98N27d+npigLBAgA8XjyL9QhoJ9ls9cp3CqE08Pn8bt26jR8/3sXFxdfXl8/nz549u+IOkZGRZ8+e9fDw0NDQWLduXXh4+JkzZ+o+p1SNHk2mTx8fF3ccUB83zjMhoYOenp60Dzp37s8mJvrXroUNGjS7b98+4o0iER4/RkwM7tzBoUNISkJZ2cqSkghAlcPx3b37opubh7SDUVQjZG1tbWLyuqxsNiF5hDSfNEnD3x88nmQaLyzE1Km4exc3b+KHH6r9ch6PV1CQBuQAchxOoqzMG15jDBTCw4cPN2/e/M8//wSgqKg4efLkmTNnsiv0Ju7Xr5+Li4v45zZt2rRt27awsFBZWbnuo0rPx48fFRX1cnNbAhAKHZKTk+ugELJYrGHDhgwb9v9GALHZaN0arVvDxwcABAI0bUoyMxUAyMurKikVSzsVRTVOHA4nOvpieHi4srJyx46dFy1C+/bYu7fm3TLz8vIOHDjCZrO7dBnl7a3arh1u34aKSk2a8veXV1Zew+O5sliiFSvm6ujo1DCTjGCgEMbExIj7oALo1q1bamrqu3fvKn7jqHil9PXr15qamg3s0igAHo+no1OclbVPKNQsLIxq3Xo104k+UVBA+/bNr1xxANqUlV3r0SOW6UQU1WBxuVzn/9W9v/6CuzvGjYOrKzZuRHW/+QuFQjs7l5cvRxJSBritW3d99uyvLyL4XTt3YtMm3L7tYmbmUrMWZI5UCmFpaWl6evqX2w0MDOTl5d+9e9eqVSvxFmVlZWVl5c8KYbm8vLyZM2f+8ccflS0LmZqaevXq1d69/xvWvWLFCmtra0m8CakLDj60bdu+vLynCQkn5s2T37ixvkxRHx//nJAzwHsWq01wcLBXdW8vyA6BQMBms2V4yvz6pKCgoLLPKVVFDg6IimLNmqVgZ8fevbu4bdtqzA6VmPj0/Xuz0tIZALS0Ivv0ScrPr8laMPv2ya1fL3fxYrGurqjeL5vxfSKRSKEKU4pIpRA+efLEw+MrN5ZOnTrVoUMHJSWl4uJPF9xEIpFAIFD52tl7YWFh//79e/fuPXXq1MoOZGBg0KZNm4ULF5Zvsba2VlVVrfU7qAuqqqpr1vwOgM9Hjx7Ytk3O15fhSLm5CAhAZqYcoAqYlpUdXLrU+v171TFjIIUuU8yTk5OjhVBSCCGy8tGrz1RVceoUDh7EwIHK8+fjl1++MweNSIQbN3DkCE6dal5Y+ATIA4Ryci/MzMxqcDspIAAbNiA8HGZmDeRWlEgkqspCEVIphFZWVi9evKjs2aZNm6ampop/TktLI4R8eTpYVFTk4eHRvHnzbdu2feNrpoKCgoGBgbOMT3WgpoYLF9C5M5o0wdixzGRIT8fOnfj7b3TqhI0bV69Z41xWJt+1a9uFC10OHoSNDVq0gLc3Ro4EnVWUoqRt7Fj06AEfH5w5g4MHYW7+lX0SE3HiBA4cgJISvLwQE6Nx797vvr792Gz2X3+tqkEV3LULq1fj2jWYmUngLcgYqY/j+EJkZKS+vv6HDx8IIcuXL3dzcxNvDwkJuXfvHiFEIBC4u7t7eXl9d/YTGR1H+FWJiURPj1y5UtfHvX+fjB1LtLTIzJmk4soiFQffFBeTc+eIlxdRVydeXuTcOVJaSlJSUsaP/2XcuDnJycl1HVpC6DhCCaLjCCVOPB2Mri758cdjJiZ25uaO584Fp6SQNWuIpSUxMSEzZ5J79yRzrIAAYmpKXryQTGv1h2QG1EvJpEmTjI2Nu3fvbmRk9PDhQ/FGFxeXJUuWEEKOHTsGwNTU1Px/UlJSvtpOQyqEhJDwcKKnR/73+5A8kUgUFHRqxozFISGhIhEJDSX9+xNDQ7J06deH1X8pO5vs3Em6dCFGRqVqajYs1kUW61LTpjZSGgssbbQQShAthFJy9WoGl+sIFADZcnLW+vqCmTMlPNuZuArK7Bfab6nXhZAQkpycHBkZWVRUVL7l48eP+fn5hBCBQJD9/wkrWTKrgRVCQsjBg8TMjKSnS6Xx9eu3qauPBYKVlV2MjS907EgOH672fBNily8nq6gMAwhAVFTGBgUl1suph76DFkIJooVQSh48eKCrO078WdPUdE9Lk/Bfh127GmwVJJKaWUZ6xKd6FbdoamqKf5CXl5f4DAuyYswYpKZiwACEh9dwANA3HD16MS9vD6BfWMiztw+8dq3mE8b07Gmsrv60oCAOYAMJS5aYzZkDDw94eqJnT0ku/ERRjVyrVq20tB7n5GzicIrNzQVGRgYSbHz3bqxcibCwr9+GbDwYWhSLqtzixWjfHiNGoAp9naohIQFv3rRisf4BypSVz7q7t65Na/Ly8sHB+52dLrG4ngAAIABJREFU1zo5rb55c8/jx0rie+wrV8LAAMOG4fDhT2s/Xb4c6uQ0cty4OXRJZ4qqATk5udjYy1u2qGzbZhgZeb72DZaVlZ07d+7s2bMBAaUrV+Lq1cZeBQGwCCFMZ6i5Q4cOBQcHHz58mOkgEiYUYvBg6OsjIEACrRUWYu1a+Ptj7lz+3bu+sbH3+vXrtWXLCi5XKtcDsrIQFobz53H+PMzNnz55MqGgIIDNfmxjs/vOnQvSOGKN0XGEEsTn89Vol+J6jxDSqZN7UlK70lK2SBSbmHjZwqIhj/4UD5/47mecsUuj1DdwODh0CD16wM8PtRxcePo0Zs9G796Ij4eenhqwXUIZK6WtDS8veHmhqAhLl9578GAg0EokapWQsOr4cfTsCX19aUegqAbl5cuXa9b4y8vLLVo03dDQsDZNPXr0+skTBT5/DQBt7RGEPAdaSCimDKOFsJ6qOLjQ0fGZurq6fjULyIsXmDEDycm1mr2wNpSUMGmS7f79WzMze3M4jwwNDU+dwrRp0NCAszO6dIGTE4yMPu2cmJh4+XKYnV378un3KIoCUFRU1L370Ddv/mCxioKDBz17dqu6M/gIhYiJQUgIQkNx/75maWkaIABYLFYKT1KTfMs4eo+w/jI0xPnzogkThtjbL2jbduiff26u4gtLS+HnBzs72Nnh4UNmqqCYhYVFUNBqd/ctkycnxMTsDwpCRgbOnYONDf79F23bonlzTJ6MP/+83b37T3Pnqnt6rt+xYz9jcSmq/nny5ElxcUdC3EWioa9f648b9/avv/DPP4iPR0HB/9vzwoXLlpbdW7fuGRoaBuDFC+zYgSFDoKuLyZORl4clS5CZqf7337MMDR0NDByWL5+ira3NzLuqZ+g9wnotIiLCze0Qn78DEGpo2F26FGtqyjYwwJffCDds2H7kyL/W1pYeHst+/VW9bVts2oSmTZkIXWVCIe7exbVr+Pvv+Wlp7kBPINvYeOTp05fbtoWiotQD0HuEEkTvEUrJ48d5bdv2KisLAgp5PO/ff49NSeEkJ+PFCyQnQ1MT5uZo3hxGRnnbt/fKywsBhEpK/QwNowoLlfr0Qd++cHaGgSS7msoSeo+wIWCxWHJyn76pFBeTGTOQlobsbBgZwdgYJiYwNoaxMdLTz2/eHFVQsP/+/fMnTswPCtrp6sps8CrhcGBvD3t7sNkmv/9+q7i4J4sVxeGYTJqEJ0/QogU6dvz0n7U1VFWRkpLi6flTRkZOx46tzpzZ02jH2FCNR2wsBg1Snzhxy/37vygoyG/evN/amlNxh7dvIS6K0dFvS0paA9oAOJwW27a969fPjM6CXkX0jLBeE4lE/fqNfPiwmJDMuXNH+/pOAyAQ4M0bpKXh1Su8fo03b3Dp0vLkZBugPyBs1qzLy5e3mQ5ePSUlJaNGTbt16765edNTp3bq6ekJBIiPx717uHcP9+8jIQEmJsjJGfHu3Sygk4LC6uXLtXx9p3y1tTdv3syd+2d6euZvv03p08fpG8d98OCBkpJSy5YtpfO2Ghd6Rihxp07h558REICBA7+/c1lZmaVl19TUn1gsobn5gaSkm+xvT9fdONAzwoaAzWaHhh5/+fKlhoZG+W1tBQWYm/+/oT/Xr/caNOjPjx91lZQu9OvXg5mstSAvL3/y5K6KWxQUYGsLW9tPD8vKkJQEF5d0oC0AgcB66dJby5bB0BCGhjAwQJMm0NeHkRH09TFrlndy8nxCzEePHhcdbW5WyRTC7u5jY2KKgEJPzza7dq2T8lukqOrZvBnr1+Pixf8+Bd/G5XLv3Lm4c2cgh8OeNOkCrYLVQguhDKjsT3m5Hj26HzhQcOBAoL1965kzJ9VNqrrE5aJdOyxa5L106fDcXFceb19YWGDz5nj3Dm/f4v17vHmDjAzcuIGMDJKSwifEBUBmplu7dnFqamZKSlBVhZwctLTA5UJNDSUlT65dKy4qOgngxInu8+fntGypyfS7pCgAKCnBpEmIi0NUVPVu8/N4vIUL50gtV0NGC2ED0b+/a//+snBjsBamT5/g4NAuISGhd+/TpqamAMzMvlwyhmVry4uLOy4UmuvpnQ8L89HSQnEx+HyUleHjR5SWIj8fr17JhYUVAQBIYaGgUycOh4N27dCuHays0K4d2rSBktKnFo8ePblu3R5jY4Nt2/5oWs87IFEyLisLQ4ZARwdRUdVepJ6qMVoIKVliZ2dnZ2f37X0uXTr0++8b0tOv+PpubtOmsgWFzRMTW54/3wMonThx0KpVam/fIiEBDx4gPBx//43Hj2FignbtoKcXf+DADj7/YFxc0uvXE+7fD5H4m6IosWfP0L8/XFywceN31uOtjcePHy9atEEoFK1a9UubNm2kdRiZQjvL1Hdr1mzdunW/srLi7t1+PXp0ZTpOg/L27Vs5OTldXd0vnyorw5MniI/H0aNH/v33vUg0BwCPZ5eVFVPnMWUD7SxTS6GhGDMGa9fCx0eKRyktLTUzc3jzZivAbtJkenJypGIdDFRiDu0s0xAkJSVt2HDxw4do4MPo0e5paXeZTtSgaGtrV9angMtFmzZo0waOjp1v3RqemdmBy00qLjbr2hUrVqBXrzpOSjVMT548Wb16h6qqUrNmMzdsMAgKQo9a9HUrKio6ceIkh8MZMmTwZ+UtNxeJiUhIQGTk68zMFkAXAPn5Vk+fvmjXrlbz7zcMtBDWa+Hhabm5HQAuYJCXJ19UVKqkREd/16lmzZoFB2/fvPlgs2aG8+btvXwZkyahWTOsWoXvXaOlqG/Jy8vr3XvU27d+LFYulzskISGyNgN5RCKRnZ1LcrIri1W2bp37oUNXHj1iJSYiKQmJiXjxAubmsLGBpWVTdfXnHz5cB9glJQn9+plPnow5c6ChIbk3JoNoIaynIiLg54fY2E7KyosKC5vIyb1SVGzesaPcH3/Ay4vpcI2MjY3NgQM24p+9vDBwII4exfDhaN4cfn7o2JHZdJSsSkpKKi7uDDgTAi2t3TzeB0Cnxq09f/48M9OouHgBgIcP47y8XrVvb9quHcaMgZUVzMzKp6OS8/I6tnjxBkKwYsXh/HzFtWvRogWmTsX06dCp+fFlGy2E9YtIhAsXsGIFBALMnYvTp9X4/JCgoFN6eo6enp7XrmHuXOzYgXXr6N9fxsjJwdsbI0Zg/354eKBDB6xcCWtrpmNRskZNrUVe3m3gCZCnoJBZy2k/dXV18/OfA4UA0dF5GROjrar69T1btGgRFLSj/GFQEJ4/x7p1+OEHjB6NX35Bs2a1CSKT6KBLJhFCysrKxD/n52PzZpibw88PS5YgLg7e3pCTA4/HmzJl4uDBgzkcjrMz7t+Hjw8GDMCwYXj5ktn4jZq8PCZNwtOncHaGqyuGDcOzZzh58h9Pz59WrvxLIBAwHZCq1x4+hJub9tixm7p2XdSv35aQkEPVXVPiM5s3a2lozDU07G5i0mvLlt9UKyuDX2NhgZ078egReDzY2WHAAMTGIi0tbfLkX318Zj979qw2wWQDkWUHDx4cNWoU0ylqaMeOQD29drq67adNW7F0KdHVJf37k9u3q/Ta/HyyZg3R1iYzZ5KcHCkHbbiKi4tLSkpq3w6fT9asIRoaofLy/YH7SkrLJ0/+tfbNypa8vDymI8iM4GBiYEDOnJFYg0uWEGtrkpEhgab4fLJpEzE2Fikq2rFY54BQI6OOhYWFEmiaCUKhsCqfcXpGyIySkpLffvsrIyMmM/Pe9u2RycnJt2/j/Hk4OFTp5Soq8PXFgwfIz4elJf7+m8yevbRZM/u+fUe9f/9eytmpz6mqwtcXI0bcKCmZCrQvKloUGnqD6VBUPRUQgIkTce5clWYQ/S5CMG8eLlzA1av42jigalNVxaxZuHkzXV7ekJABgHNubsfo6McSaLoeo4WQAR8/YuvW4pwcDUAeYGlpGc2enVNx7tAqMjLCnj24fBkBAce3bs1KTb195cqYcePmSSEy9X0uLrYaGgeBNyzWzpIS25wcpgNR9YxQiFmzsHUrIiIk0+WYEMyZg/BwhIRAsgsLNm2qr6b2GrgHJAqF9wYNshgwACdO4H93choaWgjrjlCI4GAMHw5zc9y9q+7gYMHjjdPUnGVm9rp9+/Y1brZdO7i4PBWJ+gFsQvo+efJcgpmpqhs40OPPP3vb2/88efKbAQPWODjgwQOmM1H1Rn4+Bg3Co0eIiICpqQQaJAQzZyIyEqGhEq6CADgczqVLB/r02dCjx/Jr1/xTUtS8vBAQAFNTzJqF+HgJH455dXCVVnrq7T3CzMzMhw8flpaWih8+fkyWLiXNmhEbG7JpE8nM/LRbZGTk1atXy8rKanm4mJgYbe3OwFkWa7KX18pattZ4SOoe4VedOkX09cmmTVJqvt6h9wi/4c0bYmNDJkwgkvrnJhKRqVNJ164kN1cyDVbRo0fE15cYGBAbG7JzJ8nLIy9evBg+/OcBA8bdv3+/TqNUTRXvEdIp1iTvyJGTs2evAyw1NZ/PnRt88qR6YiK8vDB+vBQ72d+9e/fo0X/V1dv4+w+5fp1Fl9irCmmvUP/0KYYORfv28PeHioqUDlJf0CnWKhMfjwED8OOPWLZMMg0KhfjpJyQn48IFMPIrFwpx7RoCAnD5MhEKHQoLVxOiamAw7fHjqxr1bGQ+nWKNMYsWrc/MDANUPnzYvHNn0B9//OTiAq6Uf9M2NjY2NjYAjIzg6Ylbt6BJlxViWsuWuHMHM2bA1hYnT4LOb9x45ObmHj9+QkNDXU1t8Lhx3M2bMWKEZFouK4OPDzIzcekSY8tTcDhwdoazM548ybKz4xHiBCA3t/PRo0mTJnWSxZUQZTByffXsGTZtQp8+eP2aDZQCkJMrmTaN07+/1KtgRRMmoF8/jBgBobDuDkpVRlERu3Zh4UI4OeHoUabTUHWiqKioQ4c+06cXeHs/HDRoxOnTEquCQiF+/BHZ2Th7tl4s0tSypba2djZwFYiRk4vav7+Vvj68vXH+PEpKmA5XHbQQVk90dPTatZuio6PFD8vKEBGBBQtga4uuXXHvHiZNwpEjS/T0euvpebZufXnkyOF1H3LDBnA4mD+/7o9MfZ23N8LCsHIlvL1RVMR0GkrK7t69m5PTubR0VknJSnX1zPbtC2rZYEJCwuDBE4cPn+bu/uLjR5w5899imcxisVhhYceHDTvZv7//9eu7b9/WjI2FjQ38/GBoCG9vnDiBggIAEAgEe/cGbt68LSsri+nUXyP9u5VSVMedZU6ePKul5Qwc1NR0mjnznJcX0dIiNjbE15fcvElEov/2LCgoePXqlajiprqVm0tatyYBAUwdXzZItbPMl/LyyPDhpGNHkpxMhELhu3fvGPwXInG0s0y5e/eS5eR6AwIg08DAupb/l/Pz8w0NOwKRwBUlJdvCQqGkckpVRgYJDCT9+xN1ddK/P2nRYoCi4koud7upqW1BQUGdxahiZxl6j7Cq3r7FihUnPn7cArTKybEJClq1bt2A7du/Pk2tsrKyMqNXLtTVcfo0evRAq1boShcxrB/U1HDsGAICYGeXzGYP53D0lZUzr18/RVe9b0ieP8fo0ea2tsNTUzspKirs2LGxlnOnPXv2rLS0I9AZgLq6UXZ2upGRkYTCSpGuLry94e2NzEwcP14QEpJXUrIYQEbGs23b7o4f303iQz5qg14a/SQhIcHXd8WBA0eE/7u3JhIhPh7+/hg7FmZmsLYGn9+MzY4EICcXMXy42Zgx9Xqy9h9+wIEDGD4cqalMR6EqmDQJ1tZrP3z46/37CykpK377bQPTiSiJCQ1F586YNg1RUZPevLmbnBzVp09t165s0sSCz78HxAI3FBTSDA0NJRK1zujqYto0ZW3tHOAdUMRi3Tt71qR5c7RsiTFjsHUroqP/u6H45s2bMWNmuLn53Lp1uy5D0jNCAHj8+HHv3uMzMxeqqNw8efJ+ly7rIiIQFQVNTXTpAltbzJ6Njh2Rn//r8OE/P3y4w9q61cqV/kyn/r6+fTF/Pjw9ERGB6szBS0mXomIxoAaAEDU+v5jpOJRkBARg6VIcPy7JdZsLCjBypGrv3nuBjUpK8n5+xypbSro+Y7FYR45snTzZSyAo+e23GT/9ZArgxQtERODuXRw/jrg4mJqia1f8+++o9HRfQgzv3p384ME5AwODOopYB1dppefgwYP9+/cXCr910XzXroMmJnYtWnQJCQmruJ3PJ/Hx5Nw5snkz6d59M3AAIABRULCdP5/8849kZrCtD376iQwaRBrQ3SiJqeN7hOXu3InR1+/I481QVOzQrl38hw91H0HyGvM9wtJSMmMGsbIiL15IstncXNKlC5k4kXzzL1xDUFBAbt4kfn4CRcXO4r/DHM7irl0v/for2bmTXLlCXr4kFecdiY6OadWqh4mJ7YYN/t9u+cWLF9nZ2d8NIPOFUF19lIODa/kcLhUVFJAHD9J4vK5AIZChrm49b16ZlxextSU6OkRZmbRpQ9zdyfTp5KefgpWVvYFSIMLRcUDdvxGpEghIt27kt9+YzlH/MFUICSHZ2dmRkZE5Obm//kpatyavXjGSQpIabSHMyiK9exNXVwmvA/PxI3FwINOmNa6vsK1adWex/gXu8nh2AQHpq1eTCRNIz56kaVOioEBatiSurmT6dKKlZQekAgItLdcbN+Ir63yzcuUmHZ1eV66kfve4Mn9ptLh45/37C/v0ieJyu/P5yM9Hfj7y8pCbCwUFKCik5eW1A5QApdJSXRWVvCFDtJo1g5kZ9PQqNuNiYPDg0KEupqZNDx7cxtR7kRJ5eZw6BQcHWFlhOAOjOaiv0NLS6ty5MwA/PzRpgs6dceEC2rVjOhZVTU+fwsMD/fph40ZI8Jpldjb69kX37tiwAbXraiNjQkOPLFy4Njs797fftjo6/r/rogIBXrxAcjKePxcVFwMwAZCb6zBgwMuSEqviYmhqQlUVqqpQUfn086VLgcXFMVUZ0SjzhZDNFnK5+X36yNvbf/otqKpCQwPq6uBwUFxs3bp1zOvXf8vJ5bZrJ79smVZl7axY4btihW9dJq9Lurr45x/06QN5+WdC4QNHR0djY2OmQ1GfzJoFbW306YOTJ9GtG9NpqCq7fBne3vDzw48/SrLZzEz06QMXF6xZI8lmZYKRkdGBA5u/+pSCAlq1QqtWANiXLlncvLmkuNi0SZNzDx/O1NKCSITcXPD5KChAfv6nn8PDUVxcwuF8fw5FxgrhkydPCCGWlpZffbagoODly5dsNtvCwkJeXv4b7cjJeTo7N1+40OGr35sUFRXv3Qs5dOiYikrTUaNOSyS5jGrXDuPH/zt0qJ+8vLua2urg4F02Nh2ZDkV9MmYMDA0xZAi2b8fQoUynoSoXFxe3Zs0uPT0tff1ZW7fqnjiB7t0l2X5GBpyd4ebWGKtg1Z0/H3j0aNCHD9mjRl3Q0tICwGZDSwta//9MJy/Pd/78zmz2CcDiOy1K9gpvVeTn5/fo0cPCwqJFixbdu3fn8/mf7RAWFqahodG+fftWrVoZGBhcvXq1sqYOHjzo5eUl5bwNh6OjJ5AKEODG6NEzmY7DPAbvEX5VbCwxNCT+37n9X081hnuE796909e3BW6yWEcVFZ1SUiTcfno6adOGLFsm4WYbs4KCgi9LzJcY6Im7c+dOkUj06NGjR48esdnsHTt2fLaDvb19RkbG/fv3k5KS5syZM2vWrG+0Jr2lAxoefX1tFisZAJudbGRUj4dANlY2NoiMxMaNWLCA6SjU18TFxRUW9gO6EjJCTa3I0FCS82m+fo3u3TFiBJYulWCrjZ2ioqKCgsJ3d2OgEB4/fnzcuHFcLpfD4YwfP/748eOf7aCiolJ+OdTc3FxIZ4+WkK1bl7VqtVxX11ZJ6YS29re+XlBMMTNDVBRu3MC4cQ12NXDZVVTUpqDgKvAMCNfSYn37rk21vHqFXr0wdSp++01STVLVwMA9wtTUVHNzc/HP5ubmr169+nIfgUCwdOnS3NzcmJiYbdsq7cZZUlLy7t27K1euiB9yOBxHR0elejIfbf3TtGnTxMRwAKmp6NQJtrbo3ZvpTNQXtLURGgovLwwZgqNH68UiAxSAgAAsXmy8bNmqy5cX6Opqbdx4qPZtrl27fePGXSwWRyRasnix54wZtW+SqgmpFMKQkJBDhz7/V8LhcPbt2wegsLBQUVFRvFFRUTE/P/+rjWhpaXG53Nzc3IcPH/aqZKqGN2/eJCYmrlq1qnzLsmXLOnToIIH30KDxeNi/nzNqlFJoaGGzZiKm4zBG2gvz1sbhw5g9W7FbtyIbG79nz577+AwYMsST6VDfUtkHuQEQCDBvnuKtW5xLl4patrSdPXu/eDufz69Ns2lpaevWnfjwIQYoVlPrMmZMdz5f5rvx1zcikUhBQYGZhXlNTU1dXFw+21g+M5C+vn52drb45+zs7K9OoqOgoODr6wtg2LBhtra248eP/+ri12ZmZk5OTvVthXqZ0LcvFi3C2LEqkZENf/H0ysjLy9fbQghg/360bTtz9+5WhIx58GC+ublJd8n2UJS0BrlC/Zs3GDIEzZrh3j2oSPSjUlBQUFb2A8AFVJWUdNBAf4HMEq9Q/93dpFIIf/jhhx9++KGyZ21tbaOiolxdXQFERkba2tp+oyl1dfWysjKRqPGetUjPzJmIj4e3N06ebFyDdmUFiwU+P5aQ7QDr48fxISGR9bwQNjw3b2LkSMyYgV9/lfxnpLS0XX5+grz8CkXFXCsrnpZWpaOcKWlj4Ex8xowZbm5ubdq0YbFYmzdvPn/+vHh7y5Yt9+zZ061bt7179xYWFlpYWOTk5Pz1119DhgzR0NCo+5yNwd9/o0cP+PnRbor1lL19x/T07aWlfdjs/RYWvzAdp3ERT6J94AD69JF84zduYOhQ+SNHQkWic8rKym5ubpI/BlVlDBTCzp07Hzt2bM+ePYSQI0eOdP3fcnkeHh56enoA2rRpExgYGBwcrKGhMW7cuPHjx9d9yEZCQQH//AN7e1hZoX9/ptNQX9i3b8OSJWsfPlzaosXYxYt72tujdWumMzUCxcWYOhX37yMqCmZmkm//4kX4+ODoUTg7KwF02kPmsQghTGeouUOHDgUHB9N7hLV06xY8PHD9eqP7I1ufO8t86dgxzJ2LkBC0acN0lK/h8/kyfYurrKwsOjqax+OpqrYaMgQWFti9Wypdds+exZQpOH0anTpJvnHqM+J7hMx0lqFkS6dOWL8egwcjOhr0InS9NWIEAPTti8uXYWXFdJqGRSAQ2Nm5vHljKRS+Fgrtf/vtd1/pTDx86BAWLMClS7C2lkr7VM3I3hqPlDT4+MDJCd7eoN2S6rMRI7BxI/r1Q3w801EaluvXr796ZZed7Z+be15O7uz8+VL5GGzfjiVLEBZGq2C9Qwsh9cmmTcjLw7JlTOegvmnYMGzaBBcXWgslqaxMubAwCwBQrKAglMYq8H5+2LABV6+iZUuJt03VFr00Sn0iJ4eTJz91nBk2jOk0VOW8vMBiwdkZFy/CxobpNLIvNhYzZ3YxNT2Qm9sZKFi3TvKznC1dipMncfMmmjSReNuUBNBCSP1HWxunTqFfv9LExEPFxWkTJgxvSb++1ktDh4LFgpsbLlzANwfiUt8REIDffsP27ayhQwMKCgoUFBS4XEn+VSQEc+YgIgLXr0OHTnRfX9FCSP0/7dvD1PTnFSuaENI+MHBkTMw/TZs2ZToU9RVDhkBREf374/x52NkxnUYG5eVh4kQkJ+P2bYgnP5bgxDEJCQmrV+/U1FTLzp719q1+WBjU1SXVNiV5tBBSn0tLe0DIbgDZ2W/Cwq75+HgznYj6Ond37N2LAQNw7hzs7ZlOI1Pu38ewYejdG1FRkNwaEp9kZ2f36fPju3cbWKwMFZXhGRnhdCGAeo52lqE+p6OjDtwHBGz2VUvLSqfKo+oDNzfs2wcPD0RHMx1Fdhw4AFdXrFyJnTslXwUBxMfHFxf3BHoQ4qWszBIKG+x05A0GLYTU506f3tmp07JmzXryeP1iYx2YjkN9h6sr9u+Hq2u4oaG9gUGHefP+YDpR/cXnY9QobNqEiAgMl9qMLpqalvn5N4GnQJSKSqGqqqq0jkRJCL00Sn3OwsIiKuosgJQUdO4MMzPQeRDrORcXcDjz3r27BGjv3z98xIjYb89l36hkZmaeOXPW0FDfzMx9+HC2jQ0iIqS4ymN8PDw89EePXpecvFBdXXXzZjrvlQyghZCqVLNmOHECgwbh2rV6Oq0XJUYIYbMB6AAoLGz17t07phPVF9nZ2R069Hv3boKcXAgQvHfv9pEjpXi4K1cwejQ2b8aIEd0BulSIzKCXRqlv6dIF69fDwwOZmUxHoSrHYrH69HHU1JwoL7+mpOQin9+D6UT1RXh4eHb2UKFwWnHxVhWVO1Ktgvv2YexYnDz5aTI8SobQM0LqO7y9kZiIIUNw5YpUehZQEnHw4NarV69mZmYaG4cNH67G5cLLi+lM9UBGhrFAcAQQAW/U1aX1544QLF+OQ4cQHo7KV2Kl6i9aCKnvW70agwdjyhTs3ct0FKoSLBbL2dlZ/POVK3BxAZ+PxryCWWkp/vwTAQH2Q4bYRkTYqqmpHDiwVRoHKinBhAl49gxRUdDTk8YRKKmjhZD6PjYbR46gRw9s3Ig5c5hOQ31P69YIDUXfvsjLw+zZTKdhQlISxoyBmRkePoSOzgJAWgtP5+Rg8GBoauLaNdDBgrKL3iOkqkRZGWfPYuNGnD/PdBSqCn74ATduYPt2/NHIBlMQgoAA9O6NKVNw6pR0ZzV78wY9e6JtW5w8SaugbKNnhFRVNWmCM2fg7o7QULRty3Qa6ntMTXHjBvr2RWEh1qxhOk2deP0aPj4oKkJEBCwspHIIkUh09+5dTU3NoqIWHh6YOxczZkjlQFRdomeEVDXY2GDrVnh6IiOD6ShUFRgY4Pp1XLuPMKC6AAAOYElEQVSGadMa/kqTJ07Azg7du0uxCpaVlTk4uLm67rCzm9mp0/K1a2kVbCBoIaSqx8sLY8Zg8GAIBExHoapASwshIXjwAD4+KCtjOo1EFRQUXLhw4f79+7m58PbG0qW4eBHLloHDkdYRo6KikpNbZGXtyc29qKBwfvDghvULbcRoIaSqbflyNG2KSZOQkJAQT9eHrfc0NBASgowMjBmD0lKm00hITk6OlVXPkSNv9ey53MTkTyUlxMaiY0fpHlReXrGoKA8AUCYnVyaN9XspRtB7hFS1sVjYuxdGRj+dOpWnpMTu1Enh3LlApkNR36KsjPPnMXIkBg3CiRMNoWfH5cuX09OHCgS+gIjHs9+5c7G0j8jn488/7TQ0FFRVe7BY/OXL59BC2GDI9v/I3NzcrKwsplM0Rjk56Wx2WkFB0IcPx6Kjc1JSUphOVBNJSUkymrwG5OURFARtbXTqdL1NGycrq96hoWESbD88PLyoqEiCDX5bXJxmaekrAMBHNTWpXQz9nxcv4OiIJk1YaWkBT5/+8+pV1JQpPtI+KFV7WVlZt2/f/u5usl0IExMTnzx5wnSKxkheXp7FygcIgNJSvoKCAtOJauLgwYNnzpxhOkXd4XCwbVvR06ezk5IOJSYeHTNmfkFBgaQa//3335OSkiTV2jc8fw4XFwQH9+3Zs0xX19bIqM+uXdLtFBsZia5d8fPP2LkTXC60tLQUFRWlekRKUiIjI9evX//d3WT70ighhOkIjZS2tvZPP7nv2eNYUsIqK3MrKzNkOlENNbZ/QhkZ6erqLYqKDAGUlbV+/fq1paUl06GqqqgIfn7w98eiRZg+ncXh7BQKhRzp9Y0BAOzZg8WLcfgwnJykehxKKqr4AZftQkgxaPXqhYsXzwCwY4eqqysiI6GhwXQm6nuaNWumpvb6w4d9hLBzc58GB1vISh08fx6zZqFNG9y9C2PjTxulWgWFQixejFOnEB4OWfktUTVDCyFVc+IVR+fNw8uXGDkS589Lsec6JRFsNvv27X///nsPIXB1vTBlCvfmTezdC01NppNVLi0Nc+bgwQPs2IG+fevooPn5GD0ahYW4cwdaWnV0UIopLJm+NDRixIhz58516dKF6SCNHSGcxMRftbXvGBpeZTpLNTx9+lRBQcHU1JTpIIwRieSfPZtUUqJhZeXHYtVqVFx0dHTr1q3V1NRqn4oQkp+fr6CgIC8vn5fX4unTqTo6t01NT9UyYbUkJc3jcgtatAhgsYR1dlBK4jIzMwHExcV9ezfZLoTPnj27cOGClZUV00EomZSVlSUnJ6eurs50kIYgNTXV2NhY2nfsKKpaBAKBsrJyr169vr2bbBdCiqIoiqol2R4+QVEURVG1RAshRVEU1ajRQkhRFEU1arQQUhRFUY2arI4jfPHixevXr8sfqqio2NraMpiHkiFXrlzR0NCws7Mr3xIRESESibp3785gKhl17949HR0dExMT8cPo6GgFBYX27dszm4qiUlNTK84krKio6ODgUNnOstprdO7cuUePHm3VqpX4oZmZ2e7du5mNRMkKNputrKz85s0bDQ0NACkpKRYWFrq6uunp6UxHkz19+/Z1dXWdM2cOAD8/P39//5CQkJYtWzKdi2rsli1b5u/vXz64Tl9f/8iRI5XtLKtnhAD69eu3b98+plNQMsnBwSEoKGjixIkAAgMDe/bsmZiYyHQo2bZs2bJjx47dvHmzadOmTGehKADo3r37iRMnqrInvUdINUY+Pj6BgYEACCGHDh0aO3Ys04lkGCFk9uzZJ0+eDAsLo1WQkkW0EFKNUZcuXbKysh4/fnz9+nUjI6PmzZsznUiG+fn5Xbt2LSIiokmTJkxnoaiaoIWQaqTGjh178ODBwMDAH3/8keksss3a2vrly5dVWf6UouonWgipRsrb2/vgwYOXLl0aOnQo01lkm6ur69GjR4cNG/bvv/8ynYWiaoIWQqqRMjY2HjRo0JQpU8SLSVG14e7ufvTo0VGjRl24cIHpLBRVbTLca5Siamnz5s1MR2g43N3d9+/fP2LEiGPHjrm7uzMdh6KqgbNs2TKmM9SEpqZm27ZtG/NKclSNGRoaOjg4KCgolG+Rl5e3sLBo164dg6lkFI/H69Chg76+PoBWrVp17dr13bt3bdu2ZbFYTEejGjV1dXUrKytzc/Oq7CyrA+opiqIoSiLoPUKKoiiqUaOFkKIoimrUaCGkKIqiGjVaCCmKoqhGjRZCiqIoqlGjhZCiKIpq1GghpCiKoho1WggpivqO06dPh4aGMp2CoqSFDqinKOo7bG1tzczMqrjGKUXJHHpGSFF158OHD8XFxTV7bU5Ozvv37yUSIzMzMzs7u7Jns7Ozs7KyqtjUx48fMzIyvr3Dx48fq5ePouoWLYQUVUPBwcE8Hu/p06fih7t37+bxeOWT9/L5fF1d3V27dgHIysoaNGiQmpqarq6uioqKpaXlqVOnxLu9ePGCx+Pt27evYsunT5/m8Xj3798XPzx79mzr1q21tLQMDAxMTEyOHTv21TwhISE8Hu/GjRsVN65du1ZfX7+87O3YscPExERPT09bW9vKyio8PLzizrt27TI3N9fW1tbR0dHW1l6/fj2A9u3bx8XFnTt3jsfj8Xi8MWPGiHe+fPlymzZteDyevr5+s2bNjhw5Ut7Opk2beDxeTEyMtbU1j8cbMmRI9X6zFFXHCEVRNZKdnc3hcLZv3y5+OHToUHl5eTs7O/FD8eJ88fHxhJCUlJSJEydevHgxKSkpIiJi8ODBXC733r174j0dHR27dOlSsWU3N7eWLVuKfz5z5gybzfbx8bl169a9e/emTp3KYrEuXbr0ZZ7S0lIDA4Nx48ZV3PjDDz94enqKf167di2bzV64cGFsbOzt27c9PT2VlJSSkpLEz65btw7A6NGjw8PD4+PjDx8+vG7dOkLIzZs3W7Ro0a1bt9DQ0NDQ0Li4OEJIdHS0nJxc165dw8PDo6KiBg4cyGKxTp06JW5q1apVAMzMzDZs2BAVFRUWFlbLXzVFSRUthBRVc7a2tkOGDCGECIVCHR2dqVOncjic7OxsQsicOXP09fVFItGXryotLTU2Nvb19RU/9Pf3B/D48WPxw3fv3nG53DVr1hBCRCJR8+bN+/btW/Hl3bt3d3Z2/mqeOXPmqKqq8vl88cPIyEgAp0+fJoTk5OSoqKhMnz69fOfi4mJTU9Off/6ZEJKbm6uqqtq/f/+vNmtjYzN06NCKWzw9PTU1NXNzc8UPy8rKWrRo0a5dO/FDcSH09/ev7PdGUfUKvTRKUTXn5OQUFhYmFArj4uKysrIWLlyooqIivt549epVJyen8tWI+Hz+rl27fH19J0+ePG3aNJFIlJycLH5q5MiRSkpKhw8fFj88dOiQSCQaPXo0gBcvXiQnJ7ds2fJKBcbGxvHx8V/NM378+Pz8/H/++Uf8MDAwUFtb283NDUBkZGRBQYGxsXF5Ozdv3jQ1NU1ISABw+/bt/Pz8CRMmVPGNx8XFubq6qqurix9yOJyhQ4fGx8dXvB3o6elZ9d8kRTGILsxLUTXn5OTk5+cXFxcXFhZmZWVlbGzcrVu3q1evduvWLSEhYdasWeLdHjx44OTkpKSk5OzsrK2tzeVy5eTkcnNzxc9qaGh4enoeOHBg2bJlbDb7wIED/fr1MzY2BiDuHbN///7yMllOJBKx2Z9/kbWysurQoUNgYOCYMWOKi4uDgoJ8fHzEKy+Km1q9evVnrzIzMwPw4cMHAOKDfpdIJEpLSzM0NKy4sUmTJoSQ7OxsLS0t8RbxIoUUVf/RQkhRNdetWzdFRcWrV6+GhYU5OTkBcHJyCggI6Natm0gk6t27t3i3TZs2qaurx8fHq6ioiLdcuHChYjs+Pj7Hjh27fv26hobGw4cPFy1aJN6uoaEBYMOGDZMmTapiJB8fn19++eXVq1dRUVE5OTk+Pj4Vm/rnn3969uz55as0NTXxv2L5XWw2W0lJSVw7y2VmZpYfpXy3KmamKGbRf6kUVXOKioqOjo7BwcERERHlhfDx48eBgYHNmzdv1qyZeLeXL19aWlqWV8Hnz58/efKkYjt9+/Zt2rRpYGBgYGCghoaGh4eHeLulpaWenl61BvCNGTNGTk7u0KFDgYGB4hNE8fbOnTvLyckFBQV99VUODg4KCgqVPauqqlpUVFRxi6Oj45UrVwQCQfmW8+fPW1hY6OjoVD0qRdUXTN+kpCjZtmLFCgBcLlfcc0QkEunp6QGYOHFi+T6zZs1SUlK6dOlScXFxbGxshw4dVFVV+/TpU7GdBQsWqKqq6ujoTJkypeL2gIAAAOPHj3/8+HFhYWFycnJgYKC4K01lBg0aZGxszOFw/vrrr4rb58+fz2azlyxZ8vLly8LCwsePH2/ZsmXfvn3iZ319fdls9uLFi1++fJmfnx8VFRUYGCh+aurUqVpaWqdOnYqNjX3+/Dkh5PLlyywWa/jw4SkpKW/fvp05cyaAgIAA8f7izjI1+nVSFAPoP1aKqpWoqCgAFcc/jBgxAsDx48fLt3z48KFLly7ir54KCgorVqxwdnb+rBCWnyPeunXrs0Ps3r3bwMCg/Murrq7uZxXuM+LOMlwuNz09veJ2oVD4xx9/lPdwAWBqahoUFCR+tqysbPHixUpKSuKnuFzuggULxE+lpaW5ubmJL58OHDhQvHHv3r3ltwOVlZVXr15dfiBaCCnZQqdYo6i6QAh5+fJldnZ2y5YtK5aiKhKJRE+ePOHz+fr6+uKzvRonKS0tffTokUAgaNKkiZGR0WfPis8U2Wy2qalpeZ2rTElJSVJSUllZWevWrZWVlWsciaKYRQshRVEU1ajRzjIURVFUo0YLIUVRFNWo0UJIURRFNWq0EFIURVGNGi2EFEVRVKNGCyFFURTVqP0fRTpUTmMZIoIAAAAASUVORK5CYII=", "text/html": [ "\n", "\n", "\n", " \n", " \n", " \n", "\n", "\n", "\n", " \n", " \n", " \n", "\n", "\n", "\n", " \n", " \n", " \n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n" ], "image/svg+xml": [ "\n", "\n", "\n", " \n", " \n", " \n", "\n", "\n", "\n", " \n", " \n", " \n", "\n", "\n", "\n", " \n", " \n", " \n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n" ] }, "metadata": {}, "execution_count": 1 } ], "cell_type": "code", "source": [ "using DFTK\n", "using Unitful\n", "using UnitfulAtomic\n", "using LinearAlgebra\n", "\n", "# Define the convergence parameters (these should be increased in production)\n", "L = 20 # height of the simulation box\n", "kgrid = [6, 6, 1]\n", "Ecut = 15\n", "temperature = 1e-3\n", "\n", "# Define the geometry and pseudopotential\n", "a = 4.66 # lattice constant\n", "a1 = a*[1/2,-sqrt(3)/2, 0]\n", "a2 = a*[1/2, sqrt(3)/2, 0]\n", "a3 = L*[0 , 0 , 1]\n", "lattice = [a1 a2 a3]\n", "C1 = [1/3,-1/3,0.0] # in reduced coordinates\n", "C2 = -C1\n", "positions = [C1, C2]\n", "C = ElementPsp(:C, psp=load_psp(\"hgh/pbe/c-q4\"))\n", "atoms = [C, C]\n", "\n", "# Run SCF\n", "model = model_PBE(lattice, atoms, positions; temperature)\n", "basis = PlaneWaveBasis(model; Ecut, kgrid)\n", "scfres = self_consistent_field(basis)\n", "\n", "# Construct 2D path through Brillouin zone\n", "sgnum = 13 # Graphene space group number\n", "kpath = irrfbz_path(model; dim=2, sgnum)\n", "plot_bandstructure(scfres, kpath; kline_density=20)" ], "metadata": {}, "execution_count": 1 } ], "nbformat_minor": 3, "metadata": { "language_info": { "file_extension": ".jl", "mimetype": "application/julia", "name": "julia", "version": "1.8.3" }, "kernelspec": { "name": "julia-1.8", "display_name": "Julia 1.8.3", "language": "julia" } }, "nbformat": 4 }