Changeset 1736

Timestamp:

Aug 20, 2005, 5:57:06 PM (20 years ago)

Author:

ole

Message:

 

Location:

inundation/validation/LWRU2

Files:

• extract_timeseries.py (modified) (1 diff)
• lwru2.py (modified) (4 diffs)
• output_ch5-7-9.txt (added)

inundation/validation/LWRU2/extract_timeseries.py

@@ -11 +11 @@
 #gauges = [[4.521, 1.196],  [4.521, 1.696],  [4.521, 2.196]] #Ch 5-7-9
 
-gauge_names = ['Boundary gauge', 'ch5', 'ch7', 'ch9']
+gauge_names = ['Boundary', 'ch5', 'ch7', 'ch9']
 #gauge_names = ['ch5', 'ch7', 'ch9']
 

inundation/validation/LWRU2/lwru2.py

@@ -9 +9 @@
 Depth at western boundary is d = 13.5 cm
 """
+
+
+read_mesh = False
 
 import sys
@@ -78 +81 @@
 from cPickle import dump, load
 
-N = 50
+N = 100
 
 try:
@@ -99 +102 @@
     #######################
     # Domain
-    #print 'Creating domain'
-    #points, vertices, boundary = rectangular(N, N/5*3,
-    #                                         len1=5.448, len2=3.402)
-    #domain = Domain(points, vertices, boundary)
+    
+    if read_mesh is True:
+        print 'Creating domain from', filenames.mesh_filename
+        domain = pmesh_to_domain_instance(filenames.mesh_filename, Domain)    
+    else:       
+        print 'Creating regular domain'
+        points, vertices, boundary = rectangular(N, N/5*3,
+                                                 len1=5.448, len2=3.402)
+        domain = Domain(points, vertices, boundary)
 
-    print 'Creating domain from', filenames.mesh_filename
-    domain = pmesh_to_domain_instance(filenames.mesh_filename, Domain)
+
 
     print "Number of triangles = ", len(domain)
@@ -159 +166 @@
 
 #Set boundary conditions
-#domain.set_boundary({'left': Bts, 'right': Br, 'bottom': Br, 'top': Br})
-domain.set_boundary({'wave': Bts, 'wall': Br})
+
+if read_mesh is True:
+    domain.set_boundary({'wave': Bts, 'wall': Br})      
+else:    
+    domain.set_boundary({'left': Bts, 'right': Br, 'bottom': Br, 'top': Br})
+
 
 #Evolve