Changeset 2349

Timestamp:

Feb 7, 2006, 4:49:08 PM (19 years ago)

Author:

ole

Message:

 

File:

• inundation/pmesh/create_mesh.py (modified) (2 diffs)

inundation/pmesh/create_mesh.py

@@ -11 +11 @@
 
 
-def convert_points_from_latlon_to_utm(polygon, refzone):
+def convert_points_from_latlon_to_utm(polygon, refzone=None):
     points = []
     for point in polygon:
         zone, easting, northing  = redfearn(point[0], point[1])
         #FIXME: Use point.latitude etc once we have a proper point set
-        assert zone == refzone
+
+        if refzone is None:
+            zone = refzone
+        else:
+            assert zone == refzone
+        
         points.append([easting, northing])
 
-    return points    
+    return points, zone    
 
 
@@ -124 +129 @@
 
     #Convert to UTM
-    bounding_polygon = convert_points_from_latlon_to_utm(bounding_polygon, refzone)    
+    bounding_polygon, refzone = convert_points_from_latlon_to_utm(bounding_polygon)