Changeset 2417 for development/momentum_sink/PCS.py

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)

File:

• development/momentum_sink/PCS.py (modified) (3 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"