Changeset 8446 for trunk/anuga_work/development/gareth/tests/channel_floodplain/channel_floodplain1.py

Timestamp:

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

Author:

davies

Message:

bal_dev: Updates

File:

• trunk/anuga_work/development/gareth/tests/channel_floodplain/channel_floodplain1.py (modified) (4 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 '$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$'