commit fa482eaba45c62b6cfa37afdc35a733a07b58757
parent 360884f20da7d8519ba7e021361855d07e182242
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date: Sat, 11 May 2019 10:15:06 -0700
Update rheology plots for consistency
Diffstat:
5 files changed, 13 insertions(+), 15 deletions(-)
diff --git a/1d_fd_simple_shear_rheology_kamb.gp b/1d_fd_simple_shear_rheology_kamb.gp
@@ -1,4 +1,5 @@
#!/usr/bin/env gnuplot
+reset
## resemble Kamb 1991 plot:
@@ -13,10 +14,8 @@ set key top left #samplen 0.9
load 'BlueSeq.plt'
-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
+plot "1d_fd_simple_shear_rheology_b0.94_kamb.txt" u ($1*1700/100000):($2*3600*24) w l lc '#666666' lw 2 dashtype '-' t "b = 0.94", \
+ "1d_fd_simple_shear_rheology_b0.40_kamb.txt" u ($1*1700/100000):($2*3600*24) w l ls 6 lw 2 t "b = 0.40", \
+ "1d_fd_simple_shear_rheology_b0.20_kamb.txt" u ($1*1700/100000):($2*3600*24) w l ls 5 lw 2 t "b = 0.20", \
+ "1d_fd_simple_shear_rheology_b0.10_kamb.txt" u ($1*1700/100000):($2*3600*24) w l ls 3 lw 2 t "b = 0.10", \
+ "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"
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
@@ -17,10 +17,9 @@ 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 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
+plot "1d_fd_simple_shear_rheology_b0.94_tulaczyk.txt" u ($1/0.502):($2/3.4e-9) w l lc '#666666' lw 2 dashtype '-' t "glass beads, b = 0.94", \
+ "1d_fd_simple_shear_rheology_b0.40_tulaczyk.txt" u ($1/0.502):($2/8.6e-9) w l ls 6 lw 2 t "b = 0.40", \
+ "1d_fd_simple_shear_rheology_b0.20_tulaczyk.txt" u ($1/0.502):($2/1.7e-8) w l ls 5 lw 2 t "b = 0.20", \
+ "1d_fd_simple_shear_rheology_b0.10_tulaczyk.txt" u ($1/0.502):($2/3.5e-8) w l ls 3 lw 2 t "b = 0.10", \
+ "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"
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) 0.94; do \
+ for b in 0.01 0.10 0.20 0.40 0.94; do \
out="$<_rheology_b$${b}_kamb.txt"; \
rm -f "$$out"; \
- for t in $$(seq 0.01 0.01 2.0) 0.94; do \
+ for t in $$(seq 0.01 0.01 2.0); do \
printf "$$t\t" >> "$$out"; \
./$< -f 1.1 -P 1.7e3 --stress-ratio $$t -b $$b | \
tail -n 1 | cut -f2 >> "$$out"; \