Changeset 3700 for anuga_core/source/anuga

Timestamp:

Oct 5, 2006, 3:42:07 PM (18 years ago)

Author:

ole

Message:

Improved runup script for Rosh

Location:

anuga_core/source/anuga

Files:

• abstract_2d_finite_volumes/util.py (modified) (1 diff)
• examples/runup_convergence.py (modified) (4 diffs)

anuga_core/source/anuga/abstract_2d_finite_volumes/util.py

@@ -810 +810 @@
 
     import pylab as p1
-    import mpl3d.mplot3d as p3
+    if surface is True:
+        import mpl3d.mplot3d as p3
 
     if report == True:    

anuga_core/source/anuga/examples/runup_convergence.py

@@ -138 +138 @@
 w0 = domain.get_maximum_inundation_elevation()
 x0, y0 = domain.get_maximum_inundation_location()
-print 'Coastline elevation = %.2f at (%.2f, %.2f)' %(w0, x0, y0)
+print
+print 'Coastline elevation = %.2f at (x,y)=(%.2f, %.2f)' %(w0, x0, y0)
 
 # Sanity check
 w_i = domain.get_quantity('stage').get_values(interpolation_points=[[x0,y0]])
-print 'Interpolated elevation at (%.2f, %.2f) is %.2f' %(x0, y0, w_i) 
+print 'Interpolated elevation at (x,y)=(%.2f, %.2f) is %.2f' %(x0, y0, w_i) 
 
 
@@ -155 +156 @@
     w = domain.get_maximum_inundation_elevation()
     x, y = domain.get_maximum_inundation_location()
-    print '  Coastline elevation = %.2f at (x,y)=(%.2f, %.2f)' %(w, x, y)    
+    print '  Coastline elevation = %.2f at (x,y)=(%.2f, %.2f)' %(w, x, y)
+    print    
 
     if w > w_max:
@@ -163 +165 @@
 
 
+print '**********************************************'
 print 'Max coastline elevation = %.2f at (%.2f, %.2f)' %(w_max, x_max, y_max)
 print 'Run up distance = %.2f' %sqrt( (x_max-x0)**2 + (y_max-y0)**2 )
-
+print '**********************************************'
 
 
@@ -172 +175 @@
 #---------------------------------------------------------------
 
-# Generate time series of "gauge" situated at right hand boundary
+# Generate time series of one "gauge" situated at right hand boundary
 from anuga.abstract_2d_finite_volumes.util import sww2timeseries
 production_dirs = {'.': 'test'}