def run_and_plot(**kwargs):
claw = setup(kwargs)
claw.run()
from clawpack.pyclaw import plot
plot.interactive_plot(setplot=setplot)
plotitem = plotaxes.new_plotitem(plot_type='2d_pcolor')
plotitem.plot_var = 1
plotitem.pcolor_cmap = colormaps.yellow_red_blue
plotitem.add_colorbar = True
plotitem.show = True # show on plot?
# Parameters used only when creating html and/or latex hardcopy
# e.g., via visclaw.frametools.printframes:
plotdata.printfigs = True # print figures
plotdata.print_format = 'png' # file format
plotdata.print_framenos = 'all' # list of frames to print
plotdata.print_fignos = 'all' # list of figures to print
plotdata.html = True # create html files of plots?
plotdata.html_homelink = '../README.html' # pointer for top of index
plotdata.latex = True # create latex file of plots?
plotdata.latex_figsperline = 2 # layout of plots
plotdata.latex_framesperline = 1 # layout of plots
plotdata.latex_makepdf = False # also run pdflatex?
return plotdata
if __name__ == '__main__':
claw = inclusion()
claw.run()
from clawpack.pyclaw import plot
plot.interactive_plot()
def run_and_plot(**kwargs):
claw = setup(kwargs)
claw.run()
from clawpack.pyclaw import plot
plot.interactive_plot(setplot=setplot)
aux[0,:,:] = rho1 * inbar + rho2 * outbar
aux[1,:,:] = lam1 * inbar + lam2 * outbar
aux[2,:,:] = mu1 * inbar + mu2 * outbar
bulk = aux[1,:,:] + 2.*aux[2,:,:]
aux[3,:,:] = np.sqrt(bulk/aux[0,:,:])
aux[4,:,:] = np.sqrt(aux[2,:,:]/aux[0,:,:])
aux[5,:,:] = 0.
aux[6,:,:] = 0.
# set initial condition
state.q[:,:,:] = 0.
claw = pyclaw.Controller()
claw.solver = solver
claw.solution = solution
claw.num_output_times = 20
claw.tfinal = 0.5
return claw
if __name__ == '__main__':
claw = inclusion()
claw.run()
from clawpack.pyclaw import plot
plot.interactive_plot()
