tadd acme.dump - 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) commit 59c213bfb870e3aade3dac55342af9a1b9323236
(DIR) parent 377476ee3e442324ad80925a0828b22025b041a5
(HTM) Author: Anders Damsgaard <anders@adamsgaard.dk>
Date: Thu, 2 Jun 2022 14:26:53 +0200
add acme.dump
Diffstat:
A acme.dump | 119 +++++++++++++++++++++++++++++++
1 file changed, 119 insertions(+), 0 deletions(-)
---
(DIR) diff --git a/acme.dump b/acme.dump
t@@ -0,0 +1,119 @@
+/home/ad/code/hansen-zoet-exp
+/lib/font/bit/lucsans/euro.8.font
+/lib/font/bit/lucm/unicode.9.font
+ 0.0000000 59.9397590
+w Newcol Kill Putall Dump Exit Dump /home/ad/code/hansen-zoet-exp/acme.dump
+c 0 New Cut Paste Snarf Sort Zerox Delcol
+c 1 New Cut Paste Snarf Sort Zerox Delcol
+f 0 2 262 262 1.2142857
+ 2 61 1905 0 0 /home/ad/code/hansen-zoet-exp/Makefile Del Snarf | Look m mc
+f 1 1 0 0 1.2142857
+ 1 52 2905 1 0 /home/ad/code/hansen-zoet-exp/ Del Snarf Get | Look
+F 1 1 4932 4932 50.7142857 4932
+ 3 55 4932 0 1 /home/ad/code/hansen-zoet-exp/+Errors Del Snarf | Look
+sh -c '\
+for v in 3.169e-6 9.506e-6; do \
+ for P in 1 1.2742749857031337 1.6237767391887217 2.0691380811147897 2.6366508987303581 3.3598182862837818 4.2813323987193934 5.4555947811685188 6.9519279617756053 8.8586679041008249 11.28837891684689 14.384498882876629 18.329807108324356 23.357214690901213 29.763514416313178 37.926901907322495 48.329302385717519 61.584821106602604 78.475997035146108 100 8 22 50; do \
+ ./cngf-pf/cngf-pf \
+ -d 530.7e-6 \
+ -m 0.6 \
+ -c 0e3 \
+ -s ${v} \
+ -L 0.11 \
+ -n ${P}e3 \
+ -A 0.33 \
+ -b 0.002 \
+ > strain_distribution_N${P}kPa_v${v}ms.txt; \
+done; done' # also close fit with lower A and higher b
+gnuplot fig-strain_distribution.gp > fig-strain_distribution.pdf
+sh -c '\
+rm -f mohr_coulomb_v*ms.txt; \
+for v in 3.169e-6 9.506e-6; do \
+ for P in 1 1.2742749857031337 1.6237767391887217 2.0691380811147897 2.6366508987303581 3.3598182862837818 4.2813323987193934 5.4555947811685188 6.9519279617756053 8.8586679041008249 11.28837891684689 14.384498882876629 18.329807108324356 23.357214690901213 29.763514416313178 37.926901907322495 48.329302385717519 61.584821106602604 78.475997035146108 100 8 22 50; do \
+ tail -n 1 strain_distribution_N${P}kPa_v${v}ms.txt >> mohr_coulomb_v${v}ms.txt; \
+done; done'
+gnuplot fig-mohr_coulomb.gp > fig-mohr_coulomb.pdf
+sh -c '\
+rm -f sediment_flux_v*ms.txt; \
+for v in 3.169e-6 9.506e-6; do \
+ for P in 1 1.2742749857031337 1.6237767391887217 2.0691380811147897 2.6366508987303581 3.3598182862837818 4.2813323987193934 5.4555947811685188 6.9519279617756053 8.8586679041008249 11.28837891684689 14.384498882876629 18.329807108324356 23.357214690901213 29.763514416313178 37.926901907322495 48.329302385717519 61.584821106602604 78.475997035146108 100 8 22 50; do \
+ (printf "%ge3\t" $P; ../cngf-pf/shear_flux < strain_distribution_N${P}kPa_v${v}ms.txt) >> sediment_flux_v${v}ms.txt; \
+ echo sediment_flux_v${v}ms.txt; \
+done; done'
+sediment_flux_v3.169e-6ms.txt
+sediment_flux_v3.169e-6ms.txt
+sediment_flux_v3.169e-6ms.txt
+sediment_flux_v3.169e-6ms.txt
+sediment_flux_v3.169e-6ms.txt
+sediment_flux_v3.169e-6ms.txt
+sediment_flux_v3.169e-6ms.txt
+sediment_flux_v3.169e-6ms.txt
+sediment_flux_v3.169e-6ms.txt
+sediment_flux_v3.169e-6ms.txt
+sediment_flux_v3.169e-6ms.txt
+sediment_flux_v3.169e-6ms.txt
+sediment_flux_v3.169e-6ms.txt
+sediment_flux_v3.169e-6ms.txt
+sediment_flux_v3.169e-6ms.txt
+sediment_flux_v3.169e-6ms.txt
+sediment_flux_v3.169e-6ms.txt
+sediment_flux_v3.169e-6ms.txt
+sediment_flux_v3.169e-6ms.txt
+sediment_flux_v3.169e-6ms.txt
+sediment_flux_v3.169e-6ms.txt
+sediment_flux_v3.169e-6ms.txt
+sediment_flux_v3.169e-6ms.txt
+sediment_flux_v9.506e-6ms.txt
+sediment_flux_v9.506e-6ms.txt
+sediment_flux_v9.506e-6ms.txt
+sediment_flux_v9.506e-6ms.txt
+sediment_flux_v9.506e-6ms.txt
+sediment_flux_v9.506e-6ms.txt
+sediment_flux_v9.506e-6ms.txt
+sediment_flux_v9.506e-6ms.txt
+sediment_flux_v9.506e-6ms.txt
+sediment_flux_v9.506e-6ms.txt
+sediment_flux_v9.506e-6ms.txt
+sediment_flux_v9.506e-6ms.txt
+sediment_flux_v9.506e-6ms.txt
+sediment_flux_v9.506e-6ms.txt
+sediment_flux_v9.506e-6ms.txt
+sediment_flux_v9.506e-6ms.txt
+sediment_flux_v9.506e-6ms.txt
+sediment_flux_v9.506e-6ms.txt
+sediment_flux_v9.506e-6ms.txt
+sediment_flux_v9.506e-6ms.txt
+sediment_flux_v9.506e-6ms.txt
+sediment_flux_v9.506e-6ms.txt
+sediment_flux_v9.506e-6ms.txt
+gnuplot fig-sediment_flux.gp > fig-sediment_flux.pdf
+iter chisq delta/lim lambda a b
+ 0 2.9070849779e+08 0.00e+00 1.09e+04 1.000000e+00 1.000000e+00
+ 1 4.3033553429e+07 -5.76e+05 1.09e+03 9.173418e-01 7.884070e-01
+ 2 3.9644105723e+06 -9.85e+05 1.09e+02 4.434906e-01 6.718084e-01
+ 3 1.5355771358e+03 -2.58e+08 1.09e+01 9.386253e-03 6.668584e-01
+ 4 8.0110915384e+00 -1.91e+07 1.09e+00 1.214331e-03 6.474419e-01
+ 5 5.1131762338e-01 -1.47e+06 1.09e-01 1.288157e-03 4.897693e-01
+ 6 2.8436049436e-03 -1.79e+07 1.09e-02 1.676210e-03 3.505343e-01
+ 7 7.3334699222e-05 -3.78e+06 1.09e-03 1.811875e-03 3.286767e-01
+ 8 3.1923482258e-05 -1.30e+05 1.09e-04 1.814273e-03 3.291044e-01
+ 9 3.1923446129e-05 -1.13e-01 1.09e-05 1.814281e-03 3.291025e-01
+iter chisq delta/lim lambda a b
+
+After 9 iterations the fit converged.
+final sum of squares of residuals : 3.19234e-05
+rel. change during last iteration : -1.13175e-06
+
+degrees of freedom (FIT_NDF) : 21
+rms of residuals (FIT_STDFIT) = sqrt(WSSR/ndf) : 0.00123295
+variance of residuals (reduced chisquare) = WSSR/ndf : 1.52016e-06
+
+Final set of parameters Asymptotic Standard Error
+======================= ==========================
+a = 0.00181428 +/- 3.279e-06 (0.1808%)
+b = 0.329102 +/- 0.0005065 (0.1539%)
+
+correlation matrix of the fit parameters:
+ a b
+a 1.000
+b -0.955 1.000