Changeset 3829 for anuga_validation/okushiri_2005/okushiri_parallel.py

Timestamp:

Oct 19, 2006, 6:21:11 PM (18 years ago)

Author:

ole

Message:

Retired modify_boundary (now rolled into set_boundary).
Got okushiri_parallel to run correctly and timed the speedup (6.5 on 8 procs)

File:

• anuga_validation/okushiri_2005/okushiri_parallel.py (modified) (4 diffs)

anuga_validation/okushiri_2005/okushiri_parallel.py

@@ -29 +29 @@
 #-------------------------
 domain = Domain(project.mesh_filename, use_cache=True, verbose=True)
-domain.set_name('okushiri_parallel')
-domain.set_minimum_storable_height(0.001)
-domain.set_default_order(2)
-print domain.statistics()
-
 
 #-------------------------
@@ -46 +41 @@
                     use_cache = True)
 
-#-------------------------
-# Boundary Conditions
-#-------------------------
-Br = Reflective_boundary(domain)
-
-domain.set_boundary({'wave': None,  # Bind this one later
-                     'wall': Br})
-
 
 #-------------------------
@@ -60 +47 @@
 domain = distribute(domain)
 
+# Parameters
+domain.set_name('okushiri_parallel')
+domain.set_minimum_storable_height(0.001)
+domain.set_default_order(2)
+domain.set_maximum_allowed_speed(0.1) # Allow a little runoff (0.1 is OK)
 
-# Bind boundary, that cannot be fully specified before distribution.
+# Set old (pre Sep 2006) defaults for limiters
+domain.beta_w      = 0.9
+domain.beta_w_dry  = 0.9
+domain.beta_uh     = 0.9
+domain.beta_uh_dry = 0.9
+domain.beta_vh     = 0.9
+domain.beta_vh_dry = 0.9
+
+
+#------------------------------------------------------------------------------
+# Setup boundary conditions
+# (MUST currently happen after domain has been distributed)
+#------------------------------------------------------------------------------
+
+Br = Reflective_boundary(domain)      # Solid reflective wall
+
 function = file_function(project.boundary_filename[:-4] + '.tms',
                          domain, 
@@ -67 +74 @@
 
 Bts = Transmissive_Momentum_Set_Stage_boundary(domain, function)
-domain.modify_boundary({'wave': Bts})
+
+domain.set_boundary({'wave': Bts,  # Bind this one later
+                     'wall': Br})