Changeset 1734

Timestamp:

Aug 19, 2005, 5:45:54 PM (19 years ago)

Author:

steve

Message:

 

Location:

inundation

Files:

• validation/LWRU2/create_mesh.py (modified) (9 diffs)
• validation/LWRU2/lwru2.py (modified) (2 diffs)
• zeus/anuga-workspace.zwi (modified) (1 diff)
• zeus/pyvolution.zbi (deleted)

inundation/validation/LWRU2/create_mesh.py

@@ -14 +14 @@
 
     #Basic geometry
-    
+
     xleft   = 0
     xright  = 5.448
@@ -20 +20 @@
     ytop    = 3.402
 
-    #Outline 
+    #Outline
     point_sw = [xleft, ybottom]
     point_se = [xright, ybottom]
-    point_nw = [xleft, ytop]    
+    point_nw = [xleft, ytop]
     point_ne = [xright, ytop]
 
@@ -29 +29 @@
     point_ml1 = [xleft + (xright-xleft)/3+1, ytop]
     point_ml2 = [xleft + (xright-xleft)/3+1, ybottom]
-    
+
     #Midway points (right)
     point_mr1 = [xleft + 2*(xright-xleft)/3+0.5, ytop]
@@ -37 +37 @@
     geo = Geo_reference(xllcorner = xleft,
                         yllcorner = ybottom)
-    
-                        
+
+
     print "***********************"
     print "geo ref", geo
     print "***********************"
-    
+
     m = Mesh(geo_reference=geo)
 
@@ -50 +50 @@
                       point_ne,
                       point_mr1,
-                      point_ml1,                      
+                      point_ml1,
                       point_nw,
                       point_sw,
-                      point_ml2,                                            
+                      point_ml2,
                       point_mr2]
 
-    
+
     dict['segments'] = [[0,1], [1,2], [2,3], [3,4],
                         [4,5], [5,6], [6,7], [7,0], #The outer border
                         [2,7], [3,6]]               #Separators
-    
+
     dict['segment_tags'] = ['wall',
                             'wall',
@@ -68 +68 @@
                             'wall',
                             'wall',
-                            'wall',                                  
-                            '',           #Interior                            
+                            'wall',
+                            '',           #Interior
                             '']           #Interior
 
-        
+
     m.addVertsSegs(dict)
 
@@ -87 +87 @@
     p2 = [xr, yt]
     p3 = [xl, yt]
-    
+
     dict['points'] = [p0, p1, p2, p3]
 
-        
+
     dict['segments'] = [[0,1], [1,2], [2,3], [3,0]]
     dict['segment_tags'] = ['', '', '', '']
-    m.addVertsSegs(dict)    
+    m.addVertsSegs(dict)
 
-    
-    base_resolution = 100 
+
+    base_resolution = 100
 
     ocean = m.addRegionEN(xleft+1, ybottom+1)
@@ -103 +103 @@
     mid = m.addRegionEN(point_ml2[0]+1, ybottom+1)
     mid.setMaxArea(0.001*base_resolution)
-    
+
 
     inner = m.addRegionEN(point_mr2[0]+1, ybottom+1)
@@ -110 +110 @@
 
     gulleys = m.addRegionEN(p0[0]+0.1, p0[1]+0.1)
-    gulleys.setMaxArea(0.00001*base_resolution)    
-    
-    
+    gulleys.setMaxArea(0.00001*base_resolution)
+
+
     m.generateMesh('pzq28.0za1000000a')
 
     import filenames
     m.export_mesh_file(filenames.mesh_filename)
-    

inundation/validation/LWRU2/lwru2.py

@@ -71 +71 @@
             File_boundary, Transmissive_Momentum_Set_Stage_boundary
 from pyvolution.mesh_factory import rectangular
-from pyvolution.pmesh2domain import pmesh_to_domain_instance 
+from pyvolution.pmesh2domain import pmesh_to_domain_instance
 from Numeric import array, zeros, Float, allclose
 import filenames
@@ -81 +81 @@
 
 try:
-    raise Exception
+    #raise Exception
     fid = open(picklefile)
     print 'Read pickled domain'

inundation/zeus/anuga-workspace.zwi

@@ -7 +7 @@
     <project name="pyvolution">pyvolution.zpi</project>
     <project name="steve">steve.zpi</project>
+    <project name="Validation">Validation.zpi</project>
 </workspace>