Changeset 5659 for anuga_work/development/Hinwood_2008/calc_norm.py

Timestamp:

Aug 15, 2008, 11:51:01 AM (16 years ago)

Author:

duncan

Message:

Current Hinwood - Removed sensor 8 from last 4 runs

File:

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

anuga_work/development/Hinwood_2008/calc_norm.py

@@ -127 +127 @@
                                points,
                                rmsd_list))
-        
+
+
+                  
 def plot_rrms_sensor_settings(run_data, outputdir_tags, quantity,
               save_as=None,
-              is_interactive=False):
+              is_interactive=False, max_rmsd=None):
     """
     For a scenario, do 
@@ -137 +139 @@
          savefig, title, axis, setp, subplot, grid, axvspan, figlegend
     from anuga.shallow_water.data_manager import csv2dict
-    
+
+    # TODO
+    # scale the plot
     plot_type = ".pdf"
     
@@ -171 +175 @@
         lines.append(plot(locations, rmsd_list))
         
+        if max_rmsd is not None:
+            #print "setting axis"
+            axis(ymin=0, ymax=max_rmsd)
 
     for break_x in run_data['break_xs']:
@@ -257 +264 @@
     return err, amount
 
-    
+
+def calc_max_rmsd(scenarios, outputdir_tags, quantity):
+    
+    max_rmsd = 0
+    
+    for run_data in scenarios:
+        id = run_data['scenario_id']
+        for outputdir_tag in outputdir_tags:
+            
+            outputdir_name = 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)
+            rmsd_list = [float(x) for x in simulation['rmsd']]
+            max_rmsd = max(max(rmsd_list), max_rmsd)
+    return max_rmsd
+
         
 def auto_plot_rrms_sensor_settings(outputdir_tags, scenarios, quantity,
                    save_as=None,
                    is_interactive=False):
-        
+    max_rmsd = calc_max_rmsd(scenarios, outputdir_tags, quantity)
+    print "max_rmsd", max_rmsd
     for run_data in scenarios:
-        plot_rrms_sensor_settings(run_data, outputdir_tags, quantity) 
+        plot_rrms_sensor_settings(run_data, outputdir_tags, quantity,
+                                  max_rmsd=max_rmsd) 
                  
 #-------------------------------------------------------------
@@ -281 +310 @@
     outputdir_tags.append("_nolmts_wdth_0.1_z_0.012_ys_0.01_mta_0.001_G")
     outputdir_tags.append("_nolmts_wdth_0.01_z_0.012_ys_0.01_mta_1e-05_G")
-    outputdir_tags = [] # _no_velocity_nolmts_wdth_0.1_z_0.012_ys_0.01_mta_0.01_G
+    outputdir_tags = []
     outputdir_tags.append("_nolmts_wdth_0.1_z_0.0_ys_0.01_mta_0.01_G")
     
@@ -289 +318 @@
     #scenarios = scenarios[4:] # !!!!!!!!!!!!!!!!!!!!!!
     #for outputdir_tag in outputdir_tags:
-     #   auto_rrms(outputdir_tag, scenarios, "stage", y_location_tag=':0.0')
+    #    auto_rrms(outputdir_tag, scenarios, "stage", y_location_tag=':0.0')
     
     #scenarios = [scenarios[0]] # !!!!!!!!!!!!!!!!!!!!!!