manus_continuum_granular1_exp

commit 3b5829e907d5977e03091ffc736399aa02e38b0d
parent 27fa480003a51e9edfeb139a15484075b5d7faf3
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Tue, 25 Jun 2019 15:21:09 +0200

Add shear-strain rate subplot and fix x-axis tics

Diffstat:
M
fig3/fig.gp
 | 
20
++++++++++++++++----

1 file changed, 16 insertions(+), 4 deletions(-)
diff --git a/fig3/fig.gp b/fig3/fig.gp
@@ -2,9 +2,9 @@
 
 #set terminal pdfcairo color size 7.5 cm, 7.5 cm
 set terminal pdfcairo color size 15.0 cm, 5.0 cm
-set multiplot layout 1,3 \
-	margins 0.1,0.98,0.23,0.95 \
-	spacing 0.04,0.04
+set multiplot layout 1,4 \
+	margins 0.09,0.97,0.23,0.95 \
+	spacing 0.03,0.03
 
 set yrange [0.0:2.0]
 
@@ -20,7 +20,7 @@ set ylabel "Vertical position, z [m]"
 set xrange [0:100]
 plot for [t = t0:tend] filename(t) u ($4/1e3):1 w l lw 1 lc "blue" t ""
 
-set xlabel "Effective normal stress [kPa]"
+set xlabel "Eff. normal stress [kPa]"
 set ylabel ""
 set xrange [50:150]
 unset ytics
@@ -29,7 +29,19 @@ plot for [t = t0:tend] filename(t) u ($3/1e3):1 w l lw 1 lc "black" t ""
 set xlabel "Shear velocity [m/s]"
 set ylabel ""
 set xrange [0.0:0.023]
+set xtics (0,0.01,0.02)
 unset ytics
 plot for [t = t0:tend] filename(t) u 2:1 w l lw 1 lc "red" t ""
 
+set xlabel "Shear strain rate [1/s]"
+set ylabel ""
+set xrange [0:3.1e-3]
+#set format x "%.0s*10^{%T}"
+set xtics (0,0.001,0.002,0.003)
+unset ytics
+prev1 = prev2 = 0
+shift(x) = (prev2 = prev1, prev1 = x)
+plot for [t = t0:tend] filename(t) u (shift($2), $0 < 1 ? 1/0 : $2 - prev2):1 w l lw 1 lc "magenta" t ""
+unset format x
+
 unset multiplot