#!/bin/sh

YEAR=2013
MON=06

cat *${YEAR}${MON}*1min.txt > tmp.txt

TYPES="Pump1 Pump2 AMBIENT1 AMBIENT2"
for TYPE in $TYPES ; do

VARS="VAP D18O DD"
for VAR in $VARS ; do

if [ $VAR == "VAP" ] ; then 
NUM=6
elif [ $VAR == "D18O" ] ; then
NUM=8 
elif [ $VAR == "DD" ] ; then
NUM=10
fi

#grep -i $TYPE tmp.txt | awk '$7<1000&&$9<0.5&&$11<2&&$12>35{print $3+'"${DOY}"'-1+$4/24+$5/60/24+9/24,$'"${NUM}"'}' > tmp_var.txt
grep -i $TYPE tmp.txt | awk '$7<1000&&$9<0.5&&$11<2&&$12>35{print $1"/"$2"/"$3"_"$4":"$5,$'"${NUM}"'}' > tmp_var.txt

cat <<EOF > draw_timeseries.gp
set terminal postscript eps enhanced color 30
set size 2
set size ratio 0.5
set xlabel 'In ${YEAR} (GMT)'
set ylabel '${VAR}'
set xdata time
set timefmt "%Y/%m/%d_%H:%M"
set format x "%m/%d"
set title '1-min Timeseries ${VAR}, ${TYPE}, ${YEAR}/${MON}'
set output '${TYPE}_${VAR}.eps'
plot 'tmp_var.txt' using 1:2 title '${VAR}' w p ps 2
EOF
#set x2data time
#set format x2 "%Y/%m/%d"
#set xtics nomirror
#set timefmt "%Y/%m/%d"
#set xrange [${DOY}:${DOY2}]

gnuplot draw_timeseries.gp
ps2img.sh ${TYPE}_${VAR}.eps ${TYPE}_${VAR}.gif
scp -p  ${TYPE}_${VAR}.gif rainbow.iis.u-tokyo.ac.jp:~/public_html/tmp/Mase/

done
done

exit