commit 287ee94897c0677e8f148e89682fdb61b76cbfbd
parent 861dda3a017d02cb4a3212f654d8abaac3b8b813
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date: Fri, 22 Apr 2022 14:21:15 +0200
add 300 m/a results and comparison
Diffstat:
6 files changed, 70 insertions(+), 45 deletions(-)
diff --git a/fig-effective_stress/Makefile b/fig-effective_stress/Makefile
@@ -1,54 +1,59 @@
BIN = ../cngf-pf/cngf-pf
FIG = $(shell basename $(CURDIR))
-stressvals = 5 8.53 10 15 20 22.2 25 30 35 40 45 50 50.4
+stressvals = 5 8 10 15 20 22 25 30 35 40 45 50
+velvals = 3.169e-6 9.506e-6
+
default: \
../$(FIG)-strain_distribution.pdf \
../$(FIG)-mohr_coulomb.pdf \
../$(FIG)-sediment_flux.pdf
-strain_distribution_N50kPa.txt: $(BIN) Makefile
+strain_distribution_N50kPa_v3.169e-6ms.txt: $(BIN) Makefile
sh -c '\
- for P in ${stressvals}; do \
- ./$(BIN) \
- -d 530.7e-6 \
- -m 0.625 \
- -c 0e3 \
- -s 3.169e-6 \
- -L 0.11 \
- -n $${P}e3 \
- -A 0.48 \
- -b 0.0003 \
- > strain_distribution_N$${P}kPa.txt; \
- done'
+ for v in ${velvals}; do \
+ for P in ${stressvals}; do \
+ ./$(BIN) \
+ -d 530.7e-6 \
+ -m 0.625 \
+ -c 0e3 \
+ -s $${v} \
+ -L 0.11 \
+ -n $${P}e3 \
+ -A 0.48 \
+ -b 0.0003 \
+ > strain_distribution_N$${P}kPa_v$${v}ms.txt; \
+ done; done'
# also close: lower A, higher b
-mohr_coulomb.txt: strain_distribution_N50kPa.txt
+mohr_coulomb_v3.169e-6ms.txt: strain_distribution_N50kPa_v3.169e-6ms.txt
sh -c '\
- rm -f $@; \
- for P in ${stressvals}; do \
- tail -n 1 strain_distribution_N$${P}kPa.txt >> $@; \
- done'
+ rm -f mohr_coulomb_v*ms.txt; \
+ for v in ${velvals}; do \
+ for P in ${stressvals}; do \
+ tail -n 1 strain_distribution_N$${P}kPa_v$${v}ms.txt >> mohr_coulomb_v$${v}ms.txt; \
+ done; done'
-sediment_flux.txt: strain_distribution_N50kPa.txt
+sediment_flux_v3.169e-6ms.txt: strain_distribution_N50kPa_v3.169e-6ms.txt
sh -c '\
- rm -f $@; \
- for P in ${stressvals}; do \
- (printf "%ge3\t" $$P; shear_flux < strain_distribution_N$${P}kPa.txt) >> $@; \
- done'
+ rm -f sediment_flux_v*ms.txt; \
+ for v in ${velvals}; do \
+ for P in ${stressvals}; do \
+ (printf "%ge3\t" $$P; shear_flux < strain_distribution_N$${P}kPa_v$${v}ms.txt) >> sediment_flux_v$${v}ms.txt; \
+ done; done'
-../$(FIG)-strain_distribution.pdf: fig-strain_distribution.gp strain_distribution_N50kPa.txt
+../$(FIG)-strain_distribution.pdf: fig-strain_distribution.gp strain_distribution_N50kPa_v3.169e-6ms.txt
gnuplot fig-strain_distribution.gp > $@
-../$(FIG)-mohr_coulomb.pdf: fig-mohr_coulomb.gp mohr_coulomb.txt
+../$(FIG)-mohr_coulomb.pdf: fig-mohr_coulomb.gp mohr_coulomb_v3.169e-6ms.txt
gnuplot fig-mohr_coulomb.gp > $@
-../$(FIG)-sediment_flux.pdf: fig-sediment_flux.gp sediment_flux.txt
+../$(FIG)-sediment_flux.pdf: fig-sediment_flux.gp sediment_flux_v3.169e-6ms.txt
gnuplot fig-sediment_flux.gp > $@
clean:
- rm -f strain_distribution_*.txt mohr_coulomb.txt sediment_flux.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-mohr_coulomb.gp b/fig-effective_stress/fig-mohr_coulomb.gp
@@ -16,8 +16,12 @@ set multiplot layout 2,1
#set xlabel "Effective normal stress, {/Symbol s}_n' [kPa]"
set ylabel "Shear stress, {/Symbol t} [kPa]"
-plot "mohr_coulomb.txt" u ($3/1e3):($3*$5/1e3) w lp title ""
+set key bottom right font ",10"
+plot "mohr_coulomb_v3.169e-6ms.txt" u ($3/1e3):($3*$5/1e3) w lp title "v = 100 m/a", \
+ "mohr_coulomb_v9.506e-6ms.txt" u ($3/1e3):($3*$5/1e3) w lp title "v = 300 m/a"
set xlabel "Effective normal stress, {/Symbol s}_n' [kPa]"
set ylabel "Shear friction, {/Symbol m} [-]"
-plot "mohr_coulomb.txt" u ($3/1e3):5 w lp title ""
+set key top right font ",10"
+plot "mohr_coulomb_v3.169e-6ms.txt" u ($3/1e3):5 w lp title "v = 100 m/a", \
+ "mohr_coulomb_v9.506e-6ms.txt" u ($3/1e3):5 w lp title "v = 300 m/a"
diff --git a/fig-effective_stress/fig-sediment_flux.gp b/fig-effective_stress/fig-sediment_flux.gp
@@ -24,9 +24,9 @@ set xlabel "Effective stress [kPa]"
set ylabel "Sediment flux [m³/a]"
set key bottom right font ",10" #samplen 0.9
set title "v = 100 m/a"
-plot "sediment_flux.txt" u ($1/1000):($2*3600*24*365.25*rswidth) w l lw 2 t "CNGF-PF", \
+plot "sediment_flux_v3.169e-6ms.txt" u ($1/1000):($2*3600*24*365.25*rswidth) w l lw 2 t "CNGF-PF", \
"rs_sediment_flux_100ma.txt" u ($1/1000):($2) w errorbars t "Hansen and Zoet 2022"
set title "v = 300 m/a"
-plot "sediment_flux.txt" u ($1/1000):($2*3600*24*365.25*rswidth) w l lw 2 t "CNGF-PF", \
- "rs_sediment_flux_100ma.txt" u ($1/1000):($2) w errorbars t "Hansen and Zoet 2022"-
\ No newline at end of file
+plot "sediment_flux_v9.506e-6ms.txt" u ($1/1000):($2*3600*24*365.25*rswidth) w l lw 2 t "CNGF-PF", \
+ "rs_sediment_flux_300ma.txt" u ($1/1000):($2) w errorbars t "Hansen and Zoet 2022"
diff --git a/fig-effective_stress/fig-strain_distribution.gp b/fig-effective_stress/fig-strain_distribution.gp
@@ -17,16 +17,30 @@ set yrange [0 : 0.11]
set xlabel "Horizontal velocity, v_x [m/a]"
set ylabel "Vertical position, z [m]"
set key bottom right font ",10" #samplen 0.9
-plot "strain_distribution_N8kPa.txt" u ($2*365*24*3600):1 w l lw 2 title "{/Symbol s}_n' = 8 kPa", \
- "strain_distribution_N22kPa.txt" u ($2*365*24*3600):1 w l lw 2 title "{/Symbol s}_n' = 22 kPa", \
- "strain_distribution_N50kPa.txt" u ($2*365*24*3600):1 w l lw 2 title "{/Symbol s}_n' = 50 kPa"
+set title "v = 100 m/a"
+plot "strain_distribution_N8kPa_v3.169e-6ms.txt" u ($2*365*24*3600):1 w l lw 2 title "{/Symbol s}_n' = 8 kPa", \
+ "strain_distribution_N22kPa_v3.169e-6ms.txt" u ($2*365*24*3600):1 w l lw 2 title "{/Symbol s}_n' = 22 kPa", \
+ "strain_distribution_N50kPa_v3.169e-6ms.txt" u ($2*365*24*3600):1 w l lw 2 title "{/Symbol s}_n' = 50 kPa"
-unset xrange
+#unset xrange
+
+##set xlabel "Normalized horizontal velocity, v_x [-]"
+#set xlabel "Shear-strain rate, dv_x/dz [1/d]"
+#set ylabel "Vertical position, z [m]"
+#set key bottom right font ",10" #samplen 0.9
+#plot "strain_distribution_N8kPa_v9.506e-6ms.txt" u ($6*24*3600):1 w l lw 2 title "{/Symbol s}_n' = 8 kPa", \
+ #"strain_distribution_N22kPa_v9.506e-6ms.txt" u ($6*24*3600):1 w l lw 2 title "{/Symbol s}_n' = 22 kPa", \
+ #"strain_distribution_N50kPa_v9.506e-6ms.txt" u ($6*24*3600):1 w l lw 2 title "{/Symbol s}_n' = 50 kPa"
+
+
+set xrange [-5 : 305]
+set yrange [0 : 0.11]
#set xlabel "Normalized horizontal velocity, v_x [-]"
-set xlabel "Shear-strain rate, dv_x/dz [1/d]"
+set xlabel "Horizontal velocity, v_x [m/a]"
set ylabel "Vertical position, z [m]"
set key bottom right font ",10" #samplen 0.9
-plot "strain_distribution_N8kPa.txt" u ($6*24*3600):1 w l lw 2 title "{/Symbol s}_n' = 8 kPa", \
- "strain_distribution_N22kPa.txt" u ($6*24*3600):1 w l lw 2 title "{/Symbol s}_n' = 22 kPa", \
- "strain_distribution_N50kPa.txt" u ($6*24*3600):1 w l lw 2 title "{/Symbol s}_n' = 50 kPa"
+set title "v = 300 m/a"
+plot "strain_distribution_N8kPa_v9.506e-6ms.txt" u ($2*365*24*3600):1 w l lw 2 title "{/Symbol s}_n' = 8 kPa", \
+ "strain_distribution_N22kPa_v9.506e-6ms.txt" u ($2*365*24*3600):1 w l lw 2 title "{/Symbol s}_n' = 22 kPa", \
+ "strain_distribution_N50kPa_v9.506e-6ms.txt" u ($2*365*24*3600):1 w l lw 2 title "{/Symbol s}_n' = 50 kPa"+
\ 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
@@ -1,3 +1,3 @@
-8.53e3 0.351 0.0238
-22.2e3 0.510 0.0434
-50.4e3 0.661 0.0434
+8.53e3 0.351 0.0053 0.0238
+22.2e3 0.510 0.0134 0.0434
+50.4e3 0.661 0.0835 0.0434
diff --git a/fig-effective_stress/rs_sediment_flux_300ma.txt b/fig-effective_stress/rs_sediment_flux_300ma.txt
@@ -0,0 +1,2 @@
+21.2e3 1.51 0.153
+50.2e3 1.99 0.0959
