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

Timestamp:

Jul 1, 2008, 3:48:00 PM (16 years ago)

Author:

duncan

Message:

Current Hinwood scenario - cropping the flume length

File:

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

anuga_work/development/Hinwood_2008/run_dam.py

@@ -66 +66 @@
     basename = 'zz_' + metadata_dic['scenario_id']
     if run_type == 1:
-        outputdir_name += '_test'
+        outputdir_name += '_test_C'
         yieldstep = 1.0
         finaltime = 15.
@@ -83 +83 @@
         maximum_triangle_area=0.001
     elif run_type == 4:
-        outputdir_name += '_good_tri_area_0.01_A'
+        outputdir_name += '_good_tri_area_0.01_C'
         # this is not a test
         # Output will go to a file
@@ -91 +91 @@
         maximum_triangle_area=0.01
     elif run_type == 5:
-        outputdir_name += '_good_tri_area_0.001_A'
+        outputdir_name += '_good_tri_area_0.001_C'
         # this is not a test
         # Output will go to a file
@@ -98 +98 @@
         finaltime = None        
         maximum_triangle_area=0.001
-    elif run_type == 6:
-        outputdir_name += '_good_tri_area_0.0001_A'
-        # 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
      
     metadata_dic = set_z_origin_to_water_depth(metadata_dic)    
@@ -231 +223 @@
 
     if run_type >= 1:
-        id = metadata_dic['scenario_id']
+        id = metadata_dic['scenario_id'] + ".csv"
         interpolate_sww2csv(pro_instance.outputdir + basename +".sww",
                             points,
-                            pro_instance.outputdir + "depth_" + id + ".csv",
-                            pro_instance.outputdir + "velocity_x_" + id + ".csv",
-                            pro_instance.outputdir + "velocity_y_" + id + ".csv")
+                            pro_instance.outputdir + "depth_" + id,
+                            pro_instance.outputdir + "velocity_x_" + id,
+                            pro_instance.outputdir + "velocity_y_" + id,
+                            pro_instance.outputdir + "stage_" + id)
   
     return pro_instance
@@ -254 +247 @@
 
 
-    #4 is 0.01 5 is 0.001
-    run_type = 5
+    # 4 is 0.01
+    # 5 is 0.001
+    run_type = 1
     #for run_data in [scenarios[5]]:
     for run_data in scenarios:
-        main( run_data['scenario_id'] + '_boundary.tsm'  , run_data,
-              run_type = run_type,
-              outputdir_name=run_data['scenario_id'])
-    gauges_for_slope()
-    #plot()
+        pro_instance = main( run_data['scenario_id'] + '_boundary.tsm'  ,
+                             run_data,
+                             run_type = run_type,
+                             outputdir_name=run_data['scenario_id'])
+        gauges_for_slope(pro_instance.outputdir,[run_data])