commit ba8bc65b74e3f4cc95c7295612747f1ad33c2d41
parent afbf73f6a49283ea015b49757b29ec515a588fff
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date: Thu, 6 Oct 2022 17:11:01 +0200
plot A_f* threshold
Diffstat:
2 files changed, 8 insertions(+), 1 deletion(-)
diff --git a/Makefile b/Makefile
@@ -130,7 +130,8 @@ fig-sediment_flux-100ma.pdf: fig-sediment_flux-100ma.gp sediment_flux_v3.169e-6m
gnuplot fig-sediment_flux-100ma.gp > $@
fig-sediment_flux_daily.pdf: fig-sediment_flux_daily.gp sediment_flux_freq1.157e-5s.txt max_depth_daily.txt
- gnuplot fig-sediment_flux_daily.gp > $@
+ skindepth=$$(awk '{print $$3; exit}' max_depth_daily.txt) &&\
+ gnuplot -e "d_s=$${skindepth}; rho_s=2600; rho_f=1000; g=9.81" fig-sediment_flux_daily.gp > $@
fig-skin_depth_diff.pdf: fig-skin_depth_diff.gp max_depth_ampl10e3.txt
gnuplot fig-skin_depth_diff.gp > $@
diff --git a/fig-sediment_flux_daily.gp b/fig-sediment_flux_daily.gp
@@ -19,6 +19,8 @@ set label "b" at screen labelx,0.48 font "{:Bold},12"
set xrange [0:*]
set yrange [0:*]
+A_f_crit = (rho_s - rho_f) * g * d_s**3.0 / 1000.0 # [kPa]
+
#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]"
@@ -26,8 +28,12 @@ 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:*]
+set arrow from A_f_crit,0 to A_f_crit,30 nohead lc "black" front
plot "sediment_flux_freq1.157e-5s.txt" u ($1):($2*3600*24*365.25*rswidth) w lp t ""
+unset arrow
+
set yrange [-0.1:1]
+set arrow from A_f_crit,-0.1 to A_f_crit,1 nohead lc "black" front
set ylabel "Max. expected deformation depth [m]"
plot "max_depth_daily.txt" u ($1/1000):($2) w l t ''
