#! /bin/bash -f
# Last edited on 2012-09-27 01:25:36 by stolfilocal

tmp="/tmp/$$"
ny=638

export GDFONTPATH=..

for prof in 'rot' 'irr' ; do

  for var in 'pressure' 'speed' 'vorticity' ; do 

    pngfile="out/plot-${prof}-${var}.png"

    gnuplot << EOF
    
    rot_omega(r) = omega
    rot_speed(r) = omega*r
    rot_pressure(r) = 0.5*rho*omega**2*r**2
    rot_vorticity(r) = 2*omega
    rot_shear(r) = 0
    
    rot_omega(r) = (1-exp(-beta*r**2))*Gamma/(2*r**r)
    irr_speed(r) = (1-exp(-beta*r**2))*Gamma/(2*r)
    irr_pressure(r) = 
    irr_vorticity(r) = 
    rot_shear(r) = 
    
    set term png font arial 14 size ${width},${height}
    set xrange [${xmin}:${xmax}]
    set output "${tmp}.png" 
    unset logscale y; unset logscale y2
    set ytics 100; set y2tics 100
    set mytics 2; set my2tics 2
    set format y "%5.0f"
    set format y2 "%-5.0f"
    set xzeroaxis linetype -1 linewidth 1.5

    set lmargin 8; set rmargin 8; set tmargin 0.5; set bmargin 2.0
    unset title
    set grid ytics nomytics
    set bars 0.0
    unset key

    set yrange [ymin:ymax]; set y2range [ymin:ymax]

    plot \

 
EOF

    convert ${tmp}.png -resize '50%' ${pngfile}

    display ${pngfile}

    rm -f ${tmp}.png

  done

done