#!/bin/sh
# see http://slackology.net/wxnotes.html for more info

TD="data"
BASE_DIR="/home/ghz/wx/particle"
DATA_DIR="$BASE_DIR/data"
PLOT_DIR="$BASE_DIR/plots"
AVG_FILTER=$BASE_DIR/avg_filter
AVG_F_64=$BASE_DIR/avg_f_64
DAY_AVG=$BASE_DIR/day_avg
NSAMP=4000
OF="$DATA_DIR/$TD"
DA_OF="$DATA_DIR/data.day_avg"

#from cronjob
FLAG="$BASE_DIR/CRONFLAG"

cd $BASE_DIR || exit 1

#DATEB=../int/date
DATEB=fbsd_date/date

while true; do

DATE=`$DATEB +%Y%m%d`
DATE8H=`$DATEB -v-8H +%Y%m%d%H`
YDATE=`$DATEB -v-1d +%Y%m%d`
YYDATE=`$DATEB -v-2d +%Y%m%d`
YDATEH=`$DATEB -v-1d +%Y%m%d%H`

TD_DUMP="cat $DATA_DIR/data.$YYDATE $DATA_DIR/data.$YDATE $DATA_DIR/data.$DATE"
TD_DUMP2="cat $DATA_DIR/data.$YDATE $DATA_DIR/data.$DATE"

if [ -e $FLAG ]; then
	(echo "$YDATE"; cat $DATA_DIR/data.$YDATE |$DAY_AVG) | paste - - >> $DA_OF
	rm $FLAG
fi
	
	$TD_DUMP | grep -A $NSAMP $YDATEH > $OF.24
	$TD_DUMP | grep -A $NSAMP -B 16 $YDATEH | $AVG_FILTER > $OF.24.avg
	$TD_DUMP | grep -A $NSAMP -B 64 $YDATEH | $AVG_F_64 > $OF.24.avg.64

	$TD_DUMP2 | grep -A $NSAMP $DATE8H > $OF.8
	$TD_DUMP2 | grep -A $NSAMP -B 16 $DATE8H | $AVG_FILTER > $OF.8.avg
	$TD_DUMP2 | grep -A $NSAMP -B 64 $DATE8H | $AVG_F_64 > $OF.8.avg.64

	cat $DA_OF |tail -n 45 > $DA_OF.45

	CPM8="\"$OF.8\" using 1:2 title 'Particle Counts (cpm)' with lines lt 2"
	AVG_CPM8="\"$OF.8.avg\" using 1:2 title '16 pt Running AVG (cpm)' with lines lt 8"
	AVG_648="\"$OF.8.avg.64\" using 1:2 title '64 pt Running AVG (cpm)' with points lt 1"
	AVG_648l="\"$OF.8.avg.64\" using 1:2 title '64 pt Running AVG (cpm)' with lines lt 1"

	CPM24="\"$OF.24\" using 1:2 title 'Particle Counts (cpm)' with lines lt 2"
	AVG_CPM24="\"$OF.24.avg\" using 1:2 title '16 pt Running AVG (cpm)' with lines lt 8"
	AVG_6424="\"$OF.24.avg.64\" using 1:2 title '64 pt Running AVG (cpm)' with points lt 1"
	AVG_6424l="\"$OF.24.avg.64\" using 1:2 title '64 pt Running AVG (cpm)' with lines lt 1"

	DA_BAR="\"$DA_OF\" using 1:2 title 'Daily Averages (cpm)' with histep lt 3"
	DA_BAR_45="\"$DA_OF.45\" using 1:2 title 'Daily Averages (cpm)' with boxes lt 1"

	$BASE_DIR/gen_index

	PROLOUGE="set title \"Radioactivity\";\
		set xtics ;\
		set y2tics ;\
		set key outside below;\
		set xlabel \"Time (UTC)\" 0.0, -1.0;\
		set xdata time;\
		set timefmt \"%Y%m%d%H%M%S\";\
		set grid;\
		set term png size 1235, 480;"

	echo "$PROLOUGE\
		set title \"Radioactivity over the last ~8 hours.\";\
		set xtics 3600 ;\
		set ylabel \"cpm\";\
		set output '$PLOT_DIR/particle_cpm_8.png';\
		plot $CPM8, $AVG_CPM8, $AVG_648;\
		set output '$PLOT_DIR/particle_cpm_8_avg.png';\
		plot $AVG_CPM8, $AVG_648l;\
		" |gnuplot

	echo "$PROLOUGE\
		set title \"Radioactivity over the last ~24 hours.\";\
		set xtics 3600 ;\
		set ylabel \"cpm\";\
		set output '$PLOT_DIR/particle_cpm_24.png';\
		plot $CPM24, $AVG_CPM24, $AVG_6424;\
		set output '$PLOT_DIR/particle_cpm_24_avg.png';\
		plot $AVG_CPM24, $AVG_6424l;\
		" |gnuplot

	echo "$PROLOUGE\
		set title \"Radioactivity: Daily averages.\";\
		set ylabel \"cpm\";\
		set xlabel \"Date (yyyy/mm, UTC)\" 0.0, -1.0;\
		set format x \"%Y\/%m\" ;\
		set output '$PLOT_DIR/particle_cpm_DA.png';\
		set mxtics 2 ;\
		set grid mxtics;\
		plot $DA_BAR;\
		set grid nomxtics;\
		set xtics 172800 ;\
		set mxtics 2 ;\
		set title \"Radioactivity: Daily averages for the last 45 days.\";\
		set xlabel \"Date (mm/dd, UTC)\" 0.0, -1.0;\
		set format x \"%m\/%d\" ;\
		set output '$PLOT_DIR/particle_cpm_DA.45.png';\
		plot $DA_BAR_45;\
		" |gnuplot

	cp $PLOT_DIR/*.png $PLOT_DIR/particles.html /wx/

	rsync -u $BASE_DIR/graph $BASE_DIR/sort_data $BASE_DIR/gen_index \
		 $AVG_FILTER $AVG_F_64 $DAY_AVG /wx/scripts/particlewx/

#	just too undercooked to post right now.
#	rsync -u $BASE_DIR/geigert.c $BASE_DIR/makefile \
#		 /wx/progs/particlewx/

	rsync -qu $BASE_DIR/data/* /wx/data/particlewx/

	echo -n .
	sleep 60
done