2016-07-11 17:14:36 +00:00
|
|
|
|
#!/usr/local/bin/bash
|
|
|
|
|
|
|
2016-07-22 21:58:25 +00:00
|
|
|
|
timestep="0.1"
|
|
|
|
|
|
lc_dir="data/STORM_III/lightcurves/Δt＝${timestep}"
|
|
|
|
|
|
ref_band='HST(𝛌＝1367)'
|
|
|
|
|
|
ref_curve="$lc_dir/${ref_band}.lc"
|
|
|
|
|
|
error_type="2"
|
|
|
|
|
|
# error types:
|
|
|
|
|
|
# 0 for covariance matrix, 1 for likelihood function, 2 for monte carlo
|
|
|
|
|
|
|
|
|
|
|
|
echo Using lightcurves in $lc_dir.
|
|
|
|
|
|
|
|
|
|
|
|
for echo_curve in $lc_dir/*
|
2016-07-11 17:14:36 +00:00
|
|
|
|
do
|
2016-07-22 21:58:25 +00:00
|
|
|
|
echo_band=$(basename $echo_curve|sed 's@\(.*\)\.lc@\1@')
|
|
|
|
|
|
echo -n $(date "+%R")\: Running psdlag for $echo_band against $ref_band.
|
|
|
|
|
|
newfilename="analyses/${echo_band}≻${ref_band}_Δt＝${timestep}"
|
2016-07-14 20:12:56 +00:00
|
|
|
|
echo "2" > tmp.psdlagargs
|
2016-07-22 21:58:25 +00:00
|
|
|
|
echo "$ref_curve 0" >> tmp.psdlagargs
|
|
|
|
|
|
echo "$echo_curve 0" >> tmp.psdlagargs
|
2016-07-14 20:12:56 +00:00
|
|
|
|
echo "8 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603" >> tmp.psdlagargs
|
2016-07-20 03:12:16 +00:00
|
|
|
|
echo "0" >> tmp.psdlagargs
|
2016-07-14 20:12:56 +00:00
|
|
|
|
echo ".9999 .9999 .9999 0.49 0.42 0.30 0.36 1.2 -0.60 -0.89 -0.72 2.39 -0.74 -1.5 -1 -3.1 -1.3 -2.3 -2 -1.9 -1.6 -2.6 -2 -1.3 -1.8 -2.6 -1.2 0 -2 -2.7 -1.3 0.63" >> tmp.psdlagargs
|
|
|
|
|
|
echo "0:0 0" >> tmp.psdlagargs
|
2016-07-22 21:58:25 +00:00
|
|
|
|
echo $error_type >> tmp.psdlagargs
|
2016-07-14 20:12:56 +00:00
|
|
|
|
echo "0" >> tmp.psdlagargs
|
|
|
|
|
|
echo "100 50 50 mcmc.dat" >> tmp.psdlagargs
|
2016-07-22 21:58:25 +00:00
|
|
|
|
time bin/psdlag tmp.psdlagargs >> $newfilename
|
|
|
|
|
|
echo ""
|
2016-07-11 17:14:36 +00:00
|
|
|
|
done
|
