2016-07-12 20:23:24 +00:00
|
|
|
# This method assumes y values are in logarithm units already.
|
2016-07-12 15:19:13 +00:00
|
|
|
|
2016-08-10 02:38:11 +00:00
|
|
|
set terminal pdf
|
2016-07-13 03:28:11 +00:00
|
|
|
set output "%OUTPUTFILE"
|
2016-08-10 02:38:11 +00:00
|
|
|
# set termoption dash
|
2016-08-04 18:07:46 +00:00
|
|
|
set termopt enhanced
|
2016-07-12 15:19:13 +00:00
|
|
|
|
2016-08-10 02:38:11 +00:00
|
|
|
set title "%TITLE\n{/*0.7 %SUBTITLE}" font "Times"
|
|
|
|
|
set xlabel "Temporal Frequency [days^{-1}]" font "Times"
|
2016-07-27 00:45:33 +00:00
|
|
|
set xrange [0.005:0.620]
|
2016-08-10 02:38:11 +00:00
|
|
|
set xtics font "Times"
|
|
|
|
|
set ytics font "Times"
|
2016-07-12 15:19:13 +00:00
|
|
|
|
2016-07-12 20:23:24 +00:00
|
|
|
# Use this if you want to show PSD directly on a log axis
|
|
|
|
|
# set logscale xy
|
2016-08-10 02:38:11 +00:00
|
|
|
# set ylabel "PSD" font "Times"
|
2016-07-13 20:16:09 +00:00
|
|
|
# set yrange [%YMIN:%YMAX]
|
2016-07-13 15:09:20 +00:00
|
|
|
# plot '%FILE1' using 1:(10**($2+2)):3:4 with xyerrorbars title "%LABEL1", \
|
|
|
|
|
# '%FILE2' using 1:(10**($2+2)):3:4 with xyerrorbars title "%LABEL2"
|
2016-07-12 15:19:13 +00:00
|
|
|
|
2016-07-12 20:23:24 +00:00
|
|
|
# Use this if you want to show the log value of the PSD
|
|
|
|
|
set logscale x
|
2016-08-10 02:38:11 +00:00
|
|
|
set ylabel "log_{10} PSD" font "Times"
|
2016-07-13 20:16:09 +00:00
|
|
|
#set yrange [%YMIN:%YMAX]
|
2016-07-27 06:11:17 +00:00
|
|
|
#set yrange [-5.5:0.5]
|
|
|
|
|
|
|
|
|
|
set key outside
|
|
|
|
|
set key center bottom
|
2016-08-10 02:38:11 +00:00
|
|
|
set key font "Times"
|
2016-07-27 06:11:17 +00:00
|
|
|
|
2016-08-10 02:38:11 +00:00
|
|
|
set yzeroaxis
|
|
|
|
|
|
|
|
|
|
plot '%FILE1' using 1:2:3:4 with xyerrorbars title "%LABEL1" pt 0 ps 1 lt rgb "green", \
|
|
|
|
|
'%FILE2' using 1:2:3:4 with xyerrorbars title "%LABEL2" pt 0 ps 1 lt rgb "red"
|
