manus_continuum_granular1_exp

commit 103219aff8e6a91f510755e753f3efac5367d2af
parent 132ce14d998edf2805b102cbb3cc527a080ac436
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Tue, 25 Jun 2019 14:39:53 +0200

Add preliminary version of figure 3

Diffstat:
M
fig3/Makefile
 | 
6
+++---
M
fig3/fig.gp
 | 
55
+++++++++++++++++++++++++++++++++++--------------------

2 files changed, 38 insertions(+), 23 deletions(-)
diff --git a/fig3/Makefile b/fig3/Makefile
@@ -3,15 +3,15 @@ FIG = fig3
 
 default: ../$(FIG).pdf
 
-unnamed.output00000.txt: $(BIN)
+unnamed.output00000.txt: $(BIN) Makefile
 	/bin/sh -c '\
 	./$< --resolution 50 --length 2.0 --normal-stress 150e3 \
 	--fluid --fluid-permeability 2e-17 --fluid-pressure-top 50e3 \
 	--fluid-pressure-ampl 50e3 \
 	--fluid-pressure-freq $$( echo "1.0/(3600*24)" | bc -l ) \
 	--time-step 1e-1 \
-	--file-interval $$( echo "60*10" | bc -l ) \
-	--time-end $$( echo "3600*24*7" | bc -l )'
+	--file-interval $$( echo "3600" | bc -l ) \
+	--time-end $$( echo "3600*24*3 + 3600" | bc -l )'
 
 ../$(FIG).pdf: fig.gp unnamed.output00000.txt
 	gnuplot $< > $@
diff --git a/fig3/fig.gp b/fig3/fig.gp
@@ -1,33 +1,48 @@
 #!/usr/bin/env gnuplot
 
-set terminal pdfcairo color size 7.5 cm, 7.5 cm
+#set terminal pdfcairo color size 7.5 cm, 7.5 cm
+set terminal pdfcairo color size 15.0 cm, 5.0 cm
 set multiplot layout 1,3
 
 set yrange [0.0:2.0]
 
 set key bottom right #samplen 0.9
 
-
-do for [t=1:10] {
+do for [t=47:71] {
 	filename = sprintf('unnamed.output%05d.txt', t)
-	set xlabel "Water pressure, p_f [Pa]"
-	set ylabel "Vertical position, z [m]" offset 2
-	#set xrange [0:150]
-	plot filename u 4:1 w l lw 2 lc "blue" t ""
-
-	set xlabel "Effective normal stress [Pa]"
-	set ylabel ""
-	#set xrange [0:200e3]
-	plot filename u 3:1 w l lw 2 lc "black" t ""
-
-	# set xlabel "Friction, mu [-]"
+	if (t == 71) {
+		set xlabel "Water pressure [kPa]"
+		set ylabel "Vertical position, z [m]"
+	} else {
+		set xlabel " "
+		set ylabel " "
+	}
+	set xrange [0:100]
+	plot filename u ($4/1e3):1 w l lw 1 lc "blue" t ""
+
+	if (t == 71) {
+		set xlabel "Effective normal stress [kPa]"
+		set ylabel ""
+	} else {
+		set xlabel " "
+		set ylabel ""
+	}
+	set xrange [50:150]
+	plot filename u ($3/1e3):1 w l lw 1 lc "black" t ""
+
+	# set xlabel "Friction [-]"
 	# set xrange [-0.05:1.5]
-	# plot filename u 5:1 w l lw 2 lc "gray" t ""
-
-	set xlabel "Horizontal velocity, v_x [m/s]"
-	set ylabel ""
-	#set xrange [0.0:0.025]
-	plot filename u 2:1 w l lw 2 lc "red" t ""
+	# plot filename u 5:1 w l lw 1 lc "gray" t ""
+
+	if (t == 71) {
+		set xlabel "Shear velocity [m/s]"
+		set ylabel ""
+	} else {
+		set xlabel " "
+		set ylabel ""
+	}
+	set xrange [0.0:0.023]
+	plot filename u 2:1 w l lw 1 lc "red" t ""
 }
 
 unset multiplot