Changeset 4044

Timestamp:

Nov 29, 2006, 2:26:23 PM (18 years ago)

Author:

duncan

Message:

make friction a parameter

File:

• anuga_work/development/dam_2006/run_dam.py (modified) (7 diffs)

anuga_work/development/dam_2006/run_dam.py

@@ -31 +31 @@
 
 
-def main():
+def main(friction):
 
     #--------------------------------------------------------------------------
@@ -47 +47 @@
     myid = 0
     numprocs = 1
-    start_screen_catcher(project.outputtimedir, myid, numprocs)
+    #start_screen_catcher(project.outputtimedir, myid, numprocs)
 
     print 'USER:    ', project.user
@@ -54 +54 @@
     #-------------------------------------------------------------------------
 
-    gate_position = 2.0
+    gate_position = 0.85
     create_mesh.generate(project.mesh_filename,
                          gate_position,
@@ -90 +90 @@
         
     domain.set_quantity('stage', elevation_tilt)
-    domain.set_quantity('friction', 0.03) 
+    domain.set_quantity('friction', friction) 
     domain.set_quantity('elevation',elevation_tilt)
 
     print 'Available boundary tags', domain.get_boundary_tags()
-    domain.set_region('dam','stage',0.15,
+    domain.set_region('dam','stage',0.20,
                                  location = 'unique vertices') 
     #domain.set_region(Set_region('dam','stage',0.2,
@@ -108 +108 @@
     t0 = time.time()
 
-    for t in domain.evolve(yieldstep = 0.05, finaltime = 10):
+    for t in domain.evolve(yieldstep = 0.01, finaltime = 30):
         domain.write_time()
     
@@ -118 +118 @@
               [gate_position - 0.45,0.2],
               [gate_position - 0.35,0.2],
-              [gate_position - 0.25,0.2],
-              [0.55,0.2]
+              [gate_position - 0.25,0.2]
               ]
 
@@ -133 +132 @@
 #-------------------------------------------------------------
 if __name__ == "__main__":
-    main()
+    main(0.000)