{"nbformat":4,"nbformat_minor":0,"metadata":{"colab":{"provenance":[]},"kernelspec":{"name":"python3","display_name":"Python 3"},"language_info":{"name":"python"}},"cells":[{"cell_type":"code","execution_count":null,"metadata":{"id":"MoIIbrsjbMJw"},"outputs":[],"source":["import numpy as np\n","import pandas as pd\n","import matplotlib.pyplot as plt\n","import xarray as xr"]},{"cell_type":"markdown","source":["# Loading Argo netCDF file"],"metadata":{"id":"CoVV9G9kjvmo"}},{"cell_type":"code","source":["filepath = '/content/5906248_Sprof.nc'\n","\n","data = xr.open_dataset(filepath)\n","display(data)"],"metadata":{"colab":{"base_uri":"https://localhost:8080/","height":360},"id":"YYsmw0Bmjyh4","executionInfo":{"status":"ok","timestamp":1692831478024,"user_tz":420,"elapsed":1218,"user":{"displayName":"Josie Adams","userId":"12688257352029271630"}},"outputId":"549b1934-576d-4f92-c48a-776b70b7c9c2"},"execution_count":null,"outputs":[{"output_type":"display_data","data":{"text/plain":["\n","Dimensions: (N_PROF: 23, N_PARAM: 8, N_CALIB: 1,\n"," N_LEVELS: 559)\n","Dimensions without coordinates: N_PROF, N_PARAM, N_CALIB, N_LEVELS\n","Data variables: (12/86)\n"," DATA_TYPE object ...\n"," FORMAT_VERSION object ...\n"," HANDBOOK_VERSION object ...\n"," REFERENCE_DATE_TIME object ...\n"," DATE_CREATION object ...\n"," DATE_UPDATE object ...\n"," ... ...\n"," NITRATE (N_PROF, N_LEVELS) float32 ...\n"," NITRATE_QC (N_PROF, N_LEVELS) object ...\n"," NITRATE_dPRES (N_PROF, N_LEVELS) float32 ...\n"," NITRATE_ADJUSTED (N_PROF, N_LEVELS) float32 ...\n"," NITRATE_ADJUSTED_QC (N_PROF, N_LEVELS) object ...\n"," NITRATE_ADJUSTED_ERROR (N_PROF, N_LEVELS) float32 ...\n","Attributes:\n"," title: Argo float vertical profile\n"," institution: AOML\n"," source: Argo float\n"," history: 2023-06-30T02:40:35Z creation (software version 1.1...\n"," references: http://www.argodatamgt.org/Documentation\n"," user_manual_version: 1.0\n"," Conventions: Argo-3.1 CF-1.6\n"," featureType: trajectoryProfile\n"," software_version: 1.15 (version 09.06.2022 for ARGO_simplified_profile)"],"text/html":["
\n","\n","\n","\n","\n","\n","\n","\n","\n","\n","\n","\n","\n","\n","\n","
<xarray.Dataset>\n","Dimensions:                          (N_PROF: 23, N_PARAM: 8, N_CALIB: 1,\n","                                      N_LEVELS: 559)\n","Dimensions without coordinates: N_PROF, N_PARAM, N_CALIB, N_LEVELS\n","Data variables: (12/86)\n","    DATA_TYPE                        object ...\n","    FORMAT_VERSION                   object ...\n","    HANDBOOK_VERSION                 object ...\n","    REFERENCE_DATE_TIME              object ...\n","    DATE_CREATION                    object ...\n","    DATE_UPDATE                      object ...\n","    ...                               ...\n","    NITRATE                          (N_PROF, N_LEVELS) float32 ...\n","    NITRATE_QC                       (N_PROF, N_LEVELS) object ...\n","    NITRATE_dPRES                    (N_PROF, N_LEVELS) float32 ...\n","    NITRATE_ADJUSTED                 (N_PROF, N_LEVELS) float32 ...\n","    NITRATE_ADJUSTED_QC              (N_PROF, N_LEVELS) object ...\n","    NITRATE_ADJUSTED_ERROR           (N_PROF, N_LEVELS) float32 ...\n","Attributes:\n","    title:                Argo float vertical profile\n","    institution:          AOML\n","    source:               Argo float\n","    history:              2023-06-30T02:40:35Z creation (software version 1.1...\n","    references:           http://www.argodatamgt.org/Documentation\n","    user_manual_version:  1.0\n","    Conventions:          Argo-3.1 CF-1.6\n","    featureType:          trajectoryProfile\n","    software_version:     1.15 (version 09.06.2022 for ARGO_simplified_profile)
"]},"metadata":{}}]},{"cell_type":"code","source":["pres_first_prof = data['PRES'].values[0,:]\n","temp_first_prof = data['TEMP'].values[0,:]"],"metadata":{"id":"47y1QKJmmTW8"},"execution_count":null,"outputs":[]},{"cell_type":"code","source":["plt.plot(temp_first_prof,pres_first_prof)\n","plt.gca().invert_yaxis()"],"metadata":{"colab":{"base_uri":"https://localhost:8080/","height":429},"id":"LTrdYyQJoWvQ","executionInfo":{"status":"ok","timestamp":1688587094315,"user_tz":420,"elapsed":664,"user":{"displayName":"Ethan C Campbell","userId":"11255944928409084259"}},"outputId":"d23b582b-5bdc-4cf5-b33f-ba0598c92a09"},"execution_count":null,"outputs":[{"output_type":"display_data","data":{"text/plain":["
"],"image/png":"iVBORw0KGgoAAAANSUhEUgAAAjAAAAGdCAYAAAAMm0nCAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAABH9UlEQVR4nO3deXhTVf4G8DdLk65JKaVNCwUKSNnLJrUKCIIUZHDfABEVwQVUwEHs/Bw2ZwTBQZ0ZRgcVdBwUl1FUXMtalLIVSqFAZSm00A0oTbpmPb8/0gYCLbSl6c3yfp7nPk3uPTf55hrJ+9x7zj0yIYQAERERkQeRS10AERERUWMxwBAREZHHYYAhIiIij8MAQ0RERB6HAYaIiIg8DgMMEREReRwGGCIiIvI4DDBERETkcZRSF+AqNpsN+fn5CAkJgUwmk7ocIiIiagAhBMrKyhAdHQ25vP7zLF4bYPLz8xETEyN1GURERNQEeXl5aNeuXb3bvTbAhISEALAfAI1GI3E1RERE1BAGgwExMTGO3/H6eG2Aqb1spNFoGGCIiIg8zLW6f7ATLxEREXkcBhgiIiLyOAwwRERE5HEYYIiIiMjjMMAQERGRx2GAISIiIo/DAENEREQehwGGiIiIPA4DDBEREXkctw4wK1asQMeOHeHv74+EhATs2rVL6pKIiIjIDbhtgPnss88we/ZszJ8/H3v37kV8fDySkpJQXFwsdWlEREQkMbcNMMuXL8fUqVPx+OOPo0ePHnj33XcRGBiIVatWSV0aERERScwtA4zJZEJ6ejpGjhzpWCeXyzFy5EikpaXVuY/RaITBYHBafIWh2oyF32Xh1fWHpC6FiIioRbhlgDl37hysVisiIyOd1kdGRqKwsLDOfRYvXgytVutYYmJiWqJUt1BlsmL1byfxwa85OHhGL3U5RERELueWAaYpkpOTodfrHUteXp7UJbWY8GC14/GX6aclrISIiKhlKKUuoC7h4eFQKBQoKipyWl9UVASdTlfnPmq1Gmq1us5t3u69bSccj/vGhEpXCBERUQtxyzMwKpUKAwYMwMaNGx3rbDYbNm7ciMTERAkrc0//3nrc8Xh0r7oDHhERkTdxywADALNnz8Z7772Hjz76CIcPH8YzzzyDiooKPP7441KX5nZu7drG8fjjtFMSVkJERNQy3PISEgA89NBDOHv2LObNm4fCwkL07dsXP/300xUde31ZudGCvacuIFB98T+jn0ImYUVEREQtQyaEEFIX4QoGgwFarRZ6vR4ajUbqcprNqfMVWLMzFztPnMfBfAOstov/+RJiw/DvSQMQGqiSsEIiIqKma+jvt9uegaErVZutGPVmKowWm2Ndu1YBGBQbhgcHxiAhNgwyGc/AEBGR92OA8RCbjhThqY/TYbbaz7j8cVRX3NO/HdqGBkhcGRERUctjgPEAn+/Jw0tfZgIAFHIZXhhxA2bcdoPEVREREUmHAcbNCSEc4WXIDeFYcl8fnnUhIiKfxwDjZmw2gQNn9Dh5vgKnL1Rh85GLs2+P7qVjeCEiIgIDjFs4W2ZEVr4emaf1WJ+Zj9+Lyp22+/vJMWVwLB4Y4DvzOxEREV0NA4yLWKw2VFtsqDJZUW22L1VmK6pM9r/VZhsqTRas3Z2HXTklTvsG+CnQu50WMa0CERseiPsHxECn9ZfokxAREbkfBphmYLUJHMo3YGfOeew4UYI9p0pQWmlu9OuM6hGJ27pFIKmnDq2CeC8XIiKi+jDAXAd9lRk/HCjA2xuOotBQXW+7AD8F/P3k9r8qRc1z+99IjT8iNGoMj4vAoNiwFqyeiIjIczHANEGBvgrvbjmOT3fnwVRzU7kQtRKDYsOQ0CkMg2Jbo0NYIAJUCqiVct5cjoiIqJkxwDRS2vHzmLx6lyO4dG4ThNG9dJg+vAsCVTycRERELYG/uI30ZfppR3h5++G+uDM+mmdYiIiIWhgDTCOF+F88ZC+szcBfvz+M0EA/hAaoEKFRo0PrQLQPC4ROG4BIjRoRIf5oFejHkENERNSMGGAa6d7+bfHh9pOO58VlRhSXGa+6j59CBm2AHzQBfggN8IM2wA+hgSqEBamgq+nEGxHiD02AEkEqJQLVCgSplAjwU0AuZ/AhIiK6nEwIIaQuwhUaOh13UwghUGGyorTSBH2VGfpKMy5UmpFfWoVTJRXILalCsaEaxWVGlFSYmvw+MhkQ6KdAoFqJQJUCSrkMSrkcCrnMsSgveXzxuRwKORxtlXIZ5Je0dX4ud3rup5DDTyGDSimHSiGHn0IOlbL2rwwqhcKxvXabSiGHn1KOQD8Fgv2V8FPIm/FoExGRL2no7zfPwDSBTCZDsFqJYLUS7Vpdva3JYsP5CiNKK832sFMTeEqrTDhfbkKhoRpFNWGnvNqCSpMVFSYLhACEACpMVlSYrC3zwZqJSilHiFqJYH/7GaVgfyXUNYGnNiQpFfYwplLa/yoVMsd2pUIOv9q/CtnFdYqLbeUyexir/auQ4+JjmT2QOW2XySCXw3l77Wtcuv2ydZe+hlwGXgokInITDDAuplLKEaUNQJS24XMYCSFQbbahwmRBhdGCCqMVlSYLLDYBm03AYhOwXvLX/tjmWHdlGxusNsBqszmtd25n399sFTBbbTBZbDBbbTDW/DVbhWOd6ZLttetNFvt6oCa0WUw4fx1nn9yVXIY6gpHz2S2lI4hdPMOlrA1icrnjsUJuD2WKmlB3cf+adnIZFAoZ/GrOpPkpnPepfS+nbfIrX682zCnqCGVOwU8uuyLgyS/bj0GOiNwFA4wbkslkCFApEKBSIDxYLXU5DWax2lBhtKLMaEaF0Ypyoxll1RaUGy0wWWywWAXMNhvMFnuQMlsFLFYbzDb7X/s6ezuLzXbldquAyWqDTdhDl80GWGsfi4thzrFd4Ip1Vpu4uE/N40tf51psArBZBQCvvPLaYFcLcs5ntwCVQg7/mps3+vvJ4a+0P1b71axX2tcHqZXo2DoIXSOD0SZEjWC1EkpejiSiejDAULNRKuTQBsqhDfSTupQms10acBoclC45k2W1P7ZYa85o1ZwBM1svnvmqDWsWq/M2s80Gq/XiPhfb2INd7Rkyq80e6qw1Qa/2vZ3bXHxuD2kXz7rZPwsu+Xx1fOZr5LOWCnIBNf2qQvyVCFErEeLvh2C1/bl9vZ/jcmVI7XN/JTSXPA7wU/CMEZEXYoAhuoRcLoMcMvgppK5EWkKIS8IP6jhrdemZLMBiqz0zBuegJ+yXGKtrJjA1WqyOx46/NesMVRYcP1uO48XlKDNaAMA+AarZirPXGOl3NTIZHJfhGnL5znHZ75LntZfmnC7vXXK58NLLeMp63ufSS37Kq+0rl0HtZ7+Lt38dfxUcmUgEgAGGiOogk9X8CEv0/iaLDRVGC8qqLSgzmlFebXFcjiyrNqPMaLlinaHmeVm12fHXVtMZ3mS1AVYAjZ9j1e0o5TLH5Ti1UmE/41Rzi4bQQPstGrSO2zXY71FV+1gbaD9jxTNS5A0YYIjI7aiUcqiUquualV0I4RjVZ6njEl5dl90sNpvjMuClzx2X+ByXCWteo55Lfpe/Tu37XP7+ta/ldOnPZoPRbO9AbzRbYbykgzwAWGwC5UYLypt4Ukohl0Hjr3QEndZB9ptwtgnxd9x8MyJEjQiNGuHBat4WgdwWAwwReSWZTIYgtRJBas//Z85qEzBarDDWXHIzOi692VBWbXa6TUNppcnxvPSS2zboq8yoNtvD0oWae1ddi0wGtApUoXWQCq2DVWgdrEZ4kAphQWq0DlYhvGadfbsaGn+e3aGW4/n/ZxMReTmFXIZAlRKBTT8hBQCoNltrQk5N0Kky43y5CcVl9ntRFRuM9scGI86WG2G1CZRUmFBSYcLR4mu/vp9ChtZBakSH+qNdq0C0axWAsCBVTcdqv4udrAP8EKiyj0BT+8mhVsoZfKjRGGCIiHxE7XD2SI3/NdvabAIllfYbbp4vN+Jchf3v+XITzlcYca5mfUmFvU2Z0QKzVaDQUI1CQzX25pY2uC6ZDI5OyrXD6u1D7RUI8JPXjDyzh5/ax8H+SrQOUiFK64/o0AC0CVZz6hUfwwBDRERXkMtlCA9W19yLKuSa7avNVpRUmHC2zIgzpVU4faESZy5UQV9V28HaDEPVxQ7XlSaLY7i+EKgZmWZDU3tay2VAkNp+pufSofaXhp+LQ+4vtgm5rA3PBnkOBhgiIrpu/n4KRIcGIDo0APExoddsL4S943LtMHpjHUPrq802VJmtNSPSakeXXRxtdr7ChPzSKhQZqmETcGy7Hn4K+1Qx2gA/6LT+iNYGoG2rACR2ao04XQg0AX7s2OwmOJkjERF5NIvVhpIK+2Wssmr7EPtyo/1MT/klQ+3LjZYr11Vb7MPyjfY56BoiSKWANsAPmprh6o5h6k7Pazo+13R4bhWo4j18GoiTORIRkU9QKuSI0Pgj4jpew2YTqDRba+4vZEZJhX2y3TOlVThWXI5tR885bqhYO8luvr664TXKZegepUGvtlp0bB2I9mGB6N+hVYP6I1HdGGCIiMjnyeX2S0fBaiV02rpDhcVqQ1m15eKQ9aqLw9cNtY9rhq1fqDQ7OjmXVplhsQkcOKPHgTN6x+sp5DLERYagXSv7Zaq2oQH2x6GBaNsqAK0C/dgf5yoYYIiIiBpAqZCjVVDjb7BosdpQoK/G/tOlyC4sw6nzlThWXI5DBQbHUpcAPwXatgpA5zZBuLFjGG7q1BrdozS8FFWDfWCIiIgkcPJcBY6fLceZ0iqcuVCF0zV/z5RW1Tv/V4i/EgmxYegaGYJ2rQIRExaAmFaBiA4NgErpHZ2LG/r7zQBDRETkZowWKwpKq3H6QhUO5uux88R57D55AeXGukdZyWRA29AAjO6pw3MjboA2wK+FK24+DDAMMERE5EUsVhsOFRiwK6cEuSWVyCupxOkLVci7UFlzDx27ntEa/OXuXujXvpWE1TYdAwwDDBER+QAhBM6Vm7DnZAmSvz6A0pp5rm7p0hozht+AxM6tJa6wcRr6+93sF8wWL16MG2+8ESEhIYiIiMDdd9+N7OxspzbDhg2DTCZzWp5++mmnNrm5uRg7diwCAwMRERGBOXPmwGK5vhsUEREReRuZTIY2IWqM6R2Fb6cPxgMD2kEpl+G3Y+cx/r0deHTVLuSXVkldZrNr9gCzdetWTJ8+HTt27EBKSgrMZjNGjRqFiooKp3ZTp05FQUGBY1m6dKljm9VqxdixY2EymbB9+3Z89NFH+PDDDzFv3rzmLpeIiMhrtG8diGUPxGPLnGF45Kb2UMplSP39LKZ8tAfedsHF5ZeQzp49i4iICGzduhVDhw4FYD8D07dvX7z11lt17vPjjz/iD3/4A/Lz8xEZGQkAePfddzF37lycPXsWKtW1h7DxEhIREfm63SdL8MC7aQCA7lEaPH5zR9zZNxr+fgqJK6ufZJeQLqfX22/aExYW5rR+zZo1CA8PR69evZCcnIzKykrHtrS0NPTu3dsRXgAgKSkJBoMBWVlZri6ZiIjIK9zYMQzPDOuMAD8FDhcY8NL/MjH49U3IyCuVurTr5tIb2dlsNsycORO33HILevXq5Vg/YcIEdOjQAdHR0cjMzMTcuXORnZ2Nr776CgBQWFjoFF4AOJ4XFhbW+V5GoxFG48Vx8wZD3TcGIiIi8iVzR3fD00M7Y+3uXHzwaw6Ky4x49r/p2PrScI+emNKlAWb69Ok4ePAgfv31V6f106ZNczzu3bs3oqKiMGLECBw/fhydO3du0nstXrwYCxcuvK56iYiIvJE20A9P3doZ9/Rvi8FLNiNfX437303DZ9NucuvLSVfjsug1Y8YMrF+/Hps3b0a7du2u2jYhIQEAcOzYMQCATqdDUVGRU5va5zqdrs7XSE5Ohl6vdyx5eXnX+xGIiIi8SkSIP958qC9C1ErszyvFjE/24UKFSeqymqTZA4wQAjNmzMDXX3+NTZs2ITY29pr7ZGRkAACioqIAAImJiThw4ACKi4sdbVJSUqDRaNCjR486X0OtVkOj0TgtRERE5GxsnyismNgfSrkMGw4X4fY3U/FNxhmYLLZr7+xGmn0U0rPPPotPPvkE33zzDeLi4hzrtVotAgICcPz4cXzyySe444470Lp1a2RmZmLWrFlo164dtm7dCsA+jLpv376Ijo7G0qVLUVhYiEmTJuHJJ5/Ea6+91qA6OAqJiIiofvvzSvHHL/bjaHE5AECn8ceTQ2IxLj4akZq6Z+RuCZLdibe+qb9Xr16Nxx57DHl5eXjkkUdw8OBBVFRUICYmBvfccw9eeeUVp0JPnTqFZ555Blu2bEFQUBAmT56MJUuWQKlsWLcdBhgiIqKrqzZb8e7W4/ho+0lcqLmDb5sQNX58YQjCg9WS1MSpBBhgiIiIGuRChQmf7MrF8pTfYbUJ9Gsfis+fSpRklJLb3AeGiIiI3FurIBWmD++Cb6bfAn8/OfblluKJD3ej0uS+U/gwwBAREREAoFdbLVZOGgi1Uo5tR8/hg205UpdULwYYIiIichjatQ0m3dQBAFBgqJa4mvoxwBAREZGT0EA/AMDpC1VuOwkkAwwRERE5GRYXAQBI/f0shizdjJe+3I+DZ/QSV+WMAYaIiIic9GqrxfxxPRCsVuL0hSp8vuc07ntnOw4XuM88gwwwREREdIXHb4nFzj+NwOrHb0R8Oy2MFhumfLgbZ8uM1965BTDAEBERUZ2C1EoMj4vAqsduRHiwCvn6asz75qDUZQFggCEiIqJraB2sRny7UACAXF73HfdbGgMMERERXZUQArkllQCAiBBpphi4HAMMERERXdWeUxcckz72jQmVtpgaDDBERER0VVbbxXvBLPnxCJb/kg2jxSphRQwwREREdA0JsWF4ZWx3aAP8UKCvxt83HcOKTcckrYkBhoiIiK5KJpPhySGdsOv/RiCxU2sAwLZj5yStiQGGiIiIGkStVGBMbx0AYF9uKQ6clu7uvAwwRERE1GB/6BPteBzsr5SsDgYYIiIiahCL1YbxK3cAAMKCVIgND5KsFgYYIiIiapDSKjOyi8oAAC+O6ippLQwwRERE1CDhwWp0aB0IAIgI8Ze0FgYYIiIiahB9lRl5NXfklfLyEcAAQ0RERA1UbrTAJgCVUo4uEcGS1sIAQ0RERA2irJnI0WSxocJokbQWBhgiIiJqkIgQNdRKe3Q4V26UtBYGGCIiImqQs2VGWGrmRVIrFZLWwgBDREREDfL1vjOw2gTiY0Kh03IUEhEREbm5arMVy1N+BwDc2rWNxNUwwBAREVEDHC0qh9FiAwDc37+dxNUwwBAREVEDrN2dCwDo3z4U7WtuZiclBhgiIiK6pl05JQCA8YPaS1yJHQMMERERXZNN2EcftQ+T/uwLwABDREREHogBhoiIiK7JUG2/86685m68UmOAISIioqs6VlyGs2X2O+92aSPtHEi1GGCIiIjoqr7PLAQADItrg1ZBKomrsWOAISIionoJIfBzlj3A9I0JlbaYSzDAEBERUb0OnjHgUIEBAPDAwBiJq7mo2QPMggULIJPJnJZu3bo5tldXV2P69Olo3bo1goODcd9996GoqMjpNXJzczF27FgEBgYiIiICc+bMgcUi7bTdREREvujrfWcAAKN76tA2NEDiai5SuuJFe/bsiQ0bNlx8E+XFt5k1axa+//57fPHFF9BqtZgxYwbuvfde/PbbbwAAq9WKsWPHQqfTYfv27SgoKMCjjz4KPz8/vPbaa64ol4iIiOqxM+c8AGCoG8x/dCmXBBilUgmdTnfFer1ejw8++ACffPIJbrvtNgDA6tWr0b17d+zYsQM33XQTfvnlFxw6dAgbNmxAZGQk+vbti1dffRVz587FggULoFK5R+chIiIib/f+thPIyjfATyHDbd0ipC7HiUv6wBw9ehTR0dHo1KkTJk6ciNxc+/wJ6enpMJvNGDlypKNtt27d0L59e6SlpQEA0tLS0Lt3b0RGRjraJCUlwWAwICsrq973NBqNMBgMTgsRERE13btbTwAAnhzSCTqtv8TVOGv2AJOQkIAPP/wQP/30E9555x3k5ORgyJAhKCsrQ2FhIVQqFUJDQ532iYyMRGGhvYdzYWGhU3ip3V67rT6LFy+GVqt1LDEx7tPRiIiIyNNUmiwwVJkBAONvdI/5jy7V7JeQxowZ43jcp08fJCQkoEOHDvj8888REOC6zj/JycmYPXu247nBYGCIISIiaqLPd+fBZLWhfVgg2rZyn867tVw+jDo0NBRdu3bFsWPHoNPpYDKZUFpa6tSmqKjI0WdGp9NdMSqp9nld/WpqqdVqaDQap4WIiIgar9psxXvbcgAAU4d2gsJNpg+4lMsDTHl5OY4fP46oqCgMGDAAfn5+2Lhxo2N7dnY2cnNzkZiYCABITEzEgQMHUFxc7GiTkpICjUaDHj16uLpcIiIin7cy9QTOlFYhPFiF+/u3k7qcOjX7JaQ//vGPGDduHDp06ID8/HzMnz8fCoUC48ePh1arxZQpUzB79myEhYVBo9HgueeeQ2JiIm666SYAwKhRo9CjRw9MmjQJS5cuRWFhIV555RVMnz4darW6ucslIiKiy+w5dQEA8MTgWASoFBJXU7dmDzCnT5/G+PHjcf78ebRp0waDBw/Gjh070KaNffz4m2++Cblcjvvuuw9GoxFJSUn417/+5dhfoVBg/fr1eOaZZ5CYmIigoCBMnjwZixYtau5SiYiI6DKGajNSfz8LAOgX00riauonE0IIqYtwBYPBAK1WC71ez/4wREREDXS0qAy3v5kKAMhamIQgtUtuGVevhv5+cy4kIiIicjh+thwAEOKvhL+fe14+AhhgiIiI6BL5pdUAgH7tW7nl6KNaDDBERETkYKi237xO49+yl44aiwGGiIiIHHJLKgEAVpt7d5FlgCEiIiIA9ukDduWUAAD6t3ffEUgAAwwRERHVSP7qAE5fqIJKKcdt3d1r9unLMcAQERERjhQa8E1GPgDg3Uf6o3ObYIkrujoGGCIiIh9nstjwj03HAAAJsWG4rVukxBVdm3t3MSYiIiKXstoEnv90H37KKgQATBkcK3FFDcMAQ0RE5MOeXZOOn7OKoFLIsfT+PhjVUyd1SQ3CS0hEREQ+qtpsxc9ZRQCA1+7tjbv7tZW4ooZjgCEiIvJRaSfOAwBkMuBeDwovAAMMERGRz9p5wn7Pl0EdwyB342kD6sIAQ0RE5KOOFdsnbhzTyzP6vVyKAYaIiMhHnSs3AgAC1Z43pocBhoiIyAdVm63IyCsFANzYMUzaYpqAAYaIiMgHbf39LAD7rNPtwwIlrqbxGGCIiIh8jM0m8LdfsgEA4we1h8LDOvACvJEdERGRTympMOG+d7Yj51wF5DLgvgHtpC6pSXgGhoiIyId8tfc0cs5VAAAW3dULXSNDJK6oaRhgiIiIfMS3+/Pxl+8PAwBeGh2HR27qIHFFTccAQ0RE5AOqzVa8+HkGAGBEtwg8cYtnTNpYHwYYIiIiH/Bx2imYrQKtg1RYMbE//P0UUpd0XRhgiIiIvFyRoRp/33gUAPDHpDiPDy8AAwwREZHX+zL9NMqMFnSP0uC+/p456uhyDDBERERebn/NHXdHdo+ASukdP/28DwwREZGXEkJgzc5c/HKoCAAwpleUxBU1HwYYIiIiL2S0WPHMf/di05FiAMDkxA7oEa2RuKrmwwBDRETkZSpNFjz/6T5sOlIMtVKOmSO74qmhnaQuq1kxwBAREXmZP6/LwobDxVAp5Fj92I24uUu41CU1O+/oyUNEREQAgEP5Bvxv72nIZMC/Jvb3yvACMMAQERF5lbW7cwEAg7uEY2SPSImrcR0GGCIiIi+x9fez+E/aKQDA5MSO0hbjYgwwREREXsBqE3j5f5kAgEduau/VZ18AduIlIiLyeFUmK17/6QgK9NVQymVIHtNd6pJcjgGGiIjIg1ltAtM+3oNtR88BAB68MQZBau//eW/2S0gdO3aETCa7Ypk+fToAYNiwYVdse/rpp51eIzc3F2PHjkVgYCAiIiIwZ84cWCyW5i6ViIjI46UePYttR89BrZTjr/f0wsI7e0pdUoto9oi2e/duWK1Wx/ODBw/i9ttvxwMPPOBYN3XqVCxatMjxPDAw0PHYarVi7Nix0Ol02L59OwoKCvDoo4/Cz88Pr732WnOXS0RE5LEyT5fi+U/2AQBu7twaExM6SFxRy2n2ANOmTRun50uWLEHnzp1x6623OtYFBgZCp9PVuf8vv/yCQ4cOYcOGDYiMjETfvn3x6quvYu7cuViwYAFUKlVzl0xEROSR3t+WgzKj/QrF8yNukLialuXSUUgmkwn//e9/8cQTT0AmkznWr1mzBuHh4ejVqxeSk5NRWVnp2JaWlobevXsjMvJi7+mkpCQYDAZkZWXV+15GoxEGg8FpISIi8lZ7Tpbg2/35AIBPnkxAv/atJK6oZbm0l8+6detQWlqKxx57zLFuwoQJ6NChA6Kjo5GZmYm5c+ciOzsbX331FQCgsLDQKbwAcDwvLCys970WL16MhQsXNv+HICIicjNWm8DLXx0AANzbv63X3m33alwaYD744AOMGTMG0dHRjnXTpk1zPO7duzeioqIwYsQIHD9+HJ07d27yeyUnJ2P27NmO5waDATExMU1+PSIiInf13rYTOFZcDo2/En+6w/uHTNfFZQHm1KlT2LBhg+PMSn0SEhIAAMeOHUPnzp2h0+mwa9cupzZFRUUAUG+/GQBQq9VQq9XXWTUREZF7M1lsWJl6AgAwd0w3hAf75m+fy/rArF69GhERERg7duxV22VkZAAAoqKiAACJiYk4cOAAiouLHW1SUlKg0WjQo0cPV5VLRETkEXbmnEdJhQmBKgUeGui7VxpccgbGZrNh9erVmDx5MpTKi29x/PhxfPLJJ7jjjjvQunVrZGZmYtasWRg6dCj69OkDABg1ahR69OiBSZMmYenSpSgsLMQrr7yC6dOn8wwLERH5tGqzFYu+OwQAuKdfWygVvjsjkEsCzIYNG5Cbm4snnnjCab1KpcKGDRvw1ltvoaKiAjExMbjvvvvwyiuvONooFAqsX78ezzzzDBITExEUFITJkyc73TeGiIjIF/1941EcLS5HaKAfXvCxYdOXkwkhhNRFuILBYIBWq4Ver4dGo5G6HCIiouvy3f58PPep/aZ1bz/cF3f1bStxRa7R0N9v3z33RERE5EE+2ZkLABjatQ3ujI++RmvvxwBDRETk5nLOVSDtxHkAwKyRNzjdHNZXMcAQERG5MSEE3t9mHzYdHxPqc3fcrQ8DDBERkZsSQuAv3x/GmprLR8/f1kXiitwHAwwREZGb2n9ajw9+zQEALLqrJ0Z0j7zGHr6DAYaIiMgNWW0CX6bnAQASYsPwaGJHaQtyMy6dC4mIiIgar7TShGkfp2NXTgkA4O5+3jlk+nowwBAREbmR7/bnI/mrAyg3WhCoUmD68C4+PWVAfRhgiIiI3ITJYsOszzJgsQm0CVFj5aQBHHVUDwYYIiIiNyCEwB+/2A+LTUAuA36eORRhQSqpy3JbDDBEREQSs1hteHX9IXy7Px8yGfDeowMZXq6BAYaIiEhChfpqTHhvB06cqwAAPH/bDRwu3QAMMERERBLamXPeEV6WPxiPezjiqEEYYIiIiCQihMD6zAIAwMjukbi3fzuJK/IcvJEdERGRRFb/dhIph4oglwEzOE1AozDAEBERScBkseGNX7IBAC+P6Ya+MaHSFuRhGGCIiIhamBACT328B5UmKwL8FJgyuJPUJXkcBhgiIqIWJITAuowz2Jx9FgDwl7t7QSGXSVyV52EnXiIiohZy6nwFXll3ENuOngMADLkhHPcNYMfdpmCAISIiagFCCExetQsnz1dCpZDjicGx7Lh7HRhgiIiIXOxcuRH/2HgUJ89XAgDWPz8YXSNDJK7KszHAEBERuVC12Yoxb2/D2TIjAOChgTEML82AAYaIiMiFfjxY4Agv/5zQD6N76iSuyDswwBAREbnQ+9tyAABP3BKLP/SJlrga78Fh1ERERC5yrtyIrHwDAGBSYgeJq/EuDDBEREQu8mbK7wCA+JhQdGwdKHE13oUBhoiIyAV+OFCAL9JPAwCeuKUjZDLerK45sQ8MERFRM6swWvD8p/tgsQn0jQll3xcX4BkYIiKiZrbxSDEsNgE/hQxrp93EqQJcgAGGiIiomW2tmefonn5t4e+nkLga78QAQ0RE1IyqzVZ8u/8MAGBUD97zxVUYYIiIiJpRkaEaZqsAANwa10biarwXAwwREVEzemfLcQDALV1aw0/Bn1lX4ZElIiJqJv9JO4m1u/MAAM8O40zTrsRh1ERERM1gZepxvPbDEQDA2D5RuKVLuMQVeTeegSEiIrpOm7OLHeHl3v5tsfS+PhJX5P0aHWBSU1Mxbtw4REdHQyaTYd26dU7bhRCYN28eoqKiEBAQgJEjR+Lo0aNObUpKSjBx4kRoNBqEhoZiypQpKC8vd2qTmZmJIUOGwN/fHzExMVi6dGnjPx0REZGL5ZVUYubaDADAhIT2WP5gXwSpeYHD1RodYCoqKhAfH48VK1bUuX3p0qX4+9//jnfffRc7d+5EUFAQkpKSUF1d7WgzceJEZGVlISUlBevXr0dqaiqmTZvm2G4wGDBq1Ch06NAB6enpWLZsGRYsWICVK1c24SMSERG5RrXZiqf/mw59lRnxMaGYP66H1CX5DnEdAIivv/7a8dxmswmdTieWLVvmWFdaWirUarX49NNPhRBCHDp0SAAQu3fvdrT58ccfhUwmE2fOnBFCCPGvf/1LtGrVShiNRkebuXPniri4uAbXptfrBQCh1+ub+vGIiIjqZbPZxB8/zxAd5q4X/Rb9Is5cqJS6JK/Q0N/vZu0Dk5OTg8LCQowcOdKxTqvVIiEhAWlpaQCAtLQ0hIaGYuDAgY42I0eOhFwux86dOx1thg4dCpVK5WiTlJSE7OxsXLhwoc73NhqNMBgMTgsREZGrfJF+Gl+kn4ZcBvxjfD9EhwZIXZJPadYAU1hYCACIjIx0Wh8ZGenYVlhYiIiICKftSqUSYWFhTm3qeo1L3+NyixcvhlardSwxMTHX/4GIiIjqcKy4DPO/yQIAzL69K0ccScBrRiElJydDr9c7lry8PKlLIiIiL1RttmLGJ/tQZbZicJdw3u9FIs0aYHQ6+5wPRUVFTuuLiooc23Q6HYqLi522WywWlJSUOLWp6zUufY/LqdVqaDQap4WIiKi5/fX7wzhSWIbwYBWWPxQPOWealkSzBpjY2FjodDps3LjRsc5gMGDnzp1ITEwEACQmJqK0tBTp6emONps2bYLNZkNCQoKjTWpqKsxms6NNSkoK4uLi0KpVq+YsmYiIqMF+OliAj3ecAgD87cG+iAjxl7gi39XoAFNeXo6MjAxkZGQAsHfczcjIQG5uLmQyGWbOnIm//OUv+Pbbb3HgwAE8+uijiI6Oxt133w0A6N69O0aPHo2pU6di165d+O233zBjxgw8/PDDiI6OBgBMmDABKpUKU6ZMQVZWFj777DO8/fbbmD17drN9cCIiosY4faESL32ZCQB46tZOuLUrJ2qUVGOHN23evFkAuGKZPHmyEMI+rOzPf/6ziIyMFGq1WowYMUJkZ2c7vcb58+fF+PHjRXBwsNBoNOLxxx8XZWVlTm32798vBg8eLNRqtWjbtq1YsmRJo+rkMGoiImouJotV3LPiV9Fh7npx1z9/FSaLVeqSvFZDf79lQgghYX5yGYPBAK1WC71ez/4wRER0XZb9fAQrNh9HiFqJH14YgpiwQKlL8loN/f32mlFIRERErvDr0XP415bjAIAl9/VheHETDDBERET1OFtmxKzPMyAEMH5Qe4ztEyV1SVSDAYaIiKgONpvAi1/sx9kyI7pGBmPeHzjPkTthgCEiIqrDe9tOIPX3s/D3k+OfE/ojQKWQuiS6BAMMERHRZfblXsCyn7MBAPPH9UTXyBCJK6LLMcAQERFdQl9lxnOf7oPFJjC2TxQevpFz67kjBhgiIqIaQgj86asDOH2hCu1aBWDxvb0hk3GqAHfEAENERFRj7e48fH+gAEq5DP8Y3w8afz+pS6J6MMAQEREB+L2oDAu+zQIAzEmKQ7/2nHvPnTHAEBGRz6syWTHjk70wWmwY2rUNpg7pJHVJdA0MMERE5NOEEFi0Pgu/F5WjTYgayx+Mh1zOfi/ujgGGiIh82psbjuLTXXmQyYA3H+yL8GC11CVRAzDAEBGRz3prw+/4+8ajAIBXxvbA4BvCJa6IGooBhoiIfNLbG47irQ328PJ/d3THlMGxEldEjcEAQ0REPucfG4/izQ2/AwD+dEc3TB3KTruehgGGiIh8yorNx/C3FHt4eXlMN0wb2lniiqgpGGCIiMhnrNh8zDHH0Uuj4/D0rQwvnooBhoiIfMI7W447wsucpDg8O6yLxBXR9WCAISIir/fvrcfx+k9HAAB/HNUV04czvHg6BhgiIvJqK1OPY/GP9vAy+/aumHHbDRJXRM2BAYaIiLzW+9tO4LUf7OFl5sgb8PwIhhdvwQBDRERe6YNfc/CX7w8DAF4YcQNmjuwqcUXUnBhgiIjI66z6NQevrj8EAHj+ti6YOZJnXrwNAwwREXmVD3/LwaKa8DJjeBfMur0rZDJOzuhtGGCIiMhr/CftJBZ8Zw8vzw7rjBdHMbx4KwYYIiLyCh+nncS8b7IAAM8M64w5SXEML16MAYaIiDzef3ecwp9rwstTt3bCSwwvXo8BhoiIPNqanafwyrqDAIBpQzvh5dHdGF58AAMMERF5rE935eL/vraHl6lDYpE8huHFVzDAEBGRR/oy/TSSvzoAAJgyOBZ/uqM7w4sPYYAhIiKPs/PEeSR/lQkAePyWjnhlLMOLr2GAISIij5JXUoln1uyF2Sowtk8U5v2hB8OLD2KAISIij1FWbcaUj3ajpMKEPu20eOP+eIYXH8UAQ0REHsFqE3hhbQZ+LypHRIgaKycNRIBKIXVZJBEGGCIi8giv/3QEm44UQ62U471HB0Kn9Ze6JJJQowNMamoqxo0bh+joaMhkMqxbt86xzWw2Y+7cuejduzeCgoIQHR2NRx99FPn5+U6v0bFjR8hkMqdlyZIlTm0yMzMxZMgQ+Pv7IyYmBkuXLm3aJyQiIo/3xZ48rEw9AQB444F4xMeESlsQSa7RAaaiogLx8fFYsWLFFdsqKyuxd+9e/PnPf8bevXvx1VdfITs7G3feeecVbRctWoSCggLH8txzzzm2GQwGjBo1Ch06dEB6ejqWLVuGBQsWYOXKlY0tl4iIPNyekyX409f24dLP39YF4+KjJa6I3IGysTuMGTMGY8aMqXObVqtFSkqK07p//vOfGDRoEHJzc9G+fXvH+pCQEOh0ujpfZ82aNTCZTFi1ahVUKhV69uyJjIwMLF++HNOmTWtsyURE5KHySirx1MfpMFsFxvTSYebIrlKXRG7C5X1g9Ho9ZDIZQkNDndYvWbIErVu3Rr9+/bBs2TJYLBbHtrS0NAwdOhQqlcqxLikpCdnZ2bhw4YKrSyYiIjdQbrRg6n/24HyFCT2jNfjbg/GQyzniiOwafQamMaqrqzF37lyMHz8eGo3Gsf75559H//79ERYWhu3btyM5ORkFBQVYvnw5AKCwsBCxsbFOrxUZGenY1qpVqyvey2g0wmg0Op4bDAZXfCQiImoBNpvAzLUZOFJYhjYharw/eSACVS79ySIP47Jvg9lsxoMPPgghBN555x2nbbNnz3Y87tOnD1QqFZ566iksXrwYarW6Se+3ePFiLFy48LpqJiIi97D052xsOFwElVKOlZMGIEobIHVJ5GZccgmpNrycOnUKKSkpTmdf6pKQkACLxYKTJ08CAHQ6HYqKipza1D6vr99McnIy9Hq9Y8nLy7v+D0JERC3uf+mn8e7W4wCAZff3Qb/2V551J2r2AFMbXo4ePYoNGzagdevW19wnIyMDcrkcERERAIDExESkpqbCbDY72qSkpCAuLq7Oy0cAoFarodFonBYiIvIs6acuOCZonD68M+7q21biishdNfoSUnl5OY4dO+Z4npOTg4yMDISFhSEqKgr3338/9u7di/Xr18NqtaKwsBAAEBYWBpVKhbS0NOzcuRPDhw9HSEgI0tLSMGvWLDzyyCOOcDJhwgQsXLgQU6ZMwdy5c3Hw4EG8/fbbePPNN5vpYxMRkbs5U1qFpz7eA5PVhqSekXjx9jipSyI3JhNCiMbssGXLFgwfPvyK9ZMnT8aCBQuu6Hxba/PmzRg2bBj27t2LZ599FkeOHIHRaERsbCwmTZqE2bNnO/V/yczMxPTp07F7926Eh4fjueeew9y5cxtcp8FggFarhV6v59kYIiI3V2G04P5303C4wIDuURp8+XQigtTstOuLGvr73egA4ykYYIiIPIPNJvD0f9Pxy6EihAer8M2MwWgbyk67vqqhv9+cC4mIiCT1t5Rs/HKoCCqFHP+eNJDhhRqEAYaIiCSzbt8ZrNhsH3G05L7eGNCBI46oYRhgiIhIEgdO6/HS/zIBAE/f2hn39m8ncUXkSRhgiIioxemrzHj2k3SYLDaM6BaBOUkccUSNwwBDREQtSgiBl77cj7ySKrRrFYDlD/aFgnMcUSMxwBARUYta/dtJ/JxVBD+FDCsm9Ic20E/qksgDMcAQEVGL2Zd7Aa/9cBgA8H93dEd8TKi0BZHHYoAhIqIWUVppwoxP9sFiE7ijtw6Tb+4odUnkwRhgiIjI5YQQ+OMX+3GmtAodWgdiyX19IJOx3ws1HQMMERG53HvbTmDD4WKolHKsmNAfGn/2e6HrwwBDREQutedkCV7/KRsAMO8PPdCrrVbiisgbMMAQEZHLlFTY+71YbQLj4qMxMaG91CWRl2CAISIil7DZBGZ9loFCQzU6hQdh8b292e+Fmg0DDBERucQ7W49j6+9noVbKsWJifwSrlVKXRF6EAYaIiJrdzhPn8bdf7P1eFt3VE92jNBJXRN6GAYaIiJrVuXIjnvt0H2wCuLdfWzw4MEbqksgLMcAQEVGzsdoEZq7NQHGZEV0igvGXe3qx3wu5BAMMERE1m39uOoZfj51DgJ8C70zsj0AV+72QazDAEBFRs9h+7Bze2vg7AOAvd/fCDZEhEldE3owBhoiIrltxWTWeX5sBIYAHB7bDfQPaSV0SeTkGGCIiui5Wm8Dzn+7DuXIj4iJDsPDOXlKXRD6AAYaIiK7Lql9zsONECQJVCqyY2B8BKoXUJZEPYIAhIqImO1duxN83HgVgn+eoS0SwxBWRr2CAISKiJvvbL9koM1rQu62W93uhFsUAQ0RETZKVr8fa3XkAgHnjekAu5/1eqOUwwBARUaMJIbDou0MQAvhDnyjc2DFM6pLIxzDAEBFRo/10sBA7c0qgVsqRfEd3qcshH8QAQ0REjVJttuKvPxwGADw1tBPahgZIXBH5IgYYIiJqlA9+zcHpC1XQafzx9LDOUpdDPooBhoiIGqzYUI0Vm48BAOaOieNcRyQZBhgiImqwpT9no9JkRd+YUNwV31bqcsiHMcAQEVGDZJ4uxZfppwEA8zlsmiTGAENERNdUO2waAO7p1xb92reSuCLydQwwRER0Td9lFmDPqQsI8FNg7uhuUpdDxABDRERXV2WyYknNsOlnhnWGTusvcUVEDDBERHQNK1NPIF9fjbahAZg2tJPU5RABaEKASU1Nxbhx4xAdHQ2ZTIZ169Y5bX/ssccgk8mcltGjRzu1KSkpwcSJE6HRaBAaGoopU6agvLzcqU1mZiaGDBkCf39/xMTEYOnSpY3/dEREdF0K9FV4d+txAMDLY7rB308hcUVEdo0OMBUVFYiPj8eKFSvqbTN69GgUFBQ4lk8//dRp+8SJE5GVlYWUlBSsX78eqampmDZtmmO7wWDAqFGj0KFDB6Snp2PZsmVYsGABVq5c2dhyiYjoOrz+4xFUma24sWMr/KFPlNTlEDk0+g5EY8aMwZgxY67aRq1WQ6fT1bnt8OHD+Omnn7B7924MHDgQAPCPf/wDd9xxB9544w1ER0djzZo1MJlMWLVqFVQqFXr27ImMjAwsX77cKegQEZHr7M29gHUZ+ZDJgHl/6AmZjMOmyX24pA/Mli1bEBERgbi4ODzzzDM4f/68Y1taWhpCQ0Md4QUARo4cCblcjp07dzraDB06FCqVytEmKSkJ2dnZuHDhQp3vaTQaYTAYnBYiImoam+3isOn7+7dD73ZaiSsictbsAWb06NH4z3/+g40bN+L111/H1q1bMWbMGFitVgBAYWEhIiIinPZRKpUICwtDYWGho01kZKRTm9rntW0ut3jxYmi1WscSExPT3B+NiMhnbDxSjIy8UgSpFJgzOk7qcoiu0OyTWDz88MOOx71790afPn3QuXNnbNmyBSNGjGjut3NITk7G7NmzHc8NBgNDDBFRE72XegIAMCmxIyJCOGya3I/Lh1F36tQJ4eHhOHbMPvmXTqdDcXGxUxuLxYKSkhJHvxmdToeioiKnNrXP6+tbo1arodFonBYiImq8jLxS7DpZAqVchsdu7ih1OUR1cnmAOX36NM6fP4+oKHvv9cTERJSWliI9Pd3RZtOmTbDZbEhISHC0SU1NhdlsdrRJSUlBXFwcWrXi7auJiFzpvW32sy939o3mTevIbTU6wJSXlyMjIwMZGRkAgJycHGRkZCA3Nxfl5eWYM2cOduzYgZMnT2Ljxo2466670KVLFyQlJQEAunfvjtGjR2Pq1KnYtWsXfvvtN8yYMQMPP/wwoqOjAQATJkyASqXClClTkJWVhc8++wxvv/220yUiIiJqfnkllfjxQAEAYOoQ3rSO3FejA8yePXvQr18/9OvXDwAwe/Zs9OvXD/PmzYNCoUBmZibuvPNOdO3aFVOmTMGAAQOwbds2qNVqx2usWbMG3bp1w4gRI3DHHXdg8ODBTvd40Wq1+OWXX5CTk4MBAwbgxRdfxLx58ziEmojIxVb9lgObAIbcEI7uUbwUT+5LJoQQUhfhCgaDAVqtFnq9nv1hiIgaQF9pRuKSjag0WfGfJwZhaNc2UpdEPqihv9+cC4mIiAAAn+zKRaXJim66EAy5IVzqcoiuigGGiIhgstiw+rccAMCTQzrxrrvk9hhgiIgI3+7PR3GZEZEaNe6Mj5a6HKJrYoAhIvJxQgi8XzN0+rGbY6FS8qeB3B+/pUREPm7b0XM4UliGQJUCEwa1l7ocogZhgCEi8nG1N6576MYYaAP9JK6GqGEYYIiIfNjhAgO2HT0HuQx44pZYqcshajAGGCIiH1Z79mVM7yjEhAVKXA1RwzHAEBH5qEJ9Nb7bnw8AmMZpA8jDMMAQEfmoD7efhNkqMKhjGOJjQqUuh6hRGGCIiHxQudGCNTtPAQCmDuXZF/I8DDBERD5o1a85KKu2oFN4EEZ0i5C6HKJGY4AhIvIxGw4V4c0NvwMApg/vArmc0waQ52GAISLyIYcLDHhh7T4IAYwf1B739m8rdUlETcIAQ0TkI86WGfHkR3tQYbLi5s6tseiunpy0kTwWAwwRkQ+oNlsx7eM9OFNahU7hQXhn4gD4KfgTQJ6L314iIi8nhMBLX2ZiX24ptAF++OCxGzllAHk8BhgiIi/3943H8O3+fCjlMrzzSH/EhgdJXRLRdWOAISLyYt/tz3eMOPrL3b1wc+dwiSsiah4MMEREXiojrxR//GI/AODJwbF4eFB7iSsiaj4MMEREXii/tApPfrQHRosNI7pFIPmO7lKXRNSsGGCIiLxMhdGCKR/twblyI7rpQvD2+H5Q8GZ15GUYYIiIvIjVJvDC2gwcLjAgPFiF9ycPRLBaKXVZRM2OAYaIyIss/ekINhwugkopx8pHB6Jdq0CpSyJyCQYYIiIv8fnuPPw79QQAYNn9fdC/fSuJKyJyHQYYIiIvsOPEefzp6wMAgOdH3IC7+nKOI/JuDDBERB7u5LkKPP3fdFhsAmP7RGHmiBukLonI5RhgiIg8mL7SjCc+2o3SSjPiY0LxtwfiIeeII/IBDDBERB7KbLVh+id7ceJsBaK0/nhv0gD4+ymkLouoRTDAEBF5ICEEFn6XhV+PnUOgSoH3Jw9EhMZf6rKIWgwDDBGRB1qfWYD/7siFTAa89VBf9IzWSl0SUYtigCEi8jDnyo2Y/20WAOC5227AqJ46iSsiankMMEREHmb+N1koqTChmy4EM4Z3kbocIkkwwBAReZAfDhTg+wMFUMhleOOBeKiU/GecfBO/+UREHqKkwoQ/rzsIAHh2WGf0ast+L+S7GGCIiDzE/G+zcL7ChLjIEMy4jZeOyLc1OsCkpqZi3LhxiI6Ohkwmw7p165y2y2SyOpdly5Y52nTs2PGK7UuWLHF6nczMTAwZMgT+/v6IiYnB0qVLm/YJiYi8wE8HC/Hd/nwo5DIse6AP1Ere74V8W6MDTEVFBeLj47FixYo6txcUFDgtq1atgkwmw3333efUbtGiRU7tnnvuOcc2g8GAUaNGoUOHDkhPT8eyZcuwYMECrFy5srHlEhF5vAsVJrxSc+noqaGd0KddqLQFEbkBZWN3GDNmDMaMGVPvdp3OeTjfN998g+HDh6NTp05O60NCQq5oW2vNmjUwmUxYtWoVVCoVevbsiYyMDCxfvhzTpk1rbMlERB5t4XdZOFduRJeIYDzPeY6IALi4D0xRURG+//57TJky5YptS5YsQevWrdGvXz8sW7YMFovFsS0tLQ1Dhw6FSqVyrEtKSkJ2djYuXLhQ53sZjUYYDAanhYjI06UcKsK6jHzIZcCy+/twqgCiGo0+A9MYH330EUJCQnDvvfc6rX/++efRv39/hIWFYfv27UhOTkZBQQGWL18OACgsLERsbKzTPpGRkY5trVq1uuK9Fi9ejIULF7rokxARtTx9pRn/9/UBAMDUIZ3Qr/2V//YR+SqXBphVq1Zh4sSJ8Pd3np9j9uzZjsd9+vSBSqXCU089hcWLF0OtVjfpvZKTk51e12AwICYmpmmFExG5gUXrD6G4zIhObYIw6/auUpdD5FZcFmC2bduG7OxsfPbZZ9dsm5CQAIvFgpMnTyIuLg46nQ5FRUVObWqf19dvRq1WNzn8EBG5m01HivC/vachkwHL7o/npSOiy7isD8wHH3yAAQMGID4+/pptMzIyIJfLERERAQBITExEamoqzGazo01KSgri4uLqvHxERORN9FVm/Okr+6ijKbfEYkAH/rtHdLlGB5jy8nJkZGQgIyMDAJCTk4OMjAzk5uY62hgMBnzxxRd48sknr9g/LS0Nb731Fvbv348TJ05gzZo1mDVrFh555BFHOJkwYQJUKhWmTJmCrKwsfPbZZ3j77bedLhEREXmrv35/CIWGasSGB+HFUXFSl0Pklhp9CWnPnj0YPny443ltqJg8eTI+/PBDAMDatWshhMD48eOv2F+tVmPt2rVYsGABjEYjYmNjMWvWLKdwotVq8csvv2D69OkYMGAAwsPDMW/ePA6hJiKvtyW7GJ/vsV86Wnp/HwSoeOmIqC4yIYSQughXMBgM0Gq10Ov10Gg0UpdDRHRNhmozkt5MRYG+Go/f0hHzx/WUuiSiFtfQ32/OhURE5CYW/3AYBfpqtA8LxJwkXjoiuhoGGCIiN/Dr0XP4dFceAPulo0CVS+9yQeTxGGCIiCRmtQm8ss5+w7pHEzvgpk6tJa6IyP0xwBARSUwhl2HFxP64vUck5o7uJnU5RB6B5yiJiNxAz2gt3nt0oNRlEHkMnoEhIiIij8MAQ0RERB6HAYaIiIg8DgMMEREReRwGGCIiIvI4DDBERETkcRhgiIiIyOMwwBAREZHHYYAhIiIij8MAQ0RERB6HAYaIiIg8DgMMEREReRwGGCIiIvI4XjsbtRACAGAwGCSuhIiIiBqq9ne79ne8Pl4bYMrKygAAMTExEldCREREjVVWVgatVlvvdpm4VsTxUDabDfn5+QgJCYFMJmvW1zYYDIiJiUFeXh40Gk2zvrYv43F1DR5X1+BxdQ0eV9fwpOMqhEBZWRmio6Mhl9ff08Vrz8DI5XK0a9fOpe+h0Wjc/ovgiXhcXYPH1TV4XF2Dx9U1POW4Xu3MSy124iUiIiKPwwBDREREHocBpgnUajXmz58PtVotdSlehcfVNXhcXYPH1TV4XF3DG4+r13biJSIiIu/FMzBERETkcRhgiIiIyOMwwBAREZHHYYAhIiIij8MA0wAlJSWYOHEiNBoNQkNDMWXKFJSXl191n2HDhkEmkzktTz/9dAtV7L5WrFiBjh07wt/fHwkJCdi1a9dV23/xxRfo1q0b/P390bt3b/zwww8tVKlnacxx/fDDD6/4bvr7+7dgte4vNTUV48aNQ3R0NGQyGdatW3fNfbZs2YL+/ftDrVajS5cu+PDDD11epydq7LHdsmXLFd9XmUyGwsLClinYAyxevBg33ngjQkJCEBERgbvvvhvZ2dnX3M/T/31lgGmAiRMnIisrCykpKVi/fj1SU1Mxbdq0a+43depUFBQUOJalS5e2QLXu67PPPsPs2bMxf/587N27F/Hx8UhKSkJxcXGd7bdv347x48djypQp2LdvH+6++27cfffdOHjwYAtX7t4ae1wB+904L/1unjp1qgUrdn8VFRWIj4/HihUrGtQ+JycHY8eOxfDhw5GRkYGZM2fiySefxM8//+ziSj1PY49trezsbKfvbEREhIsq9Dxbt27F9OnTsWPHDqSkpMBsNmPUqFGoqKiodx+v+PdV0FUdOnRIABC7d+92rPvxxx+FTCYTZ86cqXe/W2+9VbzwwgstUKHnGDRokJg+fbrjudVqFdHR0WLx4sV1tn/wwQfF2LFjndYlJCSIp556yqV1eprGHtfVq1cLrVbbQtV5PgDi66+/vmqbl156SfTs2dNp3UMPPSSSkpJcWJnna8ix3bx5swAgLly40CI1eYPi4mIBQGzdurXeNt7w7yvPwFxDWloaQkNDMXDgQMe6kSNHQi6XY+fOnVfdd82aNQgPD0evXr2QnJyMyspKV5frtkwmE9LT0zFy5EjHOrlcjpEjRyItLa3OfdLS0pzaA0BSUlK97X1RU44rAJSXl6NDhw6IiYnBXXfdhaysrJYo12vxu+p6ffv2RVRUFG6//Xb89ttvUpfj1vR6PQAgLCys3jbe8J312skcm0thYeEVpyqVSiXCwsKueg12woQJ6NChA6Kjo5GZmYm5c+ciOzsbX331latLdkvnzp2D1WpFZGSk0/rIyEgcOXKkzn0KCwvrbM9r3xc15bjGxcVh1apV6NOnD/R6Pd544w3cfPPNyMrKcvkEqN6qvu+qwWBAVVUVAgICJKrM80VFReHdd9/FwIEDYTQa8f7772PYsGHYuXMn+vfvL3V5bsdms2HmzJm45ZZb0KtXr3rbecO/rz4bYF5++WW8/vrrV21z+PDhJr/+pX1kevfujaioKIwYMQLHjx9H586dm/y6RNcrMTERiYmJjuc333wzunfvjn//+9949dVXJayM6EpxcXGIi4tzPL/55ptx/PhxvPnmm/j4448lrMw9TZ8+HQcPHsSvv/4qdSku57MB5sUXX8Rjjz121TadOnWCTqe7ojOkxWJBSUkJdDpdg98vISEBAHDs2DGfDDDh4eFQKBQoKipyWl9UVFTvcdTpdI1q74uaclwv5+fnh379+uHYsWOuKNEn1Pdd1Wg0PPviAoMGDfKJH+jGmjFjhmOgybXOpnrDv68+2wemTZs26Nat21UXlUqFxMRElJaWIj093bHvpk2bYLPZHKGkITIyMgDYT4f6IpVKhQEDBmDjxo2OdTabDRs3bnQ6G3CpxMREp/YAkJKSUm97X9SU43o5q9WKAwcO+Ox3sznwu9qyMjIy+H29hBACM2bMwNdff41NmzYhNjb2mvt4xXdW6l7EnmD06NGiX79+YufOneLXX38VN9xwgxg/frxj++nTp0VcXJzYuXOnEEKIY8eOiUWLFok9e/aInJwc8c0334hOnTqJoUOHSvUR3MLatWuFWq0WH374oTh06JCYNm2aCA0NFYWFhUIIISZNmiRefvllR/vffvtNKJVK8cYbb4jDhw+L+fPnCz8/P3HgwAGpPoJbauxxXbhwofj555/F8ePHRXp6unj44YeFv7+/yMrKkuojuJ2ysjKxb98+sW/fPgFALF++XOzbt0+cOnVKCCHEyy+/LCZNmuRof+LECREYGCjmzJkjDh8+LFasWCEUCoX46aefpPoIbquxx/bNN98U69atE0ePHhUHDhwQL7zwgpDL5WLDhg1SfQS388wzzwitViu2bNkiCgoKHEtlZaWjjTf++8oA0wDnz58X48ePF8HBwUKj0YjHH39clJWVObbn5OQIAGLz5s1CCCFyc3PF0KFDRVhYmFCr1aJLly5izpw5Qq/XS/QJ3Mc//vEP0b59e6FSqcSgQYPEjh07HNtuvfVWMXnyZKf2n3/+uejatatQqVSiZ8+e4vvvv2/hij1DY47rzJkzHW0jIyPFHXfcIfbu3StB1e6rduju5UvtcZw8ebK49dZbr9inb9++QqVSiU6dOonVq1e3eN2eoLHH9vXXXxedO3cW/v7+IiwsTAwbNkxs2rRJmuLdVF3HE4DTd9Ab/32VCSFEC5/0ISIiIrouPtsHhoiIiDwXAwwRERF5HAYYIiIi8jgMMERERORxGGCIiIjI4zDAEBERkcdhgCEiIiKPwwBDREREHocBhoiIiDwOAwwRERF5HAYYIiIi8jgMMERERORx/h/gz4b3GrNk8gAAAABJRU5ErkJggg==\n"},"metadata":{}}]},{"cell_type":"code","source":["data['PRES'].attrs['long_name']"],"metadata":{"colab":{"base_uri":"https://localhost:8080/","height":34},"id":"pQzsEoSUovux","executionInfo":{"status":"ok","timestamp":1688587163704,"user_tz":420,"elapsed":12,"user":{"displayName":"Ethan C Campbell","userId":"11255944928409084259"}},"outputId":"cfbe5d32-f01a-44bf-bc92-8dbef59bb745"},"execution_count":null,"outputs":[{"output_type":"execute_result","data":{"text/plain":["'Sea water pressure, equals 0 at sea-level'"],"application/vnd.google.colaboratory.intrinsic+json":{"type":"string"}},"metadata":{},"execution_count":68}]},{"cell_type":"markdown","source":["# Loading, processing, and plotting NB Palmer\n"," shipboard CNV file"],"metadata":{"id":"tMs6e0bNflcL"}},{"cell_type":"code","source":["filepath = 'dnbp2211_001.cnv'\n","\n","# Identify header line number (= line_idx - 1)\n","file_object = open(filepath,'r')\n","all_lines = file_object.readlines()\n","file_object.close()\n","header_names = []\n","for line_idx, line in enumerate(all_lines):\n"," if '*END*' in line:\n"," break\n"," elif ' name' in line:\n"," header_names.append(line.split('= ')[1].split(':')[0])\n","header_names.append('')\n","# print('First line of data after header:',line_idx, line)\n","\n","# Extract header names\n","for line_idx, line in enumerate(all_lines):\n"," if ' name' in line:\n"," header_names.append(line.split('= ')[1].split(':')[0])\n","header_names.append('')\n","\n","pd.read_csv(filepath,header=None,skiprows=line_idx+1,delim_whitespace=True)"],"metadata":{"colab":{"base_uri":"https://localhost:8080/","height":476},"id":"ojtjlcpHfo-V","executionInfo":{"status":"ok","timestamp":1688585652821,"user_tz":420,"elapsed":330,"user":{"displayName":"Ethan C Campbell","userId":"11255944928409084259"}},"outputId":"45f1327f-b7dc-466a-cd97-1b2e1e643d07"},"execution_count":null,"outputs":[{"output_type":"execute_result","data":{"text/plain":[" 0 1 2 3 4 5 6 7 \\\n","0 3.0 2.986 0.6714 28.697590 0.6729 28.699145 7.8876 7.9709 \n","1 4.0 3.961 0.6706 28.698075 0.6731 28.700361 7.8880 7.9623 \n","2 5.0 4.953 0.6693 28.697519 0.6709 28.699846 7.8909 7.9627 \n","3 6.0 5.945 0.6608 28.691212 0.6632 28.693173 7.8883 7.9726 \n","4 7.0 6.935 0.6685 28.698019 0.6711 28.700008 7.8871 7.9641 \n","... ... ... ... ... ... ... ... ... \n","2024 2027.0 1998.408 1.3263 30.845606 1.3290 30.847386 4.2602 4.3357 \n","2025 2028.0 1999.389 1.3256 30.845380 1.3281 30.847227 4.2630 4.3362 \n","2026 2029.0 2000.371 1.3249 30.845152 1.3273 30.846828 4.2646 4.3354 \n","2027 2030.0 2001.352 1.3245 30.845188 1.3271 30.846991 4.2660 4.3372 \n","2028 2031.0 2002.333 1.3241 30.845491 1.3269 30.847048 4.2608 4.3330 \n","\n"," 8 9 ... 15 16 17 18 19 \\\n","0 0.3025 97.2691 ... 105.930 694.320 16665 1027.1163 1027.1160 \n","1 0.3049 97.6027 ... 103.870 695.620 16696 1027.1214 1027.1215 \n","2 0.3116 97.7276 ... 104.170 696.479 16716 1027.1264 1027.1274 \n","3 0.3036 97.9193 ... 103.450 713.080 17115 1027.1322 1027.1321 \n","4 0.2961 97.9622 ... 103.960 720.073 17283 1027.1362 1027.1358 \n","... ... ... ... ... ... ... ... ... \n","2024 -0.0203 100.1124 ... 39.696 3345.650 80297 1037.1768 1037.1759 \n","2025 -0.0177 100.1098 ... 39.665 3348.430 80363 1037.1814 1037.1807 \n","2026 -0.0198 100.1064 ... 40.152 3356.138 80548 1037.1859 1037.1852 \n","2027 -0.0184 100.1108 ... 40.492 3384.184 81221 1037.1905 1037.1897 \n","2028 -0.0123 100.1082 ... 38.314 3463.743 83131 1037.1952 1037.1941 \n","\n"," 20 21 22 23 24 \n","0 33.8013 33.8016 7.93803 1 0.0 \n","1 33.8021 33.8024 7.93815 27 0.0 \n","2 33.8022 33.8036 7.93841 19 0.0 \n","3 33.8027 33.8027 7.94012 40 0.0 \n","4 33.8026 33.8023 7.93856 30 0.0 \n","... ... ... ... .. ... \n","2024 34.7218 34.7212 7.75619 34 0.0 \n","2025 34.7218 34.7213 7.75633 36 0.0 \n","2026 34.7217 34.7212 7.75646 33 0.0 \n","2027 34.7217 34.7212 7.75654 40 0.0 \n","2028 34.7220 34.7210 7.75661 27 0.0 \n","\n","[2029 rows x 25 columns]"],"text/html":["\n","
\n","
\n","
\n","\n","\n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n","
0123456789...15161718192021222324
03.02.9860.671428.6975900.672928.6991457.88767.97090.302597.2691...105.930694.320166651027.11631027.116033.801333.80167.9380310.0
14.03.9610.670628.6980750.673128.7003617.88807.96230.304997.6027...103.870695.620166961027.12141027.121533.802133.80247.93815270.0
25.04.9530.669328.6975190.670928.6998467.89097.96270.311697.7276...104.170696.479167161027.12641027.127433.802233.80367.93841190.0
36.05.9450.660828.6912120.663228.6931737.88837.97260.303697.9193...103.450713.080171151027.13221027.132133.802733.80277.94012400.0
47.06.9350.668528.6980190.671128.7000087.88717.96410.296197.9622...103.960720.073172831027.13621027.135833.802633.80237.93856300.0
..................................................................
20242027.01998.4081.326330.8456061.329030.8473864.26024.3357-0.0203100.1124...39.6963345.650802971037.17681037.175934.721834.72127.75619340.0
20252028.01999.3891.325630.8453801.328130.8472274.26304.3362-0.0177100.1098...39.6653348.430803631037.18141037.180734.721834.72137.75633360.0
20262029.02000.3711.324930.8451521.327330.8468284.26464.3354-0.0198100.1064...40.1523356.138805481037.18591037.185234.721734.72127.75646330.0
20272030.02001.3521.324530.8451881.327130.8469914.26604.3372-0.0184100.1108...40.4923384.184812211037.19051037.189734.721734.72127.75654400.0
20282031.02002.3331.324130.8454911.326930.8470484.26084.3330-0.0123100.1082...38.3143463.743831311037.19521037.194134.722034.72107.75661270.0
\n","

2029 rows × 25 columns

\n","
\n"," \n"," \n"," \n","\n"," \n","
\n","
\n"," "]},"metadata":{},"execution_count":39}]},{"cell_type":"markdown","source":["# Loading, processing, and plotting float EDF file"],"metadata":{"id":"cKHom7C6fin3"}},{"cell_type":"code","source":["filepath = '/content/1122.001.edf.txt'\n","\n","# Identify header line number (= line_idx - 1)\n","file_object = open(filepath,'r')\n","all_lines = file_object.readlines()\n","for line_idx, line in enumerate(all_lines):\n"," if line[0] != '$':\n"," break\n","print('First line of data after header:',line_idx, line)\n","\n","# Extract column names\n","col_names = all_lines[line_idx-1].split()\n","print('Header column names:',col_names,'\\n')\n","\n","# Load data\n","data = pd.read_csv(filepath,comment='$',delim_whitespace=True,na_values='*',\n"," header=None,names=col_names[2:])\n","display(data)"],"metadata":{"colab":{"base_uri":"https://localhost:8080/","height":477},"id":"cWZGx94hbZWg","executionInfo":{"status":"ok","timestamp":1688434472607,"user_tz":420,"elapsed":163,"user":{"displayName":"Ethan C Campbell","userId":"11255944928409084259"}},"outputId":"670d6bfa-41a9-40fb-8d9e-2ab204f7bc5c"},"execution_count":null,"outputs":[{"output_type":"stream","name":"stdout","text":["First line of data after header: 91 2.30 10.366 34.405 17.786 0.874579 406 1457 168 -0.902169 nan 10.366 26.4249 273.6 284.8 10.333 6.2 32\n","\n","Header column names: ['$', 'T', 'p', 't', 's', 'O2Phase', 'O2tV', 'fl', 'bb', 'cd', 'pH-Vrs', 'pH-Vk', 'theta', 'sigma', 'o2sat', '[O2]', 'T63', 'speed', 'n'] \n","\n"]},{"output_type":"display_data","data":{"text/plain":[" p t s O2Phase O2tV fl bb cd pH-Vrs \\\n","0 2.3 10.366 34.405 17.786 0.874579 406 1457 168 -0.902169 \n","1 4.0 10.247 34.407 17.779 0.877397 449 1414 175 -0.901938 \n","2 6.0 10.171 34.408 17.776 0.879281 582 1473 173 -0.901731 \n","3 7.9 10.141 34.408 17.784 0.880097 644 1505 172 -0.901716 \n","4 10.0 10.118 34.409 17.794 0.880719 946 1590 175 -0.901667 \n",".. ... ... ... ... ... ... ... ... ... \n","496 1299.8 3.092 34.401 24.207 1.087143 64 343 219 -0.918275 \n","497 1349.4 3.005 34.424 24.404 1.089901 63 343 220 -0.918504 \n","498 1399.7 2.921 34.451 24.573 1.092556 63 345 221 -0.918678 \n","499 1449.3 2.857 34.471 24.694 1.094618 65 342 225 -0.918760 \n","500 1499.4 2.768 34.496 24.096 1.097324 63 345 228 -0.918283 \n","\n"," pH-Vk theta sigma o2sat [O2] T63 speed n \n","0 NaN 10.366 26.4249 273.6 284.8 10.333 6.2 32 \n","1 NaN 10.247 26.4471 274.3 286.2 10.229 4.5 44 \n","2 NaN 10.170 26.4610 274.8 287.1 10.159 5.1 39 \n","3 NaN 10.140 26.4662 274.9 287.1 10.128 5.4 37 \n","4 NaN 10.117 26.4709 275.1 287.0 10.105 5.9 34 \n",".. ... ... ... ... ... ... ... .. \n","496 -0.828354 2.999 27.4070 321.4 177.2 3.089 NaN 1 \n","497 -0.828408 2.909 27.4336 322.0 174.0 3.002 NaN 1 \n","498 -0.828381 2.822 27.4629 322.5 171.4 2.919 NaN 1 \n","499 -0.828328 2.755 27.4849 322.9 169.6 2.854 NaN 1 \n","500 -0.827827 2.663 27.5129 323.5 183.2 2.770 NaN 1 \n","\n","[501 rows x 17 columns]"],"text/html":["\n","
\n","
\n","
\n","\n","\n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n"," \n","
ptsO2PhaseO2tVflbbcdpH-VrspH-Vkthetasigmao2sat[O2]T63speedn
02.310.36634.40517.7860.8745794061457168-0.902169NaN10.36626.4249273.6284.810.3336.232
14.010.24734.40717.7790.8773974491414175-0.901938NaN10.24726.4471274.3286.210.2294.544
26.010.17134.40817.7760.8792815821473173-0.901731NaN10.17026.4610274.8287.110.1595.139
37.910.14134.40817.7840.8800976441505172-0.901716NaN10.14026.4662274.9287.110.1285.437
410.010.11834.40917.7940.8807199461590175-0.901667NaN10.11726.4709275.1287.010.1055.934
......................................................
4961299.83.09234.40124.2071.08714364343219-0.918275-0.8283542.99927.4070321.4177.23.089NaN1
4971349.43.00534.42424.4041.08990163343220-0.918504-0.8284082.90927.4336322.0174.03.002NaN1
4981399.72.92134.45124.5731.09255663345221-0.918678-0.8283812.82227.4629322.5171.42.919NaN1
4991449.32.85734.47124.6941.09461865342225-0.918760-0.8283282.75527.4849322.9169.62.854NaN1
5001499.42.76834.49624.0961.09732463345228-0.918283-0.8278272.66327.5129323.5183.22.770NaN1
\n","

501 rows × 17 columns

\n","
\n"," \n"," \n"," \n","\n"," \n","
\n","
\n"," "]},"metadata":{}}]},{"cell_type":"code","source":["plt.plot(data['t'],data['p'],color='k')\n","plt.gca().invert_yaxis()\n","plt.grid()\n","plt.xlabel('Temperature (°C)')\n","plt.ylabel('Pressure (dbar)');"],"metadata":{"colab":{"base_uri":"https://localhost:8080/","height":448},"id":"ge-0k0Xmblj_","executionInfo":{"status":"ok","timestamp":1688434388813,"user_tz":420,"elapsed":623,"user":{"displayName":"Ethan C Campbell","userId":"11255944928409084259"}},"outputId":"930e1c1f-ab43-448b-fb6d-edeb337af7bb"},"execution_count":null,"outputs":[{"output_type":"display_data","data":{"text/plain":["
"],"image/png":"iVBORw0KGgoAAAANSUhEUgAAAkQAAAGwCAYAAABIC3rIAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAABiq0lEQVR4nO3dd1gUh/418LNUASmi0hRR1FgQewmW2Kgao8YkFlRiN8Eejb33EkWNJZbYArFdNWoigqiosYGKsWBHsYEaykpf2Hn/8GV+2WAMmF1my/k8zz6XnZldzle48Tg7RSYIggAiIiIiA2YkdQAiIiIiqbEQERERkcFjISIiIiKDx0JEREREBo+FiIiIiAweCxEREREZPBYiIiIiMngmUgfQBUqlEs+ePYO1tTVkMpnUcYiIiKgYBEHA69ev4eLiAiOjd+8DYiEqhmfPnsHV1VXqGERERPQeHj9+jMqVK79zGxaiYrC2tgbw5g/UxsZG4jRvKBQKREREwNfXF6amplLH0SjOqr8MaV7Oqr8MaV5dm1Uul8PV1VX8e/xdWIiKofBjMhsbG60qRJaWlrCxsdGJX8r/grPqL0Oal7PqL0OaV1dnLc7hLjyomoiIiAweCxEREREZPBYiIiIiMngsRERERGTwWIiIiIjI4LEQERERkcFjISIiIiKDx0JEREREBo+FiIiIiAweCxEREREZPIMqRGvWrEHVqlVRpkwZtGjRAhcvXpQ6EhEREWkBgylEu3btwrhx4zBz5kxcvnwZDRo0gJ+fH168eCF1NCIiIpKYwdzcdfny5RgyZAgGDBgAAFi/fj1+/fVX/Pjjj5g0aZLKtrm5ucjNzRWfy+VyAG9uaqdQKEov9DsU5tCWPJrEWfWXIc3LWfWXIc2r7lmVSiUyMzMhl8uRm5sLd3d3tbxvoZLklAmCIKj1u2uhvLw8WFpaYu/evejWrZu4PCgoCGlpafjll19Utp81axZmz55d5H3CwsJgaWmp6bhEREQ6SxAEJCcn4+HDh0hJSUFKSgpSU1ORmpqKzMxMZGVlISsrC9nZ2cjOzkZhDalYsSI2btyo1ixZWVno06cP0tPTYWNj885tDWIP0atXr1BQUABHR0eV5Y6Ojrh161aR7SdPnoxx48aJz+VyOVxdXeHr6/uvf6ClRaFQIDIyEj4+PjA1NZU6jkZxVv1lSPNyVv1lSPO+a9br169j48aNOHjwIJ4+fVqi9zU2NkbZsmXRqVMndcYVP+EpDoMoRCVlbm4Oc3PzIstNTU217pddGzNpCmfVX4Y0L2fVX4Y0799njYyMxMcff4y8vDxxfYMGDVClShU4OzuLj/Lly8PGxkblYW1tDQsLC8hkMo3kLC6DKEQVKlSAsbExkpOTVZYnJyfDyclJolRERES6r6CgAKNHj0ZeXh7at2+PCRMmoF27drCwsJA6WokYxFlmZmZmaNKkCaKiosRlSqUSUVFR8PLykjAZERGR7nr16hVGjRqF+Ph42NvbY//+/QgICNC5MgQYyB4iABg3bhyCgoLQtGlTNG/eHCEhIcjMzBTPOiMiIqL/U1BQgPT0dKSnpyMtLQ1paWl49eoVoqKicOrUKSQkJCAqKko8TmfBggWwtbWVOPX7M5hC1LNnT7x8+RIzZsxAUlISGjZsiPDw8CIHWhMREekSQRCQl5cnnrVV+PjrmVxve/y16Pz168LH69evi/X969Wrh1mzZqFHjx4anlSzDKYQAcCIESMwYsQIqWMQEZEeKNxjIpfLcefOHVhaWiIvL08sIoWnl/+1mOTl5SEvL0+8rt0/PUqyTX5+vkbntLS0hJ2dHezs7GBra4uCggJ8+OGHqFGjBmrWrIkOHTrAzMxMoxlKg0EVIiIiov8qISEB48aNw8GDB6FUKqWOo8LIyAgWFhbFetja2opF56+Pvy63tbVVKTsKhQK//fYbOnXqpHdn1LEQERERFYNCocC+ffvw9ddfIyUlBQBgbW0NS0tLCIKA8uXLw8rKCpaWlrC0tISFhYX4taWlJcqUKQMzMzPxlPW/fv1Pj+JsY2pqijJlysDCwgKmpqYaOX3dELAQERER/YO8vDxcuXIFv/76KzZt2oTnz58DAGrVqoV9+/ahbt26er3XxJCwEBEREf1/WVlZOHv2LE6ePIkzZ87g4sWLyM7OFtc7OjpiyJAhGD16NCpUqCBhUlI3FiIiIiIAd+/eRatWrfDy5UuV5fb29mjdujUCAwPRrVs3vTiAmIpiISIiIgJw8ODBImWoefPmaNq0KWxsbHDv3j388MMPKFu2LKytrWFra4vKlSvDxcVFosSkTixEREREALy9vVGpUiWVG5NevHgRFy9e/NfX2tjYoEaNGqhatSrc3Nzg5uaGBg0a4MMPP0SZMmU0GZvUhIWIiIgIQIMGDfD48WPcv38fsbGxSE1NxevXr9/6yMjIQGpqKh4/fgy5XA65XI7Lly/j8uXLRd537NixWLp0KYyNjSWYioqLhYiIiOj/k8lkqFGjBmrUqKGyPD8/Hy9fvsSLFy+KPO7cuYPz58/j2bNnb33PFStWwN/fH76+vqUxAr0nFiIiIiIA2dnZuHLlCi5evIjY2Fg8efJELD1//vnne7/vyJEj0apVKzUmJU1gISIiIoMgCAJSUlLw6NEjJCYmIjExUfz63r17uH79+jtvg2FkZISKFSvCwcFB5WFvb4+nT5+idevWqFChAsqVKyc+KlSoACMjo1Kckt4XCxEREem9gQMHYsuWLf+6naOjI1q0aIFmzZqhZs2aRYrP244D4oUZ9QMLERER6TWlUonQ0FCVZc2aNYObmxuqVKkinhXWqFEjuLq68tYXBoqFiIiI9JqRkRGWLFmCMWPGAABMTEzw66+//uMeHzJM/GCTiIj03sCBA2Fubg7gzRljDg4O4o1QZTIZHB0dERsbK3FKkhILERER6T1ra2tERkaiffv2sLW1BfDmIOtCL168wK+//ipVPNICLERERGQQ2rRpg+PHj2Pu3LlF1g0ePBiTJk2SIBVpCx5DREREekWhUODx48dISkrCy5cvxQsqvnz5EleuXEF0dDQAoHbt2oiLixM/SiPDxkJEREQ6JS8vD7GxsXj06BGeP38uPp48eYKEhAQ8efIESqXyne/RsGFDbNy4kWWIRCxERESk1fLy8rBnzx5cuXIF169fx9mzZ/H69et3vsbc3BwuLi7ihRQL/9fd3R3t27dHjRo1eHo9qWAhIiIirbZy5Up8++23KssqVqyIunXrwtnZWXy4uLigWrVqqFatGhwdHXmFaCoRFiIiItJqf/1Ya/To0ejXrx8aNWrEwkNqxd8mIiLSan379hW//uKLL9CkSROWIVI77iEiIiKtdvz4cfHr4OBgNG7cGJ6envDx8YGHh4eEyUifsBAREZFWa9CgAWrXro1bt24hLi4OcXFxAABjY2OcPHkSrVu3ljYg6QXucyQiIq1Ws2ZNxMfH48mTJ9izZw9mzJgBR0dHFBQUYN++fVLHIz3BQkRERDqhUqVK+OyzzzB+/Hh8+OGHAMCbs5La8CMzIiLSGnl5eZDL5SqP58+f49q1a+IjISFB3L5NmzYSpiV9wkJERESlRhAEREREICYmBnfv3sW9e/eQnJwslp/c3NxivU+VKlUwceJEdOnSRcOJyVCwEBERUam4cuUKxo4dK95L7F0sLS1hY2MDGxsbVKhQAR4eHvD09BQf5cuXL4XEZEhYiIiISKNiY2MxdepUREREAAAsLCzw2WefoVatWqhZsyYqVaoEW1tbsQCVLVsWJib864lKF3/jiIhII16/fo3o6Gj07NkTCoUCANC7d28sWrQIVapUkTgdkSoWIiIiUpuEhASMHDkSv//+O9LS0sTlDRo0wL59++Du7i5dOKJ3YCEiIqL/RKlUYuPGjTh8+DAiIiKQl5cnrqtYsSL8/PywZMkSODs7S5iS6N20+jpECxcuRLNmzWBtbQ0HBwd069YNt2/fVtkmJycHwcHBKF++PMqWLYsePXogOTlZZZvExER07twZlpaWcHBwwIQJE5Cfn1+aoxAR6a3w8HAMHz4chw8fRl5eHtq1a4d58+bh5cuXePHiBXbs2MEyRFpPqwtRdHQ0goODcf78eURGRkKhUMDX1xeZmZniNmPHjsWhQ4ewZ88eREdH49mzZ/j000/F9QUFBejcuTPy8vJw9uxZbNu2DVu3bsWMGTOkGImISK8IgoCzZ8+Kz69cuYKIiAjUq1cPtra2EiYjKhmt/sgsPDxc5fnWrVvh4OCAS5cu4aOPPkJ6ejo2b96MsLAwdOjQAQCwZcsW1KlTB+fPn8eHH36IiIgI3Lx5E8eOHYOjoyMaNmyIuXPnYuLEiZg1axbMzMyKfN/c3FyVa2HI5XIAgEKhEA8MlFphDm3Jo0mcVX8Z0rz6Nuv169cRGRmJAwcO4Ny5cwDe3HjVw8ND72b9N4Y0r67NWpKcMkEQBA1mUat79+6hZs2auHbtGurVq4fjx4+jY8eOSE1NhZ2dnbidm5sbxowZg7Fjx2LGjBk4ePCgeDNA4M1Bf+7u7rh8+TIaNWpU5PvMmjULs2fPLrI8LCwMlpaWmhiNiEgn5OTkYObMmSqHL5iZmcHf3x+BgYEwNzeXMB2RqqysLPTp0wfp6emwsbF557ZavYfor5RKJcaMGYNWrVqhXr16AICkpCSYmZmplCEAcHR0RFJSkriNo6NjkfWF695m8uTJGDdunPhcLpfD1dUVvr6+//oHWloUCgUiIyPh4+MDU1NTqeNoFGfVX4Y0ry7PWlBQgOvXr+P06dNYuXIlHj16BAAICAhA27Zt0a1bN5Wzx3R51vdhSPPq2qyFn/AUh84UouDgYFy/fh1nzpzR+PcyNzd/679yTE1Nte4XQBszaQpn1V+GNK+uzKpQKLBt2zbs3bsX586dU/mLxdraGlu2bEGPHj3e+R66Mqu6GNK8ujJrSTLqRCEaMWIEDh8+jFOnTqFy5cricicnJ+Tl5SEtLU1lL1FycjKcnJzEbS5evKjyfoVnoRVuQ0REbwiCgJMnT2Ly5Mm4cOGCuNza2hpeXl7w8/NDUFAQb51BekerC5EgCBg5ciT279+PkydPolq1airrmzRpAlNTU0RFRYn/Url9+zYSExPh5eUFAPDy8sL8+fPx4sULODg4AAAiIyNhY2ODunXrlu5ARERa6sWLF+JJKtevXwcA2NraYvLkyfDz84OnpyeMjY0lTkmkOVpdiIKDgxEWFoZffvkF1tbW4jE/tra2sLCwgK2tLQYNGoRx48bB3t4eNjY2GDlyJLy8vPDhhx8CAHx9fVG3bl3069cPS5YsQVJSEqZNm4bg4GAe/EdEhDdlqGHDhnj+/DmANzdWDQoKwoQJE4r8Q5RIX2l1IVq3bh0AoF27dirLt2zZgi+//BIAsGLFChgZGaFHjx7Izc2Fn58f1q5dK25rbGyMw4cP46uvvoKXlxesrKwQFBSEOXPmlNYYRERaKz09Hdu3bxfLEPDmv7FffPGFhKmISp9WF6LiXBGgTJkyWLNmDdasWfOP27i5ueG3335TZzQiIp1269YthIaGIiQkBBkZGeJyExMTXlCRDJJWFyIiIlKv7OxsTJkyBStXrhT/0enu7g4fHx+0b98ebdu25QknZJBYiIiIDERaWhoCAwPFPeYBAQHo168fevbsCSMjrb6TE5HGsRAREem5vLw8rFu3DnPmzEFKSgpMTU2xe/dudOvWTepoRFqDhYiISI8JgoBmzZrhjz/+AADUrVsX69evR5s2bSRORqRduI+UiEiPLVq0SCxDI0eOxNWrV1mGiN6Ce4iIiPRY4Z3oAWDlypWQyWQSpiHSXtxDRESkx5YsWSJ+vX//fgmTEGk3FiIiIj3211Pod+/eLWESIu3GQkREpKd+/fVXeHp6is9Hjx4tYRoi7cZCRESkZwRBwOLFi/Hxxx/jyZMnqFGjBk6cOCHe9JqIimIhIiLSI2FhYahZsyYmTZoEABg4cCBiY2OL3BOSiFTxLDMiIj2xcOFCTJkyBQBQtmxZDBo0CMuXL+dVqImKgYWIiEgP3Lx5EzNnzgQATJ06FVOmTIGlpaXEqYh0BwsREZEe+Pbbb6FQKNClSxfMnTuX1xsiKiHuRyUi0nFhYWH49ddfAQALFixgGSJ6DyxEREQ67Nq1a+jXrx8AoHfv3vDw8JA4EZFuYiEiItJRGzduhJeXF5RKJZo2bYqffvqJe4eI3hOPISIi0kF37tzB0KFDAQCNGjXC9u3beTYZ0X/AQkREpAP+/PNPREZG4ubNm7h165Z4X7JGjRohNjaWZYjoP2IhIiLSUunp6Th//jx++ukn7Ny5E/n5+Srr7e3tMXPmTJYhIjVgISIi0iJyuRxz5szBsWPH8Mcff0AQBHFd/fr18eGHH6J27dqoXbs22rVrBwsLCwnTEukPFiIiIi0yYsQI7NixQ3zu7u6ONm3aYMSIEWjatKmEyYj0GwsREZGWOH36NHbs2AGZTIZNmzYhICAAzs7OUsciMggsREREWmLt2rUAgI8//hgDBw6UOA2RYeGReEREWmLnzp0AgNzcXImTEBkeFiIiIi1hZ2cHAHB1dZU2CJEBYiEiItISPj4+AIDNmzejQYMGuHLlisSJiAwHCxERkZb4/vvvxa//+OMPdOzYEUqlUsJERIaDhYiISEJKpRIXLlzAhAkT4OnpqbKuZ8+evOgiUSnhWWZERKUsLy8PEREROHjwIA4dOoSkpCRxnZubG0aNGoW+ffvCwcFBwpREhoWFiIioFMXExODrr7/G1atXxWXW1tbw9/dHnz590LlzZ5iamkqYkMgwsRAREZWSLVu24JdffhGff/XVV+jWrRvatWsHMzMzCZMREQsREVEpSE9PF8tQnz59MHHiRNSvX1/iVERUiIWIiEjDBEHAkiVLAACVK1dGaGioxImI6O906vSFRYsWQSaTYcyYMeKynJwcBAcHo3z58ihbtix69OiB5ORkldclJiaic+fOsLS0hIODAyZMmID8/PxSTk9EhigzMxM9e/bE0qVLAQBjx46VOBERvY3OFKKYmBj88MMPRXYxjx07FocOHcKePXsQHR2NZ8+e4dNPPxXXFxQUoHPnzsjLy8PZs2exbds2bN26FTNmzCjtEYjIwMTHx6N9+/bYs2cPTExM8NVXX2HkyJFSxyKit9CJj8wyMjIQGBiIjRs3Yt68eeLy9PR0bN68GWFhYejQoQOANwct1qlTB+fPn8eHH36IiIgI3Lx5E8eOHYOjoyMaNmyIuXPnYuLEiZg1a9ZbD2TMzc1VuZeQXC4HACgUCigUCg1PWzyFObQljyZxVv2lr/M+ePAAY8eOxZEjRwAA9vb22LlzJ7KysvRu1rfR15/rPzGkeXVt1pLklAmCIGgwi1oEBQXB3t4eK1asQLt27dCwYUOEhITg+PHj6NixI1JTU8V7AAFvruMxZswYjB07FjNmzMDBgwcRFxcnrk9ISIC7uzsuX76MRo0aFfl+s2bNwuzZs4ssDwsLg6WlpSZGJCI98dtvv2Hr1q3Iy8uDkZERmjZtioEDB8LJyUnqaEQGJysrC3369EF6ejpsbGzeua3W7yHauXMnLl++jJiYmCLrkpKSYGZmplKGAMDR0VG80FlSUhIcHR2LrC9c9zaTJ0/GuHHjxOdyuRyurq7w9fX91z/Q0qJQKBAZGQkfHx+9v2YJZ9Vf+jbv+fPnsWHDBgBA27ZtsXbtWtSsWROA/s36LoY0K2BY8+rarIWf8BSHVheix48fY/To0YiMjESZMmVK7fuam5vD3Ny8yHJTU1Ot+wXQxkyawln1l77MO3/+fABA3759sX37dshksiLb6MusxWFIswKGNa+uzFqSjFp9UPWlS5fw4sULNG7cGCYmJjAxMUF0dDRWrVoFExMTODo6Ii8vD2lpaSqvS05OFndPOzk5FTnrrPA5d2ETkbpkZ2fj6NGjAICRI0e+tQwRkfbS6kLUsWNHXLt2DXFxceKjadOmCAwMFL82NTVFVFSU+Jrbt28jMTERXl5eAAAvLy9cu3YNL168ELeJjIyEjY0N6tatW+ozEZH+ycjIQFhYmPicd6gn0j1a/ZGZtbU16tWrp7LMysoK5cuXF5cPGjQI48aNg729PWxsbDBy5Eh4eXnhww8/BAD4+vqibt266NevH5YsWYKkpCRMmzYNwcHBb/1YjIioJKKjo/HJJ5+Ixyr07NkTLVq0kDgVEZWUVhei4lixYgWMjIzQo0cP5Obmws/PD2vXrhXXGxsb4/Dhw/jqq6/g5eUFKysrBAUFYc6cORKmJiJ9IJfLMXDgQMjlclSqVAmDBg3CtGnT+HEZkQ7SuUJ08uRJledlypTBmjVrsGbNmn98jZubG3777TcNJyMiQ/HixQvMmzcP27Ztg1wuh7OzM27cuAFbW1upoxHRe9K5QkREJKWsrCy0atUK9+7dAwBUrVoVu3btYhki0nEsRERExZSdnY3g4GDcu3cPdnZ22LVrF7y9vWFkpNXnpxBRMbAQEREVg1KpRMOGDXHnzh0AgL+/P3x9fSVORUTqwn/WEBEVQ35+vliGqlevjkWLFkmciIjUiXuIiIiKofCWHMCb44bCwsJQpUoV8eHi4qITV+4lordjISIiKobz58+LX0dFRalcEBYAjIyM0K1bN+zdu5en3RPpIH5kRkT0DwRBwLNnzxAXF4cBAwZgzJgx8PHxQaVKlYpsq1QqsW/fPixYsECCpET0X3EPERHR32zZsgUDBw58r9daWFioOQ0RlQYWIiKiv1m8eLHKcyMjI1SvXh3W1tYoW7asyv8Wfl21alV4e3ujSpUqEqUmov+ChYiIDN6LFy+wZs0aXLp0CVevXsWTJ09U1s+ZMwdTp06VKB0RlQYWIiIyOM+ePcP169dx48YNxMbGYvfu3cjPz1fZxtzcHOXLl0f9+vXRtWtXiZISUWlhISIig/Ho0SMMGjSoyBliAODo6IiJEyeiefPm+OCDD1ChQgWeLUZkQFiIiEjvZWdn4/79+xg8eDAuXLgAAKhTpw7q1q2LevXqoW3btmjZsiXMzc0lTkpEUmEhIiK9FR0djcDAQDx9+lRl+cSJE3mlaSJSwUJERHpr48aNYhkyMjJCkyZNEBgYiODgYImTEZG2YSEiIr0jl8tx8eJFlT1DN27cQO3atSVMRUTajIWIiPTG/PnzMW3atCLLP/30U9SqVUuCRESkK3jrDiLSCz/++ONbyxAA7N69m2eMEdE7cQ8REemUjIwMJCYm4vHjx0hMTMSdO3dw9uxZnD179q3bR0REwNjYuJRTEpGuYSEiIq2XnJyM06dPY8+ePdizZw8EQXjrdl988QW2bduGMmXKlHJCItJ171WIEhMT8ejRI2RlZaFixYrw8PDg9TuISK0UCgXmzZuH3bt349atWyrrbG1tUaVKFbi6uqJq1apo0aIFWrduDXd3d4nSEpGuK3YhevjwIdatW4edO3fiyZMnKv9CMzMzQ5s2bTB06FD06NEDRkY8NImI3l9ERATGjh2LmzdvAgBkMhk8PT3x0Ucf4csvv0STJk0kTkhE+qZYzWXUqFFo0KABEhISMG/ePNy8eRPp6enIy8tDUlISfvvtN7Ru3RozZsxA/fr1ERMTo+ncRKSHnj59ikqVKsHPzw83b96EnZ0dNm7ciFevXuHq1atYvXo1yxARaUSx9hBZWVnhwYMHKF++fJF1Dg4O6NChAzp06ICZM2ciPDwcjx8/RrNmzdQeloj026xZs/Ds2TPxeVpaGhYsWICffvoJrq6uaNasGUaOHMkzxohI7YpViBYuXCh+nZiYCAcHh388aNHf3189yYjI4Hz11Ve4efMmYmNjIZPJkJubi4SEBCQkJAAAfvrpJ7Rv3x6enp4SJyUifVOig6qVSiVq1KiBGzduoGbNmprKREQGIDExESdPnsSJEydw4sQJJCUlISkpCRkZGXB0dMSff/751teVK1eulJMSkSEoUSEyMjJCzZo18eeff7IQEdF7i4+PR9OmTZGVlVWs7a2srFCxYkWMHz8elStX1nA6IjJEJT7tftGiRZgwYQLWrVuHevXqaSITEem527dvFylDNWvWxLBhw+Dp6Yly5crBzs5OfJiamkqUlIgMRYnPj+/fvz8uXryIBg0awMLCAvb29ioPIqJ/07VrV6xevVrlZqt3797F+PHjoVAo0KxZM9SsWRMVK1ZkGSKiUlHiPUQhISEaiEFEhkQmk2HEiBEYNmwYNm7ciAcPHuC7774DAHz88cfo27cv1q1bh7Jly0qclIgMRYkLUVBQkCZyEJGBqlSpEvr27SsWIuDN2WQTJkxA/fr1JUxGRIbkP93LLCcnB3l5eSrLbGxs/lMgIjI89+7dU3n+8ccfIy8vD9nZ2bCwsJAoFREZkhIfQ5SZmYkRI0bAwcEBVlZWKFeunMqDiKikkpKSVJ4fPnwYzZo1g4ODA+7evStRKiIyJCUuRN9++y2OHz+OdevWwdzcHJs2bcLs2bPh4uKC7du3qz3g06dP0bdvX5QvXx4WFhbw9PREbGysuF4QBMyYMQPOzs6wsLCAt7d3kf+ApqSkIDAwEDY2NrCzs8OgQYOQkZGh9qxEVDJKpRIpKSl4/fr1W9dnZGTg/v37pZyKiAxRiT8yO3ToELZv34527dphwIABaNOmDWrUqAE3NzeEhoYiMDBQbeFSU1PRqlUrtG/fHkeOHEHFihVx9+5dlT1RS5YswapVq7Bt2zZUq1YN06dPF++DVHg17cDAQDx//hyRkZFQKBQYMGAAhg4dirCwMLVlJaJ3UyqVuHnzJqKjo3Hq1Clcu3YNCQkJyMnJUdmufv36qFSpEipVqgRPT0907NhRosREZEhKXIhSUlLg7u4O4M3xQikpKQCA1q1b46uvvlJruMWLF8PV1RVbtmwRl1WrVk38WhAEhISEYNq0aejatSsAYPv27XB0dMSBAwfQq1cvxMfHIzw8HDExMWjatCkAYPXq1ejUqROWLVsGFxcXtWYmoqKOHDmC/v3749WrV0XWGRkZoVq1aujRowfGjx+PihUrSpCQiAxdiQuRu7s7EhISUKVKFdSuXRu7d+9G8+bNcejQIdjZ2ak13MGDB+Hn54fPP/8c0dHRqFSpEr7++msMGTIEAJCQkICkpCR4e3uLr7G1tUWLFi1w7tw59OrVC+fOnYOdnZ1YhgDA29sbRkZGuHDhArp3717k++bm5iI3N1d8LpfLAQAKhQIKhUKtM76vwhzakkeTOKtuUSqVSE1NxYsXL/Dq1Su8ePECvXv3Ftd37NgRrVu3RvPmzeHm5oZbt24hICBAvN6QLs/+Lvrwsy0uQ5oVMKx5dW3WkuQscSEaMGAArl69irZt22LSpEno0qULvv/+eygUCixfvrykb/dODx48wLp16zBu3DhMmTIFMTExGDVqFMzMzBAUFCQeiOno6KjyOkdHR3FdUlISHBwcVNabmJjA3t6+yIGchRYuXIjZs2cXWR4REQFLS0t1jKY2kZGRUkcoNZxVuwiCgN9//x0XLlxAWloa0tPTIZfLIZfLoVQq3/oaW1tbjBw5EsCb/1Ddu3cPJiYmOjGvunBW/WVI8+rKrMW9PRDwHoVo7Nix4tfe3t6Ij4/H5cuXUaNGDbVfM0SpVKJp06ZYsGABAKBRo0a4fv061q9fr9HrIU2ePBnjxo0Tn8vlcri6usLX11drLiugUCgQGRkJHx8fvb+SL2fVTvPnz8eyZcv+cb2trS0qVqwIBwcHVKhQARYWFujduzc6deokbqNL8/5XnFV/GdK8ujZr4Sc8xfGfrkMEAFWrVkXVqlX/69u8lbOzM+rWrauyrE6dOvjf//4HAHBycgIAJCcnw9nZWdwmOTkZDRs2FLd58eKFynvk5+cjJSVFfP3fmZubw9zcvMhyU1NTrfsF0MZMmsJZtcvJkyffuT49PR0KhQJpaWl4+vQpypYti++//x779++Ho6Mjypcvj379+gHQjXnVhbPqL0OaV1dmLUnG9ypEUVFRWLFiBeLj4wG8KSljxoxROZZHHVq1aoXbt2+rLLtz5w7c3NwAvDnA2snJCVFRUWIBksvluHDhgniAt5eXF9LS0nDp0iU0adIEAHD8+HEolUq0aNFCrXmJDMmKFSvwzTff4Pnz58jMzERWVhYyMzNVzhrLyspCVlaWeDD1jRs3VN7j3Llz+PLLL0szNhHRW5X4OkRr166Fv78/rK2tMXr0aIwePRo2Njbo1KkT1qxZo9ZwY8eOxfnz57FgwQLcu3cPYWFh2LBhA4KDgwG8uR/SmDFjMG/ePBw8eBDXrl1D//794eLigm7dugF4U9b8/f0xZMgQXLx4Eb///jtGjBiBXr168Qwzov+gUaNGOH78OOLj45GYmIhXr14hOzsbBQUFeP36NZKTk/HgwQNcv34d0dHRKntxCw0cOFCC5ERERZV4D9GCBQuwYsUKjBgxQlw2atQotGrVCgsWLBDLijo0a9YM+/fvx+TJkzFnzhxUq1YNISEhKtc6+vbbb5GZmYmhQ4ciLS0NrVu3Rnh4uHgNIgAIDQ3FiBEj0LFjRxgZGaFHjx5YtWqV2nIS0f8xMjJC2bJlVW7MevDgQTx//lxlu/Lly2P27NlQKpXYsWMHrK2tYWVl9c5HpUqV4OHhAZlMVtpjEZGeK3EhSktLg7+/f5Hlvr6+mDhxolpC/dXHH3+Mjz/++B/Xy2QyzJkzB3PmzPnHbezt7XkRRiIJNW7cGDVr1lS5ivyff/6JP//8EwBw5cqVYr+Xi4sLpk2bpvbrnhGRYStxIfrkk0+wf/9+TJgwQWX5L7/88s7iQkSGq3Llyrhz5w4KCgqQmpoqlqHk5GScOHEC7u7uyMnJQWZm5lsfGRkZyMzMxB9//IFnz57h66+/hoeHBz766COpRyMiPVGsQvTXj5fq1q2L+fPn4+TJk/Dy8gIAnD9/Hr///ju++eYbzaQkIr1gbGyMChUqoEKFCgDenMJrZGSETp06FetskNOnT4slaNWqVSxERKQ2xSpEK1asUHlerlw53Lx5Ezdv3hSX2dnZ4ccff8S0adPUm5CI6P/767GBhZffICJSh2IVooSEBE3nICL6V40bN0b//v2xfft2WFtbSx2HiPRIiU+7JyKSirGxsfgxWe3atSVOQ0T6pFh7iP56G4t/o+77mRER/dXOnTsBAD4+PhInISJ9UqxC9PdTYi9fvoz8/HzUqlULwJurRxsbG4tXgiYi0oRff/0Vx44dg4mJCQYMGCB1HCLSI8UqRCdOnBC/Xr58OaytrbFt2zaUK1cOAJCamooBAwagTZs2mklJRAatoKAAjx49wty5cwEA7du3R40aNSRORUT6pMTXIfruu+8QEREhliHgzVln8+bNg6+vL0+9JyK1uHbtGrZv345jx47h1q1bKvdI495oIlK3EhciuVyOly9fFln+8uVLvH79Wi2hiMiwKRQKfPTRR0hLSxOXmZubo3bt2ujduzfGjh0rXTgi0kslLkTdu3fHgAED8N1336F58+YAgAsXLmDChAn49NNP1R6QiAyPsbExBEEA8Obg6XXr1qFq1aowNjaWOBkR6asSn3a/fv16BAQEoE+fPnBzc4Obmxv69OkDf39/rF27VhMZicjAGBkZ4YsvvgDw5vT66tWrswwRkUaVeA+RpaUl1q5di6VLl+L+/fsAgOrVq8PKykrt4YjIsCgUCly9ehWbN2/Gli1bAADu7u4SpyIiQ1DiQlTIysoK9evXV2cWIjJQL168wIQJE7B//36Vg6f9/PwwdOhQCZMRkaEo1kdmw4cPx5MnT4r1hrt27UJoaOh/CkVEhqGgoADR0dFo1aoVfv75Z+Tk5KBcuXL45JNPcPLkSYSHh8PS0lLqmERkAIq1h6hixYrw8PBAq1at0KVLFzRt2hQuLi4oU6YMUlNTcfPmTZw5cwY7d+6Ei4sLNmzYoOncRKTjUlJS0LlzZ5w/fx4AULlyZezZswctWrSATCaTOB0RGZpiFaK5c+dixIgR2LRpE9auXatyl3sAsLa2hre3NzZs2AB/f3+NBCUi/bJ9+3axDM2aNQtjx46FjY2NxKmIyFAV+xgiR0dHTJ06FVOnTkVqaioSExORnZ2NChUqoHr16vwXHREVW0FBATZv3gwA6NatG6ZMmQJTU1OJUxGRIXuvg6rLlSuncqVqIqLiOn36NL777jtcv34dtra2+Oyzz6SORET0/meZEREVV3Z2No4fP46ffvpJvFs9AMyePRtly5aVMBkR0RssRESkEZmZmYiMjMTvv/+On376CUlJSQAAmUyGQYMGYdiwYWjQoAF+++03iZMSEbEQEZGaKZVKbNiwAbNnzxZLEPDmLLKPP/4YgwcPFm/OqlAopIpJRKSChYiI1GrkyJHibXyqVKkCf39/tG7dGl988QXMzc0lTkdE9HbvVYjy8/Nx8uRJ3L9/H3369IG1tTWePXsGGxsbHg9AZMDu3bsnlqEhQ4bg+++/h5mZmcSpiIj+XYlv7vro0SN4enqia9euCA4OxsuXLwEAixcvxvjx49UekIh0x18/Itu4caPKAdRERNqsxIVo9OjRaNq0KVJTU2FhYSEu7969O6KiotQajoh0x9WrV9G3b1+VZbzQIhHpihJ/ZHb69GmcPXu2yG7wqlWr4unTp2oLRkTa7/Xr1zh69Ch++eUX7Nq1CwqFAo6OjhgwYAD69++POnXqSB2RiKhYSlyIlEolCgoKiix/8uQJrK2t1RKKiLSXQqFAeHg4tm/fjkOHDiE3N1dc17VrV2zYsAEODg4SJiQiKrkSFyJfX1+EhISIN3CVyWTIyMjAzJkz0alTJ7UHJCLtkJeXh++//x7fffcdnj17Ji6vXr06unXrhm7duqFVq1a8jQ8R6aQSF6Jly5bB398fdevWRU5ODvr06YO7d++iQoUK+PnnnzWRkYgklJWVhd27d2PJkiWIj48HADg4OCAwMBD9+vVDw4YNWYKISOeVuBC5urri6tWr2LVrF65evYqMjAwMGjQIgYGBKgdZE5Fuy8nJweLFi7FixQqkp6cDACpWrIgFCxagf//+PJ2eiPRKiQqRQqFA7dq1cfjwYQQGBiIwMFBTuYhIAkqlEnfu3MGFCxewcOFC3L59GwDg7u6OIUOGYNiwYbyxMxHppRIVIlNTU+Tk5GgqCxFJQBAEHD9+HFu2bEF4eDj+/PNPcZ2TkxNWrlyJzz77DEZGJb5KBxGRzijxR2bBwcFYvHgxNm3aBBMT3vmDSFclJSXh5MmT2LhxI44fPy4ut7CwQOPGjdG2bVtMmDABdnZ20oUkIiolJW40MTExiIqKQkREBDw9PWFlZaWyft++fWoLV1BQgFmzZol3ynZxccGXX36JadOmiQdxCoKAmTNnYuPGjUhLS0OrVq2wbt061KxZU3yflJQUjBw5EocOHYKRkRF69OiBlStX8jYjZDCUSiXOnz+P6OhoxMTEICYmBk+ePBHXm5ubY/DgwejZsyc+/PBDmJqaSpiWiKj0lbgQ2dnZoUePHprIUsTixYuxbt06bNu2DR4eHoiNjcWAAQNga2uLUaNGAQCWLFmCVatWYdu2bahWrRqmT58OPz8/3Lx5E2XKlAEABAYG4vnz54iMjIRCocCAAQMwdOhQhIWFlcocRFJ68uQJfH19xTPECslkMnh4eKBTp04YNGgQPvjgA4kSEhFJr8SFaMuWLZrI8VZnz55F165d0blzZwBvrob9888/4+LFiwDe7B0KCQnBtGnT0LVrVwDA9u3b4ejoiAMHDqBXr16Ij49HeHg4YmJi0LRpUwDA6tWr0alTJyxbtgwuLi6lNg+RFObPn4/4+HhYW1vDz88PLVq0QLNmzdC4cWNeTJWI6P/T6oOAWrZsiQ0bNuDOnTv44IMPcPXqVZw5cwbLly8HACQkJCApKQne3t7ia2xtbdGiRQucO3cOvXr1wrlz52BnZyeWIQDw9vaGkZERLly4gO7duxf5vrm5uSpX35XL5QDenGWnUCg0NW6JFObQljyaxFnf3/nz57Ft2zYAbz7Obtu27Vu/n1T4s9VPhjQrYFjz6tqsJclZ4kJUrVq1d16E7cGDByV9y380adIkyOVy1K5dG8bGxigoKMD8+fPF0/0L76zt6Oio8jpHR0dxXVJSUpHbCJiYmMDe3l7lztx/tXDhQsyePbvI8oiICFhaWv7nudQpMjJS6gilhrMWj0KhwNmzZ3H48GHcvXsXAFCjRg1kZGTgt99+U1dEteLPVj8Z0qyAYc2rK7NmZWUVe9sSF6IxY8aoPFcoFLhy5QrCw8MxYcKEkr7dO+3evRuhoaEICwuDh4cH4uLiMGbMGLi4uCAoKEit3+uvJk+ejHHjxonP5XI5XF1d4evrqzV371YoFIiMjISPj4/eHwDLWYvvf//7H7755hvx1hpmZmbo2bMnZs2aBVdXV3XH/c/4s9VPhjQrYFjz6tqshZ/wFEeJC9Ho0aPfunzNmjWIjY0t6du904QJEzBp0iT06tULAODp6YlHjx5h4cKFCAoKgpOTEwAgOTkZzs7O4uuSk5PRsGFDAG+uo/LixQuV983Pz0dKSor4+r8zNzeHubl5keWmpqZa9wugjZk0hbO+nSAIuHHjBkJCQrB582YAgLOzM7766isMGzZMJ260yp+tfjKkWQHDmldXZi1JRrVdaS0gIAD/+9//1PV2AN7s6vr7xeCMjY2hVCoBvPn4zsnJCVFRUeJ6uVyOCxcuwMvLCwDg5eWFtLQ0XLp0Sdzm+PHjUCqVaNGihVrzEpW2mzdvok2bNvD09MTmzZshk8kwZcoUPHz4ENOnT9eJMkREpA3UdlD13r17YW9vr663AwB06dIF8+fPR5UqVeDh4YErV65g+fLlGDhwIIA3pw2PGTMG8+bNQ82aNcXT7l1cXNCtWzcAQJ06deDv748hQ4Zg/fr1UCgUGDFiBHr16sUzzEinxcTEwN/fHykpKTAzM4OPjw+++eYbtG/fXupoREQ6p8SFqFGjRioHVQuCgKSkJLx8+RJr165Va7jVq1dj+vTp+Prrr/HixQu4uLhg2LBhmDFjhrjNt99+i8zMTAwdOhRpaWlo3bo1wsPDxWsQAUBoaChGjBiBjh07ihdmXLVqlVqzEpWWq1evYtasWThw4AAAoFmzZti/fz8qVaokbTAiIh1W4kJUuOelkJGRESpWrIh27dqhdu3a6soFALC2tkZISAhCQkL+cRuZTIY5c+Zgzpw5/7iNvb09L8JIOu/58+dYtmwZVq9eLZ5K2rJlS+zbt6/ImZZERFQyJS5EM2fO1EQOInqH6OhoBAQEIDs7GwDQtWtXzJ8/Hx4eHhInIyLSDyU+qPry5cu4du2a+PyXX35Bt27dMGXKFOTl5ak1HJGhu3z5Mr766iu0a9dOLEO//fYb9u/fzzJERKRGJS5Ew4YNw507dwC8uQhjz549YWlpiT179uDbb79Ve0AiQ3Tz5k189tlnaNKkCdavXw8A8PDwwL179xAQEPDOi6MSEVHJlbgQ3blzR7zGz549e9C2bVuEhYVh69ataj/tnsjQpKenY/Xq1WjUqBH+97//wcjICL169UJUVBT++OMPVK9eXeqIRER6qcTHEAmCIF4H6NixY/j4448BAK6urnj16pV60xEZkH379mH06NF48uQJAODTTz/F7NmzUa9ePYmTERHpvxIXoqZNm2LevHnw9vZGdHQ01q1bB+DNjVZ5pgvR+4mNjcXnn38OpVIJJycn7Ny5s8iNWImISHNKXIhCQkIQGBiIAwcOYOrUqahRowaANxdmbNmypdoDEum71NRU9O/fH0qlEv7+/hgwYAD/v0REVMpKXIjq16+vcpZZoaVLl8LY2FgtoYgMwb1797Bjxw5s27YNjx49gpWVFRYtWoSHDx9KHY2IyOCUuBA9fvwYMpkMlStXBgBcvHgRYWFhqFu3LoYOHar2gET6JiUlBWvWrMHChQvFU+mdnZ1x5MgR1K1bl4WIiEgCJT7LrE+fPjhx4gQAICkpCT4+Prh48SKmTp36zqtFExm61NRUfPXVV6hcuTJmzJiB7OxstG7dGps2bcK1a9fQoEEDqSMSERmsEhei69evo3nz5gCA3bt3o169ejh79ixCQ0OxdetWdecj0gtnz55Fo0aNsH79emRnZ6Nhw4b4+eefER0djUGDBqF8+fJSRyQiMmgl/shMoVDA3NwcwJvT7j/55BMAQO3atfH8+XP1piPSA8nJyWjTpg2USiWqVq2KzZs3o3379ry4IhGRFinxHiIPDw+sX78ep0+fRmRkJPz9/QEAz549479yif4mNzcXXbp0gVKpRLly5XDt2jV06NCBZYiISMuUeA/R4sWL0b17dyxduhRBQUHicQ8HDx4UP0ojMmQZGRm4ceMG4uPjsWXLFsTExMDOzg7Hjx9H2bJlpY5HRERvUeJC1K5dO7x69QpyuRzlypUTlw8dOhSWlpZqDUekS3JycjBq1Chs3LhRZbm1tTX2798v3vKGiIi0T4kLEfDm9h2XLl3C/fv30adPH1hbW8PMzIyFiAzawYMHVcpQhw4dULduXQwfPpx3pici0nIlLkSPHj2Cv78/EhMTkZubCx8fH1hbW2Px4sXIzc0V78xNZCiio6MxY8YMnDp1Slw2c+ZMzJo1S7pQRERUIiU+qHr06NFo2rQpUlNTYWFhIS7v3r07oqKi1BqOSJvl5eVh8ODBaNeuHU6dOgWZTIYOHTpgw4YNmDx5stTxiIioBEq8h+j06dM4e/YszMzMVJZXrVoVT58+VVswIm2Wn5+PgIAAHD9+HEZGRhg2bBimTp2KSpUqSR2NiIjeQ4kLkVKpREFBQZHlT548gbW1tVpCEWmT1NRU3L9/H/fu3cP9+/dx//59nDp1Cvfv34eVlRXCwsLE63EREZFuKnEh8vX1RUhICDZs2AAAkMlkyMjIwMyZM9GpUye1ByQqTenp6di9ezcuXbqEuLg43LlzB6mpqW/d1srKCj/99BPLEBGRHihxIVq2bBn8/f1Rt25d5OTkoE+fPrh79y4qVKiAn3/+WRMZiUpFbm4uWrRogdu3bxdZ5+TkhOrVq4uPunXrIiAgAFZWVhIkJSIidStxIXJ1dcXVq1exa9cuXL16FRkZGRg0aBACAwNVDrIm0iWCIGDOnDliGfr222/RuHFj1KlTB+7u7rygIhGRnitRIVIoFKhduzYOHz6MwMBABAYGaioXUalZu3YtvvnmG+Tk5AAAvv/+ewQHB0ucioiISlOJTrs3NTUV/9Ig0hdz584Vf6/Hjx+Pr7/+WuJERERU2kp8HaLg4GAsXrwY+fn5mshDVKrkcrnKqfKffPIJb7xKRGSASnwMUUxMDKKiohAREQFPT88iB5Xu27dPbeGING3RokW4dOkSAKBNmzbw8vKSOBEREUmhxIXIzs4OPXr00EQWolJnbGwMAOjSpQv27t0LE5P3ur0fERHpuBL/13/Lli2ayEEkiYoVKwJ4cz2tv199nYiIDEexjyFSKpVYvHgxWrVqhWbNmmHSpEnIzs7WZDYijYqJiRHvTs+rrBMRGbZiF6L58+djypQpKFu2LCpVqoSVK1fy1GTSOdnZ2di/fz8CAgLQvHlzXL9+HeXKlcPs2bOljkZERBIqdiHavn071q5di6NHj+LAgQM4dOgQQkNDoVQqNZmPSC2ePXuGfv36wcHBAZ9++inCw8NhZGSE/v374/Lly6hevbrUEYmISELFPoYoMTFR5V5l3t7ekMlkePbsGSpXrqyRcETqEB4ejmHDhiExMRHAm6ut9+rVC0OHDkWNGjUkTkdERNqg2IUoPz8fZcqUUVlmamoKhUKh9lBE6vD48WOMGTNGvBREpUqV8OOPP8LHx4fXGiIiIhXF/shMEAR8+eWX+PTTT8VHTk4Ohg8frrKsJE6dOoUuXbrAxcUFMpkMBw4cKPI9Z8yYAWdnZ1hYWMDb2xt3795V2SYlJQWBgYGwsbGBnZ0dBg0ahIyMDJVt/vjjD7Rp0wZlypSBq6srlixZUqKcpFvy8vKwePFi1K5dG/v27YOxsTHGjRuH+Ph4+Pr6sgwREVERxS5EQUFBcHBwgK2trfjo27cvXFxcVJaVRGZmJho0aIA1a9a8df2SJUuwatUqrF+/HhcuXICVlRX8/PxUbh8SGBiIGzduIDIyEocPH8apU6cwdOhQcb1cLoevry/c3Nxw6dIlLF26FLNmzcKGDRtKlJV0w8mTJ9GwYUNMmjQJWVlZaN26Na5cuYLvvvuOZ5IREdE/KvZHZpq4/lBAQAACAgLeuk4QBISEhGDatGno2rUrgDcHdjs6OuLAgQPo1asX4uPjER4ejpiYGDRt2hQAsHr1anTq1AnLli2Di4sLQkNDkZeXhx9//BFmZmbw8PBAXFwcli9frlKcSLc9f/4c48ePR1hYGIA31xdatmwZ+vXrxz1CRET0r7T2srwJCQlISkqCt7e3uMzW1hYtWrTAuXPn0KtXL5w7dw52dnZiGQLeHOxtZGSECxcuoHv37jh37hw++ugjlYvu+fn5YfHixUhNTUW5cuWKfO/c3Fzk5uaKz+VyOQBAoVBozTFThTm0JY8mvWvWnJwcrFq1CosWLUJGRgZkMhmGDRuG2bNno1y5cjp3zz1D+rkChjUvZ9VfhjSvrs1akpxaW4iSkpIAAI6OjirLHR0dxXVJSUlwcHBQWW9iYgJ7e3uVbapVq1bkPQrXva0QLVy48K3XpYmIiIClpeV7TqQZkZGRUkcoNX+dVRAEnDt3Dlu3bsWLFy8AADVr1sSwYcNQo0YNnDt3TqqYamFIP1fAsOblrPrLkObVlVmzsrKKva3WFiIpTZ48GePGjROfy+VyuLq6wtfXFzY2NhIm+z8KhQKRkZHw8fGBqamp1HE06u+zXr58GePHj8eZM2cAvDl7bP78+ejVqxeMjIp9WJxWMqSfK2BY83JW/WVI8+rarIWf8BSH1hYiJycnAEBycjKcnZ3F5cnJyWjYsKG4TeHegUL5+flISUkRX+/k5ITk5GSVbQqfF27zd+bm5jA3Ny+y3NTUVOt+AbQxk6a8fPkSs2bNwrZt2yAIAiwsLDBx4kSMHz8eVlZWUsdTK0P6uQKGNS9n1V+GNK+uzFqSjFr7z+lq1arByckJUVFR4jK5XI4LFy7Ay8sLAODl5YW0tDRcunRJ3Ob48eNQKpVo0aKFuM2pU6dUPkeMjIxErVq13vpxGWmf7Oxs7NmzBx4eHti6dSsEQUDfvn1x584dzJw5U+/KEBERlT5JC1FGRgbi4uIQFxcH4M2B1HFxcUhMTIRMJsOYMWMwb948HDx4ENeuXUP//v3h4uKCbt26AQDq1KkDf39/DBkyBBcvXsTvv/+OESNGoFevXnBxcQEA9OnTB2ZmZhg0aBBu3LiBXbt2YeXKlSofiZF2EgQBu3btgqenJ0JDQ5GZmYkPP/wQ58+fx44dO3iFdCIiUhtJPzKLjY1F+/btxeeFJSUoKAhbt27Ft99+i8zMTAwdOhRpaWlo3bo1wsPDVa6YHRoaihEjRqBjx44wMjJCjx49sGrVKnG9ra0tIiIiEBwcjCZNmqBChQqYMWMGT7nXcjExMRgzZgzOnj0LAKhQoQKWL1+Ovn378jR6IiJSO0kLUbt27SAIwj+ul8lkmDNnDubMmfOP29jb24vXnvkn9evXx+nTp987J5Wep0+fYsqUKdi+fTsAwNLSEhMmTEDdunXRvXt3liEiItIIrT2GiAxLVlYW5s6diw8++EAsQ/3798edO3cwderUtx7kTkREpC5ae5YZGQZBELBz505MnDgRjx8/BgC0bNkSISEhaNasGQDduQAYERHpLhYiksyFCxcwZswYnD9/HgBQpUoVLFmyBF988QU/GiMiolLFj8xIEosWLRLPGLOyssK8efNw69Yt9OzZk2WIiIhKHfcQUak7cuQIJk+eDODNGYULFiwQL5NAREQkBRYiKlUPHz5E3759AQBfffUV1q5dK3EiIiIifmRGpSg3Nxeff/45UlJS0KxZM6xYsULqSERERABYiKgUjRkzBrGxsbC3t8eePXt4Kj0REWkNFiIqFTt27MD69eshk8kQGhoKNzc3qSMRERGJWIhI465du4Zhw4YBAGbMmAF/f3+JExEREaliISKNSk9PR48ePZCdnQ0/Pz9Mnz5d6khERERFsBCRxgiCgIEDB+Lu3buoUqUKfvrpJxgbG0sdi4iIqAgWItKY5cuXY9++fTA1NcWePXtQoUIFqSMRERG9FQsRacTp06cxceJEAMDKlSvRvHlziRMRERH9MxYiUrvExET07NkTBQUFCAwMxPDhw6WORERE9E4sRKRWT58+RYcOHfD8+XN4eHjghx9+4L3JiIhI67EQkdo8f/4c7du3x/3791GtWjUcOXIEVlZWUsciIiL6VyxEpBbJycno0KED7t69Czc3N5w4cQKurq5SxyIiIioWFiL6z16+fImOHTvi1q1bqFy5Mo4fP84rURMRkU5hIaL/JCUlBT4+Prhx4wacnZ1x4sQJuLu7Sx2LiIioRFiI6L2lpqbCx8cHV69ehaOjI06cOIEaNWpIHYuIiKjEWIjovaSnp8PPzw+XL19GxYoVcfz4cdSqVUvqWERERO+FhYhK7PXr1wgICEBMTAzKly+PqKgo1K1bV+pYRERE742FiEokIyMDnTp1wrlz51CuXDkcO3YMnp6eUsciIiL6T1iIqNiysrLQpUsXnDlzBra2toiMjETDhg2ljkVERPSfsRBRseTm5qJr1644efIkrK2tcfToUTRp0kTqWERERGrBQkTFsmDBAhw7dgxWVlY4cuQIWrRoIXUkIiIitWEhon91+/ZtLFq0CADw448/olWrVhInIiIiUi8WInonQRAwfPhw5OXlISAgAJ9//rnUkYiIiNSOhYjeaceOHTh58iQsLCywZs0a3rmeiIj0EgsR/aM///wT33zzDQBgxowZqFatmsSJiIiINIOFiP7RxIkT8erVK3h4eIjFiIiISB+xENFbnT59Gps3bwYA/PDDDzA1NZU4ERERkeawEFEReXl5GD58OABg8ODBPKuMiIj0nqSF6NSpU+jSpQtcXFwgk8lw4MABcZ1CocDEiRPh6ekJKysruLi4oH///nj27JnKe6SkpCAwMBA2Njaws7PDoEGDkJGRobLNH3/8gTZt2qBMmTJwdXXFkiVLSmM8nfXdd9/h5s2bqFixIhYvXix1HCIiIo2TtBBlZmaiQYMGWLNmTZF1WVlZuHz5MqZPn47Lly9j3759uH37Nj755BOV7QIDA3Hjxg1ERkbi8OHDOHXqFIYOHSqul8vl8PX1hZubGy5duoSlS5di1qxZ2LBhg8bn00UPHjzAnDlzALwpRvb29hInIiIi0jwTKb95QEAAAgIC3rqu8F5Zf/X999+jefPmSExMRJUqVRAfH4/w8HDExMSgadOmAIDVq1ejU6dOWLZsGVxcXBAaGoq8vDz8+OOPMDMzg4eHB+Li4rB8+XKV4kRvrjkUHByMnJwcdOjQAX379pU6EhERUamQtBCVVHp6OmQyGezs7AAA586dg52dnViGAMDb2xtGRka4cOECunfvjnPnzuGjjz6CmZmZuI2fnx8WL16M1NRUlCtXrsj3yc3NRW5urvhcLpcDePMxnkKh0NB0JVOYQ5155s2bh/DwcJiZmWHlypXIz89X23v/F5qYVVsZ0qyAYc3LWfWXIc2ra7OWJKfOFKKcnBxMnDgRvXv3ho2NDQAgKSkJDg4OKtuZmJjA3t4eSUlJ4jZ/v36Oo6OjuO5thWjhwoWYPXt2keURERGwtLRUyzzq8ve9aO/r0KFD4lllQUFBuH//Pu7fv6+W91YXdc2qCwxpVsCw5uWs+suQ5tWVWbOysoq9rU4UIoVCgS+++AKCIGDdunUa/36TJ0/GuHHjxOdyuRyurq7w9fUVy5jUFAoFIiMj4ePj859Pid+6datYhmbMmIFp06apI6LaqHNWbWdIswKGNS9n1V+GNK+uzVr4CU9xaH0hKixDjx49wvHjx1UKiZOTE168eKGyfX5+PlJSUuDk5CRuk5ycrLJN4fPCbf7O3Nwc5ubmRZabmppq3S/Af820Z88e8RT7cePGYdasWVp7ew5t/PPXFEOaFTCseTmr/jKkeXVl1pJk1OrrEBWWobt37+LYsWMoX768ynovLy+kpaXh0qVL4rLjx49DqVSiRYsW4janTp1S+RwxMjIStWrVeuvHZYYkPDwcgYGBUCqVGDx4MJYtW6a1ZYiIiEiTJC1EGRkZiIuLQ1xcHAAgISEBcXFxSExMhEKhwGeffYbY2FiEhoaioKAASUlJSEpKQl5eHgCgTp068Pf3x5AhQ3Dx4kX8/vvvGDFiBHr16gUXFxcAQJ8+fWBmZoZBgwbhxo0b2LVrF1auXKnykZghOn36ND799FMoFAr07NkT69evZxkiIiKDJelHZrGxsWjfvr34vLCkBAUFYdasWTh48CAAoGHDhiqvO3HiBNq1awcACA0NxYgRI9CxY0cYGRmhR48eWLVqlbitra0tIiIiEBwcjCZNmqBChQqYMWOGQZ9yf/nyZXz88cfIzs5Gp06dsH37dhgbG0sdi4iISDKSFqJ27dpBEIR/XP+udYXs7e0RFhb2zm3q16+P06dPlzifPoqPj4efnx/kcjnatm2LvXv3qlySgIiIyBBp9TFEpF4JCQnw9vbGq1ev0KxZMxw8eBAWFhZSxyIiIpIcC5GBePbsGby9vfHs2TN4eHjgyJEjWnMJASIiIqmxEBmAP//8E76+vnjw4AHc3d0RGRlZ5Iw9IiIiQ8ZCpOcEQUDPnj1x48YNVKpUCceOHYOzs7PUsYiIiLQKC5GeO3LkCKKiomBubo6IiIgitzEhIiIiFiK9VlBQgEmTJgEARo0ahbp160qciIiISDuxEOmxsLAwXLt2DXZ2dmIxIiIioqJYiPRUTk6OeJPWSZMmwd7eXuJERERE2ouFSE+tW7cOiYmJqFSpEkaNGiV1HCIiIq3GQqSH0tPTMW/ePADA7NmzefFFIiKif8FCpIeWLFmClJQU1KlTB0FBQVLHISIi0nosRHrm+fPnWLFiBQBg4cKFMDGR9HZ1REREOoGFSM/Mnj0b2dnZaNmyJT755BOp4xAREekEFiI9cvv2bWzatAkAsHjxYshkMokTERER6QYWIj0ybdo0FBQUoEuXLmjdurXUcYiIiHQGC5GeuHXrFv73v/8BABYsWCBxGiIiIt3CQqQnli5dCkEQ0LVrV9SrV0/qOERERDqFhUgPPHv2DDt27AAATJw4UeI0REREuoeFSA+EhIRAoVCgTZs28PLykjoOERGRzmEh0nFpaWlYv349AO4dIiIiel8sRDpuw4YNeP36NTw8PBAQECB1HCIiIp3EQqTD8vLysHr1agDAt99+CyMj/jiJiIjeB/8G1WEnT55EcnIyXF1d0bt3b6njEBER6SwWIh1VUFCAAwcOAADGjRsHU1NTaQMRERHpMBYiHfXLL7/g2bNnKFeuHAYPHix1HCIiIp3GQqSDBEHAsmXLAADDhw9H2bJlJU5ERESk21iIdNCZM2cQGxsLMzMzBAcHSx2HiIhI57EQ6aAVK1YAANq2bQsHBweJ0xAREek+FiId8+DBA/Fg6k8++UTaMERERHqChUjHrFq1CoIgwNfXF66urlLHISIi0gssRDokPT0dmzdvBgCMGjVK4jRERET6g4VIh2zevBkZGRmoW7cufHx8pI5DRESkN1iIdER+fj5WrVoFABgzZgxkMpnEiYiIiPQHC5GOOHDgAB49eoQKFSqgb9++UschIiLSK5IWolOnTqFLly5wcXGBTCYTz556m+HDh0MmkyEkJERleUpKCgIDA2FjYwM7OzsMGjQIGRkZKtv88ccfaNOmDcqUKQNXV1csWbJEA9NoVuGp9sOHD4eFhYXEaYiIiPSLpIUoMzMTDRo0wJo1a9653f79+3H+/Hm4uLgUWRcYGIgbN24gMjIShw8fxqlTpzB06FBxvVwuh6+vL9zc3HDp0iUsXboUs2bNwoYNG9Q+j6ZcvHgRZ8+ehampKb7++mup4xAREekdEym/eUBAAAICAt65zdOnTzFy5EgcPXoUnTt3VlkXHx+P8PBwxMTEoGnTpgCA1atXo1OnTli2bBlcXFwQGhqKvLw8/PjjjzAzM4OHhwfi4uKwfPlyleKkzQr3DvXu3RvOzs4SpyEiItI/khaif6NUKtGvXz9MmDABHh4eRdafO3cOdnZ2YhkCAG9vbxgZGeHChQvo3r07zp07h48++ghmZmbiNn5+fli8eDFSU1NRrly5Iu+bm5uL3Nxc8blcLgcAKBQKKBQKdY74r169eoW9e/cCAIKDg8Xv//f/1WecVX8Z0rycVX8Z0ry6NmtJcmp1IVq8eDFMTEz+8Zo7SUlJRW5dYWJiAnt7eyQlJYnbVKtWTWUbR0dHcd3bCtHChQsxe/bsIssjIiJgaWn5XrO8r6NHjyI/Px/u7u54/vw5nj9/rrI+MjKyVPNIibPqL0Oal7PqL0OaV1dmzcrKKva2WluILl26hJUrV+Ly5culfor55MmTMW7cOPG5XC6Hq6srfH19YWNjU6pZCg8iHzx4MDp16iQuVygUiIyMhI+PD0xNTUs1U2njrPrLkOblrPrLkObVtVkLP+EpDq0tRKdPn8aLFy9QpUoVcVlBQQG++eYbhISE4OHDh3BycsKLFy9UXpefn4+UlBQ4OTkBAJycnJCcnKyyTeHzwm3+ztzcHObm5kWWm5qaluovQFJSEk6dOgXgzfFDb/vepZ1JSpxVfxnSvJxVfxnSvLoya0kyau11iPr164c//vgDcXFx4sPFxQUTJkzA0aNHAQBeXl5IS0vDpUuXxNcdP34cSqUSLVq0ELc5deqUyueIkZGRqFWr1ls/LtMme/fuFWepWrWq1HGIiIj0lqR7iDIyMnDv3j3xeUJCAuLi4mBvb48qVaqgfPnyKtubmprCyckJtWrVAgDUqVMH/v7+GDJkCNavXw+FQoERI0agV69e4in6ffr0wezZszFo0CBMnDgR169fx8qVK8Uzt7TZzp07AQC9evWSOAkREZF+k7QQxcbGon379uLzwuN2goKCsHXr1mK9R2hoKEaMGIGOHTvCyMgIPXr0EG9xAQC2traIiIhAcHAwmjRpggoVKmDGjBlaf8r948eP8fvvv0Mmk+Hzzz+XOg4REZFek7QQtWvXDoIgFHv7hw8fFllmb2+PsLCwd76ufv36OH36dEnjSWr37t0AgDZt2qBSpUoSpyEiItJvWnsMkaHbtWsXAKBnz54SJyEiItJ/LERa6MGDB4iJiRE/AiQiIiLNYiHSQoV7hzp06CBeRJKIiIg0h4VIyyiVSmzfvh0APy4jIiIqLSxEWubQoUO4desWbGxseHYZERFRKWEh0iKCIGDhwoUA3tzI1dbWVuJEREREhoGFSItER0fjwoULKFOmDEaPHi11HCIiIoPBQqRFFi1aBAAYOHAgD6YmIiIqRSxEWuLKlSs4evQojI2NMX78eKnjEBERGRQWIi1RuHeoZ8+eqFatmsRpiIiIDAsLkRa4e/cu9u7dCwCYOHGixGmIiIgMDwuRFli2bBmUSiU6deqE+vXrSx2HiIjI4LAQSez58+fYunUrAGDy5MnShiEiIjJQLEQSW7FiBfLy8tCqVSu0bt1a6jhEREQGiYVIQqmpqVi3bh0AYNKkSRKnISIiMlwmUgcwZA8fPoSTkxPKlCmDTp06SR2HiIjIYLEQSahRo0a4desWnj59CiMj7qwjIiKSCv8WlpixsTGqVKkidQwiIiKDxkJEREREBo+FiIiIiAweCxEREREZPBYiIiIiMngsRERERGTwWIiIiIjI4LEQERERkcFjISIiIiKDx0JEREREBo+FiIiIiAweCxEREREZPBYiIiIiMngsRERERGTwTKQOoAsEQQAAyOVyiZP8H4VCgaysLMjlcpiamkodR6M4q/4ypHk5q/4ypHl1bdbCv7cL/x5/FxaiYnj9+jUAwNXVVeIkREREVFKvX7+Gra3tO7eRCcWpTQZOqVTi2bNnsLa2hkwmkzoOgDet19XVFY8fP4aNjY3UcTSKs+ovQ5qXs+ovQ5pX12YVBAGvX7+Gi4sLjIzefZQQ9xAVg5GRESpXrix1jLeysbHRiV9KdeCs+suQ5uWs+suQ5tWlWf9tz1AhHlRNREREBo+FiIiIiAweC5GOMjc3x8yZM2Fubi51FI3jrPrLkOblrPrLkObV51l5UDUREREZPO4hIiIiIoPHQkREREQGj4WIiIiIDB4LERERERk8FiIdsm7dOtSvX1+8IJaXlxeOHDkidaxSsWjRIshkMowZM0bqKBoxa9YsyGQylUft2rWljqUxT58+Rd++fVG+fHlYWFjA09MTsbGxUsfSiKpVqxb52cpkMgQHB0sdTe0KCgowffp0VKtWDRYWFqhevTrmzp1brPtI6aLXr19jzJgxcHNzg4WFBVq2bImYmBipY6nFqVOn0KVLF7i4uEAmk+HAgQMq6wVBwIwZM+Ds7AwLCwt4e3vj7t270oRVExYiHVK5cmUsWrQIly5dQmxsLDp06ICuXbvixo0bUkfTqJiYGPzwww+oX7++1FE0ysPDA8+fPxcfZ86ckTqSRqSmpqJVq1YwNTXFkSNHcPPmTXz33XcoV66c1NE0IiYmRuXnGhkZCQD4/PPPJU6mfosXL8a6devw/fffIz4+HosXL8aSJUuwevVqqaNpxODBgxEZGYkdO3bg2rVr8PX1hbe3N54+fSp1tP8sMzMTDRo0wJo1a966fsmSJVi1ahXWr1+PCxcuwMrKCn5+fsjJySnlpGokkE4rV66csGnTJqljaMzr16+FmjVrCpGRkULbtm2F0aNHSx1JI2bOnCk0aNBA6hilYuLEiULr1q2ljiGZ0aNHC9WrVxeUSqXUUdSuc+fOwsCBA1WWffrpp0JgYKBEiTQnKytLMDY2Fg4fPqyyvHHjxsLUqVMlSqUZAIT9+/eLz5VKpeDk5CQsXbpUXJaWliaYm5sLP//8swQJ1YN7iHRUQUEBdu7ciczMTHh5eUkdR2OCg4PRuXNneHt7Sx1F4+7evQsXFxe4u7sjMDAQiYmJUkfSiIMHD6Jp06b4/PPP4eDggEaNGmHjxo1SxyoVeXl5+OmnnzBw4ECtuVG0OrVs2RJRUVG4c+cOAODq1as4c+YMAgICJE6mfvn5+SgoKECZMmVUlltYWOjt3t1CCQkJSEpKUvnvsq2tLVq0aIFz585JmOy/4c1ddcy1a9fg5eWFnJwclC1bFvv370fdunWljqURO3fuxOXLl/XmM/l3adGiBbZu3YpatWrh+fPnmD17Ntq0aYPr16/D2tpa6nhq9eDBA6xbtw7jxo3DlClTEBMTg1GjRsHMzAxBQUFSx9OoAwcOIC0tDV9++aXUUTRi0qRJkMvlqF27NoyNjVFQUID58+cjMDBQ6mhqZ21tDS8vL8ydOxd16tSBo6Mjfv75Z5w7dw41atSQOp5GJSUlAQAcHR1Vljs6OorrdBELkY6pVasW4uLikJ6ejr179yIoKAjR0dF6V4oeP36M0aNHIzIyssi/wPTRX/8FXb9+fbRo0QJubm7YvXs3Bg0aJGEy9VMqlWjatCkWLFgAAGjUqBGuX7+O9evX630h2rx5MwICAuDi4iJ1FI3YvXs3QkNDERYWBg8PD8TFxWHMmDFwcXHRy5/tjh07MHDgQFSqVAnGxsZo3LgxevfujUuXLkkdjd4DPzLTMWZmZqhRowaaNGmChQsXokGDBli5cqXUsdTu0qVLePHiBRo3bgwTExOYmJggOjoaq1atgomJCQoKCqSOqFF2dnb44IMPcO/ePamjqJ2zs3ORAl+nTh29/Yiw0KNHj3Ds2DEMHjxY6igaM2HCBEyaNAm9evWCp6cn+vXrh7Fjx2LhwoVSR9OI6tWrIzo6GhkZGXj8+DEuXrwIhUIBd3d3qaNplJOTEwAgOTlZZXlycrK4ThexEOk4pVKJ3NxcqWOoXceOHXHt2jXExcWJj6ZNmyIwMBBxcXEwNjaWOqJGZWRk4P79+3B2dpY6itq1atUKt2/fVll2584duLm5SZSodGzZsgUODg7o3Lmz1FE0JisrC0ZGqn+tGBsbQ6lUSpSodFhZWcHZ2Rmpqak4evQounbtKnUkjapWrRqcnJwQFRUlLpPL5bhw4YJOH9PKj8x0yOTJkxEQEIAqVarg9evXCAsLw8mTJ3H06FGpo6mdtbU16tWrp7LMysoK5cuXL7JcH4wfPx5dunSBm5sbnj17hpkzZ8LY2Bi9e/eWOprajR07Fi1btsSCBQvwxRdf4OLFi9iwYQM2bNggdTSNUSqV2LJlC4KCgmBior//2e3SpQvmz5+PKlWqwMPDA1euXMHy5csxcOBAqaNpxNGjRyEIAmrVqoV79+5hwoQJqF27NgYMGCB1tP8sIyNDZQ91QkIC4uLiYG9vjypVqmDMmDGYN28eatasiWrVqmH69OlwcXFBt27dpAv9X0l9mhsV38CBAwU3NzfBzMxMqFixotCxY0chIiJC6lilRp9Pu+/Zs6fg7OwsmJmZCZUqVRJ69uwp3Lt3T+pYGnPo0CGhXr16grm5uVC7dm1hw4YNUkfSqKNHjwoAhNu3b0sdRaPkcrkwevRooUqVKkKZMmUEd3d3YerUqUJubq7U0TRi165dgru7u2BmZiY4OTkJwcHBQlpamtSx1OLEiRMCgCKPoKAgQRDenHo/ffp0wdHRUTA3Nxc6duyo87/fMkHQ00uIEhERERUTjyEiIiIig8dCRERERAaPhYiIiIgMHgsRERERGTwWIiIiIjJ4LERERERk8FiIiIiIyOCxEBEREZHBYyEiItIi06dPx9ChQ9X2fnl5eahatSpiY2PV9p5E+oiFiIiKkMlk73zMmjVL6ohqV7VqVYSEhEiaISkpCStXrsTUqVPFZZmZmejVqxecnZ3Ru3dvZGVlFXnNyJEj4e7uDnNzc7i6uqJLly7ijTfNzMwwfvx4TJw4sVRnIdI1LEREVMTz58/FR0hICGxsbFSWjR8/XuqIxSIIAvLz80v1e+bl5b33azdt2oSWLVvCzc1NXBYSEoKyZcsiIiICFhYWKqXt4cOHaNKkCY4fP46lS5fi2rVrCA8PR/v27REcHCxuFxgYiDNnzuDGjRvvnY1I37EQEVERTk5O4sPW1hYymUxl2c6dO1GnTh2UKVMGtWvXxtq1a8XXPnz4EDKZDLt370abNm1gYWGBZs2a4c6dO4iJiUHTpk1RtmxZBAQE4OXLl+LrvvzyS3Tr1g2zZ89GxYoVYWNjg+HDh6sUDKVSiYULF6JatWqwsLBAgwYNsHfvXnH9yZMnIZPJcOTIETRp0gTm5uY4c+YM7t+/j65du8LR0RFly5ZFs2bNcOzYMfF17dq1w6NHjzB27FhxLxgAzJo1Cw0bNlT5swkJCUHVqlWL5J4/fz5cXFxQq1YtAMDjx4/xxRdfwM7ODvb29ujatSsePnz4zj/3nTt3okuXLirLUlNT8cEHH8DT0xO1a9dGWlqauO7rr7+GTCbDxYsX0aNHD3zwwQfw8PDAuHHjcP78eXG7cuXKoVWrVti5c+c7vz+RIWMhIqISCQ0NxYwZMzB//nzEx8djwYIFmD59OrZt26ay3cyZMzFt2jRcvnwZJiYm6NOnD7799lusXLkSp0+fxr179zBjxgyV10RFRSE+Ph4nT57Ezz//jH379mH27Nni+oULF2L79u1Yv349bty4gbFjx6Jv376Ijo5WeZ9JkyZh0aJFiI+PR/369ZGRkYFOnTohKioKV65cgb+/P7p06YLExEQAwL59+1C5cmXMmTNH3AtWElFRUbh9+zYiIyNx+PBhKBQK+Pn5wdraGqdPn8bvv/+OsmXLwt/f/x/3IKWkpODmzZto2rSpyvIRI0bghx9+gKmpKbZs2YLRo0eL24eHhyM4OBhWVlZF3s/Ozk7lefPmzXH69OkSzUVkUAQionfYsmWLYGtrKz6vXr26EBYWprLN3LlzBS8vL0EQBCEhIUEAIGzatElc//PPPwsAhKioKHHZwoULhVq1aonPg4KCBHt7eyEzM1Nctm7dOqFs2bJCQUGBkJOTI1haWgpnz55V+d6DBg0SevfuLQiCIJw4cUIAIBw4cOBf5/Lw8BBWr14tPndzcxNWrFihss3MmTOFBg0aqCxbsWKF4ObmppLb0dFRyM3NFZft2LFDqFWrlqBUKsVlubm5goWFhXD06NG35rly5YoAQEhMTCyyrqCgQHj+/LnK+124cEEAIOzbt+9fZxUEQVi5cqVQtWrVYm1LZIhMJG1jRKRTMjMzcf/+fQwaNAhDhgwRl+fn58PW1lZl2/r164tfOzo6AgA8PT1Vlr148ULlNQ0aNIClpaX43MvLCxkZGXj8+DEyMjKQlZUFHx8fldfk5eWhUaNGKsv+vpclIyMDs2bNwq+//ornz58jPz8f2dnZ4h6i/8rT0xNmZmbi86tXr+LevXuwtrZW2S4nJwf3799/63tkZ2cDAMqUKVNknZGREZycnFSWCYJQoowWFhZFDsgmov/DQkRExZaRkQEA2LhxI1q0aKGyztjYWOW5qamp+HXhMTl/X6ZUKkv8vX/99VdUqlRJZZ25ubnK879/hDR+/HhERkZi2bJlqFGjBiwsLPDZZ5/96wHQRkZGRYqHQqEost3fv19GRgaaNGmC0NDQIttWrFjxrd+rQoUKAN4cM/RP2/xVzZo1IZPJcOvWrX/dFnjzEVtx3pfIULEQEVGxOTo6wsXFBQ8ePEBgYKDa3//q1avIzs6GhYUFAOD8+fMoW7YsXF1dYW9vD3NzcyQmJqJt27Ylet/ff/8dX375Jbp37w7gTWH5+wHOZmZmKCgoUFlWsWJFJCUlQRAEsdTFxcX96/dr3Lgxdu3aBQcHB9jY2BQrY/Xq1WFjY4ObN2/igw8++Nft7e3t4efnhzVr1mDUqFFFSllaWprKcUTXr18vsieNiP4PD6omohKZPXs2Fi5ciFWrVuHOnTu4du0atmzZguXLl//n987Ly8OgQYNw8+ZN/Pbbb5g5cyZGjBgBIyMjWFtbY/z48Rg7diy2bduG+/fv4/Lly1i9enWRA7r/rmbNmti3bx/i4uJw9epV9OnTp8jeqapVq+LUqVN4+vQpXr16BeDN2WcvX77EkiVLcP/+faxZswZHjhz51zkCAwNRoUIFdO3aFadPn0ZCQgJOnjyJUaNG4cmTJ299jZGREby9vXHmzJli/mkBa9asQUFBAZo3b47//e9/uHv3LuLj47Fq1Sp4eXmpbHv69Gn4+voW+72JDA0LERGVyODBg7Fp0yZs2bIFnp6eaNu2LbZu3Ypq1ar95/fu2LEjatasiY8++gg9e/bEJ598onIRyLlz52L69OlYuHAh6tSpA39/f/z666//+r2XL1+OcuXKoWXLlujSpQv8/PzQuHFjlW3mzJmDhw8fonr16uJHS3Xq1MHatWuxZs0aNGjQABcvXizWNZgsLS1x6tQpVKlSBZ9++inq1KmDQYMGIScn5517jAYPHoydO3cW+6NEd3d3XL58Ge3bt8c333yDevXqwcfHB1FRUVi3bp243blz55Ceno7PPvusWO9LZIhkQkmPzCMi0oAvv/wSaWlpOHDggNRRJCMIAlq0aIGxY8eid+/eanvfnj17okGDBpgyZYra3pNI33APERGRlpDJZNiwYYNar66dl5cHT09PjB07Vm3vSaSPuIeIiLQC9xARkZRYiIiIiMjg8SMzIiIiMngsRERERGTwWIiIiIjI4LEQERERkcFjISIiIiKDx0JEREREBo+FiIiIiAweCxEREREZvP8HO4AuvoND5sEAAAAASUVORK5CYII=\n"},"metadata":{}}]},{"cell_type":"code","source":[],"metadata":{"id":"3lkcY5lYhr37"},"execution_count":null,"outputs":[]}]}