{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Read a wave2d netCDF output" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### First, import the script" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import netCDF_Reader as nrd" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Now, build your Wave object with the filename" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "filename = 'history.nc'\n", "wave = nrd.import_Wave(filename)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "#### What is a Wave object ?" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "wave?" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "#### Now, we want to know the names of variables that we can plot" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Data variables:\n", " p (time, y, x) float64 ...\n" ] } ], "source": [ "wave.variable_list()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We can choose our variable" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "variable = 'p'" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Now, we want to build the animation" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "How works the method animate?" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "wave.animate?" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "So, we just have to use the method with the variable that we choose" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "MovieWriter PillowWriter unavailable; trying to use instead.\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "...........Done!\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAArQAAAE/CAYAAAC3oVBcAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAgAElEQVR4nO3de5RkVX33//e3untmGDA4OIjcL8mESICAjBAXEgGVAPEHuWgCxgTz04cVFU18nlzgMQv9kZgQ8zyaaLxNDIK5QNTEZGLGICoqiWIYFEUgwogGhplwGwSGmenp7vr+/th7nzp1uk5duqu661R9XmvVOnWutetUddeub33Pd5u7IyIiIiJSVbXlboCIiIiIyGKoQysiIiIilaYOrYiIiIhUmjq0IiIiIlJp6tCKiIiISKWpQysiIiIilaYOrYiIiMgYMrNrzOwRM/t2yXozs/ea2RYz+5aZvSC37hIzuy/eLlm6VremDq2IiIjIeLoWOLfN+vOAdfF2KfBBADM7AHg7cBpwKvB2M1sz0JZ2MLAOrZkdbmY3m9k9ZnaXmf1Gi20q0/MXERERGSXu/mVgR5tNLgQ+5sGtwLPN7GDgp4Gb3H2Huz8B3ET7jvHATQ7w2LPA/3L3r5vZs4Dbzewmd787t02+538aoed/Wq7nvx7wuO/GeNJEREREZPAOBR7MzW+Ny8qWL5uBdWjdfTuwPd5/2szuITzZfIc26/kDt5pZ6vmfSez5A5hZ6vlfP6j2ioiIiAyrU2r7+lM+19M+W5i+C9iTW7TB3Tf0cAhrsczbLF82g4zQZszsKOBk4GuFVZXp+YuIiIgsl6d8jj+dPLKnfV4xe+8ed1+/iIfdChyemz8M2BaXn1lY/sVFPM6iDbxDa2b7AX8P/Ka7P1Vc3WKXnnr+ZnYpIVGZ1atXn/LDxxy9iNbKsvNl/YLXPWv1FhURb/nvW5aKLW+QbCRtfeghdux4Yvnf2AY21WMzZhf9qBuBy8zsBkJq6JPuvt3MbgT+MHch2DnAFYt+tEUYaIfWzKYIndm/cfd/aLHJonv+MXS+AeDEE473jZ+a/zDm9d4bP0TG6R9UVV4rt/EpEFL1Dso4vVbDoOrvl6obp8+LpXLBz/38cjcBADOjNtnfvy8zu57Q31prZlsJ1y9NAbj7h4BNwPnAFmAX8Gtx3Q4z+33gtnioq1Ka6HIZWIfWzAz4S+Aed393yWaV6fnLYFWlI5uk9qqzJCIiS8LApvr7mePuF3dY78CbStZdA1zT1wYtwiAjtKcDvwLcaWZ3xGX/GzgCqtfzFxEREVk2Rt8jtKNkkFUO/o3WubD5bSrT8xcRERFZNgvJoR0jS1LlQEREREQWbhA5tKNEHVoRERGRYacIbVvq0IqIiIgMO+XQtqUOrSyrqlU3KFK1AxERWQoG2IQ6tGXUoRUREREZdgY1dWhLjXSHtvLRPxXIliGQ3odVLZivKLqIjAbDatX8P7wURrpDKyIiIjISDGxCX8zLqEMrIiIiMuQMpRy0ow6tLIuqp4MU6WdtEREZKEMpB22oQysiIiIy9EwR2jbUoZUlNWqR2SJFakVEZBDMVLarHX3qioiIiEilKUIrIiJ9U/Uyb1WlMo/jwWqKQ5ZRh1ZERERk2OmisLbUoR1Co/hNe9RzZ4tGMZdWkTcRkeWki8LaUYdWREREZMiZIrRtqUMrAzNuUdlWRjFSKyIiy0M5tOXUoRUREREZdorQtjWSHVpFBpeXzv98+XOiaO3yULRcRKpNObTtjGSHVkRERGSUKIe2PXVoh8goVjeQ0aNqByIiy0M5tOXUoRUREREZdorQtqUOrfSNcme7o1xOERHpnalD24Y+UUVEREQqwGrW062rY5qda2bfMbMtZnZ5i/XvMbM74u1eM/tBbt1cbt3GPj7VnilCOwSqnjuryOzCVD1Sq1xaEZGlEy4K6+/nhZlNAO8HXg5sBW4zs43ufnfaxt3fmtv+zcDJuUPsdveT+tqoBRrYJ6mZXWNmj5jZt0vW/3auV//t2Ms/IK77vpndGddtHlQbRURERKqiNmE93bpwKrDF3e93973ADcCFbba/GLi+D0+l7wYZGroWOLdspbv/ibufFHv2VwBfcvcduU3OiuvXD7CNIiIiIsPPeks36DLl4FDgwdz81risxcPbkcDRwBdyi1eZ2WYzu9XMfnahT60fBpZy4O5fNrOjutx8aHv8IiIiIsNgASkHawu/dG9w9w35Q7bYpywP8iLgk+4+l1t2hLtvM7NjgC+Y2Z3u/t1eG9kPy55Da2arCZHcy3KLHfismTnw4cLJL+5/KXApwKGHHFKpfE7lzgool3apVf18i4j04LEOv3RvBQ7PzR8GbCvZ9iLgTfkF7r4tTu83sy8S8muXpUM7DP/R/x/g3wvpBqe7+wuA84A3mdlPle3s7hvcfb27rz/ggDWDbquIiIjIkksjhfU55eA2YJ2ZHW1mKwid1nnVCszsWGAN8NXcsjVmtjLeXwucDtxd3HepLHuElnDymtINcj3+R8zsU4Sk5S8vQ9tEREREhkK/69C6+6yZXQbcCEwA17j7XWZ2FbDZ3VPn9mLgBnfP/7T8fODDZlYnBEivzldHWGrL2qE1s/2BlwCvyS3bF6i5+9Px/jnAVcvURBEREZEhYAMZ+tbdNwGbCsuuLMy/o8V+XwFO6HuDFmhgHVozux44k5CQvBV4OzAF4O4fipv9HPBZd38mt+tBwKfMLLXvb939XwfVThEREZGhp6Fv2xpklYOLu9jmWkJ5r/yy+4GfGEyrhkPVLwYTaaVqF4eJiFTLYCK0o2IYcmhFREREpBNTwKCMOrQiMpZUvktEOhmm8pSmlIO21KEVERERqQClHJRTh3YJKXd2vlGJjg3Tt/jlplxaAb0Ploo+V8aIdV1bdiypQysiIiJSAYrQllOHVhZsVKKr/aBzISIig6YIbTl1aEVERESGnC4Ka08d2iWgHCcZR8qhFBHpJwOlHJRSh1ZERESkAkx1aEupQysiY031aEWkaCgr15guCmtHHVoRERGRoaeyXe2oqy8iIiIilaYIrYiIiMiwM3RRWBvq0A6Qqhss3LBeGa/XtHeqdiAi0h9KOSinDq2IiIjIkDMM08WrpdShlWVR1Whdp3YrgltdqnYgIkNZ3SAxQBHaUurQioiIiFSAynaVU4dWREREpAKUQ1tOHdo+00/ODf1MK1iun4F7/fmp7DnrfaGLw0REFsUMlBJVSh1aERERkQpQhLacOrTSN4uJvC0mArvYiF+76GmndnUbwW3VRkVtZRwoMj8Y+v8xppRDW0odWhEREZEhZ2aY6YthGXVo+2Qcvy0vJOLSbSR2YcfubR/z3iNHxde57Pl0E7ktPu64vIfyz1NROxGRHihCW0pnRkRERKQCrGY93bo6ptm5ZvYdM9tiZpe3WP9aM3vUzO6It9fn1l1iZvfF2yV9fKo9G1iE1syuAV4BPOLux7dYfybwT8D34qJ/cPer4rpzgT8DJoCPuPvVg2qndK/baFo3UdhOx+ol2uoL/F7m8SGMHioZdAiiZvmCbc5BWfR2XCO2w0YDLIiMn6EeUCEZQJUDM5sA3g+8HNgK3GZmG9397sKmf+fulxX2PQB4O7Ce8Ol4e9z3ib42skuD/I99LXBuh21ucfeT4i11ZtPJPQ84DrjYzI4bYDtFREREhl/Nert1diqwxd3vd/e9wA3AhV225qeBm9x9R+zE3kTnft/ADCxC6+5fNrOjFrBrdnIBzCyd3OK3haEwylG0pYjIlkViu4m69iv/0pkAunwtrfW3+LJ83FbHLJ6vcY7Y6gp4EZHuWf9/OToUeDA3vxU4rcV2v2BmPwXcC7zV3R8s2ffQfjewW8v9m9qLzOybZvYZM/vxuGyoTpCIiIjIsjMWEqFda2abc7dLWxy1qBg9+WfgKHc/EfgccF0P+y6Z5axy8HXgSHffaWbnA/8IrKPHExRfnEsBDj3kkEG0c+wsNjLbav9uI7Glkdwu2tRrlG9exYIuasWWRnNLIret3rlllRI65XCl9o1ipFZERDoxrPcqB4+5+/o267cCh+fmDwO25Tdw98dzs38B/HFu3zML+36x1wb2y7JFaN39KXffGe9vAqbMbC1dnNzCcTa4+3p3X3/AAWsG2mYRERGREXIbsM7MjjazFcBFwMb8BmZ2cG72AuCeeP9G4BwzW2Nma4Bz4rJlsWwRWjN7HvCwu7uZnUroXD8O/IB4coGHCCf31cvVznHSufJAdxHZdhUKyiKy845RFqn1dsdeXIS2aZ21z4mdHy0t5MWSrpRv0aZicLekMkKn3FpFakVExkyfB1Zw91kzu4zQEZ0ArnH3u8zsKmCzu28E3mJmFwCzwA7gtXHfHWb2+4ROMcBV7r6jrw3swSDLdl1PCEWvNbOthNIOUwDu/iHglcAbzGwW2A1c5O4OtDy5g2qniIiIyNAzBjKwQvyVfFNh2ZW5+1cAV5Tsew1wTd8btQCDrHJwcYf1fw78ecm6eSd32FQ9OtZLNLMYOey2YkE+GlsWgZ233Mu2K/8jXuwV8i1fS0/r6k3bZNFRax2pbWw3UVjeiLam89SpMkKn3Npucn6rYhirHagercjoq0T92Yz1PUI7SjT0rYiIiEgFLOCisLGhDu2Y6apawAJzZVMUtWWVg0JEtlMktniMesl+ZY/Xi1ZRzRSBrZVET9O3+rLIbVl1hKZ1hcoIxYhtr7m1rfYVGRbDGIWvIv1tjzGj7yOFjRJ1aEVERESGXtejf40ldWjHxCAis50qFuT3S5HVskhsWQS2XWS328oIZYqRjvx8ul+ndcQ1jdZSK6wvRm7LIrbEvZsPGvctidQm3dStVaRWRGS0GAMZKWxkqEMrIiIiMuzSSGHSkjq0PapaxGuhtWXz+3bKlS2tYJCLptZjDmlZRDZb7q2PnbarN0V9W7e3W/Ojro37NW8dYU3TWqozW4jE1gqR0ZrPpQebp7SWbXHUsaziQutILYxOzVrlWYqIlDHl0LahDq2IiIhIFahsVyl1aEVERESqQGW7SqlDO6KWMtWgmBbQSCtolKoqphjU07FK5ueKF5F5c0pC2KfQ7sI+ZYoXbCX5M1JLF2gVUg4mLF0sli4KK5bvat4vtSVLPWjz+MWWzBs+tyT1IGzT3SAMVUk9GCb5c6pBFkSqr1qDKeSYUg7aUYdWREREpAp0UVgpdWi7VJXI1nJEZhsXcjVvN9cUTW0dkU3b1AsR2GKkNi3Pf6/OIrLe3P56i8EX8mrFC73i5vmIaS1Gl1OkNkVm59LFX3F5LZ7PiewiscKgCKnFuSYVS3zViheB5VoRth+/SK0uDhMRaUER2lLq0IqIiIhUgS4KK6UO7ZjoZtCEbiOz80trxWgrtab19RYR2rKIbIrEztUnmubr2fL50dd6MUKbRX/bS63K8mPjIWv5CG28PxET8LP5bPlc03w9HqQsYps/z1nUN3sesT00R3+LLS6L1Ob3TboZfEFERCrETBeFtaEzIyIiIiKVpgjtiOg117AsXzas6y4ym6oYzHkxd3aiafuwrHldMSI7GyOws/XmCO7sXIrgxuPUcxHaejpGbHfcpusc2rjZRGxmrdaIck7UUkQ2zE9O1Jv2nYyJ+ZO11hHb1IQJmx+hTdKytE2Nudj+FBXuMlIbDhbXdZcbW5VcWhERyVHKQSl1aEVERESqQBeFlVKHtuJ6rWrQbvuFRmbrhdzZ2RSFredyaOO6mZKIbJrOzKXIbOvpXC5ome7X0zTLpW2eJo1qBinfNcynlKSJXDmUFLWdnIgR2RiqTfNTcYPJWj1OmyO2U3G+nvJlW5RaySKwqZ2pPT1GavPmVVkYkaoHIiJjTzm0balDKyIiIlIFSjkopQ5tB1WNXHWKzHqW21mbt02vkdlGXmycrzdHYwFmU4R2biJOW0dkZ2bDdO9smo/7x4G20jw0IrNzMcG2kUvrTeuT9MU2RWhTFHYiRl/zX3yn4l9GisxOTYbpislw7NnJ5shtithOxVzb9NCTrf75pMdNr0naZIGR2rBPHKHMm+u3dhupHVbDVo82nT+NGCZSPVX7/9eS/veUUodWREREZOiZIrRtqENbUWURq06Ro2JktqkObdlIYB0isyn/Nc038mQbbdmbIrOztTgfI7EzYTo9kyKzYfuZmbjfrMf5GBmdbUTMZ+L9eorQxmlW7aDeHF2vxTzWrLpBiswWorAAk/H+1FSKzLaeXzkVH3sqRYWb6+X6RIiy5s/zVIy8ZsUYssK4cVoSqW0fpUzh5+ZIbVFZpFa5tLJUhi3qXhX62xQM5dC2oQ6tiIiIyJBzWpfalEAd2orpZ73Z4vHm584WRvfqMjKb8mNTVBZgOkZmp2fiuiwyG9fvTVNvnk6HSOLMTPMUYHY25qvG5Nl6yh8t5NImKXc2BbFrqdpBqlgw2fjmOzVVa5quXBmnK6xpOq8aQ4zUrixWXGiciob0cPXCfCFSm0Y0SwHn+bm0rUYKax5NTPVpRUSqzpRD24Y6tCIiIiJVoA5tKXVoR0S39WbLas0CeMr7TJFYmiOzWZ3Zksjs3kIFgz0zuQhtjMju2ZumYfnu6RAJ3LOnOSI7PT3XNJ2NZQ5m8xHamZhbmkVmmyOZZVJJWMtyamMN2alGeyenUtQ2LFu5sjgN61etipHZlalObjyHKcrdYtSyea9NIVJrtSysGo4Zo6WpdVnkuWmgsOaE3GKN2lGpeiAiMs4GkXJgZucCf0b4mPmIu19dWP8/gdcDs8CjwP/r7v8V180Bd8ZNH3D3C/rewC4NrKtvZteY2SNm9u2S9b9sZt+Kt6+Y2U/k1n3fzO40szvMbPOg2igiIiJSCRZTDnq5dTykTQDvB84DjgMuNrPjCpt9A1jv7icCnwTelVu3291Pirdl68zCYCO01wJ/DnysZP33gJe4+xNmdh6wATgtt/4sd39sgO1ra5hyCBdzNXBZVYNsmosgltWhnTcCWEnObDEym6KxALvj/d17wvyuPSkyW4/TEG3dvTuUOdg7HaYpCpumc3OdI4rW4SrQ7Bhzzctn9jaK3E6kvNoYtd07HabTK8OfzD77hOnsbMwnjjnCc3OFKgfZ22h+Eu2891gqVBAXTxa2q1nh6vDc7vO/tXdX9aBM/j03DH8LujJeRIRBlO06Fdji7veHw9sNwIXA3WkDd785t/2twGv63Yh+GFiE1t2/DOxos/4r7v5EnL0VOGxQbRERERGpvFqtt1tnhwIP5ua3xmVlXgd8Jje/ysw2m9mtZvazvT+h/hmWHNriCXLgsxYu7/6wu28o29HMLgUuBTj0kEMG2shh1G3u7Lz9WuTQzqs3W5imEcBmC9OUO1uMzO6abhx7V4zM7o4R2V270jRERffsCeUO9u5pjswWI7L56GuqSpCiqMW811rhVKTc2mI+bhYFnm08VnrcNC1GitMxZmam4nbhT8kL56xVDm0xWpu+cKcoZBaRrTXPp1zaYn1a6L46gXJp+0MjholUx+j8f7OF5NCuLaRubij0qVodsOUHiZm9BlgPvCS3+Ah332ZmxwBfMLM73f27vTayH5a9Q2tmZxE6tC/OLT49nqDnAjeZ2X/GiO888YXZAHDiCccv/2+jIiIiIv1mLKTKwWPuvr7N+q3A4bn5w4Bt8x7a7GXA2wipotNpubtvi9P7zeyLwMnAsnRolzW8YGYnAh8BLnT3x9Py3Al6BPgUIcdDRERERPrnNmCdmR1tZiuAi4CN+Q3M7GTgw8AFsV+Wlq8xs5Xx/lrgdHK5t0tt2SK0ZnYE8A/Ar7j7vbnl+wI1d3863j8HuGqZmrmsluJiMM99p6lnQ90WB1KIpamyi8OaUw0agybEC79SysGeRnueiSkGu3eHn8mf2RlSDPbEFIPpmHKQftJP0sVZK1aFn/ZXr57K1u33rBVhul9ox36rw+PuszKsnyxchxWzBNgdv1vu3BUC+jt3hhU7n96bbbtrV0qBCNNiCkKapmF30+AO7qF99ewXrviTfu5noizFIP3cn83HFIN4AVct/kyWbVdISbDca5eVLuvzxWGgwRZERIZFv9Oc3H3WzC4DbiTkw13j7neZ2VXAZnffCPwJsB/wifhZlspzPR/4sJnVCR82V7v76HVozex64ExC/sZW4O3AFIC7fwi4EngO8IF4gmZjWPwg4FNx2STwt+7+r4Nqp4iIiMjws0FUOcDdNwGbCsuuzN1/Wcl+XwFO6HuDFmhgHVp3v7jD+tcTCvUWl98P/MT8PSSvXxeDtRxYIU7TxWGpTNdcVqYrTgvD2e4plOZKF4DB/MhsWQQ0XfS1MpbHetb+qwBYuzZEYw96TuN5H7wmRFQPWh2KZRzgjwKwek+Yn9r7TNNzn1mxb3jsVWsA2GEHAvDwrv0B2P7E6mzbhx8P7XnssfAYTz8ZntR0oZxYNphDyWgO6Rq2Wu71SkPupsjrRLrIy+LFbHG7uSxiG9fHwRJaDWdcfD11cZiIyOjRhajllv2iMBERERHpwgAitKNCHdoR02vubD0XoS3mzqbBARplusJ8GkBhbxxEYG82rG04Tho0IZXmgs6R2ZQru/pZIQH2wANDtPSIQ0LUct1zdwJwjGXp1jz7u7eFY99+OwAPfS2su/Pmh0P7Z5qjlLWp0M7DzjooTE/7UQCOPeUUAH7wwy/Mtr3/uWHdfY/sB8AD20Ju7KOP7grP4+nppvan51NkWTS2cZ4nYkh2ItYVm4zztVqK2IZjpqhuLc7XS3Jpw/1UJmxwubQiIrKM0khh0pI6tCIiIiJDzmkVrJBEHdqCYbiSexDDe5blzjYNfVvInZ039G09DXHbHJmdjsHJ3dPNw9mmQRPCsjikbUlk9ofWhIjsIYfsA8Dzjww5qifuFy6YXHtHyFe/79pPZ8fc/Nn/7v4E0IjYPhD3eyDbP5Q4PuKc52XbnvzaVwBw+EnnA/CtHwoR23v2CRHbbdtC5PipJ3Y1PZ/0/Gox+lqLzy9FZSE3IMREc4Q2i9jGSO1kHFghG9giVT8o5NJCPjLfWy5tL4ah2oGGwBWRsaYIbSl1aEVEREQqQF/my6lDWzG9VjcoW1+sOZu/X5Y7m0VoC3Vnp2Pu7J49KUI7F6eNvNLpDpHZI48I0xOP3g3AKXtuAWDnBz8KwOff9422z7MfHshFfB/47EcAeMGbQ37uGb/8awCsXncGAFOTob3/FbcvRmrT863F6OvUVH7I3rBsxWRzhHYynpOpiXAeZ+vd5dLm7xe/u3eKqqragYhIVZiqHLShDq2IiIhIFahDW0od2hFRVt0gW1+oblCsOZu/3yl3No24tTemyKbc2enpNBpYypdt5NCmuq2pzmyqZpByZlNk9oVPfRaA+//gPUBz1HQ5fD1Gho/4znYAXvh7bw0rjj4HgJnZ0P7ZeFJ2PhWqH6Tnm87B7tywZSlau3sqRW/DdGU8r+k8l+XSurfIoc1q0zZXO0ijjzW2VbUDEZFKMl0U1o46tCIiIiJDzpVy0JY6tCOutLpBi5HC5rIqB83TLKc2VTeYjfVoY4rszEyK0MaoZGE0rbw0AliqM5uqGaSc2X5EZs9413nhzkvOb17xpVAp4Zbf+UzPx2y0J7TvlKtCZHbXkWcDsHt3eD4z8Rzs2R0Si7NI7XQjWj09HaK1K1fGyHfMRd6bIrYxt7Y+2fr1SNNam5HCBlntQERElokitKXUoRURERGpAEVoy6lDO0T6WY6jU3WDVhHaYm3aRmQ25tLGCO1MDDZO723OnU0R2hSVTFf8Q6OqwbP2XwU0RgBLdWZTNYNeIrP7H7cvAP/y5n8H4JaN/wHAH90cN7i5uMfPhcl5YXrGBacC8DPvOx2AJ+9+puNjpvatPTa098Q3HAbAY4ccB8DTT4fnNxMTjNM5yEer03lKkdrpFTGHdkVzJDyd96mJmDPrnV+79Pr2Wu2gF4r6iogsB1PZrjbUoRURERGpAEVoy6lDWwGt3sBl39I6VTeYt32uykGxusFcPeVshvWzhSoHs3HFzEyIIM7sbY7Q5q1YNQXAAQesAGDdc3cCjRHAuq0ze8bX3pfdP/8doVICMTLbqxTRveXIkBe76aPTjXWnvbntvqn6wUtfFNq/7vgjANj+aMil3fl0eL67n2muegCN85TOWzECnqbpvKfXoV5rfn08Xyu28HZoVLNoXe2g+O4pq0cLqkkrIiLDTx1aERERkWFn6KKwNtShrbhi/dmO2xfqz9bzeZjZsiDl0GaR2rgi5dDOzjZHaIu5s6nmLMDq1SFiefCBYdkxdi8A91376a7anSKzWVS2jTf+zk8BcO7+/960/F+fDLmyH3jXl1vulz/2pvh4nSK1qf3H/N9Q7eC+A18AwGOPhue7Z3fzCGnQOE/ZeUvnMZ7XtGkWmS28Lu1euyxC3+F/nqsebV+k6LV+Buxd9muAcgLbUq76wo3er0vW9Wf9OFKHVkRERGTIORpYoR11aEVEREQqQL8GlVOHdkTNuzisZD7/S3OWalD8qbuQclDPUg/CzrOzYUH+Z3WAycnGH95+zwoXgx28Jgw48Ozv3gbA5g5lulJprnapBinF4Ki3rAfgibPCBWdfLR4rTt91wn4AfP+9m4HWKQjp8a6Oj19W0iuV8frR+HwOPvx4AL4Xn+9TT+4BYGZv49xkpbzieUvnsZ6lGhDnW78OcSTipteu+K297PXWz5ciItWlFJ1y6tCKiIiIDD0NfduOOrQjrttILcy/sKhTZLY+l6ap8H9z9G9yaiK7v99+4f5Bq58A4Jnbb++q/WnQhFaluVJkdv+zfgyAJ+Lysz74qtCuE05r2r5259cAuPkNn2ja7403/yfQOlKbHv/FbzipbTvT8zno2J8HYL/91gKNc5AGWoDGeUrnLZ3HRqS2eN4LAyoULwDLUWRWRGR0KYe2nDq0IiIiIkPOUcpBO+rQjqlu/igaZaLCNEUMs/kYWUxTrxcjtLkc2tXh8dbwOAAPfe3ertp5S5tBE7Kc2Tj/kn/7PwDce8CLAJj1iabtJ1/yorhdiNx+6cW/1XQcDnt36eO/uEM70/NZ8+rw/PZbfWB4zKkWg2LUm89bdv5KznM3hWf0T05EZMSZUg7aUYdWREREpAIUvCg30A6tmV0DvAJ4xIsBgZcAACAASURBVN2Pb7HegD8Dzgd2Aa9196/HdZcAvxc3/QN3v26QbR11qQh/uz+GLHfWm7epZ4MxxO3cW84XArRMTjYipPvEIgX77g4RzG/f/HDb9p7xrvMA+KObm5envFloVDNIObMpMjthYeCCZ0083bTvrnoYlnZL3C7tl3Jq3/jexrGL+bSpPbf8zmdatndrfD6HxueXnm/+HCT1DucxGzih7HXo5jWMU32XFxEZHYOI0JrZuYS+2ATwEXe/urB+JfAx4BTgceCX3P37cd0VwOuAOeAt7n5j3xvYpdIzY2abzOyoRR7/WuDcNuvPA9bF26XAB+NjHwC8HTgNOBV4u5mtWWRbRERERCrLsZ5unZjZBPB+Qn/sOOBiMzuusNnrgCfc/UeA9wB/HPc9DrgI+HFCX+8D8XjLol2E9lrgs2Z2HfAud5/p9eDu/uUOneILgY95CE3dambPNrODgTOBm9x9B4CZ3UQ4Wdf32gZpz71NtK/kwvh6dvV96w3SkLe13KFToHJyZlfYd6bDVfcvOT9MCxHa/HC2qc5sqmaQcmZTZHaqPt207+r49e2J2Wc37QefmHfsD5S1h9YR2vR80vNLzzedg/wwwPP2rTfXoS0qex3avXYiIjJafDBlu04Ftrj7/QBmdgOhb3Z3bpsLgXfE+58E/jz+wn4hcIO7TwPfM7Mt8XjFMvBdM7NVwBsJl6448G/AB919T6d9Szu07v5xM/sX4Epgs5n9FbnrU9x9/hU0vTsUeDA3vzUuK1suIiIiMpYGkEPbqr91Wtk27j5rZk8Cz4nLby3su9i+2seAp4H3xfmLgb8CXtVpx045tDPAM8BK4Fl0d8F1L1q9Mt5m+fwDmF1KSFfg0EMO6V/LRERERIbIAurQrjWzzbn5De6+ITffTX9r0X21Hhzr7j+Rm7/ZzL7ZzY6lHdqYJPxuYCPwAnfftbg2trQVODw3fxiwLS4/s7D8i60OEF+YDQAnnnC8qseLiIjISFpAqtlj7r6+zfqyflirbbaa2SRhJPkdXe7bq2+Y2U+6+60AZnYa8O8d9gHaR2jfBrzK3e9aZOPa2QhcFnM2TgOedPftZnYj8Ie5C8HOAa4YYDvGlln5d4CyL4IpHbRWa72Bx2TQfIrt7FycToVKA7WpeDV+WS7tlzbFOz/XtPhfnzw9u79/ak8cASzVmU3VDFYXUo3S8slYBSHt1+rYUBg1LGtPa+n5pOeXnm9W0SCXIFvMp03nsSzNtux1aPfaiYiIdOE2YJ2ZHQ08RLjI69WFbTYClxByY18JfMHd3cw2An9rZu8GDiFc4F9ePL47pwG/amYPxPkjgHvM7E7A3f3Esh3b5dCeschGYWbXEyKta81sK6FywVQ8/oeATYSSXVsIZbt+La7bYWa/TzjRAFelC8RERERExo/hfS7GGHNiLwNuJJTtusbd7zKzq4DN7r4R+Evgr+JFXzsInV7idh8nXEA2C7zJ3ecW2aR2lbHaGmgdWne/uMN6B95Usu4a4JpBtGsc1WJaS71NeovFdbVC5C/NWwwVNqY0zaeAbXo3z8423te7Y8GBZ/Z5DgCHnXUQAA989r9btiWr93pec4Q2Xx/2XSfsBzTqyKYRwFKd2VTNIEmR2R/Z8dWm/dbE4/xRofZsy/aUSM8nPb/0fPPnIKkVzlvxPKb50tehm9dw0WlMIiIyTAY19K27byIEGPPLrszd30PJRVnu/k7gnX1sy38tdF+NFCYiIiJSARoprJw6tGOqXZQvST9spIjhRK0wP2FNUyuEaGdnGnmjO3eFx3uCGKE97UeB8ghtcsYFpwJwy8b5aTnff2+4cHP/s34MgC+9+LeAxghgjTqz8fnEnNkUmS0ehxYR2vT4JeVnM4fG57M1Pr/0fPPnIEnnad75KznP3fzA1M3rKSIi1aYObTl1aEVERESGXnejf40rdWhHXDFy127eCjmatVqYTmRVDcJ0ajLmzE6kaS3u3/yHNjvTyB/duTPcf3hXqE1w7CmnxDXleasAP/O+UHngliPfM29dyqd9483/CcBRbwmVSRoR2E/M2wcaObMpMvuBNrmz6fGfbNtK2Dc+n/T80vPNn4MkyzmO5y2dx+y8xvPcOO/Nr0vxdWo6dofXW0REqksjRJZTh1ZERERkyA3qorBRoQ7tiEqRufTmL4vc5YOqKVczRQTTVfWNSG1zBDFFFCcnw4KJGFKciceZnc3l0D69F4DtT4Q6rT/44RcCcMQ5zwPKc2mfvPsZADZ9NJQNOP8dK+dtk0VYDwujMb/xvT8FwLn7N9diTnVms2oGbSKzm94RHu+W054p3Sbf/vR8tj+8AoCdT4dxSPLnIEnnKZ23TpHZ4utQzG2G3iLxIiJSTerQllOHVkRERKQC1KEtpw6tiIiIyNAz5dC2oQ5txZnHn5It/LTtTLTfvnjhl8+/KCxLPbCUatB8cdhUfNdMxp/Kp6biT+hT4bEn9ob5ubnGz+27doVEhO2PhmX3PzeUuTr5ta8A4IHPfqRtu2857c0AbPra+7JlrdIPoJGC8IF5a9pfgJbSDPKP18m62P5veHg+6fml55uGvE1pBtA4T9l5S+cxnteJLPWgkHIQ92+8Ti0u6OtyOFwjtMtcqQiL4dbfUXvGiSJN3SlLG5PO0t+n+fzUrypyoK6/m1Lq0IqIiIhUgL4IllOHtgLy3y6zb5yFi76ybYsXg6XIXXaM5ohSPqJXi1HeNJ2IVynFqlJMTnicWtM0RRqnVoTI48zeMM1HaPfuCRHLHTvCxWH3PRJKZx1+0vkAvODNtwPw9fd9o+U5SPKR06uP2xeAf3lzuPir1eAL7aRBE7LSYB0uAMt7wZtPBuCx2P77toXns2PHHqDxfJMUlYXGecoitNn5jNvG85zOeyNS2/z6tIvGZhFY6+7isHkXkY1IRENEZGS4yna1ow6tiIiISAUoQltOHdohUhZ17eexUt5lneZc2vz9LEfTUkS2HqfNAwCsXBGnK2txGkKMe6fDND+oQIrWPv1kiGA+sG0KgG/9UMg9PeOXfw2AI76zPazvMCQuNEp6vfgNJ4VpXH7Gu84Ld15yfvMOX9oEwC2/E8exjZNOgybkpTJd+8X23rIzDuG7LTzX9PzS881KdOUitOk8ZedtRXMObSMSniKyrQdUaPXa1bqMwC6GcvlERJaDLgprRx1aERERkSGngRXaU4d2xBWjeWVTgIlU1aAwbURqw3TFZJhOTVnTNEUep1eGt1WrCO309CwAjz4aBh64Z5+Qe7p63RkAvPD33hr3CEPddhOpLcoisCkE2wcpMntMbN9tq0J777kvPOdHH90JNJ5fkiKzK1Y2/tTSeWqcv7A8ndd0nmslr8dEi6Fvu3l9RUSk2hShLaeaMyIiIiJSaYrQjohiPdr0XaWs2oEVqiXk7xevpp+MQ95OpaFtYzroivju2WdlWD8zE9bvs0+M0M7Oj9CmqO2up0PN123bYrRyMgyJy9HnAHDKVfsAsPbYjwKdqx8MSqpmkHJmU2T2W98L7du2LUSa0/NJdWezyOyqcC7SOYFG7mw6byuy3NkwnUo5tLXm16FY3SBfwaCRX9u+uoHqz4qIVJfqz5RTh1ZERESkApRyUE4d2opJEdZO9Wiz7TtUO/Dc9720TS0eO0UI6948n1U5iHVUZ+fC/KpVYTo7G+vRzkxlx67Pxcerh2mK2D71RIhw/lfcbmY2RD53HXk2ACe+4TAAXvqiUKHgvms/nR1zIfm17aQ8WWiMAJbqzKZqBilnNkVmU/uLVQ1WrgrPfVU2bVQ5SOdp5Yowv3Iq5iRPpoh4mhYjtamSQb1pGu4vrLqB6s+KiFSDY7oorA11aEVEREQqQBHacurQFgzDuNn9rEdbPOa8qc3PoZ2Ikb+6pYht61zauRhZTBHa2ZUpQhvXzzXeXvUYwazHCG0aSasYqU15t7t3h5zaxw45DoB1xx8BwDH/9+zsmD/63dsAeOb2MMrYQ1+7F4CtNz8cHmum+TWsxaoCh511EACHnhairvuecgoAP/jhF2bbfsPDujQCWKozm6oZpJzZYmR2RRaRDc999erJON+4/jLlzq5M1Q1ShLZD7uxEuxzaJahuMAwVExSdEJFxpv+B5dShFRERERl2DvXljysMLXVoR0xZtYNsfaHaQS33bS/l03bKpZ2aSNHWGJGNEca5OD8XI7bujcd2b+TT5hUjtTufCpHPmek48tbTqwB4+PEQsb3vOS/I9j348OMBOOjYnwfggIsfBeDQPU+Edu59pumxZlbsC8CuVWsA2GoHhmPv2h+A7Q+vyLZ9+PHQnsceCyN/pRHAUp3ZVM2gGJldvTpM990vzYf1q1c1zvOq+DApMrtiIkVoe8udrTVVqChEb+dFbFXdQESkyjSwQnvq0IqIiIhUgHJoy6lDW1G9VjvI9muTa1mM7qVc2vT3M5FyaOP29bhiZTxEisim5fVchLZecvF8LR5zOkZqU53aPbv3AjCzN0REdz4dHvWRhxuR3u89K4Q699tvbZiuDhHXfVaG9ZONwgLh2CH9ld0hCMzOXaHhO3fOxcfYlW27a1doTzGCnKQ6sysLObMpMrvPPmH9PjF3dp9VjX1XrQiPu3IqVY2I04nQjolanFqats6dbTdSWCeqbiAiUj36ka3cQEcKM7Nzzew7ZrbFzC5vsf49ZnZHvN1rZj/IrZvLrds4yHaKiIiIDDej3uNtUY9mdoCZ3WRm98XpmhbbnGRmXzWzu8zsW2b2S7l115rZ93J9uZMW1aAOBhahNbMJ4P3Ay4GtwG1mttHd707buPtbc9u/GTg5d4jd7j7QJz/sFlPtoNdc2rBFzH0t5NJOxGoHk5YitnG7GFH0LEJrLefzj1+rpcePeboxB7U2Eeb37gkR2RSpTZHR3c+EsOqe3TPZEZ+Kea0pWjoZ6+JOxtBsrXDaUjJ9qqQwO9M8etnsbK4mbyGknHJliyOArcpyZ1NVgxSRTbmzYf99VjS+Vqe6sytjZHZFPI/p/E5mVQ3a585arg7tIHNnh6G6gYjIuHOWPOXgcuDz7n51DEpeDvxuYZtdwK+6+31mdghwu5nd6O4pQPnb7v7JpWjsICO0pwJb3P1+d98L3ABc2Gb7i4HrB9geERERkcpy7+22SBcC18X71wE/O789fq+73xfvbwMeAQ5c9CMvwCBzaA8FHszNbwVOa7WhmR0JHA18Ibd4lZltBmaBq939H0v2vRS4FODQQw7pQ7OrpV+5tAA15tLKpmNkkdqatVzuhVxVmLegEZFN7YzzEzEyOxWjq7tjdHXvdHOkthixBZjZW4/T2ebHqrX/nlaMvrZSjMhmkdmV4U9mn31ShHYiTpurGexTiMyuykVoV02F55JFZrORwWJktkMd2vQ6tcqh7US5s/2R/t5EZPhln48j8P9uiascHOTu2wHcfbuZPbfdxmZ2KrAC+G5u8TvN7Erg88Dl7j49qMYOskPb6qyXfepeBHzS3edyy45w921mdgzwBTO7092/W9zR3TcAGwBOPOF4/TYqIiIio2dhdWjXxuBgsiH2mwAws88Bz5u/G2/r5UHM7GDgr4BL3LNvDlcA/03o5G4gpCtc1ctxezHIDu1W4PDc/GHAtpJtLwLelF8QQ9e4+/1m9kVCfu28Dq2IiIjIqFtgDu1j7r6+9JjuLytbZ2YPm9nBMTp7MCGdoNV2PwT8C/B77n5r7tjb491pM/so8Fu9Nr4Xg+zQ3gasM7OjgYcIndZXFzcys2OBNcBXc8vWALvcfdrM1gKnA+8aYFvnGYYhcJN8G3r9uaHbi8PiwZseL/3UnR4zvVmyv6fCr67z29ZIPYgZBtQKqQaTk+lCrubUg+npiThNF2w1X8gV7qeL0mKaRfzq2ukbbEonsJhCkdIfUlpBuN98gdnKlcVpWL8qpRjE4WzToAmN0lwx5WCq8eNDGkAhmxbKdE3GaTHloJthi/t1MdgwvO/zVExcRGTJy3ZtBC4Bro7TfypuYGYrgE8BH3P3TxTWpc6wEfJvvz3Ixg4sEczdZ4HLgBuBe4CPu/tdZnaVmV2Q2/Ri4Ab3ppfp+cBmM/smcDMhh/ZuRERERGQpXA283MzuI1SsuhrAzNab2UfiNr8I/BTw2hbluf7GzO4E7gTWAn8wyMYOdGAFd98EbCosu7Iw/44W+30FOGGQbRs13V4c1ij11LxdnGk6VopXZseIk8lUD6t10LdldM+yi8HCNF1gNhkjtStihHZ6RZyubI7UzsTI7Ew+QhvLbNXn0tC8KVIb2134Kpsisel6nlqhdFiKFkMjUpymKSK7MravMQ3br4zjPaThbNOgCcXSXNAmMlu8GIz0OqQIbevhbfP386W8mp67LgYTEam8xdaW7YW7Pw68tMXyzcDr4/2/Bv66ZP+zB9rAAo0UJiIiIlIBGimsnDq0FdPrYAtZZC5NrBixnf/XkSKEKRc1DayQbVrMoS2L1Nr8+xMpyjthLadTU80R0L0xR3VmJhx8drbR3pl4vz4XpnNxmv7g64Vk2lqWMxvbMpEitPGxJxsNbuT0NkeQp1IkNv7lpBzZFJldMeHxWIVBE2qNiGhZZDbNZ5FZa47MluXStpJyZ3vNhR223FkREQkcW+qBFSpFHVoRERGRYbewsl1jQx3aiirNkS3k0s7br0PVg1aPURxwoSxSmy66b0QQ81UOUuWEOCxtrTmXNkVHZ2ZjZHaqeX52Ls032pfGR0iR2blC7mxx/ITisLsxdbYRqc2diqn4lxGLHGTzKyZTBNbj+jifTVMFg+bhbFNUNjxucwS2LDJbHEChbBrud1fVoCx3VpFZWSqqWLEww1R5R5aPUg7KqUMrIiIiUgH6QlhOHdoOqvqtuFPVg/mRWihGaxu5tLW4tn2kNjtHLYK+E+5Nq1IUcjKGRWdTpYFCPmuKzKZpbuTb7H69brGdsVkxx6j4TTblzjZq4sZpitTm2p3upwhscZoissVhaydrcX2tuZbsRC6Hdl6d2Q6R2dI6tLmKBsXI7KhUNRi2f94a8lakuqo+BK6jlIN21KEVERERqQClHJRTh7biOlU96LY+bVjXOq+220ht06hjQK3WmJ/z5hqqWR3aekhSna03R2zrWYS2FvePx6k32p0is43c2bi8w1WgtVTHNVU7SHVpc+2diPcbVRnqTfsWI7ITWeWCtF/zfC0XCU/r0muxmMhsUa+R2ar98iAiMs7UoS2nDq2IiIjIkHPvHLAZZ+rQjojF1qeFRo3aXiO1Keo6V9i+6Sr8LHcpHiv+UU6URGzn4vpUNWAuy5PNRWgLubIpd7ZTdlRqpRUitbVchLmWtcub5611RDbLh7Vi1LW5hmz+/rwRwHqMzObzZlVvVkRk9ClCW04dWhEREZEKUIe2nDq0Y6KsPm1TpK4wmlinSG2KBmcjimXrU8R2fgQxRSzrlnJj4/JCxDZFW1Oktp5ya3PtKFYzaLSnyxxamiO0+RzgWrZtqkDgTdt0ishm2xUqGDQ9bocRwDpFZtuPFFbNq3hFRKScqhyUU4dWREREZMg5aOjbNtSh7VJV6tH2WvWg1b7dRmqTLGLrhcfONWFeHdySiK3H9qXI7AQpCpuir402FGOQWcS2Qx5xIzJbqMrQ4jlZIZo7L1JbEpFtzBdGW2vx+EsZmR32928ybPVnRUSWnSvloB11aEVEREQqQCkH5dShHVHLEalNEcd5o5EBFjNsG8dq3jaL1MZj16w52tqoYNCIQBa/qfYa1ZtXszW3e61DFHVe3mshImuFyGw+GjzvcRWZHWoaHUxktOT/pnW9wehQh1ZERERkyIUc2uVuxfBSh1ZERESkAtShLacObY+qcnFY0m3qAbQp6VV8qtacelD6GPmf8H2uaRuL+6YUhHrJec1SDqw59SC/rLhtt8p++ofyVIPi+mKKQXGQhOJx2h4z7evl7Wpa3sVPZVV5nya6GExEpJxyaMupQysiIiIy7FTloC11aMdEN0Pjdhp8Yd7FXlaMEJYP2uDzLvIqlsxKkdhay3amCG6r9i82qtcqilmMwBa3LUZii/uVRWFbHSNbPsaRWRERac+Buq5hK6UOrYiIiEgFKEJbTh3aMdNLpDZJUdPSSG22Y4paxv1alPfKjlGI2GbLC8cs5ty2avYgIrSNde0jscVjlJfkmv+1et75K2lPL2VlFJmVYaX86P6o2nUc0l/q0JZTh1ZERERkyLnrorB2Blox3MzONbPvmNkWM7u8xfrXmtmjZnZHvL0+t+4SM7sv3i4ZZDsXwrFKRxwsewaNW+m2Xm+KEpbtZ+6FwRTq2a3GHDXmsn1q1Jtu2XKrU7M6Zo5ZY37CZpmw2ew4+dsEs0zQel27W7v9sscraU/H55E933rLXNnm81R2Putto7O9vIbDbhj/ntxqGlRBZMRV7e/c3Xu6jZOBRWjNbAJ4P/ByYCtwm5ltdPe7C5v+nbtfVtj3AODtwHpCHvTtcd8nBtVeERERkWG2lH3U2Bf7O+Ao4PvAL7bqh5nZHHBnnH3A3S+Iy48GbgAOAL4O/Iq77x1Uewf5teRUYIu73x+fwA3AhV3u+9PATe6+I568m4BzB9ROibqN1HaK2KYIZP7W2LZeuDVHOruN4A761jkC2zxffF6Nc9bifHSIyHYbmRURkfFSr/d2W6TLgc+7+zrg83G+ld3uflK8XZBb/sfAe+L+TwCvW3SL2hhkh/ZQ4MHc/Na4rOgXzOxbZvZJMzu8x31FRERERp5777dFuhC4Lt6/DvjZbnc0MwPOBj65kP0XYpAXhbVKiCue3n8Grnf3aTP7dcITPrvLfcODmF0KXApw6CGHLLy1kummEgJ0rlvbpOwPq1DLtqxCQqtjFtu32KhlN/u3qlYA3VcsaL1vd1+jFZUVERlvS3xR2EHuvh3A3beb2XNLtltlZpuBWeBqd/9H4DnAD9x9Nm4z8MDkIDu0W4HDc/OHAdvyG7j747nZvyCEp9O+Zxb2/WKrB3H3DcAGgBNPOF6f+CIiIjKSFhB1XRs7m8mG2G8CwMw+BzyvxX5v6+ExjnD3bWZ2DPAFM7sTeKrFdgPtow2yQ3sbsC4mBT8EXAS8Or+BmR2cev/ABcA98f6NwB+a2Zo4fw5wxQDbumCjXBOw+JzKIrZldWu7OVZZZHP+KGSt2peOtbjMmbKoa8ttO/w3UUS2N8NW2UBEZJh57yHax9x9fenx3F9Wts7MHk79NDM7GHik5Bjb4vR+M/sicDLw98CzzWwyRmnnBTX7bWA5tPEJXEbonN4DfNzd7zKzq8wsJQ2/xczuMrNvAm8BXhv33QH8PqFTfBtwVVwmIiIiMnZSHdpebou0EUhlUy8B/qm4gZmtMbOV8f5a4HTgbg81w24GXtlu/34a6MAK7r4J2FRYdmXu/hWURF7d/RrgmkG2T3qz0Iht0z6FUcfKN+ylXXPx2L1F+zpFW1s/Vnf79DK6V6/HlsGqUk1KEemP7LNpAf+7R9jVwMfN7HXAA8CrAMxsPfDr7v564PnAh82sTgiSXp0rz/q7wA1m9gfAN4C/HGRjNVKYiIiISAUsZR3aeJ3TS1ss3wy8Pt7/CnBCyf73E0q4Lgl1aPtklHNpy3QbsW3ap+Tbb1eVEkp0zMft1KZFvGaL+TY/Tu+VPOXNiogsTF1j35ZSh1ZERERkyDlLG6GtGnVopW8WErHN9u0Q6WyX1zjISGe/8qnGNRorooj8YIzjr4Jjrz+DJYwsdWhFREREhp5TV4+2lDq0fZaPRoz7N+d2z7/XqE1Vrjwd99e8HUXqREQWpyIfhctCHVoRERGRIRdyaBU0KaMOrYiIiMiwc6grQltKHVpZFp1+mh/Wn6eVUjC6NKCCiAz7AAuK0JZTh1ZERERkyDl9Gc52ZKlDO0Aqq7JwOmejY1ij7SIileLg6tGWUodWREREpAKUcVBOHVoRERGRCtDQt+XUoRUREREZcu6ui8LaUIdWRMaaqhuISNGwVjsYsuYMFXVoRURERCpAQ9+WU4d2CajagYwjVTcQEZGlog6tiIiISAUoh7acOrSyYP3ILRqV/EWdCxERGSR3VTloRx1aERERkQpQgLacOrRLSLm08w3bFaSyeMqdFdD7YKnoc2W8aKSwcurQioiIiAw5d1eVgzbUoRWRsaScZRHpZNj+TyhCW04dWhEREZEKUIe2nDq0y0A5TzKKlDMpIjJADurPlhtoLN3MzjWz75jZFjO7vMX6/2lmd5vZt8zs82Z2ZG7dnJndEW8bB9lOERERkWHmhAhtL7dxMrAOrZlNAO8HzgOOAy42s+MKm30DWO/uJwKfBN6VW7fb3U+KtwsG1U4RERGR4ee493ZbDDM7wMxuMrP74nRNi23OygUf7zCzPWb2s3HdtWb2vdy6kxbVoA4GGaE9Fdji7ve7+17gBuDC/AbufrO774qztwKHDbA9IiIiItUUB1bo5bZIlwOfd/d1wOfjfHOTQj/uJHc/CTgb2AV8NrfJb+eCk3cstkHtDLJDeyjwYG5+a1xW5nXAZ3Lzq8xss5ndmnr7rZjZpXG7zY/veGLorkhsx7FK5x261Sp1vodV1c9j1d7HVT/fIjK+ljJCSwhCXhfvXweU9sWiVwKfyQUql9Qg/6u3+oRreXbN7DXAeuBPcouPcPf1wKuBPzWzH261r7tvcPf17r7+OQfMi4aLiIiIVN4y5NAe5O7bAeL0uR22vwi4vrDsnfE6qfeY2crFNqidQVY52Aocnps/DNhW3MjMXga8DXiJu0+n5e6+LU7vN7MvAicD3x1ge0VERESGky+obNdaM9ucm9/g7hvSjJl9Dnhei/3e1suDmNnBwAnAjbnFVwD/DawANgC/C1zVy3F7McgO7W3AOjM7GniI0HN/dX4DMzsZ+DBwrrs/klu+Btjl7tNmthY4neYLxkRERETGyIJGCnss/trd+ojuLytbZ2YPm9nB7r49dlgfKdsW+EXgU+4+kzv29nh32sw+JbOG6AAACkVJREFUCvxWj23vycBSDtx9FriM0Fu/B/i4u99lZleZWapa8CfAfsAnCuW5ng9sNrNvAjcDV7v73YNqq4iIiMiwW+KUg43AJfH+JcA/tdn2YgrpBrETjJkZIf/224ttUDsDHVjB3TcBmwrLrszdb/nNwN2/Qghdj4WqD7SQLrAxry9zS6ql6hcmVelCMBGRqnPox4Vevbga+LiZvQ54AHgVgJmtB37d3V8f548ipJh+qbD/35jZgYRrqu4Afn2QjdVIYSIiIiLSxN0fB17aYvlm4PW5+e/TooqVu589yPYVqUMrfaNIbXeqHpkVEZFlEOvQSmvq0IqIiIhUwLgNZ9sLdWiHSNVzaWU8KHdWRGQ59GWwhJGlDq2IiIjIkHMHryulr8xIdmiVy7m8dP7nU97s8tNrICJVpxzaciPZoRUREREZNUo5KKcOrQxMPiI2rtFaRQVFRKQvvC+DJYwsdWhFREREhpyjKgftqEM7hEax2sG45dWOYmRW1Q1ERJZXfUw+QxdCHVoRERGRYeeK0LajDq2IiPSNIvnLYxR/2ZNmjnJo21GHVkRERKQCVOWgnDq0sqRGPZd2FHNnRURkCDjUNbBCKXVoRURERCpAKQfl1KGVZTFqkVpFZkVEZJAcx0fkM3MQ1KEVERERGXaqctCWwkoiIiIiUmkjHaGt+s/aKsMiw6DqZZiUDiIio0IR2nIj3aEVERERGQ2ukcLaUIdWllXlo+iK/omIyBJw5dC2pQ6tiIiISAW46tCWUodWREREZNgpQtuWOrQiIiIiQ091aNsZaAKgmZ1rZt8xsy1mdnmL9SvN7O/i+q+Z2VG5dVfE5d8xs58eZDtFREREhpkD9br3dBsnA4vQmtkE8H7g5cBW4DYz2+jud+c2ex3whLv/iJldBPwx8EtmdhxwEfDjwCHA58zsR919blDtFRERERlarhzadgYZoT0V2OLu97v7XuAG4MLCNhcC18X7nwReamYWl9/g7tPu/j1gSzyejCi3WqUqBlStvSIiUnWO13u7LYaZvcrM7jKzupmtb7Ndy1/jzezo+Ov7ffHX+BWLalAHg/xEPhR4MDe/NS5ruY27zwJPAs/pcl8RERGRseFe7+m2SN8Gfh74ctkGuV/jzwOOAy6Ov7JD+NX9Pe6+DniC8Kv8wAzyorBWwwsVvy6UbdPNvuEAZpcCl8bZ6aPXHfvtrls42tYCjy13I4aAzkODzkWDzkWDzkWDzkWDzkXDscvdAGDJqxy4+z0A4YfzUtmv8XHbG4ALzewe4Gzg1XG764B3AB8cVHsH2aHdChyemz8M2FayzVYzmwT2B3Z0uS8A7r4B2ABgZpvdvTQsPk50LgKdhwadiwadiwadiwadiwadiwYz27zcbQBwfBhzaFv9on4a4df2H8Rf39Pygf7SPsgO7W3AOjM7GniIcJHXqwvbbAQuAb4KvBL4gru7mW0E/tbM3k24KGwd8B8DbKuIiIjI0HrmyXtv/PdPn7m2x91WFTrkG2IgEAAz+xzwvBb7vc3d/6mL4y/6l/Z+GViH1t1nzewy4EZgArjG3e8ys6uAze6+EfhL4K/MbAshMntR3PcuM/s4cDcwC7xJFQ5ERERkXLn7uQM45ssWeYiyX9QfA55tZpMxSlv6S3u/DHRgBXffBGwqLLsyd38P8KqSfd8JvLPHh9zQeZOxoXMR6Dw06Fw06Fw06Fw06Fw06Fw06FyUa/lrfPy1/WbCr+83EH6N7ybiu2DmPl6Fd0VERESkPTP7OeB9wIHAD4A73P2nzewQ4CPufn7c7nzgT2n8Gv/OuPwYQmf2AOAbwGvcfXpg7VWHVkRERESqrHKV4RcznO6o6eJcvNbMHjWzO+Lt9cvRzqVgZteY2SNm1rJsmwXvjefqW2b2gqVu41Lp4lycaWZP5t4XV7barurM7HAzu9nM7onFwX+jxTZj8b7o8lyMy/tilZn9h5l9M56L/6/FNmPxOdLluRibzxEIdVXN7Btm9ukW68bifVFVA82h7TdbxHC6S9/aweryXAD8nbtftuQNXHrXAn8OfKxk/XmEahnrCCVFPhino+ha2p8LgFvc/RVL05xlMwv8L3f/upk9C7jdzG4q/I2My/uim3MB4/G+mAbOdvedZjYF/JuZfcbdb81tMxafI3R3LmB8PkcAfgO4B/ihFuvG5X1RSVWL0C5mON1R0825GBvu/mVCpYwyFwIf8+BWwtWXBy9N65ZWF+diLLj7dnf/erz/NOFDqlgHcSzeF12ei7EQX+udcXYq3oq5d2PxOdLluRgbZnYY8DPAR0o2GYv3RVVVrUO7mOF0R023wwP/Qvwp9ZNmdniL9eNCwyk3e1H8mfEzZvbjy92YQYs/DZ4MfK2wauzeF23OBYzJ+yL+rHwH8Ahwk7uXvi9G/HOkm3MB4/M58qfA7wBloxeMzfuiiqrWoV3McLqjppvn+c/AUe5+IvA5Gt8sx9G4vC+68XXgSHf/CcIVrP+4zO0ZKDPbD/h74Dfd/ani6ha7jOz7osO5GJv3hbvPuftJhNqYp5rZ8YVNxuZ90cW5GIvPETN7BfCIu9/ebrMWy0byfVFFVevQ9jKcLtY8nO6o6Xgu3P3xXImMvwBOWaK2DaOuh1Mede7+VPqZMdaKnjKzXkefqYSYF/j3wN+4+z+02GRs3hedzsU4vS8Sd/8B8EWgWLB+XD5HMmXnYow+R04HLjCz7xNS+M42s78ubDN274sqqVqHNivga2YrCAV8Nxa2ScPpQm443SVs41LpeC4KuYAXEPLmxtVG4FfjVe0/CTzp7tuXu1HLwcyel/K+zOxUwv+Bx5e3Vf0Xn+NfAve4+7tLNhuL90U352KM3hcHmtmz4/19gJcB/1nYbCw+R7o5F+PyOeLuV7j7Ye5+FOHz9Avu/prCZmPxvqiqSlU5WMxwuqOmy3PxFjO7gHCF8w7gtcvW4AEzs+uBM4G1ZrYVeDvhAgfc/UOEEevOB7YAu4BfW56WDl4X5+KVwBvMbBbYDVw0ov+UTwd+Bbgz5ggC/G/gCBi790U352Jc3hcHA9fFSjE14OPu/ulx/Byhu3MxNp8jrYzp+6KSNLCCiIiIiFRa1VIORERERESaqEMrIiIiIpWmDq2IiIiIVJo6tCIiIiJSaerQioiIiEilqUMrImPBzA43s++Z2QFxfk2cP3K52yYiIoujDq2IjAV3fxD4IHB1XHQ1sMHd/2v5WiUiIv2gOrQiMjbi8K+3A9cA/wM42d33Lm+rRERksSo1UpiIyGK4+4yZ/Tbwr8A56syKiIwGpRyIyLg5D9gOHL/cDRERkf5Qh1ZExoaZnQS8HPhJ4K1mdvAyN0lERPpAHVoRGQtmZoSLwn7T3R8A/gT4P8vbKhER6Qd1aEVkXPwP4AF3vynOfwD4MTN7yTK2SURE+kBVDkRERESk0hShFREREZFKU4dWRERERCpNHVoRERERqTR1aEVERESk0tShFREREZFKU4dWRERERCpNHVoRERERqTR1aEVERESk0v5/0h4L5qEufrEAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "wave.animate(variable, save=True)" ] } ], "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.7.6" } }, "nbformat": 4, "nbformat_minor": 2 }