fig-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 ""