Changeset 7702

Timestamp:

Apr 29, 2010, 3:50:10 PM (15 years ago)

Author:

steve

Message:

Using small file

File:

• anuga_core/source/anuga_parallel/run_parallel_sw_merimbula.py (modified) (5 diffs)

anuga_core/source/anuga_parallel/run_parallel_sw_merimbula.py

@@ -16 +16 @@
 #------------------------------------------------------------------------------
 
-import unittest
 import os
 import sys
 import time
 import pypar
-
 import numpy as num
 
-
-
-
+#------------------------
+# ANUGA Modules
+#------------------------
 from anuga.utilities.numerical_tools import ensure_numeric
 from anuga.utilities.util_ext        import double_precision
@@ -48 +46 @@
 #--------------------------------------------------------------------------
 
-#mesh_filename = "merimbula_10785_1.tsh"
-mesh_filename = "merimbula_43200.tsh"
-#mesh_filename = "test-100.tsh"
-yieldstep = 50
+#mesh_filename = "merimbula_10785_1.tsh" ; x0 = 756000.0 ; x1 = 756500.0
+#mesh_filename = "merimbula_43200.tsh"   ; x0 = 756000.0 ; x1 = 756500.0
+mesh_filename = "test-100.tsh" ; x0 = 0.25 ; x1 = 0.5
+yieldstep = 5
 finaltime = 200
-quantity = 'stage'
 verbose = True
 
@@ -60 +57 @@
 #--------------------------------------------------------------------------
 class Set_Stage:
-    """Set an initial condition with constant water height, for x<x0
+    """Set an initial condition with constant water height, for x0<x<x1
     """
 
@@ -76 +73 @@
 if myid == 0:
     domain = create_domain_from_file(mesh_filename)
-    domain.set_quantity('stage', Set_Stage(756000.0, 756500.0, 2.0))
+    domain.set_quantity('stage', Set_Stage(x0, x1, 2.0))
 else:
     domain = None
@@ -86 +83 @@
 if myid == 0 and verbose: print 'DISTRIBUTING DOMAIN'
 domain = distribute(domain)
+
+domain.smooth = False
 
 #------------------------------------------------------------------------------