Changeset 7492 for anuga_validation/okushiri_2005

Timestamp:

Sep 7, 2009, 4:19:23 PM (16 years ago)

Author:

steve

Message:

Changed the order of update_timestep and compute_forcing_terms.
Can now update the parameter domain.flux_timestep within
the computation of the forcing terms, so as to control stability of
the evolver.

Hopefully this will help in reducing the oscillations seen within the
culvert computation.

File:

• anuga_validation/okushiri_2005/run_okushiri.py (modified) (3 diffs)

anuga_validation/okushiri_2005/run_okushiri.py

@@ -22 +22 @@
 
 # Module imports
-from anuga.shallow_water import Domain
-from anuga.shallow_water import Reflective_boundary
-from anuga.shallow_water import Transmissive_momentum_set_stage_boundary
-from anuga.abstract_2d_finite_volumes.util import file_function
+from anuga.interface import Domain
+from anuga.interface import Reflective_boundary
+from anuga.interface import Time_boundary
+from anuga.interface import Transmissive_momentum_set_stage_boundary
+from anuga.interface import Transmissive_n_momentum_zero_t_momentum_set_stage_boundary
+from anuga.interface import file_function
 
 import project
@@ -82 +84 @@
 triangle_id = domain.get_triangle_containing_point([4.521, 1.196])
 # This should get triangle id 32833 with centroid (4.5244, 1.1972)
+bdry_id = domain.get_triangle_containing_point([0.000, 1.696])
+print 'bdry_id = ',bdry_id
 
 
@@ -92 +96 @@
 for t in domain.evolve(yieldstep = 0.05, finaltime = 22.5):
     domain.write_time()
+    print  function(domain.get_time())[0],' bdry_d = ',bdry_id,' ',\
+          domain.get_conserved_quantities(bdry_id, edge=0)[0],' ',\
+          domain.get_conserved_quantities(bdry_id, edge=1)[0],' ',\
+          domain.get_conserved_quantities(bdry_id, edge=2)[0]    
+    
 
 print 'That took %.2f seconds' %(time.time()-t0)