Changeset 1141 for inundation/ga/storm_surge/pyvolution/data_manager.py

Timestamp:

Mar 24, 2005, 1:49:49 PM (20 years ago)

Author:

ole

Message:

Verbosity in sww2asc

File:

• inundation/ga/storm_surge/pyvolution/data_manager.py (modified) (4 diffs)

inundation/ga/storm_surge/pyvolution/data_manager.py

@@ -1121 +1121 @@
 
     #Get quantity and reduce if applicable
+    if verbose: print 'Reading quantity %s' %quantity    
     q = fid.variables[quantity][:]
 
     
     if len(q.shape) == 2:
+        if verbose: print 'Reducing quantity %s' %quantity            
         q_reduced = zeros( number_of_points, Float ) 
 
@@ -1133 +1135 @@
 
     #Now q has dimension: number_of_points
+
     
     #Write prj file
+    if verbose: print 'Writing %s' %prjfile
     prjid = open(prjfile, 'w')
     prjid.write('Projection    %s\n' %'UTM')
@@ -1148 +1152 @@
 
     #Create grid
+    if verbose: print 'Creating grid'
     ncols = int((xmax-xmin)/cellsize)+1
     nrows = int((ymax-ymin)/cellsize)+1
@@ -1166 +1171 @@
 
     #Interpolate
+    
     from least_squares import Interpolation
     from util import inside_polygon
     
     interp = Interpolation(vertex_points, volumes, grid_points, alpha=0.0,
-                           precrop = False)
+                           precrop = False, verbose = verbose)
 
     #Interpolate using quantity values
+    if verbose: print 'Interpolating'
     grid_values = interp.interpolate(q).flat
 
     #Write
+    if verbose: print 'Writing %s' %ascfile    
     ascid = open(ascfile, 'w')