manus_continuum_granular1_exp

commit 4d693b1149fc9172c9bf769d3e472fa874f6d461
parent 7e2f54e17a93aafd6af63fc6ba5d7b4e05bd205c
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Wed, 25 Mar 2020 18:53:03 +0100

Add preliminary setup for variable water-pressure run

Diffstat:
M
flux_estimate/Makefile
 | 
38
++++++++++++++++++++++++++++----------

1 file changed, 28 insertions(+), 10 deletions(-)
diff --git a/flux_estimate/Makefile b/flux_estimate/Makefile
@@ -7,23 +7,41 @@ 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)
+FLAGS = -m 0.45 -c 3e3 -n 100e3 -d 1e-3 -n 100e3 -s $(VEL) -L 4.0 -p 0.35
+
+# flux per unit width under constant hydrology [m^2/s]
+SPEC_FLUX_STEADY != $(BIN) $(FLAGS) | $(SHEARFLUX)
 
 # total sediment volume passing WIDTH [m^3/s]
-FLUX != awk "BEGIN { print $(WIDTH)*$(SPEC_FLUX) }"
+FLUX_STEADY != awk "BEGIN { print $(WIDTH)*$(SPEC_FLUX_STEADY) }"
 
 # total sediment volume passing WIDTH [m^3/a]
-YEARLY_VOL != awk "BEGIN { print $(FLUX)*3600.0*24.0*365.0 }"
+YEARLY_VOL_STEADY != awk "BEGIN { print $(FLUX_STEADY)*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) }"
+YEAR_DURATION_STEADY != awk "BEGIN { print $(TARGET_VOL)/$(YEARLY_VOL_STEADY) }"
+
+default: steady variable_flux
+
+steady:
+	# steady state hydrology
+	@printf 'Specific flux: %.5g m^2/s\n' "$(SPEC_FLUX_STEADY)"
+	@printf 'Flux:          %.5g m^3/s\n' "$(FLUX_STEADY)"
+	@printf 'Yearly volume: %.5g m^3/a\n' "$(YEARLY_VOL_STEADY)"
+	@printf 'Duration:      %.5g a\n' "$(YEAR_DURATION_STEADY)"
+
+variable_flux:
+	$(BIN) $(FLAGS) \
+		-F -a 10e3 -q $$(echo 1.0/3600.0/24.0 | bc -l) -k 4.9e-17 \
+		-e $$(echo 3600.0*24.0*2.0 | bc -l) -I 1800.0 \
+		variable
+	$(SHEAFLUX) variable.output*.txt > $@
+
+clean:
+	find . -type f -name 'variable.output*.txt' -delete
+	rm -f variable_flux
 
-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)"
+.PHONY: default steady clean