{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Practice Problems\n", "### Lecture 14\n", "\n", "Rename this notebook with your last name and the lecture \n", " \n", " ex. CychB_15\n", " \n", "Turn-in this notebook on Canvas" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The file - **Datasets/greenlandSurfaceDEM5km.txt**- contains a digital elevation model of Greenland. We downloaded the dataset from the National Snow and Ice Data Center- http://nsidc.org/data/nsidc-0092#" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "data": { "image/jpeg": "/9j/4AAQSkZJRgABAQEAZABkAAD//gBeRmlsZSBzb3VyY2U6IGh0dHA6Ly9jb21tb25zLndpa2ltZWRpYS5vcmcvd2lraS9GaWxlOlRvcG9ncmFwaGljX21hcF9vZl9HcmVlbmxhbmRfYmVkcm9jay5qcGf/2wBDAAYEBQYFBAYGBQYHBwYIChAKCgkJChQODwwQFxQYGBcUFhYaHSUfGhsjHBYWICwgIyYnKSopGR8tMC0oMCUoKSj/2wBDAQcHBwoIChMKChMoGhYaKCgoKCgoKCgoKCgoKCgoKCgoKCgoKCgoKCgoKCgoKCgoKCgoKCgoKCgoKCgoKCgoKCj/wAARCAEeANwDASIAAhEBAxEB/8QAHAAAAgMBAQEBAAAAAAAAAAAABQYDBAcCAAEI/8QAUxAAAgECBAQDAwgECQkGBgMAAQIDBBEABRIhBhMxQSJRYQcUcSMyQoGRobGyFTZydCQzQ1Jic4LB0xYXNDVVkpPR8EVTVGOisyVkg6Ph8cLD0v/EABsBAAEFAQEAAAAAAAAAAAAAAAUBAgMEBgcA/8QANREAAQMCBAIHBwUBAQEAAAAAAQACAwQREiExQQVRE2FxgZGx8BQVIjKhwdEGQlPh8TMjQ//aAAwDAQACEQMRAD8A/VOPY9hd49zqqyDh81tCkLz8+GICZSy2dwp2BB6Hzw5jS9wa3UpCbC5TFj2MLq/a3xBTS1MUlFlqzJIBGpgkAaO27kmS3UHpft64uRe1LO2oIakplbI+oArTSeNgSLKOZc9P/wBYfFC6WQxM+Ydfcq4q4nGwK2jHsYpT+0nidawyVNPlSUEqLpJikZ0lPQGz20kDr2Y2ub7WZPabnS61CZU0ikrpMMigMOxYva//AOMOkp3xNLn2AHMjbzTvaY+a2LHsY0vtL4iFQY56XLImEas0bQSalJ1b35m6kC4PXrttiOj9pvE9TPO6UGVCjQIVZ0kWUhu4TVfpvY6cV3OaxjZHH4XaHbkvCoYTYX8Ctpx7GRP7Rc/EkkaQZY5VFYMsTnUd9QtzLrYaTc9b+mIm9pWec6SOOPLJOTtIVp5G3sGIW0m9lNz8QOt7SNYHAEObn1js5804zNabG/gVsWPYyn/L7PjHFIseWNFKoaN1gcq4O4IIkscVY/aPn8sPNjiysoxJjvBJ4l7H+M77n4EYtDh850H1Cb7RHzWwY9jDKn2q8UxTxRrQ5Sxcn+Sl+aBufn+ZA+vEVX7WeK6enkmbL8oKorMfkpewv/PxKOE1J2HiPyve0R81vGPY/Pye2PinQC1Bkwa1zaOW358S0Ptf4nqZHX3HJxpcptFL2t/T9cOPB6satHiPyl9oZzW+Y9jEJPapxJHOY2pMpF4w6/IyfzrH6fwxTzT2w8R0aAJSZO0tgxVopB4b2P0/XDRwmqOjfqEgqGHdb3j2PzgPblxSWCjLcm1XtbRL1/38fKf26cTSPIj5fkqst/oS/V9PEvuOsH7R4j8p/StX6Qx7H5uk9unFCoNOW5KWIuLpKAf/AF48nt14mdQVy/JrHzjluPT5+Pe46zTCPEfle6Vq/SOPYwyh9rPEVRTiR6bKVJ7CGTb/ANeK8ntf4kjdo2o8nDqLm8UtmHmPH0/DEfuiqvbD9Qme0R81vePYw/LPanxPXSKBR5QFB8ZEUmw8/n4scSe0zP6DKmmoo8tkqCrKFeBo9LWJBuZNxtuOuIZuH1EIJc3Tln5LxqYxunHOOI6t+NqTKcthZ2gdeYkiFA+oXJ1ahsI9ZGxuynY2uHfA/JlWooaKtmVJKt6dbzFAGswDEDuBe23oPLBDFJTr2Ej2xc3/ACKk92RHn96p+WrtpUtzFsCcO+E32sxiXhMIxYK1ZTA6SRf5VdrjEsAvK0dY80yX5D2L8+5gzZpxFDQUkMa1Qg5dS8lmjiJ8ZuBu5Gq3Xcntgrl707I1JFNJU1FPTIrVD76lB0WHlYjoBb1OLuYye55cBRPSRVzMkdMtRZY3sb6LD42HQDV1F8CeFMoq8pppJK7TG06iIQGzPEEc7kgmwNydI9DgpHGYKu7BcuJLjyFshrbtQQXa+w31RatiWspoqKalMtLPG915vLWVlH8UTa+97/ZivlucO/GMOR8ySGOjAq5KhtXy6Imo67gaSVADH6V7+mJBURQ00TT1VQj1RlaOIRh4wEGwDKNidz3JuNxtgfFQ0tTUcR5lOGjlqEWnhkaUMNbKCzC2y6bgd9I8zfGQ4xVtrJzgcbaWtuD256na+a0PDo+gifPK39twTb160UsbDOJZ63M0cLOXnijdSA7BrKxOxC6QNPkAD1IsaqK5aLK6yQI9PIiBGkGpi5BJW5sSNrqSVO3cbYXEr/0dU0VNJMJqvUsRe1jNdQoZewIIIK+W64K/pCnqVSeSddSHTFJHdmLbggXO+9utx388XOI0xqI29AMWEBts7gjkBrzOoBv1E0KGdkcodKc+7PtJ9dqCRcU0GdGKkg0Q1k6ukh1CNYogGNwx6lVvZiQL72AvaHJOKpYTLJLWpHS0pWOOhmCJJK/Uvq6CxAvv1Zj6YP0vDdNPLNUyKprdMiT6NMReNzpIZAoDDSTuLE26bnFeo4Agr5Zax4GR3QLJEhSCNSpsWXYjVZTtsN8Z9/8A5uwPuCPXZ2LUxPpHMc4D4T39WQGYyAvkNlI2YUtdLJTQVLoJWdnqkVmgTURqQlQFLXGkWFzrv1BwVTSzulwCnRQfo9j8COnpgVw1wnTVVLK7U+b5ZHUi0UgqQGiRd7hbbFrbk9nwcz2CHK54qlJEk5UQScS/xzKBfVcAA9ALbdcHuB8TkhnbTPJIdlztyta/ZrZAa+nhLOmp9L6EEHtz/CGrEHrZ5bC0doF+Pzm+8gf2ccV4DUkkfeUrH/vMB/zx3FI0VNHEVvNa8nqxN2+8nFdpBJLTar7S8zTaxIVTc3O2xK+uN4btYS5CVblp4pWa8aG5PbHGWU0Uc1XFEl5RO5KILsBYHcdsWVhlaCXnu9MybIaZ1a5IAN9RubXHa2/TbFHKxBTQVSVtOasVVVUKsEMm1tKlmLnxaSA27bdum+M1X/qMRAiBtyDbO9vXgrVNSOmeAQba5a27O5WXQTVkEkMbSoOZCzxlHUNYNp2bY+He+w3vbArO8oSdqhyj88AxujizKOn1j1Fxi/PHWxmmkp4E5QnWrqFKghUkuoGwBZruBcgmxO42wZr6F1yeKr5XuskULvIlLCqu7mxuxa50AA7d+u2KFH+qpWSDpgCDyFvqT3Wt3q1U8OMcfSMv2Ei/gPzptdZjScPZjU0rTyw2cki2rxagSDt8RgLUQtDJrIPyZ0vtY6fP6jv9uNNyKreoFRI4IhmdpqdnUIzILKSw2F7gHbsb49U5HQVRdpYhqc3LA+eN7S8S6aMPcMjpzQ5st89lmfMeM6XUPGCbA9sTQw/KCppXUAG7q25X+lbv6+nwwTr+HqqikljQCRI91a++g/NJ+8fEYDDVC91up/HBNrhK27CpL30TrQuzQIXZDIOujpjqrgWoI5g8QN1YdVPmP+XQ4V8rzI0jJGwPIJtGw3t/QP8Ad6fDDHFmcIhDynSNhuDc+gHc4pSRlpuoHNIKu5K0VMXapYRSBegBs4va6+fYW63+o4LyhjVVdXWRaKmmQqp8JUAoSY1ve7gLdiLbhey4C0jiomieoHKdG1QhSLwG1tV+7dfQA233JNV1VIuXzwyBIxypSxVgrOLXBRiPCA2xUnYN3FsZ3jjah8JbE3FfW9tO/wBdm0kJYHgvNgFvfDu+QZbYkj3aLc9fmDBDA/h0WyDLR5U0X5BghgEETXsJftbMn+SaCAhZGrqVQxQPpBmUE6Tsdr4dMIvtodI+B5HdymmrpitupbmrpH1mww1+LCcJsVHLkx3Ysllp4M64fDzMkct2iEoQNypAbMdN72ItsfP0GOKlY/djTvO8loBG8j6RK4IKhyo6E9vMjHOXzGCm5FXEpWJFm5sALq5c7i3Ute428u2IqZcrXO63MY5WqK0oGkp+YryIqKAQuw1MdIub+g23JqinikYyQEFzhZxJtkAdt88v6QfLdXKeGeoMYqUkeZAdCRoFRRvfwg/O7X6Dbqd8WatWiSKlWEU2WqmlSRtdV2UC2+kWJA6WtYk4F5WZquCGXM0hfdZAyHRyydRYHT0CjSL9d/PB5RDWUE6M9HHUxsHBkmCLEgB3UattvPYnqcZWrFJBIx0bQWOzzvcZ22I+99b55X4pZZ48GI9Q59SpT5dDJQa6mFKp6gh2mljFkjJ20sO/807kdrWuFDiziTPcrqInSOE5bGgipp5IleRQANTE6diQbarbD4k4uU3ETVc9RSw0FbU5Esax1jPpOsrsouHtGg26De1vPDFDHQaDFS6Y6OqdXMMT3hmAAZbX+/axsd8Dx0tQ4hoxZk8r3t9eQ6+uyNQsj4eQ6oGtha17d99e48gd1ntXxhxBVySvQT1UFOTzo2UX5cQFrBrfMFv+jhl4ezfNKyClyetlqKmV1M004c6VgO9jI19yQRYDpq7i2DUOXUdJUS1cLpHOV0yFow0UqatZ1J9HcndSOp69MQ0slxzK2FKCaqZSwa3KUXsPGAB4QSTq0m5OL0HBZpLukHRtaLkkG+l9LZ/nwUlTxenwYI2Xv2fZXqlnqYIgZGao2d+SLoT8diT03t8BgfUZfW1CwrBHC6y6pC00pVzy/EwRbb38I+LYY6uhy94446af3hROokCPqfRpNydPa1zaxsGviDhx4JaiqgknVFjLwrzjYlACz2uDvq03Gx8B67gSe34KZ7KRhAbbM/MRvkLAC9stLa3vkJdCySQGR1ib6aDvN7qzWwQSUMctHEqTusaR3iAUL4SSxI7nY/G3mcCo6Z5p9MT6ZaVBfXHfwvIFLCxNiAjHp3Hrhglyg08McVYyzQzVGlHdBe5VrAjoQbAYrUN4OKc2gelk5FUYlSdUtHrWIsUNt7kyML9LrbqcCqerkZGWt2s4doNvInYqzNTRtdrrl1aXB8UPWP8ASMslZVNBHFT1A5tiRoYAG9+hJ73tbt0F6cKZCa2WtMEss0stSgSMAomrSdLdB0JUE9sG814eM1G0FPJyIWeOoMdyxLWIa7Hc9rA+uJKDKlElXE8LSURqX1Ky6tKhU3t8cVJ3tc8lgs3YeXhoiNG3o4xjPxZ+f31Q3NBUrTzRQAwq8DFZFjiJBUbKNrnoo2874sSTPW0iGuqTGkhDpC7iRZl2JCqhAvuB029cfJ4TT5tlmX08MFZUsZWkaxjMMbFdB3Nrgr3F+mPZRTpBBRJKsLMnydPpOoxabqQAP2CWJPiO47DFiKYQtabDnoCSeWeg69ddclUkjE7nsaSbZam1siTrmdRbyQeso62klfMM5UyQRVQmjQyWlVDJYG/RbKwNv6BB6m9qgjqaiOozHN1elp6SJylM0OhWZ1a7vYm6qpUdBuSbDpiyjxy5jMZInqTykDAs/gVm0oSDspNmfa1gB3scU8tzKszOgaSGaS7JIk0yopLWUmTSbkqb3INh1Ha2L8vEp5ohCXDPUnLcm1yeSHCjY28mZAvYa3tYXOm57tNdYqdnOURVtXAEnbmO8SklVgY/NW+5tYPv8MD87yCmkif3RwNKFlIFwTa/2YKSVMFMkJe0SFVWOMAs7bbALuWOKkKzszU29LDGLoTYytGegHULpN17nYdL46JS3pmtYDoAh4vZZ1UU5ppWiYFnZb6ENiPJie3x/HHEQa5508gqFFgRbSAf5vod9+t9sOee5fBGt6aPQvVgO58ye59ThVqYRrGmwkT5t+h8wfQ4ORydIMSla66nyqsjRJLyPIymwtffBiTPEgoJY3dgTG7KPpIdJ3B7YDRgtT85P4oddxceYIv1xFXTRz0FRCSrJyX6ixHhPQ4V0bZDmLhNIubr9h8Ofq/lv7tF+QYI4HcN/q9ln7rF+QYI45wUUXsIHtykEXAMshLDRWUrXVdR2mXoMP8AhN9q8JqOFBEpKlqunsw6r8qNx64imIEbr8imSC7CByWQ09EYsvpmqVeKoEK6V1lhFuTYX6npv5bDbraneoqaDkQ0VFPNFOLmoUlWQgkcvuWbT84nv6WxBp5FKaSolmJSwh1REtYkmzMOpuRud98VY6jkAF41mVfm3vYA2Okjy9MWn8Pp+K0rPZnAPG3ncdufeqFPN0DwXC/MFWsxytYj+kaZkpY5ZEd6Vt5FkIF9LG1z4twdr79cUc3y2szqvpcpoqeSmy6ZxLVOiXfRqIBc9d7E9LaiNjpwbWvEdLZY2likeNEhkXd7sBt5tdr+W3SwxXqYahIJ4J0qFLuyToZLpMASBpI30WOw/v1Yz8dJUVMvsTzZ45nYC9uZttr4XKJxvZSO9rY3I3sOR5329ckv03CGTR09qdiIZ9SWmIEiKATclDcXK33up9DgxHKEUGNui2ckaSx7km9t9tvTEc0jzypCt0QDSFHYAf8AQwQowtCwliMSSRqX1zrqSw+cv7TA7Gx6W740o4XHwaN9ZnI4aDTXLv593PNQSVs1cWxSuAHr1355WtQqJ45JEp9Sn+UkUG5YAmy2G5JIP1K2LVPVxVQYxVNO8erQxEq7nyIP/W+GGihhjlRqvLqV6t/HO4YBUJHzVB8gbfG/njit4Vo63LWE6xzzazMt4lAC9kG17D1vgeP1Y7FnGLdv3/pNk4bKwXBBQASLw9Ry5oJJIsvhk1tGjK8YYkqdKHwkkndRY7dsDve0fM/c5TSyNOImp3imMb1EzFmZl1EgNcd2FyRbqcKXFuTR0U8FPHC6ozltABKqNiTbtc6QcCjTz1FQtPEvPDyJBGhkW7O17CxI7g72sO9r4iFSx9T7TCy1xYjY31/zvshEtXIwdE4aJ2lzyObM55bTz1VPKjCWqaQFtyGjB+YSG6EG3lgnTcbpLXx1dU1J7ozSyyRwygyP4vCXuLgCwsPS53wtfoPP8ljjr4ORO0aK8kFJIWkjG+1h8621yt7X7gXx9bK6vNsny2Wly33+ORAzuFTmCVWIYWY33N2+31xXm4e9znXjsdbC+nUOrLsT2VMrQbXB17U6/wCW4iE0rlZKVXCpIRqKi/X1PXfAqkzyrzPLeIqmCsmgjSueSBlUjVcC+/pYbYo5NwjHmVQKjPMp93plR2EOsxySuDsSF+aBYnrvttbFSi4NqKqqz/L8rUJSQVpip2fxgtZWYSaRcAKw8W1msDe9wjeDyYA8NsTt2Z/VSmqncMySM0c9n+aCkadqiS4luZDKdJP9K+GigkyozVT0lZTo6zSBeX8sfFpcWAvf5xuPO+Me4eq2mqo1NPTc1NnPKW6kdjcdb4Y8yz/MIaiopYKp6eOaKOQiO130krYX2v0387YHS0+J9t1PRcS6GEjOyPcOVlQM2zmkogI6masZ5TOiJylBIDElr9NlUruW32xytA9NWVVIaiXaBZlFNZtMugI+qRrmwHKNgFuTfFOsyDKaiojqE5pt1kSUN7yCb3du567j4dAMSVkiU1GiQRKklG6y8ynGuRYgQHTQu+nSQCL+RIxoGcFlp3dK7MZD4SAes3Ow05596RlQcBY4De2ZPdbTv6kUy+ljQTJToFZhqLXJZiP5zHc/WcVaxWjCzIpZ4rtp7lbeJfrG49QMXnJpakiOXWrH5xUC199PU9AQN98cTEc24tsbjGxacQuNCoFWlWOpg2IaNhcEdwemE3iKiFKwKbhtwcN8AEMr0/Rfnxfsk7r/AGT9zDFbOaNKugcaSZE3GLVPLgcL6JQbG6zSX5NzL9E/P/8A9fV39Phj5UsFpagnpyn3/snF2rpJqdrkXUfSXfAqr+To6iP6BifR6eE+H/l9nbBhpB00VgZhftrhv9Xss/dYvyDBHA3hr9Xss/dYvyDBLHMTqr69hT9plV7nw7DPa4Stpr/AyAYbMJPtfDHg46QSffKbYf1q4cyFs7hE/R2R78k1zi1pcNQs8zrMZKvK5abK6gwZlVwmZAF1MLHr5jVawNvLpjPeDc9zDMMxbL87q+Ygp3MRnA1Gx1HxnewFz3P34P09QIpCZFZlZSp0nS1iCOvwJxJJMamWCGgEMGXUsAp2haO8jON7q4366b3O9u24xSl/T1XTVLW0+Y1DtLG/VcjbtU8PFI/ZnxzNuToAOrXx71PlDM2qZ9QCgiFSN1U9z3BP4beeJy4RGJFmI2A7YgQGPcnc7k4+lwTYnG4bAxry8DM72CB3JFlEpVYZajrJGfm33OJDHUVPIgaWGSNijyxpqUrYkKDboLozeun1tiCd0p/lZLmIbuo7+g9ew+OC2TRZhGrzrDHNAZrVkSWMjvYWVAbXspFiD2HfAj9QTugpw9pGts9wRpfbntpqpqYNL7OBI6vO26kWsRf4VVwT1MbU4ZnSbX4B0kU2sb9Dfe49MNcOqmZnA1oikgq2rWALj7RgPQ1sNHUSQzUzw1jIsuuUkPKDvY9ywvup74OrUciklqJlbQimUggXKgXP29Mc3le2R12tDeoX+5K0bGSsacbrjZZ1xjBNmki5zlxjenU8khtrEbv95A/s4VOEYqCbivk5jFUuJWD07RvoCSpd7ncG9gQLX69MEqz9IZDNJHIyzQyOebEd1EhNyy9upwKzUvSVdNX06yNJDIkoMZudKsL2HwBGDFFJ0MjScx9lj6pwdJjIzGoWgZyqGnmNOwCzlURF+jrYLb7CcWaqsNQZDsXaRm1lFBCkkhbgbgXO/XFOvqKSu9yqMrl51BUS81JACAwQOSNxsQyjbzGPo2Fhe/ljobQ19naqdWqIAz6ZHdYgpZtJ32HQDzwt53S5nmeRzQ0JVMzTM2mlSGcqPALgKWsSd0Nm7jbtg9JOzOJFLLJa9wLb+ePmT1Pu/wCkWKKXkqZHvpHichQSx7m3fEc8JlbhvqmubiFisxq3zTNuK4Bnc/uuZVckcbNNT8sgdFJRQCQelx9Z2wRrKOaTiFYHno2kpTPTzq5K6NFiWsex3IP44eqiqWmgeVtMUUKl+aVDNGvcK3Xv0HW/rhdzMU1NlsmaSxBObVNPK1tbXZyBcfBhcfHGV4rSMpALnE5xy5261A2CxNzfdTZEjU+R3Yu0MJZIwr3GkbDQxNjcnbtcYP1kAo4ssWojRKamutRPEqxI2peW6G4vYeI3O298Bp5KWlqpRKkkNQtJ73Q6FL6ZSL3jUiyjUSdxgvXZjWnKpo6inhSSpQSO0xaTxqoLAqLAG1iQDbe4OKbqySqlp442/JzOpNrnPlsM0QhaxkMjScwBp1Hflc2VRKw1tJLUS19PVqJBGZI23RiQuhh2sx2PcHzBx6z30km42OK80j5bmFVHU1dNOKmP31Q8oURjSoZe+q506SOhbe2Lx0uAwIYWuPUdsazhNe2riyIuNRnle+57FAWlpIKrVmvlJJGuqaJtSqPpbbr9YuPjbEiuk8GtCGSRbqemxGJCt+g3xSZ/dROlgFKmWIdr/SX7Tq/tHywVGqRJucSMA6oSLdbHACrkU0FWs8YkUwuOtiPCdwfPBXN3vOem/bAav/0Gq/qX/KcHYmjCAp2DJftjhr9Xcr3J/gsXX9gYJYG8Nfq7lf7rF+QYJY5kdURXsKXtOjp5uGlSsj5lOauAyLY9BID23HTrhtwke2JynBh0u6OaylVSjFWuZl6Ed8RyRulYY2anId6QkNFyLhYpXI0YSKJ/lHVpgWO8UOvSGN9jbZb362NrXxzCfdwI0BQJ4QvcW7YirfeBHR18McE9TDPUCuSW5V45LKugHYxhe2xG5PUnE9ZTx0k5hgULEl0RQSQApsLX7bXHl0ubXwa4NxCaSd9FUNN25gnI4QQM+fO+6pz0zWQxzMNw4ZjkR/VlZ55dd7XxV5p1dccK5AxUq5+TG0libdFHVj2H1nbGkDbKmr6vzqpEBukJDt6ufmj6h4vrXBKolPJjWxbRJzAO2oC1/O9gBt2GFanqXpgoLBnvqdv5zHqf+uwGGXLqlamEOBe53BxBUUzXttI24Sq1T1c1UvIqqyKoiHhVJ4bAEtcXYC4HUE7m3TzxZzLMaurrUy2oWkpBI/MeneQmRoY7MLt83SxMdgPW+BuZkRCN0669RA2vjiugbNY2ZJuW8BWCL5t5EJHMQEkWJLCx7aDtjLcV4FTthM0DbEKdtTK1uHESOV1HxLWy5jEauXLlSlo30MhfxyP5BDvb1+zC/wARPVtWRyJRTItl1xKtzGh6A/HGlZAopqd5ajMKeepSNdSPMjcnVsdTLcE7AXHXBrMkpVo6j35IkSoBikZU3YWsoa3l2xixOGG1lM6gM7C8vsT2eur+1jnC/ECUWaxZZWLDBlRE0vOa6mF2C236EEx2tbqScPoiZAGZAYmICufEGUgEMCNrdvqwpz8NJrrJWjZeTpijLGwNl1N4f7Yti3whDT01FWxRho65uU08LAgKu4Uob7pct63PQd9dwXiBeRTHTt9XQ1rZIyGPRhjY2xQkzOky6LMGqamATq8skVO0lnlNgFAXqbsLfbi+wBB9MJ/GVOI6mlrowxmkqJqdw26aV0FSPI+Jtr7/AFYO8QnfTwGRmoTpXFjC4KOjrqqsT3vMtcyBllaHVoRmXZTo2AA7D4nftdr6marpMwq5neCkjiYQ0xTQkm+w36hnWxI7X3GA7NFV0BpYF5U00wDzSkhIFJABY9hffDRwm1JmebwSS0CPLEnLkiI2YbqFj/mrY7fb1N8YYy2JnlGI7X+/f4qKC7nBgOvr0ERZJMu4WnqsphkhzKZQRX1GkXTVsELE6EFxpuBt3wvcFV0uZQ5mGqppmVjKWq9TBpCVB1OBptYMB6Nbphq4moaLPMuy2mrKqopqySAUxnSMOAwuvj32Uspt522xmPAsk9NxCcvZ9K6ZFHQKzBSBe/bA2B//AK9LexaQfAreU9K0UL4gM7HbqTvKY6eCmjmgq11O8sUjvzDHSkaBrufAPpKbAbdLg3JKpcJquHB5bbdWAH47HFOmzKirmMcVqilkaSGYU8nL0gFQdvpbMTe/aw2xWgzRBEq1syx1CVjULFdkMo6MSfmgrvv62xr+B8RbITG8Ybacj/d9O9ZarozSEYjr4ooVdelsA+KmdKBLEaw2tD6jt9YJH14OxsJFBuDcAgjcEHvhb4xmCokd7HTe2NfAMTwFVGqS5ZRORIL2bffFSu/0Gq/qX/KcTAaJiv0X8S/tdx9fX7cR1w/gNV/Uv+U4OM2Vlftbhr9Xss/dYvyDBLA3hr9Xss/dYvyDBLHLzqry9hM9rEkcfCqNNcRGtplZghfTeVRew374c8KvtJhWp4dSCRdSSVUCsPMaxiKZ2GNzuQK8Wl3wjdYxJNSfpIU/KK1TM1MZDLezq9gFvsDdxsN999hbFOuaWZIppFXSwuCq2Budz9t8E83yeOlzqKshjjlaZg8jFBZSCoMgJ+kAF2FhffEWYzrLTukhVJQSNHUixNr+tuvrfE/6Wmlkme6T4stTrrnY8juFHxJkTAwRC2t/Xq6ByG22BdTJzKjr4Ifvcj+4H7T6YtV85ijJA1PcKq+bHoP+u2Bh+TRU1ardW8ydyfrN8b5jblDAF5mvg9kTNEAWcBD54XlFzhmy2NXpEZCLgWNvPCz5Nslcp8xqQI2lIJSMFvjboMXcqo/eMglhmiilKx2cSjwF76tTHqF1DqNxcbYpyU/Olp4D0d9b/srv+OkYJxxUy0dakhp+VoYfLMVAa2xB7tqtYbb98C6pwETuoJh0VfgSqgp8wgWpiiibloocaWLyMLnUQBfc7eXmcOVdR14ltFml1uNpUQFXv1Fh+OM7ilGTmnojTGoPOjqRMPFHIpF7r9IMNxv9eNEo62IUy1tbOGVQZgCDsQL6dXlt0645PUA4sYRLh7mlnRPJy5Zf6gOf09TVZPVrHUvPUK9QZAyBS416Rv3sFHT0wj8L1HuOaRpNUySQ1Scl2Ztle40avQG4v/Sweq87qKP3cTU5WNoVVwzDw3AJt9p+zCbxJNT1Nf7nldpWndUXR0LsQAPtIxfoJJIJWuYMwbobXyNL8bTmMlpDAgnse98UpoErKWrpJ2dY3qmbwNpYldLW9emCVSDFVyiTS7K51WFgx72Hle+K8bkQ1bwe6wVC1UvJmmHhV2VVAv1F728O5vbHRJrFliLg+SU6FBjk4/TFJT5RTkpMpeaMlmHhbwsSb+oH7ONK/RsK1Msixcmd1VVlUbC1tz64RZsyqaTiJJmoliWCJiZGlIFd4bjQo6FW12BJI1dsMFJnNHU1GW1GhjPWU/NgDtZmJuNJXpYWxznjYDap2AYWjbuRLhjGOxNaLnyA3HigvE1DS0z1VRmckUUS5hI0fNc6UY2cNp+DW+3zwt53MuX5SCtHBNFVsJG+UQyLUX3kUAfNK7WB7i4wscaZ1WVGeZhDPM+ohAyyrZkKEjbr5DpgVTRTZ1XwQx6Y3LBTIi2Bt1d/q6n0wNEJIxnILaQwlrRiKZ8i96yLiCroauAwScsTuI3bWy9SEFyCbHoNzve4FsNtVGi5jKnuhkim0VkRZQ2hgQmt16qRqTtcbnYYq5fQ0x4dpzXVMOYGFf4PJ/GMkmslWRxvY6dx0BXrucXqiIQLTz6OUkczoAQNQW1pBqBIIsA9r9QT3wY4LJIypDGuAuQDf8eXK6znG3smAIBy325b+u1TslmNjvhS4tYS1Wx8SgKcN7BlJB7dsJvEJX32QKDse+OoUvz3WdbqlqaIuulTZgbqfIjpiKsEcuXOYg2o076w30Ws1wPhbF8pqPriGphMUFUT82WJ7ejhT+IH3euCmKxBUxcv2Tw3+r2WfusX5Bgjgdw3+r2WfusX5BgjjmZ1RFCc8z+iyV6Za1mtM1iUsxjXprZb30XIBYA2LC/XA3jSaOpyGjnp5Ekieqp3V1N1KlgQfhbBnNsqp8zp3im5kZYoeZC2hxobUviHkbkeV7ix3wre0SjXLuFII8ro4RFHUwHkJZGcqVCgHoTZQLHqAN9sV6n/AIv7D5JzL4hZJHEKU1JSmaod0V276nt3soNwCSNhbc2wFZyRDK5XRKi6i1jtZrG/b5pvfbpb1tZLmUed8F1BqNM2YxM0LJMpUh+q6x84dOtvhgBSu0WXu9k56yKyRLGbaulim9lHzj6AedsCeD1D4ahpLsLbi/WO8+tVPxCBgjdl8Q/r1n/gPMorVsjI2qKMlUPm3c/UNvjfA9yScOmZ00U+VoYvEiqHjf8AnA73+u98J8yaHIx12mkD23QAFRr1wZyR3SXUCNB2IOAwG+CeXTe7hpGXUiAsR527Ykl+UrxRyCaSaqnemQkraFWPpux+02/s4sGKRMirSscjVCESKeZoj13sOYehXfoeptbe2JKGP3amjjPzgLuR3Y7sftJxDnMkqUQ5aRHmSBQsygiRiNlAINzcgkdxttgHXPDKd5IUb9CpKTLpZ83yyqzBOVBSxBGpo0KhdhYN5nfDbnUEFTw41gpCnlpy28ILMFs3+8L4AxVeXZDTinzXMFWWucuI4ojEEcjfSDci9xsNumDEmZwzUlCYQz03vah42FyukM25HqoHpfHKpXjH8JuAj9LRSNjxObbF6+yU+Mmgg0HQxcA2jtfYHClwpKicaUTsUiZRLIgZRfmctrD4+Xwxpub0NPmMjVVOzosTaJk03K37qfLCZwyHl4lqxdOb7pp0WAaQ6wenoFJNvPfBfg5ElQxvWgdfTvjnGLmmtklnjE1iwba4HlgNV0MEsVXX1EfPFDNUOqqwZNQVSdVr9CAe3rjvi3MJo+HaX3H3SfL5ZQtRzELrNG6nRvsQNQ7b9N7Xwj5VRRFKmjpRJUaZ2kEQfSrGwGwuBe23wxp+I8TZEeiLLjX8FQSSZ4QLrROF0qKnNKmsmkgqJIYhBT1UcbCF2dwfBcWvbUCepFr2wE4tqMlySpy33SKSonpHJXxFWiiDGyA9CQzEA7WA3vg1wl75lvC1SSDGE+XjiNnAHQi/bqMBs0yio4gyqjrGyuSjqhIRrZv49SwOqx6HrjA1Uz5ZS6Q3P4Wv4JCwRNkdlr9dj1ZBZpn9W1bxLW1bAl3d9iwYk3FrkbfZi/lNNItNHPRVHLzMMWihFQuuVACDZSNjfpe97nbbHfFdDHQZ1NKoRKdZ9ChAy7hSp3I6XN8feH3D5dLTyPNDHJpusDrG9hfUWLC+91AsQLdbDFxskbomg6BH3h2Elg5LSKENNGmsXqKdUMh0qg1WBIJAUtYtfp12v1x2+XwLNIHiEeoBmW5OlimhnA7FgDcnc29MCMjqRPQmNAYhFII4jK4WSZVCANYG4DAWAA6bi98FoUVKueRXiEMvhjRXZtCgDwnUTbe/nf03GDvCGt6dhwBw0vnlnkDz2sSM75G4ssTxEFrsPX6PV2Xy3Gi9SGRqVOY+qRbxyXH0l2P4X+vALiCNXrLjrpAY+uDsZEdZKmq6yKJR8R4W/wD4H6zhfzJrV0lzcajjd098SGocUEa2jjuxHzmHT4DFDNdcuX1KSGz8lgptuCFJH2HB9WWRAqg6z5YjnoBDR1FQ0SMRE4IcXA8JxcEgGqW6/UnDd/8AJ7LL2v7rF0/YGCOB3Dn6v5b+7RfkGCOOeFF17Cp7SpBFw/C5ttW09rmwvzBhrwoe1Jmj4ZjdTYrW0zX8rSjENR/yf2HySE2F0kS5dTCklnfXCklmTkvoeMkANZxvY7XXceYOEfL1kTNcxq2p/e4kLx8xkSRnsfETGB3NhqQXsCbb7NtJUu2UyNI6vUOWELyrZYzYm5JFtIG9gO1j1wvZfRxZWKt2WokRItCm3MCoRfTcAFlt0PXe/njN07T0gba/1Vt0odA57jqPWn+LnJZKetyicUJfkwTy0+lmLaTs3foDq6bWwq1qFZN+uGago8uNAtVSLylm1Te8RPpcE7N4vIFbWNxthVrmm94YqVqI77HZHt+U/djrfCi5sLQbaDTRZxmgXCLqa2C9DT6paeIjZm1t+yu/5tOBlDNC1QiM3LkJsEkGgk+l+v1Xww0Uci1FQwUDTaEXHlu33m39nF2d9xYJxRaIWxPTI09TNCHkjC05kXlJqdt9wvl06jzxQDS/RW+JYmqJCv8AB5JliOtgi3KjoTgHxWmNRTOjx4b7/bMjVOjdheDa9kp8V5W0zCsp1L1ulIonbUGAIvddyocBTfcKLsb3tgpwXl+b0dVoq54Hoo1eSZo5tRjL6QDva5O/mMXmjWpcKRrTUdMb7WRtt9iG2N+hv9mGjJcsookq6eSjelpaVIQsOtmWPZmNr+eoG3pjn3FKT2SQw2OWl7ad19b317s7LVU1e6SAYiDzy8sxa3ZzXWbSQ0tJU0tMspPhEjXsACLhh54zvK8vpa/iZ6fMTreEaoKdtlqTYhhq72BB0jrv5Ww/1+f5bWU9bl0QeauWK5jRNOlL+Eahte2M04oq6SZ4ooxyTF4lIa7IfUjuD3GG8Pf0Mgc5t7bLM8UlbI4PLr+tE71uT0Mkk8VbTaGmjRGkUkMgS5RkB8IIJ229OmFbLsnzCizeuip6JqymMzFZXK2QbEaibb9L22+3A/KeJM2yTMZqSsp6qtBF/dasuXjYsHZ9Ni3zdR2IG98FZuJwlfmIFF73QyVMvLkOqGRkZALMjX87E3FxfbpjU1EtFVtBlGA/X796HiSN/wAWhCeeFzR5hwxQJTzlYW1F/B85lO+ryAN/QnBCozZVljp6GQ1UyfJ6VUFC17dfMYzPh6GbkyxGsmhp5IVg0JJdbdbC3Tf68aNw5T0FNBJyxzKtWDygHdSeh+vGLljawlyO0lVLNG2JoA9bLLfa5mldVVEgqbRUVLLop5ETxShT4n1dCbjp2sMA0ySpqqRMwkq5lzGplZ40n25wDAAswN1vvYnY2PQb40f2q18cWW0xniV5Vl0xRmLUlz1LDy9O+M64LgqJvd555KumjllW7R04K1Fit01FhsL30i/wOHUzrC9t/G/rxstc11oAQcNvt3ePmNVanpsykzalSCgqaVXzCOXnQzNIiSElSUYCyXN77bWAONEnhnkoKWrgYPSshd2ZBexfqWAANt+mwsSb3vhO4djq6PiKQZxWtPXTB2jWJifCodtdgfAtr7EeK4FrYZ+H2jy7MJYJZ2iyuspmSMMTyIZSCN/IHYgYLMgkhb7ZAc47XGd7Hfv07lnuKzhxbA7Q53B8MzqBa3eOS5rtMCxTaSOTJ4gRuFJ0OD8Ov9nC/X/6dKdtm6dsMsgSQMsxDc5byFfmlmFnt8Tdh6NgFmcaJRRVM7pG4JikLsFBcGxtfrci/wBeOhcOqDLE2R4sSM+3f+kAIsclXDSMjcnSB8OmI6tpzSVOpy14mIHXfSe2OKed5FHu8LMCfnyfJr9/iP1D68QZnTs9FUe9SmW0Tnlr4E+ae17t9ZPwwRBzsAmr9INm1dScOZGuVQU9Q1VTiNZGl6OIC4Cr0cnQRYsovbfBzhyWpnymKatWVZ3Z2IkUoSNZsQp3QEWIU3KggEkjH3hoAcPZWFAAFLFYAWA8AwSxgCjC9hM9rWg8KIJDZTXUve2/OW332w54QfbbobgtI5mdI5K+lVnVSdAEqknbpsLA+ZGIpv8Am7sKZJk0pGymgTMhUVi3SCJOVSFrPve7uAdtyNu1lHngTCkdOXhjkjEiry+bGrKrMPq2Y7i/Q2GGbh9EFI+sDTHHZEQ30oRv9w+4YVqLMaarqxJQVbzoWA5jEi52BBvve19tup7YDcMrfZatrjpunPozNT9K0ZhfddTNTJzykcCgjUV+cLnCnV2MzFPmk4aaxpKrKrusytGNASRQNlPUAE7dbb4VpVIY7Y6vSPa9uNuhQgKMsiQyNKqvGo1FGFwbehwcyOmKUKxwySwMoAbS2tC3VvC1wN79LYCrEZ54IFFy7am/ZXf8dOGykhWniCLf1PrhZ7HXVLdSRCqjO8cE4v8AQYxN9huPvxUzOsp4475hTVKU0x0OKiMhG8VhYqSLC2q9/wC/BMCwBBtviCsggnop6WsYCmYARIsILK53LK4+adrfZjNccZUvgDacYs8xa9xy005nK2WatUT42zAy5DmNuvUetlVo6WGannNBmFEIipJh5ihafUWuoJ6L0IB+aTa5HRyy8S1EszR6UQyTNZBqjkHhRdJ7iydbW64TcwWmlMtTXxRyyeKWwCMUHe2q17Lt2ud8Q5A9NlgkNDNUszUq6XaqMSROfFpVAdtybkg9wPPGU4nTOhwxOJyGQI0BAOG97G2nPusphWB7y9wte1+vXO3rrzVTPZTR5xUvQiWOtWJmMUYvdrG1/TfpgPxhlOWw0NI9HVos9SY9RYl+ZE1yZSfogHawG9ulxiSMPU5gFkq6ybMKxByhzxHFEQG1ayUJboCANtz0tjuPhmabhmQPR1T5tFC8gLIELOG7HWQ1l6LYXsbdcXqKleI3fCCSNc8rbZDWx0ugkh6W+Ec+fq6JS8ax02XckVdJLXxnlcuJHSm0BiPAeunSBb49MXqPhfLs9ymofLYoqHMahjLSyrqTSBpvdDvy7sbm1/m9bDCxwpnVTRrTU08eWrTSTj+E1NMyOC2m5U3sW3UjyuO2454woMxfOKzOKKtjEs1S8BEIkV1ZI11C4bYEWNr9T0wQ9o6SEvls4CwIta3Xc5n11rwfiZidmpeHlrq+EpAqJCsmjpa+1i1z1Bsd/hjXqKOekysykxRyadKiMXB2Fj6dMZdwflOZTyJC8ppaGKESyNLCGAI6kXPTyHrgnmPGjxZ3VZbzn9zhOiFoIkAe+97EN/12xj6sjFhBCP8AAqR8gLwDkgHH+dVC1DwSOjyK2xAuP+r4E8MT0Ua1VVmckElZFVmWlpWnMamRjdnRACL2tbpckDFPiqoSqzYsHmqDpHgkckA9zZbb/dgZBTUcVTIpKtKCrl73Oq1uvbcYlo6cubiGXryWvlcGRhhGvL6p6kzShyo1VTmCmnqKuWSQRozTFmAuA7LcC4IsNwPqxDk/FNTm+Yx0+WZWUBivI09UsSo1t2J3st/S/bqcKWYNWJCzyVs00WoHTJMXud97E+pxpvCYWDh/L5Vip45pFMzSQrZip2Ctt1sDt6k7dcFYpqxsrYKY4S48/EnIXQmohpYoDPMMR037hrt2r6aeolRXlmMKtuFgFja+xLNc7jfbz698VauKnpXlIiBkmUOsr+N7iyt4jc/zT9ZxezCuFPKEKagVve9sA8ynecxzv4UhYGw7Kdm+43+rHQYICxrcRvbc7rJXK6hOkF2ayg3OIMyKtR1JU3HJf8pxZEQ1HbxKDivXLqpKgqoF4Xvb9k4ujVIv1Nw5+r+WfusX5Bgjgdw5+r+W/u0X5BgjjAHVGF7Gfe3UFvZ/MFAJNZSd7fy6Y0HGf+3NYn4AmWpNqc1dLzDe1l5y36YY8XaQop/+TuwrNeE81FK6UtW6CGRRaRvAjR9AwPU7hgdu2EXMIJMu49EOXlFhkk5ghMh0FlvsbdSATv0sxPS+GigoIRllBC6wPKArwVEUVpGRrOFIudxqbuRufLEmZUMEef8Av0ksDRSyKPdzDzLxm40iwuFXUPm7db9cBeIUD6Mse4ZOCucCqiGujfyyCvU0hnjp3WdpIWjBZ18QffZrjp2BPTpgBm9KYpta2MbEkEfhgtl1WJ6VFEQRCWCciz6TexACbW2t8ewG2PmdERU5aVReJC7AefljWfpniTXsFMNWi/df8lDq+nfFKXO3J9ZINkGlsyndvoryk+rdvvNvqwwBWY7AD44W6SOSmEQv8ou7H+kdz95OGejlWdAw2PceuNPKD8yor4qzKPEt0B6joMUMzr2oOShiapragh1gJKBY77kMRpJt2uOo7YZ4REMvmExspQux9Bv/AHYQ6ugzLPDJV1NO0EcvyscVzy1UKAAik3Ztxv0NrgWGAXFeI+yx2abOOiZIXWs0ZqfiF/fMtr5svnOmmpw+lZNC7vYhvPa3bfseuB6hcsqdclPPAkCK0VNOjOY7hvDpuD8CTcDcDuDHCGXS1FdDQy0sxoJJuXNrcFFdBrIt1+hbftcjA3PKQ5fmeYU/PnuXGqmWzLpbqLk6r9tVjYYyr658tQ57j8WuuX35eKZJGejbLZS8O11ZS0lRJR5VHNSSXlaUs4ZQBsq95FHYDfr3N8M2SxyfoyibMKhKqoni5hmiCiPxnYAra6qNievXrhCq83zGhip0WqqA62PMlkD6nU7aNgVSwA0m/TBnKuJ6gUMtTnlPmFQKieRqcwBeWVW2tAWa6hWIA67HBmgrYoyGl5sB1WHh4Z3UcUrQcJuiXFGcRUVHPSVcLSiSnIjLIjQlrFQLt5Df4WGEPIYYnq4TTVTAtJIGjiTZLIgDlmNiWuwIHTT5EYZMs4gqoOfVVyRzR1bSaacyrEIWjQPa8nhIs6rbrtfc4vUNFleb8Q1UgkjEc7BVamtZntdjcWBIFunW2KnEpmzRmcO6gLZgHr3y8BqvBvTuFim6aaE5BBlc8vus9TaJJmUlC9rgN5A4zzMq6GrzBcuzjL6Sinp5Bpr2laO5XYoxta1+4w4cXV0fD3CVKsppquoqFIEtQhZSl7gaR302xjWZ5jV5pUNVZlK0hLahew3JubWxlo4OmcS3bf8AC6DwqORsV5LWytzTDxnlT8O57UtHzJaOqtLFUAl9j85NfQkG/lcYHx5DPWCjqqGGWf32BzaMheVIHACkdSN1N+2/bfDRwPE9dwtMlcGkhiqFjgjkQspAOtrDoLXvc7Xt8MMdRyo/c/dLwxJJywqG2gPG0exG/dRe/wB18GaGkq6pgEQ+EXF8rXHmM9vsq1ZxBlK90ZF3i3ePsfWhS7w3wpHQrJLXRmWaaJonp5FGhbsN3PW409F39QMM1VPy4nkIAPXZQO1gLeQAA+AGOmBJUhtOkWAGwA8rYo5sb04BO5PTGw4XweGiOK+J51J+3ILO1tfLWOxSHLYbD1zQiaV6mcvIbn+7ExiV00OLqw0n1BFjj7l8HMcbX88GIqQBgT2wde8NyVNCctppWhUy7OoaJ792Xb/kfrxDNEzU1Qtt+U4+HhODOjRXSRr/ACiiYftLZW+7QcSyKTRTuEshjcEeR0nETZSvFb/w7/qDLf3aL8gwQwP4e/1Dl37tF+QYIYxJRhewle2CmFXwTPCSVDTQm4ttZwe/ww64WPaJJycihkPRauC/w1i+Iag2ieeo+S8GdJ8HPJYm7PFLlz7iSJY1e7dDpAIO3UG/bBCtqYKumjpX1ytzEQK5W7CxJFu4UgHTcXtiGuihqaipihdBGsayJGuzDpqNjvp3B/tYihiVaz3tweeF06SAU1WtrAP0rbX8jgjNSe+qanmta2o6svx9UOZLJSvLYzmMlMP4oQCmjhpYvCImjsTtsRta3/PA3NTzZIYuoZtZH9FN7fbpGCJNx5AYHBeZV1ErbhSIV+rdv/Ubf2cG6CijooxGzx39et1E97pHF7jclQpEB87vvvghla6UZulzbFSUjoMW8ukDKYj1G+CDyS1MVrNNE2UTwCQbxNzAb6Qg3IJG+4BG3niKgrkZY65qXMqdSNCmacRsJbkWS51MbkX777DFiadYEpSyxCJJw0hdSQRZgBYdbsVH3+WFOHPc0yWuqaisNIlDGfd4HlLTvG7AkeIbje5JNhtYDa2MB+ppAahrC3QcznmivDKV9QXdG4A8iAfP11psjzl46yGrgyiSaqjpXlk8barWCAsbm5uWFjuAcV6ioFTN/CIkWSVS7kDzP24s8JZ5mWaVMtS4glRYY4taI212ZjYnqCSPiAMGc8ypZaBGkiEU6kqRGLK46m/ljNtkbe1lYqqGWJubr22WU1iCkrZDBIwSdTHIHUMCtwSBceH5o3Fjj4UQ7QsUk5g06mJRS/TboL2Hxtg1xHTCJkhp43kNxd2sL/AdbYp5clbl+aQTpElVytyjxkq6j0HlfF4SEt1WadFheQdFV97SmzJWzKPnUdNflqCCiktdz07kd99vIDBWj4dnzmvrMwy2pFDRmZ5A3UQqArXt1sQDvi1TUxq1rYpg8VVqMqIYwoZTuOuwFsGuCKeSmyOVqzmywVUksbuFJZFBCgn0F8RySm1xrorlLTte8Nk+XXrGyzTjviKXibPOXTtoooXZYFkIudhc27bAYXZ4WVVdl0o8osAN+tuvpg3mGSZtSZnVtyuXTCSQiVtBQKne426EW772wS4Z4WlesjrcyeGSJCstOA5K1HQ9BZlUedh0xJEWdGImC7zoBzOXn4LfPmZA3pMQDB9d7DtT5G6xU3ISNYY7hQisGWyALtbtcE7dSScDq0mGKckbIolUjsVYNt67EYvSHXYFmewtdjvivOmu6Eka1KfC4t/fjpFFTCmp2wDYW7+awL3mR5e7Um6sy7OwA6G2B1cNfXsMT0smuhgkbYtGpI8jaxH23xwBz5dPYdcW2HdMKkymnCQayPEfwwRBhDBAxD2uSR4Rtfr+J6DHVKFDLqIVF3JPQAb74tUxaonWei96kaGZZFEOnxC5JVyxAI2I+FvjgPxWvkp8IiF3HbXTq9WT42BxzKDZqUolWeoJQRyrqbSSAreFiSOgAN79PDiWvKUkNTHNZCYZLqWBIsvUgXt1/DzxVzOVoo56JdpqttCjlrdFkazDSmylVuQbENsd74rZ7AKLJq6pqJjA7iRKssmqIyDw6wVuBrIBI6XswtfALiHGpXtaIDhJ1533Hq2Y5K3Q0wmeWu9cl+jOHr/oHLr9fdovyDBDA7hz/UGW7g/waLp+wMEcKrS9hW9o4i/QMHvAvEK2nJFr9JAemGnCZ7Waf3rhPke9pSa6unXmubAXkAsdx16fXiOWMyMcwakEeKQvMfxjbNZTTVTy1FcHiVlmCSrOT4l8NipA+ke4J/DH0KGU363wYp8kzSuoqqroqaFqWB3V9dQI2JQXO2k36nqe+AdPKZIY5WTl6wGsWDWB6bjGm4dE2CLoGuuWoUWuHxO3XFRIKeGSVtxGpa3nbtiow91p44mN5FHiPmx3Y/aTgLmed1zTJCMpqFHNLEEOdSIxsQQtuoS/qwtsCRLFWS1MXNqIHp31EaHve23mB528rg2uLYIxWcUhBAV3WdRxYyst70GHzR1xRSRWQ+eLGW1SwSFZPmE9fI4ncDY2TUdkpxU8uFnkTWws0baSCNx9Xnghw4lNU5RRLMppveQJ+QqKg19Dvbc6jffzFgBgaJVjlExLCKBGkdj0AVS34A4X+EWeGmhrEhhQSNE8gbUXFhrsWW50tcgKAOov02xH6mp2vwOHzZ/aysU1QYZLag69yfuHJImrc2jp41CNVsNIsLhAEvYdN0P14I1dYlNL443lMpa5fbQu22EXh+anynMqR2lLGspeeAu43ud/W5O2GrIJpalJpq2dHnaxUOLKq/DGHORIWrfG2wfqLL1VHQ/pXnVSJ7uH5akjw3cAj+7F6BKSWUQRxIqxoVYGLddW1w18U8xzynpG5SGUxS6lWWOPUuu3cnr64CV3EkVLDQrT1UUVesxiYRkCJ1I3Nuw3+O2FL7hV4qJ97gDPqXVRkqy5tRIsh1RLoWQtuyA3sb/Xvj5xTmCtFXUNPVpTfwyXXE8bkkWT5ugHsQR8RiQ8TGWimE8Hu+YwkxlKjxAIBbmqV30na++wN+mKdNM1FFX6SY5paljMI5i6KFVQFUm19urWF/gBgzwvhjuJOs7Jo3FtexD5X+7nkMtjyOYyt4/dU6p5JIYaGamWpEzctxKpKxDsxuL7AX38hfFuUrrLDSZXA5jqukO1uw7KOgHb68V4Z7xtI4Av0t5Y+S1KLAZV3sNwMbqi4RDRvxtzdpc8uXgAEJdLI8YXOuNbdfPtKlYhBqYgL5k2xUasgY3D7g+WBWZ1rTjrZR0GAr1LA2B2wcjgLhcpgF0aNaDHyIidMUsifVqJH3EYJ5Sb3JwpZdLetnDXYeF/tW34rg9l83LkDXIGFfHhbhCQixTKj6Tsd8SpW1YDimYCXT4Gnvy2v226nr1G3WxxSViwFrm+/wBWJEI1kdgMD5YsYw3svKmtRmVFm0dClJaaWVQa+VhJGY1BMZhTYsw5mkm22mx8scVFJG2Z5/S51UJmKyxQ1LU8K7yS2ISMd2kZV1MBYdD0BxdrZKipnjgWihjp4F1iZn3mYrbl2tdVvu1vnADzwr1XDtblNVJm9FV09VTRiTnpUrpaJHRo9QUbPpQ+YO2MXV8GlMjzE0lo05k+eWgI7t1dpK10FmjQnXkP9zX6hyD/AFHl/QfwePp+wMX8DuG05fD2WIbeGliG37AwRxY0Vtewl+1zR/kiDKiSIKymJV1DD+NXscOmEn2wEjg4kf8AjKb/AN1cTUwBmYDzHmmS/IexZhSceUWUZbmmX1NRVJOWlaySMB44VI+ie5v1wIWXRQU4hkZ0CKup7kmwHmBjOeLWrDn2Z8moVI/B4St/5FL98PuZVJgoObM4kZBquBa9hsLXPoOuDPDHtfUztH7XfdyiracxRQuP7hf6D8pKqKunmzasafMK4nmMASqsBa2pQwYC177eaDuLEll1REaEmCaWZNZs0gs3Qep67E+pOw6AZUVEkEXLjr4UdXkBuqgMdYJU+Dsb/Z9eO6CaSSmJlnSZtZ8UfT8B3ufgRueuDdPFZ+fr6qoRcI3BVKFIJ3xeFmBCnC8HFhbFiKsaM72OLjouSYWp+yCoM+WJqV30SfKPp1BV6G/2jA/Ma2ODLq2pjoJdfvFRDJqMhKEkhCbLYk6rgqVC+uB3CdfaWemLW5ouCTttg4ddJm6F5I5Ia+VXCsgujQhWk3vqKlFuPI7Yxf6qpXdBjb+037irNE5rJ24xcHLxRKlzrKqgx1VVHTiWnquQQ6AtosbAWsL33vhlp5KWohaGiRDzGJY2IsPr/wD1jBoc1eTP6qetEdqiUsQwuEN9iD2thig4qqsnzWIUk0M6ctiyltQDdhfpY45uQRktrJRn9v8AS0jMKagyyMNJLLDKh1La+gn4dMZnxpJPQ5nTCkqqMCePmPIkPycgLeo29emG3iji9KXLjT53QVIcxI6yQWCkta5RjsdN/TGTVs7Z7mUCGoIWeTlJNJqJjBO2oLe59B17Yexhcb7J1JE75nLRKOlpZ6iCvy6JhM0JjeonVxFpIKnkgjSxtewve49MTQKJYpguorz2F2NybBRv9mOMup4KVKWhp2haOmgDs6M/M1ajdB2AYs3kwF/jiv8ApGKjasQfKMap7KNgt9JO+N7+laVzYnSYdbW5m2t/X947ic7Z5yWkkDL/AD13q5FFJHIxJBv2t0wNzfMAhaCIKDazsPwxFXZzJNCyRoIh3INycL882o7Y2kUJJu5DwLqSeo7XxSdyTj5I5OI13cAb73ti+1tgpQ2yu0raKtP6URH1hgbf+o4LRyWAtgMlw9M99+Zp2/pKdvuGCkRGtfIb4rOGZUbkcoa3QvLkchextgnC6KhYNcHe+FlSSbb77DF2WXlQqinfFV8YJyTUQOZRh22Yleh88DsyzKSXLqyN9PLeGRWUDqpUgjFNn0g774qVr/wGpP8A5T/lOHthbdeX6p4VVk4ZylZCGcUkIYi9idC+eCmB3Dn6v5Z+6xfkGCOMGUYSpnPGlNlWezZdNSzOsUPNaRGXUxtqIVT1CoNRYkDsLkNaj7U5FrOCw0YkVWracDWhQ7TDsd8N9fl1HmEfLrqWGoTUrWkQMLqbj8T9p88KntZhEnDKNrlRxV04DRuV/lB1HQ/WDiWmymZ2jzTJPkPYsXXKMqmpczlrYaR6kGTxSOA20K220H8cK9fNJNHSRlEiBAdo0+aFUA+Q6m3bzx94i4c4irs0zWry+hr6qj2+WhkWNTphTV4dY3HfbfA9ayMuztzobqEVai4NhudySDufPtgxwOnMdVUvd+51/q5UXVMk4ax4yZkPXcqOZxyNUhoaNJ1ZSWJLXDFr9mA/v38hbHdAHjp9MlOtO2o+AXt28yfh8ALWG2JK+FauAwubK5FiBfob7YGUc8cGVyzQrI6RlmsxAJO3cbfH1vffGpDQx+I8vW104C4RbVtYY5LHzxSoa01RlDQ8pkttr1XBLDyHdD9oxZJxOwh4xDRetZXsuqWpqqOVGKsp2Plh7NUY5MpryxIpw8kkbgaZkfREy2IIG7g3t0Bxm6tbp1w+8P10ObxTwygBlpFppEPSzFjfb4KcCeM0oqIDGRrl+FG+7TibqhOdcNxtn4ZP4NQ1D3hlMdkAJ3AHmNxbFx8uyjh7MUpauKWshqV1MVjLMyA+VvuxPnqtPwxI/Kc5hQJGNDEgRoNjIqgWGrbcd9VxthbirM4zU0jUizzVMRKrydjY+Z7D444rJG6N5Y7K2S3dHOaqBsuLK2e2a64kqKOvzKR6U1UkjAKGSIKOWOg09AbAYLcIqxy2RloJqWeBDEghT5aoQHWbqTuym1jbYfbhUzGWtybM1SserQFPGpYASeek+V774aeBKef3CozirklqKutBhieR2DILhd/Qm3lsrC97DE0EbsTSBe5Fute4jZtI/OwtqEzVTWZo2FpDcmTSBre3UgfZ8BhIBYmqN7n3hyQfPbDZns/JpZGUgSSnSGAt8SPL/wDOE8zIWntsec34DHZaOJsbQ1gsBosALlc8zYgnrinIbE47lbxE+eIHN8E2hStC+MbnHUKsRc/Nv088RbsfCLnyxZgW2ncE9TYWw85BOKllNoGckExlX+JDAn7sEohYm3nYYocvmRtDudalftGL1I3Mp4pOupFO3nbfFV+qhKvUy3fV2XHE8mtyT0x2/wAnCEHzj1xSnfoB3xEBc3TV9dthirXvagqbn+Rf8pxIXBOKOaSWoqkD/un/ACnEzG5pbL9fcN/q9ln7rF+QYI4Vl4iy/IOHMnOYyMhlpAyqq3JCRBmPkOw+JGDOR5ic1omqTAYBzGjEbOGcFTpYNbYHUGGxIIAN97DnR1RdEMJ/tU/VZf3ym/8AdXDhhO9qptwsp/8AnKb/AN1cS03/AGZ2jzTJPkPYsrp64QZfmsf6QhhYGUiJ1j7wL3O++MzQCOgpohKJLRKGcWsxtv0263w0ZnVzQ1GbRpFTlfFvJfVvAnkp/HC1UIwijDBFYKLhOnT6sbOhiwvLuaHA6ILmERiikeBxBsCxj1A9dztt0/o44omnMGo1KVJDHd1Knt3+/p3x8rDKW1RLSllv84q5NmFt7gjf7NuuO6EMtOQ8cUZ1ttHa3r0J73H1dumCDR/6ZaW61PspTM9vlIZB6oQ4+7f7sfFqIibc1AfJvCfsNsd9euPE3FmsR5HcYnseaRdi/WxtgjwzmP6OzKSdhqidzFIO+kAD7iDgOsEOoERqh80utvsxxSK3IjZZpl1DUQSG6m/cHEb24/hcMkhFwtY90NTmNLXU1c9LSlBzmp0DSzKL+C7XVVsTfbr59oa9sygppCahKylYKXNNCICu5s1lXZhYHUdhub9sIuUZnmGXP8hUoyEWMbx+H42Bw+5HX1ctKsiJT61uH0SNGxB8vCQPrJxlq/gUDiZHNzO6gIcwFrXEA9eXholLj+iqq5MrmghSodKZYmkRmZmJBba/0VFvt8ujssS0rwxRwLDDFTkxqjKQqsR1AJJ3vudutvSiy1VBAsbUz1dOJxJeGVS9+1kIG9zc3Nja3SwxMk00MU0kmXzU41s5AkhB6Dc2a25ubdr4ynCuGztrmiVpAZc6eH9IpU8QE1IyEajX168EDz+fVVGMPqWIWH4nC0X8U52/jm/AYv11YzvI7U093JPzAev14DiUXmtDUfxrfRHkPXHTYgAAENYFIzXxwTiIyufmwSkeulfxOPgkkJ2ga/q6gfcTi1cKWynCluit8cWYgbX239b4pxNPq8KqlvOQkfcuLUYn2JaAWP8ANZvvuPwwxzk1yIUos2ryxbyyO8QQAnlO6m3ox/uxUggk0i9Q4/YjVfxBxLFTRhqmORpZNTK9pJWIIKjsCB1XyxUeSTkFEVPV1CIWaaRI79NTAbYGy1SMxMQlkFuqIbfabD78TlIY3tBFGh81UA/biGdizG97d++HsaV4BRGSY/QjiHm7aj9g2+/FXMoy9DUmSaRvkn2XwD5p8t/vxOx7ggE9upxUzGS1DUj/AMl/ynE7WZp4C/XeT5Nl1bw3kyVNHE6wwxSoLWs2gXv537g7HvfB2mp4aWBYaWGOGFfmxxqFUb32AxR4a/V7LP3WL8gwSxzU6omvYTvasL8KqD/4ym/91cOOEn2vzCm4MecjUIqqne3naVdsS03/AGZ2jzTJPkPYvzLxXUxjPs3RYpXsQpZXsLiFAR84fhjmrQLTheXMqgDwqSX+4k/fijxnSV0XEVfHDT18qaIvlIaVnR2MKaiCEN/Fq74K5o8cNM7+8EhEL6igZiF6+Cwvb4Y3ULmhoO9s9FQOVrJPnSLny6qSpc6z49fzuu/zO1yvpcdeoIZaFWlskTxDUfC5ueg9B06fUdz1xDV3inkVsxnUh2UizkA3O3X0JHoLDbFigOqmvz2n8R8bAg9Bbr5jfy322xJB8/8An2Up0U2PHHjj5i6kXMzaIJWHUKbX87bYkUBAFHRQF+zbEM4uir/OdV++/wCAOJr33wm68pYidQIwYoswmp2vC7RnyHTAVGK9MEKca1DbX74ZI0EZpjwmSmz6qJuwV9rdLYr5jmU06WdrL/NGwwNkqkhFltfucUZasvq3xXbAL3AUYaSu6ie5PlgfG12n/rW/ux0zE4hiO839a392LIbaylAspTjpFYmw645He9/qxNErgWGoE99Wxw4lKu4lOknmKD5WGLlONhZmPxGK0aPcAhAAe7Db68EKVCbMVAv31fhiF5UbircYAtfoMQVEmirJO3Mi/K3/ACbE2qw1XFug9fXA2scGop3AW2pk67m6k/ZdcQWuVGApbk773OIZwQCSW9RfH1mJIsbfXiCVjcEgG3W564nATwFG3U2xTzD/AEKp/qn/ACnFodMVcw/0Gq/qX/KcTN1UgX7Q4a/V7LP3WL8gwSwN4a/V7LP3WL8gwSxzA6ogvjsqKWchVHUk2AwucWZdQ8TQDJpa8RNFPDUTrDInNTSdabMGG7KuxHS/phhnijnheKeNJInBVkdQQwPUEHrhTn4Cy45qKynkaJWnEzwGNWQWKnwbDS3gUazchSVWwtp8CQbheIvkqaezWjSMIudZyFHTen/wsDKz2NZJVu7S5pm93vqt7vvfr1hxpcESQQRwxC0cahVFybAdOuO8WRW1DdHlR9Czksjl9g3DkrapM0zlm8z7sT1v/wBz545X2CcOoLJm+eqLk2VqYC56/wAjjXsew8cSqhmJClwN5LIv8wvD/wDtrP8A/iU/+Dj5/mF4f/2zn/8AxKf/AAca9jmaNZoXjfUFdSp0sVNj5Ebg+ow73nV/yHxS4G8li8vsQ4XNelOOIs2E0Vnkjaem1gNdV8PJ7k9T8O+Lg9g3D/8AtrP/APiU/wDg4bs24EpK6peVauaPXpDh1EhYBStyx8TOAfA7Fih6bbYbYIYqeJYqeNIo16IihQO/QYT3nV/yHxXsDeSyYewfh8f9s5//AMSn/wAHHf8AmMyICwzrPx/bp/8ABxrOPY97zq/5D4pOjbyWRt7CMgY751n/APxKf/Bx8/zDcP8A+2s//wCJT/4ONdx7HvedX/IfFewN5LIv8w/D/wDtnP8A/iU/+DirD7DuHVkmWbPs4DGUlQs9PexIUX+R66gR8ca9XUnvaKpnqIdJveF9JO1rHC/w9wZR5HmAqYKqpnCiXSs+liC8hc+O2ogE7Am17ndiTj3vOr/kPilwN5JMHsIyEdM74gH/ANSn/wAHHY9huRi3/wAc4gsOl5Kc/wD9ONZx7HvedX/IfFJ0beSylPYhkaG4znPT8Wp9/wD7OJk9i+SopVc3zu37VP8A4ONQx7CHiNUf/ofFJ0TDsswf2M5My2OcZ3a1hZqfb/7OKNf7FcijhjL57nCfLR2MslOoJ1jYfI9TewHc2xruAXEmQHPKeWCasdYXAAiMMbqhsQWBI1BrEgENt5EXBT3hVfyFe6JnJZzReyDhbMGkFFxJnNRoJDGKWnYAglSLiG1wQdsWj7C8hP8A2zn3+/T/AODjR8hyeDJaeaGnlnlErhy0xBIsixqNgNgqKPPbck4JYd7yq/5D4pejbyWS/wCYrIbW/TOff79P/g44l9gvD8sbo+c5/pdSp+Up+hFv+5xruPY97zq/5D4r3Rt5KCgpUoqKnpY2ZkhjWNS3UhQAL277Ynx7HsUU9f/Z\n", "text/plain": [ "" ] }, "execution_count": 2, "metadata": {}, "output_type": "execute_result" } ], "source": [ "from IPython.display import Image\n", "Image(filename='Figures/Topographic_map_of_Greenland_bedrock.jpg')\n", "# Image from https://en.wikipedia.org/wiki/Geography_of_Greenland" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### 1. Histograms and CDFs\n", "- Take a look at the first few lines of the dataset\n", "- Import the data into a numpy array \n", "- Determine the shape of the original dataset\n", "- Flatten the dataset\n", "- Determine the shape of the flattened dataset\n", "- Plot a histogram of the data- normalize the data and include 8 bins\n", "- Add a label to the x-axis, add a label to the y-axis, add a title\n", "- Sort the dataset from smallest to largest\n", "- Create an array of percentages\n", "- Plot a cumulative distribution function of the dataset (percent vs. elevation)\n", "- Add a label to the the y-axis, add a label to the x-axis, add a title\n", "\n", "### 2. Numpy.random\n", "- Using numpy.random, generate an array with a single random integer between 0 and 100\n", "- Generate an array of 30 random integers between 0 and 100 \n", "- Plot a normalized histogram of the random integers " ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.6.7" } }, "nbformat": 4, "nbformat_minor": 2 }