Changeset 5079

Timestamp:

Feb 27, 2008, 9:12:46 AM (17 years ago)

Author:

ole

Message:

Added some statistics to run_dam.

File:

• anuga_core/source/anuga/debug/bad_time_step/run_dam.py (modified) (6 diffs)

anuga_core/source/anuga/debug/bad_time_step/run_dam.py

@@ -27 +27 @@
 # Related major packages
 from anuga.shallow_water import Domain, Reflective_boundary, \
-                            Dirichlet_boundary,  Time_boundary, \
-                            File_boundary, \
-                            Transmissive_Momentum_Set_Stage_boundary
+     Dirichlet_boundary,  Time_boundary, \
+     File_boundary, \
+     Transmissive_Momentum_Set_Stage_boundary
 from anuga.fit_interpolate.interpolate import interpolate_sww2csv
 from anuga.abstract_2d_finite_volumes.util import start_screen_catcher, \
@@ -59 +59 @@
         finaltime = 15.
         maximum_triangle_area=0.01
+        #maximum_triangle_area=0.1 #OK
     else:
         yieldstep = 0.02
@@ -69 +70 @@
                                    outputdir_name=outputdir_name,
                                    home='.')
-    print "The output dir is", pro_instance.outputdir
+    print 'The output dir is', pro_instance.outputdir
     copy_code_files(pro_instance.outputdir,__file__,
                     dirname(project.__file__) \
@@ -90 +91 @@
 
     print 'USER:    ', pro_instance.user
+
+    
     #-------------------------------------------------------------------------
     # Create the triangular mesh
@@ -132 +135 @@
     # Create boundary function from timeseries provided in file
     function = file_function(project.boundary_file, domain, verbose=True)
+
+    #print dir(function)
+    #print function.time, function.quantities_range
+    #import sys; sys.exit()
+    
     Bts = Transmissive_Momentum_Set_Stage_boundary(domain, function)
 
@@ -145 +153 @@
 
     for t in domain.evolve(yieldstep, finaltime):
-    
-        domain.write_time()
-        print 'That took %.2f seconds' %(time.time()-t0)
-        print 'finished'
+        print domain.timestepping_statistics(track_speeds=False)
+        print domain.boundary_statistics(tags=['wave'])
+        
+    print 'That took %.2f seconds' %(time.time()-t0)
+    print 'finished'
 
     points = [[2.8,0.225],  #-1.8m from SWL