fig.gp (1990B)
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 5.93 cm, 4.45 cm font ",10" 6 set multiplot layout 1,2 \ 7 margins 0.03,0.96,0.23,0.90 \ 8 spacing 0.04,0.04 9 10 set key bottom right #samplen 0.9 11 12 file_interval=60*10 13 plot_interval=60*60 14 inc=plot_interval/file_interval 15 #set yrange [4.0:8.0] 16 L_z=8.0 17 set yrange [4.0:0.0] 18 19 t0 = 47*inc 20 tend = 71*inc 21 22 filename(t) = sprintf('sim.output%05d.txt', t) 23 24 set pm3d map 25 #set palette rgb 10,6,6 26 set palette 27 28 unset colorbox 29 #set ylabel "{/:Normal Depth [m]}" offset 1.0 30 #xmin = 0.0 31 #xmax = 200.0 32 #set xrange [xmin:xmax] 33 #set arrow from xmin,skin_depth to xmax,skin_depth nohead lc black front 34 #set arrow from xmin,max_def_depth to xmax,max_def_depth nohead front lc black dashtype 2 35 #set label "A" at screen 0.090,0.97 font ",12" 36 #set xlabel "{/:Normal Eff. normal stress [kPa]}" 37 #plot for [t = t0:tend:inc] filename(t) u ($3/1e3):(L_z - $1):(t-t0) w l lc palette lw 1 t "" 38 39 unset label 40 set label "E (stress controlled)" at screen 0.020,0.97 font ",12" 41 set xlabel "{/:Normal Shear speed [km/d]}" 42 set ylabel "" 43 unset ytics 44 unset arrow 45 xmin = -1 46 xmax = 10 47 set xrange [xmin:xmax] 48 set arrow from xmin,skin_depth to xmax,skin_depth nohead lc black front 49 set arrow from xmin,max_def_depth to xmax,max_def_depth nohead front lc black dashtype 2 50 plot for [t = t0:tend:inc] filename(t) u ($2*60*60*24/1000):(L_z - $1):(t-t0) w l lc palette lw 1 t "" 51 52 unset label 53 unset arrow 54 set label "F (stress controlled)" at screen 0.510,0.97 font ",12" 55 set xlabel "{/:Normal Shear strain rate [1/min]}" 56 set ylabel "" 57 xmin = -1.0 58 xmax = 12.0 59 set xrange [xmin:xmax] 60 set arrow from xmin,skin_depth to xmax,skin_depth nohead lc black front 61 set arrow from xmin,max_def_depth to xmax,max_def_depth nohead front lc black dashtype 2 62 unset ytics 63 plot for [t = t0:tend:inc] filename(t) u ($6*60):(L_z - $1):(t-t0) w l lc palette lw 1 t "" 64 unset format x 65 66 unset multiplot