{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Coarsen an existing MOM6 grid\n", "\n", "In this notebook, we create a coarse versions of an existing MOM6 grid, tx2_3v2.\n", "Setting CX=5 below results in a grid with ~3 degree resolution, while setting it to\n", "15 results in a grid with ~10 degree resolution. " ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "# CX: Coarsening factor\n", "\n", "# CX, grid_name = 15, 'tx10deg' # 10-degree resolution\n", "CX, grid_name = 5, 'tx3deg' # 3-degree resolution" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## 1. Import Modules" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "%%capture\n", "import xarray as xr\n", "from datetime import datetime\n", "from mom6_bathy.grid import Grid\n", "from mom6_bathy.topo import Topo" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "datestamp = datetime.now().strftime('%Y%m%d')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## 2. Create a MOM6 grid object from an existing supergrid\n", "\n", "tx2_3v2 is the current MOM6-CESM workhorse grid (0.66deg, tripolar)." ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "grid_orig = Grid.from_supergrid('/glade/p/cesmdata/inputdata/ocn/mom/tx2_3v2/ocean_hgrid_221123.nc')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## 3. Read the bathymetry" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "topo_orig = Topo.from_topo_file(grid_orig, '/glade/p/cesmdata/inputdata/ocn/mom/tx2_3v2/ocean_topo_tx2_3v2_240501.nc')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## 3. Coarsen the grid and topography:" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "topo = topo_orig[::CX, ::CX]\n", "grid = topo._grid" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "(96, 108)" ] }, "execution_count": 7, "metadata": {}, "output_type": "execute_result" } ], "source": [ "topo.depth.shape" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjsAAAGwCAYAAABPSaTdAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8ekN5oAAAACXBIWXMAAA9hAAAPYQGoP6dpAACf40lEQVR4nO2deXxU1fn/P3f2yTKTBbKRBAKCgEJRUIzaFpVKLfp1obb2Sy0q6lcLrUB/1doqtVRFbV1blLoU7bdSq/1WW3cRFbUCIoICKotsYUkCZJkks8+9vz+QOed5QkIgYSYkz/v1yut1z5xzzz333HMnZ87zOc9jWJZlQRAEQRAEoYdiS3cDBEEQBEEQjiYy2REEQRAEoUcjkx1BEARBEHo0MtkRBEEQBKFHI5MdQRAEQRB6NDLZEQRBEAShRyOTHUEQBEEQejSOdDfgaGOaJnbt2oXs7GwYhpHu5giCIAjdGMuy0NTUhJKSEthsR289IBwOIxqNdroel8sFj8fTBS3q2fT4yc6uXbtQVlaW7mYIgiAIxxBVVVUoLS09KnWHw2FU9M9CdW2i03UVFRVhy5YtMuE5BD1+spOdnQ0AOBPfgQPONLemm2GoXy2G3d7h02zZmSRdNfV4ko76lVPuRH6MnmynDrsHPK2Ot51Pn4+nmv6qyt6hncscfzf3Y7/AtEW83PVxkpW1aidJmw0NJG3LVPeXaAiQPLsvm6StaESVbW6BcPSwuVz0Awcbs9qYMENhkmXPziLpRFOzqqagL60nzv4BedR1oxUFJCuaQ9vUUqzaZGM/2utG0XqNLPpu2GrcyWPTTce3wfzcZ29S473wb2tppoN9rWsr2lYwRLOctGyiRcu3TPRG4ojhfbyS/N9xNIhGo6iuTWDbygHwZR/56lGgyUT/0VsRjUZlsnMIevxk54DpygEnHIZMdgj6ZMc4jMmOjX7B2930JbN51Dez5WX1ssmO/r1s89LnY3fTLwGHs+3JDi+rT3YcTjrZcbD2mwZN6/dnsDFjZ+da2n8hw+j8krTQNjb+/hr860ub7Bh0YmFnz1h/rnw8wMYmOzZtEuKgY910snfBpU12WOtsXjbZYe+GzaOuAw+b7LB5h92lanewe2vVLzZtssP6xWBlDUN/V3rnZOfAMEqF7CEr20BW9pFfx4RIMzpKj5/sCIIgCEJ3JGGZSHQiOmWil66+HQky2RH2w16aVmYtbRXIDDSTLJMvCuk/aW3sTbboL5F4pq2tLETy6LlWlToO59HfzVE/a4L2A7XueDrMTUc5SWe/0UTzNXOUFedmOHpdqwsEhkBrE43ZRfV2S7SxdChTieFQKzCt+oSlbe62l/GtRNvaiHhNbYfb4GK/9g1m1vLUqvtx7KXvyb4xfVib2NpPO9YMi+UV/2Nj8tg0WR+yfjF0ka2Ntt9i59o0s5ZlsveP96H8o+00JiyYOPLZTmfO7W3I1nNBEARBEHo0srIjCIIgCGnAhNkpZVTnzu5dyGRHEARBENJAwrKQsI7cFNWZc3sbMtnpzWg2d8tkWhQz3mZZexbdxss3BLj3qQ+C+Swzyq5jV/m5n9E87x76qyXuVWXtTEpjuuhLb2obW+xh2oamcioyymb6JLNZ25acn0/bULsHRwW+Xbgna3Z0jPYt6bp2rJV+ip9ra3tnihkMkrRd21acaGrixQn6deM7qNsCe109STd9Z2TyuPlU6qLB3kTHaILp2eJ+9c4Z7D2xhdm96uOF/cMz3G5aNq7qtWL0veZ9qr/bVohuU7exbeo9Wlcm9DhksiMIgiAIaUAEyqlDJjuCIAiCkAZMWEjIZCclyG4sQRAEQRB6NLKyI+yH+9lxUG+1Vlzl7710BMkzmWNbPVyElaA6iv4v0bINA5UmI8akQO4Gmg71UXXFM2ge/4FjaSM73Id5W2Yanpr/PoGkC//+uUq0owPhtA5nwDQOTDfS0bwehz7WuO6GjcPD0YWYTGOiY8ugA+ZQOp0Ok6Dt9b35hTpmWhprYD+S3nCFj6QdIdUXTta8SD7TpNU3JI9b+cSKM12O5kvHcNGXlZ+baKa+gYSji5ixUodMdgRBEAQhDchurNQhZixBEARBEHo0srIjCIIgCGnAROfCrYpLwY4jkx3hoLTnf4Nj8JA5WnRy7zaqEQiUs7J6dPIWmpdwUb1MWAst1Mqvjqvt5VwevNhkbkiChbRA9KSByWPXR5varLfVdZguxOB+SXqTLqejHCq+kpbPdTdcm6LrexzFRbTo7uo2L2HPyWHXZFobrV6uCzIjYXouT+usaiTJrG+eTtIxLb5b6Vv0ZbB9vL7N69q83ravCZAXoJW+x0nfT7LUz/VTzEeP0HkSndyN1Zlzexsy2REEQRCENJCw0Mmo513Xlp6OaHYEQRAEQejRyMqO0CHMcCR5nFlNl7MjuWw7qxZ6omg53TocLKJlQ3naEju7pncftY8FYmq42mK0dIytzlt2bfs7r5iVTXhoevdpys7lPY5uSy947nOSNrSt04nGAG0DD28gdI52zFatijKzlaOkmBaIaOc66ICI19QeWfsOk+I/LCfpqpvHJo/tTdQclmjPPJag7wkPF2FF1TjkCwGtzFr6eSb35yAKka5GNDupQyY7giAIgpAGTBhItPqZd3jnCx1DzFiCIAiCIPRoZGVHEARBENKAae3/68z5QseQyY7QISxNF+B+42OSV7yILhDqLukNG83zsH3g1pD+yePaU6n7fE7hSqXZyPxoO8kLn1BK0s4mpcnY+N9sy7LBthbzt0Broy1Gyxo5fpKOb97aXpOFLsTgW6wPI5SElZtN0qZX6Vpsm3d0ql1HDAuVMWDB1uTx3vH9SV5G6Skk7XpthUow3Y2u0QEAI1P1m6Vp7wC00vuQ8ChxUYQcbRKdNGMd7rkDBgzAtm3bWn3+4x//GPPmzUM4HMbPfvYzPPPMM4hEIpgwYQIefvhhFBYWJstu374d119/Pd5++21kZWVhypQpmDt3LhxaeJx33nkHs2bNwrp161BWVoZbbrkFV1xxxRHfZ1cgZixBEARB6AWsWLECu3fvTv4tWrQIAHDppZcCAGbOnIkXX3wRzz33HJYsWYJdu3bhkksuSZ6fSCQwceJERKNRfPDBB3jqqafw5JNPYvbs2ckyW7ZswcSJE3HWWWdh9erVmDFjBq6++mq8/vrrqb1ZhqzsCIIgCEIa6KqVnUCA7gR1u91ws115ANC3b1+SvuuuuzBo0CB885vfRGNjI5544gksXLgQZ599NgBgwYIFGDZsGJYtW4bTTjsNb7zxBj777DO8+eabKCwsxKhRo/Db3/4WN910E2677Ta4XC7Mnz8fFRUVuPfeewEAw4YNw/vvv4/7778fEyZMOOJ77SyysiMIgiAIacC0jE7/AUBZWRn8fn/yb+7cuYe8djQaxV//+ldcddVVMAwDK1euRCwWw/jx45Nlhg4divLycixduhQAsHTpUowYMYKYtSZMmIBAIIB169Yly+h1HChzoI50ISs7wsFhegLdx4Zl8rxYm2VNpqCzsRAK+Hxz8rBoN9PDlBfQap3quvu+VUHy8j+oIWkzkznP0ethP6QMs+20K0Dbn9hW1Wa9wtEl0UjDLbQ3RjmWm/p3ahimwp/k19FxZ2djNtHQ0PFGHga2rEyStjQ/Tc1ldJDGvbT9rsmVyePcf68leYbL1fZFWXgIK8beXVG8HpNUVVXB51Oax4Ot6nBeeOEFNDQ0JLU01dXVcLlcyGHhUwoLC1FdXZ0so090DuQfyGuvTCAQQCgUgvdQ4U2OEjLZEQRBEIQ00FVmLJ/PRyY7HeGJJ57Aeeedh5KSkiO+/rGEmLEEQRAEIQ0kYOv035Gwbds2vPnmm7j66quTnxUVFSEajaKBrWTW1NSgqKgoWaampqZV/oG89sr4fL60reoAMtkRBEEQhLRgdVKvY7WKhdMxFixYgIKCAkycODH52ejRo+F0OrF48eLkZ+vXr8f27dtRWbnfbFpZWYk1a9agtlaFVFm0aBF8Ph+GDx+eLKPXcaDMgTrShZixhIPTXhycQ8TIsbjvDg0zRuMbGZpGwAo0kTz750Fa74B+yeP8Rbtpxcy3iM2u5vGOZjqnT2TS9tsj9AvDt0XlZ7/8CW0DhKOJ4WCaEi2+GI/j1N4449j2sbhldm3Jv4WOM7O5pcP1dgYrRONd1X//pORx1E9HWjyDjdEvVTp20nEkz7l2a5vXbKV7EnodpmliwYIFmDJlCvGN4/f7MXXqVMyaNQt5eXnw+Xz4yU9+gsrKSpx22mkAgHPPPRfDhw/H5ZdfjnvuuQfV1dW45ZZbMG3atKRO6LrrrsMf//hH3Hjjjbjqqqvw1ltv4dlnn8XLL7+clvs9gEx2BEEQBCENpNqpIAC8+eab2L59O6666qpWeffffz9sNhsmTZpEnAoewG6346WXXsL111+PyspKZGZmYsqUKZgzZ06yTEVFBV5++WXMnDkTDz74IEpLS/H444+ndds5IJMdQRAEQUgLCcuGhHXkapLEESw1n3vuubCsg5/o8Xgwb948zJs3r83z+/fvj1deeaXda4wbNw6rVq06/MYdRUSzIwiCIAhCj0ZWdoT0ovvk4XF7GMYXyiePyXRDXL9h1dUnjwfcshk0s+MxfyQ60FFA849j2FistDjTXmk+Y7je63Awa/eStLNF6b94vCjehqOFHkMOAFqKVV94aHNhY2HAEpoblcBA6lPKWXA8SWc9vwJdwmH4NRI6hgkDZifWHExREXYYmewIgiAIQhpIh2antyJmLEEQBEEQejSysiMIgiAIaaDzAmUxY3UUmewIqYVrbUxdv8GDVNEvAaLLaVVPOy+9aAtSC3tuPB6a/qwOqY/R6+rEc7ROoL5oXAE1lrhm52hhc1NtTezkwSSdUa36xcZcCMW81FwR0kIPxTJpni2XpjMPwx9Ru8h71OXs1+wcuSmqM+f2NsSMJQiCIAhCj0ZWdgRBEAQhDZidiG+1/3wxY3UUmewI6UVbGtdNWvvpouV34eijbydnYR1aFdXyucsAm8fNCnfNMr2xbhNJe7yaSSmTBSeM0DAOXYXBzXls231mtba1nmkxGs5wkbShFXU10+sU/u+nJC1vUfdFNDupQyY7giAIgpAGTNjEz06KEM2OIAiCIAg9GlnZEQRBEIQ0kLAMJKxOOBXsxLm9DZnsCN2Hzmxt5a7shbTBQ0C0ejZaPtfomKEQSdtzctSxk4ZXSDQ2drxNbnodIyMjeRzfsbPD9RwOh9pq7gzQLe/BfqpN+tZ4AIhnUHOFPaL6sM+qFnrhLtI5CUefRCcFygkxY3UY+Q8hCIIgCEKPRlZ2BEEQBCENmJYNZid2Y5myG6vDyGRHEARBENKAmLFSh0x2hJ4B1/uIhie1tKe34hqedsIXGA6my2lo6HBZPfSETdPkAIDhpn5qkEXzuwpH3z4qYaNj0P7BWtomplfK+FzzOXXiIJLnrqf6H514Fu2H9r0cUbhPJO73SBB6CjLZEQRBEIQ0YKJzO6okWlnHkcmOIAiCIKSBzjsVlBXsjiI9JQiCIAhCj0ZWdoRjklZaA1OEet0FMxYn6Va/qHQ9FdP6HI5mxJZJdTe63x1bVibJa/z6QJJ2htR1XV90+JJw5OfTD3L9JLmvsjB5nPcC1ejomiIASDTTtD6m7bvqaFlXNkl7a9Wxq57G8jJDHY/t1V5/i57n6NP52FiyXtFRZLIjCIIgCGnAhAETndHsiAPJjpLWaWEikcCtt96KiooKeL1eDBo0CL/97W9hab4DLMvC7NmzUVxcDK/Xi/Hjx2Pjxo1pbLUgCIIgdJ4DKzud+RM6Rlp76u6778YjjzyCP/7xj/j8889x991345577sEf/vCHZJl77rkHDz30EObPn4/ly5cjMzMTEyZMQDjc8aVaQRAEQRB6L2k1Y33wwQe48MILMXHiRADAgAED8Le//Q0ffvghgP2rOg888ABuueUWXHjhhQCAv/zlLygsLMQLL7yAyy67LG1tF9IM86Nj2Jj2o4s0PFy3wH272AaUJo+3X1RA8sr/spmk47uru6RN3Z5D6HAMh00ryp7TYcRHs6JRWq/2rOKDSkheUzkdL6ZdpUvfpD5szEjbP6T2fWcISYfzqRnB3aDuJ/416ivHUR8k6Xgu1RxF/coXUMa6GpLX/7fLSTp29kkqsXYTyWultWFaofbQx/fhnCccGZ13KigrOx0lrT11+umnY/HixdiwYQMA4JNPPsH777+P8847DwCwZcsWVFdXY/z48clz/H4/xo4di6VLlx60zkgkgkAgQP4EQRAEobthWkan/4SOkdaVnV/84hcIBAIYOnQo7HY7EokE7rjjDkyePBkAUF29/5dwYWEhOa+wsDCZx5k7dy5+85vfHN2GC4IgCIJwzJDWlZ1nn30WTz/9NBYuXIiPP/4YTz31FH7/+9/jqaeeOuI6b775ZjQ2Nib/qqqqurDFgiAIgtA1mF+ZsY70T5wKdpy0ruz8/Oc/xy9+8Yuk9mbEiBHYtm0b5s6diylTpqCoqAgAUFNTg+Li4uR5NTU1GDVq1EHrdLvdcLvdB80TOkE7vlHaw+aiMYm4D5b26rJnU98iVdeNSB6XzV9D62W+RXQNj9UJ9yBc/8Dbu/kOb/I4VkX1J5umU98uORsqksd5f1tJ8kymP+lJcM0O0ekcxljimKFQm3mOuhaSjnup352EJtPZduPJJK987gqStvdXuqxgYftmA3ejup+Eh369BofnkXQ4j/6j8uzT+sJJz7UPpmPJWPJpm20w+PefpkFqL54YIL50Uk3no57LZKejpLWngsEgbDxYnt0O09z/0ldUVKCoqAiLFy9O5gcCASxfvhyVlZUpbasgCIIgCMcmaV3ZueCCC3DHHXegvLwcJ5xwAlatWoX77rsPV111FQDAMAzMmDEDt99+OwYPHoyKigrceuutKCkpwUUXXZTOpguCIAhCp0jAQKITjgE7c25vI62TnT/84Q+49dZb8eMf/xi1tbUoKSnB//zP/2D27NnJMjfeeCNaWlpw7bXXoqGhAWeeeSZee+01eDyedmoWBEEQhO6NmLFSh2Hp7op7IIFAAH6/H+NwIRyG89AnCABa2/ZbaVfaO1fTGyRaqG8RrtHYesfpyeO8z+lQrD+e/mqJZ6h8R4jmDfz9ZyRtNmuaDe735VA+eNrRkWx8cjT9IKTu1RZmJlnWRle9Ou6zhmol3G+upk0QHyedguu9dvzPCJIOayGuDDYc8j6jHzSVas+R/W/x7KXpnC8jyePGAVQ7E6VNQp91VKfl2qc0SMFSqjHyvrKKpPXxwd/Vdsf7oTRSR6jN60nErRjewb/Q2NgIn893VK5x4P/Sb5aPhyfryNccws1x/Hrsm0e1rT0FiY0lCIIgCGkggc6ZokRO3nFksiMIgiAIaUDMWKlDJjvCQTFs9NeGzZfVZlm+XdVs1LxWH2Ip3HSpJfZYJjNbZbPleC07nkVNDRt/NYyky19Xy/zuZV/Qa4YjJM1NRrpZYMMfT6J5jfTLxXKrNpoudq8GvZ9YXJ1bN5SaHopeE7NVV2JF6DP2b6ZjNJKrzLKmm5lPh7Cxr53q3kevk7GH1ls/RJmuEtTrQiscLdQNg21HbfK4eexgeh32PpIWH66Ztj16qekqXXQ2mKcEAu040lOCIAiCIPRoZGVHEARBENKABQNmJzQ7lmw97zCysiMIgiAIaeCAGaszf4fLzp078cMf/hD5+fnwer0YMWIEPvroo2S+ZVmYPXs2iouL4fV6MX78eGzcuJHUUVdXh8mTJ8Pn8yEnJwdTp05Fc3MzKfPpp5/i61//OjweD8rKynDPPfccWSd1EbKyIxwUHr7AGlii8pz0BXNu3NXhegOTqefrRIG6TqPBQkuw3e793lHHu89gGgZWdsdZqq78Arrt2PfsRyS95a7TSbrP15R2wtjLdBV+qmlIBLQ225lWIkbbmPCq/HABzav9CW1D0Z9UG3tyKImjBe+zrNdoeJGmslHqeCB9bhbd9Q1ns3pWFvvGDPahAy/iV8cJL9rFstMxUHOR0uk4m2mb+P3oriAkxIPQUerr63HGGWfgrLPOwquvvoq+ffti48aNyM3NTZa555578NBDD+Gpp55KOvKdMGECPvvss6R/u8mTJ2P37t1YtGgRYrEYrrzySlx77bVYuHAhgP1b688991yMHz8e8+fPx5o1a3DVVVchJycH1157bVruXSY7giAIgpAGTMuAaR25KerAuYFAgHzeVozIu+++G2VlZViwYEHys4oKFbPPsiw88MADuOWWW3DhhRcCAP7yl7+gsLAQL7zwAi677DJ8/vnneO2117BixQqMGTMGwH4Hwd/5znfw+9//HiUlJXj66acRjUbx5z//GS6XCyeccAJWr16N++67L22THTFjCYIgCEIa6EzE8wN/AFBWVga/35/8mzt37kGv9+9//xtjxozBpZdeioKCApx00kl47LHHkvlbtmxBdXU1xo8fn/zM7/dj7NixWLp0KQBg6dKlyMnJSU50AGD8+PGw2WxYvnx5ssw3vvENuLRA0BMmTMD69etRX695V00hsrIjCIIgCMcwVVVVxIPywVZ1AGDz5s145JFHMGvWLPzyl7/EihUr8NOf/hQulwtTpkxBdXU1AKCwsJCcV1hYmMyrrq5GQUEByXc4HMjLyyNl9BUjvc7q6mpiNksVMtkRDo7BQh/UaeEXcjNo2Rzqptyob0ge7716LMkrm/wlSe9ZOTB5nD+c+t7Pn0J/AdRcNCR5PPivjbSeMX6SLvjH58njL+YMIXmZU4tI2l7fRNKBkPqiyMii/locdqqPCDtVOtJMNUe2RqrniGdqeh8mLIxl0aXsvZerX015f15G8sQXyqFpFUKBUfKUCi+y4bahJM/MoTot+z5VV0sp1dIYTJdlaI8mkcG0QDaa3vodKuqxa7Kcgf9bQ/Joi4SeQleZsXw+X4fCRZimiTFjxuDOO+8EAJx00klYu3Yt5s+fjylTphxxO44FxIwlCIIgCGnAhK3Tf4dDcXExhg8fTj4bNmwYtm/fDgAoKtr/Q7Cmhk62a2pqknlFRUWora0l+fF4HHV1daTMwerQr5FqZLIjCIIgCL2AM844A+vXryefbdiwAf379wewX6xcVFSExYsXJ/MDgQCWL1+Oysr9O2krKyvR0NCAlStXJsu89dZbME0TY8eOTZZ59913EYspz/CLFi3C8ccfnxYTFiCTHUEQBEFICwnL6PTf4TBz5kwsW7YMd955JzZt2oSFCxfi0UcfxbRp0wAAhmFgxowZuP322/Hvf/8ba9aswY9+9COUlJTgoosuArB/Jejb3/42rrnmGnz44Yf4z3/+g+nTp+Oyyy5DScl+FyX//d//DZfLhalTp2LdunX4+9//jgcffBCzZs3q0v47HESzIxwcrgsJKF2LkUHFb5ESqpdxbldalRtmPUfy/lE9mqT7j9iZPN5SRUVvkSeo5sVl35M8/s7MD0neG7V0aTb8fRXLa7Cxk+RFEnTY+zPCJB031W8Ah432QzRB26R/1Tg8VFkR56u1QXXdRBbr3wZab8KjjlsmnUryMv/BNDy9FJvb02aeGaHPlMc/szuVDidrC/3N13gi1WWFhqm6rGamBYrTfzaZVaquCPt2jWdTzU68kPrOcfrVdfaeQQWiedt3kDSJf2Uw/06i6Tpm6CrNTkc55ZRT8Pzzz+Pmm2/GnDlzUFFRgQceeACTJ09OlrnxxhvR0tKCa6+9Fg0NDTjzzDPx2muvJX3sAMDTTz+N6dOn45xzzoHNZsOkSZPw0EMPJfP9fj/eeOMNTJs2DaNHj0afPn0we/bstG07B2SyIwiCIAhpwepk1HPrCM49//zzcf7557eZbxgG5syZgzlz5rRZJi8vL+lAsC1GjhyJ995777Dbd7QQM5YgCIIgCD0aWdkRBEEQhDSQgIFEJ4J5dubc3oZMdgQAgGMQdQAV/3ILSW/+8XFtnlu2KEjSX/sglDxeuIvqTUoyqH+c47JV2ZG5NMbW3mgWSTfHlFbog3raHp+LajQytMBaURZky8Y0Dtzubdc0O2471eFkMMlG2KHym8JUy2SxehP6ZU2aF8mji6zuOpUfLKR5LHRTr4XrcnRsGdQXlBmkY9QKqXHX74Uqkle0NIekN89S/R/PYDHZInRsxbJVWW81LdvUh2ppbK62Y1rVTaD3tqdyFEkPu0tt641v3dZmPenC5qI+pzBMva+BodQfTNbfe68GzbQOX3fDzxc6hpixBEEQBEHo0cjKjiAIgiCkAbOTAuXOnNvbkMmOIAiCIKQBEwbMTuhuOnNub0MmO70Ye3Z28jhewOKq0BBW6POJ0hd4a2m8qH88M5+kx626InncL5tqdKImHXJNceW7YfHSESRvwumrSdplU/oY/ovGbaPaGpvRtq+RbS15JJ3tovej+9I5lJ8dl6bZcTrovQVjPBif9sUUovVws32orzq20+b1Wni8K+47R883PLTvbXE6PqyEGs+bryqjZU8MkHS8WT1Xm4ONq31Um+Js1i9Ci7p20fbH+lPNTrBexcoy7PRkd58QSe94QCm3Mv9+GsnL/tuRa2AMuxqXeh8dLoFLTibp9mQp8W+NIWnHoo+O+LqC0BYy2REEQRCENHAkXpD5+ULHkMmOIAiCIKQB0eykDukpQRAEQRB6NLKy04uo+vXpJB3LVLqAgS+EeHHCjvOUVqG4tJnkXfzF90h6UO6+NusJxqnGYU1tcfLYyKfilLX1xSTdP7s+eZztoH5IuO8crg3SKfDQ9nM/F9UhpWVyMO1PFFRrE41r8a4S9LeDk8XKirZo9840GTy8kS7vaUd+1PMxVJ/qehIAaL6Yxlnb9U3ViXY2lmw22sElC9SzCJdQ7U+Okz43Z47SrgR2Z5M8Vi0yatQHLcXMxMDKOtxUE+P1q3fQ56Xj235fH5J2vbYCHcV+wpDk8e5zaD0lb+wh6fgXG1X7SvvRvB00xlx7ZD3Te33nHA4mOhkbSwTKHUYmO4IgCIKQBqxO7sayZLLTYWSyIwiCIAhpINVRz3szMtk51jGY7MpSNo+W79ItqVnb6Tq6M6jSjo07SF7kHGoiqBhQmzzmW7X7MLNQOKG22MZZqIZ9YerG/2sFKkTEe3WDSV51Hd0OX5Fdp+qJ0FASpRn1aItQgm75PVT4iCxnNHkcN2n/JpggMNerwhA47GyrLr1VuPOVeWRvCw36EI6ybckx1cZQC82jxr0eBhvP8fFqC/POb1ITqHdkHUkXu5Q5yuOgpikeuiMwXbk8MKpy6DWZewHdPOnIYeaxPfQ5eurUM7YMWs/e0+j4yHRRc5lHM5/xd6zJ23F5JTf3Yad6d007NWMFhlM3DHXfVabu0ts/IHmOvvRcOLX3fNfuDrdPENKBTHYEQRAEIQ3IbqzUIZMdQRAEQUgDYsZKHTItFARBEAShRyMrOz0Y3+dU02A0sC3X9Q3J43gwSPL2/KSApEdmqbAPxR7qTr+ZbSd3anull1eXk7zGnX7ayOPUoSeb6hROKKom6fc+Hpo8/tqIrSRvVyiHpOsiyvU+D/nQx91C0s1xGlrAY1d6D9NGfw9wfVJtO9vUMzXtDwB4HFSjoRN00T4MxZQeoplpivZeR10I9JlPtRXHMnt+THVmDWNUH+YX7CV5TvZc9zYqHdeAvtT9gcdO+17XxOQOoWO/OkC1YmHNZYCNhYcoXkHrdQRVm5xMZ2OEado06a/yxqDSEe3ZR7e4+/pRHY4XbdP4g1NJOpyrruPdS8dS5r8/Junst5QGafsv6TgzqXQMFY9uTh7bc3JIXqKhoZ0WCgeQ2FipQyY7giAIgpAGxIyVOsSMJQiCIAhCj0ZWdgRBEAQhDcjKTuqQyc4xjmGjg93SXHmY67fQvDj1PaITHz+GpIf3pefmOZWuwW3QehJ22oZP60uSx00Bqi44d/SaNtvANS4ffTGApI1MdXNrPhpI8ry7aRvGTFqbPC72NJK8hhhtU6aD+U7RNDL8yyTCwlDkuTU/O0xDEoxTkUM4os6tD1EnPFlu2oaIFobC6aT+WepPof3fZz56DPWj6L31L1U6HRuLtxCMUf2MocVuqG2mfpiOy6ManpXrKpLH7nwaKiXSRDVctjr1HPPZ8K0fQsdDTJPaRIqonievHx2H/Jk3vqjeG7OIXqeB9UuxppHh+hjfX5fSNNqGRylBXLU59ws67hoGUd3Qrkn0HdTh4U9yN6h327Hoo3Za1LuQyU7qEDOWIAiCIAg9GlnZEQRBEIQ0ICs7qUMmO4IgCIKQBix0bvt4KzOk0CYy2TnGsUw63A2H0hcYLqoZafhvqsup/aayzw96ivkhcYZJWte1BOJU81IboT5BvvxCaQ98pdQnT4mngaR1HUYj09JwvyTDTqxKHof60XtzjqL6Aj0+VyDuIXk8FtZr64fT6/RT/n0GZlGtR1WI3qv+RcX994QT9PVqDKv7aw5SXUgwTPUnUU3fkwjSegw7bf/225Q/lPLbji2fOy2XUr86Ayt2kbRP84dTFcgheYFm+lzzfEo/VeGnPqY4FYPUM65upKoWXaMDAP6N6hk3DqL1RPNZPLSEKtu3lMZrK8pqIum120tI2lauacWc9BkPfZj6AjpafmwSzcqPlG8N9WuU8c8v2zyv8UeVJJ3/ThUtEFf91La3qd6HrOykDtHsCIIgCILQo5GVHUEQBEFIA7KykzpksiMIgiAIaUAmO6lDJjvHONzPDkkPHkDy6obTsvYM5buj7gTq92V2H6r9WNKs4lLxUbNjXylJX3bGsuRxC4s71ZygOougpq3JYP5u8iuo5iEQUeeG47QRRZlUD1EbUjF+drZQTYabxUnqX0h1OQVeVVeCiQf1vP3XURqeNXuLSZ6d+d2JxZWfkniU+izhONxK42DGmbWZqRKjueo6NTNoPKN+/7uepOP76L0eDRx9+9Br7qHaD0d+fvJ4zyh6b/Hd9NwTy5WGx+emOjKvk/qeKdXit7ls9BnXRej4bo6ocRnaQ/Oyq+gzDwxSHW6V0Da43fQ6ZkLdj9tJ8+ojVJPm9ND8uF21afDTVP9lfrwOqSa+sW2NDsf/F+rb56jpcgz+LpgHLycIDJnsCIIgCEIakJWd1CGTHUEQBEFIA5ZlwOrEhKUz5/Y2ZDeWIAiCIAg9GlnZOcYx7FT7YY06PnnsvIdqJXJbqMYhx6tiAm0ZUEbyrnzhOpIeNEr5zRjuryZ5ZxRsbrN9Dhv1Q7Iz5CfpMi/V5eicU7KBpHeH1bkNUap/+HJfPkmHdigtja0v1VlYJv01lJ/TTNL/qVcxf0pyaDyjbFfbcbQyXTS2V3OE+s7RNTtuL30WfDlaT9kcVJdgRtgzd6k2NA9gZYPUP4tjsHIUcziaDN1/EwAEvkd9NgXKtRazn1D95rbt8yaWS8dHUR/ql8ll1/yzWLTiAT5a78AMOt51uGYnZmp9yMZDsJSKomxlqg/7+OlY0XVCAPWttFfTjQFAwqTtj9RT/Zozotphfdh2DLleDdPo8O8/UjSRaDOvu2DC6JRTwc6c29uQyY4gCIIgpAHR7KQOMWMJgiAIgtCjkZWdYwy+bFv1zGCSPr9CLX8/v+FrJC/aRM0qdZlqmT3Rh24WtXtp+qTcHcnjmgjdyr3stRNJuu/pu5PH+R5qRuG/RGwZymSQYaemnZ3hHJIOxlX7fS5qmhrch5owqlyq/c0huv09EaN92PJ+X5L2nKZMa1t30Ly+BdRsoYcAyGP3muMOkXSDtvU4FKNmoQQzpegmLzNGf5NwMxy08BEmfcT48rZRJF1xk9oibM/JoW1oJwRB6DsnkXTtKTR/0N/VVml7E73v5gtp4ez3NiWPh970Gcnb9mR/ku6Xrfr7xNzdJM/Ntpc3J+hzbg/d/GjEaH8mcmm9RTnqGQ/00a373I2BS3Ol4GCuB7Y15tI2sHAoA29ZkTzuyTGP7COHkXTi08+PuK72TFVVs6kbBneDOi55OvUuGQ6GCJRTh0x2BEEQBCENiBkrdYgZSxAEQRDSwIGVnc78HQ633XYbDMMgf0OHKoex4XAY06ZNQ35+PrKysjBp0iTU1NSQOrZv346JEyciIyMDBQUF+PnPf454nK5wvvPOOzj55JPhdrtx3HHH4cknnzziPuoqZLIjCIIgCL2EE044Abt3707+vf/++8m8mTNn4sUXX8Rzzz2HJUuWYNeuXbjkkkuS+YlEAhMnTkQ0GsUHH3yAp556Ck8++SRmz56dLLNlyxZMnDgRZ511FlavXo0ZM2bg6quvxuuvv57S++SIGesQ2FxKBGHLzyN58d3VvHiXwLf52o8bkDzeeAXVkAzO20rSus7lW4O+IHlZdrpt+tlPRyePPX6qgYl/mUXSz29R9m8HlWTgtPPXkvS7y4cnj6ODqS1873oaDsAxRukaXGyb+oplQ0ja2ax+xYSL6S+JgQPps2gOKv1GfCfdAuw7jm13P53eUIt2rn8VFcG4z2fhAdr5ZZXlpFvR9a3TLqb1IFuhATSFVRtcmVTLFG1mwpy4aoMtxPQ97OfMrpvUc4xl07z+t9AQITqZ728k6QFNA2mBZZ8kD81TR5Csqu9QBUp+X+UeIf8xek3//9FGfZmr9GHDr32X5O2J0jG6MaDeDX3LOgC0xGifNbWobd/2MH2GcVoteVZRk35lbmqk47m2XrU/M4O+b00B6i7BHqLXteL0OfdUuEZH/44FADNK35uuIqQ9ql0/OJ7kFT6yPHlsWCaQol3rVifNWAdWdgIB6rLB7XbD7T64hs3hcKCoqKjV542NjXjiiSewcOFCnH322QCABQsWYNiwYVi2bBlOO+00vPHGG/jss8/w5ptvorCwEKNGjcJvf/tb3HTTTbjtttvgcrkwf/58VFRU4N577wUADBs2DO+//z7uv/9+TJgw4YjvtbPIyo4gCIIgpAELgGV14u+resrKyuD3+5N/c+fObfOaGzduRElJCQYOHIjJkydj+/btAICVK1ciFoth/PjxybJDhw5FeXk5li7dv7Fh6dKlGDFiBAoLC5NlJkyYgEAggHXr1iXL6HUcKHOgjnQhKzuCIAiCcAxTVVUFn0+tgra1qjN27Fg8+eSTOP7447F792785je/wde//nWsXbsW1dXVcLlcyGG7NAsLC1FdvX/lvLq6mkx0DuQfyGuvTCAQQCgUgtdLVzhThUx2BEEQBCENmDBgdIEHZZ/PRyY7bXHeeeclj0eOHImxY8eif//+ePbZZ9M2CUkVMtlhcB8QNWco3xj2MNUeWDaqW8h7om3NQ3s4BlWQ9N4zqT11zylK11LxArVnf15cTNuQr9zZV+RQd/pRV7DNssEwtZt7d9MXMFygjiO5tB+2NFEt04TTV6uyTOOwjprKEU3YD3oMALlDqd5nX7UKF5FfSG3UWzbRPisfVJs8bnBTLUSh5hsHAHKYzx5ot/PFeKqR4uju2sNxR5t5ANXw1CXoFwvfVaH72YlH2naJv79A21+WNuZDJqpH62DOXDb+4TSSHvyTZcnjRB3VOdnfXtnmNTf8iIZmKFlM2+DbpJ5d8L9OJXnV51AtkyNT9dnfN5xM8gbdTEM32B5VZfPcdKzvbKRhSmJh9azcXDvDfEwlNOHT7iDVFPHn5nAqsUeggfaDdwPz98R+fDuK1Rg+WnrA7ogtm/apeYQ+b3b+kvrV4T6n4plaWBUaqQOh28YmjxPhMHD7P4+oDYdLuv3s5OTkYMiQIdi0aRO+9a1vIRqNoqGhgazu1NTUJDU+RUVF+PDDD0kdB3Zr6WX4Dq6amhr4fL60TqhEsyMIgiAIvZDm5mZ8+eWXKC4uxujRo+F0OrF48eJk/vr167F9+3ZUVlYCACorK7FmzRrU1qofk4sWLYLP58Pw4cOTZfQ6DpQ5UEe6kMmOIAiCIKSBA04FO/N3OPy///f/sGTJEmzduhUffPABLr74YtjtdvzgBz+A3+/H1KlTMWvWLLz99ttYuXIlrrzySlRWVuK00/av+p577rkYPnw4Lr/8cnzyySd4/fXXccstt2DatGlJndB1112HzZs348Ybb8QXX3yBhx9+GM8++yxmzpzZ5f13OIgZSxAEQRDSwIFdVZ05/3DYsWMHfvCDH2Dfvn3o27cvzjzzTCxbtgx9++43199///2w2WyYNGkSIpEIJkyYgIcffjh5vt1ux0svvYTrr78elZWVyMzMxJQpUzBnzpxkmYqKCrz88suYOXMmHnzwQZSWluLxxx9P67ZzoBdNdux+H+zGfiNuorGxzXJbL6Lxa+yanMPuprPofs9tJemE5h/nsHxmxGjZ5lI2W8/UnD6wrCFTPyXpwIsqtlCWg/r54Noa3RdJxE7j+ISYG4avn6Wuk8+0P26Dtv/p1cr+fXx/Gs8oEKTG8gxNx5LvpfUO9FPNUd0+5QClzNdA8vKHtJC0R/ONUttINQF7g9TvDtfsZGj9dnLhDpJnN+i3y+6gEgVyvzo1LfS6u3V9EovP5XGxGEua9iPmpM+GjwFTa1LCTcuaYXodQ8t2tNCFXfc+mt56h9JAmE563yX/odfx/Ev5KRk8fTnaI3SB0ulUnUdvxrWLfiVFy9R1CgqoTitamkPSdc9qY+t720mez0ufcTiq3tV4BvPzspeKaRoylc5gWJ9aklfgpvqvXdlKG7SpLp/kNZa1r73qTTodncOJSxUfP4akm8rUc7TTrw+0DKDvlJHQxlqrzUpq7JudmX10c5555pl28z0eD+bNm4d58+a1WaZ///545ZVX2q1n3LhxWLVq1RG18WjRayY7giAIgtCdSLdAuTeRds3Ozp078cMf/hD5+fnwer0YMWIEPvroo2S+ZVmYPXs2iouL4fV6MX78eGzcuLGdGgVBEASh+5Pq2Fi9mbROdurr63HGGWfA6XTi1VdfxWeffYZ7770XubnKlHTPPffgoYcewvz587F8+XJkZmZiwoQJCIfD7dQsCIIgCN2bVAuUezNpNWPdfffdKCsrw4IFC5KfVVQonzOWZeGBBx7ALbfcggsvvBAA8Je//AWFhYV44YUXcNlll3X4Wo3fOh4O536bftYzyn+Io7CAlOO+GXR9hMVM7lYODaJjd6uT45u3drht4SFUIBMaTLUf/o9Uvc632tdDjCtWq151UapNqcimGphVe/qpa26hDqlMHw0Os6FB9dOpfbeRvIRB58w5mv8eh0G1HX181DfKrqWqDVuKqI395GFbSfq047aoelxUo7NoK3Xg0z9f3evxBVRnsS9M/Z80RKmOKBinscl0mqLU2O91anolZup3svhMCe2LKdNN8yLMR4/Npirzeuh4CIaYxgTawGTffZaHXsfWrK5j2WmDeYwoS2tSPJc+G661wSUqztrgKdQHj2GnL079IFWxg8nnCj6m7c28fV3yuPGSUSQvlkXHVvZO1cYvd1P/SD4fjX+WofVpQyF93o5d9Bkbmk4rnKDPqSFG/Ybkanq2Eh+t1+1g8boi/MtGAABHgXp2TWdQX2YJppvs89ya5PEXd59A8uwZVEto6jo59q4mtPfEDPWOGGW9jbSu7Pz73//GmDFjcOmll6KgoAAnnXQSHnvssWT+li1bUF1dTeJs+P1+jB07ts04G5FIBIFAgPwJgiAIQnejU3GxOrmTq7eR1snO5s2b8cgjj2Dw4MF4/fXXcf311+OnP/0pnnrqKQAq1sbB4mwcyOPMnTuXBEQrKys7ujchCIIgCEfA/glLZzQ76b6DY4e0TnZM08TJJ5+MO++8EyeddBKuvfZaXHPNNZg/f/4R13nzzTejsbEx+VdVVdWFLRYEQRAE4VgjrZqd4uLipIvpAwwbNgz/93//B0DF2qipqUGxFgOqpqYGo0aNOmidbrf7oBFfHUETDu6vBEC8di9JJ1yDSNoyDO2Yntt4IvVbk/P2lyrBdCywWl/7AK46qifI+ZDG8Sn8T0Py2NaH+u6I11A9yp0Fyh/Oj3eOJXkmEx2FIkpTkLuO6TUubiDpmOYjpilONS51Jq3XNFVdG/dQ7UQ0wJ5NqdJOeLdQDUNwMNU86Fqa4zL3kLxvlH9J0qGEKvveRzTembuY6n1qeXwjTWtT6KN+VOIWfa56PCy+M8JgPnn0eFdOpt/gZRMJdR39vINdx+bS9AZR5lenib7inlp1bpQOM6LRAYB4hmqTEaP3bbCYW0ZY5e/8BY1R1P+ZnSQd0+RhXrZA632eatL0tyb7b8vQHo7y0uSx9a1Skhf1Ut2T7tfIl0MdtATYr+WmL9WGiS0D6X33zaIatDy3epd1X08A4HLQNNfs2FwqbUZpe3sT8Vr1bnuf39NOSeDzR5XPJkeAjtFEJn0XbJrvqla7mLRnbjja/q7uamTreepI68rOGWecgfXr15PPNmzYgP799zvGq6ioQFFREYmzEQgEsHz58rTH2RAEQRCEzmB1wZ/QMdK6sjNz5kycfvrpuPPOO/G9730PH374IR599FE8+uijAADDMDBjxgzcfvvtGDx4MCoqKnDrrbeipKQEF110UTqbLgiCIAjCMUJaJzunnHIKnn/+edx8882YM2cOKioq8MADD2Dy5MnJMjfeeCNaWlpw7bXXoqGhAWeeeSZee+01eDyedmoWBEEQhO6NmLFSR9rDRZx//vk4//zz28w3DANz5swhgcaOhFiWHZZzvw2XTJOYlobrFkxNS+FspgPLW8v8MejSeK7RaUfD01JO/fVk7aJ6jnCJ8pdjldGymauprmXwOycljy8e9gnJe2M79UUTrFX1hk+kC6KWplMAAFeDav/a9X1InvcqqsnI9qjYUhku2ke1TFOSaFHtd9eTLDQznzYuTQMRYg6R3Daqh8jU4luVD6XCkH3N1P9QSzXt00SOOjfBNDoJk6YDITWa8jOpFoivMWe6lQ6jkcUI4xqehKZ7auU4jKWJ/5Bg+/GXwoVq3HGNA7tV2CKaXs3OdELRttOttD9bqV+msjkq7SjtR8u20e6OEK/alTz2VpeTvJYM6lspkaOckvJ/Fw4Pi6kUVOMwUE/r4f9sYpp+zcYGQHOEjufYOiqaClxycvJY9wUmtI2vQGmmmkz6HjuqmR+jAez91LD099qeQuNQZ21RYsfqMGmf7AiCIAhCr6SzIR9kZafDpD02liAIgiAIwtGk16zs2GJWclnZMVhtL49vpFuWuQt9aFvPWeQDuL+soR/4tT21e+iW9lZoZq3sVbtIVuh46kTR85Fq4/rbqCnquJeoCSlncf/kcVV/aopqCdJldM9u9fhdzG2/Zx9N2+Lq5mvH0F8Trga6HG+3a6YSFjKhII9u5U7kqLr2xOlWfrooDWQ4lEnMBvowXlx6Mkl/98wPk8d9vXT5mpuFWuLZtE0RZYrgpocMF90S7NDulf9Ca+8Xm8nMYTFmv9HNZYkY23rOd8aGVL7louPXCDGTl0Plx7NpWfdeZprSLZA8lAT7mWS0s5yeOGc0STuWKPcI8R07efEjR+uY8peYTTRGx+EXv1DP3LDRxnszIyQdKtfi8DVQ00gTe8RhzZ2DzU4fVDRETc5mAX3o9d9V29az/t5x1xW9Cd29AABkaSbz5kxqGnZto//ewnHVpw4XHQ/RiMoz46lbA+isF2RxKthxes1kRxAEQRC6EyJQTh1ixhIEQRAEoUcjKzuCIAiCkA4so3MiY1nZ6TC9ZrKTuSsMx1d3Wz9GhTDI5podJzOCaroFvt02WlFA0raIZgfe1PG2JQpzaNrLtmfXKf2Bq44uxm18cgxJezarY65NMbbSbbOx45VGIBygegLbp3RoBGgUDVrPFqquCWdr/eBlYRHyqH7Go7vQz2AahhYvSX+9UD2rhhjNQzbd4h5JqPbnu+g111UXkbSjLw3XEa9Ttv9wlPaDviUcALxO1eZgjLn/Z0KWqBZyw2ajGoxwiJ5rmYfxJaYVNZjmzHLT6zgaVRsSHlrWdLJraslWGh0mITG19yZzJ62nYRDVPQXPPCV5XPabD3A0SKxZ326+0ajCDFguejMxpuewaeEDEpl0jHq+pPcWLlTjxWBj0lZLy7IhAE+xKr/hT/S9HnLthxCA+jNpYGe7oXSTPOSKnb7WMJvUOxbPpto7SwuHYsVEs9MTETOWIAiCIAg9ml6zsiMIgiAI3QpxKpgyZLIjCIIgCGlAdmOljl4z2XHs2geHbb/NPHiqciW/+7FTSDnDZJoH7TjBZCKtWKaFZ+DhIRiGTfPfs5rqC7ygvnT0Ngz8Ay27/XGmP6lTduk9IRoWIV5I7dROzQ9Igr0z0RyaNnOVnoDrSfxrqd6npUTpQuIO2p8Ne6i+J6evcveOBO0zj4vpIzTHL6ZF7+34/jQkxKo6FYZgkI86DSrLoz5YNlZRv0aOJtWOsI9qaexZ1AdLq1AOGrpGBwACzUoLZOPaGtanlu7rg1+C/ZqznOo5cp8xYP584pmqrI1rE/h1dC0Qd/PC7tuuuaLxf0nHWfVptA/N9iNadA2H8Evj26DuvaUf7YcoezZuLXxInL8oDFtY84/kZv6R8qmGx7WD9ovTprRCbi2cBQDYRwwl6cSaL9psQ/xbVO+z90R1naL7j45GKlVUn0Gfqz+i7s3FwnwE+9HnqIdHiSdo3xvaozJSqNkRUkevmewIgiAIQrdDTFEpQSY7giAIgpAGxIyVOmSyIwiCIAjpQATKKaP3THbiJvCVTTykuccxmI8NsLgolqaBsIeoDb7qHBqLpf8SrV4b02AkqO8OQItnxHVCK9a1av4BPp9zHEnnu6j+JGuLsltXN/pInnM30wgMU44oYi6quwn3ZTGWNF2Iw0nvpXE48z1TpYaV5aRDzOpDNS9Bzb+Mwfyb1NXQ9q/LKVFtsNGye4NUw1OSFUgec11NlpNqSgr6Bki6Ybvyw2TjPleYc5SEVreNffM0NFK/Rpam9+DSD661IX3BNTqtfs1pPkLizLcSuxBxRcL8klh2Vq+WbcTb1ugANJZaqzhZLJ3BQsqlg8L5K5LHO25kuj2Lvgsxt3oWXIpn0leqXQwWKyteQR3BmEQkRc/dNJnGjXM2n66Ksq+w8Am0XtsO9QAMB703K051RN0dPg6z3OpdjkTpvSWyaMdkVGl+riL0uzyqxylrL9CbcMzSeyY7giAIgtCtMNB6Z8Dhni90BJnsCIIgCEI6EDNWypA9doIgCIIg9Gh6z8qO0w7Y9t9uzN+2HgLMrq5rb5zNtGjwFGobt58wJHmcWLehw01rpedpx0eIt4jGeXIyDYnp0JY1P8kmebFCep14i4rVYzFth+mjPivsWnygQYV7Sd76faUk7dIkMJaD1hvNB00HNTs7jwfF/J0EYqq9LqbZyfMGSTqsxcYKxKi2isesymAanppCde+OOPOVwvQyYU0nEGNlDXY7RDvUvhsY4oeHD1En8z8UjWu6J67RYdcx2tHhcPRzHbR74Wqk6YxaVXjviPaFLOG8drNTgq5VsbN7Cw6h40F/jvzVTLjo03E2q8Kmh/6WtGWyWFlMpxWNtf11bBvIvny07ymnk72rzLdSU44ao80XjyZ5mc8ta/Oa3YLTR5Gk5ad9qL9Tdv7dnUXLhvuofsneyr6X+uoPOYWmIVnZSRm9Z7IjCIIgCN0JiXqeMsSMJQiCIAhCj0ZWdgRBEAQhDVjW/r/OnC90jF4z2WkaWQyH8yvthhavyeY4hHhCg7tfSLAYKnvHKkFKbtuuclpziDg+OrZl1PdMnY+m+2j+cKI5zHcL1+WE1OPn+g0jxpZHM5T9e1ejn2TZg7QfXAFNb8J8t0S4LkdbhjXYs3Bup1qbur7Kb025v4HktedLJ26x2EcsZhXP9+QpLZa5nuqeWgYy3zRan7bywcL7VB9APCwVH1x6mi9Vt4phpZXlOhxeVu//Q/jDsbRuYl1E8gAaU85GJSStfNEUrOQ+p9JLv0VUg/bFSPpOGY6OPwuHHuqNxcaKe6gfGKefOivyutQ7xseDg+lR3A7VyVyD1hJlHa75EmscSL/yqXeq7seWC6ivKn9eHUnr997SQIMXOjOoZkeP05dw8edoHfz4aCOanZQhZixBEARBEHo0vWZlRxAEQRC6FSJQThky2REEQRCENGBYnbOaSWSLjtNrJjuZO1vgsO/XCth96rZ5nKdYhHaJTbOV7xvD/Fk4qR3dcekelfhz++3h8bA6SsnvPiDpxssrab1a84fc8QXJ2zD7eJK2azZt9xfU3s1/MIS1PnNmU/9CiQymtQmpkyN+WhHXSJkJzZLKrmmjYbSIpiEcp8+puoVqa0qzmSMYjaaom6TjibZ96cRZfB03Gy/hiK6PaEcHANBYU8yHkMl89NgMPR4TLRuPsvZqaaIvwUHkYHo2j4XF+1+TPJhUboJYFk27G9W70VRG28d99HhfWI7uROKzTSRt3zeWpE2X9n1RQHU2VibttJZMNR5ce2k/2OvpmI15abp2t9L8WV46zjLzaCe67Fo+G2dcw6M/V9+2jusD04Xta8OTx7E+VHfjdrDva+3ddWVR/0gW0weaGepc08UUHLqekQeuO5qkUbNz11134eabb8YNN9yABx54AAAQDofxs5/9DM888wwikQgmTJiAhx9+GIWFhcnztm/fjuuvvx5vv/02srKyMGXKFMydOxcOhxrP77zzDmbNmoV169ahrKwMt9xyC6644oojb2wXIJodQRAEQehFrFixAn/6058wcuRI8vnMmTPx4osv4rnnnsOSJUuwa9cuXHLJJcn8RCKBiRMnIhqN4oMPPsBTTz2FJ598ErNnz06W2bJlCyZOnIizzjoLq1evxowZM3D11Vfj9ddfT9n9HQyZ7AiCIAhCOjig2enMH4BAIED+IpFIm5dsbm7G5MmT8dhjjyE3Nzf5eWNjI5544gncd999OPvsszF69GgsWLAAH3zwAZYt2+9p+4033sBnn32Gv/71rxg1ahTOO+88/Pa3v8W8efMQje5fWZs/fz4qKipw7733YtiwYZg+fTq++93v4v777z+KHXloeo0Zq3FINuyu/VuZPR4VzyCRoPM9h6vtbbGe3bS7zHy6pL2nXplScvuXkbz4tqrDa3AHyf/PbpLePKUkeZyzmC6jD5rBXMNXfi15uO08msVXwkteVvdecwndDmq56dK4PaL61NH2O7f/XH37PjPtGKMCJN0cUeanbCetuCiziaT1rehkyR9AzKT9ku2mz9Ghld/FTEY2Fp7DrS2dR5qpeaz1ErNqU2tvA+zeNdMVd4NvsjgUpm4aZPZHizvi0Jb9eSgJHsKCbzfXiVOrJ1wfKVNQVt4wkueuZ3vRuxvsYZS+Q8fL9nNVR/CwKtxUYmmm7biPPSceWoI9V0eduk6MendAi0k3iTsK1bnZHvouRFjYCUtzR9AwiD5UZo3sFuwcn5M8dvloeBxuoguEtHeOW41Z2u5V4zDhpnZZW1Trl1gK1wC6yIxVVkb/3/z617/GbbfddtBTpk2bhokTJ2L8+PG4/fbbk5+vXLkSsVgM48ePT342dOhQlJeXY+nSpTjttNOwdOlSjBgxgpi1JkyYgOuvvx7r1q3DSSedhKVLl5I6DpSZMWNGu7eiryB1lPnz56OgoKBDZXvNZEcQBEEQeiJVVVXwaT7X3G73Qcs988wz+Pjjj7FixYpWedXV1XC5XMjJySGfFxYWorq6OllGn+gcyD+Q116ZQCCAUCgEr5f9UvqKF154Ad/73vfazOcsXLgQzc3NMtkRBEEQhG5NF63s+Hw+Mtk5GFVVVbjhhhuwaNEieDyedsumi4ceeqjDk5d//OMfh1W3aHYEQRAEIR1YXfDXQVauXIna2lqcfPLJcDgccDgcWLJkCR566CE4HA4UFhYiGo2ioaGBnFdTU4OioiIAQFFREWpqalrlH8hrr4zP52t31ebtt99GXl5eh+/n1VdfRb9+/Tpcvtes7AT7GLC799uudU2Gk21ljETb7pJgP6o9yGD6DbumXTHz2Sx7G6vsMEJEtEfw+L4kTbb5ZjJn8HuoW3ws/SR52C9rDMkK9Kcu5/1LlCbD90U+ydvyXZq2h9R2UcugmpdEkNrKSQiLOH1zIyFaNhJW6cIsqtHJc9OtuXURpSva00SVCRX+fW2WBQCnTY2JHD+tN8zGB3Hjn0W1E9yNgRlTfcG34HPth67TsVrpcPgecf0i7NuPh+cg29jb32LLQ0LQemly609PUJdgsjfPXjqW6Gjpfrhe/pCkh64fmDzecC1dnjcL6HbnnEI1LjNcdNt0wqS/LZvDtF+CpZoriB107HO9T1TT5YSY9ifG3BgYmgYl6u/+jlki2gBxONoPLRLT+oGH7zETbf+WN9nXvKlprcx499+efyScc845WLNmDfnsyiuvxNChQ3HTTTehrKwMTqcTixcvxqRJkwAA69evx/bt21FZud/FSWVlJe644w7U1tYmV2AWLVoEn8+H4cOHJ8u88sor5DqLFi1K1tEW3/zmNw/rfs4888zDKt9rJjuCIAiC0K1IoQfl7OxsnHjiieSzzMxM5OfnJz+fOnUqZs2ahby8PPh8PvzkJz9BZWUlTjvtNADAueeei+HDh+Pyyy/HPffcg+rqatxyyy2YNm1aUid03XXX4Y9//CNuvPFGXHXVVXjrrbfw7LPP4uWXXz7s26utrUVtbS1Mk05A+Zb5jiCTHUEQBEFIA93Ng/L9998Pm82GSZMmEaeCB7Db7XjppZdw/fXXo7KyEpmZmZgyZQrmzJmTLFNRUYGXX34ZM2fOxIMPPojS0lI8/vjjmDBhQofbsXLlSkyZMgWff/55ckepYRiwLAuGYSCROPxgwkc02VmwYAG+//3vIyMj49CFBUEQBEHodrzzzjsk7fF4MG/ePMybN6/Nc/r379/KTMUZN24cVq1adcTtuuqqqzBkyBA88cQTKCwshGF0YvXrK45osvOLX/wCN9xwAy699FJMnToVp59+eqcbcrSJ5AO2rwTouZodmPtiCISpgCorS/lgsYWoDZhrJ9yajb7qXCq06vfxYTdZYWjXZVqfxgpq2y+5W4WTqJ1Kn0veE1w4pHAs+oiWZfmBSaclj3kIiJxNLESBQ/Nx00TzvDuZy/wslc/9vsSZy/x+RfXJY4+d6qcK3FTDUx1UPo/6ZjSTvECU7kSIsnAR0YRqo501KtNNNRqNQTVe7EzDxUOR6JIkHgKiPf9OXIRoxtvZV8DGJA9LoWukDOb+hv9K1LvFYmEoWgkjbZoPIVZvJBfHNPFNm5PHGbuKSF6Lg+pukKf8wuR5qN6rJUbLOpn/p1CL2i5sOpm2LU6fa7hZ1cX9MEX30u8wW7YasyYLb2HPow8nUVePVGNnW50TbjW4XLb2ly7MKvWD23LSsvYIvdeEFnqCj2fyXvSScBHdmc2bN+P//u//cNxxx3VZnUe0G2vnzp146qmnsHfvXowbNw5Dhw7F3XffndxnLwiCIAiCcCScc845+OSTTw5d8DA4opUdh8OBiy++GBdffDFqamrw17/+FU899RRuvfVWfPvb38bUqVNxwQUXwGaTne2CIAiCcDAMdFKz02Ut6V48/vjjmDJlCtauXYsTTzwRTrbS+V//9V+HXWenBcqFhYU488wzsWHDBmzYsAFr1qzBlClTkJubiwULFmDcuHGdvYQgCIIgCL2EpUuX4j//+Q9effXVVnkpFSgD+50E/e///i8WLFiAzZs346KLLsJLL72E8ePHo6WlBXPmzMGUKVOwbVvbOpFUkrUdsH9l5rZGqc/tzCbszaCaDI9TCRBipVT7wfF7lb6n/4XrSF71Xe2caBz5ChjXuez6edv6qebvn0bSWX9f1kbJ1jjC6kI26j4E1WPpMBrwkvI3Y9mo2/IMFl8s1FfTerDRaM+kz2Jvk/Ib9LX8nbQNYerXKMupzm1mWgmT6Vq2rish6QEn7EJbNIXp/eg6nSjzwWNjY0v3pcP97Bjs553uIyTBNDrcJw/3edNeWUeLqstFZU6Is/0GpredmFv83rRbb6W96o4BmI6Q0me3kPTnd9Kxo+soeUw2lz1E0tsaqV7GalK/Xu00XBtMB9OfZKk0H3fwUtGUR/tO42V3XEXjmBX//gOkGqMPVQgmslW/mXysM7zVKr/fokaSt/7qtj0Km0xqpb9+KV0tSeHW82OJn/zkJ/jhD3+IW2+9tVXoiSPliCY7F1xwAV5//XUMGTIE11xzDX70ox8Rz4eZmZn42c9+ht/97ndd0khBEARB6HGIQPmg7Nu3DzNnzuyyiQ5whJOdgoICLFmypF2PiH379sWWLVvazBcEQRAEQeBccsklePvttzFo0KAuq/OIJjtPPPEEFi9ejF/+8pcH9W745z//GYZhoH///l3SSEEQBEHoccjKzkEZMmQIbr75Zrz//vsYMWJEK4HyT3/608Ou84gmO3PmzMFvfvMbjBkzBsXFxV3i8OdoEyoA7F/JLfI1HY6DCQwcmW0Ln8I7aVSfvEF19BpR9UA2NfYhefaLhpK094XlbV7HsNH+tGnOG41sKoBIsOC1LcOUfd5dRQ3TLcVU+5HVjv8eDo8XpPM/cwMkvXDbt5PHmTW0P5vY/Nep6UbiLJTXyCLqyuDLetX/7+2gM/6xJVQbFk6oZ9HAfCc1tdBO81U0kHTMVA5mXMyfj8dJffJke5Q+KRCi9QaDVN9jaToci/V3e3F8+BeanfnkSYTafo3bezUNNtTjXqbD0c9lmiLu06SVHx6S2XbWsUZ8126StrvpMnuOR+lyMhxUc8b9O7VC8+/CNSWOIPOflKGEcx43FdHF6tj4LlRjdC/T7PAYUekgVuwnad3PTZzF+aquoWX7f67ez2gBE5252HdauO1gb/pYT+Vw7W4elLsLjz/+OLKysrBkyRIsWbKE5BmGkbrJziOPPIInn3wSl19++ZGcLgiCIAiCcFCOhgTmiCY70Wj0mPCaLAiCIAjdFjFjpYwj2vN89dVXY+HChV3dFkEQBEHoPVhd8NdDmDVrFlpaWg5d8Ctuvvlm1NXVHbrgVxzRyk44HMajjz6KN998EyNHjmwlHrrvvvuOpNqjSrggAdtXsZbcmg4jHKdtb45QY3lehopv42faDo67HS3QtrOozf24F7QE029YTEthZCjNydYrK0heLJuW1f2qJDz0TSj6I41/daTvie/9viT92OcjSLrig33J483fpzonG4ubFC5QrfBWUO0Pj1kV0+z3el8DQIL5m9gRULb9YJg+01iEDvso0ziUDq1KHjfF2varA1BNj43lmSyeEYKq/RaLm4X2YgCxahIRpj3Qf7JwJzdBeq+6jT9OpUytdDdEl3OIeEFcw9PmRQEYdq0fjsA5WHciEab9m+1S+hgXG+w8nluM6VGg9T/XU3EfVDotLWyM+qlWKNervsMaWuhDDxem/7+lc/te+oGlfBfFgvT7GSHaZ9snqHHZZyXNc9TTMRvX/PfwLz8SG+sQvn2Eo8ODDz6Im2++GZmZmYcuDGDevHm45ppriNub9jiiyc6nn36KUaNGAQDWrl1L8o4FsbIgCIIgpBsRKCssy8KQIUM6PIc4nFUg4AgnO2+//faRnCYIgiAIwgHEg3KSBQsWHPY5h+N0sBtsPEwNtpwobBn71/sdmrnBaaPrxRkuto1T24acn0FnkjY2rXZpdfXx0LJ9R9P0puuVwLvPJ0GS59hVT9LxrWpbdekdtego9hF0u7vlpI/7SE0Iwe/TpeXIr+iWz91nK9vaoMerSN7Oi8toG4apfonF6DJ0kJkY9W3ePOSDvtUcAFo0c6SLmbxCLbS97hraL59nFSWPjy+tIXm7WVgKS2tHrpeGA2isp8uxGTvV/QWzqFzO4abPQh9b3GU+DzWhmy55eAh7M72OjVo4aD129jPRrYUIaWBbltk2dd0MZzjbd2Ngz1fLzvHaPe2W7e64WEgTfVw2x6h5NG7RZ5HgoTwa1Pjw0mGHwCDW33oIEfb/riCPmoMzHG1vUw8xVxvpMDFG+1M3HYZm/s1cS03QzSMjJO3aqkx4TQNovYMX0O/K9der6xjMnJ65SzP/R1I4gRCBcpIpU6Yc1folLLkgCIIgCD2aXrOyIwiCIAjdCdHspA6Z7AiCIAhCOhAzVsroNZOdjMwo7BmtbbGZznZEDADqQ2qrZraP2ovDcaZj0Oz1u4JU28E1PD+Y9kby+O8PnEsvOoxqSupGFiePj7thGclr/v5pJO1qVHZ273++IHlWjBmqdQxm0WTb4fXr5HxEBQV5n1CtTcPx6g3cfT7V6ASL6duZ+4rStew5hV4zlku1B/q2793bqJ2fb6mN7VZ9aJbQvrc30/ba6WOFzat0DfwZc41X9TalP/H0oZodbxatOO7VtggzvUaMbWG2YtrzaMfNPUC3jBsxptlh+gOHJg8L92HflF7a/4a2Pd5iYTL4VvT2NlBwjdGmGYOTxwN+eWxrdpzMhUBDRI1D/s43ROkYjTCXCAm/qsu+ifY3/wUf17ZkO7KoDodrCT12le910rINTu72IvWuAByN9D2xa++y7poCALLW0G32wZPUO2fsZOE4mmj/O5oKVIL1Z1jzkJEIH6rFwrFIr5nsCIIgCEK3opNmLFnZ6Tgy2REEQRCEdCBmrJQhu7EEQRAEQejR9JqVHbvdhN2+3z6th3Lg4QACYWb31fA5qW3ZxqbV67YpN+cGc//f3LeRXieqrhM8j2pTEgk6Bx3Xf3PyeMkDVKMz5M+0XvOzTVo93Oc883+i6XQMGxVd8JAVll3l7/lGEckrfJFFqLVUSAtXE71m7sW7SNr+kh5OgvZ93eYSkg72U32a+yVtb7Aoi6SdQ5tVrcy3SMWYbSS9eckAkvZ71XPm+gc9fAgAhIuUdsJgZePsOUZLVDuMANN7ZTB/J7qmh+ljHC3MX0uG6mODh3Vgv/xI2AGus7Fz/YatzaJMuoSEdh2+LG9nmp3cUceuTqflUvr+xWJ0POR5VJr7foqaVIeTkUG/T5rqlYYnTCVpiPu4PxzVyZmZVGTSx0u1KrqW0OOg74LFx9aA/uqaW+l7crQwP/mMpHNHVSaP942iZbl/HLdH3U8on75T2yfT0Do2/dZ5CBb930D7ErmuRVZ2UkavmewIgiAIQndCtp6nDjFjCYIgCILQo5HJjiAIgiAIPZpeY8ay20zip+UAjSGqE+HxjHQfIbu92SSvL7ON69qKH45aTrL+sWlUm20r9DWRdCTB/feoOWnOZ2x+un03Teu6HK7R4ZD89g3V0Wx1b01nM/8hEwpIuu8LmiaqH23v3o9KSTpvgFqHDVRQQ7qXhQHLrGrbmYuDSicQrFG+Opz9qN7BY6eG/3gmXQvO0HwvxZjOgqdbtviTx66yZpIXCVLNBqKqL6xs2gaDx6UKaOeytWrTSdO2sOoXR5D72aHVmlq1FnetFKcfGJqvn1ZaINYmI6T1C+vPbKYpqW9SPpDycnJIXqKhAd2Znd+i75SD6bLC2rvLY7vxcReNM186Wh/H6NcQ7EF6HT00XLObfofZ8tq2bVgsppwRpfWmSqfTHg3HqzY6qSQRRQ98QNKbRp+cPDaYz6NoLu1/Xdvm3tu27s08yP+Jo4ZodlJGr5nsCIIgCEJ3QjQ7qaPbmLHuuusuGIaBGTNmJD8Lh8OYNm0a8vPzkZWVhUmTJqGmpqbtSgRBEARBEBjdYrKzYsUK/OlPf8LIkSPJ5zNnzsSLL76I5557DkuWLMGuXbtwySWXpKmVgiAIgtDFWJ34EzpM2s1Yzc3NmDx5Mh577DHcfvvtyc8bGxvxxBNPYOHChTj77LMBAAsWLMCwYcOwbNkynHbaaW1VeVDssGD/as1vT1D5ZHE7qJ23tLiOpHfsUn5g6ltozCq7QW27A8qU/5ANzYUkb+JA6ktCZ0+E+oh5b+kJJD3wtH3J4z6fUF1IopH66DmkTqcNLLP9Nydrl9IbFN5MDemf/4zea3OxmkO3lLP25NFYZJFdmt6Ax/+hoYTg0qRNfVbTfmgcTEUOBR+r60ZmUs1IhoO24dTTaQyxLxuUkxNdXwIAsRDT4fhUvwy4jda7+yyq8cqs1vRfFzCfNjGmmdLWp21xppfhP1H0fmsnRhUAxLK0wkwnZAvSNtiiqjKD+dWJ59D227Op/xadUIT2mU3z51P1P8NJXsndVJPR3bExf1p1ITVeijKpFq8uTMdSnD1zPa6Zkw5vxLPpg7W86oH4fDQmG9cK6b6iWukWuRarG5Ctue3Kf4yOB5uLxhNzutX7F4/S/jzu3vUk/flvjkseG7wbdM0OzzyaiGYnZaR9ZWfatGmYOHEixo8fTz5fuXIlYrEY+Xzo0KEoLy/H0qVL26wvEokgEAiQP0EQBEEQei9pXdl55pln8PHHH2PFihWt8qqrq+FyuZDDdmsUFhaiurq6zTrnzp2L3/zmN13dVEEQBEHoUkSgnDrStrJTVVWFG264AU8//TQ8nrZDNBwuN998MxobG5N/VVVVXVa3IAiCIHQZndHriG7nsEjbys7KlStRW1uLk09WfhISiQTeffdd/PGPf8Trr7+OaDSKhoYGsrpTU1ODoqKig9S4H7fbDbfb3erzffWZsEX2T6q8mUpbke2hjkjczBcGwmo+aGP27kCETtKG5inHMMUeaj5riFEBSktc2Z6XbxlA8gb+H7XB7/utZldvXEfbd4QanVYcoh73ix+qRP8yklf8HrX7xz3qDXSE6Hy68H3qo+eL/6fOLXqdag1aium5MU3yUDecanQsFttr51War5wvaaAh5xAqQIklqK2/tlrznZNNdTgOD9WmxOvUGFh/E2uTScdWY4t63TLX0bETySFJxAvUdcwYi19UT19bS9PecN85HP2XYIL564GDxfbyqrTBylqsTYmI6kODxcIKB6nOwlGl7j1yEh0Pm/82iqQH/mA1uhP2FuZ3yUnHbNilfX8wXzk8ZlVDC/1OiOjh0NhXmC3CtDXa+MhyMx1cnI6PDEfbeioeZy0VGHb2vl03lqRzNrXd3o13jSbpbE998jiwO4fkfTnzeJJ2K+kj7FTGB0eDapMZTmVwLCFVpG2yc84552DNmjXksyuvvBJDhw7FTTfdhLKyMjidTixevBiTJk0CAKxfvx7bt29HZWXlwaoUBEEQhGMGMWOljrRNdrKzs3HiiSeSzzIzM5Gfn5/8fOrUqZg1axby8vLg8/nwk5/8BJWVlYe9E0sQBEEQuh2yGytlpH03Vnvcf//9OP/88zFp0iR84xvfQFFREf75z3+mu1mCIAiCcMzxyCOPYOTIkfD5fPD5fKisrMSrr76azO+II9/t27dj4sSJyMjIQEFBAX7+858jHqfyj3feeQcnn3wy3G43jjvuODz55JOpuL12SbufHZ133nmHpD0eD+bNm4d58+Z1um4zbge+8mvRUqfEH+EManNvcFI7uj2gaRH60Tpddqr9qItoMX+cNFjTmx9T3zk2zS8J95MBUD1Hoq4e3Yn4Nir69u3ZR9KxscOSx8FiqteoH5VH0oOvUG4EuF4j73X6LPaMVj9jzCo6T8/cTX/i/M+J7yWPVzb2p+1n8a2+bMgnabtXvbjZGdS4X7eP+kRyF6jnzD2WxJgflURQvW6RXFrW4j5vAqqs6aQaGNPDymqxsVr5D2HafxIayWQtZnI1Q2t+65hKbZ9rGfTZuGvo10xkgNIyeVk8I67v6Q44hg5OHpsu5l/ITdtfkK1862Q5qZYmysZdFtMLRovUd1Eszr4D8unD8WWpcem00TZksuvuCCgNWuitviQvPpCea/Oqd84M8e+lrqFh8qkkbbL/Qs7XW+/OTeaVUwdEgQb1ncvjzUWz2BjVY9Mx3ZAeUy4RTqHvoRSv7JSWluKuu+7C4MGDYVkWnnrqKVx44YVYtWoVTjjhBMycORMvv/wynnvuOfj9fkyfPh2XXHIJ/vOf/wDYr6udOHEiioqK8MEHH2D37t340Y9+BKfTiTvvvBMAsGXLFkycOBHXXXcdnn76aSxevBhXX301iouLMWHChE7cbOfoVpMdQRAEQegtpFqzc8EFF5D0HXfcgUceeQTLli1DaWnpIR35vvHGG/jss8/w5ptvorCwEKNGjcJvf/tb3HTTTbjtttvgcrkwf/58VFRU4N577wUADBs2DO+//z7uv//+tE52urUZSxAEQRB6LF209Zw70o1E6KrhwUgkEnjmmWfQ0tKCysrKDjnyXbp0KUaMGIHCQuU1f8KECQgEAli3bl2yDHcSPGHChHadAaeCXrOyY9/rhM2zf5k44VNLt/EmumxueuhSqL6j2WTL/vEEnSvatDXFxdsHk7xBg3eTdEtU7S2t2UC3Rvd9f9lB76HbYtGfF1suUsvxQ56iW/CDpXQ/bvXM05PHea/TevaeQ19Y3TW/ZzU1ceWtomE+3t6jtp328dAtvx473do6KIea4XY5fcnj5ggdHxnZtE2RiHqF4nuYCbQPNYHZtG3scRd79SJ0LLn2aVthnTyMAz1VDy3gYJaHqJ+mDc06Yg/R8WyLMdPgTnUco5EvYLKoGQ6ti0MFNC+S13aIglAT22MdSv1XEt8K3XTpKSTdcJxm/nDTsZORQceDx6G+P5pjdOzsbKIPIxii+YlGlfY20mcTzKHpsBaCI2CnJi+znZghmdX0HQuexAbT8EHqeOXaNuvpSgLH0fERulN9Jwy87SOSl+mh7c3NVAM+wb6fDdYNjUHVT2E7fVcztqlxZzFz7rFAWRl1B/LrX/8at91220HLrlmzBpWVlQiHw8jKysLzzz+P4cOHY/Xq1Yd05FtdXU0mOgfyD+S1VyYQCCAUCsHrZXGAUkSvmewIgiAIQreiizQ7VVVV8PnUj7SD+Zo7wPHHH4/Vq1ejsbER//jHPzBlyhQsWbKkE404NpDJjiAIgiCkga7S7BzYXdURXC4Xjjtuf1DU0aNHY8WKFXjwwQfx/e9//5COfIuKivDhhx+S+g7s1tLL8B1cNTU18Pl8aVvVAUSzIwiCIAi9FtM0EYlEMHr06KQj3wNwR76VlZVYs2YNamtVtIBFixbB5/Nh+PDhyTJ6HQfKpNsZcK9Z2XG0GLDH9xtxTbdmg8+kWy8tZve1NHN+YC/ddmw0U1v/Hr+yydsctN4tQarL6V+kdCInj/6S5LUwm2mioQHdGb5F1dZX6Rh2fYPqFIqWU/1M9GQ1BBMe2ve+D6kWoXGU0ksYJv05tOEquqUdO9WW1PICqufpy9z2hxP0NYjG1XPNy6AuBKKs7K6mHK3BVM+R/X4GSTcM1bQJGWzcsS3NUe12bFH6m4SHhLAcqt9idIjCTW+dbLHlW375uS2l6DD6FnfTyzQ6NvbTNdS2O37eD6nArBxB0nXDWJiSLNV+3S0BAAT2UA3aJk1Lk5NFx04rzV+U9YM2JLgLATv7rok6lb6nhfVvc4iaMBwr1YP1/IDq07JZm6or1cArWIkuI3buGHWcQa/pKabbycPZqv21V40heZZFB7RHC4XBtUoNQRaOY7f2rDz0/QsVqQ43Qykcgyneen7zzTfjvPPOQ3l5OZqamrBw4UK88847eP311+H3+w/pyPfcc8/F8OHDcfnll+Oee+5BdXU1brnlFkybNi1pOrvuuuvwxz/+ETfeeCOuuuoqvPXWW3j22Wfx8ssvd+JGO0+vmewIgiAIQnci1VvPa2tr8aMf/Qi7d++G3+/HyJEj8frrr+Nb3/oWgP2OfG02GyZNmoRIJIIJEybg4YcfTp5vt9vx0ksv4frrr0dlZSUyMzMxZcoUzJkzJ1mmoqICL7/8MmbOnIkHH3wQpaWlePzxx9O67RyQyY4gCIIg9AqeeOKJdvM74si3f//+eOWVV9qtZ9y4cVi1atURtfFoIZMdQRAEQUgHEhsrZfSayU40z4TtKy1Bxk5lk28pY6PForZxV72yA9v2UucirgZ+FWVHj2ewnHEspILmOj4QpTb2fX+l+pP881tdqFsRungsSSe06BYFq6mvGcvOdDlbVf/X0biwyN5O0w3tvNgrvncfSd+0Szm1aoxS2z3X6HB/KA67stnzMAmcLJ+6P6+Lanai36ZjyahXg8Lkeg0eukFL8lASvKylSQx4ay32hut6n1ahJFjZhFe7LmuCjYWLsOvuZnhoCSqPoNf0UQ0M16ocLWwu9cy3jqcvq67RAQDTrRpV9iwdK6566vel6lylCwmeQO8tzP3qBOn3iSOsHo6dvjatx4BGpJmN3930+yRYrh6A30HbFIzSNkRz1LE9i4q4Es1UW3M41Jyq2hjJo/eS66W+ivIXqucRolLHVr7Ntu5SBewsrIpRRQd4xmDl8ysSofed0H1MtdPXXY5MdlKG7MYSBEEQBKFH02tWdgRBEAShO2Gg9Wrs4Z4vdAyZ7AiCIAhCOhAzVsroNZMdy2HBcuwfGbqexr2XaidMav6GQ3OV0TyQig+Cx7MgKtrAszmo/djBdAx+lzLKb3+hgrahbU/fh0T3Z+EIsjhfy9aRtBWnGpMjJXAFjX+VaSrraGMF9erZ56N6UJSepuQ92mfcn0z258rObkvQt3z69u+QtMvWtlAkatJnXsD87tQElVbBaaf12G20jVl+pTfYXpdL2+Ck/Z/tV/6IwlwzwLQI1i7VLwfG7QHMDKZNaFHnGpH29TIkphWXCSXaThsxWth00zYx6Ue7uPeq9oZd9L69fYO8+FGh/vujVRtK6Xtgc9OO8GYoXU7ND2l7oyy2Xpb2ijUHmCgqTMede0/bX78R5jaKjwH9uyZnBf3CaDqT+r3K96vxneGkGqOESe+naZg6t+aHI0len/kftNleWwbVPa2/m57bd6XW4JOaSJ6LvWOhvqqfEmxcBfZSv0a6fs1gWqWMYfS7JhJT/e1w0msmMrV31Za64Fip3nremxHNjiAIgiAIPZpes7IjCIIgCN0KMWOlDJnsCIIgCEK6kAlLSug9kx1N9q5rQTJ202KNx1M9RCxbqyLO/JswGzyPt0LKsgHdGFX2fNe39pK8cl8DSW/7+vHJY7+H2uOzmA2+j1PF2TKZ6CUQ60vSif+Xr9r30Zo2Wn5ouO+LoKZVcFxM7ea7/FSM0Ffzw9NcSm3u/g3Ur4ezSeVXj6VaiWCcpbXj9bUFJC/cQsseX15N0nGt37w2qufw2Kk936E5hinPo/cailPBQWNI9YvdzhzKMP1MNF+7Lhs7RhN9beN+1SbTScekZx/TAmnVRnJpxdzPjj2sGmWnww4W87MTyVP3Y2P6Hu6qKNxP60N2b/2vqiLpdlz0HBbGWKohqZuoRkihj+qE+Lua0MZDgvk4MhM03Tyinf9crCyPTebRvga4Zsceos/R1L6LAmfQh9Mnh743uk6nPkS1NV4nHd8+n6qr8eu0DX0++RpJb71A6WcSHnrfbhbvas9oVbayeCfJW7eniKSbj1d1cf9Cnh1ti8OiA6lzogzm90rXycVjbcdnE3omvWeyIwiCIAjdCBEopw6Z7AiCIAhCOhDNTsqQ3ViCIAiCIPRoes3KjuUwYX3l+8YRVPbaMIu94mhhGgfdzM7dVTI7ejxbFTCzWAyaEPOFEVB+VLhPnrXNNJaT16Ns7o1hmtccofXmeNUj9blozBmuN4nfu0ddc/vJJG/wlVTDYyWUemLjE7SsLUSVFZl+ZTs3mR+PGHW7g/rjlY6l8FUaDMusriVpR1T1Q/zs00neJ5/3J+n80gatDe37GW1hsbH0Pm0K0/51O+i99stuTB5z3z45GVRDYGqDyealeeEY1SI0a75HuC+UcIymDS0mkBGnr3TET5LEhYjF5A9cHxErVpqHRD0t7GyifeqqV22KZ7Ofmy6mT9K1K+zRmJ2IvwRD8zdkoxVvnEb7ZWiR0mmFWZ9FE/xrUT2LligdK/5sqpfZF9YGOBMccV85CS/tF1ejaj8L0YeEi96P3v/hIuafinV/WNOOFWZRHzf7gtRvTePmnOSxf2ADyQveRvWBhYbyr9X4SjHJa8pjzsKy1cCras4hWYF6qiPK3agHhqP3zb+vQ0Wqk/Pz6djJZ/6z9Pdvz95skkfGZKL974uuRMxYqaPXTHYEQRAEoVshZqyUIWYsQRAEQRB6NLKyIwiCIAhpQMxYqaPXTHbsQTtsX8VE0u3h3H+IjZql4dTMwNwPSYSGQiKagTgrHGthJ3uVrdlimhKDaSeCIaUTcLup7obHamoIKU1PnGk9uBYhllAd0SeP2vI33juapI+7YZlKsPVAt4f5s9Dam5XFOngkvU7LOs127qZ2fjNKH8be65VOh/tHitfTewttU8b96GBajy+PxcKqp0Ki4jylw+GaKBv7dqnXNFT9s5lPoSCtN9et+qI+QrVXXGfhcann3BykbdA1OgBgaRoeG/N/Y6eyLeJjytFCy0b7cpGJOjSdLBZWDk3bg5pehoVc8+xkfoE0iQavV9eGHQrD3ravlNBEOn5HVWwl6Xy3GgOf1xeSPP5O2TVfSi4vff/CzJeS3d12XKVEI9X76H0GAAltSCRYjL5IH9qmcJHWb7V0LCXeojqcwLfUOxdkQcyam2n8LsuvHp7HQe/FZPqZ3TuVMyA3+y4sLqLvgqWdu2MXdSKU8QV777XHyv0Ncd9Euuarmekiq+IsFpl2P/w7l8Qt4/7TjiZixkoZvWayIwiCIAjdCpnspAzR7AiCIAiC0KPplSs7UZ+aDtvZsr/BVtHjWerYTncLwxmgaUNbGrVsbAs7W643oZZKTQfbSszMWJZLNSrI3Jzzbesxl8pvZtumLbYMnUjoW3XpNTPKqbnJdvIJyWP/R3SNPXoW7Ri3tlXez7ZYN7A2NOWo9sc3fknybBl0S2rRq5qbeba12NxdQ9LVV5+UPM7rS++lbybdotripfezc29O8jgrk20fZ9fVQ0vsaKH7vFtt9dfK8mfBMTRzmc3GTT30XJvmSsF08Z96LHSD9sazaCIwwjy0hHYu2z5uuWk6rlkj7XVsTzu7Vd1UrIekAACbiz4LbspsD1uWelmbrqEv5wkemq5q4TZohW7e3d9IdehmzzSLuXeod6oxy10emBn03ESUuWXIVGmTmbH4s9JDcpisuXVjaZ/ZwuqhRxqp2crlY+4pstR451vydXM6AOQVqPeq36BGksdNXuu29kse579H67ElmMuDTHVuVhXNC/VhferWpAO5tL0u5iYiqj1X/r1p6pdJ4WqJaHZSR6+c7AiCIAhC2hEzVsoQM5YgCIIgCD0aWdkRBEEQhDRgWBYMvhXzMM8XOkavmezYm2ywf7VFV996zu3dFtM86DZRy2CaF+YR3V2nnRdnZT3M1uzUj5meJ4NvLda0QI729RsklADTSlgszIAVV2lHNrXzG8wY3NJfaREaRtG9xfYw1WjoW9FDbKtroIHqcIbdtVW1R9NcAED49KEk3dxP1dVcxm5uJHX/nuHZlzzO9QZJnsOg/et10PspyFVahMYg1TiEIlRvkKHpk6J2pnFgW5g92nVcTPsRCNHrhCPqXuMROkiNduQ+zuZDbD3Xmsi6odXYgp5m293BkgipNsZz2PZr5oYh4VH1ZuyiY7I9jY7hcLaZBwCGX40Bri8pclPNjs+htCnbgvkkL8q+FHQXDnluOpa4ewGPWz3jGNv6HIvQfmil49OanLmLZMFgbiT0ZxcuZv0dpNf15CiXB+Z6uk3dfQrVs+lhbPj3h8NFNTBepz6eaV5NkL7LmZ+pmwsV0Hr9W9rW7HBNStkre0l6y/eUi4kI00CF7XS8xPaod2zIX+hzrD5d3XciIlvPeyJixhIEQRAEoUfTa1Z2BEEQBKE7IbuxUodMdgRBEAQhHYgZK2X0mslO1g7A/pXZWLcZR/3Mbu5hmhjNv4iN+SFxNDN377rsgtm7uY8e3ebOXfybcaZj0J8S98Fj46Ij7ZCVdbQwu7+mK4oluA8eWm31aercimdpZtWVbb9xjUGqEQALm2GFVMcYpcUkz7uJ2ucdIeUbpe9S6iundixtf1Gm0iKEWZiMlhjVcwQiVC+ja3wsL623rolqjqKa36OGfVSn4M6ighndtw53V899suhaLKuZag8cTfQ56uPHYPIN7jdK14pF/e37ztHDUnC/JDxESMStKuY6skQ+Gx+aJsJTh45jsfbxcBHaWBp4O23fm78/nqS/WbApeVyWQUMbBBO0v7c0K03P9uYckre9lup99CZyf1kci+mgwpoPmYSbPuMMpuEJ91XH/HuJu3AK7lbj0kEjSSCxhsZjMPupMcufeYJ9LzVo73arUBLV1I+R7WtKN8THR/B41mD9XWD11p5BfVkZibbDc8SjdHzkfarSRoy+GEV/WK7Os2L4rM1ahWOVXjPZEQRBEITuhJixUodMdgRBEAQhHYgZK2XIZEcQBEEQ0oCs7KSOXjPZKXhlMxy2/VqNrVcfl/zcYq47zAxqy9VjRnHXIjFmaia+dFhhZwvTZGhJPmB5vCDTpwtxWJ6X2qwNzdeEs4HarLluyK65NLFxnyBMzxHPUf2y6xtU85KI0huIaDZ5k/kW8X/OdBa6OKipheYxhzKbrtS0Hj6qs4kzHUueprtpiVKHSF4n9eWSYHqZ3Y3Kd0qc9YuDxduJaP5wLObnI7KP6ZV0DYyb1mNn+oh4e/6SuB5MkwbZQzQvVEjTCbf2rLhfHYY+9nn7XE427rRBzPssFmY6Lc0PTKgvOo7Bg3kxbVuj8qVjM9m9PXwcSb5SrPRh9ggt699MtVabL1bj3T+wgeQl6um74MhXL5l9PRXIRAfSF9DuYbGytGYk2JdLuC9NZ29ThWNZtB/cDex9zFH5gUFMp8X8iuUtVe9K/XCm2cmkY7YlpJ5rSx3VsvGYW/4MNTCbmE+p+JdU66b7YbL8tI/c2W37YYoG6bPwsrh2dadr30uOHJJXaBuWPDYSYWBVm5cRjlF6zWRHEARBELoVYsZKGTLZEQRBEIQ0Iaao1CAelAVBEARB6NH0mpUdKxiE9ZUTkv4vKr8azQNpbJsd/8ViY2k6Cz4B59qJhEN9YIsw/UYe8xGi+bhp5WeH2dFN3f8Ju6bNQ+3oum+deJzXy/yfaNd1M38nNqY/CZWqeqN57Jos5papaXYyNlFRVOMYassv+k+5SmzYRvKMvBySzs5TOhxdAwAAhRnU7040ocVqYhqdGIt9xGMY6TqcBNObRFifEu0N+4nmqKfnOptUOtif+oExuU8WLdaUjek37NG29V8m16DxN1x/VFyE1g5mgj5ju43dqxYTynJRnQUL14WoFucpVNaJIIheT5t5oTEDSbr+ePqM9fFuZzKQ5lIW9E5/81+lfmn83NVPgabTOYHGnTLrqIbLYLot/d3mPrL4ix/1qXRmNfPXk8PeR62NmdtpXrAfvU5Q03i59zG/XHu4fyd1HKHdgrIB1STt0OLENbbQfsjcSe8t73M1lqrHUh1OsJy1SX83Mum483npd00spl6GpjNpnzUMU7ohM+RInWbHsvb/deZ8oUP0msmOIAiCIHQnZDdW6hAzliAIgiD0AubOnYtTTjkF2dnZKCgowEUXXYT169eTMuFwGNOmTUN+fj6ysrIwadIk1NTUkDLbt2/HxIkTkZGRgYKCAvz85z9HPE5X1t555x2cfPLJcLvdOO644/Dkk08e7dtrF5nsCIIgCEI6sLrg7zBYsmQJpk2bhmXLlmHRokWIxWI499xz0dKi3H7MnDkTL774Ip577jksWbIEu3btwiWXXJLMTyQSmDhxIqLRKD744AM89dRTePLJJzF79uxkmS1btmDixIk466yzsHr1asyYMQNXX301Xn/99cPuoq7CsKyebfQLBALw+/04+bI7YHftt/F796gZ6K6vM5HDcdTXi82ubLvcX4gZYwb7Fl1nweNbMS2QbjLmflSYPwtobeDxdgy2jmnXNEbxCG0fj7Gk64Z4nC8XDRcEuyYxaR5A7d1cY6TH77L7qCCCx/0a/AdV8a5vZJO8Qf/1JUnXtKh8v4dqdnSNDkB1OU1hqsGwmNjKxvQnjfXKZ4gVoH3mDNDrmE7NJwgbDjYqy4F7r7ou14nEqKsRhAtUH/OYZs4ALav7T4rm0LwEk7UkvKpe00ufo83DA2tp13TTvGFF9JfeZ7uLksfFeY0kryVKdRd1depmrSbav6WLaBO8zy9HWyTOGU3SjmbV4UaYtnf9tXRstfczzxFg70Kj6oj+j28keZ//rpykZ499KXnsYQPgl4svpfX2oWM40qCNU/ad4NzL4qNppzqpNAj+zfT7I5yraaQKaMUtQ+lA9OWp779IjOnV9jK/UZo2L69/A8kalLOPpKPa+1gVyKFt+IjGF4sUqPYb2bQP++TRm9Vj72V46L1EorT9uv8n/t1o7FN9b4bD2PaLX6GxsRE+H9V0dhUH/i+dcvHtcDjb1p4dingsjBXP33LEbd2zZw8KCgqwZMkSfOMb30BjYyP69u2LhQsX4rvf/S4A4IsvvsCwYcOwdOlSnHbaaXj11Vdx/vnnY9euXSgs3C/ymj9/Pm666Sbs2bMHLpcLN910E15++WWsXbs2ea3LLrsMDQ0NeO211474fjuDrOwIgiAIwjFMIBAgf5FI5NAnAWhs3P/DJC9vv8J85cqViMViGD9+fLLM0KFDUV5ejqVLlwIAli5dihEjRiQnOgAwYcIEBAIBrFu3LllGr+NAmQN1pAOZ7AiCIAhCOugiM1ZZWRn8fn/yb+7cuYe8tGmamDFjBs444wyceOKJAIDq6mq4XC7k5OSQsoWFhaiurk6W0Sc6B/IP5LVXJhAIIBRibt5TRK/ZjRXLMpJbr02Huu3cL6gJo556lUdCM7uYLXQp2WLmD0MzP5kshALYNlNLX6duxxQFAJaWdDJTA1vtpnkeWm+UmW+ssFrKjXpZmAy+h1lbPXaw0Ay8EXGvFmIjQpdovaV0GTp6h9oy3s9B97+bPDSGlg7FafsCYXodu7bVlZut+FbzMHMzb4Xafi0MZmF0aqE9uDmBb8fVf1rEaSSBVmYsW0TVG8uhF43lttk8gLUPTjYONXOqzcWeeTtbO+x2OibrI9SkUZirbGtOG603yswh+hi2WBt2XMDG3bdPTR7619F6Go+n59o0cxl/Tnw7PHEZwEI+cHOkoZ26acZgkufeQsves3lS8rjgY3rNjCupewS+fT+RqS5s2FiohiDbOl+vBlPUT9uw82z2nmvfJ5lFtA0V2XTQ5riVTTTAwqzsZm0KaubewixaT9xq+3d030zahqwz6CqE/i7nZ1BZgd4+AKjRwpbsa6IvlddNTWBhbXzn/4uGt9CJx4BtbeZ2LV21G6uqqoqYsdxu7j6hNdOmTcPatWvx/vvvH3kDjiF6zWRHEARBELoVXeRnx+fzHZZmZ/r06XjppZfw7rvvorS0NPl5UVERotEoGhoayOpOTU0NioqKkmU+/PBDUt+B3Vp6Gb6Dq6amBj6fD14v036lCDFjCYIgCEIvwLIsTJ8+Hc8//zzeeustVFRUkPzRo0fD6XRi8eLFyc/Wr1+P7du3o7KyEgBQWVmJNWvWoLa2Nllm0aJF8Pl8GD58eLKMXseBMgfqSAeysiMIgiAIaSDVTgWnTZuGhQsX4l//+heys7OTGhu/3w+v1wu/34+pU6di1qxZyMvLg8/nw09+8hNUVlbitNNOAwCce+65GD58OC6//HLcc889qK6uxi233IJp06YlzWfXXXcd/vjHP+LGG2/EVVddhbfeegvPPvssXn755SO/2U7SayY7LUWA/StTcMULyr4czaVaj+i71O4bHKfsy648ai/m2xf10AKGi9q37Syth1QwmO3e7mBaCi2f609M02izLNdZcL2PqdnyE3VUtxDLZtuSta30dtoNrTQOLk1vEi5l22BD9DotbpXmdv6ojQtQFI0h+txCYVov6TPWRwke+sBJrxPXtB/2MD3XEQRNa1q7RDtRBgAgVKh9wJ65ybU1WpstNnYcGUy3pWkp+Db6snzqQ2BXgxJ4uJy0nqYmurycna1uzsHGUm0j3cqtb/uN2ujXSm4m7bQGbbswb28z74ZqVTZUQPN4SBZd52QyvZrBxoBrixo/NuYGwMZ24Ovb93M/Z+3tx97HdqQSJxTQZf0VGwfQNu1U4y6az8Y+c10ROUlpWQy+TZ31aZYWNsFlZ9pB9n3i0XxM1Jt0PBRkU62Nv8/e5DHX18VN+mz0d9vG/kNzjZeu/3GxvCDT6rnt6mHxreccU2tT7cRw2+WCYeDZdqvqOlIc9fyRRx4BAIwbN458vmDBAlxxxRUAgPvvvx82mw2TJk1CJBLBhAkT8PDDDyfL2u12vPTSS7j++utRWVmJzMxMTJkyBXPmzEmWqaiowMsvv4yZM2fiwQcfRGlpKR5//HFMmDDhiG6zK+g1kx1BEARB6M10xK2ex+PBvHnzMG/evDbL9O/fH6+88kq79YwbNw6rVqUqyNihkcmOIAiCIKQBiY2VOmSyIwiCIAjpQKKep4zeM9mxK23JntGaP4IA1SL0+ZQ6PNp0iua7g/mZaKW18TJjv4au0QEAS0tbTEMSZSEV+HXaQ9efJGLMZwn7GaC/J4af+qQwmR5J94Th3kfb56bucRDTdkC6dzOX80W07J64Kpztp33PNUfhiLqfGHMFnwgz4ZAeroPrmjLpvfKQEHp4Bq7f4JodXacTZTs/Laaz0HU6iSymy7K3Xdbtpe1tNe60cZnFXeYnaD/19SndBddZZLnb1jzwkA9cH6FrQfY0UKdBeT7qK8WjaYVCUdr38RZ6HfhVWaOR5XH/Tpo/IluYjtH8FbQfbHHVh/4vqX5jz8/pOCy8X+n44lnMR1M+q1cL5dFUSvN2fEh99DgizP+TTwsRkkvbVNaHaq903YuDfS9xTUx7hNh3xM4WpekKs7xs5uOmOaaeB78m1+y01z7efp1wov1/UQntOk5WT0uE6fi069rsXEMpG5N7Or1nsiMIgiAI3QgxY6UOmewIgiAIQjpI8W6s3oys3QmCIAiC0KPpNSs7sZwEbF/Ff6rPVLbyviup3bx+CPUtYd+r6Syyqc7C5LGxoprWI4fpH1rNwLXrtjC9iYfFxtLWKrnPGI6uBWrl24dpgSw9fYhfCJbWpiiL+eQI0no9e7S8CIv/U0WHXNyr0g0jeVAi1r+6riXGdE1Mo2HEVT+4GplmJ0Gvw5eC41rYHBeVSoDJXIhmh2t0EszXi6XpdHhcKpuD+WXSNAVepo/h2pqmsGoE10O0SmsPOmYyP1FMZ6FrgTJd9Jpc66HrIzK8TDcUp89c9+2i67D2X5SNfe05x7LYg2LJIU8oQZW9geqEzEzql2nrJTnJ4/yPqBCrZCZt/5dXqAHvrSVZaBrI2quP0Uwq+LK723/mut4qx0N1Q3VBGsspL0O1metjuC8dHf7cwnGeVs/KzjQwvKzuH6ftK+7H1L7vEuy58XGnl42yseNox09QmJWNsJhs0SAba3o92ne32U5svK5GzFipo9dMdgRBEAShW2Fa+/86c77QIdJqxpo7dy5OOeUUZGdno6CgABdddBHWr19PyoTDYUybNg35+fnIysrCpEmTWgUYEwRBEIRjDqsL/oQOkdbJzpIlSzBt2jQsW7YMixYtQiwWw7nnnouWFrUEPXPmTLz44ot47rnnsGTJEuzatQuXXHJJGlstCIIgCMKxRFrNWK+99hpJP/nkkygoKMDKlSvxjW98A42NjXjiiSewcOFCnH322QD2x/AYNmwYli1blgxM1hGMqAHDvt++a2i6ln0jaLmMnTSth2uyhZhtvJ6m/ZvUNDujhuohak6lPh9MzXwcy2Z+U+rpuQ4tJE0kn8VUcjFdiFO1KR5lmhfm68eIavG5mBZIj5vVKp8JV2I0TBL5tWGy2FIxKj2AQ3PgU/4yLdtQQYdnS5nmJ4P5v+ExrAyt+U4a0qdVmsf2smlt4jbxGPOlo2t24hmsMPezo+k52tPoAFSnw+NScTJcyg8P11nY2E8/XQ/By/K0rqXgeohwtONfHQmmDYpr58aDtB6DabFs2hh1BOkzHvTQBpK2EkrPEa+jYitH3z4knbNB+ZPZ/L182r5MpnPqp358RYaxeHPsvXFqfq6cLMbdofzfZGi6KJM5EfJ7qY+bBIk1Ra8TZZq02ib1gvI2ZTAtlq5sak8fA1AfTnzscHT9jJ1pHXl7Y3GVDuylsQqz8qm+So8NaLD+jTHtjaX5DtM1fQBgi2npEPtCOIoY6KRmp8ta0vPpVruxGhsbAQB5efsFgStXrkQsFsP48eOTZYYOHYry8nIsXbr0oHVEIhEEAgHyJwiCIAjdjgMelDvzJ3SIbjPZMU0TM2bMwBlnnIETTzwRAFBdXQ2Xy4WcnBxStrCwMBmanjN37lz4/f7kX1lZ2dFuuiAIgiAI3ZhuM9mZNm0a1q5di2eeeaZT9dx8881obGxM/lVVVXVRCwVBEASh6ziw9bwzf0LH6BZbz6dPn46XXnoJ7777LkpLS5OfFxUVIRqNoqGhgazu1NTUoKio6CA1AW63G263u9XnjqANtq80CLqJm2teorn0PFeDbhWlFtLy15pI2rZ+uyrpoF1bXtOXpJsGK81Acwm1EedupLGQIrkqv5HNT50ttE12TW+Sv5batx1b6S62urMGJI+DfWk9PAaQ2baLilbEdTM7MypbrJ6o9qj2nUCvyX3aZO7Q7PPsJY8y3ZBbk2y00t1QGQAiVLJBNRvctQuLYWU52y4LHtNMa0iCxR5zsFhZpqaXYVKrVnGqdL1EjPksiYL5FNLq4v5ZuOZB1+zw1fLwVipecpUrIRTvhgTz76THMeP+kWxMs+PZoxrc/090p2Z83z6SdhSod8yRTx/qvm/TuFT1Q1W9sX4RkudisdN07YrJ+szj5mXb1q5w/zcJpveJa9oV/iy4XkZvU0uUft9Vf0nvPX+Vqnff15n/rz70O6yuVr1IHh/tF18G1Q3p44fHt+JjK7hctcl2UiPJ4/cWalb3497FvjDW+0nS0IZ3NJf1mZ/2t3uP9p3M3ik95l0iksI1APGgnDLSurJjWRamT5+O559/Hm+99RYqKipI/ujRo+F0OrF48eLkZ+vXr8f27dtRWVmZ6uYKgiAIgnAMktaVnWnTpmHhwoX417/+hezs7KQOx+/3w+v1wu/3Y+rUqZg1axby8vLg8/nwk5/8BJWVlYe1E0sQBEEQuhuGZcHohMi4M+f2NtI62XnkkUcAAOPGjSOfL1iwAFdccQUA4P7774fNZsOkSZMQiUQwYcIEPPzwwyluqSAIgiB0MeZXf505X+gQaZ3sWB2YlXo8HsybNw/z5s3r1LX8Gy3Yv9LnRP3KYFvyzCZ2Qab3adFi1LipVgIxaq+PN1JbNIHpC3yNSpvki1PbcnzXbpL2lvZLHue8xwUotA1WVEuz/uXty31Z3Vuencfnov1gBZRt38iiope9E6j5MZyj+pe73/Dso20K52sxc1jXc/83wZK29TG6tgOgvn9a1cuMt6aTVmZq/nFsLOiP5eY3pApYTJuix0oDQPx3GCxuEtctRDVfI9YhfNroPnpsTCfE/Z9Qfyc0j+tN2tPsuPtTrUcsptprJlgMNi460u7VzmOasf7u/5jypaP70TkoCa39dlpvPIP5YdKHEntusTDrb338sFvRfcLsb4KmY2EaHa6fshu0v3V/M1wDw2mOqEYFw/R7ycZ8yDQSuRJ9kHVVOSSduVm1oWV4u02ATWtjnPVDPErTFe8qvU/VCNZn7Fz3FnVvORtpe3Nf/YKeq/lTsrNdu4HxQ0m6Ufua8rCYd7a4pqeLympJT6RbCJQFQRAEobchZqzUIZMdQRAEQUgHshsrZfSayU7+O9vhsO1fHrVatH2GzDQV374jJe05nOvEd+w8dKEjINHQ0OGyZInYRbeD9nl9Cy2smdJqJh1Pr+mhS+x5Xyhzzt6RzL07a4O+6p9ZxZbqh1GTgbNBLY3Hs1gIhSgzaXCLgWbjaLXl3kOvY+hb0bmVxcttYNp5bGtxnG1FN3UzwKH2TGpV2ZzMxT8PA6LV68igprS44xBmIr1eZnaLh7WOYlt3Lb4FXzuXunYABvyJhoDQOZzxCub6IZpFs/UIFrYW1vesKt1QzMN8WDY2lrR7bfWjm/WvadDr6uEkeAgFXpfe/zwsgqM4RNKxPVoQCDYevDuY+UnvJ2ayCzV4SNruUePHbGDuEILMRYZmbrKtzaNtYCalkmeVtCBeU0vy2huhfHxk/mMZTWvHNheTJGhm/LjFtucfTTrrBVlWdjpMt3EqKAiCIAiCcDToNSs7giAIgtCd6KwXZPGg3HFksiMIgiAI6UDMWCmj90x2PB7Avl+zY2g6nfjmrWlq0LEFsYcfQjth83qTx3mfU/1Aw2Bq928qU0PQxXbuhwpo2tC2MLf0o3n2EAuTUKhvwef10GHPQ4YQrQLX3bSKf6HONVxta3RancqqMVmYBGhpg2tgsti2dW3rtMm2u/PwFro+he8IT7CvA5uHXoeUZTqRrPVKsxNjoTu47smjeWEo+xcNYcK1NsjMUMd79rbZHgCArp+J07ZzFwJ6WBXLyXRNXB+jbaVvtSGc9a+dX0gjznQ47WFjP9mbm+h7k7VKSw9gLgNiLLSElnZWszAl9PWES4uv7Gyiuhbi+gGAo1o9WIMNlSwmSUysUVvGK1oGkDxzNx0D8RBr1FHAjLatyzGtWJt5wrFL75nsCIIgCEI3wjAPskHiMM8XOoZMdgRBEAQhHYgZK2XIbixBEARBEHo0vWZlJ755K2BwpynC0cDUbO7Ge6tIXv5HXpKOna580u87gcZ18G2hv1qa+mu+RbJpXiKDrefq2cwnTMJHBQYG08tYTlWXYXE/Km1rMnhYBKOVfxntkPk74SIeb5V6NXm4CwSolsKtaWDCfWnRhJe2QQ/HYDYzXy4OWtaMa++Li/avq5a+S7p/ltK3IyTP/VkVbZSuyzHYs9lDw6pYu6vRYfR6mS+oVn2oXZaHJbEzvztkbLFvTNPidgQtzAf3ccTCaNhYeA49zUNWFL7KY56osm7mqyj3z0tpi0YOSx4bNbR/uR+bwyF08djkcfYKKtJpzzeY6CQ1xKlgyug1kx1BEARB6E5IuIjUIWYsQRAEQRB6NLKyIwiCIAjpQATKKUMmO0JKMZkPDfvilcnj4pU5JK/ugmEknfuF0ilkVYVJ3q6vZ5B0S7mW4K5xuF8dpkchxXksJKbhgebjhvtcaaXL0bN5FvONEumj6YZYnj3E9D17VcXOZlaWuRNxtqiyvg3MsVGCxf0KaSeHqQ6nFQ5N58Id1XDfOfoXtMn6N37kPk6soKYVY/Vm7aBtCvVR/WRj/Rv3Mq2NWzs3zuJoMR89lnbdhMnEQEzDZTEJof0zFb2pfBntB+frVIdzOCQ+/fyIz20P7/PLk8dte2QS2sXCQZw3Heb5QoeQyY4gCIIgpAHR7KQO0ewIgiAIgtCjkZUdQRAEQUgHFjqp2emylvR4ZLIjdBsSLOZW/rvUd0esLD953DSA+usp/z/qj2Xj1UXJ47ivbd84B0XXVvAvk0Q7i6E82BQPSazpfbhPHtceFpdKE0H4NjOfQjRMEjJqVWHvp8ynjY21V48J5GUVtVA9VXwf9cmi4yigDn2sQFPy2MjM4MVpWU3/Y+T42i17WOg+e5hmJ+/trSSdKOmTPA4X07EUKKfPIqhpbxIe5ovIxXwraU57uLzL8tA2mU1UtDPgHyr2V2LdBgi9ABEopwwxYwmCIAhCL+Ddd9/FBRdcgJKSEhiGgRdeeIHkW5aF2bNno7i4GF6vF+PHj8fGjRtJmbq6OkyePBk+nw85OTmYOnUqmpubSZlPP/0UX//61+HxeFBWVoZ77rnnaN/aIZHJjiAIgiCkA7ML/g6DlpYWfO1rX8O8efMOmn/PPffgoYcewvz587F8+XJkZmZiwoQJCIfV7tfJkydj3bp1WLRoEV566SW8++67uPbaa5P5gUAA5557Lvr374+VK1fid7/7HW677TY8+uijh9fYLkbMWIIgCIKQBlK9G+u8887Deeedd9A8y7LwwAMP4JZbbsGFF14IAPjLX/6CwsJCvPDCC7jsssvw+eef47XXXsOKFSswZswYAMAf/vAHfOc738Hvf/97lJSU4Omnn0Y0GsWf//xnuFwunHDCCVi9ejXuu+8+MilKNTLZEboNjvx8+kGQakicXypdTu5m5pfm+BKStmtuYeLM/w3a838DwNK0Hwb3ncN1ORpcotPKmY52rnsf9cHiZvKY0ue2qmtqehgAMLIySdoKaT6HuA4nzjyg6D5wmP+b9jQ6nHjtnrYz2ZK2Y+hgkk58oS2LNzJfP51A13w5+vahmTHqt8a+ZVfyOHMn1c649xWRdP1QpekJFjKfPBnMJ4/WpQ7qCgqw6EJ62ZwPSPowlWWCkCQQCJC02+2G280DwrXPli1bUF1djfHjxyc/8/v9GDt2LJYuXYrLLrsMS5cuRU5OTnKiAwDjx4+HzWbD8uXLcfHFF2Pp0qX4xje+AZdLxfCbMGEC7r77btTX1yM3N/cI77JziBlLEARBENLBAYFyZ/4AlJWVwe/3J//mzp172E2prt7/Y7KwsJB8XlhYmMyrrq5GQUEByXc4HMjLyyNlDlaHfo10ICs7giAIgpAOumg3VlVVFXw+tbPxcFd1egOysiMIgiAIxzA+n4/8Hclkp6hov/m2pqaGfF5TU5PMKyoqQm1tLcmPx+Ooq6sjZQ5Wh36NdCArO0Jasedp9lsPe0GjLLCTFruJ6FQAuPZQnYi3VtVlMm2K6Ww/RpGl+dmxmGanlWLHocdNYnGpwvS3hH+DyvfupQqN7HeYX5Ucf/Iw0UQ1O+BpnQbWvEK65GxZWswtJ+0XrnOJ79mLriD+xcZDF+oCHP003Rbzs9NKu6TDyjq20i/qvtVqLAVOol/WjiB9jpFc1afBvnQ8FD5INTqC0J387FRUVKCoqAiLFy/GqFGjAOzXAi1fvhzXX389AKCyshINDQ1YuXIlRo8eDQB46623YJomxo4dmyzzq1/9CrFYDE7n/i/XRYsW4fjjj0+bXgeQlR1BEARBSA8p3nre3NyM1atXY/Xq1QD2i5JXr16N7du3wzAMzJgxA7fffjv+/e9/Y82aNfjRj36EkpISXHTRRQCAYcOG4dvf/jauueYafPjhh/jPf/6D6dOn47LLLkNJyf4fG//93/8Nl8uFqVOnYt26dfj73/+OBx98ELNmzepER3UeWdkRBEEQhDSQ6q3nH330Ec4666xk+sAEZMqUKXjyySdx4403oqWlBddeey0aGhpw5pln4rXXXoPHo3Z5Pv3005g+fTrOOecc2Gw2TJo0CQ899FAy3+/344033sC0adMwevRo9OnTB7Nnz07rtnMAMCyrZ/ubDgQC8Pv9GIcL4TCchz5BSCmOwYOSx9YeuvXZcLtI2ooos5bhos8yUVFM0kEtBEBLMd3mHaKWHSToZYiZqz0TFwAY2nZyZxM1W1Tctw5tkejCLdddhWPgAJKOb96alnYcKY7yUpWIs43c7X3NMTOWPs4AwLBrISAsWtbw0lATuunVitLt7t3xmQutiVsxvIN/obGxkYh+u5ID/5fGD5kFh/3IxcTxRARvbrjvqLa1pyArO4IgCIKQDrqRZqenI5MdQRAEQUgHpnUwb6SHd77QIUSgLAiCIAhCj0ZWdoSUQrYHA0BLMHloZGbQPObin+h02HZy+06q98neoX7xuAdRb56xHVSIEyijdYXzlfYmTiMzEI0OABz30CZVtob6nzjW3P8faxodTnz7juQxcWkAwGjH74gVjtAP7G3/BjQ8LBwH05XpY1Y0OsIhETNWypDJjiAIgiCkhU5OdnhgP6FNxIwlCIIgCEKPRlZ2BEEQBCEdiBkrZchkR0gplqbRAQAjSxPFuNr3g2SFQuo8nsk0PHpoCdcmGmnXxbRAmSvZdfXQAk6WZ9Arc52O0D1I1NWTdKuwGbpOh4eW4JodPYxJNhNxhalPHrO55bDaKfRyTAudMkXJbqwOI2YsQRAEQRB6NLKyIwiCIAjpwDL3/3XmfKFDyGRHEARBENKBaHZShkx2hJSSaGigH2hpPU7W/g9oTCsSh4hraXSdDXB4XwLsOlY4nDxO1O7peD1Ct4Vrq2xu5S/HcNKvQT0WFgDAq8qaXjruzC82dlELhV6JaHZShmh2BEEQBEHo0cjKjiAIgiCkAzFjpQyZ7AiCIAhCOrDQyclOl7WkxyOTHaHbEN/45RGf6yjtRz/QfadwDQb3ycMgvn+YvxahZ2BGlC4LPDQW8/dkZSrNjq2R+omSvTCCcGwgkx1BEARBSAdixkoZMtkRBEEQhHRgmujU+iD3/i20iezGEgRBEAShRyMrO0KPIL5jZ7qbIPQQeFwtXbclv6OFLkXMWClDJjuCIAiCkA5kspMyxIwlCIIgCEKPRlZ2BEEQBCEdSLiIlCGTHUEQBEFIA5ZlwupE5PLOnNvbkMmOIAiCIKQDy+rc6oxodjqMaHYEQRAEQejRyMqOIAiCIKQDq5OaHVnZ6TAy2REEQRCEdGCagNEJ3Y1odjqMmLEEQRAEQejRyMqOIAiCIKQDMWOlDJnsCIIgCEIasEwTVifMWLL1vOOIGUsQBEEQhB6NrOwIgiAIQjoQM1bKkMmOIAiCIKQD0wIMmeykAjFjCYIgCILQo5GVHUEQBEFIB5YFoDN+dmRlp6McEys78+bNw4ABA+DxeDB27Fh8+OGH6W6SIAiCIHQKy7Q6/Sd0jG4/2fn73/+OWbNm4de//jU+/vhjfO1rX8OECRNQW1ub7qYJgiAIwpFjmZ3/EzpEt5/s3Hfffbjmmmtw5ZVXYvjw4Zg/fz4yMjLw5z//Od1NEwRBEIRjjt5oLenWk51oNIqVK1di/Pjxyc9sNhvGjx+PpUuXHvScSCSCQCBA/gRBEAShu5EOM1ZvtZZ068nO3r17kUgkUFhYSD4vLCxEdXX1Qc+ZO3cu/H5/8q+srCwVTRUEQRCEwyMNZqzeai3pcbuxbr75ZsyaNSuZbmxsRHl5OeKIdcp3kyAIgtDziSMGALBSsNOps/+XDrSVWzDcbjfcbner8gesJTfffHPys0NZS3oK3Xqy06dPH9jtdtTU1JDPa2pqUFRUdNBz+EM+MAjexytHr6GCIAhCj6KpqQl+v/+o1O1yuVBUVIT3qzv/fykrK6uVBePXv/41brvttlZl27OWfPHFF51uS3emW092XC4XRo8ejcWLF+Oiiy4CAJimicWLF2P69OkdqqOkpARVVVWwLAvl5eWoqqqCz+c7iq0+dgkEAigrK5M+OgTST4dG+qhjSD8dmlT3kWVZaGpqQklJyVG7hsfjwZYtWxCNRjtdl2VZMAyDfHawVZ3eTree7ADArFmzMGXKFIwZMwannnoqHnjgAbS0tODKK6/s0Pk2mw2lpaXJFR6fzydfKodA+qhjSD8dGumjjiH9dGhS2UdHa0VHx+PxwOPxHPXr6ByJtaSn0K0FygDw/e9/H7///e8xe/ZsjBo1CqtXr8Zrr73WahlOEARBEIS20a0lBzhgLamsrExjy44+3X5lBwCmT5/eYbOVIAiCIAgHp7PWkmOVY2Ky0xW43W78+te/FltmO0gfdQzpp0MjfdQxpJ8OjfRR1/L9738fe/bswezZs1FdXY1Ro0b1CmuJYaVif50gCIIgCEKa6PaaHUEQBEEQhM4gkx1BEARBEHo0MtkRBEEQBKFHI5MdQRAEQRB6NL1istMbw9m3xdy5c3HKKacgOzsbBQUFuOiii7B+/XpSJhwOY9q0acjPz0dWVhYmTZrUyglVb+Ouu+6CYRiYMWNG8jPpJ2Dnzp344Q9/iPz8fHi9XowYMQIfffRRMt+yLMyePRvFxcXwer0YP348Nm7cmMYWp55EIoFbb70VFRUV8Hq9GDRoEH7729+S2Eu9sZ/effddXHDBBSgpKYFhGHjhhRdIfkf6pK6uDpMnT4bP50NOTg6mTp2K5ubmFN6FcMxg9XCeeeYZy+VyWX/+85+tdevWWddcc42Vk5Nj1dTUpLtpaWHChAnWggULrLVr11qrV6+2vvOd71jl5eVWc3Nzssx1111nlZWVWYsXL7Y++ugj67TTTrNOP/30NLY6vXz44YfWgAEDrJEjR1o33HBD8vPe3k91dXVW//79rSuuuMJavny5tXnzZuv111+3Nm3alCxz1113WX6/33rhhResTz75xPqv//ovq6KiwgqFQmlseWq54447rPz8fOull16ytmzZYj333HNWVlaW9eCDDybL9MZ+euWVV6xf/epX1j//+U8LgPX888+T/I70ybe//W3ra1/7mrVs2TLrvffes4477jjrBz/4QYrvRDgW6PGTnVNPPdWaNm1aMp1IJKySkhJr7ty5aWxV96G2ttYCYC1ZssSyLMtqaGiwnE6n9dxzzyXLfP755xYAa+nSpelqZtpoamqyBg8ebC1atMj65je/mZzsSD9Z1k033WSdeeaZbeabpmkVFRVZv/vd75KfNTQ0WG632/rb3/6WiiZ2CyZOnGhdddVV5LNLLrnEmjx5smVZ0k+WZbWa7HSkTz777DMLgLVixYpkmVdffdUyDMPauXNnytouHBv0aDPWgXD248ePT37WW8LZd5TGxkYAQF5eHgBg5cqViMVipM+GDh2K8vLyXtln06ZNw8SJE0l/ANJPAPDvf/8bY8aMwaWXXoqCggKcdNJJeOyxx5L5W7ZsQXV1Nekjv9+PsWPH9po+AoDTTz8dixcvxoYNGwAAn3zyCd5//32cd955AKSfDkZH+mTp0qXIycnBmDFjkmXGjx8Pm82G5cuXp7zNQvemR3tQ7s3h7DuCaZqYMWMGzjjjDJx44okAgOrqarhcLuTk5JCyhYWFqK6uTkMr08czzzyDjz/+GCtWrGiVJ/0EbN68GY888ghmzZqFX/7yl1ixYgV++tOfwuVyYcqUKcl+ONj711v6CAB+8YtfIBAIYOjQobDb7UgkErjjjjswefJkAJB+Oggd6ZPq6moUFBSQfIfDgby8vF7bb0Lb9OjJjtA+06ZNw9q1a/H++++nuyndjqqqKtxwww1YtGhRyiMTHyuYpokxY8bgzjvvBACcdNJJWLt2LebPn48pU6akuXXdh2effRZPP/00Fi5ciBNOOAGrV6/GjBkzUFJSIv0kCCmiR5uxenM4+0Mxffp0vPTSS3j77bdRWlqa/LyoqAjRaBQNDQ2kfG/rs5UrV6K2thYnn3wyHA4HHA4HlixZgoceeggOhwOFhYW9vp+Ki4sxfPhw8tmwYcOwfft2AEj2Q29//37+85/jF7/4BS677DKMGDECl19+OWbOnIm5c+cCkH46GB3pk6KiItTW1pL8eDyOurq6XttvQtv06MlObw5n3xaWZWH69Ol4/vnn8dZbb6GiooLkjx49Gk6nk/TZ+vXrsX379l7VZ+eccw7WrFmD1atXJ//GjBmDyZMnJ497ez+dccYZrdwWbNiwAf379wcAVFRUoKioiPRRIBDA8uXLe00fAUAwGITNRr9q7XY7TNMEIP10MDrSJ5WVlWhoaMDKlSuTZd566y2YpomxY8emvM1CNyfdCumjzTPPPGO53W7rySeftD777DPr2muvtXJycqzq6up0Ny0tXH/99Zbf77feeecda/fu3cm/YDCYLHPddddZ5eXl1ltvvWV99NFHVmVlpVVZWZnGVncP9N1YliX99OGHH1oOh8O64447rI0bN1pPP/20lZGRYf31r39NlrnrrrusnJwc61//+pf16aefWhdeeGGP31LNmTJlitWvX7/k1vN//vOfVp8+fawbb7wxWaY39lNTU5O1atUqa9WqVRYA67777rNWrVplbdu2zbKsjvXJt7/9beukk06yli9fbr3//vvW4MGDZeu5cFB6/GTHsizrD3/4g1VeXm65XC7r1FNPtZYtW5buJqUNAAf9W7BgQbJMKBSyfvzjH1u5ublWRkaGdfHFF1u7d+9OX6O7CXyyI/1kWS+++KJ14oknWm632xo6dKj16KOPknzTNK1bb73VKiwstNxut3XOOedY69evT1Nr00MgELBuuOEGq7y83PJ4PNbAgQOtX/3qV1YkEkmW6Y399Pbbbx/0u2jKlCmWZXWsT/bt22f94Ac/sLKysiyfz2ddeeWVVlNTUxruRujuGJalufEUBEEQBEHoYfRozY4gCIIgCIJMdgRBEARB6NHIZEcQBEEQhB6NTHYEQRAEQejRyGRHEARBEIQejUx2BEEQBEHo0chkRxAEQRCEHo1MdgRBEARB6NHIZEcQBEEQhB6NTHYEQRAEQejRyGRHEARBEIQejUx2BKEHM27cOPz0pz/FjTfeiLy8PBQVFeG2224DALzzzjtwuVx47733kuXvueceFBQUoKamJk0tFgRB6HpksiMIPZynnnoKmZmZWL58Oe655x7MmTMHixYtwrhx4zBjxgxcfvnlaGxsxKpVq3Drrbfi8ccfR2FhYbqbLQiC0GVI1HNB6MGMGzcOiUSCrN6ceuqpOPvss3HXXXchGo1i7NixGDJkCNauXYszzjgDjz76aBpbLAiC0PU40t0AQRCOLiNHjiTp4uJi1NbWAgBcLheefvppjBw5Ev3798f999+fjiYKgiAcVcSMJQg9HKfTSdKGYcA0zWT6gw8+AADU1dWhrq4upW0TBEFIBTLZEYRezJdffomZM2fisccew9ixYzFlyhQyERIEQegJyGRHEHopiUQCP/zhDzFhwgRceeWVWLBgAT799FPce++96W6aIAhClyKTHUHopdxxxx3Ytm0b/vSnPwHYr+V59NFHccstt+CTTz5Jc+sEQRC6DtmNJQiCIAhCj0ZWdgRBEARB6NHIZEcQBEEQhB6NTHYEQRAEQejRyGRHEARBEIQejUx2BEEQBEHo0chkRxAEQRCEHo1MdgRBEARB6NHIZEcQBEEQhB6NTHYEQRAEQejRyGRHEARBEIQejUx2BEEQBEHo0fx/PW+ymU2ODXwAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "topo.depth.plot()" ] }, { "cell_type": "markdown", "metadata": { "tags": [] }, "source": [ "## 4. Modify Depth\n", "\n", "**Important:** After coarsening the grid, use the Topo Editor widget to erase disconnected basins.\n", "Carefully inspect the domain and make necessary adjustments to ensure physical consistency. For example, you may need to connect North and South America to prevent unintended ocean passages or isolated water bodies." ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "data": { "application/vnd.jupyter.widget-view+json": { "model_id": "d6bd814bae1b4101b1e52f0ae4c6d24e", "version_major": 2, "version_minor": 0 }, "text/plain": [ "TopoEditor(children=(VBox(children=(HTML(value='

Topo Editor

'), HTML(value='

Display

'), To…" ] }, "execution_count": 9, "metadata": {}, "output_type": "execute_result" } ], "source": [ "%matplotlib ipympl\n", "from mom6_bathy.topo_editor import TopoEditor\n", "TopoEditor(topo)" ] }, { "cell_type": "markdown", "metadata": { "tags": [] }, "source": [ "## 5. Save the grid and bathymetry files" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [], "source": [ "# MOM6 supergrid file:\n", "grid.write_supergrid(f\"./ocean_hgrid_{grid_name}_{datestamp}.nc\")\n", "\n", "# MOM6 topography file:\n", "topo.write_topo(f\"./ocean_topog_{grid_name}_{datestamp}.nc\")\n", "\n", "# CICE grid file:\n", "topo.write_cice_grid(f\"./cice_grid_{grid_name}_{datestamp}.nc\")\n", "\n", "# SCRIP grid file (for runoff remapping, if needed):\n", "topo.write_scrip_grid(f\"./scrip_grid_{grid_name}_{datestamp}.nc\")\n", "\n", "# ESMF mesh file:\n", "topo.write_esmf_mesh(f\"./ESMF_mesh_{grid_name}_{datestamp}.nc\")" ] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [], "source": [ "# Also coarsen chlorophyl and salt restoring input files:\n", "\n", "chl_orig = xr.open_dataset('/glade/p/cesmdata/inputdata/ocn/mom/tx2_3v2/seawifs-clim-1997-2010-tx2_3v2.230416.nc')\n", "chl = chl_orig.isel(LON=slice(0, -1, CX), LAT=slice(0, -1, CX))\n", "chl.to_netcdf(f\"./seawifs-clim-1997-2010-{grid_name}_{datestamp}.nc\")\n", "\n", "salt_orig = xr.open_dataset('/glade/p/cesmdata/inputdata/ocn/mom/tx2_3v2/state_restore_tx2_3_20230416.nc')\n", "salt = salt_orig.isel(LON=slice(0, -1, CX), LAT=slice(0, -1, CX)) \n", "salt.to_netcdf(f\"./state_restore_{grid_name}_{datestamp}.nc\")" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.11.14" } }, "nbformat": 4, "nbformat_minor": 4 }