commit 360884f20da7d8519ba7e021361855d07e182242
parent 80d49becbd422aa4e103662090bca17fe19652eb
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date: Sat, 11 May 2019 09:51:52 -0700
Change color sequence and include glass-bead run
Diffstat:
7 files changed, 27 insertions(+), 31 deletions(-)
diff --git a/1d_fd_simple_shear_rheology_iverson.gp b/1d_fd_simple_shear_rheology_iverson.gp
@@ -1,9 +1,13 @@
#!/usr/bin/env gnuplot
+reset
+
## resemble Iverson 2010 plot:
-set terminal pngcairo color size 18.6 cm, 11.0 cm
+#set terminal pngcairo color size 18.6 cm, 11.0 cm
+set terminal pngcairo 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
@@ -13,14 +17,10 @@ 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.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"
+load 'BlueSeq.plt'
-set xtics norotate # Restore defaults
+plot "1d_fd_simple_shear_rheology_b0.94_iverson.txt" u ($2*3600*24*365):1 w l lc '#666666' lw 2 dashtype '-' t "glass beads, b = 0.94", \
+ "1d_fd_simple_shear_rheology_b0.40_iverson.txt" u ($2*3600*24*365):1 w l ls 6 lw 2 t "b = 0.40", \
+ "1d_fd_simple_shear_rheology_b0.20_iverson.txt" u ($2*3600*24*365):1 w l ls 5 lw 2 t "b = 0.20", \
+ "1d_fd_simple_shear_rheology_b0.10_iverson.txt" u ($2*3600*24*365):1 w l ls 3 lw 2 t "b = 0.10", \
+ "1d_fd_simple_shear_rheology_b0.01_iverson.txt" u ($2*3600*24*365):1 w l ls 1 lw 2 t "b = 0.01"
diff --git a/1d_fd_simple_shear_rheology_iverson.png b/1d_fd_simple_shear_rheology_iverson.png
Binary files differ.
diff --git a/1d_fd_simple_shear_rheology_kamb.gp b/1d_fd_simple_shear_rheology_kamb.gp
@@ -11,14 +11,12 @@ set xrange [0.0:0.028]
set yrange [0.0:2.5]
set key top left #samplen 0.9
-set style line 1 linetype 1 linewidth 3 pointtype 1 pointsize 1
+load 'BlueSeq.plt'
-plot "1d_fd_simple_shear_rheology_b0.01_kamb.txt" u ($1*1700/100000):($2*3600*24) w l lw 2 t "b = 0.01", \
- "1d_fd_simple_shear_rheology_b0.10_kamb.txt" u ($1*1700/100000):($2*3600*24) w l lw 2 t "b = 0.1", \
- "1d_fd_simple_shear_rheology_b0.20_kamb.txt" u ($1*1700/100000):($2*3600*24) w l lw 2 t "b = 0.2", \
- "1d_fd_simple_shear_rheology_b0.40_kamb.txt" u ($1*1700/100000):($2*3600*24) w l lw 2 t "b = 0.4", \
- "1d_fd_simple_shear_rheology_b0.60_kamb.txt" u ($1*1700/100000):($2*3600*24) w l lw 2 t "b = 0.6", \
- "1d_fd_simple_shear_rheology_b0.80_kamb.txt" u ($1*1700/100000):($2*3600*24) w l lw 2 t "b = 0.8", \
- "1d_fd_simple_shear_rheology_b1.00_kamb.txt" u ($1*1700/100000):($2*3600*24) w l lw 2 t "b = 1.0"
+plot "1d_fd_simple_shear_rheology_b0.01_kamb.txt" u ($1*1700/100000):($2*3600*24) w l ls 1 lw 2 t "b = 0.01", \
+ "1d_fd_simple_shear_rheology_b0.10_kamb.txt" u ($1*1700/100000):($2*3600*24) w l ls 2 lw 2 t "b = 0.1", \
+ "1d_fd_simple_shear_rheology_b0.20_kamb.txt" u ($1*1700/100000):($2*3600*24) w l ls 3 lw 2 t "b = 0.2", \
+ "1d_fd_simple_shear_rheology_b0.40_kamb.txt" u ($1*1700/100000):($2*3600*24) w l ls 4 lw 2 t "b = 0.4", \
+ "1d_fd_simple_shear_rheology_b0.94_kamb.txt" u ($1*1700/100000):($2*3600*24) w l ls 5 lw 2 t "b = 0.94"
set xtics norotate # Restore defaults
diff --git a/1d_fd_simple_shear_rheology_kamb.png b/1d_fd_simple_shear_rheology_kamb.png
Binary files differ.
diff --git a/1d_fd_simple_shear_rheology_tulaczyk.gp b/1d_fd_simple_shear_rheology_tulaczyk.gp
@@ -13,16 +13,14 @@ set key top left #samplen 0.9
set logscale xy
-set style line 1 linetype 1 linewidth 3 pointtype 1 pointsize 1
+load 'BlueSeq.plt'
# divide friction by threshold friction value (mu_s)
# divide shear velocity by 1 m/h
-plot "1d_fd_simple_shear_rheology_b0.01_tulaczyk.txt" u ($1/0.502):($2/3.5e-7) w l lw 2 t "b = 0.01", \
- "1d_fd_simple_shear_rheology_b0.10_tulaczyk.txt" u ($1/0.502):($2/3.5e-8) w l lw 2 t "b = 0.1", \
- "1d_fd_simple_shear_rheology_b0.20_tulaczyk.txt" u ($1/0.502):($2/1.7e-8) w l lw 2 t "b = 0.2", \
- "1d_fd_simple_shear_rheology_b0.40_tulaczyk.txt" u ($1/0.502):($2/8.6e-9) w l lw 2 t "b = 0.4", \
- "1d_fd_simple_shear_rheology_b0.60_tulaczyk.txt" u ($1/0.502):($2/5.7e-9) w l lw 2 t "b = 0.6", \
- "1d_fd_simple_shear_rheology_b0.80_tulaczyk.txt" u ($1/0.502):($2/4.3e-9) w l lw 2 t "b = 0.8", \
- "1d_fd_simple_shear_rheology_b1.00_tulaczyk.txt" u ($1/0.502):($2/3.4e-9) w l lw 2 t "b = 1.0"
+plot "1d_fd_simple_shear_rheology_b0.01_tulaczyk.txt" u ($1/0.502):($2/3.5e-7) w l ls 1 lw 2 t "b = 0.01", \
+ "1d_fd_simple_shear_rheology_b0.10_tulaczyk.txt" u ($1/0.502):($2/3.5e-8) w l ls 2 lw 2 t "b = 0.1", \
+ "1d_fd_simple_shear_rheology_b0.20_tulaczyk.txt" u ($1/0.502):($2/1.7e-8) w l ls 3 lw 2 t "b = 0.2", \
+ "1d_fd_simple_shear_rheology_b0.40_tulaczyk.txt" u ($1/0.502):($2/8.6e-9) w l ls 4 lw 2 t "b = 0.4", \
+ "1d_fd_simple_shear_rheology_b0.94_tulaczyk.txt" u ($1/0.502):($2/3.4e-9) w l lc 'black' lw 2 t "b = 1.0"
set xtics norotate # Restore defaults
diff --git a/1d_fd_simple_shear_rheology_tulaczyk.png b/1d_fd_simple_shear_rheology_tulaczyk.png
Binary files differ.
diff --git a/Makefile b/Makefile
@@ -66,10 +66,10 @@ diurnal.gif: diurnal.mp4
# 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 '\
- for b in $$(printf "0.01\n0.10\n"; seq 0.20 0.20 1.00); do \
+ for b in $$(printf "0.01\n0.10\n"; seq 0.20 0.20 1.00) 0.94; do \
out="$<_rheology_b$${b}_kamb.txt"; \
rm -f "$$out"; \
- for t in $$(seq 0.01 0.01 2.0); do \
+ for t in $$(seq 0.01 0.01 2.0) 0.94; do \
printf "$$t\t" >> "$$out"; \
./$< -f 1.1 -P 1.7e3 --stress-ratio $$t -b $$b | \
tail -n 1 | cut -f2 >> "$$out"; \
@@ -80,7 +80,7 @@ diurnal.gif: diurnal.mp4
# 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\n0.10\n"; seq 0.20 0.20 1.00); do \
+ for b in $$(printf "0.01\n0.10\n"; seq 0.20 0.20 0.90) 0.94; do \
out="$<_rheology_b$${b}_iverson.txt"; \
rm -f "$$out"; \
for t in $$(seq 0.0001 0.002 1.0); do \
@@ -93,7 +93,7 @@ diurnal.gif: diurnal.mp4
# shear velocity rate from 0.1 m/h to 100 m/h
1d_fd_simple_shear_rheology_tulaczyk.png: 1d_fd_simple_shear 1d_fd_simple_shear_rheology_tulaczyk.gp
/bin/bash -c '\
- for b in $$(printf "0.01\n0.10\n"; seq 0.20 0.20 1.00); do \
+ for b in $$(printf "0.01\n0.10\n"; seq 0.20 0.20 1.00) 0.94; do \
out="$<_rheology_b$${b}_tulaczyk.txt"; \
rm -f "$$out"; \
for t in $$(seq 0.1 0.002 0.9); do \
