tShow skindepth and analytical solution, adjust label placement - cngf-pf-exp1 - experiments for first paper with continuum granular model
(HTM) git clone git://src.adamsgaard.dk/manus_continuum_granular1_exp
(DIR) Log
(DIR) Files
(DIR) Refs
(DIR) Submodules
(DIR) README
(DIR) LICENSE
---
(DIR) commit 8d0714ab70fd06552049eaa83855ccfc6a92ca62
(DIR) parent 8031fcda4847fb2cfc1880358f4d6d5252a5c34d
(HTM) Author: Anders Damsgaard <anders@adamsgaard.dk>
Date: Wed, 18 Dec 2019 15:37:56 +0100
Show skindepth and analytical solution, adjust label placement
Diffstat:
M fig-stick_slip_rate_depth/Makefile | 16 ++++++++++++++--
M fig-stick_slip_rate_depth/fig.gp | 49 ++++++++++++++++---------------
2 files changed, 40 insertions(+), 25 deletions(-)
---
(DIR) diff --git a/fig-stick_slip_rate_depth/Makefile b/fig-stick_slip_rate_depth/Makefile
t@@ -1,12 +1,24 @@
FIG != basename $(PWD)
+BIN = ../1d_fd_simple_shear/max_depth_simple_shear
default: ../$(FIG).pdf
sim.output00001.txt: ../fig-stick_slip_rate/sim.output00001.txt
cp ../fig-stick_slip_rate/sim.output*.txt .
-../$(FIG).pdf: fig.gp sim.output00001.txt ../skindepth.gp
- gnuplot fig.gp > $@
+../$(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
(DIR) diff --git a/fig-stick_slip_rate_depth/fig.gp b/fig-stick_slip_rate_depth/fig.gp
t@@ -1,7 +1,5 @@
#!/usr/bin/env gnuplot
-reset
-
#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,4 \
t@@ -21,47 +19,52 @@ 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)
+#load '../skindepth.gp'
+#k = 2e-17
+#f = 1.0/(3600.0*24.0)
+#f(x) = 8.0 - skindepth(f,k)
+
-set label "{/:Bold a}" at graph 0.05,0.9
+set label "{/:Bold a}" at graph 0.90,0.93
unset colorbox
-set xlabel "Water pressure [kPa]"
-set ylabel "Vertical position [m]"
+set xlabel "{/:Normal Water pressure [kPa]}"
+set ylabel "{/:Normal Vertical position [m]}"
set xrange [0:200]
-plot for [t = t0:tend] filename(t) u ($4/1e3):1:(t-t0) w l lc palette lw 1 t "", \
- f(x) title "" w l
+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] filename(t) u ($4/1e3):1:(t-t0) w l lc palette lw 1 t ""
unset label
-set label "{/:Bold b}" at graph 0.05,0.9
-set xlabel "Eff. normal stress [kPa]"
+set label "{/:Bold b}" at graph 0.90,0.93
+set xlabel "{/:Normal Eff. normal stress [kPa]}"
set ylabel ""
set xrange [0:200]
unset ytics
-plot for [t = t0:tend] filename(t) u ($3/1e3):1:(t-t0) w l lc palette lw 1 t "", \
- f(x) title "" w l
+plot for [t = t0:tend] filename(t) u ($3/1e3):1:(t-t0) w l lc palette lw 1 t ""
unset label
-set label "{/:Bold c}" at graph 0.05,0.9
-set xlabel "Shear velocity [km/a]"
+set label "{/:Bold c}" at graph 0.90,0.93
+set xlabel "{/:Normal Shear velocity [km/a]}"
set ylabel ""
set xrange [0.0:1.0]
#unset xrange
#set xtics (0,0.02,0.04,0.06)
unset ytics
-plot for [t = t0:tend] filename(t) u ($2*60*60*24*365/1000):1:(t-t0) w l lc palette lw 1 t "", \
- f(x) title "" w l
+unset arrow
+set arrow from 0,skin_depth to 1,skin_depth nohead lc black front
+set arrow from 0,max_def_depth to 1,max_def_depth nohead front lc black dashtype 2
+plot for [t = t0:tend] filename(t) u ($2*60*60*24*365/1000):1:(t-t0) w l lc palette lw 1 t ""
unset label
-set label "{/:Bold d}" at graph 0.05,0.9
-set xlabel "Shear strain rate [1/d]"
+set label "{/:Bold d}" at graph 0.90,0.93
+set xlabel "{/:Normal Shear strain rate [1/d]}"
set ylabel ""
set xrange [0:20]
unset ytics
-plot for [t = t0:tend] filename(t) u ($6*60*60*24):1:(t-t0) w l lc palette lw 1 t "", \
- f(x) title "" w l
+unset arrow
+set arrow from 0,skin_depth to 20,skin_depth nohead lc black front
+set arrow from 0,max_def_depth to 20,max_def_depth nohead front lc black dashtype 2
+plot for [t = t0:tend] filename(t) u ($6*60*60*24):1:(t-t0) w l lc palette lw 1 t ""
unset format x
unset multiplot