{ "cells": [ { "cell_type": "code", "execution_count": 2, "metadata": { "collapsed": true }, "outputs": [], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt\n", "%matplotlib inline\n", "\n", "pi = np.pi\n", "\n", "def effective_wind_speed(wind_speeds):\n", " v = np.copy(wind_speeds)\n", " v[v < 4.] = 0.\n", " v[v > 25.] = 0.\n", " v[v > 12.5] = 12.5\n", " return v\n", "\n", "def power_curve(wind_speeds):\n", " R = 287.\n", " T = 10. + 273.15\n", " c = 0.19\n", " p = 1000.\n", " A = 1.\n", " P = (c*p*A*(wind_speeds**3)) / (2.*R*T)\n", " \n", " ### scale to 100%\n", " P = P * (100./np.max(P))\n", " return P\n" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "[]" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXoAAAD8CAYAAAB5Pm/hAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAHKFJREFUeJzt3X10VfWd7/H3NwkBeSYQkOegRCmiCA34MFpt1VprW7Vz\nrw+9WtqxZWaWbXV6Z92x7V21nbXsdd3V623v6oz30uqUjg/oVHvR1s5IGanai0AIgjyI4SmAhDwQ\nCISnhJzv/eNsMGqAvXPOyT775PNay3XO2eck+e61zYed7/79ftvcHRERKVxFcRcgIiK5paAXESlw\nCnoRkQKnoBcRKXAKehGRAqegFxEpcAp6EZECp6AXESlwCnoRkQJXEncBAKNGjfKKioq4yxARSZTV\nq1c3u3v52T6XF0FfUVFBdXV13GWIiCSKmdWF+ZxaNyIiBU5BLyJS4BT0IiIFTkEvIlLgFPQiIgXu\nrEFvZk+YWaOZre+yrczMlphZbfA4ost73zGzLWa22cxuzFXhIiISTpgz+l8Cn/nQtgeBpe5eCSwN\nXmNm04E7gYuCr/lHMyvOWrUiIhLZWcfRu/trZlbxoc23ANcGzxcCy4C/C7YvcvfjwHYz2wLMBZZn\np1yR9724dg9bGg7FXYbE5Mqpo7j8vJFxl5EIPZ0wNcbd64Pne4ExwfPxwJtdPrc72PYRZjYfmA8w\nadKkHpYhfdWhYx08sGgNKQezuKuR3uYOf9q6j+f/+sq4S0mEjGfGurubWeQ7jLv7AmABQFVVle5Q\nLpGs291KyuFXfzGXT1xw1hngUmC+/MRKWo92xF1GYvR01E2DmY0FCB4bg+3vARO7fG5CsE0kq2rq\n9mMGl04aHncpEgOD9Gm9hNLToH8RmBc8nwcs7rL9TjPrb2ZTgEpgZWYlinxUzc79TC0fzNAB/eIu\nRWJgBor58M7aujGzZ0hfeB1lZruBh4BHgOfM7F6gDrgdwN03mNlzwEbgBHCfu3fmqHbpo9ydNbsO\ncOP0c+MuRWJSZKYT+gjCjLq56zRvXXeazz8MPJxJUSJnsr35MAeOdDB7sto2fZUBKSV9aJoZK4lT\ns/MAALMmjTjLJ6VQmalFH4WCXhKnZud+hgwoYWr54LhLkZiYmXr0ESjoJXHW7DzApROHU1SkAfR9\nlZG+ViPhKOglUdqOn2Dz3oPMVtumT1PrJhoFvSTKul0HSDnM0vj5Pq3IDFfzJjQFvSTKml3BhdiJ\nOqPvy8wgpZwPTUEviVJTt5+powczbKAmSvVlhqlHH4GCXhLj5ESpWRPVtunrNDM2GgW9JMaOfUdo\nOdzO7Mlq2/R1ppmxkSjoJTHW7NwPoBE3ouGVESnoJTFqdu5nSP8SKkdrolRfV6TWTSQKekmMmroD\nzNREKSHdutFaN+Ep6CURDh8/wTt7DzJb4+eFk62buKtIDgW9JMKanemJUlUVZXGXInlAF2OjUdBL\nIqzc0UKRaUaspKWXQFDSh6Wgl0RYtb2F6eOGMkR3lBKC1k3cRSSIgl7yXkdnijW79lM1WW0bSdMd\npqJR0EveW/9eK8c6UsydoqCXtPRaN0r6sBT0kveqd6QnSlVVaKKUpGkJhGgU9JL3Vu5ooWLkQEYP\nGRB3KZInNOomGgW95LVUyqne0cIcDauULrQEQjQKeslr25rb2H+kQ0EvH6DWTTQKeslrK7en+/Nz\ndCFWutB69NEo6CWvrdrRwqjBpVSMHBh3KZJHtKhZNAp6yWurgv68mRYyk/eZGSndSzA0Bb3krfrW\no+zef1T9eemWYj48Bb3krVXB+HkFvXxYka7GRqKgl7y1ansLg0qL+djYIXGXInlGM2OjUdBL3lq1\no4XZk0dQUqz/TeWDtKhZNPoNkrzUeqSDzQ2H1LaRbhUVaWZsFBkFvZn9jZltMLP1ZvaMmQ0wszIz\nW2JmtcGjFiiRyFZs34c7WshMumWodRNFj4PezMYD3wKq3H0GUAzcCTwILHX3SmBp8FokkuXb9tG/\npEg3GpHu6VpsJJm2bkqAc8ysBBgI7AFuARYG7y8Ebs3wZ0gftHzrPj4+eQT9S4rjLkXykEbdRNPj\noHf394AfAzuBeqDV3V8Bxrh7ffCxvcCYjKuUPqXlcDvv7D3EFeeNjLsUyVNq3USTSetmBOmz9ynA\nOGCQmd3d9TOeXoyi26NhZvPNrNrMqpuamnpahhSgldv3AXDF+Qp66Z5O6KPJpHVzPbDd3ZvcvQN4\nAbgSaDCzsQDBY2N3X+zuC9y9yt2rysvLMyhDCs3yrfs4p18xl0xQf166l76VoKI+rEyCfidwuZkN\ntPRCJNcBm4AXgXnBZ+YBizMrUfqa5dv2UVUxgtISjf6V7qVbN3FXkRwlPf1Cd19hZr8GaoATwBpg\nATAYeM7M7gXqgNuzUaj0Dc1tx3m3oY1bZ42PuxTJZ1rkLpIeBz2Auz8EPPShzcdJn92LRPbmtqA/\nrwuxcgZFQc67u1Y2DUF/G0teWb51H4NKi5kxfljcpUgeM9LhrvZNOAp6ySvLt+1jzpQy+ml9GzkD\n63JGL2en3ybJGw0Hj7Gt6bDaNnJWp1o38ZaRGAp6yRun+vMaPy9ncbIvr0lT4SjoJW8s37qPIQNK\nuGic+vMSjnI+HAW95I3l2/Zx2ZQyios0ikLOrEgjbSJR0Ete2HPgKHX7jnC5+vMSwsmcV+smHAW9\n5IU3tjQDcOX5o2KuRJLg5Pm8cj4cBb3khddrmxk1uL/uDyuhmEbdRKKgl9ilUs6ftjRzdeUozXKU\nUE726DWOPhwFvcRuY/1BWg63c3Wl2jYSjWbGhqOgl9i9Xpvuz181VUEv4Zh6N5Eo6CV2r9c2Me3c\nIYweOiDuUiQh3p8Zq6QPQ0EvsTra3kn1jv1q20gkJ6/kqHUTjoJeYrVi+z7aO1NcXam7jEl4poux\nkSjoJVav1zZTWlLE3CllcZciCaJFzaJR0EusXq9tYm5FGQP6FcddiiSJFjWLREEvsWk4eIx3G9rU\nn5fITs22UM6HoqCX2JwcVqn+vER1asJUzHUkhYJeYvNGbROjBvdn2rla9kCi0aJm0SjoJRaplPPG\nlmaumjqSIi1LLBFpUbNoFPQSi017D9Lc1q62jfSIWjfRKOglFss2NwFw9QW6ECs9cLJ1oxlToSjo\nJRbLNjdy8fhhjB6iZQ8kOjX7olHQS687cKSd1XX7+eSFattIz7y/THHMhSSEgl563Wu1zaQcrp02\nOu5SJKE06iYaBb30umXvNFI2qJSZE4bHXYoklFYpjkZBL72qM+Use7eJay4op1jDKqWHdIepaBT0\n0qvW7T5Ay+F2rlV/XrJAg27CUdBLr3p1cxNFBp/Q+HnJwPv3FlbSh6Ggl161bHMjsyaNYMSg0rhL\nkQQ7tUyxcj6UjILezIab2a/N7B0z22RmV5hZmZktMbPa4HFEtoqVZGs8dIx1u1v5lEbbSIaMk8sU\nx1xIQmR6Rv9T4F/dfRowE9gEPAgsdfdKYGnwWoRl76Rnw6o/L5ky3TM2kh4HvZkNAz4BPA7g7u3u\nfgC4BVgYfGwhcGumRUphWLKpgXHDBjB97NC4S5GE06Jm0WRyRj8FaAL+yczWmNkvzGwQMMbd64PP\n7AXGdPfFZjbfzKrNrLqpqSmDMiQJjrZ38nptE9dPH9PlQppIz5hmxkaSSdCXALOBx9x9FnCYD7Vp\nPD3ItdtD4e4L3L3K3avKy/WnfKH705ZmjnWkuGF6t//ui0SimbHRZBL0u4Hd7r4ieP1r0sHfYGZj\nAYLHxsxKlEKwZGMDQ/qXcNmUkXGXIgVAfxNG0+Ogd/e9wC4zuzDYdB2wEXgRmBdsmwcszqhCSbxU\nyln6TgPXXFhOaYlG9ErmtKhZNCUZfv03gafMrBTYBnyV9D8ez5nZvUAdcHuGP0MSbs2uAzS3tatt\nI1mj1k00GQW9u78FVHXz1nWZfF8pLEs2NlBSZFx7ocbPS3ZoUbNo9He05NwfNjVw2XllDDunX9yl\nSIEwLWoWiYJecmp782G2NLZxw8fUtpHsOXkxVjNjw1HQS069smEvANerPy9ZpEXNolHQS069vH4v\nF48fxoQRA+MuRQqIFjWLRkEvObPnwFHW7jrATRefG3cpUmC0qFk0CnrJmX9dn27b3DRjbMyVSKE5\nNepGp/ShKOglZ36/vp5p5w5hyqhBcZciBUbDK6NR0EtONB46RnXdfp3NS06837pR1IehoJec+LcN\nDbij/rzkhAbdRKOgl5z4/dv1nFc+iMrRg+MuRQrQqbVuYq4jKRT0knUth9tZsb2Fz84Yq7XnJSe0\n1k00CnrJulc27KUz5Xxmhto2khu6w1Q0CnrJupfW7aFi5EAuGqdbBkpumFo3kSjoJasaDx1j+dZ9\nfGHmOLVtJGfUuolGQS9Z9fK6elIOn585Lu5SpICdOoVQzoeioJesemldepJU5ZghcZciBez9UTdK\n+jAU9JI1u/cfYXXdfp3NS86dat2k4q0jKRT0kjW/XVcPwBcU9JJjJ2fG6nw+HAW9ZM1La/dw6cTh\nTCzTksSSW1rULBoFvWTF1qY2Nuw5qLN56RVa1CwaBb1kxeK39mAGN1+iRcwk9061bnRGH4qCXjLm\n7vxmzW6umjqKMUMHxF2O9AGmO0xFoqCXjFXX7WdXy1G+OHt83KVIH6FFzaJR0EvGXqjZzcDSYm68\nSGvbSO/QzNhoFPSSkWMdnfx2bT03zRjLwNKSuMuRPkKLmkWjoJeMLNnYwKHjJ/hztW2kF2lRs2gU\n9JKRF2p2M27YAC4/b2TcpUgfonH00SjopccaDx3jtdpmbp01nqIirVQpvUetm2gU9NJjL761h86U\na7SN9DotahaNgl56xN15dtUuLp04nKmjtVKl9C4tahZNxkFvZsVmtsbMfhu8LjOzJWZWGzyOyLxM\nyTc1Ow9Q29jGXXMnxl2K9EFa1CyabJzR3w9s6vL6QWCpu1cCS4PXUmAWrdzJoNJiPneJ1raR3qeL\nsdFkFPRmNgG4GfhFl823AAuD5wuBWzP5GZJ/Dh3r4Lfr6vn8zHEM6q+x89L7tARCNJme0f8E+C9A\n107ZGHevD57vBcZk+DMkz7y0tp6jHZ3cMUdtG4mH6WJsJD0OejP7HNDo7qtP9xlP/13V7ZEws/lm\nVm1m1U1NTT0tQ2Lw7KqdTDt3CJdOHB53KdJHFemMPpJMzuj/DPiCme0AFgGfMrMngQYzGwsQPDZ2\n98XuvsDdq9y9qry8PIMypDdt3HOQtbtbuWPOxFNnVSK97eTF2JSCPpQeB727f8fdJ7h7BXAn8O/u\nfjfwIjAv+Ng8YHHGVUreWLRqJ6UlRdw2S2PnJT7v33hESR9GLsbRPwLcYGa1wPXBaykAbcdP8ELN\ne3zu4rEMH1gadznSh+libDRZGTLh7suAZcHzfcB12fi+kl9+U7ObtuMnuOeKyXGXIn2c7jAVjWbG\nSijuzsLldVwyYZguwkrsdM/YaBT0EsrybfvY0tjGPZdP1kVYid2ptW6U9KEo6CWUf15ex/CB/fj8\nTM2ElfidPNXQHabCUdDLWdW3HuWVjQ3cUTWRAf2K4y5HRBdjI1LQy1k9s2InKXfuvlwXYSU/6A5T\n0Sjo5YyOdXTy5IqdXDdtNBPLBsZdjgigRc2iUtDLGb1Q8x4th9v52tXnxV2KyCm6w1Q0Cno5rVTK\n+cUb27h4/DAum1IWdzkip2hRs2gU9HJar25uZFvTYb529RQNqZS8okXNolHQy2n9/PVtjB02gM9e\nPDbuUkQ+QIuaRaOgl26tf6+VN7e18JUrK+hXrP9NJL9oUbNo9Bss3fr569sYVFrMnXMnxV2KyEdo\nHH00Cnr5iB3Nh3lp7R6+dNkkhp3TL+5yRD5Ci5pFo6CXj/jff9xKSXERX9eQSslTOqOPRkEvH7Dn\nwFGer9nNHVUTGT10QNzliHSrSDNjI1HQywcseG0b7vCX1+hsXvKXFjWLRkEvpzS3HWfRqp3cOms8\nE0ZouQPJX2rdRKOgl1Mef2M7x0+k+Otrz4+7FJEz0qJm0SjoBYB9bcdZ+P92cPPFYzm/fHDc5Yic\nlZlG3YSloBcAHlu2lWMdnTxw/QVxlyISiqHWTVgKemFv6zH++c06vjh7AlNH62xekqHITDNjQ1LQ\nCz97tZaUO/dfVxl3KSKhmWmtm7AU9H3crpYjPLtqF3fMmagbi0iiGKbWTUgK+j7up0trMTO+8Umd\nzUuymGlRs7AU9H3Yxj0Heb5mN/OumMy5wzQLVpIlPeom7iqSQUHfR7k7P3p5E8PO6aezeUmkdOtG\nSR+Ggr6PWra5iTe2NPOtT1UybKBWqJTkKdIZfWgK+j7oRGeKH728iYqRA7n78slxlyPSI2amUTch\nKej7oGerd1Hb2MaDN32M0hL9LyDJZOhibFj6Le9jWo908Ogr7zK3oowbLxoTdzkiPaaLseEp6PuY\nH7+ymf1H2vnBFy46tTCUSBKZ6WJsWD0OejObaGavmtlGM9tgZvcH28vMbImZ1QaPI7JXrmTi7d2t\nPLmiji9fUcH0cUPjLkckI+lx9BJGJmf0J4D/7O7TgcuB+8xsOvAgsNTdK4GlwWuJWSrl/NfF6xk5\nqD/f/rQWLpPk06Jm4fU46N293t1rgueHgE3AeOAWYGHwsYXArZkWKZl7tnoXa3cd4Hs3T2PoAA2n\nlOTTombhZaVHb2YVwCxgBTDG3euDt/YC3V7xM7P5ZlZtZtVNTU3ZKENOo/HQMR75/TvMnVLGrZeO\nj7sckazQombhZRz0ZjYYeB54wN0Pdn3P01dKuj0U7r7A3avcvaq8vDzTMuQMHlq8gaMdnfzotot1\nAVYKiBY1CyujoDezfqRD/il3fyHY3GBmY4P3xwKNmZUomXj57Xp+v34vf3P9BVprXgpKkYEux4aT\nyagbAx4HNrn7o13eehGYFzyfByzueXmSiZbD7Xx/8XouHj+Mr189Je5yRLLKDFKpuKtIhpIMvvbP\ngHuAt83srWDbd4FHgOfM7F6gDrg9sxKlp3740gZaj3bw5Ncuo6RYUyaksBi6GBtWj4Pe3d8gPcKp\nO9f19PtKdry0dg+L39rDA9dXMu1cjZmXwqNFzcLTaV4B2r3/CN/9zdvMmjScb3xyatzliOSEFjUL\nT0FfYDpTzrefXYs7/PSOWWrZSEFT6yacTHr0kof+8dUtrNzRwqO3z2TSSN0DVgpXUREadBOSTvcK\nyJvb9vGTpbV8fuY4bpuliVFS2AwjpSZ9KAr6ArG39RjfeLqGyWUD+dFtMzQxSgqeFjULT62bAtB+\nIsV9T9dwpL2Tp79+OUO0lo30AUWmmbFhKegLwMO/28jquv387EuzuGDMkLjLEekVBmrdhKTWTcI9\ntaKOhcvruPeqKXzuknFxlyPSe9S6CU1Bn2B/fLeJ7y/ewLUXlvOdm6bFXY5IrypSkz40BX1Cbd57\niPueqqFy9GB+9qXZGi8vfY5aN+EpHRKovvUof/HLVQzqX8w/fXUOg/vrUov0Pbo5eHhKiIRpbjvO\nf/rFClqPdrBo/uWMHXZO3CWJxEJ3mApPZ/QJ0nqkg3seX8meA0d54itzmDF+WNwlicRKa92Eo6BP\niEPHOvjqL1eytbGN/3NPFXOnlMVdkkisTOPoQ1PrJgEOHGln3hMr2bDnID/70myuuUC3XhTRHabC\nU9DnuaZDx7nn8RVsazrMY3d/nBumd3uvdZE+RzcHD09Bn8d2tRxh3hMrqW89xhNfmcNVlaPiLkkk\nbxiGq3cTioI+T63ZuZ+v/6qa9hMpfnXvXOZUqCcv0pXmS4WnoM9Dv1tXz7efe4sxQwewaP4cpo4e\nHHdJInlHF2PDU9DnkROdKf7nH97lH17dStXkESz4chVlg0rjLkskL2lmbHgK+jzReOgY33pmDW9u\na+GOqon88JaLGNCvOO6yRPKWbrkQnoI+Dyzb3Mjf/ss62o538OP/OJP/8PEJcZckkve0Hn14CvoY\nHTrWwcO/28SiVbuoHD2YJ782l2nnDo27LJFEUOsmPAV9TJZuauD7izdQ33qUv7rmfB64vlKtGpEI\ntKhZeAr6Xra9+TB//9IGXt3cxNTRg/mXv7qSj08eEXdZIoljWtQsNAV9L2k6dJzHlm3lyTfrKC0p\n4nuf/RjzrqygtETLDYn0RLp1E3cVyaCgz7Gj7Z38r3+v5Zd/2kF7Z4o/nz2ev/30hYweOiDu0kQS\nzQw8FXcVyaCgz7HH39jGY8u2csul47j/ukrOK9fkJ5FsKDKjEyV9GAr6HNtYf5Apowbx0ztnxV2K\nSEHRombhqUGcY1sbD3N++aC4yxApOFrULLycBb2ZfcbMNpvZFjN7MFc/J591ppztzYc5X2vViGSd\nFjULLydBb2bFwD8ANwHTgbvMbHouflY+29VyhPbOFOerLy+SdWam1k1IuTqjnwtscfdt7t4OLAJu\nydHPyltbm9oAFPQiOZC+wZSSPoxcXYwdD+zq8no3cFm2f8g7ew/yzafXZPvbZk3r0Q4A9ehFcqDI\nYNPeQ9zw6B/jLiUj115Yzvduzm3DI7ZRN2Y2H5gPMGnSpB59jwElxVSOye+z5fPLBzN8oJYaFsm2\nu+ZO4pzS5C8bMqYX5tRYLq5am9kVwA/c/cbg9XcA3P2/dff5qqoqr66uznodIiKFzMxWu3vV2T6X\nqx79KqDSzKaYWSlwJ/Bijn6WiIicQU5aN+5+wsy+AfwbUAw84e4bcvGzRETkzHLWo3f3l4GXc/X9\nRUQkHM2MFREpcAp6EZECp6AXESlwCnoRkQKnoBcRKXA5mTAVuQizJqAug28xCmjOUjlxKpT9AO1L\nvtK+5Kee7stkdy8/24fyIugzZWbVYWaH5btC2Q/QvuQr7Ut+yvW+qHUjIlLgFPQiIgWuUIJ+QdwF\nZEmh7AdoX/KV9iU/5XRfCqJHLyIip1coZ/QiInIaiQ76QroBuZntMLO3zewtM0vU4vxm9oSZNZrZ\n+i7bysxsiZnVBo8j4qwxrNPsyw/M7L3g2LxlZp+Ns8awzGyimb1qZhvNbIOZ3R9sT9SxOcN+JO64\nmNkAM1tpZmuDfflhsD2nxySxrZvgBuTvAjeQvlXhKuAud98Ya2E9ZGY7gCp3T9y4YDP7BNAG/Mrd\nZwTb/jvQ4u6PBP8Ij3D3v4uzzjBOsy8/ANrc/cdx1haVmY0Fxrp7jZkNAVYDtwJfIUHH5gz7cTsJ\nOy5mZsAgd28zs37AG8D9wBfJ4TFJ8hm9bkCeJ9z9NaDlQ5tvARYGzxeS/sXMe6fZl0Ry93p3rwme\nHwI2kb6fc6KOzRn2I3E8rS142S/4z8nxMUly0Hd3A/JEHvyAA38ws9XB/XSTboy71wfP9wJj4iwm\nC75pZuuC1k5etzq6Y2YVwCxgBQk+Nh/aD0jgcTGzYjN7C2gElrh7zo9JkoO+0Fzl7pcCNwH3BS2E\nguDp/mAye4RpjwHnAZcC9cD/iLecaMxsMPA88IC7H+z6XpKOTTf7kcjj4u6dwe/6BGCumc340PtZ\nPyZJDvr3gIldXk8ItiWSu78XPDYCvyHdmkqyhqC3erLH2hhzPT3m7g3BL2cK+DkJOjZBH/h54Cl3\nfyHYnLhj091+JPm4ALj7AeBV4DPk+JgkOegL5gbkZjYouMiEmQ0CPg2sP/NX5b0XgXnB83nA4hhr\nycjJX8DAbSTk2AQX/h4HNrn7o13eStSxOd1+JPG4mFm5mQ0Pnp9DejDJO+T4mCR21A1AMJzqJ7x/\nA/KHYy6pR8zsPNJn8ZC+j+/TSdoXM3sGuJb0CnwNwEPA/wWeAyaRXpn0dnfP+4ucp9mXa0m3BxzY\nAfxll35q3jKzq4DXgbeBVLD5u6T724k5NmfYj7tI2HExs0tIX2wtJn2i/Zy7/72ZjSSHxyTRQS8i\nImeX5NaNiIiEoKAXESlwCnoRkQKnoBcRKXAKehGRAqegFxEpcAp6EZECp6AXESlw/x8+TjfYCcop\npQAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "bin_width = 0.1\n", "wind_10m = np.arange(0,100,bin_width)\n", " \n", "P = power_curve(effective_wind_speed(wind_10m))\n", "\n", "ind1 = wind_10m < 30. # to limit xaxis on plots\n", "plt.plot(wind_10m[ind1], P[ind1])" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "collapsed": true }, "outputs": [], "source": [ "def rayleigh_dist_power_curve(wind_10m, P, height=10):\n", " \n", " if height == 10: adj = 1.\n", " if height == 60: adj = 1.29\n", " \n", " Rd = np.zeros(wind_10m.shape)\n", " \n", " for j in range(0, len(wind_10m)):\n", "\n", " mu = wind_10m[j] * adj\n", " si = mu / np.sqrt(pi/2.)\n", " x = wind_10m\n", "\n", " if mu == 0:\n", " F = 0.\n", " F_norm = 0.\n", " else:\n", " F = (x/si) * np.exp( (-x**2.) / (2.*si**2) )\n", " F_norm = F / np.sum(F)\n", "\n", " Rd[j] = np.sum(F_norm*P)\n", " \n", " return Rd\n" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "[]" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXoAAAD8CAYAAAB5Pm/hAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XlclWX6+PHPzb4JKG4ouIGy4wamKS65JGraqtXYMlNT\nv6aaGpupZprKaZn8zrespu3bNpPT1KhZbim5ouLOIrKIuOSGioiiiKhs9++PBxlLlHM4B845cL1f\nr/Pi8JxnuY5HLm6u516U1hohhBAtl5OtAxBCCNG0JNELIUQLJ4leCCFaOEn0QgjRwkmiF0KIFk4S\nvRBCtHCS6IUQooWTRC+EEC2cJHohhGjhXGwdAED79u11jx49bB2GEEI4lPT09GKtdYeG9rOLRN+j\nRw/S0tJsHYYQQjgUpdQhU/aT0o0QQrRwkuiFEKKFk0QvhBAtnCR6IYRo4STRCyFEC9dgoldK/UMp\nVaSUyrliWzul1Cql1N7ar22veO2PSql9Sql8pdTNTRW4EEII05jSov8CGP+zbc8Da7TWvYE1td+j\nlIoE7gaiao/5UCnlbLVohRBCmK3BfvRa6w1KqR4/2zwFGFn7fA6wDniudvtcrfUl4IBSah8wCNhi\nnXCF+K+5c+eya9cuW4chbGT06NGMGDHC1mE4hMYOmOqktT5e+7wQ6FT7vCuw9Yr9Cmq3XUUp9Qjw\nCEC3bt0aGYZorUpLS/nFL35BTU0NSilbhyOamdaaNWvWsGnTJluH4hAsHhmrtdZKKbNXGNdafwJ8\nAhAXFycrlAuzpKamUlNTw4oVKxg3bpytwxHNbPz48ZSUlNg6DIfR2F43J5RSgQC1X4tqtx8Fgq/Y\nL6h2mxBWtWXLFpRS3HDDDbYORdiAUgqtpX1oqsYm+iXAA7XPHwAWX7H9bqWUu1KqJ9Ab2G5ZiEJc\nbcuWLURERODn52frUIQNSKI3jyndK/+DcTM1TClVoJR6CJgFjFVK7QXG1H6P1joXmA/sAn4AHtda\nVzdV8KJ10lqzdetWhgwZYutQhI04OTlJojeDKb1u7rnGS6Ovsf/rwOuWBCXE9ezdu5fTp09Lom/F\nlFLU1NTYOgyHISNjhcPZssXorTt48GAbRyJsRUo35pFELxzOli1b8PPzIyIiwtahCBuR0o15JNEL\nh7N161ZuuOEGnJzkv29rJaUb88hPinAo586dIzs7W+rzrZyUbswjiV44lMsDpaQ+37pJ6cY8kuiF\nQ9m61ZhhQwZKtW5SujGPJHrhUC4PlGrbtm3DO4sWS0o35pFELxzG5YFSUrYRUroxjyR64TD27dtH\ncXGx3IgVUroxkyR64TAu1+cl0Qsp3ZhHEr1wGFu2bMHX15fIyEhbhyJsTEo35pFELxzGli1bGDRo\nkAyUElK6MZP8xAiHUFZWRlZWlpRtBCClG3NJohcOYevWrdTU1DBs2DBbhyLsgJRuzCOJXjiElJQU\nnJycpGulAKR0Yy5J9MIhpKSk0K9fP3x9fW0dirADUroxjyR6YfcqKyvZunWrlG1EHSndmEcSvbB7\nGRkZXLhwgYSEBFuHIuyElG7MI4le2L2NGzcCSIte1JHSjXkk0Qu7l5KSQmhoKJ07d7Z1KMJOSOnG\nPJLohV2rqalh48aNUrYRPyGlG/NIohd2LT8/n1OnTknZRvyElG7MI4le2LWUlBQAadGLn5BEbx5J\n9MKupaSk0LFjR0JDQ20dirAjUqM3jyR6Ydcu1+eVUrYORdgRqdGbRxK9sFsFBQUcPHhQyjbiKlK6\nMY8kemG3pP+8uBYp3ZhHEr2wWykpKfj4+NC3b19bhyLsjJRuzCOJXtitlJQUhgwZgouLi61DEXZG\nSjfmkUQv7FJJSQk5OTlSnxf1ktKNeSxK9Eqp3ymlcpVSOUqp/yilPJRS7ZRSq5RSe2u/trVWsKL1\n2LBhA1prhg8fbutQhB2S0o15Gp3olVJdgd8CcVrraMAZuBt4Hlijte4NrKn9XgizJCcn4+HhIQuN\niHpJ6cY8lpZuXABPpZQL4AUcA6YAc2pfnwPcauE1RCuUnJzMjTfeiLu7u61DEXZISjfmaXSi11of\nBd4EDgPHgbNa65VAJ6318drdCoFOFkcpWpXi4mKysrIYNWqUrUMRdkpKN+axpHTTFqP13hPoAngr\npaZfuY82fuXW+2tXKfWIUipNKZV28uTJxoYhWqANGzYASKIX1ySlG/NYUroZAxzQWp/UWlcC3wE3\nAieUUoEAtV+L6jtYa/2J1jpOax3XoUMHC8IQLU1ycjJeXl7Ex8fbOhRhp6R0Yx5LEv1hYLBSyksZ\nE5GMBvKAJcADtfs8ACy2LETR2iQnJzNs2DDc3NxsHYqwU1K6MU+jR6JorbcppRYAGUAVsAP4BPAB\n5iulHgIOAVOtEahoHYqKisjNzWX69OkN7yxaLZnkzjwWDTnUWr8MvPyzzZcwWvdCmG3dunWA1OfF\n9Tk5GcUIrbUkfRPIyFhhV5KTk/Hx8WHAgAG2DkXYscvJXco3ppFEL+xKcnIyCQkJuLq62joUYccu\nJ3q5IWsaSfTCbhw7doz8/Hwp24gGXVm6EQ2TRC/shtTnhamkdGMeSfTCbiQnJ+Pn50f//v1tHYqw\nc1K6MY8kemE3kpOTGT58OM7OzrYORdg5Kd2YRxK9sAtHjhxh//79UrYRJpHSjXkk0Qu7sHr1agBu\nuukmG0ciHIGUbswjiV7YhZUrV9KpUydiY2NtHYpwAJLozSOJXthcTU0Nq1evZuzYsTLKUZhEavTm\nkUQvbC4zM5Pi4mLGjRtn61CEg5AavXkk0QubW7VqFQBjxoyxcSTCUUjpxjyS6IXNrVy5ktjYWAID\nA20dinAQUroxjyR6YVPl5eVs3LiRsWPH2joU4UCkdGMeSfTCpjZs2EBFRYXU54VZpHRjHkn0wqZW\nrlyJu7s7CQkJtg5FOBAp3ZhHEr2wqZUrV5KQkICnp6etQxEOREo35pFEL2zm2LFj5ObmStlGmE1K\nN+aRRC9s5nK3Skn0wlxSujGPJHphM6tWraJTp07ExMTYOhThYKR0Yx5J9MImampqWLVqFWPGjKlr\nnQlhKindmEd+woRNZGVlUVRUJGUb0ShSujGPJHphE0lJSYDU50XjSOnGPJLohU0sX76cgQMH0rlz\nZ1uHIhyQlG7MI4leNLvTp0+zefNmJkyYYOtQhIOS0o15JNGLZrdy5Upqamok0YtGk9KNeSTRi2a3\nfPly2rdvT3x8vK1DEQ5KSjfmkUQvmlV1dTVJSUmMHz8eZ2dnW4cjHJSUbswjiV40q7S0NIqLi6Vs\nIywipRvzSKIXzWr58uU4OTlJt0phESndmMfF1gGI1mX58uWMHzmQgJrd8ONeKDsAFWegqhQqy0A5\ng7M7OHuAewfw7GI82vQ2Hs5uFl2/pqaGM2fOcPbsWcrKyigrK+P8+fNcvHiR6urqugeAq6tr3cPd\n3R0fH5+6h6+vLz4+PrKYuY1I6cY8FiV6pZQ/8BkQDWjgV0A+MA/oARwEpmqtSyyKUji2Cyfg2DIu\nHFrBvF+k0asjsGpY7YsKXH2Nh4sP6GqouQTVF+FSsfH9ZcrZSPYBg6DDUOgwDHzDQdX/h2l1dTWF\nhYUUFBRw9OhRTpw4walTp+oS+WXOzs54eHjg7Oxc9wCorKz8yePnXF1dadeuHQEBAQQEBBAYGEhg\nYCB+fn7yC6CJSenGPJa26N8FftBa36mUcgO8gD8Ba7TWs5RSzwPPA89ZeB3haCrL4PA8ODQXTqwF\nXQO6DTsOgnvMDLpGjDGStnd3cHKt/xw11XDpJJQXQGk+lO6CMzlw/Ac48C9jH7d20Hk0dL0Fukzg\n7EUX9u7dy759+zhw4AAVFRUA+Pj4EBgYSEhICB06dMDf3582bdrg7e2Nh4dHg4m5urqa8+fPc+7c\nOcrKyjh79iynT5/m9OnTFBYWkpeXV9e69PT0JDAwkC5dutCtWzeCg4Px8PCw1r+sQEo35mp0oldK\n+QHDgQcBtNYVQIVSagowsna3OcA6JNG3HheOQ/57sPcjqDwDPqEQ+Sfodhf3/PJFMjJ2cOjtN8GU\nFq+TM3h2Nh4Bcf/drjWc2wcnN8LJDdQcXY7T4W+o0YozF4MpLovgnBpKbOwAevToQVBQEL6+vha1\nsp2dnfH19cXX17fe1ysrKykqKuLYsWMcP36c48ePs3nzZjZu3IhSis6dO9O9e3d69OhBjx49cHd3\nb3QsQhK9uSxp0fcETgL/VEr1BdKBp4BOWuvjtfsUAp3qO1gp9QjwCEC3bt0sCEPYhYoSyH4F9n4I\nNZUQfDuEz4D2Q0ApysvLWblyFb/61a8sL2sohW4TSsFZD1KP+ZG7qwedXY/Sr2MBEb57Ge+5Aq1W\noTzHgNd08Opu2i8WC7i6utK1a1e6du1at62yspKCggIOHjzIoUOHSE1NZevWrTg5OdG9e3dCQ0Pp\n3bs37du3l1KPmaRGbx5LEr0LMAB4Umu9TSn1LkaZpo7WWiul6v0ktNafAJ8AxMXFyaflqGqqYP+n\nkPUiXDoNIb+CyOehTehPdluzZg0XLlxgypQpFl1Oa01+fj4pKSkcO3YMd3d34uIGMWDAADp1qm1T\nnM1DHfwKDn4FW+4HlzbQ8z7o/Rj4R1t0fXO4urrSs2dPevbsCUBVVRUFBQV1paVVq1axatUq/P39\n6d27N+Hh4fTo0UOmbTaB1OjNY0miLwAKtNbbar9fgJHoTyilArXWx5VSgUCRpUEKO1W6F7bcB6e2\nQccRMPAdaNuv3l0XL16Mr68vI0aMaNSltNbs3r2b9evXc+LECdq2bcuECRPo27cvbm4/64njFwF9\nX4PYV+HkJuMX0f7Pjb82OgwzEn7wHUbvnmbk4uJSV7oZO3YsZ8+erUv6mZmZpKam4unpSXh4OJGR\nkfTs2VMGlV2DlG7M0+hEr7UuVEodUUqFaa3zgdHArtrHA8Cs2q+LrRKpsB9aG8kz/Xfg5AY3fgXd\n77lmeaSmpoalS5eSmJh4dVI2wbFjx1ixYgWHDx8mICCAW2+9lZiYmIZbvkpBx2HGY8Bs+PEL2Pt/\nsPkX4DED+jxpJH33dmbHZA1+fn7ExcURFxdHZWUl+/fvZ9euXeTm5rJjxw48PDwIDw8nKiqKXr16\nSUv/ClK6MY+lvW6eBL6q7XHzI/BLjEFY85VSDwGHgKkWXkPYk6py2PpLODwfOo2GIV+AV9B1D9m2\nbRtFRUVml23Ky8tZtWoVmZmZeHl5MWnSJPr379+4hOceABHPQPjvoHA17H4bsv4MuX+FkIeM7T49\nzT+vlbi6uhIeHk54eDhVVVXs37+fvLw88vLyyMzMxMfHh+joaGJjY+ncuXOrr+lL6cY8FiV6rXUm\nEFfPS6MtOa+wU+UFsH4KlOyAvm9A5LPX7MN+pcWLF+Pi4kJiYqJJl9Fak5ubS1JSEhcvXuTGG28k\nISHBOl0UlRMEjjMeZ7Ih7y3Y93+w9wPoNhWiXmjWOn59XFxcCAsLIywsjKqqKvbu3UtWVhbbt29n\n69atdOjQgdjYWGJiYvDz87NprLYipRvzyMhYYZpTabBhMlSegxFLoOskkw9dsmQJI0aMwN/fv8F9\nz58/z9KlS8nPz6dLly5Mnjz5vzdZrc0/xviLpO/rkP+u0SX00FwIvhOi/wxt+zbNdc3g4uJCREQE\nERERlJeXk5ubS1ZWFmvWrGHNmjX06tWL/v37Ex4ejotL6/lxltKNeVrP/wzReEUpsG6iUf4Yt9Ks\nFu/evXvJy8vjsccea3DfH3/8kYULF3LhwgXGjh3L4MGDm6cu7dUV+v8NIp8zEn7+u3BkAQTdCtEv\nQrsBTR+DCby8vIiPjyc+Pp7Tp0+TlZVFZmYm3377LZ6ensTGxjJgwAA6duxo61CbnJRuzCOJXlxf\n4WqjXOMdDDetMZKiGRYtWgTA5MmTr7lPTU0Na9euZdOmTbRv357p06c3XSv+etwDIPYVo16f/3fY\n/Q4ULDJG3ca+ahct/MvatWvHyJEjGTFiBD/++CMZGRmkpqaybds2goKCGDBgAFFRUY26+e0IpHRj\nHkn04tqOJcGG28C3D4xaBZ7mJ98FCxYwcOBAunfvXu/r5eXlLFiwgAMHDtC/f38SExNxdb3GlAjN\nxa0txLwMYU/DnveMOn5SP+h+N8T8xfj3sBNKKUJCQggJCeH8+fNkZWWRkZHBkiVL+OGHH4iOjmbA\ngAE/GcjVEkjpxjyS6EX9Tm6ClNvBLxJuWmW0ds105MgRtm/fzhtvvFHv60VFRcydO5fS0lKmTJlC\nv37198G3GTc/o1bf5wnIe9PoqXP4G+j1S4h+yfgrx454e3szZMgQBg8eTEFBARkZGWRnZ5ORkUGX\nLl2Ii4sjOjra9r9IrUBKN+aRRC+udiYX1k0Cr2AYtaJRSR7gu+++A+COO+646rW9e/fyzTff4O7u\nzoMPPkhQ0PW7aNqUm78xAKvPk0Z3zH3/Bwe+hN6/gag/gkcHW0f4E0opgoODCQ4OZvz48WRlZZGa\nmsqSJUtYuXIl/fr1Iy4ujoCAxn2u9kBKN+aRRC9+6vxhSL4ZXDxh1EqLktiCBQuIiYmhd+/eP9me\nmZnJkiVL6NSpE/feey9t2rSxNOrm4dkJ4t6FiBnGvD573jUGjoX/DiJ+b0y1bGfc3d2Jj48nLi6O\nw4cPk5qaWtdNMyQkhLi4OPr06eNwg7GkdGMeSfTivyrLjN41VWUwZgP49Gj0qQoLC9m0aRMzZ86s\n26a1ZuPGjaxdu5ZevXoxdepUx5zF0bs7DP4cIv4A2S9DzqtG18zolyD0UYsXR2kKSim6d+9O9+7d\nKSsrIyMjg/T0dObNm4evry8DBw5kwIAB+Pj42DpUk0jpxjyS6IVB18DWB4w530etgLaxFp1u4cKF\naK3ryjZaa1atWsWWLVuIjY1l8uTJjj+Pi184DJsHp/4Amc9C+m+Nrpl9/wrd7mryGTMby8fHh+HD\nhzNs2DD27NlDamoqycnJrF+/noiICOLj4+nWrZtdj76V0o15JNELQ+5f4ch30P8t6DzG4tMtWLCA\nsLAwIiMj0Vrzww8/sH37duLj40lMTLTrJGK2gDij6+nxH2DHs7BpGux+C/r9DTo1bhK35uDk5FQ3\n7cKpU6dIS0sjMzOT3NxcOnbsSHx8PLGxsXbZRVNKN+ZxrMKcaBoFS41phntMN+rNFiouLmb9+vXc\neeedgLFO7Pbt2xk8eHDLS/KXKQVdEiExEwb/Ey4cgzUjYd0txs1tOxcQEMDNN9/MjBkz6v7aWrZs\nGbNnzyYpKYni4mJbh/gTUroxj7ToW7vzh4w529sNhEGfWKXcsGjRIqqrq7n99ttJSkoiLS2NoUOH\nMnr06JaZ5K/k5Ay9HoRu02DP3yH3DUiKNbpkxvzF7AFnzc3V1ZX+/fvTr18/jh49yvbt20lPT2f7\n9u307NmT+Ph4wsLCbH7zVko35pFE35rVVBlT9upqGDbf6GljBXPnziU0NJSzZ8+SmprK4MGDW0eS\nv5KLpzGlQsjDkPM67H0fDn5d20PnWaOPvh1TShEUFERQUBA333wzGRkZpKWlMX/+/LqbtwMHDsTb\n29sm8UnpxjxSumnNcl4xBkYN+hh8elnllIWFhSQnJ3PvvfeyYcMG+vXrx7hx41pXkr+SewAMnA2T\n8iHoNuNeyNJQY4qF6gpbR2cSb29vEhISeOqpp5g2bRrt27cnOTmZ2bNn891333HkyJFmT7hSujGP\ntOhbqxPrIOc1o8zQ4x6rnfabb74hNjYWJycnIiIiuOWWW1pvkr+ST08Y+pXRB3/Hs5D+1BU9dKba\nbQ+dK11587a4uJjU1FR27txJdnY2nTt3Jj4+npiYmGYZeSulG/NIi741qjhr1OXbhMLA96x66jVr\n1jBlyhR69erF7bffbvNart1pNxBuWg0jk8DFBzbdDSsGQeFaW0dmlvbt25OYmMiMGTOYOHFi3Spi\ns2fPZsWKFZw+fbpJry+lG/NIi7412vEHuHAUxm4GV+sNkMnIyCAqKgonJyemTZvWquZHN4tS0GU8\ndB5rLGCe9SKsHQ2BN0O/Wddcd9ceubm5ERcXx8CBA68aeRsaGkp8fDyhoaFW/4UvpRvzyE9ia1O4\n2hi2H/EHaH+D1U5bWlrK0qVLuXDhAvfee69d9r22O07O0Ot+6D4V9nwIua9DUn/ofi/0fdVq902a\nw5Ujb8+dO0d6ejrp6en85z//wd/fn7i4OPr374+Xl5fVrgfSojeVJPrWpPIcbH0IfMOMrn5WUlFR\nwdy5c6mqqiI7O5uoqCirnbtVcPYwavchv4Jdf4P8d+DINxD6mDF7pp1NmtaQNm3aMHLkSBISEti9\nezepqamsXr2a5ORkoqOjiY+Pt3jaZEn05pFE35pkPgflR2DsRqt1pdRas3DhQgoLC5k/fz6/+c1v\nrHLeVsnNH/r91ZgWOecvxjq2P/7T+Osr/HdWLbM1B2dnZ6KiooiKiqKoqKju5u3OnTvp0qUL8fHx\nREdHN6rEJzV688idstbi5BZj4q2wp6DDjVY7bXJyMrt376a6upq9e/dy1113We3crZZXF6PL64Qc\nCBwL2S8ZXTL3fAg1lbaOrlE6duzIxIkTeeaZZ0hMTKSiooLFixcze/ZsVq9ezZkzZ8w6n9TozSMt\n+tagpgpSHwPPrsZSeVaSn59PSkoK/fv35+mnn2bMmDF06dLFaudv9fzCIeFbKN5q/DWW9rix+Enf\n16HbnaAcr53m7u7OoEGDiI+P5+DBg2zfvp3NmzezadMm+vTpQ3x8PCEhIQ12yZXSjXkk0bcGez6A\nMzth2AJwtc7c76dOnWLhwoUEBgbi5+fHgQMHeOUV6/0SEVdoPxhGrzOWdsx8zpg0LS8O+r0BnUY7\nRB/8n1NK0bNnT3r27MnZs2dJT08nIyODPXv20K5dO+Li4ujXrx+envWXGKV0Yx7HaxII85QfM7rv\nBSZC8O1WOWVFRQXz58/HycmJqVOn8tVXX+Ht7c1tt91mlfOLeigFXSfUTpo2By4WwdqxsGYUFG20\ndXQW8fPz46abbuLpp5/m9ttvx9vbm5UrVzJ79myWLFlCYWHhVcdI6cY80qJv6TJmgK6EuPes0vLT\nWrNs2TKKioqYPn067u7uzJs3jzvvvNNm8560KnVdMqfBvk+NLpmrE4w++LGvQkC8rSNsNBcXF2Ji\nYoiJiaGwsJDt27eTnZ3Njh07CA4OJi4ujoiICFxdXaV0YyZJ9C1Z0QY4PA9iZkKbEKucMi0tjays\nLEaNGkVISAjz5s2jtLSU+++/3yrnFyZydoewJ4wumXs/hF2zjBG2QVMg5hWLF46xtc6dOzN58mTG\njh1LZmYmaWlpLFy4kKSkJGJjY2nbti0gid5Uyh7+oeLi4nRaWpqtw2hZaqphRTxcKjYm1LJCd8oT\nJ07w6aef0qtXL+655x6UUkycOJHs7GwOHjwo0x3YUuU5Y+6cvDeh8qwxTXLMTOOGbgugtebgwYNk\nZGSQl5dHdXU1BQUFDBgwgIceeqjVDtBTSqVrreMa2k9a9C3VgX9ByQ648WurJPmKigoWLFiAp6cn\nU6ZMQSlFYWEhK1as4Nlnn5Ukb2uubYzBVX0eh7zZ/x101WM6xLzsUKNs63Plzdvy8nJWr15NUVER\nRUVFvPXWW0RHRzNgwAC6dOkik+jVQ346W6LKc7DzTxAwGLrfbZVTrlixguLiYm677ba6WvzXX39N\ndXU19913n1WuIazAra0xfcLkHyF8BhyeD0vDYNsjUHbA1tFZhZeXF5GRkXz44Yd07dqVyMhIsrOz\n+eyzz/j444/Zvn07Fy9etHWYdkUSfUu0axZcLISB71jlBmxubi4ZGRkMHTqUXr2MlqHWms8//5wb\nbriBiIgIi68hrMyjA/T/XyPh9/5/cGAOLO0NWx6E0j22js5il1vtbm5uTJkypW4WTScnJ5KSknjr\nrbdYuHAhhw4dkjo+VijdKKWcgTTgqNZ6klKqHTAP6AEcBKZqrUssvY4w0fkjkPeWMTGWFSYtO3Pm\nDEuXLqVr166MGjWqbvvWrVvZtWsXn332mcXXEE3IM9DocRX5R8j7X9j3MRz80qjhR70A/o45L9HP\ne914eHgQFxdHXFwcx48fJyMjg+zsbLKysmjbti19+/alb9+++Pv72zJsm7FGi/4pIO+K758H1mit\newNrar8XzSXnL4A2Rk9aqKamhoULF6K15o477sDZ2bnutU8//RQfHx+mTZtm8XVEM/DqAgPfhskH\nIPz3cHQJLI+BlLugZKetozPb9QZMBQYGMnHiRGbMmMFtt92Gv78/69at491332XOnDlkZmZSUeEY\nq3tZi0UteqVUEDAReB2YUbt5CjCy9vkcYB3wnCXXESY6u9uYBKvPk+DTw+LTbdu2jcOHD3PrrbfW\ndWcDY0riefPmce+99+Lj41gTbbV6np2g//9A5LOw+x1jAfMjC6DrZIh+EQIa7MBhF0wZMOXm5kZs\nbCyxsbGcPXu2bkK1xYsXs3z5ciIjI+nXrx/du3dv8TdwLS3dvAM8C1w5rr6T1vp47fNCoJOF1xCm\nynoBnL2MP8ktdPLkSdasWUNYWBixsT/tkz137lzKy8t5+OGHLb6OsBH3AOOmbcQzxvq1+e8Y3XED\nxxuLmnccYddTK5g7YMrPz4/hw4eTkJDAkSNH2LlzJ7m5uezcuRN/f39iY2Pp27cv7dq1a8qwbabR\niV4pNQko0lqnK6VG1reP1lorper9JJRSjwCPAHTr1q2xYYjLirfDke+MvtMWzl9eXV3NokWLcHd3\nZ9KkSVe1dj777DNiYmIYNGiQRdcRdsDNH2JegvCnjdkx8982plUIGAQRz0LQrcZoXDvT2LlulFJ0\n69aNbt26MX78eHbv3s3OnTvZsGEDGzZsoFu3bsTExBAVFXXNeXYckSUt+qHAZKXUBMAD8FVK/Rs4\noZQK1FofV0oFAkX1Hay1/gT4BIwBUxbEIbSGnc+DewejS52FNm7cyLFjx7jrrruuKs3s3LmT1NRU\n3n333Rb/526r4uoLUc8bCf/HOcbAq413QpveEPF76Hm/sUCKnbDGXDeurq51Uy6Ulpayc+dOsrKy\nWLZsGUk7MvrLAAAcIUlEQVRJSYSGhhITE0NYWFizLHjelBqd6LXWfwT+CFDbov+91nq6Uup/gQeA\nWbVfF1shTnE9havhRDIMfNfi2SmPHz/Ohg0biImJITIy8qrXP/30U9zd3Zk+fbpF1xF2ytkDej8K\nIQ9DwULY9T+w/VHIegnCfgu9HzP66tuYtee68fX1JSEhgWHDhlFYWEh2djY5OTns2bMHV1dXIiIi\niImJoVevXg45OLApRsbOAuYrpR4CDgFTm+Aa4jKtYecfwbsHhD5q0amqqqpYuHAhXl5eJCYmXvX6\nuXPn+Ne//sXUqVNbbC1T1HJyNua8D74DitYbCX/nC5D7BoT82kj6Vrjh3+jwmmiaYqUUgYGBBAYG\nMmbMGA4fPkxWVhZ5eXlkZWXh5eVFVFQUMTExBAUFOcxftVZJ9FrrdRi9a9BanwJGW+O8wgRHv4fT\n6TD4C2OiKwusW7eOkydPcu+999Zbn/zyyy85d+4cjz/+uEXXEQ5EKeg00niUZBl98ff8Hfa8a9Tv\nw56GDsOa/cZtc0xT7OTkRI8ePejRowcTJkxg3759dbNppqam4u/vX1fP79ixo10nfZnrxpFpDbmv\ngXdP6PELi051/PhxNm/eTP/+/endu3c9l9K8//77xMXFyU3Y1qptLNz4JfT9qzFj5r6PjQ4AbQcY\nS1R2n2ZxY8NUzT1NsYuLC+Hh4YSHh3Pp0iXy8vLIzs5m48aNpKSkEBAQQGRkpN0mfUn0juzEGji1\n3Vhf1KnxH2V1dTVLlizB29ubcePG1bvPunXryMvL45///Kfd/ScWzcw72FjdKvpFOPhvoz/+1gcg\n81mjhh/6/4z++k3IlitMubu7069fP/r160dZWRl5eXns2rXLrpO+JHpHlvO6sQ5szwcsOs2WLVso\nLCxk6tSpeHjU37Pigw8+oF27djISVvyXixeEPmLU7AtXG33xs2dC7l+h+z3GTJpNtBCKvaww5ePj\nQ3x8PPHx8ZSVlbF7925yc3PtLulLondUJzdB0ToY8I5Ffy6fOnWKdevWERERcc3JyQoKCli0aBEz\nZsxoUX2LhZUoBYFjjUdpPuS/Bwe+MCZSazfQaOV3v8f4xWC1S9rfClM+Pj518+1cTvo/b+lf/jkL\nDAxs1qQvid5R5bxu9JsP/XWjT6G1ZunSpbi6utbby+ayTz75hJqaGh577LFGX0u0Er5hEP++MdfS\ngS9h30ew7WHIeMb4y7P3/wM/y2c7tffFwa9M+ufPn68r72zatImNGzfi5+dHeHg4ERERBAcHN3mX\nTUn0juh0OhxPgr5vWNRKSk9P59ChQ9xyyy20aVN///sLFy7w0UcfMWnSJHr27Nnoa4lWxs3PWOqw\nz+NwMgX2fmQk/T1/h44jjVZ+0K3g3LiVoeyldGMKb2/vuqRfXl7Onj17yMvLIy0tjW3bthEWFsbd\nd1tn3YhrkUTviHL/Cq7+0Oc3jT5FaWkpq1evpmfPnvTv3/+a+3355ZcUFxfzzDPPNPpaohVTCjoO\nNx4Xi2D/P4zeOpumgUcn6PVL6PUr8L26p9f1T2t/pRtTeHl51d3IvXTpEvv27cPdvel7KjneEK/W\n7kyu0aUt7Elj2HojLV++nOrq6nrnsrmspqaG2bNnM3DgQIYPH97oawkBgEdHY5qFW/bBiGXGfDp5\nf4Pv+8Cq4fDjF1BZZtKpHDXRX8nd3Z2oqChCQ0Ob/FqS6B3NrjfAxdvot9xI+fn55OfnM2LEiOuO\ncF2+fDn5+fnMmDHDLrqIiRbCyRm6ToARS2DKEeg3Cy6egK2/hIWBsO3XcHKLMU7kWqew8xq9vZFE\n70jO7YdD/zHqm+4BjTpFRUUFSUlJdOjQgSFDhlx337feeougoCDuuuuuRl1LiAZ5dTGmRZ60G8ak\nQLe7jP/jq26EZVHG5GoXCq86zJFq9PZAEr0j2TULlCuEN75evn79es6ePcukSZN+smLUz2VkZLBu\n3Tp++9vfOvzMfcIBKAUdh8Hgf8Btx+GGz4zJ03b8ARZ1heTxcODfdaWdllC6aU5yM9ZRnD9i9EsO\neQQ8OzfqFEVFRWzdupV+/fo1uAbAW2+9hY+PD7/+deO7bwrRKK5tIOQh43F2tzH69uC/Yct9xsI6\nwbfh3HUazk6S6E0lLXpHkfe/Rs0y8tlGHa61ZtmyZbi7uzN27Njr7rtv3z7mzp3Lo48+2moXUxZ2\nwi8c+r4Gk380Sjs974Njy3HbNJmj78No/6VwKu269Xwhid4xXDgB+z81Fn/wbtxqXJmZmRw+fJgx\nY8bg5XX9vvf/8z//g6urq3SpFPZDORmlnUH/B7cdp+rGb0jJhwG+240lEL8Ph50vwplsSfr1kETv\nCPLfhpoKiHy+UYeXl5ezatUqgoODr9tnHuDIkSPMmTOHhx56iMDAwEZdT4gm5eyOCr6Nu96F2Qef\nhUGfglcw7PorLI+FZZGQ9bLRFVkAkujt36XTsOcD6DbN7EEll61evZqLFy8yceLEBrtJvvnmm2it\nefbZxpWIhGgOl/8fX6hyh9CHYfRq4yZu/EfgGWhM3708Gr6PhKyZcHaXbQO2MUn09m7Pe1BVBlF/\natThhw8fZseOHQwePJhOna4/dWxRURGffvop06dPp3v37o26nhDNod5eNx4djbl0Rq+FW49C3AfG\ntpxXjK6ay6KMVbJOpba68o70urFnlecg/10ImgL+0WYfXl1dzbJly/D19WXkyJEN7v/2229z8eJF\nnn++cSUiIZpLg90rPTsbU4T0+Q1cOA6Hv4UjC4wuyrl/Nab3DppizLfTcUSj59xxFJLo7dnej6Ci\nBKJeaNTh27Zto6ioiGnTpuHmdv3/yCdPnuS9995j6tSphIWFNep6QjQnpZRpA6Y8A40J1sKegEun\n4OgyKFgEP/7TWCnL1Q+6TDCSfpfxFk0tYq8k0durqguw+y3oPK5RizecPXuWdevW0adPH5MS96xZ\ns7hw4QIzZ85sRLBCND+llPn96N0DoNf9xqPqgrFgSsEiOLrEGJHr5AadboIuE41pGnx6NU3wzUwS\nvb3a/5kx21/0nxt1+A8//IDWmsTExAZvwB49epQPP/yQ+++/n/Dw8EZdT4jm5uTkZNmAKRdPCLrF\neNRUQ/FmKFhsPNKfNB6+YRA4AbokGjNwNtOauNYmid4eVVcYs/p1SICOCWYfvmfPHnbv3s3o0aNN\nGvD0+uuvU11dzUsvvdSYaIWwCZNLN6ZwcjZ+1jomwIA3oXSvsebDseVGeSf/bWMywU6jjTJPl8RG\nj2mxBUn09ujAv6C8AAZ9ZvahlZWVJk9aBnDgwAE+++wzHn74YVlYRDiURpVuTOXb23iE/RaqyuFE\nspH0jy03yjwAflFGabXzGKO17+rTNLFYgSR6e1NTZfQMaBcHgePMPnzDhg2cOXOGBx988LqTll32\nyiuv4OTkxAsvNO6GrxC2YnHpxlQuXtB1ovHQGkp3w7Eko8V/ubWvXKD9EOg81kj8AfHgZD/p1X4i\nEYZD86BsPyQsNGb0M8PJkyfZvHkzffv2Nakf/M6dO5kzZw4zZsyga9eujY1YCJuwaunG9Isaa976\nRUDEDOOGbvFm46Zu4WrIfhmyXwKXNtBplJH0O48B33Czf56tSRK9PdE1xjBuv2gImmzeoVqzfPly\n3NzcGpy07PL+v//972nbtq205oVDatLSjalcPKHzaOPBG8ZI9hPJ/038l8s8Hp2M/vqdRhpffSOa\nNfFLorcnBYuNodo3fm1M4mSGrKwsDh48yKRJk/D29m5w/6SkJFavXs0777xD27ZtGxuxEDbTbKUb\nc7i3g253GA+AsoNwYg2cWGf8Ajg839ju0RE6DDcSf6fRxiydTUgSvb3QGnJeA59Q6DbVrEMvXLjA\nypUrCQoKYsCAAQ3uX1VVxR/+8AdCQ0N57LHHGhuxEDZlk9KNuXx6gE/t3PpaQ9mPULTeSPxF64zR\nusF3QsI3TRqGJHp7cXwFlGTADZ8bXb3MsGbNGi5cuGDSpGUAn3/+Obt27eK7775rcMSsEPbKLko3\n5lAK2oQYj5BfGYn//EGovtjkl5ZEbw+0Nmbb8wqGHtPNOvTo0aOkp6dzww030LlzwytPlZSU8OKL\nL5KQkMCtt97a2IiFsDm7LN2YQynwaZ4uzZLo7UHRBji5CeLeN2typZqaGr7//nvatGnDqFGjTDrm\nz3/+M6dOneK9994zqfUvhL1yiNKNnWj0NMVKqWClVLJSapdSKlcp9VTt9nZKqVVKqb21X+VOX0Ny\nXzfuyvf6lVmHpaamUlhYyM0334y7e8NDs9PT0/noo4944okn6Nu3b2OjFcIuOFzpxoYsmY++CnhG\nax0JDAYeV0pFAs8Da7TWvYE1td+LayneDoWrIPwZo6uWic6dO8fatWsJCQkhMjKywf1ramr4zW9+\nQ8eOHXnllVcsiVgIuyCJ3nSNLt1orY8Dx2ufn1NK5QFdgSnAyNrd5gDrgOcsirIly3kV3NoZCyaY\nYcWKFVRXVzNhwgSTb8Bu376dL7/8Ej8/v8ZGK4TdcPgafTOyygpTSqkeQH9gG9Cp9pcAQCFQ77JG\nSqlHlFJpSqm0kydPWiMMx3M6HY59DxHPgGsbkw/bt28fubm5JCQk0K5duwb3Lyws5LnnnmP48OH8\n4he/sCRiIeyG1OhNZ3GiV0r5AN8CT2utS698TRu/buv9lau1/kRrHae1juvQoYOlYTimnFfB1R/6\nPGHyIRUVFXz//fe0b9+eoUOHmnTME088QXl5OR9//LHcgBUthpRuTGdRoldKuWIk+a+01t/Vbj6h\nlAqsfT0QKLIsxBaqJNMYCRv+O7NWtFm7di1nz55l8uTJuLg0XHlbsGAB3377LX/5y19krnnRokjp\nxnSW9LpRwOdAntZ69hUvLQEeqH3+ALC48eG1YDmvGQk+7LcmH1JQUMC2bduIj48nODi4wf2Li4t5\n/PHHGThwIM8884wl0Qphd6R0YzpL+tEPBe4DspVSmbXb/gTMAuYrpR4CDgHmjedvDc7kwJFvIfpF\ncGt4YRAwFvpeunQpvr6+jB492qRjnnrqKUpKSli9erVJrX8hHImUbkxnSa+bjcC1Cr6mZaLWKuc1\ncPGBsKdNPmTjxo0UFRVxzz33mNRnft68eXz99dfMnDmTmJgYS6IVwi5J6cZ0Vul1I8xwNs+Ywa7P\nk8ZMdyY4efIkKSkpREdH06dPnwb3P3ToEI8++iiDBw+WKYhFiyWlG9NJom9u2TONFWvCZ5i0u9aa\npUuX4ubmxvjx4xvcv7q6mvvuu4+amhq++uorKdmIFktKN6aTRN+cTmcYrfnwGeDR3qRDtm3bxpEj\nRxg3bpxJ88y/8cYbpKSk8MEHH9CrVy9LIxbCbknpxnSS6JvTzj8bo2DDTesBc/LkSdasWUOfPn1M\nmptm/fr1zJw5k7vvvpvp082bBVMIRyOlG9NJom8uRSnGYsKRz4Nbw1MQVFdXs2jRIlxdXbnlllsa\nHOh09OhRpk6dSkhIiAyMEq2ClG5MJwXc5qA17PwTeAZCn8dNOmTjxo0cO3aMu+66Cx8fn+vuW1FR\nwdSpUzl//jxr167F19f0AVhCOCop3ZhOEn1zOP4DnNwI8R8aN2IbcOzYMTZs2EBsbKxJM1M+88wz\nbN68mXnz5hEVFWWNiIWwe1K6MZ2UbppaTRXs+D34hECvhxrcvbKykoULF+Lt7U1iYmKD+3/88ce8\n//77/O53v2PqVBmbJloPKd2YThJ9U9v3CZzdBf3fNGn1qJUrV1JcXMyUKVPw8PC47r4rVqzg8ccf\nJzExkb/97W/WilgIhyClG9NJom9KFWcg+yXoOBKCpjS4e25uLmlpadx4442EhIRcd9+cnBzuuusu\noqKimDdvnvSXF62OlG5MJ4m+KeW8CpdOw8C3jYWAr6OkpISlS5cSFBTETTfddN19CwoKmDhxIm3a\ntGHZsmW0aWP6XPZCtBRSujGdNAObSule2PMehPwK2va77q7V1dUsWLAApRR33HEHzs7O19y3qKiI\nMWPGUFJSwvr16wkKCrJ25EI4BCndmE5a9E1Ba0h/EpzcIfa1BndPSkri2LFjTJ48GX//a89mWVJS\nwrhx4zh8+DDLli2jf//+1oxaCIcipRvTSYu+KRz6DxxfAQP/Dp6dr7treno66enpDB06lIiIiGvu\nV1paysSJE8nLy2PJkiUkJCRYO2ohHIqUbkwnid7aLp2C9Kch4Abo/Zvr7nr48GGWL19OaGjodevy\np0+fZvz48ezYsYP58+dz8803WztqIRyOlG5MJ4ne2nb8ASpKYNAn4HTtWntpaSnz58/H39+f22+/\nHSen+qtoJ06cYOzYseTn5/Ptt98yefLkpopcCIcipRvTSY3emgrXwo//hIg/QNvYa+528eJFvvrq\nKyorK5k2bRqenp717nfw4EGGDx/O/v37WbZsmSR5Ia4gpRvTSaK3loozsPWX4BNqLBF4DVVVVcyb\nN4/i4mKmTZtGx44d691v27Zt3HDDDRQVFbFixQrGjBnTVJEL4ZAk0ZtOEr21pD4OF47Cjf8Gl/pb\n6FprFi1axMGDB7n11luvOV/8N998w8iRI/Hx8WHLli0MGzasKSMXwiFJjd50kuit4cC/4dDXEDMT\n2t9Q7y5aa77//ntyc3MZO3Zsveu4VlVV8cILLzB16lQGDhzItm3bCA8Pb+LghXBMUqM3ndyMtdSZ\nbNj+KHRIgMg/1rvL5eUAd+zYwbBhwxgyZMhV+xQWFnLPPfewbt06HnroId5///0G57oRojWT0o3p\nJNFbouIMbLjNWEhk2Lx6e9lcmeQTEhIYNWrUVYuC/PDDDzz44IOUlpbyxRdf8MADDzTXOxDCYUnp\nxnRSummsmkrYdDecPwTDFhiLivxMVVUV3333HTt27GD48OFXJfnS0lJ+/etfk5iYSEBAANu2bZMk\nL4SJpHRjOkn0jaE1bH/EGP0a/xF0uPGqXcrLy/nyyy/Jyclh9OjRjBw58idJ/vvvvycmJoZ//OMf\nPPfcc6Snp9dbtxdC1E9KN6aT0o25tIadL8CPX0D0SxD68FW7nD59mq+//pozZ85wxx13EB0dXffa\n3r17efrpp1m+fDkRERFs3Lix3pq9EOL6pHRjOkn05tAasv4Mu96AkF8bvWx+Jicnh6VLl+Ls7Mz9\n999Pt27dAGOE66xZs/jwww9xd3fnzTff5Mknn8TNreHFSIQQV5PSjekk0ZtK1xjTG+yeDaGPGCWb\nK0oxlZWV/PDDD2RkZBAcHMwdd9yBn58f5eXlvPrqq/z973/n0qVLPPDAA7z22msEBl5d0xdCmE4S\nvekk0Zuiqhy23AdHvoM+T8DAd0H99/bG5SkKSkpKGDp0KKNGjaqbU/7tt99m1qxZ3Hvvvbz88sv0\n6dPHVu9CiBbFycmJqqoqW4fhECTRN6Q03+hdU7ITBrwNYU/VteTPnTvHypUrycnJoV27dtx///30\n7NnzJ4dnZmbSu3dvvvrqK1tEL0SLJS1600mivxZdYyzsnfGMMaXBiKXQdSIAZWVlbNq0ibS0NLTW\njBgxgmHDhtW7buvu3btldKsQTUB63ZiuyRK9Umo88C7gDHymtZ7VVNeyupNbIP0pOJ0KnUbDkDng\n1ZUTJ06QkZHBjh07qKqqIjY2luHDh9OuXbt6T1NdXc2ePXuYMGFCM78BIVo+6XVjuiZJ9EopZ+AD\nYCxQAKQqpZZorXc1xfWsoqYSjq+C3W/BibXg0RmG/JvTvuPZm72P7Owkjh49irOzM1FRUQwfPpyA\ngIDrnvLAgQNUVFRIi16IJiClG9M1VYt+ELBPa/0jgFJqLjAFsJ9EX3Ueyg7C6TQ4mYI+shhVUUyV\na0eOtnuCnItD2b+0iJKS9wHo0KEDN998M7GxsXh5eZl0id27dwNIoheiCUjpxnRNlei7Akeu+L4A\nqH9aRwtkrJ1DlwMzcFIahUYpjROXn1O3TaF/tk8N7s6Vdee5UO3B/vIQcs6NZl95b6q1Cxcv7uTU\nqVOcOHGCwsJCysrKzI6vpKQEgLCwMKu9ZyGEwcnJiaysLKKiomwdikUSExN58803m/QaNrsZq5R6\nBHgEqBtUZC5377YUXghAo6jRykjx+ornKLTG2Fb72uXt5yo9KLnky9Hy9hw/709lVTXV1dVUVBzh\n4sWLVFdXA+Dm5tbo+MBozV+rhi+EaLxHHnnE5L+u7VnXrl2b/BqqKf70UUoNAWZqrW+u/f6PAFrr\nN+rbPy4uTqelpVk9DiGEaMmUUula67iG9muqSc1Sgd5KqZ5KKTfgbmBJE11LCCHEdTRJ6UZrXaWU\negJYgdG98h9a69ymuJYQQojra7IavdZ6ObC8qc4vhBDCNDIfvRBCtHCS6IUQooWTRC+EEC2cJHoh\nhGjhJNELIUQL1yQDpswOQqmTwCELTtEeKLZSOLbUUt4HyHuxV/Je7FNj30t3rXWHhnayi0RvKaVU\nmimjw+xdS3kfIO/FXsl7sU9N/V6kdCOEEC2cJHohhGjhWkqi/8TWAVhJS3kfIO/FXsl7sU9N+l5a\nRI1eCCHEtbWUFr0QQohrcOhEr5Qar5TKV0rtU0o9b+t4LKGUOqiUylZKZSqlHGpyfqXUP5RSRUqp\nnCu2tVNKrVJK7a392taWMZrqGu9lplLqaO1nk6mUcojV3pVSwUqpZKXULqVUrlLqqdrtDvXZXOd9\nONznopTyUEptV0rtrH0vf6nd3qSficOWbmoXIN/DFQuQA/fY9QLk16GUOgjEaa0drl+wUmo4UAb8\nS2sdXbvtb8BprfWs2l/CbbXWz9kyTlNc473MBMq01k273puVKaUCgUCtdYZSqg2QDtwKPIgDfTbX\neR9TcbDPRSmlAG+tdZlSyhXYCDwF3E4TfiaO3KKvW4Bca10BXF6AXDQzrfUG4PTPNk8B5tQ+n4Px\ng2n3rvFeHJLW+rjWOqP2+TkgD2M9Z4f6bK7zPhyONlxegNq19qFp4s/EkRN9fQuQO+SHX0sDq5VS\n6bXr6Tq6Tlrr47XPC4FOtgzGCp5USmXVlnbsutRRH6VUD6A/sA0H/mx+9j7AAT8XpZSzUioTKAJW\naa2b/DNx5ETf0gzTWvcDEoHHa0sILYI26oOOWSM0fAT0AvoBx4G3bBuOeZRSPsC3wNNa69IrX3Ok\nz6ae9+GQn4vWurr2Zz0IGKSUiv7Z61b/TBw50R8Fgq/4Pqh2m0PSWh+t/VoELMQoTTmyE7W11cs1\n1iIbx9NoWusTtT+cNcCnONBnU1sH/hb4Smv9Xe1mh/ts6nsfjvy5AGitzwDJwHia+DNx5ETfYhYg\nV0p5195kQinlDYwDcq5/lN1bAjxQ+/wBYLENY7HI5R/AWrfhIJ9N7Y2/z4E8rfXsK15yqM/mWu/D\nET8XpVQHpZR/7XNPjM4ku2niz8Rhe90A1Haneof/LkD+uo1DahSlVC+MVjwY6/h+7UjvRSn1H2Ak\nxgx8J4CXgUXAfKAbxsykU7XWdn+T8xrvZSRGeUADB4FHr6in2i2l1DAgBcgGamo3/wmjvu0wn811\n3sc9ONjnopSKxbjZ6ozR0J6vtX5FKRVAE34mDp3ohRBCNMyRSzdCCCFMIIleCCFaOEn0QgjRwkmi\nF0KIFk4SvRBCtHCS6IUQooWTRC+EEC2cJHohhGjh/j+WiW3/m1TVdwAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# My attempt at reproducing Figure 2 in Macleod et al., 2017\n", "\n", "plt.plot(wind_10m[ind1], P[ind1], color='k')\n", "\n", "Rd10 = rayleigh_dist_power_curve(wind_10m, P, height=10)\n", "plt.plot(wind_10m[ind1], Rd10[ind1], color='Grey')\n", "\n", "Rd60 = rayleigh_dist_power_curve(wind_10m, P, height=60)\n", "plt.plot(wind_10m[ind1], Rd60[ind1], color='orange')" ] }, { "cell_type": "code", "execution_count": 16, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 16, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjAAAAF8CAYAAADRg4HlAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3Xd8leX9//HXlb0ngYQACWBR3ANRsS60jmodWHHPit9u\n26oV2/r7flstrdaBHdpaB7aoBS21YEVFhghlCOIAFWUmgYTsvXOu3x93ziFRRsY5ObnPeT8fj/uR\n5OTk3B8UkneuzzWMtRYRERERN4kIdgEiIiIivaUAIyIiIq6jACMiIiKuowAjIiIirqMAIyIiIq6j\nACMiIiKuowAjIiIirqMAIyIiIq6jACMiIiKuExXsAvpjyJAhNj8/P9hliIiIiJ+sX7++3FqbdbDn\nuTrA5Ofns27dumCXISIiIn5ijNnZk+ephSQiIiKuowAjIiIirqMAIyIiIq6jACMiIiKuowAjIiIi\nrqMAIyIiIq6jACMiIiKuowAjIiIirqMAIyIiIq6jACMiIiKuE7AAY4x5xhhTaozZ2OWxDGPMImPM\n551v07t87h5jzBZjzGZjzHmBqktERETcL5BnIc0C/gj8rctj04HF1trfGmOmd358tzHmcOAq4Ahg\nOPCWMWactbYjgPWJyBd0dHRQUVER7DJEgiIrKwtjTLDLkB4KWICx1i43xuR/4eFLgDM7338OWAbc\n3fn4P6y1LcB2Y8wWYCKwKlD1iciXXXvttcyZMyfYZYgExT333MOMGTOCXYb00ECfRj3MWlvc+X4J\nMKzz/VxgdZfnFXU+9iXGmNuA2wBGjRoVoDJFwo/H4+GNN95g8uTJXH755cEuR2RA3XvvvRQVFQW7\nDOmFgQ4wPtZaa4yxffi6J4EnASZMmNDrrxeRffvss8+orq7muuuu4+abbw52OSID6uGHH8Za/Uhx\nk4FehbTHGJMD0Pm2tPPxXcDILs8b0fmYiAyQ1audQdBTTjklyJWIDDxjDB6PJ9hlSC8MdICZD9zY\n+f6NwL+7PH6VMSbWGDMa+AqwdoBrEwlrq1atIi0tjXHjxgW7FJEBZ4zRCIzLBKyFZIx5EWfC7hBj\nTBHwv8BvgbnGmG8BO4GpANbaTcaYucDHQDvwPa1AEhlYq1ev5qSTTiIiQttDSfiJiIhQgHGZQK5C\nuno/nzp7P8//NfDrQNUjIvtXV1fHxo0bueyyy4JdikhQqIXkPvpVS0R499138Xg8mv8iYUstJPdR\ngBER3wTeiRMnBrkSkeBQC8l9FGBEhNWrV3PYYYeRnp5+8CeLhCCNwLiPAoxImLPWsnr1ak4++eRg\nlyISNJoD4z4KMCJhbtu2bZSVlSnASFjTCIz7KMCIhDltYCeiOTBupAAjEuZWr15NYmIiRxxxRLBL\nEQkatZDcRwFGJMytWrWKiRMnEhkZGexSRIJGLST3UYARCWONjY188MEHmv8iYU8tJPdRgBEJY++9\n9x7t7e2a/yJhTy0k91GAEQlj3gm8J510UpArEQkutZDcRwFGJIytWrWKMWPGMHTo0GCXIhJUaiG5\njwKMSJiy1rJq1SrNfxFBLSQ3UoARCVNFRUUUFxcrwIigFpIbKcCIhKlVq1YBKMCIoBaSGynAiISp\nFStWkJCQwLHHHhvsUkSCTi0k91GAEQlTK1euZOLEiURHRwe7FJGgUwvJfRRgRMJQfX09H3zwAaee\nemqwSxEZFNRCch8FGJEwtGbNGjo6OhRgRDqpheQ+CjAiYWjlypUYY7QDr0gntZDcRwFGJAytXLmS\nI444grS0tGCXIjIoqIXkPgowImGmo6ODVatW8dWvfjXYpYgMGmohuY8CjEiY2bhxI3V1dZr/ItKF\nWkjuowAjEmZWrlwJoAAj0oVaSO6jACMSZlauXElOTg75+fnBLkVk0FALyX0UYETCzMqVKzn11FMx\nxgS7FJFBQy0k91GAEQkjRUVF7Ny5U+0jkS9QgHEfBRiRMKL5LyL7pjkw7qMAIxJGVq5cqQMcRfZB\nc2DcRwFGJIzoAEeRfVMLyX0UYETChA5wFNk/tZDcRwFGJEx4D3DUDrwiX6YWkvsowIiECR3gKLJ/\naiG5jwKMSJhYuXIlRx55JKmpqcEuRWTQUQvJfRRgRMJAe3s7q1at0vwXkf1QC8l9FGBEwsAHH3xA\nXV0dp59+erBLERmU1EJyHwUYkTCwfPlyAAUYkf1QC8l9FGBEwsDy5csZO3Ysubm5wS5FZFBSC8l9\nFGBEQpzH42H58uUafRE5ALWQ3EcBRiTEffzxx1RWVirAiByAWkjuowAjEuI0/0Xk4NRCch8FGJEQ\nt3z5ckaMGMHo0aODXYrIoKUWkvsowIiEMGutb/6LMSbY5YgMWmohuY8CjEgI27p1K8XFxWofiRyE\nWkjuowAjEsI0/0WkZ9RCch8FGJEQ9vbbb5OVlcVhhx0W7FJEBjW1kNxHAUYkhGn+i0jPqIXkPgow\nIiGqoKCAHTt2qH0k0gNqIbmPAoxIiHrnnXcAzX8R6QkFGPdRgBEJUcuXLyc1NZWjjjoq2KWIDHqa\nA+M+QQkwxpgfG2M2GWM2GmNeNMbEGWMyjDGLjDGfd75ND0ZtIqFi+fLlnHbaaURGRga7FJFBT3Ng\n3GfAA4wxJhf4ITDBWnskEAlcBUwHFltrvwIs7vxYRPqgtLSUTz/9VO0jkR5SC8l9gtVCigLijTFR\nQAKwG7gEeK7z888BlwapNhHX0/4vIr2jFpL7DHiAsdbuAh4CCoBioMZa+yYwzFpb3Pm0EmDYvr7e\nGHObMWadMWZdWVnZgNQs4jZvv/02iYmJHH/88cEuRcQV1EJyn2C0kNJxRltGA8OBRGPMdV2fY50Y\nvM8obK190lo7wVo7ISsrK+D1irjRkiVLOO2004iOjg52KSKuoBaS+wSjhXQOsN1aW2atbQPmAZOA\nPcaYHIDOt6VBqE3E9fbs2cPHH3/MWWedFexSRFxDLST3CUaAKQBONsYkGGd70LOBT4D5wI2dz7kR\n+HcQahNxvaVLlwIwefLkIFci4h5qIblP1EDf0Fq7xhjzMvAe0A5sAJ4EkoC5xphvATuBqQNdm0go\nWLp0KampqRx33HHBLkXENdRCcp8BDzAA1tr/Bf73Cw+34IzGiEg/LFmyhDPOOEP7v4j0glpI7qOd\neEVCSGFhIVu2bNH8F5FeUgvJfRRgREKId/6LAoxI76iF5D4KMCIhZOnSpWRmZur8I5FeUgvJfRRg\nREKEtZYlS5Zw5plnEhGhf9oivaEWkvvou5xIiNi+fTsFBQVaPi3SB2ohuY8CjEiIWLJkCaD5LyJ9\noRaS+yjAiISIpUuXkp2dzWGHHRbsUkRcRyMw7qMAIxICvPNfzjrrLJwNrkWkN7z/bhRi3EMBRiQE\nbN68mZKSErWPRPrIO/FdAcY9FGBEQoDOPxLpH+8IjFYiuYcCjEgIWLJkCSNHjmTMmDHBLkXEldRC\nch8FGBGX83g8LFu2jMmTJ2v+i0gfKcC4jwKMiMt9+OGHlJeXq30k0g+aA+M+CjAiLrdo0SIAzjnn\nnCBXIuJemgPjPgowIi63aNEijjjiCIYPHx7sUkRcSy0k91GAEXGx5uZm3nnnHb72ta8FuxQRV1ML\nyX0UYERcbMWKFTQ3NyvAiPSTWkjuowAj4mKLFi0iOjqaM844I9iliLiaWkjuowAj4mKLFi1i0qRJ\nJCYmBrsUEVdTC8l9FGBEXKqsrIwNGzaofSTiB2ohuY8CjIhLLV68GEABRsQP1EJyHwUYEZdatGgR\n6enpnHDCCcEuRcT11EJyHwUYERey1rJo0SImT55MZGRksMsRcT21kNxHAUbEhT777DMKCwvVPhLx\nE7WQ3EcBRsSF3nzzTUDzX0T8RS0k91GAEXGhRYsWMWbMGMaMGRPsUkRCglpI7qMAI+IybW1tLFu2\nTKMvIn6kFpL7KMCIuMyaNWuoq6tTgBHxI7WQ3EcBRsRlFi1aREREBJMnTw52KSIhQy0k91GAEXGZ\nN954gxNPPJH09PRglyISMtRCch8FGBEXKS8vZ+3atVxwwQXBLkUkpKiF5D4KMCIu8uabb2KtVYAR\n8TO1kNxHAUbERRYuXEhmZqaODxDxM7WQ3EcBRsQlPB4Pb7zxBuedd56ODxDxMwUY91GAEXGJ9957\nj7KyMrWPRAJAc2DcRwFGxCUWLlyIMYbzzjsv2KWIhBzNgXEfBRgRl1i4cCETJkwgKysr2KWIhBy1\nkNxHAUbEBSorK1mzZo3aRyIBohaS+yjAiLjAm2++icfjUYARCRC1kNxHAUbEBbzLp0888cRglyIS\nktRCch8FGJFBzuPx8Prrr3Puuedq+bRIgKiF5D4KMCKD3IYNGygtLVX7SCSA1EJyHwUYkUFu4cKF\nAFo+LRJAaiG5jwKMyCD3+uuvM2HCBIYOHRrsUkRCllpI7qMAIzKIVVVVsWrVKrWPRAJMLST3UYAR\nGcTeeOMNLZ8WGQBqIbmPAozIIDZ//nyysrKYOHFisEsRCWlqIblP1MGeYIwZAVwFnAYMB5qAjcB/\ngIXWWo23iQRAW1sbCxcu5NJLL9XyaZEAUwvJfQ4YYIwxzwK5wKvAA0ApEAeMA84Hfm6MmW6tXR7o\nQkXCzcqVK6muruYb3/hGsEsRCXlqIbnPwUZgHrbWbtzH4xuBecaYGGCU/8sSkQULFhATE8O5554b\n7FJEQp5aSO5zwDkw+wovxpixxpijOj/faq3d0tubGmPSjDEvG2M+NcZ8Yow5xRiTYYxZZIz5vPNt\nem9fVySULFiwgLPOOoukpKRglyIS8tRCcp9eTeI1xvwM+DlwuzHm7/2472PA69baw4BjgE+A6cBi\na+1XgMWdH4uEpc2bN/P555+rfSQyQNRCcp8DBhhjzA+NMV1nDx5jrb3FWnsrTvDoNWNMKnA68DT4\nRnGqgUuA5zqf9hxwaV9eXyQULFiwAICLLrooyJWIhAe1kNznYCMwFcDrxpiLOz9+0xjzujHmTeCN\nPt5zNFAGPGuM2WCMecoYkwgMs9YWdz6nBBi2ry82xtxmjFlnjFlXVlbWxxJEBrcFCxZwzDHHkJeX\nF+xSRMKCWkjuc7A5MM8D3wCONsbMB9YDU4ArrLV39fGeUcDxwBPW2uOABr7QLrJOBN5nDLbWPmmt\nnWCtnZCVldXHEkQGr4qKClasWKH2kcgAUgvJfXoyB2YsMBe4DfgezvyV+H7cswgostau6fz4ZZxA\ns8cYkwPQ+ba0H/cQca2FCxfi8XgUYEQGkFpI7nOwfWBmAW1AArDLWjvNGHMc8FdjzLvW2l/19obW\n2hJjTKEx5lBr7WbgbODjzutG4Ledb//d29cWCQULFiwgOzubCRMmBLsUkbChFpL7HGwfmOOstccA\nGGM2AFhrNwDfMMZc0o/7/gB4vnMfmW3AzTijQXONMd8CdgJT+/H6Iq7U2trK66+/zhVXXOH7jVBE\nAk8tJPc5WIB53RjzBhANvND1E9baPo+QWGvfB/b16+XZfX1NkVDwzjvvUFtbq/aRyABTgHGfAwYY\na+3dxpgUwGOtrR+gmkTC1oIFC4iNjeWcc84JdikiYUVzYNznYPvAXAfU7y+8dO7K+9WAVCYSZqy1\nzJ8/n7PPPpvExMRglyMSVjQHxn0O1kLKBDYYY9bjLKEuwznM8RDgDKAc7Zgr4hcffPAB27dv5+c/\n/3mwSxEJO2ohuc/BWkiPGWP+CEwGTgWOBppwtv6/3lpbEPgSRcLDvHnziIiI4OKLLz74k0XEr9RC\ncp+DjcBgre0AFnVeIhIg8+bN47TTTkMbNIoMPLWQ3EfrNEUGgc8++4xNmzYxZcqUYJciEpbUQnIf\nBRiRQeBf//oXAJdeqjNMRYJBLST3OWiAMcZEGGO0qZxIAM2bN48JEyYwatSoYJciEpbUQnKfgwYY\na60H+OkA1CISloqKili7dq3aRyJBpBaS+/S0hfSWMeZOY8xIY0yG9wpoZSJh4pVXXgFQgBEJIrWQ\n3Oegq5A6Xdn59ntdHrPAGP+WIxJ+5s2bx/jx4zn00EODXYpI2FILyX16FGCstaMDXYhIOCovL+ft\nt9/mnnvuCXYpImFNLST36VGAMcZEA98BTu98aBnwF2ttW4DqEgkLCxYswOPxqH0kEmRqIblPT1tI\nT+CcSP1458fXdz52ayCKEgkX8+bNIy8vj+OOOy7YpYiENbWQ3KenAeZEa+0xXT5eYoz5IBAFiYSL\nuro63nzzTb773e/6vnmKSHCoheQ+PV2F1GGMGev9wBgzBugITEki4eG1116jtbWVyy67LNiliIQ9\ntZDcp6cjMHcBS40x2wAD5AG3BKwqkTAwd+5csrOzOfXUU4NdikjYUwvJfXoaYFYAXwG86zw3B6Yc\nkfBQV1fHa6+9xrRp04iMjAx2OSJhTy0k9+lpC2mVtbbFWvth59UCrApkYSKhbMGCBTQ3NzN1qk7p\nEBkM1EJynwOOwBhjsoFcIN4YcxxO+wggBUgIcG0iIWvOnDnk5uYyadKkYJciIqiF5EYHayGdB9wE\njAAe6fJ4LfCzANUkEtKqq6t5/fXX+d73vuf7rU9EgkstJPc5YICx1j4HPGeMudxa+88BqkkkpM2f\nP5/W1la1j0QGEQUY9+npr38nGGPSvB8YY9KNMfcHqCaRkDZnzhzy8vI46aSTgl2KiHTSHBj36WmA\nucBaW+39wFpbBXw9MCWJhK7KykrefPNNpk6dqs3rRAYRzYFxn54GmEhjTKz3A2NMPBB7gOeLyD68\n8sortLe3c+WVVx78ySIyYNRCcp+e7gPzPLDYGPNs58c3A88FpiSR0DVnzhzGjBnD8ccfH+xSRKQL\ntZDcp0cBxlr7gDHmQ+Dszofus9a+EbiyREJPeXk5ixcv5qc//anaRyKDjFpI7tPTERistQuBhQGs\nRSSkzZs3j46ODrWPRAYhtZDcp0dzYIwxJxtj3jXG1BtjWo0xHcaY2kAXJxJK5syZw7hx4zj66KOD\nXYqIfIFaSO7T00m8fwSuBj4H4oFbgT8FqiiRUFNcXMyyZcu48sor1T4SGYTUQnKfHm8Daq3dAkRa\nazustc8C5weuLJHQ8uKLL+LxeLj22muDXYqI7INaSO7T0zkwjcaYGOB9Y8yDQDG9CD8i4W727Nmc\neOKJHHrooQd/sogMOLWQ3KenIeT6zud+H2gARgKXB6ookVCyadMmNmzYwPXXXx/sUkRkP9RCcp8D\nBhhjzOLOd79rrW221tZaa39prf1JZ0tJRA5i9uzZREZGavWRyCCmFpL7HKyFlGOMmQRcbIz5B9Bt\n9qG19r2AVSYSAjweD88//zznnXceQ4cODXY5IrIfaiG5z8ECzP8D7gVGAA/TPcBYYHKA6hIJCe+8\n8w6FhYU88MADwS5FRA5ALST3OWCAsda+DLxsjLnXWnvfANUkEjJmz55NUlISl1xySbBLEZEDUAvJ\nfXo0iVfhRaT3mpubeemll5gyZQoJCQnBLkdEDkAtJPfp8VECItI7r776KjU1NVx33XWBu4m10NEM\nHU1gDJhoiIgGEwURkYG7rx9Za30/NPb1wyMiIkKb/0nAqYXkPgowIgEye/ZscnJymDy5j1PFmsug\nZiPUfAyNRdC4C5p2Q9MuaKmAjkZob8SZjrYPJgpiMyAms/NtBsQPh8Q8SBjlvE3Kh/hcJ/z4ibWW\ntrY2WltbaWtro62tjY6ODtrb22lvb6ejo4OOjg48Hg8ej6dHv/EaY3xBJiIigoiICKKiooiMjOz2\n1nvFxMQQGRmp4CM9phaS+xwwwBhjMg70eWttpX/LEQkNFRUVvPbaa/zwhz8kMrIHIyHN5VC+EspW\nQOV7TnBpLt37+YhoJ3zED4fUIyF2CEQlQlQCRCZAZDxgwdMGtt1529EIrVVO2GmthIadUP5f5+Ou\nopIgZTykHg6p4yHtaMiYAHFZByzZ4/HQ0tJCS0sLzc3NNDc309LSQltb25ee6w0cUVFRREdHExcX\n5wsiXxxh6fq+d3Sma9jxeDy+ENTS0uILQ19kjCE6Otp3xcTEEBsbS0xMjC/giHipheQ+BxuBWY/z\n650BRgFVne+nAQXA6IBWJ+JSc+fOpa2tbf/tI08blP0Xdr/mXDUbnccjYiDtGBh+EaQdCWlHQcrh\nEJ8Nxk+bX7c3QENB57UNaj6B2k+gZBFsf27v8xLznSCTeSI2YwJtycfQ2BZJU1MTjY2NNDc3+77Z\nG2OIiYkhISHBFxZiYmJ84cH7wyFQvKGmvb3dN+rT9aqvr6e9vb3b10RFRREbG0tsbCxxcXHExcUR\nGxurYBOm1EJyn4OtQhoNYIz5K/Ava+1rnR9fAFwa+PJE3GnWrFkceeSRHHPMMXsf9HRA6TLY+Q8o\n/KczOhIRDVmnQf4M523mBIiMC2xxUYnOSEvq+C9/rrUaqt6HinfxVLyLLV9LZOHLGCCaSGLiD6c9\n6UQSMk4ladhpxCUP9YWAYLZrvCM50dHRxMfH7/M5Ho+H1tZW36iR9/3q6upuP7S8I0TeKz4+nujo\naLWjQpxaSO7T0zkwJ1trp3k/sNYu7DwTSUS+YOPGjaxdu5ZHHnnE+aZYvwO2PgXbnoGmYqdlM+JS\nGHkZZJ8D0SnBLhlwvnG3eOKoNYdTGzOC5vQLIB2iPTWk260kNX1IbM1q4iv+jil7Bj6LgPQTIHsy\n5JwPWac6gWyQioiI8IWSrrxzdrxtMG8rrK6uzvecyMhIEhISiI+P911RUZpCGEoUYNynp/8Cdxtj\nfgHM7vz4WmB3YEoScbdnn32W6Ogobv56Piy7EHYvdCbJ5nwdTrgRhl8IUfseJQiGpqYmampqqK2t\npbW1FYCEhASGDh1KUlIS8fHxGHPq3i9ob4Ty1c5o0p5l8Okj8PEDEJXsBLLhFzhXwoig/Hl6y9v+\niomJISVlb5j0eDw0NzfT1NTku7qGmujoaBISEkhISCAxMTHoo1DSP5oD4z49DTBXA/8L/Kvz4+Wd\nj4lIF60tzVR9+BQfP5pM2vopEDcMjrwXxn4LEkcFuzyf9vZ2qqurqa6uprm5GYDExEQyMzNJSUkh\nOvoAIylRCc6oS3bn6qq2OihZDMULnbBW1PltIu2ozpGmKc68Hpf9cI+IiPAFFK+Ojg6am5tpbGyk\nqamJhoYGampqvvT8xMRE4uPjAz73R/xHc2Dcp0cBpnO10e0BrkXEvayFXa/S9M4PeObmWhoihsOE\nB2D09YGf09JD1lrq6+uprKz0jSTEx8eTk5NDampq31si0ckw8lLnstZZ9r37Ndj9Kmz6NWy8z5kQ\nPOIyJ8wMOcU1e9R8UWRkJImJiSQmJgJ720+NjY00NDTQ2NhIaamzeswYQ1xcnO/5iYmJCjQuoBEY\n9+jRdyxjTBbwU+AIwPfd2Fqrs5BEytfAez+G8lU01CXws/lpPPavbRATG+zKAOc3yqqqKiorK2lp\naSEyMpLMzEzS09O/NB+k34yBtCOc6/C7nKXgRfOdUZnP/wSbH4W4oZB7CYz6JgybDBHunUvStf2U\nlpYGOKNbjY2NvquiooLy8nKMMcTHx5OUlKQRmkEqIiJCAcZFevqd43lgDnAR8G3gRqAsUEWJuEJT\nCXxwD2ybBXHZVI97gNGTpnPHnbcTNQjCS1tbGxUVFVRWVuLxeIiLiyM3N5fU1NSB+8EZNxQOudW5\n2mph12tOmNn5Amz9q/P5kVdA/jXOyIzL2kz7EhUVRUpKim8+jcfjoaGhgYaGBurr67uN0HhHZpKS\nkoiLi9McmiAzxqiF5CI9DTCZ1tqnjTG3W2vfBt42xrwbyMJEBi3rgS1Pwvt3O1v4j/8pHPkL/vzw\nn2hts9xyyy1BLa+trY3y8nIqKyux1pKSksKQIUM6J+MG8QdkdArkX+Vc7U3OnJkdL8K2p53RmcQ8\nyLsK8q52NtMLkR/mERERJCcnk5ycDDjzaLxhpqGhgT179rBnzx5feyopKYnk5OQDz0OSgDDGaATG\nRXoaYLxbaxYbYy7EWYF0wF16D8YYEwmsA3ZZay/q3PV3DpAP7ACmWmur+nMPEb+r/RzW3gqly2HY\n2XDi45AyDmstzzzzDKeffjqHHHJIUEpra2ujrKyMqqoqrLWkpaWRlZVFbGzwR4O+JCremQ8zcooz\nMlP4Cux8ET55yFnRlHq4E2TyrobkscGu1q8iIyO7jdC0tbX5Ak19fT21tbUAxMbG+sJMQkKC2k0D\nQC0kd+lpgLnfGJMK3AH8AUgBftzPe98OfNL5WgDTgcXW2t8aY6Z3fnx3P+8h4h/Wwmd/gvfvgohY\nOOlpGHOzb5Rg5cqVfP755/z85z8f8NI6OjooKyujosI5IsAbXGJiYga8lj6JToExNzhXcxkUvuyM\nzHx4r3NlnQZjboRRVwyaPXP8KTo6mrS0NNLS0py9eFpaqK+vp66ujsrKSioqKjDGkJSU5LsGZSgN\nAWohuUtPVyG92vluDXBWf29qjBkBXAj8GvhJ58OXAGd2vv8csAwFGBkMmkth9S2w+z+QcwGc9BQk\nDO/2lGeeeYakpCS++c1vDlhZ3sm5paWldHR0kJqayrBhw9wTXPYlLgu+8h3naiiAHS/A9lmw5lZY\n9wNnxGbMTTD0LNeuZDoQ78qluLg4hgwZ4ps/U1dX5ws14ISe5ORkX6DR6Ix/qIXkLj1dhTQCZ+Tl\nqzhnI70D3G6tLerjfWfirGpK7vLYMGttcef7JcCwPr62iP/sWQorr3a22D/h9zDu+1+am1FXV8fc\nuXO5+uqrfctrA62uro7i4mJaW1tJTEwkOzt7v1vou1biKDhiOhx+N1SsdSZL7/wH7HgeEkY6S9RH\n3wgp44JdacB8cf6Md3Smvr6e6upqKisrfaMz3kDj6gAbZGohuUtPW0jPAi8AV3R+fF3nY1/r7Q2N\nMRcBpdba9caYM/f1HGutNcbs82+RMeY24DaAUaMGz8ZgEmKshU9+56wySh4Hkxc5G7Ptw/PPP09D\nQwO33nprwMtqbW2lpKSE2tpaYmJiyMvLIykpKbRXrxgDQ05yrhMedZZlb38OPv4tbJrhrF4acxOM\nuhJiUoNdbUB5z53KzMzE4/HQ2NhIXV0dtbW1vtGZuLg4X+gJ+sRtl1ELyV1MT9KmMeZ9a+2xB3us\nRzc05jc0IzeQAAAgAElEQVTA9UA7zp4yKcA84ETgTGttsTEmB1hmrT30QK81YcIEu27dut6WIHJg\nbXWw+iYonOfMuzjpaWeztn2w1nLssccSGRnJ+vXrA/bDwlpLRUUFpaWlWGvJyspiyJAh4d06aCqG\n7bOdMFOzCSLjYdRUGHurcy5TGP3g9s6dqauro66ujsbGRsCZMOwNM0lJSTpp+yCSk5OZNm0ajzzy\nSLBLCWvGmPXW2gkHe15PR2AqjDHXAS92fnw1UNGXwqy19wD3dBZ5JnCntfY6Y8zvcPaX+W3n23/3\n5fVF+qWhAN7+hvMD8biH4LCfHPAH4apVq/jwww/5y1/+ErDw0tTUxK5du2hubiY5OZmcnBy1CQDi\nc5zN8sbfCZXrYOszTntp+3OQcqgTZEbf4Ow1E+K6zp3Jysqio6PDF2bq6uqorq4GnOMivIFGE4G/\nTC0kd+lpgLkFZw7MozhzYP4L3OTnWn4LzDXGfAvYCUz18+uLHFjFu0546WiCM1+DnHMP+iVPPPEE\nycnJXHPNNX4vx1pLWVkZpaWlREVFMWrUqG6HDUonYyDzROc6/iEoeMk5/XvDXfD+PTDiEifMZH8t\nJCf+7ktkZGS3lU3eVlNdXR0lJSWUlJT4Dq/0LtNWq0ktJLfp6SqkncDFXR8zxvwIZzJun1lrl+Gs\nNsJaWwGc3Z/XE+mzovmw8iqIy4azlzj7kBxEeXk5c+fOZdq0aSQlJfm1nK6jLqmpqeTk5PT9rKJw\nEpXozIcZcxPUfAJbn3ZGZAr/6Uz8HXMLjL3Z2TQvTHTd8Tc7O5vW1lZfmPEec+BtNaWkpJCYmBi2\nrSatQnKX/nxH/An9DDAig8K2v8GaWyDjBDhjQY9bDrNmzaK1tZVvf/vbfivFWkt5eTmlpaVERERo\n1KU/Usc7IzLHzIBd851RmY2/cq6cc51RmdyLITK82nExMTFkZmaSmZlJR0eHb/O82tpaqquru61q\nCrcdgdVCcpf+BBiNN4r7bf4DrP+hs6vu6a9AdM9GUjweD3/+85/56le/ypFHHumXUtra2ti1axf1\n9fWkpKQwfPhwjbr4Q2SMc3DkqG9Cw07Y+ixsewZWXAGxQ5yl2GO/5QSeMBMZGUlqaiqpqalYa317\nznRd1RQfH+8bnYmNjQ3pVpNaSO7Sn++OiqniXtbCxvvgo/+FEZfCqS9CZM9PZn7rrbfYunUrv/rV\nr/xSTl1dHUVFRXg8HoYPH056enpI/6AImsQ8OPr/4Mh7oeRNZ1Rm82Pw6cPOyqWx05yVZ1EJwa50\nwHXd7Tc7O5uWlhZfkCktLaW0tJTo6GjfvJnExMSQ+zuqFpK7HDDAGGPq2HdQMUCI7ZolYcNaeO8n\nsHmmM1di4l8hondZ/oknnmDIkCFcfvnl/SzFsmfPHsrLy4mNjWXkyJHExfU8SEkfRUTC8Aucq2kP\nbP+bE2ZW3+SMyOVf64SZjOOCXWlQdF3VNHToUNra2nwjM97jDSIjI0lKSiIlJSVklmirheQuB/yu\nba3d9+YXIm5lLbx3hxNeDr0djn8ETO/2UikqKmLBggXccccd/VqK2t7eTmFhIQ0NDaSnp5OTkxPe\n+7oES/ywvcuxy96BLX+Fbc/C509A+vFwyDTnUMkQ3yTvQKKjo8nIyCAjI8M3b8Y7EbimpsY3Udg7\nb8aty/zVQnIXNdglfFgLH/wMNj8K434Ixz/ap83OnnrqKTweD//zP//T51IaGxspKCigo6OD3Nxc\n0tPT+/xa4ifGwNDTnav197D9edj6V3j3O07ozevcJG/IpLDaJO+LvjhvprGx0ddqKi4upri4mLi4\nOF+rKS4uzjWtJrWQ3EUBRsLHxl85288f8j9wwsw+/RBqaWnhL3/5C+effz5jxozpUxmVlZUUFxcT\nFRXFmDFjQu8Mo1AQkw6Hfh/Gfc/ZJG/LX2Hni855TCnju2ySNyTYlQbVF5dod90NuOu8Ge8k4ISE\nhEE9yqgA4y4KMBIeNv0WPvo/GHMznPh4n3+Dnjt3LiUlJfzoRz/q9ddaaykuLqayspKkpCRGjBih\nVUaDXbdN8h6Bgjmw5SnYcIdzTtaIy+CQW2HY5F63IkPNF3cDbm9v982bqaqqorKystvhlMnJyYNu\n3ozmwLiLvntK6Nv8B+eHTf61zoTdPv6gsdby6KOPMn78eL72td6dY9rR0UFhYSH19fUMGTKEYcOG\nuWZYXTpFJznLrcd+C6o/6twk729OqEkc7Tw+5mZIGB7sSgeFqKgo0tPTSU9Px+Px+ObN1NbWUlNT\nA0BCQoKv1TQYjjbQHBh3UYCR0FbwMqy/3VkqffKsfm0lv2LFCjZs2NDrc49aW1vZuXMnLS0tDB8+\nnIyMjD7XIINE2lFOG/LY30Lhv5y5Mh/+Aj76fzD8QmcF0/ALer26LVRFRESQkpLi29+oqanJN2/G\ne7RBbGysb2QmWEcbqIXkLvrXJaGrdDn89zrImgSTXuj3D5OZM2eSkZHBdddd1+OvaWhooKCgAID8\n/Hy/HzkgQRYZB/lXO1fdFmdUZtuzsGsBxA93RmTG3gJJfZsvFYqMMSQkJJCQkNDtaIPa2tovHW0w\n0Kdoq4XkLgowEpqqN8Hbl0DSaDh9PkT1b6Ls9u3beeWVV7j77rtJSOjZJmc1NTUUFRURHR1NXl7e\noBgilwBKPgSO/Q0c/SvY9R9nX5mPfwObfg3Z5zgTf0dcCpH6e9DV/o428J6i7Q083lZTIJdoq4Xk\nLgowEnoad8Gy853fjs9cCLH9b9n88Y9/JCIigu9973s9en5FRQXFxcXEx8eTl5enybrhJCIaRl7q\nXA2FzsqlbU87h4XGDnFWL429NSyPLjiYfS3R9o7OeJdox8bG+sJMfHy8X1tNaiG5i76rSmhpq4Vl\nF0BrDXxtOSTl9/sl6+rqeOqpp5g6dSq5ubkHfK61ltLSUsrKykhOTmbkyJGDetmoBFjiSDjqXjji\nZ7BnsbMc+7M/wKePdB5dcGvn0QWJwa500NnfEu3a2lrKysooKysjKiqqW6upv//W1EJyFwUYCR2e\nDlh5LdR87Iy8pB/rl5edNWsWtbW13H777Qd8nrWW3bt3U1VVRVpaGrm5uVppJI6ISOcE7Jxzobm0\ny9EFNzuTzPOucXb8zTg+2JUOWrGxscTGxjJkyBDa29t9q5pqamqoqqryyynaaiG5iwKMhI4PfwG7\nX4UJf4Sc3i1z3h+Px8Njjz3GpEmTmDhx4gGfV1hYSF1dHVlZWQwdOlThRfYtbqhzbMFhd0DZCmdU\nZvss2PJnSDvGOZ8r/xrnebJPUVFRpKWlkZaWhsfj6dZq8p6iHRcX5wszPW01qYXkLgowEhq2P793\nl92vfNdvLzt//ny2bt3KjBkz9vscj8dDQUEB9fX15OTkkJmZ6bf7SwgzBoae5lytv4cdLzjzZd77\nMWy4C3IvhNE3OsuyI915ttBAiIiI+NIp2t7dgL2tJu/Bk95W0/7mpKmF5C7Gzf+zJkyYYNetWxfs\nMiTYKt6FRafBkJPgrEV++2ZvreWUU06hrKyMzZs37/ObXkdHBzt37qSxsVF7vIh/VG+C7c/B9r9D\ncwnEZkLetc7ITPqxYX0OU2998eDJjo4OAOLj432jM13Paho/fjxHHXUUc+fODWbZYc8Ys95aO+Fg\nz9MIjLhb425YfgnE58BX/+nX31Tffvtt1qxZwxNPPLHf8LJjxw6ampoYMWIEaWlpfru3hLG0I+C4\nB+GYGVCyyBmV2fJn+Oz3kHa0MyqTf61zirYc0BdXNTU1NfkCjfespqioKN/oTEJCgkZgXEQjMOJe\nHa3w1hlQsxHO/a+zO6ofnX/++bz//vvs2LGDuLi4bp9rb29nx44dtLS0MHLkSFJSUvx6b5FuWiqd\nIwu2zYKKtWAiYfjXnVGZ4RepxdQH3rOa6urqqK+vx+Px0N7eTmFhIaeccgpJSUnExsZqLlsQaARG\nQt+GO6FiNXz1Jb+Hlw0bNvDGG2/wm9/85oDhZdSoUSQnJ/v13iJfEpsBX/mOc9V8DNuegx1/d3b8\njc2EvKudkZmME9Ri6qGuZzV595yZOXMmJ510EiUlJQC+k7STkpJITEwcdIdPhjuNwIg77fgH/Pdq\nOPTHcMIjfn/5q666ioULF1JQUEBqaqrvcYUXGTQ87VDyljMqU/QKeFogeRzkXwejr9XxBX1w9NFH\nM3bsWObMmeNrNdXX12Ot9e0I7J0w3HXujPiXRmAkdNV8DGtvdTYCO+4Bv7/8li1beOmll/jpT3/a\nLbx457wovMigEBEFw893rtYqKPgn7JjtHCj50f+DIac4c2VGTYW4rGBX6wreVUgxMTFkZGSQkZHh\nW6btDTR79uxhz549vrkz3ku7bQ88/RcXd2mrh3e+6exceuocZ9t2P/vd735HdHR0t43rFF5kUItJ\nh0Nuda6GAtj5Iux4HtZ9H9b/CHLOc0ZmRlwMUT07yysc7Wsjuy8u025ra6O+vt4XaKqrqwFn3xnv\nZOD4+HjtwD0AFGDEPayFtdOgbrOzXDrhwNv690VxcTGzZs3illtuITs7G9gbXpqbmxk5cqTCiwxu\niaPg8Ludq+pDJ8jsfAH++x+ISoKRU5yRmWGT+31Ce6jpyUZ20dHR3ebOeFc21dfXU15eTnl5ORER\nESQmJvqCjw5yDQz97RX32PIX2PkPOObXkD05ILeYOXMm7e3t3HnnnUD3pdKjRo3SaiNxl/SjnevY\n30DpcifMFLzkHGUQN8yZ/Jt3NWSeqMm/9H4jO++8mISEBIYOHUpHRwcNDQ2+uTPeXYFjYmJ8E4ET\nExPVbvIT/VcUd6je6OxQmnMeHD49ILcoLy/n8ccfZ+rUqYwdO9a3w25TU5OWSou7mQgYdqZzTfgD\n7H7NCTOfPw6bZ0JivjNXJu9KSD8ubMNMf89CioyMJCUlhZSUFKy1tLa2+kZnqqurqaysBJyN9Lwj\nNAkJCWo39ZECjAx+7Y2w8iqIToWTn3O+GQfAww8/TENDA/feey/WWgoLC2loaGDEiBHdJvOKuFpk\nnNNGGjkFWquh6N+wc45zQvYnD0LSIU6QybsSUo8MqzDjz7OQjDG+AygzMzN9S7UbGhq6tZu6rm5K\nTEzs8blNogAjbvDeHVCzCc56I2C7j5aVlfGHP/yBq666ivHjx7Nr1y7q6urIycnRDrsSumLSYMyN\nztVSAYX/goK5zrlim34NKeP3jsykjg92tQEXyMMcjTG+FpK33eRd3dTQ0MCePXuAvZOGvSM0MTEx\nCjT7oQAjg1vBP51t1MffBTnnBuw2v/vd72hqauLee++luLiY6upqhg4dqoMZJXzEZu5dydRcBoX/\ndEZmNv4KNv7S2Sxy1JVOoEn5SrCrDYiBPMwxMjLSdx4TOHtMecNMfX09tbW1AN2WaycmJhId7f+V\nl26lACODV0MBrLkVMk6Eo+8P2G1KS0v505/+xNVXX01GRgZlZWVkZmaSlaW9MyRMxWXBV77tXE3F\nUPCyc5TBh79wrvTjYOTlThsqhEZm+jsHpj+ioqJIS0sjLS3NN3/GG2a6LteOiYnxjeSEe6BRgJHB\nydMB/70WbAec+mJAz3p58MEHaW5u5uc//zllZWWkpaWRnZ2tYVsRcA5KPfQHztVY5Kxi2jl3b5hJ\nOWzvnJr04109ZyaQLaTe1uGdP5ORkYG1lubmZhoaGmhoaKCmpoaqqiogvAONAowMTp8+BGUr4JS/\nQfLYgN2mpKSExx9/nBkzZtDR0UFKSgq5ubkKLyL7kjACDvuxczXuciYAF86Djx+ATTMgYdTeMDNk\nEkS46+yggWwh9YYxhvj4eOLj4xkyZIgCTScFGBl8qj6AD++Fkd90dg8NoAcffJATTjiBCy+8kMTE\nREaMGKHwItITCbkw7rvO1VLhHCxZOA8+f8JZmh03FEZcCiOmwLCzXHFidjBbSL2hQONQgJHBpaMF\nVl0PMZlw4hMBHY4uLi5myZIlPPPMM8TFxTFq1CjtxyDSF7GZMOYm52qrg90LnTCz4wXY8qSzBULu\nN2DkZZB9LkQnBbvifRosLaTe6m2g8W6+l5iY6OpVTgowMrh8eC9UfwRnvApxQwJ6q8cee4zHHnuM\n2NhY8vLyiIx013C3yKAUnQx5U52ro9k5MbtwntNu2jEbImKcEZncbzhX4qhgV+wTERHhihGYgzlY\noOk6KTgyMpLExERfqHHTPjQKMDJ4lC6HTx6CQ26D3AsDeqvNmzdz2mmnkZiYyCGHHBJSw6oig0Zk\nHORe5Fyedmde264FzrXu+86VdjTkXuyEmcwJAduosifc0kLqrX0FGu8qJ+/met5l2xEREb6dgr2h\nZrCOTCvAyODQVgurboSkMXDcwwG9lcfjYfv27QwbNozs7GwdtCYyECKi9h5ncPzDULsZiuY7Yebj\nGbDpfojLdn55yf0GZJ/jnDo/gNzaQuqtL65yAmhra/OFmcbGRkpLS33Pj4+P97WcEhISBs1ZToOj\nCpH1P4bGAjjnnYD2x621fPjhh+Tk5LB27VqmTZsWsHuJyAGkHAqH3+VcLRXOvJldC5xl2lufdkZv\nhk3ubDVd5KyACrDBugppIERHR5Oamuo7NsW7U7A31FRWVlJRUQE482iGDh0a9F3KFWAk+Ir+Ddue\ngcPvgaxJAbuNtZbdu3cTFRXFzJkzmTlzZsDuJSK9EJsJo69zro5WKHtn7+jM7tfg3e845zINPx9y\nzoesr0Kk/0dOQ7WF1Bdf3CnY4/HQ1NTkCzWDoa2kACPB1VwOa2+D9GPhqP8L6K0qKiqoqqri6aef\n5oQTTiApaXCuhBAJa5ExkH22c50wE2o+huKFsPt12PyYM08uKtEZnck5H4ZfAEmj/XLrcGkh9UVE\nRIRvKfZgoQAjwbX+dmitgrMWBXSfiNraWkpKSlizZg0LFy7kwQcfDNi9RMRPjIG0I5xr/J3QVg97\nlkLx63tbTgDJ45wgk3M+DD0DouL7dLtwbiG5kQKMBE/RfNj5gjPykn50wG7T1NREUVER9fX1fP/7\n32f27NladSTiRtFJMOIbzmUt1H2+N8xs+YszQhMZB0PPdMJMzrnOUQc9XBasFpK7KMBIcLRWw7vf\ndk64PfyegN2mra2NnTt3Yoxh2rRpHHvssUyZMiVg9xORAWIMpIxzrkN/CO1NzlYMxQudUPPej5zn\nxefAsHOcVU3ZZzs7CO/3JdVCchMFGAmO9+6A5lI4Y0HAWkcej4eCggI8Hg+vvvoqGzduZMWKFa7Z\npElEeiEqHoaf51wA9TucTfRK3nICzY6/O4+nHAbDznYCzbAzIWbvShq1kNxFAUYGXvGbnauOpkPG\nCQG5hbWWoqIimpqaiI6O5t577+W6667j1FNPDcj9RGSQScqHQ251LutxdvgueQtKFsO2Z+HzPzmb\n5mVM6BydOYeYSI9aSC5i3Jw2J0yYYNetWxfsMqQ32urgP0c6vy1d8L7Trw6AkpISysvLyc7O5uab\nb2bZsmV89tln5OTkBOR+IuIiHa1QsdoJMyVvQcUasB20tEewoSiBky+5x5kMnHGiKw6hDDXGmPXW\n2gkHe55GYGRgvT8dGgvhaysCFl6qqqooLy8nPT2dtWvX8uqrr/Lggw8qvIiIIzIGhp7uXEf/0tkJ\nvHQ5b/71dg5L3w0f/LzzefEw5BRnUvCwMyDzpIDsPyN9oxEYGTh73obFZ8KhP4ITHg3ILRoaGtix\nYwcJCQnk5ORw9NFHY4zho48+IiZGv0mJyP5NmTKFzz//nI/eXepsprdnGZS+DdUfAhYiYjsDzRmd\ngebkPi/Zlv3TCIwMLu2NsOZbzllHx9wfkFu0tLRQUFBAdHQ0o0aN4uGHH+bzzz9n4cKFCi8iclC+\nVUhxQ2DkZc4Fzl5Vpe84YWbPMth0H2z8pXOyduZJTqAZejoMOdk5jVsGhAKMDIwP74X6rXD2koAc\n0Nbe3s7OnTsByMvLY8+ePdx3331cfPHFnH/++X6/n4iEnv2uQopJhxEXOxdAa41zsrY30Hz8G+cw\nShPhnK495FTI6rwSRw3onyGcDHiAMcaMBP4GDAMs8KS19jFjTAYwB8gHdgBTrbVVA12fBED5avj0\nUTjk2zDsLL+/vLWWwsJC2trayM/PJzY2ljvvvJO2tjYefTQwrSoRCT093sguJrXz1OwLnY/bap3v\nc2UrnWv7LGeVEziHUHYNNGlHOydzS78F479iO3CHtfY9Y0wysN4Yswi4CVhsrf2tMWY6MB24Owj1\niT91NMPqW5x/xMc94PeX9x7Q2NDQQG5uLomJibz22mu8+OKL/PKXv2TMmDF+v6eIhKY+b2QXneLs\n+ptzrvOxp92ZN+MNNOUroWCO87moJKft5A00Q052vl56bcADjLW2GCjufL/OGPMJkAtcApzZ+bTn\ngGUowLjfxvug9hM4c2FA/pF6D2jMysoiPT2d+vp6vvOd73D44Yczffp0v99PREKX33bijYiCjOOd\n69AfOI81FHQPNJvud/anwUDqeMic6ASbzJOcHco1SnNQQf0vZIzJB44D1gDDOsMNQAlOi2lfX3Mb\ncBvAqFHqLQ5qle/Bxw/AmJtguP/noXgPaExJSWHo0KEA/OIXv6CwsJAVK1Zo4q6I9EpAd+JNHOVc\n+Vc7H7fVQvkap/VUsQZ2vQrbZjmfi4x3NvnMPAmGdIaahJE9PtMpXAQtwBhjkoB/Aj+y1tZ23d7d\nWmuNMfv8W2StfRJ4Epxl1ANRq/RBRyusvhlis+D4R/z+8k1NTRQWFhIfH8+IESMwxrBmzRp+//vf\n893vfpdJkyb5/Z4iEtoG9DDH6BTI+ZpzgXM4ZcMOJ9RUdF6f/RE+fdj5fFz23jCTOdEJOF2OQQhH\nQQkwxphonPDyvLV2XufDe4wxOdbaYmNMDlAajNrETz5+wOkBn/YvZwa/H3kPaIyKimLUqFFERETQ\n1tbGtGnTGD58ODNmzPDr/UQkPAT1MEdjIGm0c+Vf5TzW0ep8H61YszfYFP1779ckjXWCjO863u/f\nbwezYKxCMsDTwCfW2q6/ms8HbgR+2/n23/v4cnGD6o3OPgmjroSRl/r1pT0eDzt37sTj8TBmzBii\no6MB+N3vfsdHH33Ev//9b1JSNCFORHpv0B3mGBkDmROca9z3nMdaKqHiXahaD5XrnVBTMHfv1ySN\n+UKoOSFkQ00wRmBOBa4HPjLGvN/52M9wgstcY8y3gJ3A1CDUJv3laXdWHUWnwoQ/+PWlvculm5ub\nycvLIy7OOYrg008/5Ve/+hXf/OY3ufjii/16TxEJHwPaQuqr2Izup24DNJdD1XtOoKlc7wScgpf2\nfj5x9JdHamIzB752PwvGKqQVwP5mIp09kLVIAGyeCZXvwqQXIS7Lry+9Z88e6urqyMnJITnZ2e2y\nra2N6667jqSkJP7wB/8GJhEJL0FtIfVH3JDuy7gBWiqchRTeUFO5Hgpf3vv5+FxIPwbSjtn7Nvkr\nEBE58PX3kdZpif/UfubsuDviEsi70q8vXVlZSXl5ORkZGWRkZPgev++++1i/fj3//Oc/yc7O9us9\nRSS8DLoWUn/EZnafJAxO+6nqPah6H6o+gOoPoPhNsO3O5yPjIfVISD+6e7CJSQ3On+EgFGDEP6zH\nOesoIg4mPO7X5X719fXs3r2bpKQkcnJy8K5YW716NTNmzOCGG25gypQpfrufiIQnV7SQ+iM2A7LP\ncS6vjhZnr66qD/aGmqJXYOvTe5+TmNc90KQd5UwgDvJojQKM+Mdnjztng5z8LCQM99vLeg9ojI2N\nZeTIkb7w0tDQwPXXX09ubi6///3v/XY/EQlfrm0h9UdkLKQf61xe1kLT7r2Bxvt296udm+8BR98H\nR/4iODV3UoCR/qvfDh9Mh5zzYPSNfntZ7wGNxhjy8vKIjNyb9u+66y62bt3K0qVLSU0dnMObIuIu\nIdVC6g9jICHXuXK/vvfx9kao+RiqP3JWRgWZAoz0j7d1RARMfNJvrSOPx0NBQQFtbW2MHj262666\nCxcu5IknnuCOO+7gjDPO8Mv9RERCvoXUX1EJe5d1DwIKMNI/W56EPUth4l/8dmy894DGxsZGRowY\nQUJCgu9zJSUl3HzzzRx55JHcf//9frmfiAiEaQvJxRRgpO/qd8CGu5wJYWOn+e1ly8rKqK6uZujQ\noaSl7d0qu6Ojg2uuuYba2loWL17s2wdGRMQf1EJyFwUY6RtrYc2tzvsnPeW31lF1dTWlpaWkpaWR\nldV9H5n777+fpUuX8swzz3DEEUf45X4iIl5qIbmLAoz0zda/wp7FcOKfnSV2ftDQ0MCuXbtITExk\n+PDhdD3gc8mSJfzyl7/k+uuv56abbvLL/UREulILyV0igl2AuFDDTnjvDhg2GQ65zS8v6V0uHR0d\n7Tug0WvXrl1cc801jBs3jscff7xbsBER8Re1kNxFIzDSO9bCmmmAhZOe9kvryLtcGiA/P7/bcumW\nlhYuv/xy6uvrWbx4MUlJSf2+n4jIvqiF5C4KMNI7W5+GkkVw4uOQlN/vlzvQcmmAH/7wh6xZs4aX\nX35Z815EJKDUQnIXtZCk5xoK4L2fwLCz4JD/6ffLWWvZtWvXPpdLAzz11FM8+eSTTJ8+ncsvv7zf\n9xMRORC1kNxFAUZ6xnpg9c2Ap3PVUf//6uzZs4eamhqGDRv2pd10ly9fzne/+13OPfdc7fciIgNC\nLSR3UQtJembz72HPEpj4V0ga0++Xq6iooLy8nPT0dIYMGdLtc1u2bOGyyy5jzJgxzJkzp9ucGBGR\nQFELyV00AiMHV70J3p8Oud+Asd/q98vV1NRQXFxMcnLyl5ZLV1VVcdFFFwHw6quvdtvITkQkkBRg\n3EUjMHJgHa2w6jqITnFGX/q56qi+vp6ioiISEhK6nS4N0NrayhVXXMG2bdtYtGgRhxxySH+rFxHp\nMc2BcRcFGDmwj/4Pqt6H01+B+GH9eqmmpiYKCgqIiYn50l4vHo+Hm2++mcWLFzNr1iwd0igiA05z\nYLWiC+EAABR+SURBVNxFLSTZv9IV8MkDTttoxCX9eqnW1lZ27txJREQEeXl5REV1z8533XUXL7zw\nAjNmzODGG2/s171ERPpCLSR3UYCRfWutgVU3QEIeHP9ov16qvb2dHTt2YK0lPz//S3u9PPzwwzzy\nyCN8//vfZ/r06f26l4hIX6mF5C5qIcmXWQtrp0FjAZyzHKKT+/xSHR0d7Nixg7a2NvLz8790gvST\nTz7JnXfeyRVXXMHMmTN1TICIBI1aSO6iACNftuVJKHgJjv0tZE3q88t0dHSwc+dOWlpaGDVqFImJ\nid0+//e//51vf/vbfP3rX2f27NlaLi0iQaUWkruohSTdVX0I62+HnPNg/F19fhnvEQHeXXaTk7uP\n4rz00kvcdNNNnHXWWbz88stfaiuJiAw0tZDcRQFG9mqrh5VTITYDTvlbn3fb9YaXhoYGcnNzv7TL\n7ksvvcQ111zDpEmTmD9/PvHx8f6oXkSkX7wtbIUYd1CAkb3WfR9qP4NJz0Pc0D69hLWWoqIi6uvr\nGT58OOnp6d0+P3v2bK666ipOPvlk/vOf/3yprSQiEiwKMO6iACOObbNg+3Nw5L3OYY19YK2lsLCQ\n2tpasrOzycjI6Pb5p59+mhtuuIEzzzyT119/nZSUFD8ULiLiH969qRRg3EEBRqByPaz9thNcjry3\nTy/hbRt5w8sXzzd66KGHuPXWWznvvPN49dVXNfIiIoOOdwRGK5HcQQEm3DWXwfIpTsvo1DkQ0fuF\naR6Ph8LCQurq6sjJyekWXjweD3fccQd33XUXU6dO5ZVXXtGcFxEZlNRCchctow5nnnZYeSU074Fz\nV0JcVu9fonPkxTvnpWvbqKWlhVtuuYUXXniBH/zgB8ycObPb8QEiIoOJWkjuogATzt67A/YshZOf\ng4wTev3lHR0d3VYbdZ2wW1ZWxpQpU1ixYgUzZsxg+vTp2qRORAY1tZDcRQEmXH32J/js93Doj2DM\nDb3+8ra2Nnbu3ElzczMjRowgLS3N97lNmzZx0UUXUVJSwj/+8Q+uvPJKf1YuIhIQaiG5i8bzw9Gu\n12D9DyH3YjjuoV5/eUtLC9u2baO1tZW8vLxu4WXevHmccsopNDc38/bbbyu8iIhrqIXkLgow4abq\nA2feS9oxzn4vEb3bvr+xsZFt27bh8XgYPXq0b4fd1tZWfvzjH3P55Zczfvx41q5dy8SJEwPxJxAR\nCQi1kNxFLaRwUrcVlp4P0alwxgKITurVl9fW1lJYWEhUVBT5+fnExsYCUFhYyNSpU1m9ejU/+MEP\neOihh3Q0gIi4jlpI7qIAEy4ad8OSr4Ftg7MXQ0Juj7/UWktZWRmlpaXEx8eTl5dHVJTzV+e1117j\nhhtuoLW1lblz53LFFVcE6k8gIhJQaiG5i1pI4aClApZ+DVrK4MyFkHp4j7+0o6ODwsJCSktLSUtL\nY/To0URFRVFbW8ttt93GhRdeyPDhw1m3bp3Ci4i4mlpI7qIAE+paKmDJuU776IwFkHlij7+0tbWV\nbdu2+XbXzc3NJSIigrfeeoujjjqKp59+mrvvvpu1a9cybty4AP4hREQCTy0kd1ELKZQ1l/7/9u49\nOK7yvOP494ctx4rkKwYPAhvbjbFJDMFgaNxAakISIEkLadqAaQuZZMbJECAtTYrbdOKYSQfS3CCT\nACExA6khhBJuw6VACC024WIg4uILjm3A2AYtvlBJ1sWW9ukf58hZS7uyZCztHvH7zGh8zrtnz3ne\n1+/sPvu+Z/eF334sWaDxI3fCxHl9fmpjYyObN28GYMqUKdTW1tLY2MjChQu59tprOeqoo1i+fDlz\n584doODNzAaXE5hs8QjMUNX6BvxmHjStg3n3Qd2ZfXpaPp9ny5YtbNy4kaqqKqZNm0ZNTQ1Lly5l\nxowZXHfddVx66aXU19c7eTGzIcX3wGSLR2CGosaX4dEzoT2X3PMy8c/79LS2tjY2bdpEW1sbBx98\nMBMnTuT555/n4osv5vHHH+fEE0/k7rvv9tejzWxI8j0w2eIRmKEmtwwemgsdzfDR3/YpeYkIcrkc\n69evZ/fu3Rx55JHk83kuvPBC5syZw9q1a1myZAlPPvmkkxczG7I8hZQtHoEZSl65GZ76AtROhXn3\nQ+20fT6ltbWVLVu20NraypgxY6iqqmLRokVcc801dHZ2ctFFF7F48eK9fm3XzGwo8hRStjiBGQo6\nd8Hv/wnW/hgO/Qiccie8Z3zvT+nspKGhge3btzN8+HAOPfRQrrrqKq6++mra2to4//zz+eY3v8nU\nqVMHqRJmZuXlKaRscQKTdTtfh+V/A9uegpmXwnFXwkFVJQ+PCHbs2EEul6Ojo4Px48czceJEFi9e\nzBVXXMH8+fNZtGgRM2bMGMRKmJmVn6eQssUJTFZFwKtL4ZlLIDrh5Nth8md7OTxobm7mzTffpL29\nfc8v6lZXVwNQX1/PrFmzuOWWWwarBmZmFcVTSNniBCaLWrbA01+CLffCIR+GD90Io95X9NCuxCWX\ny9Ha2sqIESOYNGkSo0eP3vNpA2D16tXMnj17kCpgZlZ5PIWULU5gsqSzHV7+Eaz8NuR3wfE/gKMu\nKbqidETQ2NjI1q1baW1tpaqqirq6OsaOHbvnU0aX9vZ2NmzYwPz58werJmZmFcdTSNniBCYLIg+v\n3wn1l0Hzeqj7dJK8jJ7e49COjg527NjBtm3b6OjoYMSIESUTly7r1q0jn88zc+bMga6JmVnF8hRS\ntjiBqWT5DnjtV7DqCvi/lTDmA3Dqg3DYJ/Y+LJ+nubmZHTt20NTUBEBNTQ11dXWMGjVqr6miYtas\nWQPgBMbM3tU8hZQtTmAqUcsm2HAjrF8CO19NEpe5S+HIc+Cg5L+ss7OT5uZmGhsbaWpqIp/PM3z4\ncCZMmMDYsWMZOXJkny+3evVqAH/zyMze1TyFlC0Vl8BIOgO4GhgG/DwirixzSINj50bYch9suhve\nfDiZNpp4Khz/QzjiL+nMB60trbS0tNDc3ExLSwsAw4YNY8yYMYwePZra2tp9jrYUs2bNGiZPnkxN\nTc2BrpWZWWZ4CilbKiqBkTQM+AnwcWATsELSPRGxqryRHWC73k7uZdn+LGxbAVufSKaIgKiZRsdR\nX6O17jzahh1Ge3s7bevW097evufp1dXVTJgwgdraWmpqavYraSm0Zs0aTx+Z2buep5CypaISGOAk\nYF1EbACQdCtwFjCoCUz93ZdR01IPgAgg2JMiCBRJWRIjf9zec4Y/7os8VbQynBZGRDPVsZUR7Nxz\n5G6NomnETBrHX0TL2FPJ10xPTroNIEdTUxPbt2/nrbfeoqGhgVwux65duw5ofVetWsWCBQsO6DnN\nzLKmK4G5/PLLGTduXJmjqWxnn302p512WlljqLQE5nDg9YL9TcCfFh4gaQGwAGDy5MkDEsTInS8x\n5aCn9uxHmors2Ve3/e6PF+wHB9Ghajr0Xtp0GG8fdCytBx3KTh1Cbvfh5Fpq2bmjhZ07d9LY+Hty\nuQfJ5XI0NDSwceNGWltbB6SOhUaNGsXpp58+4NcxM6tkRx99NHV1dTzwwAPlDqXizZw50wlMf0XE\n9cD1AHPmzBmQicqZ5903EKc1M7MKNnv2bDZv3lzuMKyPiv8wSPlsBiYV7B+RlpmZmZntUWkJzApg\nuqSpkkYA5wL3lDkmMzMzqzAVNYUUER2SLgIeJPka9Q0RsbLMYZmZmVmFqagEBiAi7gfuL3ccZmZm\nVrkqbQrJzMzMbJ+cwJiZmVnmOIExMzOzzHECY2ZmZpnjBMbMzMwyxwmMmZmZZY4TGDMzM8scJzBm\nZmaWOU5gzMzMLHMUMSALOg8KSW8Brw3Q6ScAWwfo3EOR26v/3Gb94/bqP7dZ/7i9+meg2uvIiDhk\nXwdlOoEZSJKeiYg55Y4jK9xe/ec26x+3V/+5zfrH7dU/5W4vTyGZmZlZ5jiBMTMzs8xxAlPa9eUO\nIGPcXv3nNusft1f/uc36x+3VP2VtL98DY2ZmZpnjERgzMzPLHCcwZmZmljlOYIqQdIaklyWtk7Sw\n3PFUOkmvSnpRUr2kZ8odTyWSdIOknKSXCsrGS3pY0h/Sf8eVM8ZKUqK9viVpc9rP6iV9spwxVhJJ\nkyQ9KmmVpJWSvpqWu48V0Ut7uY+VIGmkpKclPZ+22eK0vGx9zPfAdCNpGLAW+DiwCVgBzI+IVWUN\nrIJJehWYExH+AagSJH0EaAZ+ERGz0rL/ALZHxJVpojwuIi4rZ5yVokR7fQtojojvlTO2SiTpMOCw\niHhO0ijgWeBs4PO4j/XQS3t9DvexoiQJqImIZklVwHLgq8BfUaY+5hGYnk4C1kXEhojYBdwKnFXm\nmCzjIuIxYHu34rOAm9Ltm0heQI2S7WUlRMQbEfFcut0ErAYOx32sqF7ay0qIRHO6W5X+BWXsY05g\nejoceL1gfxPu2PsSwG8kPStpQbmDyZCJEfFGuv0mMLGcwWTExZJeSKeYPB1ShKQpwGzgKdzH9qlb\ne4H7WEmShkmqB3LAwxFR1j7mBMYOhJMj4jjgTOAr6fC/9UMkc7mez+3dtcA04DjgDeD75Q2n8kiq\nBX4N/ENENBY+5j7WU5H2ch/rRUR0pq/1RwAnSZrV7fFB7WNOYHraDEwq2D8iLbMSImJz+m8OuJNk\nGs72rSGdi++ak8+VOZ6KFhEN6QtoHvgZ7md7Se9L+DVwc0TckRa7j5VQrL3cx/omIt4GHgXOoIx9\nzAlMTyuA6ZKmShoBnAvcU+aYKpakmvQmOCTVAJ8AXur9WZa6B7gg3b4AuLuMsVS8rhfJ1GdwP9sj\nvcFyCbA6In5Q8JD7WBGl2st9rDRJh0gam25Xk3zRZQ1l7GP+FlIR6VfnrgKGATdExL+XOaSKJWka\nyagLwHDgFrdXT5J+CcwjWX6+AVgE3AXcBkwGXgM+FxG+cZWS7TWPZGg/gFeBLxXMvb+rSToZWAa8\nCOTT4n8lua/DfaybXtprPu5jRUk6luQm3WEkgx+3RcTlkg6mTH3MCYyZmZlljqeQzMzMLHOcwJiZ\nmVnmOIExMzOzzHECY2ZmZpnjBMbMzMwyxwmMmZmZZY4TGDMzM8scJzBmZkOIpGmSlki6vdyxmA0k\nJzBmA0xSSFpasD9c0luS7i1zXL/r5/HfkvS1gYpnIOKQNFbShb08foOknKQePxkv6QxJL0taJ2nh\nO4m523mrJf2vpGEH6pyFImJDRHyx2zVHSHpM0vCBuKZZOTiBMRt4O4FZ6fohkKwhUvYFQiPiz8od\nwyAYC5RMYIAbSRak20uaXPyEZIX19wPzJb3/AMX0BeCOiOh8pyeSdIyke7v9Hdr9uIjYBTwCnPNO\nr2lWKZzAmA2O+4FPpdvzgV92PSDp7yQ9Lale0k+7PplLukvSs5JWSlqQlk2RtFrSz9LyhwoSIwrO\n+XVJl6TbP5T023T7o5JuTreb93VOSd+QtFbScmBGsYqlC3reJ+l5SS9JOic95xpJN6fnvl3Se/dR\n36Ll/YijR3sBVwJ/kp7zu92fExGPAcXWbTkJWJeOZuwCbgXO6na9rjremMZ2s6SPSXpc0h8klVrJ\n+G8pWPBO0n9J+rGk5ZJek3SypP9Mz7mk4LgL0vq9kLYDEfFiRHy621+p1YDvSq9tNiQ4gTEbHLcC\n50oaCRxLssgeko4m+VT84Yg4Dujkj28yX4iIE4A5wCXpomkA04GfRMQHgLeBzxa53jLglHR7DlAr\nqSote6zI8T3OKekEktXYjwM+CZxYom5nAFsi4oMRMQv477R8BnBNRBwNNAIXlqpvb+3QjziKtddC\nYH1EHBcRXy/xvGIOB14v2N+UlnX3PuD7wMz07zzgZOBrJIsD7kXJCvfTIuLVguJjgA0RcTLwU5JV\nkv+ZZOTnU5Leo2TF98uAuRFxLPDpUoFLOljSdcBsSf9S8NBLlG47s8zxfKjZIIiIFyRNIRl9ub/g\nodOAE4AVkgCqga5P0JdI+ky6PYkkyXgTeCUi6tPyZ4EpRS75LHCCpNFAO/AcyRv7KcAlRY4vds4J\nwJ0R0QIg6Z4S1XsR+L6k7wD3RsQySeOA1yPi8fSYpel120rUd3Qv7XBKH+Mo1V4D6ZWIeDGNayXw\nSESEpBcp/v8ygSRBJH3OSJJprqvSogCWdK2ALKkT2EWyAnA1STvfFBHPlAooIrYBXy5S3ilpl6RR\nEdHU/6qaVRYnMGaD5x7ge8A8oGs0RcBNEVH4SRlJ84CPkXzibpH0P8DI9OH2gkM7Sd7Y9hIRuyW9\nAnwe+B3wAnAqyYjB6iKx7fOcpUTEWknHk4yOfFvSI8AvSN6M9zqU0vW9uFh5X+2jvfbHZpIkqMsR\nFL9vqbDd8gX7eYq/vrZ2i+sDwHMRkU/3PwhcCyDpCJKRrQBaJM0C/gK4XtLPI+Ka/lUJgPeQJJFm\nmecpJLPBcwOwuOsTe+oR4K+7bryUNF7SkcAYYEf6ZjwT+NB+XG8ZyVTGY+n2l4Hfp2+IffEYcLaS\nb82MInnz7EFSHdASEUuB7wLHpw9NljQ33T4PWE7p+pYq72scpdqrCRjVx/oWWgFMlzQ1nfY5lyQB\nfUciYgcwLB15gWT66PmCQ44lSTYhSWZeAJA0PSJ2RsStwL3sR3KWTqltjYjd+xu/WSVxAmM2SCJi\nU0T8qFvZKuDfgIckvQA8DBxGch/JcEmrSW5EfXI/LrksPdcTEdFA8sl7WT/ifQ74Fckb7AMkb+rF\nHAM8LakeWAR8Oy1/GfhKWodxwLWl6ttLO/Q1jqLtlU6nPJ7eXNzjJl5JvwSeAGZI2iTpi+nzOoCL\ngAdJRqxui4iVfWi2vniI5D4ZSNquPo1lJFCdJjmwdzLzDSVf6X4OmArsz+jLqcB9+x21WYVR3z+M\nmZn1TXq/z73pTb1WIJ1u+8eI+PtBvu4dwMKIWDuY1zUbKB6BMTMbROmI0qMaoB+yKyadBrvLyYsN\nJR6BMTMzs8zxCIyZmZlljhMYMzMzyxwnMGZmZpY5TmDMzMwsc5zAmJmZWeY4gTEzM7PMcQJjZmZm\nmeMExszMzDLn/wEzlsZcC+U35gAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# now accounting for power losses arising from the transmission and distribution of electricity\n", "\n", "fig = plt.figure(figsize=(9,6))\n", "\n", "plt.plot(wind_10m[ind1], P[ind1], color='k')\n", "\n", "Rd10 = rayleigh_dist_power_curve(wind_10m, P, height=10)\n", "plt.plot(wind_10m[ind1], Rd10[ind1]*0.725, color='LightGrey') # <-- these seem to be too low ??\n", "\n", "Rd60 = rayleigh_dist_power_curve(wind_10m, P, height=60)\n", "plt.plot(wind_10m[ind1], Rd60[ind1]*0.725, color='orange') # <-- these seem to be too low ??\n", "\n", "plt.xlabel('Mean wind speed at 10 m ($ms^{-1}$)')\n", "plt.ylabel('Load factor (%)')" ] } ], "metadata": { "kernelspec": { "display_name": "Python 2", "language": "python", "name": "python2" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 2 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython2", "version": "2.7.13" } }, "nbformat": 4, "nbformat_minor": 2 }