hansen-zoet-exp

commit cb02243d0ed498165972a51aeeea12deb11280b8
parent 0cfe679ca76bb7c957d6667fa8b1e0cbc9121ca3
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Fri, 22 Apr 2022 12:57:03 +0200

add comparison to H&Z

Diffstat:
M
fig-effective_stress/Makefile
 | 
6
+++---
M
fig-effective_stress/fig-sediment_flux.gp
 | 
7
++++---
A
fig-effective_stress/rs_sediment_flux_100ma.txt
 | 
3
+++

3 files changed, 10 insertions(+), 6 deletions(-)
diff --git a/fig-effective_stress/Makefile b/fig-effective_stress/Makefile
@@ -14,12 +14,12 @@ strain_distribution_N50kPa.txt: $(BIN) Makefile
 		./$(BIN) \
 		-d 530.7e-6 \
 		-m 0.625 \
-		-c 30e3 \
+		-c 0e3 \
 		-s 3.169e-6 \
 		-L 0.11 \
 		-n $${P}e3 \
 		-A 0.48 \
-		-b 0.01 \
+		-b 0.001 \
 		> strain_distribution_N$${P}kPa.txt; \
 	done'
 
@@ -47,7 +47,7 @@ sediment_flux.txt: strain_distribution_N50kPa.txt
 	gnuplot fig-sediment_flux.gp > $@
 
 clean:
-	rm -f *.txt
+	rm -f strain_distribution_*.txt mohr_coulomb.txt sediment_flux.txt
 	rm -f ../$(FIG)-strain_distribution.pdf
 	rm -f ../$(FIG)-mohr_coulomb.pdf
 	rm -f ../$(FIG)-sediment_flux.pdf
diff --git a/fig-effective_stress/fig-sediment_flux.gp b/fig-effective_stress/fig-sediment_flux.gp
@@ -5,7 +5,7 @@ reset
 set terminal pdfcairo enhanced color size 7.5 cm, 5 cm
 set multiplot layout 1,1
 
-rswidth = 0.1
+rswidth = 1
 
 #set lmargin 7.0
 #set bmargin 3.5
@@ -22,4 +22,5 @@ set xlabel "Effective stress, [kPa]"
 #set ylabel "Spec. sediment flux [m²/a]"
 set ylabel "Sediment flux [m³/a]"
 set key bottom right font ",10" #samplen 0.9
-plot "sediment_flux.txt" u ($1/1000):($2*3600*24*365.25*rswidth) w p lw 2 notitle-
\ No newline at end of file
+plot "sediment_flux.txt" u ($1/1000):($2*3600*24*365.25*rswidth) w lp lw 2 t "CNGF-PF", \
+	"rs_sediment_flux_100ma.txt" u ($1/1000):($2) w p t "Hansen and Zoet 2022"+
\ No newline at end of file
diff --git a/fig-effective_stress/rs_sediment_flux_100ma.txt b/fig-effective_stress/rs_sediment_flux_100ma.txt
@@ -0,0 +1,3 @@
+8.53e3	0.351
+22.2e3	0.510
+50.4e3	0.661