hansen-zoet-exp

commit bc8eb47d791d184c1ef5aa1a4870f8759b57a8b0
parent 07937f709c4da6753ad3ce6853a82b58cb192928
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Thu,  2 Jun 2022 14:52:30 +0200

add plot with 100ma velocity only

Diffstat:
M
Makefile
 | 
13
+++++++++++--
A
fig-sediment_flux-100ma.gp
 | 
30
++++++++++++++++++++++++++++++

2 files changed, 41 insertions(+), 2 deletions(-)
diff --git a/Makefile b/Makefile
@@ -5,9 +5,10 @@ BIN = ${REPO}/cngf-pf
 PDFS = \
 	fig-strain_distribution.pdf \
 	fig-mohr_coulomb.pdf \
-	fig-sediment_flux.pdf
+	fig-sediment_flux.pdf \
+	fig-sediment_flux-100ma.pdf \
 
-stressvals_ != range -l -N 20 0 2
+stressvals_ != range -l -N 40 0.01 100
 stressvals_ += 8 22 50
 stressvals != echo ${stressvals_} | sort -n
 velvals = 3.169e-6 9.506e-6
@@ -35,6 +36,11 @@ strain_distribution_N50kPa_v3.169e-6ms.txt: ${BIN} Makefile
 			-b 0.002 \
 			> strain_distribution_N$${P}kPa_v$${v}ms.txt; \
 	done; done'  # also close fit with lower A and higher b
+	# A doubles -> scale factor doubles, stress exponent is constant at 1/3
+	# q \propto cuberoot(N)
+	# A = 0.66 -> a = 2.62
+	# A = 0.33 -> a = 1.79
+	# higher values of L do not make a difference in flux
 
 mohr_coulomb_v3.169e-6ms.txt: strain_distribution_N50kPa_v3.169e-6ms.txt
 	sh -c '\
@@ -62,6 +68,9 @@ fig-mohr_coulomb.pdf: fig-mohr_coulomb.gp mohr_coulomb_v3.169e-6ms.txt
 fig-sediment_flux.pdf: fig-sediment_flux.gp sediment_flux_v3.169e-6ms.txt
 	gnuplot fig-sediment_flux.gp > $@
 
+fig-sediment_flux-100ma.pdf: fig-sediment_flux-100ma.gp sediment_flux_v3.169e-6ms.txt
+	gnuplot fig-sediment_flux-100ma.gp > $@
+
 clean:
 	rm -f strain_distribution_*.txt mohr_coulomb_*.txt sediment_flux_*.txt
 	rm -f ${PDFS}
diff --git a/fig-sediment_flux-100ma.gp b/fig-sediment_flux-100ma.gp
@@ -0,0 +1,30 @@
+#!/usr/bin/env gnuplot
+
+reset
+
+set terminal pdfcairo enhanced color size 7.5 cm, 5 cm
+set multiplot layout 1,1
+
+rswidth = 1
+
+#set lmargin 7.0
+#set bmargin 3.5
+#set rmargin 2.0
+#set tmargin 2.0
+
+# https://sodocumentation.net/gnuplot/topic/8825/fit-data-with-gnuplot
+f(x, v) = a * x**b * v
+fit f(x, 100) "sediment_flux_v3.169e-6ms.txt" u ($1/1000):($2*3600*24*365.25*rswidth) via a,b
+fitparams = sprintf("{/:Italic q}_t = %.2e {/:Italic v N'}^{%.3g}", a, b)
+
+set xlabel "Effective stress, {/:Italic N'} [kPa]"
+set ylabel "Sediment flux, {/:Italic q}_t [m³/a]"
+set key bottom right font ",09" invert samplen 1.0
+#set title "{/:Italic v} = 100 m/a"
+set xrange [0:210]
+set xtics 0,50,200
+set yrange [0:1.1]
+plot f(x, 100) lt 1 t fitparams, \
+	"sediment_flux_v3.169e-6ms.txt" u ($1/1000):($2*3600*24*365.25*rswidth) w p lt 2 ps 0.5 t "CNGF-PF (Damsgaard et al., 2020)", \
+	"rs_sediment_flux_100ma_high-N.txt" u ($1/1000):($2) w p lt 4 ps 1 t "Hansen and Zoet (2022)", \
+	"rs_sediment_flux_100ma.txt" u ($1/1000):($2) w p lt 4 ps 1 t ""