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 Back to index

commit b859325b552d5dfa4a39a924ee47dd1089f6bb3b
parent a181459dfb92b27332b6b98501d480c665f3987c
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Tue, 25 Aug 2020 14:58:59 +0200

Reorder subplots and clean up script

Diffstat:
Mfig-stick_slip_depth/fig.gp | 69+++++++++++++++++++++++++--------------------------------------------
1 file changed, 25 insertions(+), 44 deletions(-)

diff --git a/fig-stick_slip_depth/fig.gp b/fig-stick_slip_depth/fig.gp @@ -23,35 +23,30 @@ set palette set ylabel "{/:Normal Depth [m]}" offset 0.5 -set origin -0.02, 0.50 -set size 0.36, 0.46 +set origin 0.06, 0.0 +set size 0.313, 0.46 xmin = 0.0 xmax = 200.0 set xrange [xmin:xmax] set arrow from 88,skin_depth to xmax,skin_depth nohead lc black front -#set arrow from xmin,max_def_depth to xmax,max_def_depth nohead front lc black dashtype 2 -set colorbox +set colorbox user origin 0.001, 0.09 size 0.015, 0.34 set cblabel "Time [d]" set cbrange [2:3] set cbtics 0.2 -set label "A" font ",11" at graph 0.00,1.1 +set label "D" font ",11" at graph 0.00,1.1 set xlabel "{/:Normal Effective normal stress [kPa]}" plot for [t = t0:tend:inc*2] filename(t) u ($3/1e3):(L_z - $1):(t*file_interval/(24.*3600.)) w l lc palette lw 2 t "" unset arrow unset colorbox -set origin 0.33, 0.50 -set size 0.33, 0.46 +set origin 0.373, 0.0 +set size 0.313, 0.46 unset label -set label "B (Speed controlled)" font ",11" at graph 0.00,1.1 +set label "E (Speed controlled)" font ",11" at graph 0.00,1.1 set xlabel "{/:Normal Shear strain rate [1/d]}" font ",10" -#set ylabel "" -#unset ytics xmin = -2.0 xmax = 50.0 set xrange [xmin:xmax] -#set xtics 5 -#set arrow from xmin,skin_depth to xmax,skin_depth nohead lc black front set arrow from xmin,max_def_depth to 15,max_def_depth nohead front lc black dashtype 2 plot for [t = t0:tend:inc*2] filename(t) u ($6*60*60*24):(L_z - $1):(t*file_interval/(24.*3600.)) w l lc palette lw 2 t "" unset format x @@ -59,42 +54,36 @@ unset arrow filename(t) = sprintf('stress/sim.output%05d.txt', t) -set origin 0.67, 0.50 -set size 0.33, 0.46 +set origin 0.686, 0.0 +set size 0.313, 0.46 unset label -set label "C (Stress controlled)" font ",11" at graph 0.00,1.1 +set label "F (Stress controlled)" font ",11" at graph 0.00,1.1 set xlabel "{/:Normal Shear strain rate [1/min]}" -#set ylabel "" -#unset ytics xmin = -1.0 xmax = 12.0 set xrange [xmin:xmax] -#set arrow from xmin,skin_depth to xmax,skin_depth nohead lc black front -#set arrow from xmin,max_def_depth to xmax,max_def_depth nohead front lc black dashtype 2 plot for [t = t0:tend:inc*2] filename(t) u ($6*60):(L_z - $1):(t*file_interval/(24.*3600.)) w l lc palette lw 2 t "" unset format x -set origin 0, 0 -set size 0.33, 0.46 +set origin 0.043, 0.5 +set size 0.330, 0.46 set key top right samplen 1.0 font ",10" set xrange [0:7] -#set xtics format "" set yrange [0:240] -set ylabel "{/:Normal [kPa]}" offset -0.5 +set ylabel "{/:Normal [kPa]}" set xlabel "{/:Normal Time [d]}" set xtics 1 set ytics 50 - unset label -set label "D" font ",11" at graph 0.00,1.1 +set label "A" font ",11" at graph 0.00,1.1 plot \ "rate/timeseries.txt" u ($0/1008*7):($2/1000) w l lw 2 lc "black" t "Effective normal stress", \ "" u ($0/1008*7):($3/1000) w l lw 2 lt 3 dt "-" t "Water pressure" -set origin 0.33, 0.0 -set size 0.33, 0.46 +set origin 0.373, 0.5 +set size 0.313, 0.46 unset label -set label "E (Speed controlled)" font ",11" at graph 0.00,1.1 +set label "B (Speed controlled)" font ",11" at graph 0.00,1.1 unset y2label unset y2tics set ytics 1.0 @@ -103,25 +92,22 @@ set ylabel "{/:Normal Depth of max. strain rate [m]}" set arrow from 0,max_def_depth to 7,max_def_depth nohead front lc black dashtype 2 plot "rate/max_strainrate_depth.txt" u ($0/1008*7):(8.0 - $1) w l lw 2 t "" -set origin 0.67, 0.0 -set size 0.33, 0.46 +set origin 0.686, 0.5 +set size 0.313, 0.46 unset arrow unset label -set label "F (Stress controlled)" font ",11" at graph 0.00,1.1 +set label "C (Stress controlled)" font ",11" at graph 0.00,1.1 unset y2label unset y2tics -set yrange [0.0:10.0] +set yrange [0.0:9.8] set ytics 2.0 -#set ylabel "{/:Normal Depth of max. strain rate [m]}" -#plot "stress/max_strainrate_depth.txt" u ($0/1008*7):(8.0 - $1) w l lw 2 t "" set ylabel "{/:Normal Shear speed [km/d]}" plot "stress/timeseries.txt" u ($0/1008*7):($1*3600.*24./1000.) w l lw 2 lc "red" t "" unset label filename(t) = sprintf('rate/sim.output%05d.txt', t) -set origin 0.469, 0.62 +set origin 0.498, 0.11 set size 0.17, 0.28 -#clear set xlabel "{/:Normal Shear speed [km/a]}" font ",8" offset 0,0.7 set ylabel "Depth [m]" font ",8" offset 0.7,0.0 set ytics 1 font ",8" @@ -135,11 +121,9 @@ set yrange [4.0:0.0] plot for [t = t0:tend:inc*2] filename(t) u ($2*60*60*24*365/1000):(L_z - $1):(t*file_interval/(24.*3600.)) w l lc palette lw 1 t "" filename(t) = sprintf('stress/sim.output%05d.txt', t) -set origin 0.81, 0.62 +set origin 0.81, 0.11 set size 0.17, 0.28 -#clear set xlabel "{/:Normal Shear speed [km/d]}" font ",8" -#set ylabel "Depth [m]" font ",8" unset arrow xmin = -1 xmax = 10 @@ -147,13 +131,10 @@ set xrange [xmin:xmax] set xtics 2 plot for [t = t0:tend:inc*2] filename(t) u ($2*60*60*24/1000):(L_z - $1):(t*file_interval/(24.*3600.)) w l lc palette lw 1 t "" - filename(t) = sprintf('stress/sim.output%05d.txt', t) -set origin 0.02, 0.59 -set size 0.14, 0.22 -#clear +set origin 0.095, 0.09 +set size 0.14, 0.23 set xlabel "{/:Normal Water pressure [kPa]}" font ",8" offset 0,1.3 -#set ylabel "Depth [m]" font ",8" set ylabel "" unset arrow xmin = 0