aboutsummaryrefslogtreecommitdiffstats
path: root/mass_profile/get_center_params.sh
diff options
context:
space:
mode:
Diffstat (limited to 'mass_profile/get_center_params.sh')
-rwxr-xr-xmass_profile/get_center_params.sh46
1 files changed, 46 insertions, 0 deletions
diff --git a/mass_profile/get_center_params.sh b/mass_profile/get_center_params.sh
new file mode 100755
index 0000000..a250e07
--- /dev/null
+++ b/mass_profile/get_center_params.sh
@@ -0,0 +1,46 @@
+#!/bin/bash
+origin_point=$PWD
+ZERO=`readlink -f $0`
+basedir=`dirname $ZERO`
+
+for i in `cat list.txt`
+do
+# echo $i
+ file_path=`dirname $i`
+ cd $file_path
+ mkdir -p center_results
+ cp global.cfg center_results
+ t_data_file=`grep ^t_data_file global.cfg|awk '{print $2}'`
+ t_param_file=`grep ^t_param_file global.cfg|awk '{print $2}'`
+ sbp_cfg=`grep ^sbp_cfg global.cfg|awk '{print $2}'`
+ cm_per_pixel=`grep '^cm_per_pixel' $sbp_cfg|awk '{print $2}'`
+ radius_sbp_file=`grep '^radius_sbp_file' global.cfg|awk '{print $2}'`
+ radius_file=`grep '^radius_file' $sbp_cfg|awk '{print $2}'`
+ sbp_file=`grep '^sbp_file' $sbp_cfg|awk '{print $2}'`
+ cfunc_file=`grep '^cfunc_file' $sbp_cfg|awk '{print $2}'`
+ T_file=`grep '^T_file' $sbp_cfg |awk '{print $2}'`
+
+ cp $t_data_file $t_param_file $sbp_cfg $radius_sbp_file $radius_file $sbp_file $cfunc_file center_results
+ cd center_results
+ nh=`grep '^nh' global.cfg |awk '{print $2}'`
+ z=`grep '^z' ${sbp_cfg}|awk '{print $2}'`
+ abund=`grep '^abund' global.cfg |awk '{print $2}'`
+
+ $basedir/fit_wang2012_model $t_data_file $t_param_file $cm_per_pixel
+ #echo $t_param_file
+ mv wang2012_dump.qdp wang2012_center_dump.qdp
+ cp -f wang2012_center_dump.qdp t_profile_dump.qdp
+ $basedir/coolfunc_calc.sh wang2012_center_dump.qdp $abund $nh $z cfunc.dat
+ if grep '^beta2' ${sbp_cfg}
+ then
+ $basedir/fit_dbeta_sbp $sbp_cfg
+ mv dbeta_param.txt sbp_param_center.txt
+ else
+ $basedir/fit_beta_sbp $sbp_cfg
+ mv beta_param.txt sbp_param_center.txt
+ fi
+
+ $basedir/fit_nfw_mass mass_int.dat $z >nfw_param_center.txt
+
+ cd $origin_point
+done