manus_continuum_granular1_exp

commit 4fb6f93669dbb0ca3846f5f3d6d97e9c38507d2f
parent 3ce6ca648978162c85813629a75c8190bb71e44c
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Tue, 24 Mar 2020 21:03:13 +0100

Add estimate of till flux across simplified grounding line

Diffstat:
A
flux_estimate/Makefile
 | 
29
+++++++++++++++++++++++++++++

1 file changed, 29 insertions(+), 0 deletions(-)
diff --git a/flux_estimate/Makefile b/flux_estimate/Makefile
@@ -0,0 +1,29 @@
+BIN = ../1d_fd_simple_shear/1d_fd_simple_shear
+SHEARFLUX = ../1d_fd_simple_shear/shear_flux
+
+# cross-section width [m]
+WIDTH = 100e3
+
+# average ice-flow velocity [m/s]
+VEL != echo 1000.0/365.0/24.0/3600.0 | bc -l
+
+# flux per unit width [m^2/s]
+SPEC_FLUX != $(BIN) -d 1e-3 -n 100e3 -s $(VEL) | $(SHEARFLUX)
+
+# total sediment volume passing WIDTH [m^3/s]
+FLUX != awk "BEGIN { print $(WIDTH)*$(SPEC_FLUX) }"
+
+# total sediment volume passing WIDTH [m^3/a]
+YEARLY_VOL != awk "BEGIN { print $(FLUX)*3600.0*24.0*365.0 }"
+
+# target sediment volume [m^3]
+TARGET_VOL = 5.34e11
+
+# time required to reach TARGET_VOLUME [a]
+YEAR_DURATION != awk "BEGIN { print $(TARGET_VOL)/$(YEARLY_VOL) }"
+
+default:
+	@printf 'Specific flux: %.5g m^2/s\n' "$(SPEC_FLUX)"
+	@printf 'Flux:          %.5g m^3/s\n' "$(FLUX)"
+	@printf 'Yearly volume: %.5g m^3/a\n' "$(YEARLY_VOL)"
+	@printf 'Duration:      %.5g a\n' "$(YEAR_DURATION)"