Changeset 5590

Timestamp:

Aug 1, 2008, 10:05:08 AM (14 years ago)

Author:

duncan

Message:

Current Hinwood scenario

Location:

anuga_work/development/Hinwood_2008

Files:

• calc_norm.py (modified) (4 diffs)
• run_dam.py (modified) (5 diffs)
• slope.py (modified) (1 diff)

anuga_work/development/Hinwood_2008/calc_norm.py

@@ -14 +14 @@
 from os import sep, path
 from csv import writer
+from time import localtime, strftime
+
+# Related major packages
 from Numeric import arange, array, zeros, Float, where, greater, less, \
      compress, argmin, choose, searchsorted, sqrt, sum
 
-# Related major packages
 import project          
 from os import sep
@@ -142 +144 @@
     if is_interactive:
         ion()
-        
+    time_date = strftime('plot date: %d/%m/%Y Time: %H:%M:%S',
+                                      localtime()) 
     subplot(212)    
-    plot_title = id + " Root Mean Square Deviation comparison"
+    plot_title = id + " Root Mean Square Deviation comparison" + '\n' \
+                 + time_date
+                 
     title(plot_title)
     y_label = "RMSD"
     ylabel(y_label)
+    xlabel("x location, m")
     grid(True)
 
@@ -164 +170 @@
         rmsd_list = [float(x) for x in simulation['rmsd']]
         lines.append(plot(locations, rmsd_list))
+        
+
+    for break_x in run_data['break_xs']:
+        axvspan(break_x-0.001,break_x+0.001, facecolor='g')
         
     figlegend(lines, outputdir_tags,'upper left')
@@ -273 +283 @@
     #scenarios = [scenarios[0]] # !!!!!!!!!!!!!!!!!!!!!!
     #scenarios = scenarios[4:] # !!!!!!!!!!!!!!!!!!!!!!
-    for outputdir_tag in outputdir_tags:
-        auto_rrms(outputdir_tag, scenarios, "stage")
+    #for outputdir_tag in outputdir_tags:
+     #   auto_rrms(outputdir_tag, scenarios, "stage")
     
     #scenarios = [scenarios[0]] # !!!!!!!!!!!!!!!!!!!!!!
-    #auto_plot_rrms(outputdir_tags, scenarios, "stage")
-    plot_settings(outputdir_tags, scenarios, "stage")
+    auto_plot_rrms_sensor_settings(outputdir_tags, scenarios, "stage")
+    #plot_settings(outputdir_tags, scenarios, "stage")

anuga_work/development/Hinwood_2008/run_dam.py

@@ -112 +112 @@
         finaltime = None        
         maximum_triangle_area=0.0001
+        #outputdir_name += '_good'
+        
+    elif run_type == 7:
+        # 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.00001
+        #outputdir_name += '_good'
+        
+    elif run_type == 8:
+        # this is not a test
+        # Output will go to a file
+        # The sww file will be interpolated
+        yieldstep = 0.05
+        finaltime = None        
+        maximum_triangle_area=0.00001
         #outputdir_name += '_good'
         
@@ -232 +250 @@
     Bd = Dirichlet_boundary([0.3,0,0]) 
     Bts = Time_boundary(domain, function)
+    #Bts = Transmissive_Momentum_Set_Stage_boundary(domain, function)
     domain.set_boundary( {'wall': Br, 'wave': Bts} )
     #domain.set_boundary( {'wall': Br, 'wave': Bd} )
@@ -287 +306 @@
     # 5 is 0.001
     # 6 is 0.0001
+    # 7 is 0.00001
     
     #run_type = 1
-    run_type = 5
+    run_type = 8
     #for run_data in [scenarios[5]]:
     #scenarios = scenarios[2:]
@@ -295 +315 @@
     width = 1.0
     width = 0.1
+    #width = 0.01
     for run_data in scenarios:
         pro_instance = main( run_data['scenario_id'] + '_boundary.tsm'  ,
@@ -302 +323 @@
                              outputdir_name=run_data['scenario_id'],
                              use_limits=False,
-                             friction=0.012, 
+                             friction=0.012,
                              end_tag='_G')
         #gauges_for_slope(pro_instance.outputdir,[run_data])

anuga_work/development/Hinwood_2008/slope.py

@@ -674 +674 @@
     #scenarios = [scenarios[0]]
     outputdir_tag = "_good_tri_area_0.01_limiterD"
-    outputdir_tag = "_good_lmts_wdth_0.1_z_0.012_ys_0.01_mta_0.01_F"
+    outputdir_tag = "_good_lmts_wdth_0.1_z_0.012_ys_0.01_mta_0.0001_F"
     slope_tag = ""
     #outputdir_tag = "_test_limiterC"
-    #scenarios = [scenarios[0]] # !!!!!!!!!!!!!!!!!!!!!!
+    scenarios = [scenarios[4]] # !!!!!!!!!!!!!!!!!!!!!!
     #scenarios = scenarios[4:] # !!!!!!!!!!!!!!!!!!!!!!
     
-    gauges_for_slope(slope_tag, outputdir_tag, scenarios)
+    #gauges_for_slope(slope_tag, outputdir_tag, scenarios)
     #auto_graph_slopes(outputdir_tag, scenarios) #, is_interactive=True) 
-    auto_find_min_slopes(slope_tag, outputdir_tag, scenarios)
+    #auto_find_min_slopes(slope_tag, outputdir_tag, scenarios)
     #auto_graph_froudes(outputdir_tag, scenarios) 
-    #auto_plot_froude_slopes(slope_tag, outputdir_tag, scenarios)
+    auto_plot_froude_slopes(slope_tag, outputdir_tag, scenarios)
     #g = Get_file_name(scenarios[0], outputdir_tag, slope_tag)
     #for wave_file, save_as, wave_number in Get_file_name(