{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Intro to R\n", "[NetCDF](#NetCDF) | [Time slices](#Time-slices) | [Multiple variables](#Multiple-variables) |[Multiple sites](#Multiple-sites) | [Using functions](#Using-functions) " ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "As a demo of the language, plot sin and cos on the same axis:" ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "collapsed": true }, "outputs": [], "source": [ "# set up the default plot height and width (this is just for aesthetics in the jupyter notebook, you won't need it in RStudio)\n", "library('repr')\n", "options(repr.plot.width=6, repr.plot.height=5)" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "collapsed": false }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAtAAAAJYCAYAAABLtNEpAAAACXBIWXMAABJ0AAASdAHeZh94AAAgAElEQVR4nOzdd7hdZ3Xn8e+2LcmyLcu9V8mWblFzw7iB6RDIACFM2gCBkErNJEwypEwmTwqZJKSTDISQkEwKk4SEoQaDwb1LtqV7JVmWLVdZlovkbtle88fax/fq+pZT9t7vu875fZ5Hj23p3rOXr96zz9rvu971goiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIi0rUgdgAD+93A+sDB1ICIiIiI1eQq4MnUQVVACnYezgetTByEiIiJSs3OAG1IH0av9UgfQoHOB/wKMAseWv+YD95W/rgI+A9yeILb55T8XAM8muL6IiIhIneYDzzCR80jmCuBPAZv06zFgG3BX+e+T/+zTNP9gcX55bQ0qERER6Ufz8Vzn/NSBSHs+iP+FrQX+M3DINF9zMPA24LLya3++seicEmgRERHpZ0qgg7kO2Aoc2MbXziu//opaI3opJdAiIiLSz/oqgd4ndQANGMYT4ifa+No9wLeBlbVGJCIiIiJhDUICPQZcABzQxtfOAy4Gbq0zIBERERGJaxAS6M8DS/C+g+9k+hroRXgN9Lfwbh1faiw6EREREZHMFMAfA88z0WljN96F405gF3t34fgrfCa6SaqBFhERkX7WVzXQg+QU4HeBS4FNePu6x4EteI30J4DTE8WmBFpERET6mRJoqZwSaBEREelnfZVAD0INtIiIiIhIZZRAi4iIiIh0oOkjq3N3KHB1+e9DPbxOgbfO27/Nrx/t4VozMhgB7izgyTpev3p2EHA4FNtSRyLpmW/mXQ5sDTSGDwTmQfFo6kgkPYN9gWOA3YXvu8mYFXiXqsegeC51NCISyxFMdOPoxVJgJ/Bwm7+eKK/ZzmmJbTGYZ/CowRMG/2bww1W9drXsRLCvgT0EZuWvK8FenzoyScPgJw2uNHiyHBCPGPyxwYrUsU3PVoNdPWn8Gti3wGp5MJa8Gexn8OsGdxg8O2kMf9zgoNTxTc9eDXb9pPG7C+yfwY5OHZk0z+AQg0+VOUSVE619VQMtefhJfFBVenM1ONDg7QafM3jK4A/MZ8czYW8EexDsMrD/DHYm2Eqwz4A9U/7+YamjlGYYzDf4TPnQ9+sG32twqsF7DG4weMHgI6njnGD7g/0W2LNgXwI7D+wssAvBvlz+/if962QQGJxkcLnBg+WD4CsNhg0+YrDd4AGDN6eOc4IdCfaVcqx+CmwF2Plg7ygT6h1gb08dpTTH4J3lWN1s8NqKX14JtFSulgR6MoOLDXYZ/GW5rJiY/RzYc2CfAJvmCdeOB7sB7HKwhc3HJ00yONzguwZ3G5w5w9f8pMHzBj/edHzTRFOA/QvYdn/4m/Zr3gJ2L9jflcvj0scMVhs8ZHCpwfHT/PmB5UrKkwYvTxHjlIgWgl1VJsrLp/nzeWC/DfY82H9vPj5pmsGby1WTX7P2S1A7oQRaKld7Ao1f4MxyBuRv67xOG5FcXCbPPzjH1x0FtgXsi2AZJP1SB4PC4GsGNxocO8fX/rTBcwb/pan4Zojkv4E9DrZqjq87o/y6X2wmLknB4OByxu4fZ5ugKMf65w12GJzWZIxTItkH7Atgd8xdpmE/XN6vX91MbJKCwVBZsvHxGi+jBFoq10gCjV/kjPIJc47ktbYIjgC7B+yP2vz6pWAPgP1ZvXFJKmVSvMvg1Da//iPlGE5UE22vBtsD9gNtfv33lwnIW+uNS1Ix+KcygV7UxtfON7jE4DbzjesJ2CfAHgEbafPr/wTsfrBj6o1LUihrnjcZ/EvNZZ5KoKVyjSXQ+IV+o6xxari+2Aqwfwdb11ldqF1ULiNeWF9skoLB6QaPG7yvw+/7isFlzdf023HlA93vdvh9vwq2WwlI/zH4GYOnZyo9muF7FpcJS5sTCVWyV5YPdK/p4Hvmg12Db47VamCfMfiSwYZ2HgB7pARaKtd0Ar3AYNzgc01cb9KVf6Jczu6iRaD9LdiNvvQo/cC8W8HVBv/WxfcuKWtJf7SG0Ga78mfLcdhhEmH7gN0M9hf1xCUpGJxgvkH7g11875vKlZThOmKb4ar7lOP3L7v43pPKWesPVB+XpGLwhrIsromuQUqgg9nY5a8mNZpA4xe7oHzTVL3LdqYrHgB2H9jPdfn9x4E9BvZj1cYlqRj8VLnp6qguv/9XylrShlZSbKScuXtdl9//hrL0o5ce85IRgz83uLbblZCy9v8rVcc1yxXfU95HZ91rMMv3f7y8jx9QbVySgsE+BmsNPtPQJZVAB/NJvCezdfirSY0n0PgF/6J88zSwDG6/AHY32IIeXuPn8H7Rh1cXl6Rg3qf8ToNf7uE1FpTL4A3N6tqXwb7W42t8FazBhEnqYt5i8RmDN/TwGiMGewzeWGVsM1zt4LL8qIcNrbY/vofll6qLS1Ix+BHztqEv6RpTEyXQAR0H3Iv/xeXYkzVVAn1yuYRY883bDixv3B/u8XUWgG0G++1q4pJUDN5vfrjE4h5fp9V26cSqYpvhSueBvQB2do+vc0ZZz39RNXFJKgZ/ZXB5Ba/zJ2X9ac0nA9uvgd1Fz21B7SPlREZP711Jq5yAuMPgtxq8rBLooH4DJdAvYfA3BpfWfJXW7HMFP3t7d1mHl+mJXjKXsvb5doNfq+C1CoP1Br9TQWizXem7YF+o6LU+D3ZlNa8lKRgsK2eOX1nBax1ezgK+o4rYZrjKIvx0wfdU8Fr7l4n4/+z9tSQV825GO3udxOiQEuigvg8l0C9hMGp+wltNjf1tEdhOsJ+u6PUW4IdXaCNLUOYnC+6yilp4GbyvnM2u6f1j55WzxtMcNtHV6y0rZ7PPqeb1pGkG/2DwzQpf738bfKeq15vmCh/A29DNr+j13o93lUnUhk96UdY+bzVo+oAcJdBBzQOOIKujrF+ULIHGL/zvBl+s6dU/VM4+V3TjBp/5sE2oI0c4Bvua98v9zQpfc4HB/QYfquo1p1zhb7x2udLX/AbYX1f7mtIEgxPLDdgXV/iaK8uJjNVVveakVy/Axr2Eo7LXnIdvJqzpPSd1KkvfnjY4suFLK4GWyqVOoM83PyK54nZKL964f6Xi1z0S7Cmw/1Tt60rdzFsmPWX+MFvl6/6ywRar/Jh6OwrsGbDvqfh131y+7hynwEluDP6HwcaqN1+bHwH+2Spfs3zlN5X3y4qTJfv18v6e46SUzMLgywZ/l+DSSqClckkTaPziVxj8YcWvegHe9quGHb72N2Dfqv51pU7mJ7b9Yw2ve0RZR/q2il/5F8C2Vr/aYfuA3Q7W9BKq9KBcQdlm8PM1vPb31fFw6asn9vlqXxPATizv79oQG0jZPeZ5S5PEKoGWyuWQQP+YeU/deRW+6l+Dfbm619vrtc8AM7AaljylDuVmqaeth7Zfc7z+pw2+XeEr7lMmzz20/Zr19T8Gtg2d7BaG+eEnz1iXvcvneO39yuT84xW+6vKy3r7tUxI7fP2aknOpi8EnDNYlurwSaKlcDgn0weYnu1VUFmGL8Ib931/N6017jWvBfq++15cqGXzY4O7qyyxefP2LypmVEyp6xdeDPUttx2/bkWBPg72lnteXqhn8q0FF3Vimff1fNO+PXlFZhP0vsOurea1pX//7wZ5ALe1CKPeL7DDPOVJQAi2VS55A4wH8g8E/V/RqPwG2o9rNgy+5xgfLDYraTBiA+aE9v1Hj6xdl8vGxil7xX8H+qZrXmvEafwfW8VHm0jyDY8x7jr++5ms8Z3BhBa+2T3l/rLFjkc1DXZHCMPghg0fr61g0JyXQUrlcEug3lUvsFRyNbNfXPztsR+NHI7+y3utIrwzOLLsMnFbzdX7LYG0Fr3RIucmv7kOG3lzOQmsGL3Pl7PBWg1of2A0uMfjTCl7porJGueaNqvYJsFQlAdKBsuPWXycMQQm0VC6XBHq/sh3Yz/T4SsNlfXLFXT2mvdY3wP68/utILwz+1OCyBq4zYj74Rnt8pffgp61VuCdg2uvMB3sY7N31Xkd6ZX5aYMUdhaa9zvsNtlvPJxPan4FV1qt6lusMlff7kfqvJd0yWFxOkL05YRhKoKVyWSTQeBC/b3BNj6/yq2A3VhPRnNd6V5mALGjmetKp8sFsp8F7G7reTdZzn2n7Kthnqolozmt9hsr7TEuVDIbKB7PaJwUMDis3Kr6uh1eZB/YgWCPvOZ+BrrLPtFTN4F0GD5snsakogZbK5ZRAryk/KJb18CrrwH6puqhmvdYisCfBvreZ60mnDF5V1o42cmqZwc+VS+1dbsSyw8rNg7XVuk653uvK6x3ezPWkUwa/ZLChwev9P+upJ7S9sSwNOqS6qGa93q+A3drMtaQbBl8y+FziMJRAS+WySaDhxaXK/9bldy8tl/N6SMA7vuYXwP6+uetJJwz+yOA/GrzeceVGrAu6fIUfLTfA9riE3vb19gV7AOz9zVxPOlWuavx6g9f7YfPj6btcWbPPgf17tVHNer0Gy/akU5PKN96UOBQl0FK53BLo3za4ssvv/hhYYzM15TXfjrdSyuLnJxPKzhjbeq+r7/i6lxr8UZff/RWwv6w2ojmv+edgjT1kSPsMlpSrco31nDc4yPxgoLd28d37g+0C++HqI5v1uhvAfrnZa0o7DN6dQfkGKIGWGuSWQJ9n3k+3i8MC7Aqw36o+qlmvuRDscbCKT6GTXhmcVXbfqOE0ylmv+7MGt3fxnYvLpe+Kj+6e87qvVRlHngx+3uC2BNf9J4P/08V3vqUsa1tUfVSzXvc3wG5o9prSjrIk6K9Sx4ESaKlBbgn0Pgb3Gfx4h995HNjzYGfXE9ms1/53sE83f12ZTdlW7qoE123NGq7s8Dtbm1Ibnqmx/cqyEZVxZMbgOoOGJwVo9ex9uPNuHPZXYF+qJ6pZr3sGfurhqc1fW2ZSbkqttX95B/oqgdYBFPISBbwAfBXo9IS0twF3Aw114NjLV4HvAavoBC+pyPcCjX+YF7AV2ETnNX/vAL4ExbPVRzWb4jngy/jPSzJhfqrl2UCKw26+CRwMnNv+t1iBj/kEXV2Ktfj77h3NX1tm8SbgMeDS1IGI1CGrGWgAg/9kfrT3AR181yVgn6wvqlmvfUI5+3FGmuvLVAanl7PAQ4mu/3sG3+3gOxaUpUDvrC+qWa//TrDHUEvGbBh8wOCu7ju69Hz9qzpryWhnlZv5Tqovqlmv/7tgtfd7l/YZ/K3BP6aOo9RXM9CShxwT6IUGjxt8X5vfsRg/FfDiWgObPYZ12sSSD/OT29YnvP5FZTeOI9r8jteUY7iRdnvTXP+Q8vqvTXN9mcrgKwafSnj9jxvc0sF3/DJYB19fNXt1Wct/cLoYpMVgX4MHDd6VOpZSXyXQKuGQaRXwFPAt2i/jeB3wBHBFbUHN7f+hJfCcvBkvS0jlamA37df+fQ9wFRSP1BfSbIpH8e43KU8Kk1LZseAVeClFKl8HVpSlJO14C34fTOVy/LOjh0NgpEJnAoeTdgz3LSXQMpsvAW+x9sbJ64HvlLWcqXwFOBvs6IQxCGATtZvJWrMV8ByegLT7EPjG8utT+hrpe7WKuwDYH/h2whjWAvfT1piwo4BzSFL/3FLswWtt35AuBpnkjcC6AranDkSkLtmVcAAYHFO2IDuzja++G6zRXr/TxLAPfiDFe9LGIWUN/RPdHwRRWRw/ZPDQ3J0M7OSydjRxDb2tLONYmjYOMfiE+Yxq6jj+0uCLbXzlu8B2gu1bf1SzxvHT/nkgqRlcmaKDzCxUwiGDoXxqvRWYoybThvAlxm/UH9VsihfKGLQEnt7rgMsKeCZxHF/HZ8PnumG/AbgPWFd7RLMqbgXuQjN4OXg9ye9pgM8ov7aNh9E3Af8BxfMNxDSbrwMngI0kjmOgGRyKrwKmXlXrW0qgZS7fZO4a0tcCd0DRxcEVlfs68FoaO4ZZZvBGMkg+CngEuJ65E9LXAd/ygxOTa+c9JzUyOA5Yg5fUpPYNYB5w8cxfYvviYzj5ew6KO/BDjFQHndYbgV10faqwzEUJtMzlUuB881rAmbyehLWuU3wTOATv3SoJGJwEnIZvQs3Bd5g9+dgPfwjMIVkCj+M1zR/mIpO8BngIr0FOqvDN2Vcy+0Pg2cBh5DOGv4FWUVJ7PXBpAalXJPqWEmiZy2V4/egMS+A2H3gV2STQxYN42cmrUkcywF6Hb3xK1sJuim8D5xjMdLTxufif5bJT/RJ8ub4v6gSDegNwSXmoVA4uZfZ72sXALVDsaCacOX0deCXYbBMvUpOyb/kbUPlGrZRAy6wKP8HoWmaugz4fn53OZbYR5v6wkXq9Dk8+ciiHAJ+9ewG4aIY/fwNwIxQ7mwtpNsUu4AZmnTWXupRdhzIph3jRpcAq85Zk07kYX2nJxaX4xMtM7zmp1yhwLPlMCvQlJdDSjm8xcwL9euDa8kM/F5cC52sJvHll8vEaMrpxlz3NrwZePcOXXEzaVmXT+S5KoFNZAxxJNqtqAFyHj+NXvPSPbB5wIVkd1Vw8DlzFnBvQpSYXA3cUsC11IP1MCbS04xLgzHJX71SvJqNkqfRdfFb8ZakDGUBn4LNkl6QOZIpvM20CbQvxcdLBkd+N+A5wrpbAk3g1sLHwrixZKGAPfkjVdCtrZwMLyaDl3hTfYdqEXxrwSvJakehLSqClHdfisx9Tbt62CDiLrGY+oDzR7UZ8JlSa9TpgfeE10Dn5NrB6miXw8/H7YMoTNKdzOR7XBakDGUCvIL8HKvCEaLoE+jXA2nQnaM7oUvxgq5n2HkgNyvrnV6AEunZKoGVO5ezHd3hpa63z8dPerm06pjZ8G9VBp5DjigT4GH2cl85CvxK4oVxyzkjxJN5+7+LEgQwUg33JN/m4BBg1OGbK719MdpMYgL/nnkObYZu2AjiK/MrS+o4SaGnXdBvzXgFcB0XqwzKmcylwHtgBqQMZFGXy8XK89jEr5bHeV/LShPSV5DnbCKqDTmEFsJj8ViTAD/l5jL3KIqzVrSXDBLp4Bk+iVcbRrIvw+ud7UgfS75RAS7u+AywrDxhoeQXe5i5HV+BdIM5LHcgAORM/jj7XhPQS9trUZAvxFna5xvtd4GVgB6YOZIC8CritgHtTBzJV+RD4XfYuTTsX73aRW/1zy3fxh1RpzqvQ7HMjlEBLu27Bl8DL5ThbCJxDtjfuF5fAVcbRnAuAzQVk0g7uJS5l74fAl5Fn/XNLKy49BDbnQvIdD+AJ6eQZ3YuBm/IrQXrRZcA5WglsRln/fBH5Tmz1FSXQ0pZy9uMqJvp6noPPfFyTLKi5afajWTmXQ4AvgT/MxBi+yH8v1+SjeBK4CS2BN6JMPnLvXvAtYMjghPK/c3/PXYP/XNURqRkjwNFoBroRSqClE1cxsSHkFfjO790J45lLa/ZDrcBqViYfF+B1xlkqD3a5lokZ3dyTD5jzGHKp0BBwBNmuqgF+yuqjwAVl/fN5ZJ3wF0/gHZE0kdGMV6D658YogZZOXAacUR6J/Ery/qABn/3YD++TKvUawQ+fyD0hvQK4qDxk53zyj7dVB70wdSAD4FXAtgLuSB3ITAp4Hn9IvQC/ry0gw027U0wtO5H6qP65QUqgpRPXAS/s5uAL8JmPzOusit3ABtRGqQkXAPcGOPnqamDVCBsuxJOP3B8Cr8IfAs9KHcgAuIC8659brsHvv+cD6zM7BXY6lwEv18mwjci9hr+vKIGWthXwBHDjNk5+J3AA+Scf4LM1SqDrF6EcAjz5eOF1XPKDhEg+XnwI1EbC+r2KrMshXnQFsOZQHnkF+c8+g8e7AN83IzUxWAocixLoxiiBlk5dtZCnLgI2QPFQ6mDaUNZtW5E6kD6Xdf1zS+Enat68jM0X4rPREVyNEuhaGSzBk48IkwLXA5zNDRcQIoEuduO12xemjqTPnYt3QLo9dSCDQgm0dOqyE7n71P15OspT7lV4be5pqQPpVwanACcTYwYa4PJhxpcSJ4G+CiXQdbsAeKCATakDmUsBTzzFwrFzufZQ4ozhq/FDlqQ+FwBXlZulpQFKoKVTV8/n2f1+jM9md9DA9Iqt+KEIF831ldK18/H2cOOpA2nH1Zx3x5ncNH8Rj+XcgnGyK4FjwPQQWJ+XUc7sRnAzq3dcxOXPQLEldSxtmtzBSepxLt5lSBqiBFo6UmAHbGSIj/KHkboCaAm8XucBVxbwQupA2vFR/pDF7OIhDg9ywl9xO7AdjeE6nU+Icgj3r3wfF3DlPhbnM/wa4CiwpakD6UfmJ8CuIUAZXT+J8uaTfJx3A2c/dRpbTk8dSAeuwpe3pB4vwzu0hHAdLxu5n2OfnMeeSEvKrc4LUjGDhcBKAo3hf+EdJx7IE/Pw3tUR3A48iMZwXc7CSzduTB3IIFECLZ06d4yRzcS6EV4JDIEdljqQflMmH2cQpxYT4LxtnLyJWGU9V6Ml8Lqcgx8EFKSkxw7dypJlT7P/dsJszCsM//lGemiN5HzglgIyPVW1PymBlk69fCNDlwMnmO9aj2At3n0hUtIfxRpgX8LMfNgiYOXjHPQdYo2Hq4EVYAenDqQPvQwYL9t0RnAu8Mw89lxOrIRUCXR9VP+cgBJo6YDtD5zxNd7478BuwtwMiz3ADWgGrw7n4cnHo6kDadM5gB3PvV8ATjE4LnVAbboBP4VOvXSrdz6xakfPB27cl+cvI94qyiqwA1IH0k/MV08uJNYY7gtKoKUTa4D9nmH/6/FZ3UhHZF9LmIQ/lHMI1L0An6nZMMz4jfiqRJAT/oqngFvQGK7DuQSqf8bHwDV4QrrU4IjE8bTrBjznODN1IH1mCXA4moFunBJo6cTLgY3l6W3X4R88UVwPnAWmMV+tcwlTOwr4jPnVBezBHwIjzeiqm0zFDE7EVyGCjGHbB3/PXY0/UD1DmImM4jF0qmYdzgN2FBClpWHfUDIhnWjNfIA/7Z4dqI3S9cBiYFnqQPqFwdHAqYSZ+bACH8OtDY/X4/WvUZSHUehUzQqdBzxCgANUSsuBg4Hry4fAdcR6CFQddPXOI8wDYH+JkvxIHiYn0NfhCenydOF0orgTb6MU6cMmd+fiG6/Wpw6kTUvwUylbCfR1+ENglIT0enypVr10q/Ny4LooPczx2eYdUNxV/vf1xLqn6UTC6kUrQeobSqClTXYsflzzNQAF3I2f8BfpZhjtwyZ35+EzYc+lDqRNZwEPMbHUeS2ekEY54e92fLY0yJJ9CKEOUOGlew6uI9Y9+BrgOLCTUwfSD8o2oqtQAp2EEmhp18uBx4CxSb8XLSGNFm/uziHWjfts4MayJy3AVvwI8iAJaWH4RiyN4QoYzAdWE2sT7NnsHe+NwJFlLXcEm/EOThrD1VgJ7AfclDqQQaQEWtp1Lp58PD/p96LVs10LnAG2IHUg0Zn3fm5tZopirxMTCz+561pibYa9jlh12zlbCexPmATa5uOHFk1+aN2It5AMMiaKF9BDYJVeBtxe+MqaNEwJtLTrLPzGN9kNwIpyGSmC64EFwIrUgfSBIeAgXjomMmX7MP0YvpEwM9CAx38m2H6pA+kDZwJ3FbAzdSBtWsGUhL98CIzWUvQGYsWbszPR7HMySqClDVbgycfU0+auw8dQlF66O4E7CTNbk7VzgPsLuCd1IG1azvQJ/7XAWQbzmg+pK9cDBwDDqQPpA9F6mJ8D3FnexyaL1lL0RrylaJTNuznba1VNmqUEWtqxBDiUKQl04TXRm4g1m3AdWj6swpn4zFcUZwP3Q3HvlN+/AZ/VG20+pG4U9wL3Ees9l6tos3dT659bbgTOCNRN5ga8g1OUzbtZMjgQXwmcOrElDVECLe04C6+zm65Re7ReutpIWI2ZPsxzdTbTlJsUsB3YRqxaftWQ9sj8oWkVscbwTDPm1wKHEKalKHfgNbtBVi6z1TrRMUgZXf9RAi3tOAu4aVL3gsluItbRrNcDw2AHpQ4kKvNd36uJNfNRduCYVrRVCdWQ9m4VPo6DjGFbiK+SvCSBLuAu4AHCjOHCeLGMQ3pwDrCpgMdTBzKolEBLO2ZLPtYCp5vXl0bQ+v84I2kUsS3H63CDlHDYfsAaZp6puZ54D4Gr1U2mJ2cC2wpvYxjBGvzzeqb78E3EeqjSQ2DvopUg9R0l0DIHK/Cb90xv1LX4TvAgCUjxON6LVDfv7p1NrA2EQ3jCP1PyEa2bzE14D+MgddtZilaCdBawGYrHZvjzaO0Y1wJrtJGwJ9pAmJgSaJnLEuAwZkg+yuWjLXiSHcV1hEn4s3QGsC51EB04G7gHiu0z/Pk6vK91kIS02IHXbQdZss9StNm7c5i91vUmYFVZXhXBDTstFS8AACAASURBVHjdtjYSdsEmNmEGKUHqT0qgZS6zbSBsibYcF61uOzdzfZjnZtoNhC2FH4+9lVhjONp7LhvlSsNKYs1An8ns77kb8P+vkWbC6VVxJ95/W2O4O+cAzxHrIbDvKIGWucy2gbBlLbFqitcCy8EOSB1INOUJhKuIdeOeNYEurcM3RkZxA7HeczlZiY/jKDX8C/EypBlXfQq4H99IGG0MayNhd84Exgt4OnUgg0wJtMxlugNUproBGC77UkZwM94zdVXqQAKKdgJhux1Dom3CWguM6kTCrrSOP46ygXAVnvDPVTZ1E7ES0mingOZkulNVpWGDkEBfiM84TFYA3w98C+8D+2D57z9CnGb0DZjxBMKp1uI/tyCzH8Vu4HZUxtGNs4DtgTYQDuM9f+eabbwJWBnoRMK1+P+XTiTsXLTk4wxgKxS75vi6aCuB5SqKDUIeUrU1hFlB6V+DMHAvB/5gyu99Fvi/wKuBg4Fny3//O+AfUBLdshTf6DHrh00Bu/Hm+JFu3tE+bHIRbQPhGuA+KB6Y4+vWAQsIk5AWO/ATCTWGO7eGWGP4DNpLltYBqwOdSLgW//xdmjqQSMqV3tPwlVRJaBAS6Km+H3gvvoHhrfhy9PHAUcA/Aj8AvCdZdHk5i4kNVnO5nnjLh5HizUW0DYRtdVsoTyS8l1hLyjeiVZSOlCcQriBWDX87q4CUX7OYMAlpsQ3/HNYY7kyr45VmoBMbxAT67eU/PwR8CXih/O8HgR/FZ1KVQLvVwM1zbCBsWUes2bBWDakOo2iT+f0i2gbCTmbMo20k1CpK54bwVm9BZu9sHp7wt5Ms3YGvBkZqKXozsd5zOVgNbNUJhOkNYgI9VP7zP6b5s2eAK9AbuqWTpU4dRtH/TgcWEWYG2lp1+e3O1ERrDbcOHUbRqTOBewvYkTqQNo3gpUVzjuHCD7SKtpEw2sRLDqL1MO9bg5hA31n+8+AZ/nw+3vdYOkug1+EzOyvqC6dKxUPA3ejm3YlVwCOF/9wiOAWv4e9kDK+yOPfFVg3pktSBBLKKMLPPgN+f7itr3tsRbUb3ZmLNmOdgNbFq+PtWlA+KXh0PvAn/QP3j8vfeO8PXfQ863QewY4BjafONWsBD+OlokW6GWgLvTLQb9xn4w/AdbX79WnxPxLLaIqrWNrwVm8Zw+/p1A2HLTcQaD2uBY8rPG5mDTUxSRRrDfWtQEugh4Kv4B+nXgOeBX2XvzQs/BlyLL1F/tukAM7QGL2kZ6+B7oiWka9EGlk50+mGe2hrar+Gn8IT0QcKUcRSGz+BFes8lU3aniDaGu0mgjzGfDIpgI34YiMZwe0bwjbCa5MvAICTQh+AfiD8E/ArwBTxR3gEcNunr/jtwOL658OsNx5ijNcAYFHs6+J5beWnP7ZzdBKwG2zd1IEGsxv+Oo+gmWbqZWGM42kNrSqfgJS9BSjhsHzrv99tKSIOsBBbPAesJE29yq4H7A9Xw97VBOMVqF/60Nt0T2+TNNx/BW7FpYLpuluvXAh812KeY6G6Ss7XAAfiS/XjiWLJmcDQ+qxVt9u6fO/yetcT6MF+LTw7I3M7EPw+2pA6kTUvxFdG233MFPGdwC/7/+pW6AqtYuRlW2qANhBkZhBno2Uxe2v0KSp4nW43fiDtxC37DP6XyaGpR3I3XkEbadJPKSmAPYR407Ag84e90tjHaKsrNwLFgR6UOJIBVwPpi7/t+ztbgCf+dHX5ftIT0FnQPblfZWlZyMOgJtEzLWrOync5AbwUeI9bNMNqu9VRWA+OFn9oZwRo81k5q+KFMSA2OrD6kWoyjGtJ2RdtAWK4CtlfDP8k6wnRDAjze08EOTB1Izsoafh3hnZFBKOHoxKHA1eW/D832hXNYCPwEXuzfjnN7uFYdVuIPVx192BRgNtHX84t1BFaDW/CZKZndGcRaOlwDbICi04R/HN88uwb4ZuVRVa54DqxVQ/qN1NFk7kz88KwoVtLdbOM64DSDg4IctnEznhyuYuLzV17qFDxHiXQf7mtKoPe2L7C8gtc5DG+H1+4Mf2vHdC4HIqwBtkHRTT/saEvgtwDvSB1EAKuBz6UOogOdnED4ogL2mG/EWkWIBBqIt2TfOPMN4icQa/l7Fd0l/Bvwz5JRfMN85ordYFvxMawEemar8ZMm223LKTVTAr23nVSTxN4LvKGDr/9J4C/Ipzavl6XOtcCbK4ylbjcDJ4AdBsXDqYPJkflGy2FiLR2uBj7d5fdGS0hvAX4qdRCZOwOv4Q/SRcYWAyfT+T4UCthtXk63hhAJNKATCduxBrg5UA1/31MNtEynlwT6FuAkg8UVxlOnDcBzqIxjNsP46kzHH+Zp2AJ8JanbeCOuoiwDa7dkbBCtBjYVXp4TwUq8k9H6Lr8/2hiO1v0mhRV0Px6kBkqgZQrbB7/xdptAr8dv/EES0uJpYDNh4k1iNXBPedpkBCP46lq3s423AMMG86sLqVa34v+/w6kDyVi07gWrgdugeKrL74+2t+NmYFQ9+We1ijCTGINhkBLoc4E/Ab6NbxR6FHgS7wl6GfAJvO/moFsKHEiXCXThP9PNxFqOi7Zk37Rop7etBu6GotuEfy2ePI9UF1KdioeBu9AYnk03bTlT6qYP/2TrgDWWz76audyMl4qdnjqQHJkfAHQasR4C+94gJNAF8KfANcAHgVfhm0l24TXPRwMXAb+AJ9OfZrBrw1fhDxd39fAatxBr+TDabE3ToiUfK+kh3sLvC/cTa0xoDM/AYB5e0hNp+bvXUz9vxXvyn1xNOHUr7gYeQWN4Jq3VpQ1Jo5C9DEIC/YHy1zrgB/A2MK0by0nlvy8G3g5cDvw48NEkkeZhJd7+q5eNCtGOF74FWAE2yA9O0ypnsKItf1ex1BltVUIJ9MyGgAWE6QFt++D1rr28527HOzZEGsPR6rabtAbYWvjfqWRiEBLod+NtXy4EvoDPrk61G/g34DX4cd5vayy6/PQ0e1e6GVhhcWbyb8Y/YJelDiRDJ+HLh0G6FwC9z95BvNPRlEDPbCWws4D7UgfSptPwcoZeVlEMn62MNCY0hmdWxeeyVGwQEuhh4ArgiTa+dg9eIz3IT8Gr6D35uBVPSIPUsxX3AQ+im/d0VuAn3W1JHUh77Bj8FMFeP2yizYbdAhwBdlzqQDK0glhL36uAh8qyhl5EG8O3onvwTFYSqwRpIAxCAj0GXIA/0c9lHnAxsWbbKmQHAUvo8f+/8D7YO4k1g6eb9/RWAWOFt/qLYBWe8G/q8XXWAkfaxCFHudsMPIXG8HR63ZDXtFVUkyzdTLwSjpPBDk4dSE6CltENhEFIoD+PJ4VXAu8EDpnmaxbhZRvfwrt1RDrutUqj+Ju1igeIDeXrRXEzSj6mM0qsmY9VwLgfcd2T2/CewUHGRPE8PlkQJN5GRZyBriJZuhVYYt5VKYLW506kWfMmnIjv04p0Hx4Ig5BAfwpvX7cKr4F+BK953gbciXfj2A18Ee/G8TngD1IEmoGVwF1Q7KrgtaItH6r+bnpVlPQ0qZJawcLLuVRDGpz5hMmJxKofraKGH/z/ucAfIAIoHsf3K2kM720VEy13JSODkEAb8GG8v/HvAd/BW1QdBhyB175eCfwOvonsffiH5yCq8qSjaJuwbsWP9D40dSC5KA8SGSLW0mGVtYLREtL1hEmWGrMK/wwI8hBoB+MdonqOt/DJoTuJVcaxnlgrl01YA6wv4PnUgcjeonRJqMKdwMdSB5G5lcC1Fb3WBuBkg4MKeLyi16zTOH6C4kr8YB3x3rnzCLP8ba3T+KpKlsaAH6notZpwCzAMNg+KQZ0EmGoF3v7rydSBtGkFnvBX9RAYbSXwVuAVqYPIjI7wztQgzEBL+6pcrm8tmQa5eRdP4r1TNYM3YRWwI1D7r2XA/lS3XN860nteRa9Xt1uYWDUQF23z1WpgKxTtdI1qxzpiraKUCb9FOUGxCdHG8MBQAi0lOxYvaakkgS5nne8kTAINaAl8qmitk1YBD0Jxf0WvdyuekAbpD17sALYTK2Gq20rClG8A1c82Rrun3cJE3frAs4njzSPV8A8MJdDSsgrvOtBr+6/JotWzRVvurFu0pcMqTiB8UTnz/gixxvB6NIaBF9t/jRIvga4y3g3AoQZR+oPfhreh1EOgGwb2JUwZ3WBRAi0tq4GximsnI24k1PLhhGhLh3V0DInWH1wPgRNOxU/RjDSGq54x30ysdozP4clipM+NOq0E7i+82YFkRgm0tNTR7zdaL+gNeL/NE1IHklrZ/ut4Ys18jOIb/6oUbQyPESveOq3AT6C9I3Ug7bHjgEOpcAyXByBtAkaqes0GRFu5rNMI1d/TpCJKoKWljnrXm4EjAi0fbsGXDyPVDNZlJbHafy3C239VPYYjtrI7Sae5AeWJfoV314lgBfAsPmtcpYhjWPdgV0lfe6mHEmihbP81QvVv1NbyYZAl5eI5dJpby0rg9mDtv6D6hH89cIrBQRW/bl1uxR98lIDEPUWz6haEG4g1Hm4FhsDmpw4kA9FO0RwoSqAF4DRgARV/2ExaPox2844Ub12iJR8rgG3laWZV2oDfJ4crft2aFE/g3W80huMtf1e9gbBlA96OMcrn/Xq8dWSQ7jf1MC8njFZGN1CivKGkXiuAR6G4p4bX1vJhTJV2tGhALQl/4V047kFjOJTyFM0qD9VpQl2zjbfg7dCW1vDaNSjuBR5GddCtldtIExkDRQm0QL3LRNG6AqzHT3MbpFM6pzNKrJmPOlvuRRzDA51A471zI52iuS8+Y17HGL4L2E2sMaHNsP73dWeQk3wHkhJogXqX6zcAQ4GWD2/Fy1lOTx1IKgbHUnE3gAbU+RAYrRPHBmIl/HUYAR4JdIrmEmAhNcyYF14TP06sMayHQB/DQR4AB1OUpEbqVeds43r8gyHa8uEg37xXUv2hOjWyI4GjqXcGOloJxxFgR6cOJKFop2iuwGeJ76rp9W8h1kNVtI2PdYh2iubAUQI98Kw121pXAn0XvgQVqQ/poN+8R4At5SbQCEaB54GNNb3+BuAogyNrev2qbQT2oDEcaQVlJbABCqvp9aOVRKwHloAtTB1IQtHG8MBRAi1DwH7UNFtTLh9Gu3lHq3mtWsQOHFugeLqm1x/DewkHGcPFs/iRyIOcQEc7hr7u99ytwLJyc2UE6/EjrIN0v6mWwVH4r0hjeOAogZZRYAcUO2q8xnpiJaTR4q1atOXvWpOPAp4CbidWQjqwp7mZl4ydRqzl77oT/lvxTZVBEtJiJ7CdWO+5Kq3EV5E0A50xJdDSxDLRGLFKOMaAUwdx+dCgwD9kI21eqeMI76kijuGBTKCB5fjsZZDkw+bjZXS1xVvADuBB4o3hSPFWaQQ/yOrZ1IHIzJRASxMnHd2Kd+KYV/N1qrIB/wBenjqQBE4GDkYz0FNF6wpQzkBbkTqQBFYC2wtPGCNYjt8b654xj7ayNrCrKOgI7xCUQEsTyccGvPbutJqvU5FiJ/AAg3nzHmWiZCEAOw44jGbGcLQEejFwYupAEohWwz8KPATFAzVfZwOxZnSjPbRWSUd4B6AEeqDZAXj/0bo/bO4DdhHr5j2oS+AjwMbCN81FsAJf5ryt5uuMAYeat8uLYAvwNIM7hoOUbwDN9fsdJ0wNNOA/k5PBDkodSALqAR2AEujBNoyPgVrfqGUnjmizCdFma6oSrffoKLAJij01X2cj3tYvyBgunsf7eGsM56+JGn7we/Bp5SbLCFqfSwP1EGi+arSYWKsoA0kJ9GBbAdwHxSMNXCviaW6R4q1KtCO8G0k+Cj9YZguxxsTAjWHz+v2TiTWGm5pt3IB/5g81cK0KFLuAuxmwMYz//z5DmDK6waUEerA12Ss12o7qgevEYRMfrlr+nl60JfBo77kqDOGdZIKMYVuA7w2pfQwX8BDejSPSGB64h0D872dzoIOsBpYS6MHWZPKxHlhusKCh6/Wq1YkjyGxNJZYCBxBr+bvJetdoB+yUZUgD1YljFLi78D0XESzHD7JqagxHK6UbxIfAaIcADSwl0IOtyTfqevyD4vSGrtej4iHgfmJ92PRqFNiNH78egDVdKxhtNmw9sAgvaRgU0TYQjgIPNtCBoyXaGI7W/aYKTdXES4+UQA8sW4RvVmjkjVp4W7iHiTWbEO3DpletDhyWOpA2jdBsreA4cIjBMQ1dr1db8ZaEgzSGh/G/pyia3nMQrQxpA3AC2OLUgTShPMgqWhndwFICPbhaN9Em36jRGvkPWgIdcQPhZiiaqhXciLfMCzIjVryAxxzpobVX0Za/m55tXA8sDdSJYwx/oI+U9PdCHTgCUQI9uEaBe6HY3eA1o9WzDVov6Gizd432Si1gDz7brTGcIYMDgZPQGJ7NOP65H+SU1eJxBqsTxzC+qrY1dSAyNyXQg2uY5mcbo/VWbjXyPzB1IHUL2oEjRa1gtFWJaPH2YphYHTgW4ht3m3wI3Il34og0JsYYnBnoUdSBIwwl0IMrVfJxWrBOHPswGDfvJfiybpASDitIc1pXtC4GG4DhAenEsQK4t4BHUwfSpmV4p5+m78MRx3CkhL8X0UqQBpoS6MGVYrf6RmJ14ngUP4Z8EG7ew8AThOnAwfH4oRlNj+Fos2FjeGnDIHTiiFaC1OrA8WDD1404hiOtXPYi2hgeaEqgB5IdhH+gNpp8FN4WbiexEtJoZSfdWgGMFfBC6kDa1Dqta0vD190AHGqewEcwSJ04Rok1e5diBQXizUCPASeCHZw6kDqVHThWEKsP/0BTAj2YluNv1hRPutHaKI0zGAn0MGFqRwGv197SYAeOltvwThxBxkTxArCJWO+5bkXrnztCmnvwGHBquekygnH88yrIxseunQgcRKwxPNCUQA+mEeB+KB5JcO2NxPowj7bc2a1h/O8miiQt9yZ14oh0QmXf15Can6CpDhztaXXiWJbg2l0odgH3EuahtWutvvbqwBGEEujBlPK0rmgf5mPAqeWO+b5UduAYId7yt8ZwewahhnQEH8dBxrAtwDtwND6Gy04c24k1hqO957qxAj/ISh04glACPZhSLnWOAcsM5iW6fqeC9U3tyin4DF6kpcOUCXS0hLSs4+/rThwjwH3BOnDsR9oxHGllLdp7rhupViSkS0qgB1PKzTbjePK8NNH1O1TsxI8hj7TpplPD+EazbakDaY8dDxxKujG8kXglHK2Nw/0qWveCUWAnFA8kun60UroN9Pc9GOKN4YGnBHrg2EJ8xjHJG7WAe/BZokjLcdESpk6NAuMFPJ86kDYN47XITXfgaNkAHGlwdKLrd+oOvLayn8dwtP65Q6TdcxDthMpxvBNHlI2PHSk7cGgGOhgl0INnCP97T7lcH62zRbTlzk6lLIfoxgjegWNPoutvwusUg4zh4nlgMxrDOUkd7xiwNNChVmP096FWrb72moEORAn04BkBHoDioYQxRJvR7ff6u2hLh0lb7hU+m3sHGsNZMD9B8xRidZFJ/Z4bx09BjHKo1SP4OQL9mkAPMdHhR4JQAj14Us98QLwd1ePAUrD5qQOpWtClw1T9cydbT7wx3M/Jxz6EGcO2H76JMOVD4HbiHWoVbeWyE6PA5rJNpgShBHrwJOmfO8UYsNx8BiSCDfjGxyB9UztyEvGa9+cyhiN9mEd7aO3ECLCjgJSrap04HZhP+jEc7aGqn0+FTb0iIV1QAj14cnijjgP7A6cmjqNNxXbgYfrz5j3MRElCAHYUcDjpx3C0LgbjwCFgx6YOpAZDpB8PnRgGdkFxb+I4Io7hfrwHQ/pNpdIFJdAD5cXm/ak/bLYBjxFrRqxfb94jwKZAzftH8G4hmxLHsQE4xuCIxHG06zZCbXzsiEqQuhNtVaKfD7VaQawxLCiBHjTL8LKJpMv1BRj+tB3pw7xfO3Ek3ZDXhRHgDiieThzHRuAFwoyJ4lm87V+QeDuSS0Larlzec+PEOtRqjFAbH9tjcCR5rKpJh5RAD5YR4OGEzfsni9aJI1q9YLuGSD+b24ksko9i4uAZjeGEbOJQpkjL37kk/K1DrZakDqQ9xYP4xsdIEy/tGMIfxjenDkQ6owR6sORQ/9wScflwebmDvp9EWzrMJfmAeJ04om18bMfpeBIYZAzbvmRSs13A3cAuYo3hvnsIxH/+W8uHcglECfRgySmB3oh34ihSB9KmcXznfJDZmrmVJ+kdQqzZu9zGcKQP835MPoaAXYX3CI7gFHwDdfJVlNImYo2Jfiyl0wbCoJRAD5acZu/G8fZpJ6YOpE134xsf+2kGbxjfkBdk6dAOA44ln+RjHFieOogOjAFHgx2eOpAK5XRPa8cI8ARe/pODaAlpP27mjjaGpaQEemCkb94/xe3A04S5GRYRNz7OZRS4vTxZL4JhfANqLh82G4CTDBalDqRNmwi18bEtOfQE78QwsAmKF1IHUopWSjcOnN5npXTRushISQn04FiClyBksVRU+MznFmJtwtpIrBnHuURbOhwC7oXi8dSBlDbhJUhBDtgpngTuItZ7bi4Rx3BO8W7EO3FEyQU20keldAYHA8eR15iQNkV500jvhvGlw7tSBzJJtE1N/VZDmkVHiw5kFW/hG7DuI9aY6JsxXCZ9y8lnRaIdI2Q0hvFYDsBrsyPot1K61nsx0hiWkhLowTEKjGe0dAjxEujyKFmLsvFxLqPk9WE+lxyX66MtgUeLdzZLgIXkNyZmYAWZPQQCd+ITK0Huw0WrhKtfxvAocE8Bu1MHIp3rNIE+HPhB4HPATcB24Fn8qfBa4I+BNwMHVhijVCO3pUOI1wt6Iz62T0gdSK8MFgPHkN+YmI3GcO+ixTubIaBVlhLBcfiSfTZjuPCa+NuINSb6qZRuORmNB+lMuwn0MPCX+HLlPwA/Cqwqv/8ePKk4B/gQ8OXy934T/4CWPOTU/qtlA3B42U4tgq34hrsgszWzGsU35OU0GzYLOwA4mfziHSPWbNg4cBJYP0xyjAIbyyQwgmFgD773IyfRVgKjxTubHFfVpE1zJdAHAJ/E/4LfCnwW+GH85Kd5wFH4Mtph5X+vBN4PfBX4CJ5wfKiN60itrCCT5v1TbAaeI8zNsHgOn63phxrSYeCuwpdvI1iGb9jL7cNmHFhisCB1IG0ax3+OkWYcZ5JbOcRchoEtUOxJHcgU0erix4GhPimly3FiS9o0V2K7Dp9Z/l58+eln8Bnorfjs1WTP4ydzfRb4EXxW8aeBD5bfI+mchPdczuqNWnj5zx3E+jCP9mEzkxzLIWYzDGyH4pHUgUyxEdgXOC11IO0pHsFL7zSGm5drsrSJWCUR4/iq90mpA+mF+YE6JxNrDMskcyXQHwJeAXwFX3rqxBPA3+BLFJ/pPDSp0HK89CC3pUOItxzXLwl0tN6jWSYfBTwA7CReGUfoMVyeYJpbR4u5ZDmG8fvAoeaHFEWwFT9DIPQYxj+X9yXWfVgmmSuB/gYvnWme38brHjzp358DLukkKKncCN68/7nUgUwj2klY0WpeZxLt9KucO4ZEO9K7HzpxHI8fYJPrmJhOrvWuW/DP6SBjuHgenzWPNPEynWHgwQIeSh2IdKeb2uQrgdNn+fNXArd0F47UJOdawU3EKuEYAw4DOyp1IN0yb/11Ev6zjyLX2TuI19ki/Aw0/vPeg89GBmBHAEeS4RguS+m2ojHctGglSDJFNwn0qcBa4H34MlrLQuD3gUtpb5ZampNz8rEBON7g0NSBtGkzXu8f+eY9hL/3c5wNm4bNx2uMcx3D0VYlxoElYFE2Pk5nBLit6Ly0MJUhvFtIrg+tKqVrXrQSJJmimwR6BfBtfLPgF/AOHGcDNwL/Ffg/xPowGQQ5J9Ab8TKhILMfxTP4xsfIN+8RYEcBD6cOpE2trj+5juEx/Djk/VIH0qZxPNYgGx+nlfM9bTojwDYonkodyAyiJaT9cJhKtDEsU3STQG/HW9q9D3g9/uFxDX7IytuBdwG57ZQHOBf4Ezz5HwcexZvwbwEuAz6Bf1D3GTsCOIJM36gFPI73DQ+SQAPxPmymita8fxh4BIr7Uwcyg014G7tTEsfRpmI7fo+ONOM41XLync2dTs5ldBCzDOkQsChnCOzFJjr3RLoPyxTd9mc24HK8J+7R+GD4OvDNiuKqUgH8KZ7kfxB4FX6S3C589/zRwEXAL+DJ9KeJM5PUjmF8g8jm1IHMItryYbR4p8p5Q950hsk7WdqGPwhGGhPRHwJz3ZA3k9xnG8eA44KV0gXa+PgSS/E2dpHGsEzRTQJd4P2gb8YPUflB4OP4AStrgfMqi64aHyh/rQN+AL9BLML7L55U/vtifPb8cuDHgY8mibQew8AdUDybOpBZRPswjxbvVLl/mE+VdbyFTyhsRAl0I8xXO48i4zExjazHMPFK6VobH0OOYTzu3cC9qQOR7nWTQP8H8GfAFfjJg/8E/DZwfvnnV5T/nYt34zWrF+I1249O8zW7gX8DXgNcD7ytsejql/uNG2IuHx4Ptjh1IJ0q63SXEmvpMMIYjngYRdTkYzme7OW8qjaJLQJOJOMxXJbS3YvGcFOWA5uLl7YJlkC6SaAvBD4MvIm9n56uB84EPgf8Yu+hVWYYT+rbObJ4D14jvbLWiJoVIfkYA04t26tFEGy2Zi+n4V1ygiwd2j7keQz9VNF6K48Dy8H2TR1IF0aAbYGOoV+Or9xqDFcrcgIdrQRJptFNAn0mvhnvhWn+7HHg/cA7egmqYmPABcABbXztPOBi4NY6A2pYhF6T4/hYXJY6kPYUu4H7iTVb0zIM7CriLB2eiL93cx/DG4Eh27u1Z842MnGUcDQRJgUmG8KPoZ9u9TMn0RLSjcS8B0O8MSzT6CaBbucv/V+7eN26fB6v1b4SeCdwyDRfswgv2/gW3q3jS41FVys7EK/zzjr5KHwz50PEmtGN9mHTMkSYpW/Af8ZP42VYOduI30eOSx1Im+4EniLuGM769Z81CwAAIABJREFUnjZFlGQpYhnSCWWJTDTRusjINOZKoD+Mz1L04jB8VjqVT+Ez5qvwGuhH8JrnbfiHyK7yv7+Id+P4HPAHKQKtwRAxlg4h5mEUEZOPaM37h4HN5fG9OduCn+gWZCNh0TrUQ2O4flES6A3AkmCldBBr4gXzLmAHoxKO8OZKoN+Cz1b9LHBsh689BPwWPnOU8qnW8AeBpcDvAd/Bl98Pw/sjP4jPTv8OXkLwPuKcbjWXYeA+KHalDqQN0RLSaPG2RPkwbxkmwGxjeSLe7cT6MA83hg0OxMtOoo3hCPFGLKW7l2BjGI/3GcIcQy8zmavf8RuAn8IPGWkln5cB1+I9oB/G654X48noKF4C8Rq8VvpuvM3d16oPvWN3Ah9LHUTDImy+atkIvDd1EB2YdBxy8UzqYNpR1ucuI9bS4TBeWhVBxG4yb0odRIeW4eM4+4cqZ/PwyZvs78MF7LSJUrqbU8fTpnAPgfjPd0sBua+qyRzmSqAN+HPg74GfwGdyX93G696KJ0N/jy9rShpRZj7Al7OWG+xXeIP83LWOQz4dWJ84lnadiNfpRlo6HMYPQoqgtWE5inF8dTGS1jH0D6UOpE2nkfcx9FNFOyQqYgIdrQRJZtDuiXu7gN8Ffh9Yg89Mr8HLOg7ByyDuw2emv4HXA6q/YXqRZu/G8fZqSwixya3YDvYI/jOOkkAHWzq0o/BDM6IkH+Ok3e/RqXHgULBjMz4mfarcj8SeahjYBcV9qQNpU7SEdBx4Y+ogOjQMfDd1ENK7To+sfgG4qfzVjw4Fri7/vZel2COAT9L+BswlPVxrBjYPn/2Iknzcg/d1jdQlItqHTbSlw2E81ijjYRNwtMEhxfQHNuXmNny1ZwjfFxJBtO4FIWr4J1EpXf2WA/87dRDSu04T6H63L9VseHwOTwjbbRN4cAXXnOp54NeB62p47coVYOY3wxHitBGMlkBHWzpsHUP/dOpA2jSGTzKMAFcljqUNxbNgreOQL00dTZtG8bLCKCKV0YGP4YildKcRoDStPIb+GGLdh2UG3STQ+wIfB94KHDTL10XaTNOyk2oOQngU/xm16yfxspgKFS8Av1Hta9YuWkI6BrwndRAdiJQoQbCEv4AnzTdODxMigQYCnT5nXkscaVUN/Gf796mD6ECrlO5UfIUic8UOsJ3EOdlvCJ/cirSKIjPoJoH+GD6zKVK1TcD3pg6iA+PAMj8OOfs+xeA372izdzemDqJDEQ+jOC91EG1agifRQZIP2wcfC2EeAtm7lC5AAg3EmngZAu4q/HAoCa6bkwjfi28qfDV+Mytm+CXSqQ3ASKDjkMfxOvdTEscxJ4Mj8V+RPsyjLX9DoBndUqTkYwQ/hv7u1IG06WT8UJIwY7jwEqSNaAzXJdSqmsyumwT6ZOCv8aXgCDVSLefiJxJ+G3/DPQo8iXcMuQzvdb00WXQCPrMU6Tjku/DZmgg37yH8wzHIhjxbDBxPoOSjFHEG+hiwQ1MH0oZIG4zB7wtP4WcQRLKReGM4wj0YfAwHWUGRuXSTQD+MJ55RFHgf2WuADwKvwo/S3IXXPB+NH+H9C3gy/Wm0uTKVLfiJbkFuhkUrIY3QN3UI2Fb4B3oEQ3grzEgdDMA/zE+x9jvwpLYR/zlH2LMSbUWidQz9C6kD6VCkhBT2KqXLXrQxLLPoJoH+a+DN+EaDCD5Q/loH/ADeqm4RPpN+Uvnvi4G3A5cDPw58NEmkA648Dvk2YiSkLVE+bKJssmlpHUO/O3UgHRqjum4+DSieYGLjY+6iLX9HTZbGUCld5cpj6E8h1n1YZtFNAv0/gLXAvwHnAAdUGlH13g3cAVwIfIHp+7Puxv9/XgNcD7ytsehkqojHIUdIPiL2zw2XfBS+qtU6DjmK7MfwpGPoI61IhBzD+M84Uind3cBjZD6GmTiGPtJ9WGbRTQL9LN666014j+En8CXAqb9yMQxcgcc5lz14jfTKWiOS2WT/YT5FlHrBaB/my4mVLE0WcQznnvCfgCd1kcZEtIfWlmildIavXOY+hofxY+gfTh1Ie+xksDv8kBqZTje1vtFuCGPABfhM+Vy12/OAi4Fba45JZjaGl9FEMQYcAnZcrsf1miceJxFr6XAU+HrqILo0RqwypDHybx85gh9DvyV1IO2xo/FDMyK95wAvpbOJUrpLUsfTpg3k/56LVoJ0DrA40AmPjetmBnqozV+5+DzeP/RK4J3AIdN8zSK8bONbeLeOKCfh9aONwFEGh6UOpE1bmDgOOVfBlg5tAV4rGCTel4gwozvZRuBksJw3PkY7hr51YEaQhP8loo3hCN1voq2qqWPIHLpJoKP5FN6+bhVeA/0IXvO8DW8vtKv87y/i3Tg+B/xBikAF8DfsC8RZPnwWuJ284x0Gthc+9iM4HV8di1RyMtk4sMx8M2EEG/FYl6UOZBbRSpCiHUM/VbQypAjxRhzDkeJt3FwJdKue+Q+n+b25fuXCgA/jPZ5/D/gOcD8+w3kE8CA+O/07+AfI+/D6L0mg8Fr1beS/HDdZ7kv20ZYOR4CHoNieOpAubQAWEKavfLED3/yY8+EZo8Qbw5HinSr3e9pUG4BDwY5NHch0ymPolxFrTEQfw7Wbqwa6NX2/Y5rfi+ZO/BhyyV+E5bjJNuKlP7kaQkuHTboHeJxYB3/kvhl2Ob6aGMUQ3q0qqhdL6YJsetuKT3wN4RNkuYl4DH20rjeNmyuBnq4GKlJdlMQU8TjkH00dxCxG8Pr+KEIvHRZg5vGPEmc/RbZL4JOOoY+0IW8Y+PvUQfRgHC+lG8XPR8hcsQdsC/5zvzR1NNMYJd4x9AcQ+D7chEGogZZ4sv0wn8EG4Fiw7DY+2kQpQaQbYT8sHUYcw7k+tA7jG/KCzObbYrztXqSEfy/FxBHkkcZwzmUnw8SazR3Bu5bdkTqQnHWTQM8HPom/uVreh5d5bAN+kTgnGEmeNgEnWv6H9LRsxmvtc1wCX4KvNAW5edu++CbCIPHOKGIXg9PAumltWrfWMfRRNuQtx+8HQZbrZxRtDG8m33gjltHdFvAY+kZ1k0D/V+BnmUhuTgT+HN+Udyjw28B7K4lOBtV6fGzmejOcongcuIs8Z/BGgUcKyLJH9TROBRYSePauNAYMBzoOOeeNj9FWJEaBewIeQz9VzjO60xkjz3swxBzDkeJNopsE+l34h3HrRvs2fFb6fPywhh3AT1cSnQykwo9bf4B4y4c5xjtMrJmwEXwD3l2pA+nROHAQPsEQwT3kexxytOQjdA3/JOPESqDHgWNyK6WzickgjeE+000CfQrwXfxmC/BK4Gb8WO9HgW+Sdz9RiSHa8mGuNa/RmveXCX+RUyvMbmwFniXMGC4MHyc5juFoy9/9knxsBE4wfxCMoLXxMbcxfCK+Yh9pIiNawp9ENwn00/hSH3jz/YuAKyb9eW59oCWmaMuHuSbQo8Qqh+iL5KPw0yk3ozHcE4PFeAIS6cO8L8YwXkoHmY2JmRVP4itXucU7im/KvD11IO2xY/ETm/thDNeqmwR6E540LwbeChzFRIusA4DXsPcGQ5FuZPdhPodx4CSwA1MH0lIuHS4n1o1whFjxzibiGM4t3iF8QibImLCFeB1/kHhnVvgJvfeT35iYTY5jeBjYXPjseAQj+ARA1GPoG9NNAv0XeE/Ou4F/wcs2LgFW48uWxwJ/VVWAMrA2Aktt7l7luRhnImHNxcn4hrwgy99W0F9LhxHLkIbKv4dcDAH3lclcBMvw+0D4BLoUbQznuHIZsQRpCxTPpg4kd90k0H8L/DJeA70d+LHy34/B+wb+GrFOjJI8jeGbU4PU0xcPk9/Gx1H8aPQovTxPABbRP8lHzl0BppPjxscR4pUg7YTiwdSBVCTn/uDTyXEGOlpHi34pQapdNwm0Ab8JHI/PNv9r+fvfwHvO/k98+l+kawXci69u5HYznE1uN+8RYi0dDgPPEKZWcE5jwOHlSXoR3I7//HMbw5E+zPst+cjtnjaXMbyULqeNj8Moge5LOolQcraJWMuHuX3YROzAsQWKfnkAvw1/eAkyhovn8ZhzGsMRl7/7KfnYCJxqE40DcjeO917P4j1nvjJ/CBrDfUkJtOQs2vJhbvWC0TpwREuWZlUeh7wVWJE6lg5kM4bN6/eXoDGc0hi+DyVKKd2jeGlpFmOYiQ15UY6hPwRP+vtpDNdGCbTkLGIj/9PA5qUOpDwBL9pMQrR425HrATszyWkVZTn+GRVk+dv2xRPNvhnDhe/reIh8xkQ7chrDI8DthfeEj2AEL9NVAt0GJdCSs3FgmcUZp+P4xsccjkM+FjiYWDfCaCfOtSObGd025ZR8DAE7C0/gIliClzr0TQJdilZKl1MnjohldHdB8XjqQCKIkpjIYNrAxDJuBPcCu8jj5j2Kz3rcljqQ9tgR+Ga7fks+opUhjQNHln8fqUUrQeqXY+inWk+8MZzLQ2DEMdxvkxi1UQItOduGt2HLISFtQ9E68CGHeEeALQXsSR1Im4bxDXeRjrttxxhwnMGhqQNp0ybyOQ45WknPCLCxD46hnyqXe1q7xoGlYPNTB0LMMawEuk1KoCVbhddiafmwO9E2Mw0Dd0LxVOpAKrYJH8dBxnDxFN43XGO4c9HalbVrI3B6oEOtstj4WD40R9uQpwS6A0qgJXcRl8BzSD5G8KXXKPryxl34IVPbiDWGkz8E2sQhSlr+Tm8DXtt9WupA2lNsBx4h/SrKKL6aE2RM2CL8EKVIM+ZJKYGW3EXbhDUGLC935Kc0RKxyiH5NPiDmGE6d8C8B5hFm9s72oX9noO/BHwSjjeHUExnDwF2Fn9AcwRDevakfx3AtlEBL7saA4bItWwQbgP1JuPHR4AjgKGLdCKNttulEDh/mnchh09gIsBvfmBvBqcAB9OEYLibamkUawzmsXEabFFgB3APFrtSBRKEEWnK3AVgEnJQ6kDbdhc/WpPywGeH/t3fn8XbV5b3HPyshIcwQZgghIRAyjyTMggyC2la0iq21VutUq7ZXe9tqa62ttXXqtQ7Xe6u3eh2udWjVOqCigEACJJB5JCEkQIAAYUiYE5Ln/vGsA0k4Z589rLV+69n7+369zitwztlrPZC1137W7/f8nh/sIs7o3UjgBGJ92LSiDh/mrVgNHJf/vaQyBViVJ28RTMIXPG9KHEdZVhIrga7DQ+tkYpXRdesMSmmUQEvdbQSeJcz0YVaH0ZoJwF0ZPJMwhlZMxBOlbq29WwucZHBw6kCatBav3Ux9DUe6HiYBt0O2O3UgJVlL+priVqwGxife1EodOLqcEmiptcy3QV1PvJt3yngnE+tGOAm4C7InUwdSktV4CVKUh8Cn8JHUlKPm0T7MJ9KF5Rt7WA1MCLSp1Wq8hj7JwkfzTaxOJNY1PIlYCX9yUd4M0tvqUJPZitSdOCJOHXZt8pHBY3gt75TUsbQg2UNg3i5tArGuiW5PPlbhNd5jUwfSnOxe/H2X6j7c93kVJIG2g4CTCRNvPSiBlggiTh9OzFfmp6Cpw/pZQ5gRaMATplTJx1hCbYltGd1fP3oX8DSxruE1pBt4mQjck/nOlBGcjueDkR5ak1MCLRGsAqYE6sSxGh+tObnqExv0LciLNAI9mTDJUtuiLSRMOYsyBe/AEWVL7NF4fXvXJtCZ18SvQbMozYrWVWgSsAWyR1MHEokSaIlgFXAYnhhGsBEfrUmRgEzCP+yidOA4DBhFFycfuTXEm0U5EezwBOeeBKwJ1oHjaeDO1IGULPXajlal7MQxiVgJdLfPoJRCCbREcAewkzA372w3volJipv3RLwDR5Tm/X0dOCJ92LRjDTDG4IDUgTRpDf73kuIajraF9yRgHWS7UgdSsogbAp0OlmIL8onEu4aVQLdICbTUXubJ8zpiTR+mqnmdQqzyjQnAvZBFqRVs1wpgKHEeAp8A7iHdNRzpgSpawt+uVcDkQKV0a/Ba+koXPprPlo4m1jUcLeGvBSXQEkXKRU3tSFVDGm0kIVq8bcngUWALseqgK58CN3/IiNiBo+uvYfy/8WDgpNSBNOlufBFf1ffhvofOINeE7Q+MI0y89aEEWqKIVkO6EpiWoBNHtA4c0UbMOxGthnQl1c/6jAFGEGY0zDLijZi3605gB2HKODLD33NVX8MTgPsyXwgbwUS8deSK1IFEowRaolgJTA00fdjXiaOyLcjzDhwnEish7ZXpb0iTkHYiRfvIyfio4caKz9uuE/BNMyI9tLYl39RqLTA1dSwtSFFKF21QYCLwEGRbUwcSjRJoiSJaJ44N+FbaVU7ZT+SFdlMB9Fzz/tXEK+E4CeyQCs85GVgbqAPHZHxUdn3qQCoSrR1jik4c0Up6ou1cWxtKoCWKaJ04nsM7cVSdQN8dqAPHBPwe1Cs379XE6sTRV5ZQZQIygTAPgIC/v9dBtjN1IBVJ2RquHauBCRWX0kUro4vWcq82lEBLCHknjtuJNQVe9cLHaFOHk4H7eqh5/0r8nhskAcm2A5up9iEw2gYU0eLt1ApideLoK6UbU8XJDA7JzxXtPtxL13BhlEBLJNFGP6qe7ow4ddgzN+68E8f9xJoCr+waNv88irahQ09dw/h/68Ek2GW1TZvwGbmqPjf6zhPkmrAReAeOIPHWixJoiSTl9sLtWAFMqnD6MNrU4VR6b+V3tE4cK6hu1udkfLQwyDX8fAeOXrqGN+JrO4Lch7PdeHJY1cLHScC9wTpwDCXWiHltKIGWSFYCUwJNH67CE4LSG/kbHIFviR3pRtgr7b/2FK0TR5WzKFOBJ4mzJfZofDS2Z67hDHaRpjVcJ6q8hqM9UE0GHoDs4dSBRKQEWiJZg3fiOD51IE3ahCcEVdy8J+CdC4K0hLPD8IS/Z5KPXLR+5quAE8FGVnCuSXgHjt0VnKsIk4GniZPwFyViKV1VCX+0MjotIOyAEmiJZB2+mDDIzTvbjSe0VSTQk4G7Mu+hG0GwWsHCrALGms9MRLAaT2irSPonE2sGZTKwFrJdqQOpWLRWdivxThz7VXCuaDXx0eKtFSXQEkbQThwrqebDJlo98WTgLsiiJPxFidaJ4yl8JkXX8Iv1avKxEphkcfKHVcD+wKllnsTgcHyb82jXcKQR81qJ8gYQ6RNxCryqEehIN8KeTD4yeAzvxBEkgQYquIbNFzKdTqxruBdr+KHi1nCdy+7B33dl34cnEauM7iB8fU4vXsOFUAIt0awk1layq/Dpw6Elnyfa6F20ntVF0kLCFxsPjCDMNWx9Lfd6MfnYCDxOrPtwFQsfpwAbA5XRnY7ngL16H+6YEmiJJmInjhHAKWWdwOAY/CvSjbAnR6BzVW+w06kqFmFNBh7LfOOWCMbgo7A9dw1nL4yyRquDLjveaPe0ycCWHtrIqnBKoCWaVfhuTyelDqRJd+MjEmUmTJOAvq3DA7Aj8E4qQfr9Fi5aP/PVwLEld+KI2L2grz68F0XrZ15F55Bo13C0TYtqRwm0RLMe/+CaljqQ5mSGJ0xljn5MA9ZnvsFBBJPxzg69evNeDowxODR1IE1ai/99lXkNRytBmgysyTvt9KIVxCrhWAWMBxte4jmiXcPTiDViXjtKoCWUvEfsGmLVkK6m3OQjYvuvjXmHh17U93cVZAo8ewrvdVxmvNEW5EWbri/aKmCCQRWt4YqwChiG19oXzuBo4FhiXRNTifW5UTtKoCWiaI38y17AEq0DR7RkqVD5IqO7CZNAAyUuJDQYDowj1jXc6wn0arw1XGlrO4qV3Q88THnvuYn4Lo1RyugOI17LvdpRAi0RrSBMCQcAy4CJZUwf5ospo00dRtvutgzLiTUFXmb3m4n46GCQa8L2wx/gg8RbvMwfAB8j1n24zLKTqcCGzHemjKDv/0PPXsNFUAItEa0CTg84fXhaCcc+Ea+ljTQaFm2xTRmi7eZW5qzPJGBrBg+WdPyijcM760R6z5Whyi2yi1DmNRxtQV6vbmRVKCXQEtFK/AOs1J2lipNtBh6lnNGPKcCzwB0lHLsE1tdyT8lHrORjBXA02HElHDtaT/Cp+OjrPakDSSxiP/Oy4p1CrK5C0d5ztaQEWsIJPH1Yxs17GrA68zZ2EUzFt2QPsltXaZYDx5ovPIrgdmAH5T0ERppKzpOPzFIHkthyYHrqIFqwEjgF7IASjj0NL9WLItp7rpaUQEtU0RYSljVaE60cYiqwDrJnUweS2O34Q0+QMo5sBx5zGQl03hIujGhrDsqyChhrUEZCWoaV+JbxhX5umJfRHUGsWTWNQBdACbREFa0PaVmjNdE+zKfh/y96WuZlN7cTaxZlOQXHm/fCPoV413CkeMuyHE9IozwEPoLvdFn0e24a/n4OMqtmJwBHoftwx5RAS1SriJVArwBOBits8wx7YTQl0khCtIS/TNFqSMt4aO0bDQxyTdhBxEv4S5H5uo57iXUfLqP7zVRgbaAyuil4OVaQlnv1pQRaoloFnGK+mDCCviS3yOnDsfh/f5AE2vpGq3o++citIlYZ0gpgUt7GrShTgM0ZbCvwmGWahLeO1DXsonWTKaMFarSNrKYCt0O2M3Ug0fVCAn0eL37izIDXAtcAW4CH8n/+vfxnUn/Rpg+344sfixz9mA48msFdBR6zTOPweklNHbrlwFSLcx9eQfHdb6YTa/HVVGAzZFES/rJFW0hYxixKtGtY9c8FiXLj7sSNwGf2+d6/Ad8DLsJr8Hbk//xN4N9REl17GWzF+8YGSaCB4jtxRNvRT+2/9rYaOBgYnTqQ5mT34H9/RSYgEUfvNPr8gogj0MeAFdL9Jt+L4HRi3YfVgaMgvZBA7+u1wFvwBOxV+AfYiXhv2m8Drwf+IFl00orlxKohLXo3t6nEGs2dCqxS+6/nbQCeIlYNadHXcLQFeUo+9rYSON7gyNSBNGkNxbZjPA2flQlyTdgQ4q2bqa1eTKBfnf/5XuBHwO783x8C3gxsRAl0FEvp7enD6cRLoIN80JQv83tPmVtkl6Gwazhv/3Uksaa/p6LkY08rgV2E6SaT7QTWUVy804CHszizauOAA9F9uBC9mEBPyP+8up+fPQvMI1ZS1ssiJh9HFTF9aBCxG4AS6BeLuCNhUclHXzeAdQUdr2R2DL7xja7hXAbP4DMpka7hIjtxTCZWH/4pwHbirJuptV5MoDflfw7UTmw4Xucn9bcUnz48JnUgTVqH78JXxIdN3wdAkBFoOxBP+DV6t7elwMzUQbRgFTAmb+fWqen4LppRugFMxluVqf3X3srobFGmIuu2oy0gzBN+ldEVoVcS6BOBlwNjgM/l33vLAL/3CmBRNWFJh/rq2YLcvLNn8SS6qAR6UwZPFHCsKkzCu6Zo9G5vy4HTzKdVI1iBf24UkYBEm5GYCqyH7JnUgdRMxE2tJudtNTsV7RqOtuV4rfVKAj0BuAqvb/4ZXrP1YWDWHr/zVmABcAjepUNqLvPkeT2xbt7LKCbhj9aKqK/916OpA6mZFfiDxcTUgTQnexTfza2I91y0a1gLCPu3EphocbpXrcAfWMd1chDzBgQnE6sDR7QR81orsiF+XR2O9y09Lf9zz38eucfvfRBf0PJe4OcVxyjti7aQcDnwugKOMx24voDjVEXJRz8yeNh8AdIM4sx8ddyO0WAYPisR6cN8KvCT1EHU0BJe2JJ9Q+JYmnE3vnHPVDqrv5+CPzREKaM7CM97gsRbf72QQG/DP5j6+3Da84n5T4Fb8d7CEscK4MrUQbRgGfAR380t62Tr1ynAF4sJqRJT8Q9aebFoU+ArgDM6PMZp+HqTICPQNgT/O/pY6khqaBPwOD6zFiCBzgysbwH6f3ZwoCnA3Zn/t0cQLOGvv14p4RjInoX0P0XJc0TLgMn5iFYEy/C+oePbPYDBKHy2ZGlRQVUgWq1glZbgI9BRLKPzWZ/pwEMZ3FtAPFU4Fe98E2nEvBJ5O8blxLqGi+jEMZ1YgwLTgI2QRUn4a6/XE2iJbxmwPx0kpNXK7scf1DpJQKbhG3AEGO2BvG3fccRK+Ku0nDALYQH/ezwSbFQHx4jWT3k68Ei+G6O8WBEPVVUqogXqDGI9UKn+uWBKoPd2BLA2/5IAMrgf31UyWh/STj5spgBrMl8MG8EMvF+s3lf9WwEckc8sRLAOeJrOr+FIMxJKPhorsrdyFZYC48AObufF+YLJaDvBRtt4q/aUQO9tKL6v/empA5GWLCPW9GGnU/bTiTWaOxNY2WHNdzdbh88oBLmGs+fwEbxO+lfPINb090xiveeqthQYa75oP4JleAlnuw+BY4HDCHMNW0a8kpPaUwK9t634k2WUdjziOu4KULFOpzunEK91kpKPAeQzCauJdQ233f3GvPvRScQaDdMIdGMr8YQ0yDWcPYm3QG33oXUavqPfpqIiKtlYvEVvpPdc7SmBlm4QrYZ0GXBcO1t6m9d7TyTWh7mSj8FFqyHtJN4Z+O6DQR4C7Sh8ky09BA4ggyeBO4h1H+6kBeoMYFm2dyOCOpuOdyTblDiOrtILbexSGI1v2DKiyd/v21ZcI9/tWQqMNjgyg4dTB9OEtXhN8DTgly2+dgr+vl1cdFDlsAPxBZ5KoBtbCrw7dRAtWMbzNaRZq7thzgRWZfBsCXGVYRq+adOa1IHUXLRt6ZcBr2nztTOJVQ4xHVihLbyL1UsJ9JnAG/EtaI/Pv4YD9+VfNwFfppjOBlvwnQ6bHeG/FHg7cZ5m62YN8ByeXAbYXCR7DmwNflNrNYGeDtyTwWPFx1WKKfj7QAl0YyuAUw1GZP5wVXfL8Af+afi9sxXTiDWVPANYC9mO1IHU3HLgN1MH0YKlwIfb7Mk/DfhxCTGVRVt4S1sy4At4ctr39ThwF74j0eP7/OxLVP9g8c783G2tCBYwWG7wvtRxNM++CvaNll8FXzD4QRkRlcPeARak3V46BocZ7LbONyipkG0Ae1fLr/L36n8rI6Jy2NfAvp46irozeKXJvCwOAAAgAElEQVTBU+aL8QOw48AMbHJLr/KOOWYwu6zIimcb/F6c3HA81zkndSBF6IUa6HfnX0uB1+Ot6g7B97Afnf/zYcCrgRvxkeBAN3fJLSHW9GG7reyidQPQAsImZF6feDex6qBbriE1L2ubSKxrQjX8zVkGHIDvMhlAtgWfLW51IeF0vIY/SB9zOxRfRKhruGC9kEC/CdgInAd8l/6nvrcDPwQuxrfzvqKy6KQo0RZhLQdOBxve7AvM369TiXUjnEasfr8pRVsMu4LW452Ez/AFKeGwYcAEwsSbTgabgUeI1Q+6nffcNGB9oBr+yfiob5BFu3H0QgI9EZiHrxIezE7gWmLdAMQtBiaZj4BEsByfzprYwmtOw2dMoiwg7KuRjTTamNIiQk0LswyYCtbK58hMYFPmiVYEE/DON5EeWlNaSph+5kB78c4kzD0YgFnA+jYW+8ogeiGBXg2cCxzYxO8OAy5EI2YRLcNr71qqZ0snewgfsWnl5j0deBSIsp3wWLzDjJKP5iwFplmc+/JSfN3GuBZeE20B4XTgfsgeTB1IENFmUdrZhGsqsXKEaB1Dwohyo+7E14FTgPnA6+h/p6RD8LKNa/BuHT+qLDopROaJ5Sb8aTuKJbQW7yxgSaDeozNQ79FWLMLvReNTB9KkTfjfbyulU9E+zGejB8BWLCbWLMpS4BiwE5r5ZfNZw6nEuoZnESveMHohgf4i8Hn8qfi7eKK1He/CsQn/ANiOdzY4H/gq8JkUgUrHok0fLqa1BDraYqY8XvUebUZeQ7qVMLX8meEjjk295+yFtneRruFZxJquT20JcLzBcakDadLtwNM0/54bjyfRQWZRbDg+K6truAS9kEAb8Cf4NOOngV8D9+PbyR4FPISPTn8Cf3P8IV4LLfFETKBntFBDOoNY9cQziJUs1UGrsxKptfIQOBbveBTkmrAhxKt3TW0tnpAGuYazXXg5RrOfG7OA+zJ4oLyYCjUFT/g1Al2CXtpIZRPw56mDkFItAj5gsF/mG6vU3WK8hnQ8/sEzIIMT8FGdSB/mswjVs7oWFgFzUgfRgsXA7zT5uzPxGcCN5YVTqGCLdtPL4Dnzh/wzgKtSx9OkVgZezgBuKzGWos0E7oYswg694fTCCLT0jmh9SDcDD9Jc/+rpeNuk20sNqTB2LDAKTwileS33Vk5sEXAs2IlN/O50YEWgGv5ZeLeQTYnjiCbaTOBymh8xj1bDrxKkEvV6Ar2WODdzGUReQ/oAsRaxNNu6bBawPItTXnQG3jpydepAglkEHGW+0VMEa/DdXJsZNZ+D99mPYjawSDX8LYvWjvFWYBzYyEa/lO+wOJNYI9CzUQJdml5PoKX7tLvDXyrN1pBOJ8zCFeD57gXZrtSBBHMnnpAGGcHLduMzP83Momj0rjcsAUYbHJk6kCYtx0v+BruGxwMHEeYatv3wRbu6hkuiBFq6TbQtvfN4LRvk96KNJOSjd9KKDHbjU+BBFmEBTbQuMy/nOZYw14RlaAFhu1YBOwhzDWfP4DEPNmo+C3ggg3vLj6kQE/CSRl3DJVECLd1mMTA7b5kVwWK8N/kpA/1CPpJzCrCwqqAKoAS6fYsIk3wAzcU7B28ZGqSGn1Pw96WSjxblW1yvwMu4omim7CTaAsJZwBbI7k8dSLdSAi3dZin+wRelhnQTvlCpUQIyCx/RCbL7lR0HnIgS6HZFW4S1GDge7PgGvzMLWBZoAeFs4DFgQ+pAgop2DS9i8IQ/WglStHjD6fUEegJxRiqlOevwD765qQNpTmYMXnYyF19A+Gw1MXVsFr6AcE3qQIJaBIwKtBnFGvzvu9EI3hxijd7NBJZqAWHbmklI6+Q2YOxACwnNc6XZxLqGVcNfsl5PoKXL5CNcS4i1CnywhYTRboSzgeVaQNi22/GHpSCLYbNdDN4KLFo98Sw0g9KJJcBY841zIliOdzgaaCDjVLxnf5ARXRuCzwAEiTcmJdDSjRYQZgQa8OnORiPQc4hV/6ypww7krQqXAWemjqUFA86i5C35jiFWC7t8BFratAzYRZiBjOxZfCZloIfAM4AtGdxdXUwdGQcciu7DpVICLd1oMTDL4lzfi4FjwEbt+wODo4GTiDV6pwWEnYvWS7fRQsJZeGu+O6oLpxN2Ev6+i/Seq5XMt/O+nVh10Lcx8HtuJmG2oAd8AGkrcXb9DClKgiHSilvxp+/xqQNp0jq8Q0F/o+ZnAM8AKyuNqG12NDAaJR+dupV4W3qPBjuqn5/NBpbkLfoimI3XdEfpGFJX0eqgGz20RhsUyDctUg1/mZRAS9fJvLPFA4Qp48h246Mf/cU7F08+ouxAOAtP+FelDiS4W4DjzR9GIliF/733lzBFW0A4B1isGv6O3QKclTqIFizCdyQ8fM9vBl1AGG3Xz5CUQEu3GnRzh5pZSP8JdLTFV30LCKMk/HW1Di97CDKCl+3E6y37q9uOdg2fSaw1B3V1GzDGoL9ZiTrq25Fw31KkvnriIDXxNgx/z+kaLpkSaOlWAyWkdbUQmJ2vnt7TXGKNJMwh1lRnLWW+AOs2Yi0kXMA+8RqMweuJg1zDNgR/aFHy0bml+KxEkFHorK9Ubt+Bl7PwHQij1BNPwXcgDPKei0sJtHSrRcB0g/1SB9Kkhfgox8S+bxgcj39FSkjPxBMp6VyjRU11lD+07rUtfbAFhJyOt15TAt2hvOxsOWFmUQCfKdk33mjlEHOATZA9mDqQbqcEWrrVLfhT+NTUgTQnuw9vkbTnaM0cPPlYnSSkltkYPOG/OXEg3WIhMDdQN5mb8W3nT9vje2cDtwRaQHgWcB9km1IH0iVuJtYsys3AOft872xi3dPOxD//pGRRbswiLcngIWAz8Ubw9uy80Lf9cZTk40zgYWB96kC6xG3AIYTpJpNtwhfv7pkwRSxBihRv3d0GnGFxdvxdAIwCOxHAYAQwjVjXRLT3XFhKoKWbRauDvpW94z2bWFPJZ6LWSYXZo5tMpBG8heTxGgzHr+ebkkbUmvwaloIswBcRnpo6kCatxluK9s0EzsLLAIOM6NrBeBmgruEKKIGWbhZtM4qFwFSwA/MRmznEqic+i1jxRhCtDnrPXUCn4CN4QR4Cra/kK0i8IWwAHiFMHfTzLUX7HlrnALdnXkoXwSzAiLVuJiwl0NLNbgGmGhyYOpAm3YYnzjOBCcARxBn5GI7HrQS6WPvOStTdAmA62Ag8CdmQl1NFMBMYihLowmSezEW7hm/hhQQ62qLoucAqyJ5KHUgvUAIt3WwBnpAGmQLPngDW4DfBc4F7Ml9YGMF0YH9ifdhEsACYaf7/NoJb8SnvGfhirEjlG3OBdZBtSx1Il4m2ocoC4Ayw/Yi3gHAOegCsjBJo6VqZb8e7gjAJNJAvuiHeyMccYANkj6QOpMvcBgzDFzIFkG3DF5HOId5ipjOItdtcFLfhLUWHpQ6kSQuBA1/Lf5yH9zGPdA3PJdamRaEpgZZuNw8fzY2ibxFWtNE71T+XIIOt+KzEealjacGC0dz9ErydXaRreC4avStDX0vR6akDaU72ALBpJkteBzyB97IOwE7EE/4gZX/xKYGWbrcAODNQG6WFh7FtHF4DHSkhjTZiHskt+FRyFAtns+gc4Fl8BigAOxLvFKEEumD5Q+BG9m7RWXe3jGLzucCSfFfQCM4BthPmPRefEmjpdvPxrYRPTx1Ik5ady/yndzF0F2FWUttIfLRRCXQ5bgRekjqIFiyYyZITdjB8WQY7UgfTpPOAp4ElqQPpUjcA56cOogULxrHhVGKN5p4L3AJZlIQ/PCXQ0tXyXrpbCFMHnT13OT+/7y5OfjjzEbwI5uKxLk0dSJe6GTjWfHo2ghVnccvuFUzdkjqQFpwLLIQsSsIfTX87/NXWITy+cAorD7qb0atSx9KCaGV/4SmBll4Qqg76En7Fdbw0dRitOA+4TclHadYBDxKkDtrIdp7DTfZ9XpM6lFach98npBw3ACcbnJw6kGasYeLuQ9nOm/h6kHZwdhDe+WZ+6kh6iRJo6QULCNJGySA7lTuO/RWXHJ3fFCM4HyUfpcl76Uaqg554EE8O/S5XjkodSHPsAHyzGiUf5VmLb6gS4ho+kXvn3MNJz1zPBVNSx9KkOfg6H5XRVUgJtPSC+cAUg5GpA2nCxGHsPHge5xkhkv7nN8y4PnUkXW4ecWpIL3yKA++/g1OngUXYxOhMfAMVJdAlyR8CI13DL7mTU24HLkwdSJPOB5ZBFmXHxK6gBFp6wRJ8MVOE7WTPAjZvZtQSYkzZz8E3zlDtXbluBiYZHJo6kCacux/PXUecTYzOQ8lHFSLVQZ/zKEdcC5yZDxLU3TnoAbBySqCl62XwDN7MP0Id9EuBXxOnbvt8YDlk21MH0uVuBXYSYwr8pcPZcQ2wjBgPgeeiEqQq3ABMMzg8dSCNmHcUOuFwHvsaPjhQ8/Z7NhQl0EkogZZeEaWGtK+eeD5wVn5zrDPVP1cg78iyhJpfwwanAMfh18SN1D6BVvJRocX4Q2DdZyXOAbZcxHXL8IfAureQnITPTOk+XDEl0NIrrgXOM9g/dSADyZOPk4Hr8JvhIdR6C+fnk48bU0fSIyLUkF4I3Jd555B5wNk1fwicgpKPSuQzgQup/zV8AT5aDjH6V58L3AXZfakD6TVKoKVXzAeGU+866POBB4D1+Xayd1LvEcepePKh+udq3ITvqjk8dSANnMsLo7k34Q+BU9OFM6hzgE1KPipzE/Wvg96zn/J86v8QeBa6ByehBFp6Qgbb8J39LkodSwMvBX6dr1iH+tdBnw/cCdm9qQPpEdcBI6h3d5YLeb4jS3Y/cAf1HsFT/XO1bgTOqutDoMHx+K61fV2FbsQfAqcnC2pwF6IuSEkogZZecgP1rmc7n73LIeZT7zZKqn+uUP4QuAyfYq4dgxPwMqQ9r4m6PwRqA5Vq3YyX0dV1Yd7ZwOPACv/X7AG8HKmmnxs2Di/7uzZ1JL1ICbT0kr466Nq1Jcp36DoFH2Xscy1wAtiENFEN6jxU/1y1a4BLUgcxgAuAh4Hle3yvxgsJbc81B1KBzDdTWQxcmjqWAbwEmJfBrj2+V+c66IuAzZCtTx1IL1ICLb1kPjCMetZBnw9sBda88K3sDuAu4OI0ITVip+HTnTcM9ptSqOuBuQYHpA6kH+cCN+9RggSenJ5Y04fAi4F7IFuXOpAecx31nVl7CS8eFMgTaMsSxDOYi9DoczJKoKVnZLAdH/24MHEo/bkAuHGf5AP85ljHuu2L8ZEPJR/VuhHvTVvHxaUv4UUPVNlGYAP1HDW/GCUfKVyD10HXapdKg6PxWudf7fOjG/CfTaw8qIYsw9fNaAYlESXQ0mt+TX0T6P4WgvwSuKiGq8AvA65OHUSvyR8Cb8M/OGsjX3w1BU+O9nU18LJqIxqMDcGTel3D1ZuH71JZt9Kei4HH8EGWPWR3A+up3TXMVOBYdA0nowRaes2vgbPr1A86X3x1Gv0n0NcAhwGzKg2qIdsPT+B+mTqSHnUdNUug8WR0K7C0n59dA1wINqzSiBqbAYyk/4RfSpTBk/hiwrqVpl0CXLdP/XOfq6lf3fZFwDq1YExHCbT0mnl4C6W5qQPZw6XAQ3iHhX1kDwKrqNeHzRy8tdO+U51SjevwOuiDUweyh5cBv8pgdz8/uwafrq/Te+5iYHXeZUGqV8fStEsZ+J72S/whsDYDL6j+OTkl0NJT8inwRdRrBO8y4Op+6p/71O3D5hJgOWRbUwfSo+bj10otpsDNp+MblENkfdPidRrBuxiNPqd0DTDT4IjUgQAYjAdGM3ACfS2+AL0mLRltP3zNgeqfE1ICLb3ol8ArUgcBYDAUuBz4WYNfuwY4t0ajH5ei8o1kMngK33msLrMSM/BazJ83+J1fUZt4bX+8640S6HQWAE9Qn4GBy4A7M9/4px/Z48At1KcOeja+C6wS6ISUQEsvuho4w+DI1IHgLfUOo3E5xI14zfaZlUTUkB2Ex6Gpw7SupT6zKBcDqzPY0uB3rgPmgtWh7GQu/n5SC8ZEMngOL6e7MHEofZoph6hTD/YLgVWQPZQ6kF6mBFp60c34blOXpQ4EH31ekkGDWszsUeqzDfkFePmAko+0fgbMMjgudSD4jMRgnQDm4fXRddjR7RJgcV5aIulcQw1mJczbQr6Uwdd0XA3MBDuq/KgG9XIaz1pKBZRAS8/JRz+uph5lHJfReOq7z8+AV5YcSzMuAW6C7KnUgfS4RfiIb9JrON/Q5XwGTaCzp/GykzqM4F2OSpDq4FfARPMdWFOaiy+KHqykZyGwjeTXsB2B12JflTYOUQItvepnwGV5DXIS5m205gK/aOLXfwzMBjuh3KgGpfrnGsgXnF5F+oeq8zycpmYkrib5rI8dj3eR+UnaOCSDFfhOq7+ZOJRLgKWZt2FsINuFl3mkroO+DK8fn5c4DpFaeCf+gViH+sCeYHC0wS5LuKra4HcNHm4uibcMbDPYO8qPbMAYxoIZ2Ix0MUgfg1cZPGEwImEM/2zNPQDi23mbgZ1ablQNY3g72JZ8IxVJzOBfLPF6CoPFBh9u8rffAXZPuRENGsM3wb6dNoa2DcdznXNSB1IE3USkJ2Xed3kpaUfELmXgxv37yAwfNfuNkmNq5ApgE2T9bZYh1bsWr988P2EMr6DpfuDZWmAd8KoS4xnMbwE/hay/ftVSvauAc80XUlfOYBTeRebHTb7kF8AosJnlRdWI9XVt+mma84vUj0agEzD4B4NbE507M7jX4O0tvOpysKfzThgJ2I1g/5zm3NIfg6sN/iXRuSeZDymf3sKrPg42v7yoGp77ELBnwFKXDEjOYLjBNoMrE53/jw025b3Mm33VrWD/UF5UDc99AdhOsJFpzt+xrhqBlnpQAp2AwTkGu1N0MjA4Mz/3qBZetT/YdrAryotswHMfB7YLrA5dFCRn8KcG6xOd+68NVrb4qrPz6yhB9xB7NdiTYAdUf24ZiMF3Db6R6Ny/MPhci6/6K7DV5UQ06Lk/6QMZYSmBlsIpgU7AYKjBQwZ/mODcnzLvStDqK78P9m/FRzToed8G9mA+hSg1YTAuHwWekODciww+2uKrhoDd79dT1ewrYP9V/XmlEYM3GmytekG3wWEGz1rLrfRsYl7LX/l7zhN3+0D15y2MEmgpnBLoRAz+1ZprI1fkOTODOw3e38ar3wL2QPWLoOzHaRJ3GYzBGoP/XvE5x+QzKG3UgtqX/Xqqkg3NHwATJO7SiMGRBjut4h7h+SLuR8y36G711WuqT2Tt1Dxxn1LteQulBFoKpwQ6EYOL85t3Zc3xzXdB3G0wuo1X95VSzC0+sgHPeRDYU2lKR2QwBp8wuK3ic77P4M42X/2q/Ho6sNioGp7zLLDdpG8DKf0wuN7gkxWf89vtl47YJ8DamEHshL0PbCNYC/XataMEWgqnBDqRvIxji/nfQVXn/LjBgg6OcDPYZ4qLaNDz9dWOVpjwSLMMZuRlHOMrPOcNBm0uKLUDq38gs38EW1jd+aQVBn9msKG1xXwdna9v8eJr2zzCWflARoUPZLbQr+PQlEBL4ZRAJ2TwBRt8F6oiz3dHZ1Pu9p68jrSimkH7GtgPqzmXtMNglcHfVnSuY817qJ/XwVF+APb14qJqeK4M7E4fwZM6Mjip82uqpfO9wuAZ8x0I2zlCBnYv2LuKjWzA852Wl29MreZ8pVECLYVTAp2QwfkGz1XRjcNgZl6+MbaDoxwNtgOsgh2x7ACwR8HeUP65pF0GHzJYW9G53mFwv3W0j4C9BuxxsAruefYSvPVXgs4f0iyDawz+taJzfdvg+x0e5X+CVTTwYn8Ltryac5VKCbQUTgl0QgZDDDYbvKeCc33MCuk9bT+uZgTP3gD2GGr9VWt5N47dBrMqONevDT7b4VGGg20Fe3MhQTU+1/8B+1n555FOGLw5X9S3f8nnOcLgafNNdTo5Ul9LxjGFBNb4XGvBPlj+eUqnBFoKpwQ6MYPPGJTaXzNP1Dca/HkBR3ttXpdc8jVj1/pIi9Sdwc0Gnyr5HKflifq0Ao72WbAbOj9Ow3McnI90v77c80inDA4xeLL9uuSmz/NHBg+0133jRUdbBtZiK8eWzzE3T9RPKvc8lVACLYVTAp2YwVl5DV4LG5u0fI7LzfuOHlPA0UbkpRVv7PxYA55jHN65oMKOH9Iugz8xuLuz0opBz/HxYmZQAGx2fn2NK+Z4/Z7j9/L3yYjyziFFMfiWQanrLfIHzYJ2VLX3g20udz2KfQbs+vKOXykl0EGdCXweuBZYAzwGPAXcAdwAfBwo8UbekBLoxMx7M68y+LsSz/EDg28XeMSvgP2kuOO96PgfBltT3vGlSPnivucMLijp+MPy2uc/KvCoq/06K4v9FOz/lnd8KZLBy/NBhlLaihpMyDvWFLQYz47L16NcUszxXnT8IXmC/u5yjl85JdDBZMAX8L+0vq/HgbuAu/N/3vNnXwL2qzhGJdA1YPDHeYIwvIRjn2Deb/qlBR71ovzmXcLiKMvANhB716ueY/ATg++VdOwr8in2wwo86gfy66yE9mV2XL54sMWd5iQVg/3ye3ApCaPBP1rhPdPth2DfKvaYzx+7r2d6ZfsUlEwJdDDvwf/ClgBXAof38zuHAlfgI9FGxbt6oQS6FvIavG0GhXecMPgbg9ut0D6nNgRsPdjfF3fM5499IdhzYCcWf2wpi8Gl+Sj0mBKO/RODgheu2qi8vrOEUXP7EL7xRMW7dkonDD5tsKzYeyV9Pf/vKX6xuP0W2NNgRxR7XAD7NVglnUkqogQ6mIX4jlkHNfG7w/Lfn1dqRC+mBLomDD5nML/gYw41uMtKeTCzd4E9ROGbnNh3wK4q9phStrwUaaXBpws+7onllYfYD8B+VPAxR+Bb3v9JsceVshmcbLDD4PKCj/u6fAblyCKPC7Yf3pe/4FFzm5WvEZhU7HGTUgIdzOO0NmrycWBbSbEMRAl0TRicbgW3AzN4pXnT/hKm4ezAPIEusKG/jc9Hn0uq65MyGbzN4DFre5OIfo/5ieJnUJ4/+pl5ojClwGO+E+xhsGYGTqRmDL5pBW5ulT9YLrGO2y8OeIaPgN3hyXRhx/wG2M+LO14tKIEOZgGwAWhmhG4YcAsage5pBr8w+EqBx7vW4JtFHa+fM/wNvtNaQSvB7euU3l5MymJwgMFDBu8t6HhHGTxu8OYijjfAWa4uro7U9svfDx8p5nhSNfPt6XcbzC7oeK/MFyeW1ArODs0f2N5W0PFG5/X73TaIoQQ6mHfzQg306+i/BvoQ9q6B/ovKonNKoGvE4DfMG+2PKeBYfTWpEwsIbaCzjAR7Auw1BRzr1Hz0udDpU6mWwUcN1lsBLe0M/ik/VomLq+3i/Lo7tYBjXZm/H0Z2fixJJR/IKKRrkcF8gy8XcawGZ/kw2CawAhah2yfBlnV+nNpRAh1MBnwO2AXPd9rYjnfh2ISXa+zZheMrFNJgvSVKoGskn+6bb9DRiFh+nIVFjmY3ONv/BptHx90M7PNgSzo/jqRkcHz+EPgHHR7nyHz0uaPjNHm2xWCfK+A4N3fZwqueZHBJPvgwtsPjXJgfp4CHs4ZnGgm2DexNHR7nKLx3+TuKiatWlEAHNQbfpes64Ha8NvoJvA/0PLz2+bREsSmBrhmDs803Vml7ExGD385rn08uMrYBznYKvhL89zo4xnF4y6RSdwKTahj8ncF95l2G2j3GxwzWlTv6/PzZXo/vrtlB5xd7ZT6SPb64uCSVfADi3zs8xq8M/l9RMQ1ytn8Au72zcjr7Et4fveqBvCoogZbCKYGuIYPvmpf1tPPaoQarDf6l6LganPUj+Grw/sqUmnn9d8CWo7ZfXSGvhd5o8Mk2X3+8wXaDDkfUmj7jULBFtN2Rww4GuwusoF3mJDWDM/LR45e3+frXm/ffL3CBasMzHgm2Heytbb5+Dt7W8aJi46oNJdBSOCXQNWQwLl94ckUbr/3jfOq7gG27mz7rCLwv9BfaeO3v5ItWzig+LknF4NX5NXx6i68bYnC1wQKDErcpftGZp+ObA/1uG6/9LL54UJ03uojBZ/IHwZb+Xs135nzI4GNlxTbAmf8yL+U4pcXXDQFbQGmbstSCEmgpnBLomjL4lMG9BqNaeM0Zef1pQSuyW2GX5VPYLSTCdjzYVrDStjGXdAx+ni/IanpmweDP8tHncWXGNsDZPwL2IC3tvmZn5df9y0oLS5IwONi8j/6nWnzdfxosN9i/rNgGOPMQsGvBbqKltnb2tjzxPqG82JJTAt3FjgDW5l9VUgJdUwbDDW4wWNzMCIh5y69NBl+tIr4Bovgu2BqwJke/7Sf51Hk31tz1PIPxeTLcVA9cg5n5qHUFCwf7jWA42DKwJmtfbWReM/q1cuOSVMw7I+00uKzJ339D/vuFtMFrnY3C29o1uUusXZCvP+n2jX+UQHexo+D5bhxVUgJdYwZH51OIP2w0imewXz7tvdTggCpj3CeSQ8FuyeuZG4zi2RCwT+c37opqBCUFg4vyBa0fGuT3xphvmNLRwq3O2Qx8QeGXGi/IsqPBloKtRG3ruprBJ/OZvYYtNvPF208b/HVVsQ0QyWvyWZFBFmXbGfnI8+eriSspJdBSOCXQNWcwLa9p/o7Bi7oE5KN8txg8bOm6uewZ0eFgt+It6frZutYOBPt+PkpyYeXhSeXyxOI5g/f3V9dscLnBVoPrrIPOHcWxs/LSoh+C9fNAasflifNtrZV7SFQG/5gnx/0uKszXnuy06vdyGIC9P19b8vH+HwRtNr6T7NfpjcXbSqClKQfiJSHNfL0PJdC1Z3CeeWeNJww+bN6n9LcNPmDwZF5rWtJOV+2wI/Ikels+kncu2Plgf5WPTq8Ha2lxmcRm8Jb8Wl1v8HaD2Qa/b/DFPLn+lFXSsq5ZNgHfnGIN2EfBzsE3XflfeeJxI9hhqaOU6gzX6AkAAAgUSURBVBj8bX6t/sLgTQZTDN6ZzxDusMq6xjTLLs0HKuaB/Vl+/b4M7Cp8C/tvtVYrHZoSaBnUaey9cUuzX1o9XnN5mcZ78pE6M9hmsNngreab9tSMDc+nEn+cj4TsBFuYl25o1K4H5SVJHzV4NL+G7zO43qCm/b/tOHxh4QK8xdcOvG7/DWAVLxCTOjCYkz/09V3Dmwy+YnB+6tj6Z+PyRHkF2LP5NfwtH4HuKV2VQNfwA79rjKH5Ve+zgO/hq4V3lBWQ9Do7CngGsidSRyLpmd9vhmbwVOpYmmdH+J/Zo2njkDowGAGMzOC+1LE0z4YBB0K2LXUkCQwHngXOBW5KHIt0iXPwp7LhqQMRERERKUFXjUD3QtG6iIiIiEhheqFwvd2ezhMKjUJEREREukIvJNBX4aty+2nlJSIiIiIi/TkBuBevvRmROJb+qAZaREREuplqoAO6j6RbK4uIiIhIt+iVBBpgceoARERERCS+XkqgfwwcjfcgFBERERFpSy8sIuyzE9iaOggRERERia2XRqBFRERERDrW6wn0WnxFqIiIiIhIU3qphCOCI4AdBR5vf+CAAo8nIiLSiWfyL+k9XdWqVwl0PfQlzVuSRiEiIiJSriIHCpNRAl0Pi4BZFFtScyhwLfBXwD0FHle610XAZcBfpg5EwvggXgr3g9SBSAiHAF8ErgTuTByLpLEbWJI6COlcN9dAj8T/26alDkTC+CNgTeogJJRrgY+kDkLCOBr/XJqcOhCRTvX6CPSE1AGIiIiISCy93oVDRERERKQlSqBFRERERFqgBFpEREREpAVKoEVEREREWqAEWkRERESkBUqgRURERERaoARaRERERKQFSqBFRERERFqgBLp7PYfv+NQVe85LJXag60VaszP/EmnGTvS5JCIBjE8dgIQyAhidOggJ5XjgkNRBSCj6XBIRERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERHpUAa8G7gF2A4syP89SxmU1NZhgDX4Oi9daFIz/wNY2+DnuvfIvhpdM7r3SGj7pQ5ACvdZ4L3ARuC/8JvQF4CJwHsSxiX1dGr+51rg3n5+vq3CWKS+jgd+H3i4we/o3iN7Guya0b1HRGpjDLALWAocmH/voPzfDTgzTVhSY7+DXxsvTx2I1M4w4Arg74HN+HUy0GjiGHTvkdauGd17RKQ2PoXfkK7c5/tX5t//fOURSd19CL82xiSOQ+rnKF48rT5QMqR7j0Br14zuPRLakNQBSKHOyP+8ep/v/zL/c3KFsUgMpwJPAycD/wHcBnwDeE3KoKQWHgYO2OOrEd17BFq7ZnTvEZHauB14fICfPQ6sqTAWiWE+8Bz9L+L5N7QATF7QaDRR9x7pT6NrRvceEamNx4G7BvjZXcBjFcYiMTyAf2D9PTAWOBJ4A/BQ/v3fThea1EyjZEj3HulPo2tG9x4RqY0ngE0D/Oxu4KnqQpEAMuASYFY/P/st/EPsmkojkjprlAzp3iP9Geia0b1HRGrlDgZu/fN4/nORZgzB6xO3pA5EaqNRAq17j/Sn0TUzEN17JAQtIuwuW4BDgUP2+f4hwMHA/ZVHJFHtBh4BhqYORELQvUeKonuPhKAEurvMy/+8aJ/vX5z/eVOFsUj9vRIfGXxrPz87CjgBWF5pRBKV7j3SCt17RKRWxuNTZtfzwi6Tw/J/N17Y+UkE4BjgGeAe/AOrz/7At/Br5o0J4pJ6ajQdr3uP9Gega0b3HhGpna/iN5+bgX8CFuT//uWUQUltvQe/PrYB/w58DdiQf+87qJWUvGCwelbde2Rfja4Z3XtEpFaGAR/Ep0y353/+BS+MCons63J8I4ONwKPADcBb0AeY7G2wBFr3HtnXYNeM7j0iIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiItKDJgJPA1cD2R7fPw7YCtwKDE8Ql4iIiIhIbb0XMOBd+b9nwI+AJ4DTUgUlIiIiIlJXQ/AR6CeBU4G34An1H6YMSkRERESkzk4EHgFuA7YB32Pvkg4REREREdnH6/GR5y3AyMSxiIiEMyR1ACIiUrlT8j+PBMamDEREREREpO7mADuBrwEPAKuAEUkjEhERERGpqYOBdfnXQcAVeCnHP6cMSkRERESkrr4M7ALO3uN73wB2AxckiUhEREREpKZejY82/9M+3x8J3AdsBA6tOigRERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERERGRLvf/ASWLbxzjrIgqAAAAAElFTkSuQmCC", "image/svg+xml": [ "\n", "\n" ], "text/plain": [ "plot without title" ] }, "metadata": { "image/svg+xml": { "isolated": true } }, "output_type": "display_data" } ], "source": [ "# define a sequence of points (start, end, by)\n", "x = seq(0, 6*pi, 0.1)\n", "\n", "# plot the sin of these points\n", "plot(x, sin(x), type=\"l\", col=\"blue\")\n", "\n", "# add a cos line to the plot\n", "lines(x, cos(x), col=\"red\")" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# NetCDF\n", "[NetCDF](#NetCDF) | [Summary stats](#Summary-stats) | [Time slices](#Time-slices) | [Multiple variables](#Multiple-variables) |[Multiple sites](#Multiple-sites) | [Using functions](#Using-functions)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The easiest way to load and analyze NetCDF data in R involves the [ncdf4 package](https://cran.r-project.org/web/packages/ncdf4/ncdf4.pdf \"ncdf4 Docs\"). If you know the url of the dataset that you are interested in, you can display the info about that datset using **nc_open(data_url)** and then printing the outcome. To work with the file, we will assign that dataset to a name -in this example we will use **nc**. " ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "collapsed": true }, "outputs": [], "source": [ "library('ncdf4')\n", "\n", "data_url <- 'http://hydromet-thredds.princeton.edu:9000/thredds/dodsC/MonitoringStations/butler.nc'\n", "variable <- 'AirTC_Avg'\n", "\n", "nc <- nc_open(data_url)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "If you want to check out the names of the variables in this NetCDF file, you can run:" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/html": [ "