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 d82209f391f610a27ca4367d826f76584ed56110
parent 6394add6eed5ad624b7699e13663b5fb125f82de
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Fri, 24 Jan 2020 12:53:37 +0100

Use depth as axis

Diffstat:
Mfig-experimental-setup.pdf | 0
Mfig-experimental-setup.svg | 6------
Mfig-stick_slip_rate_depth/Makefile | 2+-
Mfig-stick_slip_rate_depth/fig.gp | 18++++++++++--------
Mfig-stick_slip_stress_depth/Makefile | 2+-
Mfig-stick_slip_stress_depth/fig.gp | 18+++++++++---------
6 files changed, 21 insertions(+), 25 deletions(-)

diff --git a/fig-experimental-setup.pdf b/fig-experimental-setup.pdf Binary files differ. diff --git a/fig-experimental-setup.svg b/fig-experimental-setup.svg @@ -496,12 +496,6 @@ x="86.624481" y="79.37886" style="font-style:italic;font-variant:normal;font-weight:normal;font-stretch:normal;font-size:3.52777767px;font-family:'DejaVu Serif';-inkscape-font-specification:'DejaVu Serif Italic';stroke-width:0.26458332">x</tspan></text> - <path - style="fill:none;stroke:#000000;stroke-width:0.25570717;stroke-linecap:butt;stroke-linejoin:miter;stroke-miterlimit:4;stroke-dasharray:0.51141437, 0.25570718;stroke-dashoffset:2.64999986;stroke-opacity:1" - d="m 56.570608,52.493069 25.469573,0.175596" - id="path1593" - inkscape:connector-curvature="0" - sodipodi:nodetypes="cc" /> <text xml:space="preserve" style="font-style:normal;font-weight:normal;font-size:10.58333302px;line-height:1.25;font-family:sans-serif;letter-spacing:0px;word-spacing:0px;fill:#000000;fill-opacity:1;stroke:none;stroke-width:0.26458332" diff --git a/fig-stick_slip_rate_depth/Makefile b/fig-stick_slip_rate_depth/Makefile @@ -17,7 +17,7 @@ sim.output00001.txt: ../fig-stick_slip_rate/sim.output00001.txt ); \ max_def_depth=$$(echo $$out | cut -d" " -f1); \ skin_depth=$$(echo $$out | cut -d" " -f2); \ - gnuplot -e "skin_depth=8.0-$$skin_depth; max_def_depth=8.0-$$max_def_depth" \ + gnuplot -e "skin_depth=$$skin_depth; max_def_depth=$$max_def_depth" \ fig.gp > $@' clean: diff --git a/fig-stick_slip_rate_depth/fig.gp b/fig-stick_slip_rate_depth/fig.gp @@ -4,16 +4,18 @@ #set terminal pdfcairo color size 13.35 cm, 4.45 cm font ",10" set terminal pdfcairo color size 10.68 cm, 4.45 cm font ",10" set multiplot layout 1,3 \ - margins 0.09,0.97,0.23,0.90 \ + margins 0.09,0.98,0.23,0.90 \ spacing 0.04,0.04 -set yrange [4.0:8.0] set key bottom right #samplen 0.9 file_interval=60*10 plot_interval=60*60 inc=plot_interval/file_interval +#set yrange [4.0:8.0] +L_z=8.0 +set yrange [4.0:0.0] t0 = 47*inc tend = 71*inc @@ -32,11 +34,11 @@ set palette rgb 10,6,6 #set label "B" at screen 0.09,0.97 font ",12" unset colorbox #set xlabel "{/:Normal Water pressure [kPa]}" -set ylabel "{/:Normal Vertical position [m]}" offset 1.0 +set ylabel "{/:Normal Depth [m]}" offset 1.0 #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:inc] 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):(L_z - $1):(t-t0) w l lc palette lw 1 t "" unset label #set label "{/:Bold c}" at graph 0.90,0.93 @@ -45,7 +47,7 @@ set xlabel "{/:Normal Eff. normal stress [kPa]}" #set ylabel "" set xrange [0:200] #unset ytics -plot for [t = t0:tend:inc] 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):(L_z - $1):(t-t0) w l lc palette lw 1 t "" unset label #set label "{/:Bold d}" at graph 0.90,0.93 @@ -59,11 +61,11 @@ 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:inc] 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):(L_z - $1):(t-t0) w l lc palette lw 1 t "" unset label #set label "{/:Bold e}" at graph 0.90,0.93 -set label "C (rate controlled)" at screen 0.700,0.97 font ",12" +set label "C (rate controlled)" at screen 0.705,0.97 font ",12" set xlabel "{/:Normal Shear strain rate [1/d]}" set ylabel "" set xrange [0:20] @@ -71,7 +73,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:inc] 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):(L_z - $1):(t-t0) w l lc palette lw 1 t "" unset format x unset multiplot diff --git a/fig-stick_slip_stress_depth/Makefile b/fig-stick_slip_stress_depth/Makefile @@ -17,7 +17,7 @@ sim.output00001.txt: ../fig-stick_slip_stress/sim.output00001.txt ); \ max_def_depth=$$(echo $$out | cut -d" " -f1); \ skin_depth=$$(echo $$out | cut -d" " -f2); \ - gnuplot -e "skin_depth=8.0-$$skin_depth; max_def_depth=8.0-$$max_def_depth" \ + gnuplot -e "skin_depth=$$skin_depth; max_def_depth=$$max_def_depth" \ fig.gp > $@' clean: diff --git a/fig-stick_slip_stress_depth/fig.gp b/fig-stick_slip_stress_depth/fig.gp @@ -4,16 +4,16 @@ #set terminal pdfcairo color size 13.35 cm, 4.45 cm font ",10" set terminal pdfcairo color size 7.12 cm, 4.45 cm font ",10" set multiplot layout 1,2 \ - margins 0.05,0.95,0.23,0.90 \ + margins 0.04,0.95,0.23,0.90 \ spacing 0.05,0.05 -set yrange [4.0:8.0] - set key bottom right #samplen 0.9 file_interval=60*10 plot_interval=60*60 inc=plot_interval/file_interval +L_z=8.0 +set yrange [4.0:0.0] t0 = 47*inc tend = 71*inc @@ -38,7 +38,7 @@ unset ytics #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:inc] 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):(L_z - $1):(t-t0) w l lc palette lw 1 t "" #unset label ##set label "{/:Bold c}" at graph 0.90,0.93 @@ -47,11 +47,11 @@ unset ytics #set ylabel "" #set xrange [0:200] #unset ytics -#plot for [t = t0:tend:inc] 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):(L_z - $1):(t-t0) w l lc palette lw 1 t "" unset label #set label "{/:Bold d}" at graph 0.90,0.93 -set label "D (stress controlled)" at screen 0.045,0.97 font ",12" +set label "D (stress controlled)" at screen 0.035,0.97 font ",12" set xlabel "{/:Normal Shear velocity [km/d]}" set ylabel "" set xrange [0.0:10.0] @@ -61,11 +61,11 @@ unset ytics unset arrow set arrow from 0,skin_depth to 10,skin_depth nohead lc black front set arrow from 0,max_def_depth to 10,max_def_depth nohead front lc black dashtype 2 -plot for [t = t0:tend:inc] filename(t) u ($2*60*60*24/1000):1:(t-t0) w l lc palette lw 1 t "" +plot for [t = t0:tend:inc] filename(t) u ($2*60*60*24/1000):(L_z - $1):(t-t0) w l lc palette lw 1 t "" unset label #set label "{/:Bold e}" at graph 0.90,0.93 -set label "E (stress controlled)" at screen 0.520,0.97 font ",12" +set label "E (stress controlled)" at screen 0.515,0.97 font ",12" set xlabel "{/:Normal Shear strain rate [1/m]}" set ylabel "" set xrange [0:12] @@ -73,7 +73,7 @@ unset ytics unset arrow set arrow from 0,skin_depth to 12,skin_depth nohead lc black front set arrow from 0,max_def_depth to 12,max_def_depth nohead front lc black dashtype 2 -plot for [t = t0:tend:inc] filename(t) u ($6*60):1:(t-t0) w l lc palette lw 1 t "" +plot for [t = t0:tend:inc] filename(t) u ($6*60):(L_z - $1):(t-t0) w l lc palette lw 1 t "" unset format x unset multiplot