sww = swmodel2d.ShallowWaterModel2D(
Lx=1000.,
Ly=1000.,
nx=101,
ny=101,
h0=10.,
dt=0.8,
epsilon=0.2)
swmodel2d.set_initial_conditions_dam_break(sww,width=100.,height=1.)
swwani = animator.SWWAnimator(sww=sww,ax=ax,zlim=0.4)
anim = animation.FuncAnimation(fig,func=swwani.advance_animation,fargs=(ax, fig),
frames=500, interval=1)
if False:
anim.save('shallow_water_wave_2D.mp4', fps=1./sww.dt, extra_args=['-vcodec', 'libx264'])
#anim.save('shallow_water_wave_2D.gif', writer='imagemagick', fps=1./sww.dt)
else:
plt.show()
nx=201,
ny=101,
h0=10.,
dt=100,
epsilon=0.1,
t0=0.,
Kappa=1e-7,
f=5e-4,
)
# beta-plane
swmodel2d.set_beta_plane(sww,beta=1e-10)
# initial conditions
swmodel2d.set_initial_conditions_dam_break(sww,
width=2e4,height=3,
xfrac=0,yfrac=1,
smooth=False)
swwani = animator.SWWAnimator(sww=sww,ax=ax,zlim=0.4)
anim = animation.FuncAnimation(fig,func=swwani.advance_animation,fargs=(ax, fig),
frames=500, interval=10, blit=False)
if False:
anim.save('shallow_water_wave_2D.mp4', fps=1./sww.dt, extra_args=['-vcodec', 'libx264'])
#anim.save('shallow_water_wave_2D.gif', writer='imagemagick', fps=1./sww.dt)
else:
plt.show()
