tfig-sediment_flux_daily.gp - hansen-zoet-exp - comparison of lab and model results of sediment advection
(HTM) git clone git://src.adamsgaard.dk/hansen-zoet-exp
(DIR) Log
(DIR) Files
(DIR) Refs
(DIR) Submodules
---
tfig-sediment_flux_daily.gp (1082B)
---
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 2,1
7
8 rswidth = 1
9
10 set lmargin 8.0
11 set bmargin 3.5
12 set rmargin 2.5
13 set tmargin 2.0
14
15 A_f_crit = (rho_s - rho_f) * g * d_s**3.0 / 1000.0 # [kPa]
16
17 labelx = 0.08
18 set label "a" at screen labelx,0.98 font "{:Bold},12"
19 set label "b" at screen labelx,0.48 font "{:Bold},12"
20
21 set yrange [-0.1:1]
22 set ylabel "Max. expected shear depth [m]"
23 set arrow from A_f_crit,-0.1 to A_f_crit,1 nohead lc "black" front
24 plot "max_depth_daily.txt" u ($1/1000):($2) w l t ''
25
26 set xrange [0:*]
27 set yrange [0:*]
28
29 #set xlabel "Normalized horizontal velocity, v_x [-]"
30 set xlabel "Amplitude of water-pressure fluctuations, {/:Italic A}_f [kPa]"
31 #set ylabel "Spec. sediment flux [m²/a]"
32 set ylabel "Sediment flux, {/:Italic q}_t [m²/a]"
33 set key bottom right font ",10" invert #samplen 0.9
34 #set title "{/:Italic v} = 100 m/a"
35 set yrange [0:*]
36 set arrow from A_f_crit,0 to A_f_crit,30 nohead lc "black" front
37 plot "sediment_flux_freq1.157e-5s.txt" u ($1):($2*3600*24*365.25*rswidth) w lp t ""