ttweak parameters and plot - 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
 (DIR) LICENSE
       ---
 (DIR) commit 966ef7d8f2abe662694a68a21d709da2737d9755
 (DIR) parent cb02243d0ed498165972a51aeeea12deb11280b8
 (HTM) Author: Anders Damsgaard <anders@adamsgaard.dk>
       Date:   Fri, 22 Apr 2022 13:42:00 +0200
       
       ttweak parameters and plot
       
       Diffstat:
         M fig-effective_stress/Makefile       |       5 +++--
         M fig-effective_stress/fig-sediment_… |       7 ++++---
         M fig-effective_stress/rs_sediment_f… |       6 +++---
       
       3 files changed, 10 insertions(+), 8 deletions(-)
       ---
 (DIR) diff --git a/fig-effective_stress/Makefile b/fig-effective_stress/Makefile
       t@@ -1,7 +1,7 @@
        BIN = ../cngf-pf/cngf-pf
        FIG = $(shell basename $(CURDIR))
        
       -stressvals = 8 16 22 32 41 50
       +stressvals = 5 8.53 10 15 20 22.2 25 30 35 40 45 50 50.4
        
        default: \
                ../$(FIG)-strain_distribution.pdf \
       t@@ -19,9 +19,10 @@ strain_distribution_N50kPa.txt: $(BIN) Makefile
                        -L 0.11 \
                        -n $${P}e3 \
                        -A 0.48 \
       -                -b 0.001 \
       +                -b 0.0003 \
                        > strain_distribution_N$${P}kPa.txt; \
                done'
       +        # also close: lower A, higher b
        
        mohr_coulomb.txt: strain_distribution_N50kPa.txt
                sh -c '\
 (DIR) diff --git a/fig-effective_stress/fig-sediment_flux.gp b/fig-effective_stress/fig-sediment_flux.gp
       t@@ -15,6 +15,7 @@ rswidth = 1
        #set xrange [-5 : 105]
        #set yrange [0 : 0.11]
        
       +set xrange [0:*]
        set yrange [0:*]
        
        #set xlabel "Normalized horizontal velocity, v_x [-]"
       t@@ -22,5 +23,5 @@ set xlabel "Effective stress, [kPa]"
        #set ylabel "Spec. sediment flux [m²/a]"
        set ylabel "Sediment flux [m³/a]"
        set key bottom right font ",10" #samplen 0.9
       -plot "sediment_flux.txt" u ($1/1000):($2*3600*24*365.25*rswidth) w lp lw 2 t "CNGF-PF", \
       -        "rs_sediment_flux_100ma.txt" u ($1/1000):($2) w p t "Hansen and Zoet 2022"
       -\ No newline at end of file
       +plot "sediment_flux.txt" u ($1/1000):($2*3600*24*365.25*rswidth) w l lw 2 t "CNGF-PF", \
       +        "rs_sediment_flux_100ma.txt" u ($1/1000):($2) w errorbars t "Hansen and Zoet 2022"
       +\ No newline at end of file
 (DIR) diff --git a/fig-effective_stress/rs_sediment_flux_100ma.txt b/fig-effective_stress/rs_sediment_flux_100ma.txt
       t@@ -1,3 +1,3 @@
       -8.53e3        0.351
       -22.2e3        0.510
       -50.4e3        0.661
       +8.53e3        0.351        0.0238
       +22.2e3        0.510        0.0434
       +50.4e3        0.661        0.0434