tfig.gp - cngf-pf-exp1 - experiments for first paper with continuum granular model
 (HTM) git clone git://src.adamsgaard.dk/manus_continuum_granular1_exp
 (DIR) Log
 (DIR) Files
 (DIR) Refs
 (DIR) Submodules
 (DIR) README
 (DIR) LICENSE
       ---
       tfig.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