#!/bin/sh

if ! which gnuplot > /dev/null 2>&1
then
    echo 'gnuplot not found - skipping graph creation' >&2
    exit 1
fi

gnuplot<<EOF
    set terminal postscript eps monochrome
    set output '../numberDensity.eps'
    set decimalsign '.'

    set format xy '%g'
    set xtics 1e1 mirror
    set xlabel 'v(m^{3})'
    set ytics 1e1 mirror
    set ylabel 'n(m^{-3}m^{-3})'

    set logscale xy
    set xrange [1e-5:1.5]
    set yrange [1e-5:10]
    set key at graph 0.8,0.5

    N0 = 0.05
    v0 = 1.0

    n(x,t) = N0*exp(-t*x**2.0)*2.0*t*v0

    plot '<printf "1 0.05\n1 0\n"' w l t 'Init. Cond.',\
    n(x,10.0) ls 2 lc rgb 'black' t '[Ziff and McGrady, J. Phys. A: Math. Gen., 18: 3027-3037, 1985]',\
    '../numberDensity.transposed.dat' every ::0::19 u 1:2 w p pt 1 t 'air1',\
    '../numberDensity.transposed.dat' every ::20::24 u 1:2 w p pt 5 t 'air2',\
    '../numberDensity.transposed.dat' every ::25::28 u 1:2 w p pt 9 t 'air3'
EOF

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