mirror of
https://asciireactor.com/otho/psdlag-agn.git
synced 2024-11-25 16:45:07 +00:00
35 lines
1015 B
Gnuplot
35 lines
1015 B
Gnuplot
# This method assumes y values are in logarithm units already.
|
|
|
|
set terminal pdf
|
|
set output "%OUTPUTFILE"
|
|
# set termoption dash
|
|
set termopt enhanced
|
|
|
|
set title "%TITLE\n{/*0.7 %SUBTITLE}" font "Times"
|
|
set xlabel "Temporal Frequency [days^{-1}]" font "Times"
|
|
set xrange [0.005:0.620]
|
|
set xtics font "Times"
|
|
set ytics font "Times"
|
|
|
|
# Use this if you want to show PSD directly on a log axis
|
|
# set logscale xy
|
|
# set ylabel "PSD" font "Times"
|
|
# set yrange [%YMIN:%YMAX]
|
|
# 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"
|
|
|
|
# Use this if you want to show the log value of the PSD
|
|
set logscale x
|
|
set ylabel "log_{10} PSD" font "Times"
|
|
#set yrange [%YMIN:%YMAX]
|
|
#set yrange [-5.5:0.5]
|
|
|
|
set key outside
|
|
set key center bottom
|
|
set key font "Times"
|
|
|
|
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"
|