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 (1362B)
       ---
            1 BIN = ../cngf-pf/cngf-pf
            2 FIG != basename $(PWD)
            3 
            4 default: ../$(FIG).pdf
            5 
            6 PERM_SEQ=3.16e-16 1e-16 3.16e-17 1e-17 3.16e-18 1e-18 3.16e-19 1e-19
            7 
            8 sim-1e-18.output00000.txt: $(BIN)
            9         /bin/sh -c 'for perm in $(PERM_SEQ); do \
           10         ./$(BIN) -L 8.0 \
           11         -n 200e3 \
           12         -s $$( echo "10000/(365*24*60*60)" | bc -l ) \
           13         -F \
           14         -k $$perm \
           15         -O 0.0 \
           16         -a 180e3 \
           17         -q $$( echo "1.0/(3600*24)" | bc -l ) \
           18         -u $$( echo "2*3600*24" | bc -l ) \
           19         -S triangle \
           20         -I $$( echo "60*10" | bc -l ) \
           21         -t 0.0 \
           22         -e $$( echo "20*3600*24 + 60*10" | bc -l ) sim-$$perm; done'
           23 
           24 timeseries.txt: sim-1e-18.output00000.txt
           25         /bin/sh -c '\
           26         rm -f $@; \
           27         for perm in $(PERM_SEQ); do \
           28                 for f in sim-$$perm.output*.txt; do \
           29                         tail -n 1 "$$f" | cut -f2- >> $@; \
           30         done; done'
           31 
           32 max_strainrate_depth.txt: sim-1e-18.output00000.txt
           33         /bin/sh -c '\
           34         rm -f $@; \
           35         for perm in $(PERM_SEQ); do \
           36                 for f in sim-$$perm.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,\"\t\",$$perm \
           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         for perm in $(PERM_SEQ); do rm -f sim-$$perm.output*.txt; done
           57         rm -f *.txt
           58         rm -f ../$(FIG).pdf
           59 
           60 .PHONY: default clean