Changeset 483

Timestamp:

Nov 4, 2004, 10:18:24 AM (20 years ago)

Author:

ole

Message:

 

Location:

inundation/ga/storm_surge/pyvolution

Files:

• flatbed.py (modified) (2 diffs)
• flatbed_compare.py (modified) (3 diffs)

inundation/ga/storm_surge/pyvolution/flatbed.py

@@ -14 +14 @@
 
 #Create basic mesh
-N = 8
-points, vertices, boundary = rectangular(N, N)
+N = 150
+points, vertices, boundary = rectangular(N, N, 1000, 1000)
 
 #Create shallow water domain
 domain = Domain(points, vertices, boundary)
-domain.smooth = False
-domain.visualise = True
-domain.default_order=2
+domain.smooth = True
+domain.visualise = False
+domain.store = True
+domain.default_order=1
+
+#Set driving forces
+manning = 0.07
+manning = 0.0
+inflow_level = 10.0
+domain.set_quantity('friction', manning)
 
 
@@ -27 +34 @@
 # Boundary conditions
 Br = Reflective_boundary(domain)
-Bd = Dirichlet_boundary([0.2,0.,0.])
+Bd = Dirichlet_boundary([inflow_level,0.,0.])
 
 domain.set_boundary({'left': Bd, 'right': Br, 'top': Br, 'bottom': Br})
+domain.check_integrity()
 
-
-domain.check_integrity()
 
 ######################
 #Evolution
-#for t in domain.evolve(yieldstep = 0.05, finaltime = 1.):
-#    domain.write_time()
-
-
-#import sys; sys.exit()
-
-######################
-#Evolution
-for t in domain.evolve(yieldstep = 0.01, finaltime = 0.03):
+for t in domain.evolve(yieldstep = 10, finaltime = 500):
     domain.write_time()
 
-print domain.quantities['level'].centroid_values[:4]
-print domain.quantities['xmomentum'].centroid_values[:4]
-print domain.quantities['ymomentum'].centroid_values[:4]
-domain.distribute_to_vertices_and_edges()
-print
-print domain.quantities['level'].vertex_values[:4,0]
-print domain.quantities['xmomentum'].vertex_values[:4,0]
-print domain.quantities['ymomentum'].vertex_values[:4,0]  
+print 'Done'    
+

inundation/ga/storm_surge/pyvolution/flatbed_compare.py

@@ -1 +1 @@
 """Example of shallow water wave equation.
 
-Comparison of flatbed between mark 2(Grohm) and mark 3
+Generate slope
 
 """
@@ -10 +10 @@
 from mesh_factory import rectangular
 from shallow_water import Domain, Reflective_boundary, Dirichlet_boundary,\
-     Transmissive_boundary, Constant_height
+     Constant_height
 from Numeric import array
 
 #Create basic mesh
-N = 10
+N = 8
 points, vertices, boundary = rectangular(N, N)
 
 #Create shallow water domain
 domain = Domain(points, vertices, boundary)
-
-if N < 20:
-    domain.visualise = True
-else:
-    domain.visualise = False
-    domain.store = True
-    domain.format = 'sww'
-    domain.filename = 'compare_py3'
-
-#domain.visualise = False
 domain.smooth = False
-domain.default_order = 1
-
-
-print 'Field values'
-domain.set_quantity('elevation', 0.0)
-domain.set_quantity('friction', 1)
+domain.visualise = True
+domain.default_order=2
 
 
@@ -42 +28 @@
 Br = Reflective_boundary(domain)
 Bd = Dirichlet_boundary([0.2,0.,0.])
-Bt = Transmissive_boundary(domain)
 
-domain.set_boundary({'left': Bd, 'right': Bt, 'top': Bt, 'bottom': Bt})
+domain.set_boundary({'left': Bd, 'right': Br, 'top': Br, 'bottom': Br})
+
+
 domain.check_integrity()
-
-print domain.quantities['elevation'].centroid_values[:4]
-print domain.quantities['friction'].centroid_values[:4]
 
 ######################
 #Evolution
-for t in domain.evolve(yieldstep = 0.01, finaltime = 0.5):
-    domain.write_time()
-    #print
-
-#print domain.quantities['level'].centroid_values
-#print domain.quantities['xmomentum'].centroid_values
-#print domain.quantities['ymomentum'].centroid_values    
-
-#print 'R'
-#print domain.quantities['level'].edge_values
+#for t in domain.evolve(yieldstep = 0.05, finaltime = 1.):
+#    domain.write_time()
 
 
+#import sys; sys.exit()
+
+######################
+#Evolution
+for t in domain.evolve(yieldstep = 0.01, finaltime = 0.03):
+    domain.write_time()
+
+print domain.quantities['level'].centroid_values[:4]
+print domain.quantities['xmomentum'].centroid_values[:4]
+print domain.quantities['ymomentum'].centroid_values[:4]
+domain.distribute_to_vertices_and_edges()
+print
+print domain.quantities['level'].vertex_values[:4,0]
+print domain.quantities['xmomentum'].vertex_values[:4,0]
+print domain.quantities['ymomentum'].vertex_values[:4,0]