cngf-pf

commit 60cf1cd6cf19df4885446386216b44889551ef45
parent 78a6e567a38a425ae78900fe60689e2d6103488a
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Thu, 11 Apr 2019 20:25:04 +0200

Add comparison to Iverson 2010

Diffstat:
A
1d_fd_simple_shear_rheology_iverson.gp
 | 
25
+++++++++++++++++++++++++
A
1d_fd_simple_shear_rheology_iverson.png
 | 
0

M
Makefile
 | 
18
++++++++++++++++--
M
README.md
 | 
14
+++++++++++++-
A
iverson2010-fig2a.png
 | 
0

A
iverson2010-fig2b.png
 | 
0

6 files changed, 54 insertions(+), 3 deletions(-)
diff --git a/1d_fd_simple_shear_rheology_iverson.gp b/1d_fd_simple_shear_rheology_iverson.gp
@@ -0,0 +1,25 @@
+#!/usr/bin/env gnuplot
+
+## resemble Iverson 210 plot:
+
+set terminal pngcairo color size 18.6 cm, 11.0 cm
+set output "1d_fd_simple_shear_rheology_iverson.png"
+set xlabel "Shear strain rate [a^{-1}]"
+set ylabel "Shear friction [-]" offset 2
+
+set xrange [1e1:1e6]  # m/a
+set yrange [0.0:0.8]
+set key bottom right #samplen 0.9
+
+set logscale x
+
+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"
+
+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
@@ -6,7 +6,8 @@ HDR=$(wildcard *.h)
 
 default: 1d_fd_simple_shear.png \
 	1d_fd_simple_shear_rheology.png \
-	1d_fd_simple_shear_rheology_kamb.png
+	1d_fd_simple_shear_rheology_kamb.png \
+	1d_fd_simple_shear_rheology_iverson.png
 
 1d_fd_simple_shear: $(OBJ) $(HDR)
 	$(CC) $(LDFLAGS) $(OBJ) -o $@
@@ -32,7 +33,6 @@ default: 1d_fd_simple_shear.png \
 	done; done'
 	gnuplot $<_rheology.gp > $@
 
-
 # 1 bar is equal to 100 kPa
 1d_fd_simple_shear_rheology_kamb.png: 1d_fd_simple_shear 1d_fd_simple_shear_rheology_kamb.gp
 	/bin/bash -c '\
@@ -46,6 +46,20 @@ default: 1d_fd_simple_shear.png \
 	done; done'
 	gnuplot $<_rheology_kamb.gp > $@
 
+# shear-strain rate from 10^1 to 10^6 m/a
+# 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 \
+		out="$<_rheology_b$${b}_iverson.txt"; \
+		rm -f "$$out"; \
+		for t in $$(seq 0.0001 0.002 1.0); do \
+			printf "$$t\t" >> "$$out"; \
+			./$< -f 0.55 -P 100e3 -L 1.0 --stress-ratio $$t -b $$b | \
+			tail -n 1 | cut -f2 >> "$$out"; \
+	done; done'
+	gnuplot $<_rheology_iverson.gp > $@
+
 .PHONY: watch
 watch:
 	echo $(SRC) $(HDR) | tr ' ' '\n' | entr -s 'make && ./1d_fd_simple_shear'
diff --git a/README.md b/README.md
@@ -37,10 +37,22 @@ Glass beads have *b* = 0.94.
 
 ![1d_fd_simple_shear_rheology.png](https://gitlab.com/admesg/1d_fd_simple_shear/raw/master/1d_fd_simple_shear_rheology.png)
 
-Comparison to Kamb 1991:
+---
+
+#### Continuum model
 
 ![1d_fd_simple_shear_rheology_kamb.png](https://gitlab.com/admesg/1d_fd_simple_shear/raw/master/1d_fd_simple_shear_rheology_kamb.png)
 
 #### Kamb 1991, subglacial till from Upstream-B site
 
 ![kamb1991-fig1.png](https://gitlab.com/admesg/1d_fd_simple_shear/raw/master/kamb1991-fig1.png)
+
+---
+
+#### Continuum model
+
+![1d_fd_simple_shear_rheology_iverson.png](https://gitlab.com/admesg/1d_fd_simple_shear/raw/master/1d_fd_simple_shear_rheology_iverson.png)
+
+#### Iverson 2010, various subglacial tills
+
+![iverson2010-fig2a.png](https://gitlab.com/admesg/1d_fd_simple_shear/raw/master/iverson2010-fig2a.png)
diff --git a/iverson2010-fig2a.png b/iverson2010-fig2a.png
Binary files differ.
diff --git a/iverson2010-fig2b.png b/iverson2010-fig2b.png
Binary files differ.