def produce_report(script, args=None): import anuga if args is None: args = anuga.get_args() verbose = args.verbose #print args # Get the arguments from the calling script run_anuga_script(script, args=args) # We don't want to run plot_results in parallel args.np = 1 run_anuga_script('plot_results.py', args=args) typeset_report(verbose=verbose)
from numpy import zeros, float from time import localtime, strftime, gmtime from anuga import myid, finalize, distribute #------------------------------------------------------------------------------- # Copy scripts to time stamped output directory and capture screen # output to file #------------------------------------------------------------------------------- time = strftime('%Y%m%d_%H%M%S',localtime()) #output_dir = 'steep_island'+time output_dir = '.' output_file = 'steep_island' args = anuga.get_args() alg = args.alg verbose = args.verbose dx = 1. dy = dx L = 2000. W = 5*dx def stage_flat(x,y): w=zeros(len(x)) for i in range(len(x)): w[i]=4.5 return w def bed_elevation(x,y):
""" Runup example from the manual, slightly modified """ #--------- #Import Modules #-------- import anuga import numpy from math import sin, pi, exp from anuga import Domain from anuga import myid, finalize, distribute args = anuga.get_args() alg = args.alg verbose = args.verbose if myid == 0: #--------- #Setup computational domain #--------- points, vertices, boundary = anuga.rectangular_cross(100, 3, len1=1.0, len2=0.03) domain = Domain(points, vertices, boundary) # Create Domain domain.set_name('runup') # Output to file runup.sww domain.set_datadir('.') # Use current folder domain.set_quantities_to_be_stored({ 'stage': 2, 'xmomentum': 2, 'ymomentum': 2,