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 eb3be4b7b0e306c142e9ceacb2dc59e35dad96de
parent 103219aff8e6a91f510755e753f3efac5367d2af
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Tue, 25 Jun 2019 14:46:53 +0200

Fix iterative plotting

Diffstat:
Mfig3/fig.gp | 59+++++++++++++++++++++++------------------------------------
1 file changed, 23 insertions(+), 36 deletions(-)

diff --git a/fig3/fig.gp b/fig3/fig.gp @@ -8,41 +8,28 @@ set yrange [0.0:2.0] set key bottom right #samplen 0.9 -do for [t=47:71] { - filename = sprintf('unnamed.output%05d.txt', t) - if (t == 71) { - set xlabel "Water pressure [kPa]" - set ylabel "Vertical position, z [m]" - } else { - set xlabel " " - set ylabel " " - } - set xrange [0:100] - plot filename u ($4/1e3):1 w l lw 1 lc "blue" t "" - - if (t == 71) { - set xlabel "Effective normal stress [kPa]" - set ylabel "" - } else { - set xlabel " " - set ylabel "" - } - set xrange [50:150] - plot filename u ($3/1e3):1 w l lw 1 lc "black" t "" - - # set xlabel "Friction [-]" - # set xrange [-0.05:1.5] - # plot filename u 5:1 w l lw 1 lc "gray" t "" - - if (t == 71) { - set xlabel "Shear velocity [m/s]" - set ylabel "" - } else { - set xlabel " " - set ylabel "" - } - set xrange [0.0:0.023] - plot filename u 2:1 w l lw 1 lc "red" t "" -} +t0 = 47 +tend = 71 + +filename(t) = sprintf('unnamed.output%05d.txt', t) + +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 "" + +set xlabel "Effective normal stress [kPa]" +set ylabel "" +set xrange [50:150] +plot for [t = t0:tend] filename(t) u ($3/1e3):1 w l lw 1 lc "black" t "" + +# set xlabel "Friction [-]" +# set xrange [-0.05:1.5] +# plot filename u 5:1 w l lw 1 lc "gray" t "" + +set xlabel "Shear velocity [m/s]" +set ylabel "" +set xrange [0.0:0.023] +plot for [t = t0:tend] filename(t) u 2:1 w l lw 1 lc "red" t "" unset multiplot