cngf-pf-exp1

commit a70425a8dab4349e81422c6fd755470c668042be
parent 815d2596f662ae960930811719b6de9c532eda9e
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Wed,  4 Sep 2019 13:05:01 +0200

Increase amplitude of fluid forcing and do rate controlled shear

Diffstat:
M
fig6/Makefile
 | 
21
+++++++++++++--------
M
fig6/fig.gp
 | 
29
++++++++++++++++++-----------

2 files changed, 31 insertions(+), 19 deletions(-)
diff --git a/fig6/Makefile b/fig6/Makefile
@@ -3,15 +3,20 @@ FIG = $(shell basename $(CURDIR))
 
 default: ../$(FIG).pdf
 
+#	--stress-ratio 0.4 \
+
 sim.output00000.txt: $(BIN) Makefile ../skindepth.gp
-	/bin/sh -c '\
-	./$< --resolution 200 --length 8.0 --normal-stress 150e3 \
-	--stress-ratio 0.4 \
-	--fluid --fluid-permeability 2e-17 --fluid-pressure-top 50e3 \
-	--fluid-pressure-ampl 50e3 \
-	--fluid-pressure-freq $$( echo "1.0/(3600*24)" | bc -l ) \
-	--file-interval $$( echo "3600" | bc -l ) \
-	--time-end $$( echo "3600*24*3 + 3600" | bc -l ) sim'
+	./$< \
+	--length 8.0 \
+	--set-shear-velocity $$(echo "1000/(365*24*60*60)" | bc -l) \
+	--normal-stress 200e3 \
+	--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 "3600" | bc -l) \
+	--time-end $$(echo "3600*24*3 + 3600" | bc -l) sim
 
 ../$(FIG).pdf: fig.gp sim.output00000.txt
 	gnuplot $< > $@
diff --git a/fig6/fig.gp b/fig6/fig.gp
@@ -14,6 +14,7 @@ set key bottom right #samplen 0.9
 
 t0 = 47
 tend = 70
+#tend = 140
 
 filename(t) = sprintf('sim.output%05d.txt', t)
 
@@ -28,7 +29,7 @@ f(x) = 8.0 - skindepth(f,k)
 unset colorbox
 set xlabel "Water pressure [kPa]"
 set ylabel "Vertical position [m]"
-set xrange [0:100]
+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
 
@@ -39,23 +40,29 @@ 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
 
-set xlabel "Shear velocity [m/s]"
+set xlabel "Shear velocity [km/a]"
 set ylabel ""
-set xrange [0.0:0.06]
-set xtics (0,0.02,0.04,0.06)
+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:1:(t-t0) w l lc palette lw 1 t "", \
+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
 
-set xlabel "Shear strain rate [1/s]"
+set xlabel "Shear strain rate [1/d]"
 set ylabel ""
-set xrange [0:3.1e-3]
+#set xrange [0:3.1e-3]
+#unset xrange
+set xrange [0:20]
 #set format x "%.0s*10^{%T}"
-set xtics (0,0.001,0.002,0.003)
+#set xtics (0,0.001,0.002,0.003)
+#set xtics (100,200)
 unset ytics
-prev1 = prev2 = 0
-shift(x) = (prev2 = prev1, prev1 = x)
-plot for [t = t0:tend] filename(t) u (shift($2), $0 < 1 ? 1/0 : $2 - prev2):1:(t-t0) w l lc palette lw 1 t "", \
+#prev1 = prev2 = 0
+#shift(x) = (prev2 = prev1, prev1 = x)
+#plot for [t = t0:tend] filename(t) u (shift($2), $0 < 1 ? 1/0 : $2 - prev2):1:(t-t0) w l lc palette lw 1 t "", \
+#	f(x) title "" w l
+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 format x