#!/bin/sh
############################################################
#to   draw timeseries graph
#by   2010/03/31, hanasaki, NIES: H08ver1.0
############################################################
# GMT command preference
############################################################
PSBASEMAP=psbasemap       # default
PSTEXT=pstext             # default
PSXY=psxy                 # default
#
#PSBASEMAP="GMT psbasemap" # ubuntu
#PSTEXT="GMT pstext"       # ubuntu
#PSXY="GMT psxy"           # ubuntu
############################################################
#
############################################################
EPS=temp.eps
JPG=temp.jpg
############################################################
#
############################################################
if [ $# -lt 5 ]; then
  echo htdrawts FILE EPS RFLAG BLFAG TITLE
  exit
fi
#
FILE=$1
EPS=$2
RFLAG=$3
BFLAGARG=$4
TITLE=$5
#
if [ $# -eq 6 ]; then
  NUM=$6
  NUM=`echo $NUM | awk '{print $1+3}'`
  echo 'thick line for ' $NUM 'th column.'
else
  NUM=0
fi
############################################################
#
############################################################
LINES=`wc $FILE | awk '{print $1}'`
WORDS=`wc $FILE | awk '{print $2}'`
COLUMNS=`echo $WORDS $LINES | awk '{print $1/$2}'`
echo $LINES $WORDS $COLUMNS
if [ $COLUMNS -ge 5 ]; then
  sed -e '1d' $FILE > temp.htdrawts.txt
else
  cp $FILE temp.htdrawts.txt
fi
#
XMIN=`echo $RFLAG | sed -e 's-/- -g' | sed -e 's/-R//' | awk '{print $1}'`
XMAX=`echo $RFLAG | sed -e 's-/- -g' | sed -e 's/-R//' | awk '{print $2}'`
YMIN=`echo $RFLAG | sed -e 's-/- -g' | sed -e 's/-R//' | awk '{print $3}'`
YMAX=`echo $RFLAG | sed -e 's-/- -g' | sed -e 's/-R//' | awk '{print $4}'`
echo $XMIN $XMAX $YMIN $YMAX
#
YTIT=`echo $YMAX | awk '{print $1*1.1}'`
XTIT=`echo $XMAX | awk '{print $1*0.5}'`
XANO=`echo $XMAX | awk '{print $1*0.25}'`
YANO=`echo $YMAX | awk '{print $1*0.2}'`
JFLAG=-JX10.5/10.5
BFLAG=-Ba${XANO}::/a${YANO}::neWS
if [ "$BFLAGARG" != "." ]; then
  BFLAG=$BFLAGARG
  echo bflag changed
fi
echo $BFLAG
############################################################
#
############################################################
$PSBASEMAP $RFLAG $JFLAG $BFLAG -K > $EPS
COLUMN=4
while [ $COLUMN -le $COLUMNS ]; do
  if [ $COLUMN = $NUM ]; then
    WFLAG=-W10
    awk '{print NR,$'$COLUMN'}' temp.htdrawts.txt
  else
    WFLAG=-W1
  fi
  echo $COLUMN $WFLAG
  awk '{print NR,$'$COLUMN'}' temp.htdrawts.txt | head -10
  awk '{print NR,$'$COLUMN'}' temp.htdrawts.txt | \
  $PSXY -O $RFLAG $JFLAG $BFLAG $WFLAG -K >> $EPS
  COLUMN=`expr $COLUMN + 1`
done
############################################################
#
############################################################
$PSTEXT -O $RFLAG $JFLAG $BFLAG -N  << EOF >> $EPS
$XTIT $YTIT 20 0 0 6 $TITLE
EOF