cngf-pf-exp1

commit 37920e82855d1bcd07f08d2da7b7d92fcada8b82
parent 701ec99ae9cf8d6de043c34c2a7e3abb99577f36
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Thu, 12 Dec 2019 17:04:10 +0100

Improve plot appearance

Diffstat:
M
fig-skin_depth/Makefile
 | 
3
+--
M
fig-skin_depth/fig.gp
 | 
44
+++++++++++++++++++++++++++-----------------

2 files changed, 28 insertions(+), 19 deletions(-)
diff --git a/fig-skin_depth/Makefile b/fig-skin_depth/Makefile
@@ -3,8 +3,7 @@ BIN = ../1d_fd_simple_shear/max_depth_simple_shear
 
 default: ../$(FIG).pdf
 
-../$(FIG).pdf: fig.gp max_depth_ampl10e3.txt #triangle.gp
-	#gnuplot triangle.gp > $@
+../$(FIG).pdf: fig.gp max_depth_ampl10e3.txt
 	gnuplot fig.gp > $@
 
 max_depth_ampl10e3.txt: $(BIN) Makefile
diff --git a/fig-skin_depth/fig.gp b/fig-skin_depth/fig.gp
@@ -2,17 +2,11 @@
 
 reset
 
-#set terminal pdfcairo color size 15 cm, 5.5 cm
-set terminal pdfcairo color size 20 cm, 10 cm
+set terminal pdfcairo color size 20 cm, 5.5 cm
 
 set multiplot layout 1,3 \
-	#margins 0.19,0.72,0.10,0.98 \
-	#spacing 0.00,0.03
-
-#set lmargin at screen 0.20
-#set rmargin at screen 0.75
-#set bmargin at screen 0.25
-#set tmargin at screen 0.95
+	margin 0.07,0.92,0.25,0.95 \
+	spacing 0.18,0.0
 
 set key bottom right #samplen 0.9
 
@@ -56,35 +50,51 @@ set format y '10^{%T}'
 label_y = 3e-19
 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.03,0
+set label "hourly" at x,label_y rotate by 90 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.03,0
+set label "daily" at x,label_y rotate by 90 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.03,0
+set label "monthly" at x,label_y rotate by 90 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.03,0
+set label "yearly" at x,label_y rotate by 90 offset screen -0.01,0
 
 splot skindepth(x,y) title ""
 
 reset
 
+set xlabel "Forcing frequency [1/s]"
+set ylabel "Permeability [m^2]"
+set cblabel "Max. deformation depth [m]"
+
+set xrange [f_min:f_max]
+set yrange [k_min:k_max]
+set xtics rotate by 90
+set ytics offset 1
+set format x '10^{%T}'
+set format y '10^{%T}'
+
 set view map scale 1
 set style data pm3d
 set style function pm3d
 set xyplane relative 0
 set nomcbtics
-set xlabel "X LABEL" 
-set ylabel "Y LABEL" 
 set pm3d implicit at b
 set pm3d scansforward
 #set pm3d interpolate 10,10
 set logscale xy
 
-splot 'max_depth_ampl10e3.txt' t 'A_f = 10 kPa'
-splot 'max_depth_ampl100e3.txt' t 'A_f = 100 kPa'
+#unset colorbox
+
+set title "A_f = 10 kPa'
+splot 'max_depth_ampl10e3.txt' t ''
+
+#unset ylabel
+#set colorbox
+set title "A_f = 100 kPa'
+splot 'max_depth_ampl100e3.txt' t ''