tfig.gp - cngf-pf-exp1 - experiments for first paper with continuum granular model
 (HTM) git clone git://src.adamsgaard.dk/manus_continuum_granular1_exp
 (DIR) Log
 (DIR) Files
 (DIR) Refs
 (DIR) Submodules
 (DIR) README
 (DIR) LICENSE
       ---
       tfig.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}"