{ "cells": [ { "cell_type": "code", "execution_count": 1, "metadata": { "collapsed": false }, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "/Users/baart_f/.virtualenvs/main/lib/python2.7/site-packages/matplotlib/__init__.py:872: UserWarning: axes.color_cycle is deprecated and replaced with axes.prop_cycle; please use the latter.\n", " warnings.warn(self.msg_depr % (key, alt_key))\n" ] } ], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt\n", "import matplotlib.cm\n", "%matplotlib inline\n", "from bmi.wrapper import BMIWrapper" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "collapsed": false }, "outputs": [], "source": [ "wrapper = BMIWrapper(engine=\"xbeach\", configfile='/Users/baart_f/models/xbeach/Zwin_T01/params.txt')\n", "wrapper.initialize()" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "collapsed": false }, "outputs": [], "source": [ "zs0 = wrapper.get_var('zs').copy()\n", "uu0 = wrapper.get_var('uu').copy()\n", "vv0 = wrapper.get_var('vv').copy()" ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "collapsed": false }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAW0AAAD3CAYAAADWiwWzAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztvXuQJFd95/s5J+vR3fPoeWlmJI0QAgkhJJDQYgwIjOwx\nYcAssBG7ZxFee0HXmL0Gr3bt5fLYdRA34t4wWl8HjzV7DQhYwRrDAe8aTAALNEhcEAIEkngJNCOk\n0cxopufdM91dzzzn/pFZ3dU9XZ3ZVdVVpzS/T0SGKquzM7+qyf71r7/5O7+f8t4jCIIgjAZ62AIE\nQRCE/EjQFgRBGCEkaAuCIIwQErQFQRBGCAnagiAII4QEbUEQhBGiMGwBgiAIo4gx5qPAq4Bpa+1z\nOhzzAeAVwBzwBmvtA71eVzJtQRCE7vg48DudvmiMeQXwdGvtVcCbgb/px0UlaAuCIHSBtfbbwOlV\nDnkN8In02O8Bk8aYXb1eV4K2IAjC+nApcLBt/3D6Xk9I0BYEQRgh5EGkIAgXFN/5lvXVxva8hx/Y\nu3fvU7u81GHgsrb9Pel7PSFBWxCEC4pqYzsvvfbNuY69+2cfujzjEJVuK/EF4C3AZ4wxLwDOWGun\ncwvtgARtQRAuOGLvej6HMeZTwM3AdmPM48C7gRLgrbUfttZ+yRjzSmPMfpKSvzf2fFEkaAuCcAHi\n6L0ltbX29TmOeWvPF1qGBG1BEC44Gj4etoSukaAtCMIFRz8y7WEhQVsQhAuOWIK2IAjC6CCZ9goY\nY14OvI9kAc9HrbW3r9e1BEEQ1kI8wrNx12VFpDFGA39N0kzlWuAWY8wz1+NagiAIa8Xl3EJkvTLt\n5wP7rLUHAIwxnyZpnvKLdbqeIAhCbsTTPp/ljVIOkQRyQRCEodMY3Zg9/AeRU1NTI/zxCYIwaPbu\n3dtp2Xhu4o4rz8NnvYL2YeApbfurNkr5g298ivKeq/H1iMq+/QCMX3UlgOzLvuzL/sL+HdfdQD9w\nI5wqKr8OT1GNMRHwS2AvcAT4PnCLtfah5cdOTU359qAdIpV9+xdumtAQbd0Tsr6QtcHw9N1x3Q09\nZ9pTU1P+4mf8Qa5jjzz8ib5k9v1kXapHrLUx8Fbgq8DPgE+vFLAFQRCGQYzKtYXIunna1tqvAFfn\nOXb8mitw59ZLSe+EnPGItu4JWV/I2iB8fVk4H2ZAzsPQH0QKgiAMmlCz6DwEEbQrDz1KeU+upHwo\nhOwvirbuCVlfyNogfH1ZNHx/np9lrfw2xmwH/jtwMRABf2Wt/W+9XFNmRAqCcMHRD08758rvtwIP\nWGtvAH4T+CtjTE/JchBBe/yaK4YtYVVCzihEW/eErC9kbRC+vixir3NtGSys/LbWNoDWyu92jgKb\n0tebgJPW2mYv2oOwRwRBEAaJ60++mmfl90eAKWPME8BG4F/2etEgMu3KQ48OW8KqtAr8Q0S0dU/I\n+kLWBuHry2KAJX/vBB601l4CPBf4oDFmYy8nDCJoC4IgDJI+2SN5Vn7fBHwWwFr7CPAo0FPH0yDs\nEanT7h7R1j0h6wtZG4SvLwuXM4vOCNs/AK40xlxOsvL7dcAty455CPht4DvGmF3AM4BfrU3tmjQJ\ngiA8+aj7Qq5tNTqt/DbGvNkY80fpYX8BPM8Y8yDwNeD/sNae6kV7EJm21Gl3j2jrnpD1hawNwteX\nRZ8eRK648tta+6G21yeAf9qXi6UEEbQFQRAGSSzL2HtDPO3uEW3dE7K+kLVB+PqyiEfYGQ4iaAuC\nIAwSl10ZEixBKJc67e4Rbd0Tsr6QtUH4+rKI0bm2EJFMWxCEC45+NYwaBkEEbfG0u0e0dU/I+kLW\nBuHryyLHwplgCSJoC4IgDJK8i2tCJIhfN+Jpd49o656Q9YWsDcLXl0WflrEPBcm0BUG44Aj1IWMe\nggja4ml3j2jrnpD1hawNwteXhcyIFARBGCFGOdMOQrl42t0j2ronZH0ha4Pw9WXR8FGuLUQk0xYE\n4YJjlFdEBhG0xdPungVtfrg6VqqgCvlzg7D1hawNwteXRZ+m0mROY0+PuRl4L1AEjltrf7OXa47u\nrxtBEIQucV7n2lYjzzR2Y8wk8EHgVdba64B/0av2IIK2eNrdE4w2v2wjIG0dCFlfyNogfH1ZDHAa\n++uBv7fWHoaF/to9EYQ9IgiCMEj6tCIyzzT2ZwBFY8w3Saaxf8Ba+8leLhpEpj1+zRXDlrAqIft3\noq17QtYXsjYIX18WA1wRWQBuBF4BvBz4c2NMTx+eZNqCIFxw9KmcL8809kPACWttFagaY74FXA90\n7S8FkWmLp909oq17QtYXsjYIX18WzqtcWwYL09iNMSWSaexfWHbM54EXG2MiY8wE8OskE9q7Joig\nLQiCMEgcOte2GnmmsVtrfwH8L+DHwL3Ah621P+9Fu/K+uwJfY8we4BPALsABH7HWfsAYsxX4DHA5\n8BhgrLUznc4zNTXl33zku7hzJXw9zBVII8Gw67SXM7qtHYSAueO6G9i7d29Pd9fU1JT/4vaP5Dr2\nVSff1PP1+k0vmXYT+FNr7bXAC4G3pDWK7wC+bq29GvgG8M7eZQqCIPSPPtkjQ6HroG2tPWqtfSB9\nPUvi0+whqVO8Mz3sTuC1WecST7sL0nroyr794WXZKUF+bm2ErC9kbRC+viz6sbhmWPRFlTHmqcAN\nJJ7NLmvtNCSBHdjZj2sIgiD0ixiVawuRnoO2MWYj8DngtjTjXp73rSkPrOzbv+S3eAj7Ietbkm17\nqOzfT2X//sRTVqvst863jvvjV1059M9ntf2Q9bXqoEPRE4q+ftF0Ua4tRLp+EAlgjCkAXwS+bK19\nf/reQ8DN1tppY8xu4JvW2ms6nUMeRK6RrF+JaXBe83n6TZhJijDi9OtB5N9N5luUeMvM7z+pHkQC\nfAz4eStgp3wBeEP6+l+T1CmuinjaOWjv6xGDbkA0D/Wf7qcwB1EVVAxoD+UYNjTRmxpEm+rojXXU\nhgZMNGE8xpdjfNnhiw5f8MkWpZtONwW+9Qugy1s2iM9tFULWF7I2CF9fFrFXubYQ6XpFpDHmJuD3\ngJ8YY+4nCSfvAm4HrDHmVuAAYPohVEjxoBzouiKqQGFeUZxVxCUWAi5lhxqLiQoxSnm8VzincE7j\nY42PFTQVxCop1vSgvDqv4VPreoLwZCPUh4x56DpoW2u/A3TyM357LeeSftqr4Nv+66BQ8Ywfiymd\nbBLNNZh0O4kP12lsLVIpgS7ETExUmNw0x46xWbYW59moqozRpEgS8J1TxE7RdJq6i6i5IhVXpOJK\nzMVl5pslKnGJarNIvRnRaGpcrHFOg1MQk6ThDhSqY2APvT9FyPpC1gbh68si1HK+PEjvkVHBJ/ZH\nYc6z8bEGGw5UUWfn8BHUL1VUiKhtUejIMbmhwmVbTnP1xmmuGDvBJYXTbNcVNmhHEY8DGigqXnPO\nlTgTj3Mq3sDJeCPH6ps4Vt/MydoGTlcnmK2V8bUijXoBGhG+oZMsvdmWnSskIxdGij51+RsKQQTt\nykOPUt5z9bBldKT9afnAWGZRRHXP2AnP+BN1iodnYHoWanVOqeNs2nE9Ok6/QXu0chR1zJius1HX\n2aKbbIs8E1pToGWLK5ootnrPzqjBvJ9n3jWZK1WYjWeSbLtZpN4s0Igj4ljj4palojjTLHC8UeCJ\n6mYOVbYyXy3TqBUS6yXNwIfyua2BkPWFrA3C15eFZNrC+pHG4qgOE8ccE4/ViA6dgeNnki9N1FD1\nJqrpAA3KE2lPQcWUVZMJ3WCzdkxqGNMFCis9e45iPBWgsvyydHoK+VijxE9rY3x/Zg+zpyIaM5tp\nzGh8XaPc6P5ACBcGoZbz5SGIoC2edhsrlPSpGPRcE33kDPrQWdRcdeHLW6PdOOdR3ieFHmmmHaWb\nxqPSkypAqZUDaucwu7LvsS1qcnW5ysbNx7iq2GRusky9WoBYg4ND9QkeuXwTx86e4fjM5sQPD4yQ\nM8WQtUH4+rIQe0RYPzzopkLPxegnThMdOrXs6x7lPMoByif2iPYUtE+CtvJ0iNM9sTWK2RrFPLN0\nHDYfP+/r98zv4KvndvPjQ3s4fW4D9fOy79H9oRFGH7FHekQ87c4oB4VK8gBSN8//+unmUba63SiX\n2CNKQaQcBRVTUk1KxBQAvR6RexUuK86z5af7+OfXzHFk7DC4iIKKuP/kbu49dgmNAP48DdmXDVkb\nhK8vi34F7TzT2NPjfg24B/iX1tr/0cs1gwjaAiuvdFRJ0I4qnuI8KwbtxD7xEPsk60497ZJOgnZR\nxUTKo1FJed6AuKw4zw3j87xg9yzVXQ1KqsC4KvKxh6/n/hO7aCzYJaOb8QijSz+Cdts09r3AE8AP\njDGfT3toLz/uPSR9tXsmCKNRZkR2wIOOPcU5R+msI6qf7y9vi3aD80ngJgnarSy7rJqMqZgiigg9\n8PB484smKCjNmCpSVElmfeP2o/zxs37Ei3cfIlLDrRMMOVMMWRuEry+LPrVmzTONHeBPSPozHeuH\n9iCC9gXN8hWIy76mmlCa85TOxujGCgf6NMt2PgnyylNQjpKKGVNNyiqmpDSR0gPNtFsUVMSYToK2\nUorn7jjGm695gJfsOsRY1KSgYgoqRreWZgrCAGh6nWvLYKVp7Je2H2CMuQR4rbX2/6VPf1YGEbSl\n90hnVOzRs3UKZ2uoRnze10/FR1HOoWKHamXaOqakm5R1TEk5IsVQAvZd98x3/No/2XGEt177Q267\n7j5uu+4+Xrjr8MAVhtw/I2RtEL6+LAY4BOF9wNvb9ns+qXjaAaM8qKYnmqsTna3i6ucH7cTTToJ2\ny9MuKkdJx4ypmLLyRAz+QWQWz91xjOfuWPxrUSnPD0/spuGitFFPWHqFJxd9Csh5prE/D/i0MUYB\nO4BXGGMa1trlA4BzE0TQljrtDniSRTNzFfzsHDTOfxK5Te9KrBHnEnsEKOiYchqwiwr0kOLfzS+a\nyH3sr190hOi6+/jO0T18e3rPOqpaJGRfNmRtEL6+LPoUtBemsQNHSKax39J+gLX2aa3XxpiPA//Y\nS8CGQOwRYQUW2rDG+Pk0aDc7lI84l2y0Mu2Yso4pKygxGv/IN+6Y5o+e+SAv2PkE46nXLQjrhfcq\n17YaeaaxL79sP7QHkWlLnfbKqLSzH7GDOE4eOi7jVDzNDjfZFrQdRZV42WWlKanCEOpGEu66Z35N\n2TbAr+98goJ2/H9H9/Cd6cvWSVlCyLXGIWuD8PVl0a8VkdbarwBXL3vvQx2OvbUf1wwiaF+Q5Pmd\n60G59ky6w0EL9ohH0+o7kgTt4pCqRrrlxh3T3LhjmrrT3H9yN41Y0/DDX4gjPLnIbY8E+KMTRNAW\nT3tllFeL7U87sE3tWgzqftEeKaqYgoJhhuy1ZtntvHDnE5S141tHL+M76+Rxh5wphqwNwteXRZy3\nF06A+UIQQVtYgZY1klm+nGbZcZKJa3xqj8RECrTqztHuNDu0U8OpfnPjjmmu336Mahzx4Mmd1F1E\nPYCl78KTgyy/OmSCeEYlddoro9LRYqvdXqfi6SX2iMJTSJevJwtWusPhqdGk6htUfYMG8ZqfoqxW\np50HhedFuw7zJ9fex/MuOtLTuVYi5FrjkLVB+PqyGGCddt+RTHvQ5Ix8CvJl2p4l9ohus0d0Fw+r\nWxl2jKPuY2KfBP4iEVppWivPB5Fxa5Vk3M/edpzT9THuGVApoPDkp8MfkiNBEEFbPO0V8IuZ9qqe\ntt6Jdw7flmknQbuJ6qK3h8MT46j5mIpzNPF4oIxD06SoopUHKaxAL572IAjZlw1ZG4SvLwvppy2s\nD2ngXv0Yn5QDxovVIyW6z7QdnoaPqbmYqnfUfXKDx8oT6aR2WiuFHmDGrfCMR022lKpU4gK1WG5b\noTfE0+4R8bRXRrl0WyX2nnLT4Bw+jpMgDxR6sEec99S9o+ah6jVVX2DeF5j1mnPOU/GOBg6X49y9\netottPLctOsQb732hzxvx9G+nBPC9mVD1gbh68tCPG2h/+S0R4DEGnEOn2baC562yt/Dw5N0CXRA\nw3tq3lPzEVUf0fARDRXjdTLWrOhjtAJFNJA/MrVKepU8a+sJZhslfnFmO3PNAtW4OICrC09G3AjP\nMQ0iaIun3YEcDyK3cVFa8hcnfbW9p4ijSIwiIvfqAN+a0u5p4qmTZNlVX6DmCxSI8U6hiYlweKXY\nqPWq5++3p11QnhfvOsimYo2pJ57a84PJkH3ZkLVB+PqyGGV7JIigLbThl77OKvkDkiy7rWFUUTmK\nyqHX4H4lrU48znsaXlH3mpovUPFFar5A5COcVihXTwK3doylDz5XGxjcTyLtee6OY1y95RRPzG+S\nahKha0K1PvIgnnYOhuHfLbFGVvO0fdretPVgkFam7df4j+txOGKgwWLQrrkiFVei4svMuzKzrsw5\nV+ScU5z1DSq+2VFevzzt9SJkXzZkbRC+viy8z7eFiGTag2KtN8AaPO3279FAgaSntlrDRZMFmIk1\n0vCKGlHqaRep+hLKe2KlQYPyDu3raNdAaSgpjUqnfAyqmmRTsc7OsTlmm0Xmm6V1v6bw5ELskR4R\nT/t8FKC8z7RHtqmdy75PUVCKolJrGnzg8cTe0fBQ9xE1V0gDdpGqK6JIgnbr14BGEeEp4hlTMSjO\nq99erzrtona8ZPdBtpSqfO3wU7vuBhiyLxuyNghfXxb9CtpZ09iNMa9ncXLNOeB/t9b+pJdrBmGP\nCB1wSeBeS5aufPKbuMDaGpR5SLNsFqyRqitScUUqvkjFl6j4EvOuzJwrM+tKnPMF5ryi6h0Nv/Zl\n7t1S0J4bth/jtU99mGdMnh7QVYUnEz7nthpt09h/B7gWuMUY88xlh/0K+A1r7fXA/wV8pFftQQRt\n8bRXoLWwJsMeWfC0W98DFJSioNbWRdvhaXpP3ZNWjaSB2pWoulISvNNt3pWZd0kAn/UFzjlP1adl\ng22Ip909IWuD8PVl4Z3KtWWQOY3dWnuvtXYm3b2XZYN/u6FneyT9bXMfcMha+2pjzFbgM8DlwGOA\naRMtrIWWp70GEttCE6HWNPzA+VZ9tqbqi8ynAbviS1TjIkpBrGJ8WvsdKUfkHGUVM6bqlL1nE6DT\npzcD6U2CZ0upyiUT5zhbLzMr3raQkz7ZIytNY3/+Ksf/IfDlXi/aj0z7NuDnbfvvAL5urb0a+Abw\nzqwTjF9zRR9krB/D66e9uHWi3dNu9d9OgvbaMu0YqHnFvI+YdWXm4jHmXZlKnGbZcSvTLiVZdtyy\nSkrMuSI1r2n4pSsl17v3SCly/MbFB/k319zPDdun1/z9IfuyIWuD8PVlMejqEWPMbwJvZOlk9q7o\nKWgbY/YArwTuaHv7NcCd6es7gdf2co0LmjzGWopKj1XpOsWkj3b+sB2jqPqIOVdk1pWZjcvMxWXm\n49KywF1asEbm09dzrkglDdpxD+1g10pRO56z7Tivesp+nrb5zMCuK4w+/ZgRSb5p7BhjngN8GHi1\ntbbnhzC9ZtrvBd7G0tCyy1o7DWCtPQrsXOkb2xFPexmtAOwXg3Enlnja6feq1BpZU6btddJjxJWY\njcvMujHm4hLzrkjVJUvGkweTyevFrLvMvC8z7yLmvaPuF3Nt8bS7J2RtEL6+TLzKt63OwjR2Y0yJ\nZBr7kknrxpinAH8P/L619pF+SO86aBtjfheYttY+wOopXaAl6iPAGjLtForE19ZKrclXbnrNvC9x\nzo1xLh5bkmVX42JSAhgXF4J3dWHRTZF5V6LiIyrO0fD+vAeSghAa/bBHck5j/3NgG/BfjTH3G2O+\n36v2Xh5E3gS82hjzSmAc2GSM+SRw1Bizy1o7bYzZDRxb9SzLaP0Gb3lmoez3rO/KdH9/jn0PG/dc\nCR7OHX+UauUYOxkDFjPrdi/7lD/GBFeAh+M/O8j9s2d48auSr7Wy3Za/3Gn/qudPcjYe4wffrXOy\nMcvm526nEpc4eN80sYvY9U8upqhjjvzwKEUdc9Xzt1DUMT+9d5ZNUYWtNxWoeM2P76mxScfsvWmC\nm180kfv6vezPNRqwae3/PuNXXRnM/SX7+fb7Rp/yiqxp7NbaNwFv6s/VElSnWYBrwRjzUuDP0uqR\n/wyctNbebox5O7DVWvuOTt87NTXl33zku7hzJXz9STwDcA2rGvEQVWHDERg/UGHsxwcpPn5q1W+r\nXr6Bk6/YwzNfXOXfXXsfL951nrW28uXSzPixxjg/rG7mkfokR2uTzDTHmY9L1F2EcxqtPAUdU4qa\njOkmG6IaGws1JgsVtkZzXFyY4ZLiGXZHip1RgeKAZkkCzNRLvO+nv8Yn9j17YNcUhsMd193A3r17\ne7q5pqam/B8eyZfw3nHx83u+Xr9Zjzrt9wAvM8b8Etib7q+KeNrns3CXZAT7pZ52/lasLVpDDyoe\nZuIxZpoTnGuOMdcsUWkWqcWFhaG6jXRL9gvpFlH3BeokLVzjtuuLp909IWuD8PVl0acHkUOhL8vY\nrbV3A3enr08Bv92P817w5HgQuYjq8Hp1HJ4GjopTnI3LnG0mDyArcZF6XCD2ydL1SLXqr5Oe3Q2d\nBvK033bTRzS9IhY7WxgFRvg+ld4jORhanXaOY5b0Hsm7/raNGE/VO+acYi4uMdcsU4mL1OIijTha\naGHpdVLKp5UnUo6G0zR1Eqxjr4m9xqFJKr4TZEZk94SsDcLXl02YWXQeggjaQgcyFtacdyxrvxWd\nh7r3VJxiLi4yH6eWSBzRdHrhT8SWjEhpCkoT6zRQLwTrUf4xEC44RjjTlt4jOQjZvzuvTnuNeKDp\nFXWnqaZDcxtpwI6dJl42M295Iq/waBwRngJuyVxK8bS7J2RtEL6+TPrRMWpISKZ9geNQxCiaXieW\nh0sCdnuQVoBuZdx+MZ9XyidVJcpRVDElFVNUknEL4ZOjGVSwBJFpS++R7lneT7sbFpOKpTeyUn7x\nHdX+nidSnoKKKakmZdVgTDfZqB0l5Ra+Rzzt7glZG4SvLxPJtIVRpZU3a+XRyqVbkkF7ryAN3Bq/\n8H7yMDKZ+F7WDcZ1nQ0qZoOCMbW25fOCMBQCLefLQxCZtnja3dOrp63xlJSnrB1jUUw5iilFyeDe\nSDuiVpDWjoJ2iRWiY0o6pqybTOg6m6MKm6Mmm3SBsooWgrZ42t0TsjYIX18W7R00V9tCJIigLXQg\nV8+aRbr5i04DJQXjyjEeNRgvNChFTUo6ToK0dhQWXifBurUqclw32KhrTOoqm3STca0pqmggvbTb\niZTn0olZnrXlBNvLYf+iEAJB7JHekDrtlclzzyx42mrZlhOtFGWvmNCeTVGdjYUqdRcBChV7nE/s\njoJ2SXYdNRiPGkwUamwsVJmMqmyPqmxUHr3sdhqUp12Omtx88QF2T8zy+ceu4htHnprr+0L2ZUPW\nBuHry2SE7ZEggrbQgTUFYd/h9epEKEoqYkJ7NhdqbGpWqUXFNFh74nTKelHHlKMmY1GSYW/QNTbr\nCpNRna3aMaGHVzVS1J4rJ8+wc3yeI/MbOVbdwJH5DZyshf0gVBgigWbReQjCHhFPuzNZ99ZST3vt\nf9NpFEUVMa5hc1RlslBhY6HGRKHOWKHBWNSgHDWTTTcZ0w0mojqbohqT0TxbdJNJXWBcRefdTIP2\ntMcLTX7z4gO86ZkP8OxtxzOPD9mXDVkbhK8vE5dzCxDJtAPFk/4Fp9cQhrtYEqlUMulmTDk26yqT\nUYn5qEzTJ+PKGirCAUXlGGtZI1GdjbrKpK6wUTcYU5piAH9tFrXjyskzbCnXuPfYJcOWI4RMn+wR\nY8zLgfeRJMAftdbevsIxHwBeAcwBb0hnEHRNEJm21GmvQBqAk7UsnW+w8+q0u7wXCypmU1RhMppn\nU6HKpqjGRFRnPKozpptLAvYGXWNTVGOzrjGumugOj9mlTrt7QtYG4evLoh/VI+lQ878Gfge4FrjF\nGPPMZce8Ani6tfYq4M3A3/SqPYigLXRgrZlzeqwn7ZO9hl7pReXYqOtsjqpsjipsLFTZENWZSIN1\nK2BP6BobdC05VjcZV/GSpeuCMBL0p3rk+cA+a+0Ba20D+DTJjNx2XgN8AsBa+z1g0hizqxfpQQRt\n8bSX0ZZl+4zAvbxO25O0W3VrDKQFHBMqKeHbGFXZGFXZ0JZtJwG7ngw/iKps0g02ac+Y7ryYRuq0\nuydkbRC+vgFxKXCwbf9Q+t5qxxxe4Zg1IZ52yKwh0/YLgd7jSKair2W8b6Q8Yypmo26wSdeo6QLN\nKEIrR8NHlFTMRJRm2arGBtVkg1aUlVrjCGFBGD6hLpzJQxBBW+q0VybPUOh2T9srj1eeGIfzDq/y\n/yEV4RlTjgnVZIOuUdEF6r6AUj4N2snqxw068bQntGNMFSjQ+feKeNrdE7I2CF9fJv15EHkYeErb\n/p70veXHXJZxzJoIImgLHVDg12hgeTxxuhXWYJEkKyNVGrgbbNB1ar6I8tDwEUXVZELXmFB1JlSD\nsnIUVEQ04NWPgtAX+lPO9wPgSmPM5cAR4HXALcuO+QLwFuAzxpgXAGestdO9XFQ87RwMxb9r97Pz\netoq8bRj74i9X5OrrZWiiGJMwXjaBGpC1xlXyTaR/ndc1xnXjbSb3+pXEE+7e0LWBuHry6If1SPW\n2hh4K/BV4GfAp621Dxlj3myM+aP0mC8Bjxpj9gMfAv64V+2SaYeMAq/Xlsl6BU2giV9TMqGAgtJJ\n8ygVM6YajKs6XkPBFyiqZhKwVZNx5SgpHYyT3XSKx2c388T8JiCZzn50fuOQVQlB0ydP21r7FeDq\nZe99aNn+W/tztYQggrZ42ufjSQK21/k9bWhNV/c0fFL2lzeyKlQ6SsxRImYsDdLeQYRLPG1VZ1w3\nGVeeItmnHpSnXY0L3HXkKXz54NMBaDjNkUp20A7Zlw1ZG4SvLxN5EClk0opweW+WbhpApfZIE0XD\nJ1Np8stTRGgKaavWkoopqyZOK7T3lFTSc2RcOcaVpqTUwL015+Gx2UkOzW5a8v58s8j3jl/Cj07u\nHrAiYVSR6pEeqTz0KOU9V2cfOCQq+/YPPrNIH0J6TaanvU3tXFLbHaNoZjrO510OhaKgoKigpBxl\n1cST9NwOnRSaAAAaA0lEQVQuq2ZqjcC4iiiq7HLCu+6Z72u23XSabx+9jH947Kol78dec7SyYc3n\nG8q/a05C1gbh68tEuvwJfaN9VWMrcGd+j0p6WKsku26gaKLXlmmnVSDaKwqQetsOTxNNUsM9rhxj\nSlFUmvPbQ/Ufn2bWB85NAont8b1jF/PgqZ4WlAmC2CO9Ip52B3LYI9v0LtAaoghUkl03vKbhu1tc\nroGCgpJPLBLwFHGMKc+E8pSVIso5UqzXLNsB90xfymcfTdo5eK+Y7iKj7kTImWLI2iB8fVmoQDv4\n5SGIoC2sQE57BEBFUbK1Mm0fUfcRrouorRQUvaaoHGU8mphYqcQW0WqhamQg02nSIP2TU70PLxaE\ndkbZ05Y67RwMrSZVZy+uOcXxJZm2Q1EnWc3ouvjn1ai09E9RVskqyXEVM65gQhUoqYi8xX4912mv\n8++FkGuNQ9YG4evLRMaNCevBgqe9avBSaaatU3uklWkXcDRY652nUBTQlNKmU4X02mWlKDIYLxvg\nkbNbeOTsFh49t2Ug1xMuMAINyHkIImiLp70CObv8bdM7Qask29Zppu0jGkQ4mmu+rEahVJTe00ng\nTqpKotxedotePO17j13Cpx55Ficq61frHbIvG7I2CF9fFqNsjwQRtC8octZr+9axWYmtagXsxPz2\nKJq+QMNHuC7KmlQamAtokjZSCZHSaUAfTKnU8eoEvzizYyDXEoRRoqegbYyZBO4AriN52H8r8DDw\nGeBy4DHAWGtnVjuP1Gl3QGfbI6fcNDv005Og3fYgsuEK+B5M4dbsyEXWfq5+12n3m5BrjUPWBuHr\ny2SEM+1eDcr3A1+y1l4DXA/8AngH8HVr7dXAN4B39niNC5OWNZJVPdKeaSuF94vVI7FXuDVOsFk8\nrUIr3bYNJsveP7OFLx+8gv0zW9f9WsKFi3L5thDpOtM2xmwGXmKtfQOAtbYJzBhjXgO8ND3sTuAu\nkkDekQvS01ZkWyTaZ5b8bYtaddrJgQuZNpomSS+S/KMQ+ks3Wfb3jl/CJ/ddx8na+DooWkrImWLI\n2iB8fZmMcKbdiz1yBXDCGPNxkiz7PuDfAbta/WKttUeNMVJk24mswN0+2FepZIngSge1HkSqZAFK\n00fUnabhHbGPUSpi3evn+sSp2jj7zm5lVPQKo8l6P4g0xmwlwyY2xuwhmR+5i8Re/oi19gNZ5+7F\nHikANwIftNbeSDIe/h2cH4YyPx6p0+5Ayx7RGnS04lT2U256oXKktSKyma6IrDtPwzvW1lm7f0g/\n7e4JWRuEry+T9a/TzmMTN4E/tdZeC7wQeMvyae4r0UvQPgQctNbel+7/PUkQn25NGzbG7AaOdfj+\nBWoHjiy8ruzbv+SGCGG/fvhwT9+/6v7+/VT2r7CfduybO/wYZ2vTqLEyRBGn/LElgw/OudOcqh9O\nA7pi5qePc/BHx5Ogjeeb91S4657KwvF33TO/JJiGsv/wzFa++PjT+ca3K1T2PbJ+n7fsj/R+v+jH\nEIQMXkNiD5P+97XLD7DWHrXWPpC+ngUeIsfQX9XNQ6oWxpi7gTdZax82xrwbaJmYp6y1txtj3g5s\ntdZ29LSnpqb8m498F3euhK9HnQ578tLp4/egYiifaLD93jNseGgGf24WX60uPa5cRu/YRvXKzRy/\naYzxZ9e4/pKDPHf7QZ43fogrSnNM6hJlFfZn+9/3Xcud+67jVG2cM/UyYo8IK3HHdTewd+/enm6O\nqakp/yff+HGuY//Lbz2nq+sZY05Za7d12l/h+KeSPP+7Lg3gHem1TvvfAn9rjCkCvwLeCESANcbc\nChwATI/XuHBR4CONHy/BxBi+Uln5GK3SCTcK75NWpQ2vqXtF0yd9qEONgQ/PbOUXZ7Zz77FLePTc\nlp7KFAUhL/2oDDHGfI3Ej144LUka9p9WOLxjdmyM2Qh8DrgtK2BDj0HbWvsg8GsrfOm313IeqdNe\nGQ/4SOHHSviJMsyc72adiqfZri5aXFzjFU2vabqIhtc0/dr6aveTPHXaPzh+Mf/t4WdzsjY+cJ0h\n1xqHrA3C15dJH242a+3LOn3NGDNtjNllrZ1ezSY2xhRIAvYnrbWfz3PdIBpGXdCstkxdAZEiHisQ\njxfx0Qr/XAu13GrBB3deLzyMXOswhEHx8MxWPn/gSr577FIOzE4yUx8j2D8HhCcf6/8g8gvAG9LX\n/xroFJA/BvzcWvv+vCcOYhn7BVmnnYe0YVQ8pmlOREQFzXJnelu0O11go9KYp4hb2TZJ8F7LMIR+\nslqWfd/x3Xz84edwojpBPKQpIiFniiFrg/D1ZTGA3iO3s4JNbIy5mKS071XGmJuA3wN+Yoy5n+TX\nxLvSYcEdCSJoC53xURq0xzV6pX8txZKA7T3ELAbsGIXz4LwfXB/sNvbPbOFnZ3bQcEt/3Xz32B4e\nn91M04f9gFR4krLOQdtae4oVbGJr7RHgVenr78B5eVgmQQRt8bQ747UiLkFzTFGIVqjTjo+yTe1K\nH0SCR+HSBTbN1NOO8TgY+MrIu+6Z58TFV/PRXz6Hc43Skq/NNko0s5qFrzMh+7Iha4Pw9WUhXf6E\n/tOaFakgLkFcVqyclKZZdrqwBg/OpxZJmnE7WgtsBm9DzDWKHJ3fyNlGeeDXFoSOjHDQDuJB5Pg1\nVwxbwqoMNaNQ4IqKuLzyFJvE027dg0lrVp/2IIlb9gg+DdqDvVND7vAHYfuyIWuD8PVlMYDFNeuG\nZNqh0KHPttfgiuBKasWgvdCbRC+OcU8ybU2MxvnE0/beZ0zA6S/7Zrbyk9MX8aOTu6m7IHIDQVgk\n0ICchyB+mqT3yCoocIUkcK9kj5xuHl1sKKVanrbGeb0QuON0As0g+fGpi3jP/9zKXU88hVocZm4Q\ncv+MkLVB+PqykExbWDc84AvgCn5lSzqtHmmtiISkXViSYSdbD50K1sy+ma38+NRFfOvoUzhRjVHN\nUvY3CcKgCTQg5yGIoC112m0sb9eqwEceV1jZ095a2L1kjqT3SXvWGJVm2SpdGj6Yu/THpy7iI7+4\ngSfmN6KfVgz6ZyNkXzZkbRC+vkxCvjEzCCJoC8toj7HpAhsKCsolVLmMbzTAucUDWvZI2oDbp8vZ\nFzJtWHgU2W/mnWfWw+PntvLIzE6+M30ZT8xvZE4ybCFgQrU+8iCedg6G7d95Bb6gUWNl1NgYFBZ/\n156Ojy5ObYeF5bcuzbQdGkfrd0D/79Q5D4ebnrtPbedv9j2Hbxy9nPm4CAz/c8siZH0ha4Pw9WWy\n/svY1w3JtEcBBb6oibdM0NwRo441UPX64tdbvUdSfOpjJ6V/63P/Ha5PsL+2mQOz4xw8N8bPTlzM\nsbmNVBrFpK1sn68nCP0k1PmPeQgiaIunnYECX4yIt22guRMKZ8+h0s9ra2E3sWorB1ywQ9RiwpCW\n/PUzkh6ob+TLZy/joemLOHR0G3PVMo1mYclvh6F/bhmErC9kbRC+vixG2R4JImgLq5D6274A9c0R\nhS0FVFnTFqNTK3txRWSSZZNWjvQn2646z6x3HKxsYd/sTn5y9iJ+NrODozObOVcZpxnrNMOWHFsY\nAUY4aIunnYOh+3cKXEFR36SpbdHEpcXAeLp5dLF6ZKG8LwnWi9l2748hK95zLHbce3aSOw5dwz88\nfhUPH9rN6dMbiZt6xR+CoX9uGYSsL2RtEL6+TMTTFvrOshWSSYtWaE4kJYDJMSrp2qc1vs3Tdj5Z\nxu68oomiSTJBdC334DmnORFHHK1t5FBlKyeqJU5WNPvO7uDY6Ukq82V8MwKnlsgVhFEghGnsbcdq\n4D7gkLX21VnnDiLTlt4j2SRB2xOPe3wBIAnWW8t7FlZDohZL/lxrNaRXNDzErC1on3Wa/Y0Sd53b\nwSemn8HHHr+ev/vV8/juoas4d2YDvholq3hWIYTPbTVC1heyNghfXybrn2nnmcbe4jbg53lPLJn2\nqKCSZexurIDfPonfqVDzNYh02k9bLdxoPi35SzLuZAxZJ+YdnHaaE41xphubOVObYK5W4kStxNFa\niYPzkxyd3cLZyjjNagHf1OBALW9kIqm2MEKo9V8m/BrgpenrO0mG9p434NwYswd4JfB/A3+a58RB\nZNriaefDR+DHC/iLtsDuHagN45xuHFnopb0YtFsPIJOGUdA5ps55xcFmxH3VTXzhzGXceexqPvj4\n9Xz8VzfwxUeezf2PX8HpE5tpnitBXafZdb4IHcrn1omQ9YWsDcLXl4Vy+bYe2GmtnQaw1h4FdnY4\n7r3A21hDXi+Z9qiQVpHEZcX87iLE45QLk/iTp/HjRRYGw/jF4b5VV2TGTXCoEXMq9jg3TjUe41xj\njEqjSL1e4Gwj4mQj4mhtA09Ut3KysoHZ6hiNegHfVBCrJVUhHcP18ltOMm8hZPqQaPc6jd0Y87vA\ntLX2AWPMzeT8qQkiaEuddk4UxCXF/O4CrqRRRExuejaNiUKSbbfKAz0LQft0PE7VRzRcgZONjUzX\nNnNobpKTcxuZnRujXi3gGhrf0BBrvFN9uaEhoM+tAyHrC1kbhK8vi348iOzDNPabgFcbY14JjAOb\njDGfsNb+wWrXDSJoC6uwQhWJK3vqk4q5PRH1rRGuoGluUEnf7Vhzdn4MH22hOl/igN5GwTviOGK+\nWWa2PsZMfYz5eol6rZiU68UqqQLxbV61rGoUnsysfzlfaxr77XSYxm6tfRfwLgBjzEuBP8sK2BBI\n0JYZkWtAgytBPYK4HFF9+BE27r4yCeYlTxxH1OfHOdcsc0RNglP4psI3NaqZBmen+vF0PJOgPrcV\nCFlfyNogfH1ZhDCNvdsTBxG0hRwsn2yj0ok2RU9zInnTRSRLJGs6CdQthy3NpH1bsFYoyaaFC5d1\nDtp5prEve/9u4O485w4iaIun3QU6ic/la64kbn/fgarpVUaLpY8UB7DaK8jPrY2Q9YWsDcLXl4X0\nHhGGw4qBeVmVxwjfnIKwXig3uj8YUqedg6BqUtXSrbK/g7YA7smgPrcVCFlfyNogfH2ZSO8RIUgC\nvekEYdiMcj/tIDJt6T3SPeNXXRns08SQPzcIW1/I2iB8fZlIpi0IgjA6XLAPIo0x7wT+FUkTuZ8A\nbwQ2kLMlYQup0+4e0dY9IesLWRuEry+T9W8YtW50bY8YYy4H3gQ811r7HJJfALewtpaEgiAIA2cA\nDaPWjV487bNAHdhgjCmQrJ0/TNKS8M70mDuB12adSDzt7hFt3ROyvpC1Qfj6slA+3xYiXQdta+1p\n4K+Ax0mC9Yy19uvArpwtCQVBEIaD9/m2AOna0zbGPA349yTe9QzwWWPM73H+M9fM//MzX7mHzTfd\nDCzWf7Z+k4ewXz98mMmbXxqMnvb9mbvupnTppYxfmX49rdse+v5VVy6p5Q3l82rfD1nfco3D1hOK\nvn4RahadB+W7/G1ijDHAy6y1b0r3fx94AfBbwM1tLQm/aa29ptN5pqam/B9841OU91yNr0edDhsq\nIT90WdAW2k2owv7cIGx9IWuD4em747ob2Lt3b09FrlNTU/7/fP+Pch377ttu7Pl6/aaX6pFfAn9u\njBkDasBe4AfALBktCZcjvUe6R7R1T8j6QtYG4evLYpQz7V487QeBTwA/BB4kWeLxYZJg/TJjzC9J\nAvl7+qBTEAShfzifbwuQnuq0rbV/CfzlsrdXbEm4GlKn3T2irXtC1heyNghfXxbrXc5njNlKjvUq\nxphJ4A7gOpIJrLdaa7+32rmDWMYuCIIwUNa/eiTvepX3A19Kn/tdDzyUdeIglrGLp909Q9e2yiOa\noWvLIGR9IWuD8PVlMQBP+zXAS9PXdwJ3kQTyBYwxm4GXWGvfAGCtbZKsf1mVIIK2IAjCQFn/oL2z\nfb2KMWal9SpXACeMMR8nybLvA26z1lZWO3EQQVs87e5Z0LZ8HFknBli8FPLnBmHrC1kbhK8vC9WH\nhTPGmK8Bu9pPS/IT+J9WOHylCxaAG4G3WGvvM8a8jyQbf/dq1w0iaAuCIAyUPjyItNa+rNPXjDHT\nxphdbetVjq1w2CHgoLX2vnT/c8Dbs64bxINI6T3SPedpUxnbMLUFRsj6QtYG4evLQnmfa+uBL5Cs\nV4EO61VS++SgMeYZ6Vt7gZ9nnVgybUEQLjzWvwb7dsAaY24FDgAGwBhzMfARa21rIvu/Bf7WGFME\nfkXS3npVggja4ml3j2jrnpD1hawNwteXxXpXj1hrV1yvYq09Aryqbf9B4NfWcu4ggrYgCMJACbSD\nXx6CCNpSp909oq17QtYXsjYIX18WoQ44yEMQQVsQBGGgjHCmHUT1SOWhR4ctYVXaewiHhmjrnpD1\nhawNwteXiUxjFwRBGB36sbhmWAQRtMXT7h7R1j0h6wtZG4SvL5NYgrYgCMLIMMqZtnjaOQjZvxNt\n3ROyvpC1Qfj6MrkQB/sKgiCMLIEG5DwEEbTF0+4e0dY9IesLWRuEry8TqdMWBEEYHcTT7hHxtLtH\ntHVPyPpC1gbh68tEPG1BEIQRwo2uPxJE0BZPu3tEW/eErC9kbRC+vkzCmcb+TuBfATHwE+CN1tr6\naucOwh4RBEEYJAMYgpA5jd0YcznwJuC51trnkCTRr8s6cRBBWzzt7hFt3ROyvpC1Qfj6Mll/T/s1\nJFPYSf/72hWOOQvUgQ3GmAIwATyRdeIggrYgCMJAcT7f1j1LprED501jt9aeBv4KeBw4DJyx1n49\n68TiaecgZP9OtHVPyPpC1gbh68skgGnsxpinAf+exPeeAT5njHm9tfZTq103iKAtCIIwUPoQtPsw\njf15wHfS0WQYY/4H8CJg1aAdhD0innb3iLbuCVlfyNogfH2ZxC7f1j2Z09iBXwIvMMaMGWMUyTT2\nh7JOHETQFgRBGCje5du653bgZcaYX5IE4/dAMo3dGPNFWBjq+wngh8CDJPbKh7NOHIQ9Ip5294i2\n7glZX8jaIHx9mazzasc1TGP/S+Av13LuzKBtjPloepHptJZw1cLxtFj8VqAJ3Gat/epaBAmCIKw7\nvVWGDJU89sjHgd9Z9t6KhePGmGcBBrgGeAXwX1OvZlXE0+4e0dY9IesLWRuEry+TEe49khm0rbXf\nBk4ve7tT4firgU9ba5vW2seAfcDz+yNVEAShTzyZg3YHOhWOXwocbDvucPreqoxfc0WXMgZDyP6d\naOuekPWFrA3C15dJHOfbAqRfDyLD/JUkCIKwEoFm0XnoNtOeNsbsAlhWOH4YuKztuD3pe6ty5iv3\nLLyu7Nu/xC8LYX/mrruD0tO+P3PX3UHpad9vvQ5FzyjpW65x2HpC0dc3RtgeUT6HMGPMU4F/tNY+\nO92/HThlrb3dGPN2YKu19h3pg8i/BX6dxBb5GnCVtbbjRaampsL8ZARBCJK9e/dmFjesxtTUlP9/\nXvfZXMf+h0//i56v12/ylPx9CrgZ2G6MeRx4N0mh+GeNMbcCB0gqRrDW/twYY4GfAw3gj1cL2ND7\nP4AgCMJa8b0tnBkqmUHbWvv6Dl86r3A8Pf4vgL/oRZQgCMK6MsJ12kGsiBQEQRgogfrVeZCgLQjC\nhUeg5Xx5kKAtCMIFh5fBvt1jjHk58D6S8sOPWmtvH6KWPSRdt3aRjP78iLX2A3mHdA5IowbuAw5Z\na18dmLZJ4A7gOpLP71bg4RD0rTRAFdgwLG0h9/TpoO0/A/8UqAGPkAygPTtobX1jhO2RobZmTQPQ\nX5P0NrkWuMUY88whSmoCf2qtvRZ4IfCWVE/mkM4BchtJdU6LkLS9H/iStfYa4HrgFyHo6zBA9ZYh\na1v3nj591vZV4Fpr7Q0k7SmGpa0/rPO4MWPMPzfG/NQYExtjblzluJcbY35hjHk4LZ/OZNj9tJ8P\n7LPWHrDWNoBPk/Q1GQrW2qPW2gfS17MkDcn3kG9I57qT/iXwSpJstkUo2jYDL7HWfhwg7T8zE4i+\n5QNUx0kWfQ1NW8g9fVbSZq39urW25SncS/JzMXBtfWP9+2n/BPhnwN2dDug2aR120F7eq+QQOXqV\nDIJ0QdENJDforqwhnQPivcDbWNo2IBRtVwAnjDEfN8b8yBjzYWPMRAj6VhigOpMOUB26tmX0tafP\nOnIr8KX0dWjacuGdz7V1i7X2l9bafSSDDTrRVdI67KAdJMaYjcDnSPy5Wc7vrTJwQ8wY87skHuMD\nrH4jDMusKwA3Ah+01t4IzJH8uR/CZ9c+QPUSkoz790LQlkFoejDG/EegYa39u2Fr6Yn1z7Tz0FXS\nOuwHkYeBp7Tt5+pVsp6kfz5/DviktbY11y3PkM715ibg1caYV5L8eb/JGPNJ4GgA2iC54Q5aa+9L\n9/+eJGiH8NktH6D6P0kGqIagrZ1Oerrq6dNvjDFvILHnfqvt7SC0rRXfh5K/Vaax/0dr7T/2fIEO\nDDto/wC4Mn1QdAR4HckDomHyMeDn1tr3t73XGtJ5O52HdK4r1tp3Ae8CMMa8FPgza+3vp0/1h6ot\n1TdtjDlojHmGtfZhkrl4P0u3Yev7JfDnxpgxkuqHvST33uyQtSmW/tXU6T77AvC3xpj3kmRiVwLf\nH6S2tMrrbcBvWGtryzQPWluvHHjH1/7N5XmP7fSF1aax56SrpDVXw6j1JL0Z3s9iyd97hqjlJuBb\nJA8RfLq9i+QmtCQZxQGSUqwzQ9TZCtqvNsZsC0WbMeZ6koekReBXJGV1UQj6jDFvIwmIMXA/8IfA\npmFpa+/pA0yT9PT5B+CzK+lJy+r+N5KePutd8reStncBJeBketi91to/HrS2UcMY803gP1hrf7jC\n1yKShGIvSdL6feAWa+2qE9mHHrQFQRCebBhjXgv8F2AHcAZ4wFr7CmPMxSTrP16VHrfmpFWCtiAI\nwggh1SOCIAgjhARtQRCEEUKCtiAIwgghQVsQBGGEkKAtCIIwQkjQFgRBGCEkaAuCIIwQErQFQRBG\niP8fu+yobMI4Z8AAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "from IPython.html.widgets import interact, interactive \n", "@interactive\n", "def plot(waterlevel=(0,10)):\n", " \n", " zs = wrapper.get_var('zs')\n", " zs[:] = zs0\n", " uu = wrapper.get_var('uu')\n", " uu[:] = uu0\n", " vv = wrapper.get_var('vv')\n", " vv[:] = vv0\n", " \n", " zs[40:50, 40:50] = waterlevel\n", " for i in range(100):\n", " wrapper.update(-1)\n", " uu = wrapper.get_var('uu')\n", " zs = wrapper.get_var('zs')\n", "\n", " plt.imshow(uu, cmap=matplotlib.cm.viridis, vmin=-1.0, vmax=1.0)\n", " plt.colorbar()\n", "plot" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": true }, "outputs": [], "source": [] } ], "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.10" } }, "nbformat": 4, "nbformat_minor": 0 }