tPlot state at 1h intervals - cngf-pf-exp1 - experiments for first paper with continuum granular model
(HTM) git clone git://src.adamsgaard.dk/manus_continuum_granular1_exp
(DIR) Log
(DIR) Files
(DIR) Refs
(DIR) Submodules
(DIR) README
(DIR) LICENSE
---
(DIR) commit bdd1c0a18b64d7f774b9664c832881f2dc1e436f
(DIR) parent 8d0714ab70fd06552049eaa83855ccfc6a92ca62
(HTM) Author: Anders Damsgaard <anders@adamsgaard.dk>
Date: Wed, 18 Dec 2019 15:42:07 +0100
Plot state at 1h intervals
Diffstat:
M fig-stick_slip_rate_depth/fig.gp | 12 +++++++-----
1 file changed, 7 insertions(+), 5 deletions(-)
---
(DIR) diff --git a/fig-stick_slip_rate_depth/fig.gp b/fig-stick_slip_rate_depth/fig.gp
t@@ -12,7 +12,6 @@ set key bottom right #samplen 0.9
t0 = 47
tend = 71
-#tend = 140
filename(t) = sprintf('sim.output%05d.txt', t)
t@@ -24,6 +23,9 @@ set palette rgb 10,6,6
#f = 1.0/(3600.0*24.0)
#f(x) = 8.0 - skindepth(f,k)
+file_interval=60*10
+plot_interval=60*60
+inc=plot_interval/file_interval
set label "{/:Bold a}" at graph 0.90,0.93
unset colorbox
t@@ -32,7 +34,7 @@ set ylabel "{/:Normal Vertical position [m]}"
set xrange [0:200]
set arrow from 0,skin_depth to 200,skin_depth nohead lc black front
set arrow from 0,max_def_depth to 200,max_def_depth nohead front lc black dashtype 2
-plot for [t = t0:tend] filename(t) u ($4/1e3):1:(t-t0) w l lc palette lw 1 t ""
+plot for [t = t0:tend:inc] filename(t) u ($4/1e3):1:(t-t0) w l lc palette lw 1 t ""
unset label
set label "{/:Bold b}" at graph 0.90,0.93
t@@ -40,7 +42,7 @@ set xlabel "{/:Normal Eff. normal stress [kPa]}"
set ylabel ""
set xrange [0:200]
unset ytics
-plot for [t = t0:tend] filename(t) u ($3/1e3):1:(t-t0) w l lc palette lw 1 t ""
+plot for [t = t0:tend:inc] filename(t) u ($3/1e3):1:(t-t0) w l lc palette lw 1 t ""
unset label
set label "{/:Bold c}" at graph 0.90,0.93
t@@ -53,7 +55,7 @@ unset ytics
unset arrow
set arrow from 0,skin_depth to 1,skin_depth nohead lc black front
set arrow from 0,max_def_depth to 1,max_def_depth nohead front lc black dashtype 2
-plot for [t = t0:tend] filename(t) u ($2*60*60*24*365/1000):1:(t-t0) w l lc palette lw 1 t ""
+plot for [t = t0:tend:inc] filename(t) u ($2*60*60*24*365/1000):1:(t-t0) w l lc palette lw 1 t ""
unset label
set label "{/:Bold d}" at graph 0.90,0.93
t@@ -64,7 +66,7 @@ unset ytics
unset arrow
set arrow from 0,skin_depth to 20,skin_depth nohead lc black front
set arrow from 0,max_def_depth to 20,max_def_depth nohead front lc black dashtype 2
-plot for [t = t0:tend] filename(t) u ($6*60*60*24):1:(t-t0) w l lc palette lw 1 t ""
+plot for [t = t0:tend:inc] filename(t) u ($6*60*60*24):1:(t-t0) w l lc palette lw 1 t ""
unset format x
unset multiplot