commit 35956b56e7848746b595a3114543bad9e1011d85
parent 27160583d3c58afd6af35c00406b7965e042c16c
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date: Fri, 22 Nov 2019 13:42:35 +0100
Update plots for combining in single figure
Diffstat:
3 files changed, 33 insertions(+), 25 deletions(-)
diff --git a/fig-mohr_coulomb/fig.gp b/fig-mohr_coulomb/fig.gp
@@ -4,18 +4,18 @@ reset
## resemble Iverson 2010 fig2a
-set terminal pdfcairo size 7.5 cm, 7.5 cm
+set terminal pdfcairo size 7.5 cm, 5 cm
#set terminal pdfcairo size 15 cm, 10 cm enhanced
unset xrange
unset yrange
unset logscale x
set xtics 50
-set ytics 20
+set ytics 50
set ylabel "Shear stress {/Symbol t} [kPa]"
set xlabel "Effective normal stress {/Symbol s}_n' [kPa]"
set xrange [0:200]
-set yrange [0:120]
+set yrange [0:150]
set angles degrees
caesar(x) = (x >= 50 && x <= 5000) ? tan(28)*x - 6.5 : 1/0 # Rathbun et al., 2008, saturated Caesar till
@@ -24,15 +24,17 @@ two_rivers(x) = (x >= 30 && x <= 145) ? tan(17.8)*x + 14 : 1/0 # Iverson et
# Des Moines Lobe till properties extracted from Iverson 2010 plot, mc_des_moines_lobe.dat
wis(x) = (x >= 20 && x <= 310) ? tan(24)*x + 3 : 1/0 # Tulaczyk et al., 2000
-set key top left reverse Left font ",8" samplen 2.0
+set key top left reverse Left font ",9" samplen 2.0
+
+set label "{/:Bold c}" at screen 0.05,0.95
plot caesar(x) w l lt 3 lw 2 t "Caesar till (Rathbun et al., 2008)", \
- "mc_caesar.txt" u ($1):($2*$1 - 6.5) w l dashtype 2 lt 3 lw 3 t "Simulated Caesar till", \
+ "mc_caesar.txt" u ($1):($2*$1 - 6.5) w l dashtype 2 lt 3 lw 2 t "NGF simulated Caesar till", \
storglaciaren(x) w l lt 4 lw 2 t "Storglaciaren till (Iverson et al., 1998)", \
- "mc_storglaciaren.txt" u ($1):($2*$1 + 5) w l dashtype 2 lt 4 lw 3 t "Simulated Storglaciaren till", \
+ "mc_storglaciaren.txt" u ($1):($2*$1 + 5) w l dashtype 2 lt 4 lw 2 t "NGF simulated Storglaciaren till", \
two_rivers(x) w l lt 5 lw 2 t "Two Rivers till (Iverson et al., 1998)", \
- "mc_two_rivers.txt" u ($1):($2*$1 + 14) w l dashtype 2 lt 5 lw 3 t "Simulated Two Rivers till", \
+ "mc_two_rivers.txt" u ($1):($2*$1 + 14) w l dashtype 2 lt 5 lw 2 t "NGF simulated Two Rivers till", \
wis(x) w l lt 7 lw 2 t "WIS till (Tulaczyk et al., 2000)", \
- "mc_wis.txt" u ($1):($2*$1 + 3) w l dashtype 2 lt 7 lw 3 t "Simulated WIS till"
+ "mc_wis.txt" u ($1):($2*$1 + 3) w l dashtype 2 lt 7 lw 2 t "NGF simulated WIS till"
#"mc_des_moines_lobe.dat" w l lt 6 lw 2 t "Des Moines Lobe till (Iverson, 2010)", \
- #"mc_des_moines_lobe.txt" u ($1):($2*$1 + 4.23) w l dashtype 2 lt 6 lw 3 t "Simulated Des Moines Lobe till", \
+ #"mc_des_moines_lobe.txt" u ($1):($2*$1 + 4.23) w l dashtype 2 lt 6 lw 2 t "NGF simulated Des Moines Lobe till", \
diff --git a/fig-rate_dependence/fig.gp b/fig-rate_dependence/fig.gp
@@ -13,7 +13,7 @@ set terminal pdfcairo size 7.5 cm, 5 cm enhanced
set xrange [1e1:9e5] # m/a
set yrange [-0.0:0.7]
-set key bottom right font ",8" samplen 2.0
+set key bottom right font ",9" samplen 2.0
set logscale x
#unset xlabel
@@ -24,7 +24,7 @@ set xlabel "Shear strain rate d{/Symbol g}/dt [a^{-1}]"
#set label "{/:Bold a}" at screen 0.11,0.95
#set label "{/:Bold b}" at screen 0.11,0.44
-set label "{/:Bold b}" at screen 0.11,0.95
+set label "{/:Bold b}" at screen 0.05,0.95
plot \
"cowden.dat" w l lw 2 t "Cowden till (Tika et al., 1996)", \
diff --git a/fig-strain_distribution/fig.gp b/fig-strain_distribution/fig.gp
@@ -2,32 +2,38 @@
reset
-set terminal pdfcairo enhanced color size 7.5 cm, 7.0 cm
+set terminal pdfcairo enhanced color size 7.5 cm, 5 cm
#set lmargin 7.0
#set bmargin 3.5
#set rmargin 2.0
#set tmargin 2.0
-set xlabel "Shear displacement, {/Symbol D}x [m]"
-set ylabel "Vertical position, z [m]"
+set xlabel "{/:Normal Shear displacement {/Symbol D}x [m]}"
+set ylabel "{/:Normal Vertical position z [m]}"
set yrange [0.0:0.73]
#set xrange [-0.13:1.13]
-set xrange [-0.1:0.8]
+set xrange [-0.05:0.9]
+set xtics 0.2
+set ytics 0.2
xdisp=0.74
+set label "{/:Bold d}" at screen 0.05,0.95
+
set key bottom right font ",10" #samplen 0.9
-plot "strain_distribution_N10kPa.txt" u ($2*xdisp):1 w l lw 2 lc "#ee2430" title "{/Symbol s}_n' 10 kPa", \
- "damsgaard2013-sheardisp/sheardisp10-binned.dat" w p lw 2 lc "#ee2430" title "DEM, {/Symbol s}_n' 10 kPa", \
- "strain_distribution_N20kPa.txt" u ($2*xdisp):1 w l lw 2 lc "#53bf3a" title "{/Symbol s}_n' 20 kPa", \
- "damsgaard2013-sheardisp/sheardisp20-binned.dat" w p lw 2 lc "#53bf3a" title "DEM, {/Symbol s}_n' 20 kPa", \
- "strain_distribution_N40kPa.txt" u ($2*xdisp):1 w l lw 2 lc "#2f54ad" title "{/Symbol s}_n' 40 kPa", \
- "damsgaard2013-sheardisp/sheardisp40-binned.dat" w p lw 2 lc "#2f54ad" title "DEM, {/Symbol s}_n' 40 kPa", \
- "strain_distribution_N80kPa.txt" u ($2*xdisp):1 w l lw 2 lc "#6ccfd9" title "{/Symbol s}_n' 80 kPa", \
- "damsgaard2013-sheardisp/sheardisp80-binned.dat" w p lw 2 lc "#6ccfd9" title "DEM, {/Symbol s}_n' 80 kPa", \
- #"strain_distribution_N120kPa.txt" u ($2*xdisp):1 w l lw 2 lc "#231f20" title "{/Symbol s}_n' 120 kPa", \
- #"damsgaard2013-sheardisp/sheardisp120-binned.dat" w p lw 2 lc "#231f20" title "DEM, {/Symbol s}_n' 120 kPa"
+pointsize=0.6
+
+plot "strain_distribution_N10kPa.txt" u ($2*xdisp):1 w l lw 2 lc "#ee2430" title "NGF, {/Symbol s}_n' 10 kPa", \
+ "damsgaard2013-sheardisp/sheardisp10-binned.dat" w p ps pointsize lw 2 lc "#ee2430" title "DEM, {/Symbol s}_n' 10 kPa", \
+ "strain_distribution_N20kPa.txt" u ($2*xdisp):1 w l lw 2 lc "#53bf3a" title "NGF, {/Symbol s}_n' 20 kPa", \
+ "damsgaard2013-sheardisp/sheardisp20-binned.dat" w p ps pointsize lw 2 lc "#53bf3a" title "DEM, {/Symbol s}_n' 20 kPa", \
+ "strain_distribution_N40kPa.txt" u ($2*xdisp):1 w l lw 2 lc "#2f54ad" title "NGF, {/Symbol s}_n' 40 kPa", \
+ "damsgaard2013-sheardisp/sheardisp40-binned.dat" w p ps pointsize lw 2 lc "#2f54ad" title "DEM, {/Symbol s}_n' 40 kPa", \
+ "strain_distribution_N80kPa.txt" u ($2*xdisp):1 w l lw 2 lc "#6ccfd9" title "NGF, {/Symbol s}_n' 80 kPa", \
+ "damsgaard2013-sheardisp/sheardisp80-binned.dat" w p ps pointsize lw 2 lc "#6ccfd9" title "DEM, {/Symbol s}_n' 80 kPa", \
+ #"strain_distribution_N120kPa.txt" u ($2*xdisp):1 w l lw 2 lc "#231f20" title "NGF, {/Symbol s}_n' 120 kPa", \
+ #"damsgaard2013-sheardisp/sheardisp120-binned.dat" w p ps pointsize lw 2 lc "#231f20" title "DEM, {/Symbol s}_n' 120 kPa"
