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 6394add6eed5ad624b7699e13663b5fb125f82de
parent 66f8a43d8889cce129d3dc08193f7922d16f3858
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Wed, 22 Jan 2020 00:07:21 +0100

Update depth figures

Diffstat:
Mfig-stick_slip_rate_depth/fig.gp | 29+++++++++++++++--------------
Mfig-stick_slip_stress_depth/fig.gp | 47+++++++++++++++++++++++++----------------------
2 files changed, 40 insertions(+), 36 deletions(-)

diff --git a/fig-stick_slip_rate_depth/fig.gp b/fig-stick_slip_rate_depth/fig.gp @@ -1,10 +1,11 @@ #!/usr/bin/env gnuplot #set terminal pdfcairo color size 7.5 cm, 7.5 cm -set terminal pdfcairo color size 13.35 cm, 4.45 cm font ",10" -set multiplot layout 1,4 \ +#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 \ - spacing 0.03,0.03 + spacing 0.04,0.04 set yrange [4.0:8.0] @@ -27,28 +28,28 @@ set palette rgb 10,6,6 #f = 1.0/(3600.0*24.0) #f(x) = 8.0 - skindepth(f,k) -#set label "{/:Bold b}" at graph 0.90,0.93 -set label "B" at screen 0.09,0.97 font ",12" +##set label "{/:Bold b}" at graph 0.90,0.93 +#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]}" -set xrange [0:200] +#set xlabel "{/:Normal Water pressure [kPa]}" +set ylabel "{/:Normal Vertical position [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):1:(t-t0) w l lc palette lw 1 t "" unset label #set label "{/:Bold c}" at graph 0.90,0.93 -set label "C" at screen 0.315,0.97 font ",12" +set label "A" at screen 0.090,0.97 font ",12" set xlabel "{/:Normal Eff. normal stress [kPa]}" -set ylabel "" +#set ylabel "" set xrange [0:200] -unset ytics +#unset ytics 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 d}" at graph 0.90,0.93 -set label "D" at screen 0.545,0.97 font ",12" +set label "B (rate controlled)" at screen 0.395,0.97 font ",12" set xlabel "{/:Normal Shear velocity [km/a]}" set ylabel "" set xrange [0.0:1.0] @@ -62,7 +63,7 @@ plot for [t = t0:tend:inc] filename(t) u ($2*60*60*24*365/1000):1:(t-t0) w l lc unset label #set label "{/:Bold e}" at graph 0.90,0.93 -set label "E" at screen 0.772,0.97 font ",12" +set label "C (rate controlled)" at screen 0.700,0.97 font ",12" set xlabel "{/:Normal Shear strain rate [1/d]}" set ylabel "" set xrange [0:20] diff --git a/fig-stick_slip_stress_depth/fig.gp b/fig-stick_slip_stress_depth/fig.gp @@ -1,10 +1,11 @@ #!/usr/bin/env gnuplot #set terminal pdfcairo color size 7.5 cm, 7.5 cm -set terminal pdfcairo color size 13.35 cm, 4.45 cm font ",10" -set multiplot layout 1,4 \ - margins 0.09,0.97,0.23,0.90 \ - spacing 0.03,0.03 +#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 \ + spacing 0.05,0.05 set yrange [4.0:8.0] @@ -27,28 +28,30 @@ set palette rgb 10,6,6 #f = 1.0/(3600.0*24.0) #f(x) = 8.0 - skindepth(f,k) -#set label "{/:Bold b}" at graph 0.90,0.93 -set label "B" at screen 0.09,0.97 font ",12" +##set label "{/:Bold b}" at graph 0.90,0.93 +#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]}" -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 "" - -unset label -#set label "{/:Bold c}" at graph 0.90,0.93 -set label "C" at screen 0.315,0.97 font ",12" -set xlabel "{/:Normal Eff. normal stress [kPa]}" -set ylabel "" -set xrange [0:200] +#set xlabel "{/:Normal Water pressure [kPa]}" +##set ylabel "{/:Normal Vertical position [m]}" +unset ylabel unset ytics -plot for [t = t0:tend:inc] filename(t) u ($3/1e3):1:(t-t0) w l lc palette lw 1 t "" +#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 "" + +#unset label +##set label "{/:Bold c}" at graph 0.90,0.93 +#set label "C" at screen 0.315,0.97 font ",12" +#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 "" unset label #set label "{/:Bold d}" at graph 0.90,0.93 -set label "D" at screen 0.545,0.97 font ",12" +set label "D (stress controlled)" at screen 0.045,0.97 font ",12" set xlabel "{/:Normal Shear velocity [km/d]}" set ylabel "" set xrange [0.0:10.0] @@ -62,7 +65,7 @@ plot for [t = t0:tend:inc] filename(t) u ($2*60*60*24/1000):1:(t-t0) w l lc pale unset label #set label "{/:Bold e}" at graph 0.90,0.93 -set label "E" at screen 0.772,0.97 font ",12" +set label "E (stress controlled)" at screen 0.520,0.97 font ",12" set xlabel "{/:Normal Shear strain rate [1/m]}" set ylabel "" set xrange [0:12]