{ "cells": [ { "cell_type": "markdown", "id": "heading", "metadata": {}, "source": [ "## IR Spectrum of H2O Dimer from MD\n", "\n", "Run an MD simulation with GFN1-xTB for a gas-phase water dimer, then compute the IR spectrum from dipole-moment autocorrelation with AMS trajectory analysis.\n" ] }, { "cell_type": "markdown", "id": "runtime-note", "metadata": {}, "source": [ "This workflow can take some time to run, depending on hardware and MD length.\n" ] }, { "cell_type": "markdown", "id": "439ce65e-7d2d-46d6-977e-2bf7c979a158", "metadata": {}, "source": [ "### Initial Imports" ] }, { "cell_type": "markdown", "id": "imports-md", "metadata": {}, "source": [ "Import PLAMS classes used for MD and analysis.\n" ] }, { "cell_type": "code", "execution_count": 1, "id": "imports-code", "metadata": { "tags": [] }, "outputs": [], "source": [ "from pathlib import Path\n", "from typing import List\n", "\n", "from scm.plams import AMSAnalysisJob, AMSNVEJob, AMSNVTJob, Molecule, Atom, Settings, view" ] }, { "cell_type": "markdown", "id": "cf669b71-a1b3-4776-b677-4307ce068cbe", "metadata": {}, "source": [ "### Create Dimer" ] }, { "cell_type": "markdown", "id": "structure-md", "metadata": {}, "source": [ "Create an H2O dimer geometry and visualize it.\n" ] }, { "cell_type": "code", "execution_count": 2, "id": "structure-code", "metadata": { "tags": [] }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAASwAAAEsCAYAAAB5fY51AABPWklEQVR4Ae29CZhdVZnvvfc+p+ZKZR4qqYwkgEyCzEp3gtpiD9DXgYi0rVfsB1pE73XoT9uBJHajn32l0Vb4bO7nZ8t1YLg0j2hjCxdIGkQbZZ5DIGSsVOaqpIZUnbP39/uvc9bJrqISKpWhTtV51/O8tfa81/mfs3/1vu9ee+0gsGIKmAKmgClgCpgCpoApYAqYAqaAKWAKmAKmgClgCpgCpoApYAqYAqaAKWAKmAKmgClgCpgCpoApYAqYAqaAKWAKmAKmgClgCpgCpoApYAqYAqaAKWAKmAKmgClgCpgCpoApYAqYAqaAKWAKmAKmgClgCpgCpoApYAqYAqaAKWAKmAKmgClgCpgCpoApYAqYAqaAKWAKmAKmgClgCpgCpoApYAqYAqaAKWAKmAKmgClgCpgCpoApYAqYAqaAKWAKmAKmgClgCpgCpoApYAqYAqaAKWAKmAKmgClgCpgCpoApYAqYAqaAKWAKmAKmgClgCpgCpoApYAqYAqaAKWAKmAKmgClgCpgCpoApYAqYAqaAKWAKmAKmgClgCpgCpoApYAqYAqaAKWAKmAKmgClgCpgCpoApYAqYAqaAKWAKmAKmgClgCpgCpoApYAqYAqaAKWAKmAKmgClgCpgCpoApYAqYAqaAKWAKmAKmgClgCpgCpoApYAqYAqaAKWAKmAKmgClgCpgCx1qB8FifcCyfLwkCp+cdQRDpcz4XBMlyTNOscLWmrZgCpoApcMwVWAaYHgyCLCTKvNHJ2aa0LdP2j+KNBLP1psAgCtiFM4goB1sk2CwHUFh+oNd005w5EyFX7cQ4nt3EQTriuL0+m936uzDsXrFuXU/6uIIcnliwlOOkl9u0KWAKHFgBA9aBtem3RqBSqJcGzHdnzz55Qhy/oyFJzq4JglOrk2ROTRhW1YdhYxV79yVJvipJOuIwbEfotRzjN0w/vCWf//2fbNmyTScY7Lj9TmwzpoApUFLAgFWS4sATt+MNeVAta26eMjWK/nJCEFwKqM6ZHkWZBnbNJkmQpRaoSGAlmsbCbBi6eFExo8TuxvbG8bbeIPjl3jD84VmbNj3AIoFLea+EbZi0YgqYAoMpYMAaTJX9y0Jg5byqj7W0TFqUJF+cHAQfmRGGUxoBFB5VUB8EuVqgVM12eFdhNdNV1AIWkEpEoTBhQ4GIOp+wKAwjAaydxXuD4L7dcfw/3tLaeh+LRKsMX4qFiRLDiikwQAED1gBB/OwyALQiCGLN/01z8wfxpP7vljCcUx/HglS+ATBRR3XAScCqjaIAYAUOWOwTMd1PXDFrvyX5OI5zSRJmgFcHy3cmyU/X9PZ+9o+3b299EOfsQkDo22K1KWAKFBTod02ZKAUFLi0kxPPntbTULY7j78yKoo9NBCp4VblG1jUCKQELWAV1gEomWAlaAAiXKmVeVMEK2DloqS4a3MqT61LoGG2J49a2IPj0mZs334ZnpzBUwGRHK6aAKSAFDFgDfgceVhdPnz7/+EzmJ3Oj6LwGoDIercYBFYAVNAAmEutBfRFWAlZtJhOE1IE3Dy0df79nVQBVnohPwPI10/l8PsfSLCFisCkIrjtt06YvC1q0J+ZLMmhJRysVr4BSLVaKCiwrhIF5wWp2JnP/jCiaXxXHfYR+VS7sA0LypGqAkjeBysGKOpAJWKo9sFR7YHnPSqDKEfFpW01j7JENcbcIOeNFUfSlp5ubJ57W2voJQQtaGbTsV2oKoIB5WPt/BhFixOe2tMx6Uxw/NC8M548nBJwASJqAzji2awQw8q6cAaX6omWycF+Q8rWmBSMPrTSwioBywBK00sY6eBXsA5JJFFWtzudvOgNoPciRL7Sc1v5vyqYqVgGuKisoEBJ6hcctXFgzK45vnYZnxR3AXC2g4O6f66rAnT/d/QuqAJEzoFSNZaqq2ACrZktZDT2yvNVyhMHMr9f22legK8Iu4vh4b1UJ0FqYyVz9SHPzFwUrQcu+KVOg0hUwYPELWIx/RKfQ/PyurhXA6gKFgWAnC0ocpNSXKgtIMpibBlRVMg+awaDlQSU4+WnVaVh5YHlocUx5aoIWdxuz3EnMtUTRdauam98laCk8rPQfrH3+ylag4v9rX1qAVe7MlpZzmpLkc9VAAr9HKHIujQihO3/qpqBa0HImsAgwA6GVhg/burBQvzF/l7AIpVK4qHUqPmws1vSQDwGnkvzJxDD85/sWLDj9na++uod8Fq2wJHxBNPtbaQpUvId1ki7+xYuzDfn8d5pgE7BSp89QoBKknPlpACTvx0PL5a08tLyX5euBIaKf916Vrz3giiGhO2YRdCT4I0CWnxlF8+q6upYDKnVzqPjvrNIuUvu8+xWo6B8/oWB2BRA4afXqS8dlMufweE0eUGUkigxAFAxw0fWqYMBE3RcErhKwvNfkvS1fp2HkQZau/XYDYaVj6xyck9Awk1O3iii65pc8u0ib8hYa7v8B21RlKaDrsmLLSi5+Pry8qr+t5bEZhYESROZh1S90AyBuvggUQeV15uHl64FQ8vPeM/PbqfbHK55HkCQMxc0KkklRlK3O5b5M06yYAhWrgK7NSi16Zi+Z0dJyAT3UT80ALOa1rFCARSJw+KJ5pv0yV2u9tzRs/LRqrU/PD3W6eFx9QfL6+uI44SHr//KjGTPmLgW0ywpc9a2z2hSoCAUqFliEg45GtXH8IXUKRQh1znTFgYkp1UoaeROk3LJi7TbWH8HlQCW9zk+rHoq5Q7s8mrys/KRMpnZSGF6mUy0xYEkGKxWmQKUCK1yljpg8K4j38q5IPdABwEBQlYAFXLSF4kdvdPl08KJigi0PZDp2ep3f3u2Y+pMGmBYX59UGFSFOx6F+r2aXFJrCpBVToHIUqFRguc89Pp9/E4ntFr5uIcUBy3tTHlCuZq0DleqiMdJCEGP+IeZBa/VqF7B8remBANMxVHxdmCv91dqiZXhIOiEJf8b/O3Pm8YBLe1Tq91fSxyYqS4FK/cFzvdOHkwH4SGpLg7ygIDg5QGnBAJNH5Q1wBAIWd+8K2PAg8o/dqNYjN4PV6W3S+ws/3jiXSqlNeFvF6ZjHgqroK3au1q80YEkGKxWkQEV3HKVrwskeCgMBpXkPKFcXISVYyXqBDc/8BVkApH5ZDk4+ya4fkOAjIGmZL5qXeZh5ePnl2idlbOnAqVpGcQMCZjKZ0wqz9tcUqCwFKhVY7vonwjpOYR1wCl8HKH4H8qL68G48pFwNXAQreqE7E7A0vEwoCCnv5IvAI1gNBiyBqq+vP7gGQEsgdV5csVb7XDvV3iQ5RadZUuKY5qyYAmNfgUoFlnigBNC0QUHFOsEpB4BK0GJZL8t4ZCbIAhdn8q7YRp1IdafRFdY7L0lQUt8qv1wrtU7LBacBozSUlmsd27nzskvay/NtBbQT3bnYslhbZQpUhAKVCiz/5WrQvABfZ1DrZbl7ABqACFbuRRMAxT1TqBr4lHrAs63GeGf89gKQBCvBaSCwBCQtl3lo+XnVrFdiP+3VeWjRzlDtoJ56ZXNzfdja2sXZdAYDF/pbGfsKVDqwXG5IYPLQ0nR6XuDQvIMV03QwDTJAJQIuyk45YFHj9bi7hiTyA9dNQvBRODgYsA4ELZYrke/yYxxHcEq3h+lwX6F9k/P0yWKyC7NiClSMAhUNLDyTrJJZRQiU4OAh4cI/gCORHKiKMHJhIGARrOQleVi5fBjLa4CVhqNxQyYfCFhsV/K02D7Gcphg5SDJuVy7BoDLtS0Me7dnMmq6FVOgohSoVGBBGhdGbdW3LQgIDmlwaVriOBOomBa0VBc7mjJFYZniMQ8teUh9wKoac2NnASzniRW31fYK+xigz5lAxXjuzrNy+7Kex3ActPYx7eGp9uAFxuTV9DD0up+tW9du4aBEtVJJClQqsMQQuUivyEsi9Ep6WOChxcPQLnflYeU8KpZFAESkY+x1/haKYKU5eVfurh6g0p3FXmqNUOryXawno1XIb7GdgxzH0HDIDNJX8K6oBawStAQrzLfJ12pnZ6GpOrWaY8UUqBgFKhVY7gsGGM/pDp+A043VYYwPWki0U0scgcZ5VcDDe0qOEsBFxPCelQOWgFMElWDlPSx/DCXknWvHdv2gxXKBS9YPWBxf0BKkZLQx6eG4TL/EbLC80Dzl5K2YAhWhQKUCS6xRDupRgCNeZeTBCFppYAk0DjaCFaBwdwBZpqIDyEPSA9Glx3WYFnA05rvzzpgXtAQ6HccDy8NKsJM5L6tYOy+N6XQeqwirwmvuOc6eIHiSyoopUHEKVCqwBCm6SWVeIH+0CRjNxk2JueUWCVgONlqPea9K4aC746ea4oDFdCxgyZiWd6W7ja5TKcsUDupYOoY3Bz0dA/PAcsl65kvgY1rH8fkrgVRGKJjZDmD3hOF/MBs8X2iGJq2YAhWhQKUCS7zJbty4sXvmzJm/4m7eX8ldwpOJgMLrgOWhJcgIWgKNiOeNDgyFDp4kz/3bdRz0BC3WpcHXz8viOO5YxeOVgMU+vjuF9/xoVwy0ol1xvPGlTMaFhHcUmsDWVkyBylBA12KllmRZwfG5FQ8nRohIAFLHpsHMezndwMVNq/ZGGNglY151pzcAVprWeubT60v7aD9vOr+mU+0QRAWsLjqjdoThvb8FtLRdTBR4rZgCFaNAxQJrMb7SCiAwtbPzMUK4XcAKfyhIFIrtxdLAEKDS8wKKYOWWFUHj5lPQKsHIwytV+3WdxX1LsPLzqfOpLQIWeavMK729fb/p6PjHZNmy6HmDFapYqTQF9F+64sqZZ55Zteqxx/rO4ZX0b2pq+vGuJJn0m85OjTXlxsRSktuHge6OIPNyZWQ+DFTtQjhq/8xhFdMS1Neusynz6ZCw9B8COKmUjsu8QkuZoClTOKi2CFptdFD9k+OPjz81ZcqicMWK5x/kTT93rFrFGiumQOUo4K/HivnEVwKrm4HVO+fNWzy9puaW6bW1c6rh0N27dkXbGUFBOShBRMLUY+OwhuJ0HbWsFqtJGfsXukKwr4eVwCXzsBKoXgdBIFUCFutLuTCmfWdWeXFttOuyM89M3jF/frBt797w5b17P/G2u+66SdBasmpVnrYW6Me2VkyBsaxA6R/+WP6Qxc8WXnrppRnB6o8XLPhoS0PDL2c0NMzhdfO5cdXV0duamvpd9SKAYEEo5jwceTkKzbxpnazfPADyYZ5qZ34bP5+u/TpqHd+bzumtFVhd/ra3Be8444wwzxt3pjY15U+eOvXG//zAB/7xwlWrcnyBRIjLKul7RB0rlapARXhYy3BuVhSiueCSE0+8bnpd3RebuPgJAeMaeo6qy0I93tGvduwIHtuzJ6ine4JCPhUJJK8q7WXJw0p7WdXMy7NSLa8qbYN5WGxSKgKjCy+p0+Ggestv7ukJPnLeecFFZ58d5Pftc94azy7SfyIX9+bzmed2777te889dwUQ7rodGC+94w4dwoopMGYVGPPAWkzYtApPpAXunHvKKT+YUV//gfpMJq98Febe8OwetQFaeqbv1ra2YBteDW9dLkFL375CQA+twYDloSVweWAJVm8ELJ8TE2lkCgWVE1N4+ocLFgRXv/WtQQxcw9paXkfN0WhncViaPqarnt+164H71qy5/L8/8EBbArRCgxYKWhmrCowaYNFfiWvxjmjq1Kmva/O2bdu4Vi+lD2coh6VUlFx/TMn1OXPmL5g06SfTGxrOA0Q5vKqsclXu2UABgDt4GilBDzXvwJO5a/v2oIdp9VJPH1AgUl5rsDxW2sPStECl7dPAUsO96bjePKyUaBesdpNgnzdpUnD1mWcGTXV1wYTGxqC6AVzWgE29iFWF9jJqaR9tr3pp165Xnti27aIP3nbbK7+/8sqqs26+WYeyYgqMOQVed/GX0ycUpJYvX57B6EzeH0aDtZPtfV47vuqss7IuuX7iiYunjxt3y7T6+jmMoJADWFn3uAygUifOhAu/NLQLoAiZf62rK1hFaCjvR8BJQ0uCydsStFTL25J3lfaw3ghYbO5K2rvyg/TtpU011dXB+/CumsmrTQFW5K2CSdQ1ghaeloMWYCtCiw8QZ15rb297fvv2y//0X/7lAZJaWe4k2jOGXmirx4wCxX/XZfd5QhLJPNkSuut4xYoVwXXXXTd97ty5J9XW1i6qqqqaUl1dPZf6BSDVgaPx/I4dO55n+47UJ+n741NO+ei0ceNunFJfXydYEQI6/8RlqPGgNMSLvBt3EkChYV26AdZ4Qq8z8GaextvSOtAgtohVjpo9TCh0E7D2FeuhAovNS0UHTXtXOtY+8mdvBkxbu7sDBulzj/vE1JqeROvrq8ChN8CG15XB08rPq66eHtXU/OKBj3/8o8DqtuT22zPh0qU6vBVTYMwooIuwrIrueAEoXcvBlVdeOf6MM864rKmp6XLszZMnTx4/YcKEoI4wCWC5sdS1XWdnZ7Bv374t3d3dvyMP9b8ampvv//wll/z1zAkTrmtiO2DlkusuPANMCgHd0C7AKYdXpciql3ofdTfW2dsb9FJvo16fJHmAkRHcKPrT746cZuRRCV5paPlw0IeEjnZso9qHgh5YxG8JsIr36ji5XGYRwGrGo5qF4R06myZvC5uMNbA8RANEoDW0QG3rA7d9fdFuEvUv7djxybdef/13H8TTWjJE75RTWzEFyl6BsgKWT5Cfd955dZdccsknZsyYcc3s2bPnjh8/3kEKUMWyGrwfAQuOaMQFPY8nLoRato3807e+/OVtO15+eSoXdyKYAKxCct3DClDlAZRgJVAJTj2CFYDqwvZiHbKennxbT08mV1v7Guedp2+Tc8lrkW4lcAlAmvHwEqy8DQQWqxywBCu8N/Wsj/HWMvKgdnV0BH10YF3A513Q1BS2AKuZwGmGoAWoHLRYJ3CNY1lYT0ZN0GJfl4zv7XV3Dmh/9vm2tu+e+/Wvf1L6xNdeG+F16ZRWTIFRrYCuq7Io6iNFUj330Y9+9M3HH3/8vxx33HGnT5w4McCzyjU2NoYNDQ0RdQQ4Ig+rdMPxmJJ9PT3x9/7u78Jtq1dPnT55ckyuSncCHUh0J1BekvNqFP4VQ8B9wKsH6wZcXdhewAWskp0cq62rK9O+b9/tSVfXFS0tLWdxjuvpBHEm0BK4xCnBK4JeChVLrwoTzTzR0nWxvRrPRhazTxYPiSHi4zztuxlPcQ1DzVy/RuDi+UaAFrlwFS9KI0E4A06Cm/pkqWtGRkGuQkTV1dVhkstlG6uq8m+eO/eaJ7761am33HOPQsRu6/ZQVN+qUa1AWQCrCKv8hz/84UsWLFjwE3JV8Kkhh2UAVlbehIx551UNVFxJ8yiTCR//9a8zrz7xRDB39uyEx2Kifsl1+CJYuddnpWAlULkwEFDtJWcFoOId3d3R1s7OzN5c7kuQ6Ws6HyM7rKI6d9asWR+kvhrP5XyMd1MURlxgmeCVCFAUl+sSbIpFcFNxjhj7qXCTMukinLuL5d/ctGnTk9qAsbQ2Ma77D17t6KijzgOujHv0B/DmgVY/A1wTsKyAxTpBK9SdRPJaNb29udPnz/9AzXvfO/Wkk0++fOn3v99meS0pbGU0KzDiwFIYKM/q/e9//6WA6nbCwKC+vj4PoByoBCnBCu+qlLM6kOAKv6ZMnqxwMVQXBXVbUH6nlFwvgopOl8E+QNUjrwpQKWe1R7AiBNyGV4V1EiJ+jMPdhumwYp1gkwcsP6L+EcPSnA9wBK+3YyfieWk7VwrOV8HLcgADNsUVGkq5h7kngOzPqvL5O9Zv3fpqYWWhBwSZqNuy2ewWnlO8Zf2ePXOYzxE2ZgU/EdGNveXBVfS2JgKqau9pCV4KE2tqskrOvWnOnLdnamt//dPPfe4ikvCv/P6f/7nqrKuu4pCFopwhU9GSJUsCTGDVZw1Xrlyp5W7Z8uXLE59X1DIrpsBIKTCiwPJh4Lvf/e6zm5ubfzBlypQ8sCI1U58h9HN5K6YFsIPCCm/F6Tdn0SIHN5/bEqz84HjOs2JesBKolLPy+SrBand3d27r3r3ZHV1dr/Tl8+/jgE9h0sc7SuncVbx58+bfsE4W0fbjqRUqvgmbii2iTQ0sU75qO17aBiY3sOxxbgq82NbWtlbrikWgEygcj6iryK2tIk+3hPqXmzs7T3DQShIHLe9tKTx03lYRWrqDWJuGljytfL4q6O3NH9/SclwN0Pq3ZcsuB1YPCFqvTpwYF/uuCVAxQKIqFQ+u/QuKXUb4DNpe662YAsdcgZEEVnjSSScleFW1JNX/hQ6hDVykeeWolFT3JnDhcRxUGL1OS17NjNmzg5bjjgtaX31VXlZhnPSiVwWESl6VYCWvSiHgnp6eZFdXV75tz57sru7uX5FP+ggna8PSsPLn92DRvEAjUuZaW1tfpJYNtWg/7a+LX6BKF3k/Ge56rqU+D8/tzm3d3W/v27ZN9wGrXE6LFQJXP28LDSYDr3oPLdUFTysjaM2dNWs6ntbPH/jGN64AWvIcXaG7yIl4tWfRXeQteLFT0Pp4gJ/Dq32BDdrQ9Yndu3f/HlCtK+4irdXlZGC7/WqrTYGjpsDBSXDUThsE/HdX94U8IeFXgNVJACYHpEjHVDlAqRZ0VA+lCFjksYLF73lP8KN/+AfXSz3tVeluoPequgQqrKO7O97Z2RkKVoDruxzjU5xLUBJMvGd1oNP7C1bwkbkQqrhx2gtJr9OxZVp/sOPr2GrDbgB6MQD5Lt0VPvritm15oBWR1wrTwEp7W1PQoAFYuR5n+8NDQSueOX16feO4cbc+ctNNzXe3tvbOnzXrL8c1NZ2Fh5glV+g8WWku47t4G7B0d1I7Ojo6X3zxxYe5KfAjnh64FVghreuky5D2b9yhl89hxRQ4IgqMFLAi8lZ5EtgMmtDw3/jvzuNyWa7LjEuq+1qelaaHUnRxcXEHJ7zlLcEZF14Y/OZXvwpqCSWVq3L5qmII2Amo9tIps727O799z55MW0dHL90Z/jvn+H8wQUeA8TBi8g1LGkKDbewBNdi6gy1TG9QeBirNX4EOu+mu8OmXtm+P9bJViFHoVYun5eDF51eIqE6mk7FxaBcJWDKmE0akgD6JPLD1bW03nLBwYdA8fTod52uVH8zzPSTyZgEVefsa/dNIpD1ACvGAYWD2Ivq5XfTaa699dvv27V9n8e1qvMDFtABsxRQ46gqMCLDwqiIeSI7ptvAh/rPrYshxcfCyZN6WzAUoS08PVQXtpwT7n/3X/xq8tmZN8PIzzwTVQEu5KnlVnXSq7OCxG7ypuK29PbOto2M9ntdHOP5KTGRMe0bMjnhRewTQDND6DDqtwUv87prt20OgFQOtQrcHNnDhIZ9f0HIeF7AZL+ALVoIWENKd1F+8+GKwg460cydNcoAi9FN3kYxubihXKIDJw5L+qeKgq/Vsfzq5xttefvnli+6+++5r0Lz7dnrVL7Ve9Sm5bPJoKTAiwFq5ciXXF32kqqvfgyms0EXp8lAK7WTylvx0cfUbaqDtdGVVc3F+5G/+Jrjxq18NnnnssSDLRai7gR1Ai+R68srGjSH9rR5m0w9jazHpcLAQjdUjVvSR1LYsSfibALtuCtz12o4dddSFbg/oJWDJ0/LhYVz0ttTtQdBSz/j7APgmoD2nuTmqB1CClO6+yjywgOJgH9RBUx4XRZ13k4ULF17Bnd2TTjnllMsvuuiitXxX5mkNppwtO6IK9Ps3ekSPfOCD6Yed8J/6OC6AM9iMZ5DVbaoAKYHKG9ejg9aBD/X6NQ5aHGvytGnB57/5zeCtF10UbGLImFfWrw8AVfDU2rXJBe98Z3jzzTf/DXuvJfFPH4CyhVX6Azpoocmv8H7eTreHtRt27cq80Nqae5Xe/esYy2vjzp3B5l27nLXu3h3gRQZ4kUFIWPzw008Hv9+wIZjOKBDcMfRPDpTuxApYB4BVug2a5vSuC0ffnDlzzqOT7/0///nPZ6E7N2RdXmvg9jZvChwxBQb9d3rEjj74gQTJmNDjzQCrih8512BeIY8bj0qPy8j0fJ+34n/2wY82yFIPLV2Yn16+PPiDd787uOvOO4P/xNv6i3PPjd//vvdFjz/++AXs+uiNN94YX0jOa5QU72n9ljzTO9DnLsD0ZjytHMn4rG4yaAOZ97iqCO+efvbZ4D+3bQtaZs4MNK8clTeFgMpdHarGnII4M+ibN2/e/J6enlu//e1vv1On5vtEfkvEo4WVo6DASADLfQz+S8+QVwWceJqkP6B6yTnJuLUfcDG4i0sQOpTioaV9zmLUzjPOOSfYgSfSRQ5rA54GF/txrFKOaCS8TDVruEU88t0edCfvTm4eXAS0+uhjVpWGFhQKeBYyeIaHw2eQYNcdV2/yppSrkg31TuwgDeaFQ3EfXtYF3EFcgeZf4DtV3HgoNy0GOawtMgUGV2DELlZ+3Cd6D8oDykPKg0qwEmBUD6d4aHFbPmjbsiVoJ0TaunVrqPPi0c3RMZcsWSJojbbi7yDy2sP4j/mcN7V3dVW9tHFjfvXmzclrW7cGmwgRX+IzP0wYLG+qGljJi/ImYMkEK+k03AIws3KRp02b9rmHHnrozcznlYQf7vFsP1PgYAqMmIfFhVaXBlQaUtw+L4Us3iPgQnDLDvZhBq6TB6fjClDy4nQ+hZ4CIOeYUdxeSe3RWLx3qJj6E8BnE102rlu9aVOAt5XgbYV1fP4ZLS0OVuqjJg0H2jBCwYFa0W83CuhLl9myZcs/ovk76GM3WjUd+NlsvswUGDEPCx1eFUwEJwEkbVomz0rjXO3duzfYw+if8pK0XBAaShGY/HHS50hNdw/lOGW+jfcOMwD5a4Dj43R76H2ltTVs5W4ot2ED+oqUPkLak9K0t9IGw5wAejp/TH+tJffcc89pHNcS8MPU0nY7uAIj5mHxo94kYAlGGu9KgFICWOGLvCqFK/rvr//eKrpzKC/J9xXSNn6d/4iCmd9Ox5YVvSl3/CLAEtWse1n7LWcIZirlhUZrEcEVIuoZxO9Rr6bLwp10mm3KSDuEprjPJn28SSc/7VYe5h++i5gBFrM8pvQhDvV/rSw8PO2BephHt91NgYICIwEs5yIBpBe5wBI8qIy8qGIva5cETsNKF5surjSIfOdGQSsd0sirEtTSuTHvuXmPTV4buawEUL46xn4ELhnPYza/Rb8eADIBIKk47aRN2qSn5lUfiSIq6lh8H3+kESCWLFliifcjIawdo58CIwEsd4XgSb0ARHiUr3My4Z56XbsRQz2EvPekC85faMpBeS9Md7c82LwH4S9CbS9o+RyZoCVgyZvDIh4tCbHfSYnnn39+aDFmP9nKckauKB89/zY01B1YF5alNfE6qhbYPdyl6eEWzhNJb763ExjTrIXvZL3AZcPSHK6ytn9agZEAlgCRefXVV9sZU+pXgOVyRgPQKA2lB58FK0Eo7R3o4tIFIfgcKGxMw81fjD4sFLAERryriKFhdr7yyisOWDzTeGRcjLSqIzPt4j48nPk6PQCJ0SCSDt7Snqe09ECXdyv4D7cUdad/ak7DV9fx/S3iWOtPPvnk4d9+HG5jbL8xrcDwf6VHQBYuru9xXV2O18N7QmsH9ZhY7zwsXWyCTzocTHtYacBpH12kfh/lrBQKcp78rl27snhX//q73/1uB880upesHoGPUk6HOFEA0Wf3GnhQqU6DysNc+ujxHO+pHuqH8efBc+OJIPo15PMLOcb93Dk0YB2qmLb9QRUYKWC5fkR04HwIL2sVLVy8c+dOhinPkicuJNnT3pK/0PwF5js7DgSWPqkHnMIe7eeT7sqT4clF69evzzFq6D9pW/IsMQ9ha3LMFKAzTtp5MHkvStqlzevitVQoPpzQUHrrWP580pw2jBWvdcz8LsbKBxkpYEk/99+XH/cXmH6YH3yC51N4rqOYt9J/bpm/6HRBKRz0F5mAJcDJvHfggeX/6+tiKoaDOXlXAOsnv/3tb5/RaKcaj2usfJGpz6HcVQlO+vwDQeU1TOcB/T+KQ4GWzqPje/P/UPDYrONo6guxySOnwEgCS7DI4O38Fi/rm8Dn8/zwGVhzmx736Jd38RfCYMDSXULBajBg+f3wrmJyV8qbbWcgus8Vk8FjJdk+8NfgumvoswvUHiaqFfoJUgNB5WHlvVrpLF0PVuTB+mPquDLOF8qTJfzeq33p2nCwQ9g6U+CQFRhJYKmxCh0y9K1axo/8LC6Sd3ho+RBDF54uBt//ynsHCmF8SOgvOHfAVP5KFxQXUUzuKgGMSrZfxqM5bXhWijvHZNgCuNdKB0KzUMDyMBGE0p6pgCTz3qlg5f9RSHNt7+/Ypv8ZpLfx/xB0HuUIsYjQXqG3hlceS3dg9XGslIECIw0seTnxmjVrennv38VcND8XtIBVH6+ezwKcUBccQ9E4byHdV8sDSxecB9bAi47989wZzAApPUN4NTmz+8dool0/JQdgNHiGaY2akC1C5HU3Knwo7UHkdZPXpH8U3pNNa6wTaDtt48NtAcv/Q8GzUofckO9tG2H3K9qeZwpjfw7NWzEFDleBkQaW2i9oRbz3r5tOj5cw/VOgdYn+k+MZ5YBOVmGGH2BOHpbM//eXl+AvCl1QxYspZj95VlkuoE6g9UneVPMDjj0W7wrysVwRsEI0XEtv8zVA/ERgIg2itHclvbxX5XWT1v4mhc8Xpr0xbSdLgy19I0T/VPiO8lgG7+rxH/7wh7sVdrPPmPRii3pbNQIKlAOw9LH1w4640Lqo/5yL7gv82JdxcdUKQOSfclwMGmtcb9RxY44LWDJ/8XEx8ao/HvrFQ+BCzXARRdhjeBkf4yJ6iuPqs47mR3Bo/huWzGOPPdZHTvBf0e+L6BGjmwNW2ltKaVa6q+q9Junn/yGk9xkILL+9PCyFhPxTCPmeQm5s3Fls5ZgNu9/wW7ANjpoC7k7dUTv6oR9Y7ZHFhIin8p//80wv5QLTQH/u4mI+YV79fVTrP7/+k+styqWQhQv1Vey7hIHqvuCS+8WaakwXBwle7rEIPZ7jk2ok10jeKeOwB4yh757bVIjth0X2nqvCbR8KpoGV9siknLwxmYDlvTH+KegRq2Dt2rUdzzzzzIJ77713J5vqe5T3bMUUOGIKlIuH5T+QfuCyDCGicjEf4p1513HRXQKUNO7T6dh4rfcXjgcVy9ZijwIx/Yf/N2Alb01FF7GgVQlFnqruvL6Ml3onsLkMcOfwgLJ4ma/L9UlD1vsw2uWj0uGjv6khaKG708/rrpBQpnAQYClXmCXsvlGwUpcRvRWpEgS3z3hsFSj8Co/tOYd6NoFG7Sv98Ok5PYP//nOA1GyWT8YEN93G38nFtXrdunU9TPuSYUIXsLappCLdEsZbn4cX9BSgqZeXRR3q3YPyslTLy5J3JZN3lfaw0vmrwTwsD7liwt2FnWj/2lNPPXXKe97znm7uwkrzStO9kn5jI/ZZy83DSgsh2KjoApTF9NHaQi17FBusCFIq2rcEOrekcv44L4s7dWvJZf03OPX/ASzevZpU4QW9LmeVTp77kHAwYHn5vFemcJAbGwneW4w3q5smVz399NOdJ5xwgv9H4Xex2hQ4YgqUs4c12IdUewUvFd92/9/c14W19lfgyAOt7xAmXwNo+pjXiyNcn7a0l6U+br7LiPJXaWApT1jMFZZyhMpf4V0lAEtPD1TRXeRLjzzyyNfGcJcRyWalDBTwF30ZNMWacIQV8HDPk8+6kdDuaqClu6TOY9UdQIWDCg09sNIelvJXPoeFl+aahpfmk+3qNhKQF1Nn3BvxrK4xWB3hb88ON6gCBqxBZRkzC0vQ4s7h9Xyqz+iTAR51LNWD5m4MMgHLQ8t7WB5YA7wrveEoT6I9S492DdXzWTr93lBMsldivnDM/FBGywcxYI2Wb2r47SxBi/DwA3hL/wObLW8J432sie4Ahnhcgpfr4+ZDQrwydR1hEzq5cUeRfJdeL6ZOos9yZ/DT5K7+j8Fq+F+M7XnoChiwDl2z0bqHy2nRv20S/PkidhkwmqUPw3SpjxvscjkreVaYWFZaz93Bddh3eTvOjezWbd0XRutPYfS224A1er+74bTcQUs7zp07dwIe00XA6r1A6WwWzaHWegeoYq1uIq3Yo6y7ixzXv2ukWK2jlI5VmLW/psDRV8CAdfQ1Lrcz6DtX4r3U7WPhwoU1dE+Yh0c1Pt1YQsKthIqt5Kn2pZYLVJavSglik8dOAQPWsdO63M7kwaV2leB1gEY6z4t1BqoDCGSLj40CBqxjo3O5n0W/A2/ptvq+baqtmAKmgClgCpgCpoApYAqYAqaAKWAKmAKmgClgCpgCpoApYAqYAqaAKWAKmAKmgClgCpgCpoApYAqYAqaAKWAKmAKmgClgCpgCpoApYAqYAqaAKWAKmAKmgClgCpgCpoApYAqYAqaAKWAKmAKmgClgCpgCpoApYAqYAqaAKWAKmAKmgClgCpgCpoApYAqYAqaAKWAKmAKmgClgCpgCpoApYAqYAqaAKWAKmAKmgClgCpgCpoApYAqYAqaAKWAKmAKmgClgCpgCpoApYAqYAqaAKWAKmAKmgClgCpgCpoApYAqYAqaAKWAKmAKmgClgCpgCpoApYAocNQWSJAk5uMyKKWAKmAJlpUB46aWXZpYtW5YFVFERVq6BmtayBx98MHv77bdnWGgQK6uvzhpjChwdBcJlQRBdGgSZ2zGmszJNy5KRAUEkUA32ca+44opxixcvrh1sHQDLpKE22DbltswoW27fiLWn7BQQoGiULF6BDaWBgMsDJOYiY/boFGCUXbVqVU5HP+200xrOO++8P5g+ffoFTU1NZ2CTq6ur59XV1fU0NDRsyWQyO7HHwzC872//9m8feeyxx/q0n8DFsrymy70YsMr9G7L2jZQCIV5UhNfUDziXTp3aODubbZmWyTRxpbfUBcGUKI6DbBStzoZhVxxFG65av7413egivJRMGhLs0vseZFohXgBokjPPPHPO8ccf/6lp06a9v7m5ee6MGTMCYBWMGzcuqK2tDerr6wOgFdTU1ATZbDbo6ekJent7n+/r6/vBmjVr/ufSpUvbOZaAnOh4BznniK8yYI34V2ANKDcFFO7dEQQlj+PDs2efPDGX+7OmMPzD+iA4HZvOdKaBhteHYVBF7eBBnUuSvVxUr3DhP9mXJL/ckyT3/0Vr63Z9xiMILuft6ZhvfetbPz137twvz58/f9LkyZODCRMmxOPHj48BVtjY2ChQRVgiq6qqCjBBKQJeEXWwY8eO9Xv27PnCKaec8lPXRnJdlLKFlgFL35IVU6CgQMSVKk9IF2zmslmzPjgxCK6ekCTnTI+izDg8GiWDaqhr2IbpuJoaU/wXYSHwCqsAAd5W0MvyvUmyTeDqDIIbLtq8+UkWOXBxjhIQtewQiodV3dlnn33rokWLLpk5c2YwZcqU3MSJEzNY6L0rAYtQ0HlX3sMStKJIhwhiIBsTImbxtoLNmzf/9L777vurq666qkveFtA6kt7gIXy8g29qwDq4Pra2QhRIe1XvnjXrQ1OD4LMzwvD0ScCpHmvEecKzCgkBI0zgCoFWIFgBKWdZalCg2EoWQz2wFWQIF4NdcZzbFwS37Yiib/zRxo3PaN1ybMWhhYkiTQyIpuBV/Zww8Dxg1QesslgIrAK8q1I4KGDJFBLKBKtBisAky7a2tt7/2muvXXz++ef3aLty9LQMWPpmrFS0Ah5WJ7W0LJwbx9+aFUV/OgVI4VHlx3HdAquQ8C8EWEERVgVPCxwNBJY8K1GFBJPrNwAJElyZGHcqUw24tsdxN17XV9/S2voNNktuZ/nSoXlbulYjclREczW/WLBgwYUtLS195K2qpk6dGghWA4GlHJaH1gFgxSFLRQn4qi1bttz/0EMP/Ql3HeUBxuUGLWlrxRSoWAUW41koX3XSzJmXtMTx7wSrCXhDjUAGUClPFaVhJWDVAaW01QOitDUw35jJBI3U46IobCKcJOcVZDnuZHZfkMl8/cVZs+69u7l5jmD1IG0YwhegazUfx/ENeFEXAqJeEupVSqpzJ9CZoOTNL9P6IcBKp5f71UfC/h3nnnvuNwGVgFV2fCi7Bkk5K6bAsVBAsFpFqHfizJlXzwzDnxECTmiM4zygyuJNudDPe1QK/2qBTg3mazcNlGplAKoOU12DVWHZYl1NrXVALDuOhDbg6psXhu88KQxX3j9r1psvpA3ytA7ymbUuT3eFJXhXVxLe5ahhUrW760ceKkib7gTKPLQOctx+q8hdCVo5vLhPrl69erGgxbKDtavf/sdixoB1LFS2c5SdAh5W84HV1DC8kRAw34BXBagyHlIuwU7LFfYJTqVa4AJS3gQpWRWQiKihRcGUM/LTLA+1TSYT4oFVcVsxByTnLwyCXz/Y3HyBPK3kwNBSjkl9pW6QtwSsQtUeUkqiDzStE7SKCfYh6c/xtZ2OneDB/RM94vWRde6yKQassvkqrCHHUIGMPCs8iSsnRtGNE4EVkIrkVaUhJXfDGRdyqRa4AIRMnlQNUJAHFXkwCVIHMm3DtjJgl2XL3PQgaFgQhrc+NG3aAnCRX/b6MEweTkK+6oJMFJ0OFOPqKrDH+WWCTBE0Tj4/L1AJWul1Q9RX54vx5uiHetrb2T/R4zxD3Peob2bAOuoS2wnKTAFdkPnJM2ac3YhnNZEQsBj+ve6unyClJHqpBgJV3oCBQCUL07AiTCMW628DAVYEFyFjlsbkZ4XhrBnZ7M/IadWfrHw9xqlDngmMFs+dW5VQk8T/ED2kErLgcY6VfUUjZmNqfyGEK80MA1Z+3wTgqe/Wx7RgyZIl+w/qtxihumzIOUKf305bWQq4mGf83LkTavv6fjo+DLN1hIHkp0LFPg5M1LooRDVisMK0am8AS90UBCrBS2Ge95qIv/CPMLZ1xjEI/fBXiKpkefLYOXDjt2FePlcun++bH0WndOfzN5wWBFf9/swzqwIem1mxYoVA0ROuWBFMmTbtXfXV1eHU2tqInvbBLFaoy4UAtY9jJhyfhLyb1zJvasIwij4+Hy2z5J577uF+QdjB8ahGvkOpAWsY36btMmoV0IWYq+3t/XJDJnNcPXftgFVWsNKF4M3BinlXAx+By4dggpWS6TKFXA5WqVCvH7Q4Rj9geVip9gW4ZEl299EW7lBe+fvm5jvPeuyxe8+ZNKnp9PnzF06qqzuLW4F/2BeGLbMnTQqmTZgQTVCHUKyOvlXVdAxdw/nb+vqCcQAwP8D8aQ6xFthdfy/6d53B9Ko77rhD0ZjuHI5oMWCNqPx28mOogLvguG3/JrpxX11DKAisMgO9Km3kDEilpx20iuDy0Eon1PtNew/KfzjvXXnAaf2AwrAJric9ubEffOOCC14cX1d30rj6+hlT1REUKNXwLGC1uiioLpqfXsj8I93dQeu+fUEjMOQZQRy5nDN5XYeSeE81K6ZLBLtGi1i2ir5ezjtNrR+RSQPWiMhuJx0BBZzXQKh2LaFcHcDKAatQF4A8qRKc0tMASnkgAFcwQOM8LcDj7gYKQPKuBpqAlIZSGlhazvFKBQ9L4SJgCHPU0/v6Zr67tnZmNGeO6xpRW1eXw8MKMrW1WYFK01lq5p3xdDN3C+qDc3nY+edbtwY9QKsOYOlxm31MC17cVSydbqgTAl3Rpg51n2OxnQHrWKhs5xhpBcSj/LwZM+buC4L/wugK6nAk3PQDlTDiTKBKT2u+CK2wCK1+uSsl1QUtX3tgeTANBBbH7hcqaj3GUdR1ImnYsCFuOuGEiJ5WJNnqstkinCLqUPChhkIFI8FPhtw9L7iAMPEVPK1GPDLBSqMyyNQfS+0/lKLQUtDr6uo64VD2O9rbGrCOtsJ2/HJQQMCKyQNdhmdUqx7n8o20UDbopawLHHPgUi0DDAqvVDsPyntYqj2sVGte2+gYKh5YWqYir6roWZXWsU2I1VVVhdGuXZmq9vagacGCwnH8XUcdO318fywdH0+qmXUvsl8P0KpiH99xVH229CzhoRQPO7ys1Yey39He1oB1tBW245eDAu4xEzDxXj3jB05CoWNQUA1srYeO6qL1A5bg5M2Hhn7e7yugeIB5UOHBuDuGflvNs40uSMLVYM+GDcGEE0+EazyTqG3Sx/L7sq27O0i+KgFY1SxX3QmwsgBLoPIdTAVaPaYzlCJY4VkFe/fuDRh6ZvNQ9jlW2xiwjpXSdp6RUkBsimfNmrWQi/8MAYAFmTSs1HdARfWgBiy4qe/MQ8sBSBCSeej4WuAaCJmBwBF0tJ3uGPrjsI3apc6o2zdvDvKAJ0O+qtQdgnUur8S+ecCU4xx97CvLYzxUHfSyTyf7q9e97mKmO48qzNMwM4LXYEXaKJQEUrLM9u3bg+7u7he17bZt27xMg+16zJYZsI6Z1HaiEVJAV2fMxXoqj6pUccHL28JlKZTBAIU/5J5HoZOmhitw0w5YLH8duAQibwLBQNNptF4FIJTgNnA7baNlQEWP1MRAowtgNDY3u8R5CVBAKMd2vTL28aaOpLvxzNoxPXitbhfyBBXKCkQCne4cKi/lH4j24NJ6v67oWSXUIZDaxnAzL6jpjN4gWUa8GLBG/CuwBhwLBbg4zyyeR9h4nSdVghTrDjjNxe/X9WuzYHMg04YCkc6qbTSdNr9MtQq1S/ADrh3c9cvTraFLXRQAUBpS+9iOGwhBD8eV0Y8raKVuJ/GuO5i6s+kT7YKVvCvdMVS4p7uGym957ysNMzwqhYN5PK3srl27/uMLX/hCu96yw7EE+hEvBqwR/wqsAUdZAfFJZX6hKsDKg2eota5WbatapoMWEVOAEfMOSgKTL5r2IPLLBtZ+veritI4A1hywOqdNC3oA1j4g5wDF8h4Zx+4GRLIuYNQOaNqBm4a9UQ98XzysvAclWMmU25KH5T0wf1ewmL8KgZWGT77JH6dcagNWuXwT1o6jpYA4o7BIOSyBxo39m4aPh9BBa/ZlvHYHK/VTZ1A+3unF0bwBjn7Tgo/AMdh6LTtA0Rq1TR5SR2dn0Ao4evCMunksBkjFWLI3n8/szeXCdvJNu/CYdpAcZ1zkYDJJ9TyhpECqzzoQVspPeVj5ZLyaoW1THliesDGzfv36VZ/5zGce0POMvKRC0pRFMWCVxddgjTgGCpQocTAwCUZ+vQDljGVarmnGZw/UwVPTGUEqbXg6DlJFT8nBSh+Mbftt5/fRcm/ajuLPreNwty9+aceOuLuvL7Mnnw87qHeSg9oOyHYAqg71k8L7asCzmsWIo+oKoeOlYeVBpNyVDwcFK+XJ3KNFnNMDi20SbYNnFTNc8hfUnueff178K5tiwCqbr8IacpQV6AcsByBOOLAWMAYu86ASrJQrUt1LTQY/0Cu+3F08DysfjglE6WltB1ycaVu/n2ptWzR3bqaV3O/s6opWrV8fyZPa2dOzE2A909nX91RSXf3X5KkYYyaT8PqbSFDavXOnG5Y5DSuFgcpbCVbeu1LuygMrnXQvbtcHsKpJtF/3rW9967fKXZWTd6XfhwFLKlipBAVKv3VBSXfVXgem4nI/dIvWO0gBkHTtYSVg1XpQCUZpz0ogGggsbeuhpVrbpExDxRTPE3OHL9y9b98L923b9uNcT8+v2fpZmrMDC+bMmfMcYPpnQj6qJAZI0e7duzlUIbnuQZWGle4M+mS79648sABewj59HKeaUPD2m2+++cuEgllgJQnKqpS+xLJqlTXGFDhyCriUDodr0yE144HlwZSudYXKHDhUA6Fe6ipg0osJUvSUD7LAxz1XyHq91ut1sOJOnVumdeznwJSGlZ8WxDiePCMdvwjSOM/QN3Orqlbi8XxNfTDU7mtB4B04Gc+vX38zb8thxJvwxuKdQHiTywCt0EPKe1W666e+Vx5Y8q5kPuEO5PLsG+GBVW/duvW2Bx544GPKWzG0jWQqu2LAKruvxBp0hBXQLTNeyBytL94REz4chASiNKw07eEkCPlpwUswcaBiuXJXMmKxUvcBdyHpyDLWuTBRsPLFLxOoCNNKnlYRWAo197GvTOfqYV9ewrqH3bNX03vrOzRtBUdmXs3K8h7Bm4DWS0x/j8/mbigAnxwdPiN5XOpPpbfm+PcSClg+HHShZBTJPeNpoDjT3t6esN/XeC/hlzheAKzUcMlUdsWAVXZfiTXoaChA2MOwUft/7h5OuvrTprGx/Do9jOzoAEB0R7BkRVCJhB5JdUzLC3OwknelcNADq7hcj9E4kBW9q1CwwvqwnhSw6LYQdnK87iR5mir3Xv4ArDRA5ARmgNb9hIfn8Nm+BHz+ChCNl6emfBXQynd2diYClZ4j9B4W25D+yqjpGSXY2e4+PLO/e+GFFx6ic2iGca8ExfS5mC2fsv8bLJ82WUtMgSOpgLv48K4e1cVMiQQZXfHq1+SNsQ9K4BKosmyriwP0FIx551UBIV3tgo2HFbMupKsFOi4RLxAJVkVg6ax6sDn0YaCv2W4fyx2sqF2/Ks7TTRt3kptqj6LndOxtgwNEIZu6H+yi/hzjfMG04Ersz/G4TubzZhQeygCXuyMIqGgSnV/jmEMGD7Lu+xs2bLiXafVkF6zKMgxU+3xJa+6XWW0KjCUF9BtPeKkCL6uJVnPBzmReb2WO5E1NxJqK1kgtayhaPRd3PdPuHYR4TJqvo5a513r5GhC4N+gUa72gQkMqCw4qeuBaoeIucked6iuFp0eyKajfsSOo6ugQoNwwx11Ai/fEx1i0Lkk2PJ7JnHDDxo3d7K0jUQ1adBIx1MMm4rnJU4HSKXzeNwHpqUVQc7MzeRZ7AYA9y4PN23U05sPly5czCvMKeVZlX8zDKvuvyBp4mAroQs+0tbV18pac+7mIP8TFHHOVR/KyBAt5V+o9Ls/KeVfUzrsCMm7wPuYjTVPLUyoVlung6kSq5/vUP0sPIqunuRsDnnWMGR9sAFSv8MhMjnGqatXLnK4FenFF78yZwfjVq4PxGzcGXcwTAspierWHnXH874LVMppCW9XUAxU1QbBS82S5TZs2PUUtO2ApelSCqvbVMUZFMWCNiq/JGnmYCnDNy00Jb6H6S0wXtrtKBSpBKw0rHwZ6AghWOoA7CLXLU6kqGq9jdqMlCFq6q+iGUKauZf0qQLQWQE3E+5rAdgoZ3Z1Gjilw7Tr11KCTELGW4WQ6md+Ld9XBqXYnyR3sHpw8dJj45vimus+oYwwojrGjIfwb0G43a8AaTBVbNtYUkIcSknx+iHDoObyskwiF5Flk+vjThSk81MXgPCtqQcsBS7ACPi6sY5kvOGnOs5J3RReEQi94lmn0BI3/zqtmgl+w8ROAag5AUo/0DPBSaOgCvKIn1siyrjPPDHaRKO9ra8vvoyPolnz+uZUNDf9R3FLtPJQiIDkoHcpOo2XbA1F4tLTf2mkKDFWBzJo1axjkIPx7rOQsaWd5WZ2YPK1BDchouXvgWNNFI9cUKO9E+BZ0kkCX8QgNN/7ywa8B0L8KQuSqdgGrdoxnAN222kfJfnljGt5Yd/EyixcHe3jQeXdfX9hRXf21f6etywvcZEsrXgEDllfC6rGugLysaNGiRf8boDwNs+REOe9FsZS8LG8Hg5Zg5UHlkuTAx9WCEaaRE3YDp1t4CHkbwNKzfzuBloDVgQlofj9BS2GkPK5xjY355ne9K7Nt4sRHr1237id6LGbFwXNX7F15xYBVed95JX/icNWqVQwtlflUUQSFTq4oNNyLydOSCV4Cl6/dtIfVYB4Wy/YAL+WzngZUv+buXwe1oKVwbzfgagdcPA8Y7AFcgpa6MaiTqKAFNEM6eiZve9/7xt9+zTXz9Qzf7XQ1oAlWUgoYsFJi2OSYV0AeVXbjxo2ryGH9A7kspazEKlfk8aSh5cElaKWn5WG5edWYmwc+ezGFjS8CqD0Aqo3e5tt4NGY78zsFriK0BCyGh3HemIMWx2eI44he7sn4hoYTZkyceP/3P/7xBUvpF2XQcl9N6Y++MCumQCUp4JLt5I2u5fGVs+HEhSTQFS66a0GelC/p7LW/BaedtbFGcFCtPJQe3aliWv/9q5lW2KiiURaULnMjgJJ8j0iwO6ObQ+iNENF3laihBzoAy02eOHF+Lpv92fWf/vQ5S2+4oTvh2b7lOmDx7ubJJ5+c+CGL9UZm/5LTJUuWqJkM9DDyr5R3rT0KfwxYR0FUO2RZK+A4pAQ8/bKWAqtHucDn49z0g5aHlWpBSiRwoErVcs0EqGpqAUsXk+4ezqN7AvBxINMAe/419wKX+l9FRVi5Gk9LPeBdT3gdC6+vm6eRp02adMoZVVW/AEx/Rq9Oz1E1Y2BR8/oVPotCSYFrsO37bTvaZgxYo+0bs/YeCQV0IUcMUredHvDvIKelDqXzgY0YJPa4O4fayHtUvvbQcrBivWAl04WkHRUejgNKf0RH0V8woB7HDXYQFuqxGAFLXRtKnhbTztMCcB5aHMJBi24QudkzZrz9M3/wBz+/oLHxg6va2yfTg/1sHmY+o7GxcSojLizS3UWGjNnIdBvPST7NHcnfkaN7HlDpxmegxP1zzz2XjJZe7GrzG5V+t3ffaGNbbwqMMQXcnUI8rTlc5LcAl8XF8FDRncvvagM9ntOA1WHqDCqrKZr3rgQrmbbXjnqM598A1uM8x6eOpALXDEZPmMnooNMYzngqNon5SYBtPHUTveDHYRoKRm9u1vzmPXviNTt3Rs/v3Lm9Zty4CS0MKdPE/hqFwT/MrO39cDEaoYH2r8V+xgPNt+CdPUFT9PhN2bxEQu05nGLAOhz1bN+xoICDVktLSx0X+ncAy8e4wHWRy5nSOtBT8KLcc4VMe2h57yrtYWkHtxO1njn8N0YCXYOH5Xq/42E1AxzZNGyKxmEHPhOxCZhANWXCBNc36zfr1gWrGZSvfuLEeF5LSzSR5XhWeSyRZ4WFfsgYho1J5MEBXd6fWu2GP965c2eOURtu401dKy6++OKX+Tzi6KjPbxmw+BatVLwCuphdvocxppYy/feAa1EKXBG5LHXcCuVZeWBpOu1hKSz0wGLSTXPQ+D927oy38mIHPZfIS1LDmcBnBuagBbgmAy55Wc2TJweteGUr164Nalg/f86cYNqUKcG4pqakiW2AVYi5jqbyrDSC6CCjiOpzwN7YwYvx2fdgXzn//PO/rTYVB+dzn1Xzo60YsEbbN2btPVoK6FoQuPKTJk1qwnu5Eo/lGmxuEVx63kWpLLkqIeGfG+3Be1kKB4vAioGWcvUiYFYvM1VXhyd5KWovyXW5bvU1NaFeGjEdE7QmYzOB1VpGb3gOmzN3bjBn1qxgAusnAC7BSmGgYKUB+eRhpUcRFbQUFmq8L4WetFmnYYzAXI5OsoyUwwPYGzbcduedd15BPqtrNEPLgKVflhVTYL8CcpLcnbcFCxaMJ6y6jPmPcNGfAwy0zsHAEQkQ6ALyXpU8KE07YFBrjHa6PzzH26N/tLe9/eH2rq7bWDeTY8WNdXXRrEmTHLRmASsl5tfzuq7jFy4MGNsqAJrBRIA1HpilgSVYpUcR9bkshYIeVpy6VABWQjI+B9CqXn311fuffPLJi8lt9WgD2uI+RmnjUTBhwBoFX5I18Zgr4DmkPJYrJOZP5AI/l5nzAc6pbDCf6Vp/xbuaoVpY/hrLX+Sh6Eeh3q+3t7Y+ybwLwfCALmDfXzFfSx001ddHc3l+MIN3pOcNFx13nIPVZAAmYMm7ErC8d+U9rIFelmAlSB6s4Gz1keeqojvH/T/+8Y//ZPny5YIyLB1d0Dr4pzyYArbOFBj7Cuj68B6XZ5P71HTWxEmqUyRYKgApxotpLy0oTlzKMe4ohJuCxh+x+F5Mb54IJzY1hbWEes2MjdVCGMhxAwFroHfloSVYeZO3pXBwqIW8FsN1RVXPPvvsd0499dRPAc1Rd/fQ+mEN9du27SpRAUHKe1nKb8lUYu6+7S1Mvu6vIKftVMuz0sBW8mZkvGwnfx/Q+ASezY24RTFvxwiVj9IdQh/eKR/lTXmptGm5PCpvHHPIRbBi49y8efM++cQTT/yCNtw72qDlv4Ahf2jb0BSoUAUEH8FLpmkBaTAT5AQnv13aM9Myxu+LbwJcXxKggFVetUww8nBSNwVvWqZpv071UMJAzjVYiQgtE0LOf7r77rvVUwNmJfoco6IYsEbF12SNLEMFBKLB7I2aKmgpzPw64NgAiDKEdbGHkcAkUwI9bX65thPYND/Moms+P3v27BN4487H8bLilStXDvtgw2zDsHczYA1bOtvRFBiWAg4YPGazGOi0CBgASK/eKgGKZS6JPrD2APPbDuvshZ24oRkm5MX+msd3qpcsWSKPcFQUA9ao+JqskWNIARd+EYb9BdCAPZnYA0iAUvGg8tO+9ssFrsMsOkAybdq0hXTdWMJxEz13eJjHPCa7H/YnPyattJOYAmNHgdzixYt5qU74h0BLXpW8nX6fTsvTll7poZVeNszpWLkswtGLtT99s/o3YpgHPdq7GbCOtsJ2fFNgvwLuelu7du1sgDRHUKLwUp7+gCIp74ap8csHzhf323/U4U3xXo2E0W6yF9DzXb0FRkVYaMAa3pdte5kCw1HAeTHcIZyGp6RHEmFGAVYeSqxzsNK8ptPzfpsjBKxIx8bBm3f66ac30h61pey9LAPWcH52to8pMDwFHBAAz8JiGMjkfjB5OHlQDazprV6C2PBO328vhuDKJST8G0m+z9EajV7ab4synCn7BpahZtYkU+CwFMCTGVc8QCIoeRCpHorp1WHa73CLjgE4s0BT/bFGRTFgjYqvyRo5lhQAEup46oo8LAEobYKW5n2dXtfLiyxkWnY4RWGlPy7tOXz6HU5jDmFfA9YhiGWbmgKHqYDLsgOIVwQM6kjASkNIEPFAStd+G9WMIBH0MFa89h1uEQz98TnnqOjSoM9qwBruN277mQKHroADFrttAFhykdwdwjSMDjTtISVQyboZjkb1cIv21bn27t3bu2vXrt06Dl0bfPuGe9ijvp8B66hLbCcwBUoKOJeIRPd6vKNWPCytiAUODyLVgtOBTOsFK1kn48Vr/lCLjs347wnnDffs2bPz8ccf31A8RtkDy0ZrONRv27Y3BQ5Pgcy6det6GIr59wBrNhaT/I68xzQQVAKSnh2U6TlC9Yof2DNe4aUeni4C8KCt0/HxqgQ8gTID9J6+/vrreR9sog6sw48xD3rWI7fSgHXktLQjmQJDUcC5VQBCA/m9F1MHqJLXJHBpjCuNxpAGVRpWApbgJNO+PnHvR3wYCC6/TdGzcp4ZoEr0lp329vYH1OiVK1cq2jJgSQwrpoApUFLA3ZEDOv8bz+obwGU8lpAEDxXiCTp+rCsPrLRXpecIPZAEInVNkCmsFJC0r4ebzqht/N1GrZfHpvMQCmba2tp6AdZd2m5J4a3RmizrYh5WWX891rgxqIDyRNmNGzfuZNjlHwCgT+MhaciZrDweD6x0CChIeUvDSp6VYOXv+AlGApYHnLTzUPPbyIPjPHoFWHbLli0rP/nJT65mm1ERDurzGLCkghVT4NgqoNBL7xX8BuD4ENOTsRioRHg+pXyV95S8V3UgWKkrhLwnH0ZqP+2jImD5kNEn9wFWyKu/gu3bt39d2yxdurTsH8lRO1UMWAUd7K8pcCwVELAyPATdhpf1JTyim+VlAaRIXlJHR0cprBsKrASidBgpD2sgsORhCWryrijZ9evX//grX/nKSg0rA7BGTcdRA9ax/JnauUyB/QoIEtnW1tb/CbTeCWSWAi31zapSjsmDym/uPaV0GCgI+dyVEvUDw0jtmw4JgWGe7bObN2/exMsyPqv3Ez733HNl35XBa6DagJVWw6ZNgWOrgKCl3u5X4F2dBKROEbSATJVCw4EhYBpWvke8cl7pu4o+f6V9ZX4fQBXjYWV4hX0X+bPLvv/977fRUTTDA8+jxrvSV2PAOrY/UDubKZBWQN5NyN26Tkb//HM8pHuA1glAhvevJlmFhgKON3lUApWHlfeufDio3JXPX3nYFZPyOXlWwKqXc73/lltueVhjYPEWaCX7R1UZNcm2UaWqNdYUODQFXB+ouXPnTgBGtwKti4AUzEr0otOMXgOmF6r6F6nqfYSytHclaKWT9Dq9jgGw8spZAavdJNrfz+vq7x+tsNJnMmBJBSumwMgr4DtuhvSC/3tA9XnBSt4V4NKr5iNeoBrphaqDAcvnr4Bdwn7y3NSDPqtEO32tHqYLw1898sgjL43GMDD91Riw0mrYtCkwsgr46zEBWqcDnr+jOX9KzRtXXSd0vWEnlsdVhFao/BWWFL2rEGBlAJwLG7kjuI5+V98GVN/iOMloh5W+Gi+Qpq2YAqZAeSig4V5cMpw7iG+BV1di72bZXN88JdcFKXlWCgc1rwLYevCsniAM/OFLL710K4va2Te49tprI3JWZf/ojfsQB/ljwDqIOLbKFBhBBRQiKrSTBYCrHiidBZBOwYM6mUVN2CKmOwHSy0xvZ/oJPKrHCQHXMu/KWPCq/Gex2hQwBcpfAYGr4D4Nsa2AK9SrxNh8zDkkY+4DDfE7tc1MgdGmgK5Vwctfs/K8FOL55UyWPLJRH/rpw1gxBUwBU8AUMAVMAVPAFDAFTAFTwBQwBUwBU8AUMAVMAVPAFDAFTAFTwBQwBUwBU8AUMAVMAVPAFDAFTAFTwBQwBUwBU8AUMAVMAVPAFDAFTAFTwBQwBUwBU8AUMAVMAVPAFDAFTAFTwBQwBUwBU8AUMAVMAVPAFDAFTAFTwBQwBUwBU8AUMAVMAVPAFDAFTAFTwBQwBUwBU8AUMAVMAVPAFDAFTAFTwBQwBUwBU8AUMAVMAVPAFDAFTAFTwBQwBUwBU8AUMAVMAVPAFDAFTAFTwBQwBUwBU8AUMAVMAVPAFDAFTAFTwBQwBUwBU8AUMAVMAVPAFDAFTAFTwBQwBUwBU8AUMAVMAVPAFDAFTAFTwBQwBUwBU8AUMAVMAVPAFDAFTAFTwBQwBUwBU8AUMAVMAVPAFDAFTAFTwBQwBUwBU8AUMAVMAVPAFDAFTAFTwBQwBUwBU8AUMAVMAVPAFDAFTAFTwBQwBUwBU8AUMAVMAVPAFDAFTAFTwBQwBUwBU8AUMAVMAVPAFDAFTAFTwBQwBUwBU8AUMAVMAVPAFDAFTAFTwBQwBUwBU8AUMAVMAVPAFDAFTAFTwBQwBUwBU8AUMAVMAVNgLCjw/wOQ5MRPPC+KaAAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "execution_count": 2, "metadata": {}, "output_type": "execute_result" } ], "source": [ "molecule = Molecule()\n", "\n", "molecule.add_atom(Atom(symbol=\"O\", coords=[0.000, 0.000, 0.000]))\n", "molecule.add_atom(Atom(symbol=\"H\", coords=[0.558, 0.598, -0.513]))\n", "molecule.add_atom(Atom(symbol=\"H\", coords=[0.028, -0.815, -0.512]))\n", "molecule.add_atom(Atom(symbol=\"O\", coords=[-0.066, 0.0230, 2.724]))\n", "molecule.add_atom(Atom(symbol=\"H\", coords=[0.000, 0.000, 1.749]))\n", "molecule.add_atom(Atom(symbol=\"H\", coords=[-0.956, 0.357, 2.838]))\n", "\n", "view(molecule, width=300, height=300, guess_bonds=True, direction=\"tilt_y\")" ] }, { "cell_type": "markdown", "id": "a808bf5f-f91f-44ad-b5a9-284e6e45208c", "metadata": {}, "source": [ "### Run MD" ] }, { "cell_type": "markdown", "id": "nvt-md", "metadata": {}, "source": [ "Run NVT equilibration at 300 K with GFN1-xTB.\n" ] }, { "cell_type": "code", "execution_count": 3, "id": "nvt-code", "metadata": { "tags": [] }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[31.03|11:57:32] JOB plamsjob STARTED\n", "[31.03|11:57:32] JOB plamsjob RUNNING\n", "[31.03|11:59:56] JOB plamsjob FINISHED\n", "[31.03|11:59:56] JOB plamsjob SUCCESSFUL\n" ] } ], "source": [ "engine_settings = Settings()\n", "engine_settings.input.DFTB.Model = \"GFN1-xTB\"\n", "\n", "nvt_job = AMSNVTJob(\n", " molecule=molecule,\n", " settings=engine_settings,\n", " nsteps=20000,\n", " timestep=0.5,\n", " thermostat=\"Berendsen\",\n", ")\n", "nvt_results = nvt_job.run();" ] }, { "cell_type": "markdown", "id": "nve-md", "metadata": {}, "source": [ "Restart from the equilibrated trajectory and run NVE production while logging dipole moments.\n" ] }, { "cell_type": "code", "execution_count": 4, "id": "nve-code", "metadata": { "tags": [] }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[31.03|11:59:56] JOB plamsjob STARTED\n", "[31.03|11:59:56] Renaming job plamsjob to plamsjob.002\n", "[31.03|11:59:56] JOB plamsjob.002 RUNNING\n", "[31.03|12:02:03] JOB plamsjob.002 FINISHED\n", "[31.03|12:02:03] JOB plamsjob.002 SUCCESSFUL\n" ] } ], "source": [ "nve_job = AMSNVEJob.restart_from(nvt_job, binlog_dipolemoment=True)\n", "nve_results = nve_job.run();" ] }, { "cell_type": "markdown", "id": "analysis-md", "metadata": {}, "source": [ "Set up autocorrelation analysis of dipole-moment derivatives to obtain the IR spectrum.\n" ] }, { "cell_type": "markdown", "id": "dbba03bf-3847-4866-aad4-aeb7d6698925", "metadata": {}, "source": [ "### Run AMS Analysis" ] }, { "cell_type": "code", "execution_count": 10, "id": "analysis-code", "metadata": { "tags": [] }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[31.03|12:13:07] JOB plamsjob STARTED\n", "[31.03|12:13:07] Renaming job plamsjob to plamsjob.004\n", "[31.03|12:13:07] JOB plamsjob.004 RUNNING\n", "[31.03|12:13:08] JOB plamsjob.004 FINISHED\n", "[31.03|12:13:08] JOB plamsjob.004 SUCCESSFUL\n" ] } ], "source": [ "analysis_settings = Settings()\n", "analysis_settings.input.TrajectoryInfo.Trajectory.KFFilename = nve_results.rkfpath()\n", "analysis_settings.input.Task = \"AutoCorrelation\"\n", "analysis_settings.input.AutoCorrelation.Property = \"DipoleMomentFromBinLog\"\n", "analysis_settings.input.AutoCorrelation.UseTimeDerivative.Enabled = \"Yes\"\n", "analysis_settings.input.AutoCorrelation.WritePropertyToKF = \"Yes\"\n", "analysis_settings.input.AutoCorrelation.MaxCorrelationTime = 2000\n", "\n", "analysis_job = AMSAnalysisJob(settings=analysis_settings)\n", "analysis_results = analysis_job.run();" ] }, { "cell_type": "markdown", "id": "a081e520-41c3-4187-85a4-e1eadc99397e", "metadata": {}, "source": [ "### Extract Results" ] }, { "cell_type": "markdown", "id": "write-md", "metadata": {}, "source": [ "Write all generated analysis plots to text files and list the outputs.\n" ] }, { "cell_type": "code", "execution_count": 11, "id": "write-code", "metadata": { "tags": [] }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Wrote analysis files:\n", "- AutoCorrelation_1.txt\n", "- AutoCorrelationRaw_1.txt\n", "- Spectrum_1.txt\n", "- Integral_1.txt\n" ] } ], "source": [ "plots = analysis_results.get_all_plots()\n", "written_files = []\n", "for plot in plots:\n", " filename = f\"{plot.section}.txt\"\n", " plot.write(filename)\n", " written_files.append(filename)\n", "\n", "print(\"Wrote analysis files:\")\n", "for filename in written_files:\n", " print(f\"- {filename}\")" ] }, { "cell_type": "markdown", "id": "31cd629a-602f-458d-9bb7-b1cff7185ef2", "metadata": {}, "source": [ "Finally, plot the autocorrelation function and IR spectrum." ] }, { "cell_type": "code", "execution_count": 12, "id": "a0e2a870-a58f-47a3-ad97-7079f283d033", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkIAAAHHCAYAAABTMjf2AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/P9b71AAAACXBIWXMAAA9hAAAPYQGoP6dpAAB2UElEQVR4nO3dd3xT1fsH8E+SNunei0Kh7L0sUstGyhRcoKjIUnEAghRQcYCKghMQRVAU8KsgOAD9yVD2kLIpu4zSUqCLUjrobnJ/f7S5zW2TNm3T3qb9vF+vKr259+bcpEmePOc55ygEQRBAREREVA8p5W4AERERkVwYCBEREVG9xUCIiIiI6i0GQkRERFRvMRAiIiKieouBEBEREdVbDISIiIio3mIgRERERPUWAyEiIiKqtxgIERFZSGBgICZMmGDRc06YMAGBgYEWPScRFWMgRESib775BgqFAsHBwVU6z9atW/Hee+9VuT1arRarV69Gv3794OHhAY1Gg8DAQEycOBHHjx+v8vlri7i4OLz33nuIiIiQuylE9Q4DISISrV27FoGBgTh69CiuXr1a6fNs3boV77//fpXakp2djeHDh+O5556DIAh46623sHz5cowbNw7h4eHo3r07bt68WaX7qC3i4uLw/vvvGw2EVq5ciUuXLtV8o4jqCRu5G0BEtUN0dDQOHTqEjRs34qWXXsLatWsxb9482doze/ZsbN++HYsXL8Zrr70muW3evHlYvHixRe4nMzMTjo6ORm/LysqCg4ODRe6nsmxtbWW9f6K6jhkhIgJQmA1yd3fHQw89hFGjRmHt2rWS2/fu3QuFQoG9e/dKtsfExEChUGDNmjUACmtali1bBgBQKBTij15mZiZmzpyJgIAAaDQatG7dGp9//jkEQRD3uXnzJr799lsMHDiwVBAEACqVCrNmzUKjRo3EbadOncLQoUPh4uICJycnDBgwAIcPH5Yct2bNGigUCuzbtw+TJ0+Gj4+PeI5+/fqhQ4cOOHHiBPr06QMHBwe89dZbAIDc3FzMmzcPLVq0gEajQUBAAF5//XXk5uaW+ZimpKRg1qxZ6NixI5ycnODi4oKhQ4fi9OnTksf1/vvvBwBMnDhRfLwMH8+SNULmPIb6x3/q1KnYvHkzOnToAI1Gg/bt22P79u1ltpuoPmFGiIgAFAZCjz/+ONRqNZ5++mksX74cx44dEz+kzfXSSy8hLi4OO3bswE8//SS5TRAEPPzww9izZw+ef/55dOnSBf/88w9mz56NW7duiVmebdu2oaCgAGPHjjXrPs+fP4/evXvDxcUFr7/+OmxtbfHtt9+iX79+2LdvX6map8mTJ8Pb2xtz585FZmamuP3OnTsYOnQonnrqKTz77LPw9fWFTqfDww8/jIMHD+LFF19E27ZtcfbsWSxevBiXL1/G5s2bTbbr2rVr2Lx5M5544gk0bdoUiYmJ+Pbbb9G3b19cuHAB/v7+aNu2LT744APMnTsXL774Inr37g0A6NGjh9FzmvsY6h08eBAbN27E5MmT4ezsjKVLl2LkyJGIjY2Fp6enWY8vUZ0mEFG9d/z4cQGAsGPHDkEQBEGn0wmNGjUSpk+fLu6zZ88eAYCwZ88eybHR0dECAGH16tXitilTpgjG3l42b94sABA+/PBDyfZRo0YJCoVCuHr1qiAIgjBjxgwBgHDq1Cmz2v/oo48KarVaiIqKErfFxcUJzs7OQp8+fcRtq1evFgAIvXr1EgoKCiTn6Nu3rwBAWLFihWT7Tz/9JCiVSuHAgQOS7StWrBAACP/995+4rUmTJsL48ePF33NycgStVis5Ljo6WtBoNMIHH3wgbjt27Fipx1Bv/PjxQpMmTcTfzX0MBUEQAAhqtVqy7fTp0wIA4auvvip1X0T1EbvGiAhr166Fr68v+vfvD6CwS2X06NFYv349tFqtxe5n69atUKlUmDZtmmT7zJkzIQgCtm3bBgBIT08HADg7O5d7Tq1Wi3///RePPvoomjVrJm5v0KABnnnmGRw8eFA8n96kSZOgUqlKnUuj0WDixImSbb/99hvatm2LNm3aIDk5Wfx58MEHAQB79uwx2TaNRgOlUim2886dO3ByckLr1q1x8uTJcq/NGHMfQ73Q0FA0b95c/L1Tp05wcXHBtWvXKnX/RHUNAyGiek6r1WL9+vXo378/oqOjcfXqVVy9ehXBwcFITEzErl27LHZf169fh7+/f6kAp23btuLtAODi4gIAyMjIKPect2/fRlZWFlq3bl3qtrZt20Kn0+HGjRuS7U2bNjV6roYNG0KtVku2XblyBefPn4e3t7fkp1WrVgCApKQkk23T6XRYvHgxWrZsCY1GAy8vL3h7e+PMmTNIS0sr99qMMfcx1GvcuHGpc7i7u+Pu3buVun+iuoY1QkT13O7duxEfH4/169dj/fr1pW5fu3YtBg0aJCl4NmTJjJFemzZtAABnz55Fly5dLH5+e3t7s7frdDp07NgRixYtMnpMQECAyftZsGAB3n33XTz33HOYP38+PDw8oFQq8dprr0Gn01Wu8RVkLPMFoFRhNVF9xUCIqJ5bu3YtfHx8xJFehjZu3IhNmzZhxYoVcHd3BwCkpqZK9imZgQBgMmhq0qQJdu7ciYyMDElGIzIyUrwdAIYOHQqVSoWff/653IJpb29vODg4GJ1rJzIyEkqlssxgpTzNmzfH6dOnMWDAAJPXZcrvv/+O/v3744cffpBsT01NhZeXl/h7Rc5r7mNIROZh1xhRPZadnY2NGzdi+PDhGDVqVKmfqVOnIiMjA3/99ReaNGkClUqF/fv3S87xzTfflDqvfl6ekkHTsGHDoNVq8fXXX0u2L168GAqFAkOHDgVQmGWZNGkS/v33X3z11Velzq/T6fDFF1/g5s2bUKlUGDRoEP7880/ExMSI+yQmJmLdunXo1auX2NVWGU8++SRu3bqFlStXlrotOztbMuqsJJVKVSrz8ttvv+HWrVuSbaYeL2PMfQyJyDzMCBHVY3/99RcyMjLw8MMPG739gQcegLe3N9auXYvRo0fjiSeewFdffQWFQoHmzZvj77//NlojExQUBACYNm0aBg8eDJVKhaeeegojRoxA//798fbbbyMmJgadO3fGv//+iz///BOvvfaapKj3iy++QFRUFKZNmyYGa+7u7oiNjcVvv/2GyMhIPPXUUwCADz/8EDt27ECvXr0wefJk2NjY4Ntvv0Vubi4+/fTTKj1GY8eOxa+//oqXX34Ze/bsQc+ePaHVahEZGYlff/0V//zzD7p162b02OHDh+ODDz7AxIkT0aNHD5w9exZr166VFHUDhVknNzc3rFixAs7OznB0dERwcLDRWqaKPIZEZAZ5B60RkZxGjBgh2NnZCZmZmSb3mTBhgmBrayskJycLt2/fFkaOHCk4ODgI7u7uwksvvSScO3eu1NDvgoIC4dVXXxW8vb0FhUIhGUqfkZEhzJgxQ/D39xdsbW2Fli1bCp999pmg0+lK3XdBQYHw/fffC7179xZcXV0FW1tboUmTJsLEiRNLDa0/efKkMHjwYMHJyUlwcHAQ+vfvLxw6dEiyj374/LFjx0rdV9++fYX27dsbfQzy8vKETz75RGjfvr2g0WgEd3d3ISgoSHj//feFtLQ0cT9jw+dnzpwpNGjQQLC3txd69uwphIeHC3379hX69u0ruY8///xTaNeunWBjYyN5PEsOn6/IYwhAmDJlSqnrKdlOovpMIQismCMiIqL6iTVCREREVG8xECIiIqJ6i4EQERER1VsMhIiIiKjeYiBERERE9RYDISIiIqq3OKFiOXQ6HeLi4uDs7Fzh6fWJiIhIHoIgICMjA/7+/lAqTed9GAiVIy4urkrrFBEREZF8bty4gUaNGpm8nYFQOfSLGt64caNK6xURERFRzUlPT0dAQIBkcWJjGAiVQ98d5uLiwkCIiIjIypRX1sJiaSIiIqq3GAgRERFRvcVAiIiIiOotBkJERERUbzEQIiIionqLgRARERHVWwyEiIiIqN5iIERERET1FgMhIiIiqrcYCBEREVG9ZVWB0P79+zFixAj4+/tDoVBg8+bN5R6zd+9e3HfffdBoNGjRogXWrFlT7e0kIiIi62BVgVBmZiY6d+6MZcuWmbV/dHQ0HnroIfTv3x8RERF47bXX8MILL+Cff/6p5pYSERGRNbCqRVeHDh2KoUOHmr3/ihUr0LRpU3zxxRcAgLZt2+LgwYNYvHgxBg8eXF3NNMvdzDzcyy2Au6MaThqrehqIiIjqDKvKCFVUeHg4QkNDJdsGDx6M8PBwk8fk5uYiPT1d8lMdpqw7id6f7sGui4nVcn4iIiIqX50OhBISEuDr6yvZ5uvri/T0dGRnZxs9ZuHChXB1dRV/AgICqqVtdrYqAEBOvrZazk9ERETlq9OBUGXMmTMHaWlp4s+NGzeq5X7sbAsf+px8XbWcn4iIiMpXp4tT/Pz8kJgo7XpKTEyEi4sL7O3tjR6j0Wig0WiqvW12NswIERERya1OZ4RCQkKwa9cuybYdO3YgJCREphYVU9sUPvT5WmaEiIiI5GJVgdC9e/cQERGBiIgIAIXD4yMiIhAbGwugsFtr3Lhx4v4vv/wyrl27htdffx2RkZH45ptv8Ouvv2LGjBlyNF9CqVQAAHSCzA0hIiKqx6wqEDp+/Di6du2Krl27AgDCwsLQtWtXzJ07FwAQHx8vBkUA0LRpU2zZsgU7duxA586d8cUXX+D777+Xfeg8ABTFQdAyEiIiIpKNVdUI9evXD4JgOnAwNmt0v379cOrUqWpsVeWoFIWRUFnXQ0RERNXLqjJCdYmiKBDSMhAiIiKSDQMhmSgVrBEiIiKSGwMhmaiKHnkdM0JERESyYSAkEzEjxJQQERGRbBgIyYTD54mIiOTHQEgmHD5PREQkPwZCMuHweSIiIvkxEJIJh88TERHJj4GQTFSsESIiIpIdAyGZ6GuEOGqMiIhIPgyEZFI8aoyBEBERkVwYCMlEP4+QVidzQ4iIiOoxBkIy0XeNcdQYERGRfBgIyaR4rTEGQkRERHJhICQTsWuMcRAREZFsGAjJRMViaSIiItkxEJIJh88TERHJj4GQTDh8noiISH4MhGTC4fNERETyYyAkExVHjREREcmOgZBMFJxHiIiISHYMhGSiUHDRVSIiIrkxEJKJOLO0vM0gIiKq1xgIyURfLM2uMSIiIvkwEJKJvkaIxdJERETyYSAkE7FGiMPniYiIZMNASCbFNULMCBEREcmFgZBMlBw1RkREJDsGQjJRch4hIiIi2TEQkg0zQkRERHJjICQTZoSIiIjkx0BIJqwRIiIikh8DIZkoix55ZoSIiIjkw0BIJgrWCBEREcmOgZBMFJxHiIiISHYMhGSi5MzSREREsmMgJJPiYmlmhIiIiOTCQEgmYtcY4yAiIiLZWF0gtGzZMgQGBsLOzg7BwcE4evRomfsvWbIErVu3hr29PQICAjBjxgzk5OTUUGtNY40QERGR/KwqENqwYQPCwsIwb948nDx5Ep07d8bgwYORlJRkdP9169bhzTffxLx583Dx4kX88MMP2LBhA956660abnlpnEeIiIhIflYVCC1atAiTJk3CxIkT0a5dO6xYsQIODg5YtWqV0f0PHTqEnj174plnnkFgYCAGDRqEp59+utwsUk1gjRAREZH8rCYQysvLw4kTJxAaGipuUyqVCA0NRXh4uNFjevTogRMnToiBz7Vr17B161YMGzasRtpcFtYIERERyc9G7gaYKzk5GVqtFr6+vpLtvr6+iIyMNHrMM888g+TkZPTq1QuCIKCgoAAvv/xymV1jubm5yM3NFX9PT0+3zAWUwLXGiIiI5Gc1GaHK2Lt3LxYsWIBvvvkGJ0+exMaNG7FlyxbMnz/f5DELFy6Eq6ur+BMQEFAtbVOwRoiIiEh2VpMR8vLygkqlQmJiomR7YmIi/Pz8jB7z7rvvYuzYsXjhhRcAAB07dkRmZiZefPFFvP3221AqS8eBc+bMQVhYmPh7enp6tQRDrBEiIiKSn9VkhNRqNYKCgrBr1y5xm06nw65duxASEmL0mKysrFLBjkqlAmC6S0qj0cDFxUXyUx2KesZYI0RERCQjq8kIAUBYWBjGjx+Pbt26oXv37liyZAkyMzMxceJEAMC4cePQsGFDLFy4EAAwYsQILFq0CF27dkVwcDCuXr2Kd999FyNGjBADIrnoM0KsESIiIpKPVQVCo0ePxu3btzF37lwkJCSgS5cu2L59u1hAHRsbK8kAvfPOO1AoFHjnnXdw69YteHt7Y8SIEfjoo4/kugSRftQYa4SIiIjkoxCYkihTeno6XF1dkZaWZtFusgtx6Ri29AB8nDU4+nZo+QcQERGR2cz9/LaaGqG6hhkhIiIi+TEQkom+Rghca4yIiEg2DIRkomRGiIiISHYMhGSi4DxCREREsmMgJBOxRogpISIiItkwEJKJOI+QzO0gIiKqzxgIyUTJ1eeJiIhkx0BIJgqwRoiIiEhuDIRkUjyPEAMhIiIiuTAQkolSqV9rTOaGEBER1WMMhGTCGiEiIiL5MRCSCWuEiIiI5MdASCZK1ggRERHJjoGQTBScR4iIiEh2DIRkYlgjJDArREREJAsGQjJRiKvPs2CaiIhILgyEZKIsjoNYJ0RERCQTBkIykWSEZGwHERFRfcZASCbMCBEREcmPgZBMWCNEREQkPwZCMjHMCDEQIiIikgcDIZkoDTJC7BojIiKSBwMhmShYI0RERCQ7BkIy0a81BgA6xkFERESyYCAkE8MaIY6fJyIikgcDIZmwRoiIiEh+DIRkwhohIiIi+TEQkolCwRohIiIiuTEQkpG4Aj2LhIiIiGTBQEhG+joh9owRERHJg4GQjPSBEGuEiIiI5MFASE5FXWOsESIiIpIHAyEZiTVCzAgRERHJgoGQjFgjREREJC8GQjJijRAREZG8GAjJSD+TEGuEiIiI5MFASEYK1ggRERHJioGQjJRKfdeYzA0hIiKqpxgIyai4WJqREBERkRysLhBatmwZAgMDYWdnh+DgYBw9erTM/VNTUzFlyhQ0aNAAGo0GrVq1wtatW2uotWVjjRAREZG8bORuQEVs2LABYWFhWLFiBYKDg7FkyRIMHjwYly5dgo+PT6n98/LyMHDgQPj4+OD3339Hw4YNcf36dbi5udV8443QL7zKtcaIiIjkYVWB0KJFizBp0iRMnDgRALBixQps2bIFq1atwptvvllq/1WrViElJQWHDh2Cra0tACAwMLAmm1wm/YSKOp287SAiIqqvrKZrLC8vDydOnEBoaKi4TalUIjQ0FOHh4UaP+euvvxASEoIpU6bA19cXHTp0wIIFC6DVak3eT25uLtLT0yU/1YXzCBEREcnLagKh5ORkaLVa+Pr6Srb7+voiISHB6DHXrl3D77//Dq1Wi61bt+Ldd9/FF198gQ8//NDk/SxcuBCurq7iT0BAgEWvw1Dx8PlquwsiIiIqg9UEQpWh0+ng4+OD7777DkFBQRg9ejTefvttrFixwuQxc+bMQVpamvhz48aNamufkjVCREREsrKaGiEvLy+oVCokJiZKticmJsLPz8/oMQ0aNICtrS1UKpW4rW3btkhISEBeXh7UanWpYzQaDTQajWUbb4KCq88TERHJymoyQmq1GkFBQdi1a5e4TafTYdeuXQgJCTF6TM+ePXH16lXoDKqRL1++jAYNGhgNgmqaPiN0/U6mzC0hIiKqn6wmEAKAsLAwrFy5Ej/++CMuXryIV155BZmZmeIosnHjxmHOnDni/q+88gpSUlIwffp0XL58GVu2bMGCBQswZcoUuS5BQp8Rmr4+AlcSM+RtDBERUT1kNV1jADB69Gjcvn0bc+fORUJCArp06YLt27eLBdSxsbFQKotju4CAAPzzzz+YMWMGOnXqhIYNG2L69Ol444035LoECX1GCAC2nUtAS19nGVtDRERU/ygEru9QpvT0dLi6uiItLQ0uLi4WPfeDX+zFtduF3WKzB7fGlP4tLHp+IiKi+srcz2+r6hqrawwzQlpWTBMREdU4BkIyUhj8m4EQERFRzWMgJCODhBAiE9Lx3JpjuBhffTNZExERkZRVFUvXNYZJoH/OF86PFJ2ciT2z+snTICIionqGGSEZ6Yx0h0Unc04hIiKimsJASEYFrAsiIiKSFQMhGbFAmoiISF4MhGSk4xROREREsmIgJCNmhIiIiOTFQEhGDISIiIjkxUBIRlp2jREREcmKgZCMmBEiIiKSFwMhGRmbR4iIiIhqDgMhGXEeISIiInkxEJIRS4SIiIjkxUBIRiyWJiIikhcDIRmxWJqIiEheDISIiIio3mIgRERERPUWAyEiIiKqt2wqeoBOp8O+fftw4MABXL9+HVlZWfD29kbXrl0RGhqKgICA6mgnERERkcWZnRHKzs7Ghx9+iICAAAwbNgzbtm1DamoqVCoVrl69innz5qFp06YYNmwYDh8+XJ1tJiIiIrIIszNCrVq1QkhICFauXImBAwfC1ta21D7Xr1/HunXr8NRTT+Htt9/GpEmTLNpYIiIiIksyOxD6999/0bZt2zL3adKkCebMmYNZs2YhNja2yo0jIiIiqk5md42VFwQZsrW1RfPmzSvVoPrEVqWQuwlERET1WqVGjW3fvh0HDx4Uf1+2bBm6dOmCZ555Bnfv3rVY4+q6Sb2byd0EIiKieq1SgdDs2bORnp4OADh79ixmzpyJYcOGITo6GmFhYRZtYF0WNrAVQtv6yt0MIiKieqvCw+cBIDo6Gu3atQMA/PHHHxg+fDgWLFiAkydPYtiwYRZtYF1mo1KiZwtP7LyYKHdTiIiI6qVKZYTUajWysrIAADt37sSgQYMAAB4eHmKmiMyjVLBOiIiISC6Vygj16tULYWFh6NmzJ44ePYoNGzYAAC5fvoxGjRpZtIF1nZJxEBERkWwqlRH6+uuvYWNjg99//x3Lly9Hw4YNAQDbtm3DkCFDLNrAuk7BjBAREZFsFIIgCHI3ojZLT0+Hq6sr0tLS4OLiYvHzrzsSi7c2nZVsa+HjhLUvBMPXxc7i90dERFQfmPv5XamMkEqlQlJSUqntd+7cgUqlqswp6y1jXWNXk+5h1cHomm8MERFRPVOpQMhUEik3NxdqtbpKDapvTBVL52uZqCMiIqpuFSqWXrp0KYDCupbvv/8eTk5O4m1arRb79+9HmzZtLNvCOo4lQkRERPKpUCC0ePFiAIUZoRUrVki6wdRqNQIDA7FixQrLtrCO4/B5IiIi+VQoEIqOLqxb6d+/PzZu3Ah3d/dqaVR9ojIxfp7xERERUfWr1DxCe/bssXQ76i0GPERERPIxOxAKCwvD/Pnz4ejoWO56YosWLapyw+oLdo0RERHJx+xRY6dOnUJ+fr74b1M/ERER1dVWAIUr3QcGBsLOzg7BwcE4evSoWcetX78eCoUCjz76aLW2r6IYCBEREcnH7IyQYXeYXF1jGzZsQFhYGFasWIHg4GAsWbIEgwcPxqVLl+Dj42PyuJiYGMyaNQu9e/euwdaax9QSGwyPiIiIql+l5hGSy6JFizBp0iRMnDgR7dq1w4oVK+Dg4IBVq1aZPEar1WLMmDF4//330axZsxpsrXm4xAYREZF8KlUsnZmZiY8//hi7du1CUlISdDqd5PZr165ZpHGG8vLycOLECcyZM0fcplQqERoaivDwcJPHffDBB/Dx8cHzzz+PAwcOlHs/ubm5yM3NFX9PT0+vWsPLwUVXiYiI5FOpQOiFF17Avn37MHbsWDRo0KBGshrJycnQarXw9fWVbPf19UVkZKTRYw4ePIgffvihQnVLCxcuxPvvv1+VplaIqRohJoqIiIiqX6UCoW3btmHLli3o2bOnpdtjMRkZGRg7dixWrlwJLy8vs4+bM2eOZFRceno6AgICqqOJAAClVXVOEhER1S2VCoTc3d3h4eFh6baUycvLCyqVComJiZLtiYmJ8PPzK7V/VFQUYmJiMGLECHGbvgvPxsYGly5dQvPmzUsdp9FooNFoLNx600xl00ws50ZEREQWVKl8xPz58zF37lxkZWVZuj0mqdVqBAUFYdeuXeI2nU6HXbt2ISQkpNT+bdq0wdmzZxERESH+PPzww+jfvz8iIiKqNctTEaa6xgp0jISIiIiqW6UyQl988QWioqLg6+uLwMBA2NraSm4/efKkRRpXUlhYGMaPH49u3bqhe/fuWLJkCTIzMzFx4kQAwLhx49CwYUMsXLgQdnZ26NChg+R4Nzc3ACi1XU4sliYiIpJPpQIhuSYlHD16NG7fvo25c+ciISEBXbp0wfbt28UC6tjYWCitrOjGVEZIx74xIiKiaqcQBH7iliU9PR2urq5IS0uDi4uLxc9/KCoZz6w8Umr7sw80xoePdrT4/REREdUH5n5+W1f6pA4ylRFieEpERFT9KtU1plQqy5w7SKvVVrpB9Y3prrEabggREVE9VKlAaNOmTZLf8/PzcerUKfz44481OhlhXWCqWJo9lkRERNWvUoHQI488UmrbqFGj0L59e2zYsAHPP/98lRtWX3AeISIiIvlYtEbogQcekMzzQ+UzlRHiqDEiIqLqZ7FAKDs7G0uXLkXDhg0tdcp6gTVCRERE8qn0EhuGXTqCICAjIwMODg74+eefLda4+sD0qDFGQkRERNWtUoHQkiVLJL8rlUp4e3sjODgY7u7ulmhXvWFq8B3DICIioupXqUBo/Pjxlm5HvcWZpYmIiOTDCRVlZmpFENYIERERVT8GQjJjRoiIiEg+DIRkZnL1ecZBRERE1Y6BkMxMTajIjBAREVH1YyAkM3aNERERyceigdBbb72F5557zpKnrPNMrzVWs+0gIiKqjyo1fN6UW7du4caNG5Y8ZZ3HmaWJiIjkY9FA6Mcff7Tk6eoFkxMqMiVERERU7VgjJDPWCBEREcmn0oHQrl27MHz4cDRv3hzNmzfH8OHDsXPnTku2rV4wudZYDbeDiIioPqpUIPTNN99gyJAhcHZ2xvTp0zF9+nS4uLhg2LBhWLZsmaXbWKeZKpZmjRAREVH1q1SN0IIFC7B48WJMnTpV3DZt2jT07NkTCxYswJQpUyzWwLrO1DxCrBEiIiKqfpXKCKWmpmLIkCGltg8aNAhpaWlVblR9YpgRauzhgC+f6gKANUJEREQ1oVKB0MMPP4xNmzaV2v7nn39i+PDhVW5UfWKyRohxEBERUbUzu2ts6dKl4r/btWuHjz76CHv37kVISAgA4PDhw/jvv/8wc+ZMy7eyDisZCOl/Z0aIiIio+pkdCC1evFjyu7u7Oy5cuIALFy6I29zc3LBq1Sq88847lmthHacokZMrDoRkaAwREVE9Y3YgFB0dXZ3tqLdKZ4SK/sFAiIiIqNpVuEYoPz8fzZs3x8WLF6ujPfVOyeHz+riIXWNERETVr8KBkK2tLXJycqqjLfVSyYyQgjVCRERENaZSo8amTJmCTz75BAUFBZZuT71TctAYa4SIiIhqTqUmVDx27Bh27dqFf//9Fx07doSjo6Pk9o0bN1qkcfWBYUZIoQBYIkRERFRzKhUIubm5YeTIkZZuS71Uqli6KEfHmaWJiIiqX6UCodWrV1u6HfVW6WJp1ggRERHVlEqvPk+WUXKtMbFGSCdHa4iIiOqXSmWEAOD333/Hr7/+itjYWOTl5UluO3nyZJUbVl+xRoiIiKjmVCojtHTpUkycOBG+vr44deoUunfvDk9PT1y7dg1Dhw61dBvrFX1GiDVCRERE1a9SgdA333yD7777Dl999RXUajVef/117NixA9OmTePq81Wk5ISKRERENaZSgVBsbCx69OgBALC3t0dGRgYAYOzYsfjll18s17p6RiH+h6vPExER1YRKBUJ+fn5ISUkBADRu3BiHDx8GULgeGbt0qoarzxMREdWcSgVCDz74IP766y8AwMSJEzFjxgwMHDgQo0ePxmOPPWbRBtYnCoXCoEZI5sYQERHVA5UKhL777ju8/fbbAAqX21i1ahXatm2LDz74AMuXL7doA0tatmwZAgMDYWdnh+DgYBw9etTkvitXrkTv3r3h7u4Od3d3hIaGlrl/bcAaISIioppTqUBIqVTCxqZ45P1TTz2FpUuX4tVXX4VarbZY40rasGEDwsLCMG/ePJw8eRKdO3fG4MGDkZSUZHT/vXv34umnn8aePXsQHh6OgIAADBo0CLdu3aq2NlaFAsVrjzEMIiIiqn5mB0KxsbEVOnF1BBuLFi3CpEmTMHHiRLRr1w4rVqyAg4MDVq1aZXT/tWvXYvLkyejSpQvatGmD77//HjqdDrt27bJ42yxCwZmliYiIapLZgdD999+Pl156CceOHTO5T1paGlauXIkOHTrgjz/+sEgD9fLy8nDixAmEhoaK25RKJUJDQxEeHm7WObKyspCfnw8PDw+T++Tm5iI9PV3yU5M4szQREVHNMXtm6QsXLuCjjz7CwIEDYWdnh6CgIPj7+8POzg53797FhQsXcP78edx333349NNPMWzYMIs2NDk5GVqtFr6+vpLtvr6+iIyMNOscb7zxBvz9/SXBVEkLFy7E+++/X6W2VpYCBjNLMyNERERU7czOCHl6emLRokWIj4/H119/jZYtWyI5ORlXrlwBAIwZMwYnTpxAeHi4xYMgS/j444+xfv16bNq0CXZ2dib3mzNnDtLS0sSfGzdu1GArDWaWrtF7JSIiqp8qvNaYvb09Ro0ahVGjRlVHe0zy8vKCSqVCYmKiZHtiYiL8/PzKPPbzzz/Hxx9/jJ07d6JTp05l7qvRaKDRaKrc3spQKBRisTRrhIiIiKqf1aw+r1arERQUJCl01hc+h4SEmDzu008/xfz587F9+3Z069atJppaaYajxnSMg4iIiKpdpVefl0NYWBjGjx+Pbt26oXv37liyZAkyMzMxceJEAMC4cePQsGFDLFy4EADwySefYO7cuVi3bh0CAwORkJAAAHBycoKTk5Ns11EWLrpKRERUc6wqEBo9ejRu376NuXPnIiEhAV26dMH27dvFAurY2FgolcVJruXLlyMvL69UN968efPw3nvv1WTTzaJQgDNLExER1SCrCoQAYOrUqZg6darR2/bu3Sv5PSYmpvobZGGcWZqIiKjmWE2NUH2ggII1QkRERDWoUoHQwoULjc7mvGrVKnzyySdVblR9pmCNEBERUY2pVCD07bffok2bNqW2t2/fHitWrKhyo+or1ggRERHVrEoFQgkJCWjQoEGp7d7e3oiPj69yo+oz1ggRERHVnEoFQgEBAfjvv/9Kbf/vv//g7+9f5UbVZwpwZmkiIqKaUqlRY5MmTcJrr72G/Px8PPjggwCAXbt24fXXX8fMmTMt2sD6hDNLExER1axKBUKzZ8/GnTt3MHnyZOTl5QEA7Ozs8MYbb2DOnDkWbWB9oyzqG+OoMSIioupXqUBIoVDgk08+wbvvvouLFy/C3t4eLVu2lG2NrrqCq88TERHVrCpNqOjk5IT777/fUm2p92xVCo4aIyIiqkGcULGWMTZqjNkhIiKi6sFAqBbJ1wrihIr6GqGtZ+PR5YMd2Hf5towtIyIiqpsYCNUCL/VpBgB4+6G24qgxoDATNHntSaRl5+P5Ncdkah0REVHdZXWLrtZFc4a1xZQHW8DFzhYpmXnidsMesQIOIyMiIrI4ZoRqCRc7WwDFNUIA5xIiIiKqbgyEahkFiiMhJoGIiIiqFwOhWkZh8IwIXGiDiIioWjEQqmWUBtXS7BkjIiKqXgyEahnWCBEREdUcBkK1jGGNEOMgIiKi6sVAqJZRMCNERERUYxgI1TKGNUIcNUZERHVRTr4WS3ddwZXEDLmbwkCotlGWmFmaiIiorvkp/DoW7biMCavlXzWBgVAto+CoMSIqR1ZegdxNIKqSHRcSAQC3UrNlbgkDoVqHo8aIqCwbjsWi3dx/8MPBaLmbUq641Gzsv3yb2W0S3UjJwoNf7MXRmBRxW06+VsYWMRCqdRS1oEYoNSsPH/59Aadi78rTAKqSAq0O+Vqd3M2gavLt/msAgPl/X5C5JWXT6QQ8s/Iwxq06iv87Ey93c6iW+P3ETVy7nSnZZrjGphwYCNVC+lhIrpmlv9t/Dd8fjMbbm87Jcv9UeTn5WoQu2odHvv4PWlbbl+v0jVQcvnZH7maUKSdfizM3U8WsSkJajnhbQS0JeM/cTEWvT3Zj+d4ocVvEzVTE3MkCAPx56pZcTaNaICuvAKlZefhufxQuJZQujpY7EOLq87WQUqGAVhBqvEZIEATczcrHieuFmaAL8ek124A65M69XHx/MBoTegTC18Wuxu73atI98cPn5t0sNPF0rLH7tjYZOfkY/V04cgt0OPjGg2joZi93k4ya++c5/Hr8JmaEtoKznQ2y8oq7Ea6nZKG5t5OMrSv0y9FY3LybjU+2R+L5Xk2htlFi08ni4OdE7F3odAKUhn3/VC+cjL2L0d+GI19r+gPtbhYzQlSC/r2ipmuEtpyNx33zd+BIdHHf7c+HryNV5j9SazR9fQSW743C9PWnJNu1OgFRt+9VW81EZm5xEe3drPxquQ9LSM3Kw+u/n8ahq8mS7XGp2Rj+1QG899f5am/DzbvZyMnXQRCAyFoc9P96/CYAYPHOy/igRHfYlLUnkZFTc89zRk4+luy8jDM3UyXbb6UWZ6nO3krFxfh0/HT4urgtNSsfV5IK/+6z8+StB6Ga9dn2S2UGQQBwPOaurHVCDIRqIX2dUE3FQefj0tDrk92Yuu5Uqdve2Xyu1JtvZaRm5WH4VwewbM/VKp+rsgRBwLlbabiXW/0jbg4WfcAfvpYi2b7uaCwGfLEPi3descj9CIKAT7dHYv3RWACQZAsiYu/idkauRe7H0hbtuIxfj9/Eq7+ckgSFf5y4iXO30rHmUAySMnLKOEPlJKXniG+4ienF57+RkmXx+6oJkQkZWP1fTI3d38+HY7Fk5xW88vNJyfbLBt0dI5eH4/sDhYXcvVt6oWcLTwCFWaPJa0+i8/v/sv6wFrP0l7RsMwKcL3ddwf7Lty16vxXBQKgW0iePayojtGzPVdy8a3oI48aTVe/f/7/TcTh3Kx2f/XOpyucyl04n4MT1u2IdxdazCRj+1UG88fuZGmtDSQu2XAQALN1lmUBo3+Xb+GZvFN7ceBaCICDTYFj1e/93AU+sOFQrR+ycjyvMwNzJzMOcjWeRkZOPpIwcbDlbXFRrrJagKqJu30PPT3ZjytrCD/EkgyAxPt3yQZcl6Myo8zoUlVzuPpZysShzdis1Gzn5WuRrdfh0eyQSSjx+f5wszGINaueLRzo3BACsORSDbecSkKfV4bFvDuHJb8PNur7aKCdfi9d/P43fjt+QuykWcTUpAy/+7zgm/e842s39ByctGKiqzOwO9XTSWOw+K4qBUC2krOGMUExy+d+Gd0cmVukbumE3TU2lQJfuvoKRyw9h2Z7CAs71xwqzJoYftjXN2LejAq2u0sFKWnbx45qdr5VkhAAg5k4WUmthF5lhceT6Yzfw8bZI9P9sLyINgp+opHsWvc+/IuKQrxWwKzIJWp0gyZYlptXOQCjFRLd0SDNPLHy8I4DCoCQjJ79GiuMNv5y1eXc7pv1yCt8YFEiXrLMKae6FJ+8PELNCho5Gp+Dqbcs+xzXlp/Dr+PX4Tcz+/UydGKE5/++L+PdCInZcSER2vhbfH7hmsXOXVf/j61Ic/HgzECJDNV0jpLEt/8/guTXHMe6Ho5W+D8M3i5rqrllS1P20eOdlLNh6ERfjiz9kqzMYKxnU5BXokJqVZ7Tu5VJCBtrN+weLd1yu8v1+tfuqpEBV747MIzJKEgQB8WnSDOTaI7HILBHERZUYYltVhsFXYnoOkgyyGDfvZtfKUXZJ6cZfK0/e3wgD2vgAAG6kZKPje/9izPeHqz3Dklgi87PtXIL47yeCGuHA6/0ltzf3LizWH9jW1+j5rllhIJSek4/P/y3ObMen1s4guiL0A2T0bqRYbpLDO/eMv/+cfHcgnrq/sfi7n2vNDSopiYFQLVTTGSFzixcjEzJwt5Ifqsn3it/QqztDEZeajYVbL0q2fbf/mqQN+mDs7M00xBWl+S01FLnkB/r5uDQs3xeFNYdiSu279sh15BXosHR35WqnDDNAy/dGIdzIUHC5h6aWlJadj5z8wsf6x+e6l7p91qBWACzbNbb1bLykePdWajZuGHQHH79+F0O/3I+Pt0ViyJL9uBCXXipYq2mnYu9i3Kojkm2Bng74d0YfPNa1EbydNXBQq8TbDl9Lwe7IpGptU7yJzFlDN3u8MbQNlEoFJvQIBABM7tdcrHfs29pH3HfmwFZo28AFABBnhUHE+FVHkVtQ/F5xM9X66svuZubh/07HiTOUG/4dAbDY6D6tTkC6iWJ+V3tb9G7pBYUC6BLgBrWNfOEIA6HaqIYzQqYyBs/3alpq2+kSo0XMZZgF0nfn6HRCtXwLX7D1ojjpnCm9P92DRTsu45FlB9Hj491oN3c7wn49bZH7LznK7rFvDuHbfcbbU9XC7Uwzjr92+55k7hm56T/8PBzV6NvKG238nMXbPhvVCQ+2KcweHI1JwdmbaVW+v7TsfExeKy3ufWJFeKmg4XLiPazYF4XIhAwMW3oAD36+D7F35PuQm/vneSQbfJt20tjgw0c7opVv4eOlUChKDZ1/c+NZRCdbNpOmp9MJpTJCevMfbQ+voq6Nd4e3w78z+mD24Nbi7YGeDnikiz86NHTBuB6B6Nm8sKssrhYsr1ARyfdycSo2VbLtVhn1lTVl48mb6PXJbhwxc06sT/+JxKu/nBJrNu1LBEJZFhpQkpGTL36h79Fc2j2qUirQLdAD/ze1F1ZNuN8i91dZDIRqIX1GqKYy9enZpSN2tY0S7w5vh+CmHpLtE1Yfw9ZK1NgYBkKp2XnILdDiwS/24vHlhyweDP1t5iy2S3ddER9jnQD8dTrOIlkhczNeWp0Atar4JViRbo17uQUYuGgfPtxysdx939x4FoMW75N1GoTNp25hyrqTSMrIETMtDYpS4e2KsgMA8EAzTzT3cYRj0Rvzsz8cQV5B1Z4TUx/e5cnO1+L49ZTyd7Sw9/46j+AFO3H2VnEQOCqoEc69Pxi9WnpJ9m3s6SD+29nOBsn3cvGVhQrxDQmCgO8OXDM5DLqZV3FAplIq0MrXWTJLvkKhwJdPdcXfr/aGq70t/ItqiUxlmGorfd1aYw8HPBHUCEDx39eCrRfx+u+nZeliff33M7h5NxuLd5rXxf7L0cIi79X/xeBQVDLu5UgDH0tNvaE/j6NahXWTHsDQDn4AADuDcowODV3h4ai2yP1VFgOhWqg4K1n9L6izN9MkaV4A8HJS4/MnOgMAvJ1LF7B9uy+q1LbyGH6zTc3Kx8272Yi5k4XTN1Jx827tSS1XZQFAfRCVZiSwNCYzrwA2quIPixOxxSPcUrPyMO2XU/i/03EACmusZv92Gsv3RmHw4v0Y98MRXKlAMXF6TgFW7LsmyxBVQRDw2oYIbDkTj7WHY8Wh6o3cCz8M+xXVutgoFWjkbg+NjQoLR3YCUPhYbjsXj1+P3TA7g5hboMWKfVG4nFjYtVZeTVqDMmoT5PigXnMoBoklaoP8TUz2OKl3M3g5aTC5X3MsfborgMK/I0vbfyUZH2+LBFD4/hAxdyDa+xcHsPrn0lwNi/avDQtuluduZp7YvayvW2vu7Sj+3cSnFQb33+2/hl+P30TEjeqbGkAQBGw/F4+rBq/95Hu5KCh6XVxOrHjN1TMrj5TqFUjNyrPIaFP9ly83h8JAZ9bg1ghu6oFlz9xX5XNbEmeWroX036SW77Vc5b4xG47F4o0/zoq/dwlww8xBrdC7pbe4zdhrIT2nYmnTpPQcyRteWna+JGsSl5pjsRmQM3LyoVIqKv2t7FZqNvxc7TDpfyfg5aTGF090lnyzNWXqupM4eDUZG1/pYXZGKCtXK8l2PLEiHM/1bIq5I9rhu/3X8NfpOPx1Og4jOvtj36Xb+O3EzUpdk96KfVFYsS8KJ94JrdGhqoZDq8/HpSPmTuGHSYB7YTZjRKcGyMnTolugu/hYP9zZH78ciUX4tTuYvj4CAPD6H2cQ6OmATZN7wr2Mb5DL9kRh6a4rWHckFvtf7y8JhCb0CJTUank6qvHvjD6YseE0dl5MLHWump6HyVQRfwsf47NHdwlww7G3B0ChUOBOUQ3c9TtZuJdbACdN1d/eBUHA2VtpOGvQJS4IhR9sj3VtiPNx6eja2A02qop9p9aPLou5kwmtTjB7iHVNu5dbgAe/2AutTsD21/ogOrkw0Gjq5QTfokAoMT1HMvL23K10BDXxMHq+qvrrdBymr49AE08H7J3VD5EJGRj65QHx9nyD9xNBEMx67zLUwNUO8Wk5KNAJyMgtgIudLVKz8uBqb1vhc12MT8dj3xwCALg52AIAmns7YcNLIRU6T01gRqgW0r8n6OfiqC6GE7G52tti85SekiAIAEKK+nVb+Tph18y+ACAZbVOeu5l56L5gl2Tb2sPXccKgy+FWarZFvn2s2BeFju/9W6XUdHxqDg5fS8H+y7ex8eQtcbmKsuRrdfj7TDxSs/Lx1+k4hF8zb16Xv8/ElfpGvOq/wonoSo60sOQkkOZckyUZfkvdeTFR/Dbbyq+41uXJ+wPQrES9S2uD2iG9mDtZ+K+ceXP2FNX+xKZkQRAEsUj+4c7+kixGE08H/Phcdzjb2eL+QHej56rpEXfGRom93Lc5Brc3PuoKKP7i5OmkEet0LDX1wMoD1/Dw1//h83+Lu1xGdSvsEnr2gSaY/2gHfF2Jb/f6DFdqVj76fLpH9tXHTTkWnYK7WflIzynAU98dxt5LhRnVJp4O8CtaOichPQexKcV1WZUdCbfzQmK5NXGHrhbWAF2/k4Xrd7Kw84I0eM/ILcB7f53H17uvoMO8f/DX6TjsuZQkeU80NdzfQa1C+JwBsLct7Jae9+d5fLnzCrp8sAOrKjFpp+G6c/pAqLayukBo2bJlCAwMhJ2dHYKDg3H0aNlDun/77Te0adMGdnZ26NixI7Zu3VpDLa28ikbelWX4gnC2M/7t8YlujTD/kfb4dmw38YWfmac1+4NZ3z1hKC4tBwu2Roq/z/rtNCb973hFmm6UPnVfnrK+fc787TTGryr+mzJn5JLhrMRLdl7Bz4cL5yt6vGtDvNSnmcnjPtxyEf9dNV7cqDR4ZQqCZddoqmzNTGXFmCjeHdzOr8zj9EOvS4pKKj6fsQBaYzD6JDUrX8zqeDtr4GOw7lvYwFbo0NAVAEyuB5eSaZmMUFJ6Dvp+tkcyhcLha3ew3WD4OYBSc3U19XLEm0PbQGMjLWY1Rf+YRVlgWPrVpHv44l9pzUn3QA/MHFhYBG1nq8LYB5pUao02d4MPxlup2Vh7JBYRN1Kr1N7qYLjcUGxKltgd3djDQRzufe5WOhYZTH9xrZxi9czcglKZxpjkTLzwv+MY8fVBMSh8d/M5jPn+sGRAhOEXp3NxabhuZEb0NYdi8Pm/l5GZp8W0X05h4upjWF30BUurE8RgriT3ou4rfb3OplO3xJqjL82sPTKUalAiUFDOEhtys6pAaMOGDQgLC8O8efNw8uRJdO7cGYMHD0ZSkvEho4cOHcLTTz+N559/HqdOncKjjz6KRx99FOfO1e5V1U195llqjpD0nHws23NV8m3X1BB6jY0KY0MC0dTLEY4aGzgXpdvN/TA1d+j2zotJ1T4x2cj7GmHb9N74d0YfDGpX/A37p+e74+W+zY0ec/1O+SNwYkzs4+WsweR+LSrczu/2R0nmPMot0FW5YNhQTa5NJQgCNhrJbP47ow9cy/mWaGoxUX2wsO1sPNrP+wd/n4mT3G74BhyfloPbRRkhLyeNZAK3Nn7F2SEfg+0PdWqAqf0LnzdTc6BU1Jaz8bh+JwtrDsVApytcb2v8qqN4+ecTYrC9JzIJo1aES44zVqNXFn0X2tUqZoTO3EzFwMX7StUPjuvRxCLDnBUKBZ59oHgOmfl/F86CXtvmFToSbfyLSmNPBwR4FBeqG9Z0XStj/itBEDD6u3D0/nS3ZHRfZELxWnfXbmciPi0bPx2+jv+u3pH8fRseczQ6RRzV+OVTXSTZzpLWHSn8crZiX5TJL52u9oWvR2PZm4qUQ0TdvoePt0VKll0JamI841pbWFUgtGjRIkyaNAkTJ05Eu3btsGLFCjg4OGDVqlVG9//yyy8xZMgQzJ49G23btsX8+fNx33334euvv67hlleMAsYjoQILBULf77+Gz/65JKllyTAzw6N/Yy6rduLwtTu4lJCBlMw8JBsEQs8ENzZ5DFA8A6lWJ1QoVb7rYiKmrDtp9LaQZoVde28ObYMvnuyMtg1c0NzbCc0Nai7a+7uiicHoG0PldSNtPRuP59YYf2Pxc7GDi72NJAO1clw3dGho+g0LABZsjZR8Ox793WGsLGc6gJK+fKqL+GFeUkYFa7wAIOJGqlkzi9+8m4XwqOIPj9+O38RpI+n+Zl7l14S1Mugae6hjA/ENWr80xppDMcjK02LqulOS1L9hkP7qLydx4EphV5q3swYtvJ0woI0PHuvaEC0N/gYMM0LdAz0wpGh0S2RCBhZuvYgbKVnYfi7eIjOA376Xi7i0bDHI0H8ITlxzrNRxXk4VG02jDx6rmhHacibeaH2gn4nMWWV88HAHPN29+D0hXyvgUJR5w79rwu2MXJwp+tvdMaOP+GXJy0mNpp6OcLGzFUdBGbqVmi1+sSzQ6jD2hyOY9L/j0OkE3LybjXO30pGTr5NM3xBtUGMUm5IpCWTf+OMsvtp1BXkFOsQZzG31v/DrOBpTmLFq7OEA2zLqtK4lZ+KzfyIl3VUl6Z9uRxO1Zeb+7b+96SxW7IsSawOHdvDDFBPvRbWF1RRL5+Xl4cSJE5gzZ464TalUIjQ0FOHh4UaPCQ8PR1hYmGTb4MGDsXnz5upsapWZzAhZaF6h74xMn96/tbeRPUvzdtbgWnKmZJ0mQ6dvpOKp7w6Lv9sWjYoa3S0Awzo0EL+ZGHM3Mx9ejhoM+/IAsvO12Da9t8kXJQDcuZeLa8mZeOOPs5LJEoHCPvyfngtGgIc9svO1Yr+3nr/BSCEPR7XJWU1/ORqLQE8HvFT0JrjnUhK0WgGhRRmlr8qYCHFwBz8oFApobJTixIe9W3phzaFok8cYc7oSXQbNvZ1MpqPL69bMyMnHr8dv4uHO/vB21uBQVDKeWXkEbfycsf21PmUe++ovp3AqNhXfjQ3CoPZ+RguQ/VzszCqu9XLS4J2H2uJCfDo+fLQDDlxJxks/nRCDcMNMXMf3/sG6SQ+gta+zJNAznJ3ay0kNG5USPxiZs8Rw5Fhzbyd4GgQg3+6/Js5L9fkTnTGqaNh0RRhORHo7I1fS3bH+6A2xW6IkdQWLkC2VESq5VIueJWf/VSoV6NrYDb8cLX5PKCubUtN+P3ETWp2AzgFuaOnrjJmDWsHF3gb3B3qIXdXfjLkPH/x9odTCt2uPXMcLvZvh6u17YiAecydTMnp2/t8X0NLHCX1aeYtF2EDh4BF9dkZvxb4oPNSpgclJdpt4OuLRLv5ldi/qlxoyRf+5Y2pustSs/DIHKeiVXGj6o8c6lvk+XhtYTUYoOTkZWq0Wvr7SokFfX18kJCQYPSYhIaFC+wNAbm4u0tPTJT81zVSNkKXmp3Cxk77IRgU1wtvD2pl1bHkZoZKrSuvnHfFwUpdbMHc3Kw/pOfm4lJiB2JQscWFOU97ZfA5PrAgvFQQBhWnexp4OUCgUcFDblHpMh3fyx/2B7vjgkfYAgEZl1Dks3BaJ7DwtEtJy8MKPx/HC/46LdS/GvnkrFUDk/CFi7YRhRsjOVgVHdfW/Kbja25p88yk5Z0hJX+68gvl/X8D8vy8AKEzBA4XZkbJmIRcEQZxsbkdREad++Pn347rhiyc6w95WhTeGtjZ1ilJe6N0Mi57sAge1DXxK/O0ZfgPOytNi5f5rZWatyupmclDb4KPHOuClPs3wQDMPk/Oa/BlR8QWIdTpB8sUhJTNP7K4DgPBrdzBulfFax6DAio0+0gdC1+9kVak71dT6UD7Oll0GQT85pN7F+HSzZ7q3tNg7WZL3L/2CqmOKsla2KiUm92uB+w2eE4VCgSeCAmCjVKCFjxOGd2oAAOIs5oYDBa7dziz19zll3UkIgiAJAONSs0t90czM04qzxrdt4ILfXpaOvHJ3sMX4HoH4v6m9KnfxAPq2KvwybOp92tSXX0MlM/l2tkpJPVhtVbvDNBksXLgQ77//vqxtMFUrrbVQRqhkH79+ziBz6N8IkzJykJ2nxYs/HYcgAN+P7wY7W5XJmiBPRzWaeTvCx1lj8gV1L6dAMpFXeRMAbjtnOqAtL3vm7qjGby/3EH9v7u0kLg0Ql5qNAA8H/HCwOHOz/Xw8opIyxWD0zK00NPF0EDsxOzdyFbuA7GxVsDPIQJVMWdfEtyM3B1uTBfCmusYEQcC5W+ni8PK/TsfhxT7NEGtQkHkrNRstfJxwMvYu4lNz0L+NN1RKBTQ2Kkn3qr7bR5w80c0Ooe188VjXhpUu/G7gWhhYJqbnIF+rKzXCKj0nHytMzOANlF9vMya4ifhvU89QRYOLnHwthi09IPmgu5uVV+4UC2MfaAJHjY04aZ+5GrjawdNRjTuZeXh8+X9YNeH+SgUvpubCsvQyCC1LTAsQfu0Oen+6G/tm9y/zdXL9TiZOxabioU4NxNfXtrPx+PX4Dcwc1BodGrqiQKvDrsgk9GrhVe5rLi0rH48sO4i7WflY9sx96NTIFdeSM6FSKjC0Y9lF/e38XbBnVj+42NuioGgE6fU7WcjMLZDUycSlZZfK0mbkFGD0d4dx3GCtr7i0bPG9/qW+zfDHiVtIvpeLtzcV1rb2auEpCcaA4i/PHRu5oo2fs2TxYlPGhzTBnGFtcSXxHradi8cr/Qqz3i/1aY5j0XeRV6JmMykjx+hITkMlR8HaqpQ1NvinKqwmI+Tl5QWVSoXERGmqPTExEX5+xv9Q/fz8KrQ/AMyZMwdpaWniz40bN6re+ApSmsoIWajy3twJ/4wRM0Lpudh2Lh4HriTj4NVksb872UQg5O6ghoPaBntm9TM5FDgjN18S/NwukelZtucqpqw7aVb9kKnH0OT+SgXee7g93nu4Pb4b1w3vDm8nCRBnbDiNr/cUd4NdT85EWna++IFvODeGTYkP+pJtKbmuT3Vw0tiYnEdG3zV28Eoypv1ySgxeN526hRFfH5TUog3/6iA2GizkejsjF3kFOoz9/gimrDuJdnP/waPLDkGnEySLT97NykO+Vid2BehrS6oy+s3HWQNblQIFOgEt395W6o36wJVksZvFSWMjdsvqeZjofqqIkvdZnksJGaW6e1Kz8stdkXvagJZ4c2gbSUBtDoVCgbCitdrO3UrH1HWnKjU0Xd++leO6YVxIYYDYppwPwcpw1Njgk5EdMaKzv7gt+V6eJDAoSacT8OwPR/Dahgh8XdQ1nZ2nxezfz2DPpdt4e1Ph3GhrDsXgpZ9O4KWfTpTbjq92XxG/hL29+SxWFpUPdG7kCme78jMaAR4OcLW3haeTRsyAXL+ThfNxxfVx8Wk5RrPXR6OlXUlxqTlikO/jbCepX7RVKTD1wZYAgPdGFGbx9Wvz6X08shOeCGqEt4a1MVn72KuFF17u1xx2tip0bOSK14e0Ea+zTytvRM4fUqqcwJxlekouQFvRv1+5WE0gpFarERQUhF27iuek0el02LVrF0JCjE/QFBISItkfAHbs2GFyfwDQaDRwcXGR/NQ0U58VlsgI5Wt1koxARUeliN0T93Il9R/6OpbkomxPyfSqR1HNhaPGBg3dil+chkOkk9JzJS+29OzCdubka3HnXi4+++cStpyJxzYzilafNfh2X1mjghqZDNpu3M0SR4q4OdhKXvAtS6T78wpqLtXfvakHPnikPRQKhSQj1LeVtzji73pKJm6lZuPZH47gr9NxWFYU4JWVYdNLvpeL5Hu5koVlL8an4+bdbMkipQlpOeKHqVJRPLNsVSiVilLBnaej2mhR+P2B7tg9s5/kzbyik/75GHltmDOKLD4tW8wcpRgJeOb9dR5LdhpfBqNzI1ccenNAhV+XhsYEN8Efr/SAxkaJo9EpePybQxXOZOkzVp5OhY/v090bY9GTXSrdprKMvr8xvhwtPXdsiZGYWXkF4t/XmVtp4uro/wuPQW6BFtvOxYsB/umbadh3+TZ+L5qA9ODVZLzw4zFxVOqhqGS8+sspXE0qzJrsu3wb3xtkf1Oz8vG/8MKurV4tpEuamEM/R9KzPxzBHoOh6glpOeLfz6xBrdAlwE1ynL5rKuJGKrYULWPk66JBE4PRaUFN3MX6oXEhgdgZ1qdUIXKXADd89kRnvNinOfbN7l+qfWobJX5+IVjMsBqjVCrw7vB2UNsoxTq12b+fkQyEMCauxELF+ix7bWc1gRAAhIWFYeXKlfjxxx9x8eJFvPLKK8jMzMTEiRMBAOPGjZMUU0+fPh3bt2/HF198gcjISLz33ns4fvw4pk6dKtclmMVUNsMSw+dLpuTnDG1ToeP1o2sOXEnG1rPFH5z6+TX0Q/IXPtYRrQ0CAk+DmgsPx+IgaeW4bmL6f+G2SLxisDimftXiF386gaAPd4rbv913DU+vLC7INtTGzxmrJ9yPJ7pVvKDVGMOgzdCvx2+KS43osx3TB7SEWqXEtAEtJfv2L1pCIrDo21l1ThPw/fhuGBcSCKDwW2pTL0cEeNhj+bP34fMnCzNc526lo++ne8Rj9HVO5owKuVMUCJUUcydTshxFfFqOmGlyc1BbbObgkgFVUy9HPNU9oNR+00NbIcDDAZN6Fy4c3M/MwQCGPh7ZsdQIP2PXbujE9bsIWbgbL/5UOJLwbgUnZHR3tMxjFdTEHV8UPd8X4tNxLKZia6bp3yfc7G3h42KHhY93RLsyhmdXlVKpwJD2xZl6/Yijz/6JRN/P9qDDvH/Q65M92H4uAa/8XJzhuZuVj21nE8TARW/8qqOS7qGdF5PwZ0QcBEHAtF8i8H+n48Supu+Lsj9jH2iCTZN7SDJfD7Y1PZGlKfrC+5JlAjF3MsWuIw9HjWQEXmtfZyx6snSJgo+znWSYfnt/V/HfSqUCLXycy+16KtW9aubHyDPBjXHxgyF49oHiL5Xf7TddcJ2WlY93ih7Thzo1wIpng/BiGfOo1SZWVSM0evRo3L59G3PnzkVCQgK6dOmC7du3iwXRsbGxUBrMRNejRw+sW7cO77zzDt566y20bNkSmzdvRocOHeS6BPOY+Lu2xPB5fdeTq70twuc8CIcKFu52DnCFjVJRqi36FZj1HxRezhr4uGhwqWhCRcNRMYY1Mz4udiZTz1eT7uHm3axS62OV1f/t6aQWAw9LaOxh+lvTxlOFXUb64HDGwFaY3L95qcnv3hvRHi28nfDYfQ0BVDwzYY5Wvk4IcHeQFMLbqpTY/lpv6HSFq0u7GYxEMXz+9l66jZd/OiE+V2VJvpdnNCuSkJ4jCYTu5RZgyJLCqf8tuaDiwsc7YsLqo8jJLwwmfV3sxNmU9VaO6yZ+236xb3O4Oqgl80aZ68E2vujf2gdN5xRPwpqVp0VWXoHJ141+Nvi9l27juTXHyuwCM8YS3Xd6wzv54+/T8dh+PgEX49PR08zsRr5WJ2ZXTI1mqw7fjLkPi3dexle7r2Lf5dvo3tRTOtJJEPCyQRDk5mCL1Kx8vLYhAkBhpuOPl3tgxq8R4qi5hm72UCiAm3ezcfjaHXRu5Cq+Rx2JTkFkQjqOFI1yGt8jEC18nLBxcg/M/fM8GrrZo3Oj4sDDXCVH1vVs4Yn/rt6RrFrv6aSWzGkV4GFv9HXi46yRdG+VNVeQKW8/1BbBzTwx67fTACo2+lilVCC0nY84431Z8wl9+k+k2HXcqaGrOAWFNbCqjBAATJ06FdevX0dubi6OHDmC4OBg8ba9e/dizZo1kv2feOIJXLp0Cbm5uTh37hyGDRtWwy2uOJM1QhYIhPT94B6O6goHQQDgbGeLN4e2QVMvR/i52GH24MIRQPplMvRdY56Oakl3kWGq3/BbjZPGBk4minp3XEhEH4PMhTkqWhtUnqYGk/pN7BmI1RNLD702fEMzNgOwu6Marw5oiUZFa2s936tpqYU+t03vjW6VmHSsjZ8zvhsbhH9n9DU6LFxjo4J9UU1SWd1T288niN0NZfl6z1V8seNSqe37L9/GUhOrnlvyw/2BZp648P4Q8Xd7dWFhuovB35DhB4eTxgbP92oq+VZdEQqFolRB7x8nb5nMzhoOd98dmSR++NnbqozOMh42sBVe6lu83dKrcOuLW69UYDFOfTZIoQBc7GtuxI9SqUBg0ZqD526lS2Z4N+b3l0MQYPBF5YmgRujYyBWbp/QUu4+efaAJ5g4vrKX5KyIOAxfvl5xjyJIDyNPq0NDNXuymd1Db4PMnOmPGwFaVKvQ17HJq5u2Ijx/vVGofLye1ZJbzRu6FI1xfC5Vmk/1c7dDY4G+3ZHeaOdwc1JIpHzpVMLjr0dwLHz1WmDwoa/64/64WL33zoAW/jNYEqwuE6oPqnEeouLui8m9wL/Ruhj2z+uHwWwPwXM/Crod7uQVITC+uHfEqKmzVMwyKerbwxEePdcDGyYWjtlxMBEIAUNHYrzJvFGVp6+cMlVIBW5UCz/Vsir4tS3exVHSSuebeTgifMwDNDOqj2jZwqVQdzScjO2FQe/O+eZn7nIcNbIXHuzaUbDMM9s7dKj2twd9n4k2ez93Rsh+mSqUCY4Ibw9nORvz7MwzqG1cy6DHl8yc6Y9qAlmIR7Lubz2HTKePD6NNNzNr9Yp9mCG5Weih8t0B3PBFU3LXnUcEJFMsT6FX4WFxPySxnz2KGWeOaXgy1gVv5r6WGbva49OEQtPBxxopng+Dvaod2DVwQNrCwaNhJY4PfXg7B/Efa44XeTcUuPcNC95K1fz1beFpsdJPhl5ymno6SWcv1fJztJO8b+kD9tdBW+Gtqz8JMTFsf2Nmq0DnADUPa++Hp7gFoasZEpKYsfLwjGrjaYc6wthU+tneLwve9pIwco13ogiCIQdLvL4eUqpOs7RgI1ULVmRHSv8lZ6lu6vVolzn57umiFarWNEs4aG/FbyJASH9QKhQJjgpvgvsaFGRBTw7wrQqVUYGgHP0ws+mC0FB8XO2x48QGsfeEBBHg4QKlU4Ifx3STdMQ3dK77WElCYGTL8v5Om4iMsKjIUv+Qkbaa82KeZ5LyR84fg/YdLdye/2KeZOHLFkH5uJj1LZzmAwknaTrwz0OiHnKVHqnQOcEPYwFbwM/imvyuy9ESRgOmaIHcHW6PdTD7OGklXSsk5vqqqsUfhB6c52T49/ZelmuwW0zMniO3Q0EXMvLb3d8WhOQOwdXpveBq8Jlv5OmNsSCBsVUo0dLOXZOqmPdgCK54Nwq8vhaBLgBs0NkrJ1AlVZfh8Bno5QmOjQu+W0m7JRu72ktnMDQuiOzVyw8l3B2L5s0EACru4V4wNwsLHO1UpWHu6e2OEzxlQaui9OfTBXE6+Dhm5BdDpBHy8LVIsSE/Lzhe/BOvX77MmVlUjVN9ZIhDSj2KxxCgevYZu9ki+l4cTRUNevZ00UCgUeLCNL3aG9SlzdAIAOGmq/uY/9oEmeO/h9uXvWAndSrxxDGjri81TnNHrk8JuuxYluk7M9dT9jdElwE1c88qhEvMLOVYgeDInQHDW2MDOVoWeLbzw0+Hr6NjQFXa2KqOjmEqm94HCoPbx+xoh0NNRnCSwOgIhQDqnzRPdGuHbfdeMLnlgKXcMCqXzCgTsjkxEUGMPnI9PQ1ATd2hsVEgxMT+Qu6O6VC0TUBioGF6HqfXVKksfWBQu6aE1uXjrkWt3sGDrRcwY2AqnbxQO+ZZjIjxj7xVPdmuEVx9siQVbL2JXZFKFl2tQKBQY2M4XW87Go42fM8IGFXbnd2/qgU2TeyA7X1upMgFTDK9Bn8H5cWJ3ZOYVYMOxG+gS4AaFQiEG8WqVEp1LZLLN/dJSU+xsVXC2s0FGTgFeXXcKbfycxdnWB7b1FYvAS5ZEWAsGQrWQyYxQFbvGYpIz8en2wvoOP9fKD88tqbm3E07fTBOnkvc3SG+38Ck/RVrVjJCfix3G1/AwzUbuDnh9SGvczcxDl4DKLSioUiok9VKOZswvZGerhFYniDN2V/QN/MNHO2Dun+dMdjmqirozB7XzxeoJ96Ntg8I36w4NXRDo6SBZe83TUSMZYu5sZ4Ot03rDSWMj6fazdFeVMdMHtEQTD0eEtq2+2oTWfs7iZKA7LyZKpo8IG9gK0wa0RLqJObrcHdRo5G6PcSFNkJWnRRMPB/i72YtB0IYXH0BCeg5CmntatM1eTmo4qFXIytPi1t1sNDMRaH209SLO3EzDhNXF653JkRFSKRWY2r8F1h2NFTNT3Zt6IsDDAcueuQ95Wl2lPmjnDGuDQC8HPNZVOoJKP/O8JRl2eemLm5VKBZztbPFCb2k92D+v9UF2vrZK0yXUFB9nDTJyCrDv8m3sMxjActhgYdrKZsflxkCoFlKa6LCsakbIcE2fklPbV4V+ccyL8YW1I37lZIBKMlUsbY5XH2yBmYPMX7LBkiqzsnxZylo0Ue/gGw9iy5l4zPvrPADzgidDzz7QBAPb+SJ4QeH8Wva2KigVENPauUWjsZRKhWT0ncZGhf97tRfm/30Bvx4vTId7Oqklsxb3aekt1jr4u9qjiacD4lKzS2XUqoOD2qbcRX2r6r2H2+PRZf8ZnZl70Y7LmDagJTJM1Ah5OKqhUCjwwSPGR6wGN7NsAKSnUCjQ2MMBkQkZuJ6SZTIQOmNkYVy7Gpj405hZg1tj5qBWmL4+ApcTM8TCW6VSATtl5drUyN0BswdXbKqQyrJXq/DZqE5Iy84vt2axvJmaaxNvZ41k7T49wwkrGzEQIksxtfp8VQMh/Vw/KqUCg80ssDVHySGdTU3MZmpKWcXS5anti/lVhOHz266BCxLSc0rNReJsZyMGsQ3d7Cs1FN+wi2bW4NZ4rmegOETcvowPP2c7WzFDBBROHGdYCOpiX/xcKJUKbJrcE/dyCtC4gn8PtVVzbyf877nueOybQyb30U8CWlJVBidUVUBRIHTDYKkUc5QsmK9JCoUCS5/uKtv9V9UT3UrPbWXtzFmqJcDdOl/rLJauhQQTM15VNRCKK+rHXTXhfov24/Zo7oXRRS98B7UKQzs2qNDxhjVCvVt64dNRnaAxc02j7k2rP9tQUwzn9tk6vTde7lt6uLXGRoWQ5p5Y8WwQVo7rVqn7KQyEfdHQzR7DOzWAQqHAG0PawFalwBflrDvX0GBx2kbu9pK/Iz8X6bdBD0d1nQmC9MrqLlq256pYtH3inVBJICFn14e+EDf2jvFAyNiSOwsf74gBlZhMkOqukrV+xkYUNqqBbvDqUHe+TtchOhMTD1d1+Ly+oM3ftWLDvcujUirwyahO+Hhkx0oVHno6qeHtrMG9nAIsHt0FXk4atPRxwjd7o8RVzMc+0ARhA1uh6/wdAID+rb3xSr8W4sizusCmxNpYTwQF4M+IOJyPKz1cvaqTlX07VhpEvdKvOSb2DCw3QO7TyhuPd22IVn7O4vP8Qq+mOBx9B88+UL1dU7VBWYHQZ/8U1t8pFYX7Ga6rZqpIuSbog9Ezt9Kw40Iilu66gneKJtkDiidDNaRf7oFIb1RQI/x8+DoGtPXB6PsD0LaBC1766QQuxKWjQCfARqlATwvXuNUUBkK1kKmAp+TKxea6lJCBfZeTxNqGkjOfWkplCw9tVUpsndYb2Xlasduma2N3jAtpIgZCvi4auBt8I7m/qUedygYBhRPC/e/QdXHJCHdHNbZM641DUcl4ZuURDKjmScrMyRLa2aqwqMS6UO8MLz2Evq4yp7DfUWMDpVKBUUGNsPHkTdln2O3T0hs2SgWORqeIC3x+8PcFbJnWGwBw827pTJFvBefGorqvQ0NXnHlvEDQ2KjEbtHlyT2gFAZcTM5BXoDNZg1bbMRCqhUx1gVVm1JhOJ+CZlYfFNcDsbVUmVyWXk7GuA8MhpK5F38RnD26NHRcSxa64uqSFjzMi5g0slT0IaeaJ7a/1rpHRV1Q2pRkTDOoLzh9o5okDbzwoWWdPDoFejhgV1Ajrj90Qt+mzjInpOXjRyOrsNT2RIlmHkl90lUoFlJCOfrVGrBGqhUxlhEx1mZUlIT1HDIKAwsyKpWZQrW6GgZB+mPaU/i2weUpPyeRpdYmxLhSFQoE2fi4WH+ZL1cNwcseGbva1Yl6VgUbWWku+lytmXIHC2ZXtbVV4tx5l+IgAZoRqJVM10ZXJCMUkS4c7lpwArzYzHKVQcm0uIrm08XNGZEIGWvg44dUHW+Cn8Os4XjSZKCBd66y2MDY8/3JiBm4a1ActfzYIzhobq/miRGQpzAjVQqYyQtpKpISi70gDIWt6i7NXq/DlU13wzkNt0dEKp22nuumjxzoiuKkHFj/ZBY90aViq5m5uLcyoOGlsSi3z8MzKIzhxvbBm6K1hbeBiZ8sgiOolZoRqIZM1QpXoGiuZESo5lXtt90gX+eYyITImqIk7NrwUIv5u+Hq9tmCYWXVEclj0ZBecuJ6CC/EZWLrrCgDgWExhJquhW+3LYhHVFAZCtZDOZCBU8a6x6OTCESHDOzWAn4tdhdfpIaKyGRYW19YgCCgckDCkQwP0aumNn8JjcNdgXTR/M1Z9J6qr2DVWC5mKd9YeuW4ySDIlpqhrbPT9AXhneLtat5gfkbWb2LMpAOC5ov/Xdk4aG3z19H2Sbda6RhSRJTAjZEUOXEnGqv+iJQv3lUWrE8TZZAM9HcvZm4gqI6iJO069O9CqvmQYLoqrVABejnVzFCaROZgRsjL6Fd7NEXMnE3laHextVZKlEYjIstwd1bW6W6wkwxXSdULt7tIjqm4MhGohU2uNAYCtyvw3rEsJGQAKV4fnGx0R6Rm+H1h6yR0ia8NAqBYqa7ogG6X5T9n5uDQAQFs/56o2iYjqmNdCWwIApj7YUuaWEMmLNUJWpuTCnGWJuJEKAOjUyK16GkNEVuu10FYYFxJYalVxovqGGSErY1OU0v77TBy6fbgT4VF3jO4nCALO3CjMCHUO4GSERFQagyAiBkK1UlkD5JVFM79OXXcKyfdy8dJPx43ul5adj4zcwtXmW/hY54rARERE1Y2BkLUp0TOmD3ZK0i+06qyxMbqQJxERETEQsj4l0kWmCqvv3CsMhDydmPomIiIyhYGQlck1c8GxlMxcAKwBICIiKgtHjdVCeQWmgx1TS2zka3X4dHsk7g/0QNTtTOQWaAEAHpwxloiIyCQGQrXQOw+1xezfzxi9rcBEIPRnRBxWHojGygPRku2ezAgRERGZxK6xWuiJbgE48Hp/o7cVmOgai0vNNrrdgzVCREREJjEQqqUCPByMbjeVEXJQGx8ZxowQERGRaQyErIzWRCCkMrGWGEeNERERmcZAyMoUaAVcTswotT0n33iXGYuliYiITGMgZGWOxqRg0OL9pbZn52uN7s+uMSIiItMYCFmBZt6OuK+xW5n75JoIhDiPEBERkWkMhKyAnY0Krw5oWeY+pjJCDISIiIhMYyBkBdQ2SnHVeVNyTARCdrZcZ4yIiMgUBkK12JtD28BBrcKHj3aAjbLspyq7qFh69uDWeKVfcwBAp0au1d5GIiIia8aZpWuxl/s2xwu9msJGpcSxmBST+2l1gpgRcndQY0xwY+gEAcM7+tdUU4mIiKwSA6FazkZVmAkyNU8QULjOmD4QsrNVws1BjTlD29ZI+4iIiKyZ1XSNpaSkYMyYMXBxcYGbmxuef/553Lt3r8z9X331VbRu3Rr29vZo3Lgxpk2bhrS0tBpsteWUVSOUJwmEWBNERERkLqsJhMaMGYPz589jx44d+Pvvv7F//368+OKLJvePi4tDXFwcPv/8c5w7dw5r1qzB9u3b8fzzz9dgqy2nrIxQXoFOnFDRnoEQERGR2ayia+zixYvYvn07jh07hm7dugEAvvrqKwwbNgyff/45/P1L18J06NABf/zxh/h78+bN8dFHH+HZZ59FQUEBbGys4tJFZRVL52t14vB5ja3VxLZERESys4pPzfDwcLi5uYlBEACEhoZCqVTiyJEjZp8nLS0NLi4uZQZBubm5SE9Pl/zUBmXWCBUUF0szI0RERGQ+qwiEEhIS4OPjI9lmY2MDDw8PJCQkmHWO5ORkzJ8/v8zuNABYuHAhXF1dxZ+AgIBKt9uSbFWsESIiIrI0WQOhN998EwqFosyfyMjIKt9Peno6HnroIbRr1w7vvfdemfvOmTMHaWlp4s+NGzeqfP+WUP6oscIaIQZCRERE5pO1UGbmzJmYMGFCmfs0a9YMfn5+SEpKkmwvKChASkoK/Pz8yjw+IyMDQ4YMgbOzMzZt2gRbW9sy99doNNBoat+K7eXVCLFrjIiIqOJkDYS8vb3h7e1d7n4hISFITU3FiRMnEBQUBADYvXs3dDodgoODTR6Xnp6OwYMHQ6PR4K+//oKdnZ3F2l7TysoIPfz1f+K/7VgsTUREZDar+NRs27YthgwZgkmTJuHo0aP477//MHXqVDz11FPiiLFbt26hTZs2OHr0KIDCIGjQoEHIzMzEDz/8gPT0dCQkJCAhIQFarfF1uWqz8tYa02PXGBERkfmsZgz52rVrMXXqVAwYMABKpRIjR47E0qVLxdvz8/Nx6dIlZGVlAQBOnjwpjihr0aKF5FzR0dEIDAyssbZbgqqMYmk9hQLQ2FhFbEtERFQrWE0g5OHhgXXr1pm8PTAwEIIgiL/369dP8ru1sy1n0VUAsLNRQaEwL3NEREREVtI1RmXXCOmxPoiIiKhi+MlpJcypEWJ9EBERUcUwELISSqUC5fV6MRAiIiKqGAZCVsQwK7R8zH0IbesruV1Xh2qiiIiIagIDIStiWCfUzt+l1AixvAJdTTeJiIjIqjEQsiKGI8dsVMpS64/la5kRIiIiqggGQlbEcC4hG6UCNirp01eXpgsgIiKqCQyErIjKoFpapVRAp5MGPloGQkRERBXCQMiKGI4as1UqkV8iEPJwVNdwi4iIiKwbAyErYpjwUakUyCsoXjOtta8zlj7VVYZWERERWS+rWWKDpGyUCmTlFQdC/8zoI2NriIiIrBMzQlbEsGvMRqlAtkEgRERERBXHQMiKSLrGlAo42zGhR0REVBX8JLVSCoUC80a0x92sCIx9oInczSEiIrJKDISsWKCXIzZP6Sl3M4iIiKwWu8asCGcJIiIisiwGQkRERFRvMRCyIorydyEiIqIKYCBkRdg1RkREZFkMhKwIF1UlIiKyLAZCVoRhEBERkWUxELIiJVebJyIioqphIGRF2DNGRERkWQyErAjjICIiIstiIGRFWCxNRERkWQyErAjDICIiIstiIGRFdMwIERERWRQDISvCOIiIiMiyGAhZEQZCRERElsVAyIrk63RyN4GIiKhOYSBkRdQqPl1ERESWxE9WK7JmYnc09nDAj891l7spREREdYJC4OQ0ZUpPT4erqyvS0tLg4uIid3OIiIjIDOZ+fjMjRERERPUWAyEiIiKqtxgIERERUb3FQIiIiIjqLQZCREREVG8xECIiIqJ6y2oCoZSUFIwZMwYuLi5wc3PD888/j3v37pl1rCAIGDp0KBQKBTZv3ly9DSUiIiKrYTWB0JgxY3D+/Hns2LEDf//9N/bv348XX3zRrGOXLFkChUJRzS0kIiIia2MjdwPMcfHiRWzfvh3Hjh1Dt27dAABfffUVhg0bhs8//xz+/v4mj42IiMAXX3yB48ePo0GDBjXVZCIiIrICVpERCg8Ph5ubmxgEAUBoaCiUSiWOHDli8risrCw888wzWLZsGfz8/My6r9zcXKSnp0t+iIiIqG6yikAoISEBPj4+km02Njbw8PBAQkKCyeNmzJiBHj164JFHHjH7vhYuXAhXV1fxJyAgoNLtJiIiotpN1kDozTffhEKhKPMnMjKyUuf+66+/sHv3bixZsqRCx82ZMwdpaWniz40bNyp1/0RERFT7yVojNHPmTEyYMKHMfZo1awY/Pz8kJSVJthcUFCAlJcVkl9fu3bsRFRUFNzc3yfaRI0eid+/e2Lt3r9HjNBoNNBqNuZdAREREVkzWQMjb2xve3t7l7hcSEoLU1FScOHECQUFBAAoDHZ1Oh+DgYKPHvPnmm3jhhRck2zp27IjFixdjxIgRVW88ERERWT2rGDXWtm1bDBkyBJMmTcKKFSuQn5+PqVOn4qmnnhJHjN26dQsDBgzA//73P3Tv3h1+fn5Gs0WNGzdG06ZNzb5vQRAAgEXTREREVkT/ua3/HDfFKgIhAFi7di2mTp2KAQMGQKlUYuTIkVi6dKl4e35+Pi5duoSsrCyL3m9GRgYAsGiaiIjICmVkZMDV1dXk7QqhvFCpntPpdIiLi4Ozs7NFJ2VMT09HQEAAbty4ARcXF4udtzap69fI67N+df0a6/r1AXX/Gnl9lScIAjIyMuDv7w+l0vTYMKvJCMlFqVSiUaNG1XZ+FxeXOvnHbaiuXyOvz/rV9Wus69cH1P1r5PVVTlmZID2rmEeIiIiIqDowECIiIqJ6i4GQTDQaDebNm1en5yyq69fI67N+df0a6/r1AXX/Gnl91Y/F0kRERFRvMSNERERE9RYDISIiIqq3GAgRERFRvcVAiIiIiOotBkIyWbZsGQIDA2FnZ4fg4GAcPXpU7iaZZeHChbj//vvh7OwMHx8fPProo7h06ZJkn379+kGhUEh+Xn75Zck+sbGxeOihh+Dg4AAfHx/Mnj0bBQUFNXkpRr333nul2t6mTRvx9pycHEyZMgWenp5wcnLCyJEjkZiYKDlHbb02AAgMDCx1fQqFAlOmTAFgnc/d/v37MWLECPj7+0OhUGDz5s2S2wVBwNy5c9GgQQPY29sjNDQUV65ckeyTkpKCMWPGwMXFBW5ubnj++edx7949yT5nzpxB7969YWdnh4CAAHz66afVfWkAyr6+/Px8vPHGG+jYsSMcHR3h7++PcePGIS4uTnIOY8/7xx9/LNlHrusDyn8OJ0yYUKr9Q4YMkexjrc8hAKOvSYVCgc8++0zcpzY/h+Z8LljqvXPv3r247777oNFo0KJFC6xZs6bqFyBQjVu/fr2gVquFVatWCefPnxcmTZokuLm5CYmJiXI3rVyDBw8WVq9eLZw7d06IiIgQhg0bJjRu3Fi4d++euE/fvn2FSZMmCfHx8eJPWlqaeHtBQYHQoUMHITQ0VDh16pSwdetWwcvLS5gzZ44clyQxb948oX379pK23759W7z95ZdfFgICAoRdu3YJx48fFx544AGhR48e4u21+doEQRCSkpIk17Zjxw4BgLBnzx5BEKzzudu6davw9ttvCxs3bhQACJs2bZLc/vHHHwuurq7C5s2bhdOnTwsPP/yw0LRpUyE7O1vcZ8iQIULnzp2Fw4cPCwcOHBBatGghPP300+LtaWlpgq+vrzBmzBjh3Llzwi+//CLY29sL3377razXl5qaKoSGhgobNmwQIiMjhfDwcKF79+5CUFCQ5BxNmjQRPvjgA8nzavialfP6yrtGQRCE8ePHC0OGDJG0PyUlRbKPtT6HgiBIris+Pl5YtWqVoFAohKioKHGf2vwcmvO5YIn3zmvXrgkODg5CWFiYcOHCBeGrr74SVCqVsH379iq1n4GQDLp37y5MmTJF/F2r1Qr+/v7CwoULZWxV5SQlJQkAhH379onb+vbtK0yfPt3kMVu3bhWUSqWQkJAgblu+fLng4uIi5ObmVmdzyzVv3jyhc+fORm9LTU0VbG1thd9++03cdvHiRQGAEB4eLghC7b42Y6ZPny40b95c0Ol0giBY93MnCEKpDxmdTif4+fkJn332mbgtNTVV0Gg0wi+//CIIgiBcuHBBACAcO3ZM3Gfbtm2CQqEQbt26JQiCIHzzzTeCu7u75BrfeOMNoXXr1tV8RVLGPkRLOnr0qABAuH79uritSZMmwuLFi00eU1uuTxCMX+P48eOFRx55xOQxde05fOSRR4QHH3xQss2ansOSnwuWeu98/fXXhfbt20vua/To0cLgwYOr1F52jdWwvLw8nDhxAqGhoeI2pVKJ0NBQhIeHy9iyyklLSwMAeHh4SLavXbsWXl5e6NChA+bMmYOsrCzxtvDwcHTs2BG+vr7itsGDByM9PR3nz5+vmYaX4cqVK/D390ezZs0wZswYxMbGAgBOnDiB/Px8yXPXpk0bNG7cWHzuavu1GcrLy8PPP/+M5557TrKgsDU/dyVFR0cjISFB8py5uroiODhY8py5ubmhW7du4j6hoaFQKpU4cuSIuE+fPn2gVqvFfQYPHoxLly7h7t27NXQ15klLS4NCoYCbm5tk+8cffwxPT0907doVn332maTLwRqub+/evfDx8UHr1q3xyiuv4M6dO+Jtdek5TExMxJYtW/D888+Xus1ansOSnwuWeu8MDw+XnEO/T1U/O7noag1LTk6GVquVPNkA4Ovri8jISJlaVTk6nQ6vvfYaevbsiQ4dOojbn3nmGTRp0gT+/v44c+YM3njjDVy6dAkbN24EACQkJBi9fv1tcgoODsaaNWvQunVrxMfH4/3330fv3r1x7tw5JCQkQK1Wl/qA8fX1Fdtdm6+tpM2bNyM1NRUTJkwQt1nzc2eMvk3G2mz4nPn4+Ehut7GxgYeHh2Sfpk2bljqH/jZ3d/dqaX9F5eTk4I033sDTTz8tWcBy2rRpuO++++Dh4YFDhw5hzpw5iI+Px6JFiwDU/usbMmQIHn/8cTRt2hRRUVF46623MHToUISHh0OlUtWp5/DHH3+Es7MzHn/8ccl2a3kOjX0uWOq909Q+6enpyM7Ohr29faXazECIKm3KlCk4d+4cDh48KNn+4osviv/u2LEjGjRogAEDBiAqKgrNmzev6WZWyNChQ8V/d+rUCcHBwWjSpAl+/fXXSr/IaqsffvgBQ4cOhb+/v7jNmp+7+i4/Px9PPvkkBEHA8uXLJbeFhYWJ/+7UqRPUajVeeuklLFy40CqWbnjqqafEf3fs2BGdOnVC8+bNsXfvXgwYMEDGllneqlWrMGbMGNjZ2Um2W8tzaOpzoTZj11gN8/LygkqlKlUtn5iYCD8/P5laVXFTp07F33//jT179qBRo0Zl7hscHAwAuHr1KgDAz8/P6PXrb6tN3Nzc0KpVK1y9ehV+fn7Iy8tDamqqZB/D585aru369evYuXMnXnjhhTL3s+bnDihuU1mvNz8/PyQlJUluLygoQEpKitU8r/og6Pr169ixY4ckG2RMcHAwCgoKEBMTA6D2X19JzZo1g5eXl+Tv0tqfQwA4cOAALl26VO7rEqidz6GpzwVLvXea2sfFxaVKX1QZCNUwtVqNoKAg7Nq1S9ym0+mwa9cuhISEyNgy8wiCgKlTp2LTpk3YvXt3qVSsMREREQCABg0aAABCQkJw9uxZyRuX/s27Xbt21dLuyrp37x6ioqLQoEEDBAUFwdbWVvLcXbp0CbGxseJzZy3Xtnr1avj4+OChhx4qcz9rfu4AoGnTpvDz85M8Z+np6Thy5IjkOUtNTcWJEyfEfXbv3g2dTicGgiEhIdi/fz/y8/PFfXbs2IHWrVvL3qWiD4KuXLmCnTt3wtPTs9xjIiIioFQqxe6k2nx9xty8eRN37tyR/F1a83Oo98MPPyAoKAidO3cud9/a9ByW97lgqffOkJAQyTn0+1T5s7NKpdZUKevXrxc0Go2wZs0a4cKFC8KLL74ouLm5Sarla6tXXnlFcHV1Ffbu3SsZxpmVlSUIgiBcvXpV+OCDD4Tjx48L0dHRwp9//ik0a9ZM6NOnj3gO/TDJQYMGCREREcL27dsFb2/vWjHEfObMmcLevXuF6Oho4b///hNCQ0MFLy8vISkpSRCEwiGgjRs3Fnbv3i0cP35cCAkJEUJCQsTja/O16Wm1WqFx48bCG2+8Idlurc9dRkaGcOrUKeHUqVMCAGHRokXCqVOnxFFTH3/8seDm5ib8+eefwpkzZ4RHHnnE6PD5rl27CkeOHBEOHjwotGzZUjL0OjU1VfD19RXGjh0rnDt3Tli/fr3g4OBQI0OTy7q+vLw84eGHHxYaNWokRERESF6T+pE2hw4dEhYvXixEREQIUVFRws8//yx4e3sL48aNqxXXV941ZmRkCLNmzRLCw8OF6OhoYefOncJ9990ntGzZUsjJyRHPYa3PoV5aWprg4OAgLF++vNTxtf05LO9zQRAs896pHz4/e/Zs4eLFi8KyZcs4fN6affXVV0Ljxo0FtVotdO/eXTh8+LDcTTILAKM/q1evFgRBEGJjY4U+ffoIHh4egkajEVq0aCHMnj1bMheNIAhCTEyMMHToUMHe3l7w8vISZs6cKeTn58twRVKjR48WGjRoIKjVaqFhw4bC6NGjhatXr4q3Z2dnC5MnTxbc3d0FBwcH4bHHHhPi4+Ml56it16b3zz//CACES5cuSbZb63O3Z88eo3+T48ePFwShcAj9u+++K/j6+goajUYYMGBAqWu/c+eO8PTTTwtOTk6Ci4uLMHHiRCEjI0Oyz+nTp4VevXoJGo1GaNiwofDxxx/Lfn3R0dEmX5P6uaFOnDghBAcHC66uroKdnZ3Qtm1bYcGCBZIgQs7rK+8as7KyhEGDBgne3t6Cra2t0KRJE2HSpEmlvjha63Oo9+233wr29vZCampqqeNr+3NY3ueCIFjuvXPPnj1Cly5dBLVaLTRr1kxyH5WlKLoIIiIionqHNUJERERUbzEQIiIionqLgRARERHVWwyEiIiIqN5iIERERET1FgMhIiIiqrcYCBEREVG9xUCIiIiI6i0GQkRUq+3duxcKhaLUgo01ZdeuXWjbti20Wq247bvvvkNAQACUSiWWLFlS5vHbt29Hly5doNPpqrmlRFQZnFmaiGqNfv36oUuXLpLgIi8vDykpKfD19YVCoajxNgUFBSEsLAxjxowBULhoq5eXFxYtWoSRI0fC1dUVDg4OZZ7j/vvvx7Rp0zB27NiaaDIRVQAzQkRUq6nVavj5+ckSBB08eBBRUVEYOXKkuC02Nhb5+fl46KGH0KBBg3KDIACYMGECli5dWp1NJaJKYiBERLXChAkTsG/fPnz55ZdQKBRQKBSIiYkp1TW2Zs0auLm54e+//0br1q3h4OCAUaNGISsrCz/++CMCAwPh7u6OadOmSbqzcnNzMWvWLDRs2BCOjo4IDg7G3r17y2zT+vXrMXDgQNjZ2Yn33bFjRwBAs2bNxDaePn0a/fv3h7OzM1xcXBAUFITjx4+L5xkxYgSOHz+OqKgoyz5oRFRlNnI3gIgIAL788ktcvnwZHTp0wAcffAAA8Pb2RkxMTKl9s7KysHTpUqxfvx4ZGRl4/PHH8dhjj8HNzQ1bt27FtWvXMHLkSPTs2ROjR48GAEydOhUXLlzA+vXr4e/vj02bNmHIkCE4e/YsWrZsabRNBw4cwDPPPCP+Pnr0aAQEBCA0NBRHjx5FQEAAvL29MXz4cHTt2hXLly+HSqVCREQEbG1txeMaN24MX19fHDhwAM2bN7fgo0ZEVcVAiIhqBVdXV6jVajg4OMDPz6/MffPz87F8+XIxqBg1ahR++uknJCYmwsnJCe3atUP//v2xZ88ejB49GrGxsVi9ejViY2Ph7+8PAJg1axa2b9+O1atXY8GCBUbv5/r16+L+AGBvbw9PT08AhUGavp2xsbGYPXs22rRpAwBGAyt/f39cv369go8KEVU3BkJEZHUcHBwkmRVfX18EBgbCyclJsi0pKQkAcPbsWWi1WrRq1UpyntzcXDGwMSY7O1vsFitLWFgYXnjhBfz0008IDQ3FE088USrzY29vj6ysLLOuj4hqDgMhIrI6ht1OAKBQKIxu0w9Zv3fvHlQqFU6cOAGVSiXZzzB4KsnLywt3794ttz3vvfcennnmGWzZsgXbtm3DvHnzsH79ejz22GPiPikpKfD29i73XERUsxgIEVGtoVarJQXOltK1a1dotVokJSWhd+/eFTruwoULZu3bqlUrtGrVCjNmzMDTTz+N1atXi4FQTk4OoqKi0LVr10q1n4iqD0eNEVGtERgYiCNHjiAmJgbJyckWm4SwVatWGDNmDMaNG4eNGzciOjoaR48excKFC7FlyxaTxw0ePBgHDx4s89zZ2dmYOnUq9u7di+vXr+O///7DsWPH0LZtW3Gfw4cPQ6PRICQkxCLXQ0SWw0CIiGqNWbNmQaVSoV27dvD29kZsbKzFzr169WqMGzcOM2fOROvWrfHoo4/i2LFjaNy4scljxowZg/Pnz+PSpUsm91GpVLhz5w7GjRuHVq1a4cknn8TQoUPx/vvvi/v88ssvGDNmjFlzDhFRzeLM0kREZZg9ezbS09Px7bffVur45ORktG7dGsePH0fTpk0t3DoiqipmhIiIyvD222+jSZMmle6mi4mJwTfffMMgiKiWYkaIiIiI6i1mhIiIiKjeYiBERERE9RYDISIiIqq3GAgRERFRvcVAiIiIiOotBkJERERUbzEQIiIionqLgRARERHVWwyEiIiIqN76f3gjje/0HrVjAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "ax = plots[0].plot(title=\"AutoCorrelation\")\n", "ax;" ] }, { "cell_type": "code", "execution_count": 14, "id": "89144f53-07a4-44ef-9628-6a2962b37c31", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkQAAAHHCAYAAABeLEexAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/P9b71AAAACXBIWXMAAA9hAAAPYQGoP6dpAABrD0lEQVR4nO3deXhTZdoG8PtkbdrSdKELhRbKIqsIgkIFRLGIKIyM6KjDKIO4Dq4oKs4I7kX9xgUF3BAcR0RQ0VEHFMuiOCyyKYiWrdAKtIWWNl2znu+P5JyepEm6N2ly/66r0pycJG8Tmz553ud9XkEURRFEREREYUwV6AEQERERBRoDIiIiIgp7DIiIiIgo7DEgIiIiorDHgIiIiIjCHgMiIiIiCnsMiIiIiCjsMSAiIiKisMeAiIiIiMIeAyIiIiIKewyIiMLU8uXLIQgCdu7cKR974oknIAiC/KXVatGjRw/ce++9KCsra/R9f/HFFxg7diySkpIQGRmJnj174k9/+hPWrVvXBj9J0xw4cABPPPEEjh07FuihEFEQ0QR6AEQUfJYsWYLo6GhUVVUhJycHr732Gnbv3o0tW7Y0eNv/+7//w5w5czB27FjMnTsXkZGROHz4ML799lusXLkSV1xxRTv8BL4dOHAATz75JC655BL06NEjoGMhouDBgIiI6rn22mvRuXNnAMAdd9yBG264AR999BF27NiBCy+80OftbDYbnn76aYwfPx7ffPNNveuLi4vbbMxtQRRF1NbWwmAwBHooRNTGOGVGRA0aM2YMAODIkSN+zztz5gxMJhNGjRrl9fqkpCT5+02bNkEQBHz00Ud47LHHkJKSgqioKPzhD39AQUFBvdtu374dV1xxBYxGIyIjIzF27Fj88MMP9c47ceIEZs6cidTUVOj1emRkZOCuu+6CxWLB8uXLcd111wEALr30UnlqcNOmTQCAHj16YNKkSfj6668xfPhwGAwGvPnmmzh27BgEQcDy5cvrPZ4gCHjiiSfky9K048GDB/GXv/wFRqMRiYmJePzxxyGKIgoKCnD11VcjJiYGKSkp+Oc//+n3OSWi9sGAiIgaJNXbxMXF+T0vKSkJBoMBX3zxBUpLSxt1388++yy++uorPPLII7j33nuxfv16ZGVloaamRj5nw4YNuPjii2EymTB//nw899xzKCsrw7hx47Bjxw75vJMnT+LCCy/EypUrcf3112PhwoW46aabsHnzZlRXV+Piiy/GvffeCwB47LHH8P777+P9999H//795fvIzc3FjTfeiPHjx+PVV1/FkCFDGvksubv++uvhcDiwYMECjBgxAs888wxeeeUVjB8/Hl27dsXzzz+P3r1746GHHsJ3333XrMcgolYkElFYWrZsmQhA/PHHH+Vj8+fPFwGIubm54unTp8Vjx46J7777rmgwGMTExESxqqqqwfudN2+eCECMiooSJ06cKD777LPirl276p23ceNGEYDYtWtX0WQyycdXrVolAhBfffVVURRF0eFwiH369BEnTJggOhwO+bzq6moxIyNDHD9+vHzs5ptvFlUqldvPJJFuu3r1ahGAuHHjxnrndO/eXQQgrlu3zu14Xl6eCEBctmxZvdsAEOfPny9flp7D22+/XT5ms9nEbt26iYIgiAsWLJCPnz17VjQYDOL06dPr3S8RtS9miIionr59+yIxMRE9evTALbfcgt69e2Pt2rWIjIxs8LZPPvkkVqxYgaFDh+Lrr7/G3//+dwwbNgznn38+fv3113rn33zzzejUqZN8+dprr0WXLl3w3//+FwCwd+9eHDp0CH/+859RUlKCM2fO4MyZM6iqqsJll12G7777Dg6HAw6HA5999hkmT56M4cOH13scQRAa9bNnZGRgwoQJjTrXn1tvvVX+Xq1WY/jw4RBFETNnzpSPx8bGom/fvjh69GiLH4+IWoZF1URUzyeffIKYmBicPn0aCxcuRF5eXpMKi2+88UbceOONMJlM2L59O5YvX44VK1Zg8uTJ2L9/PyIiIuRz+/Tp43ZbQRDQu3dveZru0KFDAIDp06f7fLzy8nJYLBaYTCYMGjSoCT9pfRkZGS26vSQ9Pd3tstFoREREhFysrjxeUlLSKo9JRM3HgIiI6rn44ovlP9yTJ0/Gueeei2nTpmHXrl1QqRqfWI6JicH48eMxfvx4aLVavPfee9i+fTvGjh3b6PtwOBwAgBdffNFnPU90dHSja5Ya4i3w85VdstvtPu9HrVY36hjgXM1GRIHFgIiI/IqOjsb8+fMxY8YMrFq1CjfccEOz7mf48OF47733cOrUKbfjUgZIIooiDh8+jMGDBwMAevXqBcAZXGVlZfm8/8TERMTExGD//v1+x9HYqTMlqZjcsznl8ePHm3xfRBScWENERA2aNm0aunXrhueff97vedXV1di6davX69auXQvAWZ+k9K9//QsVFRXy5Y8//hinTp3CxIkTAQDDhg1Dr1698H//93+orKysd7+nT58GAKhUKkyZMgVffPGFW/dtiZSFiYqKAlA/uPEnJiYGnTt3rrcabPHixY2+DyIKbswQEVGDtFot7rvvPsyZMwfr1q3z2W26uroaF110EUaOHIkrrrgCaWlpKCsrw2effYbvv/8eU6ZMwdChQ91uEx8fj9GjR2PGjBkoKirCK6+8gt69e+O2224D4Ax03nnnHUycOBEDBw7EjBkz0LVrV5w4cQIbN25ETEwMvvjiCwDAc889h2+++QZjx47F7bffjv79++PUqVNYvXo1tmzZgtjYWAwZMgRqtRrPP/88ysvLodfrMW7cOLceSd7ceuutWLBgAW699VYMHz4c3333HQ4ePNgKzy4RBQMGRETUKLfffjueeeYZLFiwwGdAFBsbi7fffhtfffUVli1bhsLCQqjVavTt2xcvvvii3ANI6bHHHsPPP/+M7OxsVFRU4LLLLsPixYvdVrRdcskl2Lp1K55++mm8/vrrqKysREpKCkaMGIE77rhDPq9r167Yvn07Hn/8cXzwwQcwmUzo2rUrJk6cKN9fSkoK3njjDWRnZ2PmzJmw2+3YuHFjgwHRvHnzcPr0aXz88cdYtWoVJk6ciLVr1zZ4OyLqGASR1XxEFACbNm3CpZdeitWrV+Paa68N9HCIKMyxhoiIiIjCHgMiIiIiCnsMiIiIiCjsBTQgstvtePzxx5GRkQGDwYBevXrh6aefdmtSJooi5s2bhy5dusBgMCArK6te3xIi6nguueQSiKLI+iEiCgoBDYief/55LFmyBK+//jp+/fVXPP/883jhhRfw2muvyee88MILWLhwId544w1s374dUVFRmDBhAmprawM4ciIiIgolAV1lNmnSJCQnJ2Pp0qXysalTp8JgMODf//43RFFEamoqHnzwQTz00EMAnHsWJScnY/ny5c3umEtERESkFNA+RBdddBHeeustHDx4EOeccw5++uknbNmyBS+99BIAIC8vD4WFhW7t+o1GI0aMGIGtW7d6DYjMZjPMZrN82eFwoLS0FAkJCc1q2U9ERETtTxRFVFRUIDU1tUl7KDZXQAOiRx99FCaTCf369YNarYbdbsezzz6LadOmAQAKCwsBAMnJyW63S05Olq/zlJ2djSeffLJtB05ERETtoqCgAN26dWvzxwloQLRq1Sp88MEHWLFiBQYOHIi9e/fi/vvvR2pqKqZPn96s+5w7dy5mz54tXy4vL0d6ejoKCgoQExPTWkMnIiKiNmQymZCWloZOnTq1y+MFNCCaM2cOHn30UXnq69xzz8Xx48eRnZ2N6dOnIyUlBQBQVFSELl26yLcrKirCkCFDvN6nXq+HXq+vdzwmJoYBERERUQfTXuUuAV1lVl1dXW9eUK1Ww+FwAAAyMjKQkpKCnJwc+XqTyYTt27cjMzOzXcdKREREoSugGaLJkyfj2WefRXp6OgYOHIg9e/bgpZdewi233ALAGRXef//9eOaZZ9CnTx9kZGTg8ccfR2pqKqZMmRLIoRMREVEICWhA9Nprr+Hxxx/H3/72NxQXFyM1NRV33HEH5s2bJ5/z8MMPo6qqCrfffjvKysowevRorFu3DhEREQEcOREREYWSkN/t3mQywWg0ory8nDVEREREHUR7//3mXmZEREQU9hgQERERUdhjQERERERhjwERERERhT0GRERERBT2GBARERFR2GNARERERGGPARERERGFPQZEjVBjsQd6CERERNSGGBA14MMd+Rj0xNfI+bUo0EMhIiKiNsKAqAE/FZTB7hCx/4Qp0EMhIiKiNsKAqAF2h3OrN0dob/lGREQU1hgQNcDuCoRCfA9cIiKisMaAqAEOOUMU4IEQERFRm2FA1AC7KxDilBkREVHoYkDUAGaIiIiIQh8DogawqJqIiCj0MSBqgFRU7WCKiIiIKGQxIGoAp8yIiIhCHwOiBsgZIk6ZERERhSwGRA2QaojYh4iIiCh0MSBqgJ1TZkRERCGPAVEDuMqMiIgo9DEgaoBDZIaIiIgo1DEgagBriIiIiEIfA6IGSFt32JkiIiIiClkMiBrAPkREREShjwFRAzhlRkREFPoYEDXAwcaMREREIY8BUQPYh4iIiCj0MSBqALfuICIiCn0MiBrgkGuIAjwQIiIiajMMiBpgY6dqIiKikBfQgKhHjx4QBKHe16xZswAAtbW1mDVrFhISEhAdHY2pU6eiqKioXccoZYjYh4iIiCh0BTQg+vHHH3Hq1Cn5a/369QCA6667DgDwwAMP4IsvvsDq1auxefNmnDx5Etdcc027jtHOrTuIiIhCniaQD56YmOh2ecGCBejVqxfGjh2L8vJyLF26FCtWrMC4ceMAAMuWLUP//v2xbds2jBw5sl3GaHc4/2UfIiIiotAVNDVEFosF//73v3HLLbdAEATs2rULVqsVWVlZ8jn9+vVDeno6tm7d2m7jYh8iIiKi0BfQDJHSZ599hrKyMvz1r38FABQWFkKn0yE2NtbtvOTkZBQWFvq8H7PZDLPZLF82mUwtGhf7EBEREYW+oMkQLV26FBMnTkRqamqL7ic7OxtGo1H+SktLa9H9ObjKjIiIKOQFRUB0/PhxfPvtt7j11lvlYykpKbBYLCgrK3M7t6ioCCkpKT7va+7cuSgvL5e/CgoKWjQ2qaia8RAREVHoCoqAaNmyZUhKSsJVV10lHxs2bBi0Wi1ycnLkY7m5ucjPz0dmZqbP+9Lr9YiJiXH7agk7M0REREQhL+A1RA6HA8uWLcP06dOh0dQNx2g0YubMmZg9ezbi4+MRExODe+65B5mZme22wgxgUTUREVE4CHhA9O233yI/Px+33HJLvetefvllqFQqTJ06FWazGRMmTMDixYvbdXxyhsjRrg9LRERE7SjgAdHll1/us8dPREQEFi1ahEWLFrXzqJxEUZRXlzFDREREFLqCooYoWCmX2jMgIiIiCl0MiPywKebJ2IeIiIgodDEg8kNZN8StO4iIiEIXAyI/7IogiBkiIiKi0MWAyA+7QxkQMSIiIiIKVQyI/HA4mCEiIiIKBwyI/HCbMmNEREREFLIYEPnh4JQZERFRWGBA5Id7UTUDIiIiolDFgMgPZVE14yEiIqLQxYDID2UfImaIiIiIQhcDIj/Yh4iIiCg8MCDyg32IiIiIwgMDIj+UQRDjISIiotDFgMgPm50ZIiIionDAgMgPZRBkZxERERFRyGJA5IedW3cQERGFBQZEftjdaogYEREREYUqBkR+cOsOIiKi8MCAyA9OmREREYUHBkR+cC8zIiKi8MCAyA/l1h2Mh4iIiEIXAyI/mCEiIiIKDwyI/FAWVbMPERERUehiQOSHMghigoiIiCh0MSDyg1NmRERE4YEBkR/sQ0RERBQeGBD5YWMfIiIiorDAgMgPz6wQt+8gIiIKTQyI/PBcWcYsERERUWhiQORH/YCIEREREVEoYkDkh2cAxICIiIgoNDEg8sPucL/scHg/j4iIiDo2BkR+2D0yQuU1Vtz41jas2J4foBERERFRWwh4QHTixAn85S9/QUJCAgwGA84991zs3LlTvl4URcybNw9dunSBwWBAVlYWDh061C5jc3jUEO04VoqtR0vw4Q4GRERERKEkoAHR2bNnMWrUKGi1WqxduxYHDhzAP//5T8TFxcnnvPDCC1i4cCHeeOMNbN++HVFRUZgwYQJqa2vbfHyeRdVWm3POzOo5l0ZEREQdmiaQD/78888jLS0Ny5Ytk49lZGTI34uiiFdeeQX/+Mc/cPXVVwMA/vWvfyE5ORmfffYZbrjhhjYdn2cRtc1VRGTj+nsiIqKQEtAM0X/+8x8MHz4c1113HZKSkjB06FC8/fbb8vV5eXkoLCxEVlaWfMxoNGLEiBHYunWr1/s0m80wmUxuX83lmSGy2EWvx4mIiKhjC2hAdPToUSxZsgR9+vTB119/jbvuugv33nsv3nvvPQBAYWEhACA5OdntdsnJyfJ1nrKzs2E0GuWvtLS0Zo/Ps6jaZpcyRJwyIyIiCiUBDYgcDgfOP/98PPfccxg6dChuv/123HbbbXjjjTeafZ9z585FeXm5/FVQUNCC8XkGRK4MkZ0ZIiIiolAS0ICoS5cuGDBggNux/v37Iz/fuYorJSUFAFBUVOR2TlFRkXydJ71ej5iYGLev5vKsFbKyhoiIiCgkBTQgGjVqFHJzc92OHTx4EN27dwfgLLBOSUlBTk6OfL3JZML27duRmZnZ5uPzzBBZbc7LDIiIiIhCS0BXmT3wwAO46KKL8Nxzz+FPf/oTduzYgbfeegtvvfUWAEAQBNx///145pln0KdPH2RkZODxxx9HamoqpkyZ0ubj86whkpbb27jsnoiIKKQENCC64IILsGbNGsydOxdPPfUUMjIy8Morr2DatGnyOQ8//DCqqqpw++23o6ysDKNHj8a6desQERHR5uPzjHukKTOuMiMiIgotAQ2IAGDSpEmYNGmSz+sFQcBTTz2Fp556qh1H5VSvD5GdU2ZEREShKOBbdwQzz0yQNFXGDBEREVFoYUDkR72tOxx1GSJRZFBEREQUKhgQ+eE5ZSbtZQYwS0RERBRKGBD5UW/KTHGZdUREREShgwGRH/UyRHZmiIiIiEIRAyI/6tUQKQIiG7fvICIiChkMiPzw7EOkDIK4wSsREVHoYEDkh90j6LFwyoyIiCgkMSDyw3NWzD1DxICIiIgoVDAg8sNzc1flNBkzRERERKGDAZEf9YuqRcX3rCEiIiIKFQyI/PDc7Z4ZIiIiotDEgMgPzykzq401RERERKGIAZEfnhkiNmYkIiIKTQyI/Ki/uauiMSMDIiIiopDBgMgPz6073Jbds6iaiIgoZDAg8sPfKjNmiIiIiEIHAyI/PHfnYA0RERFRaGJA5Ee9Zfd21hARERGFIgZEftQvqmYNERERUShiQOSHZ0DEDBEREVFoYkDkh2dApLzIGiIiIqLQwYDID89l90rMEBEREYUOBkR++MsC2T2XoBEREVGHxYDID2mVmUqof52yJxERERF1bAyI/JA2d9Wo6z9NrCEiIiIKHQyI/JAyRDovARFriIiIiEIHAyI/pDIhrbr+nJmdfYiIiIhCBgMiP+x+psyYISIiIgodDIj8kKbMtF6qqhkQERERhQ4GRH6wqJqIiCg8MCDyQ8oQabzUENm47J6IiChkMCDyQ8oCaVXeMkQsqiYiIgoVAQ2InnjiCQiC4PbVr18/+fra2lrMmjULCQkJiI6OxtSpU1FUVNRu46srqq6fIbJyyoyIiChkBDxDNHDgQJw6dUr+2rJli3zdAw88gC+++AKrV6/G5s2bcfLkSVxzzTXtNjZ/q8xYQ0RERBQ6NAEfgEaDlJSUesfLy8uxdOlSrFixAuPGjQMALFu2DP3798e2bdswcuTINh+bw98qM9YQERERhYyAZ4gOHTqE1NRU9OzZE9OmTUN+fj4AYNeuXbBarcjKypLP7devH9LT07F161af92c2m2Eymdy+mkuuIfKaIWINERERUagIaEA0YsQILF++HOvWrcOSJUuQl5eHMWPGoKKiAoWFhdDpdIiNjXW7TXJyMgoLC33eZ3Z2NoxGo/yVlpbWrLGJoghpVszrKjNOmREREYWMgE6ZTZw4Uf5+8ODBGDFiBLp3745Vq1bBYDA06z7nzp2L2bNny5dNJlOzgiJlvOMtQ8QpMyIiotAR8CkzpdjYWJxzzjk4fPgwUlJSYLFYUFZW5nZOUVGR15ojiV6vR0xMjNtXcyiLpjXsVE1ERBTSgiogqqysxJEjR9ClSxcMGzYMWq0WOTk58vW5ubnIz89HZmZmm49FKqgGWENEREQU6gI6ZfbQQw9h8uTJ6N69O06ePIn58+dDrVbjxhtvhNFoxMyZMzF79mzEx8cjJiYG99xzDzIzM9tlhZlbhog1RERERCEtoAHR77//jhtvvBElJSVITEzE6NGjsW3bNiQmJgIAXn75ZahUKkydOhVmsxkTJkzA4sWL22Vs9gYyRKwhIiIiCh0BDYhWrlzp9/qIiAgsWrQIixYtaqcR1XE4lAERM0REREShLKhqiIKJe1E1a4iIiIhCGQMiH6QpM0EA1FxlRkREFNIYEPkgZYjUggCVUD8g4l5mREREoaPJNUQOhwObN2/G999/j+PHj6O6uhqJiYkYOnQosrKymt0ZOthIAY9KJcBLgohF1URERCGk0RmimpoaPPPMM0hLS8OVV16JtWvXoqysDGq1GocPH8b8+fORkZGBK6+8Etu2bWvLMbcLqURILQhQeZ0yYw0RERFRqGh0huicc85BZmYm3n77bYwfPx5arbbeOcePH8eKFStwww034O9//ztuu+22Vh1se5JqiNQqAV5mzDhlRkREFEIaHRB988036N+/v99zunfvjrlz5+Khhx6Sd63vqOQpMwFea4hYVE1ERBQ6Gj1l1lAwpKTVatGrV69mDShYOBQZItYQERERhbZmrTJbt24dtmzZIl9etGgRhgwZgj//+c84e/Zsqw0ukORVZirvq8xYQ0RERBQ6mhUQzZkzByaTCQCwb98+PPjgg7jyyiuRl5eH2bNnt+oAA6VuykyAwGX3REREIa1ZW3fk5eVhwIABAIBPPvkEkyZNwnPPPYfdu3fjyiuvbNUBBkqDU2YMiIiIiEJGszJEOp0O1dXVAIBvv/0Wl19+OQAgPj5ezhx1dMoMkZoZIiIiopDWrAzR6NGjMXv2bIwaNQo7duzARx99BAA4ePAgunXr1qoDDBS3DJGXFJGVRdVEREQho1kZotdffx0ajQYff/wxlixZgq5duwIA1q5diyuuuKJVBxgodqkxo88+RA58vvcE3vn+KESRwREREVFH1qwMUXp6Or788st6x19++eUWDyhYSKvIfPYhsou4b+VeAECvxGhc2i+pPYdHRERErahZGSK1Wo3i4uJ6x0tKSqBWq1s8qGAgrarXqlVei6orLTb5+1U7C9ppVERERNQWmhUQ+ZoiMpvN0Ol0LRpQsJAyRL76ECmfgg2/FcNUa22voREREVEra9KU2cKFCwEAgiDgnXfeQXR0tHyd3W7Hd999h379+rXuCANEWkWmUbn3IVIJgOcCM7PNgf/sPYm/jOzenkMkIiKiVtKkgEiqERJFEW+88Ybb9JhOp0OPHj3wxhtvtO4IA0TqM6Ty6EOk16hRY7XXO3/d/kIGRERERB1UkwKivLw8AMCll16KTz/9FHFxcW0yqGDgUGSIlFNmOo3Ka0B0usLcbmMjIiKi1tWsVWYbN25s7XEEHZvDex8inca97EqnVsFid+BstaVdx0dEREStp9EB0ezZs/H0008jKiqqwf3KXnrppRYPLNDqaojcV5np1O4BUdc4A/LOVKGsxgpRFL3ue0ZERETBrdEB0Z49e2C1WuXvfQmVgMDmY7d7vdY9IOrmCogsNgdqrQ4YdKHRdoCIiCicNDogUk6ThcOUmd217F7jUVTtmSFK6hQBrVqA1S7ibLUFBp2hPYdJREREraBZfYjCgTJDpMx6qVUCNIoIqVOEBkaDs/dSWTV7EREREXVEzSqqrqqqwoIFC5CTk4Pi4mI4pLbOLkePHm2VwQWSXEOkdp8yU6sEqFWCHDBF6dWIjdTiTKUZZTUsrCYiIuqImhUQ3Xrrrdi8eTNuuukmdOnSJWTqhpRsdilD5F5UrRKcGSJpkX20Xou4SC0AZoiIiIg6qmYFRGvXrsVXX32FUaNGtfZ4gobdRx8ijUqARq0C4OxFFK1Xc8qMiIiog2tWDVFcXBzi4+NbeyxBRe5ULQhQJsBUHjVE0REaxEoZIk6ZERERdUjNCoiefvppzJs3D9XV1a09nqDhEOsyRGpFAKQW3C9H6TTylFk5M0REREQdUrOmzP75z3/iyJEjSE5ORo8ePaDVat2u3717d6sMLpDkGiIvRdX1M0TOKTN2qyYiIuqYmhUQTZkypZWHEXx89SFSyTVETtF6DYwGFlUTERF1ZM0KiObPn9/a4wg6vvoQaTwzRHplDREDIiIioo4oaBozLliwAIIg4P7775eP1dbWYtasWUhISEB0dDSmTp2KoqKidhmPr1VmKo8aomi9BnGR0iozTpkRERF1RM0KiFQqFdRqtc+vpvrxxx/x5ptvYvDgwW7HH3jgAXzxxRdYvXo1Nm/ejJMnT+Kaa65pzpCbrC5D5N6HSK2Ce0AUwSkzIiKijq5ZU2Zr1qxxu2y1WrFnzx689957ePLJJ5t0X5WVlZg2bRrefvttPPPMM/Lx8vJyLF26FCtWrMC4ceMAAMuWLUP//v2xbds2jBw5sjlDbzRfGSK1SoBG7bysEgCDVu02ZcYd74mIiDqeZgVEV199db1j1157LQYOHIiPPvoIM2fObPR9zZo1C1dddRWysrLcAqJdu3bBarUiKytLPtavXz+kp6dj69atPgMis9kMs9ksXzaZTI0ei5LNVVTtrCGqO+7sVO1MrEXpNRAEQZ4ys9gcqLHaEalr1tNKREREAdKqNUQjR45ETk5Oo89fuXIldu/ejezs7HrXFRYWQqfTITY21u14cnIyCgsLfd5ndnY2jEaj/JWWltbo8Sj57VTtmjKL1jsDn0idGlpX1ojTZkRERB1PqwVENTU1WLhwIbp27dqo8wsKCnDffffhgw8+QERERGsNA3PnzkV5ebn8VVBQ0Kz7kQIitdq9iFqlaNQoBUSCIHD7DiIiog6sWXM7cXFxbnUyoiiioqICkZGR+Pe//92o+9i1axeKi4tx/vnny8fsdju+++47vP766/j6669hsVhQVlbmliUqKipCSkqKz/vV6/XQ6/VN/6E8yEXVHlt3qIW6GqIofd3TF8cd74mIiDqsZgVEr7zyittllUqFxMREjBgxAnFxcY26j8suuwz79u1zOzZjxgz069cPjzzyCNLS0qDVapGTk4OpU6cCAHJzc5Gfn4/MzMzmDLtJ7Io+RPU7VTsTa50i6p6+WO54T0RE1GE1KyCaPn16ix+4U6dOGDRokNuxqKgoJCQkyMdnzpyJ2bNnIz4+HjExMbjnnnuQmZnZ5ivMgLoMUb0+RIoaoihF8TSnzIiIiDquoF4O9fLLL0OlUmHq1Kkwm82YMGECFi9e3C6PbZf3MnPvQ6Tc7DVakSHSa51ZI4vN3i7jIyIiotYTVAHRpk2b3C5HRERg0aJFWLRoUbuPRZkhEjw6VUs1RNGKGiK16xxXHEVEREQdSNBs3RFs7Io+RO6dqgWoXTVEbgGR6ySHgxERERFRR8OAyAdfNURqlYBIrXN7EqmQWjquvB0RERF1HEE1ZRZoZyrN+HzvSVwztKvPVWYqQcBfR/WATqPCH4akyselKTOHyICIiIioo2nVgOixxx5DYWEh3n333da823az/IdjeH3jYdRYbIoMkQoqRR5NoxLQv0sMnp7ivkJO5coQ2ZkhIiIi6nBaNSA6ceJEsztDBwOpqWJ5jVWuBaqXIVJ537hVwykzIiKiDqtVA6L33nuvNe+u3VltzmDGahd91xD52MmeRdVEREQdF4uqFax258oyi93hUUNUd47axzOmkpfdMyAiIiLqaJodEOXk5GDSpEno1asXevXqhUmTJuHbb79tzbG1O4srILLZHXV7mXn2IfI1ZaZmDREREVFH1ayAaPHixbjiiivQqVMn3HfffbjvvvsQExODK6+8MiBNFFuLlCGy2kW5D5HGI0Ok8REQyRkiBkREREQdTrNqiJ577jm8/PLLuPvuu+Vj9957L0aNGoXnnnsOs2bNarUBtierXaohcs8QeS6790bDVWZEREQdVrMyRGVlZbjiiivqHb/88stRXl7e4kEFSl2GqK6GSKOu35jRGy67JyIi6riaFRD94Q9/wJo1a+od//zzzzFp0qQWDypQLDaphkiETdrcVaWC4LF1hzdqFlUTERF1WI2eMlu4cKH8/YABA/Dss89i06ZNyMzMBABs27YNP/zwAx588MHWH2U78bbKTLm7PeBnykwqquburkRERB1OowOil19+2e1yXFwcDhw4gAMHDsjHYmNj8e677+If//hH642wHSlriKRMj2cNUYNF1cwQERERdTiNDojy8vLachxBwWqvmzKzuzVmrDvH17J7qT8RGzMSERF1PE2uIbJarejVqxd+/fXXthhPQFkURdU21/eefYh8d6p2PpXcuoOIiKjjaXJApNVqUVtb2xZjCTj3PkSKzV0bVVTt/JdTZkRERB1Ps1aZzZo1C88//zxsNltrjyeg6vYyq+tDpFKhUZu7ci8zIiKijqtZjRl//PFH5OTk4JtvvsG5556LqKgot+s//fTTVhlce5NriByeGaKGi6o5ZUZERNRxNSsgio2NxdSpU1t7LAEn1RBZbB57mSnyaL6W3bOomoiIqONqVkC0bNmy1h5HUJAaM5ptdvmYRtXITtWuc5ghIiIi6niavdt9KJKmzGqtDvmYWi24rSxT+3jGpMaMDhZVExERdTjNyhABwMcff4xVq1YhPz8fFovF7brdu3e3eGDtze4QISV3aq3uGSLlLJmvKTPudk9ERNRxNStDtHDhQsyYMQPJycnYs2cPLrzwQiQkJODo0aOYOHFia4+xXUjZIcB92qt+p2rvT5mGRdVEREQdVrMCosWLF+Ott97Ca6+9Bp1Oh4cffhjr16/Hvffe22F3u7coAiIlzz5EPuIhFlUTERF1YM0KiPLz83HRRRcBAAwGAyoqKgAAN910Ez788MPWG107stq8B0QqwX2azFenau5lRkRE1HE1KyBKSUlBaWkpACA9PR3btm0D4NzvTOygAYHVyy71Gte2HUIjOlXLu90zQ0RERNThNCsgGjduHP7zn/8AAGbMmIEHHngA48ePx/XXX48//vGPrTrA9mL1MmUmdaVWBkW+OlWzqJqIiKjjatYqs7feegsOhzOAmDVrFhISEvC///0Pf/jDH3DHHXe06gDbi7caImVXapUgwC6KfjpVMyAiIiLqqJoVEKlUKqgU1cU33HADbrjhhlYbVCB4yxCp3QIiwA5/naoZEBEREXVUjZ4yy8/Pb9IdnzhxosmDCSRpY1clzwwR4G+3exZVExERdVSNDoguuOAC3HHHHfjxxx99nlNeXo63334bgwYNwieffNIqA2wv3qbM1IosWIMBETNEREREHVajp8wOHDiAZ599FuPHj0dERASGDRuG1NRURERE4OzZszhw4AB++eUXnH/++XjhhRdw5ZVXtuW4W523KTONx5QZwICIiIgoFDU6Q5SQkICXXnoJp06dwuuvv44+ffrgzJkzOHToEABg2rRp2LVrF7Zu3droYGjJkiUYPHgwYmJiEBMTg8zMTKxdu1a+vra2Vi7ajo6OxtSpU1FUVNTEH7FxGqohSjZGQKsWkBCl83p76Vw2ZiQiIup4mlxUbTAYcO211+Laa69t8YN369YNCxYsQJ8+fSCKIt577z1cffXV2LNnDwYOHIgHHngAX331FVavXg2j0Yi7774b11xzDX744YcWP7YnrxkidV1AtOLWkSivsSI20n9AxK07iIiIOp5mb+7aGiZPnux2+dlnn8WSJUuwbds2dOvWDUuXLsWKFSswbtw4AMCyZcvQv39/bNu2DSNHjmzVsVi8FFUrM0QpxgikGCN83l7OELGomoiIqMNpVmPGtmC327Fy5UpUVVUhMzMTu3btgtVqRVZWlnxOv379kJ6ejq1bt7b64zdUQ9QQNRszEhERdVgBzRABwL59+5CZmYna2lpER0djzZo1GDBgAPbu3QudTofY2Fi385OTk1FYWOjz/sxmM8xms3zZZDI1ahzea4gaHy9yyoyIiKjjCniGqG/fvti7dy+2b9+Ou+66C9OnT8eBAweafX/Z2dkwGo3yV1paWqNuZ/Gyuau6Cc8Oi6qJiIg6roAHRDqdDr1798awYcOQnZ2N8847D6+++ipSUlJgsVhQVlbmdn5RURFSUlJ83t/cuXNRXl4ufxUUFDRqHC3NEHG3eyIioo6rWQFRdnY23n333XrH3333XTz//PMtGpDD4YDZbMawYcOg1WqRk5MjX5ebm4v8/HxkZmb6vL1er5eX8UtfjWHxsdt9Y3G3eyIioo6rWQHRm2++iX79+tU7PnDgQLzxxhuNvp+5c+fiu+++w7Fjx7Bv3z7MnTsXmzZtwrRp02A0GjFz5kzMnj0bGzduxK5duzBjxgxkZma2+gozoOE+RA1hUTUREVHH1ayi6sLCQnTp0qXe8cTERJw6darR91NcXIybb74Zp06dgtFoxODBg/H1119j/PjxAICXX34ZKpUKU6dOhdlsxoQJE7B48eLmDLlBVi81RE3JEKnkZfeAKIoQfGwCS0RERMGnWQFRWloafvjhB2RkZLgd/+GHH5Camtro+1m6dKnf6yMiIrBo0SIsWrSoOcNskpZmiJTBk90hujV1JCIiouDWrIDotttuw/333w+r1So3TczJycHDDz+MBx98sFUH2F5aWkOkUgZEohj4fgZERETUaM36uz1nzhyUlJTgb3/7GywWCwBnNueRRx7B3LlzW3WA7aXFfYgE9wxRSaUZ8VE6Tp0RERF1AM0qqhYEAc8//zxOnz6Nbdu24aeffkJpaSnmzZvX2uNrNy3uVK04d29+GYY/+y2e+rL5/ZSIiIio/bRoZic6OhoXXHBBa40loLxmiJpQB6QMiA6cMkEUgdzCilYZGxEREbWtgDdmDBbeNndtzl5mAGB2rVjjNh5EREQdAwMiFylDpIyB1E2o/1GpBEin11rtANiTiIiIqKNgQOQiBUSRurpZxKYsuwfqAigpILJ5mYYjIiKi4MOAyKUuIFLLx5raS0gKoGqtnDIjIiLqSBgQuUh9iJQBUZMzRCrPDBEDIiIioo6AAZGLtHWHcspM04Q+RIBiykwuquaUGRERUUfAgMhFmjKL0rcgQ+SaYjOzqJqIiKhDYUDkIgVEBrcMUTOLql0ZIiunzIiIiDoEBkQuUg1RVAtqiFQeNUTMEBEREXUMDIhcLDZnEBPZggyRdD4bMxIREXUsDIhcrF5XmTXt6VEJ7jVELKomIiLqGBgQuch9iPTN70MknS9PmbGGiIiIqENgQOQiB0TauikzVRO27gCUnapdRdXMEBEREXUIDIhcLLb6y+6bWkMkF1XbWFRNRETUkTAgcqmrIWr+XmYaz07VDIiIiIg6BAZELt4aMza1hkguqnZlm0SRWSIiIqKOgAERAIdDlLM5Ldrt3nW+qIiBuNKMiIgo+DEggnvxs9tu980MiJSYISIiIgp+DIjgvsVGS/oQeQuIuH0HERFR8GNAhLqd7oGWdapWe1mmzwwREZF3xaZaPP7ZfuQWVgR6KEQMiIC6gmq1SoBeU/eUNLeGSMlmZw0REZE3//npJN7fdhzLfsgL9FCIGBABgMUVtGjVArSKgKg1aoi49J6IyLtqi7NFSZXrX6JAYkCEujofrVoFrSKoUTWzMaMSp8yIiLyTMujKsgWiQGFAhLopM51aBa26+Rkib+dbOWVGROSV1fWB0cL3SQoCDIhQt22HVq1ymzJrag2Rt73PmCEiIvJOzhAxIKIgwIAIihoijeCW5dE0cdm99wwRAyIiIm+k90czp8woCDAgAmB27U6v16jdpsxaY5UZM0RERN5JnfwtDIgoCDAgQl2GSK9RQa0S5MCmyXuZeV1lxl90IiJvbK4MEafMKBgwIAJgdu1OL/Ugkqa+mrvbvZK07N5qd+DzvSdQbKptyVCJiEKGNGXGDBEFg4AGRNnZ2bjgggvQqVMnJCUlYcqUKcjNzXU7p7a2FrNmzUJCQgKio6MxdepUFBUVteo4pPlrnSsg0qndA6PG8lZULX0CWre/EPet3Ivn1+XWO4eIKBxJGXRmiCgYBDQg2rx5M2bNmoVt27Zh/fr1sFqtuPzyy1FVVSWf88ADD+CLL77A6tWrsXnzZpw8eRLXXHNNq45DCoj0Guc+ZtJUWdNriOofk2qI8kurAQCnK83NHSYRUUixMUNEQUTT8CltZ926dW6Xly9fjqSkJOzatQsXX3wxysvLsXTpUqxYsQLjxo0DACxbtgz9+/fHtm3bMHLkyFYZh9nmPmXWOykaP/9ejlSjoUn3420zWKvrE1BJpcX5WFZ2ZCUiAuoyQ+xDRMEgoAGRp/LycgBAfHw8AGDXrl2wWq3IysqSz+nXrx/S09OxdevWVguIpE8neq0zQ/T+zBGoNNsQF6Vr0v14zRC5PgGVVjkzQ1xeSkTkJNVYMkNEwSBoAiKHw4H7778fo0aNwqBBgwAAhYWF0Ol0iI2NdTs3OTkZhYWFXu/HbDbDbK6bljKZTF7Pq7XasfVICUb2TFBMmTkjmgitGhGu4KgpvO12L82Rl1RZ5MclIiJmiCi4BM0qs1mzZmH//v1YuXJli+4nOzsbRqNR/kpLS/N63r+3HceM5T9i6Zajij5ELXs6vE2ZSZ+ASl0BET8JERE51S27Z782CrygCIjuvvtufPnll9i4cSO6desmH09JSYHFYkFZWZnb+UVFRUhJSfF6X3PnzkV5ebn8VVBQ4PW838/WAABOltfKNUS6FgdE9Y/ZPQIiZoiIiJykDLrdIbKJLQVcQAMiURRx9913Y82aNdiwYQMyMjLcrh82bBi0Wi1ycnLkY7m5ucjPz0dmZqbX+9Tr9YiJiXH78qbSbAMA1Frs9VaZNZe3xoxWuwhRFOUpM9YQERE5KTNDzJ5ToAW0hmjWrFlYsWIFPv/8c3Tq1EmuCzIajTAYDDAajZg5cyZmz56N+Ph4xMTE4J577kFmZmaLC6qrpIDIZq+3yqy5vPUtsjscqLLY5V92BkRERE7KTv4WuwMGtOxDKVFLBDQgWrJkCQDgkksucTu+bNky/PWvfwUAvPzyy1CpVJg6dSrMZjMmTJiAxYsXt/ix5QyR1aFYZdbCKTMvRdVWu4hS15J75+NxyoyICKirIQKYIaLAC2hAJIoNzxlHRERg0aJFWLRoUas+tpQhqmnFKTNvRdV2h4iSqrpVbzaHCJvdAY23giMiojCi7FDNbtUUaGH7V7nK7MzU1Nrs8iqztiiqtjlEuaBawiWmwaHKbMOHO/JxuoLdw4kCweZghoiCR9gGRJVuGaLWqSHyutu93SEXVEtqrfzFDwaf7P4dcz/dhyWbjgR6KERhyW3KjB8UKcDCNiCqsjgDIrPNUa8xY3P52u3eM0MkBWAUWGdctV2lVcwQEQWCcpqMGSIKtKDpVN2eRFFskxoib7vd2x0iymusbseYIQoOXPlHFFhuU2bMEFGAhWWGyGxzyP0vam321ltl5mvKrJIZomDkLyBasukIrn9zK2osfK2I2opbUTU/mFCAhWVAJGWHAOcy+LbsQ+ScMnOfkjEzQxQUpNfdWyuED3fkY3teKX76vaydR0UUPlhDRMEkTAOiuj+AtVaHPIXVJkXVXmqI2IsoOPjLEEnBUmWtrd51RNQ63BozMkNEARaWAVGl2f2PnKnWWePT4j5E3na7t4v1VpmxZiU4SJ9IvU1hSq+RVHxPRK1LFEW3rTvYh4gCLSwDIs8/cqYaKSBq/Roiu8MhZ4jio3QAGBAFC+kTqbcid+k6z+CZiFqH52aufF+kQAvLgMjzj5z0e9nyTtX1A6JKsx3VrsLcLsYIAJwyCxZmm+8MkRwQccqMqE3YPAIiZbaIKBDCMiCq8vGpvy1WmSmX3DNDFFzkGiKPDJHdIcpv1r7+XyGilvGcImMNEQUaAyKFtpgykx5Lp1EhUufMQDFDFBzqpszsXo8DzgwfEbU+m0dGyMJ2JBRgYRkQ+foj1+K9zLwUVUvTcxEalTwlxwxRcDDbva8yU06hVZrdm2oSUeuwOtx/7zynzGosdpwsq2nPIVGYC8uAyFeGSNfCHeiVGSKtWnB7LINOLWeg2JgxOJhdmSGzzQFR9L7JZBUzRERtol6GyGMK7fb3d2L08xtQUFrdnsOiMMaAyEWjEqBpxYAoQuvMBlVIGSKtWj7GrTuCg/INWPm92W3KjDVERG2h/pRZ3e+dKIrYffwsHCJw5HRlew+NwlRYBkQVXv7ItbR+CHBvzGhwBT/SKqUIDTNEwUb5BqwMUhkQEbU9zykz5YeSs9VWVLlW51ZwpSe1k7AMiLxliPTali25B9y37jC4CqhrXNMyEVqVnCHi1h3BQRkQKYNU9ykzvhkTtQV/GaJ8xTQZP5RQe2FA5NIaGSJlUbXBI8CK0DJDFGzcpsms3oMjvhkTtQ3PZffKy8q6oYpaLmyg9hGWAZH0R05Z89PSFWaA+5RZhLeAyNXniBmi4OArCLJwyoyozXy6+3fctHR7vT0efWWIOGVG7UUT6AEEgrRyKCFKh+IK5070rZEh0nipIZIop8xqmSEKCsoMka8aoiqzDaIoQvDSUoGImu79bcexJ78M/VI6uR1X/j7+fpYBEbW/sMwQSVNmCdF6+VhLt+0APDNE7k+t25QZM0QBZ3eIbnspKYMg5SdVq11k3yiiVlTjo1iaGSIKtLAMiKRpkM7ROvlYq9cQ6dwDLINWzcaMQcRzmwCzolu15+vDwmqi1iP97nmu9lX+ThaU1jVkZA0RtZewDIikP3CJygxRC/cxA7z3IVJelrJG3Loj8DwL290yRHb369ickaj1SO9/npkfqajaZnfgRJkyIOIHEmofYRcQORyi3N8iwS1D1PIpM7WfGiK9llt3BJN6GSIfRdUAUMHtO4hajfT+V+mR+ZFqiE6V17pNZ3NhA7WXsAuIqhXZmc6KDFFLt+0A/BdVGxSrzJghCrz6+5d5L6oG6meIKs02OBzuPVSIqHHkgMgj0LHanL9TBWfdt+rglBm1l7ALiKoUS+6NBq18vDWmzFReGjNKIlhDFFQ8901SBqmeGSJlDVFheS0ueOZbzFqxu20HSBSipN+1So+pMGmzZakHUXKM8wMrp8yovYRdQCT9ckXp1G5BS2sXVderIdKo2JgxiHiu9POXIVIWf+47UY4aqx2788+27QCJQpDN7oDNlV2V3oult03pg8jJsloAQL+UGLfziNpaGAZEzvRrpwitW91QW9cQcXPX4OKZITL76EMEuGeIikzON+uz1VaIIqfNiJpC+XtXaZE+nDrb4UlF1dWu46mxEfJt+CGS2kPYBUTlNc6AyGjQtn6GyM+UmUHHrTuCiee0WK3bsnvPVWZ1AZHUyNNic6DawteRqCmUHwalzxPSe6X0Oyn9XilXAecWVuDpLw/gwElTO42UwlFYB0QRiiCotZfd11tlplFs3WFzMLsQYPVXmXlvzAi4p+yLXRkiADhb7b71ABH55+3DYJRHQCQ1bozSaxCtd2aPXttwGEu35OHKhd9j9c6CdhothZuwC4hMPjJEOnXrTpnV71Rdt3WHKNafsqH2Vb8PUeOKqouUAVEVV78QNYW3cgFDvSkz5+9ipE4tB0Q//14mn//IJz+79Skiai1hFxC5ZYgUWZxWyRD5Kao2KLbuAIDvD57B3oKyFj8mNU/9KbP6NUSRroC5yqIMiMzy98wQETVNYzJEUmsUg06DThHOgEj5e+cQgZMMiKgNhG9AFKl1m9ZqjRoiVQNF1Tq1Sl5Rceu/duLmpdvZzyZA6hVVe8kQxUc5G3e6TZlVMCAiai5v+zhGurJA0u9krSJDJAVEEqlfHOv3qC2EXUBUVl2XIVJmhVpjlZmmgT5EgiC4BV6mWhtMbDoWEP4bM7o6mbsCImnKzGp3oKRKERBVMSAiagpvTWkjXR8eLXZnbWW11fn7ZtCpER2hdTu3d1I0AKDG0vBSfH7YpKYKaED03XffYfLkyUhNTYUgCPjss8/crhdFEfPmzUOXLl1gMBiQlZWFQ4cOtegxfU6ZtfYqs3oZIpXrcdyPl/KPakDU39y1flF1vBwQOd/Ez1SaoayFP1vNYJaoKbw1pY3U19VW2hyinP0xaN0zRAlROnm7pYYyRH9dtgNZL23mil5qkoAGRFVVVTjvvPOwaNEir9e/8MILWLhwId544w1s374dUVFRmDBhAmpra72e3xhuy+5bu4bIz+au0mN5FlszIAoMzzfmWpu93nVx0pSZK0NUrKhjAIAyTpkRNYnXgEiRTbfaHfIqs0idGjGKgKhbnEE+119AdLKsBptyT+PomSrkl1T7PI/Ik6bhU9rOxIkTMXHiRK/XiaKIV155Bf/4xz9w9dVXAwD+9a9/ITk5GZ999hluuOGGZj2mMiDSqlVQqwTYHWKr7WU2pk9nVNTa3PZJAwC9KyDSqBgQBQMpC6TXqGC2ObxmiKQpsyOnK/GPz/ZhWPc4t/soZYaIqEm8TZlJjRkB9/5ezhqiuimzbnGR0Lky+TV+AqIdeaXy99L7PVFjBDQg8icvLw+FhYXIysqSjxmNRowYMQJbt271GRCZzWaYzXWf5E0m90ZeymX3gDNzU2m2yQFLSwiCgH/dciEA50oIJSkz5LlclAFRYEhBT4xBi9MVZveialdx56CuRqTHRyK/tBr/3paPz/eedLsPZoiImsZbhkivVUMQ6tqRSMGOQVfXhwhwZohMrgUO/jJE2xkQUTMFbVF1YWEhACA5OdnteHJysnydN9nZ2TAajfJXWlqa2/XlHgFRXW1P6zwVgiBAEAQoZs+gEuAzA1XKP6oBYbE731ClGgW3Zfeu7ztH67HhwbHIvuZcAHWrzbrGGgBwlRlRU3nLEGlVgvz+WGOxyx9IIj1qiNymzKzKlZ+1OFRUIV/enlcif8+AiJoiaAOi5po7dy7Ky8vlr4KCuq6mVrsDVa5PFlJAFBupc7vcWgRBgFbtjIqkFWbelFa6/1H99ZQJVy38Hhtzi1t1POROCnpiXCl5bxkinUYFjVqF64Z1k6fPAKB/l04A2JiRqKm8ZYg0apUcECkDGIOXKTMpIFJOmd30zg5cufB7FFfU4nSFGUdPV8nXMSCipgjagCglJQUAUFRU5Ha8qKhIvs4bvV6PmJgYty+JSfHLEeMKgJ6dMgjzJg1Av5ROrTl8AHVF1soC6+uHOzNWo3onAKifIfrmlyL8ctKELzymZ6h1SUGPMkP04te/Yc2e32F2fYqVsoYatQoTz637f66v6/8VZoiImsbbqi+tWpBrg6QARiU4f/88M0QGj4DI7hBxqLgCVruIw8WV+PFYqdt9MyCipgjagCgjIwMpKSnIycmRj5lMJmzfvh2ZmZnNuk/pl6OTXiMHKyN6JuCW0Rk+MzgtIRVQK/dMe/aPg/D9w5fi6vO6AqhfQyT1JapqRJ8Naj5lDRHgrO1atPEInvrigFuGSDJ5cKr8fd8UZ5BdbbF7nQIgIu+8bd2hUQnQemSIDK6seidFDVHXOIPcs0jqZn222iLXaxaZahkQUYsEtKi6srIShw8fli/n5eVh7969iI+PR3p6Ou6//34888wz6NOnDzIyMvD4448jNTUVU6ZMadbjSb8cMa08PeaLRpoyUywr1ahVSIuPRG6hc87bs7lfhRQQmfmHti3JAZFHJ9yz1VavPaMu6BGP7gmR+P1sDYamxcqrE8uqrUgxtrwgnygceMsQadQq+cOH1DhX2t8s3tV3KKmTHpE6DSJdx6UM0ZnKugU0p8prceyMc7oso3MU8s5UMSCiJgloQLRz505ceuml8uXZs2cDAKZPn47ly5fj4YcfRlVVFW6//XaUlZVh9OjRWLduHSIiIpr1eJ4F1W1N6lwd4aULttTjxnPKTCrcZYaodeWdqcIjH/+MOy/piXH9kuVahpiI+v8vSJ9ilRkilUrAittG4kyFGWnxkYiL1OJMpQVnqy1IMTbv/0eicONt6w5vU2ZSrVDf5E54+Iq+6O/KyhrkPkTO98fTiq10ispr8ftZ5yreAV1ikHemyq1MgqghAQ2ILrnkEoii7/bqgiDgqaeewlNPPdUqj9f+AZFrysxL00epSNezqFqeMjMzIGpNa/efwo5jpTD8T+0WEHnulaTkuTKwa6xBXmEWG6mTAyIiahyvGSKVSq7Xk0oIpIBIEAT87ZLe8rlSg1tfGSKprcmA1Bh8te8UM0TUJEFbQ9QW2jsg8lZULZEyRFUedShyhohTZq1K+iQpLc+V6oT8TZ/6614e71qdeLbKClEUUVheC5u9/qdfIqrjLUOkUQuIc/0+nSp3BjSee0FKPDtVn6mo+0DyW2GFfFxaCcqAiJoivAKi6nbOELlqiDz3NQOctSvSlJoyy8Aps5bZdrQEOz0KKwHgjCsTd7K8FqZaKyyuT6repswk/rqXx0Y6b/dqzkEMe+ZbjMzOwexVP7Vk6EQhr9brKjOV/AHxhGvKK9JHQGTwDIgUGaL8Uuc2HYmd9Ejq5JzGZkBETRFeAZGUIYps5xoiLwGRIAh1dURVyoCIU2bNVWOxY/q7O/CXpdvrPX9nFLUGh4oq5SkzX59EAf/NOqVNJg8WVcqv37ajJT7PJyIfGSKVgHjXe7JUA2TQep/KlouqXVl1ZQ2RpFucQf7Qy4CImiI8A6J2riHyNfUiTbsoAyJTjfMPudUuyiuhbHYHNvxWJBcSkndnKs0w2xyotTqw/0S523WnK5UBUYX83Op8BD06tcpvK4YbL0zHmD6dccfFPbHsrxcAAIorzFyGT+SHt8aMygxRieu9sKEpM6mGSPl7Lekaa5CnwmutDu54T43GgKgN+ZsyA4B4jwxRrbWubT1QlyVasukIblm+E29uPtqWw+3wlFOPP//uHhApU+sHiyrdNnf1pqGtXAZ3i8X7M0dg7pX9cUnfRPmN+qTHXnVEVMfbBwaNWpDfCyWRPt4z5caMVjscDlGeClfqFheJTnoNpM8zzBJRYzEgakP+psyA+gGRVD8kqbLYIIoi1uw5AcC5rQf5dlax+/ze38vk7y02h9zfBAAOFlXIgaevwMdX5sgbQRDQLc65+kxK+RNRfWYv/b80KpVcVC3xlSFSfristdnlDzrK413jDFCpBLk+kEvvqbEYELWhulVmPqbMoqSVSs6ASFpyL6ky23GwqBJHXc3GTpXXttVQQ0KZW4aoDK/lHMJF2Tn4SREcAa6ASJoyU3t/423qZr/d4iIBMCAi8keavlLWcWrVgttegYCfompF4FNptskfJgek1m3RJH04YR0RNVVYBUTSdEZ7NdLTqKWtO7z/cnvOm3vLEK3df0q+LC1JJe+UXb8LSmuwcMMhnCyvxUc/Ojf4lXoOFVeYUewqxlRmgsb06Sx/35QMEQBFhqi6eYMnCgNS01Plh1KNooZI4isgUqkE+QPmybJa2F37dgzoogiIYhkQUfOETUBUXmOFyRVwSH+82po0ZeYr/St9KpJqXyrqZYhsWLuvUL58ptLCol0/lFNmgLMwHQC2HnGu/uqeEIlUVzAsvZHqNSq8d8uFuG1MBh4Yf45826YGRGmuDFEBM0REPkkZolhDXQCkUdWvIZK27vBGWmkmLbOPj9K5vad3ZYaImimgnarb0wnXJ/eEKJ38C9XWpAyR3kcNkeeye2mFmeRQUSVyiyqgUQlQqQRYbA4UlteiR+eoNhx1x+Wra7TUvbZztB5dYw04qZh61GlUGHtOIsaek4iC0rrsjt5HVs8XZoiIGiZniNymzFRyXy+Jr4Uoyuuk39fO0To56698f5cDomoGRNQ4YZMhkv4odouPbLfHTO6kBwA5K+EpoV5Rtfsv7qFiZ1flbnEG+Q/uSU6b+SRliDJcAaNnrVjnaD3G9El0O6bMBCk/pTZ9yow1REQNkWuI3KbMBOg1akQrdrb3NWWmvO54ibO2snO0HgNTjQCAwd2M8nkxcoaI7UqoccIoQ+T8Q5XWTtNlAPCPqwbgD0NScVGvzl6vj5P7EDn/kHvWEEkp4YRoPQxaNY6ersLJMhZW+yIVVd9+cU+crbbgknOScOXC7+XrEzvpcbGfgChSp4Zeo4LZ5mhGUbXz/6vTrl5EvlYWEoUrURTlVWaxioBI6+rXFhelRaWr1Yi/hqlSQJQvZ4j06J0UjU0PXYKkGL18HqfMqKnCLkOU1o4ZImOkFmP6JMqrzTzFK2qIHA6xXoboeEndNF8XV5bpFPvc+CRNmSXH6PG3S3pjQGqMW3auc7Qe6QmRbpkgZeAjCHW1DE3NEMVGahHleqM+wdeIqB6L3QFpL2/PDBFQ16gW8J8hipCnzOqmwgGgR+cot3IIBkTUVGETEEm1HVLxazCIi3L+wtodIipqbXLRt0RaFZcQrUdqLKfMGnLWlWmLVbyx9kqKlr9PdE1hjuwZLx/z3K9Myto1NUPk7EXEaTMiX5RdqpU1Q1JApFxp1pgpM+mDRxcfJQkMiKipwiYgOuGaakqLb78ps4Yo581LqsxyHyKt6w3CtRAKnaN1SI11/tKfLKvFN78U4sBJNmn0JE2ZKZu89VYERJ1d+4+NyEiQj3luzyHtUaZrYlE1UDdtpizOJiInaR8zQXDfVFmaMlNmbn3tZQag3qIYX4tMpIDodAXLDKhxwiggkmqIgidDBLhPm0k1RNJOzRLnlJnzj+3WoyW4/f1duOuDXX7vt6TSjMv+uQkvrz/YBqMOPhabA1Wu/Y3iFJ8+lQFRoiu1fu2wbuieEImx57jXEzlv6wqI/Ox074tUzJ1bWNHk2xKFOqlliF6jcquxa+qUmWd9UY8E7+/pg7sZoRKAn34vx88ezVmJvAmbgMhsdUAQgC6x7dOUsbHk5oyVFrnFvGcKWDllJnVYPl5S7TcVvPVoCY6crpKbEoY6KTuk8vj02Tux/pRZlF6DjQ9eIm/KqiQFqL425PXnvLRYAKjXGZuI6qbM9Bq12++XVi0VVSsyRI2YMgOc2SZfdaFp8ZGYMqQrAGBhzuHmD5zCRtgERACQEhPR5P4ybS3BS4bIs5N2gmLKTOlQke9MhLQCo9BUK28SG8qkJfexkTqoFEXs/VJiEKFVISFK51bIqXL1dvIkbQHQsxm9noa4AqJfT5nYQJPIg7TkPkKrcuszJDWwjW9kQKS8LtVo8Lui8+5xvaESgG9/LcL+E+U+zyMCwiwgCrbpMqBuiqakyoIKs/OPekqMe/DTOVqPSJ2mXl+dXD8BkbQCAwDyXHuhhTJphZlngzdjpBaf3HURVt2ZWa9eyJvrhnXD5jmXYObojCaPoVucAQlROljtIjfiJfIgNWXUa9RuQYy0CldZ++drt3vndXU1RN19TJdJeiZG42pXlujpLw9AlJa5EXkRVgFRtyAqqJbEu1aana2yyJ2q62WIXJ+cxvTpDKNBi6z+yQCcnax9URb2Hg2HgKiqfkG1ZGCqEb0UU2f+CIKA7glRjQqevN1WmjbbW1DW5NsThTIpQ6SsIdKqBfl3TW55oVbJXf69Mejqruue0HAm96EJfRGhVWF7Xim+2neqwfMpfIVVQNQ/Jabhk9pZfJSzrqWk0iI3JZMKqAFnTYy0jPy1G4di29zLcMWgFAD+i3cLFFtIHD3tO3AKFdKUWZxHhqi9ndctFgDwEwMiIjfSKrMIrRpdjBHoHK1z25RVqp2UVnr6otznzFdBtVLXWAPuGtsbAPDcV7/KdZhEnsKmUzUADOpqbPikdiZliH4vq5E3HFVmiOKjdHJKWRAEGHRqnJPszHZIW3t4sjtEuTM3ABw9HZoZor0FZbDaHbigR7xiysz/m2lbG5IeC8C5soWI6nhmiDbPuVQuqAacRdCv3TjUZ18hiXI6rTEZIgC4Y2xP/GvrMZwsr8XO46U+dw+g8BZWGaJBXYM3QyTV+ahVgtwvBwASovT1biMtJT9TaUFJpRmAc5WV1e785HOqvAY2R91ceSjWEJ2tsuCGt7biT29uxc5jpYoeRIHOEDmD7rwzVVj1YwFrFohc5FVmrhVmUXpNvY7wk89LxfAe8fVuq6RcZdajc+PqQiO0aozt62yzsTn3dKPHTOElbAKiHgmR6BQR2D+W3kgZotMVzsCmT1I0ohSbHHpLH0fqNEh3LTU9WFSJdfsLceFzObj9XzsB1K0wk3rpHD1dGXJ/mL/8+SRqrc6tAB5c/ZP8Mwc6QxQbqcMfhzqLOB/+5Ge8voHLfSl8/XffKdz5/i6UVVvklZcRLVzpq1xl1j2+8atBL+mbBADYxICIfAibgEhaTh1s4j0yQKN7d3bb9Tkhun6GCIA8bfbS+lzc8+FuWGwOfH/oDGqtdvzuWmE2rHscVAJQZbGj2BVwhYpP95wA4OxDcrykGl//UgTAfeluoPzzuvNw96XOmoWVYdIHisjTb4Um3P/RXqz7pRCrdhagyuyaMmtGjy8l6YNtF2OE3+X5ni7u0xkqwbk69yT3GyQvwiYgGhisAZFHRmNUn87Qa1SQWuQk+PgDPzQ9DgDw47GzsNqd2R+bQ8QvJ8vlbEnPxCi5admXP5/CkRAprs47U4U9+WVQCcAr1w+BTqNClE6NzJ4JuKxfUqCHB5VKwO1jewJwdkiXVsDlnanCg6t+wqbc4kAOj6jN1VrtuH/lXrmA+buDZ/DjsVIAaPSKT1/O62bETSO74+9X9W/S7WIjdXKvsM0HmSWi+sKmqHpgl+ArqAaAThEaqFUC7A4RWrWAERnxEAQBUXoNKmptbvVESreN6YmenaPwa2EFDFo1duSVYGPuaezJL5NXmKXHR6Jn5ygcL6nG018ewNNfAnMm9MXtF/eEQxTdmlQWm2qxLa8UVw5K8bvkNdBEUcTb3x8FAIzpk4irh3TFxEFd3JbvBoOYCC16JETiWEk19p8sh0oQ8LcPdqO8xopvDhRi40OX4FRZLT78MR/bjpbgrrG9cN3wtEAPm6hVrNyRj98KKxCt16DSbMOOY6VyA8bxA5JbdN8atQpPTxnUrNte0jcJu/PLsCm3GDdemN6icVDoCZuAqF+QZohUKgFxkTqcqTTj/PQ4eePCaFdA5GvKTKdRYeK5XTDx3C4AAIcoOgOigjI5HZwWH4m4KB22HD6DuEgdiivMePHrXLz4dS5UAjBxUBfcc1lv9EuJwT8+249vDhThl7E9MXdi0z55tRer3YEHV/2E//x0EgDw5xHONzTPwsxgMbCrEcdKqrH9aCmW/+8YKs02qASgotaGmct/xP6TJnll4ZLNRxgQUUgQRREfbM8H4PwA9ubmIzhZXgsLgOQYPc4N4GrfS/om4qX1B/HD4RJYbI6gfe+gwAib/xuUdTnBRiqsHt27bimotJLC15SZp6GuVPDOY6U47GrYmB4fiT8NT0Pu0xOx4+9ZeGbKIPkNwCECX+07heuWbMWp8hp8f+gMAOCd7/OCtsX9Rz8W4D8/nYRGJeC5P56LCQNTAj0kv6Q3/vdcwVD3hEj8e+YIAM5l+XaHiKz+SVAJztYIrGugULDz+FkcKq6EQavGH8/viosVmyhn9U8OaCZ3UKoRnaN1qDTbsOv42YCNg4JT2AREwWxcv2QYDVpMPi9VPjZpcCoyOkc1uARVcm43IwQBKDKZUWG2oWdiFPqldAIAec+uv4zsjr3zxmPP4+Px33vHoGdiFCrMNjzxn19Q41oBYneImPPxzzDV+t44tj3ll1Rj2Q95KK+x4r3/HQMAPDqxn5wdCmaDUp0BUYWr4ebV56Xiot6dceOF6RAE4J5xvfH2zcPl7tZbXEEpUUe2wpUdmnxeF8REaN0CopZOl7WUSiXg4j7O8Ww6yFo+cseAKAg8OrEf9s4bjx6KDUUfGH8ONj50SaNXTXWK0KJPUl2x4rxJA7zWAkXqNIiL0mFAagz+7JpDl1ZojeuXhLhILX49ZcK0t7e7bf8hiiJ+OHwGR9p4Cb8oisj5tQj/2noM3/xSiKsXbcGTXxzAH17fgkPFlYjUqfGnCzrG1JJn3ysp4H3uj4Ow9/HL8eDlfSEIAsa4MoPfH2ZARB1PeY0VC3MOYeexUvxUUIavfnZuj/HnEd0BAKN6ObccSo7RI7NXQiCHCgBu/YhqrXas3lmAvy7bgZU78kOuPQk1TfDOI4WZ1kgjD+sej4NFlcjqnyz33PBn0uBUPPvfXyG9B0w9vxtmjz8HN7+7A/tOlGPMCxtxblcjXv/zUKzdX4gFa38D4JyKW/Tn89ElNgKrdhZgSLdYjOyZ4HX3eF+KK2rxywkTkmMisP9kOf697Tg0KgEO0fs+YMdLquUxxgRhPylvYiN16BZnwO9na9AvpRP6JDszdoIgwKhoIDm6TyIWbjiMHw6fgdXugFoQmvRcErW3j37Mx2sbDuOaoV3x7a/FOHDKhFe+BWIMWljsDowfkCw3KTVGarH2vjHQqAS3hRyBcnGfRKgE4LfCCly0YANKXatAN+WexrajJbggIx5D0+KCtlULtR1BDPGQ2GQywWg0ory8HDExof0/eJGpFh/v+h1/vjAdcY3MLN3w1lZsO1oKlQDsefxyGCO1OFxcgcfW7MfOY6VwiEC/lE7IL61GtcUOjUqAzSEiJkIDg06NIpOzv1GqMQKDuhrRN6UTeidFQ69Ro6TKjN9OVaBnYhRuvDAdEVo1jp2pwkOrf8JOP/P3EVoVzusWi935ZzGmTyJuuCANf/tgNxyiiG8euBi9kzq1yvPVHu5buQef7z2Jh6/oi79d0tvrORabA0Oe+kZ+fiN1alzUqzP6d4nBwNQYjOuXhFqbHesPFCGzVwKSOvnf2oCoJWqtdrzz/VEcOV2Fwd2MuGJQitv+ipVmGy7KzoGp1iYf02tUcifqPknR+PRvFwVlI1zJHxf/gD35ZQCc712jenfGJ7t/h9TgXyUA917WB0aDFrmFFYjUadA9IRLDusfhbLUFpVUWJEbrkRQTgeQYPaL1GgiCAFOtFXqNKigCv1DQ3n+/O0RAtGjRIrz44osoLCzEeeedh9deew0XXnhho24bTgFRc3z0Yz4e+WQfLuwRj1V3ZrpdV1BajUmvbUF5jbOeaHj3OCz96wW4ZfmPckFi11gDymus8sa0vnQxRmB078745kCRfH+9EqNQUmVBhEaNW0b3QFKnCBSZajHpvFR0jTXA7hChEpwZlT35Z1FjtXe4PYiKTLXY8Fsxrh3WzW3fJk/3frhHXj3n6cpzU/D72Rr8/Hs54iK1+MdVA9C5kx4WmwMWmwNmmx0GrRrDesQhqVME7A4RFpsDGrUArVoFq92BgtJqdIuL5KqaMGe22WG2ORAToYXV7kBheS26xhqgUgkwu4Lul9cfxBHF/od6jQozR2fgnnF9YNCp8ebmI8he+xu6xhqg16gAAVg6/QJsP1qCzQdPY+7E/khvxKargbQ7/yxWbM9HVv8kZPVPhkatwv+OnMEnu07gRFk1th0tbdL9GbRqaNUCTLU2aFQCeidFw6BTQwCgEgTYRREOh7PVSYROjQiNClq1CiqVAJUAqAVn2xC1yrl9kyAIzkyxANc5AtSuf1WCv3PgOkc63/n+KQiAAOlfyJchXxbqjivOhed1HvcBt8v178Pz/gEv9+/6j+cxQRBQVWnCuMEZDIgkH330EW6++Wa88cYbGDFiBF555RWsXr0aubm5SEpqeFqIAZF/DoeIT3b/jhEZCV7fxNbtP4U7/70bggB8cfdoDOpqRKXZhnmf7Ud0hAaPXNEPggDsyS9DbmEFDhZV4OiZKtgdIqL0GvROjMba/adwqrxWvs/z0mKxeNr56BprqPd44arKbMOBUyZ0MUbgdIUZW4+W4EhxFf7z0wm58WZjSD2tAECjEtArMRonympQabYhQqtCz87Ris2CnbeRJ+cEAWrBWY+mEpyZAIcI+U1VJdS92Qry985/bQ4RNocDVpsIjdqZ5RJFwOoQYXc4t1iRbgNAvg9B+b3rTbPusvsx5Zu7WhCgUTvHolEJ8h8Etzdd6bL8c3r5gyB/7/zG25t23Zu/+9S21/vwcrzKbEd5jVVeqJDYSQ8BzjYSKpVz/BqVSv457A4RBq0aXWIjUGW2obTKCrvDAZtDhN0h1v1rF+EQnc+7dFm6XgSQ1EmPCK0Kx85Uw2xzwOZwYEdeKaotdgzuZsSJszUoqbKgf5cYXNAjDl/8dBJnq+vGeN2wbth2tAS7XZmU/l1i8NDl5+DRT/fhdIUZL147GNcNT4MoikHVA6w1rPqxAC98nYueiVEYmREPs92BX06Y8FNBGRJj9EiM1qOkyoIiUy0qav1/GKTmc5irUfDKnxgQSUaMGIELLrgAr7/+OgDA4XAgLS0N99xzDx599NEGb8+AqOW++vkU9BoVspq5QqTWasfG34rx6ykToiM0uDmzByK0TCk3xo68Uty9YjdiDFosnnY+1uw5gZxfi6BWqaDXqKDTOP89U2nBb4Um+PptlqY6ifxJiYnAdcO7YeboDMRG6iCKIr45UITHPt2HEletDeCcZto051JmHAFUW2woNplhtjmQFm9AaZUFh4or5YBVFEU5Y2O2OVBrtaPGapevtzuc/zpE5ypfh8P1vSurJP2rPEcUXcdF54dau+s2znNE1zmQ71sEXO8NzuPOy3XHpcuQL7vOU37vvLnbZc/7kN6ARK/3A9f913+cutuIiscBrDWV+GHeZAZEAGCxWBAZGYmPP/4YU6ZMkY9Pnz4dZWVl+Pzzz+vdxmw2w2yu27fLZDIhLS2NARF1WI0ttDbVWlFjsUOnVkGrUaGs2oLcwgokdtJjYKoRR09X4vezzl5Hrredujcj1782h+jMZIhAtKuLuigq32zhuizC4YDrDd+ZsteqnZkOm8OBGotd/iOgUdeNWxTrbuNQvPlJ30vXSW+MDod03HVM+kMhirDbnX8U7Io/CPLP5vVNvv6bMOD+pi7dXnk+FM+V5znStaLHY9cdFxGp1yAmQosYgwaiWLeRs06jgsOV0ZGyPHaHCI1ahYpaGwrLaxCt1yA+Sg+dpi4bplapoFG7nluV4ri67npRFFFkqkW1xY6MzlGI1GlgczhwfnocOkfrseXwGXSO1mFgqhHv/e8YTpXX4qrBKbi4T6LX1amnymvw2Kf7cKi4EoO7GXHvZX3QL4Xvp9S22juhEdSrzM6cOQO73Y7kZPfMRHJyMn777Tevt8nOzsaTTz7ZHsMjahf+ao+UYiK0bivwovUadIurmwbtk1y30o3C27XDusnfPzShb4PndzEasGxG4+o2iTqqkMt3zp07F+Xl5fJXQQF3GyciIiL/gjpD1LlzZ6jVahQVFbkdLyoqQkqK920b9Ho99Hrv+38REREReRPUGSKdTodhw4YhJydHPuZwOJCTk4PMzEw/tyQiIiJqvKDOEAHA7NmzMX36dAwfPhwXXnghXnnlFVRVVWHGjBmBHhoRERGFiKAPiK6//nqcPn0a8+bNQ2FhIYYMGYJ169bVK7QmIiIiaq6gXnbfGtiHiIiIqONp77/fQV1DRERERNQeGBARERFR2GNARERERGGPARERERGFPQZEREREFPYYEBEREVHYY0BEREREYY8BEREREYW9oO9U3VJS30mTyRTgkRAREVFjSX+326t/dMgHRCUlJQCAtLS0AI+EiIiImqqkpARGo7HNHyfkA6L4+HgAQH5+frs8oeSbyWRCWloaCgoKuI1KgPG1CC58PYIHX4vgUV5ejvT0dPnveFsL+YBIpXKWSRmNRv7PHSRiYmL4WgQJvhbBha9H8OBrETykv+Nt/jjt8ihEREREQYwBEREREYW9kA+I9Ho95s+fD71eH+ihhD2+FsGDr0Vw4esRPPhaBI/2fi0Esb3WsxEREREFqZDPEBERERE1hAERERERhT0GRERERBT2GBARERFR2AvpgGjRokXo0aMHIiIiMGLECOzYsSPQQ+rwvvvuO0yePBmpqakQBAGfffaZ2/WiKGLevHno0qULDAYDsrKycOjQIbdzSktLMW3aNMTExCA2NhYzZ85EZWWl2zk///wzxowZg4iICKSlpeGFF15o6x+tw8nOzsYFF1yATp06ISkpCVOmTEFubq7bObW1tZg1axYSEhIQHR2NqVOnoqioyO2c/Px8XHXVVYiMjERSUhLmzJkDm83mds6mTZtw/vnnQ6/Xo3fv3li+fHlb/3gdypIlSzB48GC5mV9mZibWrl0rX8/XIXAWLFgAQRBw//33y8f4erSPJ554AoIguH3169dPvj7oXgcxRK1cuVLU6XTiu+++K/7yyy/ibbfdJsbGxopFRUWBHlqH9t///lf8+9//Ln766aciAHHNmjVu1y9YsEA0Go3iZ599Jv7000/iH/7wBzEjI0OsqamRz7niiivE8847T9y2bZv4/fffi7179xZvvPFG+fry8nIxOTlZnDZtmrh//37xww8/FA0Gg/jmm2+214/ZIUyYMEFctmyZuH//fnHv3r3ilVdeKaanp4uVlZXyOXfeeaeYlpYm5uTkiDt37hRHjhwpXnTRRfL1NptNHDRokJiVlSXu2bNH/O9//yt27txZnDt3rnzO0aNHxcjISHH27NnigQMHxNdee01Uq9XiunXr2vXnDWb/+c9/xK+++ko8ePCgmJubKz722GOiVqsV9+/fL4oiX4dA2bFjh9ijRw9x8ODB4n333Scf5+vRPubPny8OHDhQPHXqlPx1+vRp+fpgex1CNiC68MILxVmzZsmX7Xa7mJqaKmZnZwdwVKHFMyByOBxiSkqK+OKLL8rHysrKRL1eL3744YeiKIrigQMHRADijz/+KJ+zdu1aURAE8cSJE6IoiuLixYvFuLg40Ww2y+c88sgjYt++fdv4J+rYiouLRQDi5s2bRVF0PvdarVZcvXq1fM6vv/4qAhC3bt0qiqIzwFWpVGJhYaF8zpIlS8SYmBj5+X/44YfFgQMHuj3W9ddfL06YMKGtf6QOLS4uTnznnXf4OgRIRUWF2KdPH3H9+vXi2LFj5YCIr0f7mT9/vnjeeed5vS4YX4eQnDKzWCzYtWsXsrKy5GMqlQpZWVnYunVrAEcW2vLy8lBYWOj2vBuNRowYMUJ+3rdu3YrY2FgMHz5cPicrKwsqlQrbt2+Xz7n44ouh0+nkcyZMmIDc3FycPXu2nX6ajqe8vBxA3YbGu3btgtVqdXs9+vXrh/T0dLfX49xzz0VycrJ8zoQJE2AymfDLL7/I5yjvQzqHv0ve2e12rFy5ElVVVcjMzOTrECCzZs3CVVddVe854+vRvg4dOoTU1FT07NkT06ZNQ35+PoDgfB1CMiA6c+YM7Ha725MIAMnJySgsLAzQqEKf9Nz6e94LCwuRlJTkdr1Go0F8fLzbOd7uQ/kY5M7hcOD+++/HqFGjMGjQIADO50qn0yE2NtbtXM/Xo6Hn2tc5JpMJNTU1bfHjdEj79u1DdHQ09Ho97rzzTqxZswYDBgzg6xAAK1euxO7du5GdnV3vOr4e7WfEiBFYvnw51q1bhyVLliAvLw9jxoxBRUVFUL4OIb/bPVE4mDVrFvbv348tW7YEeihhq2/fvti7dy/Ky8vx8ccfY/r06di8eXOghxV2CgoKcN9992H9+vWIiIgI9HDC2sSJE+XvBw8ejBEjRqB79+5YtWoVDAZDAEfmXUhmiDp37gy1Wl2vWr2oqAgpKSkBGlXok55bf897SkoKiouL3a632WwoLS11O8fbfSgfg+rcfffd+PLLL7Fx40Z069ZNPp6SkgKLxYKysjK38z1fj4aea1/nxMTEBOWbWqDodDr07t0bw4YNQ3Z2Ns477zy8+uqrfB3a2a5du1BcXIzzzz8fGo0GGo0GmzdvxsKFC6HRaJCcnMzXI0BiY2Nxzjnn4PDhw0H5exGSAZFOp8OwYcOQk5MjH3M4HMjJyUFmZmYARxbaMjIykJKS4va8m0wmbN++XX7eMzMzUVZWhl27dsnnbNiwAQ6HAyNGjJDP+e6772C1WuVz1q9fj759+yIuLq6dfprgJ4oi7r77bqxZswYbNmxARkaG2/XDhg2DVqt1ez1yc3ORn5/v9nrs27fPLUhdv349YmJiMGDAAPkc5X1I5/B3yT+HwwGz2czXoZ1ddtll2LdvH/bu3St/DR8+HNOmTZO/5+sRGJWVlThy5Ai6dOkSnL8XTS7D7iBWrlwp6vV6cfny5eKBAwfE22+/XYyNjXWrVqemq6ioEPfs2SPu2bNHBCC+9NJL4p49e8Tjx4+Louhcdh8bGyt+/vnn4s8//yxeffXVXpfdDx06VNy+fbu4ZcsWsU+fPm7L7svKysTk5GTxpptuEvfv3y+uXLlSjIyM5LJ7D3fddZdoNBrFTZs2uS1rra6uls+58847xfT0dHHDhg3izp07xczMTDEzM1O+XlrWevnll4t79+4V161bJyYmJnpd1jpnzhzx119/FRctWsTlxR4effRRcfPmzWJeXp74888/i48++qgoCIL4zTffiKLI1yHQlKvMRJGvR3t58MEHxU2bNol5eXniDz/8IGZlZYmdO3cWi4uLRVEMvtchZAMiURTF1157TUxPTxd1Op144YUXitu2bQv0kDq8jRs3igDqfU2fPl0URefS+8cff1xMTk4W9Xq9eNlll4m5ublu91FSUiLeeOONYnR0tBgTEyPOmDFDrKiocDvnp59+EkePHi3q9Xqxa9eu4oIFC9rrR+wwvL0OAMRly5bJ59TU1Ih/+9vfxLi4ODEyMlL84x//KJ46dcrtfo4dOyZOnDhRNBgMYufOncUHH3xQtFqtbuds3LhRHDJkiKjT6cSePXu6PQaJ4i233CJ2795d1Ol0YmJionjZZZfJwZAo8nUINM+AiK9H+7j++uvFLl26iDqdTuzatat4/fXXi4cPH5avD7bXQRBFUWx6XomIiIgodIRkDRERERFRUzAgIiIiorDHgIiIiIjCHgMiIiIiCnsMiIiIiCjsMSAiIiKisMeAiIiIiMIeAyIiIiIKewyIiCholJSUICkpCceOHQv0UOpZt24dhgwZAofDEeihEFEbYEBEREHj2WefxdVXX40ePXq0+2Pfe++9GDZsGPR6PYYMGVLv+iuuuAJarRYffPBBu4+NiNoeAyIiCgrV1dVYunQpZs6cGbAx3HLLLbj++ut9Xv/Xv/4VCxcubMcREVF7YUBERK1izZo10Ov1OH36dL3r5syZg6FDh/q9/X//+1/o9XqMHDnS7fgvv/yCSZMmISYmBp06dcKYMWNw5MgRAM4AZcqUKXjuueeQnJyM2NhYPPXUU7DZbJgzZw7i4+PRrVs3LFu2rMHxL1y4ELNmzULPnj19njN58mTs3LlTfnwiCh0MiIioVbz33nuYOHEiEhMT61134403Yu/evfj555993v7777/HsGHD3I6dOHECF198MfR6PTZs2IBdu3bhlltugc1mk8/ZsGEDTp48ie+++w4vvfQS5s+fj0mTJiEuLg7bt2/HnXfeiTvuuAO///57i3/G9PR0JCcn4/vvv2/xfRFRcGFAREQt5nA4sGHDBowfP97r9eeffz6MRiPWr1/v8z6OHz+O1NRUt2OLFi2C0WjEypUrMXz4cJxzzjmYMWMG+vbtK58THx+PhQsXom/fvrjlllvQt29fVFdX47HHHkOfPn0wd+5c6HQ6bNmypVV+1tTUVBw/frxV7ouIggcDIiJqseLiYlRUVKB3797ysZtuugkff/yxfLlXr144dOiQz/uoqalBRESE27G9e/dizJgx0Gq1Pm83cOBAqFR1b2XJyck499xz5ctqtRoJCQkoLi4GAEycOBHR0dGIjo7GwIEDG/9DuhgMBlRXVzf5dkQU3DSBHgARhQ7lVNbevXvd6oZUKhUEQfB5286dO+Ps2bNuxwwGQ4OP6RksCYLg9Zi0XP6dd95BTU2N19s2RmlpqddpQSLq2JghIqIWS0pKQmRkJA4ePAgAMJvNOHToEA4cOADAOaWWl5fndzn90KFD5fMlgwcPxvfffw+r1dpqY+3atSt69+6N3r17o3v37k26bW1tLY4cOdJggTgRdTwMiIioxVQqFa6++mq8/vrr2L17N5566imMGzcOa9aswdatW7FgwQKUlpZi8uTJPu9jwoQJ+OWXX9yyRHfffTdMJhNuuOEG7Ny5E4cOHcL777+P3NzcVv8ZDh8+jL1796KwsBA1NTXYu3cv9u7dC4vFIp+zbds26PV6ZGZmtvrjE1FgccqMiFrFq6++iptvvhkXXXQRzjnnHHz88cf48ssvMXHiRERERGDJkiUYMGCAz9ufe+65OP/887Fq1SrccccdAICEhARs2LABc+bMwdixY6FWqzFkyBCMGjWq1cd/6623YvPmzfJlKQukzGx9+OGHmDZtGiIjI1v98YkosARRFMVAD4KICAC++uorzJkzB/v373crlA4GZ86cQd++fbFz505kZGQEejhE1MqYISKioHHVVVfh0KFDOHHiBNLS0gI9HDfHjh3D4sWLGQwRhShmiIiIiCjsBVdOmoiIiCgAGBARERFR2GNARERERGGPARERERGFPQZEREREFPYYEBEREVHYY0BEREREYY8BEREREYU9BkREREQU9v4fvYES2rlnYRIAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "ax = plots[2].plot(xlim=(0, 5000), title=\"IR Spectrum\")\n", "ax;" ] }, { "cell_type": "code", "execution_count": null, "id": "04a03c24-fe01-41f0-89fd-e55c554e9b48", "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.8.16" } }, "nbformat": 4, "nbformat_minor": 5 }