Changeset 8446 for trunk/anuga_work/development/gareth/tests/channel_floodplain

Timestamp:

Jun 14, 2012, 12:31:32 AM (12 years ago)

Author:

davies

Message:

bal_dev: Updates

Location:

trunk/anuga_work/development/gareth/tests/channel_floodplain

Files:

• channel_floodplain1.py (modified) (4 diffs)
• plotme.py (modified) (2 diffs)

trunk/anuga_work/development/gareth/tests/channel_floodplain/channel_floodplain1.py

@@ -10 +10 @@
 import numpy
 from anuga.structures.inlet_operator import Inlet_operator
-from anuga import *
+#from anuga import *
 #from swb_domain import domain
 #from anuga import *
@@ -82 +82 @@
 
 
-domain.set_name('channel_floodplain1_bal_dev_lowvisc') # Output name
+domain.set_name('channel_floodplain1') # Output name
 #domain.set_store_vertices_uniquely(True)
 #domain.use_edge_limiter=False
@@ -142 +142 @@
 
 # Define inlet operator 
-flow_in_yval=100.0
+flow_in_yval=10.0
 if True:
     line1 = [ [floodplain_width/2. - chan_width/2., flow_in_yval],\
@@ -245 +245 @@
     if( numpy.floor(t/100.) == t/100. ):
         print '#### COMPUTING FLOW THROUGH CROSS-SECTIONS########'
+        s0 = domain.get_flow_through_cross_section([[0., 10.0], [floodplain_width, 10.]])
         s1 = domain.get_flow_through_cross_section([[0., floodplain_length-300.0], [floodplain_width, floodplain_length-300.0]])
         s2 = domain.get_flow_through_cross_section([[0., floodplain_length-1.0], [floodplain_width, floodplain_length-1.0]])
         
-        print 'Cross sectional flows: ', s1, s2 
+        print 'Cross sectional flows: ', s0, s1, s2 
         print '$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$'
 

trunk/anuga_work/development/gareth/tests/channel_floodplain/plotme.py

@@ -4 +4 @@
 
 # Time-index to plot outputs from
-index=53
+index=60
 
 #p2 = util.get_output('channel_floodplain1_bal_dev.sww', minimum_allowed_height=0.01)
@@ -10 +10 @@
 #p=util.get_centroids(p2, velocity_extrapolation=True)
 
-p2 = util.get_output('channel_floodplain1_bal_dev_lowvisc.sww', 0.01)
+#p2 = util.get_output('channel_floodplain1_bal_dev_lowvisc.sww', 0.01)
+p2 = util.get_output('channel_floodplain1.sww', 0.01)
 p=util.get_centroids(p2, velocity_extrapolation=True)