tAdd cohesion to Kamb 1991 example run - cngf-pf - continuum model for granular… | |
git clone git://src.adamsgaard.dk/cngf-pf | |
Log | |
Files | |
Refs | |
README | |
LICENSE | |
--- | |
commit 250cd96441debffc95b28f632f14203a54c91db3 | |
parent da498945df45320ea6c776126e482bda3ac416b4 | |
Author: Anders Damsgaard <[email protected]> | |
Date: Tue, 9 Jul 2019 12:10:29 +0200 | |
Add cohesion to Kamb 1991 example run | |
Diffstat: | |
M examples/1d_fd_simple_shear.png | 0 | |
M examples/1d_fd_simple_shear_rheolo… | 0 | |
M examples/1d_fd_simple_shear_rheolo… | 0 | |
M examples/Makefile | 8 ++++---- | |
4 files changed, 4 insertions(+), 4 deletions(-) | |
--- | |
diff --git a/examples/1d_fd_simple_shear.png b/examples/1d_fd_simple_shear.png | |
Binary files differ. | |
diff --git a/examples/1d_fd_simple_shear_rheology.png b/examples/1d_fd_simple_s… | |
Binary files differ. | |
diff --git a/examples/1d_fd_simple_shear_rheology_kamb.png b/examples/1d_fd_sim… | |
Binary files differ. | |
diff --git a/examples/Makefile b/examples/Makefile | |
t@@ -10,13 +10,13 @@ plots: 1d_fd_simple_shear.png \ | |
diurnal.mp4: diurnal.output00000.txt.png | |
ffmpeg -i diurnal.output%05d.txt.png -y diurnal.mp4 | |
-diurnal.output00000.txt.png: 1d_fd_simple_shear_fluid.gp diurnal.output00000.t… | |
+diurnal.output00000.txt.png: 1d_fd_simple_shear_fluid.gp diurnal.output00000.t… | |
/bin/sh -c '\ | |
for f in diurnal.output*.txt; do \ | |
gnuplot -e "filename=\"$$f\"; p_min=\"0\"; p_max=\"100e3\"" $<… | |
done' | |
-diurnal.output00000.txt: $(BIN) Makefile | |
+diurnal.output00000.txt: $(BIN) | |
/bin/sh -c '\ | |
./$< --resolution 50 --length 8.0 --normal-stress 150e3 \ | |
--fluid --fluid-permeability 2e-17 \ | |
t@@ -58,14 +58,14 @@ diurnal.gif: diurnal.mp4 | |
gnuplot 1d_fd_simple_shear_rheology.gp > $@ | |
# 1 bar is equal to 100 kPa | |
-1d_fd_simple_shear_rheology_kamb.png: $(BIN) 1d_fd_simple_shear_rheology_kamb.… | |
+1d_fd_simple_shear_rheology_kamb.png: $(BIN) 1d_fd_simple_shear_rheology_kamb.… | |
/bin/sh -c '\ | |
for b in 0.01 0.10 0.20 0.40 0.94; do \ | |
out="1d_fd_simple_shear_rheology_b$${b}_kamb.txt"; \ | |
rm -f "$$out"; \ | |
for t in $$(seq 0.01 0.01 2.0); do \ | |
printf "$$t\t" >> "$$out"; \ | |
- ./$< -f 1.1 -P 1.7e3 --stress-ratio $$t -b $$b | \ | |
+ ./$< -f 0.443 -C 1.0e3 -P 1.7e3 --stress-ratio $$t -b … | |
tail -n 1 | cut -f2 >> "$$out"; \ | |
done; done' | |
gnuplot 1d_fd_simple_shear_rheology_kamb.gp > $@ |