1d_fd_simple_shear_rheology_tulaczyk.gp (1038B)
1 #!/usr/bin/env gnuplot 2 3 ## resemble Tulaczyk 2006 plot: 4 5 set terminal pngcairo color size 18.6 cm, 17.6 cm 6 set output "1d_fd_simple_shear_rheology_tulaczyk.png" 7 set xlabel "Normalized shear stress [-]" 8 set ylabel "Normalized shear strain rate [-]" offset 2 9 10 set xrange [0.1:10.0] 11 set yrange [1.0:1000.0] 12 set key top left #samplen 0.9 13 14 set logscale xy 15 16 load 'BlueSeq.plt' 17 18 # divide friction by threshold friction value (mu_s) 19 # divide shear velocity by 1 m/h 20 21 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 "b = 0.94", \ 22 "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", \ 23 "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", \ 24 "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", \ 25 "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"