#!/usr/local/bin/gnuplot -persist


set terminal postscript portrait enhanced color dashed lw 1.5 font 'Arial-Bold, 20'  
set output "volume.ps"
set size 2,2
set origin 0,0
set key autotitle columnhead font 'Arial-Bold,20' out horiz top 
set encoding iso_8859_1
set ylabel '{\305}^3/H_{2}O' font 'Arial-Bold, 20' offset +1.5


set multiplot layout 2,2
set label 1 '(a)' at screen 0.75,1.22 font 'Arial-Bold 20'
set label 2 '(b)' at screen 1.52,1.22 font 'Arial-Bold 20'


set key autotitle columnhead font 'Arial-Bold,20' out horiz top 
set title "Without ZPE" font 'Arial-Bold, 20' 
set xrange [ -0.5 : 6.5 ]
#set yrange [ 15 : 35 ]
set size ratio 1.2
set encoding iso_8859_1
set ylabel '(V_{0}^{calc} - V_{0}^{exp})/V_{0}^{exp}' font 'Arial-Bold, 20' offset +1.5
#set ylabel '{\305}^3/H_{2}O' font 'Arial-Bold, 20' offset +1.5
set tics font 'Arial-Bold,20'
set border lw 1.5
unset key
unset logscale y
set grid
set grid lc rgb '#D3D3D3  lw 0.5
set size 0.9,2
set origin 0,0.0


plot "Figure_2_2.txt" using ($2):xticlabels(1) with linespoints  lw 2 ps 1.5 dt 2 pt 1, "" using ($3):xticlabels(1) with linespoints  lw 2 ps 1.5 dt 3 pt 2, "" using ($4):xticlabels(1) with linespoints  lw 2 ps 1.5 dt 4 pt 3, "" using ($5):xticlabels(1) with linespoints  lw 2 ps 1.5 dt 5 pt 4, "" using ($6):xticlabels(1) with linespoints  lw 2 ps 1.5 dt 7 pt 6, "" using ($7):xticlabels(1) with points  lw 2 ps 2.5 dt 8 pt 5, "" using ($8):xticlabels(1) with linespoints  lw 2 ps 1.5 dt 1 pt 7, "" using ($9):xticlabels(1) with linespoints  lw 2 ps 1.5 dt 1 pt 8, "" using ($10):xticlabels(1)  with linespoints lw 2 ps 1.5 dt 1 pt 9, "" using ($11):xticlabels(1)  with linespoints  lw 2 ps 1.5 dt 1 pt 10, "" using ($12):xticlabels(1)  with linespoints  lw 2 ps 1.5 dt 1 pt 11, "" using ($13):xticlabels(1)  with linespoints  lw 2 ps 1.5 dt 1 pt 12, 0 with lines lw 2 lc rgb "black" dt 2



set key autotitle columnhead font 'Arial-Bold,20' out horiz top 
set title "Without ZPE" font 'Arial-Bold, 20' 
set xrange [ -0.5 : 6.5 ]
#set yrange [ 15 : 35 ]
set size ratio 1.2
set ylabel '(V_{0}^{calc} - V_{0}^{exp})/V_{0}^{exp}' font 'Arial-Bold, 20' offset +1.5
#set ylabel '{\305}^3/H_{2}O' font 'Arial-Bold, 20' offset +1.5
set tics font 'Arial-Bold,20'
set border lw 1.5
unset key
unset logscale y
set grid
set grid lc rgb '#D3D3D3  lw 0.5
set size 0.9,2
set origin 0.85,0.0

set title "With ZPE" font 'Arial-Bold, 20' 
plot "Figure_2_2_zpe.txt" using ($2):xticlabels(1) with linespoints  lw 2 ps 1.5 dt 2 pt 1, "" using ($3):xticlabels(1) with linespoints  lw 2 ps 1.5 dt 3 pt 2, "" using ($4):xticlabels(1) with linespoints  lw 2 ps 1.5 dt 4 pt 3, "" using ($5):xticlabels(1) with linespoints  lw 2 ps 1.5 dt 5 pt 4, "" using ($6):xticlabels(1) with linespoints  lw 2 ps 1.5 dt 7 pt 6, "" using ($7):xticlabels(1) with points  lw 2 ps 2.5 dt 8 pt 5, "" using ($8):xticlabels(1) with linespoints  lw 2 ps 1.5 dt 1 pt 7, "" using ($9):xticlabels(1) with linespoints  lw 2 ps 1.5 dt 1 pt 8, "" using ($10):xticlabels(1)  with linespoints lw 2 ps 1.5 dt 1 pt 9, "" using ($11):xticlabels(1)  with linespoints  lw 2 ps 1.5 dt 1 pt 10, "" using ($12):xticlabels(1)  with linespoints  lw 2 ps 1.5 dt 1 pt 11, "" using ($13):xticlabels(1)  with linespoints  lw 2 ps 1.5 dt 1 pt 12, 0 with lines lw 2 lc rgb "black" dt 2

set size 2,0.5
set origin 0.45,0.12
set key outside horizontal 
set border 0 
unset title
unset tics
unset xlabel
unset ylabel
set yrange [0:0.001]
set xrange [0:0.001]

plot "Figure_2_2.txt" using ($2):xticlabels(1) with linespoints  lw 2 ps 1.5 dt 2 pt 1 lc 1, "" using ($3):xticlabels(1) with linespoints  lw 2 ps 1.5 dt 3 pt 2 lc 2, "" using ($4):xticlabels(1) with linespoints  lw 2 ps 1.5 dt 4 pt 3 lc 3, "" using ($5):xticlabels(1) with linespoints  lw 2 ps 1.5 dt 5 pt 4 lc 4, "" using ($6):xticlabels(1) with linespoints  lw 2 ps 1.5 dt 7 pt 6 lc 5, "" using ($8):xticlabels(1) with linespoints  lw 2 ps 1.5 dt 1 pt 7 lc 7, "" using ($9):xticlabels(1) with linespoints  lw 2 ps 1.5 dt 1 pt 8 lc 8, "" using ($10):xticlabels(1)  with linespoints  lw 2 ps 1.5 dt 1 pt 9 lc 9, "" using ($11):xticlabels(1)  with linespoints  lw 2 ps 1.5 dt 1 pt 10 lc 10, "" using ($12):xticlabels(1)  with linespoints  lw 2 ps 1.5 dt 1 pt 11 lc 11, "" using ($13):xticlabels(1)  with linespoints  lw 2 ps 1.5 dt 1 pt 12 lc 12



unset multiplot