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

It appears the plots are being created more-or-less correctly.

Browse Source
This commit is contained in:
othocaeS 2016-07-13 11:09:20 -04:00
parent 7f442b9f2b
commit 4a389ccbb5
3 changed files with 20 additions and 10 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
scripts/create_plots.sh
View File

@ -21,17 +21,19 @@ do
# Output curves to temporary files using perl script, move tables to # Output curves to temporary files using perl script, move tables to
# permanent location. This just assumes there are no conflicts. # permanent location. This just assumes there are no conflicts.
scripts/create_tables.pl $analysis scripts/create_tables.pl $analysis > /dev/null
mv tmp.echoPSD $echoPSD_tabfile mv tmp.echoPSD $echoPSD_tabfile
mv tmp.refPSD $refPSD_tabfile mv tmp.refPSD $refPSD_tabfile
mv tmp.timelag $timelag_tabfile mv tmp.timelag $timelag_tabfile
# Plot PSD and save using gnuplot # Plot PSD and save using gnuplot
cat scripts/templates/psd_freq.gp| cat scripts/templates/psd_freq.gp|
sed "s@%TITLE@Power Spectrum for Lightcurves $echo_band & $ref_band@"| sed "s@%TITLE@Power Spectrum for Lightcurves $echo_band \& $ref_band@"|
sed "s@%SUBTITLE@as reported by Fausnaugh et. al, STORM III, 2016@"| sed "s@%SUBTITLE@as reported by Fausnaugh et. al, STORM III, 2016@"|
sed "s@%FILE1@$refPSD_tabfile@"| sed "s@%FILE1@$refPSD_tabfile@"|
sed "s@%LABEL1@${ref_band} PSD@"|
sed "s@%FILE2@$echoPSD_tabfile@"| sed "s@%FILE2@$echoPSD_tabfile@"|
sed "s@%LABEL2@${echo_band} PSD@"|
sed "s@%OUTPUTFILE@$PSD_plotfile@" > tmp.gp sed "s@%OUTPUTFILE@$PSD_plotfile@" > tmp.gp
gnuplot tmp.gp gnuplot tmp.gp
@ -40,6 +42,7 @@ do
sed "s@%TITLE@Time Lag for Lightcurve $echo_band relative to $ref_band@"| sed "s@%TITLE@Time Lag for Lightcurve $echo_band relative to $ref_band@"|
sed "s@%SUBTITLE@as reported by Fausnaugh et. al, STORM III, 2016@"| sed "s@%SUBTITLE@as reported by Fausnaugh et. al, STORM III, 2016@"|
sed "s@%FILE1@$timelag_tabfile@"| sed "s@%FILE1@$timelag_tabfile@"|
sed "s@%LABEL1@${echo_band} Lag@"|
sed "s@%OUTPUTFILE@$timelag_plotfile@" > tmp.gp sed "s@%OUTPUTFILE@$timelag_plotfile@" > tmp.gp
gnuplot tmp.gp gnuplot tmp.gp
done done

10
scripts/templates/psd_freq.gp
View File

@ -3,7 +3,7 @@
set terminal png set terminal png
set output "%OUTPUTFILE" set output "%OUTPUTFILE"
set title "%TITLE\n{/*0.5 %SUBTITLE}" set title "%TITLE\n{/*0.7 %SUBTITLE}"
set xlabel "Temporal Frequency [days^{-1}]" set xlabel "Temporal Frequency [days^{-1}]"
set xrange [0.005:0.603] set xrange [0.005:0.603]
@ -11,12 +11,12 @@ set xrange [0.005:0.603]
# set logscale xy # set logscale xy
# set ylabel "PSD" # set ylabel "PSD"
# set yrange [0.0005:100] # set yrange [0.0005:100]
# plot '%FILE1' using 1:(10**($2+2)):3:4 with xyerrorbars, \ # plot '%FILE1' using 1:(10**($2+2)):3:4 with xyerrorbars title "%LABEL1", \
# '%FILE2' using 1:(10**($2+2)):3:4 with xyerrorbars # '%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 # Use this if you want to show the log value of the PSD
set logscale x set logscale x
set ylabel "log_{10} PSD" set ylabel "log_{10} PSD"
set yrange [-5.5:0.5] set yrange [-5.5:0.5]
plot '%FILE1' using 1:2:3:4 with xyerrorbars, \ plot '%FILE1' using 1:2:3:4 with xyerrorbars title "%LABEL1", \
'%FILE2' using 1:2:3:4 with xyerrorbars '%FILE2' using 1:2:3:4 with xyerrorbars title "%LABEL2"

13
scripts/templates/timelag_freq.gp
View File

@ -2,8 +2,15 @@ set terminal png
set termoption dash set termoption dash
set output "%OUTPUTFILE" set output "%OUTPUTFILE"
set title "%TITLE\n{/*0.7 %SUBTITLE}"
plot '%FILE1' using 1:2:3:4 with xyerrorbars\ set xlabel "Temporal Frequency [days^{-1}]"
set logscale x
set xrange [0.005:0.603] set xrange [0.005:0.603]
set logscale x
set ylabel "Lag [days]"
# Draw line at origin
set arrow from 0.005,0 to 0.603,0 nohead lt 3 lc rgb 'black' set arrow from 0.005,0 to 0.603,0 nohead lt 3 lc rgb 'black'
plot '%FILE1' using 1:2:3:4 with xyerrorbars title "%LABEL1"