#!/bin/sh

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

cp /dev/null pic.made

nx=`cat fort.1 | awk '{print $1}'`
ny=`cat fort.1 | awk '{print $2}'`
nz=`cat fort.1 | awk '{print $3}'`
xmin=`cat fort.1 | awk '{print $4}'`
dx=`cat fort.1 | awk '{print $5}'`
ymin=`cat fort.1 | awk '{print $6}'`
dy=`cat fort.1 | awk '{print $7}'`
zmin=`cat fort.1 | awk '{print $8}'`
dz=`cat fort.1 | awk '{print $9}'`

n1=$nz; f1=$zmin; d1=$dz
label1="Z (Km)"
title="Traveltime (Sec)"
#title="Traveltime (y=3000)"
style=seismic

fc=0; dc=0.1; nc=81
ccolor=red
labeleach=4; cwidth=1.2

if [ $nz -eq 1 ]; then
   FILE="fort.TTx fort.ERx"
   n1=$ny; f1=$ymin; d1=$dy
   label1="Z (Km)"
   n2=$nx; f2=$xmin; d2=$dx
   label2="X (Km)"
   fc=0; dc=0.05; nc=101
fi

for i in $FILE; do
   case $i in
   *.ps) ;;
   fort.1) ;;
   fort.50) ;;
   *)

   if [ $nz -ne 1 ]; then
      if [ "$i" = "fort.TTx" -o "$i" = "fort.ERx" ]; then
         n2=$nx; f2=$xmin; d2=$dx
         label2="X (Km)"
      else
         n2=$ny; f2=$ymin; d2=$dy
         label2="Y (Km)"
      fi
   fi

   if [ "$i" = "fort.ERx" ]; then
      fc=-20; dc=4; nc=121; labeleach=1
      title="Error (ms)"
   fi

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

   pscontour < $i.bin verbose=1 \
     n1=$n1 n2=$n2 f1=$f1 d1=$d1 f2=$f2 d2=$d2 \
     fc=$fc dc=$dc nc=$nc ccolor=$ccolor \
     labelcper=$labeleach labelcsize=15 cwidth=$cwidth \
     wbox=6.5 hbox=6.5 xbox=1.25 ybox=2.25 \
     label1="$label1" label2="$label2" title="$title" style=$style > $i.ps

   echo "OUTPUT: ${i}.ps"
   echo "${i}.ps" >> pic.made

   rm $i.bin

   ;;
   esac

done

