{ "metadata": { "name": "" }, "nbformat": 3, "nbformat_minor": 0, "worksheets": [ { "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "Basic test of the wflow BMI interface" ] }, { "cell_type": "code", "collapsed": false, "input": [ "import wflow.wflow_bmi as bmi\n", "reload(bmi)" ], "language": "python", "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "pyout", "prompt_number": 1, "text": [ "" ] } ], "prompt_number": 1 }, { "cell_type": "code", "collapsed": false, "input": [ "bmimodel = bmi.wflowbmi()\n", "bmimodel.initialize(configfile='../../examples/wflow_rhine_sbm/wflow_sbm.ini')" ], "language": "python", "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "pyout", "prompt_number": 2, "text": [ "0" ] } ], "prompt_number": 2 }, { "cell_type": "code", "collapsed": false, "input": [ "bmimodel.get_var_name(1)\n" ], "language": "python", "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "pyout", "prompt_number": 3, "text": [ "'SurfaceRunoff'" ] } ], "prompt_number": 3 }, { "cell_type": "code", "collapsed": false, "input": [ "bmimodel.get_var_shape('SurfaceRunoff')" ], "language": "python", "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "pyout", "prompt_number": 4, "text": [ "(169, 187)" ] } ], "prompt_number": 4 }, { "cell_type": "code", "collapsed": false, "input": [ "imshow(bmimodel.get_var('SurfaceRunoff'))" ], "language": "python", "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "pyout", "prompt_number": 5, "text": [ "" ] }, { "metadata": {}, "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAARoAAAD/CAYAAAAuTRk/AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzsnXd8VFX6/993SnpCQhKSkEIoEaQjSBUIVRAXRVFEVAQV\nu+uyq+66XxVdC7rruivq/iwoVkTpIhYQIghSlRo6gYQEUkjvmfL74+Zm7kwyk0mmJ/f9es0L5pZz\nz9zc87nPec5zniMYjUYjCgoKCi5E5ekKKCgotH0UoVFQUHA5itAoKCi4HEVoFBQUXI4iNAoKCi5H\nERoFBQWX43Sh+f777+nVqxcpKSm8+uqrzi5eQUHBBxGcGUej1+vp2bMnmzdvJj4+nquvvprly5dz\n5ZVXOusSCgoKPohTLZo9e/bQo0cPkpOT0Wq13Hbbbaxbt86Zl1BQUPBBnCo02dnZJCYmNnxPSEgg\nOzvbmZdQUFDwQTTOLEwQhGaPGThwIAcPHnTmZRUUFLyAAQMGcODAgSb3OVVo4uPjycrKavielZVF\nQkKC2TEHDx7kueeeY9GiRc2Wd/hwOqmpX1NY2Do30ujRAaSlPYlK5Z2Da4sWLTK7D8eOneTXX9Mb\nvvfp04VhwwZ5oGbuxfI+tFd8/T7YMjScKjRDhgzh1KlTnDt3js6dO7NixQqWL1/e6vL8/LQkJ2uI\niDBQUmKgoKDtzv/Mzc3jq692smjRuYZtc+eeJyoq3Oo54eFhREZGuqF2raesrIyiomI6d45Do3Hq\n46bgQzj1L6/RaHjrrbe49tpr0ev13HPPPQ6NOHXrlsw338zHYDCwdOkms0bY1li06Gu+/DLfbNuq\nVUX89NOnVs958MGuPP30na6umkNs3PgLS5Yc4Ouv5xEXF+vp6ih4CKe/YqZOncrUqVNtHpOammpX\nWVqtls6d4wBI6lDNTaxu2HeKFA7Tz+b5J07UcOedb1o16fr1i+Spp+6wqy7O5u23v2bHjgLuuOO/\nAGzfXkJxsfkx5eVQXm6wWkZpaY0rq+gUKipqOHiwhocf/pygIA3jxnXhnntuNDvG3uehrdOW74NH\nbNnW3NAr+qUwbsEIAH7/4ThR5wsYzH7bJ+XBiS9MX8sI5SRXNHzv3buUiAiTeCUlRTFlypgW1601\n7NyZzebNUUCRW67nCdLSdrFtWybl5bBmTRkAmZlV6PXiPZ8wYRDdu3dt0w2sJbTl++AzneZREyYw\nasIEAB566B1Uq7YyPC+tRWXk0YlCOpq+p8PT9+dRQgd0aBgzJoDk5Dizc1QqgS5dEvH393f4Nzib\ny5erOX8+k6SkRLtG/NxFXV0dmZkXePPNHQ0CI7F9ezXbtx8C4N//rmHqVB3+/n4kJSWgVqs9UV0F\nN+DUyGC7LigIOHrJ4uJi3vzPKl57/niLzosnm+msN9umQ8MKZnGROPz9BcItfK9hYSpWr55J377O\njW6eM+c/fPFFcfMH2iA4WGDs2CDWrPkjfn5+TqqZ42Rn5zBjxjKOHKmjqsr637pDBwgIEOjZ0581\naxbQsWNHq8cqeD+22rZPCg3AiROnOHjwTJP7XnppH4cO6RptD6aCZM6Z1wcjXThPIFVcIpZfuMZs\nv1YrMG5cII8/PoypU8e2uJ75+fk8//xK8vOrzLbv2lVOZqZ1/4u9REYKpKaG8uyzU+jfv7fD5TnK\ntm27+ec/t7F1ayUVFfb9ncPCYNy4UP72t3EMG3aVi2uo4CrapNDY4q23vub48YKG72vX5pOd3XSj\nVmFgOLsIp5ha/CihQ8O+YsI5RQoA//nPlfzxj7OsXjMtbRdHj2Y12l5YWMl//3uOy5dd+5t/+GEC\nkyePdqiMnTv38fvvGQ6VsXfvRT7+uLBV565cOZKbb57s0PUVPIettu0zPpqW8Mgjt5h9r6x8k717\nSxsdV15u4Nw52MlIAJLIZBrfNuy/SBzFhFNEBJculZGZmUVSUqJZGTqdjoyM87zzzk6+/rrxNdzF\nuXP5HDlyjKCgQLp27cK5c5lUVFS2qIz339/JsmWtEwlnkJVVyMWLl5Rh8DZIm7RoLCktLaWurq7R\n9q1b93HLLb82fFejx4/ahu9dyWASm1jObCqCY5g5M4Jlyx4zKyM/P5/p09/nwIFaqqtd9xuaIyxM\n7OYNGxbE+vULmTnzP2zfXt6iMiorselTcTWhoQIPPxzPK6/c67E6KLSedtd1spe8vDx++UWcd/XZ\nZ4cajZCEUE4XztOdMxzjSio6X8Hw4SFmx1RXG9i+vZKyMu/4TZGRMGZMKL/8Uk5+vnfUyV6eeCKe\nOXNGMmBAH09XRaEVKEJjB998s4UtW05RV2dg5co8cnPFOmrQMZQ9aNBRix96tBzlSioJ8nCN2x6K\nj8a3aXc+mtbwhz+M5w9/GE9VVRXFxe+xaVMxeXlGdGjYyUhGspOBHESPH6WEcoEEKhSxUVCwC++c\n1uxBAgICePfde7n11k5m2/cwlI+4m6+5ifH8zJW0LIZHQaE9o3SdrJCefoLly3fw4ouZZtvV6Eki\nixRO04l8qgngRyZQSpiHaur7JCWpePnlwYwfP1gZcfJhlK5TK+jduyezZ0N1tRj4t2VLPr/9Voce\nNRkkE0g1KowYERjIYU7Sg0vEeLjWvsfAgVrmzOnKjTemEhwc7OnqKLgIxaJpBr1ez/Hjp3jxxe/4\n8suSRvtDqOB+lvELI9nNYMDxaF/XIZ8P5bm/QUAA9OqlRasVmDmzC08+OcdjdVFwHopF4wAVFRXc\ne+9q9u2rbXK/CggD/AkAAoEKN9aupQiANHFRh6fEJjFRxerVdxATE60kw2onKM7gZjAajVRXG9E1\nnjoFQCUBfMYfiOUi17ENCK3/aN1YS+qv52flo8H0TtHVf1R44s8/dWowS5dOoXPnWIKCgrxqMqiC\n61BeJw6iQ80pkhjJISIoQ2zEYGrcBqCWxtaDHybrgvpjBKwLVC2gt9hmKRZNWShG2XbL/0tlgCu7\nfOHhMGdODGFhflx9dRKjRw912bUUvBNFaJyKHpBmaQciioYRsRHrEIVEathyK4P6Y2wJjbRfwlhf\nlrr+/3W0TCwMNC9UzulaRUaq+etfp5OQEO+U8hR8D0VoXEZV/UeF2JWqQhQKaWSlHJBPegxBbPzW\nJmaGyM4FUdTqEC0dK/26ZjFgEicN5kImXcObndsKvoLio2mG4OBg/t//m8qMGaFN76ea+1hHHh3Z\nYJHLRsSI6CDWIFofZfUfS3GoAmzlAK6SnVtWX2YNohBoZJ/WZqnTY/Lf6HGWNTNlShCffDKV6Ogo\np5Sn4JsoQtMMGo2GYcOuIimp6RgPHSrOE0sAdcTRVIoFI6YGrEO0Qprq5uhp7IOx3F8n+8gFwYBJ\nGFqb0tMo+0iWjmNiM316CAsWDGLkyCFemQpVwX0oXSc76d69I2PGFFNXZ+TQodqG7HF61JwikREc\nIYE89tHLSglND487jtzBq0IUGmfEy7S+y+TvDwMG+PHgg0MZOLAn27btNtsfFBRA//5XKiNO7Qgl\nYK+FXL58mdTU/3HkiNj16UAFT7OMHxjOFgZ7uHYgioz0/jBg20pyDUlJKrZuvYuuXbvw0UfruOce\n82VSe/ZUk5a2gNhYJZK6LaEE7LmQCgJ4lxmM5Ag3sY3VuGe5FutIXTUwBei5V2xycw3cfvsXBASo\nuHSpccIxhfaHIjQOokPNWToznt+IoGUZ7VyHZy3GmhrYvbv5ruLXX//AkSM5GI1GBEFg1qyR9O7d\n0w01VHA3itA4iRyi6EQx/RCTe18giiKaHqlyH3L/jfdQWWnkl18O8M47v5OWZsp/2r9/vCI0bRRF\naFqA0WhEp9PRVDd0I8O5hiP8H58B8B7T2MogDK0eBWq7ZGUZzHI1K7R9lOHtFrBt226mT19KRkbT\nAXKH6MZz3M1z3E0i+dzEdjfXUEHBO1EsGjvZsGEry5b9zp491n0PpQRRWp/eM4hqepDDDHayiaso\nJ8BdVfVZvvjiMP7+Wq6/fpynq6LgZBShaQadTseePQf44IP9rFtnv7P3N1LQo+Za9pNHB46RRKHH\nfTbezerVZej1vxEREcKwYYOUFBJtCCWOphlKSkpITX2LAwdaN0wbTQmL+IzVjOJn+qFr9RSB9kO/\nfhp+/PE+oqOjUKuV++Ur2Grbio/GxRQTwr+YyRVcZCY7PV0dn+DMGR3Tpn3Ajh37PF0VBSehCI0N\n0tNP8Oyzy8nJae3saKhDTQYx7KMHKozMYgfBNidPKlRWwm+/1VFU5C1xSQqOogiNDU6fzubNNy+R\nl+d4V28vKRQSylR+I9DueU+W85YUFHwTxdtmA41GRWio2NBraozUumpepFWkRFmtt6gUFLwBxaKx\nwahRg0hLm0Fa2gxuvjnc4fLkqcGbR11/hvsnRSooOBtFaGzQoUMHrrqqP1dd1Z9OnRyPgzlGAhsY\nzC38yhVctHGklJ7T8ZwwCgregCI0dtKvXyzDhjmWP6UaLcUE05FyApr008hzCnvnPCV3snv3OQ4f\nTvd0NRScgCI0dnLPPTfy5JNXO1RGfzKZxS4+YwzHSLDYK4mMipal0hRoq07jV17JYtkyZRpHW0Bx\nBruRX0mhkBD+wga+YBS/kiLbG4goFi0Z+lZj/ieU8v7KseUZkicnV1BwHYrQuJEyAjlOZ75jIAPI\nYhSnAdGGWc/VHCcekxPYXiydxfLlWuQ5gJvC+7tmW7YU8Oij/wPggQfG0aePtVSpCt6MIjRuppwA\nNjKQWewmicuEUMsIznGQJI4TQ8uExjKhuZrGvWFbQuP9HDhQx4EDuQBMmXJRERofRREaD7GCYaxg\nGMkU8i5fEoSeYKqowBGHc3MrKSgoeAbFGexhsunA/dzGSM5xL7s8XR0FBZegCI2HqUPNOToSQi3R\nXpNz2Dt5++19rF69ydPVUGgFitB4CXtJpAx/xnMKP6X70yTffVfB7t3nPF0NhVag+Gi8hA8Zxi0c\n5CF2cIFwMomgWvnzNKKsrI78/HyioqIQBPfHDhUXF1NRUWm2LSgokIiICLfXxZdQnmQv4nt6cYlQ\n3mYNLzGBNLp7ukpex5df5pOX9ymff/6wR5bZfeONNXzxxQWzbTfcEMO//rXA7XXxJRShsZPlyzey\ndOlhl16jDH8O0Zn/xwhGc5YIqlhDX5de09coKjJy4UK1x7I05udXcfq0edc2N7fKI3XxJRShscK+\nfQdJT89s+P7558f46SfXP1AlBLCKfkRQSU/ymMRJ0uhOnZIC1KNUV1fz4487OXGirNG+s2cr+eST\nbwAYNaov3bt3deq1jx49zv79pwCYOHEInTvHObV8d6AIjYyqqipyc/MBePfd7XzwQYHH6vIBw7iN\nA9zHbjKJIJNwqsyiftsvNTVGzp/PIikpgcDAQLdcs6ysjL/9bQfp6Y0d9Tt31rBz534A/vnPMmbO\nVKPRqImNjWmUYD0/Px+VSkVkZKRd183Ly+Pzz3fwyitZACxbZmTs2IEN+2NjO1FdXU1xcWnDtvDw\nMMLDHU9r4kyU5OQytm79lXnzxOHToiIDpaXNnOBiQqhhEDm8xHf8H1PYRjfPVshL8PcXiIkR+PTT\naxkzZphbrpmfn09q6v9rUmjkREZCSIiKpCQ/Vqy4m7i4WLP9jz76Pzp08OfFF+fbdd35899k7doi\niorENtOpk0BgoMkJ/tlnU/n99wxef/14w7aFC3vy2GOz7P1pTsNW21aERsbGjWlMm5bm6WqYEU4V\nqZzhKiGHSME02vGOcQRHjTEerJnnmTgxkM6drecJevrp6+jZM8Xqfnv49NNv2Lz5LNXVen74oZSS\nEvvOCwmBa68NIzjYvMu7bVsJfn4Cw4eHAfDII2O5+uqBZGVd4OWX11FZaT4pdsuWEi5csD6FZNKk\nIPLz68xW6Rg4UEv//iFmx/XtG8UTT8wBYMWK78jOLmbhwtn2/Rg7sdW2la6TjMTEGO66K4KNG4sp\nKPAOMSwmkLX0JZpKOmBapzqVswwUciggmE3GlHa59O7mzVWAdb9ZbGwas2fXMHBgyx3qNTU1bNy4\njY8/Tm+Vb668HFatsm4SnzxZBEBw8E4OHTrHhQvFfPJJPpWVVk9pkk2bGp8gzg8rMtvWp08pHTuu\nBWDlypPk5NTSocPahv1dukQzceKoll28BSgWTT2XL1+mqKiEqqoa5sxZxeHD3p2n9z+q9YwSznPK\nGMXzxklkGjsoPpwmePzxOB5+eFKj7YIgEB8fR0CAaBEZDAays3OoqRETkpWWlnP77es4caJ9BE/e\ncEMoa9f+2aEylK6THbz88qe8804GRqOR/Hwjda1bL85tdKQSf/T0VuXxhmYDD+puZIehi6er5XWE\nhUFoaOMA+IAAgZUrZzRYO+Xl5cyc+TZHjoj5gAwGIwUFUFfnfc+qK1CExoUcPnyMxYt/AODAgbJm\nHX1ehWYUqGIZFnWIzWPnMyPtVjZfVAL87EWthokTg4iMFGfL19UZ+OmnMgoLvePZdDcJCSqmTOnI\nM8/cQFJSYqvKcInQJCcnExYWhlqtRqvVsmfPHgoLC5k1axbnz58nOTmZr776qtEwm6eEZt++g+zb\nd9ps27Fj+bz55iW318UhhGBQ9QQhElT+dI3IYOHwf5NZHMrPZ7uwJyve0zVU8FHCwyEt7RYGDOjT\nqvNdIjRdu3Zl//79dOzYsWHbk08+SVRUFE8++SSvvvoqRUVFLF682O7KuIrMzCxee20Db7+d69br\nOh0hEIRE0IwH3XdgON+w64OZ33CxNIQP9g7ifJF3xVAo+AahofDJJ9cwZsxVZu3aXly29rZloevX\nr2fu3LkAzJ07l7Vr1zZ1mtt5/PGvWLYsz9PVcBxVX9CMgLovwWA+3+Yv306iWqdhyQ3foRbap/mv\n4Bjl5XDPPTtYvfpnp5fdaqERBIGJEycyZMgQ3n//fQByc3OJiRFjO2JiYsjN9awFcezYSWbN+je/\n/FJGRYWPNz7NKBBCQbcdjKVYZtJ7ZNQe/DR63tg+ol0OdSs4jtEIhYVGKiudvyRrq+NoduzYQVxc\nHPn5+UyaNIlevcxzuQqCYHUa/6JFixr+n5qaSmpqamurYRM/Py0xMUH4+ZXhC4m4zRBSQJCnHggE\nYzYYzjZ5eFFVINHBFfSJzaNPbB5rj/TiQkmYe+qq4EWocHRNsO3bs1GpviYwUMuMGWOsdqPS0tJI\nS0uzq0ynjDo9//zzhISE8P7775OWlkZsbCwXL15k3LhxHD9+3OxYd/lo8vLyyM0V5yrdc88a9u71\n8vFqOiAmEa8CIkA1EgSZY9ewEYznbJZw37DfePSa3QC89NNo0k53Jbc82FUVVvAo1qxWSWgcT0gf\nHg4ffXQNo0cPsmtultOdwZWVlej1ekJDQ6moqGDy5Mk899xzbN68mcjISJ566ikWL15McXGxx5zB\nr7/+BS+/LM54LSvzhXiIG4BK4HfgdmALcEa2v5bmEo8HanUEakVB/XT2GrZldOHVLa6L9lTwFCqs\nr9Wlw1nWuyBAWJjAv/89gPnzb7TjeCdPQcjNzWXGjBkA6HQ65syZw+TJkxkyZAi33norS5cubRje\n9hRVVTofiYkIA8YDlxD/HNcgisx5bIXXN0VVnYaqOvFPGiATHfeiAVq7Tnk1jRfAUzBHhWjNWHvp\nOO+ZNxqhpMTIm28epbCwkr/85fZWl9UqoenatSsHDhxotL1jx45s3ry51ZVxFqtXb2LnTl+Jj/ED\nUoBTgORLOgE4JhKrDl1Jp5AKZvZPZ+Wh3g7X0j60iI9Ua8126S3tSrGRlh2Gli097A1IImPAnfU+\neLAOQcggLGw1t9wyrlVpS9vEpMrTp89SWmpaQeC11/aye7fzPeeuJ7P+4zjv7LyaBcP3c1O/Y5y9\nLD4YOaWhXCoLaeZMUKuMpMRcJsivDgTQGVScuhRJVa2tx0WDKJoGxC5gawioL0dqRK5oUNIywtJy\nwM4qX/KZ2CpPcOB60tLGzusatYTz53WsXHmaa68d0n6F5qmn1rBpk0loKit96S3lOpbtG8i5onDS\nHloGwD82jeWfaSObPS/Ir5YP5n1D/y6XQAVFFYFc9885HM2OtnKGCtGZXQmyGeYtpxpRrILqv1fh\nGuvGgOjzUsu+O4LcZ6K3Up4kFL653vm4cSF8/PH9BAUFNX9wE/is0KxZs4kPPxS7b3v2VFJWpoiL\nJbU6NfuzOnP75zcDcHVSDv+4bivPbBzXcMzg5IssujHN7Lw6nZoPtg2ioFx8qPy1ep68fgcdQ6rY\ncTKRxd+MBkIxdUGMQAVid8/Rv4MkLM5MPK7C9KgbcE495WULmOotv5YcaSRIsNhvT/dN7pfxzHOu\n1aoICWneGraGzwjN/v2H2LLFlBx827YcNmyo8GCNnEEs0AM4ALgmbejlykA2pF8BQJ1BxZjumTwx\nfqe4Uw1+AXqO5USZnVOj07DpaHeyi0IBCPKro0tUMbHh5QT71/HwxL18uetqLjcMnRsRLQRnvKkl\nIZAsAGm41lHHthqTtSHVU/pXRevqLgmAPG7FmhBIQmMZIytgPlRtGQMjF3PPWEITJwbyhz9c4VAZ\nXj17u6qqimPHTqHT6fnyyz288cZFF9fO1QQB8v5tEpAAbKClI0ytZWiXHJbc/ANgAD8jv2Ym8PgX\nU+w+f0KfDJ6evp3/bB7G7lMJ5BW7Mk5H7rOppPVvcxWihVRD48baWgers7pCasyFxrIe6vrvnsss\n8PbbfXnooZnNHuezGfbOnj3P9devoajIiK5NjHpeAVwv+34AWI07H6J9mUmMXfIY4giXDoOxZdMV\n0o4lczI3kg2LvmDJN0P54IerXFJPEcln48qEXpKVoaFljlZnvZ8t//ZqGsfI+L5bwGuFZv36Lbz6\n6u765EOuvpr0VnHlm2MUEAx8KdtW4sLrNY3BKFBdJw1Dtxy9QaBGp0brp0ejcYcpLzX+QESLxBX3\ny5sasm86i5vDa4WmU6dwUlM7kZoKK1fmcPKkMx4wy/6wHOlhk/eJnfEA+gP9Ed/MWcBp24e7FMnn\nUYXDD7Pb5m1KDlTLLoarrgP2/d1dVRfvEb3AQLjjjigGD3Y8oZrXCs3w4VcxfLholhuNS9m2LY+K\nCgOHD9ehd0hzmpqwLnnzpQcazBuipcOvKToCIYhv4DzEiN8QRD9DCrAH54qMpQPR2jGC7P+a+n+V\nlRXNkSxZtcW25s6RjnEkPsZ7CQ4W+OMfJ9KnT6/mD24GrxUaOf/4x90YjUYOHDjKhAlrHFhvyYjt\nuAz5fnlfWQqUsnXuCGAwou/jY8Su0iDE7tGnQJH1U1uFFOFq6zepMfk3NIgjQ47EubR1JLGxFeLv\nyPHtF68edbLk998Pk5q6yk0LuwkW/5dER+4wSgbG1v//FOJ8JT/ErlIWkFt//EWcH3gmt1QsrTS5\nVaaXHeecIVI/jZ6rUi5y8zXHGJRykctlQSz832SyC1yVlkKNOGLXVACfFOnbHCqaHnWyRLICpXtq\ny2cijQg5q5vtPYweHcCLL45j8OC+BAfbN7Los6NOlsTGRvPUU935+ONzTvLZ2MJo8X9LH0EKEI9p\nysA5IAfRJxOOOCnSlcPx8vrJ/QvSd7AepeoYtTo1u44lEN2hkm6di7imXyZB/q4cFrSM5JVbndLQ\ndHPYG+xmGQ8jiY7lNSRBansiM2VKEHPn9nPqKqA+JTRxcbE8/fSdlJcv5auvcjhzxh3mqoDof5FG\nPS4DcYjWTB2w1eL4GmCXG+ol0dSD7h4z/ptdVxASVMvwPheaP9ghjIj31R9zK0KKY3H2sKR89FGK\ns2nKD9P2RAZg6tQu3HbbVKeW6VDOYE/x8sv38OCD3REEzD6uQQuMBCYAfRCduzcAxUCaqy7qEzSM\n0bitrdUgNuwgRNFpzm/mDKQAOsmRDqa4G/fOovZlfMqikTNr1liuvto07PbFF/t49918J18lClFc\nDiE6USOAVOA3xK5R+2bhLb/SNbaYmxfdygWX+Wcs0WHqBrmrkUuCIlk3Ci3FZ4UmISGehARTqkuN\nRo3BsIsPPyxwcPhbjh9it+kgogUjme/FiJMI2zc9EgrpEFzN7mMJbryqp4LZmnL6KtaMvfis0Fgy\ncuQQwsNDycpah04nPgDFxTr2769rpWkfiuiXycNknpciJqUKQhnOxHb8Y5uk7UXsytFoYPBgP5KS\nms8P3OKynV6iB+nduyffffckIC7a/vPPu5k48YdWCI0AdEMMuNsh2+6H2H0qwPkOSAUFzxISAu++\ne0OrV6q0hU86g+3hk0828NhjP2Fo8UtIDQxBHE49bLGvFijEZO20c+RZMRUUbNBmH5MLF4o5cqQ1\nIxIC4vQBPWJXSY4Uz1GNOBrVumxjbYp21XVSaC1e23U6fz6T9PQMm8d06BDMsGGDUKvNp9Xv3XuA\nU6ccCflvbthSh3jrrC150T44mhFNj4RCxl2VwZ70BCqqIzDPjFeO6EBXaO94TGh0Oh21tdYTiK9e\nvYOFC0/YLGPwYC3ff5/cKI/pX//6A1u2tGbioDwK1FqfS8CUH7e1SbjbBm+tGsrsiUf47Nk1pD4y\nn1MX+iNmDZQ4gTgtw1b/1ZpHuW07XtsbHhOab7/9mZde2m11f15e892e48d1TJnyLiqV+YN64kRr\nV0DoCPRFnGV92coxBsTh7SBEX41bJl55LwKg9gNhAKLw7pftjEO0cGxZplFAIuajeFWIi+cpI3tt\nBY8JTX5+KXv3OrYkSkWFkf37HR39kUaYghBFJBtRZGzNcq7DlP0tFDGmpn2+gX8/Fct/v7qae6b/\nyOq0OPYclQfuaRHvjzyfSSkgBVbGIN73fBqn5UhEDC1o31aju+jTR8MDD/Skc+dOLinfY0KTlBTF\n+PGB7NxZRbVHMheEIGa8k7pCUiCevTljpHB4KeeMq7K/eTfHz0VRUhnApk8/5UTeCPYcHSjbexlR\nQCwD+iRRCUUUkhyL/QFAF8S/Cyhi4xo0Ghgxwp/wcC3DhnXikUducd21XFZyM0yePJquXTuTmvoZ\nOTmeiLBMQFyBoA7YjSgyLaW2/rwoRMFyJIG272LQCJR2DaSmQ1OPUxHmuXhiAWnlzOOYrBs51Yj+\nnZ6IL4AONQXSAAAgAElEQVQs2qOIuxp/f3jxxfEMHz4IPz8/l16rzQ5v20cx8CtisqrWYkCMrVEh\nWkgKtrkM/F7/aW5k8Dzii6A77f5RdQFVVXDffT/y1Vc/uvxaHh3ejo6O5IUXBlBeXsPOndl89ZW7\nHas6nOPM1SF2nbSIXalKnOOzkUZkvNf/0++qPG6ec5xP3+3H/l/j7DijDvujqqsxWZpdEBOJtbQb\nJeVJltNcpsX2gcEAJ0/qyctz5EVrHx4VmvDwcO6550YArrxyO0VFv7JtWxU1Na5Yb1mekU6aietM\najFPYSA5utW0fmVEa0mXvIfuPYu44baT3JQ6k4xT4S64QgXifeiLKDrNCY1lxkFrM67lWQjbX3fX\n3XhNwN7kyaNJTIwhNXU5eXnOLFnAlDBJjnwxeWdRh9gNkze4IMS8wZJYtJ1kSUHBdWi0RoqKAjHo\nXRkiLFkg0hrXtoTXMpBSR2MLSoU4Yijtb9/WTWVlHWVlZYSGhrrsGl4jNLaRsqq19M2uxpQZz3Jo\nqxLX9PsNiMIShChyZZivUy1NbWjut8gXjfdOFiz8nbiECu67+TouXmj9uszNUwWkIzrwQ7G9fHAd\nJmsSmhZ1A+YWp5b2PEn27bczyMv7nDfffMBl1/AqD1tsbDSLFw9k0CBL/ZOvkWwPfohDpFrEB0qH\nKdpX+rjyTabHJG51iA1F+tQiOo2trb4oJUKXhNV7rZ+YzhWEhNZy5kQEtbWunI5hQHwxSCNU0ZgE\nxfJj+be2dv/k0d+W1pj0N5A+XtVMnM6lS0YyM12bX8mr7mBERATz5t1Ajx6WkxWl7obU/7b8yB8K\nLeaLr0u+E0tc3YWRHnwjouBIUxZqEEWwqYYpD8f3hS6WgNGtsyqLEO9lINaFpjUWoNyv0+6S7LgF\nH+k6gekN39SKhfIE0hrEBm3rgZP8M57qm8vfpHIxkR547/YZqFRGwjvWUKtTU1zi3/wJTkNyrJfh\nPBGWr1AqHyTw3i6rsxAE6NgRtFqB8HBXrm/uU0IDza+NLWCK2LV1TDJiN+acsyrWAvSIoycBiH4c\n+exmV6797TzCImp4/ZOf2PZDIt9+5fhyqfbTCfH+XMJ5QqPHdM9bsnyL7xMYCO+9N4o+fZIJCXFt\nDJhXCs2CBcPQaPayfHlT0bq2HrDmHr5AxMjUCkSHrKdSGOgxrVMk+ZF8YyGyfoPzuXX+MbZ9n8j2\nTYnkXXRHkKIGMbNhDeILwpmOW8v1sSy3OZf77otiyJDODd/fffcYv/3mWkf0jTeGMnVq10bbtVo1\nY8YMIioqyqXXBy8VmokTR2E0Gikt/ZVNmyqwkU2iBQQhJrRSYerre5I6xLenH97u9JXoNySfYWNz\nEATYuLI7uTnuioRWI4425eLadcNdb8mMH9/dbM2k6uoVqFSn2bfP+WIjCDBpUhDz5g1i+vTxTi+/\nJXil0ABMmnQNiYkx3HLLCoqLDZSXGyhuzXQkQPyZHRG7K/ZOmnQlch+T9wuMSmWkY3Q1N8w5hX+g\nnv97YIynq9RmeOyxWYSGruPZZw8CcPmykaqqlj0TYWEQFtZ4XEetNvLSS9cyZMgAp9TVEbxWaAC6\ndevCN9/chV6vZ9myrbz4YmbzJzVJF0Szu7XnOwN5JKof4q3XI76hvVtsgkPrePXDNHZsjmf9Fyme\nrk6b46abUhkzpj8ADz74NZs2tWyaxZ13duJPf2q8sqQgCHTuHNvEGe7Hq4XGz8+P5OQkAG67rQa1\nWsVrr51vseKLDbsG80AudyLFxcj/X4er1sZ2Nmq1kbjEcurq1BTkejopeyBipHcJJsvQe+/hffdF\nMXy4KU3G0KFXNjqmQ4cOdOggpsR46qlRREbu4csvS+y+RkREAN27N/bBeBNeLTRyBEFolEnP95CG\n4CWB8W5Lpv/QfHr0LkajMbBnWxxnjjljLlNr4lSMiPmHNYjOc8sgOmvhYO4XoFGj/EhJMUVJz5s3\nghEjBtt9/oQJozh0KLNZoenUSWDKlHBUKoFBgxJbXV934TNC8+OPv/Pcc7aTlVtHcrR5KtTcgHkQ\noXcPYatURqLiqph6SwbX3nyOshI/Hr9tHBknOjR/crNYTnq0Bz1itr0oxEdWmikvWaiW5bk/6FGj\ngdhYFY8/PpiZM6916bXCw2H06FCWLn0YjcY3mrBv1NJhziOmjeyC+53BUji7rYTn3kVAkI5/vL+D\nXT/FcUfqdRgMgpO6TK2dvyWJhZS/RoMpX3NTkd/SPXcf0dECy5dPZdCg3s0f7CC3396JZ5+d2Wj1\nD2+mnQiN3KJxN/I5S97dVZo2O4MxUy+g1wukbUhk99Y4cjKdOVnS0TgVSaAsV66zLM8y2tf1Aq/R\nCHTuHE1wsOND/lOmDKaioo6XXjrbKM3tE0/EM3v2cGJiXJPb11X4jND07ZvIjBnZbNhQRl2rej8V\niA9eNGJkrju7UN5tyag1RlKnZdH9ymJqa9TUVKtJ25DIxUxnx8k4Q2glP001tu+r1F2FxlM9vJsr\nr7yCu+4KJDe3gupqc+vvrruuoW/fxg5lb0cwGlu+MrVDFxQEWnvJ9PQT3HHHSk6erKOiVZNNQ4EU\nxKU8yvB2AXA14VE1aP0MlBVr+STte777qisf/cv56y47l2DEUaci7BMPyTfmXOe7RgPx8QJ+fqKY\nxcVp+eyzO0hMtEzE3n6w1bZ9xqIBSEnpxsaN93LXXR+3ONZApBZxnkwQohle7tT6+Ro3LzhN994l\nvLBgKAZUbp6J7S4ky0aN+Dd3jth07Cjw2WfX0a2bKCxqtcotofy+ik9ZNBLbt+8hO/syZ88W8MIL\nZ22k/tRgnjDcgNhlklaadH2uVG+jR78S7nn6GABn00OpLFHTf2gBv/3SiT1psWQcD2umBE8jH+Ku\nwn7HsmMpUcPC4NlnryA+PgKAwEA/xo8f6tKsdL5Gm7FoJEaPHgqIa2y/9NIZK0dp6z/yH65H7Ntr\nEM3pADw/58l99BlayOCx+ZQWiU7xI3siyb8QQGL3MrZ/F0/OeV9YxUFKbCUlNpNSfjQXjNkykfH3\nF7jhhlAiIwMACAnRcvvt44mL845IW1/DJy0aiSNHjjF37mpOnqyjvFEvKBgxItjakh4h9ftLcaZJ\n7a3EJVdyw/yzRHeu4h/3Xu3p6jiJYMyncjjPDxMRIZCWdgv9+7t+uLqtYKtte1WGvZbSq1cKP/xw\nP8OGNRXjYcC2WV2B2H2SgsDaLoLKyF/ePEBtrZo3/jyw+RN8hkpM6T5CcXfsjIL9+HQL02g0REVF\nNXj+zalFFJsQ2XfLpNW1iHNmghE1V8re1nbo1qeU+xals/O7WPanRVNe4olYIldhRHyEpayK7XsU\n0ZvxaYvGNnpE4ZAictWI/Xrpo6nfXoWp3y9gEh1baDHNwPZuwqNrSb3pImeOdODcsbbiuJTy+Phh\nGk2SXiwK3oj3txSHkLLngxh7ESDbJwVxGTANc/shjkjJ0ztC46heKQG6JFDWjvMOxKHrtoI0vUBb\n/39LS9VxwsKgTx8/AgJcux51e6KNC40cywcyANF6kSwaqStViJgkS35rpDzEUkJzyfEoDZ9L50tW\nlPdgREBvdbVGX0QaaZJeIM6X0NGjg/n88wWEhLhyrar2hU8LzenTZ/n739eyf789wXuWFom8ayWf\nEyMt8CZvmFJkqQHzuUuS6Bhlx/nR+iVwnY8of2ovqU1LaWqFUWnNLNf8ovnzI7n//msa8sMoOAef\nFpqiolI2bCilssVBwvJE4PI5MVJQl2XScrnPxrKrJF8axRsXG5PSMviCRSPdP3ldLf0ult3a5klK\nUjFzZoxdx95ww0CGDh3UovIVmsem0MyfP59vv/2WTp06cfjwYQAKCwuZNWsW58+fJzk5ma+++orw\ncDEh0iuvvMKHH36IWq3mzTffZPLkya7/Ba1Chym6VEpCZWupDXk+meYm6HmP7RAZV0NkXA0nfwuj\nqtwbhn6bS3qlQnwkJbGuxhm5e664wp/XX7/f4XIUWo/N1++8efP4/vvvzbYtXryYSZMmcfLkSSZM\nmMDixYsBSE9PZ8WKFaSnp/P999/z0EMPYTC0pVEAqbvU1AJ2Eu5JSWAvk+/OYfLdOTw+bhhnDnnD\niJMWsTtk7WPAfFVP704QpmA/NoVm9OjRREREmG1bv349c+fOBWDu3LmsXbsWgHXr1jF79my0Wi3J\nycn06NGDPXv2uKjankLy01jrikgOYcs3s2fQ+hnQ+hmpKldj0HtL10kaim7qI8/b4xzL8I47Injx\nxUlOKUuh9bTYR5Obm0tMjNjfjYmJITc3F4CcnByGDx/ecFxCQgLZ2dlOqmbTxMVF86c/JfHFFxfI\nyHCXJSH36VgiNRB5HhTHJvO1ltRZuYDAli89PTdH7iOSO9xdR2Kiittui0MQYPLk3gwbdpVLr6fQ\nPA45gwVBQBCsvymt7Vu0aFHD/1NTU0lNTW3V9RMS4nnxxfn89ts/ychoLkGNPI9sU8hHjpp7mzYn\nHNJ+KROcq+NrGtf5ugXZHNkezoZ3PZkfRZB9pHWzWy+6PXuqqaoykplpvYz4eIHrr4/i1Vfvtfls\nKjhOWloaaWlpdh3bYqGJiYnh0qVLxMbGcvHiRTp1ElMKxsfHk5WV1XDchQsXiI+Pb7IMudC4D+nN\nqrOyX7JUNJjiYnwBAe+ts2TZ6bB+3+1n8eIRHDt2kaeftjZjH+bOTWDRorsUkXEDlkbC888/b/XY\nFjsRpk+fzscffwzAxx9/zI033tiw/csvv6S2tpaMjAxOnTrF0KFDW1p8izh58jTTp/+LvXul8W2p\n0TX1kUaUNIixLtLaSnIMiA3CMidta5F8Nq4SAKmeJpFJ6lPJ898d4aflcWz+NM5F1/UMWq2a2bPH\n8t//XolffdDubbd14IcfJjR85s8fj1bbluZztQ1sWjSzZ8/m559/pqCggMTERF544QX++te/cuut\nt7J06dKG4W2A3r17c+utt9K7d280Gg3vvPOOy98qJSXl/PRTuUUcjbVGLU+MbZB9LJF3cxz1r7jD\nwjDvNoVE6Bg4qZhV/0wg54wnF3uThNrxexAbK3DnnZ3p0SOB5OQkbrpJTWFhFTqdgbFjr2DSpGsc\nvoaCa/GJfDTHjp3k8uXGC2odP36BRx892ChTvHNwTa5Z5yLV0eRc7X1NKS+nHebZyX05tMWT0a3S\nipwt908FBwv066dFoxFfVD17hrBkyb0EBnp6lUwFW/h8hr1nntnA2rWN0zcYjUZcF6rjCz6bxj/e\niIDBK6Yc6DHFHLVslKl7dzUbNtzfEAgqrlLqbRHXCi3Ba4Vm7drNLFmyH4DDh6vR6z3RdHxjPSaJ\nsXMKGDipmEXXXknG776QlrMxM2aE8Ne/TqBDhw4+tUCagm28Umi+/voHli07xJYtVQ6Vk5SkYt68\npIbv586V8vHHhS0sxbtFJjKhlnF3F6DWiIJ4ZFsYh7d0wL0dYueRmBiizDVqg3ip0Bxl40b7F27q\n319DebmBggIDgwYFoFKJfft+/Trw3HNzG5zSv/9+mPPnv+X336spsXD5dO+uJinJlH/k9OkasrK8\nZzqBNSITapn5f9lkHw9k45IYfvow2tNVktEytevbV0OPHh1dVBcFT+J1QmMwGFr8Nl68eDSHDmWx\nfn0WP/64ED+/phMWDRrUj82be5Oa+hq//GKaoa1SwYIF3XjyyTkN25544j1efz3HRywDgU+eTOTg\nj96W2qBlQv3cc1dz883On4hrNBrNnJTNBZoqOB+vEprz5zN59NEV/PpryxeHmzNnHFOmFLc4hkKr\nhSVL+jFt2nCz7Q89dC2Rkdv5299Ot7gu7kR0/vpKGgjbPPvsXrKzi/njH2c5tdz9+w+xcOF3Dd+H\nDu3Iv/61wKnXULCNVwlNZWU1u3dXUlBg24y4444ICgpqOHy4invuSeKKK5JISIgnIaHpSOSmmDgx\nkLFj49BoVFx//Qji4zub7e/atQu9e58BvFtowJfc1bY5dkzH55+foaDgw4ZtI0d2Z+rUsa0u86ef\ndvDRR/vYvt0UA5GTk0tgoHiNa6/txzXXtJXlZ7wXrxGaCxey2bkzndpa600mNBSGDAnk4YdHc+FC\nPmlpZ3n66dvx9/e3+zqCIDBkSEeGDk1k9uzrbB4bHd2BceMC2bevmrKyttCUvZ+9e2vZuzez4fvN\nNxej0YijT4MG9WrxsrMXLlzm2DHzRb/OnNHz4oviNcLC/BWhcQNeIzQbNuziwQcPW92vVgtceaWW\ntWsfJCxMXLZ15syWX0elUvHGG/YlQRoxYjBr1vRg0qS3OXiwFoMBdI5P2XEaKrURldqIvlaF0fv9\n1q1i1apSVq3aDMCGDXVMmzbO7nPr6uqYPXsKMTHhTJ26xVVVVLADrxGa5pg5M4y///1agoPdGx8S\nEhLC0qUzqKioIi0t3at8NiNvL2bQtBJevq4HmYd9JWpWPo/MtSuEPvbYexw4UEJJSRtVYR/Ca4Rm\n+PCeLF5sPW5myJCu9Ovn/uVJ1Wo1/fpdCUDHjmFUVdXx7ruZ5OZ6visV0bmO6C51nNkbTF2NLziD\n5dM67D2+5R6os2fP8d57m/n228vNhihs3JhFdPR67r57esO2o0ePs3LlLhYsmKKste0kvEZoBg7s\ny8CBfT1dDZtccUUP/vKXGPLzP+X8edtxPocOVXHhguvfpJ6Xu9Zgj3jI8wuLx+7bdx5BSGs4YujQ\nvk36bLKyLvHaaxfsCk1IS6smJOQEd98tfj927CQrVvzKG2+cJzr6F6ZNG0qXLkk2y1BoHq8RGl8h\nNDSUd955qNnjHnroHT75JB+j0UhVFS6JxzGi8pJ5TfYizR9TYzs/jZQoy7xrtWjROeBcw/e1a/VM\nnTq6IW6quroavV5PdXXLFpTT6YxUVIgvjnfe+Ym33hKzRj788BGqqmp55JGbWzTgoNAYn5i97Yuc\nO5dJQUERly4VsmDBNi5edP5vvv6pAgbfWMpLqV3ReaTr1Hj2uH3n2EpAZv+s+ZQUFY8/3puHHhJH\nBR566B327i2irMzAiRNN1UmaUW5uaXboACkpYvxVVpbOrFuclCSmBX311fua/2ntHJ+fve2LJCcn\nUVBQxPr1R6mocI2wih0QT/pmpN8lTX50VkoNe9KpwqlTBj755CRnz74HwLffFthM82mykswpKYF9\n+5peYTQz08DatZcICPiIhx66riFftkLLUITGyVRXV7Njx34qKmrYuvUU779f4LJrGRG8RGgEi3+b\nW/dKnsAdml6/XGXxvWl2765l9+4cO+vb3LpSTXPypJ7//jeTTp22MW3aMJKTFZ9NS1GExslcvlzI\nAw9s4fTp9rImkZSuFOzLPyN1XdSIj5+0gF9T+6UulDPupWOCXFJi5JFHjlJTo+OBB6IICgpyQp3a\nD0o2IQUn0bT/o2kkcbJ2rLRfSjgmz/vcGlSIeaJ1NC9akrO6ad544xTPP/95K+vRflGExqfxdNdJ\nrEPr0nbKE7dr6z8qi/3y9dElK6elybCkcyRhs1U/6frWxfLCBdFn88wzH5KXl9fCurRflK6Tj9Jj\nVBVqrZFjacEenn4gX7PJohF36AxB4abvch0yAiUXoKaChgYelgDakMbtXF8FZefrz5WC+FoS9KdG\nXGrXnmOlcqVuXWNhOnlSz9tvZxET8zPTpg2la9cudtal/aIIjY9y/dOFXDjsz8q/tmySoXOR1sqy\nMIwFATQB0GUIxPYybZcvPKEHTv4ABadM+xOuhpDu4uq4cqryoHoNGOoADRjVQH0+IUsdMNaJH7P6\nNWf1CRYfNeJa4FVY62oVFRl59NGj+Purue8+RWiaQxEaBQfww/Tml28OhoEzIPc4ZOwybZeLghFI\nGg7JsqVSzu+Ec9sbi4dfOCTfBqr6x1UvQGX9QXLxMgAVe6FSWvM9oH5jDbbRAqGYVtLUI4qMf/22\npoe+FexHERofZcs74XQZXMPtb+cDAlvfCuPisaYzC7aWXjONBHSEA8v8YNIIOJ0JJ84hvvWla9XS\nqBujUkNwJOjroKwIscGD2H2RNdq8o1B20fS96BxUFTWuSF0FXN4HQr3lZJAVIxclA6AKhuARUPEr\npq6QrW6WP6LQVGLy4UjnWJ/02aGDwOOPd2HUqD42ylaQUITGyQQGBjBtWjS5uVVkZFSxe3fLwuHt\n5dC3wQhqGDpbzLVy5cQq1Fq4cMhcbEI6Gek2XmfmNDYa4OwWNZUF5l2K2IEGonqJDcuIQFRf0AQJ\n9LwZzgX7U2NIhNAQRKHxB4NO/Fi2RW0QxGghU4soSFLXSnL41p9wOROQcs9IYmDp7DWCrgIu77fv\nxgRcCdqu4NcX9P5g0INRGr3KQBQUCb/6OhkRLRgJqQuloymRSkpScd11kdx/vzLp0l6UKQguZOXK\nH3j44V1cvmxE7+KwmvuW55F3RsuWt83zBscPNjDjvSoMqBr0wFAHaxcEkHvI3Lcy/HEdA+7U17/T\nVaQ9p6H0opYJb/qz6vpqLqeMg5GyFQpqacKgEcCoEhv4t7/A4aOAtCZXKGLjlm6GfDkba1MS9Bb7\n7Hh2NIkQdCvUCaLlo5PK+QbIkpUXjmhlWU6QlfxOTVs0t93WgeXL/9R8PdoZttq2IjQupLS0lN9+\nS+eOO34gO9u1vzkiQcXgu6sZdKfpjW1AxaXDGn54Wuy6NEwY0MCkV2qJ7GU0s3R+/1DDsTUmi6Ii\nX8AQ2Y3gSRMo/Xo1+ppqCApo2N/0iHYwYgPOhbJSqKnBpESWjllLp01TSIF9EjpsT8hUgeAHQrD5\nCBcqYCJwETgmO74cc6Fpfq6VIjRNo8x18hBhYWF07doZrda+uTuOUHRBS/oGI6WyCYFGoCxHTcFJ\nc8tFUMH+pRqCLabtXNit4vJJixGaYH9qcztCnQbKKsSPVUKl2iA2XstuR2vG4Q00tmhsrYluBGMN\nGC2HswXgNyAaSEQUG7mzB8zjaKz/vfbuLeeJJ97jL3+5QZn7ZCeK0LiY4OAgbrwxhvXrczl71hUB\nL9IQroGcAypyDjSfzsBogOPr7fjTB3UGVTRkZkGdPb6mgPr6NF6+uPU0FTNja7jamkAYgROIFldH\nGneXpHKbfyHo9Uaqq/U+shSPd6BEBruYqKgo3njjfqZPj6FjR2dH8Uoio0X0JzjLESSAJgTC+4Ix\nHvbtgsoamn9cmhvhaWVdGnwm8hSgrkAq1/bvHD48lCVLHiQ2VrFm7EURGjfx9NM38eST3Z1YohpR\nYKTgNSc2cJUWkqaDrhwubgBygChMXSNruEJopGFm6ePqMGhnTeJUkKN0ndzE2rXb+eabLAdLkd7u\nYD5r2gEb3j8GOo4wL95ogNLTUJ4BurL6HVI8TET9/+WogBBEX0pzwXGtwVk5bvog1jHdCeUptARF\naNzEL79ksWOHMxthS+b7WCJASAqog0EdYL7LCBj1UHIC6kpkG8sQ/Rt+NBYa+cneSCDQGVEoLyJa\naAruRBEaNxEdHUh0dDH5+Y40RrkV00oELWjDIHwI+EdD5Tm4uNbOk/WIloHlhEMDUIpo1QTgXH+R\nM+gADAU2A/kOlRQVJdCpk68sbeM9KELjJp5+eiaxsd/zxBMnPVuRwM4QdyNc2gi1efUTFe2lENGq\niQMu0XgOUCViSH8HoATvEhvn8Kc/deP++22vcKrQGEVo3MTXX29l1apzHrq6ACGjQR0FVEL+FqjO\nA305LRMDPaYQ/nDEYDd56L4B09TrYMSoW9dMwbCfSERL61dE/5EG2wF/tomICCQyMtI5VWtHKKNO\nLqa8vJwvv/yOjz8+ya5dbm50mjgIGAxBV4EQCOhAVwKlh0FfTesy1ukQY1ACMU2slKPHNAqmrT9O\niq9xN+FAEKI4nquvW0sTZ5mzY0c2W7bsdLhm7Q3FonExly8X8re/7eXcOWcPy4aJ/hYwb8Py//v1\nhIB+oDJCyWqokztB5YnFW+M30mHy11haCAZEaycE0bKR5ivZugdGi38dQUrdGY3oO8qV1UvaL/9u\nP5s2FRMdfZjx40c6Wsl2hSI0vop6IgjdTBkbpFzeYHppV++D4o/E/xstlxuuTyKFlpZ3b/SIjbcj\n4nC3NQdrBaI1ocK2VSM5ueUfRwgCrgDOYh6lXI14s4Ixzdhu2bUWLuzO/fdf72D92h+K0DiZwsJC\nnn/+K3JzxYZdUaFzcKRJYhTQWRbrkgmGk2I2uxoNCPWpEOQriujywFBupTypcUtvf8v5Pc3l17Un\nilY+69JWV81o8XGUGuA8EFb/vVh2HbllJaWx0GOv2H799XkCA3/gscdmOaGe7QdFaJxMVVU1Gzfm\nO3G5lSCgF6Kvo1LWHjPAmFefkkGLKTtcS5COb6r7ZGuNJqG+XpI/pjlspYFwBXXAZURHsKWAyC0m\neapPqSnYjm7ev78OOI1GsxKA1NR+9O7d01kVb7MoQuNVSKsBVCE2zkAgHhgC/EjDutNm7UAK3Gut\nJWCtYUn9L8typblVITROseBrSCIk/SYBe9J27t9fx8GDR0hOVtOpU6giNHagCI1XEQMkAPsRH/je\nwADgK5w7I9oerFlk0jSEAjw/dO0s5MPy9gl2ZKTAZ59NY8CA3i6rVVtCERqPEoWYG0WiEsgGBgKx\n9d/TEEdOmltsraXLkLSGYMQ3fyFNLq/iNQQgTjkoRLS67KFlv0WlEoiICCMgIKD5gxWUOBrPEYUY\n51Ej+5QgziMKQWzIOYgjJ82Jh9R1kk+6dAXyhdyCEBu0N76rNIj31h/PxO8oWOKNT4nPUFxczLlz\nF8y25eUVUVNjSxjqE3uThGiuH7HYH4howWQgvpHtxYDJsWnLkesIVfXXkfIS1yKOJlkbLXLWKFJL\nMSBag8GIgm05quR4xsO6OiPHjp0nMjJCiRS2AyVnsAOsWbOJefPMo0QNBiMVFWCwqjV+QDJi16iQ\npoWmNy0XGgn5CIor8rdYCpkG0bKRVhqQIzVydyO37KIQ63e+/rskxo7VSxAgJERgyZJBzJ073aGy\n2kbKL28AAA0SSURBVApKzmAnUlZWxtNPf0ZmZgU5OTWUlDQnmirEeI5ATEt75CFaAraGhqOxnWgq\nD/N5RhJOylNjFUsrRYdplEyOuolt7sJylrt8zW6p/mocGa0zGqGszEhtrTuH7X0XRWjsZOvWX/nt\ntwwqK+tYuTKPS5fseUClt7304EvDp2VY71boEKNum3PsBteXV9nEPncuxt1UjIzc6vE0esyHrOVT\nLxTchSI0drJu3UH++99Ldh6txhRx648YmSoJgiQ41uJX6oALTWy3JAbTPCJXZLVzBGdG+TqKDtE/\nY9lVck6X7vz5YjIyztO1axenlNdWUUadnI6U1rIjotAUYnqrGhEtnFDEh78pa8Re8hCFKhbvfDtb\nX07Wvbg2B/C//pXFSy+td1n5bQXFonEqUlepCtEHI/cBaBAFoQZRdLSI1o7l+kP2IqXX1CMG+eU7\nUFZbIhBx6F2iFnEUz3mMGRPAU08Nb/geHx/t1PLbIorQNENNTQ1ffPE9+/Y1lSdXwNRNkvwScutF\nQp6WQJprI53rj9gYWvP2r60vuxOO5lnxTaR7L91fTf3/5T4Zad1exwgKgttvjyIqKoD+/Ttz3XWp\nDpfZnlCExgKDwcDRoycoLxe7NRUV1bzyykFOnWpqETO17CMtCt+UVSEJkFxMahFFJgDHo2y9pZvi\nCiy7hfLv0r2XhCaAptfSdozISIGhQ4N4+uk/KL6YVqIIjQVVVVXcf/9a9u83OVjrmpxnJ+8KNUfb\ny53rPiRBl5BbMJbWSg2uENyxY0P48svH0GiU5tJalDvXBHV1BmptWttSPIyjAuLozGsQG1oO4giU\nCvvn9oCzgtdci+XwuWCxT47zf8fDD8ewYMFYtFqt08tuTyhCI+P06bN8+unPXLxo7YGVv12dtSqj\no29gA6LzWepCBNdvr6Pp9JnyGBfHQ/Fdj2X9bNXX+b+lZ89I+vdXZmg7SrsXGoPBwOHDxygtrWDX\nrlO89FIm/fpp6NbNj8uX60hPl4uOFNrurJnLzow1KUf0UUgjLpZlywPV5M5pdwb3eT/duqmIjzcl\nXe/cOdyDtWk72JzrNH/+fL799ls6derE4cOHAVi0aBEffPAB0dHikN7LL7/M1KlTAXjllVf48MMP\nUavVvPnmm0yePLnxBb1srlNNTQ3jxr3O7t21GI0QFCSwdesNDB7cn2++2cqNN2534dU19R/X+BZM\nOYHB5DB2d7Y73+L113vy+OOmNJ2CICAI3hin5H20eq7TvHnzePTRR7nrrrvMClu4cCELFy40OzY9\nPZ0VK1aQnp5OdnY2EydO5OTJk6hU3hkTuHTpWj777DgGAxw9WovBAMOG+fHPf06kV6/uqFQqrrlm\nIFu3il2Rt97ayapVzo3HcP2C8vJJld4j7t6MIOC1z6wvY1NoRo8ezblz5xptb0q11q1bx+zZs9Fq\ntSQnJ9OjRw/27NnD8OHDGx3rDWRkFJKWZhqKnjQpiLvv7sPo0UMbtkVGRpKaKqYAqK6upV+/s43K\nqaio46OPcigoaG1DdmWovmvKDg4WmDevE1FRgRw5UsjKldYFeNAgLWPGRPLRR5eYNCmM2loD33wj\nOqxvuSWMPn06NjonN7eSjz7Ko9qN8Yf+/gJ33x3FiBFKWk5X0CofzZIlS/jkk08YMmQIr7/+OuHh\n4eTk5JiJSkJCAtnZ2U6rqKtQqcTGMH9+P267barV46ZMGcOUKWMabS8qKiIv71OysqrIza3j6FHf\n75ZoNDBokB+hoU0HAYaHa3jiietJSkpk+fKNrFy5B4BevdR07my+qNyUKQnMmTOeixeXs2DBMKqr\n66io2A3AQw+NJDW18YvozJkMsrNXsXNnhQMCbj9RUQLDhgXx5z9PISWlu8uv1x5psdA8+OCDPPvs\nswA888wz/PnPf2bp0qVNHmutb7to0aKG/6emppKamtrSajiMIIBaLeDvb+Ttt6cybNhVrSonIiKC\nZcseA2D58o3ceefeZs/R672nGyN2Fcz/TqGh8N570xk4sK8d5wuo1eL5TzzRj/nzb2zyuBUrTF3t\nadNSbZbZvXtX1q37Czfe+Drr15c1qp/RaLSR76dlqFQwcmQw69b9xTkFtiPS0tJIS0uz69gWC02n\nTp0a/n/vvffyhz/8AYD4+HiysrIa9l24cIH4+Pgmy5ALjae4996JXHttLiqVQO/eKU4pc+LEq0lL\nsz3vJTu7gEcf3UN+vnzKgucQ/VKTzLZpNGpSUrradf6ECUNIS4sCICXF+VGzt97agUceucZsW3p6\nFo8+eqiZWCf7eOCBTjz88ATHC2qHWBoJzz//vNVjWyw0Fy9eJC4uDoA1a9bQr18/AKZPn87tt9/O\nwoULyc7O5tSpUwwdOtRWUR6lS5ckunRJcmqZ0dHRDaNxTbF//yF+/jlD5ntwNFjPMSZPDmLevH5c\nc83VrS6jud/sCHPm9CUiIqRR/bp1S6CoqBKdztys+eij85w5Y59wa7UC994bzZ13DleWS3EDNoVm\n9uzZ/PzzzxQUFJCYmMjzzz9PWloaBw4cQBAEunbtyrvvvgtA7969ufXWW+nduzcajYZ33nlHGRa0\nID+/mMOHS6itdU+ulu7d1XTv7m91//z5/Zg1y7pfytPccsu1TW7v3DmOp566o9H20tL3WbXqImfO\nWO9XpaSo6drVHz8/gT/+cRI9ezrHmlWwjZIz2M2cPn2W1NTPyMszAAJ1dY7dC5UK/Pya3vf3v3fj\n//7vrqZ3tlFeeOFjnnsuw2yb/B698EIKTzwxxwM1a/soOYO9kGee6UplZR2LF9uTTc86vXtr+N//\nrkWrbTxCFB8f41DZbYXBg/1YsmQKAImJcR6uTftEERo3ExkZwd//3pvU1P7U1tYREXGgYd+vv15i\n7VrbK1ImJKh44IHkBmHp3LkDI0ZchVrdHvPRNGbixD4EBJg/1snJka0eVVRwDkrXyYv49ts0Xntt\nF7t2VVNbC8nJKvr0CTQ7pkePUBYvvltZIVHB67DVthWh8TLS008wdeoKCgsNPPhgPK+9dp+nq6Sg\nYBeK0PgQVVVVnDx5Fp1OT6dOHUlMTPB0lRQU7EIRGgUFBZdjq20r01QVFBRcjiI0CgoKLkcRGgUF\nBZejCI2CgoLLUYRGQUHB5ShCo6Cg4HIUoVFQUHA5itAoKCi4HEVoFBQUXI4iNAoKCi7HI0Jjb0Lj\nto5yH0SU+yDSlu+DIjQeRLkPIsp9EGnL90HpOikoKLgcRWgUFBRcj9HNjB07VloCQPkoH+XThj5j\nx4612u7dno9GQUGh/aF0nRQUFFyOIjQKCgouRxEaBQUFl+N2ofn+++/p1asXKSkpvPrqq+6+vEdJ\nTk6mf//+DBo0qGFd8sLCQiZNmsQVV1zB5MmTKS4u9nAtnc/8+fOJiYlpWKcdbP/uV155hZSUFHr1\n6sWPP/7oiSq7hKbuw6JFi0hISGDQoEEMGjSI7777rmFfm7oP7hptMhqNRp1OZ+zevbsxIyPDWFtb\naxwwYIAxPT3dnVXwKMnJycbLly+bbXviiSeMr776qtFoNBoXL15sfOqppzxRNZeybds242+//Wbs\n27dvwzZrv/vo0aPGAQMGGGtra40ZGRnG7t27G/V6vUfq7Wyaug+LFi0yvv76642ObWv3wa0WzZ49\ne+jRowfJyclotVpuu+021q1b584qeByjxSDf+vXrmTt3LgBz585l7dq1nqiWSxk9ejQRERFm26z9\n7nXr1jF79my0Wi3Jycn06NGDPXv2uL3OrqCp+wCNnwloe/fBrUKTnZ1NYmJiw/eEhASys7PdWQWP\nIggCEydOZMiQIbz//vsA5ObmEhMjrpEdExNDbm6uJ6voNqz97pycHBISTGtZtYdnZMmSJQwYMIB7\n7rmnoQvZ1u6DW4VGEAR3Xs7r2LFjB7///jvfffcdb7/9Ntu3bzfbLwhCu7xHzf3utnxPHnzwQTIy\nMjhw4ABxcXH8+c9/tnqsL98HtwpNfHw8WVlZDd+zsrLMVLutExcXB0B0dDQzZsxgz549/799O0ZZ\nGAiiAPwarxAR7YKxkJjGG0jaNaWNnVcxt4iQTqt0Il4lYCyEpBEL0SLN/N0Wgn/l7kJ4X5UikJkh\nPNjdBP1+H03TAADquobneS5LtOZb35/vyO12w3A4dFKjDZ7n6aDdbDZ6edS1OVgNmvl8jrIscb1e\n0bYtDocDlFI2S3Dm/X7j+XwCAF6vF87nM8IwhFIKeZ4DAPI8R5IkLsu05lvfSins93u0bYuqqlCW\npT6h66K6rvV1URT6RKpzc7C9+3w8HiUIAvF9X9I0tf14Zy6Xi0RRJFEUyXQ61b3f73dZLBYyHo8l\njmN5PB6OK/291Wolg8FAer2ejEYj2e12//a93W7F932ZTCZyOp0cVv5bn3PIskzW67WEYSiz2UyW\ny6U0TaPv79Ic+K8TERnHL4OJyDgGDREZx6AhIuMYNERkHIOGiIxj0BCRcQwaIjLuD26iya7ZQTFk\nAAAAAElFTkSuQmCC\n", "text": [ "" ] } ], "prompt_number": 5 }, { "cell_type": "code", "collapsed": false, "input": [ "bmimodel.finalize()" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 6 }, { "cell_type": "code", "collapsed": false, "input": [], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 6 } ], "metadata": {} } ] }