Changeset 2535

Timestamp:

Mar 13, 2006, 2:46:33 PM (18 years ago)

Author:

ole

Message:

Histogram of mesh areas

Files:

• development/okushiri_2005/lwru2.py (modified) (2 diffs)
• development/stochastic_study/run_model.py (modified) (1 diff)
• inundation/pyvolution/mesh.py (modified) (2 diffs)

development/okushiri_2005/lwru2.py

@@ -115 +115 @@
 
 
+
 import sys, os
 base = os.path.basename(sys.argv[0])
@@ -124 +125 @@
 print 'Number of triangles = ', len(domain)
 print 'The extent is ', domain.get_extent()
+print domain.statistics()
 
 

development/stochastic_study/run_model.py

@@ -50 +50 @@
 print 'Number of triangles = ', len(domain)
 print 'The extent is ', domain.get_extent()
+print domain.statistics()
+
+import sys; sys.exit()
+
 
 domain.set_datadir('.')

inundation/pyvolution/mesh.py

@@ -631 +631 @@
         #Setup 10 bins for area histogram
         m = max(areas)
-        bins = arange(0, m/10, m)
+        bins = arange(0., m, m/10)
         hist = histogram(areas, bins)
 
@@ -643 +643 @@
         str += '    A in [%f, %f]\n' %(min(areas), max(areas))
         str += '    Histogram:\n'
+
+        hi = bins[0]
         for i, count in enumerate(hist):
-            str += '      >=%f: %d\n' %(bins[i], count)
+            lo = hi
+            if i+1 < len(bins):
+                #Open upper interval                
+                hi = bins[i+1]
+                str += '      [%f, %f[: %d\n' %(lo, hi, count)                
+            else:
+                #Closed upper interval
+                hi = m
+                str += '      [%f, %f]: %d\n' %(lo, hi, count)                
+                
                       
         str += 'Boundary:\n'