Changeset 1080

Timestamp:

Mar 15, 2005, 4:01:34 PM (20 years ago)

Author:

ole

Message:

More reporting in ferret2sww

File:

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

inundation/ga/storm_surge/pyvolution/data_manager.py

@@ -1442 +1442 @@
     outfile.variables['x'][:] = x - xllcorner
     outfile.variables['y'][:] = y - yllcorner
-    outfile.variables['z'][:] = 0.0
+    outfile.variables['z'][:] = z
     outfile.variables['elevation'][:] = z  #FIXME HACK
     outfile.variables['time'][:] = times
@@ -1466 +1466 @@
                 ymomentum[j,i] = vspeed[j,k,l]/100*h
                 i += 1
+
+
+    if verbose:
+        x = outfile.variables['x'][:]
+        y = outfile.variables['y'][:]        
+        
+        print '------------------------------------------------'
+        print 'Statistics of output file:'
+        print '  Name: %s' %swwname
+        print '  Reference:'
+        print '    Lower left corner: [%f, %f]'\
+              %(xllcorner, yllcorner)
+        print '  Extent:'
+        print '    x [m] in [%f, %f], len(x) == %d'\
+              %(min(x.flat), max(x.flat), len(x.flat))
+        print '    y [m] in [%f, %f], len(y) == %d'\
+              %(min(y.flat), max(y.flat), len(y.flat))
+        print '    t [s] in [%f, %f], len(t) == %d'\
+              %(min(times), max(times), len(times))
+        print '  Quantities [SI units]:'
+        for name in ['stage', 'xmomentum', 'ymomentum']:
+            q = outfile.variables[name][:].flat
+            print '    %s in [%f, %f]' %(name, min(q), max(q))
+            
+        
+
         
     outfile.close()