{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": "# Spherical grid with custom bathymetry and ridges\n\nIn this notebook, we create a spherical grid with uniform resolution. We then implement a user-defined, custom bathymetry including ridges."
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## 1. Import Modules"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "%%capture\n",
    "import numpy as np\n",
    "from mom6_forge.grid import Grid\n",
    "from mom6_forge.topo import Topo"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## 2. Create a horizontal MOM6 grid\n",
    "\n",
    "Spherical grid. x coordinates interval= [0, 360] degrees. y coordinates interval = [-80,+80] degrees"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Instantiate a MOM6 grid instance\n",
    "grid = Grid(\n",
    "    nx         = 180,         # Number of grid points in x direction\n",
    "    ny         = 80,          # Number of grid points in y direction\n",
    "    lenx       = 360.0,       # grid length in x direction, e.g., 360.0 (degrees)\n",
    "    leny       = 160,         # grid length in y direction\n",
    "    cyclic_x   = True,      # reentrant, spherical domain\n",
    "    ystart     = -80          # start/end 10 degrees above/below poles to avoid singularity \n",
    ")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## 3. Configure the bathymetry"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Instantiate a Topo object associated with the horizontal grid object (grid).\n",
    "topo = Topo(grid, min_depth=10.0)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### *Define a custom bathymetry*"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Define a sinusoidal depth array:\n",
    "i = grid.tlat.nx.data                # array of x-indices\n",
    "j = grid.tlat.ny.data[:,np.newaxis]  # array of y-indices \n",
    "custom_depth = 400.0 + 80.0 * np.sin(i*np.pi/6.) * np.cos(j*np.pi/6.)\n",
    "topo.depth = custom_depth"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.collections.QuadMesh at 0x15426d0240d0>"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjIAAAG2CAYAAABlBWwKAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8ekN5oAAAACXBIWXMAAA9hAAAPYQGoP6dpAABWsElEQVR4nO3df3wU1b0//tfm5waS3ZhINkQIUAsGFMSChi20UoxEioiSa9UHYqpcfZRPUCFWLb2AGn8E6a2gNoBaDPYqpeU+hFYUEFCxSggYpIJ4EZWvicImLZosCexmszvfP8Iss7szu5PdTWZn9/V8POZx3dkz8z7vvXPI6ZwzZwyCIAggIiIi0qEkrStAREREFC52ZIiIiEi32JEhIiIi3WJHhoiIiHSLHRkiIiLSLXZkiIiISLfYkSEiIiLdYkeGiIiIdIsdGSIiItItdmSIiIhIt9iRISIiSgCPPvooDAaDz1ZUVOT93uFwoKKiArm5ucjMzERZWRmam5t9ztHY2Ijp06ejX79+yMvLw4MPPoiurq6+TsVHiqbRiYiIqM9ceuml2Llzp/dzSsr5bsDChQvx5ptvYuPGjTCbzZg/fz5mzZqFDz/8EADgdrsxffp05OfnY8+ePTh58iTuuOMOpKam4qmnnurzXEQGvjSSiIgo/j366KPYvHkzDh48GPBdW1sbBgwYgPXr1+M//uM/AAD/93//h5EjR6Kurg4TJkzA1q1bcf311+PEiROwWCwAgDVr1uDhhx/Gv/71L6SlpfVlOl5xf0fG4/HgxIkTyMrKgsFg0Lo6REQUwwRBwOnTp1FQUICkpN6bfeFwONDZ2RnxeQRBCPjblp6ejvT0dNnyx44dQ0FBAYxGI6xWK6qrq1FYWIiGhga4XC6UlJR4yxYVFaGwsNDbkamrq8Po0aO9nRgAKC0txbx58/Dpp5/iiiuuiDifcMR9R+bEiRMYPHiw1tUgIiIdaWpqwqBBg3rl3A6HA8OGZMLW4o74XJmZmWhvb/fZ98gjj+DRRx8NKFtcXIx169bhkksuwcmTJ/HYY4/hJz/5CQ4fPgybzYa0tDRkZ2f7HGOxWGCz2QAANpvNpxMjfi9+p5W478hkZWUBAB7eNRlNhgIAQIsjy/u9o6v7JzCmdE9WyjOeBgD8sH8LAKAo/fz/cwpTvgMA9EvqvvjOeJIBAI1dOQCA/3Pme8t+0ZHnE0spjlwspThysZTiqMlJjKMmJzGOmpzEOImakxhHTU5iHDU5iXHU5NSTazwec9Jbu43HnPTabs+0u/Gfkz7z/u3oDZ2dnbC1uPF1w1CYssK/62M/7cGQcf8fmpqaYDKZvPuV7sZMmzbN+99jxoxBcXExhgwZgr/+9a/IyMgIux5ai/uOjHjLLT0zBWmG7vG7lOTz/09OOdcAUlK6L+a0DCcAwNg/FQDQz3i+4WamdF9w/ZK6pxUlec597uouY0xN9Zb1j6UURy6WUhy5WEpx1OQkxlGTkxhHTU5inETNSYyjJicxjpqcxDhqcurJNR6POemt3cZjTnputwD6ZCpCZpYBmVnhx/Gg+1iTyeTTkVErOzsbI0aMwBdffIFrr70WnZ2daG1t9bkr09zcjPz87k5kfn4+9u3b53MO8akmsYwW+Pg1ERGRBtyCJ+ItEu3t7fjyyy8xcOBAjBs3Dqmpqdi1a5f3+6NHj6KxsRFWqxUAYLVacejQIbS0nL/TtWPHDphMJowaNSqiukQi7u/IEBERxSIPBHgQ/oPDPT3217/+NWbMmIEhQ4bgxIkTeOSRR5CcnIzbbrsNZrMZc+fORWVlJXJycmAymXDvvffCarViwoQJAICpU6di1KhRmDNnDpYvXw6bzYbFixejoqJCcTirLyRUR6bT052udOzW0ZUqW0ZkNJyfVS6OqXrvBJ77LC2jFCtUHOl5lOLIxQoVJ1hO0vkJoXKSnoM5Bc8pnOtBTU4+d6BD5NSTazwec9Jbu43HnPTabt2GyCfgxqpvvvkGt912G06dOoUBAwZg0qRJ2Lt3LwYMGAAAWLFiBZKSklBWVgan04nS0lKsWrXKe3xycjK2bNmCefPmwWq1on///igvL0dVVZVWKQFIsI4MERFRrPDAg0gGh3p69IYNG4J+bzQaUVNTg5qaGsUyQ4YMwVtvvdWjuL2NHRkiIiINuAUB7gjWpI3k2HjCyb5ERESkWwl1RyYtqXuNAek6FyJxn1hG5BDOP+rnXUPBb/0BaZlQsZTiSM+jFEcuViQ5Sc8bKifpOZhT8JzCuR5U5SSz5ko0frt4zCme22085hRL7dYp9N0cmb6e7BuvEqojQ0REFCs8EOBmRyZiHFoiIiIi3eIdGSIiIg1waCk62JEhIiLSAJ9aio6E6shkp54N2CcunCRODvMvc6orM+AYceEkcXKYXBn/84SKI3ce/zhyZZTiBIsVTk7h/HZy50mEnMK5HuTKKMUJFos5yceRxorFdqsmlt5y0mu7PesOnChMsS2hOjJERESxwnNui+R4YkeGiIhIE+4In1qK5Nh4wo4MERGRBtxC9xbJ8cTHr4mIiEjHDIIQ39Oe7XY7zGYzltaXwJgZ+DZWIiIikaPdharinWhra4PJZOqVGOLfpYNH8pCVFf79hNOnPRg7qqVX66oHHFoiIiLSgAcGuGGI6Hji0BIRERHpGO/IEBERacAjdG+RHE8J1pGxpNgBALkp7d59SosiNXd1jze2dJ4fd2x1ZQBQXnwpL82uGCvYIk/+sZTiyMVSiqMmJzGOmpzEOGpyEuMkak7SxbhC5SS3wFio60FNTj25xuMxJ72123jMSa/t1tnpQl9xRzi0FMmx8YRDS0RERKRbCXVHhoiIKFbwjkx0sCNDRESkAY9ggEeI4KmlCI6NJxxaIiIiIt1KqDsy4sSwoSmnvPv6JbkBAGc8yT5lxUli0klotrPd+xxd3T+bMcX3LanSCXb+sZTiyMUKFUcaSymOmpzEOGpykk5KDJWTdIIdcwqeUzjXg5qcenKNx2NOemu38ZiTXtttpyvwd+ktHFqKjoTqyBAREcUKN5LgjmBgxB26SEJgR4aIiEgDQoRzZATOkQHAOTJERESkYwl1R0ZcNEk6dpsldmjP7ZMu4gT4Ltgljqk6ulIVyyjFChVHeh6lOHKxFOOoyEk6PyFUTtJzMKfgOYVzPajJKZzfLlFz0lu7jcec9NpuXYIn4Njewjky0ZFQHRkiIqJY4RaS4BYimCPDVxQA4NASERER6RjvyBAREWnAAwM8EdxP8IC3ZAB2ZIiIiDTBOTLRkVAdGfHNqD6LJvktpCR9cy/g+7ZX/4WTxM/SMoqxQsSRnkcpjlwspThqcpJb3EopJ+k5mFPwnMK5HqKdU0+u8XjMSXftNh5z0mu7NQTWhWJbQnVkiIiIYkXkk305tASwI0NERKSJ7jkyEbw0kkNLADR+amno0KEwGAwBW0VFBQDA4XCgoqICubm5yMzMRFlZGZqbm7WsMhEREcUQTTsy+/fvx8mTJ73bjh07AAA333wzAGDhwoV44403sHHjRuzevRsnTpzArFmztKwyERFRVHjOvWsp3C2SJ57iiaZDSwMGDPD5vGzZMlx88cW4+uqr0dbWhrVr12L9+vWYMmUKAKC2thYjR47E3r17MWHChB7HO9WVGbBPXAFSnBzmXyY79WzAMeIKkOLkMLky/ucJFUfuPP5x5MooxQkWK5ycwvnt5M6TCDmFcz3IlVGKEywWc5KPI40Vi+1WTSy95aTXdutMdQUc21s4RyY6YqY719nZiVdffRV33XUXDAYDGhoa4HK5UFJS4i1TVFSEwsJC1NXVKZ7H6XTCbrf7bERERLHGc+6uSiQbxVBHZvPmzWhtbcUvf/lLAIDNZkNaWhqys7N9ylksFthsNsXzVFdXw2w2e7fBgwf3Yq2JiIhISzHTkVm7di2mTZuGgoKCiM6zaNEitLW1ebempqYo1ZCIiCh63IIh4o0AgyBoP8j29ddf4wc/+AFef/11zJw5EwDwzjvv4JprrsH333/vc1dmyJAhWLBgARYuXKjq3Ha7HWazGUvrS2DMDHyDKxERkcjR7kJV8U60tbXBZDL1Sgzx79K6jy9Hv6zk0AcoOHPajV9e8c9erasexMQdmdraWuTl5WH69OnefePGjUNqaip27drl3Xf06FE0NjbCarVqUU0iIiKKMZoviOfxeFBbW4vy8nKkpJyvjtlsxty5c1FZWYmcnByYTCbce++9sFqtYT2xREREFEs8QhI8ETy15NF+QCUmaN6R2blzJxobG3HXXXcFfLdixQokJSWhrKwMTqcTpaWlWLVqlQa1JCIiii5xPZjwj2dHBoiBjszUqVOhNE3HaDSipqYGNTU1fVwrIiIi0gPNOzJ9qaWzezJUqyvDu09pUaS8tO71Zywp59ehyU1pB6C8+FJz1/nJVv6xgi3y5B9LKY5cLKU4anIS46jJSYyjJicxTqLmJF2MK1ROcguMhboe1OTUk2s8HnPSW7uNx5z02m7PyryJu7d4gIiePPJEryq6llAdGSIiolgR6aJ2XBCvG38FIiIi0i3ekSEiItJA5O9a4r0IgB0ZIiIiTXhggAeRzJHhyr5AgnVkxIlhtrPnJ4s5urp/AqPfBC9xkph0EtrQlFMAgH5JbgDAGY/viozSCXb+sZTiyMUKFUcaSymOmpzEOGpykk5KDJWTdIIdcwqeUzjXg5qcenKNx2NOemu38ZiTXtvtmWQ3+grvyEQHfwUiIiLSrYS6I0NERBQrIl8Qj/ciAHZkiIiINOERDPBEso4M334NIME6MuKiSdKxW0dXqmwZkXTBLnFMNUu8ds59lpZRihUqjvQ8SnHkYoWKEywn6fyEUDlJz8GcgucUzvWgJqcs6b9ZIXLqyTUejznprd3GY056bbduQ9/NkaHoSKiODBERUazwRDi0xAXxurEjQ0REpIHI337NjgzAp5aIiIgSzrJly2AwGLBgwQLvPofDgYqKCuTm5iIzMxNlZWVobm72Oa6xsRHTp09Hv379kJeXhwcffBBdXX33fio57MgQERFpwA1DxFs49u/fjxdeeAFjxozx2b9w4UK88cYb2LhxI3bv3o0TJ05g1qxZ5+vrdmP69Ono7OzEnj178Morr2DdunVYunRpRL9DpBJqaEl8M6rcQlDiPumbewHft716F07yW0hJWiZULKU40vMoxZGLFUlOPgtBhchJeg7mFDyncK4HVTnJLB4Xjd8uHnOK53YbjznFUrt1Cn032VeLoaX29nbMnj0bL730Ep544gnv/ra2Nqxduxbr16/HlClTAAC1tbUYOXIk9u7diwkTJuDtt9/GkSNHsHPnTlgsFowdOxaPP/44Hn74YTz66KNISwu8/voC78gQERHpmN1u99mcTqdi2YqKCkyfPh0lJSU++xsaGuByuXz2FxUVobCwEHV1dQCAuro6jB49GhaLxVumtLQUdrsdn376aZSzUo8dGSIiIg24EenwUrfBgwfDbDZ7t+rqatl4GzZswIEDB2S/t9lsSEtLQ3Z2ts9+i8UCm83mLSPtxIjfi99pJaGGloiIiGJFtIaWmpqaYDKdf79Wenp6QNmmpibcf//92LFjB4xGY9gxY1FCdWSyU88G7BMXThLHVP3LnOrKDDhGXDhJHFOVK+N/nlBx5M7jH0eujFKcYLHCySmc307uPImQUzjXg1wZpTjBYjEn+TjSWLHYbtXE0ltOem23Z9199wROtF4aaTKZfDoychoaGtDS0oIf/ehH5493u/H+++/jD3/4A7Zv347Ozk60trb63JVpbm5Gfn4+ACA/Px/79u3zOa/4VJNYRgscWiIiIopz11xzDQ4dOoSDBw96t/Hjx2P27Nne/05NTcWuXbu8xxw9ehSNjY2wWq0AAKvVikOHDqGlpcVbZseOHTCZTBg1alSf5yRKqDsyREREsUKAAZ4wH6EWj1crKysLl112mc++/v37Izc317t/7ty5qKysRE5ODkwmE+69915YrVZMmDABADB16lSMGjUKc+bMwfLly2Gz2bB48WJUVFTIDmf1FXZkiIiINBCtoaVoWbFiBZKSklBWVgan04nS0lKsWrXK+31ycjK2bNmCefPmwWq1on///igvL0dVVVVU69FT7MgQEREloPfee8/ns9FoRE1NDWpqahSPGTJkCN56661erlnPGARBELSuRG+y2+0wm81YWl8CY2bg21iJiIhEjnYXqop3oq2tLeQE2nCJf5ce+PB6pEfwd8nZ7sLvJ27p1brqAe/IEBERacAd4duvIzk2nvBXICIiIt3iHRkiIiINeAQDPEL4Ty1Fcmw8YUeGiIhIAx4kwRPBwEgkx8aThOrIWFLsAIDclHbvPqXVHZu7uidOtXSen0DV6soAoLyKZF6aXTFWsNUq/WMpxZGLpRRHTU5iHDU5iXHU5CTGSdScpKuKhspJbqXUUNeDmpx6co3HY056a7fxmJNe262z0wXSl4TqyBAREcUKt2CAO4LhoUiOjSfsyBAREWmAc2Sigx0ZIiIiDQgRvv1aiPLKvnqVUB0ZcTx1aMop775+SW4AwBlPsk9ZcWxVOnZrO9u9z9HV/bMZU3zfkiodl/aPpRRHLlaoONJYSnHU5CTGUZOTdCw/VE7ScWnmFDyncK4HNTn15BqPx5z01m7jMSe9tttOV+DvQrEtoToyREREscINA9wRvDQykmPjCTsyREREGvAIkc1z8cT1C4bU03yA7dtvv8Xtt9+O3NxcZGRkYPTo0fjoo4+83wuCgKVLl2LgwIHIyMhASUkJjh07pmGNiYiIKFZo2pH5/vvvMXHiRKSmpmLr1q04cuQIfv/73+OCCy7wllm+fDmee+45rFmzBvX19ejfvz9KS0vhcDg0rDkREVFkPOcm+0aykcZDS08//TQGDx6M2tpa775hw4Z5/1sQBKxcuRKLFy/GzJkzAQB/+tOfYLFYsHnzZtx66609iicumiSdhJYl3tU7t0+6iBPgu2CXODnM0ZWqWEYpVqg40vMoxZGLpRhHRU7SiZahcpKegzkFzymc60FNTuH8domak97abTzmpNd26xI8Acf2Fg8M8EQwzyWSY+OJpt25v//97xg/fjxuvvlm5OXl4YorrsBLL73k/f748eOw2WwoKSnx7jObzSguLkZdXZ3sOZ1OJ+x2u89GRERE8UnTjsxXX32F1atXY/jw4di+fTvmzZuH++67D6+88goAwGazAQAsFovPcRaLxfudv+rqapjNZu82ePDg3k2CiIgoDOLKvpFspPHQksfjwfjx4/HUU08BAK644gocPnwYa9asQXl5eVjnXLRoESorK72f7XY7OzNERBRzIp3nwjky3TT9FQYOHIhRo0b57Bs5ciQaGxsBAPn5+QCA5uZmnzLNzc3e7/ylp6fDZDL5bERERBSfNL0jM3HiRBw9etRn3+eff44hQ4YA6J74m5+fj127dmHs2LEAuu+w1NfXY968eT2OJ74Z1Wf1R78VIaVv7gV83/bqvwKk+FlaRjFWiDjS8yjFkYulFEdNTnKrdCrlJD0HcwqeUzjXQ7Rz6sk1Ho856a7dxmNOem23hsC69BYPInzXEif7AtC4I7Nw4UL8+Mc/xlNPPYVf/OIX2LdvH1588UW8+OKLAACDwYAFCxbgiSeewPDhwzFs2DAsWbIEBQUFuPHGG7WsOhERUUSECJ9aEtiRAaBxR+bKK6/Epk2bsGjRIlRVVWHYsGFYuXIlZs+e7S3z0EMPoaOjA/fccw9aW1sxadIkbNu2DUajUcOaExERRYZvv44OzV9RcP311+P6669X/N5gMKCqqgpVVVV9WCsiIiLSA807Mn3pVFdmwD5x4SRxTNW/THbq2YBjxIWTxDFVuTL+5wkVR+48/nHkyijFCRYrnJzC+e3kzpMIOYVzPciVUYoTLBZzko8jjRWL7VZNLL3lpNd260x1BRzbW/jUUnQkVEeGiIgoVnBoKTrYnSMiIiLd4h0ZIiIiDfBdS9HBjgwREZEGOLQUHQZBEAStK9Gb7HY7zGYzltaXwJgZ+AZXIiIikaPdharinWhra+u1leHFv0vTt/8nUvsHLt6nlqujE2+W/rFX66oHvCNDRESkAd6RiQ52ZIiIiDTAjkx08KklIiIi0q2EuiPT0tk9htjqyvDuU1oUKS/NDgCwpNi9ZXNT2gEoL77U3HV+jNI/VrBFnvxjKcWRi6UUR01OYhw1OYlx1OQkxknUnKSLcYXKSW6BsVDXg5qcenKNx2NOemu38ZiTXtvtWZkXWPYW3pGJjoTqyBAREcUKAZE9Qh3XT+r0ADsyREREGuAdmejgHBkiIiLSLd6RISIi0gDvyERHQnVkxIlhtrPnJ4s5urp/AqPfBC9xkph0EtrQlFMAgH5JbgDAGU+yzzHSCXb+sZTiyMUKFUcaSymOmpzEOGpykk5KDJWTdIIdcwqeUzjXg5qcenKNx2NOemu38ZiTXtvtmWQ3+go7MtHBoSUiIiLSrYS6I0NERBQreEcmOtiRISIi0oAgGCBE0BmJ5Nh4wqElIiIi0q2EuiMjrv4onYTm6EqVLSOSrjwqTg7LEjvB5z5LyyjFChVHeh6lOHKxQsUJlpN0omWonKTnYE7BcwrnelCTU5b0f3yFyKkn13g85qS3dhuPOem13boNfTjZF4aIFsSL5Nh4klAdGSIioljBOTLRwaElIiIi0i3ekSEiItIAJ/tGR0J1ZMQ3o8otBCXuk765F/B926t34SS/hZSkZULFUoojPY9SHLlYkeTksxBUiJyk52BOwXMK53pQlZPM4nHR+O3iMad4brfxmFMstVunwAXx9CahOjJERESxgndkooNzZIiIiEi3eEeGiIhIA0KEQ0u8I9ONHRkiIiINCAAEIbLjKcE6MtmpZwP2iQsniZPD/Muc6soMOEZcOEmcHCZXxv88oeLIncc/jlwZpTjBYoWTUzi/ndx5EiGncK4HuTJKcYLFYk7ycaSxYrHdqomlt5z02m7PugMnClNsS6iODBERUazwwAADV/aNGDsyREREGuBTS9HBp5aIiIhItwyCEMlUo9hnt9thNpuxtL4ExszAl5gRERGJHO0uVBXvRFtbG0wmU6/EEP8uXfbXB5HcLz3s87jPOHH4F7/r1brqAe/IEBERaUAQIt96YvXq1RgzZgxMJhNMJhOsViu2bt3q/d7hcKCiogK5ubnIzMxEWVkZmpubfc7R2NiI6dOno1+/fsjLy8ODDz6Iri5tJ0izI0NERJQABg0ahGXLlqGhoQEfffQRpkyZgpkzZ+LTTz8FACxcuBBvvPEGNm7ciN27d+PEiROYNWuW93i3243p06ejs7MTe/bswSuvvIJ169Zh6dKlWqUEQOOOzKOPPgqDweCzFRUVeb9X0zskIiLSI3GybyRbT8yYMQM///nPMXz4cIwYMQJPPvkkMjMzsXfvXrS1tWHt2rV45plnMGXKFIwbNw61tbXYs2cP9u7dCwB4++23ceTIEbz66qsYO3Yspk2bhscffxw1NTXo7OwMEb33aH5H5tJLL8XJkye92wcffOD9LlTvkIiISK/6uiMj5Xa7sWHDBnR0dMBqtaKhoQEulwslJSXeMkVFRSgsLERdXR0AoK6uDqNHj4bFYvGWKS0thd1u997V0YLmj1+npKQgPz8/YL/YO1y/fj2mTJkCAKitrcXIkSOxd+9eTJgwocexLCl2AEBuSrt3n9KiSM1d3ROnWjrPT6BqdWUAUF58KS/Nrhgr2CJP/rGU4sjFUoqjJicxjpqcxDhqchLjJGpO0sW4QuUkt8BYqOtBTU49ucbjMSe9tdt4zEmv7dbZ6UJf8QgGGKLw9mu73e6zPz09Henp8pOIDx06BKvVCofDgczMTGzatAmjRo3CwYMHkZaWhuzsbJ/yFosFNpsNAGCz2Xw6MeL34nda0fyOzLFjx1BQUIAf/OAHmD17NhobGwFAVe9QjtPphN1u99mIiIji1eDBg2E2m71bdXW1YtlLLrkEBw8eRH19PebNm4fy8nIcOXKkD2sbfZrekSkuLsa6detwySWX4OTJk3jsscfwk5/8BIcPH4bNZgvZO5RTXV2Nxx57rJdrTkREFJlwnjzyPx4AmpqafB6/VrobAwBpaWn44Q9/CAAYN24c9u/fj2effRa33HILOjs70dra6vN3t7m52Ttqkp+fj3379vmcT5y3Kjey0lc0vSMzbdo03HzzzRgzZgxKS0vx1ltvobW1FX/961/DPueiRYvQ1tbm3ZqamqJYYyIioujo7shEMkem+zzi49TiFqwj48/j8cDpdGLcuHFITU3Frl27vN8dPXoUjY2NsFqtAACr1YpDhw6hpaXFW2bHjh0wmUwYNWpUdH6UMGg+R0YqOzsbI0aMwBdffIFrr702ZO9QTrCxQSIiokS1aNEiTJs2DYWFhTh9+jTWr1+P9957D9u3b4fZbMbcuXNRWVmJnJwcmEwm3HvvvbBard45qVOnTsWoUaMwZ84cLF++HDabDYsXL0ZFRYWmf3djqiPT3t6OL7/8EnPmzPHpHZaVlQEI7B32lDgxbGjKKe++fkluAMAZT7JPWXGSmHQSmu1s9z5HV/fPZkzxXQRIOsHOP5ZSHLlYoeJIYynFUZOTGEdNTtJJiaFykk6wY07BcwrnelCTU0+u8XjMSW/tNh5z0mu77XQF/i69pa/ftdTS0oI77rgDJ0+ehNlsxpgxY7B9+3Zce+21AIAVK1YgKSkJZWVlcDqdKC0txapVq7zHJycnY8uWLZg3bx6sViv69++P8vJyVFVVhZ1DNGjakfn1r3+NGTNmYMiQIThx4gQeeeQRJCcn47bbblPVOyQiItIr4dwWyfE9sXbt2qDfG41G1NTUoKamRrHMkCFD8NZbb/Uwcu/StCPzzTff4LbbbsOpU6cwYMAATJo0CXv37sWAAQMAhO4dEhERUWLTtCOzYcOGoN+r6R0SERHpUV8PLcWrmJoj09vERZOkY7dZ4nVwbp90ESfAd8EucUzV0ZWqWEYpVqg40vMoxZGLpRhHRU7S+QmhcpKegzkFzymc60FNTuH8domak97abTzmpNd26xI8Acf2mr4eW4pTCdWRISIiihkR3pEB78gAiIGVfYmIiIjCxTsyREREGojWyr6Jjh0ZIiIiDXCyb3QkVEdGfDOqz6JJfgspSd/cC/i+7dV/4STxs7SMYqwQcaTnUYojF0spjpqc5Ba3UspJeg7mFDyncK6HaOfUk2s8HnPSXbuNx5z02m4NgXWh2JZQHRkiIqKYIRgim7DLOzIA2JEhIiLSBOfIRAefWiIiIiLd4h0ZIiIiLXBBvKhIqI7Mqa7MgH3iCpDi5DD/MtmpZwOOEVeAFCeHyZXxP0+oOHLn8Y8jV0YpTrBY4eQUzm8nd55EyCmc60GujFKcYLGYk3wcaaxYbLdqYuktJ722W2eqK+DY3sKnlqKDQ0tERESkWwl1R4aIiCimcHgoYuzIEBERaYBDS9FhEIT4foDLbrfDbDZjaX0JjJmBb3AlIiISOdpdqCreiba2NphMpl6JIf5dGrzmESRlGMM+j+esA02/eqxX66oHnCNDREREusWhJSIiIk0Yzm2RHE/syBAREWmB68hEBYeWiIiISLcS6o5MS2f3ZKhWV4Z3n9KiSHlpdgCAJcXuLZub0g5AefGl5q7zk638YwVb5Mk/llIcuVhKcdTkJMZRk5MYR01OYpxEzUm6GFeonOQWGAt1PajJqSfXeDzmpLd2G4856bXdnpV5E3ev4R2ZqEiojgwREVHM4Nuvo4JDS0RERKRbvCNDRESkAUHo3iI5nhKsIyOOp9rOnh9jdXR1/wRGv3FRcWxVOnY7NOUUAKBfkhsAcMaT7HOMdFzaP5ZSHLlYoeJIYynFUZOTGEdNTtKx/FA5ScelmVPwnMK5HtTk1JNrPB5z0lu7jcec9NpuzyS70Wc4RyYqOLREREREupVQd2SIiIhiBif7RgU7MkRERBowCN1bJMcTOzJERETa4ByZqEiojoy4aJJ0EpqjK1W2jEi6YJc4OSxLvJt37rO0jFKsUHGk51GKIxcrVJxgOUknWobKSXoO5hQ8p3CuBzU5ZUnvIofIqSfXeDzmpLd2G4856bXdug19ONmXoiKhOjJEREQxg3NkooIdGSIiIi0kyNDSrFmzenzMmjVrkJeXp6osOzJERETUazZv3oxf/OIXyMjICF0YwPr169He3s6ODBERUUxLkDsyAPDcc8+p7pj87//+b4/OnVAdGfHNqHIrWor7pG/uBXzf9updAdJvRUhpmVCxlOJIz6MURy5WJDn5rGgZIifpOZhT8JzCuR5U5SSzCm40frt4zCme22085hRL7dYpcGXfaHv33XeRk5OjuvzWrVtx0UUXqS6fUB0ZIiIi6ltXX311j8pPmjSpR+Vj5hUFy5Ytg8FgwIIFC7z7HA4HKioqkJubi8zMTJSVlaG5uVm7ShIREUWL+NRSJJtOtbS04PDhw/jkk098tnDExB2Z/fv344UXXsCYMWN89i9cuBBvvvkmNm7cCLPZjPnz52PWrFn48MMPNaopERFRdCTiyr4NDQ0oLy/HZ599BuHc67sNBgMEQYDBYIDb3fOhvbA6MrW1tbjlllvQr1+/cA730d7ejtmzZ+Oll17CE0884d3f1taGtWvXYv369ZgyZYo37siRI7F3715MmDChx7GyU88G7BMXThLHVP3LnOrKDDhGXDhJHFOVK+N/nlBx5M7jH0eujFKcYLHCySmc307uPImQUzjXg1wZpTjBYjEn+TjSWLHYbtXE0ltOem23Z92B82soeu666y6MGDECa9euhcVigcEQ+V2lsIaWfvOb3yA/Px9z587Fnj17IqpARUUFpk+fjpKSEp/9DQ0NcLlcPvuLiopQWFiIurq6iGISERFpTojCpjNfffUVli9fjuLiYgwdOhRDhgzx2cIRVkfm22+/xSuvvIJ///vfmDx5MoqKivD000/DZrP16DwbNmzAgQMHUF1dHfCdzWZDWloasrOzffZbLJagcZxOJ+x2u89GRERE2rvmmmvwz3/+M6rnDGtoKSUlBTfddBNuuukmNDc349VXX8Urr7yCJUuW4LrrrsPcuXMxY8YMJCUp95Oamppw//33Y8eOHTAajWEn4K+6uhqPPfZY1M5HRETUGwyIcI5M1GrSd/74xz+ivLwchw8fxmWXXYbUVN/3ZN1www09PmfEk30tFgsmTZqEzz//HJ9//jkOHTqE8vJyXHDBBaitrcXkyZNlj2toaEBLSwt+9KMfefe53W68//77+MMf/oDt27ejs7MTra2tPndlmpubkZ+fr1ifRYsWobKy0vvZbrdj8ODBkaZJREREEaqrq8OHH36IrVu3BnwX7mRfgyBOG+6h5uZm/M///A9qa2vx1Vdf4cYbb8TcuXNRUlKCjo4OVFVVYcOGDfj6669ljz99+nTAd3feeSeKiorw8MMPY/DgwRgwYAD+/Oc/o6ysDABw9OhRFBUVoa6uTvVkX7vdDrPZjKX1JTBmBr6NlYiISORod6GqeCfa2tpgMpl6JYb4d2nIsieRFMGIhMfhwNe/+a9erWu0DR06FNdffz2WLFkCi8USlXOGdUdmxowZ2L59O0aMGIG7774bd9xxh8+qff3798cDDzyA3/3ud4rnyMrKwmWXXeazr3///sjNzfXunzt3LiorK5GTkwOTyYR7770XVqs1rCeWiIiIYkqCrOwrderUKSxcuDBqnRggzI5MXl4edu/eDavVqlhmwIABOH78eNgVA4AVK1YgKSkJZWVlcDqdKC0txapVqyI6JxEREWlj1qxZePfdd3HxxRdH7ZxhdWTWrl2LXbt24be//S1aWlrg8Xh8vn/55ZdhMBh6/CjVe++95/PZaDSipqYGNTU14VSTiIgodiXgHZkRI0Zg0aJF+OCDDzB69OiAyb733Xdfj88ZVkemqqoKjz32GMaPH4+BAwdGZUEbIiKiRJKIK/v+8Y9/RGZmJnbv3o3du3f7fGcwGPquI7N69WqsW7cOc+bMCedwzVhSuteUyU1p9+5TWt2xuat74lRL5/kJVK2uDADKq0jmpZ1fs8Y/VrDVKv1jKcWRi6UUR01OYhw1OYlx1OQkxknUnKSriobKSW6l1FDXg5qcenKNx2NOemu38ZiTXtuts9MF6j2RTjmRE1ZHprOzEz/+8Y+jXRciIqLEkYBDS70hrJV9//M//xPr16+Pdl2IiIgSR4K8oqCyshIdHR2qyy9atAjfffed6vJh3ZFxOBx48cUXsXPnTowZMyZgss4zzzwTzmmJiIgozjz77LNYtGgR+vfvr6p8TU0N7r77bp9lXYIJqyPzySefYOzYsQCAw4cP+3wXyxN/xfHUoSmnvPv6JXWvInjGk+xTVhxblY7d2s5273N0df9sxhTft6RKx6X9YynFkYsVKo40llIcNTmJcdTkJB3LD5WTdFyaOQXPKZzrQU1OPbnG4zEnvbXbeMxJr+220xX4u/SWRJnsKwgCRowYobp/0JO7N0CYHZl33303nMOIiIhIJBi6t0iO14Ha2toeH9OTBfMiftcSERERhSFBJvuWl5f36vnDmuxLREREFAvYkSEiItKAOEcmkq0nqqurceWVVyIrKwt5eXm48cYbcfToUZ8yDocDFRUVyM3NRWZmJsrKytDc3OxTprGxEdOnT0e/fv2Ql5eHBx98EF1dgfOn+kpCDS2JiyZJJ6FliUOM5/ZJF3ECfBfsEieHObpSFcsoxQoVR3oepThysRTjqMhJOtEyVE7SczCn4DmFcz2oySmc3y5Rc9Jbu43HnPTabl2C7yt3elUfDy3t3r0bFRUVuPLKK9HV1YXf/va3mDp1Ko4cOeJ9omjhwoV48803sXHjRpjNZsyfPx+zZs3Chx9+CABwu92YPn068vPzsWfPHpw8eRJ33HEHUlNT8dRTT0WQTPgSqiNDRESUqLZt2+bzed26dcjLy0NDQwN++tOfoq2tDWvXrsX69esxZcoUAN0TdUeOHIm9e/diwoQJePvtt3HkyBHs3LkTFosFY8eOxeOPP46HH34Yjz76KNLS0uRC9yoOLREREWkh0mGlc3dk7Ha7z+Z0OlWFb2trAwDvei0NDQ1wuVwoKSnxlikqKkJhYSHq6uoAAHV1dRg9erTPU0WlpaWw2+349NNPo/Cj9Bw7MkRERFqI0sq+gwcPhtls9m7V1dUhQ3s8HixYsAATJ07EZZddBgCw2WxIS0tDdna2T1mLxQKbzeYt4/9otPhZLNPXEmpoSXyhmM+iSX4LKUlfeAf4viTNf+Ek8bO0jGKsEHGk51GKIxdLKY6anOQWt1LKSXoO5hQ8p3Cuh2jn1JNrPB5z0l27jcec9NpuDdpNWg1XU1MTTKbzi/+lp6eHPKaiogKHDx/GBx980JtV6xO8I0NERKSFKN2RMZlMPluojsz8+fOxZcsWvPvuuxg0aJB3f35+Pjo7O9Ha2upTvrm5Gfn5+d4y/k8xiZ/FMn2NHRkiIiIN9PXj14IgYP78+di0aRPeeecdDBs2zOf7cePGITU1Fbt27fLuO3r0KBobG2G1WgEAVqsVhw4dQktLi7fMjh07YDKZMGrUqPB/jAgk1NASERFRoqqoqMD69evxt7/9DVlZWd45LWazGRkZGTCbzZg7dy4qKyuRk5MDk8mEe++9F1arFRMmTAAATJ06FaNGjcKcOXOwfPly2Gw2LF68GBUVFaqGtHoDOzJEREQJYPXq1QCAyZMn++yvra3FL3/5SwDAihUrkJSUhLKyMjidTpSWlmLVqlXessnJydiyZQvmzZsHq9WK/v37o7y8HFVVVX2VRoCE6sic6soM2CcunCRODvMvk516NuAYceEkcXKYXBn/84SKI3ce/zhyZZTiBIsVTk7h/HZy50mEnMK5HuTKKMUJFos5yceRxorFdqsmlt5y0mu7daa6Ao7tNX28IJ4ghD7AaDSipqYGNTU1imWGDBmCt956q2fBe1FCdWSIiIhiRTjzXPyPJ072JSIiIh3jHRkiIiKt8K5KxNiRISIi0kIfz5GJVwZBzewfHbPb7TCbzVhaXwJjZuAbXImIiESOdheqineira3NZ7XcaBL/Lv3w4aeQnG4M+zxupwNfPP3bXq2rHvCODBERkQY42Tc62JEhIiLSAoeWooJPLREREZFuJdQdmZbO7jHEVleGd5/Sokh5aXYAgCXF7i2bm9IOQHnxpeau82OU/rGCLfLkH0spjlwspThqchLjqMlJjKMmJzFOouYkXYwrVE5yC4yFuh7U5NSTazwec9Jbu43HnPTabs/KvIm7t3BoKToSqiNDREQUMzi0FBUcWiIiIiLd4h0ZIiIiLfCOTFSwI0NERKQBzpGJjoTqyIgTw2xnz08Wc3R1/wRGvwle4iQx6SS0oSmnAAD9ktwAgDOeZJ9jpBPs/GMpxZGLFSqONJZSHDU5iXHU5CSdlBgqJ+kEO+YUPKdwrgc1OfXkGo/HnPTWbuMxJ7222zPJbvQZ3pGJCs6RISIiIt1KqDsyREREMYN3ZKJC0zsyq1evxpgxY2AymWAymWC1WrF161bv9w6HAxUVFcjNzUVmZibKysrQ3NysYY2JiIiiQ5wjE8lGGndkBg0ahGXLlqGhoQEfffQRpkyZgpkzZ+LTTz8FACxcuBBvvPEGNm7ciN27d+PEiROYNWuWllUmIiKiGKLp0NKMGTN8Pj/55JNYvXo19u7di0GDBmHt2rVYv349pkyZAgCora3FyJEjsXfvXkyYMKHH8cTVH6WT0BxdqbJlRNKVR8XJYVmGczvOfZaWUYoVKo70PEpx5GKFihMsJ+lEy1A5Sc/BnILnFM71oCYnbxwVOfXkGo/HnPTWbuMxJ722W7eBk331JmYm+7rdbmzYsAEdHR2wWq1oaGiAy+VCSUmJt0xRUREKCwtRV1eneB6n0wm73e6zERERxRoOLUWH5h2ZQ4cOITMzE+np6fjVr36FTZs2YdSoUbDZbEhLS0N2drZPeYvFApvNpni+6upqmM1m7zZ48OBezoCIiIi0onlH5pJLLsHBgwdRX1+PefPmoby8HEeOHAn7fIsWLUJbW5t3a2pqimJtiYiIokSIwkbaP36dlpaGH/7whwCAcePGYf/+/Xj22Wdxyy23oLOzE62trT53ZZqbm5Gfn694vvT0dKSnp8vHOvdmVLmFoMR90jf3Ar5ve/UunOS3kJK0TKhYSnGk51GKIxcrkpx8FoIKkZP0HMwpeE7hXA+qcpJZPC4av1085hTP7TYec4qldusUOEdGbzS/I+PP4/HA6XRi3LhxSE1Nxa5du7zfHT16FI2NjbBarRrWkIiIiGKFpndkFi1ahGnTpqGwsBCnT5/G+vXr8d5772H79u0wm82YO3cuKisrkZOTA5PJhHvvvRdWqzWsJ5aIiIhiieHcFsnxpHFHpqWlBXfccQdOnjwJs9mMMWPGYPv27bj22msBACtWrEBSUhLKysrgdDpRWlqKVatWaVllIiKi6ODQUlRo2pFZu3Zt0O+NRiNqampQU1PTRzUiIiLqG3z7dXRoPtm3L2Wnng3YJy6cJE4O8y9zqisz4Bhx4SRxcphcGf/zhIojdx7/OHJllOIEixVOTuH8dnLnSYScwrke5MooxQkWiznJx5HGisV2qyaW3nLSa7s96w6cKEyxLaE6MkRERDGDQ0tRwY4MERGRVtgZiVjMPX5NREREpJZBEIS47g/a7XaYzWYsrS+BMTPwJWZEREQiR7sLVcU70dbWBpPJ1CsxxL9Ll93zFJLTjGGfx93pwOEXf9urddUDDi0RERFpgXNkooJDS0RERKRbvCNDRESkAa4jEx3syBAREWmBQ0tRkVAdGUuKHQCQm9Lu3ae0KFJzV/fEqZbO8xOoWl0ZAJQXX8pLsyvGCrbIk38spThysZTiqMlJjKMmJzGOmpzEOImak3QxrlA5yS0wFup6UJNTT67xeMxJb+02HnPSa7t1drpA+pJQHRkiIqJYwaGl6GBHhoiISAscWooKdmSIiIi0wI5MVPDxayIiItKthLojI04MG5pyyruvX5IbAHDGk+xTVpwkJp2EZjvbvc/R1f2zGVN835IqnWDnH0spjlysUHGksZTiqMlJjKMmJ+mkxFA5SSfYMafgOYVzPajJqSfXeDzmpLd2G4856bXddroCf5fewjky0ZFQHRkiIqKYwaGlqODQEhEREekW78gQERFpwCAIMETw3uZIjo0nCdWRERdNko7dZhnO/ce5fdJFnADfBbvEMVVHV6piGaVYoeJIz6MURy6WYhwVOUnnJ4TKSXoO5hQ8p3CuBzU5hfPbJWpOemu38ZiTXtutS/AEHNtrOLQUFRxaIiIiIt1KqDsyREREsYJPLUUHOzJERERa4NBSVHBoiYiIiHQroe7IiG9G9Vk0yW8hJembewHft736L5wkfpaWUYwVIo70PEpx5GIpxVGTk9ziVko5Sc/BnILnFM71EO2cenKNx2NOumu38ZiTXtutIbAuvYVDS9GRUB0ZIiKimMGhpahgR4aIiEgDvCMTHZwjQ0RERLqVUHdkTnVlBuwTF04Sx1T9y2Snng04Rlw4SRxTlSvjf55QceTO4x9HroxSnGCxwskpnN9O7jyJkFM414NcGaU4wWIxJ/k40lix2G7VxNJbTnptt85UV8CxvYZDS1GRUB0ZIiKiWMLhochxaImIiCgBvP/++5gxYwYKCgpgMBiwefNmn+8FQcDSpUsxcOBAZGRkoKSkBMeOHfMp891332H27NkwmUzIzs7G3Llz0d7e3odZBGJHhoiISAuCEPnWAx0dHbj88stRU1Mj+/3y5cvx3HPPYc2aNaivr0f//v1RWloKh8PhLTN79mx8+umn2LFjB7Zs2YL3338f99xzT0Q/Q6Q4tERERKSBvn5qadq0aZg2bZrsd4IgYOXKlVi8eDFmzpwJAPjTn/4Ei8WCzZs349Zbb8Vnn32Gbdu2Yf/+/Rg/fjwA4Pnnn8fPf/5z/Pd//zcKCgrCTyYCBkGI7/eA2+12mM1mLK0vgTEz8A2uREREIke7C1XFO9HW1gaTydQrMcS/S+P/4wmkpBrDPk+Xy4GP/ncxmpqafOqanp6O9PT0oMcaDAZs2rQJN954IwDgq6++wsUXX4yPP/4YY8eO9Za7+uqrMXbsWDz77LN4+eWX8cADD+D7778/X4euLhiNRmzcuBE33XRT2LlEgkNLREREWhCisAEYPHgwzGazd6uuru5xVWw2GwDAYrH47LdYLN7vbDYb8vLyfL5PSUlBTk6Ot4wWOLRERESkAYOne4vkeACyd2QSiaZ3ZKqrq3HllVciKysLeXl5uPHGG3H06FGfMg6HAxUVFcjNzUVmZibKysrQ3NysUY2JiIhii8lk8tnC6cjk5+cDQMDf1+bmZu93+fn5aGlp8fm+q6sL3333nbeMFjTtyOzevRsVFRXYu3cvduzYAZfLhalTp6Kjo8NbZuHChXjjjTewceNG7N69GydOnMCsWbM0rDUREVEURGloKRqGDRuG/Px87Nq1y7vPbrejvr4eVqsVAGC1WtHa2oqGhgZvmXfeeQcejwfFxcXRq0wPaTq0tG3bNp/P69atQ15eHhoaGvDTn/4UbW1tWLt2LdavX48pU6YAAGprazFy5Ejs3bsXEyZM6FG8ls7uW2+trgzvPqXVHfPS7AAAS4rdWzY3pftZeaVVJJu7zt/a848VbLVK/1hKceRiKcVRk5MYR01OYhw1OYlxEjUn6aqioXKSWyk11PWgJqeeXOPxmJPe2m085qTXdntW5k3cvaWvn1pqb2/HF1984f18/PhxHDx4EDk5OSgsLMSCBQvwxBNPYPjw4Rg2bBiWLFmCgoIC74TgkSNH4rrrrsPdd9+NNWvWwOVyYf78+bj11ls1e2IJiLE5Mm1tbQCAnJwcAEBDQwNcLhdKSkq8ZYqKilBYWIi6ujrZjozT6YTT6fR+ttvtAWWIiIg0F8ZaMAHH98BHH32En/3sZ97PlZWVAIDy8nKsW7cODz30EDo6OnDPPfegtbUVkyZNwrZt22A0nn+y6rXXXsP8+fNxzTXXICkpCWVlZXjuuefCzyEKYqYj4/F4sGDBAkycOBGXXXYZgO4Z0mlpacjOzvYpK51F7a+6uhqPPfZYb1eXiIhIVyZPnoxgK64YDAZUVVWhqqpKsUxOTg7Wr1/fG9ULW8w8fl1RUYHDhw9jw4YNEZ1n0aJFaGtr825NTU1RqiEREVH0iENLkWwUI3dk5s+f713qeNCgQd79+fn56OzsRGtrq89dGeksan/BFgISx1NtZ8+PsTq6un8Co9+4qDi2Kh27HZpyCgDQL8kNADjjSfY5Rjou7R9LKY5crFBxpLGU4qjJSYyjJifpWH6onKTj0swpeE7hXA9qcurJNR6POemt3cZjTnptt2eS3egzfPt1VGh6R0YQBMyfPx+bNm3CO++8g2HDhvl8P27cOKSmpvrMoj569CgaGxu9s6iJiIgocWl6R6aiogLr16/H3/72N2RlZXnnvZjNZmRkZMBsNmPu3LmorKxETk4OTCYT7r33Xlit1h4/sURERBRL+vqppXilaUdm9erVALonIEnV1tbil7/8JQBgxYoV3pnRTqcTpaWlWLVqVR/XlIiIKMr6+KmleKVpR0bN+yqNRiNqamoUXztOREREiSsmJvv2FXHRJOkkNEdXqmwZkXTBLnFyWJbh3I5zn6VllGKFiiM9j1IcuVih4gTLSTrRMlRO0nMwp+A5hXM9qMnJG0dFTj25xuMxJ72123jMSa/t1m3ou8m+HFqKjoTqyBAREcUMPrUUFTGzjgwRERFRT/GODBERkQY4tBQd7MgQERFpwSN0b5EcT4nVkRHfjCq3oqW4T/rmXsD3ba/eFSD9VoSUlgkVSymO9DxKceRiRZKTz4qWIXKSnoM5Bc8pnOtBVU4yq+BG47eLx5ziud3GY06x1G6dAlf21RvOkSEiIiLdSqg7MkRERLHCgAjnyEStJvrGjgwREZEWuLJvVCRURyY79WzAPnHhJHFM1b/Mqa7MgGPEhZPEMVW5Mv7nCRVH7jz+ceTKKMUJFiucnML57eTOkwg5hXM9yJVRihMsFnOSjyONFYvtVk0sveWk13Z71h04v4ZiW0J1ZIiIiGIFH7+ODnZkiIiItMCnlqKCTy0RERGRbvGODBERkQYMggBDBBN2Izk2nhgEIb5/CbvdDrPZjKX1JTBmBr6NlYiISORod6GqeCfa2tpgMpl6JYb4d+knP30EKSnGsM/T1eXAP95/rFfrqgccWiIiIiLd4tASERGRBji0FB3syBAREWmBTy1FRUJ1ZCwpdgBAbkq7d5/SokjNXd3jjS2d58cdW10ZAJQXX8pLsyvGCrbIk38spThysZTiqMlJjKMmJzGOmpzEOImak3QxrlA5yS0wFup6UJNTT67xeMxJb+02HnPSa7t1drrQZ7iyb1RwjgwRERHpVkLdkSEiIooVXNk3OtiRISIi0gKHlqKCQ0tERESkWwl1R0acGDY05ZR3X78kNwDgjCfZp6w4SUw6Cc12tnufo6v7ZzOm+L4lVTrBzj+WUhy5WKHiSGMpxVGTkxhHTU7SSYmhcpJOsGNOwXMK53pQk1NPrvF4zElv7TYec9Jru+10Bf4uvcXg6d4iOZ4SrCNDREQUMzi0FBUcWiIiIiLd4h0ZIiIiLXBBvKhgR4aIiEgDfEVBdCRUR0Zc/VE6CS3LcO4/zu2TrkYJ+K48Kk4Oc3SlKpZRihUqjvQ8SnHkYinGUZGTdKJlqJyk52BOwXMK53pQk1M4v12i5qS3dhuPOem13boEzqDVm4TqyBAREcUMTvaNCnZkiIiItCAAiOQGEPsxANiRISIi0gTnyERHQnVkxDej+iya5LeQkvTNvYDv2179F04SP0vLKMYKEUd6HqU4crGU4qjJSW5xK6WcpOdgTsFzCud6iHZOPbnG4zEn3bXbeMxJr+3WEFgXim0J1ZEhIiKKGQIinCMTtZroGjsyREREWuBk36jQdGXf999/HzNmzEBBQQEMBgM2b97s870gCFi6dCkGDhyIjIwMlJSU4NixY9pUloiIiGKOph2Zjo4OXH755aipqZH9fvny5XjuueewZs0a1NfXo3///igtLYXD4ejjmhIREUWZJwobaTu0NG3aNEybNk32O0EQsHLlSixevBgzZ84EAPzpT3+CxWLB5s2bceutt/Y43qmuzIB94sJJ4uQw/zLZqWcDjhEXThInh8mV8T9PqDhy5/GPI1dGKU6wWOHkFM5vJ3eeRMgpnOtBroxSnGCxmJN8HGmsWGy3amLpLSe9tltnqivg2N7Cp5aiI2ZfGnn8+HHYbDaUlJR495nNZhQXF6Ourk7xOKfTCbvd7rMRERFRfIrZjozNZgMAWCwWn/0Wi8X7nZzq6mqYzWbvNnjw4F6tJxERUVjEyb6RbBS7HZlwLVq0CG1tbd6tqalJ6yoREREFYkcmKmK2I5Ofnw8AaG5u9tnf3Nzs/U5Oeno6TCaTz0ZERETxySAIsdGlMxgM2LRpE2688UYA3ZN9CwoK8Otf/xoPPPAAAMButyMvLw/r1q1TPdnXbrfDbDZjaX0JjJmBb3AlIiISOdpdqCreiba2tl77H8Li36VrRj6AlOT0sM/T5XZi12e/79W66oGmTy21t7fjiy++8H4+fvw4Dh48iJycHBQWFmLBggV44oknMHz4cAwbNgxLlixBQUGBt7NDRESkWx4AhgiPJ207Mh999BF+9rOfeT9XVlYCAMrLy7Fu3To89NBD6OjowD333IPW1lZMmjQJ27Ztg9Fo1KrKREREUcHHr6ND047M5MmTEWxky2AwoKqqClVVVX1YKyIiItKLhHrXUktn9xhiqyvDu09pUaS8tO71Zywp59ehyU1pB6C8+FJz1/kxSv9YwRZ58o+lFEcullIcNTmJcdTkJMZRk5MYJ1Fzki7GFSonuQXGQl0PanLqyTUejznprd3GY056bbdnZd7E3Wv4rqWoSKiODBERUczwCIAhgs6Ihx0ZIIYfvyYiIiIKhXdkiIiItMChpajgHRkiIiJNRLqqb3gdmZqaGgwdOhRGoxHFxcXYt29fdNPqYwl1R0acGGY7e36ymKOr+ycw+k3wEieJSSehDU05BQDol+QGAJzxJPscI51g5x9LKY5crFBxpLGU4qjJSYyjJifppMRQOUkn2DGn4DmFcz2oyakn13g85qS3dhuPOem13Z5JdiOe/eUvf0FlZSXWrFmD4uJirFy5EqWlpTh69Cjy8vK0rl5YeEeGiIhICxq8a+mZZ57B3XffjTvvvBOjRo3CmjVr0K9fP7z88su9kGDfYEeGiIhICx4h8g3drzyQbk6nUzZcZ2cnGhoaUFJS4t2XlJSEkpIS1NXV9UnKvYEdGSIiIh0bPHgwzGazd6uurpYt9+9//xtutxsWi8Vnv8Vigc1m64uq9oqEmiMjLpokHbt1dKXKlhFJF+wSx1SzxHdjnPssLaMUK1Qc6XmU4sjFChUnWE7S+QmhcpKegzkFzymc60FNTlnSd7KEyKkn13g85qS3dhuPOem13boNfThHRvB0b5EcD6CpqcnnpZHp6eG/iFKPEqojQ0REFDOi9Pi1yWRS9fbrCy+8EMnJyWhubvbZ39zcjPz8/PDroTEOLREREWkhSnNk1EpLS8O4ceOwa9eu81XweLBr1y5YrdZoZ9dneEeGiIgoQVRWVqK8vBzjx4/HVVddhZUrV6KjowN33nmn1lULGzsyREREWtBgZd9bbrkF//rXv7B06VLYbDaMHTsW27ZtC5gArCcJ1ZER34wqtxCUuE/65l7A922v3oWT/BZSkpYJFUspjvQ8SnHkYkWSk89CUCFykp6DOQXPKZzrQVVOMovHReO3i8ec4rndxmNOsdRunUJfTvZFhB2Z8A6bP38+5s+fH37cGMM5MkRERKRbCXVHhoiIKGbwpZFRwY4MERGRFjweABGsI+OJ4Ng4wqElIiIi0q2EuiOTnXo2YJ+4AqQ4Ocy/zKmuzIBjxBUgxclhcmX8zxMqjtx5/OPIlVGKEyxWODmF89vJnScRcgrnepAroxQnWCzmJB9HGisW262aWHrLSa/t9qw7cKJwr+HQUlQkVEeGiIgoZrAjExUcWiIiIiLd4h0ZIiIiLXgEhL0YjPd4MghCfN+bstvtMJvNWFpfAmNm4NtYiYiIRI52F6qKd6KtrU3VixjDIf5duuaCcqQkBS5iqFaXpxO7vn+lV+uqB7wjQ0REpAWh5y9+DDieOEeGiIiI9It3ZIiIiLQgRDhHhndkALAjQ0REpA2PBzBEsDqvwJV9gQTryFhS7ACA3JR27z6lRZGau7onTrV0np9A1erKAKC8+FJeml0xVrBFnvxjKcWRi6UUR01OYhw1OYlx1OQkxknUnKSLcYXKSW6BsVDXg5qcenKNx2NOemu38ZiTXtuts9MF0peE6sgQERHFDA4tRQU7MkRERBoQPB4IEQwtCRxaAsCnloiIiEjHEuqOjDieOjTllHdfvyQ3AOCMJ9mnrDi2Kh27tZ3t3ufo6v7ZjCm+LxeTjkv7x1KKIxcrVBxpLKU4anIS46jJSTqWHyon6bg0cwqeUzjXg5qcenKNx2NOemu38ZiTXtttpyvwd+k1HFqKioTqyBAREcUMjwAY2JGJFIeWiIiISLd00ZGpqanB0KFDYTQaUVxcjH379mldJSIiosgIQvdaMGFvvCMD6KAj85e//AWVlZV45JFHcODAAVx++eUoLS1FS0uL1lUjIiIKm+ARIt5IB3NknnnmGdx999248847AQBr1qzBm2++iZdffhm/+c1venQucdEk6SS0LMO5/zi3T7qIE+C7YJc4OczRlapYRilWqDjS8yjFkYulGEdFTtKJlqFykp6DOQXPKZzrQU1O4fx2iZqT3tptPOak13br6stHmgUPAK7sG6mYviPT2dmJhoYGlJSUePclJSWhpKQEdXV1ssc4nU7Y7XafjYiIiOJTTHdk/v3vf8PtdsNisfjst1gssNlsssdUV1fDbDZ7t8GDB/dFVYmIiHqEQ0vREdMdmXAsWrQIbW1t3q2pqUnrKhEREQWKaKKvh0NL58T0HJkLL7wQycnJaG5u9tnf3NyM/Px82WPS09ORnp7u/Sycm9XtbO9Ch7F7/LQ9VfL//HNjq+3nOrYdru4yDkf3i8NcHefHWrs6nN3/t6v7+K6Uc2Vw7oVlhvMvG+vo8oulEEcullIcuVhKcdTkJMZRk5MYR01OYpxEzUmMoyandsn/oAqVkxhHTU49ucbjMSfdtdt4zEmn7VYsJ/TBE0FdcEW0Hl4X+IJLADAIffH/rQgUFxfjqquuwvPPPw8A8Hg8KCwsxPz581VN9v3mm284vERERD3S1NSEQYMG9cq5HQ4Hhg0bpjhFoify8/Nx/PhxGI3GKNRMn2L6jgwAVFZWory8HOPHj8dVV12FlStXoqOjw/sUUygFBQU4cuQIRo0ahaamJphMptAHacxut2Pw4MGsby9hfXuX3uoL6K/OrG/vEQQBp0+fRkFBQa/FMBqNOH78ODo7A5+u6qm0tLSE7sQAOujI3HLLLfjXv/6FpUuXwmazYezYsdi2bVvABGAlSUlJuOiiiwAAJpMp5huRFOvbu1jf3qW3+gL6qzPr2zvMZnOvxzAajQnfAYmWmO/IAMD8+fMxf/58ratBREREMSbunloiIiKixJEQHZn09HQ88sgjPk8zxTLWt3exvr1Lb/UF9Fdn1pfovJh/aomIiIhISULckSEiIqL4xI4MERER6RY7MkRERKRb7MgQERGRbiVER6ampgZDhw6F0WhEcXEx9u3bp3WVUF1djSuvvBJZWVnIy8vDjTfeiKNHj/qUcTgcqKioQG5uLjIzM1FWVhbw3imtLFu2DAaDAQsWLPDui7X6fvvtt7j99tuRm5uLjIwMjB49Gh999JH3e0EQsHTpUgwcOBAZGRkoKSnBsWPHNKuv2+3GkiVLMGzYMGRkZODiiy/G448/7vPOFy3r/P7772PGjBkoKCiAwWDA5s2bfb5XU7fvvvsOs2fPhslkQnZ2NubOnYv29vY+r6/L5cLDDz+M0aNHo3///igoKMAdd9yBEydOxGR9/f3qV7+CwWDAypUrY7q+n332GW644QaYzWb0798fV155JRobG73fx9q/GaRPcd+R+ctf/oLKyko88sgjOHDgAC6//HKUlpaipaVF03rt3r0bFRUV2Lt3L3bs2AGXy4WpU6eio6PDW2bhwoV44403sHHjRuzevRsnTpzArFmzNKx1t/379+OFF17AmDFjfPbHUn2///57TJw4Eampqdi6dSuOHDmC3//+97jgggu8ZZYvX47nnnsOa9asQX19Pfr374/S0lI4HA5N6vz0009j9erV+MMf/oDPPvsMTz/9NJYvX+59z5jWde7o6MDll1+Ompoa2e/V1G327Nn49NNPsWPHDmzZsgXvv/8+7rnnnj6v75kzZ3DgwAEsWbIEBw4cwOuvv46jR4/ihhtu8CkXK/WV2rRpE/bu3Su7hH4s1ffLL7/EpEmTUFRUhPfeew+ffPIJlixZ4rOabSz9m0E6JsS5q666SqioqPB+drvdQkFBgVBdXa1hrQK1tLQIAITdu3cLgiAIra2tQmpqqrBx40Zvmc8++0wAINTV1WlVTeH06dPC8OHDhR07dghXX321cP/99wuCEHv1ffjhh4VJkyYpfu/xeIT8/Hzhd7/7nXdfa2urkJ6eLvz5z3/uiyoGmD59unDXXXf57Js1a5Ywe/ZsQRBiq84AhE2bNnk/q6nbkSNHBADC/v37vWW2bt0qGAwG4dtvv+3T+srZt2+fAED4+uuvBUGIzfp+8803wkUXXSQcPnxYGDJkiLBixQrvd7FW31tuuUW4/fbbFY+JtX8zSL/i+o5MZ2cnGhoaUFJS4t2XlJSEkpIS1NXVaVizQG1tbQCAnJwcAEBDQwNcLpdP3YuKilBYWKhp3SsqKjB9+nSfegGxV9+///3vGD9+PG6++Wbk5eXhiiuuwEsvveT9/vjx47DZbD71NZvNKC4u1uz3/fGPf4xdu3bh888/BwD885//xAcffIBp06bFbJ1FaupWV1eH7OxsjB8/3lumpKQESUlJqK+v7/M6+2tra4PBYEB2djaA2Kuvx+PBnDlz8OCDD+LSSy8N+D6W6uvxePDmm29ixIgRKC0tRV5eHoqLi32Gn2Lt3wzSr7juyPz73/+G2+0OeMGkxWKJyuvTo8Xj8WDBggWYOHEiLrvsMgCAzWZDWlqa9x9VkZZ137BhAw4cOIDq6uqA72Ktvl999RVWr16N4cOHY/v27Zg3bx7uu+8+vPLKK976ivWLhfoCwG9+8xvceuutKCoqQmpqKq644gosWLAAs2fPBhCbdRapqZvNZkNeXp7P9ykpKcjJydG8/g6HAw8//DBuu+0270sNY62+Tz/9NFJSUnDffffJfh9L9W1paUF7ezuWLVuG6667Dm+//TZuuukmzJo1C7t37/bWN5b+zSD90sVLI+NdRUUFDh8+jA8++EDrqihqamrC/fffjx07dujija0ejwfjx4/HU089BQC44oorcPjwYaxZswbl5eUa107eX//6V7z22mtYv349Lr30Uhw8eBALFixAQUFBzNY5HrhcLvziF7+AIAhYvXq11tWR1dDQgGeffRYHDhyAwWDQujoheTweAMDMmTOxcOFCAMDYsWOxZ88erFmzBldffbWW1aM4E9d3ZC688EIkJycHzIJvbm5Gfn6+RrXyNX/+fGzZsgXvvvsuBg0a5N2fn5+Pzs5OtLa2+pTXqu4NDQ1oaWnBj370I6SkpCAlJQW7d+/Gc889h5SUFFgslpiq78CBAzFq1CiffSNHjvQ+MSHWKZaujQcffNB7V2b06NGYM2cOFi5c6L0DFot1FqmpW35+fsAk+66uLnz33Xea1V/sxHz99dfYsWOH924MEFv1/cc//oGWlhYUFhZ629/XX3+NBx54AEOHDo25+l544YVISUkJ2QZj6d8M0q+47sikpaVh3Lhx2LVrl3efx+PBrl27YLVaNaxZ96Oq8+fPx6ZNm/DOO+9g2LBhPt+PGzcOqampPnU/evQoGhsbNan7Nddcg0OHDuHgwYPebfz48Zg9e7b3v2OpvhMnTgx4nP3zzz/HkCFDAADDhg1Dfn6+T33tdjvq6+s1uzbOnDmDpCTfJpmcnOz9X7exWGeRmrpZrVa0traioaHBW+add96Bx+NBcXFxn9dZ7MQcO3YMO3fuRG5urs/3sVTfOXPm4JNPPvFpfwUFBXjwwQexffv2mKtvWloarrzyyqBtMNb+jSMd03q2cW/bsGGDkJ6eLqxbt044cuSIcM899wjZ2dmCzWbTtF7z5s0TzGaz8N577wknT570bmfOnPGW+dWvfiUUFhYK77zzjvDRRx8JVqtVsFqtGtbal/SpJUGIrfru27dPSElJEZ588knh2LFjwmuvvSb069dPePXVV71lli1bJmRnZwt/+9vfhE8++USYOXOmMGzYMOHs2bOa1Lm8vFy46KKLhC1btgjHjx8XXn/9deHCCy8UHnrooZio8+nTp4WPP/5Y+PjjjwUAwjPPPCN8/PHH3qd81NTtuuuuE6644gqhvr5e+OCDD4Thw4cLt912W5/Xt7OzU7jhhhuEQYMGCQcPHvRpg06nM+bqK8f/qaVYq+/rr78upKamCi+++KJw7Ngx4fnnnxeSk5OFf/zjH95zxNK/GaRfcd+REQRBeP7554XCwkIhLS1NuOqqq4S9e/dqXSUBgOxWW1vrLXP27Fnh//2//ydccMEFQr9+/YSbbrpJOHnypHaV9uPfkYm1+r7xxhvCZZddJqSnpwtFRUXCiy++6PO9x+MRlixZIlgsFiE9PV245pprhKNHj2pUW0Gw2+3C/fffLxQWFgpGo1H4wQ9+IPzXf/2Xzx9WLev87rvvyl6z5eXlqut26tQp4bbbbhMyMzMFk8kk3HnnncLp06f7vL7Hjx9XbIPvvvtuzNVXjlxHJtbqu3btWuGHP/yhYDQahcsvv1zYvHmzzzli7d8M0ieDIEiWDSUiIiLSkbieI0NERETxjR0ZIiIi0i12ZIiIiEi32JEhIiIi3WJHhoiIiHSLHRkiIiLSLXZkiIiISLfYkSEiIiLdYkeGiIiIdIsdGSIiItItdmSI4tjkyZNx33334aGHHkJOTg7y8/Px6KOPAgDee+89pKWl4R//+Ie3/PLly5GXl4fm5maNakxE1DPsyBDFuVdeeQX9+/dHfX09li9fjqqqKuzYsQOTJ0/GggULMGfOHLS1teHjjz/GkiVL8Mc//hEWi0XrahMRqcKXRhLFscmTJ8PtdvvcdbnqqqswZcoULFu2DJ2dnSguLsaIESNw+PBhTJw4ES+++KKGNSYi6pkUrStARL1rzJgxPp8HDhyIlpYWAEBaWhpee+01jBkzBkOGDMGKFSu0qCIRUdg4tEQU51JTU30+GwwGeDwe7+c9e/YAAL777jt89913fVo3IqJIsSNDlMC+/PJLLFy4EC+99BKKi4tRXl7u08khIop17MgQJSi3243bb78dpaWluPPOO1FbW4tPPvkEv//977WuGhGRauzIECWoJ598El9//TVeeOEFAN1zZ1588UUsXrwY//znPzWuHRGROnxqiYiIiHSLd2SIiIhIt9iRISIiIt1iR4aIiIh0ix0ZIiIi0i12ZIiIiEi32JEhIiIi3WJHhoiIiHSLHRkiIiLSLXZkiIiISLfYkSEiIiLdYkeGiIiIdIsdGSIiItKt/x9TaKFqdzndXgAAAABJRU5ErkJggg==",
      "text/plain": [
       "<Figure size 640x480 with 2 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "topo.depth.plot(vmin=-50, vmax=500)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### *Apply ridges*"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [],
   "source": [
    "# apply a submerged ridge ( height of 300 m):\n",
    "topo.apply_ridge(height=400.0, width=4.0, lon=100, ilat=(10,70) )"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [],
   "source": [
    "# apply a ridge reaching above MSL ( height of 600 m):\n",
    "topo.apply_ridge(height=600.0, width=16.0, lon=240, ilat=(10,70) )"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.collections.QuadMesh at 0x15426892c5d0>"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjIAAAG2CAYAAABlBWwKAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8ekN5oAAAACXBIWXMAAA9hAAAPYQGoP6dpAABh7klEQVR4nO3dfXwU5b03/s8mmyeS7Aai2RAhgBWKKPiAGlI5p1ZTUm7qE9yt+kNMleqvnIACrVruI+ih1SDnrlB6IqjlgD2Vm0pvwSMWLaBglQQ0aAWxiIKA4CYKJiEhT7s79x+zs8zszmRnd2dnd3Y/79drX8nOzl7f67u5ZrmY65prbIIgCCAiIiKyoIxEV4CIiIgoWuzIEBERkWWxI0NERESWxY4MERERWRY7MkRERGRZ7MgQERGRZbEjQ0RERJbFjgwRERFZFjsyREREZFnsyBAREZFlsSNDRESUBh577DHYbDbFY/To0YHXu7u7UVtbi+LiYhQUFGDatGlobm5WlHHs2DFMmTIFAwYMQElJCR588EF4PB6zU1GwJzQ6ERERmeaSSy7Btm3bAs/t9nPdgHnz5uHVV1/Fhg0b4HQ6MXv2bEydOhXvvPMOAMDr9WLKlCkoLS3Frl278OWXX+Kuu+5CVlYWnnjiCdNzkdh400giIqLU99hjj2HTpk344IMPQl5ra2vD+eefj3Xr1uF//s//CQD4xz/+gYsvvhgNDQ2YMGECtmzZgh/+8Ic4efIkXC4XAGDVqlV4+OGH8dVXXyE7O9vMdAJS/oyMz+fDyZMnUVhYCJvNlujqEBFREhMEAWfOnEFZWRkyMuIz+2LTpk04evSoIWXdeuutKCoqUmzLyclBTk6O6v6HDh1CWVkZcnNzUVlZibq6OpSXl6OpqQl9fX2oqqoK7Dt69GiUl5cHOjINDQ0YO3ZsoBMDANXV1Zg1axY++ugjXHHFFYbkFKmU78icPHkSQ4cOTXQ1iIjIQo4fP44hQ4YYXu6ZM2dw6623Ysr385EV47/Af/+oBy+++CJ27dql2P7oo4/iscceC9m/oqICa9euxbe//W18+eWX+Ld/+zf80z/9E/bv3w+3243s7OyQTpHL5YLb7QYAuN1uRSdGel16LVFSviNTWFgIAHh4+3U4bisDALR0FwZe7/aIH0GuXZysVJJ7BgBwUX4LAGB0zrk/Trn9NABgQIYXAHDWlwkAOOYZBAD4R09pYN9PO0sUsbTiqMXSiqMWSyuOnpykOHpykuLoyUmKk645SXH05CTF0ZOTFEdPTpG08VTMyWrHbSrmZNXj9myHFz+d+HHg3w6jSbM5nv+P8+F0ZIbZu38P/OtXyC0Sh3zktM7GTJ48OfD7uHHjUFFRgWHDhuHFF19EXl5eTHVJpJTvyEjDSTkFdmTbxPE7e+a5P7LdfwDY7WKDys7rAQDk5mcBAAbknmtoBXbxNOOADLEhZvj8zz3iPrlZWYF9g2NpxVGLpRVHLZZWHD05SXH05CTF0ZOTFCddc5Li6MlJiqMnJymOnpwiaeOpmJPVjttUzMnKxy0Ay0xFyMzMhMPhiOq9RUVFGDVqFD799FN8//vfR29vL1pbWxVnZZqbm1FaKnYiS0tLsWfPHkUZ0lVN0j6JwMuviYiITOYVBHgFX0yPWC/V6ejowGeffYbBgwdj/PjxyMrKwvbt2wOvHzx4EMeOHUNlZSUAoLKyEvv27UNLy7kzXVu3boXD4cCYMWNiq0wMUv6MDBERUbLxQYAPsfVEhAjf/4tf/AI33ngjhg0bhpMnT+LRRx9FZmYm7rjjDjidTsycORPz58/HoEGD4HA4MGfOHFRWVmLChAkAgEmTJmHMmDGYMWMGli5dCrfbjUceeQS1tbWaw1lmSKuOTK9PTFc+dtvtyVLdR5Jr6w38Lo2pFkpnHP3P5ftoxQoXR16OVhy1WOHi9JeTfH5CuJzkZTCn/nOKpj3oyalQfqY7TE6RtPFUzMlqx20q5mTV49Zr84a8N1V88cUXuOOOO3Dq1Cmcf/75mDhxIhobG3H++ecDAJYtW4aMjAxMmzYNPT09qK6uxtNPPx14f2ZmJjZv3oxZs2ahsrIS+fn5qKmpweLFixOVEoA068gQERElAx988MVYRqRnZNavX9/v67m5uaivr0d9fb3mPsOGDcNf/vKXiOLGGzsyREREJhPnyMQ6tEQAOzJERESmS8QcmVSVVh2Z7AxxjQH5OhcSaZu0j6RbOHepX2ANhaD1B+T7hIulFUdejlYctVix5CQvN1xO8jKYU/85RdMedOWksuaKEZ9dKuaUysdtKuaUTMdtj5C6c2RSVVp1ZIiIiJKBDwK8MZ+RIYAdGSIiItMZM7READsyREREpjNksm+sK+KlCK7sS0RERJaVVmdkirK6QrZJCydJk8OC9znlKQh5j7RwkjQ5TG2f4HLCxVErJziO2j5acfqLFU1O0Xx2auWkQ07RtAe1fbTi9BeLOanHkcdKxuNWTyyr5WTV47bLGzpROB58/kcseD5GlFYdGSIiomTg5WRfw3BoiYiIiCyLZ2SIiIhM5hXERyxiHZpKFezIEBERmUwA58gYxSak+PVb7e3tcDqdWLS7CrkFoXdjJSIiknR39GFxxTa0tbXB4XAYXr70b1LTRy4UOmKb3bF4YRsGn1+LpUuXGlQ7a+IcGSIiIrIsDi0RERGZzCeIj1ik9HBKBNiRISIiMpkXNnhhi6kMIcb3p4q06si47O0AgGJ7R2Cb1qJIzR5xbLSl99wYaWtfHgDtxZdKsts1Y/W3yFNwLK04arG04ujJSYqjJycpjp6cpDjpmpN8Ma5wOaktMBauPejJKZI2noo5We24TcWcrHrc9vT2wQzGdGQI4BwZIiIisrC0OiNDRESUDHyCDT4hxjMyPCUDgB0ZIiIi03GOjHE4tERERESWlVZnZKSJYcPtpwLbBmR4AQBnfZmKfaVJYvJJaO4ucVu3R/zYcu3Ku6TKJ9gFx9KKoxYrXBx5LK04enKS4ujJST4pMVxO8gl2zKn/nKJpD3pyiqSNp2JOVjtuUzEnqx63vX2hn0s8+JABb4znEnw8IwMgzToyREREycCQOTIG1cXq2JEhIiIyGefIGIdzZIiIiMiy0uqMjLRoknzstlDq0Pq3yRdxApQLdkljqt2eLM19tGKFiyMvRyuOWizNODpyks9PCJeTvAzm1H9O0bQHPTlF89mla05WO25TMSerHrd9Qqz3pNbHK9jgFWI7l8DLr0Vp1ZEhIiJKBj5kwBfjoAiHlkTsyBAREZmMc2SMwzkyREREZFk8I0NERGQyn5AR+xwZg+pidWnVkZHujKpYNCloISX5nXsB5d1egxdOkp7L99GMFSaOvBytOGqxtOLoyUltcSutnORlMKf+c4qmPRidUyRtPBVzstxxm4o5WfW4tYXWJR58sMW8oB2HlkQcWiIiIiLLSqszMkRERMnAy1sUGCahZ2SGDx8Om80W8qitrQUAdHd3o7a2FsXFxSgoKMC0adPQ3NycyCoTERHFzAtxjkwsDyHGWxykioR2ZN599118+eWXgcfWrVsBAD/60Y8AAPPmzcMrr7yCDRs2YOfOnTh58iSmTp2ayCoTERHFTFpHJpYHJ/uKEjq0dP755yueL1myBN/61rfw3e9+F21tbVi9ejXWrVuH66+/HgCwZs0aXHzxxWhsbMSECRMijnfKUxCyTVoBUpocFrxPUVZXyHukFSClyWFq+wSXEy6OWjnBcdT20YrTX6xocorms1MrJx1yiqY9qO2jFae/WMxJPY48VjIet3piWS0nqx63PVl9Ie+l5JY0k317e3vxxz/+Effccw9sNhuamprQ19eHqqqqwD6jR49GeXk5GhoaNMvp6elBe3u74kFERJRMxFsUxPbgVUuipOnIbNq0Ca2trfjJT34CAHC73cjOzkZRUZFiP5fLBbfbrVlOXV0dnE5n4DF06NA41pqIiChyPv9k31ge7MiIkqYjs3r1akyePBllZWUxlbNgwQK0tbUFHsePHzeohkRERMbwCTb4hIyYHuzIiJLi8uujR49i27ZteOmllwLbSktL0dvbi9bWVsVZmebmZpSWlmqWlZOTg5ycHNXXmj0OxU89SrLbVX8PJ5ZYyRzHjFj/sWya4vmP5myNSxwzc0rm9mBmLLPjAMC2Y98GABSsd/h/NgIAvgjskS97p/j7mf+vEgBw9rZWAMAd33ovbCwzc0q1v1MyHbfd2ZwjYzVJcUZmzZo1KCkpwZQpUwLbxo8fj6ysLGzfvj2w7eDBgzh27BgqKysTUU0iIiJDxDqs5OVVSwEJPyPj8/mwZs0a1NTUwG4/Vx2n04mZM2di/vz5GDRoEBwOB+bMmYPKysqorlgiIiJKFtKE3Vj4uI4MgCToyGzbtg3Hjh3DPffcE/LasmXLkJGRgWnTpqGnpwfV1dV4+umnE1BLIiIiSkYJ78hMmjQJgqB+giw3Nxf19fWor683uVZERETxIy1qFwtO9hUlvCNjppZecWJYa19eYJvWokjShDCX/dzEsGJ7BwDtxZfkE8+CY/W3yFNwLK04arG04ujJSYqjJycpjp6c5JPposnJ3iV2bDN7xJ/7z1xgqZzki3GF+zupLTAWrj3oySmSNp6KOcmP21KH+PtXxUUAgNCl0kLlnBbLdX8zAADwSadLEUdPTpEct5HmxO8i6MopmjbepXIn7niQbjMQC3ZkRGnVkSEiIkoGPthivukjJ/uKkuKqJSIiIqJo8IwMERGRyXxGDC3xqiUA7MgQERGZTloLJhacIyNKq46MNDHM3XVusli3R/wIcoMmeEmTxOST0IbbTwEABmR4AQBnfZmK98gn2AXH0oqjFitcHHksrTh6cpLi6MlJPikxXE7yCXbR5CRN8s1uF/c92Xnuc7VqTka2Bz05RdLGUzEn+XFbmncGAHCiUHye6XQCALxtbSH1lGS3+1d3PZMbEltvTpF8dpHmxO8i6MopmjZ+NtMLM4i3KOAcGSNwjgwRERFZVlqdkSEiIkoGHFoyDjsyREREJpPufh0LTvYVpVVHRlo0ST522+3JUt1HIl+wSxpTLZTajv+5fB+tWOHiyMvRiqMWK1yc/nKSz08Il5O8jHjnZPO/lOERLJlTNO1BT06F8u+sMDlF0sZTMSd5GXmZ4nwXn/8ttpxshGPziHXJ8NgUseN13Iq/87tIT07xPm69NnPmyCTakiVLsGDBAjzwwANYvnw5AKC7uxs///nPsX79esVtgVwuV+B9x44dw6xZs/Dmm2+ioKAANTU1qKurU9wr0WycI0NERGQyL2wxP6Kd7Pvuu+/imWeewbhx4xTb582bh1deeQUbNmzAzp07cfLkSUydOvVcnb1eTJkyBb29vdi1axeef/55rF27FosWLYrhk4gdOzJEREQm8wkZMT+imSPT0dGB6dOn47nnnsPAgQMD29va2rB69Wo89dRTuP766zF+/HisWbMGu3btQmNjIwDgr3/9Kw4cOIA//vGPuPzyyzF58mT86le/Qn19PXp7Q88GmoUdGSIiIpP5DDgj44MNXq8X7e3tikdPT49m3NraWkyZMgVVVVWK7U1NTejr61NsHz16NMrLy9HQ0AAAaGhowNixYxVDTdXV1Whvb8dHH31k8CekHzsyREREFtXY2Ain06l41NXVqe67fv167N27V/V1t9uN7OxsFBUVKba7XC643e7APvJOjPS69FqipNVkX+nOqGoLQUnb5HfuBZR3ew0snBS0kJJ8n3CxtOLIy9GKoxYrlpwUC0GFyUleRrxzEvwv+ew2f/l9lsopmvagKyeVxeOMaA+pmJO8jC6vOBE0w9+MhJ7wp8AFu1gXn13oN448llnHbX+x+F0UexvvEcxaEC8j5quWAGDChAnYsmWLYltOTk7IfsePH8cDDzyArVu3Ijc3N+a4yYRnZIiIiEzm9d9rKZaHT7AhMzMTDodD8VDryDQ1NaGlpQVXXnkl7HY77HY7du7ciRUrVsBut8PlcqG3txetra2K9zU3N6O0tBQAUFpaiubm5pDXpdcShR0ZIiIik/n8c1xieUQy2feGG27Avn378MEHHwQeV111FaZPnx74PSsrC9u3bw+85+DBgzh27BgqKysBAJWVldi3bx9aWloC+2zduhUOhwNjxowx7sOJUFoNLREREaWjwsJCXHrppYpt+fn5KC4uDmyfOXMm5s+fj0GDBsHhcGDOnDmorKzEhAkTAACTJk3CmDFjMGPGDCxduhRutxuPPPIIamtrVc8CmSWtOjJFWV0h26SFk6Qx1eB9TnkKQt4jLZwkjamq7RNcTrg4auUEx1HbRytOf7GiySmaz06tHD05eXPE/2X0OsSx6wvzz934zQo5RdMe1PbRitNfLOYUGsfdJd4tMlu8d2S/N4sMlO/wL7BWKE6sKc1r16yLEcctwO8irVhmt/Eub+j8mniQhodiYfRNI5ctW4aMjAxMmzZNsSCeJDMzE5s3b8asWbNQWVmJ/Px81NTUYPHixQbXJDJp1ZEhIiJKBobc/TrG9+/YsUPxPDc3F/X19aivr9d8z7Bhw/CXv/wlprhG4xwZIiIisiyekSEiIjKZj3e/Ngw7MkRERCYzZGjJoLpYXVp1ZEqy2xU/9Wj2OFR/j2esZI5jRixPnk3x89LCE3GJY2ZOydwezIxldhwAcLeLvxec8ul+f88g8avRMbAVADAqv1lz30TklGp/p2Q6brs9ff2+bhTxEmqekTEC58gQERGRZaXVGRkiIqJk4BVs8MZ81ZJBlbE4dmSIiIhMZswcGQ4tAezIEBERmc6Im0ayIyNKq46Myy5ODCu2dwS2aa3uKE0Ia+k9NzGstS8PgPYqkvKJZ8Gx+lutMjiWVhy1WFpx9OQkxdGTkxRHT07yyXTR5CRVy+dfRPTKvM8tlZN8VdFwfye1lVLDtQc9OUXSxlMxJ/lxe6ZNjHXe1/pXbO0uEv+BGFb0DQBgTO5JRRw9OUVy3EaaE7+LoCunaNp4T685k33JOGnVkSEiIkoGXtjgjfGMCqfIiNiRISIiMpkPib9FQarg5ddERERkWWl1RkYaTx1uPxXYNiDDCwA468tU7CuNrcrHbt1d4rZuj/ix5dqVY+7ycengWFpx1GKFiyOPpRVHT05SHD05ycfyw+UkH5eOJidpbownX/x5YVar5XMysj3oySmSNp6KOcmPW3T650J80wkA0LMsXp+/7V0woDUktt6cIvnsIs2J30XQlVM0bby3L/RziQfBkMm+BKRZR4aIiCgZiCv7xjo0xKElgB0ZIiIi0xmxIJ6Pp2QAJMEcmRMnTuDOO+9EcXEx8vLyMHbsWLz33nuB1wVBwKJFizB48GDk5eWhqqoKhw4dSmCNiYiIKFkktCPzzTff4Nprr0VWVha2bNmCAwcO4De/+Q0GDhwY2Gfp0qVYsWIFVq1ahd27dyM/Px/V1dXo7u5OYM2JiIiiJy2IF8uDC+KJEjq09OSTT2Lo0KFYs2ZNYNuIESMCvwuCgOXLl+ORRx7BzTffDAD4wx/+AJfLhU2bNuH222+PKJ60aJJ8Elqh1A782+SLOAHKBbukyWHdnizNfbRihYsjL0crjloszTg6cpJPtAyXk7yMeOck+OfaSZutllM07UFPTtG0h3TNSV6GzSPubOsNvyCeLVNsfIK/WgV2Zex4HbcAv4skiT5u+wT9d0mPhRG3KAAvvwaQ4DMy//3f/42rrroKP/rRj1BSUoIrrrgCzz33XOD1I0eOwO12o6qqKrDN6XSioqICDQ0NqmX29PSgvb1d8SAiIkom0mTfWB6cIiNKaEfm8OHDWLlyJUaOHInXX38ds2bNwv3334/nn38eAOB2uwEALpdL8T6XyxV4LVhdXR2cTmfgMXTo0PgmQURERAmT0I6Mz+fDlVdeiSeeeAJXXHEF7rvvPtx7771YtWpV1GUuWLAAbW1tgcfx48cNrDEREVHspKGlWB48IyNKaEdm8ODBGDNmjGLbxRdfjGPHjgEASktLAQDNzc2KfZqbmwOvBcvJyYHD4VA8iIiIkonYGeFkXyMkdLLvtddei4MHDyq2ffLJJxg2bBgAceJvaWkptm/fjssvvxwA0N7ejt27d2PWrFkRx5PujKpY/TFoRUj5nXsB5d1eg1eAlJ7L99GMFSaOvBytOGqxtOLoyUltlU6tnORlxDsnmz8FafMZ2X87rJBTNO3B6JwiaeOpmJO8DMEuNiAhO/zXneAVy7X5b4Dc4VHGjttxC/C7KEjCjlub/rukU3JIaEdm3rx5+M53voMnnngCP/7xj7Fnzx48++yzePbZZwEANpsNc+fOxa9//WuMHDkSI0aMwMKFC1FWVoZbbrklkVUnIiKKmmDAVUu8aaQooR2Zq6++Ghs3bsSCBQuwePFijBgxAsuXL8f06dMD+zz00EPo7OzEfffdh9bWVkycOBGvvfYacnNzE1hzIiKi6BlxiwLOkREl/BYFP/zhD/HDH/5Q83WbzYbFixdj8eLFJtaKiIgofoxYR4ZzZEQJ78iY6ZSnIGSbtHCSNKYavE9RVlfIe6SFk6QxVbV9gssJF0etnOA4avtoxekvVjQ5RfPZqZWjJ6cMfwpS4zzcVxR4zQo5RdMe1PbRitNfLOakEidfLLd3oHgWV8+XXpZ4o2ycOFsEAPjcU6xZFyOOW4DfRVqxzG7jPVl9Ie+l5JZWHRkiIqJkYMzKvsbUxerYkSEiIjKZMUNLBLAjQ0REZDpDOjK8aglAghfEIyIiIopFWp2RafY4FD/1KMluV/09nrGSOY4ZsQJz8vw/93YNj0scM3NK5vZgZiyz4wBAoVNsSN3nie8PnQYaKrdVPGl/tHUgAOBAcVnYWGbmlGp/p2Q6bruzzZnsy8uvjZNWHRkiIqJkwMuvjcOhJSIiIrIsnpEhIiIymSG3KDCoLlaXVh2Zll5xPLW1Ly+wTWtRJGkc1WU/N55abO8AoL34kny8NjhWf4s8BcfSiqMWSyuOnpykOHpykuLoyUk+Bh1NTvYu8fDM7BF/7j9zgaVyki/GFe7vpLbAWLj2oCenSNp4KuYkP25LHeLvXxUXAdA3RybntFiu+5sBAIBPOl2KOHpyiuS4jTQnfhdBV07RtPEulRtYxoMx68hwaAlIs44MERFRMuA6MsbhHBkiIiKyLJ6RISIiMpkg2GJe0E7gKRkA7MgQERGZzph1ZDhHBkizjow0MczddW6yWLdH/AhygyZ4SZPE5JPQhttPAQAGZHgBAGd9mYr3yCfYBcfSiqMWK1wceSytOHpykuLoyUk+KTFcTvIJdtHkJE3yzW4X9z3Zee5ztWpORrYHPTlF0sZTMSf5cVuadwYAcKJQfJ7pdAIAvG1tIfWUZLf7F0U7kxsSW29OkXx2kebE7yLoyimaNn420wsz+GDAZF8CwDkyREREZGFpdUaGiIgoGRgyR4ZDSwDYkSEiIjKdMXe/NqgyFsehJSIiIrKstDojI63+KJ+E1u3JUt1HIl95VJocVih1ov3P5ftoxQoXR16OVhy1WOHi9JeTfKJluJzkZcQ7J5v/pQyPYMmcomkPenIqlP/nLUxOkbTxVMxJXkZepjhx1+d/iy0nG+HYPGJdMjw2Rex4Hbfi7/wu0pNTvI9br82cyb6CgJiHlkiUVh0ZIiKiZGDM0BI7QgA7MkRERKYTz8gkuhapgXNkiIiIyLLS6oyMdGdUtYWgpG3yO/cCyru9BhZOClpISb5PuFhaceTlaMVRixVLToqFoMLkJC8j3jkJ/pd8dpu//D5L5RRNe9CVk8ricUa0h1TMSV5Gl1ecP5Hhb0ZCT+g8kmCCXayLzy70G0cey6zjtr9Y/C6KvY33CCbNkTFkZV8C0qwjQ0RElAyMudcS58gA7MgQERGZzifYYIt5QTwCOEeGiIiILIwdGSIiIpNJVy3F8ojUypUrMW7cODgcDjgcDlRWVmLLli2B17u7u1FbW4vi4mIUFBRg2rRpaG5uVpRx7NgxTJkyBQMGDEBJSQkefPBBeDyhc5LMlFZDS0VZXSHbpIWTpMlhwfuc8hSEvEdaOEmaHKa2T3A54eKolRMcR20frTj9xYomp2g+O7Vy9OTkzRFPt/Y6xEl4F+afu4OtFXKKpj2o7aMVp79YzCk0jrtLvO11tngT7H7veh0o3+FfYK1QnCFcmteuWRcjjluA30Vascxu411ec/5RNmKODCJ8/5AhQ7BkyRKMHDkSgiDg+eefx80334z3338fl1xyCebNm4dXX30VGzZsgNPpxOzZszF16lS88847AACv14spU6agtLQUu3btwpdffom77roLWVlZeOKJJ2LLJQZp1ZEhIiJKVzfeeKPi+eOPP46VK1eisbERQ4YMwerVq7Fu3Tpcf/31AIA1a9bg4osvRmNjIyZMmIC//vWvOHDgALZt2waXy4XLL78cv/rVr/Dwww/jscceQ3Z2+FWz44FDS0RERCaTzsjE9IB4lqS9vV3x6OnpCRvf6/Vi/fr16OzsRGVlJZqamtDX14eqqqrAPqNHj0Z5eTkaGhoAAA0NDRg7dixcLldgn+rqarS3t+Ojjz4y/DPSix0ZIiIik0m3KIjlIQhAY2MjnE6n4lFXV6cZd9++fSgoKEBOTg5+9rOfYePGjRgzZgzcbjeys7NRVFSk2N/lcsHtdgMA3G63ohMjvS69lihpNbRUkt2u+KlHs8eh+ns8YyVzHDNiefJsip+XFp6ISxwzc0rm9mBmLLPjAIC7Xfy94JRP9/t7BolfjY6BrQCAUfnNmvsmIqdU+zsl03Hb7enr93WjGHWLggkTJigm7AJATk6O5v7f/va38cEHH6CtrQ1//vOfUVNTg507d8ZekQRKq44MERFRKsnMzITDob8jmJ2djYsuuggAMH78eLz77rv47W9/i9tuuw29vb1obW1VnJVpbm5GaWkpAKC0tBR79uxRlCdd1STtkwgcWiIiIjKZIXNkDFjZ1+fzoaenB+PHj0dWVha2b98eeO3gwYM4duwYKisrAQCVlZXYt28fWlpaAvts3boVDocDY8aMibku0UpoR+axxx6DzWZTPEaPHh14Xc817URERFaTiE7MggUL8NZbb+Hzzz/Hvn37sGDBAuzYsQPTp0+H0+nEzJkzMX/+fLz55ptoamrC3XffjcrKSkyYMAEAMGnSJIwZMwYzZszA3//+d7z++ut45JFHUFtb2+9wVrwlfGjpkksuwbZt2wLP7fZzVQp3TTsREZEVCTD/FgMtLS2466678OWXX8LpdGLcuHF4/fXX8f3vfx8AsGzZMmRkZGDatGno6elBdXU1nn766cD7MzMzsXnzZsyaNQuVlZXIz89HTU0NFi9ebHImSgnvyNjtdtWxtba2trDXtEfKZRcnhhXbOwLbtBZFkiaEtfSeG3ts7csDoL34knziWXCs/hZ5Co6lFUctllYcPTlJcfTkJMXRk5N8Ml00OUnV8vmXJLgy73NL5SRfjCvc30ltgbFw7UFPTpG08VTMSX7cnmkTY533tf6FzrqLxP/tDiv6BgAwJvekIo6enCI5biPNid9F0JVTNG28p9ecyb6JsHr16n5fz83NRX19Perr6zX3GTZsGP7yl78YXbWYJHyOzKFDh1BWVoYLL7wQ06dPx7FjxwBA1zXtanp6ekKuqSciIkomxqwjw7tfAwnuyFRUVGDt2rV47bXXsHLlShw5cgT/9E//hDNnzui6pl1NXV2d4nr6oUOHxjkLIiKiCAkGPSixQ0uTJ08O/D5u3DhUVFRg2LBhePHFF5GXl9fPO7UtWLAA8+fPDzxvb29nZ4aIiChFJXyOjFxRURFGjRqFTz/9FN///vfDXtOuJicnJ6Gzp4mIiMIRYMzl05RkHZmOjg589tlnmDFjhuKa9mnTpgEIvaY9UtLEsOH2U4FtAzK8AICzvkzFvtIkMfkkNHeXuK3bI35suXbl5EH5BLvgWFpx1GKFiyOPpRVHT05SHD05ySclhstJPsEumpykSb6efPHnhVmtls/JyPagJ6dI2ngq5iQ/btHpn9T5TScAQM/6vn3+tnfBgNaQ2HpziuSzizQnfhdBV07RtPHevtDPJR6MWNnXiJWBU0FCOzK/+MUvcOONN2LYsGE4efIkHn30UWRmZuKOO+5QXNM+aNAgOBwOzJkzR3FNOxERkSUZsqAdz+gACe7IfPHFF7jjjjtw6tQpnH/++Zg4cSIaGxtx/vnnAwh/TTsRERGlt4R2ZNavX9/v63quaSciIrIcwSY+YirDmKpYXVLNkYk3adEk+dhtodSO/NvkizgBygW7pDHVbk+W5j5ascLFkZejFUctlmYcHTnJ5yeEy0leRrxzEvxD1NJmq+UUTXvQk1M07SFdc5KXYfOIO9t6wy+IZ8sUG5/gr1aBXRk7XsctwO8iSaKP2z5B/13SY2HU3a8pzToyRERESYHrwBgm4Sv7EhEREUWLZ2SIiIhMFu0drJVlGFQZi2NHhoiIyGwcWjJMWnVkpDujKhZNClpISX7nXkB5t9fghZOk5/J9NGOFiSMvRyuOWiytOHpyUlvcSisneRnxzsnmT0HafEZ2sFshp2jag9E5RdLGUzEneRmCXWxAQnb4rzvBK5Zr898AucOjjB234xbgd1GQhB23Nv13SafkkFYdGSIiomRgxNBSzJdvpwh2ZIiIiMzGoSXDsCNDRERkOht4iwFj8PJrIiIisqy0OiNzylMQsk1aAVKaHBa8T1FWV8h7pBUgpclhavsElxMujlo5wXHU9tGK01+saHKK5rNTK0dPThn+FKTGebivKPCaFXKKpj2o7aMVp79YzEklTr5Ybu/AXAD6vvSyxBtl48TZIgDA555izboYcdwC/C7SimV2G+/J6gt5b1wYMbTEoSkAadaRISIiSgqcI2MYdmSIiIjMJsCAq444xwbgHBkiIiKysLQ6I9PscSh+6lGS3a76ezxjJXMcM2IFhrL9P/d2DY9LHDNzSub2YGYss+MAQKFTbEjd54nvD509ESq3VTznf7R1IADgQHFZ2Fhm5pRqf6dkOm67s82ZI2PE3a95iwJRWnVkiIiIkgLnyBiGQ0tERERkWTwjQ0REZDbBFvtkX57RAcCODBERkelsgviIqQxjqmJ5adWRaekVJ4a19uUFtmktiiRNCHPZz00MK7Z3ANBefEk+8Sw4Vn+LPAXH0oqjFksrjp6cpDh6cpLi6MlJPpkumpzsXeLRndkj/tx/5gJL5SRfjCvc30ltgbFw7UFPTpG08VTMSX7cljrE378qLgKgb7JvzmmxXPc3AwAAn3S6FHH05BTJcRtpTvwugq6comnjXSp34o4bLohnCM6RISIiIstKqzMyREREScGIOTIWGVyaOnVqxO9ZtWoVSkpKdO3LjgwREZHZ0ujy602bNuHHP/4x8vLywu8MYN26dejo6GBHRo00nuruOjfG2u0RP4LcoHFRaWxVPnY73H4KADAgwwsAOOvLVLxHPi4dHEsrjlqscHHksbTi6MlJiqMnJ/lYfric5OPS0eQkzY3Jbhf3Pdl57nO1ak5Gtgc9OUXSxlMxJ/lxW5p3BgBwolB8nul0AgC8bW0h9ZRkt/sXRTuTGxJbb06RfHaR5sTvIujKKZo2fjbTC1Ok2U0jV6xYobtj8uc//zmisjlHhoiIiOLmzTffxKBBg3Tvv2XLFlxwwQXhd/RLqzMyRERESSGNhpa++93vRrT/xIkTI9qfHRkiIiKzpfmCeC0tLWhpaYHP51NsHzduXMRlsSNDREREpmhqakJNTQ0+/vhjCP67XtpsNgiCAJvNBq838jlKadWRkRZNkk9C6/Zkqe4jkS/YJU0OK5Q60f7n8n20YoWLIy9HK45arHBx+stJPtEyXE7yMuKdk83/UoZHsGRO0bQHPTkVyv/zFianSNp4KuYkLyMvU5y46/O/xZaTjXBsHrEuGR6bIna8jlvxd34X6ckp3set12bOZF8jVva1onvuuQejRo3C6tWr4XK5YLPFfgl5WnVkiIiIkkIazZGRO3z4MP7v//2/uOiiiwwrk1ctERERkSluuOEG/P3vfze0TJ6RISIiIlP8/ve/R01NDfbv349LL70UWVnKYcGbbrop4jLZkSEiIjKZIXe/tuDQVENDA9555x1s2bIl5DVO9tVBujOq2oqW0jb5nXsB5d1eAytABq0IKd8nXCytOPJytOKoxYolJ8WKlmFykpcR75wE/0s+u81ffp+lcoqmPejKSWUVXCPaQyrmJC+jyyv+jy/D34yEntAJscEEu1gXn13oN448llnHbX+x+F0UexvvEcxa2deIey1Zz5w5c3DnnXdi4cKFcLlchpSZVh0ZIiKipJBmtyiQnDp1CvPmzTOsEwMk0WTfJUuWwGazYe7cuYFt3d3dqK2tRXFxMQoKCjBt2jQ0NzcnrpJEREQUtalTp+LNN980tMykOCPz7rvv4plnnglZ0W/evHl49dVXsWHDBjidTsyePRtTp07FO++8k6CaEhERGSBNL78eNWoUFixYgLfffhtjx44Nmex7//33R1xmVB2ZNWvW4LbbbsOAAQOiebtCR0cHpk+fjueeew6//vWvA9vb2tqwevVqrFu3Dtdff30g7sUXX4zGxkZMmDAh4lhFWV0h26SFk6Qx1eB9TnkKQt4jLZwkjamq7RNcTrg4auUEx1HbRytOf7GiySmaz06tHD05eXPEceNehzh2fWH+uTvYWiGnaNqD2j5acfqLxZxC47i7xNteZ4s3we73rteB8h3+L9dCcWJNaV67Zl2MOG4BfhdpxTK7jXd5Q+fXxIMhC+JZsCP0+9//HgUFBdi5cyd27typeM1ms0XVkYlqaOmXv/wlSktLMXPmTOzatSuaIgJqa2sxZcoUVFVVKbY3NTWhr69PsX306NEoLy9HQ0NDTDGJiIjIfEeOHNF8HD58OKoyo+rInDhxAs8//zy+/vprXHfddRg9ejSefPJJuN3uiMpZv3499u7di7q6upDX3G43srOzUVRUpNjucrn6jdPT04P29nbFg4iIKKkIBjwIQJQdGbvdjltvvRUvv/wyjh8/jnvvvRcvvPACysvLcdNNN+Hll18OuaNlsOPHj+OBBx7ACy+8gNzc3Kgqr6aurg5OpzPwGDp0qGFlExERGSKNOjLz589HZ2en7v0XLFiA06dP694/5quWXC4XJk6ciMrKSmRkZGDfvn2oqanBt771LezYsUPzfU1NTWhpacGVV14Ju90Ou92OnTt3YsWKFbDb7XC5XOjt7UVra6vifc3NzSgtLdUsd8GCBWhraws8jh8/HmuKRERExhLOzZOJ9mGVzsxvf/tbnD17Vvf+9fX1If/29yfqq5aam5vxX//1X1izZg0OHz6MW265BZs3b0ZVVRU6OzuxePFi1NTU4OjRo6rvv+GGG7Bv3z7FtrvvvhujR4/Gww8/jKFDhyIrKwvbt2/HtGnTAAAHDx7EsWPHUFlZqVmvnJwc5OTkqL5Wkt2u+KkrT49D9fdwYomVzHHMiOXJsyl+Xlp4Ii5xzMwpmduDmbHMjgMA7nbx94JT/Z8llusZJH41Oga2AgBG5Wsv+5CInFLt75RMx223p6/f1ylygiBg1KhRuu90HcnZGyDKjsyNN96I119/HaNGjcK9996Lu+66C4MGDQq8np+fj5///Of493//d80yCgsLcemllyq25efno7i4OLB95syZmD9/PgYNGgSHw4E5c+agsrIyqiuWiIiIkoYhK/taY2XgNWvWRPyeSBbMi6ojU1JSgp07d/Z7ZuT888/HkSNHoik+YNmyZcjIyMC0adPQ09OD6upqPP300zGVSURElBTS5PLrmpqauJYf1RyZ1atX4+zZs/hf/+t/4ac//SnuuecexQMQrwcfNmxYROXu2LEDy5cvDzzPzc1FfX09Tp8+jc7OTrz00kv9zo8hIiKygljnx9iEyM/H1NXV4eqrr0ZhYSFKSkpwyy234ODBg4p99Kyof+zYMUyZMgUDBgxASUkJHnzwQXg85qy/oyaqjszixYsxadIkbN++HV9//TW++eYbxYOIiIiSy86dO1FbW4vGxkZs3boVfX19mDRpkmJOyrx58/DKK69gw4YN2LlzJ06ePImpU6cGXvd6vZgyZQp6e3uxa9cuPP/881i7di0WLVqUiJQARDm0tHLlSqxduxYzZswwuj5x5bKLE8OK7R2BbVqrO0oTwlp6z00Ma+3LA6C9iqR84llwrP5WqwyOpRVHLZZWHD05SXH05CTF0ZOTfDJdNDlJ1fL5FxG9Mu9zS+UkX1U03N9JbaXUcO1BT06RtPFUzEl+3J5pE2Od97X+/zF2F4n/1x1WJP7HbEzuSUUcPTlFctxGmhO/i6Arp2jaeE+vSZN9E3DTyNdee03xfO3atSgpKUFTUxP++Z//WdeK+n/9619x4MABbNu2DS6XC5dffjl+9atf4eGHH8Zjjz2G7OzQO7DHW1RnZHp7e/Gd73zH6LoQERGlBSOGlgDxDEnwIrA9PT266tDmv12HdLGOnhX1GxoaMHbsWMVk3OrqarS3t+Ojjz4y4qOJWFQdmZ/+9KdYt26d0XUhIiKiCDQ2NioWgXU6naqr5Qfz+XyYO3curr322sCVwnpW1He73SFXFEnPI13d3yhRDS11d3fj2WefxbZt2zBu3LiQu1c+9dRThlSOiIgoJRm0oN2ECROwZcsWxTattdTkamtrsX//frz99tuxVyLBourIfPjhh7j88ssBAPv371e8pnfBm0SQxlOH208Ftg3I8AIAzvoyFftKY6vysVt3l7it2yN+bLl25Zi7fFw6OJZWHLVY4eLIY2nF0ZOTFEdPTvKx/HA5ycelo8lJmhvjyRd/XpjVavmcjGwPenKKpI2nYk7y4xad/rkQ34gTGvUsi9fnb3sXDGgNia03p0g+u0hz4ncRdOUUTRvv7Qv9XOLCoDkymZmZcDj0LygIALNnz8bmzZvx1ltvYciQIYHtpaWlgRX15Wdl5Cvql5aWYs+ePYrypKuaEnVVcVQdmTfffNPoehAREaUN+TyXqMuIcH9BEDBnzhxs3LgRO3bswIgRIxSvjx8/PuyK+pWVlXj88cfR0tKCkpISAMDWrVvhcDgwZsyY2BKKUtS3KCAiIiLrqK2txbp16/Dyyy+jsLAwMKfF6XQiLy8PTqcz7Ir6kyZNwpgxYzBjxgwsXboUbrcbjzzyCGpra3UNacUDOzJERERpYOXKlQCA6667TrF9zZo1+MlPfgIg/Ir6mZmZ2Lx5M2bNmoXKykrk5+ejpqYGixcvNiuNEOzIEBERmS0B68gIQvg3SCvq19fXa+4zbNgw/OUvf4kseBylVUdGWjRJPgmtUBpk9G+TL+IEKBfskiaHdXuyNPfRihUujrwcrThqsTTj6MhJPtEyXE7yMuKdk+CfaydttlpO0bQHPTlF0x7SNSd5GTaPuLOtN/yCeLZMsfEJ/moV2JWx43XcAvwukiT6uO0T9N8lPRZGzJEhUVTryBARERElg7Q6I0NERJQUEjC0lKrYkSEiIjKbQQviUZp1ZKQbiikWTQpaSEl+wztAeZO04IWTpOfyfTRjhYkjL0crjlosrTh6clJb3EorJ3kZ8c7J5k9B2nxGdrBbIado2oPROUXSxlMxJ3kZgl1sQEJ2+K87wSuWa/PfN7DDo4wdt+MW4HdRkIQdtzb9Nxel5JBWHRkiIqJkwMm+xmFHhoiIyGycI2MYdmSIiIhMlohbFKQqXn5NRERElpVWZ2ROeQpCtkkLJ0mTw4L3KcrqCnmPtHCSNDlMbZ/gcsLFUSsnOI7aPlpx+osVTU7RfHZq5ejJKcOfgtQ4D/cVBV6zQk7RtAe1fbTi9BeLOanEyRfL7R2YC0Dfl16WeKNsnDhbBAD43FOsWRcjjluA30Vascxu4z1ZfSHvjQsOLRkmrToyRERESYGXXxuGHRkiIiKTGXLVEjtCADhHhoiIiCyMZ2SIiIgSgWdUDJFWHZlmj0PxU4+S7HbV3+MZK5njmBErMCfP/3Nv1/C4xDEzp2RuD2bGMjsOABQ6xYbUfZ74/tBpoKFyW8V/YY62DgQAHCguCxvLzJxS7e+UTMdtd7aFJvsSAA4tERERkYWl1RkZIiKiZGDIgng8owOAHRkiIiLzcWjJMGnVkWnpFcdTW/vyAtu0FkWSxlFd9nPjqcX2DgDaiy/Jx2uDY/W3yFNwLK04arG04ujJSYqjJycpjp6c5GPQ0eRk7xKP7swe8ef+MxdYKif5Ylzh/k5qC4yFaw96coqkjadiTvLjttQh/v5VcREAfXNkck6L5bq/GQAA+KTTpYijJ6dIjttIc+J3EXTlFE0b71K5E3c88PJr43CODBEREVlWWp2RISIiSgocWjIMOzJERERmY0fGMOzIEBERmczmf8RaBqVZR0aaGObuOjdZrNsjfgS5QRO8pEli8klow+2nAAADMrwAgLO+TMV75BPsgmNpxVGLFS6OPJZWHD05SXH05CSflBguJ/kEu2hykib5ZreL+57sPPe5WjUnI9uDnpwiaeOpmJP8uC3NOwMAOFEoPs90OgEA3ra2kHpKstv9i6KdyQ2JrTenSD67SHPidxF05RRNGz+b6QVZS1p1ZIiIiJKCEUNLHJoCwI4MERGR6Qy5/JoAJPjy65UrV2LcuHFwOBxwOByorKzEli1bAq93d3ejtrYWxcXFKCgowLRp09Dc3JzAGhMREVEySWhHZsiQIViyZAmamprw3nvv4frrr8fNN9+Mjz76CAAwb948vPLKK9iwYQN27tyJkydPYurUqYmsMhERUewEgx6U2KGlG2+8UfH88ccfx8qVK9HY2IghQ4Zg9erVWLduHa6//noAwJo1a3DxxRejsbEREyZMiDietPqjfBJatydLdR+JfOVRaXJYoTRV3P9cvo9WrHBx5OVoxVGLFS5OfznJJ1qGy0leRrxzsvlfyvAIlswpmvagJ6dC+SUKYXKKpI2nYk7yMvIyxYm7Pv9bbDnZCMfmEeuS4bEpYsfruBV/53eRnpzifdx6bSZN9mVHxDBJs7Kv1+vF+vXr0dnZicrKSjQ1NaGvrw9VVVWBfUaPHo3y8nI0NDRoltPT04P29nbFg4iIKJlIc2RiebAjJEp4R2bfvn0oKChATk4Ofvazn2Hjxo0YM2YM3G43srOzUVRUpNjf5XLB7XZrlldXVwen0xl4DB06NM4ZEBERUaIkvCPz7W9/Gx988AF2796NWbNmoaamBgcOHIi6vAULFqCtrS3wOH78uIG1JSIiMgDnxxgm4ZdfZ2dn46KLLgIAjB8/Hu+++y5++9vf4rbbbkNvby9aW1sVZ2Wam5tRWlqqWV5OTg5ycnLUY/nvjKq2EJS0TX7nXkB5t9fAwklBCynJ9wkXSyuOvBytOGqxYslJsRBUmJzkZcQ7J8H/ks9u85ffZ6mcomkPunJSWTzOiPaQijnJy+jyivMnMvzNSOgJnUcSTLCLdfHZhX7jyGOZddz2F4vfRbG38R7BnDkyRlx+zZV9RQk/IxPM5/Ohp6cH48ePR1ZWFrZv3x547eDBgzh27BgqKysTWEMiIqIY8aolwyT0jMyCBQswefJklJeX48yZM1i3bh127NiB119/HU6nEzNnzsT8+fMxaNAgOBwOzJkzB5WVlVFdsURERESpJ6EdmZaWFtx111348ssv4XQ6MW7cOLz++uv4/ve/DwBYtmwZMjIyMG3aNPT09KC6uhpPP/10IqtMREQUM67sa5yEdmRWr17d7+u5ubmor69HfX29STUiIiIyAe+1ZJiET/Y1U1FWV8g2aeEkaXJY8D6nPAUh75EWTpImh6ntE1xOuDhq5QTHUdtHK05/saLJKZrPTq0cPTl5c8QpbL0OcRLehfnn1gKyQk7RtAe1fbTi9BeLOYXGcXeJt73OFm+C3e9drwPlO/wLrBWKM4RL89o162LEcQvwu0grltltvMsbOlGYkltadWSIiIiSBs+oGIIdGSIiIpMZcvk1O0IA2JEhIiIyHy+fNkxadWRKstsVP/Vo9jhUf49nrGSOY0YsT55N8fPSwhNxiWNmTsncHsyMZXYcAHC3i78XnPLpfn/PIPGr0TGwFQAwKr85bCwzc0q1v1MyHbfdnr5+X6fkk1YdGSIiomRgEwTYBF62ZAR2ZIiIiMzGy68Nw44MERGR2XivJcMk3b2WiIiIiPRKqzMyLrs4MazY3hHYprUokjQhrKX33MSw1r48ANqLL8knngXH6m+Rp+BYWnHUYmnF0ZOTFEdPTlIcPTnJJ9NFk5NULZ9/7a0r8z63VE7yxbjC/Z3UFhgL1x705BRJG0/FnOTH7Zk2MdZ5X+tf6Ky7SPy/7rCibwAAY3JPKuLoySmS4zbSnPhdBF05RdPGe3pNmuzLoSXDpFVHhoiIKBnwXkvG4dASERERWRbPyBAREZmNQ0uGYUeGiIjIZBxaMk5adWSkiWHD7acC2wZkeAEAZ32Zin2lSWLySWjuLnFbt0f82HLtysmD8gl2wbG04qjFChdHHksrjp6cpDh6cpJPSgyXk3yCXTQ5SZN8PfnizwuzWi2fk5HtQU9OkbTxVMxJftyi0z+p85tOAICe9X37/G3vggGtIbH15hTJZxdpTvwugq6comnjvX2hn0tc8BYFhuEcGSIiIrIsdmSIiIhMJg0txfqIxFtvvYUbb7wRZWVlsNls2LRpk+J1QRCwaNEiDB48GHl5eaiqqsKhQ4cU+5w+fRrTp0+Hw+FAUVERZs6ciY6ODiQSOzJERERmE4TYHxHq7OzEZZddhvr6etXXly5dihUrVmDVqlXYvXs38vPzUV1dje7u7sA+06dPx0cffYStW7di8+bNeOutt3DfffdF/TEYIa3myEiLJsnHbgulNZ792+SLOAHKBbukMdVuT5bmPlqxwsWRl6MVRy2WZhwdOcnnJ4TLSV5GvHMS/EPU0mar5RRNe9CTUzTtIV1zkpdh84g723rDL4hnyxQbn+CvVoFdGTtexy3A7yJJoo/bPkH/XdJjYchk3wjfP3nyZEyePFm9KEHA8uXL8cgjj+Dmm28GAPzhD3+Ay+XCpk2bcPvtt+Pjjz/Ga6+9hnfffRdXXXUVAOB3v/sd/sf/+B/43//7f6OsrCymdKLFMzJEREQW5fV60d7ernj09PREXM6RI0fgdrtRVVUV2OZ0OlFRUYGGhgYAQENDA4qKigKdGACoqqpCRkYGdu/eHXsyUWJHhoiIyGyCAQ8AjY2NcDqdikddXV3E1XG73QAAl8ul2O5yuQKvud1ulJSUKF632+0YNGhQYJ9ESKuhJSIiomRg84mPmMoQgAkTJmDLli2K7Tk5ObEVbDE8I0NERGRRmZmZcDgcikc0HZnS0lIAQHNzs2J7c3Nz4LXS0lK0tLQoXvd4PDh9+nRgn0RIqzMy0p1RFYsmBS2kJL9zL6C822vwwknSc/k+mrHCxJGXoxVHLZZWHD05qS1upZWTvIx452TzpyBtPiOb0GaFnKJpD0bnFEkbT8Wc5GUIdrEBCdnhv+4Er1iuzX8D5A6PMnbcjluA30VBEnbc2vTfJT1mSbQg3ogRI1BaWort27fj8ssvBwC0t7dj9+7dmDVrFgCgsrISra2taGpqwvjx4wEAb7zxBnw+HyoqKhJV9fTqyBARESWDRFy11NHRgU8//TTw/MiRI/jggw8waNAglJeXY+7cufj1r3+NkSNHYsSIEVi4cCHKyspwyy23AAAuvvhi/OAHP8C9996LVatWoa+vD7Nnz8btt9+esCuWAHZkiIiIzBflWjCxeO+99/C9730v8Hz+/PkAgJqaGqxduxYPPfQQOjs7cd9996G1tRUTJ07Ea6+9htzc3MB7XnjhBcyePRs33HADMjIyMG3aNKxYscLUPIKxI0NERJQGrrvuOgj9dJ5sNhsWL16MxYsXa+4zaNAgrFu3Lh7Vi1padWROeQpCtkkLJ0ljqsH7FGV1hbxHWjhJGlNV2ye4nHBx1MoJjqO2j1ac/mJFk1M0n51aOXpyyvCnIDXOw31FgdeskFM07UFtH604/cViTipx8sVyeweK/6vU86WXJd5fEifOFgEAPvcUa9bFiOMW4HeRViyz23hPVl/Ie+PBiKEl3j1blFYdGSIioqTAu18bhh0ZIiIikyVism+q4joyREREZFk8I0NERGQ2Q65a4ikZIM06Ms0eh+KnHiXZ7aq/xzNWMscxI1ZgTp7/596u4XGJY2ZOydwezIxldhwAKHSKDan7PPH9odNAQ+W2iv9AHG0dCAA4UKy9RkYickq1v1MyHbfd2daZ7EsiDi0RERGRZaXVGRkiIqKkYMRVSzyjAyDBZ2Tq6upw9dVXo7CwECUlJbjllltw8OBBxT7d3d2ora1FcXExCgoKMG3atJCbWhEREVmJNLQU0yPRSSSJhHZkdu7cidraWjQ2NmLr1q3o6+vDpEmT0NnZGdhn3rx5eOWVV7Bhwwbs3LkTJ0+exNSpUxNYayIiohj5hNgfPCMDIMFDS6+99pri+dq1a1FSUoKmpib88z//M9ra2rB69WqsW7cO119/PQBgzZo1uPjii9HY2IgJEyZEFK+lV5wY1tqXF9imtbqjNCHMZT83MazY3gFAexVJ+cSz4Fj9rVYZHEsrjlosrTh6cpLi6MlJiqMnJ/lkumhysneJR2dmj/hz/5kLLJWTfFXRcH8ntZVSw7UHPTlF0sZTMSf5cVvqEH//qrgIgL7JvjmnxXLd3wwAAHzS6VLE0ZNTJMdtpDnxuwi6coqmjXep3ImbkltSTfZta2sDIN7LAQCamprQ19eHqqqqwD6jR49GeXk5GhoaVMvo6elBe3u74kFERJRUBAMeBCCJOjI+nw9z587Ftddei0svvRQA4Ha7kZ2djaKiIsW+LpcLbrdbtZy6ujo4nc7AY+jQofGuOhERUUSMmCPDzowoaToytbW12L9/P9avXx9TOQsWLEBbW1vgcfz4cYNqSEREZBAB5xbFi/bBngyAJLn8evbs2di8eTPeeustDBkyJLC9tLQUvb29aG1tVZyVaW5uRmlpqWpZOTk5yMnJUX1NGk91d50bY+32iB9BbtC4qDS2Kh+7HW4/BQAYkOEFAJz1ZSreIx+XDo6lFUctVrg48lhacfTkJMXRk5N8LD9cTvJx6WhykubGZLeL+57sPPe5WjUnI9uDnpwiaeOpmJP8uC3NOwMAOFEoPs90OgEAXv9Qtprsdv+iaGdyQ2LrzSmSzy7SnPhdBF05RdPGz2Z6QdaS0DMygiBg9uzZ2LhxI9544w2MGDFC8fr48eORlZWF7du3B7YdPHgQx44dQ2VlpdnVJSIiMoQhl1/zhAyABJ+Rqa2txbp16/Dyyy+jsLAwMO/F6XQiLy8PTqcTM2fOxPz58zFo0CA4HA7MmTMHlZWVEV+xRERElDQ4MmSYhHZkVq5cCQC47rrrFNvXrFmDn/zkJwCAZcuWISMjA9OmTUNPTw+qq6vx9NNPm1xTIiIiSkYJ7cgIOu78mZubi/r6etTX15tQIyIiovizQYAt1rtf84wOgCSZ7GsWadEk+SS0bk+W6j4S+YJd0uSwQmldaP9z+T5ascLFkZejFUctVrg4/eUkn2gZLid5GfHOyeZ/KcMjWDKnaNqDnpwK5euRh8kpkjaeijnJy8jLFCfu+vxvseVkIxybR6xLhsemiB2v41b8nd9FenKK93HrtZk02dfnf1DM0qojQ0RElAxsQuxnZGI+o5MikmYdGSIiIqJI8YwMERGR2XjVkmHYkSEiIjJbYHXeWMowpipWl1YdGenOqGorWkrb5HfuBZR3ew2sABm0IqR8n3CxtOLIy9GKoxYrlpwUK1qGyUleRrxzEvwv+ew2f/l9lsopmvagKyeVVXCNaA+pmJO8jC6vOBE0w9+MhJ7QCbHBBLtYF59d6DeOPJZZx21/sfhdFHsb7xHMmezLBe2MwzkyREREZFlpdUaGiIgoKXBoyTDsyBAREZnM5hMfFLu06sgUZXWFbJMWTpLGVIP3OeUpCHmPtHCSNKaqtk9wOeHiqJUTHEdtH604/cWKJqdoPju1cvTk5M0R58b0OsSx6wvzz93B1go5RdMe1PbRitNfLOYUGsfdJd72Olu8CXa/d70OlO/wL7BWKE6sKc1r16yLEcctwO8irVhmt/Eub+j8GkpuadWRISIiSgqGDC1xbAlgR4aIiMh8XEfGMOzIEBERmc2QWxQYVBeL4+XXREREZFlpdUamJLtd8VOPZo9D9fd4xkrmOGbE8uTZFD8vLTwRlzhm5pTM7cHMWGbHAQB3u/h7wSn9l4j0DBK/Gh0DWwEAo/Kbw8YyM6dU+zsl03Hb7enr93XDGDFHhmNTANKsI0NERJQUfP5HLNiPAcCODBERkelsBsyRIRHnyBAREZFlpdUZGZddHE8ttncEtmktiiSNo7b0nhtPbe3LA6C9+JJ8vDY4Vn+LPAXH0oqjFksrjp6cpDh6cpLi6MlJPgYdTU5StXz+tbeuzPvcUjnJF+MK93dSW2AsXHvQk1MkbTwVc5Ift2faxFjnfa1/obPuInF+1rCibwAAY3JPKuLoySmS4zbSnPhdBF05RdPGe3pNmiMDriNjlLTqyBARESUFQyb7EsChJSIiIrIwnpEhIiIymw+ALcYyeEIHADsyREREpjPiqiWu7CtKq46MNDFsuP1UYNuADC8A4KwvU7GvNElMPgnN3SVu6/aIH1uuXTl5UD7BLjiWVhy1WOHiyGNpxdGTkxRHT07ySYnhcpJPsIsmJ2mSrydf/HlhVqvlczKyPejJKZI2noo5yY9bdPondX7TCUDf0h19/rZ3wYDWkNh6c4rks4s0J34XQVdO0bTx3r7QzyUuuCCeYThHhoiIiCwrrc7IEBERJQVetWQYdmSIiIjMZkRHhh0hAOzIEBERmY9XLRkmrToy0uqP8klohVJD8m+Tr0YJKFcelSaHdXuyNPfRihUujrwcrThqsTTj6MhJPtEyXE7yMuKdk+CfaydttlpO0bQHPTlF0x7SNSd5GTaPuLOtN/zKvrZMsfEJ/moV2JWx43XcAvwukiT6uO0TYr2TI5ktrToyREREyYA3jTQOOzJERERm4xwZw/DyayIiIrKstDojI90ZVbFoUtBCSvI79wLKu70GL5wkPZfvoxkrTBx5OVpx1GJpxdGTk9riVlo5ycuId042fwrS5jOy/3RYIado2oPROUXSxlMxJ3kZgl1sQEJ2+K87wSuWa/PfALnDo4wdt+MW4HdRkIQdtzb9d0mPiQ9cmtcgadWRISIiSgqGDC0ZUxWr49ASERGR2aSOTCyPKHsy9fX1GD58OHJzc1FRUYE9e/YYm5vJEtqReeutt3DjjTeirKwMNpsNmzZtUrwuCAIWLVqEwYMHIy8vD1VVVTh06FBiKktERGRxf/rTnzB//nw8+uij2Lt3Ly677DJUV1ejpaUl0VWLWkI7Mp2dnbjssstQX1+v+vrSpUuxYsUKrFq1Crt370Z+fj6qq6vR3d1tck2JiIgMZMQZmShOyDz11FO49957cffdd2PMmDFYtWoVBgwYgP/8z/80PkeTJHSOzOTJkzF58mTV1wRBwPLly/HII4/g5ptvBgD84Q9/gMvlwqZNm3D77bdHHO+UpyBkm7RwkjQ5LHifoqyukPdICydJk8PU9gkuJ1wctXKC46jtoxWnv1jR5BTNZ6dWjp6cMvwpSI3zcF9R4DUr5BRNe1DbRytOf7GYk0qcfLHc3oG5APR96WWJN8rGibNFAIDPPcWadTHiuAX4XaQVy+w23pPVF/LeuBAEwBf73a+9Xi/a29sVW3NycpCTkxOyd29vL5qamrBgwYLAtoyMDFRVVaGhoSHGuiRO0s6ROXLkCNxuN6qqqgLbnE4nKioq+v3Ae3p60N7erngQERElFcFnwANobGyE0+lUPOrq6lRDfv311/B6vXC5XIrtLpcLbrfbjKzjImmvWpI+1Eg/8Lq6Ovzbv/1bXOtGRESUDCZMmIAtW7YotqmdjUllSXtGJloLFixAW1tb4HH8+PFEV4mIiEjJoKuWMjMz4XA4FA+tjsx5552HzMxMNDc3K7Y3NzejtLTUhKTjI2k7MtKHGukHnpOTE/JHJSIiSio+IfZHhFNssrOzMX78eGzfvv1cNXw+bN++HZWVlQYnaB6bICTHzRpsNhs2btyIW265BYA42besrAy/+MUv8POf/xwA0N7ejpKSEqxdu1b3ZN/29nY4nU4s2l2F3ILQO7gSERFJujv6sLhiG9ra2uLyH2Hp36QbLvgZsjJiGwI68M0O3DZrEpYuXar7PX/6059QU1ODZ555Btdccw2WL1+OF198Ef/4xz9CpnJYRULnyHR0dODTTz8NPD9y5Ag++OADDBo0COXl5Zg7dy5+/etfY+TIkRgxYgQWLlyIsrKyQGeHiIjIkhJ008jbbrsNX331FRYtWgS3243LL78cr732mmU7MUCCOzLvvfcevve97wWez58/HwBQU1ODtWvX4qGHHkJnZyfuu+8+tLa2YuLEiXjttdeQm5ubqCoTERHFzoiOTJRmz56N2bNnJyR2PCS0I3Pdddehv5Etm82GxYsXY/HixSbWioiIKM4M6cgkxcyQhEvay6/joaVXHO9s7csLbNNaFKkkW1x/xmU/tw5Nsb0DgPbiS82ec+OpwbH6W+QpOJZWHLVYWnH05CTF0ZOTFEdPTlKcdM1JvhhXuJzUFhgL1x705BRJG0/FnKx23KZiTlY9brtU7sRNyS2tOjJERERJwecD4IutDJ6QAcCODBERkfkSNNk3FbEjQ0REZLYETvZNNUm7IB4RERFROGl1RkaaGObuOjdZrNsjfgS5QRO8pEli8klow+2nAAADMrwAgLO+TMV75BPsgmNpxVGLFS6OPJZWHD05SXH05CSflBguJ/kEO+bUf07RtAc9OUXSxlMxJ6sdt6mYk1WP27OZXpjCJ95iICY8owMgzToyREREyUAQBAhCjJN9OdsXAIeWiIiIyMJ4RoaIiMhsgv/GjzGVYUxVrC6tOjLSoknysdtuT5bqPhL5gl3SmGqhzb/B/1y+j1ascHHk5WjFUYsVLk5/OcnnJ4TLSV4Gc+o/p2jag56cAnF05BRJG0/FnKx23KZiTlY9br02k+bICEbMkTGkJpaXVh0ZIiKipODzAbZYF8RjTwbgHBkiIiKyMJ6RISIiMpsRQ0scWwLAjgwREZHpBJ8PQoxDSwKHlgCkWUdGujOq2kJQ0jb5nXsB5d1eAwsnBS2kJN8nXCytOPJytOKoxYolJ8VCUGFykpfBnPrPKZr2oCsnlcXjjPjsUjGnVD5uUzGnZDpuewQLTfYlAJwjQ0RERBaWVmdkiIiIkoJPAGy8/NoI7MgQERGZTfAB4OXXRuDQEhEREVlWWp2RKcrqCtkmrQApTQ4L3ueUpyDkPdIKkNLkMLV9gssJF0etnOA4avtoxekvVjQ5RfPZqZWTDjlF0x7U9tGK018s5qQeRx4rGY9bPbGslpNVj9sub+hE4XgQfAKEWIeWOLYEIM06MkREREmBQ0uGYUeGiIjIbAackWE/RsQ5MkRERGRZNiHFlwZsb2+H0+nEot1VyC0IvRsrERGRpLujD4srtqGtrQ0Oh8Pw8qV/k76Lm5BlC13AMBL/EN7HHQ9OxdKlSw2qnTXxjAwREZFJ8vLykINcnEZLTOUIgoDTaMHIkSMNqpl1sSNDRERkkqysLCxfuQyf4SP4hOgn+36JoxDgQ01NjYG1syZ2ZIiIiEw0c+ZMCPDBjWNRvd8neHEYB/Dsf61CdnZsw1OpgB0ZIiIiE2VlZeG5Pz6DwzgAXxQ3qTyBI8iEHXfccUccamc9aXX5tcveDgAotncEtmktitTsESd5tfSem+zV2pcHQHvxpZLsds1Y/S3yFBxLK45aLK04enKS4ujJSYqjJycpTrrmJF+MK1xOaguMhWsPenKKpI2nYk5WO25TMSerHrc9vX0ww+23347//85/wRc4gnJcpPt9XsGDI/gYf9q4HpmZmeHfkAZ4RoaIiMhkmZmZWPfyf+FzfAyPoH814WP4FLkYgJtvvjmOtbMWdmSIiIgS4MYbb0Qe8nEch3Tt3yf04igO4qVtG2Cz2eJcO+tgR4aIiCgBbDYbNr7xZxzFJ+gTesPufxQH4cBA3HDDDSbUzjrSao6MNJ463H4qsG1AhjjR6qxPOdYoja3Kx27dXeK2bo/4seXalacD5ePSwbG04qjFChdHHksrjp6cpDh6cpKP5YfLST4uzZz6zyma9qAnp0jaeCrmZLXjNhVzsupx29tn7ryT733ve3BgII7iIC7CWM39eoRuHMOn2NX4jom1swaekSEiIkqgLbtfwTF8ih6hW3Ofz/EPFMOFiooKE2tmDezIEBERJdA111yDYrhwBB+rvt4ldOIEDuONfa+bXDNrsERHpr6+HsOHD0dubi4qKiqwZ8+eRFeJiIjIMG/u/ytO4gi6hM6Q1w7jAEowBJdeemkCapb8kr4j86c//Qnz58/Ho48+ir179+Kyyy5DdXU1Wlpiu08FERFRsrjkkkvgwlAcxgHF9g6hHc04jnc+ezNBNUt+ST/Z96mnnsK9996Lu+++GwCwatUqvPrqq/jP//xP/PKXv4yoLGnRJPkktELpCjb/NvkiToBywS5pcli3J0tzH61Y4eLIy9GKoxZLM46OnOQTLcPlJC+DOfWfUzTtQU9O0Xx26ZqT1Y7bVMzJqsdtXwz3P4rVO4ffxEUXXoRhwigU2JwAgMP4CGUYjgsvvDBh9Up2SX1Gpre3F01NTaiqqgpsy8jIQFVVFRoaGlTf09PTg/b2dsWDiIgo2Y0YMQJlGIHP/Gdl2oXT+Bpu7Dn5doJrltySuiPz9ddfw+v1wuVyKba7XC643W7V99TV1cHpdAYeQ4cONaOqREREMdtz8m84BTfahNP4FB9hKL6FwYMHJ7paSS2pOzLRWLBgAdra2gKP48ePJ7pKREREugwePBhDcRH2Yw/acApNX+9KdJWSXlLPkTnvvPOQmZmJ5uZmxfbm5maUlpaqvicnJwc5OTmB54IgAAB6OjzozBXHTzuyZGOg/rHVDnE3dPaJ+3R3izcO6+s8N9bq6ewRf3rE93vs/n3gv2GZ7dzNxjo9QbE04qjF0oqjFksrjp6cpDh6cpLi6MlJipOuOUlx9OQkxdGTkxRHT06RtPFUzMlyx20q5mTR41baT/q3IxH2ntqF0uLBGIZvo7i4OGH1sAqbkMi/lg4VFRW45ppr8Lvf/Q4A4PP5UF5ejtmzZ+ua7PvFF19weImIiCJy/PhxDBkyJNHVIB2S+owMAMyfPx81NTW46qqrcM0112D58uXo7OwMXMUUTllZGQ4cOIAxY8bg+PHjcDgc4d+UYO3t7Rg6dCjrGyesb3xZrb6A9erM+saPIAg4c+YMysrKEl0V0inpOzK33XYbvvrqKyxatAhutxuXX345XnvttZAJwFoyMjJwwQUXAAAcDkfSH0RyrG98sb7xZbX6AtarM+sbH06nM9FVoAgkfUcGAGbPno3Zs2cnuhpERESUZFLuqiUiIiJKH2nRkcnJycGjjz6quJopmbG+8cX6xpfV6gtYr86sL9E5SX/VEhEREZGWtDgjQ0RERKmJHRkiIiKyLHZkiIiIyLLYkSEiIiLLSouOTH19PYYPH47c3FxUVFRgz549ia4S6urqcPXVV6OwsBAlJSW45ZZbcPDgQcU+3d3dqK2tRXFxMQoKCjBt2rSQ+04lypIlS2Cz2TB37tzAtmSr74kTJ3DnnXeiuLgYeXl5GDt2LN57773A64IgYNGiRRg8eDDy8vJQVVWFQ4cOJay+Xq8XCxcuxIgRI5CXl4dvfetb+NWvfqW450si6/zWW2/hxhtvRFlZGWw2GzZt2qR4XU/dTp8+jenTp8PhcKCoqAgzZ85ER0eH6fXt6+vDww8/jLFjxyI/Px9lZWW46667cPLkyaSsb7Cf/exnsNlsWL58eVLX9+OPP8ZNN90Ep9OJ/Px8XH311Th27Fjg9WT7ziBrSvmOzJ/+9CfMnz8fjz76KPbu3YvLLrsM1dXVaGlpSWi9du7cidraWjQ2NmLr1q3o6+vDpEmT0NnZGdhn3rx5eOWVV7Bhwwbs3LkTJ0+exNSpUxNYa9G7776LZ555BuPGjVNsT6b6fvPNN7j22muRlZWFLVu24MCBA/jNb36DgQMHBvZZunQpVqxYgVWrVmH37t3Iz89HdXU1uru7E1LnJ598EitXrsR//Md/4OOPP8aTTz6JpUuXBu4zlug6d3Z24rLLLkN9fb3q63rqNn36dHz00UfYunUrNm/ejLfeegv33Xef6fU9e/Ys9u7di4ULF2Lv3r146aWXcPDgQdx0002K/ZKlvnIbN25EY2Oj6hL6yVTfzz77DBMnTsTo0aOxY8cOfPjhh1i4cCFyc3MD+yTTdwZZmJDirrnmGqG2tjbw3Ov1CmVlZUJdXV0CaxWqpaVFACDs3LlTEARBaG1tFbKysoQNGzYE9vn4448FAEJDQ0OiqimcOXNGGDlypLB161bhu9/9rvDAAw8IgpB89X344YeFiRMnar7u8/mE0tJS4d///d8D21pbW4WcnBzh//yf/2NGFUNMmTJFuOeeexTbpk6dKkyfPl0QhOSqMwBh48aNged66nbgwAEBgPDuu+8G9tmyZYtgs9mEEydOmFpfNXv27BEACEePHhUEITnr+8UXXwgXXHCBsH//fmHYsGHCsmXLAq8lW31vu+024c4779R8T7J9Z5B1pfQZmd7eXjQ1NaGqqiqwLSMjA1VVVWhoaEhgzUK1tbUBAAYNGgQAaGpqQl9fn6Luo0ePRnl5eULrXltbiylTpijqBSRfff/7v/8bV111FX70ox+hpKQEV1xxBZ577rnA60eOHIHb7VbU1+l0oqKiImGf73e+8x1s374dn3zyCQDg73//O95++21Mnjw5aess0VO3hoYGFBUV4aqrrgrsU1VVhYyMDOzevdv0Ogdra2uDzWZDUVERgOSrr8/nw4wZM/Dggw/ikksuCXk9merr8/nw6quvYtSoUaiurkZJSQkqKioUw0/J9p1B1pXSHZmvv/4aXq835AaTLpcLbrc7QbUK5fP5MHfuXFx77bW49NJLAQButxvZ2dmBL1VJIuu+fv167N27F3V1dSGvJVt9Dx8+jJUrV2LkyJF4/fXXMWvWLNx///14/vnnA/WV6pcM9QWAX/7yl7j99tsxevRoZGVl4YorrsDcuXMxffp0AMlZZ4meurndbpSUlChet9vtGDRoUMLr393djYcffhh33HFH4KaGyVbfJ598Ena7Hffff7/q68lU35aWFnR0dGDJkiX4wQ9+gL/+9a+49dZbMXXqVOzcuTNQ32T6ziDrssRNI1NdbW0t9u/fj7fffjvRVdF0/PhxPPDAA9i6datijDtZ+Xw+XHXVVXjiiScAAFdccQX279+PVatWoaamJsG1U/fiiy/ihRdewLp163DJJZfggw8+wNy5c1FWVpa0dU4FfX19+PGPfwxBELBy5cpEV0dVU1MTfvvb32Lv3r2w2WyJrk5YPp8PAHDzzTdj3rx5AIDLL78cu3btwqpVq/Dd7343kdWjFJPSZ2TOO+88ZGZmhsyCb25uRmlpaYJqpTR79mxs3rwZb775JoYMGRLYXlpait7eXrS2tir2T1Tdm5qa0NLSgiuvvBJ2ux12ux07d+7EihUrYLfb4XK5kqq+gwcPxpgxYxTbLr744sAVE1KdkqltPPjgg4GzMmPHjsWMGTMwb968wBmwZKyzRE/dSktLQybZezwenD59OmH1lzoxR48exdatWwNnY4Dkqu/f/vY3tLS0oLy8PHD8HT16FD//+c8xfPjwpKvveeedB7vdHvYYTKbvDLKulO7IZGdnY/z48di+fXtgm8/nw/bt21FZWZnAmomXqs6ePRsbN27EG2+8gREjRiheHz9+PLKyshR1P3jwII4dO5aQut9www3Yt28fPvjgg8DjqquuwvTp0wO/J1N9r7322pDL2T/55BMMGzYMADBixAiUlpYq6tve3o7du3cnrG2cPXsWGRnKQzIzMzPwv9tkrLNET90qKyvR2tqKpqamwD5vvPEGfD4fKioqTK+z1Ik5dOgQtm3bhuLiYsXryVTfGTNm4MMPP1Qcf2VlZXjwwQfx+uuvJ119s7OzcfXVV/d7DCbbdxxZWKJnG8fb+vXrhZycHGHt2rXCgQMHhPvuu08oKioS3G53Qus1a9Yswel0Cjt27BC+/PLLwOPs2bOBfX72s58J5eXlwhtvvCG89957QmVlpVBZWZnAWivJr1oShOSq7549ewS73S48/vjjwqFDh4QXXnhBGDBggPDHP/4xsM+SJUuEoqIi4eWXXxY+/PBD4eabbxZGjBghdHV1JaTONTU1wgUXXCBs3rxZOHLkiPDSSy8J5513nvDQQw8lRZ3PnDkjvP/++8L7778vABCeeuop4f333w9c5aOnbj/4wQ+EK664Qti9e7fw9ttvCyNHjhTuuOMO0+vb29sr3HTTTcKQIUOEDz74QHEM9vT0JF191QRftZRs9X3ppZeErKws4dlnnxUOHTok/O53vxMyMzOFv/3tb4Eykuk7g6wr5TsygiAIv/vd74Ty8nIhOztbuOaaa4TGxsZEV0kAoPpYs2ZNYJ+uri7hX/7lX4SBAwcKAwYMEG699Vbhyy+/TFylgwR3ZJKtvq+88opw6aWXCjk5OcLo0aOFZ599VvG6z+cTFi5cKLhcLiEnJ0e44YYbhIMHDyaotoLQ3t4uPPDAA0J5ebmQm5srXHjhhcK//uu/Kv5hTWSd33zzTdU2W1NTo7tup06dEu644w6hoKBAcDgcwt133y2cOXPG9PoeOXJE8xh88803k66+atQ6MslW39WrVwsXXXSRkJubK1x22WXCpk2bFGUk23cGWZNNEGTLhhIRERFZSErPkSEiIqLUxo4MERERWRY7MkRERGRZ7MgQERGRZbEjQ0RERJbFjgwRERFZFjsyREREZFnsyBAREZFlsSNDRERElsWODBEREVkWOzJEKey6667D/fffj4ceegiDBg1CaWkpHnvsMQDAjh07kJ2djb/97W+B/ZcuXYqSkhI0NzcnqMZERJFhR4YoxT3//PPIz8/H7t27sXTpUixevBhbt27Fddddh7lz52LGjBloa2vD+++/j4ULF+L3v/89XC5XoqtNRKQLbxpJlMKuu+46eL1exVmXa665Btdffz2WLFmC3t5eVFRUYNSoUdi/fz+uvfZaPPvsswmsMRFRZOyJrgARxde4ceMUzwcPHoyWlhYAQHZ2Nl544QWMGzcOw4YNw7JlyxJRRSKiqHFoiSjFZWVlKZ7bbDb4fL7A8127dgEATp8+jdOnT5taNyKiWLEjQ5TGPvvsM8ybNw/PPfccKioqUFNTo+jkEBElO3ZkiNKU1+vFnXfeierqatx9991Ys2YNPvzwQ/zmN79JdNWIiHRjR4YoTT3++OM4evQonnnmGQDi3Jlnn30WjzzyCP7+978nuHZERPrwqiUiIiKyLJ6RISIiIstiR4aIiIgsix0ZIiIisix2ZIiIiMiy2JEhIiIiy2JHhoiIiCyLHRkiIiKyLHZkiIiIyLLYkSEiIiLLYkeGiIiILIsdGSIiIrIsdmSIiIjIsv4f+fzdXajvv0QAAAAASUVORK5CYII=",
      "text/plain": [
       "<Figure size 640x480 with 2 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "topo.depth.plot(vmin=-50, vmax=500)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### *Plot the ocean mask:*"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.collections.QuadMesh at 0x1542689c8ad0>"
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAhoAAAG2CAYAAAA5jy8uAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8ekN5oAAAACXBIWXMAAA9hAAAPYQGoP6dpAAA1kElEQVR4nO3dfXRU1b3/8c+EJBMkJAiYCcEE0GIDKg8mEiP4k9rU1GtRCrdFipBSikuaKCStD2kFLFaDtAK1DVAoD3ZVLpYuwfoEF6OgXgJogCpXG1GRpOAkUpsEYpPQzPn9ocx1JGCSMyd7hnm/1tprMfucfc43LJJ8+e69z3FZlmUJAADAAVGmAwAAAOcuEg0AAOAYEg0AAOAYEg0AAOAYEg0AAOAYEg0AAOAYEg0AAOAYEg0AAOAYEg0AAOAYEg0AAOAYEg0AACLAyy+/rHHjxiklJUUul0ubN2/+0jHbt2/XFVdcIbfbra985Stat25dh+9LogEAQARobGzU8OHDVVpa2q7zDx06pBtvvFFf+9rXtH//fs2ZM0c//OEPtXXr1g7d18VL1QAAiCwul0ubNm3S+PHjz3jOPffco2effVYHDhzw991yyy2qq6vTli1b2n2vaDuBhgOfz6ejR4+qZ8+ecrlcpsMBAIQwy7J0/PhxpaSkKCrKuaJ/U1OTWlpabF/HsqzTfre53W653W7b1y4vL1dOTk5AX25urubMmdOh65zzicbRo0eVmppqOgwAQBiprq7WhRde6Mi1m5qaNGhAvLy1rbavFR8frxMnTgT0zZ8/X/fff7/ta3u9Xnk8noA+j8ejhoYG/etf/1L37t3bdZ1zPtHo2bOnJOnw3oFKiGdJCgDgzBpO+DTgig/8vzuc0NLSIm9tqw5XDFRCz87/Xmo47tOAjA9UXV2thIQEf38wqhnBdM4nGqdKSgnxUUro2c1wNACAcNAVU+3xPV2K79n5+/j02e+3hISARCNYkpOTVVNTE9BXU1OjhISEdlczpAhINAAACEWtlk+tNrZjtFq+4AXThuzsbD333HMBfdu2bVN2dnaHrsNcAgAABvhk2W4dceLECe3fv1/79++X9On21f3796uqqkqSVFxcrGnTpvnPv/322/X+++/r7rvv1t/+9jctW7ZMf/rTn1RYWNih+5JoAAAQAV5//XWNHDlSI0eOlCQVFRVp5MiRmjdvniTpww8/9CcdkjRo0CA9++yz2rZtm4YPH65HHnlEv//975Wbm9uh+zJ1AgCAAT75ZGfyo6Ojx44dq7M9Oqutp36OHTtW+/bt62hoAUg0AAAwoNWy1GrjmZl2xnYlpk4AAIBjqGgAAGBAZxZ0fnF8OCDRAADAAJ8stUZAosHUCQAAcAwVDQAADGDqBAAAOIZdJwAAADZR0QAAwADfZ83O+HBAogEAgAGtNned2BnblUg0AAAwoNWSzbe3Bi8WJ7FGAwAAOIaKBgAABrBGAwAAOMYnl1rlsjU+HDB1AgAAHENFAwAAA3zWp83O+HBAogEAgAGtNqdO7IztSkydAAAAx1DRAADAgEipaJBoAABggM9yyWfZ2HViY2xXYuoEAAA4hooGAAAGMHUCAAAc06ootdqYWGgNYixOItEAAMAAy+YaDYs1GgAAINJR0QAAwADWaAAAAMe0WlFqtWys0QiTR5AzdQIAABxDRQMAAAN8csln4//7PoVHSYNEAwAAAyJljQZTJwAAwDFUNAAAMMD+YlCmTgAAwBl8ukbDxkvVmDr5cgMHDpTL5Tqt5efnS5KampqUn5+vPn36KD4+XhMnTlRNTY3JkAEAQAcYTTRee+01ffjhh/62bds2SdJ3vvMdSVJhYaGefvppbdy4UTt27NDRo0c1YcIEkyEDABAUvs/eddLZZmfHSlcyOnVywQUXBHxeuHChLr74Yl177bWqr6/X6tWrtX79el133XWSpLVr12rIkCHatWuXrrrqKhMhAwAQFJGyRiNk0qGWlhb98Y9/1A9+8AO5XC5VVFTo5MmTysnJ8Z+Tnp6utLQ0lZeXn/E6zc3NamhoCGgAAIQa32dVCTstHIRMlJs3b1ZdXZ2+//3vS5K8Xq9iY2PVq1evgPM8Ho+8Xu8Zr1NSUqLExER/S01NdTBqAABwNiGTaKxevVo33HCDUlJSbF2nuLhY9fX1/lZdXR2kCAEACJ5Wy2W7hYOQ2N56+PBhvfDCC3ryySf9fcnJyWppaVFdXV1AVaOmpkbJyclnvJbb7Zbb7XYyXABhLDdl+Jees/XoX7sgEkS6U4s6Oz+eNRrttnbtWiUlJenGG2/092VkZCgmJkZlZWX+vsrKSlVVVSk7O9tEmAAAoIOMVzR8Pp/Wrl2rvLw8RUf/XziJiYmaMWOGioqK1Lt3byUkJOiOO+5QdnY2O04AAGHPZ0XJZ2PXiS9Mdp0YTzReeOEFVVVV6Qc/+MFpx5YsWaKoqChNnDhRzc3Nys3N1bJlywxECQBAcEXK1InxROP666+XdYasLC4uTqWlpSotLe3iqAAAQDAYTzQAAIhEPsnWzhFf8EJxFIkGAAAG2H3oFg/sAgAAEY+KBgAABth/10l41ApINAAAMMAnl3yys0aDJ4MCAIAziJSKRnhECQAAwhIVDQAADLD/wK7wqBWQaAAAYIDPcsln5zkaYfL21vBIhwAAQFiiogEAgAE+m1Mn4fLALhINAAAMsP/21vBINMIjSgAAEJaoaAAAYECrXGq18dAtO2O7EokGAAAGMHUCAABgExUNAAAMaJW96Y/W4IXiKBINAAAMiJSpExINAAAM4KVqAADgnFNaWqqBAwcqLi5OWVlZ2rNnz1nPX7p0qb761a+qe/fuSk1NVWFhoZqamtp9PxINAAAMsOSSz0azOrG+44knnlBRUZHmz5+vvXv3avjw4crNzVVtbW2b569fv1733nuv5s+fr7ffflurV6/WE088oZ/+9KftvieJBgAABpyaOrHTOmrx4sWaOXOmpk+frqFDh2rFihU677zztGbNmjbP37lzp0aPHq3vfe97GjhwoK6//npNnjz5S6sgn0eiAQBAGGtoaAhozc3NbZ7X0tKiiooK5eTk+PuioqKUk5Oj8vLyNsdcffXVqqio8CcW77//vp577jn9x3/8R7vjYzEoAAAGBOs18ampqQH98+fP1/3333/a+ceOHVNra6s8Hk9Av8fj0d/+9rc27/G9731Px44d05gxY2RZlv7973/r9ttv79DUCYkGAAAGtNp8e+upsdXV1UpISPD3u91u27Gdsn37dj300ENatmyZsrKy9O6772r27Nl64IEHNHfu3HZdg0QDAIAwlpCQEJBonEnfvn3VrVs31dTUBPTX1NQoOTm5zTFz587V1KlT9cMf/lCSdPnll6uxsVG33Xabfvaznykq6ssTJdZoAABgwKmpEzutI2JjY5WRkaGysrL/i8HnU1lZmbKzs9sc88knn5yWTHTr1k2SZFlWu+5LRQMAAAN8ipLPxv/3OzO2qKhIeXl5yszM1KhRo7R06VI1NjZq+vTpkqRp06apf//+KikpkSSNGzdOixcv1siRI/1TJ3PnztW4ceP8CceXIdEAACBCTJo0SR999JHmzZsnr9erESNGaMuWLf4FolVVVQEVjPvuu08ul0v33Xefjhw5ogsuuEDjxo3Tgw8+2O57uqz21j7CVENDgxITE/XPdy5SQs/2ZV8Azl25KcO/9JytR//aBZEgFDUcb9X5l7yv+vr6dq176NQ9Pvu9NOuVCXLHx3T6Os0nTmr5NU86GmswUNEAAMCAYG1vDXUkGgAAGGDZfHurxUvVAABApKOiAQCAAa1yqbUTL0b7/PhwQKIBAIABPsveOgtfmGzlMD51cuTIEd16663q06ePunfvrssvv1yvv/66/7hlWZo3b5769eun7t27KycnRwcPHjQYMQAAaC+jicY///lPjR49WjExMXr++ef11ltv6ZFHHtH555/vP2fRokV69NFHtWLFCu3evVs9evRQbm6umpqaDEYOAIA9vs8Wg9pp4cDo1MnDDz+s1NRUrV271t83aNAg/58ty9LSpUt133336eabb5Yk/eEPf5DH49HmzZt1yy23dHnMAAAEg08u+Wyss7AztisZTYf+8pe/KDMzU9/5zneUlJSkkSNHatWqVf7jhw4dktfrVU5Ojr8vMTFRWVlZKi8vb/Oazc3NamhoCGgAAMAMo4nG+++/r+XLl2vw4MHaunWrZs2apTvvvFOPPfaYJMnr9UqS/9Gop3g8Hv+xLyopKVFiYqK/paamOvtFAADQCa2Wy3YLB0anTnw+nzIzM/XQQw9JkkaOHKkDBw5oxYoVysvL69Q1i4uLVVRU5P/c0NBAsgEACDl211mEyxoNo1H269dPQ4cODegbMmSIqqqqJEnJycmSpJqamoBzampq/Me+yO12KyEhIaABAAAzjCYao0ePVmVlZUDfO++8owEDBkj6dGFocnKyysrK/McbGhq0e/duZWdnd2msAAAEk08u//tOOtXCZDGo0amTwsJCXX311XrooYf03e9+V3v27NHKlSu1cuVKSZLL5dKcOXP0i1/8QoMHD9agQYM0d+5cpaSkaPz48SZDBwDAFsvmrhOLROPLXXnlldq0aZOKi4u1YMECDRo0SEuXLtWUKVP859x9991qbGzUbbfdprq6Oo0ZM0ZbtmxRXFycwcgBALCHt7d2kW9961v61re+dcbjLpdLCxYs0IIFC7owKgAAEAzGEw0AACJRpOw6IdEAAMCASJk6CY90CAAAhCUqGgAAGBAp7zoh0QAAwACmTgAAAGyiogEAgAGRUtEg0QAAwIBISTSYOgEAAI6hogEAgAGRUtEg0QAAwABL9raoWsELxVEkGgAAGBApFQ3WaAAAAMdQ0QAAwIBIqWiQaAAAYECkJBpMnQAAAMdQ0QAAwIBIqWiQaAAAYIBluWTZSBbsjO1KTJ0AAADHUNEAAMAAn1y2HthlZ2xXItEAAMCASFmjwdQJAABwDBUNAAAMiJTFoCQaAAAYEClTJyQaAAAYECkVDdZoAAAAx1DRAADAAMvm1Em4VDRINAAAMMCSZFn2xocDpk4AAIBjqGgAAGCATy65eDIoAABwArtOAAAAbKKiAQCAAT7LJRcP7AIAAE6wLJu7TsJk2wlTJwAAwDFGE437779fLpcroKWnp/uPNzU1KT8/X3369FF8fLwmTpyompoagxEDABAcpxaD2mnhwHhF49JLL9WHH37ob6+++qr/WGFhoZ5++mlt3LhRO3bs0NGjRzVhwgSD0QIAEByRkmgYX6MRHR2t5OTk0/rr6+u1evVqrV+/Xtddd50kae3atRoyZIh27dqlq666qqtDBQAgaCJlMajxisbBgweVkpKiiy66SFOmTFFVVZUkqaKiQidPnlROTo7/3PT0dKWlpam8vPyM12tublZDQ0NAAwAAZhhNNLKysrRu3Tpt2bJFy5cv16FDh3TNNdfo+PHj8nq9io2NVa9evQLGeDweeb3eM16zpKREiYmJ/paamurwVwEAQMed2nVip4UDo1MnN9xwg//Pw4YNU1ZWlgYMGKA//elP6t69e6euWVxcrKKiIv/nhoYGkg0AQMj5NFmw82TQIAbjIONTJ5/Xq1cvXXLJJXr33XeVnJyslpYW1dXVBZxTU1PT5pqOU9xutxISEgIaAAAwI6QSjRMnTui9995Tv379lJGRoZiYGJWVlfmPV1ZWqqqqStnZ2QajBADAPnaddIGf/OQnGjdunAYMGKCjR49q/vz56tatmyZPnqzExETNmDFDRUVF6t27txISEnTHHXcoOzubHScAgLBnfdbsjA8HRhONv//975o8ebL+8Y9/6IILLtCYMWO0a9cuXXDBBZKkJUuWKCoqShMnTlRzc7Nyc3O1bNkykyEDAIAOMJpobNiw4azH4+LiVFpaqtLS0i6KCACArhEpr4k3/sAuAAAiUoTMnYTUYlAAACKG3YWgnaxolJaWauDAgYqLi1NWVpb27Nlz1vPr6uqUn5+vfv36ye1265JLLtFzzz3X7vtR0QAAIEI88cQTKioq0ooVK5SVlaWlS5cqNzdXlZWVSkpKOu38lpYWfeMb31BSUpL+/Oc/q3///jp8+PBpD9M8GxINAAAMsPt0z86MXbx4sWbOnKnp06dLklasWKFnn31Wa9as0b333nva+WvWrNHHH3+snTt3KiYmRpI0cODADt2TqRMAAAwI1nM0vvh+r+bm5jbv19LSooqKioB3iEVFRSknJ+eM7xD7y1/+ouzsbOXn58vj8eiyyy7TQw89pNbW1nZ/nSQaAACEsdTU1IB3fJWUlLR53rFjx9Ta2iqPxxPQf7Z3iL3//vv685//rNbWVj333HOaO3euHnnkEf3iF79od3xMnQAAYIKNBZ3+8ZKqq6sDXrfhdrvtRubn8/mUlJSklStXqlu3bsrIyNCRI0f0y1/+UvPnz2/XNUg0AAAwIFhrNNr7Xq++ffuqW7duqqmpCeg/2zvE+vXrp5iYGHXr1s3fN2TIEHm9XrW0tCg2NvZL78vUCQAAESA2NlYZGRkB7xDz+XwqKys74zvERo8erXfffVc+n8/f984776hfv37tSjIkEg0AAMywgtA6qKioSKtWrdJjjz2mt99+W7NmzVJjY6N/F8q0adNUXFzsP3/WrFn6+OOPNXv2bL3zzjt69tln9dBDDyk/P7/d92TqBAAAA0w8gnzSpEn66KOPNG/ePHm9Xo0YMUJbtmzxLxCtqqpSVNT/1SBSU1O1detWFRYWatiwYerfv79mz56te+65p933JNEAACCCFBQUqKCgoM1j27dvP60vOztbu3bt6vT9SDQAADAlTN5XYgeJBgAABvD2VgAA4Bze3goAAGAPFQ0AAIxwfdbsjA99JBoAAJjA1AkAAIA9VDQAADAhQioaJBoAAJgQpLe3hjqmTgAAgGOoaAAAYECwXhMf6kg0AAAwIULWaDB1AgAAHENFAwAAEyJkMSiJBgAABrisT5ud8eGARAMAABNYowEAAGAPFQ0AAExgjQYAAHAMUycAAAD2UNEAAMCECKlokGgAAGBChCQaTJ0AAADHhEyisXDhQrlcLs2ZM8ff19TUpPz8fPXp00fx8fGaOHGiampqzAUJAECwnNp1YqeFgZBINF577TX97ne/07BhwwL6CwsL9fTTT2vjxo3asWOHjh49qgkTJhiKEgCA4Dn1ZFA7LRx0KtFYu3atPvnkk6AEcOLECU2ZMkWrVq3S+eef7++vr6/X6tWrtXjxYl133XXKyMjQ2rVrtXPnTu3atSso9wYAAM7qVKJx7733Kjk5WTNmzNDOnTttBZCfn68bb7xROTk5Af0VFRU6efJkQH96errS0tJUXl5u654AABhnBaGFgU4lGkeOHNFjjz2mY8eOaezYsUpPT9fDDz8sr9fboets2LBBe/fuVUlJyWnHvF6vYmNj1atXr4B+j8dz1vs0NzeroaEhoAEAADM6lWhER0fr29/+tp566ilVV1dr5syZevzxx5WWlqabbrpJTz31lHw+31mvUV1drdmzZ+vxxx9XXFxcp4JvS0lJiRITE/0tNTU1aNcGACBYXLK5RsP0F9BOtheDejwejRkzRtnZ2YqKitKbb76pvLw8XXzxxdq+ffsZx1VUVKi2tlZXXHGFoqOjFR0drR07dujRRx9VdHS0PB6PWlpaVFdXFzCupqZGycnJZ7xucXGx6uvr/a26utrulwgAADqp04lGTU2NfvWrX+nSSy/V2LFj1dDQoGeeeUaHDh3SkSNH9N3vfld5eXlnHP/1r39db775pvbv3+9vmZmZmjJliv/PMTExKisr84+prKxUVVWVsrOzz3hdt9uthISEgAYAQMiJkO2tnXoy6Lhx47R161ZdcsklmjlzpqZNm6bevXv7j/fo0UM//vGP9ctf/vKM1+jZs6cuu+yygL4ePXqoT58+/v4ZM2aoqKhIvXv3VkJCgu644w5lZ2frqquu6kzYAACEjgh5MminEo2kpCTt2LHjrJWFCy64QIcOHep0YJK0ZMkSRUVFaeLEiWpublZubq6WLVtm65oAAKDrdCrRWL16tcrKyvTTn/5UtbW1py38XLNmjVwulwYMGNCh635xTUdcXJxKS0tVWlramTABAAhdVDTObMGCBfr5z3+uzMxM9evXTy5XeMwTAQAQKuw+3TNcngzaqURj+fLlWrdunaZOnRrseAAAwDmkU4lGS0uLrr766mDHAgBA5IiQqZNObW/94Q9/qPXr1wc7FgAAIkeEPIK8UxWNpqYmrVy5Ui+88IKGDRummJiYgOOLFy8OSnAAACC8dSrReOONNzRixAhJ0oEDBwKOsTAUAIAvx2LQs3jppZeCHQcAAJHF7tM9z+UngwIAAJtYDAoAAGAPFQ0AAAxgjQYAAHAOUycAAAD2UNEAAMAEm1Mn4VLRINEAAMAEpk4AAADsoaIBAIAJEVLRINEAAMCASNneytQJAABwDIkGAABwDFMnAACYwBoNAADgFNZoAAAA2ERFAwAAU8KkKmEHiQYAACZEyBoNpk4AAIBjqGgAAGBApCwGJdEAAMAEpk4AAADsoaIBAIABkTJ1QkUDAAATrCC0TigtLdXAgQMVFxenrKws7dmzp13jNmzYIJfLpfHjx3fofiQaAABEiCeeeEJFRUWaP3++9u7dq+HDhys3N1e1tbVnHffBBx/oJz/5ia655poO35NEAwAAEwxUNBYvXqyZM2dq+vTpGjp0qFasWKHzzjtPa9asOeOY1tZWTZkyRT//+c910UUXdfieJBoAABhwao2GnSZJDQ0NAa25ubnN+7W0tKiiokI5OTn+vqioKOXk5Ki8vPyMcS5YsEBJSUmaMWNGp75OEg0AAEwIUkUjNTVViYmJ/lZSUtLm7Y4dO6bW1lZ5PJ6Afo/HI6/X2+aYV199VatXr9aqVas6/WWy6wQAgDBWXV2thIQE/2e32x2U6x4/flxTp07VqlWr1Ldv305fh0QDAAATgvTAroSEhIBE40z69u2rbt26qaamJqC/pqZGycnJp53/3nvv6YMPPtC4ceP8fT6fT5IUHR2tyspKXXzxxV96X6NTJ8uXL9ewYcP8f0nZ2dl6/vnn/cebmpqUn5+vPn36KD4+XhMnTjztLwgAgHAUrDUa7RUbG6uMjAyVlZX5+3w+n8rKypSdnX3a+enp6XrzzTe1f/9+f7vpppv0ta99Tfv371dqamq77mu0onHhhRdq4cKFGjx4sCzL0mOPPaabb75Z+/bt06WXXqrCwkI9++yz2rhxoxITE1VQUKAJEybof/7nf0yGDQBAWCoqKlJeXp4yMzM1atQoLV26VI2NjZo+fbokadq0aerfv79KSkoUFxenyy67LGB8r169JOm0/rMxmmh8vhwjSQ8++KCWL1+uXbt26cILL9Tq1au1fv16XXfddZKktWvXasiQIdq1a5euuuoqEyEDABAcBt51MmnSJH300UeaN2+evF6vRowYoS1btvgXiFZVVSkqKriTHSGzRqO1tVUbN25UY2OjsrOzVVFRoZMnTwZsw0lPT1daWprKy8vPmGg0NzcHbO1paGhwPHYAADrK1CPICwoKVFBQ0Oax7du3n3XsunXrOnw/49tb33zzTcXHx8vtduv222/Xpk2bNHToUHm9XsXGxvrLNKecbRuOJJWUlARs82nvHBIAAAg+44nGV7/6Ve3fv1+7d+/WrFmzlJeXp7feeqvT1ysuLlZ9fb2/VVdXBzFaAACCxNC7Trqa8amT2NhYfeUrX5EkZWRk6LXXXtOvf/1rTZo0SS0tLaqrqwuoapxpG84pbrc7aHuIAQBwjIE1GiYYr2h8kc/nU3NzszIyMhQTExOwDaeyslJVVVVtbsMBAAChx2hFo7i4WDfccIPS0tJ0/PhxrV+/Xtu3b9fWrVuVmJioGTNmqKioSL1791ZCQoLuuOMOZWdns+MEABD2XJ81O+PDgdFEo7a2VtOmTdOHH36oxMREDRs2TFu3btU3vvENSdKSJUsUFRWliRMnqrm5Wbm5uVq2bJnJkAEACI4ImToxmmisXr36rMfj4uJUWlqq0tLSLooIAICuYWp7a1cLuTUaAADg3GF81wkAABGJqRMAAOCoMEkW7GDqBAAAOIaKBgAABkTKYlASDQAATIiQNRpMnQAAAMdQ0QAAwACmTgAAgHOYOgEAALCHigYAAAYwdQIAAJwTIVMnJBoAAJgQIYkGazQAAIBjqGgAAGAAazQAAIBzmDoBAACwh4oGAAAGuCxLLqvzZQk7Y7sSiQYAACYwdQIAAGAPFQ0AAAxg1wkAAHAOUycAAAD2UNEAAMAApk4AAIBzImTqhEQDAAADIqWiwRoNAADgGCoaAACYwNQJAABwUrhMf9jB1AkAAHAMFQ0AAEywrE+bnfFhgEQDAAAD2HUCAABgExUNAABMYNcJAABwisv3abMzPhwYnTopKSnRlVdeqZ49eyopKUnjx49XZWVlwDlNTU3Kz89Xnz59FB8fr4kTJ6qmpsZQxAAAoCOMJho7duxQfn6+du3apW3btunkyZO6/vrr1djY6D+nsLBQTz/9tDZu3KgdO3bo6NGjmjBhgsGoAQAIAisILQwYnTrZsmVLwOd169YpKSlJFRUV+n//7/+pvr5eq1ev1vr163XddddJktauXashQ4Zo165duuqqq0yEDQCAbew6MaC+vl6S1Lt3b0lSRUWFTp48qZycHP856enpSktLU3l5eZvXaG5uVkNDQ0ADACDknHqOhp0WBkIm0fD5fJozZ45Gjx6tyy67TJLk9XoVGxurXr16BZzr8Xjk9XrbvE5JSYkSExP9LTU11enQAQDAGYRMopGfn68DBw5ow4YNtq5TXFys+vp6f6uurg5ShAAABM+pqRM7LRyExPbWgoICPfPMM3r55Zd14YUX+vuTk5PV0tKiurq6gKpGTU2NkpOT27yW2+2W2+12OmQAAOyJkOdoGK1oWJalgoICbdq0SS+++KIGDRoUcDwjI0MxMTEqKyvz91VWVqqqqkrZ2dldHS4AAOggoxWN/Px8rV+/Xk899ZR69uzpX3eRmJio7t27KzExUTNmzFBRUZF69+6thIQE3XHHHcrOzmbHCQAgrEXKrhOjicby5cslSWPHjg3oX7t2rb7//e9LkpYsWaKoqChNnDhRzc3Nys3N1bJly7o4UgAAgoy3tzrPasdfUlxcnEpLS1VaWtoFEQEAgGAKicWgAABEGqZOAACAc9h1AgAAYA8VDQAADGDqBAAAOMdnfdrsjA8DJBoAAJjAGg0AAAB7qGgAAGCASzbXaAQtEmeRaAAAYEKEPBmUqRMAAOAYEg0AAAw4tb3VTuuM0tJSDRw4UHFxccrKytKePXvOeO6qVat0zTXX6Pzzz9f555+vnJycs57fFhINAABMsILQOuiJJ55QUVGR5s+fr71792r48OHKzc1VbW1tm+dv375dkydP1ksvvaTy8nKlpqbq+uuv15EjR9p9TxINAAAixOLFizVz5kxNnz5dQ4cO1YoVK3TeeedpzZo1bZ7/+OOP60c/+pFGjBih9PR0/f73v5fP51NZWVm770miAQCAAS7Lst0kqaGhIaA1Nze3eb+WlhZVVFQoJyfH3xcVFaWcnByVl5e3K+ZPPvlEJ0+eVO/evdv9dZJoAABggi8ITVJqaqoSExP9raSkpM3bHTt2TK2trfJ4PAH9Ho9HXq+3XSHfc889SklJCUhWvgzbWwEACGPV1dVKSEjwf3a73Y7cZ+HChdqwYYO2b9+uuLi4do8j0QAAwIDPT390drwkJSQkBCQaZ9K3b19169ZNNTU1Af01NTVKTk4+69hf/epXWrhwoV544QUNGzasQ3EydQIAgAldvOskNjZWGRkZAQs5Ty3szM7OPuO4RYsW6YEHHtCWLVuUmZnZsZuKigYAAGYYeDJoUVGR8vLylJmZqVGjRmnp0qVqbGzU9OnTJUnTpk1T//79/es8Hn74Yc2bN0/r16/XwIED/Ws54uPjFR8f3657kmgAABAhJk2apI8++kjz5s2T1+vViBEjtGXLFv8C0aqqKkVF/d9kx/Lly9XS0qL//M//DLjO/Pnzdf/997frniQaAAAYYOfpnqfGd0ZBQYEKCgraPLZ9+/aAzx988EHnbvI5JBoAAJjAS9UAAADsoaIBAIABLt+nzc74cECiAQCACUydAAAA2ENFAwAAEzr5qveA8WGARAMAAAOC9QjyUMfUCQAAcAwVDQAATIiQxaAkGgAAmGBJsrNFNTzyDBINAABMYI0GAACATVQ0AAAwwZLNNRpBi8RRJBoAAJgQIYtBjU6dvPzyyxo3bpxSUlLkcrm0efPmgOOWZWnevHnq16+funfvrpycHB08eNBMsAAAoMOMJhqNjY0aPny4SktL2zy+aNEiPfroo1qxYoV2796tHj16KDc3V01NTV0cKQAAQeYLQgsDRqdObrjhBt1www1tHrMsS0uXLtV9992nm2++WZL0hz/8QR6PR5s3b9Ytt9zSlaECABBU7Dox7NChQ/J6vcrJyfH3JSYmKisrS+Xl5Wcc19zcrIaGhoAGAADMCNlEw+v1SpI8Hk9Av8fj8R9rS0lJiRITE/0tNTXV0TgBAOiUU4tB7bQwELKJRmcVFxervr7e36qrq02HBADA6Ug0zEpOTpYk1dTUBPTX1NT4j7XF7XYrISEhoAEAADNC9jkagwYNUnJyssrKyjRixAhJUkNDg3bv3q1Zs2aZDQ5A2Np69K+mQwA+FSHP0TCaaJw4cULvvvuu//OhQ4e0f/9+9e7dW2lpaZozZ45+8YtfaPDgwRo0aJDmzp2rlJQUjR8/3lzQAAAEg0+Sy+b4MGA00Xj99df1ta99zf+5qKhIkpSXl6d169bp7rvvVmNjo2677TbV1dVpzJgx2rJli+Li4kyFDABAUETK9lajicbYsWNlneUvyuVyacGCBVqwYEEXRgUAAIIlZNdoAABwTmONBgAAcIzPklw2kgVfeCQaIbu9FQAAhD8qGgAAmMDUCQAAcI7dp3uGR6LB1AkAAHAMFQ0AAExg6gQAADjGZ8nW9Ae7TgAAQKSjogEAgAmW79NmZ3wYINEAAMAE1mgAAADHsEYDAADAHioaAACYwNQJAABwjCWbiUbQInEUUycAAMAxVDQAADCBqRMAAOAYn0+SjWdh+MLjORpMnQAAAMdQ0QAAwASmTgAAgGMiJNFg6gQAADiGigYAACZEyCPISTQAADDAsnyybLyB1c7YrkSiAQCACZZlryrBGg0AABDpqGgAAGCCZXONRphUNEg0AAAwweeTXDbWWYTJGg2mTgAAgGOoaAAAYAJTJwAAwCmWzyfLxtRJuGxvZeoEAAA4hooGAAAmMHUCAAAc47Mk17mfaDB1AgAAHBMWiUZpaakGDhyouLg4ZWVlac+ePaZDAgDAHsv69FkYnW5UNILiiSeeUFFRkebPn6+9e/dq+PDhys3NVW1trenQAADoNMtn2W7hIOQTjcWLF2vmzJmaPn26hg4dqhUrVui8887TmjVrTIcGAEDn2apm+HgyaDC0tLSooqJCOTk5/r6oqCjl5OSovLy8zTHNzc1qaGgIaAAA4FMdXY6wceNGpaenKy4uTpdffrmee+65Dt0vpBONY8eOqbW1VR6PJ6Df4/HI6/W2OaakpESJiYn+lpqa2hWhAgDQISamTjq6HGHnzp2aPHmyZsyYoX379mn8+PEaP368Dhw40O57hnSi0RnFxcWqr6/3t+rqatMhAQBwOgNTJx1djvDrX/9a3/zmN3XXXXdpyJAheuCBB3TFFVfot7/9bbvvGdLP0ejbt6+6deummpqagP6amholJye3Ocbtdsvtdvs/W5+tym04ER5zWQAAc079rrC6YEfHv3XS1vO6/q2TknTaEoEv/h485dRyhOLiYn/fly1HKC8vV1FRUUBfbm6uNm/e3O44QzrRiI2NVUZGhsrKyjR+/HhJks/nU1lZmQoKCtp1jePHj0uSBlzxgUNRAgDONcePH1diYqIj146NjVVycrJe9XZsrUNb4uPjT1siMH/+fN1///2nnXu25Qh/+9vf2ry+1+vt0PKFtoR0oiFJRUVFysvLU2ZmpkaNGqWlS5eqsbFR06dPb9f4lJQUvfXWWxo6dKiqq6uVkJDgcMT2NTQ0KDU1lXgdQrzOCrd4pfCLmXidY1mWjh8/rpSUFMfuERcXp0OHDqmlpcX2tSzLksvlCuhrq5phUsgnGpMmTdJHH32kefPmyev1asSIEdqyZctpGdaZREVFqX///pKkhISEkP9H/nnE6yzidVa4xSuFX8zE6wynKhmfFxcXp7i4OMfv83mdWY6QnJzcofPbEhaLQQsKCnT48GE1Nzdr9+7dysrKMh0SAABh5fPLEU45tRwhOzu7zTHZ2dkB50vStm3bznh+W0K+ogEAAILjy5YjTJs2Tf3791dJSYkkafbs2br22mv1yCOP6MYbb9SGDRv0+uuva+XKle2+Z0QkGm63W/Pnzw+5easzIV5nEa+zwi1eKfxiJl501pctR6iqqlJU1P9Ndlx99dVav3697rvvPv30pz/V4MGDtXnzZl122WXtvqfL6oo9PAAAICKFxRoNAAAQnkg0AACAY0g0AACAY0g0AACAYyIi0ejoK3G7QklJia688kr17NlTSUlJGj9+vCorKwPOaWpqUn5+vvr06aP4+HhNnDjxtAenmLJw4UK5XC7NmTPH3xdq8R45ckS33nqr+vTpo+7du+vyyy/X66+/7j9uWZbmzZunfv36qXv37srJydHBgweNxdva2qq5c+dq0KBB6t69uy6++GI98MADAe9cMBnzyy+/rHHjxiklJUUul+u0dx20J7aPP/5YU6ZMUUJCgnr16qUZM2boxIkTXR7vyZMndc899+jyyy9Xjx49lJKSomnTpuno0aMhGe8X3X777XK5XFq6dGlIx/v222/rpptuUmJionr06KErr7xSVVVV/uOh9jMDzjjnE42OvhK3q+zYsUP5+fnatWuXtm3bppMnT+r6669XY2Oj/5zCwkI9/fTT2rhxo3bs2KGjR49qwoQJBqP+1Guvvabf/e53GjZsWEB/KMX7z3/+U6NHj1ZMTIyef/55vfXWW3rkkUd0/vnn+89ZtGiRHn30Ua1YsUK7d+9Wjx49lJubq6amJiMxP/zww1q+fLl++9vf6u2339bDDz+sRYsW6Te/+U1IxNzY2Kjhw4ertLS0zePtiW3KlCn63//9X23btk3PPPOMXn75Zd12221dHu8nn3yivXv3au7cudq7d6+efPJJVVZW6qabbgo4L1Ti/bxNmzZp165dbT4iO5Tife+99zRmzBilp6dr+/bteuONNzR37tyAp2GG0s8MOMg6x40aNcrKz8/3f25tbbVSUlKskpISg1Gdrra21pJk7dixw7Isy6qrq7NiYmKsjRs3+s95++23LUlWeXm5qTCt48ePW4MHD7a2bdtmXXvttdbs2bMtywq9eO+55x5rzJgxZzzu8/ms5ORk65e//KW/r66uznK73dZ//dd/dUWIp7nxxhutH/zgBwF9EyZMsKZMmWJZVmjFLMnatGmT/3N7YnvrrbcsSdZrr73mP+f555+3XC6XdeTIkS6Nty179uyxJFmHDx+2LCs04/373/9u9e/f3zpw4IA1YMAAa8mSJf5joRbvpEmTrFtvvfWMY0LtZwacc05XNE69EjcnJ8ff92WvxDWlvr5ektS7d29JUkVFhU6ePBkQe3p6utLS0ozGnp+frxtvvDEgLin04v3LX/6izMxMfec731FSUpJGjhypVatW+Y8fOnRIXq83IN7ExERlZWUZ+/u9+uqrVVZWpnfeeUeS9Ne//lWvvvqqbrjhhpCN+ZT2xFZeXq5evXopMzPTf05OTo6ioqK0e/fuLo/5i+rr6+VyudSrVy9JoRevz+fT1KlTddddd+nSSy897Xgoxevz+fTss8/qkksuUW5urpKSkpSVlRUwvRJqPzPgnHM60TjbK3E78opbp/l8Ps2ZM0ejR4/2P23N6/UqNjbW/0PvFJOxb9iwQXv37vU/mvbzQi3e999/X8uXL9fgwYO1detWzZo1S3feeacee+wxf7yn4guFeCXp3nvv1S233KL09HTFxMRo5MiRmjNnjqZMmSIpNGM+pT2xeb1eJSUlBRyPjo5W7969jcff1NSke+65R5MnT/a/9CvU4n344YcVHR2tO++8s83joRRvbW2tTpw4oYULF+qb3/ym/vu//1vf/va3NWHCBO3YscMfbyj9zIBzIuIR5KEuPz9fBw4c0Kuvvmo6lDOqrq7W7NmztW3bti5/42Bn+Hw+ZWZm6qGHHpIkjRw5UgcOHNCKFSuUl5dnOLq2/elPf9Ljjz+u9evX69JLL9X+/fs1Z84cpaSkhGzM54KTJ0/qu9/9rizL0vLly02H06aKigr9+te/1t69e097JXgo8vl8kqSbb75ZhYWFkqQRI0Zo586dWrFiha699lqT4aGLndMVjc68ErerFRQU6JlnntFLL72kCy+80N+fnJyslpYW1dXVBZxvKvaKigrV1tbqiiuuUHR0tKKjo7Vjxw49+uijio6OlsfjCal4+/Xrp6FDhwb0DRkyxL/i/VRMofRv46677vJXNS6//HJNnTpVhYWF/gpSKMZ8SntiS05OPm0R9r///W99/PHHxuI/lWQcPnxY27ZtC3iFeSjF+8orr6i2tlZpaWn+77/Dhw/rxz/+sQYOHBhy8fbt21fR0dFf+j0YSj8z4JxzOtHozCtxu4plWSooKNCmTZv04osvatCgQQHHMzIyFBMTExB7ZWWlqqqqjMT+9a9/XW+++ab279/vb5mZmZoyZYr/z6EU7+jRo0/bLvzOO+9owIABkqRBgwYpOTk5IN6Ghgbt3r3b2L+NTz75JOBlRpLUrVs3//8OQzHmU9oTW3Z2turq6lRRUeE/58UXX5TP51NWVlaXx3wqyTh48KBeeOEF9enTJ+B4KMU7depUvfHGGwHffykpKbrrrru0devWkIs3NjZWV1555Vm/B0PtZxwcZHo1qtM2bNhgud1ua926ddZbb71l3XbbbVavXr0sr9drNK5Zs2ZZiYmJ1vbt260PP/zQ3z755BP/ObfffruVlpZmvfjii9brr79uZWdnW9nZ2QajDvT5XSeWFVrx7tmzx4qOjrYefPBB6+DBg9bjjz9unXfeedYf//hH/zkLFy60evXqZT311FPWG2+8Yd18883WoEGDrH/9619GYs7Ly7P69+9vPfPMM9ahQ4esJ5980urbt6919913h0TMx48ft/bt22ft27fPkmQtXrzY2rdvn3+XRnti++Y3v2mNHDnS2r17t/Xqq69agwcPtiZPntzl8ba0tFg33XSTdeGFF1r79+8P+B5sbm4OuXjb8sVdJ6EW75NPPmnFxMRYK1eutA4ePGj95je/sbp162a98sor/muE0s8MOOecTzQsy7J+85vfWGlpaVZsbKw1atQoa9euXaZDsiS12dauXes/51//+pf1ox/9yDr//POt8847z/r2t79tffjhh+aC/oIvJhqhFu/TTz9tXXbZZZbb7bbS09OtlStXBhz3+XzW3LlzLY/HY7ndbuvrX/+6VVlZaShay2poaLBmz55tpaWlWXFxcdZFF11k/exnPwv4xWcy5pdeeqnNf7N5eXntju0f//iHNXnyZCs+Pt5KSEiwpk+fbh0/frzL4z106NAZvwdfeumlkIu3LW0lGqEW7+rVq62vfOUrVlxcnDV8+HBr8+bNAdcItZ8ZcAaviQcAAI45p9doAAAAs0g0AACAY0g0AACAY0g0AACAY0g0AACAY0g0AACAY0g0AACAY0g0AACAY0g0AACAY0g0AACAY0g0gHPY2LFjdeedd+ruu+9W7969lZycrPvvv1+StH37dsXGxuqVV17xn79o0SIlJSWd9rp3AOgsEg3gHPfYY4+pR48e2r17txYtWqQFCxZo27ZtGjt2rObMmaOpU6eqvr5e+/bt09y5c/X73/9eHo/HdNgAzhG8VA04h40dO1atra0BVYtRo0bpuuuu08KFC9XS0qKsrCxdcsklOnDggEaPHq2VK1cajBjAuSbadAAAnDVs2LCAz/369VNtba0kKTY2Vo8//riGDRumAQMGaMmSJSZCBHAOY+oEOMfFxMQEfHa5XPL5fP7PO3fulCR9/PHH+vjjj7s0NgDnPhINIIK99957Kiws1KpVq5SVlaW8vLyAJAQA7CLRACJUa2urbr31VuXm5mr69Olau3at3njjDT3yyCOmQwNwDiHRACLUgw8+qMOHD+t3v/udpE/XbqxcuVL33Xef/vrXvxqODsC5gl0nAADAMVQ0AACAY0g0AACAY0g0AACAY0g0AACAY0g0AACAY0g0AACAY0g0AACAY0g0AACAY0g0AACAY0g0AACAY0g0AACAY0g0AACAY/4/ZY7YH7KOP68AAAAASUVORK5CYII=",
      "text/plain": [
       "<Figure size 640x480 with 2 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "topo.tmask.plot()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## 4. Save the grid and bathymetry files"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [],
   "source": [
    "# MOM6 supergrid file:\n",
    "grid.write_supergrid(\"./ocean_hgrid_3.nc\")\n",
    "\n",
    "# MOM6 topography file:\n",
    "topo.write_topo(\"./ocean_topog_3.nc\")\n",
    "\n",
    "# CICE grid file:\n",
    "topo.write_cice_grid(\"./cice_grid_3.nc\")\n",
    "\n",
    "# SCRIP grid file (for runoff remapping, if needed):\n",
    "topo.write_scrip_grid(\"./scrip_grid_3.nc\")\n",
    "\n",
    "# ESMF mesh file:\n",
    "topo.write_esmf_mesh(\"./ESMF_mesh_3.nc\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python [conda env:mom6_forge]",
   "language": "python",
   "name": "conda-env-mom6_forge-py"
  },
  "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.15"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}