#! /bin/sh

WIDTH=300
WIDTHOFF=50
WIDTHOFF2=50
HEIGHT=600
HEIGHTOFF=50
HEIGHTOFF2=250
HEIGHTOFF3=500
HEIGHT2=200
WIDTH2=600
WIDTHOFF3=100
WIDTH3=800
HEIGHT3=300


vrmst=vrmst.bin
vintt=vintt.bin
vinttav=vinttav.bin
vintzav=vintzav.bin
rmspar=vrmst.par
intpar=vintt.par
vinttplot=vinttplot.par
vintzplot=vintzplot.par
unipar=unisam.par
vinttuni=vinttuni.bin
vintxz=vintxz.bin
vintzx=vintzx.bin

# plot rms velocities
ximage par=$rmspar  legend=1 title="rms velocities" <  $vrmst \
wbox=${WIDTH2} hbox=${HEIGHT2} xbox=${WIDTHOFF} ybox=${HEIGHTOFF} &


# plot interval velocities
ximage par=$rmspar  legend=1 title="interval velocities" <  $vintt \
wbox=${WIDTH2} hbox=${HEIGHT2} xbox=${WIDTHOFF} ybox=${HEIGHTOFF2} &

# plot average velocity
xgraph par=$vinttplot legend=1 title="average interval v(t)" <  $vinttav &
xgraph par=$vintzplot legend=1 title="average interval v(z)" <  $vintzav &

# plot interpolated interval velocities
ximage par=$unipar  legend=1 title="interpolated interval velocities" <  $vinttuni \
wbox=${WIDTH2} hbox=${HEIGHT2} xbox=${WIDTHOFF} ybox=${HEIGHTOFF3} &

# plot interpolated interval velocities
xcontour par=$unipar nc=50 d1=3 d2=12.5  legend=1 title="interpolated interval velocities" <  $vinttuni \
wbox=${WIDTH3} hbox=${HEIGHT3} xbox=${WIDTHOFF3} ybox=${HEIGHTOFF3} &



