Changeset 4173 for anuga_core/source/anuga/examples

Timestamp:

Jan 12, 2007, 10:18:44 AM (18 years ago)

Author:

duncan

Message:

deleting old files, updating others.

Location:

anuga_core/source/anuga/examples

Files:

• run_gippsland.py (deleted)
• run_gong.py (deleted)
• run_tsh.py (modified) (6 diffs)

anuga_core/source/anuga/examples/run_tsh.py

@@ -6 +6 @@
 
 A example of running this program is;
-python run_tsh.py visualise hill.tsh 0.05 1
+python run_tsh.py n hill.tsh 0.05 1
 """
 
@@ -13 +13 @@
 #
 
+from Numeric import array
+import time
 import sys
 from os import sep, path
-sys.path.append('..'+sep+'pyvolution')
 
-from shallow_water import Domain, Reflective_boundary, Dirichlet_boundary,\
-     Transmissive_boundary, Time_boundary
-from region import Add_value_to_region, Set_region
-from mesh_factory import rectangular
-from anuga.pyvolution.pmesh2domain import pmesh_to_domain_instance
+from anuga.shallow_water import Domain, Reflective_boundary, \
+     Dirichlet_boundary, Transmissive_boundary, Time_boundary
+from anuga.abstract_2d_finite_volumes.region import Add_value_to_region, \
+     Set_region
+from anuga.visualiser import RealtimeVisualiser
 
-from Numeric import array
 
-import time
 
 #from anuga.config import default_datadir 
@@ -50 +49 @@
     else:    
         visualise = True
+       
     filename = sys.argv[2]
     yieldstep = float(sys.argv[3])
@@ -55 +55 @@
     
     print 'Creating domain from', filename
-    domain = pmesh_to_domain_instance(filename, Domain)
+    domain = Domain(filename)
+
+    # check if the visualiser will work
+    try:
+        xx = RealtimeVisualiser(domain)
+    except:
+        print "Warning: Error in RealtimeVisualiser. Could not visualise."
+        visualise = False
+    
     print "Number of triangles = ", len(domain)
     print "domain.geo_reference",domain.geo_reference 
@@ -85 +93 @@
     inflow_stage = 10.0
 
-    ######NEW
+    
     domain.set_quantity('friction', manning)
 
@@ -93 +101 @@
     #                        domain.quantities['elevation'].vertex_values)
     #domain.set_quantity('stage', 0.0)
-    ######NEW
+    
 
     ######################