cngf-pf-exp1

commit 0b1afac2755409dde80734a7ae30ea3750c5fd34
parent 397afc1949061c21a54f1036bc0f208613cb0f68
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Thu, 10 Oct 2019 16:00:16 +0200

Annotate skin depth figure

Diffstat:
M
fig8/fig.gp
 | 
35
++++++++++++++++++++++++++++-------

1 file changed, 28 insertions(+), 7 deletions(-)
diff --git a/fig8/fig.gp b/fig8/fig.gp
@@ -2,11 +2,11 @@
 
 reset
 
-set terminal pdfcairo color size 7.5 cm, 6.0 cm
+set terminal pdfcairo color size 7.5 cm, 5.5 cm
 
-set lmargin at screen 0.24
+set lmargin at screen 0.20
 set rmargin at screen 0.75
-set bmargin at screen 0.20
+set bmargin at screen 0.25
 set tmargin at screen 0.95
 
 set key bottom right #samplen 0.9
@@ -17,14 +17,14 @@ beta_f = 3.9e-10
 mu_f = 1.787e-3
 #k = 2e-17
 k_min = 1e-19
-k_max = 1e-13
+k_max = 1e-14
 
 # forcing parameters
 #f = 1.0/(3600.0*24.0)
 #f_min = 1.0/(3600.0*24*365)  # 1/s
 #f_max = 1.0/3600.0           # 1/s
-f_min = 1.0/(24*365)          # 1/h
-f_max = 1.0                   # 1/h
+f_min = 0.5*1.0/(24*365)          # 1/h
+f_max = 2.0                   # 1/h
 
 # f: frequency, k: permeability
 #skindepth(f,k) = (k/(phi*mu_f*beta_f*3.141592654*f))**(0.5)  # f: 1/s
@@ -40,9 +40,30 @@ set samples 100
 set isosamples 60
 
 # Generate contours
+load 'BlueSeq.plt'
 set view map
 unset surface
-set contour
+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 xtics (0.0001, 0.01, 1.0)
+set xtics rotate by 90
+set ytics offset 1
+set format x '10^{%T}'
+set format y '10^{%T}'
+
+label_y = 3e-19
+set arrow from 1,k_min to 1,k_max nohead lc black
+set label "hourly" at 1,label_y rotate by 90 offset screen -0.03,0
+
+set arrow from 0.042,k_min to 0.042,k_max nohead lc black
+set label "daily" at 0.042,label_y rotate by 90 offset screen -0.03,0
+
+set arrow from 0.0014,k_min to 0.0014,k_max nohead lc black
+set label "monthly" at 0.0014,label_y rotate by 90 offset screen -0.03,0
+
+set arrow from 0.00011,k_min to 0.00011,k_max nohead lc black
+set label "yearly" at 0.00011,label_y rotate by 90 offset screen +0.02,0
+
 splot skindepth(x,y) title ""