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

Timestamp:

Aug 26, 2008, 9:18:42 AM (15 years ago)

Author:

duncan

Message:

moving Honwood files around

File:

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

anuga_work/development/Hinwood_2008/plot.py

@@ -156 +156 @@
                                      y_location_tag=':0.0',
                                      add_run_info=False):
+    """
+    Used by validation graphs to produce the final figures
+    """
     plot_type = ".pdf"
     save_as_list = []
@@ -180 +183 @@
         #save_as = pro_instance.plots_dir + sep + \
         #          outputdir_name + "_" + quantity + "_" + plot_type
-        save_as = pro_instance.plots_dir + sep + id + quantity + \
-                  '_compare' + plot_type
+        save_as = pro_instance.plots_dir + sep + 'S' + id[1:2] + \
+                  '-' + quantity + '-compare' + plot_type
         save_as_list.append(save_as)
 
@@ -210 +213 @@
                                 y_axis=None,
                                 y_location_tag=':0.0'):
-    """ 
+    """
+    Used by validation graphs to produce the final figures
     """
     from pylab import ion, plot, xlabel, ylabel, close, legend, \
          savefig, title, axis, setp, subplot, grid, figlegend, gca, \
-         text
+         text, sqrt
+    import pylab
     
     from anuga.shallow_water.data_manager import csv2dict
     
     print "file_sim", file_sim
+
+    if False:
+        fig_width_pt = 246.0  # Get this from LaTeX using \showthe\columnwidth
+        inches_per_pt = 1.0/72.27               # Convert pt to inch
+        golden_mean = (sqrt(5)-1.0)/2.0         # Aesthetic ratio
+        fig_width = fig_width_pt*inches_per_pt  # width in inches
+        fig_height = fig_width*golden_mean      # height in inches
+        fig_size =  [fig_width,fig_height]
+        print "fig_size", fig_size
+        params = {'backend': 'ps',
+                  'axes.labelsize': 10,
+                  'text.fontsize': 10,
+                  'legend.fontsize': 10,
+                  'xtick.labelsize': 8,
+                  'ytick.labelsize': 8,
+                  'text.usetex': True,
+                  'figure.figsize': [3.4, 2.1]}
+    if False:
+        params = {'backend': 'ps',
+                  'axes.labelsize': 10,
+                  'text.fontsize': 10,
+                  'legend.fontsize': 10,
+                  'xtick.labelsize': 8,
+                  'ytick.labelsize': 8,
+                  'figure.figsize': [3.4, 3.1]}
+        pylab.rcParams.update(params)
+
+
+    
     # Load in the csv files and convert info from strings to floats
     simulation, _ = csv2dict(file_sim)
@@ -231 +265 @@
     for position, i in enumerate(gauge_indexs):
         gauge_x = run_data['gauge_x'][i]
-        grid_position = (len(gauge_indexs)+1)*100 + 10 + position +1
+        
+        grid_position = (len(gauge_indexs))*100 + 10 + position +1
         subplot(grid_position)
         location_sim = str(gauge_x) + y_location_tag
@@ -241 +276 @@
     
         l_sim, l_exp = plot(time_sim, quantity_sim, time_exp, quantity_exp)
-        setp(l_sim, color='r', linestyle='--')
-        setp(l_exp, color='b')
+        setp(l_sim, color='k', linestyle='--')
+        setp(l_exp, color='k')
         grid(True)
         
@@ -282 +317 @@
     # The order defines the label
     #legend((legend_exp, legend_sim),'upper left')
-    figlegend((l_sim, l_exp),
-              (legend_sim, legend_exp),
-              'lower center')
+
+    # I couldn't get this working
+    #figlegend((l_sim, l_exp),
+     #         (legend_sim, legend_exp),
+      #        'lower center', axespad = -0.05) 
     
     if is_interactive: