fig.gp (1743B)
1 #!/usr/bin/env gnuplot 2 3 reset 4 5 ## resemble Iverson 2010 fig2a 6 7 #set terminal pdfcairo size 4.45 cm, 4.45 cm font ",10" 8 set terminal pdfcairo enhanced color size 7 cm, 4.5 cm font ",10" 9 10 set margins 8,2,3.2,1 11 12 set colorsequence podo 13 14 unset xrange 15 unset yrange 16 unset logscale x 17 set xtics 50 18 set ytics 25 19 set ylabel "Shear stress {/Symbol t} [kPa]" offset 1.0 20 set xlabel "Eff. normal stress {/Symbol s}_n' [kPa]" 21 set xrange [0:200] 22 set yrange [0:120] 23 24 set angles degrees 25 caesar(x) = (x >= 50 && x <= 5000) ? tan(28)*x - 6.5 : 1/0 # Rathbun et al., 2008, saturated Caesar till 26 storglaciaren(x) = (x >= 25 && x <= 150) ? tan(26.3)*x + 5 : 1/0 # Iverson et al., 1998 27 two_rivers(x) = (x >= 30 && x <= 145) ? tan(17.8)*x + 14 : 1/0 # Iverson et al., 1998 28 # Des Moines Lobe till properties extracted from Iverson 2010 plot, mc_des_moines_lobe.dat 29 wis(x) = (x >= 20 && x <= 310) ? tan(24)*x + 3 : 1/0 # Tulaczyk et al., 2000 30 31 #set key bottom right font ",9" samplen 1.5 maxrows 4 32 set key top left reverse Left font ",8" samplen 1.5 maxrows 4 width -3 33 34 set label "c" at screen 0.02,0.95 font ",12" 35 36 plot \ 37 caesar(x) w l lc 1 lw 2 t "Caesar till", \ 38 "mc_caesar.txt" u ($1):($2*$1 - 6.5) w l dashtype 2 lc 1 lw 2 t "{/:Bold CNGF-PF}", \ 39 storglaciaren(x) w l lc 2 lw 2 t "Storglaciaren till", \ 40 "mc_storglaciaren.txt" u ($1):($2*$1 + 5) w l dashtype 2 lc 2 lw 2 t "{/:Bold CNGF-PF}", \ 41 two_rivers(x) w l lc 3 lw 2 t "Two Rivers till", \ 42 "mc_two_rivers.txt" u ($1):($2*$1 + 14) w l dashtype 2 lc 3 lw 2 t "{/:Bold CNGF-PF}", \ 43 wis(x) w l lc 4 lw 2 t "WIS till", \ 44 "mc_wis.txt" u ($1):($2*$1 + 3) w l dashtype 2 lc 4 lw 2 t "{/:Bold CNGF-PF}"