seaice-experiments

sea ice experiments using Granular.jl
git clone git://src.adamsgaard.dk/seaice-experiments # fast
git clone https://src.adamsgaard.dk/seaice-experiments.git # slow
Log | Files | Refs | README | LICENSE Back to index

Makefile (8696B)


      1 JULIA=julia --procs 1 --optimize=3 --math-mode=fast --color=yes
      2 #JULIA=julia --procs 1 --color=yes
      3 PREFIX=thickness
      4 SCRIPT=../simulation.jl
      5 PLOTSCRIPT=../plot.jl
      6 NRUNS=8
      7 ids = \
      8 	$(PREFIX)_mu0.3_sigma_c0kPa_h01.00.pdf \
      9 	$(PREFIX)_mu0.3_sigma_c0kPa_h01.29.pdf \
     10 	$(PREFIX)_mu0.3_sigma_c0kPa_h01.67.pdf \
     11 	$(PREFIX)_mu0.3_sigma_c0kPa_h02.15.pdf \
     12 	$(PREFIX)_mu0.3_sigma_c0kPa_h02.78.pdf \
     13 	$(PREFIX)_mu0.3_sigma_c0kPa_h03.59.pdf \
     14 	$(PREFIX)_mu0.3_sigma_c0kPa_h04.64.pdf \
     15 	$(PREFIX)_mu0.3_sigma_c0kPa_h05.99.pdf \
     16 	$(PREFIX)_mu0.3_sigma_c0kPa_h07.74.pdf \
     17 	$(PREFIX)_mu0.3_sigma_c0kPa_h10.00.pdf \
     18 	$(PREFIX)_mu0.0_sigma_c200kPa_h01.00.pdf \
     19 	$(PREFIX)_mu0.0_sigma_c200kPa_h01.29.pdf \
     20 	$(PREFIX)_mu0.0_sigma_c200kPa_h01.67.pdf \
     21 	$(PREFIX)_mu0.0_sigma_c200kPa_h02.15.pdf \
     22 	$(PREFIX)_mu0.0_sigma_c200kPa_h02.78.pdf \
     23 	$(PREFIX)_mu0.0_sigma_c200kPa_h03.59.pdf \
     24 	$(PREFIX)_mu0.0_sigma_c200kPa_h04.64.pdf \
     25 	$(PREFIX)_mu0.0_sigma_c200kPa_h05.99.pdf \
     26 	$(PREFIX)_mu0.0_sigma_c200kPa_h07.74.pdf \
     27 	$(PREFIX)_mu0.0_sigma_c200kPa_h10.00.pdf
     28 
     29 .PHONY: all
     30 all: $(ids) $(PREFIX)-survived_fraction.png $(PREFIX).png
     31 
     32 all-plots: $(PREFIX)_fits.txt
     33 
     34 $(PREFIX)_fits.txt: $(PLOTSCRIPT)
     35 	$(RM) $@
     36 	for id in $(ids); do \
     37 		echo $$id; \
     38 		$(JULIA) $< \
     39 		--nruns=$(NRUNS) \
     40 		"$$id" >> $@;\
     41 	done;
     42 
     43 $(PREFIX)-survived_fraction.png: $(ids)
     44 	gm montage \
     45 		$(ids:.pdf=.pdf-survived_fraction.pdf.png) \
     46 	-geometry +4+4 \
     47 	-trim \
     48 	-tile $(shell echo "($(words $(ids))+1)/2" | bc)x2 \
     49 	$@
     50 
     51 $(PREFIX).png: $(ids)
     52 	gm montage \
     53 		$(ids:.pdf=.pdf.png) \
     54 	-geometry +4+4 \
     55 	-trim \
     56 	-tile $(shell echo "($(words $(ids))+1)/2" | bc)x2 \
     57 	$@
     58 
     59 $(PREFIX)_mu0.3_sigma_c0kPa_h01.00.pdf: $(SCRIPT)
     60 	for seed in `seq 1 $(NRUNS)`; do \
     61 		$(JULIA) $< \
     62 		--mu_s=0.3 --mu_d=0.3 \
     63 		--mu_s_wall=0.3 --mu_d_wall=0.3 \
     64 		--tensile_strength=0000e3 \
     65 		--r_min=6.0e2 --r_max=1.35e3 \
     66 		--thickness=1.00 \
     67 		--seed=$$seed \
     68 		"$@"; \
     69 	done
     70 	$(JULIA) $(PLOTSCRIPT) \
     71 		--nruns=$(NRUNS) \
     72 		"$@"
     73 
     74 $(PREFIX)_mu0.3_sigma_c0kPa_h01.29.pdf: $(SCRIPT)
     75 	for seed in `seq 1 $(NRUNS)`; do \
     76 		$(JULIA) $< \
     77 		--mu_s=0.3 --mu_d=0.3 \
     78 		--mu_s_wall=0.3 --mu_d_wall=0.3 \
     79 		--tensile_strength=0000e3 \
     80 		--r_min=6.0e2 --r_max=1.35e3 \
     81 		--thickness=1.29 \
     82 		--seed=$$seed \
     83 		"$@"; \
     84 	done
     85 	$(JULIA) $(PLOTSCRIPT) \
     86 		--nruns=$(NRUNS) \
     87 		"$@"
     88 
     89 $(PREFIX)_mu0.3_sigma_c0kPa_h01.67.pdf: $(SCRIPT)
     90 	for seed in `seq 1 $(NRUNS)`; do \
     91 		$(JULIA) $< \
     92 		--mu_s=0.3 --mu_d=0.3 \
     93 		--mu_s_wall=0.3 --mu_d_wall=0.3 \
     94 		--tensile_strength=0000e3 \
     95 		--r_min=6.0e2 --r_max=1.35e3 \
     96 		--thickness=1.67 \
     97 		--seed=$$seed \
     98 		"$@"; \
     99 	done
    100 	$(JULIA) $(PLOTSCRIPT) \
    101 		--nruns=$(NRUNS) \
    102 		"$@"
    103 
    104 $(PREFIX)_mu0.3_sigma_c0kPa_h02.15.pdf: $(SCRIPT)
    105 	for seed in `seq 1 $(NRUNS)`; do \
    106 		$(JULIA) $< \
    107 		--mu_s=0.3 --mu_d=0.3 \
    108 		--mu_s_wall=0.3 --mu_d_wall=0.3 \
    109 		--tensile_strength=0000e3 \
    110 		--r_min=6.0e2 --r_max=1.35e3 \
    111 		--thickness=2.15 \
    112 		--seed=$$seed \
    113 		"$@"; \
    114 	done
    115 	$(JULIA) $(PLOTSCRIPT) \
    116 		--nruns=$(NRUNS) \
    117 		"$@"
    118 
    119 $(PREFIX)_mu0.3_sigma_c0kPa_h02.78.pdf: $(SCRIPT)
    120 	for seed in `seq 1 $(NRUNS)`; do \
    121 		$(JULIA) $< \
    122 		--mu_s=0.3 --mu_d=0.3 \
    123 		--mu_s_wall=0.3 --mu_d_wall=0.3 \
    124 		--tensile_strength=0000e3 \
    125 		--r_min=6.0e2 --r_max=1.35e3 \
    126 		--thickness=2.78 \
    127 		--seed=$$seed \
    128 		"$@"; \
    129 	done
    130 	$(JULIA) $(PLOTSCRIPT) \
    131 		--nruns=$(NRUNS) \
    132 		"$@"
    133 
    134 $(PREFIX)_mu0.3_sigma_c0kPa_h03.59.pdf: $(SCRIPT)
    135 	for seed in `seq 1 $(NRUNS)`; do \
    136 		$(JULIA) $< \
    137 		--mu_s=0.3 --mu_d=0.3 \
    138 		--mu_s_wall=0.3 --mu_d_wall=0.3 \
    139 		--tensile_strength=0000e3 \
    140 		--r_min=6.0e2 --r_max=1.35e3 \
    141 		--thickness=3.59 \
    142 		--seed=$$seed \
    143 		"$@"; \
    144 	done
    145 	$(JULIA) $(PLOTSCRIPT) \
    146 		--nruns=$(NRUNS) \
    147 		"$@"
    148 
    149 $(PREFIX)_mu0.3_sigma_c0kPa_h04.64.pdf: $(SCRIPT)
    150 	for seed in `seq 1 $(NRUNS)`; do \
    151 		$(JULIA) $< \
    152 		--mu_s=0.3 --mu_d=0.3 \
    153 		--mu_s_wall=0.3 --mu_d_wall=0.3 \
    154 		--tensile_strength=0000e3 \
    155 		--r_min=6.0e2 --r_max=1.35e3 \
    156 		--thickness=4.64 \
    157 		--seed=$$seed \
    158 		"$@"; \
    159 	done
    160 	$(JULIA) $(PLOTSCRIPT) \
    161 		--nruns=$(NRUNS) \
    162 		"$@"
    163 
    164 $(PREFIX)_mu0.3_sigma_c0kPa_h05.99.pdf: $(SCRIPT)
    165 	for seed in `seq 1 $(NRUNS)`; do \
    166 		$(JULIA) $< \
    167 		--mu_s=0.3 --mu_d=0.3 \
    168 		--mu_s_wall=0.3 --mu_d_wall=0.3 \
    169 		--tensile_strength=0000e3 \
    170 		--r_min=6.0e2 --r_max=1.35e3 \
    171 		--thickness=5.99 \
    172 		--seed=$$seed \
    173 		"$@"; \
    174 	done
    175 	$(JULIA) $(PLOTSCRIPT) \
    176 		--nruns=$(NRUNS) \
    177 		"$@"
    178 
    179 $(PREFIX)_mu0.3_sigma_c0kPa_h07.74.pdf: $(SCRIPT)
    180 	for seed in `seq 1 $(NRUNS)`; do \
    181 		$(JULIA) $< \
    182 		--mu_s=0.3 --mu_d=0.3 \
    183 		--mu_s_wall=0.3 --mu_d_wall=0.3 \
    184 		--tensile_strength=0000e3 \
    185 		--r_min=6.0e2 --r_max=1.35e3 \
    186 		--thickness=7.74 \
    187 		--seed=$$seed \
    188 		"$@"; \
    189 	done
    190 	$(JULIA) $(PLOTSCRIPT) \
    191 		--nruns=$(NRUNS) \
    192 		"$@"
    193 
    194 $(PREFIX)_mu0.3_sigma_c0kPa_h10.00.pdf: $(SCRIPT)
    195 	for seed in `seq 1 $(NRUNS)`; do \
    196 		$(JULIA) $< \
    197 		--mu_s=0.3 --mu_d=0.3 \
    198 		--mu_s_wall=0.3 --mu_d_wall=0.3 \
    199 		--tensile_strength=0000e3 \
    200 		--r_min=6.0e2 --r_max=1.35e3 \
    201 		--thickness=10.00 \
    202 		--seed=$$seed \
    203 		"$@"; \
    204 	done
    205 	$(JULIA) $(PLOTSCRIPT) \
    206 		--nruns=$(NRUNS) \
    207 		"$@"
    208 
    209 
    210 
    211 
    212 $(PREFIX)_mu0.0_sigma_c200kPa_h01.00.pdf: $(SCRIPT)
    213 	for seed in `seq 1 $(NRUNS)`; do \
    214 		$(JULIA) $< \
    215 		--mu_s=0.0 --mu_d=0.0 \
    216 		--mu_s_wall=0.0 --mu_d_wall=0.0 \
    217 		--tensile_strength=0200e3 \
    218 		--r_min=6.0e2 --r_max=1.35e3 \
    219 		--thickness=1.00 \
    220 		--seed=$$seed \
    221 		"$@"; \
    222 	done
    223 	$(JULIA) $(PLOTSCRIPT) \
    224 		--nruns=$(NRUNS) \
    225 		"$@"
    226 
    227 $(PREFIX)_mu0.0_sigma_c200kPa_h01.29.pdf: $(SCRIPT)
    228 	for seed in `seq 1 $(NRUNS)`; do \
    229 		$(JULIA) $< \
    230 		--mu_s=0.0 --mu_d=0.0 \
    231 		--mu_s_wall=0.0 --mu_d_wall=0.0 \
    232 		--tensile_strength=0200e3 \
    233 		--r_min=6.0e2 --r_max=1.35e3 \
    234 		--thickness=1.29 \
    235 		--seed=$$seed \
    236 		"$@"; \
    237 	done
    238 	$(JULIA) $(PLOTSCRIPT) \
    239 		--nruns=$(NRUNS) \
    240 		"$@"
    241 
    242 $(PREFIX)_mu0.0_sigma_c200kPa_h01.67.pdf: $(SCRIPT)
    243 	for seed in `seq 1 $(NRUNS)`; do \
    244 		$(JULIA) $< \
    245 		--mu_s=0.0 --mu_d=0.0 \
    246 		--mu_s_wall=0.0 --mu_d_wall=0.0 \
    247 		--tensile_strength=0200e3 \
    248 		--r_min=6.0e2 --r_max=1.35e3 \
    249 		--thickness=1.67 \
    250 		--seed=$$seed \
    251 		"$@"; \
    252 	done
    253 	$(JULIA) $(PLOTSCRIPT) \
    254 		--nruns=$(NRUNS) \
    255 		"$@"
    256 
    257 $(PREFIX)_mu0.0_sigma_c200kPa_h02.15.pdf: $(SCRIPT)
    258 	for seed in `seq 1 $(NRUNS)`; do \
    259 		$(JULIA) $< \
    260 		--mu_s=0.0 --mu_d=0.0 \
    261 		--mu_s_wall=0.0 --mu_d_wall=0.0 \
    262 		--tensile_strength=0200e3 \
    263 		--r_min=6.0e2 --r_max=1.35e3 \
    264 		--thickness=2.15 \
    265 		--seed=$$seed \
    266 		"$@"; \
    267 	done
    268 	$(JULIA) $(PLOTSCRIPT) \
    269 		--nruns=$(NRUNS) \
    270 		"$@"
    271 
    272 $(PREFIX)_mu0.0_sigma_c200kPa_h02.78.pdf: $(SCRIPT)
    273 	for seed in `seq 1 $(NRUNS)`; do \
    274 		$(JULIA) $< \
    275 		--mu_s=0.0 --mu_d=0.0 \
    276 		--mu_s_wall=0.0 --mu_d_wall=0.0 \
    277 		--tensile_strength=0200e3 \
    278 		--r_min=6.0e2 --r_max=1.35e3 \
    279 		--thickness=2.78 \
    280 		--seed=$$seed \
    281 		"$@"; \
    282 	done
    283 	$(JULIA) $(PLOTSCRIPT) \
    284 		--nruns=$(NRUNS) \
    285 		"$@"
    286 
    287 $(PREFIX)_mu0.0_sigma_c200kPa_h03.59.pdf: $(SCRIPT)
    288 	for seed in `seq 1 $(NRUNS)`; do \
    289 		$(JULIA) $< \
    290 		--mu_s=0.0 --mu_d=0.0 \
    291 		--mu_s_wall=0.0 --mu_d_wall=0.0 \
    292 		--tensile_strength=0200e3 \
    293 		--r_min=6.0e2 --r_max=1.35e3 \
    294 		--thickness=3.59 \
    295 		--seed=$$seed \
    296 		"$@"; \
    297 	done
    298 	$(JULIA) $(PLOTSCRIPT) \
    299 		--nruns=$(NRUNS) \
    300 		"$@"
    301 
    302 $(PREFIX)_mu0.0_sigma_c200kPa_h04.64.pdf: $(SCRIPT)
    303 	for seed in `seq 1 $(NRUNS)`; do \
    304 		$(JULIA) $< \
    305 		--mu_s=0.0 --mu_d=0.0 \
    306 		--mu_s_wall=0.0 --mu_d_wall=0.0 \
    307 		--tensile_strength=0200e3 \
    308 		--r_min=6.0e2 --r_max=1.35e3 \
    309 		--thickness=4.64 \
    310 		--seed=$$seed \
    311 		"$@"; \
    312 	done
    313 	$(JULIA) $(PLOTSCRIPT) \
    314 		--nruns=$(NRUNS) \
    315 		"$@"
    316 
    317 $(PREFIX)_mu0.0_sigma_c200kPa_h05.99.pdf: $(SCRIPT)
    318 	for seed in `seq 1 $(NRUNS)`; do \
    319 		$(JULIA) $< \
    320 		--mu_s=0.0 --mu_d=0.0 \
    321 		--mu_s_wall=0.0 --mu_d_wall=0.0 \
    322 		--tensile_strength=0200e3 \
    323 		--r_min=6.0e2 --r_max=1.35e3 \
    324 		--thickness=5.99 \
    325 		--seed=$$seed \
    326 		"$@"; \
    327 	done
    328 	$(JULIA) $(PLOTSCRIPT) \
    329 		--nruns=$(NRUNS) \
    330 		"$@"
    331 
    332 $(PREFIX)_mu0.0_sigma_c200kPa_h07.74.pdf: $(SCRIPT)
    333 	for seed in `seq 1 $(NRUNS)`; do \
    334 		$(JULIA) $< \
    335 		--mu_s=0.0 --mu_d=0.0 \
    336 		--mu_s_wall=0.0 --mu_d_wall=0.0 \
    337 		--tensile_strength=0200e3 \
    338 		--r_min=6.0e2 --r_max=1.35e3 \
    339 		--thickness=7.74 \
    340 		--seed=$$seed \
    341 		"$@"; \
    342 	done
    343 	$(JULIA) $(PLOTSCRIPT) \
    344 		--nruns=$(NRUNS) \
    345 		"$@"
    346 
    347 $(PREFIX)_mu0.0_sigma_c200kPa_h10.00.pdf: $(SCRIPT)
    348 	for seed in `seq 1 $(NRUNS)`; do \
    349 		$(JULIA) $< \
    350 		--mu_s=0.0 --mu_d=0.0 \
    351 		--mu_s_wall=0.0 --mu_d_wall=0.0 \
    352 		--tensile_strength=0200e3 \
    353 		--r_min=6.0e2 --r_max=1.35e3 \
    354 		--thickness=10.00 \
    355 		--seed=$$seed \
    356 		"$@"; \
    357 	done
    358 	$(JULIA) $(PLOTSCRIPT) \
    359 		--nruns=$(NRUNS) \
    360 		"$@"
    361 
    362 
    363 .PHONY: clean
    364 clean:
    365 	$(RM) $(PREFIX)_*.run
    366 	$(RM) -rf $(PREFIX)_*pdf*/
    367 	$(RM) -rf test*pdf*/
    368 	for f in $(PREFIX)_*.vts; do $(RM) $$f; done
    369 	for f in $(PREFIX)_*.vtp; do $(RM) $$f; done
    370 	for f in $(PREFIX)_*.vtu; do $(RM) $$f; done
    371 	$(RM) $(PREFIX)_*.pdf
    372 	$(RM) $(PREFIX)_*.txt
    373 	$(RM) test*.pdf
    374 	$(RM) test*.txt