tfig.gp - cngf-pf - continuum model for granular flows with pore-pressure dynamics (renamed from 1d_fd_simple_shear)
 (HTM) git clone git://src.adamsgaard.dk/cngf-pf
 (DIR) Log
 (DIR) Files
 (DIR) Refs
 (DIR) README
 (DIR) LICENSE
       ---
       tfig.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 ""