1 files changed, 0 insertions, 113 deletions
diff --git a/bin/calc_coolfunc.sh b/bin/calc_coolfunc.sh
deleted file mode 100755
index d1d0b35..0000000
--- a/bin/calc_coolfunc.sh
+++ /dev/null
@@ -1,113 +0,0 @@
-#!/bin/sh
-##
-## Calculate the 'cooling function' profile with respect to the
-## given 'temperature profile' and the average abundance, redshift,
-## and column density nH, using the XSPEC model 'wabs*apec'.
-##
-## NOTE:
-## The output cooling function values should be the 'flux' values
-## with unit 'photon/s/cm^2' (different to 'calc_coolfunc_bands.sh').
-## These results will be used by 'fit_{beta,dbeta}_sbp' to derive the
-## (3D) gas density profile from (2D) surface brightness profile,
-## whose values have unit 'photon/cm^2/pixel/s'.
-##
-## Weitian LI
-## Created: 2012-08-17
-##
-## Change logs:
-## 2017-02-17, Weitian LI
-##   * Rename from 'coolfunc_calc.sh' to 'calc_coolfunc.sh'
-##   * Clean up that do not calculate and output <coolfunc_bolo>
-## 2017-02-16, Weitian LI
-##   * Do not calculate and output 'flux_cnt_ratio.txt'
-##
-
-## cmdline arguments {{{
-if [ $# -ne 5 ]; then
-    printf "usage:\n"
-    printf "    `basename $0` <tprofile> <avg_abund> <nH> <redshift> <coolfunc_outfile>\n"
-    exit 1
-fi
-TPROFILE=$1
-ABUNDANCE=$2
-N_H=$3
-REDSHIFT=$4
-COOLFUNC_DAT=$5
-NORM=`cosmo_calc ${REDSHIFT} | grep 'norm.*cooling_function' | awk -F':' '{ print $2 }'`
-
-if [ ! -r "${TPROFILE}" ]; then
-    printf "ERROR: given tprofile '${TPROFILE}' NOT accessiable\n"
-    exit 2
-fi
-[ -e "${COOLFUNC_DAT}" ] && rm -f ${COOLFUNC_DAT}
-## arguments }}}
-
-XSPEC_CF_XCM="_calc_coolfunc.xcm"
-[ -e "${XSPEC_CF_XCM}" ] && rm -f ${XSPEC_CF_XCM}
-
-## generate xspec script {{{
-cat > ${XSPEC_CF_XCM} << _EOF_
-## XSPEC Tcl script
-## Calculate the cooling function profile w.r.t the temperature profile.
-##
-## Generated by: `basename $0`
-## Date: `date`
-
-set xs_return_results 1
-set xs_echo_script 0
-# set tcl_precision 12
-## set basic data {{{
-set nh ${N_H}
-set redshift ${REDSHIFT}
-set abundance ${ABUNDANCE}
-set norm ${NORM}
-## basic }}}
-
-## xspec related {{{
-# debug settings {{{
-chatter 0
-# debug }}}
-query yes
-abund grsa
-dummyrsp 0.01 100.0 4096 linear
-# load model 'wabs*apec' to calc cooling function
-model wabs*apec & \${nh} & 1.0 & \${abundance} & \${redshift} & \${norm} & /*
-## xspec }}}
-
-## set input and output filename & open files
-set tpro_fn "${TPROFILE}"
-set cf_fn "${COOLFUNC_DAT}"
-if { [ file exists \${cf_fn} ] } {
-    exec rm -fv \${cf_fn}
-}
-
-## open files
-set tpro_fd [ open \${tpro_fn} r ]
-set cf_fd [ open \${cf_fn} w ]
-
-## read data from tprofile line by line
-while { [ gets \${tpro_fd} tpro_line ] != -1 } {
-    scan \${tpro_line} "%f %f" radius temperature
-    #puts "radius: \${radius}, temperature: \${temperature}"
-    # set temperature value
-    newpar 2 \${temperature}
-    # calc flux & tclout
-    flux 0.7 7.0
-    tclout flux 1
-    scan \${xspec_tclout} "%f %f %f %f" _ _ _ cf_photon
-    #puts "cf: \${cf_photon}"
-    puts \${cf_fd} "\${radius}    \${cf_photon}"
-}
-
-## close opened files
-close \${tpro_fd}
-close \${cf_fd}
-
-## exit
-tclexit
-_EOF_
-## xcm generation }}}
-
-## invoke xspec to calc
-printf "invoking XSPEC to calculate cooling function profile ...\n"
-xspec - ${XSPEC_CF_XCM} > /dev/null