source: anuga_validation/okushiri_2005/run_okushiri_parallel.py @ 4047

Last change on this file since 4047 was 4022, checked in by ole, 18 years ago

Removed parallelism from run_okushiri. Use run_okushiri_parallel instead.
It caused problems for the distribution which does not include parallelism.

File size: 2.8 KB
RevLine 
[3640]1"""Validation of the AnuGA implementation of the shallow water wave equation.
2
[4022]3This script sets up Okushiri Island benchmark as published at the
[3640]4
[4022]5THE THIRD INTERNATIONAL WORKSHOP ON LONG-WAVE RUNUP MODELS
6June 17-18 2004
7Wrigley Marine Science Center
8Catalina Island, California
9http://www.cee.cornell.edu/longwave/
[3640]10
[4022]11
12The validation data was downloaded and made available in this directory
13for convenience but the original data is available at
[3640]14http://www.cee.cornell.edu/longwave/index.cfm?page=benchmark&problem=2
[4022]15where a detailed description of the problem is also available.
[3640]16
17
[4022]18Run create_okushiri.py to process the boundary condition and build a the
19mesh before running this script.
20
[3640]21"""
22
23# Module imports
24from anuga.shallow_water import Domain
25from anuga.shallow_water import Reflective_boundary
26from anuga.shallow_water import Transmissive_Momentum_Set_Stage_boundary
27from anuga.abstract_2d_finite_volumes.util import file_function
28
[4022]29# from anuga_parallel.parallel_api import myid, numprocs, distribute   
[3640]30
31import project
32
33
34#-------------------------
[4022]35# Create Domain from mesh
[3640]36#-------------------------
37domain = Domain(project.mesh_filename, use_cache=True, verbose=True)
[4022]38print domain.statistics()
[3640]39
[4022]40
[3640]41#-------------------------
42# Initial Conditions
43#-------------------------
44domain.set_quantity('friction', 0.0)
45domain.set_quantity('stage', 0.0)
46domain.set_quantity('elevation',
[4022]47                    filename=project.bathymetry_filename,
48                    alpha=0.02,                   
49                    verbose=True,
50                    use_cache=True)
[3640]51
52
53#-------------------------
[4022]54# Distribute domain if run in parallel
[3640]55#-------------------------
[4022]56if numprocs > 1:
57    domain = distribute(domain)
[3640]58
[3648]59
[4022]60#-------------------------
61# Set simulation parameters
62#-------------------------
63domain.set_name(project.output_filename)  # Name of output sww file
64domain.set_default_order(2)               # Apply second order scheme
65domain.set_all_limiters(0.9)              # Max second order scheme (old lim)
66domain.set_minimum_storable_height(0.001) # Don't store w < 0.001m
67domain.set_maximum_allowed_speed(0.1)     # Allow a little runoff (0.1 is OK)
[3829]68
69
[4022]70#-------------------------
71# Boundary Conditions
72#-------------------------
[3829]73
[4022]74# Create boundary function from timeseries provided in file
75function = file_function(project.boundary_filename,
76                         domain, verbose=True)
[3829]77
[4022]78# Create and assign boundary objects
[3645]79Bts = Transmissive_Momentum_Set_Stage_boundary(domain, function)
[4022]80Br = Reflective_boundary(domain)
81domain.set_boundary({'wave': Bts, 'wall': Br})
[3645]82
[3648]83
[3640]84#-------------------------
85# Evolve through time
86#-------------------------
87import time
88t0 = time.time()
89
90for t in domain.evolve(yieldstep = 0.05, finaltime = 22.5):
91    domain.write_time()
92
[4022]93print 'That took %.2f seconds' %(time.time()-t0)
Note: See TracBrowser for help on using the repository browser.