Changeset 2395 for development/momentum_sink/run_friction.py

Timestamp:

Feb 14, 2006, 1:43:12 PM (19 years ago)

Author:

nicholas

Message:

Updated compare_sww.py. Increased distance of rear boundary and changed to dirichlet

File:

• development/momentum_sink/run_friction.py (modified) (4 diffs)

development/momentum_sink/run_friction.py

@@ -26 +26 @@
 t0 = time.time()
 
+#meshname, triagle_count =create_mesh(100,mesh_file=meshname,triangles_in_name=True)
+
 meshname, triagle_count = cache(create_mesh,(100),
                                 {'mesh_file':meshname,
@@ -50 +52 @@
 domain.set_datadir(project_friction.outputdir)
 domain.store = True
-domain.quantities_to_be_stored = ['stage']
+domain.quantities_to_be_stored = ['stage','xmomentum','ymomentum']
 
 print 'Number of triangles = ', len(domain)
@@ -58 +60 @@
 domain.set_quantity('friction', 0.01)
 domain.set_quantity('stage', 0)
-domain.set_region(Set_region('mound', 'friction', 0.8)) #, location='unique vertices'))
+domain.set_region(Set_region('mound', 'friction', 3)) #, location='unique vertices'))
 #Setup Boundary Conditions
 print domain.get_boundary_tags()
@@ -66 +68 @@
 Br = Reflective_boundary(domain)
 Bt = Transmissive_boundary(domain)
-Bd = Dirichlet_boundary([36,0,0])
+Bdw = Dirichlet_boundary([36,0,0])
+Bdb = Dirichlet_boundary([0,0,0])
 Bw = Time_boundary(domain=domain,
                    f=lambda t: [(60<t<660)*4, 0, 0])
 
-domain.set_boundary( {'wall': Br,'wave': Bd, 'back': Bt, 'exterior':Bd} ) 
+domain.set_boundary( {'wall': Br,'wave': Bdw, 'back': Bdb, 'exterior':Bdw} ) 
 
 #Evolve
 t0 = time.time()
 
-for t in domain.evolve(yieldstep = 1, finaltime = 80):
+for t in domain.evolve(yieldstep = 4, finaltime = 400):
     domain.write_time()