|
|
awk array problem
07-21-2011
, so ill try it 1st thing tomorrow.
07-21-2011
07-22-2011
CA=0 ; CB=0 ; CX=0 ; CD=0
FOUNDA=1; FOUNDB=1;
FOUNDCX=1; FOUNDD=1;
for(N=0; N<L; N++)
{
if(A[N] != C[N]) FOUNDA=0;
if(B[N] != C[N]) FOUNDB=0;
if(X[N] != C[N]) FOUNDX=0;
if(D[N] != C[N]) FOUNDD=0;# this has been changed to operate under the assumption that the 1st beam
# opens and closes before next beam opens due to spacing between beams being increased.
BEGIN { OT=0 # Time of previous measurement
MAX=0 # Max num of seconds between valid events
DAY=""; # Current day
CA=0 ; CB=0 ; CX=0 ; CD=0
# Running total of bats leaving and entering
TOTALBATS=0;
# The highest TOTALBATS has ever been
MAXBATS=0;
# Length of the patterns
L=4
# Patterns to check against
# Block 1 unBlock 0 block 0 Unblock 0
A[0]="1,1"; A[1]="1,0"; A[2]="0,1"; A[3]="0,0";
# Block 1 Block 0 unblock 1 Unblock 0
X[0]="1,1"; X[1]="0,1"; X[2]="1,0"; X[3]="0,0";
# Block 0 Unblock 0 Block 1 Unblock 1
B[0]="0,1"; B[1]="0,0"; B[2]="1,1"; B[3]="1,0";
# Block 0 block 1 unBlock 0 unblock 1
D[0]="0,1"; D[1]="1,1"; D[2]="0,0"; D[3]="1,0";
}
function print_daily(day,total,max,min,maxtime)
{
I=total; if(I<0) I=-I;
MX="no maximum"
if(maxtime > 0)
MX=sprintf("peak was at %s", strftime("%H:%M:%S",maxtime));
# printf("COUNT@%s COUNT %+d RET %d LEFT %d GUESS %d (%s)\n",
printf("Date @%s COUNT %+d IN %d OUT %d Estimate %d bats(%s)\n",
day, total, max, -min, I, MX) > "/dev/stderr";
# Reset daily counts
TOTALBATS=0; MAXBATS=0; MINBATS=0; MAXTIME=0;
MINTIME=0;
}
{ # Calculate timestamp from date string
T=mktime($1 " " $2 " " $3 " " $5 " " $6 " " $7);
T+=(60*60*16); # Add sixteen hours
$1=strftime("%Y", T); # Put these back in the strings
$2=strftime("%m", T);
$3=strftime("%d", T);
$5=strftime("%H", T);
$6=strftime("%M", T);
$7=strftime("%S", T);
# When the year, month, and/or day changes, time to print daily counts
if((DAY != $1 "-" $2 "-" $3) && (DAY != ""))
print_daily(DAY,TOTALBATS,MAXBATS,MINBATS,MAXTIME);
DAY=$1 "-" $2 "-" $3;
if($8 == "pv") # Ignore anything but PV lines.
{
# If too much time has passed since the last event, start over.
if((T-OT) > MAX) # Blank the array
for(N=0; N<(L-1); N++) C[N]="";
else # Shift elements toward the front
for(N=0; N<(L-1); N++) C[N]=C[N+1];
OT=T # Set prev time to this one.
C[L-1]=$9 "," $10; # Set the latest event in the array
# Search for events in the array.
FOUNDA=1; FOUNDB=1;
FOUNDX=1; FOUNDD=1;
for(N=0; N<L; N++)
{
if(A[N] != C[N]) FOUNDA=0;
if(B[N] != C[N]) FOUNDB=0;
if(X[N] != C[N]) FOUNDX=0;
if(D[N] != C[N]) FOUNDD=0;
}
# Count the events and mark the hour they occurred in
if(FOUNDA || FOUNDX)
{
printf("A@%s-%s-%s %s:%s:%s\n",$1,$2,$3,$5,$6,$7);
CA++; AH[$5]++;
TOTALBATS++;
}
if(FOUNDB || FOUNDD)
{
printf("B@%s-%s-%s %s:%s:%s\n",$1,$2,$3,$5,$6,$7);
CB++; BH[$5]++;
TOTALBATS--;
}
# Update our maximum daily counts
if(MAXBATS < TOTALBATS)
{
MAXBATS=TOTALBATS;
MAXTIME=T;
}
if(MINBATS > TOTALBATS)
{
MINBATS=TOTALBATS;
MINTIME=T;
}
}
}
END { # The final statistics will be printed to stderr, to easily
# seperate them from the event times printed to stdout.
# The last daily count
print_daily(DAY,TOTALBATS,MAXBATS,MINBATS,MAXTIME);
# Print the event counts
printf("A %2d\nB %2d\nT %2d\n", CA, CB, CA+CB) > "/dev/stderr";
# Print a list of hours from 1-23
STR="H";
for(N=1; N<=23; N++) STR=STR sprintf(" %2d", N);;
print STR > "/dev/stderr";
# Print hourly counts for event A
STR="A";
for(N=1; N<=23; N++)
STR=STR sprintf(" %2d", AH[sprintf("%02d", N)]);
print STR > "/dev/stderr";
# Hourly counts for event B
STR="B";
for(N=1; N<=23; N++)
STR=STR sprintf(" %2d", BH[sprintf("%02d",N)]);
print STR > "/dev/stderr";
}.\zsh.exe .\IpEther1.zsh echo ==================== echo Done. PAUSE
# because windows CMD doesn't actually expand * into arguments for you, # we need to run this in sh to do so. Bleh. exec ./zsh.exe IpEther2.zsh * > events.txt 2> stats.txt # operates on any dir name starting between year 2000 and 2999 # "2>" redirects sterr to stats.txt
#!/bin/zsh # this script is to be used with data exported from Ipether Access database ./gawk.exe -F "," -f batmon5.awk < output.txt # calls gawk using output.txt as the input to awk
|
|
|
|