{ "cells": [ { "cell_type": "markdown", "source": [ "# Cohen-Bergstresser model\n", "\n", "This example considers the Cohen-Bergstresser model[^CB1966],\n", "reproducing the results of the original paper. This model is particularly\n", "simple since its linear nature allows one to get away without any\n", "self-consistent field calculation.\n", "\n", "[^CB1966]: M. L. Cohen and T. K. Bergstresser Phys. Rev. **141**, 789 (1966) DOI [10.1103/PhysRev.141.789](https://doi.org/10.1103/PhysRev.141.789)" ], "metadata": {} }, { "cell_type": "markdown", "source": [ "We build the lattice using the tabulated lattice constant from the original paper,\n", "stored in DFTK:" ], "metadata": {} }, { "outputs": [], "cell_type": "code", "source": [ "using DFTK\n", "\n", "Si = ElementCohenBergstresser(:Si)\n", "lattice = Si.lattice_constant / 2 .* [[0 1 1.]; [1 0 1.]; [1 1 0.]]\n", "atoms = [Si => [ones(3)/8, -ones(3)/8]];" ], "metadata": {}, "execution_count": 1 }, { "cell_type": "markdown", "source": [ "Next we build the rather simple model and discretise it with moderate `Ecut`:" ], "metadata": {} }, { "outputs": [], "cell_type": "code", "source": [ "model = Model(lattice; atoms=atoms, terms=[Kinetic(), AtomicLocal()])\n", "basis = PlaneWaveBasis(model, Ecut=10.0, kgrid=(1, 1, 1));" ], "metadata": {}, "execution_count": 2 }, { "cell_type": "markdown", "source": [ "We diagonalise at the Gamma point to find a Fermi level ..." ], "metadata": {} }, { "outputs": [ { "output_type": "execute_result", "data": { "text/plain": "0.3842999767139065" }, "metadata": {}, "execution_count": 3 } ], "cell_type": "code", "source": [ "ham = Hamiltonian(basis)\n", "eigres = diagonalize_all_kblocks(DFTK.lobpcg_hyper, ham, 6)\n", "εF = DFTK.fermi_level(basis, eigres.λ)" ], "metadata": {}, "execution_count": 3 }, { "cell_type": "markdown", "source": [ "... and compute and plot 8 bands:" ], "metadata": {} }, { "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Computing bands along kpath:\n", " Γ -> X -> U and K -> Γ -> L -> W -> X\n", "\rDiagonalising Hamiltonian kblocks: 6%|█ | ETA: 0:00:07\u001b[K\rDiagonalising Hamiltonian kblocks: 19%|███ | ETA: 0:00:02\u001b[K\rDiagonalising Hamiltonian kblocks: 31%|█████ | ETA: 0:00:01\u001b[K\rDiagonalising Hamiltonian kblocks: 44%|███████ | ETA: 0:00:01\u001b[K\rDiagonalising Hamiltonian kblocks: 59%|█████████▌ | ETA: 0:00:01\u001b[K\rDiagonalising Hamiltonian kblocks: 75%|████████████ | ETA: 0:00:00\u001b[K\rDiagonalising Hamiltonian kblocks: 88%|██████████████ | ETA: 0:00:00\u001b[K\rDiagonalising Hamiltonian kblocks: 100%|████████████████| Time: 0:00:01\u001b[K\n" ] }, { "output_type": "execute_result", "data": { "text/plain": "Plot{Plots.GRBackend() n=48}", "image/png": "iVBORw0KGgoAAAANSUhEUgAAAlgAAAGQCAIAAAD9V4nPAAAABmJLR0QA/wD/AP+gvaeTAAAgAElEQVR4nOyddVgUaxvG710QQcEkVQwUVGwUEwuxsTCPgY2NLWKBjfmBjY0diIqt2GKC2IEeMbEFpVl2n+8POGIQyzI77y47v+tc59qZnX2fGxz2nreeR0REEBAQEBAQ0FTErAUICAgICAiwRDBCAQEBAQGNRjBCAQEBAQGNRjBCAQEBAQGNRjBCAQEBAQGNRjBCAQEBAQGNRjBCAQEBAQGNRjBCAQEBAQGNRjBCAQEBAQGNRjBCAQEBAQGNRr2NMD4+/tGjR6xVAEB8PM6eZS2Ca4iIYQY+qVQq55WRkXj/PsftE+Hy5QzO37qV46aUh0SCsDD+wslkMiW1/Po1Pn9WUtvMkP8WVZjERDx8qOwgmYbeuhU/7wiZLO1PQyaTxcbi0CE2qhRAni8x9TbCx48f9+vXj7UKAJg1C927sxbBNSkpKUlJSayix8fHy3nl0qVYvz7H7f/4gb59cfLkbyclEvTogatXc9yaknj9Gq1aITqaj1gymSwxMVFJjXt6YudOJbXNDPlvUYU5dQqursoOkjFfv2LTJjRsiLt3AeDzZ7Rti5cvkZCQsHgxDRzIRlVOkclkKSkp2V6m3kaoOmzfjh49WIvQSIhw4IAiTyGFC2PHDgwahA8f0k/mywc7u70tW9afOHEuhyIVpnx5dOgAb2/WOnKHRILAQDg5sdahhhw8iC5d2IQuWRKXLmHkSLRpg7FjUaAARo/GrFkAsHOnqFUrNqqUhGCEHPDgAT5/xrx5rHVoJMHBKFIE1taKfLZxYwwahAED8OvYSUjIuvj4wM2b93OlMJdMn45VqxAVxVpHLjh3DlZWKF2atQ51QyLB0aPo1ImZAJEIzs54+BCJiahcGRUqICgI166JIyIwdSozVcpAMEIOcHODlRUMDVnr0Ej27kXPnop/3NMTP35gxYr0M3PnTjA17VagwEDlTwDJRWqn8FeFaoe/P7p1Yy1CDblwAZaWMDdnLKNYMfj6YtcuLFiA/PnPtmzZQSRaVKsWY1XcIhghB5w9i/HjWYvQSGQyBATkanZWWxs7dmD+fNy5k3amW7cO799fqlRpvALzjkpixgysXKmuncKUFAQGomtX1jrUEIbjon/TpAnCwvD27RCgq0x2fF7eGgETjDC3bNwIIgwdylqHRnLxIszMYGmZq0YsLLBqFfr0wa9LH1auhKcnPn3KpUBuUOtO4fnzKFsWZcqw1qFuyGQ4dAidO7PW8Qv586NOnfLAMiC8s0opyzWCEeaW5cvRvDnEwi+SBbkcF/1Jjx6wscHEielnrK3Rvz/c3DhonBPUt1MojIsqxrVrMDKClRVrHb9z7VrQnTsHkpJeVa1albUWLhG+v3NFdDSePMHs2ax1aCQpKTh0iLNdK2vX4uxZ7NuXfsbTExcu4MIFbtrPJWraKZRKceiQsF5UEVRqXPRXKlSooK2tzVoFxwhGmCtmzEDx4qhXj7UOjeTcOZQrh7JluWlNXx87d2LMGLx+nXamQAEsXYrRoyGRcBMil6hjp/DiRZQujfLlWetQQw4eFB4g+EMwwlyxezf69mUtQlPhalz0J7a2cHVFv374uV60a1eULYuVK7mMojDq2CkUxkUVIywMRKhenbUOjUEwQsW5dQvR0fDwYK1DI5FIcPgw94/M7u749u2oqWnTOXPSNrH7+MDLC+/ecRxIMdSrU5i63ENYL6oABw/mwUxVqoxghIozdSqqVEGRIqx1aCRnzsDamvs92mIxkpOXf/my73//25x6pnx5jBjx2zoahpQvD0dHVemhZsulSzA1RYUKrHWoIQEBKjpBmFcRjFBBZDJcvowpU1jr0FT27lVWTjt392GGhp20tbv9nBqcNg1hYThxQinhcsrMmVixgqfso7lEGBdVjGfPEB0trDzgFcEIFWTVKmhpCROEbEhKwpEjylpKMGBAz8+frzdqNGvRorQz+fNjxQqMHQt2GcjTSe0Uqv5MoUwmLPdQEH9/dOkCkYi1Dk1CMEIFWbkSeSztrBpx8iRq1kSJEkoMsWYNVq5Mr4DTujWqVsXSpUqMKD9q0SkMDoaRESpVYq1DDVHZjRN5GMEIFeHTJ/z7L+aqRH0CTYTz9aJ/U6IE5szB4MHpK0h9fODjg4gI5caVB7XoFArjoorx9i0iItCkCWsdGoZghIrg7g4TE2FxMxsSE3HyJB+PzC4uKFQIPj5ph+bmmDABo0YpPa48qHinkAgHDwrrRRXhwAF07Ig8t2Fd1RGMMMccP358585d/fopq4SpQNYcPQpbWxgbKz2QSIT167FoEZ4/TzszYQJevsSRI0oPnS0q3im8ehWFC6NyZdY61BBhXJQJghHmjPDwcEfHMUlJVy9cEB532cDDuOhPypaFuzuGDEkrWKijg3XrMGYM4uJ4EpAFqtwpFMZFFePLF9y9CwcH1jo0D8EIc4ZEok2kB3zX1RUGLxgQF4egIF5T8ru6QiqFr2/aYZMmaNQICxfyJyAzypdH+/aquKeQKLe1sTSWgwfRpg10dVnr0DwEI8wZHh4W+fIFzJ5d+9y5A6y1aCKBgWjUCMWK8RdRLMbGjfDwSM9Bunw5NmzAkyf8acgM1ewU3rgBPT1YW7PWoYYI46KsEIwwByxejKtX0a+f1axZ4/Je/nW1gM9x0Z9UrIixYzFoUNoAqYkJpk3DmDF8y/ibChXQrp3KdQr9/Rn8G+UBvn9HcDDatmWtQyMRjFBezpyBtzckEkyaxFqKpvLjBy5eRMeODEJPmYLoaGzblnY4ejS+fv2tZhMrVLBTePCgMEGoCEePolkzGBiw1qGRCEYoFy9fwtkZPXqgcWNhLRwzDh1Cs2YoXJhBaG1tbN6MyZPTsm9racHXFxMm4Pt3BmJ+RdU6hTdvQlsb1aqx1qGGCOOiDBGMMHtiY9GhA2bNwrFjQneQJfv2sRxzq14dw4ZhxIi0Q1tbtGmDOXOY6fmJSnUK/f2VlQM2b5OQgKAgODqy1qGpCEaYDUQYNAh16qBYMZiZoWFD1oI0lehoBAcz/qaYORMvXmD//rTDxYuxaxfu3mUpCSrWKQwIEMZFFeHkSdjawtCQtQ5NhcGKj4iIiH///ffnYePGjfPnz//HNUR04cKFN2/eNGrUqDzT+tZz5+LtW5w/j6ZNMX06QyGaTkAAWraEvj5LDTo62LQJnTqhWTMYGaFYMcyejdGjcekS4xTJM2eiQQOMGcO4KFhoKIhQowZLDWqKMC7KFgY9wl27dg0cOHDRf8TExPx9Ta9evcaOHXvp0qV69eoFBgbyLzKVwECsXw9/f1y5gu/f0b49KyECSqy7lCPq1UOfPpgwIe1wyBBIpemLaFhRoQLatsWqVYxlCOtFFUMiwfHjvO6OFfgT4p158+aNGDEiiwtu3LhhaGj4/ft3Itq1a1eVKlUyuzIkJMTGxoZ7iURE9OQJmZjQjRtERK1b0+bNSoqjuiQnJyckJLCK/uPHj5+vP3+mwoUpNpaVlt+IiyNLSzp0KO0wNJSMjenLF6aaiJ49I0NDiopSvAWpVBoXF5cbDRUqUGhobhpQM369RXPDyZPUsCEnLfFBbGysVCplrUJepFJpcnJytpexmSP88uXLqVOnHjx4QKk7s37nyJEjrVq1KlSoEIDOnTs/ffr016FUfoiJgZMT5s9H3bq4fx8PHqBPH54lCKRz4ADatUPBgqx1AAAKFMCGDRg1ClFRAGBjgx49MHMmY1UVKqBmzQv16g04eTKIiYCwMEilqFWLSXD1RhgXZQ6DOUKxWPzixYuVK1fevn27YsWKR44c0f995ufdu3clS5ZMfa2np1esWLF3795lNlP4/fv3tWvX/jwcOHCglpZWLhXKZPjnH+0WLcjZWSqRwMtLa/RoEolkP0uWawgSiUQqleb+96lwdMl/v/E9e7THjJFJJDImSv6mYUM4OmpNnoy1a6UAPD1Rvbp2nz7SunUzeLDjjUePpkRGbhk0qN+rVzcU+LhMJvv1d55T9u3T6tYNKSnS7C/NK+Tm1/UTmQyBgfnOn0+RSFjePPIjkUi0tbXFYvVYaCmTyURyTOAzMMIpU6a4u7sDSEhIaNq06ZIlS2bPnv3rBVKp9Nffsra2dkpKSmatJSUlhYaG/jzs169f7hXOmqUdFUXz5iVLpXj3TnTypPby5clSDfoDT0P6HwyjA/j4URQWBnt7iUr9E8yfL7W1zX/qFDk4yAoWxNy5GDNG6/LlZEaPDQDg6Nhi+/ZBRYrYK/ZPJpPJFP7nHjrUbc+eaxMmjJRKeynwcTWFk7+O4GCxkZGsTJkUlbq9syD1p85wME8Fkclk8jzKMzDCn7L09PS6dOly7dq1Py4wMzP79OlT6uuUlJSvX7+WyLwYubGx8caNGzmUd/Ag9u/HzZsoVEgXwKpVGDIEJiaamAdXS0tLKpXqMsoBLJFIUkMfOYJOnVCkiGr9E+jqYsMGDBmic/8+9PUxaBB27MD27brDhzOT5Ou7cN68lOrVtZ8+VWTppkwmA6DAPzcRBQaeT0k5dPr0yEWLBuQ4sNry8xbNDceOoWtXRX7trEj9TlCjHqE8DyuMf5iwsDBzc/PU19HR0ak9v6ZNm54/fz719cWLFw0NDStUqMCPnnv34OICf38YGQFAVBT8/DB6ND/BBTKGSX5ReWjRAk2bYsaMtMNVq+Dpif8e4dhgZKTt4YHRo8Hn87pUKipQYES5ciMWL3bjL2pe4dAhODmxFqHxMDDCLl26TJs2bdmyZU5OTufOnZs4cWLq+VKlSp0/fx5A69atjYyMunXr5uPjM2TIkKlTp/KT4frbN3TtihUrULt22pnVq9GlC/5z6jTCwsJ27dqTmCgU5uWDyEg8fqy6Fdp8fBAQgCtXAMDaGs7OcGPtBS4uSExM3/XPA8uWoWbNES9enGrZsjl/UfMEoaHQ1kbVqqx1aDwMjNDV1bVAgQKfP392cHAIDw+3sLBIPb9ixQpra2sAYrH43Llz9vb279+/X7t27WheemRSKfr0Qdeu+OeftDNJSVizJn3HWCqfPn1ycBg6YEDo9OlePKgS2LsXnTtDR4e1jkwoXBirV2PIECQkAMCsWbJ9+4aWKtX40qUrrCSJxfDxwaRJPFUPfvMGy5Zh9Wo+YuU9Dh5EV6HCtwrAYI6wefPmzZtn8OQ4aNCgn6/19fVdXV15FIXJkyGVYv789DNbtsDWFlWq/HaZTKbz44dEKn2bP78tn/I0lr17MXcuaxFZ0qEDtm/HnDlYuBDfv7/X1n7z7t2CZct2Nmlix0pSw4Zo0gSLFvGRCnXcOLi64r+nWYGcERCArVtZixBgYoQqyM6dOHQIt27h5/IimQzLl2Pz5j+vnD69SIcOZ6TSlwUK1OVZpAby+jVevEBGT02qxZo1qF4dXbuidu0SXbtWCQyc/eaNV3Iyy47s0qWoXh39+sHSUolRTp/GvXvYuVOJIfIw4eH48QO2whO1CqAeK3+Uyp07mDABgYEoXjz95MGDKFoUdr8/02/ahOBg+PkZL15cd/XqtNEwAeWxZw+6doXql0A2NMSSJRg8GBKJaPPmZZ8+BZUpU4fh8lEApqaYMAFTpigxRFISXF3h4wP1WfCoWuzfj65dGWepFUhF043w40d06oRVq/6cr162DO7uv525dw/TpuHgQRgYoGJF1K2L7dv5VKqJsK27lCP69IGFBRYtAgCxGLt24cEDLFnCUtKECXj0CCdOKKv9RYtQrRratVNW+3keIaGMCqHsVG9KJZe5RpOTqWlT8vD48/yFC2RlRb+m04uKovLlaceO9DOXL/95Td6Dba7Re/dizcwoJYVV/Bzz7h2ZmND9++mHpUtTQABLSadOkaUlJSbKdXGOco2+fEmGhvTypeLa8gC5yTUaFPSoWLF3EgmHcvhg06ZNYnFZAwMriZpIV+lco6rAgwcPRoz4UbgwZs36860lSzBpEn5uGCXC4MFwdPwt3aidHYoXx+HDPKnVQPz9tbt1A8NELTmlRAnMnYv+/aNCQ8NSDw8ehIsLbt5kJqlVK1hZKaUqxejRmDIFZcpw37ImsH9/QMeO4+LjO7x8+Zy1lhyQkIDhwxfIZPNiYkrNmzePtRwu0VAjnDt3ef36s/z8mvj6xv2RIeHxY4SG4tdMbUuW4N07LF78ZyOTJsFL2EOhHJycBs+f36tYMWYVuBSjV6/YBw+aNWni5eGxGICNDXx94eSEN2+YSfLxgZcX3r/nss3Dh/HiBcaN47JNjeLataiEhNI6OkViY2NZa5GLu3cxejTMzZEvX1NgLpDo7OzMWhSXaKgR3rz5MS6uhoGBdr58f+6LX7QIrq7p8/8XL2LZMuzdm8EKwM6dEROD4GDly9Uwfvz4cejQBal0krf3NNZacoZEkqyvrx0fb3nz5rfUM05OGDMGnTqB1Tde+fIYPBjTuPtFJiRg/HisXIl8+ThrU6O4cgU7dw6YMaPVoUMza9asyVpOVsTEYNs2tGyJ9u2hr48+fWBtvenhw5ticfCHD3lrxwwPo7TKQ7E5wsePydg4ZuzYjTdSiw3+wtu3VLx4elG3Dx+oVCk6eTLTptato06dchpfbWA4R1iwYDWRqKyz8xgm0XNDaGjo7Nl+pqZxFy+mnxw2jNq1YzbfGRNDpUrRtWvZXCbnHOG0adSnDzfC1B0F5givXCFjYzp7VhlyuCQkhFxcqHBhcnCgfftIIqH166l8efrwgWJjY6tXpyZNWEuUDznnCDXOCL98IUtL2ro143cnTqSJE9Nep6SQgwPNnZtVa4mJZGZGDx/mSILawNAIHR1p69Z4JqE54fx5Mjam27fTDpOTqUULmjSJmZ4dO8jGJpu1XfIYYXg4GRnRu3dcalNfcmqEwcFkYkJBQUqSwwHR0eTrSzVrkqUleXnRhw9p548cITMzCg8nIoqNjQ0MlIrFlLsqzjwhGGEGJCRQw4Y0c2bG737/ToaG9OpV2uHkyeTgkP1T/Jw5NGSI/BLUCVZGKJFQkSL04kUM/6E55PBhMjOjp0/TDr9+JSsrWruWjRiZjJo0oY0bs7pGHiNs25Z8fLgUptbkyAivXiUTEzpzRnlyckVqF7BoUerenc6cIZks/a0bN8jIKH1EIbVCfaFC5ObGRGnOEIzwT2Qy6tOHevT47d/4VxYupP79014HBlKZMvT5c/bNfv1KxYtTZKScKtQJVkZ46RLVrp2rtekqgq8vlS9P79+nHf77L5mZMesQ3L5Npqbpw/5/k60R7t1L1aqRmiyb5wP5b9GrV8nQkI4eVaocRfj2jXx9qWpVqlSJvLwy+MZ79oxKlqQTJ9LPpBrh8OFkaMinUgURjPBPpk+nRo0osy/2xEQqUYLu3CEiioggExO6elVeGWPGkLu7vBerEayMcOZMcnfPC0ZIRHPnUrVq6fZz6RIZGdGDB2zEDB9O48Zl+m7WRvjjB5mbU3CwUoSpKXLeoiEhZGJCR44oW468bN26a+LE2QcPfuvenQoVSusCZkhkJJUrR1u2/HYy1Qi/fyexWKWHeVMRjPA3tmwhCwv6+DHTCzZsoHbtiIgSEsjGhtasyYGMiAgyMqI88b39G6yMsF49OncujxghEY0fT82apT+Bpd6Knz4xUPL1K5mY0L17Gb+btRFOnEiDBilLmJoizy0aGkqmphQYyIMcuXj69GmhQg7AElNTD1/frL61fvygWrXIy+vP86lGSER16lD9+srUygWCEaZz8SKZmtKTJ5leIJORtTWdP09ENGgQ9eqVYyW9etH//pfjT6k4TIwwKooKFaLExLxjhDIZ9e9PHTumDyq6uZGdnbwJX7hl1Sqyt8/4rSyM8MEDMjZmY96qTLa3aOpwtOq4IBHduROtpVW7cOGGe/dmlfQoOZlataLhwzN466cRnjtHYjF9/64kpdwgGGEajx+TsXE2XfiDB8nWloho+3aqWFGRvt2dO2RuTnL8wtUJJkZ44AC1bUuUu/xVqkZyMrVtS87OafPTMhn17Jk+Ic0nKSlUowb5+2fwVmZGKJNR8+a0bp3StSnAtWvXjIxq1avXXsoi22HWt2hYGJma0uHDvMnJnpgYqlaNVq6URGUxV0wkk5GzM3XsmPFSwZ9GSETFitHYscpQyhlCijUA+PoVHTtiyRK0aJHVZUuWYMoU3LuHiRMREAADgxwHqlEDlpbYt09hpQJpnDmDli1Zi+CafPng749//8WMGQAgEmHLFjx+jIUL+VaipYVVqzBhQg7K9m7fjpgYDB2qTFnZ8eMHbt3C9u3w9MSQIWjdGjY2KFcOdnbbP3+eeOMG9ex54ccPlgr/4M4dtG2LdevQsSNrKf+Rmiqydm2MHq1dpEiRLK6cPBkvX2Lv3uwTHPbvj23buBTJDB48WXlk3SNM3Swxa1Y2jVy+TBYWFBVFFSvS9u2KizlxgqpXz3RJqjrCpEdoYZGWtzov9QhT+fKFrK1p2bK0w8hIKl2adu1ioKRXrwxyzWfYI/z+nUqWpL8yT3BJ794jtLTKlCvnuG4dubtTnz7UvDlVr06lSlHhwqSjQyIRAaSlRQULkrExVahADRpQ5840YgR16LBJLC6vpVWrcOFYkYjKl6c5c/jb35bZLXrnDpmZ0cGDPMmQEw8Patgw+wH51avJ2pq+fcv0gl97hHFxpKWlWmO/f6DpQ6PZbpb4SceOtGYNOTlx0MevWZNOncptI6oD/0YYEUGmpmn/ZHnPCIno7VsqWzZ9GV5YGBkaZp/zhXPevCFDQ3rx4reTGRrhqFE0cqRyxWhpVQEigSoGBmRqStbW1KQJde9OEyfSqlV07Fg2NS5+fimHhZGTExUqRGIxVa5MS5dSUpJylWd4i6qmCx48SKVKZZ8JYe9eMjdP30udIb8aIRE1aEC5qACkdDTdCKdPJzu7TDdL/OTJEzI1pXnzqF49Dv5stm0jB4fcNqI68G+E69aRs3Pa6zxphEQUHk4lS9KxY2mHx45RyZLZfPUog3nzqGvX3878bYShoWRsTF++KEuDVEpt25JIdChfvsrt2ztn/wH5uHmTunRJc0QLC1q0SFl7H/++Re/eJRMT2r1bKeEU5tEjMjGhmzezuezCBTI2TttClgV/GOG1ayQSybXlmgkabYTZbpb4yeDBNHQomZpyU1ktOZnKlKHQUA6aUgX4N8KuXdNHp/OqERLRzZtkbEyXL6cdLl1KVapQdDSvGpKSyNLytzy6fxihVEoNGmSajDD3REZSiRJUqFB6BUfOCQ5Oc0QtLapRg1av5riA6B+36OPHVLIkm7HuLPj6lSpUyF7V/fvypn/7wwiJyMiIXFxyIVGZaK4RXrxIZmZZbZb4yYcPVKQIlSiRVVrtnLJsGfXuzVlrbOHZCFNSqFix9Bw9edgIiejs2d+evkeMoLZt+c7KffgwVa6cvtT5DyNcv54aNVLWnPfp05Q/P1WuzNN83sWL1LYt5c9P2tppjsgJv96iqS64cyc3LXOFREL29jR1ajaXvXlDZcrQvn1ytfm3Ebq7k4GBohKVjIYaYepmCTmTu0+ZQqVL05w5XEr68YOKF6d//+WyTVbwbITXr1O1aumHedsIiWj3bjI3p4gIIqLkZHJwoPHj+dbQtm36/tdfjfDrVzI1zX6UTDE8PUksZrN75NAhatqUdHRIW5tsbL4UKlSncGHrh4pmzf95iz55QqVK0Y4d3AnliDFjqHXrbB6woqOpWjXy9pa3zb+NMCmJtLVpzx5FVSoTTTTC1MoSfn5yffbhw1d6elI7O+4fw93cVH1vjZzwbIRz56aX/iANMEIiWrOGLC3Tcvx//07W1vFOTqv27Mlol59yePaMjIzSeuG/GuGQIUpxZamU7O1JS4s2beK+8RyxeTOZmq4DpgHLdHWX2NrSwIG0eTO9fZuDRlJv0VQXzM2CcyWxdStVrJhVdlkiio8nO7vsu4y/8rcRElHTpr89xaoOGmeEcm6WSKV37+EikQ1Q/eVL7vMivHtHxYtntf5YXeDZCJs0+W2MWhOMkIhmziRb27QcDpMnewNTChZsGxp6O7vPccbkyTR4MNEvRnjzJpUowf2c5c9Jwbt3OW5ZMcLDwwsUqKCjU3HQoNtt21L58lSwIAGkrU2GhlSnDvXvTxs20OvXmbbw48ePp0+pVCnato1H3fJx9SoZGWVTIU4qpa5d6Z9/cjYAnqERhoWRSKSKtQc0ywjl3yyRSvHirQFfkajafeXM1A8aRPPnK6NhXuHTCGNjycCAYmPTz2iIERKRqyvZ21NiIh0+fLRIkUb58tW0sXkXEsJT9B8/qGRJun49zQilUqpbl/t1j0FBpKtLlStTjMoX17p9mxYtoi5dqHJlKlSIRCISiahQIapcmbp0oUWL0txl+PDhJUo0MDWNknMIik8iI8ncPPtiF2PGpN14OSJDIyQiU9P0Jd+qg2YZ4bRpcm2WSGXBAhKJIkuUcBo3LicjAjkhdVdGvBpXliXi1wiPHPkzB6bmGKFUSj16UOfOlJJC79+/j4mJ9fMjMzPq14+n9J5+flS7Nkkk0ri4uJUrqUkTjtfIzJ1LYrG61rWXSunyZZozhxwdycqKChYkkYi0tO4DZQAPLa2mrAX+SXw82drSokXZXDZ3LlWvrki/PzMjnDOHChbMcWvKRoOMUP7NEvTf3+T69UpWRtShAx9RlAqfRujqSgsX/nZGc4yQiJKSqHXr33IcR0eTmxuZmJC3t9JXk8pk1Lgxbd4si4iINzamR484azl1p6CWFrOKxErC3/8FYAK0FYnadeuWvhNGFejbN5uF6+vWbXFwmFamzJeflTJzRGZGKJGQtrbKzZVqihGamrY1Nn4vz2YJ+m+52oYNSpZFRESXLpGVFcf7lniGTyOsXPnP/ZcaZYRE9OMH1ajxsXJll/HjPX5+0Tx5Qq1aUc2adOWKcqOHhlLBghONjCaPHs3Z/vmoKLKwoAIFGKTO4YFbt26NGjUqIYH8/KhKFapVi3x92Q8CLV5MtWpltQIZMZkAACAASURBVCnl/v3HBQq0BrwHDJBvPcVfZGaERNSyJVWqpFirykJTjFBLy69Ro4HyXOzhQWIxbd6sbFHpNGigcpmWcgRvRvj2LRka/vnQoGlGSESTJnkBq3V0+uzb91tV6MBAKluW+vVLW1+qDAIDA0WiBsDsBg06ctLg9etUoACVLUtfv3LSniry6y16+TI5OpKJCbm50Zs3bPScPk0lS2a1uufECbKy+l6gQN1ixeoFBCiYITQLI7x/n0SinK28VTaaUn1CT2+ng0OVbC+bNAlz52LzZgwcyIOoNCZOZFBeQB05fRotWkCs9jdjbunatZmZ2dZChd6MGFG5d2/cu5d2vkMHPHwICwtUqQJPTyQncx+6YsWKYvErsXh/vXrWuW9t5Uo0bIimTRERgWLFct+eGmBnhyNHcOUKEhJQvTp69MC1a7wKCA9H377YtQvm5hm8++wZevSAqyvmzSsUE3Pt1augLl06cK6halWULImJEzlvWPnw4MnKIyQkpEaNGtleNm4caWnR/v08KPoNqZQsLSk4mO+4XMFbj/Cff2jjxj9PamCP8Nq162ZmtpUqNXn7Ntrbm0qWJAeH3+6f8HBq144qVlRKbvevX7+GcpEesG9fEotp6dLct6TqZHaLfv9O3t5UrhzVrk1+fnyUKf3xg6pUyXhRQkwMeXhQ8eLk4cFNLegseoREtGQJ6epyEIUrNGVoNNvCvGPGkJYWHTjAj6I/WbuWOndmEzr38GOEMhkZG6clWPkVjTLC8HBaupTKll0ErBaL+4wceTUqipKSyM+PrKyoUSMKDExfyRkYSOXKkaMjNwlyf5JFhXo5iYqiChVIT0+1Fo8oj6xvUamUAgPJwYHMzMjDQ4lpqaVScnSkMWMyOP9z+TGHg+pZG6FUSjo6KrROUFOGRrNmzBisWYMDB+DkxEbAgAG4cQOPH7OJrhbcuYOiRVG2LGsdvCOTITQUnp6oUwd2drh7F9OnD+jY8W6fPpYSST1LS4wZg3r18Pgx3Nzg6YlatbBtG6TStJHS2rVRuzY8PZGUxPonAaZM8TA2tjM1DZFI8Po17OxYC1IBxGJ06IAzZ3DiBN6/h5UVnJ1x/z73gaZPx/fvWLr0t5M3bqBRI6xejYAAbNsGExPu42aIWIxWrbBoEU/hOIMHT1YeWfcI+/dXiaKRs2fT0KGMNSgGPz3CRYto9OgMzufVHmFcHAUGkotLWu09Nze6fDmDfXsfPpCHBxkakqMjBQeTTEaBgVS/PpUvT97eaWNcr19Tv35kaZle1Ck3KNwj/PhRKhKVBy5qadlyoEN9yNEt+ukTeXmRuTk1akS7diX5+Pju3x+Qew3+/lS27G/7Td+8oX79qFQp8vNTSs70rHuERBQeTiIRPX/OfWgF0PSh0X79SFtbrqoiyubrVypeXBWTD2ULP0bo4ECHD2dwPo8Z4ceP5OdHjo5UuDA5OJC3d1ar+34SE0O+vmRllTbblJKStjqxTBny9k5bJR8URJUrk6MjhYcnX79+PUnRupoKGOHr19SyJYnFJBa3FYkqN2vWXbHQaooCt2hSEu3cSWXKrBGJxunqtnN3v3nxouKjprdvk7Fxesq6uDjy8iJjY3JzU2IGn2yNkIjKllWVKSGNNsI+fUhbm86d419RxoweTdOnsxaRc3gwwoQEMjCg7xklfM0bRvjgAXl5UaNGZGhI3buTn1/GP2zWpM42NWxIFhbk7U2xsXT7NvXrR6am5OFBUVGUmEjz51O+fJ0KFBhVs2ZbxaTmyAifPqWmTUkkolKl0gZdYlQ/eRrXKHyLHjlyrEiRxvr6NXv1etmoERUrRsWLk50dDR1KS5fS8eP04kX2W5A/fqQyZcj/vwztqXtsHB0zmG7nFnmMcNUq0tFRiV3UmmuEvXqRtjZduMBEUca8eEFGRmqQZfEPeDDCU6fIzi7jt9TUCKdMmdeggdPGjQ9dXalUKbKwIFdXOnOGmyLpISHUrx8VK0aurvT2LT14QP36UdGi5OpKkZFkZdVMJFqrpdWoRQvatCmbsgN/I6cR3r1L9euTSEQWFnT6tII/SN4gN7fo27dvo3/Jb/btG4WEkJ8fublR9+5kbU26umRhQY6O5OZGfn4UEpKeiVcqlYaE3GnSROLhQUQUGkqNG5ONDV26lJufRl7kMUIi0tVNr/DFEA01wg4dSFubLl5kpShTevbMQcUvFYEHI5w0iWbPzvgtdTTCt2/fGhi0AE4ZG49YvlxZ0yTh4TR8OBUtSkOG0MOHFBFBrq5UpAjVru1VpEj58ePnBQam+aWjI/n5/ZbKPAuyNcILF6hGDRKLqUYN4i0nuCqj1Fs0NpZCQ2nnTpo+nbp1oypVSE+PypcnR0cqW7ZX/vyDihRp8+kTubqSsTEfefh+ESaXEXbrRqVL8yAnGzTRCNu3Jx0dFf0TDQkhc3M+dhRxCA9GWKMGXb2a8VvqaIQBASk6Ot1MTGyPHj2Z/dW54/NnmjOHTE2pXTs6e5YiI8nAoDbwI1++mmPG0N699ORJ2qxkoULUvTsFBmZz+2VhhAcOkIUFiURUv76qLIJQBXi+RSUSCg+nw4epWDE7YLO+fjVDQ3J15b5gVtbIaYSvX5NIlE0dKB7QuO0T7drhzBlcv47atVlLyYjatVGhAvbvZ61Dlfj4Ea9fw9aWtQ6OOHsWw4drBQfv//DhZvv2rZUdztAQM2fi1Sv07AlXV7RtiwYN2urp2XTo0KF8efj7w84OU6dCTw9ubrC0xKJFMDPDsGG4cgVE8kbx9YWpKbp3h7k5Xr3CtWsoX16ZP5VA5mhrw9ISHTuiQIE44F+JRHrtGnx8ULgwa2UZYW4OS0tMmsRah3zkESNs1w5nz+LGDdSqxVpK5tjY7BswoG6nToNZC1EVgoLQrBm0tVnr4IKrV9G7NwICUKcOr3F1dNJ2p82fj0uXTiYknDh+/NiTJ6hTB2vXYt06NG+OV69w6BCePUPNmnj/HoMGoUwZjB2LsLCsWl6zBqamGD0aDRvi82dcuJBx7i4BnklJgUjUqWDBMz17tqpQgbWaLJk6FUFBSElhrUMO8oIR2tnh3Dncvo2aNVlLyZKbN/dLJGsvXLiTpApboFWAM2fQsiVrEVwQFgYnJ+zYgUaN2AgQidC+PYYNczQ07G1v37paNXz6hN274eEBd3fs24f8+WFriwIFEBMDAF++ICAA9vYoXRru7njxIq2d5ORkmQwzZ8LAABMmoFMnxMQgIEBT8oWqBZMno1Ilj2/frvr5/Y+1lmwYOBD58mHJEtY65EBE8o+SqB6hoaF16iwRi1s8eTLU0pK1muy4ceNWv37zEhJavn49WiRirUYOJBKJVCrV1dVVUvvm5jh/Hpk91cbExBgYGCgpNIc8fowWLbBmDTp3Zi0FICLRX/fWt2948QIREYiISHvx/DnevIGuLmQyJCQAQIECd2Jj2wMGYrG3nl6bMWMwd24e6awrD/5v0V27MGMGbt1C8eJ8hv2NuLg4PT09sXw58vv2xfnzePdO2aIyRSaTSaXSfPnyZX1ZHjDCKLF4tETyRC1qFxChTh3MnKkSX5rZolQjfPgQjo6IiMj0ArUwwogING2KBQvQty9rKTlEKsW7d3jxAuHhuHwZAQF94uO/AD2BPbNnn27bFrVrC/VAsoHnW/TuXbRqhbNnUbUqbzEzIEdG+PEjzMxw6xazpRtyGqHa3+n58k2ysBjVpAmeP2ctRQ5EIsyejZkzIZOxlsKaM2fQqhVrEbkjMhItW2LaNPVzQQBaWihdGs2awcUF27cjKmqLWPwImO3i0iE+HiNGwMQEPXpg/Xq8f89aqwDw7Ru6dsXKlYxdMKeYmKByZUyZwlpHdmQ89nHr1q2zZ8/eu3fvy5cv2traxsbG9erVa9myZQXVm5ytVk0rJGTMhg1o0AAzZsDVFSo+6ujoiPnz4e+PHj1YS2HKmTO81obknC9f0LIlhg7F8OGspXCBjo6ORPIqMTGxQIECqWc+fsSpUzh6FFOnwsICDg5wcMg7i5vUC6kUffqge3e1/NKYPRu9eiE5GTo6rKVkwa97KZKTk9esWWNtbQ1AS0urdOnSNWvWrFq1qpmZGQCRSNSsWbPDGeaFZMSv+wgfPqQ6dah1a2bloeXn9GmysuIm1YhSUd4+wqQkKlw4m9rlqryP8Pt3ql2bZs5krYNTMttHmJJCISHk4UG1a6clivP1VcvcuZzD2y06aRI5OPC3ZT5r5NxH+Cv6+jRrlpLkZEOO9xHev3+/atWqM2fObN68+dmzZ6Ojo1+9ehUWFnb//v3IyMiPHz8eOHDAxMSkZ8+ezZs3j46OZmbdmWBtjWvX0Lw5bGywfj1rNVnSsiXMzLB7N2sd7Lh6FRUrqutaxPh4dOiABg0wZw5rKbygpZVW7ykkBA8ewNERQUGoWhV16qSvj09OTr5//z6p84IDlSUgAAcOYPduaGmxlqIoPXrA15e1iKz5aYnHjh1btmxZfHx81s754cOH8ePHRyg7sat8ZJhr9N49qlmTunX7rTSJqnH5MlWooOqJZpTXI5w2Lfss5KrZI0xKorZtqX9/lUgozC05SrqdlERnz9LkyVS1KhkZUaFCDQsV6tOvn6tSFaoaPNyid++SoSGFhio7Tg5QoEf49SuJxRQcrCRFWZHjHmG7du0mTJigp6eXtXGamJgsX768rArXUa1WDTdvokoV1KiBgwdZq8kEOzuULYtt21jrYISa7iCUStG3L/T0sHGjpq+o1NGBvT0WL8b9+wgJISA2JqZxUNA7iYS1sjxEVFTaAhkbG9ZSckexYqhWDVOnstaROb/9NS9fvvzp06espHBIvnzw9ERAANzd0aMHvn1jLSgjFizA3LkqUV6cZ6Ki8PQp6tdnrSOHEMHFBd+/Y9cuYcHIb5QuLQoO3jl/vqRWrTXNm+PjR9aC8gQyGfr0QZcu6NWLtRQumD8fwcGIj2etIxN+M8IVK1ZUqlSpTp0669evj0lNQaHO1K+PsDBYWKBaNRw9ylrNX9jaompVbNrEWgfvBAWhcWPkz89aR04gwsiRePoUAQFqppwfqlat6u4++uhR4w4dULcuQkNZC1J/pk9HcjIWLGCtgyPat4e29lozs3aXLl1irSUDfjPCy5cve3l5ffv2bdiwYcbGxj169AgKCiJ1ngDX04OXF/buxfjxGDYMsbGsBf3OnDmYPz8ttYfmoI7jou7uuH0bJ06gYEHWUlQYkQhubli6FO3a4cAB1mrUmUOHsGcPdu/OO2MPEREREsnaHz8G9eo1nrWWDPjNCM3Nzd3c3J4/f3758mVnZ+fjx4+3bNmyTJkyU6dOffnyJSOFHGBnl/aIWr06LlxgLOZXbGxQrx7WrWOtg1+CgtTMCOfMwfHjOHECKp/oRiXo3h2nT2PyZEydKiSOUIQnTzBsGPz9YWTEWgp3lCxZUlc3FlhsbW3PWktGZLGQJjo6et26dQ0aNACgpaXVpk2bffv2cbaahwsyXDWaBSdOUKlS5OJCci+OUzoPHpCJCankAkmlrBoND6eSJeW6UkVWja5cSZaWGrFtLkerRrPl82dq1owcHen7d66aVDmUcYt++0YVKtCWLZw3zBkKrBr9iY1NTN263MrJBg7qERYuXHjYsGFXr159/PjxiBEjTp8+3UMdExv8Qps2uH8fCQmoVu1BtWqdJ05kvxGsShU0b441a1jr4Av1Ghf188OSJTh9GmZmrKWoG4aGOHUKJiZo3DirjLICvyKToV8/dOiAAQNYS1EOy5frh4RA9XahZ5drlIguXbq0ePFiPz8/mUxWXv2LchYpgm3bYGy85sGD0Rs3Xn7HMC/6f8yZg//9Dz9+sNbBC2pkhKmrjk+dggrvFVJpdHSwcSNGjULDhjh/nrUadWDWLHz/jkWLWOtQGk2bomhRTJvGWsdfZGqEb9++XbRokZWVVdOmTXfv3t2mTZszZ848e/aMT3HKY9asLiYm0xMTi3/7ZspaCywt0aYNvL1Z61A+UikuXYK9Ss4R/MHp0xg5EseOoVIl1lLUHBcX7N2LPn2wYgVrKarN4cPYsQMHDiC7SgnqzZAh2LmTtYi/+WOoNCEhYd++fY6OjlpaWgBq167t7e39NeukkOzI6RzhH2zfTlZWFB3NoSIFiYggI6Nscm/yD+dzhMHBVLOmvBezmiNcsWKToWEdPb1JN28yic8MbucI/+D5c6pShVxcVD2bUo7g8BZ98oSMjUktbrnczBESUUICaWnRgQMcKsoKReYIJ0yYYGJi0qNHj5CQkAkTJjx+/DgkJGTs2LHF1DQpZHb07YvmzdG/P5jvEClbFp07Y/lyxjKUjVqMi/r47PjyxV9H55StLWspeYjy5XH9Oj59gr09Pn1irUbFiImBkxMWLoQm3HK6urCzg6cnax2/85sRHj58uG7duvv27Xv9+vXixYsracCo0MqV+PIFS5aw1gHMmIG1a/P4d4TqG+HXr0hImFG27Ig5c8ay1pLX0NdHQABatICtLW7fZq1GZSDCgAFo0QKDBrGWwhf/+x8ePFCtOpe/bde8fft24cKFWUlhQr582L8ftraoUQOtW7NUUro0evfGkiUq4crKICYG9+7Bzo61jsyRSNC1K4YOtff0VIdpTDVEJIKnJypXRtu2WLsWTk6sBakAs2fjwwfNqkVTqxZMTeHmpkLJln/rEf7qgi9fvpw2bVqnTp3atWuXemb37t0nTpzgVR0vmJlh+3Y4O7Nf5D1tGjZvhgqsY1UK586hfn1kl9SdJSNHwtAQHh6sdeR1evbEqVOYOFHYcY8jR7BxI/bvV+2itUpg5Ej4+7MW8QsZrxq9fv16jRo1Nm7c+O7du3v37qWefPPmzfjxqpgdJ/c0b47Jk9G1K+NsZ2ZmGDgwzy6eVvFxUS8v3L4NPz+IRKylaAA1a+LmTVy/jl69EBfHWg0jwsMxdCgCAlCiBGspvDNtGlJSsGULax3/kbERjhgxol69ev/+++/SpUt/nmzfvv3Tp08/5tHc8hMnwtISY1lPDE2dil27oM757DJFlY3w+HGsXo3Dh4VUovxhZIRTp2BggFq1brdqNTQgQPXy4iuNt2/fNmzoZGvr4umZVLcuazUsEIvh4AAvL9Y6/iMDI/z27dudO3dmz55tYGAg+uXxuEyZMgAiIyP5U8cjIhE2bcLVq9i4kaUMQ0O4uKjQ/cEVr18jKgrVq7PWkRF372LAABw4gFKlWEvRMPLnx6ZNiI2dfObM0EGDVG+XtdLw89t77ZpjUpJ25crXWWthhrc3nj3Dv/+y1gEgQyNMTk4GUPCvZ+OvX78C0M4z6dD/InVV2/TpuHWLpYwpU3DwIF68YKmBc06fRsuWqljM9uNHdOqEFSugmQ/mqsDAga2LFnVNTq41erSmVGL58qWdru7W6tXf2Kh7yd1cYGWF0qUxeTJrHQAyNEJjY2MTE5PAwEAAv/YId+7cqa+vX7FixdxHTUpKevjw4cuXL2WZzJVHR0dH/Uccj3MIVlbYsAHdu+PzZ95i/kmRIhg5EnPnMhOgDFRzXDQxEV26YPDgPFL7VE2ZP3/Khw+XPnzwk0hgY5P3d1acPYt9+yo/f37p5s0jBppd0GTiRBw7phoLpjLcZr948WIdHZ1Zs2Zt2LDBzMwsLCxs6tSp+fLlmz59eu63+vv4+BQpUsTa2rpUqVI2NjZv3rz5+5rixYsXLly4aNGiRYsWHThwYGZN5TKzTGZMnUr29iSRcN6wvHz/TkZG9PgxMwGpcJVZRiolIyN6/Tpnn1J2ZhmZjPr0oe7dSSZTahy1QamZZeRk3z4yNiYPD0pJYStELhS4RV+9IjMzunBBGXJ4IpeZZf5AV5d8fLhqLAPkzCyTsRHKZLJU5/vpl2Kx2MXFRZ4Ws+XSpUsfP34kopSUFCcnJ2dn57+vKV68eHh4eLZNKckIpVJq04bc3DhvOAcsWED//MNSAHFnhCEhZG2d408p2wjnzKHatVWoIBdzVMEIiejVK2rWjBo0oOfPWUvJjpzeogkJVKcOeXsrSQ5PcGuE3bpR6dJcNZYBuSrDJBKJFi5c+OrVq507dy5dunTLli3h4eG+vr75uEgH27hxY2NjYwBaWloNGzbMbPVNbGzsZ0YDlGIxdu7E/v3Yv59JfABwdcX58/hv64p6kzpBqFIEBGD9ehw+jAIFWEsR+J3SpXHuHHr2RIMGWL+etRpOGTVKJZamqxTLl+PNGzx4wFhGVitfzMzMevfurbzYSUlJ27dvHzx4cIbvOjo6JiYmFilSZMOGDfaZFyxISEgICgr6eWhvby/mYklGsWI4cACtW6NKFVhb5769HFOwICZPxpw5qrXtVDHOnMHEiaxF/EJYGIYNw4kTKFmStRSBjBCJMHYsWrRAnz4ICsK6dcgD2Y7XrMGtW7h2jbUOFcPcHJaWmDQJJ0+ylCGi/xJOf/r0qWjRovL0+aKionR0dP5eVpojZDLZoEGDIiMjT5w4kVrp4lfevHljbm5ORD4+PnPnzn3x4kWGud9CQ0Pbt29v/YtT+fv7c9JtTWXPnnxLluhcuBBvYMAgLXdiImrV0t+5M8HGRsp/dAASiUQqlerq6uamkYQEUfnyBZ89iytYMGe/w9jYWH19/dyEzpAPH0QtWhRYtCjJ0TGF88bVGplMlpSUpKdKuX8SE7FwYX5/f+21axObNGHzV5AF8t+it25p/fOP3unT8RYWqrAyJFfEx8fr6upy0t9Ixd9fe+hQvY8fY5SRXkcmk+XPnz/7L7Gfg6Tbt28vV66cr69vdOZ1iV6/fj1r1qwiRYrIM4GXBTKZbMSIEXZ2drGxsdleXKxYsUuXLmX4lpLmCH/FxYU6dWK2nmLVKnJ0ZBOaOJojPH6cmjZV5IPKmCOMj6d69WjhQs4bzguoyBzh35w+TaVKkasrJSaylvI7ct6iHz6QuTkdPapsOTzB7RxhKgUL0pw53DaZhpxzhOlDoz179oyKipo2bdq4ceMcHBzq169vZWVVtGhRqVT69evX+/fvX7lyJTg42NLS0s/Pz9LSMjcuPXXq1LCwsFOnTmXbrfz+/XtcXBzDVOCrVqF5cyxdyma/y9ChWLYM16+jfn0G0TlBdTZOEGHwYFhZYepU1lIEckLLlggLw9ChqFMHO3agRg3WgnJCSgp69sTgwWjfnrUUFaZrV6xZg5kz2Sn4wxhjYmJ8fX0bNWr0x8Z5fX19R0fHQ4cOpeR6XfPChQt1dXUXL17s6+vr6+u7d+/e1PP9+/dfvHgxEV2/fn327NlHjx49cOCAnZ1dkyZNMnsA4aFHSESvX5OpKZ06pew4GbN+PbVsySY0Jz3CqlXpxg1FPsh5j3DWLGrYUOV6FaqDyvYIf+LnR4aG5OVFXHdIFESeW3TsWGrVSj12g8iJMnqEnz+TWEzBwdy2SqRAj/Cn4bm4uLi4uMTFxT1+/PjTp0/58uUzNTWtXLkyVzlljIyMnJ2dnz9/nnpoamqa+qJ+/fqlSpVKPfP169e1a9fmz5+/Q4cOo0aN4nA8WgHMzbFrF3r1wo0bKFuW7+iDBmHpUly8iKZN+Q6dez58wLt3qF2btQ7A3x/bt+PaNeTPz1qKgKI4O6NJE/Trh6AgbN2qBmuddu9GYCBu3cJfqyAEfsPQENWqwd0dFy8yUsC9BfMIPz3CVBYvplq1KD6en2i/4edHdnaU+754Tsl9j9DPj7p1U/CzHPYIQ0LIxITu3uWqvbyJ6vcIU5FIyMuLjIxo1y7GSrK+Re/dI0NDun2bNzk8oYweIREdPUpiMff7enO1j1DgbyZNgqUlxo1jELp9++ibN+2MjetfunSFQfhcoAoThJGRcHLC+vUqmvJbIKdoa8PNDUePwsMDffsmb9689xbb7MAZER0NJycsXYpatVhLURPat4eBAbNqoIIRyktqeYrgYGzaxHfoFy+e6eiU+vZtwPHjl/mOnQuIcPYsHBxYakhIQOfOcHVFx44sZQhwTt26uH0bjx79b8iQay1ajHn58jVrRenIZOjbF46O6N+ftRS1om9fZhUKBSPMAfr62Ls3ecSIDkZGtseOneItbp06dUaPrlG27KO3bzNOPqCaPHiAAgVgYcFMABEGDkSVKqq1nV+AK/T10b9/kYIF3yQmJjZrpuvtjR8/WGsCAMyZg6ioPFthW3l4eSEqCqdPMwgtGGHOKFgwUk9P8uXLtK1bT/AWVCQSLVzo/uDB6tBQ4127eAubK2JiYubN22pr+4ihhhkz8O4d1q1jKEFAuYwdO+zChWlv3pwKCDC+fRvlymHYMDx5wlLS6dPYuBH790MZ28PzNvr6sLXFjBkMQgtGmDPKli3r6mpvZRVw796Y6GheQxcsiH37MG4cwsN5jasYgwZN8vePPHGiV2p5S5559uzZ6tWvdu9GQICwTDSPU7t2bRMTExsbbNuGhw9hZoYmTdChA35JvMgfL1+if3/s2YMSJRhEzwMsWoTQUPD81QrBCBVg7twpT59ub926vJMTeP6Sr1YN06ejTx++4ypEMSDcwEDr7/x5yubcuXN16gx1de25ZEmYkRHPwQVYYmoKT0+8eYPu3TF+PGrVwvr1/BX7TUiAkxNmzYKdHU8R8x5Nm6JYMbi78x1XMEIF8fZGkSIYNAjEbxZSV1eULMk0BYN8REUtcHcf9fDhRf6N8OrV6JgYswIFjIoX/85zaAFVIH9+ODvj/n2sXIkjR2BhgalTkUmRGy4ZMQJVqmDECKUHytsMGQL+J4AEI1SQ1FJNL15gzhxe44pE2LIFe/fi6FFe4+aI48fx7p3I09O2UKFCPIe+eBGrVnXx8Phn//7RzZo14zm6gEphZ4cjR3DpEhISULUqnJ2VWNfM2xt37sDXV1ntaw6zZyM+Hnv38hpUMELF0dNDYCB27cLWrbzGLVoU27fDxQXv3/MaV05kMkyfjsWLwVEmohxw4gS6dcPu3SIPj45t2rTmO7yAKIbHnQAAIABJREFUSmJpCR8fRESgdm04OsLODvv3Q8ppKYurV+HlhYAAobwlB+jooFEjzJ/Pa9B0I/z69evDhw95Da7+GBri+HFMn873kt/GjeHigt69Of575gQ/PxQqhA4d+I575AgGDcKRI2jenO/QAqpP4cIYOxb//ouxY+HtjYoV4e0tW7Zsw6JFKyQSSW5a/vABPXti61aWO4XyGN7eePCA1wf9dCMMDAzs0aNH6ut69erdvn2bPxXqTPnyOHgQzs64e5fXuLNmQSzGkiW8Bs2WxER4esLLi++4e/diyBAEBqpxmQ4BHsiXD927IzgYfn44cOD05MmXpk9/1qvX3j17cPs2YmJy3KBEgu7dMXIk2rRRglxNpWZNmJlhyhT+IqaPXunr68fFxaW+fvXqVWJiIn8q1Jy6dbFyJRwdcfUqzM15CioWY9s21KmDJk3QsCFPQbPF2xv166NBA16D7twJNzcEBaFaNV7jCqgvjRph0yYLO7uHSUnahQr1O3AAz57h2TMUKoSKFWFpmfaflRUqVMhqB87YsShWTCjsxT2jRmHePGzfzlO4dCOsUaPGmzdvhg0bVqNGjYSEhMOHDz948ODvD7i4uPAkTa3o3h2vX6NtW1y5giJFeApasiQ2bcI//yAsDMWK8RQ0C6Ki4O2Ny/ymgfP1xcKFOH8euSuRKaBxWFlZRURcTklJ+bXcaVQUXrzAw4d49Ai7d+PFCzx6BF1dWFujShVYWMDCAtbWqFgR7u4z9ux5rqW19O7dUiIRw58jbzJ1Kjw9sWULBg7kI5yIfln+v2rVqvnz53/48CGLDxDP2wWyJDQ01MXFJTQ0lLWQNMaOxb17OHWK16QSrq6IjIS/P/ctSyQSqVSqq6sr5/Xjx0MqxYoV3ESPiYkxMDDI+prVq7FsGc6cQfny3ATVZGQyWWJiYgFhvcfvpKTg1SuEhyM8HM+epf0/MvJZcnIHYLKZ2fbIyAusNfJKXFycnp4eD6Xx2rdP+4XnBplMJpVK8+XLl/Vlor+NLS4urly5ctu3b69bt+7fHyhatGiudHGKqhmhTIbu3aGrix07wNtDYlISGjTA8OHgvK+eIyN8+RJ16uDhQ5iYcBM9WyNctAhbtiAoCKVKcRNRwxGMUH6iouKNjStKpUWdnGz9/XlPw88U3owwPByVKuHZs1w95spphBn8MAULFpw1a1bVqlWLZkTqNe/fvz9xgr9km+qCWIwdOxARgdmz+QuaPz/27cOMGXyv1vmDadMwbhxnLpgtnp7w88O5c4ILCjCgaNECCQkRd+7s0DQX5BMrK5Qpg8mT+YiVsauPHj26ZJbln2/fvj1+/HjlSFJvUjcX7t7Naz2RChWweDF690Z8PH9Bf+XOHVy+DH7uCCJMnIjAQFy6JGR0FGCGtrZ2uXLlWKvI40yahGPHIJMpPZCwoZ57UjcXzpiBU/xVasKAAahdGxMm8BfxVyZPhocHChZUeiAijB2LS5cQFARDQ6WHExAQYMioURCLOVt2kAWCESqF1M2F/fvjzh3+gq5ZgwsXsHs3fxFTOXECb99iwAClB5JKMWQI7tzB2bMqsUpWQEBA2XTogP/9T+lRBCNUFnXrYsMGdOyI13yVztbXx86dGDcOERE8RQSPCdWkUgwciBcvcPw4eM9gKiAgwIZly/DmjdJ7FIIRKpEOHTBuHNq146+8Vu3amDoVPXvyV6dp2zbo6ys9oVpyMnr0wNevOH4c+vrKjSUgIKA6mJvDygpubsqNIhihcpkwAQ4O6NIFSUk8RRw3DmZm8PDgI1ZiIjw8lJ5QLSkJ3btDIkFAAPT0lBtLQEBA1ahZc9Xp0/XMzOrIlLZsRjBCpbN8OYoVw8CBPFUuFImwaRN27MCZM0qP5eODevWUm+AtPh4dOqBAAaHWvICAhnLhwhbA58MH7efPnysphGCESid1c+GrV/D05CmioSF27ED//sgyR1BuiYrC//6n3GopsbFwdISJCXbsYFDUSUBAQBVYsMA1f/4BIlGTSZOslBRCQSO0t7c/e/Yst1LyMKmbC7duvV6mTPu5c715iNi0KQYPxoABStyCM3cuundXVobPGzduuLvPa948rmJF+PmB9yr3AgICqsKgQf0TE59cu7b4xAn066eUEBkb4dy5c48cOfLHyS1btlSpUiX1tZ6eXtY77gX+oHhxmJgse/16/qJF25N4mTD09ERyMpYtU0rjL19i2zbMmKGUxhMTExs27LJ6NcLDndesgfJzOQkICKg69erh5Ens3o3Ro7lvPIPvGKlU6uXlJZVKAXz79u1nDQp7e/tHjx69ffuWexWawbhxvUxNh2lr1xkyJP9/Ba+UiJYWtm3D0qUfxo1b/fTpU24bnz5dWQnVkpKwcKFYJgPwo3jxSCGvv4CAQCotWuDoUaxbx33mkAyM8MuXL/Hx8RYWFgBOnDjh5OSUer5EiRIikSjr2hQCWdC7d9f37298+OCrqwtbWzx+rPSIpUqhcOHBPj75mzX7h8Nm797FxYtKSah2/jxq1cLt2zonT56cNcv43j3lL/gREBBQH9q0we7d8PHheMlFpisQUlJSAMTExMT8V7Y5JiaGiAoJm5lzh64uNmzAtm1o3hyrVqFbN+WGq1rVLDIy+MuXQsePo107btqcNIn7hGqRkZg6FVevYsWKVJ3VGzYspy/sGRQQEPid7t0RE4OhQ1G4MGeP4xn0CI2MjIoWLbp58+b4+Pg9e/bExsYGBwcD2L59e8GCBc15K8Gep3F2xvHjcHPD2LGQSJQYyN9//ZUrYy9ePD1qFMaORUpKbhs8eRJv3nBZLTMlBT4+qF4dJUrgwQPO3FpAQCCvMmgQVq/GpElYvZqbBjMwQrFYPHHixNWrVxsYGERFRbm7u7dq1apevXoTJ050cXHRE7Y0c4SNDW7exLNncHDA+/fKiiIWi2vWrNmwoc6tW3j6NLexZDJMm8ZlQrXLl2Fjg8BAXLkCLy/IXQNYQEBAoxk+HAsWwNUVO3Zw0FrG32fTp0+vVatWREREr169ChcuLBKJwsLC/vnnnzFjxnAQU+A/ihfHsWNYvBg2Nti5E/b2SoxlaIgTJ9Ji7diBFi0UaSQ1oVrHjhzo+fYN7u44dgwLFsDZmYMGBQQENAo3N0RHo39/GBigU6fctUXqTEhIiI2NDWsVHHDuHJUoQV5eJJMpPdbZs1SiBHl4kFSazZXJyckJCQk/DxMSqHRpCg7OrQCplPz8yNSUXF3px49ML/uRxXsCSkAqlcbFxbFWoU5o5i0aGxsrzfa7g0cmTiQtLTp2LON3pVJpcnJyto0IW7RUgubNceMGDh1Cly5Kz9Btb48bN3DmDDp3RlRUDj7o44O6dXObUO32bTRsiE2bEBQEHx8YGOSqNQEBAQ1n6VIMHIhOnXDliuKNCEaoKpQqhcuXUakS6tXD/ftKj3XxIqytUasWbt6U6yOpCdUWLFA8aHQ0xo5Fx44YORIXLuC/3AwCAgICuWLDBnTpAnt7hIUp2IJghCqEtja8vODhAXt7bN7MRyxvb3TsCB+f7K+fNw/duimYUI0I27ahcmVEReHuXTg7Q9gmLyAgwCH79qFdO9Svj/8SwOQMIZOxytG7N2xs0LUrbtzAypXQ0VFirM6dYW2N7t0REoJ16zLdGvjyJfz8FLzD7t3DyJFITMThw6hbNzdiBQQEBDLl0CE0boy6dXH/PsqXz9lnhR6hKlKpEq5fR1QUGjbEy5fKjWVlhRs3oKuLOnXw8GHG18yYgbFjYWqag2bPnDlvadmkbt0ZrVqhe3fcvCm4oICAgHK5fBnVq6NmzRxvEsvGCJOSki5evBgbG6u4NAGFMDDAvn1wcUHDhjh1SrmxUpPduLvD3h579vz57t27OHsW48bJ1dTLl9i/H1OmoEePdc+fr7t//0xISPzYsULibAEBAT64ehVlyqBSJXz6lINPZfP99PHjx2bNmj179ixX0gQUxcUFe/di8GCpjc20du0GREb+n707j6s5+/8A/rrt+6ZSQkJS9n1XTGGYUMm+FtkmZIvwzQxG2cbYyzLIWmgURQhZsmUbWRqyhWSptKruPb8/rh/Rbbndz3JvnedjHvPIp8/nnHd86t3nfM55n9fs9TV6NE6exKJFmDgRBQXfjs+ejd9+K3V6Z0YGTp/G4sVwdoaZGTp3xu7d0NLCtGme9epNGDGie+3aWuzFTFEUVZySEu7ehZkZ7OykmIFP3xHKu27dsHVrQv/+L27dchg5cufmzX42Nmz11bIlEhIwfjw6d0ZYGGrXRkyM0suX8PD4ds6nT7h7FwkJX/579QpNm6JNG7i7Y906WFl9PbHXb7/1YitQiqKoUigp4d9/YW0NW1v89x/U1cu/hCZCBdC1q62d3ds3b7YaG29wdISODgYOxMCBaN+e+emXeno4eBDr1qFz54L8/LZZWX03bJh15YrJ18yXkoJmzdCmDRwd4esLOzs6BZSiKPmipoYHD1C/PmrW3BQW1rZv345ln08ToQLQ1dW9c+fLhkSEICEB//wDT0+kp6N/f7i4wMGBycmlnz+jQweoq7ulptoCBVOn+nbosKNtWzg5Yf582NjQF34URck7LS0EBcUOHBi2ZcuIckv5l5MI69Spk52dTQttyw+BAG3bom1bLF2Kx48RHo7ff8fQoejTBwMH4uefK1OrpaAAd+/ixg3cuIGEBCQlwc4OFhaOz5+vBFTHjnXavp2Fr4SiKIpN7dvbqqg86dnzPtCl7DMFhBBuYmJDQkKCl5dXQkIC34Hw7MMHHD+OY8cQE4OmTeHuDnd31KpV6vlFRXj06Nt7vjt3ULcu2rT59p/4N58nT55kZWW1bNmSsy+kuKysLF1agY1DIpEoPz9fS4tObqqo6nmL5uTkaGpqKinIuFB+fn5eXp6hoWHZp9Gh0aqgRg2MHo3Ro5GbizNnEBaG33+HlRV++QXNmt0+cCBo1KiBDRr0lpj53N3RujUk/vSrW7euUCjk/KuhKIpihpqamrKycrmn0URYpWhpwdkZzs4oLMS5c/jnHwwbNrOwcNnhwxNtbXu3aYO2bTFkCFq2lJz5KIqiqiGaCKsmVVU4OcHJCVpaP/3996w2bVqyvSqfoihKQSnGOC9VaStXLkhJiT15cjffgVAURcmpCiXCgoKCx48fFxUVsR0NxQYNDQ2+Q6AoipJfkhPhnDlz/vzzT/HH9+7dq1evnrW1tbm5eXx8PIexURRFURTrJCRCkUi0ZcuW+vXri/84a9YsVVXVv//+u3Xr1l5eXgq93IKiKIqifiBhskx6enp2dnbjxo0BZGZmnj17dsOGDWPHjrW3t69fv/6rV69q167NeZwURVEUxQoJT4TiZz7xeskzZ84UFhb26tULgIWFBYDU1FRuI6QoiqIoFklIhDVq1DA0NIyMjASwc+dOGxubevXqAXj16hWAcpfoUxRFUZQCkZAIBQLBrFmzZs+eXbNmzcjISG9vb/Hx2NhYPT09S0tLbiOkKIqiKBZJXlC/YMGCJk2a3Lx5s3Xr1gMHDhQfzMvLW7RokYoKXYNPURRFVR2lZrWBAwd+TYFiv/76K/vxUBRFURSnSk2Eubm54eHhiYmJQqEwMDAQwPXr13V1dcWzSSmKoiiqapCcCJ8+ferk5JScnGxkZKShoSFOhOHh4adOnbp+/Tq3EVIURVEUiyRXlhk/frympmZSUtKhQ4e+HnRzc0tISEhPT+cqNoqiKIpinYREmJWVde7cuTVr1jRs2FAgEHw93qBBA0JISkoKh+FRFEVRFLskJ0KRSCRePl9cXl4eAFp6m6IoiqpKJCRCU1NTPT29CxcuACj+RBgVFaWqqmptbc1ddBRFURTFMgmTZVRUVMaNGzd//nxDQ0MDAwMAubm5R44cmTNnzsiRI3V0dDgPkqIoiqLYInnWaEBAwNOnT4cMGSKuOKqvr19UVNStW7evezNRFEVRVNUgORFqaGgcPXr0/PnzMTEx79+/19HRcXBw6NevnzgvUhRFUVSVUVa9NHt7e3t7e5Y6Fu9xUfwdZEkikYimXoqiKIpVPKQZkUg0bdo0fX19fX19b29vkUhU8pzQ0FBzc3M9PT17e3u6YIOiKIpij+REaGZmJiiF7F3u2bPn1KlTz549e/78+dmzZ3fv3v3DCW/fvvX09AwNDc3MzGzWrNnX7S8oiqIoinGSh0YXLlyYk5Pz9Y8ZGRlxcXGJiYk+Pj6yd7lz584pU6YYGRkBmDp16s6dO8eOHVv8hP3793fs2LFbt24A5s+fX69evQ8fPtSoUUP2rimKoijqB5ITocSNJry9ve/fvy97l//991+TJk3EH9vZ2T1+/LiMEywsLLS1tZ8+fVpaIrx1S/DLL0NDQjaL/0j3DaaoSrh582abNv0BQc+ejc+cOcV3OBTFKSk2F5w+fXqjRo3WrFlTsuiMVNLT078uRtTV1f3w4cMPJ2RkZNSrV+/rH/X09D5+/Fhaa4RsOn7cu379+uI/Pn78WF1dXZbwqK8KCwuFQiFftYRycnIYGYqnvsrJEbx5I0hLE6SmCt6+FaSlCd68Ebx/r/T6teDRo21AG2BkbOyZkSOLfv21oGlTCS/vqeKq5y2am5srFAoVZRqjSCSqSEaQIhEqKysTQt6/fy9jIjQxMcnMzBR/nJGRYWpq+sMJxsbGnz59+vrH9PR0ExOT0loTCPyBAC2tHmfOgO4QxSxxItTQ0OCld0IIrd4gFT+/wB07wpycBgwduigtDa9eIS0Nb94gNRWpqXj9GgIBLCxQsybMzGBujpo10aQJTE1RUIAZM357/rw1kFiv3rL4eJV9+1TU1NCqFUaNwqRJUJAfelyrnreoQCDQ1NRUoEQoFArLPa1CiVAoFCYnJ8+cOVNdXb1Ro0YyRta4cePbt287OTkBuH37dskNDm1tbffu3Sv++MmTJ4WFhV8f+Epq1SqtVase4eFo0gSzZmHFChmjoyjFk5aG0FCsWXPo8+fIgwf7pacvMjWFhQWsrdG9O8zMULMmatWCtraEa48excSJqFPHZMyYF82bFyxcqP7vvxCJsG0bQkIwcya8vVGvHgYMgK8vatbk/GujKA4QSWpKut81NDQ2bdok8XypHDlypE6dOomJiQ8ePLC0tDx06JD4uKur682bNwkh6enp+vr6+/bte/fu3ZAhQ0aNGlVaUzdu3GjdunVBAenZk/z8M1FTI7Vrk6Qk2WOkCCGkoKAgLy+Pr94/ffrEV9eKIjeXhIaSX34henrE3Z34+ka0adNv374jFby8sJD4+pK6dcmGDaR+fZKTI8zJyenZk2zf/t1p584RFxdiZEQEAmJkRFxcyLVrzH8tiqh63qLZ2dlCoZDvKCpKKBQWFBSUe1qFZo0qKyvXqlXLwcGhVq1asqdeFxeXx48fDxw4kBAyZcoUNzc38fGcnBzxM6yBgcHRo0d9fX1nz57ds2fPv/76q+wGVVURFoaOHREYiIMHYWuLefOwdKnskVKUPPr8GTExCAtDRATatcOoUdi/Hzo6AJwB5wo28u4dhg2DmhquXkWPHvjrL2hoID8fy5fD1RVDh0JL68uZ9vYQ19VIS8P69Th8GB07QlX1y8CplxdUpHjBQlFyif2UzCLxE6H444cPiakpiY0lmzcTVVVSvz55/Zrf6BQefSKUK0IhuXCBTJtGatYkXbqQtWtJWlolm7p+ndSrR3x9iVBIli0j/fuL2xfm5OQQQtzcSEBAWZd//kxCQoi9PdHSIkpKpH59YmIyVFm53tSpvpUMSDFVz1u0Sj4RVp1ESAg5f56Ym5OkJPL6NbG1JSoqZNUqHqNTeDQRyol794i/P7GyInZ2xN+fJCfL1FpQEDE2JkeOEELIixfE2Jg8eUJIsUT46BExNibv31eotdOnycCBQsAOeKms3ObzZ5liUyzV8xatkonw26DGpUuX/v7773KfILdt28bmA6pMunfH0qXo3x+XL+P+ffzxB3x9sX07zp1DiampFCXvnj/HgQPYuRN5eRg6FFFRss6Lzs/H1KlISMCVK2jQAABmzIC3N36Yi9aoEdzdsXw5Vq0qv82ffsJPPyl17944Pr6nmtoSAwOsWYNJk2SKk6I49i0Rvn37Ni4ujsdQGOHhgX//xZAhiIqCnx9GjULv3qhdG6tXg1Zqo+RcWlralCn/MzU1b9FiUUiI0qNHcHXF1q3o0gWyL1d78QJubrCxweXLX97/nTqFO3fw/xO0v7N4MezsMGXKjzmyNHFxh8UfLFoEb2+sXInISNjZyRozRXGEg4dT9vwwNComFJL+/Ymn57cj/v5EWZm0aEEyMzkNT9HRoVGOzZ4dIBBsEAhG9Ot3+dQpUlTEWMvHjxMTk+/e/H3+TGxsyLFj3458HRoV+9//SOnztcvy9i2xtydKSsTFhVTtkdJqeIuSKjo0qhiLIqWipIQ9e3D9Otat+3Jk8WLcu4d372Bqij17eA2Ookpx/Dh27nTQ1t7VqNGLPXtsHR2hrMxAs4QgMBCTJuGff+Dr++34ypWwtUW/fqVeOGcOTp/GrVtS92hqinPncPw44uJgYICtWysTNkVxqZyJzx8+fPhaBUasjLXt8kNXFxER6NwZ9evjl18AoHFjvHqFmTMxZgw2bsSpU6h+FSEoOZWXh3nzEBGBQ4c6dOt2hcGaHZ8+YcwYfPiAa9dgZvbt+MuXWLsW166Vda2ODvz84OeH6OjKdN2nD96/x6JFmDIFAQE4fpwWfqLkl+Rvudzc3MmTJxsYGBgbGzf4HsfxVZqlJY4exfjx+PffbwfXrMGtW3j2DDVqJOrpderZczB/AVIUAFy7hpYt8eYNbt6EvT0YzIK3b6NVK5ia4syZ77IggOnTMX06rKzKaWHiRDx+jDNnKh/DkiVIToauLpo0gasrCgoq3xRFsUjigKmnp6eWltbixYvt7e1dXV3Xr1/fp08fHR2d9evXMz2EKxOJ7wiLCw0lVlYkNfW7g0Ih0dbuC0QIBO3FtWwoieg7QlYVFpKAAFKr1peVDMwKCSGmpuTgQQmfOnmSNGxISv7D/vCOUOzAAdKyJZH9lVBoKNHTI1paP5atUWhV/haVqBq9IwwLC1u5cqW/v7+VlVXjxo1//fXX6OjoKVOmhIWFcZynZeTujlGj4OqK/PxvB5WU4O3dXiCYSYjx+/d0ZhvFg+RkODggNhbXr8PFhcmWi4owbx6WLMGZMxhcYsijoADTpn2pI1MRgwdDQwOyf9+7uyM9HSNHYsIENGiApCRZG6QoBklIhO/evfv06ZODgwMAVVXVr7XWZsyYERcX9/z5cy7jk93ixahXD2PGgJBvB5cv9y8qejRo0PH+/dXfv+cvOKr6IQTBwejQAe7uOHECTFQt/CIsLKJ58362tofv38fVq2jaVMI5gYFo0gR9+1a0TYEAAQGYP5+BUU0lJQQFfRkptbWFqyt42uCLon4kIRFqa2sD+Pz5MwBzc/OnT5+Kj4t3dSq5faCcEwiwfTtevMCyZd8dV1JSCgtD7dpo146nyKjq5+1b9O+P7dtx6RKmT2dgdWBxM2cu+fff7Wlpy44ehYGBhBNevMC6dVi9Wrpm7e1hY4PgYEZihKUlbt/Gnj04fRoGBti1i5lmKUoWEhKhlpZW/fr17969C6Bbt24xMTHh4eHPnz/39fVVV1dv2LAh50HKSkMD4eHYtg379//4qatXkZaGoUP5CIuqZg4fRosWaNIEFy9C5t3MJDAzG6Kt7ezl5VZafp02DT4+KLbpdUWtWIGlS1Fsk1BZDRuGjAyMGAEPD1haxhgbt5s9+3+MtU5R0pL45nDVqlWLFi0Sf9zv/5caKSkprVy5ksn3mDIrd7JMcffukZo1yZUrPx6/cIEoKZGNGxmOrQqgk2WYkplJvLxI/frk4kW2ukhMJKamZRUIPXGCNGxI8vNLPUHiZJmvRo0i//ufbCFKkpREBIK2wA2BoL7Crb6vSrdoxVXJyTLlV5YRCoWXL1/eu3dvYmIiE4ExSapESAiJiiK1a5MXL348/r//EWVlcusWk7FVATQRMuLSJdKwIfHyItnZLPbSuzcpY053fj5p1IhERZXVQtmJ8NkzUqMGKzu6dO/uJhA0VlIarKFBlixhvn32VJlbVCrVKBHml/F7ozyRNhESQlatIi1bkqysH49360YMDCTMKa/OaCKUUUEB8fcn5uYkMpLdjo4eJba2pIzv98WLyaBB5TRSdiIkhMyaRSZPrlR85RHfZoGBREODGBmRAwdY6YVxVeAWrYQqmQglL59o1qzZ6NGjT58+TYpPtawSZs1Cx44YMgRC4XfHY2Ohqoru3XkKi6pyEhPRoQPu3cPdu1/KG7GkoABz5uDPP6GqKvmEFy+wYQNWrpS1owULcOgQHjyQtZ2SNDQ0AMydi/R09O2L4cNRpw4uXmS+I4qSSHIidHNzO3HihJOTU6NGjZYuXfrixQuOw2LVunV49+6GmZmjh8esr5leRQVXr+L2bfj48BsdpdhSUlLOnIn9809hjx6YNAmHDsHYmN0e//oLNjbo3bvUE379FbNmVWaOzA8MDTFrFhYtkrWdMmhoICQEb96gWTPY26NTJ7x8yWJ3FPVFaY+KRUVFp06dcnd3V1NTU1JS6tKlS1BQUFbJIUVeVWJoVKx//0lAnL5+75SUlOLHQ0OJkhL55x+G4lNwdGhUWpmZmSYmzVVUxtep8/uzZ1z0+PYtMTEhjx6VekJ0NLG2LmuOzFflDo0SQvLySN265NIlKaOslBs3iK3tl10s5POdhSLeorKrRkOjAJSVlR0dHUNDQ1+9erV27dr8/PyJEydaWFhwmaTZM2vWCFPT+RoadczNzYsfd3eHlxcGD6a/h1KV8eyZ6ONHJSUlTReXQktLLnpcsACjR5e6GCMvD1OnYt06qKsz052GBvz9MW8eM62VrU0b3L+Pf/7BpUvQ18fMmVx0SlVTFUmq6enpW7ZssbS0rOBYmf+RAAAgAElEQVT5nKn0EyEhJD+fmJmR+/clfMrOjpiZMVBfUdHRJ0KpPHxI6tYlfn5JERERnzlZCnDrFjEzI+nppZ7g70/c3SvaWkWeCAkhRUWkSRPWp//8IDCQaGoSXV2yZQun/ZZN4W5RRlSvJ0IAQqEwKipq6NCh5ubmkydPtrS0/Pvvv7lK0KxTV8eECdi4UcKnrl9Hbi6cnTmPiVJY9+/DyQn+/li2zNrZ2VlNTY2DTmfMwG+/SS4iA+DJE6xfz8AcmR8oK2P5csydy2mBtLlz8ekThg3D1KmoUwdxcdx1TVUHkhPhgwcPFi9e3KBBg379+l26dGn69OlJSUnnz58fO3Yst+Gxa/Jk7N+P7/dbBAAtLZw+jRMnsHw5H2FRiubGDTg6YvVqeHhw1+mhQ8jIgKdnqSfMmIG5c8HGCK2zM8zMEBLCfMtlUFFBUBBSU9GsGRwc0LIlFK3sMSXHJD4n1qxZU0NDw93dPSIioqioiOmnVcbIMjQqNmwYWbNG8qfWrCFKSuTCBVmaV2x0aLQizp8nJiZcT7DKyyNWVuTUqVJP+Oefis6R+aqCQ6NiV68SCwtS4dMZ9ugRadGCKCmRn37KGDhw7MqVq3gJQ1FuUWZVyaFRyYkwNDRUIf6NZU+Ely+TevVIabm+b1+irU0yM2XpQYHRRFiuqChibFxWQmLJ0qXEza3Uz+bmEisrEh0tXZtSJUJCiKsrCQyUrgtmhYcTVdVZgB/Qet8+Hn5jVYhblHFVMhFKHhp1d3fX1dXl+NmUF506wdQUUVGSPxsZCUNDdOrEbUyUgjh6FB4eiIiAoyOn/b59i7VrERBQ6gnLl6N9e/Tpw24Yy5dj5UrwuBvNwIEYOjRHILgA5I4Y0djcHDNnIiODt3goxSUgpdSOefDgQWRk5MuXLwu+34gsKCiIk8AqJCEhwcvLKyEhQZZGQkIQEoKYGMmfffsWlpYYPhw7dsjSiUIqLCwUCoUaFdzClWlZWVny/NvY3r3w9cXx42jRguuux4xB7do/biv21ZMn6NgRt29D2rVOIpEoPz9fS0ur4pdMngwdHebn40jlzJkzjRs31tS0WLoU+/YhLQ1WVhg+HAsWVHTz4UqT81uUJTk5OZqamkpKZU20lB8ikUgoFKqWVnXpK4nPiWvXrlVWVlZRUaldu3b97zH84Cob2YdGCSGfPxNzc3L3bqknREcTJSUSEiJjP4qHDo2WZtMmYmlZ1jJ29iQkEHPzsobr+/UjK1ZUpmVph0bJ/y/n56Z0QAUlJZGRI4meHlFWJi1asPttK8+3KHuq5NCohEQoEon09fUHDBjw4cMHFgJjEiOJkBDi708mTSrrhHnziIoKPz/4eEQToUQBAcTKijx5wkPXIhHp2pXs3Cn5szduJPTsOcvK6moFvvElqEQiJIQsWkTGjKlMd2w7dozY2xNVVaKuTuztyfnzzHcht7coq6pLIkxNTQVw9epVFqJiGFOJ8O1bYmhIys777doRIyOicFumyYImwpL8/YmtLfm+MB939uwhrVuXWurByqoDcM7UtG3lGq9cIszKIubm5ObNyvXJhZCQL1NM9fSIiwtJSmKsZfm8RdlWJROhhHFeIyMjPT29TwxuRy33TE3Rr185bwEvXoRAACcnrmKi5Awh8PFBZCTi4qR+/caIvDwsWIC1a1Ha25natduoqv7euTOnLy11dDB/PhYs4LJP6Ywcidu3kZ6OadNw5w5sbGBmhpkzkZ2N+/fvyzjDgKoaJHxLqaqqLly4cOnSpdnZ2dwHxJcZM7Bhw497MxWnpoa4OFy6xFGhRUquCIUYPx43biA2lvXdJEoTEIDOndGtW6knmJpuXLnycHj4Ng6DAoBJk3D37vkBA/yePHnCcdcVp6eHJUvw5AmePkXv3vj7b+jp3WnSZEi7dp4eHtufPuU7PopXKhKPpqWlPXr0yNraukuXLjVq1Cj+KbmaNcqgNm1gbo7ISAwcWOo5dnbYtKlw4sSJq1df2rx53vjx4zgMkOJNQQFGjsT794iOho4OPzGkpGDTJty4UeoJHz/izBls21ZKvTU2qaiQrKxfIyKWvn0748qVSO4DkIqlJXbtAgBv73MbNtQnxHD3bsHOnQCgrg59fZiZoW5dNGqEZs3Qrh3s7HgNl+KE5EQYFxcnnkV969YtbuPhk7c31q8vKxECaNw4XiBIKSravGSJP02E1cHnzxg6FEIhoqJYn45fhjlz8OuvZdVL27sX/fqVWneUVQKBwMamdkLCBlvbrjx0X1nr109XVk7Jzs7dsmW0igo+fsS//+LqVfz7Lx4/xq1b2LoVeXkoLISqKnR1YWAAU1M0bIhmzdChA/r3t83KyvX1Hb6cVmJUfKWuI1QIjKwj/KqwEPXr4/hxNG9e6jlFRUUNGnR5+VKkq7vq3Tt7Tkor84auI8zJwcCBMDHBrl2l7v/Ogfh4DBmCBw+grV3qOW3bYvlymd5hV2IdYXFjx76zszOZO7fyAcintDTEx+P2bTx8iGfP8Po10tORkxMnEnkDfyopLf7wIY6X3z/4UiXXESrGF8MNVdVS96P4SkVF5fnzq9nZ11VU7Lt35yoyig8ZGXByQoMG2LOHzywoEmH6dAQGlpUFExPx9i169uQwrBJGjjQ5fJjPAFhiaooBA+Dvj/37ER+P58/x6RNyctoLBB+BSZqaw6ysMHYsLl/mO1BKBqUmwhcvXixcuNDV1bVv377iI/v27YuOjuYqMH5MnoywsPKrRmlp4do13LyJMWM4CYviXFoaHBzQti02by51liY3du2CigqGDi3rnJ07MWYMlJW5ikkSBwc8fYoXL/iMgTMaGhoi0csHD27q6EyOikLnzpg4Eba2CAzks+YcVWmSv8WvXr3avHnzoKCg58+f3717V3wwJSXFx8eHw9h4YGICZ2ds317+mQ0aICICe/Zg3Tr2w6I4dPr0uVatnFu23NmnD9atg0DAZzDZ2Vi0CGvXlhVGURH27sXIkRyGJYmKCn75BeHhPIfBJQsLsmQJfv0V48bh33+xZw+Sk2FtjcGDcfo038FRUpG4urBVq1Y//fRTZmbm2bNnLSwsxAfv3bsHIDU1lcHVjjJiakF9cTdvkrp1SWFhhU4OCCBKSiQ2ltkQ5EX1XFBvbd0DeKKt3YyX3n8wbx4ZO7accyIjSZcuDPRVuQX1xR07Rrp1YyASRSG+RX/+mSxb9u1gejoJCiJNm5LGjUlAAHn/nrfwWFJdFtR//Pjx1q1bv//+u56enqDYL6KWlpYAXr9+zVmS5kWrVqhTB0ePVuhkX1+4uuLnn/HyJcthUZxISsK7d8P19Ye4u/NfOiE5Gdu2lVpc+6tdu+RliN7REffu4c0bvuPg1rZtWLcO/z9wBgMDeHnRB0QFIyERireb0CmxYOrDhw8AVFQkr7ioSsTrKCooLAwNG6JNG3y/SweleB4+hKMj1qwZn5Fx/e+/V/MdDmbPxsyZqFWrrHM+fsTp03B35yqmMqmr4+efK/pLZJVRqxYCAjBmDAoLvzvepg2CgvDkCRwdMWMG7OwQGIiPH3mKkiqThERoampqamoaGRkJoPgT4b59+3R0dGxsbLiLjidubkhOxp07FT3/2jUIhaCTSBXa/ftwcsLSpRgnH6tDz57F7dso96X8vn28LR+UyM0NVXLuaNnGjoWVFf74Q8KnDA3h5YV79xASguRkNGiAwYNx4kTh33+HREVV8bmHikTigGlgYKCamtrixYu3bdtmbm5+584dPz8/VVVVPz8/podwZcLGO0KxpUuJp6cU5z9+TFRVyahRbMTCm+rzjvDmTWJuTsLCOOuwHEVFpEULcuhQ+We2bUtiYpjpVPZ3hISQ3Fyir0/S0hiJSN4Vv0VfvyampuT69XIuefuWBAQQE5MgJaUpurq9EhIS2A2RBVXyHaHkRCgSiebOnVt8EaJAIPD09KxIi1xiLxG+e0cMDMjbt1JcIt628K+/2AiHH9UkEd64QczMyJEj3PRWIVu2kK5diUhUzmn37hELC1JUxEynjCRCQsigQWT7dtmbUQA/3KJ79hA7O1KR75jw8Ahd3R5KSq2iopLZCo41VTIRllVZ5tWrV+fOnXv9+rWRkVH37t2tra05eUaVArOVZX7g4QFra8yfL8UlgYHw88Pp0+jRg42IuFYdKstcvAhXV2zbhv792e6qoj59gq0tIiLQpk05Z86ZAzW18mfTVJCMlWW+OnAAISE4fpyRoORayVt00CA0boylS8u/Njk5+dIlvQULjOPiUK8eSwGyokpWlpH8RKgo2HsiJITcukUsLIi0z8CDBhF1dfLiBTsxcavKPxGeP09q1mRsaJEpM2cSL6/yTyssJObm5MEDxvpl6okwK4vo6ZGPH2VvSd6VvEXT0kitWqTie7lu3kwaNpRu5Il3VfKJUPIU0OfPnwsl7UhkaGhoaGgoe5ZWCC1bon59/POPdFPywsLQtCnatEFKCqp2JVJFd/IkRo7EwYM8VyYr7tGjR/37//r0aY2kpL8BzbJPPnECVlZo3Jib0KSgowN7e0RFYcQIvkPhnIkJNm7EmDG4datCJdonTcKLF3B2RmxsWSX0KLZJfrzt0KFDA0mMjIwMDAzGjRv3oXrUEZo2TYp1FF/duEEnkcq7qCiMGYOICDnKggD27g1PShqrqqr/+nX5u77Iz/LBkqrn3FGxgQPRvDn8/St6/rJlaNoUQ4agqIjNsKgySU6Ey5cvNzQ0dHJyWrdu3YEDB9asWdO+ffvatWsHBQV5eHiEhoa6urpyHCgvXFzw8mVZm8BJpKHxpRLp6NHshEXJJjISnp44ehSdOvEdyvdatBikpra7a9ec1q1bl32mePng4MHcxCW1AQMQG4ucHL7j4MnGjQgJwcWLFTpZIEBwMFRUMG4cFHkrIAUnccDU2dl55syZxY+IRCJXV1dvb29CiLj0tjxM/GX1HaHY8uXl17iSqApMIq2S7wgPHCC1apHbt9loW1ZdupBduyp05vr1ZPhwhntn6h2hWO/ecrQchSVl3KJHjxIbG5KbW9GmcnJI585k0SJmAmNVlXxHKOGJMCsr69ixYx4eHsUPCgSCcePGHThwAECfPn1q1Kjx4MEDblI1vyZOxNGjePtW6gv79MEff8DHB2fPshAWVSn79mHmTJw4gRYt+A6lhIgIfPpU0fdq8jwuKladR0cB9O+Ptm2xYEFFz9fSwtGjCA2lRfz5ISER5ubmEkJKvgV8//59zv8Pdujr6yvzu+kLVwwN4eaGrVsrcy2tRCpXtm7FnDk4dQrNmvEdSglCIfz8sHx5hbZSSkzEmzf46Sf2w5KBiwtOnEBeHt9x8GfDBhw+jPPnK3q+sTGio7FiBQ4dYjMsShIJibBmzZo2NjY+Pj4viu0tlpiY6O/v3717dwBZWVkpKSl16tThLkxeTZ+OzZt/LCRYQbQSqZzYsgV//IHz52Fnx3cokuzeDUND9OtXoZPlYffBchkbo2VLnDrFdxz8MTDAxo3w8EB2dkUvsbJCdDSmTsWZM2xGRpUgebLMjh07Hj9+3LBhw7Zt2/bp06dFixbNmjUrKipat24dgBs3bvTo0aNt27bchsqbpk3RqFHlx3nEk0i7dWM0Jkoaq1cjMBBnzqBhQ75DkSQ/H4sXIyCgQieLdx8cNYrlmJhQzUdHAfzyC7p3l64oR7NmOHgQw4dLUeuYkp3kRNi5c+fExERfX9/atWtnZGRYW1svW7bs3r174uIyPXr0OHHihLq6Oreh8kmq/Sh+IJ5Eevs2nUTKj8BAbN2KCxdQvz7foZRiwwa0bo0uXSp0stwuHyzJxQWRkdV9LGTtWkREICZGikscHLBuHX75Bc+fsxYW9QP2p+2wiINZo2JFRcTKily7VvkWoqOJktJra+sxW7YEMxcXu6rArNGAANK4MXn1SvaW2JKRQUxNyb17FT1/0CASFMRKJMzOGhXr1ImcOMFsk3KkgrfoqVPEyopIeztv3Ejs7MiHD5UJjFXVZdYoVZKyMiZNwoYNlW+hTx9oarr/91+nyZOXZVf8pQElg0WLcOAA4uLK2dKPX4GBcHZGkyYVOlnOlw+WREdHATg6wtERc+ZId9WUKejXD337Vt/lmFz6VmItOjo6ICDAy8trxIgRLi4uH0vZQfJ8xWdBVS0TJqBhQ6Smwsyski3Y2uomJBwixCQ0VOv7xSkUkz59+jRp0oK7d9VVVZeePq1RowbfAZXuzRsEB+NW+WVkvti3D337ytHug+Vyc0P79ti0CdVgP++y/PknWrTAiRPo00eKqwID4eGBoUMRHl7d/wLZ9u2JUFVVVVtbW1ylW0tLS7sU/IXKM0NDuLsjOLjyLVy/Hn327MIlS+ImTFCqSH16qnL27z988KDJw4fCOXNOynMWBLB4MTw9UfH51/K/fPAH9eqhbl1cuMB3HHzT1kZwMCZMQHq6FFeJi84UFWHKFNYio8Q4GKVlD2fvCMUSE4m5Ofn8WdZ2duwgSkrE25uJmNikiO8IX74kzZs/1NLqYGXV6eXLl4xHxaBHj4ixsRQvgZjdfbAkNt4REkL++INMncp4q3JB2lt0yhTi4SF1Lzk5pGNH4u8v9YUsoe8Iqzs7O9jZISxM1nbGjcM//2DTJgwcyERY1P+7dAkdO2L4cJusrPjk5Mu1a9fmO6Ky+Plh9mwYGVX0fIVYPljSoEE4fBgiEd9xyIEVKxAXh6go6a7S0kJkJPbvl2mOAlW2UhPh7du3PTw8OnTo0LlzZ/GR7du3h4aGchWYnPL2xp9/MtCOszPi4hAdLXd1nxVXcDAGDcKOHfD1hZKSgO9wynH9Oq5cgbd3Rc9XoOWDP7C2hrExrlzhOw45oK2NnTsxcSJKmYBRKnHRmeXLceQIO5FVe5ITYUxMTIcOHS5cuGBgYPC1vkxeXp6vry+p3gXSnZ3x8uXhESN+e//+vYxNde6MW7dw9y6aNaM7sMjk82eMH4+NG3HpEnr14juaipk/H4sXo+K7wSvQ8sGS6NzRr7p0gZsbZsyQ+sL69RERgUmTKrqpBSUVyYnQ29t7wIAB9+/fn1+sKIKjo+OzZ8/evHnDVWzy6PXrlOzstfv3G/n5Bcremp0dkpPx5g0sLZGRIXt71dGrV7C3R3Y2Ll+W3yXzPzhxAikpGDtWiksUbppMcW5uOHKE7jH0xR9/4OpVHDiQJ3Hz8zK0aYMDB+Dmhrt3WQqt+pKQCN+9e5eUlDR37lxVVVWB4NsQk7i4aDVPhEZGRiYmn4Gwpk1bMtJgzZpIToaSEurXR/X+q62MixfRoQNcXLB/v8Js8C0SYcECLF8uxYR4hVs++INmzaCujps3+Y5DPmhpYcaMiyNG2FtadkhNTZXq2p49sXYt+vR5Nnfu6kePHrEUYTUkIRGKBz+VlH78lPjfTFNTk4Ow5JaWltaTJ/GuruECQcX2y6kAPT08eQIzM1hb4+FDplqt+oKD4e7+5aWgQN7fCX6zfz9UVKSbJ6VwywdLcnGho6PfZGUlAAPev2/45MkTaa8dNgwi0biVK8179VLAN8bySkIiNDExqVu37sGDBwEUfyIMDg42MjJq1KgRd9HJJWVl5SlTalRuY6bSqKnh3j20bYuWLem0gvIp4ktBscJC+Ptj+XLpMrdCj4uKubnR3YW+mTRp3JgxecrKna2tO1fi8latLDQ0Tnz+bMp4YNWWhEQoEAgWLly4cuXKiRMnxsfHFxYWnjhxwsPDY+XKlfPmzVOhFQ6AHj2Ql4fr15lsU0kJ587B2Rldu9K5YWVRxJeCX23eDBsb9OwpxSUKsftgudq2hVCIf//lOw75oKent2PH0lGjpm3YUJmhjKiokAsX5unqRhw7xnho1VVpCwz/+usvfX39r6dpaGj873//k7d1lBwvqC9u2TIycSIrLU+eTJSUSLAclOaWwwX1Fy4QCwsSEEBEIu4jklVWFjE3JzdvSnfV7NnEz4+dgL7H0oL6r2bOlKNV4YyQsS78kyekRg2SkVHJyy9eJObmJDVVlhAqo0ouqBeQ0udy5eTkXLly5f3797q6uh07djSq+NLfMsXHx2/ZsuXOnTva2touLi4+Pj4lN7t3c3P79OmT+GN7e/uFCxdKbCohIcHLyyshIYGRwKSSmgo7O7x4AR0d5htfvBhLluC331DK182RwsJCoVCooaHBS+9ZWVm6urrFjwQHw98fu3fDyYmXiGS1eDGSk7F7txSXFBWhbl3ExnKxcEIkEuXn52tVfEmHlC5fxsSJVeqhsOQtKq0RI9C8OXx9K3m5ry8ePcI//8gSgtRycnI0NTVLTiKRTyKRSCgUikuHlqGscU5tbe2fWBiRiY+Pb9eunY+PT3p6+oQJEwoKCvz8/H445/z586tXr7awsABQs2ZNxmOQnZkZunVDWBjGjWO+8cWLUbcuJkxAWhrWrWO+fYXz+TOmTMGNG7h0SfGGQ8XevcOGDbh2TbqrTp5U4OWDP+jYER8/4uHDKvLlMGLBAvz0E6ZNQ+XmIP7+Ozp0wM6d0i3FoSRg/9m0LH/++WePHj1KHq9Ro0ZSUlK5l/M4NEoIOXqUdOnCYvsREURZmQwYwGIXZZOTodGUFNKhAxkyhGRn8xULA379lcyYIfVV7u5s7T5YEttDo4SQqVPJH3+w2gOnGNkys39/snFj5S9PTCTGxuS//2QPpKKq5NAoz4+3N27csLW1lfipyZMn//zzz/7+/l/HSOVN37549gz377PVvrMzzp2r7mXYFHGlYElPn2L/fhSrTlEhHz/i1CkFXj5YEi0xU9KCBQgMRGFhJS+3s4OvL8aOhZSr86nvlPWOsNKys7Pj4+NLHm/Xrp1BscVQhw8fnjRp0r///mtWYou/P//8s0WLFvn5+atWrfr8+XNcXFzJ94gAEhISevfuXXyYPiEhQU1NjaGvo3y//aZWUCBYtuwze13cv6/k4KABzDA2vnbkyIbGHI4r8fiOMDs7+8KFC6mp/ZctUwsOzu/ZU7G/yz08NGxsRL6+BVJdFRysevWq8vbt+SxF9QO23xECEAphba0dG5tXr15VKMKdnZ2tw8QcgX79NEeOLBw2rJKFFkUiODtrOjkJZ8yQ7garnNzcXA0NDQV6R6iurq6url72aawkwuTk5KlTp5Y8vnr1ajs7O/HHMTExI0eOPH78eLt27cpoKisry8TE5OrVqy1atCj52YSEhHHjxv1T7GWxlZWVgMPF1U+fokMHvHyJ8v6eZbJ2bYiPTwTg0KzZubt3Zd78osL4SoQikUhdvY5Q2FJDwy4xcaWVFcf9M+zuXfTujaQkSDuvon17LF3K3UJJDhIhgAkTYGuLmTNZ7YQjsk+WETt1CtOn4949VDq5pKSgbVucPAlJPyYZViUny/DzjjA2NtbU1PTixYsVOdnQ0DAuLk7ip/h9RyjWsycJDWW3i2fPnqmqWgoEdgLB2bZtybNn7Hb3FS/vCLOzyeLFWUBtYKOJSSuOe2dD796VeQmUmMju7oMlcfCOkBASHU06d2a7E44w8o5QrFMncuSITC3s2kXs7AgH36/0HSEz4uLi3NzcgoKC7Ozs0tPTv74C3LVrV0xMDICUlBRx5SGRSBQYGKiiotKyJTOFPdkwfjy2b2e3C0tLy4KCZ7m5CY8eORQWon599O2L7Gx2O+VeQQGCg9GoERITdf744/cePWLi4vbxHZSszp9HUhLGj5f6QgXdfbBcP/2ER4/w8iXfcciZOXOwZIlMdclHj0bTpli0iLmYqhUOcvIPvLy8DItp3769+PiQIUOWLl1KCLl8+bKJiYmhoaG2tnaLFi0uXLhQWlPy8ET4+TMxNeXuKY0Qcvw4MTUlamrEx4fdjjh7IiwoIEFBxMKCODp+W2/O4K/bfBGJSMeO5MABqS8sLCS1apEHD1iIqXTcPBESQkaPJuvXc9AP6xi8RUUi0rQpOXVKpkY+fiR165LYWIZiKkWVfCLkeflEGdLT03Nzc8s+Rx4SISFk+nQeSmYEBhINDWJkRMLC2OqCg0QoFJLQUNKgAXF0JNeuffepKpAIQ0NJixakEj80jh3jYfyQs0R49ChxcOCgH9Yxe4vu2kV69pS1kZgYUqcO+fiRiYBKUSUTofy+8DQwMFCUnS7Eo6McT1+eOxfp6ejbF0OGoEED3LvHae+yE4kQFgZbW/z1F7Ztw6lTKHPWlOIRCuHvj1WrKjMDogpU2S5D7964cwdpaXzHIWeGD8fTp7h8WaZGnJzg7AwfH4ZiqjbkNxEqkKZNYWGBmBiu+9XQQEgIHj+GqSmaN4eDA+R1yeV3CEFkJNq0wZo1WLUKFy/CwYHvmFiwdStq1oSjo9QXVr3lgz9QV0evXjh6lO845IyKCmbPRqDMG36vXIkrVxAaykRM1QZNhMzw9MS2bfx0bWWF+HicPIknT1CjBiZO5CeMCjp9Gu3bY8EC+PkhPh7OznwHxI68PCxbhoCAyly7fz9+/lmxdx8sF11ZL5GHB27ckLUcq5YW9u7FjBl0o28p0ETIjGHDcPYsn3eekxNevkRAAEJCoKvLW1Yug/jJb/p0zJ2LO3fg7s53QGxaswZduqBDh8pcW7XHRcX69cOVK/j4ke845IyGBqZNY+ChsE0beHlh7FiZpqFWKzQRMkNHB66uCAnhOYxZs/DpE4YPx6RJqFNHXvb4jY+HoyNGj8bw4bh7F+7uirShfCWkp2PdOixZIvWFOTk53t6rHz/eV4kBVcWipYWePUG30ytpyhScPIn//pO1nYULkZmJLVuYiKkaoImQMZ6e2L6d/1/BVFQQFITXr9GgATp3hrX1aXPzTosW/cFxGC9fvvT2Xrhp08nBgzFsGAYPRlISvLyq4MK4HyQmJo4aFerikm9tLfW169YFb9qUkZe3//ZtHnYW4xgdHZVIVxeTJmH1alnbUVHB3r1YvBiPHjERVlVHEyFjOnWChgYuXOA7DgCAqSnOnUNcHJ488UtN/TJj/10AACAASURBVHPp0l0uLli3jvWFzEVFuH8fBw/CwcFvwwY7b2/fbt3yxSlQpawtv6qI9PR0e/tRx4/fKCyU/nkQsLOzBWL19FJq1arFeGzypn9/nD+vGNO7ODZjBg4dwuvXsrbToAH8/TFiROUrelcfNBEyaexY1qvMSKVrV7Rvby4QeGpoWDx+jIULYWkJNTXUqQMnJ/z2G27dkrWLV69w8iRWrMDo0WjVCnp6cHHBoUOoV6+Nvv42Gxv9X39V57AKOs9UVFSys4UaGqlmZpWpxZyT06d796PPn8ebm5szHpu80dVFly6IiuI7DvlTowZGjcKffzLQ1OTJMDPDH1yPBykeVopuc4bHHeol+vABDRsiORmGhnyHUswPNfJv3cKRI7hwAQ8f4t07EAJdXVhZoUMH9O2Lfv2+PbqVLLqdlYWkJCQm4v59JCbi+nUUFMDODk2awM4ObdqgVatvmyW9f/++Ro0ala6BzlRFYy7t3Ys//kjbvPlxt26dKvGF29tj2jS4ubERWvm4Kbpd3I4diI5GGHdl5BnG3i2akoLmzZGUBGNjWZtKS0OrVjhypJJTt0qqkkW3aSJk2NCh6N4dU6bwHUeFXb+OiAhcuoRHj5CWBqHwS16sWfNsTMwoQG306OV16gy5dw///ou3b2Fri+bN0awZmjVD8+YwMWErMIVLhB8+oFkzRESgbdvKXP7wIXr2xPPnKLdQPku4T4QfPqB+fbx+rag7TbJ6i44fjzp14O/PQFP//IPZs3H7NpjYM0qREmF2dnbt2u2WLFnn7e1U9pnV4L0Ntzw9MXu2IiXCdu2+K+ny9CmOHMH58zh16i/AFbA6dOjGnDlDRo1C8+Zo0KDyO8VUedOmYeTISmZBAFu2wNOTtyzIixo10K4dYmLg4sJ3KPJn/nx07oyZM6XevaukgQMRHo5587BhAxORKY6JE69nZto8ePBTuWfSn2oMc3RETg5u3uQ7jsqyssKsWYiIwJ07K5SUjgoEW6yth8+bBzc3WFvTLFiq6GjEx1f+9/e8POzdC09PRmNSBHTuaGkaNECPHggOZqa1detw/Hg1eiO7YgV0dHDkiL2BgUqLFqfKPZ/+YGOYQIAxY+RrykzlNGrUKD//cXb2nSZNWg0YgM+f+Q5IjmVlYdIkBAdXfojv4EF06IB69ZiMSiG4uuL4ceTn8x2HXJo/H6tXM/OXo6+PkBBMnIgPHxhoTZ6JU6B4umxWltKHD6EeHj3LvYomQuZ5eODAAeTm8h0HE5SUsHs39PUxfDiKiviORl7Nm4devSpTVvSroCB5r43Hkpo10awZzpzhOw651KIFWrXC7t3MtNa1K4YOxYQJzLQmh75PgQgKkmLJFk2EzLOwQKdOOHSI7zgYoqyMPXuQnw8PD4hEfEcjf65cwdGjWLGi8i3cvYtXr9C3L3MxKZSBAws3bozLysriOxB5tGgRAgMZ+x106VI8ecJ/ASxmiURYtOhLCvTyQk6OdClQjCZCVoirzFQZamo4dAgvX2LCBP5L58iVz5/h6Yl162RaMLNlCyZMqPo1d0pz/LjniROHWrfuxXcg8qhjR9SuzdhWEurqCAkp8vKa2rTpzw8ePGCmUf6IU6CuLlat+pIC16yp5DwGmghZ4eyMx4+h+HfaN5qaOHYMSUl0q7Pv/P477Ozg6lr5FrKzcfAgPDyYi0nRFBZmA5aZmQV8ByKn5s/HsmWMDcYoKT1UVv6QmDh53ToFfjAsngInTpQpBYrRRMgKFRWMGoWdO/mOg1Ha2jh2DBcu4Lff+A5FPty9i+BgrF0rUyN798LBARYWDMWkgI4e3T5smGnPnuF8ByKn+vSBpiaOH2emtUaNGnXqJNDRWfXhwyBmWuTW1xS4ejUzKfALKTe+ly83btxo3bo131FI9vgxMTMjBQV8xyGDgoKCvLy8Hw6mpRE7OxIQwHrvnz59Yr0PGRQVkbZtyY4dsrbTujU5eZKJgGQmFApzcnJ46fr1a1KjBvn8mZfOK4+zWzQ0lHTowGSDWVmkSRMSHFyZa7Ozs4VCIZPRVMCBAweePUtZuJBoahJNTeLjQyoYglAoLKjAT2H6RMiWBg1gY4PISL7jYJqJCWJisHUrNm3iOxRe/fkndHUxdqxMjVy7howMmaabVg3m5rCzw6nyl3tVU25uyMrC2bOMNaijgyNHsGCBvGwSULaOHX8ZNmxfvXq9167NnT8fubkMPQUWQxMhi8aPr1JTZr6ysEBsLFaulMftf7nx9ClWrMDWrbJurBgUBC8vWqYAAIYMwcGDfAchr5SUMGsWli9nss1GjbBnD4YNw6tXTDbLuPx83L6tS0hrgeDzw4fpixax0gv9FmTRoEG4dg0vXvAdBwvq1kVMDBYvxoEDfIfCOULg5YV589CggUztZGYiPBzjxjEUloIbNAiRkcjL4zsOeTVqFJKSGN5qu1cvTJ6MQYPkt1xGUhIsLKCjs8XF5c3evb9bsPYunSZCFmloYOjQqjZl5itra0RHw8enCg7/lm37dmRmYvp0WdvZvRu9e8PUlImYFF/NmmjdGidP8h2HvFJVxcyZWLmS4Wb9/FC3rpwWc9i/H02aoFEjpKbqHzmyZdiwYez1RRMhuyZMwLZtEAr5joMdzZohOhoTJuDECb5D4UpqKvz8EBTEwLK/4GA5/QHEFzo6Wrbx43HpEu7dY7JNgQDbt+PmTWzezGSzspsyBSNHYsYMxMdzsac3TYTsat4cNWtW5QpSLVsiPBxjxiAuju9QODF1KiZNQqtWsrZz8SKKimBvz0RMVYWrK6KjkZPDdxzySksL3t7MPxTq6CAiAkuW4Px5hluunPx8dOqEbdtw5AjzX2xpaCJkXRWrMlNSp07Ytw/u7rhxg+9QWHb4MBIT4efHQFPi4qIyzrWpYoyN0bFjNdohoRK8vREVhWfPGG62Xj3s3o3hw5GSwnDL0nr4ELVqITkZ//2HAQO465cmQtYNH45Tp5CWxnccbPrpJ+zahV9+wa1bfIfCmsxM+Phg2zZoaMja1IcPOH4co0czEVbVQkdHy6anhwkTsGoV8y07OmL6dPTvz+d8pT170KwZbGzw6hUsLTntmiZC1unpYcCAqlbotqQ+fbBpE/r1q1KF5Yrz8cHAgejalYGm/v4b/fvDyIiBpqoYFxecPo1Pn/iOQ47NmIF9+/DmDfMtz5kDGxve3ltPnIgxYzBrFkcvBX9AEyEXPD2xdWvVL1ft6orly9GnD54+5TsUpp09i9hYLFvGQFOEYOtWOk1GMgMDdOuGY8f4jkOOmZrCweFCly6jIyMZnqImEGDHDty/j3XrmG24HJ8+oWVL7NyJiAgEBHDa9Vc0EXKha1coK+PyZb7jYN+YMZg7F7164fVrvkNhTm4uJkzA+vXQ1WWgtdhYaGigUycGmqqS6OhouRIS5jx96uvpOZ/xljU1cfgwAgKYrGJTtlu3ULs2UlPx+DH69eOo05JoIuTIuHFVfMrMV1OnYupUODgUhoWdfvv2Ld/hMGDRInTqBGdnZlqrtnvwVtCAAYiLQ0YG33HIMTe3vnp6Hrm5DsnJzDduafml4gwH4zpbtqBdOzRpgpQU1KnDendlkbUYKq/kuej2D969I0ZGJDOT7zikIbHodgW1aDFDWflXM7PmhYWFlWtBTopuX7tGatYkb98y01pqKjE0lNPbgMei2z9wcSE7d/IdRAXweIsWFRVt304sLcmzZ6y0v3o1admSSLwdmCq6PXIkUVIiy5bJ3lJZaNFt+WJsjJ49sX8/33FwpWnTIoFAOz2dFBUp8KvRoiJMnIi1axmr/7J9OwYNgp4eM61VVXR0tFzKysoeHvDxgZMTUlOZb3/mTLRqhQkTmG8ZwKdPsLHBoUOIjWVmMZLsaCLkTvv2l2bO7L1oEVdrRHm1bdvKAwc6dOlyzM9Ple9YKm/5cpibY+hQZloTibBtGx0XLZ+zM+Lj8eED33HIvenTMWIEevVi5e9q0yb89x/WrGG42YQEWFggPx/Pn8tRQQmaCLkTGbk2N3f1hg0HP8ttjVvmaGhouLm5hIfXjYlBcDDf0VTKo0fYsAFbtjDW4MmTMDZGmzaMNVhVaWmhVy+E0516K8DfH/36wcmJ+beqGho4fBirVjFWQLFjR2d19Zbt21/v1g1Pn8pXlV2aCLkzdeoIfX1PU9OO6urqfMfCET09RERg8WLuJqExRSTC+PFYvJjJd/h0mkzF0dHRilu+HD16oF8/ZGcz3HKdOggLw7hxePJEpnbev8fYsc+vXn1fULBRS2t2VJTcbT0mZ+FUaUOGDLx37+qHDxuq1V4z9evjwAEMHy7r9xLHNm2CUMhk3kpJwYULjI2yVnl9++LmTVSJScdcWLUKTZrAxQX5+Qy33KUL/Pzg6lqZGrAZGZg/H5aWMDXF8eN1NTQ0lJVHTZ/uxHCITKCJkFO1a6NdO/zzD99xcKt7dyxcCGdnZGbyHUoFpKenr169+3//e7ZtG5O/t27bhuHDoa3NWINVm4YGfv4ZR47wHYeCEAiweTNMTDBkCAoLGW7c2xvt22P06IqWBCkowKZNsLODkRE2bUKbNrh7F+/eCfLyzhYVJS9dupDh+JhAEyHXqnwNbommToWDA4YOVYAdqfr395w7971I5Gpnx1ibRUXYtg2enow1WB3Q0VGpKCtj1y4AGDYMRUUMN75hA169wooVZZ1TVIRNm9CyJTQ1MWcOGjfGzZvIzMSRI2jalOF4GEcTIdf698e9e3j8mO84OLduHQoL5WW2dGny8vDff/rKysmmpjKX1i7m2DHUq4eWLRlssurr3RuJiXj1iu84FIeqKsLCkJuL8eMhEjHZsro6wsOxcaOEvUG+5j91dcycCQMDxMYiJwdHjijSDU8TIdfU1DByZJXdtr4MKioIDUV4OLZt4zuUUmRno39/9OixNSZm6K1bpxhsmU6TqQQ1NfzyCw4f5jsOhaKmhkOH8OwZpk1juGVzcxw8CA8PHDhwOSMjA8DRo3BwgKbml/x34gTy83HunBwtiqg4mgh5MGECduxgfvhC/hkZITISCxbIyxagxaWnw8kJDRti714VB4eu2sy9zXvxAtevw82NqfaqETo6WglaWoiIwLVr8PVluOVOnaCi4jR+vE+NGp1VVNKHDIGSEiIjv+Q/J3mcBFNRNBHywMYGVlaMrc5RLDY2CAnBkCFgo0xipb15A3t72Ntj82bmJ3Zv2YIxY6ClxXCz1YGjIx4/xvPnfMehaPT0cOIETpzA778z1qZIhIgIvHnzGBgDGI8YEZmfj9hY9OnDWBc8oomQH9VzyoxYr17w9YWrK/PLnirn6VN0745Ro1jZAqawELt20WkylaSiggEDcOgQ33EoICMjnD6NgwfLmeFSEdnZCA5G06b4/XcsWLBXW3tLkybCXbuq1L7SNBHyY8gQxMWxsrumQvDxQceOGDWK4Vf6lXD/PuztMXcu5sxhpf3wcDRuDAYnoFY3dHS00kxMIK7rtHlzJVt4+hTz5sHKCpGR2LoVN27g9987v3179+7dS4xGyj+aCPmhrQ03N+zezXcc/Fm/HpmZ8PfnM4br1+HoiFWr2CouDDpNRmYODnjxojrOsmaEhQViY7FihdQz1BISMHo02rZFXh4SEhAZiS5d2AlRPtBEyJtqsm19acRTvQ8cwN69/ARw7hycnbF9OwYPZquLx49x7x4GDGCr/epAWRlubnR0tPLq1kVMDBYvxoED5Z9cUICwMHTqhKFD0aYNXrzAX3+hbl32o+QbTYS86dABOjqIi+M7Dv7UqIGICMyahatXue46MhLDhuHQIfz8M4u9bN4MT09Um8qybKGjozKytkZ0NHx8EBlZ6jnv3iEwEA0a4K+/4OeHpCRMn16NCiHRRMin6rNtfWlsbbFzJwYNwsuX3HW6Zw8mTUJUFLp2ZbGXz5+xZw88PFjsopro2hXv3uHBA77jUGTNmiE6GhMmSJis/ugRpk+HjQ0SExEdjYsX4ewMgYCPKPlDEyGfRo1CZCTS0/mOg1d9+sDbGwMGVKaqbyVs3IhFixAbi1at2O0oNBStW6NhQ3Z7qQ6UlDBoEB0dlVXLlggPx5gxOHdOlJ2dLRLh9Gk4O6NHDxgaIikJu3crQC00ltBEyCcjI/TpU422rS/N3Llo0wZjxrD+xjQwEOvX4/x52Niw2xHoNBlGDRmCffv4DkLxdeqE4OBPjo6dzMycLCwiFy7EiBF4/hyLF8PYmO/geEUTIc88PREUxHcQcmDTJrx/z+T63x8QglmzsGcPzp7l4uX/gwdITka/fqx3VE107Ij8fCQm8h2H4mvc+I2WlmFu7lAHhxtXrmDoUKiq8h2THKCJkGc//YScHNy6xXccfFNVxaFDCAmp0Nw2aQmFGD8eFy/i/HmYmzPffkmbNsHLi/6IYYxAgEGD6JQZBtjY2KxcOXTKlLQNG5iuRqrIVPgOoLoTCDBmDLZvx4YNfIfCN2NjHD4MR0c0aIB27RhrtqAAI0bg40ecOQMdHcaaLUNuLvbvR0ICF31VH0OGYMQIFscMqo+JE8fyHYLcoU+E/PPwwP79yM3lOw450KIFtm6Fiwtjm+/k5MDZGUVFiIriKAsC2L8fXbvC0pKj7qqJtm0hEuH2bb7joKoimgj5Z2GBTp3oZtxfDByIyZPRvXvYyJFeT548kaWpjAz06gUzM4SFcbeYLyTk4Pz589zdq2v1PDa5u9PRUYoVNBHKhepcg7ukGTNynz71iYho26LFtIAABAXh4EHExODaNfz3H969K38Hq+zs7KdP87p3R5cu2LkTKly9AUhNTfX2XvfunU109HKOuqxOxCvrq20xJoo99B2hXHB2xq+/4r//YG3NdyhyQF1dTVlZJBRe0NU1ysjA06fIyEB6+nf/19SEoSEMDCT8PzMzYc2aidnZomnTDqxY0YjLyI2MjEQigY7Oth49pnDZbzXRogU0NHDjBpOvkCkKNBHKCRUVjByJHTuwnD5IACoqKm/f3jt27Njo0aVu9ZKV9WNqFP8/ORkJCQ+ysrqoqX1ydHwMcJoInzxR09S8cPduds2aulz2W32IR0dpIqSYRROhvPD0hL09fv+dzrkHACMjIxcXlzJO0NWFrq7kFYEFBYMDA99qaKj34XzP0D/+gI+PgGZB9gwbhl69sHJltasBRrGKJkJ50agRrK0RFUU3K5CVmpraokWzuO83ORnR0Vi/nvueq5HGjaGvj/h4dO7MdyhUFUIny8gROmVGoQUEYOpUGBjwHUdVRzejoBhHE6EcGTwY8fF4/ZrvOCjppaTg8GF4e/MdRzUwbBjCwiAU8h0HVYXQRChHNDUxaBB27uQ7Dkp6q1bB07O6Vy7mRoMGMDPDxYt8x0FVITQRyhdPT2zbBpGI7zgoabx/jz17MGMG33FUG3R0lGIWTYTypW1bGBjg3Dm+46CksXo1hg5FrVp8x1FtDB6MQ4fKr6tAURVEE6Hc8fCgU2YUSWYmtm/H7Nl8x1GdWFnByor+vkgxhiZCuTNqFE6cwMePfMdBVcxff8HZGfXq8R1HNUNHRykG8ZAIz549O6+YzMzMkue8efNmzpw5w4cPDw4OFpX+xiwpKSkpKYnNYHmgr4+ff8bevXzHAWzevPm3337jpev79+/36tWLl66lkpODTZvg68t3HEwIDw+fMkVhKsMNHozwcBQU8BZAfHz8oEGDeOueP+PHj4+KiuI7ioqKjo4eM2ZMuafxkAgvX758/vx5w/+npPRjDIWFhfb29tnZ2W5ubhs2bFi6dGlpTQmFwjLSpOLy9MTWrXwHAeTn5+fl5fHStVAozMrK4qVrqWza9H/t3WlQU9fbAPBzJWKUYDaU1QotiLJE2QxTBwMGrFSxLlVxKuCGVkdHS/tWrK3aMtRqF0VtbbGgVFoVBTccZ5RoUoKUxa0UkUUjVDCGBiGQBGKW98OZZvJHoIjCDeT5fco595xznwTjk7udg2bMQBMGdB63/qLRaFSDZzEwFxc0cSISCEgLQKvVtrW1kbZ78qhUKg2JP0BekEajUSqV/9mMnJll/Pz8tnT/K/rcuXNWVlY//PADQRDOzs5RUVFbtmwZMWDr6JiB0FCkVqPSUhQYSHYooHvt7SglBV26RHYclgqfHY2MJDsOMPiRc43w5s2bmzZt2rNnz+PHXSzbVlRUFBISQhAEQojL5SqVypqamgGPkUwEgVasgFtmzF1aGgoMRL6+ZMdhqRYtQufPo/Z2suMAg1+/HBGqVCqFQtGpkiAIe3t7hJCbm1t4eLijo+P169d37dpVXFzs8b+LD0mlUhcXF2MvFov1+PFjb2/v53fEZDL1ej2bzTbWTJ48mRgS0/FqNOyKioSqqm0kxiCVSjUaTVlZ2cDvWqVSyWQyPp8/8LvuvYqK/3N2Ps/nV5IdyKvR1NTU1NRk5p95J1TqxzzeKRpNMvC7bm1t/fvvvwfXx/VKVFdXP3r06Pvvvyc7kF5pbm6eMmXKfzYjDP2wzGV6evr27ds7VVIolIcPH3aqXLp0KYPBOHTokGnlmjVrbG1tv/32W1y0s7M7f/78m91Msnvy5Mna2lpj0c/Pb2gkQgAAAC+Py+Xa2v7HgjD9ckS4cuXKlStX9qYlh8MpLCzsVOni4mI8CmltbW1ubjYeID5vyZIlfY4TAAAAIOEaoUwmwy/a2tpOnz4d+O8NIVlZWVKpFCG0cOHCK1eu4MuHmZmZ/v7+r3W57hwAAADw0vrl1GjPfH19rays2Gx2WVlZUFDQqVOnRo0ahRCi0WhnzpyJiIhACH344YcnTpzw8fG5detWdnZ2SEjIAAcJAADAQpCQCDs6OsrLyxUKhaurq6vJhBz19fV2dnbGxyRqamrq6+unTJlCp9MHOEIAAACWg4RECAAAAJgPq507d5IdQ19IpdIbN27U/quurm78+PFkBzVEtLW1/f777+PGjbOyssI1paWlHR0dTCbz5QdvamoqLS01/WMVFhZaW1vTaDTTZmKxODc3NygoqFN3oVDIYDCoVCpCSKfTCYVCg8HwSgKzZBKJ5M8//zR+mxobG52dnckOynw1NDRIJBL8MJgluHv37tOnT41PqRUUFGi1WuOXTiQS0en0kSNHkhdgF3Q63dWrV9lsNv6/AiFUXl7e2Ng4duzYrjsYBqeff/559OjRM/4VERFBdkRDh16vj4iI2LlzJy7euHGDzWbX1ta+ksEvX77s6OhoWuPt7X38+PFOzTIzM2NjY5/vzmAw8vPzDQaDRqNZvHjx1KlT5XL5KwnMkm3ZssXBwcH4bVq+fDnZEZm1/fv3h4WFkR3FwPnqq68iIyPx69bW1uHDh0dHR+OiVColCKKhoYG86Lq1YsWKZcuW4dePHj1is9lisbi7xuRMsfZKTJo0SUDiVINDF0EQqamp/v7+c+bM8fHxiY2N3b17t1nduNvR0REdHa1QKAQCQadDSdA34eHhx44dIzsKYI5CQ0OTk5O1Wi2FQikoKODz+QUFBXiTSCTy9PR0dHQkN8Iu7d2719fXNycnZ/78+atXr46Pj582bVp3jQdxIgT9x9XVNSkpadWqVWFhYS4uLr18KnRgKJXKOXPm2NjYXLx40XjeAwDQTwICAgiCuHXrVlBQkEgkmj17tlwur6qqmjBhgkgkCg0NJTvArtHp9NTU1JUrV1ZVVdXV1Z09e7aHxrAeIeja+vXrR44cmZaWdvjwYbOarCc2NpbFYmVnZ0MWBGAAUCiUadOmCYVChJBQKOTxeDwez1g020SIEJo1a1ZYWNi2bduOHDnS87INkAhB1+rr66urqwmCaGxsJDuW/xEUFCQWix88eEB2IABYCh6PJxKJVCqVRCLx9vYOCQkRiUSNjY0VFRXm/JC3UqksKiqysbF59OhRzy0hEYKurVmzZvXq1cnJyStWrHiFy4/Z2tp2WmhQoVCMHj269yMkJibGx8dPnz69oqLiVUUFAOhBaGioWCwWiURcLnfYsGE4EQqFwokTJzo5OZEdXbcSExO9vb1zcnLWrVvX8w96SISgC+np6RKJZMeOHevXr2cymbt27XpVI3t4eHR0dBjnkq2rq5NKpV5eXi80yM6dO9euXTtjxgzIhQAMgICAAIPBkJKSwuPxEEJMJpPNZqelpZnzedGCgoKsrKzU1NTw8PDZs2cnJCT00BhulgGd1dfXf/zxx+fOncPPBh05csTf33/u3Ll+fn4vPzibzU5MTJw7d+7atWutrKzS09PXrl1rOsFQL+GnO2bOnHnt2jV3d/eXDwyA3quqqlq3bh1+TaVS9+7dS248/Q1fJrx06VJycjKu4fF4Bw4cOHnyJLmBdUepVC5fvjwlJQU/7rl3714fH5+zZ8/Omzevy/aD9YF6GxsbLy8vT09PsgMZgiorKyMiIsLDw3GRyWRyuVyNRtOHdNWlsLAwPz+/urq6Z8+excfHb9iw4fk2ZWVl9+/fnz9/fqd6BweHqVOn4kcmwsLCxo8fr1AoOq1nCV4UnU738fF5VX/fIY9Go73xxhuO/3J2du7NineDnbu7e1BQ0Ntvvz1s2DCEkKurK4fDiYqKMrdH6TGJROLp6bl06VJcpFKpPB5PqVR2lzJgijVgjn799dfLly9nZGSQHQgAYOiDa4QAAAAsGhwRAnP09OnT1tZWs5rOBgAwVEEiBAAAYNHg1CgAAACLBokQAACARYNECAAAwKJBIgQAAGDRIBECAP7DpUuXLly4QHYUAPQXuGsUAPAfIiMj1Wo1XnkHgKEHjggBGDhPnz5VKpV966tQKKRSqV6vf/kw5HL5P//8093WlpYWmUzWy6FaWlqePHnSw+9phULRw74AMAeQCAHoo+LiYhaLVVhYiItnz55lsVgbN27ERZ1ON378eLxwh1qtXrJkCZPJZLFYeKbKI0eO4GZyudze3v7rr782HVkkErFYrLy8PFy8evVqQEAAnU53dHR0cnJKTU3tMp6bN2+yWKxO2WS/4wAABrtJREFUK3GnpaWxWCyJRIKLv/32m7u7u52d3ZgxY9zd3XNzc00bnzx5ctKkSQwGw97enk6nf/rppwghPp8vEAgKCgpYLBaLxYqMjMSNr1+/HhAQwGAwHBwcnJ2dDx06ZBzn2LFjLBYrPz8/ODiYTqfPmDHjRT9bAAaUAQDQJxqNhkajff7557j4/vvvW1tbu7m54WJRURFC6OrVqwaDobm5OSYm5vz58+Xl5X/88UdcXBxBEHl5ebhlVFSUh4eHXq83jhwXF2dvb6/RaAwGg0gkGj58+Pz58/Pz8+/cubN161aEUGZmZpchTZw4MSoqyrQmODj4zTffxK/T09MRQhs2bCguLi4pKYmJiaFQKIWFhXjr0aNHEUJRUVF5eXnl5eWnTp3CS3wUFRUFBQVxOJwrV65cuXKlpKTEYDBUVFSMGjVqypQpuCY2NhYh9NNPP+GhDh8+jBAaN25cUlKSWCy+fPnyy3/aAPQfSIQA9F1kZCSPx8OvPTw81q9fjxC6f/++wWD48ssvqVSqWq1+vpder/f19V21ahUunj59GiFkTEhtbW22trYfffQRLnK53MDAQK1Wa+y+YMGCyZMndxlPcnIyhUJ5/PgxLlZVVREEgfOTRqMZO3ZsdHS0sbFOp+NwOIsWLTIYDM+ePbO3tw8ODjbNx0azZs0yvk1s9erVVCq1vr7eWMPlcp2cnHB3nAi/+OKLrj81AMwMrEcIQN/x+fxPPvlEqVQ2NTVVV1efOXPm4sWLAoHg9ddfFwgEISEhVCoVt1Sr1VlZWZWVlXK5HCGkVCrv37+PN0VFRdnZ2WVkZAQHByOEsrOzW1tbY2JiEELNzc3FxcXLli27du2acacODg65ublarZZC6fz9jYuL2759+4kTJzZv3owQOnr06IgRIxYvXowQun37tkwmc3NzM55xRQi5urr+9ddfCKGysrInT54kJSURBNGbN3779u2QkBDT1cmXLFmSkJBQU1NjXBXrnXfeeZHPEgDSQCIEoO/4fL5GoxGLxQ0NDfb29l5eXmFhYQKBICYm5vr16zt27MDNHjx4MH36dK1WO3PmTDs7O2tr6xEjRrS0tOCt1tbWS5cu/eWXX7777ruRI0dmZGQEBgZyOByEkEwmMxgMOTk5nS7m2djYNDU1jR07tlM8zs7OfD4/PT198+bNer0+MzNzwYIFDAYDIfTkyROE0MGDB3/88UfTLnQ6HSGE72dxcXHp5Ruvq6vz8vIyrcFJUS6XGxMhXhMVAPMHiRCAvps8efKYMWMEAkFDQ0N4eDhBEHw+/4MPPsjPz1er1Xw+Hzc7dOhQe3t7ZWUlm83GNcXFxc3NzcZx4uLiDhw4cOHCBS6XKxQKU1JScP3o0aMRQlu3bt22bVsvQ4qLi3vvvffu3LnT2NhYV1cXFxeH63HCS0tLW7Ro0fO9TJNlb9jY2HS6FxTfaIr3guEVXAEwf/AvFYC+IwgiNDQ0Ly9PKBTitMfn8+Vy+b59+xgMhp+fH24mkUjc3NyMWVAmk5WUlJiOExAQwOFwMjIyjh07RqFQoqOjcb2Dg4Onp2dOTo5Op+tlSPgQMCMjIyMjAx8g4np/f38ajXbq1Kkue3E4HAaD0d1WGo2mVqtNa4KDg8VisWkuz83NZbPZEyZM6GWcAJgRsi9SAjC4Gc801tbW4hp8znDevHnGNklJSRQKJSsrq729vaysLCQkhEaj+fn5mY7zzTffUCgUFxeXhQsXmtafOXOGIIh33323rKxMpVI9fPjwxIkTn332WQ8hrVmzZsyYMTY2Nlu3bjWt37NnD0Jo06ZNNTU1KpWquro6NTV1//79xgAQQhs3bqypqVEqlaWlpca7QHfs2EGlUo8fP15SUnLv3j2DwVBcXEyhUGbNmlVVVSWTyfBJ4OTkZNwe3yyDz+sCYP4gEQLwUqqrqxFCEyZMMNbgRwkPHjxorFEqlW+99RbOlxQKJSEhITo6ulMilEqlw4cPRwhduHCh0y6ysrJM1yhmMpk9J8KCggLcsqKiwrRer9fv27fPeGCKEHJycjp8+LBx6+7du21tbfGmYcOGxcfH401yuXzhwoW447Rp03Bldna28SqgtbV1YmKiTqfDmyARgsEFplgDYIDU1tbKZDJ3d3cmk/mifQ0GQ1VVVUtLi52d3Wuvvfb8/aK9p9Vq7927p1KpHB0dnZ2dO13J6+joqKio0Ol0rq6upimzSzqd7u7du2q1etKkScYMCsCgA4kQAACARYObZQAAAFg0SIQAAAAsGiRCAAAAFg0SIQAAAIsGiRAAAIBFg0QIAADAov0/bC0/FUXTT+IAAAAASUVORK5CYII=", "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" ], "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" ] }, "metadata": {}, "execution_count": 4 } ], "cell_type": "code", "source": [ "using Plots\n", "using Unitful\n", "\n", "p = plot_bandstructure(basis; n_bands=8, εF, kline_density=10, unit=u\"eV\")\n", "ylims!(p, (-5, 6))" ], "metadata": {}, "execution_count": 4 } ], "nbformat_minor": 3, "metadata": { "language_info": { "file_extension": ".jl", "mimetype": "application/julia", "name": "julia", "version": "1.6.3" }, "kernelspec": { "name": "julia-1.6", "display_name": "Julia 1.6.3", "language": "julia" } }, "nbformat": 4 }