cngf-pf-exp1

commit 276bbe01187fd68cec92c5425609ad171a3dbe64
parent e76e16e4097146f0ced3870adae4c9bc617b40c1
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Mon, 15 Jul 2019 15:02:28 +0200

Update 1d_fd_simple_shear and prepare fig1 for MC test

Diffstat:
M
fig1/Makefile
 | 
15
++++++++++++++-
M
fig1/fig.gp
 | 
23
+++++++++++------------

2 files changed, 25 insertions(+), 13 deletions(-)
diff --git a/fig1/Makefile b/fig1/Makefile
@@ -17,7 +17,20 @@ sim.output00000.txt: $(BIN) Makefile
 			tail -n 1 | cut -f2 >> "$$out"; \
 	done; done'
 
-../$(FIG).pdf: fig.gp sim.output00000.txt
+# mimic Iverson 2010 fig 2b
+mc.output00000.txt: $(BIN) Makefile
+	/bin/sh -c '\
+	for mu in $$(seq 0.30 0.05 0.50); do \
+		for C in $$(printf "0\n20e3"); do \
+			out="out_mu$${mu}_C$${C}.txt"; \
+			rm -f "$$out"; \
+			for N in $$(seq 10 20 200); do \
+				printf "$$N\t" >> "$$out"; \
+				./$(BIN) -m $$mu -f $$mu -C $$C -P $$N -L 1.0 mc | \
+				tail -n 1 | cut -f4 >> "$$out"; \
+	done; done; done'
+
+../$(FIG).pdf: fig.gp sim.output00000.txt mc.output00000.txt
 	gnuplot $< > $@
 
 clean:
diff --git a/fig1/fig.gp b/fig1/fig.gp
@@ -4,23 +4,25 @@ reset
 
 ## resemble Iverson 2010 fig2a
 
-set terminal pdfcairo size 7.5 cm, 7.5 cm
+#set terminal pdfcairo size 7.5 cm, 7.5 cm
+set terminal pdfcairo size 15 cm, 7.5 cm
 
-set multiplot layout 2,1 \
+set multiplot layout 2,2 \
 	margins 0.18,0.93,0.15,0.95 \
 	spacing 0.03,0.05
 
-set xrange [1e1:1e6]  # m/a
+set xrange [1e1:9e5]  # m/a
 set yrange [0.0:0.8]
 set key bottom right font ",7" samplen 2.0
 set logscale x
 
-unset xlabel
-unset xtics
+#unset xlabel
+#unset xtics
 set ylabel "Shear friction [-]"
+set xlabel "Shear strain rate [a^{-1}]"
 
-set label "{/:Bold a}" at screen 0.01,0.97
-set label "{/:Bold b}" at screen 0.01,0.54
+#set label "{/:Bold a}" at screen 0.01,0.97
+#set label "{/:Bold b}" at screen 0.01,0.54
 
 plot "cowden.dat"         w l lw 2 t "Cowden till (Tika et al., 1996)", \
 	"lower_cromer.dat"    w l lw 2 t "Lower Cromer till (Tika et al., 1996)", \
@@ -30,14 +32,11 @@ plot "cowden.dat"         w l lw 2 t "Cowden till (Tika et al., 1996)", \
 	"des_moines_lobe.dat" w l lw 2 t "Des Moines Lobe till (Iverson, 2010)", \
 	"wis.dat"             w l lw 2 t "WIS till (Tulaczyk et al., 2001)"
 
-#set label "a" at screen 0.01, 0.98
-#set label '{/Arial-Bold attempt to specify}' at screen 0.01 0.98
-
 set key bottom right font ",10" samplen 2.0
 
 set xtics
-set xlabel "Shear strain rate [a^{-1}]"
-set ylabel "Shear friction [-]"
+unset ylabel
+unset ytics
 
 load 'BlueSeq.plt'