SHARE
TWEET

Untitled

a guest Jul 19th, 2019 60 Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
  1. {
  2.  "cells": [
  3.   {
  4.    "cell_type": "code",
  5.    "execution_count": 1,
  6.    "metadata": {},
  7.    "outputs": [],
  8.    "source": [
  9.     "import os.path\n",
  10.     "import sys\n",
  11.     "sys.path.append(os.path.abspath('..'))"
  12.    ]
  13.   },
  14.   {
  15.    "cell_type": "code",
  16.    "execution_count": 8,
  17.    "metadata": {},
  18.    "outputs": [],
  19.    "source": [
  20.     "from importlib import reload"
  21.    ]
  22.   },
  23.   {
  24.    "cell_type": "code",
  25.    "execution_count": 48,
  26.    "metadata": {},
  27.    "outputs": [
  28.     {
  29.      "data": {
  30.       "text/plain": [
  31.        "<module 'tpisw.tp_isw_visualisation' from '/home/angusgibson/croco/PYTHON/tpisw/tp_isw_visualisation.py'>"
  32.       ]
  33.      },
  34.      "execution_count": 48,
  35.      "metadata": {},
  36.      "output_type": "execute_result"
  37.     }
  38.    ],
  39.    "source": [
  40.     "reload(iswc)\n",
  41.     "reload(visu)"
  42.    ]
  43.   },
  44.   {
  45.    "cell_type": "code",
  46.    "execution_count": 2,
  47.    "metadata": {},
  48.    "outputs": [],
  49.    "source": [
  50.     "import tpisw\n",
  51.     "from tpisw import tp_isw_class as iswc\n",
  52.     "from tpisw import tp_isw_visualisation as visu"
  53.    ]
  54.   },
  55.   {
  56.    "cell_type": "code",
  57.    "execution_count": 3,
  58.    "metadata": {},
  59.    "outputs": [
  60.     {
  61.      "name": "stdout",
  62.      "output_type": "stream",
  63.      "text": [
  64.       "Computing the Z coordinates with method  2\n",
  65.       "Done\n"
  66.      ]
  67.     }
  68.    ],
  69.    "source": [
  70.     "croco_config = iswc.CROCO_configuration('../../configs/lab_waves_ISW/ISW_his.nc')\n",
  71.     "croco_config.name = 'lab_waves'"
  72.    ]
  73.   },
  74.   {
  75.    "cell_type": "code",
  76.    "execution_count": 43,
  77.    "metadata": {},
  78.    "outputs": [
  79.     {
  80.      "name": "stdout",
  81.      "output_type": "stream",
  82.      "text": [
  83.       "Figure was saved at initial_conditions_lab_waves.png\n"
  84.      ]
  85.     },
  86.     {
  87.      "data": {
  88.       "image/png": "iVBORw0KGgoAAAANSUhEUgAAA/EAAAJQCAYAAADG9uRqAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJzs3XfYLGV9//H3h3KoSkdBVIxBjF3sigqCiFiwxi6oiSUaY2KimPyMGGNLoknswQZWNPaCDfBoQEVRwV4Qj1I99Cad7++Pe5az7Nl92j7nPGc479d17TXPztwzc8/szD77nbulqpAkSZIkSeu+DZY6A5IkSZIkaW4M4iVJkiRJ6gmDeEmSJEmSesIgXpIkSZKknjCIlyRJkiSpJwziJUmSJEnqCYN4SZIkSZJ6wiBe0hqX5OAkleTgea5XSZZPue9du+0cPs12+mDc+UpyaDd/r3lsp/fnLMl+Sb6V5ILuWD6zBHlYnqQWYTsrkqxYhCypp5JsnOTVSX6d5Mrumn7MjeFenY+FfJ9J0o2RQbykibofS1MHITNsf50OTpLcO8l7k/wyySXdj+ffJflEkj9PsuFS53GhFuMByboqya7AZ4HbAO8HXg0cuYRZkqb1UuCfgTOB/6Bd079Y0hytAQt94CtJ65uNljoDktYLnwa+A5w1z/X+DPjj4mdnZkk2Bt4CPB+4FvgG8EXgSmAX4CHA44FPAk9Y2/mbp7fRAtjfz2OdM2jn/qI1kqM1b19gU+ClVfWRpc6MtAgeCVwKPLSqrhrM7L6r+nyvztdCvs8k6UbHIF7SGldVF7GAH5lVtVQlTW8H/hL4MfDEqvrl8MKuBP4pwIFLkLd5qapzgXPnuc7V9LuUb+dueuaS5kJaPDsD5w0H8HCjuFfnZSHfZ5J0Y2R1eknzMtwGs/v7yCTnJrkiyYlJHjlmnRtUkUyyV1dN/9bArQfV9kfbdk5o471zkn9OcnySs5NcleTMJB9J8meLcHz3pwXw5wMPGw3gAarq2qr6EPD0kXU3SPL8JN9LcmmSy7q/X5Bkte/bwfEl2T7JYUnO6qrs/zTJsybkb1mSVyb5TZf2t0n+NckmE9LfoA3p4LPoFj945Nwf2qWZ2M42yU5J3t41hbgqyTlJPpXkHmPSXv+5J9m7O9ZLklyc5IvjPq8kN0vyH10ThsuSXNj9fXiSPxl3jEPrDq6rV3ezvj50bHsNpdstyQeSnDF0/XwgyW4znb8kT01yQvfZrpgpLzPkcVmSFyU5qmuacWWS85McneThs6y7VZK3dfm+IsnPkrw4SRaYlzOTnD5m/u+6Y37lyPwDuvn/MjTvdknekHbvn5NVTU4OS7LLyPpP6dZ/84T8bJLWh8HZSTYas+7Xu+VXJPl5kv837rpP8sAkn09yepefs5N8J8mr5nl+9uu2s7LbzmlJPptk35F0a+y+7677ojUNGf6uXNEtn+levV2ST3bn7LK0PiIekQlV1jNDE5tBPtKaqgzmDf8vuF2Sj3Xn6rqs+r65R5L/TnJyd51fkdau/01JthnZx3Ja8xeA9+eG3027dmkmtolPsk+SLw/t51fdtbnVmLTLu+1slOQfs6qvgdOSvDHJsnHnQZLWFZbES1qoWwPfBU4FPghsCzwJ+GySfavq6zOsu4IWaL2ke/9fQ8tOmmW/DwIOAb5Oq85+KbAbrVr7o5M8oKpOnt+h3MDzuulhVTVj9f+qunJk1geBpwKnAe8BCngs8A5gT+BpYzazNXA8cBXwCVo18CcA70tyXVUdMUiYJMDHaTUAfkOrWroMeDZw5zke30m0c/8q4HfA4UPLls+0YpLbAMfRSgWPBT4K3BJ4IvCIJI+vqi+MWfWRXZ6/BLwLuANwAHCvJHfoStdIsjntXNwW+BrweSC0a+1A2vk5dYYsruiObS/gwcAR3bzBMpLcCzgauAnwOeBnwO1pn82BSfapqhPHbPulwEO7PH0dWC0wmKNtgf8GvkU7xnOAnYBHAUcl+cuqes+Y9ZZ1+d6aVp14Ga1Jx38DuwMvXEBejgWeluT2g1ovSf4UuFW3fB/gNUPpH9JNjxma9zhas5Ovd8d0FXBH4C+ARyW5Z1Wd0aX9NK1GztOSvKyqrhnJz4Hd8b1peFmS99Ku8dOBTwEXAvft8rZPkocO0ifZn9b05WLa53sG7Zz/GfBXrHrAM6Mkr6a1Qb8U+Aztnt4ZuD/t4d3RQ8nX5H3/Gdq1O/pdeeEs+b99t/1taefjR8Cf0D6Do2Y5/Pm6LXAC8Cvgw8BmtPMP7YHoY2lNko4GNgT2AP4OeHiS+1TVJV3aw2nHdSCtT4vh/wWzHe/zgHcClwH/C6ykfQ+8nHYdPqCqxm3jI8ADad9NF9O+l14G7AiMfZAqSeuEqvLly5evsS/aj9EambfrYD7wqpFlD+vmHzUy/+Bu/sEj81cAK2bZ//KReTsCNxmT9q60H9xfmpDfw+d4zL/p0u87z3P1lG69HwBbDs3fAjixW/bUceeX9sN/w6H5dwCuAX42kv6pXfpvA5sOzd92KN+j5+vQbv5es53b2c4Z8JVu/j+NzL9/l9/zRo598LlfA+wzss7ru2UvG5r3qG7ef47J07Jxn/uE/E865gA/75Y9bWTZk7r5vwA2GLOty4C7z/OaWM7q988mwC5j0m4F/IRWA2SzMfdJ0R6gbDLhc3/QfPLWrf/sbt0XDs17Xjfvq7Q+IDYfWvZDWh8Vy4bm3WI4T0Pz96P1J/HOkfn/023/kWPW+WK37M5jrqFPjTkvg8/mb4bmfbKbd9cx299+judlv24bpwK3GLN8l6G/1/h9P3QNrJjHvXpMN/8FI/MfPrT/g8fka9J3wuHd8l3H7LuA101Y79bDxzg0/zndei8fmX/wuLyN+cz3GtnHlbQg/PYj6d/RpT9s3L0JfB/YduRzO6W7dm8+33vKly9fvtbWy+r0khbqd8C/Ds+oqq/QOhy695raaVWtrFUlN8PzT6aVLO6d1tnTQu3UTVerZjyLZ3fTQ6rq0qF8XUYrDYJWOjnqj8DfVdW1Q+v8jFaK9mdJbjKUdlAy9I9VdcVQ+vO5YYnpokurGr0f7fP9t+FlVfUtWqn8trSS2VFHVtUxI/MO66bjrpXLR2dU1VXjPvd5uj+t1P3bVfXhke1/jBYk704rPR11WFX9cMr9U1VXVtVq11a1fiPeB2wD3GvC6q+oodofI5/7QkoNB5/JPkPz9qGVYr6F9uBkT4Ak29EelB1XQ+2yq+qMWr1GClX1VeCntAd7wwYlzAcNz0xy8y7tD6vqx0OL/oYW2D67qkavi9fQHhyNK+kedw3NtS31X3fTl9aqWgTD2xn+/NbGfT8vSW5JqzVxCu2hyXDev8QNaxEshj8woYZDVf1u+BiHvI8WdI9eHwvxdNq1+rZavR+VfwIuAZ4xrukF7SHC+UP5vYxWm2AD4J6LkDdJWiMM4iUt1EkTfpydRgtE1piuXefnu7akVw/aTdJKcjcBtl+E3dQ80+8BXMf4KunfoJXs3H3Msl9X1cVj5p/WTbces4/jxqQft9/FNMj7/1XrTGvUsSPpho2rnj44vuFr5Ru06s+HdG1bX9y1qV2sofz26KbHTlg+0zF8d5HyQJI7du2IT01y+dD1+6YuyS3GrHYNrbr6qOXddFyeZ1RVv6OVNu+d1q47tCrIx9A+i2tYFeDvTavJcINzl+bpaW36z0lyzdDx3Hn0WLoHPr+iVXEe/uyfRqtqffjQtjenPTi4AHhJ1x76+hfwSloJ7HDfCoOHMyckeVeSJ2Wkbf4c3Jd2/395DmnXxn0/X3frpt+uquvGLB/3/TGNk8c9yIHrx7d/UZLjurbq13bXxnXATRl/rc/XxPu6qi6g1SDZlPYAb9Rcv5skaZ1im3hJCzWpjeI1rMEHhEleTGsHfAGtTfHvaaVaBTyG9qN/bCdvc3QWre3oLsBqndrNYCvg/BrpPRqgqq5Jci6tKcComc4jtMBmdB/jguiz55HXhRi0AZ/UT8Bg/rjgY7Vj7M4JDB1fVV2c5L60Ur1Hs6qU7twk7wD+dcKxz9U0x7Ao57c7vmNp/3+PobXbvpgW1NyN1h543PV77oSHZoN8LbSN/jG0dst7AFcDOwDHVNUlSb7HqiB+n6H0w95Ma699Fq25xRmsKgU/mFbVedQRwGuBJ9PaMUMrmb+aVqNjYBvag4MdaH04zKqqPpXWueZLaaXkzwNI8n1aTYavzWEzWwMXjCn5H2dt3PfzNbgW/jBh+aT5CzXTvfExWpv4U2nt3M+mPXiBdt1M8109sOD7usa3k1+Mz0CS1iiDeEm9kdZj9atpPwT3qJGO55LcbxF2cxwtiN+H1QOWmVwEbJtk49FAs8v39qzq7GmhJu4DuPmU257Lvmfaz04j6Rakq6r8nK5U+A60asEvpHUytgGt9HWhpjmG+dbMmOT/0Tr+2ruqlg8vSPIKJg9buH2SDccE8oNjWeh5P5YWxO9L62RtMG8wfUWSbWn3w0W0tt+D/O4IvJjWlv/+o80dkjxlwj4/SKsKfxDwziR3p5Xaf7aqzhlKNzimH1bVHsxRVX0R+GKSLYD70DpWfAHwhSR376qtz+RCYLskm80hkF8b9/18DfZ3swnLJ80vJv8unKlmwNh7I8k9aQH80cABw+cnrdf+l82wzfkYvq9/Omb5onw3SdK6xOr0kpbStcyvtGN72o/Jb40J4LdkVbXKaQzaaj83yaQfu4N9Dpci/ZD2nfqgMUkfRDvOH4xZNh8/6PYxrs32XvPc1nXM79wP2oPvmZHhvzp7d9NpjxFovcFV1U+r6q20XuGh1bSYxuAY9pqwfDB/UY5hgj+lldwuH7PswTOstxGtTf+ovbrpQtvrH0sLwvahPTA5tap+2y07hna9PYM2AsTykYcIf9It/+qYAH6Xbvlqquq0br/3SbI7q9rHHzGS7lJaUHbH7kHCvFTVZVV1bFX9HfA6WrvpGYfx63yHVgNg/zmkXRv3/XwNroX7ZcwQd4z//oBWu+mWozO75ix3Wz35rP60m35uzEPHe9MeZo0aXF8L+W7aa3RBkq1peb+C1qmlJN0oGMRLWkrnATskGfdjbpyVtKrz9+iCdqC1u6RVsZ+6LXxVHQ+8G9gO+HLGjx2+QVfK+MGh2e/rpq/v2vIO0m4OvKF7+94ps/f+bvraJJsO7WNbWgnvfJzHmB/sk3Ql5F+j9Uj9kuFlSe5D6zn/AtoQVguS5E4ZGod6yOBhyh8Xuu3O8bQmEnsmecLIvp9AC7p+xeK3GR62glZye5eR/T+H2Tv5ev3wg6ORz/3941eZWVWtpAXKD6Ad/3Dtk2/Rgp9/7N6Ptjle0U33HO63oLs3383Mtf0O76bPofXwfh4wbnjCN9OC7/d1AdkNJNkmyR5D7/eZ8H0yn2vord30TUlWa7M9Mm9t3Pfz0j0kWU4Lop83vCxtCL59x6wGrd+HWyXZb2T+/2N8s4jZrOime43kYUfg7RPWOa+b3mrC8nE+RGuK8ddpQyQOew2t7f2HJrXbl6Q+sjq9pKV0DK0n7i8n+SatreTJVfX5cYmr6rokb6GNE//jJJ+l/cDfm9Yz+tdZVSI8jRfSSoSeD/w8yXLg5C5/t6CVWO5CG995kLePJDkQ+HPgp0k+w6p2+rcBPj7aI/oCfJQ2FNqjgZ90x78xbXzp79HGa56rY4AnJ/k8bZila4BvVtU3Z1jn+bRA+N+7H/onsmqc+OuAZ03Zg/y+wJuTfIs21NtK2nk+sNv+v0+xbaqqkhxEexjxse78/YLWI/1jaL1YP3NCZ2CL5b9owfpxST5Oq+J7T1rp6Cdon+U4Z9HaD/8kyedY9bnvBLxjls9tNscAdxr6G2g96Sc5ngnt4avq7CRH0tq2n5Tkq7T2yQ+lBf8nMbkE91O0at8v6Y7lreP6O6iq9yW5B22M998kGYyAsS3tvnoQ7QHG87tV3gTs2t2zK2hNBO5Bu2d/Bxw528moqq8meQ2t6cbPu3v5NNqDgD1pJfUHd2nXxn2/EC+k3avvSHIAq8aJfzytbfrgnhr2H7Rr87NJPkYb7vD+tONYzvxr+3yvy8Pjunv6ONo5fDjtYdqZY9b5Nu1By0u6h1SD9vtv7UZwWE1VrUjyEtqDgR9099U5tJot96Pd4y8ft66kG7ckW+7CppfcjE04sS7MUudnUc1lHDpfvnytny9mHif+8AnrLB+zzsGMH5d4C1rHVqfTgsgbbJfx455vBPwd8DNaB1pn00rEb83MYxmPze8sx38fWinar2hj0F9J+zH/adqP9g1G0m9ACzZOpP0Q/SMtQH7haNpJxze0bLVj6eYvo7UPP7XLzwpaJ2GbTDhfhzJ+zPQdgY/QfiRf26U5dLZzRnuI8U5aQHQVcC7wGeBeY9KO/dwnHT+tl/E3d+fvnKHj+wStzfVcP7exxzy0fPfumjmLVoJ3Fq00b/f5bmuWfCxn5F7o5j+SFgheQmt//VVaMDr2fHXnYAUtQH47rfO4K2nVg18MZMr7/FHdfq8DdhxZ9opu2dkT1t28u/5OoQXup3V53G7S8Q+t+x5WjTN+j1ny+EhaSf3K7ro7m1Zy/K8MjQ1Ouy8/Cvyads9eTGuz/1pgh3melwNoPdSfzw3v/YcswX2/gnmME98tuz3tYcmFwGW0APkRwN936zxmzDqP7o7jClqp+JFM8d1Ke9jyji7/VwC/oTVt2HyGY9q/y+ulQ9fHrrPdj7QhML9KqxF0ZXdN/huw9VzvzW7ZwczwveXLl69+vIDNb8GmtTfb1W5sUfdi61rqPC3mK91BSpIk6UYuyYdpzV9uX1XzGYFDknohyaa3ZNPLb8sW7M6WXEdxLOeyI5vwrTr/RlEib5t4SZKkG5Gu347VRmFIsg+tSc7PDOAl3RglWXZrNrt8VzZnd1r3SRsQ9mZ7zuIK9sx2N4oSbEviJUmSbkS6ji8vofUT8gtac6U70voruArYv8aPkCBJvZVk49uw+VU7sQl35qarLb+W4qucw65sxjfqvF6XyBvES5J0I9D13v6SWRM2h1fVijWYHS2hbrSA/2JVJ5yb0/qv+Cbwhqpa6JCEkrROSrLhbdn8mu1Zxt3YamK6ayi+wkp2YwuOrXN7G8gbxEuSdCPQDc/321mSDextSawk6cYgyQa3Y4trt2Ij9mC10UhXczXX8RXO4fZsydF1Ti8DeYN4SZIkSVLvJMnt2fK6zdiAe7PNnNe7muv4Eiu5MzflK7Wyd4G8HdtJkiRJknolSe7ITa5bxgbcaw4l8MM2ZgMexo78iIt5eHbsXam2QbwkSZIkqTeS5C7c9LoA92VrwvwL0zdhA/ZnB37ARTwiN+tVIG8QL0lLIMmhSSrJXkudl7lIsluSTyc5u8v3hd38Xh3HmrIun4cky5Os0z9O+pDH9Y2fiaR12d256XVXcx33Z5sFBfADm7Ah+7Mj3+NCHp2b9+Y7zyBekjSjrqfrzwAHAF8AXg28YUkztZYlObgL0g9e6rysL5Ls1Z3zQ5c6L5Kkdce9s3VdxrU8kG2nCuAHNusC+W9zAY/pSSC/0VJnQJK0zrsNcAfg3VX13JFlbwOOBH6/1nOlG5Nn0oZBkyRpovtmmzqfq3kI2y9KAD+wORuyPzvwZc7hcdmpPlVnrdOd3VkSL0mazc7d9MzRBVV1blX9oqr+uJbzpBuRqvp9Vf1iqfMhSVp3PSDb1kquYm+2Z4NFDOAHtmAj9mMHvsl5PCE7r9Ml8gbxkjSDJNsmeW2SnyT5Y5KLkpyc5A1JthhJu1uSDyQ5I8lVSc7s3u82z33uk+TLSc5PckWSX3X722pM2uVdleNlSf45yS+TXJnk8G75siQvTvKDJBd0x7AiyWeT7DuHvBTwje7tq7p9XV/Feaa24Eme1u338iQrk3wwyc7j2trOVl29W7Z8ZN71+07y1CQnJLk0yYqRdPdJ8omuPf9VSU5L8j9JdmYOuv2+v3v7/qFzUN3Y7KPpn5Dku925Pj/JkUluMWHb2yZ5fZKfd+fpoiTHJNlvLnkb2daTk3x/9HzPss7DkhyV5NzuuvlNkn9Pslo3v911syLJ5l2a33frnJLk5UlW+0WV5NHd8ZzVpT0zyTeS/NVIuhtcE931+/Xu7fB1N/i8n9/9/c8TjuvmSa5O8uM5nLrB9ffJJKd25+/iJMcnefqE9IP7bqMk/5jk193xnZbkjUmWTVhvIff2xmn39m+6dX6R5C+H0j0/yY+7fJ+e5NVJVvt9N99jHLP+/l1+3jdh+SbddXRukk3msk1JmqsHZbs6gyvYh+3ZcA0E8AM36QL55ZzLE9fhQN7q9JI0QZLb0AKJWwPfB95Je/h5O+BvgXcBl3Vp7wUcDdwE+BzwM+D2wNOAA5PsU1UnzmGfz+v2cxnwv8BKYC/g5cCjkjygqi4cs+ongXsBX6K1X1/ZzT8ceArwE+ADwOW0kvU9gf27PM/k1cCuwEG0YH55N3/5+OTXH8c/AP8GXAAcAVwEPBQ4vvt7Mb202/bnaZ/X9QFRkmcB7waupH0upwG7AX9BO5/3rarZmgIcDlwIHAh8FjhpaNnoZ/FXwKO7fX0DuA/wJOCuSe5WVVcO5e3WtPO4K/B/wJeBLYBHAl9O8ryqevdcTkCSvwXe3OXnA930YcC3mHC+uwD41cD5tL4OVgJ3Af4eOCDJ/arq4pHVNga+SruGvgRcAzyG1kfCpt32Btt/LvA/wNm0z+ZcYMduH88C3jHDIX2mm45edwArgBOBNwJ/keS1VXXtyPrPpv3G+Z8Z9jHsnbR79pvAWcB2tD4gPphk96p65YT1PgI8kHYuLu7WeRntOJ81nHCKe/tI2nV0FHA18ATgsCRX087lQbTP7xjatffPwB9p52cxjnHgK8BvgCcl+duqGr2uHt9t803D17kkTWuvbF8r+CP7seMaDeAHbsrG7MsOHM05JLn9ulhTzCBekib7EC2A/8eqev3wgiTbA5d2f4cWON0UeHpVfXgo3ZNoP8I/lOQOVXXdpJ11Qd1buu3ee/ifRpJ3AC+gBcaj7dLp8nmnqjp3aJ2tgCfTHkDcZzTQSbLdbCegqg5NK2U/CFheVYfOtk6SPwFeRwva9qiq07r5h9CCnifPto15eghwv6r64Ug+bkcL4lYAD66qM4aWPQT4GvDfwGNn2nhVHd4VMh8IfKaqDp8h+f7Avarq+hLgJB+hPUg5EPj4UNojaJ/bU6rqyKH0W9OC1rck+VxV/WGm/KXVBngD7YHJHlW1opv/Clqw+Lgx6+xNC7i/DRwwHDym1YZ4f7f8b0dW3Rk4GXhoVV3epX818Cvgb5O8rqqu7tI+D7gKuGtVrRzeSHf/TFRVn0kbAWHidZfkg8ALgYfTgtjB/NAe0vwR+OBM+xlyp6r6zcj2l9GC80OSvGv4+hlyW+COVXV+t84/0c7PM5O8oqrO7uZPc2/fqsvfYESINwG/AP6T9rDmLoO8pdWQOQX4+yRvqqprFuEYAaiqSvIu4N+BZ9D6wxg2yPthk7YhSfO1T3aoU7iMh7EDG62FAH5gazbmIWzPsZz78yS7VdUpa23nc2B1ekkaI8k9gPvTSl1HS7QGbcGv6N7en1bq/u3hAL5L9zHgOGB3Wun3TJ4OLAPeNuap7z8BlwDPmFBV9ZXDAfxg90BopdCrPTyoqvNmyc9CPZX2kPitgwC+218BhwCjpabTOmw0gO+8gFZy/DejwUlVHUsrLX9UkpssYl7eMhzAdwal6fcezEhyV+DBwCeHA/gubxcCr6KVbD9+Dvt8Gu26eesggO+2cx3wD4z57IEXd9O/HC397R5SnNRtd5wXDwL4Lv1KWg2FrWjX+bBraKXHNzDmWl2Id3bT543M34/WGePHxpQWjzUa3HbzrgLeTruW95mw6ssHAXy3zmXAh2m/r+45lG6ae/uQ4c+oqk6lfadsDbxm+Nru0n0e2B64QROOKY5x2PuBKxg550l2p13PX6+qX81hO5I0q/2yQ/2cS3gYO7DxEoSt27KMvdme7dj41xnTfG4pWRIvSePdt5t+ZabS884e3fTYCcuPpQXwd6dVZZ33dqrqgiQ/BB5Ee2Bw8kiS745Z5+IknwceBZyU5JO0atsnrOGO6O7eTY8bk6ffJTmNVoV8sax27J37ddMHd80dRu0IbEhrHvH9RcrLuCYTgwcZ24zJ21YZP4TaDt30z+awz8F1843RBVV1ane+bz2y6H604PqJSZ44ZpvLgB2SbDfysOeiCaUR447xw8CbgJ8m+ViXv+Or6pxZj2gOquqnSb4JPDzJLYceGA1KhN81120luRWtWvs+tJLvzUaSjO3TgLl/3tPc2+P2Mehkctx1OwjqdwF+N5g5xTEO5/W8JB+n1TS4f1V9q1s073MuSTPZPzvWj7mEh7PjkgTwA9uxjAezHd/gvN8mudVw4cRSMoiXpPEGHXtNrF46ZNAG+6wJywfzV+ssbBG3c/aEdZ5E++H+VFa1V74iySeAv5+tqvYCDY5j0rb/wOIG8ZOOfdBc4B9mWX/LRczLuDbNgyrNGw7NG+Ttod1rkrnkbbbzfTarB/Hb0X4DvGqWbW8JDAfx444PxhxjVb05ybm0fgJeDLwEqCTfAP5hLn1EzME7aMHvX9A6wLs5rV34SVU16eHODXTNP75LC7r/j9bm/yJajZFdaVX6x3bUNqEN+7jPe8H39oTaBIN9zLRs48GMaY5xjHfQhgR8HvCtrvbAQbQ2/p+ZaUVJmosDcrM6iYt5ODuwbB2oOL4Dm7An23E85/8+yc5VNem7fK1Z+rMiSeumwY/zWUunWPVD+uYTlu80km7Rt9NVVV9NVV1eVYdW1e1opW9Pp5WQPx34xCz5WahBZ2g3m7B83PxBbYfVHi5nTE/pIyb1Hjs4T1tVVWZ4rVaCvRYM8vY3s+TtWTNu5YbbmnS+x11PFwEXzLLvVNXvxqw7Z1X1gaq6L+2hwSOA99KC7q8k2XGabXc+RXt48ZwkGzL/Du0A/q7L33Oqaq+qenFVvbJrh/+VRcgjLN53xEIt2jFW1QnAD4A/T7INqzq0e39XPV+SFuyRuVl9nwvZnx3Y5AbPQpfWzdmE+7ENO7DszEX6/zUVg3hJGu873fRh44ZrGjFoj73XhOWD+T9Y6HYrmnnkAAAgAElEQVS6QPZutPaoP59lO2NV1Wldm/2HAb8G9pxL53YLMDiO1foA6Dr4uuWYdS7opuOW3XPMvLkYfIYPXOD6wwbt+BfrF8Vi5m1wXT14dEFXAjvunH4H2CbJHRdh/7Oqqgur6qiq+ktab//bMvuxz3rOu0703kN72PYoWon8pbSq/HP1p930k2OWrXZOF2iN3ttzsNjH+E5anw3PpFWlL1b1/SBJC3Jgbl7f5UL2Z0c2XYcC+IGd2ZR7szU7suwPa+j305wZxEvSGFX1fdrwXHejVUe/gSTbJdm0e3s88EtaUPyEkXRPoJU8/ooxbcRHfIjWTvmvk/zpyLLX0Hq//1DNcfimJDskuc+YRVvQhsK7htZ7+GL7SLftv05yfQDZ9Rr+esYHZSfSSuOfmmTzoXW2pfXavRBvo53P/+x6qr+BJMuSzDWIHlQpv9UC83IDXVXy/wMel+TZ49IkufMcn/Z/mFXXza5D629A60l83P/6/+ym786YseSTbJHkvqPz5yNtXPFxzfYGxzRbvwxzPeeH0QL+t9E6tPtIVV0y54y20QtgJMBO8jDaQ4HFsKj39gKs6KZ7Dc+c4hg/Qqs18DLaQ4Cvjes4T5Lm6rHZqb7FBTyMHdhsHQzgB3ZhM+7B1tyMTc6dQ01BAJJsmuS7SU5O8tNuVJep2CZekiZ7Om2or9cleXz3d2jjjO9H64RqRTf00kG0Ics+luSztCGgdqeNoX0J8MzZOsirqhVJXkLrLfoHXQdS59B+JN+v2+ZqDxRmcAvgO0l+TiutPY0WLDySVq33LfMMduakqn6TNgb564CTu07NBuPEb0vruOsuI+ucleTDtKGrTkryxS6vB9A6A7w781RVv+gC5PfROlf7Mu1hysa0wPCBtPN7+zls7tu0oPMl3YOFQfvzt861B/Qxnkrr6Oy9SV4MnEBrxrEL7fzcifa5r5y4Ba6/bg6hdSL3w6Hz/TBaG+sfsfr5PqZb5/XAr5McBfyW1gb+1rRr7jjakHkLdSSt/4XjaEFkaOf8XrQO2Y6eZf1f0vqkeHKSq4Df00p8Pzhczb+qft9dL4/uZs2nKj20Nt7PAv636/zxDNq53582JOCT5rm91ayBe3u+FvUYq+qPSY5g1SgH8z3nknS9x2en+j/OY392ZIsehKe3YjOuoziZiy9IctM5/Ja6EnhIVV2aZGPguCRfqqrvzLLeROv+WZKkJVJVv02yB6206THAi2hVXlfQAqaVQ2lP6HpA/3/AvrSqvecCH6UNA/XLOe7zHUlOAf6e1tZ0c1rw/e/A6yZ0pDXJClrHZXsBe9OGnTqfFhwdQguy1oiqen2S02ltcZ9Fe5DxFdq5/Cqr2s0P+0tacPwU2vjfv6eNrf3vwJ8vMB8fSnIy8FLaOdgPuIzWu/cngI/NcTsXdA9yXtUdzxbdog+xwHbMVXV6N5ThX9M+66fRaimcDfwMeCswOlzdpG29OclZtE78DuaG5/sjE9Z5Y5LjaYHYnrRx7C+iBXiHTVpvHg6hPUjYg/Yw5gpab+kvB95Zq8aTn3RM1yZ5LPAG2ud/E9qDgOMY6nW98z5aEH9iVc3WbGV0Pz9Ksjfwr10+N6I9aHoc7aHK1EF8t5/FvLfnu+81cYzvo107Z9GGa5SkeUty1+3YmP3YgS17FJruyuZcS7EBXJxky26I0bG6fosu7d5u3L0m9eczJ5nQF5IkSYsuyU1pgfpJVXW/2dJLc9EN0/cq4C+q6r1LnJ31QpKDaePG/2tVvXKJsyOpp5Kc8GR2vvdWqwbU6JWfcDHHc8EXuWHHpYdV1WHD6brOV79P66Pk7VU1Ve2r/jzukCT1RpIdgAuHS1u79tFvonWI9emlyptuXJLcBHg+rZbJR5c4O+uF7l7+O1rfF1allzSNjR61/Vazp1pHbX9VcfzFF/yyqh45U7qquha4W9eO/tNJ7lRVP1nofg3iJUlrwuOBf0lyNK3K8La0Dv5uB5xEqyouLViSR9Cq6j+KNrze31fVbJ3laQpJ9qS1498LuDPwtqo6fUkzJUk9UlUXJllO65PEIF6StE45gdZ2+UG0MaShdZz2WuCNVXX5UmVMNxpPBA6iNc94Pat63Neasy+t2cL5tCHlXra02ZGkdV9XO/HqLoDfjPZd+sZptmkQL0ladFX1Q1qnWdIaUVUH0zrx01pSVYcChy5xNiSpb3YCjujaxW8AfLyqvjDNBg3i5yAbbVpZtuVSZ0OSJEmSlkRdft65VbXDUuejb6rqRyxgqNyZGMTPQZZtyUa7P3r2hJIkSZJ0I3T1Se8fHd5TS2SDpc6AJEmSJEmaG4N4SZIkSZJ6wiBekiRJkqSeMIiXJEmSJKknDOIlSZIkSeoJg3hJkiRJknrCIF6SJEmSpJ4wiJckSZIkqScM4iVJkiRJ6gmDeEmSJEmSesIgXpIkSZKknjCIlyRJkiSpJwziJUmSJEnqCYN4SZIkSZJ6wiBekiRJkqSeMIiXJEmSJKknDOIlSZIkSeoJg3hJkiRJknrCIF6SJEmSpJ4wiJckSZIkqScM4iVJkiRJ6gmDeEmSJEmSesIgXpIkSZKknjCIlyRJkiSpJwziJUmSJEnqCYN4SZIkSZJ6wiBekiRJkqSeMIiXJEmSJKknDOIlSZIkSeoJg3hJkiRJknrCIF6SJEmSpJ4wiJckSZIkqScM4iVJkiRJ6gmDeEmSJEmSesIgXpIkSZKknjCIlyRJkiSpJwziJUmSJEnqCYN4SZIkSZJ6wiBekiRJkqSeMIiXJEmSJKknDOIlSZIkSeoJg3hJkiRJknrCIF6SJEmSpJ4wiJckSZIkqScM4iVJkiRJ6gmDeEmSJEmSesIgXpIkSZKknjCIlyRJkiSpJwziJUmSJEnqCYN4SZIkSZJ6wiBekiRJkqSeMIiXJEmSJKknDOIlSZIkSeqJjZY6A32w0WY34WZ3etBSZ0OSJEmSlsTpJ71/qbOgjiXxkiRJkiT1hEG8JEmSJEk9YRAvSZIkSVJPGMRLkiRJktQTBvGSJEmSJPWEQbwkSZIkST1hEC9JkiRJUk8YxEuSJEmS1BMG8ZIkSZIk9YRBvCRJkiRJPWEQL0mSJElSTxjES5IkSZLUEwbxkiRJkiT1hEG8JEmSJEk9YRAvSZIkSVJPGMRLkiRJktQTBvGSJEmSJPWEQbwkSZIkST1hEC9JkiRJUk8YxEuSJEmS1BO9CuKTbJvka0l+3U23mZBu/yS/THJKkkOG5h+a5IwkJ3WvA9Ze7iVJkiRJ65Mkt0zy9SQ/T/LTJH8z7TZ7FcQDhwDHVNVuwDHd+xtIsiHwduDhwB2ApyS5w1CS/6yqu3Wvo9ZGpiVJkiRJ66VrgJdW1Z8B9wVeOBKfzlvfgvgDgSO6v48AHjMmzb2BU6rq1Kq6CjiyW0+SJEmSpLWmqs6qqh90f18C/By4xTTb3GgxMrYW3ayqzoJ2MpLsOCbNLYDTht6fDtxn6P2LkjwTOJH2ROSCNZZbSZIkSdJED37RA5c6Cwu28W/PhCN+f+ckJw7NPqyqDhuXPsmuwN2BE6bZ7zoXxCc5Grj5mEX/NNdNjJlX3fSdwGu6968B3gQ8e0I+ngs8F2DDLbab464lSZIkSeuRH1fVS2dLlGRL4JPAS6rq4ml2uM4F8VW176RlSf6QZKeuFH4nYOWYZKcDtxx6vwtwZrftPwxt693AF2bIx2HAYQDLtrtNTUonSZIkSdIkSTamBfAfrqpPTbu9vrWJ/xxwUPf3QcBnx6T5HrBbktskWQY8uVuPLvAfeCzwkzWYV0mSJEnSeixJgPcCP6+qNy/GNvsWxL8BeGiSXwMP7d6TZOckRwFU1TXAi4Cv0DoN+HhV/bRb/9+S/DjJj4C9gb9d2wcgSZIkSVpvPAB4BvCQxRrqfJ2rTj+TqjoP2GfM/DOBA4beHwWsNnxcVT1jjWZQkiRJkqROVR3H+H7bFqxvJfGSJEmSJK23DOIlSZIkSeoJg3hJkiRJknrCIF6SJEmSpJ4wiJckSZIkqScM4iVJkiRJ6gmDeEmSJEmSesIgXpIkSZKknjCIlyRJkiSpJwziJUmSJEnqCYN4SZIkSZJ6wiBekiRJkqSeMIiXJEmSJKknDOIlSZIkSeoJg3hJkiRJknrCIF6SJEmSpJ4wiJckSZIkqScM4iVJkiRJ6gmDeEmSJEmSesIgXpIkSZKknjCIlyRJkiSpJwziJUmSJEnqCYN4SZIkSZJ6wiBekiRJkqSeMIiXJEmSJKknDOIlSZIkSeoJg3hJkiRJknrCIF6SJEmSpJ4wiJckSZIkqScM4iVJkiRJ6gmDeEmSJEmSesIgXpIkSZKknjCIlyRJkiSpJzZa6gz0wdZbb86Bj9ljqbMhSZIkSUvi7R9a6hxowJJ4SZIkSZJ6wiBekiRJkqSeMIiXJEmSJKknDOIlSZIkSeoJg3hJkiRJknrCIF6SJEmSpJ4wiJckSZIkqScM4iVJkiRJ6gmDeEmSJEmSesIgXpIkSZKknjCIlyRJkiSpJwziJUmSJEnqCYN4SZIkSZJ6wiBekiRJkqSeMIiXJEmSJKknDOIlSZIkSeoJg3hJkiRJknrCIF6SJEmSpJ4wiJckSZIkqScM4iVJkiRJ6gmDeEmSJEmSesIgXpIkSZKknjCIlyRJkiSpJwziJUmSJEnqCYN4SZIkSZJ6wiBekiRJkqSeMIiXJEmSJKknDOIlSZIkSeoJg3hJkiRJknrCIF6SJEmSpJ4wiJckSZIkqScM4iVJkiRJ6gmDeEmSJEmSesIgXpIkSZKknjCIlyRJkiSpJwziJUmSJEnqCYN4SZIkSZJ6wiBekiRJkqSeMIiXJEmSJKknDOIlSZIkSeqJXgXxSbZN8rUkv+6m20xI974kK5P8ZCHrS5IkSZI0rUmx6TR6FcQDhwDHVNVuwDHd+3EOB/afYn1JkiRJkqZ1OONj0wXrWxB/IHBE9/cRwGPGJaqqbwLnL3R9SZIkSZKmNUNsumB9C+JvVlVnAXTTHdfU+kmem+TEJCdefvEFC86wJEmSJOlG686DuLF7PXdN73CjNb2D+UpyNHDzMYv+aW3mo6oOAw4D2PG2d6y1uW9JkiRJWh/c/KAXLHUWFmzb75wIR3zhx1X10rW533UuiK+qfSctS/KHJDtV1VlJdgJWznPz064vSZIkSdKS6Vt1+s8BB3V/HwR8di2vL0mSJEnSkulbEP8G4KFJfg08tHtPkp2THDVIlOSjwLeB3ZOcnuQ5M60vSZIkSdJimyE2XbB1rjr9TKrqPGCfMfPPBA4Yev+U+awvSZIkSdJimxSbTqNvJfGSJEmSJK23DOIlSZIkSeoJg3hJkiRJknrCIF6SJEmSpJ4wiJckSZIkqScM4iVJkiRJ6gmDeEmSJEmSesIgXpIkSZKknjCIlyRJkiSpJwziJUmSJEnqCYN4SZIkSZJ6wiBekiRJkqSe2GipM9AHO265jBftuetSZ0OSJEmSlsTblzoDup4l8ZIkSZIk9YRBvCRJkiRJPWEQL0mSJElSTxjES5IkSZLUEwbxkiRJkiT1hEG8JEmSJEk9YRAvSZIkSVJPGMRLkiRJktQTBvGSJEmSJPWEQbwkSZIkST1hEC9JkiRJUk8YxEuSJEmS1BMG8ZIkSZIk9YRBvCRJkiRJPWEQL0mSJElSTxjES5IkSZLUEwbxkiRJkiT1hEG8JEmSJEk9YRAvSZIkSVJPGMRLkiRJktQTBvGSJEmSJPWEQbwkSZIkST1hEC9JkiRJUk8YxEuSJEmS1BMG8ZIkSZIk9YRBvCRJkiRJPWEQL0mSJElSTxjES5IkSZLUEwbxkiRJkiT1hEG8JEmSJEk9YRAvSZIkSVJPbLQYG0myAbAbsC2w8Uxpq+qbi7FPSZIkSZLWN1MF8Ul2Al4PPAHYbA6r1LT7lCRJkiRpfbXggDrJzsAJwM5A5rraQvcnSZIkSdL6bpo28YcCtwAuBV4M3BrYuKo2mOm1CHmWJEmSJGm9NE3V9ofTqsc/p6o+sUj5kSRJkiRJE0xTMr4DcA3wmUXKiyRJkiRJmsE0QfxK4PKqumaxMiNJkiRJkiabJog/GrhJkt0WKzOSJEmSJGmyaYL41wGXAW9cpLxIkiRJkqQZLDiIr6pTgEcDD07ytSR7J9li8bImSZIkSZKGzal3+iTXzpLkId2LZMah4KuqpukRX5IkSZKk9dZcA+oZI3NJkiRJkrTmzTWI33uN5kKSJEmSJM1qTkF8VX1jTWdEkiRJkiTNbJre6SVJkiRJ0lq04CA+yalJvjOP9P+X5DcL3Z8kSZIkSeu7aXqK3xXYdB7pdwFuNcX+JEmSJElar63N6vQbAdetxf1JkiRJknSjslaC+CSbATsCl6yN/UmSJEmSdGM05+r0SW5Fq0I/bFmSBzJ5HPkAWwNPAzYGfryAPEqSJEmSJObXJv5ZwD+PzNsGWD6HdQMU8D/z2J8kSZIkSRoy347thkvci8kl8MNpLgZ+Aryrqj4yz/1JkiRJkqTOnIP4qno18OrB+yTXAWdX1c5rImOSJEmSJOmGphli7gPAhYuVEUmSJEmSNLMFB/FVdfAi5kOSJEmSJM1impL4G0hyE2AP2lByACuBH1bVxYu1D0mSJEmS1mdTB/FJ7ga8Btif1cedvy7Jl4F/rqofTrsvSZIkSZLWZ6NB97wkeRZwAnAAsCGtt/rh14bAI4ATkjx7uqxKkiRJkrR+W3AQn+TewLuBjYFTgOcCfwpsBmza/f1c4Be0Ev/DunUkSZIkSdICTFMS/4pu/eXAXavqPVV1alVdWVVXdX+/B7g78PUu7SumyWySbZN8Lcmvu+k2E9K9L8nKJD8ZmX9okjOSnNS9DpgmP5IkSZIkrU3TBPF7AgW8oKqumJSoqq4EXji0zjQOAY6pqt2AY7r34xxOa6M/zn9W1d2611FT5keSJEmSpImS7J/kl0lOSTIphp2zaTq22wK4uKp+OVvCqvpFkouAzafYH8CBwF7d30fQagG8fMz+vplk1yn3db1N60p2u/zUxdqcJEmSJGk9kGRD4O3AQ4HTge8l+VxV/Wyh25ymJP73wKZJZt1Gl/FNgdOm2B/AzarqLIBuuuMs6cd5UZIfdVXux1bHB0jy3CQnJjnxnPMuXGh+JUmSJEnrr3sDp3TNza8CjqQVTi/YNEH8p4FlwGPmkPYxwCbAJ2dLmOToJD8Z85rqQDvvBG4L3A04C3jTpIRVdVhV3bOq7rnDdlsvwq4lSZIkSeuZW3DDwuzTu3kLNk11+tcCj6X1On9+VS0flyjJg4D/ofVS//rZNlpV+05aluQPSXaqqrOS7ASsnE+Gq+oPQ9t6N/CF+awvSZIkSVo8v97sT5Y6Cwt2xrIzAO6c5MSh2YdV1WFD7zNm1Zpmv9ME8Y8D3gW8CjgmyfHAscAZXaZ2AfamdWZ3ES2Qf1yy+jFU1QfmuM/PAQcBb+imn51PhgcPALq3jwV+MlN6SZIkSZJm8OOqeukMy08Hbjn0fhfgzGl2OE0QfzirniAEeED3GhVgK+DNE7ZTwFyD+DcAH0/yHFqb/CcCJNkZeE9VHdC9/yitA7ztk5wOvKqq3gv8W5K7dftcATxvjvuVJEmSJGm+vgfsluQ2tALvJwNPnWaD0wTxv2fKagDzVVXnAfuMmX8mcMDQ+6dMWP8Zay53kiRJkiStUlXXJHkR8BVgQ+B9VfXTaba54CC+qnadZseSJEmSJN3YVdVRwFGLtb1peqeXJEmSJElrkUG8JEmSJEk9MU2b+Osl2Z7WE/2tgc2r6l8WY7uSJEmSJGmVqYL4JBsBbwT+Clg2tOhfhtJsA/wG2By4zdAQb5IkSZIkaR6mrU7/v8BLaAH8T4FrRhNU1QXAR7o0B065P0mSJEmS1lsLDuKTPIkWlK8E7llVdwHOn5D8f7vpIxe6P0mSJEmS1nfTlMQ/izZO/D9U1Q9nSfvdLu2dp9ifJEmSJEnrtWmC+D266SdnS1hVlwMXATtMsT9JkiRJktZr0wTxWwEXdQH6mt6XJEmSJEnrvWkC6wuArZJsOlvCJLsAN6W1n5ckSZIkSQswTRB/cjd98BzSPq+bnjDF/iRJkiRJWq9NE8R/FAjwmiSbT0qU5M+Bl9M6tvvgFPuTJEmSJGm9ttEU634AeD5wL+DbSd4JbAyQ5J7AXYAnA/vQgv1jquoL02VXkiRJkqT114KD+Kq6LsmjgS8A9wTePrR4uNp8uvdPWui+JEmSJEnSlD3GV9VK4AHAXwM/olWZz9Dr58BLgAdX1fnTZVWSJEmSpPXbNNXpAaiqq2ml8G9PsiVwc2BD4A9VdeG025ckSZIkSc3UQfywqroUOGUxtylJkiRJkpoFB/FJtgYeQxti7rbAtt2i84DfAMuBz1TVxVPmUZIkSZIkscAgPsnLgUOAmw7P7qYF7AkcBPxXktdV1X9MlUtJkiRJkjT/ID7JB4GnsipovxY4FTi/m7cN8Ce0dvFbA29Mcseqetai5FiSJEmSpPXUvHqnT/J84Gm0YP2HwBOBratq96q6X1Xdt6p2pwXvf96lCfDMJH+xuFmXJEmSJGn9MucgPsnGwGto1eU/Cty3qj5ZVZeNpq2qy6rqE8B9gSNpgfxrkyxqR3qSJEmSJK1P5lMS/2hgO+C3wHO6oeVm1KV5drfO9sCjFpJJSZIkSZI0vyB+b1op/Nuq6oq5rtSlfTutNH6f+WVPkiRJkiQNzCeIv3s3/doC9vOVkW1IkiRJkqR5mk8QfytaSfzPFrCfnwHXdduQJEmSJEkLMJ8g/qbAJVVV891Jt87F3HBceUmSJEmSNA/zCeK3BC6fYl9XAltMsb4kSZIkSeu1+QTxWYT9LcY2JEmSJElaL80niJckSZIkSUtoo3mmv1mSaxe4r9A6xpMkSZIkSQsw3yDe6vCSJEmSJC2R+QTxr15juZAkSZIkSbOacxBfVQbxkiRJkiQtITu2kyRJkiSpJwziJUmSJEnqCYN4SZIkSZJ6wiBekiRJkqSeMIiXJEmSJKknDOIlSZIkSeoJg3hJkiRJknrCIF6SJEmSpJ4wiJckSZIkqScM4iVJkiRJ6gmDeEmSJEmSesIgXpIkSZKknjCIlyRJkiSpJwziJUmSJEnqCYN4SZIkSZJ6wiBekiRJkqSeMIiXJEmSJKknDOIlSZIkSeoJg3hJkiRJknrCIF76/+3de7BlZ1kn4N/b6WAgIQkQQq4SR6PgjehER8TBKKhMpABFRnBgYGQmOKNTKIwOSJXgWFp4GbSmcICIlhkH8RZCooiShHC/BggkoRMScu2kc+2EdDqdyznnnT/OanPm2OfSnXN679Xneap2rb32+r5vvXtnVfr81hUAAGAkhHgAAAAYic2TLmAMHrrj9txy1tsmXQYAAAAbnCPxAAAAMBJCPAAAAIyEEA8AAAAjIcQDAADASAjxAAAAMBJCPAAAAIyEEA8AAAAjIcQDAADASAjxAAAAMBJCPAAAAIyEEA8AAAAjIcQDAADASAjxAAAAMBJCPAAAAIyEEA8AAAAjIcQDAADASIwqxFfV46vq/Kq6apg+bg9tTqyqi6pqS1VdXlWv3pv+AAAAMK1GFeKTvC7Jhd19cpILh/nFZpK8trufmuT7kvx8VX3rXvQHAACAqTS2EP/8JGcN789K8oLFDbp7W3d/fni/I8mWJMevtj8AAABMq7GF+Cd197ZkPqwnOXq5xlV1UpLvSvLpfekPAAAA66GqXjRcAj5XVaeutt/m9SxqX1TVBUmO2cOiN+zlOIclOTvJL3b3PftQxxlJzkiS4484bG+7AwAAwHIuS/KTSd6xN52mLsR397OXWlZVt1bVsd29raqOTXLbEu0OznyAf1d3v2fBolX1H+o4M8mZSfK0457Y+/JdAAAAYE+6e0uSVNVe9Rvb6fTnJXn58P7lSc5d3KDmf4E/TrKlu9+yt/0BAABglb6jqi5e8DpjvVc4dUfiV/DmJH9VVa9MckOSFyVJVR2X5J3dfXqSZyR5WZJLq+qSod+vdvffL9UfAACA/e+tH7tu0iXss5suuyVJLu3u1y7VZrnLxbt7nw4qjyrEd/edSZ61h89vTnL68P5jSfZ4PsJS/QEAAGCtLXe5+L4a2+n0AAAAsGEJ8QAAALCfVdVPVNXWJE9P8r6q+sfV9BvV6fQAAABwIOjuc5Kcs7f9HIkHAACAkRDiAQAAYCSEeAAAABgJIR4AAABGQogHAACAkRDiAQAAYCSEeAAAABgJIR4AAABGQogHAACAkRDiAQAAYCSEeAAAABgJIR4AAABGQogHAACAkRDiAQAAYCSEeAAAABgJIR4AAABGQogHAACAkRDiAQAAYCSEeAAAABgJIR4AAABGQogHAACAkRDiAQAAYCSEeAAAABgJIR4AAABGQogHAACAkRDiAQAAYCSEeAAAABgJIR4AAABGQogHAACAkRDiAQAAYCSEeAAAABgJIR4AAABGQogHAACAkRDiAQAAYCQ2T7qAMdhx+7358Fs/OukyAAAA2OAciQcAAICREOIBAABgJIR4AAAAGAkhHgAAAEZCiAcAAICREOIBAABgJIR4AAAAGAkhHgAAAEZCiAcAAICREOIBAABgJIR4AAAAGAkhHgAAAEZCiAcAAICREOIBAABgJIR4AAAAGAkhHgAAAEZCiAcAAICREOIBAABgJIR4AAAAGAkhHgAAAEZCiAcAAICREOIBAABgJIR4AAAAGAkhHgAAAEZCiAcAAICREOIBAABgJIR4AAAAGAkhHgAAAEZCiAcAAICREOIBAABgJIR4AAAAGAkhHgAAAEZCiAcAAICREOIBAABgJIR4AAAAGAkhHgAAAEZiVCG+qh5fVedX1VXD9HF7aHNiVV1UVVuq6vKqevWCZW+qqpuq6pLhdfr+/QYAAACw70YV4pO8LsmF3X1ykguH+cVmkry2u5+a5PuS/HxVfeuC5b/f3acMr79f/5IBAABgbYwtxGDQqaAAABNRSURBVD8/yVnD+7OSvGBxg+7e1t2fH97vSLIlyfH7rUIAAABYJ2ML8U/q7m3JfFhPcvRyjavqpCTfleTTCz7+har6UlX9yZ5Ox1/Q94yquriqLt4xN/PIKwcAAIBBVf1uVV0x5NNzqurI1fSbuhBfVRdU1WV7eD1/L8c5LMnZSX6xu+8ZPn5bkm9MckqSbUn+51L9u/vM7j61u0997KbN+/htAAAAYI/OT/Lt3f2dSb6S5PWr6TR16bS7n73Usqq6taqO7e5tVXVsktuWaHdw5gP8u7r7PQvGvnVBmz9K8ndrVzkAAACsTnd/YMHsp5L81Gr6Td2R+BWcl+Tlw/uXJzl3cYOqqiR/nGRLd79l0bJjF8z+RJLL1qlOAAAAWK2fTfL+1TScuiPxK3hzkr+qqlcmuSHJi5Kkqo5L8s7uPj3JM5K8LMmlVXXJ0O9XhzvR/05VnZKkk1yX5FX7uX4AAAAG577385MuYZ/df8tVSfIdVXXxgo/P7O4zd89U1QVJjtlD9zd097lDmzdk/ilr71rNekcV4rv7ziTP2sPnNyc5fXj/sSS1RP+XrWuBAAAAbCSXdvdrl1q43OXiSVJVL0/y3CTP6u5ezQpHFeIBAADgQFBVz0ny35P8YHfft9p+Y7smHgAAAA4Eb03y2CTnV9UlVfX21XRyJB4AAAD2s+7+pn3p50g8AAAAjIQQDwAAACMhxAMAAMBICPEAAAAwEkI8AAAAjIQQDwAAACMhxAPAlPuZ27886RIAgCkhxAMAAMBICPEAMMV2H4V3NB4ASIR4AAAAGA0hHgAAAEZCiAcAAICREOIBYEotvg7edfEAgBAPAAAAI7F50gWMwY6ZuXzkjvsmXQYAG8zP7OEz/x4BwMbmSDwAAACMhBAPAAAAIyHEA8AUentfN+kSAIApJMQDwIgI9wCwsQnxAAAAMBJCPAAAAIyEEA8AAAAjIcQDAADASAjxAAAAMBJCPABMGXegBwCWIsQDwMgI+QCwcQnxAAAAMBJCPAAAAIyEEA8AAAAjIcQDAADASAjxAAAAMBJCPABMEXeeBwCWI8QDwAgJ+wCwMQnxAAAAMBJCPAAAAIyEEA8AAAAjIcQDAADASAjxAAAAMBJCPAAAAIyEEA8AU2JvHxvnMXMAsPEI8QAAADASQjwAAACMhBAPAAAAIyHEAwAAwEgI8QAAADASQjwATAF3mgcAVkOIB4ARE/4BYGMR4gEAAGAkhHgAAAAYCSEeAAAARkKIBwAAgJEQ4gFgwtycDgBYLSEeAEbOTgAA2DiEeAAAABgJIR4AAABGQogHAACAkRDiAQAAYCSEeACYoLW6KZ2b2wHAxiDEAwAAwEgI8QAAADASQjwAAACMhBAPABOy1texuy4eAA58QjwAAACMhBAPAAAAIyHEA8AErNep706pB4ADmxAPAAAAIyHEA8B+5mg5ALCvhHgAOMDYSQAABy4hHgD2o/0VsAV5ADgwCfEAsJ/s72AtyAPAgUeIB4D9YFKBWpAHgAPL5kkXMAZ35MF735Hrr5x0HbAHRyW5Y9JFwBJsnwu8o2rSJfAw2ybTyrbJNPuWSRfAPCF+da7s7lMnXQQsVlUX2zaZVrZPppVtk2ll22SaVdXFk66BeU6nBwAAgJEQ4gEAAGAkhPjVOXPSBcASbJtMM9sn08q2ybSybTLNbJ9rrKp+o6q+VFWXVNUHquq41fQT4lehu22wTCXbJtPM9sm0sm0yrWybTDPb57r43e7+zu4+JcnfJfm11XQS4gEAAGA/6+57FswemqRX00+IX0FVPaeqrqyqq6vqdZOuB5Kkqv6kqm6rqssmXQssVFUnVtVFVbWlqi6vqldPuibYraoOqarPVNUXh+3z1yddEyxUVQdV1Req6u8mXQvsVlXXVdWlwynfa32H+i/t2PKBNR5y/+jZh3Lvl9+fJIdW1cULXmfszThV9ZtVdWOSf5dVHomv7lWF/Q2pqg5K8pUkP5Jka5LPJnlJd395ooWx4VXVM5Pcm+T/dPe3T7oe2K2qjk1ybHd/vqoem+RzSV7g/5tMg6qqJId2971VdXCSjyV5dXd/asKlQZKkql6T5NQkh3f3cyddDyTzIT7Jqd19xzqM/ag6/MQH6vATctBRT1nr4ddN91xmr70odfhxmb3xk7Vc26q6IMkxe1j0hu4+d0G71yc5pLvfuNL6PSd+ed+b5OruviZJquovkjw/iT9Gmaju/khVnTTpOmCx7t6WZNvwfkdVbUlyfPx/kynQ80cu7h1mDx5ejmYwFarqhCQ/nuQ3k7xmwuXAftHdD1bVIem5++dqUzY94ZsnXdKKuucye/2HU4cds2KAn2/fz17l0H+e5H1JVgzxTqdf3vFJblwwv3X4DIAVDDuavivJpydbCTxsOF35kiS3JTm/u22fTIs/SPIrSeYmXQgs0kk+UFWf29tTxVc1ePcDveOmx8zdfW3mtn91rYdfU91zmb3ho6lHPyGzN316xQC/kqo6ecHs85JcsZp+jsQvb0//YeyxB1hBVR2W5Owkv7jopi0wUd09m+SUqjoyyTlV9e3d7f4iTFRVPTfJbd39uao6bdL1wCLP6O6bq+roJOdX1RXd/ZG1XEF376qqw+Z67t5s2pRNR37DWg6/Jro7szd+IvV1h2f25osfcYAfvLmqviXzO++uT/Jzq+kkxC9va5ITF8yfkOTmCdUCMArDtcZnJ3lXd79n0vXAnnT33VX1oSTPSSLEM2nPSPK8qjo9ySFJDq+q/9vdL51wXZDuvnmY3lZV52T+kuM1DfHD+Dur6oi57q+lNmXTEU9e61Xss+7O3NZPpjYfkrlbLlmzs9m7+4X70s/p9Mv7bJKTq+obqupRSV6c5LwJ1wQwtYYbh/1xki3d/ZZJ1wMLVdUThyPwqapHJ3l2VnnqIqyn7n59d5/Q3Sdl/u/NDwrwTIOqOnS4UW2q6tAkP5p13PHZ3ff0zlsfN3frpZm7Z+t6rWavdHfmbvpMsmlz5m67dFNPwZ3hhfhldPdMkl9I8o9JtiT5q+6+fLJVQVJV707yySTfUlVbq+qVk64JBs9I8rIkPzw8iuaS4cgSTINjk1xUVV/K/I7687vbo7wAlvakJB+rqi8m+UyS93X3P6znCrv77r7v9qPmbvlC5nZM9iTo7s7ctouTnsvc7ZdPRYBPPGIOAACAKVNVR9ejj7p10/Hfk02H7ekJbetvdtsX0g/tTG+/6qDunpqbTjoSDwAAwFTp7tt61x3Hzd70mcztvG2/r3/21i+mH9yR3n7V5mkK8IkQDwAAwBTq7m3ZdefXz279ZObuu2O/rXf2tsvSu7an7/rqwcNTTaaKEA8AAMBU6u4bs2v7N8ze8PH0ru3rvr7Z27ekd96avvu6Rw33SJs6QjwAAABTq7uvy/3bT565/qPpXXet23rm7vxKesfW9NduOKS7H1q3FT1CQjwAAABTrbuvzv3bnzpz/YfT939tzcef23515u6+Nn3P1kd39wNrvoI1JMQDwJSrqj+tqq6qP510LcsZauyqOm3StQBw4OnuKzY/+ZmZue5D6QfuWbNx5+66JnPbr07vuPnQ7r5/zQZeJ0I8AKNTVW/aHRgnXcukVdUrht/jtEnXAgDr7aEr3lubv/4ZmbnuovSD9z7i8ebuvj5zd1yRvnfb4d193xqUuO6EeACYftuSXDlMF3tFkjcmOW0/1gMAE/PQlefVQSd8f2auvTD94M59Hmfuazdm7vbL0jtvPbK7d6xhietKiAeAKdfdr+/up3T36yddCwBMg5mv/G0ddPz3zQf5h/b+APrcPTdl7tYvpnfe9oTuXvuL7NeREA8AAMDozFz1vjrouO/J7LUXpmdWfyn73I5tmdv2+fR9tx/d3ev/3Lo1JsQDsKFU1WlV9ddVdVNVPVBVd1TVhVX1H6rqoCX67L4G/0PD/LOq6n1VdXtV3V9VW6rqjVV1yArrfv6wrrur6t6q+mJV/UpVHbx4HYv6/bMb2w3XwneSHxw+euOCG8vtfp00tD1p8WdL1Hfd0OYVSyx/XFX9blV9dfje24bf8l8u970XjfGCqnpvVd1cVQ9W1V1V9ZGq+rmqOni14wBAksxc/f7adMx3Z/aaC9IzK99Ufm7nrZm9+bPpXXcc292374cS15wQD8CGUVVvSXJRkp9KcmyS+5IcmeSHk/xJkg9U1WNXGOOXk5yf5N8k2ZzkUUmekuRNSf5+mR0Bv5fkvcO6jkjyUJJvTfLbSS5IsrcBdleSW4dxkmTnML/wNbuXYy5pCP9fSPLfkvyLJJXkMZn/LT9RVc9bof9hVfW3Sc5J8vzM//67Mv9b/Oskb0vy4ap63FrVDMDGMPPVf6xNT3paZq+9ID374JLt5u67PbNbP5XsuvOE7r5lP5a4poR4ADaEqvqFJL80zJ6Z5LjuflzmQ+QvJZnJfMD+o2WGeVqSNw+vo4f+Ryb5H8PyH0ry8j2s+8VJXjvM/nmSE4a+j01yRpLvTfKf9+b7dPdfdvcxST4xfPR73X3MoteNezPmUoYdE3+d5MlJ7kryb5Mc2t1HJPm2JJ9OctYKw/xZkucmuTrJzyQ5fOj/mMyH+muSPD3zO1MAYK/MXHN+bXrit80fkZ996J8t7/vuzOyNn0h2bT+pu2+aQIlrRogH4IBXVY9O8uvD7Lu7+1W798B3987u/oMkrxmW/3RVnbrEUEcm+Y3u/tXuvmPof093vzHJe4Y2L1m07srDIf/8JC/d/cdDd9/f3X+U+QA/zUegX5hk92/you7+6+6eSZLu/nKS5yS5c6nOVfXjSV6Q5JYkp3X3u3ffBXj4Dc7L/GUBO5O8oKpOWb+vAsCBaubaD9amo54yf4383Mw/fd677srMDR9Ldm3/pu6+foIlrgkhHoCN4EeSPH54/6Yl2vzvPPwIt5cs0eaBJL+3xLJzh+l3Lvr8lCQnD+9/q7v39Gz7s5LcsMS40+DFw/Tj3X3h4oXDc3V/Z5n+/3GY/tlSRz+6e2vmL3VIkh/b10IB2NhmrvtQbXr8N/5TkO/7787M9R9J7t/+lO7+6qTrWwtCPAAbwe6jyDd291f21KC7Z5N8cFH7xS7v7nuXWHbzMH38os+/e5g+lIdPfV+87k7y4SXGnQa7f48PLtNmuWU/MEzPqKpblnolefbQ7smPtGAANq6Z6z9am444KbPXXpiZ6z6czU9+Zrr7yknXtVY2T7oAANgPjh6mK10Dt3VR+8V2LNN393l7i/9tfeIwvbO7l77bzsq1TdJqfr+te/pwuOP8UcPsEcNrJY9ZfWkA8M/N3PjxqqoXJrnqoSvO+dKk61lLQjwAG8meTmV/JO1Wo1Y5Zq2wfBos9x2WWrbwbv0v7u6/XMN6AGBJ3X32pGtYD06nB2AjuG2YnrhCuxOG6Vo+N3b3uo+qqkct0+64NVznYjML3i/3LPuljpLv/g4nLLF8yWXdfX+Srw2z37FMfwBgFYR4ADaCi4fpCVX1zXtqMDxG7YeG2c+u4bo/P0wPTvL9S6y7kjxzH8ef2z3MMm3uWvB+jzsyht/lyCX67/79fmiJ5cn84/mW8vFh+qKq8rcHADwC/iEFYCM4Pw8/Au1NS7R5VR4+Gv7uNVz3JZl/NnqSvG4I7Iu9NPt+M7d7hulSATzdvTPJ7jvyvnCJZm9YZh27T4H/gao6bfHC4RF+v7xM/zOH6Tev0C5VdegKZywAwIYmxAMwalV11AqvI7t7Vx4O7y+pqrdX1ZOG/o+pqv+a5A+G5X/Z3Z9bq/qGO8+/cZj9sSRnVdVxw7oPqapXJnlH/v+j5XvjsmF6elUdv0y73Tsmfraq/ssQvFNVJ1bVO5P8dJL7luh7dh4+o+DsqnrhcOZCquqpSd6fpW8GmO4+N8k5w+ybq+ptC8+IqKpHVdW/qqrfTnL9cmMBwEYnxAMwdrev8PpQknT3W5P8/tDnVUm2VdX2zF+v/b8yf7r7RUn+01oX2N1/nod3ErwsydZh3fckeWeSTyZ5+7D8/r0c/qyhzzcluWF4XNt1w2vhdeq/neTLmf+ef5jk3qq6K/PPp//3SV6RJe4F0N0zSV6U5MbMP0Lvb5LsrKq7hzGfPoyxnJcm+Yvh/c8lubKq7h1+h11JPpXkV5I8IWt7Y0EAOKAI8QBsGN39msxfu312kluTHJb5x8ZdlORnk/xIdy/3GLlHsu5fSvKTmd+psCPJ1yXZkvnTy38syaFD07v3ctyrMn+t+nmZD+FPyPyp+U/OgqfQDM+3/4Ekb0lybeZvdvdQ5n+Lp3f3X2QZ3X1NklMW9K/M7zz4myTf393nrdD/vu5+yVDrnyW5JvN/hxyW+RvnfTDzIf7k7p7mx+0BwETV/Fl+AMAkVdXHM3/ju1/r7t+YdD0AwHRyJB4AJqyqfjAP37n+HyZZCwAw3YR4ANgPquoPq+oVVXXM7jvUV9WRVfWqJOcOzT7Y3Wv5eDsA4ADjdHoA2A+q6pIkTxtmH8j8neCPzMPPd/9ykh91PTgAsBwhHgD2g6p6XpKfSPK9SZ6U5IjM353+8iTvSXJmdy/1iDcAgCRCPAAAAIyGa+IBAABgJIR4AAAAGAkhHgAAAEZCiAcAAICREOIBAABgJIR4AAAAGIn/By6mMB4HwUMDAAAAAElFTkSuQmCC\n",
  89.       "text/plain": [
  90.        "<Figure size 1080x648 with 2 Axes>"
  91.       ]
  92.      },
  93.      "metadata": {
  94.       "needs_background": "light"
  95.      },
  96.      "output_type": "display_data"
  97.     }
  98.    ],
  99.    "source": [
  100.     "visu.show_initial_conditions(croco_config).show()"
  101.    ]
  102.   },
  103.   {
  104.    "cell_type": "code",
  105.    "execution_count": 49,
  106.    "metadata": {},
  107.    "outputs": [
  108.     {
  109.      "name": "stdout",
  110.      "output_type": "stream",
  111.      "text": [
  112.       "movie was saved in lab_waves_movie.mp4\n"
  113.      ]
  114.     },
  115.     {
  116.      "data": {
  117.       "image/png": "iVBORw0KGgoAAAANSUhEUgAAA3wAAAIMCAYAAACuQa5XAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJzt3XvMbtddJ/bfahxUD4SSNI7tXKoTo1Oom8kkwaJV4ySFxFXwjMZG8w9BjY4qJBeLIwEdqeOq/YP5o1JEC5TqjM7IlGjOaAjRSKF2xERDHZeJYwkYbDcYB096oiNXOL4lRilBHIlmWP3jPK/z+vFz2fe99tqfj3R03vd592Xt/ezbd6+110455wAAAKA+/87cBQAAAGAcAh8AAEClBD4AAIBKCXwAAACVEvgAAAAqJfABAABUSuADAAColMAHAABQKYEPAACgUgIfAABApa6buwBdvPnNb85nzpyZuxgArNif/dGX401/6z+euxgArNTjjz/+jZzzDceGW2TgO3PmTDz22GNzFwOAFfvUDbfGTzoXATCTlNL/02Q4TToBAAAqJfABQEefuuHWuYsAAAcJfAAAAJUS+AAAACol8AEAAFRK4AMAAKiUwAcAAFApgQ8AAKBSAh8AAEClBD4AAIBKCXwAAACVEvgAAAAqJfABAABUSuADgIF86oZb5y4CALyKwAcAAxD2ACiRwAcAAFApgQ8AAKBSAh8AAEClBD4AAIBKCXwA0MDpTlkOddCi8xYASiLwAcARQhwASyXwAQAAVErgAwAAqJTABwAAUCmBDwAAoFICHwAcoMMWAJZM4AMAAKiUwAcAA1MrCEApBD4AAIBKDRL4UkofTSl9JaX01ZTSfTv+nlJK/+vm70+mlN7XdFwAmIuaOgCWrnfgSym9LiL+UUT8WETcGhEfSyltnyF/LCLObv7dExEXW4wLAEUQAAFYmiFq+H44Ir6ac76Sc/6riPh0RNy1NcxdEfFP8zW/HxHfl1K6ueG4ADC50+HukW/8ZTzyjb+csTQA0M0Qge9tEfGnp35/dvNZk2GajAsAi6M2EIASXDfANNKOz3LDYZqMe20CKd0T15qDxtv+ve9pUz4AVuT8Z57c/fntZ+Ls1Suv+uyFSxdf+fkLF774qr+p0QOgBkMEvmcj4h2nfn97RDzXcJjvajBuRETknO+PiPsjIv7WW2/IX/uH/3XcdO7euHz9LXHh0WdeM/yFv/fudksBQCv7gtUSbYe9oXzqhlvjJ7/+J42H/7fPfOlVv1++/paIiPjBG7930HIBsB5DBL4/jIizKaV3RsTXIuInIuInt4b5bEScTyl9OiL+k4j4f3POz6eUvt5g3E5quhABYHina/e2DVm7t92080PnP/CaYW46d+/BafybF/98581N5uPGMjCGYxnmV3/or+OFSxePnjdO6x34cs7fTimdj4jfiYjXRcQnc85fTin99Obv/zgiPhcRd0bEVyPiLyPivzo0bt8yAUBXmnLShBvLwNTO334mYuvRhCaGqOGLnPPn4lqoO/3ZPz71c46In2k6blMvXLoY8dP/U5dRAaCVD775b8xdBACIiMOtVLYN8uL1qX3r63/xys/bD+ADQMl2PS/Y5sQNwHp1OV8sMvC94Qa9dAIAAByzyMAHAEM5XePm+T0AaiPwAcDExnoNBABsE/gAAAAKdv72M537LhH4AKAAOm4BYAwCHwAAQKUEPgBWp21tmnfwAbBUAh8AAEClBD4AAIBKCXwA0MCHzn9g0Ol5NQMAU1h84NOrGQBz+9QNt85dBADYafGBL+LaeykAYA7CHgAlqyLwAUAXmlUCUDuBDwAO8EoGAJZM4AOAgQzdsQsAnOjad4nABwAN3HTu3rmLAACtCXwAAAAFO3v1SudxBT4AmIHmnwBMQeADYBUuPPrM3EUAgMkJfABwhNo4AJZK4AOAPUp5JYPaSQC6EvgAYEBqAwEoicAHAABQKYEPAACgUgIfAKtx+fpbRpv26aacmnUCUAqBDwAmJhACMBWBDwB2KKWHTgDoY5GB7/VvvmHuIgAAABRvkYEPAMYwZK3evmabmnMCMCWBDwBGsh3uhD0AhnDTuXsbD3vdiOUAgKJ96PwH4gsXvviqz4Z+dk/IA2BOi6/ha5NuAaCt7cD2qRtunakkANDeYgOfoAcAAHDYYgMfANRk+0bmmC+JB2C52lZ8CXwArM4QrUQ8mwfAEgh8AAAAlaoi8J29eiXO335m7mIAAAAUpYrABwC1uvDoM3MXAYAFE/gAWA3hCYC1EfgAAAAqJfABwMy8WxaAsQh8AAAAherbOaXABwANefceAEsj8AGwKpevv2XuIgDAZAQ+AGhoqmfthFIAhiLwAUBLh5p2tg2FOmwBYEwCHwA01LTmTYgDoBQCHwCrpiMWAJaiyw1FgQ8AWhqqBk9NIABjE/gAYI8+tX9DhLkLjz7TexoArNuiA587owB0NcU55NA8nMMAmMKiAx8AzMFrEwBYCoEPAFpo28zzpnP3vqY271DtnjAJwJCum7sAADClC48+E7/6Q8NN73QA1IQTgNKo4QOAFgQ3AJZE4AOAI4Q8AJZK4ANg9bx8HYBaCXwArE6fjlF0qgLAkgh8ANDSWC9EPx0mvXQdgNO6Pl4g8AEwivO3n5m7CKPyXB8ASyDwAcAOnusDoAYCHwCjWVItn4AHQI0EPgBWS7NMAGon8AFAAfT+CcAYBD4AVudYD5iadwJQi+vmLgAAy3b+9jPVvEJA0AOgNr1q+FJKb0opPZRSurz5/417hvtoSukrKaWvppTuO/X5L6SUvpZS+tLm3519ygNAGZbUWUsJNOcEYCx9m3TeFxEP55zPRsTDm99fJaX0uoj4RxHxYxFxa0R8LKV066lBfiXn/J7Nv8/1LA8ALF4tNaYA9DPEDdS+ge+uiLi0+flSRNy9Y5gfjoiv5pyv5Jz/KiI+vRkPgIJ0OaksuSavS63a2atXRigJAIynb+C7Mef8fETE5v+37BjmbRHxp6d+f3bz2YnzKaUnU0qf3NckFIBxDRnclhwCAaA2RwNfSunzKaWndvxrWkuXdnyWN/9fjIjvj4j3RMTzEfFLB8pxT0rpsZTSY19/+ZsNZw2wXOdvPyM8TWDud/F5fg+AMR3tpTPn/JF9f0spvZhSujnn/HxK6eaIeGnHYM9GxDtO/f72iHhuM+0XT03r1yLitw+U4/6IuD8i4rZ335r3DQfAspyEWs+tAcDw+jbp/GxEnNv8fC4iHtwxzB9GxNmU0jtTSt8VET+xGS82IfHEj0fEUz3LA1CdpdfyLb38Y1K7B8DY+ga+T0TEHSmlyxFxx+b3SCm9NaX0uYiInPO3I+J8RPxORDwdEf885/zlzfi/mFL645TSkxHxIxHx8z3LA8DG2jphqYWaTgCG1Cvw5Zxfzjl/OOd8dvP/n20+fy7nfOep4T6Xc/4Pc87fn3P+H099/vGc89/MOb875/x3TzqAAdZFyGiu6bqyTo+bO1ip3QNgCn1r+AAGIaCwJsIeAFM52mkLwNx06lGvtQV9QQ+AqanhAyjEUOFnqhDVdT4lhbyTADbWqxkuX3/Lq/4BwNQEPoCBlRRohnCyPLUt19gEPABKIPABrFyTICfsTUOzZQCGJvABizFW6BBm9rNu9jf3fOHSxYlLAgDtCXysngtaxtB2u9oe/tD4Q26ztv92NNMEYGkEPoBYZ/BZ4zJvO92EcqyOWxiG7RWgG4EPVqrEi6cSyzS1LjWDU9UGNtFnfjqHYR/bBEB3Ah+r5gKTpta0jUwRIEtan6U009RhCwBjqCrwlXQBccixO/IwNtvfbtZLdzUGQRhL0+3c/gAMoarAx7I5sVGCGl4mXrLSL3S3n+PzXB9Dc6wApibwUQQnQJpa4rayxDLT3BcufHHuIlTt9P5jXwJoT+BbuSlPnk7U5SnhO5m7DG1eOr5v2LmeeWv6WdPpcc0QtXqlPBdIuY7tg0Pso/ZzIELgowBLPiEtuez79AkRsES7Oku56dy9kzbn1GEL+zj+An0JfCs2xknEiWnZfH/XtFkP1lkd1Mitx5z7rOMFMAeBb2BDNtFY2onhULO3pq8/WNoyRyyzzMfUuExzqHF7h5IsvZOlUsrRhB7GYbkEPlZ7AO+73GtZb95VuN9Y62QJ69rzRcNZU3POIb/ztW8/Qz/Xe2geNa/rps9x17wO+rJuyifwjWCqC6G17GClLKfOMF5rzGUduyOUqb+n7fl1mf+ati2m07Um2vY4npLON0sP6cLaNJa6jpda7rYEPhavy8G89B18CSeoqcpXajCa+yJ4aRd6fZRSjjU7tr0v4Zi1bUnl7dKbsFYszQ11g3Go/WBN635Oa1rPAt+MhmhmUapD703q8hxjk4vrKe6IDvVM1tB3TOd8NnTubXHM4DX3srVxaDtY0nJsm7LsOm55ralamwx1XJ67Q5ZSbviUss/vWidTlq2U9cA1S/o+llTWJgQ+BjHkg/Ml7mR9axDbXJCU1AnBvuUo5aKm7XAlWmLZl1jmY+Z6jm5Jz++VUovXpLax6/SGOP7OcTO3xBA1drP8kg2xDczZ/HXM64lSvvuhKibGNlRZBL6RDPUFzX23smRtTjqlHWgiXnuwaRp+l1Qz3PUgOnZ4HatGtKTtaw5LqXVgPE0u8ofYp4YKRiVtkyWVpa0Syj5FGbrMY+wbJTVcb9aslBtlAl8P21/iEna6thfch062JWzAbcxd3iGaRpUYXCN2l6vJhdqQgbbJuinlYm/O+Za27Ryya1tZ2jKUroR9opTvc6xyjL18bc4tY95069sSpo8haneHKsPU8x1ieku4fqUfga8yJe1sc9zt7xPC+9xJHvMidIoT7xDmrj0catha1bAOdtWKD71cUz7HV1JzzrkufMcqw5C6bGel3pyLmC9sDDHfqc+1YwapEraNsdZnCTeS2mp607iUcN1WdYFvqovAKb+woQ+mxw5qpTaLKGX6JVwYjTGPpgfoki7mupxUhtiflnDyOqS2GwJT1frp1OU7Sjp27Rp2afvoFOfgpZ27us5/7rIdU0r5xgrqbc/DUx2/5whqx45Lx5Z9yPVSXeA7rclGNNQdp67T77sBzVHL0Wb+pRzY5jTmOhrjjuQQSigD69P3wqGkGrc5zbn/1tTUrem5srQbMGNfFw01r2M1Mk0utptMv8nnQ52L+9x4Le3arG0t95i14tvTnrJWuhRVB75DpjyptNnxxqg9KTUUHLP0O7YRwy1DlwvZ7XGmWn9N7mBNdUdvDH0vUJayP5ZYpqaWUPbaw2WTY99cx6Suw4wxbinzmXufOX1sbHOcnLvcbe1bhq7hsMm135Chet/5fChDlrnL39sE7D7lmOMaaLWBb59DB5q+dwaGuOPSZpi2F/tNNvQSDq5T3pkZcucf2hy1u0uz5GBJP32+96ZNN79w4Yud5zGHMc4Jbee7/dmYF/HbyzvGdKcYb+hpbE+rz0X2EOfjfdvlsQDYx9zbQp95lHCDYahteeh9adc207esTY5fJ/PrOs8pvtPrRp/DApys6NN3XE9/NkTg6Rv22u54Q290Y+/4bcsx1s4xxsl0SmOesGFpzt9+ZrCatLU9v1daOJlS0/PgkNtXG2Oep0r73se45hrieqbLdeEYwWPXMLu2ydPjDrnNDlGRcWzc7fK2uQ4c8xjUtxJo6OGPSTnnQSc4hdvefWv+g89+6lWfre1kDMBxXS5uTk60Z69eeeWz7XPMhUefif/sp3/ild8/dP4DERFx07l7dw7ft0xDOH0xePLz9oXgUkMawNIdOwafPieduO6d730853zbsWlr0glAtYYIMDXcUJz7bjgAh415DBb4AKha15Po5etvqSLs7SPgAaxDNc/wnb16peoTMwDLN3fvnEIewPqo4QOgem2CztyhDACGJPABQAtLCoRq9AAQ+AAAACol8AGwCpp1ArBEu17J0IbABwA9nLyDrzRtXwwMQJ0EPgBWY87wo9YQgDkIfACww1IDmho9AE4T+ABgj0Ohr8SmnMIeANsEPgBWZahQ9IULXxxkOmMQ/AA4IfABQAWEPAB2EfgAYMuDDzzxys9LfZYPACIEPgAAgGoJfACszrHmj3fd/b5X1fLtU2LHLQBwmsAHAABQKYEPAPY4qeUr/Tk+HbYAsI/ABwAAUCmBD4BVavIcHwAsncAHAAc06bxlTppzAnCIwAcADV149JmiAqCwB8AxAh8AAEClBD4AVqtNDVnpPXUCwC4CHwDsMVTHLWOERc05AWhC4AMAAKiUwAcAR5TUUQsAtCHwAbBqXZpGak4JwFIIfAAwgJvO3TvZvAROAJoS+ADggF0dt5y9emWGkgBAewIfAIxsyBo5tXsAtCHwAbB6bUPUC5cujlMQABiYwAcAAFApgQ8AFkJzTgDaum7uAgDAEpx+F98XLnxxxpIAQHNq+AAAACrVK/CllN6UUnoopXR58/8b9wz3yZTSSymlp7qMDwBjK725ZOnlA6BMfWv47ouIh3POZyPi4c3vu/yTiPhoj/EBAABoqW/guysiLm1+vhQRd+8aKOf8SET8WdfxAWBOu16+DgBL0Dfw3Zhzfj4iYvP/WyYeHwCqpzknAF0d7aUzpfT5iLhpx5/+++GLc7Ac90TEPRER/8Fbb55y1gAAAIt0tIYv5/yRnPO7dvx7MCJeTCndHBGx+f+llvNvPH7O+f6c820559tu+Pe/r+VsAOA4NWkA1KZvk87PRsS5zc/nIuLBiccHoFBnr16ZuwhVEEIB6KNv4PtERNyRUrocEXdsfo+U0ltTSp87GSil9JsR8XsR8QMppWdTSj91aHwAAAD6O/oM3yE555cj4sM7Pn8uIu489fvH2owPAABAf31r+ABgL806+9GcE4C+BD4AOEXIAqAmAh8AAEClBD4AKJCaRgCGIPABQAt33f2+uYsAAI0JfAAAAJUS+ABgAm2aaGrOCcBQBD4A2CJwAVALgQ8ACiJsAjAkgQ+AXna9XN0L1wGgDAIfAABApQQ+ACiE5pwADE3gA4AGvH8PgCUS+ACIiG7P3Q3xrN7Zq1dW88zfoRo8tXsAbBvi/CjwAbCawAUAayPwARRsiUFsiWXu4kPnPzDYtNTuATAWgQ+gUGsJTqUSwgCogcAHABPaDpKCJQBjEvgAYCbCHgBjE/gACrf0pp1LL/+2IULa+dvPCHsATELgAwAAqJTAB1Cpud6rBwCUQ+ADYDYCJgCMS+ADZueiv7mm68o6BQAiBD6gEAIKAMDwBD6geMJgXZb+fS69/ACsi8AHLIKLbACA9gQ+AACASgl8AIUYqhZzqtpQta4AUD6BDwBaeOHSxdd8dtO5e2coCQAcJ/ABDKy2mq/algcA1kTgA+Cos1evCH4AsEACH7AYAsc4SlmvpZQDAGoi8AGrJ2jsN8e68X0AwHAEPlbPxSVjaLtdbQ9/aPwht1nbPwDUTeADoBhLCKBfuPDFuYuwSkvYNgBKJPABhItJAKBOAh+sVIkBp8QyTW3KdTDn+j42b9sCp9keALoT+Fi1k4sIFxMsldclAACHCHwzcYHGnGx/7S11nXUpd59lXep62uf87WfmLgIA9CLwAWzUFlbGsG8dNV131jHYD4BpVRX4lnIAXUo5AXi1u+5+39xFYOHaXAO4XgCGUFXgY9mc2ChB1+3Q9gsAlEjgowgulqmZ7Ru6O73/2JcA2hP4gEWZ44LPReaynfRk2vd7fPCBJwYqETiuANMR+FZuLe/8YrcSvpMSytBUlw5Lxly+Ja27OewKeV5jwVIMsZ3a1oEIgY8COCGVZdf34Tua39C9YNb+nY71YnevaQBgaQQ+6KH2i+Y18B3WpU0Nnu+eqU39XkyACIFv1cY4iTgxLZvv7xrrYZl8bzRhOwHWRuCb2BpONIeaBI7VzArm0uf5Pdt7M03Wk3W5Tl6jMh3Pv8JyCXwDG/JguLQD69LKO5Qal1s36ExhqG2rz3TG2r7tN+OrYR1PsQw1rCegH4GP1Z4M1rrcXVlfryUsjO9YDerU68p304/1N5wpOthq2jpnyZq2IKh5HfS15HWz5LK3IfCNYC0bz9rs6t69yXA1m+uVA7orn99c62+od+rxHV2bHvsOxlNSx0NrbrlEc77bsgl8hfLMynesZTm3lb7cpZevNE1vGLSZBsOwXo9rEvqsx3mtoTauiyY3iIa6wehm1LKs6bsS+Fi8Gi+cSy/flEr9fueu9VjjNrLGZZ5bSTVNQ1tSmec4Di5p/dTGumdoAt+M1rpDj3HiWtq6HOs5iznmXbMlrasxmsCWtPyllOWFSxfnLkJx5m5iXcq2oXZntyFaNwD9CHwLVfoBc6hmPm2C3hR3oocInnOGvTGUOv+5yzWH2pa5tuVZuiV8HyXUapXwXN1cz1dPWY4lzL+POZ+THmu6JX0fx64fSynrUOUQ+JhVyTvZaVOWqaTl33dntqQyRpRXnjaWWPYaagzppoSauGM3edrOp6SQV9O+NXbHWyUrYT/pM+9SbhjMZQllbEvgG8ncB/+p5136zjF1+Zo0YenazGVNNVpDLlPXWuea9oOxLXn5l1z2kgzdymHo453mheMoYT2WULM61DglTZ9+SqnZFPh6KuFLXKsuAanUAD1ViCst+A653PZFaK+EAFTKvjt3Oeae/9TGWN5S12Gp5ToxVCgpfTnXTODrYawTZQk7zLFq/VLuWLSxhPI2fS6jtGXZVa6uYW7MZ3BKuLide76lbTuUoYTtooQyjKGEGp7Tx+ixamGX8lzhEOM1ndZU11BLbUJaktrXQ3WBr/YvrK+1PYtWw/IupbarhO97n5LLBrWpbX/rcnFe6s25pRvzplWbwFaLqdblktbjsZvGY96U7jPsMdUFvqmUuvGWWq5dxmrOMfSdvqHvhpZqactQwl3zJZm7aXNt65OyLXl7W0LZ5y7j3PPvo5Syd3ksZqjhl3gOGuPm+5StbqoOfHM2LxjSWDvlnJZ8sb6k9RxRRlOXMabXdR6eU4DDSrlgm+tO+NSmuugrcR2N0RSxlGsm1yHNzPGowdoe54ioPPBFlHdSGnp+Q7WzbzPNrsOu1ZjNA+ZoglLKxeCU82G9Ljz6zKt+/8KFL+78mWGUdqNpqONzaceqY+elMac/tKbBb6iaqinPu11vepRwbda2WfOYzaC7Tru0/baP6gPfPkNvWF3vNNW0MTGeOe6Ajam25YFdStvGpwwgTS62a7kxVNr3vESnzwk1nx+mvCYcI7iP/d3MXUlTwo2Ksaw28E1pyBqzLvNY6oGz7925Esy5DNsH5rkurmo7efe94z/nhW8bJZZpbGtc5kNqWx9znyfnrG0ca9pDNJ/sM42xyz2FY+fIfdcRTc+tQ6yjQ9cyY9bIdR13jpYC2+uky3fadp5t9Ap8KaU3pZQeSild3vz/xj3DfTKl9FJK6amtz38hpfS1lNKXNv/u7FOeIRz6krpWT+/7/djwXYdpMmybjaukA2MbYzZLKWGdDLUt1Ky2sAklGKJZ1KHz6Rgtb0qrOSlh3tvfQQnHyl3fWcnHcefh7vp8r02uzYbarpscv0rah/bpW8N3X0Q8nHM+GxEPb37f5Z9ExEf3/O1Xcs7v2fz7XM/yRET7jajJ8HM3Cel7UBkq/JSwMZdSjmPmKGOXeS5hXQLTGvOGWekO1cQvdZn2GaMVypLDddML99PbwtQtd7r+fY7tt29N3aHPT5ZniptHx8rRZvihh23qup7j3xUR//nm50sR8a8i4h9sD5RzfiSldKbnvDrbXnGXr7/l6HBD1+xMXdvU9UDUZbmnrM6fO3jXaq3LDYxrznNmX3PfCC5tvkvS5xqipADXZdgxarl2TWOK9XT26pW4fP0txe6Lu5R6zOtbw3djzvn5iIjN/2/pMI3zKaUnN80+dzYJHVpJGwYAjM15b1jW5/Q1KrWb6sb60r6jpZW3VEcDX0rp8ymlp3b8u2uA+V+MiO+PiPdExPMR8UsHynFPSumxlNJjX3/5mwPMGgDW5VBLFhdWAHU62qQz5/yRfX9LKb2YUro55/x8SunmiHipzcxzzi+emtavRcRvHxj2/oi4PyLitnffmtvMBwD6uuvu90X8y7lL0c5Jk6jtzwBYj75NOj8bEec2P5+LiAfbjLwJiSd+PCKe2jcsAAAA7fQNfJ+IiDtSSpcj4o7N75FSemtK6ZUeN1NKvxkRvxcRP5BSejal9FObP/1iSumPU0pPRsSPRMTP9ywPAAAAG7166cw5vxwRH97x+XMRceep3z+2Z/yP95k/AJRsV5NKAJhS3xo+AKAgS3oZMADjE/gAoBLCHQDbBD4AAIBKCXwA0MGHzn9g7iK8yqF37AGwXgIfAABApQQ+ANg4VDN2/vYz0xUEAAYi8AHAEQ8+8ETncTWvBGBOvd7DBwDMS6AE4BA1fACwUMIeAMcIfAAAAJUS+ABggdTuAdCEwAcAIxs6nAl7ADQl8AHAggh7ALQh8AHAQgh7ALQl8AHABPqGNWEPgC4EPgAonLAHQFcCHwAUTNgDoA+BDwAm0ia8nb16RdgDoDeBDwAa+tD5D/SeRpMQJ+gBMBSBDwAmti/QqdUDYGjXzV0AAFiCs1evxAsDTw8AxqaGDwA6uOncvXMXAQCOEvgAoKUhnuUDgCkIfADQ0wuXLs5dBADYSeADgBbU7gGwJAIfAABApQQ+AACASgl8AAAAlRL4AGAAL1y6qPMWAIoj8AEAAFRK4AOABi5ff8vcRQCA1gQ+AACASgl8AAAAlRL4AAAAKiXwAcCA9NQJQEkEPgAAgEoJfADQ0U3n7p27CABwkMAHAABQKYEPAACgUlUEPg/IAwAAvFYVgQ8AAIDXEvgAAAAqJfABAABUSuADgIF5thyAofQ9pwh8AAAAlRL4AAAAFuLy9be0Gl7gAwAAqJTABwAAUCmBDwAauPDoM3MXAQBaE/gAAAAqJfABAABUSuADgBF4Fx8AJRD4AAAAKiXwAbAKL1y6qNYNgNUR+ACo3umgJ/QBsCYCHwBVE/AAWDOBD4DVEQIBWAuBD4BqCXYArN1iA9/l62+ZuwgAcJDACcDcFhv4AOCQscOWMAfAEgh8AKzSUIHtpnP3DjIdABiDwAdAddS+AcA1Ah8ANPDgA0/MXQQAaK1tY6ZvAAATzklEQVSawOduLgARzgcAcFo1gQ+Asrxw6aLwFQIoAPMS+AAAACol8AFQjblq0246d6/eOgEoUq/Al1J6U0rpoZTS5c3/b9wxzDtSSr+bUno6pfTllNLPthkfAKZ09uqVuYsAAIPpW8N3X0Q8nHM+GxEPb37f9u2I+Ps55/8oIv7TiPiZlNKtLcYHAABYnSFarvQNfHdFxKXNz5ci4u7tAXLOz+ecn9j8/K2IeDoi3tZ0fABoostJse04x5ptatYJQGn6Br4bc87PR1wLdhHxlkMDp5TORMR7I+IPuowPQJ1q78my9uUDoFxHA19K6fMppad2/LurzYxSSt8TEZ+JiJ/LOf9524KmlO5JKT2WUnrs6y9/s+3oAIxkV5hpE3CEIQAYz3XHBsg5f2Tf31JKL6aUbs45P59SujkiXtoz3OvjWtj7jZzzb536U6PxN+W4PyLuj4i47d235mPlBmA95g6Nl6+/JSJ0+AJAefo26fxsRJzb/HwuIh7cHiCllCLi1yPi6ZzzL7cdH4ByzR20AIDD+ga+T0TEHSmlyxFxx+b3SCm9NaX0uc0w74+Ij0fEj6aUvrT5d+eh8QGY3lLD25TlPqnJm2PeANDF0Sadh+ScX46ID+/4/LmIuHPz86MRkdqMD8C0hgwu29N64dLF1fReefn6W+INe/62pvUAQDn61vABMCI1SABAHwIfQKHWEPb6LmOJ60gtHgAlEfgACjdlqJliXiWGNAColcAHUKEXLl0UrEZw4dFneo3vOwFgagIfMDsXwc2tYV01XcY1rAsA6EvgA4rg4h3bAAAMT+ADiqd5YpnG+E5Opnls2nNsDy8+9Uij4bbf1QcAcxL4AAoxVIgZOgw1DWE0Yz0CMCWBD2BgNV7Ql7RMc5blwQee8NoFABZF4ANYuTkCVNd5lhQ8DxEKAShFVYFvKRcCQDdj7eNjPos2hiGnXeKzcgDAcKoKfNCFC1pKsL0d2i6/o8R1sV2Dd/bqlVbjl7hMpbPOALoR+ACinIvJKctRyjIvjV44p2dbBehO4GPV1tz7YInLXGKZSnfslRVd1mmX6Y31svTTwy9p+xAKASiFwDeTJV24UB/b325DByeGU9L6v/DoM4NMp6RlYlpj3SAB2EXgm4ED+G7WC0tm+10XNXh05VgBTE3gowhOgJRiaT12MryhwpxXMwxjqc16j6lpWYCyCXxMxsmNXfo80zWVY896TtEUdNd0StynSizT0B584Im5i8AKDLEvrWF/BI4T+FauhJNBCWVYq6Yhwne0HDU9G1RaGW981wfjxaceGXSapS0j5bGNAH0JfCum6Vp/a1vetfC9csxQHbcwvTn3b8cWYA4C38TWcLA/VEN0bPnXsH5KZv2316c5p/UN/XTdh0rZ90opRxPHXgEDlEvgG9iQB8OlHViXVt6h1LjctXaSQFlq37aadv4y93qYcv5rPkcObYom+WtYx02WUdg9bMnrpvSyD1U+gW8EpW8825ZW3qH07Syk7wuol2bqZ/vG7gilaa1zk2l0GXbp7/xbQhnH8uJTj6yq45a1HetKUNKzuEOGnbm2Gdvq+Ja4jpdY5q4EvkI1veO0BmM0i1vCuiu9jCWXb6qLoKnnrea1XqV+n017qF1iDciSnmNvc00wxM2smvS9nmp7Y0/vqpRI4JvRWnboNjUbTaex/fnS1uXQIbZPTVNfY697HSw0c+hCY+nPOZWqyXv21vAuvrkvcIe8EVJL88e5v5PT0xjiGmCNalxPS1qmJZW1ieoCX21f0D6lL2ef5ihdT95jNoFpGiznqJmd4/mbkra/fUGnpDIespRyntY1WC5xWWs3Ze3FkM0CD+3zY94wazKNKdbRnCFqiHPtVOVdYhPSObbf09Maq7a+tOuHQ/vsGOt0btUFvlIsvb1723mXsDHvc3pnm+sk0ySQ7CpjmyBTYvOkrgfRMb6vKWq5St4PGMeQNXglbD9T1JCVtE+VsM5PlFSWtkooewll2KWUC/5jllDGtRnyOxH4ehqyKcmSlLCsc94F62KIu519g9BY66BpWG06zFBlODbMEu/+9p3v3PtBLY69gH1p7+lbys3FMec19vFxLG1C9hgtYeasWZ26DG23n0PnwRJuaA5dliWfX0qqEGg6bhtVB76xL26WeHIYMlgsoflWKRcSfabRZbpzrPuSTxglbItrt9bvoOmrGU6UcPyfaz6lbiNdriXmvql0yJJvdI15XTfU+XdJhq6R3zdOCetxiJsdbbe/kq6Tqw58x/S9IzKVNhtpCTtVE2OV9WS6Q4aksQ6IUxnjju6Y5egzvSG/q5K/0y7W2hphTG3DXEn61nRNsQ31ubBa0ja+fc5aws3kudfvofP83GU7ppTyzbmdbW/vpayTQ8a6qT+V1Qa+Ie++Tf0Fdj0xHFrmNtPsOmwthm4CMvSNh7Fr2kq5GJxyPqXMl3VY8/Y1xo2huW7ulnyTpW/g7zP9XcON0czt2DS7Xr80Pe+OuR6HvI4bS9vr7DFrxbenPfeN3yY3eYaed5WBb467BYd2yjYH1jE39KVb4nIM9f122aanumu8a76HyrHvM8ri+xlOl45dpqgRH7MMTY59JW1jc4XCUpTQwuPYxXjp21AT+5ah6Q35feMdmtcQ2/bp/8e8bh2qzF1vaAx1o3vOG037VBn4tjX58g9tzLuGaTrfQ79vf96n5qTtgbDtQWUOfec79d21qQ1xB7R2YwTLMU5Ea/1+aldCxy1z3kRschE75v7Q5wKs60V2l2m2GX/Mi+uuNxaHKse+n0s7jnedzhzzHGs+Q33vQ+5L+25yj3EtuWvf6Tuvk/KP9f1eN8pUC3VsJe7aUIbqcrvNvPf9fV9ZhryT2iXgTHlBMeWFcQmBc4p5ChvQzYVHn4lf/aHxpn/ouN91nH0X14eG62OJx5d9F2+7fh/ytRxNjRH6hlRC+bpeQ3S5VpqqJmeoSoEpXyUzxLrZLu8QlSRD6FsJNLWUc561AF3c9u5b8z976Pfj7NUrETH/SgRgubYvKE53xnJSS/fgA09ExKtfx3DP//Bz8as/9NcHxzt/+5k4e/VKr/NU0wu0XRd1J6HEeRJg+U6O7ZevvyXOXr0S173zvY/nnG87Nt4qmnQCwFz69uY5VxNCAOog8AFABye1flOYqtkYAPVZ1TN8ALBU2002p362GoBlUsMHACMbunMPgQ6ApgQ+AFZr7F4WS3g1AwDrJvABQEd9O2QBgLEJfAAAAJUS+ACgoRvf9cFXfj79Tj4AKJXABwAj0/QTgLkIfACwgw5XAKiBwAcAHV149Jm4fP0tavAAKJbABwAjOqkpHPsVEACwi8AHAABQKYEPgFVS4wbAGgh8ANDRgw888ZrPdPYCQEkEPgAYgKAHQIkEPgAYyKHQpydPAOZw3dwFAIBSdAllavYAKJkaPgBWp0+HLTe+64OtxxEKAZiLwAcAAFApgQ8AOnrxqUdaDe9VEABMTeADAAColMAHwKoMXcu26118AFAKgQ8AehL6ACiVwAcABwhzACyZwAcAA2gSDL18HYCpCXwArMZQz+91eRcfAMxB4AOACXj5OgB9dWkpIvABQA9t38UHAFMS+AAgpnu+zsvXAZiSwAfAKswRtPTwCcDcBD4A6OmkWaeAB0BpegW+lNKbUkoPpZQub/5/445h3pFS+t2U0tMppS+nlH721N9+IaX0tZTSlzb/7uxTHgCYyrGeOoU/AEpwXc/x74uIh3POn0gp3bf5/R9sDfPtiPj7OecnUkpviIjHU0oP5Zz/ZPP3X8k5/889ywEAg+nTo6agB0BJ+jbpvCsiLm1+vhQRd28PkHN+Puf8xObnb0XE0xHxtp7zBYDG5uwoRQAEYE59A9+NOefnI64Fu4h4y6GBU0pnIuK9EfEHpz4+n1J6MqX0yV1NQgGgJlP1BgoAEQ0CX0rp8ymlp3b8u6vNjFJK3xMRn4mIn8s5//nm44sR8f0R8Z6IeD4ifunA+PeklB5LKT329Ze/2WbWAHCQEAZArY4+w5dz/si+v6WUXkwp3Zxzfj6ldHNEvLRnuNfHtbD3Gznn3zo17RdPDfNrEfHbB8pxf0TcHxFx27tvzcfKDQBTevGpR4525AIAU+vbpPOzEXFu8/O5iHhwe4CUUoqIX4+Ip3POv7z1t5tP/frjEfFUz/IAQLH6dAYDAF30DXyfiIg7UkqXI+KOze+RUnprSulzm2HeHxEfj4gf3fH6hV9MKf1xSunJiPiRiPj5nuUBgFcZs8MWNXoAlK7Xaxlyzi9HxId3fP5cRNy5+fnRiEh7xv94n/kDAACwX98aPgBYtDk6bJnzNREArIvABwAj8y4+AOYi8AEAAFRK4AOgWlM3nXzxqUcmnR8AHCPwAbBau57fO/3qBE0xAVi6Xr10AsDa3HX3+175+cEHnogb3/VBNXsAFEsNHwBMbI6eQQFYJ4EPgFXqErpO1+4BwBIIfADQkQAIQOkEPgBWp2+Tygt/790DlQQAxiXwAVClKV/JcOO7PvjKz8c6cDndCygAjE3gA2BVutbunTTfPKndU8sHwBIIfACwofYNgNoIfADQQ9eOW6ZscgrAegl8AKzGmJ21nH6Ob5cHH3ii17wBoAuBDwA68hwfAKUT+ABgh9M1cm2bbR7rqRMApiLwAVCdXc/H9W3OCQBLJPABQPTrobNrxy0AMDaBDwAGcqzjFgCYmsAHQPWGaM65r4MWHbcAUDKBDwAOGKu5pmcKAZiCwAcAAFApgQ+A1evTYQsAlEzgA6BqUzWdPGn6qeMWAEoi8AEAAFRK4AOAI/TECcBSCXwAVOWmc/e2Gn7X83sPPvBEu2kIhAAUSuADgD26vpLBc3wAlELgAwAAqJTAB0C1jvXQOfTrGNrWCLZtfgoAbQl8AAAAlRL4AOAAHbIAsGQCHwCrdKw5Z9vmmYIhACUS+AAAACol8AEAAFRK4ANgdQ4152z70vVtXd/dBwBjEPgAqNKxVzIAwBoIfACwRS0dALUQ+ABgDz1vArB0Ah8A1bjp3L1Hhzn2OoY+DgVEtYYAzGHRgc/zGQCUSLgDoBSLDnwAMDRhDYCaCHwAAACVEvgAYCYeTQBgbAIfAKvRpsMWPXQCUAOBDwAGdBIUPQsIwFheuHSx8bACHwCMTPgDYC4CHwAAQKUEPgAYiZo9AOYm8AHAiIQ+AOa06MB39uqVuYsAQIX69tCph08ASrHowAcAu3i/HQBcI/ABAABUSuADgBnddO7euYsAQMUEPgCY0Pnbz8xdBABWROADgFOG6nBFxy0AlEDgAwAAqJTAB8BqaE4JwFK9cOlip/EWG/hO3sHXdcEBAABqt9jABwBDG/q5u+3pqWEEYGoCHwBVaPp6g32hSxgDoEbXzV0AAJjb+dvPxA/e+L2jTh8A+jp79Upcvv6WVuOo4QNgdU4HsLHDmNczADAnNXwArJJaNwCW6OzVK/FCi+HV8AGwemM255xi+gCwzyID3//3ja/PXQQAKjFVGBP6AJjDIgPfCe/gA6APIQyA2i068AFAV3OEvR+88XuFTAAmlXLOc5ehtZTStyLiK3OXA3Z4c0R8Y+5CwB62T0pl26RUtk1K9gM55zccG2ipvXR+Jed829yFgG0ppcdsm5TK9kmpbJuUyrZJyVJKjzUZTpNOAACASgl8AAAAlVpq4Lt/7gLAHrZNSmb7pFS2TUpl26RkjbbPRXbaAgAAwHFLreEDAADgiMUFvpTSR1NKX0kpfTWldN/c5YGIiJTSJ1NKL6WUnpq7LHBaSukdKaXfTSk9nVL6ckrpZ+cuE5xIKf27KaV/nVL6o832+Q/nLhOcllJ6XUrp/0op/fbcZYETKaVnUkp/nFL6UpOeOhfVpDOl9LqI+L8j4o6IeDYi/jAiPpZz/pNZC8bqpZQ+GBF/ERH/NOf8rrnLAydSSjdHxM055ydSSm+IiMcj4m7HTUqQUkoR8d05579IKb0+Ih6NiJ/NOf/+zEWDiIhIKf03EXFbRHxvzvnvzF0eiLgW+CLitpxzo3dELq2G74cj4qs55ys557+KiE9HxF0zlwki5/xIRPzZ3OWAbTnn53POT2x+/lZEPB0Rb5u3VHBNvuYvNr++fvNvOXeiqVpK6e0R8bcj4n+buyzQx9IC39si4k9P/f5suHABaCSldCYi3hsRfzBvSeA7Nk3mvhQRL0XEQzln2yel+F8i4r+NiL+euyCwJUfE/5FSejyldM+xgZcW+NKOz9wJBDgipfQ9EfGZiPi5nPOfz10eOJFz/rc55/dExNsj4odTSprFM7uU0t+JiJdyzo/PXRbY4f055/dFxI9FxM9sHi3aa2mB79mIeMep398eEc/NVBaARdg8G/WZiPiNnPNvzV0e2CXn/M2I+FcR8dGZiwIREe+PiL+7eVbq0xHxoymlfzZvkeCanPNzm/9fioj/Pa499rbX0gLfH0bE2ZTSO1NK3xURPxERn525TADF2nSK8esR8XTO+ZfnLg+cllK6IaX0fZufr4+Ij0TEv5m3VBCRc/7vcs5vzzmfiWvXm/9nzvm/nLlYECml7950whYppe+OiP8iIg72Er+owJdz/nZEnI+I34lrHQ/885zzl+ctFUSklH4zIn4vIn4gpfRsSumn5i4TbLw/Ij4e1+5Of2nz7865CwUbN0fE76aUnoxrN3Ufyjnr/h5gvxsj4tGU0h9FxL+OiH+Rc/6Xh0ZY1GsZAAAAaG5RNXwAAAA0J/ABAABUSuADAAColMAHAABQKYEPAACgUgIfAABApQQ+AACASgl8AAAAlfr/Acwd5ttTSkxWAAAAAElFTkSuQmCC\n",
  118.       "text/plain": [
  119.        "<Figure size 1080x648 with 1 Axes>"
  120.       ]
  121.      },
  122.      "metadata": {
  123.       "needs_background": "light"
  124.      },
  125.      "output_type": "display_data"
  126.     }
  127.    ],
  128.    "source": [
  129.     "visu.show_isw_propagation(croco_config)"
  130.    ]
  131.   },
  132.   {
  133.    "cell_type": "code",
  134.    "execution_count": 13,
  135.    "metadata": {},
  136.    "outputs": [],
  137.    "source": [
  138.     "from netCDF4 import Dataset"
  139.    ]
  140.   },
  141.   {
  142.    "cell_type": "code",
  143.    "execution_count": 16,
  144.    "metadata": {},
  145.    "outputs": [],
  146.    "source": [
  147.     "%matplotlib inline\n",
  148.     "import matplotlib.pyplot as plt"
  149.    ]
  150.   },
  151.   {
  152.    "cell_type": "code",
  153.    "execution_count": 14,
  154.    "metadata": {},
  155.    "outputs": [],
  156.    "source": [
  157.     "d = Dataset('../../configs/lab_waves_ISW/ISW_his.nc.no-obc', 'r')"
  158.    ]
  159.   },
  160.   {
  161.    "cell_type": "code",
  162.    "execution_count": 32,
  163.    "metadata": {},
  164.    "outputs": [
  165.     {
  166.      "data": {
  167.       "image/png": "iVBORw0KGgoAAAANSUhEUgAAAysAAAEWCAYAAABmLAgGAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJzt3XuYbHdd5/v3pzs75AoBAhJ2dghKuI0PCicCGm8joAkDBrwcQU2UA5PDCENyREfE0aOOzyiHMziHA5qJxAeQKIeZBI0ajRGjGOWSy4SEzRYJQck2W2OI5AZkp7u/549avVO7U9W9qruqq1b1+/U86+mqWr+11q+qVlX1d/2+v98vVYUkSZIkzZqFaVdAkiRJkgYxWJEkSZI0kwxWJEmSJM0kgxVJkiRJM8lgRZIkSdJMMliRJEmSNJMMViRpBiU5I8lnktyX5GXbdMyfT/K+TW775iTvmvRxmu1/NMly89o8Y5P7eHeSLyfZv9l6SJImz2BF0lxJ8nfNP6H39S1PnHa9NuEXgXdU1XFV9bvTrsxGquo/V9VrxrGv5j184QbFPtK8Nvs2c4yq+lHgrM1sK0naPgYrkubRS5t/ZFeX29cWSHLENCq2kb56PQnYO826SJI0bQYrknaEJKcmqSSvTvJ54M+ax5+f5K+TfDHJJ5J8e982j0pycZIDSf4hyS8lWRyy/yT51SR3JLk7yU1JvrZZ9+dJXtNX9keTXNN3v5K8LslngM8k+Szw1cDvNy1Dj0jyqiT7ktyb5NYk//ua45+d5MYk9yT5bJIzR30OjSOTvLc5zt4kp/cd44lJLk3yz0k+l+QNfesOS+1Kcm6Sv0/yhSQ/O6C1ZOBxkvwWcErfc/8P69S1//n/fJL/nuR9zT5vTvLUJD/dvCe3JfnONvuSJM0OgxVJO823Ac8AvivJbuAPgV8CHgP8BHBpksc1Zd8DLAFPAZ4NfCcwLNXpO4FvBZ4KnAD8APCFEer1MuB5wDOr6muAz/NQC9EDwB3AS4BHAq8CfjXJcwCSPBd4L/CTzbG/Ffi7TTwHgO8G3t/s53LgHc0xFoDfBz4B7AZeAFyQ5LvW7iDJM4FfA34IOAl4VLPNhsepqnPWPPf/a526rvVS4LeARwP/E7iS3u/cbnppdf9thH1JkmaAwYqkefS7TUvJF5Os7e/x81V1f1V9Gfhh4IqquqKqVqrqKuA64MVJvopen4YLmvJ3AL8KvGLIMR8EjgeeDqSq9lXVgRHq/MtVdVdTr4epqj+sqs9Wz18AfwJ8S7P61cBvVtVVzfP4h6r6m008B4Brmtdjmd4//l/XPP4NwOOq6her6mBV3Qr8xpB9fR/w+1V1TVUdBH4OqJbH2Yq/rKorq2oJ+O/A44BfqaoH6QVGpyY5YQzHkSRtk5nM2ZakLXpZVf3pkHW39d1+EvD9SV7a99gu4Opm3S7gQJLVdQur2yfZ25QBOKuq/izJO4B3Aqck+SDwE1V1T8s637beyiRnAf8nvZabBeAY4OZm9R7gigGbrfschvjHvttfAo5q+tE8CXhiki/2rV8E/nLAPp7Yf4yq+lKSta1MA4/TBBqb9U99t78M3NkEQ6v3AY4DvogkqRMMViTtNP1X+G8Dfquq/u3aQklOAh4AThz0D3RV/asBj70deHuSxwMfoJeW9bPA/fSCi1VP2KBea+vyCOBS4Fzg96rqwabFaDUCuQ34mgGb3rbecxjRbcDnquq0FmUPAE9bvZPkaOCxIxxr6GshSdpZTAOTtJO9D3hpku9KspjkqCTfnuTkJoXrT4D/kuSRSRaSfE2Sbxu0oyTfkOR5SXbRC06+Aqxe1b8R+J4kxyR5Cr20rVEcCTwC+GdgqWll6e8sfjHwqiQvaOq5O8nTR30OG/g4cE+Sn0pydPN6fW2SbxhQ9n/Qe12/KcmRwC/wUGDVxj/RG2BAkrTDGaxI2rGq6jbgbODN9AKB2+i1hqx+N55LL1D4FPAv9P4JP2nI7h5Jrw/HvwB/T69z/f/drPtV4CC9f8LfA1wyYj3vBd5Ar7XmX4AfpNcpfXX9x2k63QN3A3/BQylqozyH9eqwTK8D+9cDnwPuBN5Fr/P82rJ7gX9Pr5/IAeBeegMEPNDycL8M/Memz9FPjFpXSdL8SJWt7ZKkyUmy2k/ktKr63Bj2dw69kb0OAt+4mYkhk1wMfD9wR1U9Zat1kiRNhsGKJGnsmkELPkQv/eu/0BuW+Tnlj44kaQSmgUmSJuFs4PZmOQ14hYGKJGlUtqxIkiRJmkm2rEiSJEmaSZ2bZ2Xx+GPriMc+uhmFP/RujDIipiRNUsvW6q18bfVvWw9/PAvFwkKxfHCRo/55maVjFnnG7n/mb778aBY+8wAHn3gsRx3/AAeXj2B5aaG3j/J7dCQmJUjzYYOvvoOf339nVT1ueyqzed/1r4+tL9y1vHFB4PqbHriyqs6ccJXGpnPByhGPfTQn/cz5D/24pg7/kfUHRDtNf7y+NnZf/TwM+zJeL9YfdV9dt3rtY/X2pvez+WDlsHihfz+BDFvXv1EKFopdRy9x9FEHuecfHsnTL7yHu77uBD76lgv5pk98L8efdQt//7pv4mnfciufv/sE7r7rWHhgEVYy/Ptzvddj7TajfAd3+fva4E6aD4O+s/s+3n//Yz/599tXmc37wl3LfPzKU1qVXTzpMydOuDpj1blgRZIkSdJDClhhZdrVmAiDFUmSJKnDiuLBapcG1jXdD1Zsipfam1YqVxe6ltWA25upc3966nopYSO+JlVrUsFWrUkXg7C8tMD9X3oER9y7wMHHHsOXHh/uWfkyjz36flae9hR23QP/cO+jeODgLiC9qqap0KAqD6vrTkgB638f/b2RNMNsWZEkSZI0c4pieU6nIzFYkSRJkjpupTNN1qMxWJF2gkEje/WPILZqXFku8z5y2EZW04U2SgVb47DBDUfctrdN78/K0gJUOPKB8MCjj+CBx8AXV5Z45JFf4R9P2c3CQbj73qN7x2h7iGGjww27v079OqP/fZQ0nzbzXTuDClju3JdsOwYrkiRJUsfZsiJJkiRp5hTwoH1WJEmSJM2aokwD64T5fI+k6Rk0ZO3qY10Yjnir1vtOafvc18uH3uxruHaf/X0qqrfU0gIpeOCEBZaOX+HulSM49oiD3P+EXQRY/vIRsFCwkt6u2vTLGPZ6rPc6deV7eb3XVJJmXcFyV75vRzRfwYokSZK0w/RmsJ9PC9OugCRJkqStCMstlw33lJyZ5NNJbknypgHrk+TtzfqbkjyneXxPkquT7EuyN8n543hm3WtZmdMmLmlL1g5FvJUMls18xoYNVTyOes1DutkIqWCHDV88pMyhfa59aPUQKyEP9tYfPD7UsUt8YeUYjlxY4iuPDbUAfGWROnLloTSw1eOMquspYGtfR9O/JHVQr4P91r+/kiwC7wReBOwHrk1yeVV9qq/YWcBpzfI84Nebv0vAG6vqhiTHA9cnuWrNtiOzZUWSJEnqsN48K2NpWXkucEtV3VpVB4H3A2evKXM28N7q+ShwQpKTqupAVd0AUFX3AvuA3Vt9bgYrkiRJUsetVFotwIlJrutbzuvbzW7gtr77+3l4wLFhmSSnAs8GPrbV5zWxNLAke4D3Ak+g1+fnoqr6f9aU+Ungh/rq8gzgcVV117o7bztyjaSHtJlxfJTZzNvuc73t26aIzUsqGAxOBxv2/Nakj1VBBpSrNe9HVgIrvXSvpWMgRy7zxeVjWKQ4+ChY/Aq9NLGFvhnah72nm03dmzWtRjvr+kkmaadabVlp6c6qOn3IuoG/RqOUSXIccClwQVXd07ZSw0yyz8qGeWtV9VbgrQBJXgr8HxsGKpIkSZIOKcLyeBKm9gN7+u6fDNzetkySXfQClUuq6rJxVGhiaWCbyFt7JfA7k6qPJEmSNK9GSANbz7XAaUmenORI4BXA5WvKXA6c24wK9nzg7qo6kCTAxcC+qnrbuJ7XtowGtlHeWpJjgDOB1w9Zfx5wHsDiY06YSB0lDbHVlJ4acnvQunGNaLYd+us+yiSZ640Mtrq7/hHBNih/WApYNSN7NSN81SIsHQ1HHLnMXcvHsUx48LhmIsiVUMs5/DUfeIC+5zhs3UaPTYMpXZJ2kCIcrMWt76dqKcnrgSuBReA3q2pvktc26y8ErgBeDNwCfAl4VbP5GcA5wM1Jbmwee3NVXbGVOk08WGmZt/ZS4K+GpYBV1UXARQCPeNLJs/JTKEmSJE1db1LI8SRMNcHFFWseu7DvdgGvG7DdNUzgUuNEg5UR8tZegSlgkiRJ0qaM0MG+UyY5GlirvLUkjwK+DfjhSdVFkiRJmldVYbnmc0aSSbasDMxbA06Bw5qTXg78SVXdv6WjmRwmdV+bz/Es9mfp7zey2X4rG81k36yvtq8RUItFLcLi4gpfXD6GlQrLR6+w+OXFJmcgsNDyy3PQcNPDymwn+6ZIEgArM/fjOB4TC1ba5q1V1buBd0+qHpIkSdI863Ww35Zxs7bdfD4rSZIkaYcYZwf7WWOwImk2rU2pGja08ThmXB+nMaaCtd736sz1xWGvx8oRwALsSnH30jEs1SJ11AorRy4eKr9uWlnb4YrXe3wzTO2SpJEtz+l3p8GKJEmS1GFjnMF+5hisSJIkSR234mhgkqaqfzSmQelA8zgi3kYz3s+q9WZ9P6zcxrPZty63pvm/AlksVhYhKe5ZOooHlo8gu1ZY2VUPnUfr7bt1SluLMi3qLEnanAJbViRJkiTNniI8WIvTrsZEGKxIkiRJHVaFk0JKmrIMuT3ILE6cuF1maXSwNu/D2jSsvm0eNjHksGMMkqIWAkcUCdy//AiWaoHFXSusHFGHZ35V3wFHfQ6jpICZ9iWpCzr5GxonhZQkSZI0ewpbViRJkiTNKDvYS5IkSZo5RViZ01Tb+QhWZnkIU2mc+s/1+fxOGr+1Qz7Td3+7jj/qsQb1W1nt29J2f4FaKGqx96TvXzqSgytHsLC4wvIRBUsZMPz1kL4rw47Z5rt3Tn88JWmWFPBgzce/9WvN57OSJEmSdoywPKdXMQ1WJEmSpA4rnMFe0iwY5aKJ6ZGHW/t6bOfQlBsO/9tyJvtB+x0kfX8Xiir4yvIullYWyELBYlErGX7IQfXZVDrbfF7lk6RZZMuKJEmSpJlTlbltWZnYs0qyJ8nVSfYl2Zvk/CHlvj3JjU2Zv5hUfSRJkqR51Otgv9hq6ZpJtqwsAW+sqhuSHA9cn+SqqvrUaoEkJwC/BpxZVZ9P8vhWeza1QNI49I8Uth3HGiUVbKtpagFSvbQv4CvLR7BSYSG9EcKyVNRWDrBe1prf0ZK0zeKkkKOqqgPAgeb2vUn2AbuBT/UV+0Hgsqr6fFPujknVR5IkSZpHvQ7283mhaFtCsCSnAs8GPrZm1VOBRyf58yTXJzl3yPbnJbkuyXXL990/2cpKkiRJHbPMQqulaybewT7JccClwAVVdc+A4/8vwAuAo4GPJPloVf1tf6Gqugi4COART9rjGEeSxms15WrSaWGbTO0aODFkqw17S1U4uNzLU06KLNC7VLVSo6VsrX2dBpaZzyt7knaYjv236Qz2m5RkF71A5ZKqumxAkf3AnVV1P3B/kg8DXwf87YCykiRJkgZY6WCrSRuTHA0swMXAvqp625Bivwd8S5IjkhwDPA/YN6k6SZIkSfOmCh5cWWi1dM0kW1bOAM4Bbk5yY/PYm4FTAKrqwqral+SPgZuAFeBdVfXJCdZJkiRJmiu9NLDuBSJtTHI0sGtokZ1dVW8F3rq5g2xqK0l6uO36Plmv38pmZ7Lvl9WDNLtr+qw8uNI3tv5CHerPApvotyJJmjnOYC9JkiRp5szz0MUGK5IkSVKnmQY2e0xFkDQpszKE8VZmsW9msF99Mg8uL7KwmmKW6qWCVZpjjJgKNrS+83lVT5K6YMU0MEmSJEmzpjca2OLGBTvIYEWSJEnqMCeFlKSdaCtpWJs+5gYjgq2z/tCq1To3fwtYXlmgDqWBPTRS2KEZ6beaCjanP5KS1BXzmgY2nz1xJEmSpB1idTSwNstGkpyZ5NNJbknypgHrk+Ttzfqbkjynb91vJrkjydjmTTRYkSRJkjpupRZaLetJsgi8EzgLeCbwyiTPXFPsLOC0ZjkP+PW+de8GzhzTUwIMViRpfQ8NqDX+/Y5QbtS5Ig+ledGbFHK5wvJK74Gkejtcu9OtTkgpSZqKqrBUC62WDTwXuKWqbq2qg8D7gbPXlDkbeG/1fBQ4IclJvXrUh4G7xvncDFYkSZKkjhshDezEJNf1Lef17WY3cFvf/f3NY4xYZmzsYC9JkiR12Igz2N9ZVacPWTdoJ2ub3duUGRuDFUmSJKnjxjR08X5gT9/9k4HbN1FmbEwDk6Q2JtF3Zdj+NvrBGeUHqRmaeHllgeXVXOXVIYvX3obN9Vtx2GJJmqrVeVbGMBrYtcBpSZ6c5EjgFcDla8pcDpzbjAr2fODuqjow/mfVY7AiSZIkddwKabWsp6qWgNcDVwL7gA9U1d4kr03y2qbYFcCtwC3AbwA/trp9kt8BPgI8Lcn+JK/e6vMyDUySJEnqsCpYWhlPG0RVXUEvIOl/7MK+2wW8bsi2rxxLJfp0M1hxdE1J0zKNWe1HcKiFvy+dqyqsrAQWerehb/b6QRzCWJI6Z0x9VmZON4MVSZIkScBDfVbm0cT6rCTZk+TqJPuS7E1y/oAy357k7iQ3NsvPTao+kiRJ0ryqSqulaybZsrIEvLGqbkhyPHB9kquq6lNryv1lVb1kgvWQpPEaZyrYJNPKHpYK1jyc6g1utjrLvVlfktR5G3We76qJBSvNEGYHmtv3JtlHb3bLtcGKJEmSpE2qmt8+K9sydHGSU4FnAx8bsPobk3wiyR8l+VdDtj8vyXVJrlu+774J1lSSJEnqmvTm02qxdM3EO9gnOQ64FLigqu5Zs/oG4ElVdV+SFwO/C5y2dh9VdRFwEcAjnrTHhAWpX/+FFD8d86PyUBpXkyqW2tz8i1WwvBIWF/pOkfTln63uczPnz5xeyZO0Q8zR6Idd7I/SxkTDqyS76AUql1TVZWvXV9U9VXVfc/sKYFeSEydZJ0mSJGmeFIxrBvuZM8nRwAJcDOyrqrcNKfOEphxJntvU5wuTqpMkSZI0d6rXit5m6ZpJpoGdAZwD3JzkxuaxNwOnwKGZML8P+HdJloAvA69oZsWUJEmS1JKjgY2oqq5hgwE5q+odwDsmVQdJmpiuDF/MQ2Pvr9SanOZU77gdTAuQpInp7zPYEdV0sJ9HzmAvSZIkddy85iYZrEiSJEkdN6+jgRmsSNJmTTh9ayxWRz/um8EeGFxvZ7OXtFN0LM1rI7U2zXeOGKxIkiRJHdfFYYnbMFiRJEmSOs4+K5Jm32oaj+k822dcqWCD9rPeLPYbjVbTv69K73SoHH5epG/xfJE0zLRSpsbdUjDK8+jYiGBFWHE0MEmSJEmzqDuh1WgMViRJkqQus4O9pM6Yz++q2TZro4KtSV04NErMRmlj83pZTppnHUpVGtk8P7dJmNOXy2BFkiRJ6jhbViRJkiTNnAJWVgxWJHWJaT3ba9ZSwdaq3mgxh668DavrKOfNnF7Fk7aVqU4ah2Juv5MNViRJkqSOc54VSZIkSbPJYEWSJEnS7MnO7GCf5J4Ntg9woKqeOr4qSdIOtdEs9lvYb1UIdViaQObzd03afqnOzXiuOTSnp9/CBus/W1WPXGc5Hrh/0IZJ9iS5Osm+JHuTnD/sIEm+Iclyku/bypORJEmSdpyCWkmrpWs2SgP73hb7GFZmCXhjVd2Q5Hjg+iRXVdWn+gslWQTeAlzZ4liSJEmSHqZ7gUgb67asVNWtAEmOTbLQ3H5qku9Osqu/zIBtD1TVDc3te4F9wO4BRf89cClwx6afhSTNgm1qgh8t06T343VoFvs1KiPvUNp5Vj8jqcFLfxlpWqrl0jEbpYGt+jBwVJLdwIeAVwHvbnuQJKcCzwY+tubx3cDLgQs32P68JNcluW75vvvaHlaSJEnaGXZ4sJKq+hLwPcD/W1UvB57ZasPkOHotJxdU1doO+/8V+KmqWl5vH1V1UVWdXlWnLx53XMsqS5IkSTvA6qSQbZaOaTt0cZJ8I/BDwKvbbtukil0KXFJVlw0ocjrw/vSGpDkReHGSpar63Zb1krQeZ7GfL2tHC2s5+lBVb0jLpA7dHmj1fPG80U7SNn3LNC/NuJ0+KeT5wE8DH6yqvUm+Grh6vQ3Si0AuBvZV1dsGlamqJ/eVfzfwBwYqkiRJ0og6ONJXG62Clar6ML1+K6v3bwXesMFmZwDnADcnubF57M3AKc0+1u2nIkmSJKmdeW38W7fPSpKf32gHw8pU1TVVlap6VlV9fbNcUVUXDgpUqupHq+p/tK24pJbCvI5mOJu2+mMxaPv1cozb5iBX5jZFQBpq2Ohda0fykrqubef6Fqd8kjOTfDrJLUneNGB9kry9WX9Tkue03XYzNmpZec0Gs9gHeAXw8+OojCRJkqRRjafzfDP/4TuBFwH7gWuTXL5mnsSzgNOa5XnArwPPa7ntyDYKVn4DOL5FGUmSJEnTMp6GwucCt/TNtfh+4GygP+A4G3hvVRXw0SQnJDkJOLXFtiNbN1ipql/Yys4lSZIkbYOV1iVPTHJd3/2Lquqi5vZu4La+dfvptZ70G1Rmd8ttR9Z2NDBJXedwtDvL6pj71KH3vdbr3+L5oS6yz4nUc+g7v5U7q+r0IesG7WTtB21YmTbbjsxgRZIkSeq4McXu+4E9ffdPBm5vWebIFtuOrO0M9pIkSZJm1XhGA7sWOC3Jk5McSW8grcvXlLkcOLcZFez5wN1VdaDltiNrFawkeWqSDyX5ZHP/WUn+41YPLmmbOYTx9pjE8MVr1rW+gnYoBYzRRorxXNEscdhhbacdfG5V1RLweuBKYB/wgWZC+NcmeW1T7ArgVuAWegNt/dh62261Tm3TwH4D+EngvzWVuSnJbwO/tNUKSJIkSdqaccVYVXUFvYCk/7EL+24X8Lq2225V22DlmKr6eHLYpbalcVZEkiRJ0iYUsDKfTeJtg5U7k3wNTUJBku8DDkysVpImx1GftsewcVG25dg57BJbVUjzplfLGYylbbeDU2+ksZjTj1DbYOV1wEXA05P8A/A54IcnVitJkiRJrc1rvN8qWGlmonxhkmOBhaq6d7LVkiRJktTaTg5WkpwAnAucChyx2nelqt4wsZpJmhxTwbpnTWrX6NsPecuzTr6a54kmYV4v/0rTNqcfrbZpYFcAHwVuBlYmVx1JkiRJo5jn0bzbBitHVdWPT7QmkiRJkjZnh48G9ltJ/i3wB8ADqw9W1V0TqZUkSZKk1ua1ZaXVDPbAQeCtwEeA65vluvU2SLInydVJ9iXZm+T8AWXOTnJTkhuTXJfkm0d9ApI2aT4vwMyWrfxwtNi21Q/ToRnsM3y/GXJb2ipnnZe2T7VcOqZty8qPA0+pqjtH2PcS8MaquiHJ8cD1Sa6qqk/1lfkQcHlVVZJnAR8Anj7CMSRJkqSdbY6vBbRtWdkLfGmUHVfVgaq6obl9L7AP2L2mzH1VtfrSHksn4z1JkiRpynZ4y8oycGOSqzm8z0qroYuTnAo8G/jYgHUvB34ZeDzwb4Zsfx5wHsDiY05oWWVJG1pN+engl9eOts5ow0PL9/9de1vaqnm9pKudq4PndOZ0vN62wcrvNsvIkhwHXApcUFX3rF1fVR8EPpjkW4H/BLxwQJmLgIsAHvGkPd07eyRJkiSNrO0M9u/ZzM6T7KIXqFxSVZdtcIwPJ/maJCeO2DdGkiRJ2tnm9HL+usFKkg9U1f+a5GYGvARV9ax1tg1wMbCvqt42pMxTgM82HeyfAxwJfGGUJyBpDJypfC5VQfpT/dIbFezQyGBrDToPPDe0VgfTY6S5N8cd7DdqWVkdbvglm9j3GcA5wM1JbmweezNwCkBVXQh8L3BukgeBLwM/0NfhXpIkSVIbc/of9LrBSlUdaG7+WFX9VP+6JG8BfurhWx3a9ho26AJaVW8B3tKuqpIkSZIGmtNgpe3QxS8a8NhZ46yIpClzMsDJGOfEkMPStzbU225o+pc0jJM6Sp0QeqOBtVm6ZqM+K/8O+DHgq5Pc1LfqeOCvJlkxSZIkSS3M8bWEjfqs/DbwR/TmQXlT3+P3VtVdE6uVJEmSpPZ2YrBSVXcDdwOv3J7qSJIkSRrZTgxWJO0wDlM7GaPOOD+CVLuuLFUhKd9frW9e80ikHWBeP74GK5IkSVLXGaxIkiRJmjnVzZG+2jBYkXS4/hnPNX3DUsjappZVX27f6ns66L1tm0+m+TKveSPSTjSnH2eDFUmSJKnj5vXag8GKJEmS1HUGK5J2FNPBxmtcI4JV2l0+6y83wdHI1CHzetlVUu97fk4/4gYrkiRJUoeF+b0eYbAiSZIkdZzBiqSdyYkiZ97qQF4PG9CrP/2roNrmgrV9zx1BbLbN638ukgab04+8wYokSZLUdQYrkiRJkmZOzW9jqsGKJEmS1HVzGqwsTGrHSfYkuTrJviR7k5w/oMwPJbmpWf46yddNqj6StsBuCeOx2R+SFtutXlGb1ytraiF1+CJpR8lKu2VLx0gek+SqJJ9p/j56SLkzk3w6yS1J3tT3+Pc3ccFKktPbHHNiwQqwBLyxqp4BPB94XZJnrinzOeDbqupZwH8CLppgfSRJkqS5tPZ6xbBli94EfKiqTgM+1Nw/vB7JIvBO4CzgmcAr+2KATwLfA3y47QEnFqxU1YGquqG5fS+wD9i9psxfV9W/NHc/Cpw8qfpIkiRJc6lGWLbmbOA9ze33AC8bUOa5wC1VdWtVHQTe32xHVe2rqk+PcsBJtqwckuRU4NnAx9Yp9mrgj4Zsf16S65Jct3zffeOvoKSNBdPBZsXqcMFDfnT6r5zVoR+n9JXP3OY27yimfUnq1z5YOXH1/+pmOW+Eo3xVVR2AXsME8PgBZXYDt/Xd38+aBotRTLyDfZLjgEuBC6rqniFl/jW9YOWbB62vqotoUsQe8aQ9fiNLkiRJjRFnsL+zqob2F0nyp8ATBqz6mRGqs9am/3+faLCSZBe9QOWSqrpsSJlnAe8CzqqqL0yyPpIkSdIti8upAAANb0lEQVQ8ysp4rudX1QuHHiP5pyQnVdWBJCcBdwwoth/Y03f/ZOD2zdZnkqOBBbgY2FdVbxtS5hTgMuCcqvrbSdVF0hiZDrY1024b3szxfb9ng2lfkobZvj4rlwM/0tz+EeD3BpS5FjgtyZOTHAm8otluUybZZ+UM4BzgO5Lc2CwvTvLaJK9tyvwc8Fjg15r1102wPpIkSdJc2qbRwH4FeFGSzwAvau6T5IlJrgCoqiXg9cCV9AbY+kBV7W3KvTzJfuAbgT9McuVGB5xYGlhVXcMG1+Oq6jXAayZVB0mSJGlH2IYG16bLxgsGPH478OK++1cAVwwo90Hgg6Mcs5sz2BemJUjaWYZ97/l9OL9M9ZI0gnn9yuhmsCJJkiTpIQYrkiRJkmZOQVamXYnJMFiRJEmSOmzEeVY6xWBF0ub095OY0y/IiZlGP5MCEvrfrPJ9mz3z+t+GpMmb0y91gxVJkiSp4+b1WofBiiRJktRl45nwcSYZrEjautWUpjn9opxJlY0vo62WOaxsfJ9mzbxeDpW0rexgL0mSJGkmGaxIkiRJmj2FHewlSWO0mRHBRthmNftr4iOPHTqQWjHlS9KEzOvXi8GKJEmS1HUGK5IkSZJmjZNCzqJ5fEPMpFDXDTqH5/GzOi7jStFqu5/+cqupW6ZwbY95/S9C0myoIivz+T3T3WBFkiRJUs98xioGK5IkSVLXzWsD7sKkdpxkT5Krk+xLsjfJ+QPKPD3JR5I8kOQnJlUXSZIkaW4VsFLtlo6ZZMvKEvDGqrohyfHA9UmuqqpP9ZW5C3gD8LIJ1qM7pn3+mLouzb7N9nNZ3W7a3zM7ybxe5pQ0m+b0K2diLStVdaCqbmhu3wvsA3avKXNHVV0LPDipekiSJEnzLtVu6Zpt6bOS5FTg2cDHNrn9ecB5AIuPOWFs9ZIkSZLmwbyOBjaxlpVVSY4DLgUuqKp7NrOPqrqoqk6vqtMXjztuvBXUQ2oCi5S+RQ+3lc9Jy2GHa+0x/GxORtcvX0rqrjn+32yiLStJdtELVC6pqssmeSxJkiRpJ+pNCtnBSKSFiQUrSQJcDOyrqrdN6jiSJEnSjrcy7QpMxiRbVs4AzgFuTnJj89ibgVMAqurCJE8ArgMeCawkuQB45mbTxTSDxhHkmz40Pw7Nnj7VWsyeccxk3/+abvdnJtU6JW3umO4laUbYsjKiqrqGDX4yq+ofgZMnVQdJkiRp7nW0P0obzmAvSZIkdVrN7WhgBiuafZv57O3QjJTOWPv+zOf362jGkQo2aF+VvjSt8rMxDqZ+SZpFpoFJkiRJmjkFsYO9JEmSpJlky4okSZKkmTSfsYrBiubUKB9Yc/inr+lOoZYG9UkZYOCIwpN6ned9+GL7qUiacVmZzzwwgxVJkiSpywonhZQkSZI0e0I5KaQ0t0wZmw0OZzze4Yu3Y7/zztQvSV1isCJJkiRpJs1psLIw7QpIkiRJ2oLVPittli1I8pgkVyX5TPP30UPKnZnk00luSfKmvsffmuRvktyU5INJTtjomAYr0ihqg0Xjk75lJ9nMeTRoFK5B+1ktV+ldgBv3OTsPaVOphxZJ6pCsrLRatuhNwIeq6jTgQ839w+uRLALvBM4Cngm8Mskzm9VXAV9bVc8C/hb46Y0OaLAiSZIkdVr10sDaLFtzNvCe5vZ7gJcNKPNc4JaqurWqDgLvb7ajqv6kqpaach8FTt7ogAYrkiRJUpcVowQrJya5rm85b4QjfVVVHQBo/j5+QJndwG199/c3j631vwF/tNEB7WAvjdNGFyx2WkrTuOy0SSPbjN414ghfVRDPv8FM+ZI0D9pneN1ZVacPW5nkT4EnDFj1My33P+jX5rAv2iQ/AywBl2y0M4MVSZIkqePGNc9KVb1w6DGSf0pyUlUdSHIScMeAYvuBPX33TwZu79vHjwAvAV5QtXGlJ5YGlmRPkquT7EuyN8n5A8okydubkQJuSvKcSdVHkiRJmlvb02flcuBHmts/AvzegDLXAqcleXKSI4FXNNuR5Ezgp4DvrqovtTngJPusLAFvrKpnAM8HXtc3EsCqs4DTmuU84NcnWB9JkiRp/lTB8kq7ZWt+BXhRks8AL2ruk+SJSa7oVaWWgNcDVwL7gA9U1d5m+3cAxwNXJbkxyYUbHXBiaWBNp5vVDjj3JtlHr3PNp/qKnQ28t2kC+miSE1abliZVL2mq7NOyecNeG7sbjGbQMMfjlJr8McbBfiqS5s02TApZVV8AXjDg8duBF/fdvwK4YkC5p4x6zG0ZDSzJqcCzgY+tWdV2tABJkiRJw2xPGti2m3gH+yTHAZcCF1TVPWtXD9jkYa9iM6TaeQCLj9lwoktJkiRp5yhgpXuBSBsTbVlJsoteoHJJVV02oMi6owWsqqqLqur0qjp98bjjJlNZaRbUmkU715jef7Od+jg7vaS5VVAr7ZaOmeRoYAEuBvZV1duGFLscOLcZFez5wN32V5EkSZJGUGxXB/ttN8k0sDOAc4Cbk9zYPPZm4BSAqrqQXsebFwO3AF8CXjXB+kiSJEnzqYP9UdqY5Ghg17DB2EbNKGCvm1QdpM5r+73TgcGXJmaeRwnbaJb6/vWV2UlvmqURwWblNZGkSTNYkSRJkjR7ujnSVxsGK5IkSVKXFbDSvf4obRisSPNgo3ShnWjt69HVC05beW9Xt21SxKog23WeTDMVzNQvSTuRLSuSJEmSZk91cqSvNgxWJEmSpC4rqA7OodKGwYokSZLUdXM6g73BijQvHOZ4fYeG+J1qLSZrtY/IoD4bfX1fti2teTv7rdhPRdJOZ58VSZIkSTOnytHAJEmSJM0oW1YkaQ50MR1sveGLB61bZzb7Q79l85AKZuqXJDWKWl6ediUmwmBFkiRJ6rLCDvaSJEmSZpRDF0uaC8MuvOy0UcL6n28XLkZtZSb7daSY/PNfTdfaajqYaV+SNFABZcuKJEmSpJlTZcuKJEmSpNk0rx3sUx0b5izJvcCnp10PzZwTgTunXQnNHM8LDeJ5oUE8LzTI06rq+GlXYiNJ/pjeOdzGnVV15iTrM05dDFauq6rTp10PzRbPCw3ieaFBPC80iOeFBvG8mL6FaVdAkiRJkgYxWJEkSZI0k7oYrFw07QpoJnleaBDPCw3ieaFBPC80iOfFlHWuz4okSZKknaGLLSuSJEmSdgCDFUmSJEkzqVPBSpIzk3w6yS1J3jTt+mj6kvxmkjuSfHLaddHsSLInydVJ9iXZm+T8addJ05XkqCQfT/KJ5pz4hWnXSbMjyWKS/5nkD6ZdF82GJH+X5OYkNya5btr12ck602clySLwt8CLgP3AtcArq+pTU62YpirJtwL3Ae+tqq+ddn00G5KcBJxUVTckOR64HniZ3xc7V5IAx1bVfUl2AdcA51fVR6dcNc2AJD8OnA48sqpeMu36aPqS/B1welU5UeiUdall5bnALVV1a1UdBN4PnD3lOmnKqurDwF3TrodmS1UdqKobmtv3AvuA3dOtlaapeu5r7u5qlm5crdNEJTkZ+DfAu6ZdF0kP16VgZTdwW9/9/fjPh6QNJDkVeDbwsenWRNPWpPrcCNwBXFVVnhMC+K/AfwBWpl0RzZQC/iTJ9UnOm3ZldrIuBSsZ8JhXxSQNleQ44FLggqq6Z9r10XRV1XJVfT1wMvDcJKaO7nBJXgLcUVXXT7sumjlnVNVzgLOA1zVp55qCLgUr+4E9ffdPBm6fUl0kzbimX8KlwCVVddm066PZUVVfBP4cOHPKVdH0nQF8d9M/4f3AdyR533SrpFlQVbc3f+8APkivO4KmoEvByrXAaUmenORI4BXA5VOuk6QZ1HSmvhjYV1Vvm3Z9NH1JHpfkhOb20cALgb+Zbq00bVX101V1clWdSu//ij+rqh+ecrU0ZUmObQZnIcmxwHcCjjo6JZ0JVqpqCXg9cCW9zrIfqKq9062Vpi3J7wAfAZ6WZH+SV0+7TpoJZwDn0LtKemOzvHjaldJUnQRcneQmehe/rqoqh6mVNMhXAdck+QTwceAPq+qPp1ynHaszQxdLkiRJ2lk607IiSZIkaWcxWJEkSZI0kwxWJEmSJM0kgxVJkiRJM8lgRZIkSdJMMliRJEmSNJMMViRpDiU5NcmXk9w44nY/kOSWJM5BIkmaOoMVSZpfn62qrx9lg6r6/4DXTKg+kiSNxGBFkjomyTckuSnJUUmOTbI3yddusM2pSf4mybuSfDLJJUlemOSvknwmyXO3q/6SJLV1xLQrIEkaTVVdm+Ry4JeAo4H3VdUnW2z6FOD7gfOAa4EfBL4Z+G7gzcDLJlNjSZI2x2BFkrrpF+kFHF8B3tBym89V1c0ASfYCH6qqSnIzcOpEailJ0haYBiZJ3fQY4DjgeOColts80Hd7pe/+Cl68kiTNIIMVSeqmi4CfBS4B3jLlukiSNBFeSZOkjklyLrBUVb+dZBH46yTfUVV/Nu26SZI0TqmqaddBkjRmSU4F/qCq1h0lbMi23w78RFW9ZMzVkiRpJKaBSdJ8WgYetZlJIYFfA/5lIrWSJGkEtqxIkiRJmkm2rEiSJEmaSQYrkiRJkmaSwYokSZKkmWSwIkmSJGkm/f/WGtqOZ5zI4gAAAABJRU5ErkJggg==\n",
  168.       "text/plain": [
  169.        "<Figure size 1080x288 with 2 Axes>"
  170.       ]
  171.      },
  172.      "metadata": {
  173.       "needs_background": "light"
  174.      },
  175.      "output_type": "display_data"
  176.     }
  177.    ],
  178.    "source": [
  179.     "plt.figure(figsize=(15,4))\n",
  180.     "plt.pcolormesh(d.variables['x_rho'][2,:], d.variables['time'][150:], d.variables['zeta'][150:,2,:])\n",
  181.     "plt.xlabel('x [m]')\n",
  182.     "plt.ylabel('time [s]')\n",
  183.     "plt.title('Free-surface height [m]')\n",
  184.     "plt.colorbar();"
  185.    ]
  186.   }
  187.  ],
  188.  "metadata": {
  189.   "kernelspec": {
  190.    "display_name": "Python 3",
  191.    "language": "python",
  192.    "name": "python3"
  193.   },
  194.   "language_info": {
  195.    "codemirror_mode": {
  196.     "name": "ipython",
  197.     "version": 3
  198.    },
  199.    "file_extension": ".py",
  200.    "mimetype": "text/x-python",
  201.    "name": "python",
  202.    "nbconvert_exporter": "python",
  203.    "pygments_lexer": "ipython3",
  204.    "version": "3.7.3"
  205.   }
  206.  },
  207.  "nbformat": 4,
  208.  "nbformat_minor": 2
  209. }
RAW Paste Data
We use cookies for various purposes including analytics. By continuing to use Pastebin, you agree to our use of cookies as described in the Cookies Policy. OK, I Understand
Not a member of Pastebin yet?
Sign Up, it unlocks many cool features!
 
Top