Changeset 449

Timestamp:

Oct 26, 2004, 6:01:37 PM (20 years ago)

Author:

ole

Message:

Bug huntin'

Location:

inundation/ga/storm_surge

Files:

• examples/run_tsh.py (modified) (2 diffs)
• pyvolution/config.py (modified) (1 diff)
• pyvolution/issues.txt (modified) (3 diffs)
• pyvolution/netherlands.py (modified) (4 diffs)
• pyvolution/shallow_water.py (modified) (1 diff)
• pyvolution/shallow_water_ext.c (modified) (1 diff)

inundation/ga/storm_surge/examples/run_tsh.py

@@ -50 +50 @@
 
     domain.checkpoint = False #True
-    domain.default_order = 2
+    domain.default_order = 1
     domain.visualise = visualise
 
@@ -75 +75 @@
 
     #Set friction
-    #manning = 0.07
-    manning = 0.3    
+    manning = 0.07
     inflow_stage = 10.0
     domain.set_quantity('friction', manning)
+
+    #FIXME: Make flat
+    #domain.set_quantity('elevation', 0.0)
+    #domain.set_quantity('level', 0.0)    
 
     ######################

inundation/ga/storm_surge/pyvolution/config.py

@@ -2 +2 @@
 """
 
-epsilon = 1.0e-13
+epsilon = 1.0e-12
 
 default_boundary_tag = 'exterior'

inundation/ga/storm_surge/pyvolution/issues.txt

@@ -1 +1 @@
 CURRENT PYVOLUTION ISSUES
+
 
 OPEN - issues that have not beeen resolved
@@ -18 +19 @@
 to == the bed elevation, if level < bed elevation.  This will make
 setting things up easier.
+
+
+Issue: Transmissive boundary condition can cause unphysical momentums if
+material is not moving out of the domain.
+Importance: Med
+Suggested Action: Add some sort of warning in (subclass of) Transmissive_boundary 
+alerting inflow of material.
 
 Issue:  Run-time error;
@@ -106 +114 @@
 _________________________________________________________
 
+
 ****************
 CLOSED - issues that have  beeen resolved
 ************
+
+
+New example 26/10 16:20
+
+V:\7\cit\workspaces\ws_nielsen\projects\storm_surge\examples>python run_
+ cns_ground_pts_triangulation_wave.tsh 10 500
+Creating domain from cns_ground_pts_triangulation_wave.tsh
+Number of triangles =  4962
+Output being written to C:\grohm_output\cns_ground_pts_triangulation_wav
+_ys10.0_ft500.0.sww
+Boundaries
+Available boundary tags are ['external', 'wave']
+Time = 0.0000, steps=0 (0)
+Time = 10.0000, delta t in [0.00018562, 0.47635714], steps=296 (0)
+Time = 20.0000, delta t in [0.00058530, 0.37116719], steps=267 (0)
+Time = 30.0000, delta t in [0.00014949, 0.22819244], steps=282 (0)
+Time = 40.0000, delta t in [0.00000000, 0.15756478], steps=478 (0)
+Time = 50.0000, delta t in [0.00000000, 0.32346737], steps=648 (1)
+Time = 60.0000, delta t in [0.00000000, 0.24878640], steps=664 (0)
+Time = 70.0000, delta t in [0.00000000, 0.08966319], steps=819 (0)
+Time = 80.0000, delta t in [0.00000000, 0.05706256], steps=1518 (1)
+Time = 90.0000, delta t in [0.00000000, 0.06957202], steps=2129 (3)
+Time = 100.0000, delta t in [0.00000000, 0.08955953], steps=3644 (18)
+Time = 110.0000, delta t in [0.00000000, 0.05102101], steps=2120 (0)
+Time = 120.0000, delta t in [0.00000000, 0.07261881], steps=2423 (0)
+Time = 130.0000, delta t in [0.00000000, 0.07346308], steps=2366 (1)
+Time = 140.0000, delta t in [0.00000000, 0.17995247], steps=2825 (5)
+Time = 150.0000, delta t in [0.00000000, 0.08268458], steps=1477 (0)
+Time = 160.0000, delta t in [0.00000000, 0.07423284], steps=3631 (0)
+Time = 170.0000, delta t in [0.00000000, 0.07305694], steps=4369 (31)
+Time = 180.0000, delta t in [0.00000000, 0.06342947], steps=3235 (0)
+Time = 190.0000, delta t in [0.00000000, 0.10134596], steps=2235 (0)
+Time = 200.0000, delta t in [0.00000000, 0.06387375], steps=4239 (10)

inundation/ga/storm_surge/pyvolution/netherlands.py

@@ -93 +93 @@
 N = 600 #Size = 720000
 N = 20
+#N = 150
 N = 130
 
@@ -104 +105 @@
     
 domain.check_integrity()
-domain.default_order = 1
-domain.smooth = True
+domain.default_order = 2
+
+#Output params
+domain.smooth = False
 domain.reduction = min  #Looks a lot better on top of steep slopes
 
@@ -150 +153 @@
 
 
-#Time dependent inflow
-from math import sin, pi, exp
-Bw = Time_boundary(domain=domain,
-                   f=lambda t: [0.7*(1 + sin(t*pi/2.5))*\
-                                (inflow_stage*(sin(t*pi)+0.7)),0,0])
-
 #Set boundary conditions
-#domain.set_boundary({'left': Bw, 'right': Br, 'bottom': Br, 'top': Br})
-domain.set_boundary({'left': Bd, 'right': Bt, 'bottom': Bt, 'top': Bt})
+domain.set_boundary({'left': Bd, 'right': Br, 'bottom': Br, 'top': Br})
+#domain.set_boundary({'left': Bd, 'right': Bt, 'bottom': Bt, 'top': Bt})
 
                   
@@ -172 +169 @@
 t0 = time.time()
 
-for t in domain.evolve(yieldstep = 0.1, finaltime = 2.0):
+for t in domain.evolve(yieldstep = 0.1, finaltime = 5.0):
     domain.write_time()
     

inundation/ga/storm_surge/pyvolution/shallow_water.py

@@ -450 +450 @@
     balance_deep_and_shallow(domain)
 
-    #Protect against infinitesimal depths at vertices    
-    ####dry(domain)
-    
     #Compute edge values by interpolation
     for name in domain.conserved_quantities:

inundation/ga/storm_surge/pyvolution/shallow_water_ext.c

@@ -226 +226 @@
       alpha = 1.0;  //Flat bed
 
-
+      
     //Weighted balance between stage parallel to bed elevation 
     //(wvi = zvi + hc) and stage as computed by 1st or 2nd