Skip to content

Instantly share code, notes, and snippets.

@bradyrx
Created February 26, 2020 22:33
Show Gist options
  • Save bradyrx/51edcde290a816d52b3c26afc20d849a to your computer and use it in GitHub Desktop.
Save bradyrx/51edcde290a816d52b3c26afc20d849a to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [],
"source": [
"%matplotlib inline\n",
"import numpy as np\n",
"import xarray as xr"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [],
"source": [
"ds = xr.open_dataset('/glade/u/home/sallyz/obs_data/ersstv3b.185401-201712.nc')['sst']\n",
"coslat = np.cos(np.deg2rad(ds.lat))\n",
"ds.coords['coslat'] = coslat"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Utility function to convert to -180 to 180 lon to extract S. Atlantic box\n",
"\n",
"https://github.com/bradyrx/esmtools/blob/master/esmtools/grid.py"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {},
"outputs": [],
"source": [
"import esmtools as et"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {},
"outputs": [],
"source": [
"ds = et.grid.convert_lon(ds)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Use datetime for time axis."
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {},
"outputs": [],
"source": [
"ds['time'] = np.arange('1854-01', '2018-01', dtype='datetime64[M]')"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"1. Slice to South Atlantic"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {},
"outputs": [],
"source": [
"ds = ds.sel(lat=slice(-30, -10), lon=slice(-20, 10))"
]
},
{
"cell_type": "code",
"execution_count": 8,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.collections.QuadMesh at 0x2ae1d0f81320>"
]
},
"execution_count": 8,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAZgAAAEWCAYAAABbgYH9AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3de7zlY93/8dd7xphxGMSgwcih0UGEhpJKTpG7iFJUkg5yR3F3QhK/pFtUbtWNJpJqSgoRFWNySN2pwTiMUc41mTDjMMOMmdl7v39/XNcyy5611ve712HvNbM/zx7fx17fw3V9r7Uy61rf6/C5ZJsQQgih3UYMdQFCCCGsnKKCCSGE0BFRwYQQQuiIqGBCCCF0RFQwIYQQOiIqmBBCCB0RFUwYMEmbSnpW0sihLksIoXtFBRMKSXpY0p6Vfdv/sL2m7d6hLFcjko6WNF3SYkk/rHH+vZJmSVog6R5J76o6d4qkpbkSrWxb1MhjV0mW9NWCsmwm6XpJCyXdW/1ZShov6UpJj+a8Nivx3t4v6RFJz0n6laR1+72vP+V73VCUVwidFBVMWFk9CnwV+EH/E5I2Bn4CfAZYC/g88FNJG1Rd9vNciVa2B/vlMQo4G7ilRFl+BtwOrAecCPxS0vr5XB/wO+DdZd6UpK2B7wGHAhsCC4Fzqi55Evgf4PQy+YXQSVHBhIYk/RjYFPh1/iX/hfyL3JJWydfcIOmr+Zfzs5J+LWk9SVMkzZf01+pf5pJeKWmqpCcl/U3Se9tdbtuX2f4VMK/G6U2Ap23/1snVwHPAlgO4xWeBa4F7G10kaStgB+Bk24tsXwrcRa5QbD9m+xzgryXv+wHg17Zvsv0scBJwoKSxOb/rbF9CqmBDGFJRwYSGbB8K/AN4Z/4lf0adSw8m/aremPRF/X/AhcC6wCzgZABJawBTgZ8CGwCHAOfkX+bLkXSOpKfrbHc2+bamA7Mk7SdpZG4eWwxU5/fOXAHOlPSf/cr0MuAjwFdK3Gtr4EHbC6qO3ZGPN2PrnB4A2w8AS4CtmswvhI5ZZagLEFYaF+YvOyT9Fni17evy/i+AU/N17wAetn1h3r9N0qXAe4CZ/TO1/Ungk+0sqO1eST8iVXJjSF/QB9l+Ll9yCTAZeAx4PXCppKdt/yyf/zZwku1nJRXdbk3gmX7HniFVxM2ol9/YJvMLoWPiCSa0y2NVrxfV2F8zv34Z8PrqJxFSs89LB6eYkDvZzwDeCqwK7AqcL2k7ANv32H7Udq/tP5H6Wt6T074TGGv753Xynlk1MODNwLOkfp5qawELlk+9XF5vrsqrUvk2nV8Igy2eYEIZ7Qy5/U/gRtt7lblY0nnAB+ucfsR2M01N2wE32Z6e9/8q6RZgT2BGjesNVB5V9gAmSfp33l8b6JW0je39+5cn98FsIWlsVTPZa0lPTw3Z/gPLKuaKmTl9Jf8tgNHA34vyC2GwxRNMKOMxYLlhuk26CthK0qGSRuVtR0mvqnWx7SP7jeaq3upWLpJWkTQGGAmMlDSmMiiB1KH+5soTi6TtgTeT+2Ak7S/pJUp2Aj4NXJHTnkTq79gub1cC3wcOr1P+v5MqrZNzGQ4AtgUurSrrGFIlATA679czhdQ/9Obcn/UV4LJK5ZX7lMaQfjyOyPcc1SC/EDomKphQxn8DX8pNWp9rJaP8Rfg20qCAR4F/A19n2Rdsu3yJ1DR3POkJaFE+hu0bgVNIw4UXkL7sv2b72pz2YOB+UrPTj4Cv276oUn7b/65sOd/nbD/ZoCwHA5OAp0jDh99j+4mq84tITV+QRqUtqpeR7ZnAkaSK5nFS30t1H9WhOf25pEpzEakCDGHQKRYcCyGE0AnxBBNCCKEjooIJIYTQEVHBhBBC6IioYEIIIXTEsJsHs+66IzxhQvNve2kbxkT0uPUo90tpLY9et/7boh159FE4E76Q25DHymAkfS3nMUo9LecxekQb8qC1QN2rtmkliVvvXDzX9vrFV9a2925reN6T5d7LrXcuvsb2Ps3eqxsNuwpmwoRV+O1vxjWd/t+9rX9kT/T2nzs3cI/39J/MPTBPtqEMz/Y2mq5RzsLeVVvOY2mLFXZflzzIj2ixglh7lbqjm0sbP+rplvPYbNUnii8qsMUq81tKv+kq7YmcM3L8fY+0kn7ek7385ZpNy96r4ReTpAmkYfMvJUXhnmz77HzuU8DRQA9wte0vtFLudhl2FUwIIQwWA31teLLMeoDP2r4tR8++VdJU0rIN+wPb2l7cb9mJIRUVTAghdIgxS9u0Lp/tOcCc/HqBpFmkoKkfB063vTife7wtN2yD7mgbCCGElVRfyf8B45RWYa1sR9TLM6+vtD1pwbutSKGPbpF0o6QdB+N9lRFPMCGE0CHG9JaPljLX9qSiiyStSQpvdKzt+TnG3kuANwA7ApdI2sJdEKYlKpgQQuigvjYGI8+BSy8Fpti+LB+eTQp4auAvkvqAcUDroy1aNGRNZJIOymtn9Ema1O/cCZLuz8vp7l0n/bp52d378t+XDE7JQwihHAO9uNRWRGl1uwuAWba/VXXqV8Du+ZqtSGsczW3/uxm4oeyDuRs4ELip+qCkV5Oiz24N7ENaTrfWONTjgWm2JwLT8n4IIXSVPlxqK2EXUrTs3SXNyNu+wA9Iaw7dDVwMHNYNzWMwhE1ktmcB1Fhydn/g4jwi4iFJ9wM7kdZ473/dW/Pri4AbgOM6VNwQQhgwA0vb9F1v+2aoO6u43qJ8Q6obR5FtTFr1sGI2tdcv3zAP26sM36s79lvSEZWRGfPmtW1MegghNOSSzWNlmshWRB19gpF0HbXXWj/R9hU1jkPtGrqlT9/2ZGAywGtfu+rK+f9kCKH7GHqH8TdORysY23s2kWw2MKFqfxPSyof9PSZpvO05ksaTVvcLIYSukWbyD1/d2ER2JXCwpNGSNgcmAn+pc91h+fVhLFszPYQQuoToLbmtjIZymPIBkmYDOwNXS7oGXlhz/BLgHuB3wFF2irUg6fyqIc2nA3tJug/YK++HEELXSJ38KrWtjIZyFNnlwOV1zp0GnFbj+MeqXs8D9uhYAUMIoUVpHszKWXmUETP5Qwihg/pW0qeTMqKCCSGEDoknmBBCCB1hRG9XjqUaHFHBhBBCB0UTWQghhLYzYkmLS3qvyIZdBfO8R3Df0ubXo//H0nVbLsOcpa0Hfn586VotpX9icfOfQcX8pWNazmNhz6iW81jS29o/4F633oTR2zf0eUitTxkfPbKn5TzWHv18y3lstsa8ltK/cez9LZehHdJEy2giCyGE0AHRyR9CCKHtbLXlCXlFFRVMCCF0UF88wYQQQmi31Mk/fL9mh+87DyGEDotO/hBCCB3TG/NgQgghtFvM5A8hhNAxfTGKLIQQQrulYJcrZgUj6cASlz1v+zf1TkYFE0IIHWLE0hU3VMz3SSsFN+pEegvQfRWMpIOAU4BXATvZnp6PV1anXBVYAnze9u9rpD8F+DjwRD70xUY1aQghDDa7PaGIhshvbX+k0QWSftLo/FC+87uBA4Gb+h2fC7zT9jbAYcCPG+Rxlu3t8haVSwihy4i+klthTtIESddLmiVppqRj8vFTJP1L0oy87duOktv+YKvXDOWSybMAJPU/fnvV7kxgjKTRthcPYvFCCKFlpq1PMD3AZ23fJmkscKukqfncWba/0a4bAUj6ICDbP+53/OPAc7Z/WpRHtz+7vRu4vUHlcrSkOyX9QFLdEMWSjpA0XdL0Z+b1dqakIYRQQy8jSm1FbM+xfVt+vQCYBWzcwaJ/FvhVjeMX53OFOlrBSLpO0t01tv1LpN0a+DrwiTqXnAtsCWwHzAG+WS8v25NtT7I9ae31VtgOtxDCCsaIPpfbgHGVH8J5O6JevpI2A7YHbsmHSv3YHqCRuSJ78XtKx0qts9HRJjLbezaTTtImwOXAh2w/UCfvx6qu/z5wVVOFDCGEDjGwtHwssrm2JxVdJGlN4FLgWNvzJZ0LnJpvdyrpx3bDzvmSRklaw/Zz/e4/ljQIq1DXNZFJWge4GjjB9h8bXDe+avcA0qCBEELoIqK35FYqN2kUqXKZYvsySD+2bffa7iMNLd6pTYW/APhlflqq3H8zUhPZBWUyGLIKRtIBkmYDOwNXS7omnzoaeDlwUtWoiA1ymvMlVWr4MyTdJelOYDfgvwb7PYQQQiMmzeQvsxVRGhF1ATDL9reqjnfkx3YeNHAFcKOkeZLmAjcCV9k+s0weQzmK7HJSM1j/418FvlonzceqXh/azH17PJLHe8c2kxSA2W1YMvmRReu1nMdjz7e4ZPKiNVouw4LnR7ecx6LFpZ60G+rtae13kntbD0bovnbk0eLvvba8j5azQG34LGaM3LSl9Fet8ZqWy5BMbzmHNq5ouQtwKHCXpBn52BeBQyRtR6rPHqZ+v/WA2T4POC83y6lWn4ykw2xfVCt9zOQPIYQOsdW2WGS2b6b2rPrScwAlTQT+G3g1MKYq7y0K7v1sg9PHAFHBhBDCYEqd/F01cvVC4GTgLFLXwuE0DgVTRt30XdfJH0IIKw/R6xGltkGymu1ppOauR2yfAuzeYp6udyKeYEIIoUNSJ39XLTj2vKQRwH2Sjgb+BWzQYp7xBBNCCEOhXTP52+RYYHXg08DrgA+SYj62ou50kniCCSGEDqnM5O8Wtv8KIMm2Dy+TRtKGwNeAjWy/XdKrgZ1tX5DzPLpe2niCCSGEDupjRKltMEjaWdI9pDhmSHqtpHMKkv0QuAbYKO//nfQkVCgqmBBC6BAblvaNKLUNkv8B9gbmpfL5DtKiYY2Ms30J0JfT9AClogY3bCKT9O0Secy3/aUyNwshhOEkNZF11+942//st0xKUWXxnKT1yKPFJL0BeKbMvYr6YPYHvlxwzfFAVDAhhFBDG2fyt8M/Jb0RsKRVSZ39swrSfAa4EthS0h+B9YH3lLlZUQVzVr0QABVtDA0dQggrlS4cpnwkcDZpHZnZwLXAUY0S5AXOdgVeQRqS/DfbS8vcrGEFY/t/ijIoc00IIQxP3dVEZnsu8IGBpJG0Oukp5mW2Py5poqRX2C5cIqXUO5e0vqQvSpqcF7T5gaQfDKSQIYQwHPWhUttgkHSGpLUkjZI0TdLcvDRyIxcCS0iR7yE9+dQMSNxf2XkwVwB/AK6j5OiBEEIY7tIosq6KRfY221+QdACpojgIuB74SYM0W9p+n6RDAGwvUr9RAvWUrWBWt31cyWtDCCHQfRMtWbbU8b7Az2w/WaKuWCJpNZaNItsSWFzmZmUbB6+StG/Ja0MIIWTd1EQG/FrSvcAkYJqk9YHnC9KcDPwOmCBpCjAN+EKZmzWsYCQtkDSfFO//KkmLJM2vOt40SQdJmimpr2qVSiRtlu9TWc3yvDrp15U0VdJ9+W+MZgshdJXKKLIy26CUxz6e1JcyKY8EW0iajlJTbgq7FzgQ+DDws5z2hjL3KxpF1vzSj8XuJhX6ezXOPWB7u4L0xwPTbJ8u6fi8H814IYSu0k2jyABsP1X1+jnguQbXWtKvbL8OuHqg9yo7imxamWMDYXuW7b+1kMX+LFtF7SLgXa2UJ4QQ2s0WPR5Rautif5a0YzMJi0LFjAHWAMblJqjKc9xaLAt81gmbS7odmA98yfYfalyzoe05ALbnSKq7poGkI4AjAMZt1Noa8Ev7Wg9Avai39XXon13aWh4Ll7RehoWLRrecx9JFbQjovbjFUTpLW2+eUG/reYxocS179bRcBEYsaf19jGq0uG5JY+a1ln71J9rzhf1gG/Losk7+ZuwGfELSI6SnHZEebrYtSlj0r/sTpKiZGwG3sqyCmQ/8b1Hmkq4DXlrj1Im2r6iTbA6wqe15kl4H/ErS1rab7vOxPRmYDLDlNmvUXX0thBDaqVtm8kvaodF527c1OP32Zu9b1AdztqTvAl+0fepAM7e9ZxNpFpOHwNm+VdIDwFbA9H6XPiZpfH56GQ88PtB7hRBCp3VDBQN8M/8dQxpBdgfpgWFb4BbgTQ3SLih5bDmFz5G2e0ljpgdFjhowMr/eAphI7SfVK1m2EtthpMmgIYTQNSrzYIZ6FJnt3WzvBjwC7GB7Uu643x64vyD5bcATpHVg7suvH5J0W25lqqtsQ+W1kt5ddvZmGZIOkDSbNGTuaknX5FNvAe6UdAfwS+BI20/mNOdXDWk+HdhL0n3AXnk/hBC6SpfNg3ml7bsqO7bvBopG7P4O2Nf2ONvrkZrMLgE+CTRcrKxsD+tnSJ39vZIWsayTZ62S6Zdj+3Lg8hrHLwUurZPmY1Wv5wF7NHv/EELoNBt62rSYmKQJwI9I/dp9wGTbZ1ed/xxwJrB+DmpZyyxJ55NCwxj4IMXh+ifZPrKyY/taSV+z/RlJDUf6lKpgOjwfJoQQVlptbP7qAT6bw+ePBW6VNNX2Pbny2Qv4R0EehwP/SZo8D3ATcG5BmiclHQdcnPffBzyVuzL6GiUsPUZU0n4sW1rzhjKhmkMIYThrZyyyPC2jMjVjgaRZpHVd7gHOIoVvadgXbfv5HB3lNwOYh/h+UriYX+X9m/OxkcB7GyUsVcFIOh3YEZiSDx0j6U057EAIIYQ6XL6CGSeperTs5DzFYjmSNiN10N+Sf/z/y/YdRd3k+dozgVVJ8w23A75ie7/65fdc4FOS1rTdf5ZTwwECZZ9g9gW2s92XC3kRcDspPEsIIYQ6BtCBP9f2pKKLJK1J6qc+ltRsdiLwtpL3OBnYCbgBwPaMXFk1ut8bgfOBNYFNJb0W+ITtTxbdbCC9T+tUvV57AOlCCGFYstsb7FLSKFLlMsX2ZcCWwObAHZIeBjYBbpNUa4I7QI/tZwb4Ns4C9gbmpffkO1jWXdJQ2SeY/wZul3Q9aQTZW4ATBljIEEIYZkRv+0aRCbgAmGX7WwB5yPEGVdc8TBr1VW8U2d2S3g+MlDQR+DTwp6J72/5nv+a3UgtPlnrntn8GvAG4LG872764caoQQgi2Sm0l7AIcCuxetZzJQCfBfwrYmhQt5afAMywbUVbPP3MzmSWtmodDFw1tBgYwioxUGc3NabaStJXtmwaQPoQQhpV2xiKzfTM07tCxvVlBNv9h+0RSvw2Q1uYCftEgzZHA2aQRa7OBa4GjShS59Ciyr5PGPs9k2bhnk8ZQhxBCqMWpH6aLnMDylUmtYy/IzW0faOZmZZ9g3gW8IgeiDCGEUNIghoGpS9LbSaOBN5b07apTa5FGotVK8x3Sg0RNtj9ddN+yFcyDwChylOMQQgjF3MZO/hY9SopIvx9p6ZWKBcB/1UlTmZOzC/Bq4Od5/6B+edRVtoJZCMzIq1i+UMmUqcFCCGE464Ymsjy0+A5JlwPP5Sj55HAvNeOJ2b4oX/NhYDfbS/P+eaR+mEJlK5gr8xZCCGEABjCTfzBcC+wJVGbkr5aPvbFBmo2AscCTeX9NSq5oXDbY5UWNzku61Pa7y+QVQgjDhd11FcyY6nAvtp+VtHpBmtNZNg8SYFfglDI3a8OC6ABs0aZ8Om4kfawzYmHT6ddepfm0FWut8nzLeTw5sui/icZGqmEQ1FKkNjz7t+MfX29reYxY0nobeTvWsh9Rs6u1PLWYHmBk6/9pMvrp1vNY54HWuntXufHO1gvRJl2yomXFc5J2qCyRnBcMW9Qoge0LJf0WeH0+dLztf5e5Wbt6nwb8TSPpIEkzJfVVLSKGpA9UTSKakc8vtyCOpFMk/auFCUchhNBxdrltkBwL/ELSHyT9gdRxf3StC6vDzdj+t+0r8vbvWtfU0q4nmGbcDRwIfK/6oO0p5KjNkrYBrrA9o04eZ9n+RkdLGUIITTKirztGkQFg+6+SXgm8gjRp895K530NvwF2KMiy4TXtqmAG/AxoexZAQXjpQ4CfNVmmEEIYcl0wiOwFkj7U79D2krD9oxqXv1bS/EbZAY3OD7yCkfQSYILt6kbO4waaT0nvA/ZvcP7o/IFNJ6309lStiyQdARwBsMFGQ/nQFkIYVrqvk3/HqtdjSMvO30ZaivlFbI9s9WZlQ8XcQJqgswowA3hC0o22P5MLUnNMtKTrSOtH93ei7YYrr0l6PbDQ9t11LjkXOJX0A+FU4JvAR2pdmBftmQyw1TarddMPihDCyq6LvnFsf6p6X9LawI87db+yP+fXtj1f0seAC22fLKlwmIbtPVso28E0aB6z/VjltaTvA7GEcwih63TZE0x/C4GJncq8bAWziqTxpPWXTyy6uFWSRpDCEdRd1EbS+LxGNcABpEEDIYTQNQz09XVPBSPp1yx7phpBCgFzSafuV7aC+QpwDfDHPAphC+C+Vm4s6QDgO8D6wNWSZtjeO59+CzDb9oP90pwPnGd7OnBGHr5s4GHgE62UJ4QQ2s60Z65X+1SPuu0BHrE9uyiRpDcBE/OcmPWBNW0/VJSu7Ez+X1AVzjl/8bc0c9/25cDldc7dQFrgrP/xj1W9PrSV+4cQwmDohlhk8ELcsZMG2nUh6WRgEmlo84WkwMc/IQXBbKjUAG1JW0maJunuvL+tpC8NpJAhhDAsueTW6WKkAJcLc8f+QBxAGuT1XM7nUVJsskJlm8i+D3yePCnS9p2Sfgp8dYAFDSGEYaT0csiD5XngLklTyRUGFEbGX2LbyrGhJK1R9mZlK5jVbf+l36TINkQ+CiGElVyXNJFlV+etWlEJL5H0PWAdSR8nTQf5fpmbla1g5kraslIQSe8B5jROEkIIw5zBXTSKDFjH9tnVByQd0yiB7W9I2os0a/8VwJdtTy1zs7IVzFGkiYqvlPQv4CHggyXThhDCMNZVFcxhwNn9jn24xjHghYEB1+SBAaUqlWplR5E9COyZ295G2F4w0BuFEMKw1AVNZJIOAd4PbC6pevHItYB59dLZ7pW0UNLatp8Z6H3LhorZEPgasJHtt0t6NbCz7QsGesMQQhhWuqCCAf5E6tYYRwqrVbEAKIrK0szAAKB8E9kPSeOfK7P4/05aRyAqmBBCqKeNEy0lTSAFpXwp0AdMtn22pFNJQYH7gMeBD+ehxMuKYT8CPCJpT2CR7T5JWwGvBO4quHWtgQGllK1gxtm+RNIJubA9knqbuWEIIQwnbZxo2UOKGn+bpLHArfmp4kzbJwFI+jTwZeDIOnncBLw5R8WfRopE/z7gA/XL74uaLXDZCuY5SeuxbBTZG4ABt8eFEMKw06ZRZDn24pz8eoGkWcDGtu+pumwNGjfKyfZCSR8FvmP7DEm3N7qvpIdq5Wl7i6Iyl61gPgNcCWwp6Y+k+GHvKZm2q4xSLxuMfLbp9AtG1e0PK23hmFVbzmNR36jWytDTWnqAxT2tr63Ts7TlJSfoXdriioEjWv+JqTZ8iYyot65g2fRLWi4CIxe3nkc72jbceCHC4vQ9LX6YbaQO9MFI2gzYHrgl758GfIj0w3+3xkm1M+mJ5aP5WNE/5ElVr8eQAhGvW6achf8yc2TjMcCuwBtJQSW37rfgWAghhP7KholJldA4SdOrtiNqZSlpTeBS4Fjb8wFsn2h7Amm5+aMblOgY4ATgctszc+Di6xu+BXte1fYv2/8D7F7m7Rf+BM2dQd+0vTMws0ymIYQQADSQTv65tic1ukDSKFLlMsX2ZTUu+SmpQ/7kWult30Tqh6nsPwg0HA0maYeq3RGkJ5q2xiK7VtK7gcvsbokNGkIIK4A2fWMqxeq6AJhl+1tVxyfariyfsh9wb4M8tgI+B2xG1fe/7UZPJNXDmntIE+3fW6bMA+mDWQPokfQ8aWqqba9VMn0IIQxPfW3LaRfgUNKclBn52BeBj0p6Rb7TI9QfQQZp2ZXzgPOBsr1lH62xNtfmZRKWnclf6nEohBBClTbOg7F9M7XjzvxmANn02D53gLf+JbBDjWOvK0pYdiZ//8whjVZ4xHZTUZUlnQm8E1gCPAAcbvvpfO4E0giHXuDTtq+pkX5d0mTPzUgrWr7X9lPNlCWEEDqlE6PIWvBrSZ8kLfb4wphB20/2v1DSK4GtgbUlHVh1ai3SwK9CZZvIziHVYJUZn9sAdwDrSTrS9rUl86k2FTghT9r8Omlkw3E5DM3BpDe2EXCdpK3yYjnVjgem2T5d0vF5/7gmyhFCCJ3TXRXMYfnv56uOGag1p+UVwDuAdUgPAxULgI+XuVnZCuZhUjvcTIBcCXweOBW4DBhwBdOvUvozy+bV7A9cbHsx8JCk+4GdgP/rl8X+wFvz64uAG4gKJoQQ6rJdqu8kX3sFcIWknW33//4tpWwF88pK5ZJvfI+k7W0/qBYnRGUfITV3AWxMqnAqZudj/W2YZ7Zie46kDeplnseTHwEwfuPWJ/aFEEJZ3dREloc5/yfwlnzoBuB7thvNTL1d0lGkVqUXmsZsf6TofmWnQP9N0rmSds3bOcDfJY0G6hZM0nWS7q6x7V91zYmkoW9TKodqZNXS/0W2J9ueZHvSS9ZtcdZ3CCGUZVKomDLb4DiX1Dl/Tt5el4818mNSgM29gRuBTUjNZIXKPsF8GPgkcCypAriZNJZ6KQ3CEuRFauqSdBipjW+Pqvk1s4EJVZdtAjzaPy3wmKTx+ellPCmKaAghdJcueoIBdrT92qr930u6oyDNy20fJGl/2xdJ+imw3MCrWkr9nLe9iFTbHW/7Xba/YXuh7T7bTQX2krQPqc9kP9sLq05dCRwsaXQeaz0R+EuNLK5kWYfVYcAVzZQjhBA6SS63DZJeSVu+ULYUKqZoPkylleppSa8B1iaN3i1UqoKRtB8wA/hd3t+u36pozfguKdzAVEkzJJ0HkPt6LgHuyfc7qjKCTNL5kiqhFE4H9pJ0H7BX3g8hhO5SPhbZYPg8cL2kGyTdCPwe+GxBmsk5vP+XSD/s7wG+XuZmZZvITiaN5LoBwPaMHM2zabZf3uDcacBpNY5/rOr1PGCPVsoQQggd10VNZLanSZpIGoIs4N48YremHOx4fp5jeBO1hzPXVbbHu6eZ9ZhDCGE4K9s8NlhNZHk02Gq277R9B7B6nnhZk+0+GkdnbqhsBXO3pPcDIyVNlPQd0hrPIYQQGumuUWQfr0RMAchPJkWTJqdK+pykCZLWrWxlbla2iexTwImk0AI/I40gOLVk2hBCGLa6aR4MMEKSKqN2JY0EilZArMx3OarqWITtuzEAABnsSURBVL3Z/y9SNtjlQlIFc2KZ60MIIWTdVcFcA1ySB1WZFHn5d40SDGT2f38NKxhJv6bBx2N7v2ZvPFRGYMaOaCo+JwDrjXyu5TJsuMr8lvN4dtVSsebqWrR660sm9/a1PmnVbYg02+r/I31u/bNQTzuWTG4xj3Z8kbUjjzbMZV6ydmvLcY+ZuGXxRWX8vcX0gzsEuYzjSKsS/yepk/9aUuj+uiStTlqyZVPbR1QGCdi+quhmRf8vfiP/PZA0k/Mnef8QUnyyEEIIjXRRBZNXKP4h8HvbfyuZ7ELgVuCNeX82aV2Zwgqm4W8N2zfavhHY3vb7bP86b+8H3lSycCGEMGypr9w2KGVpbk7jlrbPIE+4zBPvSz1ul32YXT/P+KwUcnNg/ZJpQwghdIfKnManIc1ppHhW/hJJq5GfxXIkgLpzZ6qVbej8L+AGSZVlMzcjRycOIYTQQBc1kZHnNA4wCv7JpCeeCZKmkJZu/nCZhGVHkf0ud+y8Mh9qOPszhBAC3djJ/6I5jcCnKZjTaHuqpNuAN5Caxo6xPbfMzRo2kVUvlWx7se078ra41jUhhBD66a5YZJ8ireuyGPgp8AwpSn6RXUmhuXYD3lz2ZkVPMBdKeiuNO3QuALYve8MQQhhWuuQJJk+q/H+2P88A5jTm9b9eTppkD/AJSXvaPqpBMqC4glmbNDytUQXzRKlShhDCMCMGb4RYEdu9kl7XRNJdgddUzf6/CLirTMKGFYztzZooTAghBOjGPpjb87DkX1A1R9n2ZQ3S/A3YFHgk708A7ixzs9amy4YQQmisuyqYdYF5wO5Vxww0qmDWA2ZJqiz8uCPwf5X5M40iukQFE0IIndSmCkbSBOBHpKgqfcBk22dLOhN4J7AEeAA4vDpi8ouKYh/exK2/3GSR2xE1qDmSzpR0r6Q7JV0uaZ18fC9Jt0q6K//dvU76UyT9K6+GOUPSvoP7DkIIoVgb14PpAT5r+1WkIcNHSXo1MJXUR7ItKXraCe0sf1VEl9tJfS93AXdVHa+r7JLJkvRBSV/O+5tK2qnFctf7UOYC77S9DXAY8OMGeZxle7u8/abF8oQQQvu1aZiy7Tm2b8uvFwCzgI1tX2u7EsH3z8Am7Sy+pCMkPUbqd5lOGvg1vUzask1k55AeyXYHvgIsAC4ltcU1xfa1Vbt/Bt6Tj99edXwmMEbS6JjYGUJY4XhAo8jGSar+4p5se3KtC/OS9dsDt/Q79RHg5wMrZKHPA1uXnVxZrWwF83rbO0i6HdIqaJKKFqkZiHofyruB2xtULkdL+hCpNv1sXp1tOZKOIIe22WjjIWsVDCEMR+X7YObanlR0kaQ1ST/wj7U9v+r4iaRmtCk10nymYRHtbzU4/QCwsKhctZStYJbmSTqVcdDrk55oGpJ0HalDqr8TbV+Rr6n5oUjaGvg68LY62Z9LWlXT+e83Wbby2ovkXwGTAbbZdlR3jekIIazU2jlMWdIoUuUypXposaTDgHcAe1Tmq/QzNv99BanlqRJB+Z3ATQW3PQH4k6RbqApyafvTReUtW8F8G7gc2EDSaaTmrC8VJbK9Z6Pz9T4USZvk+33I9gN18n6s6vrvU2JtghBCGHTtG0UmUuSUWdVPHJL2IS0ktmtefXj5Itj/L197LbBD7sNB0imkOTGNfA/4Palzf0DTRssGu5wi6VZSLBoB77I9ayA36q/eh5JHk10NnGD7jw3Sj7c9J+8eANzdSnlCCKHt2htnbBfgUOAuSTPysS+SHgBGA1NzlOQ/2z6yTh6bkoYzVyyhOFx/j+2GTWz1FC2ZvG7V7uMsi0WDpHVtP9nMTbPvUvtDOZoU9+YkSSfla99m+3FJ5wPn2Z4OnCFpO9L/fQ+TlgENIYSuIdrXRGb7ZmqH7RrICNofA3+RdDnpu/MA0tyaRq7P/di/5sVNZIXf/0VPMLfmQohU8z2VX68D/APYvOgG9dh+eZ3jXwW+Wufcx6peH9rMffsQz3lkM0kBWNJC2oqRbQhOtOYqz7eUftyqz7ZchiWrrRzzdJ9rwxdAb9+A1teora+1ASjqbb0MI9sxVrMdn+fo1t7L0peu1XohIE2gaFE3hYqxfZqk37IsIvLh/Ubu1vL+/Ld6fo2BLWpc+yJFscg2B5B0HnBlZa6JpLcDDftXQggh0G2hYgBWB+bbvlDS+pI2t/1QvYsr9UAzyv5k2rF6IqPt35IibIYQQmiki9aDkXQyqe+78jQyCvhJQZrVJX1J0uS8P1HSO8rcr2wFMzffYDNJL8tDi+eVTBtCCMNTyTAxg9iMdgCwHzmSsu1HWTaEuZ4LSYMB3pj3Z1OnG6O/shXMIcD6pKHDvwI2yMdCCCE00kVPMMCSPCWkMqdxjRJptrR9BrAUwPYiGq8R9oKyw5SfBI4pc20IIYRlumXBsewSSd8D1pH0cdLk9O8XpFkiaTWWVUpbUjWarJFSFYyk66lRx9quGek4hBBC0mWjyL4haS9gPmlW/5dtTy1IdgrwO2CCpCmk+Tilwv6XHWf6uarXY0gxwnrqXBtCCAEGu/mrlFyhFFUq1ddfmyfav4HUNHZM2cCXZZvIbu136I+SGq4DEEIIga6oYCQtoEFJbNedOCRpmu09SBFW+h9rqGwTWfWM/hHA66gdxDKEEELWzpn8rbA9FkDSV4B/k2b0C/gAdUaRSRpDmjMzTtJLWNaxvxawUZn7lm0iq57R3wM8BHy0ZNoQQhi21NcFNcwye9t+fdX+uTlK8hk1rv0EcCypMrmVZRXMfOB/y9ysbAXzKtsvik0iaXTJtCGEMDx1Xx9Mr6QPABeTSnYI0FvrQttnA2dL+pTt7zRzs7LzYP5U49j/NXPDEEIYTrpsouX7gfcCj+XtIJbFGqup2coFiqMpvxTYGFhN0va8uA1u9WZvGkIIw0YXPcHYfhjYf7DuV9REtjfwYWAToHpJzQWkdQhCCCE00A2d/BV5NeKPk9aAeeH733bN1YBbVRRN+SLgIknvtn1pJwoQQggrtS6qYIArgD8A11Gn76W/vJLmB4AtbH9F0qbAS23/pShtURPZB23/BNhM0nIrmlUv2zlQks4krQe9BHiAtC7B05I2A2YBf8uX1lydLQ+d/jmpJn4YeK/tp5otTwghtJ27LlTM6raPG2Cac0hLJe8OfIXUgnUpsGNRwqJO/kogtDVJY6WrtzUHWMj+pgKvsb0taVmf6sVsHrC9Xd7qLf15PDDN9kRgWt4PIYSuUZkH00Wd/FdJ2neAaV5v+yjgeYD8Q37VMgmLmsi+l19eZ/uP1eck7TLAQvbP+9qq3T8D7xlgFvsDb82vLwJuIK1zEEII3cNd1UZ2DPBFSUtIrUcC3GgmP7BU0kiWBbtcn/REU6jsMOVaw9SaHrpWw0eA31btby7pdkk3SnpznTQb2p4DkP9uUC9zSUdImi5p+lNPdtfzaghh5dZNTzC2x9oeYXuM7bXyftH60t8mLdWygaTTgJuBr5W5X1EfzM6kRWbW79cHsxZQuDi9pOuoHVLmRNtX5GtOJEUHmJLPzQE2tT1P0uuAX0na2vb8wndTh+3JwGSALbZZw3ct3rjZrHi6t/XR2c/0rtZyHk8tLbOMQ32Leks94TY0og3/Kkav0nrM1NVWXdpS+sWrlp1vXN+SVVv/4eJVW1uHvm+V1tIDuPBfdbF29Dmot7X/tnrHtOGNtEOXTbSs6rDf3PapkiYA4xt12NuekoNd7kF64nmX7Vll7lf0L2tVUl/LKrw4Xs18SjRp2d6z0XlJhwHvAPbIi+BgezF5rQHbt0p6ANgKmN4v+WOSxtueI2k88HhReUIIYbB1WSd/dYf9qcCzpLAvdTvsJb0BmGn7f/P+WEmvt31L0c2K+mBuBG6U9EPbj5R/D8Uk7UPqM9nV9sKq4+sDT9rulbQFMBF4sEYWVwKHAafnv1e0s3whhNAOXVbBvN72DpJuh9RhL6moOeNcYIeq/edqHKupbB/MQklnSvqNpN9XtpJp6/ku6aloqqQZks7Lx98C3CnpDuCXwJF5RU0knS9pUr7udGAvSfcBe+X9EELoHiZ18pfZCkiaIOl6SbMkzZR0TD5+UN7vq/p+rKeZDntVWpgAbPdRMo5l2cbnKaQ5J+8AjiQ9MTxRMm1Ntl9e5/ilpDHWtc59rOr1PFKbYAghdK02duD3AJ+1fZukscCtkqYCdwMHAt9rmDrp32H/HuCkgjQPSvo06akF4JPUblVaTtkKZj3bF0g6pqrZLBYcCyGEIm2qYPJo2crI2QWSZgEbV5Y8Tv33hXk002F/JKli+hLp3UwDjihT5rIVTGWYzhxJ/wE8SopPFkIIoY4BLjg2TlL1YKbJeQTs8vmmiCfbA4Ud7f3S/dj2ocC9NY7VZPtx4OCB3KeibAXzVUlrA58lzX9Zi7QQTQghhHrsgSw4Ntd2UR8KktYkdSMc28T0ja375TWStEJxo/uNIS0wuTUwpnK8TIDMUp38tq+y/Yztu23vZvt1wJZl0oYQwrDmklsJkkaRKpcpti8rWwRJJ0haAGwraX7eFpCmd1xZkPzHpPmMewM3klqvFpS5b9lRZLUsF/wyhBDCi7VrJn+eJHkBMGuggYZt/7ftscCZeQZ/ZRb/eraL4ji+3PZJwHM5wv5/ANuUuW8rFUzr04ZDCGFlZqDP5bZiuwCHArvnqR0zJO0r6QBJs4GdgaslXdMgj/urdySNlHRywX0rffBPS3oNsDYpin2hVmJkdFEAhBBC6FLtG0V2M/V/2F9eMps9JL2b1KcyDvgBqdmrkcmSXkIaRXYlKbpL0dBmoDgW2QJqfzwCWg+oFUIIK7luWtHS9vslvQ+4C1gIHNI/Un5FnpZyNqlJ7ingJmCLgdyvYRNZJdJmjW2s7dYjBIYQwkpOfS61DUpZpImkkP2XkhZqPFRSvQi+h+e/TUfOj0oihBA6pcuiKQO/Bo6yPS0PGvgM8Ff6DV/OZkl6mDTr/86q45U1ZLYtullUMCGE0CFpomVX1TA7VebO5Phi35RUc5iy7UMkvRS4BtivmZu1MooshBBCkb6SWwdJ+gKA7fmSDup3+vAaSSqeAO6y/Uj/rcx9o4IJIYQOkl1q67DqUC8n9Du3T71EtntJIWyaWqEwmshCCKFTuqcPRnVe19rv7xHgj7kp7bnKwTKTPaOCCSGEjhm8EWIFXOd1rf3+Hs3bCF68snGhYVfBPL10DS5/vHAhtrr62hDA4Pme1j/253tby2NpX+trli9pw/tY0tt6ORYtHtVS+qXPt+GfwfOttzaPWNzaf1sjelsuQnt+bbchjxE9LRahm+KMdEcn/2slzSfPYcyvyftj6icD2/+v2ZsOWQUj6UzgncAS4AHgcNtPS/oA8PmqS7cFdrA9o1/6U4CPs2zhsy/a/k3HCx5CCGW5O5ZMtt30LzlJ11PjZ4Pt3YvSDuUTzFTgBNs9kr5O6ng6zvYU0gqaSNoGuKJ/5VLlLNvfGJzihhBCE7rjCaYVn6t6PQZ4N2l1zUJDVsHYvrZq98+kpTv7OwT42eCUKIQQOmAFr19s39rv0B/LrmjcLX0wHwF+XuP4+4D9G6Q7WtKHgOmktaqfqnWRpCPIS3yuvuGaLRY1hBDKU18XtJG1QNK6VbsjSAuUvbRM2o5WMJKuq1OQE21fka85kfS4NaVf2tcDC23fXSf7c4FTSb8PTgW+SaqolpOXHZ0MsO4rN1jBf0+EEFYYpuOTKAfBraR3ItJ39UOkaMyFOlrB2N6z0XlJhwHvAPbIYQuqHUyD5jHbj1Xl833gqhaKGkIIbScGZRJlR9nevNm0QzaTX9I+wHHAfrYX9js3AjgIuLhB+vFVuwcA9Z50Qghh6Njlti4l6SBJY/PrL0m6TFKpuR5DGSrmu6RJO1PzymznVZ17CzDb9oPVCSSdL2lS3j1D0l05yuduwH8NSqlDCGEgVvAKBjjJ9gJJbwL2Bi4idVEUGspRZC9vcO4G4A01jn+s6vWhnSlZCCG0ycrRB1OZwvsfwLm2r8jzEAt1yyiyEEJYKa3oo8iAf0n6HrAn8HVJoynZ+hXRlEMIoWNKNo91dxPZe0lrwuxj+2lgXV4cbaWueIIJIYROMd1eeRTKg7Auq9qfA8wpkzYqmBBC6KQVvoWseVHBhBBCB63o82BaEX0wIYTQSW3qg5E0QdL1kmZJminpmHx8XUlTJd2X/76k4++ppKhgQgihU2zo7Su3FeshxVx8FWkax1GSXg0cD0yzPRGYlve7QlQwIYTQSW16grE9x/Zt+fUCYBawMSkg8EX5souAd3XonQxY9MGEEEInle+DGSdpetX+5ByodzmSNgO2B24BNswju7A9R9IGzRe2vaKCCSGETjHQV7qCmWt7UtFFktYELgWOtT1f6qb1oV9s2FUwzz2/Kn+ducXQFqINg0rcauPmiDYUoq/1/7C1tA159LSWx8jnWy/DKgtbz2PkwuJrGpZhcctFYGQ78lja+n9b6m0tj25YpjgxuH2FkTSKVLlMsV2Zm/KYpPH56WU88Hjbbtii6IMJIYROMW3r5Fd6VLkAmGX7W1WnrgQOy68PA65o99to1rB7ggkhhEHVvnkwuwCHAndJmpGPfRE4HbhE0keBf5CWOukKUcGEEEIntamCsX0zaVXJWvZoy03aLCqYEELomK4PZNlRUcGEEEKnGFjxw/U3bSiXTD5V0p15NctrJW1Ude4ESfdL+pukveuk79rwCCGE8IIVP1x/04ZyFNmZtre1vR1wFfBlgBz64GBga2Af4BxJI2uk79rwCCGEkLQ1VMwKZ8gqGNvzq3bXYNnskP2Bi20vtv0QcD+wU40sujY8QgghALkLpq/UtjIa0j4YSacBHwKeAXbLhzcG/lx12ex8rL/S4REkHQEcATBy3XXaUPIQQiip/Ez+lU5Hn2AkXSfp7hrb/gC2T7Q9AZgCHF1JViOrlv4fsj3Z9iTbk0aOXaOVrEIIYWCGcR9MR59gbO9Z8tKfAlcDJ5OeWCZUndsEeLRGmq4NjxBCCECqOGIU2eCTNLFqdz/g3vz6SuBgSaMlbQ5MBP5SI4uuDY8QQggviCeYIXG6pFeQVqx+BDgSwPZMSZcA95AW2DnKdi+ApPOB82xPp4vDI4QQQmLc2zvUhRgyQ1bB2H53g3OnAafVOP6xqtfz6NLwCCGEAAw0XP9KJ2byhxBCJ62kQ5DLiAomhBA6xIDjCSaEEELbub0Ljq1oooIJIYQOGs6d/PJKOjyuHklPkEat1TMOmDtIxWlFlLO9opztsyKUEcqV82W212/2BpJ+l+9Txlzb+zR7r2407CqYIpKm25401OUoEuVsryhn+6wIZYQVp5wrsqGMphxCCGElFhVMCCGEjogKZnmTh7oAJUU52yvK2T4rQhlhxSnnCiv6YEIIIXREPMGEEELoiKhgQgghdERUMJmkMyXdK+lOSZdLWqfq3AmS7pf0N0l7D3E5D5I0U1KfpElVxzeTtEjSjLyd121lzOe65rPsT9Ipkv5V9RnuO9RlqpC0T/7M7pd0/FCXpx5JD0u6K39+04e6PBWSfiDpcUl3Vx1bV9JUSfflvy8ZyjKujKKCWWYq8Brb2wJ/B04AkPRq4GBga2Af4BxJI4eslHA3cCBwU41zD9jeLm9HDnK5qtUsYxd+lrWcVfUZ/maoCwOQP6P/Bd4OvBo4JH+W3Wq3/Pl10xyTH5L+m6t2PDDN9kRgWt4PbRQVTGb7Wts9effPpJU0AfYHLra92PZDwP3ATkNRRgDbs2z/bajuX0aDMnbVZ7kC2Qm43/aDtpcAF5M+y1CS7ZuAJ/sd3h+4KL++CHjXoBZqGIgKpraPAL/NrzcG/ll1bnY+1o02l3S7pBslvXmoC1PDivBZHp2bSX/QRU0mK8LnVmHgWkm3SjpiqAtTYEPbcwDy3w2GuDwrnWEV7FLSdcBLa5w60fYV+ZoTSStpTqkkq3F9R8d2lylnDXOATW3Pk/Q64FeStrY9v4vKOOif5XIFaFBu4Fzg1FymU4Fvkn5sDLUh/9wGYBfbj0raAJgq6d789BCGoWFVwdjes9F5SYcB7wD28LIJQrOBCVWXbQI82pkSJkXlrJNmMbA4v75V0gPAVkBHOlqbKSND8Fn2V7bckr4PXNXh4pQ15J9bWbYfzX8fl3Q5qXmvWyuYxySNtz1H0njg8aEu0MommsgySfsAxwH72V5YdepK4GBJoyVtDkwE/jIUZWxE0vqVDnNJW5DK+eDQlmo5Xf1Z5i+ZigNIgxW6wV+BiZI2l7QqaaDElUNcpuVIWkPS2Mpr4G10z2dYy5XAYfn1YUC9J+/QpGH1BFPgu8Bo0mM9wJ9tH2l7pqRLgHtITWdH2R6yBR4kHQB8B1gfuFrSDNt7A28BviKpB+gFjrTdv1NzSMvYbZ9lDWdI2o7U/PQw8ImhLU5iu0fS0cA1wEjgB7ZnDnGxatkQuDz/+1kF+Knt3w1tkRJJPwPeCoyTNBs4GTgduETSR4F/AAcNXQlXThEqJoQQQkdEE1kIIYSOiAomhBBCR0QFE0IIoSOiggkhhNARUcGEEELoiKhgwrAh6dmhLkMIw0lUMCGEEDoiKpgw7Cg5U9Ldee2S9+Xjb5V0g6Rf5rWBpijPGgwhDFzM5A/D0YHAdsBrgXHAXyVV4mVtT1qv5lHgj8AuwM1DUcgQVnTxBBOGozcBP7Pda/sx4EZgx3zuL7Zn2+4DZgCbDVEZQ1jhRQUThqNGzV6Lq173Ek/5ITQtKpgwHN0EvE/SSEnrkwKFdk1U5xBWFvHrLAxHlwM7A3eQIid/wfa/Jb1yaIsVwsoloimHEELoiGgiCyGE0BFRwYQQQuiIqGBCCCF0RFQwIYQQOiIqmBBCCB0RFUwIIYSOiAomhBBCR/x/8AptBUkTBHgAAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 432x288 with 2 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"ds.isel(time=0).plot()"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"2. Create cosine of latitude mesh with same NaNs."
]
},
{
"cell_type": "code",
"execution_count": 9,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.collections.QuadMesh at 0x2ae1d1093cc0>"
]
},
"execution_count": 9,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAZUAAAEWCAYAAACufwpNAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3deZhcdZ3v8feHsMkmQRYjiwkaZLsYNaKOV0VBQa8SAVGYEbksl+GOzLjM3BFEBx+XO7iN44xLBhHEEWFQjCCgbHeAcUEIQ1jCGgJqIILAYMBgIN2f+8c5LWWlqut096mu7qrP63nO02c/36pAf/v8VtkmIiKiDuv1OoCIiOgfSSoREVGbJJWIiKhNkkpERNQmSSUiImqTpBIREbVJUokxk7STpCckzeh1LBExtSSpREeS7pO038i27V/a3sz2UC/jGo2kEyQtlrRG0jdaHH+npNslPS7pNklvbzj2MUlPl4lzZNm5xT1eJ8mSPtkhltmS/l3Sakl3NH6XkmZJulDSA+W9Zlf4bH8q6ReSfifp+5K2avpcPy2fdVWne0XULUkl+tUDwCeBM5oPSNoe+BbwQWAL4P8A35a0bcNp/1YmzpFledM9NgC+CPy8QiznADcCzwFOBr4raZvy2DDwI+CQKh9K0h7AvwBHANsBq4GvNJzyKPCPwKlV7hdRtySVGJWkfwV2An5Q/sX+t+Vf3pa0fnnOVZI+Wf6F/ISkH0h6jqSzJa2SdH3jX+CSdpV0uaRHJd0p6Z11x237e7a/DzzS4vAOwGO2f+jCxcDvgBeM4RF/DVwG3DHaSZJ2AV4KnGL7SdvnA7dQJhHbD9r+CnB9xef+GfAD29fYfgL4KHCwpM3L+11h+zyKpBox6ZJUYlS2jwB+Cbyt/Iv9M21OPYzir+ftKX45/ww4E9gKuB04BUDSpsDlwLeBbYHDga+Uf4GvQ9JXJD3WZrl5nB9rMXC7pAMlzSiLvtYAjfd7W5n0lkr6300xPR84Gvh4hWftASy3/XjDvpvK/eOxR3k9ALbvAZ4Cdhnn/SJqtX6vA4i+cWb5Cw5JPwR2t31Fuf0d4BPleW8F7rN9Zrn9n5LOB94BLG2+qe2/AP6izkBtD0n6JkVi25jil/Khtn9XnnIecBrwIPAK4HxJj9k+pzz+T8BHbT8hqdPjNgN+27TvtxTJdzza3W/zcd4volZ5U4m6PNiw/mSL7c3K9ecDr2h846Ao0nnu5IQJZUX5Z4B9gA2B1wGnS5oHYPs22w/YHrL9U4q6k3eU174N2Nz2v7W599KGyv3XAE9Q1Ns02gJ4fN2r17nXaxruNZJwx32/iMmQN5Woos6hrH8FXG37jVVOlrQQeHebw7+wPZ5ipHnANbYXl9vXS/o5sB+wpMX5BkZeSfYF5kv6dbn9bGBI0n+zvaA5nrJOZWdJmzcUgb2Y4i1pVLb/g2eS8Yil5fUj998Z2Ai4q9P9IiZD3lSiigeBdZrUjtNFwC6SjpC0Qbm8XNJurU62fXxTK6zGpW1CkbS+pI2BGcAMSRuPNCygqBR/zcibiaSXAK+hrFORtEDSTBX2Bv4KuKC89qMU9RfzyuVC4GvAUW3iv4siUZ1SxnAQsBdwfkOsG1MkBoCNyu12zqao73lNWT/1ceB7IwmrrCPamOIPxvXKZ24wyv0iapWkElX8PfCRsrjqbyZyo/KX35soKvYfAH4NfJpnfqnW5SMUxW4nUrzpPFnuw/bVwMcomvY+TvEL/v/avqy89jBgGUWR0jeBT9s+ayR+278eWcr7/s72o6PEchgwH/gviqa+77D9m4bjT1IUa0HRmuzJdjeyvRQ4niK5PERRl9JY53REef1XKRLlkxRJL2JSKJN0RUREXfKmEhERtUlSiYiI2iSpREREbZJUIiKiNgPXT2XrrWZ49o5pYRkRnd1w85qHbW/T+czW9n/9pn7k0WqDed9w85pLbR8w3mdNFQOXVGbvuAHXXbpTr8OIiGlgxqy7fzGR6x95dKjy75sZs+7eeiLPmioGLqlEREwWA8MM9zqMSZWkEhHRJcY8PXXnsuuKJJWIiC7Km0pERNTCmKEBG7UkTYojIrpoGFdaqpB0QDlb6jJJJ7Y4PlPSIkk3S7pO0p4Nxz5QTs1wq6RzOgxcOm49SyqSDi0/4LCk+U3HTiq/tDsl7d/m+q3KKWnvLn/OnJzIIyKqMTCEKy2dSJoBfBl4M7A7cLik3ZtO+zCwxPZewHso5gJC0vYUo23Pt70nxejdh9X1ORv18k3lVuBg4JrGneWXdBjFtKkHUEw1O6PF9ScCV9qeC1xZbkdETCk1vqnsDSyzvdz2U8C5wIKmc3an+H2I7TuA2ZK2K4+tDzyrnAJiE4pRwmvXs6Ri+3bbd7Y4tAA41/Ya2/dSDEG+d5vzzirXzwLe3p1IIyLGx8DTdqUF2FrS4obluKbbbU8xyd2IFaw7LfVNFH+sU84F9HxgB9v3A58DfgmsBH7bMNVDraZiRf32wLUN262+OIDtbK8EsL1S0rbtblj+4xwHsNP2U/EjR0Q/csWirdLDtuePclwt9jXf/FTgi5KWALcANwJry+qBBcAc4DHgO5LebftbVYOrqqu/YSVdQeu5x0+2fUGL/VDtixsT26cBpwHMf/HGg9UUIyJ6xzBU32+cFcCODds70FSEZXsV5SykkgTcWy77A/eOTA4n6XvAnwDTK6nY3m8cl3X84koPSppVvqXMopgFLyJiyih61NfmemCupDnA/RR1z3/aeIKkLYHVZZ3LscA1tldJ+iXwSkmbUMwGui+wuL7QnjEVmxRfCBwmaaPyy5sLXNfmvCPL9SN5Zg7xiIgpQgxVXDqxvRY4AbgUuB04z/ZSScdLOr48bTdgqaQ7KFqJva+89ufAd4H/pCgWW4+y9KZuPatgkHQQ8M/ANsDFkpbY3r/8ks4DbgPWAu+1i3EOJJ0OLLS9mKLs8DxJx1BUPh3akw8SEdFGUVHfOWFUvp99CXBJ076FDes/o/hDvNW1pwCn1BZMGz1LKrYXAYvaHPsU8KkW+49tWH+E4hUuImJKKvqp1JdUpoM0hYqI6KLhGt9UpoMklYiILsmbSkRE1MaIoSnZHqp7klQiIrooxV8REVELI55yq6EL+9fAJZU1HuaetU/0OoyIGABF58cUf0VERE1SUR8REbWwxZDzphIRETUZzptKRETUoaioH6xfs4P1aSMiJlEq6iMiolZD6acSERF1SI/6iIio1XBaf0VERB2KASWTVCIiogZGPD1gw7T0LIVKOlTSUknDkuY37H+jpBsk3VL+fEOb6z8m6X5JS8rlLZMXfUREZzYMeb1KS7/o5ZvKrcDBwL807X8YeJvtByTtSTEf8/Zt7vEF25/rYowREROgdH6cLLZvB5DUvP/Ghs2lwMaSNrK9ZhLDi4iYMENfvYVUMdU/7SHAjaMklBMk3SzpDEkz291E0nGSFkta/Oijw92JNCKihSHWq7T0i65+EklXSLq1xbKgwrV7AJ8G/rzNKV8FXgDMA1YCn293L9un2Z5ve/5WW/XPP15ETG1GDLva0i+6Wvxle7/xXCdpB2AR8B7b97S594MN538NuGhcQUZEdImBpzP2V29J2hK4GDjJ9k9GOW+W7ZXl5kEUFf8REVOIBm4+lV42KT5I0grgVcDFki4tD50AvBD4aENz4W3La05vaH78mbLZ8c3A64EPTPZniIgYjSl61FdZ+kUvW38toijiat7/SeCTba45tmH9iPE890mvzy1rthvPpRExcFZ2PqWDQXtTmXLFXxER/cJWX72FVJGkEhHRJUVFfYZpiYiIWqjWYVokHSDpTknLJJ3Y4vhMSYvK/nvXlaOSjBzbUtJ3Jd0h6XZJr6rxg/5BkkpERJcUFfX19FORNAP4MvBmYHfgcEm7N532YWCJ7b2A9wBfbDj2ReBHtncFXgzcPvFPuK4klYiILqqxR/3ewDLby20/BZwLNHck3x24EsD2HcBsSdtJ2gJ4LfD18thTth+r6zM2SlKJiOiSMfao33pkOKlyOa7pdtsDv2rYXsG6g+3eRDFQL5L2Bp4P7ADsDPwGOFPSjWX3jE278JGTVCIiummY9SotwMMjw0mVy2lNt2pVRuam7VOBmZKWAH8J3AispWiU9VLgq7ZfAvwOWKdOpg5p/RUR0SU2PD1c29/uK4AdG7Z3AB744+d5FXAUgIoh4O8tl02AFbZ/Xp76XXqRVCT9U4V7rLL9kZriiYjoG0XxV21J5XpgrqQ5wP3AYcCfNp5QDnO1uqxzORa4pkw0qyT9StKLbN8J7AvcVldgjTq9qSwA/q7DOScCSSoRES3U1aPe9lpJJ1BMXDgDOMP2UknHl8cXArsB35Q0RJE0jmm4xV8CZ0vaEFhO+UZTt05J5Qu2zxrthNHmMYmIGGQjTYpru599CXBJ076FDes/A+a2uXYJML/VsTqNmlRs/2OnG1Q5JyJiMGWYlpYkbQP8L2B24zW2j+5OWBER/SFz1Ld2AfAfwBXAUPfCiYjoH0Xrr8Ea+6tqUtnE9oe6GklERJ8Z6fw4SKoW9l0k6S1djSQiog8No0pLvxg1qUh6XNIq4H0UieVJSasa9o+bpEMlLZU03DCbI5Jml88ZmfVxYZvrt5J0uaS7y59phRYRU0qdA0pOF51af23exWffSjFGzb+0OHaP7Xkdrj8RuNL2qeUQ0CcCKaKLiCll0Fp/Vfq0kq6ssm8sbN9e9uwcrwXASB+as4C3TySeiIi62WKt16u09ItOw7RsDGxKMXrmTJ4Z0GwL4HldjGuOpBuBVcBHbP9Hi3O2s70SwPZKSdu2u1k52udxAFvMehY3rp7dhZAjov8smfAd+qloq4pOrb/+HHg/RQK5gWeSyiqKyWJGJekK4LktDp1s+4I2l60EdrL9iKSXAd+XtEc5fs24lKN9ngbwvD1mNo/qGRHRFXX3qJ8OOtWpfFHSl4AP2/7EWG9ue79xXLMGWFOu3yDpHmAXYHHTqQ9KmlW+pcwCHhrrsyIium3QkkrHgjzbQ8CkNSeWtE05bSaSdqYYx2Z5i1MvBI4s14+k6KAZETFljHGSrr5QtXboMkmHlOPz10LSQZJWAK8CLpZ0aXnotcDNkm6iGPP/eNuPltec3tD8+FTgjZLuBt5YbkdETCmD1k+lao/6D1JU2A9JepKibsW2txjvg20vAha12H8+cH6ba45tWH+EYk6AiIgpyYa19U3SNS1USipd7q8SEdG3+qloq4rK0wlLOpCiaArgKtsXdSekiIj+MIhjf1Ud+v5U4OXA2eWu90n677a7MsdxRES/cJJKS28B5tkeBpB0FnAjxdAoERHRRj9VwldRufgL2BJ4tFx/dhdiiYjoK3bqVNr5e+BGSf9O0fLrtcBJXYsqIqIviKG0/lqX7XMkXUVRryLgQ7Z/3c3AIiL6QepU2lsPeLi8ZhdJu9i+pjthRURMfxn7qw1JnwbeBSwFhsvdBpJUIiLacVGvMkiqvqm8HXhROdhjRERUlNZfrS0HNqAcPTgiIjpzKurbWg0sKWd7/ENisf1XXYkqIqJPpPirtQvLJSIixiCtv1qwfdZoxyWdb/uQekKKiOgPdr1JRdIBwBeBGcDptk9tOj4TOAN4AfB74GjbtzYcn0Ex4eH9tt9aW2ANxtKkeDQ713Sfrls9tAFLHtuh12FExICoq0lxmRC+TDF/1ArgekkX2r6t4bQPA0tsHyRp1/L8xilC3gfcDox72pJO6qpBGnOpoaRDJS2VNNww8RaS/kzSkoZlWNK8Ftd/TNL9DedN2uyUERFV2dWWCvYGltlebvsp4FxgQdM5uwNXFs/1HcBsSdsBSNoB+B/A6TV9tJZ62SzhVuBgmvq62D7b9jzb84AjgPtsL2lzjy+MnGv7ki7HGxExJkYMD69XaQG2lrS4YTmu6XbbA79q2F5R7mt0E8XvVSTtDTwfGCma+Ufgb3mmr2FX1FX8Neb3O9u3A3SYofhw4JxxxhQR0XNjKMZ52Pb8UY63+mXZfPtTgS9KWgLcQjGa/FpJbwUesn2DpH2qhzR2Y35TkTRT0l5Nuz9UUzzN3sXoSeUESTdLOqOsoGpJ0nEj2f+px56sP8qIiFbKivoqSwUrgB0btncAHvijx9mrbB9VlvS8B9gGuBd4NXCgpPsois3eIOlbNXzCdVRKKpKukrSFpK0oXq/OlPQPI8dtX9bmuisk3dpiaS4HbHXtK4DVjS0XmnyVooXDPGAl8Pl297J9mu35tudvuOWzOj06IqI+rrh0dj0wV9IcSRsCh9HU1UPSluUxgGOBa8pEc5LtHWzPLq/7f7bfPdGP1krV4q9n214l6VjgTNunSLq500W295tAbIcxyluK7QdH1iV9Dcj0xhEx5dTVpNj2WkknAJdSNCk+w/ZSSceXxxcCuwHflDQE3AYcU8vDx6BqUllf0izgncDJXYwHAEnrAYdSzNvS7pxZtleWmwdRVPxHREwZBoaH6+unUjZIuqRp38KG9Z8Bczvc4yrgqtqCalK1TuXjFNnxHtvXS9oZuHsiD5Z0kKQVwKuAiyVd2nD4tcAK28ubrjm9ofnxZyTdUr4xvR74wETiiYionQGr2tInqvao/w7wnYbt5cCEetDbXgQsanPsKuCVLfYf27B+xESeHxExGQZt7K+qFfW7SLpS0q3l9l6SPtLd0CIi+kB9FfXTQtXir69RzEn/NIDtmykq0iMioq1qzYn7adDJqhX1m9i+rqmj4touxBMR0V/66C2kiqpJ5WFJL6D8eiS9g6JvSEREtGNwja2/poOqSeW9wGnArpLup+ih2ZWOMxER/SVJZR1la6/9JG0KrGf78e6GFRHRJwas+Ktq66/tJH0d+K7txyXtLmnSe2pGREw7af3V0jcoOj8+r9y+C3h/NwKKiOgb07Tzo6RDq+xrpWpS2dr2eZTj8NteCwxVjjAiYkDVOEnXZDqp4r51VK2o/52k5/BM669XAr+teG1ExOCaRq2/JL0ZeAuwvaR/aji0BRW7kVRNKh+kGGL5BZJ+QjFG/zvGEOuU8fu1G3DXb7bpdRgRMSA09d5CRvMAsBg4ELihYf/jVBxfsWNSKUcM3hh4HfAiivZxd9p+eqzRRkQMlGlWCW/7JuAmSd8e7+/4jknF9rCkz9t+FbB0PA+JiBhMU68SvqLZkv4e2J3ipQIA2zt3urBqRf1lkg5RhwnlIyKiyfRsUnwmxey6aymmFvkm8K9VLqyaVD5IMfT9GkmrJD0uadV4Io2IGCjDFZep5Vm2rwRk+xe2Pwa8ocqFVXvUbz6B4CIiBtNIP5Xp5/dlffrd5RTG9wPbVrmwao/6l7ZYXiCpauuxVvf8rKQ7JN0saZGkLRuOnSRpmaQ7Je3f5vqtJF0u6e7y58zxxhIR0S1ytWWKeT+wCfBXwMuAI4Ajq1xYtfjrK8C1FPOqfK1cPxe4S9Kbxhpt6XJgT9t7UfTQPwlA0u4Uc7XsARwAfEXSjBbXnwhcaXsucGW5HRExtUzDOhXb19t+wvYK20fZPtj2tVWurfqmcR9wjO2l8Idf/P8H+ATwPeCycQTdeM21PNPvZQFwru01wL2SlgF7Az9rusUCYJ9y/SzgKuBDY40jIiIKkn7AKCnO9oGd7lE1qew6klDKG98m6SW2l9fUIOxo4N/K9e0pksyIFeW+ZtvZXlnGs1JS2/I+SccBxwGsv/Wz64g3IqKSKVi0NZrPTfQGVZPKnZK+SlHkBfAuiqKvjSinGG5F0hXAc1scOtn2BeU5J1M0Wzt75LIW50/on8X2aRTzwbDxC7efXv/EETF9mWk1TIvtq5v3lfXVO5bTyHdUNan8T+AvKCpvBPwY+BuKhPL6UQLcb7SbSjoSeCuwr/2HIdVWADs2nLYDxdABzR6UNKt8S5kFPFTto0RETKJp+GespKsohmpZH1gC/EbS1bY/2OnaShX1tp+kqKw/0fbbbX/O9mrbw7afGGfQB1DUgRxoe3XDoQuBwyRtJGkOMBe4rsUtLuSZ1ghHAheMJ46IiG6apq2/nm17FXAwcKbtlwGjviSMqNqk+ECKbPWjcnuepAvHGeyILwGbA5dLWiJpIUBZd3MecFv5vPfaHiqfe7qk+eX1pwJvlHQ38MZyOyJiapmGrb+A9csSoHcCF43pwornnULRAusqANtLJM0ey4Oa2X7hKMc+BXyqxf5jG9YfAfadSAwREV039RJGFR+nmJjxJ7avl7QzcHeVC6smlbW2f5uhvyIiqpuiRVsd2f4OxdBcI9vLgUOqXFu18+Otkv4UmCFprqR/Bn465kgjIgbNsKotFUg6oBxpZJmkdTp8S5pZjlBys6TrJO1Z7t9R0r9Lul3SUknv6/CcHcr7PCTpQUnnS9qhSoxVk8pfUvRwXwOcA6wic9RHRHRUV0V9ObLIl4E3UwxJf3jZEb3Rh4El5Ugl7wG+WO5fC/y17d2AVwLvbXFtozMpGkM9j6Kf4A/KfR1Vbf212vbJtl9ue365/vsq10ZEDLT6Kur3BpbZXm77KYp+gwuaztmdYtgqbN9BMS/KdrZX2v7Pcv/jwO207lQ+YhvbZ9peWy7foJjxt6NR61Tq6LI/1fip9XjqgU17HUZEDIKx1alsLWlxw/ZpZcftEdsDv2rYXgG8oukeN1E0A/6xpL2B51P09Xtw5ISykdVLgJ+PEsvDkt5NUTIFcDjwSJUP0amifqTL/sEUPeO/1fCA+6o8ICJioFVPKg/bnj/K8SqjjZwKfFHSEuAW4EaKoq/iBtJmwPnA+8t+KO0cTdHt4wvlM34KHNXxE9AhqYx02Zf0CduvbTj0A0nXVHlARMQgU30TcHUcbaRMFEcBlDP13lsuSNqAIqGcbft7HZ71CeBI2/9VXrsVxUvG0Z2CrFpRv03ZTpnyAXOoWL4WERG1uB6YK2mOpA0ppgj5o07okrYsjwEcC1xje1WZYL4O3G77Hyo8a6+RhAJg+1GKIrOOqvZT+QBwlaTl5fZsylF/IyJiFDX1U7G9tpyF8VJgBnCG7aWSji+PLwR2A74paYhiVJJjystfTTHR1i1l0RjAh21f0uZx60ma2fSmUilfVJ1O+EeS5gK7lrvuKOc7iYiIdmru/FgmgUua9i1sWP8ZxXiJzdf9mNZ1Mu18HvippO9SpMV30mKUk1Y6tf56aUMztDUULQvanhMREU2mZ4/6b5Yt0d5AkYwOtn1blWs7vamcKWkfRs9wX6diWVtExMCZhkkFiskYKYrQxqRTUnk2cAOjJ5XfjPWhERGDQNTa+mta6NSkePYkxRER0X+m6YCSE1G19VdERIxHkkpERNRmwJJK1c6PtZP0WUl3lEM0L5K0Zbn/jZJukHRL+fMNba7/mKT7y1kjl0h6y+R+goiIzqbpdMLjVnU6YUl6t6S/K7d3Kgcrm4jLgT3LIZrvAk4q9z8MvM32f6OYe/5fR7nHF2zPK5d2nXgiInpnek4nPG5V31S+AryKYiBJgMcpxvUfN9uX2R4Z6OxainFssH2j7ZHxbJYCG0vaaCLPiojoCRetv6os/aJqUnmF7fcCvwcou+5vOPolY3I08MMW+w8Bbhyl9/4JZfHZGZJmtru5pOMkLZa0eOiJ39URb0RENXlTaenpctYxA0jaBuiYWyVdIenWFsuChnNOphia+eyma/cAPg38eZvbfxV4ATAPWEkxrEBLtk8rJxebP2OzzKUSEZNn0OpUqrb++idgEbCtpE8B7wA+0uki2/uNdlzSkcBbgX1tu2H/DuXz3mP7njb3bpx05mvARRU+R0TE5OqjhFFF1QElz5Z0A7AvRSfRt9u+fSIPlnQA8CHgdbZXN+zfErgYOMn2T0a5fpbtleXmQcCtE4knIqJ2fVa0VUWnASW3ath8iGemlkTSVuUY++P1JWAj4PJiqH+utX08cALwQuCjkj5anvsm2w9JOh1YaHsx8BlJ8yj+ye6jfTFZRERPiP4q2qqi05vKDRS/tAXsBPxXub4l8EtgzngfbPuFbfZ/Evhkm2PHNqwfMZ7nrvc0bHJ/z7rnRMSAGbSkMupvV9tzbO9MMSnM22xvbfs5FPUgnaajjIiItP5q6eWNnQtt/xB4XXdCiojoIwOWVKq2/npY0keAb1F8/HcDj3QtqoiIftBnzYWrqPqmcjiwDUUz3+8D2/JM7/qIiGgnbyrrKlt5va/LsURE9J1+GoKlikpJRdK/0yKX2m45gnBERBQGrfirap3K3zSsb0wxJtfaNudGRAT0XdFWFVWLv25o2vUTSVd3IZ6IiP6SpLKupp716wEvA57blYgiIvpEetS319izfi1wL3BMt4KKiOgXGh6srFI1qexm+/eNOzJxVkREBwNYp1K1n8pPW+z7WZ2BRET0o0GbT2XUpCLpuZJeBjxL0kskvbRc9gE2mZQIIyKmsxo7P0o6QNKdkpZJOrHF8ZmSFpUz4l4nac+q19alU/HX/sD/pJg//h8a9j8OfLhLMUVE9I263kLK2Xe/DLwRWAFcL+lC27c1nPZhYIntgyTtWp6/b8VrazFqUrF9FnCWpENsn1/3wyMi+l59RVt7A8tsLweQdC6wAGhMDLsDfw9g+w5JsyVtB+xc4dpadJqk6922vwXMlvTB5uO2/6HFZZVI+izwNuAp4B7gKNuPSZoN3A7cWZ46MnlX8/VbAf8GzKaYpOudtv9rvPFERNTOtQ7Tsj3wq4btFcArms65CTgY+LGkvYHnU5Q0Vbm2Fp0q6jctf24GbN60bDbBZ18O7Gl7L+Au4KSGY/fYnlcu6ySU0onAlbbnAleW2xERU8ZIP5WKFfVbS1rcsBzX4nbNmt+DTgVmSloC/CVwI0U3kCrX1qJT8de/lKtXNM8XL+nVE3mw7csaNq8F3jHGWywA9inXzwKuopjzPiJi6nDl390P254/yvEVwI4N2zsAD/zxo7wKOApAxTzt95bLJp2urUvVJsX/XHHfeB0N/LBhe46kGyVdLek1ba7ZzvZKgPLntu1uLum4kew/tPp39UUdEdFBjU2KrwfmSpojaUPgMODCP3qWtGV5DOBY4Joy0XS8ti6d6lReBfwJsE1TncoWwIxON5d0Ba2HcznZ9gXlOSdTvJ6dXR5bCexk+5GyOfP3Je1RfjHjYvs04DSAzZ6zo4GBlb8AAAwjSURBVLf4xYCNRR0RvVFj50fbayWdQDG9+wzgDNtLJR1fHl8I7AZ8U9IQRSX8MaNdW09kf6xTk+INKepO1qeoRxmxigrFVbb3G+24pCMp5rvf1y7eEW2vAdaU6zdIugfYBVjcdPmDkmbZXilpFvBQp3giIiZbnfOplNO6X9K0b2HD+s+AuVWv7YZOdSpXA1dL+obtX9T5YEkHUNSBvM726ob92wCP2h6StDPFF7S8xS0uBI6kqJg6ErigzvgiIuqQSbpaW102Ad6DYj4VYMKTdH0J2Ai4vKhP+kPT4dcCH5e0FhgCji9nnkTS6cBC24spksl5ko4BfgkcOoFYIiLqZ8ZSUd8XqiaVsyn6hLwVOJ7izeA3E3mw7Re22X8+0LKjpe1jG9YfAfadSAwREd3WT+N6VVG19ddzbH8deNr21baPBl7ZxbgiIvpDjWN/TQdV31SeLn+ulPQ/KNo379CdkCIi+kMm6Wrvk5KeDfw1Rf+ULYD3dy2qiIh+YGeSrlZsX1Su/hZ4PYCkJJWIiE4GK6dUrlNpZZ0BJiMi4o8N2iRdVYu/Wmk1QFlERIwwkOKvygbrm4qIGI8B+03Zaeyvx2n9lQh4VlciiojoI/1UtFVFp2FaNh/teEREjC6tvyIioh591rGxiiSViIguKTo/DlZWSVKJiOimjFIcERF1yZtKRETUI3UqERFRn4z91fdm/H6ILe4c93T3ERFjM2DFXxMZ+2tCJH1W0h2Sbpa0SNKW5f4/k7SkYRmWNK/F9R+TdH/DeW+Z/E8RETEKF9MJV1n6Rc+SCnA5sKftvYC7gJMAbJ9te57tecARwH22l7S5xxdGzrV9yeSEHRExBna1pU/0LKnYvsz22nLzWlpP+nU4cM7kRRURUbMBm/mxl28qjY4Gfthi/7sYPamcUBafnSFpZruTJB0nabGkxU+tXT3RWCMiKtPwcKWlX3Q1qUi6QtKtLZYFDeecDKwFzm669hXAatu3trn9V4EXAPOAlcDn28Vh+zTb823P33D9TSb6sSIiqjFF58cqS5/oausv2/uNdlzSkcBbgX3tdQoVD2OUtxTbDzbc52vARe3OjYjoBeF0fpwskg4APgS8zvbqpmPrAYcCrx3l+lm2V5abBwHt3mgiInpnwJJKL+tUvgRsDlxeNgle2HDstcAK28sbL5B0uqT55eZnJN0i6Wbg9cAHJiXqiIixGLDWXz17U7H9wlGOXQW8ssX+YxvWj+hOZBERNRmpUxkgU6X1V0REX6qz9ZekAyTdKWmZpBNbHH+2pB9IuknSUklHNRz7QLnvVknnSNq4xo/5B0kqERFdU7Hoq0Lxl6QZwJeBNwO7A4dL2r3ptPcCt9l+MbAP8HlJG0raHvgrYL7tPYEZFI2hapekEhHRLabOOpW9gWW2l9t+CjgXWNB0joHNJQnYDHiUossGFNUdz5K0PrAJ8EANn3AdSSoREd1UvZ/K1iOdtMvluKY7bQ/8qmF7Rbmv0ZeA3SgSxi3A+2wP274f+BzwS4p+fb+1fVltn7HBwI1SHBExmcbQT+Vh2/NHOa4W+5pvvj+wBHgDRefwyyX9B0Vx1wJgDvAY8B1J77b9rarBVZU3lYiIbqqv+GsFsGPD9g6sW4R1FPA9F5YB9wK7AvsB99r+je2nge8BfzLhz9ZCkkpERLfYMDRcbensemCupDmSNqSoaL+w6ZxfAvsCSNoOeBGwvNz/SkmblPUt+wK31/Qp/0iKvyIiuqmmjo2210o6AbiUojjrDNtLJR1fHl8IfAL4hqRbKIrLPmT7YeBhSd8F/pOi4v5G4LRaAmuSpBIR0U019pYv5426pGnfwob1B4A3tbn2FOCU2oJpI0klIqJbDGSO+v7m1b9n+Mbbeh1GRAwEgwdrnJaBSyoREZPGVK2E7xtJKhER3dRHIxBXkaQSEdFNSSoREVGP/porpYoklYiIbjFQcVj7ftGzHvWSPiHp5nLWx8skPa/h2EnlfAF3Stq/zfVbSbpc0t3lz5mTF31EREUDNvNjL4dp+aztvWzPAy4C/g6gnB/gMGAP4ADgK+U8As1OBK60PRe4styOiJhCah2mZVroWVKxvaphc1OeGW1zAXCu7TW27wWWUcwj0GwBcFa5fhbw9m7FGhExLgZ7uNLSL3papyLpU8B7gN8Cry93bw9c23BaqzkDALazvRLA9kpJ247ynOOA4wA2ZpMaIo+IqGjAetR39U1F0hXlfMjNywIA2yfb3hE4Gzhh5LIWt5rQv4rt02zPtz1/AzaayK0iIsZmwOpUuvqmYnu/iqd+G7iYYrCzKnMGADwoaVb5ljILeGhCwUZE1M1O66/JImluw+aBwB3l+oXAYZI2kjQHmAtc1+IWFwJHlutHAhd0K9aIiHHLm8qkOVXSiyhmZ/4FMDInwFJJ5wG3UYz7/17bQwCSTgcW2l4MnAqcJ+kYigloDu3BZ4iIGIXx0FCvg5hUPUsqtg8Z5dingE+12H9sw/ojlDOcRURMSRn6PiIiatVHzYWrSFKJiOgSA86bSkRE1MKZpCsiImo0aBX1ch81ZatC0m8oWpu1szXw8CSFMxGJs16Jsz7TIUaoFufzbW8z3gdI+lH5nCoetn3AeJ81VQxcUulE0mLb83sdRyeJs16Jsz7TIUaYPnFON70cpTgiIvpMkkpERNQmSWVdp/U6gIoSZ70SZ32mQ4wwfeKcVlKnEhERtcmbSkRE1CZJJSIiapOkUpL0WUl3SLpZ0iJJWzYcO0nSMkl3Stq/x3EeKmmppGFJ8xv2z5b0pKQl5bJwqsVYHpsy32UzSR+TdH/Dd/iWXsc0QtIB5Xe2TNKJvY6nHUn3Sbql/P4W9zqeEZLOkPSQpFsb9m0l6XJJd5c/Z/Yyxn6RpPKMy4E9be8F3AWcBCBpd+AwYA/gAOArkmb0LEq4FTgYuKbFsXtszyuX4yc5rkYtY5yC32UrX2j4Di/pdTAA5Xf0ZeDNwO7A4eV3OVW9vvz+plIfkG9Q/DfX6ETgSttzgSvL7ZigJJWS7ctsry03r6WYcRJgAXCu7TW27wWWAXv3IkYA27fbvrNXz69ilBin1Hc5jewNLLO93PZTwLkU32VUZPsa4NGm3QuAs8r1s4C3T2pQfSpJpbWjgR+W69sDv2o4tqLcNxXNkXSjpKslvabXwbQwHb7LE8oi0DOmUHHIdPjeRhi4TNINko7rdTAdbGd7JUD5c9sex9MXBmpASUlXAM9tcehk2xeU55xMMePk2SOXtTi/q+2wq8TZwkpgJ9uPSHoZ8H1Je9heNYVinPTvcp0ARokb+CrwiTKmTwCfp/gDo9d6/r2NwattPyBpW+BySXeUbwkxIAYqqdjeb7Tjko4E3grs62c68KwAdmw4bQfgge5EWOgUZ5tr1gBryvUbJN0D7AJ0pbJ0PDHSg++yWdW4JX0NuKjL4VTV8++tKtsPlD8fkrSIouhuqiaVByXNsr1S0izgoV4H1A9S/FWSdADwIeBA26sbDl0IHCZpI0lzgLnAdb2IcTSSthmp9Ja0M0Wcy3sb1Tqm9HdZ/mIZcRBFg4Op4HpgrqQ5kjakaOxwYY9jWoekTSVtPrIOvImp8x22ciFwZLl+JNDuDTvGYKDeVDr4ErARxSs7wLW2j7e9VNJ5wG0UxWLvtd2zCRIkHQT8M7ANcLGkJbb3B14LfFzSWmAION52c8VkT2Ocat9lC5+RNI+iaOk+4M97G07B9lpJJwCXAjOAM2wv7XFYrWwHLCr//1kf+LbtH/U2pIKkc4B9gK0lrQBOAU4FzpN0DPBL4NDeRdg/MkxLRETUJsVfERFRmySViIioTZJKRETUJkklIiJqk6QSERG1SVKJgSHpiV7HENHvklQiIqI2SSoxcFT4rKRby7k/3lXu30fSVZK+W86tc7bKnnwRUU161McgOhiYB7wY2Bq4XtLI+FQvoZjv5QHgJ8CrgR/3IsiI6ShvKjGI/jtwju0h2w8CVwMvL49dZ3uF7WFgCTC7RzFGTEtJKjGIRivSWtOwPkTe5iPGJEklBtE1wLskzZC0DcVgnFNmtOSI6Sx/hcUgWgS8CriJYkTiv7X9a0m79jasiOkvoxRHRERtUvwVERG1SVKJiIjaJKlERERtklQiIqI2SSoREVGbJJWIiKhNkkpERNTm/wPsVpJEtvKIlAAAAABJRU5ErkJggg==\n",
"text/plain": [
"<Figure size 432x288 with 2 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"coslat, _ = xr.broadcast(ds.coslat, ds.isel(time=0))\n",
"coslat = coslat.where(ds.isel(time=0).notnull())\n",
"coslat.plot()"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"3. Area-weight"
]
},
{
"cell_type": "code",
"execution_count": 10,
"metadata": {},
"outputs": [],
"source": [
"area_weighted = (ds * coslat).sum(['lat', 'lon']) / coslat.sum(['lat', 'lon'])"
]
},
{
"cell_type": "code",
"execution_count": 11,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"[<matplotlib.lines.Line2D at 0x2ae1d110d860>]"
]
},
"execution_count": 11,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXAAAAESCAYAAADkJY5uAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO2dd5wV5fWHn7MVlg4ugsC6AoIFVIoK9i5KErsxMcYkJsb6s0aJJSaxEU2MGo2KPZZoFNREREVFBQtKk65UKdKR3ra8vz9m5u7euW3m7p177+ye5/OBvXfquVO+c+a85z2vGGNQFEVRwkdBrg1QFEVR0kMFXFEUJaSogCuKooQUFXBFUZSQogKuKIoSUlTAFUVRQkpRNne22267mcrKymzuUlEUJfRMnjx5rTGm3D09qwJeWVnJpEmTsrlLRVGU0CMi38abriEURVGUkKICriiKElJUwBVFUUKKCriiKEpIUQFXFEUJKSrgiqIoIUUFXFEUJaSogCuKonjgwffncd6Iz3JtRhQq4IqixPDl4vVUDhvNyo07cm1K3nDf2G/4fOH6XJsRhQq4oigxPPvpYgC+WJxfgqVEowKuKEoM1TXWUIslhZJjS5RkqIArihJDVU0tAEUFKhH5jJ4dRVFi2OUIeB574AvXbOGRDxdkfb+Tv/2eaUs3ZH2/8VABVxQlBieEUlyYvxLxk8c/5y9vz2XTjirf6w4bOZ1Bd72f1n7PeuRTTn/4k6TLfLV0A7e+PhNjTFr78Er+nh1FUXJGjS08BZIdD3xXdS0rNm73tc62XTUAmFr/+3vpy6Ws3BRchs2PR3zGc59/y/aqmsD2ASrgiqLEwZFtQ+Y8yEc+XMCHX6+OO2/YqOkMvvsDtu/yLnhB2Jgpamqz8wDM6oAOiqKEA6lTx4zxl7fnArB4+NCYee/NXgVYnnjzkkJP2xPbyICjFGlRbQt4ba5DKCLSTUTGicgcEZklIlfZ0/8oIstFZJr979RALVUUJWsIwYdOtuysbtD6BbaJ+aDfC9ds4Qf/GM/81VuAuodKbcDGeQmhVAPXGWP2BQYBl4vIfva8vxtjDrL/vRWYlYqi5ISg9Gf2d5voc9s7vDFtedrbcDzwoL1cL4ycsoyZyzfx/pxVUdNz3ohpjFlhjJlif94MzAG6BGqVoihxeXfWSiqHjY54ekHhhFCC0p8Zy600vAnz1kZNf3PGd9w5eranbTjvCEEJ+I6qGv7z5VJfIryrOrpFNR888AgiUgn0Aybak64Qkeki8pSItMuwbYqiuBg9YwUAM5dvDHQ/QSef7HLSFIuiJejm12by+PhFzFu1mXmrNkfNq601PDF+Idt2VUfZmAn9fuaTRXy+cF3UtOFj5nLDyOl8+PWalOs7Iadl32/n318siUzPuQfuICItgZHA1caYTcAjQA/gIGAF8LcE610sIpNEZNKaNakPhKIoiXE8ukwI7LuzVvLOrJVR03ZW1/D7UdNZt2UXkLkMjxqXK1ple6olCfLMT/z7x5z494+jpr0zayV3jJ7DPW9/bU+RuNtOhz/+bzbnjfg8atqazTsBb7F653y8PGkpvx81IzI9LzxwESnGEu8XjDGjAIwxq4wxNcaYWuBx4JB46xpjRhhjBhpjBpaXl2fKbkVpktRmMD/74ucm89vnJkdNe/OrFfz7i6XMXWl5v5lyIJ2u+Q5OT89iHz09t9ophpu2Wx13nEMQWAzcxyFOtGjOPXCxWgqeBOYYY+6rN71zvcXOAGZm3jxFyT0ffr2aD+auSr1gFqgNOL/YCU84ZEp+3CLrfC0o8P47ImIoUX+orQ1eKJPx6fy1bNoR30vPBw/8cOAC4DhXyuA9IjJDRKYDxwLXBGmoogTN0vXbqBw2mrGzo8X6F09/ya+emZQjq6JxhDCoHu7unoOZ8m7dQhbxnmsNM5dv9PSgqL/Mt+u2RrZx5UtT+fFjn8cs/9zn3/LdBn+9Oy9/cQqTv/3e8/Ibt1Xx0ycm8oxdftdN0A+WlB15jDETiP+GoGmDSqPCKVD0+tTlnLjf7jm2Jj51MfBgPPBqt9JmSH/cQuY43o+PX8Tj4xclXG/YyOkcXNmeswZ0jdgyaspyRk2pSz/8ylVYqtctYzimVznvzl7Fvzq2ZOy1R7N47VZuHDmdP522P18uSlzjfPT0FUxcuJ5Jt5zg6XftrE7eczTo9wLtSq/kjI3b/RchCpLIzZa/BfgiQuglhLJgzRYqh42O9HL0grsDTyY88F8/O4lnPlmcdD+JeOnLpVz3yleA9wbVXdW1vGv/5u+3WY2xf333ayYuWs+Q+8dz6xuzIsvG85DdvVDnrd7CfWO/ibuvVBblvCemogTBf7/6jgP/9C4zlnlPh/t23VYeHjc/sNdSP+KYKxwHeeP2qphOI26+sD1Nd0goHmu37KTPbe/EpCdmIob73pxV/M0lgOkc4iBOe7wMlgKxjp0jvg++P48H35+Xlk1bd9awfuuuBtuZCBVwJSdMmGellM76zruAX/bCFO5952uWrvcX1/RKpGEtj/R7844qhtz/MVOXWHFZR3Cuf+UrLnp2Ehu3JX6L2Wqnv7UoTV3y6LMF69iyszqSZ+4Q1MMynRBQOpZs3VnDL5/+ImEsPN4DatWmnZz72GeMmbkydmY9Ln1+Ms9//m3SZYY+OJ7+t4/1bK9fVMCVnOB4uX5uykJbWddt3RmARXWvu0Ho93OfLebTBWtTLudm7srNzF25mT/+z+qd6D5eVbWJa6lGUvWKUv+ikqL4UhBUFkU6xzidZ8n2qhrGfb2GKUviD8DQkBDHmJkreWjc/KTL7KxOo9atD1TAlZyQTg5vqS0yVTXe1xk/bw2Vw0azaO3WlMtGstQCCKHc+sYsfvr4xNQLunB+s5M+6MeySJzZw+EqTSDgQXng+fKWkw91VBqCCriSI/zfwekUL3pj2ncAfOlhdPWIB54n4gJ1IuzXtsc+WsC366yHlpejlcgDn/zt91z78jSWrNuW8hjW1pqYjJBEpPOQTLWKMSbyoPNKyPU7HAK+o6qGxz5aQHVNcK8j6WZE/GfSUiqHjY7EGxV/+LmBgi5e5Gw1nxsxE1n29swVjLHj15t3VHH3mLm89OVSwJsXnWjotCcmLGLU1OUcde84znn0s6TbGDF+Iac9/AkTXTVF4pGOB55qnZpakzSkFA/1wLPAA+/P4+4xc3ltavqlJ5PxzqyVHPind5n8bWovzc1TE6wc1sXrUr+iK3WI97f7CAUZKOBfW2tiunU7mDRj4PNXbw6wsc/6myi840y/5PkpXPrCFNZu2RkZasy9zFdLNzBqyrL4+/Fp1+YdVTE50HNXbAJguZfOM+l44CmsrK41vuuihFu+QyLgTlGZoB6WThWyaUv9V3hr3bwYyL+c5nwncismOKmvT10eybyIrGOv9N2G7TyXovU/Eb969kv2vnlM3Hl1Ihl/3fvf+4a/u9LhZizbyAn3fcwTSTqj1Bf3XdW1gfbOG3jHe/zooQnR+7f/nvbwJ1z7n6/irudXT/v+8V0uePKLqGkFkRBX6vXTecdJZWNNrfGtESF3wMMh4E733mYeh1ryS6E4Vc38h2jcDWu3vTGTymGjM2dcBpgwby2nPjA+plZxLknlgV/98jTO+OenvD51eWSwW0cgfvfqdG59fSarN/sflDZZadBU5Unvf28eD7jygR1vM1l8uL5X2OuWMQwfYw0tNn/1Zk5/+BOWrt/GR9/Et8uxxdmCF+FbtSk6S8ebSHmXVCfO/IWrR2PdEGe5ST00+A+J5LKGSiYIhYDvSlF6Ml2MMTzy4YKI9+w1xP7dhu3ssB8qBa6GtWc/S88zDJJho6Yze8UmVm4MbhRuvzivw+77Z/Dd73Pjq9Mj369+eRo/sct8uu/fTJQRjbIpjdTGEjtFL1FY5s3p3zHKFfpzYtN3jp7DtKUbOPKecVz41Bdxc7pT9T700jvRi6j58cATpcZ5rc/d57Z3ImUL/JDKxFrjv/httvQ7qAdFqAY1znTb0oI1WyMDrYI3D9wYw2HDP+Ck/XZnxM8H1o3Ll8dP8jpvN7s2vj51Obtqajl3YLfENrmO24qNO3h50tKoaSs37bDXib4AMi7g9l8/XpzT+JcotfGKF6fGTCuyLxp3neldcR4CEQ88URaKB1MzfW0muo4KPF5nW3ZW8+rk+LH4hpDOz8zWHVFrwEflXM+EwgMPCnerthcP3BENp9ZCJLXNte62XdUpC8E/+tECKoeNbvDgrqnIxgC18bj65WncUM+bro8fixL1kEwj4hXh9anLeWXSUmYs2xjppSdptKwWFVi30IT5a6kcNpol67alXse+k2MbJGN3HKdSR9Q3L88wT5X+fPzmRPusS3n0vq2MYvCtyNnKQsm0s+EQKg88nWM9asoyrv3PV7x7zVEUFQjdy1sC8Mf/zoopJl/jYQfuam2OqLgvhP63j2VHVS2Lhw9NuC2nG+76Lbto6aG7c0PJx5cET+Ji/3Wn9zXk5rv65WmRz3vt1oJx1x8TkUaD9QAuKSxg/bbkdSyKXNfQ5CXrqehQlnwdW/TdD6R4vyaV9+zlGHg7TN6PZb6+bRrSa8TMxu+prq2lJAB/ORQC3hD/8aEPrK6uJ9nDMzmCGq9+rxNCGfHxAu56ay4z/ngSrZoVu5aJPtmJ4qY7qry7hxu3V7F+6y7atyjxvI4f0knZCxpJIyXQfR1kSrycXpq2rjJ6+gpem7qc/hVtE3bBjtjguh4ccU5GxAP3UPnPmZIoQ8bbAzB6qcc+WhDTaOrnPCRatG60IO/byiTpOLkGk5U3hqCeEaEQ8IbgZ8QP50Q6KWrrt+6KEfBEHnhDnuI/tNO+nIfLio3bKSspok3z4mSreSbiWeaR55ROV3q3eL07exWrNn3LbT/cPzM22UfKiUUnEu+Rk5exV3kL+le0i3lrK/JwvTk1Xdxan0xIHBGOeYh5UB/3Ib57zNzYZVJuJfH23NsIqlZ5KhvTub4f+mA+u7dulp5BPggqVBMyAU9+EEZOXsYhe7WnW/u6V1g/3oBzMzix1Xg98tw3jJ/cVzfuzf/8qS/YvKOKqUs20KFFCZNvPdH/RuPuJ396FhpjWLN5p7+4fIJOLE463s8HV9KitJCOrWJvRD8/3euyTn3qxcOHxsThizxkSjnXjPv6ihsDN9F/0yHVtfnwuPn+ehLX295bM1bQa/dWvDZ1WcJMnEyRSqAN/o/Tv7KUNRaU6xQKAfdyYxljuO6Vr2hbVsypfTtzxbE92aNtc1/doh3vOvIqGEf9Yz1w//U5EvFxvdfadQHUEE5k4dcrN9OtfXPKSoK/HN6Y9h1XvzyN/hVtLZt8HLZED+Nj//ohQEx7w9L12wKP+6fjgScivq0m8n+8utKvTFoak/edaBuJuPedr5POj91a3fYue2GK392lTWoPPPuZVl4xAT3bGk0WiqOrG7ZV8eLEJdw40sp+KPRxQ9W4BTzOqrEx8Oj95wPO2I5ObYy6EErsslU1tZx8/8dc+nycGzEAnN6VM5Z77/VaFz7wfi5nfbeRI+8Zx8gEXcczhd/iSfVxOxfxMhWcczZ/9Rb63z6W6a4BMB78YH5M2mWsjWmbGBevD8U1m3fyzarNGdxxitmWguclQT1YQuGBOyS7cBK9XrkFPNlrmCPczn0Uz3t3e1wRDzyNGzkT6X0ffr2aHVW1DOnTKTLtq2VW7PZ/07/jlL6dk7YCO50yvFTrywROb1onb9rLhZ2qi3s8lq5Pnc6XaD9+cIuut4bV+A5CvCwo9xQnJ94PmRaPVFtb9v02nvlkEcPfnuurMT8Z3X8/mv4V7VLalUfNPFEE5eCFQsC9CF3C3NSY1LPEehYJodh/410MNTVuAXe26+0Mrdq0g46tSn3Hpatranlr5kq6tWvOrO82UVwofLpgXaRcav3wwU77piktcpcesGz86Js1dGrdjN6dWkUqPBYVCBu3VdG6eVGDY+aTv/2e2XZhIzdlxdGXnJfDls61n05FwXSEzi26/kJCrmsziQfeEGpq4bWpmXsTSRWLfvCD5IMcpEOtgUkpRovP6xBKrhoxRaQb8C+gE1ALjDDGPFBv/vXAvUC5Mcb/kCM+SHYIEp04d++n6trahKleta4QirPNbbuqKSwQSosKE3rgXs7P/NWbOeG+j7ll6L78+sjuqVeox2MfL/Qcq3S8aqdOizuEcuFTVhGixcOHRjIuNu2o5sA/v8tdZ/Tlp4dW+LLNzVmPfJpwnjtvOqjbzU/ozCGtnnxuAbf/vjp5WaSGT8w69t+Y0gBxGzEbfoRGTlmWsVCSSaO7erZIJw88W+TSA68GrjPGTBGRVsBkERlrjJlti/uJwJJgzLPw1ogZf7r7Rv7fVysSFnWau3Izh971Ht/bNSmcg77fH96h9+6teOeao2JuKMdbve6VryLZCYlYvNZ6rf9swTp+fWR3XyGB1T5eneuK/0f3+It3iNzHYtzXqxss4MlIyzNO4650n/fHP16YsFdo3X587yZmnbGzVzJ1yff888MFKdf1Uhog3/So1jS9MEUmyJkHboxZAaywP28WkTlAF2A28HfgBuCNQKyLscX/PPdNcn0SkXUX2Kl/0L+2G2PcF4kfPXI8LD+56Q5+1nHEK8Y7jNuImX4IwC//eH8eS79PIzbtcbnKYaMjNWrcAj5v9ZaM7SfZOv+Z5N3TdZ/SER8vZP3WXTx30aF1288zUXp18tK8LZ1sTP564HmRRigilUA/YKKI/AhYboz5KlnMVEQuBi4GqKhIz7PzUowpiET5ZGld6eCEaArT8EL9eK7uokzJijTFFvBK/vsufOoLDtmrPZcf29OzPRc8OZGDurXlH3Fio3NXbOLW12fy59P2Txp7n7Fso6dz7NSoSecYp+MlpXXZOY2yrtYYpz0jetH8UqQbR87ItQkJyVfxhjzoyCMiLYGRwNVYYZWbgZNSrWeMGQGMABg4cGCDfkVSD9zHsg3ZXyY6VHjobR2DH6fd3bnI0bJ4r+juSal+30ffrOGjb9b4EvDx89Yyfl785pHXbdG64rieCXvEGVPXW9UracXAfa/RMIH19IzJY1HKR/JVxHPalV5EirHE+wVjzCgR6QvsBTjed1dgiogcYoxZGYypia/lVZt2xHSwycz+Usck/ciEE0KJxKZ9rOvHA3fCLZu2V3H3mDlUOyl7Hh5Iubr+M32B+xXw296Y6amYmZuGlDDNow6yjYLaPA6h5MwDF0ttngTmGGPuAzDGzAA61ltmMTAwqCyUuuL/8Q/CoXe9n5bHlYp4z4SGnAgnhLJ0/TZueX2GrxxZrzHw+l2inXBCZP9x84zdMfDc3AGZDhX4Fcd0B+II4mg5DdZlpUXsDLh7emNi+fcexuLMEbn0wA8HLgBmiIhTg/MmY8xbwZiUmGTi6Q4PZEIQktWmSAfH/unLNsb0qEuFF/3eVV3L/re9w24t41c1jCvgHn/P8DFzI6MQBUG+ek6pSC9unnydQ+56H4C+Xdqol+6Dnz4xMdcmJCRnAm6MmUCKt31jTGWmDEq+n2zspY54HrjbhiCKJcXDS6PcKttzW7slfh0VL78n0SF+9KPUaXENYUdVDRtS1N72Q74/EFZv2pHSRmcwbyX85LwRM6ckaYQLltQhBz/d4dPJgY7sx8O6icYqrCP178lVLu25j32W8MGTDtn6Gencl4vXbYt42clIZ4xHJT8J6jyGqphV1uU7w1koDcGb+Cc3zpMHnqMfmEnxhoYVmfJDkBJryN7vUIKlaXvgDlm+ljfvrGbTjrpOC1U1tTHZLv5CKA3xwFMvk+oayacHUtBk42c9+tECtgU4nqkxJq3MGCX/yGkaYc4xzp/sXsxn/jO6psfeN4+JGUfTDw1pj/Kybipnze1d3/XWnJgi/OlcaFU1tbw7axWn9u2UN4NHZGOw2uFxRrbJJMbkImyoBEFQb7bhCqH4OAZB3b/urud+aEgM3AupHnDuuSM+XsjTnyxu8H7/8f48Ln9xCh/MXd3gbWWMRqB7+ZzXrPgjqOdwuAQ81wbEIU8cTiD1Q8uLV5rO7/lmlVVnZNuu4NIM/ZKP14pfatUDbzQ06QEdnB+fb97IkxMWsWmH9xhoQ17rMxEDT+ca2lldw9adyYW52q6nUlKUH/7AE+MXRkrphhmNgTceMj0qkkMoBNzBEfL356xiy85qTjuoS+Jls3Dd3/7mbF/LB21SqgeEl/27Y9inPfQJc1cmHxar2B7Id8uOamZ956+DUhDcMXpOrk3ICMZoFkpjoUl74A4T5q3l2U8XR17Zkwl4PpKrFD0HTyEU++/Y2av40/9mscxD92RHwFPVQ1f8YchOY6wSPE07C8VmzMzoOlmVw0YzqHv7HFnjn6DvxXTSCN04Dvgf3pjJio3eBpFwBFzJLLXG5PUgBYp3grr3Q3HnJfvxny/MzmC8mSBdb+rTBWs9Zb/4zUKJh+OB+8mYUf0OBnW+Gw9NuiPP4B4dYrzvVNQaw4dfr86rdAS/5/CVSUtpUVrEZS9M8bS83zzwePjN49Z6HcGh4ZPGQ5MW8KF9O/OHN2b5WmfiovVMXJRf3rnfk/i7FGM4ukkl0F52v3bLTq56aSo7q1OnBH6zajMn/f1jr+YpPlH9bjwEdSpDIeD50rsvXaprahk+Zi7tWsQv85opUmcRpr6M/JS6XeBhnEklfbSUVeMhZ4Ma5wPhlm/rbeCJCYsC308mGjF97S+zm1NcaANm46FJN2KG3AGnrKQw8H1MWfI9367bmnSZTAuCvuIHS67TTpXMEdTDOCQeeHgV3BhDWUnwh9ldeCuRLZlEG9mCRT3wxkNQ90ooPPAQ63de3YR5ZIqiNCk0hBJSjDF50xilDrOi5IYmXU42xPrNdxt28P3WqtQLZoHq2lp2pRx2zTtFXkZaVhQld2mEItIN+BfQCagFRhhjHhCR24HT7GmrgV8YY74LwsgwpxEede+4XJsQ4YoXp1KSwW6TYT4vipJNchkDrwauM8bsCwwCLheR/YB7jTEHGGMOAt4E/hCIhYA6epljV00m61pqTEZRvJCzGLgxZoUxZor9eTMwB+hijNlUb7EWBHg3hzkLpTGjMXVF8UZedKUXkUqgHzDR/n4n8HNgI3Bshm2rt9+gtqw0BNVvRfFGzrNQRKQlMBK42vG+jTE3G2O6AS8AVyRY72IRmSQik9asWZMJm5U8QT1wRfFGUJlongRcRIqxxPsFY8yoOIu8CJwVb11jzAhjzEBjzMDy8vK0jFQPPD/RjjyK4o2ghlRLKeBipRo8CcwxxtxXb/re9Rb7ETA38+bZ+9IYeF6i8q0o3shlNcLDgQuAGSIyzZ52E3CRiPTGSiP8FrgkGBPVA89HLv7XJNo0L861GYoSCnLWiGmMmUD8vjRvZd6c+Kh+5x/vzl6VaxMUJTTkvBEzl2iHEUVRwox2pVcURQkpQRW1C4eAq4IrihJicppGmGs0hKIoSphp0h64oihKmGnSMXDQMIqiKOGlSWehgDZkKooSXpp0DBw0Dq4oSnjJWVf6fEHlW1GUsNK0BzVGY+CKooSXoGqhhEfA1QdXFCWkNPksFNVvRVHCimah5NoARVGUNGnyHXk0Bq4oSljRRkz1wRVFCSlNvhGzQPVbUZSQ0uQbMbUjj6IoYUUbMXNtgKIoSpo0+Ri4KriiKGFFPfBcG6AoipImTd4D1xi4oihhJWceuIh0E5FxIjJHRGaJyFX29HtFZK6ITBeR10SkbTAmOnYEuXVFUZTgyGU52WrgOmPMvsAg4HIR2Q8YC/QxxhwAfAP8PhALbVS/FUUJKznriWmMWWGMmWJ/3gzMAboYY941xlTbi30OdA3GRAsNoSiKElbyohFTRCqBfsBE16xfAWMSrHOxiEwSkUlr1qxJx0ZrO2mvqSiKklv+8vZcxs9LX/8S4VnARaQlMBK42hizqd70m7HCLC/EW88YM8IYM9AYM7C8vDxtQ9UBVxQlzAQRRinyspCIFGOJ9wvGmFH1pl8I/AA43gTVV7TOimA3ryiKEiBBKFhKARcr+PwkMMcYc1+96UOAG4GjjTHbArDNZUfQe1AURQmOggBEzIsHfjhwATBDRKbZ024CHgRKgbF2A+PnxphLMm6hjeq3oihhJggnNKWAG2MmEF8/38q8OYlRD1xRlDAThISFpyem+uCKooSYIFKhwyPgqt+KooSYIDQsPAKeawMURVEaQNMOoagLrihKiGnSIRRFUZQwE8SwkKERcHXAFUUJM007Bq4CrihKqGnCIRRNI1QUJcw06RBKED9eURQlWzTpRkzNQlEUJcw07TTCXBugKIrSAIIoZhUaAVcFVxQlzDTtLJRcG6AoipJnhEfANQauKEqIUQ9cURQlpDTpGLg64IqihJkm7oGrgiuKEl6C0LDwCLjqt6IoIaZJ98RUFEUJM007hKIuuKIooSYHIRQR6SYi40RkjojMEpGr7Onn2N9rRWRgxi1z2xH0DhRFUQIkJ6PSA9XAdcaYKSLSCpgsImOBmcCZwGOZNysWdcAVRQkzQaQRphRwY8wKYIX9ebOIzAG6GGPGQvZCGyrgiqKEmZwXsxKRSqAfMDEAW5LvW4MoiqKEmJw2YopIS2AkcLUxZpOP9S4WkUkiMmnNmjXp2GhvJ+1VFUVRck7O8sBFpBhLvF8wxozyswNjzAhjzEBjzMDy8vJ0bLRsSHtNRVGU3FMQQM6flywUAZ4E5hhj7su8CR5RF1xRlBCTk0ZM4HDgAmCGiEyzp90ElAL/AMqB0SIyzRhzcsYttFH5VhQlzBQG0BXTSxbKBBLr52uZNScx6oArihJmmnZPzFwboCiK0gCaeDlZlXBFUcJLYVMW8CAqeSmKomSLpu2BaxBFUZQQk5M0wrxB9VtRlBDTxD1wRVGU8NK0BVwVXFFCz3MXHcKffrR/rs3ICUGEULx05MkLNAauKOHnyL3L2b6rJtdm5AT1wBVFCT1NNSW4SacRNtFzriiNjqZ6KzfxnphN9bQrSm4Y1L09pUWZl4ggYsFhIIg3j9AcSvXAFSU7lBYV0KVtc/56zoF8fccpNCvOrEyoM5Y5QtOIaUyuLVCUpkG39mW8d+3RgW1fnbHMERoPXFGU7BC0vmYilBBEadYwEhoBN6gLrhwh17oAACAASURBVCjZIGgPuSGbr+xQxqM/68/eHVs22I4HzjuIB847qMHbySWhEXBFUbKDO0ZdW5vZ7TfEFWvZrIghfTpnxI7j9unIaQd18bTsCft25Jje/oaEvPSYHpw7sGs6pnkmNAKuMXBFyQ5uD7wmwzefycD2qmsbvo0SHxk2PzxwD5755SG+tn9s747cc/aBfs3yhQq4oihJqc3wzdeQHonO20GNDwF/9Gf9KS0qoH2LkqjpJYXByl82GmtDI+CKomQHdyNjpp2nogw0QFZ7jOssHj6UIX068/Udp1DZoSxqnp/G1GbFhb7sg+x0WAqNgGsjpqJkh6CFJ1EGiZ+GyZqa+HpwYLe2DD+zb9x5jmD/7uTe3DJ036Tb79OlddT35ukIeD544CLSTUTGicgcEZklIlfZ09uLyFgRmWf/bRekoRpCUZTsELTwFBVG72CPNs0AePm3g/nmjlPoXt4i5TYSxcDfuPxwzjukIu48Z68HdG3Dr4/snnDbC+46lZ8dumfUtBalybvMPHbBAGb96WTevPIIDurW1tpfFhTciwdeDVxnjNkXGARcLiL7AcOA940xewPv298V4OwBXfnFYZW5NkNR0qIhulPeqjTlMmUl0WJ4Zv+uLB4+lPYtSigpKuDtq45i7u1DktpWnEb82om9u7NsXvzNobx55RGR7/XfEI7uVc41J/Sif0XbpNs+ef9OtCgtok+XNpFYQV6EUIwxK4wxU+zPm4E5QBfgNOBZe7FngdODMhIalnqUbc7s34Wrjt8712YoGWT8Dcfy1W0n5dqMrJBOV/cubZuzePhQeu2eOgyyb+fW3Hv2Abx37dH0q2jLBYOjvd2SooKUMee+Xdr4ttERf/cD6rAeu9HH3t6Fti1VNVaMvWu75lx1wt7+vGk7XOCsM/LSwTz1i4G+7fWCr670IlIJ9AMmArsbY1aAJfIi0jHj1tUnRApeIBJI7V8ld7RvUZLyNbqx4L50/3l+f16ZtJRxX69JuI6TGtihRXwPvGOrUnp3ahX5fs7AbgC8dtnhadn4t3MP5LRv9uDSF6Z4XseRkET35uLhQyOfd9kxdneq4XvXHsXmHdWc8c9PU+7H2cuAPdt7ttEvnq9IEWkJjASuNsZs8vpEEpGLgYsBKirix6YaGwUiSI6bh1uVFrG9qoa2ZcWs3bIrt8bkKRNvOp6tO6s57m8fpVw2nSyEsOK+s0/t25lT+3amctjohOs4onX76X04cu/d+N2r06Pmf3HzCRm1rUVpEaf09dmhxzbSi3Q5Hrg71bBnx1bxFo/ejY/9NBRPMiMixVji/YIxZpQ9eZWIdLbndwZWx1vXGDPCGDPQGDOwvNxfT6agOX6fjhzY1f+rWCoKC4IZfcMPXdo1Z/5dp9Jr99QXXFOlWXFhys4c++9hZSM0qdobaVy7jmi1aV4c8a6zyZF770a39s0j3y85ugcv/ubQqGWcTDYv92ZVtSXgiWLt71x9FB9ef0zceZcf2xOAvXZL3RjbULxkoQjwJDDHGHNfvVn/BS60P18IvJF58+oIIo3wyV8czJ9P6wNADw8t314REXJ9vzsdHTR7JzEFkjpT4NJjekS9WjcWWjUr4qz+mevm7b4///jD/fhJgmwQv7x++eHc/+PkNUtGXDCQ8TccF/k+7JR9OKzHblHL1PrwjJ1wzwEJHLzenVpRmUCgh/TpxOLhQ2nVrDj1jhqIlxDK4cAFwAwRmWZPuwkYDvxHRC4ClgDnBGOiRVA1hA/o2ob/XXEEC9du4aqXpqVeIQ77dGrF3JWbI9/zIQae6d5zjZECib2qmhUXsKOqlq7tmrPs++24s9VGXDCArbuqueblr7JmZxC8dtlh9OzYipFTlsXMK07D+3Bfbr84fC8A/v3FkrTsA/jbOQfSvmUJB3Vry+YdVUmXzfQb0kn7d+Kj3x3Dnh2SO3a3/mA/Dt0ruBh3KrxkoUwwxogx5gBjzEH2v7eMMeuMMccbY/a2/64P1NKA9FBE6Nu1TYME96he5SwePpRiO7/V8uwyZWF6ZKBUREZ46eJBPHJ+/1ybERcR65W/Pj07tmTx8KGRXF533Y6T9u/E4S7PLtc4YR5/WBfojUP24f+O6xk1x52nnSvOGtCVY3tbuRGRR22CG6vYg83OufT661KJN8BFR+wVyWDJBaHpiRl0SCKTgmt5dpk3+NS+nRjUPTNP+34Vbemdhfj4Qd3aRmUf5BMFIrQoLeLNK4/g7z+OLjq0y46BxhtSLJPentOJpSE8cF4/pv3hRF/rONf7pcf04OwB0THrdH5f0P6CMwxbIp32klQRyQ7Jj+dTRgiRgAd71BtWYMf1XYJ54PQsb8lLFw/2tGyqim8v/noQ//xZ8J6xdSzy845xzOrTpQ3tXelvOyMCHpt9UpSBQR07t2nGrD+dTNd2ZakXTkHr5kW0Lasr1NS/oi1n9k9eJrX+GXGfnuoE3dSTEXTEzmlMdO+mu4+GwrrskPy8HtMhNALu5Zg3pEhOg06pvbJzgeRDDDwVIlnqKSaSEY91wJ7tmHDjsZzat1PCZVqU+Ev1q3+OJPLX+tSitND+G9tMlIlBeWuNoUVpEfeecwBn9uvCpFtOiOoN6If2ZdFV9kZddjjnH2o1IFa0L2PSLSdEQkIO9UUsRsBTxN8uHLwnvz9lH9fUYBU80f306qWH8cbl3nLJnS76rZs1nnz+0PySVII49/YhTFy0nguf+iKt7XvV25cvHkRxUQFnJknkLxDJ+Wuac4Mmyt7JlmecqX0UCHRtVxb5Xecd3A1j4OVJSwF488ojaN+ihMOGf5CWbc5nJ63wztP70r+iHQdXxpb4SeSBD9izHft2bsXzn8dvuLvnrAMoKSrg6penYacZs2eHFtxnZ1js1jJ1N3QgUqbhwsMqmb5sA0VxUt2c31NrDLu1LOWO0/twx+jZfL7QaqqSOMs6pIon/8nO3Lp7zNzItGy1mbv3075FSUyZ2ETceXpfTjuoi6dc7rAQGgFP9drTrLgwrdZzr9svKSqga9vmHNq9A6s37Yi7TEGBQK3xlJ4WNKlCKEJ2HjIC7Kpp+JAujsgU2n8Hde/A6f26RAS8T5c2kc4X3rdZ99nJ2nGqzrVrUZKw4FGiN4pu7Zpzx+l9Ewr4uQd3Y96qzVH7S4eu7ZpHbKufa3z+oRXs3rpZlI3Obvp0acNLFw+m/+1jWb91V9S5r99o+X/H9UxYDCoZia6lCTcem1YlPzfOIW/IYBDNSwo5uld+9UVpKOERcC/LBKhI0287qa5XVqKGFPtvGDp9WCGU4O0UgWLbYxXx7qm1LStmw7a61DG3lxjvzcJvCK3+9bKjqgbw1uMy0X7iecKJlvEzIIGbROJ/5xl1ZVTre+D1adWsiPVbd0Xtv/4bxbUn9fZsx33nHki7FiX88ukv47YVABmJ8UP9N0qlPqERcC/3ZkP0O9WTvUDE8rAh5ipyhNCZ7EUETtpvd47bpyO1Bh4eN5/lG7b7NblBCNlpjRcRKjqU8eyvDmHN5p1c/0rq/OnZfz4ZgP3+8E7MvJ3VltDGEwz3A7xlaRFn9u9CaVEBj49flHSfm3ZUA5bApaIgkYB7uEidZWobIOBeXjQcr9ot4C3tmP7WnTWRaek6HGf278r81VsAKC0OtjnN8eI7tW541k5jIjSNmF5iqQ3Ro1SeYf3dJ2rkcR4C8VLP3LRsVsR5h1Tw00Mr+GTYcTHzX7lkMM9dFD0GnyMcmRBeyXKc/uhe5ZR5bGQsKymKee12PMYdVZZ6NfMgGGUlhfz5tD78/pR9GX/DsUljux3tMqh+B66tTyJhr48jrA0ZZ9JL+MUJNbkXPe9gK2WwY+u6eLuXHOpEOOfF3ZCaaXp3asXff3wg954T7BiTYSM0Au5FbBoSQnGu8x7lLbj51NjROuo/QBK9/jo3ixcPPFX44uDK9nRu0zxqmnNTZqJhMBdxej97c9tWZQ+hdUY/Kz1u385W5xX3MFlXHb93ZEQW5zgVFAjd2pcl/b1H9Spn3PXHeB6lPB5ePPD2LUo4oGsb7js3edfwZHgJvzgPE/eSPxu0J1/fMSQSK4eGhfx67d6S607sxcNZ6Kx1Rr+uMR2vmjqhCaF4EZuGhJ4dr2afTq35zVHdufOtOVHz61/k5a1KaVtWTEX7MqYv2xh5uFS0L2Ph2q2eRrv2op3u3xPxwFOv6mH/2YiAx+4zXZwGytP7deH0fnUi+/bVR0UJ2jUn9uLbdVuBOMcvxe4bWnzIS354aVEh/70ivXRBBy8jsjsPe3fDvojEhJ+KPeZF3nVGX5as3xazvSu19n3OCI2A+4mBV7Qv4+he5XyxaD1fr9qcfCWbyD3hYT/NiguZ9oeTeHjcfEvA7ekv/mYQU5d8n3S0kMICoabWeBJPt6ftfLWmJ7+JD660emwme9vOdq56Q3a3f+f43ZXjve04jyb3AyOTv/dfvzoEEbjgybq0VXcX9KIC8SS2ybjyuJ60LSvh9jdnR6a1L/PuhXppWPUS+gH46aFNoxx0mAiNgHvxFwttT6KkqIDbT+/DTx//HFZ5274Tv27ITd6pTbOUNYpLCgvYXlvjMSQU/T0yJFSKdd+95ij27JC69T/bmY6pQrfdy1uwcM3WyPdRlx1G13bNWbFhh6/u+JGyoS7tyqSAH9WrnIVrtkS+9+nSmgtdw+h9efMJ7KiuYfDd3nPT3Vx3Um82bqvi9jdn85ez+tKsuJAfHLBHyvWcUdvzpa6JEgyhEXD3zXjtib1Yv3UXz3y6ODKtKJL7at3ATthjyP6d+GT+WjbvtDINbhm6L4O6d4jaXqSbbQC2OxzeswMD92zPA+/P8yQm7mW8xsDj1QA/s18XFq3bytQlGyLTvP7WCwfvSbOSQh77aKHHNRKRXMHHXHVkVDikf4XViaZjK3+ZB/V7xNYn0w8sZ/tty4p588ojI9OLC4WqGkM7jx1MUtGmrNh3Sdsquzu81/AI0OhypJsCoRFwN+cfWkGHlqVRAl7oarhx6kOcP6iCRy8YEBlR5McHd4up1et4bf6GvvP3evzCrwfx4sQlcffjlDGtj3sZ57svG+2/p/btzAn77R41qoqXWD3AGf27sv8erX0J+J1n9GFnlbeONX//8YEMqGifMJfYLzUJ3qYynZ+f6DyMueooJi2uK87Zv6ItPcpTjxUJcM6ArsxZuYmZyzc1yLaK9mUUFQjXndTL0/JTbj0xUj5ACQ+hEXC3VjaPk5JW54Fb3+84rQ+9d28ZU/4zngebyGvLNCZmxDyL8Tccx5ad1Rz71w8j0xIJkHv6oXu1p3XzYsbOjhMvMtHr1sdri77E2WcqTtqvU8wI5Ymedx1bNaPCQ8jHK5GyoQlCUJnC2Z678bJnx5b07Fgn2KN8jPvopMklG77MCy1Ki5h/16mel/faHV3JL0In4I6nGs9bK3CFUNqUFXPFcbEt5PHu49o0QijpZFUkGhWkvFUp5a1KefRnAyI5zm7BSRQDf/m3VoXC5z5bTN+u0UWLal3hJLf9w8/sS9uyYi55PvHgsF4KX104eE9aNy/mHx/MT7i/bPWiS/Qwrs5Al/76OJtvSBE1RWkI4RFw+/a/5+wDOaxHh/iC5HFb8RpE67y29EXZE5HX+/izh/Spq7bn7i5emCKN8ILBlXFsSyzgAOcdUsFSV2qYGy91U9qWlXDNib3qBDzJW07Q1EYEPHr6jurMCrhJ8najKNkgPAJu3ywlhRJVta1vlza0bm79DK/1EuKJkUlw0/vdTirqAiipV3Y/GLq0a27v1/uOazyITKrNiavTz2XH9KCosIAH358XmeYuJBWv7SyIcU3jsXfHlvxsUAW/OGyvqOnOIA0D92zHlCXfN3g/TopgQ3oyKkpDCI+ARz5F3yz/q1dD2ZmTjqdXmyBu6gU/wmQShFDi2mQLRMvSIq47qRfH2FkCfhpPnW0kE/BUsWH37BuG7MPkb7+PEnD39pO1Mww9oDNn9e/CwjVbuWP0HLplqOBRZN8Fwh2n942ZfsWxPXlo3HxevfSwjOynxk7VC8ID/89vB2ujopKS8Ai4B+Fz5qUS1LgeuDPPfgyMuepI1m7ZSffylny/dVfc7TgilUhP9+xQxrfrosMTQ/p04qFx8/l5nHBHjE32dtu1KOaXh9d5k35CNqlCKOBBwOO8LTj1Xtq3KOGcAV255OgeUfOTxcAFOG6f3Tm2t+HsAV2jRpMJkutP7s31J3uvtpcKxwPPxAg9bg7J4UC5SngIjYA7t38yqYlUBUxVmCpuDNyeZ89yam0AdGnbPGb5VLaA9RBYsn4bQ+4fH5m2e+tmfHnzCSnWtHAq4x2/z+5R0/2UInUWTdbQ5iWE4sbpbVog8PsUtWMc3O0MIpI18Q6CvTu24ieHVHDREXulXlhRAiClgIvIU8APgNXGmD72tAOBR4GWwGLgfGNMwxJXU+B4Ock8yU5tmtGjvAW3/mC/pNuKn4WSfiNmopBGWUlRg6q0tWtRwsSbjqeDK8XLTyU753Ali7m7u3v3KG/Bgno9IuMdEieHPJH3Hu80OTHoEg/du8NAYYFw95mxoRpFyRZe7qRngCGuaU8Aw4wxfYHXgN9l2K4Ybj+9D788vDJpb7GSogLev+4YjundMem24nmjjoT5CWemCqEAETc93RTk3Vs3i6ln4cS1P7z+GD647uik6ztCu6umJuEyTnpdRfsyFg8fSr+K6GHE4ol/pFhSAjGO96CtGzShcQi4ouSalHeSMeZjYL1rcm/gY/vzWOCsDNsVQ3mrUm774f6eivOkIp6X7QiS196J1nasv8n1O7nQpYPjge/RtjndU/Twc/a7M0kKnfOb90lQbyTew8exoWWcQX+tdeIJuFPLWxvnFCUTpBsDnwn8CHgDOAfoljGLcsRZA7qwaO0W/i+N0pjJPPD2LUo4s18XfpXBOOlhPTrwyfx1ntLX3OGXeHRu05znLzqUfhVt486Pt5fKDmVcekyPyAABDm9eeQQfzF0ddzuVdrnWA7rGryyoKIo/0hXwXwEPisgfgP8C8dM0ABG5GLgYoKIif8tRlhYVcvPQ5LHzRCTLeikskMio45lixAUD+W7Ddk/x+rvP7Mt+nVszaC+reNcXNx0f9y3miL13i5nmEG8/IsKNQ/aJmd6nSxv6dIkv0CfutzvvXnNU3GJbiqL4J633emPMXGPMScaYAcC/gQVJlh1hjBlojBlYXp7bamdnD+gaGTorE4iXGHgAtCgtYm+PIti2rIQrj987UmagY+tmvute1Nfv1h7GjEyGireiZI607kYR6WiMWS0iBcAtWBkpec9fMzyeXmPtf+cecNf5na9eMpg9OzRs1BpFUTKHlzTCfwPHALuJyDLgNqCliFxuLzIKeDowC/OYSCNmtl3wgNnl6hbvNDoOrNTOJYqST6QUcGPMTxLMeiDDtoQOZ9iy4/bdPcWS4aK6JvqB5B4hXlGU/CBEPTHzjz5d2rDo7lOzPrp70LgLU8Wrva4oSu7RHhUNpLGJN0CVKwZe6iM3XlGU7KEeuBLDoXu15+Nv1jDqssMoEGmUDylFaQyogCsxXHp0D3504B50a5/ZMq+KomQWfTdWYigoEBVvRQkBKuCKoighRQVcURQlpKiAK4qihBQVcEVRlJCiAq4oihJSVMAVRVFCigq4oihKSJFsVtITkTXAt1nY1W7A2izsxw/5aBPkp11qkzfUJm80Bpv2NMbEDKiQVQHPFiIyyRgzMNd21CcfbYL8tEtt8oba5I3GbJOGUBRFUUKKCriiKEpIaawCPiLXBsQhH22C/LRLbfKG2uSNRmtTo4yBK4qiNAUaqweuKIrS6FEBVxQl9EgejjqSDZtCLeBN9aSlQz7apTZ5Q23yRHGuDYhD4DaFTsBFZH8ROQbA5EkAPx9tgvy0S23yhtrkDREZLCKvAH8Vkf1EJOcjcGfTptA0YopIAfAQcBywBJgIvGGMmSQiBcaY2qQbaCI25atdapPaFIBdHYExtm3dgC7AJGPM4yIiuXjIZNumMHng7YBWwL7A+cA64DoRaZmrCyhPbQJoC7TMM7vakn/HKh/Pn9rknQOBb4wxTwN/A0YBp4lIL2OMyVGYJ6s25bWAi8jZInKZ/bU1MBgoM8asAUYC64HL7WWzcrJEpL+I9LK/tgEOy7VN9r66i4gzkGV74PBc2yUie4pIs3o25cP5O0JEetpf2+aJTXqde7PpJyLyJxH5kT1pKjBQRHoYY7YCXwKTgN9CdsI8InK6iNwkIkPtSdOyaVNeCriItBSRkcD1wPciUmSMWQR8AlxtL7YC6+nWT0T2CPpkicheIjIaeBh4TkRONMYsBD7NlU22XZ1F5GPgeeANEelrjJkPfAhcmwu77Ljf68AzwP9s72M+8Dm5PVYHAR8DPxGR1saYBcBnubJJr3PPNomIXALcACwG7hWRXwNbgH8BV9mLbgDeA8pEpHPANpXb1/i1WA+zp0Xk7HoPuCuzYVNeCjhW7GiVMWaQMebfQI09/RngcBHpboypBlYBO4DmWbDpemCaMWYw8AbwK3v6U7ZNe2XLJpfH82PgS2PMYcD7wDAR6Y91rAZl61g5NonIPsAjwDhjzLHAV1jxQMjysYrjGXYBxmJd90fngU15cZ27bMqb69zBfkAMBobboYnLgWOA44G3gJ4icoIdzlmHdZ43BmxWD+ATY8xRxphHgeuAa+x5/wb2EZHjg7YpbwRcRDpJXWvtgUBXe/plwG0icgQwG8sT+CuAMWYmsCewM0CbnFSg7UCV/bk1MMd+Ff8E6zUpKzbZNHN9Lrb3PRxYjXVhr8JqbLonS3Y5Nm0EhhljHrC/347lfZQDXwBTcmCTw/fAPKAWOFhEmmO9qWTz/NW36QDy4Dp3bLKFfCt5cJ2LyM9F5GgRaW9PmgN0sd9S3gNmAoOANcCLwP22nccDApQEZNMxdqhyMpb3j61bs+1/ADOAl4AHgrYp5wIuIseLyHisV7ZH7MlTgO9E5CmsJ+8G4GbgNOA+oKOIPCQiM7Hqi2/MZBzOZdOD9uTxWE/6qcAQoBDrwjkauBvoHKRNtl0nishYrFfI8+zJi4B1IlJhf38Z6IsVt7wb6Coi/wjwWNW36VxjzApjzGf19tEX2GWMWWOM2QL8OYs23VPvODm2TMGqQ9EMuAk4G+s47ZFFm35iT54CrMjhde4+dwaYAOydi+vcDpV0FpFxwIVYDab/EJHWwFKgI+C0X7wE9AE6GGOeB14AhgHnATcYYzYEZNNPgceBNsaYVSJSaIypwWrgbQNgjKk1xjyDJfAZtykKY0zO/gG9sLzEs7FOzhisxrcKLC9tMlBsL3sBMML+3BGrUeVHWbDpbeBIe14PYFS9ZW8FHrI/7x6UTfb2e9p2nQb0w7qprgA6A08DP6AuLfQZ4A/2504BHiu3Tc8DN9nznPN2knOM6q0X5PmLZ9Mt9rwfAj8DumOFdjYBP8/BcXoB65W7CCtTIRfXebzr6Xp7Xu9sX+dAof23F/C8/bkI+CfwLNZb5lP28WlT7zq/s942SrJk0z+c41NvmX8B5zrXUlA2uf8VkWXEyinFWLGhg4AvjDGv2k/ZrcBiY8xyEXkLGAicg3VxfQWcJVbe6WqsUEE2bNoCLBKREvvzUhHZ1xgzBxgHXC0iYoxZhRW2yBguuw4FJhtj3rDnvY918z+LFZo40rbvQ+BNrAchxpiVwMos2fQBcJ+IPGGfI7BeHyfa828FnjbGLCO485fIpkewHna/AW4DRgPvYoV3CnNwnP4GPIkVZ+5H9q9zt03vYR2n57Aa5bJynYtIEdZbWaF9z7fGbgswxlSLyBVYjab7YR2f07HCTndjhcI+dbZljNmVJZv+DytCcLQx5iN7NUcn/gycKSJDjDHLMmVTIrIaQhGRXwLLsOKiANOBASLyOFbcqCPW69wTxpgPgQew8k1vxHplmmBvJ5OvkV5s+guWJ7AW6AD8n4hcBTyG1cKcceLYNQMre6LS/l6E1SL/F9uOZcDfRGQYcD+WkGfbpmJgAXac1D5PA7Aavz4C9sGKQ2fbpkVYov0q1vkabIy5GiuOuhkrPplNm4psm+4xxnyMdb6yfZ3HO04L7fmbsdI+A73OReRorLePdsB8e99VwLEicghEHjZ/Bv5irNj3COAIEZlor/dhDmwytk1/tNcpxGrsfRVL7I+1nZTgCdK9d72OtARex0r5mQL0tqeXA78DLrW/N8MSysH290OxcigH59imdVgxtwp7+WeBQVk6VvvY0+/HauH+BCs00BerFb7cnn8K8AfgiBzbNBqr1V3sZT8A+uXYpjHAbq71i/Pg3HWy5x+cxes81XFqgRXTvTLg6/xI4IJ63/8JXAr8AuvtACwnsxPwClBpT2sLdMkDm/6D1ZDbwz6W/YOwKam9Wd0ZVNh/hwP/rncwnsSOM9vTHgJ+kGc2PQycmqNj9bL9uRDLMzrC/t4NKw7YLA9tKsR68AV6Ufuw6Wmg1Dm/eWJTPp67Zwk4blvPpjKglLo48vnA3fbnacCV9ueBzr2ZZza9lA2bkv3LagjFGLPE/ng/0ENETjHWK9J8YISI9BaRm4AjgFl5ZtPhWKlMWcFl114icrKxWrs3GmMm2PMuAbZRl/aVTzaJMWaHMWZKnti0Hai21wm0+3fIz91W6vLRg7ZpmzFmp20HwIlYaYEAvwT2FZE3sd4SAr2O0rRpMuS4MmOunhxYr4vj633/K1br/L+BbmpTjF0f1ft+CFbjV+QVXG1Sm8JqE9abQAFW+KanPa0nVqjkCAIKl4TNpnj/clKN0G5hrxWRV7Fa2bdhxZNmGGO2Z92gPLUpjl0rsDpOvAfMM1ZXcLVJbQq7TU4nlyeA17AaBNdhhSs2qU2JyUlHHvsCKsPK8DgXWGKM+SKXQpmPNsWx6ye2XW/n6mZTm9SmAGwyWKmU52PVFnnNGHNhLoUyH22KR9bzwOtxGVZc60RjTJDdzv2QjzZBM2D5rwAAAaFJREFUftqlNnlDbfLGMqxeqPepTd7J2YAOksNC8InIR5sgP+1Sm7yhNilBEpoReRRFUZRocl7MSlEURUkPFXBFUZSQogKuKIoSUlTAFUVRQooKuNJoEZG2Yg8WLCJ72J1XFKXRoFkoSqPFLpX6pjGmT45NUZRAyGVHHkUJmuFYBcqmYY2Fua8xpo+I/AJrYIBCrBLBf8PqNn0BVtfyU40x60WkB1YVynKs0gq/McbMzf7PUJT4aAhFacwMAxYYYw7Cqu9enz5Y4xseAtwJbDPG9AM+A35uLzMCq/bFAKzR2v+ZFasVxSPqgStNlXHGmM3AZhHZCPzPnj4DOEBEWmKN/fhKvWqhpdk3U1ESowKuNFXq17aorfe9Fuu+KAA22N67ouQlGkJRGjObgVbprGhXnVskIueAVV5URA7MpHGK0lBUwJVGizFmHfCJiMwE7k1jE+cDF4nIV1gjRJ2WSfsUpaFoGqGiKEpIUQ9cURQlpKiAK4qihBQVcEVRlJCiAq4oihJSVMAVRVFCigq4oihKSFEBVxRFCSkq4IqiKCHl/wELzEjE9KCjzwAAAABJRU5ErkJggg==\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"area_weighted.plot()"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"4. Remove seasonal cycle."
]
},
{
"cell_type": "code",
"execution_count": 12,
"metadata": {},
"outputs": [],
"source": [
"climatology = area_weighted.sel(time=slice('1980', '2010')).groupby('time.month').mean('time')"
]
},
{
"cell_type": "code",
"execution_count": 13,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"[<matplotlib.lines.Line2D at 0x2ae1d11ed208>]"
]
},
"execution_count": 13,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXAAAAEKCAYAAAALoA6YAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3dd3yV9d3G8c83mySskAABAiEQNsgIiCiKA6GgKFpHWy22tlgnts5HnzqrtdU6Wu2w6qOtVFq3ggMFBLEqmzAS9iZAwgyB7N/zRw42UiAnISf3OSfX+/XKy+RwxnUcl7/c575/X3POISIioSfC6wAiIlI3KnARkRClAhcRCVEqcBGREKUCFxEJUSpwEZEQVWOBm1mamc0ysxwzW2Fmk3y3P2Bm28xsie9rTODjiojIEVbTeeBmlgqkOucWmVlTYCFwMXA5cNA590TgY4qIyNGiarqDcy4PyPN9X2hmOUD7QAcTEZETq3EF/q07m6UDc4A+wC+Aa4ADwALgNufc3hM9Pjk52aWnp9ctqYhII7Vw4cIC51zK0bf7XeBmlgjMBh5xzr1lZm2AAsABD1N1mOXHx3jcRGAiQMeOHQdt2rSp7u9CRKQRMrOFzrmso2/36ywUM4sG3gQmO+feAnDO7XTOVTjnKoG/AkOO9Vjn3PPOuSznXFZKyn/9D0REROrIn7NQDHgRyHHOPVnt9tRqdxsPLK//eCIicjw1fogJnA5cDSwzsyW+2+4Bvmdm/ak6hLIRuC4gCUVE5Jj8OQtlLmDH+KMP6j+OiIj4S1diioiEKBW4iEiIUoGLiIQoFbiHSsor+GJtAf/3xQZ2HyzxOo6IhBh/zkKReuKcY+PuQ8xZnc/s1fl8uW43h8sqAHh25loevrgPY/qm1vAsIiJVVOABdrCknH+vLWDOmqrS3rLnMACdWsXz3UEdOLNbCilNY7nv3eXcMHkRY/ul8tC43rRKjPU4uYgEOxV4PausdKzMO8Ds1fnMWZ3Pwk17Ka90xMdEMqxLK346PIMzM1NIT0741uPeun4Yf5mznqc/Xc1X63ZrNS4iNarVZlYnKysryy1YsKDBXq+hFBws4fM1+cxZXcDna/IpOFgKQK/UZpzZLYUzuyWT1SmJmKiaP3JYtaOQ219fyrJt+7UaFxHg+HuhaAVeB2UVlSzatLdqlb0mn+XbDgCQlBDD8MxkzsxMYXi3ZFo3jav1c3dv25S3b/j2avxXF/fhO1qNi8hRtAL305Y9h5hd7cPHgyXlREYYAzu24KxuKZzZLYU+7ZoTEXGsi1brRqtxEQGtwGvtUGk5X63fzZzVBcxenc+GgiIA2rdowoWntOOsbikM69qKZnHRAcug1biInIhW4EepqHTc8fpSpmbnUVpRSVx0BEMzWnFmZgpndU8hIzmBqg0aG9bRq/GHL+pDUkJMg+cQkYanFbif/vTZWt5avI3vDenImL5tGZyeRFx0pNex6N62KW/dMIy/zF7HMzPWaDUuIroSs7pFm/fy1KdrGHdKOx4d34fhmSlBUd5HREdGcNM5mUy9eTjtWjTh+smLuOkfi9hTVOp1NBHxgArcp7C4jElTFtO2WRy/Gt/Hk8Mk/jqyGr/9/G58vGIHI5+czYfL8ryOJSINTAXuc9+7K9i29zDPXNk/oB9M1pcjq/H3bz6D1BZxWo2LNEIqcOCdxdt4e/E2bjk3k6z0JK/j1EqPts14+4bTtRoXaYQafYFv3n2I/31nOVmdWnLT2V29jlMnx1qN3/zaYq3GRcJcoy7w8opKbv3nYgx46or+REWG9t+OI6vx20Z246PleZz/1Gw+Wq7VuEi4Cu3GOkm/n7mWRZv38cglfUlLivc6Tr2Ijozg5nMzee+mM2jbPI6fvarVuEi4arQFPm/DHp6duYZLB3Zg3CntvI5T73qmajUuEu4aZYHvP1zGz/+5hLSkeB68qLfXcQKm+mq8TTOtxkXCTaMrcOcc97y9jJ0HinnmygEkxob/xag9U5vxzo2n8wvfanzcs3MpLC7zOpaInKRGV+BvLNzKtOw8fj6yG/3TWngdp8FER0Zwy7mZvHrtqWzbd5jffrTK60gicpIaVYFvKCji/vdWcGrnJH52Vhev43ji1IxWXDMsnb9/tYkFG/d4HUdETkKjKfDS8komTVlMdGQET13Rn8h63Lc71Nx+fnfat2jCXW9mU1Je4XUcEamjRlPgT326muyt+3nskr60a9HE6zieSoiN4pHxfViXX8RzM9d6HUdE6qhRFPi/1xbw59nr+N6QNG2/6jOie2vGD2jPHz9bR+6OA17HEZE6CPsC31tUys//tYTOyQn88oJeXscJKr+8oBfNmkRz95vLqKhsuMEeIlI/wrrAnXPc/VY2e4pK+f2VA4iPCf9TBmsjKSGG+y/sxZIt+3jl3xu9jiMitRTWBf7avC18vGInd47qQZ/2zb2OE5TGndKOEd1TeGL6KrbsOeR1HBGphbAt8LW7Cnlo6gqGZyZz7RmdvY4TtMyMR8b3BeDed5bTkDNSReTkhGWBl5RXcMtrS4iPieJ3l51CRCM+ZdAf7Vs04c5R3ZmzOp93lmzzOo6I+KnGAjezNDObZWY5ZrbCzCYd9ee3m5kzs+TAxaydxz9axcq8A/z20n60bhbndZyQcPVp6Qzs2IKH3l/J7oMlXscRET/4swIvB25zzvUEhgI3mlkvqCp3YCSwOXARa2fO6nxemLuBH57WifN6tfE6TsiIjDAeu7QfB0vKeWjqSq/jiIgfaixw51yec26R7/tCIAdo7/vjp4A7gaA4cFpwsIRf/Gsp3dokcs+Ynl7HCTnd2jTlhhFdeXfJdmbl7vI6jojUoFbHwM0sHRgAfG1m44BtzrmlAchVa8457nwjmwPFZTxz5QDioiO9jhSSbji7C5mtE7n37WUcLCn3Oo6InIDfBW5micCbwK1UHVa5F7jPj8dNNLMFZrYgPz+/zkFr8rcvNzEzdxf3fKcHPVObBex1wl1sVCSPXdqPvAPFPPGxdiwUCWZ+FbiZRVNV3pOdc28BXYDOwFIz2wh0ABaZWdujH+uce945l+Wcy0pJSam/5NXk7jjAIx/kcHb3FCYMSw/IazQmgzq1ZMJp6bzy5UYWbtrrdRwROQ5/zkIx4EUgxzn3JIBzbplzrrVzLt05lw5sBQY653YENO0xFJdVMOm1JTSLi+bxy06hKq6crNtHdSe1WZx2LBQJYv6swE8HrgbOMbMlvq8xAc7lt19/kMOqnYU8cVk/khNjvY4TNhJjo3jkkr6s3XWQP322zus4InIMNW4O4pybC5xwWetbhTe4GTk7eeXLTVx7RmdGdG/tRYSwdnb31lzUvx3PzVrLmL6pdGvT1OtIIlJNyF6JuetAMXe8kU3P1GbcObq713HC1n0X9CIxNoq73szWjoUiQSYkC7yy0nHb60s5VFrO76/sT2yUThkMlFaJsdx3YS8Wb97H37/c6HUcEakmJAv8pS828PmaAn55QS8y9Wt9wF3cvz1ndkvhtx+vYtu+w17HERGfkCvw5dv285uPchnZqw3fH9LR6ziNgpnx6Pg+ANz79jLtWCgSJEKqwA+VljNpymKSEmL4zaX9dMpgA+rQMp7bz+/OZ6vyeW/pdq/jiAghVuAPT81hfUERT17en6SEGK/jNDoThqXTP60FD76/kj1FpV7HEWn0QqbAP1q+g9fmbWbimRmc3jVodq5tVCIjjN9c2o8Dh8t4WDsWinguJAo8b/9h7n4rm77tm3PbSJ0y6KXubZtyw4guvL14G5+t0o6FIl4KiQJ/7MNcSssreebK/sREhUTksHbjOV3pkpLAvW8vp0g7Fop4JiTa8MFxvfnrD7PISEn0OopQtWPhby7tx/b9h3liunYsFPFKSBR4i/gYHfcOMlnpSVw9tBMv/3sjizdrx0IRL4REgUtwumNUd9o2i+PuN5dRWl7pdRyRRkcFLnXWNC6aX13ch1U7C/nzbO1YKNLQVOByUs7t2YYLT2nHszPXsnZXoddxRBoVFbictPsv7EV8bCR3vbmMSu1YKNJgVOBy0pITY/nl2F4s3LSXyV9v8jqOSKOhApd6ccnA9gzPTOaxD3PZrh0LRRqEClzqRdWOhX2pdPDLd5Zrx0KRBqACl3qTlhTPbed3Y0buLt7PzvM6jkjYU4FLvfrR6Z05pUNzHnxvBXu1Y6FIQKnApV5FRhiPXdqP/YfL+NW0HK/jiIQ1FbjUu56pzfjZWV14c9FW5qzO9zqOSNhSgUtA3HROVzJSErjn7WUcKtWOhSKBoAKXgIiLjuSxS/qxde9hnvpktddxRMKSClwCZkjnJK4cnMb/fbFRl9mLBIAKXALqjlHdaRITyYPvr9S54SL1TAUuAdUqMZZfjOzG52sKmL5yp9dxRMKKClwC7qqhnejWJpGHp66kuKzC6zgiYUMFLgEXHRnBAxf2Zuvewzw/Z73XcUTChgpcGsSwrsmM6duWP362lm3a7EqkXqjApcHcM6YnAI/qCk2ReqEClwbToWU815/VlWnL8vj3ugKv44iEPBW4NKjrzsqgQ8smPPjeSsorNAhZ5GTUWOBmlmZms8wsx8xWmNkk3+0Pm1m2mS0xs+lm1i7wcSXUxUVH8r9je7FqZyGvfqXpPSInw58VeDlwm3OuJzAUuNHMegGPO+f6Oef6A1OB+wKYU8LIqN5tGJ6ZzJOfrGb3wRKv44iErBoL3DmX55xb5Pu+EMgB2jvnDlS7WwKgy+zEL2bG/Rf24lBpBU9MX+V1HJGQVatj4GaWDgwAvvb9/IiZbQF+gFbgUgtdWzflmmHpTJm/hWVb93sdRyQk+V3gZpYIvAncemT17Zy71zmXBkwGbjrO4yaa2QIzW5Cfr72h5T9uOS+TVgmx3P/ecior9QucSG35VeBmFk1VeU92zr11jLv8A7j0WI91zj3vnMtyzmWlpKTUPamEnWZx0dw1ujuLNu/jnSXbvI4jEnL8OQvFgBeBHOfck9Vuz6x2t3FAbv3Hk3B36cAO9E9rwa8/zKWwuMzrOCIhxZ8V+OnA1cA5vlMGl5jZGOAxM1tuZtnA+cCkQAaV8BQRYTw4rjf5hSU8O3Ot13FEQkpUTXdwzs0F7Bh/9EH9x5HG6JS0Flye1YGXvtjA5YPT6JKS6HUkkZCgKzElKNw5ugdxURr8IFIbKnAJCsmJsdw6shtzVufzac4ur+OIhAQVuASNH57WiczWGvwg4i8VuASN6MgI7r+wN5v3HOKFzzX4QaQmKnAJKmdkJjO6d1uem7WO7Rr8IHJCKnAJOveO7Umlczz6gQY/iJyIClyCTlpSPD87qwtTs/P4av1ur+OIBC0VuASl60d0oX2LJjzw3goNfhA5DhW4BKWqwQ89yd1RyD/mbfY6jkhQUoFL0Brdpy2nd23F76avZk9RqddxRIKOClyCVtXgh94cLCnX4AeRY1CBS1Dr1qYpE05L57V5m1m+TYMfRKpTgUvQm3ReJknxMdz/3grtkyJSjQpcgl7zJtHcNboHCzft1eAHkWpU4BISvjuoA6d0aM6vP8jlYEm513FEgoIKXEJCRITxwLje7Cos4Q8z13gdRyQoqMAlZAzo2JLvDurAS3M3sD7/oNdxRDynApeQcpdv8MNDUzX4QUQFLiElpWksk87L5LNV+czM1eAHadxU4BJyJgxLp2vrRB7S4Adp5FTgEnKqBj/0YtPuQ7w4d4PXcUQ8owKXkDQ8M4Xze7Xh2ZlryduvwQ8SvPYdKmXCS/NYsb3+ryRWgUvI+uUFvahwjl9/kOt1FJHj+tPsdcxZk09khNX7c6vAJWSlJcXzszMzeG/pduZt2ON1HJH/smN/MS9/sZHx/dvTo22zen9+FbiEtOtHdKVd8zjuf28FFZU6rVCCyzMz1lDpHD8f2S0gz68Cl5DWJCaSe8f2IifvgAY/SFBZl3+Qfy3Ywg9O7URaUnxAXkMFLiFvTN+2nJbRit9NX8VeDX6QIPHk9NXERkVw0zldA/YaKnAJeWZV+6QUFpfzu080+EG8l711H9OW5fGT4RkkJ8YG7HVU4BIWurdtytVDO/GPrzcH5HQtkdr47UerSEqI4afDOwf0dVTgEjZ+PrIbLeJjeECDH8RDc9cUMHdtATee3ZWmcdEBfS0VuISN5k2iuXNUd+Zv3Mt7S7d7HUcaIeccv/04l/YtmvCDUzsG/PVU4BJWLstKo2/75jz6QQ5FGvwgDeyj5TvI3rqfW8/LJC46MuCvpwKXsBLpG/yw80AJz85a63UcaUTKKyp5fPoqMlsncsnADg3ymjUWuJmlmdksM8sxsxVmNsl3++Nmlmtm2Wb2tpm1CHxckZoN6tSSSwa254XP17OhoMjrONJIvLFwK+vzi7hjVPeAXDZ/LP6swMuB25xzPYGhwI1m1gv4BOjjnOsHrAb+J3AxRWrn7u/0IDYqkoenrvQ6ijQCxWUVPP3pGgZ2bMHIXm0a7HVrLHDnXJ5zbpHv+0IgB2jvnJvunDtykPEroGF+ZxDxQ+umcUw6N5OZubuYmbvT6zgS5l7590Z2HCjmrtE9MGuY1TfU8hi4maUDA4Cvj/qjHwMf1k8kkfoxYVg6GSkJPPT+SkrKNfhBAmP/4TL++Nk6RnRP4dSMVg362n4XuJklAm8CtzrnDlS7/V6qDrNMPs7jJprZAjNbkJ+ff7J5RfwWExXBAxf2ZqMGP0gAPT9nHfsPl3HHqO4N/tp+FbiZRVNV3pOdc29Vu30CcAHwA3ecKyecc88757Kcc1kpKSn1kVnEb2d2S2Gkb/DDjv3FXseRMLPrQDEvzd3IuFPa0btd8wZ/fX/OQjHgRSDHOfdktdtHA3cB45xzhwIXUeTk/HJsL8orHb/+MMfrKBJmfj9zDWUVldx2fmC2i62JPyvw04GrgXPMbInvawzwLNAU+MR3258DGVSkrjq2iue6MzN4d8l25m/U4AepHxsLipgybwvfG9KRTq0SPMkQVdMdnHNzgWN9rPpB/ccRCYwbRnTlzYVbuf/dFbx/8xkNdp6uhK8nP1lNdGQEN58buO1ia6IrMaVRODL4YWXeAV7T4Ac5Scu37ee9pdu59ozOtG4a51kOFbg0GkcGPzyhwQ9ykh7/eBUt4qOZeFaGpzlU4NJomBn3j+tFYXE5T36y2us4EqK+XLeb2avzuWFEF5oFeLvYmqjApVHp0bYZVw/txOSvN7Fy+4GaHyBSzZHtYlObx/HD09K9jqMCl8bn5+dp8IPUzfSVO1m8eV+DbRdbExW4NDrN46O5Y1R35m3co8EP4reKSsfjH68iIyWBSxtou9iaqMClUbpcgx+klt5atJW1uw5yx/ndiYoMjuoMjhQiDaz64IfnNPhBalBcVsFTn6zmlA7NGd2nrddxvqECl0brP4MfNrBRgx/kBF79ahPb9zf8drE1UYFLo3b36B7EREVo8IMcV2FxGc/NWsvwzGSGdU32Os63qMClUWvdLI5bzu3KjNxdzMrd5XUcCUJ//XwDew+VceeoHl5H+S8qcGn0rhnWuWrww1QNfpBvyy8s4YXP1zO2Xyp9OzT8drE1UYFLoxcTFcH9F/ZmQ0ERL83d6HUcCSLPzVpLSXklt430ZrvYmqjARYCzfIMf/jBzDTsPaPCDwJY9h5j89SYuz0ojIyXR6zjHpAIX8flm8MMHGvwgVdvFRphx63mZXkc5LhW4iM+RwQ/vLNnOAg1+aNRy8g7wzpJt/Oj0zrRp5t12sTVRgYtUc/2ILrRrHsd9766golL7pDRWT3y8iqaxUVx/Vhevo5yQClykmviYKO4Z25OVeQeYMl+DHxqj+Rv3MCN3F9eP6ErzeG+3i62JClzkKGP7pjI0I4knPl7FvkMa/NCYOOf4zYe5tG4ayzXD0r2OUyMVuMhRzKr2STmgwQ+NzszcXSzYtJdJ52XSJMb77WJrogIXOYYjgx9e/UqDHxqLikrHbz9aRefkBC7PSvM6jl9U4CLH8c3gh/c1+KExeHfJNlbtLOS287sRHSTbxdYkNFKKeOCbwQ8b9vB+dp7XcSSASsorePKT1fRp34wxfVK9juM3FbjICVyelUaf9s14dJoGP4Sz177ezNa9h7lzVA8iIoJnu9iaqMBFTiAywnhwXB92HCjmj59p8EM4OlhSzh9mrmVYl1YMzwyu7WJrogIXqcGRwQ9/naPBD+Hoxc83sLuolDuDbFiDP1TgIn44MvjhV9M0+CGc7D5Ywl8/X8/o3m3pn9bC6zi1pgIX8cORwQ+f5uxi1ioNfggXf/xsHYdKy7l9VHevo9SJClzET98Mfnh/JaXllV7HkZO0de8h/v7lJi4blEbX1sG5XWxNVOAifoqJiuC+C3pVDX74YoPXceQkPf3pGjCYFMTbxdZEBS5SCyO6t+a8nm34wwwNfghlq3cW8tairUw4rRPtWjTxOk6dqcBFaum+C3pRVul47MNcr6NIHThX9c8uISaKG0Z09TrOSVGBi9RSx1bxTByewduLt2nwQwh6btZaZubuYtJ5mbRMiPE6zkmpscDNLM3MZplZjpmtMLNJvtsv8/1caWZZgY8qEjxuOLtq8MP/vrOc4jJNsg8V07LzeGL6asYPaM+1Z3T2Os5J82cFXg7c5pzrCQwFbjSzXsBy4BJgTgDziQSl+JgoHhnft2rzo38tpVLTe4Le0i37+MW/ljCoU0t+fUnfkLto51hqLHDnXJ5zbpHv+0IgB2jvnMtxzq0KdECRYHV2j9b8z3d6MG1ZHk99qn3Dg9n2fYf5yd8WkNI0lr9cPYi46ODf69sfUbW5s5mlAwOArwMRRiTU/HR4Buvzi/jDzLVkpCQwfkAHryPJUYpKyrn2lQUcLq1g8k9OJTkx1utI9cbvDzHNLBF4E7jVOef3DvdmNtHMFpjZgvz8/LpkFAlaZsZDF/XhtIxW3PXGMubrQ82gUlHpmDRlCat2HODZ7w+gW5umXkeqV34VuJlFU1Xek51zb9XmBZxzzzvnspxzWSkpKXXJKBLUYqIi+PNVg+jQsgnX/X0hm3cf8jqS+Pz2o1w+zdnJ/Rf2ZkT31l7HqXf+nIViwItAjnPuycBHEgk9zeOjefGawVQ6x49fmc/+w2VeR2r0/jl/M3+Zs54fntaJCSEwoLgu/FmBnw5cDZxjZkt8X2PMbLyZbQVOA6aZ2ccBTSoS5DonJ/DnqwaxaXcRN/1jEeUV2i/FK1+u2829by9neGYy913Qy+s4AePPWShznXPmnOvnnOvv+/rAOfe2c66Dcy7WOdfGOTeqIQKLBLOhGa14ZHxfPl9ToFmaHtlQUMTPXl1IenICz/1gIFEhMt+yLmp1FoqI1OzyrDTW5xfx59nr6JKSyI9OD/0LRkLFvkOlXPvyfCIMXpowmGZx0V5HCigVuEgA3DmqOxsKDvLw1JV0ahXPOT3aeB0p7JVVVHLD5EVs3XuYyT89lY6t4r2OFHDh+7uFiIciIoynruhPr3bNuPkfi8nd4feZt1IHzjnue3cF/163m19f0pfB6UleR2oQKnCRAImPieKFHw4mMS6Ka19ewK5CbT8bKC/O3cBr8zZz49lduHRQ47mYSgUuEkBtm8fx4oTB7CkqZeLfFmrjqwD4dOVOHvkgh+/0acttI0NzNFpdqcBFAqxP++Y8fWV/lm7dx+2va+Or+rRy+wFumbKYPu2a8+Tl/YmICP0NqmpDBS7SAEb1bstdo3swNTuPp2es8TpOWNhVWMxPXplPs7hoXpiQRZOY8NigqjZ0FopIA7nuzAzW5x/k9zPWkJGcwMUD2nsdKWQVl1Xw078tZO+hMl7/2Wm0aRbndSRPaAUu0kDMjF9d3JehGUnc+Ua2pvnUkXOO219fSvbWfTx9ZX/6tG/udSTPqMBFGtCRja/a+za+2rJHG1/V1tOfrmFqdh53je7BqN5tvY7jKRW4SANrER/DixOyKK90/Pjl+Rwo1sZX/np3yTaembGGywZ14LozM7yO4zkVuIgHMlIS+dNVA9lQUMSNk7XxlT8WbtrLHW9kM6RzEo+MD4+RaCdLBS7ikWFdknlkfB8+X1PAQ1NXeh0nqG3Zc4jr/r6A1OZx/OWqQcREqbpAZ6GIeOqKwR1Zn1/EX+asJyM5gWu08dV/KSwu4yevLKC0vJIpEwfTMiHG60hBQwUu4rE7R/dgfUERD01dSafkBM4Ow8kxdVVR6bjltcWszT/IKz8aQtfWiV5HCir6PUTEY5ERxtNX9KdH26qNr1btKPQ6UtB4ZFoOs1bl8+C43pyRmex1nKCjAhcJAgmxUbx4TRbxMZH8+OX55BeWeB3Jc69+tYmXvtjAj0/vzFVDO3kdJyipwEWCRGrzJrw4YTC7i0qY+PcFjXrjq7lrCrj/vRWc06M1947t6XWcoKUCFwkifTs05+kr+rN48z7ueCO7UY5kW7vrINdPXkhm60R+/70BRDayDapqQwUuEmRG90nlztHdeX/pdp7+tHFtfLW3qJRrX5lPbFQEL0zIIjFW51mciP7uiASh68/qwvr8Ip6ZsYaMlAQu6h/+G1+Vlldy3asLydtfzJSJQ+nQMvxHop0srcBFgpCZ8ej4vgzpnMQdb2SzcFN4b3zlnOPet5cxb8MenrjsFAZ2bOl1pJCgAhcJUjFREfzlqkGkNo9j4t/Cc+Mr5xxfrd/NDZMX8frCrUw6N5Nxp7TzOlbIUIGLBLGWCTG8dM1gyioqufaV8Nn4ak9RKX+ds55zfzebK5//ii/WFnDzOV259bxMr6OFFGvIT7mzsrLcggULGuz1RMLFF2sLmPDSPBJioxjVuw1j+7VjWJdWREeGzhrMOceX63fz2rwtfLx8B6UVlQxOb8n3hnRkTN9U4qIb30Qdf5nZQudc1n/drgIXCQ3zNuxhyrzNTF+5k4Ml5bSIj2Z077aM7ZfKaRmtiArSMi84WMKbC7cyZf4WNhQU0bxJNJcO7MD3hqSR2aap1/FCwvEKXGehiISIIZ2TGNI5ieKyCj5fU8C07O28v3Q7U+ZvISkhhlG923JBv1RO7ZzkeZlXVlattv8xbzPTV+ygrMIxJD2JW87tynf6aLVdX7QCFwlhxWUVzF6dz7TsPD7N2cmh0gpaJcQwuk/VyvzUzq0a9EKY/MIS3li4lSnzN7Np96VweeEAAAdfSURBVCFaxP9ntd21tVbbdaVDKCJhrrisgs9W7WJqdh4zcnZxuKyC5ERfmfdtx5DOSQEp88pKxxfrCnht3mamr9hJeaXj1M5JfP/Ujozq3Var7XqgAhdpRA6XVjBr1S6mZecxI3cnxWWVJCfGMqZvW8b2TSUr/eTLfFdhMa8vqFptb9lzmJbx0Xx3UAeuGNxR277WMxW4SCN1qLScmblVZT4zdxcl5ZW0bhrLmL6pjO2XyqCOLYnws8wrKx2fry3gta8382lO1Wp7aEYS3z+1E6N6tyE2SqvtQFCBiwhFJeXMyN3FtOztzFqVT2l5JW2aVZX5Bf1SGZB27DLfeaCY1xdsYcr8LWzde5ikhBguG9SBKwankZGi1XagqcBF5FsOlpQzI2cnU7PzmL0qn9KKSlKbx32zMu/Xvvk3q+0ZubuoqHQM69KK75/akZG9tNpuSHUucDNLA/4GtAUqgeedc8+YWRLwTyAd2Ahc7pzbe6LnUoGLBKfC4jJm5FR9ADpndVWZx0ZFUFJeSauEGL6b1YErB3ekc3KC11EbpZMp8FQg1Tm3yMyaAguBi4FrgD3OucfM7G6gpXPurhM9lwpcJPgdKC7j05U7mb9xD2d0TWFkrzaaAu+xOl/I45zLA/J83xeaWQ7QHrgIGOG72yvAZ8AJC1xEgl+zuGguGdiBSwZ28DqK1KBW/1s1s3RgAPA10MZX7kdKXqO0RUQakN8FbmaJwJvArc65A7V43EQzW2BmC/Lz8+uSUUREjsGvAjezaKrKe7Jz7i3fzTt9x8ePHCffdazHOueed85lOeeyUlJS6iOziIjgR4GbmQEvAjnOuSer/dF7wATf9xOAd+s/noiIHI8/uxGeDlwNLDOzJb7b7gEeA/5lZtcCm4HLAhNRRESOxZ+zUOYCx7vO9tz6jSMiIv7SyZ0iIiFKBS4iEqIadC8UM8sHNjXYC56cZKDA6xABEs7vDcL7/em9ha6TeX+dnHP/dRpfgxZ4KDGzBce6dDUchPN7g/B+f3pvoSsQ70+HUEREQpQKXEQkRKnAj+95rwMEUDi/Nwjv96f3Frrq/f3pGLiISIjSClxEJESpwKsxszQzm2VmOWa2wswmeZ2pvplZpJktNrOpXmepb2bWwszeMLNc3z/D07zOVF/M7Oe+fyeXm9lrZhbndaaTYWYvmdkuM1te7bYkM/vEzNb4/trSy4x1dZz39rjv38tsM3vbzFrUx2upwL+tHLjNOdcTGArcaGa9PM5U3yYBOV6HCJBngI+ccz2AUwiT92lm7YFbgCznXB8gErjS21Qn7WVg9FG33Q3McM5lAjN8P4eil/nv9/YJ0Mc51w9YDfxPfbyQCrwa51yec26R7/tCqgqgvbep6o+ZdQDGAi94naW+mVkz4Eyqds7EOVfqnNvnbap6FQU0MbMoIB7Y7nGek+KcmwPsOermi6ia7oXvrxc3aKh6cqz35pyb7pwr9/34FVAv445U4Mdx1PShcPE0cCdVw6nDTQaQD/yf7xDRC2YWFhN4nXPbgCeo2vUzD9jvnJvubaqAaCxTvn4MfFgfT6QCP4a6Th8KZmZ2AbDLObfQ6ywBEgUMBP7knBsAFBG6v4J/i+9Y8EVAZ6AdkGBmV3mbSurCzO6l6lDt5Pp4PhX4UY4zfSgcnA6MM7ONwBTgHDN71dtI9WorsNU5d+Q3pjeoKvRwcB6wwTmX75wrA94ChnmcKRD8mvIVqsxsAnAB8ANXT+dvq8CrOcH0oZDnnPsf51wH51w6VR+AzXTOhc0qzjm3A9hiZt19N50LrPQwUn3aDAw1s3jfv6PnEiYf0B4lbKd8mdlo4C5gnHPuUH09rwr8245MHzrHzJb4vsZ4HUr8djMw2cyygf7Aox7nqRe+3yreABYBy6j67zakr1o0s9eAL4HuZrbVN9nrMWCkma0BRvp+DjnHeW/PAk2BT3y98ud6eS1diSkiEpq0AhcRCVEqcBGREKUCFxEJUSpwEZEQpQIXEQlRKnCRE/DtcHhDtZ9HhONOjhKaVOAiJ9YCuKHGe4l4QAUuYcPM0n17Lr/g2zd7spmdZ2Zf+PaYHuLbc/od377MX5lZP99jH/Dt4/yZma03s1t8T/sY0MV38cXjvtsSq+07Ptl3daRIg4vyOoBIPesKXAZMBOYD3wfOAMYB9wBbgMXOuYvN7Bzgb1RdtQnQAzibqivmVpnZn6jaEKuPc64/VB1CoWqXyt5Uben6BVVX8M5tiDcnUp1W4BJuNjjnljnnKoEVVA0IcFRdgp5OVZn/HcA5NxNoZWbNfY+d5pwrcc4VULWRUpvjvMY859xW32ss8T2vSINTgUu4Kan2fWW1nyup+o3zWIc7juwnUf2xFRz/N1R/7ycSUCpwaWzmAD+Abw6HFNSw53shVYdURIKOVg7S2DxA1dSebOAQ/9m+9Jicc7t9H4Iup2qKyrTARxTxj3YjFBEJUTqEIiISolTgIiIhSgUuIhKiVOAiIiFKBS4iEqJU4CIiIUoFLiISolTgIiIh6v8Bpem+BwiBxdoAAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"climatology.plot()"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"5. Select 1979-2013 time frame"
]
},
{
"cell_type": "code",
"execution_count": 14,
"metadata": {},
"outputs": [],
"source": [
"historical = area_weighted.sel(time=slice('1979', '2013'))"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"6. Generate anomalies"
]
},
{
"cell_type": "code",
"execution_count": 15,
"metadata": {},
"outputs": [],
"source": [
"ssta = historical.groupby('time.month') - climatology"
]
},
{
"cell_type": "code",
"execution_count": 16,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"[<matplotlib.lines.Line2D at 0x2ae1d1239a20>]"
]
},
"execution_count": 16,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAESCAYAAAD+GW7gAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nOy9d5gdV30+/p4pt2xTl2w1S+4IF2yEDRhDcKE4gEmABJIQCMXxA6TBN2ACAUKCISFfSgIJPwMhEAgllC9OMCaxMdjguMjGXbYly7YkW2W10tbbppzfHzOfM2dmzsyd2/be3Z33efRo997ZO3Nnzvmc97yfxjjnyJEjR44cix9avy8gR44cOXLMD3KDnyNHjhxLBLnBz5EjR44lgtzg58iRI8cSQW7wc+TIkWOJIDf4OXLkyLFEYPT7AtKwevVqvmXLln5fRo4cOXIsGNx1111HOOdrVO8NtMHfsmULduzY0e/LyJEjR44FA8bYk0nv5ZJOjhw5ciwR5AY/R44cOZYIcoOfI0eOHEsEucHPkSNHjiWC3ODnyJEjxxJBbvBz5MiRY4kgN/g5cuTIsUSQG/wcOXLkAPDhHz6A93733n5fRk+RG/wcOXLkAPDV/30S39mxv9+X0VPkBj9Hjhw5lghyg58jR44cSwS5wc+RI8eSh+W4/b6EeUFu8HPkyLHkMVmx+n0J84Lc4OfIkWPJ41ilIX7+zX/6Jaaqi3MByA1+jhw5ljyOzgUG/+69k7jp4cN9vJreITf4OXLkWPI4Jhl8ACgai9M0Ls5vlSNHjhwt4GglbPALucHPkSNHjsWJ6aod+j03+ClgjL2MMfYIY2w3Y+yqhGN+jTF2D2PsQcbYz7tx3hw5VNh/rIKP/eghOC7v96XkWCCIhmUu1qHTscFnjOkAPg/g5QC2AXgDY2xb5JjlAP4JwKs4588E8LpOz5sjRxLe97378MVbHsfde4/1+1JyLBBEDb69SOPyu8HwzwOwm3O+h3PeAPAtAJdHjvkdAN/nnO8FAM754nSB5xgIMDAAQKXh9PlKciwUWA5P/X2xoBsGfwOAfdLv+/3XZJwKYAVj7GeMsbsYY7+f9GGMsSsYYzsYYzvGx8e7cHk5lhrKBR0AUG3YTY7MkcNDlNHbbs7wk8AUr0WXRwPAswH8OoCXAvhLxtipqg/jnF/DOd/OOd++Zs2aLlxejqWGsukbfCtn+J1g70QFtSVyD+OSTs7wk7AfwCbp940AnlYccz3nfI5zfgTAzQDO7sK5c+SIYchn+HP1pWGsegHH5Xj5Z2/G1297st+XMi+w3KikkzP8JNwJ4BTG2FbGWAHA6wFcGznmhwAuZIwZjLEhAOcD2NmFc+fIEQNJOos1PX4+UGnYmGs4GJ+p9/tS5gWW7ULXArHCXqRhOkanH8A5txlj7wLwEwA6gH/hnD/IGLvSf/8LnPOdjLHrAdwHwAXwJc75A52eO0cOFTTmTdzp3OC3DZLD5paIH8R2OUZLhiiitlijdDo2+ADAOb8OwHWR174Q+f2TAD7ZjfMtBRyaruHxI3N47omr+n0pCw4N25usOcNvH7WGdw8rS0QWsxw3ZPDzKJ0c84qXfPpmvP6a2/p9GQsSpL/mBr99VCyP2S8Vhm85rnD2A3mUTo55Rm6s2kfO8DtH1c9hWCq5DLbDYeoa7vjAxQByhp+jT+B8cQ68XqKRM/yOQQZ/rr5EGL7LYegaVg0XAeRhmTnaAOccn/jxw9h3tJJ4zM4D0/jdL92WGO88X9EC37lzH3YfnpmXc/UaJOksFWPVC5DTdqkwfMt2YWoMusbAWC7p5GgDjxyawRd+/hje+e93Jx7z4R8+iF/unsA9+yaV75M80UtwzvHe792Hy/7hF6g2nAWfbEP3bLFuy+cDSy9Kx4Wpe+bQ1LRFO3Zyg99DkBpTt5KNtuY/ATeByc9HAggN7obt4pkfvh7nX31jz8/ZS9D3WazJM/MBYvbVpcLwHQ5D98J5DZ0t2rDM3OD3EH44OHis0kQASvZwJK1eHmyNeRh4cgkCly987ZsY/mJNnpkP0C5vqWQrW07A8A2NLdqxkxv8HoISgNLGju5TfHmAyQZ3PraWi43F0SJptSCHHZqu5fXzJdCYqFrOkrgvXpSON19NXYPluGjY7qILmsgNfg9Bidppg8YfYyFJZ1I2+POg4RPDl1PLFzKI4WfdHR2ba+D8q2/Ex6/Lq30QZGftUihCZ7kuDGL4OsP+Y1Wc+sEf45t37Isde9eTR/GtO/bO9yV2BbnB7yHIhqdxBCHpyAZf6q85Hzp0xXfM6WxxGHy6Z1m35bN+NM+PHzjYs2taaJAd95UlEO1kOS4KQtLR8NCBaQDAdfcfiB37mn/+X1z1/fvn9fq6hdzg9xAitCvF7pDsIxsnub9mfR4YPk1umREnOZEXAsjgOy7P9D1cfwc2H/6ShQKZ1c8tMslPBdvhMDSSdJiow0SF+BYLcoPfQ2Sw94Lh1+1gUslGfj4YfrURP0dlAW/j5VBWK0M8tZCA5mFxXSiQJZ2lkM/gRemQpKOJOVgyc4OfIyOI4adp+BoZfEtt5OfDaVtRxFov5EnekO5ZlvtXzw1+DDLDXwrJV56k44dlSr6ssrm4TOTi+jYDBtLlUxk+I4YfGBs5y29eGL6Czc/UFrDBl3ZLWZzetLuSd1lLHbWGbPAX7ljICtsJnLYUngkgVFBtMSA3+D0E6fJpkV2GQtKRWel86MqqzNrZBczwLYdj2Ndes0g6tNguYLdF11FpOBgtGeLnxQ6vlk6QeEUoF5IryC/EkM2uGHzG2MsYY48wxnYzxq5KOe45jDGHMfbabpx30EEOQ7dFSUcu3DQfYZmqCT2bwPAbtjvwGawNx8VQ0ZuorUg6OQLUbQerhgsAFra8lxVylI6pBWaxaCSbyIVYfqFjg88Y0wF8HsDLAWwD8AbG2LaE4/4WXmesJYEsDJ+QJOnMd6YtIYnhn/Hhn+B5H/9pry+pbTguh+MGDD9Liny/tPtdh2b6bkyPzNbx8MHp2OsNx8XyIc/gL3aG77gcnHvhmECY4aex+IVYYK0bDP88ALs553s45w0A3wJwueK4PwLwPQCHu3DOBYEsGYpkkGRJpzHPUTo1FcNPMEQNx8WR2cHtc0r3a6hghH5PQz8YvuW4uPTTN+OPvvmreT+3jEs/9XO87DO3xF5v2C5WDJkAggJq+49VFk1FVRk0RgJJJzCLabkcS5LhA9gAQE5H2++/JsAY2wDgNwCE2h4udtgZJB3LPybM8GVJZz6idOIGv9/Ms13Qjmi46DH8Rob7Jy+w85V/QHHedz15rGfnePzIHF7+2VtwbK6ReMyxirpuUsN2MVoyobGgzeEL/vYmXPKpm3tyrf0EGfxA0snWzHwhFljrhsFXpWdG79JnALyPc950b8gYu4IxtoMxtmN8fLwLl9c/tMLw9x6tiPol/SyeRhh0nT4J5PMghp9l2y3vrrpdDvj+/VP43E93xV6n8hnDPUzs+fxNu7HzwDR+8mDrGcR120XR0DBcMBZ9iWTymamctmmNUBZigbVuGPz9ADZJv28E8HTkmO0AvsUYewLAawH8E2Ps1aoP45xfwznfzjnfvmbNmi5cXv/gZNDwaUD97JFxvP1rOwCEt4rzk3gVGLwh0r4X4GAGgntH3yPL/ZMZfrejk177hVvx9//9aCzkkwrkkXO5F6Bxp2UomRHd2TRsFwVDw1BRX/SNzCmSS068IjgphGEhkqJuGPw7AZzCGNvKGCsAeD2Aa+UDOOdbOedbOOdbAHwXwDs45/+vC+ceaIjEq5RIfNmw/mL3EQDhgTQfDkWZ4VMo3kLcrgIqDb+1KJ2k6KR2QWG3R2bDssrUPDB8Me4ylEiK7iTJ4C80hu+4HNc/cKClkEkaI5R4lV3SWXikqGODzzm3AbwLXvTNTgDf4Zw/yBi7kjF2Zaefv5BB2n0aWZYlB5ldk6GYF6etZPBHWghnHETQBC0XvKGdyWlrqf0n3QBFuozPhB3dpOEP95DhowWGH83FqDsBw4+Wz/77nzwysCW1r7l5D678+t0tyVhEbihKRyYAabLsQozS6cpo45xfB+C6yGtKBy3n/M3dOOdCADGANLIhG9aVftyz5bgomzpm6naoTECvIF/DcNGAxrL5HwYRNHkpQzKTpOMExqvb33v5kImnJqs4PF0LvS4knZTEnk5BhCNLDVTZyHHO0bBdFA0dQwqG/7mbdqNgaPjji0/p5uV2BTv9KpethJJaEQ3/V3uDdqNRAiBLXwuRFOWZtj1EYDyyefrJ4Nt+u7WC34ih15CZSsnUYehapgzVQQRNwnIrkk5CHaNuYLkf2jgeCWWdqhDD76Wk46HVXQ7JO57TVlcaz0GtkR8spNnvK+1uKMlqrOyNnbWjxZi0KUtfS1LSyZGMLIlXMoNYMRQwfEPXUDC0ljJtp2vttSaUB27J1GFqbEEOZiBYvFpj+Nm28O1grOQZ/MPTYYNPUTq97EAgeiq3UE8ICPxGBV3DUNHAXN2OaeKD2jmB7muzBYlzjq/e+gSOzTXw2PgsAGDL6mEAwL/+wXn49hXPxUjRiDH8UCXbBUiKcoPfQwQafrbkDWJ7lsNR0DWYOsvMOB8+OI2zPvLf+OE9T7V8nZY0qAu6BkPXcO++SVz/QLz5Qyv45h17cc++yeYHdhHRKJ0sC1cvNXx69DGG7xsmq4fSGX2yqlZSFLIhEwbfZ/hzdSfm1M3iF+gHpvzmQc0knXv2TeLD1z6ID/6/B/DIwRkYGsOJq0cAAOuXl3H+iatg6CxGAOSFcSGSotzg9xBZBoSt2CLargtDZzB1LXMc/pMTFQDAD37VusGXr8HUGUydYceTx3Dl1+9u+bNkvP/79+PVn/9lR5/RKmiCEsPPcv96yfBpwY4yfDL4vYyGIqJRs1pj+GT8i4aGYZ/h1yI9EwaV3dJ9beZUPuYvDMcqDTxycAYnrhlGIVI3R9e0GAGQo+YWYiRbbvB7iCzlkR2X4/knrcJoKdg+UvcdU9cyZYoCEA2Yo9EgWSAvTIauiWiFTtCvSoI0CUstxOGHWFuXDT4tJkfnws9lxpffeukcJzmwlqHss7woyAx/tGRitmFjNuK4HcTy2ZxzYfCbMfyDU97zKBoaHj08g9OOG4sdY2hxhh9urrPwGH4PY8JyiNIKafU4XBfrxkqoNBxhnCzHhalrcHn2TNtZPznmcBsGX2ZrhsZCzcxdl4uKnq2gXxUoaRIOmdklnV6yNksY/HAcPhmkXkZ6kI6dTdKRNHwnMPhjJQOcA0dmogvW4Bn86ZotQqCbGfynJ6sAvI5zh6brWL+8FDtG11ic4Yd2gznDzyGBttRpU5rYvMwmyOAX9OxOW6p90zHD15jYLQDhUgMy22nG4PtVi0eEZbbE8F0RodFthk8GfSJi8On+9DKWmwx9JqetiuHrmkjEixKJmTYDBHoJeVGtNkkWe8o3+E9N1tCwXSwvF2LHeHMynpBGyMMyc4QgDGmTOHxD12DoTNLwvbDMoqll7sIkG9gsckq14Yiql7IhN3QWSi2fq8tyR3bn5lyf0vGDxKvAAd4MddsVTt5eafgzNTtkLAaP4bvSz36YoqljpOhFGR2K5BFQ4lgSGraL0//yx/jBr/a3dM2dQDb4zRqvk8GnCB0Kn5Whayz2fMK7wdzg55BA7MBpUlPb9B20JK1YjgtT8xh+VklHNrBZkk7e//37sP1vbsB/3vt0yCAYmhbq6fnUZAXX3X/A/z7B92g22PuVjk/X1UpYpu24oll11zV8yUBMViSD5C/QvdTwyXFZT3DaysRA5bRNY/h3753ER659MPHcR+caqFkurr7u4fYuvg2EGX76HHj8yByA4PksL8cNvqlrqRr+Qsy0zQ1+DyE7YROPcbgwso7stCWGnyHCAggb2Cxyyk2PeJVI/+ibvwrJDbrvLCa85V934B3fuBv7j1XCRd2aDPa+STpuoD9rLJvBtxwuFohu67KW44LWT7rPnHOxKPcy0oMcsUkMX36eSWGZZPDHZ8IMHwD+9dYnEs9N47Ggz5+JoTLQK4cLqX14D8/UMD5Tx+qRonhtWQLDT9Pwc0knRwjE7NOMI4Vg6pomBpCs4Wd32trKn5Oga0y5jTV1FioPS1EPNcttkeH3R9Khe2hqXj5BlklpSQxfdTznHB+59kE88NRUW9ezbsxzCBIDbTiuMCTzIukkyILy2JIXhYYUljkaSRy78JTVmc5NTl15LPUatKBuXFFO3eXuPOA1cXnRqUE1Xkp6lKHS8OVxn4dl5gjBkWrpJEXqBCGYTAwgy+EwdYaioWdn+HWZ4acb28lKA0fnGnjNuRtj7+kRSYfAOQ9r+E0Ge6XPTttWSlPYLkfJ9KaCSmKZqdv411ufwOuvua3l67EcN2bw5XLD/ZR0ZHlCPqYeMvgewz/kM/xPvvZsXPmikwAAZ25Ylnhu0vhVY6lXODpXR8nUsHK4kJpp+9DTXr2dV5x9vHgtScOPEptQc6IFGJaZG/weIjw44pPOM6LktA2SPGzXhaF5pRWya/iBgZ2ppzvUSL887bjR2HumzpRx+LbLQ4P/vKtvjIUayuh2XfmsoHto+DuVLCzMclzh5FVp+CR1t2OcPYPvSQeU7CPLb50mMP3vYxO4+rqdMULBOW/O8O24kZdflyWd/cc8J+dIycCfv/Q0nLlhWWqDbyrzYc6jpHN0zsKq4SKG/Po/SQEPt+2ZwImrh/HsE1aI15RROopM21ZIzyAiN/g9hFxSQSWBkHEx/bBM23Xxxi/fjkcPzXoavqGhnrFI1WzdFiy1GcPfe9TLyj11XWDwX3yat7295BnrlNvwhu3GBv8jB5P7m/ar8TVNQlPT/EzlLNnOHCXD1/AVk5g+M835noSG7Yp6OsS46fkUdK3jSI/P3vgorrl5D75w82P44s17hCNWNuBJmbYNRWQOEI7DL5s6dI1hsmJh1XABI0VDyIFpJUOmq/Mv6Rydq2PFsImyaWD34Vmc9sHr8eP7w+VBZus2/vexCVx0+lrxXACIuSNDlWkblnRyhp9Dgt1E8w5aq3kyiu1w3LLLa4JS8IunZU1gmqs7Qjpo5jClZhwblpfFa+dsXoEnPvHrOHvTciUrq9tubPCbKZO53wxf1xlMLVstIqtJlA691h7D56LmPRlYYvhjZbNjlrhxxRAA4O+ufwQfu24n7vZL+4YNfnMNX8nwdQ2MMdEjYatfXAwAGGNIs3fE8LuRtZ0Vs3Ubo0UT8vC9L+J3uWfvJBqOixedFu6mxxS1gYiEyWi2ax90dOVpMMZexhh7hDG2mzF2leL932WM3ef/u5UxdnY3zjvocJpEtdBrph/7LjvwDF/Dz9rxaq5hC4PfzNgenatD1xhWDQfbWJmJqXTXuu3EjJOeos9WEhK2orh//xT+9vqHu1aKQXbamoaWWdJJM/i0aLQr6USTwEjDX1aOV2NsFVFjTs9EZuyZGL4VZ/tF/55sWukRA9ng6yw9g7xZnH4vULO8e326VCaBJCkCzQ0qRf7Vt5yHv371GcrPMzQWmsOAuvbVQkLHBp8xpgP4PICXA9gG4A2MsW2Rwx4H8CLO+VkA/hrANZ2edyEgM8Mnp61c4qBlhm9nZvhH5yysGDJDJRNMiYmpGH7NijP8tIqJsqyUlvjzys/9Av/8s8e6VopBdtqaGaN0bIeLDlkqo97uxHZdz0cjdmuOguErzlezHDz09HSqQZWPVYEM+GjRSNSyZYZfVUTpUEjlRaetBRB+3rqizowMYvg1y8Hteyaw1y/u10tULQdlU8ebn78Fuz72cgDxFqFiMfMlvBedugZvfO4Jys8z9HhYZnhOL02Gfx6A3ZzzPZzzBoBvAbhcPoBzfivn/Jj/623wGp0vesgap0paIAOv65o3gSTD0vDT/RuOm4n9ztRsrBkpgrHmDP/YXEMwHEKI4SukmrrtpDqwopAXnSzGvN2Wef9579N44d/dJIyjcNr6fpFMko7rijh8ta+lvYltSTkBRV0Txod2P8vLpvL6vn3nPlz2D7fg//zHvU3PEY1GiZZTGC4aiQuWKvOX/paxQLJ7/XmbsXK4gN8+b5M4RmMsk4ZfaTj47Wtuwws/eVPT79Ipqg0HJVOH5ueSFPQ4YZJDTptBtajJRn6pRulsALBP+n2//1oS3grgx10478AjxAZSmKOphTNtAa9WSdGM99hUoWY5qDQcrBopYLhgZJB0GrG4Y1nGUUk1dcuNGae0Sp5yM5YkFirX/am02UHpr/7zQew9WsGTviPai3BiYIx5DWQyJl4VDA2MqROv5Gd3yad+nvnahLyke9ciNHwh6ZhKY0z35SG/XV8aonIN3Ws610jJSLwH8uvyAl2zHJQMXeja65eXcfdfXopzNwdRLVkZ/nx2xqpZjtipAfCDHiJjVnJIN4OhctrmDF/Z/EY5EhhjL4Zn8N+X+GGMXcEY28EY2zE+Pt6Fy+sfZOOhmhwi3tkMnLaEmZotttTNQjOPShmGI3798tTjKwHDp126XD9HreHHo3TSmO+BqSAzM8ngU7QQ0D7DJ135wac955ztcLFgGYpaKFFwzuG4QbazirXJz2X34dnME50K35m6HzHk/04L8ljZVI4LYtuUSGQ78XtPiN43GlMkXQz7XZtUu0RRVmDIDPlcapYryEYSNI2lRi3N+olX85mPUfUXKoIX1hy+P1G5Kg1Khu8Gi/hSzbTdD2CT9PtGAE9HD2KMnQXgSwAu55xPJH0Y5/wazvl2zvn2NWvWJB22ICAbCtWEo8laNnU/IiBs8Mlp1iz5Sjb4w0VdMMit7/8Rrvy3u2LHy5IO6bKykVfN47rtpIaoRfH0ZE0ksyQ5DWVj1a7B3+RHqVAyjZe05g1rM0PiFU3agqElstboZxycjpcZUH62Gxh8Oadi/7EKxkoGlpVNpTOfjO/RuQZcl+PkD/wYl332FuU5apYTMl5RSWekmJxBTIvCyqFCKDO6bocNpwoaY6mtO+siIml+GD7lHZSlXrYqhi+TrGZQSYKi34KpL9laOncCOIUxtpUxVgDwegDXygcwxjYD+D6AN3LOH+3CORcEZI1TRdBqUjSEEWEcMzULRZ0knfRJIxv80ZIpyiFwDlz/4MHwNbkcxyqywfdelx21qnlct92YgU8ypnXbq8RJ7Dvp+uXtflrtkzQQy3zYzwlw/FIVADJJOvS+oTGYmjouPrrQ7TtazXRtYjGh/sT+uR4/MocT14zA0DRwHt/9EcN3XI73fe8+AMAjh9Q5D1XLCdVyJwMnNHy/mbvKONH9Xz1SjEg6rjIuXYbO0qOWslZ5zYK//q+H8P2706tuer4uiGgrAMrExVYYvjrxyvu9ZOpLM0qHc24DeBeAnwDYCeA7nPMHGWNXMsau9A/7EIBVAP6JMXYPY2xHp+ddCJANhWpy1GSGH3GUjpVNwUKahWZSBueKoQKOX1bC01Nqg+S6HOOzdbgcWD5Eko7P8PUmDN9yU2OSZRz05RzqEfqn375HvcORDX6bWi/tDKg+u+XLMwBiuyYV5FwIXY/XTgHiC9v+Y9kiToSkY3hlHug57hmfw4mrh8U9j36+7ED9j7vSDV3NcrB2TDL4NpVT8P4f8cMSVQyfzrN6tBDX8M0mDL+Jht/NBjhf/sXjePd30h3Y1IKxbMoMP16apGF7xeyiBEsFSrzikQRKXaOyHQvP4Hel4xXn/DoA10Ve+4L089sAvK0b51pIkCeEKqKBGH7J1ENhkWtGi7jmjdvxq71eYFOzyTPhJ1KtGi5g44oyfvrwYSWz/fiPd+KLtzwOwHMYAoAuJB2Z4au3/9EJnsSen570DP6WVZ7c8uRExUuKKYXrlchNKmptbv1p0aCMWttxRXSJrJsnQUTS6CxxgYgyOSoz0Ax0f0wpxLbSsHFgqoYT1wwLGS3O8LPvdqqWg80rh3DH40cBBPJZIOmQwVcw/IbE8BsOOOdgjKEmNYRJgt4kSidrDSgVpmsWDk3VcMq6eOmPJNA4kCUdJcN33EwOWyCQOV3u7WgAb7zoGmup7MkgIc+07SHCHaLi71clViJHxrz+OZtw3LJSSwxfY54R37C8jLrthhyihG/cvlf8TIaATttMw1fF4SdtaWnHcYFUWVGl48u6fbulGOgz6B7JTttsGj7F7Xsa/jdu34urr9sZOia6s6H2eM3QkA2+v/jQczlh1bBgmdH7WGk4ov5OGlyXo2a52LC8jDv+4mIAiigd/znTOQ5P18SCQvdu1XARjssDh6/lCP9REppF6dRsB2Ol9vjkW75yJy799M0tJeORwZelqKIRbyBUtxwRg98MNI7e+tU7xWuOw2FqXtmTLI1lBg25we8h7GYMXxqkcpkC0tMLuu+0bWLwKcxS05hItVfVuZGNOk1GTSnpeNcq51V5mbZqhr/vaEV0DgKC7z1WMvDJ154V+q4yKt2QdCzqHOVdi+XKTtsWJB0tKBp3zc17lMcQspaNiGr4DcfFVMWTnlYNF8QzjzpuKw0bm1cONf18Ghflgo61YyXoGpOidII4fO9avN/Pu/pG/OY/3eqdx3Jg6gzLyt4xJOvUbDeTpNOM4UdzPbJix5PeznaqamU2+nIABKGQEJbZKsP/2SNBtKDteoSiXNC7YvBf94Vb8dyrb8QP73kKL/n0z7uWcZ6E3OD3EDIDUoWwiW1oxGlLBqsVhk8RMRv9NPiHFQZfdsySvMJSnLby5PFq6YSvg0IYL/y7m3Dx/w3i0ykJSmMsiDRSfIdaKEqnPadtNcJoHT8On75Ts57Aclx2UqGv6KKR3eCHJR3LcTHthyuOlkyxwKictuulOkdJEKzWN2AliXUSs1VJOjQ2qg0vM5UWBdpl1S1HfGYSdJbM8DnnqNsONqxo/h1UIIfq+ExdKZvctmcCn79pd+g1WR4lFBWyS912MzdlUc1Zr0OdhpLRHYN/5xPHcHC6hnd/5148emgWMz0OY80Nfg9hu1wYVNXKLRh+QQ+xb2J+hYxROpMVSzhhN68cwmjRwD/cuCt2nGzQqMYIlVdQSTqliMFPyzoMvS6yXTVhOFSTo2o5QismeatV0KIhGL7fIxjwZJpm2ZABw9eUCSXyZxOouUczBHH4gdOWasyMlQ1xz6OfX2+ZAEoAACAASURBVG04GCkauOHdLwy9Hh1DtYhuXTR11GwHUxVLwfDVYcHlQmDwqeRDLYOkozF15BngPX+Xe7KVjKy1iIb9UNLx2bqSKLz+mtvwyZ88EnqtlpXhZ/BPEB4fn4u9Rt3oPIbfPQ1/2H+G45FWkt1GbvB7CFeSF5RhmYKhhaN0ogy/maQzVbWEE3aoYOAjr3qm8rgww/cmuXDaKhi+zPLqVlzSSdLwBcPXkMrwK75hKxgaKlZ3GL7stC3ozUsrBBo+S9xJRb/nzgPT+PT/PNp0+y00fEMTmbaUgTpWMsUzj37+XMPGcNEQkl708wiBbu0dVzI03Ld/Cmd/9L/xzTs8f00Qh6+4/5aDoYIRGHyfXdZttynD1zSWWOuHnvXGCMPPyoiHCtRWsR4z2LJBlM+vctoWDV3N8DMafPmzCJbfkrRkal3NIqad2JHc4C9cUOEsQF1ZsGa5fg/ZcNORQMNPl3Q45/jG7U9i/7GqMPhAuM49odKwQ5o8hetRWKap0PCjDD8qbUS1Z1FV0v97XWPCcKjq+lOiTNnU2068CqJ0fIPvSpm2GSQduua0dpJRKatuu/jsjbtC2cTqz45o+LYrasyMlozAaRvJyK5ZXm2fqGGK3qOobl0yddy338s4fnKiIrqmyd8z+vdlUxfscrbu+I7g5mGZOkvOtKVnTTkA4nwZDSQx/COzjdAiwTnHu79zj/LzoosfQAw/nmmbleG/92Wn4+yNy8S5gSDPo2R2R9IhDPkGf3w2N/gLFo4bsM0kDb9keDXHZYNLzC9gx+qBdf9TU/jADx4IMXwAWKuI8Nj2oZ+EEobIEBCx15tKOo6IUb902zoAcWZKcgVt3fUmGj4ZHOpQ1CooSgUIFkXLcUWIq9lE0rl//xTu8h2Ehs4Sd1KyHLJC0QovCbKGb/oLynTNwnBBFz0QgLCPgAzXUEGP9RuI3qNaxMhFF4iioQvyoHJeVy07JOm86V/uwCkf/HG2xKuUKB26j9HPyLqo61qg4cvP5PBMHbfsOiKeQSXkA4pLOioNv9ECwx8pGnjpGceFvpOX59G5wT84VcOWq34kfqdFN2f4CxiVhiMmkzrUMUgF1yWGT8y+JPRt9cCSDdGYZPBXRaIj0mQNjQUOTgLF4csTtm674nwUeRPV8Cd9g0/RG5rGxGckafjlgo5yQW9re0yOOrl3reNysWCaTSSdV37uF/ibH3khmIaWXIpa/p5y0blmUltg8P3uZbaLmZolnpUw+NJzpJDJoUKc4UcNflTGiLJy2RFt2W5ol+m6HNWGg6GCHmLijstRs5uHLqqidDjn+Jv/egg3PXIYAGKfkdVA0nGewQ/+ZtchLxKMirjJ8yK6+AGdR+kAECUm6HMcX9LpZFcKAPfunwz9TteUM/wFDNKoAXVYZlWKCVZp+CS7JNUjkQeczPCjWYRpzSjSaunIk6faCBKvDN0LI42y58kKFfsKQh3p+6n6qlYkht/O5KG/GSubcLlnmGWnramrSxd43zH8WsFI0fClv3/OlpXi52bO9IYdMHwRpVO1hf+EnrO8KNF3GioYTSUdkpTWjno7uiijLhpacA6Xh3aZsw0bFb+c8FAxbJi9EgWtR+n8754JfOkXj+MDP3hAnP/XpM5SSYs65xw3PHRI3AfyJRydCzP8XYe96KKT13oZ3LLfh3Z6UQ2/HtPws8fhA8EcoPFr+5JO2dRRs7OVLldBj/SSIKnvyExyn+huIDf4PcRcPZjcqt1v3Qq6IcmZtkLSMbxtfVIY4IxUglg2+FFMpDQbF9UypfNTSODqkUAamqxYoVrzqpDHyUoyw1dlXtYkDb+dWjpkQMbKQSSK7bowtfCuRcXyoxEWaa346Hvf/ZeX4qxNy8TrzbJJ5cJsIkqnZoleqsJpKw0OKnw3VAhnXwPxDNx9RyvQWPC8oobMM/gBw5cN9FTFQtXyGD6REhnZ4vDDC+d37twnzgt4QQdf/P3t+P/e+GwAyTvVW3Ydwdu+tgOf++nu0HFzDSd0j3cd9hj+Sb7BDxXfi4SoAhB+E/kaGy2EZcrfhXYQntPWG9eOy9surxAdbpSsmDP8BQrbcVG3XRHvrnbaOsIgygxfHpDDKeWO5fDAtKzGNF2QGL5cTuHPX3oaPvc75+AFUqbssUoj6CalqcsQkMEPafjE8FWSji8plAsGqm2EuJHEQYtdwy/wFmTaev9PV63Y+aPGM9rlK1w/JchalY1js9T6aBy+y70kOZJ0ihG5AAhCI4eLRqgjGeD5bGTsPVrB+uVlce1xhi9r+BGDX7XE/S8aGiKnyhSWCYSlSiIWoiKlf/41/g4kieFT8b89R+bAORdJeDXLCe2idh+eRUHXRIXUqMGnVqHB949HubWi4QMSw7dkyVCLMf9WEe2hSwZ/NmPIb7vIDX6PQDJMwPDVkg45mcJx+JLBT2loMp3C8OWQuChreOGpwTabTivb7pKp4xVnrQ9tO+u2K5JCdGL4CRp+KEonJbRUSAqm3lbi1VN+iQPKSm04XiSRXB4Z8LJL5cQwOreMqINUZm5yIxPZ4Ddn+IGGT9cyMdcQi7PKv0GL5nKFc/hvfrRTfGfAM/hyRm6M4UsZ3N7uJ27wS6bX6GQ4wvKz1NIBwsEIUQZPn1E2kxd9IMgFcf3yDrQwVRpOaNw8dngWxy8viSieqNM2uiuh88sLc+sGP8rwXeG0BdqvARXN+aAx1u4CkhW5we8RiJWnSTpy+Js8CGW2n9bQRGb4Q5EQuOv+5EL8+UtPA+CFtxH+8EUn4mtvOU/8/heXPQPLyiaOX1ZCFNFt59OTVRR0L6rI0L2GLTITPuIvLKFMWyPMkGTU/AWv3GaUzj6/Ls3Ja7wtfsNxQwxavo9PRerfRM8X9XvIbNT2C2YxxmKRS2loRBg+4N0j2vWpWCKVtl5eVpclOChVQt13tBoy+NHM1pCG74SdtjM124/D964hGkKZRdIBwv6RqLGiZy+S6xIMfpAlzWM9EuRFYmKugXWjJXHNcjmOmkSegvPH5cR6C2GZgMzwHXGNpOF7r7eXfJVUlqLX9Xlyg98jkGQwUvQlHSXDD2qWyJUkQwxfamgShWzwV4+GDcRYycQWP9PxiMTwo5PiJc88Dvd++CXKCU5yD02wXYdmsXrEO4+heS0ZZdZIFRtlrZ/yDKLGsWY5mKpaWDFU8KJ02jD4T05UUDI1oWFPzNZxdK6Bras9I6hqxk6Yi0k6Yc4lTzzb122BMFtsGqVj084gMPicB6GdJcViSI7vJJ8MVUatWV7PgQ1SCYbztq4MHTtSlGL9Iwx/rm7DcbkYD8MRx210nERBspkbYfjyfaTEQfruVkJLTPJF2S4Xz6VgeIlN0Xu8bMhE2V+c5F1htPkJECw4nTD8qCxE5bdpvrQTXdawXfEco+h1S8glbfAbtouXfeZm3LKr+60UyUgLhq+Kg27YYmItlyZ4VMNPc9oeN1bCjg9eguOXxeuW0GCVNfxWBjvpjHRtjxyawRq/9rrpM3xiZ2VTx6/2HsMxv0sTELBAr+5IeOLes28StstxzublnqTTxkB/csKTNMiwPOqH7dFCl+acq9Sjkk74WNngy120zlgvOW2bMHzL8XYGusZEMxsAWOaHdqoknemqBcaCcUP4zG8/C0Cgk5Pmu3o0cKxvP2FF6G/O3LBMGOCGE9fwvWsggx+c7zlbVuAFJ69GGoSkIzN8y8VJ/m4LCMYfjbloxAyBdGvHdYOSzcMFVBp23OCXTTFn5F2aiuEXjPj9rbcalhlh+JTJTb1z22Hkb/3qnfjz796nfK+b5RpU6IrBZ4y9jDH2CGNsN2PsKsX7jDH2D/779zHGzu3GeVsB5xw/+NX+EJM8NF3DwwdncNX37u/6+UiGGUuRdGalKJ6xUFhlWNJJNvg2xspGKJpGBhlCmeEnpcOrQPq+fG1r/HMZugbbdYXBv+DkVXC5V5jL4TyUyFU0tdh2f8cT3m7g2Ses8MIyLaflELd9RyvYvHJYGONH/a5QW/xOW0nF0IC40zbax1cl6QBer4Jbr7oIQDYNX5R5kIwMLaBFhbY9WfWieKIO23X+QksOTvpfzgsYLZn40Cu2iXOesWGZiPSxHTekt9P4pO9Nks6HX7kN/3Hl85VlBWQEunvwWs1y8ExpQSSGTQuvHNU1WWkIHxT5hubqjvB9rRopoma5MY18WdkMJJ2Q09aNOZrJWNNz4twjKMUWonTEoizJTl4GefsM/5ZdR5SvGxobfEmHMaYD+DyAlwPYBuANjLFtkcNeDuAU/98VAP650/O2ipt3HcGfffte/N//fiT2Xi+2UTRwR1KctjM1W7yvJzlt0zT8uhVrKiKDJpwcltmsXLCMS56xDq9+1np86reeJV6jLF6vuQgX22VaFDwmGY4zljsPXXf/Abzz3+/GjQ8fxvHLSlg+VECpoIPz1tnN4Zka1o0VhTF99NAMGAucuGmSTsxpmxLz7jmC5e+TrcZRQ8r6lQ0oOWRVDu2pqqV02I6WDAwXdCEFHJvzjGW0BPFbXrAV/+clnu/mWZuXi+9lORyO5Iim8an794gknaiWnwS6HSGnreWEMpFLEUlHllae9dH/wUV//zMAgaRzrNIQC/EqXzqcrHrfl+bHsrKJoqGBsUjiVcNBORKlRPe86sfr0/mbRSDJoDn0V9c+iOma5TltdQ2lQrojuh2sGC4MvsEHcB6A3ZzzPZzzBoBvAbg8cszlAL7GPdwGYDlj7PgunDszDvhOu2OVILKFJlqr+vFPHjyIQ00aWdPAHU3Q8C0/bHNEMcFkKaIZw49u/WWoJJ1W+nCWTB2fef05OP24oDYPMXyv1nzA8Ed9SaBhu3A5Dzl8ZYb/V//5IH503wH8au+kmJBDYoveWqRO3fZrzvj3a8/4HNaNlgSzi+rycmJVzOBrGq591wX4veduBhBefGzHDcXpk8Fo2k3LcYXBlbV2qmzqOcDjUTrLFfp9wdCwcqSAiTnvWdL/K4fjx779whNx30degrWjpaAipxsub033mnIWyOkfTcJKghbR8DnnocxxIM7wZecsEAQTUMnoyaoliqNRlNlkxZO46DOWlU0wxmIyYFUh6dDvVIm1lX62BBpLE3MN/MsvHoftl1ZoFnmUhLSKoSuHCt7C3AIpaxXdMPgbAOyTft/vv9bqMT3FbCRqBggGQCsM33E5/vDf7sKrP//LTOdLYvhzkfdlyFKE57S1lXLHxGxDaRwINFjlKJ1m1SNV0DQm2hVSTLWheRq+aKVXCgy+7fCQgSxJDF9eXMkgkLFp5TmQgSmagUN037FKqFNUlOHLOyUyeLSYGTrDWRuX4zXnbgSgcNoqGX4TDd8OdgYnrAqiaYjBM8ZiddUnq1ZIQiMUdA2rhotCyjmmkHQImsZEcpeI0rF5aAzOSiG2QKDhZ2X45NAnibDhuHC5N+ZOWuP7UPz75EU4BWMv2o2Ngg+mKhaenPDeO+24Me9+VCzB6IEgya5cMCKSTtxpW444VoXBb8VpK+0ahgq6GNvtOm0PJPSbBoAV/uLdS5bfjZ62KqE0ap2yHOMdyNgV8GQfbN68ubMrkyAMsBxH3UbMKz3gZpUSKzGnbfh9GuSqLMeopOP6coc8oGuWg6enqtiyemPiNajikH/ttLWp152Ei5+xDl/+xeOCHRl+HH7QSo8kHcdj+NITL5pBqzn5WkScNm29W9hpUc31kpRcxDlCDb2jBn+2bmOFL4GQsfjWFc/1K0tGrkV22kqx/YC32GksWy0d+jtZepMX6ZKphXYT01VL2e3KNDSsGi7gaX/cHfWZ73KFwZeh+9dqRyKqaHzS9VHxrqEm2r38uUAg6VAT8ZKp47tXPh87D06LYxgL+gFYjov3fz/ssCRJp+G4ePjgNNaNFYU0dKzSQNHQxWIRlAEP526o4vDJsRqtqNpOLR3AG+O2XxBxWOFHyIK9E/HWo4S1o97YrVlOLC+iW+gGw98PYJP0+0YAT7dxDACAc34N53w753z7mjVrVIe0BeGkkphasy25Clllh+gCE2X4qh0HISrpyMcT9ozPgfOgrogK8qB57okr8fjHL8PzTlqV6fqjeM9LTsWfXHwKXnG2p8RR+0CaRDLDJ8cWgRg+lf4Vr5vhxJxWJo/I5jS10P1aO5rcC1a+h3MNGwVDw/KhAs7etFy8ripJbTtuyKnL/PyCLBq+Sj5YFjL4EYZfaYiWgzJMnWHZkCnqIh2b83Z3uqbiUtG/9apG2iENP4HhZzQ00SgdkuzKpo4VwwU8/6RwlA81cb9tzwRu23NUvM45D9V6+tXeSWxeOSQWnmM+wyfIBl8eL3VbZfC970KOX9plthKHL8uCNcsRY5vKGVfqDm7fM4HDTSRewmRKXauz/FLMvQzN7IbBvxPAKYyxrYyxAoDXA7g2csy1AH7fj9Z5LoApzvmBLpw7M8hxWYskYbQKmYWmRZXM1m0MFYLWhdFDgwUhvn2PRukAiDluqYesHAYXhbyYjJbMWDp3KxgqGPizS08V8ouhaX6xsriGnxSlE00bDySd1g2+XB2xYATnWicxfHq+J/pRO3QPp6oWvrtjvzLWXDgYpbEhF2STv5Oqxr8MmeHLkD+r5BfhAjx5ZKpqKZOuirpX1fLQdA2Xfurn+Pmj42K30gymrsF2wpIO3QsyaOS0zcrwo1E6NC+Siq5RRdOD/g6FjNvhmTqemKjgjA2ehHNgqobNK4cD7Xy2jnJBFxIBGfxohVUqtS2jHPENtcPw5TlTtRwRokufPVu38aav3BHrg5yENLmG5NJehmZ2bPA55zaAdwH4CYCdAL7DOX+QMXYlY+xK/7DrAOwBsBvAFwG8o9Pztgo5YYXQDsOXB1laoaMjs3WsHikKaSNaD5+MX3MNX83wHxufBWPA1tXhNnIySpJDc7TLW0Sv9DCXJB3v8+t+GV6NxRm+XArCu76wjNKKdhnUa9FCnaFkhn/hyavx8d88Ex/7jTMBBOF/X/nl45iYa4hYdBmC4UvSk+26sbBNKnecBsvhMKXF6Np3XSBKS8ufQ997tmHD5eqyCqbBPA3Z5dh1eBZ7j1YSw3Fjf+uXiVYVaaOy3OvGSjB1pvQJqEBrlpB0JIavAhUyoznz2md7UuTPHx2H43K8dNtx4tgTVg0JYnF4ph7aEY0pGD7nPN1pG+mZ0IrTFghCV6sNR+z2qGzI/mNV1CwXhzPWsU8bM9GY/16gK3H4nPPrOOencs5P4px/zH/tC5zzL/g/c875O/33z+Sc7+jGeVsBxaLLur1887Pq+TILlRuKRDE+U8ea0aLY+kYlHTJ+sob/jbedj1edvV4p6cQZ/hw2rRhqmgJPLD8tmqcdmFENnyQd37DIBrLkM/yowY86bdtl+GYCw9c0hject1mELtI9TCuPq2L4co19+bhWNHwAOGvjcrxu+6bQMbKkM1WhfrcKg69rsfIZZ25YFjtOBc/fEo7+mI1InL9+5vG44d0vyrxr0CLjOmD4KQbfcXF4uo7RoiGe0093erXzL33mOnGsLOkAXtb4q8/ZIH4G4FdYDbR5l8dbEkZLOtTbcNoCwO6rL8NI0XMS2y6HTruigiF22hQ11QxJxvydLz5JLFDt+BazojeegQEEGXxZkmk44Roma0ebb2flRJC0lXh8po6T1oyILaEcaXXnE0fxJ9/yWrXJBv+Ck1fjgkiGY7TBNOGxw7MiGiINIyUDE3ON1Hj9duAlXklROnJYpstDiUNFPxIl2vw7WlyrlbBMWY+VF8hNK+MZx6KvgG/kqPrlv7x5e+xYVcy43EVL/k6ZwjKbsMmSGTTpCOroKGQ+jcXKH1AjkGYo6J78ptLwaWE2dC3WdDwN0SgdVYtBGUQQiAjRd7zt8QlsWTUU2qluXjUUYutjZQN/9apn4j0vOU18frlgiPknO4xD1+iz8Jow+N7/rRp8+uyq5Rl8ObdCGPxIqYSjcw1RJ0qGiiT84B3PxzmbV4jSJBRG2gssidIKDTvYcoU0fOnnmZoX+tgsUkRmoWkhjuOz3sAOysgGk+22xybEzypJRwY1oZ6VSgG4LseeI7Op+j2BJnSz87SKoqHh8SNzoixFyfRa9qk0/JJJ3Z7CBj2YvO1IOo5/HXooaerktfF+vpTrQOev2Q4MjeGi09fFjg3CGGWnbZzhe5JOMw2fN5UPPA3fZ/i+wVfV0WGMxYzHOZuXx45TwdCZz4LjUTppfQDSEI3SobmUlKFLUTqHZ2qewfelo8mKhRPXjKBo6CIrXcXwDV0LJZkNST0U6P6p/AdDBSPQ8KWyza2CooLkgIThgiHG1NFIz4k3fvl2/Pl37419jmqMR8tbD7ykM8i4Zdc4Tv3gj4XTVE7xDxXCslz840934xkful6p7RJCKfcJSUx128FkxcLa0aJgQvJ2WmYiQ00kmWGFpPPUpKcbnpQSoRNFtyWdt194IgDgG7fvBQAUdF3otI7Lw5m2JjF8775SATZi+O05bYMondGigT+++BT895+9UHnsaMmAxoL6M9WGm6g1G37MuDw2qoqm3tk0fDeW/BWFHIcflEYODNt5UoctOUb+fS87XRSNa4bRkmeYSMOXG7anlZ9IQ3RcBww/wWnr36/xmTrWjpVCixrVPlo9WsRwQceq4QKWDxWEYVVJXHKFVVU/W3GcqccSr1qJ0iEMFXQRik1zSU5SO1ZphEjd0bkGrn/gYKisCaB2yEZ3uv9402687at39sTwL3qD//27nwr9HiqkJN38mu3gm3d4xiutJWA45V494SnRac1oUcpIDN6XSx1Ea6ZEoTL4e47MAUiP0IlCFe/fCU47bjTEME2DCZ3Wy7SVnbbeZKf7Ss5GwfDbCssMGD5jDO++9FScui7O7gHvHq8cLoptd9VyEtPrRcy4ZPCjTeLpvM1q6TTsbJJOLSrpSE7br731PNzxFxcDCEfQ/P7zTkj9XBmrhouYmK2LJvSyoYo6o7NCj0Tp0LxKdNr6ks7hmTrWjBRFkhEAbPWlybWjRZywahiMeU5RisVXNfeR22KSPyLq4wAgJJ1f7DoiFqV2JR2q1bTJz5OQn4flcJExTL/bLscPIvZHtSuke0nz4d59k5irx0lGN7DoNfzjpDrvW1YNhVZYeVLXpFKsSbWqgbDOnNTejNLDkySdo76D50u/H9eQoyBWJ0fp0OcfNxavYZ+EsS5r+EBQZgHwJjRt2+VywoDH8DkPtr1k0IjZaJrX5LuV+GN6js16rxJWjxSk0FxHJOWoQN+DoDT4piayXZMgl1ZIguy0pbox0Th9VUXLrOGTgNfUfvfhWdBwHy4YYjfRvqTj/U+STjMNv2Bofq0cB6tGCigaOrauHsbjR+awyS+j8IHLtoXmpGfAG0qGTxFLcuTPmtG4w7lc0PHThw/jR/cfEAltrUbp0PmoLMtmYfDD5vPBp6dE/gEtrt/esQ9vu3Cr8OXVLK8Qn7zjp5/krN63Xbi15WvMgkXP8OXqkFtWD0cYfjhih95Lc8bJrfiSGD4Z9BXDhVg0A+A5eJ65fgyXbItryFHoft0OmeETUx5TJOgkoduSDhDEDQOeUSdJx8u0jZciGJ+th5pyRNPWW8m0lRl+Fqwe8VguoI7ZlkHfA/AkC68qadjorBwq4GCTZJvMGr4UpVM0tESjKevjreRUrBop4OhcI2D40udkSdxSISrpBAtwssEnPxot+F9603Zc8ox12O7LVmduXIZnSyWe6TpVZCWoie8EBGskToDKUs0dKulQzEgSZMj3LMrwaQfyO1+8XRA7Ij27D8+GSpvUbSfmlCcnMO3Cl5VNXHR6exnxzbDoDb6sx4+WzFhtbPGz5Uje/BSDn4Hh01a/ZOjSxAjen5hrxKocpsGriR+OKALQUuRNL1K1yeCPFL1KjgVDQ92JZ9qSfHJouo4VQwVhrEK+jEhtlGagZ5SV4a+SGb7d3OCTQ578DlGGf8aGZTg0XU/NsMyi4XtJad65pmvxnYSMrHVuolg5XETVcoTkIBuvZteXBHq+ZOBqTTR8U2fCMNN3PGnNCL70pu2JciONWRWxCbpe2UInjzYBAoKFQUY7DJ/G6srhgrheYvgXnhJUBKAF0Ha5CJQIhYJbbugZf+a3n4XNq4Idw3f+8Hn45VUXdZQkmYZFb/CnaxYKhobvv+P5KEfqloQ0fMsVOnta9EUWp62c0Uc75hDDn6tnTpoBvEgdmeFPVS2MloyW2Fm7Wm0ayOCvHC6E6qU4HDENH/CyKGV9WnaelUxNlLG94aFD2HLVj0T3JxVod5aV4a+SNXxF3RUZpiTpJEXOUKboffvDjcVlZArL9MM7Xb+9X5pU04qMI4NKDZPBlaWIbjH8asOBxpKNaUF6TkntG6Og76syfvReteHgyEwDwwVdqeFHSyZ719I+w98ktZGkejrnnrBCtBMlx7jjchGQISsGdTtc5I3yCwjnbV3ZdX+bjEVv8KeqFs7csAznbl4RCoEDou3q1M7cKCoZnLZyFihNjJCGP9s6ww9JOjUrsyb/9687Gy84eXVLMdZZQYsWJesUSdJxOWTiGDD8WkyfJsgM/ws/fwwAsOvwbOK55Vo6WbBqpIDZuo2av5NLM/jkfAaSDf629WPQGHD/U8kGP5vTlpJtXL/fag8Mvv98aDcix/M3u74kCIMvMfyy3xBdBXkhUGUSq/CiUz3mrPJVlSRH//hsPdT5S0Za+YxWQOfbJjV4oedx/LKS2CnRAmi5rqi3E/YVttZTt9tY9E7bqaolnIvRQlV1yxVREllr7FQtR4S5JUk6DYXBl3cPc41wo4hmiNbEn1Y4EZNwzuYV+Prbzs98rlZAuwaqXU7aN0e0PLL389G5BlYMFcQuSR745YKOpyeruD+FMcsQEkJmDd8zehNzDVQtB8enGXy/uQuQbPCHCgbGymbqLsRyeFPjIsdeU7nnJLQry63yx/+haY/hl7vA8KNROqrQVRnyfcg6dt/6gq24dNs6JVkRDN9ycGQmVZ4+2QAAIABJREFUecesOlc7kg6FZMqRaWTQj1tWEu/bLofrcnAeXGOU4fci+iYrFj3Dn67a4qGXDM+4E9tuOIGeVos4cJNQbTiCXdsJiVdy3W1RSyfi3FJpi0kYKRqhkK/pqt2Sw7ZXuPCUNXjLBVvxkVc+E0DAjF0XoQYo8gCX2Z1sBIYKOh49NItXfu4XyCJf1m0XjGXXoOmZTVctZe10GWYGhg9QHZzksdJw4jV4ohD1U2zHZ/jJU7JdZkgMn6JZhqXv3m4cfqyWjuWmG3zpPFkZPmMscWcq524cma2HIsZkRMtHFwytLX2c6tifvTEw+CuHC9CYR3joOTsuhxVxjssGv98Mf9Eb/CmpoURR2j4Dng5Mjs9oqdUkEMMHktsFyho+Y14iDy0yQXRJ9lu/fKiAnQem8b279ovvlJUl9RIFQ8OHXrlNaPmBhh9x2srsLmGyy1vvpJ2TjJrloGQkSwiqawW8ydfMOBV1TWTaphn8kqknjhXH9WrXNPMxBAy/uaTTriOPFrfZWjxevd2wzGhpBU8mSwl1lZqhdEOjLpveZ0xWGth/rIrjl6tDlKM76Vb62cr4xG+ehTc/fwtOkZIdX3X2evzgHRdg7WhJjHfbDZrF030OG/z+Mvz+08QewnV5KPJBrkZXMnU0HNcrYayxUOXLNIZft1yxtU4qrRCtyqcxFkg6bdTkfu/LTsP37t6POx4/itc8e6O3iPUgrr5TkKTjNd2QSysEAzypGqPMuCkyJq3het12Wwqvk2vk1JqEZZoGQ82v7vmBHzwAoHWGn7W7EklSXh5IPGQvCkNjotJkVtA1VKx4hdb2GX480zZt10TXsLzcWZluArHnG3ceRtVy8MJT1b0zosXg2tHvAS8q64xIsbqSqYteCiGG75DB9wmmE5aLS6aOG9/zorZ9Mp1gURv82YYNzqUa2sLgE8N3Rdzz4elsBr/huL42n15aQWNB3XONBVE6gbMx+8NeN1bCppVlsXNoFr7XLxQMbxEtci1WD5/gTfj438qDn+SrtIbrdcvNrN8DweJr2a6vN6cnXk1XbTwgOWRVx6cx/Kw7OZmE1K3mi9juqy9LfV8FUTXS38XKpbLbjd6K5pfQjisJ5BxWJVG1AxovP37gAIYKOp6f0NgnSjB6JacIhi9VJVVLOg6KhtZSlnw3saglHVFutkQMP1ycqOG4KBgaioaG/ceC1mNpjS0s/28MXRNaXRQN2w0xCcZYUGSKCj21OPAo3b9uO6g0nK5NnG4iKdNWNgTLhwoi6kKWFuSfRcu7lIW3ZjstsTU6ds4vcZslDp9qxv/nu16gZKXdYPjFFiSddkGLHX0f2f/TsaTDszF8uk9ZKrxmAe1SLIdj04qhxPsWNfjtMvxmMKQoHfLtyZIOVQolht8vdMTwGWMrAXwbwBYATwD4Lc75scgxmwB8DcBxAFwA13DOP9vJebNiqhquLy47yADP+C4rmyiZeqixcirDt722dabfxDvtGILOWFC8zWqd4QN+eVnbxY4nvNtLzbcHCVQgK5ZpK7HWFUMmPvTKbTj3hBU4f2tQGEyeBHSPGgmSGeDd41bYGk10ylJOddr6CxeV0YiWJZavOdqYhiCH5qYh7LR1esJAvZaMQekKOWGv01o69IhqlovVI8n3lBzF73jxyW2dL4qyqYsmPGm73aiDuFcGnxrJ2C4XO9OyxPCv+t79+N7d+2H4ZUT6hU7PfBWAGznnpwC40f89ChvAezjnzwDwXADvZIxt6/C8mTAdyZKUHWSA56gtmzqKphZy2qYxS2Lvhl9jXHmM44YMusYCPbodpy0QRMDcsPMQCoaGF5yyuvkfzTO8OHwnlmkrM/xt68cwVDDwW9s3hViz6n6kPYcsSU0yaAEmEtAshLBuBww/KRwyrYBa1mYbdG/qVvMonU4gX4fsNG1WvC8J8SiddGfk2194Iq7/0wsz1/BvBsaYmNdpu91ovkvWpK9WIWv4RAQp8eo9/3Evvne3F3Bhu7xlstdNdDq6LgfwVf/nrwJ4dfQAzvkBzvnd/s8z8Nogboge1wtEa87IDjLAi1oYLRliO0gFwFKdtr7BN3UGK0FjrkcYfshp22aJVmogcc++SZy7ebkyq7DfoEXJiVTLlBl+UjkIlXFP6zdQj8hmWa4NCAx+mqRT9CUdYvhJzjXq1au+vmyZwKEoHcvtmTGQr6MbdZViDVAajjKrlVAydZx+3FjH55VBUm0aw48+u5PX9UY7l6N0bFGVVH2fs5YD6QU6PfM6akbu/59a8YcxtgXAOQBu7/C8mRANqStFGm3M+AafHsDasWLTxhbktKUm3spjIkyNyU5bEaXT2sQu6Bosm6NmucrG54OAsqmjZnl6pVwPnxbaC1N2JapFthnDbyWBJmrwm5ZWcNzUsrtA0KtXhay110NO2x5JOtHr6IbBj0bp1PqQUETMPs3gR30vp7TQQ6IVhBh+xGkbRS/8NFnR9Mkzxm6Ap79H8YFWTsQYGwHwPQB/yjmfTjnuCgBXAMDmzZtbOUUM01VvwgaJV8HkshwvWmOkaIrX142VQmWSVSB93tCbaPjSBNM1JkXp+MyvxVXeNDRUqhYaPTQKnYL00mNzVkgb1jSGW977YqwdS64fpFpk0xi+5aQ7XqOgxSHQ8NOjdCzbRaXhRfMkZaMWzWRykDks0/8Os3WvgXmvDb6hsZaS/pIQLa3QrAJpL0C+lWYRa7dedRGu+LcdeOCp6VArxW4iYPiSpJNg8Aea4XPOL+Gcn6H490MAhxhjxwOA//9h1Wcwxkx4xv4bnPPvNznfNZzz7Zzz7WvWqGNrs2KqakFjQZVBeftMtWlGS4YwvseNlZo2tmgISUdLl3SMqKQTCctsOUqHwbLdTOn6/QIZ/Nm6HdOGN61MjqQA1Oy/2cLbSqVH4bStZdPwG443RtKShNIYfnanbXjn0Sv2Jyc+NeuylgX0fDn3Sgn0I/qEwbuGZU2yztcvL+NUv/Xliat7xfC9+0sJd0CybNhPht+p5bgWwJv8n98E4IfRA5i3p/oygJ2c8091eL6WQFm2NDjl7TP1ohwpGYL1nbx2xHM8pkWHUFimxlIlHVluYEoNv0VJx9eVWzV08wnZIaa3mFxz0enrsPOjL8MGqW1fWsYthcdmRRCl4z331MQrXYPlcMzW7VRfSZqG32riFTUk6RX7I9+ArsV747YDer6OK5eq7o8hS8relvHRV5+Bf3/7+aIUcbchx+FTuHZS85uBZvhN8AkAlzLGdgG41P8djLH1jLHr/GMuAPBGABcxxu7x/7WePdIGoglKKoM/VjJw995JAMCrnrXei9BIiMOn1bug614cfkp5ZFmyCUXpWG1KOr6u3GjR0M0n5BC4dqI/ygU9pC83i5ZqJUqHmkxMZQjLpPs7WbFSsyFLhg5LSrSRkXVh1zSvrHSvGX5RYvjdkI0ofN/hXGpvOL/jkvu9orJknY8UDdGNqheQ4/BpPJgJOQ6tJAx2Gx2JeZzzCQAXK15/GsBl/s+/ANAXShotQUAr66GZOv7gX+8A4EWN/P3rzsYvdx/BSWtGUptTk6ZMUTpJ5ZEbthvqwxnW8DuI0rFdf/fQvwGTBrlQVbvx3XIIZJqG3+rCp2kMps4Cp23KpCsKg99IrVBJi3bddmI7gYaTvX9q0ZQMfq8YvqThd6O0gS5F6TRrb9hrGC0s/L2CHKWjOeHXAODtF27FF295HEDvnnEWDF5sXxcRLTJGk/zbd+4T/VVHigZe++yNoj5J0dDx4NPTmKvbsckux1YbzRKvYhp+5DNaHKSerswHm+FL97rdsrvyPW+WeNXyPZSYdBrDp/cOz9RxckpUB2VL1ywX0RJBrdRMKpm6xPB7a/D1NjNro5CjdGp9MvhUErnVrPVeIJC4JIYvSa9y1c+8PHKPEO0gpGlMlFIgREPUjlUaODJbxzU374l9nqzLGrqWXDzNcUMdfkJhmX6UTassyytb4PiGbjA1/DGpTo7WJouU67w0DctsI3ktCxul0g8HpmrpTltRfTUuAcoVU5uhZGrYecALXOudpEMafnc+jyQ7lwfFwuabiHz08jPwkVduw3lSxna/oIrSkUmPXOJhIWfaDjRU2ZglQxMOMiBcORAA3vICr1v8xFwdUYhJrDNf0slWWkFjTNLw28umNHUmjNWgMnxdY0JCa9ewyGUMmjH8Vrs1yfctzWl7/LLAcdzMaQtAWU8nq9MW8Jratyv1ZUUg6XSJ4UuMlohPL9popmFZ2cSbL9jas/6vrSBUS0cwfKmkyLC609t8YzAtR5dguzxW/rUkdbEH4g6fN5y3GcvKpjLKJMTwUxKv6pHCXp6GH7zXTjalKTmJB9XgA4Hjtt2U/ZCGn8rwWw9PlUMT0yKd1ku11ZPq6ADxzG0ZrRjwj15+hvi5V/qu/N27AZnRksFPikpZCgiXVvDuRxLD76fTdlE/Ia9qY4ThS8Z2/bKSckImlU0I6twHhZtUiNZEIUnHcTmenqy1xeJk49ZOi7b5AtUuaZftjUY0/Ft2jeONX749FAnDue/LaPE+EONK670KhBN50iJAAqdtfGFqxVcj17efjygdIH2Hk/XzGPNkU5oHSVEpSwG6FIdPDF+eA3IEW+607RFsN95ijiJ1yqaOW98fCzACgET2HmP4KVE6slHX/GqZf/H9+/HzR8dDseZZIW8PB5lJrRv12HGrcfiEaJTOPXsnccuuI5it2SLeut2dDhnfZnHQ8mLwvIQ660Azhu+01E7v5LUj2H14Vhni2Q3QrpIW5Lv+8hJ0cirGGIb9xvM0VwY1P2Q+YEg7HgoXlaOH5HGdM/wewXbikg4xm6h2LyOpbIIcamckaPiCfYYMvsdWv71jHwDgqclqy99FZoqDzPDX+eUTuiHpNGxX3OPZRlCGOPCltHYfaBFuRUN9zpZkh2Aaw2/Ybkvt9D7xm2fipDXDOKlHtV7oXh2/LOhF0GmrwXJBR9WyxfMYhPDIfiGIWnKF7ZDJpkwAc4bfI6ictsR00gZ7UtkEeZtu6ppyUXD8jvVmxGl7//6p2LGtICTpDDDDp/621ZQmMmkIJV45XOyiKlLdedL2W2WUdN+yyBlfe8t5ODrXSL3XxSYafisTe/uWlbjxPb+W+fhWQd+Dnk83MFTQUZEknUEmIr2GzPBpWMpkU743/YzSWZQGf7LSEFpaXNJpbvCTyibIUkLSMSoPvcaYiPtvF+aADJhmIINyZKa97ytHxTRsJ2D4dQXDb6M8BZAeg09I6pEqIwjLVDP8QTKA1O+BYte7gbKphySddvvjLgbIeQkuC78GeBJYwdDA0H4z+m5g0Rl8y3HxnI/dgLddeKIv6YQn3eoRT8NMM/i6pnbIyiVvjYRdQEOhZ2pa8PqVLzoJZ29cFvu7ZpA/r9VwxPkEGRS5KXwrkBczywlimqkZCRA8h5YZvt66pJMGulYVw2/Yvatt3w6O+mHG0YYgnWCooKPacKQxP7jjsteg4BDb4dBY4MS+9aqLxGJb9A1+P7HoDL6pew2CHz4wDcuNFxo70S+PmhZQYOoaHIVDVnbamnoCw3fUDJ/we8/djI0rWi/gtNAknSMz7Rl8+bs1bFc4MedUGn6bYZndMviC4Sc5bQfIANIOM9rjtRMMFQxUGrY05vttzvoHXVcwfJ1h/XAQoFE0tLYTEruFwRmRXcTpx43ioQPT4DyeaLLVL486MZssOSQ5ZHcfngXgMUUvkid+jKXY3soPuV1H2UJx2p6ydhS6xvDONnuXygul5bjifs7JGn6bTtsgLLM796+Z03aQFuY/veRUnH7cKC44OTnqqFWUhYafM3xZwxeybsT2FA29rw5bYLEa/OPHcGjaY5hRXZEaIBxOYaCmpmHf0Qquf+CgeO3hg9P49A2PAvCY4ljZwFzDFvXVCWLwazLDD95PK8aVBnkyDZIhiaJc0PHY1Zfh1886vq2/3yKVr23YriTp2KHXgdbvA2mqXWP4zZy2A/ScztiwDNf/6QsTW0y2g6GCl8RoKfxWSw3hKJ144hXgMfx+hmQCi9Tgn3bcqPg56rTdstozKKTlq2DoDE9MVHDl1+8Sr1E5ZcAzNC84eTVcDtz86Hjob0USihFn+EW/cUo7kGPvF/PEWjVSxJ6rL8Nrzt2IuhSWOSc1mW+XUVKBspdsUzVwax2mzsDYwmD4vcBQQcdc3Wk7amoxgfJO7ITEK8CzG/0sqwAsQg0fAFZJjqmo03aoYOALv/dsnJXiOJX/xnG578QNJrWpazhn8wqsGDJx08PjeMVZ6wEAO544ir/+r4e8zwgxfO/BdxL3PChhXfMBTWN+yJ8twjJlhl9vk+F/4LJn4DXnbmx79xEFYwwlQ09k+N3oHTvIKJsGqtIzWsxEpBk0jUFjQVi2xuK5KJ0Qvm5hUY7IQogNx1nHy85IZ3im9KAsx4Wu6cLInL1pOcZKBhhjOHHNCA5N18SxN+86gnv9ePuQ09b/sV05Bwin+y925ggEMd6qsExL4RjPglPWjeKUdaPND2wBJVNLLJ622J/TUEFHxQri8JdyWCbgZ+j7Bl+VhHbGhmUL2+AzxlYC+DaALQCeAPBbnPNjCcfqAHYAeIpz/opOztsM8k1tpzqgPHAtx+vVSbrx1b9xhoijHSroIUMkMz1T4bTthOHLrdkWuyEBPF9A3W/4AgAVRVjmIOx0ioaeWB65n71L5wPlgg7OA7lzKdfSATzN3nG9Hr+qWlIf+40z+3BVYXT6hK4CcCPn/BQAN/q/J+FPAOzs8HyZIMsf7bAOeXUmp6Gq+iHFIRPkn1VhmZ0YfPlv+80S5gOUDTvrGxO5tMIgRYUkMfy65Sz6hZl6TUxVLegaa7ucxmIBNUWynHgNr0FBpyPycgBf9X/+KoBXqw5ijG0E8OsAvtTh+TKh2ETSaYaQpOPrk6IXrcTahgtGKD5cZvjyQkNRmWmldlvBYjckQGBMZuqeo7UbUTq9QBrDH4Tr6yXoGU3XrCXtsCXoOoPjuqhZblcaxfcCnY7IdZzzAwDg/7824bjPAHgvgOQC5z4YY1cwxnYwxnaMj483O1yJziWdOMOnZJ8Qwy/qIalBrh8j7zIoPKsTDV/GIMfhdwtlv8TCdNUz9PJ9nvGNf6fFv7qBZIY/WGGZvUDwjKwlL+cAPsN3OWq20/donCQ0fUqMsRsYYw8o/l2e5QSMsVcAOMw5v6vpwQA459dwzrdzzrevWdO8nokKzZy2zSD/jZB0rDirjDP8YOIbCkmn06iN8/1WbkvB4AuG7+c5yL6SST9rVHZk9wtJDL++BBj+sP+MJivWQJfsni+Qhl9tOB33G+gVmlogzvklSe8xxg4xxo7nnB9gjB0P4LDisAsAvIoxdhmAEoAxxtjXOee/1/ZVN0GnDF9XSTpCww8eZLmgo2a5InQzyWl7YMqL5Fk7GnRSagdf+YPn4Klj1SWhlZaFwfcMvbywTlYtDBf0gTCoRVML5WgAfolse/E7bcf8BffoXCOXdBBE6dQGrI6SjE5nzLUA3uT//CYAP4wewDl/P+d8I+d8C4DXA/hpL409EDa2ejtOW00h6Sh042F/S0tSTtjgB8dRg+rfes6mlq9FxlDB6HpY4aCCGJJIvJIkncmKheVdrAnTCUpmPA5fJf8tRlA3sCOz9a71yl3I0DQvDr9mOSgN6LPvVAT9BIDvMMbeCmAvgNcBAGNsPYAvcc4v6/Dz24JcfrQdbdGMhGUCXjEsQ2Mh9j/kO2ErdRsjRSOk4cte+q+8+Tk4MFVrq9PVUsVQxOk1V7dxYKqK45eVMVVtDIScA3hGvRHJtB2ksNFegp5BfQnkHGSBYPiW09Uidd1ERwafcz4BINYnkHP+NICYseec/wzAzzo5Z6voOCzTDcIyoxOYGD6l/Scx/BefnuTLzpGEqMGvWg6e9/Gf4htvO99n+INh8FUMv91M4IUGedEd1DDE+YSn4buoWYOr4S/uEYnOwzKpEJKqvjkZJQoZlJ22gxAjvpBRLqi5yA07D2GyOjgGv2hosVo6QbP7xT0GSqYm5lc+3r1Fz5N03KZ9k/uFwbyqLqLTsExLJF7F65tTd6aKkuHnjKcTDCUwpCeOzGGy0sCy8mBsmdMYfr9L4fYajDHB8vPxHkTp1KwFHJa50NGepCMxfClKJzqBhYbvR5CENPyc8XQEOXFFlgtuemQcR2YbA8Xwa4kMfzAnfTdBjtuc4Qdx+NXc4PcP7QxEMyFKJ8rwhyWGzzlPTLzK0TqKhiYylGXjTtqoqkJlP1Ayda9/slRNleLyF7vTFghCM5d64TQgYPh1y80Nfr/QjjMpWjwNSGD4kobfcFxwrv6MHK2DMSZknTHJOfijP34BgGxNxucDpNXKLH+QSj/0GmPlnOETDM3z5zScwdXw+5+b3mO0o+HLLQmDKB0nlkhDBr/ScFBrhLf1edRC5ygXDMw1nFA0yIlrRvD4xy8Lhd72EzQm6pYjSj0slbBMAJKGv/i/azPoGhMd8PIonT6hHabtSP1sLSlKJybpFANJpxZJrx8Ug7SQQQtqNOZ+kO6tiuEvlbBMACK3ZFAktn7C0JmI2MslnT6hHYPvStqMXB45Kul4Xeg9p61cGjlHd0AsafmAJFmpIDN8wlIy+Jdu83JMbn1sos9X0n/oGsOsnxE+qJLOYF5VF9FOpq3M8EWUjqL6IWPMK6BWjzP8HJ2jHGH4g6iSCYZvqZy2g8nyuolzNq0AADz7hBV9vpL+w9AYZv1y3oPK8Be/ht8Www9+tqTyyAXFBC77vVdzht99kKQzVjaxYXkZ77701D5fURxk1OUFfylp+JrGcPtfXBzLjF6K8Aooes8+N/h9QjvOpLCkEzRAUU3g4aLnWFTVRM/RGciImLqGX151UZ+vRg2S+epWXMNfCgYfANaNdVYFdrFADhAZVIO/6EdkO9EybkjS8as1JtS49toc2rnTqgeg8gr6IGo5Pmhiqxj+UtDwcwSQx+mgVssczKvqItoxFlFJp2Y5mKpaWDtajB1LGn41N/hdR9lnz4Octk8sXmb45PfJs62XFmRymTP8PqGdEL7Vo0GdFttxMT5TB6Deug4VPQ0/Z/jdx5Bg+IM7TGliy12vyO+T52IsLcjkclDzEgbzqvqMNzxnM/7xDecAACyX49C017Fq3TKFwS/omGvkDL8XKAsNf3ANp5Lh+wZ/UCd9jt5ADhAZ1DHb0YhkjK1kjP0PY2yX/78yNosxtpwx9l3G2MOMsZ2Msed1ct4s+MbbzscfvvDEtv5W0xheefZ6rxiS4+IgGfyxuKQzVDBQqdu507YHoNIKC03Dt10XjA32defoPuTnPahyXqdXdRWAGznnpwC40f9dhc8CuJ5zfjqAswHs7PC8TXHByavx/sue0dFnGLpX/e7QtCfpHKeQdIYLOiqWk0s6PYBg+AMs6RDDl5+/5fCBvuYcvYEcpTOocl6nYZmXA/g1/+evwutm9T75AMbYGIAXAngzAHDOGwAaHZ53XmBqGizHxeHpGgqGpmyrN1Q0UKl7Bl9jwC3vuwhcrqKWo22QwR9kpiw0/JCk4+bF85YgFoKG36nBX8c5PwAAnPMDjDFVL78TAYwD+Apj7GwAdwH4E875XIfn7jkMncF2OCZmazhurKR0AA+ZOhqOi5majZKp531ruwiKwx9k42nqmpdwE5J0+MAyvBy9gxGSdAbz+TddhhhjNzDGHlD8uzzjOQwA5wL4Z875OQDmkCz9gDF2BWNsB2Nsx/j4eMZT9AaGrsF2XTw5MYfNK4eUxwz5BdQm5hoDWyFvoaJseve2nYqn84mioYUYvuW4A8vwcvQOmszwB3TMNmX4nPNLkt5jjB1ijB3vs/vjARxWHLYfwH7O+e3+799FisHnnF8D4BoA2L59e1+1EVNjsByOx4/M4VXPWq88ZtQ3+OMztYGNvV2oWAgMH/DbHMoM3+EDf805uo9FwfCb4FoAb/J/fhOAH0YP4JwfBLCPMXaa/9LFAB7q8LzzAkPXMD5Tx3TNxpZVw8pjKGZ/39HqwFbIW6gQBn/A5ZGioYWitCzXHfhdSY7uQ18CBv8TAC5ljO0CcKn/Oxhj6xlj10nH/RGAbzDG7gPwLABXd3jeeYGhM/z8UU9W2rpabfDXjnqRO09NVkN9WHN0jmeuX4bf3r4J524e7EqMJVMX9XMAj+EPahx2jt7BWAySTho45xPwGHv09acBXCb9fg+A7Z2cqx84NFUTP5+ydlR5zBqp3EJpCZTDnU+UCzr+9rVn9fsymsJj+OE4/EGNw87RO8gZ4dqA7krzUZmCOb/k8d+95ixsXqV22q4aLohm2znDX5oomXosDn/QZagc3cdCeOa5wc+A123fmPieoWtYNeyx/KXQ8CJHHCuHC5iYDVJL7DxKZ0likPNFCIu+Hn4n+NhvnAGdsaYF2NaMFnFktp4z/CWKNSNFPPDUlPjddvMonaWIhfDMc4Ofgt89/4RMx60dLWLngcGtgZ2jt1g9WsDEXAOuy6FpzIvDH1CnXY7eYSEw/HxUdgHPP2kVAGC6ZvX5SnL0A/9/e+cfbdd45vHPN/fmXm4uiSRCoiIhQgih4kdaJipJtdaq0KXqV/0sSutHV3XG6NJMpanMoDWKLqHUoqQoVcQY1OgwUx0iM6Sp34M0hqAihPy6z/zxvMc94t7kJmefc/a5+/msddbZv859v3fvdz/7eZ/32e+7eXsrqzuMvy7zsE7k4ReTRojhh4efASfvty0d1mn4g2IxOGVqvfneCga1t7Kyw2iLGH7hyPO8DSXC4GdAnz7itP23q7eMoE5s3u4Gf/HS5eyw5SbeadsA3l6QLY3g4ef/kRQEOafTw/dhtCOkU0wihh8EBWBw+8cN/sp48aqQlDz8DZhVtWZErQyCCikNoLf0w1VAGlqhAby9IFtKHn6er3wY/CCokD59xMZ9fTJ7KE2AErdW0SiF8db13k49iVoZBBnQr7WJZWkojpUdMXiM5hLTAAAPaklEQVRaESll6eT5yofBD4IMaGtp/sjgr1odwyMXkYjhB0FBaGtp4v3lnTH8yNIpHp0x/Pxe+zD4QZABbS3lIZ0YPK2IfJSHn197X5nBlzRQ0v2SnkvfXc5UIenbkuanuXBvlrRRJeUGQd7o19rM+yvKPPzI0ikcRcjSORd40My2Bx6ki7lqJW0FnAmMN7OxQBNwRIXlBkGuaGtp4oMVqzGzNFpmePhFo9Rv05tj+FOB69Py9cAh3RzXDGwsqRloAxZVWG4Q5Ip+Le7hr+owgMjDLyClfvreHMPfwsxeA0jfQ9Y8wMz+AlwMvAK8Biwxs3+tsNwgyBVtrU0sW76aVavd4IeHXzx6hYcv6YEUe1/zM7UnBaS4/lRgJDAM6CfpmLUcf4qkxyU9vnjx4p7+H0FQV9qSh7+ywyczjzz84lGK4ffJscVf52iZZja5u32SXpc01MxekzQUeKOLwyYDL5nZ4vSb24HPADd2U94sYBbA+PHjbd3/QhDUn7aWJj5c2cGKVW7wo9O2eDQXoNP2t8Bxafk44M4ujnkF2EdSm/yd40nAggrLDYJc0a/Ffad3P/BJcCKkUzyaentaJjATmCLpOWBKWkfSMElzAMzsMeA2YC7wVCpzVoXlBkGuaGv1+YyXJIMfIZ3i0Qgv21U0AYqZvYV77GtuXwQcVLY+DZhWSVlBkGdKHv68V98BiKEVCkgR8vCDIMBj+AA/uOtPQGeKXlAcmhSjZQZBIei/cd+Prbe39u3myKC3IsXgaUFQCDbr1/LR8pF7DeeAHT/xSkrQyyn12ZZmQMsjMYl5EGTAZm2dBn/PEZs1xPymQbYMaGvhR4fuwv47bF5vKd0SBj8IMmBAW2cIZ2CZtx8Ui6P2Hl5vCWslQjpBkAHlwyEP6pffJn1QbMLgB0HGDGwPDz/IJ2HwgyBjBkVIJ8gpYfCDIGM26ttUbwlB0CVh8IMgI0YMaqu3hCBYK5GlEwQZcfeZ+/FBmtc2CPJIGPwgyIj21mbaW+OWCvJLhHSCIAgKQhj8IAiCghAGPwiCoCBUZPAlfUXSfEkdksav5bgvSHpG0vOSzq2kzCAIgmDDqNTDfxr4MvD77g6Q1ARcAXwR2Ak4UtJOFZYbBEEQrCeVzni1ANY54P9ewPNm9mI6djYwFfhTJWUHQRAE60ctYvhbAa+WrS9M24IgCIIask4PX9IDwJZd7Pqemd3ZgzK6cv9tLeWdApySVt+T9EwPythQBgNvVvHv95S86IDQ0h2hpWtCyyept45tutuxToNvZpMrLHwhsHXZ+qeARWspbxYwq8Iye4Skx82s287mWpEXHRBauiO0dE1oya+OrqhFSOe/gO0ljZTUAhwB/LYG5QZBEARlVJqWeaikhcAE4B5J96XtwyTNATCzVcC3gPuABcAtZja/MtlBEATB+lJpls4dwB1dbF8EHFS2PgeYU0lZVaImoaMekBcdEFq6I7R0TWj5JHnR8Qlk1m3/aRAEQdCLiKEVgiAICkIY/CDIIVrH24xBfWnU61MIg5+XiyOpEOd7fcnR9cnTYPZ96y0AQNLg9F33eRsl9S9brnedach7uSFF9wRJO0vaH8Dq2FEhaRdJ30k6OuqlI2nZTdLJkrp6ka7WWsZImgD1vT5JywRJVwN71lNHmZZbgYsl7VQPQyunTdLNwJ0AZla3qbwk7S3pTuAaSSdKaq1XnZG0l6QbgQvTvd1QNrShxPYESX0kXQn8GjhP0vTSSJ51ujgzgB+VHj51uoH7SroK+DkwEZghae9a60ha+ifjOhuYLmmGpFH10JL0nIxnVcwFnqynJytpCHA5ntH2JnAWcGLaVzOP1pxlaXWwpNOShprfP5J2xQdfvA24FTgAqHl9SXZlGnANcC+e4fhNYFyttVRCrzP4wGbAJsAY4GjgLeA7ktpr6WGXhQd+D/wz8ENwT6kON84uQH8z28PMjsGve71e/f4unh02DjgVGASMqJMWgOH4MCE/M7MP6+nJ4sbjWTO7DrgEuB2YKmm0mVmtjL6kZklDgdeBk4DTJA0ws4461N3S4Is3APcDGwGvlGmtyTlJtuNl4Hgz+yXuyG0D1D3UtT70CoMv6dOSRqfV/sBngDYzW4x7+m/jT+OqVpCkY3vwF87SzXEgcDXwhqSvp30d1a6oa5yT1cDhybv+MrAPMEnS7unYamsZKWnjtHo18H0AM3sBGIA/kGpC0tKalgcCY4E/SjpA0n2SzkvnqBbn5UhJP5B0cNr0JDBe0nZm9j7+lvrj+IOxaqGvMh1fSuWsMrPXgJHA/wIPA+cmXVV1msq0TE2b7gIOlTQDeAofmuUySX+XtFYttCNp4hot4dnAvBRSegtYCgytVvlVwcwa9oNXyHuA/wQeA6ak7TfgXht402sSfrGG1UjHAWX7LsK9kk8Dz+DN0k/V4ZzMBH4JvAF8DZiO30yjq6hlBN78fRB/8O5Qtq8lfV8HHFyDurKmljFp+zVp22X4sN0nAPOAcVXUIuAbuIE/IdWLr6d6cj5wWTquD7Av8DNgaI10nAD0w73XS9NxBwPv4mGvVqBvjbScUlan/wk4Nq1PTHV3QpWuzyZ46+pt4FpgYElj2TF9gf+o5v1TjU/DefhreF3nAPPMbALeuXRi2n4t8FlJI82Hdngd+BDYmIxYi47f4DcvktpwD2AkHl7aAhhiZguzjBX3RAvw9/jQFoeZN48vBV4CPpuVjm60PGZmk4CH8Jj9zmlfKXTy0fDZWYcL1qLld8APJY0EpuEtjEVmdqd5OGUObvyrgrnFmADMTOV9E9gfd0zmAKMkTTb3pt/Cz9GSGumYDOwH/BUYKeku3Gl5GHjZzJab2coaaZko6Ytm9hIet1+YDn8Cd1yWZ60jsQKvI8fgAz0eVqaxxBjgdTN7VtImkvaqkpZMaTiDj3tBpZv5faBU+TYFFqQOwEfx5vDFAGb2NO6xZFlButPRP+kYY97xtQr4I9COdzgNl7SrZRsrXpuWpyXtlCrrcuCrAOZN0q3IfiKakpZSH8b8VN7leDz2KElDzPsyRgFvm9mTqWPwfEkDaqDlCmAPfBjuxbiXf1jZ74bg3ltmSDo2hQgGpk0LgK0kNZvZA/jscfskPTcBl6bzMwn3fltqpON/8FbFaOAvwIvAHmb2JWBrSXtkoWM9tOyfOrPvA6alOn4EsDP+MMxaywAzW47XiQeAZ/Ew2+h0XKkuDQSWSToeryu71Ko/oRIaxuBLmiLpfuAiSYcnA/YIPhLnk8AX8A6Um/Am34XAUEmXS3oa73BZUulF6aGOZuA6SQfh4YPdzexUM5uLx6/fqUTDemppAq6X9HngX4ADJV0s6d/xB8OLVdKyCm8S7y5pnKRxuFHbBu+oBdgW2FPSQ3jYYLaZVXxueqhlPt5hO9zMzgNekTRT0h/wm7niAf7kDE3/33F4K++nkjbFWzVD6Mw4mY33Jwwysxvx8Nu5uHH720rOy3rquAWfinQlcLaZnWVm76V9k8zsiQ3VsQFafoU/eIaZ2VX4A+he/JycaGYvV0HLFZIGm3fgr8BDo28Ah8NHg0GC980dCfwNcLSZ/XyNFkA+qXdMqScfvAI8hjezd8eN+jlp3w7A7WXHng9cnpa3wDtwM4kRr6eOacAlZet9gD51OiffB36SlnfDOwEPraKWm4HT8Vjo+cDd+INofNJ5Zvrd0bghnlxHLd9Ov9sU2BH4fEY6mtL3aODGtNwMXAlcj8eAr8X7U/qn/b8AZpT9jZY66bgeuCAtK6t6W4GW6Wm5L7BllbX8FPj1GscemjSOwpNBwO3KV7Oqt7X65OnNwo9RiueaxzH3Bp6wNMOWfBauH0u6ATcYr6YQygI8Vny2JJnZ63j8vh46HizTYZZBdkMFWn6XtPQxs3l4p2S1tVwC3Gpm0yVta51zGj9KZ2httnmKWx60LDWzPwN/rlBLM3AB0CQfInxTUn+FeebWt4DXcC/6JuAQPPPkQqCDslCSuYdZDx2r8Ycm5tatIs81Ay1/SMeuBP6vylrOBBZJmmhmD6ftd0gag7eQ2yV9zswyDfnVilyGdCSdgHfQTE+bngKOlDQirffFQxHT8dSogcCZks4CrsJjb71GR4ZaMsnl7oGWZuAF4Cdp/aX0u1PwvO65kM3bmxlqqbg5Lmki3qG4GfB80rQS+JxSp156KF0A/KN5nHoWsK+kx9Lv/q236GhQLZa0/EPZ774CfA93JndNTlRjUu8mxpofvHPzN/hbhnOBHdP2S/Hm+aPAjXhmxb14CtkY4Ay8+bdPb9LR4FruAbZI+8/GO9L37I1a0t/dD/ha2fqVwGnA8XirA9zJ2hJPzx2Rtg0AtuptOhpcyy3AyLLf7Zellnp96i6gmwszPH3PBH6Vlptwr3XftL51MmYVxznzrqOBtfwCaE3rbQXQ0obnqZfiw0cDF6blecAZaXk8cHMVr08udISW/H1yGdIxs9Kr05fiucAHmjf/l5jZI2nfN/AUxKq9Cp8XHQ2spZSWinWOzdKbtSwzz1Mvnf8peJol+MtEYyTdjbc+5lZDQ550hJYcUu8nTg+eyqcCD5et74W/ZDWHjHrsG0lHaMm/FryF0QcPr41K20bhYYp9yThUkXcdoSU/n1xPcZiySjok3Yb34i/HOx+fMx+HpVA6QkvDaCm9KHUNPufzifhLQmeY2btF0xFackS9nzg9eBq34SNOvknK3y6yjtDSMFr2wdMsHwFOKrqO0JKPT27z8Ms4HY+nTTF/5bnoOkJLY2hZiKfy/bjOWvKiI7TkgFyHdKCzqR46OgktXZMnLUGQR3Jv8IMgCIJsyGVaZhAEQZA9YfCDIAgKQhj8IAiCghAGPwiCoCCEwQ+ChKQBkk5Py8PSS1xB0GuILJ0gSKQhle82s7F1lhIEVaERXrwKgloxE9hO0jzgOWCMmY2Vz1t6CD4Gy1h8QpUWfGam5cBBZva2pO2AK4DN8cHaTjafVCUIckGEdIKgk3OBF8xsN+C7a+wbCxyFD8g2A1hmZrvjc54em46ZhY/HsgdwDj7eehDkhvDwg6BnPGRmS4GlkpYAd6XtTwG7SmrH5zm9tWxSsdbaywyC7gmDHwQ9o3y8lY6y9Q78PuoDvJNaB0GQSyKkEwSdLAU22ZAfmg+r+1Ka/xQ547IUFwSVEgY/CBJm9hbwqKSngYs24E8cDZwk6b+B+cDULPUFQaVEWmYQBEFBCA8/CIKgIITBD4IgKAhh8IMgCApCGPwgCIKCEAY/CIKgIITBD4IgKAhh8IMgCApCGPwgCIKC8P+ZekknfSFEkgAAAABJRU5ErkJggg==\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"ssta.plot()"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"7. Convert to days to deal with differences in lengths of months."
]
},
{
"cell_type": "code",
"execution_count": 17,
"metadata": {},
"outputs": [],
"source": [
"x = (ssta['time'] - np.datetime64('1990')) / np.timedelta64(1,'D')"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"8. Compute linear trend"
]
},
{
"cell_type": "code",
"execution_count": 18,
"metadata": {},
"outputs": [],
"source": [
"y = ssta"
]
},
{
"cell_type": "code",
"execution_count": 19,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"0.00535 K/yr\n"
]
}
],
"source": [
"slope = np.polyfit(x, y, 1)[0] # in degC/day\n",
"print(f'{np.round(slope * 365, 5)} K/yr')"
]
}
],
"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.7"
}
},
"nbformat": 4,
"nbformat_minor": 4
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment