Changeset 5155 for anuga_work/development/attenuation_near_shore

Timestamp:

Mar 11, 2008, 10:16:29 AM (17 years ago)

Author:

duncan

Message:

newest

File:

• anuga_work/development/attenuation_near_shore/run_beach.py (modified) (3 diffs)

anuga_work/development/attenuation_near_shore/run_beach.py

@@ -139 +139 @@
     domain.set_default_order(2) # Use second order spatial scheme
     domain.set_timestepping_method('rk2')
-    domain.beta_w      = 1.5
     domain.beta_uh     = 1.5
     domain.beta_vh     = 1.5
-    domain.beta_h      = 1.5
+    # comment out aserts in shallow_water/shallow_water_domain.py", line 402
+    # if you want to use these values
+    #domain.beta_h      = 1.5
+    #domain.beta_w      = 1.5
     domain.set_quantities_to_be_stored(['stage', 'xmomentum', 'ymomentum'])
     domain.set_minimum_storable_height(0.001)
@@ -165 +167 @@
     Bwp = Transmissive_Momentum_Set_Stage_boundary(domain,
                    lambda t: [(1.2*sin(2*t*pi/10)), 0.0 ,0.0])
-    Bwp_velocity = Time_boundary(domain,
-                   lambda t: [(1.2*sin(2*t*pi/10)),
-                              0.24*sin(2*t*pi/10),0.0])
-    Bws = Transmissive_Momentum_Set_Stage_boundary(domain,
-                   lambda t: [1.2*sin(2*t*pi/10)+1.0*sin(2*t*pi/9.5),
-                              0.0, 0.0])
-    domain.set_boundary( {'wall': Br, 'wave': Bwp_velocity} )
+    domain.set_boundary( {'wall': Br, 'wave': Bwp} )
 
     #-------------------------------------------------------------------------
@@ -193 +189 @@
 if __name__ == "__main__":
     main( is_trial_run = False, friction=0.0,
-         outputdir_name='Velocity_wave_2nd_order')
+         outputdir_name='Primary_wave_2nd_order_rk2_wet_betas_1.5')