1 files changed, 43 insertions, 0 deletions
diff --git a/mass_profile/fit_sbp.sh b/mass_profile/fit_sbp.sh
new file mode 100755
index 0000000..76aaaaf
--- /dev/null
+++ b/mass_profile/fit_sbp.sh
@@ -0,0 +1,43 @@
+#!/bin/sh
+#
+
+if [ $# -ne 1 ]; then
+    printf "usage: $0 <mass_conf>\n"
+    exit 1
+fi
+cfg_file=$1
+if [ "$0" = `basename $0` ]; then
+    script_path=`which $0`
+    base_path=`dirname ${script_path}`
+else
+    base_path=`dirname $0`
+fi
+
+sbp_cfg=`grep '^sbp_cfg' $cfg_file | awk '{ print $2 }'`
+t_data_file=`grep '^t_data_file' $cfg_file | awk '{ print $2 }'`
+t_param_file=`grep '^t_param_file' $cfg_file | awk '{ print $2 }'`
+nh=`grep '^nh' $cfg_file | awk '{ print $2 }'`
+abund=`grep '^abund' $cfg_file | awk '{ print $2 }'`
+z=`grep '^z' $sbp_cfg | awk '{ print $2 }'`
+cm_per_pixel=`${base_path}/calc_distance ${z} | grep 'cm_per_pixel' | awk '{ print $2 }'`
+sed -i'' "s/^cm_per_pixel.*$/cm_per_pixel    ${cm_per_pixel}/" ${sbp_cfg}
+cfunc_file=`grep '^cfunc_file' $sbp_cfg | awk '{ print $2 }'`
+T_file=`grep '^T_file' $sbp_cfg | awk '{ print $2 }'`
+
+if grep -q '^beta2' $sbp_cfg; then
+    MODEL="double-beta"
+    PROG=fit_dbeta_sbp
+else
+    MODEL="single-beta"
+    PROG=fit_beta_sbp
+fi
+
+$base_path/fit_wang2012_model $t_data_file $t_param_file $cm_per_pixel 2> /dev/null
+cp wang2012_dump.qdp $T_file
+if [ ! -f ${cfunc_file} ]; then
+    $base_path/coolfunc_calc2.sh $T_file $abund $nh $z $cfunc_file
+fi
+$base_path/$PROG $sbp_cfg
+printf "## MODEL: ${MODEL}\n"
+printf "## z: ${z}\n"
+