{ "cells": [ { "cell_type": "markdown", "id": "a1db2bad-37f3-4ed7-a47e-96b9dd0209da", "metadata": {}, "source": "# Different roughness models" }, { "cell_type": "markdown", "id": "ed694b6c-5543-4e65-996d-c7a622a10f01", "metadata": {}, "source": [ "`refnx` has a variety of different interfacial roughness models available, over and above the default [Nevot-Croce][1] roughness formalism. The Nevot-Croce formalism is equivalent to a Gaussian roughness between two layers.\n", "\n", "When a different roughness model is used the calculation no longer uses the Nevot-Croce damping factor. Instead, the entire scattering length density profile is micro-sliced.\n", "\n", "[1]: L. Nevot and P. Croce, Rev. Phys. Appl. 15, 761 (1980)\n", "\n", "The definitions for all the roughness models are round in the following papers:\n", "\n", "Svechnikov, M.; Pariev, D.; Nechay, A.; Salashchenko, N.; Chkhalo, N.; Vainer, Y. & Gaman, D.,\n", "\"Extended model for the reconstruction of periodic multilayers from extreme ultraviolet and X-ray reflectivity data\",\n", "Journal of Applied Crystallography, 2017, 50, 1428-1440\n", "\n", "Stearns, D. G. J. Appl. Phys., 1989, 65, 491–506.\n", "D. Bahr; W. Press; R. Jebasinski; S. Mantl, Phys. Rev. B, 1993, 47 (8), 4385\n" ] }, { "cell_type": "code", "execution_count": 1, "id": "3429ae27-a228-4f67-b3b9-02220bf79620", "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt\n", "\n", "from refnx.reflect import Slab, SLD\n", "\n", "# the following are the different interfacial roughnesses that are available\n", "from refnx.reflect import Erf, Linear, Exponential, Tanh, Sinusoidal, Step" ] }, { "cell_type": "markdown", "id": "6baae64d-0a18-4722-b02f-034efacb3e3c", "metadata": {}, "source": [ "A comparison of the different roughness models can be achieved with only the fronting and backing slabs, with a roughness in between" ] }, { "cell_type": "code", "execution_count": 2, "id": "425947b9-e481-4d56-a121-08b78a76232b", "metadata": {}, "outputs": [], "source": [ "air = SLD(0.0)\n", "d2o = SLD(6.36)" ] }, { "cell_type": "markdown", "id": "09dda6fd-d6bf-4bd0-b166-8cb9c856845d", "metadata": {}, "source": [ "Let's create SLD plots corresponding to each of the different roughness models. In the following graph we're going to offset the SLD values vertically, to permit an easier comparison between the plots. Please note the horizontal offset we apply to the default roughness model, `rough_model=None`, which corresponds to Nevot-Croce roughness. This horizontal offset is because of the microslicing that's necessary for all the other roughness models and has no impact on the overall reflectivity calculation.\n", "Whilst not shown in this particular example it's possible to have a different type of roughness for each interface in the system." ] }, { "cell_type": "code", "execution_count": 3, "id": "7e422a82-7c77-47a5-9b4e-120c42de5626", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjEAAAG0CAYAAADD+nYCAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjgsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvwVt1zgAAAAlwSFlzAAAPYQAAD2EBqD+naQAAn1VJREFUeJztnQd8FOXWxp/t6b13IEDovSugIoiKHVGxi4q9995791659i4qWNBPUEREkN4JPRBSIL23zdbvd97Z3eyml91kd3P+944z887szEs22Xn2VJnZbDaDYRiGYRjGw5D39AQYhmEYhmE6A4sYhmEYhmE8EhYxDMMwDMN4JCxiGIZhGIbxSFjEMAzDMAzjkbCIYRiGYRjGI2ERwzAMwzCMR6KEl2IymXDixAkEBgZCJpP19HQYhmEYhmkHVL6uqqoKcXFxkMvlvVPEkIBJTEzs6WkwDMMwDNMJcnJykJCQ0DtFDFlgrD+EoKCgnp4OwzAMwzDtoLKyUhghrM/xXilirC4kEjAsYhiGYRjGs2hPKAgH9jIMwzAM45GwiGEYhmEYxiNhEcMwDMMwjEfCIoZhGIZhGI+ERQzDMAzDMB4JixiGYRiGYTwSFjEMwzAMw3gkLGIYhmEYhvFIWMQwDMMwDOORsIhhGIZhGMYjYRHDMAzDMIxHwiKGYRiGYRiPxGsbQDIM410Yq3QwG0w9PQ2GYVyMoVLb7nNZxDAM4/ZU/ZOLit8ye3oaDMN0A1X1Ne0+l0UMwzBujy6nStqQy6SFYRjvxdj+SBcWMQzDuD1mk1msQ87ph4CJsT09HYZhXEhlZSXwavvO5cBehmHcH4uI4U8shmHsYRHDMIz7Y9EwMhm7khiGaYBFDMMwHuNO4ngYhmHsYRHDMIz7Y2YRwzBMU1jEMAzj/lgsMexNYhjGHhYxDMO4PexOYhimOVjEMAzj/lgK9cq4RgzDMHawiGEYxnNiYtifxDCMHSxiGIbxIHdST8+EYRh3gj8SGIZxfywahlOsGYaxh0UMwzAelJ3Exe4YhmmARQzDMG4Pu5MYhmkOFjEMw7g/HNjLMEwzsIhhGMb94RRrhmGaQdncIMMwjDvR24vd1dfWoqqkqMEixTBeTFVVVbvPZRHDMIznBPZ6oYiprayAXltn268pL8fGHxYjP+OQ2DebzdBWt/9DnWE8Ha1e3+5zWcQwDONBMTHwGkicrPvmM2xetrRdFhafgEDI5BwBwPQCdCxiGIbxIsyWmBhPdSeZjEas+eIj5OzdbRsz6PUoyzsutlUaH9s4WZv6jByLMWefZxv3DwmFb2BQD8ycYbqfyspK3Pv59+06ly0xDMO4P2bPdif9+92X2L58WdMDMhlOv/4WDD/tjJ6YFsN4PCxiGIZxf6yBvR6iYapLS/Dtkw+isrjQZokhpl+5ABGJKbbzgqNjEBId02PzZBhPh0UMwzBuj6e5k9Z8+THKC/IaBmQyjD/nQow567yenBbDeB0sYhiG8Rx3khu2Hcjdl470v/+E2aK0jAYDDq7/RwiXuY8+i9C4eCiUKvgFBff0VBnG62ARwzCM++OmdWL0unr8+tZLqCkva3Js+KmzkDR0RI/Mi2F6CyxiGIZxe9y1d1L66pVCwASER2D07HNs4yq1BkOmn9ajc2OY3gCLGIZh3B83yk7a/PMSZGzdKLZLc3PEesK5czFy1lk9PDOG6X2wiGEYxv2xBvb2cExMeUE+1n7zmUNxusDwSAw95fQenRfD9FZYxDAM4xmuJKKHLTG7Vv4mBEzcwMEYN+cCMRbTrz+UanWPzotheissYhiGcW/srB7daYihoN0D69Y49DWiGBhi3DkXInXshO6bDMMwzcIihmEYz3AldbMlZvNP32Pj0sVNxgMjItF39NhumwfDMC3DIoZhGLdvlNjdIsZsMmHvmlViO3n4KNF8UdxeLsew02ZBLld0yzwYhmkdFjEMw7g3djEx3VXsLnd/OqqKi6Dx88d59z3GMS8M46awiGEYxr1xCOx13W2Ks49hx4pfYTQaUJh5RIwNmDiFBQzDuDEsYhiG8ZzsJBdaYlZ/9j6y03c7jA2exgXrmPb9jpro91T837qmA3buUNuYuSFW3f5862n251su4nDMtm13zO5PxPp6ey+s9XoNm40ONj7X/lCTY+aWj7dyHfu5tXWfqqpKtBcWMQzDuDfWDzeZ64rd1VZWIGdvutiedNGlUKjUCImORULaEJfcr7dCDzF9vdH2vtbXGVBTUe/w8CMxUFepg15nOc+CQWeCtloHo0ESDCQcrOKB2lZZt00GEwwGE4x6U6PzpHNoDmLfaJauZTSJB2jDMWmeYswqOOzH6bV0zNhwf8a51Olq2n0uixiGYdwb60PChVaYjC0bRQPHqJR+mDx3vsvu44nYHuRCSBhRlFWF2iqdwzl1VTqU5tUKQUDo6owoL6iF0WByuE5tha5BxDDNQ2JdrGWWtXVMOtAwJv09yBq9xm5lGbRuNvr7sT/W2p+W3cEmpzlcv5XXtXr9pnNUduB3hEUMwzBujc0E7cR4mBOHDqCyqMC2n776D1sMTG+goqhWCA0rZE3IyyhHcW61g1VEW6vHicPl0GtdIzzkChn8gtWQ21vYZDL4Baqg9lE6POAUSjl8A1RiLVPIhFVOTg96hUy8nh6UNEbHlWq5dJ7cckwsUmC4dZ/uTefQWhyXWa4hazjXuhYiweG49HpxPfH6huMO4gP24w1CpDmh4o4d2nuKyspKLHy9feeyiGEYxr2xfJl31od8UVYmvnn8vuac/eg/YTI8GbKE2IsTin0gEXJke5HNKlJRWIeS49Wdvod/iAbBkb4OY2pfJcLj/aHSSKnnJA5Cov3EuD0+/ioEhGpsokVOIsIN+mExnguLGIZh3BtbB2vnPOwOb14vBExAWDhCY+Nt4wmDhiIsLgGeZKHKPVCGmvJ6sV9dVo/dq3NQV6Vv87VkPfANdGyVEBLli4S0MChUDSYvhVKGmL7BCIqQRAvpSBImbDVg3AUWMQzDeIY7yUmWmKPbt4r15IvnY9gpM+EJnMgoR+auYofsjsJjlcjLqGjztT4BKgyaHGsTImQtSR4aLqwiDOPpsIhhGMYjLDEUd9BVasrLUHD0sNjuO2oc3BESKhlbC5F7qEzs19fohTuoOZQqOeL6h4iYCnLL9BkZiYHjY0SciBVrHAfDeCMsYhiGcW+sxodOuJNyD+xF+l8rYTZJcSLVZSViHd23P/xDQtHTaGv0OLgxHzqtwTZWcrwGR7YXNjm3/9goBIb72PaVagXSJsUiMKxhjGF6GyxiGIZxa6geh6AT1oQ/P/gPSnKzm4ynjpuI7kZXZ0BZQa1tn+qY/P31QZTlNa2JQf/UodMS4BckuXzI2hLXv+dFF8O4GyxiGIZxbyxxIB0tdFdVUiwEjEwmx0mXXimaNxJqXz8MOnk6upOinCr8/OYO1Nc0WFys+AerkTIi0rZP/8zUsdGSm4hhGM8UMUajEU8++SS+/PJL5OfnIy4uDldffTUeffRR9u8yTG9sO9DBmJis3TvEOrpfKsafexG6EyrotunnoyjLl6wshVlVQsBo/JVQaxo+doOjfHHKFWkICndMWWYYxsNFzEsvvYT33nsPn332GYYMGYKtW7fimmuuQXBwMG6//faenh7DMG4eE5O1Z6dYJw8bhe4OzP3ri/0iONeeiMQAnHf3aGga1U5hGKbzuO1f0/r163HuuefirLPOEvspKSn45ptvsHnz5p6eGsMwPZGd1I6YGArkzdq1XWwf27lNrFOGu1bE/PPtIez5O7dJ8zvKFpp8USp8/JRQqKS0ZmsxOIZhvFzETJ48Ge+//z4OHTqEAQMGYNeuXVi3bh1ef735WsT19fVisS9bzDBM73EnmYxG/PTy06ivaQiUVfv6InbAQKfPydoj6OjOYuxZndvkOBWMm3rJAAyeEuf0ezMM4wEi5sEHHxRCJC0tDQqFQsTIPPfcc5g/v/nmbC+88AKeeuqpbp8nwzAupp3F7gqOZggBQ4G7Q6adJsb6jh4HhdJ5Rd2ox9Cv7+5Czr5Sh/FRM5Mw6vQkh/RntrowTC8WMd999x2++uorfP311yImZufOnbjzzjtFgO9VV13V5PyHHnoId999t22fBFBiYmI3z5phGJf1TmojJiZn3x6xTho6HKdec6NL3oh96040ETCx/YIx4Zy+ol8QwzDdi9uKmPvuu09YYy655BKxP2zYMGRlZQmLS3MiRqPRiIVhGG/tYt2GiNm7W6wThwx3yn0ps2j9D0dg0DU0VCw4JrmpJ1+YKkr5Exo/7iXEMD2F24qY2tpaW10HK+RWMpksX8sYhukd2IrdtXKKwYDjB/aJ7cTBw5xy238WHxINFhsTmRSIEacmQK5gy0tvRwhsWoxGadtoJJ8jzPScsq6t51gW6TXi1Y7HpAtaNs2tvK7pcbP9NS3Xsa0t27YvA7Sy9eBq5nz7tf3r7E5vsmNue9vhOi1cxkpNTbXni5g5c+aIGJikpCThTtqxY4cI6r322mt7emoMw/Rwsbvs9F346ZVnodfWOZzqExiEiMTkTt7GjNoKHYxGEwqPVQkBI1fKMP2ygbbOzpQhlTAwlAVMe36eej3M9FC3ex+NZWUwlFqEodkMU00NjGWljg87o0mcZ6qtldpFGEkMGAGDsdG+Aab6epi19dK97I9R4LWJzjfDbNBLx3V6SWzQcQoWJ5FhNlnOtQgOOk5ztooSqxhpbttefDBOpdr+98ZTRcw777yDxx57DDfffDMKCwtFLMyNN96Ixx9/vKenxjBMN0LPGYGdiNm7ZlUTAUMMPvkUyBpZcNvLxp+PYvuKLIexoVPjMWhy78gwooe3sbTU4RsziYz6gwfFWpyj16P+6FGYKhy7ZxvLK6A7niuEhnhdXR0MhYX8kG8OqSOnwyJrbbyDr4HlRbZXy9pat3F+i2P2m628pqXXO/xMHHdVej2QITVq9VgRExgYiDfffFMsDMP0YprJTjp+UHIdnX3ng0gcPFQ6LJfDNzCoU7eoq9Jh96ocsU1WF7pTUKQvxs5Ogaejz8tD7ZYtMOt0Yt9QXIzaLVuF0LBhMqH+yBGYqqpcOheZSgVFWBgV0RH7ch8fKMPDKVbA7iQZFCEhkPv7QaZQAgo5ZHIFZEoFQGty41n2ZRofyH00gEolzrGeK9Z0TZkcMqUSMrUaMpVSbIt26HJZw3ESxzQf67UVtKZ70TEalzecR7+DYl9md8yyTa+zHrfftvxbuZN4+xElUoKDPVvEMAzDOBS7sxhYqktLUFGQL3oipYwYDY2fX6eaMdZVSw91Iv2fEzDoTSLeZe5DYz3igUPulqpVq2AokTpzmyqrULVyJfQnTjScRC6b2oamk+3CzpJFokMzYIAkNCzH1ImJUEZFOYhKuZ8v1ElJkFmSK+h1qrg4yH0d2ynIfH07bSljmOZgEcMwjFtj825Y3ElWK0xkcp9OCZjK4jp888xmGOqb+t1Hz0p2CwFjrK6BPkfqvm3SalHxw4+o3bzZwdVjLCmxuXlaRSaDz/BhUIaG2QSH75gxkhCxQxUbB5+0gZK1gmE8BP5tZRjGIywxlIFUWVSIY5a2AvGDBnfqctQigASMXCGzBewSsf1C0HdUQzfp7sBYXY2qP1bCZBffYywuQemXX8LUjqrjqqQk+A4fbnNp+I0fB78xYxysKYqgIOGeYRhvhEUMwzAeIWKy9+7E2r8esw3HDxzS4UsZ9Ebs35Ants+4cRj6DI9Ad0AWlKrly6Hdf8Bu0ISKX36FoaCg2dcogoNt7hmfwYMRetmlkAcG2o7LqTZWWhq7Z5heDYsYhmHcGqsLhdYKCtKUKxASE4uUEW03djy6owh71uTarqGrM6K+xoCAMI1oyOiyORuNqFj2C/S5Ul8lbXo6qtesafZcVUICfIbYCTK5DAEnT0XwuedIAaYMw7QIixiGYdwag1YKwDXDhAXvfowAS2xHe/ocrfnmIGorGwJ4rQybliC6TDvLJVS7aZNDTZSKH39C9erVjieqVAg57zwRk2JFGROL0EvmNQmAZRimfbCIYRjGraksLBANrBVqdbsFDJF7oFQIGI2/EtMuGWirRUGNGZMGt/86raHLPY7sq66C/vjxJscorTf43HOl1F6VSmyTW4hhGOfBIoZhGLemoqAAoQiBj108SHs4tEmKNek/Nhr9x0U7ZS4mnQ55Dz6Iqr9W24q/UXVXZWSkCLK1Ig/wR+RNN8F35Ein3JdhmOZhEcMwjFtTVVQoRIxvUNuF7P759hCObCsU23XVerEeOCGmU/clgVK1ejWM5eW2seq/16D6r78czqM6KokffABVtGPKMsMwrodFDMMwboNOW4elzz+B8vyGgm1JsoFICh0A3+DWRUxVqRZ7VkuBtFYiEgMQ3afjVXwpEPjEAw+g8rflTQ+qVEh48w34pKWJXWVsLGcIMUwPwSKGYRi3IWv3DpywFLOzYak+3pY7KXNXkVhHpQThlMsHiu2QKL92F6/THjqE8m+/g9logLG0DFV//AEolQiYOtVWnZayhUIuvhgBJ03pxL+OYRhnwyKGYRi3IT/jkFgPmHQyJl4wT2wbtlbCuLkC8jbSjY9sl0RM/7FRiEjoWPwMlebPveVW6HOk/klWou+/D2FXXtnBfwXDMN0FixiGYdyGPIuISR42EpFJUvPFykNZqEQFRIqSBep7tHNlDvTW1gFmM/IypNiVviMj23QVVSxdCp1dRpF2714hYJQxMQiZe5EY0/Tvj8DTT3f6v5FhGOfBIoZhGLfAZDKi4OhhsR3bX3IHEbZ2QXZuoR2/Z2PHSqm3kD1RyYEIimi95krlsmXIe7Sh8q89sc89i4Ap7CpiGE+BRQzDMG5B6fFc6OrqoNL4IDwhsZku1g0i5li61Ll5wPhom2ih4+RKamx1EWX9TSZpX69H4auviW3/aVOhtruP74jhnRIw5dpybCvYBhOkezAM0zVqqtrR2NQCixiGYXoMs8mEI9s2o66q0hYPE90vFXK5oomIsXaxriypQ1lejTDMnDxvAHz8VS1eX2QYLfulybgqMREJb78t+g91Bb1RjwV/LMDBsoNdug7DMA0Y65p2mG8JFjEMw/QYhzevxy9vvOgwFtNvgMO+te+RNSYme2+pdF7f4FYFjHbfPpuAoeq5VqjEf8yTT3RawNTqa3Gk/IjY/v3Y70LABKgCMCDUcd4Mw3QOXa0O+7G/XeeyiGEYpsfI2Zcu1qFxCQiNiYXa1w+jzpjjeJKpIfso91AFyvJrxX5SMw0cddnZqNu5U2yXL1kq1kFnn434V19x2pyvWH4FDpVJViMrj096HLP7zHbaPRimN1NZWYnFWNyuc1nEMAzTY1gDeSdddCkGTZnW/EkWd1JFcR1ycy2+chnQZ0SE42laLY5dNh/G4uKGQbkcETff7Lz51hQIASODDHEBcWJsWsI0nJFyhtPuwTBM+2ERwzBMj2A0GFB0LFNsx/RNbfE8k8ESlAtg6iUDRENHCuYNjwtwOK/il1+EgFEEB8NnyBAxFnDaqdD07eO0OaeXSJaj/qH9sfQcydLDMEzPwSKGYZgeoSQ3Gwa9Dho/f4REx7Z4Xk1FvWhArfRRYui0eFGBl+JkyHVkqquznVf2+ediHX7jjQi/9hqXzDm9WBIxwyKGueT6DMN0DBYxDMP0CPlHJFdSdN/UVnsP1ZTVg2wuQRE+thYCFT/8iLxHHmlyrszPDyEXXeiyOe8p3iPWQyOGuuweDMO0HxYxDMN0C5k7top0aisnDh8Q6+h+/R3Oyz9agYqiBgtLdYlWEjGRDUXsSr/8UqzlwcGQqaUMJZlcISwwinZ0u+4MJrMJe4v3im22xDCMe8AihmGYbqnG++tbL0NXJ2UW2RPXX+oGTZScqMbSV7ZJATAWRvkpALXcJmK0+/ejfv9+yFQqpP6+AoqQEJfMubiuGB/s/gA1eimYuN5Yj2p9NXwUPugX0s8l92QYpmOwiGEYxuWUnTguBIxSo8G4OQ3uHv+QUPQbM962n7u/TAgYv2A1wuP8xVhAmRbQGmDYuRmFB5ehbscOafy001wmYIiv9n+Frw983WR8ROQIKOX80ckw7gD/JTIM43IKMqXicFEp/TB57mUtnpd3RGriOPyUBIw5Q2oAWfLNAdTtKkLNurXQH/3Ldm7IBee7dM67i3aL9Zl9zsTAMKmXk0KmwGlJp7n0vgzDtB8WMQzDuJzCzAyxju7TshuGMo7yMirEdmy/BguLuV5nPQEh8+ZBplFDnZgE/5NPdtl8jSYj9pZI8S8Lhi0QKdUMw7gfLGIYhnE5BUeP2DKRWqKyuA61lTrIlTJEpQTaxvWFhQA0UEaGI/apJ7vl3TpWeUzEwvgqfdE3uG+33JNhmI7DIoZhGJc3eSw8ZnEnNbLEmExmZGwrgLbagOLcKjEWEalE3do1tnMMx/Mg802Bz+DB3fZOWVOpB4cPhsK+GSXDMG4FixiGYZyKtqYaO1b8Ap2lEJ2+vl5sK1VqhMcnOpxLAmblR/scxnw2/4bcb36y7ftOvAVKX8BnqOtFDHWlNsGEXUW7xD6nUjOMe8MihmEYp7Jj+S9Y//1XTcaj+qZCrnC0apw4JAXyhsf7IzTGH8bMw4hbvxqK8HCoExLEMXmQFB+jjIl26Tv1zYFv8OLmF0U9GCtc1I5h3BsWMQzDOJX8I1KH55SRYxCRmCy2qSLv4JOmNzm34FilWI87qw/6jY7Csfkvok5XgcjbnkDoJZeIY0Uf7UH94XLI5FK1Xlfxy5FfHARMlG8UJsRMcOk9GYbpGixiGIZxKoXHjor1xPPnIT6tZReQXmdEyXGpkFx0nyAYKytRt3On2A+wzzyydLF2pYjRGXU4UCpVEP7hnB9Eh2qNQsP1YJp947SAvhYwGQGz0bI2OW7T2qQHjDq742ZpTYWAmt1HG8fFCQ1rwmy/bma7u7DOo/tu2M23M3fv/aqbFsVsCRYxDMM4jdrKClSXlgAyGSJTWu8eXZRdBbPJDF+NCdofvkFV1jHAaIS6Xz+o4uNt59mMIy4UMYfKDkFv0iNUE4rUkFRbjyavoLoIOLEdMOqlfW0FULgP0Nc5PhSr8oGqPMfX0mtozFAv7ZMgMdi/jmFcQH37RROLGIZhnEahpahdaEwc1D4NvY6aoyBTciX55+5G0e8f2MYDpk5t/lugC4WFfWNHjxEwNSXA3h+AmiJpv74KOLwSqDzueB5ZTVyFTG5ZFABlcdm25YBCAyjU0rb1POpHLrZljttN9q3nNtpuskbTfftzupNu/73x4n9fLQnun9t1KosYhmGc7kqKSmlaW2XNNwdxeEuBgzuJCKo6Bp+hQ6Hu2wdyf3/RxNEBmzvJdW9UenG6e2UjkSAplrp8C4vJlg+BzDV2ZinLOSZD+64XmQZoLI0xlRogegjg06hlg38EEJwgiRArJEACYgC1f8ODzDcUUAfYiRGGcTKVlcANwe06lUUMwzBdwqDX48TBfTAaDMjavaPZejA6rQHp/xxv4sqXmU2IKNmLqKeehf+kSa274534wFydvRo5VTm2/U15m8R6SMQQdCtGA5CzEaivbhiryAFWPwfUlbX9+tgRQMI46Vs5WUISxwNxox0Vn0+QJDwYxgthEcMwTJdY+9Un2L58mcNYY0tMSW611NgxSI3z7h4lxvT5+ci7Yh7U0MF35MiWb2CxxDgrJmZ/yX7cvvr2JuMyyFybUq2rAUolS5Ut3uTPJ4DMf5o/3zcMUPlJ2/GjgUm3SGNWyKISksTWEKZXwyKGYZgukbNPiicJjY2HysdHxMMkDnF0yxRmS9V4o5IDRT0YovzfnVBTaf/RoyH3bSV+xsnZSdsLt4t1fEC86EhtZVzMOIT52IkEZ3J8O/DlhUBdadNjKn8gKq1hn9w5aWdJokWhcs18GMZLYBHDMEynMRr0KMmV3DIXPfoMgiKimj2v2CJiFP/8jIxl90uvLZMK3fmNJ3dIK9gCe50bxHte6nlYOGIhXEJtKfDrXQ1xLWXHAH2NFJeishNsIcnAOW8DUYNcMw+G8XJYxDAM02lIwJiMBmj8/REYHtnieYVZUiaSX2469CXZDscCTz211Xs4O8V6b/Fe5wbxkluo5IjdRM3Ab/cBWf86npc0GZj/HaBpaG7JMEzXYBHDMEzXs5GS+7aYmmzQGVGWL6X5hqprkPx1Q0sCJbUXSJaq+rbpTnJCYG9FfYXoUE04Jf7FZAK+vhg48lfTY+pA4Lz/SNYXpY8UgKvgj1yGcSb8F8UwTKcpsoiYSLtAXrK6/PnJPui0RlunavIIqXSVCBnSD36jR3foHmarO0nRdRGzt0SywiQFJiFY074UziZoK4GSDGn70ApJwMiVjhlAFIB75stA36atFhiGcR4sYhiG6TRFWZliHZncUJ33wIZ8m+XFnvDSffCbPLzjN7FmJ3VCw/x4+Ed8uvdTW0+kSl1l16wwVLPl/elAqVTUz8ZZrwFjru7cNRmG6TQsYhiGaRd6XT1+ef0FlOc3lKavKMxvklJdnCsF8U48ry+SBoeL7ZwF18EnOx0+D3/a8Z+2JdSkM9lJH+z5wKEejJUp8VPQKda9KQkYpS/gJ/3bMGAWMPqqzl2PYZguwSKGYZh2kbN3NzJ3bG0y7hccgvCERJvrp+S4VLgtNGcLlGUGmLT18M3eI6q/+g7tuAXE5k7qYExMubbcJmD+d/r/4KPwEdv+Kn8MCB3Q9gVytwLLbnMs21+RK63Pfw8Ycn6H5sMwjPNhEcMwTLsoOia5jlJGjMaE8y+2jYfFJ0KhlOqZVJVooaszQmYyoO71p1BPnY0taAYOhNzPUrytI3Sy2F16idRKICUoBZPjJnfsniSclt8vNUpsTL/TgMHndex6DMO4BBYxDMN0KP4lYfAwJAxq3qJSTJV5ydpRmw+/QQOhSkiQDijkCL3k0s79pM2d651k39Sx3fehjs2UMp2zCTi+TXIbXfZtQ+VcmgT1HeKeQQzjFrCIYRimXRRlS6nJUXZBvI2xupICqnMR9cC9LfZD6gi28isdFA7Wpo7tFjF/PQusfdVxbPz1QN9pHbovwzDdB4sYhmHaFdRbduJ4k0wka3PH3ANlMBpMyN5TJMYCqo9Dk2ZXSr8rdMCdRHVgDpQe6Hhn6rpyYON70jbVdKFUqNAUYMqdXZk5wzAuhkUMwzBtUpKTDbPZBN/AIPiHOvYX2vDjEaSvkQSOlWB1LZShTuqcbCt213YA8BXLr0BmheT2IpRyJQaGDWz7Hju+lNoCRA4Cbt7A7iInQTWC9CYT9EYz9AZam6AzWvZp22CCyWyG0WQWb7PZbpvGpUW6jnWbjovzzI6vscZ/mxsFhDs0TredYzlmd7DhdY7nOI7ZX8vx+g7Xst9hOkxdjZTh2B5YxDAM0yz0QUxxMHWVlchK3ynGIpNTmlTOzTtSIdYRiQGQV5ZAuX8LYpNaaejYQdpb7K6gtkAIGOpG3S+knxg7s8+Z0Cg0jifWlQGfzQHK7NofkIAhJi7stQKGfs7V9QaU1+pRUacX6/I6nW2/pt6AWp0RtTppXSe2G/Zp0eqNFpFigsFohsFqRWOYDmCqb1pnqiVYxDAM0yxZu7Zj6QtPOIw1diWZjCaU5UsCYPaNw1Dz1vMoP7QEvqc5sbGitU5MG+LCGshLlpfv53zf8ok7vwHypXMdCE4EhjVkXXmTOCmr1SO3rBb5FVoUVNWjsFLrsF1UVY/yOr2waLgShVwGlYIWuViUYl8udCMdk8tk0rZM2paLMTTZpuONX0NrK9ZN+98Y6++P45htq43XWseaOc/uRNtxhzGmo+hqq/FBO89lEcMwTLPk7JMe9L5BwQgICYXazw/DTp3lcE5FUR1MBjMUxnrkTh0HmUlKqfZJc2JXZnP7YmLalY1E1yLXEXHqY461XoLiHDtMexhkGTlcWIWD+VXILK5BVmktskpqkFVSiyqtod3XUSvlCPVTIcRXjWCxViHYVwV/jRL+GgX81Er4qmitgJ9GCT/Ltq9lUVsECl1HEisNooVEB8O0RWVlJT64Hu2CRQzDMK1mI02eOx8jZ57Z7DklxyUrjH9Nnk3AKEJC4DdurPNcSVbjQBvPv2YDeeurpKJ11no1lXlA4V6AXEzjrnPsd+RBlNXosCOnDHtyK3Egv1ISLiU1DnEZjYkK1CA22AdRQT6IDtIgJsi67SOOhfqpEeKngo9K0Z3/FIbpEixiGIZpVcREJqW0+BMqOWGpC1NzHIkfvC8K2imCgyHXNIpD6SzW9Oo22g4YTUbsLd7b1BKzdIHUpLExaWd5jIAhIXe4sBpbj5Vhe7a0HC2yxPA0IsxfjbSYQPSLDEByuB+SwvyQEuGPxFA/YSVhGG+DRQzDME2oq65CdUmx2I5oTcQcKxfrgJo8+I4aDUWAv3N/mvamBTsRU1BTgHd3vosaS0BuvbEetYZa+Cp90S9YCupFdRFw+A9pO3pYgyVHHQBMvQ/uTGGVFv9mFGPtoWKszSgWMSuN6Rvpj5GJIRgcG4SBMYFiiQzQtBk7xDDeBIsYhmknZpMZxR/tQX1W+9P/PBVKp74w+W7xQCx6fof9AYc00yEGE4YEKyEfPRd5L+50xUwaNu0ezl/t/wo/ZfzU5OxRUaOgkFssDvuXSZXy4kYBN/wNd4fiWFak52PF3nzsypHEoRUflRyjEkMxJjkUo5NDxHaov7rH5sow7gKLGIZpJ6ZqHeot6cTeDskFpVzqhwSDnU+nUWgKZYlIg8om5zkTZbQfZKqGvgO7i3eL9QX9L8CgMCmIWC6TY1r8yUDlCemk9B+ktRv3OTpRXocfdxzHsp0ncLDAURwPjQ/CSamRmNo/AmNSQqFRsjuIYRrDIoZhOlr+XiFDzH3OCVx1F3L37kFW+i7bfn7GIZQez8GoM8/B2LOlDJ7MnUVY++1hqH2VCI6UsngMBfnwP7IFwyeGIOKWm102P0Wg2hYTYzAZsK9Easx45eArbTVhhOvpk9lA9gbHFw85z+2yiH7fm48l23Lx75Fim8eM0o0n9QvHrCExmDk4WgTdMgzTOixiGKajlWPlMihDvOcBYzaZ8OtHr6K+pmmwaERaX9u/taRchzoz0HdMFKbPl1oKZC+4HjWH10FzxVPd9jM5WnEUdYY6+Kv8RYdqGyUZDQJGbvloG3qh1D7ADaA6LV9szMLizTmieJyVCX3CcOHoBMwcEo0QP3YRMYzXiJjjx4/jgQcewPLly1FbW4vU1FR88sknGDvWu74FMx6CtRCYlwVOVhQWCAGjUKkw7pwLbeP+waHoO7rhb600TxI5mpy9KP1sk9jW7pMsIprU1G6brzWVenD44Ib4F+LQ79K673Tgyp/hLmw9VooP12bij335tl+hhFBfXDQmQYiXxDBLh2yGYbxHxJSVlWHKlCk45ZRThIiJjIzE4cOHEeqsfiwM09ny915WsKsoR0qlDk9IwpSLL2/xvNJsKdjU9NMXKCg/2HBAJoOmb1+Xz7NKVyVcSdsLtjdf1M6aSj3gDLgDG4+W4K0/D2PD0RLb2JTUcFw9uQ9OTYviwm8M480i5qWXXkJiYqKwvFjp08ex5DnD9Iw7ybt+7sVZ1nowLf99GfUmVJaRC0SGYF8dgk4623bMb/w4UeDOlXy29zO8uvVVh7Fh1CRu68cNAUtWV9IAx6rCPSFe3lh5CJsyS8U+Vay9YFQCrj2pj0iDZhimF4iYZcuWYdasWZg7dy7WrFmD+Ph43Hzzzbj++uZrEdfX14vFvmwxw7gk29fbLDHZUtfniKTkFs8pL6yFGTIo9bUInTQK8c8/140zBJZnLnfYT1QGYvzfbzS4+KyE9wfCXG8Vao6jRdV4/rf9+HN/odin8vsXj0vATdNTER/iue0MGMadcVsRc/ToUbz33nu4++678fDDD2PLli24/fbboVarcdVVVzU5/4UXXsBTTz3VI3Nlek+dGG+MiSnObtsSU3rC0l6gNg+a5JbFjiugQnYHyyT31fILliM+IF5kIcno/YgfCwTGSCdSfMzY69DdlNfq8Naqw/hiQ5bo2kz9gS4dn4ibp6cijsULw/ROEWMymUQA7/PPPy/2R40ahfT0dCxatKhZEfPQQw8JwWNviSF3FMM475fS+9xJeq0WZfl5Yjsy2TGLp7KkDssX7YG2Wg+dVuo95FeTD3XyuG6d48HSgyIWJswnTAgYma4GyN0iHbzwQyCsT4/FSH2/LVdYX8prpWwjinV5+MxBSI0K6JE5MUxvw21FTGxsLAYPHuwwNmjQICxdurTZ8zUajVgYxmXYGhF6tiWGHr4r339H1IIxGAyivopfcIhY7Dm6owjFOVJvJCuh5QehTr6oW+dr351alNSn2BeTAQhJ7jEBQ92hH/phD9YfkYJ2B0YH4tGzB+Hk/pE9Mh+G6a24rYihzKSDBw86Jh8cOoTkbjZlM0wTd5KHx8RQEbs9f1l6CllIGGzX+dlCmSWlesjJcUgbE4rjl86FX10RVIlJ6E6sKdVD/ROA4sPA/l+kA32nobsxGE34aF0m3vjzELR6k2gHcPfpA3DtlD5QKrzIRMcwHoLbipi77roLkydPFu6kiy++GJs3b8b7778vFobp6WJ3nkxRljWQNwXTrrgOcrkcsf0HNjmvLL9WrOMGhCDYWIyyuiIoIiKc3+SxEb8d/Q3ZVdm2/U15Uk2aYWveBFa82HAi1YPpRnJKa3HntzuxLatM7E/uF44XLhiG5HDX/jwYhvFAETNu3Dj8+OOPItbl6aefFunVb775JubPn9/TU2PQ2y0x8GiKLIG8cQPSkDJ8VIvnWUVMaIw/9OlSE0h1kmutMHuL9+KBtQ80GZebgaH1OkAdKAXwUgZS/5noLn7dfUK4j6q0BgRqlHjs7MGYOzaBO0YzTA/jtiKGOPvss8XCMG6Btdidh8fEWLORyBLTEnVVOmhrKFjVDO2Hb6Dm0IFuETHbC6VCdkmBSRgfO14a1GsxZtOnCDGZgDt3A35h6C5qdQY8uWwvvtuaK/ZHJYXg7UtGcZVdhnET3FrEMIxbYWkA6fHuJFtKdcsipixfiofxqStB9bff2MY1/VO7JYj3nH7n4MYRN0qDFANT818gclC3ChgK3r3+8604VFAtdOst01Nxx4z+UHHsC8O4DSxiGKYXtR2or61BVXFRm5YYqyvJv7YAviNGwP/kk6EIDEDwhRd1SxDvsAi7QOMsSyXe5EnoLtYeLsKtX+8QjRojAzXC+kIdphmGcS9YxDBMe7EVu/OsH5nJZLSlhxceOyrWgeGR8PFvqGViMplRVaK17Rcekype+9XmI/D8GQhfsMDl86yor0BOVY7YHlKQAZQclw5krJTWSZO7Rah+sPYoXlx+QLzdIxND8L8rxiA6yHu6ljOMN8EihmG82J20a+Vy/PXJIpiMUrE6K40L2/38xg6cOCw1eLSHRIwqbgq6A6sVJlmvR/Cy25qe4GJLjM5gwoNLd+OHHZJ4mjsmAc+cNxQ+KrtO2QzDuBUsYhjGiwN796/7u4mAkcnkSB3fIAj09UabgFGq5bZ/n7qqCOGl+6CKa0ZQOAmyvNyx+g5UaCtQZ6gTYyILKSAaiBrUcGLKSUBwgsvmUV1vwE1fbsPaw8WibcDjZw/GlZOSOfuIYdwcFjEM46Up1uQaKc6RgnjnPfUSIhKkQpEKpRIqnwb3SHmBFP/iE6DCda+eLL1Wr8eBESPJzwRVfLzL5rgicwUOlx12GJtaWweMuwOY/iC6g6Kqelzz6WakH6+En1qB/84fjekDo7rl3gzDdA0WMQzjpcXuqkqKUV9TA7lCgZh+A6BUqZo9r6xAykQKjfGzjekLCoSAkanVUISHuzwb6dqh12J2n9nw/+oSJNbUAglj0R1kFtfgqo83I7u0FuH+anx89TiMSHRsv8AwjPvSoe+UdXV1OH7cEmxnx969e505J4Zxa2+Sp2QnWevBhMUltChgiHJLJlJItJ2IOX5CrFWxsZDJ5S6zFFlFzPTE6UjTRCKxRKomjPgxcDUH86swd9F6IWCSwvyw9KbJLGAYxsNo96fTkiVL0L9/f5x11lkYPnw4Nm2SSoETV1xxhavmxzBumJ0k86h6MK2lUhNlFndSaHRD+Xz9CYuIcaErqaC2AMV1xVDIFEgLGdDQmTpiAOAbCleyP68Sl36wEcXVOgyJCxICJiWC2wcwjNe6k5599lls27YN0dHRYn3VVVfh4YcfxmWXXdZQP4NhvBkPcye1pzKvY3sBe0uMZHFVxce5PBupv04H3+ft7pMwDq5k74kKzP9wE8pr9RieEIwvrp2AYL+WLVUMw3iBiNHr9ULAEGPGjME///yD888/HxkZGRzBz/QKbGLdTTVMVWkxqoqLbfsFRw83m07dOFi5wmKJ0RQeQe0WyZ2jTZcEhirOeSKmqLbIJlyIP45JnbSH1kn3F8iVwODz4CrSj0sChorYUezL59eOR7AvCxiG8XoRExUVhd27dwtXEhEWFoaVK1cKiwyNM0xvqRPjjjExlcVF+PiO62E0GJoci0hsKmIorbq6TIvaCh0MehNkZiNKb7sW5WbrPxJOFTEkAK9ecbVDd2orIqV6xpPA6KsApQZQu8ats+9EJS77YCMqtQbRA+mza8cjyIcFDMP0ChHzxRdfQKl0PF2tVuObb77Brbfe6oq5MYx7YXZfd9Lxg/uEgFGqNfAPbYgnSRw8HIHhEQ7n6nVGfPnYBtRW6mxjfrWFkJtNUKekAJZAXmVkJAKmTXPK/PJr8oWAofiXIRFDbOMRuTswk7KRqBqvC/siUR+kKz/eLATMaIuACWQBwzC9R8QkJDgWmsrPz0dMTIzYnjKleyp6MkxPYja6rzupOEtyAw2ZdhpmLLi51XNLT9QIAUPxyWo/JWRmM+IOr4VMpULf5b+5xD1szUIaEDoAX535lTRYmQe8nkaqEIix65XkZAoqtbj8o00orq5HWkwgPrmGBQzDeAudzp2cOXOmc2fCMO6OG1ti2puJZF/cLjY1BAtem4pLLw9C4vE1UEZFuSy+rdnGjid2SOvINFJTLrlvea0OV360GTmldUgO98Pn13EMDMN4E50udscZSUyvw41jYqwiJrIDIibEko1kKMgXa6UlcN8VWC0xQ31jgFJLLZjMf6R13CiX3LNWZ8C1n27BwYIqRAVq8OV1ExAVyI0cGcab6LSIcdU3NoZx/+wk9/rdr6uuQnWJlJUUkSS1FmiNsnxLhV5LcTtRnZeCeGNcI2KMJiP2lkgFMYf930OA/l7HE1wgYqgr952Ld2J7drnIPvriuglIDHONtYdhmJ6D2w4wjIfXibHWgwmKjIbGz7/9lhiLiDEUFIq1Mso5IuZ49XG8vvV11BgksaQ36kVzRz+TCX30ekAd0HByQBQw8Ew4m5dWHMAf+wqgVsrx0VVjMTAm0On3YBim52ERwzDtxY0aQFJn6n3//IW6qkrkZxxqtxWG6sKUF9Y1EjEFTnUnfbP/G/yRJdWAsWeMth6KkGTgTteWZPhmczb+989Rsf3KRcMxNsV1WU8Mw3ioiFEoFM6dCcO4ObbC1G7gTjq4cR1+X/SWw1hUSt82X1dVqoVRb4JcIUNQuBQfoi+0uJOio5wa/3LJwEswLFIK5JVnrMLkTZ8DAybAlaw7XIzHfpKCiO+aMQDnjnRd2wSGYTxYxOzYYcksYJjeghu5k/IPHxTryOQ+Qryoff0wcuZZzZ6r0xpgsqSHF+VUiXVwpC/kCrmjO8kJlhiDyYD9pfvF9qWDLkXfYIuw2v2b6IqN2BFwFRmFVbjpq20wmMw4f1Q8bj8t1WX3YhjGPWB3EsN0uAGk+2QjjT7zXAydPqPF83atysG676X2A/ZYXUkUrNzgTpLqPnWFI+VHRPxLgCoAKUF2mVJ5FheSi+rBVGr1uOHzbajSGjAuJRQvXjiMkw8YphfgBt59hvGw7KQetsTQPNqbUn10Z1GTMblShtSxkuvIWF4Os06q3KuKinRaPZgh4UMgpyJ2hF4LFB90mYihTKR7vtuFo8U1iAv2waLLx0CjZHc3w/QGnG6JqaiowK5du7Bz507cfvvtzr48w/R4nZiedifVlJVCW1UJmVyO8ISkdmUiXXj/GEQlWzJ0ZDKYa2uhPXQIuqwsMaQIC4NMre70nAprC6E36bGlYIvYH1pwGPj4DOmgQQuYDIBvKBDk/BiV//6dgZWWTKRFV4xBeIDG6fdgGMbDRcyRI0fw6KOPQqPR4M0330RISAgyMzOFWLGKFlpnZ2eLb4r+/v4sYhgvzU7qWRFTZGkxEBaXAGUrwkNXZ7D1RwqL9bfFwJi0Whw54wwY7TpedyUe5v3d7+OdHe84jA3LPwjUSllQNpKnOD0o+u+DhXhtpZSd9ey5QzE8IcSp12cYxktEzPz588WSnJyMoUOHorq6GpWVlQgODsbgwYPFWE5ODj766COcdtppSExMdO3MGabHit15RouBMosVxi9IDbVvw5+6PjdXEjAymWSBkcsROu/iTs/nz6w/xVolV4kGj0lGEybUaYHxNwIpJ0knyRWSiHEi2SW1uGPxTpE1dtmEJFw8jj9zGKa30W4RU1hYKIRK3759RfPHBx54ADfffDPi4xvMwx9//DHGjx/PAobxTnooO2n7bz/j4IZ1tv2Kwvx2xcM0LmpnxVBcItbqPn3Q77f/69LctAYtDpVJlpDfLvgNMf4xwKsDpHz04RcDCWPhCuoNRtz89TZU1OkxKikET8wZ7JL7MAzjJSLm7bffxk033YSIiAgsWrQIb731Fvbu3YuXX34ZAwYMcO0sGcYd6IE6MUaDAf989YlYNyY+bXC7REyopUeSFYPFjaSMiOjy/A6UHoDRbESEbwSi/aKBmmKguqChsaOLeHH5AaQfr0Sonwr/nT+aA3kZppfSbhFz9tlni8XKNddcg/feew9Tp07FhRdeiCeeeMJVc2QYt4Cq3Qq60RJTeiJXCBi1ry/OuPku27h/SCjiBgxq9bVl+S1ZYoqcJmKs2UhDI4ZKKc0FUo8khKYAGrv2Ak6Egng/+Vdyqb128QjEBvu65D4Mw3hxijVV7L311luxb98+sZ2WlgaTyQSj0ejcGTKM27mTuj+Il4ra9R8/2ba0JWCI8oKaZkWMsURyJykiwp1WnXdYhCV12ipioofCFZwor8N9S3aJ7QUn9cGpaa7rvM0wTC9IsQ4LCxOupoULF+Kuu+4SQb33338/brnlFvj68jckxhuL3cl6RMS0BmUiLXlpKyosfZGs9VOadScVWd1JkR0ObF56eCmyq7JtYxvzNor10HCLaCm0iJgo58eoGIwm3LF4B8pr9RieEIz7z3Cdu4phmF5WJ4YylH7//Xf8+uuvuPfee/Haa68hLy/PWZdnGPfpnSR3PxFzIqPc5j6yJzTWH4Hhjl8mDBZLTEfdSbuLd+OpDU81GVdCjiHfXAkYdIBOsv4g2vki5u1Vh7HlWBkCNEq8c+koUReGYZjeTbtFTF1dnagRI5dLHxyrVq0StWEGDRqEM88803Yexc2cccYZeOcdx7oRDOPx9EB2UntFjDWIN2V4BKZdOtA27hekgrzRfBsCezvmTtpZuFOsqR/SSfGW1GlqfbBvOYKrpRgVgToASJoEZ7ItqxTvrs4Q289fMAzJ4f5OvT7DMF4uYiZOnIg1a9aIIncvvPACVq5cKcQLZSr99ddfePXVVxsuqlQK1xLDeBW2wF7X3kZbU43KokLU19agtqIcMpkcEYnJrb7GaoWJSAhAQGjrFWs7G9hrDeKd028OFgxb0HDg7/9J63lfApGDgIAowCcIzqKm3oC7vt0lfvwXjIrHOSPinHZthmF6iYgxGAxCwBBLly7F+vXroVarcffdd2PkyJGunCPDuJc7yYUxMXqtFh/fcQPqqiptYyGxcVBpfDpVE6YxZqMRxtIysa3ooIixBvFSJpKNmhKgxtKfqd+pgNr5FpJn/28/sktrER/iiyfPHeL06zMM47m0+ztldHS0EC5EXFwcii0m6aqqKtfNjmF6mTupIDNDCBi5QoGA0DAEhkdi9Oxz2nxdWX5Ns0G8jaGGj6AMQpkMyrCwds+rTFuG49XHbc0dbRTtl9YhyS4RMKv2F+CbzdlCN746dwSCfFROvwfDML3AEkPVeK+++mphfaG+SCNGjMDYsWNFJd9XXnnFtbNkGHfA1nZA5vIYmJSRY3D+/Y+36zXaGj3qqvTtssRY42EUoaGQKZUddiX1Ce6DQLWlkSRRaBExUW2nfHeUkup6PLBUsv5cN6UPJvXreko4wzDeRbs/xVJSUvD333+LujCHDx/G5ZdfLloOjBs3TtSJYZjeU+zOdfcoPCaJmKiUvu1+jdWV5B+igdqn9T/pjlTrPVZxDAv/XIiK+grRodqhHoyVogMuqc5L6dyP/JiO4up6DIgOwL2zGoKVGYZhOp1iTanUtFD/pJiYmI6+nGE8l25wJxVlHRXrqOS+7W5I2ZYryVBUBLNO6matO3Kk3ZlJK46tsLmQrExPnO54kossMb/szsOKvflQKWR4Y95I+Kj4ixLDME6sEzNz5kzs3r27sy9nGI/DbLJsuEjEUHuB4pwssR3ZhiWm5EQ1fnpth3AlWQltxpVU8uGHKHz1tSbj7QnqtbqQFo5YiDl958BP4YOIf14H1r3fcNLx7U4XMeRGenKZVDTv1lP6Y0hcsNOuzTCMd9FpEWP9FsgwvQYXx8SUUZ8kvR5qXz8ER0a1em7WnhIHAUO1YJKGNrWu1GzaLG0olZBZ3L5yHx8EzZzZ5t+3NRtpStwUJAUlAcfWAZvea3oyxchEOK8J7BPL9qK0Roe0mEDcNL2f067LMIz30WkRI5q9MUxvwsm9k3IP7MXGpYthNEhiRGvJ9KPCdjJLUcmWKMuTXEijz0jGyBmJUCjlzcbDGEqkGJjE/7yLgGnT2j23vJo8lGpLoZQpkRaW5ug6ih0JjL++4eS4UYDKOS1GVqTn49fdeVDIZXjlohFclZdhmO5pO8Aw3o6z68Rs/vE7ZO3e0WQ8Pq3tkv2lluJ2kYmB8A1Qt3ie0dInqbM1YfqH9oeP0lKjpuigtO47DRh1OZxNea0Oj/0subBunNoXwxLYjcQwTOuwiGGY9mLLTnKOiCnIlIJsp86/BkGRUjdmpVqNpGEj2nT1lFuDeWNbTqk2m0wwlJZK121DxFAAr7WtALEya2XTbKTigy7JRLLyzK/7UVRVj36R/rj9tP4uuQfDMN5Fp0UMp1UzvQ4nZidVl5XaWgqMnHVWmxV57amt0EGnNQqDUEikX9uF7egPvZXCdiSKrvv9uiaZSE2q81otMRHOT3dee7gIS7fnin/TyxeN4GwkhmFcK2J27GhqBmeY3hHY2/VLFR6TrDChcfEdEjBEqcUKExTpC4Wq5dgZQ5FdYTtVy5VuSbzQQvEvY2PG2sYjfCMwK2WWtFNXBlQXSNuRzgviJbR6Ix79SXIjXTUpBWOSQ516fYZhvBd2JzFMh4vddV3FFGYe7XBROyvllniY0JjWy/wbLUG9bdWEsaZSUwDvBzM/aP6kokPSOige0NhV7HUC/1mdgaySWkQHaXDPTOcKJIZhvJtOixjqnUStCDZs2CAK3xFU/G7y5MmiPUFkZKQz58kwPY/JeZl5VktMVJ+2U4gplbo4t9q2f/xQWbv6JNlaDLQRD9NsY0crtaWSFSZb6puGSOe6kjIKq7BojfSzeHLOEARybySGYVwtYrZs2YJZs2bBz88PM2bMwIAB0rengoICvP3223jxxRfx+++/i95KDOMt2GojKTonYqiQXcbmDeI6xw/sa7cl5sfXtqP0hORCsqctS4yhuESslRGR7bLEDIts1FIgZwvw8Uy7Kn/ODeq1thbQG804NS0KZwzlCuAMw3SDiLntttswd+5cLFq0qMm3UvpgWrhwoTiHrDQM4zVY3UmdNMT89u5rKDomuZHEZWTyNkUMWWGsAiYsrkG0+AWp0WdERPv6JIW37E4ymAzYV7KveUvM0b8lASNXSXVgyI005AI4iyXbcrEpsxQ+KjmeOmcI155iGKZ7RMyuXbvw6aefNvuhQ2N33XUXRo0a1ZlLM4z7u5M6EROj12pRnHVMbA87dSbkCgXi04bANzCoXUXtAkI1uPTxCR26p6G4SKyVkRFNvmi8svUV7C3eC51RB61RiwBVAFKCUhwvYE2pPvUR4KS74EyoIu/zv0nF8+6cMQCJYa27xhiGYZwmYij2ZfPmzUhLa960TMeio6W6Fwzjde6kTsTEFGZlwmw2wT80DDNvvL3dryu1iJiw2NZdR81htLmTHEXMkfIj+GLfFw5j42LGQd64FHGxJZjXiS0FrLzw236U1epFa4HrTurj9OszDNM76JSIuffee3HDDTdg27ZtOO2002yChWJiVq1ahQ8++ACvvvqqs+fKMD1LF7KTCjMzxDq6HYG89pTlWTKROiFibIG94RHNBvJSNtKNw28U4sU+tVpgMgHFh10iYjYdLcH323LF9nPnD4VK4aQ+DgzD9Do6JWJuueUWRERE4I033sB///tfGC0FtagA3pgxY4Sr6eKLL3b2XBmmZzF3vndSwVFrNlJqp2rCtMcSY9brYaqVRI9DTEwjd5I1kHdS7CTMSJ7R/MWqTgD6WikeJrSRm6kL6I0mW2uBS8cnYUxyy0X4GIZhXJZiPW/ePLHo9XqRbk2QsFG1UlSLYTwZW5JOJ9xJBZ22xDQN6m2pOu+ROXNsvZLsaexOajWlunF13rC+gMJ5f9Ofb8jCoYJqhPmr8cAZzq/8yzBM76LLxe5ItMTGxjpnNgzjBe4kk9GI/ev+Rl1lhdinV5XkZovt6L7tt8To6gyoLqtvV02YuvS9zQoYv3HjoLBrOaA1aHG47HDTvkiNsbmSnNfDiPoivblSirO5b9ZAhPi13LiSYRimxyr25uTk4IknnhDF8Bimt/VOOrhxHVb8940m437BIQgIC28zeHj1FwdQcKwSRoNk+vEPVkPj17o1xFAgtQTwnzIFiYvea7ieQoGsyizoTXqxn1GeAYPZgHCfcMT429VlqcwDjm9t2D/yl9OL27284gCq6g0YFh+Mi8cmOu26DMP0XlwiYkpLS/HZZ5+xiGG8MzupjZgYayG7yOQ+iEiS4klI9gycMrXNWigUyLt/fZ7DWGz/kDbnZiiURIwyJtqhT9K7O97B+7vfb3I+WWEc5vLZ2UCJ5PJywElBvTuyy2zBvE+dOwQKJ3UCZximd9MpEbNs2bJWjx892lDQi2F6W9uBgqOSK2b8eXORNnlqh25RckJqLxAe74+T5vYXVp/olNZryRD6fEnEqKIdq96uyVkj1oHqQKgoSBeARqHB3IFzG06itgJWAZNItWgs/77AaGDgmegqJpMZTy7bK7YvHJ2A0Unc4JFhmB4UMeedd574ILd9M20GZ/SXYRi3wmaJafl322jQoygrU2zH9O14PIm1Oi8Jl4S09mfuWN1JSrv6TLX6WuE+In445wdH95E9xRYBExgHXPcHnM3323KwK7cCARolHpjNwbwMwziPThVooEDeH374ASaTqdll+/btTpwiw/Q8QrBbNXsr+rw4OwtGvR4af38EN7KKtIeS45IlJiw+oEOv01vdSdFRtrEDpQdgNBsR6RuJaL9Wik/aito5L4jXSkWdHi+vkDKd7pzRH1GBPk6/B8MwvZdOiRiqBUOF7lqiLStNZ6CmknTdO++806nXZZh2YdcDsbXA3oKjllTqvv07ZY20WmLC20ipboyhoFCsVXaWGPtU6lbnUuL8TCQrb6w8hJIaHVKjAnDVZOfVm2EYhum0O+m+++5DTU3TrrpWUlNTsXr1aqf9hKlr9v/+9z8MHz7caddkmA5hL8rtRMzBDWuRsWWjbb8wUypqF9Ov44JArzOiorhObIfFtd8SY9bpYCyxtBiwEzG27tStpVITLqrMe6igCl9szBLbT84ZwpV5GYZxDxFz8sknt3rc398f06ZNgzOorq7G/PnzRSuDZ599tsXz6uvrxWKlsrLSKfdnGMJsrRFDWKwaRoMBK957Ewa73zsr8WmDO/yDE4XtzIBvoEp0qW4v+kKp0aNJqcBLB/9rm9+mvE1tF7WzFzHhHasm3J7+SEaTGbOGROOk/q133GYYhuk2EfP444/j3HPPFW4lV0MtDs466yzMmDGjVRHzwgsv4KmnnnL5fJheip2IsbqTinOyhIBR+/physXzbcf9QkLRZ2SjXkTNkHOgFJuXHYXRIF27vlbfYSuMfXp1UYARiw9963BMKVdiSMSQll9sNAClR51uifk3oxirDxZBKZfhwdmDnHZdhmGYLouY3NxczJ49G2q1GnPmzME555wjGkHSvjNZvHixCBImd1JbPPTQQ7j77rsdLDGJiVxQi3ES9pYYSyRZfoYUsBqTOgCjzzy3w5fc8Uc28o82tRjGpQZ36DrWzKTSAGBA6ACcmnSq7djIyJEIUgc1n5FUtB+oKQaoEJ7SFwiKhzOglOrn/m+/2L58YjL6RHS8eSXDMIzLRAxV4qUspH///Re//PKLCLbNy8vD6aefLiw0Z599NsLsSp13turvHXfcgZUrV8LHp+2MBo1GIxaGcQUOceoWd01ehpTVE5vaubTh4lwpE+nkef0RHCm1FVCq5IjpoIjRW0VMoAyz+8zGgmELWn+BtgL431RAbxfXFpEKyJ3TTfrHHcexL68SgRolbj/N+cHCjHtCjYCplx7DtAcyesid8JnT6Yq9dHOKjaHl5ZdfxoEDB0QRPArAveGGGzB+/Hhhobn00ksRH9/xb3iU/VRYWIjRo0c7/JH8888/ePfdd0X8C3XNZpieciflW0VM/467YWordair1Il07UGT46DSdOx3uf7IEWRdcSWMpaW2sdJAYFpbQbxE4X5JwJD1JXY4IFMAk2+FM9DqjXj1D8lCdcupqaLRI+PdUCZqfn4+ysvLe3oqjAdBGqJPnz5d9uB0ue1Adna2EBZkBbngggtw//33o6ioSAgaa2Xfe++9t8PXJffUnj1SiqiVa665BmlpaXjggQdYwDDdi13LgYrCAhh0OpQczxFDMf06LmKKc6vEOiTKr8MChqheu9ZBwBjkQHqyHDeFtyOg2BrImzQRuPInOJOP1mUir0KL+BBfXM0p1b0Cq4CJioqCn58fFzpl2oQ8OSdOnBAenKSkpC79znRJxLz99tsiDsX6i0uZRGPHjhWZRNddd51YOktgYCCGDh3aJOspPDy8yTjDuBqzqcEa+OFtDb/XQZFR8A8J7bQrKbyDRe2s6E+cEOvQ+fNx6NwRuGfjQ4iL6ifaC7R9c9cUtyuprsd7fx+xdan2UbGl1NuhvwergKHPZoZpL5GRkULIGAwGqOz6vXWULjmknnnmGTz44IPil7iiogIHDx7ElClTMGnSJKxbt64rl2YYt3QnmUnNyGRQajRQ+/pi5MyzOnW5EouIiUjonIgpPyZZU14tXoxbdjwMrUbWdiq17eaWNgPhzhUxb606jOp6A4bGB+GcEXFOvTbjnlhjYOiLLMN0BKsbiYRwV+iSJYYsL1dffbUtOIeK3L3++usiqPeee+7Bpk1SnQpn8ffffzv1egzTYREDEwZMPAlz7nygU8XsirKrhGuqMKuqSyKmMjsDZHMpDKJ5ySCXyXF68unte7GtuJ3z6sIcKarG15uyxfbDZw6CnLtU9yq4Vx7TU78zXRIxVEF3w4YNQrzYc/HFF+O5557r6twYxm2wttGgdVz/tE5dY8X/0pG9V6qsayW8kyJGVSgFUZ4x8Qq8fuoC+Ch8EKBux7WMeqAs0+mWmJeWH4DBZMZpaVGY3I8L2zEM0z10ScS89tprIpiXzEIkXKzKiiww/ftzaiXjPZiNDZaYuIEdFzEmownHD5WJ7eAoX2GpiB8QioDQjpcFMFRVwafWILaHDj0FEb4dEA1lWYDJ4NS6MJszS/HHvgIo5DI8dGbnBB7DMEy3i5iTTjoJn376KRYuXIjbbrsNI0eOhE6nQ3p6Or744ouuXJph3IrKQqnBohlmRKX07fDry/JrYdSboPZRYP6TE1ttItkWJ47sFOsqX2BEwsjWTzYZgfw9gFEn7eduaWgx4IQaDWSZeu43qbDdvHGJSI1qR2Axw7gBFArx2WefiWrvFNtp5aeffsL555/v9CbGjGvocor1mWeeicOHD2PVqlUiZmXHjh1i3FrwbtiwYcLt9OabbzpjvgzTLVBjx78//1D0RyICEYJTwi+BTKGAQtnxSHpbDExiYJcEDJF5YBOiKCYtzBcaRRuWnL+eAda90XTcSfEwv+7Ow66ccvirFbhzBltfGc+CCqm+9NJLuPHGGxEa2vEsQ6bncUqJTqoRQ2KGit5Rhd3i4mIcO3ZMVPadPn06srKkTrYM4yls/20ZqktLUFdZIRa9VivGlerOpQKKgF5KK0zuuKUisyIT3x/63rbs37tWjJti2uFGOvavtA6IAUL7SEvUYGDMNegq9QYjXv79gNi+cVo/RAW2XVmbYdwJ6skXExMjrDEtsXTpUgwZMkQ851JSUkQYhT009vzzz+Paa68VpUGo7sn777/fpAI9hVyEhISIL/dU2Z6ekYwbWGJagvoW0UK9lRjGk9Dr6pF/RMrgufCRZxAQEgpTQT10Swuh7mQqqU3EJHZMxJBJ+/o/rkdBrdRagLg8W0pJ9EtIavsCJZZMpMuXADHtqObbAb7YkIWc0jpEBWqw4OQ+Tr0247nQ72ydvmtps53FV6XoUNYLVX0nAXLZZZfh9ttvR0JCQpPK8SQ+nnzyScybNw/r16/HzTffLGrikDvKCgkbKjny8MMPY8mSJbjpppswbdo0DBw4UKShz5o1S5QeWbt2LZRKpWhmfMYZZ2D37t1O7znY23CZiGEYT6Ug4zBMRgP8Q8OQPGyk+FCsN1eiCIVAO11BFUV10NZY+siYGyr0RnXQEpNbnYs+2/Jw6SEg0o+cSEB0JmUmadE3bWLrL64pAeqkYGKE9YMzKa/V4Z2/pHoz984cCD81f5QwEiRgBj/+e4/8OPY9PavDv4sU/0LxnE888QQ++ugjh2NUMoSqxz/22GNif8CAAdi3bx9eeeUVBxFDnggSNwRVlH/jjTewevVqIWK+/fZbUaH2ww8/tAmsTz75RFhlKARj5syZTviX9174k4dhGnH84D6xjh84uOFbnaVOTHviWXIPluHnN3Y0/WPTKBAc1TFLTnrRHiz8zQQ/EZeb73AsYMCg9llhghMBtXOLkb37VwYq6vRIiwnEhWMcv70yjKdBcTGnnnpqkxY5+/fvF64fe6igK8V4UpE2a/8+ivu0Qp8Z5KKi3n/Erl27kJGRIVxN9mi1Whw5IlW4ZjoPixiGsUD9kCj7KHd/utiPTxti+9mYrQ0g22GIyUqXasFQJpLaT/oTk0GGwSfFdrgI3KGMzUjRAWb6YHzwAVsHbWVkJPwnT25fUTvKRHIi2SW1+HyDFOf20JmDRGo1w9i7dMgi0lP37gxTp04VLp+HHnrIwcLSXhqXzSchQ9YXa1HYMWPG4Kuvvmq29D7TNVjEMAyAVR8vws7ff3X4WcSnDW6mAWTbD+yCzAqxPuni/qJDdVcoPLRLrA3RoQi76qqOvdhqiXFyjyQK5tUZTTi5fwSmDeAPYQZNHuCe6F588cUXhVuJXEBWBg0ahH//tQTHW6B9citZrTBtMXr0aOFSov5SQUFBTp93b8cp2UkM48kYDXqk/73SYSy6byoik1MaBiwNINtyJ1FRO2sQb3RKcJfmpTfpoT0mVdfV2M+lvRRbeiRFdLzLdkvsyC4TadVkEHpodhvuLIbxIKgcyPz580VjYyvUPofKh1DQ7qFDh0RdmXfffbeJ26k16JoRERHCLUWBvZmZmSIWhgKJc3NzXfSv6T14nlxmGCeTn3EYhvp6+AYGYcE7HwqXjUrj45DlYCt81YaIKc2rgUEnFbULjel4HArd5+mNTyO9OB0GkwHjSqQidUF92yFE6quB2uKG/aIDTnUn0dyetxS2u3B0AgbH8bdKxrt4+umnhdXE3ory3Xff4fHHHxdCJjY2VpzTEZcTNcf8559/RMAvVbivqqpCfHy8CBhmy0zXYRHD9Hpy9u0RP4PEwcOg9m1BeLQzJqYgs1Kso1KCOlXU7nD5YSw5tMS2P8eSXKROTm79hdWFwDtjgHrp/g44yZ1ErQW2HCuDj0qOe2Y6z7rDMD0BVZtvDNV8qa+vdxi78MILxdISzdV72blTqqpthQJ9yYrDOB8WMUyvJ2fvbvEzSBjSSh2VFtxJ+9adwJEdUhYCUV5QK9bRKZ2zUpAFhhgSPgS3j7odQUueolBaqJPaEDE5myUBI5NLfZGs9JnqlB5JeqMJLy6XLDsLTuqL2GC7ezAMw/QQLGKYXhf/smP5L6irluJWiBOHDtgsMS1idSfZuZgMeiP+WXwIRoNF4dgRP7BzJcz3FEtWoYmxEzEpbhIOnSgR+kmd3EZhu+JD0nroRcCFH8DZfLM5G5nFNYgIUGPhdOfWnGEYhuksLGKYXsXeNauw5suPm4z7h4QivJUKuLYUa7ljPyQSML6BKky+sCHuxC9IjYS00C5ZYoZFDIOxtBSmmhopRqdRJdEmlFiDeJ3fv6hSq8ebf0qZTnfMGIAADX9sMAzjHvCnEdOryNyxTawThwxHZJIl40cmQ/9xk1ovV24tdmd3Tl4GVc4F4lJDkDYxtstzqzPU4XDpISxYYUTij28jy1K6XRkbA7lG0z5LjAtEzKK/j6C0Rod+kf64ZFyi06/PMAzTWVjEML0Gk9Foi385+bKrEJvaUA+iLazeJCjsRMwRqR5MbGpIp+ZjNBlxsOygSKUmjpYfRVSJETN3mGHEIVi7z/gOH9H25GyF7ZwrYk6U1+GjdVKa94OzB0Gl4KoMDMO4DyximF4DNXWsr62Bj3+AqAPTIWzZSTKbeynfJmI6Vw/mrR1v4ZP0TxzGxhVL91H37Yvohx4E5HL4jRrV+oVqigEtWYVkQLhz41Ve/eMg6g0mTOgThhmDpN5NDMMw7gKLGMar0ddrcXjzBlEHJmuPlPaYNHQE5HJFp0RMaX4N1ry0FSajGfW1BtEPKSIhoFNzW3d8nVhH+kZCrZA62Q6vriHfEHyHDUXAySd3rDJvSBKgcl7WUPrxCvy447jYfuSsQR3qDswwDNMdsIhhvJpNP34nFnuSh7dh2WgGa7G7yhItCmqsjh4gMS0U8k64WGr1tThSLjV/+/bsbxHpJ5XvP77xflTiF6j7tsOiQh2qtZVSerWT42Ho3/vC8v3CU3XOiDgMT+icy4xhGMaVsIhhvJrDm9bbmjn6BgbCLzgEg06e3vELWbKo6aEe0zcYo2cliZoxnY2H2V+6HyazCdF+0TYBQ9QfkbKMNKltiJjcrcBHMwFzg6ByZnuBvw8V4d+MEqgVctw3q/2xQwzDMN0JixjGaynPz0PpiVzIFQqcd/9jIham01gsMfTffqMj0WdE1xof2qdS225hNEJ35KjY1vRrQ8QcWS0JGJkCUGoATSAw+Dw4A4PRhBcs7QWunpKCxLCOt09gGIbpDjjVgPFaju7YItbxAwd3TcCQIcZS0I5ETGdrwDRX1G5oxFDbmP7ECZjr6yFTqdquC2NNqT71EeCRPODeQ0DSBDiDJdtycaigGiF+Ktwy3Tl9lxjGW6Au1tQsUqVS4bzzGr44PPbYY7jhhhts+5dccglee+21Hppl74EtMYxXQDEcG39YjIKjGQ3F6DIlq0bf0eM6fD1tjR6Zu4pgNFjqw+wrAeUgyZVyhMd1ThBtztuMX47+ItxIG05sgMJoxoh8H9TUSjEt2nRJ2Kj79IFMqWxnXRjnunpq6g14baV07dtO7Y9gP5VTr88w7gI1cWyun9GsWbOwYsWKFl939913Y+TIkVi+fDkCAqTPgvz8fLz11lvYs0f6GyYeffRRTJ06FQsWLEBwcNc62jMtwyKG8QoKM49g/XdfNRmXyeToN7bjFop13x/GwY35tv1UjRzBvgr4Bqk61diReGrDU8iuyrbtX/2PDP4bn0XDiIS6X9/WL2QyNdSFcWIcDPHB2qMoqqpHcrgfrpjYRr8mhvFwzjjjDHzyiWOZA00LhSX1er2wvhw5cgQLFy5Egp219MMPP8TkyZORbNeodejQoejXrx++/PJL3HLLLS78V/RuWMQwXsGhTf+KdXzaYAyeepptPCwuHqGxHWuAaDKacGx3sdhOHBQKpVqBsAotUF6PiMTATs2vTFtmEzB3jL4DCpkC45d8CyALqrg4yPyk1Gi5xgdhl13W+sWqTgD6GjILAWF94CwKK7X43xrJenX/rDSolextZjoZP6aXGqF2Oyo/h/5mbUGChTpMNweVFPjvf/8rLC6rVq3C3LlzbZaba6+9ViwkgMiis3jxYtx0001NrjFnzhxxjEWM6/B6EVP43k7U+XQtHoJxf8KPh2Bm3NUI0kTBJ93u/U7XoeCP7R26lr7eiEkKQBakRBjMkOkMMBpNIkHJJ0Cq59LZQN6UoBQsGLYAZp0OB4+9IWJskj7/DOq2YmCacyWF9QUUznP3vPHnIdTpjRiVFIIzhzX/wc4wbUIC5vm4nvlBPXwCUPs77XJPPvkkXnzxRbz55ptQKBRie+DAgXj66acxb9484SYqLS3Fvn37MHbs2CavHz9+PJ577jnU19e3aOFhuobXixh9fi30Gi7S5e0EycMA+oyoBPSVVDCuawRb2gsY8hyvpQjTdEnEWAN56zMzYdbrIQ8IgCq+Y5YiV7iSDhVU4dstOWL7US5sx/QSfv31V1tci5WHH35YLMRll12Ga665pomFhsSL1YJz8OBBEZMXF9dUuNGYTqcTMTP2ribGeXi9iAm/YhCCAoJ6ehqME9n447fI3S+JAnti+w/ElLmXd+qaVIn3xOFyYQqnnkh1VXqMmpmExLQw2zkytRzqpCCnZCNp90spzD5paW1Xwq0uBI5LjSsFR/92uoihlGoqSjx7aAzGJDf8mxmmUy4dsoj01L07wCmnnIL33nvPYSwsrOH3vznrSmPq6urE2sfHp8kxX1/JTVxb20PutV6A14sYn9RQ+ASxiPEWtNXV2Ll7BYwGA/qNnQiFJYtHqVJh5Hnnwieh4+nP9C1q1ef7UVmstY1R8G789ET4BHXOfWQwGcQirg8z9pbsdagLU3/goFhrBg1q+2KfnNnQWsAeJ4mYfzOKsfpgEZRyGe4/I80p12R6MSTKnejScSX+/v5ITU1t9XhbREREiHVZWRkiIx3rR5GriWg8zjgPrxcxjHdxcMNaIWAik1Jw3n2POuWaxTnVQsAo1XIMPzVRjMX2DYZfJwXMwdKDuGrFVaih4Fs7lHIlBoZJKdHaAwdslphWqS1tEDDxdt8KA2OAtDPRVUwmM577P8kqdPnEZPSJ8IyHD8O4C5SBFBQUJOJiBgxw/GKRnp4uspisQodxPixiGLcnZ98eHFi3Rlg0stN3ibFBU0912vWPbC8U6+Sh4Zh0Xte7QP+Z/WcTAUPcljUQx8+/mFpgo/5YlhjzGdSGiCmSLDYITgSuXwVnQw0e9+VVIlCjxO2nOa/3EsN4AhRwS/Eq9iiVyg6JDrlcjhkzZmDdunUOxe+ItWvXYubMmU6bL9MUFjGMW2M06PHrmy+htqLcNkZtBNKmTO30NUtP1EBbq7ftH9lRJNZ9RznH5GuNf7lv7H24aMBFNpfV8VNmob6kxHaeIjwc6lZM2YJii4iJdH7/Iq3eiFf/kK5/8ympCPPvnOWJYTwVKmoXGxvrMEbZRwcsltL2QgXtrr/+erz88stC1BBarRY//fRTq4XzmK7DIoZx+waOJGCoceOoM+aIsdjUgQgM65x59tieYvzff3Y3GZcrZUgZ2nWTL4kVaybSmOgx8LMEGupyj8NIAkalQtL/FpESg6ZfX8jV6vZZYpxcmZf4aF0m8iq0iA/xxTVTUpx+fYZxZz799FOxtNW5vjHl5Q1fqOyL5lEm0rfffotLL71UjFENGUqxnjhxohNnzTSGRQzjVhRkHkF9TbVtf/tvy8R6xOlnYuIF87p8/b1rpawJ3yA1NL4Nv/4DJ8ZAbbffWXKrclFRXwGVXIUBoQ3+ce1uyQ3mM3Ag/CdPbv8Fi1xjiSmprsd7fx8R2/fOGgAflcKp12eY3gRlGL7//vsObQeouu8777zTo/PqDbCIYdyGgxvW4dc3X2wyLpPLMfy0WV2+fl21DtnpkjvnvDtHISzO+UGsVlfSoLBBUNkVoqvbJVl/fIcP79gFXSRi3vzzMKrrDRgaH4RzR3SwTg3DME2gfkq02LuYGNfDIoZxC8wmEzYs+VpsB0VGQe3ja0vXHHTSdASEhXfqunqdEUd3FMGgM6Iwq0pk40QkBjhNwFAm0sa8jbb9dcfXNelOTdTttoiYER0QMfVVQGWu02vCZBRW4evNUguER84cDHkne0ExDMP0NCximB5Dr9Xi4MZ1MOp1KC/IR0luNjR+/rjy5XfE2hlsWnYUu/6UKtFaGTjBOSX1yWe+8M+FKK6T+ixZkZnNmJAhQ3nuUuuJ0O7bJzZ92rLEkB9+1VNScTu9VEQL/lGAn/MK0D3/2wEYTWacPjgak/p1ThwyDMO4AyximB5j1cfvYe8ax7ThUbPnOE3AUA+k/f/mie2ENKmRI9V+GXySc/q6UENHEjAU/zIrpcHdNXRXJWLf/QzSnRtQBAdDndJGAG3pUWDdG45jCePgLKiw3V8HCkVhu4dmc2E7hmE8GxYxTLe5i3Tahoq4hZkZNgHTb+wEyGRy+AUFY9ycC5x2z8NbC6CrMyAowgfn3D5SVOF1RfzL4PDBeOHkF2zjJ5Y/ggoA6n79oE5MtLnFguec3XaLgULJYoPwVGD6QyKLCX2mOWW+ZH151q6wXd9IbozKMIxnwyKGcTk15WVY8txjKM4+1uTY0FNOx6yFd3T5HhT78sMr20T13cYMOTne6QKGsKZSW1sJWKndvFmsox+4HwFTO1jPpvBAg/VlmFRjxlks3ZaL/XmVCPJR4g4ubMcwjBfAIoZxCVUlxWIh/v3282YFTFBkNE6+9Cqn3O/A+rxmBQy5jwZNcSxm5SwaN3Uk9CdOQJ+bCygU8B09uuMXLZIsJYh0rqunpt5gK2x326n9EcqF7RiG8QJYxDBOZ9fK5Vj10Xswm022MZXGB5c++ypCYxvSeRUKhUif7gyUZSSCYC3bO1ZK2TaTL0gVNV+saPyUUCg7d4/W0Bv1OFByoIklpnbLFrH2GTIEioBOuGsKLSImqh2NITvA+/8cRWFVPZLC/HDl5GSnXpthGKanYBHDtDsThyrnUhsAK7raWmxetlTqZ2QRFNbziIDwCCFUVD6+wuJCTRudwc4/s7H+hyMwk5CxwydAhaHT46FSu6Zwm8lswkNrHxJuJKPZCJ1Jh4u3qGD89V5kWs4xWPqw+I2za9bYXox6oPiw00VMfoUW//tHKmz3wBlp0Ci5sB3DtAbFrv34449NeiH1JFOnTsXChQtx2WWXQafTiWaTS5Yswdixnfis8SK8XsR8dv+t8G2rtDvTJvU1NdBWV7X7JzX+vLk46ZIr2w5k7SC1lTps+iWziYAhxp2V4jIBQxwqO4TfMn+z7ftpzbjgLy20poYqnVYCp0/v+A0oM8mkB9QBUsNHJ0FuJK3ehDHJoThzmHPSyxnG07n66qtFCwHqb9SYvLw8hIaGwl1YtmwZCgoKcMkll4h9tVqNe++9Fw888ABWrXJ+Y1hPwutFTGVhAXSqhsqpTBeQyaBUOv4s49IGY8J5F8M3MNA2pvEPQFCEc5opEpXFdcjZXyq2s/eWwlBvRFRyIObcNhKwaCS5Qga1j2t/na0xMMMjhuPecfdC8fdmyE1vQpWUhOiHH7KdpwyPgO8wx2J3zVKRC+z6BjDoGkSMtTqvk8Rf+vEKLN0uFcx79KxBTheVDOONxMT0vNgnq7bRaBRdtd9++21cc801tuaSxPz583HPPfdg7969GDJkCHorXi9iLnr0OQR2JjaBcUCpViMsLkGsuxPKOvrpjR2oKmlIzyYmnd9PuI+6E2s20vjY8RgVNQondi4RqdSBp5zSOcvLn08Be75rOh7dDgHUzg/B5/5vv/D0zRkRh1FJ7vPNkvFe6PeuzmAp1NjN+Cp9nSLU7d1Jx44dQ58+fbB06VLRC2nTpk3o378/Fi1ahEmTJtles27dOjz00EPYunUrIiIicP755+OFF16Av79U9+qLL77AW2+9hYMHD4qxU089FW+++SaioqLE8b///hunnHIKfvvtNzz66KOiD9Mff/whBMpff/0lXmsPWYqmTJmCxYsX45lnnkFvxetFTGzqAAQFBfX0NJh2oq3RQ6c12Pb3rM4VAsY3UIWYvsFijNYJac6rYNuZbCT6oK5e+4/Y9596cucueGKHtB50DhBo+ean9AHG3+CU+f65vxAbjpZArZTj/lnO74LNMM1BAmbC1xN65Iez6bJNts7xzuaRRx7Bq6++KgQMbVO36oyMDGEpOXLkiOhk/eyzz+Ljjz9GUVERbr31VrFQN2tCr9cLsTFw4EAUFhbi7rvvFi4tEi32PPjgg+I+ffv2FUJlzZo18PPzw6BBTePkxo8fj7Vr16I34/UihvEcMncXY/l7u60xwg5Mn5+GviOd56LqKLX6WhwrycCj3xoR/+6DOAQZTFVVkPn6wm9cJyrqUkuBUinYFme+0iBinES9wYhn/08qnHfdSX2QGOaaD3aG6S1QDMpZZ50ltp966ilhISERk5aWJiwu5N658847xXESOuQCmjZtGt577z34+Pjg2muvtV2LBAodHzduHKqrqxFg5y14+umncfrpp9v2s7KyEB0d7eBKshIXFyeO92ZYxDA9CgXq6usN0GmN+OtzyfUhV8ocTML9x0Wjz4iIbp0XpVDvKtoFPQXaAjhacRSDM40YfswMM6ph1VlBs2dD3hkXW9EBKmMM+IYBAdHOnTyAj9ZlIqukFlGBGtxySqrTr88wrbl0yCLSU/d2FcPt+p7Fxkq1p8iiQiJm165d2L17N7766ivbOWStNZlMyMzMFFaUbdu24cknnxTnlpWViWNEdnY2Bg8ebHtd42yjuro6IYKaw9fXF7W1tejNsIhheoytvx3Dpl+OwqYIqNp+QgDmPjAWCpXza7t0hDe2v4Ev9n3hMLbwgDTR4PPOQ/iNN0BG6eMJCZ27QYGlvUD0EKcF8douXanFu39liO0HZ6chQMN/5kz3QV9AXOXS6UlUdgki1i9ZViFC1pQbb7wRt99+e5PXJSUloaamBrNmzRILCZ3IyEghXmif0qXtscbQWKH4GhI9zVFaWiqu1ZvhTzemW8g/WoGi7CoHCwyJGEKlUYgsI98AFU6/dnCPCxhiw4kNYp0YmAgfpQ/kRjMmZlANF6MQMZo+fbp2g0I7EeNkXlpxALU6I0YmhuC8kQ3FBRmGcQ2jR4/Gvn37kJravNWTgnRLSkrw4osvItHST40CgNvDqFGjkJ+fL4RM47Tv9PR0cbw3wyKGcTkHNuRh1WeWSrSNGHFaIk6a29+t3oUafQ0U6Ydx/wYTJkVFQi1Xw1RXh7paIxRhYZ0rZNeYgnSXiJgd2WX4Yftxsf3kOUMgd0HPKIbxFioqKrBz506HsfDw8A5fh+q1TJw4UQTyLliwQFhTSNSsXLkS7777rrDGUG0Xym6ignUkPtqbUUQihawx//77L84++2yHY2vXru3VmUkEixjG6RTlVIn4FmtadH2dlG0U1z9EZBlZiUgIxOhZSW73Duwr2YerVhqRmg/oMzajoUYxEHTmmcKN1CkO/QHs+FyqbpwjtSdAlPNEDLVfePIXycJz0ZgEYYlhGKZlKK25sSXjuuuu61S8DGURUdbSySefLOJh+vXrh3nz5onj5PL59NNP8fDDD4uAXrLcUAbSOeec0+a1qeo51YghN5S9iNmwYYMQYRdd5NxGsZ6GzEw/bS+ksrISwcHB4k3mFGvXcWRHIY7uLBIxqlay0kugswgXK4NPisP0ywa6pJu0s/nmt5cw8u5PYZLLEP/Ms5BZyvTLNBrRlVru10l//1sjgTJrgwLyo/kD9x0G1I4+8K50qb7n+13wVyuw+t7piApqPhiQYZyFVqsVgatUR6Wl4FOm65A7ibKhtm/fjuRkqfcZCaQRI0YIYeRtvzsdeX6zJYZp3y9cjR571x5HfU2DOKku0+Lw1sJmz49NDca0SweKSrpKtQKBYe7/AWesrhZWEuUvq8V+2dhUDLnwAudcvK68QcDMfgVQKIG40U4TMNX1Bry4QmpIedtp/VnAMIyXVRD+6KOPRDAwiRgKBh42bBjuuusu9HZYxDDNUpxbLWJZjHqTSB46tqsINRWOUfRWhk2LR3BUg3VC469E6ugoIV48hR13XAef39eLbWu9XL/zHP3PXaJgr7SmnkgTnFPMzh7KRiqqqkdKuB+umeKcRpsMw7gP9s0oKb6GqvoyLGJ6JUaDCQadUWxXlWqx888cVBbVOcRWFB6rbFJ0LjjKF32G29VrkcmQNCQMiT1QPdeZ1Gdm2gSMlaNxckw7Q/JnO4X83dI6ZhicTUZhNT5aJ/VdevSswdylmmGYXgNbYrwUEiLHD5VJac12YqS8sBaHNxfAoLcLYmkBKjAXkSBVkvQNVCNtcqxLu0T3FHmffijW21JlUDz/EBQyOQbHDEeAn9TmwCnkWzpdxzQUzHIGFNL2+M/p0BvNOC0tCqcNkvqwMAzD9AZYxLg5FCBL8Sj2VBTWoeBYBYzGBnVi0JmQs68UlSWSRcVsNLdLqFCgbd+REUgdQ2WtG4JuA8N9EJnU0JnaW9AdO4bcO+6EsbLSNqYvLBDNsLdNj8Mbw69wzY1dZIn5ZXce1h8pgUYpFynV3KWaYZjehNuKGOpF8cMPP+DAgQOitPLkyZPx0ksvieZZ7gx9M6bgV7KE2I9RurFNjJglN051eb3Da/V1BhQfrxaChCCXT1lBrYMlpSNo/JRIHBwGpV3xOIpToTL+0X2kiG96eMsVPV9crrsofP0N1B886DBGP4PDsYDfhPHOuYlIod4E1JY07BcecLqIqdTq8cyvUkr1raekcn8khmF6HW4rYijn/pZbbhENsgwGg0gjmzlzpigg1Lgsc2sc2V4IP58aISooDVham4WwsO5TjAgFsNJC++Icoxkms2VtMgthoa3WwWiwVxRm1Nc6Wkp0dUZxPWciRIhdZrKPv0rUXFH7Nrx99A08KiUQMX2CIbNokoAQH7eoftvTGIqLYa6vhy4nF1V//CFieRL+8x8oLeW6X9jyIv7PtBP3RjpJYBz6HfimmXgaTTAQ4ry6OG+sPCSCeftE+OOGaX2ddl2GYRhPwW1FzIoVKxz2qVBQVFSUaKI1derUJufX19eLxT7PnFj1+X74OimNtdPISFBoRFyJtU2OX5AagRG+kNv1zaHGh+Fx/vAJlBoK0rGweH/4B2t6auYeT+nnX6Dg+ecdxvxmnQ795BGiiB2J2dV7s6Crl2FYhJNEzLG10jooXloIep9HXua0Pkl7T1Tgs/VS24anzx3CwbwMw/RK3FbENIaK3hBhYWEtup+oPXpj4lKDEeAfCJlcDupkTjEgYpHJxD49VMjSQYtcJYdCIR2n+BC5bVsuBAZVm1VaCp9ZUfkq4Eeiw/JsosBX/1ANFL3IReOuGIqKUPTmm7ZCdfRem0ICcVOff5Dz7TSHc1VyFQaEDnDOjU/skNanPioJFydDlsHHfkoHeSzPGh6Lk/v37gZwDMN0jpSUFNx5551iaQ9kTKBzy8vLHcbJsEDtFC677DJRw2bAgAFYsmRJk47cvVbEUKdQ+sFNmTIFQ4daq3g48tBDD+Huu+92sMRQo62zbx3JFXt7AWaTCYWvvOoQ76LPz4epthY+w4cjZfE3Qsi+s+Md5Ox+v8nrz009F2qFZAHrEiYjkLdL2o5zTWO277flYHt2uajM+9hZg11yD4bxdq6++mp89tlnTcaps3RjT4Cn82kL4mPLli0dCs9ojmXLlqGgoACXXHKJrYbNvffeK/pJrVq1Cq7GI0QMxcZQw6x169a1eI5GoxEL0zsp+/IrlH7ySdMDMhliHnlYCBgivVhqvPjohEcxL82JdWCslGQAumpA5QdEOMmyY0dZjQ4vLpeChO86fQBigt2/EjLDuCtnnHEGPmn0udGbniORlrjArkC9oKi3E3ksrMyfPx/33HMP9u7dK9oluBK393lQV9Bff/0Vq1evRkJCQk9Ph3EDtAcP4cQjj+D4PffalsLXXxfHwq69FnGvvGxbyALjO2KEOEbxL3uKpXotQyObt+g5zZVE9WDkzq+p89xv+1FWq8fA6EBcNZkr8zJMVyDBQiX97ZfQ0FDRGJIsCtQl2srLL78s4jLJ6kBMnz5dPJ9ooT4/1Gn6scceE58zVsrKynDllVeKa/r5+WH27Nk4fPiwg4UkJCQEv//+OwYNGoSAgAAhrPLy8hzm+eGHH4rj1GMoLS0N//3vf23Hjh07JsIjKJv3lFNOEfehnkrUIJKgfwuJDArJoPNoefLJJ23upDctLnfi9ddfF+0MyDpDnoybb74Z1dSOpQWKiorw119/Yc6cOQ7j9O8lz8nixYvRay0x9Itw22234ccffxRvAjWJYhh9QSGyr70WxhJL+rIdfhMnIuree2xWl8ZkV2WjSlcFtVyNASFOspLkpwOZ/zTsH/7dZa6kdYeLsWRbrogNfv6CYVBx3BXjhojMz7qGCuDdiczX1ym1kkigkPvliiuuwK5du3D06FEhUL7//ntER0fbziN3FHW93rx5M7Zu3YobbrgBSUlJuP76620uKxIt5HKhRobkYjnzzDNFlq1KpRLn1NbWio7WX3zxhbBmXH755cIdQ12rCVo//vjjePfdd0XH7R07dojrk9C46qqrbHOhDtp0nf79+4vtSy+9FBkZGaI8CQkVusZBi7udxFJz0P3JskLPW/o3k4i5//77HUSTPeQdIdFEAqsx48ePdxCBvU7EkAvp66+/xs8//4zAwEDRxZMgxUt1Y5jeQ8lHH6H8hx8pOArGigoYS0uh6d8fIXMbWtDLVCoEzZ7dooCxdyWlhadBpZA+QLqEyQR8eQFQLX0zcyB+NJxJnc6Ih3+UrEhXTkzGmORQp16fYZwFCZiDo8f0yA904PZtkHWgyzxZ+Rs/0KmcBy3PPvssVq5cKYQJhTOQYDjnnHMcziVrxRtvvCGEE9Uw27Nnj9gnkWEVL//++68QElZBQq/56aefMHfuXDGm1+uxaNEi9OvXT+yTZefpp5+23eOJJ57Aa6+9hgsukJrRksAgEfS///3PQcSQ8DnrrLPENiW5kBuHRAxZbui5SXMkS1Nr2Af4kpWGfgYUsNuSiMnKyhKizt6VZCUuLk4c77Ui5r333rMpYnvIf0nqlvE+RKHAAwdgKG6wstTt2I7i/0q/C1YUwcFIeOdtqFPadqdU1Ffgp4yfUGuoxaa8TWLMaanUpUckAaPQAIPszKkB0cAgxw+7rvLmn4eQXVqL2GAf3HdGmlOvzTC9FXK/WJ81VqwZsOROItExfPhw0TmaxEljJk6c6GD5mTRpkhAcRqMR+/fvh1KpxIQJE2zHw8PDhdihY1bIkmEVMERsbCwKCwvFdk1NDY4cOSKsPVbrDkG100iY2EPztL8GQdchEdNe/vzzT5HpS0VmKTmG7qPVaoW1iObZmLq6OuHiag4yNtDreq2IsfcrMt4Pvd+FL72M0k8/bfZ4+E0LETBlitjWDBwIRWD7WiJ8nP6xWOwZHuGk/kW5WxqsLhd9BFeRfrwCH6yVGjw+e95QBGjc9s+WYYRLhywiPXXvjkAumdTU1BaPr18vNYYtLS0VS1czeZrD6layQqLI+vyzxqN88MEHDmKIUCgalfuwu45VWFFmb3uh2Jqzzz4bN910E5577jkh5shdRAKK0qabEzEUB0RxP81BPy9nBA63BX8aMj1G+ZIlwspiqqsTKdImSy0gDX1zsPRxkkGGoLPOFAG7nfF17yzcKdZT4qcg3j8e4b7hOD35dOeKmATX1UIwGE14YOluURPm7OGxOG1Qgz+eYdwRETzaAZeOu0IWkLvuuksIiG+//Va4bshSYe862bRJsu5a2bhxo4hJIYFBcSJkyaBzrO6kkpISEZcyeHD7SiOQq4bcMhSfQhk/nUWtVgvrUGtQIVkSPWRJsv4bv/vuu1ZfQzE6FOpBQoaCee0hFxwddzUsYhiXY6yqQtFbb6M+I8M2Rm0A6nbsaPTbqETM448h9OKLnXJfg8mAfSVSb6EHxj2APsFODg63iZhxcBUfrM3E3hOVCPZV4Yk5rk1VZJjeBlV5t8ZbWiEXED2QKcCWasZQZg9lDFHWDj3g77vvPtu52dnZoj7ZjTfeiO3bt+Odd94R5xAkZs4991zhBqL4FYrtfPDBBxEfHy/G2wvFt9x+++3CfUTzoDlTEDEJh7vtaqO1BsW3kFWH6rZQ5hJZVRpbVsgiRfE59G+gbCOK5aFYndYgkULWGDqXrDj2UFDvM888A1fDIoZxOrqcHBQvWiQq5hL1hzNgaJQyaCXillsQNPsMsa0IC4OyhYrMneFI+RFojVoEqgKRHJTsnIsa9UBdGWDQAgV7XSpiDhVUif5IxKNnDUJkYO+pX8Ew3QEVtbPGj1ihmBWqPEtBqRT4S9A577//vsj4oR5+JAQISp+muBDKxCHryx133CECge1jOGmMHvDkkqHKtr/99lsTF1JrLFiwQAiOV155RQgocmmRoGpvlV2CLEEUoDtv3jxhDaJgYWuatRX6N1GKNTVapuKxNFeKj6F/Y0vQv5lEHsUO2YsYSu+mlO6LLmpIvnAVMrOXBp9QUBIpV/pBUmob03XMOh1MuoZml6aaapR+/DFqNmyUOjXbiRizVuvwWlV8PCJuvhkyn4YHMQXm+rqwENKSQ0vw1IanMCF2Aj6c+WHXL6ivA/4zASi3i7in3kh3S9YeZ6I3mnDBf9djz/EKnJoWhY+uGuuU1FGGcSYU9JmZmSkyZloK8PRWKOlk5MiRDnVWeiP5+fkiE4osURQATZBYIlFEWV6d+d3pyPObLTFMsxhKS1Hz778wW0RL3c6dqPjpJ5j1DSKmNfzGj0fweeeJirlyjRr+U6dC0UJtAldhTal2WjZSzmY7ASOTitmNcU2m3P/WHBECJshHiRcuGMYChmEYtyQmJgYfffSRcK2RiCGLE1mKKJ6oO2AR08swGwyo270bRkuXb2N5OapW/gmDvV+YUp0PH26XYPEZMgTh1y8Qac9W5P7+8BnWvQ9eo8mIu/6+y1aRlyivl/qEDI1wUnXeLClTAUMvcmk20v68Sry1Sqrq+dS5QxAd1Lu+4TIM41mcR19Y7YKIH3300W67N4sYD8RsNEKflw8YDQ1jOh3q0vfCYKkvYDYaoN2TjvqjRxxeaywrh8kiYNqCsoRUluqU8uAghM6bB59GDTjlbtJn5GDZQazOWd1kPEAVgNFRTio8l20RMclSpoEr0BlMuOe7XdAbzTh9cDTOGxnvsnsxDNN5qJI80/OwiHEBpvp6qZqrBbPRBGNxEYxV1h4UZhiKimEoKXZ4nblOC112Nsz1WpswqT9yFMYqR9FhLCoW3Zk7iyIkBKrERLEtUyjgN3EC/CgVzs5yooyKgk8HiiT1NFbXEQmWhyc0+GFj/GMQrHEsCtUpDDogZ4vLRcy7qzOwL68SIX4qPHf+UHYjMQzD9GYRU7ZkKQwUCW4ywkzFNmhtNAmRYTYZyQ8h3CYkGMSarBtGyzlGg1iLMYMBJm29CFi1uVnMZhhrqmGqaBAZJC5MNTUu/3dRmX2ZvRVEJoOmXz+o+/a11VhRJyWL5ocyVcPbTK/xGTgQMqV3vfVWETM2ZiwGhg3s+gX1WiB7g5SNRJQdAwx1gG8YEOGE6zfD1mOlePcvyY309LlDERXIbiSGYZjW8K4nWTMUvvACahtVNuypSpIibsRi7VCGhkJJ1QztCieRMFEnJUIeYKlGK5dDnZwMZWSEFEhqQREYIMQKWVEYCWssjNOCeP94FNjyQdPxpEkO75mzqNTqcee3O0VRu/NHxeOcEXFOvwfDMIy34fUixn/6dAT6+UkPfIUCMrklq0Qhh0wml9YqtcigAVk3lEq7cxWQKRXifFrLfHwh9/WRrBgWMSL39YWcmmtZHmwyjY8QHY0FhrO6qzJNqdHXiJowTgvipXTxg79J25GDAJXFIqL0Aabc4ZK34PGf0pFbVoeEUF8RzMswDMO0jdeLmITXXuU6MV6G1qDFlvwt0Jl0Yj+rMgtmmEX8S4QvWa26SEkGUHlcaux4w2pA5dqu6T/tOI6fdp4QXsC3LhmJIB8ndNhmGIbpBXi9iGG8jze2vYGvD3zdZNxprqSjlqyDpAkuFzA5pbV47Ccpnue2U/tjTLLzKhYzDMN4OyxiGI9jY95GsU4NSYW/Suoq66PwwdVDulh4TkcZX2bgyF/Sft/pcCVUlfeOxTtQVW/A6KQQ3HZqy910GYZhmKawiGE8iipdFTIrMsU2tRKgrtRO4de7ga2NCti5WMS8vOIAtmeXI1CjxJvzRkGpcH7AMMMw3c+xY8dEOf0dO3aI1gStccUVV4iO19YS/RMnThQ9ki688MJumq1nw5+ajEdBXakp/iU+IN55AkZXA+z40nEsbhQQ2/qHT1f4fW++6FBNvDJ3OJLCHTvKMgzjOijJorWlcXNEV7Fr1y7REJK6VFuharfU7dpkV2uMaRm2xDAeWQ/Gaa0EiKNrAGM9EJIM3LReyjxT+TkU/3Mm2SW1uPf7XWL7upP64Iyhjl10GcaToZ7CBl3PPICVaso6bfvvNi8vz7b97bff4vHHH8fBgwdtYwHd1OftnXfewdy5cx3uN3v2bNG5evny5TjrrLO6ZR6eDIsYxqNwelNH4tAKaT3gDEDj2g8vrd6Im7/ehiqtFAfz4GzPqYrMMO2BBMz7d6zpkR/WDW9Ng0qjaFfTQivULZmEj3XsyJEjuPHGG7Fx40bU1NQIV88LL7yAGTNm2F6TkpKCG264ARkZGfj+++8RGhoqLCg0Zs/Ro0dFI8RNmzahf//+WLRoESZNmiSOGY1GLFmyBF999ZXDaxQKBc4880wsXryYRUw7YHcS4/asylqFm/68CQtXLsT6E+u7ZokxWwJ305c2LId+l44NmAlX88yv+5B+vBKhfiq8e9loqDgOhmHciurqaiEiVq1aJWJazjjjDMyZM0d0abbntddew9ixY8U5N998M2666SYHaw7xyCOP4N5778XOnTsxYMAAXHrppTAYpJ53u3fvRkVFhbhGY8aPH4+1a9e6+F/qHbAlhnF7XtryEvJqGsy/lJE0KGxQ5y6290dgyTVNx8l9lHwSXMk3m7Px1aZs4aV685JRiAtxbfo2w/SUS4csIj11764yYsQIsVh55pln8OOPP2LZsmW49dZbbeMkdEi8EA888ADeeOMNrF69GgMHNrQlIQFjdQk99dRTGDJkiLDepKWlISsrS1hdoqKimswhLi4OOTk5Ii5G7oIK4d4EixjGrSmuKxYCRgYZnp7yNBQyhRAwfiQ6OsPub6V1eH8g0GJSJlUx7OKGyrwuYMuxUjz+s+QKu3vGAEwbEOmyezFMT0Kumfa4dNzZEkOBvf/3f/8nYmfIclJXV9fEEjN8+HDbttUdVVhY2OI5sbFS7BudQyKGrqnRaJqN4fH19RUCpr6+XmwzLcMihvGIGJh+If1wXup5XbuYtqKhBsy8L4CoTlpzOsjx8jos/GIb9EYzzhoWi1u5HgzDuC1kPVm5ciVeffVVpKamChFx0UUXQaeTKoRbUVFjYTtIjDTOKLI/xypWrOdERESgtrZWXFetVju8rrS0FP7+/ixg2gGLGMYjGjt2KQZGdKDWAhmrAKNO6kLdTQKmVmfA9Z9tRUmNDoNjg0Q6NffQYhj35d9//8XVV1+N888/32aZobovzsZaP2bfvn1Nasmkp6dj1KhRTr+nN8IihvHubKTNHwDL73McG9JFi04HUk3vW7Ib+/IqEe6vxgdXjYWfmv/kGMadoSyiH374QQTz0heOxx57zCU1WyIjIzF69GisW7euiYihoN6ZM12faOANcMQQ41aUacuQW5UrlpyqnK7VhSErzKZF0rZPMOAXIXWlHn0luoOXfz+I/9udB5VChkVXjEE8B/IyjNvz+uuvi5TpyZMnCyEza9YsITZcAdWDaZxiffz4caxfvx7XXNNMAgLTBJmZvi56IZWVlSL/n1LYgoKCeno6TDtYk7MGt6++HSaz47cetVyNjfM3QiXvYHfn7I3Ax7MA6q907yGX14Cx54sNx/DYz3vF9isXDcfcsYnddm+G6S60Wi0yMzNFiX0fH9cFxnsrFNxL2UxUcM9aP4YyncrKyvD++++jt/7uVHbg+e31tm2DySAWxv17Ij254UkhYDQKDeQyyUhIWUkXDriwfQJm/y9AvhRDI8hc2+A+6kYB88fefDyxTBIwd58+gAUMwzDNQkHDn3/+OYqLi21jlHJ9991380+snXi9iDlp8UlQ+Hpuul9vIyUoBUvOWSKETIco2At8e3nzx0bOR3exLasMt32zAyYzcOn4RO5MzTBMq0yf7tho9p577uGfWAfwehHDeA4+Ch88e9KzHRcwxMb3pHXMMCBxYsN45EAgeTK6g4zCaiz4bAvqDSacmhaFZ84dyplIDMMwLsTrRczvF/7OMTEego/Sp30ChrpOZ/4DGPXSvkkP7P5O2j7zNSBpArqbrJIazP9wI8pq9RieEIx3LxsFJbcUYBiGcSleL2KCNcEI0nBgr1ex7Dap51Fj4kYBieO7fTq5ZbW47INNKKisx8DoQHx6zXhOpWYYhukGvF7EMF5GfrpFwMiAxAlSywBCoQZOeaRhv7umU6HF/A83iaq8fSP88eWCCQjzd6y+yTAMw7gGFjGM+5O1XmrcSNUAcjZJY0POB+Z+0qPTKqqqx2UfbkRWSS2Swvzw9fUTERnYiXgehmEYplOwiGHcm5pi4JtLpL5HNmTAtAd6cFJSP6QrPtyEo8U1iAv2wdfXT0BMMNfJYBiG6U5YxDDu5y6yFyzbPpH2IwZI1hcifgwQldZjU8wsrsHlFhcSVeH9asEEJIR2sqs2wzAM02m47QDjPmz/HFg0Bfj0zIZlz/fSsXPeAU55WFoGzOqxKR7Ir8TcRRtsMTDfL5yElAj/HpsPwzCug3on/fTTTx41p2PHjolzdu7c2e5rUsPL885z7ClXUlIiCu9Zm19So8qEhATU1NTAnWARw/QcFOOSuxU48hewbxmw3OIiCkmWLC/WZfpDQJJd7ZceYkd2Geb9byOKq+sxKDYI3944CXHcD4lhPJaioiLcdNNNSEpKgkajQUxMjOiVRJ2siby8PMyePRvuRHfN6bnnnsO5556LlJQUsT948GBMnDhR9JZyJ9idxPScgPnpJmDXN47jfaYCV/wMyN1LX/+2Jw93fbtTFLIbnRSCT64ej2C/DvZyYpheALXjM9TX98i9lRpNhwpMXnjhhdDpdPjss8/Qt29fFBQUYNWqVcIKQZCocTe6Y061tbX46KOP8PvvvzuMU1PK66+/Hg899BCUSveQD+4xC8a7oTb2m98H8nY1jNUWA4f/AGQKIGqwNOYXBpy/yK0EDH0gL1pzFC+tOCD2qRLvO5eOgr+G/3QYpjlIwLx91UU98sO5/bMlULWzEWV5eTnWrl2Lv//+G9OmTRNjycnJGD++odYUCaIff/xRuFrIrULNCpcuXYp33nkHmzZtQv/+/bFo0SJb88Ynn3xSuHrsXTlvvvmmWKxuGbrf/fffj71790KlUmHIkCH4+uuvxb2J9957D6+++ipycnLE/R599FFcccUVzc6J2Lx5M2688Ubs378fQ4cOxSOPPAJ7jEYjbrjhBvz111/Iz88XVqebb74Zd9xxB1rit99+E5YpsrzYc/rpp6O0tBRr1qzBaaedBneAP4kZ11B0CKg6IW3v+hbY9XXz5535MjBugVu+CzqDCY/9lI5vt+aI/asnp+CxswdDIe/eWjQMwzifgIAAsZDooIc1PbTbA4kEEhkkYGj70ksvRUZGRrssEwaDQYgPsmZ88803wgpEIsRqPSJxQuKCRM+MGTPw66+/CusHxaKccsopTa5XXV2Ns88+W4iLL7/8UnSFbixOTCaTeP3333+P8PBwrF+/Xoia2NhYXHzxxc3Ok8TdmDFjmoyr1WqMHDlSHGcRw3gHJiNQdkxyD1nZ8QXw75uO55HFZcrtgG9ow1j0ECB1BtyRwiotbvt6BzZlloI0yxNzhuCqyZJvmGGY1l06ZBHpqXu3+1ylEp9++qkQFGRNGT16tLDIXHLJJRg+fHiLr7v33ntx1llnie2nnnpKWFJIxKSltZ0xWVlZiYqKCiE8+vXrJ8YGDRpkO07iiIJsyVJCUDfrjRs3ivHmRMzXX38tRAq5fnx8fMRccnNzRZyPFbL20DytkHVnw4YN+O6771oUMVlZWYiLi2v2GI3TcXeBLTFM+yCRUkDpz5UNY3WlwJ9PASWHm39NZJokXlS+wMl3A2nSH767s/FoiehETcXs/NUKvHvZaJySFtXT02IYj4CsCu116fQ0FBNDgoQsCyQWli9fjpdffhkffvihEBPNYS9wyJpBFBYWtkvEhIWFietS8DBZT8jaQkLCeh1yCZGVxJ4pU6bgrbfeavZ6+/fvF/MhAWPF6tqy5z//+Q8+/vhjZGdno66uTliAyKLSEnSO/TXt8fX1FTEz7gKLGKZ58nYDB5cDRp20f2wdkLOx+XOpaaPS7hfeJxiY8QQwrGf84p3FZDLjf/8cxSu/H4DJDAyIDsB7l49Bv8iAnp4awzAugh7WJChoeeyxx7BgwQI88cQTLYoYsmxYsbqByBpCyOVyEUdnj15vaVRr4ZNPPsHtt9+OFStW4NtvvxUxLytXrmwSf+IsFi9eLKxHr732mhA4gYGBeOWVV0RMT0tERESgrKys2WMUE2O1IrkDLGK8Gb0WqMpzHKsrA3I2A/VV0r62HDi80vE8+iPUWY7bQ0IlJMluQAb0ORk49VFHN5EHQlaXB5fuxqoDhWL/glHxePb8odzIkWF6GZRK3NnaMJGRkSJ4loSMVeA0V69l1KhRYqEsHxIW5BYiEUOuJUrvvuqqq2zn0j7NqTkGDRqEL774Alqt1mY5IYuSPfT6yZMn21xUxJEjR1r9d9DcKMamOdLT03HRRe7zBdX7RUx1MSDXWmI2zIDZZLdt2bdu28ZaO97o9bYxtP81Ld2Txo31FN4PGPWA2QjUlgIGreO/iQRIdaF03IrJAJTnSCJFTMUEVOQ6ntMR5Cpg4GwgyOIX9QkBxlwNBElmT29i+Z48PPJTOkprdFAr5HjynCG4dHxih1I1GYbxLCiNeu7cubj22muFS4YsFFu3bhXuJKqP0hmmT58uas/QNehBT9YWclEFBQWJ4xR4+/777+Occ84RsSUHDx7E4cOHceWVV4rj9913n3AvkYggV9Mvv/yCH374AX/++Wez97vssstEcLE17ZkyoCh+xh4KQP78889FujTFw5Do2bJli9huCXJ30fXIGhMa2vAFla5//PhxMTd3wftFzDsjAU0vfhgpfQG5omGfuj0njGsQI3IlkDQJiB3p2AGa0p093LrSFhV1ejy5bC9+3HFc7KfFBOKNeSNFITuGYbwbykyaMGEC3njjDWGZILdPYmKiEAQPP/xwp65JlpH//ve/eP755/HMM8+ImBty5ZBwIfz8/HDgwAFRl4ZEFMXC3HLLLSJFmqDMJYp/ISFCWUYkNMj9ROKopX/DL7/8goULFwrhQxabl156SdzXCl17x44dmDdvnvhiRtlUZJUhcdUSw4YNE4HOFPxrnRtBGVUzZ860pYO7AzJzYweel0BR4MHBwah4MBBBQsTIAJnc8qCmbcu+dbvJcbRxvJOvaXK80bnkslFqJHFB+yQk1I368qj8gIBoQGFfbE0mWU1o3CpGaD8o3lGcMMLU+9uefDz9614UVNaL7KObpvfDHacNgFrpPjVqGMbdITcGWRfoYdtSICjjmfzf//2fsAyR+4hifSgYmKw65PqiYGNX/u7Ynt8VFTYrVu+1xDyUC7TxQ2B6D0eLqvHEsr1Ye7hY7PeJ8Merc0dgTLJ3W50YhmE6AmVtkauL3EdkoaLMJrJQOUPAOBPvFzEMA6Cm3oD3/j6C9/85Cp3RJCwuC6f1w83T+8FHZeduYxiGYQR33nmntAEgNTVVLO4GixjGq6k3GPHNpmy881cGSmqkdPFpAyLx1DlDuPs0wzCMh8MihvFKjCYzft55HK+vPITcsjoxlhLuhwdnD8KsIdGcecQwDOMFsIhhvAqt3oil23PxwT9HcaxEqioZFajBHTP64+KxiVApOHCXYRjGW2ARw3gF5bU6fLUpG5/8ewzF1fViLMRPhRun9hONG33VHPfCMAzjbbCIYTw6VXpbVhm+3pSN/9uTh3qDVPo7PsQX153UB5eMT+SKuwzDMF4MixjG48iv0OLX3Sfw3dYcHCqoto0Pjg3CgpP7YM6IOHYbMQzD9AJYxDAeAbUEWJ6eh2U7T2DzsVKpYwN1Q1DJMWd4HOZPTMaIhGAO2GUYxmP46KOPRBPIP/74Q+w/+OCDqKmpwTvvvNPTU/MYWMQwbusqOlhQhVX7C/HXgULsyC4TnaWtjE0Oxbkj43DOyHgE+9pXLmYYhmkf1Ofo8ccfF9VpCwoKRJ+gESNGiDEq6kZl+n/88UfRDsDZUMVa6pr9/fff28aoRUHfvn1x1113iTXTNixiGLchp7QWG46WYCMtR0pwosKx8eWQuCDhKjp7eCwSQhu1YmAYhukg1GOIyulTLyMSDSRkVq1aJfoauZolS5aIkvr2FXAjIiJE88X33nsPr7zyisvn4A2wiGF6hCqtHnuOV2B3Li3l2JVTgePlUj0XKxqlHFNSI3BqWhROSYsSAbsMw7i/FdWsl4LsuxuZSt5ul3J5eTnWrl2Lv//+G9OmTRNj1Nhw/PjxYjslJUWszz//fNsx6uJM/Pzzz3jqqaewb98+0Y36qquuEt2klUrpkUpzoEaQy5YtE9enRo/WztZWFi9ejDlz5jSZF43RtVjEtA8WMYxL0RlMOFZSg4zCahwuqEZGUTX2najA0eIaW1yL7ZdRLsPwhGBM7BuOSf3CMTY5jFOjGcbDIAFz4vH1PXLvuKcnQ9bOcgrUAZqWn376CRMnToRGo3E4vmXLFkRFRYku0meccQYUCum6JHyuvPJKvP322zj55JNFB+wbbrhBHHviiSdsrydX0Ysvvii6Un/xxRe45JJLsGfPHtHpmli3bh2uuOKKJvMiEZWbmysEk1VIMR4sYv7zn/8IRZqfny98lRTwZFXKTM9jMplRWFWP3LJaYUmh6ri0iO3SWmSV1orquc1BlhUSLcMTQkRQ7ojEEPhr3P5XkmEYL4CsJp9++imuv/56LFq0CKNHjxYWGRIbw4cPR2RkpDgvJCQEMTExtteRBYYCcMn6QpAb6plnnsH999/vIGLmzp2LBQsWiG06vnLlSvH8IgsNWYGoQzNZcRpjHcvKymIR0w7c+olBUdt33323+AWbMGEC3nzzTeEvPHjwoFDIjPOtJuTmqdQaxLrKsqZ9yg4qrqoX/YeomFxxtQ4l1fVi3NCCSLESoFGiX1QA+kcFIDUqAAOjAzEsIRgRAY7ffBiG8XzIpUMWkZ66d0djYqhbM1lXNm7ciOXLlwu3z4cffoirr7662dfs2rUL//77L5577jnbmNFoFIG6tbW18POT4vUmTZrk8Dra37lzp9iuq5Nc5z4+Pk2u7+sruc3pWoyHi5jXX39dqORrrrlG7JOYoSjyjz/+WCjh9gaLBuqVMNP/LM9as8Vv27BtPdvxHLE2i9GGbdvxhnOtx6zjDa+zHrOOWs9r7hzHc8nCQeLAaFubGvaNjuP6RvsGyzlUgl9rMEKrN0nbepNoiGjdth6v00nixVosrqMo5DLEBvsIy0p8qK8Iuk2wbPeLDEB0kIZTnxmml0DxIO116bgDJCROP/10sZALiKwnZFFpScRUV1cLa8wFF1zQ7LXaQ3h4uPg5lZWVNTlWWloq1lZLEOOhIoYixrdt24aHHnrINiaXyzFjxgxs2LChyfn19fVisVJZWSnWs99aC7mGM1k6gr9agSBfFQJ9lAj0kdZh/mphOQm3rgOkNS2RgRohZBiGYTydwYMHizgZQqVSCSuLPeR2Im9Aampqq9chyw7Fztjvjxo1Smyr1WpxHwoMnjlzpsPr0tPTxX2HDBnixH+V9+K2Iqa4uFj88kRHRzuM0/6BAweanP/CCy8IddwYKoamtHwroMesNXJd/Nfy3HUYl9mGpW8Utu2GFzieQ9vNv9aKGLfdS2a33fJ8lHK5EAZKhUxay61rueO+OC63Oy6tqdGhRiWHj1IBHxUtcrGmjB/bvlIBjUoBX5VCCJUgHxUCfJQsSBiG8XoojZriVq699loRAxMYGIitW7cKd9K5554rzqHAWkq5pjRoCvylOjJUQ+bss89GUlKSyDaiL9fkYiLx8eyzz9quT/Vfxo4di5NOOglfffUVNm/eLIrbWaHQCAruvfPOOx3mRa4tChi2upWYNjC7KcePHxeenvXr1zuM33fffebx48c3OV+r1ZorKipsS05Ojng9bTMMwzDOp66uzrxv3z6x9jTomfHggw+aR48ebQ4ODjb7+fmZBw4caH700UfNtbW14pxly5aZU1NTzUql0pycnGx77YoVK8yTJ082+/r6moOCgsQz6f3337cdp2fPf/7zH/Ppp59u1mg05pSUFPO3337rcP+9e/eK15eXlzuM0xy++eYbc2/+3amoqGj389ttLTFU9IdS2qj4kD20bx8pboVUcuMUOYZhGIZpDnpekAWflpagmi3N1XIhKwotrUFZRtZ2As1B7iQKKqZsJWvYBAUWk2XHvp4M0zodC+XuRshnOGbMGGHKs2IymcR+46hvhmEYhvE0qHwI1aqxQn2TqC6NtWge0zZu/ZOi9GrKxSe/ItWGoRRrepOt2UoMwzAM46lQzM1tt91m22cLjJeJmHnz5tkadFGxu5EjR2LFihVNgn0ZhmEYxl2wlvBgermIIW699VaxMAzDMAzDeERMDMMwDOMZsOWB6anfGRYxDMMwTKegomwEl8hnOlPQlrA21vRadxLDMAzjntADiBokFhYWin3qG2Rf6JNhmoMyjSnelX5fupqJxSKGYRiG6TTWul1WIcMw7YHq4VDV466KXhYxDMMwTKehh1BsbCyioqKg1+v5J8m0uxYcCZmuwiKGYRiGcYprqavxDQzTUTiwl2EYhmEYj4RFDMMwDMMwHgmLGIZhGIZhPBKltxfSqays7OmpMAzDMAzTTqzP7fYUxPNaEVNVVSXWiYmJPT0VhmEYhmE68RwPDg5u9RyZ2UvrRVMxnRMnTiAwMNCjiy+RIiUhlpOTg6CgoJ6eTq+G3wv3gd8L94LfD/eh0gueGSRLSMDExcW1mYbttZYY+ocnJCTAW6BfRk/9hfQ2+L1wH/i9cC/4/XAfgjz8mdGWBcYKB/YyDMMwDOORsIhhGIZhGMYjYRHj5mg0GjzxxBNizfB7wfDfhTvCn1Pug6aXPTO8NrCXYRiGYRjvhi0xDMMwDMN4JCxiGIZhGIbxSFjEMAzDMAzjkbCIYRiGYRjGI2ER4wa89957GD58uK040aRJk7B8+XLbca1Wi1tuuQXh4eEICAjAhRdeiIKCgh6dc2/hxRdfFBWf77zzTtsYvx/dx5NPPil+/vZLWloavxc9xPHjx3H55ZeLzyJfX18MGzYMW7dutR2nPJHHH38csbGx4viMGTNw+PDhnpqu15KSktLk70Imk4nnRG/7jGIR4wZQZWF6WG7btk18IJx66qk499xzsXfvXnH8rrvuwi+//ILvv/8ea9asEe0ULrjggp6ettezZcsW/O9//xMC0x5+P7qXIUOGIC8vz7asW7eO34seoKysDFOmTIFKpRJfsvbt24fXXnsNoaGhtnNefvllvP3221i0aBE2bdoEf39/zJo1SzxUGed+Ntn/TaxcuVKMz507t/d9RlGKNeN+hIaGmj/88ENzeXm5WaVSmb///nvbsf3791NavHnDhg09Okdvpqqqyty/f3/zypUrzdOmTTPfcccdYpzfj+7liSeeMI8YMaLZY/xedC8PPPCA+aSTTmrxuMlkMsfExJhfeeUVh/dIo9GYv/nmm26aZe/kjjvuMPfr10+8B73t74ItMW6G0WjE4sWLUVNTI9xKZJ3R6/XCLGuFzOlJSUnYsGFDj87VmyFT7FlnneXwcyf4/eh+yB1BjeD69u2L+fPnIzs7m9+LHmDZsmUYO3as+LYfFRWFUaNG4YMPPrAdz8zMRH5+vsPfDPW/mTBhAn9WuRCdTocvv/wS1157rXAp9bbPKBYxbsKePXuE75KqLC5cuBA//vgjBg8eLD4U1Go1QkJCHM6Pjo4WxxjnQyJy+/bteOGFF5oc4/eje6EH4KeffooVK1aI2DF6UJ588smiwy2/F93L0aNHxXvQv39//P7777jppptw++2347PPPhPHrZ9H9NlkD39WuZaffvoJ5eXluPrqq23vQ296ZnhtF2tPY+DAgdi5cycqKiqwZMkSXHXVVcKXyXQv1L7+jjvuED5mHx8f/vH3MLNnz7ZtU2wSiZrk5GR89913InCU6T5MJpOwxDz//PNinywx6enpIv6FPq+YnuGjjz4SfydkreyNsCXGTSDlnJqaijFjxggLwIgRI/DWW28hJiZGmAtJadtDkeZ0jHEuZIotLCzE6NGjoVQqxUJikoIVaZu+zfD70XPQt8sBAwYgIyOD/za6Gco4IuuwPYMGDbK596yfR42zYPizynVkZWXhzz//xIIFC2xjve2ZwSLGjb/11NfXC1FD2QCrVq2yHTt48KD44KCYGca5nHbaacK1R1Yx60LfPikWw7rN70fPUV1djSNHjogHKv9tdC+UmUSfPfYcOnRIWMaIPn36iIek/WdVZWWlyFLizyrX8Mknn4j4JIrfs9Lr/i56OrKYMZsffPBB85o1a8yZmZnm3bt3i32ZTGb+448/xI9n4cKF5qSkJPNff/1l3rp1q3nSpEliYboH++wkfj+6l3vuucf8999/i7+Nf//91zxjxgxzRESEubCwkN+Lbmbz5s1mpVJpfu6558yHDx82f/XVV2Y/Pz/zl19+aTvnxRdfNIeEhJh//vln8Vl27rnnmvv06WOuq6vr7ul6PUajUTwXKGusMb3pmcEixg249tprzcnJyWa1Wm2OjIw0n3baaTYBQ9AHwM033yzSrulD4/zzzzfn5eX16Jx7s4jh96P7mDdvnjk2Nlb8bcTHx4v9jIwMfi96iF9++cU8dOhQkTadlpZmfv/99x2OU4rvY489Zo6Ojhbn0GfZwYMHe2q6Xs3vv/8u0qab+/n2ps8oGf2np61BDMMwDMMwHYVjYhiGYRiG8UhYxDAMwzAM45GwiGEYhmEYxiNhEcMwDMMwjEfCIoZhGIZhGI+ERQzDMAzDMB4JixiGYRiGYTwSFjEMwzAMw3gkLGIYhmEYhvFIWMQwDMMwDOORsIhhGMZr+PXXX0U35fHjx+Pw4cM9PR2GYVwM905iGMZrGDhwIP7zn/9g79692LBhAxYvXtzTU2IYxoWwJYZhGK8hPDwcqampSElJgVqt7unpMAzjYpSuvgHDMIwruOaaaxAfH49nn33WYaxfv36Ijo5Geno6/+AZxsthdxLDMB6H0WhETEwM/u///k/EvxAGgwEjR47EnDlzhEupoqICMpmsp6fKMIwLYXcSwzBuxbFjx4T4aLxMnz7dds769euhUqkwbtw429iiRYvQt29f3HLLLaiqqsLRo0d76F/AMEx3we4khmHcisTEROTl5dn28/PzMWPGDEydOtU2tmzZMmFxsVpaSktL8cwzz+Dvv/9GQkICgoODsXPnTuFaYhjGe2FLDMMwboVCoRCuIlpCQkKwcOFCTJo0CU8++aTtnJ9//hnnnHOObf+JJ57A+eefj0GDBon9wYMHY9euXT0yf4Zhug+2xDAM47Zce+21wjW0cuVKyOXSd679+/fjxIkTOO2008T+vn378OWXX4pxK0OHDhWWGIZhvBsWMQzDuCWUdfT7779j8+bNCAwMdHAlnX766fDx8RH7d911F8rLy4UbyYrJZBJuKYZhvBsWMQzDuB1Lly7F008/jeXLlzeJayFX0g033GCr0Ltt2zbs2LEDSmXDx9mWLVuEFaesrAyhoaHdPn+GYboHTrFmGMatoPouEyZMwN133y0yjaxQ8TpKoyaLC7mTKHiX3EYkVh544AGHa2RnZyM5ORmrV692yGpiGMa74MBehmHciq1bt6K2tla4k2JjY23LBRdcgF9++UXUhYmIiMA777wj3Ei33nprk2uQK8nPz4/jYhjGy2FLDMMwHgNlJJ100km4//77e3oqDMO4AWyJYRjGYyABc+mll/b0NBiGcRPYEsMwDMMwjEfClhiGYRiGYTwSFjEMwzAMw3gkLGIYhmEYhvFIWMQwDMMwDOORsIhhGIZhGMYjYRHDMAzDMIxHwiKGYRiGYRiPhEUMwzAMwzAeCYsYhmEYhmHgifw/fsPzPmKS9uoAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# Create the structure, with different forms of interfacial roughness\n", "fig, ax = plt.subplots(1)\n", "\n", "for i, rough_model in enumerate(\n", " [None, Erf, Linear, Exponential, Tanh, Sinusoidal, Step]\n", "):\n", " if rough_model is not None:\n", " interface = rough_model()\n", " else:\n", " # using None for the interface roughness results in the Nevot-Croce default\n", " interface = None\n", "\n", " # Here we create the d2o slab via the Slab constructor, which allows you to specify what roughness type you want to use via the `interface` keyword.\n", " s = air(0, 0) | Slab(0, d2o, 5, interface=interface)\n", "\n", " # you could also set the interfaces attribute directly\n", " # s[-1].interfaces = interface\n", "\n", " zed, sld = s.sld_profile(z=np.linspace(-40, 100, 1001))\n", " if rough_model == None:\n", " zed += 46\n", " ax.plot(zed, sld + i * 0.5, label=repr(interface))\n", "\n", "ax.set_xlim(25, 75)\n", "ax.set_xlabel(\"z/$\\\\AA$\")\n", "ax.set_ylabel(\"$\\\\rho / 10^{-6} \\\\AA^{-2}$\")\n", "plt.legend();" ] } ], "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.14.0" } }, "nbformat": 4, "nbformat_minor": 5 }