reset
#
set terminal pslatex color rotate auxfile size 2.50,1.90 8
#
# Data sources
MountPoint = '~/NUC/'
Folder = 'Helmholtz/Results/'
#
DataFile = 'Moduli_K.txt'
#
unset key
#set grid back
#
# Line styles of plotting DEM results
set style line 1 dt 1 lw 2.5 lc rgb "#8B0000" # dark red       Spheres
set style line 2 dt 12 lw 2.5 lc rgb "#191970" # midnight blue  Ovoids_0800
set style line 3 dt 3 lw 2.0 lc rgb "#006400" # dark green     Ovoids_0625
set style line 4 dt 4 lw 2.0 lc rgb "black"   #                Ovoids_0500
#
set style line 5 dt 5 lw 2.5 lc rgb "purple"
set style line 6 dt 1 lw 3 lc rgb "blue"
set style line 7 dt 5 lw 2.0 lc rgb "orangered4"
#
# black lines
set style line 8 dt 2  lw 1 lc rgb "black"    # Thin black lines
set style line 9 dt 3  lw 1.5 lc rgb "black"    # Zero axis  lines
set style line 10 dt 2 lw 1.4 lc rgb "black"    # Zero axis  lines
set style line 11 dt 6 lw 2.0 lc rgb "#191970"
set style line 12 dt 1 lw 1.4 lc rgb "green"
#
# Experimental (Tsukamoto) lines
set style line 20 dt 3 lw 1.0 lc rgb "black"
set style line 21 dt 3 lw 1.4 lc rgb "green"
set style line 22 dt 1 lw 1.0 lc rgb "black"
set style line 23 dt 1 lw 1.0 lc rgb "black"
set style line 24 dt 2 lw 0.5 lc rgb "black"
set style line 25 dt 2 lw 2.5 lc rgb "#8B0000" # dark red       Spheres
#
set border lw 3.0
#
#set multiplot
#
set xlabel 'Moduli $\mathbf{H}^{\varepsilon\varepsilon}$, measured, $10^{4}$ kPa' offset 0,0.3
set ylabel '\parbox{2.00in}{\centering Moduli $\mathbf{H}^{\varepsilon\varepsilon}$, Eq.~(45),\\[0.4ex] $10^{4}$ kPa}' offset -1.0,0
#
set xrange[0:6]; set xtics offset 0,0 2 2; set mxtics 2; #set format x ""
set yrange[0:6]; set ytics offset 0.3,0 2; set mytics 2; #set format y "%.2f"
#
set arrow from 1,1 to 2,2 nohead ls 1
set arrow from 4,4 to 5.8,5.8 nohead ls 1
set arrow from 0,0 to 0.7,0.7 nohead ls 1
set arrow from 0,0 to 5.8,5.8 nohead ls 25
#
set arrow from 3.15,3.0 to 3.4,3.0 nohead ls 23
set arrow from 3.4,3.0 to 3.4,3.3 nohead ls 23
set label at 3.23,2.8 '\tiny $1$'
set label at 3.45,3.17 '\tiny $1$'
#
set label '\footnotesize{$H^{\varepsilon\varepsilon}_{pv}$, $(0.57)\frac{1}{9}k\ell^{2}\varrho^{\text{c}}$}' \
    at 1.90,1.6
set label '\footnotesize{$H^{\varepsilon\varepsilon}_{q\epsilon}$, $(0.39)\frac{1}{2}k\ell^{2}\varrho^{\text{c}}(1+\zeta)$}' \
    at 0.7,5.2
set label '\footnotesize{$H^{\varepsilon\varepsilon}_{p\epsilon}$, $(1.15)\frac{1}{3}k\ell^{2}\varrho^{\text{c}}\zeta$}' \
    at 0.67,0.45
#
plot MountPoint.Folder.DataFile using (0.0000001*$2):(0.57*0.0000001*$10) \
                                title '$\mathbf{E}^{\text{rev}}_{pv}$' \
                                with points pt 7 ps 0.25 lc rgb "black", \
     MountPoint.Folder.DataFile using (0.0000001*$5):(0.39*0.0000001*$13) \
                                title '$\mathbf{E}^{\text{rev}}_{pv}$' \
                                with points pt 7 ps 0.25 lc rgb "black", \
     MountPoint.Folder.DataFile using (0.0000001*$3):(1.15*0.0000001*$11) \
                                title '$\mathbf{E}^{\text{rev}}_{pv}$' \
                                with points pt 7 ps 0.25 lc rgb "black"
#
#unset multiplot
