software/JDataQuality/qaqc.sh

Go to the documentation of this file.

#!/usr/bin/env zsh
#
# \author mdejong
#
script=${0##*/}
 
if [ -z $JPP_DIR ]; then
    echo "Variable JPP_DIR undefined."
    exit
fi
 
source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null
 
set_variable:  DEBUG              QAQC_DEBUG    0
set_variable+  UTC                QAQC_UTC      UTC
 
if do_usage $*; then
    usage "$script <QA/QC file> <run range> [(parameter)]"
fi
 
if (( $# < 2 )); then
    fatal "Wrong number of arguments."
fi
 
set_variable  INPUT_FILE   $argv[1]
set_variable  RANGE        $argv[2]
set_array     PARAMETERS   $argv[3,-1]
 
RUNS=(${RANGE%%-*} ${RANGE##*-})
 
 
HEADER=(`head -1 $INPUT_FILE`)
 
if [[ "$PARAMETERS[1]" == "?" ]]; then
    echo $HEADER[*]
    exit 0
fi
 
let "INDEX_OF_RUN = ${HEADER[(ie)run]}"   # index of column "run"
 
INDICES=()
 
for PARAMETER in $PARAMETERS[*]; do
 
    INDEX=${HEADER[(ie)$PARAMETER]}
 
    if (( $INDEX > ${#HEADER} )); then
            fatal "Invalid parameter $PARAMETER."
    else
            INDICES+=($INDEX)
    fi
done
 
if (( ${#INDICES} == 0 )); then
    INDICES=({1..${#HEADER}})
fi
 
 
printf "#"
 
for i in $INDICES[*]; do
        printf "%-32s " $HEADER[$i]
done
 
printf "\n"
 
while read BUFFER; do
 
    ARRAY=(`echo $BUFFER`)
 
    if [[ ${(pj:\0:)HEADER} != ${(pj:\0:)ARRAY} ]]; then
 
        let "RUN = $ARRAY[$INDEX_OF_RUN]"
 
        if (( $RUN >= $RUNS[1]  &&  $RUN <= $RUNS[2] )); then
 
                for i in $INDICES[*]; do
 
                # convert to date
 
                if [[ $HEADER[$i] == ${UTC:-XXX}* ]]; then
                    ARRAY[$i]="$(date '+%d-%m-%Y %H:%M' -d @$ARRAY[$i])"
                fi
 
                    printf " %-32s" $ARRAY[$i]
                done
 
            printf "\n"
        fi
    fi
done < $INPUT_FILE