#! /opt/local/bin/gnuplot
set term postscript color enhanced lw 1 "Times-New-Roman,16"
set output "msv.ps"
set border 0 lw 0.5
set style data errorbars
# set bars small
unset xlabel
unset ylabel
# set key vertical
set samples 1000
LW_NUM = 1.
PS_NUM = 1.
BIG_NUM = 10.

LEFT_WIDTH = 0.59
RIGHT_WIDTH = 0.47
TOP_HEIGHT = 0.34
BOTTOM_HEIGHT = 0.44
X_L = -0.01
X_R = 0.55
Y_T = 0.7
Y_M = 0.4
Y_B = -0.05

Y_MIN = 5e-02
Y_MAX = 1e13

set style line 1 lt 1 lc rgb "red" pt 7 lw LW_NUM ps 2.0
set style line 11 lt 1 lc rgb "red" pt 5 lw LW_NUM ps 2.0
set style line 2 lt 1 lc rgb "orange" pt 7 lw LW_NUM ps 1.8
set style line 12 lt 1 lc rgb "orange" pt 5 lw LW_NUM ps 1.8
set style line 3 lt 1 lc rgb "yellow" pt 7 lw LW_NUM ps 1.6
set style line 13 lt 1 lc rgb "yellow" pt 5 lw LW_NUM ps 1.6
set style line 4 lt 1 lc rgb "green" pt 7 lw LW_NUM ps 1.4
set style line 14 lt 1 lc rgb "green" pt 5 lw LW_NUM ps 1.4
set style line 5 lt 1 lc rgb "blue" pt 7 lw LW_NUM ps 1.2
set style line 15 lt 1 lc rgb "blue" pt 5 lw LW_NUM ps 1.2
set style line 6 lt 1 lc rgb "purple" pt 7 lw LW_NUM ps 1.0
set style line 16 lt 1 lc rgb "purple" pt 5 lw LW_NUM ps 1.0


set multiplot
set logscale xy

set size LEFT_WIDTH,BOTTOM_HEIGHT
set origin X_L, Y_B
unset key
set ylabel '<v^2>, (nm/ps)^2' font "Times-New-Roman,14"
set ytics 1e03
set xlabel 'Collision rate, ps^{-1}' font "Times-New-Roman,14"
set label 'Timescale Corrected' font "Times-New-Roman,18" at 1e-01,1e-09
set label 'LJ Fluid' font "Times-New-Roman,18" at 3e-05,100. rotate by 90
plot[0.005:2000] [Y_MIN:Y_MAX] 'analysis/vvvr-True-LennardJonesFluid-0.5.out' using 2:3:(2*$4) title '{/Symbol D}t=0.5' ls 1, 'analysis/vvvr-True-LennardJonesFluid-1.0.out' using 2:3:(2*$4) title '{/Symbol D}t=1.0' ls 2, 'analysis/vvvr-True-LennardJonesFluid-2.0.out' using 2:3:(2*$4) title '{/Symbol D}t=2.0' ls 3, 'analysis/vvvr-True-LennardJonesFluid-4.0.out' using 2:3:(2*$4) title '{/Symbol D}t=4.0' ls 4, 'analysis/vvvr-True-LennardJonesFluid-8.0.out' using 2:3:(2*$4) title '{/Symbol D}t=8.0' ls 5, 'analysis/vvvr-True-LennardJonesFluid-16.0.out' using 2:3:(2*$4) title '{/Symbol D}t=16.0' ls 6, 1.87e-01 title '3kT/m' ls -1 w l

set size RIGHT_WIDTH,BOTTOM_HEIGHT
set origin X_R, Y_B
unset key
unset ylabel
unset ytics
set xlabel 'Collision rate, ps^{-1}' font "Times-New-Roman,14"
unset label
set label 'Uncorrected' font "Times-New-Roman,18" at 2e-01,1e-09
plot[0.005:2000] [Y_MIN:Y_MAX] 'analysis/vvvr-False-LennardJonesFluid-0.5.out' using 2:3:(2*$4) title '{/Symbol D}t=0.5' ls 1, 'analysis/vvvr-False-LennardJonesFluid-1.0.out' using 2:3:(2*$4) title '{/Symbol D}t=1.0' ls 2, 'analysis/vvvr-False-LennardJonesFluid-2.0.out' using 2:3:(2*$4) title '{/Symbol D}t=2.0' ls 3, 'analysis/vvvr-False-LennardJonesFluid-4.0.out' using 2:3:(2*$4) title '{/Symbol D}t=4.0' ls 4, 'analysis/vvvr-False-LennardJonesFluid-8.0.out' using 2:3:(2*$4) title '{/Symbol D}t=8.0' ls 5, 'analysis/vvvr-False-LennardJonesFluid-16.0.out' using 2:3:(2*$4) title '{/Symbol D}t=16.0' ls 6, 1.87e-01 title '3kT/m' ls -1 w l

set size LEFT_WIDTH,TOP_HEIGHT
set origin X_L, Y_M
set ytics
set ylabel '<v^2>, (nm/ps)^2' font "Times-New-Roman,14"
unset xlabel
unset xtics
unset label
set label 'Water' font "Times-New-Roman,18" at 3e-05,100. rotate by 90
plot[0.005:2000] [Y_MIN:Y_MAX] 'analysis/vvvr-True-WaterBox-0.5.out' using 2:3:(2*$4) title '{/Symbol D}t=0.5' ls 1, 'analysis/vvvr-True-WaterBox-1.0.out' using 2:3:(2*$4) title '{/Symbol D}t=1.0' ls 2, 'analysis/vvvr-True-WaterBox-2.0.out' using 2:3:(2*$4) title '{/Symbol D}t=2.0' ls 3, 'analysis/vvvr-True-WaterBox-4.0.out' using 2:3:(2*$4) title '{/Symbol D}t=4.0' ls 4, 'analysis/vvvr-True-WaterBox-8.0.out' using 2:3:(2*$4) title '{/Symbol D}t=8.0' ls 5, 'analysis/vvvr-True-WaterBox-16.0.out' using 2:3:(2*$4) title '{/Symbol D}t=16.0' ls 6, 4.67e-01 title '3kT/m' ls -1 w l

set size RIGHT_WIDTH,TOP_HEIGHT
set origin X_R, Y_M
unset ytics
unset ylabel
unset label
plot[0.005:2000] [Y_MIN:Y_MAX] 'analysis/vvvr-False-WaterBox-0.5.out' using 2:3:(2*$4) title '{/Symbol D}t=0.5' ls 1, 'analysis/vvvr-False-WaterBox-1.0.out' using 2:3:(2*$4) title '{/Symbol D}t=1.0' ls 2, 'analysis/vvvr-False-WaterBox-2.0.out' using 2:3:(2*$4) title '{/Symbol D}t=2.0' ls 3, 'analysis/vvvr-False-WaterBox-4.0.out' using 2:3:(2*$4) title '{/Symbol D}t=4.0' ls 4, 'analysis/vvvr-False-WaterBox-8.0.out' using 2:3:(2*$4) title '{/Symbol D}t=8.0' ls 5, 'analysis/vvvr-False-WaterBox-16.0.out' using 2:3:(2*$4) title '{/Symbol D}t=16.0' ls 6, 4.67e-01 title '3kT/m' ls -1 w l

set size LEFT_WIDTH,TOP_HEIGHT
set origin X_L, Y_T
set ytics
set ylabel '<v^2>, (nm/ps)^2' font "Times-New-Roman,14"
unset xlabel
unset xtics
# set xlabel 'Collision rate, ps^{-1}' font "Times-New-Roman,14"
unset label
set label 'Ideal Gas' font "Times-New-Roman,18" at 3e-05,100. rotate by 90
plot[0.005:2000] [Y_MIN:Y_MAX] 'analysis/vvvr-True-IdealGas-0.5.out' using 2:3:(2*$4) title '{/Symbol D}t=0.5' ls 1, 'analysis/vvvr-True-IdealGas-1.0.out' using 2:3:(2*$4) title '{/Symbol D}t=1.0' ls 2, 'analysis/vvvr-True-IdealGas-2.0.out' using 2:3:(2*$4) title '{/Symbol D}t=2.0' ls 3, 'analysis/vvvr-True-IdealGas-4.0.out' using 2:3:(2*$4) title '{/Symbol D}t=4.0' ls 4, 'analysis/vvvr-True-IdealGas-8.0.out' using 2:3:(2*$4) title '{/Symbol D}t=8.0' ls 5, 'analysis/vvvr-True-IdealGas-16.0.out' using 2:3:(2*$4) title '{/Symbol D}t=16.0' ls 6, 1.87e-01 title '3kT/m' ls -1 w l

set size RIGHT_WIDTH,TOP_HEIGHT
set origin X_R, Y_T
unset ytics
unset ylabel
set key top right
unset label
plot[0.005:2000] [Y_MIN:Y_MAX] 'analysis/vvvr-False-IdealGas-0.5.out' using 2:3:(2*$4) title '{/Symbol D}t=0.5' ls 1, 'analysis/vvvr-False-IdealGas-1.0.out' using 2:3:(2*$4) title '{/Symbol D}t=1.0' ls 2, 'analysis/vvvr-False-IdealGas-2.0.out' using 2:3:(2*$4) title '{/Symbol D}t=2.0' ls 3, 'analysis/vvvr-False-IdealGas-4.0.out' using 2:3:(2*$4) title '{/Symbol D}t=4.0' ls 4, 'analysis/vvvr-False-IdealGas-8.0.out' using 2:3:(2*$4) title '{/Symbol D}t=8.0' ls 5, 'analysis/vvvr-False-IdealGas-16.0.out' using 2:3:(2*$4) title '{/Symbol D}t=16.0' ls 6, 1.87e-01 title '3kT/m' ls -1 w l