Changeset 458 for inundation/ga/storm_surge/pyvolution/flatbed_compare.py

Timestamp:

Oct 28, 2004, 2:40:44 PM (20 years ago)

Author:

ole

Message:

Fixed friction bug

File:

• inundation/ga/storm_surge/pyvolution/flatbed_compare.py (modified) (3 diffs)

inundation/ga/storm_surge/pyvolution/flatbed_compare.py

@@ -14 +14 @@
 
 #Create basic mesh
-N = 50
+N = 10
 points, vertices, boundary = rectangular(N, N)
 
@@ -28 +28 @@
     domain.filename = 'compare_py3'
 
+#domain.visualise = False
 domain.smooth = False
-domain.default_order = 2
+domain.default_order = 1
 
 
 print 'Field values'
 domain.set_quantity('elevation', 0.0)
-domain.set_quantity('friction', 0.0)
+domain.set_quantity('friction', 1)
 
 
@@ -46 +47 @@
 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.2):
+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
+
+