#!/bin/sh if [ "$1" == "" ] ; then echo Usage: $0 Filename exit fi INFILE=$1 if [ ! -s $INFILE ] ; then echo $INFILE does not exist exit fi TYPES="Pump1 Pump2 AMBIENT1 AMBIENT2" for TYPE in $TYPES ; do grep $TYPE $INFILE > tmp.$$ TIMES=`awk -F'[,:/]' '{print $2 $3 $4 $5 $6}' tmp.$$ | awk '{print $1 $2}' | uniq` for TIME in $TIMES ; do YEAR=`echo $TIME | cut -c1-4 ` MON=`echo $TIME | cut -c5-6 ` DAY=`echo $TIME | cut -c7-8 ` HOUR=`echo $TIME | cut -c9-10 ` MIN=`echo $TIME | cut -c11-12 ` #echo $TIME $YEAR $MON $DAY $HOUR $MIN awk -F'[,:/]' '$2=='"${YEAR}"'&&$3=='"${MON}"'&&$4=='"${DAY}"'&&$5=='"${HOUR}"'&&$6=='"${MIN}"'{print $9,$10,$11}' tmp.$$ > awk.$$ || exit 8 awk '{SUM=SUM+$1}{SUM2=SUM2+$1*$1}{OSUM=OSUM+$2}{OSUM2=OSUM2+$2*$2}{DSUM=DSUM+$3}{DSUM2=DSUM2+$3*$3}END{print SUM/NR,(SUM2/NR-(SUM/NR)**2)**0.5,OSUM/NR,(OSUM2/NR-(OSUM/NR)**2)**0.5,DSUM/NR,(DSUM2/NR-(DSUM/NR)**2)**0.5,NR}' awk.$$ > awk2.$$ ; rm awk.$$ read VAP SVAP DO SDO DD SDD NR < awk2.$$ ; rm awk2.$$ echo $YEAR $MON $DAY $HOUR $MIN $VAP $SVAP $DO $SDO $DD $SDD $NR $TYPE done rm tmp.$$ done exit