tfig-sediment_flux.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.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