tMakefile - cngf-pf - continuum model for granular flows with pore-pressure dynamics (renamed from 1d_fd_simple_shear)
(HTM) git clone git://src.adamsgaard.dk/cngf-pf
(DIR) Log
(DIR) Files
(DIR) Refs
(DIR) README
(DIR) LICENSE
---
tMakefile (1198B)
---
1 BIN = ../cngf-pf
2 BIN2 = ../max_depth_simple_shear
3
4 PERMEABILITY = 2e-17
5 GRAIN_COMPRESSIBILITY = 1e-8
6 FLUID_PRESSURE_TOP = 100e3
7 #FLUID_PRESSURE_AMPL = 120e3
8 FLUID_PRESSURE_AMPL = 80e3
9 FLUID_PRESSURE_FREQ != echo "1.0/(3600*24)" | bc -l
10 NORMAL_STRESS_TOP = 200e3
11
12 timeseries.pdf: fig.gp timeseries.txt
13 sh -c '\
14 out=$$($(BIN2) \
15 -k $(PERMEABILITY) \
16 -P $(GRAIN_COMPRESSIBILITY) \
17 -O $(FLUID_PRESSURE_TOP) \
18 -a $(FLUID_PRESSURE_AMPL) \
19 -q $(FLUID_PRESSURE_FREQ) \
20 ); \
21 echo $$out; \
22 max_def_depth=$$(echo $$out | cut -d" " -f1); \
23 skin_depth=$$(echo $$out | cut -d" " -f2); \
24 gnuplot -e "skin_depth=$$skin_depth; max_def_depth=$$max_def_depth" \
25 fig.gp > $@'
26
27 timeseries.txt: sim.output00000.txt
28 /bin/sh -c '\
29 rm -f $@; \
30 for f in sim.output*.txt; do \
31 tail -n 1 "$$f" | cut -f2- >> $@; \
32 done'
33
34 sim.output00000.txt: $(BIN)
35 /bin/sh -c '\
36 $(BIN) -L 8.0 \
37 -n $(NORMAL_STRESS_TOP) \
38 -f 0.40 \
39 -F \
40 -k $(PERMEABILITY) \
41 -P $(GRAIN_COMPRESSIBILITY) \
42 -O $(FLUID_PRESSURE_TOP) \
43 -a $(FLUID_PRESSURE_AMPL) \
44 -q $(FLUID_PRESSURE_FREQ) \
45 -I $$( echo "60*10" | bc -l ) \
46 -e $$( echo "3600*24*7" | bc -l ) sim'
47
48 clean:
49 rm -f *.txt timeseries.pdf
50
51 .PHONY: clean