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 (4989B)
---
1 JULIA=julia --procs 1 --optimize=3 --math-mode=fast --color=yes
2 #JULIA=julia --procs 1 --color=yes
3 PREFIX=psd_width
4 SCRIPT=../simulation.jl
5 PLOTSCRIPT=../plot.jl
6 NRUNS=10
7 ids = \
8 $(PREFIX)_mu0.3_sigma_c0kPa_wider.pdf \
9 $(PREFIX)_mu0.3_sigma_c0kPa_wide.pdf \
10 $(PREFIX)_mu0.3_sigma_c0kPa_normal.pdf \
11 $(PREFIX)_mu0.3_sigma_c0kPa_narrow.pdf \
12 $(PREFIX)_mu0.3_sigma_c0kPa_narrower.pdf \
13 $(PREFIX)_mu0.0_sigma_c200kPa_wider.pdf \
14 $(PREFIX)_mu0.0_sigma_c200kPa_wide.pdf \
15 $(PREFIX)_mu0.0_sigma_c200kPa_normal.pdf \
16 $(PREFIX)_mu0.0_sigma_c200kPa_narrow.pdf \
17 $(PREFIX)_mu0.0_sigma_c200kPa_narrower.pdf
18
19 .PHONY: all
20 all: $(ids) $(PREFIX)-survived_fraction.png $(PREFIX).png
21
22 all-plots: $(PREFIX)_fits.txt
23
24 $(PREFIX)_fits.txt: $(PLOTSCRIPT)
25 $(RM) $@
26 for id in $(ids); do \
27 echo $$id; \
28 $(JULIA) $< \
29 --nruns=$(NRUNS) \
30 "$$id" >> $@;\
31 done;
32
33 $(PREFIX)-survived_fraction.png: $(ids)
34 gm montage \
35 $(ids:.pdf=.pdf-survived_fraction.pdf.png) \
36 -geometry +4+4 \
37 -trim \
38 -tile $(shell echo "($(words $(ids))+1)/2" | bc)x2 \
39 $@
40
41 $(PREFIX).png: $(ids)
42 gm montage \
43 $(ids:.pdf=.pdf.png) \
44 -geometry +4+4 \
45 -trim \
46 -tile $(shell echo "($(words $(ids))+1)/2" | bc)x2 \
47 $@
48
49 $(PREFIX)_mu0.3_sigma_c0kPa_wider.pdf: $(SCRIPT)
50 for seed in `seq 1 $(NRUNS)`; do \
51 $(JULIA) $< \
52 --mu_s=0.3 --mu_d=0.3 \
53 --mu_s_wall=0.3 --mu_d_wall=0.3 \
54 --tensile_strength=0000e3 \
55 --gamma_n=0. --gamma_t=0. \
56 --r_min=2.0e2 --r_max=1.75e3 \
57 --seed=$$seed \
58 "$@"; \
59 done
60 $(JULIA) $(PLOTSCRIPT) \
61 --nruns=$(NRUNS) \
62 "$@"
63
64 $(PREFIX)_mu0.3_sigma_c0kPa_wide.pdf: $(SCRIPT)
65 for seed in `seq 1 $(NRUNS)`; do \
66 $(JULIA) $< \
67 --mu_s=0.3 --mu_d=0.3 \
68 --mu_s_wall=0.3 --mu_d_wall=0.3 \
69 --tensile_strength=0000e3 \
70 --gamma_n=0. --gamma_t=0. \
71 --r_min=4.0e2 --r_max=1.55e3 \
72 --seed=$$seed \
73 "$@"; \
74 done
75 $(JULIA) $(PLOTSCRIPT) \
76 --nruns=$(NRUNS) \
77 "$@"
78
79 $(PREFIX)_mu0.3_sigma_c0kPa_normal.pdf: $(SCRIPT)
80 for seed in `seq 1 $(NRUNS)`; do \
81 $(JULIA) $< \
82 --mu_s=0.3 --mu_d=0.3 \
83 --mu_s_wall=0.3 --mu_d_wall=0.3 \
84 --tensile_strength=0000e3 \
85 --gamma_n=0. --gamma_t=0. \
86 --r_min=6.0e2 --r_max=1.35e3 \
87 --seed=$$seed \
88 "$@"; \
89 done
90 $(JULIA) $(PLOTSCRIPT) \
91 --nruns=$(NRUNS) \
92 "$@"
93
94 $(PREFIX)_mu0.3_sigma_c0kPa_narrow.pdf: $(SCRIPT)
95 for seed in `seq 1 $(NRUNS)`; do \
96 $(JULIA) $< \
97 --mu_s=0.3 --mu_d=0.3 \
98 --mu_s_wall=0.3 --mu_d_wall=0.3 \
99 --tensile_strength=0000e3 \
100 --gamma_n=0. --gamma_t=0. \
101 --r_min=8.0e2 --r_max=1.15e3 \
102 --seed=$$seed \
103 "$@"; \
104 done
105 $(JULIA) $(PLOTSCRIPT) \
106 --nruns=$(NRUNS) \
107 "$@"
108
109 $(PREFIX)_mu0.3_sigma_c0kPa_narrower.pdf: $(SCRIPT)
110 for seed in `seq 1 $(NRUNS)`; do \
111 $(JULIA) $< \
112 --mu_s=0.3 --mu_d=0.3 \
113 --mu_s_wall=0.3 --mu_d_wall=0.3 \
114 --tensile_strength=0000e3 \
115 --gamma_n=0. --gamma_t=0. \
116 --r_min=1.0e3 --r_max=1.0e3 \
117 --seed=$$seed \
118 "$@"; \
119 done
120 $(JULIA) $(PLOTSCRIPT) \
121 --nruns=$(NRUNS) \
122 "$@"
123
124
125 $(PREFIX)_mu0.0_sigma_c200kPa_wider.pdf: $(SCRIPT)
126 for seed in `seq 1 $(NRUNS)`; do \
127 $(JULIA) $< \
128 --mu_s=0.0 --mu_d=0.0 \
129 --mu_s_wall=0.0 --mu_d_wall=0.0 \
130 --tensile_strength=0200e3 \
131 --gamma_n=0. --gamma_t=0. \
132 --r_min=2.0e2 --r_max=1.75e3 \
133 --seed=$$seed \
134 "$@"; \
135 done
136 $(JULIA) $(PLOTSCRIPT) \
137 --nruns=$(NRUNS) \
138 "$@"
139
140 $(PREFIX)_mu0.0_sigma_c200kPa_wide.pdf: $(SCRIPT)
141 for seed in `seq 1 $(NRUNS)`; do \
142 $(JULIA) $< \
143 --mu_s=0.0 --mu_d=0.0 \
144 --mu_s_wall=0.0 --mu_d_wall=0.0 \
145 --tensile_strength=0200e3 \
146 --gamma_n=0. --gamma_t=0. \
147 --r_min=4.0e2 --r_max=1.55e3 \
148 --seed=$$seed \
149 "$@"; \
150 done
151 $(JULIA) $(PLOTSCRIPT) \
152 --nruns=$(NRUNS) \
153 "$@"
154
155 $(PREFIX)_mu0.0_sigma_c200kPa_normal.pdf: $(SCRIPT)
156 for seed in `seq 1 $(NRUNS)`; do \
157 $(JULIA) $< \
158 --mu_s=0.0 --mu_d=0.0 \
159 --mu_s_wall=0.0 --mu_d_wall=0.0 \
160 --tensile_strength=0200e3 \
161 --gamma_n=0. --gamma_t=0. \
162 --r_min=6.0e2 --r_max=1.35e3 \
163 --seed=$$seed \
164 "$@"; \
165 done
166 $(JULIA) $(PLOTSCRIPT) \
167 --nruns=$(NRUNS) \
168 "$@"
169
170 $(PREFIX)_mu0.0_sigma_c200kPa_narrow.pdf: $(SCRIPT)
171 for seed in `seq 1 $(NRUNS)`; do \
172 $(JULIA) $< \
173 --mu_s=0.0 --mu_d=0.0 \
174 --mu_s_wall=0.0 --mu_d_wall=0.0 \
175 --tensile_strength=0200e3 \
176 --gamma_n=0. --gamma_t=0. \
177 --r_min=8.0e2 --r_max=1.15e3 \
178 --seed=$$seed \
179 "$@"; \
180 done
181 $(JULIA) $(PLOTSCRIPT) \
182 --nruns=$(NRUNS) \
183 "$@"
184
185 $(PREFIX)_mu0.0_sigma_c200kPa_narrower.pdf: $(SCRIPT)
186 for seed in `seq 1 $(NRUNS)`; do \
187 $(JULIA) $< \
188 --mu_s=0.0 --mu_d=0.0 \
189 --mu_s_wall=0.0 --mu_d_wall=0.0 \
190 --tensile_strength=0200e3 \
191 --gamma_n=0. --gamma_t=0. \
192 --r_min=1.0e3 --r_max=1.0e3 \
193 --seed=$$seed \
194 "$@"; \
195 done
196 $(JULIA) $(PLOTSCRIPT) \
197 --nruns=$(NRUNS) \
198 "$@"
199
200 .PHONY: clean
201 clean:
202 $(RM) $(PREFIX)_*.run
203 $(RM) -rf $(PREFIX)_*pdf*/
204 $(RM) -rf test*pdf*/
205 for f in $(PREFIX)_*.vts; do $(RM) $$f; done
206 for f in $(PREFIX)_*.vtp; do $(RM) $$f; done
207 for f in $(PREFIX)_*.vtu; do $(RM) $$f; done
208 $(RM) $(PREFIX)_*.pdf
209 $(RM) $(PREFIX)_*.png
210 $(RM) $(PREFIX)_*.txt
211 $(RM) test*.pdf
212 $(RM) test*.txt