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 e35a1f2be8067ebdf96e269f1d5351235103a6e9
parent 6864175afb509e334ec43037eaf205ef0323c572
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Thu, 30 Jan 2020 13:43:16 +0100

Update fig3 colors

Diffstat:
Mfig-hysteresis_rate/fig.gp | 4++--
Mfig-hysteresis_stress/fig.gp | 4++--
Mfig-stick_slip_stress_depth/fig.gp | 51++++++++++++++++++++++++++++++++++-----------------
3 files changed, 38 insertions(+), 21 deletions(-)

diff --git a/fig-hysteresis_rate/fig.gp b/fig-hysteresis_rate/fig.gp @@ -22,7 +22,7 @@ unset label #set label "{/:Bold d}" at graph 0.02,0.93 set label "A" at graph 0.02,0.93 font ",12" plot \ -"timeseries_combined.txt" u ($2/1000):($4*$2/1000) w l lt 1 t "" +"timeseries_combined.txt" u ($2/1000):($4*$2/1000) w l lt 1 lc rgb "#40e0d0" t "" unset x2tics unset ytics @@ -38,7 +38,7 @@ set arrow from 170,1.4 to 170.3,0.6 head unset label #set label "{/:Bold e}" at graph 0.02,0.93 font ",12" set label "B" at graph 0.02,0.93 font ",12" -plot "" u ($2/1000):(8.0 - $6) w l lt 4 t "" +plot "" u ($2/1000):(8.0 - $6) w l lt 4 lc rgb "#87cefa" t "" set xtics format "%g" set xlabel "{/:Normal Eff. normal stress [kPa]}" diff --git a/fig-hysteresis_stress/fig.gp b/fig-hysteresis_stress/fig.gp @@ -24,7 +24,7 @@ unset label #set label "{/:Bold a}" at graph 0.02,0.93 set label "D" at graph 0.02,0.93 font ",12" plot \ -"timeseries_combined.txt" u ($2/1000):($1*60*60*24/1000) w l lt 7 t "" +"timeseries_combined.txt" u ($2/1000):($1*60*60*24/1000) w l lt 7 lc rgb "#008000" t "" unset x2tics unset ytics @@ -44,7 +44,7 @@ set arrow from 165,1.2 to 165,0.7 head unset label #set label "{/:Bold b}" at graph 0.02,0.93 set label "E" at graph 0.02,0.93 font ",12" -plot "" u ($2/1000):(8.0 - $6) w l lt 4 t "" +plot "" u ($2/1000):(8.0 - $6) w l lt 4 lc rgb "#87cefa" t "" set xtics format "%g" set xlabel "{/:Normal Eff. normal stress [kPa]}" diff --git a/fig-stick_slip_stress_depth/fig.gp b/fig-stick_slip_stress_depth/fig.gp @@ -1,15 +1,19 @@ #!/usr/bin/env gnuplot -set terminal pdfcairo color size 7.12 cm, 4.45 cm font ",10" -set multiplot layout 1,2 \ - margins 0.04,0.95,0.23,0.90 \ - spacing 0.05,0.05 +#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 terminal pdfcairo color size 10.68 cm, 4.45 cm font ",10" +set multiplot layout 1,3 \ + margins 0.09,0.98,0.23,0.90 \ + spacing 0.04,0.04 + 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] @@ -23,29 +27,42 @@ set pm3d map set palette unset colorbox -set label "D (stress controlled)" at screen 0.035,0.97 font ",12" -set xlabel "{/:Normal Shear speed [km/d]}" -set ylabel "" -xmin = -0.5 -xmax = 10.0 +set ylabel "{/:Normal Depth [m]}" offset 1.0 +xmin = 0.0 +xmax = 200.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 -unset ytics -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 "" +set label "A" at screen 0.090,0.97 font ",12" +set xlabel "{/:Normal Eff. normal stress [kPa]}" +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 "E (stress controlled)" at screen 0.515,0.97 font ",12" -set xlabel "{/:Normal Shear strain rate [1/m]}" +set label "B (rate controlled)" at screen 0.395,0.97 font ",12" +set xlabel "{/:Normal Shear speed [km/a]}" set ylabel "" -xmin = -0.8 -xmax = 12.0 -set xrange [xmin:xmax] unset ytics unset arrow +xmin = -0.1 +xmax = 1.1 +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] 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 "C (rate controlled)" at screen 0.705,0.97 font ",12" +set xlabel "{/:Normal Shear strain rate [1/d]}" +set ylabel "" +xmin = -2.0 +xmax = 22.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] filename(t) u ($6*60):(L_z - $1):(t-t0) w l lc palette lw 1 t "" +set xrange [xmin:xmax] +unset ytics +unset arrow +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