Changeset 9292 for trunk/anuga_work/development/gareth/tests/runup_sinusoid

Timestamp:

Aug 12, 2014, 12:10:51 PM (11 years ago)

Author:

davies

Message:

Adding local extrapolation and flux updating

File:

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

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

@@ -20 +20 @@
 domain.set_name('runup_sinusoid_v2')                         # Output to file runup.sww
 #domain.set_store_centroids(True)
-domain.set_flow_algorithm('DE1')
+domain.set_flow_algorithm('DE0')
 domain.set_quantities_to_be_stored({'stage': 2, 'xmomentum': 2, 
          'ymomentum': 2, 'elevation': 2})
@@ -42 +42 @@
 bumpiness=50. # Higher = shorter wavelength oscillations in topography
 tstep=0.2
-lasttime=90.
+lasttime=30.
 
 #domain.minimum_allowed_height=domain.minimum_allowed_height*scale_me # Seems needed to make the algorithms behave
@@ -110 +110 @@
     print 'Their difference is:', vol-flux_integral
     #print 'Volume less flux int', sum(dd_raw*domain.areas) - domain.boundary_flux_integral
-    
+    print 'Update_next_flux range:'
+    print domain.update_extrapolation.min(), domain.update_extrapolation.max()
+    print 'Flux_update_frequency range:'
+    print domain.flux_update_frequency.min(), domain.flux_update_frequency.max()
 
 vel_final_inds=(vv>1.0e-01).nonzero()