{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<img src=\"../_static/pymt-logo-header-text.png\">"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Coastline Evolution Model + Waves\n",
    "\n",
    "* Link to this notebook: https://github.com/csdms/pymt/blob/master/notebooks/cem_and_waves.ipynb\n",
    "* Install command: `$ conda install notebook pymt_cem`\n",
    "\n",
    "This example explores how to use a BMI implementation to couple the Waves component with the Coastline Evolution Model component.\n",
    "\n",
    "### Links\n",
    "* [CEM source code](https://github.com/csdms/cem-old): Look at the files that have *deltas* in their name.\n",
    "* [CEM description on CSDMS](http://csdms.colorado.edu/wiki/Model_help:CEM): Detailed information on the CEM model."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Interacting with the Coastline Evolution Model BMI using Python"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Some magic that allows us to view images within the notebook."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "%matplotlib inline\n",
    "import numpy as np"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Import the `Cem` class, and instantiate it. In Python, a model with a BMI will have no arguments for its constructor. Note that although the class has been instantiated, it's not yet ready to be run. We'll get to that later!"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "\u001b[33;01m➡ models: Cem, Waves\u001b[39;49;00m\n"
     ]
    }
   ],
   "source": [
    "from pymt import models"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "cem, waves = models.Cem(), models.Waves()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Even though we can't run our waves model yet, we can still get some information about it. *Just don't try to run it.* Some things we can do with our model are get the names of the input variables."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "('sea_surface_water_wave__min_of_increment_of_azimuth_angle_of_opposite_of_phase_velocity',\n",
       " 'sea_surface_water_wave__azimuth_angle_of_opposite_of_phase_velocity',\n",
       " 'sea_surface_water_wave__mean_of_increment_of_azimuth_angle_of_opposite_of_phase_velocity',\n",
       " 'sea_surface_water_wave__max_of_increment_of_azimuth_angle_of_opposite_of_phase_velocity',\n",
       " 'sea_surface_water_wave__height',\n",
       " 'sea_surface_water_wave__period')"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "waves.get_output_var_names()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "('sea_surface_water_wave__azimuth_angle_of_opposite_of_phase_velocity',\n",
       " 'land_surface_water_sediment~bedload__mass_flow_rate',\n",
       " 'sea_surface_water_wave__period',\n",
       " 'sea_surface_water_wave__height',\n",
       " 'land_surface__elevation',\n",
       " 'model__time_step')"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "cem.get_input_var_names()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We can also get information about specific variables. Here we'll look at some info about wave direction. This is the main input of the Cem model. Notice that BMI components always use [CSDMS standard names](http://csdms.colorado.edu/wiki/CSDMS_Standard_Names). The CSDMS Standard Name for wave angle is,\n",
    "\n",
    "    \"sea_surface_water_wave__azimuth_angle_of_opposite_of_phase_velocity\"\n",
    "\n",
    "Quite a mouthful, I know. With that name we can get information about that variable and the grid that it is on (it's actually not a one)."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "OK. We're finally ready to run the model. Well not quite. First we initialize the model with the BMI **initialize** method. Normally we would pass it a string that represents the name of an input file. For this example we'll pass **None**, which tells Cem to use some defaults."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [],
   "source": [
    "args = cem.setup(number_of_rows=100, number_of_cols=200, grid_spacing=200.0)\n",
    "cem.initialize(*args)\n",
    "args = waves.setup()\n",
    "waves.initialize(*args)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Here I define a convenience function for plotting the water depth and making it look pretty. You don't need to worry too much about it's internals for this tutorial. It just saves us some typing later on."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [],
   "source": [
    "def plot_coast(spacing, z):\n",
    "    import matplotlib.pyplot as plt\n",
    "\n",
    "    xmin, xmax = 0.0, z.shape[1] * spacing[1] * 1e-3\n",
    "    ymin, ymax = 0.0, z.shape[0] * spacing[0] * 1e-3\n",
    "\n",
    "    plt.imshow(z, extent=[xmin, xmax, ymin, ymax], origin=\"lower\", cmap=\"ocean\")\n",
    "    plt.colorbar().ax.set_ylabel(\"Water Depth (m)\")\n",
    "    plt.xlabel(\"Along shore (km)\")\n",
    "    plt.ylabel(\"Cross shore (km)\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "It generates plots that look like this. We begin with a flat delta (green) and a linear coastline (y = 3 km). The bathymetry drops off linearly to the top of the domain."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXkAAADuCAYAAADcF3dyAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAHiZJREFUeJzt3XmUJHWZ7vHv082mIAq2IAraiCAiKmi74sywXe1RFFBxhuFy1eEOqDCKioqcwQ10wKO4izSK4IgLigiKosgi4wZ0C9Jsigoi0Je2h0VUtu567h8RCdnZWZWRVZmRUcnzOadOVUZlRLwVp/utX73xi/cn20RExHiaM+oAIiJieJLkIyLGWJJ8RMQYS5KPiBhjSfIREWMsST4iYowlyUdEjLEk+YiIMZYkHxExxtYadQAREU2zcOFCr1ixouf7lixZ8gPbC2sIadqS5CMiOqxYsYKLL7205/vWmjNnXg3hzEiSfEREBwOrJsajr1eSfEREBxtWjUnzxiT5iIg1OCP5iIhxNWG4b9XEqMMYiCT5iIguMpKPiBhTJjX5iIgx5iT5iIhxZadcExExtjJPPiJijDmzayIixllq8hERYys1+YiIMZYplBERYy4j+YiIMWWb+yZy4zUiYixlCmVExBhLTT4iYpxldk1ExPhy5slHRIy3cRnJzxl1ABERTdNaNKTXRy+StpB0gaRrJF0l6S3l9o0lnSvpuvLzRsP6WZLkIyK6WDXhnh8VrATebvupwPOBgyVtBxwOnGd7a+C88vVQJMlHRHRo1eR7ffQ8jr3M9i/Lr+8CrgEeD+wJnFK+7RRgryH9KKnJR0SsofrsmnmSFre9XmR7Ubc3SpoP7AhcDGxqexkUvwgkbTKzgCeXJB8R0aGPefIrbC/o9SZJGwCnA4fa/rOkGUZYXZJ8RESHQT4MJWltigR/qu1vlZtvlbRZOYrfDFg+kJN1kZp8REQH24OaXSPgC8A1to9r+9ZZwGvLr18LnDnwH6KUkXxERIcB9q7ZCdgfWCrp8nLbEcAxwGmSDgBuBPYZxMm6SZKPiOhiEOUa2z8BJivA7zbjE1SQJB8R0SErQ0VEjLkk+YiIMTVR3ngdB0nyERFdpAtlRMSYyqIhERHjLDdeIyLG1zgtGjK0J16b0Ec5ImI6Wg9DDaDV8MgNs63ByPsoR0RMhwe0aEgTDC3JN6GPckTEdA2in3wT1FKTH1Uf5YiI6TCzpxzTy9CT/HT7KEs6EDgQYP3113/2tttuO7wgI2JsLFmyZIXtx8zkGHamUFYykz7K5eoqiwCevWCBf3bxJcMMNSLGxLprzf3DII4zLiP5Yc6uGXkf5YiI6Wg9DJWa/NQG1kfZNveunB13siNi9nN61/Q2yD7KE+V0poiIOgxw0ZCRmxVPvBq4byJJPiLqM1vKMb3MiiQfEVGnLBpSM9vcn3JNRNQoI/ka2eTGa0TUZsLmvjHJObMiyU/Y3L1y1ajDiIiHkIzkIyLGVGbX1MxkCmVE1CdtDWo2TovqRsRsMHueaO1lViR5Q2bXRERtUq6JiBhjHqOn7GdFkh+nPhIR0XytBmXjYJYk+fH5rRoRs4CzaEitjLl3VebJR0Q9UpOPiBhzKdfUaMJw95g8YhwRzZcGZTWzM4UyIupjPDbtzWdFkp8gs2sioj6pyUdEjLOUa+qVVsMRUafMk69ZHoaKiDqlXFOzrPEaEXUap4HlnFEHEBHRRKvsnh9VSDpJ0nJJV7Zte5+kmyVdXn68dIr950jaUdLLJO0qadN+fo6eI3lJLwD+N/B3wGbA3cCVwNnAl23f2c8Jp2PC5p77Vw77NBERwMDnyZ8MfBr4Usf2j9n+yGQ7SdoKeBewO3Ad8CdgPWAbSX8DTgBOsT3lnxxTJnlJ3wduAc4EPggsb50E2AU4U9Jxts+a6jgzVfSuGY/6WEQ03yBvvNq+SNL8aex6NHA8cJC9ejCSNgH+BdgfOGWqg/Qaye9ve0XHtr8Avyw/PippXj9RT0dWhoqIelUux8yTtLjt9SLbiyqe5BBJ/wdYDLzd9u2rRWDvO2l09nLg41VOMmWS70zwkjZs38f2bV1+CbTeexKwB7Dc9vbltvcB/0bxZwfAEba/VyXQiIi69DG7ZoXtBdM4xfHAUeWpjgI+CvxrtzdKmgu8DJjP6vn3uConqjS7RtJBwAco6vGtn9zAk6bY7WSmUYfqJsv/RUStDKuGmHNs39r6WtKJwHenePt3gHuApUDfQVWdQnkY8LTJRu3dzKAOteaxSLkmImo2xIehJG1me1n5cm+KySyT2dz2M6Z7rqpTKH8H/G26J+lwiKQrymlFG032JkkHSlosafE9d94+2dsiIgbPLtrf9vqoQNJXgZ8DT5F0k6QDgA9LWirpCopJLG+d4hDfl/Ti6f4oVUfy7wZ+Juli4N7WRttv7vN8letQ5c2LRQDznvxUZyQfEbUa3OyabjdQv9DHIX4BnCFpDnA/oOKw3rDKzlWT/AnA+UyzJtTSZx2qbb+UayKiZs1pa/BR4AXA0s6plFVUTfIrbb+t34N36rMO9YAJzN0rs/xfRNTENCnJXwdcOZ0ED9WT/AWSDqS4y9terrltsh3KOtTOFPNIbwLeC+wsaQeKS3gDcNB0go6IGC5Dc6oHy4ALy4dT2/Pv4KZQUjxZBUVt/oFzMMUUygHUoR4wkVbDEVG35rQavr78WKf86EvVJP+kLo/VrtfvyabLThfKiKiRaUySt/3+mexfNcl/gbZZMJLWB84CdpvJyasyzhqvEVGvEdfkJS0CPmV7aZfvrQ/8E3Cv7VOnOk7VJH+zpONtv7Gc2342cGK/QUdEzA5uwkj+s8CRkp5OMUml1YVya2BD4CRgygQPFZO87SMlHSvpc8CzgWNsnz7dyPuVKZQRUauiec1oQ7AvB14jaQNgAQ+2er/G9q+rHqdXq+FXtr28BDiy/GxJr7T9rb4jnwbbufEaEfVqyBRK238BLpzu/r1G8i/veH0ZsHa53UAtSX6CjOQjokbNmic/I71aDb++rkAiIpqjETX5gehVrvkP4DOdzezbvr8r8HDbldoTTJedJ14jomYPhZE8Ra+a70q6h2IlqPa7uzsAPwI+NNQIKa71fauS5COiRg0ZyUvaBngH8ERWXzRk1yr79yrXnEmxjuvWwE4Ud3f/DHwZOND23dOMuy9FP/lmXPCIeAhowOyaNt8APkcxbb3v0W7VKZTXUTTJiYh4CKjeL74GK20fP92dqz4MNVLO8n8RUacGtDWQtHH55XckvQk4g4oNItvNkiSfKZQRUbPRj+SXUPy6Ufn6HW3f67XG9gNmR5InI/mIqNmIR/K2t4SiGaTte9q/10+DyEprvEraRtJ5kq4sXz+jnF4ZETF+XPaT7/VRj59V3NZV1ZH8iRR/KpwAYPsKSV8Bjq56opmw4f7Mk4+IOo2+Jv9Y4PHAwyTtyINlmw2Bh1c9TtUk/3Dbl0hq37ay6klmrPVbNSKiLqOvyb8EeB2wOdC+CtRdwBFVD1I1ya+QtBVFsR9Jr6ZYkqo+SfIRUZcGzK6xfQpwiqRXzaTrb9UkfzCwCNhW0s0US1HtN92TRkQ03uhH8i0XSvok8CKKXz8/AT5g+3+q7NwzyUuaAyywvXu5Gskc23fNJOK+GUir4YioTaMalH0NuAh4Vfl6P+DrwO5Vdu6Z5G1PSDoEOM32X6cb5YykJh8RdWpWW4ONbR/V9vpoSXtV3blqueZcSYdR/PZ4INFXfeJqILKQd0TUqTnlmgsk/TNwWvn61RRLsFZSNcm3FvE+uG1b5SeuIiJmlQbceG1zEPA2isaQBuYCf5X0NsC2N5xq56oNyracaZQzknJNRNSqOQ3KbD9iJvtXSvKS1gbeCPx9uelC4ATb90+xz0nAHsBy29uX2zamKPnMB24AXjPZgiSrMXB/HoaKiBo1JMmreEBpP2BL20dJ2gLYzPYlVfavWq45nmJt18+Wr/cvt/3fKfY5Gfg08KW2bYcD59k+RtLh5et39T59RvIRUaNmlWs+S7HU9a7AUcBfgM8Az6myc9Uk/xzbz2x7fb6kX021g+2LJM3v2LwnsHP59SkUfxFUSPIRETVrzsDyebafJekyANu3S1qn6s5Vk/wqSVvZ/h2ApCcxjRVKgE1tLwOwvUzSJpX2atZ0pogYe42aJ3+/pLk82HHgMRQj+0qqJvl3UEzj+T1Fk5wnAq/vM9C+SDoQOBCA9TeCrPEaEXUxjanJA5+kWDBkE0kfpJhCWbkLcNXZNeeV67w+hSLJX2v73h67dXOrpM3KUfxmwPIpzrmIopUCevQTTNZ4jYg6NWQkb/tUSUuA3Sjy7162r6m6fz+LhjybYlbMWsAzJWH7S1PvsoazgNcCx5Sfz+xz/4iIejRkJC/p6cC2FIPia/pJ8FB9CuV/AVsBl/NgLd6sPnOmc5+vUtxknSfpJuC9FMn9NEkHADcC+1SKMjX5iKhTA57NkfRIioHwFsAVFKP4p0u6EdjT9p+rHKfqSH4BsJ1d/e8X2/tO8q3dqh6j7WiQRUMiok6jH8kfBSwGdrU9AVDegP1P4IPAv1c5SNUkfyXwWOruId+SkXxE1GmA8+Rn8GDo7sAzWgkewPYqSUcAS6uef8o1XiV9R9JZwDzgakk/kHRW66PqSSIiZp0J9/6o5mRgYce21oOhWwPnla873Wd7jRX4ym2VJ770Gsl/pOqBhqoB9bGIeIgZ0Eh+Bg+GrtextmuLgHWrnn/KJG/7xwDlYiF3l73lt6G40/v9qicZiCT5iKiLK4/U50la3PZ6UTn9u5cqD4YuY/W1Xdv9vyrBQfWa/EXA30naiOJPi8XAP1HXEoAZyUdE3aqtYbHC9oJhnN72LoM4zpQ1+Tay/TfglcCnbO8NPG0QAURENNLgavLd3Fo+EEqvB0NnqupIXpJeQDFyP6DcNnc4IXWRNV4jok7Db2tQ24OhVZP8W4B3A2fYvqpsUHbBsILqKuWaiKjN4BqUzeTB0LKX/Oa2/zjd81ftXXMRRV2+9fr3wJune9K+OQ9DRUTNBje7ZtoPhtq2pG9TtJWZlqo1+YiIh45WuWZ4Nfl+/EJSpQVCuumnQdno5InXiKhbc3LOLsAbJN0A/JVinrxtP6PKzrMjyeP0k4+I+lSfJ1+Hf5zJzlW7UH4YOBq4GzgHeCZwqO0vz+TklRnSTz4iatWcfvJ/kPQiYGvbXyxXhtqg6v5Va/IvLtta7gHcBGxDsVpURMR4snt/1EDSeylaHry73LQ2UHmAXbVcs3b5+aXAV23fVszsqUmeeI2IujWnXLM3sCPwSwDbt0h6RNWdqyb570i6lqJc86byz4V7+o00ImJWaNbA8r5yKmVrIe/1+9m5UrnG9uHAC4AFtu+nuMO7Z7+RRkTMGg0p11A8NHUC8ChJ/wb8CPh81Z0rJXlJ+wAry4b1/0FRD3rcdKKNiGi8Bs2Tt/0R4JvA6cBTgPfY/mTV/aveeD3S9l3lHd6XUPQ/Pr7fYCMiZo2GjOQlHWv7XNvvsH2Y7XMlHVt1/6pJvjVJ/WXA8bbPBNbpN9iIiFmjISN54H912VZ57nzVG683lzWh3YFjJa1LnS0RHnUrvOLDtZ0uImaxywZwjAY8DCXpjcCbgCdJuqLtW48Aflr1OFWT/Gso1ij8iO07yv7HmScfEeNr9LNrvkKxAt9/svoasHfZvq3qQap2ofybpN8BL5H0EuC/bf+wn2gjImaVET/xavtO4E5gX4ByicD1gA0kbWD7xirHqTq75i3AqcAm5ceXJf37dAKPiGg806Qbry+XdB1wPfBj4Ab6WGO7arnmAOB5tv9anvRY4OfAp/qKNiJiVhh9Tb7N0cDzgR/Z3lHSLpSj+yoqL//HgzNsKL+edl+DsmXmXeVxVg5rIdyIiGlrTpK/3/b/SJojaY7tC/qZQlk1yX8RuFjSGeXrvYAv9Btph11sr5jhMSIiBq9VrmmGOyRtQLE636mSlgMrq+5c9cbrcZIuBF5EMYJ/ve1BTFSKiGim0c+uadmTolfYW4H9gEcCH6i6c88kL2kOcIXt7Sm7oA2AgR+WDXdOsL2oy3kPBA4Eih8pIqIuzZgnfyjFfPjLbLfK5af0e5yeSd72hKRfSXpC1Sk7FexUtsvcBDhX0rXlYuHt510ELALQ49SYv5si4iFi9OWazYFPANuWD0P9jCLp/3zg8+SBzYCrJF1C0YESANuvqB7vg2zfUn5eXtb5n0tRb4qIaIYRj+RtHwYgaR1gAfBC4F+BEyXdYXu7KsepmuTfP60ouyh7Ic8pG56tD7yYPupLERG1GP1IvuVhwIYUhetHArcAS6vuPGWSl/RkYFPbP+7Y/vfAzX2HWtgUOKNcWWot4Cu2z5nmsSIihsAwt8KN11W93zJdkhYBT6OYbn4xRbnmONu393OcXiP5jwNHdNn+t/J7L+/nZAC2f0+xEHhERDMJmFNhJD/EJA88AVgXuI5iUH0TcEe/B+mV5OfbvqJzo+3Fkub3e7KIiFljxPM9bC9UUfJ4GkU9/u3A9pJuo7j5+t4qx+mV5Neb4nsPqxRpRMRsVGUkP2S2DVwp6Q6KZmV3AntQTFaplOR7NSi7tFxTcDWSDgCW9BduRMRs4WIk3+tjiCS9WdLXJP2RYvbhHsCvgVcCG1c9Tq+R/KEUN0n348GkvoBiVai9+446ImI2ECMv1wDzKdZ2favtZdM9yJRJ3vatwAvLrmfbl5vPtn3+dE8YETErrDXatga23zaI41TtXXMBcMEgThgR0XjNGMkPRNWHoSIiHkLciBuvg5AkHxHRTUbyERFjqurDULNAknxERDdzGtNPfkaS5CMi1pCafETE+Brg7JpRr2mdJB8R0c1gR/IjW9M6ST4iopsxKdf06l0TEfHQ0yrX9O5dM0/S4raPA7scrbWm9ZJJvj9UGclHRKyh4qIhsKJCjb3nmtbDlJF8RESn1jz5Xh8VtK9pDbTWtK5NknxERDcDaDUsaX1Jj2h9TbGm9ZVDjnw1KddERHQzmBuvI1/TOkk+ImINg1kUpAlrWifJR0R0Su+aiIgxV212TeMlyUdEdMqiIRER42x8GpSNZAqlpIWSfi3pt5IOH0UMERFTGsAUyiaoPclLmgt8BvhHYDtgX0nb1R1HRMSUxiTJj6Jc81zgt+XUIiR9DdgTuHoEsURErEnAWuNx43UU5ZrHA39se31TuS0ioiEqjOIzkp+Uumxb42qV3dpaHdvu5f31PgrcxTxgJP2gOzQhjibEAM2IowkxQDPiaEIMAE+Z8REyT35GbgK2aHu9OXBL55tsLwIWAUhaXPdqKp2aEENT4mhCDE2JowkxNCWOJsTQimMwBxqPJD+Kcs2lwNaStpS0DvDPwFkjiCMiYnID6kI5arWP5G2vlHQI8ANgLnCS7avqjiMiYlJ5GGpmbH8P+F4fuywaVix9aEIM0Iw4mhADNCOOJsQAzYijCTHAQOKovGhI48kej99WERGDog0fZhY8qfcbL7h6SRPuQ0wlbQ0iIroZk3JNo1eGakr7A0k3SFoq6fKB3bmvdt6TJC2XdGXbto0lnSvpuvLzRiOI4X2Sbi6vx+WSXjrkGLaQdIGkayRdJekt5fa6r8VkcdR2PSStJ+kSSb8qY3h/uX1LSReX1+Lr5aSGoZkijpMlXd92LXYYZhzlOedKukzSd8vXg7kWY3LjtbFJvoHtD3axvUPNf5qdDCzs2HY4cJ7trYHzytd1xwDwsfJ67FDeYxmmlcDbbT8VeD5wcPlvoe5rMVkcUN/1uBfY1fYzgR2AhZKeDxxbxrA1cDtwwBBjmCoOgHe0XYvLhxwHwFuAa9peD+BajM/DUI1N8rS1P7B9H9Bqf/CQUa7oflvH5j2BU8qvTwH2GkEMtbK9zPYvy6/vovgP/XjqvxaTxVEbF/5Svly7/DCwK/DNcnsd12KyOGolaXPgZcDny9diENdigAt5j1qTk3yT2h8Y+KGkJeWTuKO0qe1lUCQdYJMRxXGIpCvKcs5QyyTtJM0HdgQuZoTXoiMOqPF6lOWJy4HlwLnA74A7bK8s31LL/5XOOGy3rsUHy2vxMUnrDjmMjwPvBFpTYR7NoK7F3IneH7NAk5N8pfYHNdnJ9rMoSkcHS/r7EcXRFMcDW1H8mb4M+GgdJ5W0AXA6cKjtP9dxzopx1Ho9bK+yvQPF0+LPBZ7a7W3DjKFbHJK2B94NbAs8B9gYeNewzi9pD2C57SXtm7uF2v/BSbmmBpXaH9TB9i3l5+XAGRT/sUblVkmbAZSfl9cdgO1by//gE8CJ1HA9JK1NkVhPtf2tcnPt16JbHKO4HuV57wAupLg/8ChJrdlytf5faYtjYVnSsu17gS8y3GuxE/AKSTdQlHN3pRjZD+BaVCjVpFwzY41ofyBpfUmPaH0NvBhG2iztLOC15devBc6sO4BWYi3tzZCvR1ln/QJwje3j2r5V67WYLI46r4ekx0h6VPn1w4DdKe4NXAC8unxbHdeiWxzXtv3SFUUtfGjXwva7bW9uez5Ffjjf9n4M6lqMyUi+sfPkG9T+YFPgjOLfLGsBX7F9Th0nlvRVYGdgnqSbgPcCxwCnSToAuBHYZwQx7FxOjTNwA3DQMGOgGLHtDywta8AAR1DztZgijn1rvB6bAaeUs8/mAKfZ/q6kq4GvSToauIzil9EwTRbH+ZIeQ1HwuBx4w5Dj6OZdzPRajFEXyjzxGhHRQRuta3Z7bO83nn5jnniNiJh1ZlE5ppck+YiIbsakXJMkHxHRTUbyERFjaoxuvCbJR0R0MyYj+SbPk48RkrS3JEvatm3bfLV1o6wxlp1bHQZrONdekt5Tfn2ypFf32meS4zxGUi1TbWMYnLYGMfb2BX5C8ZDJrFbO5a7qncBnZ3pO238ClknaaabHihFIg7IYZ2Vvlp0oWrR2TfJlP/Evquizf5mkXcrtr5P0LUnnlP28P9y2zwGSfiPpQkknSvp0l+P+Q1sv8staTxsDG0j6pqRrJZ1aPlGJpN3K9y0tm4OtW26/QdJ7JP0E2EfSVmVMSyT9d/tfKG3n3ga41/aKLt87qhzZzymP/SFJP5e0WNKzJP1A0u8ktT/8821gv2pXPRonT7zGGNsLOMf2byTdJulZrRa7bQ4GsP30MmH+sEySUDTq2pGi5/ivJX0KWAUcCTwLuAs4H/hVl3MfBhxs+6flL5t7yu07Ak+j6EPyU2AnFQu4nAzsVsb6JeCNFP1LAO6x/SIASecBb7B9naTnUYzWd+04905A589J+YvqkcDrbbv8/fJH2y+Q9LEyhp2A9YCrgM+Vuy4Gju7yM8ZsMEtG6r1kJB/d7EvR8Iny875d3vMi4L8AbF8L/AFoJfnzbN9p+x7gauCJFI2qfmz7Ntv3A9+Y5Nw/BY6T9GbgUW0tYy+xfVPZBOxyYD7wFOB6278p33MK0N4h9OvwwF8mLwS+UbYjOIHisfxOmwF/6th2ZBnHQV798fBWH6WlwMW27ypLNPe0erpQNEx73CQ/ZzRZlVF8RvIxG0l6NMUId3tJpugbZEnv7HzrFIe5t+3rVRT/zqZ6/wNsHyPpbOClwC8k7T6DY/61/DyHosd4r6Xo7qYYsbe7FHi2pI1tty+e0opnoiO2CR78f7VeecyYjTKSjzH1auBLtp9oe77tLYDrKUbu7S6irDeXZZonAL+e4riXAP8gaSMVbWBf1e1NkrayvdT2sRTljjVq522uBeZLenL5en/gx51vKnu+Xy9pn/IckvTMLse7Bnhyx7ZzKBqhnd12f6CqbRhtx9KYLpHZNTG29qXomd/udOBfOrZ9FpgraSlFWeR1ZQ/xrmzfDHyIYiWlH1GUce7s8tZDJV0p6VcUo+DvT3HMe4DXU5RhllKMoj83ydv3Aw4oj3sV3ZeSvAjYsXVTt+0836DoE3+Wira6Ve0CnN3H+6NJxqRcky6UURtJG9j+SzmSP4OifXTnL5SRkvQJ4Du2fzSAY10E7Gn79plHFnXSJnPNazbo/cbP/LlnF0pJC4FPUJQ+P2/7mIEEWVFG8lGn95U3Pq+kKAF9e8TxdPMh4OEzPUjZU/24JPhZbAAj+fIZjc9QLB26HcXaA9sNOfLV5MZr1Mb2YaOOoRfbtzKAFcjKmTZN/CUWVQyud81zgd/a/j2ApK9RlAqvHsTBq0iSj4joptqN1Xnl8xoti2wvanv9eOCPba9vAp43gOgqS5KPiOhU/cbqih41+W7TfGu9EZokHxHRzWDKNTcBW7S93pziqe3a5MZrREQ3g5lCeSmwtaQtJa1D0Qtqxvd8+pGRfEREpwHdeLW9UtIhwA8oplCeZPuqGR+4D0nyERHdDOhhJ9vfA743kINNQ5J8RMQaPGvaFvSSJB8R0SlrvEZEjLlZ0puml/SuiYjoUK7PO6/CW1fYXjjseGYiST4iYoxlnnxExBhLko+IGGNJ8hERYyxJPiJijCXJR0SMsST5iIgxliQfETHGkuQjIsZYknxExBj7/3gDpqFZCSU/AAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 2 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "grid_id = cem.get_var_grid(\"sea_water__depth\")\n",
    "spacing = cem.get_grid_spacing(grid_id)\n",
    "shape = cem.get_grid_shape(grid_id)\n",
    "z = np.empty(shape)\n",
    "cem.get_value(\"sea_water__depth\", out=z)\n",
    "plot_coast(spacing, z)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Allocate memory for the sediment discharge array and set the discharge at the coastal cell to some value."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [],
   "source": [
    "qs = np.zeros_like(z)\n",
    "qs[0, 100] = 750"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The CSDMS Standard Name for this variable is:\n",
    "\n",
    "    \"land_surface_water_sediment~bedload__mass_flow_rate\"\n",
    "\n",
    "You can get an idea of the units based on the quantity part of the name. \"mass_flow_rate\" indicates mass per time. You can double-check this with the BMI method function **get_var_units**."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "'kg / s'"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "cem.get_var_units(\"land_surface_water_sediment~bedload__mass_flow_rate\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([ 7.])"
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "waves.set_value(\n",
    "    \"sea_shoreline_wave~incoming~deepwater__ashton_et_al_approach_angle_asymmetry_parameter\",\n",
    "    0.3,\n",
    ")\n",
    "waves.set_value(\n",
    "    \"sea_shoreline_wave~incoming~deepwater__ashton_et_al_approach_angle_highness_parameter\",\n",
    "    0.7,\n",
    ")\n",
    "cem.set_value(\"sea_surface_water_wave__height\", 2.0)\n",
    "cem.set_value(\"sea_surface_water_wave__period\", 7.0)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Set the bedload flux and run the model."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([[ -1. ,  -1. ,  -1. , ...,  -1. ,  -1. ,  -1. ],\n",
       "       [ -1. ,  -1. ,  -1. , ...,  -1. ,  -1. ,  -1. ],\n",
       "       [ -1. ,  -1. ,  -1. , ...,  -1. ,  -1. ,  -1. ],\n",
       "       ..., \n",
       "       [ 22.4,  22.4,  22.4, ...,  22.4,  22.4,  22.4],\n",
       "       [ 22.6,  22.6,  22.6, ...,  22.6,  22.6,  22.6],\n",
       "       [ 22.8,  22.8,  22.8, ...,  22.8,  22.8,  22.8]])"
      ]
     },
     "execution_count": 12,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "for time in range(3000):\n",
    "    waves.update()\n",
    "    angle = waves.get_value(\n",
    "        \"sea_surface_water_wave__azimuth_angle_of_opposite_of_phase_velocity\"\n",
    "    )\n",
    "\n",
    "    cem.set_value(\n",
    "        \"sea_surface_water_wave__azimuth_angle_of_opposite_of_phase_velocity\", angle\n",
    "    )\n",
    "    cem.set_value(\"land_surface_water_sediment~bedload__mass_flow_rate\", qs)\n",
    "    cem.update()\n",
    "\n",
    "cem.get_value(\"sea_water__depth\", out=z)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXkAAADuCAYAAADcF3dyAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAH6lJREFUeJzt3Xm4JGV99vHvPQOIgig4QlDQUQQRUUBxC8YIGINKBNeEEF81vMEFo8Zdr7iCBr0UTVyQMSJjxDVKwF1kkbgEnBFkQDAERUF4mYwIorLNnPv9o6qh6elzunqrrlNzf66rr3O6uque36k58+vnPPXU75FtIiKinZbMOoCIiJieJPmIiBZLko+IaLEk+YiIFkuSj4hosST5iIgWS5KPiGixJPmIiBZLko+IaLHNZh1ARETTHHTQQV63bt3A961evfqbtg+qIaSRJclHRPRYt24d5/7whwPft9mSJctqCGcsSfIRET0MbJhrR12vJPmIiB42bGhJ8cYk+YiIjTg9+YiItpoz3LphbtZhTESSfEREH+nJR0S0lMmYfEREizlJPiKirewM10REtFbmyUdEtJgzuyYios0yJh8R0VoZk4+IaLFMoYyIaLn05CMiWso2t87lwmtERCtlCmVERItlTD4ios0yuyYior2cefIREe3Wlp78klkHEBHRNJ1FQwY9BpG0s6SzJF0i6WJJryi3byfpdEmXlV+3ndbPkiQfEdHHhjkPfFSwHni17YcAjwWOkrQH8AbgDNu7AmeUz6ciST4iokdnTH7QY+Bx7Gts/6j8/kbgEuC+wCHAyvJtK4FDp/SjZEw+ImIj1WfXLJO0quv5Ctsr+r1R0nJgH+BcYAfb10DxQSBp+/ECnl+SfEREjyHmya+zve+gN0naGvgi8Erbv5U0ZoTVJclHRPSY5M1QkjanSPAn2/5SuflaSTuWvfgdgbUTaayPjMlHRPSwPanZNQI+Dlxi+7iul04Dnl9+/3zg1In/EKX05CMiekywds1+wPOANZIuKLe9CTgW+LykI4BfAs+ZRGP9JMlHRPQxieEa298F5huAP3DsBipIko+I6JGVoSIiWi5JPiKipebKC69tkCQfEdFHqlBGRLRUFg2JiGizXHiNiGivNi0aMrU7XptQRzkiYhSdm6EmUGp45qZZ1mDmdZQjIkbhCS0a0gRTS/JNqKMcETGqSdSTb4JaxuRnVUc5ImIUZvEMxwwy9SQ/ah1lSUcCRwJstdVWj9x9992nF2REtMbq1avX2b73OMewM4WyknHqKJerq6wAeOS++/r75543zVAjoiXustnSX0ziOG3pyU9zds3M6yhHRIyiczNUxuQXNrE6yra5Zf3iuJIdEYufU7tmsEnWUZ4rpzNFRNRhgouGzNyiuOPVwK1zSfIRUZ/FMhwzyKJI8hERdcqiITWzzW0ZromIGqUnXyObXHiNiNrM2dzakpyzKJL8nM1N6zfMOoyI2ISkJx8R0VKZXVMzkymUEVGflDWoWZsW1Y2IxWDx3NE6yKJI8obMromI2mS4JiKixdyiu+wXRZJvUx2JiGi+ToGyNlgkSb49n6oRsQg4i4bUyphbNmSefETUI2PyEREtl+GaGs0ZbmrJLcYR0XwpUFYzO1MoI6I+xq0pb74okvwcmV0TEfXJmHxERJtluKZeKTUcEXXKPPma5WaoiKhThmtqljVeI6JObepYLpl1ABERTbTBHvioQtKJktZKuqhr29sk/UrSBeXjqQvsv0TSPpKeJukASTsM83MM7MlLehzwN8CfADsCNwEXAV8FPmX7hmEaHMWczc23rZ92MxERwMTnyZ8EfAj4ZM/299t+73w7SdoFeD3wJOAy4H+BLYHdJP0BOAFYaXvBPzkWTPKSvg5cDZwKvBNY22kE2B84VdJxtk9b6DjjKmrXtGN8LCKab5IXXm2fI2n5CLseAxwPvMi+czCStgf+GngesHKhgwzqyT/P9rqebb8DflQ+3idp2TBRjyIrQ0VEvSoPxyyTtKrr+QrbKyo28jJJ/wdYBbza9m/uFIF92LzR2WuBD1RpZMEk35vgJW3TvY/t6/p8CHTeeyJwMLDW9p7ltrcBf0fxZwfAm2x/rUqgERF1GWJ2zTrb+47QxPHA0WVTRwPvA/623xslLQWeBiznzvn3uCoNVZpdI+lFwDsoxuM7P7mBBy6w20mMMA7VT5b/i4haGTZMMefYvrbzvaSPAV9Z4O1fBm4G1gBDB1V1CuVrgIfO12vvZ4xxqI2PRYZrIqJmU7wZStKOtq8pnz6DYjLLfHay/fBR26o6hfJy4A+jNtLjZZIuLKcVbTvfmyQdKWmVpFU33/Cb+d4WETF5dlH+dtCjAkmfAX4APFjSVZKOAN4jaY2kCykmsfzDAof4uqQnj/qjVO3JvxH4vqRzgVs6G22/fMj2Ko9DlRcvVgAse9BDnJ58RNRqcrNr+l1A/fgQh/gv4BRJS4DbABWH9TZVdq6a5E8AzmTEMaGOIcehuvbLcE1E1Kw5ZQ3eBzwOWNM7lbKKqkl+ve1XDXvwXkOOQ91uDnPT+iz/FxE1MU1K8pcBF42S4KF6kj9L0pEUV3m7h2uum2+HchzqiRTzSK8C3go8UdLeFKfwCuBFowQdETFdhuaMHlwDnF3enNqdfyc3hZLiziooxuZvb4MFplBOYBzqdnMpNRwRdWtOqeGfl48tysdQqib5B/a5rXbLYRsblZ0qlBFRI9OYJG/77ePsXzXJf5yuWTCStgJOAw4cp/GqjLPGa0TUa8Zj8pJWAB+0vabPa1sBfwncYvvkhY5TNcn/StLxtl9Szm3/KvCxYYOOiFgc3ISe/EeAN0t6GMUklU4Vyl2BbYATgQUTPFRM8rbfLOndkj4KPBI41vYXR418WJlCGRG1KorXzDYE+wLguZK2BvbljlLvl9j+adXjDCo1/Myup+cBby6/WtIzbX9p6MhHYDsXXiOiXg2ZQmn7d8DZo+4/qCf/Fz3Pzwc2L7cbqCXJz5GefETUqFnz5McyqNTwC+sKJCKiORoxJj8Rg4Zr/hH4cG8x+67XDwDuZrtSeYJR2bnjNSJqtin05Clq1XxF0s0UK0F1X93dG/g28K6pRkhxrm/dkCQfETVqSE9e0m7Aa4H7c+dFQw6osv+g4ZpTKdZx3RXYj+Lq7m+BTwFH2r5pxLiHUtSTb8YJj4hNQANm13T5AvBRimnrQ/d2q06hvIyiSE5ExCager34Gqy3ffyoO1e9GWqmnOX/IqJODShrIGm78tsvS3opcAoVC0R2WyRJPlMoI6Jms+/Jr6b4uFH5/LVdrw1aY/t2iyPJk558RNRsxj152w+Aohik7Zu7XxumQGSlNV4l7SbpDEkXlc8fXk6vjIhoH5f15Ac96vH9itv6qtqT/xjFnwonANi+UNKngWOqNjQOG27LPPmIqNPsx+T/CLgvcFdJ+3DHsM02wN2qHqdqkr+b7fMkdW9bX7WRsXU+VSMi6jL7Mfk/B14A7AR0rwJ1I/CmqgepmuTXSdqFYrAfSc+mWJKqPknyEVGXBsyusb0SWCnpWeNU/a2a5I8CVgC7S/oVxVJUh4/aaERE482+J99xtqR/AR5P8fHzXeAdtn9dZeeBSV7SEmBf208qVyNZYvvGcSIemoGUGo6I2jSqQNlngXOAZ5XPDwc+Bzypys4Dk7ztOUkvAz5v+/ejRjmWjMlHRJ2aVdZgO9tHdz0/RtKhVXeuOlxzuqTXUHx63J7oq95xNRFZyDsi6tSc4ZqzJP0V8Pny+bMplmCtpGqS7yzifVTXtsp3XEVELCoNuPDa5UXAqygKQxpYCvxe0qsA295moZ2rFih7wLhRjiXDNRFRq+YUKLN993H2r5TkJW0OvAR4QrnpbOAE27ctsM+JwMHAWtt7ltu2oxjyWQ5cATx3vgVJ7sTAbbkZKiJq1JAkr+IGpcOBB9g+WtLOwI62z6uyf9XhmuMp1nb9SPn8eeW2/7vAPicBHwI+2bXtDcAZto+V9Iby+esHN5+efETUqFnDNR+hWOr6AOBo4HfAh4FHVdm5apJ/lO29up6fKenHC+1g+xxJy3s2HwI8sfx+JcVfBBWSfEREzZrTsXyM7UdIOh/A9m8kbVF156pJfoOkXWxfDiDpgYywQgmwg+1rAGxfI2n7Sns1azpTRLReo+bJ3yZpKXdUHLg3Rc++kqpJ/rUU03h+RlEk5/7AC4cMdCiSjgSOBGCrbSFrvEZEXUxjxuSBf6FYMGR7Se+kmEJZuQpw1dk1Z5TrvD6YIslfavuWAbv1c62kHcte/I7A2gXaXEFRSgHd634ma7xGRJ0a0pO3fbKk1cCBFPn3UNuXVN1/mEVDHkkxK2YzYC9J2P7kwrts5DTg+cCx5ddTh9w/IqIeDenJS3oYsDtFp/iSYRI8VJ9C+W/ALsAF3DEWb+48c6Z3n89QXGRdJukq4K0Uyf3zko4Afgk8p1KUGZOPiDo14N4cSfeg6AjvDFxI0Yt/mKRfAofY/m2V41Ttye8L7GFX//vF9mHzvHRg1WN0HQ2yaEj0uvxa2GWHWUcRbTX7nvzRwCrgANtzAOUF2H8C3gn8fZWDVE3yFwF/RN015DvSk4/5XH5t8XX5vWcbR7TLBOfJj3Fj6JOAh3cSPIDtDZLeBKyp2v6Ca7xK+rKk04BlwE8kfVPSaZ1H1UYiIhadOQ9+VHMScFDPts6NobsCZ5TPe91qe6MV+MptlSe+DOrJv7fqgaaqAeNjMWNX/hp2vtf8r3d+P64s11FY6L0RVUyoJz/GjaFb9qzt2iHgLlXbXzDJ2/4OQLlYyE1lbfndKK70fr1qIxORJL9purrrL9iFfgd6X7vy13CfbacTU7SfK/fUl0la1fV8RTn9e5AqN4Zew53Xdu32/6oEB9XH5M8B/kTSthR/WqwC/pK6lgBMTz5g49+Bp78HTntd/9fm2xZRVbU1LNbZ3ncazdvefxLHWXBMvots/wF4JvBB288AHjqJACIiGmlyY/L9XFveEMqgG0PHVbUnL0mPo+i5H1FuWzqdkPrIGq+bpl/3LCXc/TvQee3p7ym+fu9oNnLtDXCvsUpxx6Zq+mUNarsxtGqSfwXwRuAU2xeXBcrOmlZQfeVP703L9X2WEx5mTL7KPhHzmlyBsnFuDC1rye9k+8pR269au+YcinH5zvOfAS8ftdGhOTdDBcXvwAHv6P/aE94KZ75l4+2dD4utt5xeXNFOk5tdM/KNobYt6T8oysqMpOqYfETEpqMzXDO9Mflh/JekSguE9DNMgbLZyR2vAYN/B0YZzomYT3N+Z/YHXizpCuD3FPPkbfvhVXZeHEkep558DP4d+LNj7phS2eumW2GL+uYKxCJXfZ58HZ4yzs5Vq1C+BzgGuAn4BrAX8Erbnxqn8coM3Nr1H3zpGKNMvZ/OS5f0/8Qep4352qp67FH3a7tbx/ygv3VDM87hoH/fDXMbx1n197b7OHXpxNGv3WnEOMy5GEdz6sn/QtLjgV1tf6JcGWrrqvtXPdNPLstaHgxcBexGsVpUREQ72YMfNZD0VoqSB28sN20OVO5gVx2u2bz8+lTgM7avK2b21KT3ZE7yU3sWU+9GPXZzxghnozMnvsp75hu2afI57I5tMV5/GKbdOv4Pj6s5wzXPAPYBfgRg+2pJlW8AqZrkvyzpUorhmpeWfy7cPGykERGLQrNKqdxaTqXsLOS91TA7Vxqusf0G4HHAvrZvo7jCe8iwkUbUpkqvP2IhDRmuobhp6gTgnpL+Dvg28K9Vd66U5CU9B1hfFqz/R4rxoPuMEm1EROM1aJ687fcC/w58EXgw8Bbb/1J1/6oXXt9s+8byCu+fU9Q/Pn7YYCNGNkrPPL35GEdDevKS3m37dNuvtf0a26dLenfV/auOyXfmrj0NON72qZLeNmywjdHvP/98F+om2U6VNuqKrV9702xnVOMm6kEXYmehu0Ry1fd3a9LPMp9pxjzod+L8CbXTnAuvf8bGC4o8pc+2vlRlbW5JXwF+RbHm4CMpLsCeZ3uvoUIdke4jc+QQO/T7hZpGr27Udnr3m1aPc1rt1JVkZtkTH+d3aNTz3r3fLH72+f5d64xl0O9WlVjezupxa7xrs2XmnhUuO/76xLHbmjcG6SXAS4EHApd3vXR34Hu2/6bKcar25J9LsUbhe21fX9Y/zjz5iGiv2c+u+TTFCnz/xJ3XgL3R9nVVD1KpJw8gaS/gT8qn/2n7x1UbGdfQPfmI2HRNqie/9cGD33jDyqn15HuVSwTeXk7V9i+r7Fd1ds0rgJOB7cvHpyT9/QhxRkQ0n2nShde/kHQZ8HPgO8AVDLHGdtXhmiOAx9j+fdnou4EfAB8cKtqIiEWhUQXKjgEeC3zb9j6S9gfmq1G/kcrL/3HHDBvK70eua1CWzLyxPM76uv7ciYiorDlJ/jbbv5a0RNIS22dNYwrlJ4BzJZ1SPj8U+PiwkfbY3/a6MY8RETF5neGaZrhe0tYUq/OdLGktsL7qzlWX/ztO0tnA4yl68C+0PanZqBERzTP72TUdh1DUCvsH4HDgHsA862BubGCSl7QEuND2npRV0CbAwLfKgjsn2F7Rp90joZxTc48JtRoRUUUDFg2R9Erge8D5tjvD5SuHPc7AJG97TtKPJd2v6pSdCvYry2VuD5wu6dJysfDudlcAK6CcQhkRUafZD9fsBPwzsLukC4HvUyT9HwwzT77qmPyOwMWSzqOoQAmA7adXj/cOtq8uv64tx/kfTTHeFBHRDDPuydt+DYCkLYB9gT8G/hb4mKTrbe9R5ThVk/zbR4qyj7IW8pKy4NlWwJMZYnwpIqIWs+/Jd9wV2IZi4PoewNXAmqo7L5jkJT0I2MH2d3q2P4Gils0odgBOKVeW2gz4tO1vjHisiIgpMCytcOF1zGWHFyJpBfBQiunm51IM1xxn+zfDHGdQT/4DwJv6bP9D+dpfDNMYgO2fUSwEHhHRTAKWVOjJTzHJA/cD7gJcRtGpvgq4ftiDDEryy21f2LvR9ipJy4dtLCJi0ZjxfA/bB6kY8ngoxXj8q4E9JV1HcfH1rVWOMyjJb7nAa3etFGlExGJUpSc/ZS4qSF4k6XrghvJxMMVklUpJflCBsh+WawreiaQjgNXDhRsRsVi46MkPekyRpJdL+qykKylmHx4M/BR4JrBd1eMM6sm/kuIi6eHckdT3BbYAnjF01BERi4GY+XANsJxibdd/sH3NqAdZMMnbvhb447Lq2Z7l5q/aPnPUBiMiFoXNZlvWwParJnGcqrVrzgLOmkSDERGN14ye/ERUvRkqImIT4kZceJ2EJPmIiH7Sk4+IaKmqN0MtAknyERH9LGlMPfmxJMlHRGwkY/IREe01wdk1s17TOkk+IqKfyfbkZ7amdZJ8REQ/LRmuGVS7JiJi09MZrhlcu2aZpFVdjyP7HK2zpvXqeV6fqvTkIyI2UnHREFhXYYx94JrW05SefEREr848+UGPCrrXtAY6a1rXJkk+IqKfCZQalrSVpLt3vqdY0/qiKUd+JxmuiYjoZzIXXme+pnWSfETERiazKEgT1rROko+I6JXaNRERLVdtdk3jJclHRPTKoiEREW3WngJlM5lCKekgST+V9D+S3jCLGCIiFjSBKZRNUHuSl7QU+DDwFGAP4DBJe9QdR0TEglqS5GcxXPNo4H/KqUVI+ixwCPCTGcQSEbExAZu148LrLIZr7gtc2fX8qnJbRERDVOjFpyc/L/XZttHZKqu1dSq23cLb670VuI9lwEzqQfdoQhxNiAGaEUcTYoBmxNGEGAAePPYRMk9+LFcBO3c93wm4uvdNtlcAKwAkrap7NZVeTYihKXE0IYamxNGEGJoSRxNi6MQxmQO1I8nPYrjmh8Cukh4gaQvgr4DTZhBHRMT8JlSFctZq78nbXi/pZcA3gaXAibYvrjuOiIh55Wao8dj+GvC1IXZZMa1YhtCEGKAZcTQhBmhGHE2IAZoRRxNigInEUXnRkMaT3Y5Pq4iISdE2dzX7PnDwG8/6yeomXIdYSMoaRET005LhmkavDNWU8geSrpC0RtIFE7tyX63dEyWtlXRR17btJJ0u6bLy67YziOFtkn5Vno8LJD11yjHsLOksSZdIuljSK8rtdZ+L+eKo7XxI2lLSeZJ+XMbw9nL7AySdW56Lz5WTGqZmgThOkvTzrnOx9zTjKNtcKul8SV8pn0/mXLTkwmtjk3wDyx/sb3vvmv80Owk4qGfbG4AzbO8KnFE+rzsGgPeX52Pv8hrLNK0HXm37IcBjgaPK34W6z8V8cUB95+MW4ADbewF7AwdJeizw7jKGXYHfAEdMMYaF4gB4bde5uGDKcQC8Arik6/kEzkV7boZqbJKnq/yB7VuBTvmDTUa5ovt1PZsPAVaW368EDp1BDLWyfY3tH5Xf30jxH/q+1H8u5oujNi78rny6efkwcADw7+X2Os7FfHHUStJOwNOAfy2fi0mciwku5D1rTU7yTSp/YOBbklaXd+LO0g62r4Ei6QDbzyiOl0m6sBzOmeowSTdJy4F9gHOZ4bnoiQNqPB/l8MQFwFrgdOBy4Hrb68u31PJ/pTcO251z8c7yXLxf0l2mHMYHgNcBnakw92JS52Lp3ODHItDkJF+p/EFN9rP9CIqho6MkPWFGcTTF8cAuFH+mXwO8r45GJW0NfBF4pe3f1tFmxThqPR+2N9jem+Ju8UcDD+n3tmnG0C8OSXsCbwR2Bx4FbAe8flrtSzoYWGt7dffmfqEOf3AyXFODSuUP6mD76vLrWuAUiv9Ys3KtpB0Byq9r6w7A9rXlf/A54GPUcD4kbU6RWE+2/aVyc+3nol8cszgfZbvXA2dTXB+4p6TObLla/690xXFQOaRl27cAn2C652I/4OmSrqAYzj2Aomc/gXNRYagmwzVja0T5A0lbSbp753vgyTDTYmmnAc8vv38+cGrdAXQSa+kZTPl8lOOsHwcusX1c10u1nov54qjzfEi6t6R7lt/fFXgSxbWBs4Bnl2+r41z0i+PSrg9dUYyFT+1c2H6j7Z1sL6fID2faPpxJnYuW9OQbO0++QeUPdgBOKX5n2Qz4tO1v1NGwpM8ATwSWSboKeCtwLPB5SUcAvwSeM4MYnlhOjTNwBfCiacZA0WN7HrCmHAMGeBM1n4sF4jisxvOxI7CynH22BPi87a9I+gnwWUnHAOdTfBhN03xxnCnp3hQDHhcAL55yHP28nnHPRYuqUOaO14iIHtr2LubAPxr8xi/+Mne8RkQsOotoOGaQJPmIiH5aMlyTJB8R0U968hERLdWiC69J8hER/bSkJ9/kefIxQ5KeIcmSdu/atlxd1ShrjOWJnQqDNbR1qKS3lN+fJOnZg/aZ5zj3llTLVNuYBqesQbTeYcB3KW4yWdTKudxVvQ74yLht2v5f4BpJ+417rJiBFCiLNitrs+xHUaK1b5Iv64l/QkWd/fMl7V9uf4GkL0n6RlnP+z1d+xwh6b8lnS3pY5I+1Oe4f9pVi/z8zt3GwNaS/l3SpZJOLu+oRNKB5fvWlMXB7lJuv0LSWyR9F3iOpF3KmFZL+s/uv1C62t4NuMX2uj6vHV327JeUx36XpB9IWiXpEZK+KelySd03//wHcHi1sx6Nkzteo8UOBb5h+78lXSfpEZ0Su12OArD9sDJhfqtMklAU6tqHoub4TyV9ENgAvBl4BHAjcCbw4z5tvwY4yvb3yg+bm8vt+wAPpahD8j1gPxULuJwEHFjG+kngJRT1SwButv14AElnAC+2fZmkx1D01g/oaXs/oPfnpPygugfwQtsuP1+utP04Se8vY9gP2BK4GPhouesq4Jg+P2MsBoukpz5IevLRz2EUBZ8ovx7W5z2PB/4NwPalwC+ATpI/w/YNtm8GfgLcn6JQ1XdsX2f7NuAL87T9PeA4SS8H7tlVMvY821eVRcAuAJYDDwZ+bvu/y/esBLorhH4Obv/L5I+BL5TlCE6guC2/147A//Zse3MZx4t859vDO3WU1gDn2r6xHKK5uVPThaJg2n3m+Tmjyar04tOTj8VI0r0oerh7SjJF3SBLel3vWxc4zC1d32+g+D1b6P23s32spK8CTwX+S9KTxjjm78uvSyhqjA9aiu4mih57tx8Cj5S0ne3uxVM68cz1xDbHHf+vtiyPGYtRevLRUs8GPmn7/raX294Z+DlFz73bOZTjzeUwzf2Any5w3POAP5W0rYoysM/q9yZJu9heY/vdFMMdG42dd7kUWC7pQeXz5wHf6X1TWfP955KeU7YhSXv1Od4lwIN6tn2DohDaV7uuD1S1G7OtWBqjEpldE611GEXN/G5fBP66Z9tHgKWS1lAMi7ygrCHel+1fAe+iWEnp2xTDODf0eesrJV0k6ccUveCvL3DMm4EXUgzDrKHoRX90nrcfDhxRHvdi+i8leQ6wT+eiblc7X6CoE3+airK6Ve0PfHWI90eTtGS4JlUoozaStrb9u7InfwpF+ejeD5SZkvTPwJdtf3sCxzoHOMT2b8aPLOqk7Zea5249+I0f/u3AKpSSDgL+mWLo819tHzuRICtKTz7q9LbywudFFENA/zHjePp5F3C3cQ9S1lQ/Lgl+EZtAT768R+PDFEuH7kGx9sAeU478TnLhNWpj+zWzjmEQ29cygRXIypk2TfwQiyomV7vm0cD/2P4ZgKTPUgwV/mQSB68iST4iop9qF1aXlfdrdKywvaLr+X2BK7ueXwU8ZgLRVZYkHxHRq/qF1XUDxuT7TfOt9UJoknxERD+TGa65Cti56/lOFHdt1yYXXiMi+pnMFMofArtKeoCkLShqQY19zWcY6clHRPSa0IVX2+slvQz4JsUUyhNtXzz2gYeQJB8R0c+Ebnay/TXgaxM52AiS5CMiNuJFU7ZgkCT5iIheWeM1IqLlFkltmkFSuyYioke5Pu+yCm9dZ/ugacczjiT5iIgWyzz5iIgWS5KPiGixJPmIiBZLko+IaLEk+YiIFkuSj4hosST5iIgWS5KPiGixJPmIiBb7/1rdStNdtfQfAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 2 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plot_coast(spacing, z)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Let's add another sediment source with a different flux and update the model."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([[ -1. ,  -1. ,  -1. , ...,  -1. ,  -1. ,  -1. ],\n",
       "       [ -1. ,  -1. ,  -1. , ...,  -1. ,  -1. ,  -1. ],\n",
       "       [ -1. ,  -1. ,  -1. , ...,  -1. ,  -1. ,  -1. ],\n",
       "       ..., \n",
       "       [ 22.4,  22.4,  22.4, ...,  22.4,  22.4,  22.4],\n",
       "       [ 22.6,  22.6,  22.6, ...,  22.6,  22.6,  22.6],\n",
       "       [ 22.8,  22.8,  22.8, ...,  22.8,  22.8,  22.8]])"
      ]
     },
     "execution_count": 14,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "qs[0, 150] = 500\n",
    "for time in range(3750):\n",
    "    waves.update()\n",
    "    angle = waves.get_value(\n",
    "        \"sea_surface_water_wave__azimuth_angle_of_opposite_of_phase_velocity\"\n",
    "    )\n",
    "    cem.set_value(\n",
    "        \"sea_surface_water_wave__azimuth_angle_of_opposite_of_phase_velocity\", angle\n",
    "    )\n",
    "\n",
    "    cem.set_value(\"land_surface_water_sediment~bedload__mass_flow_rate\", qs)\n",
    "    cem.update()\n",
    "\n",
    "cem.get_value(\"sea_water__depth\", out=z)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXkAAADuCAYAAADcF3dyAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJzt3Xm0JFWV7/Hvr4pRkMkCREELGUWU0RHbZlo2rcigYjdN89DmCdrwFGdktSOoyFOcRQpBihYHFBEURZFB2gmokqGAAhFEKKhHWTIrUFTd/f6Ik5CVlTfz5BQZN+r3WStX3ozMiNgZN+/Oc0+c2EcRgZmZ1dO0cQdgZmaj4yRvZlZjTvJmZjXmJG9mVmNO8mZmNeYkb2ZWY07yZmY15iRvZlZjTvJmZjW2yrgDMDOrmn322ScWL17c9XVz5879WUTsU0JIfXOSNzNrsXjxYq68+uqur1tl2rQZJYQzECd5M7MWASybqEddLyd5M7MWEbCsJsUbneTNzFYQbsmbmdXVRMCSZRPjDmMonOTNzNpwS97MrKYC98mbmdVYOMmbmdVVhLtrzMxqy+PkzcxqLDy6xsysztwnb2ZWW+6TNzOrMQ+hNDOrObfkzcxqKiJYMuETr2ZmteQhlGZmNeY+eTOzOvPoGjOz+gqPkzczq7e6tOSnjTsAM7OqaUwa0u3WjaTNJF0mab6kGyW9My3fQNLFkm5N9+uP6r04yZuZtbFsIrreMiwF3hMRzwdeBhwlaTvgWOCSiNgKuCQ9HgkneTOzFo0++W63rtuJWBgRv08/PwzMB54N7A/MTi+bDRwworfiPnkzsxXkj66ZIWlO0+NZETGr3QslzQR2Aq4ENo6IhVB8EUjaaLCAJ+ckb2bWoodx8osjYtduL5K0NnAucExEPCRpwAjzOcmbmbUY5sVQklalSPBnR8QP0uJ7JW2SWvGbAIuGsrM23CdvZtYiIoY1ukbA6cD8iDi56akLgMPSz4cB5w/9TSRuyZuZtRhi7ZrdgEOBeZKuTcuOA04EzpF0OHAncNAwdtaOk7yZWRvD6K6JiF8Bk3XA7zXwDjI4yZuZtfDMUGZmNeckb2ZWUxPpxGsdOMmbmbXhKpRmZjXlSUPMzOrMJ17NzOqrTpOGjOyK1yrUUTYz60fjYqghlBoeu1GWNRh7HWUzs37EkCYNqYKRJfkq1FE2M+vXMOrJV0EpffLjqqNsZtaPYOp0x3Qz8iTfbx1lSUcARwCstdZau2y77bajC9LMamPu3LmLI2LDQbYR4SGUWQapo5xmV5kFsMuuu8ZvrrxqlKGaWU2svsr0Pw9jO3VpyY9ydM3Y6yibmfWjcTGU++Q7G1od5Yjg8aVT40y2mU194do13Q2zjvJEGs5kZlaGIU4aMnZT4orXAJZMOMmbWXmmSndMN1MiyZuZlcmThpQsInjC3TVmViK35EsUgU+8mllpJiJYUpOcMyWS/EQEjy5dNu4wzGwl4pa8mVlNeXRNyQIPoTSz8risQcnqNKmumU0FU+eK1m6mRJIP8OgaMyuNu2vMzGosanSV/ZRI8nWqI2Fm1dcoUFYHUyTJ1+db1cymgPCkIaUKgseXeZy8mZXDffJmZjXn7poSTQQ8WpNLjM2s+lygrGQRHkJpZuUJojblzadEkp/Ao2vMrDzukzczqzN315TLpYbNrEweJ18yXwxlZmVyd03JPMermZWpTg3LaeMOwMysipZFdL3lkHSGpEWSbmha9lFJd0u6Nt1e02H9aZJ2kvRaSXtK2riX99G1JS/p5cC/A/8AbAI8CtwAXAh8MyIe7GWH/ZiI4LEnlo56N2ZmwNDHyZ8JfBk4q2X55yLiM5OtJGkL4APA3sCtwF+ANYCtJf0dOBWYHREd/+XomOQl/RS4Bzgf+ASwqLETYA/gfEknR8QFnbYzqKJ2TT36x8ys+oZ54jUirpA0s49VTwBOAY6MWD4YSRsB/wYcCszutJFuLflDI2Jxy7JHgN+n22clzegl6n54ZigzK1d2d8wMSXOaHs+KiFmZOzla0v8C5gDviYj7l4sg4uBJo4tYBHw+Zycdk3xrgpe0TvM6EXFfmy+BxmvPAPYFFkXE9mnZR4G3UvzbAXBcRPwkJ1Azs7L0MLpmcUTs2scuTgGOT7s6Hvgs8B/tXihpOvBaYCbL59+Tc3aUNbpG0pHAxyn64xvvPIDndVjtTProh2rH0/+ZWakClo0w50TEvY2fJZ0G/LjDy38EPAbMA3oOKncI5XuBF0zWam9ngH6oFbeFu2vMrGQjvBhK0iYRsTA9PJBiMMtkNo2IF/W7r9whlLcBf+93Jy2OlnR9Gla0/mQvknSEpDmS5jz24P2TvczMbPgiivK33W4ZJH0b+C2wjaQFkg4HTpI0T9L1FINY3tVhEz+V9Op+30puS/6DwG8kXQk83lgYEe/ocX/Z/VDp5MUsgBlbPj/ckjezUg1vdE27E6in97CJ3wHnSZoGPAGo2Gysk7NybpI/FbiUPvuEGnrsh2paz901Zlay6pQ1+CzwcmBe61DKHLlJfmlEvLvXjbfqsR/qSRMEjy719H9mVpKgSkn+VuCGfhI85Cf5yyQdQXGWt7m75r7JVkj9ULtTjCNdAHwE2F3SjhSH8A7gyH6CNjMbrYDq9B4sBC5PF6c259/hDaGkuLIKir75J/dBhyGUQ+iHetKESw2bWdmqU2r4T+m2Wrr1JDfJP6/NZbVr9LqzfkW4CqWZlSioTJKPiI8Nsn5ukj+dplEwktYCLgD2GmTnuYLwHK9mVq4x98lLmgV8KSLmtXluLeBfgMcj4uxO28lN8ndLOiUi3p7Gtl8InNZr0GZmU0NUoSX/VeBDkl5IMUilUYVyK2Ad4AygY4KHzCQfER+S9GlJXwN2AU6MiHP7jbxXHkJpZqUqiteMN4SIa4E3SVob2JWnSr3Pj4hbcrfTrdTw65seXgV8KN2HpNdHxA96jrwPEeETr2ZWrooMoYyIR4DL+12/W0v+dS2PrwFWTcsDKCXJT+CWvJmVqFrj5AfSrdTwW8oKxMysOirRJz8U3bpr/gv4Smsx+6bn9wSeFhFZ5Qn6FeErXs2sZCtDS56iVs2PJT1GMRNU89ndHYFfAJ8caYQUx3rJMid5MytRRVrykrYG3gc8l+UnDdkzZ/1u3TXnU8zjuhWwG8XZ3YeAbwJHRMSjfcbdk6KefDUOuJmtBCowuqbJ94CvUQxb77m1mzuE8laKIjlmZiuB/HrxJVgaEaf0u3LuxVBjFZ7+z8zKVIGyBpI2SD/+SNJ/AueRWSCy2RRJ8h5CaWYlG39Lfi7F143S4/c1Pddtju0nTY0kj1vyZlayMbfkI2JzKIpBRsRjzc/1UiAya45XSVtLukTSDenxi9LwSjOz+olUT77brRy/yVzWVm5L/jSKfxVOBYiI6yV9Czghd0eDiIAnPE7ezMo0/j75ZwLPBtaUtBNPddusAzwtdzu5Sf5pEXGVpOZlS3N3MrDGt6qZWVnG3yf/T8CbgU2B5lmgHgaOy91IbpJfLGkLis5+JL2RYkqq8jjJm1lZKjC6JiJmA7MlvWGQqr+5Sf4oYBawraS7KaaiOqTfnZqZVd74W/INl0v6IvBKiq+fXwEfj4i/5qzcNclLmgbsGhF7p9lIpkXEw4NE3LMAXGrYzEpTqQJl3wGuAN6QHh8CfBfYO2flrkk+IiYkHQ2cExF/6zfKgbhP3szKVK2yBhtExPFNj0+QdEDuyrndNRdLei/Ft8eTiT73iquh8ETe1uqK+fCq5487Cqur6nTXXCbpX4Fz0uM3UkzBmiU3yTcm8T6qaVn2FVdmZlNKBU68NjkSeDdFYcgApgN/k/RuICJinU4r5xYo23zQKAfi7hqbzBXzi/vdthlvHFYz1SlQFhFPH2T9rCQvaVXg7cCr0qLLgVMj4okO65wB7Assiojt07INKLp8ZgJ3AG+abEKS5QTwhC+GWqld9Ud4yZaTP+/Phw1bRZK8iguUDgE2j4jjJW0GbBIRV+Wsn9tdcwrF3K5fTY8PTcv+d4d1zgS+DJzVtOxY4JKIOFHSsenxB7rv3i35ldbc25/6udNnoPFc4/W7uCfRBlCt7pqvUkx1vSdwPPAI8BXgxTkr5yb5F0fEDk2PL5V0XacVIuIKSTNbFu8P7J5+nk3xH0FGkjczK1l1GpYvjYidJV0DEBH3S1otd+XcJL9M0hYRcRuApOfRxwwlwMYRsRAgIhZK2ihrrWoNZ7JxyWnJN8y9HXacOdJwrM4qNU7+CUnTeariwIYULfssuUn+fRTDeG6nKJLzXOAtPQbaE0lHAEcAsNb64DleVz7z7lr+cetnYL+T4IL3t39usmVmOYLK9MkDX6SYMGQjSZ+gGEKZXQU4d3TNJWme120okvzNEfF4l9XauVfSJqkVvwmwqMM+Z1GUUkDPeE7gOV6t+TOw5TuK+/1OKu7/+MXifrtNi/ubFhRfEo3HZr2qSEs+Is6WNBfYiyL/HhAR83PX72XSkF0oRsWsAuwgiYg4q/MqK7gAOAw4Md2f3+P6ZmblqEhLXtILgW0pGsXze0nwkD+E8r+BLYBreaovPlh+5EzrOt+mOMk6Q9IC4CMUyf0cSYcDdwIHZUXpPvmVzy33rLhs2QRsc0z71zc+H63r+XNj/ajAtTmS1qVoCG8GXE/Rin+hpDuB/SPioZzt5LbkdwW2i8j//yUiDp7kqb1yt9G0NfCkIdbpM3DLPbDFxr2tY9bJ+FvyxwNzgD0jYgIgnYD9FPAJ4P/kbCQ3yd8APJOya8g3uCVv0PkzsN9JcP3/XXH5bfcW9zM3HE1MVk9DHCc/wIWhewMvaiR4gIhYJuk4YF7u/jvO8SrpR5IuAGYAN0n6maQLGrfcnZiZTTkT0f2W50xgn5ZljQtDtwIuSY9bLYmIFWbgS8uyB750a8l/JndDI1WB/jGrgDv+Ai/q8Hwv4+jNuhlSS36AC0PXaJnbtUHA6rn775jkI+KXAGmykEdTbfmtKc70/jR3J0PhP1JrDJeczF0dJsrx58d6Edkt9RmS5jQ9npWGf3eTc2HoQpaf27XZ/8sJDvL75K8A/kHS+hT/WswB/oWypgB0S37lce+Dkz93wfs7J/rmi6Na3XM/bLzuYLFVxb0PFu+lcawa76v52NXlvY5T3hwWiyNi11HsPiL2GMZ2OvbJN1FE/B14PfCliDgQeMEwAjAzq6Th9cm3c2+6IJRuF4YOKrclL0kvp2i5H56WTR9NSG14jleD7t013Uz1z9Bfm6ZWbm61t3tfU/29jtvoyxqUdmFobpJ/J/BB4LyIuDEVKLtsVEG15e4ay9H4ImjXbfPXh2G9tcqNZ1ge6DC9cru/jWUTy68zVd/32AyvQNkgF4amWvKbRsRd7Z7PkVu75gqKfvnG49uBd/S7056FL4Zaaay9RnH/yGOj2X4j8TX2Uwft/jZal/nvp3fDG13T94WhERGSfkhRVqYvuX3yZmYrj0Z3zej65HvxO0lZE4S000uBsvHxFa8rj0eXDGc7nUbawNT6PE12TFrPUTS/39b398hjsGb2PBMGVfqM7AG8TdIdwN8oxslHRHS6auRJUyPJE73XBl+SXr/a9OLn1aY/tXy18s4Z2xh1SvSPLqn+52BJh898u5PQzecj2v29NL4sqv6+YcX33vr32/z8KN5P/jj5MvzzICvnVqE8CTgBeBS4CNgBOCYivjnIzrMFy/9Sp/fQy9RYr3n91m21+8butI/G6ydbt5vGtpdN9LaNXt63ddfvZ6pf3X7PzZ+LQex3ElzY7ir5ZMmywd7vZPH1Ev9k+59s3cn+fts9Hpbq1JP/s6RXAltFxDfSzFBr566f+5t+dSpruS+wANiaYrYoM7N6iuh+K4Gkj1CUPPhgWrQqkN3Azu2uWTXdvwb4dkTcV4zsKUnrwRxmX9lk28rZR79xNK/Xyzaq00dYvn7HyHcaUtmsCsc2N4acY/HaE8s/J1G3z3J1umsOBHYCfg8QEfdIenruyrkt+R9Jupmirvwl6d+FEY1xMxuBQS+kspVLo5RKt1s5lqS5PBoTefd00UNWko+IY4GXA7tGxBMUZ3j37zFQs97td9LwEvTKluiHeexWRhXprqG4aOpUYD1JbwV+AXw9d+WsJC/pIGBpKlj/XxT9Qc/qJ1ozs8qr0Dj5iPgM8H3gXGAb4MMR8cXc9XP75D8UEd9LZ3j/iaLO/CnAS3uMt3ztWjLd+metGkbRCm3e5lT7HAxyXqKu73WU76sio2skfToiPgBc3GZZ9/Vzpm2VdE1E7CTpU8C8iPhWY1nfkfdAT3tmsO1hva3Uyx9Evx+UYXyB5MQ51f5Ah6HsboaqH+NhHo9hft6Huf1e9zPZvq85ae6g5X+16obB+hk90n85feB9dY1F+n1E7Nyy7Prci6Fyk/yPgbsp5hzchWK8/FURsUPvIfdOz1LwzB4+OIN+SIa57XbbK/sPdhgxl6FK/cdVSPqjPh69vMd+YxlFo6ebjzF4kl9lRrBeRpL/6xkjS/KS3g78J/A84Lamp54O/Doi/j1nO7ndNW+imKPwMxHxQKp/7HHyZlZf4x/m+S2KGfg+xfJzwD4cEfflbiSrJQ8gaQfgH9LD/4mI63J3Mig9S8ERTQtG3Tpu1thXlVqZVTLqf/2tPJP9Lsv4z3OY+xhWS37tfbu/8MHZI++uaUhTBD5ZPjUi7sxaL7O75p3AW4EfpEUHUsxl+KXeQ+3dCknezGwyw0jy02cEa7+2+wsfOquMPvnXUcz1+iyKGaSeC8yPiKzZ+XK7aw4HXhoRf0s7/TTwW6CUJG9mVq5KFSg7AXgZ8Is0AGYPYLIa9SvInv4PaK4CtCwt60sqmflw2s7Ssv7dMTPLVp0k/0RE/FXSNEnTIuKy1NDOkpvkvwFcKem89PgA4PReI22xR0QsHnAbZmbDF1RmnDzwgKS1KWbnO1vSImBp7sq50/+dLOly4JUULfi3RMQ1fQRrZjY1jH90TcP+FLXC3gUcAqwLfDx35a5JXtI04PqI2J5UBW0IAvi5pABOjYhZbfZ7BKTTresOaa9mZjkqMGmIpGOAXwPXRESju3x2r9vpmuQjYkLSdZKekztkJ8NuqVzmRsDFkm5Ok4U373cWMAvS6BozszKNv7tmU+ALwLaSrgd+Q5H0f9vLOPncPvlNgBslXUVRgRKAiNgvP96nRMQ96X5R6ud/CUV/k5lZNYy5JR8R7wWQtBpFmfdXAP8BnCbpgYjYLmc7uUn+Y31F2UaqhTwtIh5OP7+aHvqXzMxKMf6WfMOawDoUHdfrAvcA83JX7pjkJW0JbBwRv2xZ/iqKWjb92Bg4L80stQrwrYi4qM9tmZmNQMD0nNnhRheBpFnACyiGm19J0V1zckTc38t2urXkPw8c12b539Nzr+tlZwARcTvFROBmZtUkYFpGS36ESR54DrA6cCtFo3oB8ECvG+mW5GdGxPWtCyNijqSZve7MzGzKGPN4j4jYR0WXxwso+uPfA2wv6T6Kk68fydlOtyS/Rofn1syK1MxsKsppyY9Ymtv1BkkPAA+m274Ug1Wykny36f+uTnMKLkfS4cDc3sI1M5sqomjJd7uNkKR3SPqOpLsoRh/uC9wCvB7YIHc73Vryx1CcJD2Ep5L6rsBqFJUozczqR4y9uwaYSTG367siYmG/G+mY5CPiXuAVqerZ9mnxhRFxab87NDObElYZb1mDiHj3MLaTW7vmMuCyYezQzKzyqtGSH4rci6HMzFYiUYkTr8PgJG9m1o5b8mZmNZV7MdQU4CRvZtbOtMrUkx+Ik7yZ2QrcJ29mVl9DHF0z7jmtneTNzNoZbkt+bHNaO8mbmbVTk+6abrVrzMxWPo3umu61a2ZImtN0O6LN1hpzWs+d5PmRckvezGwFmZOGwOKMPvauc1qPklvyZmatGuPku90yNM9pDTTmtC6Nk7yZWTtDKDUsaS1JT2/8TDGn9Q0jjnw57q4xM2tnOCdexz6ntZO8mdkKhjMpSBXmtHaSNzNr5do1ZmY1lze6pvKc5M3MWnnSEDOzOqtPgbKxDKGUtI+kWyT9UdKx44jBzKyjIQyhrILSk7yk6cBXgH8GtgMOlrRd2XGYmXVUkyQ/ju6alwB/TEOLkPQdYH/gpjHEYma2IgGr1OPE6zi6a54N3NX0eEFaZmZWERmteLfkJ6U2y1Y4WqlaW6Ni2+N8rNxLgduYAYylHnSLKsRRhRigGnFUIQaoRhxViAFgm4G34HHyA1kAbNb0eFPgntYXRcQsYBaApDllz6bSqgoxVCWOKsRQlTiqEENV4qhCDI04hrOheiT5cXTXXA1sJWlzSasB/wpcMIY4zMwmN6QqlONWeks+IpZKOhr4GTAdOCMibiw7DjOzSfliqMFExE+An/SwyqxRxdKDKsQA1YijCjFANeKoQgxQjTiqEAMMJY7sSUMqTxH1+LYyMxsWrbNmsOvzur/wspvmVuE8RCcua2Bm1k5NumsqPTNUVcofSLpD0jxJ1w7tzH3efs+QtEjSDU3LNpB0saRb0/36Y4jho5LuTsfjWkmvGXEMm0m6TNJ8STdKemdaXvaxmCyO0o6HpDUkXSXpuhTDx9LyzSVdmY7Fd9OghpHpEMeZkv7UdCx2HGUcaZ/TJV0j6cfp8XCORU1OvFY2yVew/MEeEbFjyf+anQns07LsWOCSiNgKuCQ9LjsGgM+l47FjOscySkuB90TE84GXAUelz0LZx2KyOKC84/E4sGdE7ADsCOwj6WXAp1MMWwH3A4ePMIZOcQC8r+lYXDviOADeCcxvejyEY1Gfi6Eqm+RpKn8QEUuARvmDlUaa0f2+lsX7A7PTz7OBA8YQQ6kiYmFE/D79/DDFH/SzKf9YTBZHaaLwSHq4aroFsCfw/bS8jGMxWRylkrQp8Frg6+mxGMaxGOJE3uNW5SRfpfIHAfxc0tx0Je44bRwRC6FIOsBGY4rjaEnXp+6ckXaTNJM0E9gJuJIxHouWOKDE45G6J64FFgEXA7cBD0TE0vSSUv5WWuOIiMax+EQ6Fp+TtPqIw/g88H6gMRTmGQzrWEyf6H6bAqqc5LPKH5Rkt4jYmaLr6ChJrxpTHFVxCrAFxb/pC4HPlrFTSWsD5wLHRMRDZewzM45Sj0dELIuIHSmuFn8J8Px2LxtlDO3ikLQ98EFgW+DFwAbAB0a1f0n7AosiYm7z4nah9r5x3F1TgqzyB2WIiHvS/SLgPIo/rHG5V9ImAOl+UdkBRMS96Q98AjiNEo6HpFUpEuvZEfGDtLj0Y9EujnEcj7TfB4DLKc4PrCepMVqu1L+Vpjj2SV1aERGPA99gtMdiN2A/SXdQdOfuSdGyH8KxyOiqcXfNwCpR/kDSWpKe3vgZeDWMtVjaBcBh6efDgPPLDqCRWJMDGfHxSP2spwPzI+LkpqdKPRaTxVHm8ZC0oaT10s9rAntTnBu4DHhjelkZx6JdHDc3femKoi98ZMciIj4YEZtGxEyK/HBpRBzCsI5FTVrylR0nX6HyBxsD5xWfWVYBvhURF5WxY0nfBnYHZkhaAHwEOBE4R9LhwJ3AQWOIYfc0NC6AO4AjRxkDRYvtUGBe6gMGOI6Sj0WHOA4u8XhsAsxOo8+mAedExI8l3QR8R9IJwDUUX0ajNFkcl0rakKLD41rgbSOOo50PMOixqFEVSl/xambWQuuvHuz1zO4vPPdOX/FqZjblTKHumG6c5M3M2qlJd42TvJlZO27Jm5nVVI1OvDrJm5m1U5OWfJXHydsYSTpQUkjatmnZTDVVoywxlt0bFQZL2NcBkj6cfj5T0hu7rTPJdjaUVMpQWxuFcFkDq72DgV9RXGQypaWx3LneD3x10H1GxF+AhZJ2G3RbNgYuUGZ1lmqz7EZRorVtkk/1xL+hos7+NZL2SMvfLOkHki5K9bxPalrncEl/kHS5pNMkfbnNdv+xqRb5NY2rjYG1JX1f0s2Szk5XVCJpr/S6eak42Opp+R2SPizpV8BBkrZIMc2V9D/N/6E07Xtr4PGIWNzmueNTy35a2vYnJf1W0hxJO0v6maTbJDVf/PND4JC8o26V4ytercYOAC6KiD9Iuk/Szo0Su02OAoiIF6aE+fOUJKEo1LUTRc3xWyR9CVgGfAjYGXgYuBS4rs2+3wscFRG/Tl82j6XlOwEvoKhD8mtgNxUTuJwJ7JViPQt4O0X9EoDHIuKVAJIuAd4WEbdKeilFa33Pln3vBrS+T9IX1brAWyIi0vfLXRHxckmfSzHsBqwB3Ah8La06BzihzXu0qWCKtNS7cUve2jmYouAT6f7gNq95JfDfABFxM/BnoJHkL4mIByPiMeAm4LkUhap+GRH3RcQTwPcm2fevgZMlvQNYr6lk7FURsSAVAbsWmAlsA/wpIv6QXjMbaK4Q+l148j+TVwDfS+UITqW4LL/VJsBfWpZ9KMVxZCx/eXijjtI84MqIeDh10TzWqOlCUTDtWZO8T6uynFa8W/I2FUl6BkULd3tJQVE3KCS9v/WlHTbzeNPPyyg+Z51e/6SIOFHShcBrgN9J2nuAbf4t3U+jqDHebSq6Ryla7M2uBnaRtEFENE+e0ohnoiW2CZ76u1ojbdOmIrfkrabeCJwVEc+NiJkRsRnwJ4qWe7MrSP3NqZvmOcAtHbZ7FfCPktZXUQb2De1eJGmLiJgXEZ+m6O5Yoe+8yc3ATElbpseHAr9sfVGq+f4nSQelfUjSDm22Nx/YsmXZRRSF0C5sOj+Qa2vGW7HU+iU8usZq62CKmvnNzgX+rWXZV4HpkuZRdIu8OdUQbysi7gY+STGT0i8ounEebPPSYyTdIOk6ilbwTzts8zHgLRTdMPMoWtFfm+TlhwCHp+3eSPupJK8Admqc1G3az/co6sRfoKKsbq49gAt7eL1VSU26a1yF0kojae2IeCS15M+jKB/d+oUyVpK+APwoIn4xhG1dAewfEfcPHpmVSRtND960dvcXfuWhrlUoJe0DfIGi6/PrEXHiUILM5Ja8lemj6cTnDRRdQD8cczztfBJ42qAbSTXVT3aCn8KG0JJP12h8hWLq0O0o5h7YbsSRL8cnXq00EfHeccfQTUTcyxBmIEsjbar4JWY5hle75iXAHyPidgBJ36HoKrxpGBvP4SRvZtbU8BKQAAAA9UlEQVRO3onVGel6jYZZETGr6fGzgbuaHi8AXjqE6LI5yZuZtco/sbq4S598u2G+pZ4IdZI3M2tnON01C4DNmh5vSnHVdml84tXMrJ3hDKG8GthK0uaSVqOoBTXwOZ9euCVvZtZqSCdeI2KppKOBn1EMoTwjIm4ceMM9cJI3M2tnSBc7RcRPgJ8MZWN9cJI3M1tBTJmyBd04yZuZtfIcr2ZmNTdFatN049o1ZmYt0vy8MzJeujgi9hl1PINwkjczqzGPkzczqzEneTOzGnOSNzOrMSd5M7Mac5I3M6sxJ3kzsxpzkjczqzEneTOzGnOSNzOrsf8PCS3ZBWUbl9UAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 2 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plot_coast(spacing, z)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Here we shut off the sediment supply completely."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([[ -1. ,  -1. ,  -1. , ...,  -1. ,  -1. ,  -1. ],\n",
       "       [ -1. ,  -1. ,  -1. , ...,  -1. ,  -1. ,  -1. ],\n",
       "       [ -1. ,  -1. ,  -1. , ...,  -1. ,  -1. ,  -1. ],\n",
       "       ..., \n",
       "       [ 22.4,  22.4,  22.4, ...,  22.4,  22.4,  22.4],\n",
       "       [ 22.6,  22.6,  22.6, ...,  22.6,  22.6,  22.6],\n",
       "       [ 22.8,  22.8,  22.8, ...,  22.8,  22.8,  22.8]])"
      ]
     },
     "execution_count": 16,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "qs.fill(0.0)\n",
    "for time in range(4000):\n",
    "    waves.update()\n",
    "    angle = waves.get_value(\n",
    "        \"sea_surface_water_wave__azimuth_angle_of_opposite_of_phase_velocity\"\n",
    "    )\n",
    "    cem.set_value(\n",
    "        \"sea_surface_water_wave__azimuth_angle_of_opposite_of_phase_velocity\", angle\n",
    "    )\n",
    "\n",
    "    cem.set_value(\"land_surface_water_sediment~bedload__mass_flow_rate\", qs)\n",
    "    cem.update()\n",
    "\n",
    "cem.get_value(\"sea_water__depth\", out=z)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXkAAADuCAYAAADcF3dyAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJzt3XnUHFWd//H3J2EVZDOAQdAgsogoW1yjDttxMooCKs4wDIMOP0GFn+KOnMENVOAo7iJBkDDigiISRVFkkREVSCQQICCCCIH8iJEdCSHP8/39UdWk0+mn+/ZWXV35vM7p83RV13Krnu5v375163sVEZiZWTVNGnYBzMxscBzkzcwqzEHezKzCHOTNzCrMQd7MrMIc5M3MKsxB3syswhzkzcwqzEHezKzC1hp2AczMymbmzJmxdOnStsvNmzfvlxExs4Aidc1B3syswdKlS7nmuuvaLrfWpElTCihOTxzkzcwaBDA2Xo28Xg7yZmYNImCsIskbHeTNzFYTrsmbmVXVeMDysfFhF6MvHOTNzJpwTd7MrKICt8mbmVVYOMibmVVVhJtrzMwqy/3kzcwqLNy7xsysytwmb2ZWWW6TNzOrMHehNDOrONfkzcwqKiJYPu4Lr2ZmleQulGZmFeY2eTOzKnPvGjOz6gr3kzczq7aq1OQnDbsAZmZlUxs0pN2jHUnbSLpC0kJJN0t6Xz5/M0mXSro9/7vpoI7FQd7MrImx8Wj7SLAC+GBEvBB4BXC0pJ2B44DLImJ74LJ8eiAc5M3MGtTa5Ns92m4nYnFE/DF//iiwEHgOcAAwO19sNnDggA7FbfJmZqtJ710zRdLcuulZETGr2YKSpgG7A9cAW0bEYsi+CCRt0VuBJ+Ygb2bWoIN+8ksjYnq7hSRtCFwAHBsRj0jqsYTpHOTNzBr082YoSWuTBfjzIuLH+ez7JU3Na/FTgSV92VkTbpM3M2sQEf3qXSPgLGBhRJxW99Ic4PD8+eHARX0/iJxr8mZmDfqYu2YGcBiwQNL8fN7xwMnA+ZKOAO4GDu7HzppxkDcza6IfzTUR8Vtgogb4fXveQQIHeTOzBh4Zysys4hzkzcwqajy/8FoFDvJmZk04C6WZWUV50BAzsyrzhVczs+qq0qAhA7vjtQx5lM3MulG7GaoPqYaHbpBpDYaeR9nMrBvRp0FDymBgQb4MeZTNzLrVj3zyZVBIm/yw8iibmXUjGJ3mmHYGHuS7zaMs6UjgSIANNthgz5122mlwhTSzypg3b97SiNi8l21EuAtlkl7yKOejq8wC2HP69PjdNdcOsqhmVhHrrjX5r/3YTlVq8oPsXTP0PMpmZt2o3QzlNvnW+pZHOSJ4csVoXMk2s9EXzl3TXj/zKI/n3ZnMzIrQx0FDhm4k7ngNYPm4g7yZFWdUmmPaGYkgb2ZWJA8aUrCI4Ck315hZgVyTL1AEvvBqZoUZj2B5RWLOSAT58QieWDE27GKY2RrENXkzs4py75qCBe5CaWbFcVqDglVpUF0zGwWjc0drOyMR5APcu8bMCuPmGjOzCosK3WU/EkG+SnkkzKz8agnKqmBEgnx1vlXNbASEBw0pVBA8OeZ+8mZWDLfJm5lVnJtrCjQe8ERFbjE2s/JzgrKCRbgLpZkVJ4jKpDcfiSA/jnvXmFlx3CZvZlZlbq4pllMNm1mR3E++YL4ZysyK5OaagnmMVzMrUpUqlpOGXQAzszIai2j7SCHpbElLJN1UN++Tku6VND9/vL7F+pMk7S7pDZL2kbRlJ8fRtiYv6ZXAfwCvAaYCTwA3ARcD34mIhzvZYTfGI1j21IpB78bMDOh7P/lzgK8B5zbM/2JEfH6ilSRtB3wU2A+4HfgbsB6wg6R/AGcAsyOi5U+OlkFe0i+A+4CLgM8AS2o7AfYGLpJ0WkTMabWdXmW5a6rRPmZm5dfPC68RcZWkaV2sehJwOnBUxKqFkbQF8O/AYcDsVhtpV5M/LCKWNsx7DPhj/viCpCmdlLobHhnKzIqV3BwzRdLcuulZETErcSfHSPpPYC7wwYh4cJUSRBwyYekilgBfStlJyyDfGOAlbVS/TkQ80ORLoLbs2cD+wJKI2CWf90ngnWQ/OwCOj4ifpxTUzKwoHfSuWRoR07vYxenAifmuTgS+APxXswUlTQbeAExj1fh7WsqOknrXSDoK+DRZe3ztyAN4fovVzqGLdqhmPPyfmRUqYGyAMSci7q89l3Qm8LMWi/8UWAYsADouVGoXyg8BL5qo1t5MD+1Qq28LN9eYWcEGeDOUpKkRsTifPIisM8tEto6Il3S7r9QulHcA/+h2Jw2OkXRj3q1o04kWknSkpLmS5i57+MGJFjMz67+ILP1tu0cCSd8Dfg/sKGmRpCOAUyUtkHQjWSeW97fYxC8kva7bQ0mtyX8M+J2ka4AnazMj4r0d7i+5HSq/eDELYMoLXhiuyZtZofrXu6bZBdSzOtjEH4ALJU0CngKUbTY2Slk5NcifAVxOl21CNR22Q9Wt5+YaMytYedIafAF4JbCgsStlitQgvyIiPtDpxht12A71tHGCJ1Z4+D8zK0hQpiB/O3BTNwEe0oP8FZKOJLvKW99c88BEK+TtUHuR9SNdBHwC2EvSbmSn8C7gqG4KbWY2WAHlaT1YDFyZ35xaH3/714WS7M4qyNrmn94HLbpQ9qEd6mnjTjVsZkUrT6rhv+SPdfJHR1KD/POb3Fa7Xqc761aEs1CaWYGC0gT5iPhUL+unBvmzqOsFI2kDYA6wby87TxWEx3g1s2INuU1e0izgqxGxoMlrGwD/CjwZEee12k5qkL9X0ukR8e68b/vFwJmdFtrMbDREGWry3wBOkPRisk4qtSyU2wMbAWcDLQM8JAb5iDhB0imSvgnsCZwcERd0W/JOuQulmRUqS14z3CJEzAfeJmlDYDorU70vjIjbUrfTLtXwm+smrwVOyP+GpDdHxI87LnkXIsIXXs2sWCXpQhkRjwFXdrt+u5r8GxumrwfWzucHUEiQH8c1eTMrULn6yfekXarhdxRVEDOz8ihFm3xftGuu+W/g643J7Ote3wd4RkQkpSfoVoTveDWzgq0JNXmyXDU/k7SMbCSo+qu7uwG/Bj470BKSnevlYw7yZlagktTkJe0AfBh4HqsOGrJPyvrtmmsuIhvHdXtgBtnV3UeA7wBHRsQTXZa7I1k++XKccDNbA5Sgd02dHwLfJOu23nFtN7UL5e1kSXLMzNYA6fniC7AiIk7vduXUm6GGKjz8n5kVqQRpDSRtlj/9qaT3ABeSmCCy3ogEeXehNLOCDb8mP4/s60b59IfrXms3xvbTRiPI45q8mRVsyDX5iNgWsmSQEbGs/rVOEkQmjfEqaQdJl0m6KZ9+Sd690syseiLPJ9/uUYzfJc5rKrUmfybZT4UzACLiRknfBU5K3VEvIuAp95M3syINv03+2cBzgPUl7c7KZpuNgGekbic1yD8jIq6VVD9vRepOelb7VjUzK8rw2+T/GXg7sDVQPwrUo8DxqRtJDfJLJW1H1tiPpLeSDUlVHAd5MytKCXrXRMRsYLakt/SS9Tc1yB8NzAJ2knQv2VBUh3a7UzOz0ht+Tb7mSklfAV5N9vXzW+DTEfH3lJXbBnlJk4DpEbFfPhrJpIh4tJcSdywApxo2s8KUKkHZ94GrgLfk04cCPwD2S1m5bZCPiHFJxwDnR8Tj3ZayJ26TN7MilSutwWYRcWLd9EmSDkxdObW55lJJHyL79ng60KfecdUXHsjbGl21EF77wmGXwqqqPM01V0j6N+D8fPqtZEOwJkkN8rVBvI+um5d8x5WZ2UgpwYXXOkcBHyBLDBnAZOBxSR8AIiI2arVyaoKybXstZU/cXGM1V+dDW87YMft71cJVp836ojwJyiLimb2snxTkJa0NvBt4bT7rSuCMiHiqxTpnA/sDSyJil3zeZmRNPtOAu4C3TTQgySoCeMo3Q63Rrv3zqtNX37b69MteUFx5rPpKEuSV3aB0KLBtRJwoaRtgakRcm7J+anPN6WRju34jnz4sn/d/WqxzDvA14Ny6eccBl0XEyZKOy6c/2n73rsmvsebdmb7s2PjK5fd0S6L1oFzNNd8gG+p6H+BE4DHg68BLU1ZODfIvjYhd66Yvl3RDqxUi4ipJ0xpmHwDslT+fTfaLICHIm5kVrDwVy5dHxB6SrgeIiAclrZO6cmqQH5O0XUTcASDp+XQxQgmwZUQsBoiIxZK2SFqrXN2ZrCjz7+ps+fpa/7w7Ybdp/SyNrVFK1U/+KUmTWZlxYHOymn2S1CD/YbJuPHeSJcl5HvCODgvaEUlHAkcCsMGm4DFe1wwL7mn9+ou3ab5sbX79PL9nrFtBadrkga+QDRiyhaTPkHWhTM4CnNq75rJ8nNcdyYL8rRHxZJvVmrlf0tS8Fj8VWNJin7PIUimgZz038Biv1XfLotav77w1T78PGpcdi9XnLbgnW8esGyWpyUfEeZLmAfuSxd8DI2Jh6vqdDBqyJ1mvmLWAXSUREee2XmU1c4DDgZPzvxd1uL6ZWTFKUpOX9GJgJ7JK8cJOAjykd6H8H2A7YD4r2+KDVXvONK7zPbKLrFMkLQI+QRbcz5d0BHA3cHBSKd0mb5C9B267r/lrE/0K8PvGulGCe3MkbUxWEd4GuJGsFv9iSXcDB0TEIynbSa3JTwd2jkj//RIRh0zw0r6p26jbGnjQkGq64/70ZScK8K34fWPdGn5N/kRgLrBPRIwD5BdgPwd8Bvi/KRtJDfI3Ac+m6BzyNa7Jr5nedOrq8+Z8ZOLlp20Od/1t1Xm1L5Fpm/evXFZ9fewn38ONofsBL6kFeICIGJN0PLAgdf8tx3iV9FNJc4ApwC2SfilpTu2RuhMzs5EzHu0fac4BZjbMq90Yuj1wWT7daHlErDYCXz4vueNLu5r851M3NFAlaB+zAbinxZgHzWrxjfMba/Wt3iN+/1in+lST7+HG0PUaxnatEbBu6v5bBvmI+A1APljIE3lu+R3IrvT+InUnfeEPafVstSnc1z510YTedOqqgb7Vl8Y9f8/2VxX3PbjyeOqfW39Eck19iqS5ddOz8u7f7aTcGLqYVcd2rff/UgoH6W3yVwGvkbQp2U+LucC/UtQQgK7JV9P9D/e+jVrNvlVbfc19D8KWG/e+z2Fodq7qvyAbvyxH9TjLJG0Mi6URMX0Qu4+IvfuxnZZt8nUUEf8A3gx8NSIOAl7UjwKYmZVS/9rkm7k/vyGUdjeG9iq1Ji9JrySruR+Rz5s8mCI14TFeDbLaequ2+pTafK1G/KyeUnQX4++JQynXjqV++RXjK6dH4VjLZvBpDQq7MTQ1yL8P+BhwYUTcnCcou2JQhWrKzTXVs8kG8NDjE0/XaxXga1IDPWQBcJMN0pYdhonOQzNj46svX/95KfuxllL/EpT1cmNonkt+64hok9RpYqm5a64ia5evTd8JvLfbnXYsfDNUJT22bOXzDddb9X9cC+qpQbumkzb6sr6n6s9LihVj2fmrX6/x2B56PFumrBrfC7XpYZa5f71rur4xNCJC0k/I0sp0JbVN3sxszVFrrhlcm3wn/iApaYCQZjpJUDY8KXe8PrE8+7v+Oiuf16atnOr/V2Pjq/7foPNafL2UXwGPLSvX+6Px+FM1q/k3+7yUtcmz8bjrj6f2vPZ/emJ5cf+z8pyvvYF3SboLeJysn3xExEtSVh6NIE+0zg2+vO61xjdMbXqdyasuu05x140ntHysHOUYptrxtwtw7drj260zUcCv3+8w/xfL2zQddZriodnnpfGzMEydfg7r/0+N75VBHE96P/ki/EsvK6dmoTwVOAl4ArgE2BU4NiK+08vOkwWrfggmT1r5LTs5scWp8UNUm67fVr3adsfGWy/T67d9uw93szJVSVG1pZSa/fKxwZ7jbo+1q7t/WwSoiY6z1/9F/WcmVSfv/0Fuo5ny5JP/q6RXA9tHxLfzkaE2TF0/9R39ujyt5f7AImAHstGizMyqKaL9owCSPkGW8uBj+ay1geQKdmpzzdr539cD34uIB7KePQVpPJn1NYVeax8TrZ+yj6Lb7MrTRlicbpppWm2rXW2+TOe4k2NvPLZ2xzGI4yzTueuH8jTXHATsDvwRICLuk5R880NqkP+ppFvJmmvek/9c6LCfl1kJdNLFctSkXIMou06vPQxKuVKpLM+7UtYG8u7opofUfvLHSToFeCTPZ/w4WRa18qnyh3hN089a/ETbrur7pJt7DIap1f96WF9eJWmTJ7tp6gxgE0nvBP4L+Fbqyklt8pIOBlbkAf6/ydqDtuqmtGZmpVeifvIR8XngR8AFwI7AxyPiK6nrpzbXnBARP8yv8P4zWZ7504GXd1jewWhWC2icV8ZaTdVrk70YZC2+2X6q+D8YlWPr9NoDFHNMJanJSzolIj4KXNpkXvv1U4ZtlXR9ROwu6XPAgoj4bm1e1yXvgJ7x7GCnw5u/2E0w6OQN0s8vi07LWvYP5yAUFdwnUqZz3u9zkXpsE+13EOemH8fYWK7rT53Xa/pfrb15sGlCi/Tfzup5X23LIv0xIvZomHdj6s1QqUH+Z8C9ZGMO7kl2AfbaiNi18yJ3ru9BHlq/YTvZZsobv58f1jIFoX4adnBvZpjnepDno9lxDbqy1I/9pZjzkf4E+bWmBJskBPm/nz2wIC/p3cB7gOcDd9S99Ezg6oj4j5TtpDbXvI1sjMLPR8RDef5j95M3s+oafu+a75KNwPc5Vh0D9tGIeCB1I0k1eQBJuwKvySf/NyJuSN1Jr7SVgmc31BoG9VO21+021m7KWENtp6ga7Ciem5qyNl10ol/v+cbtpRjksX6K/tTkN9y//YIPzx54c01NPkTg02k5I+LupPUSm2veB7wT+HE+6yCysQy/2nlRO6etFBxZxJ7M+miioDfKX27dSBkLoJ/6EeQnTwk2fEP7BR85t4g2+TeSjfW6FdkIUs8DFkZE0uh8qc01RwAvj4jH852eAvweKCTIm42kNS2YT2Qkz0OpEpSdBLwC+HXeAWZvYKIc9atJHv4PqM8CNJbP60qeMvPRfDsrivq5Y2aWrDxB/qmI+LukSZImRcQVeUU7SWqQ/zZwjaQL8+kDgbM6LWmDvSNiaY/bMDPrv6A0/eSBhyRtSDY633mSlgArUldOTWtwmqQrgVeT1eDfERHXd1FYM7PRMPzeNTUHkOUKez9wKLAx8OnUldsGeUmTgBsjYhfyLGh9EMCv8oQ7Z0TErCb7PRLyy60b92mvZmYpSjBoiKRjgauB6yOi1lw+u9PttA3yETEu6QZJz03tspNgRp4ucwvgUkm35oOF1+93FjAL8t41ZmZFGn5zzdbAl4GdJN0I/I4s6P++k37yqW3yU4GbJV1LNsYgABHxpvTyrhQR9+V/l+Tt/C8ja28yMyuHIdfkI+JDAJLWAaYDryLLQHmmpIciYueU7aQG+U91Vcom8lzIkyLi0fz56+igfcnMrBDDr8nXrA9sRNZwvTFwH7AgdeWWQV7SC4AtI+I3DfNfS5bLphtbAhfmI0utBXw3Ii7pcltmZgMQMDnhwuuAhpcFkDQLeBFZd/NryJprTouIBzvZTrua/JeA45vM/0f+2hs72RlARNxJNhC4mVk5CZiUUJMfYJAHngusC9xOVqleBDzU6UbaBflpEXFj48yImCtpWqc7MzMbGUPu7xERM5U1ebyIrD3+g8Aukh4gu/j6iZTttAvy67V4bf2kkpqZjaKUmvyARZZc7CZJDwEP54/9yTqrJAX5dsP/XZePKbgKSUcA8zorrpnZqIisJt/uMUCS3ivp+5LuIet9uD9wG/BmYLPU7bSryR9LdpH0UFYG9enAOmSZKM3MqkcMvbkGmEY2tuv7I2JxtxtpGeQj4n7gVXnWs13y2RdHxOXd7tDMbCSsNdy0BhHxgX5sJzV3zRXAFf3YoZlZ6ZWjJt8XqTdDmZmtQaIUF177wUHezKwZ1+TNzCoq9WaoEeAgb2bWzKTS5JPviYO8mdlq3CZvZlZdfexdM+wxrR3kzcya6W9NfmhjWjvIm5k1U5Hmmna5a8zM1jy15pr2uWumSJpb9ziyydZqY1rPm+D1gXJN3sxsNYmDhsDShDb2tmNaD5Jr8mZmjWr95Ns9EtSPaQ3UxrQujIO8mVkzfUg1LGkDSc+sPScb0/qmAZd8FW6uMTNrpj8XXoc+prWDvJnZavozKEgZxrR2kDcza+TcNWZmFZfWu6b0HOTNzBp50BAzsyqrToKyoXShlDRT0m2S/izpuGGUwcyspT50oSyDwoO8pMnA14F/AXYGDpG0c9HlMDNrqSJBfhjNNS8D/px3LULS94EDgFuGUBYzs9UJWKsaF16H0VzzHOCeuulF+Twzs5JIqMW7Jj8hNZm32tnKs7XVMrY9yaeKvRW4iSnAUPJBNyhDOcpQBihHOcpQBihHOcpQBoAde96C+8n3ZBGwTd301sB9jQtFxCxgFoCkuUWPptKoDGUoSznKUIaylKMMZShLOcpQhlo5+rOhagT5YTTXXAdsL2lbSesA/wbMGUI5zMwm1qcslMNWeE0+IlZIOgb4JTAZODsibi66HGZmE/LNUL2JiJ8DP+9glVmDKksHylAGKEc5ylAGKEc5ylAGKEc5ylAG6Es5kgcNKT1FVOPbysysX7TR+sH057df8Ipb5pXhOkQrTmtgZtZMRZprSj0yVFnSH0i6S9ICSfP7duU+bb9nS1oi6aa6eZtJulTS7fnfTYdQhk9Kujc/H/MlvX7AZdhG0hWSFkq6WdL78vlFn4uJylHY+ZC0nqRrJd2Ql+FT+fxtJV2Tn4sf5J0aBqZFOc6R9Je6c7HbIMuR73OypOsl/Syf7s+5qMiF19IG+RKmP9g7InYr+KfZOcDMhnnHAZdFxPbAZfl00WUA+GJ+PnbLr7EM0grggxHxQuAVwNH5e6HoczFROaC48/EksE9E7ArsBsyU9ArglLwM2wMPAkcMsAytygHw4bpzMX/A5QB4H7CwbroP56I6N0OVNshTl/4gIpYDtfQHa4x8RPcHGmYfAMzOn88GDhxCGQoVEYsj4o/580fJPtDPofhzMVE5ChOZx/LJtfNHAPsAP8rnF3EuJipHoSRtDbwB+FY+LfpxLvo4kPewlTnIlyn9QQC/kjQvvxN3mLaMiMWQBR1giyGV4xhJN+bNOQNtJqknaRqwO3ANQzwXDeWAAs9H3jwxH1gCXArcATwUESvyRQr5rDSWIyJq5+Iz+bn4oqR1B1yMLwEfAWpdYZ5Fv87F5PH2jxFQ5iCflP6gIDMiYg+ypqOjJb12SOUoi9OB7ch+pi8GvlDETiVtCFwAHBsRjxSxz8RyFHo+ImIsInYju1v8ZcALmy02yDI0K4ekXYCPATsBLwU2Az46qP1L2h9YEhHz6mc3K2rnG8fNNQVISn9QhIi4L/+7BLiQ7IM1LPdLmgqQ/11SdAEi4v78Az4OnEkB50PS2mSB9byI+HE+u/Bz0awcwzgf+X4fAq4kuz6wiaRab7lCPyt15ZiZN2lFRDwJfJvBnosZwJsk3UXWnLsPWc2+D+cioanGzTU9K0X6A0kbSHpm7TnwOhhqsrQ5wOH588OBi4ouQC2w5g5iwOcjb2c9C1gYEafVvVTouZioHEWeD0mbS9okf74+sB/ZtYErgLfmixVxLpqV49a6L12RtYUP7FxExMciYuuImEYWHy6PiEPp17moSE2+tP3kS5T+YEvgwuw9y1rAdyPikiJ2LOl7wF7AFEmLgE8AJwPnSzoCuBs4eAhl2CvvGhfAXcBRgywDWY3tMGBB3gYMcDwFn4sW5TikwPMxFZid9z6bBJwfET+TdAvwfUknAdeTfRkN0kTluFzS5mQNHvOBdw24HM18lF7PRYWyUPqOVzOzBtp03WDfZ7df8IK7fcermdnIGaHmmHYc5M3MmqlIc42DvJlZM67Jm5lVVIUuvDrIm5k1U5GafJn7ydsQSTpIUkjaqW7eNNVloyywLHvVMgwWsK8DJX08f36OpLe2W2eC7WwuqZCutjYI4bQGVnmHAL8lu8lkpOV9uVN9BPhGr/uMiL8BiyXN6HVbNgROUGZVludmmUGWorVpkM/ziX9bWZ796yXtnc9/u6QfS7okz+d9at06R0j6k6QrJZ0p6WtNtvtPdbnIr6/dbQxsKOlHkm6VdF5+RyWS9s2XW5AnB1s3n3+XpI9L+i1wsKTt8jLNk/S/9b9Q6va9A/BkRCxt8tqJec1+Ur7tz0r6vaS5kvaQ9EtJd0iqv/nnJ8ChaWfdSsd3vFqFHQhcEhF/kvSApD1qKXbrHA0QES/OA+av8iAJWaKu3clyjt8m6avAGHACsAfwKHA5cEOTfX8IODoirs6/bJbl83cHXkSWh+RqYIayAVzOAfbNy3ou8G6y/CUAyyLi1QCSLgPeFRG3S3o5WW19n4Z9zwAaj5P8i2pj4B0REfn3yz0R8UpJX8zLMANYD7gZ+Ga+6lzgpCbHaKNgRGrq7bgmb80cQpbwifzvIU2WeTXwPwARcSvwV6AW5C+LiIcjYhlwC/A8skRVv4mIByLiKeCHE+z7auA0Se8FNqlLGXttRCzKk4DNB6YBOwJ/iYg/5cvMBuozhP4Anv5l8irgh3k6gjPIbstvNBX4W8O8E/JyHBWr3h5ey6O0ALgmIh7Nm2iW1XK6kCVM22qC47QyS6nFuyZvo0jSs8hquLtICrK8QSHpI42LttjMk3XPx8jeZ62Wf1pEnCzpYuD1wB8k7dfDNh/P/04iyzHebii6J8hq7PWuA/aUtFlE1A+eUivPeEPZxln5uVov36aNItfkraLeCpwbEc+LiGkRsQ3wF7Kae72ryNub82aa5wK3tdjutcA/SdpUWRrYtzRbSNJ2EbEgIk4ha+5Yre28zq3ANEkvyKcPA37TuFCe8/0vkg7O9yFJuzbZ3kLgBQ3zLiFLhHZx3fWBVDsw3Iyl1i3h3jVWWYeQ5cyvdwHw7w3zvgFMlrSArFnk7XkO8aYi4l7gs2QjKf2arBnn4SaLHivpJkk3kNWCf9Fim8uAd5A1wywgq0V/c4LFDwWOyLd7M82HkrwK2L12UbduPz8kyxM/R1la3VR7Axd3sLyVSUWaa5yF0gojacOIeCxfWJ+/AAABmklEQVSvyV9Ilj668QtlqCR9GfhpRPy6D9u6CjggIh7svWRWJG0xOXjbhu0X/PojbbNQSpoJfJms6fNbEXFyXwqZyDV5K9In8wufN5E1Af1kyOVp5rPAM3rdSJ5T/TQH+BHWh5p8fo/G18mGDt2ZbOyBnQdc8lX4wqsVJiI+NOwytBMR99OHEcjynjZl/BKzFP3LXfMy4M8RcSeApO+TNRXe0o+Np3CQNzNrJu3C6pT8fo2aWRExq276OcA9ddOLgJf3oXTJHOTNzBqlX1hd2qZNvlk330IvhDrIm5k105/mmkXANnXTW5PdtV0YX3g1M2umP10orwO2l7StpHXIckH1fM2nE67Jm5k16tOF14hYIekY4JdkXSjPjoibe95wBxzkzcya6dPNThHxc+DnfdlYFxzkzcxWEyOTtqAdB3kzs0Ye49XMrOJGJDdNO85dY2bWIB+fd0rCoksjYuagy9MLB3kzswpzP3kzswpzkDczqzAHeTOzCnOQNzOrMAd5M7MKc5A3M6swB3kzswpzkDczqzAHeTOzCvv/c6UGEuTeqOcAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 2 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plot_coast(spacing, z)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.6.8"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}