cngf-pf-exp1

commit 50fc4c759639e4093065e6d5ebf9a8b4e49d46ef
parent c59a6698c1923b36468fe67db433aba036bb98a1
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Wed, 26 Jun 2019 14:48:18 +0200

Tweak figure appreance

Diffstat:
A
fig5/Makefile
 | 
11
+++++++++++
A
fig5/fig.gp
 | 
54
++++++++++++++++++++++++++++++++++++++++++++++++++++++

2 files changed, 65 insertions(+), 0 deletions(-)
diff --git a/fig5/Makefile b/fig5/Makefile
@@ -0,0 +1,11 @@
+FIG = fig5
+
+default: ../$(FIG).pdf
+
+../$(FIG).pdf: fig.gp
+	gnuplot $< > $@
+
+clean:
+	$(RM) ../$(FIG).pdf
+
+.PHONY: default clean
diff --git a/fig5/fig.gp b/fig5/fig.gp
@@ -0,0 +1,54 @@
+#!/usr/bin/env gnuplot
+
+reset
+
+set terminal pdfcairo color size 7.5 cm, 7.5 cm
+#set terminal pdfcairo mono size 7.5 cm, 7.5 cm
+
+set lmargin at screen 0.24
+set rmargin at screen 0.75
+set bmargin at screen 0.17
+set tmargin at screen 0.95
+
+set key bottom right #samplen 0.9
+
+# hydraulic parameters
+phi = 0.25
+beta_f = 4.5e-10
+mu_f = 0.001
+#k = 2e-17
+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 = 1.0/(24*365)          # 1/h
+f_max = 1.0                   # 1/h
+
+# 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
+
+set xlabel "Forcing frequency [1/h]"
+set ylabel "Permeability [m^2]"
+
+set logscale xyzcb
+set xrange [f_min:f_max]
+set yrange [k_min:k_max]
+set samples 100
+set isosamples 60
+#set table 'skindepth.dat'
+#splot skindepth(x,y)
+#unset table
+
+# Generate contours
+set view map
+unset surface
+set contour
+set key outside title "Skin depth\n[m]" width 1.0 samplen 2.5
+#set key outside bottom center vertical
+#set key outside top center vertical
+#set size ratio 0.35
+splot skindepth(x,y) title ""