Skip to content

Instantly share code, notes, and snippets.

@jgalazm

jgalazm/bug_in_topotools.py

Last active March 18, 2016 02:44
Show Gist options
  • Save jgalazm/6cfd1bf5c7382ec8c6d1 to your computer and use it in GitHub Desktop.
Save jgalazm/6cfd1bf5c7382ec8c6d1 to your computer and use it in GitHub Desktop.
Download ZIP
(GeoClaw) Bug behavior in topotools.py plot function
Raw
bug_in_topotools.py
{
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# Jupyter notebook for Issue in GeoClaw"
]
},
{
"cell_type": "code",
"execution_count": 1,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"from clawpack.geoclaw import topotools\n",
"import numpy as np\n",
"%matplotlib inline\n",
"import matplotlib.pyplot as plt"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Bugged case"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"#1.\n",
"def topofunction(x,y):\n",
" \"\"\"\n",
" The topography: flat floor everywhere\n",
" except for the reservoir walls\n",
" \"\"\"\n",
" z = np.zeros(x.shape)\n",
" z = np.where( (np.abs(x-2.)<2.5)*(y<-70.), 15., z)\n",
" z = np.where( (np.abs(x-2.)<2.5)*(y>5.), 15.,z)\n",
" return z\n",
"\n",
"def maketopo(tfunction):\n",
" \"\"\"\n",
" topo2 format for topography\n",
" \"\"\" \n",
" xupper = 100.\n",
" yupper = 100.\n",
" xlower = -100.\n",
" ylower = -100.\n",
" nxpoints= 100\n",
" nypoints= 100\n",
" outfile= \"channel.topotype2\"\n",
" topotools.topo2writer(outfile,tfunction,xlower,xupper,ylower,yupper,nxpoints,nypoints)\n",
"\n",
"maketopo(topofunction)"
]
},
{
"cell_type": "code",
"execution_count": 8,
"metadata": {
"collapsed": false,
"scrolled": false
},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.axes._subplots.AxesSubplot at 0x7fa42a162a50>"
]
},
"execution_count": 8,
"metadata": {},
"output_type": "execute_result"
},
{
"name": "stdout",
"output_type": "stream",
"text": [
"Error in callback <function post_execute at 0x7fa42cea7500> (for post_execute):\n"
]
},
{
"ename": "ValueError",
"evalue": "posx and posy should be finite values",
"output_type": "error",
"traceback": [
"\u001b[1;31m---------------------------------------------------------------------------\u001b[0m",
"\u001b[1;31mValueError\u001b[0m Traceback (most recent call last)",
"\u001b[1;32m/home/jose/anaconda2/lib/python2.7/site-packages/matplotlib/pyplot.pyc\u001b[0m in \u001b[0;36mpost_execute\u001b[1;34m()\u001b[0m\n\u001b[0;32m 145\u001b[0m \u001b[1;32mdef\u001b[0m \u001b[0mpost_execute\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 146\u001b[0m \u001b[1;32mif\u001b[0m \u001b[0mmatplotlib\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mis_interactive\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 147\u001b[1;33m \u001b[0mdraw_all\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 148\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 149\u001b[0m \u001b[1;31m# IPython >= 2\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
"\u001b[1;32m/home/jose/anaconda2/lib/python2.7/site-packages/matplotlib/_pylab_helpers.pyc\u001b[0m in \u001b[0;36mdraw_all\u001b[1;34m(cls, force)\u001b[0m\n\u001b[0;32m 148\u001b[0m \u001b[1;32mfor\u001b[0m \u001b[0mf_mgr\u001b[0m \u001b[1;32min\u001b[0m \u001b[0mcls\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mget_all_fig_managers\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 149\u001b[0m \u001b[1;32mif\u001b[0m \u001b[0mforce\u001b[0m \u001b[1;32mor\u001b[0m \u001b[0mf_mgr\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mcanvas\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mfigure\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mstale\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 150\u001b[1;33m \u001b[0mf_mgr\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mcanvas\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mdraw_idle\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 151\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 152\u001b[0m \u001b[0matexit\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mregister\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mGcf\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mdestroy_all\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
"\u001b[1;32m/home/jose/anaconda2/lib/python2.7/site-packages/matplotlib/backend_bases.pyc\u001b[0m in \u001b[0;36mdraw_idle\u001b[1;34m(self, *args, **kwargs)\u001b[0m\n\u001b[0;32m 2024\u001b[0m \u001b[1;32mif\u001b[0m \u001b[1;32mnot\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0m_is_idle_drawing\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 2025\u001b[0m \u001b[1;32mwith\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0m_idle_draw_cntx\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m-> 2026\u001b[1;33m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mdraw\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m*\u001b[0m\u001b[0margs\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;33m**\u001b[0m\u001b[0mkwargs\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 2027\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 2028\u001b[0m \u001b[1;32mdef\u001b[0m \u001b[0mdraw_cursor\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mself\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mevent\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
"\u001b[1;32m/home/jose/anaconda2/lib/python2.7/site-packages/matplotlib/backends/backend_agg.pyc\u001b[0m in \u001b[0;36mdraw\u001b[1;34m(self)\u001b[0m\n\u001b[0;32m 472\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 473\u001b[0m \u001b[1;32mtry\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 474\u001b[1;33m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mfigure\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mdraw\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mrenderer\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 475\u001b[0m \u001b[1;32mfinally\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 476\u001b[0m \u001b[0mRendererAgg\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mlock\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mrelease\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
"\u001b[1;32m/home/jose/anaconda2/lib/python2.7/site-packages/matplotlib/artist.pyc\u001b[0m in \u001b[0;36mdraw_wrapper\u001b[1;34m(artist, renderer, *args, **kwargs)\u001b[0m\n\u001b[0;32m 59\u001b[0m \u001b[1;32mdef\u001b[0m \u001b[0mdraw_wrapper\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0martist\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mrenderer\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;33m*\u001b[0m\u001b[0margs\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;33m**\u001b[0m\u001b[0mkwargs\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 60\u001b[0m \u001b[0mbefore\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0martist\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mrenderer\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m---> 61\u001b[1;33m \u001b[0mdraw\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0martist\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mrenderer\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;33m*\u001b[0m\u001b[0margs\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;33m**\u001b[0m\u001b[0mkwargs\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 62\u001b[0m \u001b[0mafter\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0martist\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mrenderer\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 63\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n",
"\u001b[1;32m/home/jose/anaconda2/lib/python2.7/site-packages/matplotlib/figure.pyc\u001b[0m in \u001b[0;36mdraw\u001b[1;34m(self, renderer)\u001b[0m\n\u001b[0;32m 1157\u001b[0m \u001b[0mdsu\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0msort\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mkey\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mitemgetter\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;36m0\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 1158\u001b[0m \u001b[1;32mfor\u001b[0m \u001b[0mzorder\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0ma\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mfunc\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0margs\u001b[0m \u001b[1;32min\u001b[0m \u001b[0mdsu\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m-> 1159\u001b[1;33m \u001b[0mfunc\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m*\u001b[0m\u001b[0margs\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 1160\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 1161\u001b[0m \u001b[0mrenderer\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mclose_group\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;34m'figure'\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
"\u001b[1;32m/home/jose/anaconda2/lib/python2.7/site-packages/matplotlib/artist.pyc\u001b[0m in \u001b[0;36mdraw_wrapper\u001b[1;34m(artist, renderer, *args, **kwargs)\u001b[0m\n\u001b[0;32m 59\u001b[0m \u001b[1;32mdef\u001b[0m \u001b[0mdraw_wrapper\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0martist\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mrenderer\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;33m*\u001b[0m\u001b[0margs\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;33m**\u001b[0m\u001b[0mkwargs\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 60\u001b[0m \u001b[0mbefore\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0martist\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mrenderer\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m---> 61\u001b[1;33m \u001b[0mdraw\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0martist\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mrenderer\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;33m*\u001b[0m\u001b[0margs\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;33m**\u001b[0m\u001b[0mkwargs\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 62\u001b[0m \u001b[0mafter\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0martist\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mrenderer\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 63\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n",
"\u001b[1;32m/home/jose/anaconda2/lib/python2.7/site-packages/matplotlib/axes/_base.pyc\u001b[0m in \u001b[0;36mdraw\u001b[1;34m(self, renderer, inframe)\u001b[0m\n\u001b[0;32m 2322\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 2323\u001b[0m \u001b[1;32mfor\u001b[0m \u001b[0mzorder\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0ma\u001b[0m \u001b[1;32min\u001b[0m \u001b[0mdsu\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m-> 2324\u001b[1;33m \u001b[0ma\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mdraw\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mrenderer\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 2325\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 2326\u001b[0m \u001b[0mrenderer\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mclose_group\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;34m'axes'\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
"\u001b[1;32m/home/jose/anaconda2/lib/python2.7/site-packages/matplotlib/artist.pyc\u001b[0m in \u001b[0;36mdraw_wrapper\u001b[1;34m(artist, renderer, *args, **kwargs)\u001b[0m\n\u001b[0;32m 59\u001b[0m \u001b[1;32mdef\u001b[0m \u001b[0mdraw_wrapper\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0martist\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mrenderer\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;33m*\u001b[0m\u001b[0margs\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;33m**\u001b[0m\u001b[0mkwargs\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 60\u001b[0m \u001b[0mbefore\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0martist\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mrenderer\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m---> 61\u001b[1;33m \u001b[0mdraw\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0martist\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mrenderer\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;33m*\u001b[0m\u001b[0margs\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;33m**\u001b[0m\u001b[0mkwargs\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 62\u001b[0m \u001b[0mafter\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0martist\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mrenderer\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 63\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n",
"\u001b[1;32m/home/jose/anaconda2/lib/python2.7/site-packages/matplotlib/axis.pyc\u001b[0m in \u001b[0;36mdraw\u001b[1;34m(self, renderer, *args, **kwargs)\u001b[0m\n\u001b[0;32m 1118\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0m_update_label_position\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mticklabelBoxes\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mticklabelBoxes2\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 1119\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m-> 1120\u001b[1;33m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mlabel\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mdraw\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mrenderer\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 1121\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 1122\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0m_update_offset_text_position\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mticklabelBoxes\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mticklabelBoxes2\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
"\u001b[1;32m/home/jose/anaconda2/lib/python2.7/site-packages/matplotlib/artist.pyc\u001b[0m in \u001b[0;36mdraw_wrapper\u001b[1;34m(artist, renderer, *args, **kwargs)\u001b[0m\n\u001b[0;32m 59\u001b[0m \u001b[1;32mdef\u001b[0m \u001b[0mdraw_wrapper\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0martist\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mrenderer\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;33m*\u001b[0m\u001b[0margs\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;33m**\u001b[0m\u001b[0mkwargs\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 60\u001b[0m \u001b[0mbefore\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0martist\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mrenderer\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m---> 61\u001b[1;33m \u001b[0mdraw\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0martist\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mrenderer\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;33m*\u001b[0m\u001b[0margs\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;33m**\u001b[0m\u001b[0mkwargs\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 62\u001b[0m \u001b[0mafter\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0martist\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mrenderer\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 63\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n",
"\u001b[1;32m/home/jose/anaconda2/lib/python2.7/site-packages/matplotlib/text.pyc\u001b[0m in \u001b[0;36mdraw\u001b[1;34m(self, renderer)\u001b[0m\n\u001b[0;32m 755\u001b[0m \u001b[0mposy\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mfloat\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mtextobj\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mconvert_yunits\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mtextobj\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0m_y\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 756\u001b[0m \u001b[1;32mif\u001b[0m \u001b[1;32mnot\u001b[0m \u001b[0mnp\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0misfinite\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mposx\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;32mor\u001b[0m \u001b[1;32mnot\u001b[0m \u001b[0mnp\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0misfinite\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mposy\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 757\u001b[1;33m \u001b[1;32mraise\u001b[0m \u001b[0mValueError\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;34m\"posx and posy should be finite values\"\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 758\u001b[0m \u001b[0mposx\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mposy\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mtrans\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mtransform_point\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mposx\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mposy\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 759\u001b[0m \u001b[0mcanvasw\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mcanvash\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mrenderer\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mget_canvas_width_height\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
"\u001b[1;31mValueError\u001b[0m: posx and posy should be finite values"
]
},
{
"ename": "ValueError",
"evalue": "posx and posy should be finite values",
"output_type": "error",
"traceback": [
"\u001b[1;31m---------------------------------------------------------------------------\u001b[0m",
"\u001b[1;31mValueError\u001b[0m Traceback (most recent call last)",
"\u001b[1;32m/home/jose/anaconda2/lib/python2.7/site-packages/IPython/core/formatters.pyc\u001b[0m in \u001b[0;36m__call__\u001b[1;34m(self, obj)\u001b[0m\n\u001b[0;32m 335\u001b[0m \u001b[1;32mpass\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 336\u001b[0m \u001b[1;32melse\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 337\u001b[1;33m \u001b[1;32mreturn\u001b[0m \u001b[0mprinter\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mobj\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 338\u001b[0m \u001b[1;31m# Finally look for special method names\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 339\u001b[0m \u001b[0mmethod\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0m_safe_get_formatter_method\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mobj\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mprint_method\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
"\u001b[1;32m/home/jose/anaconda2/lib/python2.7/site-packages/IPython/core/pylabtools.pyc\u001b[0m in \u001b[0;36m<lambda>\u001b[1;34m(fig)\u001b[0m\n\u001b[0;32m 205\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 206\u001b[0m \u001b[1;32mif\u001b[0m \u001b[1;34m'png'\u001b[0m \u001b[1;32min\u001b[0m \u001b[0mformats\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 207\u001b[1;33m \u001b[0mpng_formatter\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mfor_type\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mFigure\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;32mlambda\u001b[0m \u001b[0mfig\u001b[0m\u001b[1;33m:\u001b[0m \u001b[0mprint_figure\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mfig\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;34m'png'\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;33m**\u001b[0m\u001b[0mkwargs\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 208\u001b[0m \u001b[1;32mif\u001b[0m \u001b[1;34m'retina'\u001b[0m \u001b[1;32min\u001b[0m \u001b[0mformats\u001b[0m \u001b[1;32mor\u001b[0m \u001b[1;34m'png2x'\u001b[0m \u001b[1;32min\u001b[0m \u001b[0mformats\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 209\u001b[0m \u001b[0mpng_formatter\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mfor_type\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mFigure\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;32mlambda\u001b[0m \u001b[0mfig\u001b[0m\u001b[1;33m:\u001b[0m \u001b[0mretina_figure\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mfig\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;33m**\u001b[0m\u001b[0mkwargs\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
"\u001b[1;32m/home/jose/anaconda2/lib/python2.7/site-packages/IPython/core/pylabtools.pyc\u001b[0m in \u001b[0;36mprint_figure\u001b[1;34m(fig, fmt, bbox_inches, **kwargs)\u001b[0m\n\u001b[0;32m 115\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 116\u001b[0m \u001b[0mbytes_io\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mBytesIO\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 117\u001b[1;33m \u001b[0mfig\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mcanvas\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mprint_figure\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mbytes_io\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;33m**\u001b[0m\u001b[0mkw\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 118\u001b[0m \u001b[0mdata\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mbytes_io\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mgetvalue\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 119\u001b[0m \u001b[1;32mif\u001b[0m \u001b[0mfmt\u001b[0m \u001b[1;33m==\u001b[0m \u001b[1;34m'svg'\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
"\u001b[1;32m/home/jose/anaconda2/lib/python2.7/site-packages/matplotlib/backend_bases.pyc\u001b[0m in \u001b[0;36mprint_figure\u001b[1;34m(self, filename, dpi, facecolor, edgecolor, orientation, format, **kwargs)\u001b[0m\n\u001b[0;32m 2178\u001b[0m \u001b[0morientation\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0morientation\u001b[0m\u001b[1;33m,\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 2179\u001b[0m \u001b[0mdryrun\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mTrue\u001b[0m\u001b[1;33m,\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m-> 2180\u001b[1;33m **kwargs)\n\u001b[0m\u001b[0;32m 2181\u001b[0m \u001b[0mrenderer\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mfigure\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0m_cachedRenderer\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 2182\u001b[0m \u001b[0mbbox_inches\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mfigure\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mget_tightbbox\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mrenderer\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
"\u001b[1;32m/home/jose/anaconda2/lib/python2.7/site-packages/matplotlib/backends/backend_agg.pyc\u001b[0m in \u001b[0;36mprint_png\u001b[1;34m(self, filename_or_obj, *args, **kwargs)\u001b[0m\n\u001b[0;32m 525\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 526\u001b[0m \u001b[1;32mdef\u001b[0m \u001b[0mprint_png\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mself\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mfilename_or_obj\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;33m*\u001b[0m\u001b[0margs\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;33m**\u001b[0m\u001b[0mkwargs\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 527\u001b[1;33m \u001b[0mFigureCanvasAgg\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mdraw\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mself\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 528\u001b[0m \u001b[0mrenderer\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mget_renderer\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 529\u001b[0m \u001b[0moriginal_dpi\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mrenderer\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mdpi\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
"\u001b[1;32m/home/jose/anaconda2/lib/python2.7/site-packages/matplotlib/backends/backend_agg.pyc\u001b[0m in \u001b[0;36mdraw\u001b[1;34m(self)\u001b[0m\n\u001b[0;32m 472\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 473\u001b[0m \u001b[1;32mtry\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 474\u001b[1;33m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mfigure\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mdraw\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mrenderer\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 475\u001b[0m \u001b[1;32mfinally\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 476\u001b[0m \u001b[0mRendererAgg\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mlock\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mrelease\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
"\u001b[1;32m/home/jose/anaconda2/lib/python2.7/site-packages/matplotlib/artist.pyc\u001b[0m in \u001b[0;36mdraw_wrapper\u001b[1;34m(artist, renderer, *args, **kwargs)\u001b[0m\n\u001b[0;32m 59\u001b[0m \u001b[1;32mdef\u001b[0m \u001b[0mdraw_wrapper\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0martist\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mrenderer\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;33m*\u001b[0m\u001b[0margs\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;33m**\u001b[0m\u001b[0mkwargs\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 60\u001b[0m \u001b[0mbefore\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0martist\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mrenderer\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m---> 61\u001b[1;33m \u001b[0mdraw\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0martist\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mrenderer\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;33m*\u001b[0m\u001b[0margs\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;33m**\u001b[0m\u001b[0mkwargs\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 62\u001b[0m \u001b[0mafter\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0martist\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mrenderer\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 63\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n",
"\u001b[1;32m/home/jose/anaconda2/lib/python2.7/site-packages/matplotlib/figure.pyc\u001b[0m in \u001b[0;36mdraw\u001b[1;34m(self, renderer)\u001b[0m\n\u001b[0;32m 1157\u001b[0m \u001b[0mdsu\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0msort\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mkey\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mitemgetter\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;36m0\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 1158\u001b[0m \u001b[1;32mfor\u001b[0m \u001b[0mzorder\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0ma\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mfunc\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0margs\u001b[0m \u001b[1;32min\u001b[0m \u001b[0mdsu\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m-> 1159\u001b[1;33m \u001b[0mfunc\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m*\u001b[0m\u001b[0margs\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 1160\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 1161\u001b[0m \u001b[0mrenderer\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mclose_group\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;34m'figure'\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
"\u001b[1;32m/home/jose/anaconda2/lib/python2.7/site-packages/matplotlib/artist.pyc\u001b[0m in \u001b[0;36mdraw_wrapper\u001b[1;34m(artist, renderer, *args, **kwargs)\u001b[0m\n\u001b[0;32m 59\u001b[0m \u001b[1;32mdef\u001b[0m \u001b[0mdraw_wrapper\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0martist\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mrenderer\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;33m*\u001b[0m\u001b[0margs\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;33m**\u001b[0m\u001b[0mkwargs\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 60\u001b[0m \u001b[0mbefore\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0martist\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mrenderer\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m---> 61\u001b[1;33m \u001b[0mdraw\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0martist\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mrenderer\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;33m*\u001b[0m\u001b[0margs\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;33m**\u001b[0m\u001b[0mkwargs\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 62\u001b[0m \u001b[0mafter\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0martist\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mrenderer\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 63\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n",
"\u001b[1;32m/home/jose/anaconda2/lib/python2.7/site-packages/matplotlib/axes/_base.pyc\u001b[0m in \u001b[0;36mdraw\u001b[1;34m(self, renderer, inframe)\u001b[0m\n\u001b[0;32m 2322\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 2323\u001b[0m \u001b[1;32mfor\u001b[0m \u001b[0mzorder\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0ma\u001b[0m \u001b[1;32min\u001b[0m \u001b[0mdsu\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m-> 2324\u001b[1;33m \u001b[0ma\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mdraw\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mrenderer\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 2325\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 2326\u001b[0m \u001b[0mrenderer\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mclose_group\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;34m'axes'\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
"\u001b[1;32m/home/jose/anaconda2/lib/python2.7/site-packages/matplotlib/artist.pyc\u001b[0m in \u001b[0;36mdraw_wrapper\u001b[1;34m(artist, renderer, *args, **kwargs)\u001b[0m\n\u001b[0;32m 59\u001b[0m \u001b[1;32mdef\u001b[0m \u001b[0mdraw_wrapper\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0martist\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mrenderer\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;33m*\u001b[0m\u001b[0margs\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;33m**\u001b[0m\u001b[0mkwargs\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 60\u001b[0m \u001b[0mbefore\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0martist\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mrenderer\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m---> 61\u001b[1;33m \u001b[0mdraw\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0martist\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mrenderer\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;33m*\u001b[0m\u001b[0margs\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;33m**\u001b[0m\u001b[0mkwargs\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 62\u001b[0m \u001b[0mafter\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0martist\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mrenderer\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 63\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n",
"\u001b[1;32m/home/jose/anaconda2/lib/python2.7/site-packages/matplotlib/axis.pyc\u001b[0m in \u001b[0;36mdraw\u001b[1;34m(self, renderer, *args, **kwargs)\u001b[0m\n\u001b[0;32m 1118\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0m_update_label_position\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mticklabelBoxes\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mticklabelBoxes2\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 1119\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m-> 1120\u001b[1;33m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mlabel\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mdraw\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mrenderer\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 1121\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 1122\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0m_update_offset_text_position\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mticklabelBoxes\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mticklabelBoxes2\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
"\u001b[1;32m/home/jose/anaconda2/lib/python2.7/site-packages/matplotlib/artist.pyc\u001b[0m in \u001b[0;36mdraw_wrapper\u001b[1;34m(artist, renderer, *args, **kwargs)\u001b[0m\n\u001b[0;32m 59\u001b[0m \u001b[1;32mdef\u001b[0m \u001b[0mdraw_wrapper\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0martist\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mrenderer\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;33m*\u001b[0m\u001b[0margs\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;33m**\u001b[0m\u001b[0mkwargs\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 60\u001b[0m \u001b[0mbefore\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0martist\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mrenderer\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m---> 61\u001b[1;33m \u001b[0mdraw\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0martist\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mrenderer\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;33m*\u001b[0m\u001b[0margs\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;33m**\u001b[0m\u001b[0mkwargs\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 62\u001b[0m \u001b[0mafter\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0martist\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mrenderer\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 63\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n",
"\u001b[1;32m/home/jose/anaconda2/lib/python2.7/site-packages/matplotlib/text.pyc\u001b[0m in \u001b[0;36mdraw\u001b[1;34m(self, renderer)\u001b[0m\n\u001b[0;32m 755\u001b[0m \u001b[0mposy\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mfloat\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mtextobj\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mconvert_yunits\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mtextobj\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0m_y\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 756\u001b[0m \u001b[1;32mif\u001b[0m \u001b[1;32mnot\u001b[0m \u001b[0mnp\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0misfinite\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mposx\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;32mor\u001b[0m \u001b[1;32mnot\u001b[0m \u001b[0mnp\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0misfinite\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mposy\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 757\u001b[1;33m \u001b[1;32mraise\u001b[0m \u001b[0mValueError\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;34m\"posx and posy should be finite values\"\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 758\u001b[0m \u001b[0mposx\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mposy\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mtrans\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mtransform_point\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mposx\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mposy\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 759\u001b[0m \u001b[0mcanvasw\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mcanvash\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mrenderer\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mget_canvas_width_height\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
"\u001b[1;31mValueError\u001b[0m: posx and posy should be finite values"
]
},
{
"data": {
"text/plain": [
"<matplotlib.figure.Figure at 0x7fa42a654e50>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"topo = topotools.Topography()\n",
"topo.read('channel.topotype2',topo_type=2)\n",
"topo.plot(add_colorbar=True)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Solution\n",
"Move the topography a bit upwards, so $min(topo(x,y))>0$."
]
},
{
"cell_type": "code",
"execution_count": 9,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"#1.\n",
"def topofunction2(x,y):\n",
" return topofunction(x,y)+1e-10\n",
"\n",
"maketopo(topofunction2)"
]
},
{
"cell_type": "code",
"execution_count": 10,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.axes._subplots.AxesSubplot at 0x7fa429fe2850>"
]
},
"execution_count": 10,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAUgAAAEJCAYAAAD2EL76AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAGtFJREFUeJzt3X28HVV97/HP9yQkPAUEaoAkSqAYBJHHEKkIpspTtRCk\nyotSi4BKeylKqxcR1FK0atEKPqKvROACN1wE9BZ8goAUhVLIA0+BBMhVCCEJASumIJBwcn73j5mT\nbA57cuac2WvP3iffN695Zc/svWavzdnne9asWbNGEYGZmb1WT90VMDPrVA5IM7MCDkgzswIOSDOz\nAg5IM7MCDkgzswItCUhJl0paJenBhm3bSZoj6VFJN0vatuG5cyUtkbRY0pGtqIOZWau1qgV5OXDU\ngG2fBm6NiD2A24BzASTtBZwA7An8GXCJJLWoHmZmLdOSgIyIO4HnBmyeAVyRP74COC5/fCxwTUT0\nRsQTwBJgWivqYWbWSin7IMdHxCqAiHgaGJ9vnwgsa3jd8nybmVlHaedJGl/TaGZdZXTCfa+StGNE\nrJK0E/BMvn058IaG103Kt72GJIeqWU0iotK5gQmvGxsrV68t+/KlETG5yvul0MqAVL70uxE4BbgQ\n+BBwQ8P22ZIuJju03h2YW7TTnU9vYQ07yPPzYdzUumuRzvPz4fZjDnzVtmNWLqipNq030n9+K2e2\nYB+r1zL/gnL/k6aeP3+X6u/Yei0JSElXA9OBHSQ9CZwP/AtwnaTTgKVkZ66JiEWSrgUWAa8AZ4Sn\nFDIbmbp8fEpLAjIiTip46vCC138Z+HIr3tvMOliXj+BL2QdpGzFmQt01SMufzwC3IG14xo7wXzB/\nPst0d0I6IM0snVEOSDOz5twHaWZWoLvz0QFpZgl1eQvS80GaWTpSuaVp0ddOo9jw3Ccl9UnaPmX1\nHZBmlo5KLs01m0YRSZOAI8guQEnKAWlm6VRoQRZMowhwMXB2ymr3cx+kmaXT09o+SEnHAssiYmE7\n5tl2QJpZOi0MMUlbAOeRHV6v39yyN2jCAWlm6RTE1/wlv2fB/1s91L39MTAZeCC/TcskYIGkaRHx\nzEZLDpMD0szSKWhBTp2yHVOnbLd+fdZNTxbuIV+IiIeAnTbsWo8DB0REs37KlvBJGjNLp8JZ7Hwa\nxbuAKZKelHTqgJdEcenWcAvSzBIafn5tZBrF/ud3G/bOS3JAmlk6nqzCzKxAl19q6IA0s3S6Ox8d\nkGaWkFuQZmYFHJBmZgW6Ox8dkGaWkFuQZmYFWjxZRbs5IM0sHbcgzcwKdHc+OiDNLCG3IM3MCnR3\nPjogzSwhtyDNzIo4IM3MmvNsPmZmBbo7Hx2QZpaQ+yDNzAo4IM3MCnR3PjogzSwhtyDNzAp0dz76\ntq9mllBPT7mlCUmXSlol6cGGbV+RtFjS/ZJ+KGmbpNVPuXMz28RVuC82cDlw1IBtc4C3RMR+wBLg\n3JbXuYED0szSkcotTUTEncBzA7bdGhF9+erdwKSU1XdAmlk61VqQgzkN+HnlOm6ET9KYWToFrcP5\nC59lwcLfVtitPgO8EhFXD3snJTggzSydgoCcus94pu4zfv36rGseGcIudQrwHuBd1So3OAekmXWy\nVx2ESzoaOBs4LCLWpH5zB6SZpVNhNh9JVwPTgR0kPQmcD5wHjAFuUdY6vTsizqhe0eYckGaWToUr\naSLipCabLx9+ZYbOAWlm6fhSQzOzAt2dj+kDUtITwGqgj+y0/DRJ2wE/AHYBngBOiIjVqetiZm3W\n5S3IdgwU7wOmR8T+ETEt3/Zp4NaI2AO4jcSXC5lZTdIOFE+uHQGpJu8zA7gif3wFcFwb6mFm7VZh\nsopO0I6aBdkp+XmSPpJv2zEiVgFExNPA+MLSZta9urwF2Y6TNIdExEpJrwfmSHqULDQbDVw3s5Gg\ny/sgkwdkRKzM/31W0r8B04BVknaMiFWSdgKeKSr//PwNj8dMgLETElfYbBO0ZgWsXZFgx92dj2kD\nUtKWQE9EvCBpK+BI4ALgRuAU4ELgQ8ANRfsYNzVlDc0MsoZHY+PjhXtbtGO3IDdqR+D/Sor8vWZH\nxBxJ84FrJZ0GLAVOSFwPM6uDA7JYRDwO7Ndk+++Aw1O+t5l1gO7OR19JY2YJ9XR3QjogzSwdH2Kb\nmRVwQJqZFejufHRAmllCbkGamRXo7nx0QJpZQm5BmpkV8DAfM7MCbkGamRXo7nxsy3yQZrapksot\nTYvqUkmrJD3YsG07SXMkPSrpZknbpqy+A9LM0qkQkGS3eD1qwLa23q7FAWlm6VSYUTwi7gSeG7C5\nrbdrcR+kmaXT+pM04xtv1yKp1O1aJI2NiDVDfTMHpJmlUzDMZ/7cFSyY15IpzJverkWSgA8AfwUc\nDPRJGgMsB34KzMynY9woB6SZpVPQgpz6tolMfdvE9euzLik9hXnZ27Xcni8XAA9ExLqsOhoP/Clw\nsaTrImL2xt7MAWlm6VQ/wh7YS1n2di1HNjukjohngB8AP8hblBvlgDSzdCr0QUq6GpgO7CDpSeB8\n4F+A6wa7XUtjOEraBphEQ95FxIMRsXawOjggzSydCi3IiDip4KnSt2uRdD5wOvA4G/orAzisTHkH\npJmlU/+lhicBuw3nDDY4IM0spfoD8mFgHOCANLMOU/9sPl8E7ssvV1wfkhFxfJnCDkgzS6f2fOQK\n4GJgIdA31MIOSDNLp/5D7Jci4qLhFnZAmlk69QfkryR9gWz8ZOMh9oPFRTZwQJpZOrXnI9Pyf6c3\nbPMwHzPrADW3ICPi0CrlPd2ZmaVTbT7ICm+rE/MJK4qenyzp7YPtxy1IM0unvmE+E4H7Jd0DLACe\nBTYHdic73P5v4JzBduKANLN0asrHiPiapG8ARwCHkPVFvgQsBj5cZqozcECaWUo19kFGRC/w83wZ\nFgekmaVT/1nsShyQZpZO/eMgK3FAmlk6NQekpJ6IGPIlhv08zMfM0ulRuSWdJZK+LGnKcAo7IM0s\nnQq3fW2RA4Engf8t6U5Jp0naumxhB6SZpVPTQPF+EfH7iPhuREwDPgt8AVgp6VJJuw5W3n2QZpZO\nB/RBAkcDpwJTgG8As4FDgZuAPTZW3gFpZunUfxJ7CXAn8K2I+FXD9mskDTphhQPSzNKpdlfDfwA+\nTDbR7ULg1DJ3Ihxg/4j472ZPRMQZgxV2QJpZOsPMR0kTgI8Bb46ItZJ+AJwIXDnEXY2R9ClgMq++\n7evpZQo7IM0snZ5K54FHAVtJ6gO2BFYMYx83AHeTHWavG2phB6SZpTPMFmRErJD0NbIhOi8CcyLi\n1mHsaquI+OTwauFhPmaW0jCH+Uh6HTAD2AWYAGwt6aRh1ODnko4cbvXdgjSzdApakPN/tZQFv1q6\nsZKHA7+JiN8BSPoR8Hbg6lJvKz1HdmsFAedIehFYm69HRGxfZj8OSDNLp+As9tR3TmbqOyevX5/1\nxTsHvuRJ4GBJm5PdbOvdwLwhvPMfDaWaRRyQZpbOMIf5RMRcSdcD9wGv5P/OHEL59SdkJB0LvIOs\nRXlHRPyk7H5qC0hJRwNfJ+sHvTQiLqyrLpbGMSsX1F0Fq1uFgeIRcQFwQaW3l74F7AVck286S9KR\nEfHxMuVrCcj88p9vkzWbVwDzJN0QEY/UUR8zS6S+e9L0OxzYKyICQNJlwENlC9d1FnsasCQilkbE\nK2TpPqOmuphZKjVPVgE8DkxqWN8Z+HXZwnUdYk8EljWsP8WGG3yb2UhR/4zimwOLJd1N1gd5MDA3\nPytORBy/scIdf5Lm+fkbHo+ZAGMn1FcXs5FqzQpYO5zrVAZTez7yxSqF6wrI5cAbG9Yn5dteY9zU\nttTHbJM2dkDj44V7W7PfYd/roEUi4hdVytfVBzkP2F3SLpLGkF2EfmNNdTGzRNbRV2pJRdJBku6W\ntFrSy5LWSGo6u08ztbQgI2KdpDOBOWwY5rO4jrqYWTopw6+kS4APkp0IngacQnb5Yim19UFGxKCz\n+ZpZd1tLb91V6ImIRyWNzkfMzJJ0H9ntFwbV8SdpzKx79Q59hrFW+0PejfeApC8BK8mmUSvFs/mY\nWTLriFJLQqeQ5dyZZPNBvgl4f9nCbkGaWTJ19kFKGgX8U0ScDLwMfG6o+3BAmlkydQZkfjJ4N0mb\n5f2PQ+aANLNkOuAs9q+BOyTdAPyhf2NEfLNMYQekmSWzpv6z2E/my5b5MiQOSDNLpu4WZEQMud+x\nkQPSzJJZV/Mwn/5JKQZYDcwHZg12n20P8zGzZDpgmM9TQC9wVb6sJTujvQ8wa7DCbkGaWTK99Z+k\n+ZOIOKh/RdK/AXMj4iBJiwYr7IA0s2Tq7oMExkmaFBFP5esTgHH54zWDFXZAmlkyHRCQnwL+U9Ij\nZLNTTgHOlLQVMHuwwg5IM0um7skqIuJGSbeQ3bgLYFFEvJQ//tfByjsgzSyZqn2QkrYFvg/sTTb/\n7mkRcc8Qyo8GTgUOyzfdLun7EVEquR2QZpZMC4b5fAP4WUR8IA+7oQ72/g6wFXBZvv5B4ADg9DKF\nHZBmlkyVITyStgEOjYhTAPJWX+nZwHMHR8S+DetzJD1QtrDHQZpZMhVvubAr8FtJl0u6V9JMSVsM\nsQp9kib3r+SPSx/3OyDNLJmKATma7HD4OxFxAPAi8OkhVuEcsskqbpX0C+CXwNllC/sQ28ySKTpJ\ns/T237L09v8arPhTwLKI6L/58/VkgVdaRMyRNAXYM9+0uOEs9qAckGaWTNEwn52nv46dp79u/fqd\nFzz2mtdExCpJyyRNiYjHgHcDg1790kjSWOCjwDuAIGtNzoqIQQeJgwPSzBJqwUDxjwOzJW0G/IZs\nyM5QXEF2xUz/ddcnkYXliWUKOyDNLJmqw3wi4gHgoEFfWGyfiNirYf2WMtdg9/NJGjNLppcotST0\ngKTGySoOBO4rW9gtSDNLpgOuxX4rcI+k3+TruwKL83tjR352vJAD0syS6YCAnFGlsAPSzJKpOyAj\n4teS3gIcmm+6IyIeLlvefZBmlswaekstqUg6E7gOeGO+XCvpjLLl3YI0s2TqbkGSTUoxLSJeAJD0\nJeAu4JIyhR2QZpZM3TftIpskt/HGXK/k20pxQJpZMolvyFVI0uh89p+ryM5i/zB/6n1kg8dLcUCa\nWTI1HmLPBQ6IiK9Iup3s6hmAv42IeWV34oA0s2RqvKvh+sPoiJhLFphD5oA0s2RqbEG+XtInip6M\niIvK7MQBaWbJ1HjTrlHA1gzhhEwzDkgzS6bGFuTKiPh81Z04IM0smd76hvlUajn2c0CaWTJ1DfMh\nm1y3MgekmSVT1yF2RPyuFftxQJpZMh1wqWElDkgzSyblRBTt4IA0s2TcgjQzK+CANDMr0AGz+VTi\ngDSzZKoO85HUA8wHnoqIY1tSqSFwQJpZMi2YrOIsYBGwTfXaDJ1vuWBmyayjr9TSjKRJwHuA77e1\n0g3cgjSzZCpOVnExcDawbWtqM3TJAlLS+cBHgWfyTedFxE35c+cCpwG9wFkRMSdVPcysPkWtw97b\n19F7e/Hht6T3Aqsi4n5J02nRtdVDlboFedHAedck7QmcAOwJTAJulfSmiKjtok0zS6OwD3K6GD19\n1PrVtRe8pqV5CHCspPcAWwDjJF0ZESenqWlzqfsgm6X+DOCaiOiNiCeAJcC0xPUwsxqsK/nfQBFx\nXkS8MSJ2A04Ebmt3OEL6gDxT0v2Svi+pvx9hIrCs4TXL821mNsKsI0otnarSIbakW4AdGzcBAXyG\n7L6zn4+IkPTPwNeAjwz1PZ6fv+HxmAkwdkKVGptZM2tWwNoVrd9vK66kiYhfAr+sXpuhqxSQEXFE\nyZfOAn6cP14OvKHhuUn5tqbGTR1e3cysvLEDGh8v3Nua/UZf57YOy0h2iC1pp4bV44GH8sc3AidK\nGiNpV2B3hnnHMTPrbHql3NKpUp7F/oqk/YA+4AngbwAiYpGka8lGx78CnOEz2GYjk0oeYXdqAKiT\ns0lS7Hx63bUw2/SsnAkRUWnsoaToeaHca/u2rv5+KfhKGjNLp7sn83FAmlk6ZQ+xO5UD0szScQvS\nzKyAW5BmZs118hCeMhyQZpaM+yDNzIq4D9LMrDm3IM3MijggzcwK+BDbzKw5H2KbmRXwMB8zsyJu\nQZqZFXAfpJlZc+6DNDMr4hakmVlzbkGaFfjxzge+av2YlQtqqonVpsJZbEmTgCvJ7pzaB8yKiG+2\npmLlOCDNLJmKLche4BMRcb+krYEFkuZExCMtqVwJDkgzS6dCH2REPA08nT9+QdJiYCLggDSz7teq\nPkhJk4H9gHtas8dyHJBmlk5BQPY+AOseLLeL/PD6euCsiCh5n8TWcECaWToFh9ij986Wfq/Mbv46\nSaPJwvGqiLih1dUbjAPSzJJpwSH2ZcCiiPhG9doMnQPSzJKpMlmFpEOAvwIWSroPCOC8iLipNbUb\nnAPSzNKp0IKMiP8ARrWsLsPggDSzdHypoZlZc77U0MysiFuQZmbNuQVpZlbEAWlm1pzvSWNmVsR9\nkGZmzbkP0sysiAPSzKyAD7HNzJrzIbaZWRG3IM3MmlNv3TWoxgFpZum4BWlm1pz7IM3MirgFaWbW\nnFuQZmZFHJBmZgV8iG1m1ly3z+bTU6WwpPdLekjSOkkHDHjuXElLJC2WdGTD9gMkPSjpMUlfr/L+\nZtbZ1FduaVpWOlrSI3lWnNPemmcqBSSwEHgf8MvGjZL2BE4A9gT+DLhEkvKnvwt8OCKmAFMkHVWx\nDl1pzYq6a5CWP58BWR9kmWUAST3At4GjgLcAfynpzW2pc4NKARkRj0bEEkADnpoBXBMRvRHxBLAE\nmCZpJ2BcRMzLX3clcFyVOnSrtSP8F8yfz4CsD7LM8lrTgCURsTQiXgGuIcuVtqragiwyEVjWsL48\n3zYReKph+1P5NjMbgSocYg/MkFqyYtCTNJJuAXZs3AQE8JmI+HGqipnZCDDSz2JHxBHD2O9y4A0N\n65PybUXbC62cOYx37xIv3Ft3DdKaeu+CuquQ1Ej/+bXA0pXfY5eSr101YH058MaG9UGzIoVWDvNp\n7Ie8EZgt6WKyZvHuwNyICEmrJU0D5gEnA98s2mFEDOzbNLMuERGTKxSfB+wuaRdgJXAi8JetqNdQ\nVB3mc5ykZcDBwE8k/RwgIhYB1wKLgJ8BZ0RE5MX+DrgUeIysE/amKnUws5EnItYBZwJzgIfJTvou\nbnc9tCG3zMysUaqz2LaJ6x/32jD+1azrdG1ASjpM0h5116OVJO0saYu66zFckkZJOkvS9WRdKcQI\nO0RpCP6DJF0k6XhJmzU+N5JI2lrSRyTtJ2lUvm3Efc4iXReQkvaX9ATwT8BFkt5eb42qk7R7Ppzq\nLuC7kvatu07DdDjZlQ8zgePzsPyjmuvUMpK2yU80TgcuA14m+8xfgRH5x2Afsj7A44DzgfNg5H3O\njen4gGzSqpoB/GNEvAv4CXBKfla8q/6ySdpN0th8dQawMCJ2BX4NnCNpx+LSHetkYE5EzAE+C+wM\nHFtvlaqRtKWkD0n6BXB5vnl/4HsRcR7wGWC6pHfmr++a7+BAeeNjSsOmtwGPRcSfA/8MvF3SQfXU\nrh4dG5BNWlX9k2G8BXh9/viHwFLgL/qLtbeWQyNpb0nflvQQcCowJn9qM2B1/vg7ZMNrD+nCX7b/\nBCbnj+8ju8T0gP5Ds26THzo/DLwf+GpE9H/P9ie7WIKIeI5sWNvH+4u1u55V5b9rdwL/C/hHSSfn\nT/UBT0jaLCIWAIvJvpfb1FTVtuvYgOS1raqz8x/Mj8j+sgH8lixA95S0ZUR07PSc+XiuC4FdgD+N\niM9FxPP505sDz+Sf4Xdkv5R7AzvUU9th+w2wtaQdIuKlfD2Aveqt1vDk1wDfRzbEpHE42k/JWsv9\nvgcckpfp2O9gP0lbSfqThk1vBhZExL7At4CPS5oM9JJlxPj8df9BNqa5276Xw1Z7QEp6s6TrJP2r\npDENTw1sVfWS9fc8AmwpaXL+ZXyGbCDpm9pZ72FYRhbm10XEs5K2l7Rt/tyTwG5kVwsAzCX70q5t\nfzUrWQisAd6drz9Hdpnqs7XVqLrLgPPzEzL/LulzwN1kM1GNB4iIlcADkt5RZ0XLkHQe2R+unzV0\n4xwFPAoQEfeQff8+BtxGdqHHH+evuwt4K/CHdta5TrUGpKTtyA6PxwEH8eqL0Qe2qhaRTZ/2LPA4\n8JH8dS8BO/HqC9s7iqTReZgvBj4gaSHZYdnnJR0O3ED2/2AqQETcBhxKfhjXRZaThUf/4eazZAH5\nfGGJDhcRPyH7vv0XWatxb7KTFo8Cf6vM9mR/5JbWVtHy7gKOIOue+mC+7XHgvQ2v+Q5wQkQsA54m\na5j0/yHYkuz3bZNQa0Dm/TezI+Josi/g1Ib+qmVk/VmNraq9yX5gVwLvlfQhsk7y3+ch2jHyX5xR\nABHRf/v0u4E7gLPJvnT3AV8lC5AfAWdIOlzS35FdgdRVl/pHRF9EXAUslfRTss93PVmrspu9LyK+\nmAfGl4ApZIeiW5D9ofsFMCp/vtPdEREPArew4QTaZWR9i5sDRMTDwNOS9if7fu4paWZ+pdwi8tbm\nJiEialmAngH//j1ZX842+foUsr9kJzWUWdbw/H75678AbF/X5yj4bBqwfnD+WY4Htmj8fwDcDLwr\nXz+e7BLNm4GD6v4cFT7/ZsC+jZ91pCxkfci39H/nyFph+9ddr2F8jh2Ae4F98/VbgbManp8JHJY/\n3gn4H8DfAGPrrns7l7Zeapifle2J7DrL9dsiIiTtSvaX7CzgoYjok3QM8CngAmAPYB/gExHR8X0g\n+ec5mWxG9ZfJvoDfiojVDa/5C+DPgY9FxAv5tp7ogo7+TUk+HOto4K/JRlF8F7gkNhwZdCVJlwAv\nRsT/lHQYcAYwH9gOOBB4b+Pv6qaobTft6g9C8sNGSQeTfeFuk3RDRDwuaSXw1sgOAYiIH+dfztOB\nbYHPdkk4Hg58neww+aMRsbDhuS3IWsvHkJ14+l5/OEJ3nAXd1ETEGmWz4d8MfDAiXq67Ti0yE/hm\nPpzp92RjVz9N9gf9k5t6OEKbJ6sYrFUl6TjgE8AKYFlEnJ1v76pWVX5SpnfAtlFkFyH0SToBeDQi\nHqinhmYg6UTgauBFsq6qr3bT71k7tC0gB7SqrhrQquohG9v4f8h+WNcAMyPi6bZULpE8FPuinX+F\nzErILyO8kOwk2uwR1CpuqXYGZFGrqi/vg9wd2LL/8NrMrG5tnw/SrSoz6xaeMNfMrEDtlxqamXUq\nB6SZWQEHpJlZAQekmVkBB6SZWQEHpJlZAQekmVkBB6SZWQEHpJlZgf8PCqXw8mDUtMIAAAAASUVO\nRK5CYII=\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x7fa42a162390>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"topo = topotools.Topography()\n",
"topo.read('channel.topotype2',topo_type=2)\n",
"topo.plot(add_colorbar=True)"
]
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 2",
"language": "python",
"name": "python2"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 2
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython2",
"version": "2.7.11"
}
},
"nbformat": 4,
"nbformat_minor": 0
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment