manus_continuum_granular1_exp

commit 609885051096465d247666bdf09ad9506bba047c
parent 3af78234c95986555fab8e628d4bfda4683a4637
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Mon,  8 Jul 2019 17:40:21 +0200

Add simulation of different real tills

Diffstat:
A
fig7/Makefile
 | 
47
+++++++++++++++++++++++++++++++++++++++++++++++
A
fig7/fig.gp
 | 
23
+++++++++++++++++++++++

2 files changed, 70 insertions(+), 0 deletions(-)
diff --git a/fig7/Makefile b/fig7/Makefile
@@ -0,0 +1,47 @@
+BIN = ../1d_fd_simple_shear/1d_fd_simple_shear
+FIG = $(shell basename $(CURDIR))
+
+RESOLUTION := 1000
+LENGTH := 1.0
+NORMAL_STRESS := 85e3
+
+default: ../$(FIG).pdf
+
+two_rivers.txt: $(BIN) Makefile
+	/bin/sh -c '\
+	./$< --resolution $(RESOLUTION) \
+	--length $(LENGTH) \
+	--normal-stress $(NORMAL_STRESS) \
+	--friction-coefficient 0.321 --cohesion 14.0e3 \
+	--stress-ratio 0.5 \
+	--grain-size 0.1e-3 \
+	--normalize' > $@
+
+storglaciaren.txt: $(BIN) Makefile
+	/bin/sh -c '\
+	./$< --resolution $(RESOLUTION) \
+	--length $(LENGTH) \
+	--normal-stress $(NORMAL_STRESS) \
+	--friction-coefficient 0.494 --cohesion 5.0e3 \
+	--stress-ratio 0.594 \
+	--grain-size 0.1e-3 \
+	--normalize' > $@
+
+upb.txt: $(BIN) Makefile
+	/bin/sh -c '\
+	./$< --resolution $(RESOLUTION) \
+	--length $(LENGTH) \
+	--normal-stress $(NORMAL_STRESS) \
+	--friction-coefficient 0.443 --cohesion 3.0e3 \
+	--stress-ratio 0.543 \
+	--grain-size 0.1e-3 \
+	--normalize' > $@
+
+../$(FIG).pdf: fig.gp two_rivers.txt storglaciaren.txt upb.txt
+	gnuplot $< > $@
+
+clean:
+	$(RM) *.txt
+	$(RM) ../$(FIG).pdf
+
+.PHONY: default clean
diff --git a/fig7/fig.gp b/fig7/fig.gp
@@ -0,0 +1,23 @@
+#!/usr/bin/env gnuplot
+
+reset
+
+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,4 \
+#	margins 0.09,0.97,0.23,0.95 \
+#	spacing 0.03,0.03
+
+#set yrange [4.0:8.0]
+
+set key bottom right #samplen 0.9
+
+set xlabel "Normalized shear velocity [-]"
+set ylabel "Vertical position [m]"
+#set xrange [0.0:1.0]
+#set xtics (0,0.01,0.02)
+#unset ytics
+plot \
+	"two_rivers.txt" u 2:1 w l lw 1 t "simulated Two Rivers till", \
+	"storglaciaren.txt" u 2:1 w l lw 1 t "simulated Storglaciaren till", \
+	"upb.txt" u 2:1 w l lw 1 t "simulated Upstream B till"