commit 1786c984cde1baa831b14aa781c59aa681d5dbf2
parent d3c58ac566fb197debb9ce36024fee87634d4135
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date: Mon, 18 Nov 2019 11:15:21 +0100
Put Mohr-Coulomb analysis into separate figure, still needs tweaks
Diffstat:
14 files changed, 115 insertions(+), 0 deletions(-)
diff --git a/fig-mohr_coulomb/Makefile b/fig-mohr_coulomb/Makefile
@@ -0,0 +1,60 @@
+BIN = ../1d_fd_simple_shear/1d_fd_simple_shear
+FIG != basename $(PWD)
+N_SEQ != seq 1 20 220
+
+default: ../$(FIG).pdf
+
+# mimic Iverson 2010 fig 2b
+# Ceasar till (Rathbun et al., 2008)
+mc_ceasar.txt: $(BIN) Makefile
+ rm -f $@
+ for N in $(N_SEQ); do \
+ printf "$$N\t" >> $@; \
+ ./$(BIN) -m 0.5634 -f 0.5634 -C 29.6e3 -P $$N -L 1.0 mc | \
+ tail -n 1 | cut -f4 >> $@; \
+ done
+
+# Storglaciaren till (Iverson et al., 1998)
+mc_storglaciaren.txt: $(BIN) Makefile
+ rm -f $@
+ for N in $(N_SEQ); do \
+ printf "$$N\t" >> $@; \
+ ./$(BIN) -m 0.494 -f 0.494 -C 5.0e3 -P $$N -L 1.0 mc | \
+ tail -n 1 | cut -f4 >> $@; \
+ done
+
+# Two Rivers till (Iversonet al., 1998)
+mc_two_rivers.txt: $(BIN) Makefile
+ rm -f $@
+ for N in $(N_SEQ); do \
+ printf "$$N\t" >> $@; \
+ ./$(BIN) -m 0.321 -f 0.321 -C 14e3 -P $$N -L 1.0 mc | \
+ tail -n 1 | cut -f4 >> $@; \
+ done
+
+# Des Moines Lobe till (Iverson et al., 2010)
+mc_des_moines_lobe.txt: $(BIN) Makefile
+ rm -f $@
+ for N in $(N_SEQ); do \
+ printf "$$N\t" >> $@; \
+ ./$(BIN) -m 0.327 -f 0.327 -C 4.23e3 -P $$N -L 1.0 mc | \
+ tail -n 1 | cut -f4 >> $@; \
+ done
+
+# WIS till (Tulaczyk et al., 2001)
+mc_wis.txt: $(BIN) Makefile
+ rm -f $@
+ for N in $(N_SEQ); do \
+ printf "$$N\t" >> $@; \
+ ./$(BIN) -m 0.443 -f 0.443 -C 3.0e3 -P $$N -L 1.0 mc | \
+ tail -n 1 | cut -f4 >> $@; \
+ done
+
+../$(FIG).pdf: fig.gp mc_ceasar.txt mc_storglaciaren.txt mc_two_rivers.txt mc_des_moines_lobe.txt mc_wis.txt *.dat
+ gnuplot fig.gp > $@
+
+clean:
+ rm -f *.txt
+ rm -f ../$(FIG).pdf
+
+.PHONY: default clean
diff --git a/fig-mohr_coulomb/ceasar.dat b/fig-mohr_coulomb/ceasar.dat
@@ -0,0 +1,2 @@
+3.0e3 0.54
+9.0e5 0.58
diff --git a/fig-mohr_coulomb/cowden.dat b/fig-mohr_coulomb/cowden.dat
@@ -0,0 +1,2 @@
+3.2e3 0.52
+5.0e5 0.52
diff --git a/fig-mohr_coulomb/des_moines_lobe.dat b/fig-mohr_coulomb/des_moines_lobe.dat
@@ -0,0 +1,2 @@
+9.0e2 0.505
+4.0e4 0.45
diff --git a/fig-mohr_coulomb/fig.gp b/fig-mohr_coulomb/fig.gp
@@ -0,0 +1,31 @@
+#!/usr/bin/env gnuplot
+
+reset
+
+## resemble Iverson 2010 fig2a
+
+set terminal pdfcairo size 7.5 cm, 7.5 cm
+#set terminal pdfcairo size 15 cm, 10 cm enhanced
+
+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_ceasar.txt" u ($1):($2*$1 - 6.5) w l lt 3 lw 2 t "Simulated Ceasar till", \
+ "mc_storglaciaren.dat" w l lt 4 lw 2 t "Storglaciaren till (Iverson et al., 1998)", \
+ "mc_storglaciaren.txt" u ($1):($2*$1 + 5) w l lt 4 lw 2 t "Simulated Storglaciaren till", \
+ "mc_two_rivers.dat" w l lt 5 lw 2 t "Two Rivers till (Iverson et al., 1998)", \
+ "mc_two_rivers.txt" u ($1):($2*$1 + 14) w l lt 5 lw 2 t "Simulated Two Rivers till", \
+ "mc_des_moines_lobe.dat" w l lt 6 lw 2 t "Des Moines Lobe till (Iverson, 2010)", \
+ "mc_des_moines_lobe.txt" u ($1):($2*$1 + 4.23) w l lt 6 lw 2 t "Simulated Des Moines Lobe till", \
+ "mc_wis.dat" w l lt 7 lw 2 t "WIS till (Tulaczyk et al., 2001)", \
+ "mc_wis.txt" u ($1):($2*$1 + 3) w l lt 7 lw 2 t "Simulated WIS till"
diff --git a/fig-mohr_coulomb/lower_cromer.dat b/fig-mohr_coulomb/lower_cromer.dat
@@ -0,0 +1,2 @@
+5.0e1 0.60
+7.0e4 0.56
diff --git a/fig-mohr_coulomb/mc_ceasar.dat b/fig-mohr_coulomb/mc_ceasar.dat
@@ -0,0 +1,2 @@
+48 19
+200 100
diff --git a/fig-mohr_coulomb/mc_des_moines_lobe.dat b/fig-mohr_coulomb/mc_des_moines_lobe.dat
@@ -0,0 +1,2 @@
+10 7.5
+140 50
diff --git a/fig-mohr_coulomb/mc_storglaciaren.dat b/fig-mohr_coulomb/mc_storglaciaren.dat
@@ -0,0 +1,2 @@
+20 17
+135 75
diff --git a/fig-mohr_coulomb/mc_two_rivers.dat b/fig-mohr_coulomb/mc_two_rivers.dat
@@ -0,0 +1,2 @@
+30 24
+140 60
diff --git a/fig-mohr_coulomb/mc_wis.dat b/fig-mohr_coulomb/mc_wis.dat
@@ -0,0 +1,2 @@
+25 13
+200 90
diff --git a/fig-mohr_coulomb/storglaciaren.dat b/fig-mohr_coulomb/storglaciaren.dat
@@ -0,0 +1,2 @@
+5.3e2 0.57
+2.3e4 0.52
diff --git a/fig-mohr_coulomb/two_rivers.dat b/fig-mohr_coulomb/two_rivers.dat
@@ -0,0 +1,2 @@
+1.8e3 0.475
+8.0e4 0.44
diff --git a/fig-mohr_coulomb/wis.dat b/fig-mohr_coulomb/wis.dat
@@ -0,0 +1,2 @@
+1.0e1 0.504
+1.0e5 0.495
