tAdd high-amplitude scenario, calculate result for variable runs - cngf-pf-exp1 - experiments for first paper with continuum granular model
(HTM) git clone git://src.adamsgaard.dk/manus_continuum_granular1_exp
(DIR) Log
(DIR) Files
(DIR) Refs
(DIR) Submodules
(DIR) README
(DIR) LICENSE
---
(DIR) commit bb4859557c8a495cbe32920a5be9ecc08cb81691
(DIR) parent 755abffdf095b7154fa940a4dd9a40db21d2cc7b
(HTM) Author: Anders Damsgaard <anders@adamsgaard.dk>
Date: Thu, 26 Mar 2020 12:25:55 +0100
Add high-amplitude scenario, calculate result for variable runs
Diffstat:
M flux_estimate/Makefile | 33 ++++++++++++++++++++++++++++---
1 file changed, 30 insertions(+), 3 deletions(-)
---
(DIR) diff --git a/flux_estimate/Makefile b/flux_estimate/Makefile
t@@ -45,11 +45,38 @@ variable.output00048.txt:
-e $$(echo 3600.0*24.0*2.0 | bc -l) -I 1800.0 \
variable
-variable_flux: variable.output00048.txt
- $(SHEAFLUX) variable.output*.txt > $@
+variable_flux:
+ # variable water-pressure, low amplitude
+ @avg_spec_flux=$$($(SHEARFLUX) variable.output*.txt | awk '{sum+=$$1} END{print sum/NR}') && \
+ avg_flux=$$(awk "BEGIN { print $(WIDTH)*$$avg_spec_flux }" ) && \
+ yearly_volume=$$(awk "BEGIN { print $$avg_flux*3600.0*24.0*365.0 }" ) && \
+ duration=$$(awk "BEGIN { print $(TARGET_VOL)/$$yearly_volume }" ) && \
+ printf 'Specific flux: %.5g m^2/s\n' "$$avg_spec_flux" && \
+ printf 'Avg. flux: %.5g m^3/s\n' "$$avg_flux" && \
+ printf 'Yearly volume: %.5g a\n' "$$yearly_volume" && \
+ printf 'Duration: %.5g a\n' "$$duration"
+
+variable_hi.output00048.txt:
+ $(BIN) $(FLAGS) \
+ -F -a 50e3 \
+ -q $(FREQ) \
+ -k $(PERMEABILITY) \
+ -e $$(echo 3600.0*24.0*2.0 | bc -l) -I 1800.0 \
+ variable_hi
+
+variable_flux_hi:
+ # variable water-pressure, high amplitude
+ @avg_spec_flux=$$($(SHEARFLUX) variable.output*.txt | awk '{sum+=$$1} END{print sum/NR}') && \
+ avg_flux=$$(awk "BEGIN { print $(WIDTH)*$$avg_spec_flux }" ) && \
+ yearly_volume=$$(awk "BEGIN { print $$avg_flux*3600.0*24.0*365.0 }" ) && \
+ duration=$$(awk "BEGIN { print $(TARGET_VOL)/$$yearly_volume }" ) && \
+ printf 'Specific flux: %.5g m^2/s\n' "$$avg_spec_flux" && \
+ printf 'Avg. flux: %.5g m^3/s\n' "$$avg_flux" && \
+ printf 'Yearly volume: %.5g a\n' "$$yearly_volume" && \
+ printf 'Duration: %.5g a\n' "$$duration"
amplitude_dependence:
- for A in 10e3 20e3 40e3 80e3 160e3 320e3; do \
+ for A in 10e3 20e3 30e3 40e3 50e3 100e3; do \
$(MAXDEPTH) -p $(POROSITY) -k $(PERMEABILITY) -a $$A -q $(FREQ); \
done