Changeset 1295 for inundation/ga/storm_surge/pyvolution/netherlands.py

Timestamp:

May 8, 2005, 11:50:32 PM (20 years ago)

Author:

steve

Message:

Playing with coloured surfaces

File:

• inundation/ga/storm_surge/pyvolution/netherlands.py (modified) (6 diffs)

inundation/ga/storm_surge/pyvolution/netherlands.py

@@ -63 +63 @@
                 z[i] = -x[i]/20+0.4
 
+            if (x[i] - 0.72)**2 + (y[i] - 0.4)**2 < 0.05**2:# or\
+                   #x[i] > 0.75 and x[i] < 0.9 and y[i] > 0.35 and y[i] < 0.45:
+                z[i] = -x[i]/20+0.4
+
             #Wall
             if x[i] > 0.995:
@@ -79 +83 @@
 N = 130 #size = 33800
 N = 600 #Size = 720000
-N = 30
+N = 40
 
 #N = 15
@@ -100 +104 @@
 print "Number of triangles = ", len(domain)
 
+#Set bed-slope and friction
+inflow_stage = 0.4
+manning = 0.02
+Z = Weir(inflow_stage)
 
 if N > 150:
@@ -112 +120 @@
 else:
     domain.visualise = True
+    domain.visualise_color_stage = True
+    domain.visualise_timer = True
     domain.checkpoint = False
     domain.store = False
 
-
-#Set bed-slope and friction
-inflow_stage = 0.1
-manning = 0.02
-Z = Weir(inflow_stage)
 
 print 'Field values'
@@ -134 +139 @@
 
 #Constant inflow
-Bd = Dirichlet_boundary([2*inflow_stage, 0.0, 0.0])
+Bd = Dirichlet_boundary([inflow_stage, 0.0, 0.0])
 
 
@@ -151 +156 @@
 t0 = time.time()
 
-for t in domain.evolve(yieldstep = 0.2, finaltime = 10.0):
+for t in domain.evolve(yieldstep = 0.05, finaltime = 10.0):
     domain.write_time()