amoy-powerspectrum.sh

#!/bin/bash
#-----------------------------------------------
#-----------------------------------------------
### Define user-paths for convenience

export DATAFOLDER=/home/amoy/xnav/nicer-data/B1937+21
export CODESFOLDER=/home/amoy/xnav/codes

#-----------------------------------------------
#-----------------------------------------------
### Loop over the ObsIds in the data set,
### and merge the various files required for
### making the pulse profile

mapfile -t < $CODESFOLDER/obsids.txt # get list of ObsIDs as array


for i in "${MAPFILE[@]}"
do
	cd $DATAFOLDER
    export CWD=$DATAFOLDER/$i # Current Working Directory
    export CURAUX=$CWD/auxil # Current auxil Directory
    export CURXTI=$CWD/xti # Current xti Directory
 
	echo
    echo "--------------------------------------------------------"
    echo "*******Looping ... Currently processing ObsID $i*******"
    echo "--------------------------------------------------------"	
	echo
	
	nicerl3-lc $i pirange=300-1500 timebin=60.0 bkgmodeltype=scorpeon clobber=YES
	
	LOAD_LC=$(find "$CURXTI/event_cl" -type f -name "*_sr.lc")
	LC_FILE=$(find "$CURXTI/event_cl" -type f -name "*_sr.lc" -exec basename {} \; | tr '\n' ' ')
	if [[ -f "$LOAD_LC" ]]; then
		cd $CURXTI/event_cl
		powspec <<EOF
$LC_FILE
-
INDEF
INDEF
INDEF
0
default
yes
/XW
color 2 on 2
cpd $i.gif/gif
p
quit
EOF
	else
		echo "File *_sr.lc not found for ObsID $i"
    fi
done