commit 3e0a9aba9099de20853c326ead36debc1a17e410
parent ee16d066a3bded74a642e90a6864c43a6f9d6900
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date: Mon, 13 Jun 2022 13:47:06 +0200
add first flux figure with oscillations
Diffstat:
2 files changed, 38 insertions(+), 2 deletions(-)
diff --git a/Makefile b/Makefile
@@ -12,7 +12,7 @@ stressvals_ != range -l -N 40 0.01 100
stressvals_ += 8 22 50
stressvals != echo ${stressvals_} | sort -n
velvals = 3.169e-6 9.506e-6
-fluid_ampl_vals != range -N 4 1 80
+fluid_ampl_vals != range -N 4 1 40
fluid_freq_vals = 1.157e-5 # daily
grain_size = 530.7e-6
grain_friction = 0.6
@@ -62,7 +62,7 @@ strain_distribution_ampl1kPa_freq1.157e-5s.txt: ${BIN}
-c ${grain_cohesion} \
-s 3.169e-6 \
-L 0.5 \
- -n 100e3 \
+ -n 85e3 \
-A ${grain_nonlocal_ampl} \
-b ${grain_rate_dependence} \
-F \
@@ -112,6 +112,9 @@ fig-sediment_flux.pdf: fig-sediment_flux.gp sediment_flux_v3.169e-6ms.txt
fig-sediment_flux-100ma.pdf: fig-sediment_flux-100ma.gp sediment_flux_v3.169e-6ms.txt
gnuplot fig-sediment_flux-100ma.gp > $@
+fig-sediment_flux_daily.pdf: fig-sediment_flux_daily.gp sediment_flux_freq1.157e-5s.txt
+ gnuplot fig-sediment_flux_daily.gp > $@
+
clean:
rm -f strain_distribution_*.txt mohr_coulomb_*.txt sediment_flux_*.txt
rm -f ${PDFS}
diff --git a/fig-sediment_flux_daily.gp b/fig-sediment_flux_daily.gp
@@ -0,0 +1,33 @@
+#!/usr/bin/env gnuplot
+
+reset
+
+set terminal pdfcairo enhanced color size 7.5 cm, 10 cm
+set multiplot layout 2,1
+
+rswidth = 1
+
+#set lmargin 7.0
+#set bmargin 3.5
+#set rmargin 2.0
+#set tmargin 2.0
+
+#set xrange [-5 : 105]
+#set yrange [0 : 0.11]
+
+# 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 xrange [0:*]
+set yrange [0:*]
+
+#set xlabel "Normalized horizontal velocity, v_x [-]"
+set xlabel "Amplitude of water-pressure fluctuations, {/:Italic A'}_f [kPa]"
+#set ylabel "Spec. sediment flux [m²/a]"
+set ylabel "Sediment flux, {/:Italic q}_t [m³/a]"
+set key bottom right font ",10" invert #samplen 0.9
+set title "{/:Italic v} = 100 m/a"
+#set yrange [0:1]
+plot "sediment_flux_freq1.157e-5s.txt" u ($1/1000):($2*3600*24*365.25*rswidth) w p ps 0.5 t "CNGF-PF"