hansen-zoet-exp

commit 8c7aaae09ec24d954876051cff2c81b41d781c90
parent b28587911a8273ad379205b63a743e2d3e37d44e
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Tue, 14 Jun 2022 11:56:03 +0200

plot sediment flux under daily oscillations together with max depth solution

Diffstat:
M
Makefile
 | 
16
+++++++++-------
M
fig-sediment_flux_daily.gp
 | 
18
+++++++++++-------

2 files changed, 20 insertions(+), 14 deletions(-)
diff --git a/Makefile b/Makefile
@@ -13,7 +13,9 @@ 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 16 1 30
+fluid_ampl_min = 1
+fluid_ampl_max = 30
+fluid_ampl_vals != range -N 16 ${fluid_ampl_min} ${fluid_ampl_max}
 fluid_freq_vals = 1.157e-5	# daily
 grain_size = 530.7e-6
 grain_friction = 0.6
@@ -52,16 +54,15 @@ strain_distribution_N50kPa_v3.169e-6ms.txt: ${BIN}
 	done; done'  # also close fit with lower A and higher b
 	# higher values of L do not make a difference in flux
 
-.PHONY: max_depth
-
-max_depth: ${BIN}
+max_depth_daily.txt: ${BIN}
 	for freq in ${fluid_freq_vals}; do \
-		for ampl in ${fluid_ampl_vals}; do \
+		for ampl in $$(range -N 1000 ${fluid_ampl_min} ${fluid_ampl_max}); do \
+			printf '%ge3\t' $${ampl}; \
 			./${REPO}/max_depth_simple_shear \
 			-D 1.5e-5 \
 			-q $${freq} \
 			-a $${ampl}e3; \
-	done; done
+	done; done >$@
 
 strain_distribution_ampl1kPa_freq1.157e-5s.txt: ${BIN}
 	sh -c '\
@@ -126,7 +127,8 @@ 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
+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 > $@
 
 clean:
diff --git a/fig-sediment_flux_daily.gp b/fig-sediment_flux_daily.gp
@@ -2,15 +2,15 @@
 
 reset
 
-set terminal pdfcairo enhanced color size 7.5 cm, 6 cm
-set multiplot layout 1,1
+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 lmargin 8.0
+set bmargin 3.5
+set rmargin 2.5
+set tmargin 2.0
 
 #set xrange [-5 : 105]
 #set yrange [0 : 0.11]
@@ -30,4 +30,8 @@ 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:*]
-plot	"sediment_flux_freq1.157e-5s.txt" u ($1):($2*3600*24*365.25*rswidth) w lp ps 0.5 t ""
+plot	"sediment_flux_freq1.157e-5s.txt" u ($1):($2*3600*24*365.25*rswidth) w lp t ""
+
+set yrange [-0.1:1]
+set ylabel "Max. expected deformation depth [m]"
+plot "max_depth_daily.txt" u ($1/1000):($2) w l t ''