#!/usr/bin/gnuplot

reset

#set terminal postscript enhanced
#set output "vdrift_T-300K-beta-model.ps"
set terminal postscript eps enhanced 16
set output "vdrift_T-300K-beta-model.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.1f"

set grid

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

set title "{/Symbol b} = 0.1"
set xlabel "{/Symbol t}\n a)"
set ylabel "u_d({/Symbol t})"


set label 1    "{/Symbol a} = 0.01" at 6,0.13
set label 2    "{/Symbol a} = 0.02" at 6,0.25
set label 5    "{/Symbol a} = 0.05" at 6,0.38
set label 10   "{/Symbol a} = 0.10" at 6,0.5
set label 50   "{/Symbol a} = 0.50" at 6,0.60
set label 100  "{/Symbol a} = 1.00" at 6,0.75

plot [:8][-0.15:0.95] \
	"N-5e5_T-300_alpha-1_beta-0.1_speed-0.dat" using 2:6 notitle  with lines 1, \
	"N-5e5_T-300_alpha-0.5_beta-0.1_speed-0.dat"    using 2:6 notitle  with lines 1, \
	"N-5e5_T-300_alpha-0.1_beta-0.1_speed-0.dat"    using 2:6 notitle  with lines 1, \
	"N-5e5_T-300_alpha-0.05_beta-0.1_speed-0.dat"    using 2:6 notitle  with lines 1, \
	"N-5e5_T-300_alpha-0.02_beta-0.1_speed-0.dat"    using 2:6 notitle  with lines 1, \
	"N-5e5_T-300_alpha-0.01_beta-0.1_speed-0.dat" using 2:6 notitle with lines 1

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

set title "{/Symbol b} = 0.5"
set xlabel "{/Symbol t}\n c)"
set ylabel "u_d({/Symbol t})"

#set label 1    "{/Symbol a} = 0.01" at 6,-0.05
#set label 5    "{/Symbol a} = 0.05" at 6,0.20
set nolabel 1
set nolabel 2
set label 5    "{/Symbol a \243} 0.05" at 6.2,0.00
set label 10   "{/Symbol a} = 0.10" at 6,0.20
set label 50   "{/Symbol a} = 0.50" at 6,0.45
set label 100  "{/Symbol a} = 1.00" at 6,0.60

plot [:8][-0.15:0.95] \
	"N-5e5_T-300_alpha-1_beta-0.5_speed-0.dat" using 2:6 notitle  with lines 1, \
	"N-5e5_T-300_alpha-0.5_beta-0.5_speed-0.dat"    using 2:6 notitle  with lines 1, \
	"N-5e5_T-300_alpha-0.1_beta-0.5_speed-0.dat"    using 2:6 notitle  with lines 1, \
	"N-5e5_T-300_alpha-0.05_beta-0.5_speed-0.dat"    using 2:6 notitle  with lines 1, \
	"N-5e5_T-300_alpha-0.02_beta-0.5_speed-0.dat"    using 2:6 notitle  with lines 1, \
	"N-5e5_T-300_alpha-0.01_beta-0.5_speed-0.dat" using 2:6 notitle with lines 1

# ------------------------------------------------------------------------
#     +-----+  +-----+
#     |     |  |  X  |
#     +-----+  +-----+
#     +-----+  +-----+
#     |     |  |     |
#     +-----+  +-----+
set origin 0.5,0.5
set size 0.5,0.5
set lmargin 9
set rmargin 0
set key top

set title "{/Symbol b} = 0.25"
set xlabel "{/Symbol t}\n b)"
set ylabel "u_d({/Symbol t})"

set label 1    "{/Symbol a} = 0.01" at 6,-0.05
set label 2    "{/Symbol a} = 0.02" at 6,0.10
set label 5    "{/Symbol a} = 0.05" at 6,0.22
set label 10   "{/Symbol a} = 0.10" at 6,0.35
set label 50   "{/Symbol a} = 0.50" at 6,0.60
set label 100  "{/Symbol a} = 1.00" at 6,0.70

plot [:8][-0.15:0.95] \
	"N-5e5_T-300_alpha-1_beta-0.25_speed-0.dat"    using 2:6 notitle  with lines 1, \
	"N-5e5_T-300_alpha-0.5_beta-0.25_speed-0.dat"  using 2:6 notitle  with lines 1, \
	"N-5e5_T-300_alpha-0.1_beta-0.25_speed-0.dat"  using 2:6 notitle  with lines 1, \
	"N-5e5_T-300_alpha-0.05_beta-0.25_speed-0.dat"  using 2:6 notitle  with lines 1, \
	"N-5e5_T-300_alpha-0.02_beta-0.25_speed-0.dat"  using 2:6 notitle  with lines 1, \
	"N-5e5_T-300_alpha-0.01_beta-0.25_speed-0.dat"  using 2:6 notitle  with lines 1

# ------------------------------------------------------------------------
#     +-----+  +-----+
#     |     |  |     |
#     +-----+  +-----+
#     +-----+  +-----+
#     |     |  |  X  |
#     +-----+  +-----+

set origin 0.5,0.0
set size 0.5,0.5
set lmargin 9
set rmargin 0
set key top

set title "{/Symbol b} = 1.0"
set xlabel "{/Symbol t}\n d)"
set ylabel "u_d({/Symbol t})"

#set label 1    "{/Symbol a} = 0.01" at 6,-0.05
#set label 5    "{/Symbol a} = 0.05" at 6,0.20
#set label 10   "{/Symbol a} = 0.10" at 6,0.85
set nolabel 1
set nolabel 2
set nolabel 10
set label 5    "{/Symbol a \243} 0.10" at 6.2,0.00
set label 50   "{/Symbol a} = 0.50" at 6,0.3
set label 100  "{/Symbol a} = 1.00" at 6,0.50

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

set nomultiplot