aboutsummaryrefslogtreecommitdiffstats
path: root/mass_profile
diff options
context:
space:
mode:
authorAaron LI <aaronly.me@outlook.com>2016-06-08 21:13:02 +0800
committerAaron LI <aaronly.me@outlook.com>2016-06-08 21:13:02 +0800
commit1bb88e9d112c21f65a0bf74d57cc19397f9ed29e (patch)
tree9b6528043d68ce94d31e8b001679b225e857e4fd /mass_profile
parent03abfd529244c1bd7d18baaae6ad1781ce51324d (diff)
downloadchandra-acis-analysis-1bb88e9d112c21f65a0bf74d57cc19397f9ed29e.tar.bz2
Remove unused coolfunc_calc_bolo.sh
Diffstat (limited to 'mass_profile')
-rwxr-xr-xmass_profile/coolfunc_calc_bolo.sh139
1 files changed, 0 insertions, 139 deletions
diff --git a/mass_profile/coolfunc_calc_bolo.sh b/mass_profile/coolfunc_calc_bolo.sh
deleted file mode 100755
index 5c4ef5d..0000000
--- a/mass_profile/coolfunc_calc_bolo.sh
+++ /dev/null
@@ -1,139 +0,0 @@
-#!/bin/sh
-#
-# unalias -a
-#
-###########################################################
-## Task: ##
-## Calc `cooling function' data according to ##
-## given `temperature profile' ##
-## ##
-## NOTE: ##
-## given `tprofile': <radius> <temperature> ##
-## calc `cooling function' by invoking `XSPEC' ##
-## using model `wabs*apec' ##
-## ##
-## LIweitiaNux <liweitianux@gmail.com> ##
-## August 17, 2012 ##
-###########################################################
-
-## cmdline arguments {{{
-if [ $# -ne 5 ]; then
- printf "usage:\n"
- printf " `basename $0` <tprofile> <avg_abund> <nH> <redshift> <coolfunc_outfile>\n"
- exit 1
-fi
-base_path=`dirname $0`
-TPROFILE=$1
-ABUND_VAL=$2
-N_H=$3
-REDSHIFT=$4
-NORM=`$base_path/calc_distance $REDSHIFT|grep norm|awk '{print $2}'`
-
-echo $NORM
-
-
-COOLFUNC_DAT=$5
-#COOLFUNC_DAT_RATIO=flux_cnt_ratio.txt
-
-if [ ! -r "${TPROFILE}" ]; then
- printf "ERROR: given tprofile '${TPROFILE}' NOT accessiable\n"
- exit 2
-fi
-[ -e "${COOLFUNC_DAT}" ] && rm -f ${COOLFUNC_DAT}
-#[ -e "${COOLFUNC_DAT_RATIO}" ] && rm -f ${COOLFUNC_DAT_RATIO}
-## arguments }}}
-
-## specify variable name outside while loop
-## otherwise the inside vars invisible
-XSPEC_CF_XCM="_coolfunc_calc.xcm"
-[ -e "${XSPEC_CF_XCM}" ] && rm -f ${XSPEC_CF_XCM}
-
-## generate xspec script {{{
-cat >> ${XSPEC_CF_XCM} << _EOF_
-## XSPEC Tcl script
-## calc cooling function data
-##
-## generated by: `basename $0`
-## date: `date`
-
-set xs_return_results 1
-set xs_echo_script 0
-# set tcl_precision 12
-dummyrsp .01 100 4096
-## set basic data {{{
-set nh ${N_H}
-set redshift ${REDSHIFT}
-set abund_val ${ABUND_VAL}
-set norm ${NORM}
-## basic }}}
-
-## xspec related {{{
-# debug settings {{{
-chatter 0
-# debug }}}
-query yes
-abund grsa
-dummyrsp 0.3 11.0 1024
-# load model 'wabs*apec' to calc cooling function
-#model wabs*apec & \${nh} & 1.0 & \${abund_val} & \${redshift} & \${norm} &
-model wabs*apec & 0 & 1.0 & \${abund_val} & \${redshift} & \${norm} &
-## xspec }}}
-
-## set input and output filename
-set tpro_fn "${TPROFILE}"
-set cf_fn "${COOLFUNC_DAT}"
-##set cff_fn "${COOLFUNC_DAT_RATIO}"
-if { [ file exists \${cf_fn} ] } {
- exec rm -fv \${cf_fn}
-}
-
-#if { [ file exists \${cff_fn} ] } {
-# exec rm -fv \${cff_fn}
-#}
-
-## open files
-set tpro_fd [ open \${tpro_fn} r ]
-set cf_fd [ open \${cf_fn} w ]
-#set cff_fd [ open \${cff_fn} w ]
-
-## read data from tprofile line by line
-while { [ gets \${tpro_fd} tpro_line ] != -1 } {
- # gets one line
- scan \${tpro_line} "%f %f" radius temp_val
- #puts "radius: \${radius}, temperature: \${temp_val}"
- # set temperature value
- newpar 2 \${temp_val}
- # calc flux & tclout
- flux 0.01 100.0
- tclout flux 1
- scan \${xspec_tclout} "%f %f %f %f" cf_data holder holder holder
- #puts "cf_data: \${cf_data}"
- puts \${cf_fd} "\${radius} \${cf_data}"
-# flux 0.01 100.0
-# tclout flux 1
-# scan \${xspec_tclout} "%f %f %f %f" cff_data holder holder holder
-# puts \${cff_fd} "\${radius} [expr \${cff_data}/\${cf_data}]"
-}
-
-## close opened files
-close \${tpro_fd}
-close \${cf_fd}
-
-## exit
-tclexit
-_EOF_
-
-## extract xcm }}}
-
-## invoke xspec to calc
-printf "invoking XSPEC to calculate cooling function data ...\n"
-# xspec - ${XSPEC_CF_XCM}
-xspec - ${XSPEC_CF_XCM} > /dev/null
-
-## clean
-# if [ -e "${XSPEC_CF_XCM}" ]; then
-# rm -f ${XSPEC_CF_XCM}
-# fi
-
-exit 0
-