Changeset 2207

Timestamp:

Jan 13, 2006, 4:53:30 PM (19 years ago)

Author:

steve

Message:

 

Location:

inundation/pyvolution

Files:

• netherlands.py (modified) (5 diffs)
• shallow_water_vtk.py (added)

inundation/pyvolution/netherlands.py

@@ -15 +15 @@
      Transmissive_boundary, Constant_height, Constant_stage
 
-from mesh_factory import rectangular
+from mesh_factory import rectangular_cross
 from Numeric import array
 
@@ -83 +83 @@
 N = 130 #size = 33800
 N = 600 #Size = 720000
-N = 50
+N = 20
 
 
@@ -90 +90 @@
 print 'Creating domain'
 #Create basic mesh
-points, elements, boundary = rectangular(N, N)
+points, elements, boundary = rectangular_cross(N, N)
 
 #Create shallow water domain
@@ -168 +168 @@
 #V.update_quantity('elevation')
 
-for t in domain.evolve(yieldstep = 0.005, finaltime = 15.0):
+
+
+
+for t in domain.evolve(yieldstep = 0.05, finaltime = 1.0):
     domain.write_time()
     #domain.write_boundary()
@@ -174 +177 @@
     print domain.quantities['stage'].get_values(location='centroids',
                                                 indices=[0])
+    #domain.visualiser.update_quantity('elevation')
+    #time.sleep(0.1)
+    #raw_input('pause>')
     #V.update_quantity('stage')
     #rpdb.set_active()