commit 9d9525fb29cc0967cb2e574eee7c0169bac72584
parent 17c771e1335b79f45a98ef834fdb89556f79e5b7
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date: Thu, 4 Jul 2019 14:30:24 +0200
Improve fig1 makefile and add data points from Iverson 2010
Diffstat:
10 files changed, 40 insertions(+), 20 deletions(-)
diff --git a/fig1/Makefile b/fig1/Makefile
@@ -6,7 +6,7 @@ default: ../$(FIG).pdf
# mimic Iverson 2010 fig 2a
# shear-strain rate from 10^1 to 10^6 m/a
# friction around 0.55
-../$(FIG).pdf: fig.gp
+unnamed.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 \
out="out_b$${b}.txt"; \
@@ -16,6 +16,8 @@ default: ../$(FIG).pdf
./$(BIN) -f 0.55 -P 100e3 -L 1.0 --stress-ratio $$t -b $$b | \
tail -n 1 | cut -f2 >> "$$out"; \
done; done'
+
+../$(FIG).pdf: fig.gp unnamed.output00000.txt
gnuplot $< > $@
clean:
diff --git a/fig1/ceasar.dat b/fig1/ceasar.dat
@@ -0,0 +1,2 @@
+3.0e3 0.54
+9.0e5 0.58
diff --git a/fig1/cowden.dat b/fig1/cowden.dat
@@ -0,0 +1,2 @@
+3.2e3 0.52
+5.0e5 0.52
diff --git a/fig1/des_moines_lobe.dat b/fig1/des_moines_lobe.dat
@@ -0,0 +1,2 @@
+9.0e2 0.505
+4.0e4 0.45
diff --git a/fig1/fig.gp b/fig1/fig.gp
@@ -4,17 +4,35 @@ reset
## resemble Iverson 2010 fig2a
-set terminal pdfcairo size 7.5 cm, 5.0 cm
+#set terminal pdfcairo size 7.5 cm, 5.0 cm
+set terminal pdfcairo size 7.5 cm, 8.0 cm
-set xlabel "Shear strain rate [a^{-1}]"
-set ylabel "Shear friction [-]"
+
+set multiplot layout 2,1 \
+ margins 0.19,0.95,0.23,0.95 \
+ spacing 0.03,0.03
set xrange [1e1:1e6] # m/a
set yrange [0.0:0.8]
-set key bottom right #samplen 0.9
-
+set key bottom left #samplen 0.9
set logscale x
+unset xlabel
+unset xtics
+set ylabel "Shear friction [-]"
+
+plot "cowden.dat" w l t "Cowden till (Tika et al., 1996)", \
+ "lower_cromer.dat" w l t "Lower Cromer till (Tika et al., 1996)", \
+ "ceasar.dat" w l t "Ceasar till (Rathbun et al., 2008)", \
+ "storglaciaren.dat" w l t "Storglaciaren till (Iverson et al., 1998)", \
+ "two_rivers.dat" w l t "Two Rivers till (Iverson et al., 1998)", \
+ "des_moines_lobe.dat" w l t "Des Moines Lobe till (Iverson, 2010)", \
+ "wis.dat" w l t "WIS till (Tulaczyk et al., 2001)"
+
+set xtics
+set xlabel "Shear strain rate [a^{-1}]"
+set ylabel "Shear friction [-]"
+
load 'BlueSeq.plt'
plot "out_b0.94.txt" u ($2*3600*24*365):1 w l lc '#666666' lw 2 dashtype '-' t "glass beads, b = 0.94", \
diff --git a/fig1/lower_cromer.dat b/fig1/lower_cromer.dat
@@ -0,0 +1,2 @@
+5.0e1 0.60
+7.0e4 0.56
diff --git a/fig1/storglaciaren.dat b/fig1/storglaciaren.dat
@@ -0,0 +1,2 @@
+5.3e2 0.57
+2.3e4 0.52
diff --git a/fig1/two_rivers.dat b/fig1/two_rivers.dat
@@ -0,0 +1,2 @@
+1.8e3 0.475
+8.0e4 0.44
diff --git a/fig1/wis.dat b/fig1/wis.dat
@@ -0,0 +1,2 @@
+1.0e1 0.504
+1.0e5 0.495
diff --git a/iverson2010-fig2a-data.txt b/iverson2010-fig2a-data.txt
@@ -1,14 +0,0 @@
-3.2e3 0.52 Cowden till (Tika et al., 1996)
-5.0e5 0.52 Cowden till (Tika et al., 1996)
-5.0e1 0.60 Lower Cromer till (Tika et al., 1996)
-7.0e4 0.56 Lower Cromer till (Tika et al., 1996)
-3.0e3 0.53 Ceasar till (Rathbun et al., 2008)
-9.0e5 0.58 Ceasar till (Rathbun et al., 2008)
-5.2e2 0.57 Storglaciaren till (Iverson et al., 1998)
-2.3e4 0.52 Storglaciaren till (Iverson et al., 1998)
-1.9e3 0.48 Two Rivers till (Iverson et al., 1998)
-8.0e4 0.45 Two Rivers till (Iverson et al., 1998)
-1.0e3 0.50 Des Moines Lobe till (Iverson, 2010)
-4.0e4 0.46 Des Moines Lobe till (Iverson, 2010)
-1.0e1 0.51 WIS till (Tulaczyk et al., 2001)
-1.0e5 0.49 WIS till (Tulaczyk et al., 2001)
