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_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"