commit 1a08d851fee178eaa3fe67b61214a440ac16ddd7
parent 939a81fa07805ae76ed51fddf86b5097d765e3af
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date: Mon, 15 Jul 2019 16:13:05 +0200
Add MC results from Iverson2010
Diffstat:
6 files changed, 36 insertions(+), 6 deletions(-)
diff --git a/fig1/fig.gp b/fig1/fig.gp
@@ -55,17 +55,25 @@ set xlabel "Effective normal stress [kPa]"
set xrange [0:200]
set yrange [0:120]
+plot "mc_ceasar.dat" w l lw 2 t "Ceasar till (Rathbun et al., 2008)", \
+ "mc_storglaciaren.dat" w l lw 2 t "Storglaciaren till (Iverson et al., 1998)", \
+ "mc_two_rivers.dat" w l lw 2 t "Two Rivers till (Iverson et al., 1998)", \
+ "mc_des_moines_lobe.dat" w l lw 2 t "Des Moines Lobe till (Iverson, 2010)", \
+ "mc_wis.dat" w l lw 2 t "WIS till (Tulaczyk et al., 2001)"
+
+unset ylabel
+unset ytics
load 'BlueSeq.plt'
set key top left reverse Left font ",7" samplen 2.0
plot \
- "out_mu0.50_C20e3.txt" u ($1):($2*$1 + 20) w l lc '#FEE0D2' lw 2 t "mu_s = 0.30, C = 20 kPa", \
+ "out_mu0.50_C20e3.txt" u ($1):($2*$1 + 20) w l lc '#FEE0D2' lw 2 t "mu_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 "mu_s = 0.40, C = 20 kPa", \
+ "out_mu0.40_C20e3.txt" u ($1):($2*$1 + 20) w l lc '#FC9272' lw 2 t "mu_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 "mu_s = 0.50, C = 20 kPa", \
- "out_mu0.50_C0.txt" u ($1):($2*$1) w l ls 5 lw 2 t "mu_s = 0.50, C = 0 kPa", \
+ "out_mu0.30_C20e3.txt" u ($1):($2*$1 + 20) w l lc '#EF3B2C' lw 2 t "mu_s = 0.5, C = 20 kPa", \
+ "out_mu0.50_C0.txt" u ($1):($2*$1) w l ls 5 lw 2 t "mu_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 "mu_s = 0.40, C = 0 kPa", \
+ "out_mu0.40_C0.txt" u ($1):($2*$1) w l ls 3 lw 2 t "mu_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 "mu_s = 0.30, C = 0 kPa"
+ "out_mu0.30_C0.txt" u ($1):($2*$1) w l ls 1 lw 2 t "mu_s = 0.3, C = 0 kPa"
diff --git a/fig1/mc_ceasar.dat b/fig1/mc_ceasar.dat
@@ -0,0 +1,5 @@
+48 19
+75 33
+100 45
+150 73
+200 100
diff --git a/fig1/mc_des_moines_lobe.dat b/fig1/mc_des_moines_lobe.dat
@@ -0,0 +1,4 @@
+10 7.5
+50 20
+100 38
+140 50
diff --git a/fig1/mc_storglaciaren.dat b/fig1/mc_storglaciaren.dat
@@ -0,0 +1,4 @@
+20 17
+50 30
+100 53
+135 75
diff --git a/fig1/mc_two_rivers.dat b/fig1/mc_two_rivers.dat
@@ -0,0 +1,4 @@
+30 24
+50 30
+100 46
+140 60
diff --git a/fig1/mc_wis.dat b/fig1/mc_wis.dat
@@ -0,0 +1,5 @@
+25 13
+50 24
+100 46
+150 68
+200 90
