hansen-zoet-exp

comparison of lab and model results of sediment advection
git clone git://src.adamsgaard.dk/hansen-zoet-exp # fast
git clone https://src.adamsgaard.dk/hansen-zoet-exp.git # slow
Log | Files | Refs | Submodules | LICENSE Back to index

fig-strain_distribution.gp (1898B)

      1 #!/usr/bin/env gnuplot
      2 
      3 reset
      4 
      5 set terminal pdfcairo enhanced color size 10 cm, 10 cm
      6 set multiplot layout 1,2
      7 
      8 set lmargin 10.0
      9 set bmargin 3.5
     10 set rmargin 2.0
     11 set tmargin 2.0
     12 
     13 set xrange [-5 : 105]
     14 set yrange [0.09 : 0.11]
     15 
     16 #set xlabel "Normalized horizontal velocity, v_x [-]"
     17 set xlabel "Horizontal velocity, v_x [m/a]"
     18 set ylabel "Vertical position, z [m]"
     19 set key bottom right font ",10" #samplen 0.9
     20 set title "v = 100 m/a"
     21 plot "strain_distribution_N8kPa_v3.169e-6ms.txt" u ($2*365*24*3600):1 w l lw 2 title "{/Symbol s}_n' = 8 kPa", \
     22      "strain_distribution_N22kPa_v3.169e-6ms.txt" u ($2*365*24*3600):1 w l lw 2 title "{/Symbol s}_n' = 22 kPa", \
     23      "strain_distribution_N50kPa_v3.169e-6ms.txt" u ($2*365*24*3600):1 w l lw 2 title "{/Symbol s}_n' = 50 kPa"
     24 
     25 #unset xrange
     26 
     27 ##set xlabel "Normalized horizontal velocity, v_x [-]"
     28 #set xlabel "Shear-strain rate, dv_x/dz [1/d]"
     29 #set ylabel "Vertical position, z [m]"
     30 #set key bottom right font ",10" #samplen 0.9
     31 #plot "strain_distribution_N8kPa_v9.506e-6ms.txt" u ($6*24*3600):1 w l lw 2 title "{/Symbol s}_n' = 8 kPa", \
     32      #"strain_distribution_N22kPa_v9.506e-6ms.txt" u ($6*24*3600):1 w l lw 2 title "{/Symbol s}_n' = 22 kPa", \
     33      #"strain_distribution_N50kPa_v9.506e-6ms.txt" u ($6*24*3600):1 w l lw 2 title "{/Symbol s}_n' = 50 kPa"
     34 
     35 
     36 set xrange [-5 : 305]
     37 #set yrange [0 : 0.11]
     38 
     39 #set xlabel "Normalized horizontal velocity, v_x [-]"
     40 set xlabel "Horizontal velocity, v_x [m/a]"
     41 set ylabel "Vertical position, z [m]"
     42 set key bottom right font ",10" #samplen 0.9
     43 set title "v = 300 m/a"
     44 plot "strain_distribution_N8kPa_v9.506e-6ms.txt" u ($2*365*24*3600):1 w l lw 2 title "{/Symbol s}_n' = 8 kPa", \
     45      "strain_distribution_N22kPa_v9.506e-6ms.txt" u ($2*365*24*3600):1 w l lw 2 title "{/Symbol s}_n' = 22 kPa", \
     46      "strain_distribution_N50kPa_v9.506e-6ms.txt" u ($2*365*24*3600):1 w l lw 2 title "{/Symbol s}_n' = 50 kPa"