tfig.gp - cngf-pf-exp1 - experiments for first paper with continuum granular model
HTML git clone git://src.adamsgaard.dk/manus_continuum_granular1_exp
DIR Log
DIR Files
DIR Refs
DIR Submodules
DIR README
DIR LICENSE
---
tfig.gp (3280B)
---
1 #!/usr/bin/env gnuplot
2
3 reset
4
5 set terminal pdfcairo color size 17.8 cm, 4.45 cm font ",10"
6
7 set multiplot layout 1,3 \
8 margin 0.07,0.92,0.25,0.86 \
9 spacing 0.17,0.0
10
11 set key bottom right #samplen 0.9
12
13 # hydraulic parameters
14 phi = 0.25
15 alpha = 1e-8
16 beta_f = 3.9e-10
17 mu_f = 1.787e-3
18 #k = 2e-17
19 k_min = 1e-19
20 k_max = 1e-14
21
22 # forcing parameters
23 f_min = 0.5*1.0/(24*365)/36000.0
24 f_max = 2.0/3600.0
25
26 # f: frequency, k: permeability
27 skindepth(f,k) = (k/((alpha + phi*beta_f)*mu_f*3.141592654*f))**(0.5) # f: 1/s
28
29 set xlabel "Forcing frequency {/:Italic f} [s^{-1}]"
30 set ylabel "Permeability {/:Italic k} [m^2]"
31 set cblabel "Skin depth {/:Italic d}_s [m]"
32
33 set logscale xyzcb
34 set xrange [f_min:f_max]
35 set yrange [k_min:k_max]
36 set samples 100
37 set isosamples 60
38
39 # Generate contours
40 set view map
41 unset surface
42 set xtics rotate by 90
43 set ytics offset 1
44 set format x '10^{%T}'
45 set format y '10^{%T}'
46
47 set pm3d
48
49 label_y = 3e-17
50 x=1.0/3600.0
51 set arrow from x,k_min to x,k_max nohead lc "white" front
52 set label "hourly" at x,label_y rotate by 90 front offset screen -0.01,0 textcolor "white"
53 x=1.0/3600.0/24.0
54 set arrow from x,k_min to x,k_max nohead lc "white" front
55 set label "daily" at x,label_y rotate by 90 front offset screen -0.01,0 textcolor "white"
56 x=1.0/3600.0/24.0/30.0
57 set arrow from x,k_min to x,k_max nohead lc "white" front
58 set label "monthly" at x,label_y rotate by 90 front offset screen -0.01,0 textcolor "white"
59 x=1.0/3600.0/24.0/365.0
60 set arrow from x,k_min to x,k_max nohead lc "white" front
61 set label "yearly" at x,label_y rotate by 90 front offset screen -0.01,0 textcolor "white"
62
63 splot skindepth(x,y) notitle with lines palette, \
64 skindepth(x,y) with pm3d notitle
65
66 reset
67 set pm3d
68
69 x=1.0/3600.0
70 set arrow from x,k_min to x,k_max nohead lc "white" front
71 set label "hourly" at x,label_y rotate by 90 front offset screen -0.01,0 textcolor "white"
72 x=1.0/3600.0/24.0
73 set arrow from x,k_min to x,k_max nohead lc "white" front
74 set label "daily" at x,label_y rotate by 90 front offset screen -0.01,0 textcolor "white"
75 x=1.0/3600.0/24.0/30.0
76 set arrow from x,k_min to x,k_max nohead lc "white" front
77 set label "monthly" at x,label_y rotate by 90 front offset screen -0.01,0 textcolor "white"
78 x=1.0/3600.0/24.0/365.0
79 set arrow from x,k_min to x,k_max nohead lc "white" front
80 set label "yearly" at x,label_y rotate by 90 front offset screen -0.01,0 textcolor "white"
81
82 set label "a" at screen 0.01,0.95 font ",12"
83 set label "b" at screen 0.35,0.95 font ",12"
84 set label "c" at screen 0.70,0.95 font ",12"
85
86 set xlabel "Forcing frequency {/:Italic f} [s^{-1}]"
87 set ylabel "Permeability {/:Italic k} [m^2]"
88 set cblabel "Max. deformation depth [m]"
89
90 set xrange [f_min:f_max]
91 set yrange [k_min:k_max]
92 set xtics rotate by 90
93 set ytics offset 1
94 set format x '10^{%T}'
95 set format y '10^{%T}'
96
97 set view map scale 1
98 set style data pm3d
99 set style function pm3d
100 set xyplane relative 0
101 set nomcbtics
102 set pm3d implicit at b
103 set pm3d scansforward
104 #set pm3d interpolate 10,10
105 set logscale xy
106
107 #unset colorbox
108
109 set title "{/:Italic A}_f = 10 kPa"
110 splot 'max_depth_ampl10e3.txt' notitle with lines palette, \
111 '' with pm3d notitle
112
113 #unset ylabel
114 #set colorbox
115 set title "{/:Italic A}_f = 100 kPa"
116 splot 'max_depth_ampl100e3.txt' notitle with lines palette, \
117 '' with pm3d notitle