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 (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