Changeset 9257 for trunk/anuga_work/development/gareth/tests/runup_sinusoid/runup_sinusoid.py

Timestamp:

Jul 10, 2014, 5:42:24 PM (11 years ago)

Author:

davies

Message:

Cutting a boundary_flux_integral operator into the code

File:

• trunk/anuga_work/development/gareth/tests/runup_sinusoid/runup_sinusoid.py (modified) (5 diffs)

trunk/anuga_work/development/gareth/tests/runup_sinusoid/runup_sinusoid.py

@@ -9 +9 @@
 
 from math import sin, pi, exp
-#from anuga.shallow_water.shallow_water_domain import Domain as Domain
-#from anuga.shallow_water_balanced2.swb2_domain import Domain as Domain
-#path.append('/home/gareth/storage/anuga_clean/anuga_jan12/trunk/anuga_work/development/gareth/balanced_basic')
-#from swb2_domain import *
-#from balanced_basic import *
-#from balanced_dev import *
 from anuga import Domain as Domain
-#from bal_and import *
-#from anuga_tsunami import *
+from anuga.operators.boundary_flux_integral_operator import boundary_flux_integral_operator
 
 #---------
@@ -27 +20 @@
 domain.set_name('runup_sinusoid_v2')                         # Output to file runup.sww
 #domain.set_store_centroids(True)
-#domain.set_timestepping_method('euler')
-#domain.set_flow_algorithm('DE1')
-#domain.set_flow_algorithm('tsunami')
 domain.set_flow_algorithm('DE0')
 domain.set_quantities_to_be_stored({'stage': 2, 'xmomentum': 2, 
          'ymomentum': 2, 'elevation': 2})
-#domain.minimum_storable_height=1.0e-03
 domain.set_store_vertices_uniquely()
 #------------------
@@ -70 +59 @@
 
 print '#-#', domain.minimum_allowed_height
-#domain.minimum_allowed_height=0.001
-#def frict_change(x,y):
-#       return 0.2*(x>0.5)+0.1*(x<=0.5)
-#
-#domain.set_quantity('friction',frict_change)
 
 domain.set_quantity('stage', stagefun)              # Constant negative initial stage
@@ -108 +92 @@
 #Evolve the system through time
 #------------------------------
-
 for t in domain.evolve(yieldstep=tstep,finaltime=lasttime):
 #for t in domain.evolve(yieldstep=0.2,finaltime=60.):
@@ -121 +104 @@
     vv = vv*(dd>1.0e-03)
     print 'Peak velocity is: ', vv.max(), vv.argmax(), dd[vv.argmax()]
-    print 'Volume is', sum(dd_raw*domain.areas)    
+    vol=sum(dd_raw*domain.areas)
+    print 'Volume is', vol    
+    flux_integral=domain.get_boundary_flux_integral()
+    print 'Flux integral is: ', flux_integral
+    print 'Their difference is:', vol-flux_integral
     #print 'Volume less flux int', sum(dd_raw*domain.areas) - domain.boundary_flux_integral
-
+    
 
 vel_final_inds=(vv>1.0e-01).nonzero()