Makefile (1158B)
1 BIN = ../1d_fd_simple_shear/1d_fd_simple_shear 2 FIG = $(shell basename $(CURDIR)) 3 4 default: ../$(FIG).pdf 5 6 # --stress-ratio 0.40 \ 7 8 sim.output00000.txt: $(BIN) 9 /bin/sh -c '\ 10 ./$< --length 8.0 --normal-stress 150e3 \ 11 --limit-shear-velocity $$( echo "3000/(365*24*60*60)" | bc -l ) \ 12 --fluid --fluid-permeability 2e-17 --fluid-pressure-top 50e3 \ 13 --fluid-pressure-ampl 50e3 \ 14 --fluid-pressure-freq $$( echo "1.0/(3600*24)" | bc -l ) \ 15 --file-interval $$( echo "60*10" | bc -l ) \ 16 --time-end $$( echo "3600*24*7" | bc -l ) sim' 17 18 timeseries.txt: sim.output00000.txt 19 /bin/sh -c '\ 20 rm -f $@; \ 21 for f in sim.output*.txt; do \ 22 tail -n 1 "$$f" | cut -f2- >> $@; \ 23 done' 24 25 max_strainrate_depth.txt: sim.output00000.txt 26 /bin/sh -c '\ 27 rm -f $@; \ 28 for f in sim.output*.txt; do \ 29 awk " \ 30 BEGIN{getline; depth=8.0; max=0.0} \ 31 NF{ \ 32 if (\$$6 >= max && \$$6+0 > 1e-300) \ 33 { \ 34 max=\$$6; \ 35 depth=\$$1; \ 36 } \ 37 } \ 38 END{ \ 39 print depth,\"\t\",max \ 40 }" \ 41 "$$f" >> $@; \ 42 done' 43 44 ../$(FIG).pdf: fig.gp timeseries.txt max_strainrate_depth.txt 45 gnuplot $< > $@ 46 47 clean: 48 $(RM) *.txt 49 $(RM) ../$(FIG).pdf 50 51 .PHONY: default clean