#!/bin/sh

#set x

if [ -f fort.true -a -f fort.sol ]; then
    cat fort.true fort.sol > fort.TS
fi

if [ $# -eq 0 ]; then
   FILE="fort.*"
else
   FILE="$*"
fi
cat /dev/null > pic.made


scheme=`cat fort.0 | awk '{print $1}'`

nx=`cat fort.1 | awk '{print $1}'`
ax=`cat fort.1 | awk '{print $2}'`
bx=`cat fort.1 | awk '{print $3}'`
hx=`cat fort.1 | awk '{print $4}'`

nt=`cat fort.1 | awk '{print $5}'`
at=`cat fort.1 | awk '{print $6}'`
bt=`cat fort.1 | awk '{print $7}'`
ht=`cat fort.1 | awk '{print $8}'`
theta=`cat fort.1 | awk '{print $9}'`

label1="X"
label2="U"

x1beg=$ax
x1end=$bx

for i in $FILE; do
   case $i in
   *.ps|*.eps|*.out) ;;
   *.c|*.cpp|*.f|*.f90) ;;
   *.o|*.h) ;;
   fort.0|fort.1|fort.2) ;;
   xx|pic|pic.made) ;;
   *)

   ext="`echo $i | sed 's/fort\.//g'`"
   echo "pic: $ext"

   nplot=1; linecolor=red;
   linewidth=1.5; lineon=1; lineoff=0

   title="$scheme (T=$bt; Dt=$ht; Dx=$hx)"
#   title="$scheme (T=$bt; Dt=$ht; Dx=$hx; theta=$theta)"
   n=`wc -l $i | awk '{print $1}'`


      a2b < ${i} n1=2 > $i.bin


      if [ "$i" = "fort.TS" ]; then
          n=`wc -l fort.true | awk '{print $1}'`
          nplot=2; linecolor=red,blue
          linewidth=1,1.5; lineon=1,3; lineoff=0,5
      fi

      psgraph < $i.bin n=$n nplot=$nplot linecolor=$linecolor \
         linewidth=$linewidth lineon=$lineon lineoff=$lineoff \
         x2beg=-0.2 x2end=1.2 \
         wbox=6.5 hbox=6.5 xbox=1.25 ybox=2.25 \
         style=normal titlecolor=blue title="$title" titlesize=20 \
         label1="$label1" label2="$label2" > ${i}.ps

#         linewidth=2,1 mark=0,1 marksize=4,4 \
#         x1beg=$f1 x1end=3  x2beg=0 x2end=0.005 \


   echo $i.ps >> pic.made
   rm $i.bin

   ;;
   esac

done

