{ "cells": [ { "cell_type": "markdown", "id": "41e312b1", "metadata": {}, "source": [ "# PhenomSourceCalculator" ] }, { "cell_type": "markdown", "id": "d954ab9d-7605-4f31-86de-9ed93a291550", "metadata": {}, "source": [ "This example shows how to create a SASE pulse with a phenomenological model described in this paper by Trey Guest et al.: https://doi.org/10.1107/S2052252523008242" ] }, { "cell_type": "code", "execution_count": 1, "id": "f4ba1c7a-ee56-4207-9b8c-7a877b74fb09", "metadata": {}, "outputs": [], "source": [ "# import sys\n", "# sys.path.append(\"xxx\")\n", "# import wpg" ] }, { "cell_type": "code", "execution_count": 2, "id": "3b6732ec", "metadata": {}, "outputs": [], "source": [ "from SimExLite.SourceCalculators import PhenomSourceCalculator\n", "import h5py as h5\n", "from matplotlib import pyplot as plt\n", "from wpg.wpg_uti_wf import get_intensity_on_axis\n", "import numpy as np" ] }, { "cell_type": "code", "execution_count": 3, "id": "be1acfee-c960-40f8-b5d7-8c8a5e2aaf6a", "metadata": {}, "outputs": [], "source": [ "def analytical_pulse_width(photon_energy):\n", " \"\"\"\n", " Estimate analytical_pulse_width (FWHM) from radiation energy (assumes symmetrical beam)\n", "\n", " :param photon_energy: radiation energy [eV]\n", "\n", " :return sig: Radiation pulse width [m]\n", " \"\"\"\n", "\n", " sig = np.log((7.4e03/(photon_energy/1e03)))*6\n", " return sig/1e6\n", "\n", "\n", "def analytical_pulse_divergence(photon_energy):\n", " \n", " \"\"\"\n", " Estimate of analytical_pulse_divergence (half-angle) from radiation energy\n", "\n", " :param photon_energy: radiation energy [eV]\n", "\n", " :return dtheta: pulse divergence [rad]\n", " \"\"\"\n", " return ((14.1)/((photon_energy/1e03)**0.75)) / 1e06" ] }, { "cell_type": "markdown", "id": "d07e5d0f", "metadata": {}, "source": [ "## Construct Calculator" ] }, { "cell_type": "code", "execution_count": 4, "id": "fd288956", "metadata": {}, "outputs": [ { "data": { "text/plain": [ " - Parameters object -\n", "range_x [-0.0002 0.0002] [meter] The spacial mesh range in x direction. [start, end] \n", "num_x 512 Number of mesh points in x direction. \n", "range_y [-0.0002 0.0002] [meter] The spacial mesh range in y direction. [start, end] \n", "num_y 512 Number of mesh points in y direction. \n", "range_t [-2.5e-14 2.5e-14] [second] The temporal range. [start, end] \n", "num_t 250 Number of mesh points in t direction. \n", "photon_energy 10000.0 [electron_volt] The photon energy of X-ray beam. \n", "pulse_energy 0.0001 [joule] Total energy of the pulse \n", "pulse_duration 1.5e-14 [second] The length of a pulse \n", "spectral_bandwidth 1e-12 The bandwith of the beam spectrum \n", "sigma 5e-05 pulse width \n", "div 0.0025 pulse divergence " ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "psc = PhenomSourceCalculator(\"psc\",)\n", "psc.parameters['spectral_bandwidth'] = 1e-12\n", "psc.parameters" ] }, { "cell_type": "code", "execution_count": 5, "id": "473ca08f-b981-4dd5-a6da-63b917d313c5", "metadata": {}, "outputs": [ { "data": { "text/plain": [ " - Parameters object -\n", "range_x [-0.0002 0.0002] [meter] The spacial mesh range in x direction. [start, end] \n", "num_x 512 Number of mesh points in x direction. \n", "range_y [-0.0002 0.0002] [meter] The spacial mesh range in y direction. [start, end] \n", "num_y 512 Number of mesh points in y direction. \n", "range_t [-2.5e-14 2.5e-14] [second] The temporal range. [start, end] \n", "num_t 250 Number of mesh points in t direction. \n", "photon_energy 10000.0 [electron_volt] The photon energy of X-ray beam. \n", "pulse_energy 0.0001 [joule] Total energy of the pulse \n", "pulse_duration 1.5e-14 [second] The length of a pulse \n", "spectral_bandwidth 1e-12 The bandwith of the beam spectrum \n", "sigma 3.963990111718929e-05 pulse width \n", "div 2.507373968154881e-06 pulse divergence " ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# Set pulse width and pulse divergence according to photon_energy\n", "psc.parameters[\"sigma\"] = analytical_pulse_width(psc.parameters[\"photon_energy\"].value)\n", "psc.parameters[\"div\"] = analytical_pulse_divergence(psc.parameters[\"photon_energy\"].value)\n", "psc.parameters" ] }, { "cell_type": "code", "execution_count": 6, "id": "952381f2", "metadata": { "scrolled": true }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Data collection:\n", "key - mapping\n", "\n", "Phenom_wavefront - : PhenomSourceCalculator/wavefront.h5\n", "\n" ] } ], "source": [ "output = psc.backengine()\n", "print(output)" ] }, { "cell_type": "markdown", "id": "fa11c33e-3083-4953-9258-2d2f6f00ae83", "metadata": {}, "source": [ "## Get energy spectrum and temporal profile from the output" ] }, { "cell_type": "code", "execution_count": 7, "id": "18f7acaf", "metadata": {}, "outputs": [], "source": [ "from wpg.wavefront import Wavefront\n", "from wpg.wpg_uti_wf import plot_intensity_map,plot_t_wf, plot_wf, integral_intensity, get_intensity_on_axis\n", "from wpg.srw import srwlpy\n", "\n", "wfr = Wavefront()\n", "wfr.load_hdf5(output[\"Phenom_wavefront\"].filename)\n", "# srwlpy.SetRepresElecField(wfr._srwl_wf, 'time')" ] }, { "cell_type": "code", "execution_count": 8, "id": "c75cca40-f61c-4685-b482-b7149ca8557d", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjcAAAGwCAYAAABVdURTAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/SrBM8AAAACXBIWXMAAA9hAAAPYQGoP6dpAACQDUlEQVR4nO29d3gc5dX+f89WFVuSq+Qid1OMMTYYjCEESIwdIA6kEoeA6SSBJOD3DSUhEMiXHkryixOHEErgJbQQAoGYYkw3zdg0Y4O7sS3Zlq0ubZv5/bH7zDwzO7s7s9rV7kr357p0Sdppz2olzb3n3OccRdM0DYQQQgghfQRPoRdACCGEEJJLKG4IIYQQ0qeguCGEEEJIn4LihhBCCCF9CoobQgghhPQpKG4IIYQQ0qeguCGEEEJIn8JX6AX0NqqqYseOHRg4cCAURSn0cgghhBDiAE3T0NbWhpEjR8LjSR+b6XfiZseOHaivry/0MgghhBCSBdu2bcPo0aPT7tPvxM3AgQMBxH84VVVVBV4NIYQQQpzQ2tqK+vp6/T6ejn4nbkQqqqqqiuKGEEIIKTGcWEpoKCaEEEJIn4LihhBCCCF9CoobQgghhPQpKG4IIYQQ0qeguCGEEEJIn4LihhBCCCF9CoobQgghhPQpKG4IIYQQ0qeguCGEEEJIn4LihhBCCCF9ioKKm1dffRXz58/HyJEjoSgKnnzyyYzHvPzyyzj00EMRDAYxadIk3HfffXlfJyGEEEJKh4KKm46ODhxyyCFYvHixo/03bdqEk08+GccffzxWr16NSy65BOeddx6ee+65PK+UEEIIIaVCQQdnnnjiiTjxxBMd779kyRKMHz8et912GwDgwAMPxOuvv4477rgD8+bNy9cyCSFFQFc4hvKAt9DLIISUACXluVmxYgXmzJljemzevHlYsWJFymNCoRBaW1tNH4SQ0uKWpWtxyHXPY80O/v0SQjJTUuKmoaEBtbW1psdqa2vR2tqKrq4u22NuvPFGVFdX6x/19fW9sVRCSA55d/NehKMq1jZQ3BBCMlNS4iYbrrzySrS0tOgf27ZtK/SSCCEuaeuOAgBiqlbglRBCSoGCem7cUldXh8bGRtNjjY2NqKqqQnl5ue0xwWAQwWCwN5ZHCMkT7aG4uNGobQghDiipyM3s2bOxbNky02MvvPACZs+eXaAVEUJ6AyFuYlQ3hBAHFFTctLe3Y/Xq1Vi9ejWAeKn36tWrsXXrVgDxlNKZZ56p7/+jH/0IGzduxGWXXYa1a9fiT3/6Ex599FFceumlhVg+IaQX0DSNaSlCiCsKKm7ee+89zJgxAzNmzAAALFq0CDNmzMDVV18NANi5c6cudABg/PjxeOaZZ/DCCy/gkEMOwW233Ya7776bZeCE9GG6I6ouajRGbgghDiio5+a4445L+8/Krvvwcccdh1WrVuVxVYSQYqItFNG/ZuSGEOKEkvLcEEL6H+2JlBQAUNsQQpxAcUMIKWqEmRgAVKalCCEOoLghhBQ1bd0UN4QQd1DcEEKKGlncxNQCLoQQUjJQ3BBCihqmpQghbqG4IYQUNe3dRrWUSkcxIcQBFDeEkKLGlJZi5IYQ4gCKG0JIUWNOSxVwIYSQkoHihhBS1LTJ4iaP6iYSU/GvVV9gR3NX3q5BCOkdKG4IIUVNey+Vgj/09lZc+sgHOOH2V/J2DUJI70BxQwgpauS0VD49N+9t2QcA6AjH8nYNQkjvQHFDCClq2qRqqXz6ieuqgtJ1aO4hpJShuCGEFDXmJn75Ex21VWX61y1dkTR7EkKKHYobQkhRY0pL5VHc+DyK/vUX+2gqJqSUobghhBQ1srjJZ7ooKgmn7ayYIqSkobghhBQtmqaZqqXyaSiWK7G2M3JDSElDcUMIKVq6I6opopLPJn6M3BDSd6C4IYQULW0hs7E3n0385HMzckNIaUNxQwgpWuSUFJDfJn6M3BDSd6C4IYQULbKZGABiav6upVLcENJnoLghhBQtbZbITW9VS+3tCKMzHE2zNyGkmKG4IYQULVZxk89qKeu5OUCTkNKF4oYQUrQkp6V6x1AMsJEfIaUMxQ0hpGhp7zZXS+Vz5FPUIm7ouyGkdKG4IYQULYWM3DS0dOftWoSQ/EJxQwgpWtpDMQBAud8LIDel4JEUJVfWyE0kxsnghJQqFDeEkKIlpsaFSMAX/1fVU3Hzj3e2Yuo1z+H1z/ckbbOeO589dQgh+YXihhBStIhgit+rmL7Plnc370UoquLdzXuTtkUtkZp8psAIIfmF4oYQUrSI6InPE/9X1VPBIVJNLV2RpG2iFNwQUhQ3hJQqFDeEkKJF6AuvJzeCIxyNe3hau23EjSrETSIFxsgNISULxQ0hpGhRcxxNEZGb1q7k7sNWcZPPhoGEkPxCcUMIKVr0tJQ3V2mpuEHZUeSG2oaQkoXihhBStAiB4fPkxlAcjibEjZ3nJnHyoI9pKUJKHYobQkjRoumRm4S46aHgCCciN9aZVUByCozVUoSULhQ3hJCiJdHmBl5Pbvrc6Gkpm8hNlGkpQvoMFDeEkKJFj6Yk0lI9bRocicZP0BaKJkVmkj03VDeElCoUN4SQokW1lIJrOYrcAECbxVSsixtfbszLhJDCQXFDCClaNM1Snt1DwRGKGuLGWg4uzh1gEz9CSh6KG0JI0SIEhojc5KoUHEguB2daipC+A8UNIaRosc6W6qneMIkbi6k4ZokSqfbDwwkhJQDFDSGkaEmaLZWjDsVA5sgNOxQTUrpQ3BBCihZ9tlSOfDBhB54bNvEjpPShuCGEFC3WUvCeCA5N0/QmfkC6yA0NxYSUOhQ3hJCiRQgOo4lf9ueKWg5usXpuktJS2V+LEFJYKG4IIUWL1VDck2op2UwMpDEUMy1FSMlDcUMIKVqss6V60sRP9tsAQGt3qj43LAUnpNShuCGEFC25rJYKZ4rcWDw37FBMSOlCcUMIKVqEvvDpTfyyP1fEYqJhEz9C+i4UN4SQokWP3CQER0/SUhFrWipFKTinghNS+lDcEEKKFs1qKM5lWipF5CbAwZmElDwUN4SQosU6W6onFUxJhuKU4xfY54aQUofihhBStKjWeU890BuiFHxg0AcA6AjHTOXhsRg9N4T0FShuCCFFixAzeuSmJ56bhHgZPCCgP9YmlYOLyI2Rlsr6UoSQAkNxQwgpWvQ+N56el2eLtFS534vKgBeAOTUVtRqK6bkhpGShuCGEFC3WUvCeRW7i4sbv9aC63A/AbCpW2cSPkD4DxQ0hpGixloL3JJgiqqUCPg8GliXETVdyWsqYLUVxQ0ipQnFDCClakpv45SJyo+i+GvGYqmpJZedMSxFSuhRc3CxevBjjxo1DWVkZZs2ahXfeeSft/nfeeSf2339/lJeXo76+Hpdeeim6u7t7abWEkN5Es0Ru5MfcIjw3fq9HNygLsSRHafTBmdQ2hJQsBRU3jzzyCBYtWoRrrrkG77//Pg455BDMmzcPu3btst3/oYcewhVXXIFrrrkGn376Kf72t7/hkUcewS9/+cteXjkhpDdQLb1ngOyjNyJKE5DEjTARy+cUnhs28SOkdCmouLn99ttx/vnn4+yzz8aUKVOwZMkSVFRU4J577rHd/80338TRRx+NH/zgBxg3bhzmzp2LBQsWpI32hEIhtLa2mj4IIaWBmijHFmIEyD6iEpb62FhLy2Uhwz43hJQ+BRM34XAYK1euxJw5c4zFeDyYM2cOVqxYYXvMUUcdhZUrV+piZuPGjXj22Wdx0kknpbzOjTfeiOrqav2jvr4+t0+EEJIVL6xpxJdufgkrt+xNuY91Krj8mFvEbKmAzwOvYo7cRE3ihh2KCSl1CiZu9uzZg1gshtraWtPjtbW1aGhosD3mBz/4Aa677jp86Utfgt/vx8SJE3HcccelTUtdeeWVaGlp0T+2bduW0+dBCMmOl9buwhf7uvDKZ3tS7qNa+tzIj7klLJWC+yymYdk8zCZ+hJQ+BTcUu+Hll1/GDTfcgD/96U94//338cQTT+CZZ57Bb3/725THBINBVFVVmT4IIYUnlsg5RdOoCL1aKheeGz1yo8CTJnIjokSM3BBSuvgKdeGhQ4fC6/WisbHR9HhjYyPq6upsj/n1r3+NM844A+eddx4A4OCDD0ZHRwcuuOAC/OpXv4LHU1JajZB+jRAUkbTixtx7BjB8OG6Rm/gZpeWq6TpejwJxKYobQkqXgqmBQCCAww47DMuWLdMfU1UVy5Ytw+zZs22P6ezsTBIwXm+8jXq25aGEkMKg6uIm9d+u+LPOTVrK6EDs0cVNfJsQWl7FiOqwWoqQ0qVgkRsAWLRoERYuXIiZM2fiiCOOwJ133omOjg6cffbZAIAzzzwTo0aNwo033ggAmD9/Pm6//XbMmDEDs2bNwvr16/HrX/8a8+fP10UOIaQ0EIIiFM0cuTGlpbI1FIvIjc8mcqPKkZv4Nr5fIqR0Kai4Oe2007B7925cffXVaGhowPTp07F06VLdZLx161ZTpOaqq66Coii46qqrsH37dgwbNgzz58/H9ddfX6inQAjJkpiLtJRHUaAoccGRdeRGauLnsTbxk8QNIzeElD4FFTcAcPHFF+Piiy+23fbyyy+bvvf5fLjmmmtwzTXX9MLKCCH5xJG4SWzyKAq8ioKopvXYcxPwKnrkxmoo9noUQ/gwdENIyUIHLiGkIDgRN5oUuRERlZ6Wggd8yU38TIZiRaSlKG4IKVUobgghBUFERsLR1CJCZIYUBRAZ6uzHL0gdiq2l4DE5LdWz6xBCCg/FDSGkILj13PQ4chONAUjRxE9EbhQjLaVqjN4QUqpQ3BBCCoKIlqQXN/HPHg/0aEu2AZWIXAqeoomfnJYCWDFFSKlCcUMIKQgiLeXUc6P0MF0UkTw3wlCspqmWktdICCktCl4tRQjpnwhBEU7TxM9IS0HqP+NccERiKq7+98eYPXGobSl41E7cSG/5YqoGP1toEVJyUNwQQgqCPn4hbRO/+GdF7hzsQty8u3kv/vHONry1cS9qKvwA4lO/fWn63Hg9TEsRUuowLUUIKQiqW0Oxnkpyfo1drSEAwN6OsKlDsbWJn8lQzLQUISUPxQ0hpCA4GZwptIVHkQ3FzgXHnva4uGntjiAUEU38POmb+MnihuXghJQkFDeEkIIg5jqlG5xpLgU3P+aEPe1hAHGRtK8z/nXA50kSSnazpeLHUdwQUopQ3BBCCoJhKM6cloo38XM/80lEbgCgqSMubvxeD7wJ17B95CZ5jYSQ0oLihhBSEHRxk85QLM2WyqaJnyxuxGF+r5LUxE82FCty2TkjN4SUJBQ3hJCC4MRzI6elvB73TfxkcSOwa+InixtxPYDVUoSUKhQ3hJCC4K5aClk18dvTFk56zLaJn1QtJX9mWoqQ0oTihhBSEIzIjZbSuGuMX1BcV0tpmoamjuTIjX0Tv7jA0iM3PRzSSQgpLBQ3hJCCIAsHu4opWfCY0lIO+9y0dEVsz+uXSsFFxEYEj5iWIqRvQHFDCCkIslnXLjUlB03iaSl3kRs7vw2Q8NwIcROzj9x4s+iGTAgpHihuCCEFIRbLJG6M7YqiwCtSRQ4Fx24bvw1g9twYaan4NiMtRc8NIaUMxQ0hpCDIIsWu141qSkvJqaKeRW783mT/jtVQLHrdsIkfIaUJxQ0hpCBEM3pujK/lPjdpiqtMCHEjTVOAopi7EOuRm8RJvYn+N14P01KElDIUN4SQgmAyFNs08lMthmKPy1JwIW7qB1Xoj/m9HiiKTRO/xCmNyA3TUoSUMhQ3hJBeR9M0S7VUekOxiLiIY50getxMHFapPxZIGHeMJn7x6wpDsY/VUoT0CShuCCG9jjUgktlzo+jVUk5TRaLHzcRhA/THAr74vzyfpaxcXF4Yib00FBNS0lDcEEJ6nailWY2t50baxaNAMgE7u8buxETwicMNceP3mquhrJEbPS3lsjKLEFJcUNwQQnodayO+TKXg5iZ+TtNS8cjNhKFGWsrvNUduhKbSS8G91rQUxQ0hpQjFDSGk17FGbuwmg5v73BhVT06a+GmaphuKR9aUoyLgBSB5bvS0k33kxuuyMosQUlxQ3BBCeh2rl8XecxP/rCiiiZ9zH0xXJIZQQjANrgygutwPINlzIy6r97mxNPFz2g2ZEFJcUNwQQnodq0CxKwUXKSGPpTzbieCQPTx+r0cXNyItZURm4tcV/W6M2VLxY52mwAghxQXFDSGk10kSNzaGYn0iuCI+OzcUy6LE61FQpYsbczWUEDVif2spOA3FhJQmFDeEkF4nmiRukiM3QlgolpEITtJSMcvohqTIjcWcnKoUnIEbQkoTihtCSK/jyHOjirRU/Hs3Tfxi0rGKoiR5bpLGL6Ro4se0FCGlCcUNIaTXSU5L2Xlu4p+tnhtHkRuLh6ZGiBu9FDz+2Ri/YPH3sIkfISUNxQ0hpNexelnSzZbyKtYKJgfnt4gba1pKNOkzIjfm/b0uys4JIcUHxQ0hpNdxZigWnpv4924Eh1UYjUs08qurLgMgRW40e3HjpjKLEFJ8+Aq9AEJI/yMac+C5EWmpLASH7rlJHHvSwSMwuDKAQ+prAABeS+QmqRTc0geHEFJaUNwQQnodq0Cx99yk8sE4P79Xqn46etJQfbs3EbmJpSgF9zJyQ0hJw7QUIaTXcVIKntznRjzuJHIT/yxEihWvxZwcVa1Cyvm1CCHFB8UNIaTXiTmYCq5a+ty4GZxpTUtZEQMy9SZ+iWv5LIMzWS1FSGlCcUMI6XWsgZp0gzOFPlFcdA22GoqteC19bIQHyGMVUtQ2hJQkFDeEkF7HOhXcSZ8br4vxC1aDsJWk8Qua2XPDJn6ElDYUN4SQXiepQ3HayI37YZZGWsp+uyx6VFVLSmNxthQhpQ3FDSGk13HSoVjsIjJLRhM/52kpXwp1I4ubqKoZkR6rkKK4IaQkobghhPQ6bpr4WdNSTqIp8mwpO3xy5EbTkgzFbszLhJDig+KGENLrWEvB7Zr4aRZDsceF4FAdem7EWqyGYs6WIqS0obghhPQ6VoGSvs+NtUNx5vNbB2FakcVNTNVSG4qpbQgpSVx3KA6FQnj77bexZcsWdHZ2YtiwYZgxYwbGjx+fj/URQvogjpr4qebZUkKPZDMV3IpcIh6TPDcevUNxYg303BBSkjgWN2+88QZ+//vf4+mnn0YkEkF1dTXKy8uxd+9ehEIhTJgwARdccAF+9KMfYeDAgflcMyGkxEny3ETtPDfxz9beM5qbPjcpxI3Ho0BR4uXmMVVLGr/AtBQhpY2jtNQ3vvENnHbaaRg3bhyef/55tLW1oampCV988QU6Ozvx+eef46qrrsKyZcuw33774YUXXsj3ugkhJUxSKbhtWsqcWnLTxM/qobFDHsFgjdwwLUVIaeMocnPyySfjn//8J/x+v+32CRMmYMKECVi4cCHWrFmDnTt35nSRhJC+hZw2iqlaCs+NOS3lpolfpsiN2BYvA1f19XBwJiF9A0fi5sILL3R8wilTpmDKlClZL4gQ0vcR0ZcKvxdtoahDQ3HicUeem/jnVOMXALncWxJbrJYipE/AailCSK8j0kBBvxdA+j43XosPxtFUcIeRm/haVKO6ymMRUozcEFKS5EzcLFy4EF/5yldydTpCSB8mlgitlAfi/4Lsxi9Y+9x49WhK5vNn6nMjb1O1ZEMxm/gRUtq4LgVPxahRo+BJNciFEEIkRKCmXI/c2JWCxz8rWYxEsM6KssOnR25oKCakr5EzcXPDDTfk6lSEkD5OLKFcytKJG2uHYhcmXz0tlVrbSJEgm1JwDs4kpKRhqIUQ0uuISIkQN/ZTweOfrR2KnZh8HaWl7ErB9Z465vMQQkoL15Gbc845J+32e+65J+vFEEL6B0I0lKcxFGvWwZl6E7/M5880fgEAvF4pcmMZnOnGvEwIKT5cR2727dtn+ti1axdeeuklPPHEE2hubna9gMWLF2PcuHEoKyvDrFmz8M4776Tdv7m5GRdddBFGjBiBYDCI/fbbD88++6zr6xJCCocRuUkYimNqUudhETTJx/gFwD5y402KEmW8FCGkCHEdufnXv/6V9Jiqqvjxj3+MiRMnujrXI488gkWLFmHJkiWYNWsW7rzzTsybNw/r1q3D8OHDk/YPh8M44YQTMHz4cDz++OMYNWoUtmzZgpqaGrdPgxBSQGKWyA0QFzx+ySRj7VDsqhTcgaFY9txYxRCb+BFS2uTEUOzxeLBo0SIcd9xxuOyyyxwfd/vtt+P888/H2WefDQBYsmQJnnnmGdxzzz244oorkva/5557sHfvXrz55pt6t+Rx48alvUYoFEIoFNK/b21tdbw+Qkh+0MVNwBA3kZgKv9cIJuviJvGQK0OxJRJjRzpxw7QUIaVNzgzFGzZsQDQadbx/OBzGypUrMWfOHGMxHg/mzJmDFStW2B7z1FNPYfbs2bjoootQW1uLqVOn4oYbbkAsFkt5nRtvvBHV1dX6R319vfMnRQjJC0JMBH2SuLEMzxS6QvfcZDF+wZc2chP/9xeVxI21GzI7FBNSmriO3CxatMj0vaZp2LlzJ5555hksXLjQ8Xn27NmDWCyG2tpa0+O1tbVYu3at7TEbN27ESy+9hNNPPx3PPvss1q9fj5/85CeIRCK45pprbI+58sorTWtubW2lwCGkwBgdio33V9bhmcZsKfcjEcSp0qelEvtqhrgRhmKmpQgpbVyLm1WrVpm+93g8GDZsGG677baMlVQ9RVVVDB8+HHfddRe8Xi8OO+wwbN++HbfeemtKcRMMBhEMBvO6LkKIO4SY8Hs8CHg9CMfUpF43Rim4+bMTwaGPbkibloqrm1hMk/riWNJSNBQTUpK4FjfLly/PyYWHDh0Kr9eLxsZG0+ONjY2oq6uzPWbEiBHw+/3weo1Q9oEHHoiGhgaEw2EEAoGcrI0Qkl9kw6/fqyAcS27kZzUUe3NtKE5siqqqngLzsokfIX2CgjXxCwQCOOyww7Bs2TL9MVVVsWzZMsyePdv2mKOPPhrr16+HKr2d+uyzzzBixAgKG0JKiJjkifH74v+GrOLGOltKcdHEzzAIp97H5xFl6Mb59GopNvEjpKTJmbj55S9/6TottWjRIvz1r3/F/fffj08//RQ//vGP0dHRoVdPnXnmmbjyyiv1/X/84x9j7969+PnPf47PPvsMzzzzDG644QZcdNFFuXoahJBeIBYzqpNEhVTYYig2+ty4NxS7qZaSuyNbIzf03BBSmuRsttT27duxbds2V8ecdtpp2L17N66++mo0NDRg+vTpWLp0qW4y3rp1q2kYZ319PZ577jlceumlmDZtGkaNGoWf//znuPzyy3P1NAghvUBUKr0OeO0jN0YFE0yfnURT9A7FDvrcpBM3No2TCSElQM7Ezf3335/VcRdffDEuvvhi220vv/xy0mOzZ8/GW2+9ldW1CCHFgVyqHciYlnLfe0Z1FbmJJT2m+3scpqUeX/kFVm/bh+u+MTWtoCKE9A4cnEkI6XXkQZWiK3FyKTj0fQBpXIKbtJSTyI10Xa+lz43TtNSdL36GB9/aijU72SSUkGIgq8hNR0cHXnnlFWzduhXhcNi07Wc/+1lOFkYI6bvEEkUBPq/hubEOz1QtqSWPC5OvXtqdbVrKRU8dAOiOxM/RHnLeyJQQkj+y6nNz0kknobOzEx0dHRg8eDD27NmDiooKDB8+nOKGEJIRObJiGIoz9bnJIi3lYHBmKHFdj2JnXnYmboRY64qk7pZOCOk9XKelLr30UsyfPx/79u1DeXk53nrrLWzZsgWHHXYYfve73+VjjYSQPoZczZTKUJzkuXFTCm451g6v1xy5kYWQx0VlFmCk2UIUN4QUBa7FzerVq/E///M/8Hg88Hq9CIVCqK+vxy233IJf/vKX+VgjIaSPIVdL+X1xIZGqiZ/QJ0J8OAmmiFO5idyYxI3LtJTYj5EbQooD1+LG7/fr5dnDhw/H1q1bAQDV1dWuS8EJIf0TeZZT5rRUXGgIkeOka7CTtJQYqilElVxZpTfxc5iWEmKtK8x5DYQUA649NzNmzMC7776LyZMn49hjj8XVV1+NPXv24IEHHsDUqVPzsUZCSB8jZqqWymAoFpEbFz4YuRorFR6Pk7QUIzeElCKuIzc33HADRowYAQC4/vrrMWjQIPz4xz/G7t27cdddd+V8gYSQvocQH77E4EzAznMT/5zU58ZBqkgfnJl2/IK5FNxO3DhJS2nSVPFuihtCigLXkZuZM2fqXw8fPhxLly7N6YIIIX0f1VQtlcJzowrPjXUkQubzxxxEbqyl4LK4MZr4Zb5WVFpQV5jihpBigE38CCG9TtSuFDxlEz/o+wLuqqWc9LnpTAiSoM+rb3PTxE9eDyM3hBQHjsTN1772NUcjD9ra2nDzzTdj8eLFPV4YIaTvohuK5angSYMzraXg5sfToUrnT4UQN23dEQBAeUAWN6IbsnN/D5Dec/Pou9vw4wdXUgAR0gs4Skt997vfxbe//W1UV1dj/vz5mDlzJkaOHImysjLs27cPa9asweuvv45nn30WJ598Mm699dZ8r5sQUsLITfwy97lB4rNzk6+elnJQCi66Cpf7DXHjZrZULOZM3Pz5lQ3YtKcDZ85uxuyJQzKelxCSPY7Ezbnnnosf/vCHeOyxx/DII4/grrvuQktLC4B4PnzKlCmYN28e3n33XRx44IF5XTAhpPSJOfHcJDSD1XMTc+CD0Q3FDpr4tXUnxI1N5MaJvycqGXPSRWV2t4WS9ieE5AfHhuJgMIgf/vCH+OEPfwgAaGlpQVdXF4YMGQK/35+3BRJC+h6yJya158acljKa+OU4cpMQNxWyuHHj73FgKO6OxPQIUdSJYiKE9IisBmcC8aZ91dXVuVwLIaSfYEpL+ezTUsmzpRLHOhE3iV3SRW6EH6c9nCxusumpAxgDNK00dRgDhmNOxpoTQnoEq6UIIb2OSM3IkRuroVhLmgruwgejJpd3W/HokaD49+V+472eG/Ny1IHnpqk9ZKzNYWNAQkj2UNwQQnodYTvxSYZia1oqpve5iX+fVZ8bB+MXBNmmpZx4bprapcgN01KE5B2KG0JIryMEQXz8grlTsMA6W8rromuw0BtpDcUe878/U1rK1ZDOzJGbPVLkhp4bQvIPxQ0hpNcxDc7U+9ykMhTHv/e4GGYZczB+wbrNXC1lPk86nHQo3mOK3LBaipB841rcLFy4EK+++mo+1kII6SeYmvhl6HPjTRq/4LyCyRqdkUkXuXEzW8pJ5MbkuaG2ISTvuBY3LS0tmDNnDiZPnowbbrgB27dvz8e6CCF9GHlqdyDlVPD4Z8VSCu4kq+NkcKbXkrEqDxiGYjdpKTlyE3JSLcXIDSF5x7W4efLJJ7F9+3b8+Mc/xiOPPIJx48bhxBNPxOOPP45IJJKPNRJC+hhG5MbjuM+NsM9YoymNrd34zVOfYP2utqTzpx2caVE+codid5EbY93hmIqoTWimNzw3nze2oSPRS4eQ/k5Wnpthw4Zh0aJF+OCDD/D2229j0qRJOOOMMzBy5Ehceuml+Pzzz3O9TkJIH0JPG3kzdyjWB2dKQkVu5PevVdtx35ubcf+bW5LP76CJn8A2LeWyFBwAuqPJ4kaulnJSyu6WT3e24oQ7XsWlj6zO+bkJKUV6ZCjeuXMnXnjhBbzwwgvwer046aST8NFHH2HKlCm44447crVGQkgfQxcfimQoTjVbymP23MjHA8ZU707JzCufPxXWUvBy22opd4ZiwN5U3NSR38jN57vaAQBb93bm/NyElCKuxU0kEsE///lPfP3rX8fYsWPx2GOP4ZJLLsGOHTtw//3348UXX8Sjjz6K6667Lh/rJYT0AaJSZCWQoomfSEvpfW4kMSLrA5EWksVRzCKM7LBuq/DbVEs56nNjidxYTMWapuW9z01LZ/z81tQeIf0V1+MXRowYAVVVsWDBArzzzjuYPn160j7HH388ampqcrA8QkhfQ07LpKuWsva5kbWIXDElxEVYSgepDtJSyU38pA7FrmZLmddtFTctXRGTAMqHuGnujPsdUxmaCelvuBY3d9xxB7773e+irKws5T41NTXYtGlTjxZGCOmbyDd6jyddEz9znxtZqMgCQXhebCM3adJS1shNuc1sKUfVUhbPjbUcXO5xA+QnLdXcFRc3jNwQEsd1Wmr58uW2VVEdHR0455xzcrIoQkjfJeYwcqMlRW7ktFRyJES+sesdirMdv+DCUGyNxFg9N3KPG7v9c4GI3IRtzMyE9Edci5v7778fXV1dSY93dXXh73//e04WRQgpbuzKnZ0iCwbzVPBUnhsbcSNdPmrnuZGaBKbCKnzMs6XM50lHkqHYErmRe9w4PadbWhKRm1DUvokgIf0Nx+KmtbUVLS0t0DQNbW1taG1t1T/27duHZ599FsOHD8/nWgkhRcCSVzbgkGufxyc7WrI6PhYzixtjKniGUnCPfeTGSEtJ0RwHaSlrJVW5TeTG7WwpINlz0xuRm5auhKE4qjqq8CKkr+PYc1NTUwNFUaAoCvbbb7+k7Yqi4Nprr83p4gghxcdbG5vQEY5h1dZmHDSy2vXxpsiNYkRuQhma+MmBlpiNoTgSc2co9kotij0K9Kot+TgnaSlrOq0gnptEWkrV4uf3W9svE9LPcCxuli9fDk3T8JWvfAX//Oc/MXjwYH1bIBDA2LFjMXLkyLwskhBSPIjIQ2c4u264Io2kKGZDcSQWjzooetTEbCiOv7mKR1NsPTfRZENx+vELhgCoCPj068avmd1sKQDotlQsyT1u4vvn3hcjDMVA/OfgT/fECekHOBY3xx57LABg06ZNGDNmjOkfASGk/yBu5h2h7PwdVj+MiJhoWnybLyF2hAawio6Yppk8NyJyYvLcxDKnpWQ/jpySih9nfC0LLjsyNfFr7YqLwKDPg1BUzfngTE3T0NJpiJtQVEVlMLfXIKTUcCRuPvzwQ0ydOhUejwctLS346KOPUu47bdq0nC2OEFJ8iJt5qgnYmbCORpCjDJGYBl9CZ1jTUkA82hKDZkoXxfS0VLLnJl1aSi4Fr7CIG2vZuS9NmiepWsrycxHbDXGTW3XTFYmZKsVYMUWIQ3Ezffp0NDQ0YPjw4Zg+fToURbE1rSmKgliMbn1C+jJG5Ca7tJR1NIIsbsIxFeWICw1rnxvA6FYsNwKM2PW5cTA40xS58VsiN7K40bS0/ygzdSjWxY3fC3RHc+65ae40t+ZgxRQhDsXNpk2bMGzYMP1rQkj/RY/c2MxQcnO8EbkxhITJFGzpcyMfY/bcJKel1B5GbuRrZvIUWyMx1p+LeL4i/abmuJrJKm4YuSHEobgZO3as7deEkP6HuJl3ZGkoFlEXX+JmryhxU3EkptkKFDn4IqI9cvDDbvyCk6ngPpO4Mf8rlM3GmUzFmfrciJ9XMFEVZu1o3FNauqyRG4obQrJq4vfMM8/o31922WWoqanBUUcdhS1btuR0cYSQ4kPoj84eRm7k6IjfZnim0AyyQBGH2I9fiH/WNM32WCvy9ZMMxdJ/xkzl4FaxYq2WEptFyXuu+9yIHjcCihtCshA3N9xwA8rLywEAK1aswB//+EfccsstGDp0KC699NKcL5AQUlyISES24saue7AQN7IxVrMzFHvMZeLy+UTUR9YO1kZ9MrJJOG1aKoNWyOy5MUdunPTOcQM9N4Qk43pw5rZt2zBp0iQAwJNPPonvfOc7uOCCC3D00UfjuOOOy/X6CCFFRjRXhmIbcZMpLWU380n0zYmqGlRVM0VGrMMxZdIZik1pqQxiRIiXAUEf2kPRlNVSInKTc0NxFz03hFhxHbkZMGAAmpqaAADPP/88TjjhBABAWVmZ7cwpQkjfQtyssy0FtxqKASAgNfITCJ+uXeTGlJaSK6dU1WTYzT4tZT/qwQ5x/QHB+HtFq6HYKm5iOfbcJEduKG4IcR25OeGEE3DeeedhxowZ+Oyzz3DSSScBAD755BOMGzcu1+sjhBQZwmOSqyZ+AOD3pY7cZBQ3kliwDt9Mm5aSjDXWtFT8uvEUl5oh0iLEyoAyH9CaLPqE+AkmGvjkOi1l9dwwckNIFpGbxYsXY/bs2di9ezf++c9/YsiQIQCAlStXYsGCBTlfICGkuNAjN1lWS6VLS4UlQ7HQAB67tJQkOOSvI1HVkpZKvQ55m7VaSl5fRkOxJXJj9dwIcSRKwXNtKGbkhpBkXEduampq8Mc//jHpcQ7NJKR/IG72nZFYxtEEtse79tyk73MTkfrMRGIq/NJ7NqeRG6vnBhBCSkMmLSKez8Aye3ETzbPnxloKzsgNIVmIGwBobm7GO++8g127dkGV/rEoioIzzjgjZ4sjhBQf4mauafGyZ6tfJePxNg32bD03Nh2KjbRU8nqAeLWV7JdJ38TP+No+LSVmXGWK3MQXU5mI/qQbvxD/PrfiQ0RuBpb50NYdZbUUIchC3Dz99NM4/fTT0d7ejqqqKtO7KoobQvo+UUlZdISj7sVN4uZuWwoelcVN/LNbz03Aa1RZpYsqmSI3Ns/B7lp2RGXPDRwYivMUuamtKkNbdzsjN4QgC8/N//zP/+Ccc85Be3s7mpubsW/fPv1j7969+VgjIaSIkG/OnVmYioUYkCMs4sZv2+dG+i/ltfHcRC1pKREZ8qWJ2sjnAuw9N+LwTNVSyWkp1XZ7OnGzuy2EpR/vzCrq0twZNxQPHxgfBU7PDSFZRG62b9+On/3sZ6ioqMjHegghRY7sGemMuDcVqzbiw/DcJHcolqMvHhuTryktJRmK0w3NBABvmiZ+8rXcloKHYypUVdOPt1ZLyT+/UDSGO174HPe/uRldkRh+M38Kzjp6fNrryYSjKjoSkaLaqjL9MUL6O64jN/PmzcN7772Xj7UQQkoA+WafTTl42vELGUvBE9tS9bmJqXp/nHR+G8AcubFNSynJ/h47rJ4a65qs2+W1P/vRTix5ZYPu09nTbi7rzoRsJh46IACAHYoJAbKI3Jx88sn4xS9+gTVr1uDggw+G3+83bf/GN76Rs8URQooP+cadzWRwvc+NV05L2RmK459NhmK7tJTFcyPSVOkqpQCz+LGtlnIcuRHjFYxzqDaRJbtqqdYuc+Qr4tJsLH7+5X6v/hwYuSEkC3Fz/vnnAwCuu+66pG2KoiAW47sGQvoqqqpBvtdnMxncKAU3Ih12hmK72VJ2aSmr50aP+GSI3Jingts38ZPXmwqreImvKTmyFLTx3FjLwuXBoU7oTkRpygNe/fr03BCShbhRc1zGSAgpHaw342wiN/r4BUl72HtukmdLCUGipojchGOqnkbKlJbyeBRUBLzojsRQXe5P2i4iP049N3JaKmbzPOwMxVFLziuSKQdmQfz8y3wew5RNcUNIdn1uBN3d3SgrK8vVWgghRY41ipHryI18c7czBouvhaDQNM3suXFhKAaAP3x/BtpDUdRUBJK2GWmp9OcQ4kqO3JgiS4nnZNehOCly41bcJLw6ZQGvnhYLuTwHIX0R14biWCyG3/72txg1ahQGDBiAjRs3AgB+/etf429/+1vOF0gIKR6ilshtNqXgdrOl7Jr4GeMXUncotgqPSEzTt3kd/HebM6UWp84YZbvNbtSDHeJn4vd69CiT/HMShwf9yZ4bq5gJuxQmohtyuV9KS0UobghxLW6uv/563HfffbjlllsQCBjvdqZOnYq77747p4sjhBQX1qx0Zw8MxbazpRx3KI5vs4qDSMyI3GQyFGfCbtSDHeJ6fq+iCza7PjwBrzd5m2XQp/X7TAhxU+b36mkxtwKJkL6Ia3Hz97//HXfddRdOP/10eL2GCe+QQw7B2rVrc7o4QkhxkRS5ySItFbUTN2IquGSo1aulPMlpKSEQrFGVsNTET+5jkw16E7+MkRsjzWbX1ThdEz9RHRWwmYruBNEw0By5YVEHIa7Fzfbt2zFp0qSkx1VVRSQSsTkiM4sXL8a4ceNQVlaGWbNm4Z133nF03MMPPwxFUXDqqadmdV1CiDusYiKbyI1qk5Zy2ufGZ4mm2HlW1BxFbnQh5TBy4/MotqXq1j43UZvIjajWytpz4/fonhtGbgjJQtxMmTIFr732WtLjjz/+OGbMmOF6AY888ggWLVqEa665Bu+//z4OOeQQzJs3D7t27Up73ObNm/G///u/OOaYY1xfkxCSHVYxkY2hWG/i59hzYxyrl4IndkuqNpINxRmqpTKhp6UyaAUhULweJSlyo6rGVPGAzeBMsf4K0aOmB2kpem4IMXBdLXX11Vdj4cKF2L59O1RVxRNPPIF169bh73//O/7zn/+4XsDtt9+O888/H2effTYAYMmSJXjmmWdwzz334IorrrA9JhaL4fTTT8e1116L1157Dc3Nza6vSwhxjzVyk10TvzSDMzN1KNYjI6rteiIxLWeeG8VhKbgcufFZKqLkqI9tKXji67JE5MYq1jLRJRmK6bkhxMB15OaUU07B008/jRdffBGVlZW4+uqr8emnn+Lpp5/GCSec4Opc4XAYK1euxJw5c4wFeTyYM2cOVqxYkfK46667DsOHD8e5556b8RqhUAitra2mD0JIdiRHbrIRN/HPtobiaLK4UdIYiq3rMXluehy5Saw3g7gRvhmvR0kqVZeFjG0Tv0SkRnQXdu25CdtEbjh+gZDs+twcc8wxeOGFF3p88T179iAWi6G2ttb0eG1tbUpz8uuvv46//e1vWL16taNr3Hjjjbj22mt7ulRCCGw8N6Fs+twYYkDgtzHUGuMX7DoUx7+3VhfJ1VJO+tykQ2/i57BDsc+mWspO3JhKwRM/C+G5cZ2WSojB8oAUuWETP0LcR24mTJiApqampMebm5sxYcKEnCwqFW1tbTjjjDPw17/+FUOHDnV0zJVXXomWlhb9Y9u2bXldIyF9mVwYiu0iK4bnxji/3fgFa4dia/WWPH6hp5EbIy2Vfj8hsHw21VKykBGl4KpmE7kJxN9nRlwKE7lDcZDjFwjRcR252bx5s+38qFAohO3bt7s619ChQ+H1etHY2Gh6vLGxEXV1dUn7b9iwAZs3b8b8+fP1x8Q4CJ/Ph3Xr1mHixImmY4LBIILBoKt1EULsyWUpeOZqqfhnk6HYUsFk77lJ7JsjQ7HT2VKyoVg8RznqY9fET/w8K7JNS9l0KGbkhhAX4uapp57Sv37uuedQXV2tfx+LxbBs2TKMGzfO1cUDgQAOO+wwLFu2TC/nVlUVy5Ytw8UXX5y0/wEHHICPPvrI9NhVV12FtrY2/P73v0d9fb2r6xNC3JGTyE3MplrKJqVieG7kDsXmdUQsaZxwVG7i53ppJtzOlpLTUnal6kLAaVpc9Hg8ir5+kZayeogy0WXXoZjihhDn4kaID0VRsHDhQtM2v9+PcePG4bbbbnO9gEWLFmHhwoWYOXMmjjjiCNx5553o6OjQq6fOPPNMjBo1CjfeeCPKysowdepU0/E1NTUAkPQ4IST3WG++PUlLZYzc6N4Z41hrNCU5cmOkpXwe11l3Ez6b8nQ75OovIdhEusnUA0d6IjFNgweKXh0lqqXcRl1EEz+5Q3FM1RCNqXrlFiH9EcfiRqR/xo8fj3fffdex5yUTp512Gnbv3o2rr74aDQ0NmD59OpYuXaqbjLdu3QpPD/9JEUJyg7hZVwa86AjH0BmOQtM0U3TF6TnkUu2AzVRwu9lS1g7Fdp4bo8+N4yXZIntYVFXDzx5ehfFDK/E/c/c37Sd3KE4yFGtGlMokblQNfq9xbE/TUnLkBohXjVHckP6Ma8/Npk2bcr6Iiy++2DYNBQAvv/xy2mPvu+++nK+HEGKPuGlXlfvREY5B1eI3/zK/N8ORBrIYEOSuQ7Fm8sD0BNnDsm1fJ/7z4U74PAoWnbCfScwZhmKpiZ/wBEnb5EiVWHcuOxQHJDETjqqwGXROSL8hq1LwZcuWYdmyZdi1a5ce0RHcc889OVkYISQ/PPPhTnSEovje4e49akI4DAga/zo6QlFX4kaVPCoCf+LrsI2hWA4KeazVSFbPTQ5LwYUBOBRVdRERVTW0dkVRXeHX97MzFOtNBrXkbab1q+a0VE8GZ/q88WqtmKrRd0P6Pa7FzbXXXovrrrsOM2fOxIgRI1yFowkhhUXTNCx6dDXCMRXzDqoz3aSdICIOfq8HZX4PuiMqOsMxDMniHLL4sO9zk2w89maqlormrolfUGqKJ480aOoImX5uQqD4vLK4QeKz0dNHTsNZDdHl/uzmQnVJ4gaIp/e61Bgrpki/x7W4WbJkCe677z6cccYZ+VgPISSPRGLGu/q2UMS1uIlJN/KKgA/dkbBrU7FsshXonhtpKrjdbCldPCREQcSuz02Oxi+ItFQoopoiIXs7wpgwLP61PDvKZ/LcxPeXy949HgWKEn9eUX27uYmf27RUSJoKDsSjTV2RGLsUk36Pa8dZOBzGUUcdlY+1EELyjHzz7M5iwKI8JFLcUEX0wCl2nhg7z03MxnPjSeFpEURimsnE2xNkQ3G39BybOsJJawTSj18Qz9VqOLY28VO1zH11ZPRS8IARuRFrJqQ/41rcnHfeeXjooYfysRZCSJ6RxUM27+5VqYxbeFK6cyJu7Dw3NrOlLCMRrNVS4VxGbvxSWsoSubE+F0AMzrQvVRdrsVZ7idejQvIsuYneGB2KjchNfM0UN6R/4zot1d3djbvuugsvvvgipk2bBr/fHNa+/fbbc7Y4QkhuCfc0ciP5ZcQN1a24idrNlrJEbjRNsy0Ft1YjJVdLqTmvlgpFVZMQlMWNfP24adg8HFOvDPMakZuQzXaRlhLPwYlBW9M0dEdFh+L4dQM2A0gJ6Y+4Fjcffvghpk+fDgD4+OOPTdtoLiakuJH7yGQTuZGHRJbpkRt3N1Khr0yeG5+5z43cFNhrJ250w65Nh2LNvG+26GmpiGo2FLdLkZuYOXIjCsCs4xdEQ0HreAaRlioziRtnaalQVNV/TkIMyYKMkP6Ma3GzfPnyfKyDENILyIMZQz3y3HhQlgjauhVJwmzrsYncxNR4nxpNUje2kZvEOYQQCHg9CMdUs6E4h9VS3abITUj/WjY0y5Eb1RKZEUsRa7Km1QJSGbfTtJT8+gn/k90YC0L6I2xhSUg/ItxDz41c6SSiBW5Fkv3gTOPr+AgFY39F+i9leFbEehLVQnq1kWZbap4NQb+UljKVgid7brweBYqi6LOvrIZiI3JjHp4pxKLf6zF8Rw6FiTAT+zyKLg5lQUZIf8Zx5OZb3/qWo/2eeOKJrBdDCMkv8o2zJ54br9Rxt9vljVS16UPjl7vrxlRTt12PzeBMa4ficr8XLV0RhKOqdH5Xy0pCrpZKZSi2CjWfxXOTqVpKRGm8HgV+jwfdUB0Pz+y29LgBGLkhROBY3MhTwAkhpYm5FLxnQy/FTdW1oTiWXtxEoir8HlncGMcmGXZtxhfk3FBs6RtjqpaKmcVNqsGe3hTbjaaISryRYch5tZS1gZ9pzRQ3pJ/jWNzce++9+VwHIaQXMBuK3d8AYzHDL5OtoVguJxeI8QQxVUM4pqIslecm8aW1WqrMnyxuepyWMvW5MaelxLBQa+WXU3FjTUv5epCWEq+DvGZGbkh/h54bQvoRPY3cyGkYI0rgthTcMCXLiJt7NKbpAgiw9LmxGHLtPDd2aa9skHvGyM8xHFXRkegvY1SPpaiGypSWSqzfL/lmnEZu5Ing+popbggBQHFDSL/CbCjOInIj3ayNtJTbUnB7T4xIRVkNxXYdiqO6Z8Wclsrl4Eyj228s6We1N1EOnkq8qFrmtFS8Kix+Pp/XmOqdC88NDcWkv0NxQ0g/IhLNZeSmpx2KLZEbqddNylJwxRq5sUlL5SxyI82Wsgi4pkQ5eNTiudHFV+LxmCUFJ4sbOULj8xrdjSMORWe3Za4UQEMxIQKKG0L6ET323EjCpMeRG0tkRQiA5MiN8XWqDsXiBq9pxo3d11NxI3tuLJEQYSq2em58lvXFUnhyoqq5Ksrv8ehpKaeTwfXRC4HktBQNxaS/Q3FDSD+ix9VSpj432aVArKkcgew5MXtu7Jr4CUOueao2YIit3A3OjNlEbsKmdfgtnht9KnjMPm0VUzVzd2Ov7LlxmJYSoxd8xr/xAMUNIQAobgjpVxSD50YfSeC1ihvDT6PqE8HNx3pTeFpk34kQbT0fnJk8W0pc34jcWDw1lqng1sowj5yWkrob+zyK4blxGbkpD7AUnBArFDeE9CNyVS3lzUHkxmr4lSM3dkMz5e9FRMRaCg4YzytXkZuw1KG4dmAQgCFu5EgWYAzITB6/kBy5kaM6imJ4bpympXRDsc8mLZXFa0tIX4LihpB+RKSHHYpFukUuBc/WUGz1xPiktIwRuTHvY43ciCiH32tEPnIWuZHSPa3dEQDAyJpyAMbwzFTVUknjF7zJfW6E0BTHuE5LCUOxFLnRo2msliL9HIobQvoRPZ0Kbhe5yb4U3Cw+AnqfG8NQbNUnxmyp5D4yIq0l1tPz8QuGaGjrjgIARiTEjRieKYs9wBBUyU38zJ4cVdOk7sQe0+eedCjO9jUhpK9BcUNIP8LkucniBqjKhuIeN/FLFbkxJntbIzdGNVLiXPLgSVGaHs1NWsrvVXRx1doVj9yItFRL4vuIxTBsHQ9hVIbBtD0a0/Sok4jqBHxGtZgTum06FGc7EoOQvgbFDSH9iIjJUNyTyI1HN9y6j9zYl2obpeAODMW2kRtz1KKnaSlFMXr5tIXikZuBZX59jfHnYp36DdPj1m7MsucmEjMf6zYt1WXboZjihhCA4oaQfkVPp4LL3YV73sTPkpbySZGbVIbipNlNhlDSOwpHzJVNPUFOTQHAgDKfvkbAplpKRGaS0lIw7RfTNL1Hjkin+Tzu0lJ2HYqZliIkDsUNIf2IXEZusk2BpBI3uhlXjtxY9rF6WiJSmszw3ORS3Jj/RQ5MiBuR3tPFlS5Q7Dso654bqVRc7/ZsTUv1oEMxDcWExKG4IaQPsG1vJ95cvyfjfnLKoyeRG7mJX7fLniqpxiPoaRlV1ccvWPWJmNigdwAWnhev0eG3K5fixm8RN0FL5CbV+AVrtZSlVDwWUw2/UFJaKvsOxWXSyAhC+jMUN4T0AS7+xyr84O63sX5XW9r9wjmL3BhN/MJR1TQLKhNWn4pAv7lHU6elxDFWz43PxnPT08GZQJq0VNQ+MmOd+m2NUsmGaKuh2BB32XcoNtJSjNyQ/g3FDSF9gMaWbgDAZ43taffraZ8bubuw7PVw0xFXb2xnnQpu06FYSepzE/9szJYyPDfWaqn8pKWEodjsubEbjClvt3YwjqmqLmJ0YeRycKZdh+IyGooJAUBxQ0ifQNzQdzR3pd3P1KE4GnMVcQEMMeGRKokA5zdTVdX07sPWyI1PGhwpJhMkpaVS9JHxeRW9T444f0+rpYBkcTMgaPXciLRYcpM+IHn8grxdbkAIQDdEO01LCUFZZuu5YVqK9G8obgjpAwiPxfaM4sYQM5rmvOxYELOkgcTN2mkUKCaJqWTPjY2hOFWHYn1wphH98Fu69vW0zw1gTkv5vYbPSAgQa5fh5A7KZmO0T/fcyKXg5rRU2GkpeDi5FFysLyZ1QCakP0JxQ0iJo2ma7p/JFLmxzi1yW1VjTbOUuSwHj6npxI3dbCnz8R7LYEq9nFry3Bjnd7SktMiG4qDPq0dXREl9Kk+NYSg2ix898iSVgvuy7FCse25smvgBTE2R/g3FDSElTrzpXfzrHc3dGfY13zjTVdWoqoYfP7gSNy9dqz9mnZVU5nc3hVoWN8lN/JJnS1k9N+K6emTEpomfIDeGYo/pa3ENVUsMvxQjFCzjFYSosVaG2Q3O9OuGYmP8hBOEeAnaDM6Mb2fkhvRfKG4IKXHkqic3nhsg/bv7zU0d+O/HDbjvjc36Y3IaCHDf7j8qiZukqeA+ebaUvek41ewmv9eDyqC5sinXTfzK/F7dtAzEf5apPDfJ4xeSm/wZZmT3HYrj0br4aylHl+SuyozckP4MxQ0hJY4cNWnqCKe9qYkSZrtjrbQnRg7IqayYxSDr9kaaLnLjlzr0Ou1QLM92OqCuyrSv9fzZkBy5Mc4ZjqnJ1VLWwZ4xq/iJH6vaGIr9kqE6E2EpdSenouTvsyn1J6SvQHFDSIljFRY7W1KnppI8N2lESXtiEnZM1ZIiEUJ0BF1W54jjFSXZ8Cv3edFSGYoVkZYS5zM8LdNGV5v2zUlaSoqKBHweXYAB8ZJtcX2vxVCse26SqqWMyE0kKXLjfHCmLEqtFV0cwUAIxQ0hJY81+pIuNZXkuXEQuZGPs0Yq3DaNs6ZpZOQ+L0K8WHczxINlPV4Ppo6stt23J5j8LH4vPB7FNODT+vPQPUEpxKDZc5MQRvr4BWNieCbEz1tRjBJyASeDE0JxQ0jJY72JpSsHTzYUp4ncSOLGOktJ3JBF0zinhuKoJdIhI27SchO/1KXgifNJaanqCj/GDK7Q981FKXhA7v6b+Foe8Gn1ICVXc6Vu8meMXzAPznSSlhJG8KDPk2S6Nhr5uY/caJqGu1/biOXrdrk+lpBiguKGkBLHXeTGHBVIVwreIYubqP0sJbeRGyFK7PwwIuoRlg3FKSI3+mwpi3g4WEpN5bqJn0jByd4Y6/WFQBGPiwiOXdoqklQK7iYtlTwRXCBek2w8Nx9vb8X/e+ZT/OKxD1wfS0gxQXFDSIljLedOJ26ESBHRh3Sl4O0h4+YobriqpbRZpG3SRYBk9A7HtuJGpGXkPjcWQ7HVsGuZyn3wKEnc5Dot5TNXNYWjatL1U49fSC4VTyoFlyJCmRBRmTJfsrgJ9iBy89H2FgDAnvYwmjvDro8npFiguCGkxLFGX9L1uhEpj6rEAMh0kZv2UMQ4Lmr13IhScHfmVWukQyYgdSg2jMf2aSlAVByZ1zNNEjf56HMjrzNiE7mxihs1KbIjtkuDMz3ivM49NyIqY51aLj+Wjedmzc4W/etNezpcH09IsUBxQ0iJI6Iv4l7uxFAsZiSli9x02ERujI688cfdmleNpnbJ/3pkz0nKtJQkWOJdfs2RpIMkcdPWHUFPCdp0//Wn8dxY02bGkFCr+JHKyLMoBU8XuTHmS2Uhbna06l9vbqK4IaULxQ0hJY54Fz96UDmAuKE41UBMMXFaTLdOWwpu8twkz3IC3HcoNo5P3uaXqoVSpqWk4+SKI5HaqS7369tHD6pAT0mfltJsPDdG9EmsUX7cNDhTakAISNViLjw3dpEbQ3C6S0vFVA1rG9r07zftprghpYuv0AsghPQMEX0ZXVOBbXu7EIqqaA9FdQEjIwzFeuQmXSl4t121lKWJn1tDsWZOI8mIqqGIqUOx/YgGsRZr5AYAXv3F8diytwNTRpqb+mWDOS1lNhRHpCZ+VsNwUodim/EL1qGb+lTwaOa0VNrITZYdirc0daAzbByzqanT1fGEFBMUN4SUOOJdfE2FH4oSn/bdFY4liRtN03SRMlB4btKlpcKp+9xYDcVOUyB2YkQgN/EzOhSb9zFFbjQ5cmJsGDOkAmOG9DxqA1irpZI9NykNxZq9uJE7LOt+oSwGZzqJ3Dg1eQs+3dlm+n4zPTekhGFaipASR38X7/eiMhAXLfI7cIE812lAmYjcpL4BtkmRG5HOUi0dd7M1FNuJG3MTv/QdioGEodjiW8k1QanUWkRJ7PvcJEdmgOS+PqYmfpaJ4W5Kwbv1PjepS8Gddo0WCDPxzLGDAMTFTar0JiHFDsUNISWO/i7e50F5IH6zk6MuAvmmWaV7btIZio1zhETkJma+WRsN41x2KE7bxE+Vxi+Y95GPC0s371zMkbLDLnJjGH+TOxR7dM+NmApuXp9p/ELMLMzcDM4UUZmytJ4bd5EbYSb+2tQ6KArQFoqiqYPl4KQ0obghpMQRvpl45CZ+Y+uyidzIXg7huXHaxE9EbqzixK2hOF0puE8qhTbGL5j3UxRFrwqTr+mzcyjnALtScF2ERFVdQIifg3he1tlXRrUUEtslM3TSVHAHkZtousiNIW52t4Xw4RfNGc8HAGt2xsXNjDGDMLI6bk5naoqUKhQ3hJQ44gYbj9zERUuHjbiRS4wrHZSCt9mMX7B6Ztx2KNZLo9PMlkpXCg4YAiLUK5EbKS1l06E4ZBEZxvgF+47OeuQmJg3O9JrTUlFpcGgqutNEboxJ7Sp+9OBKnLL4Dazf1Z72fJ3hKBpbQwCAybUDMH5oJQBgI8UNKVEobggpcYwbrEeP3HSGUqelAl6P5MuwFyWappkjN0nVUp7ENYV51VnkRk3jkZGb2OlN/JC8nxAQsl8oF92I7ZANu3oTP5/hjbGKDPG89CZ+WmpPTsxqKJaiRJlSU1ZRJSNHbtbvaoemIWP0piExSb4y4EVVmR/jhsYN2YzckFKF4oaQEkcfouj3oiKY2lCsixufR6qosRclXZEYJP+xntJKGblxWS2VLnITiam6n8auGkhcu9c9N5ZS8LBNWirl+AUluZpKRHdECbw83TtTaipd5EaspSMcRUtXvJHhxgw9axpa4+KmtroMADBuSDxyw0Z+pFShuCGkxOmWDMUViRtbZxpDsd+rZByu2G6J/IRiqh51AeRqKbeGYnOFkIzsORHnK7cZDCmEghA3Xo+S5M3JFema+EVimqlSTV6bqsWjVDFLpMorTQ2PWCI38s8kk7hJH7mJn2+nNIZjw+70aSkRuRmREDfDBgYBAPs6et7lmZBCQHFDSIljjtwIcWPjuYkaHXEzpZPk0QtA3Dwrl5J7vebIjXNDceJ4O3HjSS0aZDwWz02+UlJAivELkggzpnMLgSI1GZT68IhIlVdPWxk9cvxSjxyh0TKNYEgbuUm8ttulMRwZxY2I3FTFxY3RUiBZJBNSClDcEFLiiMhNmc+DCr0UPHVayu/AcyN3JxbHxmRxo1ia+LmcCm4rbnyGGbfLku6RsRqK/fkUNzal4EGfHGEydwr2Sl6imJrcZFAezxCxDP1UFEUXTpmGZ8o+Kyt2FWyb93SaXj8r1siNEMl2v0eElAIUN4SUOHLkRrzj7kqTlgr4MkdurGmpcNSINAB21VI5aOInRW660kQmPJ7ktFS+kH0wRlrKqOpK8two9uJGH8+gp620pJ5B8vUypqXSiD+7n1k4puKLfanHKQhxU2eN3NgY0wkpBShuCClx5CZ+FelKwaPJnpuUkZtQ+siNPlsqyyZ+dp4bWUiIyFE6z4143vnqcQPEoylC1FgNxd3hmJ6qEz9PWajIaSnr7Cnz4EzjGKfDM+XeRlbsHgPSp6YaLWkpEQHsdNkIkJBigeKGkBJH9qdUpCkFD0tpqUyipMPGUGxKS9k08XPSfM7J+AUAaOuOG1ltxY0lcpOvSinByJpy+DyKbrIV4kYeT2Ft4gcAsVjyYE+5VNxqKJbPHc4wPFPubWTFLnIDABt2pa582qmnpeLN+/SquxDFDSlNikLcLF68GOPGjUNZWRlmzZqFd955J+W+f/3rX3HMMcdg0KBBGDRoEObMmZN2f0L6OqbITRpDsbiZ+r0evYlfqptXUuQmao5CiOqkIZUBDKkMAADe3bQ341pjmlNxE7++vaE4/lmItXyLm7+fcwSevOhoDE48TzFbqrXbqCQSIkOeYh5VNaOvjxjPoEh9btRkz5DTtFQ6w7W1gkqcM1XkJhJTsbs93sCvtjou4ES/pLBUlk9IKVFwcfPII49g0aJFuOaaa/D+++/jkEMOwbx587Br1y7b/V9++WUsWLAAy5cvx4oVK1BfX4+5c+di+/btvbxyQooDOUWhR27S9bnxelCZEEHt4ahtN1y7tJS1ZwsQv5kff8BwAMCLn9r/zcrYTfEW+KXHdHETsDMUJyq0ImIid37/jdUPrsDUUdX69yKN1JpYY9DnMZWiGyMYbCI3ibXbTQWXzy37m+yQBa0Vq+CZNjq+9lS9bna3haBp8XUPrYyLG5HeBOxHeRBS7BRc3Nx+++04//zzcfbZZ2PKlClYsmQJKioqcM8999ju/3//93/4yU9+gunTp+OAAw7A3XffDVVVsWzZsl5eOSHFgTEh2vDcpO1z41MwMBgfnKlp9kLImpYKR9WUKaU5B8bFzbK1jRnHBogbuscm2uLxKPq5RVSkzObmLQ7VPTd5jtxYsaalrGJC9tWk89xERM8fKWLlPC1lmMitWNNSR4wfDCB15EYuAxevS8Dn0YVWZ4SmYlJ6FFTchMNhrFy5EnPmzNEf83g8mDNnDlasWOHoHJ2dnYhEIhg8eLDt9lAohNbWVtMHIX0J+V280Z8knaE4XgouNIE1SgMYN26/ZHC1TsAWHDN5GAJeD7Y0dWJDhk64YhxBKkEirieuX24TuRECoTf63NhhiJuEALOICb0LcUzT03A+i7hRpciNHLHyOa2WchG5mTluEACgqSNs+1o3ikqpRBm4QPidrD2PCCkFCipu9uzZg1gshtraWtPjtbW1aGhocHSOyy+/HCNHjjQJJJkbb7wR1dXV+kd9fX2P101IMSGnpcr1Pjd2kRvDc6Moiu67sbvhicjNoIq4zyRuKE6ICctcqMqgD0dOHAIAWPZpY9q1WtM0VsSNXhcONh14PZYOxflOS1kRnpvMkRsj2uWxi9zoaSnZc+OsWiqd58bv9Zh+vvWDKvQ17+sIJ+2/01IGLtB9WWzkR0qQgqelesJNN92Ehx9+GP/6179QVlZmu8+VV16JlpYW/WPbtm29vEpC8oeqasYcJp/hpbHzScieGwAYKMRNt424CZvFTSSq6t2F7aIuIjW1fF16303MxrcjI4ZHCsHlLHKT9pI5J2CN3PjsxY3cZdine24MQ7G1QzFgCKdMJl5rZ2QrcjpvyIAgBlXE05DNncnjFKxl4AK9ISQjN6QE8WXeJX8MHToUXq8XjY3md3uNjY2oq6tLe+zvfvc73HTTTXjxxRcxbdq0lPsFg0EEg8GcrJeQYkO+gQb9XlT4E31ubG5I8mwpwHhnbvXXAEZUYlClXz82XXfhg0ZWAQB2SPOM7NDFjc1UcEA248a/t4tMiOsLAWc3XymfiLSUsUazwBDPQRYoVs+Nuaw+eThnuvELmqZJPiv7517m96IjHINHAWrK/RhUEUBjawj7OlNHbkZU20duuui5ISVIQSM3gUAAhx12mMkMLMzBs2fPTnncLbfcgt/+9rdYunQpZs6c2RtLJaQokfvUlEml4F2RmGnQJWDucwMAA8riN6+2NGkpUf4clvrc2EVdqsvjIkguj7bDruJKxm8Jw9hFJoRAENElO99JPvFbhJnV1Gs3tdxO3ERsStnFc0nVORowC5+UkZvEmgZVBODxKHoEzk7cWCeCC+w8N+9s2ounPtiRcm2EFAsFjdwAwKJFi7Bw4ULMnDkTRxxxBO688050dHTg7LPPBgCceeaZGDVqFG688UYAwM0334yrr74aDz30EMaNG6d7cwYMGIABAwYU7HkQUghkU63PaxiKgbjAEe++gXivGsBI/QxIE7lpD1nTUlJZs03UpaosIW66IlBVzbYaCjCmgqf03FjOna5DsYjcpOrImy/8PqsAs4gbJbW48Ul+HOEdkgVdMMO0dsA86iJV5EacZ8iA+OsnInB2aammRI+bYQPMEW47z81P/u997GkPYebYQRhZU55yjYQUmoKLm9NOOw27d+/G1VdfjYaGBkyfPh1Lly7VTcZbt26FRwrb/vnPf0Y4HMZ3vvMd03muueYa/OY3v+nNpRNScEJSGTgQfyevKPES745w1CxuLJ6bAWkNxfGbq13kxq5HTVUicqMmrjswIXaspPPtAMnm4HRTwQsVuQlY12i5vhB/8uBKIXi8UtpNrxzzypGb5KGXVoTw8SjJYtBYU/w84vWrSRO5aemKCx4hgARWz01XOIY9CSG0qy1EcUOKmoKLGwC4+OKLcfHFF9tue/nll03fb968Of8LIqRE0CeCJ0SAoiioSPgtrKZieXAmgLTVUuIxcVMMR1W9dNku6lLm9yLg8yAcVdHSFUkjbjJFbsxCwUnkpvfTUukFmN5kME1aynQ+SSwK4ZRW3Eh+GyVFek+kq4YkmvKlMhRrmqY/VlMeMG3Th7AmUp+720L6NiGICClWSrpaipD+jjVyAwDlAXtTcdhiKE6XlhLCoabcMBRn6lGj+266UhtQM5aCWyIR6QzFHYVKSyWt0fxvVDw18fOWx1VkGjuhp6XSDKzsTjMx3ViTOXKTynPTHorqr0lNhSVyExSRm/jruavNMIu3UtyQIofihpACsHztLlzxzw8dT9NORSiafIPXy8EtVS4Rq6E4RSm4qmr6jXlgwnQclscvpBAmVYl9072rt85asmKNithFZZKrpQrT50aQKnIjPDde02iG5LWa01IOIjdpJoJbzyM8N0ZayvzaiKhNmd+TdD7rKA9GbkgpQXFDSAH4/bLP8fC72/DKZ7t7dJ5uu8hNis6ycodiwKiWarfsJ99YRXop3ucmfUrJScVUVE09fgFIrhyy2y+pWqqXIzdJnpsU1VJCeMo/r0xpKSeem3QTwQXDB8Yrn8YMrgAgp6XMkZtUKSnAmC9lRG4McZOpKo6QQlMUnhtC+huiAZwwaGaLXRv+VJ1lRUdccXM2PDfmG5UcTTIiN1pazw1gmIrTvauPZYjcBEzpNXvRIqqMxERzu/lT+STJc2M1FFtKweXnGvB5MHZIBbY0dQKIp7A8HrvITeqInpPIzf/M2w9HThyMkw4eAcCI3Oy1dChu7gontid7pMRk8M7E74OclmLkhhQ7jNwQUgBEqH9ve3L1ihvsBiimmgxu9dwMDNp7c7qlgZRCYISjsYzCxPDcZBY3Xpv0jPXcdmbi+LHm59PbkRtrKXimPjfW6NOXJg3Vv7amqQzPTc8jN9+cMVqPBKUyFOuRGxtxUyFEciJyI6el0vmqCCkGKG4IKQAi1N9kM+vHDXaRG72E11otJdJSlmopaxM/eW5RQB/kaAyBTO25cSNu7LfLpeCpIhPW6xe6iV+y58ZcCm4VgyZxY20I6KgU3L2oE4bi9lDU1H9HpKns01Lm3yNTWoqRG1LkUNwQ0stomqbfMOz6jrhBv9FJzdz0Et6ktJTZc1NpqYYRyNU4fmlKdbo+N4ARuUmbltLSR24CjsRN+shJvkn23FiqpSyRG6sYO2qiIW6s0TUnhmInkRsrVeV+CF+zSEUB6SM3xoT55MgN01Kk2KG4IaSXCUUNoWD1QLjFriy4PMXAQ6vnZmAwfkOzVksZN0+v7oGJSgM6U3tu4jfDVptBnIJohtSWHMlIVeps7VtX6D431i7B+mypmP3Pq7rCnzLl5qRDsRPPjRWvR9HFp5yaau4S4iZ15KbTLnJjYygORWNYs6MVmqYlbSOkt6G4IaSXkd+tN/XQc2MbuUmkm97fug/H/+5lfRaQdbZU6siNuHl6TCkYIXoyVUuljdzE0ldLycIhlQCwHltocWMVYXq1VJqf16wJg23P7aZaym1/H73XjSSoReTQNnKje27ifqum9vSRm9/+Zw1O+sNrGSfDE9IbUNwQ0svIYqLHaSndUJzsuXnt8z3YtKcDT63eDiB5KrheCh6Omt5ty12P5Ru5EGU98txkaAToN0VuUqSllPSel3yT5LnxpTAU2wzGFHzr0NG253bT58atqBMCRu5106KXgtsYinXPTRRNHSHIc1jtXuNPd7YBAN7auNfVugjJBywFJ6SXMUVuOsLQNC1lG/1M6EJEusFWWEqohWk5bDEUiyZ+mhZfk3inHrIxFANGG/5M1VJOSsGdjF9IXS1V2MiNoigIeD26eMlkKLaLUs2fNgItXRHsN9w87FeP3PSwQ7EdInIj97pJn5YSnpsYdrXGozY+j4KoqqG1O5r0eys8OZ/ubHW1LkLyASM3hPQyHZLRNxxVk6qa3GAfuTG/ZxG+HuG1ECKk3O/VRwXI86WMrsfxJnriZi3SVZn63Dhp4meNvghks7LzaqnejdwA1giTdfxC6j43AkVRcMaRYzFrwhDT4+J1DKdNSyWnIp1gF7lJl5YS4xdiqoYv9nUBAMYOqdAfk39vNU2juCFFBcUNIb1Mp8Xo25NeN04iN3vb49Eh0TBw2ID4MEVFUWyHZ+qRgcQ5halYVF/1JHKjj19IMc3a73NgKPZY01K9/29M7nWTFLnxWvrcuIjKORqcGc1d5KYlXZ8b6XltaeoAANQPrtDXKL/OHeGYHtnb0x42NfwjpBBQ3BCSJ/72+ib8f8s+T6oe6bCUaO/tge/GSeSmLRTF3o6wfsMcmhA3gNzITxY35nSLSBWJm1cqM7CI3HRH1JTVPlE1/Q1fHkWQ0lCsFEPkRhY3VkOxWaCkEnJ2OKmWyjZyM0iP3MR/3zRN09NSg2zSUj6vR0/5bU6Im+EDg0aEThI3cpk4AKxN+G8IKRQUN4TkgfZQFP/vmTW47YXPkqpHrGMR9nZkP4LBfvxC/KanKMaE6s8a2+PbAl7TWINKm+GZ1j4qhrhJnWYB4kJJ6I5UHWwzdTn2OTEUW9NSBYjcyF4kq8gQT8EonXe+PsNzk/vIjXV4Zlsoqr8e1TaGYsCIAm7eEx8XMWxgUC/5b0kjbpiaIoWG4oaQPNDQ0gURsLn1uc/0dAyQ3H+mJ+XgdmXBk4YPQNDnwfH7D8eQRJTms8b4O+mhA4Om443hmcmRG9EcL2hJS6W6WXs8ih4JSuW7cWUoTjFbKjktVWjPjbVaKuGbEX1uXHjFHVVLZR25MaelWtJMBBeIKOAWPXJTZjtmg+KGFBsUN4TkgcZW45/9pztb8Z+PdurfWyM3PSkHF0JJRGAAYER1Od755RzcdcZhGFIZv6GtE+JmgEXc2HluLJEBcSPPVC0FZB6emVncJA+RtJKcliqA5yZNWsqolhI/LzeRG0MYyYJYJmvPTWX8tWmyGMztUlICEQXc0RL30IysKddL/s2Rm/h28fv0KdNSpMBQ3BCSBxoSNwNxH777tY36tqTITQ+6FIuJzRWWd97VFX74vB4MToibzxqEuDHfyAbYem7M0SBxI9/XEb+ZpYqoAJlNxdEcRG6s4qoQ4iaQxlCcafxCOuRREiLyY6UryyZ+o2vilU5f7OuCqmq6qE6VkgLM/i2/V8GREwbbvsaie/FRE+PVXxt2t6f1DRGSbyhuCMkDDa1xcXPQyCoAwPZEKS1gRG6E8OlJtZSY2CzKdq0MzhC5sRueqRuKLdVS4jlZBZJMpkZ+meZT+bLqUFw4Q7HXoyR1LPb1RNxIoilVakqP1gXctSkbWVMGn0dBOKqiobVb6nGTWtxUSr9XR04YgoFlftsxGyItNW10NarL/YiqGjbu7nC1PkJyCcUNIXlgV0IITB4+EADQJt0IRH+QuqoyAD2bLyUaAlorpAQiLSWunyot1ZGizw1g3MjFOodUms8hY+fHkBHiJlWmJuCyQ7FHSe4Y3BsIQ3GZTdTIa2ni50bc+DyKbgJPFfnoyCBoU57b60H94Hj0ZnNTB1oSkZt0aalyv/F7dcKUWgD2r/HuRJuB4QPLMCzh6+pp921CegLFDSFZ0hGK4tz73sXD72xN2iaiHJMSHWjDMVVP94hoy+hB5QB6VgouokCVKVI4QyxiJslQbFMtFbKUglunYA9JF7nJMDwzY+TGQSm4vJygz5t1d+eeIPrx2AkwryVyk86jZEVRlIwVU0LQDgi6bzAvmvBtaerUq6acRm7mHBgXN3bRORG5GTYwiKqESb0tzQBVQvINxQ0hWfLmhiYsW7sLi19en7StIWEonjisUk8/iX/2InJTPyh+o+lJ5EacqyLFjU6kpQTDLMLEaOJnRAms7f0DluiENfojk8lzI2ZLpfTcSNdKVeItp6UKUQYOGNEsO3FjTUul6guUCqPXjb24EebvVNG6dIwbUgkgHrkxOlanFqviGgeNrMLImvLE/naGYkPcDHQwY4yQfENxQ0iWNLTEfTTb93UlpRBEWmpEdTkGBMQ72fg/exFtGZ1IEWTruYnGVP0GajUUC4ZYxE1SWkovBTduRPLgTCA57ZNO3OiVNJ0pDMWxDOJGejxV5EaOhFiHVvYWQtzYiSuPJS3lJnIDyOXg9mkp8fvTk8jN5j0dUnl36tdzv9p45PGbM0bpj1VbxmzEVE03xQ+Tm/wxckMKCAdnEpIlOxMVUaoGbNvbiUkJf01M1fTqkdqqMgwo86EtFDUiN4koiUhLtYWiCEVjro2xndJwxUyGYoFVmAxO+C12SX1KjA649pEb6zllRIqjuctesKkZpoL7TCXWmTsUFypyY3hu0kRuYllGbkRayiZyE4rGEEkIRLeeG8CI3Gzc3aGnTmeOG5Ry/x8eORazxg/BgSMG6o9Zy/33dYYRUzUoSvx3w0hLMXJDCgcjN4Rkibg5AMCmRAdXAGhqDyGmavAo8cqigZZGeeKd98jqcj2C0Zwi0pEOMaPK51GSfDECqz/G6rkZP9S42YkxEXqHYkspOABUlfmSxI6MtQuulcyl4JkjN/KxhSgDB4x12vWasT43t5Eb8fO189zIc8ncVksBwLjE6/35rna0dUcxMOjDlBFVKff3ez2YMrLK5GsyDMXx32ORkhpcEYDf65HSUozckMJBcUNIloheNgCwaU+78XhC9AwbGIRP+mcv3smKyM2AMp8+7yebLsVCJJUHUptqB0uVTWV+T5LxWIiblq6I7v1JGpwpiZt0KSnAfjijTK47FBeiOzFgCBBbQ3EPmwymS0sJgVzm97iqwhKMqik3HXfE+MGmaJkThLjZ1xk2TQMXVVLCVM7IDSkkFDeEZIlZ3BiRG9GduDZR6i0iN626odiocBIpHjtTsapqiKRo5AYYVTPp3sHXlPv10uKhA4JJIqg84MXI6rLEc4h7MIzBmclpqUzipkYfztjz8Qup/DSmtFTBIjdpxI3Fo1RXXe7q3OlGMDh5zdMR8HkwqsZYz5EThrg+x/CqoL6+vR3hJHGjR24obkgBobghJAs0TbOkpZIjN4a4EZGbRFoqZFQ4iUjH3sS74BUbmnQvw81L1+LAXy/F2gb7OT16v5M0HYM9HkUXUKmEyYRhcdOoaLoWSjIUG/8m0pWBA4YfpznxfKy4GpwZSFUuLoubwhqK7dJS1uc2qqbM1bnTeW5E5KYyCzOxQJiKgezETdDn1U3IO5q7db+WHrlhKTgpAihuCMmC1u6o/i4aMKYmA0BjIqJTZ4nctHVHEI6qutG0MuDVxcLe9hDe3NCEBX99Cxf93/voCEXx9xVbEFU1vLG+yXYN+uiFDMbSzOImnprakBBo1j43bsSNEGuRmKaXqctE1YTJNkUaTaTAFCW5v47AY0pLFchQLNJSNuLK+txG1riM3IhS8Ejyz0+kItMJ2kwIU/HAoA9TRqb226RDPKftzZ34Yl/8d39EIgJYxcgNKQIobgjJApGSkkcTiBtPYyJyU1dtFTdRdEk3/IqAz5SWEpOUX1+/B3946XN9htC2vYZwktEjQBlSFOIawwbaC5MJkqkYSB6cKael0nUnBuJpLpFW2Webaot/9qXoKiweL/en9hF5lWKI3MTXEEzT50bgWtykSUuJaF02ZeACUd595MQhWfl2AGBUotLvi31d2NIU//0Uosnw3DByQwoHS8EJyQKRepowtDI+p6czgs17OjFlZJW+bbgepjcMxcJvE/B6EPB5dMNvU0fYdDP7yyvGoE3xztiK03fxQpBkSktt2tOBmKrppcaGodi4AVqrrewYVBFAQ2s39nWG9Xb/AhG5sZpuBbUDyxDwejA2caO0oxiqpQ4fNxgVgc2YPTE5reO1RJxGuvbcxH/uYVtxk75poxO+d3g9uiMqTpo2IutzjNIjN126V0tUYrGJHykGKG4IyQLRwK+uugzlAS9WbW3G5qaOuLhpSR250QVJIpU0uMKoPOm0SeMAwLa9XbaPOzWXfvXA4Xhzwx4cM3mY7XaRltrS1KGLL8A+LTU0TY8bQU2FPyFuzDc3TdOQsNykjBgMqgzgxUXH6u/+7TB3KC5M5Oa4/Yfjo9/Ms30esnAbXBlIO0XdjrSRG72BX/bPO+jz4vwvT8j6eMAQN5v3dGBH4m9Bj9yUGU38NE0zReDCURWbmzqwX+1AEJJPmJYiJAsaWuImyhHVZRif+Ke+aU8HojE1KUwvixvrROfBiWhKU3sYu9oMgzJgRH6+2Ndpa87tkErB0/GtQ0fj/V+fgCPGD7bdPrK6HEGfB5GYhvW7DGO0XRM/66wqO1KVgwszMZB6thQAjBlSoffLscNXBJEbILVAk9c30qWZGJDHLySL3Q6Hqch8I8TNu5v3QdPiaTIxLV78vsdUTU+tCu588TPMveNVLHllQ+8umPQ7KG5Iv2NLUwfm3fEqHno7eeClUxpa4+9Wa6vK9F4xG3a344t9XQjHVASlktuBweS0lEglDZE8N7sSJeRnHTUONRV+XD1/CoD4/Ci70uouPXKT+V18uuGSHo+iP4c1O+K+n4DXo0dI3BiKAWBQZSIaZfHcRCVxk0bbZKQYOhSnQ44suU1JAemrpXLhuckFwkckqrfGDqnQf8cqAl5d+Fkb+b2/dR8A4PbnP8NnjW29tVzSDym+/wyE5Jn/e3sr1jW24eala00GXzeI0Qsjqsuwf108xP7pzjZs2B2PfEwYNkC/yZnSUhbPhDD77usM60bk02eNweqr5+Lr00bq0Rs7U3Eu/BcCkZpakzA1y6LBTZ8bIHWXYqeRm0yYmvgVyFCcDnPkJhtxk7pDsVUcFwphKBYIvw0QF9IDU4xgECnWcEzFLx7/0PQ7QUguobgh/Y4XP20EEO/K+9QH27M6h+GrKdfLaT9vbNMrniYOM/7ZD5Ca+MkN/ABD3DR1hPUmf8MHGqkMYcjdZmMq7ookbnQ58J0IA+/6xrg4k5vTiZJsv1fRe5ikQ3RdTkpLSam1bKt04scaXxdj5Mbb07SUHrmxS0v1vM9NLqgq85miR+MtBnC7cvBITMXOhD8n6PPgg23NeHfz3l5YLemPFN9/BkJyiKpq+jgBIO6LESXPAHD/m1ts/SyZEBVRdVVlGFVTjupyP6Kqhv9+3AAAmJioQALM1VLCBCw8E8KfIpYQ8HlMZtr6xDtkO1NxLiM3okfJ5sSkaLl/jN8Xv1kPqUzucGzHIEvkJqZqWHDXW/jun1fo+/RE3HiKoBQ8Hd6eRm786QzFzlOR+URRFFOnY7kxIGCUg8uTwXc2d0PV4r/jwv+1tcm+EpCQnkJxQ/o0//vYB5j5/17Uu/wuS0Rtpo2uRpnfgzU7W/Heln2uztkdiemDLuuqy6Aoij588JOEZ2XScEPciBB9KKrqx1Umql0CPo++HYibiGUBISI3duXguWjoJhDdlEW3WTndE/Am/EEO/DaALG7ikZtXPtuFFRubsE7yWPRA21hmSxXfv7CeihsRKUvnuSl05AYwp6bGDzVHboTPTC4HF9HH0YPKMXpQ4ve62b4SkJCeUnz/GQjJEc2dYfz7gx1oD0Xx+xc/BwAs+3QXAODU6aNw8sEjARhpKqcIoTEg6NPTNAdZOr3KkRs5fC98NXK1y2CpvFqIDMFoEbnZl3wTMKJAPRc3InIjkNNSh48bhAPqBuK7h412dC5hKBZC7uF3tpm2+zyKowhQKrxFMH7BKaN6Ermx61Asqu2KQNzIKTdrXyK7yI3wjdUPqtB/r3dQ3JA8Ufi/EELyxIuf7tINi0s/acBznzToOf6vHjgcaiIXtN1GOKRDlHqPGWxUiMht7BXFMOgCgM/rQUXAi85wTBc3clphcGVAP+dwS5O8evEO185Q3MMhijJ1SeLGeN8zvKoMSy/5suNzCUNxvAKsG8vW7jJt9/QkbANrh+Lie38mT3gf5sCAbSXfs6Vyxaia+O+mXAYuGCilYgUiclM/uFwXRm7/9ghxSvH9ZyAkRyxN+F8CPg80DbjwgZWIqhqm19dg7JBKU5dVNwghIvsMDhpZrX89elB50rRokXoS71Rln8wQKXKTJG70tFQXVEtlSVcO01JDK4OmKh+7addOkfvcPLbyC8RUzXTuVEMzneLN0TrzRVgyAmcj5IwmfqlnSxXacwPERQoQT0lZI3G6obhLjtzEf/frB1Xowsjt3x4hTqG4IX2SjlAUr36+GwBw/alT9cdnjKnBX8+cCcDwDLgNjW9NRFHGSOJmwrBKvWRaTkkJxDvZj3ckV1PJaanhlrTUiOoyeD0KwjFV98MIcmko9ngUU0qsJ+keUS3VEY7hkXfjKamzjhpnXKsHKSnA0qG4CCM33z5sNA4fNwhXnXxgVsen61DcXiRN/ADgqwfU4gezxuAX8/ZP2mZXCm5Ebir0yM3OlmTRTkguKPxfCCF54OV1uxGOqhg3pALfOWw0uiIx7GkP4yfHTdTf7Quz5662EMJR1dTPJR2iomic5DPwez04oG4gPvyiJYW4Mbq2AsBhYwfp2walidz4vB6MqC7DF/u6sHVvpyl9lOt38bVVQf2ddE9KrKvK/PAogKrFhWDA58FPvzIZd7++CYCRWsmWYk9LDSzz47EfHZX18WKkhF2fm85wcTTxA+KdsW/45sG226rKjREMAjlyU1cVF+2RmIbd7aEkrxkhPaX4/jMQkuDtjU04/+/v4RePfZA0miATL6+L+zzmHVQHRVFw5uxxWHTCfqY0xpDKAIKJlJXoW+MEUb461jIU8oQDawEAx+6XPMNJRG6AuICRjaZD0kRuACP9taWpw/S4MBS7nV2UihFSN92eNMfzeBRUlxvP99j9hqG6wo/yHKWQTIbiIkxL9RQh2MIxs7hRVc0wkfdgtlRvYI3cdIVj2NMejzzWDy6Hz+tBXeJ3nakpkg8obkjR0R6K4rz738Npd72FF9Y04rGVX+Crt72C5z9pcHyOtQ3xsuMZY2pS7iP36nD6Dzamanp4fYylt8dPjp+E966agy/bihvjnfahYwaZPApiMjiQHLkBjEqUrZKpOKZqetoiF4ZiwGwq7mmJ9SBpNtRJB9cBAA4YkZthicVeCt5TUnluOqXqqWKI3KSjyjIZXFQYDgz6dOFLUzHJJ33vPwMpaXa1deP7d63Ai582IuD1YMER9Th4VDXauqO47J8fOkppqKqGz3fFxU2m6cMiNeXUd7OzpQuRmAa/VzFFOoD4TTfVeAK5s6+ckgLMkRu78LyIEG2RGp51StO7cxW5qauSxU3PzlmT8N0EvB58NRHROqCuKt0hjimlUvBs0KulLGkp0ePGoxRnOk6mSho5Akg9bqQKw2wN/XY8+t42/OapT0x/F6R/U9x/IaRfEY2pOP/+9/Dx9lYMqQzgsR/Nxo3fmoZ//eQoTBhaiebOCB5YsSXjebbt60R3JO6hsfbfsCLePToVNyIlVT+owlWXXTktdejYGtM24bnxexXdjCujp6X2yuIm/i7e61FydqPLR+TmmMlD9Xfx3z50FICeRx08Re656SmpOhTLDfx60ieoNzA8N/HIjeG3Md4QOH1j0RmOmsY4WGkPRXHVkx/jvjc345z73qXAIQAobkgR8aeXN+CDL1pQVebD4z8+CofU1wCIm2ovOn4SAODu1zZm/Oe1LpGSmjRsQEYB4rYkdYtNpZQTxA3d71VMZeMAMHn4AIysLsOx+w23vWmNGRwXaLLnRtzoKgLenN3o5EZ+PR1IOXPcYCgK8MMjx5oee/TC2XjuUuc9c+wo9lLwnpIqLSWq43KVhswnIg0rSsHt2ieIasV0aamOUBRf/8PrOO7Wl3XPjpXla3chnBCCb23ci4sfWtXzJ0BKHoobUhR8vL0Ff1gW7yJ83SlTk9q5nzJ9JMYMrkBTRxj/sHS8tfL5rvjwRzGtOx163t+puElhJs6ECNNPHVWddEOuDPrw6mXH469nHmZ7rLghNHdG0JLwMOSyO7GgNodpqR8dOwGrr56L4w8Ybnr8iPGDs+raK+Mr8lLwniI38ZPnnulDV4vcTAwANeXxyF1XJIbOcFSvMJQjqU7SUne++Bk27unA3o4wHnzLPmr7XMKLd/z+w6AowEtrd2F3m70QIv2HvvefgZQc0ZiKK574EFFVw4lT63DK9JFJ+/i8Hpx/zHgAwL9Xp5/kLSI3k2uTS7KtjHIYGn9z/R787B+r8Fqid86YDOkuK189sBZHjBuMC788wXa7z+tJGYGpDPp0L49Ii3XmsDuxwCxuevavQVHMFVO5pNj73PQU0ZJA04BITBI3RdSdOBPVFX69f9PG3R26uJHftGQSNx9vb8E9b2zWv3/wrS1J0azuSAzLEx2wfz5nP0xOzHRbva05J8+DlC597z8DKTn+vmILPt7eiqoyH647ZWrKm/zXpo6AogAfftGSVox8lhjQuH8GMzEghcabu/R3yS1dEdMk8a5wDD9/ZDWe+mCHPhjTbeSmfnAFHv3RbHxt6ghXxwkM3038JiFSc7kyEwPxm6oQUcVcYl3h90JR4qk+n7fv/QuTBZt8M8/luI3eYFKi39O6hjZ9rpRdWqqtO4p9HeGk429euhYxVcO8g2oxsroMe9rDeGr1DtM+b6zfg45wDCOqyzBtVDVm1MfN+qu2uhuGS/oefe8/AykpdrZ04bbn1wEALj/xAAyzKYUWDBsYxGFj4v+8XlhjP+wyElOxcXdcAGSqlAIME213RMW+zgiaO8M45uaX8JXfvYxNe+LnuX/FZuxuC2FQhR9lfg98HgUHjcpN5Y9TrBVT+YjcAEBddfznX8xelkGVAdz6nUNw52nTC72UvGAWN4ap2IjcFO9rIzMpETl99fPdiMQ0BLweU4VhRcCHcQmxs/qLZtOxezvCeHNDEwDgyhMPxJmJDtd/e32TqaPxfxMjVuZOqYXHo2B6ovUDIzeE4oYUDFXV8L+PfYCOcAyHjqnBgsPHZDxm7kHxsuLn19j3vNnS1IFwTEVFwOvI2xH0eXVBtaO5C2+sb0JrdxQ7WuIl6U+u2o4lr2wAAPzq5ClYccVX8fIvjksqA883YyyN/HRDcY5vdKfPGotD6mtw5ITBOT1vrvnOYaMxZ0ptoZeRFxRF0b1UckSjlNJSgBG5eXldPJVbP7g8yeB/aOLNyqot5kjLi582IqZqmDKiCuOGVmLB4WMwMOjD2oY2PPvxTgBAS2cEz3wY//rrh8RT2aKv1QfbmvVu4KR/QnFDCsY9b2zCG+ubUO734nffPcTRkMF5B8Ubwr21cS+aO5ND2Z81xs3Ek2sHOh5aKETQht3tWLFxD4D4ZO/G1hAueWQ1mjsjmDR8AL45YxQGVQYwepC7lFQuEKMe3t28Dxf8/T3cv2IzgNwaigFgwRFj8O+LjsbwgWyHX0gOHhWvqHtPuukbJvLSEDfC8yZM8ONsfGqHJno+vb+12fT4c4mIzNemxv/eqyv8OO+YuF/t9uc/QzSm4rGV29AVieGAuoGYmTjP5OEDURnwoiMcw/pEYQHpn1DckIKwaus+3LI0no666usHYoLNPCY7xg6pxAF1AxFTNdvUlHgnN21UddK2VMxKRCme+6QBKxKh8OtPPRjfP7weR4wfjOn1Nbj+1Kmu+trkGhG52bSnA8+vacTH2+Pen1K50RF3zBof/518e2OT/tgHiVTL4Mr8GLVzzaTh5r/pcUNtxM0YwyMjIi3toShe+zz+JkOIGwA495jxGFwZwMY9HfjLqxv16qkzZo/VfXpej4Jpo2v0czolFI3h+U8acOeLn+H3L35uO5GdlBYUN6TX2dnShQseWIlwTMXcKbX4wRGZ01EyX58WN+X+8/0vTI9vberEfxMha7m/SibmT4uHtF9cswsbdndAUYCTDx6Bm749DY9eOBtPXnQ0Zk0Y4mqNuWbS8AEo93vh8yg4c/ZYzJ1Si4DXg9kFXhfJD+L37e1Ne6FpGj7e3oJla3fBowDfOnR0gVfnjLqqMlPDxnE2vaH2rzMiLaIQ4KW1uxCOqZgwrFKvfgLiBvKfHDcRAHDrc+uwuakTA4M+nDp9lOmcVt/Nk6u245Q/vo5Pd7barnN3Wwjf+8tbuOCBlbjzxc9xx4uf4e7XNmX9vElxUBTiZvHixRg3bhzKysowa9YsvPPOO2n3f+yxx3DAAQegrKwMBx98MJ599tleWmn/ZF9HGE9/sAM3/Xct/vTyer3yIRv2doRxzn3vYXdbCPvXDsTtp0133YTum4eOhqLEU1PyWu55YxNUDfjyfsMc9bgRHDSyCuOGVOiDCqeMqEK1TafgQlJV5sdzl3wZr152PK47ZSruOnMm1v2/r+Hbh5XGjY6449Axg+DzKNjZ0o0v9nXhzhfjPaBOmT7Kdup8MaIoCiZK4sSuW7hXMgG/v3Uftjd34fpn1gAATpxal/S/4eyjx2PRCfvpfaMWzBqT5EESRQf/+XAnln68E5f/80N88EULLn1ktd7sL6ZqeOTdrfjlvz7CKX98HR9sa0Z1uR/H7x+fC7d4+XrsarUfptsdiWHjbqa8ip2Cx7QfeeQRLFq0CEuWLMGsWbNw5513Yt68eVi3bh2GDx+etP+bb76JBQsW4MYbb8TXv/51PPTQQzj11FPx/vvvY+rUqQV4Bn2bZz7ciV89+RGaO43257csXYfDxw3CNw4ZiaEDgtjXGcGqrfvQGY7hpINHYM6U4bYzf3a2dOHMv72Dz3e1Y0hlAHcvnJlVK/5RNeU4euJQvL5+Dx5f+QXOPWY8nv5gBx55N97c74Jj7HvJpEJRFHx92kj8cfl6AMCRRRoNsXZFLvYW/CR7ygNeTBtdjfe3NuOOFz/Di582wqMAF39lUqGX5orJwwfo6TQ7zw0QF3JvrG/C4yu/wD2vb0Jjawj71Q7ABV+emLSv16PgZ1+djLOPHocPtrXYGt+P238Yjhg/GO9s2osfPfi+/vjahjb88aXP8bOvTsYvHv8Q/1pl9MsaO6QC9519BMYNqcC3/vwmVm1txs1L1+F3351m+jt7b/NeLHr0A2zd24nj9h+Gb84YhfW72vHx9hZ81tiOoQMCOKS+BtNG12DGmBpMGFrJv9MCoWhyC8wCMGvWLBx++OH44x//CABQVRX19fX46U9/iiuuuCJp/9NOOw0dHR34z3/+oz925JFHYvr06ViyZEnG67W2tqK6uhotLS2oqspdOW93JGbqiil+qho06WuxTdO/N376muU48znkbTFVQ2c4hq5IDF3hKDrDMYSjKoYOCKKuugyRmIqGlm6s3taMhtZuhCIqQtEYQlE18RFDKKLC61EwcdgA7F83EPvVDsSwgcYAx8bWEB56eyteXx/PfU8YWokjJw7BlqYOvLmhCel+a4ZUBnDOl8Zj5thB6AhHsW1vF97bsg/PfdyAcExFbVUQ/3feLEwanv2U6H+v3o6fP7waA4M+RFQV3YkhgzPG1OCJHx/l+h/K2oZWfO3O1wAAd585s89W4pDS4eala/Hnlzfo3582sx43f2daAVfkniWvbMBN/10Ln0fB2t9+zbYv0fK1u3D2fe/q39dVleGJnxylz57Kht1tIXzjj69jZ0s3BlX48bOvTsa1T8cjQlVlPrR2R+HzKDjnS+MxefgAzD2oTm86+f7WffjWn94EAJw8bQQuOm4SuqMxPLhiC55cvR1uirAmDKvECQfW4qBR1RhVU46gz6OX+q9rbMPKLfvw+ud70NQRxsyxgzBrwhAcNLIK1eV+RGNa4n9b/N7S1B5GVySmD+gdNjCIoQMC8Hk8iKkaVC3+0VM8igJFiQtJr6LA41GgIPWbqVT/aYN+T86rSt3cvwsauQmHw1i5ciWuvPJK/TGPx4M5c+ZgxYoVtsesWLECixYtMj02b948PPnkk7b7h0IhhEKG6Ghttc+79pRPdrTi239+My/nzidrG9rwzEc7U273ehT85LiJ+OlXJuudUxtauvHUB9uxfO1uRFUV5QEfDh5VBU0Dnnh/Oxpau3Hrc+tszze9vgZ/+P4M17OZrMydUoeBQR/aEuWxk4cPwGmH1+N7h9dn9U5p/9qBmHNgLbbt7cTsicUZuSH9i1njB+vi5vj9h+G6Uw8q8IrcIxppjh1SkbLh4pcmD8V3DhuN5s4I9qsdgAVHjOmRsAHiPbH+tvBw/GHZ5zj76HE4YvxgbGnqxP0rNqO1O4qA14M//mAG5h5Ul3TsoWMG4VcnHYiblq7FMx/u1IsUBN8+dDTOPnoc/vraRmze04ED6qowdVQV9q+rQmNrNz7Y1owPvmjGB9tasHF3B/6ye6OjNT+/phHPp+jfVYocOqYGT/zk6IJdv6DiZs+ePYjFYqitNb9Lrq2txdq1a22PaWhosN2/ocG+78mNN96Ia6+9NjcLToPXo6A80fhM3FvFLVZRFEPd2m2zeUx8r1gOVBTAm+iDUR7wotwf/xzwetDY1o3G1hCCPg8GVQRw8OhqTBhaiaDfq79jKNO/9qI7GsPnjW1Y19COzxrb0CZN3vV5PfjqgcPxw1ljUW/pxltXXYYLvjzRNmx86Qn74ekPduDvK7agtSuCiqAXI6rLMXHYAJx88AgcPNp5FVM6ygNe/OEHM/DWxiacOHUEDhld3aPwr6IouHvhzJysjZBccOSEIThywmDUVZXhpm9Ps031FjvHTB6KC748AUelecPg93rwu+8ekvNrTxlZhSVnGPPafvONg/C/8/bHJ9tbMGRAMKmaS+b8L0/ArAmD8ZunPsGWpk5EVQ1HTxqCHx07Ua/G+v33Z9geOz/Rc6etO4KX1u7C25v24tOdrdjTHkI4qiIcVRFVNYwfWomDR1Xj6ElDUVtVhrc2NmH1tmZ8urMV3REVfq8Cv9eDgM+DYQOCGDIggMqAD5GYit3tIexpD2NPewiqqsHjUeBRElGXHv7cVA2IaRo0TUtEhJAcEbIJEFkfymX39GwoaFpqx44dGDVqFN58803Mnj1bf/yyyy7DK6+8grfffjvpmEAggPvvvx8LFizQH/vTn/6Ea6+9Fo2NyarXLnJTX1+f87QUIYQQQvJHyaSlhg4dCq/XmyRKGhsbUVeXHC4EgLq6Olf7B4NBBIOpW/oTQgghpG9R0FLwQCCAww47DMuWLdMfU1UVy5YtM0VyZGbPnm3aHwBeeOGFlPsTQgghpH9R8FLwRYsWYeHChZg5cyaOOOII3Hnnnejo6MDZZ58NADjzzDMxatQo3HjjjQCAn//85zj22GNx22234eSTT8bDDz+M9957D3fddVchnwYhhBBCioSCi5vTTjsNu3fvxtVXX42GhgZMnz4dS5cu1U3DW7duhcdjBJiOOuooPPTQQ7jqqqvwy1/+EpMnT8aTTz7JHjeEEEIIAVAEfW56m3z1uSGEEEJI/nBz/y6K8QuEEEIIIbmC4oYQQgghfQqKG0IIIYT0KShuCCGEENKnoLghhBBCSJ+C4oYQQgghfQqKG0IIIYT0KShuCCGEENKnoLghhBBCSJ+i4OMXehvRkLm1tbXAKyGEEEKIU8R928lghX4nbtra2gAA9fX1BV4JIYQQQtzS1taG6urqtPv0u9lSqqpix44dGDhwIBRF6ZVrtra2or6+Htu2beM8qxKCr1tpwtetNOHrVpr05uumaRra2towcuRI00BtO/pd5Mbj8WD06NEFuXZVVRX/aEsQvm6lCV+30oSvW2nSW69bpoiNgIZiQgghhPQpKG4IIYQQ0qeguOkFgsEgrrnmGgSDwUIvhbiAr1tpwtetNOHrVpoU6+vW7wzFhBBCCOnbMHJDCCGEkD4FxQ0hhBBC+hQUN4QQQgjpU1DcEEIIIaRPQXHjgra2NlxyySUYO3YsysvLcdRRR+Hdd9/Vtzc2NuKss87CyJEjUVFRga997Wv4/PPPTefYsGEDvvnNb2LYsGGoqqrC9773PTQ2Npr22bt3L04//XRUVVWhpqYG5557Ltrb23vlOfY1eus1GzduHBRFMX3cdNNNvfIc+wKvvvoq5s+fj5EjR0JRFDz55JOm7Zqm4eqrr8aIESNQXl6OOXPmJL1OTv5uPvzwQxxzzDEoKytDfX09brnllqS1PPbYYzjggANQVlaGgw8+GM8++2zOn29foVhet/vuuy/p76+srCwvz7kv0BuvW3d3N8466ywcfPDB8Pl8OPXUU23X8vLLL+PQQw9FMBjEpEmTcN999+XkOVLcuOC8887DCy+8gAceeAAfffQR5s6dizlz5mD79u3QNA2nnnoqNm7ciH//+99YtWoVxo4dizlz5qCjowMA0NHRgblz50JRFLz00kt44403EA6HMX/+fKiqql/n9NNPxyeffIIXXngB//nPf/Dqq6/iggsuKNTTLml66zUDgOuuuw47d+7UP376058W4imXJB0dHTjkkEOwePFi2+233HIL/vCHP2DJkiV4++23UVlZiXnz5qG7u1vfJ9PfTWtrK+bOnYuxY8di5cqVuPXWW/Gb3/wGd911l77Pm2++iQULFuDcc8/FqlWrcOqpp+LUU0/Fxx9/nL8nX8IUy+sGxDvkyn9/W7Zsyc+T7gP0xusWi8VQXl6On/3sZ5gzZ47tdTZt2oSTTz4Zxx9/PFavXo1LLrkE5513Hp577rmeP0mNOKKzs1Pzer3af/7zH9Pjhx56qParX/1KW7dunQZA+/jjj/VtsVhMGzZsmPbXv/5V0zRNe+655zSPx6O1tLTo+zQ3N2uKomgvvPCCpmmatmbNGg2A9u677+r7/Pe//9UURdG2b9+ez6fY5+it10zTNG3s2LHaHXfckd8n1E8AoP3rX//Sv1dVVaurq9NuvfVW/bHm5mYtGAxq//jHPzRNc/Z386c//UkbNGiQFgqF9H0uv/xybf/999e//973vqedfPLJpvXMmjVLu/DCC3P6HPsihXzd7r33Xq26ujpPz6xvk6/XTWbhwoXaKaeckvT4ZZddph100EGmx0477TRt3rx5PXxWmsbIjUOi0ShisVhSqLO8vByvv/46QqEQAJi2ezweBINBvP766wCAUCgERVFMzY7Kysrg8Xj0fVasWIGamhrMnDlT32fOnDnweDx4++238/b8+iK99ZoJbrrpJgwZMgQzZszArbfeimg0mq+n1q/YtGkTGhoaTO/+qqurMWvWLKxYsQKAs7+bFStW4Mtf/jICgYC+z7x587Bu3Trs27dP38f6LnPevHn6dYhzevN1A4D29naMHTsW9fX1OOWUU/DJJ5/k+yn2SXL1ujkhn39vFDcOGThwIGbPno3f/va32LFjB2KxGB588EGsWLECO3fuxAEHHIAxY8bgyiuvxL59+xAOh3HzzTfjiy++wM6dOwEARx55JCorK3H55Zejs7MTHR0d+N///V/EYjF9n4aGBgwfPtx0bZ/Ph8GDB6OhoaHXn3cp01uvGQD87Gc/w8MPP4zly5fjwgsvxA033IDLLrusUE+9TyF+72tra02P19bW6tuc/N00NDTYnkO+Rqp9+Lfnnt583fbff3/cc889+Pe//40HH3wQqqriqKOOwhdffJH7J9bHydXr5vRadtdpbW1FV1dXNsvXobhxwQMPPABN0zBq1CgEg0H84Q9/wIIFC+DxeOD3+/HEE0/gs88+w+DBg1FRUYHly5fjxBNP1EezDxs2DI899hiefvppDBgwANXV1Whubsahhx6acXw7yY7ees0WLVqE4447DtOmTcOPfvQj3Hbbbfj//r//T48OEULyx+zZs3HmmWdi+vTpOPbYY/HEE09g2LBh+Mtf/lLopZECwTuqCyZOnIhXXnkF7e3t2LZtG9555x1EIhFMmDABAHDYYYdh9erVaG5uxs6dO7F06VI0NTXp2wFg7ty52LBhA3bt2oU9e/bggQcewPbt2/V96urqsGvXLtN1o9Eo9u7di7q6ut57sn2E3njN7Jg1axai0Sg2b96c76fY5xG/99YKtcbGRn2bk7+buro623PI10i1D//23NObr5sVv9+PGTNmYP369T1/Iv2MXL1uTq9ld52qqiqUl5dns3wdipssqKysxIgRI7Bv3z4899xzOOWUU0zbq6urMWzYMHz++ed47733krYDwNChQ1FTU4OXXnoJu3btwje+8Q0A8Xcgzc3NWLlypb7vSy+9BFVVMWvWrPw+sT5MPl8zO1avXg2Px5MUuiXuGT9+POrq6rBs2TL9sdbWVrz99tuYPXs2AGd/N7Nnz8arr76KSCSi7/PCCy9g//33x6BBg/R95OuIfcR1iHN683WzEovF8NFHH2HEiBH5eGp9mly9bk7I699bjy3J/YilS5dq//3vf7WNGzdqzz//vHbIIYdos2bN0sLhsKZpmvboo49qy5cv1zZs2KA9+eST2tixY7VvfetbpnPcc8892ooVK7T169drDzzwgDZ48GBt0aJFpn2+9rWvaTNmzNDefvtt7fXXX9cmT56sLViwoNeeZ1+iN16zN998U7vjjju01atXaxs2bNAefPBBbdiwYdqZZ57Zq8+1lGlra9NWrVqlrVq1SgOg3X777dqqVau0LVu2aJqmaTfddJNWU1Oj/fvf/9Y+/PBD7ZRTTtHGjx+vdXV16efI9HfT3Nys1dbWameccYb28ccfaw8//LBWUVGh/eUvf9H3eeONNzSfz6f97ne/0z799FPtmmuu0fx+v/bRRx/13g+jhCiW1+3aa6/VnnvuOW3Dhg3aypUrte9///taWVmZ9sknn/TeD6OE6I3XTdM07ZNPPtFWrVqlzZ8/XzvuuOP0awo2btyoVVRUaL/4xS+0Tz/9VFu8eLHm9Xq1pUuX9vg5Uty44JFHHtEmTJigBQIBra6uTrvooou05uZmffvvf/97bfTo0Zrf79fGjBmjXXXVVabyRU2LlzDW1tZqfr9fmzx5snbbbbdpqqqa9mlqatIWLFigDRgwQKuqqtLOPvtsra2trVeeY1+jN16zlStXarNmzdKqq6u1srIy7cADD9RuuOEGrbu7u9eeZ6mzfPlyDUDSx8KFCzVNi5en/vrXv9Zqa2u1YDCoffWrX9XWrVtnOoeTv5sPPvhA+9KXvqQFg0Ft1KhR2k033ZS0lkcffVTbb7/9tEAgoB100EHaM888k7fnXeoUy+t2ySWXaGPGjNECgYBWW1urnXTSSdr777+f1+deyvTW6zZ27Fjb61jXMn36dC0QCGgTJkzQ7r333pw8R0XTNK3n8R9CCCGEkOKAnhtCCCGE9CkobgghhBDSp6C4IYQQQkifguKGEEIIIX0KihtCCCGE9CkobgghhBDSp6C4IYQQQkifguKGEEIIIX0KihtCCAHQ1NSE4cOH53zY6Zo1azB69Gh0dHTk9LyEkNRQ3BBCXHHWWWdBUZSkj6997WuFXlqPuP7663HKKadg3LhxjvafP39+yuf82muvQVEUfPjhh5gyZQqOPPJI3H777TlcLSEkHRy/QAhxxVlnnYXGxkbce++9pseDwWDKKc25IBwOIxAI5OXcnZ2dGDFiBJ577jkceeSRjo558skn8e1vfxtbtmzB6NGjTdvOOeccfPTRR3j33XcBAM888wzOP/98bN26FT6fL+frJ4SYYeSGEOKaYDCIuro604csbBRFwd13341vfvObqKiowOTJk/HUU0+ZzvHxxx/jxBNPxIABA1BbW4szzjgDe/bs0bcfd9xxuPjii3HJJZdg6NChmDdvHgDgqaeewuTJk1FWVobjjz8e999/PxRFQXNzMzo6OlBVVYXHH3/cdK0nn3wSlZWVaGtrs30+zz77LILBYJKwSbfGr3/96xg2bBjuu+8+0zHt7e147LHHcO655+qPnXDCCdi7dy9eeeUVhz9hQkhPoLghhOSFa6+9Ft/73vfw4Ycf4qSTTsLpp5+OvXv3AgCam5vxla98BTNmzMB7772HpUuXorGxEd/73vdM57j//vsRCATwxhtvYMmSJdi0aRO+853v4NRTT8UHH3yACy+8EL/61a/0/SsrK/H9738/Kap077334jvf+Q4GDhxou9bXXnsNhx12mOmxTGv0+Xw488wzcd9990EOgD/22GOIxWJYsGCB/lggEMD06dPx2muvZfGTJIS4JiezxQkh/YaFCxdqXq9Xq6ysNH1cf/31+j4AtKuuukr/vr29XQOg/fe//9U0TdN++9vfanPnzjWdd9u2bRoAbd26dZqmadqxxx6rzZgxw7TP5Zdfrk2dOtX02K9+9SsNgLZv3z5N0zTt7bff1rxer7Zjxw5N0zStsbFR8/l82ssvv5zyOZ1yyinaOeecY3rMyRo//fRTDYC2fPlyfZ9jjjlG++EPf5h0jW9+85vaWWedlXINhJDcweQvIcQ1xx9/PP785z+bHhs8eLDp+2nTpulfV1ZWoqqqCrt27QIAfPDBB1i+fDkGDBiQdO4NGzZgv/32A4CkaMq6detw+OGHmx474ogjkr4/6KCDcP/99+OKK67Agw8+iLFjx+LLX/5yyufT1dWFsrIy02NO1njAAQfgqKOOwj333IPjjjsO69evx2uvvYbrrrsu6Zjy8nJ0dnamXAMhJHdQ3BBCXFNZWYlJkyal3cfv95u+VxQFqqoCiPtS5s+fj5tvvjnpuBEjRpiukw3nnXceFi9ejCuuuAL33nsvzj77bCiKknL/oUOHYt++fabHnK7x3HPPxU9/+lMsXrwY9957LyZOnIhjjz026Zi9e/di4sSJWT0fQog76LkhhPQ6hx56KD755BOMGzcOkyZNMn2kEzT7778/3nvvPdNjoiJJ5oc//CG2bNmCP/zhD1izZg0WLlyYdj0zZszAmjVrslrj9773PXg8Hjz00EP4+9//jnPOOcdWSH388ceYMWNG2nUQQnIDxQ0hxDWhUAgNDQ2mD7nSKRMXXXQR9u7diwULFuDdd9/Fhg0b8Nxzz+Hss89GLBZLedyFF16ItWvX4vLLL8dnn32GRx99VK9WkgXFoEGD8K1vfQu/+MUvMHfu3KRSbSvz5s3DJ598YoreOF3jgAEDcNppp+HKK6/Ezp07cdZZZyWdf/Pmzdi+fTvmzJnj8CdECOkJFDeEENcsXboUI0aMMH186Utfcnz8yJEj8cYbbyAWi2Hu3Lk4+OCDcckll6CmpgYeT+p/S+PHj8fjjz+OJ554AtOmTcOf//xnvVoqGAya9j333HMRDodxzjnnZFzPwQcfjEMPPRSPPvpoVms899xzsW/fPsybNw8jR45MOv8//vEPzJ07F2PHjs24FkJIz2ETP0JISXP99ddjyZIl2LZtm+nxBx54AJdeeil27NjhqPnfM888g1/84hf4+OOP0wost4TDYUyePBkPPfQQjj766JydlxCSGhqKCSElxZ/+9CccfvjhGDJkCN544w3ceuutuPjii/XtnZ2d2LlzJ2666SZceOGFjrsan3zyyfj888+xfft21NfX52y9W7duxS9/+UsKG0J6EUZuCCElxaWXXopHHnkEe/fuxZgxY3DGGWfgyiuv1Mca/OY3v8H111+PL3/5y/j3v/9tW8pNCOnbUNwQQgghpE9BQzEhhBBC+hQUN4QQQgjpU1DcEEIIIaRPQXFDCCGEkD4FxQ0hhBBC+hQUN4QQQgjpU1DcEEIIIaRPQXFDCCGEkD7F/w/fZPH1Mn+3lgAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# Intensity in frequency domain\n", "srwlpy.SetRepresElecField(wfr._srwl_wf, 'frequency')\n", "data = get_intensity_on_axis(wfr)\n", "fig, ax = plt.subplots()\n", "ax.plot(data[:,0], data[:,1])\n", "ax.set_xlabel(\"Energy (eV)\")\n", "ax.set_ylabel(\"Intensity (a.u.)\")\n", "\n", "fname = \"spectrum.txt\"\n", "np.savetxt(fname, data, header = \"energy(eV) intensity\")" ] }, { "cell_type": "code", "execution_count": 15, "id": "137481bc-7ca8-4bfb-9143-1c6f7ee0dbfc", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjcAAAGwCAYAAABVdURTAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/SrBM8AAAACXBIWXMAAA9hAAAPYQGoP6dpAAB1jElEQVR4nO2deZxU5Zn9z62t927WbrYGREREZAkoQeNEExSXIWMWx5/jQkg00YiJoknUMTKaUYxGozOijiYumYkR98QlbihqIoqCuAEqCjQ2dLP2vtRy7++Pqvfe9966VXVvde2c7+fDp+nqWt4use/p85zneRRN0zQQQgghhJQInnwfgBBCCCEkk1DcEEIIIaSkoLghhBBCSElBcUMIIYSQkoLihhBCCCElBcUNIYQQQkoKihtCCCGElBS+fB8g16iqih07dqCmpgaKouT7OIQQQghxgKZp6OzsxKhRo+DxJPdmDjhxs2PHDjQ2Nub7GIQQQghJg+3bt2PMmDFJ73PAiZuamhoA0TentrY2z6chhBBCiBM6OjrQ2NioX8eTccCJG1GKqq2tpbghhBBCigwnkRIGigkhhBBSUlDcEEIIIaSkoLghhBBCSElBcUMIIYSQkoLihhBCCCElBcUNIYQQQkoKihtCCCGElBQUN4QQQggpKShuCCGEEFJSUNwQQgghpKTIq7h5/fXXsWDBAowaNQqKouCpp55K+ZhVq1bhK1/5CsrKyjBx4kQ88MADWT8nIYQQQoqHvIqb7u5uTJ8+HcuXL3d0/y1btuDUU0/F8ccfj/Xr1+OSSy7BeeedhxdeeCHLJyWEEEJIsZDXxZknn3wyTj75ZMf3v/vuu3HQQQfhlltuAQAcdthh+Pvf/47f/e53mD9/fraOSQixQVU1BCMqyv3efB+FEEJMFFXmZvXq1Zg3b57ptvnz52P16tUJH9Pf34+Ojg7TH0LIwLnwT2vx1WUr0d4TyvdRCCHERFGJm5aWFjQ0NJhua2hoQEdHB3p7e20fs2zZMtTV1el/Ghsbc3FUQkqe9dvb0NYTwrZ93fk+CiGEmCgqcZMOV155Jdrb2/U/27dvz/eRCCkJIqr4qOX3IIQQYiGvmRu3jBgxAq2trabbWltbUVtbi4qKCtvHlJWVoaysLBfHI+SAQtU000dCCCkUisq5mTt3LlauXGm67aWXXsLcuXPzdCJCDlyEYxOOUNwQQgqLvIqbrq4urF+/HuvXrwcQbfVev349mpqaAERLSueee65+/wsuuABffPEFfvGLX2DTpk2488478cgjj+DSSy/Nx/EJOaARjk2Ezg0hpMDIq7h59913MXPmTMycORMAsGTJEsycORPXXHMNAGDnzp260AGAgw46CM8++yxeeuklTJ8+Hbfccgt+//vfsw2ckDygxpwbVc3zQQghxEJeMzfHHXcctCS/9dlNHz7uuOPw3nvvZfFUhBAnROjcEEIKlKLK3BBCCgfRJKWyW4oQUmBQ3BBC0kKImjDFDSGkwKC4IYSkhV6WorghhBQYFDeEENdomgYRteGcG0JIoUFxQwhxjWzW0LkhhBQaFDeEENfIgobihhBSaFDcEEJcI5eiKG4IIYUGxQ0hxDUmccPMDSGkwKC4IYS4RnZrOOeGEFJoUNwQQlwjr1zgnBtCSKFBcUMIcY1cimIrOCGk0KC4IYS4hoFiQkghQ3FDCHGNylZwQkgBQ3FDCHFNhM4NIaSAobghhLjGNKGYmRtCSIFBcUMIcY3KVnBCSAFDcUMIcY15/UIeD0IIITZQ3BBCXGPulqK6IYQUFhQ3hBDXcP0CIaSQobghhLhGLkWxLEUIKTQobgghrjHtlqJzQwgpMChuCCGukQVNOEJxQwgpLChuCCGuUblbihBSwFDcEEJcE8ny+oWIquHp93dgR1tvxp+bEFL6UNwQQlyT7QnFj63djov//B6Ou3lVxp+bEFL6UNwQQlxjagXPQubmna37AQBBtmIRQtKA4oYQ4hpTWSoLzs2QqkDGn5MQcuBAcUMIcU22d0vJ4iZM94YQ4hKKG0KIa7KduRlU4df/vq87mPHnJ4SUNhQ3hBDXyIImnAXnRlGMv+/u6s/48xNCShuKG0KIa7JdlpKfck8XnRtCiDsobgghrjFvBc+GuDGec08nnRtCiDsobgghrsn2bilZL7EsRQhxC8UNIcQ1apYzN3Kpi84NIcQtFDeEENfI3dlZL0vRuSGEuITihhDimmwvzmSgmBAyEChuCCGuyXagWKNzQwgZABQ3hBDXmALFWRggzLIUIWQgUNwQQlwjmzXhLKgb+fn3dge5goEQ4gqKG0KIa1TT4swsPL/k3GgasK+HuRtCiHMobgghrpHXL2RlQrHlOfd0UtwQQpxDcUMIcY2cuclOK7j5c+ZuCCFuoLghhLhGy+H6BYDihhDiDoobQohrTM5NlufcABQ3hBB3UNwQQlwjh4izkbnRLIJpLwf5EUJcQHFDCHGNqSyVFefG/JyhbLRkEUJKFoobQohr5LJUOAvCw2oGZWPFAyGkdKG4IYS4JpL13VJa0s8JISQZFDeEENeo2W4FtzxnNl6DEFK6UNwQQlwja41sbwXP1msQQkoXihtCiGtMmZsszrlRlPjXI4SQVFDcEEJco2Z5iJ94ep9Hib1Gxl+CEFLCUNwQQlyjZnu3VOz5fZ7ojyjr3BtCCEkGxQ0hxDWyk5LNOTc+r5K11yCElC4UN4QQ12S7LCWe0u/1ZO01CCGlC8UNIcQ1uWoFF5kbdksRQtxAcUMIcY15iF/mMzFCzLhxboJhFRt2dDCfQwjJv7hZvnw5xo8fj/LycsyZMwdr1qxJev/bbrsNhx56KCoqKtDY2IhLL70UfX19OTotIQSIDxFn2rwRz6dnbhx0S/3m+U045b/ewMqNuzJ7GEJI0ZFXcbNixQosWbIES5cuxbp16zB9+nTMnz8fu3bZ/3B66KGHcMUVV2Dp0qXYuHEj/vCHP2DFihW46qqrcnxyQg5srAHfsJrZXm2jWyoqbpy4MV/u7wEA7GjvzehZCCHFR17Fza233orzzz8fixYtwpQpU3D33XejsrIS9913n+3933zzTRxzzDH4t3/7N4wfPx4nnngizjzzzKRuT39/Pzo6Okx/CCEDI26CcIbn0GjWQLEDcSPOlI3WdEJIcZE3cRMMBrF27VrMmzfPOIzHg3nz5mH16tW2jzn66KOxdu1aXcx88cUXeO6553DKKackfJ1ly5ahrq5O/9PY2JjZb4SQA5C43U9ZytwYZSkH4iZ2nywsKSeEFBm+fL3wnj17EIlE0NDQYLq9oaEBmzZtsn3Mv/3bv2HPnj342te+Bk3TEA6HccEFFyQtS1155ZVYsmSJ/nlHRwcFDiEDxCo2Mt0xpWduYkP8nHRLifswUEwIyXug2A2rVq3CDTfcgDvvvBPr1q3DE088gWeffRa//vWvEz6mrKwMtbW1pj+EkIFh1TKZFzeiW8q5cyMcG7aNE0Ly5twMGzYMXq8Xra2tpttbW1sxYsQI28f86le/wjnnnIPzzjsPAHDEEUegu7sbP/rRj/Dv//7v8HiKSqsRUrRYBUTGxY0+50Y4N6kfIxwbRm4IIXlTA4FAALNmzcLKlSv121RVxcqVKzF37lzbx/T09MQJGK/XC4BWNCG5xCpmMu2WWDM3TkLCqi5u+LOAkAOdvDk3ALBkyRIsXLgQs2fPxlFHHYXbbrsN3d3dWLRoEQDg3HPPxejRo7Fs2TIAwIIFC3Drrbdi5syZmDNnDjZv3oxf/epXWLBggS5yCCHZxxogzlbmxk23lDgDtQ0hJK/i5owzzsDu3btxzTXXoKWlBTNmzMDzzz+vh4ybmppMTs3VV18NRVFw9dVXo7m5GcOHD8eCBQtw/fXX5+tbIOSAxOqUZlrcaJY5N86cm+ychRBSfORV3ADA4sWLsXjxYtuvrVq1yvS5z+fD0qVLsXTp0hycjBCSiFx1S7mac6OyLEUIicIELiHENdZ1CNmac+P1OF+/oDJQTAiJQXFDCHGNtSyV6anA1t1SThoGRCs4mwsIIRQ3hBDXxAWKM+3cxNSN3+N8K7jGbilCSAyKG0KIa6xiI5zhnQf6ED9frCzlolvKSQmLEFLaUNwQQlxjdUeyNudGDPFz0S3FshQhhOKGEOIa6xbw7HVLOXdu2C1FCBFQ3BBCXGMVG5kWFJreLSWcm9SPYbcUIURAcUMIcY21TJT5zE30o3BunIinCAPFhJAYFDeEENfELc7McubGWbeU+SMh5MCF4oYQ4hqrUeOkbOQGfSu4G+dG75aiuiHkQIfihhDiGmtZKvPOTfSjHijmVnBCiAsobgghronfLZVZ6yadspTRLZXRoxBCihCKG0KIa+IyN5kuS1mcGydmDOfcEEIEFDeEENfEi5vstIL7XGwFZ7cUIURAcUMIcY1VzGRvQrHzzI3GOTeEkBgUN4QQ1wgtI8RHOGsTimND/Fx0S2V6QzkhpPiguCGEuEaUgHTxkXFxY24Fd9YtZX4sIeTAheKGEOIaITbctGq7QZ9z4xFzbpw/ptSMm67+MH7+6Pv4+2d78n0UQooGihtCiGuEORLwOQ/8usFalgJSu0OlOufmxr9txKNrv8TZf3g730chpGiguCGEuMZwbpzPoXGDaumWAlILKKMVPKNHyTvrt7fl+wiEFB0UN4QQ11gzN5lvBY9+9MfKUk5eQ5yp1NYv7OkM5vsIhBQdFDeEENeolsxN1lrB5bJUitfQSrQstaerP99HIKTooLghhLhGzbJzI57f68a5KdFAcabb7Ak5EKC4IYS4RggJPVCctTk3StxtqR7D9QuEEIobQohrrN1M2WsFd9YtJX+tlMpSslCrKffl8SSEFBcUN4QQ1xhlqdicmyxlbmTnJtlryIKmlKo4+3tC+t+HVgXyeBJCiguKG0KIa6yt4JmfUBz96PEoELGbZK8R0UrTuWnt6NP/7pHyR4SQ5FDcEEJcIwREICZuMr9bKvp8HkXRQ8XJnBv5S6UqbkqtxZ2QbEJxQwhxjbjQit1PmXZuhD7xKICipF7BIF/4VTWjR8kruzqMNnCKG0KcQ3FDCHFNXKA4S5kbj6LAq6QWUOoBUJbitnNCnENxQwhxhXyRNVrBM/waMYGiKMasm2TOhezWlJC2QWunVJYqpW+MkCxDcUMIcYV8kfV7hLjJrLrRA8WKESh23i1VOiKglWUpQtKC4oYQ4gpZPPh9wlXJ8Guo8YFip91SpeRw7GKgmJC0oLghhLhCNmn0VvCsZW6iAgc4MOfc0LkhJD0obgghrpBFRiBru6WiHz0eRZ/vkqzyZc7clIYIiKgadktLM6ltCHEOxQ0hxBWykMlGt5QsTkzdUgdY5mZvd7/pvQ6XUo87IVmG4oYQ4gpZfOjiJpI5QSE7FB6H3VKlOOemozdk+rxUvi9CcgHFDSHEFSbnxpf53VKy86IoCsTuzANtQrF16nMpBaUJyTYUN4QQV0SkGTQ+B51MblFNZSk4GuJXirulwhGxGT21c0UIMUNxQwhxhbEaQXHUyeQWufwiv0aya3spdksJ50YMSgQ4pZgQp1DcEEJcIRwErzSDJpOLM1VLoNjjwLnQStC5EYMRZXGT6QWlhJQqFDeEEFcIkeHx5KAs5YGjbqlICa5fEGWpMtm5KZVvjpAsQ3FDCHGFvNTSiavi/vmNvzt9jVJsBRffb5nPG3cbISQ5FDeEEFeI66vX4Qwat8TNuXHQLSVf9EtFANhlbtgxRYgzKG4IIa4wylKGq5LZzI3xd6fdUvI1v1Su/2Jon5gCDWR2nhAhpYzP7QP6+/vx9ttvY9u2bejp6cHw4cMxc+ZMHHTQQdk4HyGkwJD3PmWjTdk650Zx0C1Vyq3gdG4IcY9jcfOPf/wDt99+O55++mmEQiHU1dWhoqIC+/btQ39/PyZMmIAf/ehHuOCCC1BTU5PNMxNC8ogQD16PtLE7C0P8Yk/taEJxKWdu/F4FihJ1pNgKTogzHJWlvvWtb+GMM87A+PHj8eKLL6KzsxN79+7Fl19+iZ6eHnz22We4+uqrsXLlSkyaNAkvvfRSts9NCMkTellKnnOTSecm1vkkntvRbil5/UKJXP9Fqc/rUQyHrESEGyHZxpFzc+qpp+Lxxx+H3++3/fqECRMwYcIELFy4EBs2bMDOnTszekhCSOEgxIfpopuFspQQN/r6haTOjfH3UtoKDkT3d0XfC00vVRFCkuNI3Pz4xz92/IRTpkzBlClT0j4QIaSwyX4ruDFHB4Cj0lcpd0tlq/xHSCnDbilCiCsikvjw6usXMvf88noH+WOyC7t5QnHmzpJPwrHJhD6P0XJfKsKNkGyTMXGzcOFCfOMb38jU0xFCChSRb5HXL2RjQrFV3MhTiK2UZLeU7Nx46dwQ4gbXreCJGD16NDweGkGElDpyoDg7u6WiHxVLt1QyAWXO3GTsKHlFvM8+j0d3brhbihBnZEzc3HDDDZl6KkJIASOurx5Pjp0bx91SpSEAZOcmG9kmQkoZWi2EEFfoc27kVvBMzrnRnaHo515H3VKlFygWW8F9XmnNRZLSHCHEwLVz84Mf/CDp1++77760D0MIKXyEeFAUZyUjt+jOkJhz47JbqkSMG9258UkOGefcEOIM187N/v37TX927dqFV155BU888QTa2tpcH2D58uUYP348ysvLMWfOHKxZsybp/dva2nDRRRdh5MiRKCsrw6RJk/Dcc8+5fl1CSHrYTSjObObG2F0FwFi/4DBzUyplKbFHyuvxOJrSTAgxcO3cPPnkk3G3qaqKCy+8EAcffLCr51qxYgWWLFmCu+++G3PmzMFtt92G+fPn45NPPkF9fX3c/YPBIE444QTU19fjsccew+jRo7Ft2zYMGjTI7bdBCEkTO3GTnSF+0c+dtJuX4vqFkJ1zQ3FDiCMyEij2eDxYsmQJjjvuOPziF79w/Lhbb70V559/PhYtWgQAuPvuu/Hss8/ivvvuwxVXXBF3//vuuw/79u3Dm2++qU9LHj9+fNLX6O/vR39/v/55R0eH4/MRQuKJSOsRnKxGcIt1zo2zbqnSm3MjMjdej6ILPYobQpyRsUDx559/jnA47Pj+wWAQa9euxbx584zDeDyYN28eVq9ebfuYv/71r5g7dy4uuugiNDQ0YOrUqbjhhhsQiUQSvs6yZctQV1en/2lsbHT+TRFC4ohIgV8nqxHckk63lPX1S2EFg5y58cXe6FJxpQjJNq6dmyVLlpg+1zQNO3fuxLPPPouFCxc6fp49e/YgEomgoaHBdHtDQwM2bdpk+5gvvvgCr7zyCs466yw899xz2Lx5M37yk58gFAph6dKlto+58sorTWfu6OigwCFkAGhSWUpcdDMrbqIfFRfdUtZrfkTV4IsNvitW9MyNl63ghLjFtbh57733TJ97PB4MHz4ct9xyS8pOqoGiqirq6+txzz33wOv1YtasWWhubsbNN9+cUNyUlZWhrKwsq+ci5EAiIjkruvDIZCt4AucmmRtjveiXggYQzo3f43Ek8AghBq7FzauvvpqRFx42bBi8Xi9aW1tNt7e2tmLEiBG2jxk5ciT8fj+8Xq9+22GHHYaWlhYEg0EEAoGMnI0Qkhh5QrEnCzuPrHNuDNciyWM0q7gpfhEg3lMvd0sR4pq8DfELBAKYNWsWVq5cqd+mqipWrlyJuXPn2j7mmGOOwebNm6FKk6w+/fRTjBw5ksKGkBwhdIM3axOKox/1QLGDzI31SyWgbRBWpcWZnHNDiCsyJm6uuuoq12WpJUuW4N5778WDDz6IjRs34sILL0R3d7fePXXuuefiyiuv1O9/4YUXYt++ffjZz36GTz/9FM8++yxuuOEGXHTRRZn6NgghKdCdmxzNuXEioKwX/VJwbsJS5iYbIpKQUiZju6Wam5uxfft2V48544wzsHv3blxzzTVoaWnBjBkz8Pzzz+sh46amJtMyzsbGRrzwwgu49NJLMW3aNIwePRo/+9nP8Mtf/jJT3wYhJAURaQ6Nk+nBbrHOuXG0W6oExU1E6pbycHEmIa7ImLh58MEH03rc4sWLsXjxYtuvrVq1Ku62uXPn4q233krrtQghA0e4B14lO1mQ+Dk35tdNdibj84wdJ28YizM9WRGRhJQyXJxJCHGFvBXco190MzdbRlzAFUu3VLILu1X3lIIIEILR7+WEYkLckpZz093djddeew1NTU0IBoOmr/30pz/NyMEIIYWJXJbyeYxZMpmaLWMEimMfHXRLxbeCF78ICEsTiiluCHFHWnNuTjnlFPT09KC7uxtDhgzBnj17UFlZifr6eoobQkocVW5RlsRNWNXg8yZ6lPvnt3ZLJXduSnDOTUTaLcVWcEJc4bosdemll2LBggXYv38/Kioq8NZbb2Hbtm2YNWsWfvvb32bjjISQAkIesuf3Gj9CMnXhjQsUO3AtrOKmlNYveD0e4z0oge+LkFzgWtysX78el112GTweD7xeL/r7+9HY2IibbroJV111VTbOSAgpICJJnJtMIGd6AGdzbqwvXQoiQO6W8rEVnBBXuBY3fr9fb8+ur69HU1MTAKCurs51KzghpPiQnRshPIBsODciUBy9/cBbvyBtBWfmhhBXuM7czJw5E++88w4OOeQQfP3rX8c111yDPXv24H//938xderUbJyREFJAiGCvR4ledD1KVEyEM9R/raVRlrIKn1JwOEzdUpxzQ4grXDs3N9xwA0aOHAkAuP766zF48GBceOGF2L17N+65556MH5AQUlgIZ0XEbTK9GdzYCm6eUJy8W8r8eQlUpTjnhpAB4Nq5mT17tv73+vp6PP/88xk9ECGksInrZvIoQMTo7hnw81ucm/S6pYzPn/twJyY1VGNifU1GzpcrTN1SDgQeIcSAQ/wIIa7Q59zELri+DOdBrIsz0+mWEp9vaunAT/60Dpc9+kFGzpZLwnJw24HAI4QYOBI3J510kqOVB52dnfjNb36D5cuXD/hghJDCRGgMccH1ejObB4mfcxO9PZ3dUnu7okNG23uCcY8pdCLSVnAGiglxh6Oy1Omnn47vfve7qKurw4IFCzB79myMGjUK5eXl2L9/PzZs2IC///3veO6553Dqqafi5ptvzva5CSF5whAf0c8z79yI9QvRz8WFPXm3lPU5oh+DsS8UY2u4ybnxmG8jhCTHkbj54Q9/iLPPPhuPPvooVqxYgXvuuQft7e0AoqG/KVOmYP78+XjnnXdw2GGHZfXAhJD8Yi1LiTxIprqldGfIY94t5apbKvZ5KBw9UzEu0tTn3Hg9RlmK4oYQRzgOFJeVleHss8/G2WefDQBob29Hb28vhg4dCr/fn7UDEkIKC3krOJCNbimbwDJc7paK3Vc4HcVYzjEHimPvcRE6UITkg7QWZwLRoX11dXWZPAshpAgwWsGtzk1mLrzWOTfOuqXszxgq4rJUxKYsReeGEGewW4oQ4grhoCiKuVsqc63g5udPb7dU9GMwVpYqxl1Toszn8xqBYmZuCHEGxQ0hxBXWIX6Zz9xYJhQr5tuTPUYgnJqiLktJu6W4FZwQd1DcEEJcoYubuExMllrBHUznTdQKrpelLGfbuqcbz3ywo6AdnUjEmFDsc/AeEEIMKG4IIa4QQkEvS2V6zo11iJ8D1yJ+/UL0vqIsZX3oFU98gMUPvYf129sGfuAsITs3nHNDiDtci5uFCxfi9ddfz8ZZCCFFgDVQrHdLZXr9giWwnKzqFd8KHv0Yip3J6ni0tPcBAPZ1F+5wP6MVnGUpQtziWty0t7dj3rx5OOSQQ3DDDTegubk5G+cihBQoQmQY4iZbzo34GLuwJx3iZ78VPJygLNXVHwFglK0KEZFh8tK5IcQ1rsXNU089hebmZlx44YVYsWIFxo8fj5NPPhmPPfYYQqFQNs5ICCkghMgQE4QznbnREs65cd4KHrFkbqzOTXd/GAAQzJDblGlUVdO/J5+UuSnGlnZC8kFamZvhw4djyZIleP/99/H2229j4sSJOOecczBq1Chceuml+OyzzzJ9TkJIgRA3xM+bnW4pxdItlSz8m7AVPBLfLRVRNfSGYs5NuDCdG9kF80pbwTnnhhBnDChQvHPnTrz00kt46aWX4PV6ccopp+DDDz/ElClT8Lvf/S5TZySEFBARS+bGm/EJxdGPcVvBB9AtJR+tOxjW/16oZSn5vfR5FP294JwbQpzhWtyEQiE8/vjj+Od//meMGzcOjz76KC655BLs2LEDDz74IF5++WU88sgjuO6667JxXkJInrGKj8xnbuwnFLtavxD7NCw9SLgeXX2FL25kF8znVRy1wxNCDFyvXxg5ciRUVcWZZ56JNWvWYMaMGXH3Of744zFo0KAMHI8QUmhYt4LnbM5N0sWZlucQreBSpiaiafBA0fM21q8XEmbnxpPx95iQUse1uPnd736H008/HeXl5QnvM2jQIGzZsmVAByOEFCbyziMge91SinXOjYtuKc1SlhL38XuBrv5icG6M78ejOFseSggxcF2WevXVV227orq7u/GDH/wgI4cihBQuiebQRDJ05U203iGZcxO3fkFsBZfOJO7SHWsDBwo4UCxtBFcUxdHyUEKIgWtx8+CDD6K3tzfu9t7eXvzxj3/MyKEIIYWLamnVzt6cG+HcmF832Zmsn4csZSmgWJwbY8YNAC7OJMQljstSHR0d0DQNmqahs7PTVJaKRCJ47rnnUF9fn5VDEkIKh6DkKgCZ75ayzrlx1i1l/xxBS1kKQFFlbnyW0h9bwQlxhmNxM2jQIChK1CKdNGlS3NcVRcG1116b0cMRQgqPvmC0rFMZiP74yFa3lGLplko2RidRt1TIrluqKJwbsXohKhw5oZgQdzgWN6+++io0TcM3vvENPP744xgyZIj+tUAggHHjxmHUqFFZOSQhpHAQA/AqAtELr1cM8UvTBflidxfWbtuP735lDDweJa4s5WxCsX1ZSj6TWkRlKatz43UQqiaEGDgWN1//+tcBAFu2bMHYsWP1TgZCyIFFT2wIXoXf7NxE0pxQvPSvH+ONz/ZgzOBKzD14aFyruZNuqXhxE/1oKktp8WWpQhU3QpQZgxKjt9O5IcQZjsTNBx98gKlTp8Lj8aC9vR0ffvhhwvtOmzYtY4cjhBQefaGoIKgIeAEYF+B0y1LtvdHuy7ae6IZua2A5FulJvn5BtX4e3wou7mPK3IQLUyyIQHG2ck2ElDqOxM2MGTPQ0tKC+vp6zJgxA4qi2P6gURQFkUjE5hkIIaWCcG4qY+LG7x3YhVe4FKHY4/WylLUkk+T5ra6OXVnK6JaSWsEL1bkRs4S8ZueGreCEOMORuNmyZQuGDx+u/50QcuCiZ278mXFuhGgRM2es6xechGmtv2wlDxQbc7oKVdwYmZtYoNiBwCOEGDgSN+PGjbP9OyHkwEJVtbiylG+AnTyiBCM+Cp2iB4r1AXaJnyO+Wyq+FVzVMzdF4NzEZW4obghxQ1pD/J599ln981/84hcYNGgQjj76aGzbti2jhyOEFBZ9YUMYxDs36QkFw2URZSnRCu6mW8r8eaL1C4C5W6pY5txQ3BDiDtfi5oYbbkBFRQUAYPXq1bjjjjtw0003YdiwYbj00kszfkBCSOHQE4wXN5lyboQQSViWctAtZZwl9tw2reCmbqlCXb8gAsVeS+6ImRtCHOF6ceb27dsxceJEAMBTTz2F733ve/jRj36EY445Bscdd1ymz0cIKSB6Y+KmzOeRdktFf0dKd85NJGIO/wphYl2/kLRbSjPKOGFVk9YvyGWp6MdiaAU3lpPGZglxQjEhrnDt3FRXV2Pv3r0AgBdffBEnnHACAKC8vNx25xQhpHQQYWLRKQVkwrkR3VIic2N2bpx0S4mKmOjc0tcvhJOXpQpV3IQsKy64W4oQd7h2bk444QScd955mDlzJj799FOccsopAICPP/4Y48ePz/T5CCEFhHBuREkKGHi3lLVt25q5ERd2VYuKFrsBonpZymsOH5sWZ6rR3XjFlLkR7y13SxHiDtfOzfLlyzF37lzs3r0bjz/+OIYOHQoAWLt2Lc4888yMH5AQUjgI56Zcdm68GXJu9MxN9HavZc6N/DUrRubGY/pcDjmrWrTTS36OQnVurEP8nExpJoQYuHZuBg0ahDvuuCPudi7NJKT06Q3Gl6UG2i0lMjfWbilroBiICiivJ965sXYXCZcnZAoUm0tS0dcsTHGjfz9ec+amQI9LSMHhWtwAQFtbG9asWYNdu3ZBlX6gKYqCc845J2OHI4QUFtYBfkDmMjfhiP2cG1nLJJrQK27Wy1KqWdiI83VbxU3BdkvZt4JzQjEhznAtbp5++mmcddZZ6OrqQm1tran+TXFDSGkjWsErAsaPDtHRYxUTTolo1rKU/Zwb+WuJnsMniQCrk6Ra8jZA8WRudHeM1g0hjnCdubnsssvwgx/8AF1dXWhra8P+/fv1P/v27cvGGQkhBYLh3Bg/Ogbq3OjrF6y7pSxbwZO9hhEoFpkbIBRO7NyIpyzUspQQMT5L7oh5YkKc4VrcNDc346c//SkqKyuzcR5CSAHTqy/NlJ2b9DM3mqbpgkVc0EVHkMfOuUnwEuJ22bkJRhI7N4Mq/AAKWNwkcG44oZgQZ7gWN/Pnz8e7776bjbMQQgqc3mBUDJRnKHMjPyZRoFjulkrULSQeI8+5sQoXVTUCxYMrAwDSHzyYbawBaSdTmgkhBq4zN6eeeip+/vOfY8OGDTjiiCPg9/tNX//Wt76VscMRQgqLnpBwbuRW8NiE4nTEjSaLG8v6BcuFHUgsoKwZFVWLFy4RTdOXZg6qjP7cCkbUhLNz8knY2i3FreCEuMK1uDn//PMBANddd13c1xRFQSQSibudEFIa9NkM8cuUc2MM8Yt+LmdtPIrR3m2HeIxfmrkTV5aSMjfCuQGiQkI8rlDg4kxCBoZrcaOmOcuCEFL86IFiuzk3aZR4wmq8c2NdvyBeQ41oKctS4ix2ZamIamRuBlcZ4iYUUfVyVqEg3ktr5gaIijTZzSKExDOg/6P7+voydQ5CSBHQk2nnRhJE1m4pxeTcJH8Na+ZG1eLDwqomOzdGOd3aVVUIWCcUO8kdEUIMXIubSCSCX//61xg9ejSqq6vxxRdfAAB+9atf4Q9/+EPGD0gIKRz6kjk3abi68oU6bM3cSBd0Yyt2gueJm1AcP8RP1TR0x7q9asr9ujNkLV8VAka3VPRHtEf6Sc3SFCGpcS1urr/+ejzwwAO46aabEAgY1u7UqVPx+9//PqOHI4QUFj026xcGslvKLnMjbjOVpVLsVjImFCd2biIqEIy5NGU+j+7yFGI7uLF+QSzO9MR9jRCSGNfi5o9//CPuuecenHXWWfB6jR9w06dPx6ZNmzJ6OEJIYaEvzjRtBU+/W0p+TDDB+gXAGLqXcEKxrXNjETeahkjM+vF6FAQKWNwIoWe0ghtfY1mKkNSkNcRv4sSJcberqopQKJTWIZYvX47x48ejvLwcc+bMwZo1axw97uGHH4aiKDjttNPSel1CiDvsFmdmKnMjylrG+gXjfkZZyuGEYlWLCzhrmmYajidckUIUN5EkmZtE7wEhxMC1uJkyZQreeOONuNsfe+wxzJw50/UBVqxYgSVLlmDp0qVYt24dpk+fjvnz52PXrl1JH7d161ZcfvnlOPbYY12/JiEkPewWZxqZm4HNuQlbhvjJHULeFEPsdHEjzbmxZmkiqmZyeERZKliQgWJz5kZ+L9J5nwk50HDdCn7NNddg4cKFaG5uhqqqeOKJJ/DJJ5/gj3/8I5555hnXB7j11ltx/vnnY9GiRQCAu+++G88++yzuu+8+XHHFFbaPiUQiOOuss3DttdfijTfeQFtbm+vXJYS4x1icmak5N4YACeqB4ujnHlfdUuaz2JalJHHj9RR25kYvS8XcJUVR9Fk/dG4ISY1r5+Zf/uVf8PTTT+Pll19GVVUVrrnmGmzcuBFPP/00TjjhBFfPFQwGsXbtWsybN884kMeDefPmYfXq1Qkfd91116G+vh4//OEPU75Gf38/Ojo6TH8IIelhN8RPuArpiISwTaA40ZwbwEG3VLL1C5rZuQn4CljcqO7dK0KIgWvnBgCOPfZYvPTSSwN+8T179iASiaChocF0e0NDQ8Jw8t///nf84Q9/wPr16x29xrJly3DttdcO9KiEEBhlKXlxpujkSce5kXMxYYtzYzvnxubCLk8tlstS8a3ghmjweBR9KnEhtoJbMzeAeA80dksR4gDXzs2ECROwd+/euNvb2towYcKEjBwqEZ2dnTjnnHNw7733YtiwYY4ec+WVV6K9vV3/s3379qyekZBSJRhWdXFgcm68A8jcmLqlrIszpQt77CeVXbeU/LJ6W3qKspScubGKoEIgLJ1TkMq9IoQYuHZutm7dars/qr+/H83Nza6ea9iwYfB6vWhtbTXd3traihEjRsTd//PPP8fWrVuxYMEC/TaxDsLn8+GTTz7BwQcfbHpMWVkZysrKXJ2LEBKPcG2ADGZu5ECx6JZKMufGLm8iv644i6ZpCIXjy1JhqRVcFzfhwlMLejZIWgsxkGGJhBxoOBY3f/3rX/W/v/DCC6irq9M/j0QiWLlyJcaPH+/qxQOBAGbNmoWVK1fq7dyqqmLlypVYvHhx3P0nT56MDz/80HTb1Vdfjc7OTtx+++1obGx09fqEEOeINnCvVNIRnwPRC7LbDduOF2cmEVCym2O0gsc7MhFV010PX6HPuUnm3DBzQ0hKHIsbIT4URcHChQtNX/P7/Rg/fjxuueUW1wdYsmQJFi5ciNmzZ+Ooo47Cbbfdhu7ubr176txzz8Xo0aOxbNkylJeXY+rUqabHDxo0CADibieEZBa5DVwWMH7L9Fyfiw3bcuYmZF2/IF/Yk2Ru5Iu9X+6WUq3OjeF6eDwK/L7CzdyI/JHX7j0ovOMSUnA4Fjei/HPQQQfhnXfecZx5ScUZZ5yB3bt345prrkFLSwtmzJiB559/Xg8ZNzU1weNxHQ0ihGSYntheJrkkBRiZGyDqOPjMX06K7MSE4lrBpddIkjcxZ26k9QuW+TVqkWdukrlXhBAzrjM3W7ZsyfghFi9ebFuGAoBVq1YlfewDDzyQ8fMQQuLpsxngB5gvwG4vvLITI2a4aDaBYuEU2ZVk5Nf0ypkbm/ULcot1Ic+5iUjnFAwk20TIgUZareArV67EypUrsWvXLt3REdx3330ZORghpLCwW5oJDGx6bsTy8yOkqgnWL8Tun6IV3C93S6nJuqU8RZG58UuB4mTt8IQQM67FzbXXXovrrrsOs2fPxsiRI12FBwkhxYsIFJdbnBt575FbV8G6/ykU0WwDxU67pcS6AruylCYN8fN4DCEULORuKbshfnRuCEmJa3Fz991344EHHsA555yTjfMQQgoUY4CfWdx4PMZqALdtytYLdTiiSq3gTruloh8VBRDxH/v1CzA5NwPN3PQGI3H5o0zBbilCBobrpG4wGMTRRx+djbMQQgqYXpvVC4J0pxRbSyxR5ybJnJsk3VIeRdFFkKPMzQDWLzyx7ktMWfo8Hn03O0NBI6pNt5SYc1OAAWhCCg3X4ua8887DQw89lI2zEEIKGL0V3MatSPfCaxVDoYhqv35Bd27in0PfIq4outtjN+fG2i2VbuZG0zQseeR9aBpw+8rPXD3WKfriTKlTNJnAI4SYcV2W6uvrwz333IOXX34Z06ZNg9/vN3391ltvzdjhCCGFQ09S5ya9PIhVDIUl50Z2LTxSucmKeE1FMYduky3OlAcRup1zs/oLY/3M4aNqXT3WKXaLM9kKTohzXIubDz74ADNmzAAAfPTRR6avMVxMSOnSm6BbCkh/v1Scc6OqEPrFds6NbbeUcR/xGLksVebzoD+sJm4FD7s785/eatL/LnczpcPHO9px9u/fxsHDq/G9WWPwvVlj4PN69PfFPAk6+pHdUoSkxrW4efXVV7NxDkJIgdPeGwIA1FX4474mnBvXgeK4zI0Ku8WZojxj19kUkQLIellKM8pN5X4v+sNqrCxlbNtOZ87Nro4+vPBxi/55/wA7rdY1tWF/TwjvbtuPd7ftR39YxcKjx5t2YAlEJ1gkS5mbnz38HvZ2BfHgD44yvS4hxUhac24IIQcebTFxU2sjbtLN3FidHrksJRvB5f7ohb3PRkzIAWRFKl8JzVLu96C919wt5fUoCKQRKF7XtN905oGKm7DltVs6+gAYAsacuYl+zIZzo6oa/rJ+BwDgy/09GDe0KuOvQUgucSxuvvOd7zi63xNPPJH2YQghhYtwbgZVBuK+lna3lOXiLgeKZedGzNbplzaTC+RdVEb5yhAOZbF9EJnI3IjckSAYjj+PG6xiUEyBtsvcGCsoMi9u5PdgV2c/xQ0pehyLG3kLOCHkwKO9JwjAviylOzduA8VW50bVbOfclMcESp+tuImdQSpLWTM30fslWr/g/MyiY0xRolmfgTo31inKVnEjLyHN5oRi+fvY1dGf8ecnJNc4Fjf3339/Ns9BCClwDOcmcebGrXNjDQiHwqrtnBu9LBVKnLlRFEUvS0VUTRctwvWxrl8wAsXOBYp4/boKP9p6Qui3OY8brPkZEdq22wouhE42uqXkLFNrrDRGSDHDdduEEEe0JQkUe9MMFFudm5AqrV/wxJel7J0bQwx5pFkwRqDYcG4iUpt5OnNuxOsPir0H/QMsS4WkNvbo80fPIs/jEejOTRbEjfx97Oqkc0OKH4obQkhKVFVDh3Bu7LqlvOlmbmzWL9h0S5UJcWMjJoSeiraCx3dLicxNOKLpbeM+jwK/z33mRogbIfDczsixIhya6rKoiS7KXqEkmZtsOze7OunckOKH4oYQkpLO/rDuqNh1S/kylLkJRezn3CQrS5nWL5jm3ERvF5kbWYh40mwFF+JGvAcDLUuJ77+2PPp8vaEIIqqmiw15YGI2JxTLmZvddG5ICUBxQwhJiXBtyv2euK3ggOQquGwFj8vcmHZLGeqmIklZKqJ3SxmlLOucm+hzGxdw3wADxXV6WWqAgWKLc9MfiuivAQCVASMWmWwFxUAJMlBMSgyKG0JIStp6EudtgMw5N2FVTTDnRoib+Cu7ZnJu4jM3Zf74AYDpZ26MQDEw8MyNKDFVlxtlKREqVhTDsQLk0Hbm1Y3sarWyLEVKAIobQkhK9E6pivgZN0D6eZC49QthTc/QmOfcRH9U2YkJcyt47DZVLkvFZuSE7Z0bu6nHibBzbrQBlInEGWtsxE2F35tgeWgWylKSaGzrCQ1YtBGSbyhuCCkiNE3Dix+3oLmtN6ev29abeMYNYLQpu+6WspSEQqp9oDjZnBu7xZnmQHG8QyMP8XPj3PRbxI2mRd2nO175DOf84W1XQgkwAsU1scxNX0hFTygMIH6Hl1efc+PqJRwRjJjfV+ZuSLFDcUNIEbF223786H/X4orHP8jp6+p7pWxm3ADS3qMBzrkJmzI3xu3JylLyFnF5/UJ85kZe06DA73OfubGWpYCoe/PH1dvwxmd7sHFnh+PnAoyynMjc9AUjxvZ1q7jJ4oRiazCa7eCk2KG4IaSI2NEezUO0tOc2F5FsaSYwkMxNkvULkrop07ulEreCxy/OtHRLhcXSzOjn6WRurGUpIBYCjgkSN88l3782VpbqCxvPVek3z1hNdwq0E6wt7bs4yI8UORQ3hBQRXX3RkoV1x1G2ae9JPOMGyFzmxtSubbNbynbOjVTG8kqhW/Hc4rFC3Ij7DKQVvCLg1cVRf1hFT+x2t3NvIhbnJhTR0NkX0l9DJlet4ACdG1L8UNwQUkR09wtxE87p6zp2blxe3K2ZGzmzYipL+RKXpeRWcKGH5Iu1NXNjiBsl7jVTIZybcr9Xf97u/rAuUtyUuOT7i24pANjXHX2vrZmbrAaKreKG7eCkyKG4IaSI6IqJm+4cOzeiFdxurxSQfskkzrmRLrKKTbdUn833LbqV5MWZcvlK3yie0LlxfmaRTanwexGIiRuxlgJwt6cKMMpyVWU+XZjt644Ki7hAceyndbYnFAOcUkyKH4obQooI4dwEw6rrfMdAEM6N3XRiIP3FmdYN1/2JnJskZSnxNiiSuOnpN2bFCJEg3i9xViFO3HR49dk4N/u7g/rX3f43Ec6V36vo7tTe2PNVBMyZG1+aoW0niNZvIfxYliLFDsUNIUVEt1SOymXupk3fCJ5ozo0QCukN8Su3GbRnl7kJRbS4i7vcLSUEkXifKvxevZwTjCtLGc6N01k1RlnKo++7kp0bt5kbQ3B59IyNEEuVlknQ+uLMLGRuxPs+orYcANDKshQpcihuCCkiOvtkcZO73E2Hw8xNuoszjUF7hmCzG+IHxHdMidZoj2JevwDExE3saUJxZSnj+Z2UpjRNMwLFknMjwtZOn0dGvF9R5yb6fIZzY1+WykYruBA3owdXRM/QRXFDihuKG0KKCFGWAnLs3PREL7gJu6XEED+XF3erc2MqS0k/nUTJBrARN/qiTcUkiICoQBBdRomcG8BZOSm69yr69zJT5ib9spSx/duD8piY2ZdQ3KTnjjlBvO/iv69dyz0hxQTFDSFFRHe/cdHp6c/NBSgUUfUAcyLnxp/m3iNRUqqwtGsDZufGI+2C6rOEXyM2W8EFcllKuCo+G3HjZAGmvNBSdm72m5yb9CYU+6TMTaKyVLodaU4Q73tthhaCEpJvKG4IKSK6JOemO0dlqXYpU5IoUDzwzE1ycQMkHuSnmTI3iZ0b46yK/nH0oGgZ5sn3mlOeVaxe8CjRMpIopcllKffrF2JlKSlzk6gspWeEshgoFjuuBrozi5B8Q3FDSBFhDhTnVtzUlPt0YWBF7JZyP8RPbO6OiZuIfbcUIK9gMIsbebeUYvMYj+WnnPw9XPyNiQCA/37lM5NIsUOecaMoiu7cmMtSLufcqIZzU2FpWa+0dkt5c+DclBvi1W04mpBCguKGkCKiSwoUd+eoLJVqxg2Q/pwb4VyIMK2840ixKBV91o1lkF+yzE1lwBt3m1dSO9+bNQaTGqrR1hPC8lWbk561T5pxAxit5G0DKktJgWK/+cexdc6NP81ckxOC+gJPQ1CxNEWKGYobQooIuSzVm6NAcapOKSD9bimRudEH7cUusnYGkcik9CfolvJ6lDhnqcLvjbvNJ33u83pw2YmHAgCefn9H0rPKM24AY/LxQMRNRA4UWzI21rKUmHOTlbJUTLiJNRDybYQUIxQ3hBQJ4Yhq+m06V5kb8TpVljKJjOHcuJ3Qa+mW0nMt8epGXOytg/zkLeJ2geJEmRvBwcOrAZg70eyQZ9wARvu66CQDBjLnxihLCRI7N1koS0kb1IVo67cZmEhIsUBxQ0iRYC1D5aoVXOQxRBnGjrTn3CQIFNuJm0T7peRuqbhSVsBr2i4OxIsbfWt4CtEQ59zERI68CiMUTi9Q7ffGOzdWceNLY12EU4RLE/B5JHFD54YULxQ3hBQJXRanJpXTkCmEuxDwJv5xoXdLuZ1zo2duzGFaawgYSNwtlTRz47fL3NiLm1QdQlZxY/d+uJ5zI7WCW8tQFX7r+oX03DEniHJgmc+YvMyyFClmKG4IKRLkMDFQGs6Nkbkxuye2zo3f3rmRMzdxZamAF1YN4rPcSXxfmpY8EG0NFJf5By5u5Fbwcl+qQHF6AtIJohxodm5YliLFC8UNIUVCV39+nBvhpviTOjexQXkDnHOTLHOTqBVcCCRFiX9cuSPnxhARyUoxfQkyNzJuMzdCDPq8ij6hWBBflhLDCLOXuQl4WZYipQHFDSFFglXM9ORoRL7IeCR1bvQ5N+ld3K1zbmy0je5sxAeKox+9HiUuX2PXLWX9XP6+kg3hE4HiMku3lEzac25sAsWJuqWysX5BfN9lfq+054vihhQvFDeEFAlx4iZHzk3QhXPjtGSyp6sfmqbp+RHrbqmkzk3QvhXcbv2C3Zwba1nK61H025KVYhLNuZEJuRAEEVWDiPj4bAPF5sxNsm4pVdXwv6u3YuPODsevLyPe94DXo5fbrC33hBQTiXs7CSEFRVxZKkeZm5AUNk2Em8zN6s/34sx738JPjjsYwugRgWJxsZfnrQj0IX4WAaGadkvFr19IVZYCokIlHIw4cm6SlaXclIzk+8oTisV54ubzJOmWWrN1H371l48BAKuv/AZG1lU4PgcgOzcsS5HSgM4NIUWCEDdVsXJFrtYv6HmMJOLGzW6pj5rbAQAbd3ZIzo1ZKIwbWhn3uITrF6Q5N3brF+KH+MV/H04u6MLJqEhSlnKTuZHfK79liJ81bxO9T+JuqT1d/frfl6x4X3eznCIcq2jmhmUpUvxQ3BBSJIiyVENtOYDcbQU3ylL2e6UAd86N2FXVG4pIc27MP4rGDa2Ke1wicaPJmRurc+OP75ay5nIAQ7glc24SzbmRcePcRCQHxmdZv2DdCB69T+JuKfl1V3+xF39+p8nxOQDJuWG3FCkRKG4IKRK6YmJmeE0ZgBy2guudNPEXXIG+1NFBoFgsmuwNqXHdUoLxNs6NuOhaW8HFGoqAzxM3jbgy4I0b7GfN3ESfO7Vb0Ru3fsGuLOXcMQlJ75U1UGwNEwNSt5TNe2x9T95ranN8jnBE1UPZZT4v59yQkoDihpAiQTg39THnJlfrF3Tnxpcp5yZ67v4MOTft0u4ru1bwVOsXAMO5cRIoLk8WKHZTlooY83kURbGUpeIzR/4kgxKt70kyB8qKLOg4oZiUChQ3hBQJuriRnJtkE3UzhasJxemWpSwuyPhhSTI3lotuR58hbhTLEaND/FI7N+J7c1aWEoFim8yNC0Eg75USZ5XPbcWYcxP/HluFiJuSUjChuGFZihQvFDeEFAmdFnETUbWc/Had6QnF7bFFk73BiC6GyixlqbFD7MSN/foFIZZqbZybCn/8bim7zE2ZP7Vb0esgUOzKuZH2SsnPCyQIFCcp/fVJE4YBd66LuK8vtlWdgWJSClDcEFIkCOdGZG6A3ORunDk3iV0FK0KM9CUoS9XXlNmWZfT9U0nLUubHROfcmG9L17npt5SlBpq5EQMPhSMjB5TtxI0vSVlKCJG6Cr/prE6wildjzg3FDSleKG4IKRKEuKkt9+tiIBcrGPodODdVsbk0Xf2hlM9n3y1lXMzH2+Rt5PtYw7OZyNzoIVoXc24GmrkRQkgvS8mBYn+8uEu2fkE4N7q4cVOWisQmL/vM5TZxOyHFCMUNIUWCPuemzIeqmLPRm4MpsuJimmxC8eDK6EV1f3dycaOqmi5GZJdDFjd2M26i97Ffv9CRQNwoSvRCbS1DDTxzk6GylC5uPKbnBRKVpRLnmoTgqy2P/rtwU1ISjxViTS9t0bkhRQzFDSFFQnesFbym3IfKMm/stuw7N04yN0OqAgCiAiyZQOgKhmEXy5E3Yo8flsq5McSNpmmWzI1x/0p/tA3c6twky9wEk3VLhVPPuXEVKLaUpfxejy687MtSRq7JGiQXTo3h3LgoS1mGNDJzQ0oBihtCigTZuamMlS1ykbkJOsjc1JYbwqItFhi2o73H3tlx5tzEl6V6QxHdAbI6N6LjKNVuKQAo8zoIFAdjmRtfZjI3wrmRHTFRmrLvljLuZ30d4bLUplGWEo8t85kdKXZLkWKG4oaQIkDTNH2uTVWZN6fOTSiceiu4x6NgUGXUvdmXTNz02osbn1fRszCJMzfx3VLi+bweBVUBr2n9ghBD1m0LXrv1C/7UZSl9/UIgvixVE8scueuWMreCR88Rfe5k3VLyY/WzWZ0bN4Fii3jlnBtSClDcEFIERGfaRP9eLWVucurcJBE3gLPcTUcicePx4GsTh+Hg4VWYWF9te59yKfQrdid1xAYC1lX4oShK7E/0/sIFcTPnxlWgWHJShGMSVjXHe53kIX6CikB8W7hxbsm5CZtfQ7hZaZWlpKWZ0Y+cUEyKH24FJ6QIEA6NokQvfJX68swciBt9t1RycTOkKoDPd3djfxLnpi2BuPEowAOLjoSm2WdiAHPpqj+soiLgNXVKCbyKgrCm6e+R063gQOLFl+GIsSqiQneEFAS8HgQjKmor/Ghu6wUQzdKUeRKvqtCfU41/X42ylM2EYsm5sa5gEM5NbXkaZSlpaSbAshQpDejcEFIEiN/MK2IhWdF6nYvN4E4yNwD0slQycWNXlhLrBxRFSShsAHPoWJSm5DCxQIgZvSzlpBU8wQwd/fUkJ0QWWUII1FUYYsRp7kZvBZdES0WSspSiGKU766wbO+fG6fTquDk3LEuREqAgxM3y5csxfvx4lJeXY86cOVizZk3C+95777049thjMXjwYAwePBjz5s1Len9CSgHxW7S48IiLX3cONoPrQ/yS7JYCgCFC3HS7FzdO8EndRKJzyc650ctSARdlqRTOTa/kkMlZG/G4mnLj9UMORYEeKJbKTd+dNQZTR9dizkFDbB8jzm7N9gixVxsTWZrmXGQJEVNmCUpT3JBiJu/iZsWKFViyZAmWLl2KdevWYfr06Zg/fz527dple/9Vq1bhzDPPxKuvvorVq1ejsbERJ554Ipqbm3N8ckJyh/UCZJSlctgKnmQrOAAMjrWD70uSuWmz6ZayExuJEK6JKMfpzk254ZwIp0YvH1me3t65ST7bRd4rJW8ZF4+rLvNJU5odihtLKzgAnDt3PJ65+FgMrS6zfUyiWTfi34fsYDktK+mZG+uEYpalSBGTd3Fz66234vzzz8eiRYswZcoU3H333aisrMR9991ne/8//elP+MlPfoIZM2Zg8uTJ+P3vfw9VVbFy5cocn5yQ3KE7N37h3IiJwLkTN8m2ggNGoDhpK7idc6M4FzdjBlcAAJY9twkRaSBgnaksFf0onBtFUUwCJ+lW8ATCxDrAT1AmtW6LTEwi98eKPsQvRblPRgihcCLnplwWNyrCETVhh5pxP/sJxQwUk2Imr+ImGAxi7dq1mDdvnn6bx+PBvHnzsHr1akfP0dPTg1AohCFD7G3c/v5+dHR0mP4QUmz0W367rok5FdluBVdVTXcJUmVudOcmibix65byep2Lmxu+cwQCPg9e3tiKm17YZJpOLLA6N4BZ0CTL3CRqBReZFuv2cvHfo8Lv1V0Vp+Ugu1bwVIiOqbg5N2Fj75U8bfmC/1uHo65/GTvbexM+Z3zmhmUpUvzkVdzs2bMHkUgEDQ0NptsbGhrQ0tLi6Dl++ctfYtSoUSaBJLNs2TLU1dXpfxobGwd8bkJyjXW/U3VZbpwb2YVI3QqeXubGzcX9K2MH47enTwcA/OGNLdjd2Q/AIm5sdjXJpaRkmZtEF3SR8bEO1wtIGaiALm6ciQLrbiknJNoMLpfN5EDwxzva0R9W8f729oTPmThQzLIUKV7yXpYaCDfeeCMefvhhPPnkkygvL7e9z5VXXon29nb9z/bt23N8SkIGjnWKbHXMuensy664kS/UqVvBY3NuEkwhBoC23qjwGRpzeQDngWLBgmkjMbjSj7CqYe22/QCSl6UAc+nLdoifCBQnuKCLQLF1n5T4vFxybpyuYAg72NllxVieaTg3mqaZMllyZkY4e9v39SR8TqsraDzeeccVIYVGXsXNsGHD4PV60draarq9tbUVI0aMSPrY3/72t7jxxhvx4osvYtq0aQnvV1ZWhtraWtMfQooNay5CdOdkWty89cVe/MdfP9aDyvKFOmVZyoVzU19r/DLisxEbyVAUBYeOqAEAtHT0AUhQlgokKkvFP2dK5yZk79zIAW+RSXIeKI5vBU+F6KySMzf9pjZ1j9TWruqDB5sciBtrWcpNxxUhhUZexU0gEMCsWbNMYWARDp47d27Cx91000349a9/jeeffx6zZ8/OxVEJySvW366zVZa69aVP8cCbW/H6p7sBGGUpnyf5DBrAEDedSZZnit1SI+sMceNS2wAAJo8w/5JibgWPL0uZA8XxL6i3gicsS9lnbsTrDq4MpJG5MW8Fd4IeKJa6peTgb5nPq/8b6e4P62dxIm6su6WiX2NpihQneS9LLVmyBPfeey8efPBBbNy4ERdeeCG6u7uxaNEiAMC5556LK6+8Ur//b37zG/zqV7/Cfffdh/Hjx6OlpQUtLS3o6urK17dASNaxXoBq9bJU8k4Yt+zpimZYxFoDJ3ulBPJWblF+komoGjpjYqxhAM4NAEyOOTfyawv0spQsbjzJMzepBtf1Bc2rFwSXnnAIfj7/UJx4eIPrzE1YEo5OMQLFsnMTPZtHiWZyxH8ruTyYrCyVKHMTfW6Giklxkvf1C2eccQZ2796Na665Bi0tLZgxYwaef/55PWTc1NQEj/TD76677kIwGMT3vvc90/MsXboU//Ef/5HLoxOSM8TkXJGHEJmbrv4wNE0zBWYHgigp6WWpSPR1neRCvLHlmfu6g9jfHUJ9jTkH19kX0vdjNdSWmR7nlskjEzs3tmUpJVW3VCrnxr4sNbG+BhPro0LLbebGbkJxKvRAseQO6Z1csenVoj1dnhT95f5eqKpm675ZJ1ArSlQgBcMqxQ0pWvIubgBg8eLFWLx4se3XVq1aZfp869at2T8QIQWGuABZy1KhSDRMap2/kg6qNDemJyamgi6cGwAYVOmPihubdnDx3BV+r2keixvnQjCpoRqKAl0syc6N16ZbSr6o283VSdUKLgLF1rKUjOs5Nza7pVLh04f4Ga/Rl2BOjTxvKBhR0drZh5F1FXHPaRXO4jmCYTXhOgpCCp28l6UIIakRuQohMqoCPn3NQKZCxR19IYgoh7iYO90rJUi2gkEeuFduysO4FzeVAR/GDakEEF23UFNm/J5mXb8AWJwbG6ckkKL9WXdHbHY+Cfxuy1Kq+1ZwY/1CfOZGvKdlNmUpAGjaa1+asvtvzFk3pNihuCGkCLBmbjweBdUZnlIsXwzFegNjr5SzHxXJBvmJc1aX+1ARMJ7PTVlGRnRM1Zb7Tc6MdRYQYA4UJ8vcpCpLJXNuxOs6z9y4n1Dsd+TcxJelgMShYn3MgM1CUIobUqxQ3BBSBFhbwQF51k1mQsXyxVCIG2OvlENxo69giD9TT2zJZ1XAm3B6sBtEx5SctwGAn594KBYdMx7TxgzSbzOVpdIY4tebIFAsozs3YYfdUmkEiv02c27inJvYGa3uWaJQsa1zI2bdsCxFipSCyNwQQpJjdW6A6AqGne1AV4bKUvLFsNcy5ybVXimBsTwz3rkROZ7KgM/kEqSTuQGAw0dFxc2w6oDp9pOPGImTjxhpus3rMHMTVjVEVC1OAAlxWZEk2+Q2cxNKY86NnrkxBYpFZiZFWSqRuLGMGYj+nWUpUtxQ3BBSBBilA8m5iZVdOrNYlko3c7M31lIu0xM7Z1WZ1zKDJj1x843J9fjZNw/BsYcMS3lfWdDYiQm57BYMq3FdUYZzk8HMTRoTiu3WLyQqS1kXmCYSN702S0FZliLFDsUNIUWAfVkqs1OK5Ythb8hclnJ6AR4RG87X2hEvbrqDhnMji5t0Mzc+rweXnjDJ0X1l/WQ7xM+bXNwY7daJ3wfXc27SChTHDwpMFSgePagCzW292L7ffnmmEG7y98z9UqTYYeaGkCIgUVkKALqykLnpTTNQPCI2nE+sRZCRnZtyU+Ym+z+GvCmG+Pm9ii6A7C7oeqDYkXPjNHPjPlCsTyiOJHNuoh87Yv8uRPB6d2e/XsKS0VdLyM6N31jhQEgxQnFDSBFgl4sQ7c+Zcm5sy1I2r5sMMUdlZ3tv3NLFhM5NmpkbN8ilL7symKIouvNiV4pxVJaK5ZIcL85U0wgUi91SNusXrM6NMTCxXO8W67DZyt5rJ25YliJFDsUNIUWAcBNkB6WmPMOt4N3xZamQy1xIfWzycF9I1efaCMTU46qAF+VSK3i6mRs3mJybBGUwvR3cpqwkdkslDxS7K0ulM6HY2AqexLmxnLG6zKsPOeywuHyaphmZG+m/CctSpNihuCGkCLAuzgSA6rJY5iZjgWK5FTxsel2nZalyvxdDYx1TO9vNpanuWCt4ZZnP5IDk2rlJ1Hoe8CUuxfQ5cG7S3S3ld1GW89t2S5nzQFaXrSLg0ydCt/ea/630h1Xd4TE7N+yWIsUNxQ0hRYDdoDVjzk2mAsWJu6XcdPSIUHGLRdwIwVRpnXOTZqDYDakWZwKpnBsXc27cbgV349yICcVSt1S/ZcCgVdxUBryorYj+W7E6N3IGx5y5EXNuKG5IcUJxQ0gRYNctlZNAscvdUgAwMiZu4pwbKXPj93r0C3UunBtZPyQqg+mlmCShWyeBYse7pdIKFMc7N9YxAWWWKcqVAWOXlzVzI0pSAa/HdA6WpUixQ3FDSBFgV5bKZKBY0zRToDisagiGVX0ruNM5N4Ds3Jhbj/VuqVjLsRAKdkP1Mo2TzE0giXPjJlAcchko9qcxoVjulopzbizukryoNE7cJJi8zLIUKXYobggpAuxawaszGCjuCUbiunx6gxG9xOLOuREdUwmcm5go08VNDpwbxUHmxnBuEgeKMznnRry3br5/MecmaJO5KUuQuakM+KSylPnfSo/NjBv5OejckGKF4oaQIqBfH7Evl6UyN8RPlKQCPqNc1BMKu94tBSSeddMbNDs3YnlmukP83GCaUJwgwJvIuVFjLhbgtFvKaebGfZ7Jds6NpWQWV5Yq8+r7t6zOjd2MG8D4d+a0rZ2QQoPihpAiwG4Ngr5+IQOZGxEmHlzp13+L7wlG0goUO8ncAEYZJRfOTardUoAhCqwX9D7JvchO5iaNspQ850a4SokCxXJZqs8+c2P9vliWIsUOxQ0hBY6maUZZyuTcGGUp68A8twjnZnBlAJUxcdMrlarclKUSdktJE4oBoxSSi8yNaSt4isyNtRTTJ5Wpkoub+Bk0yTDWL7hwbjzxpa8+i6tnzdxEy1LCuTG7fHarF4DkJTpCigGKG0IKnFBE02eR2K1fUDUjO5EuYov3oEq/7qzIORy/C3dBiJuu/rDuKqmqZtoKDsiZmxysX5COn7IV3Jo9kjqKkrlMQhy5X5yZTqA43rkp89mXpSqStILbTSeOPkf0e+lj5oYUKRQ3hBQ4spMglxwq/F79YjvQULEoSw2pCugXup5gWL9QO12/AETFi8h4CPemLxzRBZpwboS4yUXmxtkQP/uVA1ZnJBF65ibszEVLK1AsWsHVJM6N3ZybBN1SiTI3usDtp7ghxQnFDSEFjnyxlS9ciqJIuZuBiRtRlhqUoCzlJnMDxOduxHRiRTGyIRWxi3EuMjceR5mb5OImWZgYSCNzk06gWAzxs3FuyhM4N9EhfmJCcYJWcEtZSgjQTK32ICTXUNwQUuDIKxAUy4U5U6HiVIFiN5kbwChN7YzNutGnE/u9utCoyOWcm9hrKIpZ6Mikcm6S5W0A95mbSBoTiv1JnJvyZJkbPVBsydzEMjWVfqu4EaVJihtSnFDcEFLg6G3gNgIjU8sz93T1AwCGVpXpoqM3lF6gGADGDI7Outm6tweAea+UoCJW+shJK7iDacgJu6VCqdvAgfTn3LgKFHvjnRt9zk3s/HJHnc+jIODzmFrB5fC5nrmxOjcB8e+KZSlSnPhS34UQkk/sphMLajK0X2pvV7QsNbTaXJZyuxVccGhDDQDg05ZOAOa9UoLvzRqNL/f34NQjRg7o7E4Q5pCTQLDdMEMg+QA/APDrgWKHc27SCBQLIWQ7odjGuRGiRQSKw2p0C7jI1CQquQlHsJtlKVKkUNwQUuAEbaYTCzJVltrbHXVuhlWX6Y7KQMpSk2Li5pPWqLixzrgBgFnjhuB/fzhnQOd2ihA1yUpgwvGIawUXe72cZm4czoYJqe4DxbZzbixLVWXnRojJCr8XPo+CsKqhozes/3dItFZCZG56QxFEVC0nuShCMgnLUoQUOMmcm4bYNOAv9/fGfS0R977+BW5+YZPpNtEKPqTKcG7SnVAMGOLmy/296O4Px00nzjVC1CS7SCeayuu0LOV6zk0arpjPMgVZVTVdgJb7PPp9RPlNiBhFUYxZN5IQTliWksqHzN2QYoTihpACRzgJdu7JoSOiImJTrPyTir1d/bj+uY1Y/urn2N0ZdWsiqqaLm/iylPvdUgAwuCqA4TVlAIDPdnXZZm5yiQgRJ9vAbTg39nNuUpWl3GRuVFWDMF/cbEUXSzaFMDJ10kniSwhhuQxYGythyu3g+hA/mzk3Qgh2M3dDihCKG0LyyPJXN+M/n9mQdMKwtewgI8TNJw7FzTtb9+t/b++NCpq2nqB+oR1SGTB3S6Xp3ADm3E1Pnp0bj4PMjXh/rc5Nv+NuKeeZG7mslExwWTHm3EQfLzIzgOHcAMb3YhI3Nu3giYb4KYqi/7diOzgpRihuCMkT72zdh5tf+AS///sWbNnTnfB+ycpSk0fUAgCa9vU4Cn++s3Wf/nfRFry3W6xe8MPn9ehtwb1p7pYSyLkbu8xNLnHSLSXEwRd7ukwCJ5G7YUUEip3MuZFbuV0Fii2lLzH12e9VTCJJ/FupkN5vu/1SiebcAEaei2UpUoxQ3BCSBzRNw83Pf6J/vj1JZkaUpezEzRCp/CPCu8kwiZvYb/CiU2pIVQCANJ02GE67FRwADh1RDQD4tLUzbq9UrhETij1JAsX/NGk4asp9+LS1C8v+tlG/vS/s1LkxhEeqXV+yu+MqUKx3S0Ufv7/b2Akmo5elpDPX2eyXSuTcAEYJkc4NKUYobgjJA69/tgdrJKHx5f6ehPdN5twAwGSHpanu/jA+3tGhf96pOzexGTfVUZGUqbLUIaIsVUjOTRKXpKG2HLf+6wwAwP3/2IpVn+wCYASKU4kbIRA0LfU2bbmV25/GnBvh/IjJ0kKYCkRnnbksFZ+5STZ9uUpvB2fmhhQfFDeE5IHfv/EFAEM0JOt2StYKDhji5oMv2/Ddu97EOX9429Y5eK+pTZ+KC0jiJubcDKsWzo3RBhxKsxUcAA6pjzo3rR392NEW/f7yl7lJ3S0FACdMacD3Zo0BALz+6R4AzgPF1WU+/fmtaw6siP8OniQTk+3wW4b47Uvk3PhFWUoOFDvvlgKA6pjLxlk3pBihuCEkD4iMzfGThwNILm6SlaUA4NBY7uaRd7/E2m378cZne7A7NnFYRnaKAGM2zl5pOjFgXOi6+sN6cNVNLkRQU+7H6EHRScXrmqJB5rx1SympMzcCIcraYq6I0/ULiqLoHUmpxE1IX73g7kewdYif3MIvI/6tyC3dtXZlqSR5IuGydTNzQ4oQihtCcoymadjVERUUR44fAiBFWUrvlkpelpJdmZ1tfXH3e2dLVNyIoKj4DX6P1AYOGBc1uXyRjnMDANMb6wBE3Rsgj3NuYsdPlrkRDI4JhX0xcdPvcM4NYORaUokbfcaNy+F4eqBYNWduEpWl5DPrreCOnRtOKSbFC8UNITlmX3dQ76iZOXYwgFTOTfKy1MT6alivkaIMJOgLRbA25p58fVLULRJlqX1i9UKVuSzVngFxM2vcENPn+Z9z40DcxEo8Qjg4LUsBUrt1TwrnJpKec6MvzhTOTUyADU7g3Ni1gsviJnnmRjh4zNyQ4oPihpAc09IRdVWGVQcwYVgVAGB3Z79pZolMqrJUud+LGY2D4FGAxiHRMtCOdrNz8+7W/QiGVYysK9fdlISB4tiFTu7ocRN6lZk1brDp8/xnblJ/H0OqoiJgn8uyFODCuYkFgt0M8JPvr2rRQYB6WarSb7qfcOFEJx1gZG7E2UIRVf9vnCxQ3EPnhhQh3C1FSI5pjYmbhtpyDKr0oyrgRXcwgua2Xhw8vDru/qm6pQDgnnNnY09XPx5790v8/u9bsNPi3Px9czQce/TBw/SLnJG5sXduBH6v4ir0KnP4qFqU+z3GCoM8r19wIiYM5yb6/vS6EDd2g/LsCOvOjduylPFvIKSqRqDY4twsOeFQzBo3GKdIS0kHxQSQcOpMAwAD8f+2qpi5IUUMnRtCcszOmKsyorYciqJgzOBKAIlLU8kmFAuGVZdh8ohajIwFeHdanJs3P4+Km2MmDkVNuTlYuqfL7NxY27WFG5EOfq8H08YM0j+vylMruMfB4kyByK909Ufn/DhtBQecOzf9ac4PkoPd4YimCzARBheMqCvHGUeONZ1ZhLtbO/sRjqi6aPMo9q3+VfqcG5alSPFBcUNIjmkV4qYuuvRyzODoRSdRqFjfLeUgnzEq9pw72g2h1NYTxIfN7QCAYyYOQ40ULA2GVX1SsXBurNmSS+ZNcvBdJWa2VJrK1xA/J4szBbXlfj3D1NYT1NcvZDJQrAeBLS3cqfBJZbVwRJMyN6kF6LDqMvi9CiKqhtbOflOnlGIj+tgKTooZihtCcozI3IyotYobe+dGhI8TdUvJ6M5NrFuqPxzB8x+1QNOiweOG2nJd3HT2hfUhcF6Pol+Y5QvdyVNH4OyvjnP3DVqQczf5GuInxIqTMpDHo+ilqX09QVeBYmMKcHJxI3JO1i6nVMjOTTCiJuyWssPjUTCyLpbJautN2ikFyEP8KG5I8cHMDSE5piXWFt2gOzcOy1IOShjCudnV2Ye3vtiLRfe/o1/EvjZxGABz14woSQ2pCphyNTd+5wh8sacbl594qLtvzoZZ4wbr4kIIq1yjl6UcZocGVwWwtzuIfd3BrASK9+rt92VJ72dFURR4PVH3ZX9PUJ9DZB3il4hRg8rRtK8Hzft7MW5o9N9dou+LmRtSzFDcEJJjWtsTOTeJylLJW8FlROkhFNFw92ufozcUQYXfi0kN1ThrzlgAhsDo6g9jd6cY4Ge+OP6/o8a6/bYSMqgygNv+30z0hSJ63ifXeF1kbgCjXLSro18XKoMcZI8ci5su82whN/hi4maXNDvIifACgNGDKgHsQ3NbL+prRcYqlXPDzA0pPihuCMkxelnKqXOTohVcxuNRMKKuHNv39eL1T3cDAP77zJmYN6VBv4/oltI0YPu+qKBK5yLrhm9NH5XV50+Fm8wNYGRYPviyHaoWzaXIbdWJcC5uosJkWJU75waIhrT7w6redWftlErG6EGxTFZbLw4bGR3+mChLZMy5oXNDig9mbgjJIX2hiH7ha4g5N0Lk7OnqNy1UFOjOjYPMBwA9V6Fq0Yv5nAnmQXplPo+e3diws9P0mFLFzRA/wMiwvP9lGwBg7JBK29CtFbdlKbeZG8DI3ezqdJ/bGTVIytwEk3eBibIU59yQYoTihpAc0hIrSVX4vfo4/KFVAXg9CjQN2BMrV8gYmRtnpQeRuwGAaWPq4kpB0R1I0ds27IxuCR87pNLld1JciJLSIIfZFJFh+SjWZTZ2qLP3JydlqVjXnO7cuOi4EuKm2U2gOBiBqsYvYiWkkGFZipAcIubPjKwr150Aj0dBfU0Zdrb3obWjT3dyBHoruMOZKKJjCgCOOXiY7X1qyn3Y2x3Eppi4EZONS5VTp42EBuCfJtm/H1aEGyJcs3EOxZ8Ia/eHVfSFIgldEdEtNcxloBgw9lElykslY3Qs39W8vxe9saBworJUtbQqoycUMX1OSKFD54aQHCJPJ5apj+U5xNdlnEwolpGdm6MnDrW9j3BzxHM3Di5t56bc78X3Zo1BfU156jsj3g1x6tzUlPkgqleJ2sE1TcuIc7Or033mZlSs/NgdjOhlrUTiptzv0bvc2A5Oig2KG0JyiDVMLKiPiZ3W2AVHRlxYnHbEiPxMud+Dr4wdbHsfa0t2qZel3GLNsTh9fzweJW6Hk5WO3rDewp1O5sY3gMxNRcCr33/zri4AQHmCspSiKEY7eAJxw3IVKVQobgjJISJzY3VuGmJtubstzs3uzn509IWhKM7dlTkThmB64yCcf+yEhIKoVsrhlPk8jjqBDiSsbsi4oVWOH5sqdyNKUjVlPsc5KhmxxFS0grvJ3ADRWTcA8PnuqLhJNnk5WTv4xp0dmPnrl7DsbxtdvT4huYDihpAc8mlrtDtJbAMXNMTKJa0dZudmU0s0E3PQ0CrHSydryv34y0XH4LIkA/hk52bM4ApHnUAHEvJaBI9i7GVyQmpxk35JCjCcGxEIduv+iO9FODfJxU3idvDfvfQp2ntDeOb9na5en5BcQHFDSI7QNA0bYwHeKaNqTV8TA9VaO83OzaZYq/bk2EySTCF3ULEkFY+8q2nUoApXCy7rpAnQduy1LCp1i8+yY8ytuBEdU6KiVJUkKCxCxD2WKcWftHTixQ2tAKKdV6m6wwjJNRQ3hOSIlo4+7O8JwedRMLG+2vQ1PXNjcW42xpybySPMYmig1FYYF7RGips4qst8+jyZcQ7DxALduelJ4dykkbcBjG4pgZhw7RRZzE4bU4cF00cmvK/YBWZ1bu5atdn0uXAkCSkU2NtHiERE1XDj3zZiWHUZfvRPEzJartmwIypUJtZXx2VhRFlqlyVz80lLzLkZkT3nptQ7pdJBUaLLM3d19mPsEOd5G8BoB2/vjQqC5rZefNLSga9PqofXowyoUwoAPouVkwBg8fETdSfGKd+ZOQYtHX04ctwQfPOw+qT/xu0yNzvaevHX93cAiAqlpn092LSzA0eOH2L7HG7oC0VQ5vOwTEoGDJ0bQiSe+WAH7n1jC5b9bRP+5/UvMvrcQtwcNjLehRGB4r2xRY0fNbcjFFHxWWv0QpZp50bO3NC5sUeUe9J1bj7b1Ynv378GX/vNK/jBA+/ilhc/ASCVpdJYvQAAcw6KiogJw6uw5IRJrh9fV+nHlScfhnlTGlKKiOpY5kbulnryvWaoGnDUQUNw6rSo67OpZeDOzZub92Dq0hfwn88yoEwGDp0bQmKoqob/fsWw22/82yaMHlSBBdJepOa2Xjzwjy147sMWTBhehQcXHaWP9v/gyzb86a0m/PykQ22Hs4lpwFNsxM3gygB8HgVhVcONf9uEB97cim9MrkcwoqIq4HVdekhFrUnclPYAv3QZO6QSm1o6bf97JUOIm2c+MAdt73n9CyyYPgp7Bhgo/vn8Q3Hk+CE466tjTZvcs0FlzLkR+SFN0/DY2i8BAKfPGqNnkQYqbjRNw80vfoKwquEPf9+CU6eNTDjGgBAn0LkhJMbfPmrB5l1dqC336Ru0/+OvH+thyr5QBP9692rc+8YWNLf14o3P9uDvm/cAAMIRFT97eD1WvLsdd7yy2fb5NyQIEwPGlGIAePidJgDAK5t2AQAOHVGT8YuY3ApO58ae//z2VPxh4Wwce4izqcYCOc/k8yh4evHXcNLhIxBWNVz5xIfGZOE0A8WHNNTg/H+aoOdhssmkWDbsyfeaEQyrWNfUhi17ulHh9+LkI0bqjuKnLZ3QtGhC+YWPW/DbFz5BW495lUiymThvfbEP7zW16Z9f/eRHtnvWCHEKxQ05YOgLRfDvT36IG/+2CXu6zMHdHW29uOmFTQCAH3ztIPzHtw7HuKGV2NsdxB9XbwMA/N9b29Dc1osRteX4xuR6AMBDb0eFyJPvNWPLnm4AwF/WRy8EMl39YWzbG93AbVeWAoxQcV/I/NjJLp0DJ4gdS0OqAiahQwzqa8rxzcNSl26sCOcGAL5/9HgcMaYO//Gtw1Fd5sP67W1Ys2UfgPQDxbnkjCPHor6mDF/u78VDb2/Dw2ui/95PPmIEqst8mDC8Cn6vgs7+MJrberFtbzcufug93PHqZpx02xu4a9XnuPmFTfj2nf/Awf/+HP5r5We2r3NnLKD8z9NGoq7Cjw07O/Dnd7bn7PskpUdBiJvly5dj/PjxKC8vx5w5c7BmzZqk93/00UcxefJklJeX44gjjsBzzz2Xo5OSYuY/n92AP73dhLtf+xzH/uZV/Qf1R83tOG35P7Btbw9G1JZj0dEHwe/14KffOAQA8D+vfY5dHX24a9XnAIBLTzgEvzxpMgDgpY2t+HJ/D26Xfmjv7wnhlU27sKujTx/aJ3Y4jawrT9i6K3I3ADCpoRpVsbk2h2U4TAwAh42swQVfPxjXfuvwjD/3gY4I+A6rDuCn86L/hkbUleP6b0813S/dslQuqQh48dNvRr+H657ZgEdjJanvfWUMAMDv9eDg4VF3Z9POTvzm+U0IRlR4lGh34G+e34Tlr36O95raoGnA717+FG98ttv0Gu9u3Yc3PtsDr0fBL0+arOeI7np1c9wvCYQ4Je/iZsWKFViyZAmWLl2KdevWYfr06Zg/fz527dple/8333wTZ555Jn74wx/ivffew2mnnYbTTjsNH330UY5PTvJJMKxi484OvLN1X5z9bcfzH+3E/70VFTOTR9SgNxTBFU98iF8/swFn/M9q7Orsx6SGajz+k6NRVxn9zftfZozCQcOqsL8nhK8uW4m93UGMH1qJ73xlDA4dUYPZ4wYjomo4+bY38OX+XgyvKcP3jx4PALjphU049qZXcdxvX8U/Nu/BHa9GfzM9fFRdwjPKU4u/+5Ux+N0ZM3DS4SPwremj032bEqIoCq44ebIpT0Qyw8zGQbj1X6fj4R991eSK/cuM0fjh1w7SP083UJxrzjiyEeOGVkLVotOsL/7GRMw92NhZJjr5bnnpUzz3YQs8CvD4hUfjZ988BAumj8LZXx2L6789Fd+bNQaaBly64n00t/UCiJZzf/WXjwFEMzyNQypxxpGNGF5Thh3tfXhqfbPjc2qahre+2Isn1n2JNzfv4T6sAxxFE4XSPDFnzhwceeSRuOOOOwAAqqqisbERF198Ma644oq4+59xxhno7u7GM888o9/21a9+FTNmzMDdd9+d8vU6OjpQV1eH9vZ21NZmzu7vD0f0WrqM03fX7n4a7B9sf1+7+8Xfmug49ue0ebzD1058X/vn7AtF0BdS4fcqCPg8CPg88Hs98CgKRFGgua0XG3Z04NVPduHdrfsRlGryYwZXYM5BQ+HzKNjV2YdxQ6swvbEOY4dU4h+b9+KO2G+BP/76BFxx0mT8+pmNuO8fW/THHzNxKO46e1ZcieblDa24+M/v6dNg7zzrKzjliGiHyJPvfYlLV7wPAKgMeHHrv87AxPoqzLv1ddv3o9zvwSM/notpYwbZfv2OVz7Db1/8FACw8rKv678Rk9IhHFFxzV8/hldR8OvTpqZ+QIHwxe4uvP7pbpwybWTc8tFH3tmOXzz+gf75/zuyETd+d1rcc/SFIviXO/6BT1o7MajSj6tOPgxb93bjzlWfo67Cj1cvP053Ne95/XPc8Fw00D9nwhB09oVRX1OGhtpy1NeUoSLghd8b/Rnh8yrwKgr+761t+mBBIFr2u/gbE3HEmEHoDUbwYXM7PviyDR82t6O6zIdvzxyNcUOr0NYTxP6eEHpDEYwbUomxQyvh8yjwehQoULCppQPrmvajoy8qlmaMGYRZ4wcjIH4+KYj+gQJP7O9A7PbYWRTp55iochq3xJNOJ3yhdc8HfB7Hi2qd4ub6nVdxEwwGUVlZicceewynnXaafvvChQvR1taGv/zlL3GPGTt2LJYsWYJLLrlEv23p0qV46qmn8P7778fdv7+/H/39hujo6OhAY2NjxsXNuqb9+M6db2bs+Uhqasp8qK3w678FpmLeYfW46+xZ8Hs90DQNVz35Ef68pgnfmzUGN3z7iIRTaCOqhh1tvQhGVJPgUFUN//vWNlSX+XDS1BH6TJBF96/Bqk9345JvTsK726KWu6IAd531FZw0NfHAtMfWfonLH30fE4ZV4ZXLj3P+RhCSZ97f3oZ/fL4Huzr6cem8Sbr7aWX7vh785E/r8GFzu+n2/zxtKs7+6jj9867+MI658RXXk499HgWzxg3Gtr09+pJakh++MnYQnvjJMRl9TjfiJq+t4Hv27EEkEkFDQ4Pp9oaGBmzatMn2MS0tLbb3b2lpsb3/smXLcO2112bmwEnwKArK/fYXRzuFbqey7YR3ojCj7a0Dfc4BPD7RLw32LxV/Y7nfgzKfBxFVQzCsoj+sIhRRoSHq7GiahobackwYXo2jDx6K4w4djoOGVUFRFHT0hfBeUxvWbt0Hr8eDodUBbN7VhQ07O7CjrRcBrweXnDAJC6aN1M+uKAqWfecIXH7ipJRdK16PYttR5PEoWBgrQ8ncdfYs9IdV1FX40RuM4I5XP8Pho+qSChsAOHnqCKz+fC9Om8lSESkupjcOwvTGQSnv1zikEo9feDTueOUzvLxxFzQAU0fV4syjxpruV13mw29Pn46/rG/GxPpqDK0uw+6OPrR29GN3Vz/6QhGEIipCEQ1hVUUorGHkoHL88qTJOGxkLUIRFQ+/sx1/frsJXf1h+DwKJo+swbQxgzBtdB227evBX9fvQH84gsGVAdRV+lHm8+CL3d3Y2d4HVdOgqhoimhZzj4aioaYMvSEVb32xF5+2dkbvowmHWvw96k2rquRRa4a7LbyEZI5CIrshkZOf7DH5xO/Nb+olr87Njh07MHr0aLz55puYO3eufvsvfvELvPbaa3j77bfjHhMIBPDggw/izDPP1G+78847ce2116K1tTXu/rlybgghhBCSPYrGuRk2bBi8Xm+cKGltbcWIESNsHzNixAhX9y8rK0NZWXEE9wghhBAycPLqGwUCAcyaNQsrV67Ub1NVFStXrjQ5OTJz58413R8AXnrppYT3J4QQQsiBRd7XLyxZsgQLFy7E7NmzcdRRR+G2225Dd3c3Fi1aBAA499xzMXr0aCxbtgwA8LOf/Qxf//rXccstt+DUU0/Fww8/jHfffRf33HNPPr8NQgghhBQIeRc3Z5xxBnbv3o1rrrkGLS0tmDFjBp5//nk9NNzU1ASPxzCYjj76aDz00EO4+uqrcdVVV+GQQw7BU089halTi6etkhBCCCHZI+9zbnJNtubcEEIIISR7uLl+531CMSGEEEJIJqG4IYQQQkhJQXFDCCGEkJKC4oYQQgghJQXFDSGEEEJKCoobQgghhJQUFDeEEEIIKSkobgghhBBSUlDcEEIIIaSkyPv6hVwjBjJ3dHTk+SSEEEIIcYq4bjtZrHDAiZvOzk4AQGNjY55PQgghhBC3dHZ2oq6uLul9DrjdUqqqYseOHaipqYGiKPk+Tt7p6OhAY2Mjtm/fzl1bWYTvc27g+5wb+D7nDr7XBpqmobOzE6NGjTIt1LbjgHNuPB4PxowZk+9jFBy1tbUH/P84uYDvc27g+5wb+D7nDr7XUVI5NgIGigkhhBBSUlDcEEIIIaSkoLg5wCkrK8PSpUtRVlaW76OUNHyfcwPf59zA9zl38L1OjwMuUEwIIYSQ0obODSGEEEJKCoobQgghhJQUFDeEEEIIKSkobgghhBBSUlDcEADA1q1b8cMf/hAHHXQQKioqcPDBB2Pp0qUIBoP5PlrJcf311+Poo49GZWUlBg0alO/jlBTLly/H+PHjUV5ejjlz5mDNmjX5PlLJ8frrr2PBggUYNWoUFEXBU089le8jlRzLli3DkUceiZqaGtTX1+O0007DJ598ku9jFRUUNwQAsGnTJqiqiv/5n//Bxx9/jN/97ne4++67cdVVV+X7aCVHMBjE6aefjgsvvDDfRykpVqxYgSVLlmDp0qVYt24dpk+fjvnz52PXrl35PlpJ0d3djenTp2P58uX5PkrJ8tprr+Giiy7CW2+9hZdeegmhUAgnnngiuru78320ooGt4CQhN998M+666y588cUX+T5KSfLAAw/gkksuQVtbW76PUhLMmTMHRx55JO644w4A0T1yjY2NuPjii3HFFVfk+XSliaIoePLJJ3Haaafl+yglze7du1FfX4/XXnsN//RP/5Tv4xQFdG5IQtrb2zFkyJB8H4OQlASDQaxduxbz5s3Tb/N4PJg3bx5Wr16dx5MRMnDa29sBgD+PXUBxQ2zZvHkz/vu//xs//vGP830UQlKyZ88eRCIRNDQ0mG5vaGhAS0tLnk5FyMBRVRWXXHIJjjnmGEydOjXfxykaKG5KnCuuuAKKoiT9s2nTJtNjmpubcdJJJ+H000/H+eefn6eTFxfpvM+EEJKKiy66CB999BEefvjhfB+lqPDl+wAku1x22WX4/ve/n/Q+EyZM0P++Y8cOHH/88Tj66KNxzz33ZPl0pYPb95lklmHDhsHr9aK1tdV0e2trK0aMGJGnUxEyMBYvXoxnnnkGr7/+OsaMGZPv4xQVFDclzvDhwzF8+HBH921ubsbxxx+PWbNm4f7774fHQ2PPKW7eZ5J5AoEAZs2ahZUrV+rhVlVVsXLlSixevDi/hyPEJZqm4eKLL8aTTz6JVatW4aCDDsr3kYoOihsCICpsjjvuOIwbNw6//e1vsXv3bv1r/M03szQ1NWHfvn1oampCJBLB+vXrAQATJ05EdXV1fg9XxCxZsgQLFy7E7NmzcdRRR+G2225Dd3c3Fi1alO+jlRRdXV3YvHmz/vmWLVuwfv16DBkyBGPHjs3jyUqHiy66CA899BD+8pe/oKamRs+N1dXVoaKiIs+nKw7YCk4ARNuSE10E+E8ks3z/+9/Hgw8+GHf7q6++iuOOOy73Byoh7rjjDtx8881oaWnBjBkz8F//9V+YM2dOvo9VUqxatQrHH3983O0LFy7EAw88kPsDlSCKotjefv/996csf5MoFDeEEEIIKSkYqiCEEEJISUFxQwghhJCSguKGEEIIISUFxQ0hhBBCSgqKG0IIIYSUFBQ3hBBCCCkpKG4IIYQQUlJQ3BBCCCEHIK+//joWLFiAUaNGQVEUPPXUUwX1ehdccAEURcFtt93m+rUobgghJcEnn3yCESNGoLOzM+V9N2zYgDFjxqC7uzsHJyOkMOnu7sb06dOxfPnygnu9J598Em+99RZGjRqV1mtR3BBCCpbjjjsOl1xyiaP7Xnnllbj44otRU1OT8r5TpkzBV7/6Vdx6660DPCEhxcvJJ5+M//zP/8S3v/1t26/39/fj8ssvx+jRo1FVVYU5c+Zg1apVWXs9QXNzMy6++GL86U9/gt/vT+u1KG4IIUVPU1MTnnnmGVd7dxYtWoS77roL4XA4ewcjpIhZvHgxVq9ejYcffhgffPABTj/9dJx00kn47LPPsvaaqqrinHPOwc9//nMcfvjhaT8PxQ0hpCD5/ve/j9deew233347FEWBoijYunWr7X0feeQRTJ8+HaNHj9Zv27ZtGxYsWIDBgwejqqoKhx9+OJ577jn96yeccAL27duH1157LdvfCiFFR1NTE+6//348+uijOPbYY3HwwQfj8ssvx9e+9jXcf//9WXvd3/zmN/D5fPjpT386oOfxZeg8hBCSUW6//XZ8+umnmDp1Kq677joAwPDhw23v+8Ybb2D27Nmm2y666CIEg0G8/vrrqKqqwoYNG1BdXa1/PRAIYMaMGXjjjTfwzW9+M3vfCCFFyIcffohIJIJJkyaZbu/v78fQoUMBAJs2bcJhhx2W9Hl++ctf4sYbb3T0mmvXrsXtt9+OdevWJdyM7hSKG0JIQVJXV4dAIIDKykqMGDEi6X23bdsWJ26amprw3e9+F0cccQQAYMKECXGPGzVqFLZt25a5QxNSInR1dcHr9WLt2rXwer2mr4lfEiZMmICNGzcmfR4hhJzwxhtvYNeuXRg7dqx+WyQSwWWXXYbbbrstoXNrB8UNIaTo6e3tRXl5uem2n/70p7jwwgvx4osvYt68efjud7+LadOmme5TUVGBnp6eXB6VkKJg5syZiEQi2LVrF4499ljb+wQCAUyePDljr3nOOedg3rx5ptvmz5+Pc845B4sWLXL1XBQ3hJCiZ9iwYdi/f7/ptvPOOw/z58/Hs88+ixdffBHLli3DLbfcgosvvli/z759+3DwwQfn+riEFARdXV3YvHmz/vmWLVuwfv16DBkyBJMmTcJZZ52Fc889F7fccgtmzpyJ3bt3Y+XKlZg2bRpOPfXUjL7e2LFjMXTo0Dinx+/3Y8SIETj00ENdvRYDxYSQgiUQCCASiaS838yZM7Fhw4a42xsbG3HBBRfgiSeewGWXXYZ7773X9PWPPvoIM2fOzNh5CSkm3n33XcycOVP/f2DJkiWYOXMmrrnmGgDA/fffj3PPPReXXXYZDj30UJx22ml45513TGWjTL5eJlE0TdMy/qyEEJIBfvSjH2H9+vV45JFHUF1djSFDhsDjif+d7Omnn8Z5552HHTt26PmASy65BCeffDImTZqE/fv34yc/+QnGjRuHFStWAAC2bt2KCRMmYMuWLRg3blxOvy9CSHahc0MIKVguv/xyeL1eTJkyBcOHD0dTU5Pt/U4++WT4fD68/PLL+m2RSAQXXXQRDjvsMJx00kmYNGkS7rzzTv3rf/7zn3HiiSdS2BBSgtC5IYSUBMuXL8df//pXvPDCCynvGwwGccghh+Chhx7CMccck4PTEUJyCQPFhJCS4Mc//jHa2trQ2dmZcgVDU1MTrrrqKgobQkoUOjeEEEIIKSmYuSGEEEJISUFxQwghhJCSguKGEEIIISUFxQ0hhBBCSgqKG0IIIYSUFBQ3hBBCCCkpKG4IIYQQUlJQ3BBCCCGkpKC4IYQQQkhJ8f8Bc8asEnWDP3AAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# Intensity in time domain\n", "srwlpy.SetRepresElecField(wfr._srwl_wf, 'time')\n", "[t, intensity] = get_intensity_on_axis(wfr).T\n", "data = get_intensity_on_axis(wfr)\n", "\n", "fig, ax = plt.subplots()\n", "ax.plot(data[:,0], data[:,1])\n", "ax.set_xlabel(\"t (s)\")\n", "ax.set_ylabel(\"Intensity (a.u.)\")\n", "\n", "fname = \"I_t.txt\"\n", "np.savetxt(fname, np.vstack([data[:,0], after_shift]).T, header = \"t(s) intensity\")" ] } ], "metadata": { "kernelspec": { "display_name": "simex-lite-dev3.9", "language": "python", "name": "simex-lite-dev3.9" }, "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.9.16" } }, "nbformat": 4, "nbformat_minor": 5 }