Science/psdlag-agn
1
0
Fork 0
mirror of https://asciireactor.com/otho/psdlag-agn.git synced 2024-11-21 20:25:07 +00:00
Code Issues Packages Projects Releases Wiki Activity

working on automating plots

Browse Source
This commit is contained in:
othocaeS 2016-07-12 16:23:24 -04:00
parent b6dd5fc0a9
commit 6e473cfb09
2 changed files with 30 additions and 21 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

14
plots_playground/psd_tmp.gp
View File

@ -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),

37
scripts/templates/psd_freq.gp
View File

@ -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