{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Ch-8, Hydro-Electric Plants" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## example 8.1 Page 134" ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "\n", "power developed by hydro plant is 176.52 MW\n" ] } ], "source": [ "h=100 #given height\n", "q=200 #discharge\n", "e=0.9 #efficiency\n", "p=(735.5/75)*q*h*e\n", "print \"\\npower developed by hydro plant is %0.2f MW\"%(p/1e3)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## example 8.2 page 134" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "\n", "average rate of inflow is 1000cu-m/sec\n", "\n", "power developed is 1250.35 MW\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY0AAAEZCAYAAABrUHmEAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xm8XfO9//HXRwiRhKBEIkhEDCGEyIC6TlRJW1OLalWN\nddvSGu69lGgJrTad/Mot2hLEWONVQVSaOGpMiIQQmUgqcyRSiYiMn98f37XZOc7J2eecvfYa9vv5\neJxH1l577bU+OznZn72+08fcHRERkVJsknQAIiKSHUoaIiJSMiUNEREpmZKGiIiUTElDRERKpqQh\nIiIlU9KQXDOzWWb2pbScpxLMrNbMzkk6DsknJQ3JO49+0nKeSshSrJIxShoiLWRmrcp4rk3LdS6R\nOChpSDU4wMxeN7N/m9lfzWxzM3vTzI4pHGBmm5nZYjPbP3r8XTP7V7RvcPHJzGyImT1kZneZ2YfA\nGWbW2cweM7MlZjbdzL5XdHwbMxtuZh+Y2WQzu9TMZhc9Pyva9waw3MxamdllZjbDzJaZ2VtmdkLR\n8Wea2Qtm9r/Re3rbzI6o8567mtnz0ev/bmbblfnvVKqUkobknQEnA0cD3YD9gDOB4cBpRcd9FZjr\n7q+bWU/gJuA7QGdgO6BLnfMeBzzo7lsD9wJ/Bd4DOgEnAb80s4HRsVcBu0TX/3J03brNR98CvgJ0\ncPd1wAzgi+6+FXA1cLeZdSw6vl90zHbR+R8xsw5F7/nU6H3uALQG/qfxvyqRxilpSN45cIO7L3D3\npcAIoDdwN/A1M2sXHfdd4K5o+yRghLs/7+6rgZ8B6+uc90V3fyza3h44BPiJu69299eBW4HTo+dP\nBn7p7h+6+1zgesIHe90Y57r7KgB3f8jdF0TbDwDTgf5Fr1nk7te7+7ro+anAMUXnu83dZ7j7J8AD\n0XsWaTElDakGC4q2VwJt3X0+8AJwUvQNfRBwT3RMJ2BO4QXu/jGwpM455xRtdwY+cPcVRfvei/YX\nnp9d9FzxawuKn8fMTjezCWa21MyWAvsS7ioK5tZ5/b+iuAvqvud2iJSBkoZUo8K3/EIT1cmEO4f5\n0f75wM6fHmy2JRt+YMOGzUvzgG2L7logNEcVPtg3OF+d7c+dz8x2Bf4CnA9s6+7bAG+y4d3JTnVe\nv2sUh0islDSkmj0KHAhcANxZtP8h4BgzO9TMWgPXsJH/K+4+G3gR+FXUyb4fcDahCQxC89DlZtbB\nzHYCfsTGh8S2jZ5fDGxiZmcR7jSK7WBmF0Qd+CcDewFPFj1viMRASUOqzadzGNx9JfAI0DX6k2j/\nZMK3/HsJ394/YMPmo/rmQXw7Os+86FxXuvuY6LlrCE1SM4GngQeB1Q0GGK7/e+AlQjPTvsDzdQ4b\nC/QA3gd+DpwY9dkUx7ixeEWaxeIqwmRmOxO+ve1A+IX9i7vfYGZDgO8RftkBBrv7yOg1lxO+oa0D\nLnD3p6P9fYA7gC2AJ939wliClqpjZj8Derj76Y0eXL5r/hD4prsPbPTg+l9/JnCOux9W1sBEShDn\nncYa4GJ33wcYAJxvZnsTEsh17n5A9FNIGD2BU4CehE7Jm8yscIt9M+E/SQ+gh5kNijFuqRJmti3h\nS8pfYr7OjlFT1yZmtifwX8D/xXlNkbjEljSiIY4To+2PgLf5rPOuvvbW44H73H2Nu88ijEHvb2ad\ngPbuPi467k7ghHpeL1IyMzuXMMJppLvXbfopt9bAn4BlwGhCX8pNLTifmpskMRXp0zCzrsABwMvR\nrh9HM3SHFU1I6syGQxHnEJJM3f1z+fzIEZEmcfdb3L2du59XgWu95+69out1cfdL3H1tC8433N3/\no5wxipQq9qQRDUN8CLgwuuO4mTAztjdhKOLv445BRETKI9bF0cxsM+Bh4G53fxTA3RcVPX8rYYYu\nhDuI4vHrXQh3GHPZcAmHLnx+YhNmptt1EZFmcPeSh2jHdqcRdWIPAya7+x+K9hfPWv06MCnafgz4\nlpm1NrNuhOGE46KlFJaZWf/onN8ltAl/jrvn9ueqq65KPAa9P72/antv1fD+mirOO41DCbNt3zCz\nCdG+wcC3zaw3oSNvJvB9CGPTzewBYDKwFjjPP3tH5xGG3LYhDLl9Ksa4RUSkAbElDQ8jUuq7kxm5\nkdf8EvhlPfvHA73KF52IiDSHZoRnRE1NTdIhxErvL7vy/N4g/++vqWKbEV5pZuZ5eS8iIpViZnga\nOsJFRCR/lDRERKRkKmIvkiMffwzdusG//125a+62G0yeDKbF2KuCkoZIjjz/POy+O4wZ0/ix5dKt\nG8yaFf6U/FPSEMmR0aPhyCNh880rd82DD4aXXlLSqBbq0xDJkTFj4Etfquw1C0lDqoOShkhOLF0K\nU6bAgAGVve4hhyhpVBMlDZGcqK0NH+CtW1f2ugceCG+/HTrhJf+UNERyIommKYAttoBeveDVVyt/\nbak8JQ2RnBg9OpmkAaFf48UXk7m2VJaShkgOzJsHCxZA797JXF+d4dVDSUMkB555BmpqoFWrZK5f\nSBpa/i3/lDREciDJpimAnXcOHfDvvptcDFIZShoiGeceksYRRyQbh4beVgclDZGMe/ddWLMG9tor\n2TjUr1EdlDREMq7QNJX0goFKGtVBSUMk49LQNAVwwAEwdSp89FHSkUiclDREMmz9+jByKslO8ILN\nN4f994dXXkk6EomTkoZIhr35Jmy9NeyyS9KRBGqiyj8lDZEMS0vTVIGSRv4paYhkWFLrTTXk4IPh\n5Zc1yS/PlDREMmrtWvjnP2HgwKQj+cxOO8GWW8KMGUlHInFR0hDJqFdega5dYfvtk45kQ2qiyjcl\nDZGMSlvTVIFWvM03JQ2RjEp6vamG6E4j38xz0mNlZp6X9yLSmJUrQ7PUvHmw1VZJR7Oh1ath221h\n/nxo3z7paKQxZoa7l7yegO40RDLoxRdhv/3SlzAgrHbbuzeMG5d0JBIHJQ2RDEpr01SBmqjyS0lD\nJIPSNqmvLi2Tnl/q0xDJmA8/hC5d4P33YYstko6mfgsWQM+esHgxbKKvpqmmPg2RnHv2WRgwIL0J\nA2DHHcOaWNOmJR2JlJuShkjGpL1pqkD9GvmkpCGSMWmd1FeXkkY+KWmIZMjChTBnDvTpk3QkjVPS\nyKfYkoaZ7Wxmz5jZW2b2ppldEO3f1sxGmdk0M3vazDoUveZyM5tuZlPM7Kii/X3MbFL03PVxxSyS\ndmPGwOGHQ6tWSUfSuP33h5kzQ8e95EecdxprgIvdfR9gAHC+me0NXAaMcvc9gNHRY8ysJ3AK0BMY\nBNxk9mnV45uBc9y9B9DDzAbFGLdIao0Zk43+DIDNNgt3RJrkly+xJQ13X+DuE6Ptj4C3gZ2A44Dh\n0WHDgROi7eOB+9x9jbvPAmYA/c2sE9De3Qu/encWvUakqqR9Ul9dWrwwfyrSp2FmXYEDgLFAR3df\nGD21EOgYbXcG5hS9bA4hydTdPzfaL1JVZs6Ejz8O8x+yQv0a+bNp3Bcws3bAw8CF7r78sxYncHc3\ns7LNyBsyZMin2zU1NdTU1JTr1CKJKzRNWcnTsJJ38MFw5pmwfr0m+aVFbW0ttbW1zX59rDPCzWwz\n4HFgpLv/Ido3Bahx9wVR09Mz7r6XmV0G4O5Do+OeAq4C/hUds3e0/9vA4e7+gzrX0oxwybVTTw1N\nU+eck3QkTdO9O4wYka07pGqSmhnhUSf2MGByIWFEHgPOiLbPAB4t2v8tM2ttZt2AHsA4d18ALDOz\n/tE5v1v0GpGq4J6d+Rl1qYkqX+K8YTwUOA0YaGYTop9BwFDgy2Y2DTgieoy7TwYeACYDI4Hzim4d\nzgNuBaYDM9z9qRjjFkmdyZND7e2uXZOOpOmUNPJFCxaKZMANN8CkSXDLLUlH0nQTJsBpp8FbbyUd\nidQnNc1TIlI+WRtqW6xXL3jvPVi6NOlIpByUNERSbu1a+Oc/YeDApCNpnk03hYMOgrFjk45EykFJ\nQyTlXnst1M/o2LHxY9NK/Rr5oaQhknJZbpoqUNLIDyUNkZTL0npTDRkwIKxBtX590pFISylpiKTY\nJ5/Ayy+HlW2zbPvtw8/kyUlHIi2lpCGSYi+9BPvsE0qnZt0hh6iJKg+UNERSLA9NUwVa8TYflDRE\nUiwPneAF6gzPB80IF0mpZcugc2d4/31o0ybpaFpu3TrYZhuYNQu23TbpaKRAM8JFcuK556Bfv3wk\nDAglavv2DR37kl1KGiIplaemqQI1UWWfkoZISmV1KfSNUdLIPvVpiKTQ++9Djx6weHFYuykvliyB\nbt3C4oWtWiUdjYD6NERy4Zln4LDD8pUwALbbLnTuv/lm0pFIcylpiKRQHpumCtRElW1KGiIpNHp0\nfib11aWkkW1KGiIp89578OGHsO++SUcSDyWNbFPSEEmZwtIhm+T0f2fPnrBwYejkl+zJ6a+lSHbl\nuWkKwqipfv00yS+rlDREUsQ9n5P66tKKt9mlpCGSIlOnQuvWsNtuSUcSL614m11KGiIpUmiaspKn\nWmVT//7w6quwdm3SkUhTKWmIpEg1NE1BWO12551h0qSkI5GmajRpmFk7M2tV9LiVmbWNNyyR6rNu\nHdTW5rsTvJiG3mZTKXcao4HixZm3BEbFE45I9Zo4EXbcETp1SjqSylDSyKZSksbm7v5R4YG7Lyck\nDhEpo2ppmipQ0simUpLGCjPrU3hgZgcBK+MLSaQ65Xm9qfrsvXdY9XbRoqQjkaYoJWlcBDxoZs+b\n2fPA/cCP4w1LpLqsXh2GoB5+eNKRVM4mm4RRVLrbyJZGF15291fMbE9gT8CAKe6+JvbIRKrIyy/D\nnnuGUUXVpNBEdfzxSUcipSpl9FRb4DLgQnefBHQ1s2Nij0ykilRb01SB+jWyp5TmqduB1cAh0eN5\nwLWxRSRShfK+3lRD+veH8eNhjdouMqOUpNHd3X9NSBy4+4p4QxKpLh99FIbbfvGLSUdSeVtvDV27\nwhtvJB2JlKqUpLHKzD6dp2Fm3YFV8YUkUl2efx769IEtq3Qgu5qosqWUpDEEeAroYmb3AmOAn8QZ\nlEg1qdamqQKteJstjSYNd38aOBE4C7gX6OPuz5RycjO7zcwWmtmkon1DzGyOmU2Ifr5S9NzlZjbd\nzKaY2VFF+/uY2aToueub8gZF0q7aJvXVpRVvs6WU0VNfBD5x98eBbYDBZrZriee/HRhUZ58D17n7\nAdHPyOg6PYFTgJ7Ra24y+3Stz5uBc9y9B9DDzOqeUySTliyBGTNCUaJqtcceobztggVJRyKlKKV5\n6mbgYzPbH/gv4B3gzlJO7u7PAUvreaq+hZ+PB+5z9zXuPguYAfQ3s05Ae3cfFx13J3BCKdcXSbva\n2tABvtlmSUeSnE02gQED1ESVFaUkjbXuvp7wQX2ju98ItG/hdX9sZq+b2TAz6xDt6wzMKTpmDrBT\nPfvnRvtFMq/am6YK1BmeHaUkjeVmNhg4DXg8Wia9Jd+Lbga6Ab2B+cDvW3AukUyr1kl9dSlpZEej\ny4gQ+hlOBc529wVmtgvw2+Ze0N0/XZ7MzG4FRkQP5wI7Fx3ahXCHMTfaLt4/t75zDxky5NPtmpoa\nampqmhumSOzmzoXFi2G//ZKOJHn9+sGECWENrtatk44m32pra6mtrW32683d63/C7O+EobYj3X1K\nsy9g1hUY4e69osed3H1+tH0x0NfdT406wu8F+hGan/4B7O7ubmZjgQuAccATwA3u/lSd63hD70Uk\nje68E0aMgAcfTDqSdNh/f7j1VujbN+lIqouZ4e4lFxje2J3GmYRRTEOiBQvHAiOBf5Q6K9zM7gMO\nB75gZrOBq4AaM+tNGEU1E/g+gLtPNrMHgMnAWuC8oixwHnAHoRjUk3UThkgWqWlqQ4Wht0oa6dbg\nncYGB4V+jP7AV4AjgE+Av7v7b+INr3S605AscYdddgkd4XvskXQ06TB8OIwcCX/9a9KRVJem3mmU\n0hGOu69z9xfd/Wfufiihn6PefgURadz06eHPHj2SjSNN1BmeDY12hJvZboSiS12Ljnd3Py7GuERy\nrdA0ZSV/v8u/Hj1gxQqYNw86d046GmlIKaOnHgUKo5zWR/vUDiTSAqNHw7HHJh1Fuph9NsnvxBOT\njkYa0mifhpmNc/fUL3KgPg3JivXrYYcd4PXXYSdNU93AtdfC0qXwu98lHUn1iKNP43+jRQYPNrMD\nCz8tiFGkqr3xBnzhC0oY9dGKt+lXSvPUPsB3gYF81jxF9FhEmqjal0LfmL59Q0GqVatg882Tjkbq\nU0rSOBno5u6r4w5GpBqMHg3nnJN0FOnUrl0YgjxhQujfkPQppXlqEmFJdBFpoTVrQqU+rXDTMA29\nTbdSksY2wBQze9rMRkQ/j8UdmEgejRsHu+8O222XdCTppaSRbqU0T11VtO2EWhgapiTSDFoKvXEH\nHwyDBycdhTSklHKvtYUfYKto+9n4QxPJH6031bju3eGTT2DOnMaPlcoraRmRItfEEoVIFfj4Yxg/\nHg47LOlI0s1MQ2/TrKlJQ0Sa6fnnoXdvaNs26UjSr7DiraRPU5PG92OJQqQKqGmqdOoMT69SlhE5\no84uB3D3O+MKqjm0jIikXd++cN11ap4qxYoVYamVJUtgiy2SjibfylmEqaAvn42W2gL4EvAakKqk\nIZJmS5fC1KnQv3/SkWRD27aw117w2muhf0PSo9Gk4e4/Kn5sZh2A+2OLSCSHnn02fPip/nXpCk1U\nShrp0pyO8I+BbuUORCTPND+j6dSvkU6lFGEaUfRwE6An8EBsEYnk0OjRcNddSUeRLYccApdeGkrj\nqlhVepTSEV5T9HAt8C93nx1nUM2hjnBJq/nzYd99YdEiaNUq6Wiywx06dYKxY2HXXZOOJr/i6Ah/\nD5jv7iujC7Qxs67uPquZMYpUlTFjwgKFShhNY/ZZE5WSRnqU0qfxILCu6PF64KF4whHJH9XPaD71\na6RPKUmjVXEtDXdfBWwWX0gi+bF+PYwapU7w5lLSSJ9SksZiMzu+8CDaXhxfSCL5MWZMKO26555J\nR5JNBx0Eb70FK1cmHYkUlJI0fgAMNrPZZjYbuAwtJyJSkltvhe99T6N/mqtNG+jZMyz0KOnQ6Oip\nTw80aw/g7stjjaiZNHpK0mbx4lBwaeZM2Ea1L5vtwguhSxe45JKkI8mnpo6eKnlyn7svT2vCEEmj\nu++GY49VwmgprXibLloaXSQG7p81TUnLFDrD1ZCQDkoaIjEYOxZWrYL/+I+kI8m+XXaBTTaBWbOS\njkSgtMl9mNmhQNei4z1tS6OLpIk6wMuneJJfN616l7hS1p66G9gNmMiGk/yUNETqsXw5PPwwvP12\n0pHkRyFpnHpq0pFIKXcafYCeGpokUpr77w/Lhuy4Y9KR5MfBB8Nf/5p0FAKl9Wm8CXSKOxCRvFAH\nePn16RPu3FasSDoSKeVOY3tgspmNA1ZF+9zdj4svLJFsmjQJ5syBo49OOpJ82WIL6NULXn0VDj88\n6WiqWylJY0jcQYjkxbBhcNZZsGlJQ0ykKQr9GkoaySql3GttBeIQybxPPgkT+l55JelI8ungg+Ge\ne5KOQhrs0zCzF6I/PzKz5XV+lpVycjO7zcwWmtmkon3bmtkoM5tmZk9HNccLz11uZtPNbIqZHVW0\nv4+ZTYqeu755b1UkXo8+CgccoGGhcdEkv3RoMGm4+6HRn+3cvX2dn61KPP/twKA6+y4DRrn7HsDo\n6DFm1hM4hVBOdhBwk9mno9xvBs5x9x5ADzOre06RxKkDPF477wytW8O77yYdSXWLdUa4uz8HLK2z\n+zhgeLQ9HDgh2j4euM/d10RVAWcA/c2sE9De3cdFx91Z9BqRVHj3XXj9dThBv5mxUn2N5CWxjEhH\nd18YbS8EOkbbnYE5RcfNAXaqZ//caL9Iatx2G5x2Gmy+edKR5NshhyhplNMnnzT9NYmuPRVNGMxc\nC+Udd0D//vDaa0lHImmwdi3cfjucc07SkeTfYYfB3/8e/s6l5QYPbvprkhgYuNDMdnT3BVHT06Jo\n/1xg56LjuhDuMOZG28X759Z34iFDhny6XVNTQ01NTfmiLnL33WHM+Fe+AmecAUOGwJZbxnIpyYCn\nngqL6u27b9KR5F+fPtC5cxhFdcYZSUeTTbW1tdTW1rJsGdx8c9Nf32gRJjM7ERhKaEYqdEx7qZ3h\nZtYVGOHuvaLHvwGWuPuvzewyoIO7XxZ1hN8L9CM0P/0D2N3d3czGAhcA44AngBvc/ak616nISidL\nl8Kuu8L8+WF26kUXhRVN//IX1YGuViecEOpm6E6jMp57Dk4/HaZODR3j0jzf/36o9fLrXzetCBPu\nvtEf4B1g78aOa+C19wHzgNXAbOAsYFtCQpgGPE1IGoXjBxM6wKcARxft7wNMip67oYFreSXcdZf7\nccdtuO/xx9132cX9zDPdlyypSBiSEvPmuXfo4L58edKRVJdBg9xvvDHpKLJr+nT37bYLn1fRZ2fJ\nn+ul3Gm84NHw2zSr1J3GSSfB174WZv0WW74crrgCHnwQ/vAH+OY3tSx2NRg6FN55B265JelIqsv4\n8XDccTB9upqGm+O002DPPeFnP2t6uddSksb1wI7Ao4Q7BgiZ6ZFmRxyDSiSNlSvDyqUzZsD229d/\nzMsvh7H63brBTTeFseWST+7Qo0doX+/fP+loqs9JJ4W/d9UOb5o33wxN6TNmQPv28dQI3xpYCRwF\nHBP9HNu8cLNt9Gjo3bvhhAEwYEAYVdWvHxx4INx4I6xfX7kYpXKefRbatAn/1lJ511wDv/0tLCtp\nfQopuPJKuPTSkDCao9E7jayoxJ3G974H++wDF19c2vFvvw3nnhuSxi23hNdKfpx2GvTtCxdemHQk\n1ev006F7d7jqqqQjyYZXXoGvfz0067VpE/bF0Tx1e51dhV7ns5sWbrziThrr1kGnTqH5abfdSn/d\n+vXw5z+H7H7++XD55ZoAlgdLl4YmyHfege22Szqa6vXuu+FOb+pU/TuU4uijQ9L4wQ8+2xdH89QT\nwOPRz2hCc1XVlUJ56aWQNJqSMAA22QR++EOYMCH8HHAAvPhiPDFK5dxzT5inow+qZO22G5x8Mvz6\n10lHkn7PPhvuMM5u4df9JjdPmdkmwAvufnDLLl1ecd9p/M//QNu2cPXVzT+He6gdfcEF8I1vwC9/\nCVuVuvSjpIZ76Nu67jrNzUmDuXNhv/1CAazOnZOOJp3cw2z6738fvvvdDZ+L406jrj0I1fyqhntY\n9rqli9GZhREfb70V1nzZZx8YMaI8MUrljB8fhlgPHJh0JAKw005hCPy11yYdSXo99VRoUj311Jaf\nq5Q+jY/4bH0oJywyeJm7P9zyy5dPnHcab74Z5mbMmlXeuRfPPAP/+Z9hlNUNN0DHjo2/RpL3gx+E\nodRXXJF0JFKweHGYdzB+PHTtmnQ06bJ+PRx0UPh9PfHEzz9f1juNqJ5FTy+qo+HuPdKWMOJWuMso\n92S9gQPhjTdCu2yvXmHRu5wMZsutFSvggQfgzDOTjkSKfeELYaBJS5qP8+qRR0Lf6je+UZ7zbfRO\nI0oab3i0blSaxXmncdBBYTx4nM0REyeGIb1bbx3WserePb5rSfPdcUfol1KzYvp8+GGYbPnss7D3\n3klHkw7r1oUvpNddB4MaKF1X1juN6FP4NTOr2ulLs2fDzJmhEylOvXuH4bxf/WqY5fqb32j55zRS\ndb702npr+O//1pyNYvfcE0b4HX10+c5ZSp/GVGB34F98NtTW3X2/8oXRcnHdafzxj2FCzPDhjR9b\nLu++G9rNFy8OH1IHHli5a0vD3n4bjjgC3nsPNtss6WikPitWhLuNJ54Iw9ur2erVoZ/nzjs3/qU3\njtFTRwPdgSMIy4ccSyjZWhXKMWqqqXbbLRSaueiiMBfgkkvg448rG4N83rBhoS9DCSO92rYNhYV+\n+tOkI0nesGEhaZS7lUTLiGxEce2Mtm3LeuqSLVoUlqkYN041O5K0ejV06QIvvBC+yUp6rVoVPizv\nvTeUh61GK1fC7rvD3/4W+mQ3phLzNKrGE0+Ezu+kEgbADjvAffeFIblnnx3Go3/wQXLxVKvHHgvz\napQw0m/zzcOyPYMHV+9oxBtvDIunNpYwmkNJYyOSaJpqyNe+FuaLtG8fPrzuv796/0MkQR3g2XL6\n6aGF4B//SDqSylu2LIz2vOaaeM6v5qkGlFI7Iymq2VFZ//pXGIwwZ85nK4NK+t1/P/z+96EcczUV\nRLvmmrDG1F13lXa8mqfKpJTaGUlRzY7Kuv32sPyCEka2nHxy6It67LGkI6mcJUtCU/aQIfFdQ3ca\nDWhq7YykqGZHvNatC3d0I0bA/vsnHY001eOPh3IEEydCq1ZJRxO/n/wkTHL8059Kf43uNMpg3brw\n7eT445OOpHF77w3//GdYubKmJnzDWLUq6ajyY9SosCaYEkY2fe1r0K5daKrKu/nzQ99b3MONlTTq\n0dzaGUlRzY74qAM828zC6rdXXQVr1iQdTbyuvTbMI+rSJd7rqHmqHuWonZEU1ewon0WLYI89wgxw\n/R1m25FHwimnhKbcPJo1C/r0gSlTmt4Pq+apFipX7Yyk1K3Zse++oV1Xmu7OO0NpTCWM7Lv2Wvj5\nz8P/iTy65ho477zKDNzRnUYdcdXOSEqhZkefPnD99arZUSr30F80bBgcemjS0Ug5HHdcWDvsoouS\njqS8pkwJS4VMnw4dOjT99brTaKG4amckpVCzo1u3UBJTNTtK88IL4XegWpehyKNf/AKGDoWPPko6\nkvK66qqwum9zEkZzKGnUkeWmqYa0aQO/+lVYBPHGG+HLX4Z33kk6qnQrdIDn5cuDhC9NAweGeQx5\nMXFiGD354x9X7ppqnioye3aY0LdwIWy6aZkCS5m1a0Mz1a9+FcZ0X3xxft9rc334YVioctq0sPaX\n5Me0aeHucfp02GabpKNpuWOOCbUyWpI01DzVAn/7W/hHyPOH6KabhlvZcePCHIR+/cLscvnMffeF\nuzEljPzZY4/QkvC73yUdScu99BJMmhT6LCtJSaNIHpumGlK3Zsell6pmR4HmZuTblVeGGdMLFyYd\nSfO5h1W4DRNQAAAMs0lEQVR8r7wyrOpbSUoakaVLw7fvo45KOpLKMQurgU6aFBbj22+/sOZWNZsw\nAd5/P4zrl3zaZRf4zndCE21WjR4Nc+fCGWdU/tpKGpE01M5Iyg47hII1118fanacfXb11uwYNiy8\n/2pYp6iaDR4cVoGdPTvpSJrOHa64IszNSKIpXUkjUk1NUw0p1Oxo1y5MCnzggeoanrtyZejPOOus\npCORuO24Y+gL+PnPk46k6UaMCJMUv/nNZK6v0VOku3ZGUgo1O3bbLQzTrYaaHXffDffcAyNHJh2J\nVMIHH4SO8Zdeyk5FxvXrwwjPa6+FY48tzzk1eqoZ0lw7IymFmh19+1ZPzQ51gFeXbbcNA0HirD1R\nbvffD1tuGUZ5JkV3GmSndkZSCjU73EPNjp49k46o/KZNC0sxzJ4NrVsnHY1UyvLl4S5j1Cjo1Svp\naDZuzZrwf+9Pf4Ivfal859WdRhNlqXZGUgo1O047DQ4/PKz+m7eaHbfdFkaSKWFUl/btw3Dzn/0s\n6UgaN3x4GPlVzoTRHIklDTObZWZvmNkEMxsX7dvWzEaZ2TQze9rMOhQdf7mZTTezKWZWtoGxWaud\nkZTimh2vvRaarPJSs2PNmvAf8pxzko5EkvDDH8L48WHIfVqtWhVGS117bdKRJHun4UCNux/g7v2i\nfZcBo9x9D2B09Bgz6wmcAvQEBgE3mVlZYteoqabp0iX8nV19dViC/Uc/gmXLko6qZZ54AnbfHfba\nK+lIJAlt2oRqd3FXvGuJP/85VI8cMCDpSJJvnqrbjnYcMDzaHg4UPs6PB+5z9zXuPguYAfSjhbJe\nOyMpeavZoQ5wOfvssIhnbW3SkXzeihVhImJahgcnfafxDzN71cwK9bQ6unthcv9CoFD9oTMwp+i1\nc4CdWhrAW2+FponevVt6puq0zTbhA3f48DCI4Fvfyt7SDHPmhGa2k05KOhJJ0mabhbvnK65I39yk\nG24IfYlp+ZxKcmm+Q919vpltD4wysynFT7q7m9nG/vk+99yQorFzNTU11NTUbDSAvNXOSEqhZsc1\n14SlSIYODbWKs/D3escdoQxoNa4EIBv69rfDN/qRI+GrX006muDf/4brroPnny/fOWtra6ltwS1V\nKobcmtlVwEfAuYR+jgVm1gl4xt33MrPLANx9aHT8U8BV7j626BxNHnJ70EHw29+GDz0pj4kTQ1NP\nhw6hHbZ796Qjatj69SG+hx4KlQ1FHnkkNAONHx8GfyTtpz+FefPC6L64ZGLIrZltaWbto+22wFHA\nJOAxoLAE1xnAo9H2Y8C3zKy1mXUDegAtGuswezbMnBnG5kv59O4dZpN/5SvQv39IymvXJh1V/caM\nCcntwAOTjkTS4utfD+uOPfxw0pHAokVw882hMl+aJJVLOwLPmdlEYCzwuLs/DQwFvmxm04Ajose4\n+2TgAWAyMBI4r6UVl6qhdkZSslKzQ9X5pC6zMKz1yiuT/7IzdCicemooCJYmqWieKoemNk8deSSc\nf374ZiHxcQ+riV5ySVjGeciQsAxC0hYvDsNsZ87MRwU3KR/30PF89tmhby4JhVIFkyeHdfHilInm\nqaRVY+2MpKS1Zsfdd4cF35QwpK7C3cbVV8Pq1cnE8POfh1V4404YzVGVSePJJ6u3dkZS0lSzwz00\nTZ17buPHSnU67LAw2fPWWyt/7RkzQp/KpZdW/tqlqMqkoQl9yUlDzY6xY8M3SA2CkI35xS/CT6XL\nIA8ZAhdcEFbhTaOq69P45BPo2FG1M9KgULOjW7fKjot/9NGw6Ftav8lJepx0Emy1VSgRUAkrVoQR\nhzNmhMUUK6GpfRpVlzSeeAJ+8xt49tkKBCWNWr0a/vhHmD69ctfcYouwqmlav8lJesyaBb/7XVgN\nu1KOPbayX6KUNBpx7rlhTXrVzhARUdLY6DHr1kHnzmE5dC2FLiKiIbcb9fLLYQibEoaISPNUVdLQ\nqCkRkZapmqThDv/3f0oaIiItUTVJY/Jk1c4QEWmpqkkaqp0hItJyVZc0RESk+apiyO3s2aFZauFC\nLYUuIlJMQ27r8dhjqp0hIlIOVZE01DQlIlIeuW+eWro0VL6aP19LoYuI1KXmqTpUO0NEpHxynzTU\nNCUiUj65bp5S7QwRkY1T81SR0aPDUFslDBGR8sh10lDTlIhIeeW2eUq1M0REGqfmqYhqZ4iIlF9u\nk4aapkREyi+XSUO1M0RE4pHLpKHaGSIi8chl0lDtDBGReOQ6aYiISHnlbsitameIiJSu6ofcqnaG\niEh8cpc01DQlIhKfXDVPffCBq3aGiEgTVHXzlGpniIjEK1dJQ01TIiLxykzzlJkNAv4AtAJudfdf\n13net9rKVTtDRKQJctk8ZWatgD8Cg4CewLfNbO+6x+W5dkZtbW3SIcRK7y+78vzeIP/vr6kykTSA\nfsAMd5/l7muAvwLH1z0oz01Tef/F1fvLrjy/N8j/+2uqrCSNnYDZRY/nRPs2cPzn0oiIiJRTVpJG\nSR0vqp0hIhKvTHSEm9kAYIi7D4oeXw6sL+4MN7P0vxERkRRqSkd4VpLGpsBU4EvAPGAc8G13fzvR\nwEREqkwmVmhy97Vm9iPg74Qht8OUMEREKi8TdxoiIpIOWekI3ygzG2RmU8xsupn9JOl4ysnMdjaz\nZ8zsLTN708wuSDqmcjOzVmY2wcxGJB1LuZlZBzN7yMzeNrPJUf9cbpjZ5dHv5iQzu9fMNk86ppYw\ns9vMbKGZTSrat62ZjTKzaWb2tJl1SDLGlmjg/f02+v183cweMbOtN3aOzCeNUif+Zdga4GJ33wcY\nAJyfs/cHcCEwmRJHyWXM9cCT7r43sB+Qm2ZVM+sKnAsc6O69CE3H30oypjK4nfBZUuwyYJS77wGM\njh5nVX3v72lgH3ffH5gGXL6xE2Q+aVDixL+scvcF7j4x2v6I8KHTOdmoysfMugBfBW4FclWgN/rG\ndpi73wahb87dP0w4rHJaRvhSs2U0WGVLYG6yIbWMuz8HLK2z+zhgeLQ9HMjsNOL63p+7j3L39dHD\nsUCXjZ0jD0mjpIl/eRB9szuA8A+bF/8PuARY39iBGdQNeN/Mbjez18zsFjPbMumgysXdPwB+D7xH\nGNX4b3f/R7JRxaKjuy+MthcCHZMMJmZnA09u7IA8JI08Nml8jpm1Ax4CLozuODLPzI4BFrn7BHJ2\nlxHZFDgQuMndDwRWkO2mjQ2YWXfgIqAr4e63nZl9J9GgYuZh5FAuP3PM7Apgtbvfu7Hj8pA05gI7\nFz3emXC3kRtmthnwMHC3uz+adDxldAhwnJnNBO4DjjCzOxOOqZzmAHPc/ZXo8UOEJJIXBwEvuvsS\nd18LPEL4N82bhWa2I4CZdQIWJRxP2ZnZmYRm4kaTfh6SxqtADzPramatgVOAxxKOqWzMzIBhwGR3\n/0PS8ZSTuw92953dvRuhA3WMu5+edFzl4u4LgNlmtke060jgrQRDKrcpwAAzaxP9nh5JGNCQN48B\nZ0TbZwB5+uJWKDtxCXC8u3/S2PGZTxrRN5zCxL/JwP05m/h3KHAaMDAaljoh+kfOozze9v8YuMfM\nXieMnvplwvGUjbu/DtxJ+OL2RrT7L8lF1HJmdh/wIrCnmc02s7OAocCXzWwacET0OJPqeX9nA/8L\ntANGRZ8vN230HJrcJyIipcr8nYaIiFSOkoaIiJRMSUNEREqmpCEiIiVT0hARkZIpaYiISMmUNEQq\nzMy2NrMfFj2uyeOy8JJPShoilbcNcF7SQYg0h5KGyEZEy9NMiVaqnWpm95jZUWb2QlSUp29UpOfR\nqIjNS2bWK3rtkKjozTNm9o6Z/Tg67VCgezT79jeEmfDtzOzBqBjO3UXXHxoVOXrdzH5b+b8BkQ1l\noka4SMK6AycSlql5BTjF3Q81s+OAwYSl+ce7+wlmNpCwtMYB0Wv3AAYCWwFToyUafkIoenMAhOap\n6PiewHzgBTM7lLC20wnuvld03FaVeLMiG6M7DZHGzXT3t6Jlsd8CCjUjJhFqZnwRuAvA3Z8BtjOz\n9oQ7iCfcfY27LyGsjtqR+peBH+fu86JrTAR2Bf4NfGJmw8zs68DK+N6iSGmUNEQat6poez2wOtp2\nQolTp+F6IKuLttfR8N39qjrHbebu6wiVKR8CjgGealrYIuWnpCHScs8R1SGImpred/flNJxIlgPt\nGzupmbUFOrj7SOC/gP3LEq1IC6hPQ6RxdZeC9jrbVwO3Rcufr+Cz2gv1Vnlz9yVRR/okQmnNJxu4\nRnvgb2a2BSEBXdzSNyLSUloaXURESqbmKRERKZmShoiIlExJQ0RESqakISIiJVPSEBGRkilpiIhI\nyZQ0RESkZEoaIiJSsv8PaXLtlrcgX/kAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "name": "stdout", "output_type": "stream", "text": [ "hydrograph is ploted in figure\n", "\n", "storage capacity of given plant is 3000sec-m-month\n" ] } ], "source": [ "from __future__ import division\n", "%matplotlib inline\n", "from matplotlib.pyplot import plot, title, xlabel, ylabel, show\n", "\n", "flow=[0, 1000, 800, 600 ,400 ,400 ,1200, 2400 ,2400, 1000 ,400 ,400 ,1000] #flow in matrix from in the order of months\n", "y=range(0,13)\n", "yy = [(0,1),(1,2),(2,3),(3,4),(4,5)]\n", "xx = [(0,1),(1,2),(2,3),(3,4),(4,5)]\n", "h=150\n", "e=0.85\n", "avg=sum(flow)/12\n", "print \"\\naverage rate of inflow is %dcu-m/sec\"%(avg)\n", "p=(735.5/75)*avg*h*e\n", "print \"\\npower developed is %0.2f MW\"%(p/1e3)\n", "plot(y,flow)\n", "title('hydrograph')\n", "xlabel('months')\n", "ylabel('run in cu-m/sec')\n", "show()\n", "print \"hydrograph is ploted in figure\"\n", "flow1 = range(0,12)\n", "for x in range(0,12):\n", " t=flow[x]\n", " a=avg\n", " if t" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "%matplotlib inline\n", "from matplotlib.pyplot import plot, title, xlabel, ylabel, show\n", "\n", "flow=[1500, 1000, 500, 500, 500, 1200, 2900, 2900, 1000, 400 ,600 ,1600]\n", "cod=1000#constant demand\n", "#plot2d2(flow)\n", "plot(flow)\n", "title('hydrograph for exp 8.3')\n", "xlabel('months')\n", "ylabel('run off in m**3/sec')\n", "avg=sum(flow)/12\n", "if codcod|t==avg:\n", " t=0\n", " else:\n", " t=cod-t\n", " #end \n", " flow1[x]=t \n", " #end\n", " \n", " else:\n", " flow = range(0,12)\n", " flow1 = range(0,12)\n", " for x in range(0,12):\n", " t=flow[x]\n", " a=cod\n", " if t>a|t==avg:\n", " t=0\n", " else:\n", " t=t-cod\n", " #end\n", " flow1[x]=t \n", " #end\n", "#end\n", "\n", "sto=sum(flow1)\n", "print \"storage capacity of plant is %dsec-m-month\"%(sto)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## example 8.4 Page 135" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "storage capacity of plant is -5149sec-m-month\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYQAAAEACAYAAACznAEdAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XmUFdW5/vHvEw0qhp/G6EURrmDSBDCOmGAcjxOiUXFl\nAExiuIojGofEKO1dSruMijdOaBxuHMEIihoVDSFwkXYGokFFEZFcUUBoTRyIESPD+/ujqvXcDjT0\n6dNdZ3g+a7Gos09V9VvQ67xn73fXLkUEZmZmX8g6ADMzKw1OCGZmBjghmJlZygnBzMwAJwQzM0s5\nIZiZGbCehCDpdkkNkuas5b2fS1ojaau8tlpJr0uaJ6l/XntfSXPS90bntW8i6d60fYakHYp1YWZm\n1jLr6yHcAQxo2iipG3Ao8GZeWx9gMNAnPeZGSUrfvgkYFhE1QI2kxnMOA/6Wtl8DXNGKazEzs1Zo\nNiFExJPA+2t562rgvCZtA4HxEbEyIhYCC4B+krYDOkXErHS/scAx6fbRwJh0+wHg4BZfgZmZFUWL\nawiSBgKLI+KlJm91ARbnvV4MbL+W9iVpO+nfiwAiYhXwYf4QlJmZtZ+NW7KzpI7ABSTDRZ81FzUi\nMzPLRIsSAvBVoDvwYloe6Ao8L6kfyTf/bnn7diXpGSxJt5u2k77378DbkjYGtoiI95r+UElecMnM\nrAARscFf2ls0ZBQRcyKic0T0iIgeJB/se0REAzARGCKpg6QeQA0wKyKWAcsl9UuLzMcBD6ennAgM\nTbe/D0xr5mdX7J+RI0dmHoOvz9dXbddWDdfXUuubdjoeeAboKWmRpOObfk7nfWDPBSYAc4E/AMPj\n84iGA7cCrwMLImJy2n4b8BVJrwNnAyNafAVmZlYUzQ4ZRcSx63l/xyavLwMuW8t+zwM7r6X9n8Cg\nDYrUzMzalO9ULgG5XC7rENqUr698VfK1QeVfX0upkHGm9iYpyiFOM7NSIoloq6KymZlVLicEMzMD\nnBDMzCzlhGBmZoATgpmZpZwQzMwMcEIwM7OUE4KZmQFOCGZmlnJCMDMzwAnBzMxSTghmZgY4IZiZ\nWcoJwczMACcEMzNLOSGYmRnghGBmZiknBDMzA5wQzMws5YRgZmaAE4KZmaWaTQiSbpfUIGlOXtuv\nJL0q6UVJv5O0Rd57tZJelzRPUv+89r6S5qTvjc5r30TSvWn7DEk7FPsCzcxsw6yvh3AHMKBJ2xRg\np4jYFZgP1AJI6gMMBvqkx9woSekxNwHDIqIGqJHUeM5hwN/S9muAK1p5PWZmVqBmE0JEPAm836Rt\nakSsSV/OBLqm2wOB8RGxMiIWAguAfpK2AzpFxKx0v7HAMen20cCYdPsB4OBWXIuZWcn64x/hppuy\njqJ5ra0hnABMSre7AIvz3lsMbL+W9iVpO+nfiwAiYhXwoaStWhmTmVlJGT8ehg6F3XbLOpLmbVzo\ngZL+E/g0IsYVMZ51qqur+2w7l8uRy+Xa48eambXKddfBlVfCtGmw005t+7Pq6+upr68v+HhFRPM7\nSN2BRyJi57y2/wBOAg6OiE/SthEAETEqfT0ZGAm8CUyPiN5p+7HA/hFxWrpPXUTMkLQxsDQitllL\nDLG+OM3MSkkEXHgh3HcfTJkCO2QwZUYSEaH175lo8ZBRWhD+BTCwMRmkJgJDJHWQ1AOoAWZFxDJg\nuaR+aZH5OODhvGOGptvfB6a1NB4zs1KzejWcckpSN3jqqWySQSGaHTKSNB44ANha0iKSb/y1QAdg\najqJ6NmIGB4RcyVNAOYCq4DheV/rhwN3ApsBkyJictp+G3CXpNeBvwFDinlxZmbt7ZNP4Ec/gg8/\nhMceg06dso5ow613yKgUeMjIzMrB8uUwcCBssw3cdRdsskm28bT5kJGZmf2rhgbI5aB372RWUdbJ\noBBOCGZmrfS//wv77gtHHw033AAbbZR1RIVxQjAza4WXXoL99oOzz4a6OtAGD9CUnoLvQzAzq3ZP\nPgnf+x5cfz0MHpx1NK3nhGBmVoCJE2HYMBg3Dg49NOtoisNDRmZmLXTHHXDyyTBpUuUkA3APwcys\nRf7rv+DGG6G+Hnr1yjqa4nJCMDPbAGvWwPnnJ72Cp56Crl3Xf0y5cUIwM1uPlSvhpJNg/vykkLxV\nha7J7IRgZtaMjz9OZhCtXg1Tp8Lmm2cdUdtxUdnMbB3efx/694ctt4SHH67sZABOCGZma/X227D/\n/vDNb8KYMfDFL2YdUdtzQjAza2L+fNhnn2TV0quvhi9UySelawhmZnmefx6OPBIuuQROPDHraNqX\nE4KZWeqxx2DIEPjNb+CYY7KOpv1VSUfIzKx599+fJIP77qvOZADuIZiZcfPNyRDRlCmw225ZR5Md\nJwQzq1oR8Mtfwp13whNPwFe/mnVE2XJCMLOqtHp18gyDJ55IlqLYbrusI8qeawhmVjUi4Nln4ayz\noFs3eO01ePxxJ4NG7iGYWUWLgBdegHvugXvvhc02g2OPhenT4etfzzq60uKEYGYVae7cJAncc08y\nPDRkSPJQm513Lu/HXLYlJwQzqxgLFiS9gHvvTdYhGjQI7r4b9tzTSWBDNFtDkHS7pAZJc/LatpI0\nVdJ8SVMkbZn3Xq2k1yXNk9Q/r72vpDnpe6Pz2jeRdG/aPkPSDsW+QDOrbIsWwVVXJWsO7bMPLF2a\nPMDmzTc/b3cy2DDrKyrfAQxo0jYCmBoRPYFp6Wsk9QEGA33SY26UPvtvuAkYFhE1QI2kxnMOA/6W\ntl8DXNHK6zGzKrBsGfz617DvvrD77jBvHowaBUuWfN5eLesPFVOz/2QR8STwfpPmo4Ex6fYYoPGe\nvoHA+IhYGRELgQVAP0nbAZ0iYla639i8Y/LP9QBwcIHXYWYV7r334NZb4ZBDoHdvmDkTamuTVUlv\nuQUOPhg29iB4qxTyz9c5IhrS7Qagc7rdBZiRt99iYHtgZbrdaEnaTvr3IoCIWCXpQ0lbRcR7BcRl\nZhVm+fLkOQT33JPcK3DYYTB8OBx+eDJbyIqrVfk0IkJSFCuY5tTV1X22ncvlyOVy7fFjzaydffwx\nPPpokgSmTYNcDn7846RQ/KUvZR1daauvr6e+vr7g4xXR/Oe5pO7AIxGxc/p6HpCLiGXpcND0iOgl\naQRARIxK95sMjATeTPfpnbYfC+wfEael+9RFxAxJGwNLI2KbtcQQ64vTzMpTRPL8gaefTh5R+Yc/\nQL9+yTTRY46BL3856wjLlyQiYoNL6oX0ECYCQ0kKwEOBh/Lax0m6mmQoqAaYlfYilkvqB8wCjgOu\na3KuGcD3SYrUZlbBVqyAP/0JnnkmSQLPPgudOsHee8OBB8J118E2//K10NpDsz0ESeOBA4CtSeoF\nFwEPAxOAfwcWAoMi4oN0/wuAE4BVwFkR8ce0vS9wJ7AZMCkizkzbNwHuAnYH/gYMSQvSTeNwD8Gs\nTC1d+vmH/9NPw8svwze+kSSAffZJ/u7SJesoK1NLewjrHTIqBU4IZuVh9erkA78xATzzDHz4YfKh\n35gA9twTOnbMOtLq4IRgZu3m73+HGTM+TwAzZyYLxTV+899nH+jZ0/cEZMUJwczaRERy92/jN/+n\nn06Withjj88TwLe/DVtvnXWk1sgJwcyKIiJ54PyTT36eACKSD//GBLD77tChQ9aR2ro4IZhZq730\nEvzsZ7BwIfTv/3kC6N7d6wKVk/aYdmpmFeqdd+Cii+DBB2HkSDj5ZC8HUU1c6jEzPv00WRl0p52S\nJSHmzUuWiHAyqC7+7zarYhHJQ2POPRd69UrWC/JTxKqXE4JZlWqsEyxbBjfckNQKrLp5yMisyrzz\nDpx6Khx6KHz3u8nzhp0MDJwQzKqG6wS2Pv5VMKtwrhPYhnJCMKtgrhNYS3jIyKwCuU5ghXBCMKsg\n//wnXHkl9OnjOoG1nH9NzCpA0zrB00+7TmAt54RgVuZeegnOOQcaGlwnsNbxkJFZmXrnHTjllKRO\n8L3vuU5greeEYFZm8usEHTu6TmDF418hszLyyCPJNFLXCawt+HkIZmXiqafgBz+AMWM8NGQbxg/I\nMatAq1ZB375wwQUweHDW0Vi5aGlCcA3BrAzccANssw0MGpR1JFbJ3EMwK3FLl8LOOydDRr16ZR2N\nlZN26yFIqpX0iqQ5ksZJ2kTSVpKmSpovaYqkLZvs/7qkeZL657X3Tc/xuqTRhcZjVql+8Qs48UQn\nA2t7BSUESd2Bk4A9ImJnYCNgCDACmBoRPYFp6Wsk9QEGA32AAcCN0meP6r4JGBYRNUCNpAEFX41Z\nhamvhyeegAsvzDoSqwaF9hCWAyuBjpI2BjoCbwNHA2PSfcYAx6TbA4HxEbEyIhYCC4B+krYDOkXE\nrHS/sXnHmFW1lSvh9NPhmmtg882zjsaqQUEJISLeA64C3iJJBB9ExFSgc0Q0pLs1AJ3T7S7A4rxT\nLAa2X0v7krTdrOpddx1065asVmrWHgq6MU3SV4Gzge7Ah8B9kn6cv09EhKSiVYLr6uo+287lcuRy\nuWKd2qzkLFkCl18Ozz4L2uCSoFW7+vp66uvrCz6+oFlGkgYDh0bEienr44C9gIOAAyNiWTocND0i\nekkaARARo9L9JwMjgTfTfXqn7ccCB0TEqU1+nmcZWVUZMgRqauCSS7KOxMpZe80ymgfsJWmztDh8\nCDAXeAQYmu4zFHgo3Z4IDJHUQVIPoAaYFRHLgOWS+qXnOS7vGLOqNG0azJwJtbVZR2LVpqAho4h4\nUdJY4DlgDfBn4DdAJ2CCpGHAQmBQuv9cSRNIksYqYHjeV/7hwJ3AZsCkiJhc8NWYlblPP00KyaNH\nJwvXmbUn35hmVkJGjUpuQHvkEdcOrPW8lpFZmXrrLdhjD5g1C3bcMetorBJ4LSOzMnXOOfDTnzoZ\nWHb8PASzEjB5cvLEs7vvzjoSq2buIZhl7JNPkp7B9dfDpptmHY1VMycEs4xdeSV84xtwxBFZR2LV\nzkVlswy98QZ885vw/POwww5ZR2OVxkVlszJy9tnJM5KdDKwUuKhslpFHH4V582DChKwjMUs4IZhl\nYMUKOPNMuPlm2GSTrKMxS3jIyCwDo0ZB377Qv//69zVrLy4qm7WzBQtgr72S+w66ds06GqtkLiqb\nlbCIZKjovPOcDKz0OCGYtaOHH4aFC5PZRWalxkNGZu3k44+hTx+44w448MCso7Fq4CEjsxJ16aWw\n995OBla63EMwawfz58M++8CLL0KXLllHY9XCPQSzEhMBZ5wBF1zgZGClzQnBrI098AAsXZokBbNS\n5iEjszb00UfQuzeMGwf77Zd1NFZt/AhNsxJy/vlJ72Ds2KwjsWrkhGBWIubOhQMOgJdfhs6ds47G\nqpGLymYloLGQfNFFTgZWPpwQzNrAPffA++/DaadlHYnZhis4IUjaUtL9kl6VNFdSP0lbSZoqab6k\nKZK2zNu/VtLrkuZJ6p/X3lfSnPS90a29ILOsLV8O554LN9wAG3uBeSsjrekhjAYmRURvYBdgHjAC\nmBoRPYFp6Wsk9QEGA32AAcCNkhrHtW4ChkVEDVAjaUArYjLLXF0dHHZYcleyWTkpqKgsaQtgdkTs\n2KR9HnBARDRI2haoj4hekmqBNRFxRbrfZKAOeBN4LE0qSBoC5CLi1CbndVHZysKcOXDwwfDKK7DN\nNllHY9WuvYrKPYB3Jd0h6c+SbpG0OdA5IhrSfRqAxnJaF2Bx3vGLge3X0r4kbTcrOxFw+ulw8cVO\nBlaeCh3h3BjYAzgjIv4k6VrS4aFGERGSiva1vq6u7rPtXC5HLpcr1qnNiuK3v01WND355KwjsWpV\nX19PfX19wccXOmS0LfBsRPRIX+8L1AI7AgdGxDJJ2wHT0yGjEQARMSrdfzIwkmTIaHrekNGxJENO\nHjKysvLBB8nS1g89BN/6VtbRmCXaZcgoIpYBiyT1TJsOAV4BHgGGpm1DgYfS7YnAEEkdJPUAaoBZ\n6XmWpzOUBByXd4xZ2bjoIjjqKCcDK2+tmRT3U+BuSR2AvwDHAxsBEyQNAxYCgwAiYq6kCcBcYBUw\nPO8r/3DgTmAzkllLk1sRk1m7mz0b7r03uTPZrJx56QqzVlizBvbdF044AU48MetozP4vL11h1o7G\njEmSwgknZB2JWeu5h2BWoPfeSwrJv/899O2bdTRm/8qrnZq1k+HDQUqWqDArRS1NCF5pxawAzz0H\nDz7oQrJVFtcQzFpozZqkd3D55fDlL2cdjVnxOCGYtdCtt0KHDvCTn2QdiVlxuYZg1gJ//SvstBNM\nmQK77pp1NGbNc1HZrA2ddBJsvjlce23WkZitn4vKZm1kxoxkiumrr2YdiVnbcA3BbAOsXp0sbf2r\nX8EWW2QdjVnbcEIw2wD//d/QqRP88IdZR2LWdlxDMFuPd96Bb3wDpk9PCspm5cJFZbMiO/54+MpX\n4Mors47ErGVcVDYroqefhqlTXUi26uAagtk6rFqV3JF81VVJ/cCs0jkhmK3DDTfA1lvDoEFZR2LW\nPlxDMFuLpUthl13gySehV6+sozErjIvKZkXw4x9D164walTWkZgVzkVls1Z6/HF44gkXkq36uIZg\nlmflyuSO5GuuSdYsMqsmTghmea67Lhkq+u53s47ErP25hmCWWrIkWdL62WehpibraMxar6U1BPcQ\nzFI//zmcdpqTgVWvViUESRtJmi3pkfT1VpKmSpovaYqkLfP2rZX0uqR5kvrntfeVNCd9b3Rr4jEr\n1LRpMHMm1NZmHYlZdlrbQzgLmAs0jueMAKZGRE9gWvoaSX2AwUAfYABwo6TGbsxNwLCIqAFqJA1o\nZUxmLfLpp0khefRo6Ngx62jMslNwQpDUFTgCuBVo/HA/GhiTbo8Bjkm3BwLjI2JlRCwEFgD9JG0H\ndIqIWel+Y/OOMWsXV18NX/saHHVU1pGYZas19yFcA/wC+H95bZ0joiHdbgA6p9tdgBl5+y0GtgdW\nptuNlqTtZu3irbeSVUxnzQJtcOnNrDIVlBAkHQm8ExGzJeXWtk9EhKSiTQ2qq6v7bDuXy5HLrfXH\nmrXIOefAT38KO+6YdSRmrVdfX099fX3Bxxc07VTSZcBxwCpgU5Jewu+AbwK5iFiWDgdNj4hekkYA\nRMSo9PjJwEjgzXSf3mn7scABEXFqk5/naadWdJMnJ7WDV16BTTfNOhqz4muXaacRcUFEdIuIHsAQ\n4LGIOA6YCAxNdxsKPJRuTwSGSOogqQdQA8yKiGXAckn90iLzcXnHmLWZTz5JegbXX+9kYNaoWGsZ\nNX59HwVMkDQMWAgMAoiIuZImkMxIWgUMz/vKPxy4E9gMmBQRk4sUk9k6XXll8ljMI47IOhKz0uE7\nla3qLFwIe+4Jzz8PO+yQdTRmbcd3Kputx1lnwc9+5mRg1pSXv7aq8uijybLWEyZkHYlZ6XFCsKqx\nYgWceSbcfDNssknW0ZiVHg8ZWdW44gro2xf691//vmbVyEVlqwp/+Qv06wezZ0O3bllHY9Y+XFQ2\nayIiuefgvPOcDMya44RgFe/hh5OppmefnXUkZqXNQ0ZW0T7+GPr0gdtvh4MOyjoas/blISOzPJde\nCt/+tpOB2YZwD8Eq1vz5sPfe8NJL0KVL1tGYtT/3EMxICslnnAEXXOBkYLahnBCsIj3wACxdmswu\nMrMN4yEjqzgffQS9e8Pdd8P++2cdjVl2Wjpk5IRgFef885PewdixWUdili0nBKtqc+fCAQfAnDmw\n7bZZR2OWLReVrWotWgRDh8KFFzoZmBXCCcHKXgTcdRfssQcMHJg8J9nMWs7LX1tZe/ddOOWU5J6D\nKVNg992zjsisfLmHYGXr4Ydhl13ga1+D555zMjBrLfcQrOx8+GGyUN0TTyRPPttvv6wjMqsM7iFY\nWZk+HXbdFTp0gBdecDIwKyb3EKwsrFgBtbVw//1wyy1w+OFZR2RWedxDsJI3a1ZSH2hoSBaqczIw\naxsFJQRJ3SRNl/SKpJclnZm2byVpqqT5kqZI2jLvmFpJr0uaJ6l/XntfSXPS90a3/pKsUqxcCRdd\nBEcdBRdfDOPHw1ZbZR2VWeUqtIewEjgnInYC9gJOl9QbGAFMjYiewLT0NZL6AIOBPsAA4EZJjXfP\n3QQMi4gaoEbSgIKvxirGK6/AXnsls4dmz4bBg7OOyKzyFZQQImJZRLyQbn8EvApsDxwNjEl3GwMc\nk24PBMZHxMqIWAgsAPpJ2g7oFBGz0v3G5h1jVWj1arjqKsjl4NRT4fe/9/LVZu2l1UVlSd2B3YGZ\nQOeIaEjfagA6p9tdgBl5hy0mSSAr0+1GS9J2q0JvvJEsPQEwcybsuGO28ZhVm1YlBElfAh4AzoqI\nv38+CgQREZKKtiJdXV3dZ9u5XI5cLlesU1vGIuC225JZROefD+ecAxttlHVUZuWnvr6e+vr6go8v\neLVTSV8EHgX+EBHXpm3zgFxELEuHg6ZHRC9JIwAiYlS632RgJPBmuk/vtP1Y4ICIOLXJz/JqpxVq\n6VI46SR4++1kPaKddso6IrPK0S6rnaYF4duAuY3JIDURSDv9DAUeymsfIqmDpB5ADTArIpYByyX1\nS895XN4xVuEmTIDddksWpZsxw8nALGsF9RAk7Qs8AbwENJ6gFpgFTAD+HVgIDIqID9JjLgBOAFaR\nDDH9MW3vC9wJbAZMiogz1/Lz3EOoIO+9lzzv+M9/Th5i861vZR2RWWXyA3KspE2eDCeeCN//Plx2\nGXTsmHVEZpWrpQnBS1dYu/joIzj33CQhjB0LBx2UdURm1pSXrrA29/TTSa3gn/+EF190MjArVWXT\nQ1izBr7g9FU2VqyA+np48EF45BG4+ebkaWZmVrrK5iP2hz9MPmSsdC1YANdfD0ccAZ07JzWC7t2T\nXoGTgVnpK5seAsDBBydPydpmm6wjMfi8F/CHPyR/PvooWYn0+ONh3DjYcsv1nsLMSkjZzDJavTq4\n6KJkxcvf/x569co6quq0YMHnCeCpp5KH1Rx+ePJnt91AGzyfwczaWsVPO73zzmR5g/HjXZxsD+vq\nBRx+OBx6qHsBZqWs4hMCJB9QgwfDqFHJ8IQV17p6AUcckWy7F2BWHqoiIQC89hp85zswaBD88pee\ngdQaTXsB//gHDBjgXoBZuauahADw7rtwzDHQtWsylLTZZu0fW7laWy/giCOSJOBegFllqKqEAPDJ\nJ8mw0cKFyQykf/u39o2tXLgXYFZ9qi4hQHLT2siRcPfdyQyk3r3bMbgS5l6AWXWryoTQaMwYOO+8\n6p2B5F6AmeWr6oQAn89AuvxyOOGEto2rFLgXYGbrUvUJASp7BpJ7AWa2oZwQUn/9azIDqUuXZCip\nnGcguRdgZoVwQsjzySfJsNEbb5TXDKQVK+Dxx5MEMGmSewFmVhgnhCYikhlIv/0tPPoo9OlT5OCK\nxL0AMys2J4R1GDs2eWLXuHFwyCFFCqwV3Asws7bmhNCMxx9PCs2XXQbDhhUhsBZyL8DM2pMTwnq8\n9hoceWTykPdLL23bGUjuBZhZlpwQNkBbzkByL8DMSkVZJgRJA4BrgY2AWyPiiibvFzUhQDIDadgw\n+MtfkhlInTsXdh73AsysVLU0IWR+y5akjYBfAwOAPsCxktp8NaJNN01mHh12GOy1F8ydu+HHNn12\n8KWXwrbbwn33wZIlcPvt8IMfbHgyqK+vL+gayoWvr3xV8rVB5V9fS2WeEIBvAQsiYmFErATuAdrl\nkewSXHxx8ieXg//5n7Xvt2IFTJ4MZ50FPXvC/vvD7NnJKqtvvQVPPgm1tYU/QrLSfyl9feWrkq8N\nKv/6WmrjrAMAtgcW5b1eDPRrzwB+8hPYYYdkBtKll8KJJ667FjBhgmsBZlaZSiEhZF/EAA44IPmm\n/53vwCWXwKefJnWA449P7l1wLcDMKl3mRWVJewF1ETEgfV0LrMkvLEsqiaRhZlZuymqWkaSNgdeA\ng4G3gVnAsRHxaqaBmZlVmcyHjCJilaQzgD+STDu9zcnAzKz9Zd5DMDOz0lAK006bJWmApHmSXpd0\nftbxFJOkbpKmS3pF0suSzsw6pmKTtJGk2ZIeyTqWYpO0paT7Jb0qaW5aD6sYkmrT3805ksZJ2iTr\nmFpD0u2SGiTNyWvbStJUSfMlTZFUttNH1nF9v0p/P1+U9DtJWzR3jpJOCFndtNaOVgLnRMROwF7A\n6RV2fQBnAXMpkdlkRTYamBQRvYFdgIoZ6pTUHTgJ2CMidiYZzh2SZUxFcAfJZ0m+EcDUiOgJTEtf\nl6u1Xd8UYKeI2BWYD9Q2d4KSTghkeNNae4iIZRHxQrr9EckHSpdsoyoeSV2BI4BbgYq6cyP9prVf\nRNwOSS0sIj7MOKxiWk7yhaVjOvGjI7Ak25BaJyKeBN5v0nw0MCbdHgMc065BFdHari8ipkbEmvTl\nTKBrc+co9YSwtpvWts8oljaVfiPbneQ/rVJcA/wCWLO+HctQD+BdSXdI+rOkWyR1zDqoYomI94Cr\ngLdIZv99EBHruJe/rHWOiIZ0uwEocFWzsnACMKm5HUo9IVTiMMO/kPQl4H7grLSnUPYkHQm8ExGz\nqbDeQWpjYA/gxojYA/gH5T3c8H9I+ipwNtCdpNf6JUk/yjSoNpauoFmRnzmS/hP4NCLGNbdfqSeE\nJUC3vNfdSHoJFUPSF4EHgN9GxENZx1NEewNHS3oDGA8cJGlsxjEV02JgcUT8KX19P0mCqBR7As9E\nxN8iYhXwO5L/00rTIGlbAEnbAe9kHE/RSfoPkqHb9Sb0Uk8IzwE1krpL6gAMBiZmHFPRSBJwGzA3\nIq7NOp5iiogLIqJbRPQgKUY+FhE/yTquYomIZcAiST3TpkOAVzIMqdjmAXtJ2iz9PT2EZHJApZkI\nDE23hwKV9KWs8dECvwAGRsQn69u/pBNC+s2k8aa1ucC9FXbT2j7Aj4ED06mZs9P/wEpUiV3xnwJ3\nS3qRZJbRZRnHUzQR8SIwluRL2Utp82+yi6j1JI0HngG+LmmRpOOBUcChkuYDB6Wvy9Jaru8E4Hrg\nS8DU9PPlxmbP4RvTzMwMSryHYGZm7ccJwczMACcEMzNLOSGYmRnghGBmZiknBDMzA5wQzMws5YRg\nZmYA/H9u+wc6AAAABElEQVS6oEf5Os2H4wAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "%matplotlib inline\n", "from matplotlib.pyplot import plot, title, xlabel, ylabel, show\n", "\n", "flow=[1500, 1000, 500, 500, 500, 1200, 2900, 2900, 1000, 400 ,600 ,1600]\n", "cod=1000#constant demand\n", "#[m n]=size(flow)len\n", "n=len(flow)\n", "mf = range(0,n)\n", "mf[0]=1500\n", "for i in range(2,n):\n", " mf[i]=mf[i-1]+flow[i]\n", "\n", "plot(mf)\n", "from numpy import arange\n", "dd= arange(1,mf[n-1],cod)\n", "\n", "avg=sum(flow)/12\n", "if codcod|t==avg:\n", " t=0\n", " else:\n", " t=cod-t\n", " #end\n", " flow1[x]=t \n", " #end\n", " \n", "else: \n", " for x in range(0,12):\n", " t=flow[x]\n", " a=cod\n", " if t>a|t==avg:\n", " t=0\n", " else:\n", " t=t-cod\n", " #end\n", " flow1[x]=t \n", " #end\n", "#end\n", "\n", "sto=sum(flow1)\n", "print \"storage capacity of plant is %dsec-m-month\"%(sto)\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## example 8.5 Page 154" ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "load duration data is as under\n", "[0, 20, 40, 50, 80, [], 100, [], 120, 150, 200, 250, 0, 1, 2, 3, 5, [], 7, [], 8, 9, 10, 11]\n", "(a)\n", "meanflow is 459.104938m**3-sec\n", "(b)\n", "power delevered in 36018312kW=36018.313MW\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAEZCAYAAAB8culNAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xe8U/X9x/HXu4gLRXFUHCio4B7gouK4OCqCVftrXbVa\nRW2r1lU7gA6w9mEVa6uitraVKgioqHUBAo6rVqsoIiqIgICCAiKCOJB1P78/vifemLtyk5ycJPfz\nfDzyIDk5OflEMJ9850dmhnPOOZfuG0kH4JxzrvR4cnDOOVeHJwfnnHN1eHJwzjlXhycH55xzdXhy\ncM45V4cnB1f2JM2TdHSpXKcYJFVLOi/pOFzlWq+pEyQdBgwEOqadb2a2c4xxOdccFt1K5TrFUE6x\nujLUZHIA7gAuB14F1sUbjnOlSVIrMyvIv39J65nZ2kJcy7m4ZNOttNzMxpnZYjP7KHWLPTLnmqer\npKmSlku6R9IGkt6UdELqBEmtJX0kab/o8VmS3o2ODUi/mKRBku6XNFzSJ8CPJG0n6RFJSyXNknR+\n2vkbSbpL0seSpkv6laT5ac/Pi469DnwqqZWkfpJmS1ohaZqkk9POP0fS85KGRJ/pLUlHZXzmjpL+\nG71+vKQtC/zf1LVgDSYHSQdIOgB4WtL1kr4lqVvqVsQYnWuKgFOA44BOwL7AOcBdwA/TzusNvG9m\nUyXtCdwGnAlsB2wJ7JBx3ROB0Wa2GTASuAd4D9gW+D5wjaSe0bkDgR2j9z82et/Mbp/TgeOBzaNW\nyGzgMDNrC1wF3C1pm7TzD47O2TK6/oOSNk/7zD+IPuc3gfWBXzT9n8q57DTWcrgB+DNwCHAgcE10\nLHVzrlQYcLOZLTKzZcCjwP7A3UAfSZtE550FDI/ufx941Mz+a2argd8BNRnXfcHMHonubw0cCvza\nzFab2VTgX8DZ0fOnANeY2Sdm9j5wE+ELPDPG981sFYCZ3W9mi6L79wGzCP+/pXxoZjeZ2bro+beB\nE9KuN9TMZpvZl8B90Wd2riAaHHMwsyoASTub2Zz05yT5YLQrNYvS7q8EtjWzhZKeB74v6SGgF3BJ\ndM62wILUC8zsC0lLM665IO3+dsDHZvZ52rH3gAPSnp+f9lz6a1PSn0fS2cAVhMkeAJsQWgkp72e8\n/t0o7pTMz7wJzhVINmMO99dzbHShA3GuwFK/2lNdS6cQWgILo+MLgQ5fnSxtzNe/mOHr3UIfAFuk\ntUIgdCOlvsC/dr2M+3WuJ2kn4B/AxcAWZtYOeJOvtza2z3j9TlEczsWuwZaDpD2APYHNJf0f4R+t\nAW2BDYsTnnN5ewi4FdgGuC7t+P3AS5J6AC8Df6CRH0tmNl/SC8CfJP0C2A3oS+j3h9Ct01/Sy0Ab\n4Gc0PtW0TfT8R8A3olbE3hnnfFPSpcDfgJOB3YGxac8L52LSWMuhC/AdYLPozxOiP7sBF8QfmnM5\n+2oNgJmtBB4kdN08+NUJZtMJv9pHEn6Nf8zXu33qW0dwRnSdD6Jr/d7Mnoqe+wOhK2kuMIHQul7d\nYIDh/W8A/kfoHtob+G/GaS8BnYElwNXA96IxlfQYG4vXuZypqWI/kg41sxeafWGpAzCMMJPCgH+Y\n2c2SBgHnE/7BAwwws3HRa/oTfo2tAy41swnNfV/nMkn6HdDZzM5u8uTCveeFwKlm1rPJk+t//TnA\neWZ2eEEDcy5L2SyCmy3pN9RdId23idetAa4ws9eiftrJkiYSEsVfzOwv6SdHUwtPI3RlbQ88IamL\nmWXOIHEua5K2IPzgOCvm92kP7EJoCXQGfg4MifM9nYtTNsnhYeBZYCK1U/2abL5GU/RS0/Q+k/QW\ntQNs9fWVngSMMrM1wDxJswnzvF/MIkbn6pB0AfBXYJiZZXbZFNr6wN8J6xyWA6MI6yhy5d1ELlHZ\ndCu9ZmZ5zZ+W1BF4BtgLuBI4F/gEeAW40syWSxoCvGhmI6LX/AsYZ2YP5PPezjnnmi+bqayPSeqT\n6xtEXUr3A5eZ2WeEmRedCAt2FtL4gjr/5eSccwnIplvpcmCApNWEcQQIYw5tm3qhpNbAA8DdZvZQ\n9MIP057/F2E1K4T54ulzw3eg7iIgJHnCcM65HJhZ1tOfm2w5mNkmZvYNM9vQzDaNbtkkBhF2dJ1u\nZjemHU9f4fld4I3o/iPA6ZLWl9SJMKg3qYGYKvY2cODAxGPwz+afb906Y/BgY+utjdGjw7FDDx3I\nT36SfGz+95fbrbmyaTkg6STgCEI3zzNm9mgTLwHoQViZ+rqkKdGxAcAZkvaPrjUX+AmEed+S7gOm\nA2uBiyyXT+Scy8uiRXD22fDFF/Dyy7DTTuF4jx4wdChceSV07pxsjC5+2RT7uRY4CBhBmGV0abT2\noX9jr7MwO6S+lsm4Rl5zDWGDP+dcAh5/HPr2hfPPh9//HtZL+4bYeGO44opwfNSo5GJ0xZFNy6EP\nsL9FhU4k3Qm8BjSaHFxuqqqqkg4hNpX82aC8P9/q1TBgANx7L4wcCfV9lKqqKg48MLQaXnsN9q+w\nPWDL+e8vDtlMZX0d6GlmS6PHWwJPm9m+RYivvni8t8m5Apo1C844A3bYAe64A7ZsomTQkCEwbhyM\nHdv4ea60SMIKOSAN/Al4NapydRcwGe/6ca4iDBsGhx4K554L//lP04kB4Mc/hrfegueeiz8+l5wm\nWw4AkrYjjDsYMMmiAiVJ8JaDc/lbsQIuvhgmT4Z77oF9m9kPMGwY/OMfIUHI94YtC3G0HCBUwQJo\nDRwabeHtnCtDL78M3bqFAeZXXml+YgA480xYvhzGjCl8fK40ZDPm8G9gH2AaaWUUzezceENrMB5v\nOTiXg5oauOEGuP56uO02+P7387veww/D734XBqe/ke3PTJeY5rYcskkO04G9SuUb2ZODc82XvnZh\nxIjatQv5MAvjFT/7WWhJuNIWR7fSy4RttJ1zZejxx0M3UvfuUF1dmMQAYazhT38K6x5WN1jWyJWr\nbFoOVYStLRYBq6LD5lNZnStt6WsXhg+vf+1CIRx3HJx8Mlx4YTzXd4URR7fSO8AVhOLn6WMO83KM\nMS+eHJxrWnPXLuRj8mT4znfCe7ZpE9/7uPzE0a30oZk9YmZzzGxe6pZ7iM65OOWydiEfBxwAhx0W\nFse5ypFNy+E2YHPC1tqpnkUzswcbflWjNaS3AO4FdgLmEersLo9e02QNaW85OFe/fNcu5OPtt0OC\nmDkT2rUr3vu67MXRctiYkBS+DZwQ3b6TxetSNaT3AroDF0vaA+gHTDSzLsCT0ePMGtK9gNsk+QQ5\n57JQiLUL+dhttzDuMHhwcd/XxafBloOkHwDjU3sq5f1G0kPALdHtSDNbHBVlrzaz3aNWQ42ZXRed\n/zgwyMxezLiOtxycixR67UI+FiyA/faDN9+Ebbdt+nxXXM1tOTS2K+uOwGhJ6wNPELbanpTLN3NU\nQ7or8BKwjZktjp5aDGwT3d8OSE8EC4Dtm/tezrUUDdVdSMoOO4RxjquvDonKlbcGu23M7FozOwro\nDbxOGAt4VdIoSWdL2qah16aLakg/QKgh/WnGexiN14n2JoJz9Yhr7UK++vWD++6Dd95JOhKXrybr\nOZjZCuDB6IakvYDjgeGEcYgGpdWQHm5RDWlgsaT2ZrYoKhmaqimdVQ1pgEGDBn11v6qqyvdhdy1G\nNnUXkrTVVnDppWFh3IgRSUfTslVXV1NdXZ3z6xudrSSpVVqRn7aEus6zooTR+IVDDem7gKVmdkXa\n8cHRsesk9QM2N7N+0YD0SOBgQnfSE8Cumd1YPubgWqpirl3Ix6efhoJAEyYUf2DcNaxgs5UknUb4\nlf9OVEP6dWAw8KakXllcO1VDuqekKdGtF3AtcKykmcBR0WPMbDqQqiE9Dq8h7dxXir12IR+bbgr9\n+8NvfpN0JC4fjc1Wep3QbbQRYUfWbmY2Q9JOwGgzO7h4YX4tLs8ZrsVIcu1CPlatCtNbR4yAHj2S\njsZBYdc5rDOzRWY2F5hrZjMAzOxdQl0H51yMkl67kI8NNoBBg8IAtf+WK0+NLjJLW4TWN+3Yenhy\ncC42NTVh3UKfPnDttXD77SFBlJuzzoKlS0O9aVd+GutWOhh4w8xWZhzvCBxmZnfHHl39cXm3kqtY\ncdRdSNJ//gNXXQWvvuoFgZJWsG4lM5sEdIkuum/a8XlJJQbnKlmprl3Ix8knhy6m++5LOhLXXE3l\n8r6SOgPnFSMY51qi1avhF7+ACy4Iaxf+8AdYr8kVSOVBgmuugd/+FtasSToa1xyNTWUdGD3/Unio\ngUWLyrkWYtasMEV19uxQi7nUFrUVwtFHQ6dOMHRo0pG45mhqEdyJhNXQ48zskaJF1Qgfc3CVYtgw\nuPLKMKvnoovCr+xK9fLLoYtp1qzyHFyvBIXesvsQM7sQOCi/sJxzKStWhJk8114LTz4Z1jFUcmIA\nOOgg+Na34JZbko7EZaup5JAaRhoddyDOtQTlvHYhX1dfDX/+MyxfnnQkLhs+IO1cEVTK2oV87LEH\nnHBC+O/gSl9j6xwGAlsBZwJ3EzbLu6qIsdXLxxxcuam0tQv5eO896NoVpk2D9u2TjqZlKeQ6h6uA\niYR6z0+UQmJwrtyMGxe+DCtp7UI+dtwxJMo//jHpSFxTYh2QljRU0mJJb6QdGyRpQdpOrcenPddf\n0ixJMyQ1WivCuVK2alWYifTjH8OoUZW1diFfAwaE/yZz5yYdiWtMo1NZ8764dDjwGTDMzPaJjg0E\nPjWzv2Scm6rncBC19Ry6mFlNxnnereRK2syZoe5Chw6lXXchSYMGwZw5YTqvK45CT2VF0neiX/jL\nJH0a3Zos9gNgZs8By+q7bD3HTgJGmdkaM5sHzCYU/nGuLJiFL7sePaBv39Kvu5Ckn/8cxo+HN95o\n+lyXjGy2wroR+BGwpZltGt3a5vm+l0iaKukOSZtHx7YDFqSds4DQgnCu5LXEtQv5aNs2bOf9298m\nHYlrSDa9oAuAaZndO3n4G/CH6P7VwA00PFW23v4jryHtSsmkSaEb6eijw9qFljZFNVcXXgh//Su8\n8ELYQsQVVqw1pAEkdSd8mT8NrI4OW+aYQSOv7wg8mhpzaOi5qJ40ZnZt9NzjwEAzeynjNT7m4EpC\nTU1Y1PXnP8Ntt8H3v590ROVn6NDQFff0097SilvBxxwIv+4/AzYENolum+YWHkjaNu3hd4FUr+Mj\nwOmS1pfUCegMTMr1fZyL06JF0KsXPPJIWPXsiSE3Z58NixfDhAlJR+IyZdOttK2ZHZvLxSWNAo4E\ntpI0HxgIVEnan9BlNBf4CYCZTZd0HzAdWAtc5E0EV4rGjQsDzhdcAL//vU9Rzcd664VtNfr3h2OP\n9YJApSSbbqXBwJNmNr44ITXOu5VcUlatCnP077sPhg+vzO21k2AWNub71a/g1FOTjqZyNbdbKZvk\n8BmwMWG8IVWuwwowYyknnhxcEnztQrwmTgwzvKZNg9ZeoT4WBR9zMLNNzOwbZrZhAaeyOlcWfO1C\ncRxzTEi8d96ZdCQuJdYV0nHwloMrlhUrQhGeV1+Fe+5pWdtrJ+Gll8LA/syZsNFGSUdTeeKYreRc\nizNpUtgwr02blld3ISmHHAIHHgi33pp0JA685eDc1/jahWRNmwY9e4ZyopttlnQ0laXgA9LRRVsB\n25A29dXM3sspwjx5cnBx8boLpeGcc8L4w9VXJx1JZYljttIlhPUJHwLrUsfrW/FcDJ4cXBx87ULp\nmDcPDjgA3nwTtt22ydNdlpqbHLL5X+ByYDczW5p7WM6VpvS1C6NG+dqFUtCxI1x6KXTpAocfHkqr\n9u4NnTolHVnLkk3L4Wng22a2ptETi8RbDq5QfO1CaVu+PGyrMXZsaNltsUVtojjsMFh//aQjLC9x\ndCsNBboAY8hh471C8+Tg8mUWVjhfeWUoOnPRRb7pW6mrqYHJk0OiGDMmJPajjw7J4vjjvfspG3Ek\nh0HR3dSJIiSHRGpKe3Jw+fC1C5Vh8WJ4/PGQKCZOhJ13rm1VHHQQtGqVdISlJ5bZStGFNwUws0+b\nEcxQoA/wYVqZ0C2Ae4GdgHnAqWa2PHquP9CXMPB9qZnV2avRk4PLVaruwjHHhDoCXnehMqxZE2pC\npFoVixeHHXN794bjjgvdUS6elsM+wDAg1SO7BPiRmb2ZRTD11ZAeDHxkZoMl/RpoZ2b9vIa0i4uv\nXWhZ3n03jFGMGQPPPBNah6lWxb77ttwuxDiSw/+AAWb2dPS4CrjGzLKq3ZRZ7EfSDOBIM1ssqT1Q\nbWa7R62GGjO7LjrvcWCQmb2YcT1PDi5rvnahZfvyS6iurm1VrFoVkkSfPmHMYpNNko6weOLYPmPj\nVGIAMLNqoE0OsaVsY2aLo/uLCYvrwGtIuwIbNy5sgdG9e/iC8MTQ8my4YehiuvlmmD071PfefXcY\nMiQMYn/723DjjWFFtvu6bNY5zJX0O2A4YTD6TGBOId7czExSY80AryHtcvLoo/DTn/raBVdLgt12\nC7ef/zxMTnjyydCiGDw47KOV6n468kjYYIOkI85PMWpIbwFcBfSIDj1H6O5ZltUb1N+tVGVmi6KS\noU9H3UpeQ9oVxJIlsN9+YWHbYYclHY0rB2bw2mu13U/TpoUfFampsh06JB1h/mKbrZSrepLDYGCp\nmV0XJYTNMwakD6Z2QHrXzEzgycE1xgy+9z3o3Bmuuy7paFy5+ugjGD8+JIvx42H77WvHKrp3L8/t\nVQqWHCQ92sjrzMxOzCKYr2pIE8YXfg88DNwH7EjdqawDCFNZ1wKX1Vea1JODa8zw4aGL4JVXyr9b\nwJWGdetCrYlUq+K998JYRe/eYTxj662TjjA7hUwOVY29MBqYLjpPDq4h8+eHDdsmTID99086Glep\n3n8/THYYOzaMWeyxR+1YRdeu8I0SrZJTct1KhebJwdWnpib8mjv6aOjfP+loXEuxejU891xtq+KT\nT8IYRZ8+YbFlKdWkKGTL4Y1GXmdmlsjGA54cXH2GDAkzk559tjz7g11leOedkCjGjoXnnw+V7VJj\nFbvvnuwCvEImh46NvdDM5jUnsELx5OAyzZgRtnZ+4YUwEO1cKfjiC3jqqdpWxTe+UZsoevYsfp3s\nFtGt9Pzzxre+1XKXwbtaa9fCoYfCuefChRcmHY1z9TOD6dNDkhg7Nmz8WF0N3boVL4ZCthyeN7Me\nkj6j7mI0M7O2ecSZM0m2yy7GFlvAZZfBKaf4vu4t2dVXh+b7uHH+Y8GVj+XLw6K71q2L954touWw\ndq0xdmxY9j5jRtiC+cc/Lp8pZa4wJk8Og39TpoR56M65hsWxtxKS2knaV1K31C33EPPXqhV85zth\nGtm4cTB3bigpeP758EZjw+iuYqxcCWedBTfd5InBuThks33G1cA5hP2Uvto+28x6xhpZw/HUOyC9\nZAncfnvYknnPPUOXU58+pTvn2OXnyithwYJQsMe7k5xrWhxbds8E9jaz1Y2eWCRNzVZavRpGjw5d\nTsuXh0Ll55wDm25avBhdvKqr4cwz4fXXve6zc9mKo1tpGtAu95CKa/31wxfHpElw551hgUrHjmEX\nxrlzk47O5WvFijAz6Z//9MTgXJyyaTkcRNgP6U1gVXQ4q72V4pDLOod334Vbb4WhQ8N8+MsvhyOO\n8O6IcnTeeWGR2+23Jx2Jc+Uljm6lt4C/EZJDaszBzOyZnKMM150HrCDUi15jZgc3Vl867XU5L4L7\n/HMYNiwMYm60UUgSp5/uG7SVi0cegSuugKlTW1YFL+cKIY7k8LKZHZR3ZHWvOxc4wMw+TjtWb33p\njNflvUK6piZsw3vTTWEP95/+NNzat8/rsi5GXqPBufzEkRz+QuhOeoTabiXM7NVcg4yuOxc40MyW\nph2rt750xusKun3G9OmhhOC998KJJ4ZZTsVcteia5jUanMtfHMmhmnrKdeY7lVXSHOATQrfS7Wb2\nT0nLzKxd9LyAj1OP014Xy95KH38cBjlvuQV23jkkiZNOCmsqXLK8RoNz+SubFdKStjWzhZK2BiYC\nlwCPpCcDSR+b2RYZr4t14701a+DBB8NU2IULw8B1587h1qVL+NOnxRaP12hwrjCamxyavbmxpJOB\nhZm1nZvLzBZGfy6R9B9CedDFktqn1Zf+sL7XDho06Kv7VVVVVBWwgnzr1nDaaeH22mtha4aZM+H+\n+2HWLJg9OySH9GSRuu26K2y8ccFCafFqasK01Suu8MTgXHNVV1dTXV2d8+ub3XKQ9Cdgb6C1mfXK\n6U2ljYFWZvappDbABOAq4BjqqS+d8dpEt+yuqYEPPgiJYtaskDhS9+fOha22qj9x7Lyzd4k0l9do\ncK5wCrkr6ylmNlrSzmY2p2ARhmt3Av4TPVwPGGFmf4qmstZbXzrttSVbz2HdutANkp4wUvfnz4ft\ntqvbRdW5c1ik519+X+c1GpwrrEImhylm1jX1Z8EizFMpJ4fGrFkD8+Z9PXGkbgsXwk47hYRRVRX2\nhNptt5a7SM9rNDhXeIVMDk8QZikdBDyX8XRZrZAudV9+CXPmhGm1EyeGYiCtW9cWLa+qKn7VqCR5\njQbnCq+QyWF9oBswHDgfSL9o3iukc1WJySGTWdh6PFVecOrUMGsqlSx22inpCOPjNRqci0chk8Nw\nMztL0q/MbHDBIsxTS0gOmZYtCyu6x44Nv6a/+c3aRNGjR3GrScVp5cowbfV3v4Mzzkg6GucqSyGT\nw3TC7KHHgarM59O3vSimlpgc0q1bFxaDpVoV77wDxxwTkkWvXuW9BYjXaHAuPoVMDpcCFwI7Ax9k\nPG1mtnPOUeahpSeHTIsWhdbE2LFhvKJz59Ci6NMHDjywfIodeY0G5+IVx/YZfzezn+YdWYF4cmjY\nmjVhIHfMmJAsliwJ/fe9e8O3vw3tSrQqx4oVYVO9W28NsTrnCi+W7TMk7QccQZi99JyZTc09xPx4\ncsje3LmhVTFmTFhI1rVr7VjF3nuXTteN12hwLn5xtBwuAy4AHiTMWDoZ+KeZ3ZxPoLny5JCblSvh\n6adrxyrWravtfjrqKGjTJpm4vEaDc8URR3J4A+huZp9Hj9sAL5rZPnlFmiNPDvkzCyuQU91PL78c\nZj2lksUuuxQnDq/R4FzxxJUcDjazldHjjYBJnhwqx4oVtYvvxo6Ftm1ru58OPzyePaG8RoNzxRVH\ncvg5cA5f71a608z+mkecOfPkEK+amrAbbapVMX166Hbq3TvcCrUwzWs0OFdccQ1IHwAcRu2A9JTc\nQ8yPJ4fi+ugjePzxkCjGj4cOHWpbFd2751YMyWs0OFd8zU0OWc2CN7PJZnaTmd0cd2KQ1EvSDEmz\nojrSLUo++6/HYaut4Ic/hJEjYfHiMN3UDC6+OKzU/sEPYMSIkESaUl1dXdE1Gkrt767Q/PO1LCW1\nREpSK+AWoBewJ3CGpD2Sjaq4Svkf6HrrhYHra64JXU9Tp0LPnjB6dBjEPvRQ+OMf4dVXQwLJVF1d\nza23whdfwC9/Wfz441bKf3eF4J+vZSm1KgIHA7PNbB6ApHuAk4C3kgzK1W+HHeCCC8Jt1aqwlmLs\nWDj9dPjss7AAr0+fsL1H27ahdXHrraFGg9evcK60lVTLAdgemJ/2eEF0zJW4DTaAY4+Fv/411Kx4\n5hnYd9+wsG377eHoo8OU1auu8uI9zpWDbGYrfQ+4FtiG2m27zczaFjyY8F69zOyC6PEPgUPM7JK0\nc3w02jnnctCcAelsGveDgRPMrBhdO+8DHdIedyC0Hr7SnA/nnHMuN9l0Ky0qUmIAeAXoLKljVGzo\nNOCRIr23c865SDYth1ck3Qs8BKyOjpmZPVjoYMxsraSfAeOBVsAdRUxMzjnnItmMOdwZ3f3aiWZ2\nbkwxOeecS1hWK6RLgaRewI2EFsW/zKxiduSR1AEYBnyTkIT/kdSut3GK1rG8Aiwws+8kHU8hSdoc\n+BewF+HvsK+ZvZhsVIUjqT/wQ6AGeAM418xWJRtVbiQNBfoAH6b2iJO0BXAvsBMwDzjVzJYnFmQe\nGvh81wMnEHp/3iH8/X3S2HWaHHOQ9O+M29DozYumBSyOWwNcYWZ7Ad2Biyvs86VcBkwnoxXaXJJ2\nk/SapBWSLpF0p6SrCxNi1jHUSEqvhngTMNbM9gD2pcBrcySNlXRWIa/ZjPfuSNi2v1v0ZdMKOD2J\nWArk34TvknT9gIlm1gV4Mnpcrur7fBOAvcxsP2Am0L+pi2QzID0GeCy6PQlsBnzerFDz99XiODNb\nA6QWx1UEM1tkZq9F9z8jfLFsl2xUhSVpB6A34dd1vjPOfgU8aWZtzWwIIdkk1gSWtBlwuJkNhTB2\n1tSvsiauN0jS8PRjZtbbzIY39JqYrSD8gNlY0nrAxoSZhWXJzJ4DlmUcPhG4K7p/F2GD0bJU3+cz\ns4lmVhM9fAnYoanrNDkgbWb3pz+WNBJ4PvtQC6K+xXGHFDmGooh+pXUl/AVWkr8CvwQKsT5mJ+CF\njGNJTnHuBCyR9G9gP2AycJmZfZF5oqT1zGxtsQPMh5l9LOkG4D1gJTDezJ5IOKxC28bMFkf3FxPW\ndVWqvsCopk7KZYV0F2DrHF6Xj/IYGMmTpE2A+wlfLJ8lHU+hSDqB0P85hTy/xCU9BVQBt0TdSnXW\nW0u6INq4camkhyVtGx2/StLN0f3Wkj6XNDh6vJGkL6Oxg/re95eSPpC0QFLfjKfvAA4AbjOzboTk\n9Vbaa2skXSRpFvB2dOwmSe9J+kTSK5IOi473IjT5T5P0qaQp0fFqSedF9yXpt5LmSVos6S5JbaPn\nOkbvd7akdyUtkTQgp//YtfHvAlwOdCS0aDeRdGY+1yxl0bbPFfmdI+k3wGozG9nUudmMOXwW/SP9\nVNIK4FGg2LulNrk4rtxJag08ANxtZg8lHU+BHQqcKGku4RfLUZKG5XIhMzsKeA64OOpWmpX+vKSj\ngGuAU4BtgXcJ3ZAA1YTEAnAQsJBQGx3gW8Bb9Q1CRl/YVwLHEH4cHZNxypfAMjN7OXr8MpBZ9PSk\n6D33jB5PIrQy2gEjgdGS1jezx6P47zGzTc2sa+qjU/uFdS7wo+iz7By91y0Z79cjivVo4PeSds/8\nXM1wIPAvpfT6AAAYa0lEQVSCmS2NWj0PEv5OK8liSe0Boh8THyYcT8FJOofQtZtVYm80OUgSsGf0\nj3TT6H/Gzmb2QP6hNktFL46L/jvfAUw3sxuTjqfQzGyAmXUws06EgcynzOzsPC+b2QJJfXGeSVgf\n85qZrSb8Cv+WpB2BFwn/jrYADif8N99eofTtkcAzDbzXqcBQM5sedRUNzHh+NfCxpC7R472AzC6l\nP5nZ8tQMHzMbYWbLzKzGzP4CbADslvbZGmthnQncEI3BfR59xtMlpf//fJWZrTKz14GphESUqxlA\n96h1JUJynJ7H9UrRI4SES/RnRf1Ai37g/BI4ycy+zOY12Q5IJyr6tZJaHDcduLfCFsf1IEwT7Clp\nSnTLnG1QSQrRZG/oGqnWQjgpfHkuBbaPSt2+QkgERxCSwQuE//6pxw1dM33M6716zhkJjJA0lTDY\nNz/j+a89lvQLSdMlLZe0jDDRY6sG3r++eN5Ne/weYfwwvZ98Udr9L4A2WV67DjObSphq/QrwenT4\nH7leL2mSRhH+3neTNF/SuYT9446VNBM4Knpclur5fH2BIYQW5sTo++W2pq7T6IC0mZmkVyUdbGaT\nChJ5jsxsHDAuyRjiYmb/pfR2yI2FmT1Dw1/ChfABoW8cgKhVsCW1s2ueIXS1dCV0/zxDmPZ3MPBs\nA9dcCOyY9njHjOc/Bz42s4Oi9+xHmGee7qtkJulwwq+4o8xsWnTsY9I2tmzOZ4ziWUsYSM2MrSDM\nbDBhn7WyZ2ZnNPBUZndhWWrg8zV7+UE2X0jdgf9JmiPpjej2epOvci5eyrifejwKOFfSfpI2IPTf\nv2hmqV/7zwBnA9OiadHVwPnAHDNb2sB73QecI2kPSRtTt1vpNeD/om6XXYHzmoh9U8KX+UeS1pf0\ne74+i2sR0DHqwqnPKOCKqJt1E2rHKGoaOB+Snc3lylA2yeE4YBdCU+s70e3EOINyLguWcd8AzOxJ\n4HeEwf0PCNNM0xds/Q/YkNpWwluE6ZkNtRqIBolvBJ4iLCB6MuP9/0oYd1hMWIB0dz3xpXs8us0k\nrMZdyde7qkZHfy6V9Eo9IQ0FhkcxzyF0G12S9nx9LY+KnH3j4hPb9hlqYEsISYMIv9SWRKcOiLqM\nUkv0+wLrgEvNbEIswTnnnGtUnMmhPdDezF6Lmr6TCasOTwU+jWZopJ+/J2FQ7yDCorcngC5NNJWd\nc87FILZB0Aa2hEiV/Kyv//MkYJSZrbFQQ3o2YZDQOedckRVlhoxqt4RI7VJ5iaSpku5IW5G6HV9f\n2Ob1o51zLiHZFPvJizK2hJD0N+AP0dNXAzfQ8OyOOn1e8hrSzjmXk+aUWY615VDflhBm9qFFCDt0\nprqOMrfI2IEGdn40s4q9DRw4MPEY/LP55/PPV/jbsmXGaacZ++xjTJtW/M/XXLElh4a2hIj2LUn5\nLqFwCITl66dH8747AZ0J+88451xZ+9//oGtX2GoreOkl2HPPpl+TtDi7lVJbQrye2lkSGEAo1LM/\noctoLvATADObLuk+wvYYa4GLLJd055xzJWLdOrjuOrjpJvjHP+CkMqpCE1tysIa3hGhwCwwzu4aw\n2rPFqqqqSjqE2FTyZwP/fOWu0J/v/ffhrLOgpgYmT4YdmiyvU1rKpoZ0iiRvUDjnStpjj8H558PF\nF8OAAdCqVdIRgSSsGQPSsc9Wcs65luLLL+HXv4aHH4YHHoAePZKOKHeeHJxzrgBmzIDTT4ddd4Up\nU6Bdu6Qjyk+L2CbaOefiYgZ33AGHHx66kUaPLv/EAN5ycM65nC1fDj/9KUyfDs88Ux5TVLPlLQfn\nnMtBOa5daA5vOTjnXDOU89qF5vDk4JxzWSr3tQvN4d1KzjmXhccegwMOgJ494cknKzsxgLccnHOu\nUZW0dqE5PDk451wDKm3tQnPEuStrB0lPS5om6U1Jl0bHt5A0UdJMSRPSiv0gqb+kWZJmSPp2XLE5\n51xjKnXtQnMkUUP6XOAjMxss6ddAOzPrl20Nad9byTkXp/S1C/fcUzlTVJu7t1ISNaRPBO6KTruL\nkDDAa0g75xKWWruw9dYwaVLlJIZcFGXMIa2G9EvANma2OHpqMbBNdH87amtMg9eQds4VSWrtws03\nw+23V+7aheYoVg3pBwg1pD8NBeICM7MmakJ7/5FzLlbpaxdeeaXyp6hmK9bkkFZDerhFNaSBxZLa\nm9miqGToh9HxrGtIDxo06Kv7VVVVFV+ExDkXj1TdhZ/9DPr3L426C4VSXV1NdXV1zq+Pc0BahDGF\npWZ2RdrxwdGx6yT1AzbPGJA+mNoB6V0zR599QNo5l6/0tQsjRrSMtQulVOynvhrS/YFrgfsknQfM\nA04FryHtnCuO1NqFzp1b3tqF5vAyoc65FsEMhg6Ffv3gmmtCd5Ky/h1d/kqp5eCccyWhkusuxMU3\n3nPOVbRKr7sQF285OOcqUkupuxAXTw7OuYrTkuouxMW7lZxzFaWl1V2Ii7ccnHMVoaXWXYiLJwfn\nXNlryXUX4uLdSs65suV1F+LjLQfnXFnytQvx8paDc67s+NqF+DWYHCSdEv25c/HCcc65hq1bF7a+\nOPlkuPFGuOUW2GijpKOqTA3urSRpipl1Tf1Z5Lga5HsrOdcyffBBWLuwbh3cfbdPUW2uQpYJXSpp\nItBJ0qMZt0eyDGaopMWS3kg7NkjSAklTotvxac/1lzRL0gxJ3872QzjnKttjj0G3blBV5WsXiqWx\nlsP6QDdgOHA+kJ5xzMyeafLi0uHAZ8AwM9snOjYQ+NTM/pJxbqqew0HU1nPoYmY1Ged5y8G5FmLV\nKvjVr1pW3YW4FHJX1jvM7CxJ/8wmEdTHzJ6L6kdnqi/Ak4BRZrYGmCdpNqHwz4v1nOucq3AzZsAZ\nZ8Auu/jahSQ01q10gKTtgB9K2iLzluf7XiJpqqQ7JG0eHdsOWJB2zgJCC8I514KsWAE33BDWLlx0\nka9dSEpjLYe/A08COwOTM56z6Hgu/gb8Ibp/NXADcF4D59bbf+Q1pJ2rPO+8A0OGwPDhcOyxvnYh\nX7HXkJb0dzP7ac5vELqVHk2NOTT0XFRPGjO7NnrucWCgmb2U8Rofc3CuQphBdXXYVvv55+G888JK\n5w4dko6s8hS8EpyZ/VTSfsARhF/yz5nZ1DwC3NbMFkYPvwukZjI9AoyU9BdCd1JnYFKu7+OcK11f\nfgmjRoW1CmvWwGWXhQHnNm2SjsylNJkcJF0GXAA8SBhIvjsapL45i9eOAo4EtpI0HxgIVEnan5Bo\n5gI/ATCz6ZLuA6YDa4GLvIngXGVZuBD+9je4/fawrfb114cupJZUy7lcZNOt9AbQ3cw+jx63AV6s\nr5uoGLxbybnyM3lyaCU89hj84AdwySWw++5JR9WyFHIRXLqaBu4751y91q6F++8Ps46++13Yd1+Y\nMwduvdUTQznIZlfWfwMvSUp1K50MDI01Kudc2Vq2LGyjPWRIGFi+/PKwF9J6vgd0WWmyWwlA0gHA\nYdQOSE+JO7BGYvFuJedK0Ntvw803w8iRcMIJYZD5wAOTjsqlFHy2EoCZTabuWgfnXAtnBhMnhvGE\nyZPhxz+GadNgu+2Sjszlyxt6zrlm++KLsFjt5puhVavQdfTAA759diXx5OCcy9qCBWFA+V//gkMP\nDeMKPXv6VNRK5JXgnHONWrUKJkyA008PM45WrgyV2B5+GI46yhNDpcpmEdz3gGuBbajdTdXMrG2c\ngTnnkrNgAYwdG25PPw177QWnnhoWr222WdLRuWLIZhHcO8AJZvZWcUJqnM9Wcq7w1q6FF18MyWDM\nmJAcjjsO+vQJf261VdIRunzFMVtpUakkBudc4SxZAuPHh2QwYQLsuCP07g233QaHHOLrElq6bFoO\nNwHtgYeA1dFhM7MHY46toXi85eBcDmpq4LXXQjIYOxamTw9jBn36wPHHw/ZePaWiNbflkE1yuDO6\n+7UTzezcLIIZCvQBPkwrE7oFcC+wEzAPONXMlkfP9Qf6AuuAS81sQj3X9OTgXJZWrAjrEMaMgXHj\noG3bkAx69w7bWmywQdIRumIpeHLIM5j6akgPBj4ys8GSfg20M7N+XkPaufyZhfKaqdbByy+Husu9\ne4fbrrsmHaFLSsHHHCT9O+OQAZhZ36Ze20AN6RMJ23gD3AVUA/3wGtLO5WTlylAwJ5UQ1qwJrYMr\nrgjdRl4jweUimyGnMdR2KW1EKNDzQR7vuY2ZLY7uLyZMkYVQQzo9EXgNaVc2Fi6E118Pv9yLZc6c\nkAyefRb23z8khIcfhr339rUHLn/ZVIK7P/2xpJHA84V4czMzSY397+Q1pF1JWrcOJk2qnfo5d24o\nXlPMGT7t28NZZ4VtLNq1K977uvIQew3pOi+QdgceM7Osei8za0hLmgFUmdkiSdsCT5vZ7l5D2pW6\njz+unfo5fjxsu23ox+/TB771LZ/66UpbHGMOn1H7C94IXUG/zi08INSK/hFwXfTnQ2nHvYa0Kxlm\noaso1Zf/+utQVRWSwTXXhHUBzlWqRlsOkgR0MLP3crp4Wg1pQlL5PfAwcB+wI3Wnsg4gTGVdC1xm\nZuPruaa3HFxsPvsMnniiduuIDTesnfp55JHhsXPlqKBTWaPk8HpS9aLr48nBFdrMmbVjBy++CN27\n13YXde7sg7uuMhS0WykaMH5V0sFm5l08riJ8+SU880xt6+CLL0IyuPhiePBB2HTTpCN0LnnZrJB+\nG9gVeBf4PDpsZrZvzLE1FI+3HFyzzZ9f2zqoroZ99qntLtpvP28duMoXx/YZHes7bmbzmhNYoXhy\nKE9m8M47sHp10+cWypIlYcuIsWPhgw+gV6+QDI47DrbcsnhxOFcKSmr7jDh4cigvX34J99wTagwv\nWVLcWgCbbgrHHhtaCAcfHMpZOtdSeXJwJWHRIvj738Nt//3hssvCL/ZveO1B5xLR3OTg/6u6gnr1\nVfjRj2CPPUKCeOopePzxsCW0Jwbnyof/7+rytm5dmOVzxBFw0kmw555hfOHvfw/3nXPlxxf8u5wt\nXw533AG33BK2krjsMvi//4PWrZOOzDmXL08OrtlmzYKbb4YRI8IMoHvuCWUlnXOVw7uVXFbMwrYS\nJ5wAhx4aKoq98QaMHOmJwblK5C0H16iVK+Huu+Gmm8Ljyy+H0aNho42Sjcs5F6/EkoOkecAKQr3o\nNWZ2cGP1pV1xvf8+3HYb/POfYY3AjTfC0Uf7SmLnWooku5WMUNehq5kdHB3rB0w0sy7Ak9FjV0ST\nJsEPfhCqia1YAf/9Lzz2GBxzjCcG51qSpMccMr9uTiTUlSb68+TihtMyrVkD994bxhJOPRUOPDBU\nNhsyBLp0STo651wSElshLWkO8AmhW+l2M/unpGVm1i56XsDHqcdpr/MV0gXy8ceh2+iWW6BTpzCe\ncOKJXtHMuUpU8EpwMephZgslbQ1MjMqHfiWL+tIuC+vWhR1JZ82qvc2cGf784AP43vdCUfpu3ZKO\n1DlXShJLDma2MPpziaT/AAcDiyW1T6sv/WF9rx00aNBX96uqqqiqqoo/4BJmFr7oM7/8Z82COXPC\nDqSdO4cuos6dQ6nLzp1h5529splzlaq6uprq6uqcX59It5KkjYFWZvappDbABOAq4BhgqZldJ6kf\nsLmZ9ct4bYvsVjILu5rWlwBmzYJNNqn98k/dunSBXXaBNm2Sjt45l7Sy2JVVUifgP9HD9YARZvan\naCprvfWl015b0clh2bL6v/xnzgxbTqcngPT7bdsmHblzrpSVRXLIR6Ulh2XLYMKEUJBmwgT4/PP6\nv/w7d/YCNc653HlyKHFm8OaboVzl2LEwZUrYzbRPn7BPUadOvp7AOVd4nhxK0Oefw5NP1ha0b9Wq\ntn5xz56+FYVzLn7lNJW1os2eXVvQ/oUX4KCDQjIYPx52391bB8650uYthwJZtQqee662u2jFipAM\nevcOdYx9wNg5lyTvViqi99+v7Sp66qlQ9ax379BltP/+XhbTOVc6PDnEaN06ePHF2u6i+fPhuONC\nQujVC7baKpGwnHOuSZ4cCuyjj8I4wZgx4c8OHWpbB4cc4vsQOefKgyeHPJmF6aWp7qJp08KMoj59\n4PjjYYcdYntr55yLjSeHHKxYEUpgjhkD48aFrShSU02POAI22KCgb+ecc0XnySELZvD227UziyZN\nCrUMUrOLOncuULDOOVciPDk0YOVKqK6uHUxes6a2dXDUUaG14JxzlarsF8FJ6gXcCLQC/mVm1+V6\nrXffrW0dPPtsmF7au3eoX7D33r4QzTnnGlJSM/EltQJuAXoBewJnSNoj29evWRNaB7/6Fey1Vyh3\n+eKLcNZZIVE8+yz06wf77FO6iSGf/ddLXSV/NvDPV+4q/fM1V0klB0LBn9lmNs/M1gD3ACc19oJF\ni+DOO+GUU+Cb34Rf/CLsVTR0aHhu2DA47TRo166xq5SOSv4HWsmfDfzzlbtK/3zNVWrdStsD89Me\nLwAOyTzppZdqp5rOng3HHBO6i4YMgfbtixarc85VrFJLDlmNNPftG5LB9ddDjx7QunXcYTnnXMtS\nUrOVJHUHBplZr+hxf6AmfVBaUukE7JxzZaRsp7JKWg94Gzga+ACYBJxhZm8lGphzzrUwJdWtZGZr\nJf0MGE+YynqHJwbnnCu+kmo5OOecKw2lNpW1QZJ6SZohaZakXycdTyFJ6iDpaUnTJL0p6dKkY4qD\npFaSpkh6NOlYCk3S5pLul/SWpOnR+FnFkNQ/+vf5hqSRksp2xzFJQyUtlvRG2rEtJE2UNFPSBEmb\nJxljPhr4fNdH/zanSnpQ0mZNXacskkO+i+PKwBrgCjPbC+gOXFxhny/lMmA6Wc5KKzM3AWPNbA9g\nX6BiukMldQQuALqZ2T6ELt/Tk4wpT/8mfJek6wdMNLMuwJPR43JV3+ebAOxlZvsBM4H+TV2kLJID\nOSyOKydmtsjMXovuf0b4Ytku2agKS9IOQG/gX0CJrk/PTfQr7HAzGwph7MzMPkk4rEJaQfgBs3E0\naWRj4P1kQ8qdmT0HLMs4fCJwV3T/LuDkogZVQPV9PjObaGY10cOXgCaLD5RLcqhvcdz2CcUSq+hX\nWlfCX2Al+SvwS6CmqRPLUCdgiaR/S3pV0j8lbZx0UIViZh8DNwDvEWYRLjezJ5KNquC2MbPF0f3F\nwDZJBhOzvsDYpk4ql+RQid0QdUjaBLgfuCxqQVQESScAH5rZFCqs1RBZD+gG3GZm3YDPKe9uia+R\ntAtwOdCR0KLdRNKZiQYVo2jb54r8zpH0G2C1mY1s6txySQ7vAx3SHncgtB4qhqTWwAPA3Wb2UNLx\nFNihwImS5gKjgKMkDUs4pkJaACwws5ejx/cTkkWlOBB4wcyWmtla4EHC32klWSypPYCkbYEPE46n\n4CSdQ+jazSqxl0tyeAXoLKmjpPWB04BHEo6pYCQJuAOYbmY3Jh1PoZnZADPrYGadCAOZT5nZ2UnH\nVShmtgiYL6lLdOgYYFqCIRXaDKC7pI2if6vHECYWVJJHgB9F938EVNQPtKgUwi+Bk8zsy2xeUxbJ\nIfq1klocNx24t8IWx/UAfgj0jKZ6Ton+MitVJTbZLwFGSJpKmK10TcLxFIyZTQWGEX6kvR4d/kdy\nEeVH0ijgBWA3SfMlnQtcCxwraSZwVPS4LNXz+foCQ4BNgInR98ttTV7HF8E555zLVBYtB+ecc8Xl\nycE551wdnhycc87V4cnBOedcHZ4cnHPO1eHJwTnnXB2eHJyLgaTNJF2Y9riqErcqd5XLk4Nz8WgH\nXJR0EM7lypODa/GibVlmRLuqvi1phKRvS3o+Kv5yUFQM5qGoWMr/JO0TvXZQVFzlaUnvSLokuuy1\nwC7RatTBhFXhm0gaHRVduTvt/a+NCulMlXR98f8LOFdXSdWQdi5BuwDfI2zP8jJwmpn1kHQiMICw\nZfxkMztZUk/CdhJdo9d2AXoCbYG3o60Jfk0ortIVQrdSdP6ewELgeUk9CPsWnWxmu0fntS3Gh3Wu\nKd5ycC6Ya2bTou2apwGpegVvEOo1HAYMBzCzp4EtJW1KaBGMMbM1ZraUsJvnNtS/NfkkM/sgeo/X\ngJ2A5cCXku6Q9F1gZXwf0bnseXJwLliVdr8GWB3dN0JZTKPhWhSr0+6vo+EW+aqM81qb2TpCpcP7\ngROAx5sXtnPx8OTgXHaeI9oHP+oiWmJmn9JwwvgU2LSpi0pqA2xuZuOAnwP7FSRa5/LkYw7OBZnb\nE1vG/auAodGW3J9Tu/d/vVXDzGxpNKD9BqEk49gG3mNT4GFJGxISzRX5fhDnCsG37HbOOVeHdys5\n55yrw5ODc865Ojw5OOecq8OTg3POuTo8OTjnnKvDk4Nzzrk6PDk455yrw5ODc865Ov4fj5I2Zx0G\nJg4AAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "%matplotlib inline\n", "from matplotlib.pyplot import plot, title, xlabel, ylabel, show, subplot\n", "flow=[80, 50 ,40 ,20 ,0 ,100, 150 ,200 ,250 ,120 ,100, 80]\n", "h=100; e=80\n", "\n", "subplot(211)\n", "plot(flow)\n", "title('hydrograph')\n", "xlabel('months')\n", "ylabel('run off,millon m**3/month' )\n", "fd=sorted(flow)\n", "\n", "subplot(212)\n", "plot(fd)\n", "title('flow duration')\n", "xlabel('months')\n", "ylabel('run off')\n", "\n", "t= range(0,12)\n", "for x in range(1,10):\n", " d=fd[x]\n", " ad=fd[(x-1)]\n", " if d==ad:\n", " t[x]=[]\n", " t[x-1]=t[x-1]+1\n", " fd[x]=[]\n", " #end\n", "#end\n", "ffw=fd+t\n", "print \"load duration data is as under\"\n", "print ffw\n", "mf=sum(flow)*10**6/(30*24*3600)\n", "print \"(a)\"\n", "print \"meanflow is %fm**3-sec\"%(mf)\n", "print \"(b)\"\n", "p=(735.5/75)*mf*h*e\n", "print \"power delevered in %dkW=%.3fMW\"%(p,p/1000)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## example 8.6 Page 156" ] }, { "cell_type": "code", "execution_count": 6, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "\n", "water used is \t\t31.709792m**3/sec\n", "\n", "effective head is \t200m\n", "\n", "power generated is \t53175.418569kW =\t53.175419MW\n", "\n", "peak load is \t\t70.900558Mw \n", "therefore the MW rating of station is \t70.900558MW\n", "\n", "for a head above 200m pelton turbine is suitable,\n", "francis turbine is suitable in the range of 30m-200m.,\n", "however pelton is most suitable\n" ] } ], "source": [ "mh=205#mean height\n", "a=1000*10**6#in miters\n", "r=1.25#annual rain fall\n", "er=0.8#efficiency\n", "lf=0.75#load factor\n", "hl=5#head loss\n", "et=0.9#efficiency of turbine\n", "eg=0.95#efficiency of generator\n", "wu=a*r*er/(365*24*3600)\n", "print \"\\nwater used is \\t\\t%fm**3/sec\"%(wu)\n", "eh=mh-hl\n", "print \"\\neffective head is \\t%dm\"%(eh)\n", "p=(735.5/75)*(wu*eh*et*eg)\n", "print \"\\npower generated is \\t%fkW =\\t%fMW\"%(p,p/1000)\n", "pl=p/lf\n", "print \"\\npeak load is \\t\\t%fMw \\ntherefore the MW rating of station is \\t%fMW\"%(pl/1000,pl/1000)\n", "if eh<=200:\n", " print \"\\nfor a head above 200m pelton turbine is suitable,\\nfrancis turbine is suitable in the range of 30m-200m.,\\nhowever pelton is most suitable\"\n", "else:\n", " print \"only pelton turbine is most suitable\"" ] } ], "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.9" } }, "nbformat": 4, "nbformat_minor": 0 }