#!/bin/bash

FUNC="SRP48"
gnuplot<<EOF
set output "${FUNC}.ps"
set terminal postscript enhanced color landscape
set multiplot layout 1,4 title "${FUNC} E_0({/Symbol n},J) \n two scaling methods" font ",12"

set bmargin 1.0
set tmargin 6.5
set lmargin 3.5
set rmargin 0.5

set tics font "Helvetica, 10"

set title "H_2 {/Symbol n} = 0" font ",10"
unset key
set ylabel "E_0({/Symbol n},J) [ev]
set xlabel "J"
set xrange[-0.5:11.5]
set yrange [0.25:0.9]
plot "exp.H2.v0"  lc -1 pointtype 7 pointsize 0.5, "../E0/${FUNC}.H2.v0" lc 1 pointtype 7 pointsize 0.5, "${FUNC}.H2.v0" lc 2 pointtype 7 pointsize 0.5 

#plot 2
unset ylabel
set title "H_2 {/Symbol n} = 1"
set key t
set key font "Helvetica, 7"
set xrange[-0.5:7.5]
set yrange [0.25:0.9]
plot "exp.H2.v1"  lc -1 pointtype 7 pointsize 0.5 title"exp.", "../E0/${FUNC}.H2.v1" lc 1 pointtype 7 pointsize 0.5 title"${FUNC} (a)", "${FUNC}.H2.v1" lc 2 pointtype 7 pointsize 0.5 title"${FUNC} (b)" 

#plot 3
set title"D_2 {/Symbol n} = 0"
unset key 
set xlabel"J"
set xrange [-0.5:11.5]
set yrange [0.25:0.9]
#exp
plot "exp.D2.v0"  lc -1 pointtype 7 pointsize 0.5 , "../E0/${FUNC}.D2.v0" lc 1 pointtype 7 pointsize 0.5, "${FUNC}.D2.v0" lc 2 pointtype 7 pointsize 0.5 

#plot 4
unset key
set title"D_2 {/Symbol n} = 1" 
set xlabel"J"
set xrange [-0.5:7.5]
set yrange [0.25:0.9]
plot "exp.D2.v1"  lc -1 pointtype 7 pointsize 0.5, "../E0/${FUNC}.D2.v1" lc 1 pointtype 7 pointsize 0.5, "${FUNC}.D2.v1" lc 2 pointtype 7 pointsize 0.5
EOF

ps2pdf ${FUNC}.ps
#evince ${FUNC}.pdf

FUNC="B86SRP68-DF2"
gnuplot<<EOF
set output "${FUNC}.ps"
set terminal postscript enhanced color landscape
set multiplot layout 1,4 title "${FUNC} E_0({/Symbol n},J) \n two scaling methods" font ",12"

set bmargin 1.0
set tmargin 6.5
set lmargin 3.5
set rmargin 0.5

set tics font "Helvetica, 10"

set title "H_2 {/Symbol n} = 0" font ",10"
unset key
set ylabel "E_0({/Symbol n},J) [ev]
set xlabel "J"
set xrange[-0.5:11.5]
set yrange [0.25:0.9]
plot "exp.H2.v0"  lc -1 pointtype 7 pointsize 0.5, "../E0/${FUNC}.H2.v0" lc 1 pointtype 7 pointsize 0.5, "${FUNC}.H2.v0" lc 2 pointtype 7 pointsize 0.5, "../E0/${FUNC}.QD.H2.v0" lc 3 pointtype 7 pointsize 0.5, "${FUNC}.QD.H2.v0" lc 4 pointtype 7 pointsize 0.5 

#plot 2
unset ylabel
set title "H_2 {/Symbol n} = 1"
set key t
set key font "Helvetica, 7"
set xrange[-0.5:7.5]
set yrange [0.25:0.9]
plot "exp.H2.v1"  lc -1 pointtype 7 pointsize 0.5 title"exp.", "../E0/${FUNC}.H2.v1" lc 1 pointtype 7 pointsize 0.5 title"${FUNC} (a)", "${FUNC}.H2.v1" lc 2 pointtype 7 pointsize 0.5 title"${FUNC} (b)" , "../E0/${FUNC}.QD.H2.v1" lc 3 pointtype 7 pointsize 0.5 title"${FUNC} QD (a)", "${FUNC}.QD.H2.v1" lc 4 pointtype 7 pointsize 0.5 title"${FUNC} QD (b)"

#plot 3
set title"D_2 {/Symbol n} = 0"
unset key 
set xlabel"J"
set xrange [-0.5:11.5]
set yrange [0.25:0.9]
#exp
plot "exp.D2.v0"  lc -1 pointtype 7 pointsize 0.5 , "../E0/${FUNC}.D2.v0" lc 1 pointtype 7 pointsize 0.5, "${FUNC}.D2.v0" lc 2 pointtype 7 pointsize 0.5 

#plot 4
unset key
set title"D_2 {/Symbol n} = 1" 
set xlabel"J"
set xrange [-0.5:7.5]
set yrange [0.25:0.9]
plot "exp.D2.v1"  lc -1 pointtype 7 pointsize 0.5, "../E0/${FUNC}.D2.v1" lc 1 pointtype 7 pointsize 0.5, "${FUNC}.D2.v1" lc 2 pointtype 7 pointsize 0.5
EOF

ps2pdf ${FUNC}.ps
#evince ${FUNC}.pdf

FUNC="PBEa57-DF2"
gnuplot<<EOF
set output "${FUNC}.ps"
set terminal postscript enhanced color landscape
set multiplot layout 1,4 title "${FUNC} E_0({/Symbol n},J) \n two scaling methods" font ",12"

set bmargin 1.0
set tmargin 6.5
set lmargin 3.5
set rmargin 0.5

set tics font "Helvetica, 10"

set title "H_2 {/Symbol n} = 0" font ",10"
unset key
set ylabel "E_0({/Symbol n},J) [ev]
set xlabel "J"
set xrange[-0.5:11.5]
set yrange [0.25:0.9]
plot "exp.H2.v0"  lc -1 pointtype 7 pointsize 0.5, "../E0/${FUNC}.H2.v0" lc 1 pointtype 7 pointsize 0.5, "${FUNC}.H2.v0" lc 2 pointtype 7 pointsize 0.5 

#plot 2
unset ylabel
set title "H_2 {/Symbol n} = 1"
set key t
set key font "Helvetica, 7"
set xrange[-0.5:7.5]
set yrange [0.25:0.9]
plot "exp.H2.v1"  lc -1 pointtype 7 pointsize 0.5 title"exp.", "../E0/${FUNC}.H2.v1" lc 1 pointtype 7 pointsize 0.5 title"${FUNC} (a)", "${FUNC}.H2.v1" lc 2 pointtype 7 pointsize 0.5 title"${FUNC} (b)" 

#plot 3
set title"D_2 {/Symbol n} = 0"
unset key 
set xlabel"J"
set xrange [-0.5:11.5]
set yrange [0.25:0.9]
#exp
plot "exp.D2.v0"  lc -1 pointtype 7 pointsize 0.5 , "../E0/${FUNC}.D2.v0" lc 1 pointtype 7 pointsize 0.5, "${FUNC}.D2.v0" lc 2 pointtype 7 pointsize 0.5 

#plot 4
unset key
set title"D_2 {/Symbol n} = 1" 
set xlabel"J"
set xrange [-0.5:7.5]
set yrange [0.25:0.9]
plot "exp.D2.v1"  lc -1 pointtype 7 pointsize 0.5, "../E0/${FUNC}.D2.v1" lc 1 pointtype 7 pointsize 0.5, "${FUNC}.D2.v1" lc 2 pointtype 7 pointsize 0.5
EOF

ps2pdf ${FUNC}.ps
#evince ${FUNC}.pdf

FUNC="MS.B86b"
gnuplot<<EOF
set output "${FUNC}.ps"
set terminal postscript enhanced color landscape
set multiplot layout 1,4 title "${FUNC} E_0({/Symbol n},J) \n two scaling methods" font ",12"

set bmargin 1.0
set tmargin 6.5
set lmargin 3.5
set rmargin 0.5

set tics font "Helvetica, 10"

set title "H_2 {/Symbol n} = 0" font ",10"
unset key
set ylabel "E_0({/Symbol n},J) [ev]
set xlabel "J"
set xrange[-0.5:11.5]
set yrange [0.25:0.9]
plot "exp.H2.v0"  lc -1 pointtype 7 pointsize 0.5, "../E0/${FUNC}.H2.v0" lc 1 pointtype 7 pointsize 0.5, "${FUNC}.H2.v0" lc 2 pointtype 7 pointsize 0.5 

#plot 2
unset ylabel
set title "H_2 {/Symbol n} = 1"
set key t
set key font "Helvetica, 7"
set xrange[-0.5:7.5]
set yrange [0.25:0.9]
plot "exp.H2.v1"  lc -1 pointtype 7 pointsize 0.5 title"exp.", "../E0/${FUNC}.H2.v1" lc 1 pointtype 7 pointsize 0.5 title"${FUNC} (a)", "${FUNC}.H2.v1" lc 2 pointtype 7 pointsize 0.5 title"${FUNC} (b)" 

#plot 3
set title"D_2 {/Symbol n} = 0"
unset key 
set xlabel"J"
set xrange [-0.5:11.5]
set yrange [0.25:0.9]
#exp
plot "exp.D2.v0"  lc -1 pointtype 7 pointsize 0.5 , "../E0/${FUNC}.D2.v0" lc 1 pointtype 7 pointsize 0.5, "${FUNC}.D2.v0" lc 2 pointtype 7 pointsize 0.5 

#plot 4
unset key
set title"D_2 {/Symbol n} = 1" 
set xlabel"J"
set xrange [-0.5:7.5]
set yrange [0.25:0.9]
plot "exp.D2.v1"  lc -1 pointtype 7 pointsize 0.5, "../E0/${FUNC}.D2.v1" lc 1 pointtype 7 pointsize 0.5, "${FUNC}.D2.v1" lc 2 pointtype 7 pointsize 0.5
EOF

ps2pdf ${FUNC}.ps
#evince ${FUNC}.pdf

FUNC="SRPsol63-DF2"
gnuplot<<EOF
set output "${FUNC}.ps"
set terminal postscript enhanced color landscape
set multiplot layout 1,4 title "${FUNC} E_0({/Symbol n},J) \n two scaling methods" font ",12"

set bmargin 1.0
set tmargin 6.5
set lmargin 3.5
set rmargin 0.5

set tics font "Helvetica, 10"

set title "H_2 {/Symbol n} = 0" font ",10"
unset key
set ylabel "E_0({/Symbol n},J) [ev]
set xlabel "J"
set xrange[-0.5:11.5]
set yrange [0.25:0.9]
plot "exp.H2.v0"  lc -1 pointtype 7 pointsize 0.5, "../E0/${FUNC}.H2.v0" lc 1 pointtype 7 pointsize 0.5, "${FUNC}.H2.v0" lc 2 pointtype 7 pointsize 0.5 

#plot 2
unset ylabel
set title "H_2 {/Symbol n} = 1"
set key t
set key font "Helvetica, 7"
set xrange[-0.5:7.5]
set yrange [0.25:0.9]
plot "exp.H2.v1"  lc -1 pointtype 7 pointsize 0.5 title"exp.", "../E0/${FUNC}.H2.v1" lc 1 pointtype 7 pointsize 0.5 title"${FUNC} (a)", "${FUNC}.H2.v1" lc 2 pointtype 7 pointsize 0.5 title"${FUNC} (b)" 

#plot 3
set title"D_2 {/Symbol n} = 0"
unset key 
set xlabel"J"
set xrange [-0.5:11.5]
set yrange [0.25:0.9]
#exp
plot "exp.D2.v0"  lc -1 pointtype 7 pointsize 0.5 , "../E0/${FUNC}.D2.v0" lc 1 pointtype 7 pointsize 0.5, "${FUNC}.D2.v0" lc 2 pointtype 7 pointsize 0.5 

#plot 4
unset key
set title"D_2 {/Symbol n} = 1" 
set xlabel"J"
set xrange [-0.5:7.5]
set yrange [0.25:0.9]
plot "exp.D2.v1"  lc -1 pointtype 7 pointsize 0.5, "../E0/${FUNC}.D2.v1" lc 1 pointtype 7 pointsize 0.5, "${FUNC}.D2.v1" lc 2 pointtype 7 pointsize 0.5
EOF

ps2pdf ${FUNC}.ps
#evince ${FUNC}.pdf

FUNC="optPBE-DF1"
gnuplot<<EOF
set output "${FUNC}.ps"
set terminal postscript enhanced color landscape
set multiplot layout 1,4 title "${FUNC} E_0({/Symbol n},J) \n two scaling methods" font ",12"

set bmargin 1.0
set tmargin 6.5
set lmargin 3.5
set rmargin 0.5

set tics font "Helvetica, 10"

set title "H_2 {/Symbol n} = 0" font ",10"
unset key
set ylabel "E_0({/Symbol n},J) [ev]
set xlabel "J"
set xrange[-0.5:11.5]
set yrange [0.25:0.9]
plot "exp.H2.v0"  lc -1 pointtype 7 pointsize 0.5, "../E0/${FUNC}.H2.v0" lc 1 pointtype 7 pointsize 0.5, "${FUNC}.H2.v0" lc 2 pointtype 7 pointsize 0.5 

#plot 2
unset ylabel
set title "H_2 {/Symbol n} = 1"
set key t
set key font "Helvetica, 7"
set xrange[-0.5:7.5]
set yrange [0.25:0.9]
plot "exp.H2.v1"  lc -1 pointtype 7 pointsize 0.5 title"exp.", "../E0/${FUNC}.H2.v1" lc 1 pointtype 7 pointsize 0.5 title"${FUNC} (a)", "${FUNC}.H2.v1" lc 2 pointtype 7 pointsize 0.5 title"${FUNC} (b)" 

#plot 3
set title"D_2 {/Symbol n} = 0"
unset key 
set xlabel"J"
set xrange [-0.5:11.5]
set yrange [0.25:0.9]
#exp
plot "exp.D2.v0"  lc -1 pointtype 7 pointsize 0.5 , "../E0/${FUNC}.D2.v0" lc 1 pointtype 7 pointsize 0.5, "${FUNC}.D2.v0" lc 2 pointtype 7 pointsize 0.5 

#plot 4
unset key
set title"D_2 {/Symbol n} = 1" 
set xlabel"J"
set xrange [-0.5:7.5]
set yrange [0.25:0.9]
plot "exp.D2.v1"  lc -1 pointtype 7 pointsize 0.5, "../E0/${FUNC}.D2.v1" lc 1 pointtype 7 pointsize 0.5, "${FUNC}.D2.v1" lc 2 pointtype 7 pointsize 0.5
EOF

ps2pdf ${FUNC}.ps
#evince ${FUNC}.pdf



FUNC="all_panelB"
gnuplot<<EOF
set output "${FUNC}.ps"
set terminal postscript enhanced color landscape
set multiplot layout 1,4 title "E_0({/Symbol n},J) \n panel B scaling (exp. method)" font ",12"

set bmargin 1.0
set tmargin 6.5
set lmargin 3.5
set rmargin 0.5

set tics font "Helvetica, 10"

set title "H_2 {/Symbol n} = 0" font ",10"
unset key
set ylabel "E_0({/Symbol n},J) [ev]
set xlabel "J"
set xrange[-0.5:11.5]
set yrange [0.25:0.9]
plot "exp.H2.v0"  lc -1 pointtype 7 pointsize 0.5, "SRP48.H2.v0" lc 1 pointtype 7 pointsize 0.5, "B86SRP68-DF2.H2.v0" lc 2 pointtype 7 pointsize 0.5 , "PBEa57-DF2.H2.v0" lc 3 pointtype 7 pointsize 0.5, "MS.B86b.H2.v0" lc 4 pointtype 7 pointsize 0.5

#plot 2
unset ylabel
set title "H_2 {/Symbol n} = 1"
set key t
set key font "Helvetica, 7"
set xrange[-0.5:7.5]
set yrange [0.25:0.9]
plot "exp.H2.v1"  lc -1 pointtype 7 pointsize 0.5 title"exp.", "SRP48.H2.v1" lc 1 pointtype 7 pointsize 0.5 title"SRP48", "B86SRP68-DF2.H2.v1" lc 2 pointtype 7 pointsize 0.5 title"B86SRP68-DF2", "PBEa57-DF2.H2.v1" lc 3 pointtype 7 pointsize 0.5 title"PBEa57-DF2", "MS.B86b.H2.v1" lc 4 pointtype 7 pointsize 0.5 title"MS-B86bl"

#plot 3
set title"D_2 {/Symbol n} = 0"
unset key 
set xlabel"J"
set xrange [-0.5:11.5]
set yrange [0.25:0.9]
#exp
plot "exp.D2.v0"  lc -1 pointtype 7 pointsize 0.5 , "SRP48.D2.v0" lc 1 pointtype 7 pointsize 0.5, "B86SRP68-DF2.D2.v0" lc 2 pointtype 7 pointsize 0.5 , "PBEa57-DF2.D2.v0" lc 3 pointtype 7 pointsize 0.5, "MS.B86b.D2.v0" lc 4 pointtype 7 pointsize 0.5 

#plot 4
unset key
set title"D_2 {/Symbol n} = 1" 
set xlabel"J"
set xrange [-0.5:7.5]
set yrange [0.25:0.9]
plot "exp.D2.v1"  lc -1 pointtype 7 pointsize 0.5, "SRP48.D2.v1" lc 1 pointtype 7 pointsize 0.5, "B86SRP68-DF2.D2.v1" lc 2 pointtype 7 pointsize 0.5 , "PBEa57-DF2.D2.v1" lc 3 pointtype 7 pointsize 0.5, "MS.B86b.D2.v1" lc 4 pointtype 7 pointsize 0.5
EOF

ps2pdf ${FUNC}.ps
#evince ${FUNC}.pdf

FUNC="all_panelA"
gnuplot<<EOF
set output "${FUNC}.ps"
set terminal postscript enhanced color landscape
set multiplot layout 1,4 title " E_0({/Symbol n},J) \n panel A scaling (theory method)" font ",12"

set bmargin 1.0
set tmargin 6.5
set lmargin 3.5
set rmargin 0.5

set tics font "Helvetica, 10"

set title "H_2 {/Symbol n} = 0" font ",10"
unset key
set ylabel "E_0({/Symbol n},J) [ev]
set xlabel "J"
set xrange[-0.5:11.5]
set yrange [0.25:0.9]
plot "exp.H2.v0"  lc -1 pointtype 7 pointsize 0.5, "../E0/SRP48.H2.v0" lc 1 pointtype 7 pointsize 0.5, "../E0/B86SRP68-DF2.H2.v0" lc 2 pointtype 7 pointsize 0.5 , "../E0/PBEa57-DF2.H2.v0" lc 3 pointtype 7 pointsize 0.5, "../E0/MS.B86b.H2.v0" lc 4 pointtype 7 pointsize 0.5

#plot 2
unset ylabel
set title "H_2 {/Symbol n} = 1"
set key t
set key font "Helvetica, 7"
set xrange[-0.5:7.5]
set yrange [0.25:0.9]
plot "exp.H2.v1"  lc -1 pointtype 7 pointsize 0.5 title"exp.", "../E0/SRP48.H2.v1" lc 1 pointtype 7 pointsize 0.5 title"SRP48", "../E0/B86SRP68-DF2.H2.v1" lc 2 pointtype 7 pointsize 0.5 title"B86SRP68-DF2", "../E0/PBEa57-DF2.H2.v1" lc 3 pointtype 7 pointsize 0.5 title"PBEa57-DF2", "../E0/MS.B86b.H2.v1" lc 4 pointtype 7 pointsize 0.5 title"MS-B86bl"

#plot 3
set title"D_2 {/Symbol n} = 0"
unset key 
set xlabel"J"
set xrange [-0.5:11.5]
set yrange [0.25:0.9]
#exp
plot "exp.D2.v0"  lc -1 pointtype 7 pointsize 0.5 , "../E0/SRP48.D2.v0" lc 1 pointtype 7 pointsize 0.5, "../E0/B86SRP68-DF2.D2.v0" lc 2 pointtype 7 pointsize 0.5 , "../E0/PBEa57-DF2.D2.v0" lc 3 pointtype 7 pointsize 0.5, "../E0/MS.B86b.D2.v0" lc 4 pointtype 7 pointsize 0.5 

#plot 4
unset key
set title"D_2 {/Symbol n} = 1" 
set xlabel"J"
set xrange [-0.5:7.5]
set yrange [0.25:0.9]
plot "exp.D2.v1"  lc -1 pointtype 7 pointsize 0.5, "../E0/SRP48.D2.v1" lc 1 pointtype 7 pointsize 0.5, "../E0/B86SRP68-DF2.D2.v1" lc 2 pointtype 7 pointsize 0.5 , "../E0/PBEa57-DF2.D2.v1" lc 3 pointtype 7 pointsize 0.5, "../E0/MS.B86b.D2.v1" lc 4 pointtype 7 pointsize 0.5
EOF

ps2pdf ${FUNC}.ps
#evince ${FUNC}.pdf


pdfunite all_panelB.pdf all_panelA.pdf SRP48.pdf B86SRP68-DF2.pdf PBEa57-DF2.pdf MS.B86b.pdf optPBE-DF1.pdf SRPsol63-DF2.pdf out.pdf