cngf-pf

continuum model for granular flows with pore-pressure dynamics (renamed from 1d_fd_simple_shear)
git clone git://src.adamsgaard.dk/cngf-pf # fast
git clone https://src.adamsgaard.dk/cngf-pf.git # slow
Log | Files | Refs | README | LICENSE Back to index

Makefile (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