Drop calc_distance support in favor of cosmo_calc

1 files changed, 26 insertions, 47 deletions

diff --git a/scripts/ciao_r500avgt.sh b/scripts/ciao_r500avgt.sh
index feed110..7568cdb 100755
--- a/scripts/ciao_r500avgt.sh
+++ b/scripts/ciao_r500avgt.sh
@@ -1,50 +1,45 @@
 #!/bin/sh
-#
-unalias -a
-export LC_COLLATE=C
-###########################################################
-## script to extract spectrum and prepare needed files   ##
-## for calculating the `average temperature'             ##
-## within (0.1-0.5 r500) region                          ##
-##                                                       ##
-## NOTE:                                                 ##
-## 1) r500 default in unit `pixel', if in `kpc',         ##
-## then `redshift z' and `calculator' are needed.        ##
-## 2) `background process' is same as `deproj_spectra'   ##
-## which supports `spectrum', `local' and `blanksky'     ##
-## 3) ARF/RMF files either provided or use the ARF/RMF   ##
-## of the outmost region                                 ##
-##                                                       ##
-## Weitian LI <liweitianux@gmail.com>                    ##
-## 2012/08/22                                            ##
-###########################################################
 ##
-## ChangeLogs
-## v1.1, 2012/08/26, LIweitiaNux
+## To extract the spectrum and prepare the necessary files for
+## calculating the average temperature within (0.1-0.5 r500) region.
+##
+## NOTE:
+## 1) r500 default in unit `pixel', if in `kpc',
+##    then `redshift z' and `calculator' are needed.
+## 2) `background process' is same as `deproj_spectra'
+##    which supports `spectrum', `local' and `blanksky'
+## 3) ARF/RMF files either provided or use the ARF/RMF
+##    of the outmost region
+##
+## Weitian LI <liweitianux@gmail.com>
+## Updated: 2016-06-08
+##
+## Change logs
+## v1.1, 2012/08/26, Weitian LI
 ##   modify `KPC_PER_PIX', able to use the newest version `calc_distance'
-## v1.2, 2012/08/26, LIweitiaNux
+## v1.2, 2012/08/26, Weitian LI
 ##   fix a bug with `DFT_BKGD'
-## v2.0, 2012/09/04, LIweitiaNux
+## v2.0, 2012/09/04, Weitian LI
 ##   add parameter `inner' and `outer' to adjust the region range
 ##   modify parameter `r500' to take `kpc' as the default unit
-## v2.1, 2012/10/05, LIweitiaNux
+## v2.1, 2012/10/05, Weitian LI
 ##   change `DFT_GRP_CMD' to `group 1 128 4 ...'
-## v3.0, 2013/02/09, LIweitiaNux
+## v3.0, 2013/02/09, Weitian LI
 ##   modify for new process
-## v3.1, 2015/05/27, Aaron LI
+## v3.1, 2015/05/27, Weitian LI
 ##   update 'DFT_ARF' & 'DFT_RMF' to find '*.arf' & '*.rmf' files
 ##   (specextract only use .arf & .rmf extensions since revision 2014-12)
-## v3.2, 2015/05/30, Aaron LI
+## v3.2, 2015/05/30, Weitian LI
 ##   Added options '-cmap he -bin factor 4' to ds9 command
-## v4.0, 2015/06/03, Aaron LI
+## v4.0, 2015/06/03, Weitian LI
 ##   * Copy needed pfiles to current working directory, and
 ##     set environment variable $PFILES to use these first.
 ##   * Replace 'grep' with '\grep', 'ls' with '\ls'
 ##   * replaced 'grppha' with 'dmgroup' to group spectra
 ##     (dmgroup will add history to fits file, while grppha NOT)
+## 2016-06-08, Weitian LI
+##   * Drop 'calc_distance' in favor of 'cosmo_calc'
 ##
-VERSION="v4.1"
-UPDATED="2015/06/03"
 
 ## error code {{{
 ERR_USG=1
@@ -69,26 +64,11 @@ case "$1" in
     -[hH]*|--[hH]*)
         printf "usage:\n"
         printf "    `basename $0` evt=<evt2_clean> r500=<r500_kpc> basedir=<basedir> info=<info_json> inner=<inner_val> outer=<outer_val> regin=<input_reg> regout=<output_reg> bkgd=<blank_evt|lbkg_reg|bkg_spec> nh=<nH> z=<redshift> arf=<arf_file> rmf=<rmf_file> [ grouptype=<NUM_CTS|BIN> grouptypeval=<number> binspec=<binspec> log=<log_file> ]\n"
-        printf "\nversion:\n"
-        printf "    ${VERSION}, ${UPDATED}\n"
         exit ${ERR_USG}
         ;;
 esac
 ## usage, help }}}
 
-## comology calculator {{{
-## XXX: MODIFY THIS TO YOUR OWN CASE
-## and make sure this `calc' is executable
-## NOTES: use `$HOME' instead of `~' in path
-COSCALC="`which cosmo_calc calc_distance 2>/dev/null | head -n 1`"
-# COSCALC="_path_to_calc_distance_"
-# COSCALC="$HOME/bin/mass/calc_distance"
-if [ -z "${COSCALC}" ] || [ ! -x ${COSCALC} ]; then
-    printf "ERROR: \`COSCALC: ${COSCALC}' neither specified nor executable\n"
-    exit 255
-fi
-## }}}
-
 ## default parameters {{{
 # default `event file' which used to match `blanksky' files
 #DFT_EVT="_NOT_EXIST_"
@@ -262,7 +242,7 @@ case "${R500_UNI}" in
         ;;
     *)
         printf "## units in \`kpc', convert to \`Chandra pixel'\n" | ${TOLOG}
-        KPC_PER_PIX=`${COSCALC} ${REDSHIFT} | \grep 'kpc.*pix' | tr -d 'a-zA-Z_#=(),:/ '`
+        KPC_PER_PIX=`cosmo_calc ${REDSHIFT} | \grep 'kpc.*pix' | tr -d 'a-zA-Z_#=(),:/ '`
         # convert scientific notation for `bc'
         KPC_PER_PIX_B=`echo ${KPC_PER_PIX} | sed 's/[eE]/\*10\^/' | sed 's/+//'`
         printf "## calculated \`kpc/pixel': ${KPC_PER_PIX_B}\n"
@@ -578,4 +558,3 @@ _EOF_
 ## xspec script }}}
 
 exit 0
-