{
 "metadata": {
  "name": "",
  "signature": "sha256:4f421cc7413696fe7c281427d95aff4b53c3ab081abe9cd6d31d1e2df1133e5a"
 },
 "nbformat": 3,
 "nbformat_minor": 0,
 "worksheets": [
  {
   "cells": [
    {
     "cell_type": "heading",
     "level": 1,
     "metadata": {},
     "source": [
      "    Chapter 15: Uniform Flows in Open Channels"
     ]
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 15.1, Page 516"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "from __future__ import division\n",
      "import math\n",
      "\n",
      "\n",
      " #Initializing  the  variables\n",
      "B  =4;\n",
      "D  =  1.2;\n",
      "C  =  7.6;\n",
      "n  =  0.025;\n",
      "s  =  1/1800;\n",
      "\n",
      " #Calculations\n",
      "W  =  B  +  2*1.5*D;\n",
      "A  =  D*(B+C)/2;                                                #  Area  of  parallelogram  formed\n",
      "P  =  B  +2*1.2*(D**2+(1.5)**2)**0.5;\n",
      "m  =A/P;\n",
      "i=s;\n",
      "C  =  (23+0.00155/i+1/n)/(1+(23+0.00155/i)*n/(m)**0.5);         #  By  Kutter  formula\n",
      "Q1  =  C*A*(m*i)**0.5;\n",
      "Q2  =  A*(1/n)*m**(2/3)*(i)**0.5;\n",
      "\n",
      "print \"Q using Chezy formula with C determined from the Kutter formula (m^3/s) :\",round(Q1,2)\n",
      "print \"Q using the Manning formula                                     (m^3/s) :\",round(Q2,2)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Q using Chezy formula with C determined from the Kutter formula (m^3/s) : 5.65\n",
        "Q using the Manning formula                                     (m^3/s) : 5.69\n"
       ]
      }
     ],
     "prompt_number": 1
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 15.2, Page 518"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "from __future__ import division\n",
      "import math\n",
      "\n",
      " #Initializing  the  variables\n",
      "Q  =  0.5;\n",
      "C  =  80;\n",
      "i  =  1/2000;\n",
      "\n",
      " #Calculations\n",
      "\n",
      "# A = D**2+(3/4)*D**2 = (7/4)*D**2\n",
      "D = ((4/7)*(Q/C)*(2/i)**0.5)**(2/5)\n",
      "\n",
      "#Result\n",
      "print \"Optimum depth = Optimum Width (in metres):\",round(D,3)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Optimum depth = Optimum Width (in metres): 0.552\n"
       ]
      }
     ],
     "prompt_number": 1
    }
   ],
   "metadata": {}
  }
 ]
}