tUse consistent b values for Iverson comparison - cngf-pf - continuum model for granular flows with pore-pressure dynamics (renamed from 1d_fd_simple_shear)
HTML git clone git://src.adamsgaard.dk/cngf-pf
DIR Log
DIR Files
DIR Refs
DIR README
DIR LICENSE
---
DIR commit db546c09dc4b80b452d7e84a7c31df13d7106725
DIR parent 13a67b505eeb9ac9514d31c9f3a4f9ed9b1c30ec
HTML Author: Anders Damsgaard <anders@adamsgaard.dk>
Date: Thu, 11 Apr 2019 20:40:26 +0200
Use consistent b values for Iverson comparison
Diffstat:
M 1d_fd_simple_shear_rheology_iverso… | 11 ++++++-----
M 1d_fd_simple_shear_rheology_iverso… | 0
M Makefile | 2 +-
3 files changed, 7 insertions(+), 6 deletions(-)
---
DIR diff --git a/1d_fd_simple_shear_rheology_iverson.gp b/1d_fd_simple_shear_rheology_iverson.gp
t@@ -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"
t@@ -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
DIR diff --git a/1d_fd_simple_shear_rheology_iverson.png b/1d_fd_simple_shear_rheology_iverson.png
Binary files differ.
DIR diff --git a/Makefile b/Makefile
t@@ -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 \