cngf-pf-exp1

commit dcb8d883627d6d311f6e49597ecf8eac0b6e23dd
parent a6899c4bc6bad65929b444ef98df3d24c570d7c2
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Mon,  8 Jul 2019 18:01:19 +0200

Tweak applied friction to match observed shear velocity of 320 m/a

Diffstat:
M
fig2/Makefile
 | 
15
++++++---------
M
fig2/fig.gp
 | 
8
++++----

2 files changed, 10 insertions(+), 13 deletions(-)
diff --git a/fig2/Makefile b/fig2/Makefile
@@ -13,9 +13,8 @@ two_rivers.txt: $(BIN) Makefile
 	--length $(LENGTH) \
 	--normal-stress $(NORMAL_STRESS) \
 	--friction-coefficient 0.321 --cohesion 14.0e3 \
-	--stress-ratio 0.5 \
-	--grain-size 0.1e-3 \
-	--normalize' > $@
+	--stress-ratio 0.536 \
+	--grain-size 0.1e-3' > $@
 
 storglaciaren.txt: $(BIN) Makefile
 	/bin/sh -c '\
@@ -23,9 +22,8 @@ storglaciaren.txt: $(BIN) Makefile
 	--length $(LENGTH) \
 	--normal-stress $(NORMAL_STRESS) \
 	--friction-coefficient 0.494 --cohesion 5.0e3 \
-	--stress-ratio 0.594 \
-	--grain-size 0.1e-3 \
-	--normalize' > $@
+	--stress-ratio 0.615 \
+	--grain-size 0.1e-3' > $@
 
 upb.txt: $(BIN) Makefile
 	/bin/sh -c '\
@@ -33,9 +31,8 @@ upb.txt: $(BIN) Makefile
 	--length $(LENGTH) \
 	--normal-stress $(NORMAL_STRESS) \
 	--friction-coefficient 0.443 --cohesion 3.0e3 \
-	--stress-ratio 0.543 \
-	--grain-size 0.1e-3 \
-	--normalize' > $@
+	--stress-ratio 0.537 \
+	--grain-size 0.1e-3' > $@
 
 ../$(FIG).pdf: fig.gp two_rivers.txt storglaciaren.txt upb.txt
 	gnuplot $< > $@
diff --git a/fig2/fig.gp b/fig2/fig.gp
@@ -12,12 +12,12 @@ set terminal pdfcairo color size 7.5 cm, 5.0 cm
 
 set key bottom right #samplen 0.9
 
-set xlabel "Normalized shear velocity [-]"
+set xlabel "Shear velocity [m/a]"
 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 Storglaciären till", \
-	"upb.txt" u 2:1 w l lw 1 t "simulated Upstream B till"
+	"two_rivers.txt" u ($2*365.25*24*60*60):1 w l lw 1 t "simulated Two Rivers till", \
+	"storglaciaren.txt" u ($2*365.25*24*60*60):1 w l lw 1 t "simulated Storglaciären till", \
+	"upb.txt" u ($2*365.25*24*60*60):1 w l lw 1 t "simulated Upstream B till"