Changeset 3705

Timestamp:

Oct 9, 2006, 11:08:10 AM (16 years ago)

Author:

ole

Message:

Implemented runup height detection in okushiri validation

File:

• anuga_validation/okushiri_2005/lwru2.py (modified) (4 diffs)

anuga_validation/okushiri_2005/lwru2.py

@@ -33 +33 @@
     domain = Domain(project.mesh_filename, use_cache=True, verbose=True)
 else:
-    print 'Creating regular mesh'
+    print 'Creating regular from regular mesh'
     N = 150
     points, vertices, boundary = rectangular_cross(N, N/5*3,
                                                    len1=5.448, len2=3.402)
-    print 'Creating domain'
-
-    #domain = Domain(points, vertices, boundary)
-    domain = cache(Domain, (points, vertices, boundary))
-
+    domain = Domain(points, vertices, boundary)
 
 domain.set_name('lwru2')
 domain.set_default_order(2)
 domain.set_minimum_storable_height(0.001)
-
-#domain.check_integrity()
-print 'Number of triangles = ', len(domain)
-print 'The extent is ', domain.get_extent()
+domain.check_integrity()
 print domain.statistics()
 
@@ -59 +52 @@
 domain.set_quantity('stage', 0.0)
 domain.set_quantity('elevation',
-                    filename = project.bathymetry_filename[:-4] + '.pts',
-                    alpha = 0.02,                    
-                    verbose = True,
-                    use_cache = True)
+                    filename=project.bathymetry_filename[:-4] + '.pts',
+                    alpha=0.02,                    
+                    verbose=True,
+                    use_cache=True)
 
 #-------------------------
@@ -71 +64 @@
 function = file_function(project.boundary_filename[:-4] + '.tms',
                          domain,
-                         verbose = True)
+                         verbose=True)
 
 Bts = Transmissive_Momentum_Set_Stage_boundary(domain, function)
 
-#Set boundary conditions
 if use_variable_mesh is True:
     domain.set_boundary({'wave': Bts, 'wall': Br})
@@ -88 +80 @@
 t0 = time.time()
 
+w_max = 0
 for t in domain.evolve(yieldstep = 0.05, finaltime = 22.5):
     domain.write_time()
 
+    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
+    
+    if w > w_max:
+        w_max = w
+        x_max = x
+        y_max = y
+
+
+print '**********************************************'
+print 'Max coastline elevation = %.2f at (%.2f, %.2f)' %(w_max, x_max, y_max)
+print '**********************************************'
+    
+   
+
 print 'That took %.2f seconds' %(time.time()-t0)