commit db546c09dc4b80b452d7e84a7c31df13d7106725
parent 13a67b505eeb9ac9514d31c9f3a4f9ed9b1c30ec
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date: Thu, 11 Apr 2019 20:40:26 +0200
Use consistent b values for Iverson comparison
Diffstat:
3 files changed, 7 insertions(+), 6 deletions(-)
diff --git a/1d_fd_simple_shear_rheology_iverson.gp b/1d_fd_simple_shear_rheology_iverson.gp
@@ -1,6 +1,6 @@
#!/usr/bin/env gnuplot
-## resemble Iverson 210 plot:
+## resemble Iverson 2010 plot:
set terminal pngcairo color size 18.6 cm, 11.0 cm
set output "1d_fd_simple_shear_rheology_iverson.png"
@@ -17,9 +17,10 @@ set style line 1 linetype 1 linewidth 3 pointtype 1 pointsize 1
plot "1d_fd_simple_shear_rheology_b0.01_iverson.txt" u ($2*3600*24*365):1 w l lw 2 t "b = 0.01", \
"1d_fd_simple_shear_rheology_b0.10_iverson.txt" u ($2*3600*24*365):1 w l lw 2 t "b = 0.1", \
- "1d_fd_simple_shear_rheology_b0.30_iverson.txt" u ($2*3600*24*365):1 w l lw 2 t "b = 0.3", \
- "1d_fd_simple_shear_rheology_b0.50_iverson.txt" u ($2*3600*24*365):1 w l lw 2 t "b = 0.5", \
- "1d_fd_simple_shear_rheology_b0.70_iverson.txt" u ($2*3600*24*365):1 w l lw 2 t "b = 0.7", \
- "1d_fd_simple_shear_rheology_b0.90_iverson.txt" u ($2*3600*24*365):1 w l lw 2 t "b = 0.9"
+ "1d_fd_simple_shear_rheology_b0.20_iverson.txt" u ($2*3600*24*365):1 w l lw 2 t "b = 0.2", \
+ "1d_fd_simple_shear_rheology_b0.40_iverson.txt" u ($2*3600*24*365):1 w l lw 2 t "b = 0.4", \
+ "1d_fd_simple_shear_rheology_b0.60_iverson.txt" u ($2*3600*24*365):1 w l lw 2 t "b = 0.6", \
+ "1d_fd_simple_shear_rheology_b0.80_iverson.txt" u ($2*3600*24*365):1 w l lw 2 t "b = 0.8", \
+ "1d_fd_simple_shear_rheology_b1.00_iverson.txt" u ($2*3600*24*365):1 w l lw 2 t "b = 1.0"
set xtics norotate # Restore defaults
diff --git a/1d_fd_simple_shear_rheology_iverson.png b/1d_fd_simple_shear_rheology_iverson.png
Binary files differ.
diff --git a/Makefile b/Makefile
@@ -50,7 +50,7 @@ default: 1d_fd_simple_shear.png \
# friction around 0.55
1d_fd_simple_shear_rheology_iverson.png: 1d_fd_simple_shear 1d_fd_simple_shear_rheology_iverson.gp
/bin/bash -c '\
- for b in $$(printf "0.01\n"; seq 0.10 0.20 0.90); do \
+ for b in $$(printf "0.01\n0.10\n"; seq 0.20 0.20 1.00); do \
out="$<_rheology_b$${b}_iverson.txt"; \
rm -f "$$out"; \
for t in $$(seq 0.0001 0.002 1.0); do \