aboutsummaryrefslogtreecommitdiffstats
diff options
context:
space:
mode:
-rwxr-xr-xscripts/calc_lx_simple.sh133
1 files changed, 0 insertions, 133 deletions
diff --git a/scripts/calc_lx_simple.sh b/scripts/calc_lx_simple.sh
deleted file mode 100755
index a5ae978..0000000
--- a/scripts/calc_lx_simple.sh
+++ /dev/null
@@ -1,133 +0,0 @@
-#!/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"
-