Skip to content

Instantly share code, notes, and snippets.

@tam203
Last active March 18, 2020 10:32
Show Gist options
  • Save tam203/0a1bced8416cc0e77885569c35cb10b8 to your computer and use it in GitHub Desktop.
Save tam203/0a1bced8416cc0e77885569c35cb10b8 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [],
"source": [
"import zarr\n",
"import xarray as xr\n",
"%matplotlib inline"
]
},
{
"cell_type": "code",
"execution_count": 12,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"[<matplotlib.lines.Line2D at 0x7f70f056e950>]"
]
},
"execution_count": 12,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYQAAAEiCAYAAAD5+KUgAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAgAElEQVR4nOydZ5jdxNWA37O77r3i7sW4YRtsjDHdGEwvgQRISCCU8IUOgZACCSROCC0hIRAChB5aEmISIPTeDcbGNsbYgLsN7r2X3fl+SLqrqyvpSveqrXfe59lndaXRzEiamTNz5swZUUqh0Wg0Gk1F2hnQaDQaTTbQAkGj0Wg0gBYIGo1GozHRAkGj0Wg0gBYIGo1GozHRAkGj0Wg0gBYIiMgAEZksIutF5LK085MmYvCgiKwWkQlp5ydqRORgEfk87Xxo0kdENohIH49rZ4vIuwHiUCLSt8T0TxeRl0u5N04avEAAfga8qZRqpZS6Pe3MFENEHhKR38UU/UHAEUAPpdTImNJIDGeFVUq9o5QakGaeiiEi14nINBHZISJjXa5fKiJzRWSdiEwUkYNcwrQXkeXORk1EThCRT83G8H0RGVQkL/eIyOciUisiZzuunSUik8x8LBKR34tIVZH4rhCRJSKyVkQeEJEmtmuXmM+zVUQe8osnCpRSLZVSc+JOB0BEqs2ymHs/SqnHlFJHJpF+GLRAgN7A9FJuLFYBskTAvPYG5imlNsYUv6Y4szA6Kc85L4jIvsBNwClAG+B+4L8iUukIejMww3FvP+Ax4AKgLfA/4Jki320qcBHwscu15sDlQEdgX2AM8BOviETkKOAqM1w10Af4jS3I18DvgAd88lM2upwWQSnVYP+A14EaYAuwAeiPUdEeBpYD84FrgAoz/NnAe8CtwCrgd+b5H2BUwNXAS0BvWxqDgVfM8EuBX5jnRwLjgTXAYuAOoLF5Tcw0lgFrgU+AIcB5wHZgm5nf/xV5vnnAz837twJVQDfgSfP55gKXmWHPNd9DjRn3b8zzxwNTzHy+D+xZSvxm+LHAE+b7XY8hiEfYrvcE/mPeuxK4w3bN8x17PPvbgAI2ms/zHWA0sMiR/5+a+d+I0cDuArxg5u9VoJ0t/H7mO1iD0ViOjrFsPgqMdZz7DjDB9ruF+Yxdbef2N8vVOcC7tvOXAM/ZflcAm4ExAfLyLnB2kTA/9iuPwOPADbbfY4AlLuF+BzwU8l11wBBw64CPzDjsz66Ai4Evgbm2c31t9z9j3j8BuM5+v0+69jiOAyabcSy0fztggRl2g/m3P0ZbUjSNpP9Sz0Daf8CbwP/Zfj8MPA20wujJfAGca147G9gBXIrR+DUDTsLo1e1unrsGeN8M3wqjsb8SaGr+3te8trfZwFSZ6cwALjevHQVMwujJiRl3V/PaQ5iCKMCzzcNozHuaea0w4/0V0BijlzYHOMr2fPaKNBxDKO0LVAJnmXE2KTH+sRhC51gzvhuBD8xrlRiN7K0YDV1T4CDzmuc7LvL8uQpr/h5NoUD4AEMIdDef9WNgL6AJRofh12bY7hhC6ljzOY8wf3fySPtZDMHh9vdsgLy7CYTW5vu1vselGI2Q2N7hxxhly/ktLwWet/2uNL/FjwLkJYhAeAq4yef6VOA7tt8dze/TwRGuFIHwT/OvOTAIo0F2CoRXgPZAM2fZMO99wix3Q4CvCC8QRgN7mGVjT4zO30nmtWozbJXt3rODpJH0X+oZSPsPm0AwK8lWYJDt+vkYcwzWR1zguP8FTIFh/q4ANmGoX74LTA6Yj8uB/5rHh2EIov0wRye2cA8RTiD8wPZ7X5f8Xw08aHs+e0W6C7jOEf5z4JAS4x8LvGq7NgjYbB7vjzEyqHJ5Ds93XOT5gwiE022/nwTusv2+FHjKPP458Igj/peAs2Iql24CQYBfYIwSdwArgH1s16+w8u/yLQdijIJGYwjra4Fa4OoAefEVCBijkUVAR58ws4Gjbb8bmd+n2hEulEDAqLPbgQGOOJwC4TC3smG7f6Dt2g2EFAgu1/4M3GoeV1NPBIKeQ8inI0ZlmW87Nx+jd2ix0HFPb+A2EVkjImswVENi3tMToyIUICL9ReRZc5JtHUYh7AiglHodQ4X0V2CpObnXusRnsue3N9DNyquZ319g9JDd6A1c6QjfE0MtVGr8S2zHm4Cmpl63JzBfKbXDIx9e77hcltqON7v8bmnLw6mOZzsI6BpBHoLyfxiqs8EY5fQM4FkR6SYi3YDLgF+63aiUmokxwrsDY9TaEfgMoyG3rG6sv15BMyQiJ2HMaxyjlFphnjvdFtcLZtANGCMcC+t4fdC0POiEMWq0l0NnHfU653X/fI+wnojIviLyhjmZvxZjrqZj2HjSRguEfFZg9BZ62871whhCWijHPQuB85VSbW1/zZRS75vXdvNI6y5gJtBPKdUao+GUXCJK3a6U2huj8vfH0HW7pV8Me/iFGDpUe15bKaWO9bh3IXC9I3xzpdQ/IorfmVYvj0k/v3ecFAsxRgj2PLRQSt3kFlhEXnA0shtcGsmwDMXQ03+hlKpVSr2I0bgfgDEn1RX4TESWALcBI80ORyWAUmqcUmqIUqoD8GuMcv6Rea2l7W9BkMyIyNHAvcAJSqlp1nllWNBYcR1jnp5u5t/+LEuVUitLfBcWyzFGSz1s53q6hPOqN9b99nsCC0Qbj2PMQ/RUSrUB7qauPoets6mhBYINpVQNhi7xehFpJSK9MSbLHvW57W7gahEZDCAibUTkVPPas0AXEblcRJqYce5rXmuFMQG1QUQGAhdaEYrIPmaPoxHGMN+a7AWjB+tqPx2ACcA6Efm5iDQTkUoRGSIi+3iEvxe4wMyLiEgLETlORFpFFL/z3sXATWY6TUXkQPOa3zv2o5x35eRR4AQROcp8rqYiMlpEergFVkod42hkW7o0kgWISCMRaYpRN6vMdCwroo+A40Skj/k9jsDoLHyKoVarBoaZf7/CmF8YZpZrRGRvM++dgL9hCJeZPnlpbOZFgEZmXirMa4dhWC2drJQKsmblYeBcERkkIu0w5oEesqVVZaZVCVjvt8p2XYnIaGek5rP9BxgrIs3NunRmgPx43T8IYyQVllbAKqXUFhEZCXzPdm05hnouqrIYG1ogFHIpRiM8B0N3+jg+pnBKqf9imPn901T9fAocY15bjzH5eAKGquRL4FDz1p9gFJr1GA3vv2zRtjbPrcYYvq4EbjGv3Q8MMtUWT4V5MLPwn4DRYMzFGBHdh2FZ5RZ+IvBDDDXDaoyJ3bOjit/j3r4YVhmLMKxqfN9xEcYCfzff1bcDhPfL30LgRIyR3HKMEcNPib4O3YuhqvouhvpnM/B989rDGBOgb2J0Jm7HGDnNVEptVUotsf4wrNO2m8cWt2FMan9u/v9hkby8bKZ/AHCPeTzKvHYtxnd9PsjIxxzN/B54A6NMz8cYpVhcY8Z/FYYqbLN5DlPobgCm4c4lZl6WAI8A/8CYCwzKJRiqwSUYQurBEPdaXAT8VkTWYwjjJ6wLSqlNwPXAe2ZZ3M95s4hMF5HTzeNedtWdqYIryTQ+LJZ1gkaj0WQSETkDGKyUujpg+JuBLkqpUnr6DRotEDQaTb3GVBM1xhhB7AM8j2E5GGoErTFm1zX1FHNI+ZnH5UFBJwfrKyJyMIbuvAClVEu385qdklYYaqJuGGtJ/oixlqgsGmL50iMEjUaj0QB6Ulmj0Wg0JlogaDQajQaox3MIHTt2VNXV1WlnQ6PRaOoVkyZNWqGU6uR2rd4KhOrqaiZOnJh2NjQajaZeISKerjm0ykij0Wg0gBYIGo1GozHRAkGj0Wg0gBYIGo1GozHRAkGj0Wg0gBYIGo1GozHRAkGj0Wg0gBYImoR5ZPw8jrntnbSzodFoXKi3C9M09ZNrn05knw+NRlMCeoSg0Wg0GkALBI1Go9GYaIGg0WgaPDOXrGPpui1pZyN19ByCRqNp8Bz953eoEJhz43FpZyVVAo0QRKStiIwTkZkiMkNE9jfPXyoin4vIdBH5vS381SIyy7x2lO383iIyzbx2u4iIeb6JiPzLPP+hiFRH+5gNi9paxaT5q9POhi+1tXqnPk22yGKRXLByE3tf9woLV21KJL2gKqPbgBeVUgOBocAMETkUOBHYUyk1GLgFQEQGAacBg4GjgTtFpNKM5y7gPKCf+Xe0ef5cYLVSqi9wK3BzuQ/WkLn77dmcfNf7jJ+9Mu2seFKrt27VaIoybtJCVm7cxpMfL0okvaICQURaA6OA+wGUUtuUUmuAC4GblFJbzfPLzFtOBP6plNqqlJoLzAJGikhXoLVSarwyNnJ+GDjJds/fzeNxwBhr9KAJzxdL1gOwZN3mlHPiTY0WCJng2U++5uwHJ6SdDY0HFRVGM5jU6CXICKEPsBx4UEQmi8h9ItIC6A8cbKp43hKRfczw3YGFtvsXmee6m8fO83n3KKV2AGuBDiU+U4OnPjS1Wh5kg0sen8ybny9POxsaDyrMfnFSKtYgAqEKGA7cpZTaC9gIXGWebwfsB/wUeMLs1bv17JXPeYpcyyEi54nIRBGZuHy5LsReWI2tuL7WbKAFws7J3BUbUfrjlsV7s1ZQfdVzLF67mdWbtgHGSC4JggiERcAipdSH5u9xGAJiEfAfZTABqAU6mud72u7vAXxtnu/hch77PSJSBbQBVjkzopS6Ryk1Qik1olMn1y1BNdikbHblASrFcczYZ6Zz8l3vp5b+zsqnX63l0Fve5L535qadlbL51p3vcfV/pqWS9uMTFgAwcd5qFq8xTGHnrczIpLJSagmwUEQGmKfGAJ8BTwGHAYhIf6AxsAJ4BjjNtBzaFWPyeIJSajGwXkT2M0cSZwJPm3E+A5xlHp8CvK50N6NkrFdnTcMsW7eFg3//OvNWbEwzW3mkadHx0PvzMm+FVR9ZtNpotCbOL+jL1Ts+XrCGf5gNc9JY/bhapRjV3+j4njSsWyJpB7UyuhR4TEQ+AYYBNwAPAH1E5FPgn8BZ5mhhOvAEhtB4EbhYKVVjxnMhcB/GRPNs4AXz/P1ABxGZBfwYQyVVL1m0ehMPvZduD8nZ1v7vk8UsXLWZh96fl0Z2XNHyPlts3LqDHTW1ZcVhdUCyaL4ZlG07ynsHUWBVjcoKyY3ym1RVet8QIYEWpimlpgAjXC6d4RH+euB6l/MTgSEu57cApwbJS9b5wUMf8cXSDRy7Z1c6t2qaal6snkaVaalQk6GamqGs1CtWb9xGiyZVNK6K1snA4F+/xKEDOvHgOSNLjiPpCdA4uOixSdx31j7FA8aIZZJdIVI3H5iQ+le7roiYDVt2ACn3NBz10TJdy5Kppx4hlMZe173CRY99HEvcb5RpbVRptiZZKmfFqK1V3PLS57nfE1NWJa7ZtI0XPl0CQIXUzbVpgVBPaWT23LbXpF8prEJUaR7UZCBPFvWozcgcr85YGkk8KzdsjSQeC0tlVJ++7eSFq7njjVm539/fr3fieZg4b1XuW3y2eF3ufP5SrGQkghYIEZNrfFMcNjsteKxylSWroyysVL7+uc/SzkKqbNi6I9L4rOI1cd6qejMCdGYzjSpyyt3jOfXu8QXn7SqjpNACISIWrtrE89MW5xrdLFSITK9DSDsDwL07gXlkOcRVPjZuq+GZqcnYzZdLgUOElHpNc0wLQPs3qaxI3oRcezuNiOP/8i5rN2+nX+eWQLqTpl6yKAMyKkcWRgj1jag7GVE3MvbGdd4Kf7v5mlpFZUX6HRZnHtLOkv2biEiu0iaVLT1CiIi1m7cD8OWyDUDKKiOHZUL61a4QLQ/Ck3XjHXs581t4+PWazez2i+d54qOFnmGSotIhFZMYVe+oqeWet2ezdUdNwTV76hUiiY8QtECIiUkLVrN20/ZU8zA3QwvRnGiBEJ6oRlVrN21n7DPT2VbmugM//ITXnOVGuXx66lexpR+UCkcLmETD++9Ji7jh+Znc9ebsgmv2UVYaoxUtEGLi2qc+5fT7P0glbat39gebOV3W0Cqj8Fh2/nv3bldWPH94eSYPvT+PpydH2yCHbUyVgvkrN3LpPybnmWkrpVKbg0uiDbZM0zds2VHwnPZ3mLcOQVsZ1X8+/Wpd8UD1kIWrNrHr1c/xxdL1vuG27ajlu/d8wMcLCm27k67uXy5dz6MfzE841WiprBBaN61ij+5tyopnh2l+7NaLH3DNCzwxsTRVTl6jFaBBVwqu/s80/jf1aybMXWWeU+x69fPc8PyMkvJQLhUJd8v9XpNgd0OTTH60QKjnKKUKVobaC9noP7wReZovfLoYpeDfRRqOeSs3Mn7OSr515/vc/26+RU/Sq1mPvf0drnnq00TTjIOKCil5dPX5kvW8/UXd4jO3Rmbrjlqu+19p5rgLbLt6Bfm8InV5sEa11txbqRZgUxeu4as1wfcBSdoSb3tNLbe++kXut/M15c/DFLqDfuKjhVRf9RzbY1L3aYFQzznxr+/R5xfPe163e0mMysOoVYmKVXp74b7u2fxGJmmNQBYWCkaB0Wss7d6j/vw2Zz4wIfcuvJrCUgWO8xsHwbJGspIs9yud+Nf3OPCm1wOH91qzExcfz1/Npm01Ztr5lmMT5q7KE2b272DNLdzwgjFystROUaPNTus5nyxaW3CuoNcRuXmhmU4xgeCTblrur9NcH/LS9CXs3qU1vTo0LzkOESn73eW2Y/T4QJ1bl+aDq8qmbvHLo3XtfdsWr3965QtG9e+UeEfBmV5FzBJh0/Z8yyJ7p+rbf8tfnKZU8h0nPUJIgJ+Nm5qaTjQO6rxahi+tVn1Ly4TSmWW7CiVuzn9kEkfc+lZZcZQzQgjKafv0LB7I5P5351J91XPU1ipO3KvORbPf93W7NmXhGmMyOeGOgrMMJ23Y4y84vUdMcb0lLRAS4ImJi7jn7TmJpRd3g+E20bVuy3Y2b6txhCseR9I4Ha8tWp3svtNby3R6KDbb9LgI00m+yVRh7KhVDO9VnvVTbQo9YqdwinqAMHXhGrbYRwWO9LJWR7RAKIGl67ZkwjVF2tgn5PYc+zKH/ym/9+vWExSfa2mQ5s5tpWAsXo0ozx7xhJlorZtPUoEbc6/8p2GK7MxLlCqjr9Zs5sS/vse1HsYMxQya7DnTVkYZY8v2GqZ/vZbpX69l3xtey21zl03SaeSc1h1ujW2dR8y05hD8fyfF6o3b2Of6V5nmMgfkRxIqo1IaH6XyS51vz9fjfE1tcKHixizTS0AYinVMXvx0cYm5gXWm94JpX9V9Y3udkGLO61ShG/245YIWCAH5zf+mc9zt7/LerBVA/oSYnbbNGyWZLVey0ui5pWsV6Kz0y9PKx/g5K1m+fit3vjmreGAbNvc2ZRNJNLk5oRCxeQStLXMOwTlCDZSVgoVh+U3uPyaEW5OxeO1mvg5o9ir4j1BXbtzGzS/ONMOaHalQuQmPFggBadu8MVC3lZ29INlt6p1DzuXry/c5v3LD1pL3Q969a+uy03cSdBcnt0ZCSmlAMNxwxOIfKiFpWdDwmP/DvgehfCujujyVH4f1HDWO1cVBrIycZGIOoeB6uAztf+PrHGCavbp2iPK81/k/79J1W2z3hcpGyWiBEJCOLZsA5LYurLWN5L7x13dzx87vdkgEC8MOvPl1Rt/yZuDw+UN3xVX/mVZ2HvLj97Zjn718Q86Hkv+EWfD05q3YyKG3vMmtr3xRPLALK2wbwaQ1Z1A4eVnaZjIVLo2I2+JEO5u27eA3/5tecN7rXRS4hPYhJ+Brg79Zr2euqVWJzSOs37KdBSs3ubqOsC+4LCc/bm/E3mGskOCLDMXxPy60QAiIVXDcJoLsLiqcdWnTtkKPhmHZsj2cZUpeT03hehwFbu3GmD++xaGm8HIdIdgmIYOyzBxlfTjXXU1XjF/4CMSkxEOtY0RZqvmtiBTc89Nxn/guTnzg3bk8+N684GmEy1IdyuPYGczjmmF2mgzfuvN9Rv3hjYJ3qZTxPi1qSzQKu/DRSSxZa/Tw7QLWLrgNlZE3i9fWqZ70CCGj5Hp2Hp9yxYZtSWbHFXvOPrf5G4qqUAW3JnE5GXBRmx1LCIdRGX21ZjObtplOxGw7g6U1n5InEJSy9RSDZejRD+Yz8vpXzTvy7xk3aZHvvV5eDv76RqG3TQhXTnIrjFXw0Vdck8phsNzUO0cIzo5KqftDv/DpkpxL/Lz4HOplv+gf/SB5wxUtEAJifbiKEof6OxPWo++oVazb4u3i20cehHp/dUI4OAfe9DpnPTABgDG77+IZLqnvaE+nRqnQ5rfXPPUpy9ZvNRprZTjrC+rPJqy/tjDBle1/uWanNUoxaf6qEKmXj/P9OwVCOdZwlhNB+/u3x19RbIhgI4warxy0QAjA/je+xvXmSuOsWcm4Ef/CNOP/g+/NY8+xL3uGi2pS2apQkxes4fJ/Tg5830fzVjNl4RpfHztJ6az//OqXueMnJi5i9nL3HmoQFq7exBG3vh149XusHjy91D+lRKXgBw9NLCs7odN05PSG52fmXy+jeGw39U3iIRAI4YZE75iWMF8uXc9ev305b2bfYvHaunPODTXqE1F4dty6o4YP5gTT5TttqO3sCKGctU/EPTUl3F69z08r3Y48Su5+q049c+1Tn+bMlsO2N0rBSlMtOXFeoVvx+9+dyyPj5+WdC9u5DNMbtRo0p/6/FEGXxsK0OBdIWnHb6519UJfEmpKw1OPmLVoeen8eqzdt5+XpSwC47505VF/1XME2d1lXGc1cso53zbUScXDDczN4K6D/n1/+t3Ay16ocJ981noWr/PfdtXCa8voJGoBJ8+sayixuHwp1e/mGbZCKTbxe9+xnXPt0vkVR2NW3Jc0hUH6dSGPb2TiFkJvlV02eyii4lVFSb0YLBBNnJbC2t1u3eYcjXLaVRkf/+Z3IK1ZtbZ2N+azlwVeDzl7uv3YiaFzOb3O1zWrovVkr+Mm/p+ZdP/mu9z3jyoogb97YWM8ysEurQOErcqq28L3vJOYQjOPyXm6jynibow1bd7BqY77RR7F3Wc4TWY29vfzmrV9SqqhBgNt9caIFggPl+O9sjKzKlZWGJQn6/OL5nCnee7NKM/20sDdOQRse5zd45bMluePT7/uQcZMWeVeYYv5iYviOkxes5oAbX/OdcG/V1FjRHnQ7TKuXb+9RKhQLVhYfZYV+xjAqI2WpjPxNnHfU1LJ283bfifBdSnS7HZQhv36J4de94shXfBXZ6piJyzmAxz5cwI0vzCQIbmtP4kALBJPc0vCcGV3hx7SHywJ2O+W4CdqTsfCq+KVYS1Q6urg7XEZAXoMi5/dSKPbsUd4WlMX40ytf8PXaLbwxc5lnGOuRgr6NOoFQd04pfIVOkjidTjg/x6T5qxn6m5eZMHdVpsbWcaqpcsLbVubt6a3ZFNxEXauMEqZu0xezx5M7n19lcyOEhPLlx7G3vZN2FgD33kq/X74QW3puwibMloLxO4czCsmP/jnFM0ypen37HIJS3jrw9Vu2u+5lHSitEGGDvsqqSiPWJNcaBGF7rALB+G9/nzMW160LalIVvPktGCGUkS8/9I5pJs5KoFw+JhQKjjRZvamudzj8ulc454DqyNN4afqSomHC9LLsi8SCNopOgyS35Dw1RlIYLu71sEGeKjdCCPgO8vXQ7sd2znt4EuPnrGTmdUcHit8rrWLk0nfMKjvzZX3rMx+YQPe2zVzjevNz7xFVXNSUuhQ5UNyFH6dlk8rc8cYQXgwK10eUni8/9AjBQV3vK3+kYFHKIqk4cAqkVRu38ccSff34cf4jk4qGcVPhuGEXBhC84XE24G4VLYy1iD3olyW4TI6CchYa2d+H13NPXbQGML5N6CmEEtSiCv86UWWz13a6Sbewu4xIijjnEJaZJux2jWe7Fo1Limt2CGOOctACwcTL2ZjXFntKkb8TUsI8HdIeP07OCyA0oLAhD9rwuLV55z+Sv4DJy8XAM4735DSP/EcM+1oEaeutMOXMSDkbYTcfVqXEb7n8CJ0fn7bVOQ/khmV5lSR2FY4b5WgC/j5+PuDwZVRidO98mW9KHtcoVwsEB87X7DX0feuL5Qy89sXAtvRR49XL8iOu1e9B9yV2ph94hOBS9l+avjQ/jMfI3+09xb0AKpjKyAhV6mY0xn9VVH1UypP+7rkZod22qyLmsEEEwvwAFlNRYM/n21/Gv6e2/cmzoGr2QwsEkwJdc+5//gd0rlSOYyi3o6a26OKvpLwfRslfXvsy73fQRwjSG0pjlasXQXJSqjcJ+wY5xo5bdjNUex6MX6Va0SxbX7hi349CK6P8dIMIhKRYY5t72xHCGKFclFK8EdU8iZ5DSIZcBXMuSDBxqjnOfvCjSNNfu2k7J935Hmc9MCFwz7u+cO87c/NPlDFCcBLUK6WzVx0HQdrgz5caHQm7YN//xte4MaB/IgjmLqLUHmnY24qtQ8iSQNjLthYh6PxXOVjf+NUZy8pexxM3gQSCiLQVkXEiMlNEZojI/iIyVkS+EpEp5t+xZthqEdlsO3+3LZ69RWSaiMwSkdvFVK6JSBMR+Zd5/kMRqY7jYX2fscBePf9/3fl4C9Al//g4t7/CMsew/V8fZXkf5/AEnUPw6v2//cXyAmd5QRrAqL7hjMXrmLF4XfGALrgJ+8Vrt/C3t+f43rdiwza22lx32N0juD1VCt4gXKnKkECwk4S7DGsOYYmLn7RSiSvXQUcItwEvKqUGAkMBqxtzq1JqmPln36Fjtu38BbbzdwHnAf3MP8sm7lxgtVKqL3ArcHOJz1Mydeak1n+V99si7t6l35zEz580XDbU1irWbgq/IClDWhUgjJWRO2c+MCHPHz8UbwCdk8rlcMxt73CMy1qQML3yUix68neAs6dbGLbUBi/0CAH/RmpbgqqZMNgdV7oRRVERx3/Irsq3qEAQkdbAKOB+AKXUNqXUmrAJiUhXoLVSarwyaszDwEnm5ROBv5vH44AxIsm+soJ1CLn/ydj/Wtht85Vy32/gjjdmFe1JupElPTuEmEMIojKqDT5CSPI9JDmJaC+rXhPPQQlb+4x0vBMa88e3wmdiJ2HmkvW8+fkyVtv8KPmNmIZ0L74PeprrEPoAy4EHRWSyiNwnIi3Ma5eIyCci8mvLTr4AACAASURBVICI2J2y7GqGfUtEDjbPdQfs/g8WmeesawsBlFI7gLVAhxKfqSzq3PmS998idgsV+wIk4KJHP867/u27x/PUlK9KinvrjlqueWoaKzeEsyBJmyCN6rhJi/jbW7OLziUYC9NKZ+uOGvr/8gWe9vkGQRaO5Yiw2/O3two7CWmpjLJuTZMkazdv5+wHP8pbJ+Q3p/KNod2SyJYrQQRCFTAcuEsptRewEbgKQ/2zGzAMWAz80Qy/GOhlhv0x8Lg5ynB7A1ap8buWQ0TOE5GJIjJx+fJoJ1wLVEZm8jOX5OuI4y7mFQ6JYN8CE2DCvFXMKeJF1ItnP/maRz9YkNvsJ22Cvssg4f70yhfc+MJMpi5cG13CLqzeuJ1tNbW+m9PYe+obtu2IZL2Kl2rJ3u7+ydbgWPMMNUpx/J5dw6cXeoSQPZVklvHbJz3ICv645jKDuK5YBCxSSn1o/h4HXKWUyhmCi8i9wLMASqmtwFbzeJKIzAb6m/H0sMXbA7BWDS0CegKLRKQKaAMU7KWnlLoHuAdgxIgRkb4R5wpkq3A7d3CKu+fj3F0pyvRyvcWMVNzAvuBD5Ldo41vmCMHyybN0nfcoy57fPce+THuf1amCIahveelz33S92ohi5aO2VtGlTbxeRHN5sb3Zv4+fz5tfLOcbQ7tx5ZEDEkk/DtIQcglry/MoOkJQSi0BFoqI9VXHAJ+ZcwIW3wQ+BRCRTiJSaR73wZg8nqOUWgysF5H9zPmBM4GnzfufAc4yj08BXlcJjzntK5Dt/50kOoeAYVUSOVmZ0Ar4LsOo6YLsxFZO0XLT/S5xTEw6o3f64Ldz9oMfccnjk5lX4qKsYiqyWqVCO9KDaLz6zl+5ib+8PguAxiXsdWB9pzWbtuXp3+Oid4fmsacRhCBvPm1fRpcCj4nIJxgqohuA35smpJ8AhwJXmGFHAZ+IyFSM0cQFSimrt38hcB8wC5gNWC4x7wc6iMgsDDXTVeU9VgkULEwz3rizIF/4WL5Ovz5SamX/7f+89yYuhaD67TCFv5g8eHfWirJGCG552e/G18qIsTyKWRHNW7mJfr983jeMG5GrjEoocsvNua5hv30lb+1AEF6bsTS04B/UtXAyd8rC0PYzZePWAeq/S0sg/sF9IG+nSqkpwAjH6e97hH0SeNLj2kRgiMv5LcCpQfISF7n9EByvPGlzuTSHi8V44L25xQOFIPj2gcGrQbGQFz/+MT3bu3vbjCJ+I0xyg9ti73BHTW0iE8vOlcpOSinVS9ZuoXOr0tRd5/59IjefvAff2acXEMz8ttxe95yIvBZ0aFmoYnz20oPzfqe9DmGnp3AdQkr5iDCuVk3c5X1W1ggFnlSO+FssXFX6xkKBFr4lWHaK9Ve2l+jNs5RnsN+zR/e6TYgWrNyUt5guKKWouuwsWl33nTcWcdi3S+smXHFEf9drQa3yvlgajUDo3rZQddU4xN4J5aAFgomz6AV1hRB5PiJsrLu3c+8JZ2UQEsekctxTT16xV1/1HJ8vWW/mIdYs5FGs57u+xB3Vwo5yDHPeunt62fTxP3+yNLfWYTY9csMabT/6wXyOv/1d37AvXT6KAR77W787a4XreSeJrjmJKS0tEEwKdkxLa4TgsRFKKVx8aN/ccZZ8yVgELdRhJpXT8KlvcdSf3048zWLvptT3UdJKZds9y2xuGpwOIYMyq8y9KjZt3cH37/+Qa576lAVFvBJHMYmeFTch5aAFgolVILbFuGFGmHyUy80n78Gxe9QZgtl7klnZFzpooxPmi6zdHO8ew0HynOQcQly+eMrtjHw0r7TtO+38dNwnvPrZ0uIBPZg4f3XBPgJeVFZ614kg72LL9proPJkGIG0ro50eq2d+u8NFc1r5gPIbFq9BQak9tqgJbmWUna5XkG+SZHbjWjkfXmXkHb4coeVcmBmGNs0aBQ7r50oiyLu49qlPGTdpUdFwWScjTUP62ItD9VXPZSIf5cUjPhZLWRkhRD+HEDuBRgjR4/XFSm1sR1a3970eNlpDZeR+UzkCodRV+QDvBdT9g/8E9rtfFndZPWdF6fl0Y95Nx0UaX1C0QLBIaKb16D+/7e/zPs+5XTRp9nBMLmdlOiHwCCErS6sJaHaakASbuWR96Y1tkTJQyjN43eJn6fSNod3o5rOS+smPS+91++110KdTi7zffiOEcvJQ39ACwSSpNnLmkvW+nkor8iaVy2hYzHgePXdf/nPRAY40siER6uUIIQBJZnfD1tL2Py6Gsy2tvuo5xj4z3TO8n9NAv9Xjh/TvxPtXj8k799OjEnB14chsRcK9JD8BFKR66jmEmMlKL8BeFqKYMDyoX8eCxT1ZcYMdNBdxW28sWh3MbcR/Jy/i0scnFw2X5Ov9TcSrxy3chPVD78/zu8Pzubf6OHJzI6xF3JbtNcxfGa3KJm7ev+ow+nZu6XrNr/zE7jon3ujrD/ZFLHHQ0WX1oRt2vX9cDWEpi4TCcttpw4qGCb4OId5acOStwcxFr/jXVCbMK/C5WEC5uX1y0qK8ndDSIGzqfp/oy5Dmo2F3V/vuvR9wyB/eDHVP2p2izq2b0s9DIAQhLjWqFggJMXpA50Dh7FWhnELrV6WCumP+T4mjptEDOnHisO5FwwWfQ4iXTdvKd08dJVf+e2rJe14U4/en7Am4++2xE1YgKaJrpMKOECYvCO9vKOqFp6UonLyy0K1tMt5p3dACIUJ+drS37jNogbHrD+PaADxIXVi2bgs/fmJq6Lj36tWWP546NGA+sjFCiJwI8ru6hC1SgzCyuj1PXrg/vzh2d/br421pFNrKSJX22G768rAjhFKmxAI4xY0drw5fj3bGKu8urb0Fg55DqAecsGf5Ox3ZF42VM4fg5yTPryGw2F5i2ueP2o0OLZsECrto9WbGPjM995zvz1rBtU99WhAuCXkwz8NscOK8VXy8IPgiq+01tUxdFGCTnpBE5fRQBPbu3Z7GVRW09PB1BcZk9YjfvRrYdLOYczu//DipDLlQppTyUaxuefkBi5JaBbt3bc1xHhsYjerfMfY8ONECISJuOXUoPds35/UrDwnU4Hpi3yAnLj1ygMblwJteLynqMJ27P7z0OQ+9P49J840G93v3fcgjH8wvCJfE+OCKJ6a4nj/l7vF86873A8fz34+jUfXENSqydzj8rM2+WLqeFRu2cvOLMwM79Cu2qHPXji18r1uEGSFMLdE9dTGVUdeQaptS5PXMJetYt3k7h/Tv5B6nj14hrjqhBUJENG9cCUCfTi1p3TT4CkkneVZGcXWNM66CcTZASUwAFhO+boLKjc0RbJcJhRP/kS1YtEUUxPx4e4239ZCdIGGG9mhTPBDh5hBO/Ot7gcPaWb7e34NpsTkWJ6UU0UWrN/PVms2e39bt88StPtUCISLy1g+UEY/ENEKIKn/FiCJuZ5kvtw48d9lBdCyixiomfN1UWW5ENe/zhyJbakaBX8NrqZNmLF7HpgBCLojCSER4+uIDXa+dP6pP7rjKx69QUtx08p6e16YuXJOrmxPmruLq/0wLHb+9YT9haDdO37dXQRg/ea29nWacYjrekia+yl+XlqMihhXQbkQRt7NxLjfKwd3acMBuHfzTjGiS8b+Ts7GexYu8cuhTJu2ywsvG/47X61REQb/7no5RglUurz5295zQjtsz71n79+blK0b5hmnaqNL1/MR5qzjxr+9x99uzAfj238bzjwkLQufBPpJs2qiS67+5h0uo5AWjFggREVWDa7+3nN6mUwDl588/3k+/KmdStHyJUDhCKD/OYo3MjMXreOuL5WWnc9we5RsWxEHb5oYaUyTYUNE+6eplkXPvO3ND5UGksONk/y7WpbBWRmHp0LIJ/Xdx3/ugGMtMVdO0Mg0HSt24yELPIWSc/K2XS/9c9jvnrijdH7zd9TXkezgtlrvj/+K/mYgfUYwQnOqHKOIM0sSc/8jEstP5YE5xR2hpMKxnWwCa2nbe8pubsbdXXuHs7fZXa4ov7HSbJLV3VKyjsFZGYQkrb4Z0r5tPaGK+v7CLO0XyJ9WDeGJNw8NM/LZVDYSozALtveGXppfuC9455E1MZRRFHM4RQgSxBvk+UbyXKEYZcXDH94bz6Vdr80yC/Z7XPn/lNb9iL1PnPzKpaB7cGuJe7et2VktqhBC2rrZrXudloEmVUa+cizuLFR2l8DXzdcMvl3odQsZxa3BLKddxfehKe/4iiM9LJx/JCMERRxSLiCobeElv2aSK/frkfzM/QbsjT2XkHq7cTtBrVx7CENvey9YIIu45hHKcO1qDl1Is38Im62plFDrVcDTwalIaHVoU+iWyl2GrkP/p23X+fFZv2s5nX68rGncUH/z6bw4pOGffESoKnbxXFFH05gtURmXHGKwRyLYxbvT8+AjvlfX2Bm/Fhm2uYcK2285PsFunlq7X4x4hWJ2DR8/dN1B4t7KulLFXc93v6EuP/86G8ZRWrTIqAbeekb3BuWxMP44YtEue2uaVz5bySoDtAKMoWPZhuFv+4iSKeuHskEbxToKpjBqWSPDaVB7yzV4vesxdHRS2SBVbwV43h1A84k+/WsuHc4s7GnTDqgsH9Qu2EtitkzN10Zq89Eux/2hUKXnP+vqVh7DEthe1xYF9O3DBI5M4oG8H1m+Jx925RYMXCCs2bOXdgPuuWrj1YOyVo7JCGNK9DbOXh58UjqJJchMI9hxnfw4h+knlNDcFOmZIF174dEl6GSgTr8YubCfj8sP7+V63hHaQdQjlGD6EVXXZy5/Va9/i49J7144tmBtgB7VpY4/K+92nU0v62EZNVjaPHNSFW17+nMY2gwA9hxATP3x4Ipf/y91tgRduS/Cj6oFHsa6pd4fC/EW5V7NfHPbG/MFz9uHUvXvkXb/z9OEB4g6WVhiCfJ64PIV42bSHIiGB5taZ8CJMme/WpmluQhaM3rEXWbMyCjJfYA8xeoC7KwonTRtVBiobSimE+OdWQAsElq4tHKIVo6vLln9ulaOkzxdS9Pdq35wpvzoiQMhorYyCxHHogM4cs0eXvHPH7tGVHx68q3/cjs5X2Pzu2aNNgcfVIHHEtmG9Ge8xQ+reRbGFcqVyzXG7l3V/XG2OvVf+yLkjef3K0Z5h455DCNt5sxeLIGpFL91/WRsgOvKs1yFkhCHdW/PrEwYXnHezYimlgQl7x9++vzdtbWZxrZu6awEDrkcKTDkDootG9837vWrjNrbblgo7RwRhe+4CnOwYmQQhNtdR5v+h5loAiM/GvJXH9w9KmMYyTEfeHvbgfp3o6abWNJOO38ooXHh7sfAsiy6Fp5i7lDAktSahwQuEsPrEZy89mGaNC4d5bvFs3hbeXjJso9TeYfE0rFc713CJzSE44rb3lr413H3TnAsenZSntiuYVA4pwqJaExIVbu/b34LELXzAcCGf/cRh+Surw+wtHEp4BAhrBSnn83m5ks7LSxkSIUgnL4rngPxv7oxKzyFkCLfy5PbtS9n5KHzjl//bS0cf9RyCF35xO80M7Tz3yeKc/rq5Q+CGLfx/OMXbMVkxtkfl1MiGlX17uYlLZlWGjLiLQ/0Z5u7IBYKZejmNXbvmxVcAh1UZTZi3iofHz+M3/5seTCB4nC+13imSsxJs8AIhyJJ7J269MDd3uh1aNmHeTccx76bjAsddruT3Wg0pCCeZvcE4RwjO9tReCazX5jbCAujbuSWDu7UumGgLaw7az8VPTdAYdpTpY8ZOq6ZVvPmT0blGxGrwdmkdnSrBSdj52KMHdykeyINwwiNAfFb5KGMSPsimUqVopH719HQefG+ez/qbQgT4+dEDATj/kD6h6529nXE2OXpP5QzhVqDKdVZlEeeeOFFG7VW4g9S1po0qmXTN4ezdO1+99frMZUz/eh39r3khP60S81gKUe5BsUvrplR3bJETaE1NQdivc6vY1Fphe5J79WrHyOq6DZ2CPv3GrTuYE8C00iLI81ohypnc/+/k4hsUlfPuvfKWZ5rqoTIq9bGM+/QIIbO4FajmTSIwLSR8bzioLjru4jS4m+EALGjuO7Rs4jl5uM3hOMz5Sg7sG4+FDkBNCYK9SZV7NbKezsp/++aNefTcfbnzjOGxfY9SJmTtvc2gjfEVIU21k1oH4iw7bgQVmjd9q9AltVeHbZqLh2BByurJW/k0vMTmX9NzCBlntMc2eFmiroGKYh2CwQ8O3JWJ1xzuu+q1FKqvei7nYtjZSA3uFmznLTtWFJcd1tc3XCkjBKePoFya5n/LlUn3ds04qF9HWjdtlJk5BAg3Kt20bQdKKb5Yut4zzB3f24ubT85vTINNKgfL+xyfBZ/F4mjZpIqmjYI1e6eNLNy0xtru1Q9L5WV1kiC/49bCQ2Xq5Edj+vHdkT05bZ9eCMFGP+XS4FcqR8GJw7pF6O00XPgmAQu3PX9+adz3zpxgCZtxHDl4F1/zOr8RTLE3dsId7zLvpuMK8juidzvucQn/l+/uxcAigqnY691Rgic9EejetpnnfNSFh+zGoQM6MyivgYiH0BY05Duya1RkEmLQr17ip0cN8G3gj9+zGwtXbcrPV4j6UewbHfbHtwLH5WTKr46gysfT4V692rJsnff2mne/NbtoGu1bNObfF+zP7l1b8/D4ebnz1nP96/z9A+W1TfNG3Pgtw0BimWOOUo8QMkyxor5uy/bAcYUZYv7i2IGB92/+as1mnprytZmGN797bkag+Kxeu1XRrzi8PyN3bc+Rg3cJdD8Es7SpqVUFI4R+u7TiBwcWLm4b1rOt64SynWIVyS4PNm8Ltj9yscauokLyhAEE7w1Pmr+KVz9bGjh8KdYodm+mbpvLj3KMfp+Z8nXRb1dK/+ii0bsB0KV1eOu8XLpFrhdTqf33ogN576rDSk4fjG+7T3V7Wjapcjc5zpZVdB56hBABxSprsQ297YSR/OeN2i144BLT8MKqWFb96tm+OU+49Hy8BNyc5Rv4YE5x52TL1ruvJHeLN8hz7dKmKWcfUM1D789zvW5XGV0TcB9loQTXxgHDnXzXeCDYhipQmpvvYjvztXWkHUSt5qwTny0u7un31BE9OXVEz6LhyiGJNSpuSYhE6zwxVSsjEWkrIuNEZKaIzBCR/UVkrIh8JSJTzL9jbeGvFpFZIvK5iBxlO7+3iEwzr90u5tcRkSYi8i/z/IciUh31g8ZJsaX2YXzZfLms9F3SghJFYfrTd4bxgwN3ZS+PhXDFCKoPfWT8/MBbagZ5rgqBsd8oXGluYZ9UDuqc0KjsgYLGjtsI4Z2fHep7j32EcMEhhZ0MpwCYFaCMZrgTnMet3xnKI+eOjDTOtz02SbL2kw4q3NMgaH/iNuBFpdRAYChg6RVuVUoNM/+eBxCRQcBpwGDgaOBOEbFaxLuA84B+5t/R5vlzgdVKqb7ArcDN5T1WsvziWH//Md3bNksoJ8GIovHq3rYZvzphUNEhuJunSAhuzTJj8bqCsG7OBZ1plYq98QtusZOd5s9NILi5ibBjzZu8csUo1wlytw1yio16y+2Il7pWI2y639yrBwf3i9YgxG17TQF+e+IQnr74QHq0C+5A0IvU5hBEpDUwCrjfyIjappRa43PLicA/lVJblVJzgVnASBHpCrRWSo1XRhfvYeAk2z1/N4/HAWMkibFdRLRz2TDHj7QFRBY6s8XUFBbTvlqXl9+rjhnoOeyP4rlqbJMIQS12KspQBwwocbN3LywhVmxy3Y71KbwmpN0WexX7fGFdczh5+uKDCs7NXLKOpS77BUSZbpw0bVSZ588qiwQZIfQBlgMPishkEblPRKwu2iUi8omIPCAilu6gO7DQdv8i81x389h5Pu8epdQOYC1Q0FURkfNEZKKITFy+PJt71wbBr50pZ5VmYDwar9te/TL+tE28tmV04lTHWMfuu1iVLxLsK62DrvoVcW8g/QSKdalDy3CdiSB5mf6bo3jusoMD39Ovs+FSxGuVu9uzFRvhldudc7rUADj6z++w7w2vlRdxAmRFfVgKQYp8FTAcuEsptRewEbgKQ/2zGzAMWAz80QzvVhSUz3m/e/JPKHWPUmqEUmpEp07Zt/v3wq/AlFqR7vn+3tzospDGLX6v5G999YvSEi+BoCMEIbh6yT9UsDjsZqdBVUZevVL/+41rcbjcbtGkKtQCtT99ZxiP/3BfdvGw7nETtHG5Ct9ZqS9vK4hAWAQsUkp9aP4eBwxXSi1VStUopWqBe4GRtvB2U4EewNfm+R4u5/PuEZEqoA1Q2v54CRO03o09YVDu2K8yldqxOnJwF77rspDGLX635Ddti2drPq9HDTNCKDetMNjNToOacBpCtjDxM/br7XnPbp2MQXbbZsVdl4ehFJVJyyZVHLCb93aSblZFWh54Yy8LljCNWpGV2hyCUmoJsFBErB25xwCfmXMCFt8ELBu9Z4DTTMuhXTEmjycopRYD60VkP3N+4Ezgads9Z5nHpwCvq3qywW3QRuPb+9TJyDR7VyLuy+mDLPm3ePXHo0pKe/eudXrtoCuC2zZrHFgV9KXP6lmrShaLyj5CCPptmzWuLFCrzLr+GL63r7eA/slRA3jonH04bGDn3Dm/rK3dHHwtS9S4zSEEcSKXBnF68i0Fq7zVlxnRoFZGlwKPicgnGCqiG4DfmyaknwCHAlcAKKWmA08AnwEvAhcrpawVPhcC92FMNM8GLC9m9wMdRGQW8GMMlVTs1NaqvH1KSyFoo2EfwqdZl7xGCGHm8Pt2Dj5haTX8g7q2ZvSAusYvqJfpPXq0Cfy+Fq7e5HPV7KkVeUy7sA46+vvTt4cVvNNi5aJRZQWjB3TO38gnwHN2c9Gt24mj4YmiA/OzowcUD+TgdycNCX1PFrqR/W2GAnU68Wg/TFyCL9AYVSk1BRjhOP19n/DXA9e7nJ8IFHxlpdQW4NQgeYmK5z5ZzMWPfxwobMsmVWzY6q5SCVoB7Q2EX6GNuzw75xBmL9/AmD++xV++u1cs6Y0e0ImTh/fgp0flNwiBVUZEW8mLVUy7++twjuLyMxm0XNjTiNLTKhjO2Sw/SuVQymigXfPGtG5aRe8OLZj21VrXbWeLccZ+vQMvDgxCp1bxuR23c/UxhWboO9sIYafj+U8XBw77z/P287wWWM9sO/ZTgQRtE/zy5J8PyUtjygLDgjgux1lNqir547eHFliNBJ1UrlXRTHY73+v/HVTo+gLyG+VQ/ncc8ZdiNV1sPQt4+N0X8lxYW5w2slckAiHIinInjasq+GTsUTx18YH85bt7cdIw993yosavVJXi+C8sQ3u0ydM6xLYta0zxatcVAfDrWQTtRNobCG+f6orN273957Rr3ojVm7Zzy6lDPT1sFs9I/nCzkVl4t/ikGwcH9evAkx8vKhrOa7I7yLzCfy86IG+/aTvXHD+INs0a8cdX8oVNKZPKEM3I7oz9enPqiB4MuOZF73SKmLdG0eR1bNmYFRu2RRCTMQI6YWi34gEjQinFJ2OPRCmYt2IjC1dv4pLHJ+fykhb1ZIDQcEcIYT6QX7sQ1LukPZhX57iYCuvYPYx5/HLm2wXyWi8rX0F77FFR3cF9tbETZ6M8tGdhj9ey2HGyV692uVXNQdv2UsxOobxvculhffn9yYZXy1J0zV5WTqWye9fWxQNlFAW0btqINs0aMbRnW47fs04YJaK2cSQSl64/rtracAWCS+k4a//erpud+PUU7/ie+x7Gful5jRCen7akSBzG/3IKg3MOwWqAgur0k8b+6idfe0TOPDLvGQLU9KDtdd6kcgiBYAnrUrjyyAE5KzT7owTd39m+W189WuAfC37fuZQRwreGh1R1eW4lWD++S4MVCG6IiLunQp97DillYxyPMhO4wJbRdhtzCOkvNAracNlz5eUixHptQR6hWLJ5k8oh6vAlRTbeCYo9yZtfmBki/X5UCJFsVBRIwJadSvKU4hr818d7O0IMQnxzCPFErAWCDRH3IXspBckPr8b34tH5niY7t2pCz/bNckN4Kx/lDENFYOGqzbw2Y2mguPrv0rLktKKgaVVxVx7WN4uiipQ6QoiqjNjjmbKw0GWY1/c6pH8n5tx4XCSeNMM+ycnDexQPlAFKmUKQMltIP1cMWUQLBBuC+wjBrbLfcupQ/nvRASWl46Wdae7wJXPC0G6887PDGNWvo5k/g1I6B73aN2f0gE4I8OL0JZz794l5cXnlyW+1bTkErSBBJp4topgztM+lJGGV4qRYkkkM5MI+9qj+3qucs0QpKqNyS8AZ+/Zizx5tON1nkWIpxFUMGqyVkduH9iovbr2EU/YuvVfkNUJo5NjdxBlOpPSe8NumT/whv34pd27S/NW5uLziLGf3Kj/cGp1dO7Zg7oqNge7Pc3hn5j6KXrrd5r5Uq5TbThtWcvrWNy7V/XMaZHXVspNSykfYORnnm+jcuinPXFLouTWr6BGCDUNlVEjUKiOvXl5jh9LaK1w5vUR7Cuc/MrFopElOUtrdOIShzj2Ad16DvjN741bK/sSdWzXhxDJt7vfr056l67YyMcCG7nEQ5KntOux6Ig948Jx9Qt+TVVWP3lM5AYxJZbc5hGjTCTpCsCqdtbjI8nEflSmbXT3iVantj3747sH3Sy6G21xNmB65810ZcQZJ15+aWsVf35jFnmNfivy7B+XjBd7bjWSl7bWXlzgt1DqE3GvEjZevGMXLV4yia5vw+5CE7QxmwXVGOTRYlZEbiY0QXM69NmNpwXnr9wlDuzG4W+ucv/qyCp19PUSeQPAaIdQd33eW03tJtIR5zZcf0Y8H3psLwOpNxiKqKBrwWqX4w0ufA8HmEO48PZjZcRj8VDDJzCGEe5FZcoV95v6Fc179y9iEKGzVj8sFTCHayihS3M1L3b9+VPJgzg3HculhfQsq0KT5qzj37xO54fkZeeftwfp0aplrzCOSByhVNwpJetgf5v270bppnTWNtWe1n9O9oKMq+6gpyB1eZ/QeJwAAIABJREFUm8qUQ9qOftubvfKBXVrRo13xXnXUPpjK4dLD+qWW9i+P3Z1qj+1d6wsNViC44dXDjMyksEKoEMOX0MG/fz13fvVGw7Xx+i35bhqcFS3XYJazUtn2LPb4vRqhJA1tSk3rnAN35dUfH8JB/cJZu7ilZx81BVJBOQJFY/oaQSRlMPYbg7nuxMG88KODA5X9tPNrJ+ryGqbuJ1lX9BxCAojg2gpE+Z2tArZw1eb8dF0odJhmni8jfXtaNUFURglOq4VN6Vt7GZO3jSuFvp3Dr5ewP/J1pqvlIPMqdqz3k9XJx1Jo2aSK7+9fjYiwYJWfS3GDOOcQwsYctXo3THQT5taLPb18abACwa0MexWmKAuZq6rEUyDkZ9JqtL9as9kteCDsUdYGsBSJy6Qw6HoPP1qbi7DC+GFyS6J722acOMzweVOTpzJyj/e7I+s2O8qCR4J+JQjDKMnSHELU6rYwnzeI8IwK7csoYv439euCc14fP8pKH2bi01nRtm43fNs8OzW4626/OGs9hIOdoOsCoiDse7askoIIrd6mM73OrQrXVZy0V7fcBHJNgBFC48oK9i/V22wM7NGjfBfX5ZCldQhR56Sh+YbSVkZ2PHagFxGevHB/olAMuBUwL7WMs55FsR2fvd23V+Q5y90b/ri8Nbo9c9jH6t7WmPBs5+Hi2s4Fh+zGkO5t8nxP1bkVkDrhYheYftY+1g5s5u9mjY2J7YNDzmOEpzBPSar13MjQAKHsHRCdpOgx2xe9H0IC+H38vXsXbkBSWhp1iWzYusN3c/uvVrurhsrx656V4b2rUAsp6c46oJoubZpyzJAuRcNWVoinI0IR22gjz3Ooe1wKw0PpRY99nOudt2raiLd/emjBRkBxMqJ3u9QWr9mxhGjXNk254vD+qebFbn0WBQ1thNBgVUZuJNHTsrflx972DiOvf81zyD1+zsq83z3bN+P/DtqVB84Ov+LSIisCwY2wcq6yQjh2j64lV1r7q8ipjGwna2oVVRXW6DCffarb89EvD6eVrQHq1aF55D1UPyyX2Wljlak3fjI6M3na2YnLNFmPEGxUeFgZRZtGXQLWJFTQRlpEuOb4QWWl76WG8iIu+ZGFjpdd7VNRYTg2tLu/rqk18jm8V7v8+1KVqS6qtpTfpaVai9rCR5M8eoQAfN/06ClC7L4BXG3fHWk+ffGBsaW/bUf+pivFhFGS84VJN7THDDE2tTnG3NymUoT3Zq/IXa+prU1dP1+I8ZLS3A7SiVVG4sxTozCbU8RMYxe3KUkTV1XRIwTqGsUk9IVuaTSuyj/Xq33z2PNhUawRjquOZ6GhHdClFfNuOi73u6JCcpZcAE9NMSzRsqZHfvbSg+jYsglvf7E87awAdcYJccooSaK3FoBR/TtxcN+OXO/wKrCzkL6oywBWD0cETo/J/7+FW6WprMj/DEm2P1f+e6rv9bjUAFGsQ4iaSpFA6ru4LK+CMqR7m7zJ67TFlVLK8AMW4/dzqzedWyXrIvy5yw7i4R+MTF1FB4br+jjQAoG6CRpB+E7Mk2Jujd5ZD0zI+52lHmmSWamqFH40ps4XTdJqggxpYTz5UcpWPG6MqG7PBYfsVjxgGbjVm+QcyeXnIUnDAS9ue/XLWOJN/8kygNUrrPDwdholQRqdDMkDhvZsG0u8Xlan9mcvd1+BsFRUCDOXrE80zTBcfczA3HwXpD9SsRjVvxM/P3pgrGm4CQSvPbbjzsNp+/Ti8sMLneglOQf2DXNlfdRogUC+yijuxjhI7z9D8oARvdsVD1QC7vtOSF7FT3qSOajKKi0rI2f2hvU0vs2Rg4uvw6jvFD572xTKh/G/cVUFl6c8Uotre1c9qYxRwZs3rnTddCVqgnzHtHXpdkSEpy8+kI/mReu4y22k5LT6rUpYh5Mlyx03nOViQJdWzLnh2JJ2dqsvKOVu0vqNod0SX1NTrDOXZLWN65trgQC0aFLJZ789GoAFK+N1UBWksc+CPBhZ3Z4JphAY2rNt5Kojt8bXOUJLuqHLervqOqrKeqYjwlleRLK3yDLJjlxcn73BCoSRu7Z3dVcb9zcN8iGzMEL4+TEDInPX4YbbMwr5W5gm3dYFbV/SaoYaSNvvivPRhbQXCBaS5Agzrjaiwc4h2Jd+J1mwsmRB5Efc+QziBtxybw1wy6lDY80PZMHK3Z8sdBSSxvomKzduK7yW+AfzTzDJ3dLikj0NViB4+Q+Kf4RQP1RG8VtbuamMJPfs7Vs0plXTugFs1wScxgX1D5NWz7QhjxCciEhmrKwsrj1+98TSiqvD1iAFwpSFa1ixobDHkQT1RWWUygiBumdv17xR3mrmDLyS1Kkvo8skEIGBXVpz+O67JJZmsY5AM3Nv7ySIqyg0SIFw0l/fy9vdyN7TSKMhLAgTaw6CEXdv1M11hX0diH20YFyL/61kq79ZSBY6CllBMMw/7ztrRNpZyZVTPYewk/D+7Do302moSkoJEzdx+xpyNTu1CQFxhEninQTfGzgd0dEQVUZevfI0RvheX90qmzuDlZEWCOTvFpaJhWkZqPix58FjUtkSRBUiiauM1m3x3qwoC2Sho5A0XvM642evdD0fJ17CyWutRJwc3M99s6dyCSQQRKStiIwTkZkiMkNE9rdd+4mIKBHpaP6uFpHNIjLF/LvbFnZvEZkmIrNE5HYxW0cRaSIi/zLPfygi1dE+pj8dW9YtgU+jZ+wkC7ri2AWjxxaauRGCY01C+m8kA+iXkGNHbW3xQBHjNYltnU1qBPf5745mlMfuf+USdIRwG/CiUmogMBSYASAiPYEjgAWO8LOVUsPMvwts5+8CzgP6mX9Hm+fPBVYrpfoCtwI3l/IwYWljmjX236VVEskB9aeXl4Z7aqcgtP9eu3l70tnxJD0ro/pRdqJk3z4dXM/by0a3hLYt7dzKPZ19qo31Okl9nyZV8U1eFxUIItIaGAXcD6CU2qaUWmNevhX4GQGUqiLSFWitlBqvjDHWw8BJ5uUTgb+bx+OAMRJTN3n9lrqGxfJaaF/xaE+1S+umkbtPqC9VuiJmZWKHFo05anC+hYhIXaUSyRdJXmbCUfKgz9ak54/qkzvWZqd1vPrjQ2KN//bT9uLFyw8uOG9/F+9fPYbRAzrx3ZHxeipu7+FM7/6zRvDcZQftFKvGg6xU7gMsBx4UkaHAJOBHwBjgK6XUVJe2e1cRmQysA65RSr0DdAcW2cIsMs9h/l8IoJTaISJrgQ7ACiLmp//+JHfcyPyA9rbG/iTv/vzQyKcPg8q5u88YnqpOO3bVWYXwt++PoPqq5/LStE8q219VEju3HTqwc/yJlEEWRwh9O7eMNf5mjSsZ2KV1wXln+XzonJGx5sOPVk0bMbhbm9TSj5IgAqEKGA5cqpT6UERuA8ZijBqOdAm/GOillFopInsDT4nIYNw7x1Y197uWQ0TOw1A50atXrwBZL8Rubvr12i1GQspdIlTF4OwuaCfiaHN7x7RIo+0Rm9lpRUW2GkCVd5zOECFDryN9EnwXw3u1jW0SNyx3nT481viDtHiLgEVKqQ/N3+MwBMSuwFQRmQf0AD4WkS5Kqa1KqZUASqlJwGygvxlPD1u8PYCvbWn0BBCRKqANUOBoSCl1j1JqhFJqRKdO0X2gfHkQ96Ry/ajVaYx+K6RuBGUfLWSBoKuY4ySNsnPVMfHuc1AqSb6J/1x0IFccUeju+uaT90j8/Vj7f8dFUYGglFoCLBSRAeapMcDHSqnOSqlqpVQ1RoM+XCm1REQ6iUglgIj0wZg8nqOUWgysF5H9zPmBM4GnzTifAc4yj08BXlcJ1sAkvSbGrZuPjhQmlW1CYNpXazM132IvIg1pUrmtzZ9UlshCZ+E7+/SKfae4pAnq7fRS4DERaQzMAc7xCTsK+K2I7ABqgAuUUlZv/0LgIaAZ8IL5B8aE9SMiMgtjZHBamIcIg1tdzptDyMA6hCyQyvyYON5P3rtKv4eeNmkUHb+3PrBLctZ5mmQIJBCUUlMAzzXi5ijBOn4SeNIj3ERgiMv5LcCpQfJSLvY1BxadbJt1p7AeK5OkIbgEmL1sQ+63XSilrcPNgjh65bOlHBuzysBJBjRlrtQX1Wt9o94oMKLi1BH5pml/OGVP/nDKnrnfcTeE9aUgp5HLygph+YattjzU5aJFk3S37shTGSWYrn2v4hW2d9PQqSfVKDIaJ7CbIzTADXLs5eicA6sLBEQWfBllgTTy6Uwyq68qyV6zfZSUxg5hWXMx3VB596pD2ZCAGXrDEwi2CrZXr3g2kPejvqxdSaMxrnAsRsvSu8r3iJtOHi4e3TfxNLOqMkpjJX2adG7VlM4JTNk0OJWRvSC5e9yMOf2sdnszgIjkjUyyVOmTbBi/t6/7Ghv7XFdDR1ejeGhwAqGYS+UsOLfLAmksw68Qx/vJ6LuKWzjc8M09OM90lWFPKmtWRpqdjwYnEPI3XXELEHf6GW3lHKSVS8kbIWg07tSXelTfaIACoa4gvfLZMpfr8aafpRGCX17SmFR2OrCrLxPwcWCty0z9DWR1EkETCw1PINiOZy1bn3z6KTdy3x5R5z3ELy9JZbNd80aM6t+JB84ewTFDuuZ9nwYsD3LtcP47SP6FZFUcZKljtTPRAK2M6krS6AGF3i2bxuhrHNIvyL8/ZShVlRU8/uECKsRYSu5GUtmc/Cs3/4hmHjJQ6b9jmiVnwZdR1prnXx63e9pZ0ERMgxMI9gb5iEG7FFxvXFWR14uOPv0MtHImhnB0b2TSHslANt5V59aFlj1J2OZ3MTd96djSnn4KIwSfR0179bgmehqcQMifVHavYL8/ZWgi6Vvs1qkFs237OieFX/OSgba4QXPOgbvStU0zjt2jCz9+Ympq+cjGyKgQXTzjoeHNIdhaujQ8j7oJoe01ipcuH5V6Xi47rK/ntcTI823XcKt9ZYVw3J5dG/Q78EO/l3hocALB3tBVZsA9AxiVP8msWJ0+53xG08Z18ydZqG5pz7dA3Xuw95NbpuZXKZu9dc3OQ4MTCHmuEVJZfFWYZu8OzVNpgJ15mTRvde44vQFCtlYqW02wJUS7tmma53AuCSwPvWmM2rIqgrLQWdgZadBzCGmMENwqdU2tSm3LSjszl6y3XUunxmXBZ5Adpwr9wtG7Je559bnLDuaZKV+za8cWiaabbTJQOHZCGp5AsBWkyhS6Gas2bis4Z3ixTEN9JY7f7sdpkYEs5Hqi/XcxNpPv2a554nnYpXVTfmi6skiajM4pa2Ki4QkEu5VRCgJhW01twTml0vIumv/bLiCz0BhnYuLQzMMZ+/Vmjx5tGdazbcoZ0kA2Oiw7Iw1uDsHuUz4NlVGHFoU7tlVWpKMtd6qv7O8jLSujvDmEDFV6EWmQwiCrA4QMFY2digYnEOz+ctKYmBrSvQ0H9u2Qd+6HB/dJZ8tKR5r2EVNajXHaG7Kctk9PjnRZsNhQyew6BC0RYqHBCYS9etZtipOGyghgaA+jpzmke2vm3XQco/p3SrjHY1Ry5+Mf1Ldj7jgLq4TT4KaT9+SeM+u2D2+Yb0HTUGlwAqFN80a54zRURlA3DD9mSN2G6WlkxTmpfv4h6UxcepGFzmkDlYu+XHpYX56++MBU85AFk+SdkQYnEOykNUJwI40C7hwFZEF/n5+n9CTCqXvH58+qPuEmlK88cgBDU55P0YI6Hhq0QEjD7BTcK1kW1iHkXUupB3bGfr1TSdeJ5VxO90SziRYI8dCwBULKi6/SLtR+8wRp5W1gF2Mn8caV2SiaaX+jtEl7kl+TLNmodSmRhnM7IKcJcVPR7OLibjkuNm3bkffbXvmzMKmc5hxCFuYvssQPD9417SzkoUdu8dCwBULKk8r5K4ONH0kW9BUbCldNW6Re3SQbNvCpvweNO/rDxEKDFghpqYws3FJPIktW7/f4PbsyvFfbnPO0LEwqZwWtKjHI6kipgRfP2GjQAiEtK6OsLPY5sG9H/nPRga6T62m7jRDSbYwsk+AjBjfsRWp1o1nhuD26+oZNkrTL585KgxYIaeG+gbp5LtmsAHWNX4sm8e4nHQRLOB3cr2ORkPEypHsb5t10HAO7tE41H1lBgNtOG8anvzkq7axoYkQLhBRJa2JseG9jtXa/zoYHz2uPH8TH1x5Bq6aN/G5LhEaVFbzxk9Hc8b3hObXN3r3bFblLkwRVlRUpbg6Ujx4fxEM2vm4Dw08TkoSW5NS9e3DAbh3oYbpyrqwQ2rs43UsLp9//Lq2bppQTTVbRGqN40COEFMmzMko0XckJgyyTkamWBo3+Bg0LLRBSQFeyYORek+4NpoaqWzSTKTKWnZ0GLRBSQJs0hkNXfo0TbWUUD1ogpIgu1P5kxTxXkz10zYkHLRBSIMvtXFasSDTZQLm4WdHsvAQSCCLSVkTGichMEZkhIvvbrv1ERJSIdLSdu1pEZonI5yJylO383iIyzbx2u5hdZBFpIiL/Ms9/KCLV0T1iIW7bWKZBFqvYCz86mLvP2DvtbGgyRtqD2Sm/OoKPfnl47nenBH1+NSSCdgdvA15USp0iIo2B5gAi0hM4AlhgBRSRQcBpwGCgG/CqiPRXStUAdwHnAR8AzwNHAy8A5wKrlVJ9ReQ04GbgOxE8nyuv/PgQVm3cGlf0gUm7krnRs31zerbPlgWSVq1p2jbP78RdcXj/lHKyc1N0hCAirYFRwP0ASqltSqk15uVbgZ+Rbz5/IvBPpdRWpdRcYBYwUkS6Aq2VUuOVoRx+GDjJds/fzeNxwBiJsRVo36IxfTu3iiv6omjdeDBUNg1cGiRZ+wZZcY++sxHkrfYBlgMPishkEblPRFqIyDeAr5RSUx3huwMLbb8Xmee6m8fO83n3KKV2AGuB/J3oARE5T0QmisjE5cuXB8h6Nsn5h7Gda9rIcBsxoEt6giprZGXfCE320GUiHoKojKqA4cClSqkPReQ2YCzGqOFIl/Bun0r5nPe7J/+EUvcA9wCMGDGi3naz63wZ1T12+xaN+ccP92NId+07x4mu++mR1dGsViPGQ5ARwiJgkVLqQ/P3OAwBsSswVUTmAT2Aj0Wkixm+p+3+HsDX5vkeLuex3yMiVUAbYFUJz1MvGNjVGAVUO1w07L9bh0z4E8oKGW2LGhR+jhg1Ox9FBYJSagmwUEQGmKfGAB8rpTorpaqVUtUYDfpwM+wzwGmm5dCuQD9gglJqMbBeRPYz5wfOBJ4243wGOMs8PgV4XWW1axIB3xvZi2cvPYhD+ndKOyuZZs8ebQA4YlCXlHOi0WanDYOgVkaXAo+ZFkZzgHO8AiqlpovIE8BnwA7gYtPCCOBC4CGgGYZ10Qvm+fuBR0RkFsbI4LSQz1GvEBGGdG+TdjYyT9/Orfj8d0fTpCp9t9waTUMgkEBQSk0BRvhcr3b8vh643iXcRGCIy/ktwKlB8qJpWGhhoNEkh7bd0mg0nuy0eluNK1ogaDQaT7I2qZz2Tno7O9pxjUajKUpG5AH3njmCNZu2p52NnRYtEDQaTb2haaNKurTR80pxoVVGGo3GE713R8NCCwSNRuPJ4G6GefRgbSbdINAqI41G48kRg3bh7Z8eSq8O2fKAq4kHPULQaDS+aGHQcNACQaPRaDSAFggajUajMdECQaPRaDSAFggajUajMdECQaPRaDSAFggajUajMdECQaPRaDQASH3dmExElgPzS7y9I7AiwuyUis5HPjof2coD6Hw42Rny0Vsp5bpdY70VCOUgIhOVUp4b/uh86HyknY8s5EHno+HlQ6uMNBqNRgNogaDRaDQak4YqEO5JOwMmOh/56HzUkYU8gM6Hk506Hw1yDkGj0Wg0hTTUEYJGo9FoHGiBoNFoNBpAC4REEJFM7FGeoXzocmciIqlvUiUircz/mSgfWSAr7yLpfOx0FVNEupr/U92JW0R2F5H9AVSKEzUiMkREjhKRqpTzsYeIXAmglKpNMR8jReSGtIWSiOwvIvcC+6SYh+EiMg44F9IrpyIyTER+KCJd0kjflo9BIjIKUq+zg0VkdBr5SL13EhUi0hK4CzhdRIYqpaaJSKVSqibhfLQBbgFGAstF5EPgQaXUrITz0Q64HjgAmA0cLiJ3K6VmJ5kPG9cDR4nIJKXUm0l/GxFpDdyI0QA/pJSqFRFJo+KLyA+By4A7gckpvIsOwFiMd9Ee+MA8n3Q+GgF3ACOAGcB+InKPUurDpPLgyMd+wBcisi/wulJqkohUJNWBMTspdwCHAQtEZAzwtFJqYlL52JlGCCcAC4E/YwgGkhYGJj/DsN4aCpwPdACqU8rHVqXUMOD/gMFA4sNgm0rkbeA24HdgfJuEe+m/xKjwRyql7jTzkFYvsBfwS6XUXUqpLSmU01swHn8/jNHB9yGV+rIH0EYptbdS6gyM9igNtxBDzHwMBS4AtgNXiEjzhEez7YBWwO7A6cBK4EoRaZlUPuq1QBCRI0TkSPPnc8CtSqkfA71E5DQzTOyjIEc+7gJ+BWD2xttiFPzYMfNxlPlzrFLqCvP4SIye4GBLX5xAPg4HUErtMPWgRwH3AstE5P/Ma7Vx6kgd7+MBYBnQWUROEZFbROQ0EekVV/qOfBxuHrfBEM4TROQwEXlJRH4hIt8yr8fyPhxl9AKl1GXm8XLgMxEZEEe6Hvk4wvxZA3xbRNqYz///7Z1/rJV1HcdfH1AUhsytJEDAH8DFHzfw6kwocTnURTa5ptJmGrqWTixzDRR1mjA00tj8kVpOstKVi6g2W2xtZa3EQisIdLLWNGy6VYDAhrIxPv7x+Rzvs7suO/ee5zznQd6v7ew+v849r/M53+f5PN/P82sWMNfMenLZdreNhscI4IzsIW0H3gVOI0tpbfa43MwW5egYYDYwyt3/C6wFdgA3ttujwSGZELLG9gxwOxEw3H03fXsXXwPuy+n7K/C4o+Dxb3d/s5CI3iFKNm2jXzy2p8e+nHcesBj4AXApcJeZTazA4+2c1jh2sZHowa0AlpjZGjOb2I699AHisZUojawDFgFbgSvSpbJ4uPsuYBfwNNBLlI3eIn6XmWXHY4A2uq/QO3NgArA3l29XQirGYmd6bAJWEjH4DnAvMAlYbmZdFbSNnTl5K/AC8IiZnUxslH8OnGlmH26Tx2gzW0usmztzPXkNeB64ORd7C/gZ0GNmEyrp0br7IfUi9nS3A48OML9xsd16YFkOH90BjyPy76+BnhweVrVHv2VPIxLD+VV6AKOIDeApwCpi4/j7nDe83R6NuAMjgWv6xeNJYE4V8Sh4TCb2zJcW5t0L3NWJtgH8FvhyDlvFsTDgTuC8HP8Q8BBwbcVtdBLwMPAscXxnJvDDsttn4fNO7RePxnbrfOBXwMk53g18H5jSDo/+r0Omh2Bm47NLt4PYqByV068xs4syswM0zi7qBW4ys7uBB83sI1V6eJRKpgI73P1vZnYDcKeZHVulhyXp9Apx29zXy3Bo0mOqu+/NxV8ERhMHzSab2QwvqW59MA/igPo0d3+HSIjA+/EYB2wrw6FJj6nuvo1YyS8vvHUssRPTVod+baOx/q8hSmnDPbdC7fag7zdxYD+wAMCjZHM88EoVHhmPae7+BpEILnP3h4B/EMlpZIke46zv7MeZwMScvgj4upmdS3zv9cRxHtx9C3ACsK8sj4NSRdZpMZPOBf5A1NO+m9NGApvo61I9APwdmFp43+lEjfJ3wEcr9piW8y8kykXPEaWK6VXHg9gDOwKYD/yGqKUfQ4t7gk16PEiUik4BLga6Cu+/GpjcwfZxScZjdYXxeADYApyY89cSJZM/5fzxVbfRXGYFcSZc1evs5myjp+d6cn++7xnguKo9cv7wbBsvAcsooVffz+PxnNZFlMm+BzxFlInWAdcRPeo/EmccbSF2Hlpuo025tvsDWgxkF/BnYk9qbAZsbs6bBywsLLsaWJHDEzPYn+uQxzdy+CqibntBh+NxAbAB6O2Ax5PA8sL4sDJWshbj8XHgLx2Kx2rgmzk8hkiWF3XA4Z7C+EnAvA61jcZv0kOclXdph9vGdGLD/dk2enyCKBvel+3wyFz2avoSxthsp5eU4dG0b5Uf1mQA399gEN3Ih3N4DPBTYmN/ZHH5/HsZcXCqlCzaosdjxWmdjEddfpe6eNSkfZQSjzo4lODxWF08Kmwbx+f4J4njN1fm+AzgF5S0wzQk90598ACBvBZ4k9x7Ifae1hOnK/6LOJf9aeCpnN84ELOQ6AbPr4lHWXueHxQP/S4lx6MODvIYssePgCdy/nyil3ArcfxgMZFQ2l4e+r/+nfjQAQI5msiOXwX+StbbgeOAJcANOX40cYbGrBxeRRwnOFse8jgcPOrgII+WPf4HzM7xc4hy2ewyPFr6Dp0W6BfQyfl3JfDjHB5G1PjmFJb7NvCpnHeCPORxuHnUwUEeLXt8puzPb9m/0wIDBHUccRB0Xo7fRtzrZDpxQclG8jxdecjjcPaog4M8huxxUrs9Bvuq7RPTzOx64Cp3n5Pj3wLGE9n2Fo/zhuUhj8Peow4O8qivx2CoZUKwvLOfxa15/0NcVv8TYLPHxUXykIc8auIgj/p6DJZaXqmcgRxFnIu7ANjm7huqDqQ85FF3jzo4yKO+HoOlzs9DWEQcqb/Q80Zt8pCHPGrrII/6ejRNLUtG0Nflkoc85FF/B3nU12Mw1DYhCCGEqJZaHkMQQghRPUoIQgghACUEIYQQiRKCEEIIQAlBiKYws2PzyVaY2YS84EiIDxQ6y0iIJjCzE4Ffunt3h1WEaBt1vjBNiDqxEphiZhuJ5+2e6u7d+XzgXuLRi93ELZVHEE+/2gd82t13mNkU4BHiVsh7gS+5+6vVfw0hBkYlIyGaYynwT3c/g7i3fZFu4ErgY8A9wF537wFeAL6QyzwOfMXdzyIegvJoJdZCDAL1EIRonefcfQ+wx8x2Ac/m9M3ADDMbTTwfd42ZNd5zVPWaQhwcJQQhWqd4n5oDhfECn+GUAAAAgklEQVQDxDo2DHg7exdC1BaVjIRojj3AMUN5o7vvBl4zsysALJhZppwQZaCEIEQTuPt24Hkz2wLcP4R/8Xngi2a2CXiZeLi6ELVCp50KIYQA1EMQQgiRKCEIIYQAlBCEEEIkSghCCCEAJQQhhBCJEoIQQghACUEIIUSihCCEEAKA9wC+0wI5Hr8ANgAAAABJRU5ErkJggg==\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"zarr_store = zarr.storage.ABSStore('cssp-china', prefix='zarr_hourly_1851-1859', account_name=account, account_key=key)\n",
"ds = xr.open_zarr(zarr_store)\n",
"ds.surface_air_pressure[:,100,100].rolling(time=24, center=False).mean().plot()"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3.7",
"language": "python",
"name": "py3.7"
},
"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.7.6"
}
},
"nbformat": 4,
"nbformat_minor": 4
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment