manus_continuum_granular1_exp

Experiments for first paper with continuum granular model
git clone git://src.adamsgaard.dk/.manus_continuum_granular1_exp
Log | Files | Refs | Submodules | README | LICENSE

commit 5c69b9d48594dc62fd737076595be068efd19e64
parent 3b5829e907d5977e03091ffc736399aa02e38b0d
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Tue, 25 Jun 2019 15:50:13 +0200

Update line colors for figure 3

Diffstat:
Mfig3/fig.gp | 22+++++++++++++++++-----
1 file changed, 17 insertions(+), 5 deletions(-)

diff --git a/fig3/fig.gp b/fig3/fig.gp @@ -1,5 +1,7 @@ #!/usr/bin/env gnuplot +reset + #set terminal pdfcairo color size 7.5 cm, 7.5 cm set terminal pdfcairo color size 15.0 cm, 5.0 cm set multiplot layout 1,4 \ @@ -11,27 +13,37 @@ set yrange [0.0:2.0] set key bottom right #samplen 0.9 t0 = 47 -tend = 71 +tend = 70 filename(t) = sprintf('unnamed.output%05d.txt', t) +#set palette color +set pm3d map +#set palette rainbow +#set palette rgbformulae 7,5,15 +#set palette rgb 29,30,31 +set palette rgb 10,6,6 +#set palette model RGB rgbformulae 7,5,15 +#set palette functions gray, gray, gray + +unset colorbox set xlabel "Water pressure [kPa]" set ylabel "Vertical position, z [m]" set xrange [0:100] -plot for [t = t0:tend] filename(t) u ($4/1e3):1 w l lw 1 lc "blue" t "" +plot for [t = t0:tend] filename(t) u ($4/1e3):1:(t-t0) w l lc palette lw 1 t "" set xlabel "Eff. normal stress [kPa]" set ylabel "" set xrange [50:150] unset ytics -plot for [t = t0:tend] filename(t) u ($3/1e3):1 w l lw 1 lc "black" t "" +plot for [t = t0:tend] filename(t) u ($3/1e3):1:(t-t0) w l lc palette lw 1 t "" set xlabel "Shear velocity [m/s]" set ylabel "" set xrange [0.0:0.023] set xtics (0,0.01,0.02) unset ytics -plot for [t = t0:tend] filename(t) u 2:1 w l lw 1 lc "red" t "" +plot for [t = t0:tend] filename(t) u 2:1:(t-t0) w l lc palette lw 1 t "" set xlabel "Shear strain rate [1/s]" set ylabel "" @@ -41,7 +53,7 @@ set xtics (0,0.001,0.002,0.003) unset ytics prev1 = prev2 = 0 shift(x) = (prev2 = prev1, prev1 = x) -plot for [t = t0:tend] filename(t) u (shift($2), $0 < 1 ? 1/0 : $2 - prev2):1 w l lw 1 lc "magenta" t "" +plot for [t = t0:tend] filename(t) u (shift($2), $0 < 1 ? 1/0 : $2 - prev2):1:(t-t0) w l lc palette lw 1 t "" unset format x unset multiplot