tfig.gp - cngf-pf-exp1 - experiments for first paper with continuum granular mo… | |
git clone git://src.adamsgaard.dk/manus_continuum_granular1_exp | |
Log | |
Files | |
Refs | |
Submodules | |
README | |
LICENSE | |
--- | |
tfig.gp (4357B) | |
--- | |
1 #!/usr/bin/env gnuplot | |
2 | |
3 set terminal pdfcairo color size 17.8 cm, 10.0 cm font ",10" | |
4 | |
5 set colorsequence podo | |
6 set multiplot | |
7 | |
8 file_interval=60*10 | |
9 plot_interval=60*60 | |
10 inc=plot_interval/file_interval | |
11 L_z=8.0 | |
12 set yrange [4.0:0.0] | |
13 set ytics 1.0 | |
14 | |
15 t0 = 48*inc - 1 | |
16 tend = 72*inc - 1 | |
17 | |
18 filename(t) = sprintf('rate/sim.output%05d.txt', t) | |
19 | |
20 set pm3d map | |
21 #set palette rgb 10,6,6 | |
22 set palette | |
23 | |
24 set ylabel "{/:Normal Depth [m]}" offset 0.5 | |
25 | |
26 set origin 0.06, 0.0 | |
27 set size 0.313, 0.46 | |
28 xmin = 0.0 | |
29 xmax = 200.0 | |
30 set xrange [xmin:xmax] | |
31 set arrow from 88,skin_depth to xmax,skin_depth nohead lc black front | |
32 set colorbox user origin 0.001, 0.09 size 0.015, 0.34 | |
33 set cblabel "Time [d]" | |
34 set cbrange [2:3] | |
35 set cbtics 0.2 | |
36 set label "d" font ",11" at graph 0.00,1.1 | |
37 set xlabel "{/:Normal Effective normal stress [kPa]}" | |
38 plot for [t = t0:tend:inc*2] filename(t) u ($3/1e3):(L_z - $1):(t*file_i… | |
39 unset arrow | |
40 unset colorbox | |
41 | |
42 set origin 0.373, 0.0 | |
43 set size 0.313, 0.46 | |
44 unset label | |
45 set label "e (Speed controlled)" font ",11" at graph 0.00,1.1 | |
46 set xlabel "{/:Normal Shear strain rate [1/d]}" font ",10" | |
47 xmin = -2.0 | |
48 xmax = 50.0 | |
49 set xrange [xmin:xmax] | |
50 set arrow from xmin,max_def_depth to 15,max_def_depth nohead front lc bl… | |
51 plot for [t = t0:tend:inc*2] filename(t) u ($6*60*60*24):(L_z - $1):(t*f… | |
52 unset format x | |
53 unset arrow | |
54 | |
55 filename(t) = sprintf('stress/sim.output%05d.txt', t) | |
56 | |
57 set origin 0.686, 0.0 | |
58 set size 0.313, 0.46 | |
59 unset label | |
60 set label "f (Stress controlled)" font ",11" at graph 0.00,1.1 | |
61 set xlabel "{/:Normal Shear strain rate [1/min]}" | |
62 xmin = -1.0 | |
63 xmax = 12.0 | |
64 set xrange [xmin:xmax] | |
65 plot for [t = t0:tend:inc*2] filename(t) u ($6*60):(L_z - $1):(t*file_in… | |
66 unset format x | |
67 | |
68 set origin 0.043, 0.5 | |
69 set size 0.330, 0.46 | |
70 set key top right samplen 1.0 font ",10" | |
71 set xrange [0:7] | |
72 set yrange [0:240] | |
73 set ylabel "{/:Normal [kPa]}" | |
74 set xlabel "{/:Normal Time [d]}" | |
75 set xtics 1 | |
76 set ytics 50 | |
77 unset label | |
78 set label "a" font ",11" at graph 0.00,1.1 | |
79 plot \ | |
80 "rate/timeseries.txt" u ($0/1008*7):($2/1000) w l lw 2 lc "black" t "Eff… | |
81 "" u ($0/1008*7):($3/1000) w l lw 2 lt 3 dt "-" t "Water pressure" | |
82 | |
83 set origin 0.373, 0.5 | |
84 set size 0.313, 0.46 | |
85 unset label | |
86 set label "b (Speed controlled)" font ",11" at graph 0.00,1.1 | |
87 unset y2label | |
88 unset y2tics | |
89 set ytics 1.0 | |
90 set yrange [4.0:0.0] | |
91 set ylabel "{/:Normal Depth of max. strain rate [m]}" | |
92 set arrow from 0,max_def_depth to 7,max_def_depth nohead front lc black … | |
93 plot "rate/max_strainrate_depth.txt" u ($0/1008*7):(8.0 - $1) w l lw 2 t… | |
94 | |
95 set origin 0.686, 0.5 | |
96 set size 0.313, 0.46 | |
97 unset arrow | |
98 unset label | |
99 set label "c (Stress controlled)" font ",11" at graph 0.00,1.1 | |
100 unset y2label | |
101 unset y2tics | |
102 set yrange [0.0:9.8] | |
103 set ytics 2.0 | |
104 set ylabel "{/:Normal Shear speed [km/d]}" | |
105 plot "stress/timeseries.txt" u ($0/1008*7):($1*3600.*24./1000.) w l lw 2… | |
106 | |
107 unset label | |
108 filename(t) = sprintf('rate/sim.output%05d.txt', t) | |
109 set origin 0.498, 0.11 | |
110 set size 0.17, 0.28 | |
111 set xlabel "{/:Normal Shear speed [km/a]}" font ",8" offset 0,0.7 | |
112 set ylabel "Depth [m]" font ",8" offset 0.7,0.0 | |
113 set ytics 1 font ",8" | |
114 set xtics 0.5 font ",8" | |
115 unset arrow | |
116 xmin = -0.1 | |
117 xmax = 1.1 | |
118 set arrow from xmin,max_def_depth to xmax,max_def_depth nohead front lc … | |
119 set xrange [xmin:xmax] | |
120 set yrange [4.0:0.0] | |
121 plot for [t = t0:tend:inc*2] filename(t) u ($2*60*60*24*365/1000):(L_z -… | |
122 | |
123 filename(t) = sprintf('stress/sim.output%05d.txt', t) | |
124 set origin 0.81, 0.11 | |
125 set size 0.17, 0.28 | |
126 set xlabel "{/:Normal Shear speed [km/d]}" font ",8" | |
127 unset arrow | |
128 xmin = -1 | |
129 xmax = 10 | |
130 set xrange [xmin:xmax] | |
131 set xtics 2 | |
132 plot for [t = t0:tend:inc*2] filename(t) u ($2*60*60*24/1000):(L_z - $1)… | |
133 | |
134 filename(t) = sprintf('stress/sim.output%05d.txt', t) | |
135 set origin 0.095, 0.09 | |
136 set size 0.14, 0.23 | |
137 set xlabel "{/:Normal Water pressure [kPa]}" font ",8" offset 0,1.3 | |
138 set ylabel "" | |
139 unset arrow | |
140 xmin = 0 | |
141 xmax = 200.0 | |
142 set arrow from xmin,skin_depth to xmax,skin_depth nohead lc black front | |
143 set xrange [xmin:xmax] | |
144 set xtics 100 offset 0,0.5 | |
145 set ytics 2 offset 0.5 | |
146 plot for [t = t0:tend:inc*2] filename(t) u ($4/1e3):(L_z - $1):(t*file_i… | |
147 | |
148 unset multiplot |