Makefile (2568B)
1 BIN = ../cngf-pf/cngf-pf 2 FIG != basename $(PWD) 3 4 # default parameter values 5 LENGTH = 1.0 6 NORMAL_STRESS = 100e3 7 SHEAR_VEL != echo '300/(365*24*60*60)' | bc -l 8 DIAMETER = 1e-4 9 COHESION = 0.0 10 AMPL = 0.4 11 BRATE = 0.9 12 MU_S = 0.4 13 14 default: ../$(FIG).pdf 15 16 d.txt: $(BIN) 17 /bin/sh -c '\ 18 rm -rf $@; \ 19 for d in 1e-4 2e-4 4e-4 1e-3 2e-3 4e-3 1e-2 2e-2; do \ 20 printf "$$d\t" >> $@; \ 21 ./$(BIN) \ 22 -L $(LENGTH) \ 23 -s $(SHEAR_VEL) \ 24 -n $(NORMAL_STRESS) \ 25 -c $(COHESION) \ 26 -A $(AMPL) \ 27 -b $(BRATE) \ 28 -m $(MU_S) \ 29 -d $$d \ 30 d.$$d \ 31 | tail -n 1 | cut -f5 >> $@; \ 32 done' 33 34 mu_s.txt: $(BIN) 35 /bin/sh -c '\ 36 rm -rf $@; \ 37 for mu_s in $$(awk "BEGIN{for (i=0.20; i<0.70; i+=0.05) print i}"); do \ 38 printf "$$mu_s\t" >> $@; \ 39 ./$(BIN) \ 40 -L $(LENGTH) \ 41 -s $(SHEAR_VEL) \ 42 -n $(NORMAL_STRESS) \ 43 -c $(COHESION) \ 44 -A $(AMPL) \ 45 -b $(BRATE) \ 46 -m $$mu_s \ 47 -d $(DIAMETER) \ 48 mu_s.$$mu_s \ 49 | tail -n 1 | cut -f5 >> $@; \ 50 done' 51 52 C.txt: $(BIN) 53 /bin/sh -c '\ 54 rm -rf $@; \ 55 for C in 0.0 5e3 10e3 15e3 20e3 25e3 30e3; do \ 56 printf "$$C\t" >> $@; \ 57 ./$(BIN) \ 58 -L $(LENGTH) \ 59 -s $(SHEAR_VEL) \ 60 -n $(NORMAL_STRESS) \ 61 -c $$C \ 62 -A $(AMPL) \ 63 -b $(BRATE) \ 64 -m $(MU_S) \ 65 -d $(DIAMETER) \ 66 C.$$C \ 67 | tail -n 1 | cut -f5 >> $@; \ 68 done' 69 70 A.txt: $(BIN) 71 /bin/sh -c '\ 72 rm -rf $@; \ 73 for A in 1e-10 $$(awk "BEGIN{for(i=0.1; i<=1.0; i+=0.1) print i}"); do \ 74 printf "$$A\t" >> $@; \ 75 ./$(BIN) \ 76 -L $(LENGTH) \ 77 -s $(SHEAR_VEL) \ 78 -n $(NORMAL_STRESS) \ 79 -c $(COHESION) \ 80 -A $$A \ 81 -b $(BRATE) \ 82 -m $(MU_S) \ 83 -d $(DIAMETER) \ 84 A.$$A \ 85 | tail -n 1 | cut -f5 >> $@; \ 86 done' 87 88 b.txt: $(BIN) 89 /bin/sh -c '\ 90 rm -rf $@; \ 91 for b in $$(awk "BEGIN{for(i=0.1; i<=0.9; i+=0.1) print i}"); do \ 92 printf "$$b\t" >> $@; \ 93 ./$(BIN) \ 94 -L $(LENGTH) \ 95 -s $(SHEAR_VEL) \ 96 -n $(NORMAL_STRESS) \ 97 -c $(COHESION) \ 98 -A $(AMPL) \ 99 -b $$b \ 100 -m $(MU_S) \ 101 -d $(DIAMETER) \ 102 b.$$b \ 103 | tail -n 1 | cut -f5 >> $@; \ 104 done' 105 106 L_z.txt: $(BIN) 107 /bin/sh -c '\ 108 rm -rf $@; \ 109 for L_z in 0.002 0.004 0.008 0.01 0.015 0.02 0.04 0.06 0.08 0.10 0.15 0.20 0.25 0.30 0.40 0.50 0.60 0.70 0.75 0.80 0.90 1.0; do \ 110 printf "$$L_z\t" >> $@; \ 111 ./$(BIN) \ 112 -L $$L_z \ 113 -s $(SHEAR_VEL) \ 114 -n $(NORMAL_STRESS) \ 115 -c $(COHESION) \ 116 -A $(AMPL) \ 117 -b $(BRATE) \ 118 -m $(MU_S) \ 119 -d $(DIAMETER) \ 120 L_z.$$L_z \ 121 | tail -n 1 | cut -f5 >> $@; \ 122 done' 123 124 125 ../$(FIG).pdf: fig.gp d.txt mu_s.txt C.txt A.txt b.txt L_z.txt 126 gnuplot fig.gp > $@ 127 128 clean: 129 rm -f *.txt 130 rm -f ../$(FIG).pdf 131 132 .PHONY: default clean