#!/bin/sh

ifig=2  # contour plot
ifig=1  # density plot

whls=1.0,1.0,1.0;  bhls=0.0,0.0,0.0; COL1=black # black-white (for density plot)
whls=1.0,0.97,1.0; bhls=0.0,0.5,1.0; COL1=blue  # color (for density plot)

#set x

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}'`
ny=`cat fort.1 | awk '{print $5}'`
ay=`cat fort.1 | awk '{print $6}'`
by=`cat fort.1 | awk '{print $7}'`
hy=`cat fort.1 | awk '{print $8}'`

nt=`cat fort.2 | awk '{print $1}'`
at=`cat fort.2 | awk '{print $2}'`
bt=`cat fort.2 | awk '{print $3}'`
ht=`cat fort.2 | awk '{print $4}'`
theta=`cat fort.2 | awk '{print $5}'`
sigma=`cat fort.2 | awk '{print $6}'`
epsilon=`cat fort.2 | awk '{print $7}'`

label1="X"
label2="Y"

x1beg=$ax
x1end=$bx
x2beg=$ay
x2end=$by

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 "$i.ps"

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

   if [ "$i" = "fort.U0" ]; then
      title="U0 (T=$at)"
   else
      echo " Dx=Dy=$hy, check it"
      title="$scheme (T=$bt; sigma=$sigma; Dx=Dy=$hy)"
   fi

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

   if [ $ifig -eq 1 ]; then

      psimage < $i.bin n1=$nx f1=$ax d1=$hx f2=$ay d2=$hy style=normal \
         whls=$whls bhls=$bhls bps=24 \
         legend=1 lstyle=vertright lx=6.85 ly=2.5 \
         lwidth=0.3 lheight=3.0 d1num=$d1num \
         x1beg=$x1beg x1end=$x1end x2beg=$x2beg x2end=$x2end \
         wbox=6.5 hbox=6.5 xbox=1.25 ybox=2.25 \
         labelccolor=$COL1 axescolor=$COL1 \
         title="$title" titlesize=18 titlecolor=$COL1 \
         label1="$label1" label2="$label2" > $i.ps

   else

      pscontour < $i.bin n1=$nx f1=$ax d1=$hx f2=$ay d2=$hy style=normal \
         fc=0.1 dc=0.1 \
         nc=21 labelcf=5 labelcper=5 cwidth=1.5 labelcsize=16 \
         x1beg=$x1beg x1end=$x1end x2beg=$x2beg x2end=$x2end \
         wbox=6.5 hbox=6.5 xbox=1.25 ybox=2.25 \
         grid1=dot grid2=dot d1num=0.2 d2num=0.2 brgb=1.0,0.0,0.0 \
         ccolor=red labelccolor=blue axescolor=blue \
         title="$title" titlesize=18 titlecolor=blue \
         label1="$label1" label2="$label2" > $i.ps

    fi

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

   ;;
   esac

done

