Changeset 1585 for inundation/ga/storm_surge/parallel/run_parallel_sw_merimbula.py

Timestamp:

Jul 6, 2005, 9:55:40 PM (19 years ago)

Author:

steve

Message:

Better partitioning of meribula data set

File:

• inundation/ga/storm_surge/parallel/run_parallel_sw_merimbula.py (modified) (6 diffs)

inundation/ga/storm_surge/parallel/run_parallel_sw_merimbula.py

@@ -59 +59 @@
 # mesh partition routines
 
-from pmesh_divide import pmesh_divide
+from pmesh_divide import pmesh_divide, pmesh_divide_steve
 from build_submesh import *
 from build_local import *
@@ -105 +105 @@
 
     nodes, triangles, boundary, triangles_per_proc, quantities = \
-         pmesh_divide(domain_full, nx, ny)
+         pmesh_divide_steve(domain_full, nx, ny)
 
     rect = array(domain_full.xy_extent, Float)
@@ -146 +146 @@
 
 domain.initialise_visualiser(rect=rect,scale_z =0.1)
+#domain.initialise_visualiser(rect=rect)
 domain.visualise_color_stage = True
 domain.default_order = 1
@@ -156 +157 @@
 domain.set_boundary( {'outflow': R, 'inflow': R, 'inner':R, 'exterior': R, 'open':R} )
 
-print quantities.keys()
+
 domain.set_quantity('stage', quantities['stage'])
 domain.set_quantity('elevation', quantities['elevation'])
@@ -163 +164 @@
 # Evolution
 t0 = time.time()
-domain.visualise = True
+
+print 'Processor %d No of elements %d'%(domain.processor,domain.number_of_elements)
 yieldstep = 10.0
-finaltime = 1000.0
+finaltime = 500.0
 
 #yieldstep = 10
@@ -173 +175 @@
     if myid == 0:
         domain.write_time()
-        print domain.quantities['stage'].centroid_values[0]
+        print domain.quantities['stage'].get_integral()
 
 if myid == 0: