Changeset 5664 for anuga_work/development/Hinwood_2008/plot.py

Timestamp:

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

Author:

duncan

Message:

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

File:

• anuga_work/development/Hinwood_2008/plot.py (modified) (5 diffs)

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