psdlag-agn/scripts/plot_psdlag.sh

#!/usr/bin/env bash

ywindow_of () {
    max=0
    min=-5
    while read line
#    for line in $(cat $1)
    do
        read xval yval xerr yerr <<< $line
        echo $xval $yval
        if (bc <<< "$yval > $max"); then max=$yval; fi
        if (bc <<< "$yval < $min"); then min=$yval; fi
#    done
    done < $1
    echo $(bc <<< "$min - 0.5") $(bc <<< "$max + 0.5")
}

mkdir -p analyses/tables
mkdir -p analyses/plots


for analysis in analyses/*
do
    # Grab and determine labels of analyses, skip if over the same band.
    ref_band=$(basename $analysis|sed 's@\([^≺]*\)[_ ]≺[_ ][^≺_ ]*[_ ][^_ ]*[_ ][^_ ]*@\1@')
    echo_band=$(basename $analysis|sed 's@[^≺]*[_ ]≺[_ ]\([^≺_ ]*\)[_ ][^_ ]*[_ ][^_ ]*@\1@')
    if [[ $ref_band == $echo_band ]]; then continue; fi
    err_type=$(basename $analysis|sed 's@[^≺]*[_ ]≺[_ ][^≺_ ]*[_ ][^_ ]*[_ ]\(σ∊[CLM][MFC]\)@\1@')

    # Prepare files
    echo "Plotting PSD and time lags for $echo_band, referred to ${ref_band}."
    echoPSD_tabfile=analyses/tables/${echo_band}_${err_type}_PSD.tab
    refPSD_tabfile=analyses/tables/${ref_band}_${err_type}_PSD.tab
    timelag_tabfile=analyses/tables/${ref_band}_≺_${echo_band}_${err_type}_timelag.tab
    PSD_plotfile=analyses/plots/${ref_band}_≺_${echo_band}_${err_type}_PSD.png
    timelag_plotfile=analyses/plots/${ref_band}_≺_${echo_band}_${err_type}_timelag.png

    # Output curves to temporary files using perl script, move tables to
    # permanent location. This just assumes there are no conflicts.
    scripts/create_tables.pl $analysis > /dev/null
    mv tmp.echoPSD $echoPSD_tabfile
    mv tmp.refPSD $refPSD_tabfile
    mv tmp.timelag $timelag_tabfile


    # Plot PSD and save using gnuplot
#    read ymin ymax <<< $(ywindow_of $timelag_tabfile)
#    echo $ymin $ymax
    cat scripts/templates/psd_freq.gp|
        sed "s|%TITLE|Power Spectrum for Lightcurves $echo_band and $ref_band|"|
        sed "s|%SUBTITLE|as reported by Fausnaugh et. al, STORM III, 2016|"|
        sed "s|%FILE1|$refPSD_tabfile|"|
        sed "s|%LABEL1|${ref_band} PSD|"|
        sed "s|%FILE2|$echoPSD_tabfile|"|
        sed "s|%LABEL2|${echo_band} PSD|"|
        sed "s|%YMIN|$ymin|"|sed "s|%YMAX|$ymax|"|
        sed "s|%OUTPUTFILE|$PSD_plotfile|" > tmp.gp
    gnuplot tmp.gp

    # Plot time lags and save using gnuplot
#    read ymin ymax <<< $(ywindow_of $timelag_tabfile)
    cat scripts/templates/timelag_freq.gp|
        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|%FILE1|$timelag_tabfile|"|
        sed "s|%LABEL1|${echo_band} Lag|"|
        sed "s|%YMIN|$ymin|"|sed "s|%YMAX|$ymax|"|
        sed "s|%OUTPUTFILE|$timelag_plotfile|" > tmp.gp
    gnuplot tmp.gp
done

rm tmp.*