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 4c64982dd66f0baf2645d7cf248551df233616f7
parent 9ab7accfac23f913f9fad33feda8046399f0a223
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Tue, 28 Jan 2020 16:40:44 +0100

Clean up fig 2 code and change color scheme

Diffstat:
Mfig-stick_slip_rate_depth/fig.gp | 21++-------------------
Mfig-stick_slip_stress_depth/fig.gp | 36++----------------------------------
2 files changed, 4 insertions(+), 53 deletions(-)

diff --git a/fig-stick_slip_rate_depth/fig.gp b/fig-stick_slip_rate_depth/fig.gp @@ -23,30 +23,16 @@ tend = 71*inc filename(t) = sprintf('sim.output%05d.txt', t) set pm3d map -set palette rgb 10,6,6 +#set palette rgb 10,6,6 +set palette -#load '../skindepth.gp' -#k = 2e-17 -#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" unset colorbox -#set xlabel "{/:Normal Water pressure [kPa]}" 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):(L_z - $1):(t-t0) w l lc palette lw 1 t "" - -unset label -#set label "{/:Bold c}" at graph 0.90,0.93 set label "A" at screen 0.090,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):(L_z - $1):(t-t0) w l lc palette lw 1 t "" unset label @@ -55,8 +41,6 @@ set label "B (rate controlled)" at screen 0.395,0.97 font ",12" set xlabel "{/:Normal Shear speed [km/a]}" set ylabel "" set xrange [0.0:1.0] -#unset xrange -#set xtics (0,0.02,0.04,0.06) unset ytics unset arrow set arrow from 0,skin_depth to 1,skin_depth nohead lc black front @@ -64,7 +48,6 @@ set arrow from 0,max_def_depth to 1,max_def_depth nohead front lc black dashtype 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.705,0.97 font ",12" set xlabel "{/:Normal Shear strain rate [1/d]}" set ylabel "" diff --git a/fig-stick_slip_stress_depth/fig.gp b/fig-stick_slip_stress_depth/fig.gp @@ -1,7 +1,5 @@ #!/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 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 \ @@ -21,50 +19,20 @@ tend = 71*inc filename(t) = sprintf('sim.output%05d.txt', t) set pm3d map -set palette rgb 10,6,6 +#set palette rgb 10,6,6 +set palette -#load '../skindepth.gp' -#k = 2e-17 -#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" unset colorbox -#set xlabel "{/:Normal Water pressure [kPa]}" -##set ylabel "{/:Normal Vertical position [m]}" -unset ylabel -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):(L_z - $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):(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.035,0.97 font ",12" set xlabel "{/:Normal Shear speed [km/d]}" set ylabel "" set xrange [0.0:10.0] -#unset xrange -#set xtics (0,0.02,0.04,0.06) 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):(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.515,0.97 font ",12" set xlabel "{/:Normal Shear strain rate [1/m]}" set ylabel ""