Changeset 319

Timestamp:

Sep 17, 2004, 5:03:05 PM (20 years ago)

Author:

ole

Message:

 

File:

• inundation/ga/storm_surge/pyvolution/cornell_room.py (modified) (2 diffs)

inundation/ga/storm_surge/pyvolution/cornell_room.py

@@ -18 +18 @@
 print 'Creating domain'
 points, triangles, values = from_polyfile('cornell_room_medres')
+
 
 #Create shallow water domain
@@ -62 +63 @@
 #
 print 'Initial condition'
+
+#Define water height as a lump in one corner
+def height(x, y):
+    from Numeric import zeros, Float
+    
+    N = len(x)
+    assert N == len(y)    
+
+    xmin = min(x); xmax = max(x)
+    ymin = min(y); ymax = max(y)
+
+    xrange = xmax - xmin
+    yrange = ymax - ymin    
+
+    z = zeros(N, Float)
+    for i in range(N):
+        if x[i] <= xmin + 0.25*xrange and y[i] <= ymin + 0.25*yrange:
+            z[i] = 300
+
+    return z
+
+domain.set_quantity('level', height)
+
 E = domain.quantities['elevation'].vertex_values
-print max(E)
-L = E.copy()
-N = len(L)
-L[:N/4] += 100
-domain.set_quantity('level', L)
+L = domain.quantities['level'].vertex_values
+domain.set_quantity('level', E+L)
 
 #Evolve
-for t in domain.evolve(yieldstep = 0.01, finaltime = 1.0):
+for t in domain.evolve(yieldstep = 0.05, finaltime = 5.0):
     domain.write_time()