Changeset 2417

Timestamp:

Feb 16, 2006, 3:41:21 PM (19 years ago)

Author:

nicholas

Message:

Creation of looped friction to auto generate multiple friction blocks. CCS.py created as modified copy of PCS.py ( plots cross section at various times of buildings and friction blocks)

Location:

development/momentum_sink

Files:

• CCS.py (added)
• PCS.py (modified) (3 diffs)
• create_buildings.py (modified) (3 diffs)
• loop_friction.py (added)
• project_friction.py (modified) (1 diff)
• run_friction.py (modified) (2 diffs)

development/momentum_sink/PCS.py

@@ -7 +7 @@
 import project
 from pyvolution.util import file_function
-#from pyvolution.coordinate_transforms.redfearn import degminsec2decimal_degrees, redfearn
-#from coordinate_transforms.redfearn import degminsec2decimal_degrees, redfearn
+from create_buildings import porosity, WidtH
+# from run_friction import fric
 from pylab import *
-#from compare_sww import gauge_locations
 
-    
+fric = 1    
 #swwfile = project.newoutputname + '.sww'
 #swwfile = project.outputname
-swwfile_B = project.outputdir + sep  + 'Buildings_3790.sww'
-swwfile_F = project.outputdir + sep  + 'Buildings_3662.sww'
-#swwfile_F = project.outputdir + sep  + 'friction_3135.sww'
 
-gauge_depth = Numeric.arrayrange(0, 700, 10)
+#swwfile_B = project.outputdir + sep  + 'Buildings_3790.sww'
+#swwfile_B = project.outputdir + sep  + 'Buildings_3662.sww'
+swwfile_B = project.outputdir + sep  + 'friction_n=10_3135.sww'
+swwfile_F = project.outputdir + sep  + 'friction_n=20_3135.sww'
+
+gauge_depth = Numeric.arrayrange(0, 2.5*WidtH, 5)
 gauge_breadth = 100
 gauge_locations = []
@@ -55 +56 @@
     velocity = []
     for i, GL in enumerate(gauge_depth):
+
+        # mines data from specified time values in sww file
         wB = f_B(t, point_id = i)[0]
         zB = f_B(t, point_id = i)[1]
@@ -82 +85 @@
 
     ion()
-    hold(False)
+    hold(True)
+    subplot(211)
+    
     plot(gauge_depth, stages_B, '-b',gauge_depth,stages_F,'-r')
-    title('Time_%d_ vs water depth' %t)
+    title('Time_%d_ vs water depth (n=%f, BP=%f)' %(t, fric, porosity))
     xlabel('gauge length(m)')
     ylabel('water depths (m)')
-    savefig('Comp_Time%d_B_vs_F' %t)
-    #plot(gauge_depth, diff, 'g')
-    #title('differences')
-    # raw_input('Next')
+    #savefig('Comp_Time%d_B_vs_F' %t)
+    hold(False)
+#savefig('Final_comp_n=%s' %fric)
+subplot(212)
+subplots_adjust(hspace = 0.4)
+plot(gauge_depth, diff, 'g')
+title('Final Differences in water levels')
+xlabel('gauge length(m)')
+ylabel('water differences (m)')
+savefig('Final_Water_diff_n=%s' %fric)
 print "finished"
 

development/momentum_sink/create_buildings.py

@@ -7 +7 @@
 
 
+WidtH = 200 # width of boudary in metres
+depth = 20 # depth of building side to oncoming wave
+breadth = 20 # breadth of building
+print "building footprint"
+print depth * breadth , "m^2"
+block = 625
+BL = block**0.5
+
+porosity = breadth/BL
+print porosity, " Building porosity"
 
 def create_mesh(maximum_triangle_area,
@@ -18 +28 @@
     m = Mesh()
     # Boundary of problem
-    WidtH = 200 # width of boudary in metres
+    
     #W = WidtH/8
     #L = W
@@ -26 +36 @@
 
     # inner polygons => building boundaries
-    depth = 4 # depth of building side to oncoming wave
+    
     wh = depth/2 
-    breadth = 4 # breadth of building
+    
     lh = breadth/2 
           
-    print "building footprint"
-    print depth * breadth , "m^2"
-    block = 625
-    BL = block**0.5
     ForDep = (0.2*WidtH) + (BL/2)
     RearDep = 1.2*WidtH 
-    porosity = breadth/BL
-    print porosity, " Building porosity"
         
     Breadths = Numeric.arrayrange( (BL/2), WidtH, (BL))

development/momentum_sink/project_friction.py

@@ -12 +12 @@
 
 scenario_dir_name = 'momentum_sink_2005'
-basename = 'friction'
+basename = 'friction_n='
 
 if sys.platform == 'win32':

development/momentum_sink/run_friction.py

@@ -4 +4 @@
 #           #name does not have the extension
 
+
+# set friction value for zone
+fric=10
+
+# fwojhrt
 import time
 
@@ -60 +65 @@
 domain.set_quantity('friction', 0.01)
 domain.set_quantity('stage', 0)
-domain.set_region(Set_region('mound', 'friction', 10)) #, location='unique vertices'))
+domain.set_region(Set_region('mound', 'friction', fric)) #, location='unique vertices'))
 #Setup Boundary Conditions
 print domain.get_boundary_tags()