def plot_analytic(current_data):
     import pylab 
     
     N = 2000                           # number of cells
     X = linspace(-1000.,1000.,N+1)     # vertices of cells
     t = current_data.t
     (u,h,w,z,mom) = analytical_sol(X,t)     
     pylab.plot(X,w,'r--')
     pylab.legend(['Topography', 'GeoClaw Solution', 'Analytic Solution'], loc='upper left')
from numpy import *
from matplotlib import *
from pylab import *
from analytical_solution_landslide import analytical_sol


N = 100                        # number of cells
X = linspace(-1000.,1000.,N+1) # vertices of cells
t_compare = 15
(u,h,w,z,mom) = analytical_sol(X,t_compare)

figure(1)
plot(X,w,X,z)
title('State of the slide at time t = %s' %t_compare)
legend(('stage','topography'))
show()