fig.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