fig-sediment_flux.gp (1331B)
1 #!/usr/bin/env gnuplot 2 3 reset 4 5 set terminal pdfcairo enhanced color size 7.5 cm, 10 cm 6 set multiplot layout 2,1 7 8 rswidth = 1 9 10 #set lmargin 7.0 11 #set bmargin 3.5 12 #set rmargin 2.0 13 #set tmargin 2.0 14 15 #set xrange [-5 : 105] 16 #set yrange [0 : 0.11] 17 18 # https://sodocumentation.net/gnuplot/topic/8825/fit-data-with-gnuplot 19 #f(x) = a * exp(b * x) 20 f(x, v) = a * x**b * v 21 fit f(x, 100) "sediment_flux_v3.169e-6ms.txt" u ($1/1000):($2*3600*24*365.25*rswidth) via a,b 22 fitparams = sprintf("{/:Italic q}_t = %.2e {/:Italic v N'}^{%.3g}", a, b) 23 24 set xrange [0:*] 25 set yrange [0:*] 26 27 #set xlabel "Normalized horizontal velocity, v_x [-]" 28 set xlabel "Effective stress, {/:Italic N'} [kPa]" 29 #set ylabel "Spec. sediment flux [m²/a]" 30 set ylabel "Sediment flux, {/:Italic q}_t [m³/a]" 31 set key bottom right font ",10" #samplen 0.9 32 set title "{/:Italic v} = 100 m/a" 33 set yrange [0:1] 34 plot "rs_sediment_flux_100ma.txt" u ($1/1000):($2) w p ps 2 t "Hansen and Zoet 2022", \ 35 "sediment_flux_v3.169e-6ms.txt" u ($1/1000):($2*3600*24*365.25*rswidth) w p ps 0.5 t "CNGF-PF", \ 36 f(x, 100) t fitparams 37 38 set title "{/:Italic v} = 300 m/a" 39 set yrange [0:*] 40 plot "rs_sediment_flux_300ma.txt" u ($1/1000):($2) w p ps 2 t "Hansen and Zoet 2022", \ 41 "sediment_flux_v9.506e-6ms.txt" u ($1/1000):($2*3600*24*365.25*rswidth) w p ps 0.5 t "CNGF-PF", \ 42 f(x, 300) t fitparams