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()