1d_fd_simple_shear_rheology_kamb.gp (934B)
1 #!/usr/bin/env gnuplot 2 reset 3 4 ## resemble Kamb 1991 plot: 5 6 set terminal pngcairo color size 18.6 cm, 17.6 cm 7 set output "1d_fd_simple_shear_rheology_kamb.png" 8 set xlabel "Shear stress [bar]" 9 set ylabel "Shear displacment rate [m/day]" offset 2 10 11 set xrange [0.0:0.028] 12 set yrange [0.0:2.5] 13 set key top left #samplen 0.9 14 15 load 'BlueSeq.plt' 16 17 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", \ 18 "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", \ 19 "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", \ 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", \ 21 "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"