#!/usr/bin/gnuplot

reset

set terminal postscript enhanced
set output "hotGas-drift-speed.ps"
#set terminal postscript eps enhanced 16
#set output "hotGas-drift-speed.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 7
set rmargin 1


set mytics 2
set mxtics 2

set format y "%2.2f"

lambda = 3e-7
#acc    = 1e7
#u_drift(x,y) = x / y / sqrt(lambda * acc)

# -----------------------------------------------------------------------------
# T = 10K
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 title "T = 10K"

plot [:][] \
	"N-1e5_T-10_alpha-50_beta-1.dat"	using 2:6 title "{/Symbol a} = 50 " with dots, \
	"N-1e5_T-10_alpha-5_beta-1.dat"		using 2:6 title "{/Symbol a} = 5.0" with dots, \
	"N-1e5_T-10_alpha-2_beta-1.dat"		using 2:6 title "{/Symbol a} = 2.0" with dots, \
	"N-1e5_T-10_alpha-1_beta-1.dat"		using 2:6 title "{/Symbol a} = 1.0" with dots, \
	"N-1e5_T-10_alpha-0.5_beta-1.dat"	using 2:6 title "{/Symbol a} = 0.5" with dots, \
	"N-1e5_T-10_alpha-0.1_beta-1.dat"	using 2:6 title "{/Symbol a} = 0.1" with dots

# -----------------------------------------------------------------------------
# T = 300K
set origin 0.0,0.0
set size 0.5,0.5
set key top
#set ytics 0.005
set xlabel "{/Symbol t}\n b)"
set ylabel "u_d({/Symbol t})"
set title "T = 300K"

plot [][] \
	"N-1e5_T-300_alpha-100_beta-1.dat"	using 2:6 title "{/Symbol a} = 100" with dots, \
	"N-1e5_T-300_alpha-50_beta-1.dat"	using 2:6 title "{/Symbol a} = 50 " with dots, \
	"N-1e5_T-300_alpha-5_beta-1.dat"	using 2:6 title "{/Symbol a} = 5.0" with dots, \
	"N-1e5_T-300_alpha-2_beta-1.dat"	using 2:6 title "{/Symbol a} = 2.0" with dots, \
	"N-1e5_T-300_alpha-1_beta-1.dat"	using 2:6 title "{/Symbol a} = 1.0" with dots, \
	"N-1e5_T-300_alpha-0.5_beta-1.dat"	using 2:6 title "{/Symbol a} = 0.5" with dots, \
	"N-1e5_T-300_alpha-0.1_beta-1.dat"	using 2:6 title "{/Symbol a} = 0.1" with dots



# -----------------------------------------------------------------------------
# T = 100K
set origin 0.5,0.5
set size 0.5,0.5
set lmargin 9
set rmargin 0
set key bottom
#set ytics 0.001

set xlabel "{/Symbol t}\n c)"
set ylabel "u_d({/Symbol t})"
set title "T = 100K"

plot [][] \
	"N-1e5_T-100_alpha-100_beta-1.dat"	using 2:6 title "{/Symbol a} = 100" with dots, \
	"N-1e5_T-100_alpha-50_beta-1.dat"	using 2:6 title "{/Symbol a} = 50 " with dots, \
	"N-1e5_T-100_alpha-5_beta-1.dat"	using 2:6 title "{/Symbol a} = 5.0" with dots, \
	"N-1e5_T-100_alpha-2_beta-1.dat"	using 2:6 title "{/Symbol a} = 2.0" with dots, \
	"N-1e5_T-100_alpha-1_beta-1.dat"	using 2:6 title "{/Symbol a} = 1.0" with dots, \
	"N-1e5_T-100_alpha-0.5_beta-1.dat"	using 2:6 title "{/Symbol a} = 0.5" with dots, \
	"N-1e5_T-100_alpha-0.1_beta-1.dat"	using 2:6 title "{/Symbol a} = 0.1" with dots

# -----------------------------------------------------------------------------
set origin 0.5,0.0
set size 0.5,0.5
set lmargin 9
set rmargin 0



set nomultiplot

#pause -1

#set label  "{/Symbol \341z\361 / t}" at 3,0.6
#set ytics ("0.000" 0, "0.200" 0.2, "0.400" 0.4, "0.600" 0.6, "0.800" 0.8)
#set ytics ("0.795" 0.795, "0.800" 0.800, "0.805" 0.805, "0.810" 0.810, "0.815" 0.815, "0.820" 0.820)
#set ytics ("0.795" 0.795, "0.796" 0.796, "0.797" 0.797, "0.798" 0.798, "0.799" 0.799, "0.800" 0.800, "0.801" 0.801)
#	"analytic.dat" title "teoret." with lines 2, \
#	sqrt(2/pi) title 'u_d({/Symbol \245})' with lines 4, \
#	"N-1e6-a-1E7-T-0.dat" using 2:(u_drift($5, $1)) notitle with lines 1
#	#sqrt(2/pi) title '{/Symbol \326 (2 / {/Symbol p})}' with lines 4, \