Makefile (1294B)
1 BIN = ../1d_fd_simple_shear/1d_fd_simple_shear 2 FIG != basename $(PWD) 3 4 default: ../$(FIG).pdf 5 6 PF_SEQ=1 2 3 4 5 6 7 8 9 7 8 sim-4.output00000.txt: $(BIN) Makefile 9 /bin/sh -c 'for i in $(PF_SEQ); do \ 10 ./$(BIN) -L 8.0 \ 11 -n 200e3 \ 12 -s $$( echo "10000/(365*24*60*60)" | bc -l ) \ 13 -F \ 14 -k 2e-17 \ 15 -O 0.0 \ 16 -a "$$( echo "2*$$i" | bc -l )0e3" \ 17 -q $$( echo "1.0/(3600*24)" | bc -l ) \ 18 -u $$( echo "$$i*3600*24" | bc -l ) \ 19 -S square \ 20 -I $$( echo "60*10" | bc -l ) \ 21 -t $$( echo "($$i-1)*3600*24" | bc -l ) \ 22 -e $$( echo "($$i+1)*3600*24*1.5" | bc -l ) sim-$$i; done' 23 24 timeseries.txt: sim-4.output00000.txt 25 /bin/sh -c '\ 26 rm -f $@; \ 27 for i in $(PF_SEQ); do \ 28 for f in sim-$$i.output*.txt; do \ 29 tail -n 1 "$$f" | cut -f2- >> $@; \ 30 done; done' 31 32 max_strainrate_depth.txt: sim-4.output00000.txt 33 /bin/sh -c '\ 34 rm -f $@; \ 35 for i in $(PF_SEQ); do \ 36 for f in sim-$$i.output*.txt; do \ 37 awk " \ 38 BEGIN{getline; depth=8.0; max=0.0} \ 39 NF{ \ 40 if (\$$6 >= max && \$$6+0 > 1e-300) \ 41 { \ 42 max=\$$6; \ 43 depth=\$$1; \ 44 } \ 45 } \ 46 END{ \ 47 print depth,\"\t\",max \ 48 }" \ 49 "$$f" >> $@; \ 50 done; done' 51 52 ../$(FIG).pdf: fig.gp timeseries.txt max_strainrate_depth.txt 53 gnuplot fig.gp > $@ 54 55 clean: 56 rm -f *.txt 57 rm -f ../$(FIG).pdf 58 59 .PHONY: default clean