tPlot speed in panel F and add insert plot of water pressures to panel A - cngf-pf-exp1 - experiments for first paper with continuum granular model
(HTM) git clone git://src.adamsgaard.dk/manus_continuum_granular1_exp
(DIR) Log
(DIR) Files
(DIR) Refs
(DIR) Submodules
(DIR) README
(DIR) LICENSE
---
(DIR) commit 038accc83d5b5f06ba28df5339fbf4b39f9ec779
(DIR) parent a1b10832f92cbf31f9fa1afa7a9cfd3e0e04a340
(HTM) Author: Anders Damsgaard <anders@adamsgaard.dk>
Date: Fri, 7 Aug 2020 15:00:21 +0200
Plot speed in panel F and add insert plot of water pressures to panel A
Diffstat:
M fig-stick_slip_depth/fig.gp | 32 +++++++++++++++++++++++++------
M flux_estimate/Makefile | 3 +++
2 files changed, 29 insertions(+), 6 deletions(-)
---
(DIR) diff --git a/fig-stick_slip_depth/fig.gp b/fig-stick_slip_depth/fig.gp
t@@ -23,12 +23,12 @@ set palette
set ylabel "{/:Normal Depth [m]}" offset 0.5
-set origin 0.0, 0.50
-set size 0.33, 0.46
+set origin -0.02, 0.50
+set size 0.36, 0.46
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 90,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 cblabel "Time [d]"
t@@ -110,9 +110,12 @@ unset label
set label "F (Stress controlled)" font ",11" at graph 0.00,1.1
unset y2label
unset y2tics
-set yrange [4.0:0.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 yrange [0.0:10.0]
+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)
t@@ -128,6 +131,7 @@ xmin = -0.1
xmax = 1.1
set arrow from xmin,max_def_depth to xmax,max_def_depth nohead front lc black dashtype 2
set xrange [xmin:xmax]
+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)
t@@ -144,4 +148,20 @@ 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 xlabel "{/:Normal Water pressure [kPa]}" font ",8" offset 0,1.3
+#set ylabel "Depth [m]" font ",8"
+set ylabel ""
+unset arrow
+xmin = 0
+xmax = 200.0
+set arrow from xmin,skin_depth to xmax,skin_depth nohead lc black front
+set xrange [xmin:xmax]
+set xtics 100 offset 0,0.5
+set ytics 2 offset 0.5
+plot for [t = t0:tend:inc*2] filename(t) u ($4/1e3):(L_z - $1):(t*file_interval/(24.*3600.)) w l lc palette lw 1 t ""
+
unset multiplot
(DIR) diff --git a/flux_estimate/Makefile b/flux_estimate/Makefile
t@@ -53,6 +53,9 @@ variable_flux: variable.output00095.txt
printf 'Yearly volume: %.5g m^3/a\n' "$$yearly_volume" && \
printf 'Duration: %.5g a\n' "$$duration"
+
+# -a 50e3 is now 79 m^2/a
+# -a 99e3 is now 102 m^2/a
variable_hi.output00095.txt:
$(BIN) $(FLAGS) \
-F -a 50e3 \