fig.gp (2077B)
1 #!/usr/bin/env gnuplot 2 3 #set terminal pdfcairo color size 7.5 cm, 7.5 cm 4 set terminal pdfcairo color size 15.0 cm, 5.0 cm 5 set multiplot layout 1,4 \ 6 margins 0.09,0.97,0.23,0.95 \ 7 spacing 0.03,0.03 8 9 set yrange [4.0:8.0] 10 11 set key bottom right #samplen 0.9 12 13 file_interval=60*10 14 plot_interval=60*60 15 inc=plot_interval/file_interval 16 17 #t0 = 47*inc 18 #tend = 71*inc 19 t0 = 20*inc 20 tend = 47*inc 21 22 filename(t) = sprintf('sim.output%05d.txt', t) 23 24 set pm3d map 25 set palette rgb 10,6,6 26 27 #load '../skindepth.gp' 28 #k = 2e-17 29 #f = 1.0/(3600.0*24.0) 30 #f(x) = 8.0 - skindepth(f,k) 31 32 set label "{/:Bold a}" at graph 0.90,0.93 33 unset colorbox 34 set xlabel "{/:Normal Water pressure [kPa]}" 35 set ylabel "{/:Normal Vertical position [m]}" 36 set xrange [0:200] 37 set arrow from 0,skin_depth to 200,skin_depth nohead lc black front 38 set arrow from 0,max_def_depth to 200,max_def_depth nohead front lc black dashtype 2 39 plot for [t = t0:tend:inc] filename(t) u ($4/1e3):1:(t-t0) w l lc palette lw 1 t "" 40 41 unset label 42 set label "{/:Bold b}" at graph 0.90,0.93 43 set xlabel "{/:Normal Eff. normal stress [kPa]}" 44 set ylabel "" 45 set xrange [0:200] 46 unset ytics 47 plot for [t = t0:tend:inc] filename(t) u ($3/1e3):1:(t-t0) w l lc palette lw 1 t "" 48 49 unset label 50 set label "{/:Bold c}" at graph 0.90,0.93 51 set xlabel "{/:Normal Shear velocity [km/a]}" 52 set ylabel "" 53 set xrange [0.0:1.0] 54 #unset xrange 55 #set xtics (0,0.02,0.04,0.06) 56 unset ytics 57 unset arrow 58 set arrow from 0,skin_depth to 1,skin_depth nohead lc black front 59 set arrow from 0,max_def_depth to 1,max_def_depth nohead front lc black dashtype 2 60 plot for [t = t0:tend:inc] filename(t) u ($2*60*60*24*365/1000):1:(t-t0) w l lc palette lw 1 t "" 61 62 unset label 63 set label "{/:Bold d}" at graph 0.90,0.93 64 set xlabel "{/:Normal Shear strain rate [1/d]}" 65 set ylabel "" 66 set xrange [0:20] 67 unset ytics 68 unset arrow 69 set arrow from 0,skin_depth to 20,skin_depth nohead lc black front 70 set arrow from 0,max_def_depth to 20,max_def_depth nohead front lc black dashtype 2 71 plot for [t = t0:tend:inc] filename(t) u ($6*60*60*24):1:(t-t0) w l lc palette lw 1 t "" 72 unset format x 73 74 unset multiplot