From 3c1fb2657e18a33e2a44bef03a4d8f310e1141c3 Mon Sep 17 00:00:00 2001
From: Aaron LI <aaronly.me@outlook.com>
Date: Thu, 16 Feb 2017 22:03:55 +0800
Subject: Update variable names to be more clear

---
 mass_profile/coolfunc_calc.sh     | 27 ++++++++++++---------------
 mass_profile/coolfunc_calc_erg.sh | 19 +++++++++----------
 2 files changed, 21 insertions(+), 25 deletions(-)

(limited to 'mass_profile')

diff --git a/mass_profile/coolfunc_calc.sh b/mass_profile/coolfunc_calc.sh
index 66e9123..5077b29 100755
--- a/mass_profile/coolfunc_calc.sh
+++ b/mass_profile/coolfunc_calc.sh
@@ -21,7 +21,7 @@ else
     exit 1
 fi
 TPROFILE=$1
-ABUND_VAL=$2
+ABUNDANCE=$2
 N_H=$3
 REDSHIFT=$4
 COOLFUNC_DAT=$5
@@ -55,7 +55,7 @@ set xs_echo_script 0
 ## set basic data {{{
 set nh ${N_H}
 set redshift ${REDSHIFT}
-set abund_val ${ABUND_VAL}
+set abundance ${ABUNDANCE}
 set norm ${NORM}
 ## basic }}}
 
@@ -67,7 +67,7 @@ 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 & \${abund_val} & \${redshift} & \${norm} & /*
+model wabs*apec & \${nh} & 1.0 & \${abundance} & \${redshift} & \${norm} & /*
 ## xspec }}}
 
 ## set input and output filename & open files
@@ -103,28 +103,25 @@ fi
 cat >> ${XSPEC_CF_XCM} << _EOF_
 ## 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}"
+    scan \${tpro_line} "%f %f" radius temperature
+    #puts "radius: \${radius}, temperature: \${temperature}"
     # set temperature value
-    newpar 2 \${temp_val}
+    newpar 2 \${temperature}
     # calc flux & tclout
     flux 0.7 7.0
     tclout flux 1
-    scan \${xspec_tclout} "%f %f %f %f" holder holder holder cf_data
-    #puts "cf_data: \${cf_data}"
-    puts \${cf_fd} "\${radius}    \${cf_data}"
+    scan \${xspec_tclout} "%f %f %f %f" _ _ _ cf_photon
+    #puts "cf: \${cf_photon}"
+    puts \${cf_fd} "\${radius}    \${cf_photon}"
     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}]"
+    scan \${xspec_tclout} "%f" cff_erg
+    puts \${cff_fd} "\${radius}   [expr \${cff_erg}/\${cf_photon}]"
 _EOF_
 if  [ ! -z "${COOLFUNC_BOLO}" ]; then
     cat >> ${XSPEC_CF_XCM} << _EOF_
     # coolfunc bolometric
-    set cfbolo_data \$cff_data
-    #puts "cfbolo_data: \${cfbolo_data}"
-    puts \${cfbolo_fd} "\${radius}    \${cfbolo_data}"
+    puts \${cfbolo_fd} "\${radius}    \${cff_erg}"
 _EOF_
 fi
 cat >> ${XSPEC_CF_XCM} << _EOF_
diff --git a/mass_profile/coolfunc_calc_erg.sh b/mass_profile/coolfunc_calc_erg.sh
index cd44160..02db453 100755
--- a/mass_profile/coolfunc_calc_erg.sh
+++ b/mass_profile/coolfunc_calc_erg.sh
@@ -16,7 +16,7 @@ if [ $# -ne 6 ]; then
     exit 1
 fi
 TPROFILE=$1
-ABUND_VAL=$2
+ABUNDANCE=$2
 N_H=$3
 REDSHIFT=$4
 COOLFUNC_PREFIX=$5
@@ -49,7 +49,7 @@ set xs_echo_script 0
 ## set basic data {{{
 set nh ${N_H}
 set redshift ${REDSHIFT}
-set abund_val ${ABUND_VAL}
+set abundance ${ABUNDANCE}
 set norm ${NORM}
 ## basic }}}
 
@@ -62,7 +62,7 @@ abund grsa
 dummyrsp 0.01 100.0 4096 linear
 # load model 'wabs*apec' to calc cooling function
 # (nh=0.0: do not consider aborption ???)
-model wabs*apec & 0.0 & 1.0 & \${abund_val} & \${redshift} & \${norm} &
+model wabs*apec & 0.0 & 1.0 & \${abundance} & \${redshift} & \${norm} &
 ## xspec }}}
 
 ## set input and output filename
@@ -91,17 +91,16 @@ while { [ gets \${blist_fd} blist_line ] != -1 } {
 
     ## 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}"
+        scan \${tpro_line} "%f %f" radius temperature
+        #puts "radius: \${radius}, temperature: \${temperature}"
         # set temperature value
-        newpar 2 \${temp_val}
+        newpar 2 \${temperature}
         # calc flux & tclout
         flux \${e1} \${e2}
         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}"
+        scan \${xspec_tclout} "%f" cf_erg
+        #puts "cf: \${cf_erg}"
+        puts \${cf_fd} "\${radius}    \${cf_erg}"
     }
     close \${tpro_fd}
     close \${cf_fd}
-- 
cgit v1.2.2