Changeset 5664 for anuga_work/development/Hinwood_2008

Timestamp:

Aug 18, 2008, 9:29:57 AM (16 years ago)

Author:

duncan

Message:

Current Hinwood - calc norm set up for _G end tag comparison

Location:

anuga_work/development/Hinwood_2008

Files:

• calc_norm.py (modified) (4 diffs)
• plot.py (modified) (5 diffs)
• run_dam.py (modified) (2 diffs)
• scenarios.py (modified) (8 diffs)

anuga_work/development/Hinwood_2008/calc_norm.py

@@ -131 +131 @@
                   
 def plot_rrms_sensor_settings(run_data, outputdir_tags, quantity,
-              save_as=None,
+              save_as=True,
               is_interactive=False, max_rmsd=None):
     """
@@ -140 +140 @@
     from anuga.shallow_water.data_manager import csv2dict
 
-    # TODO
-    # scale the plot
+
     plot_type = ".pdf"
     
@@ -196 +195 @@
     close()
 
+def auto_plot_test_rmsd(scenarios, outputdir_tag, quantity="stage",
+              save_as=True,
+              is_interactive=False, max_rmsd=1.0):
+    
+    tests = []
+    
+    for i in range(0,len(scenarios),2):
+        tests.append([scenarios[i], scenarios[i+1]])
+
+    for test in tests:
+        plot_test_rmsd(test, outputdir_tag, quantity,
+              save_as=save_as,
+              is_interactive=is_interactive, max_rmsd=max_rmsd)
+        
+def plot_test_rmsd(test, outputdir_tag, quantity,
+              save_as=None,
+              is_interactive=False, max_rmsd=None):
+    """
+    For a test, plot rmsd vs x location for both runs.
+    Also, add the break measurement.
+    """
+    from pylab import ion, plot, xlabel, ylabel, close, legend, \
+         savefig, title, axis, setp, subplot, grid, axvspan, figlegend
+    from anuga.shallow_water.data_manager import csv2dict
+        
+    plot_type = ".pdf"   
+    id = test[0]['scenario_id'][:2]
+    
+    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" + '\n' \
+                 + time_date
+                 
+    #title(plot_title)
+    y_label = "RMSD"
+    ylabel(y_label)
+    xlabel("x location, m")
+    grid(True)
+
+    lines = []
+    for run in test:
+            
+        outputdir_name = run['scenario_id'] + outputdir_tag
+        pro_instance = project.Project(['data','flumes','Hinwood_2008'],
+                                       outputdir_name=outputdir_name)
+        
+        file_err = pro_instance.rmsd_dir + sep + outputdir_name + "_" \
+                   + quantity + "_err.csv"
+            
+        simulation, _ = csv2dict(file_err)
+        locations = [float(x) for x in simulation['x location']]
+        rmsd_list = [float(x) for x in simulation['rmsd']]
+        lines.append(plot(locations, rmsd_list))
+        
+        for break_x in run['break_xs']:
+            axvspan(break_x-0.001,break_x+0.001, facecolor='g')
+        
+    if max_rmsd is not None:
+        #print "setting axis"
+        axis(ymin=0, ymax=max_rmsd)
+       
+    figlegend(lines, outputdir_tags,'upper left')
+    
+    if is_interactive:
+        raw_input()
+
+    save_as = pro_instance.plots_dir + sep + \
+              id + "_rmsd" + plot_type
+    if save_as is not None:
+        savefig(save_as)
+    
+    close()
+
 # Return a bunch of lists
 # The err files, for all scenarios
@@ -317 +392 @@
     #scenarios = [scenarios[0]] # !!!!!!!!!!!!!!!!!!!!!!
     #scenarios = scenarios[4:] # !!!!!!!!!!!!!!!!!!!!!!
-    #for outputdir_tag in outputdir_tags:
-    #    auto_rrms(outputdir_tag, scenarios, "stage", y_location_tag=':0.0')
+    if False:
+        for outputdir_tag in outputdir_tags:
+            auto_rrms(outputdir_tag, scenarios, "stage", y_location_tag=':0.0')
     
     #scenarios = [scenarios[0]] # !!!!!!!!!!!!!!!!!!!!!!
-    auto_plot_rrms_sensor_settings(outputdir_tags, scenarios, "stage")
+    #auto_plot_rrms_sensor_settings(outputdir_tags, scenarios, "stage")
     #compare_different_settings(outputdir_tags, scenarios, "stage")
+
+    outputdir_tag = "_nolmts_wdth_0.1_z_0.0_ys_0.01_mta_0.01_G" 
+    auto_plot_test_rmsd(scenarios, outputdir_tag)

anuga_work/development/Hinwood_2008/plot.py

@@ -15 +15 @@
                      legend_exp='Measured flume result',
                      is_interactive=False,
-                     use_axis=None):
+                     x_axis=None,
+                     y_axis=None):
     """ 
     """
@@ -47 +48 @@
     legend((legend_sim, legend_exp),'upper left')
     title(plot_title)
-    if use_axis is not None:
-        axis(use_axis)
+    if x_axis is not None:
+        axis(xmin=x_axis[0], xmax=x_axis[1])
     
+    if y_axis is not None:
+        axis(ymin=y_axis[0], ymax=y_axis[1])
+        
     if is_interactive:
         # Wait for enter pressed
@@ -109 +113 @@
         #location_exps = [location_exps[0]]
         #save_as_list = [save_as_list[0]]
-        for loc_sim, loc_exp, save_as, gauge in map(None, location_sims,
-                                                    location_exps,
-                                                    save_as_list,
-                                                    run_data['gauge_x']):
+        for loc_sim, loc_exp, save_as, gauge, gauge_elev in \
+                map(None, location_sims,
+                    location_exps,
+                    save_as_list,
+                    run_data['gauge_x'],
+                    run_data['gauge_bed_elevation']):
             time_date = strftime('plot date: %d/%m/%Y Time: %H:%M:%S',
                                       localtime())
             plot_title = "Scenario: " + outputdir_name + "\n" + \
                          "X Gauge (m):" + str(gauge) + "    " + time_date
+
+            # When the shore gauges out of the are used, don't
+            # use the standard y axis scale
+            x_axis = run_data['wave_times']
             if gauge < run_data['axis_maximum_x']:
-                use_axis = run_data['axis']
+                y_axis = [run_data['axis'][2],run_data['axis'][3]]
+        
             else:
-                use_axis = None
+                y_axis = [run_data['axis'][2] + gauge_elev,
+                          run_data['axis'][3] + gauge_elev]
             print "Doing ", plot_title
             plot_compare_csv(location_sim=loc_sim,
@@ -130 +142 @@
                              is_interactive=is_interactive,
                              save_as=save_as,
-                             use_axis=use_axis)
+                             x_axis=x_axis,
+                             y_axis=y_axis)
             if is_interactive is True:
                 break
@@ -146 +159 @@
     #outputdir_tag = "_test_C"
     #scenarios = scenarios[1:]
-    #scenarios = [scenarios[4]]
+    #scenarios = [scenarios[7]]
     is_interactive = False
     #is_interactive = True

anuga_work/development/Hinwood_2008/run_dam.py

@@ -310 +310 @@
     
     run_type = 1
-    #run_type = 4
+    run_type = 4
     #for run_data in [scenarios[5]]:
     #scenarios = scenarios[3:]
@@ -324 +324 @@
                              outputdir_name=run_data['scenario_id'],
                              use_limits=False,
-                             friction=0.012,
-                             end_tag='_H')
+                             friction=0.0,
+                             end_tag='_I')
         #gauges_for_slope(pro_instance.outputdir,[run_data])

anuga_work/development/Hinwood_2008/scenarios.py

@@ -36 +36 @@
         'ANUGA_start_time':1.83,
         'band_offset':-1.0,
-        'wave_times':[23.0,59.0]
+        'wave_times':[23.0,59.0]  # this is in Anuga time
                  }
 scenarios.append(data)
@@ -62 +62 @@
         'ANUGA_start_time':0.38,
         'band_offset':-0.5,
-        'wave_times':[20.0,59.0]
+        'wave_times':[20.0,59.0] # this is in Anuga time
         }
 scenarios.append(data)
@@ -90 +90 @@
         'ANUGA_start_time':0.05,
         'band_offset':2.0,
-        'wave_times':[15.0,60.0]
+        'wave_times':[15.0,64.0] # this is in Anuga time
                  }
 scenarios.append(data)
@@ -117 +117 @@
         'ANUGA_start_time':0.04,
         'band_offset':-2.0,
-        'wave_times':[34.0,70.0]
+        'wave_times':[34.0,74.0] # this is in Anuga time
         }
 scenarios.append(data)
@@ -151 +151 @@
         'ANUGA_start_time':12.18,
         'band_offset':-0.5,
-        'wave_times':[30.0,65.0]
+        'wave_times':[30.0,85.0] # this is in Anuga time
         }
 scenarios.append(data)
@@ -184 +184 @@
         'ANUGA_start_time':10.48,
         'band_offset':-2.0,
-        'wave_times':[30.0,65.0]
+        'wave_times':[30.0,74.0] # this is in Anuga time
                  }
 scenarios.append(data)
@@ -218 +218 @@
         'ANUGA_start_time':11.63,
         'band_offset':-0.5,
-        'wave_times':[34.0,70.0]
+        'wave_times':[34.0,75.0] # this is in Anuga time
         }
 scenarios.append(data)
@@ -252 +252 @@
         'ANUGA_start_time':12.68,
         'band_offset':-1.0,
-        'wave_times':[34.0,70.0]
-        }
-scenarios.append(data)
+        'wave_times':[34.0,75.0] # this is in Anuga time
+        }
+scenarios.append(data)