#!/usr/bin/gnuplot

reset

#set terminal postscript enhanced
#set output "start-speed-distrib-v_drift.ps"
set terminal postscript eps enhanced 16
set output "start-speed-distrib-v_drift.gpi.eps"

set line style 1 linetype 7 linewidth 1
set line style 2 linetype 1 linewidth 3

set multiplot
set size 1.0,1.0
set lmargin 8
set rmargin 1


set mytics 2
set mxtics 2

set format y "%2.3f"

#set grid

# ------------------------------------------------------------------------
#     +-----+  +-----+
#     |  X  |  |     |
#     +-----+  +-----+
#     +-----+  +-----+
#     |     |  |     |
#     +-----+  +-----+
set origin 0.0,0.5
set size 0.5,0.5
set key bottom
set ytics 0.2

set xlabel "{/Symbol t}\n a)"
set ylabel "u_d({/Symbol t})"
#set label  "{/Symbol \341z\361 / t}" at 3,0.6

set label 1     "{/Symbol a} = 0.1" at 5,0.35
set label 2     "{/Symbol a} = 0.2" at 5,-0.05
set label 5     "{/Symbol a} = 0.5" at 5,0.15
set label 10    "{/Symbol a} = 1.0"  at 5,0.52
set label 20    "{/Symbol a} = 2.0"  at 5,0.65
set label 50    "{/Symbol a} = 5.0"  at 5,0.75
set label 10000 "{/Symbol a} = 10^3" at 5,0.85

plot [:8][-0.15:0.95] \
	"N-5e5_T-300_alpha-1000_beta-1.0.dat" using 2:6 notitle  with lines 1, \
	"N-5e5_T-300_alpha-5_beta-1.0.dat"    using 2:6 notitle  with lines 1, \
	"N-5e5_T-300_alpha-2_beta-1.0.dat"    using 2:6 notitle  with lines 1, \
	"N-5e5_T-300_alpha-1_beta-1.0.dat"    using 2:6 notitle  with lines 1, \
	"N-5e5_T-300_alpha-0.5_beta-1.0.dat"  using 2:6 notitle  with lines 1, \
	"N-5e5_T-300_alpha-0.2_beta-1.0.dat"  using 2:6 notitle  with lines 1, \
	"N-5e5_T-300_alpha-0.1_beta-1.0.dat"  using 2:6 notitle  with lines 1

#	"N-5e5_T-300_alpha-50_beta-1.0.dat"   using 2:6 title "{/Symbol a} = 50" with lines 2, \
#	"N-5e5_T-300_alpha-10_beta-1.0.dat"   using 2:6 title "{/Symbol a} = 10" with lines 2, \
#	sqrt(2/pi) title 'u_d({/Symbol \245})' with lines 4

#set nolabel # remove <\dzeta> / \ tau from graphs

# ------------------------------------------------------------------------
#     +-----+  +-----+
#     |     |  |     |
#     +-----+  +-----+
#     +-----+  +-----+
#     |  X  |  |     |
#     +-----+  +-----+
set origin 0.0,0.0
set size 0.5,0.5
set key bottom
#set ytics 0.001
set xlabel "{/Symbol t}\n c)"
set ylabel "u_d({/Symbol t})"

#set nolabel 1
#set nolabel 10

plot [:8][0.795:0.801] \
	"analytic.dat" title "teoret." with lines 2, \
	sqrt(2/pi) title 'u_d({/Symbol \245})' with lines 4

# ------------------------------------------------------------------------
#     +-----+  +-----+
#     |     |  |  X  |
#     +-----+  +-----+
#     +-----+  +-----+
#     |     |  |     |
#     +-----+  +-----+
set origin 0.5,0.5
set size 0.5,0.5
set lmargin 9
set rmargin 0
set key top
#set ytics 0.005
set xlabel "{/Symbol t}\n b)"
set ylabel "u_d({/Symbol t})"

set label 1  "{/Symbol a} = 0.1" at 5,0.35
set label 2  "{/Symbol a} = 0.2" at 5,-0.05
set label 5  "{/Symbol a} = 0.5" at 5,0.15
set label 10 "{/Symbol a} = 1.0" at 5,0.52
set nolabel 20
set nolabel 50
set nolabel 10000

plot [:8][-0.15:0.95] \
	"N-5e5_T-300_alpha-1_beta-1.0.dat"    using 2:6 notitle  with lines 1, \
	"N-5e5_T-300_alpha-0.5_beta-1.0.dat"  using 2:6 notitle  with lines 1, \
	"N-5e5_T-300_alpha-0.2_beta-1.0.dat"  using 2:6 notitle  with lines 1, \
	"N-5e5_T-300_alpha-0.1_beta-1.0.dat"  using 2:6 notitle  with lines 1, \
	sqrt(2/pi) title 'u_d({/Symbol \245})' with lines 4

# ------------------------------------------------------------------------
#     +-----+  +-----+
#     |     |  |     |
#     +-----+  +-----+
#     +-----+  +-----+
#     |     |  |  X  |
#     +-----+  +-----+
# 
set origin 0.5,0.0
set size 0.5,0.5
set lmargin 9
set rmargin 0
set key top
#set xlabel "N\n d)"
#set ylabel "\n\n{/Symbol e}"
set xlabel "{/Symbol t}\n d)"
set ylabel "u_d({/Symbol t})"

set nolabel

plot [:8][:0.9] \
	"N-5e5_T-300_alpha-0.001_beta-1.0.dat" using 2:6 title "{/Symbol a} = 0.001" with lines 2



set nomultiplot

#pause -1
#"{/Symbol a} = 1.0"
#"{/Symbol a} = 0.5"
#"{/Symbol a} = 0.2"
#"{/Symbol a} = 0.1"
#	sqrt(2/pi) notitle 'u_d({/Symbol \245})' with lines 4