{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Chapter 16 - Drying"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Page 909 Example 16.1"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "     mA = 0.113 kg/sec\n",
      "\n",
      " The total drying time is then 24.0 ksec or 6.65 hr \n"
     ]
    }
   ],
   "source": [
    "from math import log\n",
    "from __future__ import division\n",
    "#For the first drying operation\n",
    "w1 = 0.25#                 #in kg/kg\n",
    "w = 0.10#                  #in kg/kg\n",
    "wc = 0.15#                 #in kg/kg\n",
    "we = 0.05#                 #in kg/kg\n",
    "f1 = w1-we#                #in kg/kg\n",
    "fc = wc - we#              #in kg/kg\n",
    "f = w - we#                #in kg/kg\n",
    "\n",
    "#x = mA\n",
    "def TotalDryingTime(t):\n",
    "    x = (1/t)*((f1-fc)/fc + log(fc/f))\n",
    "    return x\n",
    "\n",
    "print\"\\n     mA = %.3f kg/sec\"%(TotalDryingTime(15))\n",
    "\n",
    "#For the second drying operation\n",
    "w12 = 0.30#                  #in kg/kg\n",
    "w2 = 0.08#                   #in kg/kg\n",
    "wc2 = 0.15#                  #in kg/kg\n",
    "we2 = 0.05#                  #in kg/kg\n",
    "f12 = w12 - we2#             #in kg/kg\n",
    "fc2 = wc2 - we2#             #in kg/kg\n",
    "f2 = w2 - we2#               #in kg/kg\n",
    "\n",
    "t1=(1/TotalDryingTime(15))*((0.25-0.10)/0.10 + log(0.10/0.03))\n",
    "print\"\\n The total drying time is then %.1f ksec or %.2f hr \"%(t1,t1*1000/3600)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Page 909 Example 16.2"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYYAAAEZCAYAAACTsIJzAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3X2clHW9//HXR1ZCEgUzrEDFFERMORYCHg/HSfnJ6jHw\nLgtJU9M2zdvKkI7JdszU08m8N354k6VBqYlaApU6ZiaIhngDKKQYYKIoelARwf2cP77Xwswwszu7\nO9de1+y8n4/HPJy5ru98r8+Oy3z2e3uZuyMiItJsq6QDEBGRdFFiEBGRPEoMIiKSR4lBRETyKDGI\niEgeJQYREcmjxCCJMbOsmX2tgvXdYmZvmtmcStVZUP8yMzukne+dYGazKx1TiWv93Mwu7oxrSdek\nxCCxir5M3zOztWb2avTl/dHotEeP1uoYYGZNZlby99XMRgGjgU+5+8jKRL+FsuIt+kb32919TIXj\nKXk5yoyz0slZugYlBombA0e4ey/gs8Aw4MJ21mUtnNsVWObu77e5UrO6dsaTW0e3jtZRYS19Vrm0\nwlW2oMQgncbdXwFmAXsXnrPgwqiFscrMbjWz7aLTf47++1bU8hhR8N6vAVOBA6Lzk6Pjp5nZEjN7\nw8zuMbNP5rynyczOMLMlwPPF4jWzE8zsZTNbbWbfKzjXaGZ3mtkvzext4AIze9fMdsgp81kze83M\nupnZSWb2SMH1G8zsBTNbY2bX5pzbysx+Ymavm9mLZnZmSy0mM9vPzP5mZv9rZtOBHjnn+pjZ76I4\n3jSz+8ysX3TuEmAUcG30uV0dHb/KzP5hZm+b2RNm9m/FritdlxKDdAYDMLOdgcOA+UXKnAx8FcgA\nnwa2BZq/LEdF/93e3Xu5+9zcN7r7TcA3gMei8z8ws4OBHwFfBD4JvAxML7jmOGB/YMgWAZsNAa4H\nJgCfAj4G9C8oNha4w923B34CZIHjcs6fAExz9w+L/LwA/0FoQe0LHGdmzV1NXwfqgaGEVtaRlPjL\n3sy6AzOAW4E+wB3AMTnlDbgJ2CV6rCP6XN39P4FHgG9Gn9vZ0Xsej67dB/gVcEd0HakRSgwSNwNm\nmNkawpdQlvCFXWgC8BN3X+bu7wKTgC9HfyWX0y1SWGYCcJO7P+XuH0T1HWBmu+SUudTd33L39UXq\nOxa4z93/Er3/+0BTQZm/uvu9AFEX1i+Ar8CmrqUvA79sIebL3P1/3X058BDhyxhCcrnS3V9x97eA\nS4v8fM1GAnXufpW7f+judwHzmk+6+5vufre7v+/u7xA++4MK6sirOxoPWePuTe5+BfARYM8Wfg7p\nYpQYJG4OjHP3Pu4+wN3PLPFF3PxXfbN/AHXATu28bl59UbJ5A+iXU2Z5K+9fkfP+96L351pR8Poe\nYIiZDQD+H/C2uz/RwjVezXn+HqGV1Hzt3NgKr5PrU8DKgmMvs7mV1tPMpkRddG8DDwPbm1luMshr\njZjZd8xsoZm9FSX07YEdW4hBuhglBkmLV4ABOa93ATYCq2jfAGlefdFMqI+R/yXaUr3/BHbOeX/P\n6P258t4ftRruILQavkJoQbRH3rULnhcr26/g2K45sX0bGAQMj7q8DiIkjebEUJgURgHnA190997u\n3gd4m/JabdJFKDFIWkwDzoumpm5L6PKY7u5NwOuEbpzd21jfyWY21Mw+EtU3x93/Ueb77wSOMLMD\no/71/6K8fy+/IIyXjKXlbqRCuV/WvwHOMbNPmVlvYCKlk9hfgY1mdraZbW1mRxPGTZptSxhXeDsa\nGJ9c8P5V5H+uvQgJebWZdTezi4DtkJqixCBpcTPhi/TPwIuErpWzYFM3ziXAo9EMnuFF3p83d9/d\nHyCMC9xFaD3sRujzzy1fkrsvBL5JGHx9BXiT/O6domsF3P1RQhJ7Mho7KFW+8L2556cCfwCeBp4E\nfg98GCXJwuttAI4GTiJ0dR1H+JmbXQlsA6wmJJGZBde+Cjg2mrF0JWHW2CzgBWAZIamUm0yli7A4\nb9RjZjcTZl685u77lChzNWGmynvASe5ebMaKSNUwsz8Bv3L3mytU32HADe4+oBL1ibQm7hbDLYRp\nd0WZ2eHAHu4+kDBF74aY4xGJlZntT5hi+usO1NHDzA43s7pozcFk4LeVilGkNbEmBnd/BFjTQpGx\nhPnXRHPTe5tZe2ehiCTKzG4F/gicG82CandVQCOh++pvwHPARR0OUKRMHd4KoIP6seW0vP6EATGR\nquLuX61QPeuAYuMoIp0iDYPPhdPgtHeLiEiCkm4xrCR/jnZ/tlysg5kpWYiItIO7t3kNStIthnuB\nEwHMbCTwlrsX7Ub67ncd9+p8TJ48OfEYajX+ao5d8Sf/qPb42yvWxGBm0whzp/c0s+Vmdkq0o2QD\ngLvfD7xoZkuBKcAZpeq68UZY0dLGACIiUhGxdiW5+/gyypxZTl2nnQY/+AFMndrxuEREpLSku5LK\nNnEizJgBixcnHUnbZTKZpEPokGqOv5pjB8WftGqPv71iXflcKWbm7s5//zfMnQt33dX6e0REap2Z\n4e0YfK6qxLBuHQwaBHfeCSNGtP4+EZFa1t7EUDVdSQDbbAOTJ4dupSrIZyIiVamqEgPASSfBqlUw\na1bSkYiIdE1Vlxjq6uCSS2DSJGjaYhNiERHpqKpLDABHHQU9esC0aUlHIiLS9VTV4HOuhx+Gk08O\n01e7d08oMBGRFKuJwedcBx0EgwfDlClJRyIi0rVUbYsBYMECGDMGliyBXr0SCExEJMVqrsUAMHQo\njB4NV1yRdCQiIl1HVbcYAF56CYYNg0WLoG/fTg5MRCTFamLlcylnnx3+e/XVnRSQiEgVqOnE8Npr\nsNdeMG8efPrTnRiYiEiK1eQYQ7O+fUOr4SLdLl1EpMO6RIsBYO1aGDgQZs8Og9IiIrUulS0GM6s3\ns8VmtsTMJhY538fM7jazBWY218z2bu+1evWC//zPsFWGiIi0X2yJwcy6AdcC9cAQYLyZ7VVQ7HvA\n39x9KOHez1d15JoNDWEldDbbkVpERGpbnC2G4cBSd1/m7huA6cC4gjJ7AQ8BuPvzwAAz+3h7L9i9\nO1x8sbblFhHpiDgTQz9gec7rFdGxXAuAowHMbDiwK9C/IxcdPx7Wr4e77+5ILSIitasuxrrL+Zv9\nMuAqM5sPPAPMBz4sVrCxsXHT80wmU/JerFttBZddBueeC2PHhm26RURqQTabJVuBvvTYZiWZ2Uig\n0d3ro9eTgCZ3v7yF97wE7OPu7xQcb3VWUi53OPhgmDABTj21ffGLiFS71C1wM7M64HngEOAV4HFg\nvLsvyimzPbDO3T8ws9OAA939pCJ1tSkxAMydC8ceCy+8EG4JKiJSa1I3XdXdNwJnArOBhcCv3X2R\nmTWYWUNUbAjwjJktBsYA51Tq+iNGwPDhcM01lapRRKQ2dJkFbsUsXgyjRoVWQ58+MQQmIpJiqWsx\npMHgwXDkkWEwWkREytOlWwwAK1fCvvuGm/r079BEWBGR6pK6wedK6khiALjgAnjjDZg6tYJBiYik\nnBJDC9asgUGD4JFHQveSiEgt0BhDC/r0gfPPD5vsiYhIy2qixQCwbl1oNdx5Z5jKKiLS1anF0Ipt\ntoHJk7XBnohIa2omMQCcdBKsWgWzZiUdiYhIetVUYqirg0suCTfzaWpKOhoRkXSqqcQAcNRR0KMH\nTJuWdCQiIulUM4PPuR5+GE4+OWyZ0b17xaoVEUkVDT63wUEHhfUMU6YkHYmISPrUZIsBwhYZY8bA\nkiXQq1dFqxYRSQW1GNpo6FAYPRquuCLpSERE0qVmWwwAL70Ew4bBokXQt2/FqxcRSZT2Smqns88O\n/7366liqFxFJTCq7ksys3swWm9kSM5tY5PyOZjbLzJ4ys2fN7KQ44ynmwgvh9tvhxRc7+8oiIukU\n5z2fuxHu+TwaWAnMY8t7PjcCH3H3SWa2Y1R+p+i2oLl1xdZiAPjBD8Ig9G23xXYJEZFOl8YWw3Bg\nqbsvc/cNwHRgXEGZfwLbRc+3A94oTAqd4Vvfgj/9KcxUEhGpdXEmhn7A8pzXK6JjuaYCe5vZK8AC\n4JwY4ympV6+wJfekSUlcXUQkXepirLucvp/vAU+5e8bMdgf+aGZD3X1tYcHGxsZNzzOZDJlMplJx\nAtDQAD/9KWSzUOGqRUQ6RTabJZvNdrieOMcYRgKN7l4fvZ4ENLn75Tll7gcucfdHo9cPABPd/YmC\numIdY2h2++1hdtKcOWBt7pUTEUmXNI4xPAEMNLMBZtYd+BJwb0GZxYTBacxsJ2BPILH5QePHw/r1\ncPfdSUUgIpK8WNcxmNlhwJVAN+Amd7/UzBoA3H1KNBPpFmAXQpK61N1/VaSeTmkxAMycCeedB88+\nG7bpFhGpVlrgViHucPDBMGECnHpqp1xSRCQWSgwVNHcuHHssvPBCuCWoiEg1SuMYQ9UaMQKGD4dr\nrkk6EhGRzqcWQwmLF8OoUaHV0KdPp15aRKQi1GKosMGD4cgj4bLLko5ERKRzqcXQgpUrYd99w1YZ\n/ft3+uVFRDpEg88xmTgR3nwTpk5N5PIiIu2mxBCTNWtg0CB45JHQvSQiUi00xhCTPn3g/PPDJnsi\nIrVALYYyrFsXWg133hmmsoqIVAO1GGK0zTYweXIYb6iCPCoi0iFKDGU66SRYtQpmzUo6EhGReCkx\nlKmuDi65JNzMp6kp6WhEROKjxNAGRx0FPXrAtGlJRyIiEh8NPrfRww/DySeHLTO6d086GhGR0jT4\n3EkOOiisZ5gyJelIRETioRZDOyxYAGPGwJIl0KtX0tGIiBSXyhaDmdWb2WIzW2JmE4uc/46ZzY8e\nz5jZRjPrHWdMlTB0KIweDVdckXQkIiKVF1uLwcy6Ac8T7um8EpgHjHf3RSXKHwGc6+6ji5xLVYsB\n4KWXYNgwWLQI+vZNOhoRkS2lscUwHFjq7svcfQMwHRjXQvnjgaqZ77PbbuH2nz/8YdKRiIhUVpyJ\noR+wPOf1iujYFsysJzAGuCvGeCruwgvh9tvhxReTjkREpHLqYqy7LX0/XwD+4u5vlSrQ2Ni46Xkm\nkyGTybQ7sErp2xfOOgsuughuuy3paESk1mWzWbLZbIfriXOMYSTQ6O710etJQJO7X16k7N3Ar919\neom6UjfG0GztWhg4EGbPDoPSIiJpkbr7MZhZHWHw+RDgFeBxigw+m9n2wItAf3dfV6Ku1CYGgGuu\ngZkz4f77k45ERGSz1A0+u/tG4ExgNrCQ0CJYZGYNZtaQU/RIYHappFANGhrCSugKtOBERBKnBW4V\ncvvtcPXVMGcOWJvzs4hI5aWuxVBrxo+H9evh7ruTjkREpGPUYqigmTPhvPPg2WfDNt0iIklSiyEF\n6uvhE5+An/886UhERNpPLYYKmzsXjj0WXngh3BJURCQpajGkxIgRMHx4mMIqIlKN1GKIweLFMGpU\naDX06ZN0NCJSq9RiSJHBg+HII+Gyy5KORESk7dRiiMnKlbDvvuGmPv37Jx2NiNSi1G2JUUnVmBgA\nJk6EN9+EqVOTjkREapESQwqtWQODBoVFb//2b0lHIyK1RmMMKdSnD/ziF3DMMXD99VCFuU1EapBa\nDJ1g6dKQHPbdF372M/joR5OOSERqgVoMKbbHHvDYY7DVVnDAAbBkSdIRiYiUpsTQSXr2DFtlnHEG\nHHgg3HNP0hGJiBSnrqQEzJ0Lxx0Hxx8PF1+sDfdEJB6p7Eoys3ozW2xmS8xsYokyGTObb2bPmlk2\nznjSYsQIeOKJ8BgzBl57LemIREQ2iy0xmFk34FqgHhgCjDezvQrK9AauA77g7p8Bjo0rnrT5+Mdh\n1iwYORKGDQutCBGRNIizxTAcWOruy9x9AzAdGFdQ5njgLndfAeDuq2OMJ3W6dYNLLoFrr4UvfEFT\nWkUkHeJMDP2A5TmvV0THcg0EdjCzh8zsCTM7IcZ4UmvsWPjrX2HKFDjxRHj33aQjEpFaFmdiKOdv\n362BzwKHA2OA75vZwBhjSi1NaRWRtIhzPsxKYOec1zsTWg25lgOr3X0dsM7M/gwMBbb4WmxsbNz0\nPJPJkMlkKhxu8pqntE6ZEqa0Tp0K4wo730RESshms2Sz2Q7XE9t0VTOrA54HDgFeAR4Hxrv7opwy\ngwkD1GOAjwBzgS+5+8KCurrUdNVyaEqriHRU6qaruvtG4ExgNrAQ+LW7LzKzBjNriMosBmYBTxOS\nwtTCpFCrNKVVRJKiBW4p9+GHcNFF8Mtfwm9+E6a3ioiUI3UtBqmM3CmtY8fCdddpSquIxEsthiqi\nXVpFpC3UYqgBzVNazTSlVUTio8RQZXr2hFtvhdNP1y6tIhIPdSVVMU1pFZGW6J7PNer110NiaGqC\nadOgb9+kIxKRtNAYQ40q3KV1zpykIxKRalcyMZjZd3Oef7Hg3I/iDEraRlNaRaSSSnYlmdl8d9+v\n8Hmx13FTV1L5NKVVRJqpK0kATWkVkY5TYuiCNKVVRDqipa6kD4H3opfbAOtyTm/j7p02OVJdSe03\ndy588YswYYKmtIrUGk1XlZI0pVWkNmmMQUrSlFYRaQslhhqhKa0iUi51JdUgTWkVqQ2p7Eoys3oz\nW2xmS8xsYpHzGTN728zmR48L44xHAk1pFZGWxJYYzKwb4X7O9cAQYLyZ7VWk6MPuvl/0+GFc8Ug+\nTWkVkVLibDEMB5a6+zJ33wBMB8YVKdfmZo5UhllIDPfdB2edBZMmwcaNSUclIkmLMzH0A5bnvF4R\nHcvlwL+a2QIzu9/MhsQYj5QwYgQ8+STMmwdjxsBrryUdkYgkKc7EUM5o8d+And19KHANMCPGeKQF\nH/84zJ4dkoSmtIrUtjjXwa4Eds55vTOh1bCJu6/NeT7TzK43sx3c/c3CyhobGzc9z2QyZDKZSsdb\n87p1gx/9KCSHsWNh8mQ444zQ5SQi6ZfNZslmsx2uJ7bpqmZWBzwPHAK8AjwOjHf3RTlldgJec3c3\ns+HAb9x9QJG6NF21kzVPad1nH7jmGujTJ+mIRKStUjdd1d03AmcCs4GFwK/dfZGZNZhZQ1TsWOAZ\nM3sKuBL4clzxSNs0T2ndbjsYPDgsiNPAtEht0AI3adXTT8N558Grr8IVV4QBahFJP22iJ7FyD9Na\nv/1tGDQIfvKT0JIQkfRKXVeSdC1mYUD6uefgkENg1Cg45xx4c4tpAiJS7ZQYpE26d4dvfQsWLoQP\nPgithmuugQ0bko5MRCpFXUnSIc88ExLFihVh/OGww5KOSESaaYxBEuMOv/99SBB77BHGH/YqtiuW\niHQqjTFIYszgiCPg2Wfh0EPh3/8dzj4b3ngj6chEpD2UGKRiuneHc8+FRYvCbUT32guuukrjDyLV\nRl1JEpvnngvdSy+/HLqXDj9c22uIdCaNMUgqucPMmSFB7LprGKDee++koxKpDRpjkFQyCy2FZ54J\n/81k4JvfhNWrk45MREpRYpBOsfXWYUHc4sWw1VZh/OGnPw1rIUQkXZQYpFN97GNhQdzDD8Mf/hB2\nb/3d70KXk4ikg8YYJFHN4w877xzGHz7zmaQjEuk6NMYgVemww8LurWPHwsEHhxsDvf560lGJ1DYl\nBknc1lvDmWeG8Yfu3WHIkDC9VeMPIslQYpDU2GEHuPJKeOQReOihMK31nns0/iDS2WJNDGZWb2aL\nzWyJmU1sodz+ZrbRzI6OMx6pDoMHhwHpa6+F730PRo8O3U0i0jliSwxm1g24FqgHhgDjzWyLrdWi\ncpcDswCti5VNxoyBBQvg6KNDcmhogNdeSzoqka4vzhbDcGCpuy9z9w3AdGBckXJnAXcCGnKULdTV\nhQVxzz8PPXuG8Yf/+R9Yvz7pyES6rjgTQz9gec7rFdGxTcysHyFZ3BAdUm+yFNWnT1gQ9+ijYQ3E\n3nvDjBkafxCJQ12MdZfzT/ZK4AJ3dzMzWuhKamxs3PQ8k8mQyWQ6Gp9UoT33DPee/uMf4bzz4Oqr\nQ8IYOjTpyESSl81myWazHa4ntgVuZjYSaHT3+uj1JKDJ3S/PKfMim5PBjsB7wGnufm9BXVrgJlvY\nuBFuvBEaG8M6iIsvhp12SjoqkfRI4wK3J4CBZjbAzLoDXwLyvvDd/dPuvpu770YYZzi9MCmIlFJX\nB9/4Rlj/sN12oXvp8svh/feTjkykusWWGNx9I3AmMBtYCPza3ReZWYOZNcR1Xak9vXuHAenHHguP\nvfeG3/5W4w8i7aW9kqTLeeCBMP6www7w4x/DsGG6QZDUJt2oRyTHhx/CTTfBj34E774LI0dufgwf\nDr16JR2hSPyUGERKWLEC5s4N3Uxz5sD8+bD77nDAASFRHHAADBoU7hMh0pUoMYiU6YMPworqOXM2\nJ4s1a2DEiM3JYsSIMHYhUs2UGEQ64NVX81sVTz4Z7hHRnChGjgyrrrt1SzpSkfIpMYhU0MaN4T7V\nua2KVatg//3zWxU77ph0pCKlKTGIxGz16tCqaE4Wjz8On/jE5nGKkSPDrUrr4txPQKQNlBhEOtmH\nH8LChSFRNCeL5cvhc5/LTxZajS1JUWIQSYG33spvVcydGzYAzE0UQ4eGO9WJxE2JQSSFmprghRc2\nj1M89hi8+CL8y7/kJ4t+/VqvS6StlBhEqsTatTBvXn6y2Gab/HUV++0HPXokHalUOyUGkSrlDn//\n++ZEMWdO2Bhwn302T5U94ADYZRdt7SFto8Qg0oW8+25YS5HbqjDLb1V87nPhrnYipSgxiHRh7vDy\ny/nrKp59FvbaK79V8elPq1UhmykxiNSYdevCvk+5yWL9+i03DNx226QjlaQoMYgIK1bkJ4qnnoI9\n9sjvgho4UBsG1golBhHZwgcfhOSQmyzefnvLDQO33z7pSCUOqUwMZlYPXAl0A27Mvd9zdH4c8F9A\nU/Q4390fLFKPEoNIhbz6av5q7SefhF133XLDQLUqql/qEoOZdQOeB0YDK4F5wHh3X5RT5qPu/m70\nfB/gbnffo0hdSgwiMdmwYcsNA197LYxP5CaLHXZIOlJpqzQmhgOAye5eH72+AMDdL2uh/E/dfWSR\nc0oMIp1o9er8VsW8efDJT+av1v7MZ7RhYNq1NzHE+b+1H7A85/UKYERhITM7ErgU+CRwaIzxiEiZ\ndtwRjjgiPGDzhoHNLYqrrgoD3cOG5SeLvn2TjVsqI87EUNaf+O4+A5hhZqOAXwJ7FivX2Ni46Xkm\nkyGTyXQ8QhEpS7duYSX2PvvA178ejq1Zs3nDwOuvh69+NXQ3FW4YuPXWycZeS7LZLNlstsP1xNmV\nNBJozOlKmgQ0FQ5AF7zn78Bwd3+j4Li6kkRSrqkJnn8+f7X2qlXwla/A174Ge++ddIS1J41jDHWE\nwedDgFeAx9ly8Hl34EV3dzP7LHCHu+9epC4lBpEqtHQp3Hwz/PznMGBASBBf+pIW3XWW1CUGADM7\njM3TVW9y90vNrAHA3aeY2XeBE4ENwDvAt9x9XpF6lBhEqtjGjTBzJtx4I/z5z3DMMXDqqWENhbbw\niE8qE0OlKDGIdB3//CfceivcdBN85COhFXHCCbp/dhyUGESkqriH1sONN8J998Ghh4ZWxOjRWlxX\nKUoMIlK13noLfvWr0IpYvRpOOQVOPjncg0LaT4lBRLqE+fNDgpg2DfbfP3Q1jR0bup2kbZQYRKRL\nWbcO7rorJInnntO01/ZQYhCRLkvTXttHiUFEujxNe20bJQYRqSm50167dw8JQtNe8ykxiEhN0rTX\n0pQYRKTmadprPiUGEZEcmvaqxCAiUlQtT3tVYhARaUXutNdddw1jEV152qsSg4hImWpl2qsSg4hI\nO3Tlaa9KDCIiHVBs2uvXvw6HHFK9rQglBhGRCmme9vqzn8H69XDGGeGe1r17Jx1Z27Q3McS+/MPM\n6s1ssZktMbOJRc5PMLMFZva0mT1qZvvGHZOISEt69w7JYMGC0MU0dy7sthucdlqYBtvVxX1rz26E\n+z6PBlYC89jyvs8HAAvd/W0zqwca3X1kQT1qMYhIol59NSSJKVOgXz/45jfh2GOhR4+kIystlV1J\n0Zf+ZHevj15fAODul5Uo3wd4xt37FxxXYhCRVNi4EX7/e7j++tB6OOUUaGgILYq0SWtXUj9gec7r\nFdGxUr4G3B9rRCIiHVBXB+PGwezZ8Je/wAcfwLBhcMQRcP/90NSUdIQdVxdz/WX/mW9mnwdOAQ4s\ndr6xsXHT80wmQyaT6WBoIiIdM2gQXHEF/PCHMH06fP/7cOaZcPrpYY+mzp7yms1myWazHa4n7q6k\nkYQxg+aupElAk7tfXlBuX+C3QL27Ly1Sj7qSRCT13GHePLjuOrjnntCyOOMMGD48mSmvaR1jqCMM\nPh8CvAI8zpaDz7sADwJfcfc5JepRYhCRqrJ6NdxyC9xwA+ywQ0gQX/4y9OzZeTGkMjEAmNlhwJVA\nN+Amd7/UzBoA3H2Kmd0IHAX8I3rLBncfXlCHEoOIVKWmpjAecd11MGcOnHhi6GoaODD+a6c2MVSC\nEoOIdAUvvRSmu958M+y3X2hF/Md/hAHtOCgxiIhUifffhzvvDFNeV6wI011PPRV22qmy10nrdFUR\nESnQo0e4L8Rf/wozZsCyZTB4MBx/fJgCm/TfwWoxiIikwFtvhV1er78+3GXujDNgwgTo1av9daor\nSUSkC3CHBx4ICSKbDcnh9NNhyJC216WuJBGRLsAMRo+G3/42bOLXp0/Y+vvzn4c77oANGzohhmr4\nS1wtBhGpZR98AHffHVoRS5aE+0ScdlrYzK8lajGIiHRR3buHe1M//DD84Q/w+uuwzz5hd9cHH6z8\nYLVaDCIiVWjtWrjttrBw7sMPw2D1iSfC9ttvLqPBZxGRGuQOjzwSuplmz4bjjgtJYuhQJQYRkZr3\n6qvhntVTpsCuu8KjjyoxiIgI4WZC990HRx+txCAiIjk0K0lERCpCiUFERPIoMYiISJ7YE4OZ1ZvZ\nYjNbYmYTi5wfbGaPmdn7ZvbtuOMREZGWxZoYzKwbcC1QDwwBxpvZXgXF3gDOAv4nzliSVImbcyep\nmuOv5thB8Set2uNvr7hbDMOBpe6+zN03ANOBcbkF3P11d38C6IStoZJR7b9c1Rx/NccOij9p1R5/\ne8WdGPqwNFr4AAAGVUlEQVQBy3Ner4iOiYhISsWdGLT4QESkysS6wM3MRgKN7l4fvZ4ENLn75UXK\nTgbecfefFDmnBCMi0g7tWeBWF0cgOZ4ABprZAOAV4EvA+BJlSwbfnh9MRETaJ/YtMczsMOBKoBtw\nk7tfamYNAO4+xcw+AcwDtgOagLXAEHd/J9bARESkqKrYK0lERDpPqlY+t7YYLipzdXR+gZnt19kx\ntqSaF/OVEfuE6DN/2sweNbN9k4izlDLiHxfFP9/MnjSzg5OIs5Ryfvejcvub2UYzO7oz42tNGZ9/\nxszejj7/+WZ2YRJxllLmd08miv1ZM8t2cogllfHZfyfnc38m+v3p3WKl7p6KB6GraSkwANgaeArY\nq6DM4cD90fMRwJyk425j/B8HhgE/BL6ddMxtjP0AYPvoeX0VfvYfzXm+D2F9TeKxlxt/TrkHgd8B\nxyQddxs//wxwb9KxdiD+3sBzQP/o9Y5Jx92W352c8kcAf2qt3jS1GFpdDAeMBW4FcPe5QG8z26lz\nwyypmhfzlRP7Y+7+dvRyLtC/k2NsSTnxv5vzcltgdSfG15pyfvch7BBwJ/B6ZwZXhnLjT+skknLi\nPx64y91XALh7Wn5/yv3smx0PTGut0jQlhnIWwxUrk5YvqGpezNfW2L8G3B9rRG1TVvxmdqSZLQJm\nAmd3UmzlaDV+M+tH+Ad/Q3QoTYOD5Xz+Dvxr1J13v5kN6bToWldO/AOBHczsITN7wsxO6LToWlb2\nv10z6wmMAe5qrdK4p6u2Rbm/6IV/daTlH0ha4miPsmM3s88DpwAHxhdOm5UVv7vPAGaY2Sjgl8Ce\nsUZVvnLivxK4wN3dzIx0/fVdTvx/A3Z29/eimYozgEHxhlW2cuLfGvgscAjQE3jMzOa4+5JYI2td\nW753vgD8xd3faq1gmhLDSmDnnNc7E7JfS2X6R8fSoJz406qs2KMB56lAvbuv6aTYytGmz97dHzGz\nOjP7mLu/EXt0rSsn/s8B00NOYEfgMDPb4O73dk6ILWo1fndfm/N8ppldb2Y7uPubnRRjS8r5/JcD\nq919HbDOzP4MDAWSTgxt+d3/MmV0IwGpGnyuA/5OGETpTuuDzyNJ1wBoq/HnlG0kXYPP5Xz2uxAG\nuUYmHW8749+dzdOzPwv8Pem42/O7E5W/BTg66bjb+PnvlPP5DweWJR13G+MfDPyJMNjbE3iGsN4q\n9bFH5bYn7GS9TTn1pqbF4O4bzexMYDabF8Mtyl0M5+73m9nhZrYUeBc4OcGQ85QTf+FiPjM7hxQs\n5isnduAioA9wQ/RX6wZ3H55UzLnKjP8Y4EQz2wC8Q/jrKRXKjD+1yoz/WOB0M9sIvEeVff7uvtjM\nZgFPExbiTnX3hclFHbThd+dIYLaHFk+rtMBNRETypGlWkoiIpIASg4iI5FFiEBGRPEoMIiKSR4lB\nRETyKDGIiEgeJQbp0sysQ2tEzOznZnZMkePnmtk2BccuMLPjzWxyqW3VLdzCtuQxM9s52o/nuWh7\n5zTt6SQ1IjUL3ERi0tGFOqXefw5hv6XcBUOHAscRNlzLY2YTgE8CPczsfMKtbily7AHgPHd/ysy2\nBZ40sz+6+6IO/hwiZVOLQaR1DmBmF5vZLdFf8Z8CHjKzB6Jz2wHdvcR2zO5+O2Ffm/OBl9399hLH\nXnX3p6L3vAMsiq4l0mnUYpCaE22A1qvIqW+7+4PF32I/Jtzs5+TowHlAxjdvAjeasJdOqWuOJ2yH\n/GNg1+g1hcfcfVrOewYA+xHufyHSaZQYpOa4+7+3obgB3wfmuntDC+XGADe3cM1pEMYT3P3HeRco\nfmxbwk15zkl6Ly2pPUoMUnPM7BHCXdwKfcfdHyg45oSNDz9nZn289Hbjw4FvtHZtd/9Ba8fMbGvC\nzVRu83APCZFOpcQgNcfdR7XxLbMIu1f+3swOjf6CX0vYJfdNM9sbWOwV2JEyugnPTcBCd7+yo/WJ\ntIcGn6XLMrM6YH0FqnJ3v5Nwk6J7zawH8P+BWWb2IFBPuF1orgvNbHn0+EcbrnUg8BXg82Y2P3rU\nV+BnECmbtt2WLsvMhgJT3H1kzNf5A3CCu6+K8zoinUVdSdIlmdk3gLMI6w1i5e6Hxn0Nkc6kFoOI\niOTRGIOIiORRYhARkTxKDCIikkeJQURE8igxiIhIHiUGERHJ83+n3mjxfJlAtgAAAABJRU5ErkJg\ngg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7f0dd8102610>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      " k = 0.520000\n",
      "\n",
      " t1 = 1090 ksecs or 12.6 days\n"
     ]
    }
   ],
   "source": [
    "%matplotlib inline\n",
    "from matplotlib.pyplot import plot,title,show,xlabel,ylabel\n",
    "from sympy import symbols,solve\n",
    "\n",
    "E = [1, 0.64 ,0.49, 0.38, 0.295, 0.22, 0.14]\n",
    "J = [0, 0.1, 0.2, 0.3, 0.5, 0.6 ,0.7]\n",
    "plot(J,E)\n",
    "title(\"Plot for drying data\")\n",
    "xlabel(\"J = kt/L**2\")\n",
    "ylabel(\"E\")\n",
    "show()\n",
    "\n",
    "#For the 10 mm strips\n",
    "mi = (0.28 - 0.07)#        #Initial free moisture content in kg/kg\n",
    "mf = (0.13-0.07)#          #Final free moisture content in kg/kg\n",
    "#at\n",
    "t = 25#                    #time is in ksecs\n",
    "E = (0.06/0.21)\n",
    "#at E = 0.286 ,J = 0.52 from the plot of given data and J = kt/L**2\n",
    "k = symbols('k')\n",
    "k1 = solve(0.52 - (k*t)/(10/2)**2)[0]\n",
    "print\"\\n k = %f\"%(k1)\n",
    "\n",
    "#for the 60 mm planks\n",
    "m1i = (0.22 - 0.07)#       #Initial free moisture content in kg/kg\n",
    "m1f = (0.13 - 0.07)#       #Final free moisture content in kg/kg\n",
    "E = (m1f/m1i)\n",
    "#at E = 0.20 from yhe plot of the given data J = 0.63 and J = kt/L**2\n",
    "t1 = 0.63*(60/2)**(2)/k1\n",
    "print\"\\n t1 = %d ksecs or %.1f days\"%(t1,(t1*1000/(3600*24)))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Page 911 Example 16.3"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      " The heat content of this stream = 6.05 kW\n",
      "\n",
      " Heat in the water vapour = 119.3 kW\n",
      "\n",
      " The dried solids contain 0.0036 kg/s water and hence heat content of this stream is: 0.07 kW\n",
      "\n",
      " G = 2.07 kg/secs\n",
      "\n",
      " Water in the outlet stream 0.0588 kg/secs\n",
      "\n",
      " The humidity H = 0.0284 kg/kg dry air\n"
     ]
    }
   ],
   "source": [
    "from sympy import symbols,solve\n",
    "\n",
    "#(a) Air\n",
    "#G kg/s dry air enter with 0.006G kg/s water vapour and hence the heat content of this stream= \n",
    "#[(1.00G) + (0.006G × 2.01)](385 − 273) = 113.35G kW\n",
    "\n",
    "#(b) Wet solid\n",
    "#0.125 kg/s enter containing 0.40 kg water/kg wet solid, assuming the moisture is expressed on a wet basis.\n",
    "flowWater = (0.125*0.40)#               #in kg/sec\n",
    "flowDrySolid = (0.125-0.050)#           #in kg/sec\n",
    "#Hence heat content of this stream\n",
    "q = ((0.050*4.18)+(0.075*0.88))*(295-273)\n",
    "print\"\\n The heat content of this stream = %.2f kW\"%(q)\n",
    "\n",
    "#Heat out \n",
    "#(a) Air\n",
    "#Heat in exit air = [(1.00 G) + (0.006 G × 2.01)](310 − 273) = 37.45G kW.\n",
    "fd = 0.075#               #mass flow rate of dry solids in kg/sec\n",
    "w = 0.05*0.075/(1+0.05)#  #water in the dried solids leaving in kg/secs\n",
    "we = (0.050 - w)#         #The water evaporated into gas stream in kg/secs\n",
    "\n",
    "#Assuming evaporation takes place at 295 K,then:\n",
    "qout = 0.0464*(2.01*(310-295)+2449+4.18*(295-273))\n",
    "print\"\\n Heat in the water vapour = %.1f kW\"%(qout)\n",
    "\n",
    "#the total heat in this stream = (119.30 + 37.45G) kW.\n",
    "#(b) Dried solids\n",
    "\n",
    "\n",
    "#The dried solids contain 0.0036 kg/s water and hence heat conten     t of this stream is:\n",
    "q1 = ((0.075*0.88)+(0.0036*4.18)/(305-273))\n",
    "print\"\\n The dried solids contain 0.0036 kg/s water and hence heat content of this stream is: %.2f kW\"%(q1)\n",
    "   \n",
    "   \n",
    "#(c) Losses\n",
    "#These amount to 20 kJ/kg dry air or 20m kW.\n",
    "#Heat Balance\n",
    "G = symbols('G')\n",
    "G1 = solve(113.35*G + 6.05 - 119.30 - 37.45*G - 2.59 - 20*G)[0]\n",
    "print\"\\n G = %.2f kg/secs\"%(G1)\n",
    "print\"\\n Water in the outlet stream %.4f kg/secs\"%(0.006*2.07+0.0464)\n",
    "print\"\\n The humidity H = %.4f kg/kg dry air\"%(0.0588/2.07)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Page 922 Example 16.4"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      " water in the product ,b = 12.8 kg\n",
      "\n",
      " As the final moisture content is above the critical value, all the drying is at this constant rate and the time of drying is: 11700 ksecs or 3.25 hr\n"
     ]
    }
   ],
   "source": [
    "from sympy import symbols,solve\n",
    "\n",
    "#In 100 kg of feed\n",
    "\n",
    "#mass of water = \n",
    "mw = 100*30/100#                  #mass of water in kg\n",
    "#mass of dry solids = \n",
    "md = 100-30#                      #mass of dry solids\n",
    "   \n",
    "#and:\n",
    "#For b kg water in the dried solids: 100b/(b + 70) = 15.5\n",
    "b = symbols('b')\n",
    "b1 = solve(100*b - 15.5*(b+70))[0]\n",
    "print\"\\n water in the product ,b = %.1f kg\"%(b1)\n",
    "\n",
    "#Initial water content \n",
    "w1 = 30/70#            #Intial moisture content in kg/kg dry solids\n",
    "#Final moisture content\n",
    "w2 = (12.8/70)#        #Final moisture content in kg/kg dry solids\n",
    "#water to be removed\n",
    "w3 = (30-12.8)#        #water to be removed in kg\n",
    "\n",
    "#Surface of drying\n",
    "S = (0.03*70)#            #Surface for drying in m**2\n",
    "rate = (0.0007*2.1)#      #Rate of drying in kg/sec\n",
    "\n",
    "#As the final moisture content is above the critical value, all the drying is at this constant rate and the time of drying is:\n",
    "t = 17.2/0.00147\n",
    "print\"\\n As the final moisture content is above the critical value, all the drying is at this constant rate and the time of drying is: %d ksecs or %.2f hr\"%(t,t/3600)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Page 930 Example 16.5"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      " For adiabatic drying the outlet air temperature = 331.6 K\n",
      "\n",
      " Total heat = 419.5 kJ/kg or 146 kW\n",
      "\n",
      " Mass flow rate of inlet air = 2.05 kg/secs\n",
      "\n",
      " Mass flow rate of dry air ,Ga = 2.03 kg/secs\n",
      "\n",
      " the humidity of the outlet air H2 = 0.0348 kg/kg\n",
      "\n",
      " The cross sectional area of the dryer is 2.18 m**2\n",
      "\n",
      " The equivalent diameter of the dryer = 1.67 m\n",
      "\n",
      " Logarithmic mean temperature difference = 45.5 deg K\n",
      "\n",
      " The length of the dryer = 10.1 m\n",
      "\n",
      " Length/diameter ratio = 6 \n"
     ]
    }
   ],
   "source": [
    "from sympy import symbols,solve\n",
    "from math import exp,log,pi\n",
    "\n",
    "#Inlet air temperature\n",
    "Tair = 400#            #Inlet air temperature in kelvins\n",
    "H = 0.01#              #Humidity is in kg/kg dry air\n",
    "#*therefore wet bulb temperature = \n",
    "Twetbulb = 312#        #Inlet wet-bulb temperature\n",
    "NTU = 1.5#              #Number of transfer units\n",
    "\n",
    "#Then for adiabatic drying the outlet air temperature,To is given by\n",
    "To = symbols('To')\n",
    "To1 = solve(exp(1.5)*(To-312)-(400-312))[0]\n",
    "print\"\\n For adiabatic drying the outlet air temperature = %.1f K\"%(To1)\n",
    "\n",
    "#Solids outlet temperature will be taken to be maximum allowable,325K\n",
    "#From the steam tables in the Appendix, the latent heat of vaporisation of water at 312 K is2410 kJ/kg. Again from steam tables, the specific heat capacity of water vapour = 1.88 kJ/kg K and that of the solids will be taken as 2.18 kJ/kg K.\n",
    "\n",
    "#For a mass flow of solids of 0.35 kg/s and inlet and outlet moisture contents of 0.15 and 0.005 kg/kg dry solids respectively, the mass of water evaporated = 0.35(0.15 − 0.005) = 0.0508 kg/s.\n",
    "\n",
    "#For unit mass of solids ,the heat duty includes:\n",
    "       #Heat to the solids\n",
    "qsolids = 2.18*(325-300)#       #heat to solids in kJ/kg\n",
    "#Heat to raise the moisture to the dew point\n",
    "qdew = (0.15*4.187*(312-300))#  #in kJ/kg\n",
    "#Heat of vaporisation\n",
    "qvap = 2410*(0.15-0.005)#       #in kJ/kg\n",
    "#Heat to raise remaining moisture to solids outlet temperature\n",
    "qremaining = (0.05*4.187)*(325-312)\n",
    "#Heat to raise evaporated moisture to the air outlet temp.\n",
    "qevapo = (0.15-0.005)*1.88*(331.5-312)\n",
    "qtotal = qsolids + qdew + qvap + qremaining + qevapo\n",
    "print\"\\n Total heat = %.1f kJ/kg or %d kW\"%(qtotal,qtotal*0.35)\n",
    "\n",
    "#The humid heat of entering air is 1.03 kJ/kg K\n",
    "#G1 (1 +H1) = Q/Cp1(T1 − T2)\n",
    "#where: G1 (kg/s) is the mass flowrate of inlet air,\n",
    "#H1 (kg/kg) is the humidity of inlet air,\n",
    "#Q (kW) is the heat duty,\n",
    "#Cp1 (kJ/kg K) is the humid heat of inlet air\n",
    "#and: T1 and T2 (K) are the inlet and outlet air temperatures         respectively.\n",
    "G1 = symbols('G1')\n",
    "G = solve(G1*(1+0.01)-146/(1.03*(400-331.5)))[0]\n",
    "print\"\\n Mass flow rate of inlet air = %.2f kg/secs\"%(G)\n",
    "print\"\\n Mass flow rate of dry air ,Ga = %.2f kg/secs\"%(G/1.01)\n",
    "print\"\\n the humidity of the outlet air H2 = %.4f kg/kg\"%(0.01+0.0508/2.05)\n",
    "\n",
    "#At a dry bulb temperature of 331.5 K, with a humidity of 0.0347 kg/kg, the wet-bulb temperatureof the outlet air, from Figure 13.4 in Volume 1, is 312 K, the same as the inlet, which is the case for adiabatic drying.\n",
    "\n",
    "#The dryer diameter is then found from the allowable mass velocity of the air and the entering air flow and for a mass velocity of 0.95 kg/m**2.secs, the cross sectional area of the dryer is\n",
    "print\"\\n The cross sectional area of the dryer is %.2f m**2\"%(2.07/0.95)\n",
    "print\"\\n The equivalent diameter of the dryer = %.2f m\"%(((4*2.18)/pi)**0.5)\n",
    "\n",
    "#With a constant drying temperature of 312 K:\n",
    "#at the inlet\n",
    "deltaT1 = 400-312#           #inlet temperature is in deg K\n",
    "#at the outlet\n",
    "deltaT2 = 331.5-312#         #outlet temperature is in deg K\n",
    "Tlogmean = (deltaT1 - deltaT2)/log(deltaT1/deltaT2)\n",
    "print\"\\n Logarithmic mean temperature difference = %.1f deg K\"%(Tlogmean)\n",
    "#The length of the dryer, L is then: L = Q/(0.0625πDG**(0.67)*\u001bTm)\n",
    "#where D (m) is the diameter\n",
    "#and G(kg/m**2.secs) is the air mass velocity.\n",
    "L = 146/(0.0625*(pi)*1.67*(0.95)**(0.67)*45.5)\n",
    "print\"\\n The length of the dryer = %.1f m\"%(L)\n",
    "print\"\\n Length/diameter ratio = %d \"%(10.1/1.67)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Page 950 Example 16.6"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYQAAAEZCAYAAACXRVJOAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XecFeX59/HPtbt0pFsQpdhrBGVpIi5qjI0iXWwoxLKb\nmJj42JI8YowxMYmJ8begjyAISrUBKiqaLIIFFkVRLIiKiiK9ShHY6/ljZvH81i1nlz07Z89+36/X\nvJg+17l3ONeZ+565x9wdERGRtKgDEBGR5KCEICIigBKCiIiElBBERARQQhARkZASgoiIAEoIIknP\nzIaZ2byo45DUp4RQw5nZNjPbGg4FZrY9ZvqSqOOrCDNbYWZnRR1HeZhZ27D8t8YMi6OOqyRmNtLM\nJkYdh1SujKgDkGi5e8PCcTP7HBju7v+JMKRSmVm6u+8tYzUHLMHHSJTG7l4Q0bEBMLMMd98TZQwS\nDV0hSLHMLM3MbjWz5Wa2zsymmlnTcFnhr9lhZvalma03s+vMLNPMlpjZRjN7IGZfw8zsNTN7wMw2\nmdmHsb/gzayxmY01s2/MbKWZ3WVmaUW2vc/M1gF3mNkRZvafMK61ZvaYmTUO158ItAZmhb+ybzKz\nLDP7qsjn23cVEf7afcLMJprZZuDK0mIqsp9Dw6uqpjHzOoRxpZvZUWY2N/zca81sSiX8bbqZWX64\nz4Vm1jWc39PMlsSsN8fMFsZMzzOz3iXss8DMss3sE+DjcN794d93s5ktMrPu4fzzgNuAwbFXMvGW\nmSQv/bGkJL8EegM9gJbARiC3yDqdgKOAIcD9wO3AWcCJwCAz61Fk3eVAc+AO4CkzaxIuGw98DxwJ\ndADOBUYU2fZT4CDgzwS//u8O4zoeOBwYCeDulwNfAhe5+wHu/vcSPl/RPlt6A9PdvTEwKY6YCI/3\nDfAG0D9m9tBwX3uBu4AX3L0J0Ar4dwnxFCr1ysbMmgHPAf8CmgH3Ac+FCelN4Ggza2ZmtYCfAC3N\nrIGZ1QNOA0pri+gDZAInhNMLgVOApgRlMt3Marv7CwR/hylhGXcI1x9PHGUmSczdNWjA3QE+B84K\nxz8oHA+nWxL8Z08D2gIFQMuY5euAgTHTTwC/CseHAV8XOdYC4DLgYGAnUDdm2SXAf2K2/aKMuPsC\nbxf3OcLpLOCrUj7rSCAvZlmpMRVz/OHAK+G4ESSk7uH0o8BDQKsyPkNhmW6MGX4TUwbzwvHLgTeL\nbPs6cGU4/ipwMdAFeBGYAvwM6Am8W8rxC4CsMmLcAJwcU2YTK1pmGpJzUBuClKQt8LSZxdZn7yH4\nj19odcz4jmKmG8RMf11k/18AhxJU79QCVpnt+3GcRvClWqhodc/BBFck3YEDwvU3lPWByrAyZrxN\nHDHFegp4wMwOAY4FCtx9frjsZoKrhIVmthH4h7uPKyWO5l56G8KhxcRRWJYAcwkS4MpwfCNwJrAL\nyCtlv/Djcr4JuDrctwONgBYlbFveMpMkpIQgJfkSuMrd3yi6wMzaVmB/rYpMtwFmEHwJ7aL0L8Ki\n1Tt/BvYCJ7n7JjPrCzxQyvrfAfULJ8wsHTiwlGPEE9MPG7pvNLOXgMEE1S2TY5atBq4Jj3s68LKZ\nzXX3z8rabwm+BvoVmdcGmB2OzyWoRvoCuAfYBIwh+PX+P2V9lMIRMzsD+D8EV1FLw3kb+KFKq2gZ\nl6vMJDmpDUFK8iDwZzNrDWBmB5bUIFmK2Prwg8zsBjOrZWYDgeOA5939W+Al4D4zOyBszD6ySPtD\nUQ0JvuS3mFkrgi+uWKsJ6rELLQPqmtkFYd3674E6Je3c3VdVIKZJwJUEbQmT9hWA2UAzOyyc3ETw\nRbo/X5izgWPM7BIzyzCzwQRl+Wy4/HWCq5RMYKG7f0CQMDoTVCfF6wCCK8J1ZlbbzP4vwRVCoW+B\nthZeDlSwzCTJKCFISe4HZgIvmdkWgobTTjHL43mRRuw6C4CjgbUEVSj93X1juOwKoDZBu8UGYDpw\nSMw+ih7rTuBUYDMwC3iyyDr3AL+34G6n37j7ZiCb4JfySmAb/7t6pLhjlBZTcWYSNLCvcvf3YuZ3\nBN40s60EV0Q3uPuKEvZRUpnui8/d1wMXAb8laLe5iaABfUO4fDvwFrDUf7h19HVghbuvKyX+osd+\nIRyWASsIqgBjq3+mh/+uN7NF4Xh5y0ySjLkn5gU5ZlaX4PK1DsFJMsPdbwvvkphK8KtlBTDI3Tcl\nJAhJCmY2jOD5hjOijkVESpawKwR33wn0dPf2BLe/9QzvY74VmOPuxwCvhNMiIhKxhFYZhZevEFwh\npBPc8dCb4FY8wn/7JjIGSQrFVcmISJJJWJURBE+7Am8TNPCNdvebzWyjuxc+8WrAhsJpERGJTkJv\nOw1vP2sfdivwopn1LLLczUy/HEVEkkCVPIfg7pvN7DmCR+dXm9kh7v6tmbUE1hRdX0lCRKRi3L3C\nHTsmrA3BzFoU9lUT9qPyU2Axwe15V4arXQk8U9z2UT/CnSzDHXfcEXkMyTKoLFQWKovSh/2VyCuE\nlsCjYTtCGkG/J6+EPSNOM7PhhLedJjAGERGJU8ISggcP55xazPwNwDmJOq6IiFSMnlROcllZWVGH\nkDRUFj9QWfxAZVF5EnrbaUWZmSdjXCIiyczM8GRsVBYRkepFCUFERAAlBBERCSkhiIgIoIQgIiIh\nJQQREQGUEEREJKSEICIigBKCiIiElBBERARQQhARkZASgoiIACmUEJauWcrWXVujDkNEpNpKmYTw\n+HuP0+2Rbny+8fOoQxERqZZSJiHcfdbdXHPqNXR7pBtzV8yNOhwRkWon5d6HMOfTOVz29GX8qeef\n+PlpP6/kyEREktf+vg8h5RICwLL1y+g1uRfnHXke//jZP8hIS+Sro0VEkoMSQgk27tjIkCeH4O5M\nHTCVpvWaVlJ0IiLJSW9MK0HTek15buhznHDgCXQZ24Vl65dFHZKISFJL2YQAkJGWwb/O+xc3db2J\n7o90Z86nc6IOSUQkaaVslVFRc1fMZfATg/ndGb/jF51+gVmFr6pERJKS2hDK4fONn9Nrci9OP/x0\nHrjgAWqn1670Y4iIREVtCOXQrmk7Xh/+Oqu2reLcieeybvu6qEMSEUkaNSohADSq04inBz9N18O6\n0unhTry/5v2oQxIRSQo1qsqoqInvTuQ3L/2GcX3GcdExFyX8eCIiiaQ2hP305so36T+tP7/u/Gtu\n6naTGptFpNpSQqgEX23+ij5T+nDSQSfx/3r9P+pm1K2yY4uIVBY1KleCwxsfzryr5rFjzw56PtqT\nb7d9G3VIIiJVTgkh1KB2A6YOmMp5R55Hp4c7sXjV4qhDEhGpUqoyKsb0pdPJfj6b0ReOZsAJAyKL\nQ0SkPNSGkCBvr3qbvlP6MuLUEfyhxx/U2CwiSU8JIYFWbV3FxVMvpk2TNozrM476tepHHZKISInU\nqJxALQ9oSd6wPGqn1+aMcWewcsvKqEMSEUkYJYQy1M2oy4S+Exh0wiA6j+nMgpULog5JRCQhEpYQ\nzOxwM/uvmS01s/fN7IZw/kgzW2lmi8PhvETFUFnMjFu638KDFz7IRZMv4rElj0UdkohIpUtYG4KZ\nHQIc4u7vmFlD4C2gLzAI2Oru95WybVK0IRTn/TXv03tybwafOJi7z76bNNNFlogkh6RtQ3D3b939\nnXB8G/Ah0CpcXG1v2TnpoJNYMGIBr698nYunXszWXVujDklEpFJUyc9bM2sLdADeDGf90szeNbOx\nZtakKmKoTAc2OJA5l8/hoPoH0e2Rbny+8fOoQxIR2W8Jv+00rC7KA/7k7s+Y2UHA2nDxXUBLdx9e\nZBu/44479k1nZWWRlZWV0Dgrwt15YOED3DP/HqYOmEqPNj2iDklEapC8vDzy8vL2Td95553J+xyC\nmdUCngVmu/u/ilneFpjl7icXmZ+0bQjFeenTl7jsqcu45+x7GH7q8LI3EBFJgKR9MM2CR3sfBda7\n+40x81u6+6pw/EYg092HFtm2WiUEgI/XfUyvyb248OgL+du5fyMjLSPqkESkhknmhNAdeBVYAhQe\n5HbgEqB9OO9z4Fp3X11k22qXEAA27tjI4CcGk2ZpTBkwhSZ1q13ziIhUY0mbEPZHdU0IAHsK9vDb\nF3/Li5++yMxLZnJM82OiDklEaoikve20pspIy+D+8+/nN11/wxnjzuDlz16OOiQRkbjoCiGB8lbk\nMeSJIfy+x+/JycxRj6kiklCqMkpyn238jN6Te9O9dXceOP8BaqXXijokEUlRSgjVwJZdW7j0qUvZ\n9v02nhj4BM3rN486JBFJQWpDqAYa1WnEM4OfodOhneg0phNL1yyNOiQRkR/RFUIVm/DuBG566SbG\n9RnHhcdcGHU4IpJCVGVUDb3x1Rv0n9afG7vcyE3dblJjs4hUCiWEaurLzV/SZ0ofTjn4FB666CHq\nZNSJOiQRqebUhlBNtW7cmvlXzWfb99vo+WhPVm9bXfZGIiIJpIQQoQa1GzBt4DTOPfJcOo3pxOJV\ni6MOSURqMFUZJYlpS6eR83wOD174IP1P6B91OCJSDakNIYW89c1b9J3al2tOvYbf9/i9GptFpFyU\nEFLMqq2r6Du1L+2atOORPo9Qv1b9qEMSkWoi4QnBzE4EegBtCbqsXgHMc/eEPV1VkxMCwM49Oxkx\ncwQfrfuIGUNm0KpRq7I3EpEaL2F3GZnZ5Wa2EPg7cAjwGUEyaAn83czyzeyyih5YSlY3oy4TL57I\ngBMG0HlMZxZ+vTDqkESkBijttV5NgbPdfWtxC82sETAsEUFJkOlv7X4rx7c4ngsnXcj9593P0JOH\nlr2hiEgFqQ2hGnhv9Xv0ntKboScN5a6z7iLNdLewiPxYwh9MM7O/mVljM6tlZq+Y2Tozu7yiB5Ty\nO/ngk1k4YiHzvpxHv6n92Lqr2Is2EZH9Es9PzXPdfTNwEUEbwpHA/0lkUPJjBzY4kJeveJkW9Vtw\n+iOns2LTiqhDEpEUE09CKGxnuAh4IkwOqs+JQO302jzc62GGdxhO17Fdmf/l/KhDEpEUEk9CmGVm\nHwGnAa+Y2UHAzsSGJSUxM37V5VeM7zOeflP7MfbtsVGHJCIpIq5GZTNrDmxy971m1gA4wN2/TVhQ\nalSOy0frPqL35N5cdMxF3PvTe8lIK+2mMRFJdVXRqJwDFLj73nBWbaBfRQ8olee4Fsfx5og3WbJ6\nCb0m92Lzzs1RhyQi1Vg8VUbXuPvGwolw/JrEhSTl0axeM2ZfOpujmh5Fl7Fd+GT9J1GHJCLVVDwJ\nIc3shxvfzSwdqJW4kKS8aqXX4oELHuDXnX9N93HdeeWzV6IOSUSqoXgSwovAFDM728zOAaYALyQ2\nLKmIaztey9QBU7n0qUvJXZgbdTgiUs3E07ldOkEV0dnhrDnAmJg2hcoPSo3K++XTDZ/Se0pvzmxz\nJvefdz+10nVBJ1ITVEn312ZWH2jt7h9V9EDloYSw/7bs2sLQJ4eyffd2pg+cTvP6zaMOSUQSrCru\nMuoNLCasJjKzDmY2s6IHlKrRqE4jZgyZQcdDO9J5TGc+WPtB1CGJSJKLpw1hJNAZ2Ajg7ouBIxIY\nk1SS9LR07v3pvfyhxx/IGp/F8588H3VIIpLE4kkIu919U5F5BYkIRhLjyvZX8syQZxgxcwT/eP0f\nqDpORIoTT0JYamaXAhlmdrSZPQC8nuC4pJJ1O7wbb454k4lLJnL1zKvZtWdX1CGJSJKJJyH8EjgR\n2AVMBrYAv05kUJIYrRu35rWrX2PLri2cNeEsVm9bHXVIIpJEyvWCnPAW1IZhj6cJo7uMEqvACxiZ\nN5IJ705gxpAZnHLIKVGHJCKVoCruMppsZo3CTu3eAz4ws5srekCJXpql8ceef+Sv5/yVcyaew9Mf\nPh11SCKSBOKpMjrB3bcAfYHZQFugzDemmdnhZvZfM1tqZu+b2Q3h/GZmNsfMlpnZS2bWZH8+gFTc\n4JMGM/vS2dzwwg386dU/qbFZpIaL6wU5ZlaLICHMcvfdxPeCnN3Aje5+ItAFyDGz44FbgTnufgzw\nSjgtEel4aEcWjFjAzI9nMvSpoezYvSPqkEQkIvEkhIcIXp3ZEHjVzNoCZbYhuPu37v5OOL4N+BBo\nBfQGHg1Xe5Qg0UiEDj3gUOYOm0uapdFjfA++3vJ11CGJSARKbFQ2s27AG0Vbd83MgIzwSiG+gwRJ\nZC5wEvCluzeN2deGwumY9dWoHAF35y/z/0Jufi5PD36azFaZUYckIuWQyEblK4C3zWyqmQ0zs0MA\nPFCeZNAQeBL4lbtvjV0Wfuvrmz9JmBm3nXEbuRfkcsGkC5j83uSoQxKRKlTiOxfd/TqAsN7/fGB8\n2AD8H4J+jV4rq8fTsO3hSWCiuz8Tzl5tZoe4+7dm1hJYU9y2I0eO3DeelZVFVlZWvJ9J9lOf4/rQ\nrmk7ek/uzdK1S/ljzz+SZvHULopIVcrLyyMvL6/S9lfe5xDqAz0JEkRXdz+tlHWNoI1gvbvfGDP/\n3nDeX83sVqCJu99aZFtVGSWBNd+tof+0/rSo34KJF0+kYe2GUYckIqVIePfXZtasmNnb3P37Mrbr\nDrwKLOGHaqHbgIXANKA1QWP1oKJ9JSkhJI9de3aR/Vw2i1YtYuaQmbRp0ibqkESkBFWREFYQfHkX\nvle5KfBtOPzc3d+q6MFLOaYSQhJxd+5fcD/3vnYv0wdO5/TWp0cdkogUI+FPKhO8Ie18d2/u7s2B\n84BngRxgdEUPLNWHmfHrLr/mkT6PcPHUixm3eFzUIYlIAsRzhfC+u59UZN577n6ymb3j7u0rPShd\nISStD9d+SO8pvel9TG/u/em9pKelRx2SiISq4gphlZndYmZtzKxt2I/R6rCjO70XoYY5/sDjWTBi\nAe+ufpdek3uxeWdC+zkUkSoUT0IYChwOPAM8TdCecAmQDgxKXGiSrJrVa8bsS2dzRNMj6Dq2K8s3\nLI86JBGpBPFUGbVz98+LzMt09/yEBaUqo2pjdP5oRs4dyeT+kzmr3VlRhyNSo1VFldGTZnZYzAHP\nBNSqKABcn3k9U/pPYeiTQxmdr3sMRKqzeK4QMgnuJroIOBW4B7jI3b9KWFC6Qqh2Pt3wKb0m96Jn\n257867x/USu9VtQhidQ4CX8OITxIN4JeT3cQJINiu5uoLEoI1dPmnZsZ+tRQdu7ZyfSB02lWr7hn\nGkUkURKWEMxsVpFZxwOrgE0E/dL1ruhBywxKCaHa2luwl1tevoUZH89g5pCZHH/g8VGHJFJjJDIh\nZBUz2wEjSAhzK3rQMoNSQqj2xr8znpvn3MyjfR/l/KPPjzockRohkQkhzd1Lfc7AEvTNrYSQGl77\n8jUGTh/ITd1u4sYuNxL0dygiiZLIhDCXoIuKGe6+rMiyYwnedHahu/eo6MFLDEoJIWV8sekL+kzp\nw6ktT2X0haOpk1En6pBEUlYibzs9F1gP5JrZKjNbZmafmNkq4H+A1cA5FT2w1AxtmrRh/tXz2bRz\nE2dPOJs13yX0fgQR2Q/x3mWUDrQIJ9eV9WKc/Q5KVwgpp8ALuOO/dzBxyURmXjKTnxz8k6hDEkk5\nVXLbaXigVgTdVQCsKs9rNMsdlBJCypry/hR+OfuXPNzrYfoe1zfqcERSSiLbEG4Harn7neH0l8Bm\noDYw3t3vqehBywxKCSGl5X+dz8VTLyY7M5vbut+mxmaRSpLIhLAYOMPdtxVOu3uHsProVXdP2FtS\nlBBS39dbvqbv1L4c3exoxvYeS71a9aIOSaTaS2hfRoXJIHR/OG8voP+9sl9aNWrFq8NexXHOHH8m\n32z9JuqQRGq80hJCAzOrXTjh7uMBzKwOcECC45IaoF6tekzqN4m+x/Wl85jOLPpmUdQhidRopSWE\nJ4AHzaxB4Qwza0jQp9ETiQ5MagYz4/Yzbuff5/2b8x8/nynvT4k6JJEaq7Q2hAzgT8AI4Mtwdmtg\nLPA7d9+TsKDUhlAjvfvtu/SZ0ofLf3I5d/a8kzSLp3d2ESmU8NtOzaw+cFQ4udzdt1f0YHEHpYRQ\nY635bg39pvbj4IYHM6HvBBrUblD2RiICVMELctx9u7svCYeEJwOp2Q5qcBCvXPEKjes05vRHTueL\nTV9EHZJIjaFrckk6dTLqMLb3WK445Qq6ju3K61+9HnVIIjVC3E8qVyVVGUmh2Z/M5spnruTen97L\nsPbDog5HJKlVRRtCGnAp0M7d/2hmrYFD3H1hRQ9aZlBKCBLjw7Uf0mtyL/oe15e/nvNX0tPSy95I\npAaqioTwIFAAnOXux5lZM+Ald+9Y0YOWGZQSghSxfvt6Bj0xiLoZdZnUbxKN6zaOOiSRpJPwRmWg\ns7tnE7xPGXffAOgN6lKlmtdvzguXvkCbxm3oOrYrn274NOqQRFJOPAnh+7D/IgDM7ECCKwaRKlUr\nvRajLhzFLzr9gtMfOZ3/fv7fqEMSSSnxJIQHgKeBg8zsz8BrQMJ6OhUpS3ZmNpP6T2LIk0N4cNGD\nUYcjkjLifUHO8cDZ4eQr7v5hQoNSG4LEYfmG5fSa3Iuz253NP3/2T2qlqyZTaraqaFSe6O6XlzWv\nMikhSLw279zMkCeHsHvvbqYNnEazes2iDkkkMlXRqHxSkQNmAKdV9IAilalx3cY8e8mztD+kPZ3H\ndOajdR9FHZJItVViQjCz281sK3CymW0tHIA1wMwqi1CkDOlp6fz93L9ze/fb6TGuBy8sfyHqkESq\npXiqjP7i7rdWUTyFx1SVkVTI/C/nM3D6QG45/RZ+1flXej2n1CgJb0MID9IUOBqoWzjP3V+t6EHj\nOJ4SglTYF5u+oPeU3mQemsmoC0dRO7122RuJpICEtyGY2c+BV4GXgDuBF4GRcQb3iJmtNrP3YuaN\nNLOVZrY4HM6rWOgixWvTpA2vXf0a63es55wJ57D2u7VRhyRSLcTTqPwroBOwwt17Ah2AzXHufxxQ\n9AvfgfvcvUM4qMJXKl3D2g15ctCT9GjTg05jOrFk9ZKoQxJJevEkhJ3uvgPAzOq6+0fAsfHs3N3n\nARuLWaSKXUm4NEvjT2f9ibvPupuzJ5zNjI9mRB2SSFLLiGOdlWEbwjPAHDPbCKzYz+P+0syuABYB\nv3X3Tfu5P5ESDT15KEc1O4p+U/vxwdoPuLX7rWpsFilGud6HYGZZQCPgBXf/Ps5t2gKz3P3kcPog\noLBS9y6gpbsPL7KNGpWl0n295Wv6Tu3Lsc2PZUzvMdTNqFv2RiLVyP42Kpd6hRA+hPa+ux8H4O55\nFT1QIXdfE7P/McCs4tYbOXLkvvGsrCyysrL299BSw7Vq1Iq5w+Zy9YyrOXP8mTwz+BlaHtAy6rBE\nKiwvL4+8vLxK2188zyHMAG5w9wq93LaYK4SW7r4qHL8RyHT3oUW20RWCJIy7c/e8u3norYd4ZvAz\nnHaoHryX1FAVfRnNI7izaCHwXTjb3b13HMFNBs4EWgCrgTuALKA9wd1GnwPXuvvqItspIUjCPfXh\nU1z77LXkXpDLoBMHRR2OyH6rioSQVcxsd/e5FT1oWZQQpKq88+079J3SlytOuYKRWSNJs3huvBNJ\nTlXypHJVU0KQqrR622r6TevHIQ0PYULfCTSo3SDqkEQqpCp6OxVJaQc3PJj/XPEfGtVpRPdx3fly\n85dRhyQSCSUEEaBORh0e6f0Il518GV3GdOGNr96IOiSRKhdPX0bnmFm9qghGJEpmxm+7/ZYxvcfQ\nZ0ofHn3n0ahDEqlS8TQqTwC6EHRB8Wo4zHf34rqkqJyg1IYgEftg7Qf0mtyL/sf3556z7yE9LT3q\nkETKVGWNymZ2KDAAuAk41N3j6faiYkEpIUgSWL99PQOmD6BBrQZM6j+JRnUaRR2SSKmqovvry83s\nIeBJ4Bzgf4AeFT2gSHXRvH5zXrrsJQ5vdDhdx3bls42fRR2SSELFU2W0HvgUGA3kufvnCQ9KVwiS\nZEblj+KPc//IlAFTyGqbFXU4IsWqigfTDDgROCMcjgKWuftlFT1omUEpIUgSevmzl7n0qUu5q+dd\nXHPaNVGHI/IjVfEcwgFAa6AN0BZoAhRU9IAi1dU5R5zD/Kvm8883/8kNs29gT8GeqEMSqVTxXCEs\nAV4D5gGvuvvKhAelKwRJYpt2bmLIE0PY63uZNmAaTes1jTokEaBq7zI6gKAPo20VPVi8lBAk2e0p\n2MPNc27m2WXPMuuSWRzbIq6XCIokVFXcZXSymS0GlgIfmNlbZnZSRQ8okgoy0jK472f3ccvpt3DG\nuDN4cfmLUYckst/iqTJ6A7jd3f8bTmcBf3b3bgkLSlcIUo3M+2Ieg54YxM3dbuaGzjfoITaJTFXc\nZfSuu59S1rzKpIQg1c2KTSu47KnLWLllJdd1vI7hHYZzYIMDow5LapiquMvoczP7g5m1NbN2ZvZ7\nQE/oiMRo26Qt86+ezxODnmDZ+mUc8z/HcMXTV7Bg5QL040aqi3iuEJoBdwKnh7PmASPVl5FIydZv\nX8+4d8YxetFomtRtQk5mDkNOGkL9WvWjDk1SmF6QI5LECryAF5e/SG5+Lm+ufJMrT7mS6zOv56hm\nR0UdmqSghCUEM5tVynZxvVO5opQQJBV9tvEzHlr0EOPeGcdph55GTmYO5x91vhqhpdIkMiFklbKd\n3qksUkE7du9g2tJp5Obnsnb7Wq477TqGnzqcFvVbRB2aVHOJTAht3P2LCke2H5QQpKbI/zqfUYtG\n8cxHz9D72N5kd8ymU6tOBF2IiZRPIhPCYnfvEI4/6e79K3qQcgelhCA1zPrt63lk8SOMXjSaZvWa\n7WuErldLLyuU+FVVQtg3XhWUEKSmKvACXlj+Arn5uSz8emHQCN3xeo5sdmTUoUk1UBXPIYhIFUmz\nNC44+gKeG/ocC0YsIM3S6DK2Cxc8fgHPLXuOvQV7ow5RUlhpVwh7ge3hZD1gR8xid/eEvU9QVwgi\nP9ixewdTl04lNz+X9dvX73sSunn95lGHJklGzyGI1CD5X+eTm5/LjI9n0OfYPuRk5pDZKjPqsCRJ\nKCGI1EDrtq/b1wjdon4LcjJzGHziYDVC13BKCCI12N6CvfsaofO/yeeq9ldxXcfrOKLpEVGHJhFQ\nQhARAD4o+f1fAAAPPElEQVTd8CmjF41m/Dvj6XxYZ3IyczjvqPNIM907UlMoIYjI/7Jj9w6mvD+F\n3PxcNuzYwPUdr+fqDlerEboGUEIQkWK5O/nfBI3QMz+eSd/j+pLdMVuN0ClMCUFEyrRu+zrGvj2W\n0YtGc1CDg4JG6JMGUzejbtShSSVSQhCRuO0t2Mvs5bPJzc9l0TeLuKr9VVzf8XraNW0XdWhSCZQQ\nRKRClm9YzoOLHmT8O+PpclgXcjJz+NlRP1MjdDWmhCAi+2X77u37GqE37dy0rxG6Wb1mUYcm5aSE\nICKVwt1Z+PVCcvNzmbVsFhcfdzHZmdl0PLRj1KFJnJI6IZjZI8CFwBp3Pzmc1wyYCrQBVgCD3H1T\nke2UEEQitPa7tYxdPJYHFz3IwQ0PJiczh0EnDlIjdJJL9oRwBrANmBCTEO4F1rn7vWZ2C9DU3W8t\nsp0SgkgS2Fuwl+c/eZ7c/FzeXvX2vieh1QidnJI6IQCYWVtgVkxC+Ag4091Xm9khQJ67H1dkGyUE\nkSSzfMNyRueP5tF3H6Xr4V3J7pitRugkUx0TwkZ3bxqOG7ChcDpmGyUEkSS1ffd2Jr83mdz8XLbs\n2sL1Ha/nqg5XqRE6CexvQsiozGDKy93dzIr95h85cuS+8aysLLKysqooKhEpTf1a9Rl+6nCu7nA1\nC75eQG5+Lkf++0j6HdePnE45nNry1KhDrDHy8vLIy8urtP1FVWWU5e7fmllL4L+qMhKp3tZ8t2Zf\nd9wtG7YkJzOHgScOVCN0FauOVUb3Auvd/a9mdivQRI3KIqlhb8FenvvkOXLzc1m8ajFXd7ia6zpe\nR9smbaMOrUZI6oRgZpOBM4EWwGrg/wIzgGlAa3TbqUjK+mT9J4xeNJoJ706g2+HdyMnM4adH/lSN\n0AmU1AmhopQQRFLH9t3bmfTeJHLzc9n2/bagEbr9VTSt17TsjaVclBBEpFpwd95c+Sa5+bk898lz\n9D++P9mZ2WqErkRKCCJS7az5bg1j3h7Dg4sepFWjVkEj9AkDqZNRJ+rQqjUlBBGptvYU7OG5ZUEj\n9Lur3+Xq9kEjdJsmbaIOrVpSQhCRlLBs/TJG549mwpIJdG/dnZzMHM454hw1QpeDEoKIpJTvvv9u\nXyP09t3bub7j9QxrP0yN0HFQQhCRlOTuvLHyDXLzc3n+k+fpf3x/cjJz6NCyQ9ShJS0lBBFJeau3\nrd7XHfdhjQ4jJzOHAScMUCN0EUoIIlJj7CnYw7PLniU3P5clq5cwvMNwrut4Ha0bt446tKSghCAi\nNdLH6z5m9KLRTFwykTNan0F2ZnaNb4RWQhCRGu2777/j8fceJzc/lx27d5Cdmc2w9sNoUrdJ1KFV\nOSUEERGCRujXv3qd3PxcZi+fzcATBpKdmU37Q9pHHVqVUUIQESli9bbVwZPQbz1I68atycnMof/x\n/VO+EVoJQUSkBHsK9jDr41mMWjSK91a/x/AOw7m247Up2withCAiEoeP1n3E6PzRPPbeY/Ro04Oc\nzBzObnc2wZt8U4MSgkgFucPu3bB9O+zYUfZQ1no7d0JaGmRkFD/UqlXysqoYijt+Wg28IWfb99v2\nPQm9c89Osjtmc2X7K1OiEVoJQVKGO3z/fXxf0JX1JZ6eDvXqlT3Ur1/2OnXrBp9hz57ih927S16W\n6KG4Y+/eDWbRJaPok6Hz+srXyM3P5YXlLzDohEFkZ2ZzyiGnRP1focKUECRh3INfvZXxxRvvL+yM\njP3/Yi7PehkZUZdytAoKqj4RRZkEY4e9e4MfBBkZkNb4WwpOGcPunzxE+rY2NPwwhwO+6k+ttNqR\nJcMePeDoo8v391RCqEEKCuL7gq7MKpDatavmi7l+/eAXdnp61KUsNYV7kBRik8TO7/fw3PKZjF86\nio82vM/AI0cw8IhrObD24VWeCK+6Crp3L99nUkKIUEFB5f16jme9XbugTp3EfjnHrlO3bs2sYxaB\noBF6VP4oHn/vcc5scybZmdlJ3withBBj796qq97YsSOo765bN7G/mmOn69TRF7RIVdv2/TYeXxI8\nCf393u/JzszmylOupHHdxlGH9iNKCKErroDHHqv4r+GKrFenTtAoJyKpz92Z/+V8Ri0ata8ROqdT\nDj85+CdRh7aPEkJoz56g/llf0CKSaN9u+5aH33qYh956iHZN25GTmUO/4/tRO712pHEpIYiIRGT3\n3t3M/Hgmufm5fLjuQ0Z0GMG1Ha/lsEaHRRKPEoKISBL4cO2H+xqhe7brSXbHbM5qd1aVNkIrIYiI\nJJGtu7by2JLHyM3PZa/vJbtjNlecckWVNEIrIYiIJCF3Z96X88jNz2XOp3MYfOJgsjOzOfngkxN2\nTCUEEZEkt2rrKh5+O2iEPrLpkeRk5nDx8RdXeiO0EoKISDWxe+9uZnw8g1H5o/hw3Yf8/NSfc+1p\n19KqUatK2b8SgohINfTB2g8YlT+KSe9N4qx2Z5GTmUNW26z9aoRWQhARqcZiG6ELvIDszKARulGd\nRuXelxKCiEgKcHde/eLVoBH6szk8dvFjXHjMheXahxKCiEiK+WbrN9TNqEuzes3KtZ0SgoiIAPuf\nENR3poiIAEoIIiISiuwFgma2AtgC7AV2u3unqGIREZForxAcyHL3DkoGJcvLy4s6hKShsviByuIH\nKovKE3WVkd5eUAad7D9QWfxAZfEDlUXlifoK4WUzW2RmP48wDhERIcI2BOB0d19lZgcCc8zsI3ef\nF2E8IiI1WlI8h2BmdwDb3P0f4XT0QYmIVEP78xxCJFcIZlYfSHf3rWbWADgXuLNw+f58IBERqZio\nqowOBp4Oe/XLAB5395ciikVEREiSKiMREYleJHcZmVldM1tgZu+Y2Qdmdk84/29m9qGZvWtmT5lZ\n45htbjOzT8zsIzM7N4q4E6G8ZWFmbc1sh5ktDodR0X6CylNKWdwVlsM7ZvaKmR0es01NOy+KLYtU\nPS9KKoeY5b81swIzaxYzLyXPCSh/eZT7vHD3SAagfvhvBvAm0B34KZAWzv8L8Jdw/ATgHaAW0BZY\nXrheKgzlLIu2wHtRx1zFZXFAzPJfAmNq8HlRUlmk7HlRXDmE04cDLwCfA81qwjlRgfIo13kR2XMI\n7r49HK0NpAMb3H2OuxeE8xcAh4XjfYDJ7r7b3VcQ/JFT5unmcpZFSiuhLLbGrNIQWBeO18TzoqSy\nSFnFlUM4fR9wc5HVU/qcgHKXR7lElhDMLM3M3gFWA/919w+KrHI18Hw4fiiwMmbZSqByXkKaBMpZ\nFgDtwsu/PDPrXmWBVoGSysLM7jazL4FhQOFlco08L2LK4kqCq8dCKXleFFcOZtYHWOnuS4qsntLn\nBJS7PKAc50WUVwgF7t6e4JdvDzPLKlxmZr8Dvnf3SaXtIsEhVplylsU3wOHu3gH4DTDJzA6o6pgT\npaSycPffuXtrYBzwr9J2kfgoq0YcZTEe+Ge4esqeF8WUwwXAbcAdMauVdqt6ypwTUO7yKNd5EXVf\nRrj7ZuA5oCOAmQ0DLgAujVnta4L6sUKHhfNSSjxl4e7fu/vGcPxt4FPg6CoPNsGKlkWMSUBmOF4j\nz4sY+8qiJpwXMeVwKtAOeNfMPif4u79lZgdTQ84JiKs8Dir3eRFRo0gLoEk4Xg94FTgbOA9YCrQo\nsn5hQ1Ht8IN/SnjLbHUfKlAWLQge6gM4guCSuEnUnyPBZXFUzDq/BCbW4POipLJIyfOipHIosk5x\njcopd05UsDzKdV5E9WBaS+BRM0sjuEqZ6O6vmNknBH/IORY8tPaGu2d7UEc2DfgA2ANke/gJU0C5\nygI4E7jTzHYDBcC17r4potgrW0ll8YSZHUvw7oxPgesBauh5UWxZAD2AP6bgeVFsORRZZ9/fPMXP\nCShneVDO80IPpomICJAEbQgiIpIclBBERARQQhARkZASgoiIAEoIIiISUkIQERFACUEiZGbNY7rl\nXWVmK8Pxt80syvd9/4iZnWlmXROw39LK4Ggze6+Ubaea2RHh+IqYLo9PM7PPzKy9mfU2sz9UdtyS\nmpLqP53ULO6+HugA+96rvdXd74sqHjNLd/e9JSzuCWwF3ijH/jLcfU9p65RWBmbWtpR9HwU0cPfP\nCncVzv8JMB0Y5O7vmNm7BA8m/cXdd8cbu9RMukKQZGLhr9s8M1tkZi+Y2SHhgjwzu8/M8i14cVCm\nmT1tZsvM7K5wnbYWvBTlMQteHjLdzOqFy0rb7z/NLB/4lZldZGZvhr/Q55jZQeEX87XAjeH87mY2\n3sz6xwS+Lfw3y8zmmdkM4H0Leqb8m5kttODFNteUVQZFpjOK+zzAEGBmkXVPBJ4GLnP3RQDhU7pv\nELy3XKRUSgiSTAz4NzDA3TsS9Gx6d7jMgV3ungmMBmYA1wEnAcPMrGm43jFArrufAGwBssPqpweA\n/iXst5a7Z4a/zOe7exd3PxWYCtzsQb/6DwL3ufup7j6fH/egGTvdAbjB3Y8DRgCb3L0TQb/8Py/t\nl38xji36ecL5pwOLipTdM0COu79eZB8LCbowECmVqowkmdQh+IIv7L8pnaD73kKFv4jfB95399UA\nZvYZQQ+XW4Cv3L2wWucx4AaCt0idCLxcwn6nxowfHvaFcwhBX1KfxSwrrYvlWAvd/Ytw/FzgZDMb\nEE43Ao4CVsS5r+I+zz+ANsCqmPUcmEOQcF7yH16uBMFnPS/O40kNpoQgycSApe7erYTlu8J/C2LG\nC6cLz+XYX+oWTpe13+9ixh8A/u7uz5rZmcDIErbZQ3iFHXY0VruE/QH8wt3nlLCfshT3eWKn/9dx\ngIeAUQRXT4XSSLF3AkhiqMpIksku4EAz6wJgZrXM7IRy7qN14fbAUGAe8HEZ+439Ym3ED1cPw2Lm\nbwViXyyyAjgtHO9N8A7f4rzID9VWmNkxZlZ/Pz8PwBcEPV/GKgjXOc7M7oyZ3zJcX6RUSgiSTPYC\nA4C/WvCKwMVAcbd6OiX/4v0YyDGzD4DGwOjw7prS9hu7r5HAdDNbBKyNWTYLuDi8JfR04GHgzHB/\nXYBtJexvDEFXzG+Ht5COpvQrcy8y/qPPEy6bz/9+YY4DuPsuggTV28wKu8buRNBvvkip1P21pIyw\nsXaWu58ccSgJFz5/8IC7X1jGemnA20DHsm6BFdEVgqSaGvELJ3z+YKuZHVnGqhcBTygZSDx0hSAi\nIoCuEEREJKSEICIigBKCiIiElBBERARQQhARkZASgoiIAPD/AbOr0W3vHaQlAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7f0dc1cf0f90>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      " D = 4.54 m\n"
     ]
    }
   ],
   "source": [
    "from math import sqrt\n",
    "%matplotlib inline\n",
    "from matplotlib.pyplot import plot,title,show,xlabel,ylabel\n",
    "\n",
    "H = 0.036#               #Humidity is in kg/kg at 811 K\n",
    "#Taking R as 90 per cent and P as 101.3 kN/m2, then, for assumed values of Tb of 321, 333 and 344 K\n",
    "#Pw = 13,20 and 32 kN/m2, respectively\n",
    "#G = 27.8, 12.9 and 6.02 kg/s, respectively.\n",
    "     \n",
    "#for Tb = 321, 333 and 344 K,\n",
    "#G = 7.16, 7.8 and 7.54 kg/s respectively.\n",
    "Tb = [321, 333, 344]\n",
    "G1 = [27.8 ,12.9, 6.02]#           #Temperature is in kelvins\n",
    "G = [7.16 ,7.8, 7.54]#             #flow rate in kg/secs\n",
    "plot(Tb,G)\n",
    "plot(Tb,G1)\n",
    "title(\"Temperature vs Flow rate\")\n",
    "xlabel(\"Temperature Tb(K)\")\n",
    "ylabel(\"Flow rate G(kg/secs)\")\n",
    "show()\n",
    "\n",
    "\n",
    "#Plotting G against Tb for each equation on the same axis, then\n",
    "Go = 8.3#                         #Gas flow rate is in kg/secs\n",
    "Tb = 340#                         #temperature is in Kelvins\n",
    "uf = 0.61#                        #velocity is in m/secs\n",
    "\n",
    "D = sqrt(340*(8.3+(1.58*1.26))/(278*0.61))\n",
    "print\"\\n D = %.2f m\"%(D)"
   ]
  }
 ],
 "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
}