def plot_u(u, x, xv, y, yv, t, n):
#surf(xv, yv, u, show=True, zlim=[-2.4, 2.4])
#surf(xv, yv, u, zlim=[-50.0, 300.0])
#axis([0.0,400.0,0.0,430.0,-500.0,300.0])
#show()
st.surf(xv, yv, u, zlim=[zmin, zmax])
def plot_u(u, x, xv, y, yv, t, n):
#surf(xv, yv, u, show=True, zlim=[-2.4, 2.4])
#surf(xv, yv, u, zlim=[-50.0, 300.0])
#axis([0.0,400.0,0.0,430.0,-500.0,300.0])
#show()
st.surf(xv, yv, u, zlim=[zmin, zmax])
def plot_u(u, x, xv, y, yv, t, n):
#surf(xv, yv, u, show=True, zlim=[-2.4, 2.4])
#surf(xv, yv, u, zlim=[-50.0, 300.0])
#axis([0.0,400.0,0.0,430.0,-500.0,300.0])
#show()
st.surf(xv, yv, u, zlim=[zmin, zmax])
#hold('on')
#print 'xvshape = ',xv.shape[0],',yvshape = ', yv.shape[1]
b_a = zeros((xv.shape[0], yv.shape[1]))
b_a[:, :] = bottom(xv, yv)
#print 'b_a.shape=',b_a.shape
#surf(xv, yv, b_a)
st.mesh(xv, yv, b_a, zlim=[zmin, zmax])
#axis([0.0,400.0,0.0,430.0,-500.0,300.0])
#show()
#time.sleep(1.0)
filename = 'tmp_gaussian2%04d.png' % n
st.savefig(filename)
def plot_u(u, x, xv, y, yv, t, n):
#surf(xv, yv, u, show=True, zlim=[-2.4, 2.4])
#surf(xv, yv, u, zlim=[-50.0, 300.0])
#axis([0.0,400.0,0.0,430.0,-500.0,300.0])
#show()
st.surf(xv, yv, u, zlim=[zmin, zmax])
#hold('on')
#print 'xvshape = ',xv.shape[0],',yvshape = ', yv.shape[1]
b_a = zeros((xv.shape[0],yv.shape[1]))
b_a[:,:] = bottom(xv, yv)
#print 'b_a.shape=',b_a.shape
#surf(xv, yv, b_a)
st.mesh(xv, yv, b_a, zlim=[zmin, zmax])
#axis([0.0,400.0,0.0,430.0,-500.0,300.0])
#show()
#time.sleep(1.0)
filename = 'tmp_gaussian2%04d.png' % n
st.savefig(filename)
def plot_u(u, x, xv, y, yv, t, n):
#surf(xv, yv, u, show=True, zlim=[-2.4, 2.4])
#surf(xv, yv, u, zlim=[-50.0, 300.0])
#axis([0.0,400.0,0.0,430.0,-500.0,300.0])
#show()
st.surf(xv, yv, u, zlim=[zmin, zmax])
#show()
#hold('on')
#1 + 1
#print 'xvshape = ',xv.shape[0],',yvshape = ', yv.shape[1]
b_a = zeros((xv.shape[0], yv.shape[1]))
#print 'xv=',xv,',yv=',yv
b_a[:, :] = bottom(xv, yv)
#print 'b_a=',b_a
#print 'shape(b_a)=',shape(b_a)
#print 'b_a.shape=',b_a.shape
#surf(xv, yv, b_a)
#mesh(xv, yv, b_a)
st.mesh(xv, yv, b_a, zlim=[zmin, zmax])
def plot_u(u, x, xv, y, yv, t, n):
#surf(xv, yv, u, show=True, zlim=[-2.4, 2.4])
#surf(xv, yv, u, zlim=[-50.0, 300.0])
#axis([0.0,400.0,0.0,430.0,-500.0,300.0])
#show()
st.surf(xv, yv, u, zlim=[zmin, zmax])
#show()
#hold('on')
#1 + 1
#print 'xvshape = ',xv.shape[0],',yvshape = ', yv.shape[1]
b_a = zeros((xv.shape[0],yv.shape[1]))
#print 'xv=',xv,',yv=',yv
b_a[:,:] = bottom(xv, yv)
#print 'b_a=',b_a
#print 'shape(b_a)=',shape(b_a)
#print 'b_a.shape=',b_a.shape
#surf(xv, yv, b_a)
#mesh(xv, yv, b_a)
st.mesh(xv, yv, b_a, zlim=[zmin, zmax])