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 (2539B)
---
1 #!/usr/bin/env gnuplot
2
3 #set terminal pdfcairo color size 7.5 cm, 7.5 cm
4 #set terminal pdfcairo color size 13.35 cm, 4.45 cm font ",10"
5 #set terminal pdfcairo color size 10.68 cm, 4.45 cm font ",10"
6 #set terminal pdfcairo color size 11.87 cm, 4.45 cm font ",10"
7 set terminal pdfcairo color size 11.87 cm, 4.45 cm font ",10"
8 set multiplot layout 1,4 \
9 margins 0.05,0.98,0.23,0.90 \
10 spacing 0.02,0.02
11
12 set key bottom right #samplen 0.9
13
14 file_interval=60*10
15 plot_interval=60*60
16 inc=plot_interval/file_interval
17 #set yrange [4.0:8.0]
18 L_z=8.0
19 set yrange [4.0:0.0]
20
21 t0 = 47*inc
22 tend = 71*inc
23
24 filename(t) = sprintf('sim.output%05d.txt', t)
25
26 set pm3d map
27 #set palette rgb 10,6,6
28 set palette
29
30 unset colorbox
31 set ylabel "{/:Normal Depth [m]}" offset 1.8
32 xmin = 0.0
33 xmax = 200.0
34 set xrange [xmin:xmax]
35 set arrow from xmin,skin_depth to xmax,skin_depth nohead lc black front
36 set arrow from xmin,max_def_depth to xmax,max_def_depth nohead front lc black dashtype 2
37 set label "A" at screen 0.090,0.97 font ",12"
38 set xlabel "{/:Normal Water pressure [kPa]}"
39 plot for [t = t0:tend:inc] filename(t) u ($4/1e3):(L_z - $1):(t-t0) w l lc palette lw 1 t ""
40
41 unset colorbox
42 set ylabel ""
43 unset ytics
44 xmin = 0.0
45 xmax = 200.0
46 set xrange [xmin:xmax]
47 set arrow from xmin,skin_depth to xmax,skin_depth nohead lc black front
48 set arrow from xmin,max_def_depth to xmax,max_def_depth nohead front lc black dashtype 2
49 set label "B" at screen 0.220,0.97 font ",12"
50 set xlabel "{/:Normal Eff. normal stress [kPa]}"
51 plot for [t = t0:tend:inc] filename(t) u ($3/1e3):(L_z - $1):(t-t0) w l lc palette lw 1 t ""
52
53 unset label
54 set label "C (speed controlled)" at screen 0.395,0.97 font ",12"
55 set xlabel "{/:Normal Shear speed [km/a]}"
56 set ylabel ""
57 unset ytics
58 unset arrow
59 xmin = -0.1
60 xmax = 1.1
61 set xrange [xmin:xmax]
62 set arrow from xmin,skin_depth to xmax,skin_depth nohead lc black front
63 set arrow from xmin,max_def_depth to xmax,max_def_depth nohead front lc black dashtype 2
64 plot for [t = t0:tend:inc] filename(t) u ($2*60*60*24*365/1000):(L_z - $1):(t-t0) w l lc palette lw 1 t ""
65
66 unset label
67 set label "D (speed controlled)" at screen 0.705,0.97 font ",12"
68 set xlabel "{/:Normal Shear strain rate [1/d]}"
69 set ylabel ""
70 xmin = -2.0
71 xmax = 22.0
72 set xrange [xmin:xmax]
73 set arrow from xmin,skin_depth to xmax,skin_depth nohead lc black front
74 set arrow from xmin,max_def_depth to xmax,max_def_depth nohead front lc black dashtype 2
75 unset ytics
76 plot for [t = t0:tend:inc] filename(t) u ($6*60*60*24):(L_z - $1):(t-t0) w l lc palette lw 1 t ""
77 unset format x
78
79 unset multiplot