{ "cells": [ { "cell_type": "code", "execution_count": 4, "id": "f031be76-56ab-48f7-9892-afdaccbde66a", "metadata": {}, "outputs": [], "source": [ "import xarray as xr\n", "import salem" ] }, { "cell_type": "code", "execution_count": 5, "id": "5908d766-1b62-4da1-8835-7bb8c7103977", "metadata": {}, "outputs": [], "source": [ "ds = xr.open_dataset('geology-optimizedC1DF.nc')" ] }, { "cell_type": "code", "execution_count": 9, "id": "efffe4bc-8ad9-430f-8a75-93d41b8b1461", "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "/home/users/fmaussion/.miniconda3/envs/oggm_env/lib/python3.10/site-packages/salem/graphics.py:745: ShapelyDeprecationWarning: The 'type' attribute is deprecated, and will be removed in the future. You can use the 'geom_type' attribute instead.\n", " geomtype = shape.iloc[0].geometry.type\n", "/home/users/fmaussion/.miniconda3/envs/oggm_env/lib/python3.10/site-packages/salem/graphics.py:749: ShapelyDeprecationWarning: The 'type' attribute is deprecated, and will be removed in the future. You can use the 'geom_type' attribute instead.\n", " if 'Multi' in g.type:\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAASQAAAGxCAYAAADLUsscAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8o6BhiAAAACXBIWXMAAA9hAAAPYQGoP6dpAAB56ElEQVR4nO29d3xc5ZW4/9wyVRr1LktyrzLYsQlgmgFjQktiCJDgULIkyy4kgUASSsJiko2NyS4ha35kv7AsJQktBUM2JNgkGAKmuGDjAq6yLcvqbTSafu/9/THSWGONrDbSzNx7n3zuJ+idO++c1zNz5rznvOccQdM0DRMTE5MUQEy2ACYmJia9mArJxMQkZTAVkomJScpgKiQTE5OUwVRIJiYmKYOpkExMTFIGUyGZmJikDKZCMjExSRlMhWRiYpIymArJ4GzcuJEVK1bQ0dERMz5x4kQuv/zyQZ//zDPPIAgCmzdvHiMJTYyEqZAMzsaNG3nwwQf7KSQTk2RgKiQTE5OUwVRIBmbFihX84Ac/AGDSpEkIgoAgCGzYsCF6z1//+lc+97nP4XA4mDlzJv/7v/876Lz19fUsWLCAadOmsW/fvrES30SHyMkWwCR5fPOb36StrY01a9bwxz/+kdLSUgBmz54NwPbt27nrrru45557KC4u5n/+53+4+eabmTp1Kueee27cOXfu3Mmll17KhAkTeP/99ykoKBi39ZikP6ZCMjATJkygsrISgPnz5zNx4sSYx1taWnjvvfei95x77rn87W9/4/nnn4+rkN58802uuuoqli5dyq9//WvsdvuYr8FEX5hbNpMBmTdvXlQZAdjtdqZPn87hw4f73fvss89y6aWX8s1vfpOXX37ZVEYmI8JUSCYDkp+f32/MZrPh8/n6jb/44os4HA6++c1vIgjCeIhnokNMhWSSEH77298yc+ZMzjvvPLZt25ZscUzSFFMhGRybzQYQ1+oZDnl5ebz55pvMmjWL888/nw8++CAR4pkYDFMhGZy5c+cC8Mtf/pL333+fzZs309XVNaK5XC4Xf/3rXznjjDO46KKLeOuttxIpqokBMBWSwVm8eDH33nsvf/rTnzj77LM57bTT2LJly4jnczgcvPrqq1x88cVceumlvP766wmU1kTvCGbXERMTk1TBtJBMTExSBlMhmZiYpAymQjIxMUkZTIVkYmKSMpgKycTEJGUwFZKJiUnKYLhsf1VVOXbsGC6Xy8y5MhkXNE2jq6uLsrIyRDExNoDf7ycYDA74uNVqTcsEZ8MppGPHjlFRUZFsMUwMSG1tLRMmTBj1PH6/n0lVmTQ0KQPeU1JSQk1NTdopJcMpJJfLBcDZXIqMpd/jc78xiaPvNtO+zzPeoo0rs79eRdPWdlp2u5Mtypgy66uVtOzupPmTzqTJECbEu7we/eyNlmAwSEOTQs2WKrJc/S0ud5fKpAWHCQaDaaeQDOdD6t2myViQhf6XI8OBpMZ/TE+Xw+lA0uSkyzH267QjaUmWo+eHL9EugiyXOOA1HN555x2uuOIKysrKEASBtWvXDnjvLbfcgiAIPProozHjgUCA73znOxQUFJCRkcEXv/hFjh49Ouw1GU4hDUbb3i4CnQPvzfVC+/4u/O2hZIsx5rQf9OBv0+f7qWjqgNdw6O7u5tRTT+Wxxx476X1r167lww8/pKysrN9jd9xxB6+88govvvgi7777Lh6Ph8svvxxFGXhbGQ/DbdkG49D6hmSLMC4ceasp2SKMC7VvNydbhDFDRUOlfypqvLGTcckll3DJJZec9J66ujq+/e1v88Ybb3DZZZfFPNbZ2clTTz3Fr3/9a5YsWQLAb37zGyoqKnjzzTe5+OKLhyyLaSGdwORLSrFkSMkWY8yZeFExtuz+PjS9UXVBMfZca7LFGBPUk/wPwO12x1yBQGBkr6OqXH/99fzgBz9gzpw5/R7fsmULoVCIpUuXRsfKysqorq5m48aNw3otUyGdQOlpeYiy/v9ZSj6Xh2TT/zqL5uUgO/X/AxOPiooKsrOzo9eqVatGNM/q1auRZZnvfve7cR9vaGjAarWSm5sbM15cXExDw/B2HOaWrQ+CFOlLpoaGtwdPR0RZQA3rv/KMKAtoOl2nomkocaoH9Y7V1taSlZUVHe+tDjoctmzZwi9/+Uu2bt06bKe8pmnDfo7+fyKHgWiJ/OMpIX1+gKMIIMqiQRSviBrW5zp7fUjxLoCsrKyYayQK6R//+AdNTU1UVlYiyzKyLHP48GHuuuuuaNuskpISgsEg7e3tMc9tamqiuLh4WK9nKqQ+aAoceadJtx/gXgQhss6wf3gRkHSkbmMLwe5wssVIW66//no++eQTtm3bFr3Kysr4wQ9+wBtvvAHAggULsFgsrF+/Pvq8+vp6du7cyaJFi4b1euaWrQ9qSDVE9ElTjRNlO7JBv+tU0VASEGXzeDzs378/+ndNTQ3btm0jLy+PysrKfu2wLBYLJSUlzJgxA4Ds7Gxuvvlm7rrrLvLz88nLy+P73/8+c+fOjUbdhoppIfXBnmtlxlX6TyuxOCVmf60q2WKMOaJFZPZ1VYiyPnMWB9uyDZXNmzczf/585s+fD8Cdd97J/Pnz+bd/+7chz/GLX/yCL3/5y1xzzTWcddZZOJ1O/vSnPyFJwwsomBZSHyxOiawqZ7LFGHMku0TO5IxkizHmSFaRvGkuhnlO0HAsXryY4ZTWP3ToUL8xu93OmjVrWLNmzahkMRVSH0SLMRy9kkVE1bvjnp4Im6qhqfpc62BRtnTEVEh9EC2CYb6ois4d99DzA6PTkD+A2nPFG09XTB9SH7xNAV07QXvxd4Q4tE7/KTIhT5gDrx9Lthgmw8C0kPoQ6AwR6NR/wmmoO0zzzuSV4xgvwn6Fpu0dyRZjzFAGiLLFG0sXTAupDwVzspl8aWmyxRhzciZnMv3K0RcKS3Vc5Q5mf60y2WKMGYo28JWumBZSH6wuGWuG/v9JLBkStiz9J9bKDglbjj4Ta8H0Ieke0SjRJ4NEE42SHqMn9G8ODAPJYpDok4ESa/W8ThUBhf6HPtU4Y+mCqZD60LitA1FK3zdzqLTt6cJ9xJtsMcaczsPd+Nv1WS0SQNUiV7zxdMVUSH3Qa6nTEzFKNDHYFSbYZSbWphOmD6kPUy4vo/zM/MFvTHMqFxcx8aLhlYVIR0o/n8e0L5YnW4wxQ+nZssW70hXTQuqDxS7hT7YQ44BsF1H1X3kEySYh6PgndyDlk84KScdv1/ARLYL+i7Oh76JlfdG7U1uPmBZSHwwVDjfAF1XvilfVBFQtTpQtzli6MCoLadWqVQiCwB133BEdu+mmmxAEIeY644wzBp2ro6OD2267jdLSUux2O7NmzeL111/vN+9DDz0U87y1a9cmrAHfvlfraN2j706uAIf+1kDjlrZkizHm1G1s4ei7LckWY8wwfUh92LRpE0888QSnnHJKv8e+8IUv8PTTT0f/tlpPflo2GAxy0UUXUVRUxO9//3smTJhAbW1tv9bDdrud1atXc8stt/TrcJAIwj4FxQAWUtivoin6X6cSUMxaSGnGiCwkj8fD8uXLefLJJ+MqBpvNRklJSfTKy8s76Xz/+7//S1tbG2vXruWss86iqqqKs88+m1NPPTXmviVLllBSUjLidi6DMe+fp5A9Uf+Fy2Z/rZLCuTnJFmPMmXpFua6jpgrigFe6MiLJb7vtNi677LIB6+Vu2LCBoqIipk+fzre+9S2amk5e0uO1117jzDPP5LbbbqO4uJjq6mpWrlzZrw2vJEmsXLmSNWvWDLlveCAQ6NcwD2DO9ROpvn4is79aiT0vYsFZXTLF83OpPL+ICWcXIFpERItI4Sk55M/KIndqJtlVEYUlygL2PCtWl4xslxDS6ECl3n0rvejdqa31+JBOvLQ09iENe8v24osvsnXrVjZt2hT38UsuuYSrr76aqqoqampquP/++7ngggvYsmXLgG1YDh48yN///neWL1/O66+/zr59+7jtttsIh8P96vouW7aMefPm8cADD/DUU08NKu+qVat48MEH+403bG7FZrMjWUUUvwICkXYcgDVTRpRFBAEkm0jJ/FxEi4BoiYxtfXw/mWUOqm+YiChFdLqvLcCWNfsomJPNlEtLUUIqakjDXetl/2t1FMzJpmhuNkpYQw2puA9307itg7zpLjKK7ahhFSWk4anz4qn3k1nmiITnQxpKSCXQESLsV5DtEpoaGRtplQnDFKKz6Fvx6jHsPyyFVFtby+233866deuw2+1x77n22muj/11dXc3ChQupqqriz3/+M1deeWXc56iqSlFREU888QSSJLFgwQKOHTvGz3/+87iFxlevXs0FF1zAXXfdNajM9957L3feeWf0b7fbTUVFBa2fdiELsaeOBCHSpaLviW0lqLLj2Zp+87qPeNn477ujPc6Enlrm7Qe62PnrAKJFjOTGBSJfCF9rgLb9HqQexRb2R8Ylu4Qj34ooiz2KQsVT76doXg45kzIRLQKSLFKzvoGm7R3M/cYkMooi//aaqrHr+cN0HPAw/5YpCFLEIlBDKnvX1hHsCjHtyxNQQ2pkPKxy+G+NiBaR3OmZ2LItqCEVJajSstuN7JBwFtlQQ5F7lYBKoDPU00CTtLM29G4h6ZFhKaQtW7bQ1NTEggULomOKovDOO+/w2GOPEQgE+nUZKC0tpaqqin379g04b2lpKRaLJea5s2bNoqGhgWAw2M8pfu6553LxxRdz3333cdNNN51UZpvNNuQGeR/+/LPhO7W1SPskejIxFL9Kd0P/45XdDf64482fdND8SUe/8YOv18d9uW1PHECUI0pNlAVCPT3H9r9ej2QRo5ZcyBtG08BT543eK1pENA12v3CYiecXkzPFgiiLaKpGy243ziIb0780ITKHLOJrC7L9yQMUz8th6uWRE89qWKVtbxef/a6WstPzKVmQG7Hiwiqtn7o59kErJQtycU1wRpVg294uOmu6yZ+dhTVTjipNd62XQEeIzHJHVOGpYY1ARxA1rCHZIscTtBEW+Nn9/GHSuLz0oCiaiKL197oYph7ShRdeyI4dO2LGvvGNbzBz5kzuvvvuuC1PWltbqa2tpbR04MJnZ511Fs8//zyqqiKKkX/gvXv3UlpaOmCE7qGHHmLevHlMnz59OEsYGAFs2Ra8LYHEzDdGaIqGomhRy6uXrtr4ybJ177f2GxMEgT2vHO33RXcf9rL5v/b2u7/x4w5adruPK7We57Xt7cLbEohYfbKIvyNiWQbcISxtwahi6z2W4SywkVnqiMwji4R9CoGOEBMvKMZRYIuO73r+EO7DXk7//syIwtQiimrr4/tQAiqnfnMKavi41bfz2UNYs2QmLSmJjoe8Co3b2hEtInnTMqNKMNgVpm1vF7ZsC44C2/H7u8MEOkJINhGx19IMp3aDABUBNY4beLhtkFKJYSkkl8tFdXV1zFhGRgb5+flUV1fj8XhYsWIFV111FaWlpRw6dIj77ruPgoICli1bFn3ODTfcQHl5eTRa9q//+q+sWbOG22+/ne985zvs27ePlStX8t3vfndAWebOncvy5ctH3XalF9ku8blbp7HxZ7t0b+af+s3JbP+fg/iGqHw1VSPs659r4m8Pxs2mb9/noX2fp9947TvNceff+etDccc//PlnMdZg0N1jDf6pLqq8ejuLqCGVrjovQu+YojH7q5Uc29SGs9AencffEaRtbxdZlU4qzimMztF+wMO+V+uoOLeQCYsKozIc+7CVg3+tZ+rlZeROc0WVYN3GFpq2dzDp4hIceccVW/1HrXTV+ZhwTiGSJaLY/P4AfDiUf2mThJ7UliSJHTt28Nxzz9HR0UFpaSnnn38+L730UsyZoiNHjkQtIYCKigrWrVvH9773PU455RTKy8u5/fbbufvuu0/6ej/96U95+eWXEyO7JSKP3pUR9LRBSgNnrxJUUYIAscqw81B3v3tD3Uo/a7B4Xi7eRj8Nm/ofAm3e0Unzjv51xQ+tb+Tw35uiCk/tOa9V+24zjR+3Ryw+WYj6GTsOduNrDSJKAqIsEO6xXEVJQHbIiLKAPEYfKcM7teOxYcOG6H87HI5ov++hPqeXM888kw8++GDA5zzzzDP9xqqqqvD7E5MOG8ljS/0v6WgRRAFBNEqUbWRO7Xjb4kBHiEBH/5It7fu64s7Rt3tNWBubUi8D+5DS971N3xNUCUaQBEIe/dfOEWWBsF9JCwtptPT6jUzSBzO5tgdvUyCuQ1dvKEGVD1Z/mmwxxoXNv9T3+xlxauurhK1pIfVgcUpkVTqTLcaYI1oEciZnJluMcSFnSmZanaAfLuoAaSPxIm/pQvpKnmBcE5xMvkT/Pdls2VZmXVuRbDHGHFEWqP76RCSb+RFPJ8wtWw+GSaeQjbFOQY5YRpqOo6Z6dGqbCqkHwxRnsxinOBug68ipOsD2zDAHI/WMElTxtab2Ke1EoKka3Y36rxwuCETWmb7fzUFRNAElTmZ/vLF0wVRIPbTudtO6W//VIj11Pna/cDjZYow5wa4wH//3/mSLYTJMTI9fD5llDkNE2ex5VvKmuwa/Mc2xOCUK52YnW4wxxSzQpmPyZ2Xp/gMMkFXhpOwM/VZR7MWeb2PihSXJFmNMUTVxwCtdSV/JE4xRok9GqREUWad+Hdp6xfQh9SBZRILdBkgdMUo00QCKd6DtmZLGnnxTIfXgbQkQNEAuW6AzRJfkS7YYY07Yp9B5uH9VAD2hEj+ils4/N6ZC6uHYB/0LmemR1k/1H0kE6Drqo+uo/hWv3jB9SD0Uz8shozR+nXA9kTstk9yp+s9lyyxz6D5I0XswMt6VrqSv5Amm6NRcnIUGUEhTXYboPecqd1AwW98KqTd1JN6VrqSv5Ammt+OH3pEsojGiiQZJkUkE77zzDldccQVlZWUIgsDatWujj4VCIe6++27mzp1LRkYGZWVl3HDDDRw7dixmjkAgwHe+8x0KCgrIyMjgi1/84pB7J/bFVEg9GCX6JMgCigHC4UYI+/fWQ4p3DYfu7m5OPfVUHnvssX6Peb1etm7dyv3338/WrVv54x//yN69e/niF78Yc98dd9zBK6+8wosvvsi7776Lx+Ph8ssv79fsdTBMp3YPrZ+68ccpUao33Ie78bX1L8yvN7ob/HFLzuqJgbP9h2dnXHLJJVxyySVxH8vOzmb9+vUxY2vWrOHzn/88R44cobKyks7OTp566il+/etfR7tZ/+Y3v6GiooI333yTiy++eMiymAqph741kPVMw5b2ZIswLrTtjV/rWk8MfA4pMtbbNr6X4fQoPBmdnZ0IgkBOTg4Q6dcYCoVYunRp9J6ysjKqq6vZuHHjsBSSuWXrYerlZVhd+tfPFecU4ip3JFuMMafo1BzyZ2UlW4ykUlFRQXZ2dvTqbTs2Gvx+P/fccw/XXXcdWVmRf9+GhgasViu5ubkx9xYXF9PQ0DCs+fX/DRwiRfNyOPpeS7LFGHPyZrgizTDr9H1GJ7sqg0BnSNfnrlRNQI13MLJnrLa2Nqo0gFFbR6FQiK9+9auoqsrjjz8+6P2apkWbhA4V00ICEECUjOHUFuX06Mk2WozgvB+spnZWVlbMNRqFFAqFuOaaa6ipqWH9+vUxiq6kpIRgMEh7e6w7oKmpieLi4mG9jqmQON4kUs/VBXsxSqleSRZ1Xb52POlVRvv27ePNN98kPz+2WsSCBQuwWCwxzu/6+np27tzJokWLhvVa5paNSBXFQ39rNISFdOyD1mjXVT3TvKtD9+scqNTIcMuPeDwe9u8/XsyupqaGbdu2kZeXR1lZGV/5ylfYunUr//d//4eiKFG/UF5eHlarlezsbG6++Wbuuusu8vPzycvL4/vf/z5z586NRt2GiqmQiDQUPPpu/L7zeqM+TltpPdKyS7++o14S1Up78+bNnH/++dG/77zzTgBuvPFGVqxYwWuvvQbAvHnzYp731ltvsXjxYgB+8YtfIMsy11xzDT6fjwsvvJBnnnkGSZKGJYupkABbjoUpl5QZorTrnOVV7HutjmCXvisbTL60lJadnbiPeJMtSsqzePFitJN0KjnZY73Y7XbWrFnDmjVrRiWL6UMCZIdEpgFC4RDJZRtu5CMdyapwItuH9+ucbuixYqRpIdGb36V//5HY06tM79EnMEY0USH+9mx4yRqpRfqq0gQiysZQSL3NE40QZTNCxUg9YlpIgLfZz6E3G5MtxpijhTX2vnLUEMq3Zn0DvhZ999lLVJQtlTAVEpEeXm1d+s99UsMaTZ90JFuMccEIPfYSlVybSqSv5AmkYE4W075YnmwxxhxHvpV5/zwl2WKMC/P/dSq2LEuyxRhTtAFKj2jDDPunEqZCAiwZMrJD3xEZAMkmGSKBGCCjyD6kcLVJamGMT+cgGMWpbaTec4Dundp63LKZComesL/OP7zQW9bVCIo38oXU+3s6WLZ/OmIqJKBxW3v0V1XPeI752Lu2LtlijDlKUGXHszWGsHr1hqmQiDRPNAJhn4LHp+86SBBJlu48pO8mkTB4xch0JH0lTyBTryij4tzCZIsx5hSdmsPsr1UlW4wxx1Fg4/N3zUy2GGNO75Yt3pWumAqJiA9JU/TtbwCQrCIYIPIkWgTSOPJtaMwtGxFnryGKs8kCis4dvRApzmYE5/1AXWrTuXOtqZAwTlNBI/WeM0K1SEUTUOJsz+KNpQumQgLD5Hcd+6gVUUrfD+tQ6ar1svPXh5IthskIMBUSkQJURrCQUEFR9K94AcL+dC7CMTT0eA4pfTebCeSUmyaRPSkj2WKMORMvKqbinKJkizHm5M/OZvZX9R9N1AYozqaZJ7XTG6P4ViJFy/RduhZ6ayHp//1MVE3tVCJ9VWkCMU4um1HWaRZnS1dMC4meXmUG+AAbZp1GCftr8f1Fahq/xaZCAj5Y/Smk8Zs4VD596UiyRRgX6t5vMUQjAz1WjExfyROFAJmlDkOc7HXk2yKntXWONVPGkqH/+lZ6RP+fzkGQrCLzvjUl2k5bz8xYNoEcA0QTy88soHxRQbLFGHPiVYvsvdIVw2/ZREtv7Rz9+xxEi4BihAJtFoGwX//vpx5PauvfLBgEURZQFQ1N/59fwzh7RckY69QjhldIgigQ6AgmW4xxQQmqKIH0/qLKpSWD3qMqKmGfEU5qm51rdYe/LciWx/YlW4xx4eP/3p9sEUaN0tY+6D0H/lw/DpIkH5UBUkfS2IeUvqo0QVgyJHImZyZbjHEhf6Yr7Uv1aoHBmz+6Kpy6b4GkVwyvkDJLHUy6ePBtQNojwKxrqwzR7mniBcVkTdR/NFGPfdkMv2UzUjoFYIg2SIJBctn0mO1vKiSjpFMY6HiDZJACbeZJbR2iBFS6G/3JFmPMEYCuYz5DKF9fW5CQV/9VDfSI4S2ktr1dtO3tSrYYY07Iq7D9yQPJFmNc+Ox3tckWYVzQ45bN8BaSq9xhiCibxSlR/LncZIsxLhTPzzWE816PqSOGV0i501wUzMlKthhjjj3XSuVi/VeLBJh0cQnWTMMb/2mJ4d810SIYIvJklKqYYJwCbXrcspkKySKiBvX/RTXKlxTBOLlspkLSId5GvyEyw0Nehfb9nmSLMeYIokDrZ+60z9kzKoZXSA1bBs+N0gOeYz48x3zJFmPM0RTNMJUx9WghGd6pXfy5XFwTHMkWY8xxTXBQPF//UTbJJjJxSXGyxRgXehVSvCtdMbxCKpyTjbPQnmwxxpzMMgd5013JFmPMsThlys7IT7YYacU777zDFVdcQVlZGYIgsHbt2pjHNU1jxYoVlJWV4XA4WLx4Mbt27Yq5JxAI8J3vfIeCggIyMjL44he/yNGjR4cti+EVkigLKAaIPhkpZ88QznsifSniJ9cOj+7ubk499VQee+yxuI8//PDDPPLIIzz22GNs2rSJkpISLrroIrq6jh8ovuOOO3jllVd48cUXeffdd/F4PFx++eUoyvDqUhneh2SUcLhhcvZkYxzjgMT5kC655BIuueSSuI9pmsajjz7Kj370I6688koAnn32WYqLi3n++ee55ZZb6Ozs5KmnnuLXv/41S5YsAeA3v/kNFRUVvPnmm1x88cVDlsXwFlLzzk58rfqvGOk55qP9gAFSZHwKDVvbki3GuDCYD8ntdsdcgSHUkjqRmpoaGhoaWLp0aXTMZrNx3nnnsXHjRgC2bNlCKBSKuaesrIzq6uroPUNlVApp1apVCILAHXfcER276aabEAQh5jrjjDOGPOeLL76IIAh8+ctfjhnvnfehhx6KGV+7du2oenDVbWzB1zL8NyrdaN/noWWXO9lijDmBjhBH3mpKthgpQUVFBdnZ2dFr1apVw56joaEBgOLi2EBBcXFx9LGGhgasViu5ubkD3jNURrxl27RpE0888QSnnHJKv8e+8IUv8PTTT0f/tlqtQ5rz8OHDfP/73+ecc86J+7jdbmf16tXccsst/RY/UqZ9qZwjbzcR6AglZL5UpWRBLiGfQutufSulzFI7+bOyOfz3xmSLMuYMtmWrra0lK+t4WpTNZhvxa534o69p2qCGwFDuOZERWUgej4fly5fz5JNPxlUMNpuNkpKS6JWXlzfonIqisHz5ch588EEmT54c954lS5ZQUlIyIk0/EAWzsxFl/e9csyoycOSN/AOZLthyrWQboPccDL5ly8rKirlGopBKSiLVVE+0dJqamqJWU0lJCcFgkPb29gHvGSoj+ibedtttXHbZZVEH1ols2LCBoqIipk+fzre+9S2amgY3oX/yk59QWFjIzTffPOA9kiSxcuVK1qxZM+SQYiAQ6LeXBjj9BzM54+5ZiBYBW7YFBDj1m5OZe+MkZl9XxcyrK5DsIrJDYvIXSqi6sJiK8wopPzMSUrZkSBRUZ5M3w0XO5EwyyyJnmSSriC3HgiVDQrKKCCmi6yJObQM472XREMXZxotJkyZRUlLC+vXro2PBYJC3336bRYsWAbBgwQIsFkvMPfX19ezcuTN6z1AZ9pbtxRdfZOvWrWzatCnu45dccglXX301VVVV1NTUcP/993PBBRewZcuWATX0e++9x1NPPcW2bdsGff1ly5Yxb948HnjgAZ566qlB71+1ahUPPvhgv/Ft/3MQu93G/Fum4m2JFGg7+l4LokVEsgjRD7ZgEdBUkG0iYqYcNUGtLgtlp+UhWkREi0jIE2bHszXkzXAx/csTEMTIfR01HnY+d4iyM/KZcFYhalhFDam0fOrmyFtNlH4+j5xJmahhDTWk0rrHTdueLgrnZmPLsaKGIvd3HvbiawmQVelEkITo/f62IEpQRXZKaGENNazG7TEX6cmm/y+qaJDytQCaJqDF2bLFGzsZHo+H/fuPd6Spqalh27Zt5OXlUVlZyR133MHKlSuZNm0a06ZNY+XKlTidTq677joAsrOzufnmm7nrrrvIz88nLy+P73//+8ydO3dAo2UghqWQamtruf3221m3bh12e/zDhNdee230v6urq1m4cCFVVVX8+c9/joYN+9LV1cXXv/51nnzySQoKhtb+ePXq1VxwwQXcddddg9577733cuedd0b/drvdVFRUEGgPokpw8K/1hDwKaND6aX//ihpWqFnX3zHX3eDnk6dr+o037+ikeUcngiggWgR6S9M07+jAfcSLKAtRBdY7D1rk+IEoCyg9ib4Wp4wz34bYoxwDnSF8LQHKzsjHWWhH6rn/sz/U0lnTzcLvTEe2R2oAaarGx/+9H39HiM/dOhU1rCFKAtYsmabtHRGr7+IS1FBEgYX9CofebMSWY6FgTnZECYY1Qp4wbXu7sLpkHPm2iDINa4S9CgF3CFEWEERhQCWYDDzHfIboyQYMWPtouPWQNm/ezPnnnx/9u/f7cuONN/LMM8/wwx/+EJ/Px6233kp7ezunn34669atw+U6ftD2F7/4BbIsc8011+Dz+bjwwgt55plnkKTh1aUSNE0b8s/m2rVrWbZsWcyLKIqCIAiIokggEIgrwLRp0/jmN7/J3Xff3e+xbdu2MX/+/JjnqWrk0y2KInv27GHKlCncdNNNdHR0RE+RXnbZZVgsFm666SaWLVvGUJfhdrvJzs5mMV9CFvTTKkeQhKiyE2WBoDuMpmlkTXBGlZpoEWjZ7Ua2SxSdmhO9F+DIW01kFNuZcFbB8Tm6wux7rY6C2VlUXVAcHe883M1nL9dSeX4RledGaixpqkb95jYO/qWeiReVkD/DFVVs9ZvbaNreQeV5RTjyrRHrLqzSuL0DT52PkgW5SFaxZ1yjbV8XIU+YnMmZaKrWY1VqeFsCaIqGJUPqsRA1NDX1rb6wFmIDr9LZ2RnjZB4pvZ/hM1/9DnJG/11HuDvA+19ak7DXG0+GZSFdeOGF7NixI2bsG9/4BjNnzuTuu++Oq4xaW1upra2ltLQ07pwzZ87sN+ePf/xjurq6+OUvf0lFRUXc5z300EPMmzeP6dOnD2cJMdiyLEy/cgI7nulv6aQbmqKhKFq/LHd3rReA6V+eQP2WNtAg7FM49kFrvzm6G/3s+WN/31zLbjctcaJztW83UbexJarseg+YNm5to31/V3Tc1xw5VtHd5EcJqxHFKYtRZWLPs2HLtkTHPcd8hDxhJi4pRrKJkXlkgW1PHkANqpz+/VnH161pbPzZbhz5VuYsn4ggCaBpeJsC7Pz1IVwTHFScUxhVYL62ALXvNOOa4CB3miu69fW1Bmjf5yGj2I491xq1BgOdIfztQayZck83kx4FGUy+MtRjcu2wFJLL5aK6ujpmLCMjg/z8fKqrq/F4PKxYsYKrrrqK0tJSDh06xH333UdBQQHLli2LPueGG26gvLycVatWYbfb+82Zk5MD0G+8L3PnzmX58uWsWbNmOEuIQbKLZBTpP48NILPcgfRJYj+omhppknCiEvS1BuMeNo23JQY4tD7+WZVtT8SvAf7eT3f1WH09ik3R8LcH+ezlI5Seno9sE6nfFDkcGfSEadvnid7bW/xflEWsmTKiJPY8JtC+z0NWpZOiebnRseZPOjnydhMTlxRTODcn6husWVdP3futVN8wkcwyR1RRHXy9nra9Xcy6thLZEbHk/D4fG34/tH/T4ZAoH1IqkdDUEUmS2LFjB8899xwdHR2UlpZy/vnn89JLL8XsN48cOYIojj789NOf/pSXX3555PJaRBSDOED1VBlTUzWUgIYSAIj4i9SQRledj4LuMCEP0dpPgY4QDZv7n9zuPNRN56HufuP1m9qiyqwve9fWsXdtXaQAXB/rbu/aukjAo0fh+Voj1mDD1jZku4Qoi4QFfZ9xSySjVkgbNmyI/rfD4eCNN94Y1nPi8cwzzwxprKqqCr9/5C2MjJLHBpEvkRGUrygLhH1juE6NmM9M0B0iXuJR+77jxfDC2tgoJMNv2fSGt8nPgdfrky3GuLD/T3X42/Sfs9f4cYcZ9jcVUnoS8ip0HNB/WVeAtj36T6wFDFEVsxdtAAspnRVSipwjTg4Fs7OY+ZX4UTw9IYhw2vdmIDv136ts6hVlFJ2ak2wxTEaIoS0k2SEhWvWvk0VZxJZlQVP04dQ+GbYsSyT0bwA0IN7xu3R+lw2tkIzi1BYtkS+oMVJHjJEiA5ET2UICTmqnEvo3D05CpKyr/j+8oiyiaZohLKRIxUj9/8joFUNbSI3b2qOpE3om6Amz/amDyRZjXNj/f8cIuI1x7seMsumMUHc4vTfcQ0RTNDx1xog+dTeO/FxauqFqAoLOziEZess25dJSqi7Qfw8vV4WT0+6YkWwxxoWFt083RJ89vWJoCyniAA0nW4wxx0g1gkSLaAhfGUQibHGjbGm8fGMrJKNE2QzSkw0M1pdNhz4kQ2/ZJFlAMUKUzSKgGORLaqSwvx4xtIW055WjSa9pMx607+3CHSezXY9s/q+9BLvMKFu6YmiFJFlFQ5Q7FUQhrf0Kw0LTDBE5BTPKpjtmX1dFzuTMZIsx5hTPz2XqFWXJFmPMkR0Sn79zpiHOlsFxp3a8K10xtEKSLCKKAZy9eirOdjJ6FZHpQ0pfDL1lM0r0KeLoNcg6Ff2vs5eINRTPh5QEYRKEoRUSojF+TXv7uOmdSB6b/tfZi+nU1hmbHtmTbBHGhcN/03+fewBvc4APVn+abDFMRoGhfUjZVRnRLhJ6xpFvxerS/2+PZBPJKDFGFxnoqYc0wJWuGFYhibLA3JsmIdn0/09QcW4RRafkJFuMMcdV7mDGVfqvANpL75Yt3pWu6P/bOACiJbJ0Qzh7LUbxIRnDea9n9G/HD4DQGyI2gBPUKF9UwUB5bMDA+7M0/icwrkISBbzNxqidE/YrhLz6P5GOGumTZhgG2p6l8ZbNsAop2Bli6+P7ky3GuLD3j0eTLcK40PqZm9bP4rfr1iN6LD9iWB+SnCGRP9M1+I06IHdaJpZM/f/2OPKtZnG2NMewCslZYGPikpJkizEuTL64FGehLdlijDl5M7Mo+3x+ssUYN/QYZdP/z+YAGKlujnEK0RnNqS3E9xelsUIyrIUkWgRDJNaCcb6oRokm6hnDWkhKQMFT5022GOOC55jPEHWfgl0hQp70tQ6Gix6d2oZVSJ01XroPGSNEvOu3h5MtwrhQ/1FbskUYX3R4DsmwW7bMcjt5M/QfZRNEKD+rwBBFy3KnZpJZZkbZ0hnDKqSsqkwKZmUlW4zhIwgIFuuQbxctIpOWlBgiibjo1ByyJ2YkW4xxQ49RNsMqJMmS+h1HRPvxzHWpuCjyH5qGFgoOfY6enD0jOPAN6dTWU6o/BlZIoiV1P7xSVsRyU/1+BFlGsNlQGptGNJcoC5HOKmn+QR0KRokm6hnDOrW9DX4CpGbXWsV9PP1BC4dhFN11tbBG0/aOBEiV+rhrvfhaA8kWY9wwK0bqiOadnciCJdliDBmpuGhEVlLQE2bfa3VjIFHqUftOc7JFGF/MKJt+KJqfXg7QkW7Z7HlWKs8rSrA0qcmEswtwFOg/ReY4wkmuoRMOh/nxj3/MpEmTcDgcTJ48mZ/85Ceo6nGXhqZprFixgrKyMhwOB4sXL2bXrl2JW0oPhlVIuVNdhvjw2nOsFFRnJ1uMcaF4Xq4hSvUmmtWrV/Pf//3fPPbYY3z66ac8/PDD/PznP2fNmjXRex5++GEeeeQRHnvsMTZt2kRJSQkXXXQRXV1dCZXFsAop0qEiNZ3aicQo64RIgTbNSE7tBBXVfv/99/nSl77EZZddxsSJE/nKV77C0qVL2bx5c+RlNI1HH32UH/3oR1x55ZVUV1fz7LPP4vV6ef755xO6JOMqJKMknBqkfC1Ewv5KikZOx4RBFJLb7Y65AoH4Dv+zzz6bv/3tb+zduxeA7du38+6773LppZcCUFNTQ0NDA0uXLo0+x2azcd5557Fx48aELsmw9m3zzk58TfpPHfG1Bmna3p5sMcaFhs1tBLtSM3KaDCoqYhsePPDAA6xYsaLffXfffTednZ3MnDkTSZJQFIWf/exnfO1rXwOgoaEBgOLi4pjnFRcXc/hwYtOSDKuQmra2p1WUbaR0N/jpbjBGqd4jG0bm+E9bBik/UltbS1bW8WwEmy2+z/Sll17iN7/5Dc8//zxz5sxh27Zt3HHHHZSVlXHjjTdG7xOE2NfSNK3f2Ggx7JZt6uVlOPKHnoIxHvQ9mZ0o8me6KF9UkPB5Uw4BZl5TgWQ3zke6N9s/3gWQlZUVcw2kkH7wgx9wzz338NWvfpW5c+dy/fXX873vfY9Vq1YBUFISKWTYayn10tTU1M9qGi3GefdOIGdqJqKcWstX/Ym3ZByFdjKK9N88UZQFCmZlIwwz5G0CXq8XUYz9LkiSFA37T5o0iZKSEtavXx99PBgM8vbbb7No0aKEymLYLZsoiwbJ7xJSNkUmkfT+uBhhrVESdDDyiiuu4Gc/+xmVlZXMmTOHjz/+mEceeYR/+qd/AiJbtTvuuIOVK1cybdo0pk2bxsqVK3E6nVx33XWjXkZfDKyQUjv6JLpcqAk44yGlcM5eIuktr5LK72nCSVAJ2zVr1nD//fdz66230tTURFlZGbfccgv/9m//Fr3nhz/8IT6fj1tvvZX29nZOP/101q1bh8uV2BI+gqalc3254eN2u8nOzuaa+TfSsdObsh9g0elE9Y6+omVWhROI5HnpGckqUjQvJyWLtIW1EBt4lc7Ozhgn80jp/QxP+K+fIDr6b8dVn5+j3/23hL3eeGJYC6np446UjrIlQhmB/hVRL0pQTUllNJYIWuSKN56upJZXdxw55ebJyRZhXJi4pJjCU3KSLcaY4yy0Mfu6qmSLMb4k6KR2KmFYhWTPTa2Q/1jhyLchGyAULjskQ/Sei6HXhxTvSlP0/0kdACOkjUDqO+8ThVHWqXcM60NSU7x8baIwTM6ekcvXxhtPUwyrkPb/2RhFyw7/rRF/x9BrcKcrnmM+at5oGPxGPWEqJP3gPuxN6ShbojBKlC3oCRP0mIm16Y5hfUgzr6kY/CYdcOq3puDqOYukZ4pOyWH6sgnJFmN8MaNs+iGV+5RJubkJm8uaKYOaxp/QISLZRSSLwT7OZpRNP6hKajlApdxcBFlGKsgHKXFvi1EKtKVyWyuToWNchZQCUTZ5YmX0v5X2drRwGKWlFaWlNWGvYawk4uS/p+NJ70nteFe6Ylin9uG/J7+YV/jQkTF/jU+eOkigU/+VMRu3tCPI6btVGRFmlE0/hLrDSAZYvq81gKak8Sd0iIS8YTT9G4K6x7BbtsrF+u9VJtslFv1oDqIBnL2TLymj6oLEVi80GX/0byIMgCH8KpbeGkEGWKtBCtH1RWCAbP9xlyRxGFYhGcEBKspiJJqo/6Ua0qmdqAJtqYT+bfkBMEZ+l5AS0cTxwJC5bDrEsBZS49Z23Tu1vS0BNv9yb7LFGBf2rj2KZoADoDHoMMpmWAtJlNLXrB0qoiwg2YzxFktWY6wzBjN1RD9kT85ItghjTlZlBnOWG6OK4qxrKsmdmtiC8ybjj773LCfB9CHpCyM6tfVYU9u4CskAH17RIqIYxNFryFw204ekH5Sg/j+8giig+PW/zl6M8COjd0alkFatWhXtatnLTTfdhCAIMdcZZ5xx0nmefPJJzjnnHHJzc8nNzWXJkiV89NFHMff0zvvQQw/FjK9duxZBGL6Durs+8W2rU43mTzrY/cLhZIsxLmx5bB9dBilGF8V0ah9n06ZNPPHEE5xyyin9HvvCF75AfX199Hr99ddPOteGDRv42te+xltvvcX7779PZWUlS5cupa4utsys3W5n9erVtLe3j1TsKKIBEjFtORbDdOLIqnQa4j3tix6z/UekkDweD8uXL+fJJ58kN04xMZvNRklJSfTKy8s76Xy//e1vufXWW5k3bx4zZ87kySefRFVV/va3v8Xct2TJEkpKSli1atVIxI5BckijniPVKTolh/JFBckWY1yovmESVpf+SxLrnREppNtuu43LLruMJUuWxH18w4YNFBUVMX36dL71rW/R1DS8Uh9er5dQKNRPkUmSxMqVK1mzZg1Hjx4d0lyBQAC32x1zAdizLdjzrFhdMoJOPWmR08tp/HM5VITIuTLjObX1VzFy2FG2F198ka1bt7Jp06a4j19yySVcffXVVFVVUVNTw/33388FF1zAli1bsNmGtn245557KC8vj6vwli1bxrx583jggQd46qmnBp1r1apVPPjgg/3GZ15diSPDjiiLbHtiP92Nfhb9eA5qSEMNqyghlW1PHECUBGZeXYkaUiPjQZU9fziKPddK6Wl5qGEVNaQR9IRp/LgdW44FV7kDNaShhFRC3WG8TQFku4RkE3vmiTw21nt9oySc9m7VDKF8+6LDKNuwFFJtbS23334769atw263x73n2muvjf53dXU1CxcupKqqij//+c9ceeWVg77Gww8/zAsvvMCGDRsGfI3Vq1dzwQUXcNdddw0637333sudd94Z/dvtdlNRUcGmX+yJdB0RiL6Bm3+5F1EWEC1iJGTuV9BkkWMftiJaBERZjD3hLYDslJFkAUtGZAvoyLdRdnoBokVAsoh46n3s+cNRihfkMmlJSfSpjdva2fdqHZOWllA4NzuqCOs3tVG/qY3KxUVklNgj4yGVxu3tuA97KT0tD8kuRRVb+94uAu4QOVMyQaNHQap4mwOE/QqCEKk3rYY03dZFEkQBX1vAEGfL+mL4c0hbtmyhqamJBQsWRMcUReGdd97hscceIxAIIEmxvpnS0lKqqqrYt2/foPP/x3/8BytXruTNN9+M6yzv5dxzz+Xiiy/mvvvu46abbjrpnDab7eSWWZ83L15lRSWo0rKrs9+4vz0Ytw9YxwEPHQc8/cbr3muh7r0WBEmIlAXped2691to3tkZVXiBnh5q7lovoe5wRDnKAmrPMQWLU8aea4koR4uAp85LwB2i4txCLE45qlB3PFPDsQ9aOfPe2VQujtQJ0jSND1Z/ijVTpvqGSVGlFvKG2fncITJK7FSeWxRRamENf0eQ2neaySixkzvNhRZWUUIa/vYgHQc8OPKt2HKsMVZi0B1CsosIghCdZ6x/sZWAypY1g3++TFKfYSmkCy+8kB07dsSMfeMb32DmzJncfffd/ZQRQGtrK7W1tZSWlp507p///Of8+7//O2+88QYLFy4cVJaHHnqIefPmMX369OEsIeloiobSx1IJdoUJdvXvJxZRbP2ff+Tt+P64HU/X9BvLnpTBlsf3EewMRRWVElAJKCE++92RqFLrPTYR9il0HPJElaASiChB2SbhzLdG7pcFPI1+Og54yJmcSelpedHx5p2d1KxrYPLFpRTPOx7sOPT3Ro7+o5mZX6kgs9yBGo5Yg4f/3kT7vi6mXl6G7JR6xjXq3mvG1xqk4rxCUIkqtqZPOlDDGnnTXNExNazibQmQN92F+4g3MtajaPVqEUYx+pbN5XJRXV0dM5aRkUF+fj7V1dV4PB5WrFjBVVddRWlpKYcOHeK+++6joKCAZcuWRZ9zww03UF5eHo2WPfzww9x///08//zzTJw4kYaGiOWRmZlJZmZmXFnmzp3L8uXLWbNmzbAWbCTKTs+n44CH+k1tKEEABYj4WrqO+vrdH+gMUf9RW7/xzsPddB7u7jfeu708kX2v1nHgz8eiCq/3EOrhtxqxZPRYcbKItzFyFqzzUHfMuNqjSKwZMpJVROhReC2fuhEEjQlnF0TnFmWR3c8fYsplZUiyGG1vFfSE+eg/PyN3mosZV06IWnGeBh+fvVxL/qysHh+ghhbWcNd2U/d+K3kzXGRXZUSVXddRHx0HPWRPzMCWZYkqwu4mP4GOEI4CW9ShroY1Qt1h1LCGIApjX31goBC/URTSYEiSxI4dO3juuefo6OigtLSU888/n5deegmX63ji45EjRxDF46Gtxx9/nGAwyFe+8pWY+R544AFWrFgx4Ov99Kc/5eWXX07kEnRFMqNskS+0An3On/pag/ha+7f1bt7Zf0sMcOD1+rjj2//nYMzfmaV21KDGB6t2IYiRdQs9vj734W4+efpg1IrrVXbeJj9N2zuiSrC33bimahHfoENElGQCHZFtvGuCk5zJGdF5jn3UStO2DiYtKSarMgNBFhBlgT1/qKVll5sF352GPTuynfV1+9nwi+H9+xmVUSukDRs2RP/b4XDwxhtvDOs5AIcOHRr0Oc8880y/saqqKvx+/Z+4HimiRTCEo7dvHpumxqYFKUEVb1Og33MGUo7t+zy07+vvAzz6bjNH323uN777xfidY7b9vwNIFhFBFlDEMWrxbfQtm0l60d3gJ+DWfwskJajiPpJaaSNhn0LYF9kih7Uxeg9MhWSSThz8S/wtj97obvCz95WhHZQ1SW10ekbZBCJObVuW/tMpnIU2ik7JSbYY446Zy2aSVkw4qwCLS/9GcGapg+L5/XMqTdIPUyHpGNEiGsKpLRgkRcYI6P/n08AYpayrUdbZDx06tU0LKYURZBkpKwvBYo0ZF0+ZOaTnN+3ojEZ69Iy/LUjHwf6her1j+pBMxhUtHEZxu9FCx8/LhJYu5NgFecgVEwZ9/v7X6gyhkNr3e+KeGDcECaoWWVdXx9e//nXy8/NxOp3MmzePLVu2HH8ZTWPFihWUlZXhcDhYvHgxu3btSsACYjEVUoojulwIfZKDjy624KlS6TijHDFj4FZOklVk6hfLdVvrqS/5M7MomJOVbDHSlvb2ds466ywsFgt/+ctf2L17N//5n/9JTk5O9J6HH36YRx55hMcee4xNmzZRUlLCRRddRFdXV0JlMX1IKYhot6P2nEAXC/LQPF4EmxW1vQNpmoeQx4q3yEZmd//8sl4km0jJ/Fz2v1Y34D16IavKCUDLLneSJRlnEuRDWr16NRUVFTz99NPRsYkTJx6fTtN49NFH+dGPfhQtIfTss89SXFzM888/zy233DJ82QfAAL+f6YeQnYWwsBopN5fwoSNo3d20n1XBwXtPwWkPIHTLOFpUBHng3xMjtQUyUv+5vgzmQzqxUmog0D+FBuC1115j4cKFXH311RQVFTF//nyefPLJ6OM1NTU0NDSwdOnS6JjNZuO8885j48aNCV2TqZBSEM3TjWqRoKwIYcEc1OopdEwVCZaF6OjMwN4oYelWEazWAecwUuRJlIyjfIdDRUUF2dnZ0WugWvQHDx7kV7/6FdOmTeONN97gX/7lX/jud7/Lc889BxCtvlFcXBzzvOLi4uhjicLcsqUgqs+Ppb4d/AH8n6vEnyuhOCLKxfapg9L3A9jq3SjegfO3wl6F2n/0TwbVI237uwh1j1ECayozyJattraWrKzjvrWBChWqqsrChQtZuXIlAPPnz2fXrl386le/4oYbbojed2K7MU3TRtSC7GSYCimVEATQNKS8HMKHjiAV5OP88ACObi/il+fhabdStqETbcsuBoudBT1h6ja2jIvYyaZ1t8F8Rz0MVsI2KysrRiENRGlpKbNnz44ZmzVrFn/4wx8AKCmJlF5uaGiIKbTY1NTUz2oaLeaWLYWQJ5QDoHZ2IbpcaIEgSmsbqt9P7ru1THh2D9qWoYVaM8sdTLmsbCzFTRkqzy8iZ3L8Qn4mg3PWWWexZ8+emLG9e/dSVVUFwKRJkygpKWH9+vXRx4PBIG+//TaLFi1KqCymhZQiCDYb4aN1yCXFhBsaQVMRp09GChWi7K8hfHR40TKby0JmmWOMpE0tsisz4tY20j0JirJ973vfY9GiRaxcuZJrrrmGjz76iCeeeIInnngCINqdeuXKlUybNo1p06axcuVKnE4n11133aiX0RdTIaUKigKaRrihEbm8jMC0EhSbiONgG2jDd05HIk/GcPQapd1TPxKkkE477TReeeUV7r33Xn7yk58wadIkHn30UZYvXx6954c//CE+n49bb72V9vZ2Tj/9dNatWxdTCTYRmAopRdDCYQSLNXIq22bFn29BkwSs+ZmIjpkItY0ow2ghHgn7GyTKJgtoBlnrWHH55Zdz+eWXD/i4IAisWLHipCWlE4GpkFIAadY0aOlAkETCDY2ESrKxt4awdPQUwZ+VQ0a2A8teC0rj0LoAu490G6JaJEDtP5rpbjBeKWM99mUzndpJRsrPQ/AH0coK0FyRVBDh/U+wHW4lmO+geWE2tk4FS5sXXAOnipyIrzUYtz+cHmnZ7TaM8o0hXh7bKPLZUgHTQkoySmsbkqqh1dahqhrStMkoBw7jnVFE2CmS/4kHNu9GUYeXJFt6Wh6WTJkjbw3Nokpn5t44if1/PoavJf5JZN2iw/IjpkJKNoIQ9Q2J82ZDSEGuLEdo9ePc0Uq47tiIprVmWZDtxjCAM0rtJPh8nkmSMBVSshFEpJwsEAUIKaj7D6EFg3DoCKM5e2yk/C4jpcn0xfQhmSQeVUHpdKP5Ayi7eg6njSDMfyKGSq41ai6b6UMyGQtEqwW1uxu5vGzIUbTBOPZBq/572/ew67eHCHn1X4jOCJgWUhKR8vMAorWPwnXHEJ3OhMztawngb0+D08uiNOop2vd7DKN8+2KWsDVJKErr8bKrYs+JV8WdmETRmVdXpEVrIGlK1aieb82UOeOHsxIkTZqhwy2bqZBSBTWxPhDZLqElwBc1lgg2G+qh0XWcFS0iotUMsekFUyElGe+VpwMgFhUg5SbOokmHyJMWCCAVF45qjnRY55hhWkgmicb1zn6E+XMIleYkdN50aRI53CoGJ2JkhSSc5EpXzChbklFaWqGlFQEGLbo2HD596YghWiB5mwPseKYm2WKYJAhTIemUcEBBSQMLabRoqkagIw2iiWOBDlNHzC2bTllw23QyS1OrQJtUODp/UTxyprg49ebJCZ83HdBj2N+0kHRKKhVok2ZNA0DwJr5EiJF9SHq0kEyFpFNES4p8UQUB78ScSLG5jiAWT3fM+avRYqRCdEbA3LLpEaEnvysVLCRNI+PjIxFl1OZNqDICA5ev7UVHIX8wLSR9osHGlbtSJts/3NCI2NCY0ChiL03bOmje0TkGM6c+Zra/SVogiODIi98UMFUQ7faEzCPbRWSb+THWC+Y7qUOsmRbm/8vUZItxUnoTikdLycI8qi5MbLPCtEGHJ7XNLZsOEWTBEGeQAETZuE5tc8tmkhZIFtEwbYEM79TWGaaFpEdECHqM04VDCRhUIZnnkEzSge56P1sf359sMcaFmnUNyRYhaZhbNpO0wJIhk1WZmMqTqU5GsR1btiXZYpgkCFMh6ZCsSieTlpYkW4xxofK8IvJmJLa/fNpgRtlM0gEj5XcJBlprP0wfkkk6kC7F2RKBkSKKJ2L6kEzSAiWg4Gs1RlvpgDtEyDualpomqYRpIemQll1uWnYlpntJqrP3ldE1CUhrdLhlMy0kHZJZ7iCryhhRtoI5WVgyjfm7KmjagFe6YiokHZI/M4vCOdnJFmNcmLikBEeeNdlimCQIY/606BwpVYqzjQNGzmUzt2wmaYEoiwZKrh1eLpuUm4uUlTWGEo0fZk1tk7TA2xwg6DFG5KnzcPfQ2z0JAkp7+9gKZDIqTIWkQ4592JpsEcaNT188MvSb09jZGxdzy2aSDhTNyyGzLLVaII0FgihQubgIUU7nXq0jR49bNlMh6ZCiuTk4C1O7hG0ikKwilecVgWBMhaRHzC2bDjFKa6Bey8iwBdp0uGUzFZIOMUoVRUEWUBUtrb+Ao8HMZTNJC1o+deNv03+/ezWk0bA1sX3e0ooxKD+yatUqBEHgjjvuOP4ymsaKFSsoKyvD4XCwePFidu3aNRrJB8RUSDrk6D+a8TbrP7k21B3m4Ov1yRZDN2zatIknnniCU045JWb84Ycf5pFHHuGxxx5j06ZNlJSUcNFFF9HV1ZVwGUyFpEOmXFZmiCqKjnwrEy8yaAukHk4WYXO73TFXIDDwj5TH42H58uU8+eST5ObmRsc1TePRRx/lRz/6EVdeeSXV1dU8++yzeL1enn/++YSvx1RIOqRwbjaSVf9vrTXLQv4MfZy6HhGaNvAFVFRUkJ2dHb1WrVo14FS33XYbl112GUuWLIkZr6mpoaGhgaVLl0bHbDYb5513Hhs3bkz4kkyntg4xSuqIofPYhkBtbS1ZfdJkbLb4R0FefPFFtm7dyqZNm/o91tAQaaJQXBxriRYXF3P48OEEShvBVEh6QwBREgxRRdEo0cSBGCzKlpWVFaOQ4lFbW8vtt9/OunXrsJ+kvblwwlkvTdP6jSUC/dv1BkMQ4NCbDYT9Q8zvSmN8LQHqN5tRttFE2bZs2UJTUxMLFixAlmVkWebtt9/mv/7rv5BlOWoZ9VpKvTQ1NfWzmhKBqZB0hqbC0fdaDLGV8TYHaNrWkWwx0poLL7yQHTt2sG3btui1cOFCli9fzrZt25g8eTIlJSWsX78++pxgMMjbb7/NokWLEi6PuWXTGZYMmZlfqWDHszXJFmXMyZ+ZRUaxnSNvNyVblKQgqJEr3vhQcblcVFdXx4xlZGSQn58fHb/jjjtYuXIl06ZNY9q0aaxcuRKn08l11103GvHjYioknSHbRTLL9Z9YC2DPs+Io0H/O3oCMU+rID3/4Q3w+H7feeivt7e2cfvrprFu3Dpcr8f3wTIWkM0TZOC2QjO7UHis2bNgQ87cgCKxYsYIVK1aM+WubCklniKMsXyvlZKP5Awh2G0pH56hkkSsmoPkDaOWFiM0dKI1NaOHEFY4Te3PZDIoec9lMhaQz/G1BDrx+bOTP//w02qdZKXm7DQZRSPKEcpSGxrhKRsrKovbqSoLZkL1fw9GWQcYWhXBDYyQUmIBiaS2fuhElA5ce6XMIst94mjKqKFu8RLybbroJQRBirjPOOGPQuf7whz8we/ZsbDYbs2fP5pVXXol5vHfehx56KGZ87dq1Y3IeIl0JeRXa9gw9x0iQZdRz5iNNm0xo6UKaPmel49QQvioXYhwfgZSfF3mezQaCEFVGgsWKeOosBJsNwWJFcbtxNGsEisJ4iwWc+9sjyggS9oXprvfTddR3woKM81kwC7T1YaBEPIAvfOEL1NfXR6/XX3/9pHO9//77XHvttVx//fVs376d66+/nmuuuYYPP/ww5j673c7q1atpN+siD0ju1ExmXVM54OPCiad1582kdY6d+qUlNH3OinB6Bzef8Q8aTpcRnLHOcbm0BK3bS+DS01DOmA2WiIEt2u0IVgvdk1ygaohZmUiFheR90gGAagVlz36AhBbYr7qwmNLT8mIH09g6MBmhQhooEa8Xm81GSUlJ9MrLy4szy3EeffRRLrroIu69915mzpzJvffey4UXXsijjz4ac9+SJUsoKSk5aU6O0ZGdMrJDOv53eRkQsYQAlM/PRq6qQJBlpBlTaat24S8U8OdDIEcjP8PL5a7tnLl0J97PVSFYrMgTygHQMp2ofj9hp4gQUgkfPBSZe1IFDTeeij9HQpxaheB0opUXou0+gLVNwtGoIWVlIU2bjDJ7InJVRYxMotOJ6Bx+Y0tbtsUQOXsDMgblR5LNiN7NgRLxetmwYQNFRUVMnz6db33rWzQ1nfycyPvvvx+TvAdw8cUX90vekySJlStXsmbNGo4eHVoL5UAg0C/rGaDo1BwKq7PJm+FCskf+GTJL7TgLbdhyLFgy0tO91jfypJ4zn+aLqhBdLqTSErRFp9IxzY7mtBM+5xQCE7LxFQmoMigODU3WOHIsn394p/NvpX/B/5121IWzCNdFfFKCP0j4wgXIXhV5T23k9ebNxjsxB9UKjhYFIRACVUWzSAhzppK9F5ytCorbjbLvIHzwCeHaY6jnzUd0OpEKCxGcDgRZRrAMr+FjZK1p/O0bJXrcsg37W3eyRDyASy65hKuvvpqqqipqamq4//77ueCCC9iyZcuAyX0NDQ1xk/dOPK4OsGzZMubNm8cDDzzAU089Nai8q1at4sEHH+w3njczC4fTjmQR2fPKUfzBILOvm4hoEZAsIoIo8OF/fIpklVjw7WkoIRU1rKEEFLas2UdGiZ0pl5SihDTUsEqgM8TBv9STWeagsDo7cn9Iw98eoGWXG2ehDWehLTJHSCXoDuFrDWLJkHqSRCP3KyF1VL9woiygCRKi00nLbDv+PIHub88l+6CKZ4KIoEA4x4m3yErnFBHFrhHK0lAdCigClX+Q+K+GS/n0wlIun7CT1348l86tZzBhQxD+vhXllFIcRz0QCiNNn4KvJANR0Zjwu0MojU2Eex3cdcdgwRyyD/iQtx9ABaTiIsh20XROEaFMgfKdMoLdRrh2aD8u/dcqmmF/nTEshTSURLxrr702+t/V1dUsXLiQqqoq/vznP3PllVcOOPdwkvdWr17NBRdcwF133TWozPfeey933nln9G+3201FRQWfvXQEWYitGfTRf352XB4xkoYR9ilsXrMXySIiygKCHLGmgu4Qxza1IclCbA1rTUMQIyemJVlE6Nk9ZZY6KDktL3p/294uatY1UHFuEaWn5UXXeuSdJo681cSMr1SQXZXRo6hUDm9oonW3mymXl2Fzyag9ivDoxha8TQEqzi2MyC0JeOo8yHYRRI1SR4iAU0PLFMlExdcg4Z7mgDywoRJSBDQfqEEZ1arRfIqIGNZYv28m+TkeWtszkWRonWOj+G8a9jc/QcjKgpLI6zm2HELzetEyM/pF27QtuxCB0OLP4Suy4ssX8BcI2Js1CnYGEJxO1J5InlRYiNLcPOj72ZfDf2sk5DVG/7m46DDKNiyF1DcRrxdFUXjnnXd47LHHCAQCSJIU85zS0lKqqqrYt2/fgPOWlJQMK3nv3HPP5eKLL+a+++7jpptuOqnMNpttQMvsZGjq8f8PdIT6PR7yKrTs7B8W99T78dT3t+yaPumgqcfJ25eDf6nn4F/qESQhYt30vO6h9Q1YnBKiRUSURXwtkeJaHfs9WDJlJIvQcwgy8uGTHRKyTULsGZdcGqIAc7Qgkk9DUkFWYdN+K1k79jH/5uOO75AX/nrERWaZypmWbpQAqLsE3Eom22ptlKsBqjJ9KFdXoIY1PPU+GrbVkTMrn6zpmYQ7NMjNxXNAxX3Yh2tGIRYpjBoMo3T78O/cg7JgFnZRwHpYQegOIx1tR3Ha0Xr+rbTu7mG9PwDdjf5hP0dPGP4cUm8iXl++8Y1vMHPmTO6+++5+ygigtbWV2tpaSktLB5z3zDPPZP369Xzve9+Ljq1bt+6kyXsPPfQQ8+bNY/r06cNZQsqiKRpKn0N+gc4Qgc7+irD1M3fc59e8EfliV11YjGQRCLT5kTrhjZIMglkajgaRnP0KWe/sJJCTzYZ/WAkWWsh6dTOCJJA/aSrSToHdRVZECeTdB1AVlezDXpR8K20z8hBRYVIJisOD6m9FtAhYcu3Y8u0IXg9KhkpHdzc5VeXkTspERENCpe69JgJvbGbO1yeSe1pmj8TZ7F17lKZwmHn/PAVrlgU1rKIpsOePtXjqfMxZXgUQtQYPrmtACShUXVCMGtIoPCWbjoMeDvzfMSSrRO50V2SOsIYSUOk83I1kE7FlWVDDkTmUkIriN7d5qcqwFNJgiXgej4cVK1Zw1VVXUVpayqFDh7jvvvsoKChg2bJl0efccMMNlJeXR6Nlt99+O+eeey6rV6/mS1/6Eq+++ipvvvkm77777oCyzJ07l+XLl7NmzZrhLEH3SBYRVdFQu7rIe2Mf2fvKCObasNd7UD/5DBVA07C814ocCBDd8OyMPNarAvtuoXytQUK1Mq3nTiDzWBDrx61IWVm0bm+m5aOjyJMnEj5YG5Xh8Ct7OTpnBjS2oHZ2IWZFlNCu3xxCyslGUMMQ9CNk5yBNn8Jnv69FtkUsQdEi4G+NNCho+qSjx0IUov6i3q2t7BCRbRL2XBuaBrJTomR+bs+9AkpQ5ZOna8iqcDLjqoroHF3HfGx/8gCln89j8sWlkS1xWKNtTxf7Xquj7PR8Ck/JiSq2lk/dNGxuo2RhHpml9h7FptG21437sJeCOVlYMuTouPtwN4HOEFkVThBADWsEgiIMbzc6NMw2SCdHkiR27NjBc889R0dHB6WlpZx//vm89NJLMYl4R44cQRSPB/gWLVrEiy++yI9//GPuv/9+pkyZwksvvcTpp59+0tf76U9/yssvv5zIJaQ9oiwQDkRqISktrdDSihXoaxOoXu+g8/TdQkkF+aCqeCpEnM0iqs8fOTaQl0u4viEa/u+L0NqBJohooSCq2wOiBKoSm47ib4XGVgaq3NS8I/5J8V5rMHeKiyNvNYIWOaEer8JB+34PH6z+9LhcPSe7m7Z10La3K6IEexQYQPsBD4HOUFSxeXuUY8gTJtgVjo4LYmQeZ5GdjCJ7VOEFuyKWbeX5RdhzrYiSQEgJ8ZdHB1jkKNDjlk3QtDT2gI0At9tNdnY2i/lSjFNbtNtR/envk5j+5Ql4WwMc/UdifpKlgnzci6fSWi3hLw5jbZOY/LPtgyo1wWZDC4VBHbtCcZ+/awa7XziC55hv8JuTSFgLsYFX6ezsHLSC41Do/QwvuvgnyJb+waVwyM/GN/4tYa83nqTnYZsxQPX7RxTpSTUOvlGPpibuN0ZpaSV7Sya+vDIsbhlno4oWHLznm3aSDheJ4pP/rSHg7u9nMwyqFrnijacppkLqQ7oro160RGfAqyqFT29BKshDCyso4TCCxYoWSm4zyrBfSfxa0wkd+pAMfO5en8y8uoKCOdmJnTQQRCovQW3viCrtZCsjgDN+OAtrlv77zw2EwAAntZMt2CgwLSSdMRanl7VwGPVoa0JrGY0WUY587cyT2vrCtJB0higL0cOSo57LbkeaNQ3C4RH5pXqTZ/siTyiPRO1GidCjkIzQ7mlABmkUmY6YFpLOiKSxjNxqkCdVESrNQXIHUHd+hrD/MFooiFxSjNrRiRoMDTlypoXDyKUlIMuobe0ET5+JP0PC3uhDLClE3fnZ4JMMgNiTwmNkC0mPYX9TIemMrf/fwCk68RCdzkihtUAAMTMDTZbQLCLq7sg8WiiIaLcTbmhEyslGcrmG5fwP1zcgT56INnMiYkgl46NaCAQQXC6Gq0oEWY5uG0OeMBt/tiuaamOiD0yFpDMyS+14mwMxZTmkrCwUd/yUEwDV50eumgDBEOF9B5EOWUGS0AApMyP63JHW2FZqjyEcExEDAcjLjcwzgrn6+rBEWcCeZ8XbNPbHC1IWM8pmkupU3zgJe25sXaFehdKvWiREzhSpCuGaw9G6R2gqWiiIlJmB6vUOu05RFEFAtNtBU1H9fgTZgtLahlRYOLL5iPi1ABwFNubeOGnE8+gBQdMGvNIVUyHpDMkiRmoqxUPp7/uJ1sTu44DuHdM0DTE/b8QhfrmqAjE/D6lyAgBiTjZyeRmC0x6tzT1cek/Ty7muWCuwuGhE85mkFqZC0hGCCIIYW0Wxr3VzsrD9iY8JsgxV5ZA38jNNapYTtSAbenLfNK8Xrbsbrb0TpbVt4HVYrBFneByk/LyIMvN5YxzavvlVhirwD0QSFAe60hRTIekIQRIIuEMxjSJHat0IVivsP4RqtyJNmzy859psyBUTEOqaULd/Sri+AQQRITMDpaPzpP4sAKmiDM0f6xuSciKKUWltiygzAYKeHktu0akEsvuXvoHjNcWj8yTgyEGqoMctm+nU1hFqSGPTL/YkZi6vN2JxfLxrwGz8gdACgWhZWsFiRQuH0EJBtM6TK6Loa2faUU+oIKBOmYDwiS8yz1nz2PN5J5/kaoT/sxitKADN4Do8Fz74JFaW7Eyoi/y3eOos6PZDS+swV2QyXpgWko4QLSL5MxOY3a1pSMVF/ayMYU3Rx0JTh5BwK+XmInhj7xOdTrorM5EK8kAQkJvc2C0qhWEFTQDVK6NmhTl2biZyaQlScVHUElL3HYpsASdPpOn0HNTczHgvm56YXUdMUhlblsz0KyeMfqI+vhilsQmlsWl0/pneLcQAWwnRbo90IMnJRp06ATo9USe7YLGier1YO8NomU7QNJR9B5mxfg+fr3cz7e6tTPyjRvZ2K7ZWjdDEYgSLBbWjE3nyRASLjFhVTuuiEoJZAsHsEUYMUxHzpLZJKiNaxBj/0Yg54QM92hw2qbAQwWYFTUMLhyNVHy0W1Pwswi4b/mwLCCD5FBz7m1G9XqSKctSmlp7XD2F5ZzuKokTPVImygNLlRQsFsdd7KG63ECi0I3y4k3DvSXJFicxltRDMFHC0aDgOtESqZPYoWKmgIG2rPJgntU1SmlTtU6Y0N0ec0qVFqHYrojeAdvgoWkE2ngk2QpkR5SCGJOyH5Z6zTzJqn6qV0aMIPbWY+q5V/SSSgnLiKSu1pQ15UhWBokwKtnuRd9WgnmCtpasy0iumQhovekq4jiVqWKPr6ODlaZOB0tEJbk+kjG3PmNzeReZRKx3T7GQ0hBEUDf+kPCz7DkJr/HbpveeQQj4FX+vAPinR6UTt7kY95EWqORyRIaErSgGM3gbJJILodCJkZER+XQVh0A+A6HIh2KyRGtdjSHeDn89+Vzv4jeNNrzI+QSGHDx1BPHSEfP9cgtlWFIdExmctKIBgjxye7D2vJE2dRHBCLrUX2XA0CLS7NQr/PvBa1Z5jA3JJMeGqIoSAgtThIXzoSFp/YfsiqJEr3ni6YiqkEaB6vdBbU3oIH26htAi1ZuwVhSPfSkaJnZZdQwuvjxuDWIbaph1YRAlbhhOlqysylpVBcFoJ1k9lmq6YgmIXsHZqhFwq0oxuCgJhWhYUYX2rHOcrHw74muH6BmRZRuvqIjzCXDyT8cOMso0h0fSIphYQx/4UcWa5k9KFaXrwT1VQu7qQpk5CtNvpOCWfQ5fZ+OznEyj5+iG4qA1PhUD+NpHM/3MxaaOdKU3QtEBEGqSQvebpRnF7xmcd44kZZTMZDkprW48vw4eYlxMJn48hEUdv+trrcnkZWmsH7VfNg+UtfKn4IA4pxIetE/HuzKX04xAZHx9BaWnFekkRqqeAoi4LgdOmIf9ty4DzKp1uUJXIUQJJGpcGBOOCDrP9TYU0xqheL6LTidI09GiO6HKBoiDm5RI+Wjfk50kWMWHVIhOBPLESpa7hpOkrgs2GmJON0thE+Fg9Dd89E89CH3dP2sgc21F+uPcrhF4oZuKz7wNEG1uKsggHjmFnAuK7204uR1EB3nmVaCJkbD2CpqhmdC1FMRXSWNLjzBUL8xH8gSFbSKrHg1xchNLcMqyX87cHU6rCu+a0w7wZqHYZyR9Gau2KNpUU7XbE4kJQVNSWVkSnk44vnYJ7VhhJ0tjknsRrvlPpeLuEiuc3oxGxoFR3F2pXF92NfnytAeSj+wbNJQ03NCL7yvAVWMFhB0mESSUI2/emtbU0UN6amctmEp8ex6ra3Bq3FtGAaBrhhsZhv1zb3q5hP2csUXbvBUGg+5rTERUrUrED9dRi/LkiqgXyPvNjPdIW6YmXm0vT58FZ2I1lQzYbjpxCyYcqE7cfQ5VEtBDH6zUBdRuHp6ylD3aTXViA+/MTCDkEpBBkftSjjIYQKU1JzLC/yUgQrJaIH2OMyZvhQhCg9bMUUkyahuulD0AQkPLz8J4+GcUmogkCcmcAzWlHsFgJzosUW/MddVH55Fa0WVPQPt7FQGfESxbk0lXno7thaN2GtUCA8NE6MguyaTo9m/xXD6Ll56F1ewmeNQfrR3tRu1Lo382gmAppPCjIQ87LjW5XxoqcSZmoippaCqkXTUNpacX5jxDKklk4j/lQt+0GIqVFArkWrO0imbVapJHAx7uiT+17HqmXolNyCPvVISukXoTD9WRMzETLzyFYmEFXpQ1B0ZBnTUQ6UHfSOk0ph0b82kfpayCZCmk80I41onjH/gS1aBEIB1L009izLVK6ushcvzvWGinIw+IOU/ZuGOuHn6GecG4pnpIYaXcVpb0d16ajqAXZeMqt5H/YDE2tKO3taXeSW48+JPMc0jgiZWUhl5dF60LHvSc3N24/s6EgyqkVZetL9EyWpqH5fEjTp+C//PORvm+dHixvbkH++5aY/LWTMZq8vXDdMYSaOvL+/CmCPxA56JqOaAxwDml406xatYrTTjsNl8tFUVERX/7yl9mzJ7aulqZprFixgrKyMhwOB4sXL2bXrl0DzDhyTIU0DvR+4MPVk9AyHNF8rHgo7e0jzq5vP+BJyVw2uaQ4xtEqTqqk9fQi2qfLtC4cfra9VFhIy2fdBDpG3s5bcbtROjpRW9rSOtKWCN5++21uu+02PvjgA9avX084HGbp0qV09/lxePjhh3nkkUd47LHH2LRpEyUlJVx00UV0JdjvZm7ZxgkpJ5uwBprDGu0vJthsaIFAxEpobkPt7BpxyVmA5k86hnyv6HSO2DKI59MZCEGW+0cM3R5EpQhrF8h+FXliZSTHDECUkGZOQfD60XpTPVQNoSAP3+R8ZJ8Cu2o48vfEnCMaqkWWkgwSZXOfUCrYZrNhixPt/etf/xrz99NPP01RURFbtmzh3HPPRdM0Hn30UX70ox9x5ZVXAvDss89SXFzM888/zy233JKgBZkW0rihdHQivL8dsbENJAkpJxvR6USeWIng8UXaVfc0ZRwK8bp2VC4uIqvSOaTnj1gZTZ+CFgwN+f541p7S2ETu3w5S+FEHOVsa0bqOp3UIkoTqsKAUZKFOmkBgwVTarphN44WlSAEV4b1tKB2dTL2iDGuWZURr0A2DFPmvqKggOzs7eq1atWpI03Z2Rn4I8vIin7GamhoaGhpYunRp9B6bzcZ5553Hxo0bE7UawLSQxp1ea0Hp3Sa0Hy+zIVisJ93O9SWehZIzORNv8xCjTiMsh6I5rBAKIWZkjMq6UBqbkHx+lCkViC3H16KFgoif1kSd4PbcHKzNWXTOykF873i97MK5OcM+i2Q0amtryeqT5xfPOjoRTdO48847Ofvss6murgagoaEBgOLi4ph7i4uLOXz4cAIlNhXSuCHa7ZEqiSfZc49muwaRKNuQndrDUUaihJSZESmOtv8IiGJEGY3yQKE6rQKpvZvwCVuLvtab6vVC3TGypDloffpmSxYRVRWQKyYQrqsf81pTqchgUbasrKwYhTQUvv3tb/PJJ5/w7rvv9p/3hDLGmqb1Gxst5pZtnFD9/kEP3gmyHMntcrkGvOdkDRZFeeihcCk/Dyk3d0j3oioobjdCZgZaMHRcYYwyvKxt2UX44KFBo4rBL5yGZjve5kiQeipM5hdAKGRIZQQkPNv/O9/5Dq+99hpvvfUWEyYcr81eUhLpkddrKfXS1NTUz2oaLaZCSiG0cDgm4iNmZEQUx4ypiKfOikSr1IE/bEffa8bbPLSIkZDlQmlvjxTXH0LXVzEjAy0wOgtuIE4WVZQKC2mZa4muW55YCRrUrG/AX3MsugWW8vOM1ygyQWiaxre//W3++Mc/8ve//51Jk2JblE+aNImSkhLWr18fHQsGg7z99tssWrQoobKYW7YURO3qivQzC4YQc3PA40XwioSbWhBOUlepaVvHkF8jfPhoz4upqEOJmCkKKApanHbco6G3aH/cx6ZN5shVJfgLVHylDlxlpWgOG5rWP5fNP38S9u2SsbL4E5TLdtttt/H888/z6quv4nK5opZQdnY2DocDQRC44447WLlyJdOmTWPatGmsXLkSp9PJddddl4iVRDEVUorS608KH62L8dVomjCgQ3nujZPYu/Yogc4hRMF6tzmShDihFOXAIUSbDdXvR7TbETKcMY7zXme7lJMNgojSxxk/GpSurpj1aGeeSiDfRvN8GcWqoVk0rB0imfvbI9E4vx+Ly8rMZWXseLYGqbCQ4JwKbBs/HZfT8ClFghTSr371KwAWL14cM/70009z0003AfDDH/4Qn8/HrbfeSnt7O6effjrr1q3DdRL3wkgwFVIaIDqdCFZrRAlo2oDRraxK57DLj2jhMKIvgDh3RmTgk89QT5mGt9wZtzSs0tE5pC3e0AXQIjWjepRS+0wnLacrCAEVZ71IxjGNnM+6UHbtQcrJRmltw5pvJbPcgZSfR2hGOZYOf8wcJsNDG0oZZkFgxYoVrFixYkxlMRVSGqB2d0N3N4IsI2Znofn8qD5fzC+hIIIgDi+doveAY7juGOrkYqz76yEjg5bqTPx5Ak5RAk1FdDiOO7IFIeGVL0WHAzHLhTChFE+FgBAWKNwskPPr2DMugsMBHZ1ILieaKqBWlhLKsmBpigQLBEk66RZQd6jE/wFK36KhpkJKJ7RweMAT0qIciU8Mp1Fk37msR9vQXBl4F1aiWEEMQevNn8d1JIRzVz1aMBhxPo9F4qYkEW5sRijJp/CTMP5siZzffIBUXITW6T5+NstqiRwkDTWjBMKw+xDOo1lRv5FhFFEPekyuNRWSDhBsNtRwiE9fOoISHNnPY7jmMFJhIZnbAgRdE3Ad9qFaRMSgAoIwYCSsNw0mZsxiHfRMlWCzIRXkg6JEImWihPbxLpw7ZBw9851oiWkWGeXAIXwZFg68fgwtEEAz8hbNLNBmkmqILhdCcQHK/hpa93YPvxaOKCFIEogCSnMzgiyTvS8XbdMOek/+hImvZKSpk1BdzpjaRTC0A55aKByx+HqVjqbGKD4pJzvir8rNjfjOBAHB40WaPoXwnv207Ym8Rtpm6pvExTyHlM6IPSpDFHGWZbHgtinDn0NV0ELByPknUUKcWIHU3L9/WYySEYRITluGHTEYjpwN6stQzgOpSowFJNpsyOVlEasJEDIyEO12lI6OHgEiZX2VPfvJmZLJrGsr40xqMFRt4CtNMRVSOqMqqN1etCN1iIKCbB/l26mpCN0+1KYWRLsdafoU5PKymFtElwtpykTULAfa7gMQVgiX5BxXjjDkLYNgsR6ftyCf4KQiUFQEizVaP1uKE1aWHRIWp9Rv3HCYfdlMUg5VQQ2oCIJ9xAXLovQUT9MUhY6r5yMokHk0gOzuQu32IhXmIwgCWlc32v6ayHMamxEyJ6D1nGsaiv8Ijh+IFDMyEJxOwkfrsFothGZXIby3LbI0vx/iJBuPpjibSWpjWkh6QNOQszJRQ+qQy5eciDyhPOKPysjAt/RUVElA9mvIW/YgWC3IleXgyiDc2BTdaol2OxQXIh1rBSLKSMzKHNLrqT0pMoLVCj3/HT5UixAa/CT4cHL29M1A1lH6KmvTQtIJ7v2t7A0LkS/4EEuY9EXz+VC7uvCdO5P26TLl/70dtbsbobgoEknz9ESz+mwHVL8f9tVET31roSBKe3//E0SUV9/SKoLVihYKR7aJTgeipkWSjz/aMaisrZ+5cR82cHStFzPKZpKqKCEVz8E2xOzhlZvoRfP5ESxW/LkSivV4JUWlufV4mkmPs7o3AgZEH5OyshDyctC6PHHPSglTqhBraqNRMbWrKzKfIILTgTqMPnQhT5iQZ2Rlfk1SG3PLlmSknOyEzFM4N4fqW2YQnlGBNPV4tvZQG1SKWRHnsWIRkP0gzJ+DVFgYW9qj55c3qoz6oLjdKPWNCJkZ8V+gth71lKmxY5qG0tGB5o4tyzLYtrNycRFTLis76T2GwIyymSSaeF/ukSDKAmFZAg3U2mOR0PyMqYhTqob0/HBjE1JRAc4WheyDYaSG1uOWjigNKZSvBQJoNiuCzYYwfw5yacnxBx12QlnWfs8RZAtqZ1dMDajQmbP7HyU4Ya1aGn/pEoamDnylKeaWLc3p9c3IWU60sIblSAvhQADl/M8RtEk4DrYhl5YQro+UlIh3shoiCbzh+kbsf4qE28N9q0EOsQCaVFyEmmFHdGWiSgLhxuOlQFovnoLs1+itgi0V5KO0tKKFQ5E62j3F68RTZ9FeZSP3rSOx8vVJnDWd2vrFtJDSGNHlijqKhVAA2n2E6xsRXS7CTomuChk1y4HSp2Y182ZGCr31zuF0IhUWonZ3I+WObvuoNDbB7v0IrkzEmmNRRSZXVdBdJhDM7GNlKSpiRmR7p4XD0XNMzafl4CsUUM7/XKyF1beNksUM+wPmOSSTFEIQ0Hy+6J/HPmpF2NIWKbjW1YXtz5uwnD0P99RMXFuVSKKqpxtfoQNLRhlyd0/KhaJEwu6iFN2iiS7XiPvca4EA4ZrDyCXF0XmUwmwCeRqCKqCePQ/LzppoPSVhwRykY62ohTkIIYVgtoAmgbvSRkFdBtT3LNdqQaqYirJnP4febBxSyQzdow4Q4k/j7aypkNIVTYvZeomyiOhwoPgjLYUEixXLnjpsmRNBVSJhdsWNxRsmlCGjnDkD1SLgPORG2bUn2qdNkOURK6MY8VwZiIKA2u3FPTUT2SOgytA+w4E0aTbOphDOzxoJWyS6Tq9EUCDoErG1aVh8Gs7GIMreA9GDlkqnG9yRtUkF+YTqDVQZciDMsL9JSiJKVJ5fghoKc2h95EurhYIozc3Y3mxHA8K1RxHtduQPP8VaWkyoLJdgtoVgUQbWpkLUHuf6SLvmnoiy72DEMb1wNrlvH8LWUUHbTAvuqWDpEvDnWbG25yAdqMcml4MkYPEIoIG0YWt0HkES0UL0fPkUpJxsppztoG1fLg2bh9as0iR9MBWSHlAVREkj7Iv9ZRRPnYW2+7iV0etvUmsOI9QcxtkT1lda2wbt/DEStM4uZFFELc7DfqwLZuahymDxgBTQIlUpP2pGPEkd7BP71CkdnQhSLprp1I7s1uJaSOMuScIwFZJOEGWxX082dfunsVE1UULMcEa3ZL2FzaTCwjEpjq+0t0N7e6T2UW4OZX+HlgW55O06fiJbyslGyM4ifLgWGDgK2BfTqd2DuWUzSVUEIX61SC0cPt4kQFXi+ofGQhlF6xjR4+huaISGRgqP5Ued59ET333OYg2lq4komQpJr5gKSSfs+ePRgR8cp19MQZYRZ0whnGWPnBPaFFFIvRaYYLNFersJIoJVjn8odAiybv+fg4kWPT1RVeIW0FbTdztrnkPSCRkl9qTXCBKsVjRZRAyrSLVNiE5n5IGebH4tEIj0nBOFmIaYwyWjxI5kNT+6ejyHZL6rOmHKpWVkTxogj2ycUP0B1O2fIh5pRPP70YKRukgnFt8fbSRv9teqcBYNLUfPJL0wt2w6QbQI/Zza407PyexEt0k6EbNAWw+mU9skVUmH/K6hRNCGQmSt6fulSxg6PKltbtl0Qqg7TNif2gpJC4cTUm4l5A2jjrDdk0lqY1pIOmHHMzXJFqEfUmEhypRSFLuMtGErgsWakHIrm3+5d/TC6QBNU9HilBqJN5YumBaSTiiszk65yFN4ahmBfBuaFMny7y3+LxUWjnhOQRTIn50Vv4W00dAGKM6Wxj6k1PoEm4yY6csmIKdYayDh/e04/rYD+W9bjhdgEwQESUSQZaSsrGFv4SS7yKyrKxFEUyOZYX+TlEQQI5ZDKjp6o/WaSosi9Y16Gj5q4TA47DDMHDpJjnxkNSX11moyekwfkg4Qe76k8VJHUgVl7wFEuz2m+4ggioR7jwj0lsgd5NddkIWUjyaOG6oKQpx/C9OHZJJUBOg81J38c0iDoPr9CH26omi9ReJg6FsNDdxHvIPfZwTMLZtJKqIEVHY8W5Pyhe/FjAwEuy3qT1J9/iF3RenF3x5k568PjYF0JqnAqBTSqlWrEASBO+64I+7jt9xyC4Ig8Oijjw4616OPPsqMGTNwOBxUVFTwve99D3+fWjg33XQTgiDw0EMPxTxv7dq1CEPoiKFnZKdE+Zn5yRZjUNTubsKHa6M92bRQcNg5bbYsC8XzcsZGwDRDU9UBr3RlxApp06ZNPPHEE5xyyilxH1+7di0ffvghZWWD98/67W9/yz333MMDDzzAp59+ylNPPcVLL73EvffeG3Of3W5n9erVtPeUtTCJYMuyMOHskYfSxwspJxsEAbliAtqiU6OF/YeDo9BG+aKCMZAuDTG3bBE8Hg/Lly/nySefJDc3t9/jdXV1fPvb3+a3v/0tFoslzgyxvP/++5x11llcd911TJw4kaVLl/K1r32NzZs3x9y3ZMkSSkpKWLVq1UjE1i3pktuldHRGomy1RxHe2zbk9kp9SZe1moyMESmk2267jcsuu4wlS5b0e0xVVa6//np+8IMfMGfOnCHNd/bZZ7NlyxY++ugjAA4ePMjrr7/OZZddFnOfJEmsXLmSNWvWcPToSer/9CEQCOB2u2MugElfKGHSxSVUXViM1RUJNpaelkfRvBwK5mSTN8OFIAkIkkBGqR1HgQ1bjgVLRmqd9QEQLWJKR9jkqorj/z1paI0rB8JUSH3QYefaYYf9X3zxRbZu3cqmTZviPr569WpkWea73/3ukOf86le/SnNzM2effTaaphEOh/nXf/1X7rnnnn73Llu2jHnz5vHAAw/w1FNPDTr3qlWrePDBB/uNh7wKssOCbBcRRAFBhNwpmYgWMXLJAu7DXiSbSPXyiYgWAVEW0TSNjf++m+yqDKpvmIgSUlFDGr7WADueqSF3aiYV5xSihjWUkEp3g58jG5rImZxJ7tRM1HDkfk+Dj/Z9HlzlDmw51ui4ry1AoCOENcuCIIIa1lBDKmpIHTCaG+oO0/KpO/6DKUDH6eVk9pSoDdccHtVcwa4wHQc9iRAr/dE04hZoS+Mt27AUUm1tLbfffjvr1q3DHqf/+pYtW/jlL3/J1q1bh+Vo3rBhAz/72c94/PHHOf3009m/fz+33347paWl3H///f3uX716NRdccAF33XXXoHPfe++93HnnndG/3W43FRUVHH2nGVmI3U7ufvHIiU8n7Ff48D8+i/7de0LYfdTL5jV7kXqUV+9HwNsc4NimNiRZQLSIBLsj2e2apoEAslNGlARCvsh2JavSSf6sbMSe+xu2tHHsg1Ymf6GEvOlZiD1pF/v/r46GLe2ccvNkHHnWqKLa96c63Ie9ZBTamLO8CjWkoYRVDr3ZSKg7TNUFxRFl13N//aY2RFkgd6orqgSVkEpXrRfRImDNtMTcnwhrRAqoMb3eepsOjAT3ES/uI95IN96GxrT+8pn0R9CG0XFv7dq1LFu2DEk6vm1RFAVBEBBFkdWrV/ODH/wAURRjHhdFkYqKCg4dOhR33nPOOYczzjiDn//859Gx3/zmN/zzP/8zHo8HURS56aab6OjoYO3atQBcdtllWCwWbrrpJpYtWzbkxoFut5vs7GwW86V+CiklEUCyiKiKhqZo2HIsyHYposBkke5GP458KyUL83Af8SLKApJFpGFrO2pIZeKSkui9okVg7ytHsWTKzLiyomdcQFPh4//eT87kTGZfV4koRd6/7iY/H/9qP0Xzcph8cWlUgXXUeNj/p2MUz8uhcG5ORHmFVdr3e2j8uJ2izxWSUZGJ0tkN2Tk0ixn4XvmYnMkZyA4ZHBlogkTXp/UEPWHseVYAtHBEmSp+ZUBrMHtiBrZsC22tGYTrjqWFQgprITbwKp2dnWRlZQ3+hEHo/QyfL38l7mc4rIV4K/z7hL3eeDIsC+nCCy9kx44dMWPf+MY3mDlzJnfffTelpaVcfPHFMY9ffPHFXH/99XzjG98YcF6v1xujxCDiL9I0bUBF89BDDzFv3jymT58+nCWkHxoofUptBDpCBAjF3OIssmPNkGnc2j/6ePCv9f3GAh0hPvnf/nWpOw562Pjvu0GI+Gp6rcHWT914jvkiSk0WovJ46v1oWkd0ixvCHjlxbbMiOqxYHDIEfWS2KviA3KkusiZnI8oighrmcJeT1s/czLyqAmeRLXrifMezNXQe6ubMe2cBPdvWsMYnTx8kZ04RJbNtFDcHUMNVaGGVT39XizVDpuK8QtRQRDmG/QpH323BlmMhd6oLrcfqC3aH6azpxpopY8uxROcO+xVCnjCiHFmzqmip305IG6Cm9ghPaj/++OP8/Oc/p76+njlz5vDoo49yzjnnjE7GYTIsheRyuaiuro4Zy8jIID8/Pzqenx97HsZisVBSUsKMGTOiYzfccAPl5eXRaNkVV1zBI488wvz586Nbtvvvv58vfvGLMdZYX+bOncvy5ctZs2bNcJagSyRL/xZIo0KjZ77InEpAxdvU/7xQd6Of7sbjZ8Wk4iKkvFxaP/PSus+P0tIac//hd92o6xqOdxvpYduTB6L/LUhC9IDn1sf3R5WdKAsEu8Jo3m48DRpNH7dHt7maoqGpGsEehSJaRCw9ytSSIZM3LTOiBGWBQHuQzppucqZkUnVBcXTu9n0ePvt9LZWLi5hwVuQIhaqoHPuwjUPrG5j8hdIeH2BE4R37sJXmHZ1UXVh8fAsdVmnY0o7nmI/yM/MjwYawFjlPtyUxb01fNFVDE/q/7yNpM/7SSy9xxx138Pjjj3PWWWfx//7f/+OSSy5h9+7dVFZWJkLcIZGUXLYjR47EWEQ//vGPEQSBH//4x9TV1VFYWMgVV1zBz372s5PO89Of/pSXX355rMVNeURZQEmB/C6lsQkpPy/S9jpeSH9yJez87KQ1kfomzQY6Q/0eF0QBX2uA5p2ROXpPegc9YY681b90rqfOx+4X+vsGm7Z30LS9o9/44beaOPpeS3SbqwQj66jf1ErrHjeiFBn3NkcUseeYj1B3OGo99ubZyU4Za6aMKAvYNWlMFFIieeSRR7j55pv55je/CUQOKr/xxhv86le/GtdjNsPyIemBzs5OcnJyOJtLkUkDH9IQyKp0Ijsk2vb077mWcojSiM4f9ZI7NRMEaN+XPpG2MCHe5XU6OjrIzh59xcxeH9JAn+He16utrY3xIdlsNmxxUnWCwSBOp5Pf/e53LFu2LDp+++23s23bNt5+++1RyzxUDJft39UT6XmX15MsSQLpbwCkLqM15PYnRIqk0NXVlRCFZLVaKSkp4d2GgT/DmZmZVFRUxIw98MADrFixot+9LS0tKIpCcXFxzHhxcTENDQ2jlnc4GE4hlZWVUVtbi8vl0k0OXO9RhhN/EfVIOq5V0zS6urqGlEY1FOx2OzU1NQSDAx+d0DSt3+c7nnXUlxPvjzfHWGM4hSSKIhMmTEi2GGNCVlZW2nxJR0u6rTURllFf7HZ73LOAI6GgoABJkvpZQ01NTf2sprHGLD9iYmJwrFYrCxYsYP369THj69evZ9GiReMqi+EsJBMTk/7ceeedXH/99SxcuJAzzzyTJ554giNHjvAv//Iv4yqHqZB0gM1m44EHHhjUR6AHjLTW8eTaa6+ltbWVn/zkJ9TX11NdXc3rr79OVdXokqGHi+HC/iYmJqmL6UMyMTFJGUyFZGJikjKYCsnExCRlMBWSiYlJymAqJBMTk5TBVEgpSF1dHV//+tfJz8/H6XQyb948tmwZOF38j3/8IxdddBGFhYVkZWVx5pln8sYbb8Tcs2LFCgRB6HfNnDlzrJdzUoa71r689957yLLMvHnzYsZTda0mg2OeQ0ox2tvbOeusszj//PP5y1/+QlFREQcOHCAnJ2fA57zzzjtcdNFFrFy5kpycHJ5++mmuuOIKPvzwQ+bPnx+9b86cObz55psxz5Xl5H0ERrLWXjo7O7nhhhu48MILaWxs7Pd4qq3VZGiY71CKsXr1aioqKnj66aejYxMnTjzpc05sxLly5UpeffVV/vSnP8UoJFmWKSkpSaS4o2Ika+3llltu4brrrkOSpGhZ476k2lpNhoa5ZUsxXnvtNRYuXMjVV19NUVER8+fP58knnxzWHKqq0tXVRV5e3hhJmRhGutann36aAwcO8MADD4yDlCbjiamQUoyDBw/yq1/9imnTpvHGG2/wL//yL3z3u9/lueeeG/Ic//mf/0l3dzfXXHNNzPiOHTvIzMyMuXorBCaDkax137593HPPPfz2t7896RYs1dZqMjTMLVuKoaoqCxcuZOXKlQDMnz+fXbt28atf/Yobbrhh0Oe/8MILrFixgldffZWioqKYx2bMmMFrr70WM+ZyuRIn/DAZ7loVReG6667jwQcfHLS5Q6qt1WRomAopxSgtLWX27NkxY7NmzeIPf/jDoM996aWXuPnmm/nd734Xt6uw1Wpl6tSpCZN1tAx3rV1dXWzevJmPP/6Yb3/720BEqWmahizLrFu3jgsuuABIvbWaDA1TIaUYZ511Fnv27IkZ27t376BZ1y+88AL/9E//xAsvvNCvBXmqMty1ZmVl9WvD9fjjj/P3v/+d3//+90yaNGnMZDUZH0yFlGJ873vfY9GiRaxcuZJrrrmGjz76iCeeeIInnngies+9995LXV1d1NfywgsvcMMNN/DLX/6SM844I1r5z+FwxFQqDIfD/aoCCoIw7lUBexnuWkVR7NeGq6ioCLvd3m881dZqMkQ0k5TjT3/6k1ZdXa3ZbDZt5syZ2hNPPBHz+I033qidd9550b/PO++83iZqMdeNN94YveeBBx6Ie4/NZhunVcVnuGs9kQceeEA79dRT+42l4lpNBsesh2RiYpIymGF/ExOTlMFUSCYmJimDqZBMTExSBlMhmZiYpAymQjIxMUkZTIVkYmKSMpgKycTEJGUwFZKJiUnKYCokExOTlMFUSCYmJimDqZBMTExShv8fXYYF9RhiHtUAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "ds.attrs['pyproj_srs'] = 'EPSG:32632'\n", "ds['thk'].attrs['pyproj_srs'] = 'EPSG:32632'\n", "\n", "ds['thk'].salem.quick_map(interp='linear');" ] }, { "cell_type": "code", "execution_count": 8, "id": "75344953-0f73-4418-9f5c-c4643bb55e01", "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "
<xarray.Dataset>\n",
       "Dimensions:      (y: 994, x: 423)\n",
       "Coordinates:\n",
       "  * y            (y) float32 4.963e+06 4.963e+06 ... 5.062e+06 5.062e+06\n",
       "  * x            (x) float32 2.643e+05 2.644e+05 ... 3.064e+05 3.065e+05\n",
       "Data variables: (12/15)\n",
       "    thkobs       (y, x) float32 ...\n",
       "    thkinit      (y, x) float32 ...\n",
       "    lon          (y, x) float32 ...\n",
       "    lat          (y, x) float32 ...\n",
       "    usurfobs     (y, x) float32 ...\n",
       "    uvelsurfobs  (y, x) float32 nan nan nan nan nan nan ... nan nan nan nan nan\n",
       "    ...           ...\n",
       "    arrhenius    (y, x) float32 ...\n",
       "    slidingco    (y, x) float32 ...\n",
       "    thk          (y, x) float32 ...\n",
       "    usurf        (y, x) float32 ...\n",
       "    icemask      (y, x) float32 ...\n",
       "    divflux      (y, x) float32 ...\n",
       "Attributes:\n",
       "    pyproj_srs:  EPSG:32632
" ], "text/plain": [ "\n", "Dimensions: (y: 994, x: 423)\n", "Coordinates:\n", " * y (y) float32 4.963e+06 4.963e+06 ... 5.062e+06 5.062e+06\n", " * x (x) float32 2.643e+05 2.644e+05 ... 3.064e+05 3.065e+05\n", "Data variables: (12/15)\n", " thkobs (y, x) float32 ...\n", " thkinit (y, x) float32 ...\n", " lon (y, x) float32 ...\n", " lat (y, x) float32 ...\n", " usurfobs (y, x) float32 ...\n", " uvelsurfobs (y, x) float32 nan nan nan nan nan nan ... nan nan nan nan nan\n", " ... ...\n", " arrhenius (y, x) float32 ...\n", " slidingco (y, x) float32 ...\n", " thk (y, x) float32 ...\n", " usurf (y, x) float32 ...\n", " icemask (y, x) float32 ...\n", " divflux (y, x) float32 ...\n", "Attributes:\n", " pyproj_srs: EPSG:32632" ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" } ], "source": [ "ds" ] }, { "cell_type": "code", "execution_count": null, "id": "5377503c-8f4d-43d0-8681-a5e7fb8804af", "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "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.10.8" } }, "nbformat": 4, "nbformat_minor": 5 }