Skip to content

Instantly share code, notes, and snippets.

Show Gist options
  • Save gustavo-marques/3376d14da22140448006c2b4a2b12049 to your computer and use it in GitHub Desktop.
Save gustavo-marques/3376d14da22140448006c2b4a2b12049 to your computer and use it in GitHub Desktop.
CVMix convection: results from a test case with unstable initial conditions
{
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# Convective mixing due to unstable initial conditions via CVMix\n",
"\n",
"This notebook shows results from a test cases where convection occurs due to a initially unstable water column."
]
},
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Populating the interactive namespace from numpy and matplotlib\n"
]
}
],
"source": [
"import numpy as np\n",
"from netCDF4 import Dataset\n",
"import matplotlib.pyplot as plt\n",
"%pylab inline\n",
"pylab.rcParams['figure.figsize'] = (16.0, 4.0)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Initial Conditions \n",
"\n",
"The initial water column is unstable near the surface: $\\frac{\\partial \\theta_1}{\\partial z} = -0.01$; $\\frac{\\partial \\theta_2}{\\partial z} = 0.01$; salt(z) = 35."
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.text.Text at 0x7f66740cb3c8>"
]
},
"execution_count": 2,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAZoAAAEoCAYAAABhBDcdAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XmcFOW1//HPYd9B9nUYYEAFRIRRcWXRKG5BoyaYRKNG0YjL/SU3JsZcNUa92W4WQVQSTULcojFGk2gWHBBEEUFxwY2eBZhhB9lhmOX8/qia0A6z9Uz3dPfM9/1yXt39VNVTp6aAY1U9/Rxzd0RERBKlRbIDEBGRpk2JRkREEkqJRkREEkqJRkREEkqJRkREEkqJRkREEkqJRkREEkqJRkREEkqJRpoUM8s0Mzez3yVym3j3Y2YdzOwOM/vIzA6Y2Tozu8/MWjckpmQws5vN7AMz2x/+Pv6rut9NvH73ktqUaCRlhP/gJGSqilT+B83M+gFvAt8H3gF+BewCbgNmJzG0mJnZdIL4DwC/BH4ALE1qUJJ0rZIdgEicFQFHAzsTvE1cmFkb4K/AYGCyuy8J238IrAKuMbM73X1jY8dWT+dXvLr7+orG8MosKb9jST5d0UiT4u4l7v6Ru29I5DZx9N/AeOA7FUkmjGkP8BzB39HTkhBXffUHiE4y4edk/o4lyZRoJGVF3+4K3z9lZlvDZxjLzez8mraJarsLyA8/fq3iFl34c2VNt9XC5c+aWV74zGGXmS0xs6/G4fjaA98GNgBzq1hlW/jatw59Rf+ujjKzv5jZdjPba2avmtlZNaw/wsz+aGabzazczCZVWveLZrbIzHaGv4P3zOw2M2sbtc5d4W3PyeHn//yOK+8vht/PiWb2JzPbaGYHw+dWD5tZ/7r2IalBt84kHQwGlgF5wB+A7sCXgOfN7Ex3X1DL9guBbsAtBM9A/hK1bGUt2z5IcAtrEUFC6AGcC/zBzI509/+J7VA+46IwrkfcvaSK5e3C14Mx9DkEeB14D3gY6Efwu3rJzL7s7n+stP4w4A3gE+BxoD3B8yEAzOw+gmdFW4EngD3AOcB9wNlmdpa7HyT4HQNcSXC+fhBDzIcxs6sJkm8x8AKwDhgOXANcYGYT3H1tQ/Yhjcjd9aOflPgBPPgj+Z/PmRVtwJ2V1j07bH+xUnvFNr+rS3sdlg2roq0N8DJQAgyoSz/VHO/j4fpPAndV8fNGuPzcOvQV/bv6aaVl2WGsnwJdqlj/vmr6PClcvhboG9XeiuC5kgPfq7TNwuhzWJ/zAowgSK6R6N9vuOwMoAx4Ltl/XvVT9x9d0Ug6WAPcE93g7v80s7XACYncsbvnVtF20MweAKYQ/MM3r57dnxq+Tq9lvQ9i6HMncHd0g7svN7PHga8RXEX9PmrxJqq/+rg6fL3HowYjuHupmX2L4MruGoKrm3j6BtAauMXdi6IXuPvLZvYCwVVNZ3ffHed9SwIo0Ug6WOnuZVW0ryP4v+6EMbMM4DsECSWD4NZStAH17Ldj2N8qdx9dxfLOBM9oNrp7QQxdv1XNP74LCRLNcXw20bzj7sXV9DUufM2pvMDdPzGzQmCImXV193iOJqs4pxPN7PgqlvcGWhJc+ayI434lQZRoJB3sqKa9lAQOaDGzoQTPho4AFgP/IrhiKCO45fM1oG1129eiIkEVVbP8LIL/q38xxn43VdNecUXStZr2qlSsW91IsQ0EybIb8R223CN8/XYt63WK4z4lgZRoRKr3TYJ/9K5y999FLzCzywgSTX21CV+ru5q4Knx9NMZ++1TTXjFyrXJCqOkLshXr9gUOu4VIMNCgqj4bqqK/ru6+q8Y1JS1oeLM0FxW33lrGsE1W+PpsFcsmNiyc/1xJHDZ02cwmEDz/eMndl8XY77jwtltlk8LXt2Poq2LdSZUXmFkWMBDId/fqrjjrq2ImgXT6/pDUQIlGmotPCf7vPSOGbQrC10nRjWZ2NsFD8Hpz963Ah8B4MxsT1fdgglFoO4Eb6tF1V+CO6AYzywa+Evb5XAx9VVxNfd/MekX11xL4GcG/H4/UI8bazCYYJfcLMxtReaGZtTEzJaE0oltn0iy4+x4zewM4LRyB9QnBVc4LRH1vpJI5BLewnjGzPwHrgdHAVOBpgu+nNMQ9BEOcXzazx4COwBcJEuJ5MQ4CqLCIYNqaE4ElHPoeTQvgulhuRbn7a2b2E+BW4P3wd7CX4Hs0o4FXgZ/WI8ba9vtR+D2aR4FVZvYPgvPVmuB/FE4DtgBHxXvfkhhKNNKcXA78giBRXAYYUMihLxt+hru/a2aTCRLCeQR/X94BvkAwQKFBicbdnwjnALuVYEjvVoIE9oPKw3pjkA9cD/wofG0LvAXc7e7/rEeM3zGzt4EbgSsI/rHPJZgA9P88+LJm3Ln7Y2b2DvAtgtkGziJIcuuBPwGVv3gqKczcEzJZrog0IjPLJEgyv3f3K5MajEglekYjIiIJpUQjIiIJpUQjIiIJpWc0IiKSUBp1BvTs2dMzMzOTHYaISFpZsWLFVnfvVdt6SjRAZmYmy5cvT3YYIiJpxczW1GU9PaMREZGEUqIREZGEUqIREZGEUqIREZGEUqIREZGEapKJxsymmtnHZhYxs+8mOx4RkeasySWasFbGAwRTmY8ELjOzkcmNSkSk+WpyiQY4AYi4e144hflTwLQkxyTSrP3xzbU8vXxdssOQJGmKiWYAEP0nujBs+wwzm2Fmy81s+ZYtWxotOJHm6Nm3ivjzW4XJDkOSpCkmmjpx97nunu3u2b161TqDgoiI1FNTTDRFwKCozwPDNhERSYKmmGjeBIab2RAzawNMJ6gLLyIiSdDkJtV091IzuxH4J9ASeNTdVyU5LBGRZqvJJRoAd38ReDHZcYiISNO8dSYiIilEiUZERBJKiUZERBJKiUZERBJKiUZERBJKiUZERBJKiUZERBJKiUZERBJKiUZERBJKiUZERBJKiUZERBJKiUZERBJKiUZERBJKiUZERBJKiUZERBJKiUZERBIqLRONmd1lZkVmtjL8OTdq2W1mFjGzj83s7GTGKSIi6V1h8xfu/rPoBjMbCUwHRgH9gflmNsLdy5IRoIiIpOkVTQ2mAU+5e7G75wMR4IQkxyQi0qylc6K5yczeNbNHzeyIsG0AsC5qncKw7TBmNsPMlpvZ8i1btiQ6VhGRZitlE42ZzTez96v4mQY8CAwFxgIbgP+LtX93n+vu2e6e3atXrzhHLyIiFVL2GY27n1mX9czs18Dfwo9FwKCoxQPDNhERSZKUvaKpiZn1i/p4EfB++P4FYLqZtTWzIcBwYFljxyciIoek7BVNLX5iZmMBBwqA6wDcfZWZPQ18AJQCMzXiTEQkudIy0bj75TUsuxe4txHDERGRGqTlrTMREUkfSjQiIpJQSjQiIpJQSjQiIpJQSjQiIpJQSjQiIpJQSjQiIpJQSjQiIpJQSjQiIpJQSjQiIpJQSjQiIpJQSjQiIpJQSjQiIpJQSjQiIpJQSjQiIpJQKZtozOxSM1tlZuVmll1p2W1mFjGzj83s7Kj28Wb2XrjsfjOzxo9cRESipWyiISjP/AVgUXSjmY0EpgOjgKnAHDNrGS5+ELiWoITz8HC5iIgkUcpW2HT3DwGquCiZBjzl7sVAvplFgBPMrADo4u5Lw+3mARcCLzVa0CJSrfU7DvDwK7mYQYvw77WZYUALC98bGNHvLVw/eB/+R4uK5VHrVPRV0X/Fe7BD/UctJ2pfLaL74fBYWtihvir3f2j/wfsWLerQV7gcoEWLz/bVwghjq2ZfLaL6jfpdcNixR22X5Js7KZtoajAAWBr1uTBsKwnfV24XkSTL6N6BZfnb+d+XPkp2KM1WRbL7T5LGyOrdib/ffGrCE1FSE42ZzQf6VrHodnd/PsH7ngHMAMjIyEjkrkSavZ9eMoa7p43CHRxwd8odcHAcdyh3D5cdavPKyz3o77PrR73+p/9D61fVlxPdX0X/0X1Veo9/Zjmf6T/qeKL74rPLP3ssVfRVaV9EtVXeV/XHciiWzx570NfKdTtYvHorAJ3btmTa2P6NcrVTbaIxs7w49P9Ld7+/uoXufmY9+iwCBkV9Hhi2FYXvK7dXt++5wFyA7Oxsr0ccIlJHZkaHNul4A6VpWJa/nVk5q1m8eivdOrTm66cM4YqTM+navnWj7L+mM58J7AR21LPvDKBbPbetyQvAE2b2c6A/wUP/Ze5eZma7zGwC8AZwBTArAfsXEUl57s6SyDbuz1nNsvzt9OzUhu+ecxRfnTCYTm0bN+nXtrdfuPvd9enYzMrrs13U9hcRJIpewN/NbKW7n+3uq8zsaeADoBSY6e5l4WY3AL8D2hMMAtBAABFpVtydBR9vZlZOhLfX7qBPl7bccf5ILjshg/ZtWtbeQQKk7LWsuz8HPFfNsnuBe6toXw6MTnBoIiIpp7zc+dcHG5mVE2HV+l0M6Naeey4czaXZA2nbKjkJpkJNieZ4YH0D+m7o9iIiUouycufv723ggZwIH2/aTWaPDvzkkjFcdNwAWrdMja9KVpto3H1FQzpu6PYiIlK9krJynl+5njkLIuRt3cvw3p341fSxnHdMP1qlSIKpkLK3zkRE5HDFpWU8u6KIB1+JsG77fo7u14U5XxnH1FF9adEiNWfdqleiCecQ6wtUOTbO3dc2JCgREfmsAyVlPLVsLQ8vymPDzgMcO7Ard54/ijOO7p30b/7XJqZEY2aXAt8leOBe3bYea78iIlK1vcWlPP7GGuYuymfrnmKOzzyCH188htOG90z5BFOhzgnBzGYC9xMMKV5C8GXI0gTFJSLSrO06UMIfXl/Dbxbn8em+Ek7J6sHsKccxYWiPZIcWs1iuPP4fsBk42d3zExSPiEiztmPfQR5dUsDvluSz60Apk4/sxY1ThjN+8BHJDq3eYkk0A4BfK8mIiMTf1j3F/GZxPn94vYC9B8s4e1QfbpoynNEDuiY7tAaLJdGsA9omKhARkeZo064DPPxKHk8sW0NxaTnnj+nPzMnDOKpvl2SHFjexJJrfA9eZWWd3352ogEREmoPCT/fx0Cu5PP1mIWXuXDh2ADdMHsawXp2SHVrcxZJofgSMA+ab2a3AW0o4IiKxKdi6lzkLI/z5rSLM4JLxA/nGxCwyenRIdmgJU+dEE86O/ADwDJAD1VZtc3fX8GYRkSiRzbt5YEEuz68solXLFnzlxAyumziM/t3aJzu0hItlePM04E9ASyCfYB4zDW8WEanBB+t38cCCCC++v4F2rVry9VOHcO1pQ+ndpV2yQ2s0sVx53AXsA85z91cTE46ISNPwzrodzMqJMP/DTXRq24obJg3j6lOG0KNT8xtTFUuiORKYpyQjIlK95QXbuT8nwqJPttC1fWv+35kjuPLkTLp2aJxqlqkolkSzFTiYqEBERNKVu/N6blDNcmnedrp3bMOtU4/k8gmD6dyu+SaYCrEkmmeBqWbW2t1LEhVQhXBetbuAo4ETwqJmmFkm8CHwcbjqUne/Plw2nkMVNl8EbnF3T3SsItI8uTsLP9nC7JwIK9Z8Su/Obfn+eUfz5RMz6NBGY6IqxPKb+D5wIvCMmf2XuxckJqT/eB/4AvBwFcty3X1sFe0PAtcCbxAkmqmonLOIxFl5ufPvDzcxOyfCe0U76d+1HT+cNopLswfRrnVyq1mmolgSzXsEZQFOBC4wsx3AzirWc3cf1tDA3P1DqHYI9WHMrB/Qxd2Xhp/nAReiRCMicVJW7rz0/gZm50T4aONuMrp34McXH8NFxw2kTavUKjaWSmJJNC0IhjNH15qpKgs0xrzVQ8xsJUGi+767LyaYi60wap3CsK1KZjYDmAGQkZGRwFBFJN2VlpXzwjvreWBBhNwtexnWqyO/+NKxXDCmf8pVs0xFsXxhMzPeOzez+QQF1Cq73d2fr2azDUCGu28Ln8n8xcxGxbpvd58LzAXIzs7WcxwROczB0nL+/FYhcxbmsnb7Po7q25kHvjyOqaP70jJFq1mmoqQ+rXL3M+uxTTFQHL5fYWa5wAiC+jgDo1YdGLaJiMTkQEkZTy9fx0MLc1m/8wBjBnbl++eN58yj+6RsueRUlnbDIsysF7A9nBJnKDAcyHP37Wa2y8wmEAwGuAKYlcxYRSS97DtYyhNvBOWSt+wuZvzgI7jvC8cwcUSvtKlmmYqqTTRm9kXgfXf/oD4dx2H7iwgSRS/g72a20t3PBk4H7jazEqAcuN7dt4eb3cCh4c0voYEAIlIHuw+UMO/1NTzyaj7b9x7kpKE9+NX0sZw0tIcSTBxYdV8zMbNy4C53v7teHTdw+8aUnZ3ty5cvT3YYItLIdu4r4bev5fPbJQXs3F/CxBG9uGlKFtmZ3ZMdWlowsxXunl3berXdOtNDchFpcrbtKeaRV/OZ9/oa9hSX8rmRfbhpShZjBnZLdmhNUm2J5i4zu6uefStJiUhK2bzrAHMX5fH4G2s5UFrGucf048bJWRzdr+lUs0xFNSWaRTQ8WRQ0cHsRkQZbv2M/D72Sy1NvrqOs3Jl2bH9umDyMrN6dkx1as1BtonH3SY0Yh4hI3K3dto85CyM8+1Yh7nDxuIHcMHkYg3t0THZozUraDW8WEalNZPMe5iyM8PzK9bRsYUw/PoPrJw1jQDOoZpmKlGhEpMn4aOMuZudE+Pt7G2jbqgVXnpzJjNOH0qcZVbNMRUo0IpL23ivcyayc1fzrg010bNOS6ycO4+unDqFnM6xmmYqUaEQkba1Y8ymzclaz8OMtdGnXilvOGM5Vp2TSrUObZIcmUZRoRCStuDtL87YzK2c1r+Vuo3vHNnz77CO5/KTBdFE1y5SkRCMiacHdWbR6K7NzVvNmwaf07NSW2889mq9MUDXLVKezIyIpzd2Z/+FmZues5p3CnfTr2o4ffH4UXzpe1SzThRKNiKSk8nLnpfc3MntBhA837GJQ9/b87xeO4eJxqmaZbmJONGbWG8gGjgCq/N8Jd5/XwLhEpJkqLSvnb+9uYPaCCJHNexjaqyP/d+mxfH5sf1qrmmVaqnOiMbPWwEMEdV6qO9tGMG2NEo2IxORgaTl/ebuIOQsjFGzbx5F9OjPrsuM495h+qmaZ5mK5ovkhcBWQCzwOrANKExGUiDQfB0rKeGZFIQ8tzKVox35GD+jCw5eP53OqZtlkxJJovgx8Ahzn7vsTFI+INBP7D5bxxLK1zF2Uy6ZdxYzL6MY9F41mkqpZNjmxJJrewJzGSjJm9lPgAuAgwVXUVe6+I1x2G/B1oAy42d3/GbaP51CFzReBW7y6ym4ikhR7ikv5w+tr+M3iPLbtPciEod35xRfHctIwVbNsqmJJNGuBxiza8G/gNncvNbMfA7cB3zGzkcB0YBTQH5hvZiPcvQx4ELgWeIMg0UxF5ZxFUsLO/SX8/rUCHl2Sz459JZweVrM8XtUsm7xYEs3vgJlm1tXddyYonv9w939FfVwKXBK+nwY85e7FQL6ZRYATzKwA6OLuSwHMbB5wIUo0Ikm1fe9BHn01n9+/VsDu4lLOPLoPN07JYuwgVbNsLmJJND8CjiW4grgVWOHuuxIT1mGuBv4Yvh9AkHgqFIZtJeH7yu1VMrMZwAyAjIyMeMYqIsDm3Qf4zeJ8Hlu6hv0lZZwzui8zJ2cxqn/XZIcmjazaRGNm5VRdYdOA+eE6VW3q7l6nBGZm84G+VSy63d2fD9e5nWB02+N16bOu3H0uMBcgOztbz3FE4mTDzv08/EoeTy5bS0lZOZ8/tj8zJ2cxvI+qWTZXiS7lXCN3P7Om5WZ2JXA+cEbUQ/0iYFDUagPDtqLwfeV2EWkE67bvY87CXP60Yh3u8IVxA7hhUhaZPVXNsrlL2VLOZjYVuBWY6O77oha9ADxhZj8nGAwwHFjm7mVmtsvMJhAMBrgCmNXYcYs0N3lb9jBnYS7PvV1ESzO+dPwgrjt9GIO6d0h2aJIiUnmus9lAW+Df4S26pe5+vbuvMrOngQ8IbqnNDEecAdzAoeHNL6GBACIJ8/HG3TywIMLf3l1Pm1Yt+NpJQTXLvl1VzVI+K5YpaPKAX7r7/TWsMxP4lrsPbWhg7p5Vw7J7gXuraF8OjG7ovkWkeu8X7WR2ToR/rNpIxzYtufb0oVxz6lB6dVY1S6laLFc0mUBt4xG7AYPrHY2IpKy3137KrJwIOR9tpnO7Vtw8JYurThnCER1VzVJqFu9bZ50JvskvIk3EG3nbmJUT4dXIVo7o0Jr/PmsEV5ycqWqWUmc1Jhozq/wFk25VtEFQLiADuBjIi1NsIpIk7s6rka3MejnCsoLt9OzUlu+dexRfOXEwHdum8qNdSUW1/Ykp4LNDnG8Jf6pjwDcbGJOIJIm7k/PRZmblRFi5bgd9u7TjrgtGMv2EDFWzlHqrLdHMI0g0RjBc+F1gZRXrlQHbgJcrTR0jImmgvNz556qNzMqJ8MGGXQw8oj33XXQMF48fQNtWSjDSMDUmGne/suK9mV0BPOfudyc6KBFpHGXlzt/eXc8DCyJ8smkPQ3p25KeXjOHC4waomqXETZ1vtrq7/tSJNBElZRXVLHPJ37qXEX068avpYzl/TH9Vs5S4q9dTPTMbBBwHdAV2Am+7+7p4BiYi8VdcWsafVhTy4MJcCj/dz6j+XXjoq+M4a2RfVbOUhIkp0ZjZcGAOMKWKZTkE39L/JE6xiUicHCgp48lla3n4lTw27jrA2EHduHvaKCYf2VvFxiThYpkZIAt4DehBUPHyVWAjwezLpwJnAK+a2cnuHklArCISo73FpTy2dA2/XpzP1j3FnDCkOz+79FhOyVI1S2k8sVzR/C9BkrkFeMDdyysWmFkL4CbgF8B9wBfjGaSIxGbXgRLmvVbAI6/m8+m+Ek4b3pMbJx/HiUN7JDs0aYZiSTRnAC+6+2EzIodJ51dmdjZQ49T/IpI4n+49yG+X5PPb1wrYfaCUM47qzcwpWYzLOCLZoUkzFkuiaUPV36GJ9jZwWv3DEZH62LqnmF8vzuOx19ew92AZU0f15cYpWYweoGqWknyxJJp3gGpnVA5lEXypU0QawcadB3h4US5PLlvLwdJyzh/TnxunZDFC1SwlhcSSaO4DnjOzc9z9sDovZnYecBFwYbyCE5GqFX66jwcX5vLM8kLK3LnouAHcMGkYQ3t1SnZoIoeJJdH0ICgk9jcze5mg1PMmoA8wkWDI81+BnuEsAv/h7vPiE65I81awdS8PLIjw3NtFmMGl2YP4xkRVs5TUZu5e+1qAmZVzaN6zmkR3aIC7e8yTJZnZT4ELCMoO5AJXufsOM8sEPgQ+Dldd6u7Xh9uM51CFzReBW7wOB5idne3Lly+PNUSRRrN6U1DN8oV31tO6ZQsuOyGD6yYOpV/X9skOTZoxM1vh7tm1rRfLFc1VDYinPv4N3ObupWb2Y+A24Dvhslx3H1vFNg8C1wJvECSaqaics6SxVet38sCCCC+9v5H2rVtyzWlDuea0IfTurHLJkj5imevs94kMpIr9Rc8CvRS4pKb1zawf0MXdl4af5xE8L1KikbSzct0OZuesZv6Hm+ncthUzJ2Vx9alD6K5qlpKG0qWC0dXAH6M+DzGzlQTzrH3f3RcDA4DCqHUKw7YqmdkMYAZARkZVtdxEGt+bBdu5/+XVLF69lW4dWvPNz43gaydn0rW9qllK+oo50ZhZL4JKmkcDHd39mqj2IcB77r6/jn3NJ5jCprLb3f35cJ3bgVLg8XDZBiDD3beFz2T+YmajYj0Od58LzIXgGU2s24vEi7vzWu427n95NW/kb6dnpzZ895yj+OqEwXRSNUtpAmKdVPPrwP1AO8IH/cA14eI+wOsEVwmP1KU/d69xFgEzuxI4Hzij4qG+uxcDxeH7FWaWC4wAioCBUZsPDNtEUpK7s/DjLczKWc1ba3fQp0tb7jh/JJedkEH7Nio2Jk1HLJNqfo7gCuBd4E7gbOD6iuXu/r6ZrSJ4LlKnRFPL/qYCtwIT3X1fVHsvYLu7l5nZUGA4kOfu281sl5lNIBgMcAVw2HQ5IslWXu7864NNzF6wmveLdjGgW3vuuXA0l4wfqHLJ0iTFckXzHYLbVhPdfZeZHVfFOu8CJ8UlMpgNtAX+Hc4yWzGM+XTgbjMrAcqB6919e7jNDRwa3vwSGgggKaSs3HnxvQ3Mzonw8abdDO7RgZ9cPIaLxqmapTRtsSSabOApd99VwzqFVP3MJWbuXuV0N+7+LPBsNcuWA6PjsX+ReCktK+f5let5YGGEvC17yerdiV9+aSznj+lHKyUYaQZinVRzby3rdAPK6h+OSNNxsLScZ98qZM7CCOu27+fofl2Y85VxTB2lapbSvMSSaAqA8bWscyKHvrEv0iwdKCnjj2+u46FXctmw8wDHDuzKneeP4oyjVc1SmqdYEs3zwK1mdqm7P1N5oZldBYwBbo9XcCLpZN/BUh5fupa5i/PYsruY4zOP4McXj+G04T2VYKRZiyXR/ASYDjxpZpcAXQHM7EaCGjRfAFajkV7SzOw+UMK819fwyKv5bN97kFOyenD/9OOYMLS7EowIsU1B86mZTQTmAZdGLbo/fF0MfNnda3uOI9Ik7Nh3kN8uKeC3S/LZdaCUyUf24sYpwxk/WNUsRaLF9IVNd18LTDKzMQTDmHsQTAOz1N1XJCA+kZSzdU8xj7yazx9eX8Oe4lLOGtmHm6YM55iBqmYpUpV6zW/h7u+iSprSzGzedYCHF+Xx+BtrKC4t57xj+nHjlCyO6tsl2aGJpLT6zHU2GOhFMP3MlvAqR6TJKtqxn4cW5vLH5esoK3emje3PzMlZDFM1S5E6qVOiMbOewPeAy4DelZZtIpjw8n+jvqEvkvbWbNvLnAW5PPtWIWZwyfiBfGNiFhk9VM1SJBa1JhozG05QhGwQwUSapcC28H13gpkAvglcbGZnunte4sIVSbzI5j3MWRDh+XfW07KF8ZUTM5gxcRgDuqmapUh91JhozKwFwdVKBrAQuAd41d0PhsvbEgxtvh2YCDwGnJzAeEUS5sMNu5i9IMKL722gXauWXH1KJteeNpTeXVTNUqQharuiOYtgjrOngcsqpuqvEE7ZP9/MXiYoTHaxmX3O3f+dkGhFEuDdwh1O5siCAAASgUlEQVTMyonw7w820altK74xcRhfP3UIPTq1TXZoIk1CbYnmYoLaLzdVTjLR3N3DL25+nqDkshKNpLwVa7Zz/8sRXvlkC13ateK/zhzOVScPoWsHVbMUiafaEs04YIm7b6mtI3ffbGavhtuIpCR35/W8bcx6OcLredvo3rENt049kssnDKZzOyUYkUSoLdEMAl6Nob9VBCPTRFKKu/PKJ1uYnRNh+ZpP6d25Ld8/72i+fGIGHdqoXLJIItX2N6wLsCOG/nYAnesfziFm9kNgGkFxs83Ale6+Plx2G/B1gpIEN7v7P8P28RwqfPYicEtNt/yk6XN3/v3BJmYviPBu4U76d23HD6eN4tLsQapmKdJIaks0bYitvkx5uE08/NTd/wfAzG4G7gCuN7ORBJN7jgL6EwxGGOHuZcCDwLUEpZxfBKaiKpvNUlm584/3NzIrZzUfbdxNRvcO/OgLx/CFcQNp00rFxkQaU13uGSTliqBSJc+OUXFMI6j0WQzkm1kEOMHMCoAu7r4UwMzmAReiRNOslJaV89d31zM7J0Lulr0M7dWRn3/xWD5/bH9VsxRJkrokmrvM7K5EB1IVM7sXuIJg4s7JYfMAYGnUaoVhW0n4vnK7NAMHS8t57u1C5izMZc22fRzVtzOzv3wc54zuR0tVsxRJqrokmlj/ltb5CsjM5hPMLFDZ7e7+vLvfDtwePpO5Ebgzxlhq2vcMYAZARkZGvLqVRnagpIxnlq/joVfyKNqxn2MGdGXu5eM58+g+KpcskiJqTDTuntB7De5+Zh1XfZzgmcudQBHBaLgKA8O2ovB95fbq9j0XmAuQnZ2tAQNpZv/BMh5/Yw1zF+WxeXcx4wcfwb0XjWbiiF4qNiaSYlJ2XKeZDXf31eHHacBH4fsXgCfM7OcEgwGGA8vcvczMdpnZBILBAFegap9Nzp7iUua9XsAji/PZtvcgJw3twS+/NJaThvVQghFJUSmbaIAfmdmRBCPZ1gDXA7j7KjN7GviAYILPmeGIM4AbODS8+SU0EKDJ2LmvhN+9VsCjS/LZub+EiSN6cdOULLIzuyc7NBGphelrJsGts+XLlyc7DKnC9r0HeeTVPOa9tobdxaV8bmQfbpycxbGDuiU7NJFmz8xWuHt2beul8hWNNGObdx/g14vyeGzpWg6UlnHu6H7MnJzFyP6qZimSbpRoJKWs37Gfh1/J5ck311FaVs60sQOYOXkYWb3jMuGEiCSBEo2khLXb9vHgKxH+tKIQd7h43EC+MWkYmT07Jjs0EWkgJRpJqtwte5izIJe/rCyipRnTj8/guolDGXiEyiWLNBVKNJIUH2/czewFEf727nratmrBlSdnMuP0ofRRNUuRJkeJRhrV+0U7mZWzmn+u2kTHNi257vRhXHPaEHqqmqVIk6VEI43irbWfMuvl1Sz4eAud27Xi5jOGc/UpmXTrEK/JvkUkVSnRSEItzdvGrJzVLIls44gOrfn22Udy+UmD6aJqliLNhhKNxJ27s3j1VmbnRFhWsJ2endpy+7lBNcuObfVHTqS50d96iRt35+UPNzNrQYR31u2gX9d2/ODzo/jS8apmKdKcKdFIg5WXO/9YtZFZORE+3LCLgUe0576LjuHi8QNo20oJRqS5U6KReistK+fv721gdk6E1Zv3MLRnR3526bFMG9uf1qpmKSIhJRqJWUlZOc+9XcScBREKtu1jRJ9O3H/ZcZx3jKpZisjhlGikzopLy3hmeSEPLsylaMd+RvXvwkNfHc9ZI1XNUkSqp0Qjtdp/sIwnl63l4UW5bNpVzHEZ3bjnwtFMOlLVLEWkdko0Uq09xaU8tnQNv1mcx9Y9BzlxSHd+/sWxnKxqliISg5RNNGb2Q4ISzuXAZuBKd19vZpnAh8DH4apL3f36cJvxHKqw+SJwi6uyW8x27i9h3msFPLIknx37SjhteE9umjKcE4aomqWIxC5lEw3wU3f/HwAzuxm4g7CcM5Dr7mOr2OZB4FrgDYJEMxWVc66zT/ce5NEl+fxuSQG7i0s58+jezJycxXEZRyQ7NBFJYymbaNx9V9THjkCNVyZm1g/o4u5Lw8/zgAtRoqnVlt3F/GZxHn9YuoZ9B8s4Z3RfbpySxaj+XZMdmog0ASmbaADM7F7gCmAnMDlq0RAzWxm2f9/dFwMDgMKodQrDtur6ngHMAMjIyIhz5Olh484DPPRKLk8uW0tJWTkXHNufmZOzGNFH1SxFJH6SmmjMbD7Qt4pFt7v78+5+O3C7md0G3AjcCWwAMtx9W/hM5i9mNirWfbv7XGAuQHZ2drN6jrNu+z4efCWXPy0vpNydi44bwA2TsxiiapYikgBJTTTufmYdV32c4JnLne5eDBSH268ws1xgBFAEDIzaZmDYJqH8rXuZsyDCc28X0cKMS7MHcv3EYQzqrmqWIpI4KXvrzMyGu/vq8OM04KOwvRew3d3LzGwoMBzIc/ftZrbLzCYQDAa4ApiVjNhTzSebdvPAggh/fWc9rVu24PKTBjPj9KH069o+2aGJSDOQsokG+JGZHUkwvHkNh0acnQ7cbWYl4bLr3X17uOwGDg1vfolmPhDg/aKdPLAgwkvvb6RDm5Zce9pQrjltKL06q5qliDSelE007n5xNe3PAs9Ws2w5MDqRcaWDlet2MOvl1bz80WY6t23FTVOyuPqUIRzRUdUsRaTxpWyikdgty9/OrJzVLF69lW4dWvOtz43gipMz6dpe1SxFJHmUaNKcu7Mkso37c1azLH87PTu14bZzjuKrEwarmqWIpAT9S5Sm3J0FH29mVk6Et9fuoG+Xdtx5wUimH59B+zYqNiYiqUOJJs2Ulzv/+iCoZrlq/S4GdGvPvReN5pLxA1XNUkRSkhJNmigrd/7+3gYeyInw8abdZPbowE8uGcNFxw1QNUsRSWlKNCmupKyc51euZ86CCHlb9zK8dyd+NX0s5x3Tj1ZKMCKSBpRoUlRxaRnPrijiwVcirNu+n5H9uvDgV8Zx9qi+qmYpImlFiSbFHCgp46lla3l4UR4bdh7g2EHduOuCUUw5qreKjYlIWlKiSRF7i0t5/I01zF2Uz9Y9xZyQ2Z2fXDKGU7N6KsGISFpTokmyXQdK+MPrQbnkT/eVcGpWT26achwnDu2R7NBEROJCiSZJduw7yKNLCvjdknx2HShlylG9uXFKFuNUzVJEmhglmka2dU8xv16cx2Ovr2HvwTKmjgqqWY4eoGqWItI0KdE0ko07DzB3UR5PLFvDwdJyzh8TVLM8sq+qWYpI06ZEk2CFn+7joVdyefrNQsoqqllOGsbQXp2SHZqISKNQokmQgq17mbMwwp/fKsIMLhk/iBsmqZqliDQ/SjRxFtm8m9k5EV4Iq1l+dcJgrpuoapYi0nylfKIxs28BPwN6ufvWsO024OtAGXCzu/8zbB/PoQqbLwK3uLs3RpwfrN/F7AWreen9jbRv3ZJrThvKNacNoXfndo2xexGRlJXSicbMBgFnAWuj2kYC04FRQH9gvpmNcPcy4EHgWuANgkQzlQSXc35n3Q5m5USY/+EmOrdtxcxJWVx96hC6q5qliAiQ4okG+AVwK/B8VNs04Cl3LwbyzSwCnGBmBUAXd18KYGbzgAtJUKIpK3eunbecnI8207V9a775uRF8TdUsRUQOk7KJxsymAUXu/k6lKVgGAEujPheGbSXh+8rt1fU/A5gBkJGREXN8LVsYmT068p2pR3H5SYPppGqWIiJVSuq/jmY2H+hbxaLbge8R3DZLCHefC8wFyM7OrtdznDsuGBnXmEREmqKkJhp3P7OqdjM7BhgCVFzNDATeMrMTgCJgUNTqA8O2ovB95XYREUmilKyc5e7vuXtvd89090yC22Dj3H0j8AIw3czamtkQYDiwzN03ALvMbIIF2ekKPvtsR0REkiDtHiy4+yozexr4ACgFZoYjzgBu4NDw5pdI8IgzERGpXVokmvCqJvrzvcC9Vay3HBjdSGGJiEgdpOStMxERaTqUaEREJKGUaEREJKGUaEREJKGskeacTGlmtgVYU49NewJb4xxOsuhYUpOOJTXpWAKD3b1XbSsp0TSAmS139+xkxxEPOpbUpGNJTTqW2OjWmYiIJJQSjYiIJJQSTcPMTXYAcaRjSU06ltSkY4mBntGIiEhC6YpGREQSSolGREQSSommCmb2qJltNrP3o9p+aGbvmtlKM/uXmfWvZtupZvaxmUXM7LuNF3XVGngsBWb2Xrje8saLumpVHUvUsm+ZmZtZz2q2TfnzErWstmNJ+fNiZneZWVEY40ozO7eabVP+vMRwLCl/XsL2m8zsIzNbZWY/qWbb+J4Xd9dPpR/gdGAc8H5UW5eo9zcDD1WxXUsgFxgKtAHeAUam47GEywqAnsk+HzUdS9g+CPgnwZduD4s3Xc5LXY4lXc4LcBfw37VslxbnpS7HkkbnZTIwH2gbfu7dGOdFVzRVcPdFwPZKbbuiPnYEqhpFcQIQcfc8dz8IPAVMS1igddCAY0k5VR1L6BfArVR/HGlxXkK1HUvKqeFYapNO5yXtVHMs3wB+5O7F4Tqbq9g07udFiSYGZnavma0DvgLcUcUqA4B1UZ8Lw7aUU4djgeAfu/lmtsLMZjRedHVnZtOAInd/p4bV0uK81PFYIA3OS+im8Bbto2Z2RBXL0+K8hGo7FkiP8zICOM3M3jCzV8zs+CrWift5UaKJgbvf7u6DgMeBG5MdT0PU8VhOdfexwDnATDM7vdECrAMz6wB8j+oTZdqI8VhS+ryEHiS49TIW2AD8X3LDaZC6Hks6nJdWQHdgAvBt4Gkzs0TvVImmfh4HLq6ivYjgHnuFgWFbKqvuWHD3ovB1M/AcwSV1KhkGDAHeMbMCgt/3W2bWt9J66XBe6nos6XBecPdN7l7m7uXAr6k6xnQ4L3U9lrQ4LwRXJ3/2wDKgnGBSzWhxPy9KNHVkZsOjPk4DPqpitTeB4WY2xMzaANOBFxojvljU5VjMrKOZda54D5wFHDZCKpnc/T137+3umR6U+y4Exrn7xkqrpvx5qeuxpMN5ATCzflEfL6LqGFP+vEDdjiVdzgvwF4IBAZjZCIKH/ZVnbo7/eUn2yIhU/AGeJLhELiH4C/914FmCPzjvAn8FBoTr9gdejNr2XOATglEbt6frsRDcKngn/FmVqsdSaXkB4aifdDwvdTmWdDkvwB+A98I/Yy8A/dL1vNTlWNLovLQBHgv//r8FTGmM86IpaEREJKF060xERBJKiUZERBJKiUZERBJKiUZERBJKiUZERBJKiUZERBJKiUZERBJKiUYkzsysg5ndEdb8OGBm68zsPjNrXYdtJ4W1aCp+qpqBoj4x9azUr75AJ42mVbIDEGlKwulK5gPDCea7eh44H7gN6AFcV8euXgEWcvj0INH7OgqYSTClyCCgfbj+28Cfgcc8nA4e2Af8IHx/JTC4jnGINJhmBhCJk3BeqNeAo4Cz3X1J2N6JYFqSgQTT/VSeiy26j0nAAuAH7n5XDevdAdxJcFfidWA5sBvoQ1Dwajiwwt2zq9h2ITDR3RM+a68I6IpGJJ7+GxgP3FiRZADcfY+ZPQfcApwGPNOQnZjZ9wiuTtYBl7r7G1WsM5WggJpI0inRiMSBmbUnqO+xAZhbxSrbwtfDpv2PcT+ZBKWFS4Bz3b3KGYLd/R9mtqAh+xKJFw0GEImPi4BuwBPuXlLF8nbh68EG7ucqoDXwbHVJpkLU8xmRpNIVjUh8nBe+DjCzu6pYfmb4uq6KZbE4NXx9uYH9iDQaJRqR+KhIANNrWe+DBu6noghXYQP7EWk0unUm0kBhRcUMYJW7W+UfoAvBM5V17l6QzFhFkkGJRqThBoSv1dVVP4vgucqLcdjXhkr7FEl5SjQiDdcmfK3u4ftV4eujcdjXq+HrGXHoS6RRKNGINFzFFzAPG7psZhMI6q+/5O7L4rCv3xLchrvYzEbWtKKZtY3D/kQaTIlGpIHcfSvwITDezMZUtJvZYOBJYCdwQ5z2VUDwPZo2wN/N7LBv/of7Pgf4Rzz2KdJQGnUmEh/3AI8DL5vZY0BH4IuAA+fFcxCAu99nZq0IpqB508xeI5iCZg+fnYJmebz2KdIQSjQiceDuT4SzM98KfINgcsunCeYsq26QQEP2d7eZPUNwpTSZ4DlQO4IZCFYCPwYei/d+RepDiUYkTtz998DvG3F/HwI3Ndb+ROpLszeLpJCo2ZsrfOzuR8Wh357Alug2zd4sjUVXNCKppYBDdWOghno0MYquRyPSqHRFIyIiCaXhzSIiklBKNCIiklBKNCIiklBKNCIiklBKNCIiklBKNCIiklD/HzygRrJt0mkIAAAAAElFTkSuQmCC\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x7f6676780cc0>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"temp_ic = Dataset('MOM_IC.nc').variables['Temp'][0,:,1,1]\n",
"z = Dataset('MOM_IC.nc').variables['Layer'][:]\n",
"plt.figure()\n",
"plt.plot(temp_ic,-z); plt.title(r'Initial $\\theta$ profile', fontsize=20)\n",
"plt.xlabel(r'$\\theta$ [C]', fontsize=20); plt.ylabel('Depth [m]', fontsize=20)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Potential temperature ($\\theta$) drift"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.text.Text at 0x7f6673889358>"
]
},
"execution_count": 3,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAA8gAAAGTCAYAAAD0hDsfAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAIABJREFUeJzs3XmcXGWZ9//PN52FkIUkJISQhYDsoASJuIAaFATjgogLqCiKos+I2zjjuD2COig/R0VRlCcCAjMo7oKCgywiBmUJJoYtkJWQlXT2fem+fn+c01KpPtXd1V17fd+8zqu6z7nPqaurQ1dd577v61ZEYGZmZmZmZtbs+lU7ADMzMzMzM7Na4ATZzMzMzMzMDCfIZmZmZmZmZoATZDMzMzMzMzPACbKZmZmZmZkZ4ATZzMzMzMzMDHCCbGZmZmZmZgY4QTYzMzMzMzMDnCA3NUn7SvqSpHmSdkh6VtLXJA0o0P4GSc9JGlLpWHNiOFFSSPpgtWLIiaWo1y89p+qvYS3qy+9V0rT03I5tXgniGZ13zejrNc3MzPpK0gXp+9IFZX6ejvfWS8v5PMWQdGnee/PVJbru9/Oue2kprmv1ywlyk5I0DngY+CLwD+C7wCbgc8D3M9q/BDgfuDwitlYgvk+lf6Telbs/Ih4Bfgt8VdLQcsdRSLGvX3pOTbyGlSRpgqTrJK2QtFPSEknfkTQyt12Jfq9/Br5M4df/KEnfk/SYpI2SdqVx3SbpQkmDcppvS6/1ZeCZXsZjZmZ5JLVI+pCkP0taJ2l3euN4rqRrJL0545wRkr4iaY6kLen7yXJJD0j6lqQT0nb5CVR325KKvwA1QtLk9DW4vtqx9MINJO/Pv886KGlg+r5+m6SV6b+Xzem/n+9IelHeKben17uhzHFbnehf7QCs8iQNBH4HHAycGhH3p/u/CjwOfFDSJRGxKue0y0gSwB9WKMwT08dHMo59HXgQ+DjwtQrF80+9fP2gtl7DspP0AuCvwAHALcA84CTgE8CZkk6OiLU5p/T193pvRFxaIJYvAZeQ3BT8G8mb4GZgLPAq4Brg/wBTASJiG3Bpeu40kt+1mZn1gaQWkqTmTGADcBuwDBgIHAu8CzgKuDXnnIOA+4HJwCLgJqAVGEnyPvdJYDswG7g342mnAGeR3Mz+bd6xDaX4uRrUQ8DRJK91rbk+Iu7NOiDpCJLfc0fsdwJLSf6NHQ1cBHxc0lsi4laAiLgduD19v39f2aO3mucEuTn9G8mbysUdyR1ARGyR9BuSBOaVwC/gn39sTgOuiYjtFYrxRGAL8HT+gYh4KB1G+2FJl0dEe4Vi6lDU6we19xpWyA9IkuOPR8T3OnZK+jbwKZIbBh/p2F+u36ukz5PcGX4WeHtEPJjR5kzgM6V4PjMzK+g8kuT4H8CrI2Jj7kFJ+wIvzTvnKyTJ8XXAByMi8s4ZB4wDSJOme/OOX0CSIM8pdBPVOktvFPd5ylIlSRoL3A1MAL4DfD7/M5ek0SQ3zEd2voJZwkOsm4ykwcC/AyuBGRlNOnr0DszZ9wFAwM8KXHNmN0OY7isivsuVzPc8ChgKtOdc5/ycpjcDk4DTe3rtUujl6we1+RqWTdp7/DpgCXBV3uFLgK3A+eo8F7ukv1dJk0l6gncD07OSY4CI+F/g9aV4TjMzK+gV6eP1+ckxJElZRPypwDnfy0+O03NWRsTfSxznXtLpQt+XtCgdrrtW0q1Kpk7ltrs6fa89q8B1Xpoe/2Xe/nGSrlIyDWmXpDWSfi3pxKzrFLh2SLq3wLHr0+OT0+8vBRanh9+X93njgrRNwTnIkg6XdGM6zL1jutKNkg7PaNsx7H2apLdJekjSNiXD62+WNL6nP2MP/CdJcvzTiPhUVodERLRGxMdIPm+YZXIPcvM5GxgBXBsRuzOO75M+7srZdxrQBjxQ4Jq/Be7K2P9+kmTnniLi+zvJ8Nf3kQzPvTPn2L05X3f03J4O3FHE9fuqN68f1OZrWE6npo9/zO8JjojNku4nSaBfRnK3t0Opf6/vBwYAN0fEY101jIidJXg+MzMrrOMm8hG9PGdOacPpnqQXA38ERpG8L/0aGA28BZgp6ex0iC4k770fBt5LMrUoX8fw3etzrn8IMBM4iOS9/qfARODtwBsknRMRmXNt++Beks8yn6Dz0PMuX+P0psBdwDCSofBPkNyQfw9wlqTTIuLhjFP/BXhzes6fSUYKvBM4XtKUvr4Hpx0YHZ0AX+6uvd/zrStOkJvPG9LH8Vl3BUkSOUiGo5L28E0BnixUWCoivpm/T9I3SRK7H5MMj+qRiPi5pBEkbyI3RERWLy0kBbIgmT9aSUW9flDTr2Huc32S5M2yp+ZERP5crlxHpo+FhnfPJ0mQj2DvBLnUv9dT0se7u2xlZmaV8GvgP4CPSBoG/AZ4JCK6Kob4M5K/5ddImkqSrM7Oq2FRFpL6Az8nGY11akT8OefYQSTvWddKmhwROyPib5KeBt4oaVRErMtpPwg4F3gO+N+cp7maJDn+YkRcltP+B8B9wA2SDo6ILaX6uSLiXiUFyj5BEUPPJQm4ERgOvCcibso59k6SXtn/lnRMxjSpM4GXRMSjOef8hGTY/Vkkr3NfTAUGAcsj4qk+XsuanBPk5tORMJzbTbsn0sfxQAvJkOJupX88ryIpeHQV8LGsIVHdeHH6WHDIVERslLSDJIGspGJfP6jR1zDPJymuENUNdC52kmu/9LHTELq8/Xsl5WX4vY5LH5eV6HpmZtZLETFb0ntIVn54T7ohaR1JMnhdRPwu77SrSBLIT5FMcfr39JwlJD2Z34+If5Qp5DcALwC+mZscpz/LCknfIJnr+lqSSsiQvD9eRpL45U4xehPJvNcrImJP+jNMILlZvBT4Rt71/yrppySv0VtJEtNqewVJb/HfcpNjgIj4maSLST4nnULy+8x1ZW5ynPoRyet0En1PkP1+byXjOchNJO3JnAQ8HhHK30juCO4Gno2IJelp+6eP63tw/RaSYUP/B/hGRFzci8QOkuRuN5D/hzTfOpJhTt3FtaSb+b352/8UuE5vXj+o7dcQgIiYnPUzdbFd0IuYeqpHv1czM6s/EfFzkvfSM4CvklS17kcyZPlWSTekN4o72kdEfJ4kATqXJCG9L/3+g8Ajkj5UpnBfnj4enM6l3WsjSewgqY7c4Uagnc7VkDsNrwZOSB//UmDa1j157aqt4+Z7oWlfXcU7K2Nfx2g7F8yymuIe5ObSUQhheYHjryOZr3l7zr6OAgf7dG7+PEkDgJ8AbwMuiYgeDwnOu05/4IXAEz2YHzI4J76uLAR2FBHGigL7e/P6QW2/huXS0UO8X4HjHfuzltjo6e+1J1aSfHApZREQMzPrgzQZ/GO6ddwcPoekUvV7SYZe/zbvnA0kw61/lp4zBPgs8EXge5JujYjVJQ614wb327tpNzQnzmWS7gZOl3R0RDwp6QCSIcZzImJuznkd74WFRph17C9mClQ59SXerPf7PeljS1+Cyntuv99bnzlBbi4D08dCSdP708frcvY9lz7uTwGS9gF+STIU6d8i4lt9iPEYkkSyy6HBkvqR/AFe3FU7gIh4bR/iydWb1w9q9DXMe/5Sz0HumP9TqBBLR6XLveYoF/N77aGZwGtIhr9dW6JrmplZCUVEG/BzSS8kSXhfQ9fTeEhrevxfJWvXngKcTDLHuZQ6bvaeFemauT10A0mxyfeRJPHvJvnMfUOB6+evfNFhXF67rgSFP9eXKsEuZbylNovk89kESUdERLWWuLQG4AS5uaxKHzv9YZP0MmA68IeIeCjn0EpgDc8XXco/bwhJRcJTgX+JiB/2McYp6ePsbtodSbJsUiUrWvbm9YPafQ1zlXoOcscyHa+T1C+3WEdamOVkYBudq3qX+vf6Y+BzwDlp0ZAnCjWUNMhVLc3Mqmpz+qguW/X9nJ7qeI96Jcn7dE/9GtgEvEfS50kS5T0ko8RydbxPnyKpf8fc5BwdK0L05Ib3epLq13tJe+endG5OW/pYTO9tR7zTChwvJt6Siojtkv6bZNj9l0jntxfi93zriucgN5GIaAWeBE6U9KKO/ZIOJllWYCNJGf7cc4Jkrs9oSYflHpO0H8nwqFcDF3SX2On5dfgu6KJZRy/rpm5+nJelj/nrJZZNb16/9LxafQ1zYyzpHOSIWJj+XJOBj+Yd/jIwBPjvjKreJf29pnPBLyXp/b8trYDaiaTXs3dVUTMzKzFJ50k6PR0tlH/sQKBjLvF9Ofv/XdKxBa53CklStgf4WxlCvoVkmtZHJU0vEMPLJe2buy+S9Xd/TjLc91PA8cDtEfFcXrtlJEsxTia5UZ173ZcC7yJJfH/Tg1gfAiZJel3e/i+SfQN8PUmvczFFMe8nGSF2iqS35cX7NpIbCU+TjN6qhi+SFOl6t6T/UrL0014kjZH0fbovtmpNzD3Izec/gZuAu9NiVEOAd5D8kXxDXnGpDr8imRt0BrAgZ/9PSCoaPgQcWmDZo6/n3KHreEPMv0Oa65H08TJJxwFbSYpi/SKv3etI7n5mrTNYTr15/aA2X8Ny+xeSdZivlPRakpsLLyX5MPM08IWMc0r+e42Ir6Xzsi8BHpb0V5KhWFuAsSRLSh1OdgERMzMrnZeSLC20StJMnp9OcwjJFKPBJH//f5lzzruBb0iaR9Kju5LkvfdYkqHYAj4dEYXqh/RaROyW9FaS9Y9vS98/5pCMgJoIvAQ4lGRo8ba8028g6c38es73WT5Cknj+V5rczuL5dZDbgfdHxOYC5+b6JslnjFsk/Yyk4OUrSF7be8nr9Y2ILZIeBF4p6SaS9+U24Na8edK554Sk95Ek9T+TdAswj2T011tIevPfm7HEU0VExOr088ZvgX8D3ifpTpIq4QNJapJMI1kO6o5qxGh1IiK8NdlGMtTncZLCVcuAGcD4LtoPBFYDD+bs60fyhzC62FbnXWc2Sa/myG7iu5jkDuWO9DqX5R3fj6SI02/r4fWrxdewgq/VRJJhziuBXcAzJBVIO8Xf298ryZtdAJd20+5o4HvAY+lruCuN6w/AhcCgAufdSzoQwJs3b9689X5L3xM+StIj+lTe3+LbSYbF9ss75wSSnsF7SBLq7el720KSG9andPOcF6TvEdf3Ie4DgMvT949tJDdY55Mk8u8B+hc4b3763GuBgV1cfzzww/Q9chfQSpLkvaSLn+eCjGNvJkmwd6TPeTNJ7/H16TmT89ofBvwubduee92u3ltJEuL/Tn9vu9PH/wGOzGh7aXqdaRnHJhfzu+nqWnntBqbv67fz/OePzSQre1wJvLDAeQV/Zm/NtSmiNyvIWLOR9Dnga8CLI6KYua0d548g+QP8rYj4TB9j+RjJH7hXRkS1hvEUrZZew1rU299rWqDlT8CXI+LSMsR1L/DqSJbyMjMzsypIR9ldApwaEfeW4frTKOPnCasfnoNsPXUFyRCVXi09RDIvZTfw7b4Ekc4n+Rzwq3pKjlM18RrWohL9Xi9J52fPK0E8ozvWxSaZH25mZma14U/pe/TVpbiYpO+n7/cVq2tjta0h5yBLOhP4Lkllvmsi4vIqh1T3ImKHpPOBUyUNic7Flbo7/3d0sw5wD00mGdJ8fQmuVVE19BrWosn0/ve6hKTwV4fWvofDtrxrmpmZWXXdm/d9qWqH3M7enx3yn8eaTMMNsU7L2T9Nsv7cMuBh4LzoYnkXMzMzMzMzs0YcYn0SsCAiFkXELpICBWdVOSYzMzMzMzOrcY2YII8Hns35flm6z8zMzMzMzKyghpyD3BOSLgIuAhgwuOXEkZOH73V8YL+2Tuf0V+d9AAPovL9F2UvAtdB5SHtLgdq4Wbv7Ze7NblvcXjMzeGTuztaIGFPKa55x6pBYuy7772dPPDJ35x0RcWYJQ7I6MHr06Jg8eXK1wzAzswbxyCOP9OgzTiMmyMtJ1tnrMCHdt5eImEFSFIixx4yKc286Y6/jE/dZ1+nCY/tvzHzCA/pv6rRvVEv+evGJ4drVad+wAv34g9T5wD5qyWzbn877WzLOB+jXkAMHzKwUWsbNf6bU12xd18aDd0zo9fkDxi0cXcJwrE5MnjyZWbNKVYPHzMyanaQefcZpxEzpYeBwSYdIGgicC9xa5ZjMzMzMzMysxjVcD3JE7JF0MXAHyTJP10XE41UOy8ysiQVtkT3txMzMzKyWNFyCDBARt5OsaWZmZlUWQHtG/QUzMzOzWtOQCbKZmdWWdtyDbGZmZrWvEecgm5lZk5F0pqSnJC2Q9NmM4++WNFfSo5L+Kun4asRpZmZmtc09yGZmVlZB0BblG2ItqQW4CjgdWAY8LOnWiHgip9li4NURsV7S60lWMXhp2YIyMzOzuuQE2czMyq7Mc5BPAhZExCIASTcDZwH/TJAj4q857R8gWQLQzMzMbC9OkM3MrKwCaCtvgjweeDbn+2V03Tt8IfCHcgZkZmZm9ckJspmZlV0fe5BHS5qV8/2MiJjRmwtJOpUkQT6lLwGZmZlZY3KCbGZmta41IqZ2cXw5MDHn+wnpvr1IehFwDfD6iFhb2hDNzMysEThBNjOzsgooa5Eu4GHgcEmHkCTG5wLvym0gaRLwa+D8iHi6nMGYmZlZ/XKCbGZmZVfOVZAjYo+ki4E7gBbguoh4XNJH0uNXA18C9gd+IAlgTze90mZNp6092L67je270m13su3Y3UZ7BBEkG8nX7Tk3vsp6C8zMmtawQf2ZOnlURZ/TCbKZmZVVEOUu0kVE3A7cnrfv6pyvPwh8sKxBmNWRiGBx61YeXrKOh5es5+El63hm7bZqh2VmtpdjDxrObR9/ZUWf0wmymZmVV0Cbu5fMqi4iuGXOCu54fBUPL1lP65adAIwaMpCpB4/krCnjGTqohcEDWhg8sD+DB7Sw78AWBvXvR79+QoAkJOgngOTrDsp6UjOzPhg8sKXiz+kE2czMzKzBLV27jc/9Zi73L1jL+BGDeeXho3nJ5FGcdMhIXjBmKJLTWzMzcIJsZmZlFpR3DrKZFdbWHlz/1yV8846naOknLjv7OM57yST69XNCbGaWxQmymZmVmWjz4Euzilvw3GY+88u5/H3pBk49cgyXnf1CDhoxuNphmZnVNCfIZmZWVgG0ew6yWcVEBD/880K+c+d89h3UwhXvPJ63TBnvYdRmZj3gBNnMzMysgVxx59Ncec8Cpr/wQL785uMYM2xQtUMyM6sbTpDNzKzsPMTarDJ+8uBSrrxnAe+cOpHLz3mhe43NzIrkBNnMzMoqcIJsVgl3P7maL/72UaYdOYb/PPs4J8dmZr3gBNnMzMquPfxB3ayc5jy7gYt/MptjD9qPq971Yga09Kt2SGZmdckJspmZlZV7kM3Ka0nrVi68/mFGDxvIdRe8hCGD/PHOzKy3fHvRzMzMrE6t3bKTC378EO0R3PD+k1yQy8ysj3yL0czMyioQbb4fa1Zy23e18YEbZrFy4w5+8qGXceiYodUOycys7jlBNjOzsvMcZLPS+/adTzF32Qaufs+JnHjwyGqHY2bWEJwgm5lZWXkOslnpPb16M9fdv4R3Tp3IGcceWO1wzMwahhNkMzMrM9EWHmJtVioRwZdueYyhg/rzmTOPqnY4ZmYNxZ9YzMzMzOrIrf9YwQOL1vGZM49k1JCB1Q7HzKyhuAfZzMzKKoB23481K4nNO3Zz2W1P8qIJ+3HuSyZVOxwzs4bjBNnMzMrOc5DNSuO7d81nzZad/Oi9U2np5/+vzMxKzbf0zcysrCKSOci93axyJF0n6TlJj+Xsu1TScklz0m16xnkTJf1J0hOSHpf0iWLOt555atVmfvzXJZz7kokcP3FEtcMxM2tI7kE2MzOzDtcD3wduzNt/RUR8s4vz9gCfjoi/SxoGPCLpzoh4oofnWzcigv97y2MM26c//36GC3OZmZWLb82bmVnZtaNeb1Y5EXEfsK4X562MiL+nX28GngTGlzi8pnbLnBU8tHgdnznjKBfmMjMro7pMkCX9l6R5kuZK+o2kEen+yZK25wzjurrasZqZNbtkHeR+vd6sJnwsfc+9TtLIrhpKmgycADzYm/Ots807dnPZ7U9y/IT9eOdLJlY7HDOzhlavnzzuBI6LiBcBTwOfyzm2MCKmpNtHqhOemZk9z3OQ69wPgUOBKcBK4FuFGkoaCvwK+GREbOrF+RdJmiVp1po1a0oUfv276k8Lad2yk6+cdZwLc5mZlVldfvKIiD9GxJ702weACdWMx8zMCutY5qm3m1VXRKyOiLaIaAd+BJyU1U7SAJLk+KaI+HWx56dtZ0TE1IiYOmbMmNL+IHVqy8493PTAM0x/4TgX5jIzq4BG+OTxAeAPOd8fkg6v/rOkV1YrKDMzs0YgaVzOt2cDj2W0EXAt8GREfLvY862wX8x6ls079/DBUw6pdihmZk2hZqtYS7oLODDj0Bci4pa0zRdIKmfelB5bCUyKiLWSTgR+K+nYnGFeude/CLgIYNiB+5bjRzAzs1RbeFhoPZD0U2AaMFrSMuASYJqkKSSDAZYAH07bHgRcExHTgZOB84FHJc1JL/f5iLgd+EbW+da9tvbgx/cv4cWTRnDCJE/dNjOrhJpNkCPitK6OS7oAeCPw2oiI9JydwM7060ckLQSOAGZlXH8GMANg7DGjoqTBF6mNqj69mVlZBXKxrToREedl7L62QNsVwPT065mQXXI8Is4vWYBN5q4nV7N03Tb+40wv62RmVik1myB3RdKZwGeAV0fEtpz9Y4B1EdEm6VDgcGBRlcI0M7NUu4ttmRXt2pmLGT9iMGccO7baoZiZNY26TJCB7wODgDuTaU88kFasfhXwFUm7gXbgIxFR9HqOZmZWOh3LPJlZzz22fCMPLV7HF6YfTf8W//9jZlYpdZkgR8RhBfb/iqSCZtNpzxim3VKFOMzMzKzvrp25mCEDW3jnSV732MyskuoyQTYzs/oRyEW6zIqwetMOfvePFbznZQczfJ8B1Q7HzKypOEEuwPPlzMxKx+sZm/XcjX9bQlsE7z95crVDMTNrOk6Q60xbZFe87ldE50w77Z3P94dXMyuTCGjzTUezHtm+q42bHlzK6UeP5eD9h1Q7HDOzpuNPLGZmZmY14tezl7Fh224uPOWQaodiZtaU3INsZmZlJtqzl8g1sxzt7cF1Mxdz3PjhnHTIqGqHY2bWlJwgm5lZWQUeYm3WE3+ev4aFa7ZyxTuPJ13G0szMKswJspmZlZ3XQTbr3o/vX8IBwwbxhhceVO1QzMyalhNkMzMrq0C0e5knsy49t3kHM+ev4eJTD2Ngf99QMjOrFv8FNjMzM6uyPzy6ivaANx3v3mMzs2pyD7KZmZWdh1ibde33c1dw5NhhHD52WLVDMTNrak6Qy6TgcMIKjjJsi87rHQO0qPMH1ay1kc3MSiGAdhfpMito5cbtPLxkPZ8+/Yhqh2Jm1vScIJuZWZmJNi/zZFbQbXNXAvCGF42rciRmZuYE2czMyso9yGZd+/3clRwzbjiHjhla7VDMzJqeP7GYmZmZVcmz67Yx59kNvPF49x6bmdUC9yCbmVnZeYi1WbbbH02GV7/Rax+bmdUEJ8g1rJ2odghmZn0WIQ+xNivg93NXcvyE/Zi0/77VDsXMzPAQazMzq4C26NfrzaxRLWndyqPLN/LGF7n32MysVviTh5mZmVkV3Paoq1ebmdUaD7E2M7OyCqDdc5DNOvndP1Zw4sEjOWjE4GqHYmZmKfcgm5lZmclDrM3yLHhuC/NWbeaN7j02M6sp7kE2M7OyStZBdg+yWa7b5q5EgukvdIJsZlZLnCCbmVnZtXnAktlefj93BSdNHsXY4ftUOxQzM8vhTyxmZmZmFfTUqs3Mf26Lh1ebmdUgJ8hmZlZWgWiP3m89IelMSU9JWiDpsxnHJenK9PhcSS8u+Q9q1kO/n7uCfoIzj3OCbGZWazzE2szMyq69jPdjJbUAVwGnA8uAhyXdGhFP5DR7PXB4ur0U+GH6aFZREcHv567k5S/YnzHDBlU7HDMzy+Me5AprQ502M7NGFgFtoV5vPXASsCAiFkXELuBm4Ky8NmcBN0biAWCEJHffWcUteG4Li1u3ujiXmVmNcg+ymZmVXR+rWI+WNCvn+xkRMSPn+/HAsznfL6Nz73BWm/HAyr4EZlas++a3AjDtyAOqHImZmWVxgmxmZrWuNSKmVjsIs1K4f0Erh44ewvgRg6sdipmZZXCC3ITaor3TvhZ5tL2ZlUdSpKusf2OWAxNzvp+Q7iu2jVlZ7drTzgOL1nLOiydUOxQzMyvAWZGZmZVdVv2Fnm498DBwuKRDJA0EzgVuzWtzK/DetJr1y4CNEeHh1VZRs5euZ9uuNk45fHS1QzEzswLcg9wg2uncK9yPlipEYma2t6DPc5C7vn7EHkkXA3cALcB1EfG4pI+kx68GbgemAwuAbcD7yxaQWQEzF7TS0k+8/AX7VzsUMzMrwAmymZnVvYi4nSQJzt13dc7XAXy00nGZ5frL/FaOn7Afw/cZUO1QzMysgLocYi3pUknLJc1Jt+k5xz4naYGkpySdUc04zcwMSOcg93YzawQbt+1m7rINnHL4mGqHYmZmXajnHuQrIuKbuTskHUMy9+xY4CDgLklHRERbNQIsh3Yic39LH9dTzircZWZWKu1e892a3N8WtdIe8ErPPzYzq2n1nCBnOQu4OSJ2AoslLQBOAv5W3bDMzJpXBLSVcQ6yWT34y/xWhg7qz5SJI6odipmZdaGex659TNJcSddJGpnuGw88m9NmWbqvrNro12kr6vzI3vqqncjczMwqzUOsrdnNXNDKyw4dxYAW/5s2M6tlNftXWtJdkh7L2M4CfggcCkwBVgLf6sX1L5I0S9Ks7et3ljh6MzMzs8Sz67bxzNptnHyYh1ebmdW6mh1iHRGn9aSdpB8Bv0+/XQ5MzDk8Id2Xdf0ZwAyAsceMcreqmVmZBCrrMk9mtW7mglbA84/NzOpBzfYgd0XSuJxvzwYeS7++FThX0iBJhwCHAw9VOj4zM9tbO+r1ZlbvZs5v5cDh+/CCMUOrHYqZmXWjZnuQu/ENSVOAAJYAHwaIiMcl/Rx4AtgDfLSRKlibmdWjAPcgW9Nqaw/uX9jKaUeJyuyfAAAgAElEQVSPRfL/B2Zmta4uE+SIOL+LY5cBl1UwHDMzM7NMj6/YyIZtuz282sysTtRlgmxmZvXF1aitWf1lfjL/2AW6zMzqgxNkMzMrr3CRLmteM+e3cvS44YweOqjaoZiZWQ/4ln4Na4votJmZ1ZvARbqsOW3f1cYjz6z38GozszriHmQzMys79yBbM3pw8Vp2tbV7eLWZWR1xD7KZmZlZGcyc38rAln6cNHlUtUMxM7Mecg+ymZmVlZd5smY1c0ErUyePZPDAlmqHYmZmPeQeZDMzK7v2tFBXbzazerRm807mrdrMKZ5/bGZWV9yDXIT2Brmf0E7nYl/9XAjHzMokcKJrzeevC5PlnU7x/GMzs7riBNnMzMrO1ait2cxeuoF9B7ZwzLjh1Q7FzMyK0BhdomZmZmY1ZPbS9bxw/H70b/FHLTOzeuIe5DJpq3BvSdYayf2KCCFr2LWZWUmEi3TVC0nXAW8EnouI49J9lwIfAtakzT4fEbfnnTcRuBEYS1KXbUZEfDc9Ngr4GTAZWAK8IyLWl/tnqaYdu9t4YuUmLjzl0GqHYmZmRfJtTTMzK6uOKtYu0lUXrgfOzNh/RURMSbfbM47vAT4dEccALwM+KumY9Nhngbsj4nDg7vT7hvb4ik3sbgtOmDSi2qGYmVmRnCCbmVnZOUGuDxFxH7CuF+etjIi/p19vBp4ExqeHzwJuSL++AXhLCUKtabOXJh3kJ0x0gmxmVm+cINeZNiJzy9Je4D8zM7MifUzSXEnXSRrZVUNJk4ETgAfTXWMjYmX69SqSYdgNbc6zGxg/YjAHDN+n2qGYmVmRnCCbmVlZdSzz5B7kuvVD4FBgCrAS+FahhpKGAr8CPhkRm/KPR0RA4aIXki6SNEvSrDVr1hRqVvNmL93AFA+vNjOrS06QK6wNddrMzBpdhHq9WXVFxOqIaIuIduBHwElZ7SQNIEmOb4qIX+ccWi1pXNpmHPBcF881IyKmRsTUMWPGlO6HqKDnNu9g+YbtHl5tZlannCCbmVnZtaNeb1ZdHclt6mzgsYw2Aq4FnoyIb+cdvhV4X/r1+4BbyhFnrZizdAOAC3SZmdUpL/NkZmZlFV7mqW5I+ikwDRgtaRlwCTBN0hSSodFLgA+nbQ8CromI6cDJwPnAo5LmpJfrWA7qcuDnki4EngHeUbmfqPJmP7uBAS3i2IP2q3YoZmbWC06Qa4DLZpmZWS2IiPMydl9boO0KYHr69UzI7u6PiLXAa0sVY62bvXQ9x4wbzj4DWqodipmZ9YKHWJuZWdl5DrI1g7b2YO6yjUzx/GMzs7rlHmQzMyszV6O25vD06s1s29XGCZO6XAnLzMxqmBNkMzMrO/cEWzOY7QJdZmZ1zwmymZmVVeAiXdYcZi9dz6ghA5k0at9qh2JmZr3kBLmGlat4V3vGlft5OrqZmVmfzHl2A1MmjiBZ9crMzOqRsyIzMyuvSJZ66u1mVg82bt/N/Oe2cIILdJmZ1TX3IBuQ3atciHubzaxY7dkrAJk1jLnLkvnHUzz/2MysrjlBbhDtdO5maSnTB9Jikmkzs8BFuqzxzVm6AQmOdw+ymVldc1egmZmZWR/NfnYDh40ZyvB9BlQ7FDMz64OCPciSFpXg+t+JiCtLcB0zM6tbXgfZGltEMHvpek47emy1QzEzsz7qaoj1ZGAjsKGX154EeJyRmZm52JY1tKXrtrF+225OmDSy2qGYmVkfdTcH+YqI+EpvLizJE1WrrK3AJ9IWLz9hZhXmOcjWyGYvTfoSTnCBLjOzuuciXXUmqxgXlK8gl5lZXyXLNflvlDWu2UvXs+/AFo4YO6zaoZiZWR91lSC/BFjRh2v39fyCJP0MODL9dgSwISKmSJoMPAk8lR57ICI+Uo4YzMzMzCAp0PWiCfvR0s83gszM6l3BBDkiHunLhft6fjfXfmfH15K+RTJXusPCiJhSrufO0hYuBm5m1hUX6bJGtWN3G0+s2MSHXnVotUMxM7MSqOsh1pIEvAN4TbVjMTOzwlykyxrV4ys2sqc9OMHrH5uZNYReJchpYnogkLnYX0Qs7UtQRXglsDoi5ufsO0TSHJJe5S9GxF8qFIuZmRXgOcjWqDoKdE1xgS4zs4ZQVIIs6e3AZ4Hjujg3ir1ugee6iyQJz/eFiLgl/fo84Kc5x1YCkyJiraQTgd9KOjYiNmVc/yLgIoBhB+7b13DNzKyAQE6QrWHNXbaRg/bbhwOG7VPtUMzMrAR6nMhK+ihwJbAHuB9Ynn5dFhFxWjfx9AfeCpyYc85OYGf69SOSFgJHALMyrj8DmAEw9phRHvxnZmZmRZu3ahPHHDS82mGYmVmJFNPT+yngOeAVEbG4TPEU4zRgXkQs69ghaQywLiLaJB0KHA4sqkZwhQp3ZRWqaavwusSF1kfO4jWTzawUfBfSGtHOPW0sXLOV048ZW+1QzMysRIpJkMcDP6qR5BjgXPYeXg3wKuArknYD7cBHImJdxSMzM7PneR1ka1ALn9tKW3tw1IHuQTYzaxTFJMjPAoPKFUixIuKCjH2/An5V+WgaVzG9zWZmBflPiTWgeauSEidHHTisypGYmVmpFLOA7w3A6yX5XaBC2ohOWynampmZWd/NW7WZgS39OGT0kGqHYmZmJVJMgnw58DBwl6RXO1E2M7OeilCvN7NaNW/VZg4fO5T+LcV8nDIzs1rW47/oEdEGXAUcBtwDbJDUlrGVrbK1FaedyNzMzCotovdbX0gaJelOSfPTx5EZbSZK+pOkJyQ9LukTfXvW8ijwnlvs9qVq/xyNZN7KTRzp4dVmZg2lmGWezgJ+CbQAi4EVlHGZJzMzawxBVYt0fRa4OyIul/TZ9Pv/yGuzB/h0RPw9HR31iKQ7I+KJSgfbDQHPAEt6ee6rShpNk1u3dRfPbd7J0S7QZWbWUIop0nUpsA14Q0TMLE84ZmbWcAKoXoJ8FjAt/foG4F7yEuSIWAmsTL/eLOlJkpUbai1BBvhxRHylNydKai91MM3snwW6xrkH2cyskRSTIB8J3OjkuLoKVZXu5/WKzcyyjE0TYIBVQJcL1kqaDJwAPFjesKzezVu5GcBDrM3MGkwxCXIrsKtcgZiZWePq41zi0ZJm5Xw/IyJmdHwj6S7gwIzzvrB3DBGSCkYiaSjJUoGfjIhNfYq4PMaQjOSq1vmW46lVm9l/yEDGDK2ZFTDNzKwEikmQfwWcKWlAROwuV0BmZtaA+pYgt0bE1IKXjjit0DFJqyWNi4iVksYBzxVoN4Dkfe6miPh1n6Itk4hYW83zbW/zVm3iqHHDkEdwmZk1lGLWJfgisB74RToEzeqUK1ubWWX1fomnEhT3uhV4X/r1+4BbOkWXZDjXAk9GxLf7+oTW+Nrag6dWb+bIsS7QZWbWaIpJkB8FJgJvAhZKWitpUca2sDyhmplZ3Yo+bH1zOXC6pPnAaen3SDpI0u1pm5OB84HXSJqTbtP7/MxlIGmgpIck3ZX2enfV7h5JD3TVznpn6bpt7Njd7gJdZmYNqJgh1v1IlsJYmrMv69a+xxrVIfcim1kjSocVvzZj/wpgevr1TOrnves9wInA67ua7hQRuyT9F3Ab8G7g+sqE1xzmrUwrWLtAl5lZw+lxghwRk8sYh5mZNaqo6jrIjeatwPyI+GN3DSPiD2nP+dtxglxST67aTD/B4Qc4QTYzazTF9CA3vbYoZkS6mZn9kweplMoJJL3CPXUfaU+5lc5TqzYxefQQBg9sqXYoZmZWYk6QzcysAtyDXCKjgdVFtF8N7F+mWJrWvFWbOfYgF+gyM2tEBbtEJb1D0jG9vXBfzzczM7NOtgPFjOsdCuwoUyxNaevOPSxdt42jDnSCbGbWiLoaM3wz8LY+XLuv5zeN9sje2jK2gtfIWLqprcBmZlZx1ati3WieBQquCZ1hKnsX17Q+enr1ZiJcoMvMrFF1N6nWH01KrA112szMGp4T5FK5F3i5pG6TZEknAq8A/lTuoJrJvFWbAdyDbGbWoLpLkC+V1NabDX+sMTMzSBNd9X6zXN8neUV/IenoQo0kHQX8AmgDflCh2JrCU6s2M2RgCxNGDq52KGZmVgZdFem6j74nuUv6eH7dau/23oOZWfMI3zItiYh4StJXgEuB2ZJ+CdwDLEubjCdZ9/kcYBDwpYh4qhqxNqonV27iyAOH0a+fb96YmTWigglyREyrYBxmZmbWAxHxFUl7gEuAdwHn5TURsBv4QkR8vdLxNbKIYN6qzUx/4bhqh2JmZmXiZZ5KoJK9xe0F9hezEmNWoa4Wz4U2s3JyD3JJRcTXJN0EfAA4GejI2FYCM4EfR8Qz1YqvUa3etJON23dz9DgX6DIza1ROkA3ITprNzErGc4lLLk2AL6l2HM3kyVWbADhyrBNkM7NG5QTZzMzKTr4HZw3gKVewNjNreK4kZWZm5dWXJZ6cWP+TpAMk7Vut8w3mrdzEQfvtw377Dqh2KGZmViZOkM3MzOrDSuDf+nj+p0sUS1Oat2ozRx7o4dVmZo3MCbKZmZVZH9ZA9tzlXH19MVSCazStXXvaWbhmC0eN8/BqM7NG5jnIBbTV2WeItqxFRgv8CP3q7GczswbgodKlcoGkab0817+FPljUuoXdbcFR7kE2M2toTpDNzKz8nJqVyuR0swqbt9IFuszMmkHRCbKkA4CpwEgKLL8bETf2MS4zMzPLERGeFlVF81ZtZkCLOHTMkGqHYmZmZdTjBFnSAOBq4L0Unrsskn4CJ8hmZvY89yBbnZu3ahOHHTCMAS2+T2Fm1siK6UH+KvB+YCFwE/AssKccQTWbQvOd2/2J0swaQeBiW1b3nl61mZMOGVXtMMzMrMyKSZDfBTwNnBAR28sUj3WjrUDSXEzhLSfeZlZp8p8dq2Pbd7WxYuMODh0ztNqhmJlZmRUzTugA4PZKJceS3i7pcUntkqbmHfucpAWSnpJ0Rs7+EyU9mh67UpK7LMzMakH0YTOrsmfWbQXgkNGef2xm1uiKSZCXApUs3fgY8Fbgvtydko4BzgWOBc4EfiCpo1jYD4EPAYen25kVi7aH2ujXaTMzM6sFkq6T9Jykx3L2XSppuaQ56Ta9p+cWc34tW7zGCbKZWbMoJju7Hni9pP3KFMteIuLJiHgq49BZwM0RsTMiFgMLgJMkjQOGR8QDEdFRKOwtlYjVzMysQVxP9s3lKyJiSrrdXuS5PT2/Zi1qdYJsZtYsikmQLwdmAndJOlVStRYCHE9SIKzDsnTf+PTr/P1mZlZlit5vVjkRcR+wrtLn1rrFrVsZO3wQQwYVvTqmmZnVmYIJcjr3t61jA3YDbwNOBO4C1ucez9l6XNla0l2SHsvYzurzT9b9c18kaZakWdvX7yz301VFW0TmZmZWcaHeb/ZPkhZJ+ngVnvpjkuamw6hHVuH8qlrcupXJ+7v32MysGXR1K/Q+ylweJSJO68Vpy4GJOd9PSPctT7/O31/ouWcAMwDGHjPKWaOZWbm42FYpTQZGVPg5f0iy1GOkj98CPlCO8yVdBFwEMGnSpN5HXGJLWrfyumPHVjsMMzOrgIIJckRMq2AcxbgV+ImkbwMHkRTjeigi2iRtkvQy4EHgvcD3qhVkobWNy6U9Y5/Lf5mZWV9FxOqOryX9CPh9uc7PvXk9derUmritsnHbbtZu3eX5x2ZmTaJmcyhJZ0taBrwcuE3SHQAR8Tjwc+AJ4H+Bj0ZEW3ravwDXkBTuWgj8oeKBm5lZZ17mqW6lRTA7nE2yykTFzq+2xWs7CnR5DWQzs2bQ42oTkhYB34mIK7to81Hg0xFxaF8Di4jfAL8pcOwy4LKM/bOA4/r63GZmVloutlUfJP0UmAaMTm9SXwJMkzSF5HbFEuDDaduDgGsiYnqhcyPiWuAbWefXi8WtWwA4ZPS+VY7EzMwqoZhyjJPpft7TCODgXkdje8kapj2gBN0pLtRlZhXnPzul9BZJk4toHxFxYQ8bnpex+9oCbVcA03O+zzqXiDi/J89dqxav2Uo/wcRRTpDNzJpBqdcrGAbsKvE1zczM7HnHA1OKaB9AjxJk62zx2m1MGLkvg/q3VDsUMzOrgC4TZEn5JSRHZOwDaAEmAecAi0oUm2VoK9AL0+KVUMyslrkHuZRuSTergMWtW1ygy8ysiXTXg7yEvT/WfCLdChHwr32MqSbcc1znN8OWodlLPGjfwZ13DhyYfeFBGfsHZP8aon9GDbV+2XXVQhkZcs2WYDOz2vWVkl9R4TnIJTYnIm6odhDNICJYvGYrUw8eVe1QzMysQrpLkG8kSZBFsmzSXGBORrs2YC1wd0T8saQRmplZ/QsPc7H6s2bzTrbuanMPsplZE+kyQY6ICzq+lvRe4DcRUfruhRqkloy5Rln7Cu3vX+ClzZjDFIXGR2f0FhfTtpDM3mYzs3JyD7LVoUWtHUs8OUE2M2sWPS7SFRFNNWBX/Qdk7CvwcmUm0wWGQmclsgUS76xkOAol6VlP50TYzMys15Y4QTYzazq9qmItaSJwArAfsBGYHRHPljIwMzNrHJ6DXDJfBv5c7SCaxeLWrQzs34+DRmTUGjEzs4ZUVIIs6XDgB8BrMo7dA3w0Ip4uUWzVldUDXKBXOHN/oeUgMnuFs3t6M3uLC7XN6i0u1IHsnmUzqzQnyCUREV/u6rikMcDJJO8Af4uIVRUJrEEtat3K5P33paWf3zfNzJpFjxNkSYcBfwX2BxYCM4FVwIHAKcBrgZmSXhERC8oQq5mZ1SNXsS4ZSceTLKn464iYk3fsQuB7wKB01x5Jn4mI71Y4zIaxuHUrh3p4tZlZUymmB/nrJMnxJ4CrIqK944CkfsDHgCuArwHvKGWQ1aABnecgFy681Xl/5lxjsnuFC84rzuptLnAXO7MX2j3FZmaN5nzS9+HcnZJeBFwNtJDcwN4MnAp8W9LDEfHXSgda79rag6Vrt/Haow+odihmZlZBxSTIrwVuj4jv5R9Ik+XvSjoDOK1UwVVV1rDpQpWis/ZnrWEM2UOkCw6x7nkV68wVVJwgm1mtcA9yqZwMzIqI1Xn7LyZJjr8aEZcASDqFZL7yR0hGgFkRVmzYzq62dvcgm5k1mWIqUw8kew3kXLOBjK5XMzNratGHzXJNBB7J2H8asAO4vGNHRMwE7gFeXpnQGsvzSzwNrXIkZmZWScX0IP8DOKybNocBc3sfTu048s7NnfY9fc7wzLaRtbZxoSHWGT3LxQybLtg2a39TLcxlZrXMc5BLZjSwPneHpFHAZODPEbE9r/3jJL3OVqTFa7YAMHn0vlWOxMzMKqmYFOprwFslvT7roKQ3AGcDl5UiMDMzM+tkN0mSnOuE9PHvGe234n74XlncupWhg/ozZuig7hubmVnDKKYHeX/gD8DvJd0N3AesBsYCryZZ+ul3wGhJ7809MSJuLE24VVbE0k2F5iBn9Sxn9SpDgR7kUizzZGZm9WoBnWt9vI4kCX4go/04khUnrEiLWrdyyOghyPU8zMyaSjEJ8vUkb8AieXPOKsb1ZuBNOd8rPacxEuQC6yBnV6YukMj271tl6qKGWBeQWdDLzKyc3IdZKr8DvihpBvB94AiSIlw7gDsy2r8MWFy58BrHkrVbOWHiyGqHYWZmFVZMgvz+skVhZmaNy+sgl9K3SJZ6ujDdILkZ/bWI2JTbUNIxwFHA/1Q0wgawc08by9Zv560nTKh2KGZmVmE9TpAj4oZyBlIPCg2FzlyvuFBvc0ZPb1avMkB7ET3ImcOpi+hVNjMrKyfIJRERGyW9AvgqSXXqtcDPIuKqjObTSQps3lbBEBvC0rXbiIBDvMSTmVnTKaYH2aDgWsiZ1amLWd+4mKHTWddtj07Pp4DICFftnfdn7XNbt3Xb5mtbNk6QSyYiVgIf7EG7bwLfLH9Ejef5JZ6cIJuZNZuiE2RJY4BzgKOBIRHxwZz9hwCPZiwz0RhqNTku8HyFPvxm7Xdbt3VbtzWzxOI0QZ7sBNnMrOkUlSBLuhC4EtiH5wtwddzFHgv8DbgIuLaEMdaMgtWmsypTFyzo1TmRzRpKDWQP3S74ATijrUdYm1kNEJ6DXCmS3kyyqoSA+yLiV1UOqS4tad3K6KED2W/wgGqHYmZmFdbjPgRJpwMzgKdJ1jv+Ye7xiHgMeBx4SykDNDOzBhB92PpA0ihJd0qanz4WLEssqUXSbEm/79uzlo+kN0m6T9KrM479GPgN8HHgY8DPJTlB7oWOJZ7MzKz5FNOD/B/ASuDVEbFJ0gkZbeaSFA1pSEdd83Tm/sc/cVynfcUs3VRo2HR7xrLLxRTpcg+ymdWE6lax/ixwd0RcLumz6ff/UaDtJ4AngeGVCq4X3gy8GHgwd6ekNwLvA7YCVwCbSUZ0vUXSeRHx00oHWs8Wt25l2hFjqh2GmZlVQTEJ8lTg5vxlJPIsAw7sW0j1p5jK1FlDpLMS4YLXLWKIdUFOnM2seZwFTEu/vgG4l4wEWdIE4A3AZcC/Vii23jgJ+EtE7Mjb/wGS/vb3R8QvAST9N7AQeDfgBLmHNu/YzZrNOzlkjHuQzcyaUTFlWgaS3JnuygigrffhmJlZQ6rSEGtgbFr1GWAVSb2MLN8BPgO09/kZy+tAkulM+V4FbAD+OaQ6IlaRLPGUNeLLCljSug2AQz3E2sysKRXTg7wEOLGbNi8Fnup1NHWqqMrUGT3LxQzHLqZKrXuKzaxm9C3RHS1pVs73MyJiRsc3ku4ie/TSF/YKISKkzoO90+HJz0XEI5Km9SnS8hsJ7MrdIWkSMAr4XUTk/3yLSYZlWw8tat0CwCGjh1Y5EjMzq4ZiEuRbgM9IentE/CL/oKT3Ay8i7wOJmZlZH+cgt0bE1EIHI+K0gs8rrZY0LiJWShoHPJfR7GTgzZKmk6zSMFzS/0TEe/oUdXlsBibk7eu4eT27wDn5w7GtCx09yAfvv2+VIzEzs2ooJkH+BnAu8FNJbwP2A5B0MfBK4K3AfOB7pQ6y1mX29Bacg5zRttDaxllrlxYxX7mQrOJdXoLFzMqqen9jbiUpXnV5+nhLfoOI+BzwOYC0B/nfajQ5BngUeIOkoRGxJd13NskrPDOj/SEkBTathxa3bmH8iMHsM6DAG66ZmTW0HifIEbE+XVbiRuDtOYeuTB//ArwrIrqbp9xwihkKnV2ZukDbjN9OuapYOz82swZ1OclyRxcCzwDvAJB0EHBNREyvZnC9cBPw/4A/S7oBOIKkCNcq4E+5DSUJOAX4W6WDrGeLvcSTmVlTK6YHmYhYCkyT9CKS5Zz2BzYCD0TEI2WIz8zM6l1pim317qkj1gKvzdi/AuiUHEfEvSSVrmvVtSQjts4AppDcHt0NfCIi8otkvpZkbvZdFY2wjkUEi1q38pYp46sdipmZVUlRCXKHiJhLsuaxkT2cur3gsOmMtgV+C8Us8+R1kM2slnkaR2lERLukNwDnAa8A1gK/jog5Gc1HA98lGWZuPbBu6y4279jDZPcgm5k1raITZEkHA2NI+gPWpL3KJSfp7cClwNHASRExK91/OsmQuYEklTz/PSLuSY/dC4wDtqeXeV1EZBVkMTOzSnKCXDIR0U4y1PqmbtrdDNxckaAaxLPrk48Pk0a5QJeZWbPqUYIsaTTweZI71gfkHVtN8ib99YhYV8LYHiMZRvb/8va3Am+KiBWSjgPuAHLHQr27I5mulClf+XunfbMuy14RK4qYV5xVkKvgMk9Zl3APspnVCPcgWz1Ytj6pYD1h5OAqR2JmZtXSbYIs6XDgTmAiScq1h2RIl0jWXTwQ+FfgHEmnRcSiUgQWEU+mz5+/P3cZi8eBwZIGRcTOUjxvqRSuTJ1VxbrANbKqWBdMkLNKUxeKzszMzPItT3uQxztBNjNrWoVmtAIgqR9J7/Ak4M/AacDQiBgXEQcCw4DXAfcBk4H/KWu0nZ0D/D0vOb5B0hxJ/1f52XUOSRdJmiVp1vb1NZVbm5k1nujDZlYhy9ZvZ/g+/Rm+z4Bqh2JmZlXSXQ/y64CpwM+B8yJir48qaWJ6l6S7gZ+R9CKfHhF39uTJJd1F0gOd7wsR0WmtyrxzjwX+vzTGDu+OiOWShgG/As4nWZaqk4iYAcwAGHvMqJJ/BCvUK1zUMk9ZvdAFUv6CxbuyZF3DH0LNrFyc6FqdWL5hO+NHev6xmVkz6y5BPgfYCXwsPznOFREh6WLgzcDbSIZkdysiTutpoLkkTQB+A7w3IhbmXG95+rhZ0k+AkyiQIJdb4SrWPdtXcL/nIJtZnRH+c2T1Ydn6bRy8vytYm5k1s+76HV8M3B8Ra7q7UFotemZ6TtlIGgHcBnw2Iu7P2d8/LSaGpAHAG0kKfZmZmZl1KSJYvn4740d4/rGZWTPrrgd5IknS21OPk1S67jNJZwPfI1lS6jZJcyLiDOBi4DDgS5K+lDZ/HbAVuCNNjluAu4AflSKW3ijYK5wxxLpQb3PW7YuiepvNzGqFh1hbjduwbTdbd7W5grWZWZPrLkEeDmwo4nobSAp39VlE/IZkGHX+/v8E/rPAadlrK5mZWVV5mSerdcs3JBWsnSCbmTW37hLkgUBbEddrT89peoWXbsroLS5YpCvr/AJP6DnIZlbLnCBbjXt+DWQX6TIza2bdroOMP9b0yis+9nDm/r/MOKnTvqxEGFyky8waiN9JrMYt61gD2XOQzcyaWk8S5EslXVruQMzMzMyqZdn67QwZ2MKIfb0GsplZM+tJglxsP6T7CbpQ1DJPWUOsC/02ivkteR1kM6uk8Bxkq33JGsiDkTz8ysysmXWZIEe4NrKZmZWAE2SrccvWb/f8YzMz61EPspVQZuGtQnOQs25iew6ymdUh9yBbrVu+fhtTDx5Z7TDMzKzKnCBXWtYQ6wKJbOYQ62KqWJuZ1QonyFbDNu3YzaYde7zEk5mZFeqPNDMzM2sOyzsqWDtBNjNreu5BrrD2rB5kL/NkZg3OQ6ytlnUkyJ6DbGZmTpArrSmqA+QAAB0USURBVI9VrAslvS6nZmY1K/AQa6tpy9ZvA7wGspmZOUE2M7NKcIJsNWz5hu0M6t+P0UMHVjsUMzOrMifIFfaadz/Uad//397dR9ld1wcef38yCSEICWAgCUQKWlZFuqYa8alrUWmLlDXinrpSH0A9RNdqqe624rIesa67nG59aHVlRaXAKYIKpVClIlBR6REx2Ig8FEWMMkNmApJkAslMSOazf/x+KTcz9yZ3Hu79/Wbu+3XO79x7v7+nz72Eud/P/T7d9NWTmh88mdbm6XaxtvIqSepR/ZtdA1mSVDBBliR1VOAYZNXbwBbXQJYkFUyQ66BFq3CzNZNbtgo3KW+1fJQkdZ0Jsmqsf/MOnnfUkqrDkCTVgAlyDTRNhGFSXawno1nibOuOpE6K9I/MbBARlwCnA5sy88Sy7ALgHOCR8rD/npk3tHNuWX448GXgWGAD8IbM3Ny5dzE523fu4rEndroGsiQJcB1kSVKn5TQ3ddOlwKlNyj+ZmavKbUJyvJ9zzwNuyczjgVvK17Xx1BJPJsiSJFuQ66HVxFvNWpAn0cV6MpN0WQeVJGXmdyLi2Bk+dw1wcvn8MuBW4ANTuUcn9G8xQZYkPcUWZElSx0VOfVMtvDci7oqISyLisEmeuywzN5bPB4FlMxzbtPSXLchHH+okXZIkW5BroeXSTZMZg+yEXJLqzER3NrsI+CjFf8WPAh8H3j6VC2VmRrT+2SMi1gJrAY455pip3GLSBjbvYEFfcOQhC7tyP0lSvZkg10DHEuSZaHpxKmxJM8CW4NkrM4f2PI+IzwNfm+QlhiJiRWZujIgVwKZ93Oti4GKA1atXd+VfTf/m7Rx16CLmzfP7TpJkF2tJUjc4SdesVSa1e5wB3D3JS1wPnFU+Pwu4bibiminFGsiOP5YkFWxBroHXnHZH0/Kv33TShLKc16K22LEfvq2dSlKviIgrKSbUWhoR/cCHgZMjYhXFF8IG4J3lsUcBX8jM01qdm5lfBC4EvhIR7wB+Abyhm+9pf/o37+CVzz6i6jAkSTVhgixJ6iwn25o1MvPMJsVfbHHsw8Bp+zmXzPwV8OoZCXCGjTy5m0e2jbLyMCfokiQVTJBrrGlrcavxytNc5smGYkkd5d8Y1dDDW/bMYG0Xa0lSwQS5xrKvSeGk1kG2RiqpeoF/jlRPA66BLEkax0m6JElST/q3NZBNkCVJJVuQ66zpMk82w0iahdK/Xaqfgc076JsXLF98YNWhSJJqwgS5xpomwy26WDcfg9yhCqlrI0uaJLtYq44Gtuxg+eIDmd9nhzpJUsEEWZLUWa5nrJrq37zd8ceSpL2YINdYNvtBu+vrIDdjTVfS5MRY1RFIEw1s3sFLnvX0qsOQJNWIfYokSVLPeXL3GIPDI66BLEnaS21bkCPiD4ALgOcCJ2XmurL8WOA+4P7y0Nsz813lvhcClwKLgBuAczNn78wwa162bkLZdXe8sPnBzQb4uQ6ypLrwb4xqZnDrCGMJK10DWZLUoLYJMnA38Hrgc032/SwzVzUpvwg4B/g+RYJ8KvCPHYuwAtnXqou1tU9J9eWfKNXNQ5u3A66BLEnaW227WGfmfZl5//6PLETECmBxZt5ethpfDryuYwFKktqTFMs8TXWTOmDANZAlSU3UuQV5X46LiPXAVuB/ZOZ3gaOB/oZj+suyuaXFJF3hykuSaswWZNVN/+YdRMCKJSbIkqSnVJogR8TNwPImu87PzOtanLYROCYzf1WOOf77iHjeFO69FlgLcMhyJ+iQJKmXDGzZwbJDDuSA+bXtTCdJqkClCXJmnjKFc0aB0fL5nRHxM+DfAQPAyoZDV5Zlra5zMXAxwLITDp81bRvRagxy04M7F4ckTcqs+SurXuEayJKkZmZdF+uIOAJ4LDN3R8QzgeOBBzPzsYgYjoiXUEzS9Vbg01XG2hGt+imaDEuqqcAu1qqfgS07eMExh1UdhiSpZmrbrygizoiIfuClwNcj4sZy1yuAu8oxyFcD78rMx8p97wa+ADwA/Iw5NoO1JM1K05mgy0m61AG7x5KNW0ZsQZYkTVDbFuTMvBa4tkn5NcA1Lc5ZB5zY4dAqFS0m6WrWghyTaLLJtAlaktQbhoZH2DWWHH2oc5BIkvZW2wRZzb32eXc1Lf/a/b8xreuGAwQldZBdrFUn/eUST7YgS5LGM0GWJHWeCbJqZGDLdsA1kCVJE5kgzxExb2xiWYte0w7pk9RtVbUgR8ThwJeBY4ENwBsyc3OT4w6lmMPiRIp0/u2Z+b3uRapuGtw6CsDyxQdWHIkkqW5qO0mXJGmOSGAsp75Nz3nALZl5PHBL+bqZvwK+kZnPAZ4P3DfdG6u+hoZHOGThfJ620HYCSdLe/GaYI+ZNYn1kp+OS1EPWACeXzy8DbgU+0HhARCyhWCHhbIDM3Ans7FaA6r6h4RGWLbH1WJI0kQnyHDHPGXAk1Vl1f6KWZebG8vkgsKzJMccBjwB/ExHPB+4Ezs3MJ7oUo7pscHiEZYsXVh2GJKmG7GItSeq4yKlvwNKIWNewrd3r2hE3R8TdTbY1jcdlZtI8VZ8PvAC4KDN/E3iC1l2xNQdsGh5lmeOPJUlN2II8R/T1TZykazKcuEtSR03vj8yjmbm69aXzlFb7ImIoIlZk5saIWAFsanJYP9Cfmd8vX1+NCfKcNTaWRRdrE2RJUhO2IEuSOm6aLcjTcT1wVvn8LOC68Qdk5iDwUEQ8uyx6NXDvtO+sWnps+052jaUzWEuSmrIFeY7oa7LMkySJC4GvRMQ7gF8AbwCIiKOAL2TmaeVx7wWuiIgDgAeBt1URrDpvcOsIgC3IkqSmTJDniFeu/OmEsu9ufFbb54+lc1tL6pBWI3+7cevMX1G0CI8vfxg4reH1eqBlN27NHZu27UmQnaRLkjSRCbIkqaMCCCc6UE0Mbh0FYLnLPEmSmjBBnsP6ov1u1302IEvqJEeBqCaGhkeIgKUH24IsSZrISbokSVLPGBoeYenBC1nQZxVIkjSRLchz2PwWSz+l440ldZldrFUXxRJPth5LkpozQZ7DJjOztUmzpI6pcJIuabzB4VGOPtTxx5Kk5uxfJEnqsIScxibNoE3DIxzpEk+SpBZsQZ7D5k9iki5sQJbUQWGeqxoY3bWbXz2xk+UmyJKkFkyQ57AFfbvbPtZ1kCVJc90j24olnhyDLElqxQRZktR5dpVWDQwNjwCwzBZkSVILJshz2AlLBpuW/2TbkV2ORFJPS5jMiA+pUwa37mlBNkGWJDVngixJ6jxbkFUDe1qQHYMsSWrFBLkHTWryLkmS5oih4REOmD+PQw9aUHUokqSaMkHuQQc0mbzLSbokdZQNyKqBoeERli1eSITfeZKk5kyQJUkdF3axVg0MDo+w7BC7V0uSWjNB7kFNu1j7Y7qkTjJBVg1sGh7luUctrjoMSVKNzas6AEnSHJfA2DQ2aQZkpi3IkqT9sgW5Bx0wb9eEsjGbkCVJc9jjo7vYvnM3y5csrDoUSVKNmSD3oPnzJk7SJUmdEqRjkFW5PUs8uQayJGlfTJAlSZ1ngqyKDQ2PAibIkqR9M0HuQYvnj0wo27H7gAoikdQzTJBVscGtxXffchNkSdI+mCALsNu1pA7aM0mXVKGhbXaxliTtX21nsY6IP4iIeyJiLCJWN5S/KSLWN2xjEbGq3HdrRNzfsO/I6t6BJEmqi6GtIyw+cD6LDuirOhRJUo3VuQX5buD1wOcaCzPzCuAKgIj4DeDvM3N9wyFvysx1XYtyjlgQtiBL6hwn6VLVhoZHbT2WJO1XbRPkzLwPIGKfyw+dCVzVlYAkSVNngqyKDQ6PsHyJCbIkad9q28W6Tf8ZuHJc2WVl9+oPxX6ya0lSN2SRIE91k2bA0PAIRx5igixJ2rdKW5Aj4mZgeZNd52fmdfs598XA9sy8u6H4TZk5EBGHANcAbwEub3H+WmAtwCHLD5pK+JIkaRYYG0s2bRtl+ZKFVYciSaq5SluQM/OUzDyxybbP5Lj0Rsa1HmfmQPm4DfgScNI+7n1xZq7OzNWLDvMLU5I6JrEFeZaIiEsiYlNE3N1QdkFEDDRMgHlai3NPLSfKfCAizpvs+Z306BOj7B5LxyBLkvZrVnaxjoh5wBtoGH8cEfMjYmn5fAFwOsVEX5Kkqo1NY1M3XQqc2qT8k5m5qtxuGL8zIvqA/wu8BjgBODMiTmj3/E7bNDwKuMSTJGn/apsgR8QZEdEPvBT4ekTc2LD7FcBDmflgQ9lC4MaIuAtYDwwAn+9awJKkliJzypu6JzO/Azw2hVNPAh7IzAczcyfFD9hrZjS4aRjc6hrIkqT21HkW62uBa1vsuxV4ybiyJ4AXdj4ySdKkmejOdu+NiLcC64D/mpmbx+0/Gnio4XU/8OJJnN9RQ9uKBHm5CbIkaT9q24IsSZJq4SLgmcAqYCPw8U6dHxFrI2JdRKx75JFHphjuRENbR5gXsPTgA2bsmpKkuckEWZLUWQmM5dQ3VSozhzJzd2aOUQxdajYB5gDwjIbXK8uyds/fc69/m0DziCOOmLH3MDQ8ytKDFzK/z2qPJGnf/KaQJHWY6yDPZhGxouHlGTSfAPMHwPERcVxEHECx0sT1kzi/owaHRxx/LElqS23HIEuS5hAT3VkhIq4ETgaWlhNlfhg4OSJWUfQF2AC8szz2KOALmXlaZu6KiPcANwJ9wCWZeU952b9odn43DQ2PsPKwg7p9W0nSLGSCLEnqPBPkWSEzz2xS/MUWxz4MnNbw+gZgwhJOmfmWGQtwioaGR3jhrx1WdRiSpFnALtaSJGnOGt21m83bn3QGa0lSW2xBliR11p5JuqQKbBoeBVwDWZLUHhNkSVKHJeRY1UGoRw0NF2sgL1tigixJ2j8TZElS5zkGWRUZLBNku1hLktrhGGRJkjRnDW4tW5AXL6w4EknSbGALsiSpsxyDrApt2jbKwvnzWLJoQdWhSJJmARNkSVLn2cVaFRncOsKyxQcSEVWHIkmaBUyQJUmdZ4KsigwNjzj+WJLUNscgS5I6LIsEeaqbNA1DwyMc6fhjSVKbTJAlSdKclJkMDY/agixJaptdrCVJnZXAmOsgq/uGR3ax48ndLDNBliS1yQRZktR5dpVWBTaVayAvW2KCLElqjwmyJKnzTJBVgcE9CfIhjkGWJLXHBFmS1GHpOsiqxNDwKADLbUGWJLXJSbokSdKcNLSnBdkxyJKkNtmCLEnqrIRMJ+lS9w0Nj7Bk0QIOXNBXdSiSpFnCBFmS1Hl2sVYFVh62iFc958iqw5AkzSImyJKkznOSLlVg7SueVXUIkqRZxjHIkqQ5KyIOj4ibIuKn5eNhLY57X0TcExF3R8SVEeGgVUmSepAJsiSpszJhbGzq2/ScB9ySmccDt5Sv9xIRRwN/DKzOzBOBPuCN072xJEmafUyQJUmdlzn1bXrWAJeVzy8DXtfiuPnAooiYDxwEPDzdG0uSpNnHMciSpI7L6bcET9WyzNxYPh8Elo0/IDMHIuIvgV8CO4BvZuY3uxijJEmqCRNkSVKHTbsleGlErGt4fXFmXrznRUTcDCxvct75e0WRmRExIZByXPIa4DhgC/DViHhzZv7tdIKWJEmzjwmyJKnuHs3M1a12ZuYprfZFxFBErMjMjRGxAtjU5LBTgJ9n5iPlOX8HvAwwQZYkqcc4BlmS1FlJsQ7yVLfpuR44q3x+FnBdk2N+CbwkIg6KiABeDdw33RtLkqTZxwRZktR5OTb1bXouBH4nIn5K0VJ8IUBEHBURNwBk5veBq4EfAj+m+G68uPnlJEnSXGYXa0lSRyWQ028Jntq9M39F0SI8vvxh4LSG1x8GPtzF0CRJUg3VtgU5Iv5PRPxrRNwVEddGxKEN+z4YEQ9ExP0R8XsN5S+MiB+X+/667ConSapSZpUtyJIkSW2rbYIM3AScmJn/HvgJ8EGAiDgBeCPwPOBU4LMR0VeecxFwDnB8uZ3a7aAlSZIkSbNTbRPkzPxmZu4qX94OrCyfrwGuyszRzPw58ABwUjk76eLMvD0zE7gceF3XA5ckTZBjOeVNkiSpW2bLGOS3A18unx9NkTDv0V+WPVk+H18uSaqaXaUlSdIsUGmCHBE3A8ub7Do/M68rjzkf2AVcMcP3XgusLV8+/tcvuOr+mbx+zS0FHq06iBrwcyj4ORT8HArPnukLbmPzjTfn1UuncQn/u/SgO++889GI+MUkTvH/4fb5WbXPz6p9flbt87Nq30x+Vr/WzkGVJsiZecq+9kfE2cDpwKvLbtMAA8AzGg5bWZYN8FQ37MbyVve+mB5dxiMi1mXm6qrjqJqfQ8HPoeDnUIiIdTN9zcx0PghNWmYeMZnj/X+4fX5W7fOzap+fVfv8rNpXxWdV2zHIEXEq8GfAazNze8Ou64E3RsTCiDiOYjKuOzJzIzAcES8pZ69+K3Bd1wOXJEmSJM1KdR6D/BlgIXBTuVrT7Zn5rsy8JyK+AtxL0fX6jzJzd3nOu4FLgUXAP5abJEmSJEn7VdsEOTN/fR/7PgZ8rEn5OuDETsY1R/Rk1/Im/BwKfg4FP4eCn4NmK//tts/Pqn1+Vu3zs2qfn1X7uv5ZxVNDeyVJkiRJ6l21HYMsSZIkSVI3mSD3kIi4JCI2RcTdVcdSpYh4RkR8KyLujYh7IuLcqmOqQkQcGBF3RMSPys/hI1XHVKWI6IuIf4mIr1UdS1UiYkNE/Dgi1ndiNmupEyLi1Ii4PyIeiIjzqo6nzqwHtM+6QvusT0yO9Y32VVUvsYt1D4mIVwCPA5dnZs+O1Y6IFcCKzPxhRBwC3Am8LjPvrTi0ripne39aZj4eEQuA24BzM/P2ikOrRES8H1gNLM7M06uOpwoRsQFYnZmuzahZISL6gJ8AvwP0Az8Azuy1v+ftsh7QPusK7bM+MTnWN9pXVb3EFuQekpnfAR6rOo6qZebGzPxh+XwbcB9wdLVRdV8WHi9fLii3nvzFLCJWAr8PfKHqWCRNyknAA5n5YGbuBK4C1lQcU21ZD2ifdYX2WZ9on/WN2cEEWT0tIo4FfhP4frWRVKPs5rMe2ATclJk9+TkAn6JYd32s6kAqlsDNEXFnRKytOhipDUcDDzW87sckRjOs1+sK7bA+0TbrG5NTSb3EBFk9KyIOBq4B/iQzh6uOpwqZuTszVwErgZMioue63EXE6cCmzLyz6lhq4LfKfw+vAf6o7I4pST3LukJ7rE/sn/WNKamkXmKCrJ5UjpG5BrgiM/+u6niqlplbgG8Bp1YdSwVeDry2HOdyFfCqiPjbakOqRmYOlI+bgGspuq9KdTYAPKPh9cqyTJo26wqT1+P1if2xvjFJVdVLTJDVc8rJJL4I3JeZn6g6nqpExBERcWj5fBHFJDf/Wm1U3ZeZH8zMlZl5LPBG4J8y880Vh9V1EfG0ciIaIuJpwO8CznSruvsBcHxEHBcRB1D8P3x9xTFpDrCu0D7rE+2xvjE5VdZLTJB7SERcCXwPeHZE9EfEO6qOqSIvB95C8cvd+nI7reqgKrAC+FZE3EVRybwpM11yoHctA26LiB8BdwBfz8xvVByTtE+ZuQt4D3AjxSRKX8nMe6qNqr6sB0yKdYX2WZ9QJ1RWL3GZJ0mSJEmSsAVZkiRJkiTABFmSJEmSJMAEWZIkSZIkwARZkiRJkiTABFmSJEmSJMAEWXNQRNwaEbNmevaIyHHb8kmcu6FccL6WIuJr497b2VXHJElSM9Yf6sP6g6o0v+oApFam8CX1tsy8tBOxdMEvgEvL549XGMdM+xKwDlgFrKk4FklSD7D+MCdYf1BlTJBVZx9pUvYnwBLgr4At4/atLx/fChzUwbg6YUNmXlB1EDMtM78EUP7y6xecJKkbrD/MctYfVCUTZNVWsz/45R/KJcCnMnNDi/N+2dHAJElSbVl/kDQdjkHWnNNsDFFEnFyOYbkgIlZHxDciYmtEbI6IayLiGeVxz4yIqyLikYjYERHfiojnt7jPQRHxwYhYHxFPRMTjEfG9iDhzht9PRMR7IuKeiBiJiIGI+ExELGlx/JKI+NOI+KeI6I+IneX7uT4iXjru2MMiYntE/CwiosX1/qH87FY3lL02Im6JiI0RMRoRD0fEtyPi3TP53iVJ6hbrD9YfJDBBVu95EfDd8vnngTuA1wM3R8RzytcrgcuBrwO/DdwUEQc3XiQiDgVuA/4XsBu4BLgMOAL4UkT8zxmM+VPAp4HDgIuBq4BTgZuBA5oc/1zgY8BY+R4+AdwEvAr4TkScuufAzNxcXu+ZwCnjL1R+8b8GuDMz15Vla4HrgBOAfwA+DtwALALeNu13K0lS/Vh/sP6gXpGZbm6zZgM2AAkcu49jbi3+ae9VdnJ5XgJvGrfvi2X5Y8D54/Z9qNx37rjyS8vyPxtXfiDwDYovl1VtvqcEbm2x72Xl/geAw8fd53vlvg3jzlkCLG1yrZXAw8B948pXl9e5usk5F5T7zmkouxMYBY5scvyE+5blZ5fXObvqf0Nubm5ubr23WX/Y6z7WH9zc9rHZgqxec1tmXjGu7LLycStw4bh9l5ePq/YURMTTgTcD6zLzLxoPzswR4ANAAH84A/Hu+UX1Y5n52Lj7fLDZCZm5NTMfbVLeD1wNPCcijmkoX0cxU+SaaFgiIiL6gHcA24Arx11uF/Bkk3tMuK8kSXOA9QfrD+oRTtKlXrOuSdnD5eP6zNw9bt9A+biyoexFQB+QEXFBk+stKB+fO9UgG7ygfPx2k323UXTPmiAiXg6cC7wUOJKJXamOBhonI/ksRTevt1N0+wI4jeJ9X5SZjUtHXEHRLereiLiqjO2fM/ORNt+TJEmzjfWHgvUHzXkmyOo1W5uU7Wq1LzN3lXNPLGgofnr5+KJya+Xgfexr156JNIZaxDbhF9eIOIPil94RirFDPwOeoOi2dTLFuKiF4067iuJL65yIuDAzx4C15b7PjbvvJ8r7vhv4Y4qlMzIivg38afmLsiRJc4n1B+sP6hEmyNLk7fki/GRmvr9L91oGPNi4IyLmA0uB/nHnfBTYCazOzPvGnfM5ii+4vWTmjoi4FHgf8LsRcQ/F5Brfz8wfNTn+cuDycrKRlwFnUPx6fGNEPMdfgyVJmsD6g/UHzQKOQZYm7w6KX1P/Qxfu9cPyccKXEvBbFF21xvt14N4mX27zynNauYhiIox3Uowd6mPcr7/jZeaWzLwhM8+hmHjkcOAV+zpHkqQeZf2hZP1BdWaCLE1SZm6iGEezOiI+VE5GsZeIeFZEHDcDt7u0fDw/Ig5vuP6BwP9ucc4G4PiIOKrh+KCYUfKEVjfKzJ8CtwCnA+8CtlB0ndpLRLyyxZqHR5aP21vdQ5KkXmX9wfqDZge7WEtT8x7geODPgbdExG0U43yOophc40XAmcDPp3OTzPzniPg08F7g7oi4mmL2xzXAZmBjk9M+Cfw/4F8i4pry+Jfz1LqD/3Eft/wsxXqGy4BPZ+aOJsdcCzweEbdTfJkGxa/hL6JYwuHmSb5NSZJ6hfUH6w+qOVuQpSnIzGGKbkvvBR4F/hPwfuCVFMsavI9igouZcG55n60U3ZfOBG6k+CLa2SS2z1Es77AROAt4E/AQ8GKe6nLVyvUU7wdad486D/gBxQyZ7y7vtYBieYpXZuaE5RskSZL1B6w/aBaIzKw6BqmnRUQC387Mk2sQyzOBByiWXZixMVIRcTbwN8DbMvPSmbquJEm9yvqD1Bl2sZbq4bfLLzqAFZk5WFEc/42iy9NnZuJiEfE14Pdn4lqSJGkC6w/SDDNBlqr3kXGvH+/mzSPiGOAPKcZEvQ34EfDVGbr8l4DGdQ3Xz9B1JUnqddYfpA6wi7XU4yLiZOBbFLNH3gb8l8x8cJ8nSZKknmb9QXOVCbIkSZIkSTiLtSRJkiRJgAmyJEmSJEmACbIkSZIkSYAJsiRJkiRJgAmyJEmSJEmACbIkSZIkSQD8fwtWomTzQuUyAAAAAElFTkSuQmCC\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x7f6673c247f0>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"time = Dataset('prog_z.nc').variables['Time'][:]\n",
"temp = Dataset('prog_z.nc').variables['temp'][:,:,1,1]\n",
"temp_drift = temp-temp[0,:]\n",
"\n",
"plt.figure(figsize=(16,6))\n",
"plt.subplot(121)\n",
"plt.pcolormesh(time, -z, temp_drift.T)\n",
"plt.colorbar(); plt.title(r'$\\theta(z,t) - \\theta(z,t=0)$ [C]', fontsize=20)\n",
"plt.xlabel('Time [days]', fontsize=20); plt.ylabel('Depth [m]', fontsize=20); plt.ylim([-200,0])\n",
"plt.subplot(122)\n",
"plt.plot(time, temp[:,0])\n",
"plt.title('SST evolution [C]', fontsize=20)\n",
"plt.xlabel('Time [days]', fontsize=20); plt.ylabel('SST [C]', fontsize=20)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Vertical viscosities ($\\kappa_v$) and diffusivities ($\\kappa_d$)"
]
},
{
"cell_type": "code",
"execution_count": 9,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"(-200, 0)"
]
},
"execution_count": 9,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAA9EAAANJCAYAAAALBb9dAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAIABJREFUeJzs3Xu8XGV97/HPjyDiPRqwIJCCivWCNQpG66mKioJXvBf1SNVaylHSnvacai09LbW1pbXWC2pptDZyTipVKYIVi2KLVNvIpRLlIhgxcpGqAVEbFEjyO3+stc1kMrP3zFqzZtbMfN55zSt7nnV7Zs3s+e5nXZ4nMhNJkiRJkrS0PSZdAUmSJEmSpoWNaEmSJEmSBmQjWpIkSZKkAdmIliRJkiRpQDaiJUmSJEkakI1oSZIkSZIGZCNakiRJkqQB2YiWJEmSJGlANqIlqYUi4i0RcUlE/DAivhcRn4yIw1pQr6dExLkRcVNEZES8ZtJ1kiRp3Nqa02BWj4ONaElqpyOB9wNPAp4ObAMuiIgHTLJSwL2BK4DfAH484bpIkjQpR9LOnAazunE2oksR8eGI+G5E3GvC9Ti8PGL0+gnXY3NEbG56mWENs42IOLjcl+uarNM0cx8NLiKOLPfVwuNrTW4vM4/OzL/NzCsy86vAq4F9gf/W5HYHqNd5mfm7mflxYEeveSJin659lWOupqZcGzK5LXlc1mWofDWPp4/7aDTGmdVtzWkwq8fBRjQQEY+n+OCfmplbx7TN3yw/sK/sLM/My4BPAH8UEfceR10kDeXzwB8C7x3zdu9D8Z39/TFvt4rbKfbRHwLfmnBdNGXGncnmsTSTJpHV05TTYFbXsuekK9ASbwN+CPzVGLd5ePn/ZT2m/SnwJeDXgT8ZW400j24CHgH8YNIVmSIXZuYpE9juu4HLgX+fwLaHkpm3A6dAcVYA+NlJ1kdTZ9yZbB6rDczj0ZpEVk9NToNZXdfcn4mOiIcBRwEfzcxx3jNwOPBfwLXdEzLzYuBrwK9FxNy/R2pOZt6VmV/LzJsnXRf1FxF/Cfwi8JLM3D7idf9x9+VcPR5HjnKbUj8TymTzWBNnHk+3JnO6XL9Z3TJTFwgR8Yjyg3J2V/kTyh7ovhURjx1ila8DAvj7Ptv7whIf2IuGrP+p5T0HD6e46X9Hx7pe3THrmcBK4JlDrPs1EXFWRFwXET+OorfAL0bEf+8zf0TESRFxZUT8pNx/742I+y2yjSrLPCEiPh4R/xkRd0bEDRHx1xHxoFFtYykR8fCI+ERE3BoRW8v39Vld0zMi/mWRdXw1Iu6KiP2H2O7qiPj78jXcERE3R8RnIuLl5fQn9vo8d63j6nLZB0THPVMR8ZByv94SET8q13tYucy+EbG23N5Poug98mk91r3bPVjlfsqI+PUe8/9ROe1vBt0Hw+6TrnlfHhEXRcQPys/0V6PoDfPui72W8uczI2JL+fovjYjndcw71H4f8HWN9L3pWO87gVcAT8/M6wapS7ncCyLic+V27oiIb0fE5yPiDV2zvovi7Mdij4sH3a7mR4w+j2GMmRwtyuNymYghsm/Y+TuWGziTq25jie0vmscd84w0k2PEeVw+H9n3fvS5JzqmMJNjiDwu5x95Ji+ynpFndVTM6XJZs3paZeZUPYBXAQn8n46yXwZ+AlwI7Dvk+i6l6E3vXn2m/2+KSx26H98q63HKkNt7ObCuXPaLXes8qGO+o8p5/mKIdf+4fD3rKC5BWwvcWK7nj3rM/+5y2reB9wDvADYBl5Rlm+suQ/EH0TZgK/AR4M+Bs4Ht5fwrR1GvPvvj4HI9n6e4P+Wicr+sK/fVduCXOub/53L+h/VY15PKaR8f4v341fK13wF8jOJSwA9SXOpzYcd8XyvnWdFjHas7t9vxmi4EtgD/Wu6fsyg6jtgCHAp8A/gyxZfuGcCdFL8jK/vso3UdZQ+g+Hz/BHhsR/kzyn12JXDPir+/A+2Tct4/Kev2PYrLOt9O0dPkwuvfq89r+Rfgu8AG4J3Ah8vXsh14WpX9XpYdySK/86N+bzp+F/4TeMSQ+/mEsi43U3wPLOzni4FLqrx3fbbzX8BrlpjnQiBHtU0f7Xkw4jwulx9bJtOiPC6XGTZfq2T4UJlcZRt9XtvBDJHH5TIjy2QayOOu13UhNb/36ZHHZfnUZTJD5nGVfb/IazqSMWY1FXO6XNasnuLHxCtQ4YPwjvID91xgWflLmcD7gLsNua57lV8eXx1yub8ot/khYI8Kr2Hhl+aERea5XznPxUOs9yE9yvYCPgfcBRzQUb4QQJuAB3SU701xL0eye2APtQzwsPJLZ1PntstpC1/+Z9fZxhL7Y+GLMoG3d007otwn3wfuW5a9lD5/KLHzD61nDrjtR5brvxV4VI/pB3b8/JZy3Sf1mO995bTn93hNJ3fN+3/K8luB0zs/mxSd9CTwzj77aF2P9+Euissb7w38DMWX/O29Xk8D++QXynpdD+zXUb4n8Mly2u8u8n7/Qde0o8vy86rs97LsSAYL5lG9N++juC/06cB+HY97D7CvL6P4Y+SBPabtU+X961j+3sCq8nE78Pvlz7sdBCjnvxCDeSYfjDCPy/WNPZNpQR6X04bN1yoZPlQmV9nGIvuj8/txyTwuy0eSyTSUxz1eV63vffrkccd7MTWZzJB5XGXfL/K6jmRMWU2NnC6XN6un+DHxClT4UPxL+SF+FPDZ8sP3+orreli5rs8MOH9QjAeXFL39RcXtnl6u44gl5vsx8J8j2GcvLrd3fEfZB8qy1/aYf+ELaHNX+VDLsPMPquf2qdfZFH8w3adOvRZ53QtflLd1bqNj+rpy+i+Xz/ekOLK+Bbh7x3zLyy+gTYO+58Bp5bp/c4B5D6T44+WSrvK9gFuA7wB7dr2mbwLLuuZfWU7b2v16Kf7AvQv4lz77aF2Pev1OOW098Jny50q/axX2ycLnYLc/bMvf2+3AdX1ey+bufVNO/xawpcp+7/r8nbLE521U7032efTcfteyl5Xbun/V92uRdR/Zp167fYbK+S/EYJ7JByPM43J9Y89kWpDHZfmw+Volw4fK5CrbWOR1L3w/DpTHZdlIMpmG8rjrddX+3meRPC6nT00mM2QeV9n3i9R14bN5yhKfxVG8Z5VzulzerJ7ixzT2zr2K4kjqOcBDgC9l5gcrrmtF+f+SXdFHxDKKo9zHA3+emW+uuE2Ax1H8In51iflupTjaOJCIWAm8meKI8krgHl2zHNBVBygurer2BYovsm7DLvML5f9PjWLIkm4PpPhSehg7e0WtUq+l/Edm/qhH+YUUlx4+FvhwZm6LiA9QHK17CfB35XyvptiXa7P8phnAE8v/P73UjJl5Y0R8DnhmRDwyM68qJz2f4jKud2bmtq7FLs/dO674dvn/td2vNzO3R8R3KEJqUH8GPA1YGPblIzV+12CIfcLOz8E/d0/IzGsj4kbgkIi4X2Z292Taa98A3MDOz2TV/T6Ikbw3mRkVtr1gPcVZwqsi4kyK36cvZub3aqxzoV4XUjRepFp5HBEnAv8rMw8tiyaRyW3I44V6wODZVyUrh83kieUxwAgzuek8BjO5XyYPlMfleprK5H5qv2c1cxrM6qk2VY3oiHgwxRFIKO6xuB14QkQ8PzM/WWGVCz1/7r3Edu9G8eX9UorLUt5aYVsL69oTeDRwVWbescTs9+io41LrfTDFPRT3p7i/4zMUwyRspzjq9stAZ8cPC52CfKd7XWVwbemxmWGXWfiD6LeXqH7n+JtV6rWU3dZV+s+ubUJxT8rJwK+xM7BPoPhD8W+H2ObC5/SmAedfR9FpzS9T/OFF+TOUf1B02W0IjHL/9JxW2gbcbcD6kJkZEf8ALHT48q5Bl+1jmH2y8J7066X0Zoo/TJez++u9rc8y29i9M8V1DLffB9H4e7OUzPzL8nflDRRD8/xPICPi88BvZ+alo9qW5tOI8ngVxb2XC8aayS3KYxg++6pk5bCZPOk8htFkctN5DGZyv0weJo+hmUzux6xWLdPWO/fCkbCzKS5TOLl8/raoNvTEd8v/V/SbISL2Lrf3UuB/12lAlx5J8QfCfyw2U/l6lnfUcSm/RfE6fiUzj8zMX8/M/5PFGHnn95h/4QtityPr5R8W+4xgmYX575eZscjj8z2WGaZeS+l39mC/rm2SmTcB5wJPiaJ30CcBh1HcJzbMkcGF4Og+29DP2RT31fz3iFgWEQ8Eng1szMyNQ2x3ZCLiUIp7Db9P0dnGB8vfh6qG2ScL78l+fabv3zVfVa3b76OSmWdk5hMpvheeC/wN8BTg/IjYd6KV0ywYRR6voui4Z8G4M7kteQzV87VKhg+ayRPNYxhZJk99HoOZPO2Z3I9ZPb2mtRG9PjPvLI92f4niSPIuw0ZExHci4qSuskdEMZzDqrLoZopeBn+u18Yi4l7Ap4DnAG/IzHeM4DUsbPvLi85V1CnY9Sj9Yh5a/n9Wj2lP7VH2H4tM+0WKS7rqLrOh/P/JPebvp0q9lvK4iLhPj/Ijy/+734v3l///GsURb4C/HnKbC6/92YPMnMV4qB8FHkTRE+wrKa4UGfWR14FEMVzF31N09PNLFL2oPpp6R76H2ScL78mRPer2UIrLqb6Zmf2Ocg+kbfu9CZl5W2ael5m/SnGU/wEUAS3VMXAeA0TEqoj4l9g5LM7qct7OjBt3Jrclj2H47KuSlcNmchvyGOpn8lTnMZjJzFAm92NWT6FswY3Zgz6Af6K4Kf7BHWULQ098k127178AeH/X8p8A/m9X2cfL5R/aVX4/iiEvttHVAUifuq0r1/OaJeb7Tbo6z+gz32vL+XbrpbDP/Audozy/q/zo8jXs0tEB8N/KsmF69hxqGYqxN++k6E2y1xAVewFPrrONJfbJwezsSKFfb6C30dEbaDktgGso7oG7Hfhahc9qZ6+Xj+wx/cAeZQuvfT3FHy930dVjI0t3PJJ0DUvRMW1zj/e05/rY2eHIqeXzZRT3wCXwsmH3x7D7hJ29wn6TjmFyynp8opzW3avmUvvmQnp0mjHIfi/nO5LBOisZyXtT50Fx39xuHe6wsxfVZ49qWwPUped+9zHdD4bL40MpzlCdStHAfG75mU/gQV3rHVsm05I8LqcNm69VMnyoTK6yjUX2ycL341B5XE6vlck0lMddr2tdn20P/L2/2LqYskweYL9cSJ9cGHTfL/K6jsSsrlKXvu+Jj96PqbonmqLDiR9Q/BIDkJkXRMSFFL80/4NivDYoOgl57MJ85SVAx7D7Ee6zKDqsOJoiKBb8HcWXxsXAgyPilB71+dPceR/Vwln9pTo9WOhA623lgO5bgSsz82Nd8z2L4v6pc5ZY34L3UwT9xyLi4xSdIxxG8Zo/SnHk8qcy84sRcRqwBriiXOYu4FiKS4V2u+9l2GUy82sR8TqKzl+ujIh/ogjvu1HcO/NkirMOD69TrwFcBLw+Ip5A8UfY/uX+2AP4tcz8YVe9MyJOB/6yLFo77AYz86qIeAPFH1NfjohzgK9TXK7zeIrLlZ7WtcwXI2IT8DKKffTJzBz08sGRiYgXASdRnFX6vbJu2yPiFRRnYj4YEZdl5nXDrHeYfZKZ/xYRfw68iZ2fg60UR8wPo/jj4e21Xyzt2e8jdjbwXxGxgSL0g+L37fEU30EXTK5qmhHD5PF7gU9l5u+UzzdFxAuBF2Tmt9nVODO5FXkMlfK1SoYPlcltyOOyHrUyeZrzGMxkZiOT+zGrp9mkW/GDPiju2eh5dIidR8i+y86hGV4HfLdjni8A7+ix7F4UHV18qaNsD+BH7Dxq2uvxna71fJniC2fJbuopvgyvoRi0PYG3dU2/H0UHJp8Ych89iaLnxO+X9f8C8EL6HJWj+GU9CbiaYmiSb1OMeXc/+hxtq7jMoynOCnyrXOZW4AqKy7GePopt9NkfB5evex3wCIo/gL5PcST7i8DRiyx7f4o/mn4MrKjxuf0Fij8Kv0txBuDbFGdwXtpn/t/r+Iy9ZLHX1Gf5oY6gdq+P4g+pWynOCBzcYx3HlvNfTMeZpqb2CXBc+Tn+EcXvy5UU917uXWHfXEj/I9+L7vdyniOZnqPbJ1KE83Xl5/1Wiu+oN9FjeJkmH4vtdx/T+WCIPAYOKp8/vmu+9wPn91h+rJlMS/K4XGao7Bt2/o7lBs7kqtvosc2ffj8yZB6Xy9fOZEacx92vq8/0WmeimdJMHmC/XMgiuTDIvl9k2SMxq6vUZdH3xMfujyh33MyJYviGiyk6vvhFii/uh2TmrT3mfQvwJ8DjMnOpe6N6bWs5xRh278jMN9Wpd7m+NcB7KC6r+kLd9Wl4EXEkxRio/y8zXz3h6qgFOj4Tf5hFB0EaQHlm8qlZfygQTaGIeAHwD8A9M/POjvINwOezx9BUbcpk87gdzGQNyqyuxqwe3rR1LDaMKyl6L3w0RRif2qsBXXoncD1QtZfPJ1Nc4vSXS824lIi4B/AW4CwDe6IW/vB670RroTb6g4jIiPjapCvSVhGxT7mPkv4dKWk+JMV9kz8d0ikingw8gf4derUik83jVjGTNSyzeglmdT3Tdk/0wDLz9oi4jqIBfV+KI8n95v1JRLwaeFpE3Csztw65rU+yxLiWQziY4n6fdSNanwYUEY8GngccTnGfzz9m5pcmWyu1yGbgDzueVxkjdV7czq77SvPrUorLgP+ivKfy4ezsbblnb9ctyuSDMY8nxkxWRZsxqwdlVtcwk5dzR8QxFB2aHEAxHMDrMvNvJ1srtV1EvAb4W4r76M6nGELFL98+yksm/+eAs6/LzM0NVkdzLiI+RPEH93cz87Ae04MiF55D8YfDazJz0fGBNRplB0inUowxu4PivtZ7U9xHuWOSdVN7mcnDMZPVRrOczTPXiI6IZRS9TT4TuBG4BHhFZl410YpJMyYiDqajZ94lPC0zL2ysMpp7EfEU4L+AM/oE9XMoehl+DsWlxO/OzCeMt5bzy2yWmmUmq41mOZtn8XLu1cCmLLv6j4gzKXovNKilESqPYtsBhVohMy8q/4js51iKEE9gQ0Qsj4j9M7PKMD0antksNchMVhvNcjbPYsdiBwA3dDy/sSyTJM0vs2Gy3P+SpG5Tmw2zeCZ6IBFxAnACwL3uGYc//KF7TbhGkjR5l33lji2Zue8o13n00+6Vt9y6vdY6LvvKHVdSjEu6YG1mrq21UrWO2SxJuzOb22cWG9E3AQd1PD+wLNtF+QavBTjiMXvnxeevHE/tJKnFlu3/9W+Nep233Lqdut+xy/b/+k8y84gaqxgoG9QYs1mSKjKb22cWL+e+BDg0Ig6JiL2A44BzJ1wnSZpbCeyo+W8EzgWOj8ITgR9Mwz1XM8RslqQWMZvrmbkz0Zm5LSJOohgOYRnwocy8csLVkqQ5lmxveCSjiPgIcCSwT0TcCPwBcDeAzDwdOI+i989NFMNovLbRCmkXZrMktY3ZXMfMNaIBMvM8ijdFkjRhxdHuZodTzMxXLDE9gTc2WgktymyWpPYwm+uZxcu5JUmSJElqxEyeiZYktcuI7p2SJEkjYjZXZyNaktSoJNmezV4yJkmSBmc212MjWpLUuKbvu5IkScMxm6vznmhJkiRJkgbkmWhJUqMS2O7RbkmSWsNsrsdGtCSpcV4yJklSu5jN1dmIliQ1KsHOSyRJahGzuR7viZYkSZIkaUCeiZYkNc6RKCVJahezuTob0ZKkRiVp5yWSJLWI2VyPjWhJUrMStpvTkiS1h9lci/dES5IkSZI0IM9ES5IalXjflSRJbWI212MjWpLUsGA7MelKSJKknzKb6/BybkmSJEmSBuSZaElSoxLYYeclkiS1htlcj41oSVLjvGRMkqR2MZursxEtSWpUYlBLktQmZnM93hMtSZIkSdKAPBMtSWrcjvRotyRJbWI2V2cjWpLUKC8ZkySpXczmemxES5IalQTbvXtIkqTWMJvrcc9JkiRJkjQgz0RLkhrnfVeSJLWL2VydjWhJUqO870qSpHYxm+uxES1JaliwPb17SJKk9jCb63DPSZIkSZI0IM9ES5IalcAOj9lKktQaZnM9NqIlSY3zvitJktrFbK7ORrQkqVGZ3nclSVKbmM31uOckSZIkSRqQZ6IlSY3b4SVjkiS1itlcnY1oSVKjirEovfBJkqS2MJvrmcpGdES8HXg+cCfwDeC1mXlbRBwMXA1cU866ITNPnEglJUkl77uaB2azJE0Ts7mOad1znwUOy8yfB64F3tIx7RuZuap8GNKSJI2H2SxJmgtTeSY6Mz/T8XQD8NJJ1UWStDjHopwPZrMkTQ+zuZ5Z2HOvAz7d8fyQiLg8Ij4fEU+eVKUkSTttz6j10NQxmyWp5czm6lp7JjoiLgD26zHp5Mw8p5znZGAbsL6cdjOwMjNviYjDgU9ExKMy84c91n8CcALAygNauxskaeolYeclM8JslqTZYDbX09qEysyjFpseEa8Bngc8IzOzXOYO4I7y58si4hvAw4BLe6x/LbAW4IjH7J0jrbwkSTPIbJYkqcWN6MVExDHAm4CnZubtHeX7Ardm5vaIeDBwKHDdhKopSSrtsAfQmWc2S9J0MZurm8pGNPBe4O7AZyMCdg6X8RTgrRFxF7ADODEzb51cNSVJjkU5N8xmSZoSZnM9U9mIzsyH9ik/CzhrzNWRJGnumc2SpHkxlY1oSdL0SOzFU5KkNjGb67ERLUlqnGNRSpLULmZzdTaiJUmNyoTtdl4iSVJrmM31uOckSZIkSRqQZ6IlSQ0LduB9V5IktYfZXIeNaElSoxIvGZMkqU3M5npsREuSGudYlJIktYvZXJ17TpIkSZKkAXkmWpLUqCTY4ViUkiS1htlcj41oSVLjvGRMkqR2MZursxEtSWpUAjvsvESSpNYwm+txz0mSJEmSNCDPREuSGhZsdyxKSZJaxGyuw0a0JKlRXjImSVK7mM312IiWJDXOo92SJLWL2Vydhx8kSZIkSRqQZ6IlSY3KDC8ZkySpRczmetxzkqTGbc89aj0GERHHRMQ1EbEpIn6nx/T7RcQnI2JjRFwZEa8d+QuVJGlKmM3VeSZaktSoBHY0fN9VRCwD3gc8E7gRuCQizs3MqzpmeyNwVWY+PyL2Ba6JiPWZeWejlZMkqWXM5no8Ey1JmgWrgU2ZeV0ZvGcCx3bNk8B9IiKAewO3AtvGW01JkubGzGazZ6IlSQ2LgS/7quEA4IaO5zcCT+ia573AucC3gfsAv5SZO5qumCRJ7WM212EjWpLUqGIsytqXjO0TEZd2PF+bmWuHXMfRwOXA04GHAJ+NiH/NzB/WrZwkSdPEbK7HRrQkaRpsycwjFpl+E3BQx/MDy7JOrwVOzcwENkXEN4GHAxePtKaSJM2Huc1m74mWJDVuO3vUegzgEuDQiDgkIvYCjqO4PKzT9cAzACLiZ4CfA64b4cuUJGlqmM3VeSZaktSoJEZxydji28jcFhEnAecDy4APZeaVEXFiOf104I+AdRHxVSCAN2fmlkYrJklSC5nN9diIliQ1bscYLnzKzPOA87rKTu/4+dvAsxqviCRJU8Bsrs7LuSVJkiRJGpBnoiVJjcqE7Q1fMiZJkgZnNtdjI1qS1Lim77uSJEnDMZursxEtSWpU0XmJdw9JktQWZnM97jlJkiRJkgbkmWhJUuO24yVjkiS1idlcnY1oSVKjEu+7kiSpTczmemxES5Ia5n1XkiS1i9lcx1TuuYg4JSJuiojLy8dzOqa9JSI2RcQ1EXH0JOspSdK8MJslSfNims9EvzMz/6KzICIeCRwHPAp4EHBBRDwsM7dPooKSpMIO77uaF2azJE0Js7m6aW5E93IscGZm3gF8MyI2AauBf59stSRpfmXCdu+7mmdmsyS1jNlcz1Rezl1aExFfiYgPRcT9y7IDgBs65rmxLJMkTdCO3KPWQ1PDbJakKWE2V9faVx8RF0TEFT0exwJ/BTwYWAXcDLyjwvpPiIhLI+LS793iFWWSJC3FbJYkqcWXc2fmUYPMFxEfAP6xfHoTcFDH5APLsl7rXwusBTjiMXtn9ZpKkhaThMNozAizWZJmg9lcT2vPRC8mIvbvePoi4Iry53OB4yLi7hFxCHAocPG46ydJ2tUOotZD7Wc2S9J0MZura+2Z6CX8eUSsohgnfDPwawCZeWVEfBS4CtgGvNHePyVpshI82j0fzGZJmhJmcz1T2YjOzFcvMu1twNvGWB1Jkuae2SxJmhdT2YiWJE2Xee/FU5KktjGbq7MRLUlqVtp5iSRJrWI212IjWpLUqIS574BEkqQ2MZvr8Ry+JEmSJEkD8ky0JKlxXjImSVK7mM3V2YiWJDXKYTQkSWoXs7keL+eWJEmSJGlAnomWJDXOo92SJLWL2VydjWhJUqMSh9GQJKlNzOZ6bERLkhrnMBqSJLWL2Vyd90RLkiRJkjQgz0RLkpqV3nclSVKrmM212IiWJDXKYTQkSWoXs7keG9GSpMYZ1JIktYvZXJ33REuSJEmSNCDPREuSGuUwGpIktYvZXI+NaElS49KgliSpVczm6mxES5Ia51iUkiS1i9lcnfdES5IkSZI0IM9ES5IalY5FKUlSq5jN9diIliQ1zvuuJElqF7O5OhvRkqSG2QOoJEntYjbX4T3RkiRJkiQNyDPRkqTGecmYJEntYjZXZyNaktSoxM5LJElqE7O5HhvRkqRmZdELqCRJagmzuRbviZYkSZIkaUCeiZYkNW4HXjImSVKbmM3V2YiWJDUqsfMSSZLaxGyux8u5JUmSJEkakGeiJUkNC3sAlSSpVczmOmxES5IaZw+gkiS1i9lcnY1oSVLjvO9KkqR2MZur69uIjojrRrD+d2Xme0awHkmS5p7ZLEnS5C12Jvpg4AfAbRXXvRJYXnFZSdKMyPRo9wgdjNksSarJbK5nqcu535mZb62y4ojYUWW5Adf998DPlU+XA7dl5qqIOBi4GrimnLYhM09sqh6SpMHYeclImc2SpNrM5uqm8p7ozPylhZ8j4h0UR+UXfCMzV42/VpKkfuy8ZPaZzZI0Xczm6hZrRD8e+HaNddddfkkREcDLgac3uR1JklrCbJYkacL6NqIz87I6K667/ICeDHwnM7/eUXZIRFxOcQTAoFlkAAAgAElEQVT89zLzX3stGBEnACcArDxgKk/IS9LU8L6r0TCbJUmjYjZX19qEiogLgP16TDo5M88pf34F8JGOaTcDKzPzlog4HPhERDwqM3/YvZLMXAusBTjiMXt7MYMkNSQJg3pGmM2SNBvM5noqNaLLS7X2A+7Wa3pmXl+nUuU6jlqiDnsCLwYO71jmDuCO8ufLIuIbwMOAS+vWR5JUna2h5pnNkqRhmM3VDdWIjoiXAb8DHLbIsjnseis6CvhaZt64UBAR+wK3Zub2iHgwcCgwijE1JUlqJbNZkqTxGjhQI+KNwHuAbcAXgZvKnyflOHa9XAzgKcBbI+IuYAdwYmbeOvaaSZJ2cizKxpjNkqRKzOZahjkq/ZvAd4EnZeY3G6rPwDLzNT3KzgLOGn9tJEmLGsM1YxFxDPBuYBnwwcw8tcc8RwLvorjkeUtmPrX5mjXKbJYkVWM2VzZMI/oA4ANtCGlJ0nRp+mh3RCwD3gc8E7gRuCQizs3MqzrmWQ68HzgmM6+PiAc2WqnxMJslSZWYzdXtMcS8NwB3b6oikiTVsBrYlJnXZeadwJnAsV3zvBL4h4UOtjLzu2OuYxPMZklSW81sNg/TiP4w8OyIuE9TlZEkzabMeo8BHEDRoFxwY1nW6WHA/SPiwoi4LCKOH82rmyizWZJUidlc3TCXc58KPA64ICLeBPxHZv6omWpJkmZFMpJLxvaJiM4hkdaWYwoPY0+KoZeeAdwD+PeI2JCZ19at3ASZzZKkoZnN9QzciC6Hpngf8DHgnwGKISl7zZrjGEZDkjQNEqgf1Fsy84hFpt8EHNTx/MCyrNONwC2ZuRXYGhEXAY8BWh3UizGbJUmVmM21DDPE1bHAxyl6Vvsm8G0mO4yGJEkLLgEOjYhDKAL6OIr7rDqdA7w3IvYE9gKeALxzrLUcMbNZktRiM5vNwxyVPgW4HXhuZn6hmepIkmbRgPdO1Vh/bouIk4DzKRqUH8rMKyPixHL66Zl5dUT8E/AVivGKP5iZVzRbs8adgtksSarAbK5umEb0zwFnGNKSpKGNYSzKzDwPOK+r7PSu528H3t58bcbGbJYkVWM2VzZMI3oLcGdTFZEkzapofCzKOWY2S5IqMJvrGGaIq7OAZ0bE3ZqqjCRJGorZLEnSmA3TiP494PvAxyLi4EZqI0maTVnzoX7MZklSNWZzZcNczv1V4G4UPaY9PyJuA37QY77MzIeMonKSpBmQIxmLUr2ZzZKk4ZnNtQzTiN6DYtiM6zvKeu153w1JksbDbJYkacwGbkRn5sEN1kOSNMvm/LKvppjNkqTKzObKhjkTLUlSRZ4IlSSpXczmqmxES5Ka59FuSZLaxWyurG/v3BHx8oh4ZNUV111ekiTtymyWJGnyFhvi6kzgpTXWXXd5SdKscBiNUTGbJUmjYTZXttTl3HO+eyRJtSXgMBqjZDZLkuoxm2tZ7Ew0wCkRsb3KA0NeklTKrPfQLsxmSVJtZnN1i52Jvoj6Ybu55vKSJGkns1mSpAnr24jOzCPHWA9J0iyb8yPWo2I2S5JGxmyuzCGuJEnN874rSZLaxWyuzEa0JKlx4dFuSZJaxWyubqmOxSRJkiRJUskz0ZKkZjmepCRJ7WI212IjWpLUsPC+K0mSWsVsrsNGtCSpeR7tliSpXczmyrwnWpIkSZKkAXkmWpLUPI92S5LULmZzZUM3oiPigcARwP2BZb3mycwzatZLkjRLDOpGmc2SpKGZzZUN3IiOiLsBpwPH0/8y8KB4OwxqSVIhsfOShpjNkqRKzOZahjkT/UfAa4FvAOuBG4BtTVRKkiQNxGyWJGnMhmlEvxK4FnhsZv64ofpIkmZQeMlYU8xmSVIlZnN1w/TO/UDgvHGFdES8LCKujIgdEXFE17S3RMSmiLgmIo7uKD88Ir5aTntPRHiNgiS1QdZ8qB+zWZJUjdlc2TCN6OuB+zZVkR6uAF4MXNRZGBGPBI4DHgUcA7w/IhY6Ufkr4FeBQ8vHMWOrrSRJ42c2S5I0ZsM0otcBz46I+zVUl11k5tWZeU2PSccCZ2bmHZn5TWATsDoi9gfum5kbMnOhA5UXjqOukiRNyDrMZkmSxmqYRvSpwBeACyLiaRExziPfnQ6g6DhlwY1l2QHlz93lkqQJi6z3UF9msySpErO5ur4di0XEDnpf7R7ABeU8vRbNzByow7KIuADYr8ekkzPznEHWUVVEnACcALDygKGHy5YkDcNhNEbCbJYkjYzZXNliCXURDd8ynplHVVjsJuCgjucHlmU3lT93l/fb9lpgLcDKw+6bp922cpfpa5ZfX6FqkiQ1ymyWJGnC+jaiM/PIMdZjGOcCfxcRfwk8iKKTkoszc3tE/DAingh8CTgeOG2C9ZQkgb14jpDZLEkaCbO5ltZeKxURL6II2n2BT0XE5Zl5dGZeGREfBa4CtgFvzMzt5WJvoOhk5R7Ap8vHkm65816s37x6l7I1qzzaLUkjY1DPBLNZkmaI2VzZwI3oiLgOeFdmvmeRed4I/K/MfHDdimXm2cDZfaa9DXhbj/JLgcOG3da2n+zJlmtX7Fq4ati1SJL6mfcOSJpiNkuSqjKbqxumd+6DgeVLzLMc+NnKtZEkScM4GLNZkqSxGvXl3PcB7hzxOhu37Cew/OphjidIkobi0e5JMpslSbszmytbtBEdESu7ipb3KANYBqwEXgJcN6K6SZJmhUE9MmazJGkkzObKljoTvZldd+9vlI9+AvitmnUauz1v38GKjVsnXQ1JmkmR3nc1YpsxmyVJNZjN9SzViD6DIqiDYliKrwCX95hvO3AL8LnM/MxIayhJkjqZzZIkTdCijejMfM3CzxFxPHB2Zr616UqN3dYfw4aNk66FJM2ujEnXYGaYzZKkkTCbKxu4Y7HMnKvePU67bffby9Ysd3xKSarES8YaYTabzZJUmdlcWaXeuSPiIOCxwP2AHwBfzswbRlkxSdLs8L6r5pnNkqRhmM3VDdWIjohDgfcDT+8x7Z+BN2bmtSOq20St37x6t7I1qzzaLUlqF7PZbJYkjdfAjeiIeCjwb8AK4BvAF4D/BPYDfhF4BvCFiHhSZm5qoK6SpGnl0e5GmM2SpMrM5sqGORP9pxQh/RvA+zJzx8KEiNgDWAO8E/gT4OWjrOQkbLl2xe6Fq8ZfD0maeg6j0SSz2WyWpOGZzbUM04h+BnBeZp7WPaEM7XdHxNHAUaOqnCRpRhjUTTGbJUnVmM2VDdOI3ove41B2+jLw5OrVaY/lV89Vh6eSpOlkNkuSNGbDNKI3Ag9dYp6HAl+pXp32WLFx66SrIEmzw6PdTTGbJUnVmM2VDXNI90+AF0fEs3tNjIjnAi8C3jaKikmSZkdkvYf6MpslSZWYzdUNcyZ6BfBp4B8j4nPARcB3gJ8BnkoxtMYngX0i4vjOBTPzjNFUd4w2bJx0DSRJWorZLEnSmA3TiF5HcdI/KDoo6dVJyQuA53c8j3KZ6QtqSZLabx1msyRJYzVMI/q1jdViSpx228qe5WuWXz/mmkjSlJnzy74aZDabzZJUjdlc2cCN6Mz8cJMVmQbrN6/uWb5mlUEtSX1571RjzGazWZIqMZtrcawISZIkSZIGNMzl3ABExL7AS4BHAPfKzNd3lB8CfDUzfzzSWrbElmtX9J6warz1kKSp49HuRpnNPZjNkrQ4s7myoRrREfErwHuAvdnZMcnry8k/A/w7cALwNyOsoyRp2hnUjTGbJUmVmM2VDdyIjohnAmuBrwB/ABwNnLgwPTOviIgrgRcyo0G9/GqvfpekYQXed9UUs9lslqQqzOZ6hjkT/WbgZuCpmfnDiHhsj3m+AvzCSGomSZKWYjZLkjRmwzSijwDOzMwfLjLPjcB+9arUXis2bp10FSRpOnm0uylms9ksSdWYzZUN04jeC1gqqZYD26tXp+U2bJx0DSRp+oxpGI2IOAZ4N7AM+GBmntpnvsdT3Cd8XGZ+vPmaNcpsNpslaXhmcy3D3Ei0GTh8iXmeAFxTuTaSpNmUNR9LiIhlwPuAZwOPBF4REY/sM9+fAZ+p+YraYjNmsySpCrO5smHORJ8DvCkiXpaZH+ueGBGvBX4eOHlUlZsWp922creyNcuvn0BNJGlurQY2ZeZ1ABFxJnAscFXXfGuAs4DHj7d6jTGb+zCbJWniZjabh2lE/zlwHPCRiHgpcD+AiDgJeDLwYuDrwGmjrqQkaco1f8nYAcANHc9vpDgD+1MRcQDwIuBpTFFQL8FsliRVYzZXNnAjOjO/HxFPBc4AXtYx6T3l//8KvDIz566Hj/WbV+9WtmaVR7slacEI7rvaJyIu7Xi+NjPXDrmOdwFvzswdEVG7Qm1gNvdnNkvS4szm6oY5E01mXg8cGRE/TzFcxgrgB8CGzLysgfpJkmZB/aDekplHLDL9JuCgjucHlmWdjgDOLEN6H+A5EbEtMz9Ru3YTZDZLkioxmysbqhG9IDO/QjHupIAt167YvXDV+OshSXPsEuDQiDiEIqCPA17ZOUNmHrLwc0SsA/6x7SE9DLN5V2azJE3czGbz0I3oiPhZYF+KYxffK4+Az7XlVw/TybkkzZkBe/GstYnMbeV9wOdTDKPxocy8MiJOLKef3mwNJsts3p3ZLEmLMJtrGagRHRH7AL8LvAJ4YNe07wDrgT/NzFtHXkNJ0tQbx1iUmXkecF5XWc+AzszXNF+jZpnNkqQ6zObqlmxER8ShwGcprmcPYBtwS/nzA4D9gN8CXhIRRy10YV5XRLwMOAV4BLA6My8ty58JnArsBdwJ/HZm/nM57UJgf+DH5WqelZnfHUV9FrNi49z11yJJwxlDUM8Ts3lpZrMkLcFsrmzRa50iYg+KI9krgc8DRwH3zsz9M3M/4D7As4CLgIOB/zfCul1BMTTHRV3lW4DnZ+ajgV8G/m/X9Fdl5qry0XhIS5I0TmazJEmTtdSZ6GdR9Jj2UeAVmbnL8YrMvAO4ICI+B/w9xRHvZ2bmZ+tWLDOvBuju6jwzv9zx9ErgHhFx97Iuk7Fh48Q2LUnTYByXjM0Rs3kQZrMkLcpsrm6pRvRLgDuANd0h3Skzs7xp/AXASykuMRuHlwD/0RXSH46Iu4CzgD9erN5NOu22lT3L1yyf+75eJM0jg3qUzOaKzGZJ6mA2V7ZUI/pxwBcz83tLrSgzvxsRXyiXGUhEXEBx31a3kzPznCWWfRTwZxRH5Be8KjNvioj7UAT1q4Ez+ix/AnACwN7cc9AqS5KGNYYeQOeM2SxJqsdsrmWpRvRBwBeGWN+VFL2EDiQzjxpi3T8VEQcCZwPHZ+Y3OtZ3U/n/jyLi74DV9AnqzFwLrAW4bzxg5B+h9ZtX9yxfs8qj3ZKkWszmisxmSdIoLDWI4n2B24ZY320UHZo0JiKWA58Cficzv9hRvmc53AcRcTfgeRQdoEiSJihG8NAuzGZJUi1mcz1LnYneC9g+xPp2lMvUFhEvAk4D9gU+FRGXZ+bRwEnAQ4Hfj4jfL2d/FrAVOL8M6WXABcAHRlGXKrZcu6L3hFXjrYcktYKXjI2S2VyR2SxJHczmypYcJ5oJ7d7MPJvisrDu8j8G/rjPYoc3WilJktrBbJYkaUIGaUSfEhGnNF2RWbP86qWulJek+eEwGiNnNldgNkvSTmZzdYM0ooe95N23A1ixceukqyBJ7WEyjJrZXIHZLEkdTIbKFm1EZ6aHbCVJ9RnUI2M2S5JGwmyubJAz0apiw8ZJ10CSJHUymyVJI2AjWpLUrPS+K0mSWsVsrsVG9JiddtvK3crWLL9+AjWRpDEyqNViZrOkuWQ2V2YjWpLUOI92S5LULmZzdTaix2z95tW7la1Z5dFuSZImxWyWJA3DRvSYbbl2xe6Fq8ZfD0kaK492q8XMZklzyWyuzEa0JKlxXjImSVK7mM3V2Yges+VXO7ynpDmTeLRbrWY2S5o7ZnMtpoYkSZIkSQPyTPSYrdi4ddJVkKTx82i3WsxsljSXzObKbESP24aNuxUd/aDHcP63dy+XpFkQeN+VWs5sljRnzOZ6bES3gCEtaeYZ1JoyZrOkmWc2V2YjugVOu21lz/I1yx2jUpKkSTCbJUn92IiWJDUu0sPdkiS1idlcnY3oFli/eXXP8jWrPNotaQY4jIamkNksaaaZzbXYiJYkNc7OSyRJahezuTob0S2w5doVvSesGm89JElSwWyWJPVjI7oFll+9x6SrIEnN8mi3pozZLGnmmc2V2YiWJDXOS8YkSWoXs7k6G9EtsGLj1klXQZKaZVBrypjNkmae2VyZ1ypJkiRJkjQgz0S3wYaNk66BJDUnvWRMU8hsljTLzOZabERLkppnUEuS1C5mc2U2olvstNtW7la2Zvn1E6iJJEkCs1mSZCO61dZvXr1b2ZpVBrWk6RJ4yZhmh9ksaRaYzfXYiJYkNS9NakmSWsVsrsxGdIttuXbF7oWrxl8PSarLo92aFWazpFlhNlfnEFeSJEmSJA3IM9Ettvxqj3FImgGJPYBqZpjNkmaC2VyLjWhJUuNix6RrIEmSOpnN1dmIbrEVG7dOugqSNBoe7daMMJslzQyzubLWNqIj4mXAKcAjgNWZeWlZfjBwNXBNOeuGzDyxnHY4sA64B3Ae8BuZU9zt3IaNk66BJEk/ZTZjNkuS2tuIBq4AXgz8dY9p38jMXn1h/hXwq8CXKIL6GODTjdVQkjQQewCdGWazJM0Is7m61jaiM/NqgIgYaP6I2B+4b2ZuKJ+fAbyQGQvq025b2bN8zfLrx1wTSRpQ4liUM8Js7s1sljR1zOZaprWLyUMi4vKI+HxEPLksOwC4sWOeG8sySdKERdZ7aCqYzZI0Rczm6iZ6JjoiLgD26zHp5Mw8p89iNwMrM/OW8j6rT0TEoyps+wTgBIC9ueewi0/M+s2re5avWeXRbklSfWbz8MxmSZovE21EZ+ZRFZa5A7ij/PmyiPgG8DDgJuDAjlkPLMv6rWctsBbgvvGAqTmWsuXaFb0n9LoLTZLaYmq+ZWU2D89sljSVpuZbtn2m7nLuiNg3IpaVPz8YOBS4LjNvBn4YEU+M4mat44F+R8wlSWMSeMnYrDObJWm6mM31tLZjsYh4EXAasC/wqYi4PDOPBp4CvDUi7gJ2ACdm5q3lYm9g5zAan2bGOi4BWH711B33kDTvMu28ZEaYzb2ZzZKmjtlcS2sb0Zl5NnB2j/KzgLP6LHMpcFjDVZMkaS6ZzZIktbgRrd5WbNw66SpI0tDm/bIvzTazWdI0MpursxEtSWqeQS1JUruYzZXZiJ42GzZOugaSNDSPdmummc2SppDZXJ2N6Blx2m0rdytbs9zxKSVJmhSzWZJmk41oSVKzEtjh4W5JklrDbK7FRvSMWL959W5la1Z5tFtSS5jTmkNms6RWM5src2BDSZIkSZIG5JnoGbHl2hW7F64afz0kqRc7L9E8MpsltZnZXJ2NaElS89KkliSpVczmymxEz4jlV3tlvqT2GsfR7og4Bng3sAz4YGae2jX9VcCbgQB+BPyPzHRsIjXGbJbUZmZzdTaiZ8SKjVsnXQVJmpiIWAa8D3gmcCNwSUScm5lXdcz2TeCpmfn9iHg2sBZ4wvhrq3lhNkuaZ7OczR4ilSQ1K0fwWNpqYFNmXpeZdwJnAsfuUo3Mf8vM75dPNwAH1nlZkiRNLbO5Fs9Ez4oNrb/qQdKcCiCav+/qAOCGjuc3sviR7F8BPt1ojSSzWVJLmc312IiWJDVvR+017BMRl3Y8X5uZa6usKCKeRhHUv1i7VpIkTSuzuTIb0TPstNtW9ixfs/z6MddEkmrbkplHLDL9JuCgjucHlmW7iIifBz4IPDszbxltFaWlmc2SZsjcZrON6Bm2fvPqnuVrVhnUksZrDJeMXQIcGhGHUAT0ccArd6lDxErgH4BXZ+a1TVdI6sVsltQWZnN1NqIlSc0avAOS6pvI3BYRJwHnUwyj8aHMvDIiTiynnw78PrACeH9EAGxb4gi6JEmzyWyuxUb0DNty7YreE1aNtx6S5l1C80e7yczzgPO6yk7v+Pn1wOsbr4i0CLNZUjuYzXU4xJUkSZIkSQPyTPQMW361x0gktUM0f7Bbmgpms6S2MJursxEtSWreGC4ZkyRJQzCbK7MRPcNWbNw66SpIEiRE/bEopZlgNktqBbO5FhvRs2zDxknXQJIkdTKbJWnq2YiWJDXPS8YkSWoXs7kyG9Fz6LTbVu5Wtmb59ROoiaS5YU5LizKbJY2d2VyZjWhJUuPCo92SJLWK2Vydjeg5tH7z6t3K1qzyaLckSZNiNkvS9LARLUlqnke7JUlqF7O5MhvRc2jLtSt2L1w1/npImhMJOIyGtCizWdJYmc212IieQ8uv3mPSVZA0R4L0vitpCWazpHEym+vxG1uSJEmSpAF5JnoOrdi4ddJVkDRvPNotLcpsljR2ZnNlNqIlSc0zqCVJahezuTIb0fNow8ZJ10CSJHUymyVpatiIFgBHP+gxk66CpFb4+uhXaQ+gUiVms6SC2dw2re1YLCJeFhFXRsSOiDiio/xVEXF5x2NHRKwqp10YEdd0THvg5F6BJGlBZNZ6qB3MZkmaHWZzdW0+E30F8GLgrzsLM3M9sB4gIh4NfCIzL++Y5VWZeenYailJWtqch+0MMZslaVaYzZW1thGdmVcDRMRis70COHMsFZIkac6ZzZIktbgRPaBfAo7tKvtwRNwFnAX8cWbvQywRcQJwAsDe3LPRSkrSfEuPds8Xs1mSWs9srmOijeiIuADYr8ekkzPznCWWfQJwe2Ze0VH8qsy8KSLuQxHUrwbO6LV8Zq4F1gLcNx7gJ0iSmpIY1FPEbJakOWA21zLRRnRmHlVj8eOAj3St76by/x9FxN8Bq+kT1JKkMbIH0KlhNkvSnDCbK2tt79yLiYg9gJfTcc9VROwZEfuUP98NeB5FByiSJKlhZrMkaV60thEdES+KiBuBXwA+FRHnd0x+CnBDZl7XUXZ34PyI+ApwOXAT8IGxVViS1JfDaMwGs1mSZofZXF1rOxbLzLOBs/tMuxB4YlfZVuDw5msmSRranIftrDCbJWmGmM2VtbYRLUmaEQnsMKglSWoNs7mW1l7OLUmSJElS23gmWpLUMMeilCSpXczmOmxES5KaZ1BLktQuZnNlNqIlSc0zqCVJahezuTLviZYkSZIkaUCeiZYkNcseQCVJahezuRYb0ZKkhiXkjklXQpIk/ZTZXIeNaElS87zvSpKkdjGbK/OeaEmSJEmSBuSZaElSs7zvSpKkdjGba7ERLUlqnpeMSZLULmZzZV7OLUmSJEnSgDwTLUlqnke7JUlqF7O5MhvRkqSGpUEtSVKrmM112IiWJDUrgR2ORSlJUmuYzbV4T7QkSZIkSQPyTLQkqXleMiZJUruYzZXZiJYkNc+gliSpXczmymxES5IalrDDoJYkqT3M5jq8J1qSJEmSpAF5JlqS1KyETHsAlSSpNczmWmxES5Ka5yVjkiS1i9lcmY1oSVLz7LxEkqR2MZsr855oSZIkSZIG5JloSVKzMmGH911JktQaZnMtNqIlSc3zkjFJktrFbK7MRrQkqXHp0W5JklrFbK7Oe6IlSZIkSRqQZ6IlSQ1LLxmTJKlVzOY6bERLkpqVOBalJEltYjbXYiNaktS89L4rSZJaxWyuzHuiJUmSJEkakGeiJUmNSiC9ZEySpNYwm+tp7ZnoiHh7RHwtIr4SEWdHxPKOaW+JiE0RcU1EHN1RfnhEfLWc9p6IiMnUXpL0U5nFJWN1HgOIiGPKXNgUEb/TY3qU2bCpzJbHjfy1zjizWZJmhNlcS2sb0cBngcMy8+eBa4G3AETEI4HjgEcBxwDvj4hl5TJ/BfwqcGj5OGbclZYk7S53ZK3HUsoceB/wbOCRwCvKvOj0bHbmwwkUmaHhmM2SNCPM5upa24jOzM9k5rby6QbgwPLnY4EzM/OOzPwmsAlYHRH7A/fNzA2ZmcAZwAvHXnFJ0iSsBjZl5nWZeSdwJkVedDoWOCMLG4DlZXZoQGazJGkIM5vNrW1Ed3kd8Ony5wOAGzqm3ViWHVD+3F0uSZq05i8Z65cNw86jwZnNkjTNzObKJtqxWERcAOzXY9LJmXlOOc/JwDZg/Yi3fQLFJQMA/3VBfvyaUa6/5fYBtky6Ei3gfii4Hwruh8LPjXqFP+L751+QH9+n5mr2johLO56vzcy1NdepHszmifE7qOB+KLgfCu6HgtncMhNtRGfmUYtNj4jXAM8DnlFeBgZwE3BQx2wHlmU3sfOyss7yftteC8zFm9wtIi7NzCMmXY9Jcz8U3A8F90OhKwxHIjPHcQ9sv2wYdp65ZzZPht9BBfdDwf1QcD8UzOb2ae3l3BFxDPAm4AWZeXvHpHOB4yLi7hFxCMVN6Bdn5s3ADyPiiWXPn8cD54y94pKkSbgEODQiDomIvSg6uTq3a55zgePLnkCfCPygzA4NyGyWJA1hZrO5zeNEvxe4O/DZcjSMDZl5YmZeGREfBa6iuJTsjZm5vVzmDcA64B4U92l9ere1SpJmTmZui4iTgPOBZcCHyrw4sZx+OnAe8ByKTq9uB147qfpOMbNZkjSQWc7m2HklluZFRJwwL/crLMb9UHA/FNwPBfeDNBn+7hXcDwX3Q8H9UHA/tI+NaEmSJEmSBtTae6IlSZIkSWobG9FzJP4/e3ceLVtZ3nn8++MyOAYUFJAhYMQBTURFNIkDjgyJonEI2B2GGNFW0MSkFZpkBU0wxDiPhDY0kEYRpY2oGAIoIonIoIACAle8yqQICIgDcO99+o+9D9Q9t85Q06k653w/rr2qau/97v1UeRfPed/9DslxSW5O8t1xxzJOSbZL8tUkVyS5PMlbxh3TOCR5QJILklza/g7vGHdM45RkRZJvJ/niuGMZlySrknwnySWjmAlU0vrMzQ1zc8PcvC5zs7l5UtmdexlJ8hzgLuDEqnrSuOMZl2wPR0AAACAASURBVCRbA1tX1beSPBS4GHhZVV0x5tAWVDtT7oOr6q4kGwHnAW+pqvPHHNpYJHkrsCvwG1X1h+OOZxySrAJ2rSrX5JQWiLm5YW5umJvXZW42N08qn0QvI1V1LnDbuOMYt6q6qaq+1b7/OXAlsM14o1p41bir/bhRuy3LVrUk2wJ/AHxi3LFIWl7MzQ1zc8PcfD9zsyaZlWgta0l2AJ4CfHO8kYxH203qEuBm4MyqWpa/A/ABmrVv1447kDEr4KwkFyc5eNzBSFqezM3m5pa5uWFunkBWorVsJXkIcCrw51V157jjGYeqWlNVuwDbArslWXZdCZP8IXBzVV087lgmwLPafw97AW9qu5lK0oIxN5ubwdw8jbl5AlmJ1rLUjjM6FTipqv7fuOMZt6q6HfgqsOe4YxmD3wde2o45Ohl4fpL/O96QxqOqbmhfbwY+B+w23ogkLSfm5nWZm83NYG6eVFaitey0k3b8C3BlVb1v3PGMS5JHJNmsff9A4EXA98Yb1cKrqsOratuq2gHYF/hKVf33MYe14JI8uJ3MhyQPBl4MLOvZgiUtHHNzw9zcMDc3zM2Ty0r0MpLkU8A3gMcluT7Ja8cd05j8PvAnNK2al7Tb3uMOagy2Br6a5DLgQppxV8t2CQmxJXBekkuBC4AvVdW/jzkmackzN9/H3NwwN6uTuXlCucSVJEmSJEnz5JNoSZIkSZLmyUq0JEmSJEnzZCVakiRJkqR5shItSZIkSdI8WYmWJEmSJGmerERryUlyTpJFM+18kpq2bdVD2VVJVo0wvIEk+eK073bguGOSJC08c/PkMDdLg9tw3AFIM+kj2R5UVcePIpYF8EPg+Pb9XWOMY9g+CVwE7ALsM+ZYJEkDMjcvCeZmaUBWojXJ3tFl358DmwIfBG6fduyS9nV/4EEjjGsUVlXVkeMOYtiq6pMAbSu3iVqSFj9z8yJnbpYGZyVaE6tb4mr/g78p8IGqWjVDuR+NNDBJkpYpc7MkOSZaS1C3cVdJdm/H/RyZZNck/57kjiQ/S3Jqku3a8x6d5OQkP03yqyRfTfLkGe7zoCSHJ7kkyS+S3JXkG0n2G/L3SZJDklye5NdJbkjykSSbznD+pkn+Z5KvJLk+yT3t9zktye9OO/dhSX6Z5PtJMsP1vtD+drt27HtpkrOT3JTk7iQ3JvlakjcO87tLkpYGc7O5WVpKrERruXk68PX2/f8GLgD+CDgryePbz9sCJwJfAp4LnJnkIZ0XSbIZcB7wLmANcBxwAvAI4JNJ/n6IMX8A+DDwMOBY4GRgT+AsYOMu5z8BOApY236H9wFnAs8Hzk2y59SJVfWz9nqPBl44/ULtHzB7ARdX1UXtvoOBzwM7A18A3gucDjwQOGjgbytJWm7MzeZmaXGpKje3RbMBq4ACdpjlnHOaf9rr7Nu9LVfAf5t27F/a/bcBR0w79jftsbdM2398u/9t0/Y/APh3miS5yzy/UwHnzHDs99rjK4GHT7vPN9pjq6aV2RTYosu1tgVuBK6ctn/X9jqf7VLmyPbY6zr2XQzcDTyyy/nr3bfdf2B7nQPH/W/Izc3NzW24m7l5nfuYm93clsHmk2gtN+dV1UnT9p3Qvt4BHD3t2Int6y5TO5JsDvx34KKqenfnyVX1a+DtQIDXDCHeqdbjo6rqtmn3Obxbgaq6o6pu6bL/euCzwOOTbN+x/yKaWTr3SccSHklWAK8Ffg58atrlVgP3drnHeveVJGkO5mZzs7SoOLGYlpuLuuy7sX29pKrWTDt2Q/u6bce+pwMrgEpyZJfrbdS+PqHfIDs8tX39Wpdj59F0V1tPkt8H3gL8LvBI1u9atg3QOcnLx2i6vf0pTTc4gL1pvvfHq6pzaY+TaLqJXZHk5Da2/6yqn87zO0mS1Mnc3DA3S4uElWgtN3d02bd6pmNVtbqd02Ojjt2bt69Pb7eZPGSWY/M1NUHJT2aIbb3W5SQvp2nV/jXNeKvvA7+g6ca2O81Ysk2mFTuZJvm+LsnRVbUWOLg99s/T7vu+9r5vBN5Ms7RJJfka8D/b1nNJkubL3GxulhYVK9FS76YS+vur6q0LdK8tgWs7DyTZENgCuH5amb8D7gF2raorp5X5Z5pEvY6q+lWS44G/AF6c5HKaSUu+WVWXdjn/RODEdhKX3wNeTtNSfkaSx9vyLUlaYOZmc7O0YBwTLfXuApqW42cvwL2+1b6ul1yBZ9F0XZvuMcAVXZL0Bm2ZmXycZoKR19OMt1rBtJbu6arq9qo6vapeRzOhy8OB58xWRpKkETA3t8zN0uhZiZZ6VFU304w92jXJ37STfKwjyW8l2XEItzu+fT0iycM7rv8A4B9mKLMK2CnJozrOD81snjvPdKOqugY4G/hD4A3A7TRdydaR5HkzrFv5yPb1lzPdQ5KkUTA3m5ulhWR3bqk/hwA7Ae8E/iTJeTRjox5FM2nJ04H9gB8McpOq+s8kHwYOBb6b5LM0M2/uA/wMuKlLsfcDxwDfTnJqe/7vc//akS+Z5ZYfo1mTckvgw1X1qy7nfA64K8n5NH8UhKbl/+k0S2yc1ePXlCRpGMzN5mZpQfgkWupDVd1J043rUOAW4BXAW4Hn0Sw78Rc0E4cMw1va+9xB051rP+AMmoR6T5fY/plm+Y2bgAOA/wZcBzyD+7ugzeQ0mu8DM3cXOwy4kGZ20je299qIZvmQ51XVestrSJI0auZmc7O0UFJV445BWtaSFPC1qtp9AmJ5NLCSZlmMoY0rS3Ig8H+Ag6rq+GFdV5KkUTA3S5qN3bmlyfDcNmEDbF1VPx5THH9F0wXsI8O4WJIvAn8wjGtJkrTAzM2SurISLY3fO6Z9vmshb55ke+A1NOPIDgIuBT4zpMt/Euhcm/KSIV1XkqRRMjdLmpHduaVlLsnuwFdpZu48D/gfVXXtrIUkSdLImJulyWYlWpIkSZKkeXJ2bkmSJEmS5slKtCRJkiRJ82QlWpIkSZKkebISLUmSJEnSPFmJljRUSQ5PcmGSO5P8NMkXkjxpAuJ6TpLTktyQpJIcOO6YJEmStPhYiZY0bLsDHwN+D3g+sBo4K8nDxxkU8BDgu8BbgF+NORZJWpYmtaF1qfF3lkbLSrRGIskJSW5O8uAxxvC09onjn40rhlFIskP7vY6f5/m7t+dPbd8bZXxVtUdV/Z+q+m5VfQf4E+ARwO+P8r7ziOv0qvpfVfVZYG23c5JsMe23cg1ASRqu3ZnMhtalZnf8naWRsRKtoUvydJqK09FV9YsFuN9ftBWe13Tur6qLgX8D/i7JQ3q4Xk+V1IW+3gC+BrwD+MgC3/ehNP+t+dkC37cfv6T5jd4B/HDMsUjSktNrQ+skNMq3cfTVML/QDdlT5vM723As9c9KtEbhKOBO4OMLdL+nta8Xdzn2D8BWwJsXKJZJdk5VHVlVC12J/iBwCfCNBb5vz6rql+1vdCSwaszhSNLE6qh4/TDJA2Y4Z1V7zoazXGrGhtaFbpRv7znUhvkO42rIntLtd7bhWOqTlWgNVZLHAi8ETqmqhRp3+jTgLuDq6Qeq6gLge8Drk/jvfYEleR/wLOAVVbVmyNf+++kt6F223Yd5T0nSerYH/nyA8rM1tC50ozyMrmF+XA3ZU9b7nW04lvpnpWIZS/KEtqLxuWn7n9HOYPzDJE/p8bJ/CgT49Az3PG+OSs+5PcR/dNv16PE0k0at7bjOn3ScejJNkn/RPK55JPCD9uMB02I7cNq5r05ybpI7kvwqyXfaiTw26fV6SQ5McmqSa9tr3ZnkP5P89/n+Hv3q7G6e5LeSfDbJrUl+nuQ/0k5EkuQRSY5NclOSX6eZsOR5s1z3/cB+wPOr6toeY3ppkrPbe92d5MYkX0vyxo7TPgA8YY7tgl7uK0nqyc+A24DDkmzRa+HZGlrH1CgPE9owP0jD8SgbtKXlarbuNVr6ntq+fmtqR5IDgH8GzgdeVVU/7fGaLwTWtOW7+TfgrC77D6Kp6H6lh3t9CzgBOAD4L+DMjmPndLz/z/b1RcAZc1zzHGAzmhmcL23jnXLJ1Jsk7wIOB24BPkmTcPcC3gXskeTFVXXPfK9H08p+OXAucBOwObA38K9JHldVfzNH3MOwA/BN4Erg+Pbzy4Fzkvwu8O80TwQ+DTwc2Bf4cpLHVtWPOi+U5IPAHwPPq6qexn8lOZjm3+CPgS/Q/MaPBH6H5t/JxwCq6pb2mCRpPH4JvAd4P/C3wKHzLdg2tO5Lkye6NbTO2SjP7BNWfr2qntNDPEcDb+/YtTbJ1Pv9q+pf2/cnA0cyv78p5nPfl9L8jbAzTW69FbgG+HRVfazj1A8A/3eOy/1o+o55/M6S+mAlenm7rxKdZAVNIvxzmkrKn1fVvb1cLM2kH7sAV840dqmq3tOl3HtoKtD/B3jnfO9XVack2YymEn1CVR07w6kXtq9zJtOqOifJKpqEdknbxWl6vL9LU4G+Dtitqn7c7j8c+Bzwh8BfAe+az/VaT6qq70+7z8bAl2la+I+pqhvmin9AzwX+uqqO6ojhb2j+P/kmcArwxqpa2x47EzgR+It2myrzUZoxbC8DfpZkq/bQXVV11zzieD1wD/Dkqrq580A/Tzo6yj4EeEz7cQNg+yS7ALdNbwSQJM3bR4FDaJ7OfqiqrpmrwDwbWheyUR5G0zA/q/k2GkN/DceDNGhLmp2V6OVtqhK9iuYp43OA11XVJ/q83jbACponqXNK08T7UeB/tK+HVlWvM0Ou9zR9uqq6I8mvaZLqMPxp+/r3UxXo9j6rk/wlzRPkP6N5Kj0v0yvQ7b572grp84EX0FRYR2kVcPS0fSfQVKI3Af7nVAW69UngOJqGk05TXa7Pnrb/HTSt9/OxGlivEaf9I6JfuwJfnRbPO2i+44EDXFeSlq2qujfJYcBngH8E/mi28+fT0LrQjfLt9YbeMD8PI2k0bssP2qAtaRaOiV7edqH5j/fnaVp8v91LBTrJG5J0tjhv3r7OuZRR++T7eJoK9Lur6pA+KtDQVKLvBb4zx3m3AQMlpGn3hC6t3FV1NXA9sGOSTed7wSTbJ/loku8l+eXU+Cbg1PaUbQaOem6XdBkrdWP7enVV/bzzQHvuT4Btp+3PDNuR84zjJOBBwBVJ3p/kZUke0fvXWVdVnTNDXAcOem1JWs6q6rM0E1a9PMmz5jj9jTQzRZ9N0+g+tf1Vxzk9N8on+RjwlzSN8q+d1ug7X/NqmAeG2TA/ikZjmN/vLKlPPoleppI8mmasLsDNNOOanpHkJVX1hXleZhfWHdc7NfFH16UuOu69Ec1TzFcCf1tVPbUWd1xnQ+C3gSuq6u45Tn9gR3yDmqocz5Tcb6JJrpsBd8x1sfb/iwuAhwFfB/6jLbeGZlzyATRPgkdtvVjbp+tdj7VWAxsNM4iqel+SW2j+AHgzzRCDSvI1mqfhFw3zfpKkofhLmm7Q7wGeOdNJVZWZjnXotVH+OGB/mkb5t89RZDa9NMxvOcB9ppwEvJem0fhkmmWw/rOP+WjWM8/fWVKffBK9fE21tn4O2B04ov18VA8zTu4CfLvj81RXpM27nAtAmrUkP0dTgf6rfivQrZ1pKuwzthi399yApkJ782zn9WCqQrnVDMe3nnbeXN5K85u9tqp2r6o3V9XftE9uB560ZDGqqhOr6pk0v8sfAP9C03XujGE8lZYkDVdVfQP4LE2D/B8PeLleGuVPpqlA/+0gFehxNMxX1ftoGsp/SNNo/DngJ0m+mmTXQa8vaXSsRC9fU5Xok6rqnvbp8zdpEsh6yyol2aX9j/rUUk67ted2Pom+Cfgp8LhuN2zHOH2JZszwG6vqvQN+h6mxuN+e9awmnrBurLOZ6tK8YobjU/fbffqBJI+h6d78g6q6fZ7Xm5rs6tQux547a6RLXFXdXlWnV9XraLr/P5zhjEOTJA3f4TRPcv+hnRyzXwvdKA9japi30VhanKxEL19TlejOCuhft6/v6Ex+SXai6WI0Vck+jGam5gfRUTFtxzSfC2zRVibpuMamNN2UnwscWFUfny24NGsWr7c28zRTyfXO2a7F/d3KvjrrWff7GVDMPN7puPb1rzsTXMcM5xvQJMH5Xm9V+7p7584ke9BMULasJHleOtYV6fDI9vWXCxmPJGl+qmolzYzSO9LDclddLHSjPIyuYX5ebDSWFhfHRC9fT6HpbvyDqR1VdVaSc2gqc/8D+GB76CPAl6rqsPbzyiQvA15aVTeyrlOBVwB7ACs79n8S+D2asb+PTnJkl5j+oaML1VQDz+pZvsPF7etRSZ4E/AK4vKo+M+28F9M8Df78LNe6T1XdleSbwLOTnARc3ZY/raouq6r/SvJu4G3Ad5N8tr33XsCTgPOAf5rv9Wj+4DgI+Ex7rRvb6+xJ01gxaLe4xeZzwF1JzqdpYAjwbODpNP+fd1vSRJI0Gd5J00X5CKCfyb2oqkpyLvCKJI9pK+fAfY3ypwPPoGmUn3XliiTHt/EcVFXHz3LqqBrmZ5TkecA5XSZWtdFYmnBWopehJNvQ/Af6a13+w30EzRqIRyQ5jqbL0ouB3aaddy/dW2FPpenitD/NDJlTXZ+mWlN363ItgJunzd7828DPaVqau6qqc5McStPafSjN5FvvollmY+q7bkqzvMMXq+q6ma7VxZ8A76epyO5HU5G7Hrisvffbk3ybZm3M/Wkm1/o+zdP891bVPfO9XlVd1ibSv6fpyrUhcCnNMiG3s/wq0YfRNMI8leYpw69pxou9Hfh49bh+uSRp4VTVbUneBbx7wEstZKM8jKhhfg42GkuLVPpbVUjLRZKXAv8PeFBnxbD9D/7Xuk3ikeRwmsrsU6tqrm5R3e65GXArTWX0bX0H31zrUOBDwLOr6rxBrrVYJdmdpsX8HT0sM7Xstb0ynusMp5K0vnYZxhuqatsuxzYBvkezwgTARlU1VyV2+jU2Bq4DVlXVM9p9G9D0onvILEVvrqr7Zs5uG7x/C/jNqpp1tu8kh9A0yv8mbcN8VR3RcXxT4MfAGVX1snl+j92ZIQcneQNNI8GTaSYrnWo0/hRNo/E6S0uOkjlP6o1jojWXopkQ674llpI8m6Yb1UwV5PcDP6Lp0tWPZ9M86X5fn+UBSPJAmklOTl2uFehp/rYdZ/69cQcyqZJs0bFG97Ke1E2SZlNV6VaBbo/dXVU7tuek1wp0e417gA8AuyV5SrtvbVU9tOO63bbOCvRmwO8Ax8xVgW6v/5GqelxVPaC91hHTTtmfZvKx9/T6fWa43zFV9fKqenRVPaiqHl5VT6mqdy9kBVpS73wSrVkl2Zpm3PQJNF2zHk8zhnd74AlV1bVCluQ5wPOA91TVLxYo3OkxPIGmK/TxVbVqHDFMgiQ7AAd27Lqlqj4ylmAmXJIH0Yx1v49P7yVpPNoZuK+iGfr0kj7Kv4RmiNcOVfXjAWN5IM2wrf+qqlf2UG531h0/fVVVPX6QWIYlyRY0E7jdxyfR0vwsyUp0kj1pJsVaAXyiqo4ec0iLWpL9gKNpxlF/GzgT+CvgoVXV16QhkjRM7RwOf0jTlfNJXY6HJi/sTTNZz4FVNetSNhouc7P6MQmN8m0cfTXMT3JDtg3HUv+WXCW6XWboauBFNBM3XQjsV1VXjDUwSdLItH9o3wWcOEMlem+asY570wxH+eDUOEuNnrlZkrSULMUx0bsBK6vq2nY8zcnAPmOOSZI0QlV1LnDbLKfsQ1PBrqo6H9isHa6ihWFuliQtGUuxEr0NzWyOU65v90mSli9zw3j5+0uSloxlu050koOBgwEe/KA87fGP2XjMEUnS+F182d23VNUjhnnNPZ734Lr1tjUDXePiy+6+nGb5lynHVtWxA11UE6czNz8Anrb9mOORpElwdTOWfqi5+TFJ/XLAa9zULPe251ACWmSWYiX6BmC7js/btvvW0f7xdSzArk9+QF1whqlaklZsfc0Ph33NW29bw6D/jV2x9TW/rqpdB7jEvHKDRqbn3Py4pP55YWKTpIn2vGb98KH6JfD6Aa9xJGwxhFAWpaXYnftCYKckOybZGNgXOG3MMUnSslXA2gH/NwSnAfun8Uzgjqq6aRgX1ryYmyVpgoTmaeog23K25L5/Va1OcghwBs0yGsdV1eVjDkuSlrFizYhXw0vyKWB3YIsk1wN/C2wEUFXHAKfTzMy9kqYB/qCRBqR1mJslabKENkmqL0uuEg1QVafT/MEkSRqz5kn0aJdTrKr95jhewJtGGoRmZW6WpMkx9SRa/VmK3bklSZIkSRoJGyAkSSM3pHHNkiRpCOzOPRgr0ZKkkSqKNTXa7tySJGn+7M49GH87SdLIjXpMtCRJmj+fRA/GMdGSJEmSJM2TT6IlSSNVwBqfREuSNDHszj0YfztJ0sjZnVuSpMlhd+7B2J1bkiRJkqR58km0JGmkCpydW5KkCWJ37sH420mSRs5VoiVJmhx25x6MlWhJ0kgV5cRikiRNECvRg3FMtCRJkiRp6JLsmeSqJCuTHNbleJJ8qD1+WZKnzlU2ycOTnJnkmvb1YR3HDm/PvyrJHh37n5bkO+2xDyXJIN/LSrQkabQK1gy4SZKk4dpwwG0uSVYAHwX2AnYG9kuy87TT9gJ2areDgY/Po+xhwNlVtRNwdvuZ9vi+wBOBPYGPtdehve7rOu615zy+woysREuSRqpoxkQPskmSpOGZ6s49yDYPuwErq+raqroHOBnYZ9o5+wAnVuN8YLMkW89Rdh/ghPb9CcDLOvafXFV3V9UPgJXAbu31fqOqzq+qAk7sKNMXx0RLkkYsrGGgXlOSJGmIhjQ79xZJLur4fGxVHdvxeRvguo7P1wPPmHaNbudsM0fZLavqpvb9j4EtO651fpdr3du+n76/b1aiJUmSJEm9uqWqdh1nAFVVSRZ84JeVaEnSSBWw1nHNkiRNjAWanfsGYLuOz9u2++ZzzkazlP1Jkq2r6qa2q/bNc1zrhvb9bHH0xDHRkqSRW9N26e53kyRJwzPVnXuUE4sBFwI7JdkxycY0k36dNu2c04D921m6nwnc0XbVnq3sacAB7fsDgM937N83ySZJdqSZQOyC9np3JnlmOyv3/h1l+uKTaEnSSBVYEZYkaYIsxJPoqlqd5BDgDGAFcFxVXZ7kDe3xY4DTgb1pJgH7JXDQbGXbSx8NnJLktcAPgVe3ZS5PcgpwBbAaeFNVrWnLvBE4Hngg8OV265uVaEmSJEnS0FXV6TQV5c59x3S8L+BN8y3b7r8VeMEMZY4Cjuqy/yLgSb3EPhsr0ZKkkVtbPomWJGlSDGl27mXL306SNFJ255YkabIs0MRiS5aVaEnSSBVhjfNYSpI0MXwSPRj/qpEkSZIkaZ5sgJAkjZxjoiVJmhx25x6MlWhJ0kg5JlqSpMliJXowVqIlSSMW1pSjhyRJmiRWBPvnXzWSJEmSJM2TDRCSpJEqYK1ttpIkTYwAGw1aE1w9jEgWJyvRkqSRc0y0JEmTI4ENrUT3zUcDkiRJkiTNk0+iJUkjVeXEYpIkTZIENlox7igWLyvRkqSRW2t3bkmSJsZQunMvY/50kqSRataJ9km0JEmTYigTiy1ji/KvmiT/lOR7SS5L8rkkm7X7d0jyqySXtNsx445VkqTlwNwsSVouFmUlGjgTeFJV/Q5wNXB4x7HvV9Uu7faG8YQnSbpfMyZ6kE2LgrlZkhaLACsG3JaxRfkQv6r+o+Pj+cArxxWLJGl2rhO9PJibJWkRCYu0JjgZlsJfNX8KfLnj845td7GvJXn2uIKSJN1vTWWgTYuOuVmSJtlUJXqQbRmb2K+f5Cxgqy6Hjqiqz7fnHEGzzPdJ7bGbgO2r6tYkTwP+LckTq+rOLtc/GDgYYPttJvZnkCRpYixkbt5yFF9AkqQhmNjaY1W9cLbjSQ4E/hB4QVVVW+Zu4O72/cVJvg88Frioy/WPBY4F2PXJD6ihBi9Juk8RZ+deIhYyNz8uMTdL0ihNbE1w8i3Kny7JnsDbgOdW1S879j8CuK2q1iR5NLATcO2YwpQktdY6OdiSZ26WpEVkamIx9WVRVqKBjwCbAGcmATi/ne3zOcA7k9wLrAXeUFW3jS9MSZLrRC8b5mZJWiycWGwgi/Knq6rHzLD/VODUBQ5HkqRlz9wsSVouFmUlWpK0eBTOsC1J0kTxSfRA/OkkSSPnOtGSJE0Yx0T3zUq0JGmkqmCNE4tJkjQ5fBI9EP+qkSRJkiRpnmx/kCSNWFiLY6IlSZoYPokeiD+dJGmkCrtzS5I0cRwT3Tcr0ZKkkXOdaEmSJohPogfiXzWSJEmSpAWT5OFJzkxyTfv6sBnO2zPJVUlWJjlsPuWTHN6ef1WSPTr2Py3Jd9pjH0qSdv+BSX6a5JJ2+7O54rcSLUkaqSKsrcE2SZI0RFNPogfZBnMYcHZV7QSc3X5eN8RkBfBRYC9gZ2C/JDvPVr49vi/wRGBP4GPtdQA+DrwO2Knd9uy43aerapd2+8RcwVuJliSN3Bo2GGiTJElDNP5K9D7ACe37E4CXdTlnN2BlVV1bVfcAJ7flZiu/D3ByVd1dVT8AVgK7Jdka+I2qOr+qCjhxhnvOiz3hJUkjVcBaJxaTJGmyjHdisS2r6qb2/Y+BLbucsw1wXcfn64FnzFF+G+D8aWW2Ae5t30/fP+UVSZ4LXAX8RVV13nc9VqIlSZIkSb3aIslFHZ+Prapjpz4kOQvYqku5Izo/VFUlqX6DGLQ88AXgU1V1d5LX0zzZfv5sBaxES5JGLKxxnWhJkibHcGbnvqWqdp3pYFW9cMbbJz9JsnVV3dR2tb65y2k3ANt1fN623QcwU/mZytzQvl/vWlV1a8f+TwDvninuKfavkySN1FR37kE2SZI0ROMfE30acED7/gDg813OuRDYKcmOSTammTDstDnKnwbsm2STJDvSTCB2Qdv1+84kz2xn5d5/qkxbCZ/yUuDKuYL3SbQkSZIkaSEdDZyS5LXAN9jlogAAIABJREFUD4FXAyR5FPCJqtq7qlYnOQQ4g2YE93FVdfls5avq8iSnAFcAq4E3VdWatswbgeOBBwJfbjeANyd5aXv+bcCBcwVvJVqSNHJ255YkaYKEsU4s1nahfkGX/TcCe3d8Ph04fb7l22NHAUd12X8R8KQu+w8HDu8hfCvRkqTRqopdsiVJmiTDGRO9bPlXjSRp5NbUBgNt85FkzyRXJVmZ5LAuxzdN8oUklya5PMlBQ/+ikiQtFuMdE72oWYmWJC16SVYAHwX2AnYG9kuy87TT3gRcUVVPBnYH3ttOVCJJkjRvy7wNQZI0agWsHf2Y6N2AlVV1LUCSk4F9aCYW6Qzloe2snA+hmTxk9agDkyRp4ox5TPRiZyVakjRimXeX7AFsA1zX8fl64BnTzvkIzdIXNwIPBf64qtaOOjBJkiaOY6IH4k8nSRqpZp3ogZ9Eb5Hkoo7Px1bVsT1eYw/gEuD5wG8BZyb5elXdOWhwkiQtKlaiB+JPJ0laDG6pql1nOX4DsF3H523bfZ0OAo6uqgJWJvkB8HjggqFGKkmSljQr0ZKkkVsz+nksLwR2SrIjTeV5X+A10875Ec2akl9PsiXwOODaUQcmSdJEsibYN386SdJIFRlGd+7Z71G1OskhwBk0U6UcV1WXJ3lDe/wY4O+A45N8h6Yj29ur6paRBiZJ0iRyYrGBWImWJI3c2gVYUbGqTgdOn7bvmI73NwIvHnkgkiRNOsdED8R1oiVJkiRJmifbHyRJI1UFa0bcnVuSJPXAJ9ED8aeTJI3cqMdES5KkHjkmum9WoiVJI9VMLOboIUmSJoZPogfiXzWSJEmSJM2T7Q+SpJFbg925JUmaGD6JHog/nSRppArHREuSNFGsRA/En06SNGKOiZYkaeI4sVjfFuVfNUmOTHJDkkvabe+OY4cnWZnkqiR7jDNOSZKWC3OzJGm5WMxPot9fVe/p3JFkZ2Bf4InAo4Czkjy2qtaMI0BJUmOtY6KXC3OzJC0GduceyFL76fYBTq6qu4EfJFkJ7AZ8Y7xhSdLyVQVrHBO9nJmbJWnSWIkeyKLszt06NMllSY5L8rB23zbAdR3nXN/ukyRJo2duliQteRNbiU5yVpLvdtn2AT4OPBrYBbgJeG8f1z84yUVJLvrprfYok6RRWlsbDLRpMixkbr5jyLFLkqZZMeC2jE3sQ/yqeuF8zkvyv4Evth9vALbrOLxtu6/b9Y8FjgXY9ckPqP4jlSTNpohLXC0RC5mbH5eYmyVpVOzOPZBF2byfZOuOjy8Hvtu+Pw3YN8kmSXYEdgIuWOj4JEnrWksG2jT5zM2StIhMVaIH2Zaxxfr1351kF6CAVcDrAarq8iSnAFcAq4E3OfunJEkLYuS5eaMu++7tL1ZJkvq2KCvRVfUnsxw7CjhqAcORJM2iwO7cy4C5WZIWkbDsxzUPYlFWoiVJi4uTg0mSNEEcEz0QfzpJ0miVE4tJkjRxrAn2zUcDkiRJkiTNk5VoSdJIFc7OLUnSRBnz7NxJHp7kzCTXtK8Pm+G8PZNclWRlksPmUz7J4e35VyXZo2P/UUmuS3LXtHtskuTTbZlvJtlhrvitREuSRm5t26W7302SJA3R1MRig2yDOQw4u6p2As5uP68bYrIC+CiwF7AzsF+SnWcr3x7fF3gisCfwsfY6AF8AdusSy2uBn1XVY4D3A/84V/BWoiVJIzU1O7eVaEmSJsT414neBzihfX8C8LIu5+wGrKyqa6vqHuDkttxs5fcBTq6qu6vqB8DK9jpU1flVddMcsXwWeEGSWf/4sBItSZIkSVpIW3ZUaH8MbNnlnG2A6zo+X9/um638bGVmcl+ZqloN3AFsPlsB52STJI2cT5MlSZowg9cEt0hyUcfnY6vq2KkPSc4CtupS7ojOD1VVSarfIAYt3w8r0ZKkkSrski1J0kSZGhM9mFuqateZDlbVC2e8ffKTJFtX1U1JtgZu7nLaDcB2HZ+3bfcBzFR+tjIzmSpzfZINgU2BW2crYHduSdLIOTu3JEkTZPxjok8DDmjfHwB8vss5FwI7JdkxycY0E4adNkf504B92xm3dwR2Ai7oIZZXAl+pqlmfbFuJliRJkiQtpKOBFyW5Bnhh+5kkj0pyOtw3PvkQ4AzgSuCUqrp8tvLt8VOAK4B/B95UVWvaa787yfXAg5Jcn+TI9lr/AmyeZCXwVrrMFD6d3bklSaNVjomWJGmiTD2JHpOquhV4QZf9NwJ7d3w+HTh9vuXbY0cBR3XZ/zbgbV32/xp4VQ/hW4mWJI3W1BJXkiRpggw+JnrZshItSRo5K9GSJE2QMT+JXuwcEy1JkiRJ0jzZ/iBJGimXuJIkacL4JHog/nSSpJErK9EakY1m2H/vgkYhSYuMleiB+NNJkkbOtZ4lSZowTizWN8dES5IkSZI0Tz6JliSNVLlOtPoQZu6qLUkakN25B+JPJ0kaOcdES5I0QaxED8Tu3JIkSZIkzZPtD5KkEXOJK0mSJo4Ti/XNSrQkaeTszi1J0gSxO/dA/OkkSSNVOLGYJEkTxUr0QBwTLUmSJEnSPNn+IEkarWqWuZIkSRPCJ9ED8aeTJI3cWuzOLUnSJCknFuublWhJ0kgVTiwmSdIkqcAaa4J9c0y0JEmSJEnzZPuDJGnEXCdakqSJ4pPogfjTSZJGzonFJEmaHBVYvWLQTslrhxLLYmQlWpI0co6JliRpclTCmg0HrQreM5RYFqMZf7kk1w7h+h+oqg8N4TqSJC175mZJksZvtuaHHYA7gNv7vPb2wGZ9lpUkLRFVPokeoh0wN0uShmDNCte46tdcz/DfX1Xv7OfCSUbWST7Jp4HHtR83A26vql2S7ABcCVzVHju/qt4wqjgkSfPjxGJDZW6WJA2kCGuwEt2vRTkmuqr+eOp9kvfStMpP+X5V7bLwUUmSZuLEYkufuVmSFo8irLYS3bfZKtFPB24c4NqDlp9TkgCvBp4/yvtIkjQhzM2SJI3ZjJXoqrp4kAsPWn6eng38pKqu6di3Y5JLaFrA/7qqvr4AcUiSZuGY6OEwN8/fRl323Tvqm0rSIrJmcXZKnggT+8slOQvYqsuhI6rq8+37/YBPdRy7Cdi+qm5N8jTg35I8saru7HL9g4GDAbbfZmJ/Bkla9IpYiV4iFjI3bznc0CVJHRwTPZi+ao9tV62t6N7QS1X9aJCg2mu8cI4YNgT+CHhaR5m7gbvb9xcn+T7wWOCiLtc/FjgWYNcnP8DRepI0Qv5HdvSWWm5+fOI/G0kaESvRg+mpEp3kVcBhwJNmKVu9XrdPLwS+V1XXT+1I8gjgtqpak+TRwE7AMNbUlCRpIpmbJUlaWPNOqEneBHwIWA38J3BD+35c9mXd7mIAzwHemeReYC3whqq6bcEjkyTdz3WiR8bcLEnq1zifRCd5OPBpYAdgFfDqqvpZl/P2BD4IrAA+UVVHz1U+yeHAa4E1wJur6ox2/1HA/sDDquohHfc4EPgnmhwK8JGq+sRs8ffSKv0XwM3A71XVD3ooNxJVdWCXfacCpy58NJKkWS1Ax9yZEu20c3YHPkDT5fmWqnru6CMbKXOzJKlnE7DE1WHA2VV1dJLD2s9v7zwhyQrgo8CLgOuBC5OcVlVXzFQ+yc40DbpPBB4FnJXksVW1BvgC8BGgc+LLKZ+uqkPmG/wGPXzRbYDPTEKSliQtLlUZaJtLR6LdC9gZ2K9NpJ3nbAZ8DHhpVT0ReNXwv+mCMzdLknrWjInecKBtQPsAJ7TvTwBe1uWc3YCVVXVtVd0DnNyWm638PsDJVXV3mxtXttehqs6vqpsGDRx6q0RfB2wyjJtKkjRksyXaKa8B/t/UBFtVdfMCxzgK5mZJ0mK0ZUeF9sd0X5RhG5o8N+X6dt9s5WcrM5tXJPlOks8m2W6uk3upRJ8A7JXkoT2UkSSJqsG2eZhP0nws8LAk5yS5OMn+w/l2Y2VuliT1ZQ0rBtqALZJc1LEd3Hn9JGcl+W6XbZ1G7qoqBhj4NWh5mm7eO1TVbwNncv8T7hn18hz+aOCpNP3K3wZ8q6p+3leYkqRloxjKxGJbJOlcEunYdjmkXmxIs/TSC4AHAt9Icn5VXT1ocGNkbpYk9WxIS1zdUlW7zniPWZZFTPKTJFtX1U1JtqaZ32O6G4DOp8Lbcv/kXzOVn63MTHHe2vHxE8C7ZzsfengS3Q7G/ijwGOArwO1J1nTZxjkrqCRpabqlqnbt2KZXoOeTNK8HzqiqX1TVLcC5wJNHF/LomZslSYvUacAB7fsDgM93OedCYKckOybZmGbCsNPmKH8asG+STZLsSLOs4gWzBdJWwqe8FLhyruB7WeJqH+CzNLOe/gC4kfEuoyFJWgwKGP0SV/clWprK8740Y6A7fR74SJINgY2BZwDvH3Vgo2RuliT1o2Dcs3MfDZyS5LXAD4FXAyR5FM0KG3tX1eokhwBn0OS546rq8tnKV9XlSU4BrqDJh29qG5xJ8m6avw0elOT69j5HAm9O8tL2/NuAA+cKvpfu3EcCvwT+oKrO66GcJGmZm+e45gGu3z3RJnlDe/yYqroyyb8Dl9GsV/yJqvruaCMbuSMxN0uSepZhzLDdt7YL9Qu67L8R2Lvj8+nA6fMt3x47Cjiqy/63AW/rsv9w4PAewu/pl3sccKJJWpLUswVYJ7pboq2qY6Z9/ifgn0YfzYIxN0uSejakMdHLVi+zc98C3DOqQCRJUs/MzZIkLbBenkSfCuyZZKOqundUAUmSlpoMY3ZudWduliT1xSfR/evlSfRfAz8DPpNkh5FEI0lammrATTMxN0uSejbVnXvAdaKXrV6eRH8H2IhmNtOXJLkduKPLeVVVvzWM4CRJS0ANZZ1odWduliT1rMi4Z+de1HqpRG9AM+33jzr2dfuryL+UJElaGOZmSZIW2Lwr0VW1wwjjkCQtZXbJHglzsySpX+Nc4mqx85eTJC0AH4RKkjQpXOJqMFaiJUmj55NoSZImhpXowcw4O3eSVyfZud8LD1pekiSty9wsSdL4zbbE1cnAKwe49qDlJUlLhUtcDYu5WZI0FKtZMdC2nM3Vnds/XSRJgynAJa6GydwsSRpI053bkb39muuXOzLJkX1e2yQvSQKgzAjDZG6WJA3EMdGDma0SfS6DJ9tVA5aXJEn3MzdLkjRmM1aiq2r3BYxDkrSU+fxzKMzNkqRh8Ul0/+wIL0kaPcdES5I0MYos+8nBBmElWpI0cvFJtCRJE8OJxQbjLydJkpaFjWbYf++CRiFJWuysREuSRsu1niVJmjiOie6flWhJ0ojFMdGSJE0Ql7gazAbjDkCSJEmSpMXCJ9GSpNGzO7ckSRPDJ9GDsRItSRo9K9GSJE0Ul7jqX8+V6CSPBHYFHgbdf/mqOnHAuCRJS4mV6JEyN0uSeuESV4OZ9y+XZCPgGGB/Zh5LHZo/lUzUkiSNmLlZkqSF10vzw98BBwHfB04CrgNWjyIoSdISUjg79+iYmyVJPXNM9GB6qUS/BrgaeEpV/WpE8UiSlqDYnXtUzM2SpL5Yie5fL0tcPRI4faGSdJJXJbk8ydoku047dniSlUmuSrJHx/6nJflOe+xDSXz0IUmToAbcNBNzsySpZ0VYzYqBtuWsl0r0j4DfGFUgXXwX+CPg3M6dSXYG9gWeCOwJfCzJ1P+LHwdeB+zUbnsuWLSSJC08c7MkSQusl0r08cBeSTYdUSzrqKorq+qqLof2AU6uqrur6gfASmC3JFsDv1FV51fV1AQqL1uIWCVJGpPjMTdLkno0NTv3INty1ksl+mjgPOCsJM9LspAt3522oZk4Zcr17b5t2vfT90uSxiw12KYZLdncfBVb8yyOXGeTJA3PGlYMtA0iycOTnJnkmvb1YTOct2c7TGhlksPmU77b8KIkD0rypSTfa4clHd1x/iZJPt2W+WaSHeaKf8YmhCRr6T4SLcBZ7TndilZVzatpIslZwFZdDh1RVZ+fzzX6leRg4GCA7bdZ3i0pkjRyzs49FMspN8OCPFyXpGVpAmbnPgw4u6qObivHhwFv7zyhHRb0UeBFNI2wFyY5raqumKn8tOFFj6JpZH5se8n3VNVXk2wMnJ1kr6r6MvBa4GdV9Zgk+wL/CPzxbMHPllDPZcTTuVTVC/sodgOwXcfnbdt9N7Tvp++f6d7HAscCbP+k36gP3779OscP3exHfYQmSdJILZvcnOxU8NQ+QpEkLQL7ALu3708AzmFaJRrYDVhZVdcCJDm5LXfFLOXvG14E/CDJSmC3qvoG8FWAqronybe4Pz/tA/d1d/os8JEkaYchdTVjJbqqdp/p2JidBnwyyftoWhd2Ai6oqjVJ7kzyTOCbwP7Ah8cYpyQJnGF7iMzNkqRhGfOT6C2r6qb2/Y+BLbuc022o0DPmKL8NcP60MusMI0qyGfAS4IPT71NVq5PcAWwO3DJT8BPbjznJy2kS7SOALyW5pKr2qKrLk5xC0wKxGnhTVa1pi72RZpKVBwJfbrc53XrPgzlp1W7r7Dt0F59ES9LQWIleEhYyN8MDgMcPNX5JUmNqiasBbZHkoo7Px7Y9ioDZhwetE0tVJf3PgNJL+SQbAp8CPjT1hLsf865EJ7kW+EBVfWiWc94E/GVVPbrfgKZU1eeAz81w7CjgqC77LwKeNOi9JUnD5eRgo2FuliT1Y2p27gHdUlW7zniPWYYHJflJkq2r6qZ2JYebu5w201AhgJnKz1YGmiFD11TVB7rc5/q2kr0pcOtMsUNvs3PvAGw2xzmbAb/ZwzUnwupfb8gtV2++ziZJ0iKwA0s0N8Mm3L+09NQmSVoiTgMOaN8fAHSbuPJCYKckO7aTge3blput/GnAvu2M2zvSDi8CSPL3NBXkP58lllcCX5ltPDQMvzv3Q4F7hnzNkVvxa9jsyl7aEyRJPfFJ9DgtytzMBsCDp+37+TgCkaSlacxjoo8GTknyWuCHwKsBkjwK+ERV7d2OTz4EOANYARxXVZfPVn6m4UVJtqXpRv494FvtShYfqapPAP8C/Gs7CdltNJX1Wc1aiU6y/bRdm3XZR/ultgdeAfTdt1yStERZiR4ac7MkaVDjXuKqqm4FXtBl/43A3h2fTwdOn2/59th6w4uq6nqa5SC7nf9r4FU9hD/nk+hVrPunz1vabSYB3tpLAJNgw1+uZfNLfzHuMCRpSUo5JnrIVrEMcjMbs+7iWABXjiMQSVp6hjSx2LI1VyX6RJpEHZplKS4DLuly3hqawddnV9V/DDVCSZLUydwsSdIYzVqJrqoDp94n2R/4XFW9c9RBLbhf/ArOv3TcUUjS0lVde1CpD8slN2/964t5/ZXr/rt51gPX79Jw3q/8tyVJ/RjC7NzL1rx/uapaVjNvffj29YeXHbqZa0dLUl/szj0Syy03s02XfSsXPApJWvTGPSZ6seur+SHJdsBTaKYIvwP4dlVdN8zAJElLh2OiR8/cLEmaLyvRg+mpEp1kJ+BjwPO7HPsKzRTiVw8ptrE6adVu6+07dBefREuSJstyys1s3WWfT6IlSQts3pXoJI8B/gvYHPg+cB7wY2Ar4Fk0U4yfl+T3qsqUJkm6n0+iR8LcLEnql0+i+9fLk+h/oEnSbwE+WlVrpw4k2QA4FHg/8C7axa4Xs1uu3nz9nbssfByStOi5xNUoLavcvN6SV5KkvrjE1WB6qUS/ADi9qj48/UCbtD+YZA/ghcMKTpIkzcrcLEnSAuulEr0x3deh7PRt4Nn9hzM5NrtyeU14Kkkj5ZPoUVlWubnr7NySpJ41E4u5xFW/evnlLgUeM8c5jwEu6z+cybH5pb8YdwiStHRYiR6VZZWb2W7cAUjS0uGY6P71Uol+F/C5JHtV1ZenH0zyB8DLgZcNKzhJ0tLgmOiRMTdLknrmEleD6aUSvTnwZeCLSc4GzgV+AmwJPJdmaY0vAFsk2b+zYFWdOJxwF9D5l447AkmS5rKscvORb8l6+551aPcWmvM+vP65kiQNQy+V6ONpOuSFZoKSbpOUvBR4ScfntGUWXaKWJGkROB5zsySpR87OPZheKtEHjSyKReLDt2/fdf+hm/1ogSORpEXG7tyjsuxzM1uNOwBJWpycWKx/8/7lquqEUQYiSVqiXCd6ZMzNkqR+OCZ6MDY/9OCkVbt13X/oLj6JliRpLLYddwCSpOWm50p0kkcArwCeADy4qv6sY/+OwHeq6ldDjXJC3HL15t0P7LKwcUjSouOT6JFazrl5GN25N+qy797BLytJE8sn0YPpqRKd5LXAh4AHcP/EJH/WHt4S+AZwMPAvQ4xRkrTYWYkeGXOzJKkfTizWv3lXopO8CDgWuAz4W2AP4A1Tx6vqu0kup1mLckkm6s2u3GDcIUjSohMcEz0q5mbszi1JfWieRDuyt1+9/HJvB24CnltVdyZ5SpdzLgN+dyiRSZKkuZibJUlaYL1UoncFTq6qO2c553qW8GITm1/6i3GHIEmLk0+iR2XZ5+aNt5rtq0uSunFM9GB6qURvDMxVi9wMWNN/OBPu/EvHHYEkLT4ucTVKyz43/6/NN+26/1nnr/+P7rxnZtThSNKiYSW6f70M8l0FPG2Oc54BXNV3NJKkpakG3OYhyZ5JrkqyMslhs5z39CSrk7yy7+8zOVZhbpYk9WjqSfQg23LWy5PozwNvS/KqqvrM9INJDgJ+BzhiWMEtFh++ffv19h26mWtHS9JCSbIC+CjwIpruyxcmOa2qruhy3j8C/7HwUY6EuXkm29097ggkSUtUL5XodwP7Ap9qW+83BUhyCPBs4I+Aa4APDztISdIiN/ru3LsBK6vqWoAkJwP7AFdMO+9Q4FTg6SOPaGGYmyVJPStc4moQ865EV9XPkjwXOBF4VcehD7WvXwdeU1XLbvatk1bttt6+Q3fxSbQkTRnCmOgtklzU8fnYqjq24/M2wHUdn6+n6cZ8fwzJNsDLgeexRCrR5uaZbb7VLeMOQZImmEtcDaKnX66qfgTsnuR3aJbL2By4Azi/qi4eQXySpKVg8Er0LVW164DX+ADw9qpamyydCabMzZKkXjk792D6an6oqsto1p0UcMvVm6+/c5eFj0OSlrEbgO06Pm/b7uu0K3ByW4HeAtg7yeqq+reFCXG0zM3r2nyDW8cdgiRpieq5Ep3kN4FH0DxX+GnbAr6sbXZlL5OcS9Iy08MM2wO4ENgpyY40led9gdesE0bVjlPvkxwPfHGpVKDNzevbHCvRkjQbn0T3b16V6CRbAP8L2A945LRjPwFOAv6hqm4beoSSpEVv1OtEV9XqdjKtM4AVwHFVdXmSN7THjxltBAvP3CxJ6lcRJxYbwJyPUJPsBFwEvAXYElgD3Az8tH2/FfBW4KIkjx5WYEleleTyJGuT7Nqx/0VJLk7ynfb1+R3HzmnXCL2k3R7Z/erDtfmlv1hvkyQtrKo6vaoeW1W/VVVHtfuO6VaBrqoDq+qzCx/lcJib57YFt6y3SZImQ5KHJzkzyTXt68NmOG/PNoesTHLYfMonObw9/6oke7T7HpTkS0m+1+axozvOPzDJTzvy1J/NFf+slegkG9C0ZG8PfA14IfCQqtq6qrYCHgq8GDgX2AH4v3PdsAffpVma49xp+28BXlJVvw0cAPzrtOP/rap2abebhxjP/2/v3sPtqus7j78/BgUvlUhQ5CKPtM1U0bapYqyj470KaMVLdaDWW6nUKqmtz7TCw7RiK1PaTq94TdEROmhksJS0ogygtKVjBLSkFaIYMWoigpGLgoqEfOePtSI7J/tccvbZe+1z9vv1POs5e/3W+q39PdvI93z377d+S5I0XzXgph8xN0uSBlXt6tyDbAM6Bbi8qlYCl7f7u0myDHgXcAxwJHBCkiNn6t8ePx54HHA08O72OgD/s6oeA/wc8NQkx/S83Ud68tTZswU/22//PJqFWM4HTqiq3f6Uqaq7gcuSXA58BHhZkl+oqktne+PZVNUmgKkrqFbVv/XsXgc8MMm+bSzd2LCxs7eWpMVg2NO5J4y5eQ6emFfs0XZyfbrvue/MU4YdjiSNnY7viT4OeGb7+hzgCuCtU85ZDWyuqhsBkqxr+10/Q//jgHVt/vlKks3A6qr6NPApgKr6YZLP0SxCOi+zTed+GXA3sGZqku7VHjsZuAf4pfkGMw8vAz43JUmf0w7D/16W0jNMJGkxcyR6IZmbJUkD2fWIq0E24MAk1/RsJ+1FCAdV1U3t62/S3Jo01aHA13v2t7ZtM/WfqQ8ASZYDv0gzgr3Ly9pbki5I0vu0j75mG4l+AvCvVfWt2S5UVbckubLtMydJLqO5b2uq06rqoln6Pg74Y5pv5Hd5ZVVtS/JjwEeBVwHnTtP/JOAkgP140FxDnrOzbj+8b/ua5RO/YKokaTATkZv3n2vAe8EVuyVpQW2vqqOmOzhTPundqapK5j9nbW/6J9kH+DDw17tGuIF/AD5cVXcn+XWake1nT3cNmL2IfhRw5VwCal1Hs0ronFTVc/fi2j+S5DDgQuDVVfXlnutta39+N8mHaKYA9E3UVbUWWAvw0Byw4OMc521Z3bd9zSqLaEkTxtHkhTYRufmQAf6gms4KFxeTJKAdid453OncM+WTJDcnObiqbkpyMM3imFNto8l5uxzWtgFM13+mPtDkmC9V1V/2xNn7DevZwJ/M8qvNOp37ocDts12kx+00C5oMTTv8/jHglKr61572fdrHfZDk/sALaRZAkSR1KAuwaTfmZknSYAp27Fg20Dag9TQLUdL+7DfT6WpgZZIjkjyAZsGw9bP0Xw8cn2TfJEcAK4GrAJK8g2ai02/1vklbhO/yImDTbMHPNhL9AJpHZczVzrbPwJK8BDgLeDjwsSTXVtXzae7v+kng95P8fnv684C7gEvaJL0MuAz4m4WIZT6237Ci/4FVo41DksaCI9ELydw8Twc6nVuSAKgK9+4YeIXtQZwJnJ/kROCrwCsAkhwCnF1Vx1bVjiQEavzyAAAf3klEQVQnA5fQ5JAPVNV1M/WvquuSnE+z+NgO4E1VdW87W+o04AvA59rlOd7ZrsT9m0le1J5/K/Da2YKfyyfXyZ8+VXUhzbSwqe3vAN4xTbcnDjUoSZLGg7lZkrRotVOon9On/RvAsT37FwMXz7V/e+wM4IwpbVuZZnJbVZ0KnLoX4c+piD49yel7c1HB8k2zzZSXpMnhI64WnLl5HlxYTJIazUh0p4+4WtTmUkTv7e1o/qnUWrHxrq5DkKTxYGZYaObmeXJxMUkCCovoAcxYRFeVw6nzZAEtST0s4RaMuXn+LKAlqVEVdtxjET1fnd5NvqRt2Nh1BJIkqcftWdu3/YL6oz3ajste3R4nSZogFtGSpOEq74mWJGm8hJ33WgrOl5/ciJ11++F7tK1Z/rUOIpGkEbKI1hhzmrekiVOA90TPm0W0JGnoHImWJGmMVCyiB2ARPWLnbVm9R9uaVY5ES5LUlQN99JUkaS9YRI/Y9htW7Nm4avRxSNJIORKtMebzoyVNnAJ27O3TErWLRbQkaeiczi1J0pjZ0XUAi5dF9Igt3+TjPSVNmMKRaI215dzWdQiSNFqFRfQArOgkSZIkSZojR6JHbMXGu7oOQZJGz5FojbF+C4vd0UEckjQyjkQPxCJakjRUwXuiJUkaKwXc03UQi5dF9Kht2Nh1BJI0ehbRGmPrs2mPtp+uY/qee0A+PuxwJGn4Cri36yAWL4voMXDW7Yf3bV+z3OdHS5LUhemeHb1zxHFIksaPRbQkaehSDkVLkjRWvCd63iyix8B5W1b3bV+zypFoSUuAj7jSIrRimpHob404DkkaChcWG4iPuJIkSZIkaY4ciR4D229Y0f/AqtHGIUnD4urcWmyWc3vfdkeiJS0JjkQPxCJakjR8FtGSJI0Pi+iBWESPgeWbnFUvaWlzJFqLzYHf639P9JdGHIckDYVF9EAsosfAio13dR2CJEnqsW//GlqSJItoSdIIOBItSdJ4cSR63iyix8GGjV1HIEnDU07n1uLz9sP7t7+pHrJH2+dz55CjkaQFVsA9XQexeFlES5KGzyJakqTxUcC9XQexeFlEj7Gzbt/za/A1y7/WQSSSJAlg+XccdZakSWcRPcbO27J6j7Y1qyyiJS0uwencWjr2ccExSUuBq3MPxCJakjR8ZRUtSdLYsIgeiEX0GNt+w4o9G1eNPg5JGpQj0Voybuk6AElaABbRA7lf1wFIkiRJkiZHkgOSXJrkS+3Ph01z3tFJvphkc5JT5tI/yant+V9M8vye9k8k2ZjkuiTvTbKsbd83yUfaPp9J8ujZ4reIHmPLN91vj02SFp1agE0aF7f22SRpMdox4DaYU4DLq2olcHm7v5u2yH0XcAxwJHBCkiNn6t8ePx54HHA08O5dxTLwiqr6WeDxwMOBl7ftJwK3VdVPAn8B/PFswVuVSZKGLjsH2yRJ0gLaNZ27uyL6OOCc9vU5wIv7nLMa2FxVN1bVD4F1bb+Z+h8HrKuqu6vqK8Dm9jpU1Xfac/YBHsB9X9P3XusC4DlJMlPw3hM9xlZsvKvrECRpYTiarKXC1bklLQXd3xN9UFXd1L7+JnBQn3MOBb7es78VePIs/Q8FNkzpc+iunSSX0BTVH6cpmHd7n6rakeQOYAWwfbrgx7aITvJy4HTgscDqqrqmbX80sAn4Ynvqhqp6Q3vsicAHgQcCFwNvrlrES8Ju2Nh1BJIk/Yi5Gd7+qj3b3nZl/3OveNpwY5Gkjh2Y5Jqe/bVVtXbXTpLLgEf26Xda705VVTL/JUj3pn9VPT/JfsB5wLOBS+fznmNbRAOfB14KvK/PsS9XVb91qt8DvB74DE2iPprmWwZJUodcnXvJMDdL0lJQwD0DX2V7VR017VtUPXe6Y0luTnJwVd2U5GD6P/tgG/Conv3D2jaA6frP1GdXXD9IchHNNO5Le/psTbIPsD+zzDsa2yK6qjYBzDId/UfaD++hVbWh3T+XZm78kkrUZ91+eN/2Ncu/NuJIJGmOCp8TvUSYm6fh4mKSFpsC7u00gvXAa4Az258X9TnnamBlkiNoCt3jgV+epf964ENJ/hw4BFgJXJXkIcCPtUX3PsALgH+Zcq1PA78EfHK2GVNjW0TP4ogk1wJ3AP+9qv6FZi771p5zdpv/LknqjiPRE8HcLEmLSbf3RJ8JnJ/kROCrwCsAkhwCnF1Vx7b3J58MXAIsAz5QVdfN1L+qrktyPnA9zW/4pqq6N8mDgfVJ9qVZXPtTwHvba70f+Nskm2m+Fj1+tuA7LaJnmidfVf2+jQC4CTi8qr7d3mf190keN4/3Pgk4CWA/HrS33Ttz3pbVfdvXrHIkWpI0uHHJzfvvbecuudiYJO2Vqvo28Jw+7d8Aju3Zv5jmVqA59W+PnQGcMaXtZuBJ05z/A+573NWcdFpEzzRPfoY+dwN3t68/m+TLwH+iGeI/rOfUPea/T7nOWmAtwENzgGMkkjRM/ld20RiX3HzIAIvMSJJm0f3q3IvaopvOneThwK3tsPyP08xzv7Gqbk3ynSQ/T7N4yauBs7qMdRi237Ci/4F+S7lI0hgITude6iY9N3tPtKRFxyJ6IGNbRCd5CU2ifTjwsSTXVtXzgacDf5DkHmAn8Iaq2pW+3sh9j9H4OEtt4RJg+ab7dR2CJGlCmZun4XRuSZooY1tEV9WFwIV92j8KfHSaPtcAjx9yaJKkvVHl6txLhLlZkpaIhXnE1cQa2yJa/a3YeFfXIUjSXnM6t5a0fk83laRx1v0jrhY1i2hJ0vBZREuSNF68J3reLKIXmw0bu45AkiT1ePvZ/dvf9sd7tl3x1uHGIkkaPovoJeKs2w/fo23Ncp8dLWk8OJ1bE8lVuyWNK1fnHohFtCRpuArYaRUtSdLYcGGxgVhELxHnbVm9R9uaVY5ESxoT1tCaRD76StK4cmGxgfjQYUmSJEmS5siR6CVi+w0r9mxcNfo4JKkf74nWRPKeaEnjynuiB2IRLUkavrKKliRprFhEz5tF9BKxfJMz8yWNr1GMRCc5GvgrYBlwdlWdOeX4K4G3AgG+C/xGVfncQA2P90RLGlcuLDYQi+glYsXGu7oOQZI6k2QZ8C7gF4CtwNVJ1lfV9T2nfQV4RlXdluQYYC3w5NFHq4nhdG5JWpIsoiVJw1WMYnXu1cDmqroRIMk64DjgR0V0Vf2/nvM3AIcNPSpJksaRq3MPxCJ6qdjgjERJ4ylAhn9P9KHA13v2tzLzKPOJwMeHGpEm3tv/Y8+2t720/7lX/N1wY5Gk3biw2EAsoiVJw7dz4CscmOSanv21VbV2PhdK8iyaIvppA0clSdJiZBE9EIvoJeys2w/v275m+ddGHIkkDWx7VR01w/FtwKN69g9r23aT5GeAs4FjqsplnzR6/quTpEXPIlqSNHQjmM59NbAyyRE0xfPxwC/vFkNyOPB3wKuq6oZhByRJ0thyde6BWEQvYedtWd23fc0qR6IljdAIFharqh1JTgYuoXnE1Qeq6rokb2iPvxf4fWAF8O4kADtmGd2WFp4rdksaFy4sNm8W0UvY9htW9D+warRxSJp0BcMfiaaqLgYuntL23p7Xvwb82tADkWbidG5J48B7ogdyv64DkCRJkiRpsXAkeglbvsnvSCSNhwx/IFpaFL7jdG5J48CR6IFYREuShm8E07klSdIcubDYQCyil7AVG+/qOgRJgoIM/pxoaUm4+QddRyBJNEW0C4vNm0X0UrZhY9cRSJKkHh+apv1th/Rvv+IbQwtFkjRPFtGSpOFzOre01yygJQ2V90TPm0X0BDrr9sP3aFuz3GdHSxoia2hpRi44JmmkXFhsIC7fLEmSJEnSHDkSPYHO27J6j7Y1qxyJljQ8cTq3NCMXHJM0Uh2vzp3kAOAjwKOBLcArquq2PucdDfwVsAw4u6rOnK1/klOBE2mWTvvNqrqkbf8EcDBNDfwvwJuq6t4krwX+FNjWvu07q+rsmeJ3JFqSNHxVg22SJGnh7Fqde5BtMKcAl1fVSuDydn83SZYB7wKOAY4ETkhy5Ez92+PHA48Djgbe3V4HmkL7Z4HHAw8HXt7zdh+pqlXtNmMBDY5ET6TtN6zYs3HV6OOQNCEK8BFX0oxu6ToASZOl+3uijwOe2b4+B7gCeOuUc1YDm6vqRoAk69p+18/Q/zhgXVXdDXwlyeb2Op+uqu+05+8DPIABVmyxiJ5Ayzc5AUGSpHFyc5+2A0YehSSNzEFVdVP7+pvAQX3OORT4es/+VuDJs/Q/FNgwpc+hu3aSXEJTVH8cuKDnvJcleQbwReC3q6r3ffdgES1JGqpQ3hMtSdK4GXwk+sAk1/Tsr62qtbt2klwGPLJPv9N6d6qqksz7D4W96V9Vz0+yH3Ae8GzgUuAfgA9X1d1Jfp1mZPvZM13HInoCrdh4V9chSJo0FtHSjPo94cqRaElDszALi22vqqOmfYuq5053LMnNSQ6uqpuSHEz/u1q2AY/q2T+M+xb/mq7/TH12xfWDJBfRTP2+tKq+3XP4bOBPpot7F+f1SpKGz4XFJEkaH90vLLYeeE37+jXARX3OuRpYmeSIJA+gWTBs/Sz91wPHJ9k3yRHASuCqJA9pi22S7AO8APhCu39wz3u+CNg0W/CORE+iDRu7jkCSJPXY1qftsdOc2+FTaSRpoZwJnJ/kROCrwCsAkhxC8yirY6tqR5KTgUtoHnH1gaq6bqb+VXVdkvNpFh/bwX2PsXowsD7JvjQDyZ8C3tte6zeTvKg9/1bgtbMFbxEtSRouV+eWJGm8dLw6dzuF+jl92r8BHNuzfzFw8Vz7t8fOAM6Y0nYz8KRpzj8VOHUvwh/fIjrJy4HTab6IXV1V17TtrwR+p+fUnwGeUFXXJrmC5gHa32+PPa+qfGrEHDz/kJ/tOgRJY+FLQ7mqC4stDebm0Zru/40PHWkUkpak7h9xtaiNbRENfB54KfC+3saqOo9mNTWS/DTw91V1bc8pr9yV1CVJY8IieqkwN4/QV7sOQNLStTALi02ssS2iq2oTQJKZTjsBWDeSgCRJmnDmZkmSxriInqP/SrM0ea9zktwDfBR4R5XDH5LULVfYnjDmZklaDAZfYXtidVpEz/QA7qrqt8x5b98nA9+rqs/3NL+yqrYl+TGaRP0q4Nxp+p8EnASwHw+aT/iSpLkoLKIXkXHJzfvPJ3hJ0tyZmuet0yJ6pgdwz8HxwIenXG9b+/O7ST4ErGaaRF1Va4G1AA/NAf4TkqRhcnXuRWNccvMhiblZkjSW7td1APOR5H40zwJb19O2T5ID29f3B15IswCKJEkaMnOzJGlSjG0RneQlSbYCTwE+luSSnsNPB75eVTf2tO0LXJLk34FrgW3A34wsYEnStFI10KbxYG6WJGmMFxarqguBC6c5dgXw81Pa7gKeOPzIJEl7zUJ4STA3S5I0xkW0JGmJKGCnRbQkSePDB0UPYmync0uSJEmSNG4ciZYkDZnPiZYkabwUsKPrIBYti2hJ0vBZREuSNEaczj0Ip3NLkiRJkjRHjkRLkobPkWhJksaI07kHYREtSRouV+eWJGnMOJ17EBbRkqQhK6idXQchSZJ+xCJ6EN4TLUmSJEnSHDkSLUkaPu+JliRpzHhP9HxZREuShst7oiVJGjNO5x6ERbQkafgciZYkaYy4OvcgvCdakiRJkqQ5ciRakjR8jkRLkjRGnM49CItoSdKQlUW0JEljxencg7CIliQNVwE7fU60JEnjw5HoQXhPtCRJkiRJc+RItCRp+JzOLUnSGHE69yAsoiVJw2cRLUnSGHE69yAsoiVJQ1aw0yJakqTx4Uj0ILwnWpIkSZKkObKIliQNV0HVzoE2SZK0kHZN5x5km78kByS5NMmX2p8Pm+a8o5N8McnmJKfMpX+SU9vzv5jk+X2uuT7J53v2903ykbbPZ5I8erb4LaIlScO3swbbJEnSAtsx4DaQU4DLq2olcHm7v5sky4B3AccARwInJDlypv7t8eOBxwFHA+9ur7Prmi8F7pzyVicCt1XVTwJ/AfzxbMFbREuShq9qsE2SJC2gbkeigeOAc9rX5wAv7nPOamBzVd1YVT8E1rX9Zup/HLCuqu6uqq8Am9vrkOQhwFuAd8wQywXAc5JkpuAtoiVJkiRJo3RQVd3Uvv4mcFCfcw4Fvt6zv7Vtm6n/TH3+EPgz4HvTvU9V7QDuAFbMFLyrc0uShqsKdnpfsyRJ42NBHnF1YJJrevbXVtXaXTtJLgMe2affabtFUlVJ5j3tbC79k6wCfqKqfnsu9zzPxiJakjR8TsmWJGmMLMgjrrZX1VHTvkPVc6c7luTmJAdX1U1JDgZu6XPaNuBRPfuHtW0A0/Wfrs9TgKOSbKGpgR+R5IqqemZPn61J9gH2B749XezgdG5J0gjUzp0DbZIkaSF1fk/0euA17evXABf1OedqYGWSI5I8gGbBsPWz9F8PHN+uuH0EsBK4qqreU1WHVNWjgacBN7QF9NRr/RLwyaqZv/13JFqSJEmSNEpnAucnORH4KvAKgCSHAGdX1bFVtSPJycAlwDLgA1V13Uz9q+q6JOcD19MMtb+pqu6dJZb3A3+bZDNwK02xPiOLaEnSkLnCtiRJ42VBpnPP/92rvg08p0/7N4Bje/YvBi6ea//22BnAGTO89xbg8T37PwBePvfoLaIlScNW+KxnSZLGyoIsLDaxvCdakiRJkqQ5ciRakjR85eJgkiSNj26ncy92FtGSpKEqoJzOLUnSGHE69yDGdjp3kj9N8oUk/57kwiTLe46dmmRzki8meX5P+xOT/Ed77K+TpJvoJUk/UtWMRA+yzUGSo9u8sDnJKX2Op80Nm9vc8oQF/12XOHOzJC0Vu0aiB9km19gW0cClwOOr6meAG4BTAZIcSbPs+OOAo4F3J1nW9nkP8Hqa54GtbI9Lkpa4Ng+8CzgGOBI4oc0XvY7hvvxwEk3O0N4xN0uSJt7YFtFV9X+ratdXHBuAw9rXxwHrquruqvoKsBlYneRg4KFVtaF9OPa5wItHHrgkaQ+1swba5mA1sLmqbqyqHwLraPJFr+OAc6uxAVje5g7NkblZkpaKXdO5B9km12K5J/pXgY+0rw+lSdy7bG3b7mlfT22XJHVt+AuLHQp8vWd/K/DkOZxzKHDTcENbsszNkrRoubDYIDotopNcBjyyz6HTquqi9pzTaP4XPm+B3/skmul8AHdeVhd8cSGvP+YOBLZ3HcQY8HNo+Dk0/BwaP7XQF/wut11yWV1w4ICX2S/JNT37a6tq7YDXVB/jlJtPB3Pz5PFzaPg5NPwcGguem+GmS+D0QXPzxP5v02kRXVXPnel4ktcCLwSe004DA9gGPKrntMPatm3cN62st326914LTOQfYEmuqaqjuo6ja34ODT+Hhp9DY0qhuiCqahT3wE6XG/b2nIlnbu6G/w1q+Dk0/Bwafg6NRZybl6yxvSc6ydHA7wIvqqrv9RxaDxyfZN8kR9AsUnJVVd0EfCfJz7crf74auGjkgUuSunA1sDLJEUkeQLPI1fop56wHXt2u0v3zwB1t7tAcmZslSRrve6LfCewLXNo+DWNDVb2hqq5Lcj5wPc1UsjdV1b1tnzcCHwQeCHy83SRJS1xV7UhyMnAJsAz4QJsv3tAefy9wMXAszaJX3wNe11W8i5i5WZI08XLfTCxNiiQneS+hn8Mufg4NP4eGn4PUDf+/1/BzaPg5NPwcGn4O48ciWpIkSZKkORrbe6IlSZIkSRo3FtETJMkHktyS5PNdx9KlJI9K8qkk1ye5Lsmbu46pC0n2S3JVko3t5/D2rmPqUpJlSf4tyT92HUtXkmxJ8h9Jrh3GSqCS9mRubpibG+bm3Zmbzc3jyuncEyTJ04E7gXOr6vFdx9OVJAcDB1fV55L8GPBZ4MVVdX3HoY1Uu1Lug6vqziT3B64E3lxVGzoOrRNJ3gIcBTy0ql7YdTxdSLIFOKqqJva5j9KomZsb5uaGuXl35mZz87hyJHqCVNU/A7d2HUfXquqmqvpc+/q7wCbg0G6jGr1q3Nnu3r/dJvJbtSSHAS8Azu46FkmTxdzcMDc3zM33MTdrnFlEa6IleTTwc8Bnuo2kG+00qWuBW4BLq2oiPwfgL2mefbuz60A6VsBlST6b5KSug5E0mczN5uaWublhbh5DFtGaWEkeAnwU+K2q+k7X8XShqu6tqlXAYcDqJBM3lTDJC4FbquqzXccyBp7W/ns4BnhTO81UkkbG3GxuBnPzFObmMWQRrYnU3mf0UeC8qvq7ruPpWlXdDnwKOLrrWDrwVOBF7T1H64BnJ/nf3YbUjara1v68BbgQWN1tRJImibl5d+ZmczOYm8eVRbQmTrtox/uBTVX1513H05UkD0+yvH39QOAXgC90G9XoVdWpVXVYVT0aOB74ZFX9SsdhjVySB7eL+ZDkwcDzgIleLVjS6JibG+bmhrm5YW4eXxbREyTJh4FPAz+VZGuSE7uOqSNPBV5F863mte12bNdBdeBg4FNJ/h24mua+q4l9hIQ4CLgyyUbgKuBjVfWJjmOSljxz84+YmxvmZvUyN48pH3ElSZIkSdIcORItSZIkSdIcWURLkiRJkjRHFtGSJEmSJM2RRbQkSZIkSXNkES1JkiRJ0hxZRGvJSXJFkkWz7HySmrI9ci/6bkmyZYjhDSTJP0753V7bdUySpNEzN48Pc7M0uH26DkCazjyS7euq6oPDiGUEvgp8sH19Z4dxLLQPAdcAq4DjOo5FkjQgc/OSYG6WBmQRrXH29j5tvwXsD/wVcPuUY9e2P18NPGiIcQ3Dlqo6vesgFlpVfQig/ZbbRC1Ji5+5eZEzN0uDs4jW2OqXuNr/4O8P/GVVbZmm39eGGpgkSRPK3CxJ3hOtJajffVdJntne93N6kqOSfCLJHUluS/LRJI9qz/vxJOuSfCvJ95N8KsnPTvM+D0pyapJrk9yV5M4kn05ywgL/PklycpLrkvwgybYk70yy/zTn75/kd5J8MsnWJD9sf5/1SZ4y5dyHJfleki8nyTTX+4f2szuqp+1FSS5PclOSu5N8I8k/JXnjQv7ukqSlwdxsbpaWEotoTZonAf/Svv4b4CrgpcBlSR7T7h8GnAt8DHgGcGmSh/ReJMly4ErgfwD3Ah8AzgEeDnwoyTsWMOa/BM4CHgasBdYBRwOXAQ/oc/5jgTOAne3v8OfApcCzgX9OcvSuE6vqtvZ6Pw48d+qF2j9gjgE+W1XXtG0nARcBRwL/APwZcDHwQOB1A/+2kqRJY242N0uLS1W5uS2aDdgCFPDoGc65ovmnvVvbM9t+BbxyyrH3t+23AqdNOfZ77bE3T2n/YNv+u1Pa9wM+QZMkV83xdyrgimmO/ef2+GbggCnv8+n22JYpffYHDuxzrcOAbwCbprQf1V7ngj59Tm+Pvb6n7bPA3cAj+py/x/u27a9tr/Parv8Nubm5ubkt7GZu3u19zM1ubhOwORKtSXNlVZ03pe2c9ucdwJlTjp3b/ly1qyHJCuBXgGuq6k96T66qHwBvBQL88gLEu+vb4zOq6tYp73Nqvw5VdUdVbe/TvhW4AHhMksN72q+hWaXzuPQ8wiPJMuBE4LvAh6dcbgdwT5/32ON9JUmahbnZ3CwtKi4spklzTZ+2b7Q/r62qe6cc29b+PKyn7UnAMqCSnN7nevdvfz52vkH2eEL785/6HLuSZrraHpI8FXgz8BTgEew5texQoHeRl3fTTHv7VZppcADH0vze76mq3kd7nEczTez6JOva2P61qr41x99JkqRe5uaGuVlaJCyiNWnu6NO2Y7pjVbWjXdPj/j3NK9qfT2q36TxkhmNztWuBkpuniW2Pb5eTvITmW+0f0Nxv9WXgLpppbM+kuZds3ynd1tEk39cnObOqdgIntcfeN+V9/7x93zcCv0nzaJNK8k/A77TfnkuSNFfmZnOztKhYREt7b1dC/4uqesuI3usg4MbeA0n2AQ4Etk7p84fAD4GjqmrTlD7vo0nUu6mq7yf5IPDbwPOSXEezaMlnqmpjn/PPBc5tF3H5z8BLaL4pvyTJY/zmW5I0YuZmc7M0Mt4TLe29q2i+Of4vI3ivz7U/90iuwNNopq5N9ZPA9X2S9P3aPtN5D80CI79Oc7/VMqZ80z1VVd1eVRdX1etpFnQ5AHj6TH0kSRoCc3PL3CwNn0W0tJeq6haae4+OSvJ77SIfu0nyE0mOWIC3+2D787QkB/Rcfz/gj6bpswVYmeSQnvNDs5rnkdO9UVV9CbgceCHwBuB2mqlku0nyrGmeW/mI9uf3pnsPSZKGwdxsbpZGyenc0vycDKwE/gB4VZIrae6NOoRm0ZInAScAXxnkTarqX5OcBawBPp/kApqVN48DbgNu6tPtL4D3Av+W5KPt+U/lvmdH/uIMb/lummdSHgScVVXf73POhcCdSTbQ/FEQmm/+n0TziI3L9vLXlCRpIZibzc3SSDgSLc1DVX2HZhrXGmA78DLgLcCzaB478ds0C4cshDe373MHzXSuE4BLaBLqD/vE9j6ax2/cBLwGeCXwdeDJ3DcFbTrraX4fmH662CnA1TSrk76xfa/70zw+5FlVtcfjNSRJGjZzs7lZGpVUVdcxSBMtSQH/VFXPHINYfhzYTPNYjAW7ryzJa4H/Bbyuqj64UNeVJGkYzM2SZuJ0bmk8PKNN2AAHV9U3O4rjv9FMAXvnQlwsyT8CL1iIa0mSNGLmZkl9WURL3Xv7lP07R/nmSQ4HfpnmPrLXARuB/7NAl/8Q0PtsymsX6LqSJA2TuVnStJzOLU24JM8EPkWzcueVwG9U1Y0zdpIkSUNjbpbGm0W0JEmSJElz5OrckiRJkiTNkUW0JEmSJElzZBEtSZIkSdIcWURLkiRJkjRHFtGSJEmSJM2RRbQkSZIkSXP0/wFxK4h0k2PopAAAAABJRU5ErkJggg==\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x7f667334c6a0>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"kv = Dataset('visc.nc').variables['conv_kv'][:,:,1,1]\n",
"kd = Dataset('visc.nc').variables['conv_kd'][:,:,1,1]\n",
"n2 = Dataset('visc.nc').variables['conv_N2'][:,:,1,1]\n",
"kd_int = Dataset('visc.nc').variables['Kd_interface'][:,:,1,1]\n",
"zi = Dataset('visc.nc').variables['zi'][:]\n",
"\n",
"plt.figure(figsize=(16,14))\n",
"plt.subplot(221)\n",
"plt.pcolormesh(time, -zi, kv.T)\n",
"plt.colorbar(); plt.title(r'$\\kappa_v(z,t)$ added by cvmix_conv[m$^2$ s$^{-1}$]', fontsize=20)\n",
"plt.xlabel('Time [days]', fontsize=20); plt.ylabel('Depth [m]', fontsize=20); plt.ylim([-200,0])\n",
"\n",
"plt.subplot(222)\n",
"plt.pcolormesh(time, -zi, kd.T)\n",
"plt.colorbar(); plt.title(r'$\\kappa_d(z,t)$ added by cvmix_conv [m$^2$ s$^{-1}$]', fontsize=20)\n",
"plt.xlabel('Time [days]', fontsize=20); plt.ylabel('Depth [m]', fontsize=20); plt.ylim([-200,0])\n",
"\n",
"plt.subplot(223)\n",
"plt.pcolormesh(time, -zi, kd_int.T)\n",
"plt.colorbar(); plt.title(r'$\\kappa_d (z,t)$ total [m$^2$ s$^{-1}$]', fontsize=20)\n",
"plt.xlabel('Time [days]', fontsize=20); plt.ylabel('Depth [m]', fontsize=20); plt.ylim([-200,0])\n",
"\n",
"plt.subplot(224)\n",
"plt.pcolormesh(time, -zi, n2.T, cmap=plt.cm.jet)\n",
"plt.colorbar(); plt.title(r'N$^2(z,t)$ [s$^{-2}$]', fontsize=20)\n",
"plt.xlabel('Time [days]', fontsize=20); plt.ylabel('Depth [m]', fontsize=20); plt.ylim([-200,0])"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3",
"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.6.1"
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment