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