1d_fd_simple_shear

continuum model for granular flows with pore-pressure dynamics
git clone git://src.adamsgaard.dk/1d_fd_simple_shear
Log | Files | Refs | README | LICENSE Back to index

1d_fd_simple_shear_rheology_iverson.gp (981B)

      1 #!/usr/bin/env gnuplot
      2 
      3 reset
      4 
      5 ## resemble Iverson 2010 plot:
      6 
      7 #set terminal pngcairo color size 18.6 cm, 11.0 cm
      8 set terminal pngcairo size 18.6 cm, 11.0 cm
      9 set output "1d_fd_simple_shear_rheology_iverson.png"
     10 
     11 set xlabel "Shear strain rate [a^{-1}]"
     12 set ylabel "Shear friction [-]" offset 2
     13 
     14 set xrange [1e1:1e6]  # m/a
     15 set yrange [0.0:0.8]
     16 set key bottom right #samplen 0.9
     17 
     18 set logscale x
     19 
     20 load 'BlueSeq.plt'
     21 
     22 plot "1d_fd_simple_shear_rheology_b0.94_iverson.txt" u ($2*3600*24*365):1 w l lc '#666666' lw 2 dashtype '-' t "glass beads, b = 0.94", \
     23          "1d_fd_simple_shear_rheology_b0.40_iverson.txt" u ($2*3600*24*365):1 w l ls 6 lw 2 t "b = 0.40", \
     24          "1d_fd_simple_shear_rheology_b0.20_iverson.txt" u ($2*3600*24*365):1 w l ls 5 lw 2 t "b = 0.20", \
     25          "1d_fd_simple_shear_rheology_b0.10_iverson.txt" u ($2*3600*24*365):1 w l ls 3 lw 2 t "b = 0.10", \
     26          "1d_fd_simple_shear_rheology_b0.01_iverson.txt" u ($2*3600*24*365):1 w l ls 1 lw 2 t "b = 0.01"