1d_fd_simple_shear

commit 426077b91107ec1f7c67880a0fc3812045adeaab
parent b3e24e64943968daacc934d3549a983991b95207
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Thu, 11 Apr 2019 17:10:58 +0200

Add missing gnuplot script

Diffstat:
A
1d_fd_simple_shear_rheology_kamb.gp
 | 
24
++++++++++++++++++++++++

1 file changed, 24 insertions(+), 0 deletions(-)
diff --git a/1d_fd_simple_shear_rheology_kamb.gp b/1d_fd_simple_shear_rheology_kamb.gp
@@ -0,0 +1,24 @@
+#!/usr/bin/env gnuplot
+
+## resemble Kamb 1991 plot:
+
+set terminal pngcairo color size 18.6 cm, 17.6 cm
+set output "1d_fd_simple_shear_rheology_kamb.png"
+set xlabel "Shear stress [bar]"
+set ylabel "Shear displacment rate [m/day]" offset 2
+
+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
+
+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"
+
+set xtics norotate # Restore defaults