cngf-pf-exp1

commit 755abffdf095b7154fa940a4dd9a40db21d2cc7b
parent e8da30a680fa43cddf55379aafa7bec6ab833fb8
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Thu, 26 Mar 2020 11:18:24 +0100

Add analysis of amplitude effect

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

1 file changed, 13 insertions(+), 2 deletions(-)
diff --git a/flux_estimate/Makefile b/flux_estimate/Makefile
@@ -1,5 +1,6 @@
 BIN = ../1d_fd_simple_shear/1d_fd_simple_shear
 SHEARFLUX = ../1d_fd_simple_shear/shear_flux
+MAXDEPTH = ../1d_fd_simple_shear/max_depth_simple_shear
 
 # grounding-line width [m]
 WIDTH = 100e3
@@ -7,7 +8,10 @@ WIDTH = 100e3
 # average ice-flow velocity [m/s]
 VEL != echo 1000.0/365.0/24.0/3600.0 | bc -l
 
-FLAGS = -m 0.45 -c 3e3 -n 100e3 -d 1e-3 -n 100e3 -s $(VEL) -L 4.0 -p 0.35
+POROSITY = 0.35
+PERMEABILITY = 4.9e-17
+FREQ != echo 1.0/3600.0/24.0 | bc -l
+FLAGS = -m 0.45 -c 3e3 -n 100e3 -d 1e-3 -n 100e3 -s $(VEL) -L 4.0 -p $(POROSITY)
 
 # flux per unit width under constant hydrology [m^2/s]
 SPEC_FLUX_STEADY != $(BIN) $(FLAGS) | $(SHEARFLUX)
@@ -35,13 +39,20 @@ steady:
 
 variable.output00048.txt:
 	$(BIN) $(FLAGS) \
-		-F -a 10e3 -q $$(echo 1.0/3600.0/24.0 | bc -l) -k 4.9e-17 \
+		-F -a 10e3 \
+		-q $(FREQ) \
+		-k $(PERMEABILITY) \
 		-e $$(echo 3600.0*24.0*2.0 | bc -l) -I 1800.0 \
 		variable
 
 variable_flux: variable.output00048.txt
 	$(SHEAFLUX) variable.output*.txt > $@
 
+amplitude_dependence:
+	for A in 10e3 20e3 40e3 80e3 160e3 320e3; do \
+		$(MAXDEPTH) -p $(POROSITY) -k $(PERMEABILITY) -a $$A -q $(FREQ); \
+	done
+
 clean:
 	find . -type f -name 'variable.output*.txt' -delete
 	rm -f variable_flux