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
 
if do_usage $*; then
    usage "$script <QA/QC file> [(parameter)]"
fi
 
if (( $# == 0 )); then
    fatal "Wrong number of arguments."
fi
 
set_variable  INPUT_FILE   $argv[1]
set_array     PARAMETERS   $argv[2,-1]
 
HEADER=(`head -1 $INPUT_FILE`)
 
if [[ "$PARAMETERS[1]" == "?" ]]; then
    echo $HEADER[*]
    exit 0
fi
    
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
 
while read BUFFER; do
 
    ARRAY=(`echo $BUFFER`)
 
    if [[ ${(pj:\0:)HEADER} != ${(pj:\0:)ARRAY} ]]; then
 
        for i in $INDICES[*]; do
            printf "%-24s %-32s\n" $HEADER[$i] $ARRAY[$i]
        done
    fi
done < $INPUT_FILE