Changeset 2446

Timestamp:

Feb 24, 2006, 4:09:02 PM (19 years ago)

Author:

nicholas

Message:

Attempted loop_buildings.py did not work.

Location:

development/momentum_sink

Files:

• CCS.py (modified) (7 diffs)
• PCS.py (modified) (2 diffs)
• create_buildings.py (modified) (6 diffs)
• loop_buildings.py (added)
• loop_friction.py (modified) (1 diff)
• run_buildings.py (modified) (3 diffs)

development/momentum_sink/CCS.py

@@ -12 +12 @@
 from numerical_tools import norm, corr, err  
 
-swwfile_B = project.outputdir + sep  + 'Buildings_3260.sww'
 
 
-gauge_depth = Numeric.arrayrange(0, 1.5*WidtH, 20)
+
+gauge_depth = Numeric.arrayrange(46, 1.5*WidtH, 25)
 print len(gauge_depth),"gauge_depth length"
 gauge_breadth = 100
@@ -24 +24 @@
     gauge_locations.append(gauge_location)
 
+
+swwfile_B = project.outputdir + sep  + 'Straight_offset_11.sww'
 
 #Read model output
@@ -38 +40 @@
 least_norm=[1]
 norm_OK=[]
-Friction = Numeric.arrayrange(1,51,1)
+Friction = Numeric.arrayrange(1,80,1)
 for fric in Friction:
     swwfile_F = project.outputdir + sep  + 'friction_n=' + repr(fric) + '_3135.sww'
@@ -47 +49 @@
                       use_cache = True)
 
-    t=500
+    t=1000 # specifies times to look at cross sections
     
     stages_B = []
@@ -54 +56 @@
 
         # mines data from specified time values in sww file
+        # Building file
         wB = f_B(t, point_id = i)[0]
         zB = f_B(t, point_id = i)[1]
         uhB = f_B(t, point_id = i)[2]
         vhB = f_B(t, point_id = i)[3]
-
+        # Friction file
         wF = f_F(t, point_id = i)[0]        
         uhF = f_B(t, point_id = i)[2]
         vhF = f_B(t, point_id = i)[3]
+
         #m = sqrt(uh*uh + vh*vh)   #Absolute momentum
         velB = sqrt(uhB*uhB + vhB*vhB) / (wB + 1.e-30) #Absolute velocity
         velF = sqrt(uhF*uhF + vhF*vhF) / (wF + 1.e-30) #Absolute velocity                   
-
         stages_B.append(wB)
         stages_F.append(wF)
@@ -73 +76 @@
     Ar_F = Numeric.array(stages_F)
     diff = abs(Ar_B - Ar_F)
-    norm_F = err(Ar_B, Ar_F, 2, relative = True)
+    norm_F = err(Ar_B, Ar_F, 2, relative = True)# 2nd norm (RMS) test
      
     print " "
     print norm_F, "norm N_F"
     print " "
-    G = corr(Ar_F, Ar_B)
     norm_OK.append(norm_F)
 
@@ -91 +93 @@
 print " "
 print "_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_"
-print fric_OK, " <- This Manning's n is OK"
+print fric_OK, " <- This Manning's n is best fit from 2nd norm test"
 print "_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_"
 print least_norm, "2nd normal of Friction vs Buildings  "

development/momentum_sink/PCS.py

@@ -11 +11 @@
 from pylab import *
 
-fric = 40    
+fric = 25    
 #swwfile = project.newoutputname + '.sww'
 #swwfile = project.outputname
 
 #swwfile_B = project.outputdir + sep  + 'Buildings_3790.sww'
-swwfile_B = project.outputdir + sep  + 'Buildings_3662.sww'
+swwfile_B = project.outputdir + sep  + 'Straight_offset.sww'
 #swwfile_B = project.outputdir + sep  + 'friction_n=10_3135.sww'
-swwfile_F = project.outputdir + sep  + 'friction_n=40_3135.sww'
+swwfile_F = project.outputdir + sep  + 'friction_n=51_3135.sww'
 
-gauge_depth = Numeric.arrayrange(0, 1.5*WidtH, 5)
+gauge_depth = Numeric.arrayrange(46, 1.5*WidtH, 25) # Random special
+#gauge_depth = Numeric.arrayrange(0, 1.5*WidtH, 25)
 gauge_breadth = 100
 gauge_locations = []
@@ -45 +46 @@
 
 
-T = Numeric.arrayrange(0,1000,20)
+T = Numeric.arrayrange(0,1000,10)
 #T = [ 20, 50, 100, 150, 200 ]
 

development/momentum_sink/create_buildings.py

@@ -2 +2 @@
 import Numeric
 import math
+import random
 # add inundation dir to your pythonpath
 from pmesh.mesh import Mesh
 from coordinate_transforms.geo_reference import Geo_reference
 
+depth = 11 # depth of building side to oncoming wave
+breadth = 11 # breadth of building, width of building to oncoming wave
+#from loop_buildings import depth
 
 WidtH = 200 # width of boudary in metres
-depth = 15 # depth of building side to oncoming wave
-breadth = 15 # breadth of building
+#breadth = depth
+
 print "building footprint"
 print depth * breadth , "m^2"
@@ -15 +19 @@
 BL = block**0.5
 
-porosity = breadth/BL
+porosity = breadth/BL                           
 print porosity, " Building porosity"
+
+#random.uniform(0,1) #proper implemantation of random generator
 
 def create_mesh(maximum_triangle_area,
@@ -25 +31 @@
     triangles in the mesh to the mesh file name.
     """
+
+    
     # create a mesh instance of class Mesh
     m = Mesh()
+
     # Boundary of problem
-    
-    #W = WidtH/8
-    #L = W
     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,2], 'wave':[3], 'back':[1]})
 
@@ -41 +46 @@
     #Th = (45 *(3.142/180)) # sets an initial rotation      
     Th = 0
-    ForDep = (0.2*WidtH) + (BL/2)
+    ForDep = (0.2*WidtH) + (BL-whs)
     RearDep = 1.2*WidtH 
         
@@ -61 +66 @@
             lh4 = (+lhs) * math.cos(Th) - (-whs) * math.sin(Th)
             polygon = [[D+wh1,B+lh1],[D+wh2,B+lh2],[D+wh3,B+lh3],[D+wh4,B+lh4]]
-            #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.        
-            Th = Th + (37.3 *(3.142/180)) # keeps rotating individual buildings.
+            m.add_hole_from_polygon(polygon, tags={'wall':[0,1,2,3]})        
+            #Th = Th + (90 *(3.14159/180)) # keeps rotating individual buildings.
 
-            #print polygon
     m.generate_mesh(maximum_triangle_area=maximum_triangle_area)
     triangle_count = m.get_triangle_count()
@@ -73 +76 @@
     else:
         if triangles_in_name is True:
-            mesh_file = mesh_file[:-4] + '_' + str(triangle_count) \
+            mesh_file = mesh_file[:-4] + '_Ran_' + str(depth) + '_' + str(triangle_count) \
                         + mesh_file[-4:]
         m.export_mesh_file(mesh_file)

development/momentum_sink/loop_friction.py

@@ -24 +24 @@
 from pmesh.mesh import importMeshFromFile
 
-Friction = Numeric.arrayrange(1,50,1)
+Friction = Numeric.arrayrange(50,80,1)
 for fric in Friction:
     meshname = project_friction.meshname

development/momentum_sink/run_buildings.py

@@ -15 +15 @@
 
 
-from create_buildings import create_mesh
+from create_buildings import create_mesh, depth
 #from building_generator import create_mesh
 
@@ -46 +46 @@
 
 
-domain.set_name(project.basename + '_%d' %triagle_count)
+domain.set_name(project.basename + '_Br=%d_%d' %(depth, triagle_count))
 domain.set_datadir(project.outputdir)
 domain.store = True
@@ -75 +75 @@
 t0 = time.time()
 
-for t in domain.evolve(yieldstep = 10, finaltime = 500):
+for t in domain.evolve(yieldstep = 10, finaltime = 1000):
     domain.write_time()