Changeset 3115 for production/wollongong_2006/run_flagstaff.py

Timestamp:

Jun 8, 2006, 2:10:15 PM (19 years ago)

Author:

ole

Message:

Got parallel wollongong going and flagged two issues with FIXME

File:

• production/wollongong_2006/run_flagstaff.py (modified) (7 diffs)

production/wollongong_2006/run_flagstaff.py

@@ -54 +54 @@
 print 'I am processor %d of %d on node %s' %(myid, numprocs, processor_name)
 
-bounding_rectangle = zeros(4, Float) # Buffer for results
-
 
 #------------------------------------------------------------------------------
@@ -94 +92 @@
                                     interior_regions=project.interior_regions)
     
-    ## Add buildings
-    #dict = importUngenerateFile(project.buildings_filename)
-    #Segment.set_default_tag('wall') # This should bind to a Reflective boundary
-    #mesh.addVertsSegs(dict)
-
-  
     # Add buildings
     # This should bind to a Reflective boundary
@@ -117 +109 @@
     print domain.statistics()
 
-    domain.set_name(project.basename)
+    #domain.set_name(project.basename)
     domain.set_datadir(project.outputdir)
-    domain.set_quantities_to_be_stored(['stage', 'xmomentum', 'ymomentum'])
-
-    bounding_rectangle[:] = domain.get_extent(absolute=True)
+    #domain.set_quantities_to_be_stored(['stage', 'xmomentum', 'ymomentum'])
+    domain.set_quantities_to_be_stored(None)
+
+    domain.set_quantity('elevation', 
+                        filename=project.demname + '.pts',
+                        use_cache=True,
+                        verbose=True)    
+
 
     # Subdivide the mesh
@@ -161 +158 @@
 
 
-pypar.broadcast(bounding_rectangle, 0)
-
 # Build the domain for this processor
 domain = Parallel_Domain(points, vertices, boundary,
@@ -168 +163 @@
                          ghost_recv_dict = ghost_recv_dict)
 
+# FIXME (Ole): Name currently has to be set here to get the processor number
+# right. It would be easy to build into Parallel_Domain
+domain.set_name(project.basename)
 
 #------------------------------------------------------------------------------
@@ -175 +173 @@
 domain.set_quantity('stage', project.initial_sealevel)
 domain.set_quantity('friction', 0.03)
-domain.set_quantity('elevation', 
-                    filename=project.demname + '.pts',
-                    use_cache=True,
-                    verbose=True)
+
+#
+# FIXME (Ole): This one segfaults which is bad, because set_quantity is
+# time consuming and should be done here rather than on processor 0
+#
+#domain.set_quantity('elevation', 
+#                    filename=project.demname + '.pts',
+#                    use_cache=True,
+#                    verbose=True)
 
 
@@ -196 +199 @@
                      'ghost': None})
 
-                            
+
+
+print 'P%d: Ready to evolve. Value of store is %s' %(myid, str(domain.store))
+
+
 #------------------------------------------------------------------------------
 # Evolve system through time