 Timestamp:
 Mar 11, 2008, 10:16:29 AM (15 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

anuga_work/development/attenuation_near_shore/run_beach.py
r5102 r5155 139 139 domain.set_default_order(2) # Use second order spatial scheme 140 140 domain.set_timestepping_method('rk2') 141 domain.beta_w = 1.5142 141 domain.beta_uh = 1.5 143 142 domain.beta_vh = 1.5 144 domain.beta_h = 1.5 143 # comment out aserts in shallow_water/shallow_water_domain.py", line 402 144 # if you want to use these values 145 #domain.beta_h = 1.5 146 #domain.beta_w = 1.5 145 147 domain.set_quantities_to_be_stored(['stage', 'xmomentum', 'ymomentum']) 146 148 domain.set_minimum_storable_height(0.001) … … 165 167 Bwp = Transmissive_Momentum_Set_Stage_boundary(domain, 166 168 lambda t: [(1.2*sin(2*t*pi/10)), 0.0 ,0.0]) 167 Bwp_velocity = Time_boundary(domain, 168 lambda t: [(1.2*sin(2*t*pi/10)), 169 0.24*sin(2*t*pi/10),0.0]) 170 Bws = Transmissive_Momentum_Set_Stage_boundary(domain, 171 lambda t: [1.2*sin(2*t*pi/10)+1.0*sin(2*t*pi/9.5), 172 0.0, 0.0]) 173 domain.set_boundary( {'wall': Br, 'wave': Bwp_velocity} ) 169 domain.set_boundary( {'wall': Br, 'wave': Bwp} ) 174 170 175 171 # … … 193 189 if __name__ == "__main__": 194 190 main( is_trial_run = False, friction=0.0, 195 outputdir_name=' Velocity_wave_2nd_order')191 outputdir_name='Primary_wave_2nd_order_rk2_wet_betas_1.5')
Note: See TracChangeset
for help on using the changeset viewer.