reset
#
set terminal pslatex color rotate auxfile size 3.00,1.80 8
#
# Data sources
MountPoint = '~/NUC/'
Folder = 'Helmholtz/Results/'
#
DataFile1 = 'Moduli_K.txt'
DataFile2 = 'Moduli_K_spline.txt'
#
D=0.000165
kn=6000
#
unset key
#set grid back
#
# Line styles of plotting DEM results
set style line 1 dt 1 lw 2.0 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 0.5 lc rgb "black"
set style line 24 dt 2 lw 0.5 lc rgb "black"
#
#set multiplot
#
set border lw 3.0
#
set xlabel 'Strain, $-\varepsilon_{11}$, percent' offset 0,0.3
set ylabel 'Moduli, $H_{\!\ast\!\ast}^{\,\varepsilon\varepsilon}\!\cdot (\ell /k^{\text{n}})$' offset 0.3,0
set y2label 'Moduli, $H_{\!\ast\!\ast}^{\,\varepsilon\varepsilon}\: /\;p$' offset 0.0,0
#
set xrange[0:11];  set xtics offset 0,0 5; set mxtics 5; #set format x ""
set yrange[0:1.80]; set ytics nomirror offset 0.3,0 0.5; set mytics 5; set format y "%.1f"
set y2range[0:364]; set y2tics nomirror offset -0.3,0 100; set my2tics 2; #set format y "%.1f"
#
#set xzeroaxis ls 23
#
#set key at 0.70,1.20 left top Left reverse samplen 3.5 spacing 1.50 nobox
unset key
#
set label '\footnotesize $H_{q\epsilon}^{\varepsilon\varepsilon}$' at 2.60,1.30
set label '\footnotesize $H_{pv}^{\varepsilon\varepsilon}$' at 1.60,0.50
set label '\footnotesize $H_{p\epsilon}^{\varepsilon\varepsilon}$' at 0.55,0.26
#
plot MountPoint.Folder.DataFile1 using (100*$1):($2*D/kn) \
                                 notitle \
                                 with points pt 6 ps 0.4 lc rgb "black", \
     MountPoint.Folder.DataFile2 using (100*$1):($2*D/kn) \
                title '$K_{pv}^{\varepsilon}\cdot (\overline{D}/k^{\text{n}})$' \
                                 with lines ls 1, \
     MountPoint.Folder.DataFile1 using (100*$1):($5*D/kn) \
                                 notitle \
                                 with points pt 6 ps 0.4 lc rgb "black", \
     MountPoint.Folder.DataFile2 using (100*$1):($5*D/kn) \
       title '$K_{q\epsilon}^{\varepsilon}\cdot (\overline{D}/k^{\text{n}})$' \
                                 with lines ls 1, \
     MountPoint.Folder.DataFile1 using (100*$1):($3*D/kn) \
                                 notitle \
                                 with points pt 6 ps 0.4 lc rgb "black", \
     MountPoint.Folder.DataFile2 using (100*$1):($3*D/kn) \
       title '$K_{p\epsilon}^{\varepsilon}\cdot (\overline{D}/k^{\text{n}})$' \
                                 with lines ls 1
#
#set size 0.45,0.4
#set origin 0.45,0.3
#unset label
#unset xlabel
#unset ylabel
#set xrange[0:0.6];  set xtics offset 0,0 0.5; set mxtics 5; #set format x ""
#set yrange[0:1.7];  set ytics offset 0.3,0 1; set mytics 10; set format y "%.1f"
#plot MountPoint.Folder.DataFile1 using (100*$1):($2*D/kn) \
#                                 notitle \
#                                 with points pt 6 ps 0.4 lc rgb "black", \
#     MountPoint.Folder.DataFile2 using (100*$1):($2*D/kn) \
#                title '$K_{pv}^{\text{rev}}\cdot (\overline{D}/k^{\text{n}})$' \
#                                 with lines ls 1, \
#     MountPoint.Folder.DataFile1 using (100*$1):($5*D/kn) \
#                                 notitle \
#                                 with points pt 6 ps 0.4 lc rgb "black", \
#     MountPoint.Folder.DataFile2 using (100*$1):($5*D/kn) \
#       title '$K_{q\epsilon}^{\text{rev}}\cdot (\overline{D}/k^{\text{n}})$' \
#                                 with lines ls 1, \
#     MountPoint.Folder.DataFile1 using (100*$1):($3*D/kn) \
#                                 notitle \
#                                 with points pt 6 ps 0.4 lc rgb "black", \
#     MountPoint.Folder.DataFile2 using (100*$1):($3*D/kn) \
#       title '$K_{p\epsilon}^{\text{rev}}\cdot (\overline{D}/k^{\text{n}})$' \
#                                 with lines ls 1
#unset multiplot
