commit d3c58ac566fb197debb9ce36024fee87634d4135
parent 2f19733c517f53aab139d18d20b062176cefd854
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date: Mon, 18 Nov 2019 10:49:30 +0100
Remove Mohr-Coulomb analysis from rate-dependence figure
Diffstat:
8 files changed, 10 insertions(+), 70 deletions(-)
diff --git a/damsgaard2013-fig8.pdf b/damsgaard2013-fig8.pdf
Binary files differ.
diff --git a/fig-rate_dependence/Makefile b/fig-rate_dependence/Makefile
@@ -1,5 +1,5 @@
BIN = ../1d_fd_simple_shear/1d_fd_simple_shear
-FIG = $(shell basename $(CURDIR))
+FIG != basename $(PWD)
default: ../$(FIG).pdf
@@ -8,7 +8,7 @@ default: ../$(FIG).pdf
# friction around 0.5
sim.output00000.txt: $(BIN) Makefile
/bin/sh -c '\
- for b in $$(printf "0.01\n0.10\n"; seq 0.20 0.20 0.90) 0.94; do \
+ for b in 0.01 0.10 0.20 0.40 0.94; do \
out="out_b$${b}.txt"; \
rm -f "$$out"; \
for t in $$(seq 0.0001 0.002 1.0); do \
@@ -17,24 +17,11 @@ sim.output00000.txt: $(BIN) Makefile
tail -n 1 | cut -f2 >> "$$out"; \
done; done'
-# mimic Iverson 2010 fig 2b
-mc.output00000.txt: $(BIN) Makefile
- /bin/sh -c '\
- for mu in $$(seq 0.30 0.05 0.50); do \
- for C in $$(printf "0\n20e3"); do \
- out="out_mu$${mu}_C$${C}.txt"; \
- rm -f "$$out"; \
- for N in $$(seq 1 20 220); do \
- printf "$$N\t" >> "$$out"; \
- ./$(BIN) -m $$mu -f $$mu -C $$C -P $$N -L 1.0 mc | \
- tail -n 1 | cut -f4 >> "$$out"; \
- done; done; done'
-
-../$(FIG).pdf: fig.gp sim.output00000.txt mc.output00000.txt *.dat
- gnuplot $< > $@
+../$(FIG).pdf: fig.gp sim.output00000.txt *.dat
+ gnuplot fig.gp > $@
clean:
- $(RM) *.txt
- $(RM) ../$(FIG).pdf
+ rm -f *.txt
+ rm -f ../$(FIG).pdf
.PHONY: default clean
diff --git a/fig-rate_dependence/fig.gp b/fig-rate_dependence/fig.gp
@@ -5,9 +5,9 @@ reset
## resemble Iverson 2010 fig2a
#set terminal pdfcairo size 7.5 cm, 7.5 cm
-set terminal pdfcairo size 15 cm, 10 cm enhanced
+set terminal pdfcairo size 7.5 cm, 10 cm enhanced
-set multiplot layout 2,2 \
+set multiplot layout 2,1 \
margins 0.09,0.97,0.11,0.97 \
spacing 0.03,0.15
@@ -21,10 +21,8 @@ set logscale x
set ylabel "Shear friction {/Symbol m} [-]"
set xlabel "Shear strain rate d{/Symbol g}/dt [a^{-1}]"
-set label "{/:Bold a}" at screen 0.01,0.96
-set label "{/:Bold b}" at screen 0.525,0.96
-set label "{/:Bold c}" at screen 0.01,0.45
-set label "{/:Bold d}" at screen 0.525,0.45
+set label "{/:Bold a}" at screen 0.11,0.95
+set label "{/:Bold b}" at screen 0.11,0.44
plot "cowden.dat" w l lw 2 t "Cowden till (Tika et al., 1996)", \
"lower_cromer.dat" w l lw 2 t "Lower Cromer till (Tika et al., 1996)", \
@@ -46,38 +44,3 @@ plot "out_b0.94.txt" u ($2*3600*24*365):1 w l lc '#666666' lw 2 dashtype '-' t "
"out_b0.20.txt" u ($2*3600*24*365):1 w l ls 5 lw 2 t "b = 0.20", \
"out_b0.10.txt" u ($2*3600*24*365):1 w l ls 3 lw 2 t "b = 0.10", \
"out_b0.01.txt" u ($2*3600*24*365):1 w l ls 1 lw 2 t "b = 0.01"
-
-unset xrange
-unset yrange
-unset logscale x
-set xtics 50
-set ytics 20
-set ylabel "Shear stress {/Symbol t} [kPa]"
-set xlabel "Effective normal stress {/Symbol s}_n' [kPa]"
-set xrange [0:200]
-set yrange [0:120]
-
-set key top left reverse Left font ",8" samplen 2.0
-
-plot "mc_ceasar.dat" w l lt 3 lw 2 t "Ceasar till (Rathbun et al., 2008)", \
- "mc_storglaciaren.dat" w l lt 4 lw 2 t "Storglaciaren till (Iverson et al., 1998)", \
- "mc_two_rivers.dat" w l lt 5 lw 2 t "Two Rivers till (Iverson et al., 1998)", \
- "mc_des_moines_lobe.dat" w l lt 6 lw 2 t "Des Moines Lobe till (Iverson, 2010)", \
- "mc_wis.dat" w l lt 7 lw 2 t "WIS till (Tulaczyk et al., 2001)"
-
-unset ylabel
-unset ytics
-load 'BlueSeq.plt'
-set key top left reverse Left font ",8" samplen 2.0
-
-plot \
- "out_mu0.50_C20e3.txt" u ($1):($2*$1 + 20) w l lc '#FEE0D2' lw 2 t "{/Symbol m}_s = 0.3, C = 20 kPa", \
- "out_mu0.45_C20e3.txt" u ($1):($2*$1 + 20) w l lc '#FCBBA1' lw 2 t "", \
- "out_mu0.40_C20e3.txt" u ($1):($2*$1 + 20) w l lc '#FC9272' lw 2 t "{/Symbol m}_s = 0.4, C = 20 kPa", \
- "out_mu0.35_C20e3.txt" u ($1):($2*$1 + 20) w l lc '#FB6A4A' lw 2 t "", \
- "out_mu0.30_C20e3.txt" u ($1):($2*$1 + 20) w l lc '#EF3B2C' lw 2 t "{/Symbol m}_s = 0.5, C = 20 kPa", \
- "out_mu0.50_C0.txt" u ($1):($2*$1) w l ls 5 lw 2 t "{/Symbol m}_s = 0.5, C = 0 kPa", \
- "out_mu0.45_C0.txt" u ($1):($2*$1) w l ls 4 lw 2 t "", \
- "out_mu0.40_C0.txt" u ($1):($2*$1) w l ls 3 lw 2 t "{/Symbol m}_s = 0.4, C = 0 kPa", \
- "out_mu0.35_C0.txt" u ($1):($2*$1) w l ls 2 lw 2 t "", \
- "out_mu0.30_C0.txt" u ($1):($2*$1) w l ls 1 lw 2 t "{/Symbol m}_s = 0.3, C = 0 kPa"
diff --git a/fig-rate_dependence/mc_ceasar.dat b/fig-rate_dependence/mc_ceasar.dat
@@ -1,2 +0,0 @@
-48 19
-200 100
diff --git a/fig-rate_dependence/mc_des_moines_lobe.dat b/fig-rate_dependence/mc_des_moines_lobe.dat
@@ -1,2 +0,0 @@
-10 7.5
-140 50
diff --git a/fig-rate_dependence/mc_storglaciaren.dat b/fig-rate_dependence/mc_storglaciaren.dat
@@ -1,2 +0,0 @@
-20 17
-135 75
diff --git a/fig-rate_dependence/mc_two_rivers.dat b/fig-rate_dependence/mc_two_rivers.dat
@@ -1,2 +0,0 @@
-30 24
-140 60
diff --git a/fig-rate_dependence/mc_wis.dat b/fig-rate_dependence/mc_wis.dat
@@ -1,2 +0,0 @@
-25 13
-200 90
