reset
#
set terminal pslatex color rotate auxfile size 2.70,1.70 8
#
# Data sources
MountPoint = '~/SSHFS/'
Folder = 'Helmholtz/Results/'
#
DataFile = 'Results_vs_strain_CIDC_const_p_A_Prefix1_100kPa_kn_12000_mu0.50_DEquil_Template_Frict_0.25_DBumpies-10648-6-06-08-064-c_Frict_0.05_LRforce_-0.20e-4-Int10-1-Nevada-D50.txt'
#
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 2.0 lc rgb "black"
set style line 24 dt 2 lw 0.5 lc rgb "black"
#
set border lw 3.0
#
#set multiplot
#
set xlabel 'Strain, $-\varepsilon_{11}$, percent' offset 0,0.3
set ylabel 'Stress ratio, $q/p$' offset 0.5,0
set y2label '\parbox{1.30in}{\centering Helmholtz energy,\\[0.2ex] $\psi \,/\, p, 10^{-2}$}'
#
set xrange[0:11];  set xtics offset 0,0 2 5; set mxtics 5; #set format x ""
set yrange[0:1.30]; set ytics offset 0.3,0 0.5; set mytics 5; set format y "%.1f"
set ytics nomirror
set y2range[0:1.30]; set y2tics offset -0.3,0 0.5; set my2tics 5; set format y2 "%.1f"
#
set arrow from 1.40,0.3163 to 0.2,0.3163 head ls 23
set label '\footnotesize $\psi_{\text{o}}$' at 1.60,0.3163
#
set arrow from 1.00,0.14 to 0.08,0.010 head ls 23
set label '\footnotesize $q_{\text{o}}$' at 1.20,0.16
#
set key at 4.00,0.500 left top Left reverse samplen 2.5 spacing 1.30 nobox
#
plot MountPoint.Folder.DataFile using (100*$1):($3/$2) \
                                title '\footnotesize $q/p$' \
                                with lines ls 1, \
     MountPoint.Folder.DataFile using (100*$1):(100*$46/$2) \
                                title '\footnotesize $\psi \,/\, p, 10^{-2}$' \
                                with lines ls 2
#
#set size 0.55, 0.55
#set origin 0.30, 0.20
#
#set border lw 2.0
#unset xlabel
#unset ylabel
#set xrange[0:0.3]; set xtics offset 0,0 2  0.2; set mxtics 2; set format x "%.1f"
#set yrange[0:0.8]; set ytics offset 0.3,0  0.5; set mytics 5; set format y "%.1f"
#
#plot MountPoint.Folder.DataFile using (100*$1):($3/$2) \
#                                title '$-\varepsilon_{11}$' \
#                                with lines ls 1
#
#unset multiplot
