tMakefile - cngf-pf-exp1 - experiments for first paper with continuum granular model
 (HTM) git clone git://src.adamsgaard.dk/manus_continuum_granular1_exp
 (DIR) Log
 (DIR) Files
 (DIR) Refs
 (DIR) Submodules
 (DIR) README
 (DIR) LICENSE
       ---
       tMakefile (1274B)
       ---
            1 BIN = ../cngf-pf/cngf-pf
            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 triangle \
           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