{
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# Swift transit model example III - LDTk-based limb darkening\n",
"\n",
"*Author:* Hannu Parviainen
\n",
"*Last modified:* 7.7.2020\n",
"\n",
"The *LDTk limb darkening model*, ``pytransit.LDTkLDModel``, works as an example of a more complex limb darkening model that is best implemented as a subclass of ``pytransit.LDModel``. The *LDTk limb darkening model* uses *LDTk* to create a set of stellar limb darkening profile samples given the stellar $T_\\mathrm{Eff}$, $\\log g$, and $z$ with their uncertainties, and uses the profiles directly to calculate the transit. The profiles are created from the PHOENIX-calculated *specific intensity spectra* by Husser (2013), and the model completely avoids approximating the limb darkening profile with an analytical function.\n",
"\n",
"The model is parameter free after the stellar parameters have been given. The model can be *frozen* for model optimisation, and *thawn* for MCMC posterior estimation. When *frozen*, the model returns the average limb darkening profile interpolated from the profile at the given $\\mu$ locations. When *thawn*, each model evaluation chooses a random limb darkening profile from the sample and uses interpolation to evaluate the model at the wanted $\\mu$ values. "
]
},
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Populating the interactive namespace from numpy and matplotlib\n"
]
}
],
"source": [
"%pylab inline\n",
"rc('figure', figsize=(13,6))"
]
},
{
"cell_type": "code",
"execution_count": 74,
"metadata": {},
"outputs": [],
"source": [
"def plot_lc(time, flux, c=None, ylim=(0.9865, 1.0025), ax=None, alpha=1):\n",
" if ax is None:\n",
" fig, ax = subplots()\n",
" else:\n",
" fig, ax = None, ax\n",
" ax.plot(time, flux, c=c, alpha=alpha)\n",
" ax.autoscale(axis='x', tight=True)\n",
" setp(ax, xlabel='Time [d]', ylabel='Flux', xlim=time[[0,-1]], ylim=ylim)\n",
" \n",
" if fig is not None:\n",
" fig.tight_layout()\n",
" return ax"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Import the model\n",
"\n",
"First, we import the ``SwiftModel`` and ``LDTkLDModel`` and some simple transmission functions from *LDTk*."
]
},
{
"cell_type": "code",
"execution_count": 29,
"metadata": {},
"outputs": [],
"source": [
"from pytransit import SwiftModel, LDTkLDModel\n",
"from ldtk import sdss_g, sdss_i, sdss_z"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {},
"outputs": [],
"source": [
"time = linspace(-0.05, 0.05, 1500)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Example 1: single passband\n",
"\n",
"The *LDTkLDModel* is initialised by giving it the stellar parameters and passband transmission functions,"
]
},
{
"cell_type": "code",
"execution_count": 80,
"metadata": {},
"outputs": [],
"source": [
"ldm = LDTkLDModel(teff=(5500, 150), logg=(4.5, 0.1), z=(0.0, 0.1), pbs=[sdss_i], frozen=True,\n",
" cache='/Users/hannu/Dropbox/data/Husser2013_bin50/')"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"and given to the ``SwiftModel`` as any other limb darkening model."
]
},
{
"cell_type": "code",
"execution_count": 81,
"metadata": {},
"outputs": [],
"source": [
"tm = SwiftModel(ldm)\n",
"tm.set_data(time)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"after which the transit model evaluation goes as usual"
]
},
{
"cell_type": "code",
"execution_count": 82,
"metadata": {},
"outputs": [],
"source": [
"flux1 = tm.evaluate(k=0.1, ldc=[None], t0=0.0, p=1.0, a=4.2, i=0.5*pi, e=0.0, w=0.0)"
]
},
{
"cell_type": "code",
"execution_count": 83,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAA6AAAAGoCAYAAABYCYj1AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nOzdeXTV5b3v8c83M4QhIyGEQICEIYxiQGSmakWOFbW11Q5ah1o92tbeM9S267RnuL3Vtuf29nSynqOttlZrrVZbUcEJCPM8DwlzCIQMkBBChp08949sbMQAAbL3bw/v11qs7P0bdj4bswyf/Ty/32POOQEAAAAAEGgxXgcAAAAAAEQHCigAAAAAICgooAAAAACAoKCAAgAAAACCggIKAAAAAAiKOK8DBENGRobLy8vzOgYAAAAARIV169ZVOecyz94eFQU0Ly9Pa9eu9ToGAAAAAEQFMzvQ2Xam4AIAAAAAgoICCgAAAAAICgooAAAAACAoKKAAAAAAgKCggAIAAAAAgoICCgAAAAAICgooAAAAACAoAlZAzexpMztmZlvPsd/M7L/MrNTMNpvZxA775prZLv++Rzts/6GZ7fQf/4qZpQQqPwAAAACgewVyBPQ3kuaeZ/8Nkgr8f+6X9EtJMrNYST/37y+UdIeZFfrPWSRpjHNunKTdkr4ZkOQAAAAAgG4XsALqnFsiqeY8h8yX9Kxrt1JSipllS5osqdQ5t9c51yzpBf+xcs4tdM75/OevlDQwUPkBAAAAAN3Ly2tAcyQd6vC8zL/tXNvPdo+kN8714mZ2v5mtNbO1lZWV3RAXAAAAAHA5vCyg1sk2d57tfzvR7NuSfJKeO9eLO+eedM4VOeeKMjMzLysoAAAAAODyxXn4vcsk5XZ4PlBSuaSEc2yXJJnZXZJulHSNc+5DxRQAAAAAELq8HAF9TdKd/rvhTpFU65w7ImmNpAIzG2JmCZJu9x8rM5sr6RuSbnLONXgVHAAAAABw8QI2Ampmz0uaLSnDzMokfVdSvCQ5556QtEDSPEmlkhok3e3f5zOzhyW9JSlW0tPOuW3+l/2ZpERJi8xMklY65x4I1HsAAAAAAHSfgBVQ59wdF9jvJD10jn0L1F5Qz96e3z3pAAAAAADB5uUUXAAAAABAFKGAAgAAAACCggIKAAAAAAgKCigAAAAAICgooAAAAACAoKCAAgAAAACCggIKAAAAAAgKCigAAAAAICgooAAAAACAoKCAAgAAAACCggIKAAAAAAgKCigAAAAAICgooAAAAACAoKCAAgAAAACCggIKAAAAAAgKCigAAAAAICgooAAAAACAoKCAAgAAAACCggIKAAAAAAgKCigAAAAAICgooAAAAACAoKCAAgAAAACCggIKAAAAAAgKCigAAAAAICgooAAAAACAoKCAAgAAAACCggIKAAAAAAgKCigAAAAAICgooAAAAACAoKCAAgAAAACCggIKAAAAAAgKCigAAAAAICgooAAAAACAoKCAAgAAAACCggIKAAAAAAgKCigAAAAAICgooAAAAACAoKCAAgAAAACCImAF1MyeNrNjZrb1HPvNzP7LzErNbLOZTeywb66Z7fLve7TD9jQzW2RmJf6vqYHKDwAAAADoXoEcAf2NpLnn2X+DpAL/n/sl/VKSzCxW0s/9+wsl3WFmhf5zHpX0jnOuQNI7/ucAAAAAgDAQF6gXds4tMbO88xwyX9KzzjknaaWZpZhZtqQ8SaXOub2SZGYv+I/d7v8623/+M5Lel/SNC2VpbXM60dAsSTLZ33Z0fGgf2SQz6+zQDsfaR7apk+POd6x9sO0C36uzbwAAAADAc8451Tf51NrmvI4S8gJWQLsgR9KhDs/L/Ns6236V/3GWc+6IJDnnjphZv658o+1H6jTh3xddfuIQdClldmT/PvrmvJGaOiwj4PkAAACASNXY0qofv71bv1txQKeaW72OExa8LKCdDem582y/uBc3u1/tU3uVkTNE3/1EoVyHV+n4gs599OU/fKzrZFvXj/3w617ea33oZc+81gXO77i91Tm9vvmIPvc/q/To3JH68qxhnQcFAAAAcE6VJ5t019Ortf1InW4aP0BjcvooPpZ7vJ5xz+Odb/eygJZJyu3wfKCkckkJ59guSRVmlu0f/cyWdOxcL+6ce1LSk5JUVFTk7p42pDuzh7VHrhmuf3xpk77/xk45SQ9QQgEAAIAuqznVrM88uUJHTjTqqbuKdM2oLK8jhZx7zrHdy4r+mqQ7/XfDnSKp1j+9do2kAjMbYmYJkm73H3vmnLv8j++S9GqwQ0eCHgmx+untV+gT4wfosTd2asGWI15HAgAAAMKCr7VND/9+vcpqTus3d0+ifF6kgI2Amtnzar9hUIaZlUn6rqR4SXLOPSFpgaR5kkolNUi627/PZ2YPS3pLUqykp51z2/wv+5ikF83sXkkHJd0WqPyRLibG9J+3jdehmgZ946XNGjOgrwal9/Q6FgAAABDSHntjp5bvqdYPPzVOVw1N9zpO2LHOrn+MNEVFRW7t2rVexwhJh2oa9Hf/tVRDMpL1pwenKo556wAAAECnVu2t1meeXKkvTBms/7h5jNdxQpqZrXPOFZ29nbYR5XLTeup7t4zVprJa/W7lAa/jAAAAACGpsaVVj768RYPSeuqb80Z6HSdsUUChG8dla+bwTP3nwt06VtfodRwAAAAg5Pzs3VLtqzqlx24dq54JXt7LNbxRQCEz07/fNFpNrW36Pwt2eB0HAAAACCnlJ07rv5fu1c0TBmhqfobXccIaBRSSpLyMZN03fYj+vLFcWw/Xeh0HAAAACBk/WrhLTtI/Xj/C6yhhjwKKD3x51jCl9IzXD97a5XUUAAAAICRsPVyrVzYc1j3ThmhgKqtGXC4KKD7Qt0e8Hp6TryW7K7WstMrrOAAAAIDnvv/GDqX0iNffzxnmdZSIQAHFh3x+ymAN6Jukx9/cqWhYogcAAAA4l7X7a7SstFoPzclXn6R4r+NEBAooPiQpPlZfu7ZAm8tqtXh3pddxAAAAAM/89N1SpScn6LNXDfI6SsSggOIjbrlioLL7JumX7+/xOgoAAADgic1lJ7R4d6XunTGEZVe6EQUUH5EQF6P7ZgzVqn01WnfguNdxAAAAgKD76bul6tsjXl+YMtjrKBGFAopO3T4pVyk94/XEYkZBAQAAEF12V5zUou0V+uLUPPXm2s9uRQFFp5IT4/TFqXlatL1CJRUnvY4DAAAABM2vl+1XYlyMvjg1z+soEYcCinO66+o8JcbF6JkV+72OAgAAAATFiYZmvbKhTLdckaPU5ASv40QcCijOKTU5QZ8YP0Avrz+susYWr+MAAAAAAfeHNYfU2NKmuxj9DAgKKM7rzqsHq6G5VS+vK/M6CgAAABBQvtY2PbvigKYMTdOo7D5ex4lIFFCc17iBKZqQm6JnVx6Qc87rOAAAAEDAvL2jQodPnNbd04Z4HSViUUBxQXdePVh7K09pWWm111EAAACAgPn96kMa0DdJ147K8jpKxKKA4oLmjc1WWnKCfrfygNdRAAAAgIAoO96gpSWVuq0oV7Ex5nWciEUBxQUlxcfq1ity9PaOClXXN3kdBwAAAOh2f1zbfs+T24oGepwkslFA0SW3FeXK1+b0yobDXkcBAAAAulVrm9Mf1x7SjIJMDUzt6XWciEYBRZeM6N9b43NT9OLaQ9yMCAAAABGluLRK5bWNun1SrtdRIh4FFF32maJc7a6o1+ayWq+jAAAAAN3mD2sOKi05gZsPBQEFFF124/hsJcXH6MW1h7yOAgAAAHSLqvomLdpeoVuvyFFCHPUo0PgbRpf1SYrXvDHZem1juU43t3odBwAAALhsr24sV0ur02eYfhsUFFBclNuKcnWyyac3tx3xOgoAAABw2f684bDG5vRVQVZvr6NEBQooLsqUoWkamNpDf95Q7nUUAAAA4LKUHqvXlsO1uvmKHK+jRA0KKC6KmWn+hAEqLq1SFWuCAgAAIIy9uvGwYkz6xPhsr6NEDQooLtr8CTlqbXN6fTPTcAEAABCenGtf435afob69U7yOk7UoIDiog3P6q2R/Xvr1Y2HvY4CAAAAXJJ1B46r7Php3cL026CigOKSzJ+Qo/UHT+hgdYPXUQAAAICL9sqGw+oRH6vrR/f3OkpUoYDikpyZJ/+XzdyMCAAAAOGl2dem17cc0XWFWUpOjPM6TlShgOKSDEztqUl5qfrzhsNyznkdBwAAAOiypSWVOtHQopuvGOB1lKhDAcUlu2lCjkqO1Wv7kTqvowAAAABd9vqWI+qTFKfp+ZleR4k6FFBcsr8bm63YGNOCLdwNFwAAAOGhydeqRdsr9PHR/ZUQRx0KNv7GccnSkhM0ZWia3thylGm4AAAACAvLSqt0stGneWO5+ZAXKKC4LHPHZGtv1Sntrqj3OgoAAABwQQu2HFVvpt96hgKKy3L96CyZSW9sZRouAAAAQluzr00Ltx3VdYVZTL/1CH/ruCz9eiepaHCq3tx61OsoAAAAwHkt21Olukaf5o3J9jpK1KKA4rLdMCZbO4+e1L6qU15HAQAAAM5pweYj6p0YpxnDM7yOErUCWkDNbK6Z7TKzUjN7tJP9qWb2ipltNrPVZjamw76vmdlWM9tmZo902D7BzFaa2UYzW2tmkwP5HnBhc8e0X8DNNFwAAACEqpbWNi3cXqFrC7OUGBfrdZyoFbACamaxkn4u6QZJhZLuMLPCsw77lqSNzrlxku6U9BP/uWMkfUnSZEnjJd1oZgX+c34g6d+ccxMkfcf/HB4akNJD43NTmIYLAACAkLVqb41qT7d8MHgCbwRyBHSypFLn3F7nXLOkFyTNP+uYQknvSJJzbqekPDPLkjRK0krnXINzzidpsaRb/Oc4SX38j/tKKg/ge0AX3TCmvzaX1epQTYPXUQAAAICPWLj9qJLiYzSzgLvfeimQBTRH0qEOz8v82zraJOlWSfJPpR0saaCkrZJmmlm6mfWUNE9Srv+cRyT90MwOSfqRpG929s3N7H7/FN21lZWV3fSWcC43+D9Jemsbo6AAAAAILc45LdxWoZkFmeqRwPRbLwWygFon29xZzx+TlGpmGyV9RdIGST7n3A5Jj0taJOlNtRdVn/+cByV93TmXK+nrkp7q7Js75550zhU554oyM/mUI9AGpyerMLuP3mAaLgAAAELMlsO1OlrXqI+PZvqt1wJZQMv0t1FLqX1k80PTZZ1zdc65u/3Xc94pKVPSPv++p5xzE51zMyXVSCrxn3aXpJf9j/+o9qm+CAHXj+6v9QePq6q+yesoAAAAwAcWbqtQjEnXjOzndZSoF8gCukZSgZkNMbMESbdLeq3jAWaW4t8nSfdJWuKcq/Pv6+f/Okjt03Sf9x9XLmmW//HH9LdiCo9dW9hPzknv7jjmdRQAAADgAwu3H9XkIWlKTU648MEIqLhAvbBzzmdmD0t6S1KspKedc9vM7AH//ifUfrOhZ82sVdJ2Sfd2eIk/mVm6pBZJDznnjvu3f0nST8wsTlKjpPsD9R5wcQqz+ygnpYcW7ajQpyflXvgEAAAAIMD2V53S7op6fefGsxfkgBcCVkAlyTm3QNKCs7Y90eHxCkkFZ5/n3zfjHNuLJV3ZjTHRTcxM147qpz+sPaTGllYlxXOBNwAAALy1aHuFJOm6wiyPk0AK7BRcRKFrC7PU2NKm4pIqr6MAAAAAWrj9qEZl91FuWk+vo0AUUHSzq4akq3dinN7eUeF1FAAAAES5qvomrT1wXB9n9DNkUEDRrRLiYjRrRKbe3nFMbW1nr7oDAAAABM87OyrknPTx0RTQUEEBRbe7rjBLVfVN2lR2wusoAAAAiGKLth9TTkoPFWb38ToK/Cig6Hazh/dTbIx9cME3AAAAEGxNvlYtK63SnJGZMjOv48CPAopu17dnvCbnpXEdKAAAADyzel+NTre06mMj+3kdBR1QQBEQ1xVmaXdFvQ5Un/I6CgAAAKLQuzuPKTEuRlcPzfA6CjqggCIgrh3VfqH3uzuPeZwEAAAA0ej9XZW6eli6eiSwNn0ooYAiIAal99TQzGS9t6vS6ygAAACIMvuqTmlf1SnNGcH021BDAUXAzBnRTyv3Vut0c6vXUQAAABBF3vPPwqOAhh4KKAJmzoh+ava1acXeKq+jAAAAIIq8t+uYhmUma1B6T6+j4CwUUATMpCGp6pkQq/d2Mg0XAAAAwdHQ7NOqvTWMfoYoCigCJjEuVlOHZei9XcfknPM6DgAAAKLAstJqNbe2aQ7Lr4QkCigCas7ITJUdP609lfVeRwEAAEAUeG/XMSUnxGpSXprXUdAJCigCarZ/6sP73A0XAAAAAeac03s7j2l6QYYS4qg6oYj/KgionJQeGp7VS+/tYj1QAAAABNauipM6UtuojzH9NmRRQBFwc0b00+p9Napv8nkdBQAAABHszM0vZ3MDopBFAUXAzR7RTy2tTstKWY4FAAAAgfPezmMqzO6jrD5JXkfBOVBAEXBFeanqlRin95mGCwAAgACpbWjRuoPHmX4b4iigCLj42BhNz8/Q+7sqWY4FAAAAAbG0tFKtbU5zRmZ6HQXnQQFFUMwZmakjtY3aefSk11EAAAAQgd7deUwpPeM1ITfV6yg4DwoogmLW8PapEIt3sxwLAAAAupdzTktLqjSjIFOxMeZ1HJwHBRRB0b9vkkZk9dbSEgooAAAAuteuipOqPNmkGQUZXkfBBVBAETQzh2dozb7jamhmORYAAAB0n+KS9tUWpudTQEMdBRRBM6MgU82tbVq1t8brKAAAAIggxaVVGpqZrAEpPbyOgguggCJoJg9JU2JcjJYwDRcAAADdpMnXqlV7azSD0c+wQAFF0CTFx+qqoelawo2IAAAA0E02HDyh0y2tmkYBDQsUUATVzIIM7ak8pcMnTnsdBQAAABGguKRKsTGmKcPSvY6CLqCAIqhmDm9fGHgpo6AAAADoBktLqzQhN0V9kuK9joIuoIAiqAr69VL/PklcBwoAAIDLVtvQoi1lJ5h+G0YooAgqM9OMggwVl1TJ19rmdRwAAACEsRV7q9TmxPqfYYQCiqCbOTxTdY0+bSqr9ToKAAAAwtjSkiolJ8RqQm6K11HQRRRQBN30/AyZSUuZhgsAAIDLUFxapSlD0xUfS60JF/yXQtClJido3MAUlmMBAADAJTtU06AD1Q2azvTbsEIBhSdmFmRo46ETqm1o8ToKAAAAwlBxaZUkrv8MNxRQeGLm8Ey1OWn5niqvowAAACAMFZdUKatPooZl9vI6Ci4CBRSemJCbot6JcSzHAgAAgIvW1ua0bE+Vpudnysy8joOLQAGFJ+JjY3T1sHQt2V0l55zXcQAAABBGtpXX6URDC9Nvw1BAC6iZzTWzXWZWamaPdrI/1cxeMbPNZrbazMZ02Pc1M9tqZtvM7JGzzvuK/3W3mdkPAvkeEDgzhmfq8InT2l/d4HUUAAAAhJGlpe2z6Kbmp3ucBBcrYAXUzGIl/VzSDZIKJd1hZoVnHfYtSRudc+Mk3SnpJ/5zx0j6kqTJksZLutHMCvz75kiaL2mcc260pB8F6j0gsKbnt39ideYCcgAAAKArlpVWaWT/3urXO8nrKLhIgRwBnSyp1Dm31znXLOkFtRfHjgolvSNJzrmdkvLMLEvSKEkrnXMNzjmfpMWSbvGf86Ckx5xzTf7zjgXwPSCA8tJ7Kielh5aVUEABAADQNY0trVqz//gHgxkIL4EsoDmSDnV4Xubf1tEmSbdKkplNljRY0kBJWyXNNLN0M+spaZ6kXP85wyXNMLNVZrbYzCYF8D0ggMxM0/LTtXxPlVrbuA4UAAAAF7Z6X42afW2s/xmmAllAO7sd1dkt4zFJqWa2UdJXJG2Q5HPO7ZD0uKRFkt5Ue1H1+c+Jk5QqaYqkf5L0onVy6yszu9/M1prZ2spK7rQaqqYXZKqu0aeth2u9jgIAAIAwsKy0SgmxMZo8JM3rKLgEgSygZfrbqKXUPrJZ3vEA51ydc+5u59wEtV8Dmilpn3/fU865ic65mZJqJJV0eN2XXbvVktokfeTjD+fck865IudcUWZmZne/N3STqcPaLxznOlAAAAB0xdKSKk0cnKKeCXFeR8ElCGQBXSOpwMyGmFmCpNslvdbxADNL8e+TpPskLXHO1fn39fN/HaT2abrP+4/7s6SP+fcNl5QgifYSpjJ6JWpUdh8to4ACAADgAqrqm7T9SJ1mFDDAFK4C9rGBc85nZg9LektSrKSnnXPbzOwB//4n1H6zoWfNrFXSdkn3dniJP5lZuqQWSQ855477tz8t6Wkz2yqpWdJdjoUkw9r0/HQ9s/yATje3qkdCrNdxAAAAEKKW76mWJE3jBkRhK6Dj1s65BZIWnLXtiQ6PV0gqOMe5M86xvVnS57sxJjw2LT9D/710n9YeqOHTLAAAAJxTcUml+vaI19icvl5HwSUK5BRcoEsmD0lTfKxxHSgAAADOyTmn4pIqTR2WrtiYzu53inBAAYXneibEaeKgVK4DBQAAwDntqzql8tpGpt+GOQooQsL0/AxtK69Tzalmr6MAAAAgBJ2ZLTeD9T/DGgUUIWFaQYack5bvYRQUAAAAH7W0pEq5aT00OD3Z6yi4DBRQhIRxOX3VOzGOabgAAAD4CF9rm1buqdZ0pt+GPQooQkJcbIymDEvnRkQAAAD4iE1ltTrZ5NP0fFZMCHcUUISM6fkZOlRzWgerG7yOAgAAgBBSXFIlM2nqsHSvo+AyUUARMs7c0YxRUAAAAHS0rLRKYwb0VWpygtdRcJkooAgZwzKT1b9PEteBAgAA4AP1TT6tP3hc07n7bUSggCJkmJmm5Wdo2Z4qtbU5r+MAAAAgBKzaWy1fm9MMbkAUESigCCkzCjJ0oqFF24/UeR0FAAAAIaC4tEpJ8TGaODjV6yjoBhRQhJSp+e0XlnMdKAAAAKT2GxBNyktTUnys11HQDSigCCn9eidpRFZvrgMFAACAjtY2quRYvWZw/WfEoIAi5EzLz9DqfTVqbGn1OgoAAAA8dGZQgvU/IwcFFCFnekG6mnxtWn/guNdRAAAA4KHi0iqlJydoZP/eXkdBN6GAIuRMHpKuuBjjOlAAAIAo5pxTcWmVpuVnKCbGvI6DbkIBRcjplRinKwalcB0oAABAFNtdUa/Kk02s/xlhKKAISdPyM7T5cK1qG1q8jgIAAAAPLC2plCRNZ/3PiEIBRUianp8h56QVexkFBQAAiEbFpVUampmsASk9vI6CbkQBRUgan5ui5IRYLS2hgAIAAESbZl+bVu2t0QxGPyMOBRQhKT42RlOGpnMdKAAAQBRaf/C4Tre0ahoFNOJQQBGypuVnaH91gw7VNHgdBQAAAEFUXFKl2BjTlGHpXkdBN+tSATWzwk62ze72NEAHM/x3PFu+h1FQAACAaFJcWqUJuSnqkxTvdRR0s66OgL5oZt+wdj3M7KeSvh/IYEB+v17q1zuR60ABAACiSG1DizaXnWD6bYTqagG9SlKupOWS1kgqlzQtUKEASTIzTc/P0PI91Wprc17HAQAAQBCs2FulNve32XCILF0toC2STkvqISlJ0j7nXFvAUgF+0wsyVHOqWTuO1nkdBQAAAEFQXFqlXolxmpCb4nUUBEBXC+gatRfQSZKmS7rDzF4KWCrA78zUi2Km4QIAAESF4pIqTRmapvhY7pcaibr6X/Ve59x3nHMtzrmjzrn5kl4NZDBAkrL6JGl4Vi8VsxwLAABAxDtU06D91Q1c/xnB4rp43DEzG3TWtsXdHQbozLT8DP1+1UE1trQqKT7W6zgAAAAIkDODDlz/Gbm6OgL6uqS/+r++I2mvpDcCFQroaHp+hpp8bVp/4LjXUQAAABBAxaVV6t8nScMye3kdBQHSpQLqnBvrnBvn/1ogabKk4sBGA9pdNTRdcTHGNFwAAIAI1tbmtLy0StPyM2RmXsdBgFzSlb3OufVqvyEREHC9EuN0xaAUCigAAEAE21Zep+MNLUy/jXBdugbUzP5Xh6cxkiZKqgxIIqAT0/Iz9JN3SnSioVkpPRO8jgMAAIBudmawgRsQRbaujoD27vAnUe3Xgs4PVCjgbDMKMuSctHxPtddRAAAAEADFpZUa2b+3Mnsneh0FAdSlEVDn3L8FOghwPuMGpqhXYpyKS6s0b2y213EAAADQjRpbWrVm/3HdOWWw11EQYOctoGb2F0nuXPudczd1eyKgE/GxMZoyNF3LuA4UAAAg4qzZX6NmX5umc/1nxLvQCOiPgpIC6ILp+el6e0eFDtU0KDetp9dxAAAA0E2KS6qUEBujyUPSvI6CALtQAd3nnDsYlCTABZz5RKy4tEp3TB7kcRoAAAB0l6UlVZo4OEU9E7p0hSDC2IVuQvTnMw/M7E8BzgKc17DMXurfJ4nlWAAAACJIdX2Tth+p03TufhsVLlRAO64AO/RiX9zM5prZLjMrNbNHO9mfamavmNlmM1ttZmM67PuamW01s21m9kgn5/6jmTkz4yc1SpiZpuVnaHlpldraznlpMgAAAMLImcGFGQWZHidBMFyogLpzPL4gM4uV9HNJN0gqlHSHmRWeddi3JG10zo2TdKekn/jPHSPpS5ImSxov6UYzK+jw2rmSrpPE9OAoM6MgQ8cbWrT9SJ3XUQAAANANlpZUKaVnvMbk9PU6CoLgQgV0vJnVmdlJSeP8j+vM7KSZXagBTJZU6pzb65xrlvSCPrp2aKGkdyTJObdTUp6ZZUkaJWmlc67BOeeTtFjSLR3O+7Gkf9ZFlmKEv6n56ZLa/0cFAACA8OacU3FJlablZyg2xi58AsLeeQuocy7WOdfHOdfbORfnf3zmeZ8LvHaOpEMdnpf5t3W0SdKtkmRmkyUNljRQ0lZJM80s3cx6SponKdd/3E2SDjvnNp3vm5vZ/Wa21szWVlZWXiAqwkW/3kkakdWb5VgAAAAiQOmxeh2ta9QMrv+MGhcaAb0cnX2EcfaI5WOSUs1so6SvSNogyeec2yHpcUmLJL2p9qLq85fRb0v6zoW+uXPuSedckXOuKDOT+eSRZHpBhlbvr1FjS6vXUQAAAHAZlvhntbH+Z/QIZAEtk3/U0m+gpPKOBzjn6pxzd2xT7EQAACAASURBVDvnJqj9GtBMSfv8+55yzk10zs2UVCOpRNIwSUMkbTKz/f7XXG9m/QP4PhBipudnqNnXprX7j3sdBQAAAJdhaUmlhmYma2Aqa7xHi0AW0DWSCsxsiJklSLpd0msdDzCzFP8+SbpP0hLnXJ1/Xz//10Fqn6b7vHNui3Oun3MuzzmXp/aSO9E5dzSA7wMhZvKQNMXHGsuxAAAAhLEmX6tW7a3RTO5+G1UCttKrc85nZg9LektSrKSnnXPbzOwB//4n1H6zoWfNrFXSdkn3dniJP5lZuqQWSQ855xjugiQpOTFOVwxKVXFppaSRXscBAADAJVh34LhOt7Sy/meUCVgBlSTn3AJJC87a9kSHxyskFZx9nn/fjC68ft5lRkSYmp6foR+/vVs1p5qVlpxw4RMAAAAQUpaWVCkuxjRlWLrXURBEgZyCCwTM9IIMOSet2FPtdRQAAABcgqUllZo4OFW9EgM6JoYQQwFFWBqX01e9k+L803ABAAAQTqrrm7StvI7lV6IQBRRhKS42RlcPTedGRAAAAGFo2Z5qOSfNGM4NiKINBRRha3pBhg7VnNb+qlNeRwEAAMBFWLq7Un17xGtsTl+voyDIKKAIW7P8n5gtKWEaLgAAQLhwzmlpSZWm5acrNsa8joMgo4AibA1OT1Zeek8t3kUBBQAACBd7Kut1tK5RM1j/MypRQBHWZg3P1PI91WrytXodBQAAAF2wZHf7PTxY/zM6UUAR1maNyNTpllat3X/c6ygAAADogqUllRqakazctJ5eR4EHKKAIa1OGpishNkaLdzMNFwAAINQ1+Vq1cm+Nphcw+hmtKKAIaz0T4jR5SJre33XM6ygAAAC4gHUHjut0SyvTb6MYBRRhb9bwTO2uqFf5idNeRwEAAMB5LN5VqfhY01QKaNSigCLszRrhX46FabgAAAAh7f1dlZqUl6ZeiXFeR4FHKKAIewX9eim7bxLXgQIAAISw8hOntavipGaPYPmVaEYBRdgzM80anqnikiq1tLZ5HQcAAACdODNYMHtEP4+TwEsUUESEWcMzdbLJp42HTngdBQAAAJ14f9cxDeibpIJ+vbyOAg9RQBERpuZnKDbGtHgX03ABAABCTbOvTcUlVZo1op/MzOs48BAFFBGhb494TRyUwnWgAAAAIWjtgRqdam7l+k9QQBE5Zg3P1JbDtaqqb/I6CgAAADo4s/zKNJZfiXoUUESMWcPbL2hfWsIoKAAAQChh+RWcQQFFxBg9oI/SkxP0PteBAgAAhAyWX0FHFFBEjJgY06wRmXp/V6V8LMcCAAAQElh+BR1RQBFRrh2VpdrTLVp74LjXUQAAACDpvZ0sv4K/oYAioswcnqmE2Bi9s6PC6ygAAABRr7GlVUtLqvSxUSy/gnYUUESUXolxumpomt7ZcczrKAAAAFGvuKRKp1tadf3o/l5HQYiggCLiXDsqS3urTmlPZb3XUQAAAKLawu1H1TsxTlcNSfc6CkIEBRQR55pR7Re4Mw0XAADAO61tTm/vOKY5I/spIY7agXb8JCDiDEztqZH9e+vt7UzDBQAA8Mq6A8dVc6pZHx+d5XUUhBAKKCLSdYVZWnugRsdPNXsdBQAAICot2n5UCbExmjWc9T/xNxRQRKRrRmWpzUnv72YUFAAAINicc1q4vUJXD0tX76R4r+MghFBAEZHG5fRVv96Jemsr14ECAAAE2+6Keh2obmD6LT6CAoqIFBNjun50f72/+5gamn1exwEAAIgqC7cdlSRdN4oCig+jgCJi3TC2vxpb2rR4V6XXUQAAAKLKwu0VumJQivr1SfI6CkIMBRQRa3JemtKTE7Rg61GvowAAAESN8hOnteVwra4rZPQTH0UBRcSKi43Rx0f317s7KtTY0up1HAAAgKjwtn8t9o8X9vc4CUIRBRQR7YYx/XWquVVLdjMNFwAAIBje2nZUQzOTld+vl9dREIIooIhoVw9LV98e8XqDabgAAAABV3OqWSv31mjuaEY/0TkKKCJafGyMPl6Ypbe3V6jJxzRcAACAQHpr21G1tjnNG5vtdRSEKAooIt68sdk62eTTstIqr6MAAABEtAVbjigvvadGD+jjdRSEqIAWUDOba2a7zKzUzB7tZH+qmb1iZpvNbLWZjemw72tmttXMtpnZIx22/9DMdvrPecXMUgL5HhD+puanq3dSnN7YwjRcAACAQKk51azle6o1b2y2zMzrOAhRASugZhYr6eeSbpBUKOkOMys867BvSdronBsn6U5JP/GfO0bSlyRNljRe0o1mVuA/Z5GkMf5zdkv6ZqDeAyJDYlysrhuVpYXbK9TS2uZ1HAAAgIi0kOm36IJAjoBOllTqnNvrnGuW9IKk+WcdUyjpHUlyzu2UlGdmWZJGSVrpnGtwzvkkLZZ0i/+4hf5tkrRS0sAAvgdEiBvGZqv2dItW7Kn2OgoAAEBEep3pt+iCQBbQHEmHOjwv82/raJOkWyXJzCZLGqz2QrlV0kwzSzeznpLmScrt5HvcI+mNzr65md1vZmvNbG1lJUtwRLsZBRlKTojV65uPeB0FAAAg4lTXNzH9Fl0SyALa2U+eO+v5Y5JSzWyjpK9I2iDJ55zbIelxtU+3fVPtRdXX8UQz+7Z/23OdfXPn3JPOuSLnXFFmZuZlvRGEv6T4WM0dk63XtxxRQ7PvwicAAACgy/68sVytbU43X3H2eBPwYYEsoGX68KjlQEnlHQ9wztU55+52zk1Q+zWgmZL2+fc95Zyb6JybKalGUsmZ88zsLkk3Svqcc+7sUgt06tNFA1Xf5NObrAkKAADQbZxz+uPaQxo/sK+GZ/X2Og5CXCAL6BpJBWY2xMwSJN0u6bWOB5hZin+fJN0naYlzrs6/r5//6yC1T9N93v98rqRvSLrJOdcQwPyIMJOHpGlwek/9cW2Z11EAAAAixrbyOu08elKfKursijngwwJWQP03CnpY0luSdkh60Tm3zcweMLMH/IeNkrTNzHaq/W65X+vwEn8ys+2S/iLpIefccf/2n0nqLWmRmW00sycC9R4QWcxMn5o4UCv2VutgNZ9dAAAAdIcX1x5SYlyMbho/wOsoCANxgXxx59wCSQvO2vZEh8crJBWcfZ5/34xzbM/vzoyILp+8cqB+/PZuPbf6gL55wyiv4wAAAIS1+iafXll/WHPH9FffHvFex0EYCOQUXCDkDEjpoRvGZOv5VQd1qombEQEAAFyOP649pJNNPt09bYjXURAmKKCIOvfOGKK6Rp9eWHPowgcDAACgU82+Nj29bJ+KBqdqQm6K13EQJiigiDoTB6Vq6rB0/ezdEtU2tHgdBwAAICz9duUBHao5rYc+xhVy6DoKKKLSt/9ulE6cbtFjb+4QK/kAAABcnLLjDfrJ27s1oyBDs4dneh0HYYQCiqg0ekBf3T9zqJ5ffUg/f6/U6zgAAABh4/ipZt319Go5Sf9602iZmdeREEYCehdcIJR94/qROlbXpB8t3K2Ungn6/JTBXkcCAAAIaY0trbrv2bU6VHNaz947WcMye3kdCWGGAoqoFRNj+sGnxqn2dIv+5dWtyuiVoLljsr2OBQAAEJJ8rW16+PcbtP7gcf38sxM1ZWi615EQhpiCi6gWHxujn392oibkpuirL2zU6n01XkcCAAAIOc45/cur2/T2jgr96ydGa95YPrTHpaGAIur1SIjVU3dN0sCUHrrvmTUqqTjpdSQAAICQ8l/vlOr51Qf14OxhumtqntdxEMYooICktOQEPXPPZCXGx+qup1frSO1pryMBAACEhOdXH9SP396tWyfm6J+vH+F1HIQ5Cijgl5vWU7+5e5LqGn364tNrVHuaNUIBAEB0W7S9Qt9+ZYtmDc/U458cxx1vcdkooEAHowf01a++cKX2VtXr/mfXqrGl1etIAAAAnlh34Li+8vx6jcnpq198bqLiY6kOuHz8FAFnmZafoR/dNl6r9tXoH17cpLY253UkAACAoNpWXqt7n1mj/n2S9PQXJyk5kcUz0D34SQI6MX9Cjo7VNel7C3YoNsb0o9vGKyGOz2sAAEDkW7O/Rvf8eo16J8XpmXsmK6NXoteREEEooMA5fGnmUPnanB5/c6eONzTric9fyad/AAAgor2385gefG6dBvTtod/ed5VyUnp4HQkRhiEd4DwenD1MP/jUOC3fU63P/vdKVdc3eR0JAAAgIP6yqVxfenathmX20osPXE35REBQQIEL+HRRrn71+Su18+hJ3fbECh2qafA6EgAAQLd6btUBffWFDZo4OFXP3z+FabcIGAoo0AXXFmbpufuuUlV9kz75y+XacaTO60gAAADd4hfvl+rbr2zVnBH99Ow9k9UnKd7rSIhgFFCgi4ry0vTHB6bKTPr0r1Zo9b4aryMBAABcsrY2p+8v2KEfvLlL8ycM0K++cKWS4mO9joUIRwEFLsKI/r31pwenKrN3or7w1Cot3HbU60gAAAAXrb7JpwefW6dfLdmrL0wZrB9/egLrfCIo+CkDLtLA1J566YGpGpndRw/8bp2eLt4n51grFAAAhIcD1ad06y+WadH2Cv3LjYX69/mjFRNjXsdClKCAApcgLTlBv7/vKl0zKkv//tft+uoLG3Wqyed1LAAAgPMqLqnSTT9bpoq6Jj17z1W6d/oQmVE+ETwUUOASJSfG6Vefv1L/dP0Ivb65XLf8Ypn2VtZ7HQsAAOAjnHP6n6V7defTq5TVJ1GvPTxN0wsyvI6FKEQBBS5DTIzpoTn5euaeyao82aSbfrZMb3FdKAAACCGNLa36hz9u0v9+fYeuK8zSy38/TYPTk72OhShFAQW6wYyCTP31qzM0NDNZX/7tOj3+5k75Wtu8jgUAAKLcoZoGfeZXK/Ty+sN65NoC/fJzV6pXYpzXsRDFKKBAN8lJ6aEXv3y17pg8SL98f4/u+vVqVdc3eR0LAABEqb9sKte8nyzV3spTeuLzV+qRa4dzsyF4jgIKdKOk+Fh9/9ax+sGnxmnN/uO68afFWneA9UIBAEDwNDT79I2XNusrz29QflYvLfjaDM0d09/rWIAkCigQEJ8uytXLD05VXKzptidW6Edv7VILU3IBAECAbS+v0yd+WqwX1x3SQ3OG6cUvX63ctJ5exwI+QAEFAmRMTl8t+OoMfXLiQP3svVLd+ovlKj3GXXIBAED3c87pmeX7dfMvlulko0+/u/cq/dP1IxUfyz/3EVr4iQQCqHdSvH5423g98fmJKjveoBt/ulTPrtgv55zX0QAAQISoOdWs+3+7Tt99bZumDUvXG1+boWn5LLGC0MQtsIAgmDsmWxMHpeqf/7RZ33l1m97ecUw//NQ4ZfVJ8joaAAAIU845LdhyVN99bavqTvv0nRsLdfe0PJlxoyGELkZAgSDp1ydJv/7iJP3H/NFava9a1/+/JVqw5YjXsQAAQBg6VteoB363Tg/9fr2y+/bQqw9P0z3Th1A+EfIYAQWCyMz0havzNDU/Q1//w0b9/XPrde2oLP3b/NHKSenhdTwAABDinHN6aV2Z/uOv29Xka9OjN4zUfdOHKI5rPREmLBquRSsqKnJr1671OgbwIS2tbfr1sn368aISmUlfv3a47p6Wxy8QAADQqbLjDfrWK1u1ZHelJuWl6vFPjtPQzF5exwI6ZWbrnHNFH9lOAQW8VXa8Qf/6Wvt1oaOy++h7t4zRxEGpXscCAAAhoqW1Tc8s368fL9otSfrGDSP1+asGKyaG6bYIXRRQCihCmHNOb22r0L++tk0VJxv12cmD9M9zR6pvj3ivowEAAA8tL63Sd1/bppJj9Zo9IlP/++YxGpjKup4IfecqoFwDCoQAM9PcMf01vSBD/7lwl55Zvl9vbavQv9w4SjeNH8ANBQAAiDLlJ07re6/v0Otbjig3rYf+584iXTOqH/8mQNhjBBQIQVsP1+pbr2zR5rJaTRyUon+5sVBXMC0XAICI1+Rr1f8s3aefvVuqNuf00Jx83T9zqJLiY72OBlwUpuBSQBFmWtucXlp3SD9auFuVJ5t00/gB+ue5I5h2AwBABHLOaeH2Cn1/wQ7tr27Q3NH99e2/G6XcNH7vIzydq4AG9HabZjbXzHaZWamZPdrJ/lQze8XMNpvZajMb02Hf18xsq5ltM7NHOmxPM7NFZlbi/8qwECJSbIzpM5MG6f1/nK2vfixfC7cf1cf+c7F+8OZOnWxs8ToeAADoJusPHtenf7VCX/7tOsXGmJ69Z7Ke+MKVlE9EpICNgJpZrKTdkq6TVCZpjaQ7nHPbOxzzQ0n1zrl/M7ORkn7unLvGX0RfkDRZUrOkNyU96JwrMbMfSKpxzj3mL7WpzrlvnC8LI6CIBOUnTutHb+3SyxsOK6NXgr5+3XB9piiXZVsAAAhT+6pO6Ydv7dSCLUeV0StRX7+ugN/tiBhejIBOllTqnNvrnGtWe6Gcf9YxhZLekSTn3E5JeWaWJWmUpJXOuQbnnE/SYkm3+M+ZL+kZ/+NnJN0cwPcAhIwBKT30fz8zQa89PE1DM3rp269s1bz/Wqq3th1VNEylBwAgUlTXN+lfX9um6/7vYr2/q1KPXFugxf80W5+7ajDlExEvkHfBzZF0qMPzMklXnXXMJkm3Sio2s8mSBksaKGmrpO+ZWbqk05LmSTozhJnlnDsiSc65I2bWr7Nvbmb3S7pfkgYNGtQtbwgIBeMGpugPX56it7Yd1Q/e3KUv/3adxg3sq/913XDNGp7J3fEAAAhR9U0+/WbZPj2xeK9Ot7TqM5Ny9cg1BerXJ8nraEDQBLKAdvav4LOHaR6T9BMz2yhpi6QNknzOuR1m9rikRZLq1V5UfRfzzZ1zT0p6UmqfgnuR2YGQ1r5sS7auHZWlVzYc1k/eKdEXf71Gk/JS9Q8fH6EpQ9O9jggAAPxONfn0zIr9+u8le3W8oUXXjsrSozeMUH6/3l5HA4IukAW0TFJuh+cDJZV3PMA5Vyfpbkmy9mGbff4/cs49Jekp/77/4389Saows2z/6Ge2pGMBfA9ASIuLjdFtRbmaPyFHf1hzUD99t1S3P7lS0/Mz9A8fH87SLQAAeKih2adnVxzQk0v2quZUs2aPyNQj1w7XhNwUr6MBnglkAV0jqcDMhkg6LOl2SZ/teICZpUhq8F8jep+kJf5SKjPr55w7ZmaD1D5N92r/aa9Jukvto6d3SXo1gO8BCAsJcTH6wtV5uq0oV79dcUC/XLxHt/xiuWaPyNTfz87X5CFpXkcEACBqnG5u1e9WHtATi/eo+lSzZg7P1CPXFmgiHwwDgV0H1MzmSfp/kmIlPe2c+56ZPSBJzrknzOxqSc9KapW0XdK9zrnj/nOXSkqX1CLpfznn3vFvT5f0oqRBkg5Kus05V3O+HNwFF9HmzDUmTy/br5pTzZqUl6q/n52v2SO4RhQAgECpbWjRb1fu16+X7Vf1qWbNKMjQI9cO15WDKZ6IPue6C25AC2iooIAiWjU0+/SHNYf030v2qry2UaOy++jB2cP0d2OzFRtDEQUAoDscrW3UU8V79ftVB3WquVVzRmTqoTn5KspjBhKiFwWUAooo1uxr06sbD+uXi/dob+Up5aX31JdnDdOtE3OUGBfrdTwAAMLSnsp6Pbl4r17eUKY2J31iXLa+PGuYRmX38Toa4DkKKAUUUGub08JtR/WL9/doy+FaZfZO1BemDNZnrxqkjF6JXscDACDkOee0al+Nfr1snxZur1BCbIw+MylXX5oxVLlpPb2OB4QMCigFFPiAc07FpVV6qnif3t9VqYS4GN08YYDumT5EI/vzqS0AAGdrbGnVXzaV6+ll+7XjSJ1Se8brc1cN1hen5fEhLtAJCigFFOhU6bGT+vWy/frT+jI1trRpWn667pk2RHNG9FMM14kCAKLcsbpG/W7lAT236qCqTzVreFYv3TNtiG6+IkdJ8VzGApwLBZQCCpzXiYZmPb/6kJ5Zvl9H6xo1JCNZn7tqkD45caBSkxO8jgcAQFBtOnRCv1m+X3/dXC5fm9M1I/vp7mlDNHVYOneUB7qAAkoBBbqkpbVNb2w9qt8s26f1B08oIS5GN47N1mevGqQrB6fySxcAELEamn16bWO5nlt1UFsO1yo5IVa3FeXqi1PzlJeR7HU8IKxQQCmgwEXbcaROv191UK9sOKz6Jp9GZPXW56YM0s1X5KhPUrzX8QAA6BY7j/p/360/rJP+33ef9/++683vO+CSUEApoMAlO9Xk0182let3qw5o6+E69YiP1Y3jsvWpKwdqUl4a14oCAMJO7ekWvbapXC+tPaRNZbUfzPj53JRBmjiIGT/A5aKAUkCBbrG57ISeW3lQf91crlPNrcpN66FPThyoT04cyO3nAQAhra3NadmeKv1xbZne3HZUzb42jezfW58uytUtV+RwzwOgG1FAKaBAt2po9umtbUf10royLd9TLeekKUPT9MmJAzVvbLaSE+O8jggAgCRpT2W9/rzhsF5ef1iHT5xW3x7xmj9hgD5dlKvRA/ow2gkEAAWUAgoETNnxBr2y/rBeWl+mA9UN6pkQqxvGZOuTV+ZoypB0pugCAIKuoq5Rf9lUrj9vPKyth+tkJk3Pz9Cni3J1XWEWS6gAAUYBpYACAeec07oDx/XSujL9dfMR1Tf5NKBvkj4xfoBumjBAhdl8ygwACJy6xha9ueWoXt10+IPZOWNz+mr+hAG6afwA9euT5HVEIGpQQCmgQFCdbm7Vwu1H9erGci3ZXSlfm9OwzGTNn5Cjm8YP4Hb2AIBu0eRr1Xs7K/XqxsN6Z+cxNfvaNDi9p+ZPyNH8CQM0LLOX1xGBqEQBpYACnqk51awFW47otU3lWr2vRpI0fmBffWL8AN0wNls5KT08TggACCetbU6r9lXr1Q3lWrD1iE42+pTRK0E3jhugm6/I0fiBfZlxA3iMAkoBxf9v785j7LoO+45/zyxv9n2lhhySGlEiKdpabMuS7DiKLNvyErto3cJp07hOWsNtAjhpgzRp/mlRFHWQIInbpjEMx3BSpwiyFTEM25IV2/EqW7YoiZREUlwkkkPOcBbOPvPmLad/3DujIcVNnHlvOJzvB3i4991z73vnEYd33u+dc8/VDeHMxDxffu4Mf/fMGZ4/MwXAG7e28Oi+Xt67bws77RmVJF3CQq7AD46N8fgLQ3z9hXOMzmRpyFTynn29fOjuPt420EFVZcV6V1NSygBqAJVuOCdGZ/nawSG+dvAsz56eBGB3bxPv2tvDI3t6eENfixMYSdImNjmf41uHz/H488N86/A5ZhcLNNZU8dAdXTy6r5d37u6hLuNkQtKNyABqAJVuaIMT83zt4BCPHRzix6+MU4zQ3VTDO/f08K693Tw40OmMhZK0CZydnOeJF4Z5/IVhfnBsjHwx0tVUw7v29vDuvT08MNBBTZV/D6QbnQHUACptGOOzi3zz0DmeeHGYbx8ZYXaxQF11Je+4vZNH9vTw8O5uOhpr1ruakqQ1UCxGnhuc5JuHzvHNw+d4Lh0Rc2tnA+++s5d339nD3VtbHREjbTAGUAOotCFl8wWePD7OEy8M88SLw5ydXCAEeFN/G4/s7eGRPd0MdDU62YQkbSCT8zm+89II3zh0jn84PMLY7CIhwN3bWnlkTw/vubOH27qb1ruaklbBAGoAlTa8GCPPn5niiReTMHpwMJnEqK+1jp/a1cnbd3XytoFO2hoy61xTSdJKMUaODM/wjbSX8yevnKdQjLTWV/PTt3fxM3d0847bu2j3/C3dNAygBlDppnNmYp5vHDrHd14a4fvHxpheyBNCctPxn9rVydtv6+JN29vIVDkroiSV29hMlh8cH+N7R8f49pERBifmAdi7pZmf2d3Fw7u7uXtbG5UOrZVuSgZQA6h0U8sXijx7epLvvjTKd14aYf+pCQrFSF11Jfff2s7bd3Xxjl2d3NbtcF1JKoXphRw/OjHO946O8f1joxwamgagsaaKBwc6eHh3Nw/d0U1vS+0611RSORhADaDSpjK1kOPJY2N89+go33lplBOjs0Ays+79t3akj3Z2djYYSCXpOizkCvzklfN87+go3z82xoHBSQrFSE1VBW/e0caDA508ONDBG/pavD+ntAkZQA2g0qZ2anyO7x4d5QfHxnjy+BjnprMAdC0H0nbeurODgS4DqSRdyvxigf0nz/Ojl8d58vgYT78ywWKhSGVF4O5trTw40MEDAx3c29/mbbMkGUANoJKWxBg5MTrLk8fH+eGJJJAOTyWBtLOxJgmjt3bwwK3tzrAradOamFvkxy8ngfNHJ8Y5ODhJvhgJAfb0NvPgQAdvu62Tt+xsp7Gmar2rK+kGYwA1gEq6jBgjL4/N8cPjSRh98vg4Q1MLALTWV3P3tlbu2dbGPf2t3N3fSnNt9TrXWJLWVoyR0+fnefrkeZ56eZynTpzn8HByDWemsoI3bm3hLTvbuW9HO/dub6OlzvOgpCu7XAD15ypJm14IgZ2dDezsbOAj9/UTY+Tk+BxPHh9j/8kJ9p+c4B+OHCFGCAFu62rknv5W7ulv497+Nm7rbnQWR0kbyuR8judOT/DMyQmeOTXBs6cnGJ1ZBJJJg+7d3sbP3rWFt+xo565trQ6plbRm7AGVpGswtZDjuVOT7D95nv2nJth/8jzn53JA8mXtrm0ty72kd21rpbOxZp1rLEmJXKHI4aFp9p9aCpznOTYyu1x+W3cjd21NRnjcs62V3b1NThokadUcgmsAlbSGlobt7j95nv0nJ3j65HkODU1TKCbn1N7mWvb1tbCvr5k39LWwr6+FnmZvPSCptPKFIkdHZjg4OMXBwUkODE5ycHCSbL4IQGdjhru3taaPNt64rcXLCiSVhAHUACqpxOYW8xw4/eoXvgODkxwfnWXpNNvVVMO+W5rTYJo8bmmpdZIjSdclmy9wZGiGg2eSc87BM1McOju1HDbrM5Xs3dLMXcuBs5WtbXWecySVhdeAbJ2JLwAAEGtJREFUSlKJ1WeqeOutHbz11o7lbbPZPC+enUpD6RTPn5nk2y+NLveUtjdkuPOWZvZsaeaOnibu6G3itu5Gr7eStCzGyMh0lkND0xwZnubQ0DQvnJniyPA0+fRc0lRbxb5bWviFB7azr6+FO29pYWdng9enS7rhGEAlqYQaaqp484523ryjfXnb/GKBQ0PJ8LiDg1McPDPJF77/Motpr0VFgB2dDezubeL2niZ29zZxR28z/e31fpmUbnKz2TyHh6c5PJQ8Dg1NcXhoevmac0hGU+zZ0sxDd3QloyluaWFbuz2bkjYGA6gklVldppJ7+tu4p79teVu+UOTlsbn0S+cUh4eTHo6vHhxaHsJbW13Bru4mdnU3MtDdyEBXAwNdjWzvaCBT5YQh0kYyPrvIsZEZjp6b4di5mWR9ZIZT4/PL+9RnKrm9p4n33NnLHb3JCIndvc20N2TWseaStDpeAypJN7C5xTwvDc8kwTTtFTk2MsPZyYXlfSorAv3t9cuBdKArCai3dTXSUu/kItJ6KRQjp8/PcWxkhmPnZpOwOZI8VvZo1lRVcGtX8qPS0lD83b3NbG2ro8JRD5I2KK8BlaQNqD5TxV3bklu7rDSTzXNiZHb5y+zSF9xvHxllsVBc3q+jIcOOzga2d9SzvT1ddtSzvaOBtvpqh+xJq5TNFzh9fp6TY3O8MjbLy2NznBxP1k+dn18eWg/J/8eB7kYe3bcl+cEo/aGor9WgKWnzsAdUkm4il+pxeWV8llfG5i7oNYVk0pKlMLq9PQmmW9vq6WutY0trLTVVToQkAUwv5HhlOVgm4XLp+ZnJeVZ+larPVNKf/n/a0ZGOSuhu4NbORtocOitpE7EHVJI2gcqKkATKjgYe3n1h2UKuwKn0C/TLY7OcHJ/j5bE5nh+c5LGDQ8uzaQKEAF2NNfS11dHXWkdfWx1b02Vfaz19bXU01vgnRBtfNl9geDLL4MQ8ZyfnOTMxz5nJBc5MzHN2IllOZ/MXHNPRkKG/o5637Ghje8fW5ZEF/e0NdDZmHFkgSVfgtwdJ2iRqqyvZ1dPErp6m15TlC0XOTCxwemKOwfPzDE7MLy8PDE7y2PND5AoXjphpqatmS0stPc219DTX0NtcS3dzLb3N6baWGjoaapy5V+smmy8wOrPIuakFzqah8kwaKs9OzjM4scDoTPY1x7U3ZLiltZb+jnruv7WdW1rr6G+vp7+jnv72eppqvbZakq5XSQNoCOFR4NNAJfC5GOOnLipvAz4PDAALwC/GGA+mZb8G/GsgAgeAj8UYF0IIdwOfAWqBPPDvYow/KuXnkKSbXVVlRfLluqP+kuXFYmRkJsvpC8LpHEOTWYanFnjx7BSjM1mKF13VUVkR6G6qobu5lp6mGnrTwNrZmKGjoYaOxgydjTV0NtZQl3HIr66uUIyMzy4yMp1lZCabLJceM1lGphcYmc4yOrPI5HzuNcfXZyq5pbWOLS217NnSzJaWZMh5X7ptS0udbVGSSqhk14CGECqBI8C7gNPAU8DPxRhfWLHP7wIzMcb/EkLYDfxRjPGdIYQ+4LvA3hjjfAjhL4GvxBi/EEJ4HPiDGONXQwjvA34jxvjQleriNaCSVHr5QpHRmUWGpxYYmlrgXLocnkpC6nC6fqlQAEkw6EiDaWdjhvaGDG31GVrqq2mrz9BaV/3qerqsrTYobGQxRmayeSbmcpyfW2R8dpGJuVy6XOT8XI7xuWR9fDbHyHSW8dnX/tAB0JCppKup5tVHY7LsTJdbWpLh5M11VQ6RlaQyWI9rQO8DjsYYj6cV+AvgQ8ALK/bZC/x3gBjjoRDCjhBCz4q61YUQckA9cCbdHoHmdL1lxXZJ0jqqqqygt6WW3pZa7rrCfgu5AmOzi4xOZxmbTXqqxmYWGZvJJttnspyZWODA4CTn53IXzCJ6sZqqiuUw2lJXvbzeVFtFQ00Vjemj4aJlsl5JY22Vky2twmK+yEw2z/RCjumFPFMLOWYW8kwvvLptOnvR83Q5MZ9jYm7xNUO7l4QArXXVtKU/RPS11nL3tpblYPlq0KylsylDfcariiRpIyjl2boPOLXi+WngrRft8yzwj4HvhhDuA7YDW2OMPwkh/B5wEpgHHo8xPp4e86vAY2l5BfDgpd48hPBx4OMA/f39a/OJJEmrVltdmUxs1Fp3TfvPLxaYmF9c7iWbnMsxMf/q+vm5pGxiPseJ0VmenptgZiHPfK5wTa9fXRloqKmiIVO1HFwbaqqoraqgproyXVZQU1VJTVUFtdXJsmapfEVZTVUlmaoKKisCVRUhWVYurVe8um15WUFlZWDlZbIXD0xa+fTiUUvFIuSKRfKFSK5QJF+M5JeXcbksXyiSS8tyhUg+3Z7NF5hfLDCfK7KQK7CQKzC/vCwyv1i4YPt8rkA2V2Q+V2A2myd7hR8HltRUVdBUW01zbfLv21hbRXdTbfJjQUOGtvpqWusztNdnaGtIfkBoq8/QXFft9cOSdBMqZQC91F+Ni3/m/BTw6RDCMyTXee4H8um1oR8CdgITwF+FEH4+xvhF4N8CvxZj/JsQwj8D/gR45DVvFONngc9CMgR3jT6TJKnM6jKV1GXq2NJybYF1SaEYmV3MM7OQZzabZyabZzZbYCabYyZbWN6WbL9wOTmf41yuwGI+CWbZfJFsup6/1PjPm0RtdRKw69JHbRqw6zKVtNZXp8+TsvpMZRIoa6poqq2mqXblMllvrKkiU1Wx3h9LknQDKWUAPQ1sW/F8KxcNl40xTgEfAwjJBRkn0sd7gBMxxpG07G9Jejq/CHwU+GT6En8FfK50H0GStFFVVgSaa6tpXuMZS/OFIouFItlckYV80iO4FE6XAmq+GCmkvYyF5ecrti89T8sLMV7wq+3FlyiGFaUry0IIVFcmPalVlSvWKwJVlem2i8tWLOtWBMqaqgoq7HGUJJVYKQPoU8CuEMJOYBD4CPDPV+4QQmgF5mKMiyQz3n47xjgVQjgJ3B9CqCcZgvtOYGkWoTPATwPfAh4GXirhZ5Ak6QJJsKugPrPeNZEkaeMpWQCNMeZDCL8CPEZyG5bPxxifDyF8Ii3/DLAH+LMQQoFkcqJfSst+GEL4a+Bpklut7CcdTgv8G5Jhu1Ukt275eKk+gyRJkiRp7ZTsNiw3Em/DIkmSJEnlc7nbsDgzgCRJkiSpLAygkiRJkqSyMIBKkiRJksrCACpJkiRJKgsDqCRJkiSpLAygkiRJkqSyMIBKkiRJksrCACpJkiRJKgsDqCRJkiSpLAygkiRJkqSyMIBKkiRJksrCACpJkiRJKgsDqCRJkiSpLAygkiRJkqSyMIBKkiRJksrCACpJkiRJKgsDqCRJkiSpLAygkiRJkqSyMIBKkiRJksrCACpJkiRJKgsDqCRJkiSpLAygkiRJkqSyMIBKkiRJksrCACpJkiRJKgsDqCRJkiSpLAygkiRJkqSyMIBKkiRJksrCACpJkiRJKgsDqCRJkiSpLAygkiRJkqSyCDHG9a5DyYUQRoBX1rseN6FOYHS9KyG9DrZZbTS2WW00tlltNLbZ0tkeY+y6eOOmCKAqjRDCj2OMb17vekjXyjarjcY2q43GNquNxjZbfg7BlSRJkiSVhQFUkiRJklQWBlCtxmfXuwLS62Sb1UZjm9VGY5vVRmObLTOvAZUkSZIklYU9oJIkSZKksjCASpIkSZLKwgCqKwohtIcQvh5CeCldtl1mv0dDCIdDCEdDCL95ifJfDyHEEEJn6WutzWy1bTaE8LshhEMhhOdCCP8vhNBavtprM7mG82YIIfyPtPy5EMK913qsVArX22ZDCNtCCN8MIbwYQng+hPDJ8tdem9FqzrNpeWUIYX8I4cvlq/XNzwCqq/lN4O9jjLuAv0+fXyCEUAn8EfBeYC/wcyGEvSvKtwHvAk6Wpcba7FbbZr8O7IsxvhE4AvxWWWqtTeVq583Ue4Fd6ePjwB+/jmOlNbWaNgvkgf8QY9wD3A/8sm1WpbbKNrvkk8CLJa7qpmMA1dV8CPjTdP1PgX90iX3uA47GGI/HGBeBv0iPW/IHwG8AznilclhVm40xPh5jzKf7PQlsLXF9tTld7bxJ+vzPYuJJoDWEsOUaj5XW2nW32Rjj2Rjj0wAxxmmSL/R95ay8NqXVnGcJIWwF3g98rpyV3gwMoLqanhjjWYB02X2JffqAUyuen063EUL4IDAYY3y21BWVUqtqsxf5ReCra15D6dra4OX2udb2K62l1bTZZSGEHcA9wA/XvIbShVbbZv+QpAOlWKoKblZV610Brb8QwhNA7yWKfvtaX+IS22IIoT59jXdfb92kSylVm73oPX6bZNjYn7++2knX5Kpt8Ar7XMux0lpbTZtNCkNoBP4G+NUY49Qa1k26lOtusyGEDwDnYow/CSE8tOY12+QMoCLG+MjlykIIw0vDZ9IhCecusdtpYNuK51uBM8AAsBN4NoSwtP3pEMJ9McahNfsA2nRK2GaXXuOjwAeAd0ZvlqzSuGIbvMo+mWs4Vlprq2mzhBCqScLnn8cY/7aE9ZSWrKbNfhj4YAjhfUAt0BxC+GKM8edLWN9NwyG4upovAR9N1z8K/N0l9nkK2BVC2BlCyAAfAb4UYzwQY+yOMe6IMe4g+U9+r+FTJXbdbRaSGfOA/wh8MMY4V4b6anO6bBtc4UvAL6SzNN4PTKbDyq/lWGmtXXebDcmv0H8CvBhj/P3yVlub2HW32Rjjb8UYt6bfXz8CfMPwuXbsAdXVfAr4yxDCL5HMYvtPAUIItwCfizG+L8aYDyH8CvAYUAl8Psb4/LrVWJvdatvs/wJqgK+nPfdPxhg/Ue4PoZvb5dpgCOETaflngK8A7wOOAnPAx6507Dp8DG0iq2mzwNuAfwkcCCE8k277TzHGr5TzM2hzWWWbVQkFR5dJkiRJksrBIbiSJEmSpLIwgEqSJEmSysIAKkmSJEkqCwOoJEmSJKksDKCSJEmSpLIwgEqSJEmSysIAKknSGgghdIQQnkkfQyGEwXR9JoTwv0vwfl8IIZxYuqfdJcpn0uXAUj3Wug6SJL1e3gdUkqQ1FkL4z8BMjPH3SvgeXwC+HGP868uUz8QYGy/3XJKk9VC13hWQJOlmFkJ4CPj1GOMH0mC6E9gC3A78e+B+4L3AIPCzMcZcCOFNwO8DjcAo8K9ijGev8j47gf9L8rf9a6X5NJIkrY5DcCVJKq8B4P3Ah4AvAt+MMb4BmAfeH0KoBv4n8OEY45uAzwP/7Rpe99PAH8cY3wIMlaTmkiStkj2gkiSV11fTXs4DQCWv9lYeAHYAdwD7gK+HEEj3uWLvZ+ptwD9J1/8P8DtrWGdJktaEAVSSpPLKAsQYiyGEXHx1MoYiyd/lADwfY3zgOl7biR0kSTc0h+BKknRjOQx0hRAeAAghVIcQ7ryG474HfCRd/xelqpwkSathAJUk6QYSY1wEPgz8TgjhWeAZ4MFrOPSTwC+HEJ4CWkpYRUmSrpu3YZEkaQO62m1YLrG/t2GRJK07e0AlSdqYJoH/GkL4xJV2CiEMhBCeAYbLUy1Jki7PHlBJkiRJUlnYAypJkiRJKgsDqCRJkiSpLAygkiRJkqSyMIBKkiRJksri/wOVZMVm8dry+wAAAABJRU5ErkJggg==\n",
"text/plain": [
""
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"plot_lc(time, flux1);"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Example 2: multiple passbands "
]
},
{
"cell_type": "code",
"execution_count": 84,
"metadata": {},
"outputs": [],
"source": [
"ldm = LDTkLDModel((3500, 150), (4.5, 0.1), (0.0, 0.1), [sdss_g, sdss_z], frozen=True,\n",
" cache='/Users/hannu/Dropbox/data/Husser2013_bin50/')"
]
},
{
"cell_type": "code",
"execution_count": 85,
"metadata": {},
"outputs": [],
"source": [
"lcids = zeros(time.size, int)\n",
"lcids[time.size//2:] = 1"
]
},
{
"cell_type": "code",
"execution_count": 86,
"metadata": {},
"outputs": [],
"source": [
"tm = SwiftModel(ldm)\n",
"tm.set_data(time, lcids=lcids, pbids=[0,1])"
]
},
{
"cell_type": "code",
"execution_count": 87,
"metadata": {},
"outputs": [],
"source": [
"flux1 = tm.evaluate(k=0.1, ldc=[None], t0=0.0, p=1.0, a=4.2, i=0.5*pi, e=0.0, w=0.0)"
]
},
{
"cell_type": "code",
"execution_count": 69,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAA6AAAAGoCAYAAABYCYj1AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nOzdeXjV5Z3//9c7JxsJS8gKhEACJoGwIyCgLIpWXOo22tY61dqqtV91bGe62GXa6UwXW/U33+6ObZ2p/Xa0rVu1rkgFVEAW2RIgkLCGLSsEErKe+/dHjm3EAAFyzucsz8d1cSXns5y8jtCrvLjvz32bc04AAAAAAARbnNcBAAAAAACxgQIKAAAAAAgJCigAAAAAICQooAAAAACAkKCAAgAAAABCIt7rAKGQmZnp8vPzvY4BAAAAADFh7dq1tc65rBOPx0QBzc/P15o1a7yOAQAAAAAxwcx293ScKbgAAAAAgJCggAIAAAAAQoICCgAAAAAICQooAAAAACAkKKAAAAAAgJCggAIAAAAAQoICCgAAAAAIiaAVUDN73Myqzaz0JOfNzH5iZhVmttHMpnY7t9DMygPnHuh2/CEz2xq4/jkzSwtWfgAAAABA3wrmCOj/SFp4ivNXSCoM/LpL0i8lycx8kn4eOF8i6WYzKwncs0jSeOfcREnbJH0tKMkBAAAAAH0uaAXUObdMUv0pLrlW0hOuy0pJaWY2VNIMSRXOuR3OuTZJTwWulXPudedcR+D+lZKGBys/AAAAAKBvefkMaK6kvd1eVwWOnez4iT4j6ZWTvbmZ3WVma8xsTU1NTR/EBQAAAACcCy8LqPVwzJ3i+N9vNPuGpA5Jvz/ZmzvnHnPOTXPOTcvKyjqnoAAAAACAcxfv4c+ukpTX7fVwSfslJZ7kuCTJzG6TdLWkBc65DxRTAAAAAED48nIE9AVJtwZWw50p6Yhz7oCk1ZIKzazAzBIlfSJwrcxsoaSvSrrGOdfsVXAAAAAAwJkL2giomT0pab6kTDOrkvRtSQmS5Jx7VNLLkq6UVCGpWdLtgXMdZnavpNck+SQ97pwrC7ztzyQlSVpkZpK00jl3d7A+AwAAAACg7wStgDrnbj7NeSfpnpOce1ldBfXE4+f1TToAAAAAQKh5OQUXAAAAABBDKKAAAAAAgJCggAIAAAAAQoICCgAAAAAICQooAAAAACAkKKAAAAAAgJCggAIAAAAAQoICCgAAAAAICQooAAAAACAkKKAAAAAAgJCggAIAAAAAQoICCgAAAAAICQooAAAAACAkKKAAAAAAgJCggAIAAAAAQoICCgAAAAAICQooAAAAACAkKKAAAAAAgJCggAIAAAAAQoICCgAAAAAICQooAAAAACAkKKAAAAAAgJCggAIAAAAAQoICCgAAAAAICQooAAAAACAkKKAAAAAAgJCggAIAAAAAQoICCgAAAAAICQooAAAAACAkKKAAAAAAgJCggAIAAAAAQoICCgAAAAAICQooAAAAACAkKKAAAAAAgJCggAIAAAAAQoICCgAAAAAICQooAAAAACAkKKAAAAAAgJAIWgE1s8fNrNrMSk9y3szsJ2ZWYWYbzWxqt3MLzaw8cO6BbsfTzWyRmW0PfB0crPwAAAAAgL4VzBHQ/5G08BTnr5BUGPh1l6RfSpKZ+ST9PHC+RNLNZlYSuOcBSYudc4WSFgdeAwAAAAAiQHyw3tg5t8zM8k9xybWSnnDOOUkrzSzNzIZKypdU4ZzbIUlm9lTg2s2Br/MD9/9W0hJJXz1dlk6/0+Hmtg8dN9mHL+7pUO8uC1z74TM9Xdvze/b8rj1d29vrenrPBJ/1mBMAAADAmXPO6Whrh/x+53WUsBe0AtoLuZL2dntdFTjW0/ELAt/nOOcOSJJz7oCZZffmB20+0KjJ/77o3BNHiawBSbpqwlDdOXeUctP6eR0HAAAAiEhNrR16+PVyPbVqr463d3odJyJ4WUB7GoJzpzh+Zm9udpe6pvYqM7dA3/5oyQfOux7esacf4nq68CR6fs8PH+ztzz7X9+z5/Zw2H2jU79/drSdX7dGXLy/WZy4sUFwcI6IAAABAb+2tb9atj6/SrromXT8lV2OGDFCCjzVe3/eZH/Z83MsCWiUpr9vr4ZL2S0o8yXFJOmRmQwOjn0MlVZ/szZ1zj0l6TJKmTZvmbr+woC+zR7x9h4/r238u1Xdf2qINVUf00I0TlZzg8zoWAAAAEPaqGpp146PL1dLu15N3ztTMURleRwo7nznJcS8r+guSbg2shjtT0pHA9NrVkgrNrMDMEiV9InDt+/fcFvj+Nkl/DnXoaJGb1k+/unWaHrhijF7csF/3/u86dXT6vY4FAAAAhLWm1g7d8ds1am7r1FN3UT7PVNBGQM3sSXUtGJRpZlWSvi0pQZKcc49KelnSlZIqJDVLuj1wrsPM7pX0miSfpMedc2WBt31Q0h/N7LOS9ki6KVj5Y4GZ6e55o5WS6NO3/lymrz6zSQ/fNJEFigAAAICT+O5LW1R+6Kie+MwMjR060Os4ESeYq+DefJrzTtI9Jzn3sroK6onH6yQt6JOA+JtbZ+WrvqlN//eN7Ro7dIDumDPK60gAAABA2HmzvFpPrtqjz80bpTmFWV7HiUg8JQtJ0v0LCnX5uBw9+MpWrd1d73UcAAAAIKw0tXboa89sUlFOf33x0iKv40QsCigkdU3H/dGNk5Q7uJ/u+f061R1r9ToSAAAAEDYeW7ZDBxtb9IMbWLzzXFBA8TeD+iXoF7dMVX1Tm77xXOkZbUEDAAAARKtDjS16bNkOXTVxqM4fOdjrOBGNAooPGDdskL54WZFeLTuolzYd8DoOAAAA4LlHXi9Xh9+vr14+xusoEY8Cig+5c06BJg0fpH99vlT1TW1exwEAAAA8s/Vgo/60tkq3zcrXiIwUr+NEPAooPiTeF6cf3ThJjS0devj1cq/jAAAAAJ756V8rlJoYr3svOc/rKFGBAooeFQ8ZoFtnjdSTq/aodN8Rr+MAAAAAIVdRfVQvbzqgW2eNVFpKotdxogIFFCf1hUuLlJ6SqO+8WMaCRAAAAIg5P3+zUsnxPn32ogKvo0QNCihOalC/BH3xsiKt3tWgJeU1XscBAAAAQmZXbZP+vH6fbrlghDL6J3kdJ2pQQHFKH5uWp7z0fnr49XJGQQEAABAz/mvZDsX74nTX3FFeR4kqFFCcUmJ8nL6woEhl+xv1WtlBr+MAAAAAQXe4uU3PravS9ZNzlT0w2es4UYUCitO6bkquRmel6pHXt6nTzygoAAAAotsf1+xVS7tft83O9zpK1KGA4rR8caYvXlak7dXH9ErpAa/jAAAAAEHT6Xd6YsVuzchPV8mwgV7HiToUUPTKleOHalRWqn7xZiXPggIAACBq/XVrtaoajjP6GSQUUPRKXJzp7rmjtflAo5Ztr/U6DgAAABAUv12+S0MGJusj43K8jhKVKKDoteum5GrIwGT94s0Kr6MAAAAAfW5PXbPerqjVzTNGKMFHVQoG/qui1xLj43THnAK9u7Nea3c3eB0HAAAA6FNPv1clM+nGacO9jhK1KKA4IzfPGKG0lAT9cgmjoAAAAIgefr/TM2urdNF5mcpN6+d1nKhFAcUZSU2K16dn5+uNLdWqqD7mdRwAAACgTyyvrNO+w8d107Q8r6NENQooztg/zhypRF+cfrdil9dRAAAAgD7xp7V7NTA5Xh8pYfGhYKKA4oxl9k/S1ROH6um1VTra0u51HAAAAOCcHDnerldLD+q6KblKTvB5HSeqUUBxVm6bna+mtk49s7bK6ygAAADAOXll0wG1dvh14/ksPhRsFFCclUl5aZqcl6YnVuyW3++8jgMAAACctefX79OozFRNyB3kdZSoRwHFWfv07HztqG3SWxW1XkcBAAAAzsr+w8f17s56XTclV2bmdZyoRwHFWbtiwhClpybqD6v3eB0FAAAAOCsvbNgv56RrJw/zOkpMoIDirCXF+3T9lFwt2nxIdcdavY4DAAAAnLHn1+3TlBFpGpmR6nWUmEABxTn5+PQ8tXc6Pbdun9dRAAAAgDNSfvCoth48qusm53odJWZQQHFOinIGaMqINP1h9V45x2JEAAAAiBzPr98nX5zpqolDvY4SMyigOGcfn5an7dXH9N6ew15HAQAAAHrF73d6Yf1+zSnMVGb/JK/jxAwKKM7Z1ZOGKSXRpz+u3ut1FAAAAKBX1u09rH2Hj7P4UIhRQHHO+ifF6+qJQ/Xixv061trhdRwAAADgtF4tPaAEn2nB2Byvo8QUCij6xMen56m5rVMvbzrgdRQAAADglJxzennTQc0pzNLA5ASv48QUCij6xNQRg5WfkaLn3mM1XAAAAIS30n2N2nf4uBaOH+J1lJhDAUWfMDPdMHW4Vuyo077Dx72OAwAAAJzUK6UH5IszXcb025CjgKLPXD+la/+k59kTFAAAAGHKOadXSg9q1qgMDU5N9DpOzKGAos/kpadoRn66nn2vij1BAQAAEJbKDx3Vztompt96hAKKPnXD1FxV1jRpY9URr6MAAAAAH/LKpoMyky4fRwH1AgUUferKiUOVGB+nZ9+r8joKAAAA8CGvlh7U9Px0ZQ1I8jpKTKKAok8NTE7QZSU5emHDfrV1+L2OAwAAAPzN7romlR86yuinh4JaQM1soZmVm1mFmT3Qw/nBZvacmW00s1VmNr7bufvNrNTMyszsC92OTzazlWa23szWmNmMYH4GnLkbpuSqobldb1fUeB0FAAAA+Js3tlRLEqvfeihoBdTMfJJ+LukKSSWSbjazkhMu+7qk9c65iZJulfTjwL3jJd0paYakSZKuNrPCwD0/kvQd59xkSd8KvEYY6drQN15/2XDA6ygAAADA3yzeckiF2f01IiPF6ygxK5gjoDMkVTjndjjn2iQ9JenaE64pkbRYkpxzWyXlm1mOpLGSVjrnmp1zHZKWSro+cI+TNDDw/SBJ+4P4GXAWEuPjtHD8EL2++ZBa2ju9jgMAAACosaVdq3bW69ISRj+9FMwCmitpb7fXVYFj3W2QdIMkBabSjpQ0XFKppLlmlmFmKZKulJQXuOcLkh4ys72SHpb0tZ5+uJndFZiiu6amhqmgoXb1xGE61tqhpdv4bw8AAADvLS2vUYff6dKx2V5HiWnBLKDWw7ETN4d8UNJgM1sv6T5J6yR1OOe2SPqhpEWSXlVXUe0I3PN5SV90zuVJ+qKk3/T0w51zjznnpjnnpmVlZZ3zh8GZmT06Q+mpiXpxAwPUAAAA8N4bWw4pPTVRk/MGex0lpgWzgFbp76OWUtfI5gfaiHOu0Tl3e+B5zlslZUnaGTj3G+fcVOfcXEn1krYHbrtN0rOB7/+krqm+CDPxvq5puIu3VKu5reP0NwAAAABB0tHp15LyGl1cnC1fXE/jZAiVYBbQ1ZIKzazAzBIlfULSC90vMLO0wDlJukPSMudcY+BcduDrCHVN030ycN1+SfMC31+ivxdThJmrJw7V8fZO/XVrtddRAAAAEMPW7G7QkePtTL8NA/HBemPnXIeZ3SvpNUk+SY8758rM7O7A+UfVtdjQE2bWKWmzpM92e4tnzCxDUruke5xzDYHjd0r6sZnFS2qRdFewPgPOzQUFGcoakKS/bDigqycO8zoOAAAAYtTiLYeU6IvTnCIezfNa0AqoJDnnXpb08gnHHu32/QpJhSfeFzg35yTH35Z0fh/GRJD44kxXjh+ip1bv1bHWDvVPCuofNwAAAKBHb2yp1gWj0vn7aBgI5hRcQFdMGKrWDr+WlrMaLgAAAEKvsuaYdtY26TK2XwkLFFAE1fT8dGWkJurVsoNeRwEAAEAMWrzlkCTpkjE8/xkOKKAIKl+c6bKSHL25tVqtHZ1exwEAAECMeWNLtcYMGaDhg1O8jgJRQBECl48fomOtHVpeUed1FAAAAMSQI83tWru7QZeOZfptuKCAIuhmj85Q/6R4vVrKNFwAAACEztsVter0O108htVvwwUFFEGXFO/TJWOytWjLIXX6nddxAAAAECOWbqvWwOR4TRqe5nUUBFBAERILxw9RfVObVu+q9zoKAAAAYoBzTku31WhOYZbifdSecMHvBEJiXlGWEuPjmIYLAACAkCg/dFSHGls1r4jpt+GEAoqQSE2K19zCLL1edlDOMQ0XAAAAwfX+PvRzKaBhhQKKkFk4foj2H2nRpn1HvI4CAACAKLd0W43GDBmgIYOSvY6CbiigCJlLx2bLF2dMwwUAAEBQNbV2aPWueqbfhiEKKEImLSVRM0el67UyCigAAACCZ0Vlndo7HQU0DFFAEVILxw1RZU2TKqqPeh0FAAAAUWrpthqlJPp0fv5gr6PgBBRQhNRlJUMkiWm4AAAACArnnJZsq9bs0RlKivd5HQcnoIAipIYMStbkvDQt2nzI6ygAAACIQrvqmrW3/jjTb8MUBRQhd+nYbG2oOqLqxhavowAAACDKLC2vliTNK8r2OAl6QgFFyC0YmyNJ+uvWao+TAAAAINos2VajgsxUjchI8ToKekABRciNGTJAuWn9tJgCCgAAgD7U0t6plTvqmH4bxiigCDkz0yVjsvX29lq1tHd6HQcAAABRYt2ew2pp9+ui8zK9joKToIDCEwvGZut4e6dWVNZ5HQUAAABRYkVlrXxxpgtGpXsdBSdBAYUnZo7KUEqiT29sYTVcAAAA9I3llXWakDtIA5ITvI6Ck6CAwhPJCT7NKczUX7dWyznndRwAAABEuKbWDq3fe1izR2d4HQWnQAGFZxaMzdGBIy0q29/odRQAAABEuFW76tXhd5o9muc/wxkFFJ65uDhbZtLiLayGCwAAgHOzorJOib44nT9ysNdRcAoUUHgma0CSJg1P0+KtPAcKAACAc7O8slZTRqSpX6LP6yg4BQooPHXp2GxtrDqi6sYWr6MAAAAgQh1ublPZ/kam30YACig8tWBsjiTpr1uZhgsAAICzs3JHvZyTZp/HAkThjgIKT40ZMkC5af30Bs+BAgAA4CytqKxVvwSfJg1P8zoKToMCCk+ZmeYXZ2l5Za1aOzq9jgMAAIAItLyyTtML0pUYT70Jd/wOwXMXF2erua1Ta3Y1eB0FAAAAEab6aIu2Vx9j/88IQQGF52afl6FEX5ze5DlQAAAAnKEVlXWSRAGNEBRQeC4lMV4zCtK1ZFuN11EAAAAQYVZU1mlAcrzGDRvkdRT0AgUUYWF+cZYqqo9pb32z11EAAAAQQZZX1mnmqAz54szrKOgFCijCwvzibEliFBQAAAC9tre+WXvqm5l+G0EooAgLo7NSlZfeT0vLeQ4UAAAAvbNix/vPf2Z6nAS9RQFFWDAzzS/K1jsVdWzHAgAAgF5ZUVmnjNREFeX09zoKeokCirAxvzhLx9s7tWpnvddRAAAAEOacc1peWatZozNkxvOfkYICirAxa3SGEuPjtKSc50ABAABwajtqm3SosZXptxEmqAXUzBaaWbmZVZjZAz2cH2xmz5nZRjNbZWbju52738xKzazMzL5wwn33Bd63zMx+FMzPgNBJSYzXBQXpWsJzoAAAADiN5ez/GZGCVkDNzCfp55KukFQi6WYzKznhsq9LWu+cmyjpVkk/Dtw7XtKdkmZImiTpajMrDJy7WNK1kiY658ZJejhYnwGhd3FxtiprmtiOBQAAAKe0orJWwwYla2RGitdRcAaCOQI6Q1KFc26Hc65N0lPqKo7dlUhaLEnOua2S8s0sR9JYSSudc83OuQ5JSyVdH7jn85IedM61Bu5juCyKzC/OkiRGQQEAAHBSfr/Tiso6zRqdyfOfESaYBTRX0t5ur6sCx7rbIOkGSTKzGZJGShouqVTSXDPLMLMUSVdKygvcUyRpjpm9a2ZLzWx6ED8DQqwgM1Uj0lN4DhQAAAAntfXgUTU0tzP9NgIFs4D29E8R7oTXD0oabGbrJd0naZ2kDufcFkk/lLRI0qvqKqodgXviJQ2WNFPSlyX90Xr4Zw8zu8vM1pjZmpoaykykMDNdXJyldypr1dLOdiwAAAD4sOWVtZK6FrFEZAlmAa3S30ctpa6Rzf3dL3DONTrnbnfOTVbXM6BZknYGzv3GOTfVOTdXUr2k7d3e91nXZZUkv6QPLX3lnHvMOTfNOTctKyurrz8bgmh+cbZa2v1sxwIAAIAeraisU0Fmqoal9fM6Cs5QMAvoakmFZlZgZomSPiHphe4XmFla4Jwk3SFpmXOuMXAuO/B1hLqm6T4ZuO55SZcEzhVJSpRUG8TPgRCbOYrtWAAAANCzjk6/3t1Zz+hnhApaAQ0sHnSvpNckbZH0R+dcmZndbWZ3By4bK6nMzLaqa7Xc+7u9xTNmtlnSi5Lucc41BI4/LmmUmZWqa2Gj25xzJ07tRQTrl+jTrFEZLEQEAACAD9m074iOtXbw/GeEig/mmzvnXpb08gnHHu32/QpJhSe5d85JjrdJ+sc+jIkwNL84S995cbP21DVrBEtrAwAAIOD9/T9njqKARqJgTsEFztrFxdmSpCXbGAUFAADA362orNOYIQOU2T/J6yg4CxRQhKX8zFTlZ6Toza0UUAAAAHRp7ejU6l08/xnJKKAIW/OLs7ViRx3bsQAAAECStG7PYbV2+DV79Ic2wUCEoIAibM0vzlJLe9cqZwAAAMDyyjrFmTSjIN3rKDhLFFCErZmjMpScEMc0XAAAAEiSVlTWakLuIA3ql+B1FJwlCijCVnKCTzNHZWjZNvYDBQAAiHXNbR1at+ewZjH9NqJRQBHW5hVlaUdtk/bUNXsdBQAAAB5avatBHX7H/p8RjgKKsDavKEuStHQ7o6AAAACxbHllrRJ8pmn5g72OgnNAAUVYK8hMVV56Py0tp4ACAADEshWVdZqSN1gpifFeR8E5oIAirJmZ5hdla3llrdo6/F7HAQAAgAeONLerdN8R9v+MAhRQhL15RVlqbuvUmt1sxwIAABCL3t1ZJ78Tz39GAQoowt6s0RlK8BnTcAEAAGLU8so6JSfEafKINK+j4BxRQBH2UpPiNT0/XUvZjgUAACAmrais0/T8dCXF+7yOgnNEAUVEmFeUpa0Hj+rgkRavowAAACCEao62qvzQUZ7/jBIUUESEecVd27EsYxQUAAAgpqzcUSdJmj060+Mk6AsUUESE4pwByhmYxDRcAACAGLO8sk4DkuI1fthAr6OgD1BAERHMTPOKsvTW9hp1dLIdCwAAQKxYUVmrC0alK95HdYkG/C4iYswrylZjS4c2VB32OgoAAABCYN/h49pV16xZTL+NGhRQRIyLCjPlizMtYTsWAACAmLCi8v3nP1mAKFpQQBExBvVL0JS8NJ4DBQAAiBHLK2uVnpqo4pwBXkdBH6GAIqLMK8rSxqojqj3W6nUUAAAABJFzTisq6zRrVIbi4szrOOgjFFBElPe3Y3l7e63HSQAAABBMu+qadeBIC/t/RpleFVAzK+nh2Pw+TwOcxvhhg5Semsg0XAAAgCi3vLJrwIHnP6NLb0dA/2hmX7Uu/czsp5J+EMxgQE/i4kxzCzO1bFuN/H7ndRwAAAAEyfLKOg0ZmKyCzFSvo6AP9baAXiApT9JySasl7Zd0YbBCAacyrzhLdU1tKtvf6HUUAAAABIHf77Sysk6zR2fIjOc/o0lvC2i7pOOS+klKlrTTOecPWirgFOYUdj0HuqS82uMkAAAACIZt1UdV19TG859RqLcFdLW6Cuh0SRdJutnMng5aKuAUMvsnaeLwQTwHCgAAEKWWV3Tt/0kBjT69LaCfdc59yznX7pw76Jy7VtKfgxkMOJV5RVl6b0+DjjS3ex0FAAAAfWx5ZZ1GZqRo+OAUr6Ogj/W2gFab2YjuvyQtDWYw4FTmFWXJ76R3KtmOBQAAIJp0dPr17o46Vr+NUvG9vO4lSU6SqesZ0AJJ5ZLGBSkXcEqT89I0IDleS8trdOWEoV7HAQAAQB8p29+oo60dmjU60+soCIJeFVDn3ITur81sqqTPBSUR0AvxvjjNKczU0m01cs6xOhoAAECUWF4ZeP5zFCOg0ai3U3A/wDn3nroWJAI8M68oSwcbW7Tt0DGvowAAAKCPLK+sVVFOf2UNSPI6CoKgVyOgZvbP3V7GSZoqiSVI4am5RX/fjqV4yACP0wAAAOBctXX4tXpXvT4xfYTXURAkvR0BHdDtV5K6ngm9NlihgN4YOqifinMGsB0LAABAlFi/97Ba2v1svxLFevsM6HeCHQQ4G/OLs/T4OzvV1Nqh1KTerqkFAACAcLS8slZm0swCCmi0OuXf2M3sRXWtftsj59w1fZ4IOAPzirL0X8t2aEVlnS4tyfE6DgAAAM7B8so6jR82SINSEryOgiA53ZDRwyFJAZyl8/MHKyXRp6XbaiigAAAAEex4W6fW7WnQZy4s8DoKguh0BXSnc25PSJIAZyEp3qfZozO0ZFs127EAAABEsDW769Xe6Xj+M8qdbhGi59//xsyeCXIW4KzMK8rS3vrj2lXX7HUUAAAAnKXllXWKjzNNz0/3OgqC6HQFtPtw0qgzfXMzW2hm5WZWYWYP9HB+sJk9Z2YbzWyVmY3vdu5+Mys1szIz+0IP937JzJyZZZ5pLkSXeUXZkrq2YwEAAEBkWl5Zp8l5aSwsGeVOV0DdSb4/LTPzSfq5pCsklUi62cxKTrjs65LWO+cmSrpV0o8D946XdKekGZImSbrazAq7vXeepMskMT0YGpGRooLMVLZjAQAAiFBHjrdrU9VhzWb6bdQ7XQGdZGaNZnZU0sTA941mdtTMGk9z7wxJFc65Hc65NklP6cN7h5ZIWixJzrmtkvLNLEfSWEkrnXPNzrkOSUslXd/tvv+U9BWdYSlG9JpXlKWVO+rU0t7pdRQAAACcoVU76+V30uzzmNwY7U5ZQJ1zPufcQOfcAOdcfOD7918PPM1750ra2+11VeBYdxsk3SBJZjZD0khJwyWVSpprZhlmliLpSkl5geuukbTPObfhVD/czO4yszVmtqamhpGxaDevOEst7X6t2lnvdRQAAACcoeWVtUqKj9OUEWleR0GQnW4E9Fz0tBzpiSOWD0oabGbrJd0naZ2kDufcFkk/lLRI0qvqKqodgTL6DUnfOt0Pd8495pyb5pyblpWVdQ4fA6hotI0AACAASURBVJFgZkGGEuPjmIYLAAAQgVZU1ml6frqS4n1eR0GQBbOAVikwahkwXNL+7hc45xqdc7c75yar6xnQLEk7A+d+45yb6pybK6le0nZJoyUVSNpgZrsC7/memQ0J4udABOiX6NMFBekUUAAAgAhTd6xVWw8eZfuVGBHMArpaUqGZFZhZoqRPSHqh+wVmlhY4J0l3SFrmnGsMnMsOfB2hrmm6TzrnNjnnsp1z+c65fHWV3KnOuYNB/ByIEPOKslRRfUxVDWzHAgAAECmWV9ZJEgsQxYigFdDA4kH3SnpN0hZJf3TOlZnZ3WZ2d+CysZLKzGyrulbLvb/bWzxjZpslvSjpHudcQ7CyIjrML+6aas0oKAAAQORYXlmrAUnxmpA7yOsoCIGgbrLjnHtZ0ssnHHu02/crJBWeeF/g3JxevH/+OUZEFBmd1V+5af20tLxGt1ww0us4AAAA6IV3Kup0wagMxfuCOTkT4YLfZUQNM9Pcoiwtr6xTW4ff6zgAAAA4jb31zdpT36wLz2P6bayggCKqzCvK0rHWDr23hxnbAAAA4e6dilpJ0oXs/xkzKKCIKheel6H4OOM5UAAAgAjwTmWdsgYkqTC7v9dRECIUUESVAckJOn/kYC0pp4ACAACEM+ecVlTW6sLRGTIzr+MgRCigiDoXj8nWlgON2n/4uNdRAAAAcBLlh46q9libZjP9NqZQQBF1Lh2bI0lavOWQx0kAAABwMu9UdO3/yfOfsYUCiqgzOitV+RkpemNLtddRAAAAcBLLK2qVn5Gi3LR+XkdBCFFAEXXMTJeOzdGKyjoda+3wOg4AAABO0NHp17s765l+G4MooIhKC8bmqK3Tr7e3sxgRAABAuNlQdUTHWjt04WgKaKyhgCIqTcsfrEH9ErRoM9NwAQAAws37+3/OGp3hcRKEGgUUUSnBF6f5xVl6s7xanX7ndRwAAAB0805FrcYNG6j01ESvoyDEKKCIWpeOzVF9U5vW7WnwOgoAAAACjrd1at2ew6x+G6MooIha84qzFB9nrIYLAAAQRlbvqldbp1+zmX4bkyigiFoDkxN0wah0vcF+oAAAAGHjncpaJfhMMwrSvY4CD1BAEdUWjMlRRfUx7apt8joKAAAAJC2vqNOUvMFKSYz3Ogo8QAFFVLt0bI4kMQoKAAAQBg43t6l0/xHNPo/pt7GKAoqoNiIjRcU5A7SY50ABAAA8t6KyTs6JBYhiGAUUUW/B2Gyt2lWvI83tXkcBAACIacu212hAUrwm56V5HQUeoYAi6l1akqNOv9OSbYyCAgAAeMU5p2XbajVrdIYSfNSQWMXvPKLe5OFpyuyfyHYsAAAAHtpR26R9h49rblGW11HgIQoool5cnOmSMdlaUl6ttg6/13EAAABi0rJtNZKkeRTQmEYBRUz4SMkQHW3p0IoddV5HAQAAiEnLttUoPyNFeekpXkeBhyigiAkXFWaqf1K8Xt54wOsoAAAAMae1o1Mrd9Qz/RYUUMSG5ASfFozN1mubD6q9k2m4AAAAobR2V4OOt3dqbiEFNNZRQBEzrpwwVIeb2/XujnqvowAAAMSUpdtrlOAzzRqd4XUUeIwCipgxryhLqYk+vbSJabgAAAChtGxbraaOGKzUpHivo8BjFFDEjOQEny4Zm6PXyw6qg2m4AAAAIVF9tEVbDjTy/CckUUARY64cP0R1TW16dyfTcAEAAELh7e21kth+BV0ooIgpF4/J1oDkeD29tsrrKAAAADFh2bYaZaQmqmToQK+jIAxQQBFTkhN8um5yrl7edEBHmtu9jgMAABDVOv1OS7fVaG5RluLizOs4CAMUUMScj0/PU2uHX8+v3+d1FAAAgKj23p4GNTS3a8HYbK+jIExQQBFzxucO0vjcgXpy1R4557yOAwAAELXe2HJICT5jASL8DQUUMenmGSO09eBRFiMCAAAIojc2H9IFBRkamJzgdRSECQooYtI/TB2uzP5Jeui1cvn9jIICAAD0tV21TaqsaWL6LT6AAoqYlJzg01cWFmvt7gZWxAUAAAiCN7YckiRdOjbH4yQIJxRQxKwbpw7X9PzB+sErW1Tf1OZ1HAAAgKiyeEu1inL6Ky89xesoCCMUUMSsuDjTd6+boMaWDj30WrnXcQAAAKLGkeZ2rdpVrwWMfuIEFFDEtOIhA3TbrHw9tXqPSvcd8ToOAABAVFiyrVqdfqdLef4TJwhqATWzhWZWbmYVZvZAD+cHm9lzZrbRzFaZ2fhu5+43s1IzKzOzL3Q7/pCZbQ3c85yZpQXzMyD63X9podJTEvVvL5SxLQsAAEAfWLylWhmpiZqcN9jrKAgzQSugZuaT9HNJV0gqkXSzmZWccNnXJa13zk2UdKukHwfuHS/pTkkzJE2SdLWZFQbuWSRpfOCebZK+FqzPgNgwqF+Cvnx5sdbsbtALG/Z7HQcAACCitXf6taS8WhePyZYvzryOgzATzBHQGZIqnHM7nHNtkp6SdO0J15RIWixJzrmtkvLNLEfSWEkrnXPNzrkOSUslXR+47vXAMUlaKWl4ED8DYsRN0/I0IXeQfvDyVjW1dpz+BgAAAPRoza4GNbZ0MP0WPQpmAc2VtLfb66rAse42SLpBksxshqSR6iqUpZLmmlmGmaVIulJSXg8/4zOSXunph5vZXWa2xszW1NTUnNMHQfTzxZn+7ZoSHWxs0S+WVHgdBwAAIGK9seWQEn1xmlOY5XUUhKFgFtCexttPfMDuQUmDzWy9pPskrZPU4ZzbIumH6ppu+6q6iuoHhqXM7BuBY7/v6Yc75x5zzk1zzk3LyuIPP07v/JHpun5Krn61bKd21zV5HQcAACDiOOe0aPMhzRydodSkeK/jIAwFs4BW6YOjlsMlfeABO+dco3PudufcZHU9A5olaWfg3G+cc1Odc3Ml1Uva/v59ZnabpKsl3eJYNQZ96IErxijeZ/ruS1u8jgIAABBxthw4qj31zVo4bojXURCmgllAV0sqNLMCM0uU9AlJL3S/wMzSAuck6Q5Jy5xzjYFz2YGvI9Q1TffJwOuFkr4q6RrnXHMQ8yMG5QxM1r2XnKdFmw9p2TambgMAAJyJV0sPKM6kj4xj/0/0LGgFNLBQ0L2SXpO0RdIfnXNlZna3md0duGyspDIz26qu1XLv7/YWz5jZZkkvSrrHOdcQOP4zSQMkLTKz9Wb2aLA+A2LTZy8q0MiMFH3nxTK1d/q9jgMAABAxXik9qOn56crsn+R1FISpoE7Mds69LOnlE4492u37FZIKT7wvcG7OSY6f15cZgRMlxfv0r1eV6I4n1ui3y3fpjjmjvI4EAAAQ9iqqj2l79TH920dP3HkR+LtgTsEFItaCsdmaV5SlH7+xXTVHW72OAwAAEPZeKzsoSbp8PM9/4uQooEAPzEzf+miJWjv9+vLTG+T3s9YVAADAqbxSekCT89I0dFA/r6MgjFFAgZMYndVf/3rVWC0pr9Hj7+z0Og4AAEDY2l3XpNJ9jbqC0U+cBgUUOIV/nDlSl4/L0Q9f3aoNew97HQcAACAsPfvePplJ10we5nUUhDkKKHAKZqYf/sNEZfVP0n1PrtPRlnavIwEAAIQVv9/p2XVVunB0JtNvcVoUUOA00lIS9ZObp2jf4eP6+nOlco7nQQEAAN63ZneD9tYf1w1Tc72OgghAAQV6YVp+ur54aaFe3LBff1pT5XUcAACAsPHM2iqlJPq0kOc/0QsUUKCXPj//PM0enaFvvVCqiuqjXscBAADw3JHj7Xpx435dNWGoUhLjvY6DCEABBXrJF2f6z49PVmpivO7933Vqae/0OhIAAICn/rRmr5rbOnXb7HyvoyBCUECBM5AzMFkPf2ySth48qu+9tMXrOAAAAJ5p7/Trf5bv0oz8dI3PHeR1HEQICihwhi4uztZdc0fpdyt369XSA17HAQAA8MTvV+5WVcNx3T1/lNdREEEooMBZ+NJHijVp+CB95emN2lvf7HUcAACAkDp4pEU/Xrxds0dn6OLibK/jIIJQQIGzkBgfp5/ePFVO0h2/XaNG9gcFAAAxoqGpTZ/6zbtq6/Dr2x8dJzPzOhIiCAUUOEsjMlL0y1vOV2XNMd3z+/fU3un3OhIAAEBQNbV26NP/s1q765v1q9umqXjIAK8jIcJQQIFzcFFhpr5//QS9tb1W3/pzqZxzXkcCAAAIitaOTt31uzUq3XdEP7t5imaPzvQ6EiIQm/UA5+hj0/O0u75JP3+zUvkZqfrcvNFeRwIAAOhTHZ1+3f/ker1TUaeHb5qkj4wb4nUkRCgKKNAH/uWyYu2ua9YPXtmqvPQUXTlhqNeRAAAA+oRzTt94rlSvlh3Uv15dohvPH+51JEQwCijQB+LiTA/fNEkHjrToi39Yr6GDkjVlxGCvYwEAAJwT55x+8MpW/WHNXt13yXn67EUFXkdChOMZUKCPJCf49NinzlfOwGTd+cQatmcBAAAR75dLK/XYsh361MyR+ufLiryOgyhAAQX6UEb/JP337dPV3ul0+/+s1pHjbM8CAAAi0/++u0c/erVc10wapu9cw3Yr6BsUUKCPjc7qr0f/8XztrmvS3b9bq5b2Tq8jAQAA9JpzTr9+a4e+8fwmXVycpUc+NklxcZRP9A0KKBAEs0Zn6KEbJ2nlzjp9/v+tVWsHJRQAAIS/1o5OfeXpjfruS1t0eckQ/eKW85XgozKg7/CnCQiS66bk6nvXTdCb5TW673/Xqb3T73UkAACAk6o91qpbfvWu/rS2Sv+0oFC/uGWq+iX6vI6FKEMBBYLokxeM0L99tESvbz6kL/5hvTr9zutIAAAAH7J5f6Ou/dk72rTviH568xT982VFTLtFULANCxBkn76wQG2dfn3/5a1KjI/TwzfyHAUAAAgfr5Ud1Bf/sF4DkxP09N2zNWH4IK8jIYpRQIEQuGvuaLW2+/XIom1Kio/T96+fwEpyAADAU36/0y+WVOjh17dpUl6afvWp85U9MNnrWIhyFFAgRO5bUKjWDr9+9maFkuJ9+vZHSyihAADAE1UNzfrynzZqxY46XTd5mB78h4lKTuB5TwQfBRQIoX/5SJFa2jv167d3KjE+Tl+7YgwlFAAAhIxzTk+vrdJ3Xtws55x++A8T9LFpefx9BCFDAQVCyMz0javGqrXDr8eW7ZDf7/T1K8fyTCgAAAi62mOt+tqzm7Ro8yHNKEjXIzdNUl56itexEGMooECImZm+c804+eJMv357p+qa2vSjGyeyxxYAAAiaV0sP6uvPbdKx1g5986qx+syFBfwDODxBAQU8EBdn+vZHS5TZP1EPv75NDc1t+sUtU5WSyP8kAQBA3zlyvF3febFMz763T+NzB+r/+9hkFeUM8DoWYhh/2wU8Yma695JCZfRP0jee26RP/upd/fenp2twaqLX0QAAQITz+52eea9KP3y1XA3NbfqnBYW675LzmHEFz1FAAY/dPGOE0lMTdd+T63Tjo8v1xGcvUG5aP69jAQCACLVh72F9+4Uyrd97WFNGpOm/Pz2dvT0RNvgnECAMXD5uiJ74zAxVN7bqxl8u1/ZDR72OBAAAIkztsVZ99emNuu4X76iq4bgeuWmSnrl7NuUTYYUCCoSJmaMy9IfPzVKH3+nGR1do7e4GryMBAIAI0NHp13+/s1MXP7xEz7xXpTsuKtCbX5qnfzh/OAsNIexQQIEwUjJsoJ79/GwNTknQLb9eqdfLDnodCQAAhCm/3+nP6/fpsv9cpu+8uFmT89L06hfm6BtXlWhAcoLX8YAeUUCBMJOXnqKnPz9bxTkDdNfv1uqni7fLOed1LAAAECacc3q19KCu+PFbuv+p9UqKj9Ovbp2mJz4zQ+dls8ItwhuLEAFhKLN/kv7wuVl64JmNemTRNm09eFQP3TSRbVoAAIhhzjkt3VajR17fpk37jmhUZqp+cvMUXT1hKFNtETH42ywQppITfPrPj09WybCBevCVrdpR26Rf3Xq+hg9O8ToaAAAIIeec3tpeq5/+dbtW72rQ8MH99NCNE3X9lFzFs60KIkxQ/8Sa2UIzKzezCjN7oIfzg83sOTPbaGarzGx8t3P3m1mpmZWZ2Re6HU83s0Vmtj3wdXAwPwPgJTPTXXNH6/FPT1dVQ7Ou/dk7endHndexAABACLR3+vX8un268idv69bHV2lPfbP+47rx+uu/zNdN0/Ion4hIQftTa2Y+ST+XdIWkEkk3m1nJCZd9XdJ659xESbdK+nHg3vGS7pQ0Q9IkSVebWWHgngckLXbOFUpaHHgNRLX5xdl6/p4LNSglQbf8+l39v5W7vY4EAACCpKm1Q795e6fmP7REX/jDerV3+vWjGydq2Vcu1qdmjlRiPMUTkSuYU3BnSKpwzu2QJDN7StK1kjZ3u6ZE0g8kyTm31czyzSxH0lhJK51zzYF7l0q6XtKPAu8xP3D/byUtkfTVIH4OICyMzuqv5++5UPc/uU7ffL5UWw826tsfHacE/vUTAICocPBIi363cpd+t2K3Gls6NKMgXf9+7ThdXJzNM56IGsEsoLmS9nZ7XSXpghOu2SDpBklvm9kMSSMlDZdUKul7ZpYh6bikKyWtCdyT45w7IEnOuQNmlt3TDzezuyTdJUkjRozokw8EeG1gcoJ+fdt0PfRauR5dWqnN+xv1s09O1bC0fl5HAwAAZ8E5pxWVdfrdyt16ffMh+Z3TwnFDdNfcUZoygifNEH2CWUB7+meaE/eSeFDSj81svaRNktZJ6nDObTGzH0paJOmYuopqx5n8cOfcY5Iek6Rp06axhwWihi/O9MAVYzQ+d6AeeGaTrvzJW3rkpklaMDbH62gAAKCXGlva9ezaKv1u5W5V1jQpLSVBd1xUoE9eMEIjM1K9jgcETTALaJWkvG6vh0va3/0C51yjpNslycxM0s7ALznnfiPpN4Fz3w+8nyQdMrOhgdHPoZKqg/gZgLB19cRhGj9skP7P79/TZ3+7Rp+bO0pfuryYKbkAAISxsv1H9Pt39+j5dfvU3NapyXlpeuSmSbpq4lAlJ/i8jgcEXTAL6GpJhWZWIGmfpE9I+mT3C8wsTVKzc65N0h2SlgVKqcws2zlXbWYj1DVNd1bgthck3aau0dPbJP05iJ8BCGv5mal69v/M1vde2qL/WrZDq3fV66efnKpcpuQCABA26pva9Of1+/SnNVXafKBRSfFxunbyMH1qZr4mDB/kdTwgpMy54M1ONbMrJf1fST5Jjzvnvmdmd0uSc+5RM5sl6QlJnepanOizzrmGwL1vScqQ1C7pn51ziwPHMyT9UdIISXsk3eScqz9VjmnTprk1a9ac6hIg4v1l43498MwmxfuMKbkAAHisvdOvpeU1enptlRZvPaT2TqcJuYN04/nDde3kYUpLSfQ6IhBUZrbWOTftQ8eDWUDDBQUUsWJXbZPu+d/3VLa/UXfOKdCXLi9WUjzTeQAACAW/32nd3ga9uOGA/rLxgGqPtSqzf6Kum5yrG6cN15ghA72OCITMyQpoMKfgAgix/MxUPfP52fr+y1v0q7d26q3ttXrkY5M0bhjTewAACAbnnMr2N+rFjfv1lw0HtO/wcSXGx+mS4mzdMDVXF4/JZn0GoBtGQIEo9ebWan31mY2qb2rT/QsK9fn5oxXP/wECAHDOnHPaXn1Mf9l4QH/ZsF87apsUH2eaU5ipj04apstKcjQgOcHrmICnmIJLAUUMOtzcpm/9uUwvbNivScMH6ZGPTdZ52f29jgUAQMTp9Du9t6dBizYf0utlB7Wrrllm0qxRGfropGFaOG6IBqfyXCfwPgooBRQx7KWNB/TN5zepua1TX1k4RrfPzldcXE9b9QIAgPe1tHfq7e21en3zQS3eUq26pjYl+EyzR2fqspIcfaQkR9kDk72OCYQlngEFYthVE4dqesFgff3ZTfqPv2zW62UH9fBNk5SXnuJ1NAAAwkpDU5v+urVar28+qGXbanW8vVMDkuJ18ZhsXVaSo/nFWUyvBc4BI6BADHHO6em1Vfr3FzfL75y+dHmxbp2VLx+joQCAGNXpd9q074je2lajt7bXau2eBnX6nXIGJukjJUN0WUmOZo7KUGI86ygAZ4IpuBRQ4G/2HT6urz+7SUu31WjS8EH6/g0TWCkXABAz9h8+rre212jZ9lq9U1Grw83tkqTxuQM1ryhLHykZogm5g3hcBTgHFFAKKPABzjm9uPGA/v3FMjU0t+uOiwp0/6WFSklkZj4AILrUHWvVqp31WrmjTm9X1KqypkmSlDMwSXMKszSnMFMXnZepjP5JHicFogfPgAL4ADPTNZOGaV5hlh58dYv+a9kOvbTpgL573XjNL872Oh4AAGettlvhXLmjTtsOHZMk9UvwaXpBum6eMUJzi7JUmN1fZoxyAqHECCgASdKqnfX62rMbVVnTpGsmDdM3rx6r7AGs7AcACH81Rz9YOLdXdxXOlESfpuWna+aodM0claEJuYOUwJ7YQEgwBZcCCpxWa0enfrmkUr94s1JJ8XH6pwWFum12PgsvAADChnNO+w4f13t7DmvVzjqt3FGvikDhTP1b4czQzFHpGk/hBDxDAaWAAr22o+aY/uMvm/VmeY1GZaXqW1eXMC0XAOCJ5rYObaw6onV7Dmvdngat23tYNUdbJXUVzukF7xfODI0fNlDxFE4gLFBAKaDAGfvr1kP6j79s0c7aJl06NlvfvKpE+ZmpXscCAESpTr9TZc2xQOFs0Lo9h1V+6Kg6/V1/Xy3ITNWUvDRNGZGmyXmDNXboAAonEKYooBRQ4Ky0dfj13+/s1E8Wb1d7p9P/396dR8d91vcef39n04xmtEvWblve4y1OYrIC2SDN1iRlK7QFynK5SctSbrmFlnPvaS/cc5NbDm1pe4GQBgqh6QmEtGHJRhay4Bgcx47jLYmleJFly5Jt7Zr1uX/MT7Lk2LFjaWas0ed1zu/89plnksfSfPQsv4+/vY1PX7WIWInmMBMRkTOXzjjaDw2ypbOPl/b18XJnH1v39zOSTANQVhJgzdxKL3BWcW5rJdXRUIFLLSKnSwFUAVRkSrr7R7nj4Z3cv3EfNdEQn7lqEX9w0TyNDxURkVNKZxwdPW8Mm8OJbNiMBP2saCpnZXMFq5orWN1SwcK6mJ7DKTKDKYAqgIpMi817j3L7QztY195La3WEL1yzlN9d3aQvCSIignOOQ4Nxdh4YOLYcHOCVgwOMJjNANmwubypnlRc2V3lh06/fIyJFRQFUAVRk2jjnePrVHm5/aAfbu/pZ3ljOX1y7lMuX1Ol5aiIis8RgPMXOA9lwufPAADsO9LPzwABHhpPj19TGSljWUMbShjLOacyGzoV1UY3bFJkFFEAVQEWmXSbj+OlL+/naozvZe3iESxbU8MXrlrGmtbLQRRMRkWkST6V5vWd4PGCOtWruOzIyfk1pyM+S+rLxsLm0oYyl9WXUxEoKWHIRKSQFUAVQkZxJpDL82/rd/OMTr9E7lOD6VQ38+TVLWVgXK3TRRETkNKQzjv1HR+joGRpf2nuG6OgZpPPICN4ktPh9xoLaKEsbxsJmOcsaymiujGgohohMogCqACqSc4PxFHc90853nm5nJJnmljXNfObqxbTp0S0iIgXnnKN3KJENmIeOBcyOniFe7x0mkcqMXxsN+Wmri9JWG6OtppQFdTGWNpSxoC5KScBfwE8hIjOFAqgCqEje9AzGufPpdr6/7nWSaZcNolct0jNERURyLJXO0NU3yt7Dw+w5bunoGWJgNDV+bdBvzK0upa02xoK6KG212WVBbZS6shKN6ReRKVEAVQAVybtDA3G+/atd3LN+N8m04/fOa+azVy1mbk1poYsmIjIjpdIZugfidPWN0Hl0lM4jI+w9MjweODuPjJDKHPtuF/AZLVURWqtLxwPm2NJcGdFkQCKSMwqgCqAiBdM9MMq3nmrnh+t3k8o43nt+M7ddsUhdc0VEJnDOcXQ4SefREbr6Rtl/dIT9fSPsPzpK19ER9h8d4eBAnHRm8ne36miI1upS5laXMrc6wtzq0vH9hvKwQqaIFIQCqAKoSMF194/yzV/t4ofr95BKZ7huVSO3Xb6Qlc0VhS6aiEhOjSbTHBqI0z0wSnd/nO4J2119o17QHBl/VuaYkN9HY2WYxoowTZURmioiNFVGaKwM01wZobEiTFk4WKBPJSJycgqgCqAiZ43ugVHufvZ17nl+N4PxFJcvqeO2KxZyUVu1xhyJyIzhnGMwnsqGyf5soMyGzDjd/aNeyMxu908YeznG7zNqYyEaKyI0VYZpqojQWBmhuTLsHYtQEw1pdlkRmZEUQBVARc46fSNJ7nl+N3c/20HvUILz5lbyibe3ce2KBnUZE5G8c84xlEhzeDBBz1Ccw4MJDg8dv53g8FCc3sEEvUOJSTPHjgkFfMwpK/GWMHPKj23XTdiujobwK1yKSJFSAFUAFTlrjSbT3LdhL3c908Gew8M0VYT58CXz+dCFrVSWhgpdPBGZYVLpDP2jKY4OJ+gbSXJ0JEn/SJKjw9kleywxfmzsmr7hJIn0GwMlQCTopzoaojYWojoaojpaQo23PTFo1peFKY8E1JtDRGY9BVAFUJGzXjrjeGJHN3c/28G69l7CQR/vOb+Fj106n8X1ZYUunojkmHOO0WSGgXiSwdEUg3FvmbA9MJpiaMLxAW/dP5oNkn3DSQbib+zuOlGsJEBFJEhlaXDSujwSpGYsXEaz4bImFqImWkIkpGdfioi8FQqgCqAiM8r2rn6++1wH/7FpP4lUhncsruXjl7Vx+ZI6jYcSyaN0xpFIZUikM+Pr5HH7iVSG0WSakUSakaS3JNKMJtMMe8fGzk/aH7snkWYokWYwnnrDDK8nEvAZZeEAsXCAaChAWThAWThIZSRIRWmQykiIikiAytIQFePHjoXMoLr4i4jknAKoAqjIjNQ7GOfe3+zh++t20z0QZ0FtlD++bD7vPb+FeppzygAAEzJJREFUaEmg0MUTyZl0xvGdZ9o5MpzAuex+xjmcg4xz3n621XDStstuZ5wj492TcUzedo54KkNyLERO3J4YNNPutALhm/H7jNKgn3DIT2nITyToJxzMriOh7BIO+ImV+ImFA8RKgsTCAcpKAsRKAt6xydslAZ+6uIqInOUUQBVARWa0RCrDQy93cfezHWze10dZOMDvr23lDy+ep+eJSlHaur+PG77xLEG/EfD58Bn4fIbPDJ9lg52NbZu37ctu+8wwA5/Z5OsmbIf8PkIB3/g6OLY/4VjIf+x40G+UBI671u8jGPBR4vcR9sLlWLAMB7OBU62NIiKz08kCqJoPRGRGCAV83LymmZvObWLjnqN897kOvvvr17nr2Q4uXVjDhy6cyzUr6ikJaJyWFIdkOvsH4js/vJYrl80pcGlERESmhwKoiMwoZsYF86q4YF4V3f2j3LdhL/f+Zi+fufdFqqMh3n9BC7//tlYW1MUKXVSRKUlnsrOx6jEdIiJSTBRARWTGmlMe5tNXLeZPrljEM6/1cO/6Pdz1bAfffrqdtfOqeP/aFm5Y3URMY0VlBhprAQ0ogIqISBHRtzIRmfF8PuPyJXVcvqSO7v5R7t/YyY9e2MsX79/CXz+4jetWNfD+C1q5qK1aM+jKjDE2+U9AYyhFRKSIKICKSFGZUx7mtisWcuvlC9i45yg/fmEvP93cxU82dtJUEeamNc3ccl4TyxrKC11UkTeV8gKouuCKiEgxUQAVkaI0cazo/7xxBY9uO8B/btrPd55p51u/2sWyhjJuOS87qVFTZaTQxRV5g1Q6OwZUXXBFRKSYKICKSNGLhPzcvKaZm9c00zsY5+dbunjgxU5uf2gHdzy8g4vaqrllTTPXrWqkIhIsdHFFALWAiohIccrpwBIzu9bMdprZa2b2pROcrzKzB8zsJTP7jZmtnHDu82a21cxeNrN7zSzsHV9jZs+b2SYz22BmF+byM4hIcamJlfCRS+bzwJ9cxq/++xV8/l1L6O6P86WfbOFtX/0lt/7gBR5++QDxVLrQRZVZbmwMqJ6jKSIixSRnLaBm5gf+GXg3sA/4rZk96JzbNuGyvwI2Oed+z8yWeddfbWbNwGeB5c65ETO7D/gg8D3g/wJ/45x7yMyu9/avyNXnEJHiNa8mymevXsxnrlrEls4+Hnixk59u7uLhrQcoCwe4etkcfmdFA+9cUkdUM+lKniXTegyLiIgUn1x+o7oQeM051w5gZv8O3AxMDKDLgf8D4JzbYWbzzax+QtkiZpYESoH93nEHjM0eUjHhuIjIGTEzVrdUsrqlki9ffw6/3tXLg5v38/j2g/zHpv2EAj7esaiWa1bU865z6qmJlRS6yDILjM+CqwAqIiJFJJcBtBnYO2F/H3DRcddsBt4DPOt1pZ0HtDjnXjCzrwF7gBHgUefco949fwY84p33AZee6M3N7FPApwDmzp07PZ9IRIpewO/jnUvqeOeSOlLpDBt2H+HRrQd5ZOsBHt/Rjc+2sHZeNdesqOea5Q3MrSktdJGlSKXGH8OiACoiIsUjlwH0RL8x3XH7twP/YGabgC3Ai0DKzKrItpa2AUeBH5nZHznn7gFuAz7vnLvfzD4A/Avwrje8kXN3AncCrF279vj3FRE5pYDfx8ULarh4QQ3/48Zz2NbVPx5Gv/rz7Xz159tZ1lDGNSsauGZ5PSuayjFTWJDpkUqPtYBqDKiIiBSPXAbQfUDrhP0Wjusu65zrBz4GYNlvbR3e8jtAh3PukHfuJ2RbOu8BPgp8znuJHwF35e4jiIhkmRkrmipY0VTB59+9hD29wzy67QCPbjvIPz3xKt94/FUaysNcuayOK5bO4e2LajVuVKYkndEYUBERKT65/Hb0W2CxmbUBnWQnEfqDiReYWSUw7JxLAJ8EnnbO9ZvZHuBiMysl2wX3amCDd9t+4HLgKeAq4NUcfgYRkROaW1PKJ9+xgE++YwG9g3Ee39HNkzu6+enmLu79zV5Cfh8XtlVz5bI5XLVsDm210UIXWWaY1PgsuAqgIiJSPHIWQJ1zKTP7NPAI4Afuds5tNbNbvfPfAs4Bvm9mabKTE33CO7fezH4MbARSZLvm3um99H8h2203AIzijfMUESmUmlgJH1jbygfWtpJIZdiw+zBP7ujmyZ2H+MrPtvGVn21jfk0pVy6bw5VL53BhWzXhoL/QxZazXFrPARURkSJkzhX/8Mi1a9e6DRs2nPpCEZFptvfwME/u7OaJHd2s29VLPJWhJODjgnlVXLyghksW1nBuSyWhgMb5yWTffGoXdzy8g+3/61oiIf3BQkREZhYze8E5t/b44xqgJCKSQ63VpXzkkvl85JL5jCTS/HpXD8+91su69l6+/tgr8BhEgn7Wzj8WSFc3VxDwK5DOdmNjQDULroiIFBMFUBGRPImE/Fx9Tj1Xn5N93PGRoQTrO3pZtysbSP/2kZ0AREN+3tZWzSVeIF3RVKFumLPQ2BhQv2ZWFhGRIqIAKiJSIFXRENeubOTalY0A9AzGWd9+mHXtPazb1ctTOw8BUBYOcFFbNRe2VXPBvGpWNpdTElCXzGKXSjt8Bj798UFERIqIAqiIyFmiNlbCDasbuWF1NpB294+yrr2X59uzraS/3N4NQCjgY3VzBRfMqxpfamIlhSy65EAq4/QMUBERKToKoCIiZ6k55WFuXtPMzWuaATg0EOeF3UfYuOcIG14/zHefe51vP90OQFttlPPnVrF2fhXntlSyuD5GUONIZ7R0JqPxnyIiUnQUQEVEZoi6shKuXdnAtSsbABhNpnm5s48Xdh9hw+4jPLWzm/s37gOyraTLG8tZ1VzBqpYKVrdUsKgupsmNZpBk2mnsr4iIFB0FUBGRGSoc9LN2fjVr51fzXwHnHLt7h3mps48t+47y0r4+Hnixkx88v9u7PhtKV7dUsrI5G0oX1sUUcs5S6YwjoP83IiJSZBRARUSKhJkxvzbK/NooN53bBEAm4+joHWLLvj62dPaxZV8f923Yy/d+/TqQfQTMiqZyljeVs7ShjGUNZSyuL6M8HCzgJxHwxoCqxVpERIqMAqiISBHz+YyFdTEW1sW45bzsWNJ0xtF+aJAtnX285AXT+1/Yx1AiPX5fc2WEpQ1lLKnPhtIl9WUsnBPV7Lt5lEpn1AIqIiJFRwFURGSW8fuMxfXZls73nN8CZLvv7jsywisHB9hxYICdBwZ45eAAz7x6iGQ6+zzKgM9oq42ypKGMJXPKWFAXZWFdjLbaKJGQgul0S2c0BlRERIqPAqiIiGBmtFaX0lpdytXn1I8fT6QyvN475IXSfnYeGOSlfUf5xZYunDt2f3NlZDyQLqyLsqAuxoK6KA3lYcwUos5EKuM0k7GIiBQdBVARETmpUMDHkvpsF1y8caUAI4k0HT1DtPcM0n5oiF2Hsusfbdg7qStvachPS1WE1qrS8YDbWhUZ346V6NfQyagFVEREipF+84uIyFsWCflZ7k1eNJFzjoP9cdoPDbKrZ4iOQ0PsPTLM3sPDrO84zGA8Nen6qtKgF0pLaak+FlSbKsI0VIQpm8WTISU1BlRERIqQAqiIiEwbM6PBC4+XLqqddM45x9HhJHsOD3uhdGQ8nG7r6uexbQdJpDOT7omVBGioCNNYEaa+PLueuF9fHqaqNFSULYVqARURkWKkACoiInlhZlRFQ1RFQ5zbWvmG85mM4+DAKHsPj9DVN8LB/lG6+kY50Jddv3qwh+6BUTJu8n0+g+poiJpoCTWxELWxY+va2OTjtbGSGTNhkh7DIiIixUgBVEREzgo+n9FYEaGxInLSa1LpDD2DCbr6RujqG+XQQJzewTg9Qwl6BuL0DiXYvO8ovYOJN3T3HRMK+KiIBCkPB7LrSNDb99aRwHH7QcrCASIhP6WhAKVBP748tEymMuqCKyIixUcBVEREZoyA3zfexfe8U1w7mkzTMxindzBB71CcnoEEPUNx+kaS9I8k6R9J0TeS5PBQgo6eoeyx0RTp45tYTyAc9BENjYXSbDCNBP0EAz5Cfh+hgBH0Z7ePHfMR9HvHvWNBvw+fz/AZ+M3wmY3vH+yPUx0NTc9/OBERkbOEAqiIiBSlcNBPS1UpLVWlp32Pc46hRHo8pPZ5y1A8xXAizXAiux5JpBnytofjaYaTaUYTaUZGkiTTGRKpzPg6kXaTjqVOI+COuWFV45l8dBERkbOWAqiIiIjHzIiVBIiVBGiuPHlX4KlIZ7KBdCyUZhxknPOW7FjYjHOkM47mqtyUQUREpFAUQEVERPLI7zP8Pj/h4MyYDElERGQ6aXo9ERERERERyQsFUBEREREREckLBVARERERERHJCwVQERERERERyQsFUBEREREREckLBVARERERERHJCwVQERERERERyQsFUBEREREREckLBVARERERERHJCwVQERERERERyQtzzhW6DDlnZoeA3YUuRxGqBXoKXQiRt0B1VmYa1VmZaVRnZaZRnc2dec65uuMPzooAKrlhZhucc2sLXQ6R06U6KzON6qzMNKqzMtOozuafuuCKiIiIiIhIXiiAioiIiIiISF4ogMpU3FnoAoi8RaqzMtOozspMozorM43qbJ5pDKiIiIiIiIjkhVpARUREREREJC8UQEVERERERCQvFEDlTZlZtZk9Zmaveuuqk1x3rZntNLPXzOxLJzj/BTNzZlab+1LLbDbVOmtmf2tmO8zsJTN7wMwq81d6mU1O4+emmdk3vPMvmdn5p3uvSC6caZ01s1Yze9LMtpvZVjP7XP5LL7PRVH7Oeuf9Zvaimf0sf6UufgqgcipfAh53zi0GHvf2JzEzP/DPwHXAcuBDZrZ8wvlW4N3AnryUWGa7qdbZx4CVzrnVwCvAX+al1DKrnOrnpuc6YLG3fAr45lu4V2RaTaXOAingz51z5wAXA3+qOiu5NsU6O+ZzwPYcF3XWUQCVU7kZ+Fdv+1+BW05wzYXAa865dudcAvh3774xfwf8BaAZryQfplRnnXOPOudS3nXPAy05Lq/MTqf6uYm3/32X9TxQaWaNp3mvyHQ74zrrnOtyzm0EcM4NkP1C35zPwsusNJWfs5hZC3ADcFc+Cz0bKIDKqdQ757oAvPWcE1zTDOydsL/PO4aZ3QR0Ouc257qgIp4p1dnjfBx4aNpLKHJ6dfBk15xu/RWZTlOps+PMbD5wHrB+2ksoMtlU6+zfk21AyeSqgLNVoNAFkMIzs18CDSc49eXTfYkTHHNmVuq9xjVnWjaRE8lVnT3uPb5MttvYD99a6UROyynr4Jtcczr3iky3qdTZ7EmzGHA/8GfOuf5pLJvIiZxxnTWzG4Fu59wLZnbFtJdsllMAFZxz7zrZOTM7ONZ9xuuS0H2Cy/YBrRP2W4D9wEKgDdhsZmPHN5rZhc65A9P2AWTWyWGdHXuNjwI3Alc7PSxZcuNN6+Aprgmdxr0i020qdRYzC5INnz90zv0kh+UUGTOVOvs+4CYzux4IA+Vmdo9z7o9yWN5ZQ11w5VQeBD7qbX8U+M8TXPNbYLGZtZlZCPgg8KBzbotzbo5zbr5zbj7Zf+TnK3xKjp1xnYXsjHnAF4GbnHPDeSivzE4nrYMTPAh8xJul8WKgz+tWfjr3iky3M66zlv0r9L8A251zX89vsWUWO+M665z7S+dci/f99YPAEwqf00ctoHIqtwP3mdknyM5i+34AM2sC7nLOXe+cS5nZp4FHAD9wt3Nua8FKLLPdVOvsPwElwGNey/3zzrlb8/0hpLidrA6a2a3e+W8BvwCuB14DhoGPvdm9BfgYMotMpc4ClwEfBraY2Sbv2F85536Rz88gs8sU66zkkKl3mYiIiIiIiOSDuuCKiIiIiIhIXiiAioiIiIiISF4ogIqIiIiIiEheKICKiIiIiIhIXiiAioiIiIiISF4ogIqIiIiIiEheKICKiIhMAzOrMbNN3nLAzDq97UEz+385eL/vmVnH2DPtTnB+0FsvHCvHdJdBRETkrdJzQEVERKaZmf01MOic+1oO3+N7wM+ccz8+yflB51zsZPsiIiKFECh0AURERIqZmV0BfME5d6MXTNuARmAJ8N+Ai4HrgE7gd51zSTO7APg6EAN6gD92znWd4n3agH8j+7v94dx8GhERkalRF1wREZH8WgjcANwM3AM86ZxbBYwAN5hZEPhH4H3OuQuAu4H/fRqv+w/AN51zbwMO5KTkIiIiU6QWUBERkfx6yGvl3AL4OdZauQWYDywFVgKPmRneNW/a+um5DHivt/0D4I5pLLOIiMi0UAAVERHJrziAcy5jZkl3bDKGDNnfywZsdc5dcgavrYkdRETkrKYuuCIiImeXnUCdmV0CYGZBM1txGvc9B3zQ2/7DXBVORERkKhRARUREziLOuQTwPuAOM9sMbAIuPY1bPwf8qZn9FqjIYRFFRETOmB7DIiIiMgOd6jEsJ7hej2EREZGCUwuoiIjIzNQHfMXMbn2zi8xsoZltAg7mp1giIiInpxZQERERERERyQu1gIqIiIiIiEheKICKiIiIiIhIXiiAioiIiIiISF4ogIqIiIiIiEhe/H+dI/6nVaZpfgAAAABJRU5ErkJggg==\n",
"text/plain": [
""
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"plot_lc(time, flux1);"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"### Thawing the model\n",
"\n",
"After thawing, the model takes a random sample from the limb darkening profile sample set every time it is evaluated. We don't want this behaviour when fitting a model to observations, since this sort of randomness can easily confuse even the best optimiser, but is exactly what we want then doing MCMC for parameter posterior estimation."
]
},
{
"cell_type": "code",
"execution_count": 70,
"metadata": {},
"outputs": [],
"source": [
"ldm.frozen = False"
]
},
{
"cell_type": "code",
"execution_count": 78,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAA6AAAAGoCAYAAABYCYj1AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nOzdeZjdd2Hf+/d3ZrRakrWNFmuxvMiLLNuyLWTAmDWAjQEvLDGkxSULcQJc6G16Q5LbJL1tnyZp2iwNDSGBpLnNhdALJJC6gCEE18Y2trFlSxrt68xo9uXs6+/bP84xkYVkSZbO/GZ5v55Hz8z5LWc+55Hh0Wd+3yXEGJEkSZIkqdXa0g4gSZIkSZoZLKCSJEmSpAlhAZUkSZIkTQgLqCRJkiRpQlhAJUmSJEkToiPtABNh+fLlccOGDWnHkCRJkqQZ4ZlnnhmKMXaefHxGFNANGzbw9NNPpx1DkiRJkmaEEMKRUx13CK4kSZIkaUJYQCVJkiRJE8ICKkmSJEmaEBZQSZIkSdKEsIBKkiRJkiaEBVSSJEmSNCEsoJIkSZKkCdGyAhpC+HwIYSCEsOM050MI4Q9DCPtDCM+HEG4+4dwdIYQ9zXOfOuH4fwgh7G5e/9UQwuJW5ZckSZIkXVitfAL6F8AdL3P+TmBj889HgD8GCCG0A59unt8EfCCEsKl5z8PA5hjjDcBe4FdaklySJEmSdMG1rIDGGB8BRl7mkruBv4wNTwCLQwirgW3A/hjjwRhjBfhi81pijN+KMdaa9z8BrG1VfkmSJEnShZXmHNA1wLETXnc3j53u+Ml+Gvifp3vzEMJHQghPhxCeHhwcvABxJUmSJEnnI80CGk5xLL7M8X+8MYRfA2rAX53uzWOMn40xbo0xbu3s7DyvoJIkSZKk89eR4s/uBtad8Hot0AvMPs1xAEIIDwDvBN4SY3xJMZUkSZIkTV5pPgH9GvCh5mq4rwbGY4zHgaeAjSGEy0IIs4H7m9cSQrgD+GXg3THGQlrBJUmSJEnnrmVPQEMIXwDeCCwPIXQDvwHMAogxfgZ4CHgHsB8oAB9unquFED4GfBNoBz4fY9zZfNs/AuYAD4cQAJ6IMT7Yqs8gSZIkSbpwWlZAY4wfOMP5CHz0NOceolFQTz5+5YVJJ0mSJEmaaGkOwZUkSZIkzSAWUEmSJEnShLCASpIkSZImhAVUkiRJkjQhLKCSJEmSpAlhAZUkSZIkTQgLqCRJkiRpQlhAJUmSJEkTwgIqSZIkSZoQFlBJkiRJ0oSwgEqSJEmSJoQFVJIkSZI0ISygkiRJkqQJYQGVJEmSJE0IC6gkSZIkaUJYQCVJkiRJE8ICKkmSJEmaEBZQSZIkSdKEsIBKkiRJkiaEBVSSJEmSNCEsoJIkSZKkCWEBlSRJkiRNCAuoJEmSJGlCWEAlSZIkSRPCAipJkiRJmhAWUEmSJEnShLCASpIkSZImhAVUkiRJkjQhLKCSJEmSpAlhAZUkSZIkTQgLqCRJkiRpQlhAJUmSJEkTwgIqSZIkSZoQFlBJkiRJ0oSwgEqSJEmSJoQFVJIkSZI0ISygkiRJkqQJYQGVJEmSJE2IlhXQEMLnQwgDIYQdpzkfQgh/GELYH0J4PoRw8wnn7ggh7Gme+9QJx5eGEB4OIexrfl3SqvySJEmSpAurlU9A/wK442XO3wlsbP75CPDHACGEduDTzfObgA+EEDY17/kU8J0Y40bgO83XkiRJkqQpoKNVbxxjfCSEsOFlLrkb+MsYYwSeCCEsDiGsBjYA+2OMBwFCCF9sXrur+fWNzfv/K/APwC+fKUu+VOXxvcdfcqyjrZ22We2EcMLB8NL7Tnp56mPNN3jxeAftL7kmhPDSn/Gj93npwROv6Tjpb2VWx8vlCSf9vB+/ZtZJr5ddfBHhVBdKkiRJOmdJkvB8zyjFYiXtKJNeywroWVgDHDvhdXfz2KmO39r8fmWM8ThAjPF4CGHF2fygg8MFPvD5H55/4mlk2fxZ/Ku7ruGeW9anHUWSJEmakp7vHuZn/vwZBvPVtKNMGWkW0FM9gosvc/zc3jyEj9AY2su8VZezsXP+y15/qh9wyh8aIZ585hQXnk3gl1wT4znnOPOh+OLbvuSi4XyF4UKVT/73F/jM/zrAV3/hdcybc/JzUkmSJEmn88UnD/IrX+0iAhfNbmPFwjl0tDnK8EVHTnM8zQLaDaw74fVaoBeYfZrjAP0hhNXNp5+rgYHTvXmM8bPAZwG2bt0aH/4Xb7qQ2ae8/+/xQ/zG13exu6/Arf/+7/m7j97K+s7FaceSJEmSJr0vP3X4R+Xzni2r+P37b0k70qQTfunUx9PchuVrwIeaq+G+GhhvDq99CtgYQrgshDAbuL957Yv3PND8/gHgbyc69HTxwddcxgu/+TZWLpxFplTjrj96gt6RXNqxJEmSpElt9/FRfvkrO4nAR27fYPk8R63chuULwOPA1SGE7hDCz4QQHgwhPNi85CHgILAf+FPgFwFijDXgY8A3gS7gSzHGnc17fgt4awhhH/DW5mu9QnNnz+KJX30raxbPJluuc89/eZysE6clSZKkU6rXEz74Z09Ri/CWq5bxq3ddl3akKSfEH5skOP1s3bo1Pv3002nHmLSq9YSt/+Zhxks1brhkIV/56OvoaE/z4bgkSZI0+fwfX3iar23vZ8m8Dp79jbenHWdSCyE8E2PcevJxW4aY1d7Gt/757bQFeL43y7/5+i5mwi8mJEmSpLPV1TvC323vB+Cvf/7VKaeZuiygAmDlxfP5vfddD8BfPXmEb7/Qe4Y7JEmSpJnjwf/2LAnwzutWcNWqi9OOM2VZQPUjd9+8njdds5xahF/+m50MZYtpR5IkSZJS9/VnuzkyUmJWG/zhT/3YqFKdAwuoXuJzH9rGojntjBSqfOrLL1BPHIorSZKkmateT/jNv9sFwC++8XLa3OvzvFhA9RJtbYE/faCxlPQ/7Bnk4Z0OxZUkSdLM9aePHGA4X+WiWW188q3XpB1nyrOA6sfcenknb9jYGIr7b/9HF5liNe1IkiRJ0oSr1hL++B8OAPCv795ECD79PF8WUJ3SH/+Tm+log+6xMn/y93vSjiNJkiRNuM99bx/j5ToXz23nvVsvTTvOtGAB1SnNnzOLj77xcgD+21PdHBrKpZxIkiRJmjjVWsKfPHoYgE/dcXW6YaYRC6hO65NvvYZFc9oZL9X5/W/vIXFBIkmSJM0QX37mCKPFGhfNbuMDr74s7TjThgVUpxVC4DffvQmA7+waYE9fNuVEkiRJUuvV6wm//53G3M9ffMMVKaeZXiygeln33ryOpfM6yFUS/vC7u92WRZIkSdPeP+zrpz9TpqMNfuFNG9OOM61YQPWyQgj8+juvBeCRPcNsPzaSciJJkiSptf7Tt/YRgXdct8J9Py8wC6jO6O7mU9B8JeGzjxzwKagkSZKmrcOD2R9NPfvNe65POc30YwHVGYUQ+LV3NDbdfXT/MDt6x1JOJEmSJLXG7317L7UErl5xEUsvmpt2nGnHAqqzct/W9Syc006unPBfH/UpqCRJkqafbLHK33cNAvBv7rku5TTTkwVUZyWEwM+/obEv6Hf3DLP7+HjKiSRJkqQL62+fO0a2UmfBnDa2Xd6ZdpxpyQKqs/bzb7iS2e0wXqzxhSePuC+oJEmSpo0YI3/x/aMAvOeWtSmnmb4soDprs9rbeP/WtSTAw10DHBx0X1BJkiRND/v7sxwezgPwS2+9OuU005cFVOfkX779GgIwlKvwtWd70o4jSZIkXRB//ugBaglcvmw+C+fNTjvOtGUB1Tm5eP4c3nx1J/UIf7ejj4FMKe1IkiRJ0nkpVWp8c+cAAJ94y+Upp5neLKA6Z7/yzmsB6Bkr8N3d/SmnkSRJks7PY/sHGS3WmN0G775pfdpxpjULqM7ZlZ0LuWzpPMo1+Oqz3WRL1bQjSZIkSa/Y5x87TALcetlSQghpx5nWLKB6RT76lisB2NWb4YeHR1JOI0mSJL0yA5ki24+NAfBLd7j4UKtZQPWK3HfTOuZ2BHLlhL95tpu6W7JIkiRpCvrGzuPkKwkLZrdx47qlaceZ9iygekXa2gL33byGBHjy0CiHh3JpR5IkSZLOSZIk/PUPuonA2zatTDvOjGAB1Sv2iZ+4CoChfJnvdrkYkSRJkqaWo8MFDg42HqR88q1XpZxmZrCA6hVbuWgeV3bOp1KHb3X1MVaopB1JkiRJOmt/91w3xVpk2fwO1i9bkHacGcECqvPykdsb+yTt6cvx3LHRlNNIkiRJZ6deT/jK9l4A7rlpTcppZg4LqM7LvbesY1Yb5Cp1vvFCn4sRSZIkaUrYN5Cje6wEwMff4vDbiWIB1XmZ1d7GG65aTj2Bp48Mc2Qon3YkSZIk6Yz+9ofdVGqRlQtnsXj+7LTjzBgWUJ23j7258RujnrESj+0fSjmNJEmS9PIqtYRv7OoD4J6b1qacZmaxgOq8bVm/hGUXdVCqRr7TdZx8uZZ2JEmSJOm0dvWM0TveGH77kTdckXKamcUCqgvi/VvXEYFdx7Ps7sukHUeSJEk6rYde6KFci6y4qINlF81JO86MYgHVBfEzr2ushjtWrPLdrgFidDEiSZIkTT61WsK3dzemjd114yUpp5l5LKC6IJYvnPujPUEf3z/IYK6cdiRJkiTpx+zpz3B8rAjAg2/cmHKamccCqgvmg9suBeDQcJ7njronqCRJkiafh3cep1iLLJ3XzspFc9OOM+NYQHXB3L9tPW1Atlzn73f3U6snaUeSJEmSfqReT/jGzgEA3nrdqpTTzEwWUF0w8+d0sGXdIqoJPHtsjMPDubQjSZIkST9yaDjP0ZECAL/4Rle/TYMFVBfUh29rLEZ0fLTIM4cchitJkqTJ45E9/RSrCQtmt3Hp8oVpx5mRWlpAQwh3hBD2hBD2hxA+dYrzS0IIXw0hPB9C+EEIYfMJ5z4RQtgRQtgZQvjkCce3hBCeCCE8F0J4OoSwrZWfQefmzutXM6sNCtWER/YNUay4J6gkSZLSlyQJD73QRwReffnStOPMWC0roCGEduDTwJ3AJuADIYRNJ132q8BzMcYbgA8Bf9C8dzPwc8A24EbgnSGEF5eo+h3gX8cYtwC/3nytSaKjvY3bNy6nlkBX3ziHhvNpR5IkSZLoHi2yr78xReyf3XZZymlmrlY+Ad0G7I8xHowxVoAvAnefdM0m4DsAMcbdwIYQwkrgWuCJGGMhxlgDvgfc27wnAoua318M9LbwM+gVeHFP0IFMmacPjaScRpIkSYKnDg+Tq9TpCHDblZ1px5mxWllA1wDHTnjd3Tx2ou3AfQDNobSXAmuBHcDrQwjLQgjzgXcA65r3fBL4DyGEY8DvAr9yqh8eQvhIc4ju04ODgxfoI+lsvPbK5czrgGIt4fsHhsmXHYYrSZKkdH3jhePUI2xavZAQQtpxZqxWFtBT/a3Gk17/FrAkhPAc8HHgWaAWY+wCfht4GPgGjaL6Yov5BeCfxxjXAf8c+NypfniM8bMxxq0xxq2dnf6GYyKFELhtYyf1BPYOjLNvwNVwJUmSlJ6RXJntx8YBeN/WdWe4Wq3UygLazT8+tYTGk82XDJeNMWZijB9uzuf8ENAJHGqe+1yM8eYY4+uBEWBf87YHgK80v//vNIb6apJ54DWNcfWD4xWeOjSUchpJkiTNZM8eG2W0WAXgPbesTTnNzNbKAvoUsDGEcFkIYTZwP/C1Ey8IISxungP4WeCRGGOmeW5F8+t6GsN0v9C8rhd4Q/P7N/OPxVSTyG1XLmducxjuk4dGyTT/By9JkiRNtO/u6qeawJrFc5g/Z1bacWa0jla9cYyxFkL4GPBNoB34fIxxZwjhweb5z9BYbOgvQwh1YBfwMye8xZdDCMuAKvDRGOOLm0r+HPAHIYQOoAR8pFWfQa9cW1vgdRs7+XbXIPsHMuztz7B1w7K0Y0mSJGmGKVZqPHZwGIA7Nq1KOY1aVkABYowPAQ+ddOwzJ3z/OLDx5Pua524/zfFHgVsuYEy1yD999Qa+3TXIcLbCM4dHLaCSJEmacLv7MvRnygA8cNuGdMOopUNwNcPdvrGTOe2NYbhPHx0l6zBcSZIkTbDv7xukWE1YOLud9csWpB1nxrOAqmXa2gK3XbmcWgL7B7LsG8imHUmSJEkzSL2e8O2uASKw7fKlaccRFlC12AdvvRSA4WyZZ4+OpJxGkiRJM8mh4TwHB/MAPPDaDemGEWABVYu96ZqVzGqDQjXhmUNj5EoOw5UkSdLE+OGRUfKVOu2hMT1M6bOAqqXa2wJbNyyhlsDBoSwHmr+BkiRJklrte1191CJsXHERIYS04wgLqCbAe7c2Nvvtz5R57tjoGa6WJEmSzt9IrswPezIA3HXD6pTT6EUWULXcXdevoY3GMNynD4+QLzsMV5IkSa21q3ecsUIFgA801yVR+iygarm5s9q5dtUCyvXIgYEsh4YKaUeSJEnSNPfEwWFK1cjiee0sXzA37ThqsoBqQrxryxoABjJlnncYriRJklqoVk94dN9gY/uVDW6/MplYQDUh3ndLYx5orlLnmcOjlKr1lBNJkiRpujo6UuDwcGPU3f3bHH47mVhANSGWLZzLusVzqdQiBwazHBrKpR1JkiRJ09Qzh0fIleu0AW+4ekXacXQCC6gmzJ3XryYB+jJldvSMpx1HkiRJ09Sj+waoRdiwbD7tbW6/MplYQDVh7t+2HoBMqcYzh0ep1hyGK0mSpAtrvFDhh0cbDzvuvH5Vyml0MguoJszlnQtYNr+DSj3hwGCOoyPFtCNJkiRpmuk6nmU4Xwbgg7duSDeMfowFVBPqTdesoJZA73iRvf0Ow5UkSdKF9dShQUrVyMLZ7axZMi/tODqJBVQT6idf1ViFbLxQ44dHxkiSmHIiSZIkTRf1esJjB4ZJgBvXX5x2HJ2CBVQTauuGJczrCJRrdfb0Z+nPlNKOJEmSpGmiZ6zI/v7Gbgv33rQm5TQ6FQuoJlQIgVdfsZxqAt0jefb2Z9OOJEmSpGlie/cY4+UaAHduXp1yGp2KBVQT7r6bG7+NGslXea57NOU0kiRJmi6eODhMtQ4rF85m/pxZacfRKVhANeHedt0q2oBiNWFXzzijzVXKJEmSpFcqX67y7KERAF5z+bKU0+h0LKCacHM62rlm1QIq9cjRkQL7BnNpR5IkSdIUt7cvS09zfZH3bl2bchqdjgVUqXjbdauIwHCuQldPJu04kiRJmuK2d4+TL9dpA157RWfacXQaFlCl4j23rAMgX66z/dgY+eZkcUmSJOlcxRh5+uAQtQgbls+nrS2kHUmnYQFVKtYtnc/y+R2UawlHRwocGc6nHUmSJElT1HCuzI6+xqi6N17t08/JzAKq1Lzuqk5qEfqyRfb0uR2LJEmSXpk9/VmGshUAPvCq9Smn0cuxgCo177m5MTk8U6jx3LEx6vUk5USSJEmaip4/NkqxmjC3I7Bx1aK04+hlWECVmtde2cmsNijV6hwcyNIzVkw7kiRJkqaYJIk8dWiEeoSrVi5MO47OwAKq1LS3BbasW0ylDn2ZIvvdjkWSJEnnqG+8SFdf49+Rd16/KuU0OhMLqFL1jutXAzCSq7GjeyzlNJIkSZpq9vbnGCs25n++t7nTgiYvC6hSdfeWSwAo1Ors6csymi+nnEiSJElTybNHRyhXIwvntNO5cG7acXQGFlClaumCuaxZPJdKLeHYSIFDQ4W0I0mSJGmKqNUTnjw0RALcsPbitOPoLFhAlbo3Xd1JPcJwoUzX8fG040iSJGmK6Bkrsn+g8QDjvpvXpJxGZ8MCqtTdd1NjrH6mUGdnb4ZytZ5yIkmSJE0Fu/syZMo1AO7cvDrlNDobFlCl7qZLFzO3A8r1OkeG8xwdyacdSZIkSVPAM4dGqdQinRfNYv6cWWnH0VmwgCp1IQRuWr+USh0GMkX2DbgdiyRJkl5euVrniUPDALzq8qUpp9HZsoBqUrhj80oARvM1dvSME2NMOZEkSZIms2OjhR+NnHv/zWtTTqOzZQHVpPCuGxqTxgvVOocGc/RnSiknkiRJ0mTW1TtOrtJYO+R1V61IOY3OlgVUk8LSBXO45OI5VOoJPaNFDg45D1SSJEmn94ODw9TqcMnFc+lot9ZMFS39mwoh3BFC2BNC2B9C+NQpzi8JIXw1hPB8COEHIYTNJ5z7RAhhRwhhZwjhkyfd9/Hm++4MIfxOKz+DJs7tGzupJTBarLCnN5N2HEmSJE1ShUqNp4+OAvCay5eknEbnomUFNITQDnwauBPYBHwghLDppMt+FXguxngD8CHgD5r3bgZ+DtgG3Ai8M4SwsXnuTcDdwA0xxuuA323VZ9DEuvemxjDcbLFOV1+GXKmWciJJkiRNRkeGChwfa0zZuu8W539OJa18AroN2B9jPBhjrABfpFEcT7QJ+A5AjHE3sCGEsBK4FngixliIMdaA7wH3Nu/5BeC3Yozl5n0DLfwMmkDbLlvG7DYo1Wr0jBQ5MuJquJIkSfpxO5vzPwPw6ss7046jc9DKAroGOHbC6+7msRNtB+4DCCFsAy4F1gI7gNeHEJaFEOYD7wDWNe+5Crg9hPBkCOF7IYRXtfAzaAK1tQVuXLeEcg0G80X29ltAJUmS9ON+cHCYWgKXLJ5Le1tIO47OQSsL6Kn+Szh5b43fApaEEJ4DPg48C9RijF3AbwMPA9+gUVRfHI/ZASwBXg38S+BLIYQf+1khhI+EEJ4OITw9ODh4IT6PJsDbrltJBMbyNXYdH6deT9KOJEmSpEkkW6zyXPcYAK+7clnKaXSuWllAu/nHp5bQeLLZe+IFMcZMjPHDMcYtNOaAdgKHmuc+F2O8Ocb4emAE2HfC+34lNvwASIDlJ//wGONnY4xbY4xbOzt9LD9V3L3lEgCK1YSjgwWOux2LJEmSTnBgKE//eOPfiO9x/88pp5UF9ClgYwjhshDCbOB+4GsnXhBCWNw8B/CzwCMxxkzz3Irm1/U0hul+oXnd3wBvbp67CpgNDLXwc2gCrVg0j1WLZlOq1enLljg44DBcSZIk/aOdPaPkK3XagK0bfAI61bSsgDYXD/oY8E2gC/hSjHFnCOHBEMKDzcuuBXaGEHbTWC33Eye8xZdDCLuArwMfjTGONo9/Hrg8hLCDxsJGD8QYTx7aqynstVcsp5rASL7C7j63Y5EkSVJDjJEnD45Qi7B+6TzanP855XS08s1jjA8BD5107DMnfP84sPE0995+muMV4J9cwJiaZO7esoavPNtLtlxjb3+OXKnKgrmz0o4lSZKklI3kK+zsGQfgto0+/ZyKWjkEV3pFbrtyOR1tUKok9I4WOTySTzuSJEmSJoGDg1kGchUA3nfzujNcrcnIAqpJp6O9jc1rLqZcSxjKl9jvdiySJEkCdvSMU6zUaQ9w4/olacfRK2AB1aT0tk2N7VjGC1W6esdJErdjkSRJmsmSJPKDQ435nxuWzucUOzFqCrCAalK6e8saAArVhENDBY6Pux2LJEnSTDaUK9HVlwXg9Ve7zeJUZQHVpLRmyXyWL5hNqVqnP1vi4KDzQCVJkmayff05hrNlAN5z85qU0+iVsoBq0rr18qVUExjLV9jrdiySJEkz2vM9oxRrCR1tcN2axWnH0StkAdWk9e4bLgEgU66zpz9LrlxNOZEkSZLSUE8iPzw8Si2By5df5PzPKcwCqknrTdespA0oVer0jBY5OlJIO5IkSZJS0DdWZG9zZ4Q3Of9zSrOAatKa3dHG1asWUq4nDOVK7O93HqgkSdJMtKc/y3C+sf/ne29Zn3IanQ8LqCa1N1+7giRCtlSn6/gYMca0I0mSJGmCvdA9RrGaMLs9cOXKBWnH0XmwgGpSe3EeaK5S59BQnv6M27FIkiTNJJVanWePjVCPcEWn8z+nOguoJrWrVy9i4Zx2ytUa/Zmy27FIkiTNMD2jRQ4MNNYCefM1zv+c6iygmvS2XrqUch3GC2X29GfTjiNJkqQJtHcgy0ihMf/zvpvXpZxG58sCqknvzutXAZApJuzuy1Cs1FNOJEmSpIny/LExytWEOR2BK1YsTDuOzpMFVJPeO65fDUCpVqNnrMjREYfhSpIkzQSlao3nj41Si7BxhYsPTQcWUE16C+bO4rJlF1GqJQxnyuwbyKUdSZIkSRPg2HCBA4ON+Z8/cZ3zP6cDC6imhNdduYxaAtlKja7ecbdjkSRJmgG6+rKMFxvzP+/d4v6f04EFVFPC3VvWAJAv1Tg8UmAwW045kSRJklrthe4xSrXIvI42Ll12UdpxdAFYQDUl3LJhCXM7AsVqnb6xIgeHHIYrSZI0nRXKNXb0jFGPcNVK539OFxZQTQkhBG5Yu5hyLTJeqLCnzwIqSZI0nR0aznN4uDH/8+3Xr0w5jS4UC6imjLddt5II5Ep1dh/PUHI7FkmSpGmrqzfDeLEKwN1b1qacRheKBVRTxrtvbM4DrdQ5Plbg2KjbsUiSJE1HMUaePzpCuRaZP7uNNYvnpx1JF8hZFdAQwqZTHHvjBU8jvYwVi+ayatFsSrU6w/kK+wcsoJIkSdNRtlSjqz9LPcK1qxamHUcX0Nk+Af1SCOGXQ8O8EMJ/Bv59K4NJp/Kay5dTTSBbqtJ13O1YJEmSpqMDAzmOjjTmf96xeVXKaXQhnW0BvRVYB3wfeAroBW5rVSjpdN615RIAsqU6B4fyDOcqKSeSJEnShba7b5xsqQbAu5rTsDQ9nG0BrQJFYB4wFzgUY0xalko6jddv7KSjDYrVOv1jJQ4OuhquJEnSdBJj5PnuMcq1yMI57ay6eF7akXQBnW0BfYpGAX0V8DrgAyGE/79lqaTT6Ghv49rViyjVEjKlCnv7s2lHkiRJ0gU0ki+z53hj/uem1YvSjqML7GwL6M/EGH89xliNMfbFGO8G/raVwaTTefM1K0kiZMtuxyJJkjTd7O/P0T1WBODtzv+cds62gA6EENaf+Af4XiuDSadzz02NeaCFcp3usSLdo4WUE0mSJOlC2d2XJVtuzP+8u7n+h6aPjrO87n8AEQg05oBeBuwBrmtRLum0Llu+gCXzOihUaozkK+wfyHHlSpfnliRJmurqSeSFnjEqtciiOe0sWzA37Ui6wLidOOUAACAASURBVM7qCWiM8foY4w3NrxuBbcCjrY0mnd7WDcso1yFXrLLL7VgkSZKmhcFsiX0Djfmfm9c6/3M6OtshuC8RY/whjQWJpFTcdUNjPkC2UuPwcIGhXDnlRJIkSTpf+/qz9I425n/eed3qlNOoFc5qCG4I4f884WUbcDMw2JJE0ll4+3WrCGynVEnoyxQ5NJSnc6FDNCRJkqay3ccz5JoLTN51o/M/p6OzfQK68IQ/c2jMCb27VaGkM5k3u4MrVlxEsVYnV6yxt8/tWCRJkqayWj1hZ2+GSi2yeF47Sy+ak3YktcBZPQGNMf7rVgeRztXtVyxn/0CebLnGnv4sxUqNebPPdl0tSZIkTSZ94yUODDbmf96wbnHacdQiL/uv9RDC12msfntKMcZ3X/BE0lm675a1/PnjR8iX6vSMFukeKbJxlavhSpIkTUV7+7McHysB8A7nf05bZ3pc9LsTkkJ6BTavuZiLZrVRrNYZLVTYP5i1gEqSJE1RXb3j5KuN+Z93bF6Vchq1ypkK6KEY49EJSSKdoxACN126hEf3D5MvVtl9PMPbr1tNW1tIO5okSZLOQaWW0NXXmP+5ZF4Hi53/OW2daRGiv3nxmxDCl1ucRTpnb920EoBsucaR4SKDbsciSZI05XSPFTg0lKce4cZ1S9KOoxY6UwE98VHS5ef65iGEO0IIe0II+0MInzrF+SUhhK+GEJ4PIfwghLD5hHOfCCHsCCHsDCF88hT3/lIIIYYQlp9rLk0f77yhsTx3oVpnIFPk8FA+5USSJEk6V3v7svSPN+d/3uDw2+nsTAU0nub7MwohtAOfBu4ENgEfCCFsOumyXwWeizHeAHwI+IPmvZuBnwO2ATcC7wwhbDzhvdcBbwUcHjzDLVswh3WL51Gs1smW6+xxOxZJkqQpZ1fvOPnm/p93XmcBnc7OVEBvDCFkQghZ4Ibm95kQQjaEkDnDvduA/THGgzHGCvBFfnzv0E3AdwBijLuBDSGElcC1wBMxxkKMsQZ8D7j3hPt+D/i/OMdSrOnpNVcspVKHXLnK/oEsxXIt7UiSJEk6S8VKjb39GSr1yLL5HSycNzvtSGqhly2gMcb2GOOiGOPCGGNH8/sXXy86w3uvAY6d8Lq7eexE24H7AEII24BLgbXADuD1IYRlIYT5wDuAdc3r3g30xBi3v9wPDyF8JITwdAjh6cHBwTNE1VT2rhsb/1nlyzV6x0scHSmknEiSJEln6+hIkSPDBeoRblrv/M/p7kxPQM/HqZYiPfmJ5W8BS0IIzwEfB54FajHGLuC3gYeBb9AoqrVmGf014NfP9MNjjJ+NMW6NMW7t7Ow8j4+hye61Vy5nVhsUKgmjuTL7B3NpR5IkSdJZ2n18nIHxxkKSdzXX99D01coC2k3zqWXTWqD3xAtijJkY44djjFtozAHtBA41z30uxnhzjPH1wAiwD7gCuAzYHkI43HzPH4YQHCg+g7W3Ba5fczGFSp18ucbevhxJ4uhsSZKkyS7GyK6+zI/2/3xbc4cDTV+tLKBPARtDCJeFEGYD9wNfO/GCEMLi5jmAnwUeiTFmmudWNL+upzFM9wsxxhdijCtijBtijBtolNybY4x9LfwcmgLedO1KEiBfrnNkJM9grpR2JEmSJJ1BrlTjQH+OSi3SuWAWF82dlXYktVjLCmhz8aCPAd8EuoAvxRh3hhAeDCE82LzsWmBnCGE3jdVyP3HCW3w5hLAL+Drw0RjjaKuyauq7d0tjuEauWmMwW+LAgNuxSJIkTXYHhnIcG8mTALc4/3NG6Gjlm8cYHwIeOunYZ074/nFg48n3Nc/dfhbvv+E8I2qaWLv0IjoXzCZbqjaH4WZ57ZVuEStJkjSZ7TmeZTBbAeDOG1annEYToZVDcKUJtXXDEkq1SL5S4+BQjnzJ7VgkSZImqySJ7O7LUPzR/E+XdZkJLKCaNu5q/tYsV6rTO17i2JjbsUiSJE1Wo/kyB4dylGuRVRfPYd7slg7O1CRhAdW08dZrV9FGYz/QTL7CnuOZtCNJkiTpNA4M5ekZLZIAr1rn/M+ZwgKqaWPOrHauXrWAQrVOoVpn30Ceaj1JO5YkSZJOoas3w3C2sf/nO5z/OWNYQDWtvPGaFdSSxlPQYyN5+sbdjkWSJGmyqdWTxvzPWkIA3nyt+3/OFBZQTSv33rQGgEy5xnCuwp4+h+FKkiRNNgPZMkeG8lRqkTWL5zJnVnvakTRBLKCaVq5auYjF8zooVhNylSp7+7MkSUw7liRJkk6wtz9L73hj/udtVyxNO44mkAVU0862DUsoVOoUy3WODhcYzpfTjiRJkqQT7OnLMJqvAnDvLWtTTqOJZAHVtHPPTY3/E8uUaxwfL7FvIJdyIkmSJL2oUkvY25+lWKvT0QbbNixPO5ImkAVU085bNq2kncZCRLlylb19GWJ0GK4kSdJk0DNW4MhwjkodLu9cQFtbSDuSJpAFVNPOnI52rl+7iHy5TqlS49BQnmyplnYsSZIkAXv7svSNNaZIvemqzpTTaKJZQDUt3Xn9auoRcpWE42MlDg/l044kSZIkYHdfhvFic/7nzWtSTqOJZgHVtHRvcx5otlhlvFily+1YJEmSUles1NjXn6NUS5jb0cbVqxalHUkTzAKqaWnForlcunQeuUqdUrXO/sEcpWo97ViSJEkz2pGhPMdG8lQTuGbVQkJw/udMYwHVtPWma1ZSqUfylTo9IwW6RwppR5IkSZrRuvqy9I+XAHj75lUpp1EaLKCatt7b3FMqU6yQKVbZ3Z9NOZEkSdLMlSSRruPjZMuNUWl333RJyomUBguopq3rLlnE4nkd5CoJuXKNvX1ZqvUk7ViSJEkz0lihwqHBLOV6wqK57Vxy8fy0IykFFlBNWyEEbrtiGYVKnWIl4ehIgYFMKe1YkiRJM9K+wRzdYyVqCdywdnHacZQSC6imtftuWQdAtlRhJFdmj8NwJUmSUrGrZ5zBbGP/z3dcvzLlNEqLBVTT2uuv6mRWG+TKdXLlGnv6MiRJTDuWJEnSjFKpJeztz1GoNOZ/3nW9+3/OVBZQTWuz2tu45dLF5Cp1KrU6R4cLjBYqaceSJEmaUXrHCxwdzVOuRVYsnM3F82enHUkpsYBq2nvXjWtIImRKVQYyZfb0OQxXkiRpInX1Zjk+WqQeYdulS9KOoxRZQDXt3XVDY4nvTKlOplRlV+84dYfhSpIkTYgYI3v6swznG6PQ7mtulaeZyQKqaW/x/NlcvXIB+XKNSrXGwcE8w7ly2rEkSZJmhHy5xv7+DKVqQnuAN1ztAkQzmQVUM8K7bryEatJYjGgwW6Lr+HjakSRJkmaEg0N5ekYLlOuRjSsW0N4W0o6kFFlANSO8b2tjO5axUpVctcau3izVepJyKkmSpOlvV+84/ZnG6LOf2OTTz5nOAqoZYeWiuWxYNp9sqU61mnB4OMfx0ULasSRJkqa1Wj1hd1+W8VINgJ/cti7lREqbBVQzxh3XraJSj+TKdcYKNZ49Nk6MLkYkSZLUKn3jRQ4P5ajUEi6e18G6JRelHUkps4BqxvjgresBGC1WKFZr7O7LMFaoppxKkiRp+tp1PMvx8SLVBG5avzjtOJoELKCaMdYvu4g1i+cwXqxRryf0ZUps7x7zKagkSVILxBjZ2TvOYK65/cqWNSkn0mRgAdWMctf1l1CqJmRLdUqVGtuPjf1oUrwkSZIunPFChYMDOYrlOgF42+ZVaUfSJGAB1Yzy4ds2ANCfKxECDGbKPHNkxBVxJUmSLrA9fTl6x4uUa5ENS+czd1ZH2pE0CVhANaOsXjyfzZcsZDhXoVZNqCUJz/eMcWQ4n3Y0SZKkaeX5njEGMyUS4I3XdKYdR5OEBVQzzk/deim1BI6Nl2lrC+SKNZ48NEKhUks7miRJ0rRQrNTZ159lrNhY8PF9t7j9ihosoJpx7rt5LfM62ujLFKnVEpIE9vZl2defTTuaJEnStHB4KE/PaJFiLWHh3HY2rbk47UiaJCygmnHmzGrn3pvWMFqocTxXJhIpVmr84NAoY4VK2vEkSZKmvBd6xhjIlqjW4dWXLU07jiYRC6hmpE++dSPtbbCvL0stSQghcngox44et2WRJEk6H9V6QtfxcYbyjV/s/+Sr1qecSJOJBVQz0opF87h/63r6MmUODeVJEihW6zxzZIzBrNuySJIkvVL940WOjhQpVOrMboM3XbMy7UiaRFpaQEMId4QQ9oQQ9ocQPnWK80tCCF8NITwfQvhBCGHzCec+EULYEULYGUL45AnH/0MIYXfznq+GEBa38jNo+vpX79zEigWz2NGTIVeuEdqgd7TIM0dGqSc+BZUkSXoldvRm6BtrbL9y3dpFtLeFtCNpEmlZAQ0htAOfBu4ENgEfCCFsOumyXwWeizHeAHwI+IPmvZuBnwO2ATcC7wwhbGze8zCwuXnPXuBXWvUZNL3Nnd3O77xvC+V65IXeLNQjlXrCc8dG6R4ppB1PkiRpykmSyM6ecQZyjRFl9964JuVEmmxa+QR0G7A/xngwxlgBvgjcfdI1m4DvAMQYdwMbQggrgWuBJ2KMhRhjDfgecG/zum81jwE8Aaxt4WfQNPfGq1fwlms66R4t0pspEYChXJkfHB6hVK2nHU+SJGlKGc6VOTSUI1dq/HP9vpvdfkUv1coCugY4dsLr7uaxE20H7gMIIWwDLqVRKHcArw8hLAshzAfeAZzqv96fBv7nqX54COEjIYSnQwhPDw4OntcH0fT22++5kbkdbezoyVCPkXoCO3rGOTiQSzuaJEnSlNLVl6F3rEipFrl82XwWzJuVdiRNMq0soKca7H3yxLrfApaEEJ4DPg48C9RijF3Ab9MYbvsNGkW1duKNIYRfax77q1P98BjjZ2OMW2OMWzs7O8/rg2h6W75wDp94y0ay5Tp7+jLEGBkvVnjy0DDZUjXteJIkSVNCjJHtx8boz5SJwNs3r0o7kiahVhbQbl761HIt0HviBTHGTIzxwzHGLTTmgHYCh5rnPhdjvDnG+HpgBNj34n0hhAeAdwI/Fd0zQxfAz7/hcjYsm8++gTz5co0kiewZyLGrN5N2NEmSpClhvFDhwGCeTHP47U+92u1X9ONaWUCfAjaGEC4LIcwG7ge+duIFIYTFzXMAPws8EmPMNM+taH5dT2OY7hear+8Afhl4d4zRlWJ0QbS1tfEf338j9Qgv9GQIIVAo1Xjq0IjbskiSJJ2F3f1Zjo/lKVbrLL9oNmuXXJR2JE1CLSugzYWCPgZ8E+gCvhRj3BlCeDCE8GDzsmuBnSGE3TRWy/3ECW/x5RDCLuDrwEdjjKPN438ELAQeDiE8F0L4TKs+g2aWWy5dyjtvWEV/tkzPWJEkJhwZzfPU4WEXJJIkSTqD7cfGGMiUqSXwlmtXpB1Hk1RHK988xvgQ8NBJxz5zwvePAxtPvq957vbTHL/yQmaUTvTv7rme7+0eYGdvls4Fs6nWIs8eGWP1xfPYsm4xIbiPlSRJ0snGixX29ecYyjfWz/hnr92QbiBNWq0cgitNORfPn80v3XEtxWqd3X05IDJWqPDkwWF6x4ppx5MkSZqUuo5n6R0rUKzWWTyvg2svuTjtSJqkLKDSST70mku5ZtVCDg8XyJaqtLcFDg3lefLQCIVK7cxvIEmSNIMkSeSZI6P0j5eoJfCGjcvTjqRJzAIqnSSEwH96/40E4Llj4wQiSYxsPzbKc8fGSBIXXpYkSXrRcK7MgYEsI4XG8Nufft3lKSfSZGYBlU5h0yUX895b1jCcr7KrL0NbgHy5xpMHhzk0nE87niRJ0qSxvXuMvvEi+UqNhXPbuXH9krQjaRKzgEqn8Zvv3szqRXN4oTvLQLZMjNAzWuTxA0OMF6tpx5MkSUpduVbnh0dHGciWqNTh9isdfquXZwGVTmP+nA7+7IGtzGoPPHlolGotIYkJO3szPHVomFo9STuiJElSqg4O5jg0mGc411gn4+dff0XKiTTZWUCll3HdmsX8xrs3UajUeeLwCBAoVes8dWiEvf3ZtONJkiSlpp5EnjgwzFCuRLZcpXPhbIff6owsoNIZfPDWDbz7xtUMZCvs6s1AhP5MiUf3DzGYKaUdT5IkKRVHh/N0Hc8wmK1SS+Dem9akHUlTgAVUOgu/+74tXLF8Prv7c/SNF0liZG9flscPDlOu1dOOJ0mSNKFK1TqP7BtkKF9htFABHH6rs2MBlc7C7I42/vynt7FgdhvPHB2jXK1TrtV55sgozx5xaxZJkjRz1OoJTx8eZmdvhnypynipxvWXLGTZgjlpR9MUYAGVztL6pRfxu+/fQqUeefLwKElMGCtUePzAELv7MmnHkyRJarkYIy/0jPEPe4aoJQnHRooAfOInrk45maYKC6h0Du7YvJoPv/ZSxoo1dvZkSSJ0jxV4dP8QPaOFtONJkiS11J6+LN/uGiBXrlJPEgZyZVYtnM1PbFqZdjRNERZQ6Rz933dt4sY1izg8UqRntECSNOaDfm/PIOPNORCSJEnTzbGRAg939TGYKbNo7iyeO5qhlsDH37Ix7WiaQiyg0jlqb2/jzz/8KpbM72B7zziZUo1aPWF79xj/sGeAUtVFiSRJ0vQylC3xzR19HB4qsHBOO3v7MhwbLbBh2Tw+eOulacfTFGIBlV6BpQvm8l9+6maI8Nj+YQqVGqVajWeOjvLkwWHqLkokSZKmiWypyjd39dHVN86C2W0cHM7z2IFh2gP8/k9uIYSQdkRNIRZQ6RV6zRWd/Lt7N1OqJXz/4AjlasJ4ocpjB4Z4vnss7XiSJEnnrVSt852ufp49MsbcjjaOjBZ4dN8QMcJ//eltbFm/NO2ImmIsoNJ5uH/bpfyLt24kV67zxIERytU6x8eKPLJ3gENDubTjSZIkvWK1esKjewf5fvNp5/GxEo/tGyaJ8OcffhWvvbIz7Yiagiyg0nn6+Fuu4p/eup6xUmMIbrWesL8/x7d39TOULaUdT5Ik6ZwlSeTpwyP8/Z4BYpIwkCnxv/YPkUT40we2cvtVK9KOqCnKAipdAP/PPZu547qVDOaqPN89Ti2J7OzN8PCufnKlatrxJEmSzlqtnvDU4RG+sbOPQqXOYK7MYwcaTz7/+J/cwpuvccsVvXIWUOkCCCHwRx+8mVdduoTe8TK7jo9TrtZ49ugo39rVT6lSSzuiJEnSGVVqCY/uH+KhF46TLdYYyZd5bP8wtQQ+/VM38bbrVqUdUVOcBVS6QDra2/h/f/ZWNq64iCMjJfYN5MmXazx1eIRvdw1QrSdpR5QkSTqtUrXOt3f18c0dfRSrNUaLFZ44OEw9gT+4fwt3bL4k7YiaBiyg0gU0d1Y7f/1zr2b1ojnsH8xzZLRAtlDh+weH+O7uAbdnkSRJk1KmWOXr23v5zu4BarU6A5kyj+4bpFaH//j+G3nXljVpR9Q0YQGVLrClC+fypQdfw5L5HezoyXJ4uMBovsKj+wZ5dP8giSVUkiRNIkPZEl/9YTeP7R8iJgl7BnI8sneIWh3+/Xs3c+/Na9OOqGmkI+0A0nS0bulF/PVHXsP7/+RxdhzPUo8JMYl8d/cAczva2XbZUjdtliRJqesZK/C153rp6s2QJHWe7R6nZ6zMxXM7+LMHtvKqy5alHVHTjAVUapGrVi3ivz/4Wt7/J9+nqy9PkkQikYd39TGno40t65ekHVGSJM1QMTZW7P/Wjj72D+bIFSv88FiGXKXG5ksW8pc/vY2lC+amHVPTkAVUaqGNKxfy5Qdv471/8n32DBSoJ5BE+J87jjN3djvXrFqUdkRJkjTDlCo1/tf+IR7bP0T/eJHj40Ve6M1ST+CDt67j3969mbY2Z+qpNSygUotdvmIBX/mF1/K+z3yf/UMF6jESga8/18usmwNXrFiYdkRJkjRD9GWKfPOFPp7vGSNfqtHVl+HISIk57YHfu/8G3r3F+Z5qLQuoNAE2LF/AV37hdbznM49xaLhIPUIg8JVne7h7yxquWmkJlSRJrRNj5PnuMb65o4/Dw3mGMmV29mXIVxJWXzyHv/rZW7m803+PqPUsoNIEWbdsPl/96G2859Pf5+hIkRgjAfjKM8e456Y1XLP64rQjSpKkaahYqfGdrgEe2z9E33iBPf05esfLBODuLZfwO++5gTmz2tOOqRnCAipNoDWL5/O3H7+Ne/7o+xwbLVFPIvUk4cvP9HD3TbB5jSVUkiRdOL1jBb62vZftR8foHimwbzBHqRZZtWg2//n+m3jV5f+7vTsPsyyry3z/XXs6U8w5RM5jZWVVUSMkVRQlBVKCigKCoDiBiFRji62iV237ud3t020LttcLt1W4NCIq+qg4XJEZikEZiqKKmsfMqpynmIcz7mndP9Y+EZFZWWNmRGRGvp/HIOLss88564S7MuKN31q/tXq5hygXGQVQkSU23Ffhn3/ppfzYh77FgfEWrTglznKy71qSLOc6dccVERGRs9RJMr7x+Bhfe3SUvSenefREnfFmiu/Bz9ywhf/6uhfg+2o0JEtPAVRkGaztrfCpd7+Mn/3Yd7jz4CR3HpgkyyxZZknznBdv055bIiIi8vw8PjLLZx84zgNHpnn05AyHJtpkFravrvLHP3ktV2zQH7tl+SiAiiyTnkrI3976Et7zd/fwz/ce59sHJmnECZnNyXLLSzQlRkRERJ6DeifhtodO8vV9Yzx4dJonRhu0M0vJN7zr5h38+qt3Y4xZ7mHKRU4BVGQZ+b7H+99yHVuHavyvr+zjvqN1OklOlkOc5tx86drlHqKIiIic56y13H9kms/cf4zvPDHGIyMNGnGOAV6xazXve/NVDPdVl3uYIoACqMiyM8bwnu/fzZZVFX77n+7nkZEmjSQjyy2dJOeWy9dqM2gRERE5o+NTLT7/4An+7bER7jw4xXQ7BeCy4R7+4M1Xc+UmTbeV84sCqMh54k17trB5qMI7/vwuDk92aD96gk6aMt1KeN21G4kChVARERFxphodvvzoCLc9fIJv759krJ4AMNwX8p9/+AX80NUbl3mEImemACpyHrlhxxo+/e6X8eMf/hYnZjvc9vAI9XbKTDvhx/Zsoaes/2RFREQuZp0k4+t7x/j0vUf42r4xxhuu4tlT8rn15u28+3t3aeaUnNf026zIeWbrmhqf/5WbeetHb+feo7N87bExJhoxM+2Un7p+C2v6yss9RBEREVliWZbz3UNTfOLOg3zxoREmW/PB86eu38J7Xr2bUugv8yhFnpkCqMh5qL8W8Q+/8D38l08+wF/fcZi7j8ww0YyptxN+6oatbF/Ts9xDFBERkSVgreW+o1P87R0H+cz9J5hqZQD0lX1+5sat/Motuwm1TEcuIIsaQI0xPwB8APCBj1hr33va/YPAR4GdQBv4OWvtA8V9vwy8EzDA/7bWvr84PgT8LbANOAD8mLV2cjHfh8hyCAKP333j1dywYxW/+ff3cXCizSfuPMxMM+anbtzGNZvVVEBERGSlynPLvYcn+eBX9/H1x8doxhaA/rLPz920nXe/che+r+ApF55FC6DGGB/4Y+BVwBHgO8aYT1prH1pw2m8D91hr32CMuaw4/xZjzJW48Hk9EAOfM8Z82lq7F/gt4DZr7XuNMb9V3P7NxXofIsvtdddu5KqNfbz1o3dweLLNP95zjLFmzM+8ZBuv2L1W+3mJiIisIHluuf2JMX7/c4/wwNEZUpc7WdMT8o6XbuedL9+p4CkXtMWsgF4P7LPWPgFgjPkb4PXAwgB6BfB7ANbaR4wx24wxw8DlwO3W2mbx2K8BbwB+v3iOVxSP/3PgqyiAygq3fU0vX3zPK/h3f/kdvvbYOF9+ZIyR6Q5jszE/ct0GwkBrPkRERC5knTjlU/cd5QO37ePQZBsAD9i9tsZv/MBuXnn5Ov3RWVaExQygG4HDC24fAW447Zx7gTcCXzfGXA9sBTYBDwC/a4xZBbSA1wB3Fo8ZttYeB7DWHjfGrD3TixtjbgVuBdiyZcs5eUMiy6kc+nzs7TfwR7c9xgdu28cDx2f5/c89zNHpJm+9cRtDtdJyD1FERESeo/HZNv/jMw/xuQdO0khyAMIAbtqxiv/y2ivYvqZvmUcocm4tZgA9059o7Gm33wt8wBhzD3A/cDeQWmsfNsa8D/giUMcF1fS5vLi19sPAhwH27Nlz+uuKXJCMMfzS9+3mRdsG+cW/upvRRsIffXkf+07O8B9u2c2l6/RDSkRE5EJw20PHee9nHmbvWGvu2GA14Eeu3cCv3rKLvpq63svKtJgB9AiwecHtTcCxhSdYa2eAtwMYN6dgf/GBtfZPgT8t7vsfxfMBnDTGrC+qn+uBkUV8DyLnpZdespYv/9rNvPMv7+LOg9N86v4RHjvZ4D++5nKtCxURETlPTTbavPezD/Mv956g2a12erB7uJdfffWlvOLStVrfKSuesXZxioPGmAB4DLgFOAp8B/hJa+2DC84ZAJrW2tgY807gZdbatxb3rbXWjhhjtgBfAG601k4aY/4nML6gCdGQtfY3nm4se/bssXfeeefTnSJyQbLW8kdf3ssffXkvnQxqkeHd37uLt9+0nXKkXZZERESWWzvJ+LNvPM6ff+MgJ2bjueMDZZ/XXrWOX7hlFxsGass4QpHFYYy5y1q750nHFyuAFi/6GuD9uG1YPmqt/V1jzLsArLUfMsbcCPwFkOGaE72ju6WKMebfgFVAArzHWntbcXwV8HfAFuAQ8GZr7cTTjUMBVFa6+w5P8osfv4vD0x0AXrl7Fe/70WtZ06fpOyIiIkstzXI+c+9RPvjVfTw80pw7XgoMl6yu8a6X7+D7XrCeiv5YLCvYsgTQ84UCqFwM6u2E3/jEvXz+oZNkFoZ7I/7gzVfzskuHl3toIiIiK14nzfj0vUf4s28d4qEjM2TF8cCD4b4yr796PT/xki1sGKjhe1oqIyufAqgCqFwErLX8zR0H+L3PPsZMO8UD3vDC9fz3cTfNfwAAIABJREFUH7laf2UVERE5x5pxyj/ffZiPf/swDx+bJS+O+x70l0Nu2DHEW1+yhWu2DFHVz2G5yCiAKoDKReTxkVl+8+/v4a5DM1hgbU/Ih37mRbxw66rlHpqIiMgF7eB4nU/ceZjP3n+cxxd0sA086CuFXLmhh9deu4nv2bWG4b4ynqqdcpFSAFUAlYtMJ8n42Ncf50++9gTTbTcR6E3Xred9b7oG3/eXeXQiIiIXhiTL+bfHRvmHuw7z9cfHmW7N7wwY+tAbhexeX+OWy4a5+dK1bF1doxTo56yIAqgCqFyErLU8dGya3/mXB7nz4BS5haGKz0fetocXblu93MMTERE5L52YbvLPdx/l0/cd58Fjs/PrOQ1EgaGvFLFzbZU92wZ52aVr2T3cR28lXNYxi5xvFEAVQOUiNttO+OtvHeCDX3uCqbb7y+3rr17H//Xj1xFovzEREbnItZOMf31slM8/eIKv7x3lZLFdigEi31AJPYZqETvW9HD1pn5evH0Vu4Z7WVWLtPe2yFNQAFUAlYtcllvuPzLF7376Ib57aIrMQjU0vO9N1/DaazYu9/BERESWTCfN+NbjY3z+wRN8c984BydOXctZCjyqUcDa3pAtgzUu39DPdVsHuXS4l7W9JYVOkWdBAVQBVASA0dk2f3v7Af7s9kOMNxIAdq+t8b/ftoctq3qWeXQiIiLnXivO+NYTY3z1kZN8c984T4w15zrWdgNnJfDprYQMVAK2rqpxyXAPV24cYPe6Xtb1lRU6RZ4jBVAFUJE5nTTjviOTfPArj/PNx8dpp+7fgZ988Ub+8+uuohyqeYKIiFy4xmbbfOXRUf71sRG+e2iKo1PtuftCD8qhTzn06S0F9JY9Bmsl1vdXuWRtD1es7+OS4R7W9JaX8R2IXPgUQBVARZ7k5EyLLzx4nD/9xgEOT7TIcqgEht9749W8/rqN+muviIic96y17Bup8+VHTvLNvWPcd3SayQWdakuBoRL6VCOfniigFhl6yhFDPWWG+0pcsraXS9f1sG1VjYFqtIzvRGRlUQBVABU5o06a8fDxGf7qW/v54kOjTLdTLHDpmhrve/M1XLdlcLmHKCIiMufkTItv7hvn9ifGufvwFPtH6yTFfFoDVEKvCJs+/dWQcuhRDnyq5Yihasj6gQqXrO3h0uFeNg5UKGnWj8iieKoAGizHYETk/FEKfK7dPMjGgSo37DzJX91+kEdO1HlstMEb/uSb3LRzkN9749VaHyoiIktuphVzx4EJvrlvnLsPTbBvpMFsJ5u7PwoM1SigVvLoL4f0lwJKUUAQGEphQF8pYKinxNZVNXasqbJ1VY3VtRKepxk+IstFFVARmROnOY+emOH/u+coX3rwJCdm23SK9aFvum49/+drr6K/qn3ORETk3JtsdLjzwCR3HZzgvsNTPHqyzngzmbs/9KASBdQin4FqQH8lohZ6BIGPMYbQ9yiHHoPViA0DFXasqbF9dQ/DfWWiQFuOiSw1TcFVABV51iYaMXfsH+df7j7CnYcmmWgkJLn74X/rzTv4pVsuVaMiERF5Xqy1HBhvcNeBSe48OMkDx6Y4NN5ipj2/bjPwcOs2Sz795YjVvRE9JZ9K4ON5kOXgGQ/f9+ivBKztK7F1qMbONTXWD1ToKQXqYyCyzBRAFUBFnpMstxwYr/O1R0f43P3HeXy0yUzbBdGeyOPfv+IS3nHzDkqBgqiIiJzZWL3DfYenuPfIFA8dn2HfSJ2jE03ifP6cKDDUQp9K5DNQCVnVG9EXhZRLPiXPI84htzmBZ7DG0FsKGKpFrB+osHNNjc2DVforIYGvKqfI+UQBVAFU5Hlpxin3HZ7itodP8I19Yxyf7lDvpCQ59JY8fuHlCqIiIhe7yUbMfUenuOfwFA8fm2HvyTrHplu0kvmk6QHlyKMaBVRDn6FayOqeMgPVgFLoE/keubXEWU6WWcDgGeivuvPWD5TZOlRlXX+Zvkqonzsi5zkFUAVQkbNyYqbFtx+f4CuPnOCh4zOMNRLq7YQ4g76yz797+SX8/Mu26xcCEZEVKs8thyYaPHhshoePTfPYyToHxpscmWzSPEPQrIQ+tSigt+KzqlpisBbSUwqpRj7GWIzxaCcZSZZjLVgg9D2GaiGDtRKbBqpsGaqwurdMbzkgVIVT5IKiAKoAKnLWkixn78lZ7tg/wbceH2X/WIOJRspsO6GTWQVREZELnLWW8XrM3pFZHj4+wyPHZ3h8tM7hyRZj9Zh8wa+NHrgtTqIiaJYDVtciVveGVKKQgXJEKfRIc4u1liTLaKeWOHVh1RjoL0cM1EJW1SI2D1VZ31+mtxzSUwo0pVbkAqcAqgAqcs4045T7j0xzx/5xvntwkmNTTSZbKTOthHZq6S/73HrzTt7xsh1qViQicp7pJBkHxhvsPVln38gs+8caHBxvcny6xUQjZkExE4BSYKiEvmsKFAX0VALW1EoMVn3KYUBfJaKvFIAHaZaT5ZBkGc04p1OETWstlVLAYDVkoBKxrr/MlsEK/dWInnJAJfTVNEhkhVEAVQAVOecmGzHfOTDOXQeneOjoFGONNlONlKlWSivNqUYeP33DNn7xlZfQX9H2LSIiS8Fay1i9w96Ts+wbqbNvtM6hsQZHptqMznaYXtBttqscepQDj3LgmgG5hkABa3vL9JZCosijJ/LpKYdEvkea5+QW0szSiFOasZtK21UNA/qqAX3liFU9Iev7y6zucVNpVd0UuTgogCqAiiwKay1HJ1t8e/8Y9x6e4fHRGaYaMdPtjOlmwmycEXrwhhdu4tdefSnDfZXlHrKIyAVttp1wcLzJoXFXuTw40eDYVIuTM23G6zFTrYT0tCpm4EG5W8UsPtfKPv3ViFXVkFoppBS6dZvdiqRnDGmWk+SWLLc044xmnBEnGd2n9z2ohAF9lZC+SsCqWon1/WUGqhHVyKdWCjQTRuQipQCqACqyqLLc8uiJGe7YP8HDx6c5Otmi3k6YbqdMNRMmWyke8Oor1/J/vPoydq7tXe4hi4icd2ZbCcemWxyaaHJwvMGh8SZHp5qcnIkZnW0z1XRr7k8XeFAqKpjl0KcUevSUAvrKAQPViP5KSK3kE/k+5dCjWg6oRgGhMVgoKpqWTmppdFJacUontWS5axDkeVAqmgr1lAJ6yiGraiHDfW7NZq0UUCv5Wv8vInMUQBVARZZEO0m55/A0dx+a5PGROqOzLeqdjHo7ZbyRMNaIAbhh2xC/9H2XcNPO1Vr3IyIrmrWWmXbC8ak2RyebHJ1qcXSqzbGZJqMzHcZmY6aaMdOt5EnrLwF846bIlgJvvopZTJPtLQf0lSKqJZ9q6FEK/eKcgErkUwp9Ag9yCxZXyYyLkNlMMpI0cxXOItR6HgSeTyk09JZCesoB/ZWQ1b0lBioRtZJPNQyolnx1pRWRp6UAqgAqsqRmWjHfPTjF/UenOThWZ6LRoRlnzMYZE/WYkdkOmYWtqyrcevMOfmzPFv0yIyIXjFacMjrb4eRMh5OzLU5Mdzg502ZktpgG20yYbMbMthManYwzFC0BKAeGKDBEfkAp8KiEHqXIoxL49JQC+qsu9JUDFzyjorpZLsJo6HtEgY9vwPMMeW5J85xWktHouI8kz8lyS5q5iqZvDIFnCAKPyPeplfyiqhnQX4kYqAbUotCF2shXgyAReV4UQBVARZbFeL3jpuUem+bIVJvZdkw7yWjFGRONmGMzHdpJTm/J5y3Xb+GdL9vB2r7ycg9bRC4C1lrSLGeimTDZiBlvdBivx+6j0Wa8njDVcmFytpVSjxNm2yn1jpue+lRC3xD5hijwKfkeUeA+ypFHNfColiJ6yy70hb6rbJYCnzBwazDLoUclCoqA6cKiZ4wbb25Jspwks3Mhs51mJGk+d3+aWwyW0PfxPTeOKDDF9Fm3LrO/GtJfiYrXc0GzHLrzRUTOBQVQBVCRZXV8qsW390+w9+QsI7NtGp2UTpLTTjKm2gnHplpMNlMMcOPOIX7+ZTv43t1r9Vd3EXlKee72lJztuGA400qY6bjP082E2barQrqPhJlmykw7oR4nNNquoU779G49p/GAKDAEviH0PUKvCJOhW+9YCT3XQTZ0wa4a+QS+R+gZSoFHEHiUfJ8odCGzJ/KIwu45HqHvwqVnILeWNIdOmtFO3PhaiftIixKq20/T7avpeS7oGmPwjSH0DZVS6Jr/RG7LlO7az+7U3UrkqqmegqaILDIFUAVQkfPC/rEG335inMMTTcbqHZpxQpxZOkVV9PhMh2NTLdIc1vZG/PiLN/NzN21nsFZa7qGLyPOU57YIVzntJKfZyWgm6fznOKfZSZlqJcy0YmY7KdPNhHonpdFJqbczmsVWH+3UBbJOktFJcp6mEHkKz8NVJT0XCiPPFJVJf64CGfkuTJaKQFmOXAXTeBB4HoFnCD1vLoxGgUfoe1SKsFcKinO6n32Dbzwsdv7fuSQjTt3+mO0ko53kdNKMLLeAIc9z0txNlQ18t37TL0KqC6welcijGgZFqPXc2sxy4BoQRfPrRLWsQUSWkwKoAqjIeSPPc54Ya3LngXEOjbeYasV0UvcLZZpDnOaMzHY4PNlkuuiee/Olq7n15h3cqKZFchHJiu0vwDWQAbAWbG7JrCW37pzcFrdzi11wLLXudp6fdr+1xFlOJ8lJspzOglDUOeV2RiezJIm7L07dXo9xZt3tLCfNiimh3fOL43HinivJMpLMVe3y5/ErR9ANf373czcIdoOeC2Vh0K1QGiLfw/c9yoGbvlot9p10j8WtgSyexzPucb7vHhf4bn1lt5lP6Ju5abC+7xEYU/wbVEyHTXPaRZiMUxceO4n7PrTTjDjJOf1tp1mOMRB6Hp5v8A1gDH4RcivhfMCMAo9aKaAnCiiFRVgOuxVYj8j39G+iiJyXniqABssxGBG5uHmexyVre9i5psa+kTp3Hpjk0ESDRjslyd0vrtXIZ8tghdlOyqHxJv+2d4yvPjbGxv4SP3H9Ft760m30VaLlfisii+bQeIMf/MC/0ozdFNHz6c/FRV7CM+Bh8Dw3hdT3zCkfpcCjFrl1hYHnuc9FmHOBrjhuIAw8Is8jDF0jntDzCYNieqoHofHwPDCe547NjcGdExTPFfhmfgps5KqaJd8jLIJb4Jm5MQceGGMwGHJrSTL3R7AktbRSt0ygFedzfwjohvBOmhf7bNri+2GwuNBvczCeC5W+ZyhHPtbi3nMResuhX1RP3bFq5NNb7JdZCvyi461rNhT5nqbLisiKogAqIsvGGMOu4V52rulh78gsdx+e4vBEk0Ynpcda4iQjCjwGqhFXbern8EST/eNN/uCLe3n/bXt55WXD/OxLt6oqKivSvYenaMQ5w70l+ish4AJXN/y54ERx3AXA7jnzX7tw1X0c3cdTTOvE4BXTO7udUY2Zb3rjpn+6Dqt+UanzPA+veCLP2Lkn9Siex1v4unbuv02DC4tufC5M+p57XW8uoLoqp1808Ql8F8CCYrps6Bfn+O75u4/vvleMqw4nWbHdSJYTFxXZdpIx005IMxces9xVbztZRpLOV5q7Y+3estbOfe1Ctgu5UeDGQPG99IuqbHecbgquR7lYj1kKXRjuNiQq+f7c12r8IyIXEwVQEVl2nmfYva6PXWt7OTjR5N7Dkzwx1mC2nVIru7VTrTjlkjU97Fzbw1QzZt9Igy89dJIvPHSSNb0Rr716A2+7aRtbh2rL/XZEzomxegeAXWtr7Nk2NBd2uiFyodP//mKM1z31lPDpPrkvPO/U6uXCYOjhgmY3NGKYq3DOhcwiAHZvd6ePBnMhsbhdPD4sXsQtS3Sj6O5NWfyfm0pcVBuTzAXIpAiQjU5GnLkAmeWWzLoQ2T0/7QbPLC++R2a+amy7n6wLrMYW35tupdZgPM89JrdQBE0XzMEvKqpBETS7FdZK5Pbb7FYqu4Gy29lW1UsRkSdTABWR84bnGbavrrFtVZVjUy3uPTLF3pN1ZloJlWqJLM9pJTlhb4lVtYjrdwxxYKzB/rEGH/3GAT76jQNcsb6XN+/ZzI++cKOm6MoFrdVxU2+3rq7yums3zlc/i0BjbBE8DWBNcZ87jjs094XFUixbxNUsuzU+F8asPS3UGk45Zq0lt4bcuqphWmxfkuWQFIEwzd3WIO00cefk7v40y+fCpbXuc3eqas6pITKz1r2X+ayNtcWbtca9J1tM/S2Cs6U45hkqgcFEwfz7s+41TFHN9YxXVI/dOtCFU4C7IbNcbIMSBcUemwvXl3YrskUTIhERee4UQEXkvGOMYeNglQ0DFV68rcN9R6Z55MQMk82E3rIrzySJpdnJ2DXcwwvW9zHVStg3Umf/WIPf+ZeH+O+feogXbxvizXs285qr1lGJ9M+dXFiaSQrA6EzMA0dn6LYhyhdMFbULPtsnl/sWTCPtHlrQyAgXCLshzT2HO9ZtUDgfGhecz/wU1fmguCDQzn1mLux5Bqwtmu3gwrHppskiHJZ9H+PPr6d0Q3V7X3arvp5n5iqXpniB7lpL33PV3O6UXt+4abLlwCcK/PltVBasEQ26YbK7dtQzaDq/iMji0m9kInLeMsawprfMLZeXeeGWQR48Ps1Dx2YYm+0Q+h7D/QFJltOIU/rKIS/aNshLdgxxZLLN46N1vntoktv3T/Db/3QfL7tkDW/es4lbLh8mUOVCLgAz7RiAVpwwOtsuQpedy3gLg1I38i2sHMJ8t1aLm1bbLYN2w1u3ahosCI7uf7251Np9Dop1ptalR+aH0o2MT66mep5xt4sqrS2mvnYrlt6CJkJnut1d+xn53fWfC5sYeUUFs5gWO/fZm7ut6a8iIucfBVARuSAM1iK+55I1XLNpgEeOz/DgsRlOznTIrWWwGoE1dLKMZidj3UCZzYMVjLE8Md7i8ZFZvvrYCF96ZIRq6PE9u9bwxus28srL1xIF/nK/NZEzqncSwK3V7K+ECxZ6ntospxsaT+2TO3+uedKxbiQ9vWrppqd2/zzTbRjUbTDU7TjbnQrsFdVCz8yvIe3uldndssQvOtie2myou2504RrSBR/FY7pBUkREVhYFUBG5oPSWQ168fRVXbOjnibE6Dx2d4chkkzjLqAQ+fQMhaWaZbackWc6O1TV2D/eQpJa9o7PsH2nw5Udc86JK6HHTJat5w3UbeNUV6xRG5bzS3X6lrxyxa7gX6DYUmp/e2rVwKuzc12c6Nnd8vtrYDZYwX5U0xf0sqEiaBfd3g2N3LWW3YikiIvJMFEBF5IJUKwVctXGAS4d7OTbZ4pETs+wfqzPVTAk8WFVz21Z00pzZdoIxHles7+PqjQMkac6jI7M8PtLgq4+O8KWHRyiHHjftXMUbrtvIq65YRylUGJXllWYugK7tj7hm88Ayj0ZEROTcUAAVkQtaKfDZvqaHLatqjNU77B2ps390lhPTHSxQDTzWD1SwFurtlEac4hnDlRv7uHZTP0lueeT4LHtH6nztsVFue2SUcmC4cedqfuDKdbz6BevcFF+RJZakLoCWPP2oFhGRlUM/1URkRfA9w3BfmeG+Mldt6OfoVJO9I3UOjTcZnY0JfUOt5DNUq5LlOZOtlHaS43lw1aZ+rtsyQJzmPHx8hsdO1vn63lG+8ugov/UP97N7XQ/fu3str7lqPVdu7FeXTFkSWdHtNgrVNEtERFaORQ2gxpgfAD4A+MBHrLXvPe3+QeCjwE6gDfyctfaB4r5fBX4e1xnhfuDt1tq2MeZa4ENAGUiBf2+tvWMx34eIXFj6qyH91X52ru1hZKbD/rE6B8aajNbbzLQyosAwVI0oBYZWkjPVSmjGOb6BqzcN8MItAyS55bGTdR4frXN0osUHv/YEH/zaEwxVQ16yYxXf/4Jhbrl8HT1l/R1PFkdaBNCKAqiIiKwgi/abkzHGB/4YeBVwBPiOMeaT1tqHFpz228A91to3GGMuK86/xRizEfgPwBXW2pYx5u+AtwAfA34f+B1r7WeNMa8pbr9isd6HiFy4SoHP5qEqGwcqXLEh5uR0m/1jDY5ONTk52yb0PMqhx3BviSgwzLRSZjspSQKeZ7lyYz/XbBrAM4bDkw0ePVHnwHiTzz14gs88cALf3MsLNvbxyt3D/PDV67mkaBQjci7ERQAtqzmWiIisIIv5p/vrgX3W2icAjDF/A7weWBhArwB+D8Ba+4gxZpsxZnjB2CrGmASoAseK4xboK77uX3BcROSMPM+wuqfE6p7SXFX02HSLQ+MNTsx0ODbdpuQbaqWQ9f0VAs8w1YxpxBm5tXgG1vdX2DxUI/AM062YR47N8vhYgydG69x3ZIb337aXNb0R128b4pWXreWVl61lsFZa7rcuF7A0KyqgkQKoiIisHIsZQDcChxfcPgLccNo59wJvBL5ujLke2ApsstbeZYz5A+AQ0AK+YK39QvGYXwE+X9zvAS8904sbY24FbgXYsmXLuXlHInLBK4c+W1ZV2TxU4dLhXk7OtDky2eToZIupVsxMKyYKPfpKEauGInxjmGgmNDopcZphjCX0fa7fOcSNl6wizS2PjzZ47MQsh8YbfPaBE3z6/hMYYOeaKjdesppXXb6Ol2wfIlJnXXkOFEBFRGQlWswAeqYuHfa02+8FPmCMuQe3zvNuIC3Whr4e2A5MAZ8wxvy0tfbjwC8Av2qt/QdjzI8Bfwp835NeyNoPAx8G2LNnz+mvKyIXOWMMQ7WIoVrE9tU1Rmc7jMy0OTrd4sR0m5F6h6mWIfA9essBGwYrlH3DRNN10m20M/c8nuGStTUuX9eL8WCyHvPoyTr7Rxscnmryl986xF9+6xClwHDVxn5u3rWGV1y2lis39OF5WtsnTy213TWgWmcsIiIrx2L+VDsCbF5wexOnTZe11s4Abwcwrq3k/uLj+4H91trR4r5/xFU6Pw68Dfjl4ik+AXxk8d6CiFwMyqFbK7ppsMKOZg8jsx1GZ9ucmOkw0YgZnW0z3fTwjKG3EjLcV6Ia+Ey2UxrtlHac0SYnx1Iphdy4cxXfs2s1nTTnxFSLR0/OcmCsycPHZrjz4BR/+KW9VCOPy9f3sWfrIDddspo9W4eolhQ0ZF5WVEBrkf5QISIiK8di/rbzHWCXMWY7cBTXROgnF55gjBkAmtbaGNfx9l+ttTPGmEPAS4wxVdwU3FuAO4uHHQNeDnwVeCWwdxHfg4hcRIwxDNYiBmsRO9bUGK/HjNU7TDViTsy2mWrGjNVjZpoxnu/RUwoYqoX0lCrMtFOaSUqznTGb5eR5ju95bByosm11DWOg3krZP9lk/0idI5MtHjsxw10Hp/h//3U/gQc719R40dYhbty5mht3DLG6t7zc3xJZRnnRhKha0hRcERFZORYtgFprU2PMu4HP47Zh+ai19kFjzLuK+z8EXA78hTEmwzUnekdx37eNMX8PfBe31crdFNNpgXfipu0GuK1bbl2s9yAiF6/Q91jXX2Zdf5l2kjHRcGF0pplwcrbDdCthvBEz20owxtBXCegph2zoq9BOM5pxRiPOaCUZ9U4OxlAKDFdu6OO6TQPk1lLvpJyYbvPEaJ1DEy1OTLf46zsO89d3uOXzGwfKXLupnxt2rOJ7dq1h++qa9iC9iGQ2B6CnFC3zSERERM4dY+3KXx65Z88ee+eddz7ziSIiz6AVZ4zVO4w3YuqthJHZNtPtlGYnJfQM1hh6SwE95YD+Sugek2Q0OhmNOCXNcizgGUM59KmGHoHvYXMXSCcaMXtH6hyabHJ8qsXJ2Q6ZyyEMVEIuW9fLNZv7eeHWIV60dZDVPeq0u1Ld9N4vcXSqw9d/4+VsGupZ7uGIiIg8J8aYu6y1e04/rgVHIiLPQSVy60U3D1Wpd1LG6x3G6jGtTsJoI6bezpjpJK66OdOmHPr0lQL6KgHDfSU6aU6S5jTjlFaSM9FIsFg8z1AJfDatqrFtTQ1wYXeqmfDEWJ2DY02OTDW5/+g0t++fwC2Xh3V9JS5b18t1mwd44dYhrt0yQG85XMbvkJwrefGHh5KvNaAiIrJyKICKiDxPPaWAnlLAliKMTjRixhsx7Thjshkz005odDJOxm38ukcYGHrLAdUoYE2xvjO3lk43kMY54/UOWIvveZRDj+H+MpuGKthd0Ely6p2E0dmYg2MNjky1ODHT4vYnJvjqY2OA25tq01CFK9b3cfWmfq7eNMiVG/vor4SavnuByYsZSlGgH9UiIrJy6KeaiMhZMsbQWw7pLYdsXVWjsSCMtuKUmVbCTCel2U4Zr8dMELuAGfnUIp9K6LOmt+wCh4U4y2nEKa1OxuisC6RR4BP6HgPViOH+Mlds6AULs3FKK844Od3mwHiTI5NNTky3ue3hET734Mm5Ma7rLbFzbQ+Xb+jjyo39XLNpgC1DVXxPofR8lRVNiEJf/z8SEZGVQwFUROQcq5UCaqWAzUNVWnHGeMNt59LopDTjlEYno5lktDsZzU4CxhB4hkroUysFlEKf1bUStgYGS5xa6p2EVpIzXo9J85wo8Ah9j3Lgs7Y3ZH1/heu2DJJkObPtjFaScnSqxZEJVyUdmWlzz6FJvvH4+Nw4eyKfraur7B7u44oNvVy2ro/L1vUxVIvwFEyXXbcCqj8SiIjISqIAKiKyiCqRz6aoyqbB6lw33clmzGzbNSRqdlLaSU69k9JKMmY7Kb7x8D1LTymkHPmUAp9VPSUw4BtDmlmacUa9k9CIM8YbMWHgEfqGyPfoKfus6onYNFjl+m2WLLfMdhLiNGeqmXBgrMGxadfg6MR0m4ePzfKPd8+PebAasHmwyo41PexaW2N3EUyH+8uEWo+4ZLprQAN9z0VEZAVRABURWSLl0GfDQIUNAxWy3DLdSphqxky1EjpJTpxmtNOMVpzT6KTMdlImmwmh7yqktVJAqah81ko+veUA3zNYLJ0kZ7adUO9kjM66Kmkp8Ag9QxB49JRCoqphXX+F3ev6sNbSijPSYO34AAASyklEQVSaaUYzTjk+2ebYdIuR6TajjZgjUy0eODZDvqBRek/ks2GgzLY1NXat6WHn2l52rqmxfXWNWilUpe4cs7hvvqrRIiKykiiAiogsA98zDNUihmpuj8dWnDHViplqJsy0ErLc0k4ykiynEWe0k5zpVkKaW0Lfw/cM1cin5LttXCqRz1CtxLp+D2MMeW5pximz7ZRGnDEy2yHN3NTdyHchthr59FQCrC2xbagHjCVJLa0ko5NmdJKckzNtjk01OTnbZmymw0Qz5SsPj/KFB0fm3osBVtVC1vVX2DxUZeuqKttX19ixpocda2r0lUNVTp+HfLkHICIisggUQEVEzgOVyKcSVVjf76qjs+2EqWbCVCuhFWckWUaSWbI8pxnnJGlGvZ0ymedEvodnoBK6CqnxDOXQI/J9hvtKBL6HMZDl0IpTppoJzSKUJllOKXSh1PcMpcCntxTQXwkAw+ahKh6DxJkLxHGW00kzxhsJJ6ZajNZjxuttJhoJJ6bbPHpiliQ/dX/p/nLAqp6IdX1l1hcV4C1Dblry1lVVhmqRew+q9J0izy0G1L1YRERWFAVQEZHzjO8ZBqoRA1VXHW0nGdOthOmWq472li1JlpHnYK1bDxpnOc0kZaqVuwqpMfi+a2zkewbfGMqhTxh4bBiYX8uZW2h2UqZaLpTOtFNGZzsYA6XAJwqKx0Y+PeXAhSEL6/vhmo39pNYSpzlxmpOkljjLmG6lHJtqMlbvMFrvMNVImGmn3H90im89MYE97f2WAsOqWom1vSXW9ZfZMFBh42CV4b4Sw31l1veVGaiGRSdgc9EEMnv6N0pERGQFUAAVETnPlUOfcugz3FeeC5zdQDrbTilHgQukFsDSTnLaSU5mLVPNBICgCG6hbygFHuAaFpVDjyhwazsD38NgAEuaWWba7jXaRffdOMvxjKEcuD1NfeMeXwrcdF6LC03rB+DKjX3EmQulaZaTZDlp7sLqRCNmrN5mrB4z2XSvMd1MODTZ5MFjM0+qoAIEBvqrIf2VkFW1iNW9Zdb0lljbV2K4t8xwX5l1fWXW9pUoBS50X+hh1SqBiojICqQAKiJyATHGzG3zsmGgQp5b6sVeo7PtlHonpRRYeko5WQ6mx5Ln0Eoy0jynk+RMxgkeEBaBMwo9As+QA5Fn3JTcwKcSBQxV57dkMRg6Wea6+LZSOmnORCMhznIslsgv1pcG7vlKoUc58jD4xegt1rppvcbgQmlmSbKc3EKSu6A620oZr3eYbCXMtGJmWikz7ZR6J6Heztg7Uueew9NnDKoAldDQUwrpLQf0lkMGqiED1YjBashQrUR/JaCvEtJX3NdfCRmsRlQiH8+YuYrxck8JtoBmJYuIyEqjACoicgHzPENf2YUpcFWzVpIx206LD1fBLIU+SZZjqu5xaWZpJylJ7rrh5lg8oO15mDb4viHw3XRbz7iqaSn0KQUeQ9US6/oq+AaKgilJljPTdkGxFbtOvtMtt9UMuLAb+oYo8PCMIQoMke9TDj0qYQDGYm1RQe0DU3TqTXMXUnMgtzlZasmsJc2hnaSMN2ImGx2mW6mrCLcSmrGbjtyMM6Ynm+w7mdNKM54ir87xTbEWN3Thuxb51Mo+tSigpxRSDT3KJZ9K4FMtBVTCYvyRTzUMKIc+1eL+UlEZDk4L5cYYjHHfU0Px+bSQaYwL+880XhERkQuRAqiIyApijKEaBVSjgOE+dyxO3T6js+35KmngQRR4pJmlv2LxjCFJc5I8J82gk2bEaT4XktppBq2U0DdzoTPwXdOiUuARBR7l0GewGhL4Hh4GC+TW0kpSGh3XNKmVuGZKzTgjyWJya3HTgY0LqYGHbyDyfcLAEHgeYbeAig+Rqwy66mCZravcPVnupvimmcVai8Xtf5rl7utud9+ZVkIzyWh03PYzrSSbC8ytJKWdZLRTVymeaMQcn3bfh06ak52DQOgZN53Y81zA9/EwHnOhtBtMjTHEmXUhX0REZAVRABURWeGiwGMomN/yJc8tzcQFwnrHfbTijCDyiPIcC/gmLDrnWpJimmy3C25Xmhs6qbvtFWU8a902MVHgFcHUVQn7KiFre8t4xW4sxrhqZ5zmtOKMRuzG0UlyOqml0YlJivDY7VrkqomGwHf7m3qeq6wGRQffwPMgdIEWLFiDNcVaSvvkbrK5dV2Fc+u+Nrg9VfPcVVmtNeTWkhfh1lropDlx6qYfd1K3prWdZsRJTpxlxKmbShynmXvuHJIiCOdF5TbLcheOrXutvKj8Fu/WBWgL080O5Ug/pkVEZGXRTzYRkYuM5xl6SgE9pfkfAWmWuyplnM4F0zjN8T2DMZbQN3PnWyDLLJnNi0CWk+f2lOm4aZ7T7BggdmkTF6q6Vc7Q7zZBcg2W1vSW2TR4atOgbkCNs5xOktGIXdUySS2dtKh45ilZnhcvYYrw7DoJ+74h9NzUV88wf7voAOyqkL4bdhFyLd0g2h3DwrLn05Uju0HSBU07n5ux1p7yvgyAMeR5DtaQY4tAbufeB8D7b3uEONE8XBERWVkUQEVEhMD36K969FfDuWOdNKPZyWgmGc1OSiN201UN4AXg5a7COZetLK6qt2BrFluEqW4IzawlizPaRXUVmAuA3c613e1WSqFP5HuUQo+ecshQT4nINwvWSFqy3AVe9zFfqW0nmZs2m9u5CmTWKaboPmkjGDcG33MNiLqNiDwDvuf2J/WMxTduerDnLQjKphtLDQFguyHVGjCnv063xuqqsda6bsQUwdN2z7BgjaUchgRPeg4REZELmwKoiIickVvf6TO44FiWW7d2Ms5oFh+tJCVOXQU0MJBbQxh51ErgwdxaUIObjpqkLih6xTTabljLLWRFx95Gx5LmneJVXajzDHPTb4OiqVE59IsGRx5RYCiHAUO1yE3N9VzcM8Zgc0tq3XTiLCuaGxVrRt0046yothbHiumxcZaTJ921pG76cg5zobGbHLsVTOYiZvHGjQuy3VMNRYBesNaTuftOebtEvqcAKiIiK44CqIiIPGu+Z+gth/SWw1OOJ1nuwmgxTdYF04y06NzjGUNu82KrFlcN7U43zYtprnkREOM0ww88ItxU2e560e7ayNxSdPq15LYDRZhzzYfANwbf8/C9hYF1Pqy6Sqv73O1UG3im+PA4ZU6uMXPNjFzF9dSvu+tHu8e778N2jy9Y55nPrfss1nvaU6f+dnWPVaOg+A6IiIisHAqgIiJy1kLfo7/i0V85NZjGaU4ryegkGe3ENexpF51n87mMZ7DW4vsuNFYj38Wxop+Qze3cFN7u3qG+NUQBC/YYnTtlrlppcGtVO0lOmuUuGHbPYz7oecW0Wt+41+/uA9qdktsNqQun5y78OjIenregi615cldbzxhMcQ5FwyNjIC+aJbk3ML9w1BjDv+0dJVnQ9ElERGQlUAAVEZFFExVbtHBaMAW3xrSduAZDrW5ALdZvzvX+6SZJXCgrhT6lcH66qrXdEAcLjpLl1nWuzXFB1fhzKzDnPzM/rbboQktRYU27nWrz+em38+s051+7+3rGnLavpwGP+X0/u1Ns588xxTl03xzdd9ANyEmWu4qsiIjICqIAKiIiy6K7xvRM4dTtvZkV2524kLrw61P35OyGN1uEOg8/gGjhtNbuZ7swOM6zdn5qrQuApuhWuzDWumrpXGOl4lm7a1ztgqm1cy2Bi//tdsV1XXK704nt3H0Lx9F9TDn0njTVWURE5EKnACoiIuedbuW09ynuTzLXZTfJ3DYt7ms7dzzOcpLUzk3zPZWdr0Iu2HYlDMz8FizWTfu1uetISzc0Lkykp2/L4i94BbtgT0/mI3K3Wrrw0d01ru5rNx0ZA5sGa/RX9GNaRERWFv1kExGRC0638+0zmduiJbV0MtcUKc0sSV5s0VKsKXUdcfO5hkiOPeUTzHe0tVB08C1qpd35tXPnn7pNi13wVOa0vUYp7nc9mSzWuiA8VIsYqkXP+XsjIiJyPlMAFRGRFWsuqEYAzzyd1a0dzeeCane7loXrQbudbTPrzlnYEbc7xbbbCfepPdP2KpbhvjKDNU3BFRGRlUUBVEREpOA63PqUzsFPx/ltY04NpZbTQmq3ARKnrhUFqJ2LgYiIiJxH9JNNRERkEXQ74HqnrxUVERG5iKm/u4iIiIiIiCwJBVARERERERFZEgqgIiIiIiIisiQUQEVERERERGRJGPv0feJXBGPMKHBwucexAq0GxpZ7ECLPga5ZudDompULja5ZudDoml08W621a04/eFEEUFkcxpg7rbV7lnscIs+Wrlm50OialQuNrlm50OiaXXqagisiIiIiIiJLQgFUREREREREloQCqJyNDy/3AESeI12zcqHRNSsXGl2zcqHRNbvEtAZUREREREREloQqoCIiIiIiIrIkFEBFRERERERkSSiAytMyxgwZY75ojNlbfB58ivN+wBjzqDFmnzHmt85w/68bY6wxZvXij1ouZmd7zRpj/qcx5hFjzH3GmH8yxgws3ejlYvIs/t00xpj/p7j/PmPMC5/tY0UWw/O9Zo0xm40xXzHGPGyMedAY88tLP3q5GJ3Nv7PF/b4x5m5jzKeWbtQrnwKoPJPfAm6z1u4Cbitun8IY4wN/DPwgcAXwE8aYKxbcvxl4FXBoSUYsF7uzvWa/CFxprb0aeAz4j0syarmoPNO/m4UfBHYVH7cCH3wOjxU5p87mmgVS4NestZcDLwF+UdesLLazvGa7fhl4eJGHetFRAJVn8nrgz4uv/xz4kTOccz2wz1r7hLU2Bv6meFzX/w38BqCOV7IUzuqatdZ+wVqbFufdDmxa5PHKxemZ/t2kuP0X1rkdGDDGrH+WjxU51573NWutPW6t/S6AtXYW9wv9xqUcvFyUzubfWYwxm4AfAj6ylIO+GCiAyjMZttYeByg+rz3DORuBwwtuHymOYYx5HXDUWnvvYg9UpHBW1+xpfg747Dkfocizuwaf6pxne/2KnEtnc83OMcZsA64Dvn3ORyhyqrO9Zt+PK6DkizXAi1Ww3AOQ5WeM+RKw7gx3/adn+xRnOGaNMdXiOV79fMcmciaLdc2e9hr/CTdt7K+e2+hEnpVnvAaf5pxn81iRc+1srll3pzE9wD8Av2KtnTmHYxM5k+d9zRpjfhgYsdbeZYx5xTkf2UVOAVSw1n7fU91njDnZnT5TTEkYOcNpR4DNC25vAo4BO4HtwL3GmO7x7xpjrrfWnjhnb0AuOot4zXaf423ADwO3WG2WLIvjaa/BZzgnehaPFTnXzuaaxRgT4sLnX1lr/3ERxynSdTbX7JuA1xljXgOUgT5jzMettT+9iOO9aGgKrjyTTwJvK75+G/DPZzjnO8AuY8x2Y0wEvAX4pLX2fmvtWmvtNmvtNtx/5C9U+JRF9ryvWXAd84DfBF5nrW0uwXjl4vSU1+ACnwTeWnRpfAkwXUwrfzaPFTnXnvc1a9xfof8UeNha+4dLO2y5iD3va9Za+x+ttZuK31/fAnxZ4fPcUQVUnsl7gb8zxrwD18X2zQDGmA3AR6y1r7HWpsaYdwOfB3zgo9baB5dtxHKxO9tr9o+AEvDFonJ/u7X2XUv9JmRle6pr0BjzruL+DwGfAV4D7AOawNuf7rHL8DbkInI21yxwE/AzwP3GmHuKY79trf3MUr4Hubic5TUri8hodpmIiIiIiIgsBU3BFRERERERkSWhACoiIiIiIiJLQgFUREREREREloQCqIiIiIiIiCwJBVARERERERFZEgqgIiIiIiIisiQUQEVERM4BY8wqY8w9xccJY8zR4uu6MeZPFuH1PmaM2d/d0+4M99eLzzu74zjXYxAREXmutA+oiIjIOWaM+a9A3Vr7B4v4Gh8DPmWt/funuL9ure15qtsiIiLLIVjuAYiIiKxkxphXAL9urf3hIphuB9YDlwLvAV4C/CBwFHittTYxxrwI+EOgBxgDftZae/wZXmc78Ne4n+2fW5x3IyIicnY0BVdERGRp7QR+CHg98HHgK9baq4AW8EPGmBD4X8CbrLUvAj4K/O6zeN4PAB+01r4YOLEoIxcRETlLqoCKiIgsrc8WVc77AZ/5auX9wDZgN3Al8EVjDMU5T1v9LNwE/Gjx9V8C7zuHYxYRETknFEBFRESWVgfAWpsbYxI734whx/1cNsCD1tobn8dzq7GDiIic1zQFV0RE5PzyKLDGGHMjgDEmNMa84Fk87hvAW4qvf2qxBiciInI2FEBFRETOI9baGHgT8D5jzL3APcBLn8VDfxn4RWPMd4D+RRyiiIjI86ZtWERERC5Az7QNyxnO1zYsIiKy7FQBFRERuTBNA//NGPOupzvJGLPTGHMPcHJphiUiIvLUVAEVERERERGRJaEKqIiIiIiIiCwJBVARERERERFZEgqgIiIiIiIisiQUQEVERERERGRJ/P9D6Te9bi99gAAAAABJRU5ErkJggg==\n",
"text/plain": [
""
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"flux1 = tm.evaluate(k=0.1, ldc=[None], t0=0.0, p=1.0, a=4.2, i=0.5*pi, e=0.0, w=0.0)\n",
"ax = plot_lc(time, flux1);\n",
"\n",
"for i in range(10):\n",
" flux1 = tm.evaluate(k=0.1, ldc=[None], t0=0.0, p=1.0, a=4.2, i=0.5*pi, e=0.0, w=0.0)\n",
" ax = plot_lc(time, flux1, ax=ax, c='C0', alpha=0.25); "
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"---\n",
"\n",
"©2020 Hannu Parviainen"
]
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.7.7"
}
},
"nbformat": 4,
"nbformat_minor": 4
}