fig.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_interval/(24.*3600.)) w l lc palette lw 2 t "" 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 black dashtype 2 51 plot for [t = t0:tend:inc*2] filename(t) u ($6*60*60*24):(L_z - $1):(t*file_interval/(24.*3600.)) w l lc palette lw 2 t "" 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_interval/(24.*3600.)) w l lc palette lw 2 t "" 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 "Effective normal stress", \ 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 dashtype 2 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 lc "red" t "" 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 black dashtype 2 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 - $1):(t*file_interval/(24.*3600.)) w l lc palette lw 1 t "" 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):(t*file_interval/(24.*3600.)) w l lc palette lw 1 t "" 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_interval/(24.*3600.)) w l lc palette lw 1 t "" 147 148 unset multiplot