Changeset 8359 for trunk/anuga_work/development/gareth/tests

Timestamp:

Mar 13, 2012, 3:46:12 PM (13 years ago)

Author:

davies

Message:

balanced_dev: important change to 'protect'

Location:

trunk/anuga_work/development/gareth/tests

Files:

• channel_floodplain/plotme.py (modified) (1 diff)
• dam_break/plotme.py (modified) (1 diff)
• runup_sinusoid/runup_sinusoid.py (modified) (5 diffs)
• shallow_steep_slope/channel_SU_sparse.py (modified) (1 diff)

trunk/anuga_work/development/gareth/tests/channel_floodplain/plotme.py

@@ -4 +4 @@
 
 # Time-index to plot outputs from
-index=1100
+index=75
 
 #p2 = util.get_output('channel_floodplain1_bal_dev.sww', minimum_allowed_height=0.01)

trunk/anuga_work/development/gareth/tests/dam_break/plotme.py

@@ -9 +9 @@
 p2_st=util.get_centroids(p_st)
 
-p_dev = util.get_output('dam_break_20120308_153245/dam_break.sww', 0.001)
+p_dev = util.get_output('dam_break_20120313_153101/dam_break.sww', 0.001)
 p2_dev=util.get_centroids(p_dev, velocity_extrapolation=True)
 

trunk/anuga_work/development/gareth/tests/runup_sinusoid/runup_sinusoid.py

@@ -18 +18 @@
 #Setup computational domain
 #---------
-points, vertices, boundary = anuga.rectangular_cross(40,40)
+points, vertices, boundary = anuga.rectangular_cross(40,40, len1=1., len2=1.)
 
 domain=Domain(points,vertices,boundary)    # Create Domain
@@ -25 +25 @@
 domain.set_quantities_to_be_stored({'stage': 2, 'xmomentum': 2, 'ymomentum': 2, 'elevation': 1})
 #domain.set_store_vertices_uniquely(True)
+domain.minimum_allowed_height=0.001
 
 #------------------
@@ -30 +31 @@
 #------------------
 scale_me=1.0
+
+#domain.minimum_allowed_height=domain.minimum_allowed_height*scale_me # Seems needed to make the algorithms behave
+
 def topography(x,y):
         return (-x/2.0 +0.05*numpy.sin((x+y)*50.0))*scale_me
@@ -67 +71 @@
 # Associate boundary tags with boundary objects
 #----------------------------------------------
-domain.set_boundary({'left': Br, 'right': Bd, 'top': Br, 'bottom':Br})
+domain.set_boundary({'left': Br, 'right': Br, 'top': Br, 'bottom':Br})
 
 #------------------------------
@@ -75 +79 @@
 for t in domain.evolve(yieldstep=0.1,finaltime=20.0):
     print domain.timestepping_statistics()
+    xx = domain.quantities['xmomentum'].centroid_values
+    yy = domain.quantities['ymomentum'].centroid_values
+    dd = domain.quantities['stage'].centroid_values - domain.quantities['elevation'].centroid_values 
+    dd = (dd)*(dd>1.0e-03)+1.0e-06
+    vv = ( (xx/dd)**2 + (yy/dd)**2 )**0.5
+    vv = vv*(dd>1.0e-03)
+    print 'Peak velocity is: ', vv.max(), vv.argmax()
 
 print 'Finished'

trunk/anuga_work/development/gareth/tests/shallow_steep_slope/channel_SU_sparse.py

@@ -57 +57 @@
 
 Br = anuga.Reflective_boundary(domain) # Solid reflective wall
-Bo= anuga.Dirichlet_boundary([-9.9369037,0.0,0]) # Outflow for steady uniform flow with a depth integrated velocity of 0.1
+#Bo= anuga.Dirichlet_boundary([-9.9369037,0.0,0]) # Outflow for steady uniform flow with a depth integrated velocity of 0.1
+
+def constant_water(t):
+    return -9.936903
+
+#Bo= anuga.shallow_water.boundaries.Transmissive_n_momentum_zero_t_momentum_set_stage_boundary(domain, constant_water) # Outflow for steady uniform flow with a depth integrated velocity of 0.1
+Bo = anuga.Transmissive_boundary(domain)
 domain.set_boundary({'left': Br, 'right': Bo, 'top': Br, 'bottom': Br})
 #------------------------------------------------------------------------------