Changeset 2341

Timestamp:

Feb 6, 2006, 4:26:06 PM (19 years ago)

Author:

nicholas

Message:

Added friction comparision for building, Fiddled with loops fro create buildings

Location:

development/momentum_sink

Files:

• create_buildings.py (modified) (1 diff)
• friction_block.py (added)
• run_buildings.py (modified) (3 diffs)
• run_friction.py (added)

development/momentum_sink/create_buildings.py

@@ -21 +21 @@
     #W = WidtH/8
     #L = W
-    outer_polygon = [[0,0],[WidtH,0],[WidtH,WidtH],[0,WidtH]]
+    outer_polygon = [[0,0],[2*WidtH,0],[2*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 = 15
-    wh = depth/2 # depth of building side to oncoming wave
-    breadth = 22
-    lh = breadth/2 # breadth of building
+    depth = 15 # depth of building side to oncoming wave
+    wh = depth/2 
+    breadth = 15 # breadth of building
+    lh = breadth/2 
     print "building footprint"
     print depth * breadth , "m^2"
+    block = 625
+    BL = block**0.5
+    print BL
     #[12.5,37.5,62.5,87.5,112.5,137.5,162.5,187.5]
-    
-    for W in [12.5,37.5,62.5,87.5,112.5,137.5,162.5,187.5]:
-        for L in [12.5,37.5,62.5,87.5,112.5,137.5,162.5,187.5]:
-            polygon = [[W-wh,L-lh],[W+wh,L-lh],[W+wh,L+lh],[W-wh,L+lh]]
+    D = WidtH/2 + BL/2
+    B = BL/2
+    for D in [112.5,137.5,162.5,187.5]:
+        for B in [12.5,37.5,62.5,87.5,112.5,137.5,162.5,187.5]:
+            polygon = [[D-wh,B-lh],[D+wh,B-lh],[D+wh,B+lh],[D-wh,B+lh]]
             m.add_hole_from_polygon(polygon, tags={'wall':[0,1,2,3]})# Adds holes with reflective boundaries.        
             #print polygon
+            
+            print B
 
-    m.add_circle([20,70], 5, hole=True, tag='wall')
+        print D
 
     m.generate_mesh(maximum_triangle_area=maximum_triangle_area)

development/momentum_sink/run_buildings.py

@@ -16 +16 @@
 
 from create_buildings import create_mesh
+#from building_generator import create_mesh
 
 from pmesh.mesh import importMeshFromFile
@@ -24 +25 @@
 t0 = time.time()
 
-meshname, triagle_count = cache(create_mesh,(100),
+meshname, triagle_count = cache(create_mesh,(10),
                                 {'mesh_file':meshname,
                                  'triangles_in_name':True}
@@ -69 +70 @@
 
 domain.set_boundary( {'wall': Br,'wave': Bd, 'back': Bt, 'exterior':Bd} ) 
-
+domain.set_region(Set_region('mound', 'friction', 3, location='unique vertices'))
 #Evolve
 t0 = time.time()