Changeset 6441 for branches/numpy/anuga/coordinate_transforms/geo_reference.py

Timestamp:

Mar 4, 2009, 8:26:22 AM (15 years ago)

Author:

rwilson

Message:

After changes to get_absolute, ensure_numeric, etc.

File:

• branches/numpy/anuga/coordinate_transforms/geo_reference.py (modified) (5 diffs)

branches/numpy/anuga/coordinate_transforms/geo_reference.py

@@ -8 +8 @@
 
 import types, sys
+import copy
+
 from anuga.utilities.numerical_tools import ensure_numeric
 from anuga.utilities.anuga_exceptions import ANUGAError, TitleError, \
@@ -240 +242 @@
     # @return The changed points.
     # @note If 'points' is a list then a changed list is returned.
-    # @note The input points data is changed.
     def change_points_geo_ref(self, points, points_geo_ref=None):
         """Change the geo reference of a list or numeric array of points to
@@ -250 +251 @@
         is_list = isinstance(points, list)
 
-        points = ensure_numeric(points, num.float)
+        points = ensure_numeric(copy.copy(points), num.float)
 
         # sanity checks
@@ -308 +309 @@
         is_list = isinstance(points, list)
 
-        # convert to numeric array
-        points = ensure_numeric(points, num.float)
+        # convert to numeric array, force a copy
+        points = ensure_numeric(copy.copy(points), num.float)
 
         # sanity checks
@@ -328 +329 @@
             points[:,0] += self.xllcorner
             points[:,1] += self.yllcorner
-            #self.is_absolute = True
 
         # if we got a list, return a list