Changeset 2505

Timestamp:

Mar 9, 2006, 11:27:05 AM (18 years ago)

Author:

sexton

Message:

adding some diagnostics

File:

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

inundation/pyvolution/data_manager.py

@@ -1197 +1197 @@
 
     dem_elevation_r = reshape(dem_elevation, (nrows, ncols))
+    totalnopoints = nrows*ncols
 
     # calculating number of NODATA_values for each row
     nn = 0
     no_counts = []
+    print 'I think NODATA_values is:    ', NODATA_value
     for i in range(nrows):
         v = [dem_elevation_r[i,:]]
@@ -1207 +1209 @@
         no_counts.append(this_count)
 
-    if verbose and nn > 0: print 'There are %d NODATA_values in the elevation' %nn
+    if verbose and nn > 0:
+        print 'There are %d values in the elevation', %totalnopoints
+        print 'There are %d NODATA_values in the elevation' %nn
     
     # dimension definitions
@@ -1214 +1218 @@
     
     nopoints = nrows_in_bounding_box*ncols_in_bounding_box-nn
+    #nopoints = nrows_in_bounding_box*ncols_in_bounding_box
     outfile.createDimension('number_of_points', nopoints)
     outfile.createDimension('number_of_dimensions', 2) #This is 2d data
@@ -1239 +1244 @@
         else:
             newcols = ncols_in_bounding_box
-
+        
         telev = zeros(newcols, Float)
         tpoints = zeros((newcols, 2), Float)