commit b6947c9d3614320a252e75677be9ce0869a2488b
parent e2a4b60861f06f1fd2fd4ab5648835a02f0a41f0
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date: Wed, 8 Jan 2020 12:20:56 +0100
Add figure only showing stick-slip experiment forcing
Diffstat:
2 files changed, 48 insertions(+), 0 deletions(-)
diff --git a/fig-stick_slip_forcing/Makefile b/fig-stick_slip_forcing/Makefile
@@ -0,0 +1,26 @@
+BIN = ../1d_fd_simple_shear/1d_fd_simple_shear
+FIG != basename $(PWD)
+
+default: ../$(FIG).pdf
+
+sim.output00000.txt: $(BIN)
+ /bin/sh -c '\
+ ./$(BIN) --length 8.0 \
+ --normal-stress 200e3 \
+ --stress-ratio 0.40 \
+ --fluid \
+ --fluid-permeability 2e-17 \
+ --fluid-pressure-top 100e3 \
+ --fluid-pressure-ampl 80e3 \
+ --fluid-pressure-freq $$( echo "1.0/(3600*24)" | bc -l ) \
+ --file-interval $$( echo "60*10" | bc -l ) \
+ --time-end $$( echo "3600*24*7" | bc -l ) sim'
+
+../$(FIG).pdf: fig.gp sim.output00000.txt
+ gnuplot fig.gp > $@
+
+clean:
+ rm -f *.txt
+ rm -f ../$(FIG).pdf
+
+.PHONY: default clean
diff --git a/fig-stick_slip_forcing/fig.gp b/fig-stick_slip_forcing/fig.gp
@@ -0,0 +1,22 @@
+#!/usr/bin/env gnuplot
+reset
+
+set terminal pdfcairo color size 7.5 cm, 3.5 cm
+
+set key top right samplen 1.0 font ",10"
+set xtics 1 format "%g"
+set xlabel "{/:Normal Time [d]}"
+set yrange [0:320]
+set ylabel "{/:Normal [kPa]}"
+set ytics 100
+
+#set label "{/:Bold a}" at screen 0.06,0.99
+#set label "{/:Bold b}" at screen 0.06,0.75
+#set label "{/:Bold c}" at screen 0.06,0.54
+#set label "{/:Bold d}" at screen 0.06,0.30
+
+set label "{/:Bold a}" at graph 0.02,0.9
+
+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"