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 (3184B)
---
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 D_min = 1e-9
15 D_max = 1e-3
16
17 # forcing parameters
18 f_min = 0.5*1.0/(24*365)/36000.0
19 f_max = 2.0/3600.0
20
21 # f: frequency, D: diffusivity
22 skindepth(f,D) = (D/(3.141592654*f))**(0.5) # f: 1/s
23
24 set xlabel "Forcing frequency {/:Italic f} [s^{-1}]"
25 set ylabel "Diffusivity {/:Italic D} [m^2/s]"
26 set cblabel "Skin depth {/:Italic d}_s [m]"
27
28 set logscale xyzcb
29 set xrange [f_min:f_max]
30 set yrange [D_min:D_max]
31 set samples 100
32 set isosamples 60
33
34 # Generate contours
35 set view map
36 unset surface
37 set xtics rotate by 90
38 set ytics offset 1
39 set format x '10^{%T}'
40 set format y '10^{%T}'
41
42 set pm3d
43
44 label_y = 3e-9
45 x=1.0/3600.0
46 set arrow from x,D_min to x,D_max nohead lc "white" front
47 set label "hourly" at x,label_y rotate by 90 front offset screen -0.01,0 textcolor "white"
48 x=1.0/3600.0/24.0
49 set arrow from x,D_min to x,D_max nohead lc "white" front
50 set label "daily" at x,label_y rotate by 90 front offset screen -0.01,0 textcolor "white"
51 x=1.0/3600.0/24.0/30.0
52 set arrow from x,D_min to x,D_max nohead lc "white" front
53 set label "monthly" at x,label_y rotate by 90 front offset screen -0.01,0 textcolor "white"
54 x=1.0/3600.0/24.0/365.0
55 set arrow from x,D_min to x,D_max nohead lc "white" front
56 set label "yearly" at x,label_y rotate by 90 front offset screen -0.01,0 textcolor "white"
57
58 splot skindepth(x,y) notitle with lines palette, \
59 skindepth(x,y) with pm3d notitle
60
61 reset
62 set pm3d
63
64 x=1.0/3600.0
65 set arrow from x,D_min to x,D_max nohead lc "white" front
66 set label "hourly" at x,label_y rotate by 90 front offset screen -0.01,0 textcolor "white"
67 x=1.0/3600.0/24.0
68 set arrow from x,D_min to x,D_max nohead lc "white" front
69 set label "daily" at x,label_y rotate by 90 front offset screen -0.01,0 textcolor "white"
70 x=1.0/3600.0/24.0/30.0
71 set arrow from x,D_min to x,D_max nohead lc "white" front
72 set label "monthly" at x,label_y rotate by 90 front offset screen -0.01,0 textcolor "white"
73 x=1.0/3600.0/24.0/365.0
74 set arrow from x,D_min to x,D_max nohead lc "white" front
75 set label "yearly" at x,label_y rotate by 90 front offset screen -0.01,0 textcolor "white"
76
77 set label "a" at screen 0.01,0.95 font ",12"
78 set label "b" at screen 0.35,0.95 font ",12"
79 set label "c" at screen 0.70,0.95 font ",12"
80
81 set xlabel "Forcing frequency {/:Italic f} [s^{-1}]"
82 set ylabel "Diffusivity {/:Italic D} [m^2/s]"
83 set cblabel "Max. deformation depth [m]"
84
85 set xrange [f_min:f_max]
86 set yrange [D_min:D_max]
87 set xtics rotate by 90
88 set ytics offset 1
89 set format x '10^{%T}'
90 set format y '10^{%T}'
91
92 set view map scale 1
93 set style data pm3d
94 set style function pm3d
95 set xyplane relative 0
96 set nomcbtics
97 set pm3d implicit at b
98 set pm3d scansforward
99 #set pm3d interpolate 10,10
100 set logscale xy
101
102 #unset colorbox
103
104 set title "{/:Italic A}_f = 10 kPa"
105 splot 'max_depth_ampl10e3.txt' notitle with lines palette, \
106 '' with pm3d notitle
107
108 #unset ylabel
109 #set colorbox
110 set title "{/:Italic A}_f = 100 kPa"
111 splot 'max_depth_ampl100e3.txt' notitle with lines palette, \
112 '' with pm3d notitle