#!/usr/bin/gnuplot
#
# NAME
#   vortStatsPlot: gnuplot script for plotting vorticity evolution
#
#
# DESCRIPTION
#   This file can be executed by running gnuplot and typing
#     load vortStatsPlot
#
#
set terminal postscript portrait 8
set output "vortStatsPlot.ps"
set size 0.4,0.3

#set terminal pbm color
#set output "vortStatsPlot.ppm"

set title "Vorticity evolution"
set xlabel "time      (L/U)"
set ylabel "vorticity (U/L)"

set grid
set logscale xy

g1(x) = 946.598 * x**(-1.73481)
g2(x) = 24.2745 * x**(-1.01816)
g3(x) = 469.551 * x**(-1.36965)

plot [5:10000]\
     "turbTowed06-Re10k-Fr1g/stats.dat" every 2 using 1:2  title 'Fr=inf    |vort|' , \
     g3(x) title 't^-1.37', \
     "turbTowed-Re10000-Fr10/stats.dat" every 1 using 1:11 title 'Fr=10  vort_y' , \
     "turbTowed-Re10000-Fr10/stats.dat" every 1 using 1:13 title 'Fr=10  vort_z' , \
     ((x>  5) && (x< 200)) ? g1(x) : 1/0 title 't^-1.73', \
     ((x>200) && (x<8000)) ? g2(x) : 1/0 title 't^-1.02'

# \
     "turbTowed-Re10000-Fr10/stats.dat" using 1:2  title 'Fr=10  |vort|' , \
     "turbTowed-Re10000-Fr10/stats.dat" using 1:9  title 'Fr=10  vort_x' , \
     "turbTowed-Re10000-Fr10/stats.dat" using 1:11 title 'Fr=10  vort_y' , \
     "turbTowed-Re10000-Fr10/stats.dat" using 1:13 title 'Fr=10  vort_z' , \
     444 * x ** -1.33 with lines , \
     "turbTowed06-Re10k-Fr1g/stats.dat" using 1:2  title 'Fr=inf |vort|' , \
     "turbTowed06-Re10k-Fr1g/stats.dat" using 1:9  title 'Fr=inf vort_x' , \
     "turbTowed06-Re10k-Fr1g/stats.dat" using 1:11 title 'Fr=inf vort_y' , \
     "turbTowed06-Re10k-Fr1g/stats.dat" using 1:13 title 'Fr=inf vort_z' , \
     435 * x ** -1.35 with lines

pause -1 "Press return."