Skip to content

Instantly share code, notes, and snippets.

@bradyrx
Created February 26, 2020 22:20
Show Gist options
  • Save bradyrx/937d96c6ba56fd64923a5ac8767ef2c7 to your computer and use it in GitHub Desktop.
Save bradyrx/937d96c6ba56fd64923a5ac8767ef2c7 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": [
"Use datetime for time axis."
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {},
"outputs": [],
"source": [
"ds['time'] = np.arange('1854-01', '2018-01', dtype='datetime64[M]')"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"1. Slice to North Atlantic"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {},
"outputs": [],
"source": [
"ds = ds.sel(lat=slice(0, 20), lon=slice(310, 340))"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.collections.QuadMesh at 0x2b339cb81cc0>"
]
},
"execution_count": 5,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAZAAAAEXCAYAAACDChKsAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3deZhcVZ3/8fcnTUIISUAI+yLLBBSVNbKIDKKogAouoKAyuPxEFBQHcUFQGNQRlxnFhSWKgBpQGBYBkUWGRVD2HQKKLEMIskMWEpLu/vz+OKdI0anlVt3q7kr39/U89+m69557zqlKp07fe5avbBNCCCG0asxwVyCEEMKyKRqQEEIIbYkGJIQQQluiAQkhhNCWaEBCCCG0JRqQEEIIbYkGJCxF0vqS5knqGe66hBC6VzQgAUkPS9q1sm/7/2xPtN03nPVqRNIhkm6W9JKk02qc/6CkmZLmSrpX0nurzh0jaXFuJCvbRjXy2FmSJX2rSV02kHSlpBcl3Vf9WUpaS9IFkmbnvDYo8N4+LOkRSfMlnS9plQHv6y+5rKua5RXCYIoGJCyrZgPfAn458ISkdYDfAIcBk4EvAWdIWr0q2e9yI1nZHhyQx1jgeOCGAnU5E7gNWBU4EvgfSavlc/3AJcAHirwpSa8DTgb2B9YAXgROqEryLPAj4Lgi+YUwmKIBGeUk/RpYH7gw/yX+5fwXtSUtl9NcJelb+S/feZIulLSqpBmS5ki6qfova0mvkXS5pGcl3S/pg52ut+1zbZ8PPFPj9LrA87b/6OQPwHxg4xaK+CJwGXBfo0SSNgG2Bo62vcD2OcBd5AbD9hO2TwBuKljuR4ALbV9jex7wdeD9kibl/P5k+yxSAxrCsIoGZJSzvT/wf8B78l/i36uTdF/SX8XrkL6I/wqcCqwCzASOBpC0InA5cAawOrAfcEL+y3opkk6Q9Hyd7c4239bNwExJe0rqyY+vXgKq83tPbuDukfSZAXV6NfAJ4NgCZb0OeND23Kpjd+Tj7Xhdvh4A2/8AFgGbtJlfCINmueGuQFhmnJq/zJD0R2Az23/K+2cD38zp3g08bPvUvH+rpHOAvYF7BmZq+7PAZztZUdt9kn5FasTGk76A97E9Pyc5C5gOPAFsB5wj6XnbZ+bzPwa+bnuepGbFTQReGHDsBVJD2456+U1qM78QBk3cgYSinqh6vaDG/sT8+tXAdtV3EqTHMmsOTTUhd2J/D3gLMA7YGfiFpC0BbN9re7btPtt/IfV17J2vfQ8wyfbv6uR9T1XH+07APFI/S7XJwNylr14qr52q8qo0rm3nF8JQizuQANDJJZkfBa62/fYiiSWdBHy0zulHbLfzKGhL4BrbN+f9myTdAOwK3F4jvYHKrcbbgGmS/pn3VwL6JL3B9l4D65P7QDaSNKnqMdYWpLufhmz/mSUNb8U9+fpK/hsBywN/a5ZfCEMt7kACpLuJpYaxtukiYBNJ+0sam7c3SnptrcS2DxowGqp6q9t4SFpO0nigB+iRNL7S6U/qsN6pcschaStgJ3IfiKS9JL1KybbA54Hf52u/Tupv2DJvFwA/Bz5ep/5/IzVKR+c6vA/YHDinqq7jSY0AwPJ5v54ZpP6ZnXJ/0rHAuZXGKffpjCf98Tcmlzm2QX4hDJpoQALAd4Cj8iOnw8tklL/o3kHqdJ8N/BP4Lku+QDvlKNKjs6+S7mAW5GPYvho4hjScdi7py/w/bV+Wr90XeID0WOhXwHdtn16pv+1/Vrac73zbzzaoy77ANOA50vDavW0/VXV+AenRFKRRXQvqZWT7HuAgUkPyJKnvo7qPaP98/YmkRnEBqYELYcgpAkqFEEJoR9yBhBBCaEs0ICGEENoSDUgIIYS2RAMSQgihLSNqHkjPpBW93Gorl8ukv+nM48Zc8npIy+91AZUdX9El4zNKvw8o/W/SDXUAGNOB9ZXH9JZ7M1pcvhJe+FLpPOby3NO2V2uesr537rKin3m22Pu55c6XLrW9W5nyus2IakCWW21l1vn2waXy6F9Y8iNZUP6mbsyi8nmot3QWjFlcrjFUB76s1IkvzcUdyGNRuet7Sl4PsNyL5fNYfk75D3SFp8t9oONml59U3zez/LzKP/l/HimbxzPP9nHjpesXStuz1t+nlC2v24yoBiSEEIaSgf5ueWQwDKIBCSGENhmzuHvjrg26aEBCCKGEuAMJIYTQMmP6RvFqHtGAhBBCCf3dMtxwGAzqPBBJ60m6UtLMHEfh0Hx8lRzy9O/556vqXL9bDon6gKSvDmZdQwihVQb6cKFtJBrsiYS9wBdtvxbYHjhY0makFVSvsD0VuCLvv4KkHuBnwO7AZsB++doQQuga/bjQNhINagNi+3Hbt+bXc0mxs9cB9gJOz8lOB95b4/JtgQdsP2h7EfDbfF0IIXQFA4vtQttINGR9IJI2ALYCbgDWsP04pEZG0uo1LlmHFN2uYhYpfvXAfA8EDgTombJSZysdQggNeAQ/nipiSNbCkjSRFNTnC7bnFL2sxrGl/qVsT7c9zfa0nkkrlqlmCCG0xtBXcBuJBr0ByeE2zwFm2D43H35C0lr5/FqkyGsDzQLWq9pflxThLoQQukKaiV5sG4kGexSWgFOAmbb/u+rUBcAB+fUBLIlHXe0mYKqkDSWNI4UNvWAw6xtCCK0RfQW3kWiw70B2JMVwfquk2/O2Bylu9Nsl/R14e95H0tqSLgaw3QscAlxK6nw/K8eLDiGErmCg38W2kWhQO9FtX0vtvgyAt9VIPxvYo2r/YuDiwaldCCGUY2DRKA6rFDPRQwihhP5OxABaRkUDEkIIbUoz0aMBCSGE0CIj+jr0CEvSesCvgDVJA7em2z5e0u+ATXOylYHnbW/ZkUJLigYkhBBK6OAjrMrST7dKmgTcIuly2x+qJJD0X8ALnSqwrGhAQgihTUYsck9n8kqrc1RW6JgrqbL0073w8rSIDwJv7UiBHTCiGpCeHjN58oJSeSxcfmyp6xeNL/+R9i/uwC/kS+Vvq/t6S8ZEL3k9gPrK59GRmOgvlauHF5avQycetY/pwL9Jz+Ryv+PqnVi6DuNeXK95omYeLp9FmkhY+P/aFEk3V+1Ptz29VsIBSz9V7AQ8Yfvvrdd0cIyoBiSEEIZaC53oT9ue1ixRg6Wf9gPObL2GgycakBBCaJMt+ty5eSB1ln5C0nLA+4FtOlZYB0QDEkIIJfR3aBhvg6WfAHYF7rM9qyOFdcjonUIZQgglpXkgYwptBdRb+gnSWoBd9fgK4g4khBDaZsRid+ZrtNHST7Y/1pFCOiwakBBCKKEvljIJIYTQqk7ORF8WRQMSQggl9HdwFNayJhqQEEJoU6UTfVkj6f0Fki3MITXqigYkhBDaZLSs9oH8nBQJtlHl/5Um8ZgGtQGR9Evg3cCTtl+fjxVaWVLSw8BcoA/oLTKDM4QQhpJNx0ZhDbE/2v5EowSSftMsk8F+56cBPyUtUQxAiytL7mL76UGrXQghlKKOTSQcSrY/2ok0gx3S9pq8KNhSunFlyRBCaIWho0uZDBVJHwVk+9cDjn8KmG/7jCL5DOc7b7aypIHLJN0i6cB6mUg6UNLNkm7ufWH+oFQ0hBDq6eBM9KH0ReD8Gsd/m88VMpwP75qtLLmj7dmSVgcul3Sf7WsGJsrLIU8HmDB1bQ9OVUMIYWlGy2pM9B7bcwcezHFICse0GJYGpMjKkrZn559PSjoP2BZYqgEJIYThYpbZTvSxkla0/YrHNjkS4riimQzXfVXDlSUlrZjfCJJWBN4B3D2E9QshhAJEX8Gty5wC/E91H3V+/dt8rpDBHsZ7JvAWUiSuWcDRtk+hxsqSktYGfmF7D2AN4LzUz85ywBm2LxnMuoYQQqvMsjkT3fYPJM0Drs4BrAzMB46zfWLRfAZ7FNZ+dY5/rMax2cAe+fWDwBatljd2TB9rTFzqsV5L5o1bvtT1C8aXC4kLsHBR+X+WRR3Io/elcnm4t/x/LC8u/5db/+Ly9RhT8uP0ch14Hx3432p1oh7lQi73ji9fB49ZvXQenQhpCy1FJOwqtk8CTsoNiGr1iUg6wPbp9fJYJh/ehRBCN7DVVXcgkqYC3wE2A8ZXjtveqN41tuc1yPJQIBqQEEIYDF02D+RU4Gjgh8AuwMdpvFxJMw2v7ap3HkIIy5IUUKqn0DZEVrB9BemR1CO2j6HcZO2GUyPiDiSEENqUOtG7qg9koaQxwN8lHQI8BpTpMIo7kBBCGCxdNhP9C8AE4POkeXYfBQ4okd91jU7GHUgIIbSp22ai274JQJJtf7xZeklrAP8JrG17d0mbATvk6RbYPqTR9XEHEkIIJfQzptA2FCTtIOleYGbe30LSCQ0uOQ24FFg77/+NdBdTSDQgIYTQJhv6rELbEPkR8E7gmVQ/30EKDFXPFNtnAf05fS8pBlMhDR9hSfpxgTzm2D6qaIEhhDBSGNHbP2QjrAqx/aheOWG0UYMwX9Kq5NFWkrancYymV2jWB7IX8I0mab4KRAMSQhiVumwm+qOS3gRY0jhSZ/rMBukPAy4ANpZ0HbAasHfRwpo1ID9sNI0dQNKrihYWQggjSRcO4z0IOB5YB5gFXAYcXC+x7Vsl7UwKMy7gftuLixbWsAGx/aNmGRRJE0III1PnljKRtB4p/PeapD6J6baPz+c+BxwC9AJ/sP3lWnnkEOAfaaHMCaS7kFfb/pSkqZI2tX1RkesLDeOVtBrwKWCD6muaBWUPIYSRroMx0XuBL+a7gknALZIuJ61Ovhewue2XcpC9miR9D/gWsAC4hLQo7Rds/6bOJacCtwA75P1ZwNlA5xoQ4PfAn4E/0UIPfQghjGQ2LO5QJ7rtx4HH8+u5kmaSHkV9irTM+kv53JMNsnmH7S9Leh+pMdgHuBKo14BsbPtDkvbLeS+Qii/ZXLQBmWD7K0UzDSGE0aDFiYRTJN1ctT89h+ReSg7utBVwA/B9YCdJ3wYWAodXJgzWUIknsQdwpu1nm7QHiyStwJJRWBsDLxV7O8UbkIsk7WH74qIZhxDCaNDCI6ynbU9rlijH5ziH9OhpTg4B/ipge+CNwFmSNrJda6HDCyXdR3qE9dnc/bCwQXFHkx51rSdpBrAj8LGib6hh74+kuZLmkNaEv0jSAklzqo43JOmXkp6UdHfVsWMkPSbp9rztUefa3STdL+kBSV8t+oZCCGGoVEZhFdmKkDSW1HjMsH1uPjwLONfJjaQO9ik162N/ldSfMS2PpnqR1H9SqywB9wHvJzUaZ+brripUWZqPwppUNKM6TgN+ShpZUO2Htn9Q7yJJPcDPgLeTPrybJF1g+96S9QkhhI7q4CgskeKRz7T931WnzictyX6VpE2AccDT9fKx/VzV6/mkULW10lnS+ba3Af7QTp0LvXNJVxQ5NpDta4Bn26jXtsADth+0vYgU6L1mKxpCCMOm4N1HwTuQHYH9gbcOeELzS2Cj/CTnt8ABdR5fteN6SW9s9+JmS5mMB1Ykdf68iiVrw09myeJb7ThE0r8BN5OGrT034Pw6wKNV+7OA7erU8UDgQIAV11yRtScUnoVf07yx45snamDO4nIx1QHmjxtXOo8XF5XPY8HYcvHdFy8uPzqlvwNx1ft7y9ejb1G5evQvX36oZ9/48p9FXwfq0Tuh3PWLJ5Z/H73jy/9+d4KB3g7dgdi+lvrxNz7akUKWtgvwaUmPkO5UlKrizYtc3KwT/dOklRnXJo0Vrry5OaRHTO04Efgm6bP/JvBfwMD5JLU+xJotbh7FMB1gymundKpVDiGEprplJrqkrRudt31rnVO7lym3WR/I8ZJ+CnzN9jfLFFSV5xOV15J+Tu0JK7OA9ar21wVmd6L8EELopG5oQEh/iAOMB6YBd5D+EN+cNBT4zXWum1vwWE1N771s95HGFHeEpLWqdt8H3F0j2U3AVEkb5gXB9iUt+BVCCF2jMg+kU6Ow2q6HvYvtXYBHgK1tT8ud41sBDzS49FbgKVIckL/n1w9JulXSNs3KLfrw7jJJH2hlhiKApDOBvwKbSpol6ZPA9yTdJelO0vO3f89p15Z0Mby8Jv0hpEAnM4GzbN/TStkhhDAU+lGhbYi8xvZdlR3bdwNbNkh/CbCH7Sm2VyU90joL+CzQKBAVUHwi4WGkzvQ+SQtY0tEyudFFtvercfiUOmlnU3WnkyctxsTFEEL3ctc8wqqYKekXpKVLTOp8b7Sc+zTbB1V2bF8m6T9tHyap6YigQg1IB+aDhBDCiGOgt7+rArt+HPgMafI3wDWkgUv1PCvpK6ThwQAfAp7Lc/H6mxVW9A4ESXuyJDTiVUWX+w0hhJGqxbWwBp3thZJOAi62fX+BSz5MWs7k/Lx/bT7WA3yw2cVFl3M/jrQGy4x86FBJb87T5kMIYdRyFzUg+Q/975Nmq28oaUvgWNt71kqf44d8TtJE2/MGnG7U+Q4U70TfA3i77V/a/iWwGx0cmRVCCMuqLutEP5q0ksfzALZvJ8VxqknSmyTdC9yb97eQ1LTzvKKVh3crV71eqYXrQghhRLI7u5hiB/TabmU5jh8C7wSeAbB9B0u6Kpoq2gfyHeA2SVeSRmD9K3BEC5UMIYQRSPR1Vyf63ZI+DPRImgp8HvhLowtsPzpghkbhoIGF3rntM0lr0Z+btx1s/7bxVSGEMPLZKrQNkc8BryMFhToDeIElI7JqeVTSmwBLGifpcBoP+32FwqOwSI3N0/maTSRtklfbDSGEUalb1sKq8i7bRwJHVg5I2ocU57yWg4DjSQvYzgIuAw4uWljRUVjfJY0PvoclY4NNGmMcQgijk1M/SBc5gqUbi1rHgJdHYX2k3cKK3oG8F9i0EtQ9hBBCMoQjrOqStDtpZOw6kn5cdWoy0Fsj/U+os8I5gO3PFym3aO/PgywJ1h5CCIH0DdwlfSCzSfGVFpJCb1S2C0ijrAa6OZ8fD2xNWkjx76R1swp3ohe9A3kRuD1HIXz5LqRoKxVCCCOT6Osf/juQPPz2DknnAfPzKuqV8OBLrWll+/R8/mPALjl+OnkW+2VFyy3agFxALKceQghL6aaZ6KQv/12ByqzyFfKxN9VJvzYwiSWhxyfSQrTZoospnt7ovKRzbH+gaKEhhDAS2F3XgIyvXpLE9jxJjYIQH8eSOX4AOwPHFC2slWG8jWzUoXxKWX5MLxuu8EypPJ4bWy7g87OLVix1PcC8seXjqs8bWz5mdNnY7C8tLv/rtbi/A/HM+8pP9FpUMq5676Lyn0Xfog7EmF/YgXjkJWOaq6/8F+7CKd3zpd1lw3jnS9q6EsI2B4VaUC+x7VMl/RHYLh/6qu1/Fi2sUw1Izd58Sb8E3g08afv1+dj3gfcAi4B/AB+3/XyNax8mhVbsI03Pn9ahuoYQQsd02TDeLwBnS6qEAF+LNAXjFSStWWko8s/fN0pTz2DPwT+NtPBitcuB19venBRGsdGSKLvY3jIajxBCNzKiv39MoW1I6mPfBLyGFBPks8Brbd9SI2mRYH1N03TqDqTmPZztayRtMOBYdQ//9cDeHapDCCEMuW66AZH0bwMObSUJ278acHwLSXMaZQU0Og+00YBIehWwnu07qw5/pdV8sk8Av6tzzqRY7AZOtj29Tn0OBA4EWGmtFdqsRgghtKH7OtHfWPV6PPA24FbgFQ2I7fIdahRfyuQqYM+c/nbgKUlX2z4sV6bwuOGqPI8kzZCcUSfJjrZnS1oduFzSfbXW3soNy3SAdV63cjf9MRBCGA266FvH9ueq9yWtBPx6sMor+mBuJdtzgPcDp9rehjTWuC2SDiB1rn/Ert0FZXt2/vkkcB4pSEoIIXSVTs1El7SepCslzZR0j6RD8/FjJD0m6fa8tRLM70Vgaptvramij7CWk7QWKUbukc0SNyJpN9Ijr51tv1gnzYrAGNtz8+t3AMeWKTeEEAZDB0dh9QJftH2rpEnALZIuz+d+aPsHzTKQdCFL7onGAJsBZ3WshgMUbUCOBS4FrrN9k6SNSOumNCTpTOAtwBRJs0jhFo8gTa2/PAcxud72QZLWBn5hew9gDeC8fH454Azbl7T0zkIIYZDZ4A6NsLL9OPB4fj1X0kzSMuutqG5keoFHbM9qdIGkNwNT85yQ1YCJth8qUljRmehnU7UcsO0HgaYzz23vV+PwKXXSzibHWc/5b1GkbiGEMJwGYx5IHr26FXADsCNwSB5hdTPpLuW5Gtf0AF+3Xbh7QdLRwDRgU+BU0qK5v8llNlWo6ZS0iaQrJN2d9zeXdFTRSoYQwojlglt6EnNz1XZgrewkTQTOAb6Q+55PBDYmrZT7OPBfNauRFlB8MXecF/U+0gCp+TmP2aS1sQop+gjr58CXgJNzIXdKOgP4VgsVDSGEEaalpdqfbjYpWtJYUuMxw/a5ALafqDr/c+CiBlksBO7KfSfzKwcbrJy+yLbzdIlK/3NhRRuQCbZvHBB4fakgJSGEMOp06BGW0hfsKcBM2/9ddXyt3D8C6Y7h7gbZ/CFvRWt4lqSTgZUlfYo0N+/nRetctAF5WtLGlYpI2pvc2RNCCKNWZycS7gjsT7qDuD0f+xqwn6QtU2k8DHy6QR4r2z6++kBlOHAttn8g6e2kWeebAt+wfXm99AMVbUAOJk3We42kx4CHgI8WLSSEEEasDjUgtq+l9rJQRdatqjgAOH7AsY/VOFbpdL80d7oXbjSqFR2F9SCwa/X8jHYKCyGEEacLZqJL2g/4MLChpOrgf5OBmjEubPdJelHSSrZfaKfcokuZrAH8J7C27d0lbQbsYLvmkNwQQhg1uqABAf5C6laYwitHac0F7qx5RdJqp/srFH2EdRppjHBlFvrfSIsgRgMSQhi9TMceYZWqhv0I8IikXYEFtvslbUJa2v2uBpfW6nQvrGgDMsX2WZKOyJXtldTXbqEhhDBSdFlAqWuAnfKq6VeQJh5+CPhIrcTNwpU3U7QBmS9pVZaMwtoeaOuZWQghjCj9w38HUkW2X5T0SeAntr8n6ba6iaWHqPEQznahMOVFG5DDgAuAjSVdB6xGFwaCGjeml1cv/3SpPCb1TCx1/cSel0pdDzCnd3zpPBaUjGcOML+vXB6L+srHK+v10ERya2ZRX7nwCXMXl49zv2DR2NJ5LFxcPo++kl+Y6qrv2/LUXXcgkrQD6Y7jk/lYo/+I1RMbxwP7AKsULazp/3BJY3LGO5PGCQu43/biooWEEMKItGSZkm5xKGnB2vNs35MXvr2yXmLbA0do/UjStcA3ihTWtAHJnTH/ZXsH4J4imYYQwuigruhEr8hB966p2n8QqDuiStLWVbtjSHckHV8L6zJJHwDOrRcAKoQQRqUu+kbMI68OBzag6vvd9lvrXFI95LeXNEn8g0XLa6UPZEWgV9JC0mMs255ctKAQQhiRuqgBIYXdOAn4BVBkpOwn813KyyRtWLSwojPRC9/ShBDCqGG6bRRWr+0TW0j/P8DWNY5tU+TiojPRBxYAaRjvI7brrsor6Zek2OdP2n59PrYKaRLiBqSFwT5YJzjKbqT1W3pIkQqPK1LXEEIYSl02CutCSZ8FzgNeHhJq+9nqRJJeA7wOWEnS+6tOTSYNmiqk6COsE0itVGVG4xuAO4BVJR1k+7I6150G/BT4VdWxrwJX2D5O0lfz/leqL8qLfP0MeDswC7hJ0gW27y1Y3xBCGBrd1YAckH9+qeqYgYHzOjYl/XG/MvCequNzgU8VLaxoA/Iw6VnZPQB5LawvAd8EzgVqNiC2r8mhGavtRYqTDnA6cBUDGhBgW+CByrM5Sb/N10UDEkIIddgu1H9h+/fA7yXtYPuv7ZZXtAF5TaXxyIXfK2kr2w+q9VlBa1SCo9h+XNLqNdKsAzxatT8L2K5WZjks5IEAq65dfvJcCCG0opseYeWIhp8B/jUfugo4ucG8vdskHUx6nPXyoyvbnyhSXtFpvvdLOlHSznk7AfibpOWBwZhQWKtVqvnPZHu67Wm2p01apfws2xBCaIlVbBsaJ5I6wE/I2zb5WD2/BtYE3glcDaxLeoxVSNE7kI8BnwW+QPpyv5Y01ngxsEvRwrInKiEaJa0FPFkjzSxgvar9dYHZLZYTQgiDy0D/cFfiFd5oe4uq/f+VdEeD9P9iex9Je9k+XdIZwKVFCys6jHdBvuu4yPb9A07PK1pYdgGpo+e4/PP3NdLcBEzN45EfA/YlBUsJIYSu0k2PsIA+SRvb/gdAXsqk0XyQyhOk5yW9HvgnaYRsIYUeYUnaE7gduCTvbzkg6lW9684E/gpsKmlWXiHyOODtkv5OGmV1XE67tqSLIS0XDxxCaglnAmdV98GEEELXcMFtaHwJuFLSVZKuBv4X+GKD9NPz0u9Hkf64vxf4btHCij7COpo0MuoqANu31xhdtRTb+9U59bYaaWcDe1TtX0xrsYBDCGHoddEdiO0rJE1lycK399muuUR4Xih3Tp6Hdw1LD/Vtqmgnem+7MXNDCGGkkotvQ1KfNKJqBdt32r4DmJAnFi7Fdj/pSU/bijYgd0v6MNAjaaqkn5Bi8IYQwujWr2Lb0PiU7ecrO/nuotHEwMslHS5pPUmrVLaihRV9hPU5Ujz0l4AzSX0T3yxaSAghjFRd1ok+RpIqq6bnVT0aTZCrzPc4uOpYrZnrNRUdhfUiqQE5skj6EEIYNbqrAbkUOEvSSaSaHUQe/FRL0Znr9TRsQCRdSIOPx/aeZQoPIYRl2hD2bxT0FeDTpNnoIi0z9Yt6iSVNIIXrWN/2gZUOeNsXFSms2R3ID/LP95NmK/4m7+9HWh+rq4yljzWXe755wgYmjCkX03ylnhdLXQ/w4tjyS7IsdPlZ+S/1l8tjscvFEe+UMR34E7Hse3mpv3x8+Hl95eOqL+pAPcpasafc/zGAiR3Io2ML63VRA5IjyJ4G/G+NOXu1nArcArwp788ixRQp1IA07ES3fbXtq4GtbH/I9oV5+zDw5iIFhBDCiNaheSC5I/tKSTMl3SPp0AHnD5dkSVMa5NHqnL2NbX+PPKHQ9gJqLyVVU9FRWKvlGY2VSm4IrFa0kBBCGKk6OIy3F/ii7dcC2wMH55XPkbQeaeL1/zXJozJn73lIc/ZoPLN8kaQVyE2cpI2piiPSTNH72X8HrsTlcvoAABsuSURBVJJUCX24AXkF3BBCGNU69Agrr1JeWal8rqSZpJXJ7wV+CHyZ2ks/Veu1/UILq6QfTbpbWU/SDGBH0tqHhRQdhXVJ7lx5TT5Ud3ZjCCGMGq11ok+RdHPV/nTb02slzCt9bAXckB9LPWb7jgINwyvm7AGfp8GcPduXS7qVdMcj4FDbTxd9Q81GYW1t+9Zc0EukKIR104QQwqhTvAF52va0ZokkTQTOIa1+3kuaPvGOgmVUz9mrrKz7rSbX7Ezq0zYwlhQOt5BmdyCnSnoLjTtVTiG1lCGEMPp0cBRWDgh1DjDD9rmS3gBsCFTuPtYFbpW0re1/Dri2B/gP21+i4Jy9vMr6v5AmiAN8WtKutg9ucNnLmjUgK5GGeDVqQJ4qUlAIIYw0onPzQJRaiFOAmbb/G8D2XcDqVWkeBqbVesxku0/SNi0WuzPw+qqZ66cDdxW9uGEDYnuDFisTQgijh0GdCyi1I7A/cJek2/Oxr+WVyYu6LQ/bPRuY/3I17XPrpL8fWB94JO+vB9xZtLDhn1UUQgjLss6NwrqWJnMwCvxRvwrwDPDW6suAeg3IqsBMSTfm/TcCf63MHWm22kg0ICGEUEZ3zUT/eIuXfKNMecPSgEjaFPhd1aGNgG/Y/lFVmreQxjw/lA+da/vYIatkCCEU0GVrYbUkrzSCpMlUtQe2ny1yfaEGJHfufATYyPaxktYH1rR9Y5NLa8prtGyZ8+4hxT2vNXTsz7bf3U4ZIYQwJJbhBkTSgaTQHAuAftIjtM4u5w6ckDN/K3AsMJc01OyNLda3lrcB/7D9SNOUIYTQTTrbiT4cvgS8rpXJg9WKNiDb2d5a0m2QolxJKr9kbLIvS8YgD7SDpDuA2cDhtu8ZmCC3oAcCrL52dOmEEIZYF9yBSDqs0fnKsOAa/gG0vYR40W/cxflRU2Ws8GqkO5JSciO0J3BEjdO3Aq+2PU/SHsD5wNSBifJSANMBNnnDCl3wTxlCGE26pA9kUv65KenJUGUF3vcA1zS47gjgL5JuoGoRRdufL1Jo0Qbkx6Q+itUlfRvYGziq4LWN7A7cavuJgSdsz6l6fbGkEyRNafdWK4QQBkUXNCC2/wNA0mXA1rbn5v1jSHNC6jkZ+F/S5MGWbwqKLqY4Q9ItpP4KAe+1PbPVwmrYjzqPryStCTxh25K2JS09/0wHygwhhM4oGOtjCK0PLKraX0Tj5dx7bTd8/NVIs8UUV6nafZKqL3tJqxQd6lUn7wmk9e0/XXXsIADbJ5Hucj4jqZc0QmDfynT7EELoBqKF6EtD49fAjZLOIzVt7wN+1SD9lbkf+UJe+QirI8N4b8mVEKlley6/XpkU2KTtgOy2XyTNgqw+dlLV658CP20lzzEyE7SoecIGxvb0lbp+vBaXur5TejowNKTPReONhWb6C8duq6/P5b+qFrv8QJMxJX+3Vh5TPuzzmsvNLZ3Hf5bOIemmUVi2vy3pj8BO+dDHbd/W4JIP55/V/dCdGcZre0MASScBF1TWZJG0O7BrkQJCCGFE677nIhOAObZPlbSapA1tP1QrYeU7vl1F/yx6Y/WCXrb/SFrFMYQQRrcOxUTvBElHA19hyR3FWOA3DdJPkHSUpOl5f6qkwpO3izYgT+dCNpD0aklHEh3aIYTRrmA89CEc6vs+0tSI+QC2Z7NkiG8tp5I62t+U92fRPADVy4o2IPsBq5GG8p5PWp9+v6KFhBDCiNVFdyDAojzYqDJnb8Um6Te2/T1gMYDtBbQwLqDoMN5ngUOLZhpCCKNFN3WiA2dJOhlYWdKngE8AP2+QfpGkFVjS4GxM1WisZooupnglNdpQ22+tkTyEEEaNLpmJDoDtH0h6OzCHNCv9G7Yvb3DJMcAlwHqSZpCCWhVeEr7omL7Dq16PBz5ACvYeQgijV/dNJCQ3GI0ajeq0l+VJ4tuTHl0d2spqH0UfYd0y4NB1kq4uWkgIIYxYXdCASJpLg5rYnlznuitsvw34Q41jTRV9hFU9I30MsA2wZpFrQwhhpBLd8QjL9iQASccC/yTNSK/EcVpqFJak8aT5IlMkvYolHeeTgbWLllv0EVb1jPReUpTATxYtJIQQRqwuaECqvNP2dlX7J+aVdr83IN2ngS+QGotbWNKAzAF+VrSwog3Ia20vrD4gafmihYQQwohkUH9XtSB9kj4C/JbUtO0HLLU+k+3jgeMlfc72T9otrOg8kL/UOPbXdgsNIYSRossmEn4Y+CDwRN72Ycl6V0sp03hA89V41wTWAVaQtBWvfE42oUzBIYQwInTRDYjth4G9hqq8Zo+w3gl8DFgXqA6JOBf42iDVKYQQlhnd0IlekaPFfooUA+Tl73fbnxiM8pqtxns6cLqkD9g+ZzAqEEIIy7QuakCA3wN/Bv5Ejb6PgSRVRmptZPtYSesDa9q+sUhhzR5hfdT2b4ANagVtbxCovSlJD5PuZPpIUbGmDTgv4HhgD1LQ94/ZvrXd8kIIoeM62L8haT1S8Kc1SeFlp9s+XtI3SY+l+kmB/T6WF0msZYLtr7RQ7Ak537cCx5K+k88hxVVvqlknemUhromkscTV28QWKlnPLra3HNh4ZLsDU/N2IHBiB8oLIYSOEWktrCJbAb3AF22/ljQz/GBJmwHft7257S2Bi4BvNMjjIkl7tPAWtrN9MLAQwPZzwLiiFzd7hHVyfvkn29dVn5O0YwuVbMdewK/yypLXS1pZ0lq2Hx/kckMIobgORdrO322P59dzJc0E1rF9b1WyFWn80OxQ4GuSFpGWaVfKrvZMdGCxpJ5KnrkPpfDykEWH8dYa6lVq+BepwpdJuiXH5B1oHeDRqv1Z+VgIIXSNFobxTpF0c9VW63sv5SltAGwF3JD3vy3pUVJ/Rd07ENuTbI+xPd725Lxfr/EA+DEpTMfqkr4NXEsL0X6b9YHsQAo0stqAPpDJQE/RQurY0fZsSasDl0u6z/Y11cXXuGapljf/IxwIsNY6PUweU3gl4kExaZjLrxjXvP+sqTElH+6OK/6HTNfrKx4ioaaeLulpLfs+AMaW/L2YVL4KTOnpklkErS2m+HSdx/WvIGkiqR/iC7bnANg+EjhS0hHAIcDRda6tdIpvaPubuV9lrXqd4rZn5MUU30b6zn2v7ZlF31CzO5BxpL6O5Xhl/8ccYO+ihdRS6QSy/SSpBdx2QJJZwHpV++sCS3Uc2Z5ue5rtaa9apegNVQghdEYH+0CQNJbUeMywfW6NJGeQVkOv5wRgB5ZMHpxHg6VJJG0PPGb7Z7Z/CsyStF299AM16wO5Grha0mm2HymaaTM5StaY/JxvReAdpBEA1S4ADpH0W2A74IXo/wghdJtOBZTKdw+nADOrR7hKmmr773l3T+C+BtlsZ3trSbdB6hSX1KhT/ERg66r9+TWO1VV0LawXJX0feB0pHgi5cu0GlFoDOC99XiwHnGH7EkkH5XxPAi4mDeF9gDSMt3CQkxBCGBKmY53opGBO+wN3Sbo9H/sa8ElJm5I6tx8BDmqQR6ud4soDlQCw3S+paLtQuAGZAfwOeDep8gcATxUtZCDbDwJb1Dh+UtVrAwe3W0YIIQyFTs0DsX0ttft+L24hm4Gd4nsDX2+Q/kFJn2fJNInPAg8WLaxop8Gqtk8BFtu+Ok+L375oISGEMGK54DYUVbFnAF8GvkMaEvxe22c1uOQg0kCpx0j9ztuRByUVUfQOZHH++bikd5E6s9ctWkgIIYxE3RJQqkLSr23vT1U/SdWxpeRBTPu2W17RBuRbklYCvkia/zGZFIwkhBBGL7uTfSCd8Lrqndwfsk29xDky4SdZun+70OKLhR5h2b7I9gu277a9i+1tgI2LXBtCCCNZJ4fxtl0H6YgcF31zSXPyNpe0dtYFDS79NWntrXcCV5OeLM0tWm6ZiRNLLa4YQgijTTcElLL9nRwX/ft5BnplFvqqtr/a4NJ/sf11YH5eff1dwBuKllumAenAfNIQQliGGeh3sW1oPFC9I6lHUs1Z61mlf/t5Sa8HViLFEimkTAPSVQ/+QghhWHTRKCzgbZIulrSWpDcA15NWD6lnuqRXAUeRHnXdC3y3aGHN1sKaS+23LmCFooWEEMJI1U2jsGx/WNKHgLtIE7D3G7iSOoCkQ20fT5r1/hxwDbBRq+U1vAOprORYY5tku/BsxRBCGLEqI7GabUNA0lTSku7nAA8D+0uqtfJkZWWPUquqRyMQQgjt8uCPsGrRhcDBtq/Ia2sdBtzEgOG9wMwcFXZ1SXdWHa/ED9m8SGHRgIQQQpvSRMIueoYF21YtAW/gvyQtNYzX9n6S1gQuJS3Q2JZY/zyEEMroL7gNIklfBrA9R9I+A07XW4j2KeAu248M3IqWGw1ICCGUILvQNsiqlyM5YsC53WpdYLuPFCWxcAz0geIRVgghtGtoh+g2ojqva+1XewS4Lj/mml85WB2PpJFoQEIIoW1GQzdJsBHXeV1rv9rsvI2h8XyRmkZUA7LCuM3ZYv2bS+XxzOxyiwy/6PIx0Re6Oyb5jy35p1VPB95GT5cseNBX8iH24g58x3QinnlfB+qx0GWffJfvEFjcN795oqHSHZ3oW0iaQ56jl1+T98fXu8j2f5QpdFgakBzo/VekRbz6gel5Ukt1mrcAvwceyofOtT0w7G0IIQyfLhnGa7unneskXUmNO5Si0WaH6w6kF/ii7VslTQJukXS57XsHpPuz7XcPQ/1CCKGY7rgDadfhVa/HAx8gfT8XMiwNiO3HSdGysD1X0kxgHdI6LCGEsOxYhtsP27cMOHSdpKuLXj/sfSCSNgC2Am6ocXoHSXeQOnkOt31PjesPJIdgXH/99QevoiGEUEOXTSRsiaRVqnbHkIJPrVn0+mFtQCRNJK3Z8oXK7MkqtwKvtj1P0h7A+cDUgXnYng5MB5g2bdqy+y8ZQlj2mM6MTBg+t5DehUiPrh4iRSgsZNgaEEljSY3HDNvnDjxf3aDYvljSCZKm2H56KOsZQgj1iCGZJDhobG9Y5vphmYmeF/k6hbSUcM0JK5LWzOmQtC2prs8MXS1DCKGALlqNt1WS9skDmZB0lKRzJW1d9PrhugPZEdgfuEvS7fnY14D1AWyfBOwNfEZSL7AA2DcvDhZCCN1j2f5a+rrtsyW9mRQX/QfAicB2RS4erlFY19IkJK7tnwI/HZoahRBCG8ygL5Q4yPryz3cBJ9r+vaRjil4ciymGEEIJnVpMUdJ6kq6UNFPSPZIOzce/L+k+SXdKOk/Syh2s/mOSTgY+CFwsaXlaaBeiAQkhhLYZ+vuLbc1VJli/FtgeOFjSZsDlwOtzkKe/sfRqu2V8kBQTZDfbzwOrAF8qevGwzwMJIYRllulYH0i9Cda2L6tKdj2pf7hTZb4InFu1/3IdiogGJIQQyijeBzJFUvVqr9PzPLalNJhg/Qngd61VcPBEAxJCCCW0MA/kadvTmuZXZ4K1pCNJj7lmtFPPwRANSAghlNHBYbz1JlhLOgB4N/C2bprOEA1ICCG0y4a+zozjrTfBWtJuwFeAnXOfRdeIBiSEEMro3A1BvQnWPwaWBy7Pi3Ncb/ugThVaRjQgIYRQRudGYdWbYH1xRwoYBNGAhBBCuwx0R0z0YRENyACvGjOh1PXjOxAT/dn+RaXz6MQK06Wf7Hbi/5W64z9n2Zjm89uLOPoKc/vHdSCPuuGxC5vfv3yp6/s6MH+5v3RcdoB/diAPg5fttUzKiAYkhBDK6J5BUUMuGpAQQmiX6dgorGVRNCAhhFBG3IGEEEJoXfcGixoK0YCEEEK7TNGVdkekYVvOXdJuku6X9ICkr9Y4L0k/zufvbCXMYgghDJllOKRtWcMVE70H+BmwO7AZsF9e977a7sDUvB1ICrMYQgjdJRqQIbct8IDtB20vAn4L7DUgzV7Ar5xcD6wsaa2hrmgIIdRl476+QttINFwNyDrAo1X7s/KxVtMg6UBJN0u6+amnnup4RUMIoaF+F9tGoOFqQGqt9zLwEy6SBtvTbU+zPW211VbrSOVCCKGwUfwIa7hGYc0C1qvaXxeY3UaaEEIYPnaMwhoGNwFTJW0oaRywL3DBgDQXAP+WR2NtD7yQ4/WGEEL3iDuQoWW7V9IhwKVAD/BL2/dIOiifP4m0hPEewAPAi8DHh6OuIYRQn0dsB3kRwzaR0PbFDFjnPjccldcGDh7qeoUQQmGxnHsIIYS2xXLuIYQQWmXAcQcSQgihZY6AUiGEENo0mu9A5BE0vEzSU8AjTZJNAZ4eguqUFfXsnGWhjhD17LRm9Xy17VKzjyVdkssp4mnbu5Upr9uMqAakCEk325423PVoJurZOctCHSHq2WnLSj2XZcO2nHsIIYRlWzQgIYQQ2jIaG5Dpw12BgqKenbMs1BGinp22rNRzmTXq+kBCCCF0xmi8AwkhhNAB0YCEEEJoy4hqQCSNl3SjpDsk3SPpP/LxffJ+v6RpA645QtIDku6X9M5urKekDSQtkHR73k6qn/uQ1PP7ku6TdKek8yStXHVNN32eNes5HJ9ngzp+M9fvdkmXSVq76ppu+ixr1rPbfjerzh8uyZKmVB0b8s9zxLM9YjZSFMOJ+fVY4AZge+C1wKbAVcC0qvSbAXcAywMbAv8AerqwnhsAd3fR5/kOYLl8/LvAd7v086xXzyH/PBvUcXJVms8DJ3XpZ1mvnl31u5n31yOFingEmDKcn+dI30bUHYiTeXl3bN5se6bt+2tcshfwW9sv2X6IFHtk2y6s57BoUM/LbPfm49eTokVC932e9eo55BrUcU5VshVZEra52z7LevUcFvXqmfd/CHyZV9ZxWD7PkW5ENSAAknok3Q48CVxu+4YGydcBHq3an5WPDboW6wmwoaTbJF0taachqCJQqJ6fAP6YX3fz51ldTxiGz7NeHSV9W9KjwEeAb+TkXfdZ1qkndNHvpqQ9gcds3zEg+bB9niPZiGtAbPfZ3pL01+a2kl7fILlqZTE4NRtQSGv1fBxY3/ZWwGHAGZImD3c9JR0J9AIzKodqZTH4tWy5nsPyedaro+0jba+X63dIpdq1shjsOrZRz2763dwcOJJXNm4Vw/Z5jmQjrgGpsP08qS+h0eJls0jPSyvWBWYPYrWWUqSe+bb7mfz6FtLz202GpIJL6vCKeko6AHg38BHblf+IXfd51qrncH+eDf7NzwA+kF933WdZ5eV6Dvdnmcut1HMvUv/GHZIeJn1mt0paky74PEeiEdWASFqtaqTNCsCuwH0NLrkA2FfS8pI2BKYCN3ZbPXP6nvx6o1zPB4ernpJ2A74C7Gn7xapLuurzrFfP4fg8G9RxalWyPVnye9Btn2XNenbZ7+Zttle3vYHtDUiNxta2/8kwfZ4j3UiLB7IWcHr+hR4DnGX7IknvA34CrAb8QdLttt9p+x5JZwH3kh5xHGy7r9vqCfwrcKykXqAPOMj2s8NYzwdIo1kulwRwve2DuvDzrFlPhufzrFfHcyRtCvSTRg0dBNCFn2XNetJlv5v1Eg/j5zmixVImIYQQ2jKiHmGFEEIYOtGAhBBCaEs0ICGEENoSDUgIIYS2RAMSQgihLdGAhBBCaEs0IKOIpHnNU5XK/xeSNsuvv9bG9RtIurvF9Avyeki1zh8j6fBW67EskrSypM9W7W+stLz6oP6bh9EtGpDQMbb/n+17827LDUib/pHXQxo0lZnWXW5l4OUGxPagfy4hRAMyykl6taQrlIIFXSFp/Xz8NEk/lvQXSQ9K2jsfHyPpBKUgPhdJurjq3FWSpkk6Dlgh/wU8Y+CdhVKwn2Py622UggL9FTi4Kk2PUkCom3LdPl3w/RypFDDoT6TYKpXjG0u6RNItkv4s6TVVx6/P5Rxb+Ytd0lskXSnpDOCufOyjSkGMbpd0ctUSHu+Q9FdJt0o6W9LEfPw4Sffm+v+gQZ1XyzO9b8rbjvn4tvnzvy3/3DQff11VPe5UWmbkOKBy1/H9Ip9VCKW5C4KSxDY0GzCvxrELgQPy608A5+fXpwFnk/7I2Ax4IB/fG7g4H18TeA7YO5+7ihwIq7osBgQdAg4Hjsmv7wR2zq+/X0kHHAgclV8vD9wMbDig7gPz3Yb0ZT8BmEyK+XB4PncFMDW/3g743/z6ImC//PqgSr2BtwDzK2WSgn1dCIzN+ycA/wZMAa4BVszHv0JaDXYV4H6WrPawcoN/lzOAN+fX6wMz8+vJLAmItStwTn79E9LikADjgBUGfhaN/s1ji61T20hbCyu0bgfg/fn1r4HvVZ0733Y/cK+kNfKxNwNn5+P/lHRluwVLWon0xXp1Vfm759fvADav3N0AK5EWwHuoQZY7Aec5L5wo6YL8cyLwJuDsvCYWpEYJ0vt/b359BlB9p3CjU/AhgLeRGqibch4rkOJQbE9qYK/Lx8cBfwXmAAuBX0j6A6mhqmdXYLOquk2WNCm/59PzHYZJQZPI+R8paV3gXNt/r7o2hCETDUgYqHpxtJeqXmvAz1b08srHpeOr8qq3GJuAz9m+tMWyauU3BnjerfcJzB9Qn9NtH1GdQNJ7SMGM9ht4saRtSQ3PvqT4GW+tU84YYAfbCwZc/xPgStvvk7QB6Q4P22dIugF4F3CppP/HEKyAG8JA0QcS/kL6goMUae7aJumvBT6Q+0LWID3qqWWxpMpfzE8Aq0taVdLypPgcOMVxeEHSm6vKr7gU+EwlD0mbSFqxSd2uAd4naYX8F/x7cjlzgIck7ZPzkqQt8jXXsyQGx74DM6xyBbC3pNVzHqtIenW+fkdJ/5KPT8h1nQisZPti4AtAo8brMpYEaEJSJe1KwGP59ceqzm8EPGj7x6RlyjcH5gKTGpQRQsdFAzK6TJA0q2o7DPg88HFJdwL7A4c2yeMcUpyFu4GTgRuAF2qkmw7cKWmG7cXAsTntRbwy9snHgZ/lTvTqv8B/QVp6+9bcAX8yTe6Ybd8K/A64Pdfzz1WnPwJ8UtIdwD2k4EOQvtwPk3QjaYnwWu8Fp9FlRwGX5c/qcmAt20+RvtzPzMevB15D+jK/KB+7Gvj3BlX/PDAtd4jfy5Kl0r8HfEfSdUD1SLAPAXcrDV9+DfArp6BO10m6OzrRw1CJ5dxDyyRNtD1P0qqkoDw7OgXtGep6bABcZLtROOBmeUwAFti2pH1JHep7NbtuWSFpnu2Jw12PMDJFH0hox0VK0eDGAd8cjsYj6wNWUgq81e6ch22Anyr1Qj9PGom2zJO0Meku7InhrksYueIOJIQhIulIYJ8Bh8+2/e3hqE8IZUUDEkIIoS3RiR5CCKEt0YCEEEJoSzQgIYQQ2hINSAghhLb8f+Nbzrb2DmiJAAAAAElFTkSuQmCC\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": 6,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.collections.QuadMesh at 0x2b339cc94780>"
]
},
"execution_count": 6,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAY0AAAEXCAYAAABRWhj0AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3deZwdVZ338c83gQgEwr5vAScCkQdRM0GHcd+AQVBEh7gjPhlmwN3nkRFn9MGZEXXGGVQQM7KOgBvgIEYBGSEubAHCEhaNQTQGWURJEATS+T5/1Glyudzuru6+t5d7v+/Xq15dt+pU1e8WoX996tQ5R7aJiIioY8p4BxAREZNHkkZERNSWpBEREbUlaURERG1JGhERUVuSRkRE1JakEU8jaRdJD0uaOt6xRMTEkqQRSPqlpFf2f7b9K9sb2+4bz7gGI+lYSYslPSbpzBb73yTpdkmrJd0m6XUN+z4h6YmSGPuX3Vuc4yWSLOmfhohlpqQfSnpE0h2N91LS9pIukrSynGtmje/2Zkl3S/qjpG9L2qLpe/20XOuKoc4V0W5JGjFZrQT+CTi9eYekHYGvAh8EZgD/BzhX0jYNxb5eEmP/srzpHOsDJwHX1IjlPOBGYEvgeOBbkrYu+9YC3wfeUOdLSXo28GXgbcC2wCPAKQ1FHgT+Azixzvki2i1Jo8dJ+i9gF+A75S/u/1v+crak9UqZKyT9U/kL92FJ35G0paRzJK2SdF3jX9CS9pR0maQHJd0p6U3tjtv2Bba/Dfyuxe6dgD/Y/p4r3wX+CDxzGJf4EHApcMdghSQ9C3ge8HHbj9o+H7iFkiRs32v7FOC6mtd9C/Ad24tsPwz8A3CYpE3K+X5g+xtUSTNizCVp9DjbbwN+Bby2/MX9mQGKHkH11++OVL98rwLOALYAbgc+DiBpOnAZcC6wDTAPOKX8Bf00kk6R9IcBlptH+LUWA7dLOkTS1PJo6jGg8XyvLUltqaS/bYppV+BdwAk1rvVsYLnt1Q3bbirbR+LZ5XgAbP8CeBx41gjPF9FW6413ADFpnFF+gSHpe8Bs2z8on78JfLKUOxj4pe0zyucbJJ0PHA4sbT6p7b8D/q6dgdruk3Q2VeLagOqX7htt/7EU+QawALgX2A84X9IfbJ9X9n8e+AfbD0sa6nIbAw81bXuIKrmOxEDn22SE54toq9Q0oq57G9YfbfF547K+K7BfY42B6pHLdmMTJpSG6M8ALwWmAS8BviJpXwDbt9leabvP9k+p2i4OL8e+FtjE9tcHOPfShsbzFwEPU7WbNJoBrH760U8714saztWfUEd8voixkJpGALRzqONfA1faflWdwpJOBd46wO67bY/kMc++wCLbi8vn6yRdA7wSWNKivIH+KsUrgDmSfls+bwr0Sfpftg9tjqe0aewuaZOGR1TPoarlDMr2j1iXbPstLcf3n3934BnAz4Y6X8RYSE0joKo1PO2V0xG6GHiWpLdJWr8sfy5pr1aFbR/d9BZT4zJgwpC0nqQNgKnAVEkb9DfcUzU6v6i/ZiHpucCLKG0akg6VtLkqc4H3Av9djv0HqvaDfctyEfCfwJEDxP8zqkT08RLD64F9gPMbYt2A6hc/wDPK54GcQ9Xe8qLSPnQCcEF/QiptNBtQ/cE3pVxz/UHOF9FWSRoB8CngY+Vx0odHc6Lyy+3VVA3nK4HfAp9m3S/NdvkY1WOx46hqKo+Wbdi+EvgE1auvq6l+gf+L7UvLsUcAy6ge+ZwNfNr2Wf3x2/5t/1LO+0fbDw4SyxHAHOD3VK/CHm77/ob9j1I9doLqbaxHBzqR7aXA0VTJ4z6qtozGNp+3leO/RJUIH6VKahFjQpmEKSIi6kpNIyIiakvSiIiYACSdLuk+SbcOsF+SPi9pmaSbJT2vYd8BpSPtMknHdTLOJI2IiInhTOCAQfYfCMwqy3yqdi1UDSx6ctk/G5gnaXangkzSiIiYAGwvohpbbCCHAmeXoXGuBjaTtD0wF1hme7ntx4GvlbId0VX9NKbpGd6A6eMdRkRMAqv5/QO2tx665MBe87Lp/t2D9QaDvv7mx5YCf2rYtMD2gmFcbkeqflD9VpRtrbbvN4zzDktXJY0NmM5+esV4hxERk8AP/K27R3uO3z3Yx7WX7FKr7NTtf/4n23NGcblWY9p4kO0d0VVJIyJiLBlYy9qxutwKYOeGzztR9YWaNsD2jkibRkTECBnzhPtqLW1wEfD28hbVC4CHbN9DNQLCLEm7SZpG1dn0onZcsJXUNCIiRqFdNQ1J51ENsrmVpBVU0w2sD2D7VGAhcBDVaAaPUIa2sb1G0rHAJVTD6pxeRhboiCSNiIgRMqavTaNq2J43xH4DxwywbyFVUum4JI2IiFFY27k25wmpo20aknaW9ENJt5d5CN5Xtm9RpgP9efm5+QDHj1kvx4iI4TLQh2st3aLTDeFrgA/Z3gt4AXBM6al4HHC57VnA5eXzU4x1L8eIiJFYi2st3aKjScP2PbZvKOurqeaS3pGqt+JZpdhZwOtaHD6mvRwjIobLwBN2raVbjFmbhqSZwHOBa4Bty6ti2L5H0jYtDqnVy1HSfKpxWNiAjdobdETEINxlj57qGJOkIWljqolw3m97ldSqA+PTD2ux7Wn/dUo3/AUAM7RFb/3Xi4jxZejrsd86He/cV6aiPB84x/YFZfO9ZaAtys/7Whw6UO/HiIgJoeoRXm/pFp1+e0rAacDttj/XsOsi4B1l/R2sm5+50Zj2coyIGD7RV3PpFp2uaexPNafxyyUtKctBVPMov0rSz4FXlc9I2kHSQqh6OQL9vRxvB77RyV6OERHDZWCt6y3doqNtGrZ/TOu2CYCnDUdreyVVN/n+z2PWyzEiYrgMPN5jQ/ilR3hExCisdfc8eqojSSMiYoSqHuFJGhERUYMRfXk8FRERdeXxVERE1GLE45463mGMqe5KGhtvxNo5zxvvKCJiMrjiW6M+RdW5L4+nIiKipjSER0RELbboc2oaERFR09rUNCIioo6qn0ZqGhERUYMRT7i3fo321reNiGizvvTTiIiIOtIjPCIihmVt3p6KiIg60hAeERG1GaVNo50knQ4cDNxne++y7evAHqXIZsAfbO/b4thfAquBPmCN7TmdjDUiYrhseu7tqU7Xq84EDmjcYPuvbe9bEsX5wAWDHP+yUjYJIyImILG25lLrbNIBku6UtEzScS32by7pQkk3S7pW0t4N+94n6VZJSyW9v41f8ik6mjRsLwIebLVPkoA3Aed1MoaIiE4x0OcptZahSJoKnAwcCMwG5kma3VTso8AS2/sAbwdOKsfuDfxvYC7wHOBgSbPa9T0bjWcLzouAe23/fID9Bi6VdL2k+QOdRNJ8SYslLX78iT92JNCIiIH0MaXWUsNcYJnt5bYfB74GHNpUZjZwOYDtO4CZkrYF9gKutv2I7TXAlcDr2/UdG41n0pjH4LWM/W0/jyrrHiPpxa0K2V5ge47tOdPWn96JOCMiWjJirestNewI/Lrh84qyrdFNwGEAkuYCuwI7AbcCL5a0paSNgIOAnUf59VoalxYcSetRffHnD1TG9sry8z5JF1Jl4UVjE2FExNDMsBrCt5K0uOHzAtsLGj63yixu+nwicJKkJcAtwI1ULwrdLunTwGXAw1TJZU3dwIZjvJr9XwncYXtFq52SpgNTbK8u668GThjLACMihqbhzKfxwBAv9azgqbWDnYCVjQVsrwKOhCfbhe8qC7ZPA04r+/6lnK/tOvp4StJ5wFXAHpJWSDqq7DqCpkdTknaQtLB83Bb4saSbgGuB79r+fidjjYgYLlP1CK+z1HAdMEvSbpKmUf2evKixgKTNyj6AdwOLSiJB0jbl5y5UT3I68pJRR2satucNsP2dLbatpHoOh+3lVG8ADMuajcQD+2ww3MMiohdd0Z7TtGvmPttrJB0LXAJMBU63vVTS0WX/qVQN3mdL6gNuA45qOMX5krYEngCOsf37tgTWpLd6pUREtJGtto49ZXshsLBp26kN61cBLV+ltf2itgUyiCSNiIhRyHSvERFRSzUJ09TxDmNMJWlERIxQ1RCeAQsjIqKmDI0eERG19PcI7yVJGhERo7A2NY2IiKjDJpMwNZL0+RrnWGX7Y22KJyJi0jBizdq8PdXoUOAfhyhzHJCkERE9qV09wieLoZLGv9s+a7ACkjZvYzwREZNGXrltYvs/hjpBnTIREd2pvcOITAa1GsIlbU01leDMxmNsv6szYUVETA515//uFnXfnvpv4EfAD4C+zoUTETF52PBEGsJb2sj2RzoaSUTEJNOLnfvqPoy7WNJBHY0kImISWotqLd1i0KQhabWkVcD7qBLHo5JWNWwflKTTJd0n6daGbZ+Q9BtJS8rSMhlJOkDSnZKWSTpuuF8sIqLT+t+eqrN0i6HentpklOc/E/gicHbT9n+3/a8DHSRpKnAy8CqqeW6vk3SR7dtGGU9ERFv12ttTtb6tpMvrbGtmexHw4Ajimgsss73c9uPA16g6GkZETBw1axk9U9OQtAEwHdiqdOLr/+YzgB1Gcd1jJb0dWAx8qMVctjsCv274vALYb4AY5wPzAaZuvjmr/mztKMKKiKjPwJrUNJ7ib6h+se8JXN+w/DfV46OR+BLwTGBf4B7g31qUaZWW3epkthfYnmN7ztSNp48wpIiI4UubRhPbJ0n6IvBR259sxwVt39u/Luk/gYtbFFsB7NzweSdgZTuuHxHRTt2UEOoYsl5luw9o2+u2krZv+Ph64NYWxa4DZknaTdI04AjgonbFEBHRDv39NHqpplH3Ydylkt4gaVjfXNJ5wFXAHpJWSDoK+IykWyTdDLwM+EApu4OkhQC21wDHApcAtwPfsL10ONeOiBgLvdZPo26P8A9SNYj3SXqUqs3BtmcMdpDteS02nzZA2ZU01GhsLwQW1owvImLsufceT9VKGm3orxER0XUMrFnbW29P1Z7uVdIhwIvLxytst2rAjojoGb049lTdodFPBP4cOKdsep+kv7Sd4T0ioqe5x5JG3XrVQcCrbJ9u+3TgANr4RlVExGTVzobwocbck7S5pAsl3SzpWkl7N+z7gKSlkm6VdF7pnN12w3kYt1nD+qbtDiQiYrKx29e5r2HMvQOB2cA8SbObin0UWGJ7H+DtwEnl2B2B9wJzbO8NTKXqqtB2dds0PgXcKOmHVG9OvRj4+04EFBExeYi+9jWEPznmHoCk/jH3GgdqnU31+xjbd0iaKWnbsm89YENJTwAb0aEO0bW+re3zgBcAF5Tlhba/1omAIiImE1u1lhpajbm3Y1OZm4DDACTNBXYFdrL9G+BfgV9RDc/0kO1LR/nVWhpOipwCPAD8HniWpBcPUT4ioqsNc+yprSQtbljmN52uzph7JwKbS1oCvAe4EVhTBpQ9FNiNajDZ6ZLe2s7v2q/u21OfBv4aWAr0DyNrYFEngoqImBRctWvU9IDtOYPsH3LMPdurgCMByggdd5XlNcBdtu8v+y4A/gL4au3oaqrbpvE6YA/bj7U7gIiIyayNQ4Q8OeYe8Buqhuw3NxaQtBnwSJln6N3AIturJP0KeIGkjYBHgVdQjVDednWTxnJgfSBJIyKiMO3rp2F7jaT+MfemAqfbXirp6LL/VGAv4GxJfVQN5EeVfddI+hZwA7CG6rHVgrYE1qRu0ngEWFJm63sycdh+byeCioiYHETf2vZ17ms15l5JFv3rVwGzBjj248DH2xbMAOomjYvI0OQREU/Taz3C6w5YeNZg+yWdb/sN7QkpImJysJM0Rmr3Np1nVKZM62PjXVeNdxgR0UN6bcDCdnVlbPnSmaTTJd0n6daGbZ+VdEcZO+XC8jZAq2N/WSZrWiKpI28BRESMll1v6RadHgj+TKrBDRtdBuxdxk75GYMPR/Iy2/sO8W5zRMS4MGLt2im1lm7Rrm/Ssn5mexHwYNO2S8t0rgBXU3VgiYiYlFxz6RbDThplaN59mjZ/ZITXfxfwvQH2mWpu8utbdLdvjGd+f7f8vlWPjDCMiIgRcFvHnpoU6g4jcgVwSCm/BLhf0pW2PwhV7WG4F5Z0PFUnlHMGKLK/7ZWStgEuk3RHqbk8he0FlE4sG/7ZDt2U0CNiMuix3zp1axqbljFPDgPOsP184JUjvaikdwAHA2+xWzcR2V5Zft4HXEg1bHBExITSazWNukljPUnbA28CRjU3uKQDqB5nHWK75fMkSdMlbdK/DrwauLVV2YiI8ZS3p1o7gWo8lF/Yvk7S7sDPhzpI0nnAVcAeklZIOgr4IrAJ1SOnJZJOLWV3kNTffX5b4MeSbgKuBb5r+/vD+mYRER1mg9dOqbV0i7o9wr8JfLPh83JgyB7gtue12HzaAGVXUuYdL+d/Tp3YIiLGUzfVIuqolf4kPUvS5f2d9CTtI+ljnQ0tImIS6LF3buvWmf6TqhPeEwC2b6ZDk5ZHREwe9RrBu6khvO7YUxvZvraaKOpJawYqHBHRM7qoFlFH3aTxgKRnUm6PpMOpJi+PiOhdGeV2QMdQdaDbU9JvqOak7cik5RERk0qSxtOVt5leWfpMTLG9urNhRURMEj32eKru21PbSjoN+Jbt1ZJmlz4XERG9LW9PtXQmVee+HcrnnwHv70RAERGThqkeT9VZJhBJb6yzrZW6SWMr298A1gKUoc37akcYEdGlJukwIq3mMRpsbqMn1W0I/6OkLVn39tQLgIdqHhsR0b3WTqxaxGAkHUg18saOkj7fsGsGNbtR1E0aHwQuAp4p6SfA1sDhw4h1TGy0/hPsu+1vxjuMiJgE2jUCqiZeLWIwK4HFVFNdXN+wfTXwgTonGDJpSJoCbAC8BNiDapa+O20/MdxoIyK6yiRr5LZ9E3CTpHNH+jt8yKRhe62kf7P9QmDpSC4SEdGdJl4jd00zJX0KmE1VKQDA9u5DHVi3IfxSSW9Q0zgiERE9b3K+cnsG8CWqdoyXAWcD/1XnwLpJ44NUQ6M/JmmVpNWSVo0k0oiIrjI5k8aGti8HZPtu258AXl7nwFpJw/YmtqfYnmZ7Rvk8YxQBR0RMfqZ6e6rOUoOkAyTdKWmZpONa7N9c0oWSbpZ0raS9y/Y9yqR2/csqSYP1pftTaa/+uaRjJb0e2KZOjHV7hD+vxfJMSYO2iUg6XdJ9/fNwlG1bSLpM0s/Lz80HOHbQmxcRMRHI9ZYhzyNNBU4GDqRqa5gnaXZTsY8CS2zvA7wdOAnA9p2297W9L/B84BHgwkEu935gI+C9pfzbgHfU+b51H0+dAlxNNa/Gf5b1rwE/k/TqQY47EzigadtxwOW2ZwGXl89PUfPmRUSMv/Y9npoLLLO93PbjVL9jD20qM5vq9ya276Bq0N62qcwrqKbmvnvAkO3rbD9se4XtI20fZvvqOkHW7afxS+Ao20sByi/w/wN8ErgAuHSAwBZJmtm0+VDgpWX9LOAK4CNNZZ68eeV6/TfvtprxRkRMNjsCv274vALYr6nMTcBhwI8lzQV2BXYC7m0ocwRwXqsLSPoOg6Qw24cMFWTdpLFnf8IoJ75N0nNtLx/BC1Xb2r6nnOceSa2eo9W5eQBImg/MB5i+3fThxhIRMSrD6Ny3laTFDZ8X2F7QeKoWxzSf/UTgJElLgFuAG2noyS1pGlXHvYGGBPnX2tEOoG7SuFPSl6iqSwB/TfVo6hmUKWDbrM7NqzZWN30BwJZ7bT3x3lGIiO5Wv5/GA7bnDLJ/BbBzw+edqHpwr7uUvQo4EqB0gbirLP0OBG6w3VjzaDz+yuZtpV155zKN95Dqtmm8E1hG1XjyAWB52fYE1Tu+w3GvpO1LsNsD97UoM+TNi4gYd6YaxrXOMrTrgFmSdis1hiOohm96kqTNyj6AdwOLSiLpN48BHk01necKSTMkbUH1yOsMSZ+rE2TdSZgelXQKcLHtO5t2P1znHA0uomqlP7H8/O8WZZ68ecBvqG7em4d5nYiIjmvX2FO210g6lmoaiqnA6baXSjq67D8V2As4W1IfVRvvk/MaSdoIeBXwNzUut6ntVZLeDZxh++OSatU0aiUNSYcAnwWmAbtJ2hc4YahGE0nnUTV6byVpBfBxqmTxjTKJ06+AN5ayOwBfsX3QQDevTqwREWOqjQ/FbS8EFjZtO7Vh/Spg1gDHPgJsWfNS65UnPW8Cjh9OjHXbND5O9UbTFSW4JS3einoa2/MG2PWKFmVXUg3Z2//5aTcvImLCmZwtqSdQ/VH+E9vXSdod+HmdA+smjTW2H8rQUxER69TtuDfR2P4m1dBQ/Z+XA2+oc2zdhvBbJb0ZmCpplqQvAD8ddqQREd2mjcOIjBVJO5XhSO6TdK+k8yXtVOfYuknjPcCzgceoWuZXkTnCIyLaNozIGDuD6qWkHaj6xX2nbBtS3benHqFqLBlWg0lERNebeAmhjq1tNyaJM4cY4PBJQw04OOou5xERXWti1iLqeEDSW1nXp2Me8Ls6Bw5V0+jvcn4YsB3w1YYL/HJ4MXbe9KmPsd+my8c7jIiYBL46dJF6JmfSeBfwReDfqb7BTyk9zYcyaNLo73Iu6ZO2X9yw6zuSFo0s1oiILjI5k8YngXfY/j1UU1ZQVRLeNdSBdRvCty7v8VIusBuw9QgCjYjoKpO0IXyf/oQBYPtB4Ll1DqzbT+MDwBWS+p/9zKSMLBsR0dMmXkKoY4qkzZtqGrXyQd23p74vaRawZ9l0h+3HRhRqRES3mJi1iDr+DfippG9Rpb03Af9c58Ch3p56nu0bAEqSuGmwMhERPWcSJg3bZ5e5PV5ONRXFYbZrTXI3VE3jDEkvpfX8Fv1Oo+azsIiIrjMJkwZUk+kxgtlQh0oamwLXM3jSuH+4F42I6AZi0j6eGrGhXrmdOUZxRERMPgbVm2Cpa9R9eyoiIlpJTSMiImrrsaRRt3NfW0naQ9KShmVV82BZkl4q6aGGMv84HrFGRAxmknbuG7G6070KeAuwu+0TJO0CbGf72pFctMwzvm8591SqecAvbFH0R7YPHsk1IiLGRBclhDrq1jROAV5INVAhwGrg5DbF8ArgF7bvbtP5IiLGRmkIr7N0i7pJYz/bxwB/Aihdz6e1KYYjWDc8b7MXSrpJ0vckPbtVAUnzJS2WtPjh3z/eppAiImpyzaVL1E0aT5THSAaQtDUw6twpaRpwCA1z1Ta4AdjV9nOALwDfbnUO2wtsz7E9Z+PN25XHIiLq6bU2jbpJ4/NUbQ7bSPpn4MfAv7Th+gcCN9i+t3mH7VW2Hy7rC4H1JW3VhmtGRLRPj9U06g5YeI6k66naHwS8zvbtbbj+PAZ4NCVpO+Be25Y0lyrB1ZpZKiJiTHRZQqhjqAELt2j4eB8Nv+AlbVHGYB8RSRsBrwL+pmHb0QC2TwUOB/5W0hrgUeAI2z32nyciJjIx+BhL3Wiomsb1VHlUwC7A78v6ZsCvgN1GemHbjwBbNm07tWH9i1TTEda20ZTHec4GvxppSBERw9ZNb0bVMWibhu3dbO8OXAK81vZWtrcEDgYuGIsAIyImtB5r06jbEP7npTEaANvfA17SmZAiIiaRNiYNSQdIulPSMknHtdi/uaQLJd0s6VpJezfs20zStyTdIel2SS8c9XdroW7SeEDSxyTNlLSrpONJo3RE9Lqar9vWeeW2dGs4meqt0tnAPEmzm4p9FFhiex/g7cBJDftOAr5ve0/gOUA7XlZ6mrpJYx6wNdVrt98GtmFd7/CIiN7VvprGXGCZ7eW2Hwe+BhzaVGY2cDmA7TuAmZK2lTQDeDHVpHjYftz2H0b3xVqr+8rtg8D7OhFARMRk1saG8B2BXzd8XgHs11TmJuAw4MelK8KuwE5AH9WEeGdIeg7VS0zvs/3HtkVX1B2w8Ie0yJW2X97ugCIiJpNh9PbeqszL3W+B7QWNp2pxTPPZTwROkrQEuAW4EVgDrA88D3iP7WsknQQcB/xD7ehqqjufxocb1jcA3kAVaERE7xrem1EP2J4zyP4VwM4Nn3cCVj7lcvYq4Eh4cvTxu8qyEbDC9jWl6Leokkbb1X08dX3Tpp9IurID8URETC7te532OmCWpN2opos4AnhzYwFJmwGPlDaPdwOLSiJZJenXkvYoU0+8AritbZE1qPt4qrFn+BTg+cB2nQgoImKyEO0bjND2GknHUvWLmwqcbntp00gZewFnS+qjSgpHNZziPcA5ZSDY5ZQaSbvVfTzV2DN8DVV16KhBj4iI6AVt7LhX+sMtbNrWOFLGVcCsAY5dAgz2+Kst6iaNvWz/qXGDpGd0IJ6IiMnDoLVd1N27hrr9NH7aYttV7QwkImIy6rX5NIYa5XY7qneHN5T0XNa9EjaDqrU+IqK3dVFCqGOox1OvAd5J9erX5xq2r6bqzh4R0dO6qRZRx6BJw/ZZwFmS3mD7/DGKKSJi8kjSWEfSW21/lWp8kw8277f9uRaH1SLpl1Q1lj5gTXOnl9Jx5STgIOAR4J22bxjp9SIi2q7L2ivqGOrx1PTyc+MW+9pxq15m+4EB9h1I9WrZLKrxV77E08dhiYgYN6L3JmEa6vHUl8vqD2z/pHGfpP07FlXlUODsMsXr1WWs+O1t39Ph60ZE1Ndjs1DXfeX2CzW3DYeBSyVdL2l+i/2tRnzccZTXjIhoq7xy26DM/PQXwNZNbRozqLq5j8b+tldK2ga4TNIdthc1Xr7FMU+79SXhzAfYccepPHP9R0YZVkRETV02lWsdQ9U0plG1Z6wHbNKwrAIOH82Fba8sP++jmtxpblORIUd8LMcvsD3H9pwtt6xbcYqIaA+trbd0i6HaNK4ErpR0pu2723VRSdOBKbZXl/VXAyc0FbsIOFbS16gawB9Ke0ZETDTdlBDqqDv21COSPgs8m2o+DWBUkzBtC1xYvVXLesC5tr/fNJrjQqrXbZdRvXLbkREbIyJGzPRcQ3jdpHEO8HXgYOBo4B1UUwuOiO3lVBOfN29vHM3RwDEjvUZExFjopkbuOuo2Amxp+zTgCdtX2n4X8IIOxhURMTm45tIl6tY0nig/75H0V1QN0jt1JqSIiMmhnZMwTRZ1k8Y/SdoU+BBV/4wZwPs7FlVExGRgp02jFdsXl9WHgJcBSErSiIie12tvT42mY8PTBjCMiOg16YT/U5AAAA7VSURBVBFeX6se2xERvcNAj033Opqk0Vt3KiKilR77TTjU2FOraX1LBGzYkYgiIiaRbnr0VMdQw4hsMlaBRERMSnl7KiIianHvvT2VpBERMUJV577UNCIioq7UNCIioq5eq2lk1qKIiJGqO1hhzbwi6QBJd0paJum4Fvs3l3ShpJslXStp74Z9v5R0i6QlkhaP9qsNJDWNiIgRM2pT5z5JU4GTgVdRzVx6naSLbN/WUOyjwBLbr5e0Zyn/iob9L7P9QFsCGkBXJY3119+HHXYcXYJ97J7d2xRNRPSE9j2emgssK/MNUWYtPRRoTBqzgU9Vl/UdkmZK2tb2ve0KYijj8nhK0s6SfijpdklLJb2vRZmXSnqoVLWWSPrH8Yg1ImJAbusc4TsCv274vKJsa3QTcBiApLnArqybpsLApZKulzR/NF9rMONV01gDfMj2DZI2Aa6XdFlTNQzgR7YPHof4IiLqqV/T2KqprWGB7QUNn1uN59d88hOBkyQtAW4BbqT6fQqwv+2VkrYBLpN0h+1FdYOra1yShu17gHvK+mpJt1Nl1OakERExsdV/OvWA7TmD7F8B7NzweSeqCe/WXcpeBRwJIEnAXWXB9sry8z5JF1I97mp70hj3t6ckzQSeC1zTYvcLJd0k6XuSnj3A8fMlLZa0+P77RzxteUTEiMiutdRwHTBL0m6SpgFHABc95VrSZmUfwLuBRbZXSZpentogaTrwauDWtn3JBuPaEC5pY+B84P0lgza6AdjV9sOSDgK+DcxqPkep3i0AmDNnTm+9MB0R48tAX3t+7dheI+lY4BJgKnC67aWSji77TwX2As6W1Ef1ZOaocvi2wIVV5YP1gHNtf78tgTUZt6QhaX2qhHGO7Qua9zcmEdsLJZ0iaatOv04WEVGXqF2LqMX2QmBh07ZTG9avovUfz8uB57QtkEGM19tTAk4Dbrf9uQHKbFfK9b8lMAX43dhFGRFRQ/884UMtXWK8ahr7A28DbilvAUDVaWUXeDKzHg78raQ1wKPAEXYX3fmI6A499mtpvN6e+jFDTBdr+4vAF8cmooiIETAZsDAiIurrtQELkzQiIkbMsLa3qhpJGhERI2XSphEREcPQWxWNJI2IiNFIm0ZERNSXpBEREbXY0Ndbz6eSNCIiRiM1jYiIqC1JIyIiajHQpjnCJ4skjSbrK7ckIuoyOG0aERFRVx5PRURELSZvT0VExDCkphEREfV01wRLdSRpRESMlOm5UW7HZbpXAEkHSLpT0jJJx7XYL0mfL/tvlvS88YgzImJQPTbd63jNET4VOBk4EJgNzJM0u6nYgVQTqM8C5gNfGtMgIyLqSNIYE3OBZbaX234c+BpwaFOZQ4GzXbka2EzS9mMdaETEgGzc11dr6RbjlTR2BH7d8HlF2TbcMkiaL2mxpMX3339/2wONiBjUWtdbusR4JQ212NZ8V+uUwfYC23Nsz9l6663bElxERG099nhqvN6eWgHs3PB5J2DlCMpERIwf994c4eNV07gOmCVpN0nTgCOAi5rKXAS8vbxF9QLgIdv3jHWgERGDSk2j82yvkXQscAkwFTjd9lJJR5f9pwILgYOAZcAjwJHjEWtExMDcVY3cdYxb5z7bC6kSQ+O2UxvWDRwz1nFFRNTWg0Ojj1vnvoiIruC19ZYaanR63lzShaXD87WS9m7aP1XSjZIubtO3e5okjYiIETLgta61DKVmp+ePAkts7wO8HTipaf/7gNtH+70Gk6QRETFSdjtrGnU6Pc8GLq8u7TuAmZK2BZC0E/BXwFfa9fVaSdKIiBiFYdQ0turviFyW+U2nqtOh+SbgMABJc4FdqbojAPwH8H+Bjr4D3FWj3F5//fUPSLp7iGJbAQ+MRTyjlDjbZzLECImz3YaKc9fRXmA1v7/kB2u/sVXN4g/YPmCQ/XU6NJ8InCRpCXALcCOwRtLBwH22r5f00prxjEhXJQ3bQ3YJl7TY9pyxiGc0Emf7TIYYIXG221jEOUQSGK4hOzTbXkXpfiBJwF1lOQI4RNJBwAbADElftf3WNsYH5PFURMREMWSnZ0mblX0A7wYW2V5l++9t72R7ZjnufzqRMKDLahoREZNVzU7PewFnS+oDbgOOGus4ezFpLBjvAGpKnO0zGWKExNlukyXOJ9Xo9HwV1RxDg53jCuCKDoQHgNxFY6JERERnpU0jIiJqS9KIiIjauippSNqgjMdyk6Slkv5f2f7G8nmtpDlNx/x9GeflTkmvmYhxSpop6VFJS8py6sBnH5M4PyvpjjL+zYWSNms4ZiLdz5Zxjsf9HCTGT5b4lki6VNIODcdMpHvZMs6J9m+zYf+HJVnSVg3bxvx+diXbXbNQdY7ZuKyvD1wDvIDqjYM9qBqH5jSUn03Vw/IZwG7AL4CpEzDOmcCtE+h+vhpYr2z/NPDpCXo/B4pzzO/nIDHOaCjzXuDUCXovB4pzQv3bLJ93pnoD6W5gq/G8n924dFVNw5WHy8f1y2Lbt9u+s8UhhwJfs/2Y7buo5u6YOwHjHBeDxHmp7TVl+9WsG8Zgot3PgeIcc4PEuKqh2HTW9QCeaPdyoDjHxUBxls//TjWcRmOM43I/u1FXJQ14cmjgJcB9wGW2rxmkeJ2xXjpimHEC7KZqyOMrJb1oDEIEasX5LuB7ZX0i38/GOGEc7udAMUr6Z0m/Bt4C/GMpPuHu5QBxwgT6tynpEOA3tm9qKj5u97PbdF3SsN1ne1+qvyrnqmm8+SZ1xnrpiGHGeQ+wi+3nAh8EzpU0Y7zjlHQ8sAY4p39Tq1N0Psphxzku93OgGG0fb3vnEt+x/WG3OkWnYxxBnBPp3+Y+wPE8NaH1G7f72W26Lmn0s/0HqraBwcaGGXKsl06rE2epUv+urF9P9Tz2WWMS4LoYnhKnpHcABwNvsd3/P9+Eu5+t4hzv+znIf/NzgTeU9Ql3Lxs8Ged438ty3f44D6Vqr7hJ0i+p7tkNkrZjAtzPbtFVSUPS1g1vyGwIvBK4Y5BDLgKOkPQMSbtR9bS8dqLFWcpPLeu7lziXj1eckg4APgIcYvuRhkMm1P0cKM7xuJ+DxNjYu/cQ1v07mGj3smWcE+zf5o22t7E909UYTCuA59n+LeN0P7tRtw0jsj1wVvlHPAX4hu2LJb0e+AKwNfBdSUtsv8bVuC7foBrDZQ1wjO2xmCV+WHECLwZOkLQG6AOOtv3gOMa5jOotlMskAVxt++gJeD9bxsn43M+BYjxf0h5UcyDcDfSPMzTR7mXLOJlg/zYHKjyO97PrZBiRiIioraseT0VERGclaURERG1JGhERUVuSRkRE1JakERERtSVpREREbUkaPUTSw0OXGtX5vyJpdln/6AiOnynp1mGWf7SMP9Rq/yckfXi4cUxGkjaT9HcNn5+paqjyjv43j96TpBFtY/vdtm8rH4edNEboF2X8oY7p7/E8wW0GPJk0bHf8vkRvStLocZJ2lXS5qgl2Lpe0S9l+pqTPS/qppOWSDi/bp0g6RdXENxdLWtiw7wpJcySdCGxY/tI9p7kGoWqCnE+U9eermkjnKuCYhjJTVU2idF2J7W9qfp/jVU2y8wOquUn6tz9T0vclXS/pR5L2bNh+dbnOCf1/mUt6qaQfSjoXuKVse6uqiX+WSPpyw/AZr5Z0laQbJH1T0sZl+4mSbivx/+sgMW9delxfV5b9y/a55f7fWH7uUbY/uyGOm1UN8XEi0F+7+GydexUxIp4Ak3pkGZsFeLjFtu8A7yjr7wK+XdbPBL5J9YfFbGBZ2X44sLBs3w74PXB42XcFZfKoxmvRNFEP8GHgE2X9ZuAlZf2z/eWA+cDHyvozgMXAbk2xN5/3+VS/4DcCZlDNmfDhsu9yYFZZ3w/4n7J+MTCvrB/dHzfwUuCP/dekmiDrO8D65fMpwNuBrYBFwPSy/SNUo6xuAdzJulEXNhvkv8u5wF+W9V2A28v6DNZNIvVK4Pyy/gWqARgBpgEbNt+Lwf6bZ8kymqXbxp6K4XshcFhZ/y/gMw37vm17LXCbpG3Ltr8Evlm2/1bSD0d6YUmbUv0yvbLh+geW9VcD+/TXYoBNqQaZu2uQU74IuNBlcEJJF5WfGwN/AXyzjEEFVSKC6vu/rqyfCzTWCK51NWEPwCuoktJ15RwbUs3j8AKqpPqTsn0acBWwCvgT8BVJ36VKTgN5JTC7IbYZkjYp3/msUpMw1URDlPMfL2kn4ALbP284NqKjkjSiWeNgZI81rKvp53Cs4amPQjdoONdAg58JeI/tS4Z5rVbnmwL8wcN/xv/HpnjOsv33jQUkvZZqAqB5zQdLmkuVbI6gmn/i5QNcZwrwQtuPNh3/BeCHtl8vaSZVTQ7b50q6Bvgr4BJJ72YMRpaNgLRpBPyU6pcaVDOy/XiI8j8G3lDaNraleozTyhOS+v8yvhfYRtKWkp5BNb8FruZBeEjSXzZcv98lwN/2n0PSsyRNHyK2RcDrJW1Y/lJ/bbnOKuAuSW8s55Kk55RjrmbdHBZHNJ+wweXA4ZK2KefYQtKu5fj9Jf1Z2b5RiXVjYFPbC4H3A4MlrEtZN6kRkvrLbgr8pqy/s2H/7sBy25+nGvJ7H2A1sMkg14hoiySN3rKRpBUNyweB9wJHSroZeBvwviHOcT7VPAW3Al8GrgEealFuAXCzpHNsPwGcUMpezFPnDjkSOLk0hDf+pf0VqmGsbyiN6F9miJqx7RuArwNLSpw/atj9FuAoSTcBS6km7IHqF/oHJV1LNdx2q++Cq7fCPgZcWu7VZcD2tu+n+oV+Xtl+NbAn1S/wi8u2K4EPDBL6e4E5pVH7NtYNO/4Z4FOSfgI0vsH118Ctql413hM429VESD+RdGsawqOTMjR6DJukjW0/LGlLqols9nc10c1YxzETuNj2YFPlDnWOjYBHbVvSEVSN4ocOddxkIelh2xuPdxzRPdKmESNxsapZ06YBnxyPhFH0AZuqmqxqpH0Sng98UVVL8h+o3iCb9CQ9k6q2de94xxLdJTWNiDEi6XjgjU2bv2n7n8cjnoiRSNKIiIja0hAeERG1JWlERERtSRoREVFbkkZERNT2/wEqac4li0dYhwAAAABJRU5ErkJggg==\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": 7,
"metadata": {},
"outputs": [],
"source": [
"area_weighted = (ds * coslat).sum(['lat', 'lon']) / coslat.sum(['lat', 'lon'])"
]
},
{
"cell_type": "code",
"execution_count": 8,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"[<matplotlib.lines.Line2D at 0x2b339cce96a0>]"
]
},
"execution_count": 8,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXAAAAESCAYAAADkJY5uAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO2dd5gUVdaHf6d7EhPIQw5DDipBBwQBCQIKuqurrmmVXXfNEeMqa1oj664on7rr4mLO2VVQJAnCCjggSUBAMgwwpGFgmNj3+6OqerqrK3ZXdXf1nPd55pnuCrdO37p16txzzz2XhBBgGIZhvIcv0QIwDMMw0cEKnGEYxqOwAmcYhvEorMAZhmE8CitwhmEYj8IKnGEYxqOkxfNizZs3FwUFBfG8JMMwjOdZvnz5ASFEvnp7XBV4QUEBioqK4nlJhmEYz0NE27W2swuFYRjGo7ACZxiG8SiswBmGYTwKK3CGYRiPwgqcYRjGo7ACZxiG8SiswBmGYTwKK3CGYeols37aixveXJ5oMWKCFTjDMPWSFTsOY+6GfYkWIyZYgTMMUy8RQvqzw9tLt2P0lAXuCBQFcZ1KzzAMkywEAgIBmxr8L5+udUma6GALnGGYeklAAF5fEZgVOMMw9ZKAELZdKArJshg8K3CGYeolihKORhkHkkN/swJnGKZ+oijhaJRxbZJocFbgDMPUSwIxWeDG56zdXRoXNwsrcIZh6iWKEa1Ws7e/9yO+XrvX5Fx95fzl6j047/lF+GJ1cYwSmsMKnGGYeoliIa/Yfhhz1tVN6Pl85R7c8JbxDE0jF8qmfccAAL/sP+aAlMZwHDjDMPUSRQlfOm0JAGDb5HMtnxsI6O+Lp3ecLXCGYeolsYxD1srW+4KNJVix43D4TnkfUfTlW4UVOMMwCaW6NoCXF25BVY2BWesCVgcZtQYkFR/4719Zhgv/+b/g9oUbS1BeVQsAILivwVmBMwyTUN5ash1PzFyP6Yu2xvW6WgORakX99dpinPf8Iny+ck/4uRrm+8FjlZjwyjLMXCMNXrIFzjBMynOsogYAcLyyJq7XVevgb37ai52HToRtUwYkN+0vC9teq6H8K+UeRIX8Pw7623wQk4jaA3gDQCsAAQDThBBTiagfgJcAZAGoAXCTEGKZm8IyDJN6KKowHhZrKGoL/DqN3OBB2VTqWCsKRSlP2ReP32MlCqUGwF1CiBVElAdgORHNBvA0gL8KIb4iovHy9xHuicowTCqiKL44629LeVCUY77bVIKXFvwS3K4VhaJsC30xlJZXo1F2eixiGmLqQhFCFAshVsifywCsB9AW0supoXxYIwB7tEtgGIaxQJxNcCupZIVsg6/aVYqaEKtbaAQLKuUp/vFpC7eg76PfYM+RExHHOoWtOHAiKgDQH8BSABMBzCKif0B6EZzhtHAMw6Q+enr0/+ZuwpaSY3jusv6Wyyo9UY2qmgDy8zJNj7WkwHUOUXtQHvxsLbrk54TtOyr79ncfOYE0P6FFXpbp9exieRCTiHIBfAxgohDiKIAbAdwhhGgP4A4A03XOu46IioioqKSkxAmZGYZJIZTID5/KAJ8yeyM+W2mvY3/GU3Mx4Ik5lo61EgeuF2qo3v5+0U4s3HRALjd83+9fWYaBT8y1JJNdLClwIkqHpLzfFkJ8osgFQPn8IYCBWucKIaYJIQqFEIX5+fmxysswTIoSTdz07iMn8OgX64IDh8flGGwjjlfWYOaaYktx4HpKXr09EBBBF4tagStx4W4kt7IShUKQrOv1QogpIbv2ABgO4FsAowBsclw6hmFSnljU2p3vr8TSrYcw/pRWKCxoaumcSZ+uwecr96B1I3OXhpavGwB2Hi4PU9Q1AYFaeRRTT+nXBATS/c76+a34wIcAuArAGiJaKW+bBOBaAFOJKA1ABYDrHJUsBSguPYHDx6vRu01D84MZpp6i6MFoxjBrg1av9XN2HioHUBd/boSe0Xz1qz9EbKup1bbAFaprA0j3Ozv1xlSBCyEWQT/C5zRHpUkxBj81D4C9JDlMcvLTnlLM37Aft4zqlmhRUg7Fyo3GNlWUvto9ccGLi9GxWTamGgyAWtH5dnoHdfnFtfdX1QSQnWGjQAvwTEyGCeF4ZQ0K7puB/64KHzz79QuL8Y9vNiZIqtQmFgtcz2++cueRiOnvwevJ/2OJQtGiutb44Kpa53O9sAJnmBB2Hpa618/PDR/SUbrqRgNRRdsOYfBTc1FWUe2egClI3UzMqDR4WBmWrifC/xsea6Nks2XW3EjWxQqcYUJQ/JhpOr5Ko4f0H9/8jOLSCqzeVeqKbKlKLMEZisq3Yk0v334Y5zy3ECVlldJ1rShnG7LVmChwoxzi0cILOjBMCMpDqBctYPSMKgNU1S50lVOZoA/c5YmYn6zYhQ1765JSmQ18rtp5JJigygpaGQpDsWPNW4UVOMOEoISC+XS0iZGlpyjwGhNfKBOOUl/pPvsOgeBtslDlFdXhytgsLvv8FxfbkqXGxMR2Y41jVuAuUFUTML2ZTHKiKBM9a9BIgStKXyvVKKNPRbU00SUrPQoFLjtRrNS4+p46fZvMfOButApW4C5w0b/+hzW72Q/qRZRnUK87bPSMKi/tNPWccMYQxTLOTPMDkBRhZY35jEqgTinr+Z+nztmEIV2bobCgaUS8itMK1dQH7sKLnQcxXYCVt3cJ5nTWediMHkLFetdzvzDaqJX13R+uQu+HZmkeGwgIPP7luuBknLqIEu378uycjbj4pe8BRN4XpxWqqQXuggnOCpxJCZZsOYhPf9wVcznBmX06HjCjgapgqGFc1yVPLmoDArvkUEyrKIpNUaif/rhb99j1e4/iP4u24uZ3VoSdY6XG1S52pxWqmQXuhhOFFTiTElw2bQnueH9VzOWYWWVGz6iiuJ0Y/vjmp714dXF814h0gufnbcLQv83HjoPWlXiw3izoN6Vu6+Lype/VNQEcKa8yOdvdnlEiLHD2gTOe5r1lO/B/c83zqB06XoUm2ekor6pFTqZ+sw+uDqPzrB86XoWqmgBaaSRCUp5fJ55TZXmvq4d0cqC0+LB2dykWbJRSRu8+cgIdmmXbOt9Kz0U5RnGHKPfr+reWmypItz1bNSbho3bytViFFTjjae77ZI3uvuraAK57owiDuzTDkzM3YHj3fCzYWIKF94zUVS6Khaf3sI+esgCATn4bE3+sEUIIVNYEkJXut31uMrD7yAmc9/yi4Hc7sfB1LhTrxyr3xyz/SChujy2bR6GwC4VhLLP9YDnm/1yCJ2duAICgdTh90Rb86bXIbHJAbCGAdlwBat5augM9H/wau3WW33rgszWYs25f1LK5jdp9YSeM1mwgMuxY+X/dDEzLl4kq37gdzHzgPIgZR+Zt2Icr/7PUlSTs8UQIgdnr9plaB6mInsX1+vfbMXfD/oh7++TM9fj2Z0nJW33YS8oqg3HMdcXZr+sv5MRLuw5p+47fWrID17xRZLtcN9h64DhW7jwStk39mJgldtLCyqMWUJngVp/Pgvtm4MPlO23LZAeOQkkirn1jORZtPmBJ8QUCAu8s3WE5djWezF63D9e+URS2ojYA7C+rwOb9xxIkVXwwC+dTP1DTFm7Bu8t2hG3bX1aBxZsP6JYx4Ik5uGr6Uqk8eZtek/l85W7M26BtRSuZ6vRysNhBCBGz4XHwWKVuxM3If3yLC0xmKUaTTsBOdkBf0IVivXz1TEyn4ThwB/luUwlKT0SfFU5pIEcrarDGJDnRjDXFmPTpGkyd486iRDW1AVw+bQmWbDlo+9zDctd224HjYdsHPTk36M9NVcwUuNEDpZx66b+X4Hf/WWpYzg/bDoeVp1fs7e+txB9fC7eiAwGB4tITQYWX4YAC7/3QLIyb+p2lY695/Qec/8KisG0HjlXitMfnYMps6+lz1b85Gl31+co9FuSWB5nlb24oRS+Rkgr88PEqXDV9GW5+e0XUZShd6Ktf+wG/emGRoUVxrFJa2ePQcbMwpugoLq3A91sO4q4P7IfJpcnBr2rroD54VMyiDoz83cqpW1UvPiOCvlwbLpTpi7Zi8FPz8LOcZMlvMtL27OyNun5yhRPVtWFJmyTZBBZuLImwques349VKgPlwDEpW9836/Za+g1A5G9esuUg+jwyC0d1Uuve9u6P+GpNcdi5K3cewfrio8bXCXpQlCgUyyImHHahWKRCdmVs2l9mcqQ+ysO/Svb1GVpr8v9ktAbS5Kx6yZAhr7o2gAmvLMOPOw4nWhQA1h4oOxaxmQtFi1W7pPZVozMJSO0KmTp3E26SDZPagMCmfdba+H9X7cGEV5bhgyLrfmA7g37qunx76Q4crajBhmJt+f67ag9ujMLAUuo2mMMqCZ85PTgKxSaxjDrbmQ5dt6yT9n47jey9ZTvw4Gdrg9+rawOmsclGKBny9hw5gQc/W5tQRb7zUDkWbizBnVH0JKwSamGaT8qp2x9xj+TKzkyz9ogU3Dcj+LK3c79zMsIjedXBG1pFKQsDPDt7I8Y8uxCbZUNlz5ET+HK19io02w5Ig6M7TWZJ3vLOCsxYXWxF9CDzN+zHrsPavQKzHoVdhOpZSEajSQ+OAwewZlcprn2jCJ/efAbKKmrQvWVexDGxLNGkYCdzmVl3LiAAq4tRK3HNj11wMgCg21++Qg/5N6plqqiuRWV1AI2y03XLUx6gFTuOYMWOIxjVswVG9mxhTRiHyZRjnE9UuTfYG3oLzB6Y0P1ag9Ub95UFezBOc6S8Cnd/uDqYvjYoh6qhBYQAVEpd6RQUbT8EANh/tBJdW+Th/BcXBxcrUJORpu1KU/NliPLeebgcD3++FpPO7RVMNKXF1TohmUBdYq/PV+5Gk+wMnNk9P6Ku7ejgujBCQkV1rcdcKGyB4/l5m7D3aAUGPzUPY59dqBn54UQ1qS1wo7r3mYQ0aVkJ1bUBPDt7I45Xmq+M/bPcTVb3KO79aDX6PvqNYcNQ73PaIoqGChejdUJ/rx0LXK04V+08grHPLsThcvsD4Vaswm/W7cOc9fswXw5b1DtXaGyra2/S99W7S/HMNz/rKm8tvli1B7996X+Gx5RX1eL177fj67XWfeF63P7eSkx4ZRmA2Nx5ym9evuMwej74taciqdx413jOAldTUyugNzM6FlUVYYEbVL9y6Pq9ZTj72YUR+7We54+W78LUuZtQWRPAfeN6RiWjsvCukYWvvnYiU50q7o1qF9YGVAj9uWYWT6i7xW6c/I6D5WiSo93zsWJoZWdoW7TqQcaAEBEvblIp8MlfbTC/oIpb3/0x7Pv9n6xBfl6m5rGVMYTfaU3oUS/ua7Xm9xw5gfIqyeDx4rwGzoUCLcUaiRNdlcjUkwbHyv0YvRF0LYtMCe/TezEs2nQgYlkvPd0r/d7wnRXVtej54NcY07tl2HYn4oyjJR7+ytBLmD3jN7y1HL+UHMcPfxltIZNcOGf+fT56top031m5LqD/MGu5F9RtRGkHdgbFgsuW6Zg16vh3K3zww040bKDvvgMiJ/QUl54I+413f7jKcrjvGZPn2ZYxmXDDheI5Ba5GqZTS8mqAgEYN0iNCjcy4+tVl8BFh+h8GBLeplaVR5ZsNlmqdqviBs9O1b8GV0yNjj/V+j/I8lJRVwu8jNM3JCHanZ6umX6t9urUBAb+PMHf9Pgzs1BR5WcYPZCwoD25Uq49bJFSpmVlpS7YcCn42W89QC3WoXlAGG1PC1Xyxeg+eCYm/1lLg/jiH0Om9KO79eLXpuWp3yeCn5qFjSB6aj5abpwD+eW8ZWjeOTB7mNdiFgkhlqVRK30e/ASAlGVKen0PHqzB6ygK8PKEQnZrn6Jap9kMCkUrGqPLN9JG2D1za5sQgmfKADXhiDgCpDvRkUq87WBsQ2Hu0An96vQhn9WwR9hJzmvhb4NauN27qd466lmL5lW8tCbeEBUTEi8jnC8/EZ0mmGAb2Y7ltWj2b7TZSzQLA2c8tRGeD59crcBy4BlqWk9KwT1TXYvP+Y5i+aIvtciMscAM3oJlFqQ5Xe+abn7H9oPUJIsHr6GzXahh6YZBavn0lLO2XEncHhBRjzE0vfGhdWH1g1hcfdXQVJUsWuEXhVu8qxUZVrHfQhRKFQoim7mPSOw4prS02JlQlK+xC0UCrG6m2TLIMQqD0iLTA9SvfzHgLPXPnoRN4ft7mun12bqqNhXatWlpC1A1s2vUD20XPpVFSVolXF2/F3WN7BK3LaAm9T4mKEXbyspdNWxKxze8jbD943LJuLNp2KJhwKxpi+T21ARGcqVzf4ThwDbSUgvrBbaAz4m8Hw0FMs6RJIdZ7eXV4Y3ZAf2uWoVeuum5Cv0fjB7ZDXSa58O33f7IGc9bvw+AuzTCsW35M17AziOkWbl928eaDGP73by0du+tweXBNyERw/VvLPRkx4gZuzMT0vAJXW7CDnpyLk9o0DNumNx26vKrG8mzNWLo/oUqyXDWJxU6p+oOY5i+x4PVE5Hfl2ERZ4EosvxOXDy0iURa4+rql5dU4obKA4yWa2vqNygcu1+rr/9uGF+dvNjk6HFbeIbAFHtkA1e1j79EK7D1aYXiOQu+HZgVnqAHAP7/djF/1aYP2TbMj1LrhWogmN+bAsUocKq/CjoPlKC4Nl80RC9yGTFoWeHAhX5e1ijJZxl0fePx6FPoyhH8/Y/JcHFe9uBP1cokmvYQi6sP//clhaeoXHIWCSGUcywoqQF1eCQB4+uuf8VHRLsy7e0TEcbF0f8ZoTO5xotxgGRoDrHoKQq3Trn9zOdo3kcK63NYpegpVx7MSFeEWuAMFRsF/vtuCx75ch58fHwcAEcr7hjeXo6wy+lTH8YZtaGfgiTwaOG1lHa3QGXBxqRXbuamHjlfhnOcW4uUJhWjftC6WVusloN91Dd/+v18OApDyjLu56OtzczYGJzrpuYKcuP6Dn61Fx2Y5uHNM94RZudvkMLnq2kAwmVgoX/8U+9T0WLDt1vBQwqhkxo326H0F7nClKMmFzFw1ocRiRds58+DxKhw8XoW3lmzH/eN7GcpmlFhLDzcn2DxnYbGLDcVlWL2rFDeP7Br1dT6XlyZ7Yd4m9FaNhcSb45U1aJydkVAZ1Lz3ww68YNOPzerbGdyoR0/FgVdp5NBwupusN5CnVtIF982IyCcRFVG8gNRKQWuA1eogZijxypKi9554YuZ6/H3Wz47EywYEsHa38eIAbpOM4XMHjtlfdIQHIp2hXmcj3LivDN0f+AqzfgqfGu68Ba4MtIVrmeOVNREr7nwhJ5OKRYRong11r9xKLLzCniMndPM9u+lCCaU2IHDT28uDS9VFLmIQHzncpiaKhX2TkVS5H4nkpStPxZkxhshq4RkXyqZ90ixBtTXgtA9cz9oY++xCBIQ0TV1Notv32t2lERneNJLAAQAmvr9Stxw7i1jEQllFDWau2YvCjk1xSrtGEfsDQsAHwgdFO5GbmYbxp7SOi1xOE+sAe6woBkasBITQXRqNsUZmuj/mSWpamCpwImoP4A0ArSCllp8mhJhKRO8D6CEf1hjAESFEP8cllGnYQFvUWDK/aaGnw9zqRQZjsGsDEdEKeqh7B1oJ9aPpmRCA5dsPo1frPGRnRNb3nHX7sO3gcVwzrLPtsrVQJFSLqtT1vR9JyZK0XppeYM2uUsy0ubqNkzji4gPw0oJf8PiM9Y6UVV9xyziyYoHXALhLCLGCiPIALCei2UKIS5UDiOgZAM4lk9BAL37Vaf+ckunNqL6dXJZMuc5fPl2L922sV2hGNAp879EKXPSv/+E3/dvi2Usj38XXvCGtqP7G99vRPDcDn9w0JCYZ7SyAoWbnoXI8O8f6qumJwKi34yWi8Zsz4biVht9UgQshigEUy5/LiGg9gLYA1gEASaELlwAY5Y6IErqDcg47MKy8KT9dsTtcBge6yh+vME+raYdoXmzKKWZJrXYcKseOQ/YyymmhV21WqnPSp2vw3aYDMcvAMPHALQvc1iAmERUA6A8gNFn1MAD7hBDmcWIxoBvVbOFht1N3VvxU1SoHsxOuTqc9NLF0TGoDAle/ugz7VTNanUZ5Kavrz4rv2M2QR4ZxGreaq2UFTkS5AD4GMFEIERqfdTmAdw3Ou46IioioqKQkMu+2VfQs8C9W7cFXa5zzMxIBHxbtjFhJJBR1bvFYImG2HyzHO0t3OB5iFItMP+05ivk/l+DfC62n4RVCYMbqYtTYcC/pSWhFdlbfjJeIJoWBFSxFoRBROiTl/bYQ4pOQ7WkALgRwmt65QohpAKYBQGFhYdRaRS/axI6SscKR8mrc85HxSiN+1es0Fmv3y9XFYSuBW8HK29yJ6JzQPDFm/HfVHtz+3krcP64nrh/exdI5um4xC+8ANsAZL+GWD9z0CZV93NMBrBdCTFHtHg1ggxDCWQeuBrFkyhNCWoTWKdSSuBGgb0R5Va3pUlROjO1qTQMP5btNJdgp+8L3ye4WOyuj61VbWWV1RMy9GtbfjJdwI4QQsGaBDwFwFYA1RKQMq08SQswEcBkM3CdOEku0ydS5m/DM7I1YcM8IdGwW+9JMoZbjtIW/oJHJwq5OM2W2efSFExOc0uVGV1pejc0lkes/XjV9Gfw+wi9Pjg++YP02lohTXnzqgejJX23AzsMnIo4/UVWLrHQfiIh94IynSGQUyiLoGDxCiD84LZAesVjgyrkV1Q6F/4WI8uTMDbiwf1tnynUQJxS4UuW/m75Ed1q68mKtVdb4tNFS9W7pgWOVOKyywEvLq9H30W9w55juuO2sbmyBM57CLYPDMzMxnfDphtZhLG4P9ZnqRRoSzbKth7BLw4K1i/ISMMspcsrDs4LZEdXjA49+sQ7LdxzWPE+5BepbUVMrIl5Ah8olhT5l9ka8tWR72MrmDJPsuGVweEaBO7FaTKhOiMUlE6Fw9OatJ4hL/u3MElpWrfiyyhqsk1PFqkcXX1m81bR89WWqawNh24pLT4RNntpfVon9NnztDJNoEjkTMymodUBJhiqkaPNUjH12ATbuC5/oYhRy6GWiccPY6SkFp9Kr+jQ1gXALfPBT89A8N9O2LAyTLCTFRJ5E4sTsdScscLXyBpLPAneKaOq8sqYWRdsOWTpW6FjgWi6UA8fY4ma8S8In8iQaJyzwUEvPyQV8U9UCjyZNwcvfbcXFL32PtbvNU+McPF6FC15cHOGvrwkEErYcGsO4Qb13oTjtA3cyDa2d2YdeIpbZYwdN4rgB4JMVuzQjg2oCIu6x9QzjJm5Z4J5R4E5mHfygaCfmrt9nfqBFnLTmk4lYGt3Xa/fi/R92GB6j16mqqRUpW6dM/aTeW+BOKHDFqLvXZKq8XVLVhRJLk3t3mbHyBvQHkmsDImVWs2EYIIFT6ZMFJywyt1YpT1UXitvovZRrAsLRnOsMk2jq/UQeJyzwH7YdwpYDxrmuoyFVlU2iZqvXBgLsQmFSinrvA3figXZrWahUVTYLNpbgxfm/xP26NQHBK6EzKUW9jwN3evFiJ0lVf63ZFHq3YOXNpBr13gee6BW+jUjViTyJghU4k2qwBZ7ED3WqRqEkClbgTKrBMzGT+KHmKBRnSebeFsNEg1tRKJ5R4Mk8UMgWuLOw/mZSjXrvA3crhtsJqtgCZxjGgHrvA09mFwrDMIwR9d4HnswWOMMwjBFsgbMFzjCMR3FrUrOHFHiiJWAYhomOem+BswuFYRivUu8VOLtQGIZJJDeN6IIl958V1bnkkqb1jgJnC5xhmATSunEDtGqUFdW59d4CT+ap9AzDpD6xTMap94OYyTwTk2HiRYbfM49syhGLFc0WOCtwhkG6P0GrbDCxWeD1fSIP+8AZBkhjCzxhUAyOkHpvgXMUCsOwBZ5IYtHB9d4C5zhwhgHS2QJPGOwDjwG2wBmGFXgi8cVQ9ZxOlqfSMwzS2IViyp/P6YmLT2vneLmx+MDr/YIOPIjJMECaW6ZcCnHjiC44s3u+4+W65ceOBe8ocHahMAzSYunH1yPc0LVuWdGx4JnWwIOYDANkZ/gTLYIncEPX+r2owImoPRHNJ6L1RPQTEd0esu9WIvpZ3v60m4LW8LqTDIMmORmJFsETxOKv1iMZvVdpFo6pAXCXEGIFEeUBWE5EswG0BHA+gD5CiEoiauGmoGyBM/Wd7i1z0Tw3M9FieAI3lK1PLvSO0d2xdk8pZq/b5/xFbGJqgQshioUQK+TPZQDWA2gL4EYAk4UQlfK+/W4Kyj5wpj7j9xG+uWM4T+SxiBveDiWW+/bR3XDLyK7OXyAKbPnAiagAQH8ASwF0BzCMiJYS0QIiGuC8eHVwFApTn1EsSrcmhKQe9urp0fNPwnOX9jM8JtSqT5b7YFmBE1EugI8BTBRCHIXkfmkCYBCAewB8QBrDtER0HREVEVFRSUlJ1IJyMiumPuNWBMRpHZu4Um6isVtdEwYXoF2TBgD0QzVDlXZo+Q2z0vDHIZ1sy+gElhQ4EaVDUt5vCyE+kTfvAvCJkFgGIACgufpcIcQ0IUShEKIwPz/62Ey2wJn6jKIvnLT87hjdHS9c0d+x8pKJaOpJOaVLfi62TT4X7Zs2CC/Tp63AVz9yNk5q0zAqOWPFShQKAZgOYL0QYkrIrs8AjJKP6Q4gA8ABN4QEeCYm4w38LoUqKArJyeL7d2iMnEwrcQzew0415WUpdSCdpShntc4JrXt1lEuiwvOt3L0hAK4CsIaIVsrbJgF4BcArRLQWQBWA3wvhnpnMg5iMF/AToRbOt1VFeTjpSUnlWZ169XTrqK749ucSrNldCgD4/v5RaJIthWaqq0OtzkKterXCdiNs0QqmClwIsQj6L7QrnRVHH3ahMF4g3U+oqnW+XEV5OPkY+HzkaHnJhJ4C79EqD3eN7YGC+2YAALIz0pCV7pfPCT9JbTOSgQWeqDFN78zEZAuciROxxFq7tuCC0q138DFI8xFc6CwkBXqDvllp4TNZw90iyjbpk3ruSehMTLW1rva53z22Oz6+8QwbEkeHZxQ4W+BMvLh6SAG+njgsqnOra90ZrNFTKjGV6SM0bJCGSwvb4+aRXRwrNxnQM4jz88JfzmFuEVL5wFVVrTeIqfW9Q7OcuET4eEeB81R6Jk74fYSeraKLKih3w38Cd2YW+olARH7OHDcAACAASURBVPjbxX1wagfvhhM+fXEfPHb+SWHb9CxwdS4Zo2iVSB+4fvnqcuI1vOCZIWi2wJl4kYxJi8gFC9xvYFF6icGdm6F902zkZKYFXVh6P8en0qyhM1uF7E+qs8DD6zpUaavLN3OpuIV3FDj7wJkYePhXvdEkOwMT319pemwyKjOfjlKJhXAFnoQ/2iKKUr7w1LpFHPQUqHp7aB0E3VSyF0ytcvwa7pY6tC3waVedho37ygzljwVPuFCWbDmIyhoOBGeip0+7xrigf1tLx8YSy53nWly181EooWGE3lXf2u4KvfeRuncV+uLKSpfUYUWN5AZTvyz1ZmJqyaDkbR97UivcMqqbruyx4gkFPmN1MQD3JknEm87NcxItQr3DTtOJpvs7aXxPLL5vFBo2SLd9rhUU+Z3sh/o0rE8voiW7vgtFv5xMOUKlsloyFtUvy9DLRPq8VZZ9nJKOeUKBK4o7FSYerHv0bJzctlGixah32Hn5R9PMMvw+tG3cwLUZeXVx4M6p8LQU8YFryq5ngRvc3Ga50oSeC/q3AQDcpIrMMXrJqXfFS1d5QoErlZORAityZ2ekcW7zBGDHwlT7gy8b0B5TLzPOVKec45YlG/SBm3gSP77xjIiIDP0yU8MC1xp01vs9RgPU2Rlp2PDYObhrTA8AwE0jumLb5HPrzg3tsfhMLPA41acnNKIy6ylVVuRm/R1/7Cgo9bG/7tsG5/cz9p+7ne5VeUEIEyfKaR2bWJ6IFPo8OSF1olycdlwoZoO1Wen+COVcdx3tz1K59q7jFJ5Q4EplpaeABQ7w6kKJwI5rIyIkzIpiciHZlEbxll7+1RYjtkKtRD2FY+X3tMjLxNy7hmN0L1cX5dJFU4HrWeAx3KDwMEL9wdBYr2MHT2hEpW7OObkV7h7bPbHC2OCMLs00t7MCjz9WLOMceZJHh2bZts9VjnDPhaJY4OYM7iy1u7aNGxgeZzSzsG67+e8JCIEu+bkJM7BI47J6+jMW14aRBa7+XuPSjFw1nlDgCm0bN3A1JMcpmudmYOZtw/DyhELN/RzSHn/MFOvTF/XBDw+Mxhe3DMUZXcLT2tswwF2zvOzonfy8TGybfC6uHNTR8Lg0C1EoVi47tKtUX38+pyeGdYtYEsB1tJSyXn3FMsgcdm9NJu6cqHZnRq4aTyhwsmF9WGX67wsx5ZK+DpYo8eENg1H0wBj0btNQN9eyi1l3GR2UZ++fvztV875fMqA9sjPScEq7yAghK1YoBXNJu2uB28GspxceRqh9jNFl/T7CwntG4m8X9wEAtG+ajTf/dLptOWNFu26kbZ2a52DzE+NMjrV/HXU5yrfO+VKIcMuGWVFfx5ZMcblKjCiV46Tea92oATLSpJ/fuXkOXrt6ABqk+03OMseKBcYWePxRHrjxp7TGBSYDkpHnmh+jPM/J4ANXMDMUrIQRGuW5Jkjupsy02J+bWNCSXdlGJGWIbCgv2hCbDzzkc8Q+acu1wzpj/t0j4hYq7AkFrtSW2Qi8HdL9FGycfh9hRI8WjuRbseJjSxYfeNEDo3XdPKlGLHk/kskHrsfUy/rhs5uHhG1TDIXerRvirjGRY0dpIf4E3Z6DwWXt/NR7zu5h/WAAV5zeAc/81loPOTMtUo0FMwvK3/97y1D8ZXyvYO7vaDCywJXm5SdCpzhO1POEAicXphH7fRT2lgacGXiw8oZPEv2NZjkZGN2rBR674GTcMdo7g8OxYtfNoTysa/96Nr66XTvNrFLk8aoa2/I8fXEfPPKr3obHmEl8fr+26Ne+cdg2xVAY3asFbj0rcuzIylR6o+ZsdRWabZPPxZWnG/vj1b3fm0Z0wYge1tbQ1bqfpNpX0DwH157ZObi/RV4m2jSy5+YwMgKU68Q7nN4bCtyFSkn3++pusvzJCdeGNReK+YX6avhinYbkdKJXDeqIBhmeaApRE0sbUs7NzUxD42ztqfJKGzq9U1MA9iadFTTLwR9MVjWPRn6lPeuGCFoaxDRwodiRyeTYGbcNRdEDo8PkMXvRXn9mZzx7qbaVbubSWjrpLCy+b5SxUDplSp91LPA4x8Kn9lNrQJo/0gJ3pFwHFHjRA6Nx7zk9nRLJEvHqFSRqskcsg4uhyi1NL4xBPuSRX5+EhfeMtJUTRW+C2sbHx+G5S/uFyWDHjaj4wK3UuX4Yof1zojk2NzMtbAISkXmvo3+HxvhN/3aa+5QXj+6LycILQo1RMiu3Z+Lq4QkFrig8J+tG8lU7HzlgZdKH2XRoyfqQPndunoM7Rnc3jemNlXh5dRKVziaWy4YqwHQdZdu1RS4AKSFSh2bZtn6n3ks/I80XbE/BNmrjRinPjRVZogkjtLOQr5li01ogwewcO7lJnCD8RS59biGv8OOGMWhJpvheLkrkRuvk241ClKSTda5+GL+540y8+aeBYdtMw7uo7uHIyUzD7aO7uZYkSSFeFniicm7EctnQW6plzX5378iIFW3sXM/IQq4bHLVenoKZC8UKTi2kYiZBpEVrfpK1noVz7S00QCE7Iw1PX9QH/71lKIC6ds0WuAZBC9zBMn0hXTQnlaP6BnZvmYch8sSQ+8dJbhGzR4KIgg9sNOFj0eBkhI8RCXOhxNB6QpWAP2gR1+1v3zRbfYqt6xnNYFRbdnbukp2eq57iqTFYytCOrlJ+48CCppauTyDTl5aRslRChJ2cEame8XnJgPZoJQ+Eup0LRw9PrMgjXLPA5YdRfthuHtkFHy/fjb1HK6IuV8uf6fNRWFYzs/hcH4WOakv/8/MysevwiajlMiNuPvAEWeCxvDdCz83OSMNTF56Czs1zcOm0JY5cz9gC17bs8rLSUFZhHPFi57nRM2JqDEb2zUr97OYhwWMy0nzY+tR4fL/lIK54eWnk9Snyu5n1bOSuVEILnVwIxtBlEwxJduxylvCIBS79V+rv6Yv64N5z7MWVqvFRZKO55+yeWDLprJjKtRYHbrzfF2qBy9umXVWIpy/qoxnzasaUS/oGB8P0iNfs0IRFUMakwMNPvnxgB02rO+xyNl5U6SrtmZPhx4WntpXLCZfh92cUAABmTTwTqx4aa1huIBDuAzd6qUSTv1rvN867azje/NNA9GvfGH1DQhvJwK+tLiu0FwoA71xzOi46NXzA0khkZXKRk0sxGl5PXs0nI4rnMxa8YYGrHvtLBrQHADz99c9Rl0mg4BvcSfeBpUFMUws8crQ7Py8Tlwxojwc+W2tbpj7tGqFrizzD9SDLKu3HL0dDtYNd2ltHdUX7Jtlo2SgLB8oqcdeHq3SPzfTbn8DhI+llq6V0nOwNqldv+enRcyKuo1jp/do3DuvNGRFQWeBpfh+qdCxS3eiaKOicn4vO+bma+/QVuPq4cNfSGV2bY+ba4rBjjIyllg0zcf3wzhFKPxaM7nm3Frl48YpTMaybtdh1p/CEBX6hHCp0Vq+WjpVJvjqrw611BvUwMwqI6ixVdaNRv2xe+UMh3r12kIlM5re5g2xRur2SiJMKvF/7xrhkQHsM756PkT21U5mOO7kV3rnmdDTSid82IjibT6NKTP2zBlXeV1bCirVmVOeK4o5uEFPxgUsnpxsUEq9c+7o5VyK+U8TYQHWNao1KI9cTEe4f1wvdW+ZFI2YYN46QVuYxiu0nIpzbp3XcMzJ6QoGf0q4Rtk0+19EpqoS6N7iTCtyKBW7uA6egTOrS1Mp/VM+W6N2moWF5TXIyTGW6YmAH/PeWIZZnv0WLk3lgwnoqOsc0zErHGV2jy5AXjCzQuKdmLhLFJ/rrvm3w2AUnax6jKG4jBR7L4JhQhRE2zdVvB/FSPFZ/hlaKWPXLP14Dhvee3QMbHx9nLS98nPGEAncDH1HQunHSL+vERB4fhT584eVpnat3yd6tG2LVQ2PRyMKkEiJCn3aNXV9JpLmBErFLqKhWu+b2LqCUbb9cZf/ff9sHvxvYIWyfMo7hV7lHtKiLA7cgrwrlZamU/841g/DQedpT9uO1hqNe+1K36tCj2jWR5kBcKLtDmsoGSbyMXSKKu2/bKskpVRwIU+AOmuBWrAIrE3mClqqqOC1RjZSuXdeBm8/xyxMKMVEn58qfz+mJ164eYKu8sN/tgtyhCYoi95lESCh+Z58vwnLLldMMW7k3yrVjSSernNm+aTb+OLST9nVUMo7okY8L+rVBYccmmseHlmsHvd+hfoEobr/v7h2JGbdJ+WeGdmuObZPPDSp0L6/j6RQpr8ALmmlHCxBZc3eYoU7C44QFLvnA7cyii+462mW591CM6d1SdybjaR2bYIBOjLAeYfn1XVHgivWrpcCNz1V2ax13QX8pwuTta07H/eN6onG2ea9Ez0qfNfFMfD1RO8GWnYk86rj1FnmZeO6y/sHoCqdQS3Jh/7b44PrByM4Ij6dQfm/7ptkRPUg3woq9iieiUGLhmUv6IjPNj/OeXxS2nciZQcycTH/Y6htOZCOkMB+4eXl6DTnZFDig7wOX4n7tlWXkA2+Rl4n9ZZX2CtQpX9uFYq7BtWKZNz0xLuhv7tgsB9cP7xLc98IV/bFyx5Gw49UDkWp6tNIfpNNzw2nRKDsdD5zbC4fLq/Di/F9cmxeglqVRdjoGdop8cRsZQnXGDSvwlLfAs9L9msnVpZleig88+tZa0Cx8YNXqGoJmBBV4DG1UrSwfOLeXZl7oUPSu1yU/B1cPKYhemKBM2r/dKEZYoVPzHFwfkhI0zIMS8mXOncNxm5w+1az+zujSTDdXtToGOxT1TNmIc6H9MjcaLDyvTxs8oPJR28lnosbuudcM64yOcntW2o7TijwyXFBbOKPesRPPRqqQ8gpczyL2UV3oVCyNtFfrhro5ovWwpMBtWBm6FrhKg18zrLNmXuhQ9F5AH1w/GHeaKH8r6P10K0qmW4tc3D++V/B7aO8k9PSuLXItP9zvXDsIN4/sqiOTfiFKPelZilaSMVlBGS+JZgarOg7cCraMmih+nla8t12UyTmswOuBC8Vo5pfycMTSEIgkJW4HKy8M9exTI/QegmgSEemV5ZRC0hswDk0uZhUjH3jdLMToZVbK15K4zset176cyfti5kKxdm749i9vHap7jnKosGCBR/Pr1D2QaMahFJdlTkbKqy9TUr4GdLtoBJRVVAMAmlgYRNIjmkZszYVi/cHVO6aDyXRvLfTry76C1ULvl/uJTP39EdZb2Ioy4Tur5SRMscQ3K3Whdb+ULXpKmkCO5H2ps6Ltn6ssrNtUNQ/AaL1GZQKSW6kV1KkgoqmjymrJ8tJbNLw+kfI1oBcDTUQ4tWMTXD6wPW4dZexWMCI6y8j8GL2JPFr4CDipTUNc0K8tnpi5HgDwxh8Homdr+7PQ9K5Hvtgy+inozQr1+SzEVquub2SB18jdq1jim6dNKMSri7eiqcYLPktWRHeN1fefOxHlZGcgUs0do7ujZ6s8jNKZparFyB4t0Dk/BzfpuJViRR1PHU0v5dWrB+CzH3c7OqfAq5gqcCJqD+ANAK0ABABME0JMJaJHAFwLoEQ+dJIQYqZbglqloFk2th0sx62juuJXfdsgPy9T99h0vw9PXdgnpuspz9U9Z/fAZz/utnROXpb5e9NOGCERYcZtw7DvaEVQgZ/ZPboZlXovJL9DFrheGKGPzF8Pkdev26AoOMXaVCzwtBgs8NM6NsFpOnHQaX6fYU4SCplnoBDNLFfFDRZNqpKMNB/O79fW1jmNszMw764R9i9mEfULPJoXU6/WDW27LVMVKxZ4DYC7hBAriCgPwHIimi3ve1YI8Q/3xLPHtsnnYvJXG/DSgl/QIMPvSB4EM5QGePPIrrqDYWpenlCImWuK8fiM9brH1PnnrTdwJxSskQ/cifL1ZrT5feY+dqMBsIw0HyaN74nTOkohaTVBF0piRrrUUShbnxofVTnRDES6QWaaLyI1azS9T3VceaLyw6cKpu91IUSxEGKF/LkMwHoA9l7rcaRW1nxO+B9fuKI/HjZZLTya9temcQNcM6yz4TFDuzXH2N4tTa8fLkvsv9loUM6J8vV6H9HEgasVyHVndglazHUulMQEWvl84W0wmjUYAWBw52bISveZthe3UHqCL08otJwF0YiGWel4eUJhMHRTPdNz1sQz8Z8JhTFfp75gywdORAUA+gNYCmAIgFuIaAKAIkhW+mGnBbSLkoDeiTf7aR2boHWjBvjrF+t0j3Eyd8j4U1oFk/RnpfsxzWZDdkKSe87pgeraABZuOoADx+omwviIHMmXMbJHCzx0Xm+8u2wHNu0/Fla+1QRRCkYvaaUdxCvLnhpCpAslGvLzMrHhsXEOSBQbTvYAxvRuibN6tsCY3i0jesk9WuUZTk5iwrFsnhBRLoCPAUwUQhwF8C8AXQD0A1AM4Bmd864joiIiKiopKdE6xFHO6imlnB3UuVnMZVlptE72bKde1h9v/un0qM934iFrnpuJKZf2Q4ZK8fl90VmQaogIfxzaKSKCwJKyUx1iJI6ylFbCXCgUnd862RAxRMEY4fNRXFycqY6lJkZE6ZCU99tCiE8AQAixTwhRK4QIAHgZwECtc4UQ04QQhUKIwvx891KV5mRIOUmUhDdGoVJWsaKvnIjMUIjV7ePky0S9lJbbrkpLL0vVdyOlHxzETJAWJXImjDBZ0Hp5xyuDIaOPlSgUAjAdwHohxJSQ7a2FEMoSGb8BYH+pGId45Fe9MdSFlTCszYJ08HoxFuakO0cd+2xU9tkntcSsn/bFdD0rv10tg9H9qU60BQ5nwggTTTAppuqn3DSiSzApF5M4rPjAhwC4CsAaIlLW5JoE4HIi6gfpHm8DcL0rElrgD0M6uVJuvF0odhjWrTmOqZZBc1IWvbUE/zK+FzbtL8MHRbuC2wo7NsW/rypE4eOzceBYlaXy1aVH4UExPKc2OBaSKAs8cQs4u4H6l9x7Ts+EyMGEY6rAhRCLoD0+lvCYb7exYr0lKrxLy1fupCx6q5Ffe2ZnfL22OEyBO3HZaJSdUa9gTO+WeO+Hnbpx3G6jDiP0LKpmMLCgaVRZLhl3SPmZmLGgrGzdIN2PUzs2xuLNBwEA6x49Gy99+wv+b95m11ewsYOTkhit5q1+UURTB7UB+zHF6kOMFORZvVpiy5Pjdd0YGX5fMNTQDZzKHZMsKPfngxsGJ1gSJhRW4AYoFvj6x87BscoanPzwLACQks/LDTqZHlE7CqNLvvH6onoWOBAZmlenI61fX71ArRVr1Y4LBTD2Qa96eKzp9WLBqWRWiSaWVMuM+7ACNyDUKox4FmPIUaHw3nWDsHLnEUz+aoPlc+4Y3R0fr9iluc+qKKsfGWu4wjZgzwKPpg4iF6gN33/P2T3QtnEDTHx/ZXBblmr1o1jqvkGG3/ygGHAqDjxZcCu5FRMbrMA1uOfsHvhoebiSVIcL2kn3qsegzs0wqHMzbNp3DPvLKiydc/vobrh9tHbyLauyNMwyX4vRSIGrFZOZnvry1qERGfGqVe4L9fVuHtk1YpBWrcCT2UORKhZ4n3aN8cO2w2yHJymswDXQymuiVhZK19KJR/SZS/o6UEr8BlTVisnMf60Vk9+oQTp24gRuGtEF//z2F83UoErGv0YN0lF6ojoiFWky+5hzM9MM3VBe4b5xPTGgoAlO11j2jEk8rMAtEqHAlRlqSWRlKZK4rdfUESOKYrVz3ZcnFGLu+v24clBH3ZA0JePfs7M3YurcTREWeDJbuI//5mTX1pWMJ+l+H845uXWixWB0YAVuEbW154QLxWl8REj3Ex4413oCrGhQFGfbxg1w0WntoprQ0bpRA1w5qKOlYytqpBVYvORCad2oQaJFYOoBnlbg+XmZKIlx5XGrRI5hKi6U5NEiPh9h0xPRpS3V48HzeuPnvUfDtikKPDczzZF1Ms1QVmDxkguFYeKBpxX44j+PitukArWyEMHtcbl8wrj41HZolB0+6Km/jJg7KPc4O8O5KBSGSQU8rcD1FgdwA7Wv24lFc72AVjpWs5XLBxQ0QemJamzcd0xzv13uHNMdPiL85tRwV00qTVVnmGhIgYSXiSEZfeBuoKXAzXJsv3DFqZg18UzHZGicnYFHfn1ScGasAnHrZeo5/AhESX2ZoZaukQxKSdGqFyYnhBRaaDZZKFZSvffDMGZ42oWSSOoS3ae2EtEKk+zUPAfDujXHDcO7hG1v1CAd+8sqg+MCKx8e42ooHeejZuo79UaBP3b+SejUPBdXTl/qSHkdmmYDANo0znKkPC/h95FmNsRXrx6AWT/tQ4uGUp1kZ7jbvFiBM/WdeqPArxpc4Gh5fzijAN1a5mJo1+aOlutl2jXJxp+GupObXYtknsjDMPGg3ihwJ3jjjwORn5cJQHItDHNhFSDGOsmUypdhEgErcBuc2Z0VNsMwyQMrcEaTBfeMwM5DJxItBsMwBrACZzTp2CwHHZsZL/rAMExiYQXOeI6vJw7DlpLjiRaDYRIOK3Am6RhuMtbQs1VD9GzVME7SMEzywgqcSSq2TT6Xl+9iGIvwVHom6eDwQIaxBitwhmEYj8IKnGEYxqOwAmcYhvEorMAZhmE8CitwhmEYj8IKnGEYxqOwAmcYhvEorMAZhmE8CitwhmEYj8IKnGEYxqOwAmcYhvEorMAZhmE8iqkCJ6L2RDSfiNYT0U9EdLtq/91EJIiIV/dlGIaJI1bSydYAuEsIsYKI8gAsJ6LZQoh1RNQewBgAO1yVkmEYhonA1AIXQhQLIVbIn8sArAfQVt79LIB7AXACZ4ZhmDhjywdORAUA+gNYSkS/BrBbCLHKBbkYhmEYEyyvyENEuQA+BjARklvlLwDGWjjvOgDXAUCHDh2ik9JhWjbMTLQIDMMwMUNWlq8ionQAXwKYJYSYQkSnAJgLoFw+pB2APQAGCiH26pVTWFgoioqKYpc6Bg4cq0RWuh+5mbyaHMMw3oCIlgshCtXbTbUYSetbTQewXggxBQCEEGsAtAg5ZhuAQiHEAcckdonmuWx9MwyTGljxgQ8BcBWAUUS0Uv4b77JcDMMwjAmmFrgQYhEAw1VmhRAFTgnEMAzDWINnYjIMw3gUVuAMwzAehRU4wzCMR2EFzjAM41FYgTMMw3gUSxN5HLsYUQmA7XG4VHMAyRaTnowyAckpF8tkDZbJGqkgU0chRL56Y1wVeLwgoiKtWUuJJBllApJTLpbJGiyTNVJZJnahMAzDeBRW4AzDMB4lVRX4tEQLoEEyygQkp1wskzVYJmukrEwp6QNnGIapD6SqBc4wDJPysAJnGMbzyGmvk4p4yORpBV5fb1o0JKNcLJM1WCZLpCdaAA1cl8lzCpyITiKiEQAgksSBn4wyAckpF8tkDZbJGkQ0mIg+BPAPIupNRP76JJNnBjGJyAfgBQCjAOwAsBTA50KIIiLyCSECLFPyysUysUwuyNUCwFeybO0BtAVQJIR4mYgoES+ZeMvkJQu8CYA8AL0A/A7AQQB3EVFuohpQksoEAI0B5CaZXI2RfHWVjPePZbJOXwAbhRCvAngGwCcAziei7kIIkSA3T1xlSmoFTkQXE9FN8teGAAYDyBZClAD4GMAhADfLx8blZhHRqUTUXf7aCMAZiZZJvlZnIsqWvzaFtBReouuqIxFlhciUDPdvKBF1lb82ThKZuJ1bk+lyIvorEf1a3vQjgEIi6iKEOA7gBwBFAK4H4uPmIaILiGgSEZ0rb1oZT5mSUoETUS4RfQzgbgCHiShNCLEVwGIAE+XDiiG93foTURu3bxYRdSKiGQBeBPAmEY0RQmwB8L9EySTL1ZqIFgJ4C8DnRHSKEGIzgG8B3JkIuWS/32cAXgPwhWx9bAawBImtq34AFgK4nIgaCiF+AfB9omTidm5ZJiKiGwDcC2AbgL8T0TUAjgF4A8Dt8qFHAMwBkE1ErV2WKV9u43dCepm9SkQXh7zgbo2HTEmpwCH5jvYJIQYJId4FUCtvfw3AECLqLISoAbAPQAWABnGQ6W4AK4UQgwF8DuCP8vZXZJk6xUsmlcVzKYAfhBBnAJgL4D4iOhVSXQ2KV10pMhFRTwD/AjBfCDESwCpI/kAgznWlYRm2BTAbUrsfngQyJUU7V8mUNO1cQX5BDAYwWXZN3AxgBICzAMwE0JWIRsvunIOQ7nOpy2J1AbBYCHGmEOIlAHcBuEPe9y6AnkR0ltsyJY0CJ6JWVDda2xdAO3n7TQAeJqKhANZBsgT+AQBCiLUAOgKodFEmJRToBIBq+XNDAOvlrvhiSN2kuMgkk6X6nC5fezKA/ZAa9j5Ig01Px0kuRaZSAPcJIabK3x+DZH3kA1gGYEUCZFI4DGATgACAAUTUAFJPJZ73L1SmPkiCdq7IJCvy40iCdk5EE4hoOBE1lTetB9BW7qXMAbAWwCAAJQDeAfCcLOdZkBZhz3BJphGyq3I5JOsfst5aJ/8BwBoA7wGY6rZMCVfgRHQWEX0Hqcv2L3nzCgB7iOgVSG/eIwD+AuB8AFMAtCCiF4hoLaT84qVO+uFUMv2fvPk7SG/6HwGcA8APqeEMB/AUgNZuyiTLNYaIZkPqQl4mb94K4CARdZC/vw/gFEh+y6cAtCOi512sq1CZLhFCFAshvg+5xikAqoQQJUKIYwAejaNMT4fUkyLLCkh5KLIATAJwMaR6ahNHmS6XN68AUJzAdq6+dwLAIgDdEtHOZVdJayKaD+D3kAZMnyeihgB2AmgBQBm/eA/AyQCaCSHeAvA2gPsAXAbgXiHEEZdkugLAywAaCSH2EZFfCFELaYC3EQAIIQJCiNcgKXjHZQpDCJGwPwDdIVmJF0O6OV9BGnzrAMlKWw4gXT72KgDT5M8tIA2q/DoOMn0NYJi8rwuAT0KOfRDAC/Lnlm7JJJffVZbrfAD9IT1UtwBoDeBVAOehLiz0NQAPyZ9buVhXapneAjBJ3qfct7FKHYWc5+b905LpAXnfrwBcCaAzJNfOUQATElBPb0PqcqdBilRIRDvXak93y/t697sK4AAABhRJREFUxLudA/DL/7sDeEv+nAbgnwBeh9TLfEWun0Yh7fyJkDIy4iTT80r9hBzzBoBLlLbklkzqvzTEGZJiSiEk31A/AMuEEB/Jb9njALYJIXYT0UwAhQB+C6lxrQJwEUlxp/shuQriIdMxAFuJKEP+vJOIegkh1gOYD2AiEZEQYh8kt4VjqOQ6HcByIcTn8r65kB7+1yG5JobJ8n0L4EtIL0IIIfYC2BsnmeYBmEJE/5HvESB1H5fK+x8E8KoQYhfcu396Mv0L0svuWgAPA5gB4BtI7h1/AurpGQDTIfmZ+yP+7Vwt0xxI9fQmpEG5uLRzIkqD1Cvzy898Q8hjAUKIGiK6BdKgaW9I9XMBJLfTU5BcYf9TyhJCVMVJptsgeQiGCyEWyKcpeuJRABcS0TlCiF1OyaRHXF0oRHQ1gF2Q/KIAsBrAaUT0MiS/UQtI3bn/CCG+BTAVUrzpnyF1mRbJ5TjZjbQi098gWQIHADQDcBsR3Q7g35BGmB1HQ641kKInCuTvaZBG5P8my7ELwDNEdB+A5yAp8njLlA7gF8h+Uvk+nQZp8GsBgJ6Q/NDxlmkrJKX9EaT7NVgIMRGSH7UMkn8ynjKlyTI9LYRYCOl+xbuda9XTFnl/GaSwT1fbORENh9T7aAJgs3ztagAjiWggEHzZPArgb0LyfU8DMJSIlsrnfZsAmYQs0yPyOX5Ig70fQVL2I2UjxX3cNO9V3ZFcAJ9BCvlZAaCHvD0fwD0AbpS/Z0FSlIPl76dDiqEcnGCZDkLyuXWQj38dwKA41VVPeftzkEa4F0NyDZwCaRQ+X94/DsBDAIYmWKYZkEbdST52HoD+CZbpKwDNVeenJ8G9ayXvHxDHdm5WTzmQfLq3utzOhwG4KuT7PwHcCOAPkHoHgGRktgLwIYACeVtjAG2TQKYPIA3kdpHr8lQ3ZDKUN64XAzrI/ycDeDekMqZD9jPL214AcF6SyfQigPEJqqv35c9+SJbRUPl7e0h+wKwklMkP6cXnaqO2IdOrADKV+5skMiXjvXsdLvttQ2TKBpCJOj/y7wA8JX9eCeBW+XOh8mwmmUzvxUMmo7+4ulCEEDvkj88B6EJE44TURdoMYBoR9SCiSQCGAvgpyWQaAimUKS6o5OpERGcLabS7VAixSN53A4By1IV9JZNMJISoEEKsSBKZTgCokc9xdfq3x+/dcdTFo7stU7kQolKWAwDGQAoLBICrAfQioi8h9RJcbUdRyrQcSHBmxkS9OSB1F78L+f4PSKPz7wJozzJFyLUg5PtASINfwS44y8QyeVUmSD0BHyT3TVd5W1dIrpKhcMld4jWZtP4Sko1QHmEPENFHkEbZyyH5k9YIIU7EXaAklUlDrmJIEyfmANgkpKngLBPL5HWZlEku/wHwKaQBwYOQ3BVHWSZ9EjKRR25A2ZAiPC4BsEMIsSyRijIZZdKQ63JZrq8T9bCxTCyTCzIJSKGUv4OUW+RTIcTvE6kok1EmLeIeBx7CTZD8WmOEEG5OO7dDMsoEJKdcLJM1WCZr7II0C3UKy2SdhC3oQAlMBK9HMsoEJKdcLJM1WCbGTTyzIg/DMAwTTsKTWTEMwzDRwQqcYRjGo7ACZxiG8SiswBmGYTwKK3AmZSGixiQvFkxEbeTJKwyTMnAUCpOyyKlSvxRCnJxgURjGFRI5kYdh3GYypARlKyGthdlLCHEyEf0B0sIAfkgpgp+BNG36KkhTy8cLIQ4RURdIWSjzIaVWuFYIsSH+P4NhtGEXCpPK3AfgFyFEP0j53UM5GdL6hgMBPAGgXAjRH8D3ACbIx0yDlPviNEirtf8zLlIzjEXYAmfqK/OFEGUAyoioFMAX8vY1APoQUS6ktR8/DMkWmhl/MRlGH1bgTH0lNLdFIOR7ANJz4QNwRLbeGSYpYRcKk8qUAciL5kQ569xWIvotIKUXJaK+TgrHMLHCCpxJWYQQBwEsJqK1AP4eRRG/A/AnIloFaYWo852Uj2FihcMIGYZhPApb4AzDMB6FFTjDMIxHYQXOMAzjUViBMwzDeBRW4AzDMB6FFTjDMIxHYQXOMAzjUViBMwzDeJT/BwWgjJz0gY7lAAAAAElFTkSuQmCC\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": 9,
"metadata": {},
"outputs": [],
"source": [
"climatology = area_weighted.sel(time=slice('1980', '2010')).groupby('time.month').mean('time')"
]
},
{
"cell_type": "code",
"execution_count": 10,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"[<matplotlib.lines.Line2D at 0x2b339cde0978>]"
]
},
"execution_count": 10,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXoAAAEGCAYAAABrQF4qAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3dd3xUZd7+8c+XFEJCh9BDAtKlE5ogay/YG6CgoCLsIvZdly3Ps/vbXVfXtupaVgUUaYKKvaKyKkpLIDRDk04oCSUEQkJC7t8fGZ7NsgmEMJMzM7nerxcvMuXMXMeEy5N7zn0fc84hIiLhq5rXAUREJLBU9CIiYU5FLyIS5lT0IiJhTkUvIhLmIr0OUJqGDRu6pKQkr2OIiISM1NTULOdcfGmPBWXRJyUlkZKS4nUMEZGQYWZbynpMQzciImFORS8iEuZU9CIiYU5FLyIS5lT0IiJhTkUvIhLmVPQiImFORS8iISE7t4Dpi7awfneO11FCTlBOmBIROW51RjZTF2zhvbQd5BUU0axODB/fey714qK9jhYyVPQiEnSOFhbx6aqdvLFgC6lb9hMTVY1rujXnnDYN+NVbK3hgdhqTR/amWjXzOmpIUNGLSNDYmX2EGYu2MnPxNrIO5ZPYIJbfX9GRm3olUCc2CoCDeYX8z3ureOmbn7j7/DYeJw4NKnoR8ZRzjgUb9/LGD1uYm76bIue4oH0jbu2fyKC28f911D6ib0uWbNrHU1+spUfLupxzVkOPkocOFb2IeCInr4B3l+1g6oItrN9ziLqxUYwe2IoR/RJJqB9b5nZmxqPXd2F1Rjb3zkzjk3sH0qh2TCUmDz0qehGpVOt35/DGgi3MWbqdw0eP0bVFHZ64sStXdWtGTFREuV4jrnokL43oxTXPf889M5cxfXRfIiN0EmFZVPQiEnCFx4qY++NupizYzMKN+4iOqMaV3ZpyW/8kuifUrdBrtmtci0eu68yDs5fz9Nx1PHxZB/+GDiMqehEJmD05eby5eBszFm1l18E8mtetwcOXtWdocgINalY/49e/vmcLlmzex4v/+onkpHpc0KGxH1KHHxW9iPiVc46ULft5Y8EWPlu1k4JjjnPbNuTP13bmgg6NiPDzKZF/uOpslm/L5oFZy/n43oG0qFf2+H5VpaIXEb/IPVrI+2kZvLFgC+k7D1IrJpJb+yUxol9LWsfXDNj7xkRF8NKInlz53HzunrGMt8b2JzpS4/UlqehF5IxsyjrM1AVbeCt1Gzl5hXRoUou/XteFa3s0Iza6ciomsUEcT9zUlZ9PW8pfP0nnj1efXSnvGypU9CJSYU9+vpbn520gsppxeZem3NY/keTEephV/ozVyzo35Y4BrZj8/SZ6J9Xniq5NKz1DsFLRi0iFfLA8g+fnbeD6Hs2ZMLgDjWp5fy77hMs7sGzbfn79zgo6Nq0V0CGjUKKBLBE5bWt35fDrt1eQnFiPx27oGhQlDxAdWY0XbulJVIQxbvpSjhw95nWkoKCiF5HTkn2kgLFTU6gZE8mLw3sG3QefzerW4O9Du7N2dw7/+/4qr+MEheD6DolIUCsqcjw0O43t+4/w4vCeQbv0wHntGzH+/Da8lbqd2SnbvI7jORW9iJTbC/M28GX6Hn53RUd6J9X3Os5J3X9RO845qwH/894q0nce9DqOp1T0IlIu/1q7h6e/XMc13Zsx6pwkr+OcUkQ149lhPahTI4px05eSk1fgdSTPqOhF5JS27cvlvjfTaN+4Fo9e38WT0ycrIr5Wdf5xcw+27stlwjsrcc55HckTKnoROam8gmOMnZqKc46Xb+1VaZOg/KVv6wb88pL2fLxyJ1N+2Ox1HE+o6EWkTM45fvfuKn7ceZBnhnUnsUGc15EqZOyg1lzYoRGPfJLOsq37vY5T6VT0IlKmaYu28s7S7dx3YduQXhmyWjXjqSHdaFQrhvEzlrH/8FGvI1WqUxa9mSWY2TwzSzez1WZ2n+/+WWaW5vuz2czSyth+s5mt9D0vxd87ICKBkbplP3/6cDXnt4/nvgvbeh3njNWNjebF4T3Zk5PHg7PTKCqqOuP15TmiLwQecs51BPoBd5tZJ+fcUOdcd+dcd+AdYM5JXuN833OT/ZBZRAJsT04e46an0rRODZ4Z2uO/rtsaqrol1OX3V3Ri3tpMXvrmJ6/jVJpTFr1zbqdzbqnv6xwgHWh+/HEr/vh9CDAzUCFFpPIUHCti/IxlZB8p4J8jelEnNsrrSH51W/9ErujalKe+WMuCn/Z6HadSnNYYvZklAT2ARSXuPhfY7ZxbX8ZmDvjCzFLNbMxJXnuMmaWYWUpmZubpxBIRP3rs0zUs3rSPR6/vQqdmtb2O43dmxt9u6EpSgzjumbmMPTl5XkcKuHIXvZnVpHiI5n7nXMlpZjdz8qP5Ac65nsDlFA/7DCrtSc65V5xzyc655Pj4+PLGEhE/+mB5BpPmb2Jk/0Su69HC6zgBU7N6JC+O6Mmh/ALunbmMwmNFXkcKqHIVvZlFUVzy051zc0rcHwlcD8wqa1vnXIbv7z3Au0CfMwksIoFxfEXKXon1+N0VnbyOE3AdmtTmz9d0ZuHGfTzzZVkDEuGhPGfdGDAJSHfOPX3CwxcBa5xz28vYNs7Mah3/GrgE0HJyIkHmYF4BP5+WGrQrUgbKTckJDEluwfPzNjBv7R6v4wRMeb6bA4BbgQtKnE452PfYME4YtjGzZmb2ie9mY2C+mS0HFgMfO+c+81N2EfGDoiLHg7OWs21fLi8O70njIF2RMlD+dE1nOjSpxQOz0thx4IjXcQLCgnHth+TkZJeSolPuRSrD81+v58kv1vGHqzpx+4BWXsfxxKasw1z1j/m0aVST2SF6cXEzSy3rFPbQ2xsR8Ztv1mXy1NzQWZEyUFo1jOPxG7uStu0Aj36a7nUcv1PRi1RR2/blcu/MZSG3ImWgDO7SlFHnJPHa95v5ZOVOr+P4lYpepAo6viJlkXP8c0TorUgZKL8d3JHuCXV5+O0VbMo67HUcv1HRi1Qx/7Ei5dDuJDUMzRUpAyE6shovDO9JZITxi2mp5BWEx8XFVfQiVczxFSnvvbAtF3YM3RUpA6V53Rr8fUh31uzK4Q/vr/Y6jl+o6EWqkOMrUp7XPp77w2BFykA5v0Mjxp13FrNStjF/fZbXcc6Yil6kisjMyWfc9FSa1InhmaHdw2ZFykC576K2NKxZnVe/2+h1lDOmohepAgqPFTF+xlKyjxTw8ohk6sZGex0p6FWPjGBk/0S+WZfJut05Xsc5Iyp6kSrgsU/XsCiMV6QMlOH9EqkeWY3J8zd5HeWMqOhFwtyHyzOYWAVWpAyE+nHR3NCrBXOW7SDrUL7XcSpMRS8SxtbuyuHhKrQiZSDcMaAVRwuLmLZwi9dRKkxFLxKmquqKlP7WplFNLujQiKkLtoTsefX6zouEoZIrUr5wS9VbkdLfRg9sxd7DR3k/bYfXUSpERS8Shl781wa+TN/Nbwd3pE+r+l7HCXn9z2pAx6a1mfjdJoJxxd9TUdGLhJmSK1LePiDJ6zhhwcwYPbAV6/cc4tsQnEClohcJIz9lHuK+N7UiZSBc1a0ZjWpVZ2IITqBS0YuEibW7chj68kIiq5lWpAyA6MhqjDwnie/WZ7F2V2hNoFLRi4SBVTuyGfbKAiKqwZtj+mtFygC5pU9LYqKqMWl+aB3Vq+hFQtzSrfu5+dWFxEZHMntsf9o0qul1pLBVLy6aG3u14L1lGWTmhM4EKhW9SAhbtHEvt05cRP24aGb/vD+JDXQkH2h3DGjF0WNFTA2hCVQqepEQNX99FiNfW0yTOjHMHtuf5nVreB2pSmgdX5OLOjZi2sLQmUClohcJQV+v2c0dU5aQ1CCOWWP7a0JUJbtzYGv2HT7Ku8tCYwKVil4kxHy2aidjp6bSvnEt3hzTj4Y1q3sdqcrp17o+ZzerzaT5mygqCv4JVCp6kRDyftoO7p6xjC7N6zD9rr5aV94jZsboc1uxYc8hvlmf6XWcU1LRi4SI2Uu2cf+sNHon1WPqnX2pHRPldaQq7YouzWhcuzqTvgv+tepV9CIhYOqCzTz8zgoGtmnIa6P6EFddk6G8dnwC1fwNWaTvPOh1nJNS0YsEuYnfbeR/3l/NRR0bMXFkMjWiI7yOJD639GlJjagIJgX5FahU9CJB7Pmv1/OXj9O5oktTXhzei+qRKvlgUjc2mpuSW/B+2g72HMzzOk6ZVPQiQcg5x5Ofr+XJL9ZxfY/mPDusuy4cEqRuH9CKwiIX1BOo9JMjEmScczzycTrPz9vAsN4JPHlTNyIj9E81WLVqGMdFHRsH9QQq/fSIBJGiIsf/vr/6/y7m/dfrulCtmpYaDnajB7Zif24Bc5YG5wQqFb1IkDhW5JgwZwVTF25h7KDW/PHqs1XyIaJPq/p0aV6HSfM3BuUEKhW9SBAoPFbEg7PTmJ2ynfsubMuEyzvooiEh5PgEqp8yD/PNuuCbQKWiF/HY0cIixs9YxvtpGTx8WXseuLidSj4EDe7SlCa1Y5gYhGvVq+hFPJRXcIyfT0vls9W7+N8rOzHuvDZeR5IKioqoxqgBSXy/YS8/ZgTXBCoVvYhHco8WMnpKCl+v2cMj13XmjoGtvI4kZ+jm3i2JjQ6+CVQqehEPHMovZNTkJfzwUxZP3tSN4X0TvY4kflAnNoohyQl8sDy4JlCdsujNLMHM5plZupmtNrP7fPfPMrM035/NZpZWxvaXmdlaM9tgZhP8vQMioSb7SAEjJi4idet+nh3Wgxt7tfA6kvjR7QOSKCxyvLEgeCZQleeIvhB4yDnXEegH3G1mnZxzQ51z3Z1z3YF3gDknbmhmEcALwOVAJ+BmM+vkv/gioWXf4aPc8upCVmdk8+LwnlzVrZnXkcTPEhvEcUmnxkxbtIUjR4NjAtUpi945t9M5t9T3dQ6QDjQ//rgVnx4wBJhZyuZ9gA3OuY3OuaPAm8A1/gguEmoyc/K5+ZWFbNhziFduS+bSs5t4HUkCZPS5rTmQW8A7S7d7HQU4zTF6M0sCegCLStx9LrDbObe+lE2aA9tK3N5Oif9JnPDaY8wsxcxSMjOD7zxUkTOxM/sIQ19ewNZ9ubw2qjfnt2/kdSQJoOTEenRtUYfJQXIFqnIXvZnVpHiI5n7nXMlzh26m9KN5gNJOBi51r51zrzjnkp1zyfHx8eWNJRL0tu3LZcjLC9iTk88bd/bhnDYNvY4kAWZm3DmwFRuzDjNv7R6v45Sv6M0siuKSn+6cm1Pi/kjgemBWGZtuBxJK3G4BZFQsqkjoWbJ5Hzf+8weycwuYNrovvZPqex1JKsngLk1pWieGiUFwBarynHVjwCQg3Tn39AkPXwSscc6VNRC1BGhrZq3MLBoYBnxwJoFFQkFRkePFf21g2CsLqREVwayx/emeUNfrWFKJoiKqMeqcJBZs3MuqHdmeZinPEf0A4FbgghKnUw72PTaME4ZtzKyZmX0C4JwrBMYDn1P8Ie5s59xqv6UXCUJ7D+Vz++tLePyztVzeuQkf3jOQjk1rex1LPDCsT/EEqskeT6A65YUnnXPzKX2sHefcqFLuywAGl7j9CfBJxSOKhI7Fm/Zx78xl7Ms9yiPXdeaWPi21bk0VVqdG8QSqaQu38PBlHWhSJ8aTHJoZK+IHRUWOF+Zt4OZXF1IjOoJ3x53D8L6JKnnhjgGtOOYcbyzY7FkGFb3IGdp7KJ9Rry/hic/XMrhLUz4YP4Czm9XxOpYEiZYNYrm0UxOmL9pK7tFCTzKo6EXOwKKNexn83Hcs3LiXv17XheeGdadWTJTXsSTIjD63FdlHCngn1ZsJVCp6kQooOVQTGx3Je+MGcEtfjcdL6Xol1qNbQl0meTSBSkUvcpqyDuUz8rXFPPH5Wq7s2owP7xlIp2Y6q0bKZmaMHtiKzXtz+WpN5U+gOuVZNyLybws37uXemcs4cKSAR6/vwrDeCTqKl3K5vHMTmtetwcTvNnJxp8aV+t46ohcph6Iix/Nfr+eWVxdSs3rxUM3NOnVSTkOkbwLVok37WLm9cidQqehFTuH4UM2TX6zjyq7N+EBDNVJBQ/skEBcdwaRKvq6sil7kJBZu3MvgZ79j8aZ9PHp9F54d1p2a1TXiKRVTOyaKob1b8tGKnezMPlJp76uiFynFsSLHP74qMVRzt4ZqxD9uH5BEkXNM+aHyrkClohc5QdahfEa9tpin5q7jqm7FQzVaq0b8JaF+LJd1bsKMRVs4nF85E6hU9CIlLPjp30M1j13fhWeGaqhG/O/Oga05mFfI25U0gUpFL0LxUM1zX61n+MSF1IwpHqoZpqEaCZBeifXo0bIur32/iWOVMIFKRS9VXmZOPiMnL+bpueu4ulszPhyvoRoJvNEDWxdPoErfHfD30u+kUqX98FMW972ZxsEjBfzthi4MSdYEKKkcl57duHgC1fxNXBLgC8XriF6qpGNFjme/XM+IiYuoHRPJ++MHMLS3hmqk8kRGVOP2AUks3rSPFdsPBPS9VPRS5WTm5HPb5EX8/ct1XNO9OR+MH0iHJhqqkco3tHcCNatHMinAV6BS0UuVkpNXwIiJi0jZvJ/Hb+jK00O6EaezasQjtWKiGNY7gY9X7CTjQOAmUKnopco4VuS4d+YyNmQeYuLIZIZoQTIJAqOOT6BasDlg76GilyrjkY/Tmbc2k/939dmc2zbe6zgiALSoF8vlXZoyY9HWgE2gUtFLlTBt4RYmf7+J2wckMaJfotdxRP7D6IGtyMkr5K2UbQF5fRW9hL3567P4wwerOb99PL+/opPXcUT+S4+W9eiVWI/J328OyAQqfQolYW3DnkP8YnoqbeJr8tzNPYiopjF5CU73X9SWrftyKXKOCPz7c6qil7C1//BR7pyyhOiIakwcmayLdktQC+TnRip6CUtHC4sYOy2Vndl5zLyrLwn1Y72OJOIZjdFL2HHO8bt3V7J40z6euLErvRLrex1JxFMqegk7L3+7kbdSt3PvhW25pntzr+OIeE5FL2Hls1W7+Ntna7iya1MeuKit13FEgoKKXsLGqh3ZPDArja4t6vLkTd0061XER0UvYWFXdh53TllCvdgoXr2tFzFREV5HEgkaOutGQl7u0UJGv7GEQ3mFvP2Lc2hUK8brSCJBRUUvIa2oyPHgrOWszjjIxNuSdWUokVJo6EZC2lNz1/LZ6l38bnBHLuzY2Os4IkFJRS8h653U7bww7ydu7pPAnQNbeR1HJGip6CUkLd60jwlzVnDOWQ340zWddYaNyEmo6CXkbNl7mLFTU0ioF8tLw3sRFaEfY5GTOeW/EDNLMLN5ZpZuZqvN7L4Sj91jZmt99z9exvabzWylmaWZWYo/w0vVk32kgDunpFDkYNKo3tSJ1UJlIqdSnrNuCoGHnHNLzawWkGpmc4HGwDVAV+dcvpk1OslrnO+cy/JDXqnCCo8VMX7GUjZnHWbqnX1p1TDO60giIeGURe+c2wns9H2dY2bpQHPgLuAx51y+77E9gQwqVZtzjj9+uJrv1mfxtxu60P+sBl5HEgkZpzW4aWZJQA9gEdAOONfMFpnZN2bWu4zNHPCFmaWa2ZiTvPYYM0sxs5TMzMzTiSVVwJQfNjNt4VbGDmrN0N4tvY4jElLKPWHKzGoC7wD3O+cOmlkkUA/oB/QGZptZa+fcidfBGuCcy/AN7cw1szXOuW9PfH3n3CvAKwDJycn+v5aWhKx5a/fwp49+5KKOjXn4sg5exxEJOeU6ojezKIpLfrpzbo7v7u3AHFdsMVAENDxxW+dchu/vPcC7QB9/BJeqYe2uHO6ZsYz2TWrz7LDuuhSgSAWU56wbAyYB6c65p0s89B5wge857YBoIOuEbeN8H+BiZnHAJcAq/0SXcJd1KJ87Xl9CjegIJo1MJq66VuwQqYjy/MsZANwKrDSzNN99vwUmA5PNbBVwFBjpnHNm1gyY6JwbTPGZOe/6JrNEAjOcc5/5eyck/OQVHGPs1FSyDuUze2x/mtWt4XUkkZBVnrNu5kOZlyQfUcrzM4DBvq83At3OJKBUPc45JryzgtQt+3nhlp50S6jrdSSRkKYphRJ0nv96A++lZfDLS9pxRdemXscRCXkqegkqH63I4Km567iuR3PuPr+N13FEwoKKXoJG2rYDPDR7OcmJ9Xjshi5aqEzET1T0EhR2HDjC6CkpNKpdnZdv7UX1SF0KUMRfdL6aeO5wfiGjp6SQX3CMGXf1pUHN6l5HEgkrKnrxlHOO3727krW7DjJ5VG/aNa7ldSSRsKOhG/HUrCXbeC8tg/subMd57U+2AKqIVJSKXjyTvvMgf/hgNQPbNGT8BTrDRiRQVPTiiUP5hdw9fSm1a0Tx96Faw0YkkFT0Uumcc/x2zko27z3Mc8N6EF9LH76KBJKKXirdzMXb+GB5Bg9c1E4XEBGpBCp6qVQ/Zhzkjx+u5ty2DTXzVaSSqOil0uTkFXD3jKXUiy0el6+mcXmRSqHz6KVSOOf4zZyVbNl7mJl39aOhJkWJVBod0UulmL5oKx+t2MlDl7Snb2uNy4tUJhW9BNyqHdn86aMfGdQunl/87Cyv44hUOSp6CaicvALGHx+XH9JN4/IiHtAYvQSMc44Jc1aybf8RZt7VT4uViXhER/QSMNMWbuHjFTt56JJ29GlV3+s4IlWWil4CYtWObP78UTrntY/n54M0Li/iJRW9+N3BvALGTV9K/bhonh6i8+VFvKYxevEr5xwT3lnBjgNHmDWmH/Xjor2OJFLl6Yhe/Grqwi18snIXv7q0PclJGpcXCQYqevGblduz+ctH6ZzfPp4x57b2Oo6I+KjoxS+yjxQwbkYqDWpG85TG5UWCisbo5Yw55/j12yvIOJDH7LEalxcJNjqilzM25YfNfLZ6Fw9f2p5eiRqXFwk2Kno5I8u3HeCRT9K5sEMj7tK4vEhQUtFLhWUfKV5fPr5mdZ68SevYiAQrjdFLhTjnePjt5ezKzmPW2P7U07i8SNDSEb1UyGvfb+bz1bv59WUd6JVYz+s4InISKno5bWnbDvDop+lc1LERo89t5XUcETkFFb2cluzcAu6evpRGtWJ48qZumGlcXiTYaYxeys05xy/fXs7ug3nM/nl/6sZqXF4kFOiIXspt0vxNzP1xNxMu70DPlhqXFwkVKnopl2Vb9/PYp2u4uFNj7hyocXmRUKKil1M6kHuU8TOW0bh2DE/eqHF5kVBzyqI3swQzm2dm6Wa22szuK/HYPWa21nf/42Vsf5nvORvMbII/w0vgOef45VvL2ZOTxwvDe1InNsrrSCJymsrzYWwh8JBzbqmZ1QJSzWwu0Bi4BujqnMs3s0YnbmhmEcALwMXAdmCJmX3gnPvRf7sggTTxu018mb6H/72yE90T6nodR0Qq4JRH9M65nc65pb6vc4B0oDnwC+Ax51y+77E9pWzeB9jgnNvonDsKvEnx/xwkBKRu2c/fPlvDpWc35vYBSV7HEZEKOq0xejNLAnoAi4B2wLlmtsjMvjGz3qVs0hzYVuL2dt99pb32GDNLMbOUzMzM04klAXAg9yj3zFhKkzoxPK5xeZGQVu6iN7OawDvA/c65gxQP+9QD+gG/Ambbf7dBae3gSnt959wrzrlk51xyfHx8eWNJABwtLOL+WWlkHsrnhVt6UqeGxuVFQlm5it7Moigu+enOuTm+u7cDc1yxxUAR0PCETbcDCSVutwAyzixy2TIOHOFoYVGgXr5KyCs4xtipKfxrbSZ/uOpsumlcXiTkleesGwMmAenOuadLPPQecIHvOe2AaCDrhM2XAG3NrJWZRQPDgA/8EfxE+w8f5ern5/P791biXKm/NMgpHMovZOTkxfxrXSaPXNeZEf0SvY4kIn5QniP6AcCtwAVmlub7MxiYDLQ2s1UUf8g60jnnzKyZmX0C4JwrBMYDn1P8Ie5s59zqQOxIvbhobunTktkp23npm58C8RZh7UDuUYZPXETKlv08M7Q7w/uq5EXCxSlPr3TOzaf0sXaAEaU8PwMYXOL2J8AnFQ14Oh64uB2b9+by+GdraVk/liu7NquMtw15e3LyuHXiYjZlHeal4T255OwmXkcSET8Kq0XNzIzHb+xKxoEjPDh7OU3r1NBa6aewfX8uIyYuYvfBfCaP6s3Atid+zCIioS7slkCIiYrglduSaVonhrveSGHr3lyvIwWtjZmHGPLPBew9fJRpo/uo5EXCVNgVPUD9uGheG9WbY0WOUa8vJju3wOtIQSd950GGvLyA/MIi3hzTj16J9b2OJCIBEpZFD9A6viav3NqLbftyGTstRaddlrBs636GvryAqIhqzBrbn7Ob1fE6kogEUNgWPUDf1g14/MauLNy4j9/M0WmXAD/8lMXwiYuoFxfN7LH9adOopteRRCTAwurD2NJc16MFW/bm8syX62nVMJbxF7T1OpJnvkrfzS+mLyWpQSzT7uxLo9oxXkcSkUoQ9kUPcN+FbdmyN5cnv1hHQv1Yrule6nI7Ye3D5Rk8MCuNTs1qM+X2PtSL02UARaqKKlH0ZsZjN3Rhx/4j/OqtFTSrW4PeSVXnw8c3F2/lN++upHdSfSaNTKZWjNauEalKwnqMvqTqkRG8fGsvmterwZg3UticddjrSJVi4ncbmTBnJYPaxjPl9j4qeZEqqMoUPRQvk/DaqOLVlG9/fQn7Dx/1OFHgOOd45st1/OXjdC7v3IRXb0umRnSE17FExANVqugBkhrG8cptyezYf4Sx01LJLzzmdSS/c87xyMfpPPPlem7s1YJ/3NyD6Mgq960WEZ8q+a+/d1J9nripK4s37eM374TXaZfHihy/mbOSifM3MeqcJB6/oSuREVXy2ywiPlXiw9jSXNO9OVv35vLU3HUkNojjvotC/7TLgmNFPDArjY9W7GT8+W146JJ2ujKUiFTdogcYf0EbNu/N5e9frqNlgxpc16OF15EqLK/gGHdPX8pXa/Yw4fIO/PxnZ3kdSUSCRJUuejPj0eu7sONALr9+eyXN6tSgb+sGXsc6bYfyC7lrSgoLN+3lz9d25lZdMERESqjyg7fRkdX454hetKhfg7HTUtmYecjrSKflQO5RRkxcxOLN+3h6SDeVvIj8lypf9AB1Y4tPu6xmxtQbw2UAAAkYSURBVB2vL2FfiJx2mZmTz7BXFvJjxkFeHN4zpIeeRCRwVPQ+iQ3iePW2XmRk5zF2akrQn3a548ARhry8gC17c5k0KplLdVUoESmDir6EXon1eeqmbizZvJ+H314RtKddbso6zJB/LiArJ5+pd/bh3LbxXkcSkSBWpT+MLc1V3ZqxdV8uT3y+lsQGcTx4cTuvI/2HNbsOMmLiYoqcY+aYfnRurrXkReTkVPSlGHfeWWzZe5jnvlpPYv1YbugVHGPfadsOMHLyYmpERTBtdF/aNKrldSQRCQEq+lKYGX+5tgvb9x9hwpzi1S77n+XdaZe7svP4as1u/vpxOg1qVmf66L4k1I/1LI+IhBYLxnHo5ORkl5KS4nUMso8UcMNLP7DnYB5zxg2otKsx5RUcY8nmfXy7LpNv1mWybnfxKZ8dm9bm9dt701gXDBGRE5hZqnMuudTHVPQnt21fLte+8D1x1SN5d9w5NKhZ3e/v4Zzjp8xDfLMui2/XZbJo017yCoqIjqhGn1b1GdSuIYPaxdO+cS0taSAipTpZ0Wvo5hQS6sfy6shkbn5lIWOmpjJ9dF9ios58ud/sIwX8sCGLb9dn8u26LHYcOAJA6/g4hvVuyc/axdO3dX1io/UtEpEzoxYph54t6/H3od0ZN30pv3xrOc8N60G1aqd3ZH2syLFyR/b/DcekbTvAsSJHreqRnNOmAePOP4tBbeM19i4ifqeiL6fBXZoy4fIOPPbpGpIaxPHLS9ufcpvdB/P4Zl0m367LZP6GLA7kFmAGXZrXYdx5ZzGoXTzdE+oSpWWERSSAVPSnYeyg1mzOOszz8zbQskEsQ5IT/uPxvIJjpGzez7frM/lmbSZrd+cAEF+rOhd2aMygdg05t2089XVhbhGpRCr602Bm/Pnazmzff4TfzllJi7o1aFwnhm/WZvLt+kwWbvz3h6jJSfWY0LMDg9rG07GpPkQVEe/orJsKOJhXwA0v/sCGzEMc/8/XqmEcg9oWnx3Tr3UD4qrr/6EiUnl01o2f1Y6J4rXbe/P81xvo3LwOP2unD1FFJHip6CuoRb1YHruhq9cxREROSad7iIiEORW9iEiYU9GLiIQ5Fb2ISJhT0YuIhLlTFr2ZJZjZPDNLN7PVZnaf7/4/mtkOM0vz/RlcxvabzWyl7znBe3K8iEiYKs/plYXAQ865pWZWC0g1s7m+x/7unHuyHK9xvnMuq8IpRUSkwk5Z9M65ncBO39c5ZpYONA90MBER8Y/TWgLBzJKAb4HOwIPAKOAgkELxUf/+UrbZBOwHHPCyc+6VMl57DDDGd7M9sLbcwbzTEAjn31TCef+0b6ErnPfvTPYt0TkXX9oD5S56M6sJfAM84pybY2aNfYEc8GegqXPujlK2a+acyzCzRsBc4B7n3LcV3JGgYmYpZa0tEQ7Cef+0b6ErnPcvUPtWrrNuzCwKeAeY7pybA+Cc2+2cO+acKwJeBfqUtq1zLsP39x7g3bKeJyIigVGes24MmASkO+eeLnF/0xJPuw5YVcq2cb4PcDGzOOCS0p4nIiKBU56zbgYAtwIrzSzNd99vgZvNrDvFQzebgbFQPFQDTHTODQYaA+/61mKPBGY45z7z6x54q9TPG8JIOO+f9i10hfP+BWTfgnI9ehER8R/NjBURCXMqehGRMKeir4CyloUIJ2YWYWbLzOwjr7P4m5nVNbO3zWyN73vY3+tM/mJmD/h+JleZ2Uwzi/E605kws8lmtsfMVpW4r76ZzTWz9b6/63mZsaLK2LcnfD+XK8zsXTOr64/3UtFXzPFlIToC/YC7zayTx5n87T4g3esQAfIs8JlzrgPQjTDZTzNrDtwLJDvnOgMRwDBvU52x14HLTrhvAvCVc64t8JXvdih6nf/et7lAZ+dcV2Ad8Bt/vJGKvgKcczudc0t9X+dQXBRhsyyEmbUArgAmep3F38ysNjCI4lOGcc4ddc4d8DaVX0UCNcwsEogFMjzOc0Z8kyv3nXD3NcAU39dTgGsrNZSflLZvzrkvnHOFvpsLgRb+eC8V/RnyLQvRA1jkbRK/egZ4GCjyOkgAtAYygdd8Q1MTfXM8Qp5zbgfwJLCV4vWpsp1zX3ibKiAa+9bgOr4WVyOP8wTKHcCn/nghFf0Z8C0L8Q5wv3PuoNd5/MHMrgT2OOdSvc4SIJFAT+Al51wP4DCh+6v/f/CNVV8DtAKaAXFmNsLbVFIRZvY7ioeIp/vj9VT0FVTashBhYgBwtZltBt4ELjCzad5G8qvtwHbn3PHfwN6muPjDwUXAJudcpnOuAJgDnONxpkDYfXxmvu/vPR7n8SszGwlcCQx3fpropKKvgLKWhQgHzrnfOOdaOOeSKP4g72vnXNgcFTrndgHbzKy9764LgR89jORPW4F+Zhbr+xm9kDD5oPkEHwAjfV+PBN73MItfmdllwK+Bq51zuf56XRV9xRxfFuKCU11hS4LSPcB0M1sBdAf+6nEev/D9lvI2sBRYSfG/75BeLsDMZgILgPZmtt3M7gQeAy42s/XAxb7bIaeMfXseqAXM9fXKP/3yXloCQUQkvOmIXkQkzKnoRUTCnIpeRCTMqehFRMKcil5EJMyp6EX8wLci5rgSt88Lx5U/JTSp6EX8oy4w7pTPEvGAil6qHDNL8q35PdG3bvt0M7vIzL73rXHex7fm+Xu+dcEXmllX37Z/9K0j/i8z22hm9/pe9jHgLN8klyd899Usse79dN9sVZFKV56Lg4uEozbATcAYYAlwCzAQuBr4LbANWOacu9bMLgDeoHgWLUAH4HyKZzCuNbOXKF4YrbNzrjsUD91QvKrp2RQvFfw9xTOq51fGzomUpCN6qao2OedWOueKgNUUX8jCUbx0QBLFpT8VwDn3NdDAzOr4tv3YOZfvnMuieEGtxmW8x2Ln3Hbfe6T5Xlek0qnoparKL/F1UYnbRRT/plvaMMvx9UJKbnuMsn8zLu/zRAJKRS9Sum+B4fB/wzBZp7jmQA7FQzkiQUdHGCKl+yPFV6FaAeTy72VxS+Wc2+v7MHcVxVcF+jjwEUXKR6tXioiEOQ3diIiEORW9iEiYU9GLiIQ5Fb2ISJhT0YuIhDkVvYhImFPRi4iEuf8PVDCskScmyOYAAAAASUVORK5CYII=\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": 11,
"metadata": {},
"outputs": [],
"source": [
"historical = area_weighted.sel(time=slice('1979', '2013'))"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"6. Generate anomalies"
]
},
{
"cell_type": "code",
"execution_count": 12,
"metadata": {},
"outputs": [],
"source": [
"ssta = historical.groupby('time.month') - climatology"
]
},
{
"cell_type": "code",
"execution_count": 13,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"[<matplotlib.lines.Line2D at 0x2b339ce204e0>]"
]
},
"execution_count": 13,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAESCAYAAADwnNLKAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nOydeZwkZX3/P09VV1/TPffsvcseLKcs13IoCIKoKCreShKvaJBfJBpfJhETTdQkhmhMDGokRFE0KsYTVJRLBJZDWJCbXViWhZ09Z2dn5+yjjuf3R9W3+qnqqu6q7p7pnp3n/Xrta3uqq6qfrq56vs/3ZpxzSCQSiWThorR7ABKJRCJpL1IQSCQSyQJHCgKJRCJZ4EhBIJFIJAscKQgkEolkgSMFgUQikSxwEu0eQCMMDg7y1atXt3sYEolEMq946KGHDnDOh/zb56UgWL16NTZv3tzuYUgkEsm8gjH2QtD2lpiGGGPXMsb2M8aeCHn/jxljjzn/7mWMnSi8t4Mx9jhj7BHGmJzdJRKJZI5plY/g2wAurPH+8wDO5ZxvAPCPAK7xvX8e5/wkzvnGFo1HIpFIJBFpiWmIc34XY2x1jffvFf68H8CKVnyuRCKRSJqnHVFDHwDwa+FvDuAWxthDjLFL2zAeiUQiWdDMqbOYMXYebEFwtrD5LM75bsbYIgC3Msa2cM7vCjj2UgCXAsCqVavmZLwSiUSyEJgzjYAxtgHANwBczDkfpe2c893O//sB/AzA6UHHc86v4Zxv5JxvHBqqin6SSCQSSYPMiSBgjK0C8FMA7+acPyNs72KM5ek1gFcDCIw8kkgkEsns0Krw0R8AuA/A0YyxYcbYBxhjlzHGLnN2+XsAAwD+yxcmuhjAJsbYowAeAPArzvlvWjEmiUQiaSUjkyWc+fnb8dTuiXYPpeW0KmrokjrvfxDABwO2bwdwYvUREolE0lnct30UeyeK+M/bn8F/v/vwinSXtYYkEokkAj0ZDQCwd6LU5pG0HikIJBKJJAJF3QQA7BsvtnkkrUcKAolEIolAoWwLgr0TUhBIJBLJgqTgaAQAYFq8jSNpPVIQSCQSSQRmyhVBMF022jiS1iMFgUQikUSgKGgE0yUpCCQSiWTBMSNoAVIQSCQSyQKkULbc11Mls8ae8w8pCCQSiSQCBV1qBBKJRLKgKQjO4ikpCCQSiWThMVM2kUoozmspCCQSiWTBUdBNDOVTAKSPQCKRSBYkhbKJwZwtCKSPQCKRSBYgBd1Ef1cSCpOCQCKRSBYkhbKJbFJFVzIhncUSiUSyECnoJjKaiq5UQmoEQTDGrmWM7WeMBbaZZDZXMca2McYeY4ydIrx3IWNsq/PeFa0Yj0QikbSambKJTFJFV0rFtHQWB/JtABfWeP+1ANY7/y4F8HUAYIypAL7mvH8cgEsYY8e1aEwSiUTSMkqGibSmIpeSpqFAOOd3AThYY5eLAXyH29wPoJcxthTA6QC2cc63c87LAK539pVIJJKOgXOOsmEhqSrozmg4VNDbPaSWMlc+guUAdgp/DzvbwrZXwRi7lDG2mTG2eWRkZNYGKpFIJH4Mi8PiQCqhYHlvBrvGCu0eUkuZK0HAArbxGturN3J+Ded8I+d849DQUEsHJ5FIJLUoG3bBuWRCwYq+DA5MlTxlqec7cyUIhgGsFP5eAWB3je0SiUTSMXgFQRYAMHwYaQVzJQhuBPAeJ3roTADjnPM9AB4EsJ4xtoYxlgTwLmdfiUQi6RhKjiBIJVSs6MsAAIbHZto5pJaSaMVJGGM/APAKAIOMsWEA/wBAAwDO+dUAbgLwOgDbAMwAeL/znsEYuxzAzQBUANdyzp9sxZgkEomkVRzuGkFLBAHn/JI673MAHw557ybYgkIikUg6krJp+wOSCQWDuSQAYHSq3M4htRSZWSyRSCR1INNQUlWQUBUkVQUF6SyWSCSShYPrI9DsKTOtKTJqSCKRSBYS5CNIqfaUmUmqno5l8x0pCCQSiaQOorMYALLJhDQNSSQSyUJCDB8FgLSmSkEgkUgkCwm/RpCRPgKJRCJZWIjho4D0EUgkEsmCo6STaYg0AmkakkgkkgVF2fSahqSPQCKRSBYY1T4CFUVpGpJIJJKFg5hZDDg+AqkRSCQSycKhEj5a0QhmpEYgkUgkCwdqU8mY3UsrrakoGRYsK7CP1rxDCgKJRCKpQ9mwXP8AYJuGAKBoHB5agRQEEolEUoeSYbpmIcA2DQFwcwm+fc/z+PxNT7dlbK2gJYKAMXYhY2wrY2wbY+yKgPf/mjH2iPPvCcaYyRjrd97bwRh73HlvcyvGI5FIJK2kSiMgQaCbeGLXOD7zi6dwzV3b2zW8pmm6MQ1jTAXwNQCvgt2D+EHG2I2c86doH875FwF80dn/DQA+xjk/KJzmPM75gWbHIpFIJLNB2QwxDekmHnphrF3Dahmt0AhOB7CNc76dc14GcD2Ai2vsfwmAH7TgcyUSyWHA5h0HMV0y2j2MmpCzmKiYhiyMTtudyrKOcJiPtEIQLAewU/h72NlWBWMsC+BCAD8RNnMAtzDGHmKMXdqC8UgkknnC3vEi3nb1ffi7nz3e7qHUJMxZXNBNHJwuAQAMc/5GELWiZzEL2BZ2Rd4A4B6fWegszvluxtgiALcyxrZwzu+q+hBbSFwKAKtWrWp2zBKJpAPYPW43gN82MtXmkdTGbxpKa6IgKLv7WBaHogRNiZ1NKzSCYQArhb9XANgdsu+74DMLcc53O//vB/Az2KamKjjn13DON3LONw4NDTU9aIlE0n4OTNqr6d5Mss0jqU0p1DRUEQRApSbRfKMVguBBAOsZY2sYY0nYk/2N/p0YYz0AzgVwg7CtizGWp9cAXg3giRaMSSKRzAP2jBcBAD0Zrc0jqU1oHoHuFQRUpXS+0bRpiHNuMMYuB3AzABXAtZzzJxljlznvX+3s+mYAt3DOp4XDFwP4mZOtlwDwfc75b5odk0QimR+QIMh0uKO1bFjBeQR+QWCYADpbqAXRCh8BOOc3AbjJt+1q39/fBvBt37btAE5sxRgkEsn8Y4/jI+j0Am5V4aOOIJguGRib0bGkO429E0W3JtF8Q2YWSySStrHnkK0RdHq3L3/4aDppv947XoRpcSztTQMgjWD+IQWBRCJpGzvHZgAAM+V5kEcgaARJVYHCgBcP2uNf0ZcFABTnqY9ACgKJRNIWirqJvROORtDhE6jfNMQYQ0ZT8cKoLQhW9WcAQJqGJBLJwuaov/s1Pnjdg5H3Hx6bAXcyjgrzQSNQvQ7tTFLFjlE79mVVv60RSNOQRCJZ0JRNC7c9vT/y/mRWWTPY1fFNXvymIaDSkwAAVvaRIJAagUQikUSGzCrHLMl3tLOYc15lGgIqkUP5dAK9WTshbr7mEUhBIJF0AE/tngDn87dWTSPsnyxBUxlW9GU6WiOgbOGk6i0dQUXmFnenkdLsqVSahiQSSUPc/ewIXnfV3fj+Ay+2eygNYzbQsrGkW0gnVGSSCRR0s2PbPpapcX2AaQgAFnen3GQzaRqSSCQNMTxmJ1U9tnO8zSNpnGIDCWElw0RKU9yVdae2fXQFgeozDTnjXtaTQSphv5aCQCKRNERlNdmZE2EUGskMLuoWUgnVFQSdah5yTUMJb9QQCYgjF+UqpqEOz5AOQwoCiQfDtPDd+1+APk+rKM5H5vtqEoifGfy1O7bhzmdGkNKUqv6/nUaYaWjXIVuTO3JRTpqGJIcX3/v9i/j0z5/AdffuaPdQFgzJeT6JAPFMQ5xzfPHmrTgwVUIqoaIrZZc8myjqszW8pggTBLsdQbBuKOeajebrbygFgcTDoRn7YRwvdOZDeThCjtb5bBqKY9YR761UQsHaoS4AwLb9ndmcphTiI3jJ8h4AwMr+LBhjSCaUlvyGJcOc83uhJdVHJRJJ4xiWPdHM1zo1QDwfwYjTjAawBcGRQ7Zp5Yld47j4pMAut22FfAQpn0bwrfedhj3jRahOR7JUQmlJHsGJn70F3WkND/zdBU2fKypSEEgCmX/N9uYv5I+ZzxpBVEGwf6KIT/280nsqpalIqAqOWZLHE7smZmt4TRFmGurNJt1EMsD29bTiNyzqFop6qf6OLUSahiQeeGi7aclsoTtNz8sdYF8uGWZDiW3FiKahK3+9Bb9/vtKynFbZRy/J47kO7VscJgj8ZJLKvNXqWiIIGGMXMsa2Msa2McauCHj/FYyxccbYI86/v496rERyuGOY5CNo7yQyPqPj6E/9Bv991/bYx0bVCMKSsnqzSUwWO7PwXFgegZ+MprY08mkuE+yaFgSMMRXA1wC8FsBxAC5hjB0XsOvdnPOTnH+fi3msRHLYQj6CdtepGZuxWy5+974XYh8bVRB0+3oTk0aQS9nZxUYHhi1X8ggiCIIW5hEcmsOAjVZoBKcD2MY53845LwO4HsDFc3CsRHJYQKahdmfWGs4KtJEmMVFXwhO+yU0UBAAwVYr32R/+3sP4qx89GuuYuEQ1DaVbIAhEs5zoVJ9tWiEIlgPYKfw97Gzz81LG2KOMsV8zxo6PeaxE0nFYFm/Jw+o6i9usEVAuwHSpkSzhyjG1fAwUOrqiz27kQsl0+bQtCMg8tPqKX+Ffbnq67uf+6vE9+PFDw7HHG4fIpqGk2lCpDRFaFADAgan5JQiCAkz8d8LDAI7gnJ8I4CsAfh7jWHtHxi5ljG1mjG0eGRlpeLCS2lBMe4fW/+oovnXvDpz2z7dhx4Hpps5D5pCCbuK/fretFUNrCFrNlhswz4gr4VoF6CaKOk49og9vO3UFAEBzKnr6BQGAhnwVs0EpJHzUD/kIyoaF1Vf8qqHfUryO800jGAawUvh7BYDd4g6c8wnO+ZTz+iYAGmNsMMqxwjmu4Zxv5JxvHBoaasGwJUHQ6qeRyWCh8djwIQDAvc+NNnUeWgWesLwHd2yJ3til1TTj6BQTysRVrZ/xgo6ejOY6iS1He8inbd/BVMnouPIm0aOGbNMQmdb+87ZnY3+W+Bv88MGdc+YzaYUgeBDAesbYGsZYEsC7ANwo7sAYW8IYY87r053PHY1y7OHMcyNTeONXN+GebQfaPRQXilzphFDGTucIpz3hU3uaqxpqWBZUhWHdUJfbw7cdNGPfPjBVdl/rVvi9M1Ew0J1OIOEkYdE8Rz6CyaIe2bwyV/doZEGg2aYhEoqNRIHRb7B2sAv3bR/Flr2Tsc/RCE0LAs65AeByADcDeBrA/3HOn2SMXcYYu8zZ7W0AnmCMPQrgKgDv4jaBxzY7pvnCs/sm8djwOP74G79v91Bc6Ob99r078MHrNrd5NJ0NXavHhpsUBCZHQmFY3J3GvolS2xrUiBNw3JXo8NiM+1qvMQGSRqC6gsDel0xDUyWj7gT61d8+i4u/ugmPOhrZbBM3fLQZgUoawVlHDgKYu4qsLcksdsw9N/m2XS28/iqAr0Y9dqHQYRowAO8q67an97VxJJ0PPaTN1mXSTQ5NVbC4O42yYeHQjI6+rmT9A1uMKAjGZnQM5VORj6WeCkAl+siPZXFMFHV0i4LAEXo5RxD84tE9OH5ZT83PuvWpfXh0eBz/dcfc+FPKpgmFAYkIzuIZ3WzKxFbQbbPSYM6+9tMxo6gaRWYWtxGjhgrdLqRvIDrTji241go4CrppQVMZlvSkAaBt5iFxAjs0U66xp5eibmJksoRVjqlseGwG2wOyhKfKBjiHTyOw3+t2fAS3Pb0PX7pla83PIzkjCp+4/OLR3XhxdKb+jghuXB9EWlPBeSUfoxEKZfuCDObthUDccNpGkYKgjRiCU61T+tXO18Ya7WCmRFE2zf12hmUhoSpY3G2vAn/44M46R8wOBSF8NY55g8oxrx60q4i+9ev34fwv3Vm1H028fdmk4COojsihCrhh0OQoCoI4zw/nHH/xgz/gdVfdHWn/smHVNQsBlWb2B6ejCwLT4vj0z59wI8/oupNG0EhORyNIQdBGxDC7Rnq+zgZSI4jOjPPQNhvlopscmsKwpMeOrf/2vTuaDkltBHHyj2PeoAYtqweynu3+6/Ljh4aRVBWcd8winLl2AADwJqfaqBNLAsDrlA3yVVCIqTjeWpFKfui4qKvtsmlVdScLglpXioLgvjoRZS+MTuO797+A933rAQCViZ8EwVQDOR2NIAVBGxFtqWF21blGRgtFZ8aZSJoVBIZpawTLezO46ISlANrTD0LUBuNoBBMF+zos8vkUSFMg7npmBOccNYj+riSOGOjCjisvwsscp6jI6HQlfj6o/tBUqfraxMnKnopZ06hkWHVzCIBgjeCS/7m/5jEUPrvD0ZaKrkZgm4ZmpGno8McUfASdshIPiti4dtPzuPtZmcTnZ7rcIo3A4kg4iVV/fOYqAO3p3ytO/nEyZMumvS91GiNe8Nngx2bKrh8kiIc//SrkUgkcmKxMpH5BoJtWYIXPOFnZkzEn1zg+AgAYjWEa8n8X0sTyaQ3JhIIpaRo6/PFoBE3amVtFkEbwuV8+hXd/84E2jKazKZCz2ORN+Xh0wQadTdqTKUWPzCWiOSiORkD3TFfSJwgOVgQB5xwTRQM9vqJzIv1dSSzKpzylFfztK8OiaOIIrrgaQVQfQdYxDY3FEATidZ4uGRh3tKuulIpcKiGjhhYC4uTfKVUX/YLAb6LYeXBmzhxYnc50xGzaehiCRkCTSSP1fpqloJvu51P0ShRIi/RrBDsFQTBVMmBavKYgAOwwUnGB5BcEpCF0Jb02+zgNYUTfQBQBYvsI6k+VlAtBYdd/cf6RAGqXkxY///kD09g3WUR/VxKphIpsUnUDEmYbKQjaiHjDd4ppyD+OL968xfP3y79wB97/rQfnckgdyfn/9jtPLZhmzEO6aSGhkEZAE/HcC4KibqLP6bjViEaQTXknZ9FWTguKeoLAr1X4BSIJAr+JKUpDGM45rn/gRc/vFqWeT1TT0IkrenHBsYvcRQFpd7WebfF33jE6jf0TRSzutr9bLpWQ4aMLAdFH0CmmIX/46P/e/yIA23lFk53YYWqhst2J6qHSCM0IAsPkbvE1mjym26B1FXUL/U4iWxxTC2kEOZ9GMCms5qMKAkosI/zaJ02MS50IKxKcUTSCe7aN4oqfPo5P/vRxd9v+qIIggmlIURgu2rDU/ZuER61M6aLw3o4D09g7UXTDiLNJdc58RVIQtBFv1FBnagTEonzao6ZaVnN28cMFMgc0o9EZVrVG0C5ncXcmAYXF00hcjUAw16zoy3gcvRRZRIljYeR9wsQvkChiiFbNpMFEcRaTsBa1nV2H6ielRTUNAcDy3koILR1TKxJPbPH5/IEZ7JsoYYnz3bqkRrAwEHMHykZnTKphq5dcOuGJYDjtn2/Dp294InDfhQRlkTbjIyibHJozaaQSSuyJuFUUyiYymuoWT4tK2bSQUJjbWwAA1i/KeQQBaQT+DmV+qjWCYNMQCeDerH2+KOGjilJd9V70Y4QR1TQEAMudPgsAkHK0iFqLBBr3+kU5bNs/iQNTJSwiQZCUzuIFQSdqBGGCwDAtT0zz6HQZN/whsGL4YY/o2CdzSCPOftPi2DNegGFa0JxJijGGbDLRFtPQTNlAJplwyylHhSZKmmczmop8WvM4eieimoZS0QSB4iSgxdEICsI1TSYUDOZSkcpMxBEEi4VcikgagXOdzz1qCI8Oj4NzeDQCaRrqAL7+u+fw5O7gypLPjUzhmX3NlYgVJ49mVpStomTYTTWOWZL3bD/nqCHoJq9SU49b1j2Xw+sYaJI8YXkP/vo1RwNozEfwn7c9g5f+y28xPFZwo4YA28Ty++0HsWe88Vo6jTBVMpFLJWK3XCw7CVfLejN42boB/O8HT0d3JhGoEfRko2kEyYQCFqAZ7Z8sgTFgwwq7MN1JK3sBRNMIRMdzPpXAqv4MXoygEZQMy13d10MsTBdFEFB01sdedRROWG5/p6Od5y+fTmCioOO3W/bNumYgBUEIRd3Ev/5mC95+9X2B77/yS3fi1f9xV1OfIWoEndCMg2rKv/dlq/HZNx7vbs9oCnTTqorg6JRs6LmGJsl3nrYSPRl7RdqIae+OrXaS3nhB90wg2aSKp/ZM4KKrNrVgtNGZKunIpRowDTkr5rSm4vt/diZOPaIf+bSGyaLu+pHGCzoUBuSStQse0wq/bFjIBgikPYcKGMqlcPFJy/Cjy17qJuBFCXcVHc+5dAKr+rORBEEcHwEAHDGQxcYj+qCpETQCw0RSVdCVSuAXf3E2tvzjhTj1iD4AwEBXEpMlA3/67c141b/fOas+OSkIQqCwstmc7MwOSyij7zyUS3lufE1VUDasKnNFJwivdkCr1IymIpmwV/KNXAuxbIEm2K9p4ohTvKxZDCdjN5fSbNNQHGdxwESZTyegm9w1NU4UdeTTWqCdXuScoyrdBzPJatPInvEilvVmwBjDaav7MdCVQkJhkSZ0Me+jK5nA8r4M9owXasb5A/FMQwBw51+fhx//v5dVNAIz/FoWyiZSWuXclJ0MAAO5iplp93gx0ndslJYIAsbYhYyxrYyxbYyxKwLe/2PG2GPOv3sZYycK7+1gjD3OGHuEMdYxnVBGnOxGv82ylXg0gg7wEZAgWNSdclczgN2Qo2xaVeppJ5iz2gFNTtmk6l6ni792Dx56YSzWecQJQNQI9jm/Q6LOpNlKSNvrSqkNmYb84ZXUepL8BNOO2akey3srztZsUvXY9QFg93gBy3orOQTJhIIjF+Wwde9E3XOLPq7T1/SjN5OExcNDdXcfKuDOZ0Yih4/6oWNqhY+WDNOtUeSn39eTYs/47JUnb1oQMMZUAF8D8FoAxwG4hDF2nG+35wGcyznfAOAfAVzje/88zvlJnPONzY6nVeyfsB/GrlT9qoONYpodqhHkqzUCPVAQtF94zQWjUyX8yTd+j/2T9oNIk2RGEAQA8Nav3xvrvIpQcVM8D2mKS3vD6/K0GooIy6cTdqetGLV7SkZ1dc5uXzN62xEd7Vm654rzccvHzkFG88bRc86x+1DBzSEgjl6Sx9YILR1FjeCS01e5kUcTISUn3vCVTXjvtQ/ENg0R0XwEpkcLEKHCc8TeThYEAE4HsI1zvp1zXgZwPYCLxR045/dyzmm5dD/sJvUdDWkE/kzHViJqAZ0wqZIgGOhKISk4L7UEg25yz4MEdE5ZjNnm+gd3YtO2A/jWPTsAeE1Dmm+lGMekIpbv0NTq1f9QLnqHsGYhId+VSiCtKZ749nqEmYYAURCYnjyDWizvzeCoxfmq6KVDMzqKuoWlvqzio5fksXu8WLdi63TJQD6VwI2Xn4Wjl+TdUNaJkOOoeJxp8YYEQSpS1JAVqhEM+H7/2WxY1ApBsByA2Elj2NkWxgcA/Fr4mwO4hTH2EGPs0haMpyXQpBgmrVuB2WHO4pGpIvqydtVDcYLTVAW6sXBNQ/RA0yRfcE1DiSqTAWkNQWzZO4HnhT4Dog+AEsoA4LyjbTv5XPaooAk7l0ogmVBj3Y9lw6yKqunvsiexF0adhisxBAGRdXwVnHPsPDjjLs4omYwgc9JIjWsP2BrBYD6FDSvsSCO/sKpFUxpBnTyCtBZ8btE0lE8lOl4jCDJkBt7BjLHzYAuCTwibz+KcnwLbtPRhxtg5IcdeyhjbzBjbPDIy+yWR6aaq10i7GU++YXF3IukU0xA1xBBv/IqPwLtKrDdZjE2X8W83b+2YpjuNQosBKmMw45qGFGgJ7+1fq2TBhV++G+f92+/cv6kT13In7JL45ntPw/nHLIpUP6dVkJDPpRLQVBYrU7oU4Ew9YXkPlvdm8IMH7BIlM7rhls+ICpVY+MEDO/HyL9yBO7bsB1CdlEa/T73rNVMyPMKIspzDNAKRZnwE9UxDqZDFZreQXLe4J93xgmAYwErh7xUAqjKNGGMbAHwDwMWcc7dtD+d8t/P/fgA/g21qqoJzfg3nfCPnfOPQ0FDQLi2FQin9zio/zZQWME3uOgw7QSMYL+hupiZpBApr3EfwqRuewFfv2Ia76vQyGJks4czP346ndtd3+LUDmmiuf3Anbnx0t3tPZJKJKtMQ+ZZEtu6drIoAKhsWpkoG/urVR+GeK87HBcctdt9TFIZcKhGr2UqzUI5ILm1rOfE0gmpBoCoMF5+0DPdvPwjdtDBTMiP7CIhMMoGCbuI3T+4FAFx1+7P2GH1O54ogqH29psuGx9TragQBjW78NGIZiJRQZoSbhsSubUt70h1vGnoQwHrG2BrGWBLAuwDcKO7AGFsF4KcA3s05f0bY3sUYy9NrAK8G0BF1C+im8tvFq/drpsYMd28CvQNWzZNFw10l0U2sKgyaqsDi1SWB62kx486KV2W1o19ufnIv9k4Ucd29Oxoc+exCo+cc+MgP/oBP/MQuWpbR1KqV4j7fw8o5x2u+fBfe8d/efBRqDt+b9ToEibSmxIrlbxYSBF2OcIvTqS4sqoYWFWXDwkzZrCodXY+MpmCyqONhJxqLnsW8rwxF2rlX62oEZdNTIbXiI6he7Pk1/YGu4N+pFpFMQ+Vw0xAA3PDhs3D335yHnoxW9fy1kqY9oZxzgzF2OYCbAagAruWcP8kYu8x5/2oAfw9gAMB/OVLOcCKEFgP4mbMtAeD7nPPfNDumVkAPQr1WcSXdBOqkzYdhWpa70ugEx+tk0cD6RU5mp6sRMNf84XfG1dOGaCKr1+YvatZpuwhbHWeTqmfC1FRWZRoih+O2/VOe7RO+mjl+0ppa1yzZSqaKomlIieX/CYuqodpDJcPCTLkR01DC1czzqYTbWaxhjaBkYGVfpShcxUdQPcH6hcNgPr7jPopGMFUyqvo4iJzoZE6nEmqsLmxxaUlIDOf8JgA3+bZdLbz+IIAPBhy3HcCJ/u2dAD38M7rtrGIhq9pmNYJ0B5mGJp2kH0A0DTFXKIimoTWDXXUTXMi0kQiIiBEhQTCbEVrNIAq8D5y9Bt/c9DwAW8CJ0T5DuVSVs9jfrpGg3ztMSKZjZvc2ixg1lEwosUyeVGLCT6UMs4mZciOmocr+S3vTmNxnC1N/YTpXENQxpfkjl1IJFamEEhg+KvZNBirN5OMQJY9goqDXrb8EzL6GKDOLQ6AVEee1J/s4iTd+DFMwDbXZWUytBPNCrRegYhoCbNX8jJxPHzoAACAASURBVDX9uPkvz8HFJy2DafE63Zfs61bPhETmlJk2tGeMgriiWzPY5b5mjLlVQwFgUXe1Q+/Fg9PwUzYs95xh0SjphIKibmF4bAbbR6YC92klU2UDyYRi/1MZdNOKHAgRlnlLwmGqaMCwOLIx7ezi/mKkkH/BQIupsOeUc47RqRJGJksY8q3suzNaoEbg7zs8kGvANFTHWWxaHJMlo25pbmD2FwYLUhAUdRNP7AouJkeIK/RalSCb+XFMiwumofYKgoJuwrS4oBHYK11yFgO2mSyZUHD0kry7rVZGNEXZ1CvTMTxmF1eLEsbXDsTV8WAu5Vn9irbxIwayVRpAkEYwWdTd+8vvbCYokuTsf70D53/pzsYHH4HP3/Q0Nj17wO0FoKkKOPeGr27ZO4GfPDQceHy4IHB6+Dq+omzMLH2PRiDkDqi+jOt6pqEfPzSMU//pNhgWx2pBkANU2K36vhPDfIHqPglRUBRWMwIrakVWwBaqs2kqXJCC4IcP7sTrv7IJ1zoqfhCiFK/VN7QZQWBYlfK97TYN+eu8uz4C52YGbDMZTVy0rZYAI5tmve+2x2kOErep+FyhCwXlhvJJbPrE+fjVR84G4C0DsXqgC7vHC557Ish8NlE0KhpBiCDwR6nUsjM/svMQVl/xq0jZtX5Kholr7tqOJ3dPuLZq0nJELfX1V23Cx3/0aKAGWAr1EdjbqGdD3DwCURAs6Q7Psk4naguCTdsOuK/XVAmCaifsdffuwN/8+DHPdwozDdejluM9atc2wL4fDIvPmi9xQQoCCuX7ws1bQsu7ls1KWFdNjaAJKW1a3DW9tLvWEKnHrvPSue9VxtwHYqZkugKAEqBqTfL0YNaNLnIeiCAVvRMQi4YN5lIYyqdw/DK7ZLC4Ol0z2AXOvc1OglabEwXdXSVqIaYhv82dErOC+OWjdrT277bur/dVMDpVwjV3PeeafUQtLCdoBIBXEyKtzm8y4ZzbPoIAgUah0YcaFATi/ot7wgUBfU7Yill0yB8xkPW8l0upVXPAP9z4JABgZZ+3lEUjJBMKdo0FF7aLJwgc89csaQULUhDQ5FXULfzi0eDmKrppueFvtZpDNKcRcCQcQdBu0xA5zMheSf9/4OVrPBODqxEErBr90IPpNw09s28SP3zQTjQyrUrpirlqyxcXcUXndxrSSvHl6wdds4NoVgjqpTtR1N3rFlUjeGZfuJ+Arq4SYdX68R89is/ftAWPO6bRqQBBQOVFxO/d5zwLu32tHd3vEeQsVkkjcExDMYMBMlplf39ZCZGU07sg6Fl8/7cecPtuA9VlO7oCKpwSFAEWV4CJHJrR8Zsn9+L/Nu+seo80kSjRclEjoxqlM8M0ZpmZsol8KoGyablNyIn3fesBLM6noZsci/JJ7Bkv1qwf05QgMG2NIKGytoeP+k1DaU3FjisvAgDc+tQ+dz96uKlsci2NoCIIvPu88aubUNQtvPWUFW6WrjiGToMmu/e9bHVgqN8Df/dKdKc11xQm+gWCVqmilhDqLPbFlj9Xw2FMPt0o1gvKZqbv5NEIfIECus83MjajY894wQ1pBCpaQ6BpSCMfQQs0ghqmIcYYUongqBrq+aAw4DNvPL7KxBPUF5j6FHz2jcfjvKMXQa0T9RaFFwJMhHE0glSittbTLAtSEBR1O5RNMVjVzfO7rZUs2NWOGllzsmsqfNRuWq6pCspt1ggqpqHqm1IMkSSNgMomR9Fk/PtQdMfodNnjkOxUQVAyLAzmUviM0KxHZFHenqTSmop8OoHhMUEQ6CYUBohK0URRd82OYc7itK+aZy3zpBWjzAntSfOhaI7r8pmGxPt+KJ/Cs/unsPuQNyqKnp+g7FiavA5N25/RTPhovQJ8dlRN+LP40nUDeM9LV1dt7wowDRmmhbedugJvOaV1tTGDkuni+giA2dMIFqRpiOKJ/bG5/otMD0aU5tONYFocCZVBUzpHI+jOVK8NRPMFJZeRcIgSbx4mSPdNFN3PHcylOtZHoJuWpxprLVb0Zd0oKMAWIuIKGrDNMWR2Cao6ClSbhvQIHdCimIZIfSDZMVkSTUNe4ST+bvQs+E1DpC0H1cshQXDQ0Qji5omIwqVexFE6UTu8kjQhP12pRFUNrYIe3iMgLrd//FwA3utMxNMIpCBoOXZyS6Iqe9NfD4Zu3CgO0UYwLDINxavtMhvQqijIjis6NCtRQ45GEOLkFmPQ/T4CirTZN1HClFPnZVlvGlMlY1bb8TVKnA5VK/oy2HXIKwhW9HkdlGXTqmlSAbymod6sVvP+oGsWRzOghUdNZ7EgfGj/XT5BUFMjcLaNuz6Cxk1D6YSCb753I279WGBNSqS14PBKGtfn33xC4HG5pG0iFv0hBT1+8lsY64Zy6O9KBgaljBd0t8VnPerlSjTLghQERd1ERlOqVhEHprzZhLQKqi0ImosaSigMfV1JN5W+XZCJJmiF6i9JLf4ftlIVHXB+bYdWiqJGsLQn7XSLmrts2qjogpO8Hst7MxgeK7iTc1E3kXacmURRN917KsxZTCtA6oIWJnCBirknzqKEJs0g01BQ+00S5jt8eRF0/wdNZv7w0biTKy1KkqqChKrglccuxvrF+cB9wxKudNPCn79iXZVWRtB3ponasjiKuhVpco6KbX6qHptd2yuallQp1yE1gpZBdU9SmuKZyEd9kzF1J/PHAYur1q374sduE7ZGoGDdYFdNZ+BcQA96kHnB6yNwwked/8PCXsWM6yqNwJn89guCgJyBM3WqvbaDuBrBVMlw1f6SYSGlKbjj46/Al95+ortyrZiGgs9LJrjXb1hqlwGvYRoSo+DqQWchQRAUNRRkGqLx7jgw7bn/CzU0gqQrCBorIULnrFWUjUhpalVopW5ansKOQdB3JocxXZdmIoWqP0ML9H/ZZbGjXRO6BrNVb2hBCoKCbiGTVOtqBPQj+R254rx2x5b9DZuHDNNCQmFYtyiHPePFtoZPUpxzUJ/cZIBGkHQ1guAb07OaFK4f59z9nrZpyH69yEn9j9Pha64ox9AIVjix5+QnKOkmUgkVqwe78NZTVzjFwwSNIETAHLOkG99632n4pzedAE1lNTUCt2FOhPuQ5nCa2EXbNeVEBJVGIGFe0E18+oYn3PvFNQ0layeUMRZtQhchDSKKJpEOiBqqjC38eNIIZnzXsFU+AiA4V4E+M6rAkc7iFvP13z2Hp/dMIKOpSPnsin7zDDnP/JMdmVFOWtmLmbKJJ3fXLlcRBvkI1g058ecj4UlDsw096P70fSDYNEQCI6x8hDj5i1rDVMlwr9/odNldkVLkTa2cjXYRRyOg9oJkDikZlmcCpFIBtLio1aD+vGMWIZlQ6vqQ6JrNlA380y+f8iS0hUEmhiAHPfmExEAA22Fub//f+190wyFp4kwlAjQCZ3/TWZXHzc5NJhQkFBZp1ZzWbAEr4o6txqROWj8tSGZDEHSlEoFRXzNls2blUZGohfUaZcEJgn/9zRYAcKKGvBrB6FTJcwOE+QjIKUfNpYPsf1EwLQ5NZVg7lAMAbD/QPvOQaXEoLDiVXly1uHkEdWqth2kEYinrkmG6K9LBvH0tO1EQ6GZwdc0g+pz+AgenyzAc04Q4SaYE01BSVSJNjvXKQtPkdd9zo/jGpudxxU8fC92XO8YhWgCJlTdJW3C1PVGYmxzHLet2/6YVbq1VN8X3A42bWjLOc1qPtM/MCwDFsv13FNMQfR9qOpRuoWmoK5UILJ8yXTYiXxe6jtJZ3GLoBvNEDc2UPX1Cyabp9xHQipZuIlElPzhdxobP3IxNzx5APchHQJNHlJZ5s4XJuadvrogYukc+Ak2pnUcgagqis1gM4yvqJiaLOnKphHutW636/vKx3fifu7Y3dY44piFqYHJwuuzeW6IQSSVUlAzbNBRVy9CcaqB+xgs6bn1qn2sa2u1UPo0yVrqnRwVzKN37QT4C3bSwvDeDL75tA4CKwKbfK2yyrgiCxlKWskkVmQgmpbTmbXQ/UdTxtz+rNBAKP79fENQXHnHJBySt2Z8V3zR09Z3PVYXwtoIFlVAmOrkymlqVjThRMNCT0dwQOVoV+B9C0zkPZWKKdu2ndk9gomjgEz95DPdccX7N8VDUUFCtlHf8930wTAs//fOzYn/PRjAtjhA54HkoSBPQAiJLRILsy4DYe8BOABov2IKAHshWawSXf/8PAIA/O2dtw+cI68AVRHdGg8Lsfs3BgkBBSbc1grAcAj9BJUhuf3ofPv6jR3FoRnebt9M1rxWXTo9AyTDx0AsHcf/2g/jQOWtx7NJuvPHEZc7nBUQNmRYSqu3PAioJbnTvh02cyYQKIPrK109GUyP5CHK+yfbfb3nGLTYX5L8QjwPmwDQUIAimYzTrIfPiC6Mz2D4yjWW9zddBEmmJRsAYu5AxtpUxto0xdkXA+4wxdpXz/mOMsVOiHttKxIk2o6lVsccTBd2TUJVU7frsYzM6Vl/xK9zwyC4Adq9hoGI6Elcihwq2bXjXoULNJDHOuVt0Lih9/IHnD+LhFw81/F3jYpjhGoHoN6j4CGoXnfNoBAGCYHFPGgXdxOPD4zhqSd592Dsxakg3eWhxOD+qwtCbTWJ0uhy4WiYfQVyNwG+C+8B1m13tyh/bH6W+fUm3cOtT+6GpDB+9YD3edPJyKOQsDuispZscmqq4mhtV5KVInbCJk+7tRuPyB3MpDHTVbwqTT2se88sfdlaenVqmJfIRTPsFQQ3hERfbR2BWFZ6LoxGI5sUTV/a0bGxE09+WMaYC+BqA1wI4DsAljLHjfLu9FsB659+lAL4e49iWIa4YSoZVFTU0UdTRndbciS+VUJBUFdzjrCy+5xSvIo0gn6pexY4IrQrHQrIZgcrkmFCYYytGlbNrLrE4D3QU+0n6o4bCTEM+s8J4QcctT+51J68l3WnsPlTAs/uncMaafnei6MiooRgaAWCbWMZmBI1AqzYNxTE3aTETDmtFGJFZs2xa2HWogKU9mapVabCPwNZgsj6B7WYWh1VRdb57o93nvvJHJ+Mf3lB/Ssin7cQwep637p1w36u1uichQb8VfZ9W5hFQ0InfYTxdii4IxGczqAxMs7RC7J0OYBvnfDvnvAzgegAX+/a5GMB3uM39AHoZY0sjHtsyRPVsqmS4zmIyGdkageY2W9dUuxUhFaZ7yXJbEpNkJ41AFCZiz9payR/0QKoqqxTNmsMetX4My4okCCq1hqgfQZizWPQRcHz0+j/g0u8+5D6gS7rTrgA9c22/242qmY5vtWgmc9vuyRs94qU/m8ToVNn9/Ws5i6PgNw3Vayxfy7xGE17JsLBrbMY1K/k/z/6cynkoqa4iCCo+glRCcbUJP/TdG9UIlvZk3EisWlCxxKmSXb5DdKrWmtT9/pBamdKNkktp7tgIy+Io6GZs38mKFpTGDqIVgmA5ALHG6rCzLco+UY5tGWJSx3TJQCqhwOKVSWuiaLeNIwuJpjJv6KQzGZBGkNZUqArzrGJHPIIg/IEljYCcrhRf3i5MKzh01I/rLA5wKIp4s1ItbHdCY/dNlJBMKJ7Suyv7s4JpaHauQVjfiSjE1Qj6ujSMzZSFrNtqH0Ec01BC8TqL6X7781esC9w/Stn0km5i16FC4MQSVGKcTIduJm65EjVUa6Kl797KBK0gyNY/WTSqfutak3qlZlYlTwJo3LkdBEUXinND5XOiX5fbP34ubvroy1s2LpFWCIKg2cNvLwjbJ8qx9gkYu5QxtpkxtnlkZCRol7qIN8glZ6wS1EIThmlhqmSgO5Nws2u1hOIRBLQSo9WZqjBkNDXUNFQrC9AUzgHMfiu6epiWVTOmncbpOotdh2JY1FBw+OiBqRJ6M5rn4cynNHtVyWbPNNRMsl6cEhMA0N+VsqOGAmLsXdOQEcM0lPCahmgSXh6yOqzlZ3FLS5QM7J8sBZ4jqKBg2bSgJZj7O5GPoF6BtjUDdo5MKyfWIMhcMlU0qn7rWtoIY7Zplq5vPed3Iyztsa/xHqGfNf2Gcdp3rhvKRfL/NEIrBMEwgJXC3ysA+Lu9hO0T5VgAAOf8Gs75Rs75xqGhoYYGSjfIzz98Fk5Z1ecp5ETvdacrpqGkqnhsn/QQUR6ByhgySW/Y2v7Jkjuh1kr+oImSTCz+5La5xrB4zeqV/rDRdB1Tjmc1aXE3fn3PeBE9PkGQ1ux4+myNJiGNIDrnGs31cFX4GA9sT0bDeEF3TX3iPUQBCrbzNZq5KenLI6BrlAsYU0JhtU1Dzu9ll4pAsGkoIBBANy1oiv07dSUTgkZg1ZxoKfdgtoMAKhqBXiUI6tn7NZW5SaP0zKZiZkHXYmmvnSy5VxAEJHCCylO3g1Z82wcBrGeMrWGMJQG8C8CNvn1uBPAeJ3roTADjnPM9EY9tGXSD0E1DGYdF3XSbhXRnNE/0hF8j4Jy7vY5JIygIN/nYdNn94WtqBL5M3npldGcbyymJHYa/V3FaU5HRVIxNBxfLM3yORmLfRBG9Wc19OMUkNluotm7CEIVUoxrBZNEA54hcHAywBYFucrdFozdqSHUyi+OZhowA05DogKUWjKev6a+pVdFi4+k9do0sfzN3oNKnmn430+KweOUeyKZUj0ZQK9mOWnr6m8G3GvIRTJYCNIJ6gkDQuOiZjZpAGIX+bBJJVcHu8Up0Fy1MZttkFpWm9TXOucEYuxzAzQBUANdyzp9kjF3mvH81gJsAvA7ANgAzAN5f69hmxxRGlSBwwzZN157bnU5UzCCq4pkcS4aFB3eM4br7XgBgPzAZXyJLyTAxlM9iJwo1ncW6r7ZPmEbAOW+4cXYcDIu7mlAQ/oxiwI6OoVrzfsTqmobJ3fh1w+LoyWhujoYYsmoL1dYJQ39wQCO47QQj1IwnaF8yE1bnEdimoXxE4aIlvI2LXLOCMIn87wfOwBO7xnHLU/uwcyy4xARlOgMVIXn0kuBqnmI2M/2W9Cx0JRO4//lRjDkhsrU0gmOX2uc/KaT6Z6sgk8lk0agSsPU0L7ExVClGxndUFIVhcU/KqxHo4WXf20FLRsE5vwn2ZC9uu1p4zQF8OOqxswVNDDmhHSNgq7cTQpMIMpX7b6CyYXpWtwmFTEOVbSXDcnMRapl6Kj4CwVkcIDgMK7oJoRkopyEMf7E5wHaKhjX8oOuU1hT85sm9nvd6MkmknYdV/MxsUm2paUic/Bt1FlPeQ3cDgmDfhP3gB0UNibV76qGFOIszSRW//ujLYVocK/uzWNmfxd3bDoQKU//9uLw3E2pzTjoCC6gENiSFiLHt+6bxzmvuQ1G3sLI/PJKlN5vE3X9zHhZ114/8aQZ6pqeKuit4v/DWDdg7Uaw7qYs+gpJRW8NplKXdGdzwyG58/FVHY9VAtuM0ggVVYoISTihUkVTGaaFscD6tubZy1fcAlgxviKXCqk1DJcNyH65aph7yEdAkT9EkfubKb1BPECRdJ7EgCLLJqmY+BE0eQavF3qzmbhc1Lr+/pVlEoUJNbz73i6fw2HD0RD3SCOI46UgQUChx2pdHYDh+h+gJZYrHNETfqyuVwLFLu92wZsC+t8OEqT/slFbrQYiaLtnP6V7d47SrfGbfFHaOzeD9L1tTc/wr+7OBRelaSVDU0NnrB/GRV66ve6xoBis7ZcNbzSlH9AEArrn7OQAVbVEsadNOFpYgKJnIpRKuD+DYpd1QGHD3swfcEtSD+aRb+wfwPjxlnyBQFTvBxn1gTAumY/oA6mgEfh+BprqOKrEUxlyFlNbXCOz3xMUVJU4FQRNX0ATQk9HcnryJOdIIpooGJgoGrr3nebz32gcin6PiO4rnIwAqE6boaKbV5lTRiBk1VG0aCrJ9Z1O2w92fxQpU944QBYgfWxBYnuOoj4RYunowl8IFxy2O9D1mk2TCDuwQfQRRK3tqHo3AmhWhdcVrj8ERA1mMOf2bnz8wDVVhWNmfrXPk3LCgBMHSnjROX9Pv/j2UT+Fl6wZx46O7sX+yBIUBA10pfPN9G/HJ1x6DJd1pjyAoGZbnAVMVux4RqeI08XdHEASG30cgaARiSYYoPYFbgeHUPQojKG+glkZANlczYEISncWq4CPoyWiug7UVTPtMQ6PTtrCPuhIHmtMIqPRD1uMstj97dLrsljeoh6Yw6JblLhDofgsyK9C2oIg1f6jvy9cPhn5mWvDX0HFBpqx6TeXnErp//L7AemhC459SjJLjcRnMpdwSNM8fmMaq/myssOTZpDNGMUf82Tlrce37TvNse+OJy/DiwRnc9vR+DORSUBWGFX1ZfOjcdWCMeerGl52wP0JVFI+Dk1bvrkZQYzVf0QjIR1BxFosRN3E7Eu2fKOL1X7k7Uk16kXolJt50sp3nN5SvPPj9XUlMFo3ApDJDWGH56clobi0XUfj0d4ULlkbwaARlA6POuaOuFKdKBr59zw4AjfkI9owXkE2qnqxbilQzLR7ZLKCpCjiv3DMzriCo/h7+zF8Rf1+NE1eEO3AzyUoUm+FzFv/7O0509xvMd44gWNmfxY7RGUwVDWScZM8oeKOGZsdHAAC9mYpPbfuBaawJiNhqFwtKEATxmpcsQVJV8PSeicDVDd0gfV1Jt3wwoTKvacjVCNL1ncVuJIarEVScxWWfXyIO3/v9i3hi1wR+8MCL2HlwJnJIqmHWFgQfOmctnvjsa9wGMoB9TQAEmodIqykHrEx7MpqrfqseQZDCoYIeqEU0gjgZzpRMtxVp1JXiJ378GJ7aY5fEyMfII8inE2DM7mTnn6zFSUY0QdaCTDKGKwiM0I5fWV9ROBG65/7i/CPxnT893T1vEB4fgUk+Anv/t5yyApc61Vwpa7YTWDvYhecPTGO6bLjO4ygkRR9BjN4TcenJ2oLAsjh2SEHQWfRkNJyx1jYXDQWsbsg01J9NouzEfxOKYj941HWLJm26CaNpBJXwUQph9Tqo49nMKURtMJfCy79wBz72w0ciHWfW0QgYY1UTaL8zkZHdU8St3RIgyHqzSdc0JEZEDXQlwXmwYGkEEoLJhIKyYbmmoXoF0MYLOp4bmcJvt+x3t4XV0glCUZgrOHI+849YMCyqIPBn+s6UTWRDOn5RSOpEQOcxOv4ly3twzlG1kzKDTEPib0Vjb7SY3GywdiiHkckS9owXIwt7gExDlTyC2XJs92aSGC/omCjqKOhmy0tJN8OCFwQAcPIq26MfZLMtuxqB5ob9ESpjWL8oB93k2D4yVSk7nFDrloxwfQRCghZN+kFNw6OyxwlZpHHf/vT+Wru72L0R4t0O5EANaneomxYUFjz+3owGqiSi+kxDAFpmHqrkhmgoGRWNoJ5t/u1X34tXfunOpiKYemmi9E1IS3sqGlVU0xDZrMm0M1M2kQmZgMmXMVHUcWCq5JZOB4In9DBE05BfIxDPEUM+zjprnZavjw2PR87RAOzr6wkfnYWoIcD2jU2VDLclbm8Mc+NsIwUBgBNX2NET+yZKVe/Ryn2gK+X4CARBoDBscI59bHjcU3Y4ral4fNd4aGq96TqLvT4Czjl0Q/ARxBQEe53sReo8FTUvxrB4rFUv4A3ZqzqfU78+iJ6MhiHHxPS+syqhh9Tda3SqtiDgnHsiq8Kgiaw7Y5copmsS1meZeGZf8y1DaZL3r5hFQdCbjTYR0D1C4x6dKoVGMdH2iYKB93zzAXz0+kcq3ztgQg8joymCaYgWLZXj6P6Ne8/MJtT7++B0Gct6oq+2gxLKZgP6vV8Ynfb83QlIQQDgBGcyP2d9uLrsagSG6Cy2+w1nk/akXxI0As457n1uFBddtQk/eWi46nz+ZvGphO0Q1E3u8xHEW5VSYStadUQVBFadqKEgxLR+P3oNQdCd0dCT0bDjyovw7jOPcLf356JpBP9w45NY88n6OYhFw47Vz2gqSrqFA85545Ty+NRFx+Kmj8Sv+Ei2c7/2IWoB0Z3FjmnIWRQ8vmscL1kWHPpJjuqJou76N2iBU44lCIJ8BJX747xjFgGA29WsE1jZn3U1lFUD0cMyxYSy2cojACq/zY7RGc/fnYAUBAAW5dPY/KkLcPn5R1a9d8oqO7Iil9KqfASqwqAqDOsX5/Hs/kmPRkBNwZ8/MI2P/+jRqhWsSbHZgrMYQJVDOq5piMZAeRG1CsmJ1Cs6F4RY8bH6fFZg7SKxhIefimmoWjMT+Y5T4qOeU7mkW0gn7C5zZdNy6yLFaQB+wvIeT9P2qFBXLX+xOtGu3xcjagiwf6N9E0XsGS+Glmyg6Kb9E5VyBvsn7dcV01D9xz6dVHFoRse9zx1wo9jE445anMeOKy9yzaqdQCqhunH5ceLzxYSy2cojACp+lRelRtC5DDqho36u+9PTcevHzkHKmUxEUw1NnCt6M9h9qCj0qK2+kfzhfLqvDLVYCdWfzRwHUv9djSDicfXKUAdR6fca5COwfQ5vcFaMl527Dscu7fb0IfBDD8poDY1gj1C4q97KnmrlJ1Xb7EZ5BXE0gqihpn4GSCOoUUIgaiSSmMPxiNOC8cQQQZBL2hFLv3/+oLttv6MRkI8hiumDktX+6H9+j7u32WXfozq328laJxJnVSxBoLjXZrZKTACVif95RyOIE5I823SOy79Dyac15NMakol9AIAZwQxCk/jyvgxufXqfO8EE3UgHp8ueScX0OYvzgpPPIwhirF6pSiRQ8RFEXeWbFocas6ZRNqmCsWAfgV1Lh+Erl5yMq951Ehhj4Jxj2NdfV0RTFfRktJqmoS1O1UzAFq61JmoSBClNxXhBx7Qv8S8KjQoC0m6C8gEHuuyexlELm5FJZrJouGUrwur7KApDd1rDk7srrRorGoFj4onQbU3MWv7vO7fjog1LcaTTuL6TWTuUwx1bR7AyRicvsaifrRHMjiCge2L7iO2D6iTTkBQEEXFLAwiCgFbQy3rsDOTdziQXpBEcnC571FXD5yym0NWRyRJEK1IpRmaxKEDINBRVJTDrVB8NgkJKjMvb+AAAIABJREFUf/DATqxfnMcbNix1JzfDtFznIm375OuOrXtOmiTDEMMi61UqLeoW0prdd7psWK4Qj6cRNGYmII0g6LNu+dg5sSKjSCN469fvxRWvPcYeV42wze5MAjsP2vdiMqE4gQxmPB+BT5M5c+1A5PG2kwuOXYxn9k3GMg35fQSzmVkMAMNjBWQ0ddbrL8VBmoYiEiQIFFcjsG86ascY5GzyP/h+HwHdJCOTJc+EfuVNT0cuzSz6L0gziOMjiGsaAuxwxQNTJXzkB3/whKrqdfobhNHflcTBGlFDE4L2US+8s2iQRqCgZJjubxdHI4gTjy5CE3WQIBjIpbB+cXjBNz+iyZLuo1o19imEdDCXAucctzy1D1/77baapSL8+Ju5dNLqtRYvXTeA737gjFilG+bKR5DWVDfAopP8A4AUBJGhm8MjCJxJdpnTiIaab6RDNAIRw+cjII3gG3dv91THnC6b2LTtQKQxGma18zTq3G7VKToXhR9u3omHXxxzxmK5na7iUK/MhJizUK/rVVE37ZwO0giEhutRabRlYb0ObnEQhfmByRK6fGUr/JAgWNmfcTPBtx+YrioVUQv/9+6kmPdWQ0XnDKdo5GyZhoDKc95pglUKgogkhaqRBD1Yq/qzSCgMD+ywHXQpTcG/vOUEz/H+bFm/j4AetEeHx/Fvtzzj2Xc4pNGIn7BG8lEwGhQE4qR961P78Jb/utcZC49ki/YzkKtjGirE0Ah0OxQwpSmYLBpOZJQtCKLkIQBouEHJSat6sbQnjb+84KiGjhc568gBvH7DUgDAyFSpbttMyiVY2ZfF9//sDAD2QiYoMSwM/73QaRNXK6EmPEUh6m+2IM2/0xzvTX1jxlg/Y+xWxtizzv9VsWSMsZWMsTsYY08zxp5kjH1UeO8zjLFdjLFHnH+va2Y8swmpdK7tHRVTTD6t4S2nLHe3J1UFl5y+Cm8+ubLNP7npvjyCoBXerR87B11JFS+MRhMEQaGmUVek9YrOhUHn99ec0U0rdqYyUCltHVRGGfBqBPV9BJWoIeo3MZBLweLVlThFWtEIqDut4b5PvtJT7bZRGGN4zfFLANimw3p9bt9+6kq88cRl+MDZa3DEQBeOWZLHZFF3HaJRBIE/W7zTTBmthBZ5FFU2WwllQEVT3Li6c8JugeY1gisA3M45Xw/gdudvPwaAj3POjwVwJoAPM8aOE97/D875Sc6/OelU1giLu20Ve+dYJepFbPz9tlNXuq9pUhfLBPt7+5q0OqsxWSYTClYNdOHFiJVEgzSCom6FTqoijfoIiLecssJ7vhjN2UX6u1IwLR5YKwfwRijV611QMixbEAiqPmUvB5VpJuLmU8wFdC8dmCrXjWS64LjFuOqSk90Q01wq4akSG2WiWz3gLYjWm+msFWwroetBZt9Ug+bAKDy3344YOuvI8BLg7aBZQXAxgOuc19cBeJN/B875Hs75w87rSQBPA1ju36/TodIAdvp6GjuuvMitJwMAxwckHYlOpyofAWkEdRrGH9GfjSwI3JaCPhtnFK3ANOOXmADs5j4A8DevORqnHtHnXifDakwjcMtMhJiHJou6G4YX9r2Kuol9E0XHR6B4fgc6tpafoDW1T1tLxhUEpdiF3vLpBCZLelWnsVqcsXYAmz5xnucchyt0PcjsO5s+gsvOtau2ntJBiXhA84JgMed8D2BP+AAW1dqZMbYawMkAfi9svpwx9hhj7Nog05Jw7KWMsc2Msc0jIyNNDjs+/V3JwAbuRL1VGpkmiEqtocpDee37Nnr20VQFS3vTnqbXtSDTULfvoY3S9cvkjWkEP/zQmdj0ifOQUBUctThfKT9t8sDrVI96hecmiwYWOQ63MNPQ5d9/GGd8/vaKaUjUCBwbbVh+Buc8djb3XCCWs87GDGnNpzVMORoBY9X2/zBW9FVCMDupplCr0XwRgbMZ1vnul67GjisvmrUQ1UapOxrG2G2MsScC/l0c54MYYzkAPwHwl5xzynb5OoB1AE4CsAfAl8KO55xfwznfyDnfODRUu4TubMAYwxJntRtmYz1qcXjCjV8Q+GsNAcD5xyzGmWsrNuWkqiCf1jBdNiKZd0j193fTCpowOeeeczZSdI4+iyaMhMJcAWdHDcU/X95tQh4cETRR1F0zXZhGcJsTxjo2oyOtKZ4V3kCIRrDrUAE3PLJrznpEx0U0M8ZNcsunbdNQ2an/1KgD/HCFnmd6RoP6PBzu1L2jOOcXhL3HGNvHGFvKOd/DGFsKILDmMWNMgy0Evsc5/6lw7n3CPv8D4JdxBj/XLOlO48WDM6GC4MbLz/ZMJFwwMoRrBN5zifHbWoIhl1LBOTCjm25MO+c88GEmByhNckctzuGZfVOY0asn1Zdd+VuoCsOmT5zvjqcZHwFgCzWKpDLMxvIIaOUbNslPFg1sWGGv6sM0Heo/AKBKI6iYhrwT/j/98in8+om9+M6fnu5ue/9Zq2OPf7YQwznrOYv95BxBYJizV1lzPkPXhBJCBzqo/eZc0exdcSOA9zqv3wvgBv8OzJ6xvgngac75v/veWyr8+WYATzQ5nlmFNIJkyASX1tTQMLsqjcCZMP1zL61eGbPzEXIp+3w3Pb4HRd3E5d9/GEd96teeXAP/OU0nNJIcuEEawZ7xIoYdxzfn3Gle39ztIGoEumXV7IAVBk14YZP8ZNFAd1pzWoQGaw1iDR+/IKA47mnfsfS7ff13zwEA/vnNL8E/vOH42OOfLUSNIKhFZS260xrKpoWpkhFbOB+3tBsnrwpvaXk4oLmCgJo6Hb6O8TCa9QBdCeD/GGMfAPAigLcDAGNsGYBvcM5fB+AsAO8G8DhjjNpl/a0TIfQFxthJsP1zOwB8qMnxzCpL65iG/DCnvkNPRsN4QfekrxuWHVXjX9mTfTKr2UlDVOLgb378GB7ZeQi/fGwPADt5bYOv5yyFs37u4pfAsrhbv6iezZssRHFLTPhRVeaavPQGV5/kFA3SCCyLY6pkIJ9OIJtUsWXvJD7x48fw+hOXYu1Qzo3iyqUTrrM55XMW0z7UO5Ygs9izTlRHJ6X/A97JP262M5nbRqfLsZul3/TR+CW45xv0jFGfgMEFqBE0JQg456MAXhmwfTeA1zmvNyGk4g3n/N3NfP5cU89H4OfD563D/skijl3ajS/evBXjBd1dkZohCVxkn6SkITFaY+veSsG1IGcnmYaW9WRwwooePOBUoKwVM09jAaJlnNbC6yNozNREgqAoaAS7DhVw1pW/xf+8x3amd6VUpDUVdz9rZ1z/cPNOdCVVPPm5C+33hUkzk/RqBMudYmTjBa8zmj6P2lnOZuRII4jjie8stq/HwemyNA0FQCXDt+ydRD6dqCqvsRCQd0UMXI0g4iQxkEvhq390Cla4k09lFWqEtIaklSjZgcVJTTSXBDWs8VeXpAmwbNaOGnrDVzYBiB5NEobKbI2Ac26XoW6RaWizk7H9jbu3A7BXx/5wxmlhf38fZHESpd9wzKcR0OdRwnGnTQaiIz9u+CiZF8emyy1JljvcoIZIuw4VMLQAtQFACoJYULRKmI8gDLI/i4IgTCOgSYtMAbm0KAgqdu2g5ir+EgK0+qtnGtq6z9Y0mjYNOYLN4nYeQSOTjqowJBOKxzRE34eSybJJ1ZPM50c89uz1Q65AzKcSyDg+A3/JD78pqtM0ApG4jXKaMQ0tBAaEBkEL0SwESEEQi6VOH9S4D5PbPlAQBHYJhgBB4JiGyG4p2oOnS5XJKighyu1A5UzINAHWCokUhUSzseJkWjItDt1oLKEMsCd60RFMAnPSaYCTTao4wpf5KkKr+3OPGkIulXCvx+rBLjDG0JvRMO7TCPyCoNM0ApG4yUikQYwXdCkIAkhrquuMJ9PtQuPwTRecBYbydhezRgXBIcEuHeojcExDZCISBYGoEQRN7n7TEK1q/RqBWHBNPOfoVO0WkfWg72NaHLrVWNE5wNsvF6iUfKDJO5tM4IgaPWkLZRN/dMYqfP7NduE/Kldx7FK79HNfNlmlEfgFaydqBKet7kNfNhnbhCf6FBpJ8lsI5FIJzJRNt5LwQkMKghioCsMR/dnYBbgqGoFQudQKbu5OGgGtzrtSwT6CYI3AZxpyfQSWb7+KIBBt61RGu1FIwzEsq+Ey1IDt4BW/K32vCcE0VKs5eUE3kRVW9K88dhEuPWctPnye3ZO6J6tV+Qj8IbadWG3zR5e9rKHjRJ9CXLPmQoE6v82XBjytRgqCmHz3g2cgF9dZR9myQi+DcB+BPYFVmtoHT6bBGoG3umSYj0AsuDYtjOn8Y2pWCKkLfZ+yYcHijUchZTTVI+j8Qi+bTITG0nPObUEgxN2nEir+VuiO1pfVsOOAt37TTNmEwiqhtIOHkYlA7DbWqLluoXDGAhUE8q6IyfLeTM0G7EGkEnYp5NGpsmumCKv2SeGjNKkGZRD3ZrWaGkHSrxH4BIEYejpdMpBUFVx6zlq8feNKNAN9n6Lh1UziktFsjcC0OLbtn6oSetmkimOX5vG5i70JX2XDQlG3wDmQrpF9G2YaEh2FcbN3OxlRKErTUDD/8c4T8eHz1jXckW6+I++KOSKXTuDae57Hhs/cAsBuVRlFI/CjqQx92aQ72YpQdUlaiWthGoEuagQmdMtCugUTBEUN0fkbDVXMJG0fwc1P7sWr/+POqn4M2aQKxhje89LVnu3/efsz+PJtdlOfbA1nL9XeESnopsdReDjV49FUxV0UZBZgHZ0ovPnkFfjr1xzT7mG0jYUp/tpALpXwVNTUzWDTEM0/4ntnru3H/dvtWPrutIZUQkGphkZAQoQmYn+fAnGFPVnUwTkaivn342oEztgaNUNkNBUjkyXsnyjC4sD2kSnP+2Edur52x3OVfWqY7zLJBAq6CcsptEfmpMM5YqQrqaJsWG6tJYlERC4P5ghR5aTaPkE2dOo7LGoE11/6UlzuODq7MxpSmhqsEVgcSaG6JGN2TH7JDNcI/t/3HgbQuBlHRHUFQfS690FkHY2AHNnDQjMgIFof4UwN0w6ZSshXUjJsc9LhHENOgrG3w1okSjoDKQjmCG9imAndtAKLvFHBOP975Ds4dmke6YQS7CMwqpO4qHG7SFBWcisyTkmwkbbSqJaRcdpz/m6rXcxW7NmcSigebemOv3oFPv4quy/w0p60mxxUS1iQIKDIJIoYOpw1AgoZ7peCQBKAFARzhFgRc7JooGxYgRFBVKvHP4fuciojnrl2AGlNDc0j8E++YklmwNYGrvz1lqpjW6kRFFwfQWPnpGJ6D+4YA1AJGwWqa/GvGezCxSfZDe/+6tVHu+0ZZ2p0ICMhUSib+PTPn8DrnRIbpBEcjtUn6X45nHsPSxpH+gjmiKxHEOgom1ZghMLZTi/Td562yrP9gy9fg+mSgbefuhKbnj0Q6COgxiMifkHwf5t3uhPsOzauwP9tHgbQfME5oFKiolnT0CWnr8LP/rDLLZonErTSXzWQxdOfuxCZpIqTV/XiwFQJL1sXHgZIZpKZsonv3v+Cu30wl8QnX3sMLjhucUPj7mRIEEgfgSQIqRHMEZaQzTtZCtcIVvZnsePKi3DqEd4yAuuGcrjqkpORSdqVN4u6iamSgS/f9oybcWs3HvFOvsmE4kkoE4XCaasr3dBa6yNozlkMVHdZ83+GH/IJrB3K4cbLz65p76+YhryRQxlNxYfOXYd1Q+Gd5uYrZA6UPgJJEFIQzBGGKUbqGCgJvQniktYUlAwLv92yH1++7Vm89ep7AdimIX+ceNLnIxDDIvPCZNtKHwE5YZvRMvx9l4m+Fpg2MsnqCqfA4ZVE5oeSDaVGIAmiKdMQY6wfwA8BrIbdWOYdnPOxgP12AJgEYAIwOOcb4xx/OGAIZR0mi06TmgZX4amErRFs3Wu3ft62f6pS+lkJ0AgEQTAmhLB2Zyo/f2s0Asoj8Ca2NUK3r8TD0Yvz+OzFx7fEfu93FhPHLY1X1XM+0gpBKjn8aPbpvwLA7Zzz9QBud/4O4zzO+UkkBBo4fl6jW6IgMFAyzIa7YKU1BUXdwtN7Ko1qyqaFsmkF+wgEbWRUFAQejWAW8gia0Aj8/QbSmoIz1w7gyEX5xgfoQIJgquStN9TJFUebhYIVwkxukoVNs0//xQCuc15fB+BNc3z8vGHD8h73tasRNGwaUlEyTDy9Z8LdVtTtQm/+c2pquEYgFlZrhWmI7PelWfARtLJ1ZMZxFu8ZL7rbjlnSvIDpZG64/Cx8+Z0nNV1qXHJ40qwgWMw53wMAzv9hVcs4gFsYYw8xxi5t4Ph5z19esB43Xn4WGAOmaoSPRiGVUGBxeyJb4jTLKekm9ICooZRPIxCzm8VVd0s1AqO5qCGgWiNItbA0ApWf+NnDuwAAl527Dj/80Etbdv5OZO1QDm86eXm7hyHpUOr6CBhjtwFYEvDW38X4nLM457sZY4sA3MoY28I5vyvG8XAEyKUAsGrVqjp7dx4JVcGGFb3IpRKYaNJZ3CNEfqzsz2DvRBFF3TYNVfkIfBrB6HQJp6/ux5+ft85jCmlFVUo3j6DcXB4BUHFkq04f5Fb2ByBnMTWqP2Ntf0eWnZZI5oq6Txfn/ALO+UsC/t0AYB9jbCkAOP/vDznHbuf//QB+BuB0561IxzvHXsM538g53zg0NBTnO3YU+VQCEwUdhsUbFgRiX1Vq2Vg07Gxl/zmrnMUzOtYvzuEVRy/yTK7JBpvIiCR8Reea8RFQJnafI/RaaRryC5VeKQQkC5xml1k3Aniv8/q9AG7w78AY62KM5ek1gFcDeCLq8Ycb+bSGA455ptHJTSyFsMwRBIWyCSMsocwxDZkWx9hM2S3DIIaStlIjaLYMNQBozrlorK3UCPyVRaU2IFnoNPt0XQngVYyxZwG8yvkbjLFljLGbnH0WA9jEGHsUwAMAfsU5/02t4w9n8umE2xKyUY1gkSAIVvTZnbqKuq0R+O3yomno0EwZnAN9AbHkLfERqP6EssY1AiqV0ZPVoKmspT4CP4u6F2Z7QomEaCqPgHM+CuCVAdt3A3id83o7gBPjHH84k08n3GqajQoCMWt2eR+ZhhwfQY0SE9SMJSipqJVRQ8Umaw0BwCmrenHJ6atw2blrce2m5/HSGiUjmmHHlRfNynklkvmErDU0x+TTGkanbY2gUXOHWGKZTCekEfiTuNKaiumyAdPiGJ0qO8dUZ9DORh5BM+dMqAr+5S128/nPXvySpsfm54w1/VjSIzUBiQSQgmDOyaUTbrp/K+zeFPlT1MlH4F3Zb1jRg29usvDY8CE3dLSvq9om3oqic4qv6FwrzjlbHO7hohJJHGStoTlGjI9vjSCwz1HSLcdH4D3nOeuHoDDgd1tHcHAmXCNophwE4fcRaLJRukQyL5AawRwjZsw26iMAgF9cfjbKplXRCAwTZaNaEPR1JXHMkm48svMQNjoVTYM0gpZWH21B0TmJRDJ3SEEwx4gaQVJtPDb+hBV2yYqpkl1KuehmFldPvst6Mxgem8HodBm5VCIwbLUVk3bCV3SuFcJFIpHMPvJJnWM8pqEWhESmE5XJ17CqNQIAWNKTwt6JIsZmyqFliGejH0ErIpEkEsnsIwXBHJNPCaahltjlFSQU5vRBrk4oA4Al3WkcmtGx+1BhVgVBwi06Z0FVWFXilkQi6UykIJhjci3WCAA7cohKKgetwhc7CVNP75l0w039hHX+igOdoxwQxiqRSDoX+bTOMR5ncYsmy7SmYMpp8B5sGrIFwVTJCMwqbhViJrHshCWRzB+kIJhj1gx2ua+biRoSSWsqJmoJAqGEQphG0ApUKQgkknmJFARzjJgV3KqKmna2sp0j4O9ZDFTqEQGzO0EzxtwCblIQSCTzBxk+2kZaZRrqy2p4YXQGQKVqp4gofPymoe/86el4ZOehlowDsENVxwu6FAQSyTxCagRt4MbLz8I7N66s6sLVKL1ZDSOTdv2isOgfMkP5TUPnHDWEj7xyfUvGAQDLe20zlBQEEsn8QWoEbWDDil5seFtvy87Xk0m6PQeCTEOA3RBn1AjPI2gVdH4pCCSS+YPUCA4DerNiJFJwGCjtk0vNruynkhdiG0yJRNLZSEFwGNAnCIKwTmNX/8mpeMfGFZ6opdlArIYqkUjmB00JAsZYP2PsVsbYs87/fQH7HM0Ye0T4N8EY+0vnvc8wxnYJ772umfEsVHozFTNMmGlo/eI8vvC2E6sa17Sat56yAgDwmuMXz+rnSCSS1tHsrHAFgNs55+sB3O787YFzvpVzfhLn/CQApwKYgd3AnvgPep9z/v/bO/Nou8YzDj+/jJUmEhEhhhAiXEJEIqRoEKG0hC5VEWPUWFMXWmWhlabSmlLjEmOWKcSU0mhXDdXSVkukRCk1dZBGRJGyhMjbP97vyF5Xbm5uztlnn3PP+6y1193jOc/dZ+/9ft/7fXvvWc23D1qnZ6ZGUPTzfTZbpwevT/4qA/v2KNQjCIKVp9xAMBaYlsanAfu1sv5o4BUze6PM7w0yrNFtWY0gHu0QBEFbKfeqsbaZzQNIf/u2sv5BwO3N5p0o6VlJNywvtRS0TraxOO/UTxAE7Y9WrxqSHpI0dznD2LZ8kaQuwL7AjMzsq4FNgG2AecDFK9j+GElPSXpqwYIFbfnqdk/fHsveOFZ0aigIgvqj1b6EZrZ7S8skzZfUz8zmSeoHvLWCj9oLmG1m8zOf/dm4pGuBB1bgMRWYCjB8+HBrzbuR6BWpoSAIyqDcq8bPgcPT+OHAzBWsO45maaEUPErsD8wt06fhibeCBUHQVsq9akwGxkh6GRiTppG0rqTPegBJ6paW39Ns+59Kek7Ss8CuwHfK9GlYuqZuo/Ge4CAI2kpZt5ma2UK8J1Dz+W8Ce2emPwTWXM56h5bz/cEyxo3oz02/fz3u6A2CoM3Es4baCed+bQsm7DiAPt27tr5yEARBhkgotxM6dBD91+zW+opBEATNiEAQBEHQ4EQgCIIgaHAiEARBEDQ4EQiCIAganAgEQRAEDU4EgiAIggYnAkEQBEGDI7P6e36bpAVAnu806AO8nePnt4Vw+Ty14gHh0hLhsnyKdtnQzNZqPrMuA0HeSHrKzIYX7QHhUsseEC4tES7Lp5ZcskRqKAiCoMGJQBAEQdDgRCBYPlOLFsgQLp+nVjwgXFoiXJZPLbl8RrQRBEEQNDhRIwiCIGhwIhAEQR0hKV5BV8PU6+/T0IGgln40SQ39W7RELfxGkmrpBU6dixYoIalP+lv4a/Ek9cyMF3nM1OV5XJfS5SBpS0m7AFjBDSSStpJ0WnJZWrDLNpKOlrROkR7JpUnSSCj2N5I0UtK1wHZFOTRzmQFcJGmLoi6+crpJuh2YCWBmnxbhkny2lzQTuE7SBEldizhmJI2QdAtwQTqv6+raWley5SCpg6SrgLuBsyRNlDS8tKwgrUnAj0uBqYiTW1JnSdcA1wOjgEmStq+2R3LpmS6804GJkiZJGliQy9F4D4/ZwDNFlnol9QWuAGbhd6WeAkxIy6pa+jXnwzTZR9LxyaPq55CkrYErgbuAGcBuQFWPl3RdOQ+4DngQf/3vt4Eh1fQol4YJBMAaQA+gCRgPLAROk9S92qXxTKrht8DPgB+Bl6wKOKG2Anqa2TAzOwQ/Joq6Bf4MvCfbEOBYYE1go4Jc+gNnm9nVZvZRkaVe/KLykpndCFwM3AOMlTTIzKyawUBSJ0n9gPnAUcDxknqZ2dICjt0RwN/N7Gbg18AXgH9kXHPfL+na8QZwhJndihfuNgQKT5e1hXYdCCRtK2lQmuwJfAnoZmYL8JrBO3j0zv2gSS6bApjZknTS7AlcC7wl6Vtp2dIquZT2y6fAgak0/nVgB2C0pKFp3bxdBkhaLU1eC5wLYGavAL3wQJU7yaNrGu8NDAb+JGk3Sb+SdFbaP9XYJ+Mk/VDSvmnWM8BwSZuY2QfAn4Gn8GCZa/os47JP+q4lZjYPGAC8DjwGnJncci1QZVzGpln3A/tLmgQ8B6wPXCbpe8k1l/0iaVSzWvN0YE5KSy0EFgH98vju3DCzdjfgB+kvgD8ATwJj0vyb8VIeeBVuNP4jrltFl90yyy7ESzHbAn/Dq7frF7BfJgO3Am8BhwIT8ZNsUI4uG+FV6YfxoLxZZlmX9PdGYN+cj5XmHk1p/nVp3mXAWOBIYA4wJEcXAcfhF/4j0zHxrXSMnANcltbrAOwEXA30q6LLkcAX8RLvlLTevsD7eAqtK9C5Si7HZI7pnwKHpelR6dgdmYNHD7w29g5wA9C75JdZpzPw+zzPnTyGdlMjaFZKOx2YY2Yj8QatCWn+DcCOkgaY2RK8evsRsBoVZAUu9+EnNpK64aWGAXiqam2gr5n9q5L56JVxAb4PvAAcYF7NngK8BuxYKY8WXJ40s9HAo3ibwJZpWSkNsx7wz7RtxY7VFXg8AvxI0gDgPLw28qaZzTRPy8zCg0IumF9JRgKT0/d9G9gFL7DMAgZK2t285L0Q3z/vVdFld2Bn4L/AAEn344WZx4A3zGyxmX1SJZdRkvYys9fwdoF/pdWfxgs0iyvtAXyMHyOHAG8CB2T8SjQB883sJUk9JI3IwaPitJtAgJeaSif5B0DpgFwdeCE1Oj6BV6svAjCzuXjpptIHTUsuPZNLk3mD2xLgT0B3vKGrv6StrbL56BW5zJW0RTqQFwPfBDCv3q4H/LWCHlmXUhvJ8+n7rsDzvQdL6mveVjIQeMfMnkkNkudI6pWzx5XAMOAYYAFeKzggs11fvLRXMSQdllINvdOsF4D1JHUys4eAuXi6bgFwGzAl7ZvReEm5SxVdnsVrIoOAfwOvAsPMbB9gA0nDquyyS2pI/xVwXjrGDwK2xANlJT16mdli/Jh4CHgJT9cNSuuVjqXewIeSjsCPla2q2YazqtR9IJA0RtKvgQslHZguao8Dm0p6BvgK3nBzG15tvADoJ+kKSXPxhp73KvFjraSJ10e6AAAF2klEQVRLJ+BGSXvjqYihZnasmc3G8+PvluvRBpeOwDRJewC/BPaUdJGk3+EB49WcXJbg1euhkoZIGoJf8DbEG4gBNga2k/Qonn6YbmZl7ZuV9Hgebyjub2ZnAf+QNFnSH/GT/PlyHJKHJPVL/9vheI3wckmr4zWgvizr/TIdb69Y08xuwVN4Z+IXvO9WYJ+0xeVOYAv82DjVzE4xs/+lZaPN7OkqutyBB6R1zewaPDA9iO+XCWa2yu8racHjSkl9zDsOfIynV98CDgRvO0mb7wmMA74MjDez65vVGGqTonNT5Qz4QfEkXl0fil/sT0/LNgPuyax7DnBFGl8bbziuWP65jS7nARdnpjsAHQpyORe4NI1vgzdA7p+jy+3ACXi+9RzgATxADU+eJ6ftxuMX6d0L8vhO2m51YHNgjwp5dEx/BwG3pPFOwFXANDzHfAPeVtMzLb8JmJT5jC4FukwDzk/jqtRxW4bLxDTeGVgnR4/Lgbubrbt/8huId0IBv658s1LnT7WGWrpjcqUo5YrN86TbA0+b2cy07CHgEkk34xeRf6Y0zAt4HvpUSTKz+Xj7QFEuD2dczCrQ26IMl0eSSwczm4M3iObtcjEww8wmStrYzF5Ny55gWZpuunl3vKI9FpnZi8CLZbp0As4HOkqahQeYT5PfEkknAvPwEvdtwH54L5gLgKVkUlLmJdKiXD7FAyrmV76ySrsVcPljWvcT4D85epwMvClplJk9lubfK6kJr013l7SrmVU0dVgt6io1JOlIvFFoYpr1HDBO0kZpujOezpiId+HqDZws6RTgGjy3Fy7Ld6lIf/SVcOkEvAJcmqZfS9sdg/dLnw3l361aQY+yq/WSRuGNmGsAf09OnwC7KjUmpmB1PvAT8xz4VGAnSU+m7X5Trke4lOVhyeMHme2+AZyNFzK3TgWr+qToKsnKDniD6n34XZWzgc3T/Cl4Nf8J4Ba8p8eDeDe3JuAkvAq5Q7jUlMsvgLXT8lPxRvzt2pNHxmdn4NDM9FXA8cAReC0FvFC2Dt6FeKM0rxewXrjk69JGjzuBAZntdq7kPilqKFygjT9Y//R3MnBHGu+Il3B3StMbpAtcRfKo4ZKby01A1zTdrb16lD4X72Nfyj+PBy5I43OAk9L4cOD2nH+fcKlRjyKHukoNmVnp9vEpeD/mPc1TCO+Z2eNp2XF4N8lcHwkQLmW7lLrPYsueXdPuPEqfa97HvrTvx+DdQcFvkGqS9ABeW5mdh0O41L5HoRQdiVZ1wHu3PJaZHoHfPDaLCvQeCJf241JDHh3xFMODwMA0byCe6tiJCqdewqU+PYoY6vJVlamHy1JJd+E9ChbjDZ4vmz+jJlzCpaY8kkvpBrDrgHvxO94X4qmH98OlWJda8SiEoiNRGdG7G/70zrdJfc/DJVxq2SO57IB3B30cOCpcasulVjyqPdTdfQQZTsDzdWPMb/0Ol3CpdQ/wLq1nA5eES0261IpHVanL1BAsq/IX7QHh0hK14lIrHkFQq9RtIAiCIAgqQ111Hw2CIAgqTwSCIAiCBicCQRAEQYMTgSAIgqDBiUAQBK0gqZekE9L4uunmtCBoN0SvoSBohfTo6gfMbHDBKkGQC/V8Q1kQVIvJwCaS5gAvA01mNlj+Xtr98GfUDMZfdNMFf4vWYmBvM3tH0ibAlcBa+EPujjZ/2U0Q1ASRGgqC1jkTeMXMtgHOaLZsMHAw/iC7ScCHZjYUf6ftYWmdqfjzaoYBp+PPuw+CmiFqBEFQHo+a2SJgkaT3gPvT/OeArSV1x99jOyPzAriu1dcMgpaJQBAE5ZF9Hs3SzPRS/PzqALybahNBUJNEaigIWmcR0GNVNjR/fPFr6f22yBlSSbkgKJcIBEHQCma2EHhC0lzgwlX4iPHAUZL+AjwPjK2kXxCUS3QfDYIgaHCiRhAEQdDgRCAIgiBocCIQBEEQNDgRCIIgCBqcCARBEAQNTgSCIAiCBicCQRAEQYMTgSAIgqDB+T/xiwvjSH2hdwAAAABJRU5ErkJggg==\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. Compute linear trend"
]
},
{
"cell_type": "code",
"execution_count": 14,
"metadata": {},
"outputs": [],
"source": [
"x = np.arange(len(ssta))\n",
"y = ssta"
]
},
{
"cell_type": "code",
"execution_count": 17,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"0.01729 K/yr\n"
]
}
],
"source": [
"slope = np.polyfit(x, y, 1)[0] # in degC/month\n",
"print(f'{np.round(slope * 12, 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