source: trunk/DataCheck/Archive/fixAuxKeyWords.sh@ 16217

#!/bin/bash

if [ "$#" != 3 ]
then
        echo "Error: fixAuxKeyWords.sh should be called with 2 arguments. Aborting."
        exit
fi

file=$1
tempFile="temp.txt"
reportFile=$2
processErrors=$3

if [ -a $tempFile ]
then
        rm $tempFile 2>/dev/null
fi

#get current keywords value
result=`/home/isdc/lyard/FACT++/fitsdump $file -h -o $tempFile 2>/dev/null`

if [ -a $tempFile ]
then
        timesys=`grep 'TIMESYS' $tempFile | grep -E -o 'UTC'`
        mjdref=`grep 'MJDREF' $tempFile | grep -E -o '[-]*[0-9]*[.]*[0-9]+'`
        tstarti=`grep 'TSTARTI' $tempFile | grep -E -o '[-]*[0-9]*[.]*[0-9]+'`
        tstartf=`grep 'TSTARTF' $tempFile | grep -E -o '[-]*[0-9]*[.]*[0-9]+'`
        tstopi=`grep 'TSTOPI' $tempFile | grep -E -o '[-]*[0-9]*[.]*[0-9]+'`
        tstopf=`grep 'TSTOPF' $tempFile | grep -E -o '[-]*[0-9]*[.]*[0-9]+'`
        date_obs=`grep 'DATE-OBS' $tempFile`
        date_end=`grep 'DATE-END' $tempFile`
        telescope=`grep 'TELESCOP' $tempFile`
        package=`grep 'PACKAGE' $tempFile`
        origin=`grep 'ORIGIN' $tempFile`
        timeunit=`grep 'TIMEUNIT' $tempFile`

else
        echo "Could not list keywords in $file" >> $processErrors
        exit
fi

rm $tempFile

#retrieve the start and stop time from the data itself
result=`/home/isdc/lyard/FACT++/fitsdump $file -c Time --minmax --nozero -o $tempFile 2>/dev/null`

if [ -a $tempFile ]
then

        tstart=`grep 'min' $tempFile | grep -E -o '[-]*[0-9]*[.]*[0-9]+'`
        tstop=`grep 'max' $tempFile | grep -E -o '[-]*[0-9]*[.]*[0-9]+'`
        tstarti2=`echo $tstart | grep -E -o '[-]*[0-9]*[.]' | grep -E -o '[-]*[0-9]*'`
        tstartf2=`echo $tstart | grep -E -o '[.][0-9]+'`
        if [ "$tstarti2" == "" ] #no decimal part
        then
                tstarti2=$tstart
        fi
        tstartf2="0"$tstartf2
        tstopi2=`echo $tstop | grep -E -o '[-]*[0-9]*[.]' | grep -E -o '[-]*[0-9]*'`
        tstopf2=`echo $tstop | grep -E -o '[.][0-9]+'`
        if [ "$tstopi2" == "" ] #no decimal part
        then
                tstopi2=$tstop
        fi
        tstopf2="0"$tstopf2

else
#       echo "Could not minmax $file" >> $processErrors
        exit
fi

rm $tempFile
#output the values to be added/updated to temp text files
modified="false"
if [ "$telescope" == "" ]
then
        echo "TELESCOP FACT / Telescope that acquired this data" >> $tempFile
        modified="true"
        echo "TELESCOP in $file" >> $2
fi
if [ "$package" == "" ]
then
        echo "PACKAGE FACT++ / Package name" >> $tempFile
        modified="true"
        echo "PACKAGE in $file" >> $2
fi
if [ "$origin" == "" ]
then
        echo "ORIGIN FACT / Institution that wrote the file" >> $tempFile
        modified="true"
        echo "ORIGIN in $file" >> $2
fi
if [ "$timeunit" == "" ]
then
        echo "TIMEUNIT d / Time given in days w.r.t. to MJDREF" >> $tempFile
        modified="true"
        echo "TIMEUNIT in $file" >> $2
fi
if [ "$mjdref" != "40587" ]
then
        mjdref="40587"
        echo "MJDREF "$mjdref" / Store times in UNIX time (sec from 1970ish)" >> $tempFile
        modified="true"
        echo "MJDREF in $file" >> $2
fi
if [ "$timesys" != "UTC" ]
then
        timesys="UTC"
        echo "TIMESYS "$timesys" / Time system" >> $tempFile
        modified="true"
        echo "TIMESYS in $file" >> $2
fi
if [ "$tstarti2" != "$tstarti" ] || [ "$tstartf2" != "$tstartf" ]
then
        echo "TSTARTI "$tstarti2" / Time when first event received (integral part)" >> $tempFile
        echo "TSTARTF "$tstartf2" / Time when first event received (fractional part)" >> $tempFile
        date_obs2=`echo "$tstarti2 + $tstartf2 + $mjdref" | bc -l`
        date_obs2=`/home/isdc/lyard/FACT++/MjDtoISO $date_obs2`
        echo "DATE-OBS "$date_obs2" / Time when first event was received" >> $tempFile
        modified="true"
        echo "TSTART in $file" >> $2
fi

if [ "$tstopi2" != "$tstopi" ] || [ "$tstopf2" != "$tstopf" ]
then
        echo "TSTOPI "$tstopi2" / Time when last event received (integral part)" >> $tempFile
        echo "TSTOPF "$tstopf2" / Time when last event received (fractional part)" >> $tempFile
        date_end2=`echo "$tstopi2 + $tstopf2 + $mjdref" | bc -l`
        date_end2=`/home/isdc/lyard/FACT++/MjDtoISO $date_end2`
        echo "DATE-END "$date_end2" / Time when last event was received" >> $tempFile
        modified="true"
        echo "TSTOP  in $file" >> $2
fi

if [ "$modified" == "true" ]
then
        echo "INGEST v0.1 Version of Etienne ingest script" >> $tempFile
fi

if [ -a $tempFile ]
then
        fmodhead $file $tempFile 2>&1 1>/dev/null
fi

fchecksum $file update+ 2>&1 1>/dev/null

rm $tempFile 2>/dev/null