Changeset 5577 for anuga_work/development/Hinwood_2008/run_dam.py

Timestamp:

Jul 25, 2008, 4:36:05 PM (15 years ago)

Author:

duncan

Message:

Current Hinwood scenario calculating RMSD's

File:

• anuga_work/development/Hinwood_2008/run_dam.py (modified) (8 diffs)

anuga_work/development/Hinwood_2008/run_dam.py

@@ -62 +62 @@
          outputdir_name=None,
          run_type=0,
+         width=1.0,
+         use_limits=True,
          end_tag = '_limiterD'):
 
@@ -68 +70 @@
     
     if run_type == 1:
-        outputdir_name += '_test' + end_tag
         yieldstep = 1.0
         finaltime = 15.
         maximum_triangle_area=0.1
+        outputdir_name += '_test'
         
     elif run_type == 2:
-        outputdir_name += '_test_long_time' + end_tag
         yieldstep = 0.5
         finaltime = None
         maximum_triangle_area=0.01
+        outputdir_name += '_test_long_time'
         
     elif run_type == 3:
-        outputdir_name += '_yieldstep_0.1_tri_area_0.01' + end_tag
         yieldstep = 0.1
         finaltime = None       
         maximum_triangle_area=0.01
+        #outputdir_name += '_yieldstep_0.1'
+        
     elif run_type == 4:
-        outputdir_name += '_good_tri_area_0.01' + end_tag
         # this is not a test
         # Output will go to a file
@@ -92 +94 @@
         finaltime = None        
         maximum_triangle_area=0.01
+        #outputdir_name += '_good'
+        
     elif run_type == 5:
-        outputdir_name += '_good_tri_area_0.001' + end_tag
         # this is not a test
         # Output will go to a file
@@ -100 +103 @@
         finaltime = None        
         maximum_triangle_area=0.001
-     
+        #outputdir_name += '_good'
+        
+    elif run_type == 6:
+        # this is not a test
+        # Output will go to a file
+        # The sww file will be interpolated
+        yieldstep = 0.01
+        finaltime = None        
+        maximum_triangle_area=0.0001
+        #outputdir_name += '_good'
+        
+    if use_limits is True:
+        outputdir_name += '_lmts'
+    else:
+        outputdir_name += '_nolmts'
+    outputdir_name += '_wdth_' + str(width)
+    outputdir_name += '_z_' + str(friction)
+    outputdir_name += '_ys_' + str(yieldstep)
+    outputdir_name += '_mta_' + str(maximum_triangle_area)
+    outputdir_name += end_tag
+    
     metadata_dic = set_z_origin_to_water_depth(metadata_dic)    
         
@@ -150 +173 @@
     # this creates the mesh
     #gate_position = 12.0
-    create_mesh.generate(mesh_filename, metadata_dic,
+    create_mesh.generate(mesh_filename, metadata_dic, width=width,
                          maximum_triangle_area=maximum_triangle_area)
 
@@ -172 +195 @@
     domain.set_minimum_storable_height(0.0001)
 
-    domain.set_default_order(2) # Use second order spatial scheme
-    domain.set_timestepping_method('rk2')
-    domain.use_edge_limiter = True
-    domain.tight_slope_limiters = True
-    
-    domain.beta_w      = 0.6
-    domain.beta_uh     = 0.6
-    domain.beta_vh     = 0.6
+    if use_limits is True:
+        domain.set_default_order(2) # Use second order spatial scheme
+        domain.set_timestepping_method('rk2')
+        domain.use_edge_limiter = True
+        domain.tight_slope_limiters = True
+        
+        domain.beta_w      = 0.6
+        domain.beta_uh     = 0.6
+        domain.beta_vh     = 0.6
     
 
@@ -223 +247 @@
     print 'finished'
 
-    flume_y_middle = 0.5
+    flume_y_middle = 0.0
     points = []
     for gauge_x in metadata_dic['gauge_x']:
@@ -259 +283 @@
     #from plot import plot
 
-
+    # 1 is fast and dirty
     # 4 is 0.01
     # 5 is 0.001
+    # 6 is 0.0001
+    
+    #run_type = 1
     run_type = 5
     #for run_data in [scenarios[5]]:
     #scenarios = scenarios[2:]
     #scenarios = [scenarios[0]]
+    width = 1.0
+    width = 0.1
     for run_data in scenarios:
         pro_instance = main( run_data['scenario_id'] + '_boundary.tsm'  ,
                              run_data,
-                             run_type = run_type,
+                             width=width,
+                             run_type=run_type,
                              outputdir_name=run_data['scenario_id'],
-                             end_tag='_limiterD')
+                             use_limits=False,
+                             friction=0.012, 
+                             end_tag='_G')
         #gauges_for_slope(pro_instance.outputdir,[run_data])