Remove 'collect_infodata.sh' and 'json2csv.py'

Superseded by 'manifest.py' and 'results.py'

2 files changed, 0 insertions, 648 deletions

diff --git a/scripts/collect_infodata.sh b/scripts/collect_infodata.sh
deleted file mode 100755
index bdda7f7..0000000
--- a/scripts/collect_infodata.sh
+++ /dev/null
@@ -1,574 +0,0 @@
-#!/bin/sh
-##
-## Collect the calculated mass data and update them to the INFO json file.
-##
-## JSON parser: http://json.parser.online.fr/
-##
-## Weitian LI <liweitianux@gmail.com>
-## August 31, 2012
-##
-## Change logs:
-## 2017-02-17, Weitian LI
-##   * Remove cooling time and radius results calculated by 'cooling_time'
-## 2017-02-16, Weitian LI
-##   * Remove luminosity results calculated by 'calc_lx'
-## 2017-02-09, Weitian LI
-##   * Update to use the new style config files
-##   * Some cleanups
-## v3.4, 2015/06/03, Weitian LI
-##   * Copy needed pfiles to current working directory, and
-##     set environment variable $PFILES to use these first.
-##   * Replaced 'grep' with '\grep', 'ls' with '\ls'
-## v3.3, 2013/10/14, Weitian LI
-##   add key `Unified Name'
-## v3.2, 2013/05/29, Weitian LI
-##   add key `XCNTRD_RA, XCNTRD_DEC'
-## v3.1, 2013/05/18, Weitian LI
-##   add key `Feature'
-## v3.0, 2013/02/09, Weitian LI
-##   modified for `new sample info format'
-## v2.2, 2012/12/18, Weitian LI
-##   add `beta' and `cooling time' parameters
-## v2.1, 2012/11/07, Weitian LI
-##   account for `fitting_dbeta'
-## v2.0, 2012/10/14, Weitian LI
-##   add parameters
-## v1.1, 2012/09/05, Weitian LI
-##   add `T_avg(0.2-0.5 R500)' and `T_err'
-##
-
-## error code {{{
-ERR_USG=1
-ERR_DIR=11
-ERR_CFG=12
-ERR_RES=13
-ERR_COSC=14
-ERR_BETA=101
-ERR_BETA2=102
-ERR_JSON=201
-ERR_MLOG=202
-ERR_BLOG=203
-ERR_CLOG=204
-## error code }}}
-
-## usage, help {{{
-case "$1" in
-    -[hH]*|--[hH]*)
-        printf "usage:\n"
-        printf "    `basename $0` json=<info_json> cfg=<main_cfg> res=<final_result> basedir=<repro_dir> massdir=<mass_dir>\n"
-        exit ${ERR_USG}
-        ;;
-esac
-## usage, help }}}
-
-## default parameters {{{
-# default basedir
-DFT_BASEDIR=".."
-# default mass dir
-DFT_MASSDIR="mass"
-# default pattern for json info file
-DFT_JSON_PAT="*_INFO.json"
-# main config file
-DFT_CFG_PAT="mass.conf"
-# default result file
-DFT_RES_FINAL_PAT="final_result.txt"
-# default sbprofile region file
-DFT_SBPROFILE_REG="sbprofile.reg"
-# default radial spectra region file
-DFT_RSPEC_REG="rspec.reg"
-## default parameters }}}
-
-## functions {{{
-# process commandline arguments
-# cmdline arg format: `KEY=VALUE'
-getopt_keyval() {
-    until [ -z "$1" ]
-    do
-        key=${1%%=*}                    # extract key
-        val=${1#*=}                     # extract value
-        keyval="${key}=\"${val}\""
-        echo "## getopt: eval '${keyval}'"
-        eval ${keyval}
-        shift                           # shift, process next one
-    done
-}
-## functions }}}
-
-## parameters {{{
-# process cmdline args using `getopt_keyval'
-getopt_keyval "$@"
-
-# init directory
-INIT_DIR=`pwd -P`
-
-# check given parameters
-# check given dir
-if [ -d "${basedir}" ] && \ls ${basedir}/*repro_evt2.fits > /dev/null 2>&1; then
-    BASEDIR=${basedir}
-elif [ -d "${DFT_BASEDIR}" ] && \ls ${DFT_BASEDIR}/*repro_evt2.fits > /dev/null 2>&1; then
-    BASEDIR=${DFT_BASEDIR}
-else
-    read -p "> basedir (contains info json): " BASEDIR
-    if [ ! -d "${BASEDIR}" ]; then
-        printf "ERROR: given \`${BASEDIR}' NOT a directory\n"
-        exit ${ERR_DIR}
-    elif ! \ls ${BASEDIR}/*repro_evt2.fits > /dev/null 2>&1; then
-        printf "ERROR: given \`${BASEDIR}' NOT contains needed evt files\n"
-        exit ${ERR_DIR}
-    fi
-fi
-# remove the trailing '/'
-BASEDIR=`( cd ${INIT_DIR}/${BASEDIR} && pwd -P )`
-printf "## use basedir: \`${BASEDIR}'\n"
-# mass dir
-if [ ! -z "${massdir}" ] && [ -d "${BASEDIR}/${massdir}" ]; then
-    MASS_DIR=`( cd ${BASEDIR}/${massdir} && pwd -P )`
-elif [ -d "${BASEDIR}/${DFT_MASSDIR}" ]; then
-    MASS_DIR=`( cd ${BASEDIR}/${DFT_MASSDIR} && pwd -P )`
-else
-    read -p "> mass dir (relative to basedir): " MASS_DIR
-    if [ ! -d "${BASEDIR}/${MASS_DIR}" ]; then
-        printf "ERROR: given \`${BASEDIR}/${MASS_DIR}' NOT a directory\n"
-        exit ${ERR_DIR}
-    fi
-fi
-# remove the trailing '/'
-MASS_DIR=`echo ${MASS_DIR} | sed 's/\/*$//'`
-printf "## use massdir: \`${MASS_DIR}'\n"
-# check json file
-if [ ! -z "${json}" ] && [ -r "${BASEDIR}/${json}" ]; then
-    JSON_FILE="${json}"
-elif [ "`\ls ${BASEDIR}/${DFT_JSON_PAT} | wc -l`" -eq 1 ]; then
-    JSON_FILE=`( cd ${BASEDIR} && \ls ${DFT_JSON_PAT} )`
-else
-    read -p "> info json file: " JSON_FILE
-    if ! [ -r "${BASEDIR}/${JSON_FILE}" ]; then
-        printf "ERROR: cannot access given \`${BASEDIR}/${INFO_JSON}' file\n"
-        exit ${ERR_JSON}
-    fi
-fi
-printf "## use info json file: \`${JSON_FILE}'\n"
-# check main config file
-if [ ! -z "${cfg}" ] && [ -r "${cfg}" ]; then
-    CFG_FILE="${cfg}"
-elif [ -r "${DFT_CFG_PAT}" ]; then
-    CFG_FILE="${DFT_CFG_PAT}"
-else
-    read -p "> main config file: " CFG_FILE
-    if [ ! -r "${CFG_FILE}" ]; then
-        printf "ERROR: cannot access given \`${CFG_JSON}' file\n"
-        exit ${ERR_CFG}
-    fi
-fi
-printf "## use main config file: \`${CFG_FILE}'\n"
-SBP_CFG=`\grep '^sbp_cfg' ${CFG_FILE} | awk '{ print $2 }'`
-printf "## sbp config file: \`${SBP_CFG}'\n"
-# check final result file
-if [ ! -z "${res}" ] && [ -r "${res}" ]; then
-    RES_FINAL="${res}"
-elif [ -r "${DFT_RES_FINAL_PAT}" ]; then
-    RES_FINAL="${DFT_RES_FINAL_PAT}"
-else
-    read -p "> final result file: " RES_FINAL
-    if [ ! -r "${RES_FINAL}" ]; then
-        printf "ERROR: cannot access given \`${RES_FINAL}' file\n"
-        exit ${ERR_RES}
-    fi
-fi
-printf "## use final result file: \`${RES_FINAL}'\n"
-## parameters }}}
-
-## directory & file {{{
-BASE_PATH=`dirname $0`
-printf "## BASE_PATH: ${BASE_PATH}\n"
-EVT_DIR="${BASEDIR}/evt"
-SPEC_DIR="${BASEDIR}/spc/profile"
-## dir & file }}}
-
-## prepare parameter files (pfiles) {{{
-CIAO_TOOLS="dmkeypar"
-
-# Copy necessary pfiles for localized usage
-for tool in ${CIAO_TOOLS}; do
-    pfile=`paccess ${tool}`
-    [ -n "${pfile}" ] && punlearn ${tool} && cp -Lvf ${pfile} .
-done
-
-# Modify environment variable 'PFILES' to use local pfiles first
-export PFILES="./:${PFILES}"
-## pfiles }}}
-
-cd ${BASEDIR}
-printf "## enter directory: `pwd -P`\n"
-
-## in dir `repro' {{{
-punlearn dmkeypar
-EVT_RAW=`\ls *repro_evt2.fits`
-OBS_ID=`dmkeypar ${EVT_RAW} OBS_ID echo=yes`
-DATE_OBS=`dmkeypar ${EVT_RAW} DATE-OBS echo=yes`
-EXPOSURE_RAW=`dmkeypar ${EVT_RAW} EXPOSURE echo=yes | awk '{ print $1/1000 }'`
-## ACIS_TYPE
-DETNAM=`dmkeypar ${EVT_RAW} DETNAM echo=yes`
-if echo ${DETNAM} | \grep -q 'ACIS-0123'; then
-    ACIS_TYPE="ACIS-I"
-elif echo ${DETNAM} | \grep -q 'ACIS-[0-6]*7'; then
-    ACIS_TYPE="ACIS-S"
-else
-    printf "*** ERROR: unknown detector type: ${DETNAM}\n"
-    ACIS_TYPE="UNKNOWN"
-    exit 1
-fi
-## dir `repro' }}}
-
-## in dir `repro/evt' {{{
-cd ${EVT_DIR}
-EVT_CLEAN=`\ls evt2_c*_clean.fits`
-EXPOSURE_CLEAN=`dmkeypar ${EVT_CLEAN} EXPOSURE echo=yes | awk '{ print $1/1000 }'`
-## dir `repro/evt' }}}
-
-## in dir `repro/mass' {{{
-cd ${MASS_DIR}
-
-# misc {{{
-N_H=`\grep '^nh' ${CFG_FILE} | awk '{ print $2 }'`
-ABUND=`\grep '^abund' ${CFG_FILE} | awk '{ print $2 }'`
-TPROFILE_DATA=`\grep '^tprofile_data' ${CFG_FILE} | awk '{ print $2 }'`
-NFW_RMIN_KPC=`\grep '^nfw_rmin_kpc' ${CFG_FILE} | awk '{ print $2 }'`
-Z=`\grep '^z' ${SBP_CFG} | awk '{ print $2 }'`
-KPC_PER_PIXEL=`cosmo_calc -b --kpc-per-pix ${Z}`
-SBP_DATA=`\grep '^sbp_data' ${SBP_CFG} | awk '{ print $2 }'`
-RMAX_SBP_PIX=`tail -n 1 ${SBP_DATA} | awk '{ print $1+$2 }'`
-RMAX_SBP_KPC=`echo "${RMAX_SBP_PIX} ${KPC_PER_PIXEL}" | awk '{ printf("%.2f", $1*$2) }'`
-SPC_DIR="$(dirname $(readlink ${TPROFILE_DATA}))"
-if [ -f "${SPC_DIR}/${DFT_RSPEC_REG}" ]; then
-    RMAX_TPRO_PIX=`\grep -iE '(pie|annulus)' ${SPC_DIR}/${DFT_RSPEC_REG} | tail -n 1 | awk -F',' '{ print $4 }'`
-    RMAX_TPRO_KPC=`echo "${RMAX_TPRO_PIX} ${KPC_PER_PIXEL}" | awk '{ printf("%.2f", $1*$2) }'`
-fi
-[ -z "${NFW_RMIN_KPC}" ] && NFW_RMIN_KPC="null"
-[ -z "${RMAX_SBP_PIX}" ] && RMAX_SBP_PIX="null"
-[ -z "${RMAX_SBP_KPC}" ] && RMAX_SBP_KPC="null"
-# misc }}}
-
-## determine single/double beta {{{
-if \grep -q '^beta2' ${SBP_CFG}; then
-    MODEL_SBP="double-beta"
-    n01=`\grep '^n01' ${RES_FINAL} | awk '{ print $2 }'`
-    beta1=`\grep '^beta1' ${RES_FINAL} | awk '{ print $2 }'`
-    rc1=`\grep -E '^rc1\s' ${RES_FINAL} | awk '{ print $2 }'`
-    rc1_kpc=`\grep '^rc1_kpc' ${RES_FINAL} | awk '{ print $2 }'`
-    n02=`\grep '^n02' ${RES_FINAL} | awk '{ print $2 }'`
-    beta2=`\grep '^beta2' ${RES_FINAL} | awk '{ print $2 }'`
-    rc2=`\grep -E '^rc2\s' ${RES_FINAL} | awk '{ print $2 }'`
-    rc2_kpc=`\grep '^rc2_kpc' ${RES_FINAL} | awk '{ print $2 }'`
-    BKG=`\grep '^bkg' ${RES_FINAL} | awk '{ print $2 }'`
-    # beta1 -> smaller rc; beta2 -> bigger rc
-    if [ `echo "${rc1} < ${rc2}" | bc -l` -eq 1 ]; then
-        N01=${n01}
-        BETA1=${beta1}
-        RC1=${rc1}
-        RC1_KPC=${rc1_kpc}
-        N02=${n02}
-        BETA2=${beta2}
-        RC2=${rc2}
-        RC2_KPC=${rc2_kpc}
-    else
-        # beta1 <-> beta2 (swap)
-        N01=${n02}
-        BETA1=${beta2}
-        RC1=${rc2}
-        RC1_KPC=${rc2_kpc}
-        N02=${n01}
-        BETA2=${beta1}
-        RC2=${rc1}
-        RC2_KPC=${rc1_kpc}
-    fi
-else
-    MODEL_SBP="single-beta"
-    N01="null"
-    BETA1="null"
-    RC1="null"
-    RC1_KPC="null"
-    N02=`\grep '^n0' ${RES_FINAL} | awk '{ print $2 }'`
-    BETA2=`\grep '^beta' ${RES_FINAL} | awk '{ print $2 }'`
-    RC2=`\grep -E '^rc\s' ${RES_FINAL} | awk '{ print $2 }'`
-    RC2_KPC=`\grep '^rc_kpc' ${RES_FINAL} | awk '{ print $2 }'`
-    BKG=`\grep '^bkg' ${RES_FINAL} | awk '{ print $2 }'`
-fi
-## single/double beta }}}
-
-## get `mass/virial_radius' {{{
-# 200 data {{{
-R200_VAL=`\grep '^r200' ${RES_FINAL} | awk '{ print $2 }'`
-R200_ERR_L=`\grep '^r200' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $1 }'`
-R200_ERR_U=`\grep '^r200' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $2 }' | tr -d '+'`
-M200_VAL=`\grep '^m200' ${RES_FINAL} | awk '{ print $2 }'`
-M200_ERR_L=`\grep '^m200' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $1 }'`
-M200_ERR_U=`\grep '^m200' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $2 }' | tr -d '+'`
-MGAS200_VAL=`\grep '^gas_m200' ${RES_FINAL} | awk '{ print $2 }'`
-MGAS200_ERR_L=`\grep '^gas_m200' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $1 }'`
-MGAS200_ERR_U=`\grep '^gas_m200' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $2 }' | tr -d '+'`
-FGAS200_VAL=`\grep '^gas_fraction200' ${RES_FINAL} | awk '{ print $2 }'`
-FGAS200_ERR_L=`\grep '^gas_fraction200' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $1 }'`
-FGAS200_ERR_U=`\grep '^gas_fraction200' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $2 }' | tr -d '+'`
-[ -z "${R200_VAL}" ]       && R200_VAL="null"
-[ -z "${R200_ERR_L}" ]     && R200_ERR_L="null"
-[ -z "${R200_ERR_U}" ]     && R200_ERR_U="null"
-[ -z "${M200_VAL}" ]       && M200_VAL="null"
-[ -z "${M200_ERR_L}" ]     && M200_ERR_L="null"
-[ -z "${M200_ERR_U}" ]     && M200_ERR_U="null"
-[ -z "${MGAS200_VAL}" ]    && MGAS200_VAL="null"
-[ -z "${MGAS200_ERR_L}" ]  && MGAS200_ERR_L="null"
-[ -z "${MGAS200_ERR_U}" ]  && MGAS200_ERR_U="null"
-[ -z "${FGAS200_VAL}" ]    && FGAS200_VAL="null"
-[ -z "${FGAS200_ERR_L}" ]  && FGAS200_ERR_L="null"
-[ -z "${FGAS200_ERR_U}" ]  && FGAS200_ERR_U="null"
-# 200 }}}
-# 500 data {{{
-R500_VAL=`\grep '^r500' ${RES_FINAL} | awk '{ print $2 }'`
-R500_ERR_L=`\grep '^r500' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $1 }'`
-R500_ERR_U=`\grep '^r500' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $2 }' | tr -d '+'`
-M500_VAL=`\grep '^m500' ${RES_FINAL} | awk '{ print $2 }'`
-M500_ERR_L=`\grep '^m500' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $1 }'`
-M500_ERR_U=`\grep '^m500' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $2 }' | tr -d '+'`
-MGAS500_VAL=`\grep '^gas_m500' ${RES_FINAL} | awk '{ print $2 }'`
-MGAS500_ERR_L=`\grep '^gas_m500' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $1 }'`
-MGAS500_ERR_U=`\grep '^gas_m500' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $2 }' | tr -d '+'`
-FGAS500_VAL=`\grep '^gas_fraction500' ${RES_FINAL} | awk '{ print $2 }'`
-FGAS500_ERR_L=`\grep '^gas_fraction500' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $1 }'`
-FGAS500_ERR_U=`\grep '^gas_fraction500' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $2 }' | tr -d '+'`
-[ -z "${R500_VAL}" ]       && R500_VAL="null"
-[ -z "${R500_ERR_L}" ]     && R500_ERR_L="null"
-[ -z "${R500_ERR_U}" ]     && R500_ERR_U="null"
-[ -z "${M500_VAL}" ]       && M500_VAL="null"
-[ -z "${M500_ERR_L}" ]     && M500_ERR_L="null"
-[ -z "${M500_ERR_U}" ]     && M500_ERR_U="null"
-[ -z "${MGAS500_VAL}" ]    && MGAS500_VAL="null"
-[ -z "${MGAS500_ERR_L}" ]  && MGAS500_ERR_L="null"
-[ -z "${MGAS500_ERR_U}" ]  && MGAS500_ERR_U="null"
-[ -z "${FGAS500_VAL}" ]    && FGAS500_VAL="null"
-[ -z "${FGAS500_ERR_L}" ]  && FGAS500_ERR_L="null"
-[ -z "${FGAS500_ERR_U}" ]  && FGAS500_ERR_U="null"
-# 500 }}}
-# 1500 data {{{
-R1500_VAL=`\grep '^r1500' ${RES_FINAL} | awk '{ print $2 }'`
-R1500_ERR_L=`\grep '^r1500' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $1 }'`
-R1500_ERR_U=`\grep '^r1500' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $2 }' | tr -d '+'`
-M1500_VAL=`\grep '^m1500' ${RES_FINAL} | awk '{ print $2 }'`
-M1500_ERR_L=`\grep '^m1500' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $1 }'`
-M1500_ERR_U=`\grep '^m1500' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $2 }' | tr -d '+'`
-MGAS1500_VAL=`\grep '^gas_m1500' ${RES_FINAL} | awk '{ print $2 }'`
-MGAS1500_ERR_L=`\grep '^gas_m1500' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $1 }'`
-MGAS1500_ERR_U=`\grep '^gas_m1500' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $2 }' | tr -d '+'`
-FGAS1500_VAL=`\grep '^gas_fraction1500' ${RES_FINAL} | awk '{ print $2 }'`
-FGAS1500_ERR_L=`\grep '^gas_fraction1500' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $1 }'`
-FGAS1500_ERR_U=`\grep '^gas_fraction1500' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $2 }' | tr -d '+'`
-[ -z "${R1500_VAL}" ]      && R1500_VAL="null"
-[ -z "${R1500_ERR_L}" ]    && R1500_ERR_L="null"
-[ -z "${R1500_ERR_U}" ]    && R1500_ERR_U="null"
-[ -z "${M1500_VAL}" ]      && M1500_VAL="null"
-[ -z "${M1500_ERR_L}" ]    && M1500_ERR_L="null"
-[ -z "${M1500_ERR_U}" ]    && M1500_ERR_U="null"
-[ -z "${MGAS1500_VAL}" ]   && MGAS1500_VAL="null"
-[ -z "${MGAS1500_ERR_L}" ] && MGAS1500_ERR_L="null"
-[ -z "${MGAS1500_ERR_U}" ] && MGAS1500_ERR_U="null"
-[ -z "${FGAS1500_VAL}" ]   && FGAS1500_VAL="null"
-[ -z "${FGAS1500_ERR_L}" ] && FGAS1500_ERR_L="null"
-[ -z "${FGAS1500_ERR_U}" ] && FGAS1500_ERR_U="null"
-# 1500 }}}
-# 2500 data {{{
-R2500_VAL=`\grep '^r2500' ${RES_FINAL} | awk '{ print $2 }'`
-R2500_ERR_L=`\grep '^r2500' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $1 }'`
-R2500_ERR_U=`\grep '^r2500' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $2 }' | tr -d '+'`
-M2500_VAL=`\grep '^m2500' ${RES_FINAL} | awk '{ print $2 }'`
-M2500_ERR_L=`\grep '^m2500' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $1 }'`
-M2500_ERR_U=`\grep '^m2500' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $2 }' | tr -d '+'`
-MGAS2500_VAL=`\grep '^gas_m2500' ${RES_FINAL} | awk '{ print $2 }'`
-MGAS2500_ERR_L=`\grep '^gas_m2500' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $1 }'`
-MGAS2500_ERR_U=`\grep '^gas_m2500' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $2 }' | tr -d '+'`
-FGAS2500_VAL=`\grep '^gas_fraction2500' ${RES_FINAL} | awk '{ print $2 }'`
-FGAS2500_ERR_L=`\grep '^gas_fraction2500' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $1 }'`
-FGAS2500_ERR_U=`\grep '^gas_fraction2500' ${RES_FINAL} | awk '{ print $3 }' | awk -F'/' '{ print $2 }' | tr -d '+'`
-[ -z "${R2500_VAL}" ]      && R2500_VAL="null"
-[ -z "${R2500_ERR_L}" ]    && R2500_ERR_L="null"
-[ -z "${R2500_ERR_U}" ]    && R2500_ERR_U="null"
-[ -z "${M2500_VAL}" ]      && M2500_VAL="null"
-[ -z "${M2500_ERR_L}" ]    && M2500_ERR_L="null"
-[ -z "${M2500_ERR_U}" ]    && M2500_ERR_U="null"
-[ -z "${MGAS2500_VAL}" ]   && MGAS2500_VAL="null"
-[ -z "${MGAS2500_ERR_L}" ] && MGAS2500_ERR_L="null"
-[ -z "${MGAS2500_ERR_U}" ] && MGAS2500_ERR_U="null"
-[ -z "${FGAS2500_VAL}" ]   && FGAS2500_VAL="null"
-[ -z "${FGAS2500_ERR_L}" ] && FGAS2500_ERR_L="null"
-[ -z "${FGAS2500_ERR_U}" ] && FGAS2500_ERR_U="null"
-# 2500 }}}
-FGRR=`\grep '^gas_fraction.*r2500.*r500=' ${RES_FINAL} | sed 's/^.*r500=//' | awk '{ print $1 }'`
-FGRR_ERR_L=`\grep '^gas_fraction.*r2500.*r500=' ${RES_FINAL} | sed 's/^.*r500=//' | awk '{ print $2 }' | awk -F'/' '{ print $1 }'`
-FGRR_ERR_U=`\grep '^gas_fraction.*r2500.*r500=' ${RES_FINAL} | sed 's/^.*r500=//' | awk '{ print $2 }' | awk -F'/' '{ print $2 }' | tr -d '+'`
-[ -z "${FGRR}" ]       && FGRR="null"
-[ -z "${FGRR_ERR_L}" ] && FGRR_ERR_L="null"
-[ -z "${FGRR_ERR_U}" ] && FGRR_ERR_U="null"
-## mrl }}}
-## repro/mass }}}
-
-cd ${BASEDIR}
-## orig json file {{{
-printf "## collect data from original info file ...\n"
-OBJ_NAME=`\grep '"Source\ Name' ${JSON_FILE} | sed 's/.*"Source.*":\ //' | sed 's/^"//' | sed 's/"\ *,$//'`
-OBJ_UNAME=`\grep '"Unified\ Name' ${JSON_FILE} | sed 's/.*"Unified.*":\ //' | sed 's/^"//' | sed 's/"\ *,$//'`
-OBJ_RA=`\grep '"R\.\ A\.' ${JSON_FILE} | sed 's/.*"R\.\ A\.":\ //' | sed 's/^"//' | sed 's/"\ *,$//'`
-OBJ_DEC=`\grep '"Dec\.' ${JSON_FILE} | sed 's/.*"Dec\.":\ //' | sed 's/^"//' | sed 's/"\ *,$//'`
-OBJ_XCRA=`\grep '"XCNTRD_RA' ${JSON_FILE} | sed 's/.*"XCNTRD_RA":\ //' | sed 's/^"//' | sed 's/"\ *,$//'`
-OBJ_XCDEC=`\grep '"XCNTRD_DEC' ${JSON_FILE} | sed 's/.*"XCNTRD_DEC":\ //' | sed 's/^"//' | sed 's/"\ *,$//'`
-REDSHIFT=`\grep '"redshift' ${JSON_FILE} | sed 's/.*"redshift.*":\ //' | sed 's/\ *,$//'`
-COOLCORE=`\grep -i '"cool.*core' ${JSON_FILE} | sed 's/.*"[cC]ool.*":\ //' | sed 's/\ *,$//'`
-[ -z "${COOLCORE}" ] && COOLCORE="null"
-OBJ_FEATURE=`\grep '"Feature' ${JSON_FILE} | sed 's/.*"Feature":\ //' | sed 's/^"//' | sed 's/"\ *,*$//'`
-OBJ_NOTE=`\grep '"NOTE' ${JSON_FILE} | sed 's/.*"NOTE":\ //' | sed 's/^"//' | sed 's/"\ *,*$//'`
-
-# T & Z {{{
-T_1R500=`\grep '"T(0\.1.*R500' ${JSON_FILE} | sed 's/.*"T.*":\ //' | sed 's/\ *,$//'`
-T_1ERR=`\grep '"T_err(.*0\.1.*R500' ${JSON_FILE} | sed 's/.*"T.*":\ //' | sed 's/\ *,$//'`
-T_1ERR_L=`\grep '"T_err_l.*0\.1.*R500' ${JSON_FILE} | sed 's/.*"T.*":\ //' | sed 's/\ *,$//'`
-T_1ERR_U=`\grep '"T_err_u.*0\.1.*R500' ${JSON_FILE} | sed 's/.*"T.*":\ //' | sed 's/\ *,$//'`
-Z_1R500=`\grep '"Z(0\.1.*R500' ${JSON_FILE} | sed 's/.*"Z.*":\ //' | sed 's/\ *,$//'`
-Z_1ERR=`\grep '"Z_err(.*0\.1.*R500' ${JSON_FILE} | sed 's/.*"Z.*":\ //' | sed 's/\ *,$//'`
-Z_1ERR_L=`\grep '"Z_err_l.*0\.1.*R500' ${JSON_FILE} | sed 's/.*"Z.*":\ //' | sed 's/\ *,$//'`
-Z_1ERR_U=`\grep '"Z_err_u.*0\.1.*R500' ${JSON_FILE} | sed 's/.*"Z.*":\ //' | sed 's/\ *,$//'`
-T_2R500=`\grep '"T(0\.2.*R500' ${JSON_FILE} | sed 's/.*"T.*":\ //' | sed 's/\ *,$//'`
-T_2ERR=`\grep '"T_err(.*0\.2.*R500' ${JSON_FILE} | sed 's/.*"T.*":\ //' | sed 's/\ *,$//'`
-T_2ERR_L=`\grep '"T_err_l.*0\.2.*R500' ${JSON_FILE} | sed 's/.*"T.*":\ //' | sed 's/\ *,$//'`
-T_2ERR_U=`\grep '"T_err_u.*0\.2.*R500' ${JSON_FILE} | sed 's/.*"T.*":\ //' | sed 's/\ *,$//'`
-Z_2R500=`\grep '"Z(0\.2.*R500' ${JSON_FILE} | sed 's/.*"Z.*":\ //' | sed 's/\ *,$//'`
-Z_2ERR=`\grep '"Z_err(.*0\.2.*R500' ${JSON_FILE} | sed 's/.*"Z.*":\ //' | sed 's/\ *,$//'`
-Z_2ERR_L=`\grep '"Z_err_l.*0\.2.*R500' ${JSON_FILE} | sed 's/.*"Z.*":\ //' | sed 's/\ *,$//'`
-Z_2ERR_U=`\grep '"Z_err_u.*0\.2.*R500' ${JSON_FILE} | sed 's/.*"Z.*":\ //' | sed 's/\ *,$//'`
-[ -z "${T_1R500}" ]  && T_1R500="null"
-[ -z "${T_1ERR}" ]   && T_1ERR="null"
-[ -z "${T_1ERR_L}" ] && T_1ERR_L="null"
-[ -z "${T_1ERR_U}" ] && T_1ERR_U="null"
-[ -z "${Z_1R500}" ]  && Z_1R500="null"
-[ -z "${Z_1ERR}" ]   && Z_1ERR="null"
-[ -z "${Z_1ERR_L}" ] && Z_1ERR_L="null"
-[ -z "${Z_1ERR_U}" ] && Z_1ERR_U="null"
-[ -z "${T_2R500}" ]  && T_2R500="null"
-[ -z "${T_2ERR}" ]   && T_2ERR="null"
-[ -z "${T_2ERR_L}" ] && T_2ERR_L="null"
-[ -z "${T_2ERR_U}" ] && T_2ERR_U="null"
-[ -z "${Z_2R500}" ]  && Z_2R500="null"
-[ -z "${Z_2ERR}" ]   && Z_2ERR="null"
-[ -z "${Z_2ERR_L}" ] && Z_2ERR_L="null"
-[ -z "${Z_2ERR_U}" ] && Z_2ERR_U="null"
-# T & Z }}}
-## json data }}}
-
-## output data to JSON file {{{
-printf "## save collected data into file \`${JSON_FILE}' ...\n"
-mv -fv ${JSON_FILE} ${JSON_FILE}_bak
-cat > ${JSON_FILE} << _EOF_
-{
-    "Obs. ID": ${OBS_ID},
-    "Source Name": "${OBJ_NAME}",
-    "Unified Name": "${OBJ_UNAME}",
-    "Obs. Date": "${DATE_OBS}",
-    "Detector": "${ACIS_TYPE}",
-    "Exposure (ks)": ${EXPOSURE_RAW},
-    "Clean Exposure (ks)": ${EXPOSURE_CLEAN},
-    "R. A.": "${OBJ_RA}",
-    "Dec.": "${OBJ_DEC}",
-    "XCNTRD_RA": "${OBJ_XCRA}",
-    "XCNTRD_DEC": "${OBJ_XCDEC}",
-    "nH (10^22 cm^-2)": ${N_H},
-    "redshift": ${REDSHIFT},
-    "E(z)": null,
-    "T_ref (keV)": null,
-    "Z_ref (solar)": ${ABUND},
-    "Rmax_SBP (pixel)": ${RMAX_SBP_PIX},
-    "Rmax_Tpro (pixel)": ${RMAX_TPRO_PIX},
-    "Rmax_SBP (kpc)": ${RMAX_SBP_KPC},
-    "Rmax_Tpro (kpc)": ${RMAX_TPRO_KPC},
-    "NFW_Rmin (kpc)": ${NFW_RMIN_KPC},
-    "Model_SBP": "${MODEL_SBP}",
-    "n01": ${N01},
-    "beta1": ${BETA1},
-    "rc1": ${RC1},
-    "rc1_kpc": ${RC1_KPC},
-    "n02": ${N02},
-    "beta2": ${BETA2},
-    "rc2": ${RC2},
-    "rc2_kpc": ${RC2_KPC},
-    "bkg": ${BKG},
-    "R200 (kpc)": ${R200_VAL},
-    "R200_err_lower (1sigma)": ${R200_ERR_L},
-    "R200_err_upper (1sigma)": ${R200_ERR_U},
-    "M200 (M_sun)": ${M200_VAL},
-    "M200_err_lower (1sigma)": ${M200_ERR_L},
-    "M200_err_upper (1sigma)": ${M200_ERR_U},
-    "M_gas200 (M_sun)": ${MGAS200_VAL},
-    "M_gas200_err_lower (1sigma)": ${MGAS200_ERR_L},
-    "M_gas200_err_upper (1sigma)": ${MGAS200_ERR_U},
-    "F_gas200": ${FGAS200_VAL},
-    "F_gas200_err_lower (1sigma)": ${FGAS200_ERR_L},
-    "F_gas200_err_upper (1sigma)": ${FGAS200_ERR_U},
-    "R500 (kpc)": ${R500_VAL},
-    "R500_err_lower (1sigma)": ${R500_ERR_L},
-    "R500_err_upper (1sigma)": ${R500_ERR_U},
-    "M500 (M_sun)": ${M500_VAL},
-    "M500_err_lower (1sigma)": ${M500_ERR_L},
-    "M500_err_upper (1sigma)": ${M500_ERR_U},
-    "M_gas500 (M_sun)": ${MGAS500_VAL},
-    "M_gas500_err_lower (1sigma)": ${MGAS500_ERR_L},
-    "M_gas500_err_upper (1sigma)": ${MGAS500_ERR_U},
-    "F_gas500": ${FGAS500_VAL},
-    "F_gas500_err_lower (1sigma)": ${FGAS500_ERR_L},
-    "F_gas500_err_upper (1sigma)": ${FGAS500_ERR_U},
-    "R1500": ${R1500_VAL},
-    "R1500_err_lower": ${R1500_ERR_L},
-    "R1500_err_upper": ${R1500_ERR_U},
-    "M1500": ${M1500_VAL},
-    "M1500_err_lower": ${M1500_ERR_L},
-    "M1500_err_upper": ${M1500_ERR_U},
-    "M_gas1500": ${MGAS1500_VAL},
-    "M_gas1500_err_lower": ${MGAS1500_ERR_L},
-    "M_gas1500_err_upper": ${MGAS1500_ERR_U},
-    "F_gas1500": ${FGAS1500_VAL},
-    "F_gas1500_err_lower": ${FGAS1500_ERR_L},
-    "F_gas1500_err_upper": ${FGAS1500_ERR_U},
-    "R2500": ${R2500_VAL},
-    "R2500_err_lower": ${R2500_ERR_L},
-    "R2500_err_upper": ${R2500_ERR_U},
-    "M2500": ${M2500_VAL},
-    "M2500_err_lower": ${M2500_ERR_L},
-    "M2500_err_upper": ${M2500_ERR_U},
-    "M_gas2500": ${MGAS2500_VAL},
-    "M_gas2500_err_lower": ${MGAS2500_ERR_L},
-    "M_gas2500_err_upper": ${MGAS2500_ERR_U},
-    "F_gas2500": ${FGAS2500_VAL},
-    "F_gas2500_err_lower": ${FGAS2500_ERR_L},
-    "F_gas2500_err_upper": ${FGAS2500_ERR_U},
-    "T(0.1-0.5 R500)": ${T_1R500},
-    "T_err(0.1-0.5 R500)": ${T_1ERR},
-    "T_err_l(0.1-0.5 R500)": ${T_1ERR_L},
-    "T_err_u(0.1-0.5 R500)": ${T_1ERR_U},
-    "Z(0.1-0.5 R500)": ${Z_1R500},
-    "Z_err(0.1-0.5 R500)": ${Z_1ERR},
-    "Z_err_l(0.1-0.5 R500)": ${Z_1ERR_L},
-    "Z_err_u(0.1-0.5 R500)": ${Z_1ERR_U},
-    "T(0.2-0.5 R500)": ${T_2R500},
-    "T_err(0.2-0.5 R500)": ${T_2ERR},
-    "T_err_l(0.2-0.5 R500)": ${T_2ERR_L},
-    "T_err_u(0.2-0.5 R500)": ${T_2ERR_U},
-    "Z(0.2-0.5 R500)": ${Z_2R500},
-    "Z_err(0.2-0.5 R500)": ${Z_2ERR},
-    "Z_err_l(0.2-0.5 R500)": ${Z_2ERR_L},
-    "Z_err_u(0.2-0.5 R500)": ${Z_2ERR_U},
-    "F_gas(R2500-R500)": ${FGRR},
-    "F_gas_err_l(R2500-R500)": ${FGRR_ERR_L},
-    "F_gas_err_u(R2500-R500)": ${FGRR_ERR_U},
-    "Cool_core": ${COOLCORE},
-    "Feature": "${OBJ_FEATURE}",
-    "NOTE": "${OBJ_NOTE}"
-},
-_EOF_
-## output JSON }}}
-
-exit 0
diff --git a/scripts/json2csv.py b/scripts/json2csv.py
deleted file mode 100755
index a06a135..0000000
--- a/scripts/json2csv.py
+++ /dev/null
@@ -1,74 +0,0 @@
-#!/usr/bin/env python3
-#
-# Copyright (c) 2012-2017 Weitian LI <liweitianux@live.com>
-# MIT license
-#
-# Convert JSON to CSV format
-#
-# References:
-# [1] https://docs.python.org/3/library/json.html
-# [2] https://docs.python.org/3/library/csv.html
-# [3] http://stackoverflow.com/questions/1871524/convert-from-json-to-csv-using-python
-# [4] http://json.parser.online.fr/
-#
-#
-# Weitian LI <liweitianux@live.com>
-# Created: 2012-08-31
-#
-# Change logs:
-# 2017-02-06, Weitian LI
-#   * Use `argparse`
-#   * Cleanups
-# v3.5, 2015-11-09
-#   * Get column keys from the first input json block
-#   * Use python3
-# v3.4, 2015-11-09
-#   * Add keys 'XPEAK_RA', 'XPEAK_DEC', and 'XPEAK_XCNTRD_dist (pix)'
-#   * Add 'colkeys' to record CSV column keys
-#   * Update file header description
-# v3.3, 2013-10-14
-#   * add key `Unified Name'
-# v3.2, 2013-05-29
-#   * add key `XCNTRD_RA, XCNTRD_DEC'
-# v3.1, 2013-05-18
-#   * add key `Feature', corresponding to `collectdata_v3.1'
-
-import os
-import argparse
-import csv
-import json
-from collections import OrderedDict
-
-
-def main():
-    parser = argparse.ArgumentParser(description="Convert JSON to CSV")
-    parser.add_argument("infile", help="Input JSON file")
-    parser.add_argument("outfile", nargs="?",
-                        help="Output CSV file (default to use the same " +
-                        "basename as the input file)")
-    args = parser.parse_args()
-
-    with open(args.infile) as f:
-        # Use `OrderedDict` to keep the orders of keys of all json files
-        data = json.load(f, object_pairs_hook=OrderedDict)
-
-    # Column keys
-    colkeys = list(data[0].keys())
-
-    if args.outfile:
-        outfile = args.outfile
-    else:
-        outfile = os.path.splitext(args.infile)[0] + ".csv"
-    with open(outfile, "w") as csvfile:
-        outfile_writer = csv.writer(csvfile)
-        # CSV header row
-        outfile_writer.writerow(colkeys)
-        # CSV data rows
-        for row in data:
-            outfile_writer.writerow([row.get(key) for key in colkeys])
-
-
-if __name__ == "__main__":
-    main()
-
-#  vim: set ts=4 sw=4 tw=0 fenc=utf-8 ft=python: #