tfig-strain_distribution.gp - hansen-zoet-exp - comparison of lab and model results of sediment advection
HTML git clone git://src.adamsgaard.dk/hansen-zoet-exp
DIR Log
DIR Files
DIR Refs
DIR Submodules
---
tfig-strain_distribution.gp (2185B)
---
1 #!/usr/bin/env gnuplot
2
3 reset
4
5 set terminal pdfcairo enhanced color size 8.3 cm, 8.3 cm font ",10"
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.02 : 0]
15
16 labely = 0.98
17 set label "c" at screen 0.08,labely font "{:Bold},12"
18 set label "d" at screen 0.58,labely font "{:Bold},12"
19
20 #set xlabel "Normalized horizontal velocity, {/:Italic v}_x [-]"
21 set xlabel "Horizontal velocity, {/:Italic v} [m/a]"
22 set ylabel "Depth beneath ice-bed interface, {/:Italic z} [m]"
23 set key bottom right font ",10" #samplen 0.9
24 set title "{/:Italic v} = 100 m/a"
25 plot "strain_distribution_N8kPa_v3.169e-6ms.txt" u ($2*365*24*3600):(-1*$1+0.11) w l lw 2 title "{/:Italic N}' = 8 kPa", \
26 "strain_distribution_N22kPa_v3.169e-6ms.txt" u ($2*365*24*3600):(-1*$1+0.11) w l lw 2 title "{/:Italic N}' = 22 kPa", \
27 "strain_distribution_N50kPa_v3.169e-6ms.txt" u ($2*365*24*3600):(-1*$1+0.11) w l lw 2 title "{/:Italic N}' = 50 kPa"
28
29 #unset xrange
30
31 ##set xlabel "Normalized horizontal velocity, v_x [-]"
32 #set xlabel "Shear-strain rate, dv_x/dz [1/d]"
33 #set ylabel "Vertical position, z [m]"
34 #set key bottom right font ",10" #samplen 0.9
35 #plot "strain_distribution_N8kPa_v9.506e-6ms.txt" u ($6*24*3600):($1-0.11) w l lw 2 title "{/Symbol s}_n' = 8 kPa", \
36 #"strain_distribution_N22kPa_v9.506e-6ms.txt" u ($6*24*3600):($1-0.11) w l lw 2 title "{/Symbol s}_n' = 22 kPa", \
37 #"strain_distribution_N50kPa_v9.506e-6ms.txt" u ($6*24*3600):($1-0.11) w l lw 2 title "{/Symbol s}_n' = 50 kPa"
38
39
40 set xrange [-15 : 315]
41 #set yrange [0 : 0.11]
42
43 #set xlabel "Normalized horizontal velocity, {/:Italic v}_x [-]"
44 #set xlabel "Depth, {/:Italic v} [m/a]"
45 unset ylabel
46 set xtics 0,100,300
47 set key bottom right font ",10" #samplen 0.9
48 set title "{/:Italic v} = 300 m/a"
49 plot "strain_distribution_N8kPa_v9.506e-6ms.txt" u ($2*365*24*3600):(-1*$1+0.11) w l lw 2 title "{/:Italic N}' = 8 kPa", \
50 "strain_distribution_N22kPa_v9.506e-6ms.txt" u ($2*365*24*3600):(-1*$1+0.11) w l lw 2 title "{/:Italic N}' = 22 kPa", \
51 "strain_distribution_N50kPa_v9.506e-6ms.txt" u ($2*365*24*3600):(-1*$1+0.11) w l lw 2 title "{/:Italic N}' = 50 kPa"