Changeset 5449

Timestamp:

Jun 30, 2008, 9:35:07 AM (17 years ago)

Author:

duncan

Message:

Current Hinwood scenario

Location:

anuga_work/development/Hinwood_2008

Files:

• plot.py (modified) (1 diff)
• run_dam.py (modified) (2 diffs)
• slope.py (modified) (2 diffs)

anuga_work/development/Hinwood_2008/plot.py

@@ -77 +77 @@
     return file_exp, file_sim, location_sims, location_exps, outputdir_name, \
            save_as_list
-
+def plot():
+    from scenarios import scenarios
+    outputdir_tag = "_good_tri_area_0.01_A"
+    outputdir_tag = "_test"
+    #scenarios = scenarios[1] # !!!!!!!!!!!!!!!!!!!!!! 
+    #scenarios = [scenarios[5]] # !!!!!!!!!!!!!!!!!!!!!! 
+    # Supported formats: emf, eps, pdf, png, ps, raw, rgba, svg, svgz.
+    plot_type = ".pdf"
+    for run_data in scenarios:
+        
+        temp = Hinwood_files_locations(run_data, outputdir_tag, plot_type)
+        file_exp, file_sim, location_sims, location_exps, outputdir_name, \
+                  save_as_list = temp
+        print "run_data['scenario_id']", run_data['scenario_id']
+        #location_sims = [location_sims[0]]
+        #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']):
+            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
+                         
+            print "Doing ", plot_title
+            plot_compare_csv(location_sim=loc_sim,
+                             file_sim=file_sim,
+                             location_exp=loc_exp,
+                             file_exp=file_exp,
+                             plot_title=plot_title,
+                             y_label='Water height (m)',
+                             is_interactive=False,
+                             save_as=save_as)
+    
 #-------------------------------------------------------------
 if __name__ == "__main__":

anuga_work/development/Hinwood_2008/run_dam.py

@@ -250 +250 @@
     
     from scenarios import scenarios
+    from slope import gauges_for_slope
+    #from plot import plot
 
 
     #4 is 0.01 5 is 0.001
-    run_type = 1
+    run_type = 5
     #for run_data in [scenarios[5]]:
     for run_data in scenarios:
@@ -259 +261 @@
               run_type = run_type,
               outputdir_name=run_data['scenario_id'])
-    
+    gauges_for_slope()
+    #plot()

anuga_work/development/Hinwood_2008/slope.py

@@ -104 +104 @@
     
     outputdir_tag = "_good_tri_area_0.01_A"
+    outputdir_tag = "_test"
     scenarios = [scenarios[1]] # !!!!!!!!!!!!!!!!!!!!!!
     for run_data in scenarios:
@@ -117 +118 @@
 def gauges_for_slope():
     from scenarios import scenarios
-    from scenarios import scenarios
     
     outputdir_tag = "_good_tri_area_0.01_A"
-    scenarios = [scenarios[1]] # !!!!!!!!!!!!!!!!!!!!!!
+    outputdir_tag = "_test"
+    #scenarios = [scenarios[1]] # !!!!!!!!!!!!!!!!!!!!!!
     dx = 0.005
     for run_data in scenarios: