fig.gp (900B)
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 layout 2,1 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('sim.output%05d.txt', t) 19 #set origin 0.095, 0.09 20 #set size 0.14, 0.23 21 set xlabel "{/:Normal Effective normal stress [kPa]}" 22 set ylabel "" 23 unset arrow 24 xmin = 0 25 xmax = 200.0 26 set arrow from xmin,skin_depth to xmax,skin_depth nohead lc black front 27 set xrange [xmin:xmax] 28 set xtics 100 offset 0,0.5 29 set ytics 2 offset 0.5 30 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 1 t "" 31 32 reset 33 set xlabel "{/:Normal Time [d]}" 34 set ylabel "{/:Normal Shear velocity [m/s]}" 35 plot "timeseries.txt" u ($0*file_interval/(3600*24)):($1) w l lw 1 t ""