Changeset 2395

Timestamp:

Feb 14, 2006, 1:43:12 PM (19 years ago)

Author:

nicholas

Message:

Updated compare_sww.py. Increased distance of rear boundary and changed to dirichlet

Location:

development/momentum_sink

Files:

• compare_sww.py (modified) (2 diffs)
• create_buildings.py (modified) (2 diffs)
• friction_block.py (modified) (5 diffs)
• process_gauges.py (added)
• run_buildings.py (modified) (1 diff)
• run_friction.py (modified) (4 diffs)
• run_rotate.py (modified) (1 diff)

development/momentum_sink/compare_sww.py

@@ -5 +5 @@
 import project
 from pyvolution.util import file_function
+import Numeric
+#from friction_block #import test
+
+#print "test", test
+#import sys
+#sys.exit() 
 
 
 #swwfile = project.outputname + '.sww'
-swwfile = project.outputdir + sep  + 'friction_12330.sww'
-gauge_locations = [[10,10],[20,20]]
+#swwfile = project.outputdir + sep  + 'FrictionBlock_0.01_3.sww'
+swwfile = project.outputdir + sep  + 'Buildings_3662.sww'
+gauge_depth = Numeric.arrayrange(0, 700, 20)
+gauge_breadth = 100
+gauge_locations = []
 
-quantities = ['stage','elevation']
+for GD in gauge_depth:
+    gauge_location = [GD,gauge_breadth]
+    gauge_locations.append(gauge_location)
+
+quantities = ['stage','elevation','xmomentum']
 gauges = file_function(swwfile,
                   quantities = quantities,
@@ -18 +31 @@
                   use_cache = True)
 
-times = [0,1,2,3,4]
+times = [150,500]
 
 for gauge, _ in enumerate(gauge_locations):

development/momentum_sink/create_buildings.py

@@ -18 +18 @@
 
     # Boundary of problem
-    WidtH = 500 # width of boudary in metres
+    WidtH = 200 # width of boudary in metres
     #W = WidtH/8
     #L = W
-    outer_polygon = [[0,0],[1.5*WidtH,0],[1.5*WidtH,WidtH],[0,WidtH]]
+    outer_polygon = [[0,0],[5*WidtH,0],[5*WidtH,WidtH],[0,WidtH]]
     print outer_polygon
     m.add_region_from_polygon(outer_polygon, tags={'wall':[0,1,2], 'wave':[3]})
 
     # inner polygons => building boundaries
-    depth = 18 # depth of building side to oncoming wave
+    depth = 20 # depth of building side to oncoming wave
     wh = depth/2 
-    breadth = 18 # breadth of building
+    breadth = 20 # breadth of building
     lh = breadth/2 
       
@@ -46 +46 @@
 
     for D in Depths:
-        Breadths = Breadths + BL/2
+        #Breadths = Breadths + BL/2 #Used to offset buildings
         for B in Breadths:
             polygon = [[D-wh,B-lh],[D+wh,B-lh],[D+wh,B+lh],[D-wh,B+lh]]

development/momentum_sink/friction_block.py

@@ -6 +6 @@
 from coordinate_transforms.geo_reference import Geo_reference
 
-
+#test = 300
 def create_mesh(maximum_triangle_area,
                 mesh_file=None,
@@ -16 +16 @@
     # create a mesh instance of class Mesh
     m = Mesh()
-
+    
+    #print "test", test  
     # Boundary of problem
     WidtH = 200 # width of boudary in metres
@@ -28 +29 @@
     #Boundary
     dict = {}
-    dict['points'] = [[0,0],[0.2*WidtH,0],[1.2*WidtH,0],[3.5*WidtH,0],[3.5*WidtH,WidtH] \
+    dict['points'] = [[0,0],[0.2*WidtH,0],[1.2*WidtH,0],[5*WidtH,0],[5*WidtH,WidtH] \
                       ,[1.2*WidtH,WidtH],[0.2*WidtH,WidtH],[0,WidtH]]
     
@@ -37 +38 @@
                             'wall',                           
                             'wall',
-                            'wall',
+                            'back',
                             'wall',
                             'wall',
@@ -50 +51 @@
     m.generate_mesh(maximum_triangle_area=maximum_triangle_area)
     triangle_count = m.get_triangle_count()
+  
+
     
     if mesh_file is None:    
-        return m, triangle_count
+        return m, triangle_count, Width
     else:
         if triangles_in_name is True:

development/momentum_sink/run_buildings.py

@@ -65 +65 @@
 Br = Reflective_boundary(domain)
 Bt = Transmissive_boundary(domain)
-Bd = Dirichlet_boundary([10,0,0])
+Bdw = Dirichlet_boundary([36,0,0])
+Bdb = Dirichlet_boundary([0,0,0])
 Bw = Time_boundary(domain=domain,
                    f=lambda t: [(60<t<660)*4, 0, 0])
 
-domain.set_boundary( {'wall': Br,'wave': Bd, 'back': Bt, 'exterior':Bd} ) 
+domain.set_boundary( {'wall': Br,'wave': Bdw, 'back': Bdb, 'exterior':Bdw} ) 
 
 #Evolve
 t0 = time.time()
 
-for t in domain.evolve(yieldstep = 1, finaltime = 20):
+for t in domain.evolve(yieldstep = 10, finaltime = 1000):
     domain.write_time()      
 

development/momentum_sink/run_friction.py

@@ -26 +26 @@
 t0 = time.time()
 
+#meshname, triagle_count =create_mesh(100,mesh_file=meshname,triangles_in_name=True)
+
 meshname, triagle_count = cache(create_mesh,(100),
                                 {'mesh_file':meshname,
@@ -50 +52 @@
 domain.set_datadir(project_friction.outputdir)
 domain.store = True
-domain.quantities_to_be_stored = ['stage']
+domain.quantities_to_be_stored = ['stage','xmomentum','ymomentum']
 
 print 'Number of triangles = ', len(domain)
@@ -58 +60 @@
 domain.set_quantity('friction', 0.01)
 domain.set_quantity('stage', 0)
-domain.set_region(Set_region('mound', 'friction', 0.8)) #, location='unique vertices'))
+domain.set_region(Set_region('mound', 'friction', 3)) #, location='unique vertices'))
 #Setup Boundary Conditions
 print domain.get_boundary_tags()
@@ -66 +68 @@
 Br = Reflective_boundary(domain)
 Bt = Transmissive_boundary(domain)
-Bd = Dirichlet_boundary([36,0,0])
+Bdw = Dirichlet_boundary([36,0,0])
+Bdb = Dirichlet_boundary([0,0,0])
 Bw = Time_boundary(domain=domain,
                    f=lambda t: [(60<t<660)*4, 0, 0])
 
-domain.set_boundary( {'wall': Br,'wave': Bd, 'back': Bt, 'exterior':Bd} ) 
+domain.set_boundary( {'wall': Br,'wave': Bdw, 'back': Bdb, 'exterior':Bdw} ) 
 
 #Evolve
 t0 = time.time()
 
-for t in domain.evolve(yieldstep = 1, finaltime = 80):
+for t in domain.evolve(yieldstep = 4, finaltime = 400):
     domain.write_time()      
 

development/momentum_sink/run_rotate.py

@@ -74 +74 @@
 t0 = time.time()
 
-for t in domain.evolve(yieldstep = 1, finaltime = 20):
+for t in domain.evolve(yieldstep = 1, finaltime = 60):
     domain.write_time()