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

Timestamp:

Aug 6, 2008, 9:36:38 AM (16 years ago)

Author:

duncan

Message:

Current Hinwood scenario

File:

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

anuga_work/development/Hinwood_2008/calc_norm.py

@@ -38 +38 @@
     return band_condition
     
-def auto_rrms(outputdir_tag, scenarios, quantity):
+def auto_rrms(outputdir_tag, scenarios, quantity, y_location_tag=':0.0'):
     """
     Given a bunch of scenarios that have CSV guage files, calc the
@@ -56 +56 @@
         for gauge_x in run_data['gauge_x']:
             gauge_x = str(gauge_x)
-            location_sims.append(gauge_x + ':0.0')
+            location_sims.append(gauge_x + y_location_tag)
             location_exps.append(gauge_x)
         
@@ -236 +236 @@
         sim['rmsd'] = sim['err']/sqrt(sim['amount'])
         settings[outputdir_tag] = sim
-    print "settings", settings
+    #print "settings", settings
 
     aux = [(settings[k]['rmsd'], k) for k in settings.keys()]
@@ -273 +273 @@
     
     outputdir_tags = []
+    #outputdir_tags.append("_nolmts_wdth_0.1_z_0.012_ys_0.01_mta_0.01_H")
+    outputdir_tags.append("_nolmts_wdth_0.1_z_0.0_ys_0.01_mta_0.01_G")
+    #outputdir_tag = "_good_tri_area_0.01_D"
+    outputdir_tags.append("_nolmts_wdth_0.1_z_0.012_ys_0.01_mta_0.01_G")
     outputdir_tags.append("_lmts_wdth_0.1_z_0.012_ys_0.01_mta_0.01_G")
+    outputdir_tags.append("_nolmts_wdth_0.1_z_0.012_ys_0.05_mta_1e-05_G")
+    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.append("_nolmts_wdth_0.1_z_0.0_ys_0.01_mta_0.01_G")
-    outputdir_tags.append("_nolmts_wdth_0.1_z_0.012_ys_0.01_mta_0.001_G")
-    outputdir_tags.append("_nolmts_wdth_0.1_z_0.012_ys_0.01_mta_0.01_G")
-    outputdir_tags.append("_nolmts_wdth_0.01_z_0.012_ys_0.01_mta_1e-05_G")
-    outputdir_tags.append("_nolmts_wdth_0.1_z_0.012_ys_0.05_mta_1e-05_G")
-    outputdir_tags.append("_nolmts_wdth_0.1_z_0.012_ys_0.05_mta_1e-05_G")
-    outputdir_tags = [] # _no_velocity_nolmts_wdth_0.1_z_0.012_ys_0.01_mta_0.01_G
-    outputdir_tags.append("_nolmts_wdth_0.1_z_0.012_ys_0.01_mta_0.01_G")
-    outputdir_tags.append("_no_velocity_nolmts_wdth_0.1_z_0.012_ys_0.01_mta_0.01_G")
-    outputdir_tag = "_good_lmts_wdth_1.0_z_0.012_ys_0.01_mta_0.001_F"
+    
+    #outputdir_tag = "_good_lmts_wdth_1.0_z_0.012_ys_0.01_mta_0.001_H"
     #outputdir_tag = "_test_limiterC"
     #scenarios = [scenarios[0]] # !!!!!!!!!!!!!!!!!!!!!!
     #scenarios = scenarios[4:] # !!!!!!!!!!!!!!!!!!!!!!
     #for outputdir_tag in outputdir_tags:
-    #    auto_rrms(outputdir_tag, scenarios, "stage")
+     #   auto_rrms(outputdir_tag, scenarios, "stage", y_location_tag=':0.0')
     
     #scenarios = [scenarios[0]] # !!!!!!!!!!!!!!!!!!!!!!
-    #auto_plot_rrms_sensor_settings(outputdir_tags, scenarios, "stage")
-    compare_different_settings(outputdir_tags, scenarios, "stage")
+    auto_plot_rrms_sensor_settings(outputdir_tags, scenarios, "stage")
+    #compare_different_settings(outputdir_tags, scenarios, "stage")