commit 4afb091c0cf8494a8ea5aaae3ac2704f58001457
parent d47e95b4ebcc3323cd9a824252dd6ecde5c2931d
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date: Tue, 14 Jun 2022 15:01:05 +0200
add labels to subplots
Diffstat:
3 files changed, 26 insertions(+), 29 deletions(-)
diff --git a/fig-sediment_flux.gp b/fig-sediment_flux.gp
@@ -7,13 +7,14 @@ set multiplot layout 2,1
rswidth = 1
-#set lmargin 7.0
-#set bmargin 3.5
-#set rmargin 2.0
-#set tmargin 2.0
+set lmargin 8.0
+set bmargin 3.5
+set rmargin 2.0
+set tmargin 2.0
-#set xrange [-5 : 105]
-#set yrange [0 : 0.11]
+labelx = 0.05
+set label "a" at screen labelx,0.98 font "{:Bold},12"
+set label "b" at screen labelx,0.48 font "{:Bold},12"
# https://sodocumentation.net/gnuplot/topic/8825/fit-data-with-gnuplot
#f(x) = a * exp(b * x)
@@ -21,22 +22,22 @@ f(x, v) = a * x**b * v
fit f(x, 100) "sediment_flux_v3.169e-6ms.txt" u ($1/1000):($2*3600*24*365.25*rswidth) via a,b
fitparams = sprintf("{/:Italic q}_t = %.2e {/:Italic v N'}^{%.3g}", a, b)
-set xrange [0:*]
-set yrange [0:*]
-
-#set xlabel "Normalized horizontal velocity, v_x [-]"
set xlabel "Effective stress, {/:Italic N'} [kPa]"
-#set ylabel "Spec. sediment flux [m²/a]"
set ylabel "Sediment flux, {/:Italic q}_t [m³/a]"
-set key bottom right font ",10" invert #samplen 0.9
+set key bottom right font ",09" invert samplen 1.0
+#set title "{/:Italic v} = 100 m/a"
+set xrange [0:210]
+set xtics 0,50,200
+
set title "{/:Italic v} = 100 m/a"
-set yrange [0:1]
-plot f(x, 100) t fitparams, \
- "sediment_flux_v3.169e-6ms.txt" u ($1/1000):($2*3600*24*365.25*rswidth) w p ps 0.5 t "CNGF-PF", \
- "rs_sediment_flux_100ma.txt" u ($1/1000):($2) w p ps 1 t "Hansen and Zoet (2022)"
+set yrange [0:1.1]
+plot f(x, 100) lt 1 t fitparams, \
+ "sediment_flux_v3.169e-6ms.txt" u ($1/1000):($2*3600*24*365.25*rswidth) w p lt 2 ps 0.5 t "CNGF-PF (Damsgaard et al., 2020)", \
+ "rs_sediment_flux_100ma_high-N.txt" u ($1/1000):($2) w p lt 4 ps 1 t "Hansen and Zoet (2022)", \
+ "rs_sediment_flux_100ma.txt" u ($1/1000):($2) w p lt 4 ps 1 t ""
set title "{/:Italic v} = 300 m/a"
set yrange [0:*]
-plot f(x, 300) t fitparams, \
- "sediment_flux_v9.506e-6ms.txt" u ($1/1000):($2*3600*24*365.25*rswidth) w p ps 0.5 t "CNGF-PF", \
- "rs_sediment_flux_300ma.txt" u ($1/1000):($2) w p ps 1 t "Hansen and Zoet (2022)"
+plot f(x, 300) lt 1 t fitparams, \
+ "sediment_flux_v9.506e-6ms.txt" u ($1/1000):($2*3600*24*365.25*rswidth) w p lt 2 ps 0.5 t "CNGF-PF (Damsgaard et al., 2020)", \
+ "rs_sediment_flux_300ma.txt" u ($1/1000):($2) w p lt 4 ps 1 t "Hansen and Zoet (2022)"
diff --git a/fig-sediment_flux_daily.gp b/fig-sediment_flux_daily.gp
@@ -12,13 +12,9 @@ set bmargin 3.5
set rmargin 2.5
set tmargin 2.0
-#set xrange [-5 : 105]
-#set yrange [0 : 0.11]
-
-# https://sodocumentation.net/gnuplot/topic/8825/fit-data-with-gnuplot
-#f(x, v) = a * x**b * v
-#fit f(x, 100) "sediment_flux_v3.169e-6ms.txt" u ($1/1000):($2*3600*24*365.25*rswidth) via a,b
-#fitparams = sprintf("{/:Italic q}_t = %.2e {/:Italic v N'}^{%.3g}", a, b)
+labelx = 0.08
+set label "a" at screen labelx,0.98 font "{:Bold},12"
+set label "b" at screen labelx,0.48 font "{:Bold},12"
set xrange [0:*]
set yrange [0:*]
diff --git a/fig-skin_depth_diff.gp b/fig-skin_depth_diff.gp
@@ -74,9 +74,9 @@ x=1.0/3600.0/24.0/365.0
set arrow from x,D_min to x,D_max nohead lc "white" front
set label "yearly" at x,label_y rotate by 90 front offset screen -0.01,0 textcolor "white"
-set label "a" at screen 0.01,0.95 font ",12"
-set label "b" at screen 0.35,0.95 font ",12"
-set label "c" at screen 0.70,0.95 font ",12"
+set label "a" at screen 0.01,0.95 font "{:Bold},12"
+set label "b" at screen 0.35,0.95 font "{:Bold},12"
+set label "c" at screen 0.70,0.95 font "{:Bold},12"
set xlabel "Forcing frequency {/:Italic f} [s^{-1}]"
set ylabel "Diffusivity {/:Italic D} [m^2/s]"
