{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "Tracking density 1023.3 through various wind events. " ] }, { "cell_type": "code", "execution_count": 16, "metadata": { "collapsed": true }, "outputs": [], "source": [ "import numpy as np\n", "import netCDF4 as nc\n", "import matplotlib.pyplot as plt\n", "from salishsea_tools import (nc_tools, viz_tools, gsw_calls)\n", "import numpy.ma as ma\n", "from matplotlib import animation, rc\n", "%matplotlib inline" ] }, { "cell_type": "code", "execution_count": 17, "metadata": { "collapsed": true }, "outputs": [], "source": [ "mesh_mask = nc.Dataset('/home/vdo/MEOPAR/NEMO-forcing/grid/mesh_mask_downbyone2.nc')\n", "new_domain = mesh_mask.variables['tmask'][0,:,334:898,114:398] \n", "def calc_rho(Sal, TempC, P):\n", " sqrSal = np.sqrt(Sal)\n", " R1 = ((((6.536332e-9 * TempC - 1.120083e-6) * TempC + 1.001685e-4)\n", " * TempC - 9.095290e-3) * TempC + 6.793952e-2) * TempC - 28.263737\n", " R2 = (((5.3875e-9 * TempC - 8.2467e-7) * TempC + 7.6438e-5)\n", " * TempC - 4.0899e-3) * TempC + 8.24493e-1\n", " R3 = (-1.6546e-6 * TempC + 1.0227e-4) * TempC - 5.72466e-3\n", " SIG = (4.8314e-4 * Sal + R3 * sqrSal + R2) * Sal + R1\n", " V350P = 1.0 / 1028.1063\n", " SVA = -SIG * V350P / (1028.1063 + SIG)\n", " rho = 28.106331 - SVA / (V350P * (V350P + SVA)) + 1000\n", " return rho" ] }, { "cell_type": "code", "execution_count": 18, "metadata": { "collapsed": true }, "outputs": [], "source": [ "Jan04 = nc.Dataset('/ocean/vdo/MEOPAR/completed-runs/SalishSeaLake/LakeJan0.4/SalishSea_1h_20170101_20170107_grid_T.nc')" ] }, { "cell_type": "code", "execution_count": 19, "metadata": { "collapsed": true }, "outputs": [], "source": [ "rc('animation', html='html5')" ] }, { "cell_type": "code", "execution_count": 33, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 33, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXMAAAD8CAYAAACFK0QrAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztvX+UHOV15/25QlYsW7IyhKCgH2NGCswG69UQR1EIcRYR\ny7JQTEjkjcEbE97X7OF4l1gyRq8F8Rs5Sw57kIOsoMR513pt1gQHAwkiRo4QkjgBnxyMsSCMkMAT\nEIJBww8FM5ZHXoIsdN8/qp6ep6urqqu6q6ere+7nnD7dXfVU1e0e6Vu373Ofe0VVMQzDMDqbKe02\nwDAMw2geE3PDMIwuwMTcMAyjCzAxNwzD6AJMzA3DMLoAE3PDMIwuwMTcMAyjSUTkBRF5SkSeFJG9\n4bZzReRRt01EliYcu1JEhkTkORG5ztt+qojsFpFnw+eeNBtMzA3DMIrhQlU9V1WXhO+/CPx3VT0X\n2BC+r0JETgG+DFwEnAN8XETOCXdfBzyoqmcBD4bvEzExNwzDaA0KvCd8PQt4OWbMUuA5VX1eVY8D\ndwKXhPsuAW4LX98G/G7axaY2bW4BnHbaaXrmmWe22wzDSGXo2VcB6D/rF9psyeTm8ccff11Vf76Z\nc/zGsnfqj944mWns00/99ADw796mraq6NTJMgT0i8jbwlXD/Z4AHRORmAsf5/JjTzwVe8t4fBn4t\nfD1bVV8JX78KzE6zsxRifuaZZ7J37952m2EYsawbvBSAZSysbLt54K52mTPpEZEXmz3Hj944yR3f\nTtXGCue+9/C/e6GTJD6gqiMicjqwW0R+APwn4BpVvUdEPgZ8DVjeiL2qqiKSWnvFwiyGkYIT8rjt\ni7dvYGDt5gm2yCgjqjoSPh8B7iUIn1wBbAuH/F24LcoIMN97Py/cBvCaiJwBED4fSbOhFJ65YbSL\nJLE2jKyIyLuBKao6Fr5eAdxAECO/AHgI+C3g2ZjDvw+cJSJ9BCJ+GfCfw333EdwQbgqfv5Vmh4m5\nMSkpSsQHb7mmkPMYHc1s4F4RgUBT71DVnSJyDLhFRKYSxNyvAhCROcBXVXWVqp4QkT8CHgBOAW5V\n1QPheW8C7haRK4EXgY+lGWFibnQscYIcF8tuRLgzxcQHcp/W6EJU9Xli/jWo6j8DvxKz/WVglfd+\nB7AjZtwPgQ9mtcPE3OgIsgpysx63TWwanYqJuWFgIm50PibmRkdw88BdhU5Wmngb3YaJudE1mEAb\nkxnLMzc6gnpeuQm5MdkxMTdKT5bwyrrBSy1n3JjUmJgbpefmgbsye94m6MZkxWLmRsfgC7qJtmFU\nY5650XWY0BuTEfPMjY7BRNowkjExN0qLibdhZMfEvAkG1m5mbGFQ4P75Nde22ZruohkhtzRFYzJi\nYt4EunyUGe02oksowgs3ETcmMybmBXDs0Kx2mzCpMNE2jFpMzHOyePuGyutjh2ZZeKVAdg3312xb\n0TtU9d6E3DDiMTFPYfH2DRUxiRMaE/JiWdE7FPs9g4m4YdRDVFN7hAaDRF4AxoC3gROqukRETgXu\nAs4EXgA+pqqj4fjrgSvD8WtU9YG0879r9nz9hc9/BoAZfUer9u27+IY8n6dl+B65T1ns61SyxMpN\nyA0fEXk8Q4PlVN63eJrmaOjc9PUmgjye+YWq+rr3/jrgQVW9SUSuC9+vF5FzCPrYvQ+YA+wRkbNV\n9e1GDIwT0XYIqH9N36bF2zek2pN0E4ieczKxbvBSdg3314RQfEzADSMfzYRZLgGWha9vI2hauj7c\nfqeqvgUcEpHnCLpSfzfpROfMn83eNdeWRrjr4Wxy9rrnfRffkCrecefoBup9Zl+0AxGPH2cCbhiN\nk1XMlcDDfhv4iqpuBWar6ivh/lcJmpoCzAUe9Y49HG6rQkSuImxw2tsb/O+OE0NfKMtEnID526Je\n567h/tJ9hrxkvVHF4WLh7ntJ88oNw8hPVjH/gKqOiMjpwG4R+YG/U1VVROoH36uP2QpsBViyZEnl\n2KjX6/Dft1sUGxG1ThevuM+c5zNZSMUwWksmMVfVkfD5iIjcSxA2eU1EzlDVV0TkDOBIOHwEmO8d\nPi/clou0kEW9OHURZBHsNBviJvbWDV7aUuFqxnPOQ1E3JhNxwyiOumIuIu8GpqjqWPh6BXADcB9w\nBXBT+Pyt8JD7gDtE5EsEE6BnAY81YlySl+5vK0LU84hg1us5oZqI+iILtmyqyQLyafWvAhNlY7Ij\nIqcAe4ERVf2IiAwA/xOYQZDt9weq+uPIMf0EGYGOBcAGVf0LEfl94E+BXwKWqureejZk8cxnA/eK\niBt/h6ruFJHvA3eLyJXAi8DHAFT1gIjcDTwNnACubjSTxZGUSRL3viiK8vxbKXQLtmyqvD52aFas\noLdKyE3ADaOKtcAzwHvC918F1qnqwyLySeD/Bv7EP0BVh4BzoXIzGAHuDXfvB1YDX8lqQF0xV9Xn\ngYGY7T8EPphwzI3AjVmNyEOasE9mokLui3ge4Y3+kjDRNox0RGQe8NsEmvfZcPPZwHfC17uBB4iI\neYQPAgdV9UUAVX0mPHdmOzp6BWi9VMCs3nWZJlez4Hvkad54I0Js4m0YNZwmIn6YY2uYwOH4C+Bz\nwExv2wGCNO1/AH6f6nnEOC4DvtmMkR0t5pNhQY4v3FHS4uSGYSTzw7dncPvo+RlH3/160gpQEfkI\ncERVHxeRZd6uTwJbRORPCOYRjyedXUSmAb8DXJ/RoFg6VsyLXF7fLcLfLZ/DMDqI3wB+R0RWAe8E\n3iMi31DVTxAkiyAiZxOEYZK4CHhCVV9rxpCOFPPJHCt33vixQ7MqpXet4JdhtAdVvZ7Qow4983Wq\n+gkROT1M5Z4C/D8EmS1JfJwmQyzQoQ2djx2alZilsXj7ho4Q+6wpi3FCvaJ3iBl9R3l+zbUm5IZR\nTj4uIv8K/AB4GfhfACIyR0R2uEFhuveHgG3+wSLyeyJyGPh14B9FJLVYIWSsmthqlixZonv31k2j\nrKKTqxjmyRiJ5pA3M7lpGN1AEVUTf+F9p+on7vhQprGbzr2766omdgRlreXiE61Tkualr76gdpsJ\nuWEYUTpSzLOEUQbWbmbwlmsmwJrGqVcG1scE3DCMNDpOzLOEVxZs2QQLT1bS+jotrmzCbRhGXjpq\nAjSp3nk0pPL8mmuZebC8H823d9dwf2KrtDjWDV46IfVeDMPoLDrOMwcypeSVOcQSd1PKEnLxRdx/\nbZ68YRgdI+ZByGRWu82YEMzzNgwjL6UX84G1m4MXC9trR9kwb9wwDJ/Si/ngLdeMC3qXEa3X7oda\nTKwNw8hDqcU8KuKdlpWSRJYc+FZ3JTIMo7sob8oHMLbwZOXRjbjMlGh2i7/fMAwjC6X1zAfWbq7E\nybvFI09i3eClrOhN3mceumEY9SidmPuhlZkHp5Q6xbAZsnrdJuSGYWShdGLuM+at4oT2eejRvPBm\n677kCZ+YZ24YRhZKJ+a+Jz6wdnNVvHzBlk0tEfR1g5dWYtXRDJOsLN6+odB2bRYvNwwjD6UTcx8n\n5G5pfismQp1oRisYrujNJ8r+6s00b9rCK4ZhtIJSi7mPE/K8xbMWbNnE6gseq9p288BdVZ50I9QT\nZVtubxjGRFJqMS+iWFZUyCE9ewTqi2/eEIgbb6JuGEarKLWYu9WfE51n3qp4tcXBDcNoFaUWcxif\nEPWzWtz7LKGWvevfzxtX/6SpkEqryOKpFx2usfCPYXQnpRfzZnlo5/rCGjzfPHBXJu+6VSKZp3do\nFgaH5wMw0PtSU+cxDKP9dISYL1u5EVZVm5p1ArRefDwPZcv5TrqxZLXx9tHzubznEQaH55ugG0aH\n0xFi3iitiFG3Ktc9L3GC7T7v4u0bOHZoVqyd0V8XrRR0/zqX9zwC2K8Aw2gVpS60lUY0hh6llZON\ny1ZubNm5s5Dked88cBc3D9zFvotvSL3huHGO20fP5/bR81v2nTkhN4xuRUROEZF/EZFvh+9/X0QO\niMhJEVmS59hw25+JyD4ReVJEdonInHo2ZPbMReQUYC8woqofEZFTgbuAM4EXgI+p6mg49nrgSuBt\nYI2qPpD1OnGM9k8Dsme0FCFKu4b7Y5ft3zwAC9jUEg+9nt1Fh3iiXnoRouvOd3nPI1zeM77dPHKj\ny1kLPAO8J3y/H1gNfKWBYwH+XFX/BEBE1gAbgE+lnSRPmCV6weuAB1X1JhG5Lny/XkTOAS4D3gfM\nAfaIyNmq+naOa1VYsGUTM5YfZXXvELuG+yv9P1vZsHnbw0sThXrx9g3M6Bt/7Wi2Xks7KFLIfRF3\nmIAbkwERmQf8NnAj8FkAVX0m3Jf72PD4H3vD3g1oPTsyiXnCBS8BloWvbwMeAtaH2+9U1beAQyLy\nHLAU+G6Wa0WZ0Xe08npF7xC7iPeYfbJmncSxa7ifmQensGzlRkb7p6HLRzMdt2DLJnp3nOChnesb\num4WexuZgE3KWBkcnl/lOec51+2j51dtNy/c6HJOE5G93vutqrrVe/8XwOeAmQ2cO/FYEbkR+EPg\nKHBhvRNl9czjLjhbVV8JX78KzA5fzwUe9cYdDrdFDb0KuAqgtzc+3WRg7WZYXr1tRe9QZYJv5sEp\n9Awdr+xz4ttsTvnYwpPo8p8AQX663zDCndttc3b0Dh1vWMjzkCboTmzBF9z4Cc48guufF+K9eBNw\no5P48fF3Vv2/rsPrqhob9xaRjwBHVPVxEVmWx4Z6x6rq54HPh2HrPwK+kHa+umKexVhVVRGp+zMg\ncsxWYCvAkiVLYo/tGTrO8MJZVd45hILaC1wQvN813I/sCVzDY4dmBfsapHLu6LYY2hVaSfbiAwEP\nPOVqwU0T23r55gO9L9UIer1zGsYk4TeA3xGRVcA7gfeIyDdU9RMFHvu3wA6aFfOkCwKvicgZqvqK\niJwBHAnHjwD+//x54bbcTN8/wsz+PsaYxS76E0V1Re8QfLKRK+SnMjE6ULtv2cqNDK+aGlsPxpGW\nUpiFrLHtooXWhNswalHV64HrAUJnd11GIU89VkTOUtVnw6GXAD+od766Yp50QRH5c+AK4Kbw+Vvh\nIfcBd4jIlwgmQM8CktUthTcXBdEZ55n73ev9161m28NLx+PhA0H4Jy6cs2QjDD+8NPVcWYW7nmgX\nKa4m1IZRLCLye8BfAj8P/KOIPKmqHw5TDL+qqqvqnOImEeknSON7kTqZLACimj064on5R0Tk54C7\nCYISLxKkJr4Rjvs8ga98AviMqt6fdt4lS5bo3r17a7a7fO7hVVMrgu6yWWJDLyXA97wbmYRNE/Fm\nRTcaUy/TalbDyIqIPJ4Uw87Ku86ao7/4pf+SaexTv/NnTV9vIsi1AlRVHyLIWkFVfwh8MGHcjQSZ\nLw1z0bw1sGguo/3TmHkQxqhOSRyjOpYeN5nRDoEvauFNFuHOU1vFz0SpzDHc0pyNhmGUh45Zzu9E\nfDx7ZVqNoHcDLgtlIMMkbh5P3Y0d6CWI91/cgHGGYZSWjhFzCIR8+v6RSiwdksMurSRpdWizXrnV\nLzEMo1FKW5vl/sNbgHFP3M8nn75/hJ6h4/R/8UV6d5ygd8cJjh2aVRF2R4480twMrN1csy1a8yQv\nt4+eb0JuGEZDlFbMgcoiHF/InVc+ff9I5Xn6/hH6v/giMw9OmRBBX9E7FNv9aGDtZhZs2RR7zSx2\n2ISkYRiNkiubpVUkZbP4uMwWJ+JJvPzRPsYWnpyQbBeXsjjaP40PffLR+gdkPCdkr9duGJMNy2aJ\np/Qx84vmrQleLKqpCBBL4MXXTo5GPeMixH1G31FG+3sK7VF67+/+RfjKxNwwjOyUXsyBqgnPegSe\n+1zqZbvUC3tkEfsVvUPsWt7PzD097D54HmMLT6au/jQMw2gVpY6Zw/hEqCOPsDfDruH+qkcSK3qH\n+NAnH62EWdwkaLOx+lY21zAMo/voCM/8oZ3rq7r7RCdB4+gZOs7YwqkcO1RMLnpc5cQoqy94jMXb\nNyB7ejh16N0c2B8Ukhz63Hurxs3oO5p4Dpdnbr05DcPIQ+k9c8dDO9dnLjHrRN4tNIpmuDRLmrfu\nZ7q8uWguL3+0r2r/zINTmPtff8SB357ddCqjYRiGoyM8cx/fS8/ioc88OIWxhSerBL2oBUZO0KNe\n9uoLHquU563hAuCTwbEDazeHE7bvT6y2mGfJvmEYk5eOE3OoDbvE4SZCg9ouYT2X0GPOs2o0KRzi\ne+aNVHB0ZXv9Vnj+efxwC2AhF8MwUulIMYdA0Ctpiwn4gg4xfUMP9tAzdLzSoSgaX08TaH9fUtgl\na0YMveP55VEshm4YRhY6VswhW5hlPFWRiqjDeJ2XYAy8TB9E8sWzetx+jXX/hpBl0tRRL6XRifrt\no0GWi8XaDcPw6Wgxr7ca1B/35qK5NfVdHNHCXdHFRllDKHEt55ygF9VMw4VdFm/f0La2dYZhlI+O\nFvP7D2+pG2pxZBX+OJoR4qLLCDgPfUXvUGpzZ8MwJhcdk5qYRHRRUTPUxNQ9WlmB0TAMo1k6XsyL\nwg/BJOWl11sN2g7WDV7KusFLY0vyGoYxeejoMIsjT7ilCCaimXRS+CRpmf+HPvkoC7ZssmqLhjFJ\n6RrP/P7DWwoNuRS9ajQLu4b7E1eFLtiyqW69ltUXPMaCLZtaZZ5hGCWmKzzzIvDTFrOQtPqzGdyk\nZhyrk1aU1owL6sMcOzTLvHTDmER0jWfumGjvvGwxdAhuCt3W6NowyoqIvFNEHhORQRE5ICL/Pdz+\nZyKyT0SeFJFdIjIn4fhbReSIiOxP2H+tiKiInJZmR9eJOeQT9CJK6pZV0A3DmBDeAn5LVQeAc4GV\nInIe8OequlhVzwW+DWxIOP7rwMq4HSIyH1gBDNczomvDLFkmRZ2Q5w2xxNFs2KXRcgCGYbQXDXpv\nHgvfviN8qKr+2Bv2biC2R6eqfkdEzkw4/Wbgc8C36tlRSjF3ItxsyCRN0IsU8kbxBdyt5lywZVNN\nOYBGRH3XcD8MFGCkYRh1EZFTgMeBXwS+rKrfC7ffCPwhcBS4MOc5LwFGVHVQROqOL6WYO+KEOKvA\nR4/NUscFxkvmQmOlcvN42Ct6h6qqJkJtI+dGOg7dPHCXCblhpPD2W6fkyVg7TUT8jvNbVXWrP0BV\n3wbOFZGfBe4VkUWqul9VPw98XkSuB/4I+EKWC4rIu4A/JgixZKKUYu571NGY9rKVGyuC7IQ9VrgX\nza2M889RL0ZeZHNmR5rAr+gd4uaLi1uSb8v7J56sN1z723Qsr6vqkiwDVfVHIvJPBDFwf0Lzb4Ed\nZBRzYCHQBzivfB7whIgsVdVX4w4opZhDtaC7ErWOY4fey8yDU6qE/U1PvJupwzIR+J560f/BnbCY\ncBjGxCAiPw/8NBTy6cCHgI0icpaqPhsOuwT4QdZzqupTwOneNV4Alqjq60nHlFbMIRD0uCYUM/qO\non0wvHAWvdTPRonz0OOIeuV5mlhkZaIqHfre4q7hfquwOMHYzXRScQZwWxg3nwLcrarfFpF7RKQf\nOAm8CHwKIExR/KqqrgrffxNYRhDOOQx8QVW/lteIUot5lKTmzK7ELVR76GnjoqSFV4pqCt0INw/c\n1VDc3HF5zyNc3vMI6waHTGBaQPRvY9/x5ENV9wG/HLP9ownjXwZWee8/nuEaZ9YbUzfPPCUh/lQR\n2S0iz4bPPd4x14vIcyIyJCIfrneNLKQJahYhr0daxUR3/SKW+C/enpRqmoxb4t+IULiSudDYZKqR\nHRNyo51kWTSUlBB/HfCgqp4FPBi+R0TOAS4D3kcwCfDX4c+PhhlbeLJGyI8dmkXvjhPNnLaGeoLu\nrtssi7dvaEjUoTHBiAp6o9fudtzfJev3YzdHo0zUDbMkJcQTBPSXhdtvAx4C1ofb71TVt4BDIvIc\nsBT4bpGGQzAx6peuzUK9+Hm0+XMcRcXSnWjkjWc7QW9UTFb0DjE4PB9gUvcUjX5/K8IuUVlW9JqQ\nG2UjU8w8LiFeRGar6ivhkFeB2eHrucCj3uGHoXaWUkSuAq4C6O3tje4GgjTE0f5pzOgbrdp+7NAs\nZh6cUunjGU09zBNqibsZjPZPyyTqRdGsqDsaEZjoMVEhkz09DN5yTe7zlpV631EzOfrW+cloJxI4\n3hkHhwnxwKeBf1bVn/X2japqj4j8FfCoqn4j3P414H5V/fuk8y5ZskT37t1bs31g7eaqlESoFvKH\ndq4HgjzzqKedN3YePT66MjRN1IucHG222mFRHqPrNepCNNseXtrRVRgziXgJzmnUR0Qez5r3ncTP\n9M7Xues+k2nsobXrmr7eRJCr0Jaq/ghwCfGvicgZAOHzkXDYCDDfO2xeuC0XA2s31wio7Omhd8cJ\nBm+5piLkEL8qNG8Braj4Rz32tHh6UZOjENwYmolpNzpRGuX20fOrYu2rL3iso0ML9b6ThlfaJnB5\nzyOVUJYj+t4wiqRumCUpIR64D7gCuCl8doVg7gPuEJEvAXOAs4DH8hjlC3mVp3px8jEP7Vxf46Hn\nDbnUjq2u31Iv9NKKvPRGiROaIsR43eCl7Brujw2/ZDl/q7zVpGu32jv25y/crxnDaAdZYuZJCfHf\nBe4WkSsJEuI/BqCqB0TkbuBp4ARwdVi3oOW4VaNZlu9nEfnxMXNzFeQqIi/deecurFRE3LrZidMi\naIXopn2eiVpqf/PAXVWed9zE8mSebDZaT66YeauIi5n77c/yxmrjYuhJ5PHc31w0LupZJkbLvnLU\nFx8/pFIG8ohr3ptTXA/XZsW8npAbxWEx83hKuwLUr144sHZzZs902cqNkCNenicUE4ybW8l2mYhM\nl1Yy0PtSS+O4zaxebeWvByfkzf4asLCKUSZKK+a+ULZaNLOWx3VjRvv7Mp23yDIAC7ZsKk02SVTE\n4rx6J5RZJh4bEcW4a0YzcJohbpl+dFv0OibuRjspbZgFxkMtfiaJn5IYJa4olyNOqBtNZ8wTbilK\nzIuMnTucV55X/JLEPE5Mo3nrjX6O6C+Igd6XqsQ1yw2mSLFPwtIRW4+FWeIprWcOQazcLRxyojm2\ncGouL3X6/pFMDS2WrdxYJe5Jwp435bEo73xG31HGKCb90dGoqNUTSifgfnaPa8TR6K+LuDh0mnAO\neOvQWhWyMeE2ykSpxRyoEvIZfUdT87kf2rm+yjtP8uCzHBsXeokT8iyxc9/mZoS91WGWqBfdSLu6\n1BWUbep+5Ivu4PD8TOEQ/4Zlom10AqUW82UrN0L/tIqIOy/32KFZmSZFXfOKrK3mfPF3vwhiKhFU\npSnmjefH3YyyCPyxQ7MqYaeiRN1ldbjWde7G1OgNpxOyONykb9TWaBjHBNzoNHKtAJ1opu8foWfo\neJUAOkEfW3iSZSs3ctG8NVVt4x7auZ7p+0cqj7xhEf88UC3co/3TmhLyJJJWkLrt0X1+2mYz7Lv4\nhiohj5Kl4JSjE4TckZQD7j8Mo9MotWfuPOpo3DzKm4vmsmzlxooA+y3npu8fiW0MHT0+jjmRm4Fb\n4p9mSzPkKQlQlJde75pxOdlxFNWuLqmUgXVKMox0Sp3NUo/o4qBmen/6+ebudZzIt0LIXVij0fou\nzQq6P1cQ/XzRkEueOHpeYTchN7Jg2SzxlNozr4fz3POs+EzCvxEk3RRa5ZE76k3wJlFkLD349eGF\nkqiuN5PVU4fs3npUxE28DSM/HS3mDj8c04x37uPfHPIs4c9L1PNt1ktvBFekzDF9/3gufRErXZ2o\nxzWW9oXcRNwwGqcrxNxRlJD7tFLIi6aZejb+PIOjZ+h45fP7+fJ5vHOH31i6mmCS1YTcMJqj1Nks\nebn/8Jamwy0Od2PoGTqeuzVdHvyMlaTslUZoJOPFT+F02UA9Q8crmS6+XXkyXaC2Prpj38U3mJAb\nRgF0lZhDvoVCeejdcSJTw+dmmHlwCjMPTqlcyz3SxqfRrKDnuVajFJVmaRjtQkTmi8g/icjTInJA\nRNZG9l8rIioipyUcf0143H4R+aaIvNPb92kR+UG4/4tpdnRVmMURFzIoiiKrJcYJZNpCKL9phzvW\nP0ecXUVMjvrhlqJwnn1ZiocZRhOcAK5V1SdEZCbwuIjsVtWnRWQ+sAIYjjtQROYCa4BzVPXNsBfE\nZcDXReRC4BJgQFXfEpHT04zoOs+8SOq1kmsFA2s3JxYMG7zlmhpvfWzhycojjQVbNmX2gp13Hhey\naibU4tj28FJO/fK7LbxidAWq+oqqPhG+HgOeYXzp+Gbgc0BaDvhUYLqITAXeBbwcbv+vwE2q+lZ4\n7iMJxwNdLOb3H96SeRl/GnGCniUE0ioGb7mm6pHXs80j6O6zp3nljQj66gseY8nGJyzEYnQdInIm\n8MvA90TkEmBEVQeTxqvqCHAzgef+CnBUVXeFu88GflNEviciD4vIr6ZduyvDLEUTXUDkrwSt1xc0\njaRjBm/JF/d/fs21NQ2wfZGPimbW0Mv9h7eklhVultUXPMbi7RvMQzcmnFPeyjUPdJqI+Ksat6rq\n1uggEZkB3AN8hiD08scEIZZERKSHIJTSB/wI+DsR+YSqfoNAn08FzgN+laBN5wJNWOnZ9WLue+fN\nxNHjVoTGiTo0n8bYSJw7Kdae5v020vCimRWhhtGhvF5vBaiIvINAyP9WVbeJyP9BINCDIgIwD3hC\nRJaq6qveocuBQ6r6b+F5tgHnA98ADgPbQvF+TEROAqcB/xZnQ9eGWeJoNuzi0vXqUVQYZmDt5qaO\nj/L8mmsrjyyMV45sXZ69eeVGpyOBWn8NeEZVvwSgqk+p6umqeqaqnkkgzO+PCDkE4ZXzRORd4Xk+\nSBBzB/gH4MLwGmcTLM1+PcmOrvfMo0QFvRFv3a/hAtUTo9H4crNhGCfoSSGUerjr+577gi2bKtvz\n9FeNYl65YQDwG8DlwFMi8mS47Y9VdUfcYBGZA3xVVVep6vdE5O+BJwhCM/8CuBDOrcCtIrIfOA5c\nkRRigQ4vtFUkrr5LI6tIkwpyJdHK1aRpQu+HXLKkRUZvJK0IsVjdcCMvRRTaetfs+XrWZZ/NNHbf\nls92RKEtE/MEGo2v5xH2VtV6SQtdxPVVzWJPK2PlJuhGHkzM45l0YZas+BUZ8xD17N9cNDcxDFPk\npKmPK164VPOJAAAaj0lEQVSVJupZrxfXdShNyHcN9yN7eugZOs7wqqmsvuCx1PPvGu5vWzs5w+gm\nTMzrEDdpmkfgoz1EowuPnLgXsbI0KrzR0rJpNV+KaDoNgdDvWt7PKD08v2Y8bJPUVHkiq0MaRjdj\nYZYCaCblMakwWN7a6VnF2K9+2AxZwizRkrd+do4uH63pP9roRKwxubAwSzwm5gVTVE0YJ/JZRb0o\nzzorzXYcSvLU42qeG4aPiXk8FmYpmKKKfLnwjN8oApJj3UV53Fnxl/HnnQxNEvJGzmUYRoCJeQto\ndPI0iWCxUu12J/LtbpxRT9jTxNvHsloMo3HqLlFMqtUrIqeKyG4ReTZ87vGOuV5EnhORIRH5cCs/\nQJkpotBXGq55RP8XX6R3xwlkT0+hDS4aYddwf8PVFLOKvmEYtWTxzGNr9QL/J/Cgqt4kItcB1wHr\nReQcgnq87wPmAHtE5GxVfbs1H6HcNJsNU4/xcMwIMJeeIS/98WAPYwtPTng8HRovj7t4uH5apWEY\ntdT1zFNq9V4C3BYOuw343fD1JcCdqvqWqh4CngOWFm14J9Mqj91vdefTael/K3qHWDd4KYu3b2DB\nlk0s3r6hJs3SMIxqcsXM/Vq9wGxVfSXc9SowO3w9F3jUO+ww44XajRYTre7o8tcneoK0GXyv3rfZ\nCbqlMhpGLZnF3K/Vq6o/Dss6AqCqKiK5chxF5CrgKoDe3t48h3YFRZXmjcOFXCBSc51xD71ThD2O\nGX1H0b5qcfcxoe8cBofnc/vo+ey+9Tz7mzVJJjGP1uoNN78mImeo6isicgbgWhqNAPO9w+eF26oI\ni7tvhSDPvEH7u4KksEsRIu/37/QrOMaFXjpV4KN2O6G3uHtncHnPI+zmvHab0fFkyWapqdUbch9w\nRfj6CuBb3vbLRORnRKQPOAtIL9BhxNJM67u06o9Jddb9TJhOi7MbnclA70sM9L5kXnkBZOme4Gr1\n/paIPBk+VgE3AR8SkWcJumXcBKCqB4C7gaeBncDVkzWTpSiaFfRGG1FHxb3TRL7o5h6GUWbqhllU\n9Z8BSdj9wYRjbgRubMIuI0IzVRxdga+0Gut58AW9zKEZXT7K4u0bkD09HeH5DQ7Pr3o/0PtSzf7o\nNsNw2ArQkuGaZDge2lnd3LmRcgHjIZfaSVForvxuJ2TJOFE/dmgWvTtO1HynZSAq5En7Wyno9W4m\nRrmxQlslZNnKjTXbGumAlIR/s2hVj88yC/yxQ7NyN7JuFb6A3j56fuV1O0sb1Lth+Da3Q/Ct0FY8\n5pmXkCTPschaL7X9S8dDMEUIe5mzZWb0HS1F+CUoX3B+3XETTT2BNo+9nJiYdxBFlgaIa0odTWGE\nYj12J/BlEXU//DJRuen16s+UpdhYu71vIz8m5h1Os4uP4rz0Vre2K1ucPboIaSLz08si3nGYiHcW\nJuZdRNRzzyruaV66T5HCXjYv3acdol42TMg7DxPzLiav117PS/eJLjxqVNzL5qU7XLGvZj3nLGV9\n07ouufrweTswxWXHmEB3Nybmk4SseerRQl3+gqO0PPVmxL2MeetOSPMKet6a7Fm7LuXpwOSE/PbR\n87m855Fc9mTF/Ttqdc3+TkBEbgU+AhxR1UXe9k8DVwNvA/+oqp+LOfYFYCwcc8JlzYjIXYCrOPez\nwI9U9dw0O0zMJxlZRD0adnFkFXaorgOThzIJuxP0xds3VImp85IXbNkEwOoLHuPynkeqUgvbwbrB\nSyvi7Wxxz5f3PFJ4jroJeYWvA38F/I3bICIXEpQDH1DVt0Tk9JTjL1TV1/0Nqlq5y4vIJqDufwbL\nM5/EZAm9RAU9K1GxbzQM025Bd+TtTZqniXWjRIU7y9huCLWUMc88LA/+beeZi8jdwFZV3VPnuBeA\nJVEx9/YLMAz8lqo+m3Yu88yNVBpdrOR6lo7fDPLlsZdFxB1+/DoLE9ECL0nE/dCKH2rpBiEvilP+\nXfPULDpNRHxvc2tY9TWNs4HfFJEbgX8H1qnq92PGKUE3treBr8Sc9zeB1+oJOZiYGy1mvK3d+LaX\nP9pXed3uZtR5SBL0vEIfh/Pk/ZuA791nuTnExcdNyAvh9QZ+CUwFTgXOA34VuFtEFmhtKOQDqjoS\nhmF2i8gPVPU73v6PA9/MekHDmFDm3HOo8vrNRXOrQjJO3Mua5RIn3FmFPMtEatKYesLuC7kJdyk4\nDGwLxfsxETkJnAb8mz9IVUfC5yMici9Bi83vAIjIVGA18CtZLmhiPkkouptRUUzfP8L0/eOi7trc\nQWcJetK4onLV63nmfjjFqiuWgn8ALgT+SUTOJogzVsXFReTdwBRVHQtfrwD8fzDLgR+o6uEsFzQx\nNxJpdBFSI/iiDtM6XtB333oe+wooD5An9h4VdDAvfSIQkW8Cywhi64eBLwC3AreKyH7gOHBF2F5z\nDvBVVV1F0Df53rAF51TgDlXd6Z36MjKGWMCyWSYNWYQ4T6pZK4Xdeelx8fQyCnuaqOdddFTExGmZ\nSwQUQRHZLDNnzdNf+fVPZxr78APXWdVEw2iE8YbUtYJeRk89zUtvVVZLtwu2kR8T80lCloqLeVb1\nNdr5KCvT948w2t9XFUN3lGlhkWPXcH/V+2azW+IwATfSyNID1OhSmmkYXeQ5upGouBtGqzHP3Cit\nGM+551BVTnocZQy7tIrF2zdM6kqORjom5kYhTGTmS5TJIuiu8FeUbQ8vLU0bPKN9mJgbNSxbuZHR\n/mno8tGq7a4hsk9Si7uiYuo9Q8cZW1j/n+lkEfQ4Vl/wWFXHpJ6h46VsWm20FhNzo4plKzcyvGoq\nM/pGE8f4Qu8aOUBto+SivPO4SdA4yiboRSzzz8qK3iHoBS4IrjuwdrOJ+iTDxNyoIhDyZEF8aOd6\nBtZurnueosMsnSroaUQnSaNNsP3PkVQTJm7fit4h+GSw33nsvTtOmLB3OSbmRhVpQui62u+7Jajl\n7Y+NTszdf3iLCTrpxblkTw89Q8crxcjuP1w7ubls5UYAdvefV/XZZx6cgoSvd3Ne1TGdVLzMKA4T\ncyM3QWhlVpVoxmVaFCno3ShQYwtPMuceJ+TxGUW+N33RvDWVksKuOJlfxtUvV1xzvnKW5jEKxMTc\nqGLfxTdUxcGTmNF3tCYsEKWdxb3K5J1H2TXcHzuZXI92ZgwZ5cfE3GgIf7LTif+CLZssRS7Cit6h\nSmzbCfjc/SPAj6rGNdJTs6zrA4z2YGI+SXGx2Lha4jPS1+mEY47WePBl84TdL4d22uV74XP3v1h3\n/EXz1phIGw1hYj6JeePqnwA/qbyfUfD5i4yZZ538bDf+ZOe2h5dmFnEfE3SjEerWZhGRW0XkSFiX\n1207VUR2i8iz4XOPt+96EXlORIZE5MOtMtyoz0Xz1nDRvDUVL9wnS5qa82zrxcYdcUvNJ5so7Rru\nZ9dwP9seXkr/F19suIeq+9tZXNzIShbP/OvAXwF/4227DnhQVW8SkevC9+tF5ByCgurvA+YQNCo9\nW1XfLtZsIwt+SzaXGz54yzUMrN1cs7rTxxdvt6qQgz2MLTxZCVn4sWBH2WLmExlecd9HnpBKHJPt\n5mcUR10xV9XviMiZkc2XEHTWALgNeAhYH26/U1XfAg6JyHMEPe2+W4y5Rh5cGVldPloRm8XbNwTN\nqBJwQh63kjO4OfSk3ghaQc/Q8arYfhbaESeXPT30e/1N82JCbjRDozHz2ar6Svj6VYL2RxB0FHjU\nG3c43FaDiFwFXAXQ29vboBlGGk4cFm/fkLkk64y+o7H54j4Dazezm/M45nnq7thWVvab6Lh5vaX4\nfpbK3i+/nzn7Gxdyw2iWpidAw752uXvPqepWYCsEbeOatcNIZt/FtSs20xhYu5k59xyK9RSXrdwI\nXks32VMdfgGqij4B9Cyqvp/niSNP3z9SWSgDEyPoWZs1u8/Yf88h3lw0t8rOPJ/RPHKjCBoV89dE\n5AxVfUVEzgCOhNtHgPneuHnhNqPNPL/m2kyLgSBIUXxz0VyWrdxYM1Hqhztm9B1F+4CUCVJ/hWIz\n+KGWVgl6Fk/cCfipQ8crcfE3F9X++Hxz0dyaG1EUq5ViFEmjYn4fcAVwU/j8LW/7HSLyJYIJ0LOA\nx5o10iiOuInLOIZXTY1doegmUGcenBIIOfHx6ehEqMuoaSS7w4liI7HzInB1VARiY+KNfCYTcqNo\n6oq5iHyTYLLzNBE5DHyBQMTvFpErgReBjwGo6gERuRt4GjgBXG2ZLOXh2KFZ7CJb7HzmwSk1guPq\nnDtkT5CRGg2ztIJ6Xm6UvMv5k7xylys+vcF4eNRuJ/yWS24UTd08c1X9uKqeoarvUNV5qvo1Vf2h\nqn5QVc9S1eWq+oY3/kZVXaiq/ap6f2vNN/Lw/JprK9kqK3qHKg8ft39s4UkWbNlU2T6wdnNFyMcW\nnqyp4Bd3Dh93Y8gjyHHkCdtkzY9PI/rZGsH33F1s/c1Fcy2PvEsQkWtE5ICI7BeRb4rIOyP7/0BE\n9onIUyLyiIgMePvWhscdEJHPNGWHavvnHpcsWaJ79+5ttxmTBteAYvUF4xEwP/Tii2BWMXPeuX9s\nXM55s+LlVw3MGjfP4qEneea7bz2vqkxtoyTdxJzn3s6wiwuBRW3wqzTG0S6bReRxVV3SzDlmzpqn\nv/Lrn8409uEHrku8nojMBf4ZOEdV3wwjEztU9evemPOBZ1R1VEQuAv5UVX9NRBYBdxKkbx8HdgKf\nUtXnGvlMtpx/EjJ9/wi9zGUbS6sEPQ/R+LVbWIQnsHGLiNzEYLMEHvr4L4U0soRc4uqO7xru59jC\nk/SEm5uxPa6UrTvX9P0jVbn8STjhT7IhS9gmbjVw1L7KmIgtFudPZCowXUR+CrwLeNnfqaqPeG8f\nJUgMAfgl4Huq+r8BRORhYDXwxUaMMM98EuILh++hu2wNR5JXPnjLNVXvnacfxQloNO+8mclQZ7fD\nD/3Uo9G4vt9b09GI7WklbNPSGuuVvm009h7nkSd56WWiCM981rTZev7ptc2x49g58pcvAq97m7aG\nqdXOnrXAjcCbwC5V/YOkc4nIOuA/qOp/EZFfIkge+fXw2AeBvaqa7SdD9Nwm5pMT5x26cIUTOifm\nWYXc4fLYk8IVNw/cFWtDozQi6M1O0kZrkMdNyuYV+SzFyGyitJo2iHlamKUHuAe4lKCu8d8Bf6+q\n34gZeyHw18AHVPWH4bYrgf9GUPHuAPCWqjYUOzcxn+S4mi1xE5rOEx3tn5Yo4j4LtmxKDdu4FD//\nXJ0m6ECl3Rs0/uvCBLpxSibmvw+sVNUrw/d/CJynqv8tMm4xcC9wkar+a8K5/gdwWFX/OvsnGaf5\nqXqjYxlYu7mSmTLz4JTKA6gSqzkZ6408v+Za9q5/P9seXhq7f0XvELp8tGrx0v2Ht1QeeYkTUv8z\nxHHs0Kyms1x0+WglrNRsdo7R8QwD54nIu0REgA8Cz/gDRKQX2AZcHhVyETndG7MauKNRQ0zMJzGD\nt1zD82uurRG/uPS/rB70QzvXM/PgFLY9vLRmcZL/fmDt5sTJuEaI2lwvC6dZUZ/Rd5Q3rv4Jw6um\n8vJH+3KLuqUldgeq+j3g74EngKcINHWriHxKRD4VDtsA/Bzw1yLypIj4YYh7RORpYDvBupzqFlQ5\nsDCLAYyHWyBezKfvH8nlPTuv34+jx6U/uonFaLZHHqJNjn2KTF9Mw78xuLh6lhBM3l8kRU1+djJl\nCrOUCRNzowpf1B1Z4uVJ53LlcuMEHcYnFeOyJ7KKelzsPMpEibqPi63Xy06JMpEeeyfeDEzM4zEx\nN1qKHx/3M118UY9OikbJIm5ZBB2yizoU463PPDilMueQRTjLEHopu8CbmMdjMXOjpey7+IZKDRcn\n4NFwSz3PP+8EadqS/zzL85udKJ3Rd7RSgTJrTL0MQmrx/M7ExNwohKSeldH3fkMHJ5Z+DZg00kQ9\nGsooUtAbEXV3zIy+o7ELqspOGW4qRj4671+ZUUr8xS9VNT3CZ39JfZw4FtE/1Al6Fi84b010X5yz\njnWfeUbfUR7aWb/7UpI3HF0d6t5bjrvhY2LeRTgRbddSbCcSrlTuePGtqTVFuNx7f7vfdDrufV5a\nUf88+hmyjvd/fbgMnrwZQlB7o8pbL8aEvHuxCdAuwF+aH6VRIWwWPzXRx00KOruiIhdlbOHJisee\nJ46bdULUXaMd+PVekuq1+9Uh/e8ymtrpiP46iRP6Thd0mwCNxzzzLsD954xLK2zWu22GmQenMMas\nGq/cxwl1Utw8KfRSr/iUL45+WYJEO9sg6GMLT6LLfxKkZ0b6no+LeK2Q+++1D95YPr792KH31gh8\nlLibYqcLvGGeedcSJ+w+EyHuzjv3iXrfrg1dEnnsjBOpqLdbVi/dEf1+kuxJCvFEb5r++aIC36le\nu3nm8ZiYdzlxQukEotkJx6zXT+tKVI+8N50sgg7lFXX3Kybt10Jc6CqJaK2dJDpJ2E3M4zExnwTE\nVUZ0TISgw3gYJYuY1yuze+zQrIa6GCVluRS1yGgiiIubO+K+2zxt9jpF0E3M47GY+SShaFFyHnfW\nm0EzIu6ux8La7JA8N6Ok1MW0mHq74ulJxHnhvvfdaPZOEd2fjPZiYj4JSBOjBVs2VWWX5CXrgh8W\nFtMcOcrA2s3jYvzRPoC6JXujDZYdrUhlnAhcCeNoX9R633c9AXdNp8vonRu1mJgbQPxkJaQLQlav\ntVmv3N/v3zxm9B1ljFlAdeaHE58sqYxRbz3OSy+bdx6HX5Mekic7k4Q57rtyx5igdwYWM+9iMnvN\nE0zelnQ+rga6E9u4mHAR5XT9a/iUXdSj2Su+991Ioa8yirjFzOMxMe8S/Bi2q1SYpaZIXMw5KQMl\nScjiYrb+QpcsE3WONEHP0lk+jXaU1J1o4lIRW1FXvZ2YmMdjYt4luEyPKGmCnnXy0Img3wTCFzNf\nsNPEO088N07U0zoT5Slh0K6Suq0mLZulmwTdxDweE/Muwq8d7qiXgxzt8hNHVMwBhldNTS2clZZC\nWG8Jv08rPXVHmrindTHyKYOom5hnx8S8RZiYF0NUzJOE3A+L1Cv2lFT3ZfCWa2JvHnH49Vji4vhZ\ns1za5a1naUtXlknSLGGWThHtJEzM47Fsli4hq7DmXSTkC7nLmBhbeDLz9RxjC082PSHrpyFmHZ9E\ntDIjYVqjv9+/UdSbDC2bkENgc8/Q8Y4XbyMbJuYdTpKo+l55XgH365FHJzTTRMtfhp5WXCuKu0lk\noag88Hq1awbWboYW5Jz7n3MibgCj/dNYtnJj28oiGxOHiXkHUy+skncxkPNEhz/3Xu8c4+fyxccJ\nN8R7qHm78/jhim4gy+fIMiar4KedywR9ctCy/zkislJEhkTkORG5rlXXmawULeQQ/KePtjgbW3iy\n8vDP3bvjRE3Io1EhdotdJqpOTCvxF+4UQe+OE5muWY/R/ml1f40YjVFP6yRgS7h/n4i8P+uxeWiJ\nZy4ipwBfBj4EHAa+LyL3qerTrbhet7Fs5caKqMb9LE/LIsmaoRJHUrnaOO9wfFJtLjCt5hj/uDSx\nid4oBm+5JrXZRj2K8GQbuUYrflG4m2XvjhM16aBGOciodRcBZ4WPXwP+X+DXitbJVoVZlgLPqerz\nACJyJ3AJYGJeBxfqcB5ZdFk5AAeDbvfPF1yT3A+l+J69v4jI3Sz8SbWox5ck4ml55jMPTmFg7WYG\nb7mm0k90+n4aFvWJotVhoaQORHnJM3Fs5CKL1l0C/I0GqYOPisjPisgZwJkZjs1Mq8R8LvCS9/4w\nwR2pgohcBVwVvn1LRPa3yJaJ4DTg9Zac+YHkXbLls0VcIdH2tPO3KnJWc80R0r6D1n3vE0M2+9O/\ng1wU+Hdr53f/3vpD0vnxT488sHPkL0/LOPydIuLnTm9V1a3h67palzBmbsZjM9O2CdDwy9gKICJ7\nOyGPM4lOtt9sbx+dbH8n2w6gqivbbUPRtErMR4D53vt54TbDMIxuIovWJY15R4ZjM9OqgN/3gbNE\npE9EpgGXAfe16FqGYRjtIovW3Qf8YZjVch5wVFVfyXhsZlrimavqCRH5I4JI3ynArap6IOWQrSn7\nOoFOtt9sbx+dbH8n214YSVonIp8K9/9PYAewCngO+N/A/5V2bKO2lKI2i2EYhtEc3bHczjAMY5Jj\nYm4YhtEFtF3MO23Zv4i8ICJPiciTLvdURE4Vkd0i8mz43NNuOx0icquIHPHz+NPsFZHrw7/FkIh8\nuD1WV2yJs/1PRWQk/P6fFJFV3r4y2T5fRP5JRJ4WkQMisjbcXvrvPsX2jvjuJy2q2rYHQdD/ILCA\nYE34IHBOO23KYPMLwGmRbV8ErgtfXwdsbLednm3/EXg/sL+evcA54d/gZ4C+8G9zSsls/1NgXczY\nstl+BvD+8PVM4F9DG0v/3afY3hHf/WR9tNszryyFVdXjgFvO2mlcAtwWvr4N+N022lKFqn4HeCOy\nOcneS4A7VfUtVT1EMPu+dEIMjSHB9iTKZvsrqvpE+HoMeIZgxV/pv/sU25Moje2TmXaLedIy1zKj\nwB4ReTwsSQAwW4O8UYBXgdntMS0zSfZ2yt/j02H1uVu9MEVpbReRM4FfBr5Hh333Eduhw777yUS7\nxbwT+YCqnktQCe1qEfmP/k4Nfnd2TL5np9lLUHFuAXAu8ArQXPuiFiMiM4B7gM+o6o/9fWX/7mNs\n76jvfrLRbjHvuGX/qjoSPh8B7iX4OflaWAWN8PlI+yzMRJK9pf97qOprqvq2qp4E/j/Gf86XznYR\neQeBGP6tqm4LN3fEdx9neyd995ORdot5Ry37F5F3i8hM9xpYAewnsPmKcNgVwLfaY2Fmkuy9D7hM\nRH5GRPoI6i8/1gb7EnFCGPJ7BN8/lMx2ERHga8Azqvolb1fpv/sk2zvlu5+0tHsGlmCZ678SzIB/\nvt321LF1AcGs/SBwwNkL/BzwIPAssAc4td22ejZ/k+An8U8JYplXptkLfD78WwwBF5XQ9tuBp4B9\nBCJyRklt/wBBCGUf8GT4WNUJ332K7R3x3U/Why3nNwzD6ALaHWYxDMMwCsDE3DAMowswMTcMw+gC\nTMwNwzC6ABNzwzCMLsDE3DAMowswMTcMw+gC/n+ZDJVk9XqirAAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig, ax = plt.subplots()\n", "diff = ax.contourf(ma.masked_array(list_of_depths[40,:,:], mask = 1-new_domain[0,:,:]),\n", " levels = np.linspace(0,58,8))\n", "fig.colorbar(diff,ax=ax)" ] }, { "cell_type": "code", "execution_count": 38, "metadata": { "collapsed": true }, "outputs": [], "source": [ "JanW04 = nc.Dataset('/ocean/vdo/MEOPAR/completed-runs/SalishSeaLake/LakeJan-0.4/SalishSea_1h_20170101_20170107_grid_T.nc')" ] }, { "cell_type": "code", "execution_count": 23, "metadata": { "collapsed": false }, "outputs": [ { "ename": "RuntimeError", "evalue": "No mappable was found to use for colorbar creation. First define a mappable such as an image (with imshow) or a contour set (with contourf).", "output_type": "error", "traceback": [ "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m", "\u001b[0;31mRuntimeError\u001b[0m Traceback (most recent call last)", "\u001b[0;32m\u001b[0m in \u001b[0;36m\u001b[0;34m()\u001b[0m\n\u001b[1;32m 9\u001b[0m \u001b[0minterval\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0;36m0.25\u001b[0m\u001b[0;31m#in seconds\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 10\u001b[0m \u001b[0mani40\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0manimation\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mFuncAnimation\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfig\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0manimate30\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0mframes\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;36m10\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0minterval\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0minterval\u001b[0m\u001b[0;34m*\u001b[0m\u001b[0;36m1e+3\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mrepeat\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;32mFalse\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m---> 11\u001b[0;31m \u001b[0mplt\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mcolorbar\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 12\u001b[0m \u001b[0mani40\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n", "\u001b[0;32m/home/vdo/anaconda3/lib/python3.6/site-packages/matplotlib/pyplot.py\u001b[0m in \u001b[0;36mcolorbar\u001b[0;34m(mappable, cax, ax, **kw)\u001b[0m\n\u001b[1;32m 2252\u001b[0m \u001b[0mmappable\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mgci\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 2253\u001b[0m \u001b[0;32mif\u001b[0m \u001b[0mmappable\u001b[0m \u001b[0;32mis\u001b[0m \u001b[0;32mNone\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m-> 2254\u001b[0;31m raise RuntimeError('No mappable was found to use for colorbar '\n\u001b[0m\u001b[1;32m 2255\u001b[0m \u001b[0;34m'creation. First define a mappable such as '\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 2256\u001b[0m \u001b[0;34m'an image (with imshow) or a contour set ('\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n", "\u001b[0;31mRuntimeError\u001b[0m: No mappable was found to use for colorbar creation. First define a mappable such as an image (with imshow) or a contour set (with contourf)." ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAeQAAAHiCAYAAAA597/kAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAFqxJREFUeJzt3X/M7ndd3/HXuz3t6KBKEawHRBCDJuDkQE46omiYgKlk\nCriFwRgrk1nMlEDGtjDJBLNsk40fiUSZJTRUgygbIJWgW2lYiJvW1a4trWUUSZmU054BSqkQsPS9\nP+7vcbc3933O3XPu+5z3ua7HI7lyX9f38/3e1+fb733neV/X9e33VHcHADizzjnTEwAABBkARhBk\nABhAkAFgAEEGgAEEGQAGEGTYZ1V1R1U980zPA5hNkIGvUxteX1WfW26vr6raNP7YqvpwVX2pqj62\n9Q+Oqvr7VfWpqvrzqvrNqnrY6d8LOLsIMqyIqjqwh9/u8iTPTfKkJN+T5EeSvGzT+LuS/K8k35Tk\nNUn+c1U9YpnHE5P8cpIXJ7k4yZeS/NIezg1WkiDD6XGoqm6uqi9U1W9U1YOODVTVT1TVJ6rq81V1\ndVU9cln+2KrqzaGtqv9WVf94uf+SqvrvVfXmqvpcktft4XwvS/LG7v50d9+Z5A1JXrI873cmeUqS\n13b3l7v7PUluTvJ3lm1flOS3uvsj3X1vkn+V5Meq6sI9nB+sHEGG0+P5SS5N8u3ZeMX5kiSpqh9M\n8u+W8YNJPpXk1x/A9/2bST6ZjVei/2br4PLW8Z8d5/ZtO3zfJya5adPjm5Zlx8Y+2d1fPM74X27b\n3X+c5CtJvvMB7Besnb18iwvY2S9092eSpKp+K8mhZfmLklzZ3TcsY/8yyZ9W1WN3+X0/091vWe7f\nt3Wwu38tya+dxHwfkuQLmx7fk+Qhy+fIW8eOjT9qh22PjXuFDMfhFTKcHndtuv+lbEQrSR6ZjVfF\nSZLlLd7P5f/H7UT+ZE9m9/XuTfINmx5/Y5J7e+Nfo9k6dmz8iztsu3Uc2IYgw5n1mSSPOfagqh6c\njROl7kzy58viv75p/W/Zsv1x/7m2qnpRVd17nNtOb1nfmo0Tuo550rLs2NjjtnwmvHX8L7etqu9I\ncn6Sjx9vrrDuBBnOrHcl+UdVdaiq/lqSf5vkuu6+o7v/bzbC/A+q6tyq+vEk3/FAvnl3v7O7H3Kc\n2//ZYdNfSfJPq+pRVfWoJK9K8o7le348yY1JXltVD6qqH0vyN5K8Z9n2nUl+pKq+f/kD418nee+W\nz5yBLQQZzqDu/lA2zkJ+T5Ij2QjuCzat8hNJ/nk23sZ+YpL/cZqm9stJfivJR5fbB5Zlx7wgyeEk\nf5qNk9L+7vIHRLr71iQ/mY0wH03y4CT/5DTNG85atfGREABwJnmFDAADCDIADCDIADCAIAPAAIIM\nAAOc1ktnnn/OBX3Bga0X8AGA1XTPXxz9bHc/YjfrntYgX3DgG/K93/z3TudTAsAZ8zt3vuVTJ15r\ng7esAWAAQQaAAQQZAAYQZAAYQJABYABBBoABBBkABhBkABhAkAFgAEEGgAEEGQAGEGQAGECQAWAA\nQQaAAQQZAAYQZAAYQJABYABBBoABBBkABhBkABhAkAFgAEEGgAEEGQAGEGQAGECQAWAAQQaAAQQZ\nAAYQZAAYQJABYABBBoABBBkABhBkABhAkAFgAEEGgAEEGQAGEGQAGECQAWAAQQaAAQQZAAY4YZCr\n6kFV9QdVdVNV3VpVP7csf11V3VlVNy63Z+//dAFgNR3YxTpfSfKD3X1vVZ2X5Her6reXsTd39xv2\nb3oAsB5OGOTu7iT3Lg/PW269n5MCgHWzq8+Qq+rcqroxydEk13T3dcvQy6vq5qq6sqou2rdZAsCK\n21WQu/tr3X0oybcmuaSqvjvJW5M8LsmhJEeSvHG7bavq8qq6vqqu/+r9X96jaQPAanlAZ1l3958l\n+XCSS7v77iXU9yd5W5JLdtjmiu4+3N2Hzz/nglOfMQCsoN2cZf2Iqnrocv+CJM9K8rGqOrhptecl\nuWV/pggAq283Z1kfTHJVVZ2bjYC/u7s/UFW/WlWHsnGC1x1JXrZ/0wSA1babs6xvTvLkbZa/eF9m\nBABryJW6AGAAQQaAAQQZAAYQZAAYQJABYABBBoABBBkABhBkABhAkAFgAEEGgAEEGQAGEGQAGECQ\nAWAAQQaAAQQZAAYQZAAYQJABYABBBoABBBkABhBkABhAkAFgAEEGgAEEGQAGEGQAGECQAWAAQQaA\nAQQZAAYQZAAYQJABYABBBoABBBkABhBkABhAkAFgAEEGgAEEGQAGEGQAGECQAWAAQQaAAQQZAAYQ\nZAAYQJABYABBBoABBBkABhBkABhAkAFgAEEGgAEEGQAGOGGQq+pBVfUHVXVTVd1aVT+3LH9YVV1T\nVbcvXy/a/+kCwGrazSvkryT5we5+UpJDSS6tqqcmeXWSa7v78UmuXR4DACfhhEHuDfcuD89bbp3k\nOUmuWpZfleS5+zJDAFgDu/oMuarOraobkxxNck13X5fk4u4+sqxyV5KLd9j28qq6vqqu/+r9X96T\nSQPAqtlVkLv7a919KMm3Jrmkqr57y3hn41Xzdtte0d2Hu/vw+edccMoTBoBV9IDOsu7uP0vy4SSX\nJrm7qg4myfL16N5PDwDWw27Osn5EVT10uX9Bkmcl+ViSq5Nctqx2WZL379ckAWDVHdjFOgeTXFVV\n52Yj4O/u7g9U1e8leXdVvTTJp5I8fx/nCQAr7YRB7u6bkzx5m+WfS/KM/ZgUAKwbV+oCgAEEGQAG\nEGQAGECQAWAAQQaAAQQZAAYQZAAYQJABYABBBoABBBkABhBkABhAkAFgAEEGgAEEGQAGEGQAGECQ\nAWAAQQaAAQQZAAYQZAAYQJABYABBBoABBBkABhBkABhAkAFgAEEGgAEEGQAGEGQAGECQAWAAQQaA\nAQQZAAYQZAAYQJABYABBBoABBBkABhBkABhAkAFgAEEGgAEEGQAGEGQAGECQAWAAQQaAAQQZAAYQ\nZAAYQJABYABBBoABBBkABjhhkKvq0VX14ar6o6q6tapesSx/XVXdWVU3Lrdn7/90AWA1HdjFOvcl\neVV331BVFyb5w6q6Zhl7c3e/Yf+mBwDr4YRB7u4jSY4s979YVbcledR+TwwA1skD+gy5qh6b5MlJ\nrlsWvbyqbq6qK6vqoj2eGwCsjV0HuaoekuQ9SV7Z3fckeWuSxyU5lI1X0G/cYbvLq+r6qrr+q/d/\neQ+mDACrZ1dBrqrzshHjd3b3e5Oku+/u7q919/1J3pbkku227e4ruvtwdx8+/5wL9mreALBSdnOW\ndSV5e5LbuvtNm5Yf3LTa85LcsvfTA4D1sJuzrL8vyYuTfLSqblyW/UySF1bVoSSd5I4kL9uXGQLA\nGtjNWda/m6S2Gfrg3k8HANaTK3UBwACCDAADCDIADCDIADCAIAPAAIIMAAMIMgAMIMgAMIAgA8AA\nggwAAwgyAAwgyAAwgCADwACCDAADCDIADCDIADCAIAPAAIIMAAMIMgAMIMgAMIAgA8AAggwAAwgy\nAAwgyAAwgCADwACCDAADCDIADCDIADCAIAPAAIIMAAMIMgAMIMgAMIAgA8AAggwAAwgyAAwgyAAw\ngCADwACCDAADCDIADCDIADCAIAPAAIIMAAMIMgAMIMgAMIAgA8AAggwAAwgyAAxwwiBX1aOr6sNV\n9UdVdWtVvWJZ/rCquqaqbl++XrT/0wWA1bSbV8j3JXlVdz8hyVOT/FRVPSHJq5Nc292PT3Lt8hgA\nOAknDHJ3H+nuG5b7X0xyW5JHJXlOkquW1a5K8tz9miQArLoH9BlyVT02yZOTXJfk4u4+sgzdleTi\nPZ0ZAKyRXQe5qh6S5D1JXtnd92we6+5O0jtsd3lVXV9V13/1/i+f0mQBYFXtKshVdV42YvzO7n7v\nsvjuqjq4jB9McnS7bbv7iu4+3N2Hzz/ngr2YMwCsnN2cZV1J3p7ktu5+06ahq5Nctty/LMn79356\nALAeDuxine9L8uIkH62qG5dlP5Pk55O8u6pemuRTSZ6/P1MEgNV3wiB39+8mqR2Gn7G30wGA9eRK\nXQAwgCADwACCDAADCDIADCDIADCAIAPAAIIMAAMIMgAMIMgAMIAgA8AAggwAAwgyAAwgyAAwgCAD\nwACCDAADCDIADCDIADCAIAPAAIIMAAMIMgAMIMgAMIAgA8AAggwAAwgyAAwgyAAwgCADwACCDAAD\nCDIADCDIADCAIAPAAIIMAAMIMgAMIMgAMIAgA8AAggwAAwgyAAwgyAAwgCADwACCDAADCDIADCDI\nADCAIAPAAIIMAAMIMgAMIMgAMIAgA8AAggwAA5wwyFV1ZVUdrapbNi17XVXdWVU3Lrdn7+80AWC1\n7eYV8juSXLrN8jd396Hl9sG9nRYArJcTBrm7P5Lk86dhLgCwtk7lM+SXV9XNy1vaF+3ZjABgDZ1s\nkN+a5HFJDiU5kuSNO61YVZdX1fVVdf1X7//yST4dAKy2kwpyd9/d3V/r7vuTvC3JJcdZ94ruPtzd\nh88/54KTnScArLSTCnJVHdz08HlJbtlpXQDgxA6caIWqeleSpyd5eFV9Oslrkzy9qg4l6SR3JHnZ\nPs4RAFbeCYPc3S/cZvHb92EuALC2XKkLAAYQZAAYQJABYABBBoABBBkABhBkABhAkAFgAEEGgAEE\nGQAGEGQAGECQAWAAQQaAAQQZAAYQZAAYQJABYABBBoABBBkABhBkABhAkAFgAEEGgAEEGQAGEGQA\nGECQAWAAQQaAAQQZAAYQZAAYQJABYABBBoABBBkABhBkABhAkAFgAEEGgAEEGQAGEGQAGECQAWAA\nQQaAAQQZAAYQZAAYQJABYABBBoABBBkABhBkABhAkAFgAEEGgAEEGQAGEGQAGECQAWCAEwa5qq6s\nqqNVdcumZQ+rqmuq6vbl60X7O00AWG27eYX8jiSXbln26iTXdvfjk1y7PAYATtIJg9zdH0ny+S2L\nn5PkquX+VUmeu8fzAoC1cuAkt7u4u48s9+9KcvFOK1bV5UkuT5IHnXvhST4dAKy2Uz6pq7s7SR9n\n/IruPtzdh88/54JTfToAWEknG+S7q+pgkixfj+7dlABg/ZxskK9Octly/7Ik79+b6QDAetrN//b0\nriS/l+S7qurTVfXSJD+f5FlVdXuSZy6PAYCTdMKTurr7hTsMPWOP5wIAa8uVugBgAEEGgAEEGQAG\nEGQAGECQAWAAQQaAAQQZAAYQZAAYQJABYABBBoABBBkABhBkABhAkAFgAEEGgAEEGQAGEGQAGECQ\nAWAAQQaAAQQZAAYQZAAYQJABYABBBoABBBkABhBkABhAkAFgAEEGgAEEGQAGEGQAGECQAWAAQQaA\nAQQZAAYQZAAYQJABYABBBoABBBkABhBkABhAkAFgAEEGgAEEGQAGEGQAGECQAWAAQQaAAQQZAAYQ\nZAAYQJABYABBBoABDpzKxlV1R5IvJvlakvu6+/BeTAoA1s0pBXnxt7r7s3vwfQBgbXnLGgAGONUg\nd5IPVdUfVtXlezEhAFhHp/qW9dO6+86q+uYk11TVx7r7I5tXWEJ9eZI86NwLT/HpAGA1ndIr5O6+\nc/l6NMn7klyyzTpXdPfh7j58/jkXnMrTAcDKOukgV9WDq+rCY/eT/FCSW/ZqYgCwTk7lLeuLk7yv\nqo59n1/r7t/Zk1kBwJo56SB39yeTPGkP5wIAa8v/9gQAAwgyAAwgyAAwgCADwACCDAADCDIADCDI\nADCAIAPAAIIMAAMIMgAMIMgAMIAgA8AAggwAAwgyAAwgyAAwgCADwACCDAADCDIADCDIADCAIAPA\nAIIMAAMIMgAMIMgAMIAgA8AAggwAAwgyAAwgyAAwgCADwACCDAADCDIADCDIADCAIAPAAIIMAAMI\nMgAMIMgAMIAgA8AAggwAAwgyAAwgyAAwgCADwACCDAADCDIADCDIADCAIAPAAIIMAAMIMgAMIMgA\nMMApBbmqLq2q/11Vn6iqV+/VpABg3Zx0kKvq3CS/mOSHkzwhyQur6gl7NTEAWCen8gr5kiSf6O5P\ndvdXk/x6kufszbQAYL2cSpAfleRPNj3+9LIMAHiADuz3E1TV5UkuXx5+5XfufMst+/2cQzw8yWfP\n9CROg3XZz8S+rqp12dd12c9k1r4+ZrcrnkqQ70zy6E2Pv3VZ9ld09xVJrkiSqrq+uw+fwnOeNdZl\nX9dlPxP7uqrWZV/XZT+Ts3dfT+Ut6/+Z5PFV9e1VdX6SFyS5em+mBQDr5aRfIXf3fVX100n+S5Jz\nk1zZ3bfu2cwAYI2c0mfI3f3BJB98AJtccSrPd5ZZl31dl/1M7OuqWpd9XZf9TM7Sfa3uPtNzAIC1\n59KZADDAvgT5RJfUrA2/sIzfXFVP2Y957KeqenRVfbiq/qiqbq2qV2yzztOr6gtVdeNy+9kzMde9\nUFV3VNVHl/24fpvxs/6YJklVfdem43VjVd1TVa/css5Ze1yr6sqqOlpVt2xa9rCquqaqbl++XrTD\ntmfVpXJ32Nf/UFUfW35G31dVD91h2+P+vE+yw36+rqru3PQz+uwdtl2FY/obm/bzjqq6cYdt5x/T\n7t7TWzZO8PrjJI9Lcn6Sm5I8Ycs6z07y20kqyVOTXLfX89jvW5KDSZ6y3L8wyce32c+nJ/nAmZ7r\nHu3vHUkefpzxs/6YbrNP5ya5K8ljVuW4JvmBJE9JcsumZf8+yauX+69O8vod/lsc9/d62m2Hff2h\nJAeW+6/fbl+XseP+vE+67bCfr0vyz06w3Uoc0y3jb0zys2frMd2PV8i7uaTmc5L8Sm/4/SQPraqD\n+zCXfdPdR7r7huX+F5PclvW+UtlZf0y38Ywkf9zdnzrTE9kr3f2RJJ/fsvg5Sa5a7l+V5LnbbHrW\nXSp3u33t7v/a3fctD38/G9dPOKvtcEx3YyWO6TFVVUmen+Rdp3VSe2g/grybS2qu1GU3q+qxSZ6c\n5Lpthr93eXvst6vqiad1Ynurk3yoqv5wufraVit1TBcvyM6/3KtyXJPk4u4+sty/K8nF26yzisf3\nx7Pxrs52TvTzfjZ4+fIzeuUOH0Os2jH9/iR3d/ftO4yPP6ZO6jpFVfWQJO9J8sruvmfL8A1Jvq27\nvyfJW5L85ume3x56Wncfysa/7vVTVfUDZ3pC+2m52M2PJvlP2wyv0nH9K3rjvb2V/18vquo1Se5L\n8s4dVjnbf97fmo23og8lOZKNt3JX3Qtz/FfH44/pfgR5N5fU3NVlN6erqvOyEeN3dvd7t4539z3d\nfe9y/4NJzquqh5/mae6J7r5z+Xo0yfuy8XbXZitxTDf54SQ3dPfdWwdW6bgu7j728cLy9eg266zM\n8a2qlyT520letPwB8nV28fM+Wnff3d1f6+77k7wt289/lY7pgSQ/luQ3dlrnbDim+xHk3VxS8+ok\n/3A5M/epSb6w6S2zs8LyecXbk9zW3W/aYZ1vWdZLVV2Sjf/enzt9s9wbVfXgqrrw2P1snBiz9R8J\nOeuP6RY7/rW9Ksd1k6uTXLbcvyzJ+7dZZyUulVtVlyb5F0l+tLu/tMM6u/l5H23L+RvPy/bzX4lj\nunhmko9196e3Gzxrjul+nCmWjTNuP56NM/hesyz7ySQ/udyvJL+4jH80yeEzfXbbSezj07Lx1t7N\nSW5cbs/esp8/neTWbJy9+PtJvvdMz/sk9/Vxyz7ctOzPSh7TTfv74GwE9hs3LVuJ45qNPzKOJPmL\nbHxm+NIk35Tk2iS3J/lQkoct6z4yyQc3bft1v9eTbzvs6yey8bnpsd/Z/7h1X3f6eZ9622E/f3X5\nPbw5G5E9uKrHdFn+jmO/n5vWPeuOqSt1AcAATuoCgAEEGQAGEGQAGECQAWAAQQaAAQQZAAYQZAAY\nQJABYID/B9fYaggai3nhAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig,ax = plt.subplots(figsize=(8,8))\n", "def animate30(i):\n", " ax.clear()\n", " den = ax.contourf(np.ones((40,20)), \n", " levels = np.linspace(0, 58, 8) #, colors='black'\n", " )\n", " ax.set_title('hour = %03d'%(i)) \n", " return ax\n", "interval = 0.25#in seconds\n", "ani40 = animation.FuncAnimation(fig,animate30,frames=10,interval=interval*1e+3, repeat=False)\n", "plt.colorbar()\n", "ani40" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": true }, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": 46, "metadata": { "collapsed": false }, "outputs": [], "source": [ "Jun04 = nc.Dataset('/ocean/vdo/MEOPAR/completed-runs/SalishSeaLake/LakeJun0.4/SalishSea_1h_20160601_20160607_grid_T.nc')" ] }, { "cell_type": "code", "execution_count": 48, "metadata": { "collapsed": false }, "outputs": [], "source": [ "def animate_one_density(file, rho):\n", " depth3 = file.variables['deptht'][:]\n", " saline3 = file.variables['vosaline'][:]\n", " temp3 = file.variables['votemper'][:]\n", " a = np.expand_dims(depth, 1).T\n", " b = np.expand_dims(np.expand_dims(a,2), 3)\n", " c = np.ones((144,40,564, 284)) * b\n", " density3 = calc_rho(saline3, temp3, c)\n", " index3 = (np.abs(density3 - rho)).argmin(axis=1)\n", " list_of_depths3 = np.zeros((144,564,284))\n", " for t in range(144):\n", " for y in range(564):\n", " for x in range(284):\n", " list_of_depths3[t,y,x] = depth3[index3[t,y,x]]\n", " mask_land = 1 - new_domain[0,:,:]\n", " fig,ax = plt.subplots(figsize=(8,8))\n", " def animate30(i):\n", " ax.clear()\n", " masked_depths = ma.masked_array(list_of_depths3[i,:,:], mask = mask_land)\n", " den = ax.contourf(masked_depths, \n", " levels = np.linspace(0, 58, 8) #, colors='black'\n", " )\n", " #ax.set_ylim(ylims)\n", " #ax.set_xlim(xlims)\n", " ax.set_title('hour = %03d'%(i)) \n", " return ax\n", "\n", " interval = 0.25#in seconds\n", " ani40 = animation.FuncAnimation(fig,animate30,frames=144,interval=interval*1e+3, repeat=False)\n", " return ani40" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# June, vtau = 0.4" ] }, { "cell_type": "code", "execution_count": 49, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/html": [ "" ], "text/plain": [ "" ] }, "execution_count": 49, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAeoAAAHiCAYAAAAnPo9XAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3X+wHedd3/HPNwElahLgukmFLdkooY5mYrUy7Z3QmfSH\nhrQgUagpZTymkHGKwXQakCJouTZMm8x0XCx+RJULdHBLBjPEcTyEEMNgpTjDnUzb/LpOIyI7veAm\ndpGQfyQWkBCPhd2nf5zdq72r3bO/nt19nt33a8Zj6dxzz9mzZ7Wf/T77/DDnnAAAQJheMvYGAACA\ncgQ1AAABI6gBAAgYQQ0AQMAIagAAAkZQAwAQMIIaGImZPW5m/3Ds7QAQNoIaQG22cNzMvpj8d9zM\nLPPzvWb2+2b2FTP73/kLETP752b2hJn9hZn9lpldMfynAOJCUAMTZ2Zf5fHlbpX0XZIOSPqbkr5T\n0g9nfv5eSf9L0l+V9FOSfsPMXpNsx3WSflnSWyTtkvQVSb/kcduASSKogXFdb2Z/YGZ/ZmbvM7OX\npz8wsx8ys8fM7Fkze8DMrkoe32tmLhvAZrZuZj+Y/PmtZvY/zOyEmX1R0js9bu/Nkn7eOXfWOXdO\n0s9Jemvyvq+X9LckvcM595xz7v2S/kDSP0t+9/sk/bZz7iPOuS9L+reSvtvMXuVx+4DJIaiBcd0o\n6ZCk12pRob5VkszsWyT9dPLzKyU9Iem+Bq/7zZI+p0Xlekf+h0kT9J8u+e+akte9TtLpzN9PJ4+l\nP/ucc+5LS36+9bvOuf8j6XlJr2/wuYDZ8dkkBqC5u5xzfyJJZvbbkq5PHv8+Se92zn0q+dntki6Y\n2d6ar/snzrn/lPz5hfwPnXP3Srq3xfa+UtKfZf7+55Jemdynzv8s/fnukt9Nf05FDSxBRQ2M68nM\nn7+iRZhJ0lVaVNGSpKSp+Iu6FHpV/tjL1l3uy5K+JvP3r5X0ZbdY3Sf/s/TnXyr53fzPARQgqIEw\n/Ymkb0j/Ymav0KKD1jlJf5E8/Fcyz//63O8vXRbPzL7PzL685L+ypu9HtOhIljqQPJb+7HW5e875\nn2/9rpl9o6Qdkv5w2bYCc0dQA2F6r6R/YWbXm9nLJP0HSR93zj3unHtGi8D+fjN7qZn9gKRvbPLi\nzrn3OOdeueS//1vyq78m6cfMbLeZ7Zb045J+NXnNP5T0aUnvMLOXm9l3S/obkt6f/O57JH2nmf29\n5MLj30v6zdw9bQA5BDUQIOfcQ1r0in6/pPNaBPFNmaf8kKR/o0Vz+HWS/udAm/bLkn5b0meS/34n\neSx1k6RVSRe06Az3PcmFhZxzj0j6l1oE9tOSXiHpXw203UC0bHFrCQAAhIiKGgCAgBHUAAAEjKAG\nACBgBDUAAAEjqAEACFgQU4i++tWvdnv37h17MwBM3OYfXZoIbt+1+TligGE9/PDDX3DOvabqeUEE\n9d69e7WxsTH2ZgCYsANHT+h1b7r097+UdPrksdG2BzCzJ6qfRdM3gBk4cPREo8eBkBDUACaNMEbs\ngmj6BgDf6gZ0/nk0hyM0BDWAyfBRPR88dFzrp9Y8bA3gB0ENIHg0X2POCGoAQRornKmmERo6kwFA\ngvvTCBEVNYBZKArhbNVOSCNUBDWAVsqapusEXh/N2m2ClnBGDAhqAKXaBOqQ95YJWswBQQ3gMvSy\nBsJBUAOIDpU05oSgBhA0QhlzR1ADCBIBDSwQ1AAuc/rksd7uUxPAQDMENYBCaaDWDWwCGOgHM5MB\nKEVIA+MjqAEUYogWEAaavgFsQ0ADYaGiBrCFkAbCQ1AD6IyAB/pjzrmxt0Grq6tuY2Nj7M0AkNEl\nfOlcBlQzs4edc6uVzyOoASxDYAP9qBvUNH0DWOr0yWMELjAighpALYQ1MA6CGgCAgBHUAGprWlXT\nGxzojqAG0AhN4MCwmJkMQCGqYSAMBDUASQQzECqCGpg5AhoIG0ENzNQQAc39bKA7gjpS+ZMsJ0SE\nhOMR8IegBmaor2qagAb8I6iBiSOUgbgR1AAqEcrAeJjwJELcn0ZdPqppji9gXFTUESk66XISRRlC\nGpgGgjoCBDQAzBdBHQgmnUBIuBAEwkFQj4yARl9OnzzW+PgioIHw0JkMmChCGpgGc86NvQ161dfu\nca9764/Xeu7cTibLTrZz2xeop00rDccSMDwze9g5t1r1vFpN32b2uKQvSXpR0gvOuVUzu0LS+yTt\nlfS4pBudcxeS598u6Zbk+Ueccx9q8RkKVZ2EpnbCST/PkE3k2fea2v6csibHCN8rEI9aFXUS1KvO\nuS9kHvsZSc865+40s9skrTjn1szsDZLeK+mNkq6S9JCk1zvnXix7/SYVdR1TPQn56v3NCX1a+D6B\nOHmtqEvcIOlg8ud7JK1LWksev88597ykz5vZY1qE9kfLXmjftV+vjQaVY9HJZg6dsrKfO/286f+X\nnYDnsG/moMv3SEAD8aob1E7SQ2b2oqRfds7dLWmXc+588vMnJe1K/rxb0scyv3s2eaxXczsR5Xv0\nHjh64rJ9QGei8LS9f8zFFjBfdZu+dzvnzpnZX5P0e5J+VNIDzrmvyzzngnNuxcx+QdLHnHO/njz+\nK5IedM79Ru41b5V0qyRdc801f/uJJ57Y9p5zuxfdlM8T99z3pQ8hBynfLxAmr03fzrlzyf+fNrMP\naNGU/ZSZXemcO29mV0p6Onn6OUlXZ359T/JY/jXvlnS3JK2url52tVBVRRRVkHNAQIcn1JDm+wWm\noXIctZm9wsxelf5Z0rdKOiPpAUk3J0+7WdIHkz8/IOkmM3uZmb1W0rWSPtFm406fPFZ57zXUk2Qf\nfH/WOe27voS6DwlpYDrqVNS7JH3AzNLn3+ucO2Vmn5R0v5ndIukJSTdKknPuETO7X9Kjkl6Q9LZl\nPb7rqFNd558/NV0+Y6hhAv+meOwDc1cZ1M65z0k6UPD4FyW9ueR37pB0R+etyxhjPPGYhvqcdXqN\nh6Dr/uh6YVPU4z4koX9/ANpjru+R1T3pdzkR1wmZWAK7ra7h6iucfQ8tnOr3BeCS6IK6TmUdysxa\ny5qrhwjoKQixem1j2fc49+8YwHJBzPW9urrqNjY2Or1GaHNiD9lU20Ro+2mZOYQ0gPkaYmayoCyr\ntIcKzZgnGAlpW6YS0ADgw2SCuk9lwdFlxqixgjGkQC4yhZAOfR8DiMvkgrrrNJpNTCFUpsZXSMZ0\niwDAtE0uqPNCG1bDSb4/Pvct3xOAUEw+qFNjLAUYwoVBTGK7jQAAQ5hNUA/N19rRUxdq73gACMUs\ngpoxy+OgRQEAupt8UNcJiz4COj9cjIuA7tiHAOZo0kEdQkU353DxNT/7nPchAEw6qOsgBPrj40KJ\n7wfA3M0+qBEmAhoAFiYZ1CE0eaMdAhoAtnvJ2BswtgNHTxDsPWkSuqdPHiOkAaDAJCvqNg4cPUFQ\nZPjqrT7klK4AMEWTDOomvY0J5+V8XsAQ0gDQ3CSDOlW1uhUhXU+X6pqLJQDoZtJBXRUS+Z8TGMv1\nURGzzwFguUkHNcJFQANAPZMNau6HhodwBoDmJhnUVSFddu+ant/tsM8AoD+TDOo6yoYNsYhGM+wn\nAOiXOefG3gatrq66jY2Nzq9Tp5Ju+hpzDaKmtw7mup8AoC0ze9g5t1r1vMnMTOYjpHEJ+wsAwjD5\npu+mgUNAXcI4dAAY3+SDGu302Wu+6LUJfQAoNomgZijWOHx2vMt+h4Q2AFwymXvUGE/dFciaPA8A\nsBB9Rc1JvR9V96eL+PouqKgB4JJog5qAngZCGQCWizKo64Y0IdDesvHkPi6S+G4AoJ4og7oupgT1\no0lIs78BwK8ogzoNg6KqjyZxP9oELiENAP5FGdRSedMsYQEAmJLohmfVHeKDfiy7EOK7AQD/ogtq\nDC8fwLRaAMBwogrqsqknCY5h1K2WqawBwJ9ogrronjQBPZ46+5+wBoDuouxMRkCPh/AFgGEFG9QE\nQhi6fA9cUAFAd8EGdRFO/MMipAFgfNHco5bopDSkrvuZ7woA/Iiqok6V9f6egrJwm8rnAwA0E2xQ\nlwXTFKu0GD9T0fcT4+cAgNAFG9RNhbgAh68VqOp+rmXv18SybVv2mmVzsAMA2ptMUIemKKz6COc6\n29DktVjCEgDCEk1QHzx0XBf27Rh7MyoNHcYAgGmLqtd3lbGbXMd+/zI+tyvUzwgAU2XOubG3Qaur\nq25jY6PWc+sGhc+m4zrvE9p2ddmGqtenJQAAujOzh51zq1XPi6bpO2RDBZuvanZZxzsqZgAIC0Hd\nUZ1ga9Opq837NKnui5atJKQBIDwE9YCqhpD5mA1syN8DAPQvuqBuOxY5FDFucxb3pwFgWNEFdVbf\nTbanTx6LPli7mkow+7j9AABjiDqolwlxprKQzGnfZC+2CGwAsZlsUEt+wnpqHa2GDKi2U5EOtQ0A\nEIPog3qqJ+Jl9+KnUA2O/ZmorAHEIroJT7JiG7/cRsxB0sd+a7M/mMAFQIjqTngS9RSiK5sXS3/m\nawawsSv2sd8/NG32R52Z2NjPAEIVdVDHsEgH/EuD1Xe4hhTWIW0LgHFFf4+6LyGdKKfag72PZuxl\n7xXSd1pH2fZO8VgAUI6gLhDiCf3A0RNa2byo9VNrY2+KF23DpktIxRjWRYqmfwUwXVE3fS+7R819\nx3CNGS6nTx4j3ABEJeqKuo971GOE+9SCI4Ym27Lx8SEM2+ICE0BWtEEde0eirvdnQwg9AsWPIdcy\nBxCfKJu+657YljWNx6SoGZ+Q7C6EfUhIA6hSu6I2s5dK2pB0zjn3HWZ2haT3Sdor6XFJNzrnLiTP\nvV3SLZJelHTEOfchz9tdaQ4nthCaaZsIqff62CHd5P1D2WcAxtGk6fuopM9K+prk77dJ+rBz7k4z\nuy35+5qZvUHSTZKuk3SVpIfM7PXOuRc9brd3Q8/pffDQcW89uMeYjrPL2tdjB88YM5V1GVYGYN5q\nBbWZ7ZH0jyXdIenHkodvkHQw+fM9ktYlrSWP3+ece17S583sMUlvlPRRb1tdIZaT28FDxyXpssAe\nu9qbsqEWCun6HcZyDAPoX92K+j9K+glJr8o8tss5dz7585OSdiV/3i3pY5nnnU0e86LOCbDJmOMx\nQnFl8+K2Huu+tmEq9+SLDPE9VXXWG+pYIaQBZFUGtZl9h6SnnXMPm9nBouc455yZNVrdw8xulXSr\nJF1zzTVNfrWWC/t2bAW2tL1qTSvZsaYg9fm+MZ7Um95bH+NiimF6AEJRp6J+k6R/YmbfLunlkr7G\nzH5d0lNmdqVz7ryZXSnp6eT55yRdnfn9Pclj2zjn7pZ0t7RYPavDZ1gqDcX0xJutOvOVbUymcFJf\nFthza/6fwvcJoB+VQe2cu13S7ZKUVNT/2jn3/Wb2s5JulnRn8v8PJr/ygKR7zexdWnQmu1bSJ3xt\ncNdwTX831mbiKZ7Q5xbKqSl+lwD86zLhyZ2S7jezWyQ9IelGSXLOPWJm90t6VNILkt7mu8e3j0o4\ntko6lJN67KFath/7/lyhfH8A4tMoqJ1z61r07pZz7ouS3lzyvDu06CGOjpqe4Lvcf/cVJtnX6TMA\n89tb9F51P1Od53V5fQBoy5zr7fZwbaurq25jY6PWcw8eOh5dNdxGlylGfd57bxu6faxURZgCmBIz\ne9g5t1r5PII6LEW91FNNgi87PC325mofCGsAoZl0UEvx3WNuq01lOtZ92BgQ2ABCUTeoo109ay5i\nC9cmTeVD3csGgJhFE9SH9xxZ/GG/t0nOJmusxTry77ds/vQmzwWAOYum6fvwniN6LgnpuTR7x4Qm\nZQBopm7TdxTrUWdDGmGiEgaAfkTT9I3tmHYTAOYhiooal5w+eay0mXnM5meavgGgH1TUkQg1CEPd\nLgCYiiiC+sGzd83+PnXdntwMeQKAaYmq6XvnmcVqmbGufOXDgaMntiZ9qbKsmRwAEIcoKmppUVVL\nl3qAx7yW9FBYEQoA4hdNUKfSwK5bVU5R1QXKEE3ehDQADCOaCU+KTH2Bjiormxe188y50SaCIawB\noL1ZzPXN6lAafba27L4nuAHAv6g6kxWhCTyMkAYA9CPKipp5vwEAcxF1RU1IL6xsXpz1kDUAmLIo\nK2pckr1YScOaCxgAmI6oe31zj3S7bFWddrSTht9PdCoDgGqTWuayDM29213Yt6OwmmaGMgCIV9RN\n34spRelUlndh346tKjoday1Jp5PJYqR+q+y685IDAKpFHdTMUlZtcQGzuJg5vOfI1uPZ0C7DrQUA\nGF/UQZ0a835srMoubrL7Mq2I2acAMJ6oO5OVIVgul7+fX3aroKxDWqrptK00fwNAsVlMIYr66oZr\n+ryVzYtbVXd6j7vNeuAHjp4grAGgg6h7fZdhApDu0h7kF/bt0HP7d7cK6RQtHADQ3iSDev3UWmGz\nLdopG/YFAOjfJO9R51HRhYNmcABYmMWEJ3UwdCssXDQBQDOTD2oAAGI2i6CmY1lYDhw9QWUNADXN\nIqgRJsIaAKpNPqjp/Q0AiNksen1LVG8xoEc4gDmh13cO96kBADGaTVCvn1ojrANHqwcAXG42QS1x\nvzoGhDUAbDeroJZoAgcAxGV2QY3wUVUDwCWzDGqq6vHU7dmdTopCaAOYu9kMz8pK5/9mRah+1Anj\nNgHM8C0AU1J3eNZXDbExoUk7lVGt+dWkWgYA1DPLpu8UFZofp08eY18CQE9mWVHDjzbh3KWaTn+X\niwIAc0JQo7Wy4KRpGwD8IajR2dDBnH0/qmsAUzfre9QS91cBAGGbfVCnCGwAQIgI6pyVzYtMiAIA\nCAZBnbN+ak07z5wjsAEAQSCoS+w8c27sTQAAgKAu8uDZu/Tg2bskMS84AGBcDM9a4lJVvVsSc4Nn\npR3vGDMNAP0KOqgP7zmy7e9plTu0nWfO6bn9u0d57xBle8ePFdj00AcwF0EHdd7hPUcGC+v8RYK0\naAanqi5WFJxU2wDQXdBB/eDZuy4LzL6r7KKAxiVjV7Jjvz8ADC3ooJYuD+t8E3RataWdvnaeOVc7\nvNPX3faayZ+zvb5p9g4DIY0xFLUM5W/5cGyiT8EHtVRcWZc1QT+3f/fWc8sCOx/QRa91Yd9r6fGd\nw8kIc7Ls1k3+Z/m/828FPkUR1HWkYbuyeXErgA8eOr5VGRdV5lX3m7M/n3Noh3LSOXD0RDDbgmnr\n2r+CYxU+RRPU+aAt69h1+WOL0K5TRS8z105koZ1saGpEyDgu0Ydoglq61JR98NDx2r9zKWB3FzyG\nvFhONDQ1IhQce+hbVEHdxYV9O2bdfF2lj5PN6ZPHBhuiRaWNIXGcYUjmnBt7G7S6uuo2NjZqPz+t\nqJtWxtmgpqq+XF8nnzHGU3MiRVtlxyvHFHwzs4edc6tVz5vNXN9U09X6CtR0rW9OdIgVxy7GNJug\nzlbQVNPjIbABoJnZBHUW1XW5oZqp+w5rpi+FL1xYYmyVncnM7OWSPiLpZcnzf8M59w4zu0LS+yTt\nlfS4pBudcxeS37ld0i2SXpR0xDn3Id8b3qQqJpibGapjVt+dzehgNn30/scc1Kmon5f0Lc65A5Ku\nl3TIzP6OpNskfdg5d62kDyd/l5m9QdJNkq6TdEjSL5nZS/vY+DoI6bANcWI9cPQEFTYqcZwgVJUV\ntVt0C/9y8tevTv5zkm6QdDB5/B5J65LWksfvc849L+nzZvaYpDdK+qjPDa+LYVntDVlZZ9+vL8wW\nFT+CFHNUaxx1UhE/LOmvS/pF59zHzWyXc+588pQnJe1K/rxb0scyv35W2dlGRsS60u1MpSlcojk8\nRoQz5q5WUDvnXpR0vZl9naQPmNn+3M+dmTUakG1mt0q6VZKuueaa2r938NDxVr22s6thoZ0hQm6o\nSVLqvMfK5kWtn1rrfVtwuTbHQNvjsuq9aInB2BpPeGJm/07SVyT9kKSDzrnzZnalpHXn3L6kI5mc\ncz+dPP9Dkt7pnCtt+q474UmbkM6ebA/vOdJ6rm9cbsiTV0hVFSft/jT9nn19F1Xvy3eOPtSd8KQy\nqM3sNZL+0jn3p2a2U9J/k3Rc0j+Q9EXn3J1mdpukK5xzP2Fm10m6V4v70ldp0dHs2qQqL1QnqJvO\nRpbel85XRIS1X3MNa4mTtw9DVs5VCGsMrW5Q12n6vlLSPcl96pdIut859ztm9lFJ95vZLZKekHSj\nJDnnHjGz+yU9KukFSW9bFtJNdA1pabGwR5NFPbDckPd8h+p0Vhf3u4c35r5mKBjGEs1c33Waveve\nU2w7VziWG/rEFUpgc8Luru53OdRwvjY4DtCUt6bvIVQFdVVIt/kHQlj7F8KJaqzwDuGzx6zJ9zbU\nvvZ1LHFsoMykFuUoC9Mu80anlTdjrKdlrLnE08kyQqnyp2wq09wCdQW/HnWfS86tn1rjfrVHIQ1j\nGfN+dpN711NqZp1zZ6wpfzaML+im72VN3j7/YbQdm41iIZ60Ql0Tu+t2hbCvh9y3Y33eOV+EoD8+\ne31PXlpZE9Z+hFRZj2XM+6hz3/d9GGoiHqBI0EGd3kce6h/IyuZFwtqT0IayZN9/qJnP8p95qOOY\n8PaLgMbYgg7qIaXzgBPW0zfkIiChCO3CqY05T66DeYsiqIcITxbrQEiGWqs7dDGMzY/xogdxiSKo\ni+4h+5oVivHUCEn2eA5tJrYpI6ARsiiC+sDRExJBCk9iC74mgRDbZyvTNQSz+2HZa421CAjQRBRB\n3acu96bzk6VQlZdjXuxivveHz1AaU5uRA2WfL3/sTWnsOuYh6HHUWWX/uLqsGVy32bts9rL8Gtes\nylXP2Ce8EAOrbcAOtcxj9v2G3n9DjEf3sQ1AU5MbR11WvV7Yt6NTtbYsVNP3TAP5wbN3LX2tw3uO\nJH+61DGN0A5PiPd+Q9qWVNm/p6GHnS0bbjbEfiOkMbZogjqtmtP1pIuaqn02ry5bKrNMPsiZnhR9\n6TM8mr52m23pGrB9BTShjBBF0/SdlYa1tHzBjirL/rF3aVJv8j5zE9KJMLTvJaR9M7YQm9cB3ya1\nelaer3vBy5r2fIX0sveZG/ZDsbFW/ArZUPuDfY8YRNP0nRfbvV/mCmYO8Lw57IuqGdFowgaqRRXU\nW5Oe9BzSfQ0lGnq+6RAxTGve6J0NNBfNPer0vnTTSrruP9w+171u875TF8IJdcx9H8LnH5LPfV12\nwTu3fYr4TW54lpQOk2oW1mnPa5/3nNHdHJvB5/Z5s7r0DF/2u3Pep5iPaCrqvLpX6EWTlSwL7bFW\nGZprZS2Ne7IdY78TLgCk+hV1tEEt6bKFOoo0DWqp+3q+2eFj0vIhZHMO6KyxwovmbwBjmUVQLxtP\nXTbtZ6pJVV1H2fvF1jt9TGMEVwgXSqGts8wFBDCMWQR1lUtTeqpWhYswzC2wQ7m9QkADwyKoM9LA\nZtGMeAwdGiFU1infn32sEQ0Alptkr++20jm4l80TDoSi7yFHBDQQl1kEdSoNbBbLCB8Toyy0HcbW\ntUMkgHBEOdd3V+un1rSyebGywxnGN1ST9JRCjJAGpmVWFXXW+qk1HTx0nGZwBK/J2H5m6gKmZ7ZB\njXiMNQlNqIqaw0PqDAfAr1n0+l6mzqQpCAvzr1eb+8UMEINJr0ftG/eq43Lg6IlewzTmkGN9ZWB6\nZh/U66fWtjqXIS5TqHwBoMrsgxooEmtVysULMD2zv0edlY6v5p51c3V7Ig/93l3FGHyxXmQAc8M9\n6hbShTpoBq9vZfNiZTD0fd+073vWMSGkgemhos6hqq5vZfNi5ZKheUMEah9hFeqFAMEMxIuKuqWd\nZ85Joqqu48K+HY0DbIheyaGGKgC0QUVdILvaFpV1PW2q66w+wnWoVajGREUNxIvVszrIrrYlEdZD\nSAPHZxi2XdAiJtn91fViCUCYCOoamA98OZ9hmH0tH6HtcxWu0yePBVlVpzhGgWmi6bsCncuqjTU8\nqs181z62NeSwlmgOB2JB07cnaVMic4KX67OJucnrhl7xAkAb9Pquaf3UGpXKErGMZY5lOwEgRUXd\n0OmTx6iulxizA1fVWszZn8+hoxmAaaCibiFdxCP9D+NrUyVPsbLm4gOYHirqloqGwUzxxN+Gz57W\nIeH7BTAGKmqPphZMiAvHHzBNDM/qCdXXJSEESJNhXm1eIwQh7GcA9dUdnkVQ9yz0k3uoxhjvXPWe\nMXyXhDUQD8ZRB8L3TFtzkd9XBBCAuaKiHgmh3U6TwJ7rPuaiBogDy1wGjqFd7TQJ3yGW1ASAvlFR\nB2Ku1V8XU6qusxdtXSfT4eIEiAMVdWTS6o+TbH11wzf0kJa2hzMtLQCy6EwWoKrpL3HJsqlAY9tn\nF/bt0MrmRe08c07rp+7a9rPYPgsAf2j6jhQn7unqeuFBqwwQB4ZnTVx6Mma9bEiL5vK0IgcwLVTU\nE8PKXtNRVRkv+66pqoHwMTPZzBHY0+F7xjRCHAgDTd8zl67uRWADQNwYnjUD3LecLjoVAtNHRT1x\n66fWtjqcYVqaztIGIE5U1DOwfmqNE/XEUEkD80FFPSPpEB7EhVAG5o2KekbWT62xGMhMEfZAvKio\nZ2bnmXPJn3ZLYqIUAAgd46hxGaqv6aKvAhAOVs9Ca5zMp+vA0RNciAGRqQxqM7vazH7fzB41s0fM\n7Gjy+BVm9ntm9kfJ/1cyv3O7mT1mZptm9m19fgD04/TJY9zLnrA0sNPQZggfEK7Kpm8zu1LSlc65\nT5nZqyQ9LOm7JL1V0rPOuTvN7DZJK865NTN7g6T3SnqjpKskPSTp9c65F8veg6bvMDGr2bzRsgL0\ny9sUos6585LOJ3/+kpl9VoueSDdIOpg87R5J65LWksfvc849L+nzZvaYFqH90eYfA2NjSNd85ZvI\nCW5gHI3uUZvZXknfJOnjknYlIS5JT0ralfx5t6Q/zvzaWaVdjLe/1q1mtmFmG88880zDzcYQ0h7i\nNIFD4v42MJbaw7PM7JWS3i/p7c65PzezrZ8555yZNeo+7py7W9Ld0qLpu8nvYhgPnr1LEvcvsV02\nrKmyERJfF5Jlx3X2XJgufDSEWkFtZl+tRUi/xzn3m8nDT5nZlc6588l97KeTx89Jujrz63uSxxCp\n7AFJRYUBQ9kgAAAWMklEQVSsOscDYY6++T4v5S9Gt/6e3AYcupWxTmcy0+Ie9LPOubdnHv9ZSV/M\ndCa7wjn3E2Z2naR7dakz2YclXUtnsvkgzFGEwEYfxjjfrGxe9FJR+1yP+k2S3iLpM2b26eSxn5R0\np6T7zewWSU9IulGSnHOPmNn9kh6V9IKkty0LaUxP9oRMaAOYirEuNuv0+v7vkqzkx28u+Z07JN3R\nYbswEduajQDAo7ncemFmMgAAAkZQo3dTuKIFEJY5tdQR1BjE6ZPHCGwAXswppCWWucTA0rCe2z80\nXPrOuWBDW3M9b7DMJUYz1390WCCwi4357yK076Ttvgjtc5SpOzyLoEZnh/cckXRpJrMmCGvEclL1\nLdRjf9n30XWb637XQ73P2AhqDCadVi+/eEfdfyys0gUpnpOrL03CqGyYY36fhRr+fYr5uCGoMaiy\nsM4q+gd14OgJVujCNr5mfQrdEGOApxzcMQd0yufMZECl7Im1rEIuO2kQ0sjiePBnyGbsIS8KphDS\nTVBRY1BTvsKHX3OorKv+PQwVSE3/XfZ5AdD0/WJGRY0gcU8NdU29sg7p2K871e/QgTnVgG6KoMao\nyu5bA9I01r5uczwP/VlDC+lYv+u+ENQIDpOioMgUQruOED/bEMOqQvzcoSCoESwCG2WmMstZ7Nvv\nA/ugGkGN4HFfG2UOHD0R7Yk+pu2u2s9Nx4SjGYIa0aHSBoY3l9nCQkRQI1rZf/iE9nzlv3sCIQx8\nD/4Q1JgEqmykQrh/3eU4rNtpzlfnOsY9h4+gxqTUHQ8K9GXZ8df0Pm/ZRQfH+LwQ1JgcqmsMqclx\n1vaY7PNY5uI2fAQ1Jot72PPWR4/wKR5HsX2mdAGgqU8vm0VQYxaosuep6/1qjpdmhro3PaeQlliU\nAzPEyRdZbS/ixr5vHNoSmHQga45FOQCgBl+TdQxxq6VpGJYNXfO5fQR0/6ioMTtU1OiiSTDN4Vgj\nqNujogaAHhSF71zDaq6fe2gENQB0NIfKGeMhqAEgAFX3uENbu51qejgENQDU4GNikLrh1vR5sfU4\nRzMENQAUyIdg2zAcItSGrLYJ6eER1ABQoOu46rH0NR48lM83RwQ1AHQUQoiFOI84/CCoMTmH9xy5\n7LHn9u/WhX07RtgaTN3YIT1E0I79GeeOoMakPbd/t5fXqTpRUZXM09ABRqexeSKoMTkPnr2rsKpu\nq87JipW65mXqATb1zxcbphDF5ORDOq2q6zZ9+zhJEdbTN6VqmmAeB1OIAi1wwkJdXZfQrPv6feA4\njwtBjcnbeeact3vVdfmYHANxKFuhytfrtUEQTwtN35i0tBm8bvO37xMcYY28Pi7iCOY40fQNaNGx\nLG/I8KSyRh7HA5qiosZsVZ0wfVYpBw8dZxw3ekNFHScqaqDCkAsarJ9a2/ozFRWAJghqzF7duY59\nvh8VNoC6CGogZ6hmxJXNi4Q1gErcowYCQHM4uuI+dXzq3qMmqIGAENjwjQAPF0ENRIzARl8I7nAQ\n1IBnRR3A+j7pEdjoG8E9HoZnAR4cPHRcUjKjWUHHr3yQrmxe3Pb37LCsNtLXo9MZ+pIe412PVfSH\noAYqNAnJ/HO7LtyQzlNOYKMvHFPho+kbKLCtkq6pzXCrsgBP3z/rwr4dBDZ6lx5jVNj94x410EGb\ne8NtJzLJhjUToSAkK5sXCewecY8aGFAatulJrU7QF1XThPR8pS0y+X4OqTGOjQv7dvS+7jaqEdRA\ngaarXjU9mZU9j9W2pq8siHeeOaf1U5ev9pYu1brzjLatq85F3XwQ1EBgmAt82tLvdWXzonaeOSep\neDnW1LKfNe1LUXaRkKKZO0wvGXsDgBC1rWp9NRNywkQd66fWtH5qTSubF7f9V6To8Z1nziWV/BrH\nXMCoqAHPfDVd0ww+fc/t371VVXeRhuzhPUe2Deer8/7p7y2r3DEuKmqgQAgzjhHS05UGqY+Qznrw\n7F1b1XFaLaf/lfG9DfCPihoo4aOiTQO/6esQ0tPVV0jnLauQ0w5qdZ6L8RHUQE+yVXk+9BnqMi9D\nhXNdBHNcCGqgJweOnigNZMamTl/2PnFVQGcrXEIUedyjBgocOHrCS/Nz9nWKQpkm7mlqEtJ5+WZp\ngKAGcvoIz6bBT6Udp/zwqLZN3YQ1smj6xmy0WWhjTAzPiouPgM5iyBRSVNSYjdAmdDh98hiV8wT4\nqqKLUFlDYvUszFAIVWqTgA5he1Gsr4Behip7OuqunkVFjagd3nNEh/cc0cFDxwvXcC4SQhVL+Mav\n7uxfQFfco0b0sisKHTx0fGkTd0gBuWz4FsLWZyVNxYy8yqA2s3dL+g5JTzvn9iePXSHpfZL2Snpc\n0o3OuQvJz26XdIukFyUdcc59qJctBzKyHcSKhkN1XY0qH6jZ6r3L6xLWcSGgMYY6FfWvSvoFSb+W\neew2SR92zt1pZrclf18zszdIuknSdZKukvSQmb3eOfei380Gqm2rnj2GtLT9JJ2uExxLb3K0Q0hj\nLJVB7Zz7iJntzT18g6SDyZ/vkbQuaS15/D7n3POSPm9mj0l6o6SP+tlcYLv8CW6opu2iEytrSE9X\n3/ej097dBDaKtL1Hvcs5dz7585OSdiV/3i3pY5nnnU0eA6LUpFk6vTce0n1wdDdGz24gq3Ovb7cY\n39V4jJeZ3WpmG2a28cwzz3TdDMC7oe4dE+xh6nN8dBGqaZRpW1E/ZWZXOufOm9mVkp5OHj8n6erM\n8/Ykj13GOXe3pLulxTjqltsBbONzNq8DR09oZfNi44lSVjYvNm4Cz29zGhA0pY+DgEZI2gb1A5Ju\nlnRn8v8PZh6/18zepUVnsmslfaLrRgJNtAnKZaqGfOWf6+O9CehxFN2Lzg7/y+oa4AQ06qqcmczM\n3qtFx7FXS3pK0jsk/Zak+yVdI+kJLYZnPZs8/6ck/YCkFyS93Tn3YNVGMDMZ+uCjss5XVstOrr47\nk/m+4EC5IScvqTqOMB91Zyar0+v7e0t+9OaS598h6Y6q1wVicGHfjkYncV/hmr0/zj1s/8acVYyQ\nRlPMTAZUuBS8ywcw+Oj1XdSBjdD2a8yQDm1hGMSBub6BmsY+yRLS3S0L6Z1nzvXacWzs4wfxoqIG\nMHksoIGYEdSYLF9DteqOpy7rTJYNCTqHDSsf0GNOWJLOD09ljaZo+gYqHDh6ojLw6/T4JqTHVTek\nmX0MoSGoMWk+ZxcrCuuDh47rwNETS0P4wr4djXuPF2GVrWmou246kKocRz0ExlGjb7GMca4bxnQs\nq6ftDGNlk5z0gabw7epcyExln3kbRw1Mwfqpta0TADOHISRTvHedDdtln6tp68KU9lETVNSYhcN7\njmxVSUOGbNsObVWVNRV1fW2q6i4VdVGYhNbc3SXwyj5L+pp1Q7ru601Z3YqaoMYspOv99h3WywK2\nabj6fK05K+obUBXYbYO6yZzwY/ERgEXb3/R1868xh2DOo+kbKLDzzLlB7z9m+VwRy+cqYVPnoyPf\n2MpCbKym4+zrpNtQp5JuWm1jgYoas5BW1FJ/VXXbjmDLQpfOZf7UHVPt40KuKoS6VNQE3HTQ9A3k\n5MPaZ1A3HTrlsxm87WvOUZ1mcJ8tLn0ENkE9HTR9AznpqkVpYPex0hXCln7fQzWFLwvi9VNrhC5q\noaLGbPkartU2qNtWwPQI96NNJ7MuWN4SeTR9AzX4XpJyqG2gKdyPoeYCJ6RRpG5QM4UoZivmIKsz\n//jpk8dolq+QTu8KhIygBkbWd6AS2NXSsH5u/+5ehu8d3nNkW2dGoAk6k2G20vCKubKGP9sr60VY\ns5IWQkBFjdnr0hls7JBvsg1U1fVlK2yfqKrRBp3JAIXRqazrdqQdo+oM+Rn7AiMmbVfgqkIHMzCO\nGmgglCk5u0wzmv7OgaMntLJ5kTG6nmS/iwv7Xrv1HbUN7b4Cuqxa54IgflTUQIE6Par74mvtbIZw\n9a9JaHcNzK7N5gR2eKiogUillfAQIRpKS0Ks0gU/0nvZ2cAeO5jbvl5+Br+6uBDoDxU1ELAh750T\n2O2U3cPuElxT6XRGeC9HRQ3MUJcmeYartVO0jCYBtZC/4GC/tENFDQSu6T3rMac2nbP8/erQmr5D\nQmAvMIUoMFO+ApYZzZrpa+z1FDFTWzMENTCC9ERV52TVZpiVz2qYsK6PecPro6quj6AGRrYsrA8e\nOr50TeNlaLoeDyFU7sGzd7F/GuIeNTCibEhnm0zXT63p8J4jlzWjtqnYfFTEhH59XfZ3nRaWoqb1\nseYkJ3C7YT1qIALpibnqvmbXJlXCehhd93NZUDe97913cBPQfjA8C4hAesJLm7dDvsfJ5CjV8vun\nyfzry2SDt05o55/DKmBxo6JGFLLNwFOew7rPICyr9tL3XFYN+prWdM7qDN+q2xN6WVjvPHOuMsy7\nBjcVtR9U1Ijetmbh/bu3LToh0RtZunwfLAv6LhcBhHR3l/bhIkTb3o8ufs3sY6+VpMsmYSl63aaB\nTUCPg6BG0J7LBHQegX25Ns3Ty/Zf/rWaXBigWD6w82G5LJybXDAVPbcovKuq86zsxQWhPRyavhGF\nuoEQe2j3GbJdf7/re6HYssq375aMZe+dqlt1E9zN0esbk9Q1fGJRN/iWrTvdJDyL9ll+/PayvgHc\nw+6mKDDH2J9VwU1TuV8ENSZtWQjNKail5WFd97WYHzwM+aAc++JnWXAT2t0R1JiFrlVjyLo0g/ts\nQm+KsO4uDcixg3qZsuU9yxDUl6PXN2ZhykszNu0YNsV9MFchB3Rq+zZWj9s+vOcIYd0SQY3ZOHD0\nRHRVdV+TjMS2HxCXsglXWJ+6HRblwCRMOXh8frYhlq5keUzkPbd/d+EwMJa6rId71JiUKQ/j6lpZ\nt/nM2fes8/vpDHIxNN2iWp3hW03lm8XnXFXXvUdNRY1ZOnD0xKzu6XYN6fTvTYaNIW58h+HgHjXQ\ngyaThYSoairSss+TrY4YWx0v3yFd1LlszpV0UzR9Y1LaVMl9hKjPsctjN3l3eW3COm5tArvJ+Orn\n9u9eukjJ1DE8Cwhc0f1f383xbS9C6vY2r3rO6WQiljndZpiS9CKrycQn2cCls5gfBDVmL4TFPUIc\nguVzaNiUx7tPTZ1pRH1Vv9mQZ5x1OYIaSEwpRHxddPQ1jhthKgrpbJg2DdIHz95FVe0Bvb6BCfLZ\nq31l82Kre5VFvxNbp7q5y94/blvt1v1dqulyVNSIno/OVlSN7bQJXprBw3Vh345tF1jpJCVpVdwl\nTNPfZU3r5uj1jaj12SO6zyBJ3zddSrLPntFdq9iy/cCKW9PluwkcxVg9C5PX5QTfJmTy6zNnNQna\n7HsPOXyp6Weus22sujU9TdakLqqS8z9DOYIak9DXybtrwJSFduhjhpt87roXEVXrYTdFYI+rSVDX\nQWCXI6gxGaGMLc7KBnXRWFOfgV0UhENOgjJ0ZS0R1iEoC2yC2h+CGpPi+8TtqwrMh1h6cvMV1H3e\nQ28brsvet89e3YT3sAjq/jEzGbDEhX07tqriNoGd3pPbKUm6fPk+H/oeytTH+tx9Th6TfU1COw6E\ntB9U1IhCHxV1qmlQpyGdDl1Jq+eicPKx3UP3TPc5B/mQ86gzzK67tlOFoh0qaqCEj/DIh7TUXzVZ\n9Lp9BlLdSrtOMA49PStjtLtJx1H7nCYU3RHUmLQ+AuK5/buX3oPuKySGDJ86Adtke3xV30ULmcCv\nC/t2aP0UIR0Smr4RtLbhNMT93bkaomptcxuBe9h++R52h8vR6xvRC2Vt6awhZhIDQpDeryas+8M9\nakQttJAuGjcNTFl+3m+Mh6BGUEIJ6MN7jmx1GJMIZ8xTOoyRqnpcLHOJqPUR0gcPHd/qMJb+B8zV\nhX07uOc/MipqBCOUk4GPYC5qMhwj8Olg5Q9DvzAWghpBmNLJL5SQzsu3Pkxpnw9h7P3V11zydQ09\nJh6X9BbUZnZI0klJL5X0X51zd/b1XghH0cmsj5m1YjhZpDM5Pbd/t1Y2L452ci3aVwcPHZcCuHhA\nO77nlG+ij6lnsVwvQW1mL5X0i5L+kaSzkj5pZg845x7t4/0QhrLQ7RrGY1cyqaZTVGZndjp46HhQ\nYb24iNi99eeqSVwQpjGOKXqCD6+vivqNkh5zzn1OkszsPkk3SCKoJ8hnkBZdqYdw9Z6enOpePOR7\nya6fWitdw3oI+bB+8OxdW3OWp8a6kEB7hPQ89BXUuyX9cebvZyV9c0/vhcCEEKxdZAOr7MRUFmpV\nJ7JQwzA/t3MorRgxGLPlJ9TjCX71MjOZmX2PpEPOuR9M/v4WSd/snPuRzHNulXRr8tf9ks5435B5\neLWkL4y9EZFi37XHvmuPfdfe1PbdNzjnXlP1pL4q6nOSrs78fU/y2Bbn3N2S7pYkM9uoM40aLse+\na4991x77rj32XXtz3Xd9TXjySUnXmtlrzWyHpJskPdDTewEAMFm9VNTOuRfM7EckfUiL4Vnvds49\n0sd7AQAwZb2No3bO/a6k36359Lv72o4ZYN+1x75rj33XHvuuvVnuuyCWuQQAAMVYlAMAgICNHtRm\ndsjMNs3sMTO7beztCZ2ZPW5mnzGzT5vZRvLYFWb2e2b2R8n/V8bezhCY2bvN7GkzO5N5rHRfmdnt\nyXG4aWbfNs5Wj69kv73TzM4lx92nzezbMz9jvyXM7Goz+30ze9TMHjGzo8njHHcVluy72R97ozZ9\nJ1ON/qEyU41K+l6mGi1nZo9LWnXOfSHz2M9IetY5d2dysbPinJv9ArJm9vclfVnSrznn9iePFe4r\nM3uDpPdqMaveVZIekvR659yLI23+aEr22zslfdk593O557LfMszsSklXOuc+ZWavkvSwpO+S9FZx\n3C21ZN/dqJkfe2NX1FtTjTrnLkpKpxpFMzdIuif58z1aHNyz55z7iKRncw+X7asbJN3nnHveOfd5\nSY9pcXzOTsl+K8N+y3DOnXfOfSr585ckfVaLmRo57ios2XdlZrPvxg7qoqlGl30xkJykh8zs4WR2\nN0na5Zw7n/z5SUm7xtm0KJTtK47Faj9qZn+QNI2nTbfstxJmtlfSN0n6uDjuGsntO2nmx97YQY3m\n/q5z7npJhyW9LWmm3OIW9zLoyl8D+6qR/yzpdZKul3Re0s+PuzlhM7NXSnq/pLc75/48+zOOu+UK\n9t3sj72xg7pyqlFs55w7l/z/aUkf0KKp56nk/k56n+fp8bYweGX7imNxCefcU865F51z/0/Sf9Gl\nJkb2W46ZfbUWQfMe59xvJg9z3NVQtO849sYPaqYabcDMXpF0spCZvULSt2qxmMkDkm5OnnazpA+O\ns4VRKNtXD0i6ycxeZmavlXStpE+MsH1BSkMm8U91aREd9luGmZmkX5H0WefcuzI/4rirULbvOPZ6\nnJmsDqYabWyXpA8sjmd9laR7nXOnzOyTku43s1skPaFFL8nZM7P3Sjoo6dVmdlbSOyTdqYJ95Zx7\nxMzu12LN9BckvW2KvUfrKNlvB83sei2abB+X9MMS+63AmyS9RdJnzOzTyWM/KY67Osr23ffO/dhj\nZjIAAAI2dtM3AABYgqAGACBgBDUAAAEjqAEACBhBDQBAwAhqAAACRlADABAwghoAgID9fwhJGkmE\nkcfBAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "animate_one_density(Jun04, 1023.3)" ] }, { "cell_type": "code", "execution_count": 50, "metadata": { "collapsed": true }, "outputs": [], "source": [ "JunW04 = nc.Dataset('/ocean/vdo/MEOPAR/completed-runs/SalishSeaLake/LakeJun-0.4/SalishSea_1h_20160601_20160607_grid_T.nc')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# June, vtau = -0.4" ] }, { "cell_type": "code", "execution_count": 51, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/html": [ "" ], "text/plain": [ "" ] }, "execution_count": 51, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAeoAAAHiCAYAAAAnPo9XAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3X+wHedd3/HPNwElahLgukmFLdkooY5mYrUy7Z3QmfSH\nhrQgUagpZTymkHGKwXQakCJouTZMm8x0XCx+RJULdHBLBjPEcTyEEMNgpTjDnUzb/LpOIyI7veAm\ndpGQfyQWkBCPhd2nf5zdq72r3bO/nt19nt33a8Zj6dxzz9mzZ7Wf/T77/DDnnAAAQJheMvYGAACA\ncgQ1AAABI6gBAAgYQQ0AQMAIagAAAkZQAwAQMIIaGImZPW5m/3Ds7QAQNoIaQG22cNzMvpj8d9zM\nLPPzvWb2+2b2FTP73/kLETP752b2hJn9hZn9lpldMfynAOJCUAMTZ2Zf5fHlbpX0XZIOSPqbkr5T\n0g9nfv5eSf9L0l+V9FOSfsPMXpNsx3WSflnSWyTtkvQVSb/kcduASSKogXFdb2Z/YGZ/ZmbvM7OX\npz8wsx8ys8fM7Fkze8DMrkoe32tmLhvAZrZuZj+Y/PmtZvY/zOyEmX1R0js9bu/Nkn7eOXfWOXdO\n0s9Jemvyvq+X9LckvcM595xz7v2S/kDSP0t+9/sk/bZz7iPOuS9L+reSvtvMXuVx+4DJIaiBcd0o\n6ZCk12pRob5VkszsWyT9dPLzKyU9Iem+Bq/7zZI+p0Xlekf+h0kT9J8u+e+akte9TtLpzN9PJ4+l\nP/ucc+5LS36+9bvOuf8j6XlJr2/wuYDZ8dkkBqC5u5xzfyJJZvbbkq5PHv8+Se92zn0q+dntki6Y\n2d6ar/snzrn/lPz5hfwPnXP3Srq3xfa+UtKfZf7+55Jemdynzv8s/fnukt9Nf05FDSxBRQ2M68nM\nn7+iRZhJ0lVaVNGSpKSp+Iu6FHpV/tjL1l3uy5K+JvP3r5X0ZbdY3Sf/s/TnXyr53fzPARQgqIEw\n/Ymkb0j/Ymav0KKD1jlJf5E8/Fcyz//63O8vXRbPzL7PzL685L+ypu9HtOhIljqQPJb+7HW5e875\nn2/9rpl9o6Qdkv5w2bYCc0dQA2F6r6R/YWbXm9nLJP0HSR93zj3unHtGi8D+fjN7qZn9gKRvbPLi\nzrn3OOdeueS//1vyq78m6cfMbLeZ7Zb045J+NXnNP5T0aUnvMLOXm9l3S/obkt6f/O57JH2nmf29\n5MLj30v6zdw9bQA5BDUQIOfcQ1r0in6/pPNaBPFNmaf8kKR/o0Vz+HWS/udAm/bLkn5b0meS/34n\neSx1k6RVSRe06Az3PcmFhZxzj0j6l1oE9tOSXiHpXw203UC0bHFrCQAAhIiKGgCAgBHUAAAEjKAG\nACBgBDUAAAEjqAEACFgQU4i++tWvdnv37h17MwBM3OYfXZoIbt+1+TligGE9/PDDX3DOvabqeUEE\n9d69e7WxsTH2ZgCYsANHT+h1b7r097+UdPrksdG2BzCzJ6qfRdM3gBk4cPREo8eBkBDUACaNMEbs\ngmj6BgDf6gZ0/nk0hyM0BDWAyfBRPR88dFzrp9Y8bA3gB0ENIHg0X2POCGoAQRornKmmERo6kwFA\ngvvTCBEVNYBZKArhbNVOSCNUBDWAVsqapusEXh/N2m2ClnBGDAhqAKXaBOqQ95YJWswBQQ3gMvSy\nBsJBUAOIDpU05oSgBhA0QhlzR1ADCBIBDSwQ1AAuc/rksd7uUxPAQDMENYBCaaDWDWwCGOgHM5MB\nKEVIA+MjqAEUYogWEAaavgFsQ0ADYaGiBrCFkAbCQ1AD6IyAB/pjzrmxt0Grq6tuY2Nj7M0AkNEl\nfOlcBlQzs4edc6uVzyOoASxDYAP9qBvUNH0DWOr0yWMELjAighpALYQ1MA6CGgCAgBHUAGprWlXT\nGxzojqAG0AhN4MCwmJkMQCGqYSAMBDUASQQzECqCGpg5AhoIG0ENzNQQAc39bKA7gjpS+ZMsJ0SE\nhOMR8IegBmaor2qagAb8I6iBiSOUgbgR1AAqEcrAeJjwJELcn0ZdPqppji9gXFTUESk66XISRRlC\nGpgGgjoCBDQAzBdBHQgmnUBIuBAEwkFQj4yARl9OnzzW+PgioIHw0JkMmChCGpgGc86NvQ161dfu\nca9764/Xeu7cTibLTrZz2xeop00rDccSMDwze9g5t1r1vFpN32b2uKQvSXpR0gvOuVUzu0LS+yTt\nlfS4pBudcxeS598u6Zbk+Ueccx9q8RkKVZ2EpnbCST/PkE3k2fea2v6csibHCN8rEI9aFXUS1KvO\nuS9kHvsZSc865+40s9skrTjn1szsDZLeK+mNkq6S9JCk1zvnXix7/SYVdR1TPQn56v3NCX1a+D6B\nOHmtqEvcIOlg8ud7JK1LWksev88597ykz5vZY1qE9kfLXmjftV+vjQaVY9HJZg6dsrKfO/286f+X\nnYDnsG/moMv3SEAD8aob1E7SQ2b2oqRfds7dLWmXc+588vMnJe1K/rxb0scyv3s2eaxXczsR5Xv0\nHjh64rJ9QGei8LS9f8zFFjBfdZu+dzvnzpnZX5P0e5J+VNIDzrmvyzzngnNuxcx+QdLHnHO/njz+\nK5IedM79Ru41b5V0qyRdc801f/uJJ57Y9p5zuxfdlM8T99z3pQ8hBynfLxAmr03fzrlzyf+fNrMP\naNGU/ZSZXemcO29mV0p6Onn6OUlXZ359T/JY/jXvlnS3JK2url52tVBVRRRVkHNAQIcn1JDm+wWm\noXIctZm9wsxelf5Z0rdKOiPpAUk3J0+7WdIHkz8/IOkmM3uZmb1W0rWSPtFm406fPFZ57zXUk2Qf\nfH/WOe27voS6DwlpYDrqVNS7JH3AzNLn3+ucO2Vmn5R0v5ndIukJSTdKknPuETO7X9Kjkl6Q9LZl\nPb7rqFNd558/NV0+Y6hhAv+meOwDc1cZ1M65z0k6UPD4FyW9ueR37pB0R+etyxhjPPGYhvqcdXqN\nh6Dr/uh6YVPU4z4koX9/ANpjru+R1T3pdzkR1wmZWAK7ra7h6iucfQ8tnOr3BeCS6IK6TmUdysxa\ny5qrhwjoKQixem1j2fc49+8YwHJBzPW9urrqNjY2Or1GaHNiD9lU20Ro+2mZOYQ0gPkaYmayoCyr\ntIcKzZgnGAlpW6YS0ADgw2SCuk9lwdFlxqixgjGkQC4yhZAOfR8DiMvkgrrrNJpNTCFUpsZXSMZ0\niwDAtE0uqPNCG1bDSb4/Pvct3xOAUEw+qFNjLAUYwoVBTGK7jQAAQ5hNUA/N19rRUxdq73gACMUs\ngpoxy+OgRQEAupt8UNcJiz4COj9cjIuA7tiHAOZo0kEdQkU353DxNT/7nPchAEw6qOsgBPrj40KJ\n7wfA3M0+qBEmAhoAFiYZ1CE0eaMdAhoAtnvJ2BswtgNHTxDsPWkSuqdPHiOkAaDAJCvqNg4cPUFQ\nZPjqrT7klK4AMEWTDOomvY0J5+V8XsAQ0gDQ3CSDOlW1uhUhXU+X6pqLJQDoZtJBXRUS+Z8TGMv1\nURGzzwFguUkHNcJFQANAPZMNau6HhodwBoDmJhnUVSFddu+ant/tsM8AoD+TDOo6yoYNsYhGM+wn\nAOiXOefG3gatrq66jY2Nzq9Tp5Ju+hpzDaKmtw7mup8AoC0ze9g5t1r1vMnMTOYjpHEJ+wsAwjD5\npu+mgUNAXcI4dAAY3+SDGu302Wu+6LUJfQAoNomgZijWOHx2vMt+h4Q2AFwymXvUGE/dFciaPA8A\nsBB9Rc1JvR9V96eL+PouqKgB4JJog5qAngZCGQCWizKo64Y0IdDesvHkPi6S+G4AoJ4og7oupgT1\no0lIs78BwK8ogzoNg6KqjyZxP9oELiENAP5FGdRSedMsYQEAmJLohmfVHeKDfiy7EOK7AQD/ogtq\nDC8fwLRaAMBwogrqsqknCY5h1K2WqawBwJ9ogrronjQBPZ46+5+wBoDuouxMRkCPh/AFgGEFG9QE\nQhi6fA9cUAFAd8EGdRFO/MMipAFgfNHco5bopDSkrvuZ7woA/Iiqok6V9f6egrJwm8rnAwA0E2xQ\nlwXTFKu0GD9T0fcT4+cAgNAFG9RNhbgAh68VqOp+rmXv18SybVv2mmVzsAMA2ptMUIemKKz6COc6\n29DktVjCEgDCEk1QHzx0XBf27Rh7MyoNHcYAgGmLqtd3lbGbXMd+/zI+tyvUzwgAU2XOubG3Qaur\nq25jY6PWc+sGhc+m4zrvE9p2ddmGqtenJQAAujOzh51zq1XPi6bpO2RDBZuvanZZxzsqZgAIC0Hd\nUZ1ga9Opq837NKnui5atJKQBIDwE9YCqhpD5mA1syN8DAPQvuqBuOxY5FDFucxb3pwFgWNEFdVbf\nTbanTx6LPli7mkow+7j9AABjiDqolwlxprKQzGnfZC+2CGwAsZlsUEt+wnpqHa2GDKi2U5EOtQ0A\nEIPog3qqJ+Jl9+KnUA2O/ZmorAHEIroJT7JiG7/cRsxB0sd+a7M/mMAFQIjqTngS9RSiK5sXS3/m\nawawsSv2sd8/NG32R52Z2NjPAEIVdVDHsEgH/EuD1Xe4hhTWIW0LgHFFf4+6LyGdKKfag72PZuxl\n7xXSd1pH2fZO8VgAUI6gLhDiCf3A0RNa2byo9VNrY2+KF23DpktIxRjWRYqmfwUwXVE3fS+7R819\nx3CNGS6nTx4j3ABEJeqKuo971GOE+9SCI4Ym27Lx8SEM2+ICE0BWtEEde0eirvdnQwg9AsWPIdcy\nBxCfKJu+657YljWNx6SoGZ+Q7C6EfUhIA6hSu6I2s5dK2pB0zjn3HWZ2haT3Sdor6XFJNzrnLiTP\nvV3SLZJelHTEOfchz9tdaQ4nthCaaZsIqff62CHd5P1D2WcAxtGk6fuopM9K+prk77dJ+rBz7k4z\nuy35+5qZvUHSTZKuk3SVpIfM7PXOuRc9brd3Q8/pffDQcW89uMeYjrPL2tdjB88YM5V1GVYGYN5q\nBbWZ7ZH0jyXdIenHkodvkHQw+fM9ktYlrSWP3+ece17S583sMUlvlPRRb1tdIZaT28FDxyXpssAe\nu9qbsqEWCun6HcZyDAPoX92K+j9K+glJr8o8tss5dz7585OSdiV/3i3pY5nnnU0e86LOCbDJmOMx\nQnFl8+K2Huu+tmEq9+SLDPE9VXXWG+pYIaQBZFUGtZl9h6SnnXMPm9nBouc455yZNVrdw8xulXSr\nJF1zzTVNfrWWC/t2bAW2tL1qTSvZsaYg9fm+MZ7Um95bH+NiimF6AEJRp6J+k6R/YmbfLunlkr7G\nzH5d0lNmdqVz7ryZXSnp6eT55yRdnfn9Pclj2zjn7pZ0t7RYPavDZ1gqDcX0xJutOvOVbUymcFJf\nFthza/6fwvcJoB+VQe2cu13S7ZKUVNT/2jn3/Wb2s5JulnRn8v8PJr/ygKR7zexdWnQmu1bSJ3xt\ncNdwTX831mbiKZ7Q5xbKqSl+lwD86zLhyZ2S7jezWyQ9IelGSXLOPWJm90t6VNILkt7mu8e3j0o4\ntko6lJN67KFath/7/lyhfH8A4tMoqJ1z61r07pZz7ouS3lzyvDu06CGOjpqe4Lvcf/cVJtnX6TMA\n89tb9F51P1Od53V5fQBoy5zr7fZwbaurq25jY6PWcw8eOh5dNdxGlylGfd57bxu6faxURZgCmBIz\ne9g5t1r5PII6LEW91FNNgi87PC325mofCGsAoZl0UEvx3WNuq01lOtZ92BgQ2ABCUTeoo109ay5i\nC9cmTeVD3csGgJhFE9SH9xxZ/GG/t0nOJmusxTry77ds/vQmzwWAOYum6fvwniN6LgnpuTR7x4Qm\nZQBopm7TdxTrUWdDGmGiEgaAfkTT9I3tmHYTAOYhiooal5w+eay0mXnM5meavgGgH1TUkQg1CEPd\nLgCYiiiC+sGzd83+PnXdntwMeQKAaYmq6XvnmcVqmbGufOXDgaMntiZ9qbKsmRwAEIcoKmppUVVL\nl3qAx7yW9FBYEQoA4hdNUKfSwK5bVU5R1QXKEE3ehDQADCOaCU+KTH2Bjiormxe188y50SaCIawB\noL1ZzPXN6lAafba27L4nuAHAv6g6kxWhCTyMkAYA9CPKipp5vwEAcxF1RU1IL6xsXpz1kDUAmLIo\nK2pckr1YScOaCxgAmI6oe31zj3S7bFWddrSTht9PdCoDgGqTWuayDM29213Yt6OwmmaGMgCIV9RN\n34spRelUlndh346tKjoday1Jp5PJYqR+q+y685IDAKpFHdTMUlZtcQGzuJg5vOfI1uPZ0C7DrQUA\nGF/UQZ0a835srMoubrL7Mq2I2acAMJ6oO5OVIVgul7+fX3aroKxDWqrptK00fwNAsVlMIYr66oZr\n+ryVzYtbVXd6j7vNeuAHjp4grAGgg6h7fZdhApDu0h7kF/bt0HP7d7cK6RQtHADQ3iSDev3UWmGz\nLdopG/YFAOjfJO9R51HRhYNmcABYmMWEJ3UwdCssXDQBQDOTD2oAAGI2i6CmY1lYDhw9QWUNADXN\nIqgRJsIaAKpNPqjp/Q0AiNksen1LVG8xoEc4gDmh13cO96kBADGaTVCvn1ojrANHqwcAXG42QS1x\nvzoGhDUAbDeroJZoAgcAxGV2QY3wUVUDwCWzDGqq6vHU7dmdTopCaAOYu9kMz8pK5/9mRah+1Anj\nNgHM8C0AU1J3eNZXDbExoUk7lVGt+dWkWgYA1DPLpu8UFZofp08eY18CQE9mWVHDjzbh3KWaTn+X\niwIAc0JQo7Wy4KRpGwD8IajR2dDBnH0/qmsAUzfre9QS91cBAGGbfVCnCGwAQIgI6pyVzYtMiAIA\nCAZBnbN+ak07z5wjsAEAQSCoS+w8c27sTQAAgKAu8uDZu/Tg2bskMS84AGBcDM9a4lJVvVsSc4Nn\npR3vGDMNAP0KOqgP7zmy7e9plTu0nWfO6bn9u0d57xBle8ePFdj00AcwF0EHdd7hPUcGC+v8RYK0\naAanqi5WFJxU2wDQXdBB/eDZuy4LzL6r7KKAxiVjV7Jjvz8ADC3ooJYuD+t8E3RataWdvnaeOVc7\nvNPX3faayZ+zvb5p9g4DIY0xFLUM5W/5cGyiT8EHtVRcWZc1QT+3f/fWc8sCOx/QRa91Yd9r6fGd\nw8kIc7Ls1k3+Z/m/828FPkUR1HWkYbuyeXErgA8eOr5VGRdV5lX3m7M/n3Noh3LSOXD0RDDbgmnr\n2r+CYxU+RRPU+aAt69h1+WOL0K5TRS8z105koZ1saGpEyDgu0Ydoglq61JR98NDx2r9zKWB3FzyG\nvFhONDQ1IhQce+hbVEHdxYV9O2bdfF2lj5PN6ZPHBhuiRaWNIXGcYUjmnBt7G7S6uuo2NjZqPz+t\nqJtWxtmgpqq+XF8nnzHGU3MiRVtlxyvHFHwzs4edc6tVz5vNXN9U09X6CtR0rW9OdIgVxy7GNJug\nzlbQVNPjIbABoJnZBHUW1XW5oZqp+w5rpi+FL1xYYmyVncnM7OWSPiLpZcnzf8M59w4zu0LS+yTt\nlfS4pBudcxeS37ld0i2SXpR0xDn3Id8b3qQqJpibGapjVt+dzehgNn30/scc1Kmon5f0Lc65A5Ku\nl3TIzP6OpNskfdg5d62kDyd/l5m9QdJNkq6TdEjSL5nZS/vY+DoI6bANcWI9cPQEFTYqcZwgVJUV\ntVt0C/9y8tevTv5zkm6QdDB5/B5J65LWksfvc849L+nzZvaYpDdK+qjPDa+LYVntDVlZZ9+vL8wW\nFT+CFHNUaxx1UhE/LOmvS/pF59zHzWyXc+588pQnJe1K/rxb0scyv35W2dlGRsS60u1MpSlcojk8\nRoQz5q5WUDvnXpR0vZl9naQPmNn+3M+dmTUakG1mt0q6VZKuueaa2r938NDxVr22s6thoZ0hQm6o\nSVLqvMfK5kWtn1rrfVtwuTbHQNvjsuq9aInB2BpPeGJm/07SVyT9kKSDzrnzZnalpHXn3L6kI5mc\ncz+dPP9Dkt7pnCtt+q474UmbkM6ebA/vOdJ6rm9cbsiTV0hVFSft/jT9nn19F1Xvy3eOPtSd8KQy\nqM3sNZL+0jn3p2a2U9J/k3Rc0j+Q9EXn3J1mdpukK5xzP2Fm10m6V4v70ldp0dHs2qQqL1QnqJvO\nRpbel85XRIS1X3MNa4mTtw9DVs5VCGsMrW5Q12n6vlLSPcl96pdIut859ztm9lFJ95vZLZKekHSj\nJDnnHjGz+yU9KukFSW9bFtJNdA1pabGwR5NFPbDckPd8h+p0Vhf3u4c35r5mKBjGEs1c33Waveve\nU2w7VziWG/rEFUpgc8Luru53OdRwvjY4DtCUt6bvIVQFdVVIt/kHQlj7F8KJaqzwDuGzx6zJ9zbU\nvvZ1LHFsoMykFuUoC9Mu80anlTdjrKdlrLnE08kyQqnyp2wq09wCdQW/HnWfS86tn1rjfrVHIQ1j\nGfN+dpN711NqZp1zZ6wpfzaML+im72VN3j7/YbQdm41iIZ60Ql0Tu+t2hbCvh9y3Y33eOV+EoD8+\ne31PXlpZE9Z+hFRZj2XM+6hz3/d9GGoiHqBI0EGd3kce6h/IyuZFwtqT0IayZN9/qJnP8p95qOOY\n8PaLgMbYgg7qIaXzgBPW0zfkIiChCO3CqY05T66DeYsiqIcITxbrQEiGWqs7dDGMzY/xogdxiSKo\ni+4h+5oVivHUCEn2eA5tJrYpI6ARsiiC+sDRExJBCk9iC74mgRDbZyvTNQSz+2HZa421CAjQRBRB\n3acu96bzk6VQlZdjXuxivveHz1AaU5uRA2WfL3/sTWnsOuYh6HHUWWX/uLqsGVy32bts9rL8Gtes\nylXP2Ce8EAOrbcAOtcxj9v2G3n9DjEf3sQ1AU5MbR11WvV7Yt6NTtbYsVNP3TAP5wbN3LX2tw3uO\nJH+61DGN0A5PiPd+Q9qWVNm/p6GHnS0bbjbEfiOkMbZogjqtmtP1pIuaqn02ry5bKrNMPsiZnhR9\n6TM8mr52m23pGrB9BTShjBBF0/SdlYa1tHzBjirL/rF3aVJv8j5zE9KJMLTvJaR9M7YQm9cB3ya1\nelaer3vBy5r2fIX0sveZG/ZDsbFW/ArZUPuDfY8YRNP0nRfbvV/mCmYO8Lw57IuqGdFowgaqRRXU\nW5Oe9BzSfQ0lGnq+6RAxTGve6J0NNBfNPer0vnTTSrruP9w+171u875TF8IJdcx9H8LnH5LPfV12\nwTu3fYr4TW54lpQOk2oW1mnPa5/3nNHdHJvB5/Z5s7r0DF/2u3Pep5iPaCrqvLpX6EWTlSwL7bFW\nGZprZS2Ne7IdY78TLgCk+hV1tEEt6bKFOoo0DWqp+3q+2eFj0vIhZHMO6KyxwovmbwBjmUVQLxtP\nXTbtZ6pJVV1H2fvF1jt9TGMEVwgXSqGts8wFBDCMWQR1lUtTeqpWhYswzC2wQ7m9QkADwyKoM9LA\nZtGMeAwdGiFU1infn32sEQ0Alptkr++20jm4l80TDoSi7yFHBDQQl1kEdSoNbBbLCB8Toyy0HcbW\ntUMkgHBEOdd3V+un1rSyebGywxnGN1ST9JRCjJAGpmVWFXXW+qk1HTx0nGZwBK/J2H5m6gKmZ7ZB\njXiMNQlNqIqaw0PqDAfAr1n0+l6mzqQpCAvzr1eb+8UMEINJr0ftG/eq43Lg6IlewzTmkGN9ZWB6\nZh/U66fWtjqXIS5TqHwBoMrsgxooEmtVysULMD2zv0edlY6v5p51c3V7Ig/93l3FGHyxXmQAc8M9\n6hbShTpoBq9vZfNiZTD0fd+073vWMSGkgemhos6hqq5vZfNi5ZKheUMEah9hFeqFAMEMxIuKuqWd\nZ85Joqqu48K+HY0DbIheyaGGKgC0QUVdILvaFpV1PW2q66w+wnWoVajGREUNxIvVszrIrrYlEdZD\nSAPHZxi2XdAiJtn91fViCUCYCOoamA98OZ9hmH0tH6HtcxWu0yePBVlVpzhGgWmi6bsCncuqjTU8\nqs181z62NeSwlmgOB2JB07cnaVMic4KX67OJucnrhl7xAkAb9Pquaf3UGpXKErGMZY5lOwEgRUXd\n0OmTx6iulxizA1fVWszZn8+hoxmAaaCibiFdxCP9D+NrUyVPsbLm4gOYHirqloqGwUzxxN+Gz57W\nIeH7BTAGKmqPphZMiAvHHzBNDM/qCdXXJSEESJNhXm1eIwQh7GcA9dUdnkVQ9yz0k3uoxhjvXPWe\nMXyXhDUQD8ZRB8L3TFtzkd9XBBCAuaKiHgmh3U6TwJ7rPuaiBogDy1wGjqFd7TQJ3yGW1ASAvlFR\nB2Ku1V8XU6qusxdtXSfT4eIEiAMVdWTS6o+TbH11wzf0kJa2hzMtLQCy6EwWoKrpL3HJsqlAY9tn\nF/bt0MrmRe08c07rp+7a9rPYPgsAf2j6jhQn7unqeuFBqwwQB4ZnTVx6Mma9bEiL5vK0IgcwLVTU\nE8PKXtNRVRkv+66pqoHwMTPZzBHY0+F7xjRCHAgDTd8zl67uRWADQNwYnjUD3LecLjoVAtNHRT1x\n66fWtjqcYVqaztIGIE5U1DOwfmqNE/XEUEkD80FFPSPpEB7EhVAG5o2KekbWT62xGMhMEfZAvKio\nZ2bnmXPJn3ZLYqIUAAgd46hxGaqv6aKvAhAOVs9Ca5zMp+vA0RNciAGRqQxqM7vazH7fzB41s0fM\n7Gjy+BVm9ntm9kfJ/1cyv3O7mT1mZptm9m19fgD04/TJY9zLnrA0sNPQZggfEK7Kpm8zu1LSlc65\nT5nZqyQ9LOm7JL1V0rPOuTvN7DZJK865NTN7g6T3SnqjpKskPSTp9c65F8veg6bvMDGr2bzRsgL0\ny9sUos6585LOJ3/+kpl9VoueSDdIOpg87R5J65LWksfvc849L+nzZvaYFqH90eYfA2NjSNd85ZvI\nCW5gHI3uUZvZXknfJOnjknYlIS5JT0ralfx5t6Q/zvzaWaVdjLe/1q1mtmFmG88880zDzcYQ0h7i\nNIFD4v42MJbaw7PM7JWS3i/p7c65PzezrZ8555yZNeo+7py7W9Ld0qLpu8nvYhgPnr1LEvcvsV02\nrKmyERJfF5Jlx3X2XJgufDSEWkFtZl+tRUi/xzn3m8nDT5nZlc6588l97KeTx89Jujrz63uSxxCp\n7AFJRYUBQ9kgAAAWMklEQVSsOscDYY6++T4v5S9Gt/6e3AYcupWxTmcy0+Ie9LPOubdnHv9ZSV/M\ndCa7wjn3E2Z2naR7dakz2YclXUtnsvkgzFGEwEYfxjjfrGxe9FJR+1yP+k2S3iLpM2b26eSxn5R0\np6T7zewWSU9IulGSnHOPmNn9kh6V9IKkty0LaUxP9oRMaAOYirEuNuv0+v7vkqzkx28u+Z07JN3R\nYbswEduajQDAo7ncemFmMgAAAkZQo3dTuKIFEJY5tdQR1BjE6ZPHCGwAXswppCWWucTA0rCe2z80\nXPrOuWBDW3M9b7DMJUYz1390WCCwi4357yK076Ttvgjtc5SpOzyLoEZnh/cckXRpJrMmCGvEclL1\nLdRjf9n30XWb637XQ73P2AhqDCadVi+/eEfdfyys0gUpnpOrL03CqGyYY36fhRr+fYr5uCGoMaiy\nsM4q+gd14OgJVujCNr5mfQrdEGOApxzcMQd0yufMZECl7Im1rEIuO2kQ0sjiePBnyGbsIS8KphDS\nTVBRY1BTvsKHX3OorKv+PQwVSE3/XfZ5AdD0/WJGRY0gcU8NdU29sg7p2K871e/QgTnVgG6KoMao\nyu5bA9I01r5uczwP/VlDC+lYv+u+ENQIDpOioMgUQruOED/bEMOqQvzcoSCoESwCG2WmMstZ7Nvv\nA/ugGkGN4HFfG2UOHD0R7Yk+pu2u2s9Nx4SjGYIa0aHSBoY3l9nCQkRQI1rZf/iE9nzlv3sCIQx8\nD/4Q1JgEqmykQrh/3eU4rNtpzlfnOsY9h4+gxqTUHQ8K9GXZ8df0Pm/ZRQfH+LwQ1JgcqmsMqclx\n1vaY7PNY5uI2fAQ1Jot72PPWR4/wKR5HsX2mdAGgqU8vm0VQYxaosuep6/1qjpdmhro3PaeQlliU\nAzPEyRdZbS/ixr5vHNoSmHQga45FOQCgBl+TdQxxq6VpGJYNXfO5fQR0/6ioMTtU1OiiSTDN4Vgj\nqNujogaAHhSF71zDaq6fe2gENQB0NIfKGeMhqAEgAFX3uENbu51qejgENQDU4GNikLrh1vR5sfU4\nRzMENQAUyIdg2zAcItSGrLYJ6eER1ABQoOu46rH0NR48lM83RwQ1AHQUQoiFOI84/CCoMTmH9xy5\n7LHn9u/WhX07RtgaTN3YIT1E0I79GeeOoMakPbd/t5fXqTpRUZXM09ABRqexeSKoMTkPnr2rsKpu\nq87JipW65mXqATb1zxcbphDF5ORDOq2q6zZ9+zhJEdbTN6VqmmAeB1OIAi1wwkJdXZfQrPv6feA4\njwtBjcnbeeact3vVdfmYHANxKFuhytfrtUEQTwtN35i0tBm8bvO37xMcYY28Pi7iCOY40fQNaNGx\nLG/I8KSyRh7HA5qiosZsVZ0wfVYpBw8dZxw3ekNFHScqaqDCkAsarJ9a2/ozFRWAJghqzF7duY59\nvh8VNoC6CGogZ6hmxJXNi4Q1gErcowYCQHM4uuI+dXzq3qMmqIGAENjwjQAPF0ENRIzARl8I7nAQ\n1IBnRR3A+j7pEdjoG8E9HoZnAR4cPHRcUjKjWUHHr3yQrmxe3Pb37LCsNtLXo9MZ+pIe412PVfSH\noAYqNAnJ/HO7LtyQzlNOYKMvHFPho+kbKLCtkq6pzXCrsgBP3z/rwr4dBDZ6lx5jVNj94x410EGb\ne8NtJzLJhjUToSAkK5sXCewecY8aGFAatulJrU7QF1XThPR8pS0y+X4OqTGOjQv7dvS+7jaqEdRA\ngaarXjU9mZU9j9W2pq8siHeeOaf1U5ev9pYu1brzjLatq85F3XwQ1EBgmAt82tLvdWXzonaeOSep\neDnW1LKfNe1LUXaRkKKZO0wvGXsDgBC1rWp9NRNywkQd66fWtH5qTSubF7f9V6To8Z1nziWV/BrH\nXMCoqAHPfDVd0ww+fc/t371VVXeRhuzhPUe2Deer8/7p7y2r3DEuKmqgQAgzjhHS05UGqY+Qznrw\n7F1b1XFaLaf/lfG9DfCPihoo4aOiTQO/6esQ0tPVV0jnLauQ0w5qdZ6L8RHUQE+yVXk+9BnqMi9D\nhXNdBHNcCGqgJweOnigNZMamTl/2PnFVQGcrXEIUedyjBgocOHrCS/Nz9nWKQpkm7mlqEtJ5+WZp\ngKAGcvoIz6bBT6Udp/zwqLZN3YQ1smj6xmy0WWhjTAzPiouPgM5iyBRSVNSYjdAmdDh98hiV8wT4\nqqKLUFlDYvUszFAIVWqTgA5he1Gsr4Behip7OuqunkVFjagd3nNEh/cc0cFDxwvXcC4SQhVL+Mav\n7uxfQFfco0b0sisKHTx0fGkTd0gBuWz4FsLWZyVNxYy8yqA2s3dL+g5JTzvn9iePXSHpfZL2Snpc\n0o3OuQvJz26XdIukFyUdcc59qJctBzKyHcSKhkN1XY0qH6jZ6r3L6xLWcSGgMYY6FfWvSvoFSb+W\neew2SR92zt1pZrclf18zszdIuknSdZKukvSQmb3eOfei380Gqm2rnj2GtLT9JJ2uExxLb3K0Q0hj\nLJVB7Zz7iJntzT18g6SDyZ/vkbQuaS15/D7n3POSPm9mj0l6o6SP+tlcYLv8CW6opu2iEytrSE9X\n3/ej097dBDaKtL1Hvcs5dz7585OSdiV/3i3pY5nnnU0eA6LUpFk6vTce0n1wdDdGz24gq3Ovb7cY\n39V4jJeZ3WpmG2a28cwzz3TdDMC7oe4dE+xh6nN8dBGqaZRpW1E/ZWZXOufOm9mVkp5OHj8n6erM\n8/Ykj13GOXe3pLulxTjqltsBbONzNq8DR09oZfNi44lSVjYvNm4Cz29zGhA0pY+DgEZI2gb1A5Ju\nlnRn8v8PZh6/18zepUVnsmslfaLrRgJNtAnKZaqGfOWf6+O9CehxFN2Lzg7/y+oa4AQ06qqcmczM\n3qtFx7FXS3pK0jsk/Zak+yVdI+kJLYZnPZs8/6ck/YCkFyS93Tn3YNVGMDMZ+uCjss5XVstOrr47\nk/m+4EC5IScvqTqOMB91Zyar0+v7e0t+9OaS598h6Y6q1wVicGHfjkYncV/hmr0/zj1s/8acVYyQ\nRlPMTAZUuBS8ywcw+Oj1XdSBjdD2a8yQDm1hGMSBub6BmsY+yRLS3S0L6Z1nzvXacWzs4wfxoqIG\nMHksoIGYEdSYLF9DteqOpy7rTJYNCTqHDSsf0GNOWJLOD09ljaZo+gYqHDh6ojLw6/T4JqTHVTek\nmX0MoSGoMWk+ZxcrCuuDh47rwNETS0P4wr4djXuPF2GVrWmou246kKocRz0ExlGjb7GMca4bxnQs\nq6ftDGNlk5z0gabw7epcyExln3kbRw1Mwfqpta0TADOHISRTvHedDdtln6tp68KU9lETVNSYhcN7\njmxVSUOGbNsObVWVNRV1fW2q6i4VdVGYhNbc3SXwyj5L+pp1Q7ru601Z3YqaoMYspOv99h3WywK2\nabj6fK05K+obUBXYbYO6yZzwY/ERgEXb3/R1868xh2DOo+kbKLDzzLlB7z9m+VwRy+cqYVPnoyPf\n2MpCbKym4+zrpNtQp5JuWm1jgYoas5BW1FJ/VXXbjmDLQpfOZf7UHVPt40KuKoS6VNQE3HTQ9A3k\n5MPaZ1A3HTrlsxm87WvOUZ1mcJ8tLn0ENkE9HTR9AznpqkVpYPex0hXCln7fQzWFLwvi9VNrhC5q\noaLGbPkartU2qNtWwPQI96NNJ7MuWN4SeTR9AzX4XpJyqG2gKdyPoeYCJ6RRpG5QM4UoZivmIKsz\n//jpk8dolq+QTu8KhIygBkbWd6AS2NXSsH5u/+5ehu8d3nNkW2dGoAk6k2G20vCKubKGP9sr60VY\ns5IWQkBFjdnr0hls7JBvsg1U1fVlK2yfqKrRBp3JAIXRqazrdqQdo+oM+Rn7AiMmbVfgqkIHMzCO\nGmgglCk5u0wzmv7OgaMntLJ5kTG6nmS/iwv7Xrv1HbUN7b4Cuqxa54IgflTUQIE6Par74mvtbIZw\n9a9JaHcNzK7N5gR2eKiogUillfAQIRpKS0Ks0gU/0nvZ2cAeO5jbvl5+Br+6uBDoDxU1ELAh750T\n2O2U3cPuElxT6XRGeC9HRQ3MUJcmeYartVO0jCYBtZC/4GC/tENFDQSu6T3rMac2nbP8/erQmr5D\nQmAvMIUoMFO+ApYZzZrpa+z1FDFTWzMENTCC9ERV52TVZpiVz2qYsK6PecPro6quj6AGRrYsrA8e\nOr50TeNlaLoeDyFU7sGzd7F/GuIeNTCibEhnm0zXT63p8J4jlzWjtqnYfFTEhH59XfZ3nRaWoqb1\nseYkJ3C7YT1qIALpibnqvmbXJlXCehhd93NZUDe97913cBPQfjA8C4hAesJLm7dDvsfJ5CjV8vun\nyfzry2SDt05o55/DKmBxo6JGFLLNwFOew7rPICyr9tL3XFYN+prWdM7qDN+q2xN6WVjvPHOuMsy7\nBjcVtR9U1Ijetmbh/bu3LToh0RtZunwfLAv6LhcBhHR3l/bhIkTb3o8ufs3sY6+VpMsmYSl63aaB\nTUCPg6BG0J7LBHQegX25Ns3Ty/Zf/rWaXBigWD6w82G5LJybXDAVPbcovKuq86zsxQWhPRyavhGF\nuoEQe2j3GbJdf7/re6HYssq375aMZe+dqlt1E9zN0esbk9Q1fGJRN/iWrTvdJDyL9ll+/PayvgHc\nw+6mKDDH2J9VwU1TuV8ENSZtWQjNKail5WFd97WYHzwM+aAc++JnWXAT2t0R1JiFrlVjyLo0g/ts\nQm+KsO4uDcixg3qZsuU9yxDUl6PXN2ZhykszNu0YNsV9MFchB3Rq+zZWj9s+vOcIYd0SQY3ZOHD0\nRHRVdV+TjMS2HxCXsglXWJ+6HRblwCRMOXh8frYhlq5keUzkPbd/d+EwMJa6rId71JiUKQ/j6lpZ\nt/nM2fes8/vpDHIxNN2iWp3hW03lm8XnXFXXvUdNRY1ZOnD0xKzu6XYN6fTvTYaNIW58h+HgHjXQ\ngyaThYSoairSss+TrY4YWx0v3yFd1LlszpV0UzR9Y1LaVMl9hKjPsctjN3l3eW3COm5tArvJ+Orn\n9u9eukjJ1DE8Cwhc0f1f383xbS9C6vY2r3rO6WQiljndZpiS9CKrycQn2cCls5gfBDVmL4TFPUIc\nguVzaNiUx7tPTZ1pRH1Vv9mQZ5x1OYIaSEwpRHxddPQ1jhthKgrpbJg2DdIHz95FVe0Bvb6BCfLZ\nq31l82Kre5VFvxNbp7q5y94/blvt1v1dqulyVNSIno/OVlSN7bQJXprBw3Vh345tF1jpJCVpVdwl\nTNPfZU3r5uj1jaj12SO6zyBJ3zddSrLPntFdq9iy/cCKW9PluwkcxVg9C5PX5QTfJmTy6zNnNQna\n7HsPOXyp6Weus22sujU9TdakLqqS8z9DOYIak9DXybtrwJSFduhjhpt87roXEVXrYTdFYI+rSVDX\nQWCXI6gxGaGMLc7KBnXRWFOfgV0UhENOgjJ0ZS0R1iEoC2yC2h+CGpPi+8TtqwrMh1h6cvMV1H3e\nQ28brsvet89e3YT3sAjq/jEzGbDEhX07tqriNoGd3pPbKUm6fPk+H/oeytTH+tx9Th6TfU1COw6E\ntB9U1IhCHxV1qmlQpyGdDl1Jq+eicPKx3UP3TPc5B/mQ86gzzK67tlOFoh0qaqCEj/DIh7TUXzVZ\n9Lp9BlLdSrtOMA49PStjtLtJx1H7nCYU3RHUmLQ+AuK5/buX3oPuKySGDJ86Adtke3xV30ULmcCv\nC/t2aP0UIR0Smr4RtLbhNMT93bkaomptcxuBe9h++R52h8vR6xvRC2Vt6awhZhIDQpDeryas+8M9\nakQttJAuGjcNTFl+3m+Mh6BGUEIJ6MN7jmx1GJMIZ8xTOoyRqnpcLHOJqPUR0gcPHd/qMJb+B8zV\nhX07uOc/MipqBCOUk4GPYC5qMhwj8Olg5Q9DvzAWghpBmNLJL5SQzsu3Pkxpnw9h7P3V11zydQ09\nJh6X9BbUZnZI0klJL5X0X51zd/b1XghH0cmsj5m1YjhZpDM5Pbd/t1Y2L452ci3aVwcPHZcCuHhA\nO77nlG+ij6lnsVwvQW1mL5X0i5L+kaSzkj5pZg845x7t4/0QhrLQ7RrGY1cyqaZTVGZndjp46HhQ\nYb24iNi99eeqSVwQpjGOKXqCD6+vivqNkh5zzn1OkszsPkk3SCKoJ8hnkBZdqYdw9Z6enOpePOR7\nya6fWitdw3oI+bB+8OxdW3OWp8a6kEB7hPQ89BXUuyX9cebvZyV9c0/vhcCEEKxdZAOr7MRUFmpV\nJ7JQwzA/t3MorRgxGLPlJ9TjCX71MjOZmX2PpEPOuR9M/v4WSd/snPuRzHNulXRr8tf9ks5435B5\neLWkL4y9EZFi37XHvmuPfdfe1PbdNzjnXlP1pL4q6nOSrs78fU/y2Bbn3N2S7pYkM9uoM40aLse+\na4991x77rj32XXtz3Xd9TXjySUnXmtlrzWyHpJskPdDTewEAMFm9VNTOuRfM7EckfUiL4Vnvds49\n0sd7AQAwZb2No3bO/a6k36359Lv72o4ZYN+1x75rj33XHvuuvVnuuyCWuQQAAMVYlAMAgICNHtRm\ndsjMNs3sMTO7beztCZ2ZPW5mnzGzT5vZRvLYFWb2e2b2R8n/V8bezhCY2bvN7GkzO5N5rHRfmdnt\nyXG4aWbfNs5Wj69kv73TzM4lx92nzezbMz9jvyXM7Goz+30ze9TMHjGzo8njHHcVluy72R97ozZ9\nJ1ON/qEyU41K+l6mGi1nZo9LWnXOfSHz2M9IetY5d2dysbPinJv9ArJm9vclfVnSrznn9iePFe4r\nM3uDpPdqMaveVZIekvR659yLI23+aEr22zslfdk593O557LfMszsSklXOuc+ZWavkvSwpO+S9FZx\n3C21ZN/dqJkfe2NX1FtTjTrnLkpKpxpFMzdIuif58z1aHNyz55z7iKRncw+X7asbJN3nnHveOfd5\nSY9pcXzOTsl+K8N+y3DOnXfOfSr585ckfVaLmRo57ios2XdlZrPvxg7qoqlGl30xkJykh8zs4WR2\nN0na5Zw7n/z5SUm7xtm0KJTtK47Faj9qZn+QNI2nTbfstxJmtlfSN0n6uDjuGsntO2nmx97YQY3m\n/q5z7npJhyW9LWmm3OIW9zLoyl8D+6qR/yzpdZKul3Re0s+PuzlhM7NXSnq/pLc75/48+zOOu+UK\n9t3sj72xg7pyqlFs55w7l/z/aUkf0KKp56nk/k56n+fp8bYweGX7imNxCefcU865F51z/0/Sf9Gl\nJkb2W46ZfbUWQfMe59xvJg9z3NVQtO849sYPaqYabcDMXpF0spCZvULSt2qxmMkDkm5OnnazpA+O\ns4VRKNtXD0i6ycxeZmavlXStpE+MsH1BSkMm8U91aREd9luGmZmkX5H0WefcuzI/4rirULbvOPZ6\nnJmsDqYabWyXpA8sjmd9laR7nXOnzOyTku43s1skPaFFL8nZM7P3Sjoo6dVmdlbSOyTdqYJ95Zx7\nxMzu12LN9BckvW2KvUfrKNlvB83sei2abB+X9MMS+63AmyS9RdJnzOzTyWM/KY67Osr23ffO/dhj\nZjIAAAI2dtM3AABYgqAGACBgBDUAAAEjqAEACBhBDQBAwAhqAAACRlADABAwghoAgID9fwhJGkmE\nkcfBAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "animate_one_density(JunW04, 1023.3)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# January, vtau = 0.4" ] }, { "cell_type": "code", "execution_count": 52, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/html": [ "" ], "text/plain": [ "" ] }, "execution_count": 52, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAeoAAAHiCAYAAAAnPo9XAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3X+wHedd3/HPNwElahLgukmFLdkooY5mYrUy7Z3QmfSH\nhrQgUagpZTymkHGKwXQakCJouTZMm8x0XCx+RJULdHBLBjPEcTyEEMNgpTjDnUzb/LpOIyI7veAm\ndpGQfyQWkBCPhd2nf5zdq72r3bO/nt19nt33a8Zj6dxzz9mzZ7Wf/T77/DDnnAAAQJheMvYGAACA\ncgQ1AAABI6gBAAgYQQ0AQMAIagAAAkZQAwAQMIIaGImZPW5m/3Ds7QAQNoIaQG22cNzMvpj8d9zM\nLPPzvWb2+2b2FTP73/kLETP752b2hJn9hZn9lpldMfynAOJCUAMTZ2Zf5fHlbpX0XZIOSPqbkr5T\n0g9nfv5eSf9L0l+V9FOSfsPMXpNsx3WSflnSWyTtkvQVSb/kcduASSKogXFdb2Z/YGZ/ZmbvM7OX\npz8wsx8ys8fM7Fkze8DMrkoe32tmLhvAZrZuZj+Y/PmtZvY/zOyEmX1R0js9bu/Nkn7eOXfWOXdO\n0s9Jemvyvq+X9LckvcM595xz7v2S/kDSP0t+9/sk/bZz7iPOuS9L+reSvtvMXuVx+4DJIaiBcd0o\n6ZCk12pRob5VkszsWyT9dPLzKyU9Iem+Bq/7zZI+p0Xlekf+h0kT9J8u+e+akte9TtLpzN9PJ4+l\nP/ucc+5LS36+9bvOuf8j6XlJr2/wuYDZ8dkkBqC5u5xzfyJJZvbbkq5PHv8+Se92zn0q+dntki6Y\n2d6ar/snzrn/lPz5hfwPnXP3Srq3xfa+UtKfZf7+55Jemdynzv8s/fnukt9Nf05FDSxBRQ2M68nM\nn7+iRZhJ0lVaVNGSpKSp+Iu6FHpV/tjL1l3uy5K+JvP3r5X0ZbdY3Sf/s/TnXyr53fzPARQgqIEw\n/Ymkb0j/Ymav0KKD1jlJf5E8/Fcyz//63O8vXRbPzL7PzL685L+ypu9HtOhIljqQPJb+7HW5e875\nn2/9rpl9o6Qdkv5w2bYCc0dQA2F6r6R/YWbXm9nLJP0HSR93zj3unHtGi8D+fjN7qZn9gKRvbPLi\nzrn3OOdeueS//1vyq78m6cfMbLeZ7Zb045J+NXnNP5T0aUnvMLOXm9l3S/obkt6f/O57JH2nmf29\n5MLj30v6zdw9bQA5BDUQIOfcQ1r0in6/pPNaBPFNmaf8kKR/o0Vz+HWS/udAm/bLkn5b0meS/34n\neSx1k6RVSRe06Az3PcmFhZxzj0j6l1oE9tOSXiHpXw203UC0bHFrCQAAhIiKGgCAgBHUAAAEjKAG\nACBgBDUAAAEjqAEACFgQU4i++tWvdnv37h17MwBM3OYfXZoIbt+1+TligGE9/PDDX3DOvabqeUEE\n9d69e7WxsTH2ZgCYsANHT+h1b7r097+UdPrksdG2BzCzJ6qfRdM3gBk4cPREo8eBkBDUACaNMEbs\ngmj6BgDf6gZ0/nk0hyM0BDWAyfBRPR88dFzrp9Y8bA3gB0ENIHg0X2POCGoAQRornKmmERo6kwFA\ngvvTCBEVNYBZKArhbNVOSCNUBDWAVsqapusEXh/N2m2ClnBGDAhqAKXaBOqQ95YJWswBQQ3gMvSy\nBsJBUAOIDpU05oSgBhA0QhlzR1ADCBIBDSwQ1AAuc/rksd7uUxPAQDMENYBCaaDWDWwCGOgHM5MB\nKEVIA+MjqAEUYogWEAaavgFsQ0ADYaGiBrCFkAbCQ1AD6IyAB/pjzrmxt0Grq6tuY2Nj7M0AkNEl\nfOlcBlQzs4edc6uVzyOoASxDYAP9qBvUNH0DWOr0yWMELjAighpALYQ1MA6CGgCAgBHUAGprWlXT\nGxzojqAG0AhN4MCwmJkMQCGqYSAMBDUASQQzECqCGpg5AhoIG0ENzNQQAc39bKA7gjpS+ZMsJ0SE\nhOMR8IegBmaor2qagAb8I6iBiSOUgbgR1AAqEcrAeJjwJELcn0ZdPqppji9gXFTUESk66XISRRlC\nGpgGgjoCBDQAzBdBHQgmnUBIuBAEwkFQj4yARl9OnzzW+PgioIHw0JkMmChCGpgGc86NvQ161dfu\nca9764/Xeu7cTibLTrZz2xeop00rDccSMDwze9g5t1r1vFpN32b2uKQvSXpR0gvOuVUzu0LS+yTt\nlfS4pBudcxeS598u6Zbk+Ueccx9q8RkKVZ2EpnbCST/PkE3k2fea2v6csibHCN8rEI9aFXUS1KvO\nuS9kHvsZSc865+40s9skrTjn1szsDZLeK+mNkq6S9JCk1zvnXix7/SYVdR1TPQn56v3NCX1a+D6B\nOHmtqEvcIOlg8ud7JK1LWksev88597ykz5vZY1qE9kfLXmjftV+vjQaVY9HJZg6dsrKfO/286f+X\nnYDnsG/moMv3SEAD8aob1E7SQ2b2oqRfds7dLWmXc+588vMnJe1K/rxb0scyv3s2eaxXczsR5Xv0\nHjh64rJ9QGei8LS9f8zFFjBfdZu+dzvnzpnZX5P0e5J+VNIDzrmvyzzngnNuxcx+QdLHnHO/njz+\nK5IedM79Ru41b5V0qyRdc801f/uJJ57Y9p5zuxfdlM8T99z3pQ8hBynfLxAmr03fzrlzyf+fNrMP\naNGU/ZSZXemcO29mV0p6Onn6OUlXZ359T/JY/jXvlnS3JK2url52tVBVRRRVkHNAQIcn1JDm+wWm\noXIctZm9wsxelf5Z0rdKOiPpAUk3J0+7WdIHkz8/IOkmM3uZmb1W0rWSPtFm406fPFZ57zXUk2Qf\nfH/WOe27voS6DwlpYDrqVNS7JH3AzNLn3+ucO2Vmn5R0v5ndIukJSTdKknPuETO7X9Kjkl6Q9LZl\nPb7rqFNd558/NV0+Y6hhAv+meOwDc1cZ1M65z0k6UPD4FyW9ueR37pB0R+etyxhjPPGYhvqcdXqN\nh6Dr/uh6YVPU4z4koX9/ANpjru+R1T3pdzkR1wmZWAK7ra7h6iucfQ8tnOr3BeCS6IK6TmUdysxa\ny5qrhwjoKQixem1j2fc49+8YwHJBzPW9urrqNjY2Or1GaHNiD9lU20Ro+2mZOYQ0gPkaYmayoCyr\ntIcKzZgnGAlpW6YS0ADgw2SCuk9lwdFlxqixgjGkQC4yhZAOfR8DiMvkgrrrNJpNTCFUpsZXSMZ0\niwDAtE0uqPNCG1bDSb4/Pvct3xOAUEw+qFNjLAUYwoVBTGK7jQAAQ5hNUA/N19rRUxdq73gACMUs\ngpoxy+OgRQEAupt8UNcJiz4COj9cjIuA7tiHAOZo0kEdQkU353DxNT/7nPchAEw6qOsgBPrj40KJ\n7wfA3M0+qBEmAhoAFiYZ1CE0eaMdAhoAtnvJ2BswtgNHTxDsPWkSuqdPHiOkAaDAJCvqNg4cPUFQ\nZPjqrT7klK4AMEWTDOomvY0J5+V8XsAQ0gDQ3CSDOlW1uhUhXU+X6pqLJQDoZtJBXRUS+Z8TGMv1\nURGzzwFguUkHNcJFQANAPZMNau6HhodwBoDmJhnUVSFddu+ant/tsM8AoD+TDOo6yoYNsYhGM+wn\nAOiXOefG3gatrq66jY2Nzq9Tp5Ju+hpzDaKmtw7mup8AoC0ze9g5t1r1vMnMTOYjpHEJ+wsAwjD5\npu+mgUNAXcI4dAAY3+SDGu302Wu+6LUJfQAoNomgZijWOHx2vMt+h4Q2AFwymXvUGE/dFciaPA8A\nsBB9Rc1JvR9V96eL+PouqKgB4JJog5qAngZCGQCWizKo64Y0IdDesvHkPi6S+G4AoJ4og7oupgT1\no0lIs78BwK8ogzoNg6KqjyZxP9oELiENAP5FGdRSedMsYQEAmJLohmfVHeKDfiy7EOK7AQD/ogtq\nDC8fwLRaAMBwogrqsqknCY5h1K2WqawBwJ9ogrronjQBPZ46+5+wBoDuouxMRkCPh/AFgGEFG9QE\nQhi6fA9cUAFAd8EGdRFO/MMipAFgfNHco5bopDSkrvuZ7woA/Iiqok6V9f6egrJwm8rnAwA0E2xQ\nlwXTFKu0GD9T0fcT4+cAgNAFG9RNhbgAh68VqOp+rmXv18SybVv2mmVzsAMA2ptMUIemKKz6COc6\n29DktVjCEgDCEk1QHzx0XBf27Rh7MyoNHcYAgGmLqtd3lbGbXMd+/zI+tyvUzwgAU2XOubG3Qaur\nq25jY6PWc+sGhc+m4zrvE9p2ddmGqtenJQAAujOzh51zq1XPi6bpO2RDBZuvanZZxzsqZgAIC0Hd\nUZ1ga9Opq837NKnui5atJKQBIDwE9YCqhpD5mA1syN8DAPQvuqBuOxY5FDFucxb3pwFgWNEFdVbf\nTbanTx6LPli7mkow+7j9AABjiDqolwlxprKQzGnfZC+2CGwAsZlsUEt+wnpqHa2GDKi2U5EOtQ0A\nEIPog3qqJ+Jl9+KnUA2O/ZmorAHEIroJT7JiG7/cRsxB0sd+a7M/mMAFQIjqTngS9RSiK5sXS3/m\nawawsSv2sd8/NG32R52Z2NjPAEIVdVDHsEgH/EuD1Xe4hhTWIW0LgHFFf4+6LyGdKKfag72PZuxl\n7xXSd1pH2fZO8VgAUI6gLhDiCf3A0RNa2byo9VNrY2+KF23DpktIxRjWRYqmfwUwXVE3fS+7R819\nx3CNGS6nTx4j3ABEJeqKuo971GOE+9SCI4Ym27Lx8SEM2+ICE0BWtEEde0eirvdnQwg9AsWPIdcy\nBxCfKJu+657YljWNx6SoGZ+Q7C6EfUhIA6hSu6I2s5dK2pB0zjn3HWZ2haT3Sdor6XFJNzrnLiTP\nvV3SLZJelHTEOfchz9tdaQ4nthCaaZsIqff62CHd5P1D2WcAxtGk6fuopM9K+prk77dJ+rBz7k4z\nuy35+5qZvUHSTZKuk3SVpIfM7PXOuRc9brd3Q8/pffDQcW89uMeYjrPL2tdjB88YM5V1GVYGYN5q\nBbWZ7ZH0jyXdIenHkodvkHQw+fM9ktYlrSWP3+ece17S583sMUlvlPRRb1tdIZaT28FDxyXpssAe\nu9qbsqEWCun6HcZyDAPoX92K+j9K+glJr8o8tss5dz7585OSdiV/3i3pY5nnnU0e86LOCbDJmOMx\nQnFl8+K2Huu+tmEq9+SLDPE9VXXWG+pYIaQBZFUGtZl9h6SnnXMPm9nBouc455yZNVrdw8xulXSr\nJF1zzTVNfrWWC/t2bAW2tL1qTSvZsaYg9fm+MZ7Um95bH+NiimF6AEJRp6J+k6R/YmbfLunlkr7G\nzH5d0lNmdqVz7ryZXSnp6eT55yRdnfn9Pclj2zjn7pZ0t7RYPavDZ1gqDcX0xJutOvOVbUymcFJf\nFthza/6fwvcJoB+VQe2cu13S7ZKUVNT/2jn3/Wb2s5JulnRn8v8PJr/ygKR7zexdWnQmu1bSJ3xt\ncNdwTX831mbiKZ7Q5xbKqSl+lwD86zLhyZ2S7jezWyQ9IelGSXLOPWJm90t6VNILkt7mu8e3j0o4\ntko6lJN67KFath/7/lyhfH8A4tMoqJ1z61r07pZz7ouS3lzyvDu06CGOjpqe4Lvcf/cVJtnX6TMA\n89tb9F51P1Od53V5fQBoy5zr7fZwbaurq25jY6PWcw8eOh5dNdxGlylGfd57bxu6faxURZgCmBIz\ne9g5t1r5PII6LEW91FNNgi87PC325mofCGsAoZl0UEvx3WNuq01lOtZ92BgQ2ABCUTeoo109ay5i\nC9cmTeVD3csGgJhFE9SH9xxZ/GG/t0nOJmusxTry77ds/vQmzwWAOYum6fvwniN6LgnpuTR7x4Qm\nZQBopm7TdxTrUWdDGmGiEgaAfkTT9I3tmHYTAOYhiooal5w+eay0mXnM5meavgGgH1TUkQg1CEPd\nLgCYiiiC+sGzd83+PnXdntwMeQKAaYmq6XvnmcVqmbGufOXDgaMntiZ9qbKsmRwAEIcoKmppUVVL\nl3qAx7yW9FBYEQoA4hdNUKfSwK5bVU5R1QXKEE3ehDQADCOaCU+KTH2Bjiormxe188y50SaCIawB\noL1ZzPXN6lAafba27L4nuAHAv6g6kxWhCTyMkAYA9CPKipp5vwEAcxF1RU1IL6xsXpz1kDUAmLIo\nK2pckr1YScOaCxgAmI6oe31zj3S7bFWddrSTht9PdCoDgGqTWuayDM29213Yt6OwmmaGMgCIV9RN\n34spRelUlndh346tKjoday1Jp5PJYqR+q+y685IDAKpFHdTMUlZtcQGzuJg5vOfI1uPZ0C7DrQUA\nGF/UQZ0a835srMoubrL7Mq2I2acAMJ6oO5OVIVgul7+fX3aroKxDWqrptK00fwNAsVlMIYr66oZr\n+ryVzYtbVXd6j7vNeuAHjp4grAGgg6h7fZdhApDu0h7kF/bt0HP7d7cK6RQtHADQ3iSDev3UWmGz\nLdopG/YFAOjfJO9R51HRhYNmcABYmMWEJ3UwdCssXDQBQDOTD2oAAGI2i6CmY1lYDhw9QWUNADXN\nIqgRJsIaAKpNPqjp/Q0AiNksen1LVG8xoEc4gDmh13cO96kBADGaTVCvn1ojrANHqwcAXG42QS1x\nvzoGhDUAbDeroJZoAgcAxGV2QY3wUVUDwCWzDGqq6vHU7dmdTopCaAOYu9kMz8pK5/9mRah+1Anj\nNgHM8C0AU1J3eNZXDbExoUk7lVGt+dWkWgYA1DPLpu8UFZofp08eY18CQE9mWVHDjzbh3KWaTn+X\niwIAc0JQo7Wy4KRpGwD8IajR2dDBnH0/qmsAUzfre9QS91cBAGGbfVCnCGwAQIgI6pyVzYtMiAIA\nCAZBnbN+ak07z5wjsAEAQSCoS+w8c27sTQAAgKAu8uDZu/Tg2bskMS84AGBcDM9a4lJVvVsSc4Nn\npR3vGDMNAP0KOqgP7zmy7e9plTu0nWfO6bn9u0d57xBle8ePFdj00AcwF0EHdd7hPUcGC+v8RYK0\naAanqi5WFJxU2wDQXdBB/eDZuy4LzL6r7KKAxiVjV7Jjvz8ADC3ooJYuD+t8E3RataWdvnaeOVc7\nvNPX3faayZ+zvb5p9g4DIY0xFLUM5W/5cGyiT8EHtVRcWZc1QT+3f/fWc8sCOx/QRa91Yd9r6fGd\nw8kIc7Ls1k3+Z/m/828FPkUR1HWkYbuyeXErgA8eOr5VGRdV5lX3m7M/n3Noh3LSOXD0RDDbgmnr\n2r+CYxU+RRPU+aAt69h1+WOL0K5TRS8z105koZ1saGpEyDgu0Ydoglq61JR98NDx2r9zKWB3FzyG\nvFhONDQ1IhQce+hbVEHdxYV9O2bdfF2lj5PN6ZPHBhuiRaWNIXGcYUjmnBt7G7S6uuo2NjZqPz+t\nqJtWxtmgpqq+XF8nnzHGU3MiRVtlxyvHFHwzs4edc6tVz5vNXN9U09X6CtR0rW9OdIgVxy7GNJug\nzlbQVNPjIbABoJnZBHUW1XW5oZqp+w5rpi+FL1xYYmyVncnM7OWSPiLpZcnzf8M59w4zu0LS+yTt\nlfS4pBudcxeS37ld0i2SXpR0xDn3Id8b3qQqJpibGapjVt+dzehgNn30/scc1Kmon5f0Lc65A5Ku\nl3TIzP6OpNskfdg5d62kDyd/l5m9QdJNkq6TdEjSL5nZS/vY+DoI6bANcWI9cPQEFTYqcZwgVJUV\ntVt0C/9y8tevTv5zkm6QdDB5/B5J65LWksfvc849L+nzZvaYpDdK+qjPDa+LYVntDVlZZ9+vL8wW\nFT+CFHNUaxx1UhE/LOmvS/pF59zHzWyXc+588pQnJe1K/rxb0scyv35W2dlGRsS60u1MpSlcojk8\nRoQz5q5WUDvnXpR0vZl9naQPmNn+3M+dmTUakG1mt0q6VZKuueaa2r938NDxVr22s6thoZ0hQm6o\nSVLqvMfK5kWtn1rrfVtwuTbHQNvjsuq9aInB2BpPeGJm/07SVyT9kKSDzrnzZnalpHXn3L6kI5mc\ncz+dPP9Dkt7pnCtt+q474UmbkM6ebA/vOdJ6rm9cbsiTV0hVFSft/jT9nn19F1Xvy3eOPtSd8KQy\nqM3sNZL+0jn3p2a2U9J/k3Rc0j+Q9EXn3J1mdpukK5xzP2Fm10m6V4v70ldp0dHs2qQqL1QnqJvO\nRpbel85XRIS1X3MNa4mTtw9DVs5VCGsMrW5Q12n6vlLSPcl96pdIut859ztm9lFJ95vZLZKekHSj\nJDnnHjGz+yU9KukFSW9bFtJNdA1pabGwR5NFPbDckPd8h+p0Vhf3u4c35r5mKBjGEs1c33Waveve\nU2w7VziWG/rEFUpgc8Luru53OdRwvjY4DtCUt6bvIVQFdVVIt/kHQlj7F8KJaqzwDuGzx6zJ9zbU\nvvZ1LHFsoMykFuUoC9Mu80anlTdjrKdlrLnE08kyQqnyp2wq09wCdQW/HnWfS86tn1rjfrVHIQ1j\nGfN+dpN711NqZp1zZ6wpfzaML+im72VN3j7/YbQdm41iIZ60Ql0Tu+t2hbCvh9y3Y33eOV+EoD8+\ne31PXlpZE9Z+hFRZj2XM+6hz3/d9GGoiHqBI0EGd3kce6h/IyuZFwtqT0IayZN9/qJnP8p95qOOY\n8PaLgMbYgg7qIaXzgBPW0zfkIiChCO3CqY05T66DeYsiqIcITxbrQEiGWqs7dDGMzY/xogdxiSKo\ni+4h+5oVivHUCEn2eA5tJrYpI6ARsiiC+sDRExJBCk9iC74mgRDbZyvTNQSz+2HZa421CAjQRBRB\n3acu96bzk6VQlZdjXuxivveHz1AaU5uRA2WfL3/sTWnsOuYh6HHUWWX/uLqsGVy32bts9rL8Gtes\nylXP2Ce8EAOrbcAOtcxj9v2G3n9DjEf3sQ1AU5MbR11WvV7Yt6NTtbYsVNP3TAP5wbN3LX2tw3uO\nJH+61DGN0A5PiPd+Q9qWVNm/p6GHnS0bbjbEfiOkMbZogjqtmtP1pIuaqn02ry5bKrNMPsiZnhR9\n6TM8mr52m23pGrB9BTShjBBF0/SdlYa1tHzBjirL/rF3aVJv8j5zE9KJMLTvJaR9M7YQm9cB3ya1\nelaer3vBy5r2fIX0sveZG/ZDsbFW/ArZUPuDfY8YRNP0nRfbvV/mCmYO8Lw57IuqGdFowgaqRRXU\nW5Oe9BzSfQ0lGnq+6RAxTGve6J0NNBfNPer0vnTTSrruP9w+171u875TF8IJdcx9H8LnH5LPfV12\nwTu3fYr4TW54lpQOk2oW1mnPa5/3nNHdHJvB5/Z5s7r0DF/2u3Pep5iPaCrqvLpX6EWTlSwL7bFW\nGZprZS2Ne7IdY78TLgCk+hV1tEEt6bKFOoo0DWqp+3q+2eFj0vIhZHMO6KyxwovmbwBjmUVQLxtP\nXTbtZ6pJVV1H2fvF1jt9TGMEVwgXSqGts8wFBDCMWQR1lUtTeqpWhYswzC2wQ7m9QkADwyKoM9LA\nZtGMeAwdGiFU1infn32sEQ0Alptkr++20jm4l80TDoSi7yFHBDQQl1kEdSoNbBbLCB8Toyy0HcbW\ntUMkgHBEOdd3V+un1rSyebGywxnGN1ST9JRCjJAGpmVWFXXW+qk1HTx0nGZwBK/J2H5m6gKmZ7ZB\njXiMNQlNqIqaw0PqDAfAr1n0+l6mzqQpCAvzr1eb+8UMEINJr0ftG/eq43Lg6IlewzTmkGN9ZWB6\nZh/U66fWtjqXIS5TqHwBoMrsgxooEmtVysULMD2zv0edlY6v5p51c3V7Ig/93l3FGHyxXmQAc8M9\n6hbShTpoBq9vZfNiZTD0fd+073vWMSGkgemhos6hqq5vZfNi5ZKheUMEah9hFeqFAMEMxIuKuqWd\nZ85Joqqu48K+HY0DbIheyaGGKgC0QUVdILvaFpV1PW2q66w+wnWoVajGREUNxIvVszrIrrYlEdZD\nSAPHZxi2XdAiJtn91fViCUCYCOoamA98OZ9hmH0tH6HtcxWu0yePBVlVpzhGgWmi6bsCncuqjTU8\nqs181z62NeSwlmgOB2JB07cnaVMic4KX67OJucnrhl7xAkAb9Pquaf3UGpXKErGMZY5lOwEgRUXd\n0OmTx6iulxizA1fVWszZn8+hoxmAaaCibiFdxCP9D+NrUyVPsbLm4gOYHirqloqGwUzxxN+Gz57W\nIeH7BTAGKmqPphZMiAvHHzBNDM/qCdXXJSEESJNhXm1eIwQh7GcA9dUdnkVQ9yz0k3uoxhjvXPWe\nMXyXhDUQD8ZRB8L3TFtzkd9XBBCAuaKiHgmh3U6TwJ7rPuaiBogDy1wGjqFd7TQJ3yGW1ASAvlFR\nB2Ku1V8XU6qusxdtXSfT4eIEiAMVdWTS6o+TbH11wzf0kJa2hzMtLQCy6EwWoKrpL3HJsqlAY9tn\nF/bt0MrmRe08c07rp+7a9rPYPgsAf2j6jhQn7unqeuFBqwwQB4ZnTVx6Mma9bEiL5vK0IgcwLVTU\nE8PKXtNRVRkv+66pqoHwMTPZzBHY0+F7xjRCHAgDTd8zl67uRWADQNwYnjUD3LecLjoVAtNHRT1x\n66fWtjqcYVqaztIGIE5U1DOwfmqNE/XEUEkD80FFPSPpEB7EhVAG5o2KekbWT62xGMhMEfZAvKio\nZ2bnmXPJn3ZLYqIUAAgd46hxGaqv6aKvAhAOVs9Ca5zMp+vA0RNciAGRqQxqM7vazH7fzB41s0fM\n7Gjy+BVm9ntm9kfJ/1cyv3O7mT1mZptm9m19fgD04/TJY9zLnrA0sNPQZggfEK7Kpm8zu1LSlc65\nT5nZqyQ9LOm7JL1V0rPOuTvN7DZJK865NTN7g6T3SnqjpKskPSTp9c65F8veg6bvMDGr2bzRsgL0\ny9sUos6585LOJ3/+kpl9VoueSDdIOpg87R5J65LWksfvc849L+nzZvaYFqH90eYfA2NjSNd85ZvI\nCW5gHI3uUZvZXknfJOnjknYlIS5JT0ralfx5t6Q/zvzaWaVdjLe/1q1mtmFmG88880zDzcYQ0h7i\nNIFD4v42MJbaw7PM7JWS3i/p7c65PzezrZ8555yZNeo+7py7W9Ld0qLpu8nvYhgPnr1LEvcvsV02\nrKmyERJfF5Jlx3X2XJgufDSEWkFtZl+tRUi/xzn3m8nDT5nZlc6588l97KeTx89Jujrz63uSxxCp\n7AFJRYUBQ9kgAAAWMklEQVSsOscDYY6++T4v5S9Gt/6e3AYcupWxTmcy0+Ie9LPOubdnHv9ZSV/M\ndCa7wjn3E2Z2naR7dakz2YclXUtnsvkgzFGEwEYfxjjfrGxe9FJR+1yP+k2S3iLpM2b26eSxn5R0\np6T7zewWSU9IulGSnHOPmNn9kh6V9IKkty0LaUxP9oRMaAOYirEuNuv0+v7vkqzkx28u+Z07JN3R\nYbswEduajQDAo7ncemFmMgAAAkZQo3dTuKIFEJY5tdQR1BjE6ZPHCGwAXswppCWWucTA0rCe2z80\nXPrOuWBDW3M9b7DMJUYz1390WCCwi4357yK076Ttvgjtc5SpOzyLoEZnh/cckXRpJrMmCGvEclL1\nLdRjf9n30XWb637XQ73P2AhqDCadVi+/eEfdfyys0gUpnpOrL03CqGyYY36fhRr+fYr5uCGoMaiy\nsM4q+gd14OgJVujCNr5mfQrdEGOApxzcMQd0yufMZECl7Im1rEIuO2kQ0sjiePBnyGbsIS8KphDS\nTVBRY1BTvsKHX3OorKv+PQwVSE3/XfZ5AdD0/WJGRY0gcU8NdU29sg7p2K871e/QgTnVgG6KoMao\nyu5bA9I01r5uczwP/VlDC+lYv+u+ENQIDpOioMgUQruOED/bEMOqQvzcoSCoESwCG2WmMstZ7Nvv\nA/ugGkGN4HFfG2UOHD0R7Yk+pu2u2s9Nx4SjGYIa0aHSBoY3l9nCQkRQI1rZf/iE9nzlv3sCIQx8\nD/4Q1JgEqmykQrh/3eU4rNtpzlfnOsY9h4+gxqTUHQ8K9GXZ8df0Pm/ZRQfH+LwQ1JgcqmsMqclx\n1vaY7PNY5uI2fAQ1Jot72PPWR4/wKR5HsX2mdAGgqU8vm0VQYxaosuep6/1qjpdmhro3PaeQlliU\nAzPEyRdZbS/ixr5vHNoSmHQga45FOQCgBl+TdQxxq6VpGJYNXfO5fQR0/6ioMTtU1OiiSTDN4Vgj\nqNujogaAHhSF71zDaq6fe2gENQB0NIfKGeMhqAEgAFX3uENbu51qejgENQDU4GNikLrh1vR5sfU4\nRzMENQAUyIdg2zAcItSGrLYJ6eER1ABQoOu46rH0NR48lM83RwQ1AHQUQoiFOI84/CCoMTmH9xy5\n7LHn9u/WhX07RtgaTN3YIT1E0I79GeeOoMakPbd/t5fXqTpRUZXM09ABRqexeSKoMTkPnr2rsKpu\nq87JipW65mXqATb1zxcbphDF5ORDOq2q6zZ9+zhJEdbTN6VqmmAeB1OIAi1wwkJdXZfQrPv6feA4\njwtBjcnbeeact3vVdfmYHANxKFuhytfrtUEQTwtN35i0tBm8bvO37xMcYY28Pi7iCOY40fQNaNGx\nLG/I8KSyRh7HA5qiosZsVZ0wfVYpBw8dZxw3ekNFHScqaqDCkAsarJ9a2/ozFRWAJghqzF7duY59\nvh8VNoC6CGogZ6hmxJXNi4Q1gErcowYCQHM4uuI+dXzq3qMmqIGAENjwjQAPF0ENRIzARl8I7nAQ\n1IBnRR3A+j7pEdjoG8E9HoZnAR4cPHRcUjKjWUHHr3yQrmxe3Pb37LCsNtLXo9MZ+pIe412PVfSH\noAYqNAnJ/HO7LtyQzlNOYKMvHFPho+kbKLCtkq6pzXCrsgBP3z/rwr4dBDZ6lx5jVNj94x410EGb\ne8NtJzLJhjUToSAkK5sXCewecY8aGFAatulJrU7QF1XThPR8pS0y+X4OqTGOjQv7dvS+7jaqEdRA\ngaarXjU9mZU9j9W2pq8siHeeOaf1U5ev9pYu1brzjLatq85F3XwQ1EBgmAt82tLvdWXzonaeOSep\neDnW1LKfNe1LUXaRkKKZO0wvGXsDgBC1rWp9NRNywkQd66fWtH5qTSubF7f9V6To8Z1nziWV/BrH\nXMCoqAHPfDVd0ww+fc/t371VVXeRhuzhPUe2Deer8/7p7y2r3DEuKmqgQAgzjhHS05UGqY+Qznrw\n7F1b1XFaLaf/lfG9DfCPihoo4aOiTQO/6esQ0tPVV0jnLauQ0w5qdZ6L8RHUQE+yVXk+9BnqMi9D\nhXNdBHNcCGqgJweOnigNZMamTl/2PnFVQGcrXEIUedyjBgocOHrCS/Nz9nWKQpkm7mlqEtJ5+WZp\ngKAGcvoIz6bBT6Udp/zwqLZN3YQ1smj6xmy0WWhjTAzPiouPgM5iyBRSVNSYjdAmdDh98hiV8wT4\nqqKLUFlDYvUszFAIVWqTgA5he1Gsr4Behip7OuqunkVFjagd3nNEh/cc0cFDxwvXcC4SQhVL+Mav\n7uxfQFfco0b0sisKHTx0fGkTd0gBuWz4FsLWZyVNxYy8yqA2s3dL+g5JTzvn9iePXSHpfZL2Snpc\n0o3OuQvJz26XdIukFyUdcc59qJctBzKyHcSKhkN1XY0qH6jZ6r3L6xLWcSGgMYY6FfWvSvoFSb+W\neew2SR92zt1pZrclf18zszdIuknSdZKukvSQmb3eOfei380Gqm2rnj2GtLT9JJ2uExxLb3K0Q0hj\nLJVB7Zz7iJntzT18g6SDyZ/vkbQuaS15/D7n3POSPm9mj0l6o6SP+tlcYLv8CW6opu2iEytrSE9X\n3/ej097dBDaKtL1Hvcs5dz7585OSdiV/3i3pY5nnnU0eA6LUpFk6vTce0n1wdDdGz24gq3Ovb7cY\n39V4jJeZ3WpmG2a28cwzz3TdDMC7oe4dE+xh6nN8dBGqaZRpW1E/ZWZXOufOm9mVkp5OHj8n6erM\n8/Ykj13GOXe3pLulxTjqltsBbONzNq8DR09oZfNi44lSVjYvNm4Cz29zGhA0pY+DgEZI2gb1A5Ju\nlnRn8v8PZh6/18zepUVnsmslfaLrRgJNtAnKZaqGfOWf6+O9CehxFN2Lzg7/y+oa4AQ06qqcmczM\n3qtFx7FXS3pK0jsk/Zak+yVdI+kJLYZnPZs8/6ck/YCkFyS93Tn3YNVGMDMZ+uCjss5XVstOrr47\nk/m+4EC5IScvqTqOMB91Zyar0+v7e0t+9OaS598h6Y6q1wVicGHfjkYncV/hmr0/zj1s/8acVYyQ\nRlPMTAZUuBS8ywcw+Oj1XdSBjdD2a8yQDm1hGMSBub6BmsY+yRLS3S0L6Z1nzvXacWzs4wfxoqIG\nMHksoIGYEdSYLF9DteqOpy7rTJYNCTqHDSsf0GNOWJLOD09ljaZo+gYqHDh6ojLw6/T4JqTHVTek\nmX0MoSGoMWk+ZxcrCuuDh47rwNETS0P4wr4djXuPF2GVrWmou246kKocRz0ExlGjb7GMca4bxnQs\nq6ftDGNlk5z0gabw7epcyExln3kbRw1Mwfqpta0TADOHISRTvHedDdtln6tp68KU9lETVNSYhcN7\njmxVSUOGbNsObVWVNRV1fW2q6i4VdVGYhNbc3SXwyj5L+pp1Q7ru601Z3YqaoMYspOv99h3WywK2\nabj6fK05K+obUBXYbYO6yZzwY/ERgEXb3/R1868xh2DOo+kbKLDzzLlB7z9m+VwRy+cqYVPnoyPf\n2MpCbKym4+zrpNtQp5JuWm1jgYoas5BW1FJ/VXXbjmDLQpfOZf7UHVPt40KuKoS6VNQE3HTQ9A3k\n5MPaZ1A3HTrlsxm87WvOUZ1mcJ8tLn0ENkE9HTR9AznpqkVpYPex0hXCln7fQzWFLwvi9VNrhC5q\noaLGbPkartU2qNtWwPQI96NNJ7MuWN4SeTR9AzX4XpJyqG2gKdyPoeYCJ6RRpG5QM4UoZivmIKsz\n//jpk8dolq+QTu8KhIygBkbWd6AS2NXSsH5u/+5ehu8d3nNkW2dGoAk6k2G20vCKubKGP9sr60VY\ns5IWQkBFjdnr0hls7JBvsg1U1fVlK2yfqKrRBp3JAIXRqazrdqQdo+oM+Rn7AiMmbVfgqkIHMzCO\nGmgglCk5u0wzmv7OgaMntLJ5kTG6nmS/iwv7Xrv1HbUN7b4Cuqxa54IgflTUQIE6Par74mvtbIZw\n9a9JaHcNzK7N5gR2eKiogUillfAQIRpKS0Ks0gU/0nvZ2cAeO5jbvl5+Br+6uBDoDxU1ELAh750T\n2O2U3cPuElxT6XRGeC9HRQ3MUJcmeYartVO0jCYBtZC/4GC/tENFDQSu6T3rMac2nbP8/erQmr5D\nQmAvMIUoMFO+ApYZzZrpa+z1FDFTWzMENTCC9ERV52TVZpiVz2qYsK6PecPro6quj6AGRrYsrA8e\nOr50TeNlaLoeDyFU7sGzd7F/GuIeNTCibEhnm0zXT63p8J4jlzWjtqnYfFTEhH59XfZ3nRaWoqb1\nseYkJ3C7YT1qIALpibnqvmbXJlXCehhd93NZUDe97913cBPQfjA8C4hAesJLm7dDvsfJ5CjV8vun\nyfzry2SDt05o55/DKmBxo6JGFLLNwFOew7rPICyr9tL3XFYN+prWdM7qDN+q2xN6WVjvPHOuMsy7\nBjcVtR9U1Ijetmbh/bu3LToh0RtZunwfLAv6LhcBhHR3l/bhIkTb3o8ufs3sY6+VpMsmYSl63aaB\nTUCPg6BG0J7LBHQegX25Ns3Ty/Zf/rWaXBigWD6w82G5LJybXDAVPbcovKuq86zsxQWhPRyavhGF\nuoEQe2j3GbJdf7/re6HYssq375aMZe+dqlt1E9zN0esbk9Q1fGJRN/iWrTvdJDyL9ll+/PayvgHc\nw+6mKDDH2J9VwU1TuV8ENSZtWQjNKail5WFd97WYHzwM+aAc++JnWXAT2t0R1JiFrlVjyLo0g/ts\nQm+KsO4uDcixg3qZsuU9yxDUl6PXN2ZhykszNu0YNsV9MFchB3Rq+zZWj9s+vOcIYd0SQY3ZOHD0\nRHRVdV+TjMS2HxCXsglXWJ+6HRblwCRMOXh8frYhlq5keUzkPbd/d+EwMJa6rId71JiUKQ/j6lpZ\nt/nM2fes8/vpDHIxNN2iWp3hW03lm8XnXFXXvUdNRY1ZOnD0xKzu6XYN6fTvTYaNIW58h+HgHjXQ\ngyaThYSoairSss+TrY4YWx0v3yFd1LlszpV0UzR9Y1LaVMl9hKjPsctjN3l3eW3COm5tArvJ+Orn\n9u9eukjJ1DE8Cwhc0f1f383xbS9C6vY2r3rO6WQiljndZpiS9CKrycQn2cCls5gfBDVmL4TFPUIc\nguVzaNiUx7tPTZ1pRH1Vv9mQZ5x1OYIaSEwpRHxddPQ1jhthKgrpbJg2DdIHz95FVe0Bvb6BCfLZ\nq31l82Kre5VFvxNbp7q5y94/blvt1v1dqulyVNSIno/OVlSN7bQJXprBw3Vh345tF1jpJCVpVdwl\nTNPfZU3r5uj1jaj12SO6zyBJ3zddSrLPntFdq9iy/cCKW9PluwkcxVg9C5PX5QTfJmTy6zNnNQna\n7HsPOXyp6Weus22sujU9TdakLqqS8z9DOYIak9DXybtrwJSFduhjhpt87roXEVXrYTdFYI+rSVDX\nQWCXI6gxGaGMLc7KBnXRWFOfgV0UhENOgjJ0ZS0R1iEoC2yC2h+CGpPi+8TtqwrMh1h6cvMV1H3e\nQ28brsvet89e3YT3sAjq/jEzGbDEhX07tqriNoGd3pPbKUm6fPk+H/oeytTH+tx9Th6TfU1COw6E\ntB9U1IhCHxV1qmlQpyGdDl1Jq+eicPKx3UP3TPc5B/mQ86gzzK67tlOFoh0qaqCEj/DIh7TUXzVZ\n9Lp9BlLdSrtOMA49PStjtLtJx1H7nCYU3RHUmLQ+AuK5/buX3oPuKySGDJ86Adtke3xV30ULmcCv\nC/t2aP0UIR0Smr4RtLbhNMT93bkaomptcxuBe9h++R52h8vR6xvRC2Vt6awhZhIDQpDeryas+8M9\nakQttJAuGjcNTFl+3m+Mh6BGUEIJ6MN7jmx1GJMIZ8xTOoyRqnpcLHOJqPUR0gcPHd/qMJb+B8zV\nhX07uOc/MipqBCOUk4GPYC5qMhwj8Olg5Q9DvzAWghpBmNLJL5SQzsu3Pkxpnw9h7P3V11zydQ09\nJh6X9BbUZnZI0klJL5X0X51zd/b1XghH0cmsj5m1YjhZpDM5Pbd/t1Y2L452ci3aVwcPHZcCuHhA\nO77nlG+ij6lnsVwvQW1mL5X0i5L+kaSzkj5pZg845x7t4/0QhrLQ7RrGY1cyqaZTVGZndjp46HhQ\nYb24iNi99eeqSVwQpjGOKXqCD6+vivqNkh5zzn1OkszsPkk3SCKoJ8hnkBZdqYdw9Z6enOpePOR7\nya6fWitdw3oI+bB+8OxdW3OWp8a6kEB7hPQ89BXUuyX9cebvZyV9c0/vhcCEEKxdZAOr7MRUFmpV\nJ7JQwzA/t3MorRgxGLPlJ9TjCX71MjOZmX2PpEPOuR9M/v4WSd/snPuRzHNulXRr8tf9ks5435B5\neLWkL4y9EZFi37XHvmuPfdfe1PbdNzjnXlP1pL4q6nOSrs78fU/y2Bbn3N2S7pYkM9uoM40aLse+\na4991x77rj32XXtz3Xd9TXjySUnXmtlrzWyHpJskPdDTewEAMFm9VNTOuRfM7EckfUiL4Vnvds49\n0sd7AQAwZb2No3bO/a6k36359Lv72o4ZYN+1x75rj33XHvuuvVnuuyCWuQQAAMVYlAMAgICNHtRm\ndsjMNs3sMTO7beztCZ2ZPW5mnzGzT5vZRvLYFWb2e2b2R8n/V8bezhCY2bvN7GkzO5N5rHRfmdnt\nyXG4aWbfNs5Wj69kv73TzM4lx92nzezbMz9jvyXM7Goz+30ze9TMHjGzo8njHHcVluy72R97ozZ9\nJ1ON/qEyU41K+l6mGi1nZo9LWnXOfSHz2M9IetY5d2dysbPinJv9ArJm9vclfVnSrznn9iePFe4r\nM3uDpPdqMaveVZIekvR659yLI23+aEr22zslfdk593O557LfMszsSklXOuc+ZWavkvSwpO+S9FZx\n3C21ZN/dqJkfe2NX1FtTjTrnLkpKpxpFMzdIuif58z1aHNyz55z7iKRncw+X7asbJN3nnHveOfd5\nSY9pcXzOTsl+K8N+y3DOnXfOfSr585ckfVaLmRo57ios2XdlZrPvxg7qoqlGl30xkJykh8zs4WR2\nN0na5Zw7n/z5SUm7xtm0KJTtK47Faj9qZn+QNI2nTbfstxJmtlfSN0n6uDjuGsntO2nmx97YQY3m\n/q5z7npJhyW9LWmm3OIW9zLoyl8D+6qR/yzpdZKul3Re0s+PuzlhM7NXSnq/pLc75/48+zOOu+UK\n9t3sj72xg7pyqlFs55w7l/z/aUkf0KKp56nk/k56n+fp8bYweGX7imNxCefcU865F51z/0/Sf9Gl\nJkb2W46ZfbUWQfMe59xvJg9z3NVQtO849sYPaqYabcDMXpF0spCZvULSt2qxmMkDkm5OnnazpA+O\ns4VRKNtXD0i6ycxeZmavlXStpE+MsH1BSkMm8U91aREd9luGmZmkX5H0WefcuzI/4rirULbvOPZ6\nnJmsDqYabWyXpA8sjmd9laR7nXOnzOyTku43s1skPaFFL8nZM7P3Sjoo6dVmdlbSOyTdqYJ95Zx7\nxMzu12LN9BckvW2KvUfrKNlvB83sei2abB+X9MMS+63AmyS9RdJnzOzTyWM/KY67Osr23ffO/dhj\nZjIAAAI2dtM3AABYgqAGACBgBDUAAAEjqAEACBhBDQBAwAhqAAACRlADABAwghoAgID9fwhJGkmE\nkcfBAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "animate_one_density(Jan04, 1023.3)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# January, vtau = -0.4" ] }, { "cell_type": "code", "execution_count": 53, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/html": [ "" ], "text/plain": [ "" ] }, "execution_count": 53, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAeoAAAHiCAYAAAAnPo9XAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3X+wHedd3/HPNwElahLgukmFLdkooY5mYrUy7Z3QmfSH\nhrQgUagpZTymkHGKwXQakCJouTZMm8x0XCx+RJULdHBLBjPEcTyEEMNgpTjDnUzb/LpOIyI7veAm\ndpGQfyQWkBCPhd2nf5zdq72r3bO/nt19nt33a8Zj6dxzz9mzZ7Wf/T77/DDnnAAAQJheMvYGAACA\ncgQ1AAABI6gBAAgYQQ0AQMAIagAAAkZQAwAQMIIaGImZPW5m/3Ds7QAQNoIaQG22cNzMvpj8d9zM\nLPPzvWb2+2b2FTP73/kLETP752b2hJn9hZn9lpldMfynAOJCUAMTZ2Zf5fHlbpX0XZIOSPqbkr5T\n0g9nfv5eSf9L0l+V9FOSfsPMXpNsx3WSflnSWyTtkvQVSb/kcduASSKogXFdb2Z/YGZ/ZmbvM7OX\npz8wsx8ys8fM7Fkze8DMrkoe32tmLhvAZrZuZj+Y/PmtZvY/zOyEmX1R0js9bu/Nkn7eOXfWOXdO\n0s9Jemvyvq+X9LckvcM595xz7v2S/kDSP0t+9/sk/bZz7iPOuS9L+reSvtvMXuVx+4DJIaiBcd0o\n6ZCk12pRob5VkszsWyT9dPLzKyU9Iem+Bq/7zZI+p0Xlekf+h0kT9J8u+e+akte9TtLpzN9PJ4+l\nP/ucc+5LS36+9bvOuf8j6XlJr2/wuYDZ8dkkBqC5u5xzfyJJZvbbkq5PHv8+Se92zn0q+dntki6Y\n2d6ar/snzrn/lPz5hfwPnXP3Srq3xfa+UtKfZf7+55Jemdynzv8s/fnukt9Nf05FDSxBRQ2M68nM\nn7+iRZhJ0lVaVNGSpKSp+Iu6FHpV/tjL1l3uy5K+JvP3r5X0ZbdY3Sf/s/TnXyr53fzPARQgqIEw\n/Ymkb0j/Ymav0KKD1jlJf5E8/Fcyz//63O8vXRbPzL7PzL685L+ypu9HtOhIljqQPJb+7HW5e875\nn2/9rpl9o6Qdkv5w2bYCc0dQA2F6r6R/YWbXm9nLJP0HSR93zj3unHtGi8D+fjN7qZn9gKRvbPLi\nzrn3OOdeueS//1vyq78m6cfMbLeZ7Zb045J+NXnNP5T0aUnvMLOXm9l3S/obkt6f/O57JH2nmf29\n5MLj30v6zdw9bQA5BDUQIOfcQ1r0in6/pPNaBPFNmaf8kKR/o0Vz+HWS/udAm/bLkn5b0meS/34n\neSx1k6RVSRe06Az3PcmFhZxzj0j6l1oE9tOSXiHpXw203UC0bHFrCQAAhIiKGgCAgBHUAAAEjKAG\nACBgBDUAAAEjqAEACFgQU4i++tWvdnv37h17MwBM3OYfXZoIbt+1+TligGE9/PDDX3DOvabqeUEE\n9d69e7WxsTH2ZgCYsANHT+h1b7r097+UdPrksdG2BzCzJ6qfRdM3gBk4cPREo8eBkBDUACaNMEbs\ngmj6BgDf6gZ0/nk0hyM0BDWAyfBRPR88dFzrp9Y8bA3gB0ENIHg0X2POCGoAQRornKmmERo6kwFA\ngvvTCBEVNYBZKArhbNVOSCNUBDWAVsqapusEXh/N2m2ClnBGDAhqAKXaBOqQ95YJWswBQQ3gMvSy\nBsJBUAOIDpU05oSgBhA0QhlzR1ADCBIBDSwQ1AAuc/rksd7uUxPAQDMENYBCaaDWDWwCGOgHM5MB\nKEVIA+MjqAEUYogWEAaavgFsQ0ADYaGiBrCFkAbCQ1AD6IyAB/pjzrmxt0Grq6tuY2Nj7M0AkNEl\nfOlcBlQzs4edc6uVzyOoASxDYAP9qBvUNH0DWOr0yWMELjAighpALYQ1MA6CGgCAgBHUAGprWlXT\nGxzojqAG0AhN4MCwmJkMQCGqYSAMBDUASQQzECqCGpg5AhoIG0ENzNQQAc39bKA7gjpS+ZMsJ0SE\nhOMR8IegBmaor2qagAb8I6iBiSOUgbgR1AAqEcrAeJjwJELcn0ZdPqppji9gXFTUESk66XISRRlC\nGpgGgjoCBDQAzBdBHQgmnUBIuBAEwkFQj4yARl9OnzzW+PgioIHw0JkMmChCGpgGc86NvQ161dfu\nca9764/Xeu7cTibLTrZz2xeop00rDccSMDwze9g5t1r1vFpN32b2uKQvSXpR0gvOuVUzu0LS+yTt\nlfS4pBudcxeS598u6Zbk+Ueccx9q8RkKVZ2EpnbCST/PkE3k2fea2v6csibHCN8rEI9aFXUS1KvO\nuS9kHvsZSc865+40s9skrTjn1szsDZLeK+mNkq6S9JCk1zvnXix7/SYVdR1TPQn56v3NCX1a+D6B\nOHmtqEvcIOlg8ud7JK1LWksev88597ykz5vZY1qE9kfLXmjftV+vjQaVY9HJZg6dsrKfO/286f+X\nnYDnsG/moMv3SEAD8aob1E7SQ2b2oqRfds7dLWmXc+588vMnJe1K/rxb0scyv3s2eaxXczsR5Xv0\nHjh64rJ9QGei8LS9f8zFFjBfdZu+dzvnzpnZX5P0e5J+VNIDzrmvyzzngnNuxcx+QdLHnHO/njz+\nK5IedM79Ru41b5V0qyRdc801f/uJJ57Y9p5zuxfdlM8T99z3pQ8hBynfLxAmr03fzrlzyf+fNrMP\naNGU/ZSZXemcO29mV0p6Onn6OUlXZ359T/JY/jXvlnS3JK2url52tVBVRRRVkHNAQIcn1JDm+wWm\noXIctZm9wsxelf5Z0rdKOiPpAUk3J0+7WdIHkz8/IOkmM3uZmb1W0rWSPtFm406fPFZ57zXUk2Qf\nfH/WOe27voS6DwlpYDrqVNS7JH3AzNLn3+ucO2Vmn5R0v5ndIukJSTdKknPuETO7X9Kjkl6Q9LZl\nPb7rqFNd558/NV0+Y6hhAv+meOwDc1cZ1M65z0k6UPD4FyW9ueR37pB0R+etyxhjPPGYhvqcdXqN\nh6Dr/uh6YVPU4z4koX9/ANpjru+R1T3pdzkR1wmZWAK7ra7h6iucfQ8tnOr3BeCS6IK6TmUdysxa\ny5qrhwjoKQixem1j2fc49+8YwHJBzPW9urrqNjY2Or1GaHNiD9lU20Ro+2mZOYQ0gPkaYmayoCyr\ntIcKzZgnGAlpW6YS0ADgw2SCuk9lwdFlxqixgjGkQC4yhZAOfR8DiMvkgrrrNJpNTCFUpsZXSMZ0\niwDAtE0uqPNCG1bDSb4/Pvct3xOAUEw+qFNjLAUYwoVBTGK7jQAAQ5hNUA/N19rRUxdq73gACMUs\ngpoxy+OgRQEAupt8UNcJiz4COj9cjIuA7tiHAOZo0kEdQkU353DxNT/7nPchAEw6qOsgBPrj40KJ\n7wfA3M0+qBEmAhoAFiYZ1CE0eaMdAhoAtnvJ2BswtgNHTxDsPWkSuqdPHiOkAaDAJCvqNg4cPUFQ\nZPjqrT7klK4AMEWTDOomvY0J5+V8XsAQ0gDQ3CSDOlW1uhUhXU+X6pqLJQDoZtJBXRUS+Z8TGMv1\nURGzzwFguUkHNcJFQANAPZMNau6HhodwBoDmJhnUVSFddu+ant/tsM8AoD+TDOo6yoYNsYhGM+wn\nAOiXOefG3gatrq66jY2Nzq9Tp5Ju+hpzDaKmtw7mup8AoC0ze9g5t1r1vMnMTOYjpHEJ+wsAwjD5\npu+mgUNAXcI4dAAY3+SDGu302Wu+6LUJfQAoNomgZijWOHx2vMt+h4Q2AFwymXvUGE/dFciaPA8A\nsBB9Rc1JvR9V96eL+PouqKgB4JJog5qAngZCGQCWizKo64Y0IdDesvHkPi6S+G4AoJ4og7oupgT1\no0lIs78BwK8ogzoNg6KqjyZxP9oELiENAP5FGdRSedMsYQEAmJLohmfVHeKDfiy7EOK7AQD/ogtq\nDC8fwLRaAMBwogrqsqknCY5h1K2WqawBwJ9ogrronjQBPZ46+5+wBoDuouxMRkCPh/AFgGEFG9QE\nQhi6fA9cUAFAd8EGdRFO/MMipAFgfNHco5bopDSkrvuZ7woA/Iiqok6V9f6egrJwm8rnAwA0E2xQ\nlwXTFKu0GD9T0fcT4+cAgNAFG9RNhbgAh68VqOp+rmXv18SybVv2mmVzsAMA2ptMUIemKKz6COc6\n29DktVjCEgDCEk1QHzx0XBf27Rh7MyoNHcYAgGmLqtd3lbGbXMd+/zI+tyvUzwgAU2XOubG3Qaur\nq25jY6PWc+sGhc+m4zrvE9p2ddmGqtenJQAAujOzh51zq1XPi6bpO2RDBZuvanZZxzsqZgAIC0Hd\nUZ1ga9Opq837NKnui5atJKQBIDwE9YCqhpD5mA1syN8DAPQvuqBuOxY5FDFucxb3pwFgWNEFdVbf\nTbanTx6LPli7mkow+7j9AABjiDqolwlxprKQzGnfZC+2CGwAsZlsUEt+wnpqHa2GDKi2U5EOtQ0A\nEIPog3qqJ+Jl9+KnUA2O/ZmorAHEIroJT7JiG7/cRsxB0sd+a7M/mMAFQIjqTngS9RSiK5sXS3/m\nawawsSv2sd8/NG32R52Z2NjPAEIVdVDHsEgH/EuD1Xe4hhTWIW0LgHFFf4+6LyGdKKfag72PZuxl\n7xXSd1pH2fZO8VgAUI6gLhDiCf3A0RNa2byo9VNrY2+KF23DpktIxRjWRYqmfwUwXVE3fS+7R819\nx3CNGS6nTx4j3ABEJeqKuo971GOE+9SCI4Ym27Lx8SEM2+ICE0BWtEEde0eirvdnQwg9AsWPIdcy\nBxCfKJu+657YljWNx6SoGZ+Q7C6EfUhIA6hSu6I2s5dK2pB0zjn3HWZ2haT3Sdor6XFJNzrnLiTP\nvV3SLZJelHTEOfchz9tdaQ4nthCaaZsIqff62CHd5P1D2WcAxtGk6fuopM9K+prk77dJ+rBz7k4z\nuy35+5qZvUHSTZKuk3SVpIfM7PXOuRc9brd3Q8/pffDQcW89uMeYjrPL2tdjB88YM5V1GVYGYN5q\nBbWZ7ZH0jyXdIenHkodvkHQw+fM9ktYlrSWP3+ece17S583sMUlvlPRRb1tdIZaT28FDxyXpssAe\nu9qbsqEWCun6HcZyDAPoX92K+j9K+glJr8o8tss5dz7585OSdiV/3i3pY5nnnU0e86LOCbDJmOMx\nQnFl8+K2Huu+tmEq9+SLDPE9VXXWG+pYIaQBZFUGtZl9h6SnnXMPm9nBouc455yZNVrdw8xulXSr\nJF1zzTVNfrWWC/t2bAW2tL1qTSvZsaYg9fm+MZ7Um95bH+NiimF6AEJRp6J+k6R/YmbfLunlkr7G\nzH5d0lNmdqVz7ryZXSnp6eT55yRdnfn9Pclj2zjn7pZ0t7RYPavDZ1gqDcX0xJutOvOVbUymcFJf\nFthza/6fwvcJoB+VQe2cu13S7ZKUVNT/2jn3/Wb2s5JulnRn8v8PJr/ygKR7zexdWnQmu1bSJ3xt\ncNdwTX831mbiKZ7Q5xbKqSl+lwD86zLhyZ2S7jezWyQ9IelGSXLOPWJm90t6VNILkt7mu8e3j0o4\ntko6lJN67KFath/7/lyhfH8A4tMoqJ1z61r07pZz7ouS3lzyvDu06CGOjpqe4Lvcf/cVJtnX6TMA\n89tb9F51P1Od53V5fQBoy5zr7fZwbaurq25jY6PWcw8eOh5dNdxGlylGfd57bxu6faxURZgCmBIz\ne9g5t1r5PII6LEW91FNNgi87PC325mofCGsAoZl0UEvx3WNuq01lOtZ92BgQ2ABCUTeoo109ay5i\nC9cmTeVD3csGgJhFE9SH9xxZ/GG/t0nOJmusxTry77ds/vQmzwWAOYum6fvwniN6LgnpuTR7x4Qm\nZQBopm7TdxTrUWdDGmGiEgaAfkTT9I3tmHYTAOYhiooal5w+eay0mXnM5meavgGgH1TUkQg1CEPd\nLgCYiiiC+sGzd83+PnXdntwMeQKAaYmq6XvnmcVqmbGufOXDgaMntiZ9qbKsmRwAEIcoKmppUVVL\nl3qAx7yW9FBYEQoA4hdNUKfSwK5bVU5R1QXKEE3ehDQADCOaCU+KTH2Bjiormxe188y50SaCIawB\noL1ZzPXN6lAafba27L4nuAHAv6g6kxWhCTyMkAYA9CPKipp5vwEAcxF1RU1IL6xsXpz1kDUAmLIo\nK2pckr1YScOaCxgAmI6oe31zj3S7bFWddrSTht9PdCoDgGqTWuayDM29213Yt6OwmmaGMgCIV9RN\n34spRelUlndh346tKjoday1Jp5PJYqR+q+y685IDAKpFHdTMUlZtcQGzuJg5vOfI1uPZ0C7DrQUA\nGF/UQZ0a835srMoubrL7Mq2I2acAMJ6oO5OVIVgul7+fX3aroKxDWqrptK00fwNAsVlMIYr66oZr\n+ryVzYtbVXd6j7vNeuAHjp4grAGgg6h7fZdhApDu0h7kF/bt0HP7d7cK6RQtHADQ3iSDev3UWmGz\nLdopG/YFAOjfJO9R51HRhYNmcABYmMWEJ3UwdCssXDQBQDOTD2oAAGI2i6CmY1lYDhw9QWUNADXN\nIqgRJsIaAKpNPqjp/Q0AiNksen1LVG8xoEc4gDmh13cO96kBADGaTVCvn1ojrANHqwcAXG42QS1x\nvzoGhDUAbDeroJZoAgcAxGV2QY3wUVUDwCWzDGqq6vHU7dmdTopCaAOYu9kMz8pK5/9mRah+1Anj\nNgHM8C0AU1J3eNZXDbExoUk7lVGt+dWkWgYA1DPLpu8UFZofp08eY18CQE9mWVHDjzbh3KWaTn+X\niwIAc0JQo7Wy4KRpGwD8IajR2dDBnH0/qmsAUzfre9QS91cBAGGbfVCnCGwAQIgI6pyVzYtMiAIA\nCAZBnbN+ak07z5wjsAEAQSCoS+w8c27sTQAAgKAu8uDZu/Tg2bskMS84AGBcDM9a4lJVvVsSc4Nn\npR3vGDMNAP0KOqgP7zmy7e9plTu0nWfO6bn9u0d57xBle8ePFdj00AcwF0EHdd7hPUcGC+v8RYK0\naAanqi5WFJxU2wDQXdBB/eDZuy4LzL6r7KKAxiVjV7Jjvz8ADC3ooJYuD+t8E3RataWdvnaeOVc7\nvNPX3faayZ+zvb5p9g4DIY0xFLUM5W/5cGyiT8EHtVRcWZc1QT+3f/fWc8sCOx/QRa91Yd9r6fGd\nw8kIc7Ls1k3+Z/m/828FPkUR1HWkYbuyeXErgA8eOr5VGRdV5lX3m7M/n3Noh3LSOXD0RDDbgmnr\n2r+CYxU+RRPU+aAt69h1+WOL0K5TRS8z105koZ1saGpEyDgu0Ydoglq61JR98NDx2r9zKWB3FzyG\nvFhONDQ1IhQce+hbVEHdxYV9O2bdfF2lj5PN6ZPHBhuiRaWNIXGcYUjmnBt7G7S6uuo2NjZqPz+t\nqJtWxtmgpqq+XF8nnzHGU3MiRVtlxyvHFHwzs4edc6tVz5vNXN9U09X6CtR0rW9OdIgVxy7GNJug\nzlbQVNPjIbABoJnZBHUW1XW5oZqp+w5rpi+FL1xYYmyVncnM7OWSPiLpZcnzf8M59w4zu0LS+yTt\nlfS4pBudcxeS37ld0i2SXpR0xDn3Id8b3qQqJpibGapjVt+dzehgNn30/scc1Kmon5f0Lc65A5Ku\nl3TIzP6OpNskfdg5d62kDyd/l5m9QdJNkq6TdEjSL5nZS/vY+DoI6bANcWI9cPQEFTYqcZwgVJUV\ntVt0C/9y8tevTv5zkm6QdDB5/B5J65LWksfvc849L+nzZvaYpDdK+qjPDa+LYVntDVlZZ9+vL8wW\nFT+CFHNUaxx1UhE/LOmvS/pF59zHzWyXc+588pQnJe1K/rxb0scyv35W2dlGRsS60u1MpSlcojk8\nRoQz5q5WUDvnXpR0vZl9naQPmNn+3M+dmTUakG1mt0q6VZKuueaa2r938NDxVr22s6thoZ0hQm6o\nSVLqvMfK5kWtn1rrfVtwuTbHQNvjsuq9aInB2BpPeGJm/07SVyT9kKSDzrnzZnalpHXn3L6kI5mc\ncz+dPP9Dkt7pnCtt+q474UmbkM6ebA/vOdJ6rm9cbsiTV0hVFSft/jT9nn19F1Xvy3eOPtSd8KQy\nqM3sNZL+0jn3p2a2U9J/k3Rc0j+Q9EXn3J1mdpukK5xzP2Fm10m6V4v70ldp0dHs2qQqL1QnqJvO\nRpbel85XRIS1X3MNa4mTtw9DVs5VCGsMrW5Q12n6vlLSPcl96pdIut859ztm9lFJ95vZLZKekHSj\nJDnnHjGz+yU9KukFSW9bFtJNdA1pabGwR5NFPbDckPd8h+p0Vhf3u4c35r5mKBjGEs1c33Waveve\nU2w7VziWG/rEFUpgc8Luru53OdRwvjY4DtCUt6bvIVQFdVVIt/kHQlj7F8KJaqzwDuGzx6zJ9zbU\nvvZ1LHFsoMykFuUoC9Mu80anlTdjrKdlrLnE08kyQqnyp2wq09wCdQW/HnWfS86tn1rjfrVHIQ1j\nGfN+dpN711NqZp1zZ6wpfzaML+im72VN3j7/YbQdm41iIZ60Ql0Tu+t2hbCvh9y3Y33eOV+EoD8+\ne31PXlpZE9Z+hFRZj2XM+6hz3/d9GGoiHqBI0EGd3kce6h/IyuZFwtqT0IayZN9/qJnP8p95qOOY\n8PaLgMbYgg7qIaXzgBPW0zfkIiChCO3CqY05T66DeYsiqIcITxbrQEiGWqs7dDGMzY/xogdxiSKo\ni+4h+5oVivHUCEn2eA5tJrYpI6ARsiiC+sDRExJBCk9iC74mgRDbZyvTNQSz+2HZa421CAjQRBRB\n3acu96bzk6VQlZdjXuxivveHz1AaU5uRA2WfL3/sTWnsOuYh6HHUWWX/uLqsGVy32bts9rL8Gtes\nylXP2Ce8EAOrbcAOtcxj9v2G3n9DjEf3sQ1AU5MbR11WvV7Yt6NTtbYsVNP3TAP5wbN3LX2tw3uO\nJH+61DGN0A5PiPd+Q9qWVNm/p6GHnS0bbjbEfiOkMbZogjqtmtP1pIuaqn02ry5bKrNMPsiZnhR9\n6TM8mr52m23pGrB9BTShjBBF0/SdlYa1tHzBjirL/rF3aVJv8j5zE9KJMLTvJaR9M7YQm9cB3ya1\nelaer3vBy5r2fIX0sveZG/ZDsbFW/ArZUPuDfY8YRNP0nRfbvV/mCmYO8Lw57IuqGdFowgaqRRXU\nW5Oe9BzSfQ0lGnq+6RAxTGve6J0NNBfNPer0vnTTSrruP9w+171u875TF8IJdcx9H8LnH5LPfV12\nwTu3fYr4TW54lpQOk2oW1mnPa5/3nNHdHJvB5/Z5s7r0DF/2u3Pep5iPaCrqvLpX6EWTlSwL7bFW\nGZprZS2Ne7IdY78TLgCk+hV1tEEt6bKFOoo0DWqp+3q+2eFj0vIhZHMO6KyxwovmbwBjmUVQLxtP\nXTbtZ6pJVV1H2fvF1jt9TGMEVwgXSqGts8wFBDCMWQR1lUtTeqpWhYswzC2wQ7m9QkADwyKoM9LA\nZtGMeAwdGiFU1infn32sEQ0Alptkr++20jm4l80TDoSi7yFHBDQQl1kEdSoNbBbLCB8Toyy0HcbW\ntUMkgHBEOdd3V+un1rSyebGywxnGN1ST9JRCjJAGpmVWFXXW+qk1HTx0nGZwBK/J2H5m6gKmZ7ZB\njXiMNQlNqIqaw0PqDAfAr1n0+l6mzqQpCAvzr1eb+8UMEINJr0ftG/eq43Lg6IlewzTmkGN9ZWB6\nZh/U66fWtjqXIS5TqHwBoMrsgxooEmtVysULMD2zv0edlY6v5p51c3V7Ig/93l3FGHyxXmQAc8M9\n6hbShTpoBq9vZfNiZTD0fd+073vWMSGkgemhos6hqq5vZfNi5ZKheUMEah9hFeqFAMEMxIuKuqWd\nZ85Joqqu48K+HY0DbIheyaGGKgC0QUVdILvaFpV1PW2q66w+wnWoVajGREUNxIvVszrIrrYlEdZD\nSAPHZxi2XdAiJtn91fViCUCYCOoamA98OZ9hmH0tH6HtcxWu0yePBVlVpzhGgWmi6bsCncuqjTU8\nqs181z62NeSwlmgOB2JB07cnaVMic4KX67OJucnrhl7xAkAb9Pquaf3UGpXKErGMZY5lOwEgRUXd\n0OmTx6iulxizA1fVWszZn8+hoxmAaaCibiFdxCP9D+NrUyVPsbLm4gOYHirqloqGwUzxxN+Gz57W\nIeH7BTAGKmqPphZMiAvHHzBNDM/qCdXXJSEESJNhXm1eIwQh7GcA9dUdnkVQ9yz0k3uoxhjvXPWe\nMXyXhDUQD8ZRB8L3TFtzkd9XBBCAuaKiHgmh3U6TwJ7rPuaiBogDy1wGjqFd7TQJ3yGW1ASAvlFR\nB2Ku1V8XU6qusxdtXSfT4eIEiAMVdWTS6o+TbH11wzf0kJa2hzMtLQCy6EwWoKrpL3HJsqlAY9tn\nF/bt0MrmRe08c07rp+7a9rPYPgsAf2j6jhQn7unqeuFBqwwQB4ZnTVx6Mma9bEiL5vK0IgcwLVTU\nE8PKXtNRVRkv+66pqoHwMTPZzBHY0+F7xjRCHAgDTd8zl67uRWADQNwYnjUD3LecLjoVAtNHRT1x\n66fWtjqcYVqaztIGIE5U1DOwfmqNE/XEUEkD80FFPSPpEB7EhVAG5o2KekbWT62xGMhMEfZAvKio\nZ2bnmXPJn3ZLYqIUAAgd46hxGaqv6aKvAhAOVs9Ca5zMp+vA0RNciAGRqQxqM7vazH7fzB41s0fM\n7Gjy+BVm9ntm9kfJ/1cyv3O7mT1mZptm9m19fgD04/TJY9zLnrA0sNPQZggfEK7Kpm8zu1LSlc65\nT5nZqyQ9LOm7JL1V0rPOuTvN7DZJK865NTN7g6T3SnqjpKskPSTp9c65F8veg6bvMDGr2bzRsgL0\ny9sUos6585LOJ3/+kpl9VoueSDdIOpg87R5J65LWksfvc849L+nzZvaYFqH90eYfA2NjSNd85ZvI\nCW5gHI3uUZvZXknfJOnjknYlIS5JT0ralfx5t6Q/zvzaWaVdjLe/1q1mtmFmG88880zDzcYQ0h7i\nNIFD4v42MJbaw7PM7JWS3i/p7c65PzezrZ8555yZNeo+7py7W9Ld0qLpu8nvYhgPnr1LEvcvsV02\nrKmyERJfF5Jlx3X2XJgufDSEWkFtZl+tRUi/xzn3m8nDT5nZlc6588l97KeTx89Jujrz63uSxxCp\n7AFJRYUBQ9kgAAAWMklEQVSsOscDYY6++T4v5S9Gt/6e3AYcupWxTmcy0+Ie9LPOubdnHv9ZSV/M\ndCa7wjn3E2Z2naR7dakz2YclXUtnsvkgzFGEwEYfxjjfrGxe9FJR+1yP+k2S3iLpM2b26eSxn5R0\np6T7zewWSU9IulGSnHOPmNn9kh6V9IKkty0LaUxP9oRMaAOYirEuNuv0+v7vkqzkx28u+Z07JN3R\nYbswEduajQDAo7ncemFmMgAAAkZQo3dTuKIFEJY5tdQR1BjE6ZPHCGwAXswppCWWucTA0rCe2z80\nXPrOuWBDW3M9b7DMJUYz1390WCCwi4357yK076Ttvgjtc5SpOzyLoEZnh/cckXRpJrMmCGvEclL1\nLdRjf9n30XWb637XQ73P2AhqDCadVi+/eEfdfyys0gUpnpOrL03CqGyYY36fhRr+fYr5uCGoMaiy\nsM4q+gd14OgJVujCNr5mfQrdEGOApxzcMQd0yufMZECl7Im1rEIuO2kQ0sjiePBnyGbsIS8KphDS\nTVBRY1BTvsKHX3OorKv+PQwVSE3/XfZ5AdD0/WJGRY0gcU8NdU29sg7p2K871e/QgTnVgG6KoMao\nyu5bA9I01r5uczwP/VlDC+lYv+u+ENQIDpOioMgUQruOED/bEMOqQvzcoSCoESwCG2WmMstZ7Nvv\nA/ugGkGN4HFfG2UOHD0R7Yk+pu2u2s9Nx4SjGYIa0aHSBoY3l9nCQkRQI1rZf/iE9nzlv3sCIQx8\nD/4Q1JgEqmykQrh/3eU4rNtpzlfnOsY9h4+gxqTUHQ8K9GXZ8df0Pm/ZRQfH+LwQ1JgcqmsMqclx\n1vaY7PNY5uI2fAQ1Jot72PPWR4/wKR5HsX2mdAGgqU8vm0VQYxaosuep6/1qjpdmhro3PaeQlliU\nAzPEyRdZbS/ixr5vHNoSmHQga45FOQCgBl+TdQxxq6VpGJYNXfO5fQR0/6ioMTtU1OiiSTDN4Vgj\nqNujogaAHhSF71zDaq6fe2gENQB0NIfKGeMhqAEgAFX3uENbu51qejgENQDU4GNikLrh1vR5sfU4\nRzMENQAUyIdg2zAcItSGrLYJ6eER1ABQoOu46rH0NR48lM83RwQ1AHQUQoiFOI84/CCoMTmH9xy5\n7LHn9u/WhX07RtgaTN3YIT1E0I79GeeOoMakPbd/t5fXqTpRUZXM09ABRqexeSKoMTkPnr2rsKpu\nq87JipW65mXqATb1zxcbphDF5ORDOq2q6zZ9+zhJEdbTN6VqmmAeB1OIAi1wwkJdXZfQrPv6feA4\njwtBjcnbeeact3vVdfmYHANxKFuhytfrtUEQTwtN35i0tBm8bvO37xMcYY28Pi7iCOY40fQNaNGx\nLG/I8KSyRh7HA5qiosZsVZ0wfVYpBw8dZxw3ekNFHScqaqDCkAsarJ9a2/ozFRWAJghqzF7duY59\nvh8VNoC6CGogZ6hmxJXNi4Q1gErcowYCQHM4uuI+dXzq3qMmqIGAENjwjQAPF0ENRIzARl8I7nAQ\n1IBnRR3A+j7pEdjoG8E9HoZnAR4cPHRcUjKjWUHHr3yQrmxe3Pb37LCsNtLXo9MZ+pIe412PVfSH\noAYqNAnJ/HO7LtyQzlNOYKMvHFPho+kbKLCtkq6pzXCrsgBP3z/rwr4dBDZ6lx5jVNj94x410EGb\ne8NtJzLJhjUToSAkK5sXCewecY8aGFAatulJrU7QF1XThPR8pS0y+X4OqTGOjQv7dvS+7jaqEdRA\ngaarXjU9mZU9j9W2pq8siHeeOaf1U5ev9pYu1brzjLatq85F3XwQ1EBgmAt82tLvdWXzonaeOSep\neDnW1LKfNe1LUXaRkKKZO0wvGXsDgBC1rWp9NRNywkQd66fWtH5qTSubF7f9V6To8Z1nziWV/BrH\nXMCoqAHPfDVd0ww+fc/t371VVXeRhuzhPUe2Deer8/7p7y2r3DEuKmqgQAgzjhHS05UGqY+Qznrw\n7F1b1XFaLaf/lfG9DfCPihoo4aOiTQO/6esQ0tPVV0jnLauQ0w5qdZ6L8RHUQE+yVXk+9BnqMi9D\nhXNdBHNcCGqgJweOnigNZMamTl/2PnFVQGcrXEIUedyjBgocOHrCS/Nz9nWKQpkm7mlqEtJ5+WZp\ngKAGcvoIz6bBT6Udp/zwqLZN3YQ1smj6xmy0WWhjTAzPiouPgM5iyBRSVNSYjdAmdDh98hiV8wT4\nqqKLUFlDYvUszFAIVWqTgA5he1Gsr4Behip7OuqunkVFjagd3nNEh/cc0cFDxwvXcC4SQhVL+Mav\n7uxfQFfco0b0sisKHTx0fGkTd0gBuWz4FsLWZyVNxYy8yqA2s3dL+g5JTzvn9iePXSHpfZL2Snpc\n0o3OuQvJz26XdIukFyUdcc59qJctBzKyHcSKhkN1XY0qH6jZ6r3L6xLWcSGgMYY6FfWvSvoFSb+W\neew2SR92zt1pZrclf18zszdIuknSdZKukvSQmb3eOfei380Gqm2rnj2GtLT9JJ2uExxLb3K0Q0hj\nLJVB7Zz7iJntzT18g6SDyZ/vkbQuaS15/D7n3POSPm9mj0l6o6SP+tlcYLv8CW6opu2iEytrSE9X\n3/ej097dBDaKtL1Hvcs5dz7585OSdiV/3i3pY5nnnU0eA6LUpFk6vTce0n1wdDdGz24gq3Ovb7cY\n39V4jJeZ3WpmG2a28cwzz3TdDMC7oe4dE+xh6nN8dBGqaZRpW1E/ZWZXOufOm9mVkp5OHj8n6erM\n8/Ykj13GOXe3pLulxTjqltsBbONzNq8DR09oZfNi44lSVjYvNm4Cz29zGhA0pY+DgEZI2gb1A5Ju\nlnRn8v8PZh6/18zepUVnsmslfaLrRgJNtAnKZaqGfOWf6+O9CehxFN2Lzg7/y+oa4AQ06qqcmczM\n3qtFx7FXS3pK0jsk/Zak+yVdI+kJLYZnPZs8/6ck/YCkFyS93Tn3YNVGMDMZ+uCjss5XVstOrr47\nk/m+4EC5IScvqTqOMB91Zyar0+v7e0t+9OaS598h6Y6q1wVicGHfjkYncV/hmr0/zj1s/8acVYyQ\nRlPMTAZUuBS8ywcw+Oj1XdSBjdD2a8yQDm1hGMSBub6BmsY+yRLS3S0L6Z1nzvXacWzs4wfxoqIG\nMHksoIGYEdSYLF9DteqOpy7rTJYNCTqHDSsf0GNOWJLOD09ljaZo+gYqHDh6ojLw6/T4JqTHVTek\nmX0MoSGoMWk+ZxcrCuuDh47rwNETS0P4wr4djXuPF2GVrWmou246kKocRz0ExlGjb7GMca4bxnQs\nq6ftDGNlk5z0gabw7epcyExln3kbRw1Mwfqpta0TADOHISRTvHedDdtln6tp68KU9lETVNSYhcN7\njmxVSUOGbNsObVWVNRV1fW2q6i4VdVGYhNbc3SXwyj5L+pp1Q7ru601Z3YqaoMYspOv99h3WywK2\nabj6fK05K+obUBXYbYO6yZzwY/ERgEXb3/R1868xh2DOo+kbKLDzzLlB7z9m+VwRy+cqYVPnoyPf\n2MpCbKym4+zrpNtQp5JuWm1jgYoas5BW1FJ/VXXbjmDLQpfOZf7UHVPt40KuKoS6VNQE3HTQ9A3k\n5MPaZ1A3HTrlsxm87WvOUZ1mcJ8tLn0ENkE9HTR9AznpqkVpYPex0hXCln7fQzWFLwvi9VNrhC5q\noaLGbPkartU2qNtWwPQI96NNJ7MuWN4SeTR9AzX4XpJyqG2gKdyPoeYCJ6RRpG5QM4UoZivmIKsz\n//jpk8dolq+QTu8KhIygBkbWd6AS2NXSsH5u/+5ehu8d3nNkW2dGoAk6k2G20vCKubKGP9sr60VY\ns5IWQkBFjdnr0hls7JBvsg1U1fVlK2yfqKrRBp3JAIXRqazrdqQdo+oM+Rn7AiMmbVfgqkIHMzCO\nGmgglCk5u0wzmv7OgaMntLJ5kTG6nmS/iwv7Xrv1HbUN7b4Cuqxa54IgflTUQIE6Par74mvtbIZw\n9a9JaHcNzK7N5gR2eKiogUillfAQIRpKS0Ks0gU/0nvZ2cAeO5jbvl5+Br+6uBDoDxU1ELAh750T\n2O2U3cPuElxT6XRGeC9HRQ3MUJcmeYartVO0jCYBtZC/4GC/tENFDQSu6T3rMac2nbP8/erQmr5D\nQmAvMIUoMFO+ApYZzZrpa+z1FDFTWzMENTCC9ERV52TVZpiVz2qYsK6PecPro6quj6AGRrYsrA8e\nOr50TeNlaLoeDyFU7sGzd7F/GuIeNTCibEhnm0zXT63p8J4jlzWjtqnYfFTEhH59XfZ3nRaWoqb1\nseYkJ3C7YT1qIALpibnqvmbXJlXCehhd93NZUDe97913cBPQfjA8C4hAesJLm7dDvsfJ5CjV8vun\nyfzry2SDt05o55/DKmBxo6JGFLLNwFOew7rPICyr9tL3XFYN+prWdM7qDN+q2xN6WVjvPHOuMsy7\nBjcVtR9U1Ijetmbh/bu3LToh0RtZunwfLAv6LhcBhHR3l/bhIkTb3o8ufs3sY6+VpMsmYSl63aaB\nTUCPg6BG0J7LBHQegX25Ns3Ty/Zf/rWaXBigWD6w82G5LJybXDAVPbcovKuq86zsxQWhPRyavhGF\nuoEQe2j3GbJdf7/re6HYssq375aMZe+dqlt1E9zN0esbk9Q1fGJRN/iWrTvdJDyL9ll+/PayvgHc\nw+6mKDDH2J9VwU1TuV8ENSZtWQjNKail5WFd97WYHzwM+aAc++JnWXAT2t0R1JiFrlVjyLo0g/ts\nQm+KsO4uDcixg3qZsuU9yxDUl6PXN2ZhykszNu0YNsV9MFchB3Rq+zZWj9s+vOcIYd0SQY3ZOHD0\nRHRVdV+TjMS2HxCXsglXWJ+6HRblwCRMOXh8frYhlq5keUzkPbd/d+EwMJa6rId71JiUKQ/j6lpZ\nt/nM2fes8/vpDHIxNN2iWp3hW03lm8XnXFXXvUdNRY1ZOnD0xKzu6XYN6fTvTYaNIW58h+HgHjXQ\ngyaThYSoairSss+TrY4YWx0v3yFd1LlszpV0UzR9Y1LaVMl9hKjPsctjN3l3eW3COm5tArvJ+Orn\n9u9eukjJ1DE8Cwhc0f1f383xbS9C6vY2r3rO6WQiljndZpiS9CKrycQn2cCls5gfBDVmL4TFPUIc\nguVzaNiUx7tPTZ1pRH1Vv9mQZ5x1OYIaSEwpRHxddPQ1jhthKgrpbJg2DdIHz95FVe0Bvb6BCfLZ\nq31l82Kre5VFvxNbp7q5y94/blvt1v1dqulyVNSIno/OVlSN7bQJXprBw3Vh345tF1jpJCVpVdwl\nTNPfZU3r5uj1jaj12SO6zyBJ3zddSrLPntFdq9iy/cCKW9PluwkcxVg9C5PX5QTfJmTy6zNnNQna\n7HsPOXyp6Weus22sujU9TdakLqqS8z9DOYIak9DXybtrwJSFduhjhpt87roXEVXrYTdFYI+rSVDX\nQWCXI6gxGaGMLc7KBnXRWFOfgV0UhENOgjJ0ZS0R1iEoC2yC2h+CGpPi+8TtqwrMh1h6cvMV1H3e\nQ28brsvet89e3YT3sAjq/jEzGbDEhX07tqriNoGd3pPbKUm6fPk+H/oeytTH+tx9Th6TfU1COw6E\ntB9U1IhCHxV1qmlQpyGdDl1Jq+eicPKx3UP3TPc5B/mQ86gzzK67tlOFoh0qaqCEj/DIh7TUXzVZ\n9Lp9BlLdSrtOMA49PStjtLtJx1H7nCYU3RHUmLQ+AuK5/buX3oPuKySGDJ86Adtke3xV30ULmcCv\nC/t2aP0UIR0Smr4RtLbhNMT93bkaomptcxuBe9h++R52h8vR6xvRC2Vt6awhZhIDQpDeryas+8M9\nakQttJAuGjcNTFl+3m+Mh6BGUEIJ6MN7jmx1GJMIZ8xTOoyRqnpcLHOJqPUR0gcPHd/qMJb+B8zV\nhX07uOc/MipqBCOUk4GPYC5qMhwj8Olg5Q9DvzAWghpBmNLJL5SQzsu3Pkxpnw9h7P3V11zydQ09\nJh6X9BbUZnZI0klJL5X0X51zd/b1XghH0cmsj5m1YjhZpDM5Pbd/t1Y2L452ci3aVwcPHZcCuHhA\nO77nlG+ij6lnsVwvQW1mL5X0i5L+kaSzkj5pZg845x7t4/0QhrLQ7RrGY1cyqaZTVGZndjp46HhQ\nYb24iNi99eeqSVwQpjGOKXqCD6+vivqNkh5zzn1OkszsPkk3SCKoJ8hnkBZdqYdw9Z6enOpePOR7\nya6fWitdw3oI+bB+8OxdW3OWp8a6kEB7hPQ89BXUuyX9cebvZyV9c0/vhcCEEKxdZAOr7MRUFmpV\nJ7JQwzA/t3MorRgxGLPlJ9TjCX71MjOZmX2PpEPOuR9M/v4WSd/snPuRzHNulXRr8tf9ks5435B5\neLWkL4y9EZFi37XHvmuPfdfe1PbdNzjnXlP1pL4q6nOSrs78fU/y2Bbn3N2S7pYkM9uoM40aLse+\na4991x77rj32XXtz3Xd9TXjySUnXmtlrzWyHpJskPdDTewEAMFm9VNTOuRfM7EckfUiL4Vnvds49\n0sd7AQAwZb2No3bO/a6k36359Lv72o4ZYN+1x75rj33XHvuuvVnuuyCWuQQAAMVYlAMAgICNHtRm\ndsjMNs3sMTO7beztCZ2ZPW5mnzGzT5vZRvLYFWb2e2b2R8n/V8bezhCY2bvN7GkzO5N5rHRfmdnt\nyXG4aWbfNs5Wj69kv73TzM4lx92nzezbMz9jvyXM7Goz+30ze9TMHjGzo8njHHcVluy72R97ozZ9\nJ1ON/qEyU41K+l6mGi1nZo9LWnXOfSHz2M9IetY5d2dysbPinJv9ArJm9vclfVnSrznn9iePFe4r\nM3uDpPdqMaveVZIekvR659yLI23+aEr22zslfdk593O557LfMszsSklXOuc+ZWavkvSwpO+S9FZx\n3C21ZN/dqJkfe2NX1FtTjTrnLkpKpxpFMzdIuif58z1aHNyz55z7iKRncw+X7asbJN3nnHveOfd5\nSY9pcXzOTsl+K8N+y3DOnXfOfSr585ckfVaLmRo57ios2XdlZrPvxg7qoqlGl30xkJykh8zs4WR2\nN0na5Zw7n/z5SUm7xtm0KJTtK47Faj9qZn+QNI2nTbfstxJmtlfSN0n6uDjuGsntO2nmx97YQY3m\n/q5z7npJhyW9LWmm3OIW9zLoyl8D+6qR/yzpdZKul3Re0s+PuzlhM7NXSnq/pLc75/48+zOOu+UK\n9t3sj72xg7pyqlFs55w7l/z/aUkf0KKp56nk/k56n+fp8bYweGX7imNxCefcU865F51z/0/Sf9Gl\nJkb2W46ZfbUWQfMe59xvJg9z3NVQtO849sYPaqYabcDMXpF0spCZvULSt2qxmMkDkm5OnnazpA+O\ns4VRKNtXD0i6ycxeZmavlXStpE+MsH1BSkMm8U91aREd9luGmZmkX5H0WefcuzI/4rirULbvOPZ6\nnJmsDqYabWyXpA8sjmd9laR7nXOnzOyTku43s1skPaFFL8nZM7P3Sjoo6dVmdlbSOyTdqYJ95Zx7\nxMzu12LN9BckvW2KvUfrKNlvB83sei2abB+X9MMS+63AmyS9RdJnzOzTyWM/KY67Osr23ffO/dhj\nZjIAAAI2dtM3AABYgqAGACBgBDUAAAEjqAEACBhBDQBAwAhqAAACRlADABAwghoAgID9fwhJGkmE\nkcfBAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "animate_one_density(JanW04, 1023.3)" ] }, { "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.6.0" } }, "nbformat": 4, "nbformat_minor": 2 }