blob: a5ae9787a04a2243153b44fcb627aa3b85b14e1f (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
|
#!/bin/sh
##
## Run calc_lx in BATCH mode
##
## Weitian LI <liweitianux@gmail.com>
## Created: 2012/08/31
##
UPDATED="2014/06/18"
##
## Changelogs:
## 2015/06/03:
## * Replaced 'grep' with '\grep', 'ls' with '\ls'
## 2014/06/18:
## * use env variable 'MASS_PROFILE_DIR'
##
unalias -a
export LC_COLLATE=C
# fix path for python
export PATH="/usr/bin:/usr/local/bin:$PATH"
## usage, `path_conffile' is the configuration file
## which contains the `path' to each `repro/mass' directory
if [ $# -ne 1 ]; then
printf "usage:\n"
printf " `basename $0` <mass_dir>\n"
printf "\nNOTE:\n"
printf " script cannot handle \`~' in path\n"
printf "\nVersion: ${UPDATED}\n"
exit 1
fi
## set the path to the script {{{
if [ -n "${MASS_PROFILE_DIR}" ]; then
CALCLX_SCRIPT="${MASS_PROFILE_DIR}/calc_lx"
CFCALC_SCRIPT="${MASS_PROFILE_DIR}/coolfunc_calc.sh"
FIT_TPRO="${MASS_PROFILE_DIR}/fit_wang2012_model"
else
printf "ERROR: environment variable 'MASS_PROFILE_DIR' not set.\n"
exit 2
fi
if [ -z "${CALCLX_SCRIPT}" ]; then
printf "ERROR: \`CALCLX_SCRIPT' not set\n"
exit 250
elif [ ! -r ${CALCLX_SCRIPT} ]; then
printf "ERROR: CANNOT access script \`${CALCLX_SCRIPT}'\n"
exit 251
fi
## script path }}}
# result lines
RES_LINE=100
# process dir
MDIR="$1"
# mass fitting conf
MCONF="fitting_mass.conf"
# process
cd ${MDIR}
printf "Entered dir: \``pwd`'\n"
# conf file
if [ ! -r "${MCONF}" ]; then
printf "ERROR: configuration file \`${MCONF}' not accessiable\n"
else
LOGFILE="calclx_`date '+%Y%m%d%H'`.log"
[ -e "${LOGFILE}" ] && mv -fv ${LOGFILE} ${LOGFILE}_bak
TOLOG="tee -a ${LOGFILE}"
# fitting_mass logfile, get R500 from it
MLOG=`\ls ${MCONF%.[confgCONFG]*}*.log | tail -n 1`
R500_VAL=`tail -n ${RES_LINE} ${MLOG} | \grep '^r500' | awk '{ print $2 }'`
R200_VAL=`tail -n ${RES_LINE} ${MLOG} | \grep '^r200' | awk '{ print $2 }'`
# radius_sbp_file {{{
RSBP=`\grep '^radius_sbp_file' ${MCONF} | awk '{ print $2 }'`
TMP_RSBP="_tmp_rsbp.txt"
[ -e "${TMP_RSBP}" ] && rm -f ${TMP_RSBP}
cat ${RSBP} | sed 's/#.*$//' | \grep -Ev '^\s*$' > ${TMP_RSBP}
RSBP="${TMP_RSBP}"
# rsbp }}}
TPRO_TYPE=`\grep '^t_profile' ${MCONF} | awk '{ print $2 }'`
TPRO_DATA=`\grep '^t_data_file' ${MCONF} | awk '{ print $2 }'`
TPRO_PARA=`\grep '^t_param_file' ${MCONF} | awk '{ print $2 }'`
SBP_CONF=`\grep '^sbp_cfg' ${MCONF} | awk '{ print $2 }'`
ABUND=`\grep '^abund' ${MCONF} | awk '{ print $2 }'`
NH=`\grep '^nh' ${MCONF} | awk '{ print $2 }'`
Z=`\grep '^z' ${SBP_CONF} | awk '{ print $2 }'`
cm_per_pixel=`\grep '^cm_per_pixel' ${SBP_CONF} | awk '{ print $2 }'`
CF_FILE=`\grep '^cfunc_file' ${SBP_CONF} | awk '{ print $2 }'`
printf "## use logfile: \`${LOGFILE}'\n"
printf "## working directory: \``pwd -P`'\n" | ${TOLOG}
printf "## use configuration files: \`${MCONF}, ${SBP_CONF}'\n" | ${TOLOG}
printf "## use radius_sbp_file: \`${RSBP}'\n" | ${TOLOG}
printf "## R500 (kpc): \`${R500_VAL}'\n" | ${TOLOG}
printf "## R200 (kpc): \`${R200_VAL}'\n" | ${TOLOG}
printf "## redshift: \`${Z}'\n" | ${TOLOG}
printf "## abund: \`${ABUND}'\n" | ${TOLOG}
printf "## nh: \`${NH}'\n" | ${TOLOG}
printf "## T_profile type: \`${TPRO_TYPE}'\n" | ${TOLOG}
printf "## cfunc_file: \`${CF_FILE}'\n" | ${TOLOG}
## fit temperature profile {{{
T_FILE="_tpro_dump.qdp"
if [ "${TPRO_TYPE}" = "wang2012" ]; then
printf "fitting temperature profile (wang2012) ...\n"
[ -e "wang2012_dump.qdp" ] && mv -fv wang2012_dump.qdp wang2012_dump.qdp_bak
[ -e "fit_result.qdp" ] && mv -fv fit_result.qdp fit_result.qdp_bak
${FIT_TPRO} ${TPRO_DATA} ${TPRO_PARA} ${cm_per_pixel} 2> /dev/null
mv -fv wang2012_dump.qdp ${T_FILE}
[ -e "wang2012_dump.qdp_bak" ] && mv -fv wang2012_dump.qdp_bak wang2012_dump.qdp
[ -e "fit_result.qdp_bak" ] && mv -fv fit_result.qdp_bak fit_result.qdp
else
printf "ERROR: invalid tprofile type: \`${TPRO_TYPE}'\n"
exit 10
fi
## tprofile }}}
## calc `flux_ratio' {{{
printf "calc flux_ratio ...\n"
CF_FILE="_cf_data.txt"
FLUX_RATIO="__flux_cnt_ratio.txt"
[ -e "flux_cnt_ratio.txt" ] && mv -fv flux_cnt_ratio.txt flux_cnt_ratio.txt_bak
printf "## CMD: sh ${CFCALC_SCRIPT} ${T_FILE} ${ABUND} ${NH} ${Z} ${CF_FILE}\n" | ${TOLOG}
sh ${CFCALC_SCRIPT} ${T_FILE} ${ABUND} ${NH} ${Z} ${CF_FILE}
mv -fv flux_cnt_ratio.txt ${FLUX_RATIO}
[ -e "flux_cnt_ratio.txt_bak" ] && mv -fv flux_cnt_ratio.txt_bak flux_cnt_ratio.txt
## flux_ratio }}}
printf "## CMD: ${CALCLX_SCRIPT} ${RSBP} ${FLUX_RATIO} ${Z} ${R500_VAL} ${TPRO_DATA}\n" | ${TOLOG}
printf "## CMD: ${CALCLX_SCRIPT} ${RSBP} ${FLUX_RATIO} ${Z} ${R200_VAL} ${TPRO_DATA}\n" | ${TOLOG}
L500=`${CALCLX_SCRIPT} ${RSBP} ${FLUX_RATIO} ${Z} ${R500_VAL} ${TPRO_DATA} | \grep '^Lx' | awk '{ print $2,$3,$4 }'`
L200=`${CALCLX_SCRIPT} ${RSBP} ${FLUX_RATIO} ${Z} ${R200_VAL} ${TPRO_DATA} | \grep '^Lx' | awk '{ print $2,$3,$4 }'`
printf "L500= ${L500} erg/s\n" | ${TOLOG}
printf "L200= ${L200} erg/s\n" | ${TOLOG}
fi
printf "\n++++++++++++++++++++++++++++++++++++++\n"
|
