Changeset 1900

Timestamp:

Oct 12, 2005, 12:20:22 PM (19 years ago)

Author:

ole

Message:

Changed interpolation points to absolute UTM coordinates according to new
file_function interface

File:

• inundation/pyvolution/generic_boundary_conditions.py (modified) (2 diffs)

inundation/pyvolution/generic_boundary_conditions.py

@@ -197 +197 @@
         boundary_keys = domain.boundary.keys()
 
-        #Record ordering #FIXME: should this happen in domain.py
+
+        xllcorner = domain.geo_reference.get_xllcorner()
+        yllcorner = domain.geo_reference.get_yllcorner()        
+        
+
+        #Record ordering #FIXME: should this happen in domain.py?
         self.boundary_indices = {}
 
@@ -206 +211 @@
             x2 = V[vol_id, 4]; y2 = V[vol_id, 5]
 
-            #Midpoints
+            #Compute midpoints
             if edge_id == 0: m = array([(x1 + x2)/2, (y1 + y2)/2])
             if edge_id == 1: m = array([(x0 + x2)/2, (y0 + y2)/2])
             if edge_id == 2: m = array([(x1 + x0)/2, (y1 + y0)/2])
+
+            #Convert to absolute UTM coordinates
+            m[0] += xllcorner
+            m[1] += yllcorner
+            
+            #Register point and index
             self.midpoint_coordinates[i,:] = m
-
             self.boundary_indices[(vol_id, edge_id)] = i