#! /bin/tcsh -f

# t_lcvs 4.0

# Creates light curves using Xronos
# Uses code from do_lcvs (v3.2) script by D.Turcan

# Adaptation by Paul O'Neill

# Calls the scripts: t_lcvs_nbint
# Uses the script  : do_html_flc.xco
# Uses the windows : win_full2.wi win_full1.wi win_10per.wi

# NOTE: the script t_lcvs_nbint is called.  This script contains the 
# following lines:
#
#   #! /bin/tcsh -f
#
#   set diff = `echo $1 $2 | awk '{print int(($2-$1))}'`
#   echo $diff $3 | awk '{print int($1/$2+1)}' 
#

# Date in final form: 17 December 2004


########################################################################################
# Welcome 
########################################################################################

echo "gnt_infrm (t_lcvs)"
echo "gnt_infrm (t_lcvs) ++++++++++++++++++++++++++++++++++++++++++++++++++++++"
echo "gnt_infrm (t_lcvs)  t_lcvs 4.0  STARTED "`date`
echo "gnt_infrm (t_lcvs) ++++++++++++++++++++++++++++++++++++++++++++++++++++++"
echo "gnt_infrm (t_lcvs)"


########################################################################################
# Get parameters
########################################################################################

set obsid = $1           # observation sequence number
set SEQ = $2             # root directory for current sequence
set SCRIPTS = $3         # directory containing scripts
set XCO = $4             # directory containing xco scripts


########################################################################################
# Set directory variables
########################################################################################

# NOTE: aux, screened, unscreened, and work directories MUST be sub-directories
# of the sequence directory (SEQ). This structure is hard wired.

set AUX = $SEQ/aux            # contains house keeping data
set SCR = $SEQ/screened       # contains screened data
set UNSCR = $SEQ/unscreened   # contains unscreened data
set WORK = $SEQ/work          # a work area


########################################################################################
# Check for necessary files
########################################################################################

cd $WORK


#### look for GIS event files ####

foreach file ( gis2_src.evt gis3_src.evt )
    if ( ! -e $WORK/${obsid}_${file} ) then
	echo "gnt_fatal (t_lcvs) File $file does not exist."
	exit 1
    endif
end


#### look for SIS event files ####

foreach inst ( sis0 sis1 )
    foreach ext ( src src_soft src_soft2 src_hard bgd bgd_soft bgd_soft2 bgd_hard )
	set file = ${obsid}_${inst}_${ext}.evt
	if ( ! -e $WORK/$file ) then
	    echo "gnt_fatal (t_lcvs) File $file does not exist."
	    exit 1
	endif
    end
end


########################################################################################
# Check time order
########################################################################################

echo "gnt_infrm (t_lcvs) Time ordering the event files"


#### get start and stop times ####

set s0s = `fkeyprint ${obsid}_sis0_src.evt+1 TSTART | grep "TSTART  =" | awk '{print $3}'`
set s0t = `fkeyprint ${obsid}_sis0_src.evt+1 TSTOP | grep "TSTOP   =" | awk '{print $3}'`
set s1s = `fkeyprint ${obsid}_sis1_src.evt+1 TSTART | grep "TSTART  =" | awk '{print $3}'`
set s1t = `fkeyprint ${obsid}_sis1_src.evt+1 TSTOP | grep "TSTOP   =" | awk '{print $3}'`
set g2s = `fkeyprint ${obsid}_gis2_src.evt+1 TSTART | grep "TSTART  =" | awk '{print $3}'`
set g2t = `fkeyprint ${obsid}_gis2_src.evt+1 TSTOP | grep "TSTOP   =" | awk '{print $3}'`
set g3s = `fkeyprint ${obsid}_gis3_src.evt+1 TSTART | grep "TSTART  =" | awk '{print $3}'`
set g3t = `fkeyprint ${obsid}_gis3_src.evt+1 TSTOP | grep "TSTOP   =" | awk '{print $3}'`


#### determine SIS0 and SIS1 relative START times ####

set diff = `echo $s0s $s1s | awk '{print ($1 > $2)}'`  # = 1 if SIS0 start > SIS1 start 
                                                       # ie SIS1 start first
if ( $diff == "1" ) then             # SIS1 started observing before SIS0
    set sis = `echo sis1 sis0`
    set sbeg = $s1s    
else                              # SIS0 started before or same as SIS1
    set sis = `echo sis0 sis1`
    set sbeg = $s0s
endif


#### determine SIS0 and SIS1 relative STOP times ####

set diff = `echo $s1t $s0t | awk '{print ($1 > $2)}'`  # = 1 if SIS1 stop > SIS0 stop 
                                                       # ie SIS1 stopped last

if ( $diff == "1" ) then         #  SIS1 stoppped later than SIS0
    set send = $s1t
else                             #  SIS0 stopped after or same as SIS1
    set send = $s0t
endif


#### determine GIS2 and GIS3 relative START times ####

set diff = `echo $g2s $g3s | awk '{print ($1 > $2)}'`  # = 1 if GIS2 start > GIS3 start 
                                                       # ie GIS3 start first

if ( $diff == "1" ) then             # GIS3 started observing before GIS2
    set gis = `echo gis3 gis2`
    set gbeg = $g3s    
else                                 # GIS2 started before or same as GIS3
    set gis = `echo gis2 gis3`
    set gbeg = $g2s
endif


#### determine GIS3 and GIS2 relative STOP times ####

set diff = `echo $g3t $g2t | awk '{print ($1 > $2)}'`  # = 1 if GIS3 stop > GIS2 stop 
                                                       # ie GIS3 stopped last

if ( $diff == "1" ) then         #  GIS3 stoppped later than GIS2
    set gend = $g3t
else                             #  SIS0 stopped after or same as SIS1
    set gend = $g2t
endif


########################################################################################
# Create light curve  -  source, full-band, 16s bin, SIS
########################################################################################

echo "gnt_infrm (t_lcvs)"
echo "gnt_infrm (t_lcvs) Source, full-band, 16s bin, SIS"

echo $WORK/${obsid}_$sis[1]_src.evt >! $WORK/lcvs.tmp
echo $WORK/${obsid}_$sis[2]_src.evt >> $WORK/lcvs.tmp

set nbint = `$SCRIPTS/t_lcvs_nbint $sbeg $send 16`     # number of 16s bins in time-series
echo "nbint = $nbint"
if ( $nbint < 2 ) echo "gnt_warng (t_lcvs) Interval ($nbint) is smaller than requested bin size (16)."


#### create light curve using Xronos lcurve task ####

echo "gnt_infrm (t_lcvs) Xronos command being used is:"

echo "gnt_infrm (t_lcvs) lcurve nser=1 cfile1=@$WORK/lcvs.tmp"
echo "gnt_infrm (t_lcvs) window=$SCRIPTS/win_full2.wi dtnb=16 nbint=$nbint tunits=1 clobber=yes"
echo "gnt_infrm (t_lcvs) outfile=$WORK/${obsid}_sis_src_16s.flc rescale=2.0 plot=no"

lcurve nser=1 cfile1=@$WORK/lcvs.tmp \
    window=$SCRIPTS/win_full2.wi dtnb=16 nbint=$nbint tunits=1 clobber=yes \
    outfile=$WORK/${obsid}_sis_src_16s.flc rescale=2.0 plot=no

fselect $WORK/${obsid}_sis_src_16s.flc+1 $WORK/${obsid}_sis_src_16s.flc \
    expr='\!ISNULL(RATE1)' clobber=yes     # remove null bins


########################################################################################
# Create light curve  -  source, full-band, 16s bin, GIS
########################################################################################

echo "gnt_infrm (t_lcvs)"
echo "gnt_infrm (t_lcvs) Source, full-band, 16s bin, GIS"

echo $WORK/${obsid}_$gis[1]_src.evt >! $WORK/lcvs.tmp
echo $WORK/${obsid}_$gis[2]_src.evt >> $WORK/lcvs.tmp

set nbint = `$SCRIPTS/t_lcvs_nbint $gbeg $gend 16`     # number of 16s bins in time-series
echo "nbint = $nbint"
if ( $nbint < 2 ) echo "gnt_warng (t_lcvs): Interval ($nbint) is smaller than requested bin size (16)."


#### create light curve using Xronos lcurve task ####

echo "gnt_infrm (t_lcvs) Xronos command being used is:"

echo "gnt_infrm (t_lcvs) lcurve nser=1 cfile1=@$WORK/lcvs.tmp "
echo "gnt_infrm (t_lcvs) window=$SCRIPTS/win_full2.wi dtnb=16 nbint=$nbint tunits=1 clobber=yes "
echo "gnt_infrm (t_lcvs) outfile=$WORK/${obsid}_gis_src_16s.flc rescale=2.0 plot=no"

lcurve nser=1 cfile1=@$WORK/lcvs.tmp \
    window=$SCRIPTS/win_full2.wi dtnb=16 nbint=$nbint tunits=1 clobber=yes \
    outfile=$WORK/${obsid}_gis_src_16s.flc rescale=2.0 plot=no

fselect $WORK/${obsid}_gis_src_16s.flc+1 $WORK/${obsid}_gis_src_16s.flc \
    expr='\!ISNULL(RATE1)' clobber=yes      # remove null bins


########################################################################################
# Create light curve  -  source, full-band, 256s bin, SIS
########################################################################################

echo "gnt_infrm (t_lcvs)"
echo "gnt_infrm (t_lcvs) Source, full-band, 256s bin, SIS"

echo $WORK/${obsid}_sis_src_16s.flc >! $WORK/lcvs.tmp  # use the 16s light curve

set nbint = `$SCRIPTS/t_lcvs_nbint $sbeg $send 256`     # number of 256s bins in time-series
echo "nbint = $nbint"
if ( $nbint < 2 ) echo "gnt_warng (t_lcvs) Interval ($nbint) is smaller than requested bin size (256)."


#### create light curve using Xronos lcurve task ####

echo "gnt_infrm (t_lcvs) Xronos command being used is:"

echo "gnt_infrm (t_lcvs) lcurve nser=1 cfile1=@$WORK/lcvs.tmp"
echo "gnt_infrm (t_lcvs) window=$SCRIPTS/win_full1.wi dtnb=256 nbint=$nbint tunits=1 clobber=yes"
echo "gnt_infrm (t_lcvs) outfile=$WORK/${obsid}_sis_src_256s.flc plot=no rescale=1.0"

lcurve nser=1 cfile1=@$WORK/lcvs.tmp \
    window=$SCRIPTS/win_full1.wi dtnb=256 nbint=$nbint tunits=1 clobber=yes \
    outfile=$WORK/${obsid}_sis_src_256s.flc plot=no rescale=1.0

fselect $WORK/${obsid}_sis_src_256s.flc+1 $WORK/${obsid}_sis_src_256s.flc \
    expr='\!ISNULL(RATE1)' clobber=yes         # remove null bins


########################################################################################
# Create light curve  -  source, full-band, 256s bin, GIS
########################################################################################

echo "gnt_infrm (t_lcvs)"
echo "gnt_infrm (t_lcvs) Source, full-band, 256s bin, GIS"

echo $WORK/${obsid}_gis_src_16s.flc >! $WORK/lcvs.tmp   # use the 16s light curve

set nbint = `$SCRIPTS/t_lcvs_nbint $gbeg $gend 256`    # number of  256s bins in time-series
echo nbint = $nbint
if ( $nbint < 2 ) echo "gnt_warng (t_lcvs) Interval ($nbint) is smaller than requested bin size (256)."

echo "gnt_infrm (t_lcvs) Xronos command being used is:"

echo "gnt_infrm (t_lcvs) lcurve nser=1 cfile1=@$WORK/lcvs.tmp"
echo "gnt_infrm (t_lcvs) window=$SCRIPTS/win_full1.wi dtnb=256 nbint=$nbint tunits=1 clobber=yes"
echo "gnt_infrm (t_lcvs) outfile=$WORK/${obsid}_gis_src_256s.flc plot=no rescale=1.0"

lcurve nser=1 cfile1=@$WORK/lcvs.tmp \
    window=$SCRIPTS/win_full1.wi dtnb=256 nbint=$nbint tunits=1 clobber=yes \
    outfile=$WORK/${obsid}_gis_src_256s.flc plot=no rescale=1.0

fselect $WORK/${obsid}_gis_src_256s.flc+1 $WORK/${obsid}_gis_src_256s.flc \
   expr='\!ISNULL(RATE1)' clobber=yes         # remove null bins


########################################################################################
# Create light curve  -  source, full-band, 5760s bin, SIS
########################################################################################

echo "gnt_infrm (t_lcvs)"
echo "gnt_infrm (t_lcvs) Source, full-band, 5760s bin, SIS"

echo $WORK/${obsid}_sis_src_16s.flc >! $WORK/lcvs.tmp

set nbint = `$SCRIPTS/t_lcvs_nbint $sbeg $send 5760`
echo nbint = $nbint
if ( $nbint < 2 ) echo "gnt_warng (t_lcvs) Interval ($nbint) is smaller than requested bin size (5760)."

echo "gnt_infrm (t_lcvs) Xronos command being used is:"

echo "gnt_infrm (t_lcvs) lcurve nser=1 cfile1=@$WORK/lcvs.tmp"
echo "gnt_infrm (t_lcvs) window=$SCRIPTS/win_10per.wi dtnb=5760 nbint=$nbint tunits=1 clobber=yes"
echo "gnt_infrm (t_lcvs) outfile=$WORK/${obsid}_sis_src_5760s.flc plot=no rescale=1.0"

lcurve nser=1 cfile1=@$WORK/lcvs.tmp \
    window=$SCRIPTS/win_10per.wi dtnb=5760 nbint=$nbint tunits=1 clobber=yes \
    outfile=$WORK/${obsid}_sis_src_5760s.flc plot=no rescale=1.0

fselect $WORK/${obsid}_sis_src_5760s.flc+1 $WORK/${obsid}_sis_src_5760s.flc \
   expr='\!ISNULL(RATE1)' clobber=yes


########################################################################################
# Create light curve  -  source, full-band, 5760s bin, GIS
########################################################################################

echo "gnt_infrm (t_lcvs)"
echo "gnt_infrm (t_lcvs) Source, full-band, 5760s bin, GIS"

echo $WORK/${obsid}_gis_src_16s.flc >! $WORK/lcvs.tmp

set nbint = `$SCRIPTS/t_lcvs_nbint $gbeg $gend 5760`
echo nbint = $nbint
if ( $nbint < 2 ) echo "gnt_warng (t_lcvs): Interval ($nbint) is smaller than requested bin size (5760)."

echo "gnt_infrm (t_lcvs) Xronos command being used is:"

echo "gnt_infrm (t_lcvs) lcurve nser=1 cfile1=@$WORK/lcvs.tmp"
echo "gnt_infrm (t_lcvs) window=$SCRIPTS/win_10per.wi dtnb=5760 nbint=$nbint tunits=1 clobber=yes"
echo "gnt_infrm (t_lcvs) outfile=$WORK/${obsid}_gis_src_5760s.flc plot=no rescale=1.0"

lcurve nser=1 cfile1=@$WORK/lcvs.tmp \
    window=$SCRIPTS/win_10per.wi dtnb=5760 nbint=$nbint tunits=1 clobber=yes \
    outfile=$WORK/${obsid}_gis_src_5760s.flc plot=no rescale=1.0

fselect $WORK/${obsid}_gis_src_5760s.flc+1 $WORK/${obsid}_gis_src_5760s.flc \
    expr='\!ISNULL(RATE1)' clobber=yes


########################################################################################
# Create light curve  -  source, full2-band, 16s bin, SIS
########################################################################################

echo "gnt_infrm (t_lcvs)"
echo "gnt_infrm (t_lcvs) Source, full2-band, 16s bin, SIS"

echo $WORK/${obsid}_$sis[1]_src_full2.evt >! $WORK/lcvs.tmp
echo $WORK/${obsid}_$sis[2]_src_full2.evt >> $WORK/lcvs.tmp

set nbint = `$SCRIPTS/t_lcvs_nbint $sbeg $send 16`
echo nbint = $nbint
if ( $nbint < 2 ) echo "gnt_warng (t_lcvs): Interval ($nbint) is smaller than requested bin size (16)."

echo "gnt_infrm (t_lcvs) Xronos command being used is:"

echo "gnt_infrm (t_lcvs) lcurve nser=1 cfile1=@$WORK/lcvs.tmp"
echo "gnt_infrm (t_lcvs) window=$SCRIPTS/win_full2.wi dtnb=16 nbint=$nbint tunits=1 clobber=yes"
echo "gnt_infrm (t_lcvs) outfile=$WORK/${obsid}_sis_src_16s_full2.flc rescale=2.0 plot=no "

lcurve nser=1 cfile1=@$WORK/lcvs.tmp \
    window=$SCRIPTS/win_full2.wi dtnb=16 nbint=$nbint tunits=1 clobber=yes \
    outfile=$WORK/${obsid}_sis_src_16s_full2.flc rescale=2.0 plot=no

fselect $WORK/${obsid}_sis_src_16s_full2.flc+1 $WORK/${obsid}_sis_src_16s_full2.flc \
   expr='\!ISNULL(RATE1)' clobber=yes


########################################################################################
# Create light curve  -  source, full2-band, 256s bin, SIS
########################################################################################

echo "gnt_infrm (t_lcvs)"
echo "gnt_infrm (t_lcvs) Source, full2-band, 256s bin, SIS"

echo $WORK/${obsid}_sis_src_16s_full2.flc >! $WORK/lcvs.tmp  # use the 16s light curve

set nbint = `$SCRIPTS/t_lcvs_nbint $sbeg $send 256`     # number of 256s bins in time-series
echo "nbint = $nbint"
if ( $nbint < 2 ) echo "gnt_warng (t_lcvs) Interval ($nbint) is smaller than requested bin size (256)."


#### create light curve using Xronos lcurve task ####

echo "gnt_infrm (t_lcvs) Xronos command being used is:"

echo "gnt_infrm (t_lcvs) lcurve nser=1 cfile1=@$WORK/lcvs.tmp"
echo "gnt_infrm (t_lcvs) window=$SCRIPTS/win_full1.wi dtnb=256 nbint=$nbint tunits=1 clobber=yes"
echo "gnt_infrm (t_lcvs) outfile=$WORK/${obsid}_sis_src_256s_full2.flc plot=no rescale=1.0"

lcurve nser=1 cfile1=@$WORK/lcvs.tmp \
    window=$SCRIPTS/win_full1.wi dtnb=256 nbint=$nbint tunits=1 clobber=yes \
    outfile=$WORK/${obsid}_sis_src_256s_full2.flc plot=no rescale=1.0

fselect $WORK/${obsid}_sis_src_256s_full2.flc+1 $WORK/${obsid}_sis_src_256s_full2.flc \
    expr='\!ISNULL(RATE1)' clobber=yes         # remove null bins



########################################################################################
# Create light curve  -  source, full2-band, 5760s bin, SIS
########################################################################################

echo "gnt_infrm (t_lcvs)"
echo "gnt_infrm (t_lcvs) Source, full2-band, 5760s bin, SIS"

echo $WORK/${obsid}_sis_src_16s_full2.flc >! $WORK/lcvs.tmp

set nbint = `$SCRIPTS/t_lcvs_nbint $sbeg $send 5760`
echo nbint = $nbint
if ( $nbint < 2 ) echo "gnt_warng (t_lcvs) Interval ($nbint) is smaller than requested bin size (5760)."

echo "gnt_infrm (t_lcvs) Xronos command being used is:"

echo "gnt_infrm (t_lcvs) lcurve nser=1 cfile1=@$WORK/lcvs.tmp"
echo "gnt_infrm (t_lcvs) window=$SCRIPTS/win_10per.wi dtnb=5760 nbint=$nbint tunits=1 clobber=yes"
echo "gnt_infrm (t_lcvs) outfile=$WORK/${obsid}_sis_src_5760s_full2.flc plot=no rescale=1.0"

lcurve nser=1 cfile1=@$WORK/lcvs.tmp \
    window=$SCRIPTS/win_10per.wi dtnb=5760 nbint=$nbint tunits=1 clobber=yes \
    outfile=$WORK/${obsid}_sis_src_5760s_full2.flc plot=no rescale=1.0

fselect $WORK/${obsid}_sis_src_5760s_full2.flc+1 $WORK/${obsid}_sis_src_5760s_full2.flc \
   expr='\!ISNULL(RATE1)' clobber=yes


########################################################################################
# Create light curve  -  source, soft-band, 16s bin, SIS
########################################################################################

echo "gnt_infrm (t_lcvs)"
echo "gnt_infrm (t_lcvs) Source, soft-band, 16s bin, SIS"

echo $WORK/${obsid}_$sis[1]_src_soft.evt >! $WORK/lcvs.tmp
echo $WORK/${obsid}_$sis[2]_src_soft.evt >> $WORK/lcvs.tmp

set nbint = `$SCRIPTS/t_lcvs_nbint $sbeg $send 16`
echo nbint = $nbint
if ( $nbint < 2 ) echo "gnt_warng (t_lcvs): Interval ($nbint) is smaller than requested bin size (16)."

echo "gnt_infrm (t_lcvs) Xronos command being used is:"

echo "gnt_infrm (t_lcvs) lcurve nser=1 cfile1=@$WORK/lcvs.tmp"
echo "gnt_infrm (t_lcvs) window=$SCRIPTS/win_full2.wi dtnb=16 nbint=$nbint tunits=1 clobber=yes"
echo "gnt_infrm (t_lcvs) outfile=$WORK/${obsid}_sis_src_16s_soft.flc rescale=2.0 plot=no "

lcurve nser=1 cfile1=@$WORK/lcvs.tmp \
    window=$SCRIPTS/win_full2.wi dtnb=16 nbint=$nbint tunits=1 clobber=yes \
    outfile=$WORK/${obsid}_sis_src_16s_soft.flc rescale=2.0 plot=no

fselect $WORK/${obsid}_sis_src_16s_soft.flc+1 $WORK/${obsid}_sis_src_16s_soft.flc \
   expr='\!ISNULL(RATE1)' clobber=yes


########################################################################################
# Create light curve  -  source, soft-band, 5760s bin, SIS
########################################################################################

echo "gnt_infrm (t_lcvs)"
echo "gnt_infrm (t_lcvs) Source, soft-band, 5760s bin, SIS"

echo $WORK/${obsid}_sis_src_16s_soft.flc >! $WORK/lcvs.tmp

set nbint = `$SCRIPTS/t_lcvs_nbint $sbeg $send 5760`
echo nbint = $nbint
if ( $nbint < 2 ) echo "gnt_warng (t_lcvs): Interval ($nbint) is smaller than requested bin size (5760)."

echo "gnt_infrm (t_lcvs) Xronos command being used is:"

echo "gnt_infrm (t_lcvs) lcurve nser=1 cfile1=@$WORK/lcvs.tmp"
echo "gnt_infrm (t_lcvs) window=$SCRIPTS/win_10per.wi dtnb=5760 nbint=$nbint tunits=1 clobber=yes"
echo "gnt_infrm (t_lcvs) outfile=$WORK/${obsid}_sis_src_5760s_soft.flc plot=no rescale=1.0"

lcurve nser=1 cfile1=@$WORK/lcvs.tmp \
    window=$SCRIPTS/win_10per.wi dtnb=5760 nbint=$nbint tunits=1 clobber=yes \
    outfile=$WORK/${obsid}_sis_src_5760s_soft.flc plot=no rescale=1.0

fselect $WORK/${obsid}_sis_src_5760s_soft.flc+1 $WORK/${obsid}_sis_src_5760s_soft.flc \
    expr='\!ISNULL(RATE1)' clobber=yes


########################################################################################
# Create light curve  -  source, soft2-band, 16s bin, SIS
########################################################################################

echo "gnt_infrm (t_lcvs)"
echo "gnt_infrm (t_lcvs) Source, soft2-band, 16s bin, SIS"

echo $WORK/${obsid}_$sis[1]_src_soft2.evt >! $WORK/lcvs.tmp
echo $WORK/${obsid}_$sis[2]_src_soft2.evt >> $WORK/lcvs.tmp

set nbint = `$SCRIPTS/t_lcvs_nbint $sbeg $send 16`
echo nbint = $nbint
if ( $nbint < 2 ) echo "gnt_warng (t_lcvs): Interval ($nbint) is smaller than requested bin size (16)."

echo "gnt_infrm (t_lcvs) Xronos command being used is:"

echo "gnt_infrm (t_lcvs) lcurve nser=1 cfile1=@$WORK/lcvs.tmp"
echo "gnt_infrm (t_lcvs) window=$SCRIPTS/win_full2.wi dtnb=16 nbint=$nbint tunits=1 clobber=yes"
echo "gnt_infrm (t_lcvs) outfile=$WORK/${obsid}_sis_src_16s_soft2.flc rescale=2.0 plot=no "

lcurve nser=1 cfile1=@$WORK/lcvs.tmp \
    window=$SCRIPTS/win_full2.wi dtnb=16 nbint=$nbint tunits=1 clobber=yes \
    outfile=$WORK/${obsid}_sis_src_16s_soft2.flc rescale=2.0 plot=no

fselect $WORK/${obsid}_sis_src_16s_soft2.flc+1 $WORK/${obsid}_sis_src_16s_soft2.flc \
   expr='\!ISNULL(RATE1)' clobber=yes


########################################################################################
# Create light curve  -  source, soft2-band, 5760s bin, SIS
########################################################################################

echo "gnt_infrm (t_lcvs)"
echo "gnt_infrm (t_lcvs) Source, soft2-band, 5760s bin, SIS"

echo $WORK/${obsid}_sis_src_16s_soft2.flc >! $WORK/lcvs.tmp

set nbint = `$SCRIPTS/t_lcvs_nbint $sbeg $send 5760`
echo nbint = $nbint
if ( $nbint < 2 ) echo "gnt_warng (t_lcvs): Interval ($nbint) is smaller than requested bin size (5760)."

echo "gnt_infrm (t_lcvs) Xronos command being used is:"

echo "gnt_infrm (t_lcvs) lcurve nser=1 cfile1=@$WORK/lcvs.tmp"
echo "gnt_infrm (t_lcvs) window=$SCRIPTS/win_10per.wi dtnb=5760 nbint=$nbint tunits=1 clobber=yes"
echo "gnt_infrm (t_lcvs) outfile=$WORK/${obsid}_sis_src_5760s_soft2.flc plot=no rescale=1.0"

lcurve nser=1 cfile1=@$WORK/lcvs.tmp \
    window=$SCRIPTS/win_10per.wi dtnb=5760 nbint=$nbint tunits=1 clobber=yes \
    outfile=$WORK/${obsid}_sis_src_5760s_soft2.flc plot=no rescale=1.0

fselect $WORK/${obsid}_sis_src_5760s_soft2.flc+1 $WORK/${obsid}_sis_src_5760s_soft2.flc \
    expr='\!ISNULL(RATE1)' clobber=yes


########################################################################################
# Create light curve  -  source, hard-band, 16s bin, SIS
########################################################################################

echo "gnt_infrm (t_lcvs)"
echo "gnt_infrm (t_lcvs) Source, hard-band, 16s bin, SIS"

echo $WORK/${obsid}_$sis[1]_src_hard.evt >! $WORK/lcvs.tmp
echo $WORK/${obsid}_$sis[2]_src_hard.evt >> $WORK/lcvs.tmp

set nbint = `$SCRIPTS/t_lcvs_nbint $sbeg $send 16`
echo nbint = $nbint
if ( $nbint < 2 ) echo "gnt_warng (t_lcvs): Interval ($nbint) is smaller than requested bin size (16)."

echo "gnt_infrm (t_lcvs) Xronos command being used is:"

echo "gnt_infrm (t_lcvs) lcurve nser=1 cfile1=@$WORK/lcvs.tmp"
echo "gnt_infrm (t_lcvs) window=$SCRIPTS/win_full2.wi dtnb=16 nbint=$nbint tunits=1 clobber=yes"
echo "gnt_infrm (t_lcvs) outfile=$WORK/${obsid}_sis_src_16s_hard.flc rescale=2.0 plot=no"

lcurve nser=1 cfile1=@$WORK/lcvs.tmp \
    window=$SCRIPTS/win_full2.wi dtnb=16 nbint=$nbint tunits=1 clobber=yes \
    outfile=$WORK/${obsid}_sis_src_16s_hard.flc rescale=2.0 plot=no

fselect $WORK/${obsid}_sis_src_16s_hard.flc+1 $WORK/${obsid}_sis_src_16s_hard.flc \
    expr='\!ISNULL(RATE1)' clobber=yes


########################################################################################
# Create light curve  -  source, hard-band, 5760s bin, SIS
########################################################################################

echo "gnt_infrm (t_lcvs)"
echo "gnt_infrm (t_lcvs) Source, hard-band, 5760s bin, SIS"

echo $WORK/${obsid}_sis_src_16s_hard.flc >! $WORK/lcvs.tmp

set nbint = `$SCRIPTS/t_lcvs_nbint $sbeg $send 5760`
echo nbint = $nbint
if ( $nbint < 2 ) echo "gnt_warng (t_lcvs): Interval ($nbint) is smaller than requested bin size (5760)."

echo "gnt_infrm (t_lcvs) Xronos command being used is:"

echo "gnt_infrm (t_lcvs) lcurve nser=1 cfile1=@$WORK/lcvs.tmp"
echo "gnt_infrm (t_lcvs) window=$SCRIPTS/win_10per.wi dtnb=5760 nbint=$nbint tunits=1 clobber=yes"
echo "gnt_infrm (t_lcvs) outfile=$WORK/${obsid}_sis_src_5760s_hard.flc plot=no rescale=1.0"

lcurve nser=1 cfile1=@$WORK/lcvs.tmp \
    window=$SCRIPTS/win_10per.wi dtnb=5760 nbint=$nbint tunits=1 clobber=yes \
    outfile=$WORK/${obsid}_sis_src_5760s_hard.flc plot=no rescale=1.0

fselect $WORK/${obsid}_sis_src_5760s_hard.flc+1 $WORK/${obsid}_sis_src_5760s_hard.flc \
    expr='\!ISNULL(RATE1)' clobber=yes


########################################################################################
# Create light curve  -  source, hard-band, 16s bin, GIS
########################################################################################

echo "gnt_infrm (t_lcvs)"
echo "gnt_infrm (t_lcvs) Source, hard-band, 16s bin, GIS"

echo $WORK/${obsid}_$gis[1]_src_hard.evt >! $WORK/lcvs.tmp
echo $WORK/${obsid}_$gis[2]_src_hard.evt >> $WORK/lcvs.tmp

set nbint = `$SCRIPTS/t_lcvs_nbint $gbeg $gend 16`
echo nbint = $nbint
if ( $nbint < 2 ) echo "gnt_warng (t_lcvs): Interval ($nbint) is smaller than requested bin size (16)."

echo "gnt_infrm (t_lcvs) Xronos command being used is:"

echo "gnt_infrm (t_lcvs) lcurve nser=1 cfile1=@$WORK/lcvs.tmp"
echo "gnt_infrm (t_lcvs) window=$SCRIPTS/win_full2.wi dtnb=16 nbint=$nbint tunits=1 clobber=yes"
echo "gnt_infrm (t_lcvs) outfile=$WORK/${obsid}_gis_src_16s_hard.flc rescale=2.0 plot=no"

lcurve nser=1 cfile1=@$WORK/lcvs.tmp \
    window=$SCRIPTS/win_full2.wi dtnb=16 nbint=$nbint tunits=1 clobber=yes \
    outfile=$WORK/${obsid}_gis_src_16s_hard.flc rescale=2.0 plot=no

fselect $WORK/${obsid}_gis_src_16s_hard.flc+1 $WORK/${obsid}_gis_src_16s_hard.flc \
    expr='\!ISNULL(RATE1)' clobber=yes


########################################################################################
# Create light curve  -  source, hard-band, 5760s bin, GIS
########################################################################################

echo "gnt_infrm (t_lcvs)"
echo "gnt_infrm (t_lcvs) Source, hard-band, 16s bin, GIS"

echo $WORK/${obsid}_gis_src_16s_hard.flc >! $WORK/lcvs.tmp

set nbint = `$SCRIPTS/t_lcvs_nbint $gbeg $gend 5760`
echo nbint = $nbint
if ( $nbint < 2 ) echo "gnt_warng (t_lcvs): Interval ($nbint) is smaller than requested bin size (5760)."

echo "gnt_infrm (t_lcvs) Xronos command being used is:"

echo "gnt_infrm (t_lcvs) lcurve nser=1 cfile1=@$WORK/lcvs.tmp"
echo "gnt_infrm (t_lcvs) window=$SCRIPTS/win_10per.wi dtnb=5760 nbint=$nbint tunits=1 clobber=yes"
echo "gnt_infrm (t_lcvs) outfile=$WORK/${obsid}_gis_src_5760s_hard.flc plot=no rescale=1.0"

lcurve nser=1 cfile1=@$WORK/lcvs.tmp \
    window=$SCRIPTS/win_10per.wi dtnb=5760 nbint=$nbint tunits=1 clobber=yes \
    outfile=$WORK/${obsid}_gis_src_5760s_hard.flc plot=no rescale=1.0

fselect $WORK/${obsid}_gis_src_5760s_hard.flc+1 $WORK/${obsid}_gis_src_5760s_hard.flc \
    expr='\!ISNULL(RATE1)' clobber=yes


########################################################################################
# Create light curve  -  background, full-band, 16s, SIS
########################################################################################

echo "gnt_infrm (t_lcvs)"
echo "gnt_infrm (t_lcvs) Background, full-band, 16s bin, SIS"

echo $WORK/${obsid}_$sis[1]_bgd.evt >! $WORK/lcvs.tmp
echo $WORK/${obsid}_$sis[2]_bgd.evt >> $WORK/lcvs.tmp

set nbint = `$SCRIPTS/t_lcvs_nbint $sbeg $send 16`
echo nbint = $nbint
if ( $nbint < 2 ) echo "gnt_warng (t_lcvs): Interval ($nbint) is smaller than requested bin size (16)."

echo "gnt_infrm (t_lcvs) Xronos command being used is:"

echo "gnt_infrm (t_lcvs) lcurve nser=1 cfile1=@$WORK/lcvs.tmp"
echo "gnt_infrm (t_lcvs) window=$SCRIPTS/win_full2.wi dtnb=16 nbint=$nbint tunits=1 clobber=yes"
echo "gnt_infrm (t_lcvs) outfile=$WORK/${obsid}_sis_bgd_16s.flc rescale=2.0 plot=no"

lcurve nser=1 cfile1=@$WORK/lcvs.tmp \
    window=$SCRIPTS/win_full2.wi dtnb=16 nbint=$nbint tunits=1 clobber=yes \
    outfile=$WORK/${obsid}_sis_bgd_16s.flc rescale=2.0 plot=no

fselect $WORK/${obsid}_sis_bgd_16s.flc+1 $WORK/${obsid}_sis_bgd_16s.flc \
    expr='\!ISNULL(RATE1)' clobber=yes


########################################################################################
# Create light curve  -  background, full-band, 5760s, SIS
########################################################################################

echo "gnt_infrm (t_lcvs)"
echo "gnt_infrm (t_lcvs) Background, full-band, 5760s bin, SIS"

echo $WORK/${obsid}_sis_bgd_16s.flc >! $WORK/lcvs.tmp

set nbint = `$SCRIPTS/t_lcvs_nbint $sbeg $send 5760`
echo nbint = $nbint
if ( $nbint < 2 ) echo "gnt_warng (t_lcvs): Interval ($nbint) is smaller than requested bin size (5760)."

echo "gnt_infrm (t_lcvs) Xronos command being used is:"

echo "gnt_infrm (t_lcvs) lcurve nser=1 cfile1=@$WORK/lcvs.tmp"
echo "gnt_infrm (t_lcvs) window=$SCRIPTS/win_10per.wi dtnb=5760 nbint=$nbint tunits=1 clobber=yes"
echo "gnt_infrm (t_lcvs) outfile=$WORK/${obsid}_sis_bgd_5760s.flc plot=no rescale=1.0"

lcurve nser=1 cfile1=@$WORK/lcvs.tmp \
    window=$SCRIPTS/win_10per.wi dtnb=5760 nbint=$nbint tunits=1 clobber=yes \
    outfile=$WORK/${obsid}_sis_bgd_5760s.flc plot=no rescale=1.0

fselect $WORK/${obsid}_sis_bgd_5760s.flc+1 $WORK/${obsid}_sis_bgd_5760s.flc \
    expr='\!ISNULL(RATE1)' clobber=yes


########################################################################################
# Create light curve  -  background, full-band, 16s, GIS
########################################################################################

echo "gnt_infrm (t_lcvs)"
echo "gnt_infrm (t_lcvs) Background, full-band, 16s bin, GIS"

echo $WORK/${obsid}_$gis[1]_bgd.evt >! $WORK/lcvs.tmp
echo $WORK/${obsid}_$gis[2]_bgd.evt >> $WORK/lcvs.tmp

set nbint = `$SCRIPTS/t_lcvs_nbint $gbeg $gend 16`
echo nbint = $nbint
if ( $nbint < 2 ) echo "gnt_warng (t_lcvs): Interval ($nbint) is smaller than requested bin size (16)."

echo "gnt_infrm (t_lcvs) Xronos command being used is:"

echo "gnt_infrm (t_lcvs) lcurve nser=1 cfile1=@$WORK/lcvs.tmp"
echo "gnt_infrm (t_lcvs) window=$SCRIPTS/win_full2.wi dtnb=16 nbint=$nbint tunits=1 clobber=yes"
echo "gnt_infrm (t_lcvs) outfile=$WORK/${obsid}_gis_bgd_16s.flc rescale=2.0 plot=no"

lcurve nser=1 cfile1=@$WORK/lcvs.tmp \
    window=$SCRIPTS/win_full2.wi dtnb=16 nbint=$nbint tunits=1 clobber=yes \
    outfile=$WORK/${obsid}_gis_bgd_16s.flc rescale=2.0 plot=no

fselect $WORK/${obsid}_gis_bgd_16s.flc+1 $WORK/${obsid}_gis_bgd_16s.flc \
    expr='\!ISNULL(RATE1)' clobber=yes


########################################################################################
# Create light curve  -  background, full-band, 5760s, GIS
########################################################################################

echo "gnt_infrm (t_lcvs)"
echo "gnt_infrm (t_lcvs) Background, full-band, 5760s bin, GIS"

echo $WORK/${obsid}_gis_bgd_16s.flc >! $WORK/lcvs.tmp

set nbint = `$SCRIPTS/t_lcvs_nbint $gbeg $gend 5760`
echo nbint = $nbint
if ( $nbint < 2 ) echo "gnt_warng (t_lcvs): Interval ($nbint) is smaller than requested bin size (5760)."

echo "gnt_infrm (t_lcvs) Xronos command being used is:"

echo "gnt_infrm (t_lcvs) lcurve nser=1 cfile1=@$WORK/lcvs.tmp"
echo "gnt_infrm (t_lcvs) window=$SCRIPTS/win_10per.wi dtnb=5760 nbint=$nbint tunits=1 clobber=yes"
echo "gnt_infrm (t_lcvs) outfile=$WORK/${obsid}_gis_bgd_5760s.flc plot=no rescale=1.0"

lcurve nser=1 cfile1=@$WORK/lcvs.tmp \
    window=$SCRIPTS/win_10per.wi dtnb=5760 nbint=$nbint tunits=1 clobber=yes \
    outfile=$WORK/${obsid}_gis_bgd_5760s.flc plot=no rescale=1.0

fselect $WORK/${obsid}_gis_bgd_5760s.flc+1 $WORK/${obsid}_gis_bgd_5760s.flc \
    expr='\!ISNULL(RATE1)' clobber=yes


########################################################################################
# Create light curve  -  background, full2-band, 16s, SIS
########################################################################################

echo "gnt_infrm (t_lcvs)"
echo "gnt_infrm (t_lcvs) Background, full2-band, 16s bin, SIS"

echo $WORK/${obsid}_$sis[1]_bgd_full2.evt >! $WORK/lcvs.tmp
echo $WORK/${obsid}_$sis[2]_bgd_full2.evt >> $WORK/lcvs.tmp

set nbint = `$SCRIPTS/t_lcvs_nbint $sbeg $send 16`
echo nbint = $nbint
if ( $nbint < 2 ) echo "gnt_warng (t_lcvs): Interval ($nbint) is smaller than requested bin size (16)."

echo "gnt_infrm (t_lcvs) Xronos command being used is:"

echo "gnt_infrm (t_lcvs) lcurve nser=1 cfile1=@$WORK/lcvs.tmp"
echo "gnt_infrm (t_lcvs) window=$SCRIPTS/win_full2.wi dtnb=16 nbint=$nbint tunits=1 clobber=yes"
echo "gnt_infrm (t_lcvs) outfile=$WORK/${obsid}_sis_bgd_16s_full2.flc rescale=2.0 plot=no"

lcurve nser=1 cfile1=@$WORK/lcvs.tmp \
    window=$SCRIPTS/win_full2.wi dtnb=16 nbint=$nbint tunits=1 clobber=yes \
    outfile=$WORK/${obsid}_sis_bgd_16s_full2.flc rescale=2.0 plot=no

fselect $WORK/${obsid}_sis_bgd_16s_full2.flc+1 $WORK/${obsid}_sis_bgd_16s_full2.flc \
    expr='\!ISNULL(RATE1)' clobber=yes


########################################################################################
# Create light curve  -  background, full2-band, 5760s, SIS
########################################################################################

echo "gnt_infrm (t_lcvs)"
echo "gnt_infrm (t_lcvs) Background, full2-band, 5760s bin, SIS"

echo $WORK/${obsid}_sis_bgd_16s_full2.flc >! $WORK/lcvs.tmp

set nbint = `$SCRIPTS/t_lcvs_nbint $sbeg $send 5760`
echo nbint = $nbint
if ( $nbint < 2 ) echo "gnt_warng (t_lcvs): Interval ($nbint) is smaller than requested bin size (5760)."

echo "gnt_infrm (t_lcvs) Xronos command being used is:"

echo "gnt_infrm (t_lcvs) lcurve nser=1 cfile1=@$WORK/lcvs.tmp"
echo "gnt_infrm (t_lcvs) window=$SCRIPTS/win_10per.wi dtnb=5760 nbint=$nbint tunits=1 clobber=yes"
echo "gnt_infrm (t_lcvs) outfile=$WORK/${obsid}_sis_bgd_5760s_full2.flc plot=no rescale=1.0"

lcurve nser=1 cfile1=@$WORK/lcvs.tmp \
    window=$SCRIPTS/win_10per.wi dtnb=5760 nbint=$nbint tunits=1 clobber=yes \
    outfile=$WORK/${obsid}_sis_bgd_5760s_full2.flc plot=no rescale=1.0

fselect $WORK/${obsid}_sis_bgd_5760s_full2.flc+1 $WORK/${obsid}_sis_bgd_5760s_full2.flc \
    expr='\!ISNULL(RATE1)' clobber=yes


########################################################################################
# Create light curve  -  background, soft-band, 16s, SIS
########################################################################################

echo "gnt_infrm (t_lcvs)"
echo "gnt_infrm (t_lcvs) Background, soft-band, 16s bin, SIS"

echo $WORK/${obsid}_$sis[1]_bgd_soft.evt >! $WORK/lcvs.tmp
echo $WORK/${obsid}_$sis[2]_bgd_soft.evt >> $WORK/lcvs.tmp

set nbint = `$SCRIPTS/t_lcvs_nbint $sbeg $send 16`
echo nbint = $nbint
if ( $nbint < 2 ) echo "gnt_warng (t_lcvs): Interval ($nbint) is smaller than requested bin size (16)."

echo "gnt_infrm (t_lcvs) Xronos command being used is:"

echo "gnt_infrm (t_lcvs) lcurve nser=1 cfile1=@$WORK/lcvs.tmp"
echo "gnt_infrm (t_lcvs) window=$SCRIPTS/win_full2.wi dtnb=16 nbint=$nbint"
echo "gnt_infrm (t_lcvs) tunits=1 clobber=yes outfile=$WORK/${obsid}_sis_bgd_16s_soft.flc rescale=2.0 plot=no"

lcurve nser=1 cfile1=@$WORK/lcvs.tmp \
    window=$SCRIPTS/win_full2.wi dtnb=16 nbint=$nbint \
    tunits=1 clobber=yes outfile=$WORK/${obsid}_sis_bgd_16s_soft.flc rescale=2.0 plot=no

fselect $WORK/${obsid}_sis_bgd_16s_soft.flc+1 $WORK/${obsid}_sis_bgd_16s_soft.flc \
    expr='\!ISNULL(RATE1)' clobber=yes


########################################################################################
# Create light curve  -  background, soft2-band, 16s, SIS
########################################################################################

echo "gnt_infrm (t_lcvs)"
echo "gnt_infrm (t_lcvs) Background, soft2-band, 16s bin, SIS"

echo $WORK/${obsid}_$sis[1]_bgd_soft2.evt >! $WORK/lcvs.tmp
echo $WORK/${obsid}_$sis[2]_bgd_soft2.evt >> $WORK/lcvs.tmp

set nbint = `$SCRIPTS/t_lcvs_nbint $sbeg $send 16`
echo nbint = $nbint
if ( $nbint < 2 ) echo "gnt_warng (t_lcvs): Interval ($nbint) is smaller than requested bin size (16)."

echo "gnt_infrm (t_lcvs) Xronos command being used is:"

echo "gnt_infrm (t_lcvs) lcurve nser=1 cfile1=@$WORK/lcvs.tmp"
echo "gnt_infrm (t_lcvs) window=$SCRIPTS/win_full2.wi dtnb=16 nbint=$nbint"
echo "gnt_infrm (t_lcvs) tunits=1 clobber=yes outfile=$WORK/${obsid}_sis_bgd_16s_soft2.flc rescale=2.0 plot=no"

lcurve nser=1 cfile1=@$WORK/lcvs.tmp \
    window=$SCRIPTS/win_full2.wi dtnb=16 nbint=$nbint \
    tunits=1 clobber=yes outfile=$WORK/${obsid}_sis_bgd_16s_soft2.flc rescale=2.0 plot=no

fselect $WORK/${obsid}_sis_bgd_16s_soft2.flc+1 $WORK/${obsid}_sis_bgd_16s_soft2.flc \
    expr='\!ISNULL(RATE1)' clobber=yes


########################################################################################
# Create light curve  -  background, hard-band, 16s, SIS
########################################################################################

echo "gnt_infrm (t_lcvs)"
echo "gnt_infrm (t_lcvs) Background, hard-band, 16s bin, SIS"

echo $WORK/${obsid}_$sis[1]_bgd_hard.evt >! $WORK/lcvs.tmp
echo $WORK/${obsid}_$sis[2]_bgd_hard.evt >> $WORK/lcvs.tmp

set nbint = `$SCRIPTS/t_lcvs_nbint $sbeg $send 16`
echo nbint = $nbint
if ( $nbint < 2 ) echo "gnt_warng (t_lcvs): Interval ($nbint) is smaller than requested bin size (16)."

echo "gnt_infrm (t_lcvs) Xronos command being used is:"

echo "gnt_infrm (t_lcvs) lcurve nser=1 cfile1=@$WORK/lcvs.tmp"
echo "gnt_infrm (t_lcvs) window=$SCRIPTS/win_full2.wi dtnb=16 nbint=$nbint"
echo "gnt_infrm (t_lcvs) tunits=1 clobber=yes outfile=$WORK/${obsid}_sis_bgd_16s_hard.flc rescale=2.0 plot=no"

lcurve nser=1 cfile1=@$WORK/lcvs.tmp \
    window=$SCRIPTS/win_full2.wi dtnb=16 nbint=$nbint \
    tunits=1 clobber=yes outfile=$WORK/${obsid}_sis_bgd_16s_hard.flc rescale=2.0 plot=no

fselect $WORK/${obsid}_sis_bgd_16s_hard.flc+1 $WORK/${obsid}_sis_bgd_16s_hard.flc \
    expr='\!ISNULL(RATE1)' clobber=yes


########################################################################################
# Create light curve  -  background, hard-band, 16s, GIS
########################################################################################

echo "gnt_infrm (t_lcvs)"
echo "gnt_infrm (t_lcvs) Background, hard-band, 16s bin, GIS"

echo $WORK/${obsid}_$gis[1]_bgd_hard.evt >! $WORK/lcvs.tmp
echo $WORK/${obsid}_$gis[2]_bgd_hard.evt >> $WORK/lcvs.tmp

set nbint = `$SCRIPTS/t_lcvs_nbint $gbeg $gend 16`
echo nbint = $nbint
if ( $nbint < 2 ) echo "gnt_warng (t_lcvs): Interval ($nbint) is smaller than requested bin size (16)."

echo "gnt_infrm (t_lcvs) Xronos command being used is:"

echo "gnt_infrm (t_lcvs) lcurve nser=1 cfile1=@$WORK/lcvs.tmp"
echo "gnt_infrm (t_lcvs) window=$SCRIPTS/win_full2.wi dtnb=16 nbint=$nbint"
echo "gnt_infrm (t_lcvs) tunits=1 clobber=yes outfile=$WORK/${obsid}_gis_bgd_16s_hard.flc rescale=2.0 plot=no"

lcurve nser=1 cfile1=@$WORK/lcvs.tmp \
    window=$SCRIPTS/win_full2.wi dtnb=16 nbint=$nbint \
    tunits=1 clobber=yes outfile=$WORK/${obsid}_gis_bgd_16s_hard.flc rescale=2.0 plot=no

fselect $WORK/${obsid}_gis_bgd_16s_hard.flc+1 $WORK/${obsid}_gis_bgd_16s_hard.flc \
    expr='\!ISNULL(RATE1)' clobber=yes


########################################################################################
# Check if all files were created
########################################################################################

echo "gnt_infrm (t_lcvs)"
echo "gnt_infrm (t_lcvs) checking to see if all light curves were created"
echo "gnt_infrm (t_lcvs)"


#### source region light curves ####

set file = "${obsid}_sis_src_16s.flc"
if ( ! -e $WORK/$file ) echo "gnt_warng (t_lcvs) **** SEQUENCE $obsid - $file was not created ****"

set file = "${obsid}_gis_src_16s.flc"
if ( ! -e $WORK/$file ) echo "gnt_warng (t_lcvs) **** SEQUENCE $obsid - $file was not created ****"

set file = "${obsid}_sis_src_256s.flc"
if ( ! -e $WORK/$file ) echo "gnt_warng (t_lcvs) **** SEQUENCE $obsid - $file was not created ****"

set file = "${obsid}_gis_src_256s.flc"
if ( ! -e $WORK/$file ) echo "gnt_warng (t_lcvs) **** SEQUENCE $obsid - $file was not created ****"

set file = "${obsid}_sis_src_5760s.flc"
if ( ! -e $WORK/$file ) echo "gnt_warng (t_lcvs) **** SEQUENCE $obsid - $file was not created ****"

set file = "${obsid}_gis_src_5760s.flc"
if ( ! -e $WORK/$file ) echo "gnt_warng (t_lcvs) **** SEQUENCE $obsid - $file was not created ****"

set file = "${obsid}_sis_src_16s_full2.flc"
if ( ! -e $WORK/$file ) echo "gnt_warng (t_lcvs) **** SEQUENCE $obsid - $file was not created ****"

set file = "${obsid}_sis_src_256s_full2.flc"
if ( ! -e $WORK/$file ) echo "gnt_warng (t_lcvs) **** SEQUENCE $obsid - $file was not created ****"

set file = "${obsid}_sis_src_5760s_full2.flc"
if ( ! -e $WORK/$file ) echo "gnt_warng (t_lcvs) **** SEQUENCE $obsid - $file was not created ****"

set file = "${obsid}_sis_src_16s_soft.flc"
if ( ! -e $WORK/$file ) echo "gnt_warng (t_lcvs) **** SEQUENCE $obsid - $file was not created ****"

set file = "${obsid}_sis_src_5760s_soft.flc"
if ( ! -e $WORK/$file ) echo "gnt_warng (t_lcvs) **** SEQUENCE $obsid - $file was not created ****"

set file = "${obsid}_sis_src_16s_soft2.flc"
if ( ! -e $WORK/$file ) echo "gnt_warng (t_lcvs) **** SEQUENCE $obsid - $file was not created ****"

set file = "${obsid}_sis_src_5760s_soft2.flc"
if ( ! -e $WORK/$file ) echo "gnt_warng (t_lcvs) **** SEQUENCE $obsid - $file was not created ****"

set file = "${obsid}_sis_src_16s_hard.flc"
if ( ! -e $WORK/$file ) echo "gnt_warng (t_lcvs) **** SEQUENCE $obsid - $file was not created ****"

set file = "${obsid}_sis_src_5760s_hard.flc"
if ( ! -e $WORK/$file ) echo "gnt_warng (t_lcvs) **** SEQUENCE $obsid - $file was not created ****"

set file = "${obsid}_gis_src_16s_hard.flc"
if ( ! -e $WORK/$file ) echo "gnt_warng (t_lcvs) **** SEQUENCE $obsid - $file was not created ****"

set file = "${obsid}_gis_src_5760s_hard.flc"
if ( ! -e $WORK/$file ) echo "gnt_warng (t_lcvs) **** SEQUENCE $obsid - $file was not created ****"


#### background region light curves ####

set file = "${obsid}_sis_bgd_16s.flc"
if ( ! -e $WORK/$file ) echo "gnt_warng (t_lcvs) **** SEQUENCE $obsid - $file was not created ****"

set file = "${obsid}_sis_bgd_5760s.flc"
if ( ! -e $WORK/$file ) echo "gnt_warng (t_lcvs) **** SEQUENCE $obsid - $file was not created ****"

set file = "${obsid}_gis_bgd_16s.flc"
if ( ! -e $WORK/$file ) echo "gnt_warng (t_lcvs) **** SEQUENCE $obsid - $file was not created ****"

set file = "${obsid}_sis_bgd_5760s.flc"
if ( ! -e $WORK/$file ) echo "gnt_warng (t_lcvs) **** SEQUENCE $obsid - $file was not created ****"

set file = "${obsid}_sis_bgd_16s_full2.flc"
if ( ! -e $WORK/$file ) echo "gnt_warng (t_lcvs) **** SEQUENCE $obsid - $file was not created ****"

set file = "${obsid}_sis_bgd_5760s_full2.flc"
if ( ! -e $WORK/$file ) echo "gnt_warng (t_lcvs) **** SEQUENCE $obsid - $file was not created ****"

set file = "${obsid}_sis_bgd_16s_soft.flc"
if ( ! -e $WORK/$file ) echo "gnt_warng (t_lcvs) **** SEQUENCE $obsid - $file was not created ****"

set file = "${obsid}_sis_bgd_16s_soft2.flc"
if ( ! -e $WORK/$file ) echo "gnt_warng (t_lcvs) **** SEQUENCE $obsid - $file was not created ****"

set file = "${obsid}_sis_bgd_16s_hard.flc"
if ( ! -e $WORK/$file ) echo "gnt_warng (t_lcvs) **** SEQUENCE $obsid - $file was not created ****"

set file = "${obsid}_gis_bgd_16s_hard.flc"
if ( ! -e $WORK/$file ) echo "gnt_warng (t_lcvs) **** SEQUENCE $obsid - $file was not created ****"


########################################################################################
# Create full band gif files
########################################################################################

cd $WORK


#### look for flc file ####

if ( -e $WORK/${obsid}_sis_src_5760s.flc ) then

   set flc = $WORK/${obsid}_sis_src_5760s.flc

else

   echo "gnt_warng (t_lcvs) **** SEQUENCE $obsid - file ${obsid}_sis_src_5760s.flc does not exist ****"

   if ( -e $WORK/${obsid}_sis_src_16s.flc ) then
      set flc = $WORK/${obsid}_sis_src_16s.flc
   else
      set flc = ""
      echo "gnt_warng (t_lcvs) **** SEQUENCE $obsid - file ${obsid}_sis_src_16s.flc does not exist."
   endif

endif


#### plot flc file ####

if ( $flc != "" ) then

    echo "gnt_infrm (t_lcvs) creating plot for file $flc"


# get min and max TIME for the plot

    fstatistic infile=$flc rows="-" colname=TIME
    set timemin = `pget fstatistic min`
    set timemax = `pget fstatistic max`
    set timedif = `echo $timemin $timemax | awk '{print 0.2*($2-$1)}'`
    set timemin = `echo $timemin $timedif | awk '{print $1-$2}'`
    set timemax = `echo $timemax $timedif | awk '{print $1+$2}'`
    echo "r x $timemin $timemax" >! $WORK/${obsid}_do_html_flc.tmp

# add error to value to get max column

    fcalc infile=${flc}+1 outfile=${flc}.tmp clname=MIN expr="RATE1-ERROR1" copycol=yes copyall=yes rowrange="-" clobber=yes
    fcalc infile=${flc}.tmp+1 outfile=${flc}.tmp clname=MAX expr="RATE1+ERROR1" copycol=yes copyall=yes rowrange="-" clobber=yes

# get min and max RATE1 for the plot

    fstatistic infile=${flc}.tmp  rows="-" colname=MIN
    set rate1min = `pget fstatistic min`
    fstatistic infile=${flc}.tmp  rows="-" colname=MAX
    set rate1max = `pget fstatistic max`
    set rate1dif = `echo $rate1min $rate1max | awk '{print 0.2*($2-$1)}'`
    set rate1min = `echo $rate1min $rate1dif | awk '{print $1-$2}'`
    set rate1max = `echo $rate1max $rate1dif | awk '{print $1+$2}'`
    echo "r y $rate1min $rate1max" >> $WORK/${obsid}_do_html_flc.tmp

# create plot - do_html_flc.xco has generic commands

    cat $WORK/${obsid}_do_html_flc.tmp $XCO/do_html_flc.xco >! $WORK/${obsid}_do_html_flc.xco
    fplot $flc xparm="TIME[XAX_E]" yparm="RATE1[ERROR1]" rows="-" device="/null" pltcmd="@$WORK/${obsid}_do_html_flc.xco"

endif


# create jpeg file for web page

set file = web_flc

if ( -e ${file}.gif ) then
    echo "gnt_infrm (t_lcvs) creating thumbnail for ${file}.gif"
    convert -resize 170x136 gif:${file}.gif gif:${file}_thumb.gif
    mv ${file}.gif ${file}_full.gif
    mv ${file}_thumb.gif ${file}_full_thumb.gif
else
    echo "gnt_warng (t_lcvs) **** light curve image for web was not produced ****"
endif


########################################################################################
# Create full2 gif files
########################################################################################

cd $WORK


#### remove previous files used in full ####

if ( -e web_flc.gif ) rm web_flc.gif
if ( -e ${obsid}_do_html_flc.tmp ) rm ${obsid}_do_html_flc.tmp
if ( -e ${obsid}_do_html_flc.xco ) rm ${obsid}_do_html_flc.xco


#### look for flc file ####

if ( -e $WORK/${obsid}_sis_src_5760s_full2.flc ) then

   set flc = $WORK/${obsid}_sis_src_5760s_full2.flc

else

   echo "gnt_warng (t_lcvs) **** SEQUENCE $obsid - file ${obsid}_sis_src_5760s_full2.flc does not exist ****"

   if ( -e $WORK/${obsid}_sis_src_16s_full2.flc ) then
      set flc = $WORK/${obsid}_sis_src_16s_full2.flc
   else
      set flc = ""
      echo "gnt_warng (t_lcvs) **** SEQUENCE $obsid - file ${obsid}_sis_src_16s_full2.flc does not exist."
   endif

endif


#### plot flc file ####

if ( $flc != "" ) then

    echo "gnt_infrm (t_lcvs) creating plot for file $flc"


# get min and max TIME for the plot

    fstatistic infile=$flc rows="-" colname=TIME
    set timemin = `pget fstatistic min`
    set timemax = `pget fstatistic max`
    set timedif = `echo $timemin $timemax | awk '{print 0.2*($2-$1)}'`
    set timemin = `echo $timemin $timedif | awk '{print $1-$2}'`
    set timemax = `echo $timemax $timedif | awk '{print $1+$2}'`
    echo "r x $timemin $timemax" >! $WORK/${obsid}_do_html_flc.tmp

# add error to value to get max column

    fcalc infile=${flc}+1 outfile=${flc}.tmp clname=MIN expr="RATE1-ERROR1" copycol=yes copyall=yes rowrange="-" clobber=yes
    fcalc infile=${flc}.tmp+1 outfile=${flc}.tmp clname=MAX expr="RATE1+ERROR1" copycol=yes copyall=yes rowrange="-" clobber=yes

# get min and max RATE1 for the plot

    fstatistic infile=${flc}.tmp  rows="-" colname=MIN
    set rate1min = `pget fstatistic min`
    fstatistic infile=${flc}.tmp  rows="-" colname=MAX
    set rate1max = `pget fstatistic max`
    set rate1dif = `echo $rate1min $rate1max | awk '{print 0.2*($2-$1)}'`
    set rate1min = `echo $rate1min $rate1dif | awk '{print $1-$2}'`
    set rate1max = `echo $rate1max $rate1dif | awk '{print $1+$2}'`
    echo "r y $rate1min $rate1max" >> $WORK/${obsid}_do_html_flc.tmp

# create plot

    cat $WORK/${obsid}_do_html_flc.tmp $XCO/do_html_flc.xco >! $WORK/${obsid}_do_html_flc.xco
    fplot $flc xparm="TIME[XAX_E]" yparm="RATE1[ERROR1]" rows="-" device="/null" pltcmd="@$WORK/${obsid}_do_html_flc.xco"

endif


# create jpeg file for web page

set file = web_flc

if ( -e ${file}.gif ) then
    echo "gnt_infrm (t_lcvs) creating thumbnail for ${file}.gif"
    convert -resize 170x136 gif:${file}.gif gif:${file}_thumb.gif
    mv ${file}.gif ${file}_full2.gif
    mv ${file}_thumb.gif ${file}_full2_thumb.gif
else
    echo "gnt_warng (t_lcvs) **** light curve image for web was not produced ****"
endif


########################################################################################
# All done 
########################################################################################

echo "gnt_infrm (t_lcvs)"
echo "gnt_infrm (t_lcvs) ++++++++++++++++++++++++++++++++++++++++++++++++++++++"
echo "gnt_infrm (t_lcvs)  t_lcvs 4.0  FINISHED "`date`
echo "gnt_infrm (t_lcvs) ++++++++++++++++++++++++++++++++++++++++++++++++++++++"
echo "gnt_infrm (t_lcvs)"