Changeset 5736 for anuga_core/source/anuga/coordinate_transforms

Timestamp:

Sep 5, 2008, 11:30:19 AM (17 years ago)

Author:

ole

Message:

Implemented run-time version of get_energy_through_cross_section and tested it.

File:

• anuga_core/source/anuga/coordinate_transforms/geo_reference.py (modified) (1 diff)

anuga_core/source/anuga/coordinate_transforms/geo_reference.py

@@ -212 +212 @@
         assert points.shape[1] == 2, msg                
 
-            
+        
+        # FIXME (Ole): Could also check if zone, xllcorner, yllcorner 
+        # are identical in the two geo refs.    
         if points_geo_ref is not self:
-            #add point geo ref to points
+            # If georeferences are different
+        
             if not points_geo_ref is None:
+                # Convert points to absolute coordinates
                 points[:,0] += points_geo_ref.xllcorner 
                 points[:,1] += points_geo_ref.yllcorner 
         
-            #subtract primary geo ref from points
+            # Make points relative to primary geo reference
             points[:,0] -= self.xllcorner 
             points[:,1] -= self.yllcorner