tMakefile - seaice-experiments - sea ice experiments using Granular.jl
(HTM) git clone git://src.adamsgaard.dk/seaice-experiments
(DIR) Log
(DIR) Files
(DIR) Refs
(DIR) README
(DIR) LICENSE
---
tMakefile (5160B)
---
1 JULIA=julia --procs 1 --optimize=3 --math-mode=fast --color=yes
2 #JULIA=julia --procs 1 --color=yes
3 PREFIX=width
4 SCRIPT=../simulation.jl
5 PLOTSCRIPT=../plot.jl
6 NRUNS=1
7 R_MIN=5.0e2
8 R_MAX=1.0e3
9 MU=0.5
10 SIGMA_C=200e3
11
12 ids = \
13 $(PREFIX)_mu$(MU)_sigma_c0kPa_w6e3.pdf \
14 $(PREFIX)_mu$(MU)_sigma_c0kPa_w7e3.pdf \
15 $(PREFIX)_mu$(MU)_sigma_c0kPa_w8e3.pdf \
16 $(PREFIX)_mu$(MU)_sigma_c0kPa_w9e3.pdf \
17 $(PREFIX)_mu$(MU)_sigma_c0kPa_w1e4.pdf \
18 $(PREFIX)_mu0.0_sigma_c$(SIGMA_C)Pa_w6e3.pdf \
19 $(PREFIX)_mu0.0_sigma_c$(SIGMA_C)Pa_w7e3.pdf \
20 $(PREFIX)_mu0.0_sigma_c$(SIGMA_C)Pa_w8e3.pdf \
21 $(PREFIX)_mu0.0_sigma_c$(SIGMA_C)Pa_w9e3.pdf \
22 $(PREFIX)_mu0.0_sigma_c$(SIGMA_C)Pa_w1e4.pdf
23
24
25 .PHONY: all
26 all: $(ids) $(PREFIX)-survived_fraction.png $(PREFIX).png
27
28 all-plots: $(PREFIX)_fits.txt
29
30 $(PREFIX)_fits.txt: $(PLOTSCRIPT)
31 $(RM) $@
32 for id in $(ids); do \
33 echo $$id; \
34 $(JULIA) $< \
35 --nruns=$(NRUNS) \
36 "$$id" >> $@;\
37 done;
38
39 $(PREFIX)-survived_fraction.png: $(ids)
40 gm montage \
41 $(ids:.pdf=.pdf-survived_fraction.pdf.png) \
42 -geometry +4+4 \
43 -trim \
44 -tile $(shell echo "($(words $(ids))+1)/2" | bc)x2 \
45 $@
46
47 $(PREFIX).png: $(ids)
48 gm montage \
49 $(ids:.pdf=.pdf.png) \
50 -geometry +4+4 \
51 -trim \
52 -tile $(shell echo "($(words $(ids))+1)/2" | bc)x2 \
53 $@
54
55 $(PREFIX)_mu$(MU)_sigma_c0kPa_w6e3.pdf: $(SCRIPT)
56 for seed in `seq 1 $(NRUNS)`; do \
57 $(JULIA) $< \
58 --width=6e3 \
59 --r_min=$(R_MIN) --r_max=$(R_MAX) \
60 --mu_s=$(MU) --mu_d=$(MU) \
61 --mu_s_wall=$(MU) --mu_d_wall=$(MU) \
62 --tensile_strength=0000e3 \
63 --seed=$$seed \
64 "$@"; \
65 done
66 $(JULIA) $(PLOTSCRIPT) \
67 --nruns=$(NRUNS) \
68 "$@"
69
70 $(PREFIX)_mu$(MU)_sigma_c0kPa_w7e3.pdf: $(SCRIPT)
71 for seed in `seq 1 $(NRUNS)`; do \
72 $(JULIA) $< \
73 --width=7e3 \
74 --r_min=$(R_MIN) --r_max=$(R_MAX) \
75 --mu_s=$(MU) --mu_d=$(MU) \
76 --mu_s_wall=$(MU) --mu_d_wall=$(MU) \
77 --tensile_strength=0000e3 \
78 --seed=$$seed \
79 "$@"; \
80 done
81 $(JULIA) $(PLOTSCRIPT) \
82 --nruns=$(NRUNS) \
83 "$@"
84
85 $(PREFIX)_mu$(MU)_sigma_c0kPa_w8e3.pdf: $(SCRIPT)
86 for seed in `seq 1 $(NRUNS)`; do \
87 $(JULIA) $< \
88 --width=8e3 \
89 --r_min=$(R_MIN) --r_max=$(R_MAX) \
90 --mu_s=$(MU) --mu_d=$(MU) \
91 --mu_s_wall=$(MU) --mu_d_wall=$(MU) \
92 --tensile_strength=0000e3 \
93 --seed=$$seed \
94 "$@"; \
95 done
96 $(JULIA) $(PLOTSCRIPT) \
97 --nruns=$(NRUNS) \
98 "$@"
99
100 $(PREFIX)_mu$(MU)_sigma_c0kPa_w9e3.pdf: $(SCRIPT)
101 for seed in `seq 1 $(NRUNS)`; do \
102 $(JULIA) $< \
103 --width=9e3 \
104 --r_min=$(R_MIN) --r_max=$(R_MAX) \
105 --mu_s=$(MU) --mu_d=$(MU) \
106 --mu_s_wall=$(MU) --mu_d_wall=$(MU) \
107 --tensile_strength=0000e3 \
108 --seed=$$seed \
109 "$@"; \
110 done
111 $(JULIA) $(PLOTSCRIPT) \
112 --nruns=$(NRUNS) \
113 "$@"
114
115 $(PREFIX)_mu$(MU)_sigma_c0kPa_w1e4.pdf: $(SCRIPT)
116 for seed in `seq 1 $(NRUNS)`; do \
117 $(JULIA) $< \
118 --width=1e4 \
119 --r_min=$(R_MIN) --r_max=$(R_MAX) \
120 --mu_s=$(MU) --mu_d=$(MU) \
121 --mu_s_wall=$(MU) --mu_d_wall=$(MU) \
122 --tensile_strength=0000e3 \
123 --seed=$$seed \
124 "$@"; \
125 done
126 $(JULIA) $(PLOTSCRIPT) \
127 --nruns=$(NRUNS) \
128 "$@"
129
130 $(PREFIX)_mu0.0_sigma_c200kPa_w6e3.pdf: $(SCRIPT)
131 for seed in `seq 1 $(NRUNS)`; do \
132 $(JULIA) $< \
133 --width=6e3 \
134 --r_min=$(R_MIN) --r_max=$(R_MAX) \
135 --mu_s=0.0 --mu_d=0.0 \
136 --mu_s_wall=0.0 --mu_d_wall=0.0 \
137 --tensile_strength=$(SIGMA_C) \
138 --seed=$$seed \
139 "$@"; \
140 done
141 $(JULIA) $(PLOTSCRIPT) \
142 --nruns=$(NRUNS) \
143 "$@"
144
145 $(PREFIX)_mu0.0_sigma_c200kPa_w7e3.pdf: $(SCRIPT)
146 for seed in `seq 1 $(NRUNS)`; do \
147 $(JULIA) $< \
148 --width=7e3 \
149 --r_min=$(R_MIN) --r_max=$(R_MAX) \
150 --mu_s=0.0 --mu_d=0.0 \
151 --mu_s_wall=0.0 --mu_d_wall=0.0 \
152 --tensile_strength=$(SIGMA_C) \
153 --seed=$$seed \
154 "$@"; \
155 done
156 $(JULIA) $(PLOTSCRIPT) \
157 --nruns=$(NRUNS) \
158 "$@"
159
160 $(PREFIX)_mu0.0_sigma_c200kPa_w8e3.pdf: $(SCRIPT)
161 for seed in `seq 1 $(NRUNS)`; do \
162 $(JULIA) $< \
163 --width=8e3 \
164 --r_min=$(R_MIN) --r_max=$(R_MAX) \
165 --mu_s=0.0 --mu_d=0.0 \
166 --mu_s_wall=0.0 --mu_d_wall=0.0 \
167 --tensile_strength=$(SIGMA_C) \
168 --seed=$$seed \
169 "$@"; \
170 done
171 $(JULIA) $(PLOTSCRIPT) \
172 --nruns=$(NRUNS) \
173 "$@"
174
175 $(PREFIX)_mu0.0_sigma_c200kPa_w9e3.pdf: $(SCRIPT)
176 for seed in `seq 1 $(NRUNS)`; do \
177 $(JULIA) $< \
178 --width=9e3 \
179 --r_min=$(R_MIN) --r_max=$(R_MAX) \
180 --mu_s=0.0 --mu_d=0.0 \
181 --mu_s_wall=0.0 --mu_d_wall=0.0 \
182 --tensile_strength=$(SIGMA_C) \
183 --seed=$$seed \
184 "$@"; \
185 done
186 $(JULIA) $(PLOTSCRIPT) \
187 --nruns=$(NRUNS) \
188 "$@"
189
190 $(PREFIX)_mu0.0_sigma_c200kPa_w1e4.pdf: $(SCRIPT)
191 for seed in `seq 1 $(NRUNS)`; do \
192 $(JULIA) $< \
193 --width=1e4 \
194 --r_min=$(R_MIN) --r_max=$(R_MAX) \
195 --mu_s=0.0 --mu_d=0.0 \
196 --mu_s_wall=0.0 --mu_d_wall=0.0 \
197 --tensile_strength=$(SIGMA_C) \
198 --seed=$$seed \
199 "$@"; \
200 done
201 $(JULIA) $(PLOTSCRIPT) \
202 --nruns=$(NRUNS) \
203 "$@"
204
205 compress:
206 tar cvfz seed1-vtp.tar.gz $(PREFIX)*seed1/*.vtp
207 tar cvfz seed1-vtu.tar.gz $(PREFIX)*seed1/*.vtu
208 tar cvfz plots.tar.gz *.png *.pdf
209
210 .PHONY: clean
211 clean:
212 $(RM) $(PREFIX)_*.run
213 $(RM) -rf $(PREFIX)_*pdf*/
214 $(RM) -rf test*pdf*/
215 for f in $(PREFIX)_*.vts; do $(RM) $$f; done
216 for f in $(PREFIX)_*.vtp; do $(RM) $$f; done
217 for f in $(PREFIX)_*.vtu; do $(RM) $$f; done
218 $(RM) $(PREFIX)_*.pdf
219 $(RM) $(PREFIX)_*.png
220 $(RM) $(PREFIX)_*.txt
221 $(RM) test*.pdf
222 $(RM) test*.txt