commit e88be62300cf88187b45798ba8a1d6e0b3c04724
parent 50440e41d31728921faf47fd9cbc5cb53352b054
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date: Tue, 7 Jan 2020 15:34:59 +0100
Add local experiments for comparison
Diffstat:
7 files changed, 299 insertions(+), 1 deletion(-)
diff --git a/fig-rate_dependence/fig.gp b/fig-rate_dependence/fig.gp
@@ -30,7 +30,7 @@ plot \
"caesar.dat" w l lw 2 t "Caesar till (Rathbun et al., 2008)", \
"storglaciaren.dat" w l lw 2 t "Storglaciaren till (Iverson et al., 1998)", \
"two_rivers.dat" w l lw 2 t "Two Rivers till (Iverson et al., 1998)", \
- "wis.dat" w l lw 2 t "WIS till (Tulaczyk et al., 2001)", \
+ "wis.dat" w l lw 2 t "WIS till (Tulaczyk et al., 2000)", \
"out_b0.94.txt" u ($2*3600*24*365):1 w l lc '#666666' lw 2 dashtype 2 t "{/:Bold current study: NGF model}, b = 0.94", \
"out_b0.10.txt" u ($2*3600*24*365):1 w l lc '#666666' lw 2 dashtype 3 t "{/:Bold current study: NGF model}, b = 0.10"
diff --git a/fig-stick_slip_rate_depth_local/Makefile b/fig-stick_slip_rate_depth_local/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_rate_nonlocality/sim.output00001.txt
+ cp ../fig-stick_slip_rate_nonlocality/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_rate_depth_local/fig.gp b/fig-stick_slip_rate_depth_local/fig.gp
@@ -0,0 +1,74 @@
+#!/usr/bin/env gnuplot
+
+#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 \
+ margins 0.09,0.97,0.23,0.95 \
+ 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
+t0 = 20*inc
+tend = 47*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 a}" at graph 0.90,0.93
+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 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:inc] filename(t) u ($3/1e3):1:(t-t0) w l lc palette lw 1 t ""
+
+unset label
+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
+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:inc] 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.90,0.93
+set xlabel "{/:Normal Shear strain rate [1/d]}"
+set ylabel ""
+set xrange [0:20]
+unset ytics
+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:inc] filename(t) u ($6*60*60*24):1:(t-t0) w l lc palette lw 1 t ""
+unset format x
+
+unset multiplot
diff --git a/fig-stick_slip_rate_local/Makefile b/fig-stick_slip_rate_local/Makefile
@@ -0,0 +1,75 @@
+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 \
+ --set-shear-velocity $$( echo "1000/(365*24*60*60)" | bc -l ) \
+ --grain-size 1e-3 \
+ -A 1e-3 \
+ --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*2" | bc -l ) sim'
+
+timeseries.txt: sim.output00000.txt
+ /bin/sh -c '\
+ rm -f $@; \
+ for f in sim.output*.txt; do \
+ tail -n 1 "$$f" | cut -f2- >> $@; \
+ done'
+
+max_strainrate_depth.txt: sim.output00000.txt
+ /bin/sh -c '\
+ rm -f $@; \
+ for f in sim.output*.txt; do \
+ awk " \
+ BEGIN{getline; depth=8.0; max=0.0} \
+ NF{ \
+ if (\$$6 >= max && \$$6+0 > 1e-300) \
+ { \
+ max=\$$6; \
+ depth=\$$1; \
+ } \
+ } \
+ END{ \
+ print depth,\"\t\",max \
+ }" \
+ "$$f" >> $@; \
+ done'
+
+# central finite differences
+sediment_flux.txt: sim.output00000.txt
+ /bin/sh -c '\
+ rm -f $@; \
+ for f in sim.output*.txt; do \
+ awk " \
+ BEGIN{getline; integral=0.0} \
+ NF{ \
+ if (NR > 1) \
+ { \
+ integral+=(\$$2+v_prev)/2.0*(\$$1-z_prev); \
+ } \
+ z_prev=\$$1; \
+ v_prev=\$$2; \
+ } \
+ END{ \
+ print integral \
+ }" \
+ "$$f" >> $@; \
+ done'
+
+../$(FIG).pdf: fig.gp timeseries.txt max_strainrate_depth.txt sediment_flux.txt
+ gnuplot fig.gp > $@
+
+clean:
+ rm -f *.txt
+ rm -f ../$(FIG).pdf
+
+.PHONY: default clean
diff --git a/fig-stick_slip_rate_local/fig.gp b/fig-stick_slip_rate_local/fig.gp
@@ -0,0 +1,52 @@
+#!/usr/bin/env gnuplot
+reset
+
+set terminal pdfcairo color size 7.5 cm, 10.5 cm
+
+set multiplot layout 4,1 \
+ margins 0.19,0.72,0.10,0.98 \
+ spacing 0.00,0.03
+
+set key top right samplen 1.0 font ",10"
+unset xlabel
+set xtics format ""
+set yrange [0:320]
+set ylabel "{/:Normal [kPa]}"
+set ytics 100
+
+set label "{/:Bold e}" at graph 0.05,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"
+
+unset label
+set label "{/:Bold f}" at graph 0.05,0.9
+unset ylabel
+unset ytics
+set yrange [0:80]
+set y2range [0:80]
+set y2tics 20
+set y2label "{/:Normal Shear stress [kPa]}"
+plot "" u ($0/1008*7):($4*$2/1000) w l lw 2 lt 1 t ""
+
+unset label
+set label "{/:Bold g}" at graph 0.05,0.9
+unset y2label
+unset y2tics
+set ytics 0.5
+set yrange [2.5:0.0]
+set ylabel "{/:Normal Depth of max. strain rate [m]}"
+plot "max_strainrate_depth.txt" u ($0/1008*7):(8.0 - $1) w l lw 2 lt 4 t ""
+
+unset label
+set label "{/:Bold h}" at graph 0.05,0.9
+unset ylabel
+unset ytics
+set yrange [0.0:1e-4]
+set y2range [0.0:1e-4]
+set y2tics 3e-5
+set y2label "{/:Normal Specific sediment flux [m^2/s]}"
+set xtics 1 format "%g"
+set xlabel "{/:Normal Time [d]}"
+plot "sediment_flux.txt" u ($0/1008*7):1 w l lw 2 lt 6 t ""
diff --git a/fig-strain_distribution_local/Makefile b/fig-strain_distribution_local/Makefile
@@ -0,0 +1,29 @@
+BIN = ../1d_fd_simple_shear/1d_fd_simple_shear
+FIG != basename $(PWD)
+
+default: ../$(FIG).pdf
+
+strain_distribution_N10kPa.txt: $(BIN) Makefile
+ /bin/sh -c '\
+ for P in 10 20 40 60 80 120; do \
+ ./$(BIN) \
+ --friction-coefficient 0.4 \
+ --set-shear-velocity 0.0369 \
+ --origo 0.05 --length 0.64 --normal-stress $${P}e3 \
+ --porosity 0.275 \
+ --grain-size 0.04 \
+ -A 1e-3 \
+ -b 0.022 \
+ --normalize > strain_distribution_N$${P}kPa.txt; \
+ done'
+
+../$(FIG).pdf: fig.gp strain_distribution_N10kPa.txt damsgaard2013-sheardisp/sheardisp*-binned.dat
+ gnuplot fig.gp > $@
+
+clean:
+ rm -f *.txt
+ rm -f ../$(FIG).pdf
+
+.PHONY: default clean
+
+
diff --git a/fig-strain_distribution_local/fig.gp b/fig-strain_distribution_local/fig.gp
@@ -0,0 +1,41 @@
+#!/usr/bin/env gnuplot
+
+reset
+
+set terminal pdfcairo enhanced color size 7.5 cm, 5 cm
+
+set margins 8,2,3.2,1
+
+#set lmargin 7.0
+#set bmargin 3.5
+#set rmargin 2.0
+#set tmargin 2.0
+
+set xlabel "{/:Normal Shear displacement {/Symbol D}x [m]}"
+set ylabel "{/:Normal Vertical position z [m]}"
+
+set yrange [0.0:0.73]
+#set xrange [-0.13:1.13]
+set xrange [-0.05:0.9]
+set xtics 0.2
+set ytics 0.2
+
+xdisp=0.74
+
+set label "{/:Bold d}" at screen 0.05,0.95
+
+set key bottom right font ",10" #samplen 0.9
+
+pointsize=0.6
+
+plot "strain_distribution_N10kPa.txt" u ($2*xdisp):1 w l lw 2 lc "#ee2430" title "{/:Bold NGF}, {/Symbol s}_n' 10 kPa", \
+ "damsgaard2013-sheardisp/sheardisp10-binned.dat" w p ps pointsize lw 2 lc "#ee2430" title "DEM, {/Symbol s}_n' 10 kPa", \
+ "strain_distribution_N20kPa.txt" u ($2*xdisp):1 w l lw 2 lc "#53bf3a" title "{/:Bold NGF}, {/Symbol s}_n' 20 kPa", \
+ "damsgaard2013-sheardisp/sheardisp20-binned.dat" w p ps pointsize lw 2 lc "#53bf3a" title "DEM, {/Symbol s}_n' 20 kPa", \
+ "strain_distribution_N40kPa.txt" u ($2*xdisp):1 w l lw 2 lc "#2f54ad" title "{/:Bold NGF}, {/Symbol s}_n' 40 kPa", \
+ "damsgaard2013-sheardisp/sheardisp40-binned.dat" w p ps pointsize lw 2 lc "#2f54ad" title "DEM, {/Symbol s}_n' 40 kPa", \
+ "strain_distribution_N80kPa.txt" u ($2*xdisp):1 w l lw 2 lc "#6ccfd9" title "{/:Bold NGF}, {/Symbol s}_n' 80 kPa", \
+ "damsgaard2013-sheardisp/sheardisp80-binned.dat" w p ps pointsize lw 2 lc "#6ccfd9" title "DEM, {/Symbol s}_n' 80 kPa", \
+ #"strain_distribution_N120kPa.txt" u ($2*xdisp):1 w l lw 2 lc "#231f20" title "{/:Bold NGF}, {/Symbol s}_n' 120 kPa", \
+ #"damsgaard2013-sheardisp/sheardisp120-binned.dat" w p ps pointsize lw 2 lc "#231f20" title "DEM, {/Symbol s}_n' 120 kPa"
+
