{ "cells": [ { "cell_type": "markdown", "id": "ab7a3f04f2cbbdc4", "metadata": {}, "source": [ "# Usage\n", "\n", "**Authors:** Judith Korth and Hannu Parviainen
\n", "**Edited:** 4 March 2025" ] }, { "cell_type": "code", "execution_count": 1, "id": "72228aed66a65c83", "metadata": { "ExecuteTime": { "end_time": "2025-03-04T13:41:58.856609Z", "start_time": "2025-03-04T13:41:58.850256Z" } }, "outputs": [], "source": [ "%matplotlib inline" ] }, { "cell_type": "code", "execution_count": 2, "id": "initial_id", "metadata": { "ExecuteTime": { "end_time": "2025-03-04T13:43:46.766800Z", "start_time": "2025-03-04T13:43:46.764964Z" } }, "outputs": [], "source": [ "import matplotlib.pyplot as plt\n", "\n", "from numpy import pi\n", "from cheopslcs import LCSim" ] }, { "cell_type": "markdown", "id": "cb40c22a-be7f-485a-b59a-c99aeecddcbc", "metadata": {}, "source": [ "## Simulator setup\n", "\n", "The light curve simulation is done using `cheopslcs.LCSim` class, which is initialised by giving the observing window width in hours, exposure time in seconds, and the white noise standard deviation. Optionally, we can also give the limb darkening model used by the transit model (one of: `constant`, `linear`, `quadratic`, `nonlinear`, `general`, `power-2`, and `power-2-pm`, see the [PyTransit documentation](https://github.com/hpparvi/PyTransit/blob/9a39e7162b93e9f61f8a79409606b8c1649fc750/doc/source/notebooks/models/roadrunner/roadrunner_model_example_1.ipynb) for details)." ] }, { "cell_type": "code", "execution_count": 3, "id": "4408f366-4eeb-476c-a309-adc82d9e08be", "metadata": {}, "outputs": [], "source": [ "lcs = LCSim(window_width=5.0, exp_time=60, white_noise=1e-4)" ] }, { "cell_type": "markdown", "id": "d50080b7d2326a59", "metadata": {}, "source": [ "## How to simulate a transit\n", "\n", "After the simulator has been set up, we can create light curves by calling the simulator as a function. A simple transit can be created by giving the \n", "\n", "- planet-star radius ratio,\n", "- zero epoch,\n", "- orbital period in days,\n", "- semi-major axis divided by the stellar radius,\n", "- impact parameter, and\n", "- limb darkening coefficients specific to the used limb-darkening model\n", "\n", "as arguments, and the call will return two NumPy arrays containing mid-exposure times and fluxes." ] }, { "cell_type": "code", "execution_count": 4, "id": "c97f1a43534cc798", "metadata": { "ExecuteTime": { "end_time": "2025-03-04T14:06:12.665109Z", "start_time": "2025-03-04T14:06:12.544227Z" } }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAA/MAAAGbCAYAAACIxMC9AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8hTgPZAAAACXBIWXMAAA9hAAAPYQGoP6dpAABYk0lEQVR4nO3de3xU9Z3/8fdkSMhAyHAJhiCEpBQJBeQSKZcoSrsFQSyB9me8LCt0l5/8hNbIbpEotMhWQuvKatWgoFiUKuyWqO2K2KzWCw0lEKOCKCgQQ+NQOpROKgkhJOf3B51x7pncmJzJ6/l4jDJnvufM90zO7fO9WgzDMAQAAAAAAEwjLtoZAAAAAAAALUMwDwAAAACAyRDMAwAAAABgMgTzAAAAAACYDME8AAAAAAAmQzAPAAAAAIDJEMwDAAAAAGAy3aKdAbNqamrS559/rl69eslisUQ7OwAAAACAGGAYhv72t79p4MCBiosLXf9OMN9Kn3/+uQYPHhztbAAAAAAAYtCJEyc0aNCgkJ8TzLdSr169JF38gZOTk6OcGwAAAABALKipqdHgwYM9MWcoBPOt5G5an5ycTDAPAAAAAGhXzXXnZgA8AAAAAABMhmAeAAAAAACTIZgHAAAAAMBkCOYBAAAAADAZgnkAAAAAAEyGYB4AAAAAAJMhmAcAAAAAwGQI5gEAAAAAMBmCeQAAAAAATIZgHgAAAAAAkyGYBy4Bh6tOpUedcrjqop0VAAAAADGgW7QzAMS67fuqVFB8QE2GFGeRCueNVt6E9Ghnq0txuOp03HlWmSk9lWa3hVwGAAB8RXK/5J7aefG3iW0E810IJ/Ol53DVeQJ5SWoypHuLD2rqFf35G1wiwQpTJLWogIVzBwDQFUVSIUGlRefF3yb2Ecx3EbFwMpsxoDruPOsJ5N0aDUOVzlrT7ENnFOmxEKwwpWDHAcmiiAtYIj13zHh8AgDMobW14y29N3mnlxS2QsLhqtP+yr9QaRGhS/2cEIsVSjxrBSKY7wKidTI3d8K15IQ0a2FEZkpPxXkFjpJktViUkdIjepkyMYerTpt3H9fTu49HdCwEK0xpkqQIC1giPXfMenwCADo//3vMPddnafQgu8/zU3u0QvPehkXSDaMHBK2QeOUDhxoNQz999eOAz91pqLTw5f/bLromUwuvzmxRAUvPBKvOnm+MOJANVaH0ygcO3XBlmidNZwyMg8UIPGsFF9Vg/u2339aDDz6o8vJyORwOvfjii8rNzQ27zltvvaVly5bpww8/1MCBA7V8+XItXrzYJ82OHTu0atUqHT16VEOHDtUDDzyguXPn+qQpKirSgw8+KIfDoZEjR+rhhx/WNddc0967GHUOV53+54PPI6odDhVc+5fSRlLqG+6Eay4g899eqNrVrAG9NGZwn2bzGWme2/o7h/q+wnmjdW/xQTUahqwWi9bOGxXxd1IC+aXt+6q0YscBnzg8XMGUw1Wn01/UBxSmBBMn6fTZes8Ahe7fPJKWFR1RWMbfHQDMI5Lnp7bcD/zvMYWvfizpy+enqVf0D0izYscBSV+WXTd3b/L/HkPS/xw4GTRPP3nlo7B5ptLiS8FaLxiSNr5zXE/tPu5TMCMFPrseqHYFFJpEGsgGq1CSLv79Hvj739BQ6AKiaAkWQwQ7xoMdz13x+SmqwfzZs2c1ZswYLVy4UN/5zneaTX/8+HHNmjVLixYt0tatW/X73/9ed955p/r37+9Zf8+ePcrLy9O///u/a+7cuXrxxRd10003affu3Zo4caIkafv27crPz1dRUZFycnL05JNPaubMmTp06JDS02OnhMf7ZPDnHbyk2W0hg2//kkTJ98R3flEfEJSHO+HePvLnsAHZ20f+HFD6HBcXeCFqkpRbVKp1zeTTP88tKalurrmaFHih9f9ud8nr7hXTVOmsVUZKj4ibx1EC+SX3Q0awmNy7lDlYYZIlyDreLLr491r6fEXQYzxcy4pwhWX+efLfn1A3G/7uANB5RXqvjqQ2PZLAo/yzMyELpN2VG0u/+dWANJHeL92CFV63RksrLWJZuOdwybdgJtSza6j1wlVkhKpQ8ub/HO5fQOT/3NGSAqu2DJgYqoLkkVvGhq1cCVZRGKqAItYCfothGO1w6radxWJptmb+nnvu0a9//Wt99NGXJYKLFy/W+++/rz179kiS8vLyVFNTo1dffdWT5vrrr1efPn30wgsvSJImTpyo8ePHa8OGDZ40I0aMUG5urgoLCyPKb01Njex2u1wul5KTk1uyq5eEw1WnnHVvBL2A+F8wbvn6YG3bd8K31E/S/bkj9eOXP2zRxd1qseiRW8Zq6fMVAZ89dss4/WBbRcjtNfd5qO8rvnOy5haVRrSe5e//MfyCs90rpjXbbFpS0AKD5rS0aZv7AvTTXR8HBJH++Ywl4S6upUedunXT3rDrexcm+R/77oDd313f+Koe/d2nIY8dq8Wi5TOH62evHvZpWeH/oNZcnrz/9s21WvHPu/ffPdZuQABgBu5rr3/hfbDnp+aeS1pSsRCsRVp7CBbohHtulKT5E9P13N6q0NuU9Oit4zR+SB/uTwr/HN5eXlg0SZOH9vO8D1WA1DPBqn2VZ5ptUeEWJ+nnt47T4D42nT3f6HPce2/Xf3lzx3Ww88g7jbuC5IFXPg7IU7AYwf18FKyi0Gd/whSydeYKk0hjTVP1md+zZ4+mT5/us2zGjBl6+umn1dDQoPj4eO3Zs0d33313QJqHH35YknT+/HmVl5drxYoVPmmmT5+u0tLSkN9dX1+v+vp6z/uampo27k3HClXC6h+8GJKeLzsRkK5J0qqXPmzx9zYaho6c/FvQ2kxZQjd3bu7zcN/33/v/GPF6huc/vtvwLtkLNpiL/6BpLclmsAFjmutGsO7VjwO+I5b7oDVXgxGqqZi3cKW3hiRLkEKcYQN6hd2mO4AvvnOyas83eWrkf/N+dbOBvHeepl7RX5KaHSgoXLN+/1YrnfkGBACxIlTBbajnp0bD0L7K5mvTmxuINVyLtLbyr4l133PvmZmln+78+OLYMl6sFou+e9Ug/bKsKuh+uQu6b7hyYAfk1nxCtdprT/7dA4M9X/j/jSPpdihdjAGCVcr5b9d/ebCuHQU7Dqhn926q/mtd0DEW3Mf+X+saQo7BECdJFumemVkBlSuSmj1P3Pk439jkU0nZHt0iOwNTBfMnT55Uamqqz7LU1FRduHBBTqdTaWlpIdOcPHmx34/T6VRjY2PYNMEUFhbq/vvvb6c96XihBl5rLnhpDz9/41NZ9GXwZLVYtPz64TIMI+iFJM4irZ03StlD+jR7oQlWw7o1TElxJNzNpsPVtAYbNK0lwgVkg/v2iKh5XKz2QYukP2DehHSfpmJxkmZdOUD/84HvOdtoGNLff1v/Yz9YDXskx9xPXvnIk4+qv5wNG8QHq7loNAw9s7tST+0+FnQ970KangnWgEKHOEmf/vlvMXkDAoDOzP/+FAmrxaIJGeHvLZEMxBqqUubWr6dr+74TAc2m3YK1OAvVOk0KEvTNzNLpL8577lnu++WYwX0CxgBafv1wXTmod9guhF1FsFpnf3GSbpmYrhdCFIpEKlT3wFCaDOlnuw77BMPBWqq2VbBNhSsYcGs0DK3bGViJJfnuq7tQwn3MSYq40CRUJWUsVJSZKpiXLjbH9+buJeC9PFga/2WRpPFWUFCgZcuWed7X1NRo8ODBLcv8JRRq4LVIgpdQWnLiG5LiDOmxW8fpj16lcd5Bfpykf5maqYU5X47mGapvjxQYkEWSz3A3MPc2vUv2Qv0u7lLB5n63UN9ntVjUIyEuaI3//bkjm/2bxHIftHD99LyD1rwJ6Zp6RX/P2AOStPPAyYDuIaFKb/MmpOvbYwYGjF3gfcyFOsaD1aT4C1VzESeFDOTd63kXJnl/t/t4itUbEAB0Zi3tRx4q8PUX7JnCv6Y12ACuVotF3//mV/X9b35V5ZVngjY7vnliugb2sQUE3f5d94JxB327V0zTwqszAu6X/vdh7j8XRdLtzvtZZOk3vqpKZ60+qP5r0ODa+9/ehSY9EuJ04i91Pn/3SA/PRsPQlZf39hnDSVJAPqIl2Df7F0x5H5/elWNtEQsVZaYK5gcMGBBQe37q1Cl169ZN/fr1C5vGXROfkpIiq9UaNk0w3bt3V/fu3dtjNy6ZUBfdwnmjVbDjQEAzqjhJa3JH6kd+/eS9+0FJgSd+qFrSJl28EHmXUHoH+d79qtwlmlOv6O+50Hh/h39A9soHjqB9f1bdMEKz/j7dhnu/T9WcU25RaUBtp3ffrtKjzrDBljvgDxb0+ZdOS/KpiXWvf/Z8Y2DzIl0M1MIN0ubep1i9aTbXhN47aHW/3PwD8VClt+51/NeXAs8TSUGPr3CtM0I9wMVJmjV6QMgRgd3Hjn/zOKn5wrNYuAEBQGcWSRcvKXh/ce97i/8z079ck6l+vbr7BHLBalr9Wzl6F+rPHmPT2fMXgs6WE+z5r3eP+JCFC97c99zJQ/sFfe4Idh/tavwHQ24uqPR/jnP/hpOH9vOpZJAU9N/ev3ewZ8lIuJ8Z/P9+3vkIVkAULaFaEzcahsorz4Qd4PtfpmaqX1L3ZgsoYqWizFTB/OTJk/Wb3/zGZ9lvf/tbXXXVVYqPj/ekKSkp8ek3/9vf/lZTpkyRJCUkJCg7O1slJSU+09WVlJRozpw5l2AvLq1wwUuwgDNvQrrirXEBNwfvflChLkD+taSh+sE3Serbs3uz09j5f4f3RfCGK9O0dudHAd/nf7F0/39dkFYK7n0KNYWZ++Y86O8DgGSm9AxaohmsdPreG0YElGo7XHUhHwrcN233v0PtUywKN9qqFDjzgjf3sex/A/IuvY3kt/M/T4IdX6FqUkI9wLnPr2CBvHu9P4boQyYFH9/BLVZuQADQmQVr5bj8+uFBm6F7P1N4j4vj/czkvi9sfOe4p9D58t62kDWtoSpA3MLVlPvf10IVLvijoDg8/9lybhg9oNka+XDPccGC62D/dmuugMn7+cK/QizcM0Oa3RZQQOS9D+5KK/+KtuXXDw84hqXwFRLu9axxlqAVc+G634YbY8u/0CRUAUWsDdYY1dHsv/jiC3366aeSpHHjxmn9+vWaNm2a+vbtq/T0dBUUFKi6ulrPPvuspItT040aNUp33HGHFi1apD179mjx4sV64YUXPFPTlZaWaurUqXrggQc0Z84cvfzyy1q5cmXA1HTz58/XE088ocmTJ2vjxo3atGmTPvzwQw0ZMiSivHf20ewj5XDVBb0JhFoezvZ9VQEBc7CRxf1H5w73eUu/L9yAYMH2KWAKM78ScCmyaewi5Z3nYP7vNV/R07uPR7xPscT99wlodibfqQWD/R6hRrv3H+m1JYIdX5IiOubCjWIb7tzwFmnhAQCgYwV7fmjumSLS2UpCzQDkrS33sub2KVQrSARq6Qj1HfV7ej+fBHt29T7mWtMlwr1ej4Q4z+C/zcUIkTwz+Y+xEOz3jJP04pIpGjO4T8jtNhdbhPu9zHSMRxprRjWYf/PNNzVt2rSA5bfffrt+8YtfaMGCBaqsrNSbb77p+eytt97S3XffrQ8//FADBw7UPffco8WLF/us/6tf/UorV67UsWPHNHToUD3wwAOaN2+eT5qioiL97Gc/k8Ph0KhRo/Sf//mfmjp1asR5j5Vgvr1FeoK7T6K2BmGtvVC51w12EfHuUtDagobmvjdUX7fdKy6eD129P1pzv1GweVU76m8VyQOcv1DHtXepcbjp9lpaeAAAiK7m7kOhrvnNTct7KaalbcuzVFcSyTS5l2pgQO+/mdQ5nhtb88wUSaDd0tgi0rx1dqaYmu66665TuLKEX/ziFwHLrr32Wr377rtht/vd735X3/3ud8OmufPOO3XnnXdGlE9ELpI+yd6fhxp1P9ImXm3puxVsYBvvLgDB+tG3x6BjwZoy+TeBMsuFpqOk2W3qmxR6mjb/3yfUgI9t/R2DHV+RHHOhjmvv5l/B0gSreQ917jDnPAB0nJZeY8NNLZpmDz7FqtViUXZGn5ADsV6qLlX0g49Mc03cL+UYR+Ga50dLa56ZIhlUsaWxRaR5ixWm6jMP8wp1EnVUEBaJ5goS2lrQ0BxGhQ2vpb9/Z/o9IzmuQ6Xxn6c32LkTriknAKBtWnONbe6eFe6+EGwg1s5wL4Mv998w2CDSXWGMo47S2kA7lgP0lohqM3szo5l9+4pW85fmmumYtZ9NrDD77x/Jcd3SY7+juhQAANp2jW1tk2GYi8NVF3IQaaC9mKLPvJkRzMeO5m6s3Hiji9/fV0cM9gcAuCiaY/nAXPhboyOZos880Bk010yHZjzRxe/vq6O7fwBAVxbNsXxgLvyt0RnERTsDAIDIufvsWS0XJ+5jznkAaD9cYwGYCc3sW4lm9gCiieZ9ANBxuMYCiCaa2QNADKN5HwB0HK6xAMyAZvYAAAAAAJgMwTwAAAAAACZDMA8AAAAAgMkQzAMAAAAAYDIE8wAAAAAAmAzBPAAAAAAAJkMwDwAAAACAyRDMAwAAAABgMgTzAAAAAACYDME8AAAAAAAmQzAPAAAAAIDJEMwDgMk5XHUqPeqUw1UX7awAgGlxLQVgNt2inQEAQOtt31elguIDajKkOItUOG+08iakRztbAGAqXEsBmBE18wBgUg5XnefhU5KaDOne4oPUKgFAC3AtBWBWBPMAYFLHnWc9D59ujYahSmdtdDIEACbEtRSAWRHMA4BJZab0VJzFd5nVYlFGSo/oZAgATIhrKQCzIpgHAJNKs9tUOG+0rJaLT6FWi0Vr541Smt0W5ZwBgHlwLQVgVhbDMIzmk8FfTU2N7Ha7XC6XkpOTo50dAF2Yw1WnSmetMlJ68PAJAK3EtRRAZxFprMlo9gBgcml2Gw+eANBGXEsBmA3N7AEAAAAAMBmCeQAAAAAATIZgHgBiiMNVp9KjTuZHBgAAiHH0mQeAGLF9X5UKig+oyZDiLFLhvNHKm5Ae7WwBAACgA1AzDwAxwOGq8wTyktRkSPcWH6SGHgDCoDUTADOjZh4AYsBx51lPIO/WaBiqdNYyOjMABEFrJgBmR808AMSAzJSeirP4LrNaLMpI6RGdDAFAJ0ZrJgCxgGAeAGJAmt2mwnmjZbVcjOitFovWzhtFrTwABBGuNRMAmAXN7AEgRuRNSNfUK/qr0lmrjJQeBPIAEIK7NZN3QE9rJgBmQ808AMSQNLtNk4f2I5AHgDBozQQgFlAzDwAAgC6H1kwAzI5gHgAAAF1Smt1GEA/AtGhmDwAAAACAyRDMAwAAAABgMgTzAAAAAACYDME8AAAAAAAmQzAPAAAAAIDJRD2YLyoqUmZmphITE5Wdna133nknbPrHH39cI0aMkM1m0/Dhw/Xss8/6fN7Q0KA1a9Zo6NChSkxM1JgxY7Rr1y6fNBcuXNDKlSuVmZkpm82mr3zlK1qzZo2ampraff8AAAAAAGhvUZ2abvv27crPz1dRUZFycnL05JNPaubMmTp06JDS09MD0m/YsEEFBQXatGmTJkyYoLKyMi1atEh9+vTRjTfeKElauXKltm7dqk2bNikrK0uvvfaa5s6dq9LSUo0bN06S9NOf/lRPPPGEtmzZopEjR2r//v1auHCh7Ha77rrrrkv6GwAAAAAA0FIWwzCMaH35xIkTNX78eG3YsMGzbMSIEcrNzVVhYWFA+ilTpignJ0cPPvigZ1l+fr7279+v3bt3S5IGDhyo++67T0uWLPGkyc3NVVJSkrZu3SpJmj17tlJTU/X000970nznO99Rjx499Nxzz0WU95qaGtntdrlcLiUnJ7dsxwEAAAAACCLSWDNqzezPnz+v8vJyTZ8+3Wf59OnTVVpaGnSd+vp6JSYm+iyz2WwqKytTQ0ND2DTuYF+Srr76ar3++us6cuSIJOn999/X7t27NWvWrJD5ra+vV01Njc8LAAAAAIBoiFow73Q61djYqNTUVJ/lqampOnnyZNB1ZsyYoaeeekrl5eUyDEP79+/X5s2b1dDQIKfT6Umzfv16ffLJJ2pqalJJSYlefvllORwOz3buuece3XLLLcrKylJ8fLzGjRun/Px83XLLLSHzW1hYKLvd7nkNHjy4HX4FAAAAAABaLuoD4FksFp/3hmEELHNbtWqVZs6cqUmTJik+Pl5z5szRggULJElWq1WS9Mgjj2jYsGHKyspSQkKCli5dqoULF3o+ly721d+6dauef/55vfvuu9qyZYv+4z/+Q1u2bAmZz4KCArlcLs/rxIkTbdxzAAAAAABaJ2rBfEpKiqxWa0At/KlTpwJq691sNps2b96s2tpaVVZWqqqqShkZGerVq5dSUlIkSf3799dLL72ks2fP6rPPPtPHH3+spKQkZWZmerbzwx/+UCtWrNDNN9+s0aNHa/78+br77ruD9tN36969u5KTk31eAAAAAABEQ9SC+YSEBGVnZ6ukpMRneUlJiaZMmRJ23fj4eA0aNEhWq1Xbtm3T7NmzFRfnuyuJiYm6/PLLdeHCBe3YsUNz5szxfFZbWxuQ3mq1MjUdAAAAAMAUojo13bJlyzR//nxdddVVmjx5sjZu3KiqqiotXrxY0sWm7dXV1Z655I8cOaKysjJNnDhRZ86c0fr163Xw4EGf5vF79+5VdXW1xo4dq+rqaq1evVpNTU1avny5J82NN96oBx54QOnp6Ro5cqQqKiq0fv16fe9737u0PwAAdBCHq07HnWeVmdJTaXZbtLMDAACAdhbVYD4vL0+nT5/WmjVr5HA4NGrUKO3cuVNDhgyRJDkcDlVVVXnSNzY26qGHHtLhw4cVHx+vadOmqbS0VBkZGZ40586d08qVK3Xs2DElJSVp1qxZeu6559S7d29PmkcffVSrVq3SnXfeqVOnTmngwIG644479KMf/ehS7ToAdJjt+6pUUHxATYYUZ5EK541W3oT0aGcLAAAA7Siq88ybGfPMA+iMHK465ax7Q01eV3arxaLdK6ZRQw8AAGACnX6eeQBA+zvuPOsTyEtSo2Go0lkbnQwBAACgQxDMA0AMyUzpqTi/2T2tFosyUnpEJ0MAAADoEATzABBD0uw2Fc4bLavlYkRvtVi0dt4omtgDAADEmKgOgAcAaH95E9I19Yr+qnTWKiOlB4E8AHhhtg8AsYJgHgBiUJrdxkMqAPhhtg8AsYRm9gAAAIh5DledJ5CXpCZDurf4oByuuuhmDABaiWAeAAAAMY/ZPgDEGoJ5AAAAxDxm+wAQawjmAQAAEPOY7QNArGEAPAAAAHQJzPYBIJYQzAMAAKDLYLYPALGCZvYAAAAAAJgMwTwAAAAAACZDMA8AAAAAgMkQzAMAAAAAYDIE8wAAAAAAmAzBPAAAAAAAJkMwDwAAAACAyRDMAwAAAABgMgTzAAAAAACYDME8AAAAAAAmQzAPAACAmOdw1an0qFMOV120swIA7aJbtDMAAAAAdKTt+6pUUHxATYYUZ5EK541W3oT0aGcLANqEmnkAAADELIerzhPIS1KTId1bfJAaegCmRzAPAACAmHXcedYTyLs1GoYqnbXRyRAAtBOCeQAAAMSszJSeirP4LrNaLMpI6RGdDAFAOyGYBwAAQMxKs9tUOG+0rJaLEb3VYtHaeaOUZrdFOWcA0DYMgAcAAICYljchXVOv6K9KZ60yUnoQyAOICQTzAAAAiHlpdhtBPICYQjN7AIhxzK0MAAAQe6iZB4AYxtzKAAAAsYmaeQCIUcytDAAAELsI5gEgRjG3MgAAQOwimAeAGMXcygAAALGLYB4AYhRzKwMAAMQuBsADgBjG3MoAAACxiWAeAGIccysDAADEHprZAwAAAABgMgTzAAAAAACYDME8AAAAAAAmQzAPAAAAAIDJEMwDAAAAAGAyBPMAAAAAAJhM1IP5oqIiZWZmKjExUdnZ2XrnnXfCpn/88cc1YsQI2Ww2DR8+XM8++6zP5w0NDVqzZo2GDh2qxMREjRkzRrt27QrYTnV1tf7xH/9R/fr1U48ePTR27FiVl5e3674BAAAAANARojrP/Pbt25Wfn6+ioiLl5OToySef1MyZM3Xo0CGlp6cHpN+wYYMKCgq0adMmTZgwQWVlZVq0aJH69OmjG2+8UZK0cuVKbd26VZs2bVJWVpZee+01zZ07V6WlpRo3bpwk6cyZM8rJydG0adP06quv6rLLLtPRo0fVu3fvS7n7AAAAAAC0isUwDCNaXz5x4kSNHz9eGzZs8CwbMWKEcnNzVVhYGJB+ypQpysnJ0YMPPuhZlp+fr/3792v37t2SpIEDB+q+++7TkiVLPGlyc3OVlJSkrVu3SpJWrFih3//+9822AginpqZGdrtdLpdLycnJrd4OAAAAAABukcaaUWtmf/78eZWXl2v69Ok+y6dPn67S0tKg69TX1ysxMdFnmc1mU1lZmRoaGsKmcQf7kvTrX/9aV111lf7P//k/uuyyyzRu3Dht2rQpbH7r6+tVU1Pj8wIAAAAAIBqiFsw7nU41NjYqNTXVZ3lqaqpOnjwZdJ0ZM2boqaeeUnl5uQzD0P79+7V582Y1NDTI6XR60qxfv16ffPKJmpqaVFJSopdfflkOh8OznWPHjmnDhg0aNmyYXnvtNS1evFg/+MEPAvrfeyssLJTdbve8Bg8e3A6/AgAAAAAALRf1AfAsFovPe8MwApa5rVq1SjNnztSkSZMUHx+vOXPmaMGCBZIkq9UqSXrkkUc0bNgwZWVlKSEhQUuXLtXChQs9n0tSU1OTxo8fr7Vr12rcuHG64447tGjRIp/m/v4KCgrkcrk8rxMnTrRxzwEAAAAAaJ2oBfMpKSmyWq0BtfCnTp0KqK13s9ls2rx5s2pra1VZWamqqiplZGSoV69eSklJkST1799fL730ks6ePavPPvtMH3/8sZKSkpSZmenZTlpamr72ta/5bHvEiBGqqqoKmd/u3bsrOTnZ5wUAAAAAQDRELZhPSEhQdna2SkpKfJaXlJRoypQpYdeNj4/XoEGDZLVatW3bNs2ePVtxcb67kpiYqMsvv1wXLlzQjh07NGfOHM9nOTk5Onz4sE/6I0eOaMiQIW3cKwAAAAAAOl5Up6ZbtmyZ5s+fr6uuukqTJ0/Wxo0bVVVVpcWLF0u62LS9urra05f9yJEjKisr08SJE3XmzBmtX79eBw8e1JYtWzzb3Lt3r6qrqzV27FhVV1dr9erVampq0vLlyz1p7r77bk2ZMkVr167VTTfdpLKyMm3cuFEbN268tD8AAAAAAACtENVgPi8vT6dPn9aaNWvkcDg0atQo7dy501ND7nA4fJq+NzY26qGHHtLhw4cVHx+vadOmqbS0VBkZGZ40586d08qVK3Xs2DElJSVp1qxZeu6553zmkJ8wYYJefPFFFRQUaM2aNcrMzNTDDz+s22677VLtOgAAAAAArRbVeebNjHnmAQAAAADtrdPPMw8AuPQcrjqVHnXK4aqLdlYAAADQBlFtZg8AuHS276tSQfEBNRlSnEUqnDdaeRPSo50tAAAAtAI18wDQBThcdZ5AXpKaDOne4oPU0AOIabRGAhDLqJkHgC7guPOsJ5B3azQMVTprlWa3RSdTANCBaI0EINZRMw8AXUBmSk/FWXyXWS0WZaT0iE6GAKAD0RoJQFdAMA8AXUCa3abCeaNltVyM6K0Wi9bOG0WtPICYFK41EgDECprZA0AXkTchXVOv6K9KZ60yUnoQyAOIWe7WSN4BPa2RAMQaauYBoAtJs9s0eWg/AnkAMY3WSAC6AmrmAQAAEHNojQQg1hHMAwAAICal2W0E8QBiFs3sAQAAAAAwGYJ5AAAAAABMhmAeAAAAAACTIZgHAAAAAMBkCOYBAAAAADAZgnkAAAAAAEyGYB4AAAAAAJMhmAcAAAAAwGQI5gEAAAAAMBmCeQAAAAAATIZgHgAAAAAAkyGYBwAAAADAZAjmAQAAAAAwGYJ5AAAAAABMhmAeAAAAMcPhqlPpUaccrrpoZwUAOlS3aGcAAAAAaA/b91WpoPiAmgwpziIVzhutvAnp0c4WAHQIauYBAABgeg5XnSeQl6QmQ7q3+CA19ABiVquC+bq60BdFh8PR6swAAC4dmqICiCXHnWc9gbxbo2Go0lkbnQwBQAdrVTA/btw4vfvuuwHLf/WrX+nKK69sc6YAAB1r+74q5ax7Q7du2qucdW9o+76qaGcJANokM6Wn4iy+y6wWizJSekQnQwDQwVoVzH/rW9/SlClTtG7dOhmGoS+++EILFizQ7bffrh/96EftnUcAQDuiKSqAWJRmt6lw3mhZLRcjeqvForXzRinNbotyzgCgY7RqALxHH31UN9xwgxYuXKhXXnlFn3/+uZKTk7Vv3z597Wtfa+88AgDaUbimqDz0AjCzvAnpmnpFf1U6a5WR0oNrGoCY1urR7KdPn6558+Zpw4YN6tatm37zm98QyAOACbibonoH9DRFBRAr0uw2gngAXUKrmtkfPXpUkydP1v/8z//otdde0/LlyzVnzhwtX75cDQ0N7Z1HAEA7oikqAACA+VkMwzCaT+arV69euuGGG/TEE0+od+/ekqTS0lL90z/9k3r16qWKior2zmenU1NTI7vdLpfLpeTk5GhnBwBazOGqoykqAABAJxNprNmqmvmioiJt27bNE8hL0pQpU1RRUaHx48e3ZpMAgEsszW7T5KH9COQBAABMqFU186BmHgAAAADQ/iKNNVs1AN6zzz4b8jOLxaL58+e3ZrMAAAAAACACraqZ79Onj8/7hoYG1dbWKiEhQT169NBf/vKXdstgZ0XNPAAAAACgvXVon/kzZ874vL744gsdPnxYV199tV544YVWZxoAAAAAADSvVcF8MMOGDdO6det01113tdcmAQAAAABAEO0WzEuS1WrV559/3p6bBAAAAAAAflo1AN6vf/1rn/eGYcjhcOixxx5TTk5Ou2QMAAAAAAAE16pgPjc31+e9xWJR//799Y1vfEMPPfRQe+QLAAAAAACE0Kpm9k1NTT6vxsZGnTx5Us8//7zS0tJatK2ioiJlZmYqMTFR2dnZeuedd8Kmf/zxxzVixAjZbDYNHz48YJq8hoYGrVmzRkOHDlViYqLGjBmjXbt2hdxeYWGhLBaL8vPzW5RvAAAAAACipV37zLfU9u3blZ+fr/vuu08VFRW65pprNHPmTFVVVQVNv2HDBhUUFGj16tX68MMPdf/992vJkiX6zW9+40mzcuVKPfnkk3r00Ud16NAhLV68WHPnzlVFRUXA9vbt26eNGzfqyiuv7LB9BAAAAACgvUU8z/yyZcsi3uj69esjSjdx4kSNHz9eGzZs8CwbMWKEcnNzVVhYGJB+ypQpysnJ0YMPPuhZlp+fr/3792v37t2SpIEDB+q+++7TkiVLPGlyc3OVlJSkrVu3epZ98cUXGj9+vIqKivSTn/xEY8eO1cMPPxzxPjLPPAAAAACgvUUaa0bcZz5YzXYwFoslonTnz59XeXm5VqxY4bN8+vTpKi0tDbpOfX29EhMTfZbZbDaVlZWpoaFB8fHxIdO4g323JUuW6IYbbtA//MM/6Cc/+Umz+a2vr1d9fb3nfU1NTbPrAAAAAADQESIO5n/3u9/p2LFjysjIUFxc21vnO51ONTY2KjU11Wd5amqqTp48GXSdGTNm6KmnnlJubq7Gjx+v8vJybd68WQ0NDXI6nUpLS9OMGTO0fv16TZ06VUOHDtXrr7+ul19+WY2NjZ7tbNu2Te+++6727dsXcX4LCwt1//33t25nAQAAAABoRy2KyocNGyan0+l5n5eXpz/96U9tyoB/Tb5hGCFr91etWqWZM2dq0qRJio+P15w5c7RgwQJJF+e4l6RHHnlEw4YNU1ZWlhISErR06VItXLjQ8/mJEyd01113aevWrQE1+OEUFBTI5XJ5XidOnGjF3gIAAAAA0HYtCub9u9fv3LlTZ8+ebdUXp6SkyGq1BtTCnzp1KqC23s1ms2nz5s2qra1VZWWlqqqqlJGRoV69eiklJUWS1L9/f7300ks6e/asPvvsM3388cdKSkpSZmamJKm8vFynTp1Sdna2unXrpm7duumtt97Sz3/+c3Xr1s2nBt9b9+7dlZyc7PMCgFjhcNWp9KhTDlddtLMCAACACLRqnvn2kJCQoOzsbJWUlGju3Lme5SUlJZozZ07YdePj4zVo0CBJF5vMz549O6Dpf2Jioi6//HI1NDRox44duummmyRJ3/zmN3XgwAGftAsXLlRWVpbuueceTw0+AHQV2/dVqaD4gJoMKc4iFc4brbwJ6dHOFgAAAMJoUTBvsVgCmsBHOuBdMMuWLdP8+fN11VVXafLkydq4caOqqqq0ePFiSRebtldXV3vmkj9y5IjKyso0ceJEnTlzRuvXr9fBgwe1ZcsWzzb37t2r6upqjR07VtXV1Vq9erWampq0fPlySVKvXr00atQon3z07NlT/fr1C1gOALHO4arzBPKS1GRI9xYf1NQr+ivNbotu5gAAABBSi4J5wzC0YMECde/eXZJ07tw5LV68WD179vRJV1xcHNH28vLydPr0aa1Zs0YOh0OjRo3Szp07NWTIEEmSw+HwmXO+sbFRDz30kA4fPqz4+HhNmzZNpaWlysjI8KQ5d+6cVq5cqWPHjikpKUmzZs3Sc889p969e7dkVwGgSzjuPOsJ5N0aDUOVzlqCeQAAgE4s4nnmpYvN0SPxzDPPtDpDZsE88wBigcNVp5x1b/gE9FaLRbtXTCOYB2AaDledjjvPKjOlJ9cuAKbX7vPMS10jSAeAriTNblPhvNG6t/igGg1DVotFa+eN4mEYgGkw7geArqpFNfP4EjXzAGKJw1WnSmetMlJ6EMgDMA1aFwGIRR1SMw8AiE1pdhsPvgBMh3E/AHRlLZpnHgAAAOgsMlN6Ks5vYiWrxaKMlB7RyRAAXEIE8wAAADAl97gf1r9Plcy4HwC6EprZAwAAwLTyJqRr6hX9GfcDQJdDMA8AAABTY9wPAF0RzewBAAAAADAZgnkAAAAAAEyGYB4AAAAAAJMhmAcAAAAAwGQI5gEAAAAAMBmCeQAAAAAATIZgHgAAAAAAkyGYBwAAAADAZAjmAQAAAAAwGYJ5AAAAAABMhmAeAAAAAACTIZgHAAAAAMBkCOYBAAAAADAZgnkAAAAAAEyGYB4AAACm43DVqfSoUw5XXbSzAgBR0S3aGQAAdB4OV52OO88qM6Wn0uy2aGcHAILavq9KBcUH1GRIcRapcN5o5U1Ij3a2AOCSIpgHAEji4RiAOThcdZ5rlSQ1GdK9xQc19Yr+FEIC6FJoZg8ACPlwTPNVAJ3NcedZz7XKrdEwVOmsjU6GACBKCOYBADwcAzCNzJSeirP4LrNaLMpI6RGdDAFAlBDMAwB4OAZgGml2mwrnjZbVcvGiZbVYtHbeKJrYA+hy6DMPAPA8HN9bfFCNhsHDMYBOLW9CuqZe0V+VzlplpPTgWgWgSyKYBwBI4uEYgLmk2W1cpwB0aQTzAAAPHo4BAADMgT7zAAAAAACYDME8AAAAAAAmQzAPAAAAAIDJEMwDAAAAAGAyBPMAAAAAAJgMwTwAAAAAACZDMA8AAAAAgMkQzAMAAAAAYDIE8wAAAAAAmAzBPAAAAAAAJkMwDwAIyuGqU+lRpxyuumhnBQAAAH66RTsDAIDOZ/u+KhUUH1CTIcVZpMJ5o5U3IT3a2QIAAMDfUTMPAPDhcNV5AnlJajKke4sPUkMPAADQiUQ9mC8qKlJmZqYSExOVnZ2td955J2z6xx9/XCNGjJDNZtPw4cP17LPP+nze0NCgNWvWaOjQoUpMTNSYMWO0a9cunzSFhYWaMGGCevXqpcsuu0y5ubk6fPhwu+8bAJjRcedZTyDv1mgYqnTWRidDAAAACBDVYH779u3Kz8/Xfffdp4qKCl1zzTWaOXOmqqqqgqbfsGGDCgoKtHr1an344Ye6//77tWTJEv3mN7/xpFm5cqWefPJJPfroozp06JAWL16suXPnqqKiwpPmrbfe0pIlS/SHP/xBJSUlunDhgqZPn66zZ892+D4DQGeXmdJTcRbfZVaLRRkpPaKTIQAAAASwGIZhNJ+sY0ycOFHjx4/Xhg0bPMtGjBih3NxcFRYWBqSfMmWKcnJy9OCDD3qW5efna//+/dq9e7ckaeDAgbrvvvu0ZMkST5rc3FwlJSVp69atQfPx5z//WZdddpneeustTZ06NaK819TUyG63y+VyKTk5OaJ1AMAstu+r0r3FB9VoGLJaLFo7bxR95gEAAC6BSGPNqA2Ad/78eZWXl2vFihU+y6dPn67S0tKg69TX1ysxMdFnmc1mU1lZmRoaGhQfHx8yjTvYD8blckmS+vbtGzJNfX296uvrPe9rampCpgUAs8ubkK6pV/RXpbNWGSk9lGa3RTtLAAAA8BK1ZvZOp1ONjY1KTU31WZ6amqqTJ08GXWfGjBl66qmnVF5eLsMwtH//fm3evFkNDQ1yOp2eNOvXr9cnn3yipqYmlZSU6OWXX5bD4Qi6TcMwtGzZMl199dUaNWpUyPwWFhbKbrd7XoMHD27lngOAOaTZbZo8tB+BPAAAQCcU9QHwLBbfjpmGYQQsc1u1apVmzpypSZMmKT4+XnPmzNGCBQskSVarVZL0yCOPaNiwYcrKylJCQoKWLl2qhQsXej73t3TpUn3wwQd64YUXwuazoKBALpfL8zpx4kQL9xQAAABt5XDVqfSokxk2AHR5UQvmU1JSZLVaA2rhT506FVBb72az2bR582bV1taqsrJSVVVVysjIUK9evZSSkiJJ6t+/v1566SWdPXtWn332mT7++GMlJSUpMzMzYHvf//739etf/1q/+93vNGjQoLD57d69u5KTk31eAAAAuHS276tSzro3dOumvcpZ94a27ws+aDIAdAVRC+YTEhKUnZ2tkpISn+UlJSWaMmVK2HXj4+M1aNAgWa1Wbdu2TbNnz1ZcnO+uJCYm6vLLL9eFCxe0Y8cOzZkzx/OZYRhaunSpiouL9cYbbwQN9AEAANB5OFx1Kig+4Jk6s8mQ7i0+SA09gC4ragPgSdKyZcs0f/58XXXVVZo8ebI2btyoqqoqLV68WNLFpu3V1dWeueSPHDmisrIyTZw4UWfOnNH69et18OBBbdmyxbPNvXv3qrq6WmPHjlV1dbVWr16tpqYmLV++3JNmyZIlev755/Xyyy+rV69entYBdrtdNht9QwEAADqb486znkDerdEwVOmsZWwPAF1SVIP5vLw8nT59WmvWrJHD4dCoUaO0c+dODRkyRJLkcDh85pxvbGzUQw89pMOHDys+Pl7Tpk1TaWmpMjIyPGnOnTunlStX6tixY0pKStKsWbP03HPPqXfv3p407qnwrrvuOp/8PPPMM54++AAAAOg8MlN6Ks4in4DearEoI6VH9DIFAFEU1XnmzYx55gEAAC6t7fuqdG/xQTUahqwWi9bOG6W8CenRzhYAtKtOP888AAAA0BJ5E9I19Yr+qnTWKiOlB83rAXRpBPMAAAAwjTS7jSAeANQJ5pkHAAAAAAAtQzAPAAAAAIDJEMwDAAAAAGAyBPMAAAAAAJgMwTwAoFkOV51KjzrlcNVFOysAuiCuQQAQiNHsAQBhbd9XpYLiA2oypDiLVDhvNPM6A7hkuAYBQHDUzAMAQnK46jwP0ZLUZEj3Fh+kdgzAJcE1CABCI5gHAIR03HnW8xDt1mgYqnTWRidDALoUrkEAEBrBPAAgpMyUnoqz+C6zWizKSOkRnQwB6FK4BgFAaATzAICQ0uw2Fc4bLavl4tO01WLR2nmjlGa3RTlnALoCrkEAEJrFMAyj+WTwV1NTI7vdLpfLpeTk5GhnBwA6lMNVp0pnrTJSevAQDeCS4xoEoCuJNNZkNHsAQLPS7DYeoAFEDdcgAAhEM3sAAAAAAEyGYB4AAAAAAJMhmAcAAAAAwGQI5gEAAAAAMBmCeQAAAAAATIZgHgAAAJ2Ow1Wn0qNOOVx10c4KAHRKTE0HAGgRh6tOx51nlZnSk6miAHSI7fuqVFB8QE2GFGeRCueNVt6E9GhnCwA6FYJ5AEDEeMAG0JEcrjrtr/yL5zojSU2GdG/xQU29oj8FiADghWAeABARh6uOB2wAHca7sNBfo2Go0lnLtQYAvNBnHgAQkePOswEP2Y2GoVc+cNCnFUCb+BcW+rNaLMpI6XFpMwUAnRzBPAAgIpkpPRVnCVz+k1c+Us66N7R9X9WlzxSAmBCssNDNarFo7bxR1MoDgB+CeQBARNLsNhXOGy2rJTCidze5p4YeQGsEKyyMk/T4reO0e8U0xuYAgCAI5gEAEcubkK7dK6Zp5Q0jAj5rNAyVV56JQq4AmJ1/YaHVYlHhd0brhisHUiMPACFYDMMI0agJ4dTU1Mhut8vlcik5OTna2QGAS8rhqlPOujcCmsUywj2AtnC46lTprFVGSg+CeABdVqSxJjXzAIAWc9ei+d9EaG4PoC3S7DZNHtqPQB4AIkAwDwBolbwJ6fr5reMClrunkAKASDhcdSo96qQQEABaiHnmAQCtlj2kj+Is8mluzxRSACLlPbc83XQAoGWomQcAtFqwQauWXz9cx51nqWUDEJb/3PJ00wGAlqFmHgDQJnkT0jX1iv6qdNbqg+q/6qevfkwtG4BmBZtb3t1Nhz7zANA8auYBAG2WZrcpI6WHJ5CXqGUDEF6wueXppgMAkSOYBwC0i3C1bBKDXAHwFaybztp5o6iVB4AI0cweANAu3LVswQbDY5ArANLFQr3jzrPKTOkpSRrct4eK75ys2vNNzC0PAC1EMA8AaBfuWrZ7iw+q0TA8g+Htr/xL0EGupl7Rnwd3oAvxLtRzt6439GUB3+Sh/aKZPQAwHYJ5AEC7CTUYnj8GuQK6Fv+R670vCxTwAUDr0GceANCugg2G549BroCuJdiYGt68x9cAAESGYB4A0O7CPbgzyBXQ9QQbud4bBXwA0HI0swcAtLtgg+HFSXr01nEaP6QPgTzQxfiPqWGRJItkGBTwAUBrWQzDCNPoCaHU1NTIbrfL5XIpOTk52tkBgE5n+74qn8Hw1s4bxQj2QBfkP4J9pbPWUwvv/jeBPAB8KdJYk5p5AECH8B4Mj4d1IPZ5B+3u8725aSm5LgBA60W9z3xRUZEyMzOVmJio7OxsvfPOO2HTP/744xoxYoRsNpuGDx+uZ5991ufzhoYGrVmzRkOHDlViYqLGjBmjXbt2tfl7AQAtl2a3afLQfgEP7A5XnUqPOuVw1UUpZwDa0/Z9VcpZ94Zu3bRXOeve0JNvHdVv3q8OOi0l5z0AtI+o1sxv375d+fn5KioqUk5Ojp588knNnDlThw4dUnp6YFPMDRs2qKCgQJs2bdKECRNUVlamRYsWqU+fPrrxxhslSStXrtTWrVu1adMmZWVl6bXXXtPcuXNVWlqqcePGtep7AQDtp7maOgCdn3/Tef+gvfDVj4Oux7SUANB+otpnfuLEiRo/frw2bNjgWTZixAjl5uaqsLAwIP2UKVOUk5OjBx980LMsPz9f+/fv1+7duyVJAwcO1H333aclS5Z40uTm5iopKUlbt25t1fcGQ595AGg5h6tOOeve8BkYz2qxaPeKaZIU0EQXQOfjXyB384TBer7sRETrus93znEACK3T95k/f/68ysvLtWLFCp/l06dPV2lpadB16uvrlZiY6LPMZrOprKxMDQ0Nio+PD5nGHey35nvd311fX+95X1NT0/xOAgB8BJuyrtEw9MzuSj21+xi19UAn53DVBdTCtySQZ9R6AGg/Uesz73Q61djYqNTUVJ/lqampOnnyZNB1ZsyYoaeeekrl5eUyDEP79+/X5s2b1dDQIKfT6Umzfv16ffLJJ2pqalJJSYlefvllORyOVn+vJBUWFsput3tegwcPbsvuA0CXFGyu6TjJE8hL9KsFOrNgBXLNiZP0+K3jtHvFNArpAKAdRX0APIvF96nOMIyAZW6rVq3SzJkzNWnSJMXHx2vOnDlasGCBJMlqtUqSHnnkEQ0bNkxZWVlKSEjQ0qVLtXDhQs/nrfleSSooKJDL5fK8TpyIrBQaAPAl91zT1r9fb60Wi/7lmsygtfWVztoo5BBAOMEK5Pzd+vV0n3O88DujdcOVA6mRB4B2FrVm9ikpKbJarQG14adOnQqoNXez2WzavHmznnzySf3pT39SWlqaNm7cqF69eiklJUWS1L9/f7300ks6d+6cTp8+rYEDB2rFihXKzMxs9fdKUvfu3dW9e/e27DIAQL5T1vVIiFPVX2oVZ1FAP3r3PNQAOgf3oHf3zMzSz149rMYgwy5ZLRZ9/5tf1fe/+VWmpQSADha1YD4hIUHZ2dkqKSnR3LlzPctLSko0Z86csOvGx8dr0KBBkqRt27Zp9uzZiovzbWSQmJioyy+/XA0NDdqxY4duuummNn8vAKB9pNltevvInz19by2SLBbJMOhXC3RG/oPe3XN9lq4c1FsfVP/VE9j7n7ucwwDQsaI6Nd2yZcs0f/58XXXVVZo8ebI2btyoqqoqLV68WNLFpu3V1dWeueSPHDmisrIyTZw4UWfOnNH69et18OBBbdmyxbPNvXv3qrq6WmPHjlV1dbVWr16tpqYmLV++POLvBQB0LP9BtAxJcYb02K3jNH5IH0lS6VFnyJHtvafFImAAOlawQe9+tuuwdq+YpslD++nbYwZSCw8AURDVYD4vL0+nT5/WmjVr5HA4NGrUKO3cuVNDhgyRJDkcDlVVVXnSNzY26qGHHtLhw4cVHx+vadOmqbS0VBkZGZ40586d08qVK3Xs2DElJSVp1qxZeu6559S7d++IvxcA0LGCDaLVJKlvz+4+NfbBRrZnnnrg0go1C4V7vnj3CwBwaUV1nnkzY555AGi9UPPNF985WXOLSn2Wx0n6+a3jlD2kj07VnFNuUamMIPPUE0wAHSPU+cp5BwAdI9JYM+qj2QMAup5go9qvnTdKZ883Bq2xX/p8haYUvqE5j/sG8hIj3wMdyXvQO//zlUAeAKIrqs3sAQBdl/eo9u6+tg5XXcDI9m6hmpEx8j3QMUINekffeADoHKiZBwBETZrdpslD+/mMfu1dY9+cOIu0/PrhOu48K4erriOzCpiOw1Wn0qPOVp0boQa9I5AHgM6DmnkAQKfirrEvrzyjH2yrCFpLL10sjb7zuqH66a6PGQwP8NOWgSIdrjr9zwefhx30DgAQfdTMAwA6nTS7TbPHDPSppXfPRS9dbFp/z8wsFb151Kfm8N7ig9TQo8sLVqse6tzwr73fvq9KOeve0AOvfByQli4tANC5UDMPAOi0/PvVS/L8u7npsoCuKtJzI1ifeHdLF38MegcAnQ/BPACgU/Ofw9r73/6D5VFzCEiZKT2bPTeC1d7/9NWP1RRke6tuGKFZV6YRyANAJ0MzewCAKQWb3i7UYHhtGQgMMJtQUz+6g/FQfeKbdLE7izerxUIgDwCdlMUw/GfsRSRqampkt9vlcrmUnJwc7ewAQJflcNWp0lmrD6r/erFm0W/Ar7YMBAaYmfvc8B6B3vt8aI67EIDzBQAurUhjTZrZAwBMzR2k3PbUHwIG/Moa0CvoQGBTr+hPTSNinn8XFf+m9eHESSq+c7LGDO7TcRkEALQJzewBAKYXasCv/97/x5ADgXmjGT7MqKXHbbDzRJLmTwyseW+SVHs+WA96AEBnQc08AMD0gg34JUlb91YFpHUPBOZw1em486wOVLuCNs8HOrOWdh9xuOp0+ov6oAPjffeqQfplWRWDSQKAyVAzDwAwPf8Bv0KJk/TPV2fo1+9/rpx1b+jWTXtVuPPjFs1VTy0+oq25eeRDzR3//Rfek2FI7tPE3Sd+zOA+YQfMAwB0TtTMAwBigntO+lc+cOgnr3wU8Pns0WnaedChje8cD7udRsPQKx84dEOQEbwZTA/R4G5FkpnSU2l2W8huJa984FCjYfi0NPGfO96QFGdIj906TuOH9PEc4+7zx3/APABA58Vo9q3EaPYA0Dk5XHXKWfeGT7ATJ0lBmuGH4x+sB9uu1WLR7hXTCHzQYYIVIE29on/AsRhKnBR07vgXFk3S5KH92ju7AIB2EGmsSTN7AEBMCTbH9r9ck9miQF4KbLocqjbUfzA9oKVCdd0I1ZxeUkTdSqS/zx3vl4z+8AAQG2hmDwCIOf5NhiXpqd3Hgwb0VotFy68fLmucJaB5vneT+2CD7MVJOn22Xg5XHbXzaJVwXTfCNae/4co07V4xLWS3EjerxaLlM4frZ68eVqNh0B8eAGIIzexbiWb2AGAu2/dV6d7ig56AZvn1w3XloN6e/sHBmtG7uYMsSZ5tuCs7DYXuP+/f1zmYSNIgNoXqEvLzW8cpe0gfnao5p9yiUgV7Uoukyb07cM+bkC6Hq47+8ABgEpHGmgTzrUQwDwDm01xA4x3w+4uT9OKSKbosOVHllWf0g20VYfvPRzJYHgPqmUdbCl1CrVt61KlbN+0Nuo53YVEo7mPu7SN/DltQBQAwl0hjTZrZAwC6jDS7LWxwE25E/CZJuUWlWjdvtAb37RGy/3ya3ab3T5zRiuIDnhpVd1/nqVf093x/qP7Q3mnc6ai5j662FLqEWjfUvO9u/ossQZa5jzlGogeArolgHgAAL2l2m264Mk1rd34UEGQZfw+4i++cHBCEuQcV276vSit2HAgaeHlPeRduQL1Qtfv3XJ+l0YPsAYF9cwE/BQIt4/17SYqo0CXYNvZX/iVg3YIdB3TiL7UqevOomoyLQbrFoqBN6b0ZQdJ5D2TXXEEVACD2EMwDAODHPSJ+wY4DAdN6NRqGas83qXDeaJ+mzWvnjZJ0MfALFZf95JWPtHbnR7rn+iwN7J0YskBACl5zX/jqx5J8A/sD1a6AecW9A36a8ocWrJDD//f6l6sDZ0LwLnRpbhv+miQ99rujnvfued//PXekfvTyhyFnXWAgOwCAP/rMtxJ95gEg9r1/4kzAAGTefeP9++CH6wMdjOXv/zEM38HKpPD9qSPhDux/uuvjsH37u6pgrR4G9k7UXdveCxiQTkEKXdx91dsyB7y3FxZNUtVfzvoOsBjk2GAgOwCIffSZBwCgjcYM7qN1QWrg3UGUd9PmUH2gg/V1dnPXyj526ziNH9LHJzgLNhVeSzQZulhj77fcvyl/VxSu1YO/Jkn/9+qv6Ondx30GmAvVhH7pN7/a4r+Zu0XG5KH9AqZU9A/caU4PAHAjmAcAIIxIBhfzruX17gPtDvz8a8e9NUnq27N70O3+89WZenr38dYH9Arfz7qrCjZeQShxkq4cbFfxnZNVe75JH1T/1dOtwV+TpJ+//mnA8nAj04crIHK/BwAgGIJ5AACaEa421L+WN1hte+8e8SGnvPPvJ3/cedanH7xF0v+dmql+Sd09/aUjFayf9fLrh+u486xnv7qiSFs9uFtVLH2+ImS3heZ4j6fA9HEAgPZEn/lWos88AEAK3bf9hUWTNHloP897d1/nD6r/GjCIWd6E9LCDprn7aEsKug13YOi9PE7Sv1yTqYVXZ/qs5z1YXlcaDM9/oLrt+6pCFrDESVoTZEC6OCmg20I4q24YoVl/n73AnQf6uwMAmhNprEkw30oE8wAA6WKA5j/gWXODzPkHdcG24S9U4YB/YOhw1emZ3ZV6avexZgdmi/XB8IK1dPCf6z1UAcvgvj2CFtL4d1sIFfjH+m8LAOg4DIAHAMAl4J7GLtQgeaHW8f68uT7cwfq5h2v67w7kpS/nRX/klrHNzmsfS0K1dPCfJz7NbtPkof307TEDAwpYgk0dGGx6uLwJ6Yq3xrXoGAAAoK0I5gEAaKNIBskLJ1wf7pYGhsEKBhoNQ/p7rXSoee1jif84Bv6CFWIEG3guWCFN3oT0gMBfavsxAABASxHMAwDQDtoyZViwwLG1A6QFKxiwWizKzujT4hYEZtWalg7BhArQQ/2tmTYOAHApEcwDANAJtFfNbrhm/12l9jhcS4c4Sf98dUbE2yJABwB0VgyA10oMgAcA6My6+sjp3qPVu1s6nP7ifMDAgF1lNH8AgHkwmn0HI5gHAKBz8y7QkNTlRvMHAJgTo9kDAIAuzbuJfOlRZ5cazR8AEPviop0BAACAjubuR+8tVkfzBwB0DQTzAAAg5rkHBrRaLkb0sTyaPwCga6CZPQAA6BK6ymj+AICugWAeAAB0GUw1BwCIFTSzBwAAAADAZAjmAQAAAAAwGYJ5AAAAAABMhmAeAAAAAACTIZgHAAAAAMBkoh7MFxUVKTMzU4mJicrOztY777wTNv3jjz+uESNGyGazafjw4Xr22WcD0jz88MMaPny4bDabBg8erLvvvlvnzp3zfH7hwgWtXLlSmZmZstls+spXvqI1a9aoqamp3fcPAAAAAID2FtWp6bZv3678/HwVFRUpJydHTz75pGbOnKlDhw4pPT09IP2GDRtUUFCgTZs2acKECSorK9OiRYvUp08f3XjjjZKkX/7yl1qxYoU2b96sKVOm6MiRI1qwYIEk6T//8z8lST/96U/1xBNPaMuWLRo5cqT279+vhQsXym6366677rpk+w8AAAAAQGtYDMMwovXlEydO1Pjx47VhwwbPshEjRig3N1eFhYUB6adMmaKcnBw9+OCDnmX5+fnav3+/du/eLUlaunSpPvroI73++uueNP/6r/+qsrIyT63/7NmzlZqaqqefftqT5jvf+Y569Oih5557Lmhe6+vrVV9f73lfU1OjwYMHy+VyKTk5uZW/AAAAAAAAX6qpqZHdbm821oxaM/vz58+rvLxc06dP91k+ffp0lZaWBl2nvr5eiYmJPstsNpvKysrU0NAgSbr66qtVXl6usrIySdKxY8e0c+dO3XDDDZ51rr76ar3++us6cuSIJOn999/X7t27NWvWrJD5LSwslN1u97wGDx7c8p0GAAAAAKAdRK2ZvdPpVGNjo1JTU32Wp6am6uTJk0HXmTFjhp566inl5uZq/PjxKi8v1+bNm9XQ0CCn06m0tDTdfPPN+vOf/6yrr75ahmHowoUL+n//7/9pxYoVnu3cc889crlcysrKktVqVWNjox544AHdcsstIfNbUFCgZcuWed67a+YBAAAAALjUotpnXpIsFovPe8MwApa5rVq1SidPntSkSZNkGIZSU1O1YMEC/exnP5PVapUkvfnmm3rggQdUVFSkiRMn6tNPP9Vdd92ltLQ0rVq1StLFvvpbt27V888/r5EjR+q9995Tfn6+Bg4cqNtvvz3od3fv3l3du3f3yad0MagHAAAAAKA9uGPMZnvEG1FSX19vWK1Wo7i42Gf5D37wA2Pq1Klh1z1//rxx4sQJ48KFC0ZRUZHRq1cvo7Gx0TAMw7j66quNf/u3f/NJ/9xzzxk2m82TZtCgQcZjjz3mk+bf//3fjeHDh0ec/xMnThiSePHixYsXL168ePHixYsXr3Z/nThxImxMGrWa+YSEBGVnZ6ukpERz5871LC8pKdGcOXPCrhsfH69BgwZJkrZt26bZs2crLu5i9//a2lrPv92sVqsMw/CUbIRK05Kp6QYOHKgTJ06oV69eIVsSxAp3l4ITJ04w2B8QBOcIEB7nCNA8zhMgvK50jhiGob/97W8aOHBg2HRRbWa/bNkyzZ8/X1dddZUmT56sjRs3qqqqSosXL5Z0sZ96dXW1Zy75I0eOqKysTBMnTtSZM2e0fv16HTx4UFu2bPFs88Ybb9T69es1btw4TzP7VatW6dvf/ranKf6NN96oBx54QOnp6Ro5cqQqKiq0fv16fe9734s473FxcZ4Cha4iOTk55k8coC04R4DwOEeA5nGeAOF1lXPEbrc3myaqwXxeXp5Onz6tNWvWyOFwaNSoUdq5c6eGDBkiSXI4HKqqqvKkb2xs1EMPPaTDhw8rPj5e06ZNU2lpqTIyMjxpVq5cKYvFopUrV6q6ulr9+/f3BO9ujz76qFatWqU777xTp06d0sCBA3XHHXfoRz/60SXbdwAAAAAAWiuq88zDHCKd5xDoqjhHgPA4R4DmcZ4A4XGOBIraPPMwj+7du+vHP/6xz2j+AL7EOQKExzkCNI/zBAiPcyQQNfMAAAAAAJgMNfMAAAAAAJgMwTwAAAAAACZDMA8AAAAAgMkQzAMAAAAAYDIE89CZM2c0f/582e122e12zZ8/X3/961/DrlNcXKwZM2YoJSVFFotF7733XkCa+vp6ff/731dKSop69uypb3/72/rjH//YMTsBdKDWnCOGYWj16tUaOHCgbDabrrvuOn344Yc+aa677jpZLBaf180339yBewK0n6KiImVmZioxMVHZ2dl65513wqZ/6623lJ2drcTERH3lK1/RE088EZBmx44d+trXvqbu3bvra1/7ml588cWOyj7Q4dr7HPnFL34RcM+wWCw6d+5cR+4G0GFaco44HA7deuutGj58uOLi4pSfnx80XVe7jxDMQ7feeqvee+897dq1S7t27dJ7772n+fPnh13n7NmzysnJ0bp160Kmyc/P14svvqht27Zp9+7d+uKLLzR79mw1Nja29y4AHao158jPfvYzrV+/Xo899pj27dunAQMG6Fvf+pb+9re/+aRbtGiRHA6H5/Xkk0925K4A7WL79u3Kz8/Xfffdp4qKCl1zzTWaOXOmqqqqgqY/fvy4Zs2apWuuuUYVFRW699579YMf/EA7duzwpNmzZ4/y8vI0f/58vf/++5o/f75uuukm7d2791LtFtBuOuIckaTk5GSfe4bD4VBiYuKl2CWgXbX0HKmvr1f//v113333acyYMUHTdMn7iIEu7dChQ4Yk4w9/+INn2Z49ewxJxscff9zs+sePHzckGRUVFT7L//rXvxrx8fHGtm3bPMuqq6uNuLg4Y9euXe2Wf6CjteYcaWpqMgYMGGCsW7fOs+zcuXOG3W43nnjiCc+ya6+91rjrrrs6LO9AR/n6179uLF682GdZVlaWsWLFiqDply9fbmRlZfksu+OOO4xJkyZ53t90003G9ddf75NmxowZxs0339xOuQYunY44R5555hnDbre3e16BaGjpOeIt1PNTV7yPUDPfxe3Zs0d2u10TJ070LJs0aZLsdrtKS0tbvd3y8nI1NDRo+vTpnmUDBw7UqFGj2rRd4FJrzTly/PhxnTx50uf47969u6699tqAdX75y18qJSVFI0eO1L/9278F1NwDnc358+dVXl7uc3xL0vTp00OeE3v27AlIP2PGDO3fv18NDQ1h03DPgNl01DkiSV988YWGDBmiQYMGafbs2aqoqGj/HQA6WGvOkUh0xftIt2hnANF18uRJXXbZZQHLL7vsMp08ebJN201ISFCfPn18lqemprZpu8Cl1ppzxL08NTXVZ3lqaqo+++wzz/vbbrtNmZmZGjBggA4ePKiCggK9//77Kikpacc9ANqX0+lUY2Nj0OM73DkRLP2FCxfkdDqVlpYWMg33DJhNR50jWVlZ+sUvfqHRo0erpqZGjzzyiHJycvT+++9r2LBhHbY/QHtrzTkSia54H6FmPkatXr066CAp3q/9+/dLkiwWS8D6hmEEXd5WHbVdoKUuxTni/7n/OosWLdI//MM/aNSoUbr55pv1q1/9Sv/7v/+rd999tx32EOhYzR3fkaT3X97SbQKdWXufI5MmTdI//uM/asyYMbrmmmv0X//1X7riiiv06KOPtnPOgUujI675Xe0+Qs18jFq6dGmzo2JnZGTogw8+0J/+9KeAz/785z8HlGy1xIABA3T+/HmdOXPGp3b+1KlTmjJlSqu3C7SXjjxHBgwYIOliCXFaWppn+alTp8KeV+PHj1d8fLw++eQTjR8/PpLdAC65lJQUWa3WgJqOcMf3gAEDgqbv1q2b+vXrFzZNW+5FQDR01DniLy4uThMmTNAnn3zSPhkHLpHWnCOR6Ir3EWrmY1RKSoqysrLCvhITEzV58mS5XC6VlZV51t27d69cLlebgu7s7GzFx8f7NBd2OBw6ePAgwTw6hY48R9xN572P//Pnz+utt94Ke/x/+OGHamho8CkAADqbhIQEZWdnB3QHKSkpCXl8T548OSD9b3/7W1111VWKj48Pm4Z7Bsymo84Rf4Zh6L333uOeAdNpzTkSiS55H4nOuHvoTK6//nrjyiuvNPbs2WPs2bPHGD16tDF79myfNMOHDzeKi4s970+fPm1UVFQYr7zyiiHJ2LZtm1FRUWE4HA5PmsWLFxuDBg0y/vd//9d49913jW984xvGmDFjjAsXLlyyfQPaQ2vOkXXr1hl2u90oLi42Dhw4YNxyyy1GWlqaUVNTYxiGYXz66afG/fffb+zbt884fvy48corrxhZWVnGuHHjOEfQ6W3bts2Ij483nn76aePQoUNGfn6+0bNnT6OystIwDMNYsWKFMX/+fE/6Y8eOGT169DDuvvtu49ChQ8bTTz9txMfHG7/61a88aX7/+98bVqvVWLdunfHRRx8Z69atM7p16+YzkwRgFh1xjqxevdrYtWuXcfToUaOiosJYuHCh0a1bN2Pv3r2XfP+AtmrpOWIYhlFRUWFUVFQY2dnZxq233mpUVFQYH374oefzrngfIZiHcfr0aeO2224zevXqZfTq1cu47bbbjDNnzvikkWQ888wznvfPPPOMISng9eMf/9iTpq6uzli6dKnRt29fw2azGbNnzzaqqqouzU4B7ag150hTU5Px4x//2BgwYIDRvXt3Y+rUqcaBAwc8n1dVVRlTp041+vbtayQkJBhDhw41fvCDHxinT5++RHsFtM3jjz9uDBkyxEhISDDGjx9vvPXWW57Pbr/9duPaa6/1Sf/mm28a48aNMxISEoyMjAxjw4YNAdv87//+b2P48OFGfHy8kZWVZezYsaOjdwPoMO19juTn5xvp6elGQkKC0b9/f2P69OlGaWnppdgVoEO09BwJFnsMGTLEJ01Xu49YDOPvo2sAAAAAAABToM88AAAAAAAmQzAPAAAAAIDJEMwDAAAAAGAyBPMAAAAAAJgMwTwAAAAAACZDMA8AAAAAgMkQzAMAAAAAYDIE8wAAAAAAmAzBPAAACLB69WqNHTv2kn/vm2++KYvFIovFotzc3LBpr7vuOuXn5/u8d6/73nvvdWg+AQCItm7RzgAAALi0LBZL2M9vv/12PfbYY/r+979/iXIU6PDhw7rssstatE5xcbGOHj2qr3/96x2UKwAAOg+CeQAAuhiHw+H59/bt2/WjH/1Ihw8f9iyz2WxKSkpSUlJSNLInSbrsssvUu3fvFq3Tt29f1dTUdEyGAADoZGhmDwBAFzNgwADPy263y2KxBCzzb2a/YMEC5ebmau3atUpNTVXv3r11//3368KFC/rhD3+ovn37atCgQdq8ebPPd1VXVysvL099+vRRv379NGfOHFVWVrY4z2fPntU//dM/KSkpSWlpaXrooYfa+CsAAGBuBPMAACAib7zxhj7//HO9/fbbWr9+vVavXq3Zs2erT58+2rt3rxYvXqzFixfrxIkTkqTa2lpNmzZNSUlJevvtt7V7924lJSXp+uuv1/nz51v03T/84Q/1u9/9Ti+++KJ++9vf6s0331R5eXlH7CYAAKZAMA8AACLSt29f/fznP9fw4cP1ve99T8OHD1dtba3uvfdeDRs2TAUFBUpISNDvf/97SdK2bdsUFxenp556SqNHj9aIESP0zDPPqKqqSm+++WbE3/vFF1/o6aef1n/8x3/oW9/6lkaPHq0tW7aosbGxg/YUAIDOjz7zAAAgIiNHjlRc3Jf1AKmpqRo1apTnvdVqVb9+/XTq1ClJUnl5uT799FP16tXLZzvnzp3T0aNHI/7eo0eP6vz585o8ebJnWd++fTV8+PDW7goAAKZHMA8AACISHx/v895isQRd1tTUJElqampSdna2fvnLXwZsq3///hF/r2EYrcgtAACxjWAeAAB0iPHjx2v79u267LLLlJyc3OrtfPWrX1V8fLz+8Ic/KD09XZJ05swZHTlyRNdee217ZRcAAFOhzzwAAOgQt912m1JSUjRnzhy98847On78uN566y3ddddd+uMf/xjxdpKSkvTP//zP+uEPf6jXX39dBw8e1IIFC3ya/AMA0NVQMw8AADpEjx499Pbbb+uee+7RvHnz9Le//U2XX365vvnNb7a4pv7BBx/UF198oW9/+9vq1auX/vVf/1Uul6uDcg4AQOdnMeiIBgAAOok333xT06ZN05kzZ9S7d+8Wr19ZWanMzExVVFRo7Nix7Z4/AAA6C9qnAQCATmfQoEG65ZZbWrTOzJkzNXLkyA7KEQAAnQs18wAAoNOoq6tTdXW1pIt95QcMGBDxutXV1aqrq5MkpaenKyEhoUPyCABAZ0AwDwAAAACAydDMHgAAAAAAkyGYBwAAAADAZAjmAQAAAAAwGYJ5AAAAAABMhmAeAAAAAACTIZgHAAAAAMBkCOYBAAAAADAZgnkAAAAAAEzm/wOoE9pNEXhNQgAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "time, flux = lcs(radius_ratio=0.1, zero_epoch=0.0, period=2.0,\n", " scaled_semi_major_axis=8.0, impact_parameter=0.1,\n", " limb_darkening=[0.2, 0.3])\n", "\n", "fig, ax = plt.subplots(1, 1, figsize=(10,4), constrained_layout=True)\n", "ax.plot(time, flux, '.')\n", "plt.setp(ax, xlabel='Time [d]', ylabel='Flux');" ] }, { "cell_type": "markdown", "id": "ee07ff8a-dca8-49f2-9aba-08802cdfd188", "metadata": {}, "source": [ "The resulting light curve contains gaps with a period of CHEOPS's rotation period of 99 minutes. The fractional width and phase of the gaps can be changed with the `efficiency` and `eff_phase` arguments. Both `efficiency` and `eff_phase` range from 0 to 1; `efficiency` has a default value of 0.8, and `eff_phase` defaults to 0.0." ] }, { "cell_type": "code", "execution_count": 5, "id": "01b4f7df-6cbf-482f-a431-0e82991e0d92", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAA/MAAAGbCAYAAACIxMC9AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8hTgPZAAAACXBIWXMAAA9hAAAPYQGoP6dpAABPkElEQVR4nO39fXxU9YH3/79PhoQMhAw3gRCEkJRLCQXkJlBuorTsriCIJdLf1XizXEK7/MpVqEZ2i0ShRVYJamW1ahBQLEotubZEbQtiU61IGkogRgVRUCCGxqF0WDqpSQghOd8/6EwzycxkEpJMTvJ6Ph7jg5z5nJnPiTlz5n0+d4ZpmqYAAAAAAIBlRIS7AgAAAAAAoGUI8wAAAAAAWAxhHgAAAAAAiyHMAwAAAABgMYR5AAAAAAAshjAPAAAAAIDFEOYBAAAAALCYHuGugFXV19friy++UJ8+fWQYRrirAwAAAADoAkzT1N/+9jcNGTJEERGB298J8630xRdfaNiwYeGuBgAAAACgCzpz5oyGDh0a8HnCfCv16dNH0pVfcGxsbJhrAwAAAADoCioqKjRs2DBv5gyEMN9Knq71sbGxhHkAAAAAQJtqbjg3E+ABAAAAAGAxhHkAAAAAACyGMA8AAAAAgMUQ5gEAAAAAsBjCPAAAAAAAFkOYBwAAAADAYgjzAAAAAABYDGEeAAAAAACLIcwDAAAAAGAxhHkAAAAAACyGMI9OwemuVuFJl5zu6nBXBQAAAAA6vR7hrgC6L6e7WqddlTpS7tajb3yielOKMKTsBWOVMTkx3NWzNM/vNjmutxIc9nBXBwAAhEFHfR/gewcQHoR5eHXkB3HuoTJl5R1Rvem7vd6UHsg7qhnXDbyqOrT2WDr6YtSS92tYVlLA/Rr+brk5AgBA1xLqd4eO+j7A9w4gfAjz3UiwD/9gH8RtHXCd7mq/Qd6jzjRV6qpq9Xs1Ppb7b07R2KGOFl/0Qt3Pc0wt/R215OLXsKzx922m/P+/avi7baubIwAAoP209Q37tvw+EOw7Dt870Bl0554hYQ3z7777rh5//HEVFxfL6XTq1VdfVXp6etB99u3bpxUrVuijjz7SkCFDtHLlSi1dutSnzK5du7RmzRqdPHlSI0aM0COPPKLbbrvNp0xOTo4ef/xxOZ1OjR49Wk8++aRuvPHGtj7ETqO5sB7og/jdE39p9d3WQBem067KgEFekmyGoaS4Xq06Tn/Hkv3GJ5JaftHzt5+/D4vW3ARoycWvcdmGv7p6U8radUS9e/bQsH52FZX+T5Pf7dXeHAEAAO2nPW7Y+/uu1ZLvA6EMhXS6q/WbD7/gewfCqrv3DAlrmK+srNS4ceO0ePFifetb32q2/OnTpzV37lwtWbJEO3bs0B/+8Ad9//vf18CBA737HzhwQBkZGfrP//xP3XbbbXr11Vf17W9/WwUFBZoyZYokKTc3V5mZmcrJyVFaWpo2b96sOXPm6NixY0pM7Hr/85v78A/0gV9ceqHFd1v9ffg3vjDdf3OKIgz5DfQ2w9D6BWNafQEIdqOgpRe9xvv9tbrW54J2/80pGtI3ulU3D0K9yAa6UPrUT9LyV0oCPn81N0cAAMA/+LupH6hVMJSyod6wTx3er9nvDg1fOzmud5PvWhGSzlfWeCcbDqX1v6FAjT2N8b0DHSVYxpH8/413tVb8sIb5OXPmaM6cOSGXf+6555SYmKgnn3xSkjRq1CgdPnxYP/nJT7xh/sknn9RNN92krKwsSVJWVpb27dunJ598Ur/4xS8kSRs3btR3v/td/du//Zt3nzfffFObNm1SdnZ2Gx5h+AW7a7r7Q6cmJ/XT+S9rmnzg2wxD8hO4g91tDfTh3/jC9Nje47p/Tooee+O46kxTNsPQyptH6vqhfZUU1yvgCScFvvB4+Lt4hVL/UPbb8MYnMv2E9kACfaBICvg7T4rr5feGSGtd7c0RAAC6u2Ct1JL89s4LpWz2grEa1r9XSDfsAzWGeAL65ndPNnm/7AVj9UDeUdWZpgxd+T62/JUSv40snh6FDevoj7/GnoYiJH33hqSAv8OuEqDQOQS6wfViQameLzjVpLW+K7biW2rM/IEDBzRr1iyfbbNnz9YLL7yg2tpaRUZG6sCBA7rvvvualPHcALh06ZKKi4u1atUqnzKzZs1SYWFhwPeuqalRTU2N9+eKioqrPJr2Fyhcezy8+2Pvvw1JhiGZ5j8CYOrwfgEDZ0NOd7UOl/5P0PdqqM40df01fVWwaqZKXVVNAry/+gfrdtZQgsPuc/FqLNDd4ub2i1DgoB9MnWnq6bc+085DZU2Oo/HvfOXNI7Wt4LReKDgd9L2Mv//HTzV9rLlllOZen8BFEwCAVgrWSr1q1xGf63GgG/3esvpHA4fnhn/e96cFbUxo+BqPvvGJvj9zhDa9c6pJQG9c9oG8oypYNVMFq2aquPSC7tlZErD1v2GPwtsnD2t2KKS/xh5Jmjc2QXuOOrVl/2k9X3A6YIBqyZxEQDD+GuMMSVsLTvmcl1m7juhSXb1+/PpHXW5+B0uF+bNnzyo+Pt5nW3x8vC5fviyXy6WEhISAZc6ePStJcrlcqqurC1rGn+zsbD300ENtdCTtr7lJ5hozJUWY0jN3TtDE4f28f9QNA66/Vt7mbhj447mLnBTXS9NGDPBb98Y3BxpfeAK1eJ92VWrGdQO9Nwo+LP+rTw+AYK3UGZMTNeO6gX73W3nzSD26t3Wt5K8UlXn/3XD3hr/zP/21Whv2fKLmXt4T0CU1uTg3ZDMMgjwAAFehue9Spvc/zfNXrM40daj0gk9vxWA37OslPfvOSa26OUXX9LUH/A7gee3dHzp1y/UJ6h8TFdL3l3pTeqXoTMDngzX2REjac9TZJCilDO7TquGIgD+Ne3j4a4zzd17WS1rz2kdNXq8rzO9gqTAvSYZh+Pxs/v1/XMPt/so03hZKmYaysrK0YsUK788VFRUaNmxYyyrfgQKNAV84JVEvHyxr+oSu/KH3793T5w+6YcBt3IL+wZkLWpV3pNkW4oYXpoZ3kf19iId6c6BxF5pALffTRgzQN8cNUamrSr2iIlR5qa7JWLGG//Z8MPjbr+HFtrEISU//PZQ/uucT1Qevvle9rvxeHn2j+SDfOKDPG2dX5aXLTXoT0LUeAICr19yEvW3h4d0fe1urPcMNpcA37M2/D1d86o7xzdbt4d0fa/2ej4POVdQcf3U77apsMlzyuzckacv+0z771pmm/vvwn4LOSdRwTgC+tyCYQF3kMyYnKmVwH6XnFDabSRrrCvM7WCrMDx48uEnr+blz59SjRw8NGDAgaBlPS3xcXJxsNlvQMv707NlTPXv2bIvD6BD+up3YDEP/v0lD9fOismYnLPF356vh9iPl7qAtyY3HwUtNL0yNu7e0pDdBhOQN8lLglntP3RtO1NK4q7vn341vLjTer+EFzV+L/y3XD5EkXdPXHnRCusa/p0Dd1RqX8xfQG95s6RUVoapL9QGHLQAAgNA1N59OW/HMJ1SwamaTG/ZZu440aSCoM03p799LQume33iuolCH60nST2+foHnjrny/8dddvuH3vOf9DBPcEaAByVs/BW7gATyam8y78lJdq4J8V2j8slSYnzZtmn7961/7bPvtb3+rSZMmKTIy0lsmPz/fZ9z8b3/7W02fPl2SFBUVpdTUVOXn5/ssV5efn6/58+d3wFF0jMbdTjx/sOOG9fM7NrzhH3SgO1+htJp7Wqgn+rnD2j8m+Cyswe6AN7zwREiaO3awfnMk8LCIxrO7BuqyH+wmgL8PDs/FtmHLfePw7K/7mb/jCNZdTWp6oQz0YdPwZgsAAGgb/r5LLf36V5Sz76RPcGjYO6/xED1/3eE9vRQb8tfdN1CLo80wlJrUr0ndVt48UrYIw2dOJM9rN56rSFKTYYWNed5H8h+mGt+ACDb3UHO6yvhltI9AE915hpIEWr1hXfpo/ajBOHnP9kBZxYrCGua//PJLffbZZ96fT58+rffff1/9+/dXYmKisrKyVF5erpdeekmStHTpUj3zzDNasWKFlixZogMHDuiFF17wzlIvSffee69mzJihRx99VPPnz9frr7+u3/3udyooKPCWWbFihRYuXKhJkyZp2rRp2rJli8rKypqsV291gbrIB2vNDXTnq/GYJ388wd/TQt1YoN4CnotKoBPRc8JJ8natDxbkG79uS7rJNXdzoeHzgUK0v4t/414Kjf+fNCwfIenfZiRrcVpyl/iQAQDAqvx9l0oc0KtJY4nnu4+/G/0Nh8QFmocnUHffccP6aUOA+Yv81c3prtb6PR/7fe3G31saDysMNs9QKEvqeuqz+0NnkxsK0pV5fy6bZsCbB11h/DLaR6BeMp6hJI1Xb/D8/WZMTlSkLSLg+doVGKbZittnbeSdd97RzJkzm2y/++679bOf/UyLFi1SaWmp3nnnHe9z+/bt03333aePPvpIQ4YM0f33398khP/yl7/U6tWrderUKY0YMUKPPPKIFixY4FMmJydHjz32mJxOp8aMGaP/+q//0owZM0Kue0VFhRwOh9xut2JjY1t24J1Y4UmX7tx6sMn21beM8vvB7BEh6dVl0zVuWL+gr597qMzviRbK8053tdI2vB1yi3co+zVmMwzvXWZ/+zV8vjlOd3XQ2fqvtjwAAAiPq73GN/d96Grer6WvHcr7tOQ7UXNlne5qv3MCtOQ7Frqfhn/XjXn+dqSmDWaSNb9jh5o1wxrmrayrhvlAH8B535+m23IKA461b4sLRXPPB7rR0HCG90Cv2/ADoPGEfIFuAjTer6XHCQAArKkj1kVvz4DRHq/dku9EoZTlOxZayumuDtjz4xdLpvpdJcuqCPPtrKuGeSnwh2vj7Q27jnfEXa62bCmX5PffLbm5AAAAup5AcwehZd+JQinLdyy01NXmAasgzLezrhzmpeDdrML5octdXAAA0F66S1AArKw75IFQs6alZrNHxwk2uVs4L2bB1r0HAAC4GqFM9AYgvMgD/0CYh+WE+4YCAADomppbeQdA50AeuCIi3BUAAAAAOgPP8rI2w5CkJku0AUBnQss8AAAA8Hd04QVgFYR5AAAAoAG68AKwArrZAwAAAABgMYR5AAAAAAAshjAPAAAAAIDFEOYBAAAAALAYwjwAAAC6Pae7WoUnXXK6q8NdFQAICbPZAwAAoFvLPVSmrLwjqjelCEPKXjBWGZMTw10tAAiKlnkAAAB0W053tTfIS1K9KT2Qd5QWegCdHmEeAAAA3dZpV6U3yHvUmaZKXVXhqRAAhIgwDwAAgG4rOa63IgzfbTbDUFJcr/BUCABCRJgHAABAt5XgsCt7wVjZjCuJ3mYYWr9gjBIc9jDXDACCYwI8AAAAdGsZkxM147qBKnVVKSmuF0EegCUQ5gEAANDtJTjshHgAlkI3ewAAAAAALIYwDwAAAACAxRDmAQAAAACwGMI8AAAAAAAWQ5gHAAAAAMBiCPMAAAAAAFgMYR4AAAAAAIshzAMAAAAAYDGEeQAAAAAALIYwDwAAAACAxRDmAQAAAACwGMI8AAAAAAAWQ5gHAAAAAMBiCPMAAAAAAFgMYR4AAAAAAIshzAMAAKBbcrqrVXjSJae7OtxVAYAW6xHuCgAAAAAdLfdQmbLyjqjelCIMKXvBWGVMTgx3tQAgZLTMAwAAoFtxuqu9QV6S6k3pgbyjtNADsBTCPAAAALqV065Kb5D3qDNNlbqqwlMhAGgFwjwAAAC6leS43oowfLfZDENJcb3CUyEAaAXCPAAAALqVBIdd2QvGymZcSfQ2w9D6BWOU4LCHuWYAEDomwAMAAEC3kzE5UTOuG6hSV5WS4noR5AFYDmEeAAAA3VKCw06IB2BZdLMHAAAAAMBiwh7mc3JylJycrOjoaKWmpmr//v1Byz/77LMaNWqU7Ha7Ro4cqZdeesnn+draWq1bt04jRoxQdHS0xo0bp7179/qUuXz5slavXq3k5GTZ7XZ95Stf0bp161RfX9/mxwcAAAAAQFsLazf73NxcZWZmKicnR2lpadq8ebPmzJmjY8eOKTExsUn5TZs2KSsrS1u3btXkyZNVVFSkJUuWqF+/frr11lslSatXr9aOHTu0detWpaSk6M0339Rtt92mwsJCTZgwQZL06KOP6rnnntP27ds1evRoHT58WIsXL5bD4dC9997bob8DAAAAAABayjBN02y+WPuYMmWKJk6cqE2bNnm3jRo1Sunp6crOzm5Sfvr06UpLS9Pjjz/u3ZaZmanDhw+roKBAkjRkyBA9+OCDWrZsmbdMenq6YmJitGPHDknSvHnzFB8frxdeeMFb5lvf+pZ69eqll19+OaS6V1RUyOFwyO12KzY2tmUHDgAAAACAH6FmzbB1s7906ZKKi4s1a9Ysn+2zZs1SYWGh331qamoUHR3ts81ut6uoqEi1tbVBy3jCviTdcMMNeuutt3TixAlJ0gcffKCCggLNnTs3YH1rampUUVHh8wAAAAAAIBzCFuZdLpfq6uoUHx/vsz0+Pl5nz571u8/s2bP1/PPPq7i4WKZp6vDhw9q2bZtqa2vlcrm8ZTZu3KhPP/1U9fX1ys/P1+uvvy6n0+l9nfvvv1933HGHUlJSFBkZqQkTJigzM1N33HFHwPpmZ2fL4XB4H8OGDWuD3wIAAAAAAC0X9gnwDMPw+dk0zSbbPNasWaM5c+Zo6tSpioyM1Pz587Vo0SJJks1mkyQ99dRTuvbaa5WSkqKoqCgtX75cixcv9j4vXRmrv2PHDr3yyit67733tH37dv3kJz/R9u3bA9YzKytLbrfb+zhz5sxVHjkAAAAAAK0TtjAfFxcnm83WpBX+3LlzTVrrPex2u7Zt26aqqiqVlpaqrKxMSUlJ6tOnj+Li4iRJAwcO1GuvvabKykp9/vnn+uSTTxQTE6Pk5GTv6/zwhz/UqlWrdPvtt2vs2LFauHCh7rvvPr/j9D169uyp2NhYnwcAAAAAAOEQtjAfFRWl1NRU5efn+2zPz8/X9OnTg+4bGRmpoUOHymazaefOnZo3b54iInwPJTo6Wtdcc40uX76sXbt2af78+d7nqqqqmpS32WwsTQcAAAAAsISwLk23YsUKLVy4UJMmTdK0adO0ZcsWlZWVaenSpZKudG0vLy/3riV/4sQJFRUVacqUKbpw4YI2btyoo0eP+nSPP3jwoMrLyzV+/HiVl5dr7dq1qq+v18qVK71lbr31Vj3yyCNKTEzU6NGjVVJSoo0bN+o73/lOx/4CAAAAAABohbCG+YyMDJ0/f17r1q2T0+nUmDFjtGfPHg0fPlyS5HQ6VVZW5i1fV1enJ554QsePH1dkZKRmzpypwsJCJSUlectcvHhRq1ev1qlTpxQTE6O5c+fq5ZdfVt++fb1lnn76aa1Zs0bf//73de7cOQ0ZMkTf+9739KMf/aijDh0AAAAAgFYL6zrzVsY68wAAANbjdFfrtKtSyXG9leCwh7s6ANBEqFkzrC3zAAAAQEfJPVSmrLwjqjelCEPKXjBWGZMTw10tAGiVsC9NBwAAALQ3p7vaG+Qlqd6UHsg7Kqe7OrwVA4BWIswDAACgyzvtqvQGeY8601Spqyo8FQKAq0SYBwAAQJeXHNdbEYbvNpthKCmuV3gqBABXiTAPAACALi/BYVf2grGyGVcSvc0wtH7BGCbBA2BZTIAHAACAbiFjcqJmXDdQpa4qJcX1IsgDsDTCPAAAALqNBIedEA+gS6CbPQAAAAAAFkOYBwAAAADAYgjzAAAAAABYDGEeAAAAAACLIcwDAAAAAGAxhHkAAAAAACyGMA8AAAAAgMUQ5gEAAAAAsBjCPAAAAAAAFkOYBwAAAADAYgjzAAAAAABYDGEeAAAAAACLIcwDAAAAAGAxhHkAAAAAACyGMA8AAAAAgMUQ5gEAAAAAsBjCPAAAALo8p7tahSddcrqrw10VAGgTPcJdAQAAAKA95R4qU1beEdWbUoQhZS8Yq4zJieGuFgBcFVrmAQAA0GU53dXeIC9J9ab0QN5RWugBWB5hHgAAAF3WaVelN8h71JmmSl1V4akQALQRwjwAAAC6rOS43oowfLfZDENJcb3CUyEAaCOEeQAAAHRZCQ67sheMlc24kuhthqH1C8YowWEPc80A4OowAR4AAAC6tIzJiZpx3UCVuqqUFNeLIA+gSyDMAwAAoMtLcNgJ8QC6FLrZAwAAAABgMYR5AAAAAAAshjAPAAAAAIDFEOYBAAAAALAYwjwAAAAAABZDmAcAAAAAwGII8wAAAAAAWAxhHgAAAAAAiyHMAwAAAABgMYR5AAAAAAAshjAPAAAAAIDFEOYBAAAAALCYsIf5nJwcJScnKzo6Wqmpqdq/f3/Q8s8++6xGjRolu92ukSNH6qWXXvJ5vra2VuvWrdOIESMUHR2tcePGae/evU1ep7y8XP/6r/+qAQMGqFevXho/fryKi4vb9NgAAAAAAGgPPcL55rm5ucrMzFROTo7S0tK0efNmzZkzR8eOHVNiYmKT8ps2bVJWVpa2bt2qyZMnq6ioSEuWLFG/fv106623SpJWr16tHTt2aOvWrUpJSdGbb76p2267TYWFhZowYYIk6cKFC0pLS9PMmTP1xhtvaNCgQTp58qT69u3bkYcPAAAAAECrGKZpmuF68ylTpmjixInatGmTd9uoUaOUnp6u7OzsJuWnT5+utLQ0Pf74495tmZmZOnz4sAoKCiRJQ4YM0YMPPqhly5Z5y6SnpysmJkY7duyQJK1atUp/+MMfmu0FEExFRYUcDofcbrdiY2Nb/ToAAAAAAHiEmjXD1s3+0qVLKi4u1qxZs3y2z5o1S4WFhX73qampUXR0tM82u92uoqIi1dbWBi3jCfuS9Ktf/UqTJk3S//7f/1uDBg3ShAkTtHXr1qD1rampUUVFhc8DAAAAAIBwCFuYd7lcqqurU3x8vM/2+Ph4nT171u8+s2fP1vPPP6/i4mKZpqnDhw9r27Ztqq2tlcvl8pbZuHGjPv30U9XX1ys/P1+vv/66nE6n93VOnTqlTZs26dprr9Wbb76ppUuX6p577mky/r6h7OxsORwO72PYsGFt8FsAAAAAAKDlwj4BnmEYPj+bptlkm8eaNWs0Z84cTZ06VZGRkZo/f74WLVokSbLZbJKkp556Stdee61SUlIUFRWl5cuXa/Hixd7nJam+vl4TJ07U+vXrNWHCBH3ve9/TkiVLfLr7N5aVlSW32+19nDlz5iqPHAAAAACA1glbmI+Li5PNZmvSCn/u3LkmrfUedrtd27ZtU1VVlUpLS1VWVqakpCT16dNHcXFxkqSBAwfqtddeU2VlpT7//HN98skniomJUXJysvd1EhIS9NWvftXntUeNGqWysrKA9e3Zs6diY2N9HgAAAAAAhEPYwnxUVJRSU1OVn5/vsz0/P1/Tp08Pum9kZKSGDh0qm82mnTt3at68eYqI8D2U6OhoXXPNNbp8+bJ27dql+fPne59LS0vT8ePHfcqfOHFCw4cPv8qjAgAAAACg/YV1aboVK1Zo4cKFmjRpkqZNm6YtW7aorKxMS5culXSla3t5ebl3LPuJEydUVFSkKVOm6MKFC9q4caOOHj2q7du3e1/z4MGDKi8v1/jx41VeXq61a9eqvr5eK1eu9Ja57777NH36dK1fv17f/va3VVRUpC1btmjLli0d+wsAAAAAAKAVwhrmMzIydP78ea1bt05Op1NjxozRnj17vC3kTqfTp+t7XV2dnnjiCR0/flyRkZGaOXOmCgsLlZSU5C1z8eJFrV69WqdOnVJMTIzmzp2rl19+2WcN+cmTJ+vVV19VVlaW1q1bp+TkZD355JO66667OurQAQAAAABotbCuM29lrDMPAAAAAGhrnX6deQAAAAAA0DqEeQAAAHRJTne1Ck+65HRXh7sqANDmwjpmHgAAAGgPuYfKlJV3RPWmFGFI2QvGKmNyYrirBQBthpZ5AAAAdClOd7U3yEtSvSk9kHeUFnoAXQphHgAAAF3KaVelN8h71JmmSl1V4akQALQDwjwAAAC6lOS43oowfLfZDENJcb3CUyEAaAeEeQAAAHQpCQ67sheMlc24kuhthqH1C8YowWEPc80AoO0wAR4AAAC6nIzJiZpx3UCVuqqUFNeLIA+gyyHMAwAAoEtKcNgJ8QC6LLrZAwAAAABgMYR5AAAAAAAshjAPAAAAAIDFEOYBAAAAALAYwjwAAAAAABZDmAcAAAAAwGII8wAAAAAAWAxhHgAAAAAAiyHMAwAAAABgMYR5AAAAAAAshjAPAAAAAIDFEOYBAAAAALAYwjwAAAAAABZDmAcAAECX4XRXq/CkS053dbirAgDtqke4KwAAAAC0hdxDZcrKO6J6U4owpOwFY5UxOTHc1QKAdkHLPAAAACzP6a72BnlJqjelB/KO0kIPoMsizAMAAMDyTrsqvUHeo840VeqqCk+FAKCdEeYBAABgeclxvRVh+G6zGYaS4nqFp0IA0M4I8wAAALC8BIdd2QvGymZcSfQ2w9D6BWOU4LCHuWYA0D5aNQFedXW17Hb/H4xOp1MJCQlXVSkAAACgpTImJ2rGdQNV6qpSUlwvgjyALq1VLfMTJkzQe++912T7L3/5S11//fVXXSkAAACgNRIcdk0bMYAgD6DLa1WYv+mmmzR9+nRt2LBBpmnqyy+/1KJFi3T33XfrRz/6UVvXEQAAAAAANNCqbvZPP/20brnlFi1evFi7d+/WF198odjYWB06dEhf/epX27qOAAAAAACggVaFeUmaNWuWFixYoE2bNqlHjx769a9/TZAHAAAAAKADtKqb/cmTJzVt2jT95je/0ZtvvqmVK1dq/vz5WrlypWpra9u6jgAAAAAAoIFWhfnx48crOTlZH3zwgW666SY9/PDDevvtt5WXl6evfe1rbV1HAAAAAADQQKvCfE5Ojnbu3Km+fft6t02fPl0lJSWaOHFiW9UNAAAAAAD4YZimaYa7ElZUUVEhh8Mht9ut2NjYcFcHAAAAANAFhJo1WzUB3ksvvRTwOcMwtHDhwta8LAAAAAAACEGrWub79evn83Ntba2qqqoUFRWlXr166X/+53/arIKdFS3zAAAAAIC2FmrWbNWY+QsXLvg8vvzySx0/flw33HCDfvGLX7S60gAAAAAAoHmtCvP+XHvttdqwYYPuvffetnpJAAAAAADgR5uFeUmy2Wz64osv2vIlAQAAAABAI62aAO9Xv/qVz8+macrpdOqZZ55RWlpam1QMAAAAAAD416own56e7vOzYRgaOHCg/umf/klPPPFEW9QLAAAAAAAE0Kpu9vX19T6Puro6nT17Vq+88ooSEhJa9Fo5OTlKTk5WdHS0UlNTtX///qDln332WY0aNUp2u10jR45sskxebW2t1q1bpxEjRig6Olrjxo3T3r17A75edna2DMNQZmZmi+oNAAAAAEC4tOmY+ZbKzc1VZmamHnzwQZWUlOjGG2/UnDlzVFZW5rf8pk2blJWVpbVr1+qjjz7SQw89pGXLlunXv/61t8zq1au1efNmPf300zp27JiWLl2q2267TSUlJU1e79ChQ9qyZYuuv/76djtGAAAAAADaWsjrzK9YsSLkF924cWNI5aZMmaKJEydq06ZN3m2jRo1Senq6srOzm5SfPn260tLS9Pjjj3u3ZWZm6vDhwyooKJAkDRkyRA8++KCWLVvmLZOenq6YmBjt2LHDu+3LL7/UxIkTlZOTo4cffljjx4/Xk08+GfIxss48AAAAAKCthZo1Qx4z769l2x/DMEIqd+nSJRUXF2vVqlU+22fNmqXCwkK/+9TU1Cg6Otpnm91uV1FRkWpraxUZGRmwjCfseyxbtky33HKL/uVf/kUPP/xws/WtqalRTU2N9+eKiopm9wEAAAAAoD2EHOZ///vf69SpU0pKSlJExNX3zne5XKqrq1N8fLzP9vj4eJ09e9bvPrNnz9bzzz+v9PR0TZw4UcXFxdq2bZtqa2vlcrmUkJCg2bNna+PGjZoxY4ZGjBiht956S6+//rrq6uq8r7Nz50699957OnToUMj1zc7O1kMPPdS6gwUAAAAAoA21KJVfe+21crlc3p8zMjL05z//+aoq0Lgl3zTNgK37a9as0Zw5czR16lRFRkZq/vz5WrRokaQra9xL0lNPPaVrr71WKSkpioqK0vLly7V48WLv82fOnNG9996rHTt2NGnBDyYrK0tut9v7OHPmTCuOFgAAAACAq9eiMN94eP2ePXtUWVnZqjeOi4uTzWZr0gp/7ty5Jq31Hna7Xdu2bVNVVZVKS0tVVlampKQk9enTR3FxcZKkgQMH6rXXXlNlZaU+//xzffLJJ4qJiVFycrIkqbi4WOfOnVNqaqp69OihHj16aN++ffrpT3+qHj16+LTgN9SzZ0/Fxsb6PAAAAAAACIewzWYfFRWl1NRU5efn+2zPz8/X9OnTg+4bGRmpoUOHymazaefOnZo3b16Trv/R0dG65pprdPnyZe3atUvz58+XJP3zP/+zjhw5ovfff9/7mDRpku666y69//773hZ8AAAAAAA6q5DHzEtXusQ37gIf6oR3/qxYsUILFy7UpEmTNG3aNG3ZskVlZWVaunSppCtd28vLy71ryZ84cUJFRUWaMmWKLly4oI0bN+ro0aPavn279zUPHjyo8vJyjR8/XuXl5Vq7dq3q6+u1cuVKSVKfPn00ZswYn3r07t1bAwYMaLIdAAAAAIDOqEVh3jRNLVq0SD179pQkXbx4UUuXLlXv3r19yuXl5YX0ehkZGTp//rzWrVsnp9OpMWPGaM+ePRo+fLgkyel0+qw5X1dXpyeeeELHjx9XZGSkZs6cqcLCQiUlJXnLXLx4UatXr9apU6cUExOjuXPn6uWXX1bfvn1bcqgAAAAAAHRaIa8zL0mLFy8OqdyLL77Y6gpZBevMAwAAdA5Od7VOuyqVHNdbCQ57uKsDAFelzdeZl7pHSAcAAIB15B4qU1beEdWbUoQhZS8Yq4zJieGuFgC0u7BNgAcAAABcDae72hvkJanelB7IOyqnuzq8FQOADkCYBwAAgCWddlV6g7xHnWmq1FUVngoBQAcizAMAAMCSkuN6K6LRwko2w1BSXK/wVAgAOhBhHgAAAJaU4LAre8FY2f6+VLLNMLR+wRgmwQPQLbRoAjwAAACgM8mYnKgZ1w1UqatKSXG9CPIAug3CPAAAACwtwWEnxAPoduhmDwAAAACAxRDmAQAAAACwGMI8AAAAAAAWQ5gHAAAAAMBiCPMAAAAAAFgMYR4AAAAAAIshzAMAAAAAYDGEeQAAAAAALIYwDwAAAACAxRDmAQAAYDlOd7UKT7rkdFeHuyoAEBY9wl0BAAAAoCVyD5UpK++I6k0pwpCyF4xVxuTEcFcLADoULfMAAACwDKe72hvkJanelB7IO0oLPYBuhzAPAAAAyzjtqvQGeY8601Spqyo8FQKAMCHMAwAAwDKS43orwvDdZjMMJcX1Ck+FACBMCPMAAC8mlALQ2SU47MpeMFY240qitxmG1i8YowSHPcw1A4COxQR4AABJTCgFwDoyJidqxnUDVeqqUlJcL4I8gG6JlnkAABNKAbCcBIdd00YMIMgD6LYI8wAAJpQCAACwGMI8AIAJpQAAACyGMA8AYEIpAAAAi2ECPACAJCaUAgAAsBLCPADAK8FhJ8QDAABYAN3sAQAAAACwGMI8AAAAAAAWQ5gHAAAAAMBiCPMAAAAAAFgMYR4AAAAAAIshzAMAAAAAYDGEeQAAAAAALIYwDwAAAACAxRDmAQAAAACwGMI8AAAAAAAWQ5gHAAAAAMBiCPMAAACwDKe7WoUnXXK6q8NdFQAIqx7hrgAAAAAQitxDZcrKO6J6U4owpOwFY5UxOTHc1QKAsKBlHgDgF61fADoTp7vaG+Qlqd6UHsg7ymcUgG4r7GE+JydHycnJio6OVmpqqvbv3x+0/LPPPqtRo0bJbrdr5MiReumll3yer62t1bp16zRixAhFR0dr3Lhx2rt3r0+Z7OxsTZ48WX369NGgQYOUnp6u48ePt/mxAYBV5R4qU9qGt3Xn1oNK2/C2cg+VhbtKALq5065Kb5D3qDNNlbqqwlMhAAizsIb53NxcZWZm6sEHH1RJSYluvPFGzZkzR2Vl/r80btq0SVlZWVq7dq0++ugjPfTQQ1q2bJl+/etfe8usXr1amzdv1tNPP61jx45p6dKluu2221RSUuIts2/fPi1btkx//OMflZ+fr8uXL2vWrFmqrKxs92MGgM6O1i8AnVFyXG9FGL7bbIahpLhe4akQAISZYZqm2Xyx9jFlyhRNnDhRmzZt8m4bNWqU0tPTlZ2d3aT89OnTlZaWpscff9y7LTMzU4cPH1ZBQYEkaciQIXrwwQe1bNkyb5n09HTFxMRox44dfuvxl7/8RYMGDdK+ffs0Y8aMkOpeUVEhh8Mht9ut2NjYkPYBACsoPOnSnVsPNtn+iyVTNW3EgDDUCACuyD1UpgfyjqrONGUzDK1fMIYx8wC6nFCzZtgmwLt06ZKKi4u1atUqn+2zZs1SYWGh331qamoUHR3ts81ut6uoqEi1tbWKjIwMWMYT9v1xu92SpP79+wcsU1NTo5qaGu/PFRUVAcsCgJV5Wr8admel9QtAZ5AxOVEzrhuoUleVkuJ6KcFhD3eVACBswtbN3uVyqa6uTvHx8T7b4+PjdfbsWb/7zJ49W88//7yKi4tlmqYOHz6sbdu2qba2Vi6Xy1tm48aN+vTTT1VfX6/8/Hy9/vrrcjqdfl/TNE2tWLFCN9xwg8aMGROwvtnZ2XI4HN7HsGHDWnnkANC5JTjsyl4wVjbjSn9WT+sXX5oBdAYJDrumjRjAZxKAbi/sS9MZhu/gJ9M0m2zzWLNmjc6ePaupU6fKNE3Fx8dr0aJFeuyxx2Sz2SRJTz31lJYsWaKUlBQZhqERI0Zo8eLFevHFF/2+5vLly/Xhhx8GbbmXpKysLK1YscL7c0VFBYEeQJdF6xcAAEDnFraW+bi4ONlstiat8OfOnWvSWu9ht9u1bds2VVVVqbS0VGVlZUpKSlKfPn0UFxcnSRo4cKBee+01VVZW6vPPP9cnn3yimJgYJScnN3m9H/zgB/rVr36l3//+9xo6dGjQ+vbs2VOxsbE+DwDoymj9AgAA6LzCFuajoqKUmpqq/Px8n+35+fmaPn160H0jIyM1dOhQ2Ww27dy5U/PmzVNEhO+hREdH65prrtHly5e1a9cuzZ8/3/ucaZpavny58vLy9Pbbb/sN+gAAAAAAdFZh7Wa/YsUKLVy4UJMmTdK0adO0ZcsWlZWVaenSpZKudG0vLy/3riV/4sQJFRUVacqUKbpw4YI2btyoo0ePavv27d7XPHjwoMrLyzV+/HiVl5dr7dq1qq+v18qVK71lli1bpldeeUWvv/66+vTp4+0d4HA4ZLfTAgUAAAAA6NzCGuYzMjJ0/vx5rVu3Tk6nU2PGjNGePXs0fPhwSZLT6fRZc76urk5PPPGEjh8/rsjISM2cOVOFhYVKSkrylrl48aJWr16tU6dOKSYmRnPnztXLL7+svn37est4lsL7xje+4VOfF198UYsWLWqvwwUAAAAAoE2EdZ15K2OdeQAAAABAWws1a4ZtzDwAAAAAAGgdwjwAAAAAABZDmAcAAAAAwGII8wAAAAAAWAxhHgAAAJ2a012twpMuOd3V4a4KAHQaYV2aDgAAAAgm91CZsvKOqN6UIgwpe8FYZUxODHe1ACDsaJkHAABAp+R0V3uDvCTVm9IDeUdpoQcAEeYBAADQSZ12VXqDvEedaarUVRWeCgFAJ0KYBwA0i/GqAMIhOa63IgzfbTbDUFJcr/BUCAA6EcI8ACCo3ENlStvwtu7celBpG95W7qGycFcJQDeR4LAre8FY2Ywrid5mGFq/YIwSHPYw1wwAws8wTdNsvhgaq6iokMPhkNvtVmxsbLirAwDtwumuVtqGt326udoMQwWrZvJlGkCHcbqrVeqqUlJcLz57AHR5oWZNZrMHAAQUbLwqX6gBdJQEh53PHABohG72AICAGK8KAADQORHmAQABMV4VAACgc6KbPQAgqIzJiZpx3UDGqwIAAHQihHkAQLMajld1uqt12lWp5LjeBHsAAIAwIcwDAEKWe6hMWXlHVG9KEYaUvWCsMiYnhrtaAAAA3Q5j5gEAIXG6q71BXpLqTemBvKNyuqvDWzEAAIBuiDAPAAhJsGXqAKCtOd3VKjzp4oYhAARAN3sAQEg8y9Q1DPQsUwegPTCkBwCaR8s8ACAkLFMHoCMwpAcAQkPLPAAgZCxTB6C9BRvSw2cOAPwDYR4A0CINl6kDgLbkdFfr/Jc1DOkBgBAQ5gEAABB2DcfJG5IMQzJNhvQAQCCEeQAAAIRV43HypqQIU3rmzgmaOLwfQR4A/GACPABAq7F0FIC24G+cfL2k/r17EuQBIABa5gEArcLSUQDaCktfAkDL0TIPAGgxlo4C0Fac7mqddlXq/jkpLH0JAC1AyzwAoMVYOgpAW2jcw+f+m1N0/dC+LH0JACGgZR4A0GKeLrEN0SUWQEv46+Hz2N7jBHkACBFhHgDQYgkOu7IXjPV2iY2Q9N0bksJaJwDWEqyHDwCgeYR5AECrZExOVMGqmfr/3/gVyZC27D+ttA1vK/dQWbirBsAC6OEDAFeHMA8AuCrPF5xiIjwALda4hw+T3gFAyzABHgCg1QJ1k939oVO3XJ/Al3IAQWVMTtSM6waq1FXFWHkAaCFa5gEAreavm6wkPbz7Y7rcAwjK6a5W4UmXJGnaiAEEeQBoIVrmAQCt5ukm+0DeUdWZvk30ni73M64byJd0AJL+sab8kXK3Hn3jE++SdNkLxipjcmK4qwcAlkKYBwBcFU832d0fOvXw7o99nmPteQAeDdeUb4gbfwDQOnSzBwBctQSHXbdcn9Cky32EpPOVNUyIB3RzjdeUb4wl6QCg5QjzAIA20XhmakOSKWn5KyWMnwe6OX+TZTbEknQA0HJ0swcAtBlPl/vi0gu6Z2eJz5J1WbuOqHfPHkod3k/SlS/3yXG96VYLdAOeyTL9BXqWpAOA1iHMAwDaVILDrv4xTVvh6nWlld7TE98UE18B3YFn0rv756TosTeOq840ZTMMrbx5pK4f2pcl6QCglQjzAIA2F6wVruEmJr4CuraGk95FGNL9N6cQ4AGgjTBmHgDQ5hqPnw+Gia+ArqnxpHf1pvTY3uMEeQBoI4R5AEC7yJicqIJVM/XMHROazHLfEBNfAV2P012t33z4RZPeOdy8A4C2E/Ywn5OTo+TkZEVHRys1NVX79+8PWv7ZZ5/VqFGjZLfbNXLkSL300ks+z9fW1mrdunUaMWKEoqOjNW7cOO3du/eq3xcA0HIJDrvmjRvSZJZ7T4N9hKTv3pAUruoBaAe5h8qUtuFtPbL7kybPcfMOANpOWMfM5+bmKjMzUzk5OUpLS9PmzZs1Z84cHTt2TImJTSdD2rRpk7KysrR161ZNnjxZRUVFWrJkifr166dbb71VkrR69Wrt2LFDW7duVUpKit58803ddtttKiws1IQJE1r1vgCAq+OZ5b7UVeX9Iv9iQameLzilLftP6/mC00yEB3QBwdaTZ9Z6AGhbhmmaQVb9bF9TpkzRxIkTtWnTJu+2UaNGKT09XdnZ2U3KT58+XWlpaXr88ce92zIzM3X48GEVFBRIkoYMGaIHH3xQy5Yt85ZJT09XTEyMduzY0ar3laSamhrV1NR4f66oqNCwYcPkdrsVGxvbyt8AAHRPTne10ja87fOF32YYKlg1ky/6gIUVnnTpzq0Hm2xfc8sozb0+gfMbAEJQUVEhh8PRbNYMWzf7S5cuqbi4WLNmzfLZPmvWLBUWFvrdp6amRtHR0T7b7Ha7ioqKVFtbG7SMJ+y35n0lKTs7Ww6Hw/sYNmxYaAcKAGjitKvp0nWMpQWsz7OSRUM2wyDIA0A7CFuYd7lcqqurU3x8vM/2+Ph4nT171u8+s2fP1vPPP6/i4mKZpqnDhw9r27Ztqq2tlcvl8pbZuHGjPv30U9XX1ys/P1+vv/66nE5nq99XkrKysuR2u72PM2fOXM3hA0C3FugLP2NpAWtrvJIFXesBoP2EfZ15o9GyRaZpNtnmsWbNGp09e1ZTp06VaZqKj4/XokWL9Nhjj8lms0mSnnrqKS1ZskQpKSkyDEMjRozQ4sWL9eKLL7b6fSWpZ8+e6tmzZ2sOEQDQiOcL/wN5R1VnmrIZhlbePFKnXZXeMqddlUqO600IACym8RwZnMMA0D7CFubj4uJks9matIafO3euSau5h91u17Zt27R582b9+c9/VkJCgrZs2aI+ffooLi5OkjRw4EC99tprunjxos6fP68hQ4Zo1apVSk5ObvX7AgDaXsMv/B+W/1WPvvGJ6s0rs91LkikpwhAT4wEWlOCwE+IBoJ2FrZt9VFSUUlNTlZ+f77M9Pz9f06dPD7pvZGSkhg4dKpvNpp07d2revHmKiPA9lOjoaF1zzTW6fPmydu3apfnz51/1+wIA2laCw66kuF7eIC9dCfGe4fT1pvRA3lE53dXhqiIAAECnFNZu9itWrNDChQs1adIkTZs2TVu2bFFZWZmWLl0q6co49fLycu9a8idOnFBRUZGmTJmiCxcuaOPGjTp69Ki2b9/ufc2DBw+qvLxc48ePV3l5udauXav6+nqtXLky5PcFAHQcf5PhNVRnmtr9oVO3MIEWAACAV1jDfEZGhs6fP69169bJ6XRqzJgx2rNnj4YPHy5JcjqdKisr85avq6vTE088oePHjysyMlIzZ85UYWGhkpKSvGUuXryo1atX69SpU4qJidHcuXP18ssvq2/fviG/LwCg43gmwwsW6B/e/bHW7/mYLvcAAAB/F9Z15q0s1LX/AADNyz1U5p0Mz5AkQ/J3dWItegAA0NWFmjXDPps9AACNZ7+WpN0fOvXw7o99ynnWoifMAwCA7i5sE+ABANBQgsOuaSMGeGfBvuX6hIBr0Tvd1So86Qp5YryWlgcAAOjsaJkHAHRKgdai31ZwWi8UnFa9GdrSdbmHypSVdyTk8gAAAFbAmPlWYsw8AHQMp7vauxb9hj2fqPFFK0LST++coNTh/Zp0v3e6q5W24W2fyfUYdw8AADqzULMm3ewBAJ1aw7Xo/d19rpe0/JUSpW14W7mHynye87fsnWfcPQAAgJUR5gEAnV5za9FLV5a2eyDvqHdcvNNdrfNf1gQcdw8AAGBljJkHAHR6oaxFL11pdd/9oVN1pqlH3/hE9aZkSDL+vtSdzTC0fsEYutgDAADLY8x8KzFmHgA6VsO16CMk3T5lmHYWnWk24EtXuqE9fecETfQzrh4AAKAzYZ15AECX0ngt+gSHXeOG9vUG/GDqJfXv3ZMgDwAAugzCPADAMjxr0Ht4Av7uD516ePfHAffzN07e6a7WaVelkuN6S5L33wR+AABgBYR5AIClJTjsuuX6BK3f87HfLvf+xsk3XHveMz+eKdahBwAA1kGYBwBYXoLDruwFY71d7m2GoZU3j9T1Q/t6u+R7ON3V3iAvyWe5u3pTytp1RL179vC7bj0AAEBnQZgHAHQJ/sbUN+Z0V+s3H34RdNI8z7r1tNIDAIDOjDAPAOgyGo+pb6hh1/pQ0EoPAAA6s4hwVwAAgPbWuGt9Q5516P3xtNKnbXhbuYfK2rOKAAAALULLPACgyzvtqvQb5NfcMkpzr0+QJBWXXtA9O0v8lvO00qcM7qNBsdHMfA8AAMKOMA8A6PKS43orwpBPULcZhuZen+AN5PPG2VV56XLAdevrJc1/tlCG/jHz/f03p2jsUAfBHgAAdDjDNP18Y0GzKioq5HA45Ha7FRsbG+7qAACakXuozGe2+/ULxvid3M7prg7aSu9PSyfLa7jGfWe7CdCZ6wYAQHcQatakZR4A0C2EMtu9dGUSPU8rfdauI6oP4bVbMllew4n4OlvrfuO6MZs/AACdFy3zrUTLPAB0fR+cuaD0nEK15ErZMKD3jrKp8lKdN6g73dVK2/C23xb/cIdnf3WzGYYKVs0M+00GAAC6E1rmAQC4SuOG9dOGBWO93fMNSTIUNNzXm1L2G5/4bPME/IgIBey6X29KD+Qd1YzrBoYlPPubJLDONFXqqiLMAwDQCRHmAQAIonH3fEkqdVXpw/K/6rE3jvudLK8xfwHfnzrT1O4PnbqlwcR8DcewS2q38eyBJgn0HDMAAOhc6GbfSnSzBwC0ZrK8UHi63EvyjmE3/v5cczPp+5vALtRJ7UKdJBAAALSfULMmYb6VCPMAAI+GIThUa24ZpcumGbB1v2F4D6bhWHt/E9hJCmlSO0/g7x1lU9Wl+qCTBAIAgPZDmG9nhHkAQENOd3XI3e8bTizndFdr94dOPbz741a/t80wlPf9abotp9Cnh4C/Mf4Rkn565wSfWfeZxR4AgM6DMN/OCPMAgEA8wb5XVISqLtX7BHx/3deDzXIfqn+dkqgdB8tCLu8J7TOuGxjSLPasPw8AQMcgzLczwjwAoCU8AT9Q9/XmuuqHMpN+S9kMQ0/dMV7LXylp8twvlkzVtBEDvHWj5R4AgI7B0nQAAHQiCQ570BZtz6z5/ibUi5D06rLpGhQb3eKZ9IOpM02dOPu3oLPYO93V3iAvhX8JPQAAcEVEuCsAAACuSHDYNW/cEGUvGCubcWUKPJthKPtbYzVu2JUx7tNGDND3ZoxQwaqZWn3LKL+vYzT6OULSw+mjFdH4CUk/ffszmab097fzDgPwBPVg688DAIDwoWUeAIBOpvHa9v5awBMcdt1yfYLW7/m4Sav6yjkj/Y7Rj7RF+O3Kb0qKMKVn7pygiQ0mxpNYfx4AgM6KMA8AQCfUXLd8T5nsBWP9rg3/zXFDmtwM8Nwk8Dd7fr2k/r17+p307v45KU1uDtDFHgCA8CLMAwBgYYFa8QPdDAjWot+wtb3xpHf335yi64f2Zf15AAA6CcbMAwBgcZ6x9KGGbE+LfsNx+Q1b2/1NevfY3uMEeQAAOhFa5gEA6IaCjcsPNukdYR4AgM6BMA8AQDcVqCs+k94BAND50c0eAAD4aK4bPgAACD9a5gEAQBOhLI8HAADChzAPAAD8CmV5PAAAEB50swcAAAAAwGII8wAAAAAAWAxhHgAAAAAAiyHMAwAAAABgMYR5AAAAAAAsJuxhPicnR8nJyYqOjlZqaqr2798ftPyzzz6rUaNGyW63a+TIkXrppZealHnyySc1cuRI2e12DRs2TPfdd58uXrzoff7y5ctavXq1kpOTZbfb9ZWvfEXr1q1TfX19mx8fAAAAAABtLaxL0+Xm5iozM1M5OTlKS0vT5s2bNWfOHB07dkyJiYlNym/atElZWVnaunWrJk+erKKiIi1ZskT9+vXTrbfeKkn6+c9/rlWrVmnbtm2aPn26Tpw4oUWLFkmS/uu//kuS9Oijj+q5557T9u3bNXr0aB0+fFiLFy+Ww+HQvffe22HHDwAAAABAaximaZrhevMpU6Zo4sSJ2rRpk3fbqFGjlJ6eruzs7Cblp0+frrS0ND3++OPebZmZmTp8+LAKCgokScuXL9fHH3+st956y1vm3//931VUVORt9Z83b57i4+P1wgsveMt861vfUq9evfTyyy+HVPeKigo5HA653W7Fxsa27MABAAAAAPAj1KwZtm72ly5dUnFxsWbNmuWzfdasWSosLPS7T01NjaKjo3222e12FRUVqba2VpJ0ww03qLi4WEVFRZKkU6dOac+ePbrlllu8+9xwww166623dOLECUnSBx98oIKCAs2dOzdgfWtqalRRUeHzAAAAAAAgHMLWzd7lcqmurk7x8fE+2+Pj43X27Fm/+8yePVvPP/+80tPTNXHiRBUXF2vbtm2qra2Vy+VSQkKCbr/9dv3lL3/RDTfcINM0dfnyZf3f//t/tWrVKu/r3H///XK73UpJSZHNZlNdXZ0eeeQR3XHHHQHrm52drYceeqhtDh4AAAAAgKsQ1jHzkmQYhs/Ppmk22eaxZs0anT17VlOnTpVpmoqPj9eiRYv02GOPyWazSZLeeecdPfLII8rJydGUKVP02Wef6d5771VCQoLWrFkj6cpY/R07duiVV17R6NGj9f777yszM1NDhgzR3Xff7fe9s7KytGLFCu/PbrdbiYmJtNADAAAAANqMJ2M2OyLeDJOamhrTZrOZeXl5Ptvvuecec8aMGUH3vXTpknnmzBnz8uXLZk5OjtmnTx+zrq7ONE3TvOGGG8z/+I//8Cn/8ssvm3a73Vtm6NCh5jPPPONT5j//8z/NkSNHhlz/M2fOmJJ48ODBgwcPHjx48ODBgwePNn+cOXMmaCYNW8t8VFSUUlNTlZ+fr9tuu827PT8/X/Pnzw+6b2RkpIYOHSpJ2rlzp+bNm6eIiCvD/6uqqrz/9rDZbDJN03tnI1CZlixNN2TIEJ05c0Z9+vQJ2JOgq6ioqNCwYcN05swZJvsD/OAcAYLjHAGax3kCBNedzhHTNPW3v/1NQ4YMCVourN3sV6xYoYULF2rSpEmaNm2atmzZorKyMi1dulTSla7t5eXl3rXkT5w4oaKiIk2ZMkUXLlzQxo0bdfToUW3fvt37mrfeeqs2btyoCRMmeLvZr1mzRt/85je9XfFvvfVWPfLII0pMTNTo0aNVUlKijRs36jvf+U7IdY+IiPDeUOguYmNju/yJA1wNzhEgOM4RoHmcJ0Bw3eUccTgczZYJa5jPyMjQ+fPntW7dOjmdTo0ZM0Z79uzR8OHDJUlOp1NlZWXe8nV1dXriiSd0/PhxRUZGaubMmSosLFRSUpK3zOrVq2UYhlavXq3y8nINHDjQG949nn76aa1Zs0bf//73de7cOQ0ZMkTf+9739KMf/ajDjh0AAAAAgNYK6zrzsIZQ1zkEuivOESA4zhGgeZwnQHCcI02FbZ15WEfPnj314x//WD179gx3VYBOiXMECI5zBGge5wkQHOdIU7TMAwAAAABgMbTMAwAAAABgMYR5AAAAAAAshjAPAAAAAIDFEOYBAAAAALAYwjx04cIFLVy4UA6HQw6HQwsXLtRf//rXoPvk5eVp9uzZiouLk2EYev/995uUqamp0Q9+8APFxcWpd+/e+uY3v6k//elP7XMQQDtqzTlimqbWrl2rIUOGyG636xvf+IY++ugjnzLf+MY3ZBiGz+P2229vxyMB2k5OTo6Sk5MVHR2t1NRU7d+/P2j5ffv2KTU1VdHR0frKV76i5557rkmZXbt26atf/ap69uypr371q3r11Vfbq/pAu2vrc+RnP/tZk2uGYRi6ePFiex4G0G5aco44nU7deeedGjlypCIiIpSZmem3XHe7jhDmoTvvvFPvv/++9u7dq7179+r999/XwoULg+5TWVmptLQ0bdiwIWCZzMxMvfrqq9q5c6cKCgr05Zdfat68eaqrq2vrQwDaVWvOkccee0wbN27UM888o0OHDmnw4MG66aab9Le//c2n3JIlS+R0Or2PzZs3t+ehAG0iNzdXmZmZevDBB1VSUqIbb7xRc+bMUVlZmd/yp0+f1ty5c3XjjTeqpKREDzzwgO655x7t2rXLW+bAgQPKyMjQwoUL9cEHH2jhwoX69re/rYMHD3bUYQFtpj3OEUmKjY31uWY4nU5FR0d3xCEBbaql50hNTY0GDhyoBx98UOPGjfNbplteR0x0a8eOHTMlmX/84x+92w4cOGBKMj/55JNm9z99+rQpySwpKfHZ/te//tWMjIw0d+7c6d1WXl5uRkREmHv37m2z+gPtrTXnSH19vTl48GBzw4YN3m0XL140HQ6H+dxzz3m3ff3rXzfvvffedqs70F6+9rWvmUuXLvXZlpKSYq5atcpv+ZUrV5opKSk+2773ve+ZU6dO9f787W9/27z55pt9ysyePdu8/fbb26jWQMdpj3PkxRdfNB0OR5vXFQiHlp4jDQX6/tQdryO0zHdzBw4ckMPh0JQpU7zbpk6dKofDocLCwla/bnFxsWprazVr1izvtiFDhmjMmDFX9bpAR2vNOXL69GmdPXvW5++/Z8+e+vrXv95kn5///OeKi4vT6NGj9R//8R9NWu6BzubSpUsqLi72+fuWpFmzZgU8Jw4cONCk/OzZs3X48GHV1tYGLcM1A1bTXueIJH355ZcaPny4hg4dqnnz5qmkpKTtDwBoZ605R0LRHa8jPcJdAYTX2bNnNWjQoCbbBw0apLNnz17V60ZFRalfv34+2+Pj46/qdYGO1ppzxLM9Pj7eZ3t8fLw+//xz78933XWXkpOTNXjwYB09elRZWVn64IMPlJ+f34ZHALQtl8uluro6v3/fwc4Jf+UvX74sl8ulhISEgGW4ZsBq2uscSUlJ0c9+9jONHTtWFRUVeuqpp5SWlqYPPvhA1157bbsdD9DWWnOOhKI7Xkdome+i1q5d63eSlIaPw4cPS5IMw2iyv2mafrdfrfZ6XaClOuIcafx8432WLFmif/mXf9GYMWN0++2365e//KV+97vf6b333muDIwTaV3N/36GUb7y9pa8JdGZtfY5MnTpV//qv/6px48bpxhtv1P/7f/9P1113nZ5++uk2rjnQMdrjM7+7XUdome+ili9f3uys2ElJSfrwww/15z//uclzf/nLX5rc2WqJwYMH69KlS7pw4YJP6/y5c+c0ffr0Vr8u0Fba8xwZPHiwpCt3iBMSErzbz507F/S8mjhxoiIjI/Xpp59q4sSJoRwG0OHi4uJks9matHQE+/sePHiw3/I9evTQgAEDgpa5mmsREA7tdY40FhERocmTJ+vTTz9tm4oDHaQ150gouuN1hJb5LiouLk4pKSlBH9HR0Zo2bZrcbreKioq8+x48eFBut/uqQndqaqoiIyN9ugs7nU4dPXqUMI9OoT3PEU/X+YZ//5cuXdK+ffuC/v1/9NFHqq2t9bkBAHQ2UVFRSk1NbTIcJD8/P+Df97Rp05qU/+1vf6tJkyYpMjIyaBmuGbCa9jpHGjNNU++//z7XDFhOa86RUHTL60h45t1DZ3LzzTeb119/vXngwAHzwIED5tixY8158+b5lBk5cqSZl5fn/fn8+fNmSUmJuXv3blOSuXPnTrOkpMR0Op3eMkuXLjWHDh1q/u53vzPfe+8985/+6Z/McePGmZcvX+6wYwPaQmvOkQ0bNpgOh8PMy8szjxw5Yt5xxx1mQkKCWVFRYZqmaX722WfmQw89ZB46dMg8ffq0uXv3bjMlJcWcMGEC5wg6vZ07d5qRkZHmCy+8YB47dszMzMw0e/fubZaWlpqmaZqrVq0yFy5c6C1/6tQps1evXuZ9991nHjt2zHzhhRfMyMhI85e//KW3zB/+8AfTZrOZGzZsMD/++GNzw4YNZo8ePXxWkgCsoj3OkbVr15p79+41T548aZaUlJiLFy82e/ToYR48eLDDjw+4Wi09R0zTNEtKSsySkhIzNTXVvPPOO82SkhLzo48+8j7fHa8jhHmY58+fN++66y6zT58+Zp8+fcy77rrLvHDhgk8ZSeaLL77o/fnFF180JTV5/PjHP/aWqa6uNpcvX27279/ftNvt5rx588yysrKOOSigDbXmHKmvrzd//OMfm4MHDzZ79uxpzpgxwzxy5Ij3+bKyMnPGjBlm//79zaioKHPEiBHmPffcY54/f76Djgq4Os8++6w5fPhwMyoqypw4caK5b98+73N33323+fWvf92n/DvvvGNOmDDBjIqKMpOSksxNmzY1ec3//u//NkeOHGlGRkaaKSkp5q5du9r7MIB209bnSGZmppmYmGhGRUWZAwcONGfNmmUWFhZ2xKEA7aKl54i/7DF8+HCfMt3tOmKY5t9n1wAAAAAAAJbAmHkAAAAAACyGMA8AAAAAgMUQ5gEAAAAAsBjCPAAAAAAAFkOYBwAAAADAYgjzAAAAAABYDGEeAAAAAACLIcwDAAAAAGAxhHkAANDE2rVrNX78+A5/33feeUeGYcgwDKWnpwct+41vfEOZmZk+P3v2ff/999u1ngAAhFuPcFcAAAB0LMMwgj5/991365lnntEPfvCDDqpRU8ePH9egQYNatE9eXp5Onjypr33ta+1UKwAAOg/CPAAA3YzT6fT+Ozc3Vz/60Y90/Phx7za73a6YmBjFxMSEo3qSpEGDBqlv374t2qd///6qqKhonwoBANDJ0M0eAIBuZvDgwd6Hw+GQYRhNtjXuZr9o0SKlp6dr/fr1io+PV9++ffXQQw/p8uXL+uEPf6j+/ftr6NCh2rZtm897lZeXKyMjQ/369dOAAQM0f/58lZaWtrjOlZWV+j//5/8oJiZGCQkJeuKJJ67ytwAAgLUR5gEAQEjefvttffHFF3r33Xe1ceNGrV27VvPmzVO/fv108OBBLV26VEuXLtWZM2ckSVVVVZo5c6ZiYmL07rvvqqCgQDExMbr55pt16dKlFr33D3/4Q/3+97/Xq6++qt/+9rd65513VFxc3B6HCQCAJRDmAQBASPr376+f/vSnGjlypL7zne9o5MiRqqqq0gMPPKBrr71WWVlZioqK0h/+8AdJ0s6dOxUREaHnn39eY8eO1ahRo/Tiiy+qrKxM77zzTsjv++WXX+qFF17QT37yE910000aO3astm/frrq6unY6UgAAOj/GzAMAgJCMHj1aERH/aAeIj4/XmDFjvD/bbDYNGDBA586dkyQVFxfrs88+U58+fXxe5+LFizp58mTI73vy5EldunRJ06ZN827r37+/Ro4c2dpDAQDA8gjzAAAgJJGRkT4/G4bhd1t9fb0kqb6+Xqmpqfr5z3/e5LUGDhwY8vuaptmK2gIA0LUR5gEAQLuYOHGicnNzNWjQIMXGxrb6df7X//pfioyM1B//+EclJiZKki5cuKATJ07o61//eltVFwAAS2HMPAAAaBd33XWX4uLiNH/+fO3fv1+nT5/Wvn37dO+99+pPf/pTyK8TExOj7373u/rhD3+ot956S0ePHtWiRYt8uvwDANDd0DIPAADaRa9evfTuu+/q/vvv14IFC/S3v/1N11xzjf75n/+5xS31jz/+uL788kt985vfVJ8+ffTv//7vcrvd7VRzAAA6P8NkIBoAAOgk3nnnHc2cOVMXLlxQ3759W7x/aWmpkpOTVVJSovHjx7d5/QAA6CzonwYAADqdoUOH6o477mjRPnPmzNHo0aPbqUYAAHQutMwDAIBOo7q6WuXl5ZKujJUfPHhwyPuWl5erurpakpSYmKioqKh2qSMAAJ0BYR4AAAAAAIuhmz0AAAAAABZDmAcAAAAAwGII8wAAAAAAWAxhHgAAAAAAiyHMAwAAAABgMYR5AAAAAAAshjAPAAAAAIDFEOYBAAAAALCY/w/Rrvun73y7NAAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "time, flux = lcs(radius_ratio=0.1, zero_epoch=0.0, period=2.0,\n", " scaled_semi_major_axis=8.0, impact_parameter=0.1,\n", " efficiency=0.5, eff_phase=0.2,\n", " limb_darkening=[0.2, 0.3])\n", "\n", "fig, ax = plt.subplots(1, 1, figsize=(10,4), constrained_layout=True)\n", "ax.plot(time, flux, '.')\n", "plt.setp(ax, xlabel='Time [d]', ylabel='Flux');" ] }, { "cell_type": "markdown", "id": "b0f3e9fed1eab871", "metadata": {}, "source": [ "## How to simulate a full phase curve\n", "\n", "A phase curve with a secondary eclipse can be simulated by giving the simulator a non-zero `geometric_albedo` argument. Let's create a long observation of 1.5 days with observing efficiency of 0.6, and also set a non-zero orbital ``eccentricity`` and ``argument_of_periastron`` to create a more interesting light curve (both eccentricity and argument of periastron can be set for simple transit simulations too, but their effect is very small)." ] }, { "cell_type": "code", "execution_count": 6, "id": "a65a637ab74612d1", "metadata": { "ExecuteTime": { "end_time": "2025-03-04T14:07:58.873660Z", "start_time": "2025-03-04T14:07:58.749365Z" } }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAA/MAAAGbCAYAAACIxMC9AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8hTgPZAAAACXBIWXMAAA9hAAAPYQGoP6dpAACMLElEQVR4nO3dd3gUVdsG8HtLeiWEJISEDoHQCRCKNKkKggICorHziagIKCpiwfKKWBEUeFUEfUVBKVakSYfQQu8ttJBOek92vj+Snewmu5tsyM7MZu/fdeWCnZ3ZOZvJzs4z5znPUQmCIICIiIiIiIiI7IZa7gYQERERERERkXUYzBMRERERERHZGQbzRERERERERHaGwTwRERERERGRnWEwT0RERERERGRnGMwTERERERER2RkG80RERERERER2hsE8ERERERERkZ1hME9ERERERERkZxjMExEREREREdkZBvN27uuvv8aAAQPg7e0NlUqF9PR0uZtERERERERENsZg3g4MGDAAK1asMPlcbm4uhg8fjtdff13aRhEREREREZFstHI3gO7M9OnTAQA7duyQtR1EREREREQkHfbMExEREREREdkZBvNEREREREREdobBvAJ98MEH8PT0FH92796NKVOmVFpGREREREREjolj5hVoypQpGD9+vPj44YcfxtixYzFmzBhxWaNGjeRoGhERERERESkAg3kF8vPzg5+fn/jYzc0NAQEBaNmypYytIiIiIiIiIqVgMG/nEhISkJCQgEuXLgEATp48CS8vLzRu3NjohgARERERERHVHRwzb+eWLl2KLl26YPLkyQCAfv36oUuXLvjjjz9kbhkRERERERHZikoQBEHuRhARERERERFR9bFnnoiIiIiIiMjOMJgnIiIiIiIisjMsgKcgOp0Ot27dgpeXF1QqldzNISIiIiIiIokJgoCsrCwEBwdDrTbf/85gXkFu3bqF0NBQuZtBREREREREMrtx4wZCQkLMPs9gXkG8vLwAlB40b29vmVtDREREREREUsvMzERoaKgYH5rDYF5B9Kn13t7eDOaJiIiIiIgcWFVDr1kAj4iIiIiIiMjOMJgnIiIiIiIisjMM5omIiIiIiIjsDIN5IiIiIiIiIjvDYJ6IiIiIiIjIzjCYJyIiIiIiIrIzDOaJiIiIiIiI7AyDeSIiIiIiIiI7w2CeiIiIiIiIyM4wmCciIiIiIiKyMwzmiYiIiIiIiOwMg3kiIiIiIiIiO8NgnoiIiIiIiMjOMJgnIiIiIiIisjMM5omIiIiIiIjsDIN5IiIiIiIiIjvDYJ6IiIiIiIjIzjCYJyIiIiIiIrIzDOaJiIiIiIiI7AyDeSIiIiIiIiI7w2CeiIiIiIiIyM4wmCciIiIiIiKyMwzmiYiIiIiIiOwMg3kiIiIiIiIiO8NgnoiIiIiIiMjOMJgnIiIiIiIisjN2F8wvXrwYzZo1g6urKyIiIrB7926L6+/cuRMRERFwdXVF8+bNsXTp0krrrF27FuHh4XBxcUF4eDjWr19/R/t95plnoFKpsGDBAqvfHxEREREREVFV7CqYX716NaZPn445c+bg6NGj6Nu3L+655x5cv37d5PqxsbG499570bdvXxw9ehSvv/46pk2bhrVr14rrREdHY8KECYiKisLx48cRFRWF8ePH48CBAzXa72+//YYDBw4gODi49n8BRERERERERABUgiAIcjeiuiIjI9G1a1csWbJEXNa2bVvcf//9mDdvXqX1X331Vfzxxx84e/asuGzKlCk4fvw4oqOjAQATJkxAZmYm/vnnH3Gd4cOHo169evj555+t2m9cXBwiIyOxadMmjBgxAtOnT8f06dOr/f4yMzPh4+ODjIwMeHt7V3s7IiIiIiIiqhuqGxfaTc98YWEhYmJiMHToUKPlQ4cOxb59+0xuEx0dXWn9YcOG4fDhwygqKrK4jv41q7tfnU6HqKgozJo1C+3atavWeyooKEBmZqbRDxEREREREVFV7CaYT0lJQUlJCQIDA42WBwYGIiEhweQ2CQkJJtcvLi5GSkqKxXX0r1nd/c6fPx9arRbTpk2r9nuaN28efHx8xJ/Q0NBqb0tERERERESOy26CeT2VSmX0WBCESsuqWr/i8uq8pqV1YmJi8MUXX2DFihUW21LR7NmzkZGRIf7cuHGj2tsSERERERGR47KbYN7f3x8ajaZSL3xSUlKlXnO9oKAgk+trtVrUr1/f4jr616zOfnfv3o2kpCQ0btwYWq0WWq0W165dw0svvYSmTZuafU8uLi7w9vY2+iEiIiIiIiKqit0E887OzoiIiMCWLVuMlm/ZsgW9e/c2uU2vXr0qrb9582Z069YNTk5OFtfRv2Z19hsVFYUTJ07g2LFj4k9wcDBmzZqFTZs21fxNExEREREREZmglbsB1pg5cyaioqLQrVs39OrVC19//TWuX7+OKVOmAChNW4+Li8MPP/wAoLRy/ZdffomZM2di8uTJiI6OxrJly8Qq9QDw4osvol+/fpg/fz5Gjx6N33//HVu3bsWePXuqvd/69euLPf16Tk5OCAoKQlhYmK1/LURERERERORg7CqYnzBhAlJTU/Huu+8iPj4e7du3x4YNG9CkSRMAQHx8vNHc782aNcOGDRswY8YMfPXVVwgODsbChQsxduxYcZ3evXtj1apVeOONN/Dmm2+iRYsWWL16NSIjI6u9XyIiIiIiIiIp2dU883Ud55knIiIiIiJybHVunnkiIiIiIiIiKsVgnoiIiIiIiMjOMJgnIiIiIiIisjMM5omIiIiIiIjsDIN5IiIiIiIiIjvDYJ6IiIiIiIjIzjCYJyIiIiIiIrIzDOaJiIiIiIiI7AyDeSIiIiIiIiI7w2CeiIiIiIiIyM4wmCciIiIiIiKyMwzmiYiIiIiIiOwMg3kiIiIiIiIiO8NgnoiIiIiIiMjOMJgnIiIiIiIisjMM5omIiIiIiIjsDIN5IiIiIiIiIjvDYJ6IiGpdUYkOeYUlcjeDiIiIqM5iME9ERLVu5MI96PzuZuQUFMvdFCIiIqI6icE8ERHVuvOJWSgo1uH4zXS5m0JERERUJzGYJyKiWlVcohP/X6ITZGwJEREpjU4nYPWh67iUlCV3U4jsHoN5IiKqVfnF5cF8sZlgPrugGL8evoG0nEKpmkVEZPf+ORmPNTE3AQBn4zPx1IpDOH0rQ+ZWWefvk/F4de1JDP5sl9xNIbJ7DOaJiKhW5ReVF757YvkhFBn01Os9t/IIZq05gU82n5eyaUREdiu/qATPrjyCl389jts5hXh17Qn8ey4JIxbuMcqIUqr9V1Kx9UwizsZnistmrzuJg7G3ZWwVkX1jME9EJDFBEKCrw+nnhsE8ALy65oTR41NxGdh5IRkAsPLAdcnaRURkz26m5Yr/7/reFpy4Wd4j/+KqYygsVl5AfyExCwu2XkBWfhEmfr0fT/9wGIt3XBaf//ngdYz/bzQA1OnvRSJb0crdACIiR/PIsgNIzS7EXy/cBa2m7t1TzS8yvqBcdzQOnq5aPNa7KVo08MTIRXuMni8s1sFZW/d+D0REtel8QrbZ5/4+GQ93Zw0+frCThC2q2tDPS1PpY1NyLK539yc7kJiZj0mRjTFnRLgUTSOqE3j1REQkkYSMfLz1+ynsvZSKcwlZuH47t+qN7FDFnnkA+CH6GgZ9utNkz9GpWxnYeCoegsBeGSIiU7afS8Ira45bXOfXsrH0SrTpdILF56+k5CCnsATf7I6VqEVEdQN75omIJDJ99VHsv1L3xwYWFFcO5vVe+PlIpWVjFu8DALw5MhxP3dXMZu0iIrI3giBgxb6reO+vM7DnLPSKGVtEVDvYM09EZEM6nSAWJqoYyGfkFeFiYpZdFC6yhqWLtk2nE80+995fZ7D+qHJ7loiIpFRUosObv5/CO3+WBvITuoVaXN/DWSNRy2xLP3Y+5tptzFh9DNdT62YWG1FtYM88EZGNCIKAB5bsQ1Z+ETZN71fp+XkbzuHg1dsI9nHF1IEtMS4iBK5O9n8xZirNvrpmrD4OnQ4Y0bEh/hd9DYPDA9HM36MWW0dEpHwZeUV4/qcj2H0xBSoV8Po9bfF032ZYffiG2W2K7Lnr3sCibZdQWFKCk3GZ2HUhGRtOxuP8+/fI3SwiRWIwT0RkI4mZBTh+Ix0AsOdSSqXnD14t7am/lZGPN347hYX/XsT/9WuOSZGN4e5sv6fngjusqPzSr8dxITEL/911BZ9tuYCz7w2vpZYRESnftdQcPLniEC4n58DNSYMvJnbG0HZBVW5XWKxDiU6ARq2SoJW28/nWC0aPC4p1yC8qqRM3u4lqG9PsiYhs5FJSeeXhJ5YfMrveoDYBaOjjiqSsArz/91n0+XAbFv57ERl5RVI0s9YlZeZXe90RHRqaXP5D9DUAQN4d9PITEckhKTMfD3+7H3+fiLd624Oxt3H/V3txOTkHDX1c8euUXlUG8ntfu1v8v6WaJVKrzbZcTDRfyZ/IkTGYJyKywqtrTmDckn0oqsY4958OXqu0rGl990rLZgxpjZ2zBmL+2A5oWt8dablF+GzLBfT5cBvmbzyHlOyCWmm7VC4lG190Pd67qdl1Q/0q/z4A00H8rgvJmL3upN39Poio7kvJLsCbv53C6VsZeP7no9h7KRXP/XQEMdfKa6XkFBRjxMLdmPvHaZOvsTbmJh7+dj/ScovQMcQHvz/XB+0b+Vjcb2M/dzT0dhUf5xUqJ5jPzi82emzpu6AqPO8TmWa/eZxERBLLLyoRxyueuZWJTqG+KCguwdIdVzAgrAE6hfqK6+67lIINJytPxRMe7I2rFYr51Pd0hrNWjQndG2Ns1xD8fTIei7dfxvnELCzZcRnL98ZiYvfG+L9+zRHs62bT91gb9D0oHRr5QKtRYcaQ1ni6bzMkZhZg7JJ9RusObhuApTsvV+t1P950HifjMnAmPhO/P9en1ttNRFRdOp0AlQpQqUpT2uesP4lNpxPxv/3GN3HHLokGAIT6ueHG7TwAwOlbmZg7qp3Ra3265Ty+2l56Lry3QxA+fbAz3CwUtGvg5YIvH+qClgGeUKtVcNGqUVCsU1Q2U3ZBeTD/7aPdMCCsAcZFhODQ1dt4588zVr0Wg3ki09gzT0RUTVeSc8T/q8su4P4XfQ2fb72A0V/tFedJPxWXgUnfHjD5GoXFlQsU1XN3Fv+v1agxunMj/PNiX3wdFYFOIT7IL9Jhxb6r6P/xdry65gSupuRUeg0luVzWM/+fB9pj/dQ+8HFzQkg9d0Q0qYfNM/oZVWTu1tSvyvGd647cRFGJDifjMgBArENARCQ1QRBQVKLDqK/2YNzSaBSV6JBbWIzjNzIsbqcP5A1fByjtSX/upyNiIP/8wJb48qGuZgP5LyZ2hr+nM5Y+0hWRzeujvqcLAIjjyZU0BVxWWc98oLcLBocHQqtRo30jHzzRp5lRltrgtoFVvlZqTqHN2klkz9gzT0RUDRcSs/DeX+U9CbmFxcgvKsEXWy+Ky3p88C+eH9gSb5tJoQRQKT3f39PFZFEftVqFoe2CMCQ8EHsvpeLL7Rex/8ptrD58A7/G3MDIjsF4bmBLhAV51cK7qz1pOYVIyS696GrRwLPS860DveDv5Wy0bOkjEZj8w2Gzrznzl+N463fzv1MiIqlM/Ho/TsZlILcsnb3Hf7YiLdf6+ia5hSXIKSjG0z8cxombGXDWqPHh2A4Y0zXE4najOzfCqE7BYkaAnpuTBhl5RXc0m0ht0/fMe7pUDjdyDIYDfPtYNwiCgGazN4jLHurRGD8fvC4+/udkPO7rFIxGdpCdRiQl9swTEVXD0M93IfpKqvj43b/OoM2bG5FlkEaYnFVgMZAHSqsNdzZIx58U2dji+iqVCne18seq/+uFtc/2wsCwBtAJwB/Hb2HYgl2Y/MNh3ErPs/gaUtKPl2/k6wYPExdwAKCC8UXo4LYBeKBLI4uva5iuSUQkhxKdgAOxt8VAHkCNAnkA2H8lFaO/2osTNzPg5+GMlZMjqwzk9SoG8gDEnvxt55Jq1B5b0I+Z93R1qvRcToVzukqlwhcTOwMozep6Z1Q7dAr1haeLFi5aNY7fzMCwz3fhl0M3xKwGImIwT0RUI6dvZVZ73Yd6lKeVFxSX4NvHuomPrZlAKKKJH5Y/0QN/vXAXRnRoCJUK2HImEd/ujrXiVWxLP16+RUDlXnm9AWENAJT2JAGlF3H6ZQAQ4OVS5X54MUdEUqvN8eiTfziM+Ix8tAzwxG9T+6B7U787er2QeqU91leSlVP1XX8T1svEjd1cE4X6RnduhJNzh+LhyCZw1qrx29TeOPXOMPzzYl90beyL7IJivLL2BJ5ccQiJVsyaQlSXMZgnIrIxQQDuaukPAHikZxP4e7qgY0hpheL7OgVb/XrtG/ngq4e74oW7WwEA8oqU02utn46vlYVgvltTP/z+XB+j6ZRGdGiIyX2bYdlj3ar1O4lTUDYCETmGzacrFzWtiolOdACATgD6tvLH2md7o7GJWU6spR93XlSinBudWRbS7PXt7dLY12i5l0Evvj4DoXkDT/w6pTdeu6cNnDVqbD+fjKGf78Lvx+J4Y5ccHsfMExFV4U4vFnSCgG8e7YaLSVnoUDbN0K9TeiEjtwgBBlMKWcvVqfR+bLGCLt4uJmUBAFpaCOYBGFX+B0oL/80ZEQ6gtOL/ygPXLBZy+s/fZ7HkkYg7aywRURWyC4rx9a4rOB2XgX+tTGFvUt8dCyZ0xgOLS2fxCAv0QkNfV+w4n4yonk3w9n3h0Gpqp1/Nqex1CqsxbapU9Kn07i6V68J8PK4j1h65ifurGGKlp1GrMKV/C9zdJgAv/XIcJ+My8OKqY/jnZALef6A9/D2rzugiqosYzBMRVeFOL450Qul4xo4hvuIyF60GAd7mpx2qDm1ZFfgSnXKC+cvV6JmvSkMfNxyYPRid3t1sdp1/TiVg46kEDG8fVOP9EBFZkpJdgPn/nMOvMTdrtP2OlwcYjW9//u6W6NvKH9dScyvd0LxTztqyYL5YOcF8UVlbXLSVb1jU83DG032bW/2arQO9sG5qbyzefhmLtl3ExtMJOHT1Nv7zQHsMb9/wjttMZG+YZk9EVIX8wqovjj4e1xHjIkLQoZEPejY3Hvs4pmv1eh6spVGX9cwrJJjPLijGrYzScYxV9cxXxcfdCdGz78auWQNx+YN7jZ6bOqAFAOCt308hI69mxaeIiCw5cCUV3d7fWuNA/u42AWIgP2tYGPq1boCh7QLh6+5c64E8ADhpSvdVccYUOenb4lRL2Qd6Tho1XhzcCr891wdhgV5IzSnElB+P4MVVR5GeyynsyLGwZ56IqApVFT36aFxHPNgtFA8azJ9+Nj4TAV4uSMstuuPA1pyyazeUKGTMoL5X3t/TBb7uzlWsXbWGPqanIJo2qBX+OZWA2JQczN94Dh880OGO90VEZGjeP+eqve6yx7rhqe/Lp9fcNWug0Tj45wa2xHMDa7V5leh7v5UUzBeWDQGr7WBer30jH/zxQh98sfUilu68jN+P3UL05VTMH9sRA9sE2GSfRErDnnkiIjMEQcDOC8m4/6u9ZteJbOaHMSbG/LVt6I36ni42C+QBQFN2gVSikDHzF8uC+ZYBHrX+2s+W9cZ/HRUBVycN5o0pDeB/OnAdBwymDCQiqg06K26SDmobiJg3BouPvd2k7ysTx8wrKM2+2EY984ZctBq8MrwN1j7bG80beCApqwBPrDiEV9YcR1Y+M7eo7mMwT0R1yplbmbidc2dpdiU6AX+duIWRi/bgse8OIsFgCpxn+hmP8ftpcs9aK2BkLf2YeaWk2V8Sg/nav4Hx6vA2uPrhCAxtVzpGvmfz+uKUf7PXnUR+LU4ZRURkeE4J9jFdqLSeuxO+jiotxFnf0wUTu4fi4cjGtZKZZK3yAnjK+D4ADNPsrZmEtWa6NK6HDdP64um7mkGlAn45fBPDF+zG3ksp4jpHr6fh+Z+O4EJiFl74+Sg21WB2AiKlYZo9EdUZZ25l4t6Fu+GsUePCf+6xevv8ohKsOxKH/+66jGupuQAAd2eN0Xy44cHe4v+XPhIBjdr2Fynm6PdtTQ+SLZVPS+clyf5eu6ct/j2bhCspOfhy2yW8PCxMkv0SUd1WohNwtew74OfJPbFo20WxHoihA68PFgvPAcCHYztK1saKnB0wzb4iVycN3hgZjqHtgvDyr8dx/XYuHv72AF4a0hovDGolzirw14l4AMCfx2/h6ocjJGkbka2wZ56I6ow9l5IBWF99Piu/CEt3Xkbfj7bj9fUncS01F/XcnTBjcGvsffVufPpgJ3FdH7fyOXCDzPTWSEWjUlrPfPWmpastPm5OeHd0ewDA0p2XcTY+U5L9ElHdduN2LgqLdXDRqtGjmR9cnUzPPOJsokq7XBw1zd6UHs388M+LffFwZGMAwIJ/LyrqJgdRbWLPPBE5rOSsAqzYF4sfoq8hK790PtyGPq6Y3Lc5JvYIhbtz6SlyTNdGyMovQtcm9ZCWWz4Gr02QND3Q5mg1+qnp5L9IyS8qwfXbpT1ZdzItnbWGtw/C8HZB2Hg6Aa+uPYH1U/vImi1BRPZPn2XUvIEnNGqVyanVlMZZo7yeeSnT7CvycNHi3dHt8dPB6yjRCfjnFFPqqW5iME9EdZKlOchv3M7F17uu4JfDN1BQ1ovRMsATU/q3wKhOwZV6W1QqFR7v0wwAkFtYjKb13dGlcT2zvTVS0ShonvnYlBzoBMDLVYsGXi6S7vud0e2w93IKTtzMwPK9sRbnLo5Lz8Nbv53CU3c1Q++W/hK2koiUTBAEXEzKxq4LyfjtWByA8iwjU+d6f0/px8Vbosh55iVOs69Io1bBz90ZqTmFmPbzUZPrpOcWylLjgKi2MJgnojppyo8xeHNkOFbuv4a5o9qhfSMfJGXlY8mOy/jrRLwYAHcO9cXUAS0wuG0g1NXo0XV31mL7ywPE+YPlpE+zV0IwXz5e3lPy302gtytev7ctZq87iU83X8CwdkEI9XM3ue6sX49j3+VU/HsuiWMliRzc7ZxC7LmUgt0XkrH7YopRsVMA6Ft2w++hHo2x/mgcOoX64uWhrfHRxvPijBpKoe/9tnaYmS3Zap55a1SVVRH5wb84/771NXaIlILBPBHVGSoYB5Hv/XUGAPDodwcrrduvdQM8278Fejb3szr4VEIgD5T3zCthzPxFG1ayr44J3ULx29E4HIi9jdfXn8QPT/YweZz0Nx2IyPEUlehw5Foadl9Mwa6LyTgZlwHD+qEuWjUim9dHv1b+6Ne6AVoHlg6l6tHMD/++1B+NfN3g6qRB31YNZHoH5imxAJ6+LVoZ0uz1sgqKLT5fUKxDYbFOUfUPiKzBYJ6I6gwBloNatQq4p0NDPNu/Bdo38pGoVbajv0DSKSCYvyxxJfuK1GoV5o3pgOFf7MbuiylYfzQOY7qGVFpPSb1WRFT7BEEwupF3NSUHuy8mY+eFFERfTkFOofE0lm2CvNC3LHjv3tR8sbsWDeS5UVldzgosgKdPs3eWsWf+vk7B+OnAdYvrJGbmm83mIlI6BvNEZPeyC4rh6aJFVTHtIz2biNXP6wKNuvQCSRk989JWsjeleQNPTB/cCh9tPI93/zqDfq0bwN/TePx+QZFyLnSJqHat2BuLjzedx8rJPXH8RjqW7YkVC3Pq+Xk4466WpcF731b+CPSWd1aS2qJPZdcJpVXktTIG0HpKSLN/+75wtArwxDt/njG7zi+Hb2DaoFaytpOophjME5HdScrKx+vrTiKqV1Psu5yC/+68Uq3t3GQuWFfblDJmvrhEh9iUHADyBvMAMLlvc/x5PB5n4zPx3l9n8MXELkbPs2eeyP6V6AQcvnobHUJ8xFlHikt0mFsWsP2w7yr+OhGPwhIdtGoVIprUQ7/WDdCvVQO0C/auVn0Ue2OYUVBQrKxgXs40exetBk/0aWYxmF+07RIEAXh5WJiELSOqHQzmiciuCIKAHv/5FwCw9WySVdsqZax7bVHKmPlrt3NRVCLAzUmDRr5usrbFSaPG/LEdcP9Xe/H7sVsY3TkYd7cJFJ83vPFx4EoqvFydEB7sLUdTiaiGlu+Nxft/n0X7Rt54/Z626NqkHtYdiROf33ImEYUlOjhr1Ih5czC8XJ1kbK00DAu95RWVwMNF/kt8JaTZVzTn3rYY0bEhen+4zWj5sj2xDObJLinn00VEVA0X76CAWV3rjFHKmHl9UbkWAR6K6PHqGOKLp+4qnUrwjfWnkG2mANKkbw/g3oW7cbUsq4CI7MP/9l8DAJyKy8Skbw/goW/24/X1J8Xn9UXPQv3cHCKQB0rrhrg6lV7W51WoCyAXJaTZ6619them3d0Sj/dpimATN53zikqw71KK7N+nRNaS/9NFRFRNey6mYMuZxBpv7+kqf09FbVKrlNEzrw/mWyqoQNSMIa0R6ueGWxn5+GTTeZPr6HvpB3yyA6dvZUjZPCK6A+oKWVZHr6ebXE/pRetqmz7VvqBYWcG8nGn2ehFN/DBzaJh4Y+Hpshu+hiZ9ewADPtmhqBkBiKpid8H84sWL0axZM7i6uiIiIgK7d++2uP7OnTsREREBV1dXNG/eHEuXLq20ztq1axEeHg4XFxeEh4dj/fr1Vu23qKgIr776Kjp06AAPDw8EBwfj0Ucfxa1bt+78DRMRgNL5gB9ZdgAfmwnMTBkSXp5e3SnUF4/0bGKLpslGq1bGmHlxjvlAeSrZm+LurMW8BzoCAL6PvoqYa2kW13/k2wNSNIuIakF1R0w93rupTduhNPq6MHmFyghG88sKjlY117scXr2nDZ4b2KLS8uu3c3EyrvTm7tYziZj4dTRupuVWWo9IKZT36bJg9erVmD59OubMmYOjR4+ib9++uOeee3D9uukpJ2JjY3Hvvfeib9++OHr0KF5//XVMmzYNa9euFdeJjo7GhAkTEBUVhePHjyMqKgrjx4/HgQPlF3ZV7Tc3NxdHjhzBm2++iSNHjmDdunW4cOECRo0aZdtfCFEdp9MJeOO3k/hx/zXcuF39L9OPx3XEqXeGGaXL/f5cH3jXsXRLjcKCeaX1gt3Vyh/jIkIgCMBra09YnLIpLbcIxSU69Jr3L5q+9jeW7rwsYUuJyJLiEp3RcBiNmWh+bNcQfB0Vge0vD8C+1+5G75b+UjVREfQ98/kK6ZkvKDvnmpvuT05OGjXuad/Q5HNjFu/DV9sv4ekfDmP/lduYve6kyfWIlMCugvnPPvsMTz31FJ5++mm0bdsWCxYsQGhoKJYsWWJy/aVLl6Jx48ZYsGAB2rZti6effhpPPvkkPvnkE3GdBQsWYMiQIZg9ezbatGmD2bNnY9CgQViwYEG19+vj44MtW7Zg/PjxCAsLQ8+ePbFo0SLExMSYvdEAAAUFBcjMzDT6IaJy/55Lwo/7r+ON306JYySrsn5qbzzYLRSeLto6X7lcCQXwdDrBoGdeWcE8UFrsyN/TGReTsrFo20WL677063HEZ+QDAD7855wUzSOiKuQWFmP66mMY8MkO/H6stMhdxTR7vQZeLhjaLgjN/D1Mjouu61zFnnmFBPNFpe1QYjAPAEE+5qclNMwC3H0xBbPXnZCiSURWs5tgvrCwEDExMRg6dKjR8qFDh2Lfvn0mt4mOjq60/rBhw3D48GEUFRVZXEf/mjXZLwBkZGRApVLB19fX7Drz5s2Dj4+P+BMaGmp2XSJHZNgTsybmpsl12jY0rkRu+OXcLtjHNg1TiPI0e/luWtzKyENeUQmcNCo08XOXrR3m1PNwxtv3tQMALN5hubf992McGkUkl+ISHV5ffxLrjpSf65Oy8hH+1ib8dSIeALC0bBpSc4U2R3Y03dPqKNz0BfCKlBHM6zME9IX5lKa+h3O11/354A0btoSo5pT56TIhJSUFJSUlCAwMNFoeGBiIhIQEk9skJCSYXL+4uBgpKSkW19G/Zk32m5+fj9deew2TJk2Ct7f5KY9mz56NjIwM8efGDZ4oiPQ2nIzHfzacrXI95wqFdVy05T0A0wa1xAt3t8RfL9xV6+1TAiWk2etnF2ha30MR8xqbMrJjQwxqEyD+nrRqFT4a21HmVhGRoT2XUvDTgeuY+ctxFJdlVW0oC+L1tGoV9l9JRUZuYaXtG3i5oH2jun0Dtypimr0CgvkSnSBOTeeqVWbPvEqlwvG3yzvrJvetXBTPEAvjkRIp88rLgorzRAuCYHHuaFPrV1xendes7n6LioowceJE6HQ6LF682MI7AVxcXODt7W30Q+TorqXm4MttFzF15RGTzztr1fjfUz3Ex8/f3croecNCO+7OWrw0NKzOXuBp1aXvVc5g/rKCU+z1VCoV3ru/PTycSy8o3Z01GN89FOum9pa5ZUSkZ3geOxufBQCY++cZo3VOxmVg4tf7catsOIyh1go+B0nFTUHBvGEblJpmDwA+bk54Y0Rb9G3lj5lDwvBAl0Zm180tkP/3SlSR3QTz/v7+0Gg0lXrDk5KSKvWa6wUFBZlcX6vVon79+hbX0b+mNfstKirC+PHjERsbiy1btjA4J6qBd/88g082XzD7/JT+LdCxka/4uLGfO76OihAfK7Fqrq2UxfKyjpm/mKi8aelMCfZ1w6v3tAEA1Pd0AVD6t2PJcyuPcMo6IokYpobfSMvFkeuWZ6AAAO+y6UZ7t6iPN0eG26xt9sK57PuvsET+udINg3mlfy8/3bc5/vdUJNycNcgpKDa7Xnah+eeI5KLsT5cBZ2dnREREYMuWLUbLt2zZgt69Tfeu9OrVq9L6mzdvRrdu3eDk5GRxHf1rVne/+kD+4sWL2Lp1q3izgIiss/tSisXn3Z01cDEYf6dSGWfOKDXV2xacy95rUYlOzDqS2qXksmBeQdPSmfNIZBN88EAHzC9Lsfdw1lpc/++T8RixcA8A4GJiFrItXOQR0Z0xDKISM/MxZrH5ukRAacr97lfuxtUPR+CnyT3RJogdKOLQKwWkg+eXVbJ31qjN1jhQolwLxQNz+R1ACmRXV70zZ87Et99+i++++w5nz57FjBkzcP36dUyZMgVA6Rj0Rx99VFx/ypQpuHbtGmbOnImzZ8/iu+++w7Jly/Dyyy+L67z44ovYvHkz5s+fj3PnzmH+/PnYunUrpk+fXu39FhcXY9y4cTh8+DBWrlyJkpISJCQkICEhAYWFlcd1EZF5YSaCwm8e7Sb+391ZY3SX31I12rpOXx9AJ8jTOy8IAi4mlqbDKr1nHigtmjUpsjF6NPMDUP2iTH0+3IYhn+/CG+tPIiEjHyv2xjKwJ6plqw6V1w1KzCyocv3T7w6Dj3vdmm70TmkVMMOJnr5n3kWhxe/M+b9+zQEA9Uz8bf0acxOvrrE8zSmR1Cx3SyjMhAkTkJqainfffRfx8fFo3749NmzYgCZNmgAA4uPjjaaCa9asGTZs2IAZM2bgq6++QnBwMBYuXIixY8eK6/Tu3RurVq3CG2+8gTfffBMtWrTA6tWrERkZWe393rx5E3/88QcAoHPnzkZt3r59OwYMGGCj3whR3aDTCdh8JhFLd17GybjKac0NDQJ2NycNVCoVds0aiMKSEni7OuGulv5o5OuGtg2V3ztcmwwvkgqKdXCSOCshObsAmfnFUKuA5g08JN13bVCpVHB1UiO/yPKFWVx6HgDgt2O3cOR6Oq7fzsX5xCzMG8MiekS1oUQn4Oj1dPHx1ZQcBHi5ICnLdFA/unOwUbFTKqUpG3ulpGBeyePlTenXugF2vzIQqTmFuP+rvUbPfb2rdDaFrk18MaF7YzmaR1SJXQXzADB16lRMnTrV5HMrVqyotKx///44csR0IS29cePGYdy4cTXeb9OmTWVLcSWyZ4XFOvx2LA5Ld17GleQcs+sFeLmI/9ePCWxcv3y8s5uzBjtnDRBTDB2Fs0HwXlBUAk8XaU/pl8rGy4f6udvdBZve4TeGIL+oBBP+G43LyTl4qEdjtGjggff/Nj2TwvXbuQCAbeeSpGwmUZ1UVKLDJ5vPo1OIr9Hy6CupFiuHy1n0U8m0CpjhRE9/k1Sp09JZEurnjgKD3vfmDTyMrlGy8pmZRcphd8E8Edm/nIJi/HzwOr7dHYuEzNKqxN6uWrg7a8XHeu+Nbgc/g7lgC8yktznSWHk9tVoFZ60ahcU6s78XW9KPl28VoPwUe3M8XbTwdNHih6cisfFUAiZ2D63WfPNSZ0EQ1UWrDt3Af8vmjjeUkVdkcTv2n5imKZuqtVgBBfD0N2Oc7fRcGernBo1aBXdnDYJ93IyC+ff/Potvd8diw4t9ja5PiOTAYJ6IJHM7pxDf77uK76OvIj239GItwMsFT/dthkmRTTD5+8NiMD9vTAcMDQ8UK4/rFShgyh0lcZEzmC+blq6FHQfzeo183fDUXaVzDLs5V33xeTMtD8/+GIN3RrdDgJfj1m0guhOnbtZstggl9DwrkZPYMy//mG59MG+vNz5dtBocf3so1Cpg5urjlZ5PyMzHH8fi8Hgfy3PTE9kag3kisrm49Dx8u/sKVh28IU4/1MzfA8/0a44HujYSxz72alEf0VdSUc/dCQ/1MB6P9mBECLafT8LIjsGSt1/JXLQaZKEYBcXS3+TQT0vXKqBu1SpwqzBkoJGvmzhu3tA/pxLQKsATM4eGSdU0ojohLacQ7/11BuuOxhktj2hSDyU6AcdupJvcTl/j4p4OQRK00v4oacy8PjvAXoN5AOLQNXNDBU6YqPFDJDUG80RkMxcTs7B05xX8fixOvLho38gbz/ZvieHtgyqNcZ/ctzk8XLQY0jaw0mt9/GAnlOgEhxsXXxV9Zf+CKoq42YI4LV0d6Jk35GYwZV3rQE/882I/tHh9g8l17WnKJSIl2HAyHm/9fgop2YVQqwAPF604BjmknhvuaR+EKT+W1joyvJHWtL471k/tgzPxmejdgtP/mqLVKGfMvL5nXt8me2bu5sjeSykQBMFoelwiqTGYJ6I7IggCVh26gS6NfcV5fo9eT8OSHZex+UyiuF7vFvXx7IAWuKulv9kvPjdnjZjqbAoD+cr0Fe3zJR5+kJFbhOSyStN1Lpg36Jl/eWiYxb87d2f7LPxHJLWkrHy89dtpbDydAKC01sZH4zoi1M8dgz/biez8Ynw4piPcnDVY/kR3xCbnYFJkYyzbE4uPN53HA11CUM/DGX1a+sv8TpRLf64qUsSYefvvmdczd3MkMbMAl5Nz6tx3INkXBvNEdEd+PHAdb/52CmoVsOPlgXhl7XHsv3JbfH5Yu0BM6d8CXRrXk7GVdZd+iILUY+YvJZfOL9/Qx1XyKvq2ZtiT1Laht8V19RkRG0/Fo2WAJ1rWsSEHRHdKEASsOxKHd/86g4y8ImjVKjw7oAWev7uleP7a+fJAaDUquJXdHBsYFoCBZaNXnhvYEuO7hcLfk4XGqqJV0Jj5Yp1+zLz934Q3dXOkfSNvnIrLxO6LyQzmSVZ16wqMiCT31/HSyt86AXht3Qnsv3IbGrUKY7o0wjP9mzO4sTExzV7iYF4/Xr4uXsS0C/ZGRJN6aBPkhVA/d4vr5hWVIPpyqpgWfPXDEVI0kcguxKXn4fV1J7HzQjKA0gDoo7GdEB5sfJPMx93J4us08HKx+DyV0vfMK2HMfF3qmTc1TeLIjsE4FZeJ7/bG4u42AWhS30OGlhExmCeiGigq0WHl/mto6u8hFrQDgH2XUwEAa6b0Yk+8RMqDeWnT7PWV7OtiMO+i1WDts72Nlo3p2gjrjsRVWnfxjstibxgRldLpBKw8eB0fbjiLnMISOGvVmD64Ff6vb3OHnEZUKvrAWVFj5tX2f7xD/dwqLRsXEYIf91/Djdt5GLtkH757vDs6hvhK3zhyeAzmichqm08nYu6fZwCUFigy5OqkRudQXxla5ZhcysZ3S10A72JS3axkb84HD3RAbkGJON7X0MJtl8T/5xQUw6OODTsgMlRVwa+rKTl4de0JHIgtHW4V0aQe5o/tWCdv/CmNknrmi0vqTpr9S0PCkJVfjN+P3RKX+Xu6YN2zvfH48kM4E5+JiV/vx5JHItC/dQMZW0qOyP5vlxGR5OIzyqfpqjhlV7CvGyu7SshVpjT7utwzb4qrkwZLoyLwx/N94ONmPiU4qawoIFFddD01F73mbcPSnZcrPVeiE/DNrisY/sUuHIi9DTcnDd6+Lxy/PNPLYc4TcisfMy9/MF9Yh9Ls63k444uJXTAgzDhQD/B2xepneuKulv7ILSzBUysOYW3MTZlaSY7K/j9hRGRzM1Yfw2PfHYSu7AIhMTPf7LperpbHPlLtEnvmJUyzzykoFm/itHKwi/SOIb5Y8UR3s89vPFW5556orvji34tIyMzHh/+cM1p+ITELY5bsw382nEV+kQ59WtbH5hn98ESfZpyFRELl1ewVUACvDk1Np/fe6PYYGNYAP02OFJd5uTrhu8e7Y3TnYBTrBLz063F8tf0SBEH+GyrkGJgLSERmFZfocOhqGtYfLR0rfCUlG0E+bvhmd6zZbbowxV5SchTAu5KcAwCo7+GMeh6OV2G6S+N6ePquZvh2T+XPwedbLuDuNgEIC3KM4QfkWAw7WXecT8LKA9fRyNcNKw9cQ1GJAC9XLd4Y0Rbju4UyQ0sGSuqZ16f6O9WBMfN6oX7uWP5Ej0rLnbVqfD6+M4K8XfHfXVfw8abzSMjIx9xR7Xgzi2yOwTwRmfXRpvP4etcV8XFmfjFupadZ3Gbm0Na2bhYZ0AfzUs4zfzGpdFq6Fg7WK2+o0ETPV+tAT1xIzMaM1cfw23N94KytOxex5Nh0OgFqtcqoHsTjyw8ZrTO4bQDev78DgnxcpW4eldGUBc5KGDNfWHaD2UnrGMGsWq3C7HvbItDbFe/9fQb/238NyVkFWDCxM1zLMuiIbIFXGkRkRBAE/HzwOk7cTDcK5AFgzOJ9uJycbXF7b6bZS0qOeeYvicXvHDiYN/h9DwxrgHdGtcOPT0einrsTzsRnYtG2izK2jqj2vP/XGXT7z1YkZubD00xxx4FhDfDNo90YyMtMn9KujJ75ulPN3hpP3tUMix7qAmeNGhtPJyBq2QFk5BbJ3SyqwxzrE0ZEVdpxIRmz153EqC/3mnz+nbIq9oae7NMMTeq7Y7mFscRkGy5OZWn2Elazd7Tid6a0bVg+T/byJ3rgsd5NEeDlivfv7wCgdMq6o9ctZ7EQ2YNv98Tidk4hIj/412yQ2L6RD9PqFUCrqGr2pW1wxAylkR2D8f2TPeDlqsWhq2kYt3RfpWLBRLXF8T5hRGTR5STLPe+mvHVfOHbOGoiBYQE2aBFZIsc885ccbFo6UyZFNsbr97bB39PuMlo+omNDjOoUjBKdgJd+OY68QumOC5GtLd5RuYo9ANRzd7zaGUqkEcfMy18Ar1CcZ94xb/L0alEfv07phSBvV1xMysbYxftwLiFT7mZRHcRgnshB6XQCYq7dRk5BMQRBwJIdlzHpm/1YfeiGyfXv6xRs9NirLN2yDQt9yUrqNPuC4hJcu50LwLF75p00avxfvxZoF+xT6bl3R7dDoLcLrqTkYP7Gcya2JlIOnU7A6+tPVppuLrewWLxxVxVL0zWSdPQp7UUlyumZ19aBqelqqk2QN9ZN7Y1WAZ5IyMzHg0uiEX05Ve5mUR3DAnhEDurPE7fw4qpjuKulP+p5OOPP47csrj+pR2P0a+WPWWtOAADef6A9Cop0uLdjQymaS2ZIXc3+akouSnQCvFy0CPR2kWSf9sbX3Rnzx3bE48sPYcW+qxgaHojeLf3lbhYRikt0ePS7g2gV4Il3RrfH/6Kv4s3fT4vPDwhrgDZBpUNI3vnjDFYfNn1ztyIlTIVGhj3z8gfz+r8J5zo0NV1NBPu64dcpvTD5h8M4dDUNj313EJ9N6ISRHYOr3pioGhz3dhmRg/v54HUAwJ5LKVUG8gDQwMsZA9sEwNVJDSeNCj2b18f47qFmCyKRNPRVcgskqmav76lrEeDJMbIWDAgLwKTIxgCAWWtOIDOfBZBIfgev3sa+y6n4PvoaABgF8gAwfMFu7LqQjMe+O1hlID99cCvoTwFdGtezSXvJOk4a5cwzr0+zd3Lgnnk9X3dn/O+pSAxvF4TCEh1e+PkovjMxtSlRTfAqnMhBBXpbV3W4gacrfNydsPuVu3E7p9Dq7ck2pO6Z109L58iV7Ktrzr1tsediCq7fzsW7f57BJw92krtJ5OCKDdKvzd1gevS7g1W+zgNdGuGFu1thYvfGiM/IQxiHWymCW9nNXSmnKjVHX5SV07KVcnXS4KuHu+KdP0/jh+hrePevM7iVnocm9d0xICwAoX7ucjeR7BRvlxE5KI2VRWm83Urv/TXwcuGFm4Loq9lLdfHGSvbV5+GixafjO0GlAtbE3MSWM4lyN4kcnNogm6bLu1tq9Bpv3xeOD8d2gEatQpCPK3vlFcTNuTRwzlVA4U39d5KLA1azN0ejVuGdUe3wyvAwAKUzRbz5+2mM/2+0zC0je8ZPGJGDSssprNZ6DX1c8VCPxkypViipC+CJlewDGcxXR/emfvi/vs0BALPXnUBqdoHMLSJHZngar+m46if6NBPPO6Qs+mBeCbNo6L+T2DNvTKVSYeqAlmjgVV5zJj4jX8YWkb1jME/kgI7fSMf288nVWnf3KwMxb0wHG7eIakrKNPviEh2upOQAAFo2YHZGdc0Y0hqtAz2Rkl2IOetPQRDkL05FjudUXAbe/fOM3M0gG3J3Ks2gy1NAmr2+Z97ViaGGKU4VsiNZV4VqimPmiRzQpG/2W3zeSaPClP4tMCAswKGnlbEHUs4zfyMtD4XFOrg6qdGonpvN91dXuDpp8Nn4zrj/q73YeDoBvx2LwwNdQuRuFjmQguISjFy0x+rtVv9fT1xLzcV9nYKxJuYGU+oVztW59Psgr6gEgiDImlGXX6xPs2fPvCm3KvTGz1x9DF9HdcO127loWt+d2ZBUbbxKJ3JAOVWk4J16ZxheGhqGiCa8cFM6F7Gave175vUp9s39Pa2uueDo2jfywYuDWgEA3vr9NOIz8oyeLyguwZYziUjMZLol1a6U7AKMWGh9IA8A3Zr6YXz3ULg5axDVqynaN/Kp5dZRbXJ3Lu2jEwQgX4LvBEv030ku7Jmvlq1nkzDnt1MY+MkOzP3jdNUbEJXhJ4zIwehMjJOseAOYd9Lth5Rp9hwvf2eeHdACnUJ9kZVfjFfWnDBKt5+z/hQm/3AYL/1yXMYWUl2TkVeER5cdFD+71nh+YEvetLMzbgbj0+VOtdf3zHPMfPXppwzWTx1JVB0M5onsTEJGPk7ezKjRtgdjb+P+xXsrLfd2dcLyJ7ojyNsVK57ofqdNJAlJmWavn5auZQMG8zWh1ajx2fhOcNGqsftiCn48cF18btOpBADAnkspcjWP6pi8whI8/f0hnInPNLvOwdcHYfOMftjx8gBx2QcPdEDsvHvx8rAwCVpJtUmjVsG57Dsht7BY1rboMwNYzZ7ItvgJI7IT284l4lJSNnrO+xf3fbnHqp6Wa6k5ePbHGIz/bzROmLgRoFGrMDAsAPtfH4QBYQG12WyyMX2vhxQ985fZM3/HWjTwxGv3tAEAfPD3WVwtKyhYwqJ4VIsKi3WY8mMMDl1Ng5erFi0aeFRap7m/BwK8XdE60AtN/T3QJsgLzho17uvUkON17Zi7QiraFxSxZ/5OFJfIO0yC7AcL4BHZgWM30vHkisNGy47fSK9yru+M3CIs2nYR30dfRVGJALUKeKhHY8wY0hrnE7Lw8LcHAADOLHJnt/S9HoXFOuh0AtQ2SosVBIFzzNeSx3o1xebTiYi+koqXfj2OX57pBR2DeaolJToBM1Yfw84LyXB1UmP5493xffQ1XE7OMVqvQ4jx+PffnuuDohIdvFydpGwu1TJXrQZAkexj5vOL2TN/J/KKSlCQVwQntRo+7vxMknn8hBHZgXMm0iQtjWUsKtFhxd5Y9P9kO77dE4uiEgH9WzfAxun98J8HOsDf0wV9Wvrjo3Ed4e2qxaJJXWzZfLIhF4Nej0Ib3smPz8hHTmEJtGoVmtSv3MtH1adWq/DJ+E7wctEi5loavt51BYax/KqD17HxVII4tRNRdQmCgNfXncTfJ+PhpFHhv1Hd0K2pHzxdKveOzqqQRu/qpGEgXwfop4LLl2DolTk6nYBCzjNv0YgODY0ePxhhPMPJzF+Oo9v7W9H1/S0max0R6TGYJ1KotJxCZOSVzjtq6svwrxO3Ki0TBAFbzyRi2Oe7MPfPM0jPLULrQE+seKI7vn+yB1oHGs8NPr5bKI6/PRTdm/rZ5k2QzRn2etiyov3Fsl75pv4ecGImxx1r5OuGt+4LBwB8tuW80TCJ19adxJQfY/DxpvNyNY/skCAI+GDDWaw+fANqFfDFxC7o37oBAMDDuXIiZkg9d6mbSBLQXy/IeTMwx2C8vqcLk4BNmT+uI94d3U58PKarcTC/5UwigNJMm9ScQiRk5GPZnlhkcT56qoBXZEQKVFBcgi7vbUGndzbj0NXblarNA6XTmBhWwz59KwMPf3sAT/9wGFdScuDv6YwPHuiADdP6WhwHz7GR9k2rVkGfpGHLInhiij2L39WacREhGNw2EEUlpntdVh+6IXGLyJ59tf0SvtkdCwD4cExH3GvQ8+dRIaC6uw1ro9RVLmIwL1+afXZBaTCvVauYZm+Gp4sW93dpJD5u4OVidt34jDy88PMRvPfXGcxed1KK5pEd4e0yIgUyLG734NJos+tl5BWhoFiHTzefx68xNyEIgLNWjafuaoapA1owZdIBqFQquGg1pePrbFgE71JZJXsWv6s9KpUK88Z0wJEFabidU1jp+fqezjK0iuzR9/uu4pPNFwAAb4xoi/HdQ42e93Itv9z7v37N8Wz/FpK2j6TjWhY8y9kzn51fGsx7umrZYWCBYb0iVyc1mtR3x7XU3ErrnU/IwqGraQCAv07E48tJkjWR7ACDeSKFSczMx4iFe6q17vyN5/D7sVvILatae1+nYLwyLAyhfkyfdCQuTuqyYF6CnnkWv6tVDbxc8J/72+PZlUcqPVffg8E8VW3dkZt4+4/TAIBpg1rh6b7NK63j71ne6/d476aox7+tOkvKGU7M0ffMM8XeMlcnDZ7o0xS5BSUIqeeOzTP6IeyNjZXW+/NEvAytI3vBTxmRwmw6nVDtdX8+WJqG27WxL94YGY6ujevZqlmkYC5iT4xtLt4EQRDHzDOYr333dGiITx7shJd/PW60/Mj1dAiCwJ4tByEIAj7fehHhDb0wvH15inzMtTQEeLmYvEm7+XQCZq05AaA0SJ8xuJXJ1/YzCN6Z9ly3iQXw5OyZZzBfbW/fVz5u3kVruljgrgvJRo8FQUDMtTRcSc7BmK6NoGUdG4fGTxmRQqRmF+DP47cw988z1d4mpJ4bXrunDUZ04LzAjszWPTGpOYVIzy2CSlU6TzrVvnERIcgrKsGbv50yWr7tXBJ+iL6GYF83zBvTQabWkRT2XU7Fwn8vAgD6t26A8GBvLNlxWXz+6ocjjNbfeykFz/90FCU6AWO7huCtkeFmvwcM0+ydGczXaUoogKdPszf8u6PqWTOlF8ZZGF4JAPd/tRfHb2YAABp4u2CghbpIVPfxU0Ykg5hrtxGXno9RnYLFZQM+2YGs/GILW1W2dWZ/TvtCYk9bgY0u3i4mlvbKh9Zz59+bDUX1bFIpmH/q+8Pi/98d3Y4zCdRhqQZ1E3ZeSMbOCr1xKdkFYrr80etpmPzDYRSW6DCsXSDmj+0AtYXpSg3T7BnM122uWvnT7LPKeuYrFl6kqnVr6oc597ZFVkExHowIQd+PtldaRx/IA0ByZoGUzSMF4hmdSALXU3Px+7E46HQCCopLMHZJNKb9fBTnE7KQml2AA1dSrQ7kAc7fSqVcbHzxdimZKfZK0P7tTYi+nCp3M6iaTsVl4MGl+xBz7Xa11nfWWM6uuvuTHcgpKMa5hEw8vvwQcgtLcFdLfyx8qEuVabahfu6YdndLzBoWZjaVl+oGJaTZ6+eY55COmpncrzlmDmkNb7eqixjrb5yQ4+KnjEgC/T7ejhdXHcNr604YFTe5cTsXQz/fhQlf7ze53YCwBkaP3Z01RtVPiQCDnnkbFcC7XDZevhWDeZubWzb3vCkFxTo8tvyghK2hO/HQN/tx6GqayRlJNp9OQJ8Pt2HbuURxWV4VwVdmfjF+PngdUcsOIiOvCF0a++K/URHVDs5nDg3DcwNbWvcmyO64KCDNvqikNJhnJtGd8apGZgPnnSd+yohs7FRceTrUL4dvGj2XnldklFpZ0UdjOxo9/nJSF6x4ojsCvFzwdVRE7TaU7JaLkz6Yt03P/MWyaelaMJi3ucf7NMOl/9yD/q0bmHy+sFiHXw/fkPVCnapHn22lEyo/98W/FxGXnocnV5QPo8iuRnbW+3+fRXJWAdoEeWHF4z2YxkyVuNq4IGp1FJeU/tEzmL8zlobO6GXlF2PZnlh8tycW6bnmryep7uK3AJGNjVxkfpq5JTsumX3Ox80JAd6uWP5Ed/i5OyMsyEtMqz84Z3Ctt5Psl5hmb6OLt0vsmZeUVqPG90/2wNw/TmPFvquVnp+15gRO38rEmyPDMfOXY+jQyMfkdGSkTOm5heJnCijN0PJ00WLDyerNZNK0vjt+eKoHfNyrTsElx6OEnvlCsWeehXlt7frtXCzbEwsAePevMzg0ZzAaeLlUsRXVJbxlRmRDgmCiS8bA5eQcs8+VlHXnDAwLQKdQX46PJ7NsmWafmV+ExLICO+yZl5bGQq/Min1Xse7ITfx+7Bbe//ssikvk64Wjqq0/ehMnb2bg92Nx6PzuFqMsmrFL9iHquwOIvlJaD2FEh4bmXgYNvFzw49ORCPBytXmbyT6J1exlLICn75nnlGm1y7BuzWv3tAEAbDmTaLTOjvNJkraJ5MeeeSIbupPCJFG9mtRiS6guKw/ma//iTd+DGOTtCm9X9gRKyTCU9/d0QUq2cdXiiwa9uy3n/IOHeoRi3pjyoTnbzyfhyLU0zBjculrpmmQ7M1YfBwAMblt5CqmkrAIkZZUf29aBpVlYa4+UDsuaMbg1inU6DG4biCb13eHr7lzpNYj0lFAAr1hX1jPP806tMjymwb5uJtf5+eB1DAgLYO+8A+EtMyIbysg1X5jE39P8Bdl797fHjMGtbdEkqoNsWc3+UiIr2cvFcMrwrTP74ffn+iDUr/wC7utdV4zW//ngDaNsoCeWH8KibZfw54lbNm8rldt0OgFv/37K5HPxGflVbv9U32aYFNlYfPzcwBZ4aWgYOoX6MpCnKumnplNGmj3DjDulvznj5aLF2K4hAIBuTeqhWX0Pk+sfuZ6O2etOSNY+kh975olsKN1CMP/y0DC8tu4kAOB/T/WAVq3GQ9+UVrUf27UR5wKmatN/2ecV1v7FG6elk0+HEF/x/77uzvB1d8buV+7G2CX7EHMtzeQ2yVkFCPA2TsGOTTE/nIdq3zP/izH73PmELIvbDmoTAE8XLdoFe8PTRYsAbxemKpNV9Gn2cs4zzzT72vPT5J54/68zmDMiHO2CvdGkvjv6t24g1kYwZevZJKTnFvLmn4NgME9khiAIEITqVRM1Jz3PfGXRZv7ld1UbeLmIX34A4O7MjyZVn34u2kwbTFFzMbE0+GAwL72RHRqWTkEW6mu03M3CRdyfJ+Lx1F3NjJbJWdXa0aRZmJ0EAIpNlbYvs/zx7ujZvD6A0oDs0JzB0LKAGFlJf3O3QAFT0znz7/eOdW1cD+um9hEfjynrnQeAp+9qhj9P3EKflv5YdyTOaLt+H23HZ+M7Y3B4oGRtJXkwYiAykJSVj1NxGejXqgHGLtkHqFRY/2zvGgX02QXF+CH6GgCgXbA3Tt/KNHo+wNsV93cORlpuEVoFeEGjVmHWsDCE+rnXynshx+FTFsynWcgEqSl9zzwr2UtPrVYhqmfl2hmWimG+99cZJGTk4Yk+5QG9LQojkmmHrt6u8bYD2xiPp3dzZtFTsp5YAE/Gm3hF7JmXxBsjw/HGyHAAwJC2gXh25RHxucz8Yjz9w2E8O6AFXhrSmseiDmMwT2RgxMI9SM4qwMwhrXH8Zun88Fn5xdWeAkgQBJy4mYGfD17HH8dvIbcs7TnY1w0P9WiMkzczsPrwDQCAh4sGCyZ2Mdr+uYEta/HdkKOoV5ZKV9tzzOYVluBmWh4A9swrScUg78MxHXAiLgM/HbgOAPhmd6w4VRHAnnkp5BeVYNmeWCzZcdmq7YJ9XHErIx87Xh5gm4aRwxEL4Ml4E6+IY+Yl5+Jk/LtWqQBBAJbsuIyj19Ow8KEunAWjjqpRMJ+Xlwc3N9NVFOPj49GwoflpVYiULLmsorC+ijCg79WyHMxn5BXh92Nx+PngDZyNz6z0vCAAj/RsgktJ2WIw7+XCyuBUO+p56HvmazeYv5ycDUEA6rk7ob4nK+MqhVuFi7ZBbQNRz8NZDOYBwDCbW850WzntuZiClgGeCPKx3QWsTifg9+Nx+HjjedyqRnE7Q+MiQvDWfeFQAfDiTBFUS1wUUACvmPPMSy6yWX2jx5+M6wRXJw1eWXMc+6/cxoiFe/DlQ10Q2by+mVcge1WjW2ZdunTBkSNHKi1fs2YNOnbsaGILIvkdv5GO6MupRssEQcDXuy7j+Z+OGH3xXUvNFf+fZ+YLURAEHL56Gy/9chyRH2zFW7+fxtn4TDhr1bi/czBW/V9PcV39NC2NDKYScXXiHWuqHfoiN2k5tZtmr5+WrlWAV62+Lt0ZwzHz3zzaDQ28XODrZj4YPBmXIUWzFGXnhWQ8suwA+n60rdJz+6+k4utdlzHo0x04fAdp8fuvpGL0V3sxY/Vx3MrIR7CPKxZM6IxPHuxU5bbbXx6AeWM6wNvViYE81aryqenkT7Nnz7x0PFy0+Pel/uJjbzcnjOjYEH+8cBdaB3oiOasAk749gKU7LxvNemJOSnaBmGFBylajnvkhQ4agd+/emDt3Ll599VXk5OTg+eefx6+//ooPP/ywtttIdMcEQcDor/YCAA6/MRj+Zb2My/bE4oMN5wCU9pKYUjGYT8spxLqjcVh18LrRPM+tAz3xUI/GeKBLo0oVRPXF7dycNTj8xmA4qdVQqXjHmmpHQNl8sslZBdDphFqbU1wfzLdgir1idW3sC8Byz+7FpGyk5RSinkfdq2ycX1SC4zfSodWoEdGknrh814VkAOVBhV5qdgEmfr1ffPz48kM49c4wq/Z5JTkb8/45hy1nEgEAni5aPDugBZ66qxlcnTTivs3x83A2KoBKVJuU0DOvDwJZwFFafgbXnh4upX8HLRp44rfn+uCN9aew7mgcPvznHGKupeGTBzuJ9XYqupycjUGf7kSPZn745ZlekrSdaq5GwfyiRYswYsQIPPHEE/j7779x69YteHt749ChQwgPD6/tNhJZbeeFZOQVlmB4+yAAxlO0HLuejtScAoyLCMWfx8vnX358+SGTr5VXWAJBEBB9JRWrDt7AxlMJ4hyqbk4ajOzYEBN7NEbXxr5mA/SGBmme/kxXploW5O0KjVqFwhIdkrMLEFg2NVlGbhFcnNQWC6aZEn05FbPXnUBKdmnaPovfKYvhBZhfWYDu5Wr563z90Tg80rPJHU95efxGOgQAnStU2JdDanYBIt7fCgDwcNbg2NtD4aRRIzEz36hmwOXkbGTlF+P7fVdx/Ea60WtkFxRX+wbY7ZxCLPz3In7cfw3FOgEatQoP9QjF9MGtjc7rps7xw9sF4fm7W+LojXSM6MChiGQ7+poaBcU6CIIgS8cBx8zLw7CeiuHv3t1Zi0/Hd0K3pn6Y+8dpbDmTiPsW7cHih7uifSOfSq/z6+HSoaYHY81nLmXkFsHVWQ0XrQaCIGDr2SS0b+SNhj6mh2GT7dS4AN7QoUMxZswYLFmyBFqtFn/++ScDeVKEEp2Ax747CAA4+PogBHi7Gs2//fQPhwEAt9LzcepW5fHtFS3YehHXb+cazdXcLtgbE3s0xujOwfC20CO24onu+OnAdbwyvE1N3w5RlbQaNYK8XRGXnof9V1IxrF0Q0nOLMOjTHegU6oufJves+kUMTPp2Pwyz8Fj8Tlke6BqCcwlZeKRnE/FCPdTPHbPvaYNTtzKNblLqvfvXGXy1/RIe6NIIE7qHolWg9UMn0nIKxQynC+/fA7UKOJ+YhbZB3rWWDWKN346Vv8+cwhLkFpTAx12NWWtOGK036NOdVbxOnNF0TxXlF5Xg+31X8eX2S8jKLy59zTYBeO2eNiZ/jwHelYP51oGeaN/Ix+SFM1FtMrx5W1Css3gzd92Rm1h96AYWPdQFAd61V1tCPwUjx8xLy8XgZm39CplYKpUKkyIbo0MjHzy7MgbXb+dizJJ9eG90O4zvFmp006ew2Di9Pi49D59sOg+NWoWPx3VERl4ROr+7BfXcnfDq8DZw1qox85fj8HV3wrG3htr2TZpRVKLDi6uOom2QN14Y1MrsOmqVChoZvq9sqUbB/OXLlzFp0iQkJCRg06ZN2LlzJ0aPHo1p06bhP//5D5ycOP6L5JOZVz5uOCmrAAHersgpLK603hf/XqzW6+0sS5n0cNZgVOdGmNSjMTqEVO+CbEBYAAaEBVS9ItEdCvR2QVx6Hl5cdQwNfVwRX1aMa9/l1Gr1zmTmF2FdzE0cu5GOisPpWgUymFeSZv4e+PrRbpWWP9O/BRIy8isF888OaIG1MTeRlFWAb/fE4ts9sejS2BcTuoViZKdgeLpU71LAcNq1nIJifLrlPH7cfx2v3dMGU/q3uLM3VQMVx3PmFZXAB044bWWNgKPX08VgPjW7ANdu56Jr43oQBAF/nYjH/I3nxFkdwht6Y86ItujT0t/s6/l7uuDhyMb46eB18bPkbaGmAVFtcjUI6PKLSswG84IgYOYvxwEAqw7dwDQzAZA1MvOL4O3qJAaD7JmXlkqlwhcTOyM1uxDNG5j+3u4Q4oO/XrgLL/1yHP+eS8Kra0/i0NU0vDe6vdizX1hS3gF2KSkLgz/bJT5+fmBLDP289HFabhFeW3dSfC49t0i2bJA/j9/ChpMJ2HAyAWq1Cl0a+6J3i/LzdFGJDkM/3wUPFw3+fP6uOjXUtUbBfOfOnTFixAhs2rQJvr6+GDJkCO699148+uij2LJlC44ePVrb7SSqtnSDYD4zrwh7LqZg69nEO3rND8d0wH2dguFRzYteIqkZ1mmIr1BVOzPP8vSKBcUlGPb5rkrbAaXjgYNqsceGbCvIxxX/vNgX93yxW1z26vA2eGlIa+y8kIzVh25g27kkHL2ejqPX0/HuX2cwokNDjO8eim5N6pm9wMnML8KH/5wTH3d5b4v4/0X/XpQ8mF/070V8uuWC0bKMvCIE+biaLVpqjk4QkJFXBLUKeGz5QZyKy8Qrw8Ow+XQijpWl5Qd6u+DloWEY0zWkWr06/3mgAx7t1RTDFpRe9DKYJ6loNWpo1SoU6wSjIni/HY1DgJcLejavj2M3042G69TG2PZ/zybiqe8PY9awsPIx82oG81Ib3blRlev4ujvjm0e7YcnOy/h083msibmJU3EZWPJIBJr5e6DA4O/GMJAHgCk/xohDTU35cf81RPVqWuP219SpuPJM2483nQcATLu7JWYODcNnWy4g5tptMcM2I6+oUm0re1ajT9nixYuxatUq+Pr6ist69+6No0ePomvXrrXVNrP7btasGVxdXREREYHdu3dbXH/nzp2IiIiAq6srmjdvjqVLl1ZaZ+3atQgPD4eLiwvCw8Oxfv16q/crCALmzp2L4OBguLm5YcCAATh9+vSdvVmqtouJWfh+31Ucvnobj3x7QFz+3E9H8MiyA1ix76rF7Qe3DbT4/MQejRnIk6KZK2QDAMnZBSaXH7+Rjt0Xk7H7QorJQB4oLX5Xl+5gO4K2Db3x7uh2AIAJ3UIBlF7gD2obiK8f7Ybo2YMw+542aN7AA7mFJfg15iYeXBqNQZ/txNKdl5GUVflvYeqPR3DFYKiRoQZe0tcBqRjIA8CwBbswfMEu5BZaF8zHXEvDXR9uQ4e5m8ULwo82nsexG+lwd9Zg5pDW2P7yADzYLdSq9EzDdZv4uVvVJqI7oe+Nzy8qwc4Lyfhf9FVMX30Mk74tvR4as3if0fCTjzaex47zSXe0zzd/OwWgNJA6cj0dAOBdRS0Pko9arcJzA1vix6cj4e/pjHMJWbhv0R58tuUCfo25aXa7cwlZFl/3zd/liX0SMyt/by3cdgk5BcVY+O9F7L1UPptVUpbpayJ7VaNgPioqyuRyLy8vLFu27I4aZMnq1asxffp0zJkzB0ePHkXfvn1xzz334Pr16ybXj42Nxb333ou+ffvi6NGjeP311zFt2jSsXbtWXCc6OhoTJkxAVFQUjh8/jqioKIwfPx4HDpQHhNXZ70cffYTPPvsMX375JQ4dOoSgoCAMGTIEWVmW/+jtiU4nID23EKlmAgM5Dfl8F97+4zTGLY1GXHqeuDwtt3pTdX37WDccnDMIX0zsXOm57k3rVd6ASGEsBfMpJj6zm04nYPRXexG17KDFm10tzaTqkbI92qspTr0zDB+O7VDpuQZeLnimfwv8O7M/1kzphfHdQuDurMGV5Bx8+M859Jq3DZN/OIytZxJRXKLDyZsZ2HMpxey+lFTUs6oLTXPbZBUYD8VSqYCJ3UOx4+UBmDaoFdydrQ9KQuq5IcDLBR0a+aB7Uz+rtyeqKf30dAmZ+Xjsu4NGAda7f50xuc0LP99ZVm19E+cBZqQoX+8W/vh7Wl/0aOqH7LLA9054OFtXcLe2ZOabvt6/kZZbadlHG89Bp6t6ej57oRKqM9lgBT/88IP5F1SpzAb7dyoyMhJdu3bFkiVLxGVt27bF/fffj3nz5lVa/9VXX8Uff/yBs2fPisumTJmC48ePIzo6GgAwYcIEZGZm4p9//hHXGT58OOrVq4eff/65WvsVBAHBwcGYPn06Xn31VQBAQUEBAgMDMX/+fDzzzDPVen+ZmZnw8fFBRkYGvL29rfjNSGP53li88+cZ+Lg5lU6vpqCxUE1f+/uOtr/64Qjx/9dSc9DQxw2Hr93Gst2xeO/+9gj2ZXVOUrYFWy9gwVbTX8JvjgzHp5vPY1hZRe1hn+8SCxRVRa7x0CSt7IJi/H3iFlYfuiH2qgGlgX9yFb0Yg9sG4NvHutu4hcbu9JxflU3T+yEsyPoigRUVleigEwRxujAiKfT5cJtRx0Z1GV4LWeupFYfw7znj3v3drwxEKLNS7EJRiQ6fbDqP/+66ckev4+Wqxcm51k33eScEQcCUH2Ow6bTp4bRfR0Xg//4XU2n5/LEdMKF7Y1s3745UNy6sUf7Liy++aPS4qKgIubm5cHZ2hru7u02C+cLCQsTExOC1114zWj506FDs27fP5DbR0dEYOtS4quKwYcOwbNkyFBUVwcnJCdHR0ZgxY0aldRYsWFDt/cbGxiIhIcFoXy4uLujfvz/27dtnNpgvKChAQUH5RVJmZtWV1eWknwIpI68IS3ZcrpViKbWhtudSbVK/dP7f3i38jYpnECmZr4UekPfKemLWH43D+qNxVr0ue+Ydg6eLFhO6N8aE7o1xMTELvxy+gXVH4qoM5AFYPfWhEtzdJgDbzplOK1arUCuBPMACYCSPmgTyALDjfBL6tmqAEp0AnSCgRCegRBCg0xn+HyaWCSaHc1mq1ULK4qRRY/a9bdGrRX2zUzVXi8Qd3ldScswG8gBMBvIA8Orak2ha3wORzevbqmmSqVEwn5aWVmnZxYsX8eyzz2LWrFl33ChTUlJSUFJSgsBA47HNgYGBSEhIMLlNQkKCyfWLi4uRkpKChg0bml1H/5rV2a/+X1PrXLt2zex7mjdvHt555x2zzytNPYNiEWuP3FRMMG/qYjOknptYfbiiJ/o0xfK9V23cKiJp1ea0QoZYyd7xtAr0wpwR4Zg1rA3e/O0UVh++YXF9tZ3VVPBw1qBfK3+jYH7DtL64d2FpLZw6lH1JZJU7CuJM8GKtIbszICwAVz8cgeELduFcQhZeGtLaZI0Sc6Q+fTpZUWTR1UltVBTyo03nsfbZ3rZolqRq7ZZxq1at8OGHH1bqta9tFQsxVTUFgqn1Ky6vzmvW1jqGZs+ejYyMDPHnxg3LF0xy8zOYs7Li/JVySs0prLTM3JzBG6b1xdv3tbN1k4gkZ6uhICH1mCLpqJy1agxrb7k4KAAU68xXNlaaSZGNcWDOYNSr8B1mam54IjKmUgFatQrOWjXcnDTwdNHC21ULLzOF7lg81X59Or4TvpjYGc8NbGn2+Jqis3709h1x0lb/b2zRQ12x/InyIWF9Wth/rzxQw555czQaDW7dulX1ijXg7+8PjUZTqRc+KSmpUo+4XlBQkMn1tVot6tevb3Ed/WtWZ79BQUEASnvoGzZsWK22AaWp+C4u9nMBYXjxY6rQiVxMpdkPaRtYqYjHR2M7Ijy4dMzJk32a4bu9sQBY4I7qhmCfyj3zzlq1ON9vdQwMa4Dt55PFx88PbGlV9W6qexpXY7xrYbH0XdmdQ31x7EY63hoZjo82nTPqbbFEo1LB00WLitebhjeoLRWTJLIHSx7uimdXHqm0vEtjXxw1qIkBAIsf7oqpZet2b1oPix+OgLNWDY1aBY1KBbW69HOjUavMBucJGfnoOe9f8fFjvZrIMj0Z1Z52wT5oF+wDoDQzNyu/2Oy6kyIbY1KPxhi5aE+lc6utVXd/Izo2xJDwQJyKyxCX3V3FTFb2okY983/88YfRz++//46lS5ciKioKffr0qe02AgCcnZ0RERGBLVu2GC3fsmULevc2nSLRq1evSutv3rwZ3bp1g5OTk8V19K9Znf02a9YMQUFBRusUFhZi586dZttmj/wM0uyVlDpl6oPcIcQHHzxQXsX5p8mRGN89VHz81n3hiJ13L/564S788GSkFM0ksqmKFcUbeLmgVYB1KfKBBqn6nzzYCS8PC6uVtpH9ahnghY/GdsSzA1rgm0e7Yc69baFRq4y+A4oszDlsK/rTfqifO/56oa/Z9fq1bmD0WB+LGPYeLX2kK1QqFVb9X0+0C/bGd49LW8yPqLbd06Ehlj4SUWn5sHZBlZYNahuAXs3rY9qgVvh1Sm808HKBj5sTPF20cHPWwEWrgVajttjLbpi5+fG4jnhndHu0tPL7h5TrvfvbAwCGhpcHv78/1wcn5w5F7Lx78cEDHcSboILEifYV99bLzBh4z7IZSdoEeWFw20A83rspOof62rZxEqlRRHb//fcbPVapVGjQoAHuvvtufPrpp7XRLpNmzpyJqKgodOvWDb169cLXX3+N69evY8qUKQBK09bj4uLEavtTpkzBl19+iZkzZ2Ly5MmIjo7GsmXLxCr1QGkxv379+mH+/PkYPXo0fv/9d2zduhV79uyp9n5VKhWmT5+ODz74AK1atUKrVq3wwQcfwN3dHZMmTbLZ70NqbgbTTTTz95CxJcbMnTh6G6TPhJpIFVapVGjfyMdm7SKSktqgB71Py/r48alITPx6v8VtnLVqzL2vHVQqYE3MTbwyvA36tW4ANycNBrYJsHWTyU4Y3ggFgKheTbDywHWxsKIcwbz+Lq4KsBg0fB0VgcTMfPT/eIe4PgAMCQ9EfQ9ndGnsi+HtSzPqejavj7+nmb8xQGRPhrULxIzBrfH51vLxzr1NpBW7aDX4+f963tG+nLXlfYO5hbVblJjk1791Axx5cwhctGq0e3sTAKBjiE+FIcul/0rfM1++Q193J/z8fz3xx/FbmFY21WKonxsSMwswuV9zAIBWo8a3j3WTtpE2VqNgXifT+LgJEyYgNTUV7777LuLj49G+fXts2LABTZo0AQDEx8cbzf3erFkzbNiwATNmzMBXX32F4OBgLFy4EGPHjhXX6d27N1atWoU33ngDb775Jlq0aIHVq1cjMjKy2vsFgFdeeQV5eXmYOnUq0tLSEBkZic2bN8PLq3Yq4irFQz1C8fNBhY3tN3Pi0GrKTzL1PZUzxp/I1lKzC6FSqfDUXc1wIPY2AKB1oCcuJGaL67x9Xzge6dlErLb9UI/SKVru7dCw8gsSGXB10uDRXk2w80Iydl1IlrVnXn8BueKJ7thxPhkr9l0FUDrtUOfQenB10ogzlJSuX7qBl6sT9r8+CFoOI6E6SqVS4cXBrYyC+fbBPnjqrmZYtqd0mOHzA1vW+n67NPat9dck+emzL6Jn3w0nE5ka+sdSD7oyvHngWZYxNqpTMO5q6Q9fNyeoVEBeUQncnZWTUVzb7O6dTZ06FVOnTjX53IoVKyot69+/P44cqTxuyNC4ceMwbty4Gu8XKP0jnjt3LubOnWvxdeyfPB9WS8y1pZGvGx7o0gherto6/SEm0gsL9ML5xCwxFW5ouyD8+1J/hNZzx5qYm3h9/UkAwPTBrYwCeSJrOWnUeKxXE+y6kIzCEum/EfQXcPrryQFhAejapJ4YzLdt6G1yerlGBoUi+fdPjmDaoFZY+O9FLH+iO9RqFd4cGY43RrQFULsF6qJn340bt/PQMcS31l6TlKehj+liu+JfkowBgqfB8C/DoR91PQao9rubOXNmtV/0s88+q1FjSPnkSqOxpGJbwhuWFrlTqVT4fEJn6RtEJJOfJkdiz6UUo3GRLcrmiX+oRyg8XDTo2rgeQqtR1IyoKvpg2Joii7VFP7xKVX4JKY6JBICSCvPLLXusG/49l4RHezcBkSOZMbgVnujd1KiIsS2qzDf0cTMb6FHdJ8YHUo+ZN9jdXS39Jd23UlQ7mD969Gi11uM0FHWb/uhK/WG1RN+WNkFe+HhcJzT1Z6BCjqm+pwtGd25k8jmVSmX2OaKa0AfzsqTZi3n25cvUahWGhAfiZlpepXoog9oGYlAdqVxMZA2VSlVpKkai2qa/sSp1Z59hMdOXhjpm0d5qB/Pbt2/HlStX0LRpU6jVTE1zVErumVepVOgQwoJ2RERScC6b31fOYF5doQPhm0e7QRAEdiwQEUmovGdeWvr96WdfcERWReWtWrVCSkqK+HjChAlITEys9UaRcqkUOGZef1eOdYyIiKQjb5p9KVOnfQbyRETSEjN3Je7tEwxmNnFUVgXzFQ/Qhg0bkJOTU6sNImVTlX9aZW2HoYpVjYmIyPb0veJyfB2IF3A87xMRyU/mnnlHjuaZL09WKR8zryD6NHtH/iQTETkQged9IiLFkGvMfPl3geOyKphXqVRm5xUkx6CSsSfGHLGqMf8UiYgkJ0dBVJ73iYiUQ75zsf67wHG/DKyaeE8QBDz++ONwcXEBAOTn52PKlCnw8PAwWm/dunW110JSJEVVs+ddOSIiycl57cTzPhGRchiei6UsQlpeDFWS3SmSVcH8Y489ZvT4kUceqdXGkPIp8caXwEHzRESykWXMvP4/PO0TEcnOMHgXBOkuyctDAMf9MrAqmF++fLmt2kF2Rllp9qUc92NMRCQ9Ocerl1cw5pmfiEhuhmdinSBALdG5Wcdq9iyAR9ZR4tR0rGpMRCQfOb4PmJBFRKQchudiKb8TxCFXDvxdwGCerKL/sCipZ17HsZNERJKT9eKJ530iIsUwzJKSMkYQOOaKwTxZp3xqOgVF82VtUTvybTkiIpnIOWbekcdJEhEphlHPvHRfCpzZhME8WUlVHs0rBlNsiIikJ281e17AEREphVGavQw98478VcBgnqwizjMvczsMlRfAc+SPMhGRXOSYZ74Uz/pERPKTOztW7v3LicE8WUXsmFfQoHmBV3VERA6lPCOLJ34iIrkZzzMv3X6ZnctgnqylwAJ44ngZmdtBRORIZJ2ajuMkiYgUQyX3mHnJ9qg8DObJKsqcmq70X17UERFJT5YCeBwnSUSkGHJVs9cxS4vBPFlHiVPTccw8EZH05C2Ap28Dz/tERHKTb555BQUkMmEwT1ZR4mUTqxoTEclHjksp8bwvw76JiMg8KQPs8mlKJdul4jCYpxpR0jzzTLMnIpKenKdcXsARESmHfD3zlffvaBjMk1WUmWZf2hhHnpaCiEgucqQ5lo+Z53mfiEhuco2ZB2MABvNkHSVeOCnpxgIRkaOQdcw8q9kTESmG0blYjqnppNul4jCYJ6uU98wrJ4JmISQiIvnIM2Zehp0SEZFJxrG8HGPmHTcGYDBPVtF/VJR0HVVezZ6IiBwBx8wTESmHYTCtkzBI0OlYDJXBPFmn7MOqpF4RVrMnIpKDfCddjpknIlIOo555GarZO/JXAYN5skp5z7xyonn2zBMRyUeem7u8iUtEpBSyV7OXcJ9Kw2CerKLEavbgmHkiIsnJWgCP0xERESmG4TW4lDFCeTFUx/0yYDBPVtGnNCorlud4GSIiucgyNV3Zv448HRERkRJJmr1btiu1A38VMJgnqyixZ17HHhoiIsnJecoVa6XI2AYiIiqnkqFKdvlQW8f9NmAwT1YxSKKRsRXGBFa/ICKSjSxT05X9y5u4RETKIMeMVxxyxWCerKTED4s+nceRU2yIiKQm5xhF3sQlIlIWtQwzXimpILdcGMxTjSgpzZ535YiIZCTD94GOU5ISESmKOBRXwi8FHYtgM5gn66iUOM982b+OPF6GiMihcDoiIiJFEYtkS9kzz/opDOapZhSV1sIeGiIiyclaAE/fBp74iYiUQeyZlw7rpzCYJyspsZo9P8hERPKRpQAee2OIiBRFLIAn7aB5AI49TSmDebKKIueZF9MtHfeDTEQkNTmvnXgTl4hIWeTo8NNnCjvydwGDebKKInvmBQ6eJCKSi6S9MOI+S//lTVwiImWQ43zMEIDBPFmpfA5J5UTzOn6QiYgkJ2cgzd4YIiJlkaVnnmlaDObJOqryaF4xWAiJiEg+8oyZl2GnRERklv4qXCfhCbp8RivHxWCerKLMMfOlrVE78ieZiMiBsDOGiEhZxOmrJdynjjNaMZgn65Sn0CgpnC/lwJ9jIiLJyVpDhRWMiYgURY5q9hwzz2Ce6gDxg8yLOiIih8Ax80RECiPDPPP6vTnyjV0G81QjSuqXFy/qZG4HEZEjkqMgKqvZExEpS3nPvHT7LO/Qk26fSsNgnqwijodRUDQvsPoFEZFD4Zh5IiJlUclQJbs8BHDcLwMG82QVBRaz5weZiEgGcgbS+jGZPOsTESmDvFPTSbdPpWEwT1ZRYgE8VrIkIpKPHF8HOl7AEREpihwdfhxqy2CerKTInnlWsiQikpwSio4yI4uISBnkGIqr45h5BvNkHfHiTUnRfBlH/iATEclF6q8Dw8wwnveJiJRBLYYIUk5Np++Zd9wvAwbzZBWVDB/Uqug/yI48LQURkaMw7PXhWZ+ISCnkK5KtduCI1oHfOtWEHNNOVIXTUhARSU885Ur8fWC4OyWk+hMRkbwF8NgzT1RdSpyaTvyf436QiYgchVGavYztICKicuV1taScmo5FsBnMk1Xk+KBWhT3zRETSk2vYlXHPvKS7JiIiM2Sdms6BMZgnq8jxQa0Kp6UgInIcxmPmeeYnIlICOc7H5R16jvtdwGCerKLECyf2zBMRSU+u7wPDTAAVr2KIiBRBfx2uk7DHT78nRw4B+DVINaKgjnlOS0FEJCOpM7VYzZ6ISHnkKJKtv3HgyB16dhPMp6WlISoqCj4+PvDx8UFUVBTS09MtbiMIAubOnYvg4GC4ublhwIABOH36tNE6BQUFeOGFF+Dv7w8PDw+MGjUKN2/etGrfx48fx0MPPYTQ0FC4ubmhbdu2+OKLL2rrrSuKMtPsSznyB5mISGpKOOc6cmolEZGS6M/HkoYIZTtz5Omp7SaYnzRpEo4dO4aNGzdi48aNOHbsGKKioixu89FHH+Gzzz7Dl19+iUOHDiEoKAhDhgxBVlaWuM706dOxfv16rFq1Cnv27EF2djZGjhyJkpKSau87JiYGDRo0wI8//ojTp09jzpw5mD17Nr788sva/0XIrPyjopxovnxaCiIikprU3wbsmSciUi5B0jR71s3Syt2A6jh79iw2btyI/fv3IzIyEgDwzTffoFevXjh//jzCwsIqbSMIAhYsWIA5c+ZgzJgxAIDvv/8egYGB+Omnn/DMM88gIyMDy5Ytw//+9z8MHjwYAPDjjz8iNDQUW7duxbBhw6q17yeffNJo382bN0d0dDTWrVuH559/3pa/Gskps2den2LjyB9lIiLHYDRmnqd9IiJFKJ/hRDqsm2UnPfPR0dHw8fERg2kA6NmzJ3x8fLBv3z6T28TGxiIhIQFDhw4Vl7m4uKB///7iNjExMSgqKjJaJzg4GO3btxfXqcm+ASAjIwN+fn4W31dBQQEyMzONfpROPy5dQbE8P8hERDIoHx8p8dR0rGZPRKQ4skxNV7536XaqMHYRzCckJCAgIKDS8oCAACQkJJjdBgACAwONlgcGBorPJSQkwNnZGfXq1bO4jrX7jo6Oxi+//IJnnnnG4vuaN2+eOA7fx8cHoaGhFtdXBPGDqpxwvrySpeN+kImIHAXnmSciUp7y63AJ0+zZoSdvMD937lyoVCqLP4cPHwZgOoVaEIQqU6srPl+dbSquY82+T58+jdGjR+Ott97CkCFDLO5n9uzZyMjIEH9u3LhhcX0lkP5jWjV+kImIZCBDSiUg7bRHRERUPfL0zHPMvKxj5p9//nlMnDjR4jpNmzbFiRMnkJiYWOm55OTkSj3vekFBQQBKe9YbNmwoLk9KShK3CQoKQmFhIdLS0ox655OSktC7d29xneru+8yZM7j77rsxefJkvPHGGxbfF1Ca9u/i4lLlekoiVqpU0LUUP8hERI7DKM2eJ34iIkWQo8OPHXoy98z7+/ujTZs2Fn9cXV3Rq1cvZGRk4ODBg+K2Bw4cQEZGhhh0V9SsWTMEBQVhy5Yt4rLCwkLs3LlT3CYiIgJOTk5G68THx+PUqVPiOtXd9+nTpzFw4EA89thj+M9//lM7vyAFUmLPPPhBJiKSnGxDmzhmnohIceTo8NMP++XUdArXtm1bDB8+HJMnT8b+/fuxf/9+TJ48GSNHjjSqZN+mTRusX78eQOkf1PTp0/HBBx9g/fr1OHXqFB5//HG4u7tj0qRJAAAfHx889dRTeOmll/Dvv//i6NGjeOSRR9ChQwexun119q0P5IcMGYKZM2ciISEBCQkJSE5Olvg3ZXsqBY6Z16dcspo9EZH0pP46YDV7IiLlkaMoqlg3y4G/C+xiajoAWLlyJaZNmyZWnh81alSledzPnz+PjIwM8fErr7yCvLw8TJ06FWlpaYiMjMTmzZvh5eUlrvP5559Dq9Vi/PjxyMvLw6BBg7BixQpoNJpq7/vXX39FcnIyVq5ciZUrV4rLmzRpgqtXr9bq70FuSvywcJ55IiLpyfV9wHnmiYgUSM6p6Rz428Bugnk/Pz/8+OOPFtepeCdIpVJh7ty5mDt3rtltXF1dsWjRIixatKjG+65qH2Rb4lF33M8xEZHDMK5mzxM/EZESqGVMs3fkGMAu0uxJOcR55pWTZS+2xZHHyxAROQrDG/dqnvaJiBShvK6WDGn2ku1ReRjMk1XEMfMKKoHHavZERNIzPOfKMUYSYM88EZFSqGSokl1ezd5xvwsYzFONKLFn3oE/x0REDkNJ3z9ERFRKn72rk3Seef2+HReDebKKEueZ13Pk4hdERFIz7AmRdIykPhuLp3wiIsWQI3u3fGo6yXapOAzmySpyjIepiiDwwo6IyGGw3hERkWJJWwCv9F+m2RNVU/k88/K2wxBTbIiIpGc0Zl7C/ZbPK8yzPhGRUojZuxLuk3WzGMyTlcRq9jK3w5COg+aJiBwGZyIiIlIeMXtXyoKo7NFjME/WkaNSZVV4YUdEJD257p9yzDwRkfKUj5mXTnks77hfCAzmySqKHDNf9i8v7IiI5CFlT4xOvIHLkz4RkVLIOzWddPtUGgbzZBVFjpnnhR0RkcMQmI5FRKQ45UNxJUyz55h5BvNkLSV+XDgtBRGR1AxvoEqaVslYnohIceTo8NPvS+3AXfMM5qlGFNQxzxQbIiIHxHM+EZFylBfAk26fnJ6awTxZqfyum3LCec4xSUQkA4NTrizzCrNvnohIOWSYmq7Crh0Sg3myigKL2SuqGB8REdkWq9kTESmPHFPTEYN5spK+91tJn1Om2RMRSc/wnCtpwSOOkSQiUhxZpqZTUDwiFwbzZBUl9sxzmiIiIsdRPq8wEREphRxj5ivv3fEwmCerqOT9pJrElEsiIsfBqemIiJRHJcdE88RgnqwjRwpNlXhdR0QkOcNzrqQF8Ezsn4iI5CXHOVlR8YhMGMyTVfSp7ArqmC+/sOOVHRFRnccZTIiIyJAjfx0wmCfriD3zyonm9SmXLIZERCQd+YJpDq0iIlKa8umr5W2Ho2EwT1ZR4JB5Bd1WICIiW+OQeSIiApQVj8iFwTxZRdlT0/HSjohIKrKPmec5n4iI4Ng3dxnMk91jMSQiIsehK7tzwHM+EZFyiHW1ZG6Ho2EwT1ZR4qQT4oUdr+yIiCRjeM6Vso5KeTaWZLskIiIFUlINL7kwmCerlBe3UNCHh+MniYgcRvnXD8/6RETk2Dd3GcyTVVQKvHgSxMrGymsbERHVLoHV7ImIlIfV7GXBYJ6sosRpJ5hySUQkPcObu5IWwGM2FhERQVnxiFwYzJNVysfMK+fTwws7IiLHo+YdXCIigjIzh6XCYJ6so8SeeabZExFJzrgAnnSYjUVEpDxK7PBzBAzmySpKnHaCF3ZERI5DvIErczuIiEheSopH5MJgnqyixGr25fPM89KOiEgOUn4nlN/A5TmfiIgcu0OPwTxZRYnzzLNnnojIcSjp+4eIiEopsUi2I2AwT1YRe0IU9EHV9wgxliciko58Y+Y5NR0REYF3DsBgnqykxIsnMc1egW0jIqLaxXM+EREZcuSvAwbzVCNKug9W3jPvyB9lIiL5yDPPPM/5RERKocQi2Y6AwTxZRYmXTuJJQ4mNIyKqo+QLpplmT0REvHEAMJinOqC8l4aIiOo6Hc/5RERkwJFnN2EwTzWixKnp1A78QSYikprRKVeONHue84mIFEOO6asVFI7IhsE8WUWJ106sbExE5Dg4gwkREVEpBvNUZzCYJyKSjnHHvIQ9MaYaQERE5IAYzFONKCmrhZWNiYgch/6cz6FVRETKIccpWcobyUrFYJ6spLyLJx3T7ImIJGc4Zl3SqenANHsiIirnyDEAg3myeyx+QUTkQMQCePI2g4iISG4M5qlGlBRAi700vLIjIpKFlF8J+n1xaBURkXLoz8mSZmopKB6RC4N5sooS42XOM09EJD25zrkCe+aJiMiAI9/cZTBPdk/spXHczzERkcNgwSMiIqJSDOapRhR1McXKxkREkjM85QoS5jqW98zznE9EpBT6U7IsU5U6MAbzZBUlXjqxsjERkeMoHzNPRETk2Nm5DObJ7nH8JBGR9IymppNwvwKnIyUiIgLAYJ5qSEnVI8ubwis7IqK6jnVSiIiUi9XspcVgnqyixDGKOvbSEBHJStILKnEGE570iYjIsbvzGMyT3ePUdEREjoM3cImIiEoxmKcaUVJaS3nKJa/siIjkIGn1Yt7AJSJSHP11uKRp9qxnz2CerKPIiyeB1eyJiOQgxz1U3sAlIiJDjvx1wGCe7B6LIREROQ5WsyciIirFYJ7sHqemIyKSh3jalTStssK+iYhIdrKck5llbz/BfFpaGqKiouDj4wMfHx9ERUUhPT3d4jaCIGDu3LkIDg6Gm5sbBgwYgNOnTxutU1BQgBdeeAH+/v7w8PDAqFGjcPPmzRrvOzU1FSEhIVCpVFW2zx4pMWDWj5dhyiURUd1XfgOX53wiIqWRI7525O8DuwnmJ02ahGPHjmHjxo3YuHEjjh07hqioKIvbfPTRR/jss8/w5Zdf4tChQwgKCsKQIUOQlZUlrjN9+nSsX78eq1atwp49e5CdnY2RI0eipKSkRvt+6qmn0LFjx9p501QtLIZERCQPseCRpHtlnRQiIiIA0MrdgOo4e/YsNm7ciP379yMyMhIA8M0336BXr144f/48wsLCKm0jCAIWLFiAOXPmYMyYMQCA77//HoGBgfjpp5/wzDPPICMjA8uWLcP//vc/DB48GADw448/IjQ0FFu3bsWwYcOs2veSJUuQnp6Ot956C//880+V76ugoAAFBQXi48zMzJr/kiQmKKicvY69NEREDoNDq4iIlEfOgqiOzC565qOjo+Hj4yMG0wDQs2dP+Pj4YN++fSa3iY2NRUJCAoYOHSouc3FxQf/+/cVtYmJiUFRUZLROcHAw2rdvL65T3X2fOXMG7777Ln744Qeo1dX7tc6bN09M3ffx8UFoaGi1tpOTSoF9IQKr2RMRyUraqYhKKfH7iIjI0cnR4efI3wZ2EcwnJCQgICCg0vKAgAAkJCSY3QYAAgMDjZYHBgaKzyUkJMDZ2Rn16tWzuE5V+y4oKMBDDz2Ejz/+GI0bN672+5o9ezYyMjLEnxs3blR7W6qMvTRERNKS47QrsAIeERERAJmD+blz50KlUln8OXz4MADTKdSCIFSZWl3x+epsU3GdqvY9e/ZstG3bFo888ojF163IxcUF3t7eRj/2QklpLeVj5nllR0QkB0HCbwWBY+aJiBRHnpu7SopI5CHrmPnnn38eEydOtLhO06ZNceLECSQmJlZ6Ljk5uVLPu15QUBCA0p71hg0bisuTkpLEbYKCglBYWIi0tDSj3vmkpCT07t1bXKeqfW/btg0nT57EmjVrAJT/Yfn7+2POnDl45513LL5He6LE3u/yavYyN4SIyMHIMkaSY+aJiBRLlvDagb8PZA3m/f394e/vX+V6vXr1QkZGBg4ePIgePXoAAA4cOICMjAwx6K6oWbNmCAoKwpYtW9ClSxcAQGFhIXbu3In58+cDACIiIuDk5IQtW7Zg/PjxAID4+HicOnUKH330UbX3vXbtWuTl5Yn7PnToEJ588kns3r0bLVq0qMmvhqzAavZERPLimHkiIiLp2UU1+7Zt22L48OGYPHky/vvf/wIA/u///g8jR440qibfpk0bzJs3Dw888ABUKhWmT5+ODz74AK1atUKrVq3wwQcfwN3dHZMmTQIA+Pj44KmnnsJLL72E+vXrw8/PDy+//DI6dOggVrevzr4rBuwpKSnitr6+vjb93ZDBHUBe1xERSao0oJZ4YjqB2VhEREojx6xSzLK3k2AeAFauXIlp06aJledHjRqFL7/80mid8+fPIyMjQ3z8yiuvIC8vD1OnTkVaWhoiIyOxefNmeHl5iet8/vnn0Gq1GD9+PPLy8jBo0CCsWLECGo3Gqn07GiV9ePQXdmpe2RER1XlMsyciUjAZYgRHztSym2Dez88PP/74o8V1KhZBUKlUmDt3LubOnWt2G1dXVyxatAiLFi26o30bGjBgAAsySIiFjYmIZFJ24pXyG09fJ4U3cImIHBujLTuZmo7IkvJeGl7YERHVdbxXTkSkPHJehTtyCMBgnmpEymmIqsLxk0RE8pIyG403cImIlEtJMYIjYDBPVlHitRPT7ImI5CHLvMIy7puIiJSDmVoM5qkOYDEkIiJ5STo1HbOxiIgUR85zsiN/HTCYpxpR0p0wgf00RESykOPijWd8IiLlUlKM4AgYzJNVlDj1A3vmiYgcCMfMExEROD4fYDBPdYAYzMvbDCIihyPHDV79xRvP+URESiLfWdmR7+0ymKcaUeJ9MM45TERU9zEbi4hIuZQYI9RlDObJKkq8eGIxJCIieejPu5IWwCvfu3Q7JSIixeH4fAbzVAfoxDR7XtgREdV17JknIlIeeavZO+4XAoN5qhEl3QkTx0867ueYiEhWUhYh4ph5IiLlUlKM4AgYzJNVlBgw86RBRCQPOb4S2DNPRERUisE82T1xzmFe2BERyULSMfNlO2PRUyIi5ZDzjOzIXwcM5qmGlNMdLnDMPBGRLOSY6503cImIlItzv0uLwTxZRZkBM8fMExHJScpLN97AJSIioDxTy5ExmCe7x/GTRETykGfMvBjNExGRQshbzd5xMZinGlHSjTB9Uzh+koio7hPT7GVtBRERmaKkGMERMJgnqygxXtb30iiwaUREdVvZiVfKVMfybCye9YmIHBnvGzCYpzpAxzR7IiKHwZ55IiLlkbWOiQMHAQzmqUaUdCesvEfIcT/IRERykrYAHoueEhEplZJiBEfAYJ6sosRrJ05TREQkD1nnFZZx30REJD+Oz2cwT3UBCxsTEclKygsqjpknIlIeVrOXB4N5qhElzetY3jPvyB9lIiLpyXHeFcCip0REiqWgGMERMJgnqygxXmY1eyIiuUlfzZ4nfSIixyZwhD6DebJ/HDNPRCQPOc675bE8T/pEREoha5q9A38daOVuAJXT9zBnZmbK3BLzsrKyoCvIRbFaq5h2FufnQFekQ3ZWFjKdSuRuDhGRwyjJz4GuoBhZmZnIdJWmhyQ3u/R7qDAvWzHfQ0REjq4wNxu6glzk5mRJdm4uKNtnfk7d+z7Qv5+qhjarBCUNfnZwN2/eRGhoqNzNICIiIiIiIpnduHEDISEhZp9nMK8gOp0O58+fR3h4OG7cuAFvb2+5m0RWyMzMRGhoKI+dHeKxs188dvaLx85+8djZLx47+8TjZr9qeuwEQUBWVhaCg4OhVpsfGc80ewVRq9Vo1KgRAMDb25sfVjvFY2e/eOzsF4+d/eKxs188dvaLx84+8bjZr5ocOx8fnyrXYQE8IiIiIiIiIjvDYJ6IiIiIiIjIzjCYVxgXFxe8/fbbcHFxkbspZCUeO/vFY2e/eOzsF4+d/eKxs188dvaJx81+2frYsQAeERERERERkZ1hzzwRERERERGRnWEwT0RERERERGRnGMwTERERERER2RkG80RERERERER2hsG8AqSlpSEqKgo+Pj7w8fFBVFQU0tPTLW6TnZ2N559/HiEhIXBzc0Pbtm2xZMkSaRpMopocOwA4e/YsRo0aBR8fH3h5eaFnz564fv267RtMAGp+3PSeeeYZqFQqLFiwwGZtJNOsPXZFRUV49dVX0aFDB3h4eCA4OBiPPvoobt26JV2jHdTixYvRrFkzuLq6IiIiArt377a4/s6dOxEREQFXV1c0b94cS5culailVJE1x27dunUYMmQIGjRoAG9vb/Tq1QubNm2SsLVkyNrPnd7evXuh1WrRuXNn2zaQzLL22BUUFGDOnDlo0qQJXFxc0KJFC3z33XcStZYMWXvsVq5ciU6dOsHd3R0NGzbEE088gdTU1JrtXCDZDR8+XGjfvr2wb98+Yd++fUL79u2FkSNHWtzm6aefFlq0aCFs375diI2NFf773/8KGo1G+O233yRqNQlCzY7dpUuXBD8/P2HWrFnCkSNHhMuXLwt//fWXkJiYKFGrqSbHTW/9+vVCp06dhODgYOHzzz+3bUOpEmuPXXp6ujB48GBh9erVwrlz54To6GghMjJSiIiIkLDVjmfVqlWCk5OT8M033whnzpwRXnzxRcHDw0O4du2ayfWvXLkiuLu7Cy+++KJw5swZ4ZtvvhGcnJyENWvWSNxysvbYvfjii8L8+fOFgwcPChcuXBBmz54tODk5CUeOHJG45WTtsdNLT08XmjdvLgwdOlTo1KmTNI0lIzU5dqNGjRIiIyOFLVu2CLGxscKBAweEvXv3SthqEgTrj93u3bsFtVotfPHFF8KVK1eE3bt3C+3atRPuv//+Gu2fwbzMzpw5IwAQ9u/fLy6Ljo4WAAjnzp0zu127du2Ed99912hZ165dhTfeeMNmbSVjNT12EyZMEB555BEpmkgm1PS4CYIg3Lx5U2jUqJFw6tQpoUmTJgzmJXYnx87QwYMHBQBVXuBSzfXo0UOYMmWK0bI2bdoIr732msn1X3nlFaFNmzZGy5555hmhZ8+eNmsjmWbtsTMlPDxceOedd2q7aVSFmh67CRMmCG+88Ybw9ttvM5iXibXH7p9//hF8fHyE1NRUKZpHFlh77D7++GOhefPmRssWLlwohISE1Gj/TLOXWXR0NHx8fBAZGSku69mzJ3x8fLBv3z6z29111134448/EBcXB0EQsH37dly4cAHDhg2TotmEmh07nU6Hv//+G61bt8awYcMQEBCAyMhI/PbbbxK1mmr6mdPpdIiKisKsWbPQrl07KZpKFdT02FWUkZEBlUoFX19fG7SSCgsLERMTg6FDhxotHzp0qNnjFB0dXWn9YcOG4fDhwygqKrJZW8lYTY5dRTqdDllZWfDz87NFE8mMmh675cuX4/Lly3j77bdt3UQyoybH7o8//kC3bt3w0UcfoVGjRmjdujVefvll5OXlSdFkKlOTY9e7d2/cvHkTGzZsgCAISExMxJo1azBixIgatYHBvMwSEhIQEBBQaXlAQAASEhLMbrdw4UKEh4cjJCQEzs7OGD58OBYvXoy77rrLls0lAzU5dklJScjOzsaHH36I4cOHY/PmzXjggQcwZswY7Ny509ZNJtT8Mzd//nxotVpMmzbNls0jC2p67Azl5+fjtddew6RJk+Dt7V3bTSQAKSkpKCkpQWBgoNHywMBAs8cpISHB5PrFxcVISUmxWVvJWE2OXUWffvopcnJyMH78eFs0kcyoybG7ePEiXnvtNaxcuRJarVaKZpIJNTl2V65cwZ49e3Dq1CmsX78eCxYswJo1a/Dcc89J0WQqU5Nj17t3b6xcuRITJkyAs7MzgoKC4Ovri0WLFtWoDQzmbWTu3LlQqVQWfw4fPgwAUKlUlbYXBMHkcr2FCxdi//79+OOPPxATE4NPP/0UU6dOxdatW232nhyFLY+dTqcDAIwePRozZsxA586d8dprr2HkyJEs9nSHbHncYmJi8MUXX2DFihUWP5dUM7Y+X+oVFRVh4sSJ0Ol0WLx4ca2/DzJW8ZhUdZxMrW9qOdmetcdO7+eff8bcuXOxevVqkzfeyPaqe+xKSkowadIkvPPOO2jdurVUzSMLrPnc6XQ6qFQqrFy5Ej169MC9996Lzz77DCtWrGDvvAysOXZnzpzBtGnT8NZbbyEmJgYbN25EbGwspkyZUqN98zacjTz//POYOHGixXWaNm2KEydOIDExsdJzycnJle7y6OXl5eH111/H+vXrxZSMjh074tixY/jkk08wePDgO38DDsyWx87f3x9arRbh4eFGy9u2bYs9e/bUvNFk0+O2e/duJCUloXHjxuKykpISvPTSS1iwYAGuXr16R213dLY8dnpFRUUYP348YmNjsW3bNvbK25C/vz80Gk2lXomkpCSzxykoKMjk+lqtFvXr17dZW8lYTY6d3urVq/HUU0/h119/5XWIDKw9dllZWTh8+DCOHj2K559/HkBpgCgIArRaLTZv3oy7775bkrY7upp87ho2bIhGjRrBx8dHXNa2bVsIgoCbN2+iVatWNm0zlarJsZs3bx769OmDWbNmASiN4Tw8PNC3b1+8//77aNiwoVVtYDBvI/7+/vD3969yvV69eiEjIwMHDx5Ejx49AAAHDhxARkYGevfubXKboqIiFBUVQa02TqzQaDRizy/VnC2PnbOzM7p3747z588bLb9w4QKaNGly5413YLY8blFRUZUuTocNG4aoqCg88cQTd954B2fLYweUB/IXL17E9u3bGRzamLOzMyIiIrBlyxY88MAD4vItW7Zg9OjRJrfp1asX/vzzT6NlmzdvRrdu3eDk5GTT9lK5mhw7oLRH/sknn8TPP/9c43GfdGesPXbe3t44efKk0bLFixdj27ZtWLNmDZo1a2bzNlOpmnzu+vTpg19//RXZ2dnw9PQEUHotqVarERISIkm7qWbHLjc3t9KwFo1GA6A8I80qNSqbR7Vq+PDhQseOHYXo6GghOjpa6NChQ6WplsLCwoR169aJj/v37y+0a9dO2L59u3DlyhVh+fLlgqurq7B48WKpm+/QanLs1q1bJzg5OQlff/21cPHiRWHRokWCRqMRdu/eLXXzHVZNjltFrGYvD2uPXVFRkTBq1CghJCREOHbsmBAfHy/+FBQUyPEWHIJ+qp5ly5YJZ86cEaZPny54eHgIV69eFQRBEF577TUhKipKXF8/Nd2MGTOEM2fOCMuWLePUdDKx9tj99NNPglarFb766iujz1d6erpcb8FhWXvsKmI1e/lYe+yysrKEkJAQYdy4ccLp06eFnTt3Cq1atRKefvppud6Cw7L22C1fvlzQarXC4sWLhcuXLwt79uwRunXrJvTo0aNG+2cwrwCpqanCww8/LHh5eQleXl7Cww8/LKSlpRmtA0BYvny5+Dg+Pl54/PHHheDgYMHV1VUICwsTPv30U0Gn00nbeAdXk2MnCIKwbNkyoWXLloKrq6vQqVMn4bfffpOu0VTj42aIwbw8rD12sbGxAgCTP9u3b5e8/Y7kq6++Epo0aSI4OzsLXbt2FXbu3Ck+99hjjwn9+/c3Wn/Hjh1Cly5dBGdnZ6Fp06bCkiVLJG4x6Vlz7Pr372/y8/XYY49J33Cy+nNniMG8vKw9dmfPnhUGDx4suLm5CSEhIcLMmTOF3NxciVtNgmD9sVu4cKEQHh4uuLm5CQ0bNhQefvhh4ebNmzXat0oQatKfT0RERERERERyYTV7IiIiIiIiIjvDYJ6IiIiIiIjIzjCYJyIiIiIiIrIzDOaJiIiIiIiI7AyDeSIiIiIiIiI7w2CeiIiIiIiIyM4wmCciIiIiIiKyMwzmiYiIiIiIiOwMg3kiIiKqsblz56Jz586S73fHjh1QqVRQqVS4//77La47YMAATJ8+3eixfttjx47ZtJ1ERES2opW7AURERKRMKpXK4vOPPfYYvvzyS7zwwgsStaiy8+fPIyAgwKpt1q1bh8uXL6NHjx42ahUREZHtMZgnIiIik+Lj48X/r169Gm+99RbOnz8vLnNzc4Onpyc8PT3laB4AICAgAL6+vlZt4+fnh8zMTNs0iIiISCJMsyciIiKTgoKCxB8fHx+oVKpKyyqm2T/++OO4//778cEHHyAwMBC+vr545513UFxcjFmzZsHPzw8hISH47rvvjPYVFxeHCRMmoF69eqhfvz5Gjx6Nq1evWt3mnJwcPProo/D09ETDhg3x6aef3uFvgYiISJkYzBMREVGt2rZtG27duoVdu3bhs88+w9y5czFy5EjUq1cPBw4cwJQpUzBlyhTcuHEDAJCbm4uBAwfC09MTu3btwp49e+Dp6Ynhw4ejsLDQqn3PmjUL27dvx/r167F582bs2LEDMTExtnibREREsmIwT0RERLXKz88PCxcuRFhYGJ588kmEhYUhNzcXr7/+Olq1aoXZs2fD2dkZe/fuBQCsWrUKarUa3377LTp06IC2bdti+fLluH79Onbs2FHt/WZnZ2PZsmX45JNPMGTIEHTo0AHff/89SkpKbPROiYiI5MMx80RERFSr2rVrB7W6vL8gMDAQ7du3Fx9rNBrUr18fSUlJAICYmBhcunQJXl5eRq+Tn5+Py5cvV3u/ly9fRmFhIXr16iUu8/PzQ1hYWE3fChERkWIxmCciIqJa5eTkZPRYpVKZXKbT6QAAOp0OERERWLlyZaXXatCgQbX3KwhCDVpLRERknxjMExERkay6du2K1atXIyAgAN7e3jV+nZYtW8LJyQn79+9H48aNAQBpaWm4cOEC+vfvX1vNJSIiUgSOmSciIiJZPfzww/D398fo0aOxe/duxMbGYufOnXjxxRdx8+bNar+Op6cnnnrqKcyaNQv//vsvTp06hccff9wo5Z+IiKiuYM88ERERycrd3R27du3Cq6++ijFjxiArKwuNGjXCoEGDrO6p//jjj5GdnY1Ro0bBy8sLL730EjIyMmzUciIiIvmoBA4wIyIiIjuzY8cODBw4EGlpafD19bV6+6tXr6JZs2Y4evQoOnfuXOvtIyIisjXmnREREZHdCgkJwUMPPWTVNvfccw/atWtnoxYRERFJgz3zREREZHfy8vIQFxcHoHSsfFBQULW3jYuLQ15eHgCgcePGcHZ2tkkbiYiIbInBPBEREREREZGdYZo9ERERERERkZ1hME9ERERERERkZxjMExEREREREdkZBvNEREREREREdobBPBEREREREZGdYTBPREREREREZGcYzBMRERERERHZGQbzRERERERERHbm/wHtfZWjlBFahwAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "lcs = LCSim(window_width=24*1.5, exp_time=60, white_noise=1e-5)\n", "\n", "time, flux = lcs(radius_ratio=0.1, zero_epoch=0.0, period=0.7,\n", " scaled_semi_major_axis=4.0, impact_parameter=0.1,\n", " eccentricity=0.2, argument_of_periastron=0.23*pi,\n", " geometric_albedo=0.5, limb_darkening=[0.2, 0.3],\n", " efficiency=0.6, eff_phase=0.2)\n", "\n", "fig, ax = plt.subplots(1, 1, figsize=(10,4), constrained_layout=True)\n", "ax.plot(time, flux, '-')\n", "plt.setp(ax, ylim=(0.9995, 1.0005), xlabel='Time [d]', ylabel='Flux');" ] }, { "cell_type": "markdown", "id": "fb6f7b9a4e448215", "metadata": {}, "source": [ "---\n", "
© 2025 Judith Korth and Hannu Parviainen
" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "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.12.3" } }, "nbformat": 4, "nbformat_minor": 5 }