commit 9cd65d9f9940142b7eff3a53053172c530f6479f
parent 37920e82855d1bcd07f08d2da7b7d92fcada8b82
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date: Thu, 12 Dec 2019 20:10:45 +0100
Improve figure appearance
Diffstat:
1 file changed, 31 insertions(+), 26 deletions(-)
diff --git a/fig-skin_depth/fig.gp b/fig-skin_depth/fig.gp
@@ -3,10 +3,9 @@
reset
set terminal pdfcairo color size 20 cm, 5.5 cm
-
set multiplot layout 1,3 \
- margin 0.07,0.92,0.25,0.95 \
- spacing 0.18,0.0
+ margin 0.07,0.92,0.25,0.87 \
+ spacing 0.17,0.0
set key bottom right #samplen 0.9
@@ -25,8 +24,9 @@ f_max = 2.0/3600.0
# f: frequency, k: permeability
skindepth(f,k) = (k/(phi*mu_f*beta_f*3.141592654*f))**(0.5) # f: 1/s
-set xlabel "Forcing frequency [1/s]"
-set ylabel "Permeability [m^2]"
+set xlabel "Forcing frequency {/:Italic f} [1/s]"
+set ylabel "Permeability {/:Italic k} [m^2]"
+set cblabel "Skin depth {/:Italic d}_s [m]"
set logscale xyzcb
set xrange [f_min:f_max]
@@ -35,41 +35,44 @@ set samples 100
set isosamples 60
# Generate contours
-load 'BlueSeq.plt'
+#load 'BlueSeq.plt'
set view map
unset surface
-set contour base
-set cntrlabel format '%8.3g' font ',7' start 0.1 interval -1
-set cntrparam levels discrete 0.1, 0.3, 1, 3, 10, 30, 100, 300
-set key outside title "Skin depth\n[m]" width 1.0 samplen 2.5
+#set contour base
+#set cntrlabel format '%8.3g' font ',7' start 0.1 interval -1
+#set cntrparam levels discrete 0.1, 0.3, 1, 3, 10, 30, 100, 300
+#set key outside title "Skin depth\n[m]" width 0.5 samplen 2.5
set xtics rotate by 90
set ytics offset 1
set format x '10^{%T}'
set format y '10^{%T}'
-label_y = 3e-19
+set pm3d
+
+label_y = 3e-18
x=1.0/3600.0
-set arrow from x,k_min to x,k_max nohead lc black
-set label "hourly" at x,label_y rotate by 90 offset screen -0.01,0
+set arrow from x,k_min to x,k_max nohead lc black front
+set label "hourly" at x,label_y rotate by 90 front offset screen -0.01,0
x=1.0/3600.0/24.0
-set arrow from x,k_min to x,k_max nohead lc black
-set label "daily" at x,label_y rotate by 90 offset screen -0.01,0
+set arrow from x,k_min to x,k_max nohead lc black front
+set label "daily" at x,label_y rotate by 90 front offset screen -0.01,0
x=1.0/3600.0/24.0/30.0
-set arrow from x,k_min to x,k_max nohead lc black
-set label "monthly" at x,label_y rotate by 90 offset screen -0.01,0
+set arrow from x,k_min to x,k_max nohead lc black front
+set label "monthly" at x,label_y rotate by 90 front offset screen -0.01,0
x=1.0/3600.0/24.0/365.0
-set arrow from x,k_min to x,k_max nohead lc black
-set label "yearly" at x,label_y rotate by 90 offset screen -0.01,0
+set arrow from x,k_min to x,k_max nohead lc black front
+set label "yearly" at x,label_y rotate by 90 front offset screen -0.01,0
-splot skindepth(x,y) title ""
+splot skindepth(x,y) notitle with lines palette, \
+ skindepth(x,y) with pm3d notitle
reset
-set xlabel "Forcing frequency [1/s]"
-set ylabel "Permeability [m^2]"
+set xlabel "Forcing frequency {/:Italic f} [1/s]"
+set ylabel "Permeability {/:Italic k} [m^2]"
set cblabel "Max. deformation depth [m]"
set xrange [f_min:f_max]
@@ -91,10 +94,12 @@ set logscale xy
#unset colorbox
-set title "A_f = 10 kPa'
-splot 'max_depth_ampl10e3.txt' t ''
+set title "{/:Italic A}_f = 10 kPa"
+splot 'max_depth_ampl10e3.txt' notitle with lines palette, \
+ '' with pm3d notitle
#unset ylabel
#set colorbox
-set title "A_f = 100 kPa'
-splot 'max_depth_ampl100e3.txt' t ''
+set title "{/:Italic A}_f = 100 kPa"
+splot 'max_depth_ampl100e3.txt' notitle with lines palette, \
+ '' with pm3d notitle
