cngf-pf-exp1

commit 0e035e180e5a00886eb101441709b74ac4db3869
parent 592a9ca5274d53ca20e1af52618cbbdb1c3c3e38
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Fri, 14 Aug 2020 09:48:12 +0200

Analyze WIS with analytical solution

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

1 file changed, 20 insertions(+), 7 deletions(-)
diff --git a/flux_estimate/Makefile b/flux_estimate/Makefile
@@ -9,20 +9,32 @@ WIDTH = 100e3
 VEL != echo 1000.0/365.0/24.0/3600.0 | bc -l
 
 POROSITY = 0.35
-PERMEABILITY = 4.9e-17
+DIFFUSIVITY = 1e-8
 FREQ_DAILY != echo 1.0/3600.0/24.0 | bc -l
+FREQ_WEEKLY != echo 1.0/3600.0/24.0/7.0 | bc -l
 FREQ_MONTHLY != echo 1.0/3600.0/24.0/30.0 | bc -l
-FREQ_YEARLY != echo 1.0/3600.0/24.0/30.0 | bc -l
-FREQ=$(FREQ_DAILY)
-FLAGS = -m 0.45 -c 3e3 -n 100e3 -d 1e-3 -s $(VEL) -L 4.0 -p $(POROSITY) -D 1e-8 -q $(FREQ)
+FREQ_YEARLY != echo 1.0/3600.0/24.0/365.0 | bc -l
+FREQ_DECADAL != echo 1.0/3600.0/24.0/365.0/10.0 | bc -l
+FREQ_CENTENNIAL != echo 1.0/3600.0/24.0/365.0/100.0 | bc -l
+FREQ=$(FREQ_WEEKLY)
+T_END != echo 3600.0*24.0*14.0 | bc -l
+FLAGS = -m 0.45 -c 3e3 -n 100e3 -d 1e-3 -s $(VEL) -L 4.0 -p $(POROSITY) -D $(DIFFUSIVITY) -q $(FREQ) -I 1800.0
 
 TARGET_VOL = 5.34e11
 
-default: steady_flux variable_flux variable_flux_hi
+default: \
+	steady_flux \
+	def_depth
 
 steady.txt:
 	$(BIN) $(FLAGS) > $@
 
+def_depth:
+	for q in $(FREQ_DAILY) $(FREQ_WEEKLY) $(FREQ_MONTHLY) $(FREQ_YEARLY) $(FREQ_DECADAL) $(FREQ_CENTENNIAL); do \
+		printf '%s\t' "$$q"; \
+		$(MAXDEPTH) -p $(POROSITY) -D $(DIFFUSIVITY) -q $$q -a 50e3;\
+	done
+
 steady_flux: steady.txt
 	# steady-state hydrology
 	@avg_spec_flux=$$($(SHEARFLUX) steady.txt) && \
@@ -38,7 +50,7 @@ steady_flux: steady.txt
 variable.output00095.txt:
 	$(BIN) $(FLAGS) \
 		-F -a 10e3 \
-		-e $$(echo 3600.0*24.0*2.0 | bc -l) -I 1800.0 \
+		-e $(T_END) \
 		variable
 
 variable_flux: variable.output00095.txt
@@ -59,7 +71,7 @@ variable_flux: variable.output00095.txt
 variable_hi.output00095.txt:
 	$(BIN) $(FLAGS) \
 		-F -a 50e3 \
-		-e $$(echo 3600.0*24.0*2.0 | bc -l) -I 1800.0 \
+		-e $(T_END) \
 		variable_hi
 
 variable_flux_hi: variable_hi.output00095.txt
@@ -74,6 +86,7 @@ variable_flux_hi: variable_hi.output00095.txt
 		printf 'Yearly volume: %.5g m^3/a\n' "$$yearly_volume" && \
 		printf 'Duration:      %.5g a\n' "$$duration"
 
+
 truffer_comparison:
 	for A in 10e3 20e3 40e3 80e3 160e3 320e3 640e3 1000e3; do \
 		printf '%s\t' "$$A"; \