Changeset 5192

Timestamp:

Apr 7, 2008, 3:55:51 PM (17 years ago)

Author:

duncan

Message:

near shore work, new limiters primary function, with build-up phase

File:

• anuga_work/development/near_shore_PMD/run_beach.py (modified) (6 diffs)

anuga_work/development/near_shore_PMD/run_beach.py

@@ -77 +77 @@
 
     # Structured mesh variables
-    dx = 5.
+    dx = 1.0
     dy = dx
     L = 2080.
@@ -87 +87 @@
         outputdir_name += '_test'
         yieldstep = 1
-        finaltime = 100.
+        finaltime = 200.
         maximum_triangle_area=3000
         thinner=True
     else:
         yieldstep = 1.
-        finaltime = 1600
+        finaltime = 1700 #1600
         maximum_triangle_area = 100
         thinner=True
@@ -193 +193 @@
     Br = Reflective_boundary(domain)
     Bd = Dirichlet_boundary([10.,0,0])
+
+    
+    def wave_func(t):
+        if t < 100:
+            wave = [(1.2*t/100.0*sin(2*t*pi/10)), 0.0 ,0.0]
+        else:
+            wave = [(1.2*sin(2*t*pi/10)), 0.0 ,0.0]           
+        return wave
+    
     Bwp = Transmissive_Momentum_Set_Stage_boundary(domain,
-                   lambda t: [(1.2*sin(2*t*pi/10)), 0.0 ,0.0])
+                   wave_func)
     Bwp_velocity = Time_boundary(domain,
                    lambda t: [(1.2*sin(2*t*pi/10)),
@@ -206 +215 @@
                               'bottom': Br ,
                               'left' : Br ,
-                              'right' : Bwp } )
+                              'right' : Bwp} )
     else:
         domain.set_boundary( {'wall': Br, 'wave': Bwp} )
@@ -218 +227 @@
     for t in domain.evolve(yieldstep, finaltime):    
         domain.write_time()
-
+        
     print 'That took %.2f seconds' %(time.time()-t0)
     print 'finished'
@@ -247 +256 @@
         friction=0.0,
         is_structured=True,
-        outputdir_name='structured_old_limiters')
+        outputdir_name='structured_ramp_up2')