{ "cells": [ { "cell_type": "markdown", "metadata": { "deletable": true, "editable": true }, "source": [ "## Tutorial on how to implement periodic boundaries" ] }, { "cell_type": "markdown", "metadata": { "deletable": true, "editable": true }, "source": [ "This tutorial will show how to implement Periodic boundary conditions (where particles that leave the domain on one side enter again on the other side) can be implemented in Parcels" ] }, { "cell_type": "markdown", "metadata": { "deletable": true, "editable": true }, "source": [ "The idea in Parcels is to do two things:\n", "1) Extend the fieldset with a small 'halo'\n", "2) Add a periodic boundary kernel to the `.execute`" ] }, { "cell_type": "markdown", "metadata": { "deletable": true, "editable": true }, "source": [ "We'll start by importing the relevant modules" ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "collapsed": true, "deletable": true, "editable": true }, "outputs": [], "source": [ "from parcels import FieldSet, ParticleSet, JITParticle, plotTrajectoriesFile\n", "from parcels import AdvectionRK4\n", "from datetime import timedelta as delta\n", "import math" ] }, { "cell_type": "markdown", "metadata": { "deletable": true, "editable": true }, "source": [ "We import the Peninsula fieldset; note that we need to set `allow_time_extrapolation` because the Peninsula fieldset has only one time snapshot." ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "collapsed": false, "deletable": true, "editable": true }, "outputs": [], "source": [ "fieldset = FieldSet.from_parcels('Peninsula_data/peninsula', allow_time_extrapolation=True)" ] }, { "cell_type": "markdown", "metadata": { "deletable": true, "editable": true }, "source": [ "Extending the fieldset with a halo is very simply done using the `add_periodic_halo()` method. Halos can be added either in the zonal direction, the meridional direction, or both, by setting `zonal` and/or `meridional` to `True`" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "collapsed": false, "deletable": true, "editable": true }, "outputs": [], "source": [ "fieldset.add_periodic_halo(zonal=True)" ] }, { "cell_type": "markdown", "metadata": { "deletable": true, "editable": true }, "source": [ "The other item we need is a custom Kernel that can move the particle from one side of the domain to the other." ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "collapsed": false, "deletable": true, "editable": true }, "outputs": [], "source": [ "def periodicBC(particle, fieldset, time, dt):\n", " if particle.lon < fieldset.halo_west:\n", " particle.lon += fieldset.halo_east - fieldset.halo_west\n", " elif particle.lon > fieldset.halo_east:\n", " particle.lon -= fieldset.halo_east - fieldset.halo_west\n" ] }, { "cell_type": "markdown", "metadata": { "deletable": true, "editable": true }, "source": [ "The fieldset constants `halo_east` and `halo_west` in the kernel above were set during the call to `add_periodic_halo()`, and store the original zonal extent of the fieldset. If we would have also used `meridional=True` in `add_periodic_halo()`, `fieldset.halo_north` and `fieldset.halo_south` would also have been added.\n", "\n", "***Note that some hydrodynamic data, such as the global ORCA grid used in NEMO, already has a halo.*** In these cases, **do not** extent the fieldset with the halo but only add the periodic boundary kernel, where you use the explicit values for halo_east and halo_west" ] }, { "cell_type": "markdown", "metadata": { "deletable": true, "editable": true }, "source": [ "Now define a particle set and execute it as usual" ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "collapsed": false, "deletable": true, "editable": true }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Compiled JITParticleAdvectionRK4periodicBC ==> /var/folders/r2/8593q8z93kd7t4j9kbb_f7p00000gn/T/parcels-501/e314053eec65b2ef80b5495b2134f745.so\n" ] } ], "source": [ "pset = ParticleSet.from_line(fieldset, pclass=JITParticle, \n", " size=10, start=(0.7, 0.1), finish=(0.7, 0.4))\n", "output_file = pset.ParticleFile(name=\"PeriodicParticle\", outputdt=delta(hours=1))\n", "pset.execute(AdvectionRK4 + pset.Kernel(periodicBC),\n", " runtime=delta(hours=24), dt=delta(minutes=5),\n", " output_file=output_file)" ] }, { "cell_type": "markdown", "metadata": { "deletable": true, "editable": true }, "source": [ "And finally plot the particle trajectories" ] }, { "cell_type": "code", "execution_count": 6, "metadata": { "collapsed": false, "deletable": true, "editable": true }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAikAAAF5CAYAAABa9bp0AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAAPYQAAD2EBqD+naQAAIABJREFUeJzsnXl4HEeZ/z81l46Z0WVJlmTJlizft+WQOHEukmXDkmVJ\n7EBiyMkRWEJgw0J+u7AQrgQ22QQIy7LALglnuOKQEAgECEtCYuewbPlObEs+ZFmyZEsa3XPV74+a\nq0cjW5IlayS9n+epp7urq0vVo5nub7311ltKa40gCIIgCEK6YZvoBgiCIAiCIKRCRIogCIIgCGmJ\niBRBEARBENISESmCIAiCIKQlIlIEQRAEQUhLRKQIgiAIgpCWiEgRBEEQBCEtEZEiCIIgCEJaIiJF\nEARBEIS0RESKIAiCIAhpSdqIFKXUHUqpBqVUn1Jqi1LqTcO87galVFgptSkp/5FIfmL67fi0XhAE\nQRCEsSYtRIpS6nrgQeAeYDVQB/xeKVV4huvmAA8Azw9R5BlgJlASSRvHqs2CIAiCIIwvaSFSgLuA\nb2utf6C13gd8COgF3jvUBUopG/Aj4LNAwxDFBrTWrVrrE5HUOdYNFwRBEARhfJhwkaKUcgJrgD9F\n87RZmvmPwIWnufQe4ITW+pHTlLlcKdWilNqnlPovpVTBmDRaEARBEIRxxzHRDQAKATvQkpTfAixM\ndYFSah1wG7DyNPU+AzyOsbJUA18GfquUujAiggRBEARBSGPSQaQMhQIGiQmllAf4IfABrXX7UBdr\nrX+ecLhbKbUTOAhcDvw5Rb0zgKuAQ0D/2TRcEARBEKYZmUAl8Hut9cmxqjQdREobEMI4uCZSzGDr\nChiryBzg10opFcmzASil/MBCrfUgHxWtdYNSqg2YRwqRghEoPx7VHQiCIAiCAPAe4CdjVdmEixSt\ndUAptRW4EngKICI+rgQeTnHJXmB5Ut69gAf4KHA01d9RSpUDM4DjQzTlEMCPfvQjFi9ePLKbSFPu\nuusuvvrVr050M8aMqXQ/U+leQO4nnZlK9wJyP+nK3r17ufHGGyHyLh0rJlykRHgI+H5ErLyCme2T\nDTwKoJT6AdCotf6U1toP7Em8WCnVgfG33Rs5dmMcax8HmjHWk38H3gB+P0Qb+gEWL15MTU3NmN7c\nRJGbmztl7gWm1v1MpXsBuZ90ZirdC8j9TALG1F0iLUSK1vrnkZgoX8AM+2wHrtJat0aKlAPBEVQZ\nAlYANwN5QBNGnHxWax0Ys4YLgiAIgjBupIVIAdBa/xfwX0Ocu+IM196WdNwPvHXsWicIgiAIwrlm\nwuOkCIIgCIIgpEJEyhRm48aptQrAVLqfqXQvIPeTzkylewG5n+mGkrhmBqVUDbB169atU82JSRAE\nQRDGldraWtasWQOwRmtdO1b1iiVFEARBEIS0RESKIAiCIAhpiYgUQRAEQRDSEhEpgiAIgiCkJSJS\nBEEQBEFIS0SkCIIgCIKQlohIEQRBEAQhLRGRIgiCIAhCWiIiRRAEQRCEtEREiiAIgiAIaYmIFEEQ\nBEEQ0hIRKYIgCIIgpCUiUgRBEARBSEtEpAiCIAiCkJaISBEEQRAEIS0RkSIIgiAIQloiIkUQBEEQ\nhLRERIogCIIgCGmJiBRBENKblha4+GKorjbbEycmukWCIJwjHBPdAEEQ0oiWFtiwAY4fh9JS2LQJ\niotHXo/fDw0NcNNN0NwMRUXwjW9Abi4EgxAKmW00ne74nnvgjTdMvfX1cMkl8PnPg8MBdvvItu3t\n8LGPQWsrlJTAL34BFRWg1MR9VoIgDImIFEGY7gwMmJd3eztcfz3s3Gny6+th7Vr4wAegpyeeenut\nx6nyg0Hr3zh6FNatG5v2vvEGbNx49vUcOQJz5oDNBtnZJrnd8W3ifqq8b33LfEZgtlddBU8+Cfn5\n4PGMTvgIgmBBRIogTHZaWmD9ejh2DGbMgC98wbwgo8Ijmk6dGpzX3m7ExVAcOgQPPWR9SUdf1Lm5\nUFY2+CUeTXffbR2aKSuDX/7SatlITEPl/e3fwubN8XouvBCefTZubRnJ9rrroKkpXldREXzpS4NF\nVnS/t9d8bkePDj7f2Wn9rLZvN6IHTNvz8oxgSUyp8rQ2bWhvh1mz4Fe/EouMIEQQkSIIE8WZhgt6\ne81L/sQJUzZxm7i/b1/ccnH4MPz938frcLniL8OCArOdMwdWrbK+KKPnPvpR2LYtfv1FF8Ff/zq6\n+/vud60iparKCIyR8qtfGRGW+Dl5PKNrU1WVVaQsWAC33z66ui6+GF58MX68YgXcf39qIdjebj6L\n11+PHyeLHDDWndJSqKyEmTPN9yExJefNmGEsQVFkCEqYYohIEYRzSTBoXiTHjsEttxiBAWa4YPFi\n89KMipDubuu1SpmXUuLLaulSU1d7e7xceTls2WJER1bWyIYdfve7wYJgtGzaNDZ1FRePXiiNV5uG\nqmskgiAUMkJl1SpjqYmSlwfvfGdciG7bFt8fGLDWYbMZa1BUtNTVQVubOVdfD1dcAY8/bqxYXu/o\n71UQJggRKYIwUlL1VouKzLDAsWOmpx7dJu4fO2au1Tp1vX19sGgRXHbZ4F7zzJlQWGiGEZKprbX2\n6OfMMcMGo2EsBcFY1jVWpNP92e3GgjV7tlWkLF4MX/nK4PJaQ1dXXMQmW9ZSCdvdu813CoxIKSsz\nadas1PulpZCRIRYZIW1IG5GilLoD+ARQAtQBd2qtXx3GdTcAPwF+pbVen3TuC8D7gTzgReAftdYH\nxrrtwhQnHDYP7SNHzMvk7rvNzBUwvdXZs80LxO+3XldcHH8J1NTA299ufTF88IPwyivx8jU18Mgj\nI2/fWFoHhHPPcP9/SkFOjknz5qUukzwEVVNjfIqSBfPhw/DSS2a/v99aR2GhGWqM+irV15v2pZvg\nFKYFaSFSlFLXAw8CtwOvAHcBv1dKLdBat53mujnAA8DzKc79P+AjwC1AA/ClSJ2Ltdb+5PLCNOB0\nFpCjR02KCpHEdOwYBALxepKHTzIz4YtftAqQkhLjD3I6fv3r9BsOEc494z2cdToLiNbQ0THY8nf/\n/dZyx4+PTfsEYYQoPZTp+Vw2QqktwMta649FjhVwFHhYa33/ENfYgL8A3wMuBXITLSlKqSbgAa31\nVyPHOUALcIvW+ucp6qsBtm7dupWampoxvT9hAunvj4uPD30IDiQY0jIzjeDo64vnORzGp6OiInWa\nPRve8Q5rb3XdOhEJwtQi2SIj33HhDNTW1rJmzRqANVrr2rGqd8ItKUopJ7AGuC+ap7XWSqk/Aqeb\nCnAPcEJr/YhS6tKkOqsww0Z/SqjTp5R6OVLnIJEipDlDWUHa240AOXw49ba5eeg6MzJMoLBEATJz\npnW2RCpkeEWY6qTjd1z8ZKYlEy5SgELAjrFyJNICLEx1gVJqHXAbsHKIOksAPUSdJaNuqXBuCYWM\n+fnIEXjve61RRysrjeNhoqOgy2WExpw5sGwZvO1tZj+ad9NN1ngby5bBXXeNvF0yvCJMddLxO75h\nQ9y6U18Pf/d3xq8mI2Ni2yWMK+kgUoZCYYSGNVMpD/BD4ANa6/ZBV42izkTuuusucnNzLXkbN25k\n41hEuJyOnK73Ex2KOXw4dWpsHBy5NIrTCZ/9bFyAzJ5t6j2dFSRVvA1BECYHyX4xtbVmyLakxDwD\nhko5OfFrxBozJjz22GM89thjlrzOVHF/xoAJ90mJDPf0Ahu01k8l5D+K8TO5Nqn8SqAWCGFEB8QX\nSgwRt74cBFZprXckXPt/wDat9aDus/ikjAMdHfDmN5tInFEKC81CcYcPW4dilDIPjaEeNO9/v4n9\nEUXGyAVhepHsJ7N0KXz84/Hh3Wg6etTq6J6XF+/MvPaaVey86U3muXKmIV7hjExZnxStdUAptRW4\nEngKYo6zVwIPp7hkL7A8Ke9ewAN8FDiqtQ4qpZojdeyI1JkDXAB8czzuY0pyul5H4lBMKl+QI0fA\n5xtcZ3e3iQPx1rdaRUh5+enNtk8+KVYQQZjODHfmUihkOkBR0ZIoYtqSJou++qoJeDh7djxFLbPR\n/YoKY7GJItaYc8qEW1IAlFLvAr4PfJD4FOTrgEVa61al1A+ARq31p4a4/hEGz+65G/h/wK3AIeCL\nwFJgaaopyGJJScLng8svt4ZILyoyAiM6FBMKxc/l5w/+kc+ZA/feCzt2xMuJBUQQhIki2RqzeDF8\n+MODO1nJQ0szZ8afb5s3W5dWOP98Y42Z5gtKTllLCoDW+udKqULgC8BMYDtwlda6NVKkHBjCOWHI\nOu9XSmUD38YEc3sB+LtpEyPlTL4gjY2DY4IkHqeygnR1GYvHRRdZhUhFhXXcN5HLLxcLiCAI6cFw\nrTEDA/FnZLKluLXVWvaVV4yl5XRhCyoqzIKcUSEj1phhkxaWlHRgylhSBgaMyn/HO2Dnznh+QYGZ\nEXP06OAfWVFR6h/VffdZ6xAriCAI051ka8yiRcYak9zha2qyWps9nvgztq7OCJUoq1fDn/9sOnuT\n1CIzpS0pQgpSKe3sbKPuGxtNVMhU+8kCJEpfH6xZA9deaxUk5eVmTDYVV1whVhBBEIRERuIbc/z4\n4AjWR49aFwQFM6yel2eEzKxZ5rk81LaoyDj6ThNrjFhSIkyoJSUQMF+4aFjqpib48peN6Ihit1tV\nOZiZMuXl1i9wNN19t9WfRKwggiAI6UGyNWb5cvj0p62dzui2qckaisHpNMtvnDplhuCjLFgA//3f\n8eU5zvGq12JJmUxEFW5TkxES995r/EASRUhTk1HATU1m9dJEsehwDF4pNzcXvvGNuAgpK7N6nCfz\nu9+JFUQQBCEdGckaS+GweUckC5hvfMNa7o03jPU7isdj6k5c4Tq6Hz12OODmm9PaGiOWlAgxS8qz\nz1Lz+c8P/U+LLpfe3Dx0ev556OkZ/EdsNhN4KPmLknhcVmaEzaWXytoZwvig9emTlJEyk6VMOrbp\nXJU5fNi6grXTaSwu0WtHw1m8Z8SScq745CeNUxOY0MsrV8KFF1pFSOKCdMPB44E77jDiw2Yb/GU7\nedKkHTvieZdcYsYuu7qM2e7CC+Hzn0+fH4iUmZxlhKmDUqnT6c5N1jLnsi1Kmed0urRnqDJ9ffDU\nUyb2lMcD11xjtiOp52tfs/rHpOFq1yJSkkkO9tPcDE88kbrsUF+GQMD6Uujthf/8z9F9STMywO+H\nH/84/X80Y1Um+oBIh7ZImfQrk05tmagyggDwX/91dtc/+6zVYl9aenb1jQMiUpIpLDTjfVEuusgM\n3yQ/KE7HiRPDH28UBOG0tHS3sOHnGzjefZxSTymbrt9EsXt0v6exrEsQJj3puNp1EuKTEiHmk/KH\nP1Dzuc+JwBCEUaC1pi/Yx8FTB7n1yVtp6W6hMLuQe6+4F4/LQzAcJBAOEAgFCIQD5jiyHwgFBp0P\nhAJ8t/a7HPUdjf2NipwKbl9zO06bE4fNgdNutg6bI5aXmJ+Y90+//yd2ndgVq+tNZW/iTzf/CbfL\njU2NbP0WETyCEGe8fFJEpESYMsHcBGGUHO44zHW/uI7jXcfJz8rnX9b9CzZlo8vfhW/Ah2/AR9dA\nZN+fsD/gi5XpGugipENn/mNDEBUVTrsztm3tabXUaVd2CrMLY4ImGA7GxI4+/SLnp8Xj8uB1efFm\neMnJyInte12Dj70ZXu5/8X5eP/l67PoLZl3A5vdtRslwjDANEcdZQRAGcbre/EBwgNbeVlp7WjnR\nc4LW3sg24Tgxr8sfj7lwrOsYNz5xIwAKFXtxR1/WORk5eDO8lHpLyXHlWF7sORk5/POz/0xLTzyi\nZnlOOX+6+U8W60eyGLEre8oX/MXfu5gXj8bHzdeWr+Wv7009AyGswzGLTKKIieZd94vrqD0ef34u\nKVzCZy/7rEVkdfm7jACLCLGjvqMxQRY91xcc7Dz/8rGXybw3k8LsQoqyiyhyF1GUXUSxu9hynLjN\ny8zDpmxilRGEIRCRIgiTCK01Hf0dNPoaafQ18tFnPsqB9gMA1LfXM+/heRS7i2ntbcU3MHj9JbfT\nbV6akRfl4sLFXDr7Uordxfz7i/9Oa288YvHs3Nns/vBu3E73iK0D33rtWxaRMid3DgtmLBjVPW+6\nfhPrf7be8gIfCpuykeHIIIPUK2o/855nBtU1GjEQDAe5+HsX8/Kxl2N5C2cs5CPnf4TWntaYAGzq\naqKupY7WnlbaetsGWXocNgeF2YX4Bnz0BnoB839c+z9ruf8t91OeU84s7yxKvaU4bPK4FqYf8q0X\nhAkgVc+5MLuQ1p7WmABp9DVyrOvYoOPoyywVGs36xetjvfdEQVLkLiLbmT3ktU/se8IiUipyKvC4\nPKO6v5EIizNR7C4e0nIyUXU5bA6e2vjUiARPKBziVN+pmHUrcfvgSw9ayjZ0NPDOX7wzdmxTNma6\nZxrRkjOLcm9kGxEx0fxsZ7ZYZYQphfikRBCfFGG86fH3cLjzMIc6DnHnM3dS314fO5dhzzBDFeFA\nLM9pczIrZ1bsJRRNicc3/PIGXmp8KXbNuop1o34Jn+g5MSZWBmHkJA9pratYx1MbnzLC1HfMIlhj\nW98x2vuta8DkZebhD/ktQnbBjAV8/5rvU5lXyUz3TPGZEcYFcZwdZ0SkCMPhdL3UvkBfTIQ0tDdw\nqOMQhzoPxY4TrRTJ5GXmce8V91pESJG76IwzTkRYTA1G+3/sDfRyzHfMIlzue+E+fP7BQ30AmY5M\n5uTOoTKvMpaq8qpi+8Xu4piIEYuMMBJEpIwzIlKEMxEIBbjwfy9k6/GtsbzC7ELmFcyjob3B4oNh\nV3Zm584e9BKozKukKr+K6395PS8dHRsLiCAkkmyVuWDWBXzn7d8xojkqoCPi+VDHITr6O2JlMx2Z\nse/p9uPbae5pjp2T76hwOmR2jyCcA3r8PRxsP8jBUwet2/aDHO44PGh6bbe/m0WFi3hr9VstIqTM\nW3ZaR8cnrn9izHw2BCGRVP5Axe5iVsxckbJ8R38HhzsOx0TLoY5DNHQ0cKrvlKXc8e70C5kuTH1E\npAjTgpjpuus4he5C7rnsHtr72mMC5MCpAxw8ddBiDXE73VQXVFOdX836ReupLqjmP1/5T3a37o6V\nWVO6hkfe8ciI2zOWzqCCkMhIv1t5mXnkleSxsmSlJT/ZIlPqSb+Q6cLUR0SKMOXQWtPc3cyBUwfY\nf2o/B04d4Dtbv8PJvpMA1HfUc/VPrgbMA706v5rqgmreMvctsf3q/GrL+HyU9YvXiwVEmBaM5Qwt\nQRgtIlKEtGeo6bpNXU0cOHXAiJGT+znQfiB2HJ3doFBU5FbQE+ix1FnuLWfPHXvwZnhH1BaxgAjT\nhXT7rosj7/RERIqQtoTCIRp9jVz9k6tjQyz17fVUfa2KMGH6g/2AiSExO3c28wvms65iHbesvIV5\nBfOYVzCPuflzyXRkDjJdz8mbM2KBIgjCxLHh5xtiv+H69nqu/P6V/OY9v6E8p3zE6y4JkwcRKcKE\nEgwHOdxxOGYBOXDqQMwiUt9ejz/kH3SNw+bgC2/+QkyIVOZVkuFIHWE0ipiuBWFyk+y4u6t1F3O+\nNocMewZz8+fGngeJaXbubInUO8mR/54wbkTNs01dTeRn5vPxCz9OW2+bRYgc6jhEMBwETPCy6MPm\nquqrYg+aTz/3act6K8tnLudjaz82orakm+laEISRUeoptQRArCmt4Ytv/iIHTx2MPVN+s/831LfX\nx54pDpuDyrxK8yzJn8dM90x+susndPu7Kc8p51c3/EqGjNIciZMSQeKknB2d/Z2Dpu7+Ys8v6Bzo\ntJTLdGRSnV+dstdTkVOB3WYfVLcELBMEYbjPgWA4yNHOoymts3tb91rWT3LZXVxYfqHFYb66oJq5\n+XMpyCo4l7c36ZE4KcI5JdlJ7Zfv+iXBcJCDpw5S314fm7obFSWJMRVyM3KpLqiO9WaizM6dTcPH\nGkY8fixWEEEQhvsccNgcVOVXUZVfxVuq32I5V/31auo74taYbGc2Mz0zqWup44l9T1iWGcjLzIuL\nlry5MRGTk5HDR3/3UZq7m6XTdA4QkSLE8A34aGhvoKGjgY///uM0dDQAxkmt7MEySw+kPKec6vxq\nlhcv5x0L32HphUR7IMnOqhU5FeLgJgjChFHqLbWIlKVFS/nZdT+LHcdiJyV1xjYf3Uyjr3HQKtb1\n7fWs/NZK/mntPxlhlGfE0YysGbJG0hghwz0RptJwz1BT9fqD/bGw2A0dDbHw2NHjRGuIQll+kDOy\nZvDoNY9SnV9NVX4VmY7MM7ZDhmkEQUgnzuaZNBAc4FDHIS5/9HLLcgEum4tMZya+gfh6SR6XJyZY\nqvKqYktjRI97A71Tbjq1rN0zzkwFkdIX6ONI5xE2/HyDJSqq1+XF4/JYvOMdNgezc2fHfkCJvYCq\nvCrLdD+QdTsEQRAg9YrVL9z2Au397ZYOYENHQ2yJgUMdh2IhE8A8fxOHw6vyqvjqVV9lTt4c5uTO\nIS8zb9JZYsQnZRqRyhJSlF3Eqb5THO48zJHOIxzuOBzf7zzM4Y7DQ66yq7XmAzUfoCq/Krbg3ayc\nWaedmidTdgVBEAaT6tmolKIgq4CCrALWlK0ZdE1Yh2npbokJmDufudPi/9LQ0cA1P7smduxxeZiT\nO4fZubOZkzuHOXnW/VJPKW29bVPOGpOKtLGkKKXuAD4BlAB1wJ1a61eHKHst8ClgHuAE9gMPaq1/\nlFDmEeCWpEt/p7V+2xB1TqglJbrkeqOvkQ8+/UH2n9ofO5flyMKmbJaoqZmOzNiX1rLNm8Mnnv2E\nZaVesYIIgiCkD8nWmIsqLuLxdz3O4Q5rx/OI70gsL1HUOGwO7MrOQGggljc3by4PXfUQFbkVlOeU\nU5hdeE59AKe0JUUpdT3wIHA78ApwF/B7pdQCrXVbiktOAl8C9gF+4O3AI0qpFq31HxLKPQPcCkTt\nZgOMM6msINnObBp9jZZ0tPMojV3x4+QVRxNx2p3cc9k9FiVdlF00pDnwt+/5rVhBhClDSwts2ADH\nj0NpKWzaBMVTr8MoTCOGWqm6xFPCBeUXpLzGN+DjSOeRmCX9X//0rxaRUt9Rb7HGuOwuZnlnUZ5T\nHksVORWx/Qx7Bv/4239M+1lKaWFJUUptAV7WWn8scqyAo8DDWuv7h1nHVuBprfU9keNHgFyt9fph\nXj8mlpRkhWxXdkI6dMbrHDYHs7yzqMitYF/bPtp649pspnsm/7DwH1ARraWUQqEsW2BQ3kjLj0Ud\nQ5Ufz3ZPl785nu0ey7/ZekJx442K5mZFyUzFj3+syMsFv18RCCj8A5Gtn/j+gCLgV/j9Jn9gAHPe\nD/fdB/XxCRnMmwef+Qw4nSa5XPH9ofKix+3t8L73wYkTUFYmgkeYvKSyxmx616bBHWLfUctxorBJ\npCCrgI3LNlLmLaPUU0qZt8zse0uHnK2U2Cn3nvRS94U6mGqOs0opJ9ALbNBaP5WQ/yhGZFw7jDqu\nBH4FvENr/Vwk7xHgHUAAaAeeA/5Na53SZDFWIqX64WpLVMRMRyYXVVxE9HPWaLTWli1gyfOH/Lxx\n8g38IT9Ou5N5+fNw2B0jqmOovImoY6jyI22HMM3QyiSiWxL2E7bRsoPOpSo/OM9uV9hsYFMKm03F\ntzaTb4/s222mrD2Sr5QiHDJiLBRU2O0wq0zhcEwNUTwuf3Oq3c8Edix8Az4e2vIQHf0d5Gfm84kL\nP0FeVt5p6wA41XeKY13HeODFB/D54zOSTofL7qLUU0qpNyJePEa8PLr90bhrQhPwHWAKDvcUAnag\nJSm/BVg41EVKqRzgGJABBIEPRwVKhGeAx4EGoBr4MvBbpdSFehyVWXLo5jWla/jTzX8arz837UhX\nEZXudSSX7+7RtLVp2tqgtU3zwx9qWls1KA1oMjIhN1fT1a3p6zN55hyxfYdD4/FqPB7weDTuSHr1\n1eg1gNJ4czTvudGUt9s1TifY7Rq70xw7HBqHA+yR89GtzW7yH/gPTf3BeNvmVmvuvBOCIU0opAkE\nzTYY1IRCJj96LhjUhMLE9n/5uKa7O34vmVma1as1fb2a/n5NX7+mvw/6Bsy+3x+598i9xD4HpcnI\n0ASDpm6UJoDmRL3mgrUajxvzebg1mVk6+h85u+/wOfiehXU4bTs7U/W3Phac6DnBB57+wJjUlQp/\nyG/8ZDoPj9vfGIp0sKSUYsTGhVrrlxPy7wcu1lpfNMR1CqgCPMCVwGcxlpTnhyhfBRwErtRa/znF\n+Rpg66WXXkpubq7l3MaNG9m4ceOw7kdigwjnilS+GjNmmKGMxkY4dmxwiuZ3d1vrstkgHI4f5+TA\nRz4CBQWQn2/dRvezsiCFBZiLL4YX41Zo1q2Dv56F3/aJE7B+/dj4pIy0bcEgdHZCR4cZKkrcdnSY\noaiOjnh5ux20tn6WDodpd1mZSbNmxfcTj/v74brrxPdmOjIRQqy1p5X3PfU+WnpamOmeyXfe/h1m\nZM0YVh2//uWv+fWmX1N7vJbO/sjSJ/3AEUCGe4as57tAudb6705T5gTwaa31d1Ocm/RxUoSpz8AA\nHDoEBw/CnXdafTVcLgiFTIridMZfhLNmQXl5fD+aysrgb/5m7ITFWIqKsWas25ZK9PzlL0ZANjWZ\ndOxY6v2TJ611JQvFOXPggQdg7lyT8vNH305BGA8SO+VT1icFhnScPYJxnH1gmHX8L1Cltb5iiPPl\nwGGMteXpFOdFpAhpQXu7ESHRVF8f329sND31VOTlwZe/bBUgRUXm5Xcm0llYpDNn87n195vrouLl\nwx+2ChelrP/rvDwjVqqr48IlmioqjCAFmQ0lTAxTegoy8BDwfWVm6ESnIGcDjwIopX4ANGqtPxU5\n/hfgNczwTQZwNXAj8KHIeTdwD8YnpRkTT+XfgTeA35+rmxKEZFpazEutsdEMqdx2G7S2WoVI4vBB\nfn78pXThhfH96mrYuBFeeiledulS+NCHRteu4uKzG5KZrpzN55aZCVVVJgE8/LDVKnPRRfDUU+a7\nkZxeew2OHIlbzex2mD3bfDd27TLfMzBlr7nG+j0RhMlEWogUrfXPlVKFwBeAmcB24CqtdTSEajnG\nOTaKG/gzSMS+AAAgAElEQVRmJL8PEy/lPVrrX0bOh4AVwM1AHsbv+PfAZ7XWgXG+HUEAjN/HG2/A\nvn3w+utm+/TT0NsbL/PP/2xeLtXVsGqV6QEnCpHTmfifeGJwL16YvGzaNPj/GfUBOu+8weUDATh6\ndLCAaW+3ltu8GSorYeFCWLTIbKP7ZWWp/YoEIV1Ii+GedECGe4TTMZQJXWtjqk8UItH9o0fj15eU\nmJfC1q3Q1RXPr6qy+pUIwtmS7Cczb5757ka/lwcOGIdgAI8HFiwYLF7mz4fsbBk6EobPVB/uEYS0\n5tprTY8UjKhYutRYQN54Iz5Txuk0L4RFi+DGG60P/jwTvmDQC6Ss7NzehzD1SWWRSRQWgQA0NBjB\nEhXWr78Of/iDGXqMMns2+Hzx4cf6elOvDAsK5xIRKYKQQH+/eWDv2WPS7t1m+/rr1nJdXbB6Ndxw\ngxEjixYZq4jjDL+oVC8QQRhLzuQn43Qa68mCBfD2t1vPnToVFy+vv278ZBI5fhxBOKeISBGmNEOZ\nq/v6zEM4KkKiguTgwfg00NJSWLIErrrKmMcPHozXe9558D//M/L2iIOqkM4UFBgH7QsvNMcvvGC1\n/JWWTky7hOmLiBRhSpM8TLN4sXFGra+PT++cNcsM31x9tRElS5fGy0X59KfFAiJMP8TyJ0w0IlKE\nKUFfnxlb373bmpKdUvv64NZbjRBZssSIkaQAwykRC4gwHZHvvTDRiEgR0pZUQzVerxEj0eGZRDES\ntYxUVBgRcu218OSTZjZDlJoaePDBibkfQRBGj8w0mp6ISBHSkt5e4wtSV2eO6+tNmHC/P+4zUl5u\nxMg73hG3jCxZYoKkRbn7bjFXC8JUYMOGuH9Mfb2ZKffYY8Zp3e2e2LYJ44eIFGFC6eqCvXvjzqvR\n/YaGweHfMzLgG9+ICxIZphGE6UPyzKL9++NB7ubMiXdSomm4Q7lCeiMiRRgXkk2z//u/JgZDoiDZ\ns8eEh49SWWkeLuvXmwfM178OO3bEzy9bBu9//zm/FUEQ0oDSUquP2dq1ptOS+Dz51a/goYesTvFR\nwbJkiQmqeN990NYmQ0aTBYk4G0Eizp494bARHfv2mTVkGhoGl7HZTMCz6EMjmhYuHGyylUXvBEGI\nMtznQTS8QKJ42bPH+KYlrhAOMHMmfPKT8VhHlZVmHSRh5EjEWSFt6OkxkVYTo1VGt319qa8pLoY/\n/tGE287MHN7fkaEaQRCiDPd5kJVl1sFatcqaPzBgOkiJ1tv2dvjsZ+Prablc5hkVFS2J6x1Ffd3E\ngffcIiJFGET0R3j0qJlNc+ONZn2aqBhJXJNm5kzzAz7/fLjppviP+6abrCuvzp8Py5ef+3sRBEEA\n49M2Z45VpLzpTfD88/H1txI7XY8+avKjlJWZZ93evdDcbPLq680swsSAd8LYIiJlGqO16Q0cOGCc\n0KLbZ56xrtT7qU/Fxcd73mPtZUTXpElGVugVBCHdSBWczmYzYQsqKuAtb7GW7+qKr2QeTSdPWsu8\n9JJ5Fs6fPzhVVMjw0dkiImWKkWyKfPxx4yuSKEKi+wcOmKEbMMu1V1SYH1by+jOVlWZMdyTIUI1w\ntrS0tLBhwwaOHz9OaWkpmzZtovgs7OpjXZ8w+Rjpc8nrhTVrTIqSvEjo3LnwtreZ5+pvf2usK9FV\npl0ucz5ZvOTlwUc/aiwyMmR0ekSkTAH8fjhyxDiq/uM/xteYqa83JspoXJGoEJk3z3jG33ST2Z8/\n3/yQor4islKvMBqOHTvGtddey/HjxykoKOCee+4hMzOT3t5e+vr66OvrG9H+3r176Ys4OdXX1zNn\nzhzmzp2L0+nE5XINSmfK/+EPf8iRI0di9b35zW/m0UcfpaioiMLCQtxuN0qpifwIhUnAmVaZDgbh\n8OF4hzCannwSDh0a7LxbX2/8Z+6+2zyHq6pM8njO6W2lLTK7J0K6zO5J5ZRVVGTyGxrMFzp529gY\nFyLJFBTA975nxMjcucap7EzIrJrpw5msC6FQiLa2NlpaWmhubj7ttrW19bR/SylFVlYWWVlZZGdn\nn3H/pz/9KZ2dnbHrc3Nzue222wgEAvj9fksaTl5jYyOh5DdEAhkZGTHBUlhYeMb9YDDIDTfcIJYZ\nYdj4/UaoXHZZ3K8FzMrUdrtZhT1KcbF5ZkdTVVV83+GAd70rvZ7R4zW7R0RKhIkUKVqbeftHjhgn\n1X374ueysowFJNFHpLAw/oVN3FZVwc03Wx1W162TYRdhMFprTp48yZVXXsmOhGA0RUVFrF692iI8\nwkkKOCcnh5kzZ1JSUmLZfv3rX6etrS1WrqKigi1btsRER0ZGxogsFRdffDEvJpj01q1bx1/P4suc\nXN95553Hd7/7XVpbW2lra6OtrW3I/ba2NoJRG/4QlJSUcNddd1FVVUVlZSVVVVXMmDFDrDPCIJKt\n1evWGQfelhbT8UzshEZTohOvUtZgl7Nnwxe/aByD58wx8WGcznN3PyAiZdwZC5Ey1NS0gQFj7Thy\nxKTDhwfvJyroRHJy4HOfs4oRr3foNogVRIgSDodpamri4MGDHDx4kAMHDli2iVaKKBkZGVx99dUp\nRUh0mzWEOW6sRcWJEydYv379mFkqzqY+rTWdnZ0W4fLe977XIspcLheZmZn4fL5YnsfjoaqqKpai\n4iWavAk/ZvGZmT6M5jnd32/eF/X1Zqg+0YHXZrNa0202M0w/e3ZcuET3o9ve3rGdSi0iZZwZC5GS\nrI49HiMompsHh3iP/13zZYp+kaLT4aJUVcE//VO8bKo01LlzdU06tGGqXTPUueQX2c9+9jP6+vos\n4iO6X19fT39E/SqlKC8vp7q6mnnz5lFdXU11dTX33nsvddEFkjg7YTHWoiLdSSXKXnjhBdrb22lo\naIilQ4cOWfb7E3okM2bMiImXzZs3cyzhx3+2Ik+YuqSyxDz7rLXzG+0AR7eNjXGHXjDDS4mjn7Nn\nmw7xrFlmXbRZs0wnebiGQBEp48xYiJTqamvYZpvN+JNEP2KtU6fEc6GQCYimtflyZGSkvnao+oTp\nQhjQg5JSKiZobDYVSTZsNmU5FxdCYbq7uwmHQ9hsitzcHGw226QVcefymoGBfjZv3kx/fy+ZmZlc\nfPE6srIyT3uN1pqBgT66u7vp6emiq6uL7u4uurp8NDc3o3Uo4f8JpaUl5OfnU1CQz4wZBcyYkU92\ndnbk/zn1PtOpds14teHUKbjjDrPUSHExfOtbxg3gdNeEw6bDfPSoES2f+xykMKZacLvjgiW6je5n\nZsJnPmMsOqWl8LnP1fKWt4hIGTfGw5Iykf4gZxJDIxE9cs25uUZrCAQCNDUdp7HxWCwdO9bEyZNt\ngBqU3G4v73vfeyksLCYvLx+lbJPiXs/lZzrRbRjuNXv27KG7uyfyvwWnMwOPx0NfXz8DA/5IeYXd\nbicjIzOWXK4MXK4MHA4HWiu0hmAwSFPTcYLBMA6Hg5kzS7DZ7Gn9fxAmNytX1lJXJ2Hx05pUU9Mm\nikQVLUw8ycM0v/jFL+jq6mLXrl3s3LmTXbt2sWvXLvbv3x+bgVJRUcGyZct461uXs2zZ5Sxbtow7\n7riDzZs3x+pdtWodX//6AxN1W8IYcuJEIevX355yuCwYDHLw4EH27NnDnj172L17N3v27GHfvn0M\nDAwAkJeXx5IlS1i6dCl//OMf6e9viFwL5eWTY+goHcTiVL4m8Vx7O3z5y2abn2+mQOfmDr++u+4y\nlpwoCe5ZY4pYUiKkyxRkYerR1tbG5Zdfzu7du2N5Simiv73CwkKWL1/OsmXLYmnp0qXkplhnfrr5\nfQinJxQKUV9fbxEue/bsYdu2bZZyWVlZ3HXXXaxevZpVq1Yxd+5cbDbbBLVamAokjxyMlyVFREoE\nESnC2aK1prGxkdraWrZt28a2bduora2lMXGxkAgzZszgZz/7GcuXLxeRIYw569at46WEWAT5+flk\nZWXR1NQEgNfrZeXKlaxatSomXJYuXUpG1AlOEM5A8gyl8fJJkeEeQRgF4XCYAwcODBIkJyPzAgsL\nC1m9ejXvfve7qamp4f7776e2Nv67XbRoEVdeeeVENV+Y4jzxxBMpLW4tLS3U1dWxbds2tm/fzh/+\n8Ae++c1vorXG4XCwZMkSi3ApLS3lfe97n1juhEEkLzFQO2ayxIpYUiKIJUVIJupH0tTURG5uLrfe\neisHDx6ktraWuro6uru7AeM7snr1ampqali9ejWrV6+mvLycxCBeMkwjpCs9PT3s2LGD7du3x8TL\nzp07LVOlo8i0aGEoZAryOCMiRQDw+/3s2rWLrVu38pnPfIaWlhbL+QULFljEyOrVqyksLJyg1grC\n+BAMBnn99de54oorOHHiRCx/7ty5HIwuDiYICYyXSEmb4R6l1B3AJ4ASoA64U2v96hBlrwU+BcwD\nnMB+4EGt9Y+Syn0BeD+QB7wI/KPW+sC43YQwqUgUJNG0Y8cO/H4/NpsNe9Ia65WVlbz++usT1FpB\nOHc4HA6WLl3K/PnzLSKltLR0AlslTEfSQqQopa4HHgRuB14B7gJ+r5RaoLVONbHpJPAlYB/gB94O\nPKKUatFa/yFS5/8DPgLcAjREyv9eKbVYa+0f73sSJpbkKb8//elPaW1t5bXXXkspSJYsWcKaNWu4\n+eabWbNmDStXruSqq66yRBSdNWvWBN6RIJx7Nm3aNGiYUhDOJWkx3KOU2gK8rLX+WORYAUeBh7XW\n9w+zjq3A01rreyLHTcADWuuvRo5zgBbgFq31z1NcL8M9UwS/38/5559vCfcenfJrs9lYvHgx5513\nHmvWrIkJErfbPage8SMRBEEYHlN2uEcp5QTWAPdF87TWWin1R+DCYdZxJbAA+EvkuAozbPSnhDp9\nSqmXI3UOEinC5KS/v5+dO3dSW1sbs5Ds3LmTQCBgKVdQUMBTTz01pCBJRXFxsTgJCoIgTCATLlKA\nQsCOsXIk0gIsHOqiiGXkGJABBIEPa62fi5wuAfQQdZaMQZuFCaCvr4+6urqYIKmtrWXXrl0Eg0Hs\ndjtLly6lpqaG2267je985zvs3Lkzdu2iRYu46KKLJrD1giCcDS0tLaxfv57m5maxbE4j0kGkDIUC\nTjcW1QWsBDzAlcBXlVL1Wuvnz6JOYYKJ+pIcO3YMr9fLDTfcwBtvvEFtbS179uwhFArhdDpZtmwZ\na9as4fbbb6empoYVK1aQlZUVq+dd73qXjKULwhRiw4YNsQB19fX1nH/++Xz7299m9erVIlamMBPu\nkxIZ7ukFNmitn0rIfxTI1VpfO8x6vguUa63/LjLccxBYpbXekVDm/4BtWuu7UlxfA2y99NJLB4Uj\n37hxIxs3bhzxvQnDo6Wlhe3bt7N9+3b+4z/+g7aERSCUUqxZs4aampqYD8myZcskMqYgTDOqq6up\nT1hmPnFpidLSUktYgGjofyULmI0Ljz32GI899pglr7Ozk+effx6mYpyUIRxnj2AcZ4e1eppS6n+B\nKq31FZHjoRxnb9Za/yLF9eI4O85Eo7RGBUk0cFRzczNgQnUHAgFLEKmqqirLg0mY+vj9LezatQG/\n/zguVynLlm3C5ZKe8nTn4osvtsy2u+iii/jBD34Qi/gcTdHnSU5OjiV67urVq5kxYwY33HCDOMOP\nA1PWcTbCQ8D3IzN0olOQs4FHAZRSPwAatdafihz/C/AaxlqSAVwN3Ah8KKHOrwH/ppQ6ABwCvgg0\nAk+O/+1MXxKHa3Jycrjpppuor69n+/bt7Nixg56eHgDKy8tZtWoV73//+1m1ahWrVq2iqqqKSy+9\n1PIgKisrm6hbEc4RWmvC4T6CQR+hUCe7d19PT4+ZmdXfX09d3d+wcOF3UcqFUk5sNudp9p0oNXjh\nPBE+k59U06GLi4uprq7muuuui5Vrbm6OdYK2bdvGb37zG77+9a8DVutLfX09V1xxBc8995wIlTQm\nLSwpAEqpDwN3AzOB7Zhgbq9Fzj0HHNJavzdy/EXgXUA50IeJl/I1rfUvk+r8HCb2Sh7wAnDHUMHc\nxJIyOsLhMPX19ezatYudO3fy8MMPW4ZrAJYuXRoTItE0VJRWmfY7uTAv/2sZGGjE4Shgzpx7sNns\nBIOdhEI+gsHOiPjwpciLn9M6OIatsmOzuSKCxYnN5iIQOElieCSXq4Sysg/jcpXgcpXicpWQkVGK\n01mMzeYcw7YI6UBXVxd1dXVcc801sfW1EikuLmbFihUsX748tl2yZInFz004PRIWf5wRkXJ6tNY0\nNzfHxEh0u2fPHnp7ewEzzbe3t1eGa6YY4bCfgYGj9PU10N9/iP7++Lar6zW0DqS8TiknDkcudnsu\nDkcOdnsODkd0f+i8/fs/Rk/P9lg9Hk8NixY9gtYBwuEAWvtHvN/Y+FWCwY6EtrlwOPIJBE5g9aVX\nOJ2FFvGSLGRcrhLAzr59t4plZpKRPGRUU1PDv/7rv7Jz50527NjBzp07qa+vj8VUmj9/vkW4LF++\nnOzsbN75zndKRyoJESnjjIiUwTNrbrrpJo4cORITJdEeSFZWFkuXLmXZsmUsX748ti0pKeGSSy6x\nPARkQbL0ItWwh8NRwMBA4yABEt0ODBwj/iJXZGSUk5lZRWZmJSdP/ppgsD1Wf0ZGBWvWvIbdnoPd\nnjnKNp5g1671YyoAamsvxueLfy9zctZRU/NXwuEggUArfv9x/P7m2HZgwHrs9x8nHO4bsn6Xq5y5\nc+8jO3sx2dmLcDg8Z9VeYXwYjqW2u7ub3bt3s3PnTot4iT7/bDYb4XA4Vn7x4sU8/fTTVFZWYrMN\nHmqcLohIGWemm0iJWkb27t3Lvn372LdvHz/84Q/p6OiwlFu8ePEgMVJVVTVoXZsoMlyTvgwMNFNX\ndyW9vXtieUplRIZaQrE8l6uMzMzKmBDJyqqKHFeSkVGBzeaKlR3q5Z9unK3w0VoTCnXFBMvu3e+K\nWGGi2En8DDMyyiOCxSS324gXp7NYZpxMQqLPyx07dvDud7+bU6dODSqTmZnJwoULWbJkCYsXL46l\n+fPn43K5UtQ6tZjqjrPCWZK8Vk1UHAQCAQ4ePBgTIomixOfzAWYxsfnz5+P3W5c0qqqqYs+ePan+\n3JBIlNaJJxwO0Nu7j+7uOnp66ujuNsn6UjXYbFlUV38lJkoyMmaPyAKybNmmQS//dMTlKj4r8aSU\nwuHIweHIITt7AVlZ8y2fZ07OWlaseIbe3tfp7d0bS+3tz9LU9F8xnxuHIz9BvCyKiJfF2GyZ7N79\nLhk+SlOUUpSWllJaWsrixYst1uLzzjuPL33pS+zduzeWnn322ZjlxW63M2/evJhoKSsr43/+53/o\n7Oxk1qxZ0pE7A2JJiTCZLSlaa9auXcsrr7wSyysoKKC4uJgDBw4QDJoHZG5uLosXL2bRokWxtHjx\nYqqqqnA6nYPGa2WoJv3x+9ssQqSnZwc9PXtiTqKZmZW43SvxeEw6fPg+uru3xq5PV8tHujMSy0w4\n7Kev72BEuOyjt3cvPT1mPxzuiZSyxpl0u1dQU/PyqIfMhPFjuNbi1tZW9u7dy549eywCprGx0VIu\nJyeH6667jvnz58fSvHnzyM7OPle3NCbIcM84k+4iJWpuPHDgQMoUtYpEycrK4gMf+IBFlMycOfO0\npmYZqklPjB/Jevr7j2C3Z1NQcBV9fQfo7q7D728CjEXE7V6Gx7MyQZSswOHITapr7P09hNGhdZiB\ngWP09u5lz553EwxaZ50o5cDtXo7Xex5e75vwet+E271UZh9Ncqqqqjh06FDs2O12s2jRIvbv3295\njs+aNcsiXKLJ4/Hwnve8J+2e0yJSxpmJFinRdSkaGxvJycnh1ltv5cSJExYhEp1Fo5SioqKCefPm\nxdL3vvc99u3bF6tPrCCTG7+/hc7Ozfh8m2lq+m9CofjDSykX+fl/E7OOuN0ryc6ej1Kp/YSE9CfZ\nt8ftXklZ2Yfo6nqVrq5X6enZDYSx2TLxeFbFRIvXex7Z2QtTxoYR0pOhLNZaa1pbW9m/f3/KFI0x\nlczs2bP5yle+QlVVFVVVVRQXT4zfk4iUceZciJRQKMTx48c5dOjQoPTiiy9apu4CzJ071yJEoqmq\nqorMTKsZWKwgk5dw2E93dx0+32Z8vi34fJvp7z8EGCfWYLAzYVgAMjPnsnbtwQlqrTAenMnCFQr1\n0N29HZ/v1Zhw6evbD4Dd7sXrXRMTLZmZlRw48HGxlqUpo3lWRy3p+/fvZ8OGDZZYVMmzjbKzs6ms\nrIyJlsRUWVlJXl7ekD6MZ4OIlHHmbEVKoiUkNzeXD37wg7S3t1uEyJEjRwgE4jElioqKqKyspLKy\nkmeffZbOzs7YOYkvMnUZGGiKCZLOzs10d28lHO5HKRdebw05OReSk7OWnJwLycgoZ9u2SybFDBrh\n3BIIdNDdvdUiXAYGjg4q53Yvp6bmFfFvmSKkssQ8/fTTHDp0iIaGBkuK5kWt8AB5eXkEAgGLZWbh\nwoU8+uijVFRUUFJSMuTszdMhImWcOVuRkvzFGeHfJvn/YGYTOFBKxUx30f3k4+GeG4s6pP7hndN6\nAJ9vM+FwPzabg8zMCkKhTgKBk4TDvSgFdnt2JHBYEU5nES5XATabY1Cd4XAfHR3PEQ73RXxS/ga7\nPXvKfnZS/+jrDwbb2bfvVoLBUygFkVMY/5bFeL2r8HhW4vWuxOnMn1T3NlH1pxsjtcREh5ESxcuX\nv/xlurq6UpZ3OByUlZVRUVExZAqHw1x33XWWNjQ2NqaXSFFKVQO3AdXAx7TWJ5RSfwcc0VrvHqsG\nnivUWYqU5BU6i4qKeOCBB9BaxxIw5H5nZyff/va36ezsJCcnh9tvvx2Px2Mpd7o6zlR/OtQxGdp4\ntnWEwwECgXZ6enYTDg8Q/XlpTSSiai52ew42mxelXJP2PidbGwXhbJjMIivVuaNHj1rcC5xOp8XK\nP5zPI/G3tW7dOh5++OH0iZOilLoMeAZ4EbgU+DRwAlgJvA+4buirpyalpaUWkbJgwQJuueWWEdXx\nyU9+cqybJYwzAwPNdHa+QGfnC3R0PE9Pzw7MVFKrI2NmZhVr18rw3USSrkJqLOpIPPb729i//2P4\n/a04HIXMn/81HI6CiIgOMzDQRFfXNrq7t9PdvY3e3oNoDQ7HjIgj9irc7hU4HDOor/80gUAbTmch\nVVX34XDkp+1nkC51TIY29vT08OSTT9Ld3Y3b7eYf/uEfyMrKGnYdP/rRjywzkY4fP854Mdpgbl8B\n/k1r/ZBSKtFm9Bxw59k3a/KRaoVOYWqhtaa/vyEmSDo7X4g5L2ZmziUv71LKyz9Kbu6l7Nt3q8WP\nxOWS1ZwnmsSe5FRn5crTdWRXAW+LHQUC7fh8L9HZ+SKdnX/F5/smWg8QCtmYMyfqkNmEx/Nl8YWa\nQnzzm98c9bV1dXUW94bS0tKxaFJKRitSlgPvTpF/Apgx+uZMXiTS6tTAurZNCVVV99HbuzsiTF7A\n7z8GGGfE/Py3UFX1RXJzLyEjwypCJkskVkFwOvOZMeNqZsy4GoBweICurq3s3Pl2gsF4+He/f/x6\ny8LkIlWnPDlI3VgxWpHSAZQCDUn5q4FjZ9UiQZggwuEAdXV/Gxmygf7+eurqLkcpBx7PGmbOfDe5\nuZeQm7sOp7PgtHWdbRh2QZgobLYMcnMvIjt7cZI1cPx6y8LkIlWnPN1Eyk+Bf1dKvZPIALxSah3w\nH8APxqpxgjCehEL9dHW9Ehm6eZ7Ozpcs8UjAPJgvuGA/drt7glopCBODWAOFdGC0IuVTwDeBo5jl\nP/dEtj8BvjQ2TROEsSUY7MLneynmT+LzvYzWfuz2XHJzL6ay8rO0tPw4ZkkB42siAkWYjog1UEgH\nRiVStFm97ANKqS8CywAPsE1rvX8sGycIw8XqS2J6fUrZ6ez8a8xS0tW1DQjhdBaTm3sJ1dUPkJt7\nKR7PcqIh5UtKbpXeoyAIQpowWksKAFrrI8CRMWqLIIyaXbs2xMbP+/vr2bKlinDYRFnMyJhNXt6l\nlJbeTl7epWRlLRhyhof0Hs89/hY/uzbswn/cj6vUxbJNy3AVu9KuTkEQzj3DFilKqYeGW1Zr/fHR\nNUcQhofWIbq7d+LzvUhn54v4fFss55Wys2jRD8nLu4TMzDkT1Eohig5rAqcCBFoC+Jv9JrWYbfOj\nzQRaTSCp/vp+Xl74Mrlrc8EGyqbAHtlGjpVdnfFc2xNt+I/7Y3Vuv2w7Cx9ZSEZ5Bq4SFzaHLMgn\nCJOBkVhSVicdr8H4obweOV4AhICtY9AuQbAQDHbT1fVyJJbDi/h8mwmFulDKicdTg8s1E7+/KVbe\n7V5BScmNE9jiqUmyhWLh9xZCGPzN/rgAiYiPxP3AiQA6aI3+anPbcJW4CPqClnwd0NiybRA24kYH\ntdmGdSyPEJZjHUo4F4ZAmzV6Zu++XrZduC3yh8FV6iKjPIOMWRlmG02z4ltbxmAhIxaaicPvb2Hn\nzr8nEDiJy1UmCydOE4YtUrTWb47uK6U+DnQBt2it2yN5+cAjwAtj3UhhapPKn0TrQCy4VGfni3R3\n1wEhHI48cnLWMXv2v5Kbuw6v903Y7VkpV5EVxobAqQDdO7rp2dHD4XsPEzgRt3q8uvBVS1lbphEe\nzplOXCUuvOd5cZW4TJoZ3zpnOnF4zOOn9uJafC/Go1d6VnlY9viys2pzcp3eC7ws/M5CBhoHTDo2\nENtv/2M7A40DhHwhSx3OIucgIdP8/Wb69vfF7n/X+l3U/HV8Vk0XrOzatYGurtcA6O9v4LXXVjNv\n3kORhTgrpkWQvunIqNbuUUodA/5WJ63Ro5RaBjyrtZ504TXPdu0eYfTU1q7D53spdqxUBloPAJCZ\nWU1u7rpYys5ejFJiqh8PwsEwfW/0GUFS1xMTJgON5n+hMhRo0P74M8NZ7GTpL5bGhIjdax/xy8J/\nws+u9WPskzKKOoNdQYt4GWgcwH/MbzlOttDYsm3Mvns2ntUePDUeMmZlyMtynNiypZr+/sSlJRyA\nscP7zvAAACAASURBVMK5XKWR1cPNCuJe7xrs9qwJaed0ZbxWQR6t42wOUJQivwjwjr45wlQnGlre\n53uFrq6X8flesQgUAJstk0WLfkxOzjoyMkomqKVTG3+bn54dPXTXGSHSvaObnt096AEjQDLKM3Cv\ncDPzxpm4V7rxrPCQtSCL7Zdvt1gosuZnkXdp3lm1xVXsGnNrxGjqdHgdOBY5cC8aesp57UW1+DbH\n71+5FI3faCR40rwsnUVOPDUevKu9ZlvjJXNupgiXMcDlKrWIlJycC1i27HF8vi10dm7G59vCoUP3\nRFYZd+DxrLIIl8zMSgKBE4OstjJklN6M1pLyA+AS4J+BVzAB3dYCDwAvaK1HtrJeGiCWlPEhEDhF\nV9er+Hwv4/O9TFfXKwQCbYCJQZKTcz4+38v098eDF+fkrJMZNmOAv8VYE/qP9GPPtpP/lnz6D/bT\nXdcdcyq1ZdpwL3PjXuHGs9Jjtss9OGc4U9c5DlaPyUSq+3cWORloHKC7tpuubV1mW9uF/5j5jO05\ndjyrjWCJbrMWZmFz2MTHZQSkGtJNFhjhcJCenp34fJsjaQt9fQcAcDpnonXAEupfnjVjx3hZUkYr\nUrIx0WXfC0SfZkHgf4FPaq17hro2XRGRMnriPiVN2O1eioqup7d3D11dr8QW4HM4CsjJOR+v93xy\nci7A630TLldR5PozP3yE4RH2h+l6rYuO/+vgyANHCHXE/SyUS1HwtwUxy4h7hZvs+dlmRoww5vhP\n+OneZgRLdNt/sB8AW5YN9wo3A4cH8Df7Y9fkrMsRH5cxxu9vjXSSNtPY+BDhcH/snFIOCgs34PWu\nwes9D6+3BocjdwJbO3lJK5ESu1gpN1ANKODAZBQnUUSkjIxQqJfu7h10d2/j0KHPEwi0JJxVMTFi\nhMkFZGVVi8l7HAgPhPG94qPj/zro+EsHvpd8hPvC2L12dFAT7gvHymbOzWTtwbUT2Foh2Bmke3tE\nuNR2c+JnJ9CB+DPY5rYx78F55L05j6z5WfKbGWNqay8etB5RVlY1XV3bYktiZGXNw+s9D49nTUS8\niHAZDunmkwJARJTsOGNBYVITCJykq2sb3d3x1Nv7BhBGKQdgdWTNzKxkzZotKesSzo5QXwjfFh8d\nf+mg8y+d+Lb4CPeHsefaybskj8ovVJJ3WR6e1Z5B/iOuUhlGmGgcuQ7yLssj7zLjx9PX0Gf5Hymb\n4o073oAQuGa5yH9zPnlX5BnRUimOoGdLqvWIXK5itA7R2/s6XV1bI+k12tqeigWEzMqaj9e7hqys\nhbS2Pk443I3LNUusvueA0Q73/Bnjh5ISrfUVZ9OoiWA6W1ISpwA7nQXMmvUR+voaYoJkYOAoADab\nG49nJR7Parze1Xg8q3G7l7J9+5WW3omM846OVP4Jdo+dzpc66fxLp7GUvOxD+zWOfAe5l+bGXnie\nlZ5BwzbT3X9kMpDqf2TLstH5Qicdf+6g/bl2urd1g4bMqkzy3pxH/hX55L05j4yyjIlu/pQmLlxe\ni4kXn28zELdOOhz5lJV9CLd7BR7PCrKyFmCznVXff9KSVsM9SqmvJmU5gVWYdXy+r7X+2CjqvAP4\nBFAC1AF3aq1fHaLs+4GbI38PTAC5TyWWV0o9AiQ78P5Oa/22IeqcViIlGOykp2c3PT27OXToHvz+\n45bzTmchHs/qWPJ6V5OVNS+2xk0i4lMyNiTH9rB77YT7w+iAxlnotIgS93K3ibQqTHkC7QE6/tJB\nx59N6tkZGZZYmBW3tFyeB2HECXec2bJlrsXJ32bLwuEowO8/BoBSLtzuJbjdy2PCxe1ejstVglIq\nZUyoqfKsTKvhHq31XanylVKfwyw2OCKUUtcDDwK3Y2YL3QX8Xim1QGvdluKSyzArLr8E9AP/Ajyr\nlFqitU582z4D3IrxmQEYGGnbJjvBYBe9vXtigqSnZxc9PbtjPyqwDRIeGRmzWbv20LDHw2W9m7Oj\nd38vp545RXdttyVfhzXzvjqP3MtycS8RUTJdceY7KbqmiKJrIo7mJ/xGtDxnLC1N/20iLduybYR7\nTS9fAs2NDy5XmUWkeDw11NT8lUDgFD09O+nu3klPzw56enbS2rop5ufidBbidi+nt3dfrEPY31/P\nrl3r5dl5Bs7KcXZQZUrNA17RWheM8LotwMtRC4wyb8ejwMNa6/uHcb0NaAfu0Fr/KJL3CJCrtV4/\nzDZMSktKXJkfw27PpbT0fQwMHI2JkoGBw5GSiszMubjdSxPSMrKyFlJX9zcyXHMOCfWF6Pi/Dk79\n9hQnnzlJ/8F+lFNhy7YR6ozPxpGZHsJwGDg2QPuf23njw28Q7ooPRdgybVTdV8WMt80ga4E44Y4F\nI7Eaax2mv7/BIlza2p5E68SAgHby8i4jO3uRJWVklE+6/1daDfcMWZlSNwH/PpKIs0opJ9ALbNBa\nP5WQ/yhGZFw7jDq8QAtwndb6t5G8R4B3AAGMgHkO+Det9akh6kh7kRIOB+nvP0Rf3xv09r5BX99+\nWlp+TCjUaSmXkTEnJkKigiQ7ezF2e3bKemW4Zvz5/+ydd3xb1d3/31eStWzL247tLDs7cZpNJiMk\nrDAKCaOhjEIpBVqgodDSPjwthae07NlCKVACgdAAAcqvoQkJhBSy9yAJie0s25H3kC1ZlnV+fxxb\nlmQ7cRLbkuzzfr3u617de+7Rsa909bnfdZqtJeWflVO5uhKvy4tpgImkS5JIvCSR+PPj8dZ5VQyJ\n4rQJdhca4g00OhsR9QJztpmkOUkkzkkk/rx49JbWbltF1xOcXWQy9cNmm0Jd3T7q6r7zVdrW6aKx\nWoe1Ei96vY29e38YlvfqsHL3aJoWPDGKBqQDE4FHT7G7ZOREhfag/XZgWAf7eBwoAFb67fsM+BDI\nR6ZJ/wlYpmnaVNGZyuwMCfZRjhr1AUJ4cDqlCJFiRIoSlysPIWRlS53OjMUy2Pe6GbN5IFOm5Lf1\nVu2i3DWdT2Ndk7Xks0BrSdw5cWT9XxaJcxKxDrcGPi3FoCwnitMmZ2lOK5Grj9FT+WUlZcvKKP20\nlIKXCtBZdMSfHy9FyyWJWLJU1lB30V52EchAXZfrUJNgaVkqKlb4CmD643LlsXXrFAYPfh6LZRBm\nc1aPnArgdANn3yQwu8cLlABfCCFWnGJf6UiBMVUIscFv/xPADCHEtJOc/yAy4Pbc4LmEgtplAbnA\nLCHEl20c7zZLitdbj8t1GJcrn++++xkuV67fUR0t0eN6LJYsLJahWK1DsViG+LalOVDXSpkrV033\n05yVU3+kHs2gYRpgomZ9TZvWkuZJ9RSK7kYIQd2+OulmXFZG1ZoqhEdgHWElcU4iSXOSiJsRh86o\n5sYKN9zuUpzO/ezefWWbgqUZozEDi2VQk2gZhMWS7duOikrq0uDdiHD3nNYAzsDdo2na/cBvkcJj\nWwfeqxj4HyHE39s4Nh7Ycs455xAXF1i4Z/78+cyfP7/NPtu64AZDIvX1x3C5DuFy5fsWpzMfl+sQ\nbnch7WVwGwyJDB++EKt1KGZzFjpd2+XJW95fuWpChdftpXJNJftu2ucrMw+gj9Mz8HcDSbykDWuJ\nQhEmeKo9VKyq8IkWd6EbfYycPsE2w0bJeyU0lDUo12MY0fqhdBqjRi3B6czD6czF5crF6cxt2s4L\nEDR6vQ2LJRuX6ygeT5lvf0zMWMaM+RKDIa7D96rFixezePHigH1VVVWsWbMGwkGkaJqWB0wSQpQF\n7Y8Htgohsk+xv7YCZ48gA2efbOecB5AC5cL2UpWD2vcFDgPfF0L8vzaOn5YlJfhDo2kmoLGVG8Yf\ngyEBszkbiyWLqqq1TaJFYjRmkpZ2Q/OYmnv1W2snOEY7x090rL1zO9Jv5517an9PZ5x7ev9HT5WH\nmq011GyqoXZ7razoqmnS+CVk26g0I8P+Plxdg5Oe212f5TM5t/uuQSjFrBACxw6HT7D4x7aACuIO\nF071odTjqfITMHJtt78VMDVAM3p9LGbzAEym/pjN/VutjcZMdDpDu5aYsLKkaJrmBfoIIYqD9qcB\nR4QQp1RlSNO0a4GFwE9pSUG+GhguhChpmtDwmBDit03tfwU8AsxHpiE34xBC1DaV6/89MiblODAY\nGbcSDXxPBIZXN4/htERK8PThmhaFxTK06ZVoc93yPxcI4cHtLkQID5pmICoqramK68nPbVmf6FhX\nnNuRfjvnXIWid9LVYu/E53oqG6El0Qx0OqISDQHn9ESh2L2C+0zOPf1rUFm5JsCSomnGpofqlsyw\nttFhMmXi8VTS2Fjj29scYhAWgbOapl3h9/IiTdP800r0wCzg0KkOQgixRNO0ZKTwSAO2AxcJIUqa\nmvRFTmDYzJ3IAnIfBHX1h6Y+GoHvIQu+xQOFwHLgd20JlDMhePrw2NizVExIJ9OekOmYODr1cxud\nHqrWVFG+vIzy5eW4i+rRReuJnxVH4oUJJMxKICrF6DvHXVrPvlv20VDsJirNwPA3hvvNIty7RWbz\nNajwNPBQXh4lDW4SDQbu7dsXq06jQUCDaKTBK2gQXhqEF4/XS4MQuP22ffvx4vG1FXi8XtZVV1HV\nKG8PGoJMo5GfZ2aSHGUgyWAgSqfr8Ji77xqcybldf/0K/n2M+kP1oAnQBKYBJjJ+mtGhc8/0+3d6\n/8czObfnXb9mwdF8DWQzeSwm5ns4HNvweuvRNBMxMWPQ6Uzt9N26/+rqwClPgguBdjanZElpsqCA\nHHGwbbIBKVB+2ZY7Jdw5XUuKigmJXPzL0EclRpFybQpVa6qoWFWB1+mVaZuXJ5F0WRLx58SrgMIO\n4mxs5IDTyd66OvY1LZ+WllLrPdmTWtsYNQ2jTodJ0zDpdAHbJp2OnQ4H9Se4jyVHRdHXZKKvyUSm\n0ejb9l9iDK2f1+xuN/N276bI7SbdaGRpTg6pxt4Rl6GmVFC0R3vJGuHm7slHxqS0H2YcYURCnRRF\n5+H1eNkyfouvxHgzcefGkXSZFCbWYSro9USUuN0+EbKvrs4nSg65XL7nrtSoKIZbrWxzOKhpbPEf\npBuNLBk5UgqOdgSIUafDqGknvQYztm7lm+qWGIopNhsLhw/nWH19u0tJQ6BBNU6vbxEyTet37HZy\nXS2+++k2G1+re4Oil9Peg3lYuHuaEUJkddYAFIruwm13U/4fWbekYkUFnorA4GbTQBPjVo8L0ejC\ng2DrwfujRlHn9bK3tjZAkOyrq6PM01SzB8i2WBhutXJ1SgrDrVaGW60Ms1pJipJur2AhkW02MyM+\nvlPGvDQnh7ltWDyGWtsuXghQ7/VS2I6A2V1by3/Kyyl0uwPO2e5w8OzRo0yLi2NcTAxGnbKsKXof\n3V1Xq8MiRdO0e4BXhRCupu12EUK8cMYjUyjOENEoqN5YLQuqLSvDsUXOjRM7KZbMuzMp/biU2p0t\nlhRTZu+eVVYIwSU7d7LNIf9PeS4XmevW+awi0TqdT4BcnJjIiOhohlutDLZYMJ3kB7stIdFZpBqN\np2zhMOl0ZFksZFnaL341fetW1voJK72m8dv8fFxeL2adjkmxsUyz2ZgeF8dUm43kXuIKUii6kw67\ne5pcPBOFEGVN2+0hTjUFORxQ7p6egbvETfnycsqXlVO+vBxPuQdDgoHEixJJnJNI4kWJPt+68rvD\nUZeLlRUVrKqoYGVFBfYgN0iSwcDikSMZbrXS12TqVe6vYre7lbCKNxjY5nCwtqqKtdXVfFNVRVGT\nxWWoxcL0uDim2WxMi4tjuNWKrhf9vxS9m7CKSemJKJESWfgHvepj9CRcmEDVmipqNtWAgJjxMb6y\n37FnxaIzKNM8QEVDA6srK1nZJEq+czrRgPExMcxOSOD/lZWxp67O117FYZwYIQSHXS6fYFlbXc1O\nhwMvkGAwMLVJsEy32Zhks+FobOy1wbiKnk1YxaRomvY74CkhRF3QfgvwgBDikc4YnELRFq7DLnbM\n3oHzoNO3r3ZPLSnzUsi4M4PEixMx9endrptmXI2NrKuu9omSzTU1eIHBFguzExJ4LDubmfHxJDbF\njtzXr1+XuWV6IpqmMdBiYaDFwvVpaQDUeDxsqK5mbXU1a6uqeOLIEaobG9EDZp3Ol+WU53Ixd/du\nJQI7iP+DSW+1fPZGTje7pxFIb6OYWxJQLISIuCk2lSUlfGkoa6DiywoqVsrFldu6WqI5y8yUvCkh\nGF144RWC7Q6HT5R8XVWF0+slJSqKWQkJzE5IYFZ8PANPEIuh6FwaheDb2lrWVldz/8GDOPxSsa06\nHS8OGcLlSUmkKIvKCQme5VlVwQ0vwsqSgqyR0pa6GQOUn/5wFAo5g3DV11VSlKyqwLHNAQIswywk\nXpxIwqwEjjx+hJoNLVUPjRm98wZvd7u5dOdODtfXI4TAKwQVjY1YdTrOiY/n0awsZickMDo6WsVH\nhAi9pjE6JobRMTG8ffx4QJaTXtO4bf9+NGBGXBxXJSdzZXKyEpFt4D8/FoBju4PC1wpJmJmAOdvc\nq+KlehOnWnG2gpayd99pmuYvVPRADPBK5w1P0VMJMN32MdL/f/pTu62WipUVVK2tQrgFxj5GEmYn\n0PeevsTPisfc1+w7P256XKug195EntPJkuJi/u/w4YAiaf1MJj4ePZopNptKkQ1D2spyEsCnpaV8\nVFrKr/PyWJCby9iYGK5MTuaq5GRGR0erH2DAmG7ElddiRdU0je9++h14wdTfRPzMeBJmJhA/Mx5z\nf/MJelJEEqdacfZmpBXlDeAXgH9ZfDdwSAixrlNH2E0od0/3IRoFmydspnZHYCE1faxe3mhmJZAw\nOwHrCFVMzZ9DTifvl5SwpKSEzTU1WHQ6dBAgUrLNZnKnKLdXpFLj8fBZeTkfl5by77IyqhsbyTab\nubLJwjItLg59L/1OtJWNpzPpqFxTSeWXcnHskFZXc7Y5QLSYMlSMWlcTVtk9mqadC6zt7HlwQokS\nKV1HY10j1Ruqqfq6iqqvq6heV01jTWNAG2OGkSmHp6gsnCCOuFxSmBQXs7GmBrNOx6WJiVybmsql\nSUlctGNHgPtAZeP0HNxeL19WVvJRSQmflJVx3O0mJSqKK5KSuColhVnx8Zj1ERf+16U0lDdQ+VWL\naKndLR+ELMMs2CbbqN5QjdflxdTXpAJvO5mwEikBHciMnij/fUKI6naahy1KpHQe7hK3T5BUfV2F\nY6sD4RHo4/TETY8jbkYcxYuLA0rSqyC4Fo65XHzQZDFZV12NSdO4JCmJa1NSuCwpiVi/eWbaquWh\nUlp7Hl4h2FBdzcdNbqEDTicxej0z4+PZX1eHW8iJFdX1D8Rd7Jai5YtKji88jtfZYnWMSoli4MMD\nsU2zEZ0TrR6QzpCwEimaplmBJ4BrgaTg4yq7p+fjH1NiiDeQdnMatTtrqfq6Cud+mRps6m8ibkac\nb4keFY2mk6ZqVUgtkML6eilMiov5proao6ZxcZPF5PKkJGxtTICn6J0IIdhbV8fHpaX86cgRHH5z\nIo2wWtk+caKKR2qD9YPWB8a0GDXwgvAIdNE6bGfZsE2zETc1DtsUm99M5oqOEG7ZPU8CM4E7gbeB\nnwGZwE+BBztnaIpwpKGyAccWB/t+vI/6w/W+/Y6tDqJHR5MwK4GBvx9I3PS4EwavGVONvdJy4j83\nTlJUFFclJfFZU6qwQdO4MCGBt4YP54rkZOKUMFG0gaZpjIyOZmR0NK8XFQWIlL11dWSsXct1qanc\nkJbGFJtNxXU1ERx4GzspljGfj6FmSw3Va6upXldN0d+LOPLHI4B0EcVNjcM2VYqXqKQo9lyzRz1Y\ndTOnexe8HLhJCLFa07R/AP8VQhzUNO0w8EPgnU4boSJkeGo8OLY5qNlUQ81mufgKqAXd90wDTUza\nOan7BxlhfH/XLjbUyNTpPJeLTTU1XJKYyBvDhvH95GQSotTTWwB2O8ybB0VFkJ4OS5dCamqoRxU2\npBuN5PnN1Dw2JoYLExJ4x27nr4WFDDKb+WFaGjekpTHkBBMu9gZylua0st7qLXriZ8QTP0NOdimE\nwJXvonpdNVVrZfzc8beOgxeZv9qkB115LnZduovxG8crEdjFnK67xwGMEkIc1jTtGDBXCLFR07Qs\nYJcQIqazB9rV9DZ3T3D1xhGLRuAucgcIkrp9dSBAZ9ERMy6G2ImxcpkUy/4f76d6rSqs1BEavF5W\nVFTw1vHjLCkpCTg20GQif+rUEI0sjKmvh8JC+P73Ydeulv3jx8OKFZCQAMql0W5MUqMQrKmsZJHd\nzgclJVQ3NjI5NpYb0tK4LjVVFY47BTwODzWbatgzb0+rmdMNSQZix8cSMz6G2HFybRlk8bm1exPh\nFpOyE7hbCPGVpmkrge1CiPubZkf+tRAis7MG2F30JpHSUNnA9nO2BwSuNqMZNWLGBAoS6whrq6Ay\nFVNyYoQQbHM4eNtu5127neKGBnKio3F4PByqb3GT9bpsHCGgogIKCtpfjh2D0tIT96PTSaGSlCSX\n5OSW7bZeNy/+P869xErjbGzk/5WV8bbdzmfl5QghuDgxkRvS0rgiORmryhDqEMEVb60jrKRcm4Jj\nmwPHVgf1x+T3Wh+rJ2ZsjBQu42OJGReDdYQVT5mnR5f1DzeRsgBoFEK8oGnabOBTQId0H90nhHi+\nswbYXfREkeKt91K3rw7HLge1u2rlsruW+qP1rdpGJUfxveXfk1HuRvWEeroU1Nfzjt3OW8ePs6eu\njtSoKH6YlsZNaWmMiYmhpKGhZ2Xj+P/Q9+kDL74orSD+gsNfgBQWgrNlziU0TQqDzMy2lwcfhG3b\nWtrn5MDDD0NZmRQyZWUti//rigopiIKJjW0RLLm5UFnZcmzyZFi/vsv+VeFAqdvNkpISFtntrKuu\nJkavZ15yMpcmJvJcQQHHe8rnsgs42YOZu8Qt3eNba3zCpdk9rjPr0AwajY6W+KHYybFMWD+h2/+O\nriKsREqrTjRtADABKAVuEELcfsaddjORJlKCK7YOfnYw7iK3FCS7pSCp21/n86Ga+puIzokmenQ0\nMaNjOPrUURzbHb7+lLvm9HF4PHxUWsrbdjsrKyow6XRcmZzMjWlpXJiQgKEnuiWEgG+/hcsvh/z8\nttuYzYGCo2/f1iIkPR1OFIdTXAxz5566taOxUQqQEwmZd9+F2iBr4ujRcO65cM45cmmaNLAnkut0\n8o7dziK7nQP+wpFeaOHrIjxVHhw7pHDJ/998vA5vwPGotCiiR0YTPSoa60gr0SOjsY6yYkyOPIEY\n1iLF15mmjQG2qhTkrsFb78WZ66TuuzoO/uJgQHZNM4Z4A9Gjo31LzOgYonOiMcQFxkgrd82Z0SgE\nqysreev4cT4sKaHW6+WcuDhu6tOHq1NSemZmjt0OK1fKmJCVK6VVJJg+feTxzEzpjgnnoMIZM+Cb\nb1peDxkCZ58Na9bAwYNy37BhLaLl3HOl0OphCCHou24dhe6WuXFU5eLOJ9hdZBlmIfW6VOq+raN2\nTy3OA06ER/4eR6VEtRIuUclRfHfHd2F7zw63FGTFGdLetOPCK6g/Wk/dd3U4v5OCpG6/3HYddsko\nc2iVXRPVJ4qJmydizDB2KNq8t6YAny7NqcOHXS4E8jIUud0MsVj4df/+3JCWRlZPmxSurg7++1/4\n/HO57Nwp948ZA9dfDxdcIF0v6/xmwhg0SFojIoGlS9u30hQUyL/9q6/k8uqrcn9WVqBoycoKbyHW\nATRNI8tsDhAp6crV0+m0lV3kLzK8bi/Og05q99T6hEvlV5UU/b0I0RBoTHDludgyaQuDnhiEZbAF\ny2BLqwfRnoKypDTRnZYU0SjYOmUrNZtbZvE1JBkwpZtwHnTidUklokVp8gM41IJ1qDVgvefqPWra\n8m7ieH0907ZtI98v1bOP0chHo0YxuSfVofB6Yft2KUhWrJBWhvp6yMiQguSCC2D27EAXyOm6YyKN\n4mIpWtaskaJl507p8srMDBQtw4ZFpGhRlYvDF2+DFC/bZ26nwe43E42OlodWpPWlWbAEL1GJUe0+\nGHcWyt3TxZyqSDnRBRdegdvuxnXIhSvfFbg+5MJ12NVKGevMOvrc2gfrsBYxYupvardUs3LXdC0O\nj4ePS0tZZLfzeUUF3qDjPcYcfuRIi6Vk1SoZsxEdDeed1yJMRoyIyB/eLqWiAr7+ukW0bN0q42BS\nU6VgGTsWPvgAqqt7tnhTdBvB7iLbdBuj/zUaZ64T58Gg5YCThpIWQWNIMCAaREDgrnWklVEfjMI8\nwIzeeuY/2WEhUjRNW3qSJvHAub1BpAR/YIx9jESPicZ1yEX94XqfNcQffYwec5YZc5aZmk01uIta\nzKvWUVYGPzs4wI0T8ITu/xvRznab7U+l7Zm078q+/ba7su9Gr2BtTTWflJayvKIcp1cwKTaW76ck\n835xMdtraxFNbSfGxPDB6JwzH0t3XZ+SErjpJigsRLOYYdIkmcmSexA0HYwfB+efD+fPhElngSmq\n4337bXfF9YkIS1VNDaxd2yJa1q4NzC6aOBE2bQrd+BQRz6k+mHqqPAEC5vCfDrcK3G0mKjkKU38T\n5gFmzAPMvu3mdVRyFA3FDSe0xISLSPlHR9oJIW457RGFiFMVKcHzQCgUim4mnIV4TQ0Ir99uAQYD\nGE1NAlBr99yuFOIB7buyb//94Xh9zqB9pF6f6g3VNFa2WFJ0Fh3owFvbtnDxR2eRFn3/CRqDQwzC\nInA2EsVHV9HWPBA5Hzc9WfvpvgAR6K8H29k+lfZd2XdA+67s239/N/Z93O1mVXk5K8orOFxfT4Je\nz6z4BC5ITGSIxYzm/y0/lbGE07U/chRWfQFfrUbU1Pgd0CA1DV56sVde+07tu533ES+9BIcONb3S\nZHE5Wxzk5YHHCN/7HkyYCAMHgqap730k9t3O+4T8e99O33HT4qjeUI3X5UVn0mGbYkNn0nX42les\nqsAff09AV9Izw4G7gZNFaivCh+bMnIL6evSaRkpUFOtrarCm67hqdDJ/Tktjdk+pZ1JbC0uWn5lv\n2AAAIABJREFUwGuvSZdDcjL85Cb48svAomhDpsM1Kkaiy7jxl0EBxS/ImJT8fHj7bXjzL/C3fMjO\nhptvlq64gQNDPWqFol1ahTikd8/vXacGzkYykVAnRXHqVHk8TNqyJaBYVbxezwtDhnBlcjKxPaGe\niRCwZYsUJu++Cw6HDHi97Ta44gowmXpPFk6k4PXKwNs335SisrZWBiv/6Eeygm9MxE1/pujhnCwm\nJixiUroSTdN+BtwP9AF2IOcGajPSTNO024CbgObIxS3Ab4Pba5r2CHAbMqD3G+BOIcTBdvpUIqWH\nUNnQwCdlZbxfXMyKigoagj7jPSYzp6IC3nlHipMdO2ShsVtvhVtuUU/lkYTDIUXjwoXwxRcyu+qa\na6RgOftsNZGiIiLoKpESFp9+TdOuA54Gfg+MQ4qU5ZqmJbdzyrnAu8B5wBTgKLBC07R0vz5/Dfwc\n+ClwFlDb1KfyyfRAKhoa+EdREZfu3Enq2rX8aN8+Kj0enhw0iImxsQFtI7pQlRCwejXccIO0iCxY\nIAuoLVsmYyD+8AclUCKNmBjp7lm1SrqDfv1rmSV03nny2j78MGzcKCvkDhok18XFoR61QtEthIUl\nRdO09cAGIcS9Ta81pPB4QQjxRAfO1wEVwM+EEIua9hUCTwohnm16bQPswM1CiCVt9KEsKRFGWUMD\nn5SW8n5JCSsrKmgUghlxcVyTksLclBQyTSaghxSqOn5cPmm/9pos2T5kiHTn3HSTLEWvwO6wM2/J\nPIocRaTHpLP0uqWkRkeoS0sIWUyv2R0UEPgMTJ8u3UUKRZgQFtk9XYGmaVHIyQkfa94nhBCapq0E\npnawm2ggCihv6jML6TZa5ddntaZpG5r6bCVSFOFLc+BrkdtNclQU16WksLyigi8qK2kUgnPi4nh2\n0CDmpqSQ0SRM/Ek1GiNvsjS7XcaQ5OWB2y0nyzMapRvg9delGyAS6oc00VkCotHbiMvj8i1Oj9O3\nfesnt7KreBcAeRV5zP3nXL6+NUJ/yDVNWkxmzIDnn5fl90tKWo5/+60svJfcnrG55+F/H4jYhw3F\nKRNykQIkA3qklcMfOzCsg308DhQAK5te90EmTrXVp3rsjDAu3bmTLQ45Y3Oey8XGmhpmxsfz/ODB\nzE1Opk8bwiSiqamBqVMDZxfOzpbBsfHxoRvXaSKEYM67c9haJB+u8iryGPe3cVwx9ApcjX6Co8HZ\npvjw39/gbTjJu7WwuXAzdy+7m0mZk5iUMYmhSUPR6yKuzqSMURk6NFCkVFZC//4yM+gXv5Dl+Hs4\n83bv5ptqmV2S53Ixd/fuyHv4UJwy4SBS2kMjMFu77Uaa9iBwLbLS7ckSt0/a54IFC4iLiwvYN3/+\nfObPn3+yoSg6CbfXy9dVVSwrK2NZeTl76+oCjg8wmfhi7NgQja4LKSiAF16Av/0NqqpaH48ggWJ3\n2FmVv4qVeStZmbeSo9VHA46X1JawsXAjZoMZi8GC2WAmxhhDsjXZ99psMGOJatn2b9vWsR9/8mO2\nHm+xMttMNlbkreClTS8BEGOMYUL6BCZlTGJixkQmZU4iKz4rMiraBk+G+Pe/w4cfwksvwSuvwGWX\nwS9/KecPioS/5zQocgfe3rc4HDx79ChzkpIYarFExnXsISxevJjFixcH7Ktq657VCYQ8JqXJ3VMH\nzBNC/Mtv/5tAnBDiqhOcez/wW2CWEGKb3/4sIBcYK4TY6bd/NbBNCLGgjb5UTEoIKayv57PycpaV\nlfF5RQU1jY2kG43MSUzkm+pq9vkJlek2W896gtqxA55+GhYvBqsV7rhDZnls3tzSJsxjEBxuB2sO\nr/GJkma3y+jU0czOns3yg8v5tvRbX/vp/aZ3uiumuLaYuf+c28qlVOmqZEvhFjYXbmZT4SY2FW7i\nSNURABItiVKwZEzyiZdMW2anjqtLqa+XaefPPAO7d8O4cXDffXDttdI92IOYsXWrz5ICspSA0+ul\nXgiyzWYuSUxkTlIS58XHY9VHoMUswunRKcjtBM4eQQbOPtnOOQ8gBcqFbaUqnyBw9iYhxPtttFci\npRtpFIIN1dU+a8k2hwMdMMVmY05SEnMSExkbE4OmaT0j8DUYIeSkfk89Jdf9+0uz/Y9/DDZb2Nc1\naWhsYGPBRlbmrWRV/irWHVuHx+uhn60fs7NnMzt7NudnnU+fGOldbU9AhIri2mIpWgo2sblIru21\n0jucHpPOpMxJTEyX1paJGRNJtiaHd2CuELBypRQr//mPnJ357rvh9tshISHUo+sU2roPROv1rK6s\nZFlZGf8uK+NwfT1mnY6Z8fHMSUzkkqQkBlksoR56r6Cni5RrgYXIdOGNwALgamC4EKJE07S3gGNC\niN82tf8V8AgwH1jr15VDCFHr1+bXwI+AQ8CjwChgVFtuISVSuo7mgLdj9fWYdDpyoqNZXVlJucdD\nksHAxU1PQBclJpIUFXXyDiMZtxvee0+Kk127YPx4uP9+uPpqCOO/XQjBtyXfSktJ/kpWH1qNw+0g\n3hzPzIEzfcJkSOKQiDS7CyE4Vn0swNqyuXAzla5KALLis6iur6bMWeY7pyusQZ3Cnj3w7LOysq3B\nIGvn/OIXMn25ByOEYF9dHcvKy/msrIw1VVU0CMFQi4U5SUlMjY3luYIC7D3pYSeM6NEiBUDTtLuA\nXwFpwHZkMbfNTce+AA4JIW5tep0P9G+jmz8IIR7x6/Nh4HZkMbf/IlOUVTG3bqK8oYH/VlVx13ff\nUejnT47W6VjQrx9zEhM5y2ZDH4E/aqdMZaWMNXnhBSgshEsvleIkjGMIjlUfY1XeKlbmSxfOccdx\njHojM/rPYHaWFCXj08dHZjBqBxBCkFuRK60thZv566a/4mpsma8ryZLEgbsPkGAJU0uF3Q4vvwx/\n+QuUlcGVV0pX0PTpYfuZ60xqPB5WVVSwrMmNXBAU0zI6OprNEyZgVMXyOoUeL1JCjRIpZ06p282a\nqiq+qqxkdWUlu2prEcjobI9fux5T8bUjHDokU0hfe01aUW68Uf5QjBwZ6pEFYHfYufK9K8mvzEen\n6Yg2RnOw/CAaGuPSx/lEyfT+07FGWUM93JAw440ZfHP0G99rDQ2TwcQPcn7AnRPvZFLGpPC0Ijmd\nsGiRdAXt2wdjx8p9bjdkZISdK7ErEELQf906jgUJFYtOx1SbjXPj4zk3Pp7JsbGYVTzLaaFEShej\nRMqpU+x281VlpVyqqthdWwtIEdL8pT83Lo4b9u4NCHjrcYGv/tjtcu6VQ4dkUGN5OcTFwV13wc9/\nHpaF145WHWXq61MpqCnw7UuLTuPFS15kZtZMkq29pxbHiQiOq3nlslf41/5/8bctf+NI1RHGp4/n\nzol3Mj9nPtHG6FAPtzVer4xXuf76wOyxMA/K7iyCA2/HREdzQ1oaa6qq+G9VFZUeD0ZNY3KzaImL\nY2pcHNFKtHQIJVK6GCVSTozd7ebyXbs47HKh0zRi9XrfpH2DLRbOjYvjvCZh0s9sDji3Rwa+toXX\nCzk5sHdvy77sbNi5U9a6CDO2H9/O0+ue5r3d79HobUT4ZednJ2STe09uCEcXOTR6G/ns4Ge8vPll\nPjvwGTaTjZvG3MQdE+9gZEp4WcwAGZuSl9fyOjsbcnv+tT7RfahRCHY5HHzVZAleU1lJmceDQdOY\nFBvLxNhYVlVUUOf1ktmT72FnQI+tOKsITyoaGthcU8OmpmVZWRluP0EroqJ4Z8QIzo2P95Wfb4+I\nrPh6KtTXy4n+nn46UKA0E0YCRQjB53mf8+TaJ1mZt5IBcQN46oKneG/3e6wvWO9rlx6TfoJeFP7o\ndXouG3oZlw29jEOVh3h1y6u8vu11Xtz4IucOOJc7J97JVSOuwqgPkx+19PRAkZLeO671ie5Dek1j\nbGwsY2NjubdvX7xCsLeuzmcpfqWw0DdR6SGXi9GbNvG/AwYwyWZjTHS0chF1IcqS0kRvtqTUNjay\nzU+QbKqp4WCTlcSm1zMxNpYtNTVUNTb6zulVcSXtUVHREgxbVARXXCGrxO7a1dImTEzp7kY37+1+\nj6fWPsWu4l2MTx/PA9Me4OqRV2PQGcIuRTjScTe6Wbp3KS9vfpk1h9eQGp3KbeNu4/YJtzMgfkBo\nBxfm6e3hyKD168lztQRNGzUNATQIgUHTGB0dzaTYWJ/VZVR0NFG9LCBXuXu6mJ4uUprTgAvdbmx6\nPdenprLf6WRTTQ17amvxAmadjnExMb4v2ySbjSEWCzpNa+XP7dFxJSfj8GF47rmWYNibbpLVPocP\nD7sfgEpXJa9ueZXnNzxPYU0hlw65lPun3c+5A84NzyDPHsie4j28svkV3tr5FjX1NcwZMoc7J97J\nxYMv7rGZUT2Ntu5/q8aOZafDEWBx/jboXjqx6V6aZTbz67w8jvdgl7cSKV1MTxMpQgiO1tezu7aW\n3bW1PHX0KCUNgfOejPUXJCdR/70mruREbNsGTz4pZ6W12cI6GPZI1RGeX/88f9/6d+ob67lh9A38\nctovwy5GojnOOEw0XZficDtYvGsxL29+mW3HtzEwfiC3j7+dy4ddzh3/7w5lxQpjOnr/87dKN4uX\n5tg9f7LNZl4fNoyc6GiSe8h9VImULiYSREp7s4CWuN0+MeK/VDe5Z2L0ehqaykc3k2U2k9fb3TUd\nQQhYsUKKk1WrYOBAWLBAFsiKiQn16FqxrWgbT617in/u/ic2k407J97Jz8/6Oemx4Rd3IARMmCC1\nXzODB8Ojj8owHqs1cO2/fTohAOEiiIQQbCrcxMubX+a93e9R76kPCFoO2yJxitOisqGBkZs2tZp7\nqJm0qChGRUeT47eMio7GZpAho5Ey+7MSKV1MuIsUIQRTtm5lY02Nb1+cXo9Jp6O4yUJi1DRGWK2+\nD/romBhyoqPpbzJx9rZtyl1zKgRXhp0wAR54QP7KGcIr3lwIwfLc5Ty59km+yP+CgfEDuW/Kfdwy\n7hZijOElpI4fl9XbP/9cLkVFp9eP0di+gGlv/dprgRNLh0O4ULmznKEvDg2oZNvX1pejC46e4CxF\npBHsLppms/H6sGGtHiwPOJ14m9r0M5nIiY5mp8MRUIguXO/dKrunFyCEoKyhgQNOJwecTg42rQ/U\n1XHQ6QwIXAXwCMG9GRlSkERHM9hiwdCOu2ZpTk4rc6UiCLsdvv99OHAAHA4pVC65RBZjO++8sKnS\n2TyHTGFNIQadAb1Oz77SfUzMmMg/r/4nc0fMxaALj6+20wn//a8UJCtWyGxskPXEbrwR/v1vWcW9\nmWnTYPlyqK2Fujq59t8OXrd3rKKi9b7S0sCx7d8PJSWQktJ9/49gEi2JDE8eHlAk7lj1MX708Y94\ndOaj9IvrF7rBKTqNtu6/qUYjw6Ojudqvnauxkf1OZ4BwsQdZYNZVVzNt61YGWywMsVgC1vFhPLXG\n6aIsKU10hyXFfw6bOIOB29PTsTc0BIgRfyGSbjQGfAgX2u09ezbgULJnD8yeLR/1mxk3DrZ22gNB\npzHt9WmsO7bO9zrBnMDHP/iYs/ufHfJgWK9XCpFmUfLf/8oM7YwMuPBCuOAC+W9udrN0Z5zxjBnw\nzTeB+/R6Oab582XVeJuta977RPhnVvWJ7sPlQy/nmfXPUF1fzb2T7+XBGQ+Gb+l9RZcTbIXpbzJx\nXny870HWP9YwyWBgiNUaIFwSDAb+cOgQJQ0NXeouUu6eLqazRIrH66XI7eawy8WR+noOu1y+7dWV\nlTi93oD26UZjK0U8xGplkNlMTJBbQQWvdjINDfDJJ3Juk9WrQaeTv7LNhFmRKyEES/Ys4YalN+AR\nLRMNhLrwWmFhiyhZuVIKD6tVTkt04YVyGTEi9IaoYEH06qvw1VeweLEUU2aznFJp/ny5DqpJ2K3U\n1Nfw1NqneGrdU5j0Jv7n7P/hZ2f9DLMhhINShIST3ferPB4ONgmWg36W9wNOpy8UwJ8Eg4H5qan0\nN5vpbzLR32xmgMlEusl0RvOoKZHSxXRUpDRXXj3qchGt13NZUhJlHo8UIi4Xx+rraQw6x6Bp9DOZ\nKKqvx+X3/84wGvnse9/rmj9I0T4lJfKR/f335fa4cfCDH8C778KOHS3txo6FhQtDMsTgW8UO+06e\nXPskO+07sJniqK5vKWs+Nm0cb13VfeOsc8KWLbBuLaxbD7kH5YBHjoCp02DqVBg7RsaNRArHj8N/\nlsNny2Q9Pms0zJolvX1TpoAhRJnCJXWlvLz5ZZZ+u5TUmFTuPutuLh1yKbo2fkzCwxmpCCeqGxu5\ncvfuAGuLAUDT8AT99uuBviYTA8xmn4CJNxh48/hxHI2N9DOZTvhgrERKF9NRkRJselMoFAqFojdw\nohADFTgbJgSnkWUYjXzUwSDU8oYGfp2XR2lDA8lRUTyenU1iDwx0Civq6uSkah98IF03/fvLDJ1L\nL4XY2FCPrl1q3XUs3PEm7+56l1iTjbsm3smcIZei76Yqlj/5SUuQazMWC0yYCJPPgsmT5b+yu104\n3flIJYQMrv38c1jxOZQUy3iZCy6QLqyhQ7v/799p38lLG19ip30nEzMm8vNJP2d4yvDuHYQioihv\naOBBv9+dP5/C787cpgKgzbSXRt2VKJFyiqQbjQHlkbPMZs46hWi7i5OSumJYCggshBEfL904778v\nM3WuuAJefBHOP1/GnoQpjd5G3tj2Bg99+RA19TU8OO1+fjX9V92aSrxjBxz4FPCrQZWRAYd2Qm/T\n1FPOgpvPAu9vZNDt4sXw/nPwzkMwbJiMX5k/XwqW7mCybQa3DZ7Ov/b/iwdXPciPFk1mfs58/nj+\nH8lKyOqeQSgijktO83cny2wOECnpIfDhKndPEx1196jg1TBm2jRY15L1gsEAv/oV/PSn8rE/zPk8\n93N+ueKX7CrexQ3fu4HHzn+sW1NQDx+Ghx6ScyWaTOCnxcOipki40NAg6/q9+y589JHUwOPHw2WX\nwbJlUF7ePcXiPF4P/9j2D36/+veU1pVy16S7eOich0i2Jnfdmyp6Fafye6diUrqYcC/mpmgHr1em\nZixaBK+/Lm30zWRlBc72GqbsLdnL/Z/fz7IDy5jRfwbPXPgMkzInddv7l5fDH/8IL70ECQnw8MNw\n+eVw3XWhr84a7jidstbL4sVSsPh//LpL2NW6a3lu/XM8/s3jCASJlkR0mo7M2ExVYl/RbSiR0sUo\nkRJh7Nkjhck778DRo1KQeDxyu5kwf/wvrSvl4dUP88rmVxgQP4AnZj/B3BFzu63WidMpJ3D+05+g\nsVEanRYsCMtq/xFBVhYcOtTy2maTFW4TE7vn/UtqSxjzyhiKHC1lfFWJfUV30VUiJXyd8wpFMIWF\n8PTTMtYkJ0cWurjsMhkskJsLmzdLYZKdLddLl4Z6xG1S76nnqbVPMfiFwby9823+PPvPfHvXt8wb\nOa9bBEpjI/zjHzKO4qGHZOXX3Fz43/9VAuVMyMwMfF1bK+NWXnstsPxOV5ESnYIlyhKwz1+wRDx2\nu6zIN2iQXBcXh3pEim5ABc4qwpuaGmlHX7RIBgJERUlfxMMPyyIW/v7R1NSwtZw0l7I/WH6QSlcl\nDY0N3DnpTn5/7u9Jie6euuxCyJiJBx+E3bvh2mulm2fw4G55+x7P0qWBxeJeeQWeeEJmSr32mqwZ\nKB80u470mHTyKvICXvcY5s1rKRmclyf/2WH6fVd0HsqSoggv7HZpBUlPl5OqpKbCzTfLaMVXX5VV\nt95/X86xE0EByxcvuphvjn6DvdZOfWM9Y/qM4aU5L3WbQNmwAWbOlIan5GTYuBH++U8lUDqTZo2c\nmyvXOTnw1luwZo3MhJ80Ce64A8rKTt7X6bL0uqVM7zed7IRspvebztLrwtOaeFoEz0a5cSP87ndy\nGm0VttBjUSJFER5UV8taJt/7HqxdK8VIaSmkpcm0ky+/hB//WKYWRxCldaXc9e+72G7fHrC/yq9i\nbFdy4ABcc42smlpWJoM8v/hC/mAquoezz5ZTQD33nAywHTYM/v73rnEBpUan8vWtX5N7Ty5f3/p1\nzwqaTQ+yCsXHy7IC48dLF+9998kg+sbgmt+KSEaJFEXoyMuTMwxfcIF8vL/mGjl9rT96fUSkDwfj\n8Xp4ccOLDH1xKO/uepes+MAaFl1thrfb4Wc/g5EjpRXlzTdh+3aYMyf0c+j0RgwGuOce+O47ac26\n/XYpHDdtCvXIIoilSwNjznbvlnEpK1bAxRdLBXjOObKoz+23yyKOISg+puhclEhRdB8ej3zS+fWv\n5a/noEEypUSvh2eekakQZ50VeE7w01MEsCpvFWNfGcu9/7mXq0dezXd3f8f629Z3ixne4YA//EG6\ncd59Fx57TFZNvflm+W9WhJa0NCkYv/5a/n5OnizL+HSlC6jHEOxPS02VMWoXXAAvvwwFBdIKe9NN\n0lx4ySXSZXz99TIoaOpUFXQbgagU5CZUCnIn4l/5NSUFbrlFipPPPpNFOVJT5ePkZZfB7NmB5emD\np6qNoAId+RX5/HLFL/lo30dM7zedFy55gfHpXftZav5XFxZKC0l1tYw1vvtu+M1vui/9VXHqeDzy\nt/V//1cKyMceg9tuU2KyUxACdu2SQfcffRQ4cSjIyUO3bAnr6tORhqqT0sUokdJJNDbKFIa2bgqX\nXSYzcyZO7FE3h1p3LX/6+k88tfYpkq3JPHHBE8zPmd8t6cQzZrQkPIDUc5s2RaSHrNdit8uMqzff\nlF+Nv/yltUFRcYb07x9YQwkgKUlGk8+aJZfBg5Uv9AxQEwwqwhO3Wz6RfPWVTGP45hv5OO9P//4y\nAr+HIYRg8e7F/OrzX1FaV8oD0x7gwRkPEm2M7pb3r6+Xbnl/YmKUQIk00tJk3Zqf/ETGEU2ZIi0q\njz0mQ7UUnUCwSMnJgSuvlGUNfv5z+XDVt68UK+efL9fBhW8UIUGJFMWp4XTC+vVSkKxZI+fKcTrl\nr+P06TLe5L33pKm1mX7dN/9Md7G1aCv3fHYP3xz9hrkj5vLUBU916wRv+/ZJV3uwHozAEB5FE9Om\nSSvYK6/IInsffii/Tp98IpPdIsz7GV4EF7Fp/kc++qj0j65ZI+NYVq2ChQvlOcOGScEyYYIsf1Ba\nqi5CKBBChMUC/AzIR869uh6YdIK2I4EPmtp7gXvaaPP7pmP+y7cn6HM8ILZs2SIUQojjx4WYPl2I\ngQOFGDlSiHvuka+jooQAIRIShLjiCiGeflqITZuEaGhoOddul22zs+Xabg/d39FJHK85Lqa/Pl0M\neHaASHsyTfAwYtRfRomVuSu7dRxerxB/+5sQFosQw4YJ8fnnPe5frRDyOt5yi/yq+S/Tp4d6ZL2A\nkhIhliwR4o47hBgypPVFGDRIiG3bAu95CrFlyxYBCGC86ERtEBYxKZqmXQcsBG4HNgILgGuAoUKI\n0jbaT2w6vgV4FnhcCPFCUJvfA/OAWUCzo9EjhChvZwwqJqXZf7B1q4zms9tbjkVFSfPouefKNL9R\no3pUXMnJmPb6NNYda5lhOSs+i+/u/g6DrvuMkWVl0iXw0Ucyw/KZZyC6ezxLihCRmSmDopvJzpbJ\nLYpuZMAAOHKk9X6LRdZoOeusliUrq9fGtfT0mJQFwN+EEG8BaJp2B3ApcCvwRHBjIcRmYHNT28dP\n0K9HCFHS+cPtAdTVwc6dUpBs3SrjSnbvlikHen1rAdKvHyxZEpqxhhCP18PC7QvZWLAxYL+mad0q\nUFatkpmVLpe0Nl91Vbe9tSKEZGUFihTlzgsB/foFipQpU+DJJ6VvbuNG6Y979ll5LClJipVJk+R6\n4ECZYx6B2YrhQshFiqZpUcAE4LHmfUIIoWnaSmDqGXY/RNO0AsAFrAN+I4Q4epJzeg7N+akFBTLN\n99prZQnSrVth714ZLBYVJYPIJkyQj+fjx8Po0XDhhYFpI73s7iiE4MO9H/LQFw+xv2w/ydZkSuta\njHrdNSeK2y3jE556SiYivPWWiufrTbQVSqHoZtqLZ5kxo6VNaWmLaNm4Ef76V3jkkcB+8vKkFfqj\nj2QmUVRU9/4dEUrIRQqQDOgBe9B+OzDsDPpdD/wI2A+kAw8DazRNyxFC1J5Bv+GLwyEjKvfuletX\nXpF1SZrZvVvmOM6YIctfjh8vBYrJ1LqvXnx3XJm3kt+s+g2bCzdz0aCLeHfeu/S19WXuP+dS5Cgi\nPSa9W+ZE+e47GRy7Ywf8+c9w//29ysOmIKznzOw9dOQiJCfL4nGXXCJfCyGn85g8ObBw3P79spBl\nVBQMHy7d5jk5LeusLFUoJ4hwECntoSGDcE4LIcRyv5e7NU3bCBwGrgX+0d55CxYsIC4uLmDf/Pnz\nmT9//ukO5fTxL4rWLBRSUuT+ZjHSLEj27oVjx1rO7dtXunT8GThQqvyO0AvvjhsLNvKbVb/hi/wv\nmNJ3Cl/e/CXnDTzPd/zrW7vn/yEEvP463HuvvIzr13f97LkKhaIT0TR5vx0yJFCknHWWnBp7927Y\ns0euly9vmQ7EbJYiZtQoGQvz4Yfy4bNfP2mBCRNX0eLFi1m8eHHAvqqqrpmPLOSBs03unjpgnhDi\nX3773wTihBAn9L5rmpYPPBscONtO243A50KI/2njWHgFzlZUSPu+f1G02Fipsisr5WuDQZoNhw+H\nESPkMny4XGJjW1f6mj691wmPjrC3ZC8PffkQS/cuZVTKKP54/h+5YtgV3VKMLZjycul1+/BDWSvj\nuedUcKxCEbF0pIK2EDLH3F+47Nkj3Uf+kyWazXDeebK0/6BB8t4/aJC0vlgs3fpntUWPDZwVQjRo\nmrYFmYXzLwBN/jrMAk4qPDqKpmkxwCDgrc7qs0O0ZQ1JTZVWjkOH5Hw1bS1tqVIh4IEHWsTIyfya\nvdhl0xGOVB3h4dUPs3DHQvrZ+rHwyoX8cPQP0etCY2798ku48Ub50fjgA/mxUSgUEUxHLNKaJu/P\n6elyHqJmBg2ScSzNmExy+eoreOMNWZ+qmczMFtEyaJAM4H35Zfk7kpkZ0QG7IRcpTTy2J2CJAAAg\nAElEQVQDLGwSK80pyFbgTQBN094Cjgkhftv0OgpZK0UDjECmpmljAIcQIrepzZPAp0gXTybwB8AD\nBNqougKvF0pKZMDqDTdIVwzID9zgwWC1tk7vHTBAKuKzzoLrrpPbjz0mM3CaGTMGfvvbjo+jF7ps\nToTdYWfeknkcqz5Gg7eBktoS4s3xPHfRc9w+4XZMhjZic7oBtxt+/3t4/HH5oPTWW9LNo1AoejHp\n6YEiJScHPv5YbgshHz5zc1uWgwdlEc2PPw6MRTx0SLqeZsyQbqO+feXaf9tmC3zv9h6uQ0BYiBQh\nxBJN05KBR4A0YDtwkV/6cF+kwGgmA9hGS8zK/U3LV8D5fue8CyQBJcDXwBQhxKnPN+p/wdLSZLqZ\nyyVFSFtLURE0NLT3x8Idd0gR0rxkZLQdLDVzprKEdCJXLL6CjYUtMTn9bf3ZfdduYk2xJzirA/i7\nTIPdp+0cs9vh2msER49Kz57LKXj8UbjvvqaPQv2p93nax7q6fzXmkx8L13Gd7rFwHVckjfnHP5bF\nkSoqICFBvl66tPV5MTHyAXbMmJZjd90lM46acbth5crW5zZjswWKlhUrWtKu8/Lk71CIHnhDHpMS\nLpwwJiU4tuPUOg78YJhMLem8PeGLdKbHuuG9vcIrF6+X4CiTKE3f7nknfD+FQqHoLXSgimCPjUmJ\nCIqKAl/HxckCPf6BlcFBls2va2ul+nU4pOKdN0+uT3ZeZx3r6v7D9JgQgr1l+1iVv4q9JfuIN8eh\n0+kpc0ozqAAGJWbzwLQHun3M9W6Ne+4BVz1NY9FITZG1UMLl/xdW7x0px8J1XOF0LFzH1d3Hurr/\nkhIZ4Hb8OPTpA++8IzNDO9rn7NlyXrZmQlgnS4mUjtCWb/DxExW6DeL55zt/TIo2cXlcvLPzHZ5d\n/yx72MP4aeO5b8oirhl1DZWuyoBaJ49ctxSiu9fP+u23sqbe3gY5mVQz04cCN3XrUBQKRU8lLi5Q\nZJwqH38cNqEGSqR0BJUlE/YU1xbz8qaX+cumv1BaV8oVw67gr5f+lbP7n01zKnFqdGq31Tppi7ff\nluFIAwfC6tXwm9+oj5RCoQhDwijpQomUjhBGF0wRyJ7iPTy7/lkW7VyEXqfnlrG3cO/kexmSNCTU\nQ/NRVwd33y2zBm++Gf7yF1n7RH2kFAqF4sQokaKIKOwOO3OXzCWvIg+Xx0Wlq5KM2AwePu9hbp9w\nO4mWxFAPMYB9++Caa2TM2T/+AT/6UahHpFAoFJGDEimKiGF38W4uXnQxBTUFvn1DE4ey665dGPXG\nEI6sbRYtku6d/v1l8chRo0I9IoUicmmuc+Q/f1ZqN8eUKbofJVIUYc2x6mMs3rWYRbsWsdO+Ex2B\nM+x5hCfsBIrTKedvfO01GWD/178GJnQpFIpTZ96SeXxzVJaCyKvI48r3rmTtj9eGeFSKrkbNqaoI\nOypdlby29TVmLpxJ/2f787vVv2N48nA++cEnTO47OaBtekzoUuPaYt8+OfHpokVyksCFC5VAUSg6\ngyJHYCmIDcc2cOsnt/LZgc9wN7pDNCpFV6MsKYqwwOVxsezAMhbtXMS/D/wbj9fD+Vnn88b332Du\niLnYTLJs85S+UwLSiJdeFz5pMe++KycH7NdPTjY9enSoR6RQ9BzSY9LJq2gpBZFhy+Cbo9/wj+3/\nIMGcwJXDr+TaUdcyK2sWUfoTzGmmiCiUSFGEhOYA2PyKfLzCS11DHTXuGiakT+DPs/7MD3J+QHps\naytJqNOI28LphHvvhb//HX74Q3jlFWU9USg6m6XXLW31gJJiTWFX8S6W7FnCkj1LfILlquFXcc2o\na5Rg6QGosvhNnLAsvqLTsDvsrMhdwS9X/JKSuhLf/r62vnx+4+cMTx4ewtGdGnY7zJkj5/TyeODp\np+EXv2hdwFGhUHQ9Qgh22nfy/rfvs2TPEg6UHyDRkshF2Rexo3gHzgYnGbEZKuC2i1Bl8RURibvR\nzdqja1l+cDnLc5ez7fg2gFbBrka9MaIEitcL06YFFiL+8ENYsCB0Y1IoejOapjGmzxjG9BnDozMf\nZYd9B+/veZ9n1j+Dy+MCIL8ynymvTeHT+Z8yMmWkr9CjInxRIkXR6eRV5LH84HL+k/sfvsj/Aofb\nQYo1hYsGX8R9U+/jguwLAiL1IfwCYE/EoUNw662BAgVaT/GkUChCg6ZpjO0zlrF9xvLe7vfIq2z5\nsuZX5pPzcg59Yvpwftb5zMqaxaysWQyIHxDCESvaQ4kUxRlhd9i58p9Xkl+Rj07TYTKYOFR5CIPO\nwLR+0/jtjN9y0eCLGNtnLDqtJZmsLf9yuCMEvPoq3H8/JCbKuid79rQcD+EcXAqFoh3SY9MDRMrU\nvlN5+LyHWZW3ii8OfcHiXYsRCAYlDJKCJXsWMwfOJCU6RdVmCQNUTEoTKialYwghOFB+gA3HNrD+\n2HoW7lhIbUOt73if6D68fNnLnJ91vi8jpydw5Ajcdht8/jn85CdyxmKXq/WUTqnq/qVQhBXFtcWt\nHoj8hUaFs4LVh1azKn8Vq/JXsa90HwBj0sZQUldCYU2hr+30ftPDLnA/XFAxKYqQUOGsYGPBRtYf\nW8+Ggg1sKNhAubMcgOHJw1v5dK1GK1cOvzIUQ+0ShJDl7BcsgNhY+M9/4KKL5DGbTc2/o1CEOyfL\nCEywJHDViKu4asRVABRUF/BF/hesyl/Fop2LAtpuP76dJ755gokZExmfPp54c3yXjl2hRIqClnLT\nhTWF2Ew2rh99Pd+WfMv6Y+vZX7YfgERLIpMzJ3Pv5HuZ0ncKkzImkWBJYMYbMyI2tuRkFBRIq8ln\nn8k5d559FuLVPUmh6NFk2jK5ccyN3DjmRg6WHwy4v+l1eh5d8ygOtwOAIYlDmJgx0beM6zOOWFOs\nchN1Isrd00RvcvcIISisKWRPyR52F+/m8W8ep7i2OKDNhPQJTOk7hcmZk5nSdwqDEwe3GQl/MlNq\nJCIEvP22LG1vtco4lMsuC/WoFApFd9PW/S3Zmsz+0v1sLtwsl6LNbCvahtPjRENjRMoISmpLAkos\n9AY3UVe5e5RIaaIniJS21DvAnmIpRppFyZ6SPVS6KgGwGCw0isaAstID4weSf29+SP6GUFNUBD/9\nKXz6KdxwAzz/vAySVSgUivbweD3sLdnrEy6vbXst4J6q03RM7zedEckjGJEyguHJwxmRPIJ+cf0C\nEgoiGRWTomiXhsYGjlQd4fvvfZ89JTLdJK8ij8xnMvF4PUBLHZKc1BzmDJlDTmoOo1JGMTB+IOe+\neW6ASTMzNjMkf0eosNth3jw4cADKyiAhAT76CK7sOaE1CoWiCzHoDIxOG83otNHcMu4Wdth3BNxT\n+8b2JdOWyYaCDby18y1f3RZrlNUnWJoFTLI1mQdXPoi91t5jrNNnghIpYUZ7vswKZwV5FXnkVeSR\nW5EbsH2k6ghe4W3VV6wxllcvf5Wc1BwGJw7GoGv7ckdiOnBncsklsG1by+vsbCVQFArF6dPWPbVZ\naDR6GzlSdYS9pXvZW7JXrkv3suzAMipcFQH95FXkMfrl0dx91t0MjB/IwPiBZMVnkR6b3mMsMCdD\nuXuaCKW7RwhBaV0pBTUFXP/h9ewt3es7Fh0VjVFvDPjwxpniGJQ4iOyEbLLjs8lOyGZQ4iAeXPkg\nW4q2+Nr1Bj/omVBYCA8/LOfc8Sc7G3JzQzIkhULRSxFCUFxbzLi/jQuY8dmoNxJniguIcTHqjfSP\n6+8TLc0Cxmay8chXj1BaV9rtUwAod08E4W8NSYtO48VLXqS+sZ5j1ccoqC6goKZAbtcUUFBdQGFN\nIfWN9W32pdN0PDDtAZ8QyU7IJsGc0GYQ67IfLuvVFpGOUlUFTzwhs3WsVsjKgny/EBxVlE2hUHQ3\nmqaRFpNGdkJ2gEiZlDGJr2/9GofbweHKwxyqPORb8ivz2Vq0laV7l1LmLAvoL78yn0EvDOL8rPPJ\niMkgIzaD9Nh0MmIzfEuyNdlnkQnXjCRlSWki2JJyogvmbHBSUldCcW1xm8un333qC0wNxhplJTM2\nk0xbJn1tfeW23+ufL/s5mwo3+dora0jnUV8PL78M//d/UFcna5/86ldyvyrK1nk031MEImC7+Zj/\ntmp36u3CcUw9uV13j8nhdvDe7veocdcQY4zhulHXEW2MPml/1fXVvLH9Deoa6ugoBp2BPjF9yIjN\n4GD5QV8NLDj13x6V3dPFBIuU4PofZ0KULoqzB5xNZmwmcaY43/62Pmx17jr+k/sf6hrqsEZZuTD7\nQsxR5hN/yMP0yx1O7crK4FiBwO0WJCdDerogyhgZYw/3dv6vFYpIQUPzWaQ1mtaaFrDdGe26872K\naooCrPJmvZm+cX071N/B8oM0eBt852YnZJN7T8f93srd080U1XRstrgYYwwp1hRMBhMgL/ahykM4\nPU5fmyh9FMW1xZTUlnToA9jX1te3va9sX7vtwunL02wy1LTO/zKeUh9B+44VwKaNGmWlGgMGwuQZ\nGgkJXXPT6Krr1B3twnFMob7h94b/ZySP/XT76MmcSd2qcC3MqURKOwRPSnUqpq+eWOAsUmhOJz50\nCBwOGX8ybZqMQZk+PdSjUygUiq7jZFMAnIhwzfIMG5GiadrPgPuBPsAO4G4hxKZ22o4EHgEmAAOA\nXwghXjiTPoM5kwt2Jh8UxekjBFx4Iezc2bJv+HA5v04veIhSKBSK0yZcf7fCQqRomnYd8DRwO7AR\nWAAs1zRtqBCitI1TrEAusAR4tpP6DCBcL5iiNXV1sHixDIr1FygAbrcSKAqFQhGphEs1mAXA34QQ\nbwkh9gF3AHXArW01FkJsFkL8WgixBHC31eZU+1REHvv3ywydzEw5EWBaGowYEdhGpRMrFApF5BJy\nkaJpWhTSbbOqeZ+QKQMrganh0qciPGhogA8/hNmzpStn0SI5105uLvz737B6tYw9yc6W66Xh4VZV\nKBQKxWkQDu6eZEAP2IP224FhYdSnIoQUFsrKsK++KrenTZMC5eqrwWRqaZeaKmNQFApFz8Juh6uu\nkmtVz6j3EA4ipT006PTiC13Rp6ILsNtlgbX8fBlXUlEBFoucmfjOO2HMmFCPUKFQdCfz5sG6dXI7\nL0/OubVly4nPUUQ+4SBSSoFGIC1ofyqtLSFd3ueCBQuIi4sL2Dd//nzmz59/mkNRnApOJ6xcKWNM\n7H5XKjsbtm6FoEujUCh6CUVBpau2boXRo+XDzLx5clsFyXcPixcvZvHixQH7qqqquuS9wqLirKZp\n64ENQoh7m15rwBHgBSH+f3v3HiVFeeZx/PuAjKAIiIgoSMTBiEhAI2qI9wveo2sSL4nG+7oe1yOJ\nUReTqJi4OUtUTNyjJ8azqLhK1nUJSFaOJosRA3g0ghARNDqgMl6AcPECijDv/vFUn6lpema6qy9V\nw/w+59SZ6eq3qt5nqqf76fd9661wRzvbLgfuzr8EudR9Woo3GOzsVq+G3/8eZsyAZ57xRKVbNx9/\nkqOb/ol0bkceCXNjk4APGwajR8PMmT4f0tChzQnLoYcqYam1as04m/rA2cgk4Eozu8jMhgG/xi8z\nfgjAzKaY2c9zhc2sm5mNMrODgDpgYPS4vth9SrreeAPuuAOOOgoGDIDLL/dkZcIEWLYMDjusZXld\npSPSuU2b1nJQ/HPPwSOPwKpVMGsWHHssTJ4Mhx8OgwfDuHEwfbqXra/3JGfVqrSjkFJloiUFwMyu\nBm7Eu2hewSde+0v03GxgRQjhsujxl4DlbDu+5LkQwvHF7LPA8dWSUiW58SUrVvi3mx494M03/efY\nsXDWWXD66X4Jcc6qVbrpn4iUZssWHzg/bZovjY0tnx89Gl4qajpPKZVuMFhlSlIq69NP/c1g7ly4\n6y4f+JrTv79fpTN2LOy0U3p1FJHtV1MT7L23Xw0YN3gwHH20t+IefTTsv7+6hipBNxiUTGts9IRk\n3jz/+cor/q2md2+/OieuZ09vPRERqZYuXWDIkJZJyrBhcNppMGeOz1K9dSvsvrsnLLmkZdQoWLPG\nx7aoJTd9SlKkaLmb9733HvTqBeef79PQz50L77zjZerrfQ6TK67wvuDhw/0fPz7gTeNLRKQWpk1r\nvdv444/9kuY5c+D552H8ePj8c9hlF+jaFdav93INDb4Pzb+UDiUp0qY1a+DVV325/faWlwUvXuyD\n1M45xxOTr3/dB8HmK/RGISJSbW1N7rjLLn5D0pNO8seffQZ/+YsnLbff3rLsCy/AKafAyJG+jBrl\n3UR1ddWtvyhJ6ZRyLSLxpKFHD3jtteaEJLd88IFvU1fndxmO22ef5smV2qJZYEUk67p39yuAjjwS\nnnpq29bfujp4/HG/KhF8moQDDmhOXEaOhL328skm1U1UOUpSOpGtW33syOmnewIC3pQ5eLA3c4IP\nIBs61CdGuvJKGDHCl6FD4bjjWv7j7rVX7WMQEam2trqJ1q/398/Fi2HRIv/5u9/5xQJxDQ0+X8vE\nid4NXl8PffvWPpaOTklKB1aoRaSuzv85Ghp8Svnc7w0N8PbbLSdIy+ne3a+2GTHCvxn06FH4eOq2\nEZHOoK3W3z59mltccpqa/P32iCNadomvXAnxycr79PEvfLmkJbf06gXXXOMt12qBaUlJSgeycaMP\nWs0tN93kc4+AJyF77eWtJTm77OITH+27L5x5pv8cMgRuvtmnlM4ZMQIuuqj946vbRkRkW126eLIx\ndGjLJGXMGJ9N+623tl3mzvWW7fxu9IYGf0++4goYNAgGDmxe+vf3Qb1xhb6sbk8JjpKUlOXu7Lly\npWfZ48Z5s2E8Gckt+bdGyL+2v3dvuO++5sSkb9/C1/+PHq0WERGRSivU2tynDxxyiC/5PvvMW2CO\nO65lcvPRR/Doo/6+v2VL8/quXX2/Awc2JzAzZngrOXiCc8YZngB161bdWGtFSUoFxDPZAQPggQc8\nOVi92q+Oyf8Z/33lSm8qBHj3Xc+ed9zRW0Vyy4gRLR/nltNOazlG5IAD4Lzz2q+vWkRERCqv1PfW\n7t39fTu/BWb0aN9PU5N/VjQ2+rJyZfPvjY2wdGnz9A85L73k3f59+nh9dt+95RJf16WLt8ivWeOf\nKVlshdGMs5FyZpzNv/FVefXwWRLNmltB8n/mft+yxV/YW7bADjt4hp3LntvbNslPbdtxjq9ts79t\n2sfXttnZdu1auO46Txb69YNf/hJ22624fVxwQcvu+yFDfA6r1atbLmvWtJz5uzU9e8Kuu7Zc+vbd\n9ifAL37hA4n33BMmTFjA2LGaFr9qyklS6uu9mS2nrs5vdtWvny+77uoZa/xPnft98mRvQckZNAgu\nuaT5+fyfhdYV+1Pb1mbbtI+vbbO/baWOL5IVo0YtYNEiTYufSXvu2TJJOfRQePrp4ra9+mrdSE9E\nytNRkrPOtm3ax2+vzLhxsGRJc9kDD4RJk1rfdvNmb41Ztw5uu615Vl7wlppqUJJSAeVcmqvxISJS\nrkJdDCLtmT07+ZfkJ55oOcyhX79t7zpdCUpSKkCJhoiIdDTlfHblfzmfMMHvbF9pSlJERESkJPkJ\nzoKKjUJpqUt1disiIiJSHiUpIiIikklKUkRERCSTlKSIiIhIJilJERERkUxSkiIiIiKZpCRFRERE\nMklJioiIiGSSkhQRERHJJCUpIiIikklKUkRERCSTlKSIiIhIJilJERERkUzKTJJiZv9sZsvNbJOZ\nvWBmh7ZT/hwzWxqVX2Rmp+Y9/6CZNeUtT1U3imyZOnVq2lWoqO0pnu0pFlA8WbY9xQKKp7PJRJJi\nZucBdwG3AgcDi4CnzaxfK+XHAI8BDwAHAdOB6WY2PK/oLGAPYEC0fKcqAWTU9vbi357i2Z5iAcWT\nZdtTLKB4OptMJCnAD4D7QwhTQgjLgKuAjcBlrZQfB8wKIUwKIbweQrgVWABck1fu8xDC6hDCqmjZ\nULUIREREpKJST1LMrBtwCPB/uXUhhAD8ERjTymZjoufjni5Q/lgz+9DMlpnZfWbWt0LVFhERkSpL\nPUkB+gFdgQ/z1n+Id9EUMqCI8rOAi4DjgRuBY4CnzMzKrbCIiIhU3w5pV6ANBoSk5UMIj8eeW2Jm\nfwXeAo4Fni2wfXeApUuXllzRrNqwYQMLFixIuxoVsz3Fsz3FAoony7anWEDxZFXss7N7Jfdr3rOS\nnqi7ZyPwrRDCk7H1DwG9QwhnF9jmbeCuEMI9sXUTgLNCCAe3caxVwI9DCA8UeO67wKNlhCIiItLZ\nXRBCeKxSO0u9JSWE8IWZvQycADwJEHXJnADc08pm8ws8PzZaX5CZDQJ2A95vpcjTwAXACuCz4iMQ\nERHp9LoD++CfpRWTeksKgJmdCzwM/BPwIn61z7eBYSGE1WY2BVgZQvhRVH4M8BwwHvhf/NLi8cBX\nQwivmdnO+OXM/wN8AAwFJgI7AyNDCF/UMj4REREpXeotKeDjR6I5UX6Kz2vyCnByCGF1VGQQsCVW\nfr6ZfQf412j5G97V81pUZCswEh842wd4D8/ublGCIiIi0jFkoiVFREREJF8WLkEWERER2YaSFBER\nEcmkTpWkVPomhmkrJR4zG25mT0Tlm8zs2lrWtRglxnOFmc0xs7XR8of2zmctlRjL2Wb2kpmtM7NP\nzGyhmV1Yy/q2p9T/ndh250evt2nVrmMpSjw/F0cxbI3drHRjLevblgTva73N7F4zey/aZpmZnVKr\n+ranxHPzbIEbyTaZ2cxa1rk1Cc7N96PzsdHM3jGzSWa2Y63q254Sz80OZnaLmb0ZlV9oZieXfNAQ\nQqdYgPPwS4svAoYB9wNrgX6tlB8DfAFcB+wP3AZ8DgxPO5aE8YzGr3A6F2gErk07hjLjeQS/x9NI\n4MvAZGAdsGcHjOVo4KzodTYEuDZ67Y1NO5Yk8cS2+xLwLvAnYFracZRxfi6OXlu7A/2jZfe040gY\nSzfgJWAm8DVgMHAU8JW0Y0kYT5/YOekPDI/+d77XAWP5LrAp2m4wcGL0Xn1n2rEkjGdi9P9/Mn5p\ncu6efKNKOm7agdfwD/wC8KvYYwNWAje2Uv63wJN56+YD96UdS5J48rZdTvaSlMTxROW7ABuACzt6\nLNE2LwO3pR1L0nii8/E8cCnwINlKUkp9L7gYWJt2vSsUy1X41ZBd0657JeIpsP33gfVAj44WC/Dv\nwB/y1t0JzEk7loTxNAJX5a17AphSynE7RXePVfcmhjWXMJ7MqlA8O+PfEtdWvIIlqEQsZnYC3jr0\nXDXqWIoy4rkVWBVCeLC6NSxNGfH0NLMVURP8dDMbXuWqtithLN8g+rJlZh+Y2V/N7CYzS/2zoELv\nA5cBU0MImypfw+IljGUecEiuC8XM9gVOw+cCS1XCeHbEex/iNgFHlnLsTMyTUgNt3cRw/1a2KeYm\nhmlJEk+WVSKeiXjmnp9Y1lqiWMysF17/HfE5ga4OIcyuViVLUHI8ZnYE3oIyqrpVSyTJ+Xkd//Bb\nDPQGbgDmmdmBIYTGalW0CEli2Re/6ep/AqcC+wH3Rfu5vTrVLFpZ7wNmdhhwIP7aS1vJsYQQpprP\nF/ZnM7No+1+HECZWtabFSXJungauM7Pn8fvmnQh8kxLHwnaWJKU1Zd3EMIOyXr9SFRWPmY3Hx9oc\nE0LYXPVaJdNeLB/jH+o98Vs+3G1mDSGEObWoXAIF4zGznvh4oX8MIayrea2Sa/X8hBBewJu6vaDZ\nfGApcCXeYpQ1bb3WuuAfLFdG34QXmtlA4HrST1JaU+z72uXAqyGEl6tcn3K0GouZHQv8CO+SexGf\nKf0eM3s/hNARz8044DfAMqAJT1QmU2IS2VmSlDX4LLR75K3vz7aZYc4HJZavpSTxZFnieMzseuBG\n4IQQwpLqVK8kiWKJPjAaooeLo+6Em4C0k5RS46nHB8zOjL4NQvTNycw2A/uHEJZXqa7FKPt/J4Sw\nxcwW4h8iaUoSy/vA5uj1lrMUGGBmO4QQtrSyXS2U8z7QAx/Y+ZPqVK1kSWL5KT5eI9dFuiRK+u8n\n/QSy5HhCCGuAb5pZHbBbCOF9M/s3fExk0VLvh6yF4FPh525iCLS4ieG8VjabHy8fafMmhrWSMJ7M\nShqPmd0A/Bi/hcLCatezGBU8N13wrp9UJYhnKfAV4CC8ZWgUfuPQ2dHv71a5ym2qxPmJxm+MoPWb\nldZEwljmsm1ytT/wfsoJSrnn5jygjozcyT5hLDvhLQ5xTdGmVqB8zZRzbkIIm6MEpRvwLWB6qQfv\nFAveHbCJlpdP/Z3oUkJgCvDzWPkxwGaaL0GegF9+lZVLkEuNpxv+IXEQPvZhYvS4Pu1YEsZzY3Q+\nzsaz+9yycweMZTzeXzskKv9DfMDZpWnHkiSeAttn7eqeUs/PzfgXlCHAwcBU4FP8BqgdLZZB+FVw\nv8LHo5yOtxqPTzuWcl5r+JVkj6Vd/zLPza34lUnn4ZfsjsWvxMpEXAniOSx6fx6CX+b+R+BNoFdJ\nx0078Br/ka8GVkR/6PnA6Nhzs4HJeeW/hfenbcIHzZ2cdgxJ48Gb4JvwJrv4MjvtOBLGs7xALFvx\nm0h2tFh+hg/O/BRvVv0z8O20Y0gaT4FtM5WkJDg/k6LX2yb8ZqUz8buppx5HknMDHI5/+92Ifwj+\nC9F93LKwJIhnv+h///i0617m66wLnhC/Eb0XrADuocQP9QzFczSwJHqdrYreBwaUekzdYFBEREQy\nqVOMSREREZGOR0mKiIiIZJKSFBEREckkJSkiIiKSSUpSREREJJOUpIiIiEgmKUkRERGRTFKSIiIi\nIpmkJEVEREQySUmKiHQIZrbczK6t0r6fNbNJ1di3iCSnJEVE2mVmD5rZtJSrMRr4Te6BmTWZ2Zkp\n1kdEqmyHtCsgIlKMEMLf066DiNSWWlJEpCxmtreZzTCzj81sg5n9l5n1jz1/q5ktNLMLoy6b9WY2\n1cx2jpXpaWaPmtknZtZoZt/P74KJd/eY2XIgANOjFpWGaP1D+S0+Zna3mT0be59Q/LAAAAMVSURB\nVLyTmU2J6ttoZtcViKnOzO40s5VRneab2TGV/LuJSPuUpIhIuWYAfYCjgBOBeuC3eWXqgbOA04DT\ngWOA8bHn7wbGAGcAY6N9fbWNYx4KGHAxMCB6DJ64FBJff2e0/28AJwHHAofklb8XOBw4F/gK8N/A\nLDOrb6NOIlJh6u4RkcTMbCwwAtgnhPBetO57wBIzOySE8HKuKHBxCGFjVOYR4ATgZjPrCVwEnB9C\n+FP0/KXAe60dN4SwxswANoQQVpVQ352By4Dvxo51MbAyVmYwcAmwdwjhg2j1JDM7FbgU+EmxxxOR\n8ihJEZFyDAPezSUoACGEpWa2HjgAyCUpK3IJSuR9INcltC/+XvRSbB8fmdnrVahvPdANeDF2rHV5\nxxoBdAXesCgTitQBa6pQJxFphZIUESmHUbiLJX/9F3nPB5q7my22Ln8fpWoqsF23AvtsrVsIoCew\nBe9uasp77pMEdRKRhDQmRUTK8Row2MwG5laY2XCgd/RcMd7Ck4LDYvvoBezXznZf4C0ecauBPfPW\nHRT7/c3oWF+LHWtX4MuxMguj/e4RQmjIW4ruWhKR8qklRUSK1cfMRuWtWwYsBh41sx/grRb3As+G\nEBYWs9MQwidm9jBwp5mtwxONCcBW2m7xWAGcYGbzgM9DCOuB2cD10biY+cCFePfNguhYn5rZfwB3\nmNna6Fi3R8fK1edvZvYYMMXMrseTlv7A8cCiEMKsYuISkfKpJUVEinUM/mEfX24B/gFYDzwHPIO3\nVpxf4r5/AMwDZkb7+DOeAH0WK5OfsPwQvxLoHZqTkGeAnwET8XEnPYGH87a7AXgeeDI61vM0j53J\nuQSYgl8JtAz4HT6Z3DslxiUiZbAQ2vqiIiJSe2a2E9AIXBdCeDDt+ohIOtTdIyKpM7OD8CuFXsTn\nXLkFbzmZkWa9RCRdSlJEJCuuxwewbsa7X44MIaxNt0oikiZ194iIiEgmaeCsiIiIZJKSFBEREckk\nJSkiIiKSSUpSREREJJOUpIiIiEgmKUkRERGRTFKSIiIiIpmkJEVEREQy6f8BZyUl4LSBdUMAAAAA\nSUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plotTrajectoriesFile('PeriodicParticle.nc')" ] }, { "cell_type": "markdown", "metadata": { "deletable": true, "editable": true }, "source": [ "We can see that the particles start at 0.7E, move eastward, and once they hit the boundary at 0.895428E, they jump to the other side of the domain (the horizontal lines). So we have periodic boundary conditions!" ] }, { "cell_type": "markdown", "metadata": { "deletable": true, "editable": true }, "source": [ "*As a note, one may ask why we need the halo. Why can't we use simply the `PeriodicBC` kernel? This is because, if the particle is close to the edge of the fieldset (but still in it), `AdvectionRK4` will need to interpolate velocities that may lay outside the fieldset domain. With the halo, we make sure `AdvectionRK4` can access these values.*" ] } ], "metadata": { "kernelspec": { "display_name": "Python 2", "language": "python", "name": "python2" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 2 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython2", "version": "2.7.13" } }, "nbformat": 4, "nbformat_minor": 1 }