commit f08311efbe3f5173a4021fecfd95a650e1f43b0b
parent 6d14ae3970a33b1a68a60e9a2774fc664f0d90d4
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date: Fri, 31 Jan 2020 09:45:39 +0100
Produce fig2 as single figure
Diffstat:
4 files changed, 184 insertions(+), 33 deletions(-)
diff --git a/fig-stick_slip_depth/Makefile b/fig-stick_slip_depth/Makefile
@@ -0,0 +1,33 @@
+FIG != basename $(PWD)
+BIN = ../1d_fd_simple_shear/max_depth_simple_shear
+
+default: ../$(FIG).pdf
+
+rate/sim.output00001.txt: ../fig-stick_slip_rate/sim.output00001.txt
+ mkdir -p rate
+ cp ../fig-stick_slip_rate/sim.output*.txt rate/
+
+stress/sim.output00001.txt: ../fig-stick_slip_stress/sim.output00001.txt
+ mkdir -p stress
+ cp ../fig-stick_slip_stress/sim.output*.txt stress/
+
+../$(FIG).pdf: fig.gp rate/sim.output00001.txt stress/sim.output00001.txt $(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=$$skin_depth; max_def_depth=$$max_def_depth" \
+ fig.gp > $@'
+
+clean:
+ rm -f rate/*.txt
+ rm -f stress/*.txt
+ rm -f ../$(FIG).pdf
+
+.PHONY: default clean
diff --git a/fig-stick_slip_depth/fig.gp b/fig-stick_slip_depth/fig.gp
@@ -0,0 +1,107 @@
+#!/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 terminal pdfcairo color size 10.68 cm, 4.45 cm font ",10"
+#set terminal pdfcairo color size 11.87 cm, 4.45 cm font ",10"
+set terminal pdfcairo color size 17.8 cm, 4.0 cm font ",10"
+set multiplot layout 1,6 \
+ margins 0.05,0.98,0.23,0.90 \
+ spacing 0.02,0.02
+
+set key bottom right #samplen 0.9
+
+file_interval=60*10
+plot_interval=60*60
+inc=plot_interval/file_interval
+#set yrange [4.0:8.0]
+L_z=8.0
+set yrange [4.0:0.0]
+
+t0 = 47*inc
+tend = 71*inc
+
+filename(t) = sprintf('rate/sim.output%05d.txt', t)
+
+set pm3d map
+#set palette rgb 10,6,6
+set palette
+
+unset colorbox
+set ylabel "{/:Normal Depth [m]}" offset 1.8
+xmin = 0.0
+xmax = 200.0
+set xrange [xmin:xmax]
+set arrow from xmin,skin_depth to xmax,skin_depth nohead lc black front
+set arrow from xmin,max_def_depth to xmax,max_def_depth nohead front lc black dashtype 2
+set title "A" font ",11" offset -7,-0.6
+set xlabel "{/:Normal Water pressure [kPa]}"
+plot for [t = t0:tend:inc] filename(t) u ($4/1e3):(L_z - $1):(t-t0) w l lc palette lw 1 t ""
+
+unset colorbox
+set ylabel ""
+unset ytics
+xmin = 0.0
+xmax = 200.0
+set xrange [xmin:xmax]
+set arrow from xmin,skin_depth to xmax,skin_depth nohead lc black front
+set arrow from xmin,max_def_depth to xmax,max_def_depth nohead front lc black dashtype 2
+set title "B" font ",11" offset -7,-0.6
+set xlabel "{/:Normal Eff. normal stress [kPa]}"
+plot for [t = t0:tend:inc] filename(t) u ($3/1e3):(L_z - $1):(t-t0) w l lc palette lw 1 t ""
+
+set title "C (rate controlled)" font ",11" offset 0,-0.6
+set xlabel "{/:Normal Shear speed [km/a]}"
+set ylabel ""
+unset ytics
+unset arrow
+xmin = -0.1
+xmax = 1.1
+set xtics 0.3
+set xrange [xmin:xmax]
+set arrow from xmin,skin_depth to xmax,skin_depth nohead lc black front
+set arrow from xmin,max_def_depth to xmax,max_def_depth nohead front lc black dashtype 2
+plot for [t = t0:tend:inc] filename(t) u ($2*60*60*24*365/1000):(L_z - $1):(t-t0) w l lc palette lw 1 t ""
+
+set title "D (rate controlled)" font ",11" offset 0,-0.6
+set xlabel "{/:Normal Shear strain rate [1/d]}"
+set ylabel ""
+xmin = -2.0
+xmax = 22.0
+set xrange [xmin:xmax]
+set xtics 5
+set arrow from xmin,skin_depth to xmax,skin_depth nohead lc black front
+set arrow from xmin,max_def_depth to xmax,max_def_depth nohead front lc black dashtype 2
+unset ytics
+plot for [t = t0:tend:inc] filename(t) u ($6*60*60*24):(L_z - $1):(t-t0) w l lc palette lw 1 t ""
+unset format x
+
+filename(t) = sprintf('stress/sim.output%05d.txt', t)
+
+set title "E (stress controlled)" font ",11" offset 0,-0.6
+set xlabel "{/:Normal Shear speed [km/d]}"
+set ylabel ""
+unset ytics
+unset arrow
+xmin = -1
+xmax = 10
+set xrange [xmin:xmax]
+set xtics 2
+set arrow from xmin,skin_depth to xmax,skin_depth nohead lc black front
+set arrow from xmin,max_def_depth to xmax,max_def_depth nohead front lc black dashtype 2
+plot for [t = t0:tend:inc] filename(t) u ($2*60*60*24/1000):(L_z - $1):(t-t0) w l lc palette lw 1 t ""
+
+unset arrow
+set title "F (stress controlled)" font ",11" offset 0,-0.6
+set xlabel "{/:Normal Shear strain rate [1/min]}"
+set ylabel ""
+xmin = -1.0
+xmax = 12.0
+set xrange [xmin:xmax]
+set arrow from xmin,skin_depth to xmax,skin_depth nohead lc black front
+set arrow from xmin,max_def_depth to xmax,max_def_depth nohead front lc black dashtype 2
+unset ytics
+plot for [t = t0:tend:inc] filename(t) u ($6*60):(L_z - $1):(t-t0) w l lc palette lw 1 t ""
+unset format x
+
+unset multiplot
diff --git a/fig-stick_slip_rate_depth/fig.gp b/fig-stick_slip_rate_depth/fig.gp
@@ -2,11 +2,12 @@
#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 terminal pdfcairo color size 10.68 cm, 4.45 cm font ",10"
-set multiplot layout 1,3 \
- margins 0.09,0.98,0.23,0.90 \
- spacing 0.04,0.04
-
+#set terminal pdfcairo color size 10.68 cm, 4.45 cm font ",10"
+#set terminal pdfcairo color size 11.87 cm, 4.45 cm font ",10"
+set terminal pdfcairo color size 11.87 cm, 4.45 cm font ",10"
+set multiplot layout 1,4 \
+ margins 0.05,0.98,0.23,0.90 \
+ spacing 0.02,0.02
set key bottom right #samplen 0.9
@@ -27,18 +28,30 @@ set pm3d map
set palette
unset colorbox
-set ylabel "{/:Normal Depth [m]}" offset 1.0
+set ylabel "{/:Normal Depth [m]}" offset 1.8
xmin = 0.0
xmax = 200.0
set xrange [xmin:xmax]
set arrow from xmin,skin_depth to xmax,skin_depth nohead lc black front
set arrow from xmin,max_def_depth to xmax,max_def_depth nohead front lc black dashtype 2
set label "A" at screen 0.090,0.97 font ",12"
+set xlabel "{/:Normal Water pressure [kPa]}"
+plot for [t = t0:tend:inc] filename(t) u ($4/1e3):(L_z - $1):(t-t0) w l lc palette lw 1 t ""
+
+unset colorbox
+set ylabel ""
+unset ytics
+xmin = 0.0
+xmax = 200.0
+set xrange [xmin:xmax]
+set arrow from xmin,skin_depth to xmax,skin_depth nohead lc black front
+set arrow from xmin,max_def_depth to xmax,max_def_depth nohead front lc black dashtype 2
+set label "B" at screen 0.220,0.97 font ",12"
set xlabel "{/:Normal Eff. normal stress [kPa]}"
plot for [t = t0:tend:inc] filename(t) u ($3/1e3):(L_z - $1):(t-t0) w l lc palette lw 1 t ""
unset label
-set label "B (rate controlled)" at screen 0.395,0.97 font ",12"
+set label "C (rate controlled)" at screen 0.395,0.97 font ",12"
set xlabel "{/:Normal Shear speed [km/a]}"
set ylabel ""
unset ytics
@@ -51,7 +64,7 @@ set arrow from xmin,max_def_depth to xmax,max_def_depth nohead front lc black da
plot for [t = t0:tend:inc] filename(t) u ($2*60*60*24*365/1000):(L_z - $1):(t-t0) w l lc palette lw 1 t ""
unset label
-set label "C (rate controlled)" at screen 0.705,0.97 font ",12"
+set label "D (rate controlled)" at screen 0.705,0.97 font ",12"
set xlabel "{/:Normal Shear strain rate [1/d]}"
set ylabel ""
xmin = -2.0
diff --git a/fig-stick_slip_stress_depth/fig.gp b/fig-stick_slip_stress_depth/fig.gp
@@ -2,12 +2,11 @@
#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 terminal pdfcairo color size 10.68 cm, 4.45 cm font ",10"
-set multiplot layout 1,3 \
- margins 0.09,0.98,0.23,0.90 \
+set terminal pdfcairo color size 5.93 cm, 4.45 cm font ",10"
+set multiplot layout 1,2 \
+ margins 0.03,0.96,0.23,0.90 \
spacing 0.04,0.04
-
set key bottom right #samplen 0.9
file_interval=60*10
@@ -27,42 +26,41 @@ set pm3d map
set palette
unset colorbox
-set ylabel "{/:Normal Depth [m]}" offset 1.0
-xmin = 0.0
-xmax = 200.0
-set xrange [xmin:xmax]
-set arrow from xmin,skin_depth to xmax,skin_depth nohead lc black front
-set arrow from xmin,max_def_depth to xmax,max_def_depth nohead front lc black dashtype 2
-set label "A" at screen 0.090,0.97 font ",12"
-set xlabel "{/:Normal Eff. normal stress [kPa]}"
-plot for [t = t0:tend:inc] filename(t) u ($3/1e3):(L_z - $1):(t-t0) w l lc palette lw 1 t ""
+#set ylabel "{/:Normal Depth [m]}" offset 1.0
+#xmin = 0.0
+#xmax = 200.0
+#set xrange [xmin:xmax]
+#set arrow from xmin,skin_depth to xmax,skin_depth nohead lc black front
+#set arrow from xmin,max_def_depth to xmax,max_def_depth nohead front lc black dashtype 2
+#set label "A" at screen 0.090,0.97 font ",12"
+#set xlabel "{/:Normal Eff. normal stress [kPa]}"
+#plot for [t = t0:tend:inc] filename(t) u ($3/1e3):(L_z - $1):(t-t0) w l lc palette lw 1 t ""
unset label
-set label "B (rate controlled)" at screen 0.395,0.97 font ",12"
-set xlabel "{/:Normal Shear speed [km/a]}"
+set label "E (stress controlled)" at screen 0.020,0.97 font ",12"
+set xlabel "{/:Normal Shear speed [km/d]}"
set ylabel ""
unset ytics
unset arrow
-xmin = -0.1
-xmax = 1.1
+xmin = -1
+xmax = 10
set xrange [xmin:xmax]
set arrow from xmin,skin_depth to xmax,skin_depth nohead lc black front
set arrow from xmin,max_def_depth to xmax,max_def_depth nohead front lc black dashtype 2
-plot for [t = t0:tend:inc] filename(t) u ($2*60*60*24*365/1000):(L_z - $1):(t-t0) w l lc palette lw 1 t ""
+plot for [t = t0:tend:inc] filename(t) u ($2*60*60*24/1000):(L_z - $1):(t-t0) w l lc palette lw 1 t ""
unset label
-set label "C (rate controlled)" at screen 0.705,0.97 font ",12"
-set xlabel "{/:Normal Shear strain rate [1/d]}"
+unset arrow
+set label "F (stress controlled)" at screen 0.510,0.97 font ",12"
+set xlabel "{/:Normal Shear strain rate [1/min]}"
set ylabel ""
-xmin = -2.0
-xmax = 22.0
+xmin = -1.0
+xmax = 12.0
set xrange [xmin:xmax]
set arrow from xmin,skin_depth to xmax,skin_depth nohead lc black front
set arrow from xmin,max_def_depth to xmax,max_def_depth nohead front lc black dashtype 2
-set xrange [xmin:xmax]
unset ytics
-unset arrow
-plot for [t = t0:tend:inc] filename(t) u ($6*60*60*24):(L_z - $1):(t-t0) w l lc palette lw 1 t ""
+plot for [t = t0:tend:inc] filename(t) u ($6*60):(L_z - $1):(t-t0) w l lc palette lw 1 t ""
unset format x
unset multiplot
