commit ef1faed8b3b9c7328e9ccb034c6df5695bbe1cd3
parent a7c3437271b1f22e65531d07d23657da57ad015e
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date: Tue, 21 Jan 2020 23:04:59 +0100
Add fig-stick_slip_stress_depth
Diffstat:
2 files changed, 103 insertions(+), 0 deletions(-)
diff --git a/fig-stick_slip_stress_depth/Makefile b/fig-stick_slip_stress_depth/Makefile
@@ -0,0 +1,27 @@
+FIG != basename $(PWD)
+BIN = ../1d_fd_simple_shear/max_depth_simple_shear
+
+default: ../$(FIG).pdf
+
+sim.output00001.txt: ../fig-stick_slip_stress/sim.output00001.txt
+ cp ../fig-stick_slip_stress/sim.output*.txt .
+
+../$(FIG).pdf: fig.gp sim.output00001.txt ../skindepth.gp $(BIN) Makefile
+ sh -c '\
+ out=$$($(BIN) \
+ --normal-stress 200e3 \
+ --fluid-permeability 2e-17 \
+ --fluid-pressure-top 100e3 \
+ --fluid-pressure-ampl 80e3 \
+ --fluid-pressure-freq $$( echo "1.0/(3600*24)" | bc -l )\
+ ); \
+ max_def_depth=$$(echo $$out | cut -d" " -f1); \
+ skin_depth=$$(echo $$out | cut -d" " -f2); \
+ gnuplot -e "skin_depth=8.0-$$skin_depth; max_def_depth=8.0-$$max_def_depth" \
+ fig.gp > $@'
+
+clean:
+ rm -f *.txt
+ rm -f ../$(FIG).pdf
+
+.PHONY: default clean
diff --git a/fig-stick_slip_stress_depth/fig.gp b/fig-stick_slip_stress_depth/fig.gp
@@ -0,0 +1,76 @@
+#!/usr/bin/env gnuplot
+
+#set terminal pdfcairo color size 7.5 cm, 7.5 cm
+set terminal pdfcairo color size 13.35 cm, 4.45 cm font ",10"
+set multiplot layout 1,4 \
+ margins 0.09,0.97,0.23,0.90 \
+ spacing 0.03,0.03
+
+set yrange [4.0:8.0]
+
+set key bottom right #samplen 0.9
+
+file_interval=60*10
+plot_interval=60*60
+inc=plot_interval/file_interval
+
+t0 = 47*inc
+tend = 71*inc
+
+filename(t) = sprintf('sim.output%05d.txt', t)
+
+set pm3d map
+set palette rgb 10,6,6
+
+#load '../skindepth.gp'
+#k = 2e-17
+#f = 1.0/(3600.0*24.0)
+#f(x) = 8.0 - skindepth(f,k)
+
+#set label "{/:Bold b}" at graph 0.90,0.93
+set label "B" at screen 0.09,0.97 font ",12"
+unset colorbox
+set xlabel "{/:Normal Water pressure [kPa]}"
+set ylabel "{/:Normal Vertical position [m]}"
+set xrange [0:200]
+set arrow from 0,skin_depth to 200,skin_depth nohead lc black front
+set arrow from 0,max_def_depth to 200,max_def_depth nohead front lc black dashtype 2
+plot for [t = t0:tend:inc] filename(t) u ($4/1e3):1:(t-t0) w l lc palette lw 1 t ""
+
+unset label
+#set label "{/:Bold c}" at graph 0.90,0.93
+set label "C" at screen 0.315,0.97 font ",12"
+set xlabel "{/:Normal Eff. normal stress [kPa]}"
+set ylabel ""
+set xrange [0:200]
+unset ytics
+plot for [t = t0:tend:inc] filename(t) u ($3/1e3):1:(t-t0) w l lc palette lw 1 t ""
+
+unset label
+#set label "{/:Bold d}" at graph 0.90,0.93
+set label "D" at screen 0.545,0.97 font ",12"
+set xlabel "{/:Normal Shear velocity [km/d]}"
+set ylabel ""
+set xrange [0.0:10.0]
+#unset xrange
+#set xtics (0,0.02,0.04,0.06)
+unset ytics
+unset arrow
+set arrow from 0,skin_depth to 10,skin_depth nohead lc black front
+set arrow from 0,max_def_depth to 10,max_def_depth nohead front lc black dashtype 2
+plot for [t = t0:tend:inc] filename(t) u ($2*60*60*24/1000):1:(t-t0) w l lc palette lw 1 t ""
+
+unset label
+#set label "{/:Bold e}" at graph 0.90,0.93
+set label "E" at screen 0.772,0.97 font ",12"
+set xlabel "{/:Normal Shear strain rate [1/m]}"
+set ylabel ""
+set xrange [0:12]
+unset ytics
+unset arrow
+set arrow from 0,skin_depth to 12,skin_depth nohead lc black front
+set arrow from 0,max_def_depth to 12,max_def_depth nohead front lc black dashtype 2
+plot for [t = t0:tend:inc] filename(t) u ($6*60):1:(t-t0) w l lc palette lw 1 t ""
+unset format x
+
+unset multiplot
