commit 6fc849a1b22a0e5876667e698b10d0248836347c
parent 0a7a86a9b0a6dc84d7c32393fd7743437b223850
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date: Wed, 11 Dec 2019 13:26:03 +0100
Use second as time unit in skin-depth figure
Diffstat:
1 file changed, 17 insertions(+), 16 deletions(-)
diff --git a/fig-skin_depth/fig.gp b/fig-skin_depth/fig.gp
@@ -20,17 +20,14 @@ k_min = 1e-19
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 = 0.5*1.0/(24*365) # 1/h
-f_max = 2.0 # 1/h
+f_min = 0.5*1.0/(24*365)/36000.0
+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
-skindepth(f,k) = (k/(phi*mu_f*beta_f*3.141592654*(f/3600)))**(0.5) # f: 1/h
+skindepth(f,k) = (k/(phi*mu_f*beta_f*3.141592654*f))**(0.5) # f: 1/s
+#skindepth(f,k) = (k/(phi*mu_f*beta_f*3.141592654*(f/3600)))**(0.5) # f: 1/h
-set xlabel "Forcing frequency [1/h]"
+set xlabel "Forcing frequency [1/s]"
set ylabel "Permeability [m^2]"
set logscale xyzcb
@@ -54,16 +51,20 @@ 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
+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 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
+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 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
+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 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
+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
splot skindepth(x,y) title ""
