#!/bin/bash
today=`date -u +%Y%m%d`
echo doing $today
cd $HOME/spectro_radiometer
file=$1$today-spec.csv
tail -10 $file >tmp$$
awk '
BEGIN {FS=", *"}
/./ {nrec++
	for (k = 10; k <= 2048; k++)
	{
		specdata[k-10] += $k
	}
  }
END {
	ind = 1
	for (s = -BW/2; s <= BW/2; s += (BW/2048))
	{
		d = s/FC
		d *= -1
		d *= 299792
		printf ("%f %f\n", d, specdata[ind]/nrec)
		ind += 1
	}
		
}
' FC=1420.4058e6 BW=3.0e6 <tmp$$ |sort -n >graphdata$$
curd=`date -u +%Y-%m-%d/%H:%M:%S`
gnuplot <<!EOF!
set xrange [-130:130]
set yrange [-1:1.5]
set grid
set xlabel 'Red Shift (km/sec)'
set ylabel 'Relative Brighness (dB)'
set title  'Neutral Hydrogen Spectrum: $curd UTC'
set terminal png large size 700,500
set output 'live_spec.png'
plot 'graphdata$$' with lines title 'H1'
!EOF!

rm -f tmp$$ graphdata$$

ftp -n <<!EOF!
open www.ccera.ca
user ccera $2
cd public_html/files/live_data/h1
put live_spec.png
bye
!EOF!