Changeset 2415 for development/momentum_sink/create_buildings.py

Timestamp:

Feb 15, 2006, 4:47:01 PM (19 years ago)

Author:

nicholas

Message:

 

File:

• development/momentum_sink/create_buildings.py (modified) (4 diffs)

development/momentum_sink/create_buildings.py

@@ -5 +5 @@
 from pmesh.mesh import Mesh
 from coordinate_transforms.geo_reference import Geo_reference
+
 
 
@@ -16 +17 @@
     # create a mesh instance of class Mesh
     m = Mesh()
-
     # Boundary of problem
     WidtH = 200 # width of boudary in metres
@@ -23 +23 @@
     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]})
+    m.add_region_from_polygon(outer_polygon, tags={'wall':[0,2], 'wave':[3], 'back':[1]})
 
     # inner polygons => building boundaries
-    depth = 20 # depth of building side to oncoming wave
+    depth = 4 # depth of building side to oncoming wave
     wh = depth/2 
-    breadth = 20 # breadth of building
+    breadth = 4 # breadth of building
     lh = breadth/2 
-      
+          
     print "building footprint"
     print depth * breadth , "m^2"
@@ -39 +39 @@
     porosity = breadth/BL
     print porosity, " Building porosity"
-    
+        
     Breadths = Numeric.arrayrange( (BL/2), WidtH, (BL))
     print Breadths, "Breadths"