Changeset 4386 for anuga_work

Timestamp:

Apr 17, 2007, 2:23:21 PM (18 years ago)

Author:

duncan

Message:

new UA friction scenario - elevation calculated by a function, initial stage of 0.06

Location:

anuga_work/development/friction_UA_flume_2006

Files:

• create_mesh.py (modified) (2 diffs)
• parallel_run_dam.py (modified) (4 diffs)

anuga_work/development/friction_UA_flume_2006/create_mesh.py

@@ -7 +7 @@
 
 
-def generate(mesh_filename, is_coarse = False):
+def generate(mesh_filename, maximum_triangle_area=0.01):
     """
     Generate mesh for University of Aberbeen dam break flume.
@@ -68 +68 @@
     # this is the location of the slope region.
     slope.setTag("slope")
+    #slope.setTag("dam")
     
-    if is_coarse:
-        m.generate_mesh(maximum_triangle_area=0.01)
-    else:
-        m.generate_mesh(maximum_triangle_area=0.0001)
+    m.generate_mesh(maximum_triangle_area=maximum_triangle_area)
 
     m.export_mesh_file(mesh_filename)

anuga_work/development/friction_UA_flume_2006/parallel_run_dam.py

@@ -65 +65 @@
         yieldstep = 1
         finaltime = 15.1
+        maximum_triangle_area=0.01
     else:
         yieldstep = 0.02
         finaltime = 15.1
+        maximum_triangle_area=0.0001
+        
+        maximum_triangle_area=0.001
         
 
@@ -107 +111 @@
     # Create the triangular mesh
     #-------------------------------------------------------------------------
-    
+
+    # this creates the mesh
     #gate_position = 12.0
     create_mesh.generate(mesh_filename,
-                         is_coarse=is_trial_run) # this creates the mesh
+                         maximum_triangle_area=maximum_triangle_area)
 
     head,tail = path.split(mesh_filename)
@@ -157 +162 @@
 
     for t in domain.evolve(yieldstep, finaltime):
+    
         domain.write_time()
-    
         print 'That took %.2f seconds' %(time.time()-t0)
         print 'finished'
@@ -199 +204 @@
                  0.019, 0.02,
                  0.0225, 0.025, 0.03, 0.035, 0.04] # 17 processes
-    #frictions = [0.0]
+    frictions = [0.0, 0.014]
     for i, friction in enumerate(frictions):
         if i%size == rank:
             print "I am processor %d of %d on node %s" %(rank, size, node)
-            main(friction, is_trial_run = True, outputdir_name='friction_UA_new_meshII')
+            main(friction, is_trial_run = True,
+                 outputdir_name='friction_UA_new_meshII_0.06_stageIII')