mirror of
https://asciireactor.com/otho/industry-website.git
synced 2024-11-22 20:55:07 +00:00
47 lines
1.7 KiB
Bash
47 lines
1.7 KiB
Bash
|
#!/usr/bin/env bash
|
|||
|
|
|||
|
# This script is meant to encapsulate all functions in proper order to produce
|
|||
|
# the analyses. Parameters that this script should control: Δt, σ type, input
|
|||
|
# dataset(lightcurve directory), more?
|
|||
|
|
|||
|
ref_band=$(cat ref_band)
|
|||
|
ref_curve="lightcurves/${ref_band}.lc"
|
|||
|
err_str=$(cat err_type)
|
|||
|
|
|||
|
mkdir -p analyses/tables
|
|||
|
|
|||
|
for echo_curve in lightcurves/*
|
|||
|
do
|
|||
|
# Grab and determine labels of analyses, skip if over the same band.
|
|||
|
echo_band=$(basename $echo_curve|sed 's|\(.*\)\.lc|\1|')
|
|||
|
if [[ $ref_band == $echo_band ]]; then continue; fi
|
|||
|
|
|||
|
echo "Analysing $echo_band ≺ $ref_band."
|
|||
|
|
|||
|
if [[ -e "analyses/${ref_band}_≺_${echo_band}/" ]]; then
|
|||
|
echo "Results already exists. Create tables from stored results."
|
|||
|
cp analyses/${ref_band}_≺_${echo_band}/*.out .
|
|||
|
else
|
|||
|
time python scripts/analyze_lightcurve.py $ref_curve $echo_curve >> ${echo_curve}.log
|
|||
|
mkdir -p "analyses/${ref_band}_≺_${echo_band}"
|
|||
|
cp *.out analyses/${ref_band}_≺_${echo_band}/
|
|||
|
fi
|
|||
|
|
|||
|
echo "Tabling PSD and time lags referred to ${ref_band} for $echo_band{${err_str}}."
|
|||
|
|
|||
|
# Propagate tables into analyses/tables
|
|||
|
echoPSD_tabfile=analyses/tables/PSD_${echo_band}_\{${err_str}\}.tab
|
|||
|
refPSD_tabfile=analyses/tables/PSD_${ref_band}_\{${err_str}\}.tab
|
|||
|
timelag_tabfile=analyses/tables/timelag_${ref_band}_≺_${echo_band}_\{${err_str}\}.tab
|
|||
|
|
|||
|
# Output curves to temporary files using perl script, move tables to
|
|||
|
# permanent location. This just assumes there are no conflicts.
|
|||
|
scripts/extract_tables.pl
|
|||
|
mv tmp.echoPSD $echoPSD_tabfile
|
|||
|
mv tmp.refPSD $refPSD_tabfile
|
|||
|
mv tmp.timelag $timelag_tabfile
|
|||
|
rm *.out
|
|||
|
done
|
|||
|
|
|||
|
rm tmp.*
|