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 (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