Changeset 6192 for anuga_core/source/anuga/abstract_2d_finite_volumes/util.py

Timestamp:

Jan 19, 2009, 6:49:40 AM (16 years ago)

Author:

ole

Message:

Work on Patong and some diagnostic output

File:

• anuga_core/source/anuga/abstract_2d_finite_volumes/util.py (modified) (4 diffs)

anuga_core/source/anuga/abstract_2d_finite_volumes/util.py

@@ -392 +392 @@
 
         if boundary_polygon is not None:
-            #removes sts points that do not lie on boundary
+            # Remove sts points that do not lie on boundary
             # FIXME(Ole): Why don't we just remove such points from the list of points and associated data?
             # I am actually convinced we can get rid of neighbour_gauge_id altogether as the sts file is produced using the ordering file.
@@ -404 +404 @@
             for i in range(len(boundary_polygon)):
                 for j in range(len(x)):
-                    if num.allclose(vertex_coordinates[j],boundary_polygon[i],1e-4):
+                    if num.allclose(vertex_coordinates[j],
+                                    boundary_polygon[i], 1e-4):
                         #FIXME:
                         #currently gauges lat and long is stored as float and
@@ -475 +476 @@
         if boundary_polygon is not None:
             #removes sts points that do not lie on boundary
-            quantities[name] = num.take(quantities[name],gauge_id,1)
+            quantities[name] = num.take(quantities[name], gauge_id, 1)
             
     # Close sww, tms or sts netcdf file         
@@ -488 +489 @@
         #vertex coordinates is position of urs gauges
 
+    if verbose:
+        print 'Call interpolation function'
     # Return Interpolation_function instance as well as
     # starttime for use to possible modify that of domain