Changeset 3193 for inundation/pmesh/mesh_interface.py

Timestamp:

Jun 21, 2006, 3:33:49 PM (18 years ago)

Author:

duncan

Message:

bug fixing

File:

• inundation/pmesh/mesh_interface.py (modified) (6 diffs)

inundation/pmesh/mesh_interface.py

@@ -22 +22 @@
 def create_mesh_from_regions(bounding_polygon,
                              boundary_tags,
-                             maximum_triangle_area,
+                             maximum_triangle_area=None,
                              filename=None,
                              interior_regions=None,
                              poly_geo_reference=None,
                              mesh_geo_reference=None,
-                             minimum_triangle_angle=28.0):
+                             minimum_triangle_angle=28.0,
+                             verbose=True):
     """Create mesh from bounding polygons, and resolutions.
 
@@ -36 +37 @@
     is a list of indices referring to segments associated with that tag
 
-    Resolution is the maximal area per triangle for the bounding polygon
-    (excluding interior regions, see later)
+    maximum_triangle_area is the maximal area per triangle
+    for the bounding polygon, excluding the  interior regions.
 
     Interior_regions is a list of tuples consisting of (polygon, resolution)
@@ -45 +46 @@
     pmesh functionality for that
 
-    poly_geo_reference is the geo_reference of the polygons.
-    The bounding polygon and the interior polygons
+    poly_geo_reference is the geo_reference of the bounding polygon and
+    the interior polygons.
     If none, assume absolute.  Please pass one though, since absolute
     references have a zone.
@@ -62 +63 @@
     """
     #FIXME (OLE-DSG)
-    # To do make maximum_triangle_area optional?
     # check the segment indexes - throw an error if they are out of bounds
     #(DSG) Yes!
@@ -117 +117 @@
         excluded_polygons = None
 
-    from Numeric import array
-
-
     # Convert bounding poly to absolute values
     # this sort of thing can be fixed with the geo_points class
@@ -157 +154 @@
     else:
         m.generate_mesh(minimum_triangle_angle=minimum_triangle_angle,
-                        maximum_triangle_area=maximum_triangle_area)        
+                             verbose=verbose)
         m.export_mesh_file(filename)