commit 88deaea32af656582c48c70399e7035a08c440af
parent 7ec3eef89b0a01f91498a9e648b03e817b4dee67
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date: Sun, 7 Jul 2019 11:17:30 +0200
Add depth of max deformation to hysteresis plot
Diffstat:
3 files changed, 23 insertions(+), 13 deletions(-)
diff --git a/fig2/fig.gp b/fig2/fig.gp
@@ -3,11 +3,6 @@ reset
set terminal pdfcairo color size 7.5 cm, 8.0 cm
-# set xrange [0.0:0.7]
-# set yrange [0.0:0.01]
-# set key top left #samplen 0.9
-
-
set multiplot layout 3,1 \
margins 0.19,0.85,0.23,0.95 \
spacing 0.03,0.03
@@ -20,10 +15,8 @@ plot \
"timeseries.txt" u ($0/1008*7):($2/1000) w l lw 2 lc "black" t "effective normal stress", \
"" u ($0/1008*7):($3/1000) w l lw 2 lt 3 dt "-" t "water pressure"
-#set xtics
unset xtics
unset ylabel
-#unset yrange
unset ytics
set yrange [0:6.0]
set y2range [0:6.0]
@@ -32,11 +25,9 @@ set y2label "Shear velocity [km/d]"
plot "" u ($0/1008*7):($1*60*60*24/1000) w l lw 2 lt 7 t ""
set xtics
-#unset yrange
unset y2label
unset y2tics
set ytics
-#set yrange [0:2]
set yrange [1.4:0.0]
set xlabel "Time [d]"
set ylabel "Depth of max. strain rate [m]"
diff --git a/fig3/Makefile b/fig3/Makefile
@@ -33,7 +33,10 @@ max_strainrate_depth.txt: sim.output00000.txt
END{print depth,\"\t\",max}" "$$f" >> $@; \
done'
-../$(FIG).pdf: fig.gp timeseries.txt max_strainrate_depth.txt
+timeseries_combined.txt: timeseries.txt max_strainrate_depth.txt
+ paste $^ > $@
+
+../$(FIG).pdf: fig.gp timeseries_combined.txt
gnuplot $< > $@
clean:
diff --git a/fig3/fig.gp b/fig3/fig.gp
@@ -1,7 +1,11 @@
#!/usr/bin/env gnuplot
reset
-set terminal pdfcairo color size 7.5 cm, 5.0 cm
+set terminal pdfcairo color size 7.5 cm, 7.5 cm
+
+set multiplot layout 2,1 \
+ margins 0.19,0.75,0.23,0.85 \
+ spacing 0.03,0.03
# set xrange [0.0:0.7]
# set yrange [0.0:0.01]
@@ -10,9 +14,21 @@ set terminal pdfcairo color size 7.5 cm, 5.0 cm
set xrange [40:160]
set x2range [110:-10]
set x2tics 20 nomirror tc lt 3
+unset xtics
-set xlabel "Effective normal stress [kPa]"
set x2label "Water pressure [kPa]" tc lt 3
set ylabel "Shear velocity [km/d]"
plot \
-"timeseries.txt" u ($2/1000):($1*60*60*24/1000) w l lt 7 t ""
+"timeseries_combined.txt" u ($2/1000):($1*60*60*24/1000) w l lt 7 t ""
+
+set xtics
+set xlabel "Effective normal stress [kPa]"
+unset x2tics
+unset ytics
+unset x2label
+unset ylabel
+set y2label "Depth of max. shear strain [m]"
+set yrange [1.4:0.0]
+set y2tics
+set y2range [1.4:0.0]
+plot "" u ($2/1000):(8.0 - $6) w l lt 4 t ""
