From 6e473cfb09b6aa5a0c463aef132cbe5a4f586594 Mon Sep 17 00:00:00 2001
From: othocaeS <xaw5719@wmich.edu>
Date: Tue, 12 Jul 2016 16:23:24 -0400
Subject: [PATCH] working on automating plots

---
 plots_playground/psd_tmp.gp   | 14 +++++++++----
 scripts/templates/psd_freq.gp | 37 +++++++++++++++++++----------------
 2 files changed, 30 insertions(+), 21 deletions(-)

diff --git a/plots_playground/psd_tmp.gp b/plots_playground/psd_tmp.gp
index b851571..de25dc6 100644
--- a/plots_playground/psd_tmp.gp
+++ b/plots_playground/psd_tmp.gp
@@ -3,14 +3,20 @@
 # plot 'file' with errorbars
 # plot 'file' using 1:2:(sqrt($1)) with xerrorbars
 # plot 'file' using 1:2:($1-$3):($1+$3):4:5 with xyerrorbars
+set terminal png
+set output "tmp.png"
 
-plot 'tmp.sourcePSD' using 1:($2):3:4 with xyerrorbars, \
-'tmp.reprocPSD' using 1:($2):3:4 with xyerrorbars,
-#'tmp.ccPSD' using 1:($2):3:4 with xyerrorbars
-#plot 'tmp.timelag' using 1:2:3:4 with xyerrorbars
+set title "Yo"
+set xlabel "Temporal Frequency [1/days]"
+set ylabel "log_{10} PSD"
 set logscale x
 set xrange [0.005:0.603]
 set yrange [-5.5:0.5]
+
+plot 'tmp.sourcePSD' using 1:($2):3:4 with xyerrorbars title "Driving Emission", \
+'tmp.reprocPSD' using 1:($2):3:4 with xyerrorbars title "Reprocessed Emission",
+#'tmp.ccPSD' using 1:($2):3:4 with xyerrorbars
+#plot 'tmp.timelag' using 1:2:3:4 with xyerrorbars
 # set yrange [:1]
 #    (x, y, ydelta),
 #    (x, y, ylow, yhigh),
diff --git a/scripts/templates/psd_freq.gp b/scripts/templates/psd_freq.gp
index d0fed30..fb37deb 100644
--- a/scripts/templates/psd_freq.gp
+++ b/scripts/templates/psd_freq.gp
@@ -1,20 +1,23 @@
+# This method assumes y values are in logarithm units already.
 
+set title "%TITLE"
+set terminal png
+set output "tmp.png"
 
-
-# This method assumes data is in logarithm units already.
-
-plot 'tables/PSD_lcdrive.dat' using 1:(10**($2+2)):3:4 with xyerrorbars, \
-'tables/PSD_lc1.dat' using 1:(10**($2+2)):3:4 with xyerrorbars, \
-'tables/PSD_lc2.dat' using 1:(10**($2+2)):3:4 with xyerrorbars
-#plot 'tmp.timelag' using 1:2:3:4 with xyerrorbars
-set logscale xy
+set title "Yo"
+set xlabel "Temporal Frequency [1/days]"
 set xrange [0.005:0.603]
-set yrange [0.0005:100]
-# set yrange [:1]
-#    (x, y, ydelta),
-#    (x, y, ylow, yhigh),
-#    (x, y, xdelta),
-#    (x, y, xlow, xhigh),
-#    (x, y, xdelta, ydelta), or
-#    (x, y, xlow, xhigh, ylow, yhigh).
- pause -1
+
+# Use this if you want to show PSD directly on a log axis
+# set logscale xy
+# set ylabel "PSD"
+# set yrange [0.0005:100]
+# plot '%FILE1' using 1:(10**($2+2)):3:4 with xyerrorbars, \
+# '%FILE2' using 1:(10**($2+2)):3:4 with xyerrorbars
+
+# Use this if you want to show the log value of the PSD
+set logscale x
+set ylabel "log_{10} PSD"
+set yrange [-5.5:0.5]
+plot '%FILE1' using 1:2:3:4 with xyerrorbars, \
+'%FILE2' using 1:2:3:4 with xyerrorbars
\ No newline at end of file