cngf-pf-exp1

commit 879eda3514b0f30b2bfdc5ef873c856573559b92
parent 276bbe01187fd68cec92c5425609ad171a3dbe64
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Mon, 15 Jul 2019 15:44:15 +0200

Add Mohr-Coulomb plot from model

Diffstat:
A
fig1/RedSeq.plt
 | 
26
++++++++++++++++++++++++++
M
fig1/fig.gp
 | 
27
++++++++++++++++++++++++---

2 files changed, 50 insertions(+), 3 deletions(-)
diff --git a/fig1/RedSeq.plt b/fig1/RedSeq.plt
@@ -0,0 +1,25 @@
+# line styles for ColorBrewer Reds
+# for use with sequential data
+# provides 8 red colors of increasing saturation
+# compatible with gnuplot >=4.2
+# author: Anna Schneider
+
+# line styles
+set style line 1 lc rgb '#FFF5F0' # very light red
+set style line 2 lc rgb '#FEE0D2' # 
+set style line 3 lc rgb '#FCBBA1' # 
+set style line 4 lc rgb '#FC9272' # light red
+set style line 5 lc rgb '#FB6A4A' # 
+set style line 6 lc rgb '#EF3B2C' # medium red
+set style line 7 lc rgb '#CB181D' #
+set style line 8 lc rgb '#99000D' # dark red
+
+# palette
+set palette defined ( 0 '#FFF5F0',\
+    	    	      1 '#FEE0D2',\
+		      2 '#FCBBA1',\
+		      3 '#FC9272',\
+		      4 '#FB6A4A',\
+		      5 '#EF3B2C',\
+		      6 '#CB181D',\
+		      7 '#99000D' )+
\ No newline at end of file
diff --git a/fig1/fig.gp b/fig1/fig.gp
@@ -5,11 +5,11 @@ reset
 ## resemble Iverson 2010 fig2a
 
 #set terminal pdfcairo size 7.5 cm, 7.5 cm
-set terminal pdfcairo size 15 cm, 7.5 cm
+set terminal pdfcairo size 15 cm, 10 cm
 
 set multiplot layout 2,2 \
 	margins 0.18,0.93,0.15,0.95 \
-	spacing 0.03,0.05
+	spacing 0.03,0.15
 
 set xrange [1e1:9e5]  # m/a
 set yrange [0.0:0.8]
@@ -34,7 +34,6 @@ plot "cowden.dat"         w l lw 2 t "Cowden till (Tika et al., 1996)", \
 
 set key bottom right font ",10" samplen 2.0
 
-set xtics
 unset ylabel
 unset ytics
 
@@ -45,3 +44,25 @@ plot "out_b0.94.txt" u ($2*3600*24*365):1 w l lc '#666666' lw 2 dashtype '-' t "
          "out_b0.20.txt" u ($2*3600*24*365):1 w l ls 5 lw 2 t "b = 0.20", \
          "out_b0.10.txt" u ($2*3600*24*365):1 w l ls 3 lw 2 t "b = 0.10", \
          "out_b0.01.txt" u ($2*3600*24*365):1 w l ls 1 lw 2 t "b = 0.01"
+
+unset xrange
+unset yrange
+unset logscale x
+set xtics 50
+set ytics 20
+set ylabel "Shear stress [kPa]"
+set xlabel "Effective normal stress [kPa]"
+
+load 'BlueSeq.plt'
+set key top left reverse Left font ",7" samplen 2.0
+
+plot "out_mu0.50_C0.txt" u ($1):($2*$1) w l ls 5 lw 2 t "mu_s = 0.50, C = 0 kPa", \
+     "out_mu0.45_C0.txt" u ($1):($2*$1) w l ls 4 lw 2 t "", \
+     "out_mu0.40_C0.txt" u ($1):($2*$1) w l ls 3 lw 2 t "mu_s = 0.40, C = 0 kPa", \
+     "out_mu0.35_C0.txt" u ($1):($2*$1) w l ls 2 lw 2 t "", \
+     "out_mu0.30_C0.txt" u ($1):($2*$1) w l ls 1 lw 2 t "mu_s = 0.30, C = 0 kPa", \
+     "out_mu0.50_C20e3.txt" u ($1):($2*$1 + 20) w l lc '#FEE0D2' lw 2 t "mu_s = 0.30, C = 20 kPa", \
+     "out_mu0.45_C20e3.txt" u ($1):($2*$1 + 20) w l lc '#FCBBA1' lw 2 t "", \
+     "out_mu0.40_C20e3.txt" u ($1):($2*$1 + 20) w l lc '#FC9272' lw 2 t "mu_s = 0.40, C = 20 kPa", \
+     "out_mu0.35_C20e3.txt" u ($1):($2*$1 + 20) w l lc '#FB6A4A' lw 2 t "", \
+     "out_mu0.30_C20e3.txt" u ($1):($2*$1 + 20) w l lc '#EF3B2C' lw 2 t "mu_s = 0.50, C = 20 kPa"