diff options
Diffstat (limited to 'scripts')
| -rwxr-xr-x | scripts/clean_massdir.sh | 219 |
1 files changed, 0 insertions, 219 deletions
diff --git a/scripts/clean_massdir.sh b/scripts/clean_massdir.sh deleted file mode 100755 index dc9c4d7..0000000 --- a/scripts/clean_massdir.sh +++ /dev/null @@ -1,219 +0,0 @@ -#!/bin/sh -# -# clean files in `mass' dir -# -# v2, 2013/05/04, LIweitiaNux -# available for `empty mass dir' -# - -IMG_DIR=${1:-../img} -SPC_DIR=${2:-../spc/profile} - -if [ ! -f fitting_mass.conf ] || [ ! -f fitting_dbeta_mass.conf ]; then - # empty mass dir - printf "## EMPTY mass dir ...\n" - rm -rf *.* - - # expcorr_conf - EXPCORR_CONF=`ls ${IMG_DIR}/*_expcorr.conf 2> /dev/null` - if [ ! -z "${EXPCORR_CONF}" ]; then - Z=`grep '^z' ${EXPCORR_CONF} | awk '{ print $2 }'` - N_H=`grep '^nh' ${EXPCORR_CONF} | awk '{ print $2 }'` - ABUND=`grep '^abund' ${EXPCORR_CONF} | awk '{ print $2 }'` - printf "# redshift: ${Z}\n" - printf "# nh: ${N_H}\n" - printf "# abund: ${ABUND}\n" - fi - - # cosmo_calc - COSMO_CALC=`which cosmo_calc 2> /dev/null` - if [ ! -z "${COSMO_CALC}" ]; then - CM_PER_PIXEL=`${COSMO_CALC} ${Z} | grep 'cm/pixel' | awk -F':' '{ print $2 }' | tr -d ' '` - printf "# cm_per_pixel: ${CM_PER_PIXEL}\n" - fi - - cat > wang2012_param.txt << _EOF_ -A 5.0 1.0 500 T -n 5.0 0.1 10 T -xi 0.3 0.1 1.0 T -a2 2000 1000 1e+05 T -a3 1000 600 3000 T -beta 0.5 0.1 0.5 T -T0 1.0 1.0 2.0 T -_EOF_ - - cat > fitting_mass.conf << _EOF_ -t_profile wang2012 -t_data_file tcl_temp_profile.txt -t_param_file wang2012_param.txt -sbp_cfg fitting_sbp.conf -nh ${N_H} -abund ${ABUND} -radius_sbp_file sbprofile.txt -# nfw_rmin_kpc 1 -_EOF_ - - cat > fitting_dbeta_mass.conf << _EOF_ -t_profile wang2012 -t_data_file tcl_temp_profile.txt -t_param_file wang2012_param.txt -sbp_cfg fitting_dbeta_sbp.conf -nh ${N_H} -abund ${ABUND} -radius_sbp_file sbprofile.txt -# nfw_rmin_kpc 1 -_EOF_ - - cat > fitting_sbp.conf << _EOF_ -radius_file radius_sbp.txt -sbp_file flux_sbp.txt - -cfunc_file coolfunc_calc_data.txt -T_file t_profile_dump.qdp - -n0 0.005 -rc 30 -beta 0.7 -bkg 0 - -cm_per_pixel ${CM_PER_PIXEL} -z ${Z} -_EOF_ - - cat > fitting_dbeta_sbp.conf << _EOF_ -radius_file radius_sbp.txt -sbp_file flux_sbp.txt - -cfunc_file coolfunc_calc_data.txt -T_file t_profile_dump.qdp - -n01 0.05 -rc1 30 -beta1 0.7 -n02 0.005 -rc2 300 -beta2 0.7 -bkg 0 - -cm_per_pixel ${CM_PER_PIXEL} -z ${Z} -_EOF_ - - # link files - [ -f ${IMG_DIR}/flux_sbp.txt ] && ln -svf ${IMG_DIR}/flux_sbp.txt . - [ -f ${IMG_DIR}/radius_sbp.txt ] && ln -svf ${IMG_DIR}/radius_sbp.txt . - [ -f ${IMG_DIR}/sbprofile.txt ] && ln -svf ${IMG_DIR}/sbprofile.txt . - [ -f ${SPC_DIR}/tcl_temp_profile.txt ] && ln -svf ${SPC_DIR}/tcl_temp_profile.txt . - exit 0 -fi - -######################################################################## -rm -rf _* -rm -rf *?backup* -rm -rf *backup?* -rm *_bak *.log -rm global.cfg flux_sbp.txt radius_sbp.txt sbprofile.txt -rm tcl_temp_profile.qdp tcl_temp_profile.txt - -mkdir backup -cp fitting_mass.conf fitting_sbp.conf backup/ -cp fitting_dbeta_mass.conf fitting_dbeta_sbp.conf backup/ -cp beta_param_center.txt dbeta_param_center.txt backup/ -cp wang2012_param.txt backup/ -cp results_mrl.txt final_result.txt backup/ - -rm *.* - -[ -f backup/fitting_sbp.conf ] && cp backup/fitting_sbp.conf . -[ -f backup/wang2012_param.txt ] && cp backup/wang2012_param.txt . -[ -f fitting_mass.conf ] && cp fitting_mass.conf fitting_dbeta_mass.conf -[ -f fitting_sbp.conf ] && cp fitting_sbp.conf fitting_dbeta_sbp.conf - -SBP_CONF=`ls backup/fitting_sbp.conf backup/fitting_dbeta_sbp.conf 2>/dev/null | head -n 1` -EXPCORR_CONF=`ls ${IMG_DIR}/*_expcorr.conf` -N_H=`grep '^nh' ${EXPCORR_CONF} | awk '{ print $2 }'` - -sed -i'' "s/^beta.*$/beta 0.5 0.1 0.7 T/" wang2012_param.txt - -if [ -f backup/fitting_mass.conf ]; then - # mass_conf - cp backup/fitting_mass.conf . - sed -i'' "s/^t_data_file.*$/t_data_file tcl_temp_profile.txt/" fitting_mass.conf - sed -i'' "s/^nh.*$/nh ${N_H}/" fitting_mass.conf - cp fitting_mass.conf fitting_dbeta_mass.conf - sed -i'' "s/fitting_sbp/fitting_dbeta_sbp/" fitting_dbeta_mass.conf - # sbp_conf - cp backup/fitting_sbp.conf . - SBP_CONF="fitting_sbp.conf" - SBP2_CONF="fitting_dbeta_sbp.conf" -elif [ -f backup/fitting_dbeta_mass.conf ]; then - # mass_conf - cp backup/fitting_mass.conf . - cp backup/fitting_dbeta_mass.conf . - sed -i'' "s/^t_data_file.*$/t_data_file tcl_temp_profile.txt/" fitting_dbeta_mass.conf - sed -i'' "s/^nh.*$/nh ${N_H}/" fitting_dbeta_mass.conf - cp fitting_dbeta_mass.conf fitting_mass.conf - sed -i'' "s/fitting_dbeta_sbp/fitting_sbp/" fitting_mass.conf - # sbp_conf - cp backup/fitting_dbeta_sbp.conf . - SBP_CONF="fitting_dbeta_sbp.conf" - SBP2_CONF="fitting_sbp.conf" -else - # - printf "*** ERROR: fitting_mass.conf & fitting_dbeta_mass.conf not exists ***\n" -fi - -radius_file=`grep '^radius_file' $SBP_CONF | awk '{ print $2 }'` -sbp_file=`grep '^sbp_file' $SBP_CONF | awk '{ print $2 }'` -cfunc_file=`grep '^cfunc_file' $SBP_CONF | awk '{ print $2 }'` -T_file=`grep '^T_file' $SBP_CONF | awk '{ print $2 }'` -cm_per_pixel=`grep '^cm_per_pixel' $SBP_CONF | awk '{ print $2 }'` -z=`grep '^z' $SBP_CONF | awk '{ print $2 }'` - -if [ "x${SBP2_CONF}" = "xfitting_sbp.conf" ]; then - rm -f ${SBP2_CONF} - cat > ${SBP2_CONF} << _EOF_ -radius_file $radius_file -sbp_file $sbp_file - -cfunc_file $cfunc_file -T_file $T_file - -n0 0.005 -rc 30 -beta 0.7 -bkg 0 - -cm_per_pixel $cm_per_pixel -z $z -_EOF_ -elif [ "x${SBP2_CONF}" = "xfitting_dbeta_sbp.conf" ]; then - rm -f ${SBP2_CONF} - cat > ${SBP2_CONF} << _EOF_ -radius_file $radius_file -sbp_file $sbp_file - -cfunc_file $cfunc_file -T_file $T_file - -n01 0.05 -rc1 30 -beta1 0.7 -n02 0.005 -rc2 300 -beta2 0.7 -bkg 0 - -cm_per_pixel $cm_per_pixel -z $z -_EOF_ -else - # - printf "*** ERROR ***\n" -fi - -[ -f ${IMG_DIR}/flux_sbp.txt ] && ln -sf ${IMG_DIR}/flux_sbp.txt . -[ -f ${IMG_DIR}/radius_sbp.txt ] && ln -sf ${IMG_DIR}/radius_sbp.txt . -[ -f ${IMG_DIR}/sbprofile.txt ] && ln -sf ${IMG_DIR}/sbprofile.txt . -[ -f ${SPC_DIR}/tcl_temp_profile.txt ] && ln -sf ${SPC_DIR}/tcl_temp_profile.txt . - |