Changeset 1761 for inundation/validation/LWRU2

Timestamp:

Aug 25, 2005, 10:06:22 AM (20 years ago)

Author:

steve

Message:

Added a more symmetric rectangle_cross procedure. Can be used in place of rectangle from mesh_factory

File:

• inundation/validation/LWRU2/lwru2.py (modified) (5 diffs)

inundation/validation/LWRU2/lwru2.py

@@ -74 +74 @@
 from pyvolution.shallow_water import Domain, Reflective_boundary,\
             File_boundary, Transmissive_Momentum_Set_Stage_boundary
-from pyvolution.mesh_factory import rectangular
+from pyvolution.mesh_factory import rectangular_cross
 from pyvolution.pmesh2domain import pmesh_to_domain_instance
 from Numeric import array, zeros, Float, allclose
@@ -91 +91 @@
 #######################
 # Domain
-    
+
 if read_mesh is True:
     print 'Creating domain from', filenames.mesh_filename
@@ -101 +101 @@
 
 
-    
-else:   
+
+else:
     print 'Creating regular mesh'
     N = 100
-    points, vertices, boundary = rectangular(N, N/5*3,
+    points, vertices, boundary = rectangular_cross(N, N/5*3,
                                              len1=5.448, len2=3.402)
     print 'Creating domain'
-    
+
     #domain = Domain(points, vertices, boundary)
     domain = cache(Domain, (points, vertices, boundary))
@@ -135 +135 @@
                     verbose = True,
                     use_cache = True)
-                    
+
 domain.set_quantity('friction', 0.0)
 domain.set_quantity('stage', 0.0)
@@ -157 +157 @@
 #Set boundary conditions
 if read_mesh is True:
-    domain.set_boundary({'wave': Bts, 'wall': Br})      
-else:    
+    domain.set_boundary({'wave': Bts, 'wall': Br})
+else:
     domain.set_boundary({'left': Bts, 'right': Br, 'bottom': Br, 'top': Br})