#!/bin/sh

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

legend=0
legend=1

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

#set x

cat /dev/null > pic.made

#------------------------------------------

if [ $# -eq 0 ]; then
   FILE="fort.*"
else
   FILE="$*"
fi

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}'`

freq=`cat fort.2 | awk '{print $1}'`
Qual=`cat fort.2 | awk '{print $2}'`

label1="X"
label2="Y"

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

#----------------------------
wbox=6.5
hbox=`bc -l <<END
      $wbox * ( $nx - 1 ) / ( $ny - 1 )
END`
   ybox=`bc -l <<END
      ( 11 - $hbox ) / 2
END`
   ly=`bc -l <<END
      $ybox + 0.4
END`
#----------------------------

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

   title="$ext (Freq=$freq; Q=$Qual)"

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

#---------------------------
   if [ $i2D -eq 1 ]; then
#---------------------------

   if [ $ifig -eq 1 ]; then

         psimage < $i.bin n1=$nx f1=$ax d1=$hx f2=$ay d2=$hy style=$style \
            whls=$whls bhls=$bhls bps=24 \
            legend=$legend lstyle=vertright lx=1.3 ly=2.3 \
            lwidth=0.2 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 titlecolor=$COL1 \
            title="$title" titlesize=18 \
            label1="$label1" label2="$label2" > $i.ps
   else

      pscontour < $i.bin n1=$nx f1=$ax d1=$hx f2=$ay d2=$hy style=$style \
         labelcf=1 labelcper=1 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 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

#         fc=0.1 dc=0.1 nc=21 \

   fi

#---------------------------
   else
#---------------------------

      nplot=1; n=`wc $i | awk '{print $1}'`; linecolor=black;
      linewidth=1.5; lineon=1; lineoff=0
      title="Trace @($Obsx,$Obsy)"

      if [ "$ext" = "7273" ]; then
         n1=`wc fort.72 | awk '{print $1}'`
         n2=`wc fort.73 | awk '{print $1}'`
         nplot=2; n="$n1,$n2"; linecolor="black,black";
         linewidth="1,2"; lineon="1,4"; lineoff="0,4"
      elif [ "$ext" = "8372" ]; then
         n1=`wc fort.83 | awk '{print $1}'`
         n2=`wc fort.72 | awk '{print $1}'`
         nplot=2; n="$n1,$n2"; linecolor="black,blue";
         linewidth="1.5,2"; lineon="1,4"; lineoff="0,4"
      elif [ "$ext" = "8373" ]; then
         n1=`wc fort.83 | awk '{print $1}'`
         n2=`wc fort.73 | awk '{print $1}'`
         nplot=2; n="$n1,$n2"; linecolor="black,blue";
         linewidth="1.5,2"; lineon="1,4"; lineoff="0,4"
      elif [ "$ext" = "827273" ]; then
         n1=`wc fort.82 | awk '{print $1}'`
         n2=`wc fort.72 | awk '{print $1}'`
         n3=`wc fort.73 | awk '{print $1}'`
         nplot=3; n="$n1,$n2,$n3"; linecolor="black,blue,red";
         echo "Color=$linecolor"
         linewidth="1.5,2,2"; lineon="1,4,2"; lineoff="0,4,3"
      elif [ "$ext" = "837273" ]; then
         n1=`wc fort.83 | awk '{print $1}'`
         n2=`wc fort.72 | awk '{print $1}'`
         n3=`wc fort.73 | awk '{print $1}'`
         nplot=3; n="$n1,$n2,$n3"; linecolor="black,blue,red";
         echo "Color=$linecolor"
         linewidth="1.5,2,2"; lineon="1,4,2"; lineoff="0,4,3"
      fi

      psgraph < $i.bin n=$n nplot=$nplot linecolor=$linecolor \
         linewidth=$linewidth lineon=$lineon lineoff=$lineoff \
         wbox=6.5 hbox=6.5 xbox=1.25 ybox=2.25 \
         style=normal label1="t" label2="U" \
         title="$title" titlesize=20 > ${i}.ps


#---------------------------
   fi
#---------------------------

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

   ;;
   esac

done

