Changeset 2625 for inundation/parallel/run_parallel_merimbula.py

Timestamp:

Mar 29, 2006, 9:41:58 AM (19 years ago)

Author:

linda

Message:

Added ghost boundary layers

File:

• inundation/parallel/run_parallel_merimbula.py (modified) (5 diffs)

inundation/parallel/run_parallel_merimbula.py

@@ -87 +87 @@
     # read in the test files
 
-    filename = 'test-100.tsh'
-#    filename = 'merimbula_10785.tsh'
+#    filename = 'test-100.tsh'
+    filename = 'merimbula_10785.tsh'
     nx = numprocs
     ny = 1
@@ -95 +95 @@
 
     domain_full = pmesh_to_domain_instance(filename, Advection_Domain)
-    domain_full.set_quantity('stage', Set_Stage(200.0,300.0,1.0))
-#    domain.set_quantity('stage', Set_Stage(756000.0,756500.0,4.0))
+#    domain_full.set_quantity('stage', Set_Stage(200.0,300.0,1.0))
+    domain_full.set_quantity('stage', Set_Stage(756000.0,756500.0,4.0))
 
     nodes, triangles, boundary, triangles_per_proc, quantities  =\
             pmesh_divide(domain_full, nx, ny)
-
+    
     # subdivide the mesh
     rect = array(domain_full.xy_extent, Float)
-
-    print rect
 
     rect = array(rect, Float)
@@ -128 +126 @@
              rec_submesh(0)
 
-#if myid == 0:
-#    print 'ghost'
-#    print ghost_recv_dict
-#
-#if myid == 0:
-#    print 'full'
-#    print full_send_dict
-
-
 pypar.broadcast(rect,0)
-print rect
 
 domain = Parallel_Domain(points, vertices, boundary,
@@ -150 +138 @@
 
 T = Transmissive_boundary(domain)
-R = Reflective_boundary(domain)
-domain.set_boundary( {'outflow': R, 'inflow': R, 'inner':R, 'exterior': R, 'open':R} )
+#R = Reflective_boundary(domain)
+domain.set_boundary( {'outflow': T, 'inflow': T, 'inner':T, 'exterior': T, 'open':T, 'ghost':None} )
 
 
@@ -170 +158 @@
 
 
+from norms import linf_norm
 
 #yieldstep = 1
+#finaltime = 4000
+
 yieldstep = 1000
 finaltime = 50000
-# finaltime = 4000
+
 for t in domain.evolve(yieldstep = yieldstep, finaltime = finaltime):
     if myid == 0:
         domain.write_time()
+##    Stage = domain.quantities['stage']
+##    max_stage = max(max(Stage.edge_values))
+##    print max_stage
+##    if max_stage > 4.0:
+##        print "######################", max_stage
 
 if myid == 0: