Changeset 4350

Timestamp:

Apr 3, 2007, 3:51:30 PM (18 years ago)

Author:

sexton

Message:

update of Rosh' example

File:

• anuga_validation/solitary_waves/solitary_wave_runup.py (modified) (5 diffs)

anuga_validation/solitary_waves/solitary_wave_runup.py

@@ -22 +22 @@
 from anuga.shallow_water import Transmissive_Momentum_Set_Stage_boundary
 from abstract_2d_finite_volumes.util import file_function
-from pylab import plot, xlabel, ylabel, title, ion, close, savefig,\
-     figure, axis, legend, grid, hold
+#from pylab import plot, xlabel, ylabel, title, ion, close, savefig,\
+#     figure, axis, legend, grid, hold
 
 
@@ -39 +39 @@
 amplitude = 0.5     # Solitary wave height H
 normalized_amplitude = amplitude/offshore_depth 
-simulation_name = 'runup_convergence'   
+ 
 coastline_x = -highest_point/slope
 
@@ -54 +54 @@
 
 # Structured mesh
-dx = 5           # Resolution: Length of subdivisions on x axis (length)
-dy = 5           # Resolution: Length of subdivisions on y axis (width)
+dx = 40.           # Resolution: Length of subdivisions on x axis (length)
+dy = 40.          # Resolution: Length of subdivisions on y axis (width)
 
 length = east-west
 width = north-south
-points, vertices, boundary = rectangular_cross(length/dx, width/dy,
-                                               len1=length, len2=width,
-                                               origin = (west, south)) 
-
-domain = Domain(points, vertices, boundary) # Create domain
+##points, vertices, boundary = rectangular_cross(length/dx, width/dy,
+##                                               len1=length, len2=width,
+##                                               origin = (west, south)) 
+##
+##domain = Domain(points, vertices, boundary) # Create domain
 
 
@@ -70 +70 @@
 interior_polygon = [[400,north-10],[west+10,north-10],
                     [west+10,south+10],[400,south+10]]
+res = 0.065
+simulation_name = 'runup_convergence' + str(res)
 meshname = simulation_name + '.msh'
 create_mesh_from_regions(polygon,
                          boundary_tags={'top': [0], 'left': [1],
                                         'bottom': [2], 'right': [3]},
-                         maximum_triangle_area=dx*dy/4,
+                         maximum_triangle_area=10.,#dx*dy/4.,
                          filename=meshname,
-                         interior_regions=[[interior_polygon,dx*dy/32]])
+                         #interior_regions=[[interior_polygon,dx*dy/32.]])
+                         interior_regions=[[interior_polygon,res]])
 
 
@@ -174 +177 @@
 print '**********************************************'
 
-
+import sys; sys.exit()
 #-----------------------------------------------------------------------------
 # Interrogate further