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