diff --git a/plots/psd_freq.gp b/plots/psd_freq.gp
new file mode 100644
index 0000000..b3972fa
--- /dev/null
+++ b/plots/psd_freq.gp
@@ -0,0 +1,20 @@
+
+
+
+# This method assumes data is in logarithm units already.
+
+plot '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, \
+'tables/PSD_lcdrive.dat' using 1:(10**($2+2)):3:4 with xyerrorbars
+#plot 'tmp.timelag' using 1:2:3:4 with xyerrorbars
+set logscale xy
+set xrange [0.005:0.603]
+set yrange [1:10000]
+# 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
\ No newline at end of file
diff --git a/plots/psd_tmp.gp b/plots/psd_tmp.gp
new file mode 100644
index 0000000..98b89c9
--- /dev/null
+++ b/plots/psd_tmp.gp
@@ -0,0 +1,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
+
+plot 'tmp.sourcePSD' using 1:($2+2):3:4 with xyerrorbars, \
+'tmp.reprocPSD' using 1:($2+2):3:4 with xyerrorbars, \
+'tmp.ccPSD' using 1:($2+2):3:4 with xyerrorbars
+#plot 'tmp.timelag' using 1:2:3:4 with xyerrorbars
+set logscale x
+set xrange [0.005:0.603]
+# 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
\ No newline at end of file
diff --git a/plots/timelag_freq.gp b/plots/timelag_freq.gp
new file mode 100644
index 0000000..3b93d6d
--- /dev/null
+++ b/plots/timelag_freq.gp
@@ -0,0 +1,8 @@
+
+
+
+plot 'tables/timelag_lc1.dat' using 1:2:3:4 with xyerrorbars, \
+'tables/timelag_lc2.dat' using 1:2:3:4 with xyerrorbars
+set logscale x
+set xrange [0.005:0.603]
+pause -1
\ No newline at end of file
diff --git a/plots/timelag_tmp.gp b/plots/timelag_tmp.gp
new file mode 100644
index 0000000..63edb2f
--- /dev/null
+++ b/plots/timelag_tmp.gp
@@ -0,0 +1,18 @@
+
+
+# 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
+
+plot 'tmp.timelag' using 1:2:3:4 with xyerrorbars
+#plot 'tmp.timelag' using 1:2:3:4 with xyerrorbars
+set logscale x
+set xrange [0.005:0.603]
+# 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
\ No newline at end of file