#! /opt/local/bin/gnuplot
set term postscript color enhanced lw 1 "Times-New-Roman,16"
set output "vac_True.ps"
set border 3 lw 0.5
set style data points
# set bars small
set samples 1000
LW_NUM = 1.
PS_NUM = 1.
BIG_NUM = 10.

set style line 1 lt 1 lc rgb "red" pt 7 lw LW_NUM ps PS_NUM
set style line 2 lt 1 lc rgb "orange" pt 7 lw LW_NUM ps PS_NUM
set style line 3 lt 1 lc rgb "yellow" pt 7 lw LW_NUM ps PS_NUM
set style line 4 lt 1 lc rgb "green" pt 7 lw LW_NUM ps PS_NUM
set style line 5 lt 1 lc rgb "blue" pt 7 lw LW_NUM ps PS_NUM
set style line 6 lt 1 lc rgb "purple" pt 7 lw LW_NUM ps PS_NUM

X1 = 0.
X2 = 0.25
X3 = 0.44
X4 = 0.63
X5 = 0.82
Y3 = 0.71
Y2 = 0.35
Y1 = 0.

Y_MIN_IG = 1e-05
Y_MIN_WA = 1e-05
Y_MIN_LJ = 1e-05
Y_MAX_IG = 1e00
Y_MAX_WA = 5e08
Y_MAX_LJ = 1e10
Y_GAMMA = 1e-13
X_TYPE = -150
X_GAMMA = 40.

LEFTWIDTH = 0.3
MIDWIDTH = 0.23
RIGHTWIDTH = 0.23
TOPHEIGHT = 0.34
MIDHEIGHT = 0.34
BOTTOMHEIGHT = 0.39

set multiplot

set size LEFTWIDTH,BOTTOMHEIGHT
set origin X1,Y1
# set key top right
unset key
set logscale y
set xtics nomirror 50
set ytics nomirror
set xlabel '{/Symbol D}{/Times-Italic t}, ps' font "Times-New-Roman,14"
set ylabel '<v(t) v(t+{/Symbol D}t)>, (nm/ps)^2' font "Times-New-Roman,14"
set label '{/Symbol g}=0.01' at X_GAMMA,Y_GAMMA
set label 'LJ Fluid' font "Times-New-Roman,18" at X_TYPE,1e04 rotate by 90
plot[0:100] [Y_MIN_LJ:Y_MAX_LJ] 'analysis/vvvr-True-LennardJonesFluid-0.5-0.01.vac' title '{/Symbol D}t=0.5' ls 1 ps 2.0, 'analysis/vvvr-True-LennardJonesFluid-1.0-0.01.vac' title '{/Symbol D}t=1' ls 2 ps 1.8, 'analysis/vvvr-True-LennardJonesFluid-2.0-0.01.vac' title '{/Symbol D}t=2' ls 3 ps 1.6, 'analysis/vvvr-True-LennardJonesFluid-4.0-0.01.vac' title '{/Symbol D}t=4' ls 4 ps 1.4, 'analysis/vvvr-True-LennardJonesFluid-8.0-0.01.vac' title '{/Symbol D}t=8' ls 5 ps 1.2, 'analysis/vvvr-True-LennardJonesFluid-16.0-0.01.vac' title '{/Symbol D}t=16' ls 6

set size MIDWIDTH,BOTTOMHEIGHT
set origin X2,Y1
unset key
unset ylabel
unset label
unset ytics
set label '{/Symbol g}=0.1' at X_GAMMA,Y_GAMMA
plot[0.5:100] [Y_MIN_LJ:Y_MAX_LJ] 'analysis/vvvr-True-LennardJonesFluid-0.5-0.1.vac' title '{/Symbol D}t=0.5' ls 1 ps 2.0, 'analysis/vvvr-True-LennardJonesFluid-1.0-0.1.vac' title '{/Symbol D}t=1' ls 2 ps 1.8, 'analysis/vvvr-True-LennardJonesFluid-2.0-0.1.vac' title '{/Symbol D}t=2' ls 3 ps 1.6, 'analysis/vvvr-True-LennardJonesFluid-4.0-0.1.vac' title '{/Symbol D}t=4' ls 4 ps 1.4, 'analysis/vvvr-True-LennardJonesFluid-8.0-0.1.vac' title '{/Symbol D}t=8' ls 5 ps 1.2, 'analysis/vvvr-True-LennardJonesFluid-16.0-0.1.vac' title '{/Symbol D}t=16' ls 6

set origin X3,Y1
unset label
set label '{/Symbol g}=1' at X_GAMMA,Y_GAMMA
plot[0.5:100] [Y_MIN_LJ:Y_MAX_LJ] 'analysis/vvvr-True-LennardJonesFluid-0.5-1.0.vac' title '{/Symbol D}t=0.5' ls 1 ps 2.0, 'analysis/vvvr-True-LennardJonesFluid-1.0-1.0.vac' title '{/Symbol D}t=1' ls 2 ps 1.8, 'analysis/vvvr-True-LennardJonesFluid-2.0-1.0.vac' title '{/Symbol D}t=2' ls 3 ps 1.6, 'analysis/vvvr-True-LennardJonesFluid-4.0-1.0.vac' title '{/Symbol D}t=4' ls 4 ps 1.4, 'analysis/vvvr-True-LennardJonesFluid-8.0-1.0.vac' title '{/Symbol D}t=8' ls 5 ps 1.2, 'analysis/vvvr-True-LennardJonesFluid-16.0-1.0.vac' title '{/Symbol D}t=16' ls 6

set origin X4,Y1
unset label
set label '{/Symbol g}=10' at X_GAMMA,Y_GAMMA
plot[0.5:100] [Y_MIN_LJ:Y_MAX_LJ] 'analysis/vvvr-True-LennardJonesFluid-0.5-10.0.vac' title '{/Symbol D}t=0.5' ls 1 ps 2.0, 'analysis/vvvr-True-LennardJonesFluid-1.0-10.0.vac' title '{/Symbol D}t=1' ls 2 ps 1.8, 'analysis/vvvr-True-LennardJonesFluid-2.0-10.0.vac' title '{/Symbol D}t=2' ls 3 ps 1.6, 'analysis/vvvr-True-LennardJonesFluid-4.0-10.0.vac' title '{/Symbol D}t=4' ls 4 ps 1.4, 'analysis/vvvr-True-LennardJonesFluid-8.0-10.0.vac' title '{/Symbol D}t=8' ls 5 ps 1.2, 'analysis/vvvr-True-LennardJonesFluid-16.0-10.0.vac' title '{/Symbol D}t=16' ls 6

set origin X5,Y1
unset label
set label '{/Symbol g}=100' at X_GAMMA,Y_GAMMA
# plot[0.5:100] [Y_MIN_LJ:Y_MAX_LJ] 'analysis/vvvr-True-LennardJonesFluid-0.5-100.0,msd' title '{/Symbol D}t=0.5' ls 1 ps 2.0, 'analysis/vvvr-True-LennardJonesFluid-1.0-100.0,msd' title '{/Symbol D}t=1' ls 2 ps 1.8, 'analysis/vvvr-True-LennardJonesFluid-2.0-100.0,msd' title '{/Symbol D}t=2' ls 3 ps 1.6, 'analysis/vvvr-True-LennardJonesFluid-4.0-100.0,msd' title '{/Symbol D}t=4' ls 4 ps 1.4, 'analysis/vvvr-True-LennardJonesFluid-8.0-100.0,msd' title '{/Symbol D}t=8' ls 5 ps 1.2, 'analysis/vvvr-True-LennardJonesFluid-16.0-100.0,msd' title '{/Symbol D}t=16' ls 6


set size LEFTWIDTH,MIDHEIGHT
set origin X1,Y2
unset xlabel
unset xtics
set ytics
set ylabel '<[r(t+{/Symbol D}t) r(t)]^2>, nm^2/ps' font "Times-New-Roman,14"
unset label
set label 'Water Box' font "Times-New-Roman,18" at 2.*X_TYPE,1e04 rotate by 90
# plot [0.5:100] [Y_MIN_WA:Y_MAX_WA] 'analysis/vvvr-True-WaterBox-0.5-0.01.vac' title '{/Symbol D}t=0.5' ls 1 ps 2.0, 'analysis/vvvr-True-WaterBox-1.0-0.01.vac' title '{/Symbol D}t=1' ls 2 ps 1.8, 'analysis/vvvr-True-WaterBox-2.0-0.01.vac' title '{/Symbol D}t=2' ls 3 ps 1.6, 'analysis/vvvr-True-WaterBox-4.0-0.01.vac' title '{/Symbol D}t=4' ls 4 ps 1.4, 'analysis/vvvr-True-WaterBox-8.0-0.01.vac' title '{/Symbol D}t=8' ls 5 ps 1.2, 'analysis/vvvr-True-WaterBox-16.0-0.01.vac' title '{/Symbol D}t=16' ls 6

set size MIDWIDTH,MIDHEIGHT
set origin X2,Y2
unset key
unset ylabel
unset ytics
unset label
# plot[0.5:100] [Y_MIN_WA:Y_MAX_WA] 'analysis/vvvr-True-WaterBox-0.5-0.1.vac' title '{/Symbol D}t=0.5' ls 1 ps 2.0, 'analysis/vvvr-True-WaterBox-1.0-0.1.vac' title '{/Symbol D}t=1' ls 2 ps 1.8, 'analysis/vvvr-True-WaterBox-2.0-0.1.vac' title '{/Symbol D}t=2' ls 3 ps 1.6, 'analysis/vvvr-True-WaterBox-4.0-0.1.vac' title '{/Symbol D}t=4' ls 4 ps 1.4, 'analysis/vvvr-True-WaterBox-8.0-0.1.vac' title '{/Symbol D}t=8' ls 5 ps 1.2, 'analysis/vvvr-True-WaterBox-16.0-0.1.vac' title '{/Symbol D}t=16' ls 6

set origin X3,Y2
# plot[0.5:100] [Y_MIN_WA:Y_MAX_WA] 'analysis/vvvr-True-WaterBox-0.5-1.0.vac' title '{/Symbol D}t=0.5' ls 1 ps 2.0, 'analysis/vvvr-True-WaterBox-1.0-1.0.vac' title '{/Symbol D}t=1' ls 2 ps 1.8, 'analysis/vvvr-True-WaterBox-2.0-1.0.vac' title '{/Symbol D}t=2' ls 3 ps 1.6, 'analysis/vvvr-True-WaterBox-4.0-1.0.vac' title '{/Symbol D}t=4' ls 4 ps 1.4, 'analysis/vvvr-True-WaterBox-8.0-1.0.vac' title '{/Symbol D}t=8' ls 5 ps 1.2, 'analysis/vvvr-True-WaterBox-16.0-1.0.vac' title '{/Symbol D}t=16' ls 6

set origin X4,Y2
# plot[0.5:100] [Y_MIN_WA:Y_MAX_WA] 'analysis/vvvr-True-WaterBox-0.5-10.0.vac' title '{/Symbol D}t=0.5' ls 1 ps 2.0, 'analysis/vvvr-True-WaterBox-1.0-10.0.vac' title '{/Symbol D}t=1' ls 2 ps 1.8, 'analysis/vvvr-True-WaterBox-2.0-10.0.vac' title '{/Symbol D}t=2' ls 3 ps 1.6, 'analysis/vvvr-True-WaterBox-4.0-10.0.vac' title '{/Symbol D}t=4' ls 4 ps 1.4, 'analysis/vvvr-True-WaterBox-8.0-10.0.vac' title '{/Symbol D}t=8' ls 5 ps 1.2, 'analysis/vvvr-True-WaterBox-16.0-10.0.vac' title '{/Symbol D}t=16' ls 6

set origin X5,Y2
# plot[0.5:100] [Y_MIN_WA:Y_MAX_WA] 'analysis/vvvr-True-WaterBox-0.5-100.0.vac' title '{/Symbol D}t=0.5' ls 1 ps 2.0, 'analysis/vvvr-True-WaterBox-1.0-100.0.vac' title '{/Symbol D}t=1' ls 2 ps 1.8, 'analysis/vvvr-True-WaterBox-2.0-100.0.vac' title '{/Symbol D}t=2' ls 3 ps 1.6, 'analysis/vvvr-True-WaterBox-4.0-100.0.vac' title '{/Symbol D}t=4' ls 4 ps 1.4, 'analysis/vvvr-True-WaterBox-8.0-100.0.vac' title '{/Symbol D}t=8' ls 5 ps 1.2, 'analysis/vvvr-True-WaterBox-16.0-100.0.vac' title '{/Symbol D}t=16' ls 6


set size LEFTWIDTH,TOPHEIGHT
set origin X1,Y3
unset xlabel
set ytics
set ylabel '<[r(t+{/Symbol D}t) r(t)]^2>, nm^2/ps' font "Times-New-Roman,14"
set label 'Ideal Gas' font "Times-New-Roman,18" at 10.*X_TYPE,1e-03 rotate by 90
plot[0.5:100] [Y_MIN_IG:Y_MAX_IG] 'analysis/vvvr-True-IdealGas-0.5-0.01.vac' title '{/Symbol D}t=0.5' ls 1 ps 2.0, 'analysis/vvvr-True-IdealGas-1.0-0.01.vac' title '{/Symbol D}t=1' ls 2 ps 1.8, 'analysis/vvvr-True-IdealGas-2.0-0.01.vac' title '{/Symbol D}t=2' ls 3 ps 1.6, 'analysis/vvvr-True-IdealGas-4.0-0.01.vac' title '{/Symbol D}t=4' ls 4 ps 1.4, 'analysis/vvvr-True-IdealGas-8.0-0.01.vac' title '{/Symbol D}t=8' ls 5 ps 1.2, 'analysis/vvvr-True-IdealGas-16.0-0.01.vac' title '{/Symbol D}t=16' ls 6

set size MIDWIDTH,TOPHEIGHT
set origin X2,Y3
unset key
unset ylabel
unset ytics
unset label
plot[0.5:100] [Y_MIN_IG:Y_MAX_IG] 'analysis/vvvr-True-IdealGas-0.5-0.1.vac' title '{/Symbol D}t=0.5' ls 1 ps 2.0, 'analysis/vvvr-True-IdealGas-1.0-0.1.vac' title '{/Symbol D}t=1' ls 2 ps 1.8, 'analysis/vvvr-True-IdealGas-2.0-0.1.vac' title '{/Symbol D}t=2' ls 3 ps 1.6, 'analysis/vvvr-True-IdealGas-4.0-0.1.vac' title '{/Symbol D}t=4' ls 4 ps 1.4, 'analysis/vvvr-True-IdealGas-8.0-0.1.vac' title '{/Symbol D}t=8' ls 5 ps 1.2, 'analysis/vvvr-True-IdealGas-16.0-0.1.vac' title '{/Symbol D}t=16' ls 6

set origin X3,Y3
plot[0.5:100] [Y_MIN_IG:Y_MAX_IG] 'analysis/vvvr-True-IdealGas-0.5-1.0.vac' title '{/Symbol D}t=0.5' ls 1 ps 2.0, 'analysis/vvvr-True-IdealGas-1.0-1.0.vac' title '{/Symbol D}t=1' ls 2 ps 1.8, 'analysis/vvvr-True-IdealGas-2.0-1.0.vac' title '{/Symbol D}t=2' ls 3 ps 1.6, 'analysis/vvvr-True-IdealGas-4.0-1.0.vac' title '{/Symbol D}t=4' ls 4 ps 1.4, 'analysis/vvvr-True-IdealGas-8.0-1.0.vac' title '{/Symbol D}t=8' ls 5 ps 1.2, 'analysis/vvvr-True-IdealGas-16.0-1.0.vac' title '{/Symbol D}t=16' ls 6

set origin X4,Y3
plot[0.5:100] [Y_MIN_IG:Y_MAX_IG] 'analysis/vvvr-True-IdealGas-0.5-10.0.vac' title '{/Symbol D}t=0.5' ls 1 ps 2.0, 'analysis/vvvr-True-IdealGas-1.0-10.0.vac' title '{/Symbol D}t=1' ls 2 ps 1.8, 'analysis/vvvr-True-IdealGas-2.0-10.0.vac' title '{/Symbol D}t=2' ls 3 ps 1.6, 'analysis/vvvr-True-IdealGas-4.0-10.0.vac' title '{/Symbol D}t=4' ls 4 ps 1.4, 'analysis/vvvr-True-IdealGas-8.0-10.0.vac' title '{/Symbol D}t=8' ls 5 ps 1.2, 'analysis/vvvr-True-IdealGas-16.0-10.0.vac' title '{/Symbol D}t=16' ls 6

set origin X5,Y3
plot[0.5:100] [Y_MIN_IG:Y_MAX_IG] 'analysis/vvvr-True-IdealGas-0.5-100.0.vac' title '{/Symbol D}t=0.5' ls 1 ps 2.0, 'analysis/vvvr-True-IdealGas-1.0-100.0.vac' title '{/Symbol D}t=1' ls 2 ps 1.8, 'analysis/vvvr-True-IdealGas-2.0-100.0.vac' title '{/Symbol D}t=2' ls 3 ps 1.6, 'analysis/vvvr-True-IdealGas-4.0-100.0.vac' title '{/Symbol D}t=4' ls 4 ps 1.4, 'analysis/vvvr-True-IdealGas-8.0-100.0.vac' title '{/Symbol D}t=8' ls 5 ps 1.2, 'analysis/vvvr-True-IdealGas-16.0-100.0.vac' title '{/Symbol D}t=16' ls 6