{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "Notebook to quickly examine one simulation" ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "collapsed": false }, "outputs": [], "source": [ "import netCDF4 as nc\n", "from salishsea_tools import nc_tools\n", "import matplotlib.pyplot as plt\n", "import numpy as np\n", "import os\n", "\n", "%matplotlib inline" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "/data/nsoontie/MEOPAR/SalishSea/results/2Ddomain/3.6/timestep_check/diff1e-6_river\n" ] } ], "source": [ "results = '/data/nsoontie/MEOPAR/SalishSea/results/2Ddomain/3.6/'\n", "\n", "run = 'timestep_check/diff1e-6_river'\n", "\n", "path = os.path.join(results,run)\n", "\n", "print(path)\n", "\n", "filepart = 'output.abort_0000.nc'" ] }, { "cell_type": "code", "execution_count": 17, "metadata": { "collapsed": false }, "outputs": [], "source": [ "y = 3\n", "f = nc.Dataset(os.path.join(path,filepart))\n", "Us = f.variables['vozocrtx'][:,:,y,:]\n", "Ws = f.variables['vovecrtz'][:,:,y,:]\n", "Ts = f.variables['votemper'][:,:,y,:]\n", "Ss = f.variables['vosaline'][:,:,y,:]\n", "ssh = f.variables['sossheig'][:,y,:]\n", " \n", "depths = f.variables['deptht'][:]\n", "x = f.variables['nav_lon'][y,:]\n", "\n", "xx,zz = np.meshgrid(x,-depths[:])" ] }, { "cell_type": "code", "execution_count": 13, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "(-50, 0)" ] }, "execution_count": 13, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAwcAAAEKCAYAAABHWSm4AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xu8JFV56P3fwwwDjKDc4nBVUERFo+hJ8P4yJop4iWiO\nNxKNQY+5GC+vmhjRnDhoAppE4jmJ+uZ4iRwUEo4XovE6EneiiaLkiIIjKoSJjMAgIDeBuWye949a\ntXftmurevXv3nu7Z8/t+Pv2p7qpVa61a1TNda9d6akVmIkmSJEl7jLsCkiRJkiaDnQNJkiRJgJ0D\nSZIkSYWdA0mSJEmAnQNJkiRJhZ0DSZIkSYCdA0kdImJjRPzyuOuxWBHx2Yh4ybjrsTNFxIcj4u0j\nzvOoiLgnIob6zSjfp18q798cEe9vbHtuRFwTEbdHxCMj4sERcWlE3BYRrxrVMYxaRLwvIv5o3PWQ\npFGzcyAtUxHxxIj4t4i4JSJuioivRsQvDLh7ltcg5Qx9MRoRB0bEJyPijnIBeeow+ZS81kXEuc11\nmfmMzDy31z5LoXkhPCYDn7tRiohTI+KjPTbP1Cczz8zMVzS2/QXwyszcLzO/DbwRuCgz752Zfx0R\nT46IL5fv8dUjqOfqiHhvRPyk5PnPw+STmb+bmX+y2PpI0qRZOe4KSBq9iLg38I/AbwMXAHsBTwK2\njLNeHd4D3A3cF3gU8JmI+HZmbhhvtRYlgRh3JcbgmcBnFrJDRARwP6B5vu8P/Fvj8x3AB4DVwJsX\nWUeA/0X1h7GHADcDxy80g4jYIzPvGUFd2vkGQDo7qaQx8s6BtDwdS3WN8fdZuTsz12fmZQAR8cCI\n+KeIuLH8BfUjEXGfhRYSEb8F/BrwxjIs5B8WsO+9gF8F/ntm3pmZ/wr8A7DgYUARcTJwOvDCUo9v\nlfVTEfHy8v43I+JfI+LsiPhpRFwZEY+PiNMi4kcRsTkifqOR514R8RcR8Z8RcX0ZRrJ32XZwRPxj\nyeemiPiXqJxLdbH76VKP3y/pH1vu4vy0DJk5sVHOVEScFREXR8StEXFhRBwwwDE/sZHnj5p1Bw4s\n9bstIr4eEQ9o7PeQiFhf6n1FRDy/sW2fiHhXuftxS0R8JSL26ij7v0bE1RFxXPm8B/AU4PPl80tK\nu90YEW9u7bsuIs6NiFXA7cAK4NvlfFwErAX+utT9mMz8ZmZ+FOi8a9DveLrSAr8C/FZm3lT+bXyr\nf0vP3B17X1TD1O4AnhyNO2YR8b2IeGYj/cry7+r48nm+8/8nEfGvwM+Ao+erjyQtJTsH0vL0fWC6\nXMCc3ONi80+BQ4GHAkcC67oyKhehP+3alpn/C/go8M4yLOSUsk994dz1+lTZ/Vhge2Ze2cjy28DD\netTjfmX/Izrq8XngTODvSj0eVW9i7hCbE0oZBwLnU91VeTTwQODFVBelq0vadwDHAI8sy8OBPy7b\n3gBcAxxMddfj9HKh+RLgR8CzSj3+IiIOp7qL87bMPAD4feDjEXFQo14vAU6jOh/bgf/Z1QaNtrg/\n8Fngf5Q6HF+OC6q7Fi+iOp8HAFdSneu6Q7Ye+AjwcyXdeyPioWXfv6C6g/O40kZ/0Gq/iIjTStv8\ncuMOzwnAf2TmzaXD8F7g14HDgIOA5jlLgMzcmpn7lnWPyMxjMvOXga8Av1eGFTW/G13tMN/xtJ0A\n/CfwtnLx/p2I+NV+ZTScCry91PmrzP1unVe2154G3JCZlw54/l8M/DdgX6rvjySNjZ0DaRnKzNuB\nJ1JdvLwfuCEi/iEi7lu2X5WZF2Xmtsy8EfhL4MQeeX21XNT0M2cYTWY+KzMP6PF6dkm2L3BbK5/b\ngf161ONHZf9Nfeow33CeqzPznDJs4wKqi9e3lXZYD2wFjomIAF4BvD4zb8nMO4CzqC4+KekOBY7K\nzOly16OXFwOfLR0YMvNLwCVUw3CgOkf/OzM3ZOadwH8HXlDq0MuvAevLnaHpzLy5jNev8/tEZl6S\nmdNUnbd66MyzGm1wT2ZeCnwCeH756/9pwGsz87qy/euZubVR7uuoLm5PzMz/aKxvDil6HvDp8r3Z\nWo6nOQRnkCFXgw7L6nk8PdIfATwcuIXq/L0KOKfcUegngQsz82sAmVkPz6vreT7w7PrOEtX5Ob+8\nH+T8fzgzv1eOYfsgBy5JS8XOgbRMZeYVmXlaZh5JdUF0GPBugIhYExF/FxGbIuJW4Fyqv/DuTHcA\n926tuw9VB2GpbG68vwsgM3/SWrcv1V+hVwP/Xt/xAD5H9Vd6gD+n+ov8FyPiqoj4wz5l3p/q4vun\njbyeABzSSHNN4/2PgD0bZXU5AviPPtvbx1n/hf7+wGNadfk1YA3V+d8buKpPvm8A3pOZ17bWP53q\nTgZUF90zHbjS4bmpT55dBh1z3/N4IuLIqALdb4+IuhN6F7AN+JPM3J6Z/wJ8GThpgLKu6bWh3OH4\nHlUHYTXV0KXzGnVcyPmXpLEyIFnaDWTm9yPiHOC3yqozgWng4Zl5S0Q8B/irYbNvr4iIz1Hduejy\nL5n5TOAHwMoyrrwePvJI4PIh6zHKANEbqS4kj8vM69oby52E3wd+PyIeBvxTRHwjM7/Mju3xI+Dc\nzPytdj4N92u931bq0Ms1VENkFupHwD9n5g4Xw+XOwd1UQ6i+02P/k4AvRMT1mfmJst8hwKGNsfvX\nUQ1Vq/NdzdJ1PHseT7Fv63N9XO07E6MIAD6famjRCmBD487KIOffAGRJE8M7B9IyFNWz4l9fxjsT\nEUdSXbh8rSTZlyr48baS5g8WUdxm4AHNFZn59DLmvuv1zJLmZ1RDQN4W1eMln0j1F9dhHz26GThq\nnuE4A8nqSTTvB94dET8HEBGHR8RJ5f0zI6IefnQbVUer7pxspophqH0E+JWIOCkiVkTE3hGxtj43\nVBeqL46Ih5YL6bcB/6cMferlo8BTIuL5Jfj1oIh4ZCO/Xj4DHBsRL46IPcvrFyPiIeWYPwScHRGH\nlro+LqrA4dp3gZOB90TEr5R1T6e6q1L7GPCsiHhC2fdtLPy3ZuYYorI31d2UiCpQvK7TP/Y6nh75\n/jPVxfrppd2eQBUA/YVB69Nn3d9RxRr8DtX5qc13/nvlL0ljYedAWp5uBx4DXBzV01W+RvVX0zeU\n7WdQBeLeCnwa+Dg9/noZEU+KiH5DfT4IHFeGTHxigfV8JbAPcAPVRdTvZOb3etTjfmWIyA4BycX/\nKcubIuKSju1dz//vdwH+h1RDh75ehl6tpwqiBnhQ+Xw71WM335OZ9fPyzwL+qLTH60uMxClUj+G8\ngeri9A3MXhAmVYfow1R/dV8FvKZPvcjMa4BnlHxuAr4FPGK+4yyxKCdRxU78uJR3VikTqrshlwHf\nLPme1aonmfkdqrH+74/qKVHPZHZIESVI+feohtVcS/W40OawmXb9us5Bc92JwJ1UHZsjqe7o1OP3\n75jneOZmWo3nP4Wq7W4B/gZ4SWb+oCt9nzrvsC4zr6f6LjwO+PvG+vnOf/t4JWmsov8fpxaRcfWj\n8W6qW6wfyMx3LklBkrQLi4gvUw07+dC467JQEbGS6oL86HKhLknaxS3JnYOIWAH8NdXt5+OAU/s8\nWk6Sdne76rCSA4A/smMgScvHUg0rOgG4MjM3ZuY2qrGYpyxRWZK0q+sK6v71Moyq/bpsHBXskpk/\nycy/GXc9RiEivtujvU+df29JWj6W6mlFhzN3jOkmqvHPkqSGzHxyj/UfZW5gq5ZQZnZOvidJu5ul\nunNgcJUkSZK0i1mqOwc/pnqqRO1IGpPiAESEHQhJkiTtFJk58fFdC70+XopjWqrOwSXAgyLiKKpH\n2b2Q6hnrLe9douKX2p4DpKmbdp8F5Lt6gPwW6iPAi4fcd1S2D5CmfXxdbTxMmkHav12/hbb13wCv\nanzeVpZ1/Zp1qPPe3uNzVz068jm4NbHwlrKsp3yaM23XbTPvMtsTEu8a1q1bx7p168ZdjWXPdl56\ntvHSs42X3q7YxiOYAmen+ZMB0/3REpW/JJ2DzNweEa+imlhmBfDBXs8ulyRJklQZ5E/QS2mp7hyQ\nmZ9j7qyZkiRJkvpY7MV5mVLgEmBTZv5KRBxINTnj/YGNwAsy85Ze+ztD8rL3iPmTaJFOGHcFlr21\na9eOuwq7Bdt56dnGS882Xnq28dLaZ8BXH68FNjD7gKA3Aesz81jgovK5pyWbIXk+VcCFMQdzLUXM\nwSRY7jEH7XoYcyBJ0iSJiF0mIPn9A6Z9BTsGJEfEEcCHgT8FXl/uHFwBnJiZmyPiEGAqMx/SK99d\n+YpzjLa1PnddpNYXe7e31jebvH3hemdHPnXe7QvOST11g3QEBtE+3mHV7d/VXqMa1XdXnzKgeXHe\n3WFo79tuw7otFtLRlCRJu6JFXuH9JfAHQPOvgWsyc3N5vxlYs4TlS5IkSRqVXn+63AD0e7pPRDwL\nuCEzvxURa7vSZGbO97hUOweSJEnShOh1cf4I5kaSfnLHJI8Hnh0RzwD2Bu4dEecCmyPikMy8PiIO\nBW7oV74ByZIkSdKE2HPAV1tmvjkzj8zMo4EXAf+UmS8BPgW8tCR7KXBhv/K9czAS7RgEGGw8+12t\nz/1iF+pT1TUWvx2XUI9N7zf+f5BTP+on7baPZdj9e31uapfRL4ah33H2ivHoV/YgsQF1fephfzd1\n7F8PF6yHCc7GLjz8J98H4PLX/2K14sNlw8PK8gON7L5kELIkSbuKEV591cOH3gFcEBEvpzzKtN9O\ndg4kSZKkCTGKx49k5j8D/1ze3ww8ZdB97RxIkiRJE2LcF+fjLl+SJElSMepB3Qtl50CSJEmaEOO+\nOB9z+c2+0UKrUgeE9pscrN33as9Ge1fHtkH6a/VosK5A5F66gnHbgcTttF32bC27ttW6jqm9rl+g\nbr9z0m7vdptu65P23q3Pzf23tT43R97V7/sFB7ePp6s+tboeB7bWN8sc5PvQqnO0RgvOeZrwzXPL\n/IXy8dIjZpPsV5YPLMtLjmrtO4BnnjH7/nVvHXw/SZI0Vt45kCRJkgSM/+J83OVLkiRJKrxzIEmS\nJAkYzaNMF2PMnYNjmO0fNZui1/j/ZoxAe2x6bb/G+zrP9vjz9uRjTfu0ll1p+8UctOvTbuKuU17v\nc1vHtvZ4/a797+qRth5T38z3oFaa2zvybccRDBKPUbdFs726zi2zp2ivxrotrd23t9ICHFyWd/co\nGuBnrXX9wjcOKMt6bP8tHdU9uLVt77JsntbrynHuU5Z1+EBdh+lG2isPnJvmRWV5fCPNoWVZxyP8\nVVn++2xsxAs5E4Cr1h0DwF2/Ux3MR459HgC3NxruWt7SyPxPkSRJk8s7B1qmxv3VliRJ2vWM++J8\n3OVLkiRJKvYc9Oq83+iIRbBzIEmSJE2IlXYOJEmSJAHsuWK85Y+3c7ByLexf3jeHqO/dSlcHoF7X\nXFlPCHVTWa6pFgffezZJHdhZB7leVZb5w1YeAEdVi/3WzMmOjY0k7Xm56no2A2TbwbJ1mv1bS5gN\nWK2Dcevj6+oJ1sGz+5ZlM+63Dnjdu7WsY49vabTJka18Nh64Y73qbbeUk1IH5R7cSNOOY95U0t5Y\nls1vVl2POkb2qLJsfvnbbVCf1kMaaY6p68XcMm5vpKm31fndwVzNtq3rsbYsL23lC/DwstxYlvcq\ny+b39cqyrNunDjb+aVk2A5LrtnhQWZ5clg9ppCnHfP9HXAHAfz7kIXPLAZ7GFwD48L1PA+Cq6w+g\n6XfjHTPv35VvQ5Ik7RoGvnOwVOWPt3hJkiRJtT33mj/NUrJzIEmSJE0K7xxIkiRJAsZ+dT7e4i9i\ndpz2vrODwVfteycA+5Tl7T+tBqvf8/17ze57SRkrv6ksy9jwfZ7305kkR917IwC3lAH11118dLXh\n78uA76lGXZ5Xli9KAA466loAbvrm4bNpflyWdVxB16Rc7XnD6vH7x1Qb9j34lpmkd95eDdi/59Zy\nXBv75HdUWdbxBDfSW12vvatj4Y6Y2bTHmmqWsP0OqAbq33pjNVZ91d5bZtLU7X7r5WXw+yHVtgce\nNjvofUuZwWyvMrj/mpurYIatVzTiG+q6Hjz32I9cfc0OVd7Kqu58DsmZNA98wAYArr3tsDl1vvX6\nRjDEjaXh6/Zpx4A0v/ElhmHtsZ8H4N9ufjwA27fNBkMct6Yqc+Od1Xdnr1Lmij1mAwlu+M79qvoc\nUU02d9CBVcDE5s33BeCen62eLXP/mFP2w4/9JgBXHXLMTJKD7l3tf3T5Qqw+tgow2Xrsqpk0Gzl6\nzmEdf+zXAbgyPl4d9s8+PrPtC5yEJEnaRYy5c7DHeIvXstUOKpckSdL8Vgz4aomIvSPi4oi4NCI2\nRMRZZf26iNgUEd8qr5N33HuWw4okSZKkSTHk1Xlm3h0RT87MOyNiJfDViHgikMDZmXn2EhYvSZIk\naeQW8bSizLyzvF1FdX+hHm8f3XvsyGFFkiRJ0qRYOeCrQ0TsERGXApuBL2fmd8umV0fEtyPigxGx\nf/feJY/M7Ld9yUREnpX/L/uVGaz2aczqVQe5riizR02XgVXXcthMmo0lQvfGMqvUkVRBrj/PZTNp\n1nADAHeWmbp+wIMB+DeqwNNL+IWZtE/hSwA8hosBuG/Z9wbuO5Nmc5kZ7fYym9f+Zcatu2ZmApsN\n1F3NnXOOYQ2bAdi3MWPX1pK2Dsa9sUQSN/OrHTQzKxjlmGaDXOv2qet1UInG3YutpU6zgayrSzuv\nKJHO0+XbVde36UoeCMB+ZSaxuo2b6nNVt821Zea51Y3zWeddt1d9zrd3DJi7oeRzTZmtbb9Gex1V\nAnRvKue8bttbGjO41e1yC93f+xWNGcn2L53p+jtzGT8PwF2Ntq2/B/Vx1e21ojGbWh0cXJ/jfcrx\n3lTO59ZG+//L1icBcNiqara3Z/AZADZw3EyaA0o71cfb/N7XnsFn55T5Wv4HAM+JxwHw4HzOTNoH\n8/2Z949iww55SZK03EUEmTnwX8/HJSIyHztg2q/T85gi4j7AF4A3ARuAn5RNbwcOzcyX98rXYUWS\nJEnSpOgINgaYurV6DSIzb42IzwC/kJlT9fqI+ADw6X772jmQJEmSJkWPq/O1B1Wv2hmtAR0RcTCw\nPTNviYh9gKcCZ0TEIZl5fUn2XGgMsxm8eEmSJEk73fBX54cC50TEHlRxxedm5kUR8b8j4niqpxZd\nDfx2v0zGGnPwjXz4zLj4pnpceHtcfHPsdj2mvB7jX49NP6gxO9jKks/21pj868oY7npcO8CxZVx2\nHadQ16FZZj2evV5Xp2mO6a/r2o6bqMehN4+3uV+V716lnNmYg/oYVrXaabrjnlNdv3bZTe11dT77\ndMQc1G28cuYYZuMImmPum3Wv27iZZuXM+Zy77FK3bTOmolaf4/a25jna0lEP6G6vuh4Hl+9MHTfR\nzK8+5l4xDM1tdf1WlfavYxeasRWX8ihgNlagjndoxhXU+9ffxTr/Zrv93p3vBeDi1ScAcHWJwXl2\nuVO4+s7Zc3Xn6tnv0wEd51mSpOVul4o5eNqAab/QO+ZgMbxzoCXR7jxIkiRpAIt4lOko2DmQJEmS\nJsWYr86HnucgIp4fEd+NiOmIeHRr2+kR8cOIuCIiTlp8NSVJkqTdwIoBX0tkMX2Ty6ginv+muTIi\njgNeCBwHHA58KSKOzcx7FlGWJEmStPyN+c7B0MVn5hVQBXi0nAKcn5nbgI0RcSVwAvD1dsI64BL6\nB6nWmkGldfBnHYBaB+E2A3fbk6jVwcD18r6N8usJsept7X2bZbYn75ru04y9gnIBVrfyqfPdvyO/\nQcbw10HBvQK6u9bVx9fV/itb6/qdox2Djrf32bZj27br0VVWOyi7Vp97mA2sbk6e1pV/s4715Gz1\nvs1J0PaZmcDtpyWfHSdBq4Of64nf6vrU383mcdZtUdevnuisGRBefwfretVB9s3zcc3qI+YcXz15\n3MyxNb4uB1wzG5zciMGXJEmTaFcdVtTHYcCmxudNVHcQJEmSJPUzycOKImI9cEjHpjdnZt/Z1Vo6\nn5f6l+vumHn/+LUrefzaPReQpSRJkrSjqakppqamxl2N4UzysKLMfOoQef6YuYMXjijrdvC6dfvO\nvB9kWJEkSZI0n7Vr17J27dqZz2ecccb4KrNQe4+3+FH1TZqBB58CzouIs6mGEz0I+EbXTvtP38KK\n7b3H0q+crmKYt6+oRj9Nr5yt7vSK6n7KDmPVpxtjykve9X5bV1STW9Vj15vj0tuTlNX51OU0dU2o\nVavjBhYyXr+dXzumYVDt2Id+MQe99mmm2dp60G7XeP12Pu0J4rr276pXrd6/bpNmW9Rj+re2Jo/b\n0vFA4HaaXucFZr8H9bYtzN7RquMHVjcmgOu1f52mawK9Wr1tdfm+1ZOhNeMm6jR1zEF7cjuAI++s\nRu5du7qaPO1pfGFOOSuvbXy4V8+qS5KkSbOEQ4YGMXTnICKeC/xP4GDgMxHxrcx8emZuiIgLgA3A\nduCVOa5pmCVJkqRdySQPK+onMz8JfLLHtjOBM4fNW5IkSdot7aqdA0mSJEkjtqsOKxqF1T/rfm59\n24rtZf60LbPpp0vN63iEOj5hRUcIw/R0iSMoG1etbORTYgrqGIM6TqHOr6kd+9COaQB2GGXeL6Zi\nh3qWfHYcqT7g/q1j6dq2w/o5cQRz99urNa9AV77t/Jvj4tv59tu/nc+WjlaoYwJWtcromv+gHcfR\n/Nw+zjrepLm+Lr8us07TFSdRxw+smpkjo3esx4rWMRzETTPLer6EOu1+JfahjmnoOs7H31aF82y4\n97FzjvPGh+zLwVdU+8f91s2kd4CfJEkTzjsH49Prgnm5GeQ4d5enRQ1ynF0dk8WVOX8Hse4YjErd\nMZAkSbsYOweSJEmSADoewrhT2TmQJEmSJoV3DiRJkiQBY786H2vxXcHDC8+jdyByu5wVK6u0dYBy\nHWAMsNeWezr3mV7ZXDd3/9rK6cECq+ezsqNeO6bpHSjdK/h5IUHRTe3A6/nq0axL1/b2OZru+PbN\nHMvKjknUOgLAAfZasWP7twOK+02Ctqq0+8oVdbDw7P28fhO21erJztoByf3UE6btP11NdLZqxWyQ\n9UwQ9JYSBN1xHvb6WbXc3rr1ePA1JQj5oevmrYMkSZpAYw6J7X0VKkmSJGnnWjngqyUi9o6IiyPi\n0ojYEBFnlfUHRsT6iPhBRHwxIvbvV7ydA0mSJGlSDNk5yMy7gSdn5vHAI4AnR8QTgTcB6zPzWOCi\n8rknOweSJEnSpFgx4KtDZt5Z3q4qqX4KPBs4p6w/B3hOv+INSJYkSZImxd7D7xoRewD/F3gg8L7M\n/G5ErMnMzSXJZmBNvzzG2jmIaZiJ3WzWpGtdS0e86rzqTlY7QLm5bk7d6O6Y9Qt+7gqy7bVP77Td\nwb699ukVHDxoPXqlmS4zSS8kcLxu04WU07V/O+h7TvoyU3Z7luy5+cyd8bpfcHUdjL5irx1n0K71\nm0iunhW6Dmxul9E8PytXT89Ju9+tZZ97NYKNS31WlhjlLAHX0fzO13OclYDkR1zxgyqNgciSJO3a\nFnF1npn3AMdHxH2AL0TEk1vbMyJyiYqXJEmSNFI9/h45dXn1GkRm3hoRnwH+C7A5Ig7JzOsj4lDg\nhn772jmQJEmSJkWPq/O1x1ev2hkXzN0eEQcD2zPzlojYB3gqcAbwKeClwDvL8sIhipckSZK00w1/\ndX4ocE6JO9gDODczL4qIbwEXRMTLgY3AC5am+FHY3uN9v3W1XjXvs0+0Pjfv2kSPGIZe63tZyLwV\nvcbl94pFaO4zyJj+hda9bZi4jr77DFDnev/Z8faz29rHU6dZtWLHmIvpdtzAlt4xDDPxJfUkdx0x\nB/3U52vV3d31nFvmHXPLvrVa3qsZN1K3U53fytZ6oJ7LbeW15fNNC6qyJEmaVENOgpaZlwGP7lh/\nM/CUQfPxzoEkSZI0KRbxtKJRsHMgSZIkTYoh7xyMip0DSZIkaVKM+erczoEkSZI0KewcDKlXcOsC\nJuyaE6Dca78FttB8QcA5wK2irmDjHfJdwHGOzCBldgXP9tq/q21bE+DFfEHrNAJ2m1nX7dWeUK8r\njzrwd+/WvgPuv7Lednf/egKsbAUbl/nT5mrXfUvrc1MJSI6nretdqCRJ2nXYOZAkSZIEGHMgSZIk\nqfDOgSRJkiQA9hpv8XYOJEmSpEmxW985WOqg2l7BpP2CU3vl0dxvIWW3xKiOeWcGJHeV1Svwevs8\n2wfJY9gymxZyjruCgufTrMuK1rp2mR0zGw8UbNyvfe4eII0kSdr17NadA0mSJEmz7BxIkiRJgsEe\ne7+U7BxIkiRJE2LaOweSJEmSwM5Bt4XMxDvIPv0CkXu1wEKDlReStleA9KjKqo0qWHWxgduLDZ7u\nFVzcL+i41wzJXem7goIXEtDcLrNXOV31Wmi71etuqxbxjHUDV0+SJE2+LXutGjDl1iUpfzI7BzvL\nOI5+ubf4qI9v1PkN0/HcGZaqgyhJknYp0yvGG3Sw3C9VJUmSpF3GNHYOJEmSJAHbx9w52GOspUuS\nJEmaMc3KgV5tEXFkRHw5Ir4bEZdHxGvK+nURsSkivlVeJ/crf+g7BxHx58CzqKIhrgJOy8xby7bT\ngZdRhV++JjO/OGw5kiRJ0u5iEcOKtgGvy8xLI2Jf4N8jYj2QwNmZefYgmSzmzsEXgYdl5iOBHwCn\nA0TEccALgeOAk4H3RoR3KCRJkqR5TLNioFdbZl6fmZeW93cA3wMOL5tj0PKHvnOQmesbHy8G/mt5\nfwpwfmZuAzZGxJXACcDXd8hkO8M/brPeb0Xrc9cTXdpH2e8xlYM8ErVXXQbRlXaYNljIIzAHaZN+\nx91Vv16Ph53uSDPfvr3W9cpvZWtb899Hr+9BV/1qd7fy3dKRpq3fY0/bZTfLbP9brsvaq0f6Xvn6\n5CJJkpalLQz6KNPeIuIo4FFU199PAF4dEb8BXAK8ITNv6bXvqP6i/zLgs+X9YcCmxrZNzPZatLvw\n4lWSJGnBho05qJUhRR8DXlvuILwPOBo4HrgOeFe/8vveOSjjlA7p2PTmzPx0SfMWYGtmntcnq+xX\njiRJkqQZGoQGAAASfUlEQVTeMQffnLqTS6bu7LtvROwJfBz4SGZeCJCZNzS2fwD4dL88+nYOMvOp\n81TgN4FnAL/cWP1j4MjG5yPKuh2sezcz3Ya1j6lekiRJ0mJMTU0xNTU17moMpVfn4NFr9+PRa/eb\n+fz/nXHznO0REcAHgQ2Z+e7G+kMz87ry8bnAZf3Kj8zh/qhfHoP0LuDEzLyxsf444DyqOIPDgS8B\nx2SroIjI/A+GjzmoLSTmYEXrc5dhhsMs9hiGMakxB/3iOXrt22vdoGm6Yg56BfoPEnMwiHHHHJT9\n4tnr+tezJfOtC0ovSdJyEBFk5sBBueMSEfmNfPhAaU+Iy+ccU0Q8EfgX4DvMjtp5M3Aq1ZCiBK4G\nfjszN/fKdzGToP0VsApYX3VU+FpmvjIzN0TEBcAGqsuZV7Y7BpIkSZJ21C+eoJ/M/Crd8cSfW0g+\ni3la0YP6bDsTOHPYvCVJkqTd0SLmORiJxdw5kCRJkjRCW0fwKNPFsHMgSZIkTYjtu/Wdgy2M/nn4\ngwQHdx11HRjaK6i0ud9Cgm+XSnOirl7BxcPWa5Dj69XO7QDxrnxrwwYbd21r61X3Yes1iEHKate5\nXyB9e1sjTfzqugVUTJIk7SqGjTkYFe8cSJIkSRPCmANJkiRJgJ0DSZIkScXuHXMgSZIkacbWOTOj\n7nx2DiRJkqQJ4bAiSZIkSYDDiiRJkiQVPspUkiRJEuCwIkmSJEmFnQNJkiRJgJ0DSZIkScWW3fpR\npjf1WD/dY/0gHanmvnX6O8py77LsOuq7y3J7nzospCPX6xhG5e7G+/bxbG99HuRYBqlvO98udV36\npa23DfLta6bpledCvsVdeSykPoO0wai0zkk8f91OLFySJI2Ddw4kSZIkAXYOJEmSJBXjnudgj7GW\nLkmSJGnGNCsHerVFxJER8eWI+G5EXB4RrynrD4yI9RHxg4j4YkTs3698OweSJEnShJhmxUCvDtuA\n12Xmw4DHAr8XEQ8F3gSsz8xjgYvK554cViRJkiRNiGFjDjLzeuD68v6OiPgecDjwbODEkuwcYIo+\nHQQ7B5IkSdKE2MKqRecREUcBjwIuBtZk5uayaTOwpt++dg4kSZKkCdEVT7AQEbEv8HHgtZl5e0TM\nbMvMjIjst7+dA0mSJGlC9BpWtGnqKn48dVXffSNiT6qOwbmZeWFZvTkiDsnM6yPiUOCGfnnYOZAk\nSZImRK/OwaFrj+XQtcfOfP7GGV+asz2qWwQfBDZk5rsbmz4FvBR4Z1leSB92DiRJkqQJsYh5Dp4A\nvBj4TkR8q6w7HXgHcEFEvBzYCLygXyZ2DiRJkqQJMWzMQWZ+ld7TFDxl0HzsHEiSJEkTYthHmY6K\nnQNJkiRpQmwdwaNMF8POgSRJkjQhFhFzMBJ2DiRJkqQJsdh5DhbLzoEkSZI0IYw5kCRJkgTs7p2D\nm3us395j/SC17bVvc/+ufLa3lv3275embbr1ebHnu51fP4McyyBpR6VdxqR1TbeMuwLd4tR1466C\nJEnaSYw5kCRJkgQYcyBJkiSp8FGmkiRJkoDxDyvqNcXyvCLi7RHx7Yi4NCIuiogjG9tOj4gfRsQV\nEXHSaKoqSZIkLW/TrBzotVSG7hwAf5aZj8zM44ELgbcCRMRxwAuB44CTgfdGxGLKkSRJknYL06wY\n6LVUhr5oz8zbGx/3BW4s708Bzs/MbZm5EbgSOGHoGkqSJEm7iXF3DhZ1TyIi/hR4CXAXsx2Aw4Cv\nN5JtAg5fTDmSJEnS7mCi5zmIiPXAIR2b3pyZn87MtwBviYg3Ae8GTuuRVS6umpIkSdLyt4W9xlp+\n385BZj51wHzOAz5b3v8YOLKx7Yiybgfrzpt9v/bnq5ckSZK0GFNTU0xNTY27GkMZ952DyBzuj/oR\n8aDM/GF5/2rghMx8SQlIPo9qmNHhwJeAY7JVUERkfqpH5s6QPHh+/ThD8rIwyhmSM986srwkSdpV\nRASZGeOux3wiIh+Ylw+U9qp4+JIc02Iuz86KiAdTXa5eBfwuQGZuiIgLgA1Ul4OvbHcMJEmSJO1o\n3PMcDN05yMzn9dl2JnDmsHlLkiRJu6Nh5zCIiA8BzwRuyMyfL+vWAf8N+ElJdnpmfr5fPs4/IEmS\nJE2IRTzK9G+p5hhrSuDszHxUefXtGICjviVJkqSJMWxAcmZ+JSKO6ti0oLgEOweSJEnShNiyddWo\ns3x1RPwGcAnwhsy8pV9ihxVJkiRJE2J6+8qBXgN6H3A0cDxwHfCu+XbwzoEkSZI0Iaa3dw8rmv7K\nV7nnq19dUF6ZeUP9PiI+AHx6vn3sHEiSJEkTolfngMedyB6PO3H28zv+bN68IuLQzLyufHwucNl8\n+9g5kCRJkibE9m3DBSRHxPnAicDBEXEN8FZgbUQcT/XUoquB354vHzsHkiRJ0oS4Z3q4y/PMPLVj\n9YcWmo+dA0mSJGlS9BpWtJPYOZAkSZImxd3jvTy3cyBJkiRNiu3jLd7OgSRJkjQp7BxIkiRJAuwc\nSJIkSSq2jbd4OweSJEnSpJgeb/F2DiRJkqRJsVsPK7p1rKX3V5+YQVpoZ5zEuh5j/sKM1HI6lkEN\n8H2Kl6xb8mpIkqQJdfd4i/fOgSRJkjQpdus7B5IkSZJm2TmQJEmSBNg5kCRJklT4KFNJkiRJgI8y\nlSRJklQ4rEiSJEkS4KNMJUmSJBXeOZAkSZIEjL1zsMd4i5ckSZI0Y/uAr5aI+FBEbI6IyxrrDoyI\n9RHxg4j4YkTsP1/xdg4kSZKkSbFtwNeO/hY4ubXuTcD6zDwWuKh87svOgSRJkjQppgd8tWTmV4Cf\ntlY/GzinvD8HeM58xRtzIEmSJE2K0T6taE1mbi7vNwNr5tvBzoEkSZI0KZYoIDkzMyJyvnR2DiRJ\nkqRJ0R1PAD+agmumFprb5og4JDOvj4hDgRvm28HOgSRJkjQpOuIJADh8bfWqfe2MQXL7FPBS4J1l\neeF8O9g5kCRJkibFkMOKIuJ84ETg4Ii4Bvhj4B3ABRHxcmAj8IL58rFzIEmSJE2KITsHmXlqj01P\nWUg+dg4kSZKkSdEr5mAnWfQ8BxHxhoi4JyIObKw7PSJ+GBFXRMRJiy1DkiRJ2i1sGfC1RBZ15yAi\njgSeCvxnY91xwAuB44DDgS9FxLGZec9iypIkSZKWvSV6lOmgFnvn4Gzgja11pwDnZ+a2zNwIXAmc\nsMhyJEmSpOVv24CvJTJ05yAiTgE2ZeZ3WpsOAzY1Pm+iuoMgSZIkqZ/pAV9LpO+woohYDxzSsekt\nwOlAM54g+mTVORvbuk/Mvl/70OolSZIkLcbU1BRTU1PjrsZwxjysKDLnnUV5x50iHg5cBNxZVh0B\n/Bh4DHAaQGa+o6T9PPDWzLy4lUfmucNXfMnVJ2aQqIydcRLreoz5C6NFGuD7FC9Zt2TFZ751yfKW\nJGlSRQSZ2e8P2RMhIpKnD3ht/rmlOaahApIz83JgTf05Iq4G/ktm3hwRnwLOi4izqYYTPQj4xigq\nK0mSJC1rY36U6ajmOZjp4mTmhoi4ANhA9XfuV+YwtyckSZKk3c0SPqZ0ECPpHGTmA1qfzwTOHEXe\nkiRJ0m5jzEPInSFZkiRJmhTLZFiRJEmSpMVawseUDsLOgSRJkjQpHFYkSZIkCbBzIEmSJKkw5kCS\nJEkSsDweZSpJkiRpBBxWJEmSJAlwWJEkSZKkYhGPMo2IjcBtJZdtmXnCQvOwcyBJkiRNisUNK0pg\nbWbePGwGdg4kSZKkSbH4mINYzM57LLp4SZIkSaOxbcBXtwS+FBGXRMQrhineOweSJEnSpFjcnYMn\nZOZ1EfFzwPqIuCIzv7KQDOwcSJIkSRNvqrx6y8zryvInEfFJ4ATAzoEkSZK0vKwtr9oZc7ZGxGpg\nRWbeHhH3Ak7aIdEAxts5GMUMcCvKchGPfeprqfJdqEmpx3K0Yv4kI1NuFcZp63ZioZIkaTewBvhk\nREB1jf/RzPziQjPxzoEkSZI0MYabBS0zrwaOX2zpdg4kSZKkibH4Z5kuhp0DSZIkaWIMd+dgVOwc\nSJIkSRPjrrGWbudAkiRJmhjeOZAkSZIEGHMgSZIkqfDOgSRJkiTAOweSJEmSCu8cSJIkSQJ8WpEk\nSZKkwmFFkiRJkgCHFUmSJEkqvHMgSZIkCfDOgSRJkqTCOweSJEmSAO8cSJIkSSp8lKkkSZIkwDsH\nkiRJkorxxhzsMeyOEbEuIjZFxLfK6+mNbadHxA8j4oqIOGk0VZUkSZKWu20DvnYUESeX6+8fRsQf\nDlP60J0DIIGzM/NR5fW5UqnjgBcCxwEnA++NiMWUo0WYumLcNVj+bOOlNzU1Ne4q7BZs56VnGy89\n23jp2cZLbfuAr7kiYgXw11TX38cBp0bEQxda+mIv2qNj3SnA+Zm5LTM3AlcCJyyyHA1p6vvjrsHy\nZ+dg6flDtHPYzkvPNl56tvHSs42X2tB3Dk4ArszMjZm5Dfg7quvyBVls5+DVEfHtiPhgROxf1h0G\nbGqk2QQcvshyJEmSpN3AcHcOqK63r2l8HuoavG9AckSsBw7p2PQW4H3A28rntwPvAl7eI6tcaMUk\nSZKk3c/QjzIdyfV2ZC4+n4g4Cvh0Zv58RLwJIDPfUbZ9HnhrZl7c2scOgyRJknaKzOwaDj9RFnp9\n3DymiHgssC4zTy6fTwfuycx3LiTPoR9lGhGHZuZ15eNzgcvK+08B50XE2VS3Mh4EfKO9/65wgiRJ\nkqSdZZHXx5cADyp/tL+W6gFBpy40k8XMc/DOiDie6hbG1cBvA2Tmhoi4ANhANSDqlTmK2xOSJEmS\nOmXm9oh4FfAFYAXwwcz83kLzGcmwIkmSJEm7viWdfyAiNkbEd8okad8o6w6MiPUR8YOI+GLjKUdO\nnjaEHm385xHxvfIkqU9ExH0a6W3jBepq48a2N0TEPRFxYGOdbTyEXu0cEa8u3+fLI+KdjfW28wL1\n+P/ihIj4Rln3zYj4xUZ623iBImL/iPhY+c5uiIjH+Ls3Wh1t/Fh/90avq50b2/ztW84yc8leVMON\nDmyt+zPgjeX9HwLvKO+PAy4F9gSOopofYY+lrN9yePVo46fWbQe8wzYefRuX9UcCn29ut41H287A\nk4H1wJ7l88/ZziNv4yngaeX904Ev28aLauNzgJeV9yuB+/i7t1Pa2N+9ndDO5b2/fcv8tTNmLm4H\nVjy7fOEoy+eU906eNrw5bZyZ6zPznvLxYuCI8t42Hl5XgNDZwBtb62zjxWm38+8CZ2U1mQuZ+ZOy\n3nYeXruNr6O6uALYH/hxeW8bL1D5a/WTMvNDUI3/zcxb8XdvZHq1sb97o9Xnuwz+9i17S905SOBL\nEXFJRLyirFuTmZvL+83AmvLeydOG09XGTS8DPlve28bD2aGNI+IUYFNmfqeV1jYeXtd3+UHA/xMR\nX4+IqYj4hbLedh5OVxu/CXhXRPwI+HPg9LLeNl64o4GfRMTfRsT/jYj3R8S98HdvlLraeHUrjb97\ni9fZzv727R4W87SiQTwhM6+LiJ8D1kfEFc2NmZnR/3muRkvPb4c2zsyvAETEW4CtmXlen/1t4/l1\nfY9PB5pjKvs9esw2HkxXO68EDsjMx5ax8BcAD+ixv+08v642fivwmsz8ZEQ8H/gQ1RCNLrZxfyuB\nRwOvysxvRsS7qTpfM/zdW7RebfzH4O/eCHW18xnAk/C3b9lb0jsHWeZBKEMBPkl1i2lzRBwC1VwJ\nwA0l+Y+pxrHVjmD29rZ66NHGRMRvAs8Afr2R3DYeQkcbn0j1V5VvR8TVVO347xGxBtt4aD2+y5uA\nT5T13wTuiYiDsZ2H0qONT8jMT5YkH2N2KIBtvHCbqP6q+s3y+WNUF1jX+7s3Mr3a2N+90epq50dR\nxRP427fMLVnnoNx+2q+8vxdVT/MyqknSXlqSvRS4sLz/FPCiiFgVEUfTY/I0zerVxhFxMvAHwCmZ\neXdjF9t4gXq08Tcyc01mHp2ZR1P9J/roMmzANh5Cn/8vLgR+qaw/FliVmTdiOy9Yjza+HLgyIk4s\nyX4J+EF5bxsvUGZeD1xTvqsATwG+C3waf/dGolcb+7s3Wj3a+d8z8xB/+5a/pRxWtAb4ZETU5Xw0\nM78YEZcAF0TEy4GNwAvAydOG1KuNfwisoho2APC1zHylbTyUzjZupZlpQ9t4aL2+y3sCH4qIy4Ct\nwG+A7Tykrjb+QkTcBLwnIvYC7gJ+C2zjRXg18NGIWAVcBZxGNRmRv3uj027jlwHfxN+9Uev6Ljf5\n27dMOQmaJEmSJGDpn1YkSZIkaRdh50CSJEkSYOdAkiRJUmHnQJIkSRJg50CSJElSYedAkiRJEmDn\nQJIkSVJh50CSJEkSAP8/lTiNCGS3Os0AAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig, ax = plt.subplots(1,1,figsize=(15,4))\n", "t=0\n", "mesh=ax.pcolormesh(xx,zz,Ss[t,:,:])\n", "ax.set_title('Sal: t={} : {}'.format(t, run))\n", "fig.colorbar(mesh,ax=ax)\n", "ax.set_xlim([500,650])\n", "ax.set_ylim([-50,0])" ] }, { "cell_type": "code", "execution_count": 14, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "(-50, 0)" ] }, "execution_count": 14, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAw8AAAEKCAYAAABUjmlMAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xu4JFV56P/vywwDcpHhYsBhxoABkjNBI0ZHoz9ljhdE\no6DPiaIejReOufBTczExIDkRYoLXEI05/JITxR8mASUKiklQB4870UQuEuTigIJhhAEZBAVHuc3l\nPX9U1ezaNd29a3fvnq6Z/f08Tz/dXbVq1Vure2bX6lpvrchMJEmSJGk2u006AEmSJEk7BzsPkiRJ\nklqx8yBJkiSpFTsPkiRJklqx8yBJkiSpFTsPkiRJklqx8yBpQYqIjRFx2KTj2JEiYl1EPHee63x9\nRHxlyG0Pi4itEbFb+f6fI+K1tfV/EhHfj4g7y/cvi4jby8/uF+bnCOZfRNwQEc+edBySNA52HqSd\nVET8uDyJ2liegD1Qe/+qSccH205WnzPktk+KiKsj4icR8fVRThYjYioiTq4vy8x9M3PdsHUOEcOM\nE+UJyfKxQ0XEX0fEm2Yrl5kvysy/Lbd5HPC7wM9l5rKyyAeAU8rP7tqIeHP53XgoIj42D3E+PiL+\nMSJ+VHZa3jtMPZl5dGb+66jxSFIX2XmQdlKZuU95ErUv8F3gxdX7zLxg0vGVEoi5bhQRS4DPAh8H\nlgLnAZ+NiN1HiKMr5tweu4DjgX+a4zaPA+7NzHsBIiLKZWtrZe4A3gWcO2qA5XduDXAZcDBwKPB3\nc6xj8ahxTKJuSZoLOw/SLiYidouIUyPiloi4JyI+GRH7l+uqX79fHxG3RcS9EfEbEfHUiLguIn4Y\nER+u1fX6iPi3iPhwRNwXETe2vZIQEX9LcbL3ufJqyO/N4TBWA4sy80OZuSkzP0xx0j3nqxgR8afA\ns4C/LOP4i3L51oh4fPn6/4+Ic8phMxsj4isRcUhEfKhskxsj4km1OpdFxKcj4u6I+M+IeEtt3ary\n1/D7I+KuiPhAuar6Jfq+ch9PK8u/MSLWRsQPIuLz5S/u1GJ8S0R8p/wl/H3lSfRsx/ymss4fRcQ3\n67EDx0TEteXn+YmI2KO23Ysj4hvlMf9bRDyhtm5FRFxUHvM99e9JY9/vL9tv3/L9E4H7MvPOiFgU\nER8oj+U7wC83tp2KiJOjGFr1RWBZ2VbnAz8CFgHXRsTNAJl5cWZ+Fri3Tyx9j6eH1wPrM/ODmflg\nZj6SmdcPKF/tY11EvD0irgM2lse4LiKeU35PHqj+/ZXljymPf1H5frbP/5TyeL81WyyStCPYeZB2\nPW8BTgCeDTwW+CHwvxplVgFHAK8EPgS8g+LE/OeBV8TM8dqrgFuAA4F3AhdFxFKAspPyuV5BZOZr\ngduYviLygXKb+8qTuV6Pt5eb/zxwXaPKa8vl24mIV0fEtX3iOB34CvD/lnG8tVc54OXA6cBBwCPA\n5cBVwAHAp4Czy33tBnwOuAZYBjwX+O2IOK6s50PAn2fmfsDjgX8olz+rfN6vjOOKiDgROA14Wbnf\nrwDNq0YvBX4ReDJwIvDGPvFXbfFyis/ptZn5aIrvwg+q1eVxvgA4HHgixUkzEXEM8FHgTeUx/zVw\nSUTsXp7o/iNwK/DTFL/Kz4gzCn8DHA08PzM3lqteVG5LWfcvA08CngL8CjOvCiWQmfkl4IXAnWVb\nvbq8wgbwxMw8snnYPdqh3/Es6dN0Twe+W3Ygvx8RX46Io/uUbXplGe/SzNxSHVNm3gl8DfhvtbKv\nBv4hM7e0/PxPBJ4KrGwZiySNlZ0Hadfz68AfZuadmbkJOBP4lZg51v5d5S+ra4CNwPmZeU95svMV\n4Jha2bvLKwBbMvNCil9AXwyQme/JzJfMJbjMXJqZ+/d5vK8stg9wf2PTHwH70kNmnp+Zs+VEDPrF\nPoGLMvOazHwYuBj4SWb+XWYmcCHTbfJU4KDM/JPM3JyZtwIfoTiBhKLjcWREHJSZD2TmFQP2/xvA\nuzPzW5m5FXg38KSIWFEr897MvC8zbwc+CMyWz/I/ym2uBsjM72TmbbXj/IvMvCszf0jRCaquSvwa\n8NeZeVUWPg48DPwSRQfyscDvl7/KP5yZ/17b5+7AJyiGmL0kMx+qrXsR8M/l61dQdKzuKPd/Vp92\nYcDyXnoNS+t3PE/vU8dypjvTj6UYZtVmqFzVpneU352m8yk/s/Kq0UnlMmj3+b+7/Px71S1JO5yd\nB2nXcxhwcfVrPsUY8c0U47grG2qvH+zxfu/a+zsa9X+X4hf3cdoIPLqxbD+KDsSwZst7uLv2+qHG\n+wcpOjRQ/PK+rH7FhOLX458q158MHAXcGBFXRsSMoTkNPw18qFZPNfzm0FqZ22uvb2P2tl8OfGfA\n+rtqr5vH9bbGcS2nOJFeAXy3PMHt5QjgJcAfZ+bmamF5herngKqj8dgexzMfenU0+h5PeaWqurlA\nlYvxIPCVzPxC2Sn8AMXVtp9rsf/bB6y7CPiliDiE4mrg1sz8ai3GuXz+kjRxdh6kXc9twPGNX/T3\nyszvDVnfoY33P832HYp+tjthj5l3iWo+Ti2LraUYUlP3ROCbc4p8QBwjuB24tdG+j87M6mrMLeUw\nm8cA7wU+FRGP6hPDbcCvNeraOzMvr5V5XOP1bG1/O8XJfFtVXLcBf9qIZZ/M/GRZ5+Oqcfo93Egx\nnOrSiDiqtvwFwJfKqzcA3+txPPOhX9v2PJ7ySlV1c4Gqczdj2Ft5lWCU/RcriissX6S44vBqZg5L\navP5dynZX5LsPEi7oL8CzqoSLyPiMRFxwhzrqJ84/VREvLUc+/5y4GeZHoYymw3Az9QXZO0uUT0e\n7ymLfRnYUu53j4h4K7AV+D9zPI6+cTTM5UTxSorE2LdHxKPKBNmjI+IpABHxmoh4TFn2foqTv63A\n98vnehx/BbwjIlaW2+5XtnHd70XE0nIoy1uBT84S30fKbZ5c5iEcUU/C7aE69r8BfiOKhO+IiL0j\n4pcjYh/gCooT//dExF4RsWdEPKNeSWZ+giJ35rIoE9EphizV77J0IfDWiDg0iiTiUxlB2fZ7AouB\nReV3pergDDqeXv4OeHpEPLes47cpPrMbR4mxdD7wOorch/Nry9t8/pLUKXYepF3Ph4BLgC9GxI8o\nEjZX1da3+SWzXuYK4EiKE6l3Ab9S/ppKRLwjIgZ1JN4N/GE5LON32x5AmavxUuBXKRK+fxV4aX1I\nTF1E/PeIuGFAlR+iyPv4QUR8sNcu6ZG426MMZULsiylyBf6Tol3+N9PDrF4A3BARG4E/B15Z5gg8\nAPwp8G9le6zKzM9QXJ34RETcD1xfbl/3WeBqigTtf2SW25Jm5qfK/VR3KLoI2L9f8dpxXU2RXPyX\nFAnWN1O0O+VwpZdQXNG4jeJKxCt61PFx4I+BL0UxAd9xwOdr+/sb4AsUv/J/Hfg0g7+PPT+Dmv8J\nPAD8AfAaiqFHp892PD13lPntso6/Ksu/BDih33duji6haLvvZe0OTi0+f686SOqcmL6aPM8VRxxP\nkdy3CPhIZg412Y6kyYmI1wMnZ+azZiur+RcRW4EjMvM/Jx3LXEXEKopE4n4JypKkndBYrjyUl3z/\nkmJioJXAqyLiv4xjX5KkTkqKW8ZKknYh4xq2tAq4JTPXlcMPPkFxr2pJO5dew3e04/Rs+4j4qz4J\n5+fs6AD7KW+R+oVJxzGqiHhcn7b+UUQsn3R8krSjjWXYUkT8CvCCzHxT+f41wNMy8y2Dt5QkSZLU\nVeO68uAvlZIkSdIuZvGY6r2DYlKhygpgfb1ARNjBkCRJ0g6RmXO5LfdEzPX8eBLHNK7Ow9eBI8tb\n9d1JMTnOq5qFPp0vbFXZHjwCwMMs2bZsMVsA2JeNACwq39dtZF8AlvAwAA+y14yyB26bzBMWMfNu\nfPdy0Hb1VWWWch8AW3o0330snRHXwTMm7p1pQznh7xaK25L/VI+yj7DHjGOoVMcG022xpGynPWpl\n/+KM+3nDGdPDcqu4fljGWanauF6mvo+m6rOobwewFw9se121c7Ns9b5X+zXrqX/mVfmqvXqp9lmV\nWcadM9Zvrm1b1ffItngWbVdm8bZj2GNGmfrnce4Zd/L2M/bggfL71YzvwC3T37O9flK0wYZHHzhj\n34+qtVvzOKtj2nfLxm1lzll0CgCn3l7edfT+4un0o/9nEV+8a1vZp+Rzt73+ZS5jZ3TGGWdwxhln\nTDqMXZ7tPH628fjZxuO3M7bx3OZ9nKw/aVnuD8caRX9j6Txk5uaIeDPF/bwXAR/NzPmYaEc7iV6d\nOUmSJA22+6QDmMW4rjyQmZcCl46rfkmSJGlXM7aT83nS9fg0oqet3mPSIezynrm6/zAqzY/Vq1dP\nOoQFwXYeP9t4/Gzj8bONx+tRkw5gFnYednFPW70nG2cvphE8c7X/jMbNP1Q7hu08frbx+NnG42cb\nj9eCHbY0nx7Ylug8ndRcJbA2k3rrY+2bib+LGwm83+Ko7fZVJfVWSdG9fI9l2+2ruf8qcfrOsuz0\n+ulj2L/cR7XNOg4H4IFan7M6hsWNfe3To0tQ7bMZVz2B+p4yEbyZIF5PTGZAovRsHmgkpdf3NXMf\ngzWT3duqkoyrfd3KYQDsxYPAdHs2Y5yrKpEdpj/jKuG9SsyvrFt02LbX+z662H+VRH3vogNb7Kt9\nu0mSpJ1b10/Oux6fJEmStGB45UGSJElSK10/Oe96fJIkSdKC0fUrD5E5mYmeIyI/mq9uNZ67Gpte\nL1stq0+u1VQflw7bT961pDHBWVHv5rJs/wnJ9mrssz6pWDXevaq7GmPfnPCreTz1sr0mMmvGVc+L\naE6+Vo3vr/Z5L9Pj6qtcibvLCeoeaIzPBziIe4DpvINeZap1zbjq+RXNnIKqTR4sY9/co22bx1RX\nxVEde6/PqMqRqD77FdwOTH9m9YkBq+2qXJFm/TCdazK9rniu505UEwFW38VmHk71edS3q/bxDY4B\nZubYVN+LKpeiiv0Z/Pu2Mm/+o48A8Nw//icA3s2pAKx64vUA5E9NT4az7yXf3/Z6417bT34oSdKu\nLiJ2mhmmL2lZ9gR2rRmmtcA5SZwkSdLceatWSZIkSa10/eR8t0kHIEmSJKmwe8tHU0SsiIgvR8Q3\nI+KGiHhrY/3bImJrRBxQW3ZaRNwcETdFxHFt4ut650aSJElaMEY4Od8E/E5mfiMi9gGujog1mXlj\nRKwAng98tyocESuBk4CVwKHAZRFxVGZuHVN8o/sML92WXPpwI7kZtp/IrJ5EWyWeNidzqyerNhNX\nK0v5IQAHc/d266qE2I3ss90+qwTbKtG3WldP2G2O9a+SaKt9Vcm5MJ3gXCXlVtvWj6mq+54y6bkq\nW598rkoSPqhMBt5QJkNX9dT3uVefpN76MVTbVfFt6fE1aU7GViUFV9vWE6erNm0eS30ytSqZutqu\narf72H9bmfpx1NXbvJnMXh1ntbye6Fy1U5VQXtVTb5sqsblfUjrAgWWCeVWmqrc6lvo+D2PdjH1e\nUyZMH8EttX0W7VJ9xj/LtwC4jif0PP5ennHZ/5l+0/+eApIkqWOGvdtSZt4F3FW+/nFE3AgsA24E\nzgbeDny2tsmJwAWZuQlYFxG3AKuAywftxysPkiRJUkfMx8l5RBwGHANcEREnAusz87qIGTdnWsbM\njsJ6iisQY49PkiRJ0jzod+XhyvIxm3LI0qeA3wK2Au+gGLK0rciAzWedw8HOgyRJktQR/W7Vemz5\nqJzTo0xE7A58Gvi7zPxMRDwBOAy4trzqsJwiF+JpwB3Aitrmy8tlA010kjj+PuHH5YJ6esND5XOz\na7NP7fVBjedqm7tqZWamQ8Ce5fPy8vmQ2roqjrsa7+v7bNa7tFEW+nfHqjjrZX/cWFfVt2etzEON\n52rf9f0091kNsf9J+bx3bV01R9qejeeH2F61XZVSUP+Mqn1W21Wx19u0SqO4v1G2auN6fVU7L208\n31Mrc0+jbFV/PX1j6cwyS5b/qHjes8hH2Guv6QSAu28r80ZuWjwz9ul0hun2OpCZaj8LPOqIH85Y\n9eBdRZ7GbvsVH8C++0/nPDxhSTGJWzVJ37f/9YlFuE+ZPtAHNhb/bWz9TvEBPPYZtwKwjDu3lbn6\nyc8sXry5XPCk8vkXPwnAKTn9b/+cs3532+t8B5IkLTg70yRx61uWXc7MSeKi6B2cB9ybmb/Tp/5b\ngV/MzB+UCdPnU+Q5HApcBhyRs3QOvPKg8eg/ebQkSZL6GOHk/JnAa4DrIuKactk7MvPSWpltHYPM\nXBsRFwJrKX7SPWW2jsOI8UmSJEmaT7u3PTufeVNSMvOrzDKHW2Y+vvH+LOCsOYRn50GSJEnqisVD\ndh52FDsPkiRJUkfs3vGh35PtPHyVYi68pi2N91Uj1u9dVSXzVinpVZLrT2plmnVX21fzce1XW1cl\n9VYJydW2vVLeq331upfW4sZzpZn4DNsf196N94O2H5QwXfVEe8VZHVfVxvs0ytbX1RO36/uG6dir\nueCqOf6q9qrHVG1X7bvaplcCdrOeelzVZ1vFVW1fn1+wSp4u43tk30cXz2X9P64fU5WjXCVwN5Oi\nYfvk+Kptav+wHzy8nMiukai+de/iA71//+mM9a/+P4fM3OdlZVzrDtpWZttxlWW+d9/hxfPyw6fL\nVPF8tXx+EjOcE9NJ2nwUSZK0k2h95WFCOh6eJEmStHDsvsfsZSbJzoMkSZLUFR0/O+94eJIkSdIC\n0vGz88mGd1ef5f2yx+vRNvMFmmP5e9VTbd+cEA62zymotu21z2ZOBj3KNA2aaG22bWfb51w02+S+\nPsuhfy5Fr3VVbkJzYr62qs+v2sdPBpRpqsdV5VNU8VWTvPVq26pNHxpQtirTK57KrY2yVTzV96w+\n0V1zsr/vl8+9jq3aZ1XfTbV1hzTKXlC92Kt4+pd3Tq/7RK3cG3vsR5IkdYedB0mSJEmteLclSZIk\nSa10/Oy84+FJkiRJC4h3W5IkSZLUSsfPzicb3kOzF5mhV8LuwwPK9Nu+VzJ0c12b+pqxDCrfa/lc\n6h7XFORt2msu69p8owYlZzcTnnttN0zC+iBtktEHtUWvZP3ZYqmOs9qmnmjevAnAPT3KVKq54JoT\n3B178/Tr3ztyQCCSJKlT7DxIkiRJasWEaUmSJEmtdPzsvOPhSZIkSQtIx8/OOx7eAKPkAPTLbxi1\nvnEYZ93jMGy8w+SatKlnLjkYVQ5Or8uFVf7B7j3W9dt3rxiqvIXmBHX1Ms08nl45GVU9ezfeH/OS\n4rk+idyGPnFJkqTu6fjZecfDkyRJkhYQb9UqSZIkqZWOn53vNuyGEfHyiPhmRGyJiCc31p0WETdH\nxE0RcdzoYUqSJEkLwKKWjwkZuvMAXA+8DPjX+sKIWAmcBKwEjgfOiYhR9iNJkiQtDItbPhoiYkVE\nfLn8cf+GiHhrufyAiFgTEd+OiC9GxNLaNnP+wX/oCyOZeVO50+aqE4ELMnMTsC4ibgFWAZcPuy9p\nhkHJ1M1v9FwSr4dN0t40h/XNfQzadlCS9t49lvVbv7Ml3UuStJANP2xpE/A7mfmNiNgHuDoi1gBv\nANZk5vsi4g+AU4FTGz/4HwpcFhFHZebWQTsZxxWBZcD62vv1ZUCSJEmSBhly2FJm3pWZ3yhf/xi4\nkeIc/ATgvLLYecBLy9fbfvDPzHVA9YP/QAP7NmVv5ZAeq96RmZ+brfKa7Ln05jOmXx+wGg5cPYcq\nJUmSpO1NTU0xNTU16TCGMw8J0xFxGHAMcAVwcGZWN27fABxcvl7GzJFBrX7wHxheZj5/jrEC3AGs\nqL1fXi7b3pFnDFG9JEmS1N/q1atZvXr1tvdnnnnm5IKZqz1H27wcsvRp4Lcyc2M9xSAzMyJ6/6hf\nFpmt/vm6GVQ98eES4PyIOJui93IkcOU87UcabJhJ4maro15PrwnbRjGovkHrflI+V7kN9w7YdtDE\ndpIkqVv63Elpan3xGCQidqfoOPxtZn6mXLwhIg7JzLsi4rHA3eXy9j/41wzdeYiIlwF/ARwE/FNE\nXJOZL8zMtRFxIbCW4jTslMyctRcjSZIkLXh9zs5XH1Y8KmdeMXN9FJcYPgqszcwP1lZdArwOeG/5\n/Jna8jn/4D/K3ZYuBi7us+4s4Kxh65YkSZIWpOHHBT0TeA1wXURcUy47DXgPcGFEnAysA14BMOwP\n/h2fw06SJElaQIacAC4zv0r/O6k+r882c/7B387DzqrX/AC7D1jXz64+Hn5Q/sKo9TTrG2Y+hc21\n535xtcl9qJKrDiyfNwL7lq+/OkRckiRpMjp+dt7x8KR51tVv/HzHte/sRSRJUgd19Vyl1PHwJEmS\npAVkj0kHMJidB0mSJKkrOn523vHwJEmSpAWk42fnkw1vLom9Tc1E3/lKEh4lpkkbJvZqm109cbrS\nJql5rv8q5pIo3a/sMMnWvVSf49fnqT5JkrRjDXm3pR2l430bSZIkaQHp+Nl5x8OTJEmSFpCOn513\nPDxJkiRpAXHYkiRJkqRW9py9yCTtvJ2HURKbd+akaI3fXJOX5zLDdPNf3OY+y3utq9ffTHTfWD5/\nq0UMkiSpuzp+dt7x8CRJkqQFxGFLkiRJklrp+Nl5x8OTJEmSFpCOn513PDxpJzDuSeIGrTN/R5Kk\nXYvDliRJkiS14t2WJEmSJLXilQdJkiRJrXT87Lzj4UmSJEkLSMfPzjsennaYXom3u/dYNqh8fZth\nE3kH7XNcmhOujav+XvsYtO82bbilfL5lrkFJkqRO6vjZ+W6TDkCSJElSaVHLR0NEnBsRGyLi+sby\nt0TEjRFxQ0S8t7b8tIi4OSJuiojj2obX8b6NJEmStIAMf3b+MeDDwMerBRHxX4ETgCdm5qaIeEy5\nfCVwErASOBS4LCKOysyts+3EKw+SJElSV+zR8tGQmV8BfthY/JvAuzNzU1nm++XyE4ELMnNTZq6j\nGAC9qk14dh4kSZKkrljc8tHOkcCzI+LyiJiKiKeUy5cB62vl1lNcgWgVntTbMEnPo854PMkZk3fE\nvvvtY6773jJ7EUmStBPqc3Y+9e/FY4ja9s/Mp0fEU4ELgcf3KZsjhCdJkiRph+tzdr762cWjcubZ\nrWpbD1wEkJlXRcTWiDgIuANYUSu3vFw2K4ctSZIkSR2Ri9o9WvoM8ByAiDgKWJKZ9wCXAK+MiCUR\ncTjF8KYr21TolQdJkiSpI7YMeXYeERcAxwIHRsTtwB8B5wLnlrdvfQT4VYDMXBsRFwJrgc3AKZnZ\nathStCw37yIied5k9i3t9Jo5Dy0nicvb5j0SSZI6LyLIzJh0HLOJiHzoJ+3K7rk3EzkmrzxIO5Mh\nOw2SJGnn8PAeS1qWfGSscfRj50HaWXiHJUmSdnlbFrVPaJgEOw+SJElSR2zBzoMkSZKkFjbbeZAk\nSZLUxpaOn54PPc9DRLw/Im6MiGsj4qKI2K+27rSIuDkiboqI4+YnVEmSJGnXtoVFrR6TMsokcV8E\nfj4zfwH4NnAaQESsBE4CVgLHA+dEhJPRSZIkSbPoeudh6Osimbmm9vYK4L+Vr08ELsjMTcC6iLgF\nWAVcPnSU0s5k0F2RFrUoI0mSFqyHaXur1smYr0FVbwQuKF8vY2ZHYT1w6DztR5IkSdpldT3nYWB0\nEbEGOKTHqndk5ufKMqcDj2Tm+QOqcippSZIkaRY79a1aM/P5g9ZHxOuBFwHPrS2+A1hRe7+8XLa9\n75wx/Xr/1XDA6kG7kyRJkmY1NTXF1NTUpMMYStc7D5E53EWBiDge+DPg2My8p7Z8JXA+RZ7DocBl\nwBHZ2FFEJM/zgoR2QTsy5+GWuRXP2+Zpv5Ik7UQigsyMSccxm4jIK/PoVmVXxQ0TOaZRBlV9GFgC\nrIkIgK9l5imZuTYiLgTWApuBU5odB0mSJEnb26lzHgbJzCMHrDsLOGvYuiVJkqSFqOvDlrrdtZEk\nSZIWkEcWyK1aJUmSJI1os1cetMvo6sRm9X9jWxrLuhbzKPFsrr1eN2IckiSpk3bZnAdJkiRJ88uc\nB0mSJEmtdL3zsNukA5AkSZJU2MyiVo+miDg3IjZExPW1Ze+PiBsj4tqIuCgi9qutOy0ibo6ImyLi\nuLbx2XmQJEmSOuIR9mj16OFjwPGNZV8Efj4zfwH4NnAabJvU+SRgZbnNORHRql9g50GSJEnqiC0s\navVoysyvAD9sLFuTmVvLt1cAy8vXJwIXZOamzFwH3AKsahOfOQ+SJElSR4zxVq1vBC4oXy8DLq+t\nWw8c2qYSOw+SJElSR4zjVq0RcTrwSGaeP6BYtqnLzoMkSZLUEf3utvStqbv41tSGOdcXEa8HXgQ8\nt7b4DmBF7f3yctms7DxIkiRJHdGv83DE6kM5YvX0yKJ/PPO6WeuKiOOB3weOzcyHaqsuAc6PiLMp\nhisdCVzZJj47D5IkSVJHDDvPQ0RcABwLHBQRtwPvpLi70hJgTUQAfC0zT8nMtRFxIbAW2AyckpkO\nW5IkSZJ2Jg/3vg3rrDLzVT0Wnzug/FnAWXPdz2Q7D1t20H7m0oHbUTEtRJsb7+fr29frM5vL51jF\n1SueZsyTtG7SAUiSpHHr+gzTXnmQJEmSOsLOgyRJkqRWxjjPw7yw8yBJkiR1xDjmeZhP3Y5OkiRJ\nWkActiRJkiSpFTsPkiRJklp5mCWTDmEgOw+SJElSR5jzIEmSJKkVhy1JkiRJasXOgyRJkqRWnOdB\nkiRJUivmPEiSJElqxWFLkiRJklp5xFu1SpIkSWrDnAdJkiRJrZjzIEmSJKkVcx4kSZIktWLnYZDN\nI2zbjHxQXdW6QUc7SiwaTtfavGvxVNZNOgBJkrSjmPMgSZIkqRVzHiRJkiS1MsqtWiPiNOA1wFbg\neuANwN7AJ4GfphjP8IrMvG/Yfew2dHSSJEmS5tVmFrV6NEXEYcCbgCdn5hOARcArgVOBNZl5FPCl\n8v3Qhu48RMS7IuLaiPhGRHwpIlbU1p0WETdHxE0RcdwoAUqSJEkLxRYWt3r08CNgE7BXRCwG9gLu\nBE4AzivLnAe8dJT4Rrny8L7M/IXMfBLwGeCdABGxEjgJWAkcD5wTEV7hkCRJkmaxhUWtHk2Z+QPg\nz4DbKDpTXKzZAAAORklEQVQN92XmGuDgzNxQFtsAHDxKfEPnPGTmxtrbfYB7ytcnAhdk5iZgXUTc\nAqwCLh86SkmSJGkB6Her1o1T/8HGqWv6bhcRPwP8NnAYcD/wDxHxmnqZzMyIyFHiGylhOiL+FHgt\n8CBFBwFgGTM7CuuBQ0fZjyRJkrQQ9Os87LX6qey1+qnb3n/vzI81izwF+PfMvBcgIi4Cfgm4KyIO\nycy7IuKxwN2jxDdwOFFErImI63s8XgKQmadn5uOAjwEfHFDVSD0cSZIkaSF4mD1aPXq4CXh6RDwq\nIgJ4HrAW+BzwurLM6yjSDYY28MpDZj6/ZT3nA/9cvr4DWFFbt7xctr3vnjH9er/VsHR1y91JkiRJ\nvU1NTTE1NTXpMIYy7AzTmXltRHwc+DrFrVr/A/jfwL7AhRFxMuWtWkeJLzKHuygQEUdm5s3l67cA\nqzLztWXC9PkUw5gOBS4DjsjGjiIiedYIFyTmMsN0v23mur00Cevmr6q8bf7qkiRpZxERZGZMOo7Z\nRET+TN7Qqux34uiJHNMoOQ/vjoifBbYA3wF+EyAz10bEhRSXSTYDpzQ7DpIkSZK212sOhy4Z+srD\nyDv2yoPUzrr5q8orD5KkhWhnuvKwvBjYM6v1ceROd+VBkiRJ0jwaNudhR7HzIEmSJHWEnQdJkiRJ\nrTz8yJJJhzCQnQdJkiSpI7Zs7vbpebejkyRJkhaQLZsdtiRJkiSpBTsPkiRJklrZvMnOgyRJkqQW\ntm7p9ul5t6OTJEmSFhKHLUmSJElq5aFun553OzpJkiRpIdk86QAGs/MgSZIkdYWdB0mSJEmt2HmQ\nJEmS1MqmSQcwmJ0HSZIkqSu2TDqAwew8SJIkSV3hsKUxGaZhO/5hSACsm3QAkiRpYh6adACD7Tbp\nACRJkiSVNrd89BERiyLimoj4XPn+gIhYExHfjogvRsTSUcKz8yBJkiR1xYidB+C3gLVAlu9PBdZk\n5lHAl8r3Q7PzIEmSJHXFCJ2HiFgOvAj4CBDl4hOA88rX5wEvHSW8nTfnQZIkSdrVjHar1j8Hfh94\ndG3ZwZm5oXy9ATh4lB145UGSJEnqii0tHw0R8WLg7sy8humrDjNkZjI9nGkoXnmQJEmSuqJfPsON\nU3DT1KAtnwGcEBEvAvYEHh0RfwtsiIhDMvOuiHgscPco4UXRAdnxIiJ51mT2LXXauvFVnbeNr25J\nkroqIsjMnr/Gd0lEJH/d8vz41/sfU0QcC/xeZr4kIt4H3JuZ742IU4GlmTl00rRXHiRJkqSumL95\nyapeyHuACyPiZIqfKF8xSqV2HiRJkqSumIfOQ2b+C/Av5esfAM8bvdaCnQdJkiSpK+bvysNY2HmQ\nJEmSumK0W7WOnZ0HSZIkqSt63Ia1S+w8SJIkSV3x0KQDGMzOgyRJktQV5jxIkiRJasWcB0mSJEmt\nmPMgSZIkqRWHLUmSJElqxc6DJEmSpFY6nvOw26gVRMTbImJrRBxQW3ZaRNwcETdFxHGj7kOSJEla\nEB5u+ZiQka48RMQK4PnAd2vLVgInASuBQ4HLIuKozNw6yr4kSZKkXV7Hhy2NeuXhbODtjWUnAhdk\n5qbMXAfcAqwacT+SJEnSrm9Ty8eEDN15iIgTgfWZeV1j1TJgfe39eoorEJIkSZIG2dLyMSEDhy1F\nxBrgkB6rTgdOA+r5DDGgquy59LtnTL/ebzUsXT0oHEmSJGlWU1NTTE1NTTqM4XR82FJk9j6vH7hR\nxNHAl4AHykXLgTuApwFvAMjM95RlPw+8MzOvaNSRPGvu+5Z2eevGV3XeNr66JUnqqoggMwf90N0J\nEZG8sOX58aWTOaahEqYz8wbg4Op9RNwK/GJm/iAiLgHOj4izKYYrHQlcOR/BSpIkSbu0jt+qdb7m\nedjWRcrMtRFxIbCW4sLLKTnM5Q1JkiRpoZngbVjbGGrY0rzs2GFLUm/rxle1w5YkSQvRTjVs6Zda\nnh9/bScatiRJkiRpDDo+bGnkGaYlSZIkzZMRbtUaEcdHxE0RcXNE/ME4wvPKgyRJktQVQ96qNSIW\nAX8JPI/iLqhXRcQlmXnj/AXnlQdJkiSpOza3fGxvFXBLZq7LzE3AJ4AT5zs8rzxIkiRJXTF8zsOh\nwO219+sp5mCbV3YeJEmSpK7od6vWh6bg4alBW+6Q25jaeZAkSZK6ol/Ow+LVxaPyozObJe4AVtTe\nr6C4+jCv7DxIkiRJXTH8sKWvA0dGxGHAncBJwKvmJaYaOw+SJElSV/S5DetsMnNzRLwZ+AKwCPjo\nfN9pCew8SJIkSd0x5K1aATLzUuDSeYulBzsPkiRJUleM0HnYEew8SJIkSV0xfM7DDmHnQZIkSeqK\njl95cIZpSZIkSa3YeZAkSZLUisOWpElYN+kAJEmS5s7OgyRJktQZ3c6YtvMgSZIkdUa3M6btPEiS\nJEmd4ZUHSZIkSa08OOkABrLzIEmSJHWGVx4kSZIktWLOgyRJkqRWvPIgSZIkqRWvPEiSJElqxSsP\nkiRJklrxbkuSJEmSWnHYkiRJkqRWHLYkSZIkqZVuX3nYbdIBSJIkSapsavmYu4h4W0RsjYgDastO\ni4ibI+KmiDhutjq88iBJkiR1xniuPETECuD5wHdry1YCJwErgUOByyLiqMzc2q8erzxIkiRJnTG2\nKw9nA29vLDsRuCAzN2XmOuAWYNWgSrzyIEmSJHXG/N+qNSJOBNZn5nURUV+1DLi89n49xRWIvuw8\nSJIkSZ0xdD7DGuCQHqtOB04D6vkM0aNcJQftx86DJEmS1Bn9ch6+VT56y8zn91oeEUcDhwPXllcd\nlgNXR8TTgDuAFbXiy8tlfQ2d8xARZ0TE+oi4pny8sLZuTlnbkiRJkqB/jsPjgRfWHu1k5g2ZeXBm\nHp6Zh1MMTXpyZm4ALgFeGRFLIuJw4EjgykH1jZIwncDZmXlM+bgUtsvaPh44JyJMzJ6U+6YmHcGu\nzzYeu6mpqUmHsCDYzuNnG4+fbTx+tvG4bW75GNq2YUmZuRa4EFgLXAqckpkDhy2NelLfa7zUnLO2\nNUb3T006gl2fbTx2/qHaMWzn8bONx882Hj/beNzGN88DQGY+PjN/UHt/VmYekZk/l5lfmG37UTsP\nb4mIayPioxGxtFy2jOJySGXWrG1JkiRJsAOuPIxkYML0LFnb/x/wx+X7dwF/Bpzcp6qBlz8kSZIk\nwThu1TqfYpZhTe0qiTgM+FxmPiEiTgXIzPeU6z4PvDMzr2hsY4dCkiRJO0RmDro9aSfM9fx4Esc0\n9K1aI+Kxmfm98u3LgOvL15cA50fE2RTDlXpmbe8MH6AkSZK0o+wM58ejzPPw3oh4EsWQpFuBX4ci\nazsiqqztzbTI2pYkSZLUffMybEmSJEnSrm+s8y9ExLqIuK6cRO7KctkBEbEmIr4dEV+s3aXJyeWG\n0KeN3x8RN5Z3wrooIvarlbeN56hXG9fWvS0itkbEAbVltvEQ+rVzRLyl/D7fEBHvrS23neeoz/8X\nqyLiynLZVRHx1Fp523iOImJpRHyq/M6ujYin+XdvfvVo46f7d2/+9Wrn2jr/9i1kmTm2B8VwpgMa\ny94HvL18/QfAe8rXK4FvALsDh1HMD7HbOOPbFR592vj5VdsB77GN57+Ny+UrgM/X19vG89vOwH8F\n1gC7l+8fYzvPextPAS8oX78Q+LJtPFIbnwe8sXy9GNjPv3s7pI39u7cD2rl87d++Bf7YETM/NxM/\nTii/kJTPLy1fO7nc8Ga0cWauycyt5dsrgOXla9t4eL0SmM4G3t5YZhuPptnOvwm8OzM3AWTm98vl\ntvPwmm38PYqTL4ClwB3la9t4jspfu5+VmecCZObmzLwf/+7Nm35t7N+9+TXguwz+7Vvwxt15SOCy\niPh6RLypXHZwZm4oX28ADi5fO7nccHq1cd0bgX8uX9vGw9mujSPiRGB9Zl7XKGsbD6/Xd/lI4NkR\ncXlETEXEU8rltvNwerXxqcCfRcRtwPuB08rltvHcHQ58PyI+FhH/ERF/ExF749+9+dSrjfdqlPHv\n3uh6trN/+wSj3W2pjWdm5vci4jHAmoi4qb4yMzMG38/WbO7ZbdfGmfkVgIg4HXgkM88fsL1tPLte\n3+PTgPqYzkG3VrON2+nVzouB/TPz6eVY/AuBx/fZ3naeXa82fifw1sy8OCJeDpxLMQSkF9t4sMXA\nk4E3Z+ZVEfFBis7ZNv7dG1m/Nv4j8O/ePOrVzmcCz8K/fQveWK88ZDkPRDnU4GKKS1gbIuIQKOaK\nAO4ui99BMY6uspzpy+fqo08bExGvB14E/Pdacdt4CD3a+FiKX2WujYhbKdrx6og4GNt4aH2+y+uB\ni8rlVwFbI+IgbOeh9GnjVZl5cVnkU0wPNbCN5249xa+yV5XvP0VxAnaXf/fmTb829u/e/OrVzsdQ\n5DP4t2+BG1vnoby8tW/5em+Knur1FJPIva4s9jrgM+XrS4BXRsSSiDicPpPLaVq/No6I44HfB07M\nzIdqm9jGc9Snja/MzIMz8/DMPJziP9knl8MSbOMhDPj/4jPAc8rlRwFLMvMebOc569PGNwC3RMSx\nZbHnAN8uX9vGc5SZdwG3l99VgOcB3wQ+h3/35kW/Nvbv3vzq085XZ+Yh/u3TOIctHQxcHBHVfv4+\nM78YEV8HLoyIk4F1wCvAyeWG1K+NbwaWUAxLAPhaZp5iGw+lZxs3ymxrQ9t4aP2+y7sD50bE9cAj\nwK+C7TykXm38hYi4F/hfEbEH8CDwa2Abj+AtwN9HxBLgO8AbgEX4d28+Ndv4jcBV+HdvvvX6Ltf5\nt2+BcpI4SZIkSa3siFu1SpIkSdoF2HmQJEmS1IqdB0mSJEmt2HmQJEmS1IqdB0mSJEmt2HmQJEmS\n1IqdB0mSJEmt2HmQJEmS1Mr/BUAdcI3g0vUxAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig, ax = plt.subplots(1,1,figsize=(15,4))\n", "t=0\n", "mesh=ax.pcolormesh(xx,zz,Ts[t,:,:])\n", "ax.set_title('Temp: t={} : {}'.format(t, run))\n", "fig.colorbar(mesh,ax=ax)\n", "ax.set_xlim([500,650])\n", "ax.set_ylim([-50,0])" ] }, { "cell_type": "code", "execution_count": 15, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "(-50, 0)" ] }, "execution_count": 15, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAwkAAAEKCAYAAABZkBkLAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XuYJVV57/Hfr7unZxhmABEFYbip4MloFNGMJEbtaEDk\nGImJBjXRKF5iOF6SGD0OnBxmNAHUSDReEhM10RwHQhRQoqCjsRWjiBhQcEQFGWG4KRJgcC493fOe\nP6qqu6q69u597V3d/f08z35676pVq95ae8/sWrvWW8sRIQAAAADIDA06AAAAAAD1QicBAAAAQAGd\nBAAAAAAFdBIAAAAAFNBJAAAAAFBAJwEAAABAAZ0EAIuC7fW2/3HQccwn2y+3fWUf6t1q+1kdbvvP\ntt+ePn+a7Rtz6x5j+zrbD9h+ne0Vti+zfZ/tf+1V/L1m+/dtf37QcQDAfKKTACwhtvfafmRp2Qbb\n/9Li9l2dlNp+h+170sd5XdQzZvu2/LKIODciXt1pnR3GMX1CvMhE+qhk+9By+1dtGxFXRsT/yK17\ni6QvRcR+EfF+SS+U9HBJB0bEabYPsf0Z27enn9Ujuj0Q22+0/WPbD9reYvuYduuIiE9ExLO7jQUA\nFhI6CQDmZUZF238k6VRJj08fv5Uuw8JziqTLm6x3g+VHStpSev3DiNibvt4r6XOSfrfrCCXZfpWk\n0yWdEhGrJP1PSfe0WcdwL2KZ77oBoFt0EgA0OqErFrJ/SdLfSfpV29tt39vmfv5Q0l9HxB0RcYek\nv5b08jbrkO19lZygHprG8YDtR+SviNg+Kv0l+uW2b7X9c9uvtf0rtr9r+79tv69U7+npL8332r4i\n/yu27b+xfbft+9PtH2v7NZJeIuktaRyfTsseavtTtn+a/oL9+lw9G2x/0vaFadzftv34Fo75cNsX\np3XeUxH7u9K4f2z75Nzy/W1/xPYdtrfZfrvtodz6V6fH/IDt79k+rmLfv5TWe1pu8SlKTuZl+4m2\n/yut40JJK3LbTl/xsf0fksYkvT9tr02S/kLSaenrV0TETyPi7yVd06Admh5PqeyQpLMl/UlE3ChJ\nEXFLRPx3k6bOrpb9p+3zbd8jaUP+Cprtv7P9rtI2n7b9p+nzVt7/f7F9v5J/EwBQS3QSABSkJ9C/\nVl4eEd+X9FpJ34iI1RFxYFr+rek2VY98R2KtpO/kXn9X0mObxPEd2y+qiOMXkk6WdEcax34Rcaeq\nr4isk/RoSS+S9F5JZ0p6Zrrf37P99HRfp0paL+n5kg6SdKWkC9J1z5b0NEnHRMT+SobI/Dwi/kHS\nJyS9I43j1PTE9DJJ10o6VNKzJP2J7ZNyMT1P0kWSHiJpk6RLbY80aYdhSf8u6RYlv7wflsWWeoqk\nGyU9VNI7JX0kt+6fJU1IepSkJ0o6SdKr0npfqOQk+qURsV8aV6HjZ/t4SVdIel1E/Gu6bFnaHptt\nj0q6VNLH0uP5NyVXAWa9FxHxzLRd/1faXi+RdI6kC9PX/9SoDVo5ngprlLTVL6cdxR+nJ+mtdIrX\nSbpZyVCovyqt2yRpusNk+yGSTpR0QRvv/7+ln6VNLcQCAANBJwFAQUQ8JCK+3mD1rBOsiDgv3abq\ncWCu6CpJ9+deP5AuaxTHEyLiwlbjaLDs7RExERGbJW2XtCki7kmvZFwpKfvl/LWSzo2IH6RDX86V\ndFx6NWFC0mpJv2R7KC1zV4P9/oqkgyLiLyNiMiJukfRhJZ2UzDURcXFETEk6X8kv7yc0agclJ6yP\nkPTmiNgZEbtL789PIuIjERGSPi7pEbYfbvtgSc+R9Kfpdj+T9J5cLK9S0sH5tiRFxM0RcWuu3mdI\n+rSSTsTncsufLum6tLN2gqSRiHhvRExFxKckfavJsZTby6p+32ZvNPfxlK1J/54o6XGSfkPSiyW9\nsoXd3RERH4iIvRGxq7Tua5LC9tPS1y+Q9PX0M9HK+//1iPiMJFXUDQC10fDXKwCL0pSkZaVlyyTt\nmYd9Pyhpv9zr/dNl/XR37vnOitdZJ+VISe+1/e7S9odGxJdtv1/SByQdaftiSX8eEdsr9nekkmFQ\n+SEtw5K+mnu9LXsSEWF7m5JOQCOHK+kI7G2wfrrDEhE70h/KVym5IrJM0p25H8+HJGUdgTVKfi2v\nYkl/JGk8Ir5aWjc91EjJr+W3l9b/pMmxSJ3nwBypJsdj+3uSsiFiJ0vakT5/Z0Q8IOkB2x9K4//w\nHPtqlJSdvWcXKulwXKlkyNnHczG2/P4DQJ1xJQFYWm6VdHRp2dGStra4/awTPNtnpmPKqx4P5Ip+\nTzO/3EvSEyTd0EbsTeNosKxVt0p6TekqyL4RcZUkRcT7IuLJSoZMHSvpzQ32eaukW0r17BcRz82V\nOTx7kg5PWSPpjiax3SbpCLef5HqbpN2SHpqLZf+I+OXc+kc32DaUdBKOtH1+ad1zNNNJuFPJkJ68\nI9uIsZ33rOnxRMRj02FLqyPiPyX9QMlVoE72OVeZCyS9wPaRSq70fCpdPtf73/SuUQBQJ3QSgKXl\nXyX9H9uH2R6y/ZuSnivpky1uf5ekNem4dElSRJyTOzkrP/JXDj4u6c/SxM7DJP2ZkjHmnbhb0kNt\n5+tvadhKSbbN30s60/ZaaTpB9oXp8yfbfkp6zDsk7VJyRSaLI39L2aslbbf9Ftv72B62/TjbT86V\neZLt56d5CH+S1ndVkxi/qeRk/DzbK53MLTArZ6QszdP4gqTzba9O3+9HZXkYSn5N/3PbxzvxaBdv\nObpdyS/yT7d9btoWR0taHhE/SMt8XdKk7TfYXmb7d5QMuWmmPNyouNJeoZnk5xXp61aOp3z8O5R8\n3t9ie5XtNZJerSS/oysRcZ2SuyR9WNIV6ZUKae73v5PPKAAMBJ0EYGl5m5ITu68pSVI9T9JLImL6\ntpTpFYCnNtj+P5RcEbjL9k/b2XFEfEhJUuf1SpKWL0uTfyvZvsH2ixvUdaOSX3N/7OSuPo/Q7F9p\nW/7FOCIulfQOSRemd525XlJ2X/z9JP2DkvbaquTkMLu7zUckrU2TtC9OhwQ9V8kVkx9L+lm6bdaZ\nCSXj/E9L6/t9Sb+T5idUB5jU+VtKfvW/Vckv6r+Xq698nPnXL5M0quS2o/cqSSw+JK33k0qScjcp\nyQ+5WEnycX7f9ysZ0/8c229TcgvRz+bW75H0O0ruUvXzNK5PqahZfFXx70jjCSUJ2b9o5XgaeJ2S\nIW13KPncf6KFBOlGbVpetklJEvx08nGL7z9XEgAsCE5y3fpQcXIbvvcoGY/54Yh4R192BAALhO2z\nJT06Il466Fg6Yfuzkt4XEVcMOhYAQH/15UpCOnb2/UouVa+V9GIn91gHgKVsoQ83GU8fAIBFrl/D\njdZJuikitqaXoy9UMtMqACxllcNNbB/RKPE7HUtfCxHxrsVw207bf9+gvT846NgAoC76dQvUw1S8\nhdw2JRP+AMCSFREbGyy/VclcDJgHEfFaJXNjAAAa6NeVBBKzAAAAgAWqX1cSblfuXuDp88IEMrbp\nSAAAAGBeRETt88LaPT/u5zH1q5NwjaRjbB+l5NZzpymZnbLoByGtmEwjmX0HwKHhZN3eqTTMydxc\nQruyZenrqiNZtbs6ugeXF7fNb5+Ntl1V8R5lMWZxVMQ8q0wWZzav7Ip82exvsq+hfZMJQvfev+9M\nmSyeLNZsXtz8wIQVaaz3pZ+T+3LrPrZBeuOG2XGuKL3Ob5PtK20m7VOxTfa+TR+ni9vmn2fbZcfb\n5D2fdk+68wNmKhxaXnw/l++TzJM0NTl7jqnJPcmyvb9YmSy4ybOPIbtxYrn9p2Y+TPusSpbt3jkq\nSVr9kGSS3fvvyd0p8i/OlV65QToqiW90RfJ3OD2+nXflyl6b1vucZELW1fsl9f30xzN96ocelcyr\nNTqUHN+db0vnPvv1mWq+88zHSJJO10clSd8+MLljaTw+OU4fMfP5ffjHb51+frfyt8FfODZs2KAN\nGzYMOoxFj3buP9q4/2jj/luIbZybqb32/rLFcv+nr1H0qZMQEZO2Xyfp80pugfqRiPh+P/aFmpqc\nuwgAAACKls1dpCHbW5XMNTMlaU9ErOu0rn5dSVBEXC7p8n7VDwAAACw2XZ6ch6SxiLh3wHGg9p4w\nNugIFr8njg06gkVvbGxs0CEsCbRz/9HG/Ucb9x9t3F/7zF1kLj0ZW9WvuxuhLo4bG3QEi9/xY4OO\nYNHjC2l+0M79Rxv3H23cf7Rxfy1r8dFASPqi7Wtsv7qbOBbulYRszHuW2FtOjm2njvLzObfLOmhN\nEph3lQJZlS3PLSsVmU5YHsklTY+k+8oSn2fn6M7I2qBZMncjq3LP72tYam7ZPgeUk5AlLGdJzTuz\nJORDKpLVu5Al1UvS3lXLm5TUTKK2JD1s4f6TAwAA/dflmcJTI+JO2w+TtNn2jRFx5QDiAAAAANAr\nja4SbJE0112AIuLO9O/PbF8iaZ0kOgkAAADAQtbo5Pzx6SNzSWm97ZWShiNiu+19JZ0kaWOv4wAA\nAAAwz7q4BerBki5J54QYkfSJiPhCp5UNvpNQMRFWZm950rL8OP/t6d9yS+bHf5frzra/p2JnB6V/\nswnYypOhSdKD6fPyRGT5nIT70vHp5ZZN8wxG12yfXjSxKy2bTu42etADyd8VEzO73HZQ8/gkjaYT\nfk1s2y9ZcFS6ImuKfI7BQSrK8h8ezCXCZ5OMZfkTVZ+Scs5FJp8PkbXTdN5IaYK6ByvG8u8qv57Z\nz95S7sfOu/Yt7jO3/U6l67Jci6xMLvdiaP9fJPVmE65ly3P5BtkkatkEa9v/O5nFbu/uXOwjs7eT\npJ33rZ51DDog+fPQ/X4uSVqu5P2874CZz0U2idpKJe+rnlzcjyTtUBLzb6e/I3z735PJ1KZzV559\n90zhjwsAACwQnXYSIuIWScf1Ko7BdxKwOJU7UgAAAJhTD26B2hN0EgAAAICaqMvJeV3iAAAAAJa8\nLnISeopOAgAAAFATdTk5H2wcI1E9EVl5WZakm88LzbpZ2RGkSbFZEm/edILw5EixnvzRZ4nFK5JE\n0mxSrr2FuEqJ1NnffIL0AbsrY98nl5iayfY1mq7bsT0ZhTY8MlNfObk226ZwfA8m64YOLiXirsja\nseJtLrfBqvwEbtlxZYna6brJXHLzSGlZub78hHArXKx3+m9F0no2eVzVJ7NJkvsse7JtmsQ1Hd9E\n4eXIspnP38xnJ913lpxc+IwWg80+O5Wf27To9okkg3pidLS4H0nbdyQJz/dNHjB7ewAAsKhxJQEA\nAABAQV1OzusSBwAAALDkcSUBAAAAQAG3QJU0tO9M/kB+HPhwKSdhekKr/MKRNPQslyAd15+fiGx4\nZLJQ387pMeLZtjPVlbcfTe/zP5GLa3qseSnm/PJsWbbPmb9JLFOTjZt8+T4Ts5ZNt0sppyHfRrPy\nHdJ6snYbyeVpTOdapJOBDS3fXdxPru7dy0cL68rHL81MMqYVVfkFxZyN8r5mH61m533klCcrm/48\nFHIV0pyGbJ6G8gRuuc9H1t5Tk8V9ja6YHVl27NNtkY9jVXFdZmI69yL3nmcfvdLx5Y8t+6ysXJm8\nbztXPSRZcedM+cfoB5Kkt+kvkgV/kK44OCtxb672fWcdDwAAqCeuJGBxq0pIBwAAQFN1OTmvSxwA\nAADAkres1bPzPt/9kE4CAAAAUBMjdBIAAAAA5C1rY1qofhpoJ2Hl6p3TSZp55aTObAKxiV0zCaXT\nSblpsujqNHl3eCiXgKupQn1TaYLrhNLJxnL7qdpekpbnJi+bbDCZV1UycnZcy0cnCrFMjc7UMbW3\nur5G+8krt1F+u2xd1cRr2bod6QRsWZLy8oqy2fbZNlVxZceetW32t5x8nl+WtU0+Qbid7TMTy2aX\nyU9Klj+GqjpWlibey44ve8+K2xXjyk94N7Usa+8sETpZl03sNzFSjEmq+Jzlktan378stTudU013\nzZTfrdFZdRZtzT1/7BxlAQBAXbR8JaHPahIGAAAAgGWzf1scCDoJAAAAQF3U5Oy8JmEAAAAAqMvZ\n+UDDWL1ye+XybPx+JhvHnx8rProizUlIx2/vM7Qz3XYmx2EkrWd6/Hc6tntHxZj9qu0laWpopony\n+QQFTYaHj6o4Jn4q3+RD1dsU8haUjoXPchpKrwvrRocrX+dl64b3K+V9VExt1kp95WWN8iykmfeh\nWX5GOa+iaT0V+SzlidCyMll9+VyA8jGXj1eSsvSV8nFV5WeUczeqc1Wqjyu/fKad0uNbUbVFyS0b\nJUkbb2mhLAAAqK+adBIanKYCAAAAmHfDLT4asD1s+1rbl3UTRk36KgAAAAB6cHb+RklbJK3uphKu\nJAAAAAB1sbzFRwXbaySdIunDktxNGFxJAAAAAOqiu7Pzv5H0Zkn7DTaMLo1q93RycSsmcwOwRksT\nXi1PE4SrEnCnk3WHkkTQLGk3X3aldhTKZqY0O4m4WVzl4ynXp6oE4SYDyyYbrGun3aqUE7Sb1VdO\nms4rx5dP9J5rn1XbVSYPN6pvdHaZciL19HtfUd/yUlJ5dizLc+9RFvPE0PJCPROjo7kypfc8jWFq\nZHaC9XQCdWmb/GR205MAZmU6+le6s5ONAADAoDX43h+/Xxp/oPFmtp8r6acRca3tsT6FAQAAAGDe\nNfjteOzA5JHZePusIr8m6Xm2T1Fyb8T9bH88Il7WSRjkJAAAAAB1MdLioyQizoyIwyPiaEkvkvQf\nnXYQsjAAAAAA1EHvzs6jm40H2klY2ea46eHC9ZfSpFkqTUClmTHd5XHzIxXj1LNJz8pj8ycLOQnV\nzdVkLrWm4/DnqneuurtRrjcfZzmebF1VnFk9M/kKs3MuWjFVeq9aybmoytdoFHtVfeX8lezzlf9c\nLG9wPPs0y69IY5gaSuvL5UlkdWeft6xsfpK3WXkUI9m/8WY3KUiuP74kfbVBe5qUBQAAtdWDs/OI\n+Iqkrww4DAAAAAA90eD2pvONTgIAAABQFzU5O+84cdn2C21/z/aU7eNL69bb/pHtG22f1H2YAAAA\nwBIw3OKjz7rpq1wv6fmSPpRfaHutpNMkrZV0mKQv2j42IvZ2sS8AAABg8avJlYSOw4iIGyXJnpVM\neaqkCyJij6Sttm+StE7SVeWCw5rseDK1sqpk5NlKyaIVSc7Ntu9mArOqeqfaSNJtRaOk36rJ3hol\nczeLp2r55HSybzG5uSoRupzEnU8ynknWbTyZ2lRp3XDF52GqtF2zydQaJZVXLS8nMDebAK9cT/E4\nJyvL5I9lZlmWuJz+nWznn+tpuee3trEdAAAYqJp0EvoxT8KhkrblXm9TckUBAAAAQDMLYbiR7c2S\nDqlYdWZEXNbGfirv03r3hn/QULpq1djxWj12fFUxAAAAoGXj4+MaHx8fdBidqcmVhKZhRMSJHdR5\nu6TDc6/XpMtmOXjDa3o21AYAAACQpLGxMY2NjU2/3rhx4+CCadeKQQeQ6FVfJZ+Y8BlJm2yfr2SY\n0TGSrq7aqNFEVZ1oNva8rKpj0s72ncTVSplWxrn3al/tbNOOdjp9zSaaa2USuuo8g+r9V9XXTu5F\nppwX0Uw5X6Oq7sk26gMAAEvAPAwlakXHnQTbz5f0t5IOkvRZ29dGxHMiYovtiyRtkTQp6YyI6Gpa\naAAAAGBJWAjDjZqJiEskXdJg3TmSzum0bgAAAGBJWuidBAAAAAA9ttCHG/VCu+OwhzU1a9x+r8bY\nN8pJqNpnN/VLjce1dzsuvVGcrY7d71VuRFmjuRmqdJunMFX6SLeyTTmuSQ33PaE+22c53qZ2zTzd\nqZVNi/5Q1q/HTzoJDQAADFJNfsKvSRit6fXJ6yD2OYhjaAWJszMWwx236CAAALBA1eTsvCZhAAAA\nANDyQQeQoJMAAAAA1EVNzs5rEgYAAACAupyd1ySM1jQbN9/punbKNCqbzzOYq55mCcyDzFdod9/t\nlG8lYXmm3uQjWZVw3CjJt6psJwnLVcuzfbZyvN28f8XP0GRx2WT2d+56yEUAAGCB6/B0wvYKSV9R\nMmBpVNKnI2J9p2EsqE4CAAAAsKh1eHYeEbts/0ZE7LA9Iulrtn89Ir42j2EAAAAA6Lkuzs4jYkf6\ndFTJNYl7O61rqPMwAAAAAPTUcIuPCraHbF8n6W5JX46ILZ2GwZUEAAAAoC5WVC8e35I8momIvZKO\ns72/pM/bHouI8U7CWFKdhLomCHcbz1zbVyVWN5u5upV4OklGbkfVNo1mqq4u25uPdrktqtqmHE8r\nbdMsIXpmXVrvpOes7+Fx6pxlAADAAtDgFGbs8ckjs/FTjauIiPttf1bSkyWNdxIGw40AAACAuuhw\nuJHtg2wfkD7fR9KJkq7tNIwldSUBAAAAqLXOz84fIeljtoeUXAj4l4j40vyHAQAAAKC3Or8F6vWS\njh9wGL3Xyljvdrfvdv+dlO2knmaTs3UTT9XY+FYmEGu8n84+Lp1MRFb13mdlGrVbp9rJz2jWTpPT\n8c3dTuVjaUluMrV9lNzhbKuObn17AABQf4NLnS2oTScBAAAAWPIa3N1ovtFJAAAAAOqCKwkAAAAA\nCmpydl6TMAAAAADU5ex8oGG0MwlYL3Sb6FpOmO1XsnSzevuZYNy4vu4mgGtn4rVm22bHNZxm8FbV\nOzLHe9xuLOXyzdq22QRps8umict7K8oOlfY9ObvIIT+7X5L084c9dM59AQCABYROAgAAAIACchIA\nAAAAFNTk7LwmYQAAAADQ8kEHkKCTAAAAANRFTc7Oa5O43EpCabOk1E6SWDtVTkztNtG4WZlyomwn\nsxd3q534mmlvNurGyeHNZlzuRdJ28+TkuWfbblZmOuk9TVienEz+jow0OZYscTk3ucpnHnZSw30A\nAIAFjE4CAAAAgIKanJ3XJAwAAAAAwd2NAAAAAORN1eTsvCZhAAAAAKCTIGl3g3s8NUpMbTcptdEM\nve0mQM9Vf/W63lwrKtezEGZcbjQjdeXswo3qGKpI5E23z9ZVxdDOrNqdHEMr9VUdZ6OYpyZHpv8O\nj6QZytnbN1Sq5ARJNyRPX61/bDkuAACwcOxePtpiyYm+xlGTvgoWi05OruumnY5GS/VVdHhmlck6\nCM3c0INgAABArU0N1+McqvxbJQAAAIABmdJwS48y24fb/rLt79m+wfYbuomDKwkAAABATbQz9L1k\nj6Q/jYjrbK+S9G3bmyPi+51URicBAAAAqIlOc0sj4i5Jd6XPH7T9fUmHSuqok9DxcCPb77L9fdvf\nsX2x7f1z69bb/pHtG20zNSwAAADQgk6HG+XZPkrSEyV9s9M4uslJ+IKkx0bEEyT9UNL6NKi1kk6T\ntFbSyZI+aJvcBwAAAGAO3XYS0qFGn5T0xoh4sNM4Oh5uFBGbcy+/Kel30+enSrogIvZI2mr7Jknr\nJF1VrmNCM7d4anf8Vfk2plXbZ2WGS7dUrbqMM337yg7GgTWLvVe3I20lrm5uOdqKycm56xsZmXlf\nyuWz231m8nf0ydZly5q1aXldfp+9uKtS23WkXeBye08f/8jssnO1jSRNjaZlsmZa1V5YAABg4dmt\n6lugXj2+S1eP72q6re1lkj4l6f9FxKXdxNGrnITTJV2QPj9UxQ7BNkmH9Wg/WCBa6VAAAACgqNEP\nzE8aW6Unjc38YviBjfcX1tu2pI9I2hIR7+k2jqadBNubJR1SserMiLgsLXOWpImI2NSkqug8RAAA\nAGBp6GJUxFMl/YGk79q+Nl22PiKu6KSypp2EiDix2XrbL5d0iqRn5RbfLunw3Os16bJZ7t3wgenn\no2MnaMXYuubRAgAAAHMYHx/X+Pj4oMPoSKedhIj4mno4B5ojOvuR3/bJkt4t6RkRcU9u+VpJm5Tk\nIRwm6YuSHh2lHdmOR8XMFLL9zUko/q1CTsLc2slJqCrbTk5CO/I5CYOQzajcKCchH19WdvdEMt6w\nKhcha4OVozslSXf+19HJim0zZR7+vFs7ivVuHdHRdgAALGS2FREedBxzsR1Xx+NaKrvON/T1mLo5\ng32fpFFJm5MhUPpGRJwREVtsXyRpi5KUyzPKHQQAAAAAs/XqB+ZudXN3o2OarDtH0jmd1g0AAAAs\nRb24U2Mv1KOrAgAAAKAwRcAg0UkAAAAAaqLdPN1+GWgnYbtWVy4vJ4BmyZ7tKicsD6ezUmUJzfk3\nYULLk3336I1pJWm4KrG1UZlZ9VckvM4u0/qxDLeQ/NtKfRMt77G6vmxZVTyN1rXSFlkycDtlq1Ql\nWDf6x1y5r5HiumbHOy29JXKnycoAAGDhWPA5CQAAAAB6i5wEAAAAAAV0EgAAAAAUkJMAAAAAoCDL\nkx00OgkAAABATTDcCAAAAEABw40AAAAAFHALVAAAAAAFDDcCAAAAUEAnAQAAAEABnQQAAAAABbu5\nBaq0/YHVlcunJpMe1PDI1Jx1lMtk27ZStmq7TrdvZf9lk3uSsiPL5q43KzsfWomn35q1Y3ldr9+X\nVsrmy8y1/6nJ2f/Mmn3ehkcmJUmj2i1Jevixt84ZDwAAWBzqciVhaNABAAAAAEhMabilR5ntj9q+\n2/b1vYiDTgIAAABQE5MabulR4Z8kndyrOMhJAAAAAGqi03kSIuJK20f1Kg46CQAAAEBN1CUngU4C\nAAAAUBONOgk/Gd+qn4z/ZN7ioJMAAAAA1MRujVYuP2TsWB0yduz0669t/Gpf46CTAAAAANREpzkJ\nvcbdjQAAAICa6OIWqBdI+rqkY23fZvsV3cRRj64KAAAAgI4TlyPixb2Mg04CAAAAUBMN5kCYd3QS\nAAAAgJqoS05CPaIAAAAAwDwJAAAAAIomGtwCdb7RSQAAAABqgpwEAAAAAAXkJAAAAAAoICcBAAAA\nQAGdBEk7H1zZ3x1M1qOR5zKxu8cVjkx1tfnEVPKxGBqe7EU0bRlZ1n7sU7n3ebjLY59vWbyjK2Y+\nBMtHJyRJO9Xnfx8AAKB2yEkAAAAAUEBOAgAAAIACboEKAAAAoKAuw42GOt3Q9tttf8f2dba/ZPvw\n3Lr1tn9k+0bbJ/UmVAAAAGBxm9JIS49+67iTIOmdEfGEiDhO0qWSzpYk22slnSZpraSTJX3Qdjf7\nAQAAAJa7G8u6AAAONElEQVSEKQ239Oi3jk/eI2J77uUqSfekz0+VdEFE7ImIrZJukrSu4wgBAACA\nJaIunYSurlXY/itJL5W0UzMdgUMlXZUrtk3SYd3sBwAAAFgKFsQ8CbY3SzqkYtWZEXFZRJwl6Szb\nb5X0HkmvaFBVdBcmAAAAsPjt1vJBhyBpjk5CRJzYYj2bJH0ufX67pMNz69aky2Z799tnnv/q06Vf\ne0aLuwMAAACqjY+Pa3x8fNBhdKSbKwm2T1byw/2wpA9HxDs6riuisx/5bR8TET9Kn79e0rqIeGma\nuLxJyfCjwyR9UdKjo7Qj26Hbd3Uad2sWyIzLPdejWYcXyozLeYOYcbmdfQ6PJG06sav4K0HVjMv9\ncLeO6FvdAADUlW1FhAcdx1xsx6PihpbK3uzHFY7J9rCkH0j6TSU/0H9L0osj4vudxNJNTsK5th8j\naUrSzZL+WJIiYovtiyRtkTQp6YxyBwEAAADAbF3Mk7BO0k3pjYNk+0IlNxSa305CRLygybpzJJ3T\nad0AAADAUtTFHAiHSbot93qbpKd0WhkzLgMAAAA10SgnYff4VZoYv6pyXaqnI3foJAAAAAA10aiT\nMDL2VI2MPXX69YMb/7ZcpHzzoMOVXE3oCJ0EAAAAoCZ2T4x2uuk1ko6xfZSkOySdJunFnVZGJwEA\nAACoianJzk7PI2LS9uskfV7JLVA/0umdjSQ6CQAAAEBtTHVxC/+IuFzS5b2Ig04CAAAAUBPddBJ6\niU4CAAAAUBOTe+gkAAAAAMjZO1WP0/N6RAEAAABAYrgRAAAAgIJd9Tg9r0cUAAAAAKTJQQeQoJMA\nAAAA1AWdBAAAAAAFdBIAAAAAFOwZdAAJOgkAAABAXUwNOoAEnQQAAACgLhhuJGnXaOfbjjToZlXd\nW3bSXewnOt+2FxodZx8NDc/Pp3NkWXvHNjyAtijvO5sqPR/L8Mhkuq76n1O2XpKWj05Ulmm0HAAA\nLDG7Bh1AgisJAAAAQF1wJQEAAABAAZ0EAAAAAAV0EgAAAAAUcAtUAAAAAAU1uQXq0KADAAAAAJCa\nbPHRBtsvtP0921O2j29lG64kAAAAAHXRn1ugXi/p+ZI+1OoGdBIAAACAuuhD4nJE3ChJdutzh9FJ\nAAAAAOqCuxsBAAAAKOiwk2B7s6RDKladGRGXtVsfnQQAAACgLhrdAvXmcenH4w03i4gTexkGnQQA\nAACgLhrdAvWoseSR+dLGTvfQUmICt0AFAAAA6mJXi4822H6+7dsknSDps7Yvn2sbriQAAAAAddGf\nuxtdIumSdrahkwAAAADURaOchHlGJwEAAACoi0Y5CfOMTgIAAABQF8yTAAAAAKCATgIAAACAgprk\nJHR9C1Tbb7K91/aBuWXrbf/I9o22T+p2HwAAAMCSsLvFR591dSXB9uGSTpT0k9yytZJOk7RW0mGS\nvmj72IjY282+AAAAgEWvJsONur2ScL6kt5SWnSrpgojYExFbJd0kaV2X+wEAAAAWvz0tPvqs406C\n7VMlbYuI75ZWHSppW+71NiVXFAAAAAA0M9Xio8+aDjeyvVnSIRWrzpK0XlI+38BNqorKpe/dMPP8\nKWPSCWPNwgEAAADmND4+rvHx8UGH0ZmaDDdyRPX5e9ON7MdJ+pKkHemiNZJul/QUSa+QpIg4Ly17\nhaSzI+KbpTpCN3eRpjDSoAs1OVyxrFn/Za79tN8+PdXoOHu9Tc7Q8Px8OkeWtRfncJfH1Y1s31Pp\n5ysfy/DIZLquus+drZek5aMTkqTdE6OFMtnyfrtbR8zLfgAAqBPbioguTgjnh+3Qc1o897y8v8fU\nUeJyRNwg6eDste1bJD0pIu61/RlJm2yfr2SY0TGSru5FsAAAAMCiVpNboPZqnoTpLk9EbLF9kaQt\nSi6YnBGdXK4AAAAAlpp5uL1pK3rSSYiIR5ZenyPpnF7UDQAAACwZNclJYMZlAAAAoC4W2XAjAAAA\nAN0a3H1aCugkAAAAAHXBcCMAAAAABXQSAAAAABT0ISfB9rskPVfShKSbJb0iIu5vts1Q78MAAAAA\n0JHdLT7a8wVJj42IJ0j6oaT1c21AJwEAAACoi8kWH22IiM0RsTd9+U1Ja+bahuFGAAAAQF30/xao\np0u6YK5CdBIAAACAumh0C9Q949LkeMPNbG+WdEjFqjMj4rK0zFmSJiJi01xh0EkAAAAA6qLRUCKP\nScvGZl7v2lhYHREnNqvW9sslnSLpWa2EQScBAAAAqIs+3ALV9smS3izpGRGxq5Vt6CQAAAAAddGf\nnIT3SRqVtNm2JH0jIs5otgGdBAAAAKAu+nAlISKOaXcbboEKAAAAoIBOAgAAAIAChhuVjcQ87qvR\nPa66LJsaGu7D9aqSkWXtx1VluIPj68W2vdzX8Eh1ey8fnWhpGQAAQF3QSQAAAABqo/+zqbWCTgIA\nAABQG/0fCdIKOgkAAABAbXAlAQAAAEDBzkEHIIlOAgAAAFAjXEkAAAAAUEBOAgAAAIACriQAAAAA\nKOBKAgAAAIACriQAAAAAKODuRgAAAAAKGG4EAAAAoIDhRgAAAAAKuJIAAAAAoIArCQAAAAAKen8l\nwfbbJT1PUkj6uaSXR8RtzbYZ6nkUAAAAADq0p8VHW94ZEU+IiOMkXSrp7Lk24EoCAAAAUBu9vwVq\nRGzPvVwl6Z65tqGTAAAAANRGf3ISbP+VpJdK2iHphLnKM9wIAAAAqI3JFh9Ftjfbvr7i8VuSFBFn\nRcQRkv5Z0t/MFUXHVxJsb5D0Kkk/SxedGRGXp+vWSzpd0pSkN0TEFzrdDwAAALB0NLqS8MP0US0i\nTmxxB5skfW6uQt1cSQhJ50fEE9NH1kFYK+k0SWslnSzpg7a5YjEg8Z9fHXQIi97kV74+6BAWvfHx\n8UGHsCTQzv1HG/cfbdx/tHG/Nbpy8Eglp9bZo3W2j8m9PFXStXNt0+3JuyuWnSrpgojYExFbJd0k\naV2X+0GH4utXDjqERW/yq3QS+o0vpPlBO/cfbdx/tHH/0cb91pe7G52bDj26TtKYpDfNtUG3icuv\nt/0ySddIelNE3CfpUElX5cpsk3RYl/sBAAAAloDez5MQES9od5umnQTbmyUdUrHqLEl/J+lt6eu3\nS3q3pFc2iq3dwAAAAIClp/e3QO2EI7o/f7d9lKTLIuKXbb9VkiLivHTdFZLOjohvlrah4wAAAIB5\nERFVw+Rrpd3z434eUzd3N3pERNyZvny+pOvT55+RtMn2+UqGGR0j6ery9gvhjQIAAADmS53Oj7vJ\nSXiH7eOUDCW6RdIfSVJEbLF9kaQtSgZVnRG9uFwBAAAAYF70ZLgRAAAAgMWjr/MX2N5q+7u2r7V9\ndbrswHRGuB/a/oLtA3Ll19v+ke0bbZ/Uz9gWiwZt/C7b37f9HdsX294/V542blNVG+fWvcn2XtsH\n5pbRxh1o1M62X59+nm+w/Y7cctq5TQ3+v1hn++p02bds/0quPG3cJtsH2P5k+pndYvspfO/1VkUb\nn8D3Xu9VtXNuHd99S0FE9O2hZBjSgaVl75T0lvT5/5Z0Xvp8raTrJC2TdJSS+RWG+hnfYng0aOMT\ns7aTdB5t3Ps2TpcfLumK/HrauLftLOk3JG2WtCx9/TDauedtPC7p2enz50j6Mm3cVRt/TNLp6fMR\nSfvzvTcvbcz33jy0c/qc774l8piPmZDLCRjPSz94Sv/+dvqcSdg6V2jjiNgcEXvTl9+UtCZ9Tht3\nriqR6HxJbykto427U27nP5Z0bkTskaSI+Fm6nHbuXLmN71RykiVJB0i6PX1OG7cp/fX6aRHxUUmK\niMmIuF987/VMozbme6+3mnyWJb77lox+dxJC0hdtX2P71emygyPi7vT53ZIOTp8fqmTitQyTsLWm\nqo3zTpf0ufQ5bdyZWW1s+1RJ2yLiu6WytHHnqj7Lx0h6uu2rbI/bfnK6nHbuTFUbv1XSu23fKuld\nktany2nj9h0t6We2/8n2f9n+R9v7iu+9Xqpq45WlMnzvda+ynfnuW1q6nXF5Lk+NiDttP0zSZts3\n5ldGRLj5/WDJqp7brDaOiCslyfZZkiYiYlOT7WnjuVV9jtdLyo+5bHbLMtq4NVXtPCLpIRFxQjpW\n/iJJj2ywPe08t6o2PlvSGyLiEtsvlPRRJUM3qtDGzY1IOl7S6yLiW7bfo6QTNo3vva41auP/K/G9\n10NV7bxR0tPEd9+S0dcrCZHOo5AOEbhEyaWnu20fIiVzLUj6aVr8diXj3DJrNHPZGw00aGPZfrmk\nUyT9fq44bdyBijZ+hpJfWb5j+xYl7fht2weLNu5Yg8/yNkkXp8u/JWmv7YNEO3ekQRuvi4hL0iKf\n1MwQAdq4fduU/Mr6rfT1J5WcaN3F917PNGpjvvd6q6qdn6gk34DvviWib52E9LLU6vT5vkp6ntcr\nmWztD9Nifyjp0vT5ZyS9yPao7aPVYBI2zGjUxrZPlvRmSadGxK7cJrRxmxq08dURcXBEHB0RRyv5\nz/T4dDgBbdyBJv9fXCrpmenyYyWNRsQ9op3b1qCNb5B0k+1npMWeKemH6XPauE0RcZek29LPqiT9\npqTvSbpMfO/1RKM25nuvtxq087cj4hC++5aOfg43OljSJbaz/XwiIr5g+xpJF9l+paStkn5PYhK2\nDjVq4x9JGlUynECSvhERZ9DGHals41KZ6TakjTvW6LO8TNJHbV8vaULSyyTauUNVbfx52z+X9AHb\nyyXtlPQaiTbuwuslfcL2qKSbJb1C0rD43uulchufLulb4nuv16o+y3l89y1yTKYGAAAAoGA+boEK\nAAAAYAGhkwAAAACggE4CAAAAgAI6CQAAAAAK6CQAAAAAKKCTAAAAAKCATgIAAACAAjoJAAAAAAr+\nP2QjgdqWWo/TAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig, ax = plt.subplots(1,1,figsize=(15,4))\n", "t=0\n", "mesh=ax.pcolormesh(xx,zz,Us[t,:,:])\n", "ax.set_title('U: t={} : {}'.format(t, run))\n", "fig.colorbar(mesh,ax=ax)\n", "ax.set_xlim([500,650])\n", "ax.set_ylim([-50,0])" ] }, { "cell_type": "code", "execution_count": 16, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "(-50, 0)" ] }, "execution_count": 16, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAx8AAAEKCAYAAABzIOikAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xu8ZFV95/3Pt0/TgEJslci1HTS0GUmiYBIgGrXjLYhR\ndJJ4ySQS9VESH3TmlUxEZJ6JxAgYI1Hj6PAoZsgFCfGWNkKAMGljLoIYRAwQLrEjjdBeUeTWfU7/\n5o/a1V1dXVWnTp1bnerP+/XanL3XXnvttVYVffY6e11SVUiSJEnSYlu13BmQJEmStHew8SFJkiRp\nSdj4kCRJkrQkbHxIkiRJWhI2PiRJkiQtCRsfkiRJkpaEjQ9JEyPJvUmOXO58LKUkm5M8e4HT/NUk\nnx3x2iOT7Eiyqjm+NMmvdJz/3STfSPK15vglSe5oPrsnL0wJFl6SLyd5xnLnQ5JWOhsf0oRLckaS\nS7vCbu0T9tIh0tuc5Fkj5uWYJF9Icl+Sa+fzsJlkU5LXdIZV1YFVtXnUNEfIw24P2sukmm1JJTk/\nyWtni1dVJ1XVnzTXPBb4DeA/VtVhTZTfB17ffHbXJzmt+W48mOSPFiCfj0/yV0m+1zR63jFKOlX1\no1X1d/PNjyTt7Wx8SJPvM8BTkwQgyaHAauCYjr9OHwr8EDDMw1UBmWsmkqwB/hL4Y2AtcCHwl0n2\nmWtaHfkYF3OujwlwIvDpOV7zWOBbVfUtgOY7+Vjgxo44dwJvAz483ww237krgb8BDgYOB/50jmms\nnm8+liNtSRpXNj6kyXctsA9wTHP8dOBvgVu6wm6vqrsHJZTkT2g9LH6q6Sbz3+aQjw3AVFW9p6q2\nV9Uf0npon/NblCRvb/L8viYf723CdyR5fLP/v5O8v+n2c2+SzyY5JMl7knwnyU1JjulI87AkH0vy\n9ST/luQNHeeOa/4a/90kdyf5/eZUu7F2T3OP45v4r05yY5JvJ/nr5i/+dOTxDUlub/4S/3vthuEs\nZX5tk+b3kvxLZ96BY5Ncn+SeJBcn2bfjup9L8sWmzP+Q5Mc6zq1L8vGmzN9M8od97v3Opv4ObI6f\nBNxTVV9LMpXk95uy3A68oOvaTUlek1bXsCuAw5q6ugj4HjAFXJ/kVoCq+kRV/SXwrT556VueHn4V\n2FJV766qB6pqW1XdMCB++x6bk7wpyZeAe5sybk7yrOZ7cn+SR3bEP7Yp/1RzPNvn//qmvP86W14k\nadLY+JAmXFVtA64GntkEPQP4LPD3zX477DMASd6c5FN90voV4KvAzzXdZH6/ueae5mGw1/am5vIf\nAb7UleT1TfgekvxSkuv75OPMpgz/b5OPN/Yp/i8CZwIHAduAzwGfBx4FfBQ4r7nXKuBTwHXAYcCz\ngf+a5HlNOu8B/qCqHgE8HviLJvzpzc9HNPm4OsnJwBnAS5r7fhb4SFe+Xgz8OPAU4GTg1X3y366L\nXwR+G/iVqvoB4EXAt9unm3L+LPA44Em0HrpJcixwAfDapsznAxuT7NM8KP8V8BXgP9B6K7BbPtPy\nQeBHgedW1b3NqZOaa2nSfgGthuxPAL/A7m+lCqiqugp4PvC1pq5+qaoObOI8qarWdxe7Rz30K8+a\nPlV3AvDvTQP0G0n+NsmP9onb7eVNftdW1Uy7TFX1NeCfgJ/viPtLwF9U1cyQn//JwE8CRw+ZF0ma\nGDY+pL3DZ9jV0PhpWn+x/2xH2NObOFTVuVX1wrkkXlVrq+qRfbbfa6IdAHy369LvAQfSQ1VdVFWz\njQkZ9MaggI9X1XVV9RDwCeC+qvrTqirgEuDYJu5PAgdV1e9W1XRVfQX4EK0HUGg1XNYnOaiq7q+q\nqwfc/9eAc6rqX6tqB3AOrS5u6zrivKOq7qmqO4B3A6+YpZz/T3PNFwCq6vaq+mpHOd9bVXdX1Xdo\nNaLab0VeB5xfVZ+vlj8GHgJ+CjgOOBT4reatwENV9Y8d99wHuJhWF7kXVtWDHedOAtpjhl5Kq2F2\nZ3P/s/vUCwPCe+nVra5feU7ok8YRtD7D99Aq66cZrqtfu07vbL473S6i+cyat1Yva8JguM//nObz\n75W2JE00Gx/S3uHvgJ9uuor8YFXdTuuvt09twn6E4cZ7zMe9wA90hT2CVgNkVLON+/h6x/6DXccP\n0GoQQesv/4d1vrGh9dfrxzTnXwM8AbgpyTVJduta1OU/AO/pSKfdfejwjjh3dOx/ldbblkGOAG4f\ncL6zu1x3uX6zq1xH0HoQXwf8e/OA3MtRwAuB36mq6XZgkrXAfwTaDZVDe5RnIfRqqPQtT/Om7N5m\na49FeQD4bFVd3jQqfx94dJP/2dwx4NzHgZ9KcgitBvyOqvr7jjzO5fOXpL2Kg92kvcPnaD3ovxb4\nB4Cq+l5a052+jlZXmH8fMq09HviTfL9XeOPtVXUurUHFv9l17klAz3EGo+RjHu4AvlJVT+h5o6rb\naHWtIcnPAx9N8qg+efgq8Laq6u5q0+mxwE0d+3cOkb+jZomzW5Y78vL2qjq7O0KSnwIem2Sq6VbU\n7SbgfwKXJXlWVd3ShP8scFXz9gjgrqYMbZ3789GvbnuWp3FR1/H1wNPaB8OMrZnl/q0TVd9JcgWt\nNx5Hs3u3qmE+/3GaLEGSlpRvPqS9QFU9QGvg+W+w+xuOv2/CPjOH5LbSmhmrM/0Dmn78vbZzm2h/\nC8wkeWOSfZO8EdgB/J8Ri7VHPrrM5UHzGloDi9+UZP9mgPGPJvkJgCS/nOQHm7jfpfXwuAP4RvOz\nMx//C3hLkqObax/RjNno9N+SrG264rwR+PNZ8veh5pqnNOMwjuocxNxDu+wfBH4trQHzSfLwJC9I\ncgCtcUB3AecmeViS/ZI8tTORqroYeAvwN2kG8tPqctU5y9UlwBuTHN68RXvzLGUZqKn7/Wj9cWyq\n+a5MDVGeXv4UOCHJs5s0/iutz+ymPvHn4iLgFFpjPzobPcN8/pK017LxIe09PgP8IK0GR9tnaQ2K\n3dkgSfKWdK0B0uUc4L833Up+Y9ibV9V2WgOtXwl8p/n54s4uPZ2S/OckXx6Q5HuAX0hrRqF397ol\nPQY+94hD85f/n6M1VuLfaD2g/v/s6ib2s8CXk9wL/AHw8maMxP3A24F/aOrjuKr6JPAO4OIk3wVu\naK7v9JfAF2gNcP8rZplWtqo+2tynPUPUx4FH9oveUa4v0Hrb9T5aA9RvpVXvNN2tXkjrjcpXab1d\neWmPNP4Y+B3gqrQWcHwe8Ncd9/sgcDmttwzXAh9j8F/2e34GHf4/4H7gdOCXaXWdOnO28vS8Uett\nzS/TahB8uynvi/p95+ZoI626u6tzBq0hPn/fekjaq2XXm/MFTjg5kdZAyingQ1U10sJOkjRJkuwA\njqqqf1vuvMxVkuNoDcTuN8BbkqSBFuXNR/N6+320FqE6GnhFkicuxr0kSUumaE35K0nSSBar29Vx\nwG1VtbnpanExrXnNJWlv1/N1c5L/1TFbU+f2/qXOYD/NFLeXL3c+5ivJY/vU9feSHLHc+ZOkSbZY\ns10dzu5TCW4Bjl+ke0nSilFVU33Cf43WGhFaZM0aKT3Xl5EkLa7FevPhgDpJkiRJu1msNx930lrA\nqm0drbcfOyWxgSJJkqQlUVVzmYJ9Wcz1+XgllKnbYjU+rgXWN9Myfo3WQkyv6I606u7vL9Lt1bbj\nnW9n1W+dudzZmGgrqY5nDn74cmdhJG9961t561vfutzZmHjW8+Kzjhefdbz4VmIdz22N0eX1u0PG\n+++LmovFsyiNj6qaTnIarbnfp4ALqmohFnWSJEmSJtY+y52BRbZYbz6oqsuAyxYrfUmSJGnSLNrD\n+ZiY9PLt9fLUpy93Fiaedbz4NmzYsNxZ2CtYz4vPOl581vHis44X1/7LnYFFtmgrnM9646Qc8yEt\nrZU65kOSpPlIsiIGZyepDw4Z97U44HzR7JhpZXPV1PQeYW2d54ZNr9f1ve41TH6GuUe/e45y7Zr9\nHtq5P719qud1o9bJXK5brPTa13eWs61d3tnu24o8tfvP/bbNK18979OkveaA+3vmc7e4D64BYP+D\n7gFgprl2qDLdse/OsOcf/3EALrviP7UCNrV+fOzskwD4+ddeujPuER+8rSO1o/reR5IkLb8V8XA+\nD5NePkmSJGnFmPQB54u1yKAkSZKkOVo95NZLkhOT3Jzk1iSn94nz3ub89UmObcJ+OMl1Hdt3k7yx\nOffWJFs6zp043/JJkiRJGgOjvvlIMgW8D3gOrQW/P59kY+dyF0lOAo6qqvVJjgc+AJxQVf8KtBsi\nq5rrP9FcVsB5VXXeiFnbzbI2PlbvMzOwr/tCGWaMRrfV+8wAg/viL7SBYz3aYxZWz+xxqp3Xbe3r\n23Gb8rbPA2x7sDVuoD2Wolf5+o23mEtdrJqa7lueudTtVI/ydl83sN4ebNdJ++Yd1/aon+70u/O6\ns64f2jX+on2PqbV90pnuuOeDrXFh7XK1f05v3zW3xR71075+665kfpbLAbjsy82Yj2Z02n86sJnd\n+kPf2xl33Qfv6MiVYz4kSRpn8+h2dRxwW1VtBkhyMXAy0LnW3ouACwGq6uoka5McXFUdTxk8B7i9\nqjofIBZsYLvdrrQohhk0L0mSpN3tP+TWw+FAZ4NhSxM2W5wjuuK8HLioK+wNTTetC5KsHbIoPdn4\nkCRJksZEvzEeXwb+tGPrYdj1M7rfYuy8Lska4IXAX3Sc/wDwOOAY4C7gXUPepyf/PC1JkiSNiX7d\nrn6y2dou3DPKncC6juN1tN5sDIpzRBPW9nzgC1X1jXZAVX29vZ/kQ8Cn+mZ+CL75kCRJksbEPGa7\nuhZYn+TI5g3Gy4CNXXE2Aq8ESHICcE/XeI9XAB/pvCDJoR2HLwFuGKFYO62INx/dCwDO+/oeg7d3\nDvJtjgcNhp7P4oLDLJTYawG7dl4HDprvMyi9sywLsYDgIKN+Rv0Goc90DNbuNfh89oSbN4nTzRvG\nHp95d/qDBuGzT488TO9+fWee+znwYfcCcP+2/XdPv989uhzF7a2dbzYB3/xS6+e72zH69AaVJElj\nbdQB51U1neQ04HJgCrigqm5Kcmpz/vyqujTJSUluA+4DXtW+PsnDaQ02f21X0u9Icgyt7llfAU4d\nMYvACml8SJIkSXuD+TycV9VlwGVdYed3HZ/W59r7gIN6hL9yHlnag40PSZIkaUxM+grnNj4kSZKk\nMTHpHadXbuOjz/iGzn78c1kUr19//15jGEYZN9HrmlHSGWqBvq5FBntGGWJMSzvOvvtvm/WaYeqp\ne4xFr3EOe4zLmasBizF252kuYzT63gc6xnz0qe/d8tIs4EhrAcdtq9fsmV5jUF0cyebWzs65KFqL\nkJ61cz6KSf+7iSRJk2nSf4Ov3MaHxtpiD2yXJEmaRJP+cD7p5ZMkSZJWjH2GfTpfoX/ntfEhSZIk\njYnVNj4kSZIkLYV9Rhz2ulIsa+NjavXMUIPChxmE3D2YeeD1PQYltwcfDzOgu53eoHENoyyMOEzc\noQbRN+XbsXrPfA6T94XQ6/MYZYB3Z3m7Fxkcqm7b1/T4M8Iw+dl5/1EHwO+Rn9aPh/EAAA+sethI\n6c/QxN85JcZRHf8FfmxX3Acmft4MSZImx9BvPlaoCS+eJEmStHLss+9y52Bx2fiQJEmSxsWEP51P\nePEkSZKkFWTCn84npnjzWjBunuYyrqPXdfNZbBDmtphiv3sPykP3WIte6XTb9uC+rNnvoTmnt9g6\n623XooBr5pdok8zqnd/BAd+H1dW6d3PRFHOoi45k7+XA1s4B7ZCuMR8/vCvu/Txs+HtIkqTlNTFP\n571NePG0XPo1PCRJkjSAs11JkiRJWhIT/nQ+4cWTJEmSVhBnu5IkSZK0JCb86XziitdrMPO2PRYZ\nTI/r5j74eNSB5mOhXReLvNhgp0GTAuyxmF8Tt7OOZ6b7DM7uTK97Acmd34c9P6vVfQa+dw5K3/bg\nvnveA3p+h+ZiX7bNPb2OItzD2tbOfu2QJwJw/NOawyN3xXXAuSRJK8gKfrwcxqrlzoAkSZKkxtSQ\nWw9JTkxyc5Jbk5zeJ857m/PXJzm2I3xtko8muSnJjUlOaMIfleTKJLckuSLJ2vkUz8aHJEmSNC5W\nD7l1STIFvA84ETgaeEWSJ3bFOQk4qqrWA68DPtBx+j3ApVX1ROBJwE1N+JuBK6vqCcBVzfHIbHxI\nkiRJ42LExgdwHHBbVW2uqu3AxcDJXXFeBFwIUFVXA2uTHJzkEcDTq+rDzbnpqvpu9zXNzxfPt3gr\nU3ff/hF0LqzX7v/fHpfQc+G+BbjnYowTaY9R2DZGS2sMWmRwWTSL+83ZCItWtse07PxcOj/z+Syw\n2GvMR/vF54E/0Pr50j3jPrBj/10H/rlBkqTxNvqj4uHAHR3HW4Djh4hzBDADfCPJHwFPBr4A/Jeq\nuh84uKq2NvG3AgePnEN8FNEiGauGhyRJ0kqx75Dbnob9S2v3LDdFq8nzFOD9VfUU4D56dK+qqprD\nfXpauW8+JEmSpEnT5+l80zda2wB3Aus6jtfRerMxKM4RTViALVX1+Sb8Y0B7wPrWJIdU1d1JDgW+\nPnsh+hv5zUeSX0zyL0lmkjyl69wZzSj6m5M8bz4ZlCRJkvYafWa32nAIvPXHdm09XAusT3JkkjXA\ny4CNXXE2Aq8EaGazuqeqtlbV3cAdSZ7QxHs28C8d15zS7J8CfHI+xZvPm48bgJcA53cGJjmaVmGP\nptWv7G+SPKGqdszjXpIkSdLkG/HpvKqmk5wGXE6ryXJBVd2U5NTm/PlVdWmSk5LcRqtr1as6kngD\n8GdNw+X2jnPnApckeQ2wmV0jTEcycuOjqm4GSPZYHO1k4CPNKPvNTeGOAz436r1WslVLsIhfz8Hx\n87CiF0/s0q7/Hc3ikb0WoWzrtfjhTvNYVLDzO7CjuccUcxh43mOhxJ0Dzg9oAn649eO+U5uXmZ/e\ndfn93+9YZPAHhr+tJElaBvN4DKuqy4DLusLO7zo+rc+11wM/2SP828BzRs/V7hZjwPlh7N6/bAut\nNyCSJEmSBpnHIoMrwcC2VZIrgUN6nHpLVX1qDvfpOSp++znnsqP5q32e+nTytGfMIUlJkiRpT5s2\nbWLTpk3LnY3RTE4HlJ4GFq+qnjtCmv1G0e9hnzPezI4H1oxwC0mSJKm3DRs2sGHDhp3HZ5111vJl\nZq72W+4MLK6Falt1dojfCFyU5Dxa3a3WA9cs0H36avflH9hvv23xh2GsOO0F8aBrUTxgavUiVVjn\nZzXK2JgRFgDsNLVq7gv+zXt8Tb8xJ8MUv+Nj+SYHtXbaiww+sfXjz/d9eWun433lQ50NfMd8SJI0\n3lZwl6phjNz4SPIS4L3AQcCnk1xXVc+vqhuTXALcSOuR6vXNgiSSJEmSBtmbu10NUlWfAD7R59zZ\nwNmjpi1JkiTtlWx8SJIkSVoSE97tajGm2l08PfrLr5qaHnktjR3feviscXqlPZ97Drqune5Crw0y\nzLodC71WyLYH92Vq9czAdTX2X3vvgt5zGA98c+3skTpNTw09tmR6eorpHnE7x9P0c+jRX5n9Bqt3\nbd/i0XyLR7fW+TiA1rQOR8BF/BIX8Us8+6mfZtUP3ceqH7pvqLxLkqQxsXrIbYVaWVmf5wDjbqse\nPfuD2UIvuDdMegt9z2EaM8M8IM/F/gfcP2ucB+45cEHvOYz9D7pnye85TMPurhsft6D3/NutGxY0\nPUmStERW1tP5nE148SRJkqQVZN/lzsDisvEhSZIkjYsJfzqf8OJJkiRJK8iEP51PePEa7bEi0xkc\nbwH0G1/RHlPRvYDfUui5gGC7Tubwaq9zLMoog+LbixUOXAiyfa7HQPV5DYofMPC9W89B8t3FXeaF\nKrfymNZOs9YgR7R+XPW15wCjfT6SJGkMTPhsV3tH40OSJElaCSb86XzCiydJkiStIBP+dD7hxZMk\nSZJWELtdSZIkSVoS+y13BhbXymp8DLHIYOdg4YEDmxfwnlKnQd+79kDwmVH+rNFxybfaI83XVuvn\nEdktfUmStEKtrKfzOVu13BmQJEmS1JgacushyYlJbk5ya5LT+8R5b3P++iTHdp2bSnJdkk91hL01\nyZYm/LokJ86neBPetpIkSZJWkBGfzpNMAe8DngPcCXw+ycaquqkjzknAUVW1PsnxwAeAEzqS+S/A\njcCBHWEFnFdV542Ws9355kOSJEkaF6uH3PZ0HHBbVW2uqu3AxcDJXXFeBFwIUFVXA2uTHAyQ5Ajg\nJOBDQPfieAu2WN7KaHxMTzn2Ylz1+Wwe+P7DliEzi2A6Iy9OOb19ataFEWeYGn78xz67tq08hq08\nhv0Puof9D7oHjqS1SZKklW30bleHA3d0HG9pwoaN8wfAbwE7eqT9hqab1gVJ1s6hNHuw25UWxf4H\n3L/cWZAkSVp5+sx2tek62PTFgVfWkHfY461Gkp8Dvl5V1yXZ0HX+A8DvNPtvA94FvGbIe+3Bxock\nSZI0Lvp0iNjwE62t7az/vUeUO4F1HcfraL3ZGBTniCbs54EXNWNC9gN+IMkfV9Urq+rr7chJPgR8\ninlYGd2uJEmSpL3B6GM+rgXWJzkyyRrgZcDGrjgbgVcCJDkBuKeq7q6qt1TVuqp6HPBy4P9UVTve\noR3XvwS4Yb7FkyRJkjQORnw6r6rpJKcBl9N6f3JBVd2U5NTm/PlVdWmSk5LcBtwHvKpfch3770hy\nTBP2FeDU0XLYslc3PlbvMzN7pAHaC7rtmFmCamwP6nYRuYkw0iKDHe7hkQA8rBlb89DBa+adJ0mS\nNAbm8VhZVZcBl3WFnd91fNosaXwG+EzH8StHz9Ge9urGhyRJkjRWJnyCVxsfkiRJ0riY8KfzCS+e\nJEmStILsu9wZWFw2PiRJkqRxMeFP5xNevL3MBK0Cv6SD+RfIzPTc8zqzo/nMhrm0I849989rcVFJ\nkjSuVs6jz0gmvHiSJEnSCjLhT+cTXjxJkiRp5ajJ6cjSk40PSZIkaUysoB7nI5nw4kmSJEkrx6Q3\nPlYtdwbmZPX8ViTvtu3uH5g1zqp9H5o9To9Vx3fMrO45WHrNfkOmNz01cAD59Pbh38lte3D2OduG\nSe8RB31n6Hs+8P2H7ZZ2r/T3X3vvrOl01kW7TlfvM7Pn6vSrZ4b6fjxwz4GzxpkZYuD+qkff13Fv\n5tyMn2aK6Y5VhB591J27zvWpL9bu+p5te3AN2x5cw/337s/99+4/t5tLkqSx9dC+a4baVqqV1bZa\n4Nmc1hzyvVnj7Hho9gf3uczINExDYJh77vHwPcAwDZ5h0vvuNx859D33P+D+WeMM0xAYpi7m0igd\npsEzNUR6O7718KHvOYxv3Xb47JHuWVn/u0qSpLmbmZrsQR8+zUiSJEljYgYbH5IkSZKWwLSND0mS\nJElLYWbCH89HHnCe5J1JbkpyfZKPJ3lEx7kzktya5OYkz1uYrEqSJEmTbYapobaVaj6zXV0B/EhV\nPRm4BTgDIMnRwMuAo4ETgfcnWVmzakmSJEnLYNIbHyO/16mqKzsOrwZ+vtk/GfhIVW0HNie5DTgO\n+Fx3GsNMabqbzlmI5nLtzusm+zXWOHng+w8basartrnMGLbTELNStWfx2jazsNM0z1f7ler0nL7H\nu3bbs6b1muZZkiStXA+xcqfRHcZCvZF4NXBps38YsKXj3BZgiHlENUnm0vCQJElSywyrh9p6SXJi\nM+zh1iSn94nz3ub89UmObcL2S3J1ki8muTHJOR3xH5XkyiS3JLkiydr5lG9g46O50Q09thd2xDkT\n2FZVFw1IquaTSUmSJGlvMGq3qyRTwPtoDXs4GnhFkid2xTkJOKqq1gOvAz4AUFUPAj9TVccATwJ+\nJsnTmsveDFxZVU8ArmqORzawr0tVPXfQ+SS/CpwEPLsj+E5gXcfxEU3YHqbfcQ47ZtJK66lPJ097\nxuw5liRJkgbYtGkTmzZtWu5sjGQe4zmOA26rqs0ASS6mNRzipo44LwIuBKiqq5OsTXJwVW2tqna3\nlTXAFPCdjmue2exfCGxiHg2QkQdBJDkR+C3gmU1rqW0jcFGS82h1t1oPXNPz5qefwfT2lTtgRpIk\nSeNnw4YNbNiwYefxWWedtXyZmaN5rPNxOHBHx/EW4Pgh4hwBbG3enHwB+CHgA1V1YxPn4Kra2uxv\nBQ4eNYMwvxHYf0irZXRlEoB/qqrXV9WNSS4BbgSmgddXld2uJEmSpFn0G8/xz5vu5Z83fX/QpcM+\nb6fXdVU1AxzTLJ9xeZINVbVpt4hVlWRez/Xzme1q/YBzZwNnj5q2JEmStDfq1+3qyRvW8uQNu8Z6\nX3DW1u4o3UMf1rH7JFC94uwxPKKqvpvk08CP0+pitTXJIVV1d5JDga8PW5ZeXH9DkiRJGhPbWDPU\n1sO1wPokRyZZQ2vdvY1dcTYCrwRIcgJwT1VtTXJQexarJPsDzwW+2HHNKc3+KcAn51M+F76QJEmS\nxsSoYz6qajrJacDltAaMX1BVNyU5tTl/flVdmuSkZh2++4BXNZcfClzYLAy+CviTqrqqOXcucEmS\n1wCbgZeOWDRgb2t8rG66qE13d3UbUntBuH0XJjvLob0o3Y4JX8BmObQXzRxmEoU5zWTRXkzxgF1B\nLi4oSdJk6jfmYxhVdRlwWVfY+V3Hp/W47gbgKX3S/DbwnJEz1WXvanxIkiRJY2weU+2uCDY+JEmS\npDFh40OSJEnSkpjHOh8rgo0PSZIkaUxsW8mDi4dg40OSJEkaE3a7kiRJkrQk7HYlSZIkaUnMZ6rd\nlWCySydJkiStIHa7kiRJkrQkbHxIkiRJWhI2PiRJkiQtiYecancMrJ5p/ZxeoJbgqKVu8rF6n9bP\n6e3j0TLdMbP0H+Ns95xZqM9qXHUWf3rulz/EGgBmppuEVtfOc/3qdtW+D839RpIkaUXxzYckSZKk\nJWHjQ5IkSdKScJ0PSZIkSUvCdT4kSZIkLQm7XUmSJElaEpPe+Fi13BmQJEmS1PIQa4baeklyYpKb\nk9ya5PQ+cd7bnL8+ybEd4R9OsjXJDV3x35pkS5Lrmu3E+ZTPxockSZI0JmZYPdTWLckU8D7gROBo\n4BVJnthxwO+1AAAShUlEQVQV5yTgqKpaD7wO+EDH6T9qru1WwHlVdWyz/fV8ymfjQ5IkSRoTM0wN\ntfVwHHBbVW2uqu3AxcDJXXFeBFwIUFVXA2uTHNIcfxb4Tp9sZSHKBjY+JEmSpLExj8bH4cAdHcdb\nmrC5xunlDU03rQuSrJ1LebrZ+JAkSZLGxDRTQ2091JC36H6LMdt1HwAeBxwD3AW8a8j79ORsV5Ik\nSdKY6LfOx12bbuGuTbcOuvROYF3H8TpabzYGxTmiCeurqr7e3k/yIeBTg+LPxsaHJEmSNCb6TbX7\nmA1P5DEbdo0fv+6sy7qjXAusT3Ik8DXgZcAruuJsBE4DLk5yAnBPVW0dlJ8kh1bVXc3hS4AbBsWf\njY0PSZIkaUxs6zON7myqajrJacDlwBRwQVXdlOTU5vz5VXVpkpOS3AbcB7yqfX2SjwDPBB6d5A7g\nf1TVHwHvSHIMre5ZXwFOnUfxbHxIkiRJ46LPeI6hVNVlwGVdYed3HZ/W59rutyTt8FeOnKEebHxI\nkiRJY6LfmI9JMdmlkyRJklaQfmM+JoWND0mSJGlM2PgYV6tnFjS56enF+aB3zDRVvE8rv6umpvc8\nNwamhqnPzjqaJf62B/dlzX4P7RY2Mz0+5V1u29gXgJkB37vO74okSdo7zGfMx0rg06AWRXfDQ5Ik\nSbNzzIckSZKkJTHqVLsrhY0PSZIkaUxMererVaNemORtSa5P8sUkVyVZ13HujCS3Jrk5yfMWJquS\nJEnSZJth9VDbSjVy4wP4vap6clUdA3wS+G2AJEfTWs79aOBE4P1J5nMfSZIkaa8ww9RQ20o1cqOg\nqu7tODwA+GazfzLwkaraXlWbgduA40bOoSRJkrSXmPTGx7ze2SR5O/ArwAPsamAcBnyuI9oW4PD5\n3EeSJEnaG6zkhsUwBjY+klwJHNLj1Fuq6lNVdSZwZpI3A+8GXtUnqZpfNiVJkqTJ91CzFtikGtj4\nqKrnDpnORcClzf6dwLqOc0c0YXuYfsc57JgJAHnq08nTnjHk7SRJkqTeNm3axKZNm5Y7GyPZq998\nDJJkfVXd2hyeDFzX7G8ELkpyHq3uVuuBa3re/PQzmN4+2RUsSZKkpbVhwwY2bNiw8/iss85avszM\nkY2P/s5J8sPADHA78OsAVXVjkkuAG4Fp4PVVZbcrSZIkaRaTvs7HyI2PqvqFAefOBs4eNW1JkiRp\nb7SS1/AYxmSXTpIkSVpB7HYlSZIkaUlMeuPDlcclSZKkMfHQtjVDbb0kOTHJzUluTXJ6nzjvbc5f\nn+TY2a5N8qgkVya5JckVSdbOp3w2PiRJkqQxMTO9eqitW5Ip4H3AicDRwCuSPLErzknAUVW1Hngd\n8IEhrn0zcGVVPQG4qjkemY0PSZIkaUzMTE8NtfVwHHBbVW2uqu3AxbSWw+j0IuBCgKq6Glib5JBZ\nrt15TfPzxfMpn2M+JEmSpDHRp2ExjMOBOzqOtwDHDxHncOCwAdceXFVbm/2twMGjZhBsfEiSJElj\no98C3PUPf0f942cHXTrsunoZMs4e6VVVJZnX+n02PiRJkqQxsWOmz+P5Cc9qbW3vOqc7xp3Auo7j\ndbTeYAyKc0QTZ58e4Xc2+1uTHFJVdyc5FPj6EMXoyzEfkiRJ0riYnhpu29O1wPokRyZZA7wM2NgV\nZyPwSoAkJwD3NF2qBl27ETil2T8F+OR8iuebD0mSJGlcPDja43lVTSc5DbgcmAIuqKqbkpzanD+/\nqi5NclKS24D7gFcNurZJ+lzgkiSvATYDLx29cDY+JEmSpPExPfqlVXUZcFlX2Pldx6cNe20T/m3g\nOaPnanc2PiRJkqRxMY/Gx0pg40OSJEkaFzY+JEmSJC2J7cudgcVl40OSJEkaFzPLnYHFZeNDkiRJ\nGhd2u5oAoy9Trwk2s2P370V7RdGp1f3/5NBz4Z8maKbP96zzmlVTrX9R7t+2PwAPPbCmdWK/bXvE\nkSRJe6EHlzsDi2vvaHxIkiRJK8GE/w3SxockSZI0Lmx8SJIkSVoSNj4kSZIkLQmn2pUkSZK0JJxq\nV5IkSdKSsNuVJEmSpCXhVLuSJEmSloRvPiRJkiQtCRsfkiRJkpaEjQ9JkiRJS2LCp9pdtdwZkCRJ\nktSYGXKbgySPSnJlkluSXJFkbZ94Jya5OcmtSU7vCP/FJP+SZCbJUzrCj0zyQJLrmu39s+XFxock\nSZI0Lh4ccpubNwNXVtUTgKua490kmQLeB5wIHA28IskTm9M3AC8B/q5H2rdV1bHN9vrZMmLjQ5Ik\nSRoX00Nuc/Mi4MJm/0LgxT3iHEerIbG5qrYDFwMnA1TVzVV1y5zv2oOND0mSJGlcbB9ym5uDq2pr\ns78VOLhHnMOBOzqOtzRhs3lc0+VqU5Kfni2yA84lSZKkcTHH8RxtSa4EDulx6szOg6qqJNUjXq+w\n2XwNWFdV32nGgnwyyY9U1b39LrDxIUmSJI2Lfl2qtmyCOzf1vayqntvvXJKtSQ6pqruTHAp8vUe0\nO4F1HcfraL39GHTPbcC2Zv+fk9wOrAf+ud81Nj4kSZKkcdGv8XHIhtbWds1Zc0l1I3AK8I7m5yd7\nxLkWWJ/kSFpvNF4GvKJHvOzcSQ4CvlNVM0keT6vh8W+DMuKYD0mSJGlcLM6Yj3OB5ya5BXhWc0yS\nw5J8GqCqpoHTgMuBG4E/r6qbmngvSXIHcALw6SSXNek+E7g+yXXAXwCnVtU9gzIy7zcfSX4TeCdw\nUFV9uwk7A3g1rV5rb6yqK+Z7H0mSJGniPbTwSTbP6M/pEf414AUdx5cBl/WI9wngEz3CPwZ8bC55\nmVfjI8k64LnAv3eEHU3rNc3RtEbI/02SJ1TVjvncS5IkSZp4c59Gd0WZb7er84A3dYWdDHykqrZX\n1WbgNlrzBkuSJEkaZHG6XY2NkRsfSU4GtlTVl7pOHcbuI+OHnSNYkiRJ2rvNDLmtUAO7Xc0yX/AZ\nwPM6ow9Ique8wdPvOIcdM63L8tSnk6c9Y2BmJUmSpNls2rSJTZs2LXc2RjPh3a4GNj76zRec5EeB\nx9Ea3Q5wBPCFJMez5xzBRzRhe9789DOY3j41QrYlSZKk3jZs2MCGDRt2Hp911pympV1ee3Pjo5+q\n+jIdy7In+Qrw41X17SQbgYuSnEeru9V64JqFyKwkSZI00VbweI5hLNQigzu7VVXVjUkuoTU/8DTw\n+qoaZbl2SZIkae+yCFPtjpMFaXxU1eO7js8Gzl6ItCVJkqS9ht2uJEmSJC0Ju11JkiRJWhIreBrd\nYdj4kCRJksaF3a4kSZIkLQkbH5IkSZKWhGM+JEmSJC0Jp9qVJEmStCTsdiVJkiRpSdjtSpIkSdKS\ncKpdSZIkSUtiwrtdrVruDEiSJElqTA+5zUGSRyW5MsktSa5IsrZPvBOT3Jzk1iSnd4S/Lcn1Sb6Y\n5Kok6zrOndHEvznJ82bLi40PSZIkaVxsH3KbmzcDV1bVE4CrmuPdJJkC3gecCBwNvCLJE5vTv1dV\nT66qY4BPAr/dXHM08LIm/onA+5MMbF/Y+JAkSZLGxSK8+QBeBFzY7F8IvLhHnOOA26pqc1VtBy4G\nTgaoqns74h0AfLPZPxn4SFVtr6rNwG1NOn055kOSJEmabAdX1dZmfytwcI84hwN3dBxvAY5vHyR5\nO/ArwAPsamAcBnyu65rDB2XExockSZI09jY1W29JrgQO6XHqzM6Dqqok1SNer7DO684EzkzyZuDd\nwKv6RR2UzsQ1Pmamp3buT2+fGhBzjK3umGOtXZ4mbPU+Czv/2qD0Vk213unNTC/h12T1zODjuV4/\nwHRTt93fk87vEKvb//+k6z6d/19lt3R2zMxeX/d+58Ddjtt1LUmS1NuGZms7a7ezVfXcflcm2Zrk\nkKq6O8mhwNd7RLsTWNdxvI7Wm4xuFwGX9rnmiCasL8d8SJIkSWNjUUacbwROafZPoTVovNu1wPok\nRyZZQ2sg+UaAJOs74p0MXNeR7suTrEnyOGA9cM2gjEzcmw9JkiRp5VqU3hDnApckeQ2wGXgpQJLD\ngA9W1QuqajrJacDlwBRwQVXd1Fx/TpIfprUE4u3ArwNU1Y1JLgFubDL++qrau7pdSZIkSSvX3OfR\nnU1VfRt4To/wrwEv6Di+DLisR7xfGJD22cDZw+bFxockSZI0Nh5Y7gwsKhsfkiRJ0thY+Dcf48TG\nhyRJkjQ2JnsGTBsfkiRJ0tjwzYckSZKkJeGbD0mSJElLwjcfkiRJkpaEs11JkiRJWhJ2u5IkSZK0\nJOx2JUmSJGlJ+OZDkiRJ0pLwzYckSZKkJeGbD0mSJElLwjcfkiRJkpaEU+1KkiRJWhK++ZAkSZK0\nJCZ7zMeqUS9M8tYkW5Jc12zP7zh3RpJbk9yc5HkLk1VJkiRp0m0fcluZRm58AAWcV1XHNttlAEmO\nBl4GHA2cCLw/yXzuo/n4x88sdw4m39WbljsHE2/Tpk3LnYW9gvW8+KzjxWcdLz7reLFND7mtTPNt\nFKRH2MnAR6pqe1VtBm4DjpvnfTSqf/q75c7B5Ltm03LnYOL5i25pWM+LzzpefNbx4rOOF5tvPgZ5\nQ5Lrk1yQZG0TdhiwpSPOFuDwed5HkiRJ2gtM9puPgQPOk1wJHNLj1JnAB4DfaY7fBrwLeE2fpGrU\nDEqSJEl7j8meajdV828XJDkS+FRV/ViSNwNU1bnNub8Gfruqru66xgaJJEmSlkRV9RouMFbm+ny8\nEsrUbeSpdpMcWlV3NYcvAW5o9jcCFyU5j1Z3q/XANd3Xr8TKkiRJkhbL3vB8PJ91Pt6R5BhaXaq+\nApwKUFU3JrkEuJFWh7TX10K8XpEkSZK0oi1ItytJkiRJms2irr+RZHOSLzWLEF7ThD0qyZVJbkly\nRccsWS5OOII+dfzOJDc1M5F9PMkjOuJbx3PUq447zv1mkh1JHtURZh2PoF89J3lD833+cpJ3dIRb\nz3PU59+L45Jc04R9PslPdsS3jucoydokH22+szcmOd7fewurRx2f4O+9hdernjvO+btPo6uqRdto\ndcd6VFfY7wFvavZPB85t9o8GvgjsAxxJa32QVYuZv0nY+tTxc9t1B5xrHS98HTfh64C/7jxvHS9s\nPQM/A1wJ7NMc/6D1vOB1vAn42Wb/+cDfWsfzquMLgVc3+6uBR/h7b0nq2N97S1DPzb6/+9zmtS3F\nyuPdA2de1HyhaX6+uNl3ccLR7VbHVXVlVe1oDq8Gjmj2rePR9RoAdh7wpq4w63h+uuv514Fzqmo7\nQFV9owm3nkfXXcd30Xp4A1gL3NnsW8dz1Py1/elV9WGAqpququ/i770F06+O/b23sAZ8l8HffZqn\nxW58FPA3Sa5N8tom7OCq2trsbwUObvZdnHA0veq406uBS5t963g0e9RxkpOBLVX1pa641vHoen2X\n1wPPSPK5JJuS/EQTbj2Pplcdvxl4V5KvAu8EzmjCreO5exzwjSR/lOSfk3wwycPx995C6lXHD+uK\n4++9+etZz/7u00KYz2xXw3haVd2V5AeBK5Pc3HmyqiqD5zN2NPzs9qjjqvosQJIzgW1VddGA663j\n2fX6Hp8BdPZpHTQ1nnU8nF71vBp4ZFWd0IxFuAR4fJ/rrefZ9arj3wbeWFWfSPKLwIdpdWHpxToe\nbDXwFOC0qvp8knfTatzt5O+9eetXx/8D/L23gHrV81nA0/F3n+ZpUd98VLMOSNNV4hO0XsFtTXII\ntNYKAb7eRL+TVj/CtiPY9fpfffSpY5L8KnAS8J87olvHI+hRx8+k9Veh65N8hVY9fiHJwVjHI+vz\nXd4CfLwJ/zywI8lBWM8j6VPHx1XVJ5ooH2VXVwnreO620Pqr8Oeb44/SeoC72997C6ZfHft7b2H1\nqudjaY3n8Hef5mXRGh/N67kDm/2H02op30BrEcJTmminAJ9s9jcCL0+yJsnj6LM4oXbpV8dJTgR+\nCzi5qh7suMQ6nqM+dXxNVR1cVY+rqsfR+kf6KU23Cut4BAP+vfgk8Kwm/AnAmqr6JtbznPWp4y8D\ntyV5ZhPtWcAtzb51PEdVdTdwR/NdBXgO8C/Ap/D33oLoV8f+3ltYfer5C1V1iL/7NF+L2e3qYOAT\nSdr3+bOquiLJtcAlSV4DbAZeCi5OOKJ+dXwrsIZWtwqAf6qq11vHI+lZx11xdtahdTyyft/lfYAP\nJ7kB2Aa8EqznEfWq48uTfAv4n0n2BR4AXgfW8Ty8AfizJGuA24FXAVP4e28hddfxq4HP4++9hdbr\nu9zJ330aiYsMSpIkSVoSSzHVriRJkiTZ+JAkSZK0NGx8SJIkSVoSNj4kSZIkLQkbH5IkSZKWhI0P\nSZIkSUvCxockSZKkJWHjQ5IkSdKS+L8Q49KHBgGMlAAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig, ax = plt.subplots(1,1,figsize=(15,4))\n", "t=0\n", "mesh=ax.pcolormesh(xx,zz,Ws[t,:,:])\n", "ax.set_title('W: t={} : {}'.format(t, run))\n", "fig.colorbar(mesh,ax=ax)\n", "ax.set_xlim([500,650])\n", "ax.set_ylim([-50,0])" ] }, { "cell_type": "code", "execution_count": 20, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 20, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAA3AAAAEKCAYAAABTzLdZAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XmcHHWd//H3Zyb3TMhBYggQCEIIl8uhcggsAx6LILAq\n7IrCKoruKh6ouIq4a0Rd1hUV/SneJyt4LSIIKJFlAA9ANOEIsAgIJCEHOSaZSSbHZD6/Pz5V6ZpO\nz2T6mOrpzOv5eNSjj+qu+nZ1JTPv+Xzr+zV3FwAAAABg+GuqdwMAAAAAAINDgAMAAACABkGAAwAA\nAIAGQYADAAAAgAZBgAMAAACABkGAAwAAAIAGQYADgAqY2Wwz6zWzXe7/UTO7xczOr3c78mRm3zOz\nT9Z4m1WdI2b2tJmdktz/qJl9M7PutWa22Mw6zexwM5trZgvNbL2ZvbtWn6HWzOyrZvaxercDABrZ\nLveLBwA0gmoCg5lNNbOfm1lX8kv+uVW0Y56ZXZN9zt1Pc/dr+nvPUMiGlTrxZMmVmZ1rZj/sZ/X2\n9rj7f7j72zPrrpT0Lnef6O4PSPpXSbe7+27u/mUzO9nM7jCzDjP7aw3aOcHMrjaz55Nt3lnJdtz9\nne7+qWrbAwAj2ah6NwAAULavSNok6QWSjpR0s5k94O6P1LdZVXFJVu9G1MHpkm4u5w1mZpL2kZT9\nvveV9PvM4y5J35I0QdJHq2yjJH1D8UffgyStkXREuRswsyZ3761BW4q3a5Lk7rkHcACoBypwACDJ\nzD5sZkuSLmiPZbquHW1m95vZOjNbbmafK3rreWb2TFKZGNQvymb2DklvlPSvSRe4X5TRzhZJr5P0\nb+6+0d1/J+kXksru8mhmp0q6VNI/Ju1YkDzfbmZvS+6/xcx+Z2afN7O1ZvaEmb3MzC4ws2fNbIWZ\n/VNmm2PN7MrkmCxPusyNS9ZNM7NfJttZbWZ3WbhGEUhuStpxSfL6Y83s98nrF5rZSZn9tJvZFWZ2\nb/Ld3GBmUwbxmU/IbPPZbNslTU3at97M7jGzF2bed5CZzU/a/ZiZnZNZN97MPpdUETvM7G4zG1ti\n3683s7+a2SHJ4yZJr5D0q+Tx+clxW1V8LqWVUjMbI6lTUrOkB5Lv43ZJbZK+nLT9AHf/o7v/UFLJ\n6ttAn6fUayWdIekd7r7aw4KBj/T2KvNXLbrkdkk62TKVZzN71MxOz7x+VPLv6Ijk8c6+/0+Z2e8k\nbZC0387aAwC7CgIcgBHPzOZKukjSS9x9N0mvkvR0svqLkr7g7pMkvVDST4refrykAyW9XNK/J9tK\ng8LaUvtz929I+qGkzyRd4M5K3pOGm1LLjcnbD5TU4+5PZDb5gKRD+/ls+yTv37tEO34l6T8k/Shp\nx5HpKvXtTnh0so+pkq5LjsFRkvaXdJ4iOExIXvufkg6QdHhyu5ekf0/WfVDSYknTFNXDS5MwcL6k\nZyW9JmnHlWa2l6RfSrrc3adIukTS/5jZ7pl2nS/pAkkzJfVI+lKpY5A5FvtKukXxnU5TVJEeSFdL\neoOkeZKmSHpC0qeT97VImi/pvyVNT153tZkdnLz3SkUl9LjkGH2o6PiZmV2QHJuXZyqlR0t6yt3X\nJKHuaklvkrSnpN0lZb8zlyR33+Lurclzf+PuB7j7yyXdLemipAtl9twodRx29nmKHS3pGUmXJwHr\nQTN73UD7yDhX0ieTNv9Wfc+ta5P1qb+TtNLdFw7y+z9P0oWSWhXnDwCMCAQ4AJC2SRor6VAzG+3u\nz7r7U8m6LZLmmNm0pOJ1b9F7P+Hum939QUUYOEKS3P23yS+eA+nTZdDdX+PuU/pZzkxe1ippfdF2\nOiVNLLWD5LNMcfclA7RhZ10X/+ru30+6qP1EETAud/et7j5fcYwOMDOT9HZJH3D3DnfvknSFIiAo\ned1MSbPdfVtSPezPeZJuSUKm3P03ku5XdDmUIgT8wN0fcfeNkv5N0j8kbejPGyXNd/cfJ/tfk1w/\nlm7vene/3923KQJ22k3wNZlj0OvuCyVdL+mcpIp2gaT3ufuyZP097r4ls9/3KwLISZnzSurbffJs\nSTcl582W5PNkuxsOpnvpYLug9vt5+nn93pIOk9Sh+P7eLen7SWVuIC7pBnf/gyS5++aidl4n6cy0\nQqv4fq5L7g/m+/+euz+afIaewXxwANgVEOAAjHhJxeJiRfVlhZldZ2Yzk9VvU1S9HjWz+7JdvhLL\nM/c3SmoZ4uZ2Sdqt6LlJihA3VFZk7ndLkrs/X/Rcq6KaM0HSn9LKoaRbFdUuSfqsorJ1m5k9aWYf\nHmCf+yoC0trMto6XtEfmNYsz95+VNDqzr1L2lvTUAOuLP2da6dpX0jFFbXmjpBmKStk4SU8OsN0P\nSvqKuz9X9PyrFRVBKYLR9pCdhNLVA2yzlMFeA9bv5zGzWRaD43SaWfqHgm5JWyV9yt173P0uSXco\nKtU7s7i/Fcm/u0cVIW6CopvmtZk2lvP9A8CIwSAmACDJ3a+TdJ2ZTZT0dUmfkfRPyS+Zb5TiGiZJ\nPzOzqbXYZfETZnarpBP6ef1d7n66pMcljUquc0q7yh0u6eEK21HLQSVWKX7ZP8TdlxWvTCpyl0i6\nxMwOlfS/Znafu9+hHY/Hs5Kucfd3DLC/fYrub03a0J/Fiu6A5XpW0p3uvkNgSSpwmxTdRR/s5/2v\nkvRrM1vu7tcn79tD0szMtWTLJG3vwpgEmt132FJt9Pt5Eq1Fj9PPVVzhq8WgIdcpulE2S3okU6Ec\nzPfPoCUARiQqcABGPDM70MxOSQae2Kz4hXxbsu48M5uevHSd4pfGgULPYLuxrVBcU7edu786uQas\n1HJ68poNiu5ul1sM7X6ConJR6bD/KyTN3knXw0HxGGHwm5KuSo+Zme1lZq9K7p9uZmlXy/WKY5we\nyxWKa+pS/y3pDDN7lZk1m9k4M2tLro2S4jifZ2YHJ2Hnckk/Tbp59ueHkl5hZuckA2bsbmaHZ7bX\nn5slHZicC6OT5aVmdlDymb8j6fNmNjNp63EWg42kFkk6VdJXzOyM5LlXK6qTqZ9Jeo2ZHZ+893KV\n/zN6+2ewME5RlTSLwWXSNv2yv8/Tz3bvVASqS5Pjdrxi0JRfD7Y9Azz3I8W1b/+i+H5SO/v++9s+\nAOzyCHAAENe/XSHpeUUlZJpidEYpfrl82Mw6JX1B0hsy1/KUCgsuSWZ2YvKe/nxb0iFJ97Dry2zv\nuySNl7RS8Yvuv7j7o6VeaDGISaeVGMQk8dPkdrWZ3V9ifan50QYKSR9WdJO8x8zWKQbLODBZNyd5\n3KkY8v4r7p7OJ3aFpI8lx+MDyTV7ZymGwF+pCBAfVOGXdleE1u8pvrMxkt47QLvk7oslnZZsZ7Wk\nBZL+Zmef0907FVW0N0hamuzvimSfUlQVH5L0x2S7VxS1U8k1kq+R9E2L0T9PV6H7pJKBTS5SdCF8\nTjFUf7aLYHH7+j33EicpuvTeLGmWojKaXk/WtZPP03ejcX3ZWYpj16GoUJ/v7o+Xev0Abd7hOXdf\nrjgXjpP048zzO/v+iz8vAIwYNvAfKwe5EbNmxcXFS9z9jJ29HgCAapjZHYoudt+pd1vKZWajFKFp\nvyRMAQAwaLWqwL1PMaEofw0DAOSlUbvQTZH0McIbAKASVQe4pFvOaZK+pcb9YQoAaDylBoJ5U9Jl\ntHh5qB4NLMXdn3f3r9e7HbVgZov6Od7n7vzdAIBKVN2F0sx+qpgIdjdJl9CFEgAAAACGRlUVODN7\njaSVyTDIVN8AAAAAYAhVVYEzs/+QdL6kHsVEprtJ+h93/6fMa7guDgAAAMCI5u41KXjVZBRKSTKz\nk1SiC6WZ7WRaHqBg3rx5mjdvXr2bgQbAuYJycL5gsDhXUA7OFwyWmdUswNV6HjiSGgAAAAAMkVG1\n2lAyGeudO30hAAAAAKAita7AAVVpa2urdxPQIDhXUA7OFwwW5wrKwfmCeqjZNXD97oBr4AAAAACM\nYMP5GjgAAAAAwBAhwAEAAABAgyDAAQAAAECDIMABAAAAQIMgwAEAAABAgyDAAQAAAECDIMABAAAA\nQIMgwAEAAABAgyDAAQAAAECDIMABAAAAQIMgwAEAAABAgyDAAQAAAECDIMABAAAAQIMgwAEAAABA\ngyDAAQAAAECDIMABAAAAQIMgwAEAAABAg6gqwJnZODO718wWmtkjZnZFrRoGAAAAYHjr6JC6uurd\nipHF3L26DZhNcPeNZjZK0m8lXeLuv82s92r3AQAAAGD4mT5deuELpXvvrXdLhjczk7tbLbY1qtoN\nuPvG5O4YSc2S1lS7TQAAAADD36pVUm9vvVsxslR9DZyZNZnZQkkrJN3h7o9U3ywAAAAAjYAAl6+q\nA5y797r7EZL2lvS3ZtZWdasAAAAANASulspX1V0oU+6+zsxulvQSSe3ZdfPmzdt+v62tTW1tbbXa\nLQAAAIA6IsDtqL29Xe3t7UOy7aoGMTGzaZJ63L3DzMZL+rWkT7j77ZnXMIgJAAAAsAsyk1pbpc7O\nerdkeBtOg5jMlPR9M2tSdMe8JhveAAAAAOzauAYuX1VPI7DTHVCBAwAAAHZJZtL48dLGjTt/7UhW\nywpc1YOYAAAAABi5qMDliwAHAAAAoGJ0tssXAQ4AAABAxajA5YsABwAAAKBiVODyRYADAAAAUDEC\nXL4IcAAAAAAqRhfKfBHgAAAAAKBBEOAAAAAAVMRqMrMZykGAAwAAAFCRJtJE7jjkAAAAACpCBS5/\nBDgAAAAAFaEClz8OOQAAAICKUIHLHwEOAAAAQEUIcPkjwAEAAACoCF0o88chBwAAAFARKnD5I8AB\nAAAAqAgVuPxxyAEAAABUhACXPw45AAAAgIrQhTJ/o+rdAAAAAACN5WMfk044gQpcPVR9yM1slpnd\nYWaLzOxhM3tvLRoGAAAAYHh68klp2TIqcPVQiwrcVknvd/eFZtYq6U9mNt/dH63BtgEAAAAMM1u2\nSNu2EeDqoeoKnLsvd/eFyf0uSY9K2rPa7QIAAAAYnrZujQBHF8r81fSQm9lsSUdKureW2wUAAAAw\nfGzZIvX2UoGrh5oNYpJ0n/yZpPcllbjt5s2bt/1+W1ub2traarVbAAAAADlLu1BSgSutvb1d7e3t\nQ7Jtc/fqN2I2WtIvJd3q7lcVrfNa7AMAAADA8HDiidLZZ0uf/ay0dKnEr/sDMzO5e03qlbUYhdIk\nfVvSI8XhDQAAAMCuhy6U9VOLoufxks6TdLKZLUiWU2uwXQAAAADDEF0o66fqa+Dc/beq8WAoAAAA\nAIavdBRKKnD5I3gBAAAAKEvahZIKXP445AAAAADKQgWufghwAAAAAMqSXgNHgMsfAQ4AAABAWehC\nWT8ccgAAAABlSbtQEuDyxyEHAAAAUBa6UNYPAQ4AAABAWehCWT8ccgAAAACD1tsb1Te6UNYHhxwA\nAADAoG3dGrd0oawPAhwAAACAQUsDXG8vAa4eCHAAAAAABm3LlrilAlcfBDgAAAAAg5YNcFwDlz8O\nOQAAAIBBy3ahJMDlj0MOAAAAYNDoQllfBDgAAAAAg0YXyvrikAMAAAAYNLpQ1heHHAAAAMCg0YWy\nvghwAAAAAAYtO5F3qre3Pm0ZiQhwAAAAAAYtrcD19kruhfvIR9UBzsy+Y2YrzOyhWjQIAAAAwPCV\n7UKZBrdsNQ5DqxYVuO9KOrUG2wEAAAAwzGW7UKYVOAJcfqoOcO5+t6S1NWgLAAAAgGHm05+WPvCB\nwmO6UNYX18ABAAAA6NePfiR94QvSmjXxeNMmadw4KnD1MiqPncybN2/7/ba2NrW1teWxWwAAAABV\nWrtWmjVLuuEG6bnnYuqAKVMIcANpb29Xe3v7kGzbPD3q1WzEbLakm9z9RSXWeS32AQAAAKB2enoi\nkO2zT/+vef55ac4c6YwzpLvvlp55RjrhhKjG7b13rF+wIG6nTcuv7Y3GzOTuNZk1jy6UAAAAwAh0\n9dXSIYdEACv2+OPSjTdKv/mNdNxx0uzZEd7MpGefpQJXT7WYRuA6Sb+XdKCZLTazC6pvFgAAAICh\nsmFDBLjTTpM+85m+6zZvll75Suntb5e+8Q3pzW+W9t031h1zjLR06Y4B7oYb4rUYerUYhfJcd9/T\n3ce6+yx3/24tGgYAAADs6tylJUsiNK1cKX3wg1H1eu97o8viAw9IF14ovexl0u23S4sXx/t+9Svp\nySf7bmvBAumii6TLLpOeeiqe6+2V1q/fcb8XXii99KXS174W2/rlLwvr5s+PbpWzZ0fXyL//+0KA\nO/roCG5TpvQdhfJd75L++Z9remjQD7pQAgAAAHXylrdIhx8ela0vflG6807p4oujq+IZZ0innipN\nmhQVsVe8IqpinZ3Sm94kHX+8dP31EahWrYrXtrRIjz4qnXtujB55+unS9OkxFUAatrq6pJtvlr70\nJWnqVOmnP41At//+0uWXx/Ovf7303e9Kt90WI04ecIC0557SfvvFNtIKXDp9QHNzXQ7fiFSTQUwG\n3AGDmAAAAGCE2LZNamqKACZJd9whffSj0tveFstDD0W1as89pfe8R3rta6Wnn5bOPlv63e+kP/5R\nOuigeO+tt0oHHxyVMPeovp1wQoSnQw+NMPe+98VgIk89FWHtc5+LUPWJT0RF7uUvj32ccUYEtJkz\npR/+MCp6t9xSaPfdd0tPPCH9+tcR+K64Qmpt7fvZ1qyJ6+IuuED6+MfjtZ2d0qJF0vjxUnd3ISSi\nr1oOYkKAAwAAAAbpueekr341rh077ri+637zmwhiL32p9POfS2PHSsceK51zjvSDH0QQu+++CEeL\nF0dV7POfjxDW0xMTZE+YMPD+V66MYHbyydKYMVFNe+c7oyL3d39XCI7F1q6NNnR0RNA74ABp8uTy\nP/8vfhFdKq+6Srr22tj/I49E2OvqIsD1hwAHAAAA5OzZZ6VTTpFOPDGuGbvsshjFcf36uFbta1+T\nrrtO+vGPpZ/9TNq6Nbo9Xn+9tHp1hJ+TT45KmBQVralT6/uZynXnnVJbm/T970dXy40bo8vmpEnS\nunUEuP7UMsDlMpE3AAAAMBy5S3/6U1w7dvDB8dzKldKf/yztvnsMp79xYyxXXRXVsosvjirYVVdF\nkJswQTrssBgM5MUvjm6Ll1wSFa7p06NL5fTpcZ1ZVqOFNym6b6a32VEouQYuPwQ4AAAAjEiLF0vv\nfndcl9bVFYFr770jnB1ySFzTdeCBEdAmTIjry173unjvkUdGFaoUM2nu3Pw+R57SbpeTJ/cdhXLM\nmPq1aaQhwAEAAGCXsW5dTDh9zz0xsEdLS1yfNXeutGlTDIu/cmV0eWxvj+H6f/KTCHDPPRdznB17\nbGXXh40EkydHQJ00qVCBa26WRpEqcsM1cAAAAKhKV5d0//0xV1g6zPxQ6u2NwTief15avjwqZvfc\nI/3hD3Gt2ezZMUrjEUfEhNXr1sVAGxMnSi94QXRnPPjg6Oo4c+bQt3dX87vfRRfK178+QtzTT0t7\n7RW3/NpfGtfAAQAAoC7Wro1h4xctkh5+OCpZCxZIL3pRDEN/7LHS+98fYW7MmLh27L77Ysj6Cy+M\nQUBS3d0RuFavjjD229/GJNJdXTHQx5w5fatm6e3q1dGlcfr0CGSHHx6jQl5+ebynv5EYURvHHy89\n9lgE6d7eqMCNHl3vVo0cBDgAAABsl843tnBhhLMnn4ygtXZtVLHWrYvrww47LJYzz4zh9Ftb43Xf\n+14M4LF2rbR5cwStOXNiIuo3vzkeb9gQIzBu2xYDhaTL0UfH/GW77y795S8RCGfOjOfTsDZ9ujRt\nWgzRj/ppbqYLZb3QhRIAAGCE6OmJ8PTMM9K998a8ZXfdFcPgjx0b14utWSPttltUtY44IsJXS0s8\nd8gh0j77xKiKldi0KQLhxIkxAmNLC9WyRvXUU9EF1azQbfXBB+lC2R+6UAIAAIxQ7lEFW7Mmqlxr\n1sSyenXhfrqsXdt32bw5KmUzZ8Zk069+tfRf/xUVr82bo+vilCkR1obCuHFxbRoaX3NzdJ9saoqF\nLpT5IcABAAAMA+5RGXv44bjWa/36CGqrV8foiOkIicuWRRCaMiWW3XePatbUqXF/r73ierSpUwuv\nmTw5bgeqeLW2xvuBwWhqii6UZnGfLpT54VADAADkpLs7QtiSJYXbZ56JecgefDBC1GGHSXvsEcO0\nT5oUXdOOP17ac88IZzNnRoAD6im9Bq65mWvg8sahBgAAuyT36BK4fHksHR3xC+eWLXEd2IYNUmdn\nVLnWr49wtXlzLO7RjXDatBg4Y/fd45dUKW5bWwvXiy1YELf77x9Lc7O0YkXsc9myvoGtqytC2F57\nxYTRe+8dw9mfc05UzaZNq+8xAwYr7ULp3rcCl3arxNAhwAEAgIbR0xMhaP36GE4+DWfFy7JlcStF\nxWqPPaIb4ahRca1OS0uMhjhxYlS5Zs2Kx2PGxGAeZrGPVatiP+mQ6WkbNmyIdkycKB11VMx99tRT\nMVy+uzRjRuzzqKNilMY0sE2bxqAd2DWkXSjTUSjTfx9btlAhHmoEOAAAMKDe3qhOdXUVgkt6W/zc\nxo0x0mB/y4YNEZROOSWqWOk8UqWWrq6obC1bFsPaL14c14a1tERwSkNSusydK510Ut/nWlsJTMBQ\nyE4jkIY5KSrYBLihRYADACAHvb2FitCee1b2/q6u6PK3cWMs3d1xm3b7S5eOjgg+q1bF+tbWCDzj\nxsUvWmnXwTSAZZdt2wqDZ2zdGn9N7+mJ96bdBltb+97P3k6YEMvUqfGe7DJ+fCyrV8fQ9Vu2FEaw\nK17MYptz50p/+7cxdP2sWRHK0q6MAOonrbr19kblOq3AbdoUVW0MHQIcAAAl9PZG+MmGm46OHYdl\nL146OgrXUm3aFLc9PRFKpkyJ0HLyyTGfllQIY8VLGrLWr49l48ZCUEpD0oQJhXA0dmxhmTw5AtTc\nuRGY0m1t2hSfq7W1UJ0qXpqa4tqvSZPil7J0qXUV6+yza7s9APkq7kKZrcBhaFUd4MzsVElXSWqW\n9C13/0zVrQIADHs9PYWQsmVL/PDOLulfZWfMiIrLggXxngkTCl3uspWZceMKkwVPnhzLYOcVSq+L\nquXS3R3tylab0qHYs0OzH3TQjs+lwWrs2LhN/1I9blxs+9pr+1bjsoEsXcaPj2ORLmm4AoDhYKAu\nlBhaVQU4M2uW9GVJr5C0VNIfzexGd3+0Fo0DgF1Fb28EnbSy0tER1/MsWRLd3IqXtWujy9s++0gH\nHlhY9t03BkMYM2bo27x1a8w79eyz0iOPxNxUDz0kLVoU7evtjZAxdmy0Jx1KOl2amuIzr1wZoebI\nIyPAZCtJ6XVR3d1xm3YR7OiIZcKEeK9ZtKenJ7a9226FSYe7uuL5UtWkUsvkyTGYxM5eN2HC0ASm\n1lbpHe+o/XYBIE/ZUSizg5hs2lTfdo0E1Vbgjpb0hLs/LUlm9iNJZ0kiwAHY5XV3R9jaY4+oFK1a\nFeHmqaeke+6J4LJkifTEE9EFLq00pQMwzJoVQeIFL4hg9uIXxwh106ZFJWf9+pgf6vHHY3vXXBNh\natmy6B6Xvn/atPjhmVaLxo+P9owaVVjM+nbp6++2uzuGP1+yJKpme+xRGOb8RS+SzjpLOvTQ2Ge6\n3aHiHmFu7dp4nI4euHVrPJ9ek9XaWhg1EACQj2wXSipw+ao2wO0laXHm8RJJxxS/aN06LmYE0Njc\noyK0fHkME37TTdJtt0V4WL06AllPT0zAu88+0nHHRaUsHRlv8uTKBl446qgdn9u2LdqRVvBWry5c\nr7VhQwS/bduiPT09EXjc+3bpa2kpDDKRPpfezpgRoW3GjPpOzGpW6D4IABhesl0ouQYuX9X+aPbB\nvOiNb5ynl7407re1tamtra3K3QIYKdyjKpROtLt+fYSV0aP7Lml1ZvTo6FL3f/8Xy5o1sQ0pfrhk\nR9/r7o5rt9KR9tKl1OPOzgg4M2ZIRxwhnXaa9M1vRqVs69YIUTNm5FMFam4uTAQMAEA9pF3M04m7\nt24tjHgLqb29Xe3t7UOybXMfVAYr/WazYyXNc/dTk8eXSurNDmRiZn7FFa6PfKTqtgJoMO6F65RW\nrIiq0Zo1hfmiOjsLway/2/XrI5xNmlSoxqQ/KNIlrTKly/jxMbDE3LmFSXPTJf3hkl7jNHp035H2\nxowp/VxLS2wXAACE5ub4mbrffvGz/eKLpYsuqnerhiczk7vX5M+81Vbg7pc0x8xmS3pO0j9KOrf4\nRZ2dVe4FwLDT2ystXRrXZ/3lL3GbLsuXR3DbsiXCUEtLdCWcOTO67aUj+k2cGMFs1qxCQCu+3W23\nfAbsAAAA5Wlujj+cpoOYcC1yPqoKcO7eY2bvlvRrxTQC3y41AiUBDhj+3GMgi3S+qTSArV4dg2Ys\nWxbB7OmnI6Q98USErOwIiW1t0pw5MSz6uHER3hj2HACAXVNzc2F04G3b+Jmfl6ovT3f3WyXdOtBr\n1q+vdi/AyNLbG/8hdnfHwBldXYW5tVpaYtTCSZMG/ktXb29hQIu0O+LKldJjjxXC2LJlMUz8mjWx\nPp3Ad8KEwtDwU6dG5SxdjjoquibOmUM/dwAARrL0koY0wFGBy0cu44tRgcNI1NMT4Wvt2vKWjo74\nN9PUFNdcTZkSoS2dW2vDhghi3d1R4UqHj08nDk5DW2dnPJ/tkrj77hG+Zs2SXvKSQiibNi1eM3Zs\nvY8aAABoFNnKGwEuP7kEOCpwGI7co4tgd3d0GezsLAya0d/S2VlY0gmH0/tplSzV1BShacqU0su0\naVHFKrVut912PuT8pk0RErdti8+QLi0t8f6JE+s7BDwAANi1pb+r0IUyX1TgMGysWxchyH3Hpacn\nAkt3d99h4NMQlb3dsCEC2caNfYNNqWXUqEL1Kh0wo9QydWpMtDxxYmEEw+L7ra1REcsO7DqUf4ka\nN65wn3kWAQBA3tLA1tTEICZ5IsChIps2FYaFX7pUeugh6eGH4y8wU6ZEoEmDVmdn/KMeMya66E2Y\nEPNXzZoVgWfhQumOO2Ikw3TC3uyw72bxunHjYmlt7Rua0vvTp8cwti0tsY80mI0f3/9SycTKO8N/\nXgAAYCSoos4RAAAN6ElEQVRIf4+iC2W+6EI5jPT29l8pyo4KmE4ynL1Ny9ZmcTt6tDR5clSOpk6N\nqtCGDYWqVbZilVar0orVYO739MSkxenQ8IcdJr3udfE51q6Nbe+zT6FK1dxcaP+GDRH67rornjvs\nMOlLX5KOOSbaDQAAgOGPLpT1MWIrcM8/H0Ohd3dHNWnTpggTPT3SqlXRnW/SpAhBLS3R1S6d46Kn\npzBx8ObNhWCS3i/1ePPmviGoVCjasqVQoSpVLRo3bseJhtPbpqYIab29hWu7OjpidME1ayLYtbQU\n5t5Kq1YtLXF/xozYb7ZyVVzFyj43fjz/SAEAAEaybBdKKnD5yS3AuZf+Ut0jCKVBqlT1KX1+06ak\n0UmYyt6m99PH27YVhk9Pl46OCG0LF8bjgw6KADNuXASSMWPivdOmRXhbulRatCgqRtu2RWhLq1vp\nPseOLSxpF8EpUwr3s+taWvoGoeKANG4cJz4AAAAaQ7YCxzVw+cklwI0ZE9dH3XGH9MlPRljJBrR0\nuPQ0SJWqPGUDTloBS0PVtm2FboTpYhYhrHgI9fPPl668Utp/fypIAAAAQKXoQlkfuQS4WbOktjbp\nxBOlG2+U9tyzb0BjqHMAAACgsaSBjYm885VLdHrsscJoggAAAAAaX3YUyt5eKnB5ySXA8WUCAAAA\nu5ZsgOtvvAvUHtEKAAAAQNmyXSglAlxeCHAAAAAAypatwGVvMbQ4zAAAAADKlh2FUqIClxcCHAAA\nAICyFQc3Alw+CHAAAAAAylbcdZIulPngMAMAAAAoW3NzhDYqcPkiwAEAAAAoW3NzhLa08kaAywcB\nDgAAAEDZ0upbGtzoQpmPig+zmZ1jZovMbJuZHVXLRgEAAAAY3qjA1Uc1OfkhSa+VdFeN2gIAAACg\nQaQBjmvg8jWq0je6+2OSZHxTAAAAwIhDF8r64DADAAAAKBtdKOtjwAqcmc2XtEeJVR9195sGu5N5\n8+Ztv9/W1qa2trbBvhUAAADAMFTchZIKXEF7e7va29uHZNvm7tVtwOwOSR909z/3s96r3QcAAACA\n4eXMM6Xbb5fe+lbpy1+WbrtNeuUr692q4cnM5O41qVHWKidTMAUAAABGECbyro9qphF4rZktlnSs\npJvN7NbaNQsAAADAcFZ8DRxdKPNRzSiUP5f08xq2BQAAAECDKB6FkgpcPsjJAAAAAMrGPHD1QYAD\nAAAAUDa6UNYHhxkAAABA2ehCWR8EOAAAAABlYyLv+iDAAQAAACgbE3nXB4cZAAAAQNnoQlkfBDgA\nAAAAZUsn8qYLZb4IcAAAAADKRhfK+uAwAwAAAChb2oWSCly+CHAAAAAAysZE3vVBgAMAAABQNiby\nrg8OMwAAAICyMQplfRDgAAAAAJSNQUzqg8MMAAAAoGzFXSipwOWDAAcAAACgbOkccHShzBcBDgAA\nAEDZGMSkPjjMAAAAAMrGNAL1QYADAAAAUDaugasPAhwAAACAshVPI0AXynxUdZjN7LNm9qiZPWBm\n15vZpFo1DAAAAMDwRRfK+qg2J98m6VB3P1zS45Iurb5JAAAAAIY7ulDWR1UBzt3nu3tv8vBeSXtX\n3yQAAAAAwx1dKOujlof5rZJuqeH2AAAAAAxTVODqY9TOXmBm8yXtUWLVR939puQ1l0na4u7XltrG\nvHnztt9va2tTW1tbJW0FAAAAMEw0NzORd3/a29vV3t4+JNs2d69uA2ZvkfR2SS93900l1nu1+wAA\nAAAwvHz969KXviS9853Se94jPf20tO++9W7V8GRmcveaRNydVuB20pBTJX1I0kmlwhsAAACAXRNd\nKOuj2mvg/p+kVknzzWyBmV1dgzYBAAAAGOaYRqA+qqrAufucWjUEAAAAQONIr39LK3CMQpkPDjMA\nAACAslGBqw8CHAAAAICyFV8DRwUuHxxmAAAAAGUrnsibClw+CHAAAAAAykYXyvogwAEAAAAoWzqR\nN10o88VhBgAAAFA2ulDWBwEOAAAAQNmYyLs+CHAAAAAAylZ8DRxdKPPBYQYAAABQtuKJvKnA5YMA\nBwAAAKBsjEJZHwQ4AAAAAGWjC2V9cJgBAAAAlI0ulPVBgAMAAABQtoMPli64gC6UeSPAAQAAACjb\nC14gXXghE3nnjcMMAAAAoGJU4PJFgAMAAABQMQYxyReHGQAAAEDFGMQkXwQ4AAAAABWjC2W+Kg5w\nZvZJM3vAzBaa2e1mNquWDQMAAAAw/DGISb6qOcz/5e6Hu/sRkm6Q9PEatQkAAABAg6ACl6+KA5y7\nd2YetkpaVX1zAAAAADQSKm/5GlXNm83s05LOl7RR0rE1aREAAACAhmFGiMvTgIfazOab2UMlljMk\nyd0vc/d9JH1P0hdyaC8AAACAYcSM7pN5GrAC5+6vHOR2rpV0S38r582bt/1+W1ub2traBrlZAAAA\nAMNZUxMBrlh7e7va29uHZNvm7pW90WyOu/8luf8eSUe7+/klXueV7gMAAADA8Pab30inny5t3lzv\nlgxfZiZ3r0nMreYauCvMbK6kbZKelPTOWjQIAAAAQOOgApevigOcu59dy4YAAAAAaDxcA5cvxosB\nAAAAULGmJkahzBOHGgAAAEDFqMDliwAHAAAAoGIEuHwR4AAAAABUjC6U+eJQAwAAAKgYFbh8EeAA\nAAAAVIwKXL441AAAAAAqRgUuXwQ4AAAAABUjwOWLAAcAAACgYnShzBeHGgAAAEDFqMDliwAHAAAA\noGJNTQS4PBHgAAAAAFTMjC6UeeJQAwAAAKgYFbh8EeAAAAAAVIxr4PJFgAMAAABQMbpQ5otDDQAA\nAKBidKHMFwEOAAAAQMXoQpkvAhwAAACAijGRd7441AAAAAAqRgUuX1UHODP7oJn1mtnUWjQIAAAA\nQOOYMUN67Wvr3YqRo6oAZ2azJL1S0jO1aQ5Guvb29no3AQ2CcwXl4HzBYHGuoBycL2HSJOnKK+vd\nipGj2grc5yX9ay0aAkj8R4jB41xBOThfMFicKygH5wvqoeIAZ2ZnSVri7g/WsD0AAAAAgH6MGmil\nmc2XtEeJVZdJulTSq7Ivr2G7AAAAAABFzN3Lf5PZYZJul7QxeWpvSUslHe3uK4teW/4OAAAAAGAX\n4u41KXhVFOB22IjZXyW92N3XVN8kAAAAAEAptZoHjiobAAAAAAyxmlTgAAAAAABDr1YVuJLM7FQz\ne8zM/mJmHx7KfWH4M7NZZnaHmS0ys4fN7L3J81PNbL6ZPW5mt5nZ5Mx7Lk3On8fM7FX9bx27IjNr\nNrMFZnZT8phzBSWZ2WQz+5mZPWpmj5jZMZwvKCX57heZ2UNmdq2ZjeVcQcrMvmNmK8zsocxzZZ8f\nZvbi5Bz7i5l9Me/PgaHXz7ny2eTn0ANmdr2ZTcqsq9m5MmQBzsyaJX1Z0qmSDpF0rpkdPFT7Q0PY\nKun97n6opGMlXZScEx+RNN/dD1QMjvMRSTKzQyT9o+L8OVXS1WY2pH90wLDzPkmPqNBNm3MF/fmi\npFvc/WBJfyPpMXG+oIiZzZb0dklHufuLJDVLeoM4V1DwXcV3nVXO+ZEOUvFVSW9z9zmS5phZ8TbR\n+EqdK7dJOtTdD5f0uGLU/pqfK0P5n9DRkp5w96fdfaukH0k6awj3h2HO3Ze7+8LkfpekRyXtJelM\nSd9PXvZ9SX+f3D9L0nXuvtXdn5b0hOK8wghgZntLOk3St1SYpoRzBTtI/sJ5ort/R5Lcvcfd14nz\nBTtar/hj4gQzGyVpgqTnxLmChLvfLWlt0dPlnB/HmNlMSRPd/b7kdT/IvAe7iFLnirvPd/fe5OG9\nipH6pRqfK0MZ4PaStDjzeEnyHJD+FfRIxck9w91XJKtWSJqR3N9Tcd6kOIdGli9I+pCk3sxznCso\nZT9Jz5vZd83sz2b2TTNrEecLiiSjZX9O0rOK4Nbh7vPFuYKBlXt+FD+/VJw3I9FbJd2S3K/puTKU\nAY7RUVCSmbVK+h9J73P3zuw6j1F1Bjp3OK9GADN7jaSV7r5AhepbH5wryBgl6ShJV7v7UZI2KOni\nlOJ8gSSZ2f6SLpY0W/GLU6uZnZd9DecKBjKI8wOQmV0maYu7XzsU2x/KALdU0qzM41nqmzAxApnZ\naEV4u8bdb0ieXmFmeyTrZ0pKJ4MvPofSCeOx63uZpDMt5pi8TtIpZnaNOFdQ2hJJS9z9j8njnykC\n3XLOFxR5iaTfu/tqd++RdL2k48S5goGV87NnSfL83kXPc96MEGb2FsUlIG/KPF3Tc2UoA9z9igvx\nZpvZGMWFezcO4f4wzCUXa35b0iPuflVm1Y2S3pzcf7OkGzLPv8HMxpjZfpLmSLpP2OW5+0fdfZa7\n76cYYOB/3f18ca6gBHdfLmmxmR2YPPUKSYsk3STOF/T1mKRjzWx88jPpFYqBkjhXMJCyfvYk/yet\nT0bDNUnnZ96DXVgyAMmHJJ3l7psyq2p6rowagrZLiovIzezdkn6tGOXp2+7+6FDtDw3heEnnSXrQ\nzBYkz10q6T8l/cTM3ibpaUn/IEnu/oiZ/UTxw7VH0ruciQtHqvR751xBf94j6YfJHwyflHSB4mcP\n5wu2c/cHzOwHij8y90r6s6RvSJoozhVIMrPrJJ0kaZqZLZb076rsZ8+7JH1P0njFCLm/yvNzYOiV\nOFc+rvi9doyk+ckgk39w93fV+lxhIm8AAAAAaBDMZQIAAAAADYIABwAAAAANggAHAAAAAA2CAAcA\nAAAADYIABwAAAAANggAHAAAAAA2CAAcAAAAADYIABwAAAAAN4v8DJOd1rdnMj+sAAAAASUVORK5C\nYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig, ax = plt.subplots(1,1,figsize=(15,4))\n", "t=0\n", "ax.plot(ssh[t,:])\n", "ax.set_title('ssh: t={} : {}'.format(t, run))\n", "#ax.set_xlim([500,650])\n", "#ax.set_ylim([-50,0])" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": true }, "outputs": [], "source": [ "#" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.4.3" } }, "nbformat": 4, "nbformat_minor": 0 }