inp=$1
fig=$2
# ***********************************
# Basic condition for GMT drawing
# ***********************************
cat << EOT > _bc.awk
BEGIN{FS=","}
NR==2{xmin=\$1;xmax=\$2;dx=\$3;}
NR==3{ymin=\$1;ymax=\$2;dy=\$3;}
END{
    printf "range=%d/%d/%d/%d\n",xmin,xmax,ymin,ymax
    printf "scale=12/8\n"
    printf "xga=a%.1fg%.1f\n",dx,dx
    printf "yga=a%.1fg%.1f\n",dy,dy
}
EOT
awk -f _bc.awk _draw.txt > _bc.txt
# ***********************************
# Data for pstext
# ***********************************
cat << EOT > _text.awk
BEGIN{FS=",";i=-1}
NR==2{xmin=\$1;xmax=\$2;dx=\$3;}
NR==3{ymin=\$1;ymax=\$2;dy=\$3;}
NR==4{tmax=\$1;wmax=\$2;}
NR==5{tmin=\$1;wmin=\$2;}
6<=NR{i=i+1;y[i]=\$1;s[i]=\$2;}
END{
    n=i
    for(i=0;i<=n;i++){
        gsub(/[_]+/ ," ",s[i])
        printf "%10.3f %10.3f 10p,Helvetica 0 BR %s\n",xmax,y[i],s[i]
    }
    printf "%10.3f %10.3f 10p,Helvetica 0 BL max:%.3fm(%.1fsec)\n",tmax,wmax,wmax,tmax
    printf "%10.3f %10.3f 10p,Helvetica 0 TL min:%.3fm(%.1fsec)\n",tmin,wmin,wmin,tmin
}
EOT
awk -f _text.awk _draw.txt > _d.txt
# ***********************************
# Data for horizontal line of RWL
# ***********************************
cat << EOT > _l1.awk
BEGIN{FS=",";i=-1}
NR==2{xmin=\$1;xmax=\$2;dx=\$3;}
NR==3{ymin=\$1;ymax=\$2;dy=\$3;}
6<=NR{i=i+1;y[i]=\$1;s[i]=\$2;}
END{
    n=i
    printf "%10.3f %10.3f %10.3f %10.3f\n",xmin,y[0],xmax,y[0]
}
EOT
awk -f _l1.awk _draw.txt > _l1.txt
# ***********************************
# Data for horizontal line of the structure
# ***********************************
cat << EOT > _l2.awk
BEGIN{FS=",";i=-1}
NR==2{xmin=\$1;xmax=\$2;dx=\$3;}
NR==3{ymin=\$1;ymax=\$2;dy=\$3;}
6<=NR{i=i+1;y[i]=\$1;s[i]=\$2;}
END{
    n=i
    for(i=1;i<=n;i++){
        printf "%10.3f %10.3f %10.3f %10.3f\n",xmin,y[i],xmax,y[i]
    }
}
EOT
awk -f _l2.awk _draw.txt > _l2.txt
# ***********************************
# GMT drawing
# ***********************************
gmt set FONT_LABEL 12p
gmt set MAP_TICK_LENGTH_PRIMARY -0.1c
source _bc.txt
xlabel="Time (sec)"
ylabel="Water Level (EL.m)"
gmt psbasemap -R$range -JX$scale -Bx$xga+l"$xlabel" -By$yga+l"$ylabel" -BSWne+g#ffffff -P -K > $fig
awk 'BEGIN{FS=","}2<=NR{print $1,$2}'  $inp | gmt psxy -R -J -W1.5p,#000000 -O -K >> $fig
gmt psxy _l1.txt -R -J -Sv0.05+s -W1.0,#777777 -O -K >> $fig
gmt psxy _l2.txt -R -J -Sv0.05+s -W1.0,#777777,10_3_3_3:0 -O -K >> $fig
gmt pstext _d.txt  -R -J -F+f+a+j -N -O >> $fig