Changeset 4884

Timestamp:

Dec 12, 2007, 3:25:15 PM (17 years ago)

Author:

ole

Message:

Shark bay embayment study.
Reduced verbose output further (ticket:238)

Files:

• anuga_core/source/anuga/fit_interpolate/fit.py (modified) (3 diffs)
• anuga_core/source/anuga_parallel/parallel_shallow_water.py (modified) (2 diffs)
• anuga_work/production/shark_bay_2007/project.py (modified) (1 diff)
• anuga_work/production/shark_bay_2007/run_shark_bay_embayment_study.py (added)
• anuga_work/production/shark_bay_2007/run_shark_bay_frequency_sweep.py (modified) (1 diff)

anuga_core/source/anuga/fit_interpolate/fit.py

@@ -268 +268 @@
                      in_and_outside_polygon(point_coordinates,
                                             self.mesh_boundary_polygon,
-                                            closed = True, verbose = verbose)
+                                            closed = True,
+                                            verbose = False) # There's too much output if True
         #print "self.inside_poly_indices",self.inside_poly_indices
         #print "self.outside_poly_indices",self.outside_poly_indices
@@ -274 +275 @@
         
         n = len(inside_poly_indices)
-        if verbose: print 'Building fitting matrix from %d points' %n        
+        #if verbose: print 'Building fitting matrix from %d points' %n        
+
         #Compute matrix elements for points inside the mesh
         triangles = self.mesh.triangles #Did this for speed, did ~nothing
@@ -347 +349 @@
                     # I suggest deleting this verbose output and make
                     # Geospatial_data more informative for txt files.
+                    #
+                    # I still think so (12/12/7, Ole).
             
 

anuga_core/source/anuga_parallel/parallel_shallow_water.py

@@ -152 +152 @@
             t0 = time.time()
             pypar.reduce(self.local_timestep, pypar.MIN, 0,
-                         buffer=self.global_timestep,
-                         bypass=True)
+                         buffer=self.global_timestep)#,
+                         #bypass=True)
 
         else:
@@ -178 +178 @@
         #Broadcast minimal timestep to all
         t0 = time.time()
-        pypar.broadcast(self.global_timestep, 0,
-                        bypass=True)
+        pypar.broadcast(self.global_timestep, 0)#,
+                        #bypass=True)
 
         self.communication_broadcast_time += time.time()-t0

anuga_work/production/shark_bay_2007/project.py

@@ -51 +51 @@
 source='shark_bay'
 
-boundary_event = 'july2006'
+#boundary_event = 'july2006'
 #boundary_event = '10000'
-#boundary_event = 'experimental'
+boundary_event = 'experimental'
 
 if setup =='trial':

anuga_work/production/shark_bay_2007/run_shark_bay_frequency_sweep.py

@@ -163 +163 @@
     print 'domain id', id(domain)
     domain.beta_h = 0
-    #domain.tight_slope_limiters = 1
+    domain.tight_slope_limiters = 1