fig.gp (1162B)
1 #!/usr/bin/env gnuplot 2 reset 3 4 set terminal pdfcairo color size 15 cm, 10 cm 5 6 set multiplot layout 4,1 \ 7 margins 0.15,0.88,0.13,0.95 \ 8 spacing 0.0,0.03 9 10 file_interval_seconds=600.0 11 conv=1.0/(24.0*3600.0/file_interval_seconds) * 140.0/138.5 12 13 set key top right samplen 1.0 14 unset xlabel 15 set xtics format "" 16 set xrange [0:140] 17 set yrange [0:380] 18 set ylabel "[kPa]" 19 set ytics 100 20 21 plot \ 22 "timeseries.txt" u ($0*conv):($2/1000) w l lw 2 lc "black" t "effective normal stress", \ 23 "" u ($0*conv):($3/1000) w l lw 2 lt 3 dt "-" t "water pressure" 24 25 unset yrange 26 set yrange [5e-20:1e-15] 27 set ytics (1e-19,1e-18,1e-17,1e-16) 28 set ylabel "Permeability [m^2]" 29 set logscale y 30 plot "max_strainrate_depth.txt" u ($0*conv):3 w l lw 2 lt 5 t "" 31 unset logscale y 32 33 unset ylabel 34 unset ytics 35 set yrange [0:100] 36 set y2range [0:100] 37 set y2tics 20 38 set y2label "Shear stress [kPa]" 39 plot "timeseries.txt" u ($0*conv):($4*$2/1000) w l lw 2 lt 1 t "" 40 41 set xtics format "%g" 42 unset y2label 43 unset y2tics 44 set ytics 0.5 45 set yrange [3.0:0.0] 46 set xlabel "Time [d]" 47 set ylabel "Depth of max. strain rate [m]" 48 plot "max_strainrate_depth.txt" u ($0*conv):(8.0 - $1) w l lw 2 lt 4 t "" 49 50 unset multiplot