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 (4957B)


      1 JULIA=julia --procs 1 --optimize=3 --math-mode=fast --color=yes
      2 #JULIA=julia --procs 1 --color=yes
      3 PREFIX=width_monodisperse
      4 SCRIPT=../simulation.jl
      5 PLOTSCRIPT=../plot.jl
      6 NRUNS=10
      7 ids = \
      8 	$(PREFIX)_mu0.3_sigma_c0kPa_w6e3.pdf \
      9 	$(PREFIX)_mu0.3_sigma_c0kPa_w7e3.pdf \
     10 	$(PREFIX)_mu0.3_sigma_c0kPa_w8e3.pdf \
     11 	$(PREFIX)_mu0.3_sigma_c0kPa_w9e3.pdf \
     12 	$(PREFIX)_mu0.3_sigma_c0kPa_w1e4.pdf \
     13 	$(PREFIX)_mu0.0_sigma_c200kPa_w6e3.pdf \
     14 	$(PREFIX)_mu0.0_sigma_c200kPa_w7e3.pdf \
     15 	$(PREFIX)_mu0.0_sigma_c200kPa_w8e3.pdf \
     16 	$(PREFIX)_mu0.0_sigma_c200kPa_w9e3.pdf \
     17 	$(PREFIX)_mu0.0_sigma_c200kPa_w1e4.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_w6e3.pdf: $(SCRIPT)
     50 	for seed in `seq 1 $(NRUNS)`; do \
     51 		$(JULIA) $< \
     52 		--width=6e3 \
     53 		--r_min=1.0e3 --r_max=1.0e3 \
     54 		--mu_s=0.3 --mu_d=0.3 \
     55 		--mu_s_wall=0.3 --mu_d_wall=0.3 \
     56 		--tensile_strength=0000e3 \
     57 		--seed=$$seed \
     58 		"$@"; \
     59 	done
     60 	$(JULIA) $(PLOTSCRIPT) \
     61 		--nruns=$(NRUNS) \
     62 		"$@"
     63 
     64 $(PREFIX)_mu0.3_sigma_c0kPa_w7e3.pdf: $(SCRIPT)
     65 	for seed in `seq 1 $(NRUNS)`; do \
     66 		$(JULIA) $< \
     67 		--width=7e3 \
     68 		--r_min=1.0e3 --r_max=1.0e3 \
     69 		--mu_s=0.3 --mu_d=0.3 \
     70 		--mu_s_wall=0.3 --mu_d_wall=0.3 \
     71 		--tensile_strength=0000e3 \
     72 		--seed=$$seed \
     73 		"$@"; \
     74 	done
     75 	$(JULIA) $(PLOTSCRIPT) \
     76 		--nruns=$(NRUNS) \
     77 		"$@"
     78 
     79 $(PREFIX)_mu0.3_sigma_c0kPa_w8e3.pdf: $(SCRIPT)
     80 	for seed in `seq 1 $(NRUNS)`; do \
     81 		$(JULIA) $< \
     82 		--width=8e3 \
     83 		--r_min=1.0e3 --r_max=1.0e3 \
     84 		--mu_s=0.3 --mu_d=0.3 \
     85 		--mu_s_wall=0.3 --mu_d_wall=0.3 \
     86 		--tensile_strength=0000e3 \
     87 		--seed=$$seed \
     88 		"$@"; \
     89 	done
     90 	$(JULIA) $(PLOTSCRIPT) \
     91 		--nruns=$(NRUNS) \
     92 		"$@"
     93 
     94 $(PREFIX)_mu0.3_sigma_c0kPa_w9e3.pdf: $(SCRIPT)
     95 	for seed in `seq 1 $(NRUNS)`; do \
     96 		$(JULIA) $< \
     97 		--width=9e3 \
     98 		--r_min=1.0e3 --r_max=1.0e3 \
     99 		--mu_s=0.3 --mu_d=0.3 \
    100 		--mu_s_wall=0.3 --mu_d_wall=0.3 \
    101 		--tensile_strength=0000e3 \
    102 		--seed=$$seed \
    103 		"$@"; \
    104 	done
    105 	$(JULIA) $(PLOTSCRIPT) \
    106 		--nruns=$(NRUNS) \
    107 		"$@"
    108 
    109 $(PREFIX)_mu0.3_sigma_c0kPa_w1e4.pdf: $(SCRIPT)
    110 	for seed in `seq 1 $(NRUNS)`; do \
    111 		$(JULIA) $< \
    112 		--width=1e4 \
    113 		--r_min=1.0e3 --r_max=1.0e3 \
    114 		--mu_s=0.3 --mu_d=0.3 \
    115 		--mu_s_wall=0.3 --mu_d_wall=0.3 \
    116 		--tensile_strength=0000e3 \
    117 		--seed=$$seed \
    118 		"$@"; \
    119 	done
    120 	$(JULIA) $(PLOTSCRIPT) \
    121 		--nruns=$(NRUNS) \
    122 		"$@"
    123 
    124 $(PREFIX)_mu0.0_sigma_c200kPa_w6e3.pdf: $(SCRIPT)
    125 	for seed in `seq 1 $(NRUNS)`; do \
    126 		$(JULIA) $< \
    127 		--width=6e3 \
    128 		--r_min=1.0e3 --r_max=1.0e3 \
    129 		--mu_s=0.0 --mu_d=0.0 \
    130 		--mu_s_wall=0.0 --mu_d_wall=0.0 \
    131 		--tensile_strength=0200e3 \
    132 		--seed=$$seed \
    133 		"$@"; \
    134 	done
    135 	$(JULIA) $(PLOTSCRIPT) \
    136 		--nruns=$(NRUNS) \
    137 		"$@"
    138 
    139 $(PREFIX)_mu0.0_sigma_c200kPa_w7e3.pdf: $(SCRIPT)
    140 	for seed in `seq 1 $(NRUNS)`; do \
    141 		$(JULIA) $< \
    142 		--width=7e3 \
    143 		--r_min=1.0e3 --r_max=1.0e3 \
    144 		--mu_s=0.0 --mu_d=0.0 \
    145 		--mu_s_wall=0.0 --mu_d_wall=0.0 \
    146 		--tensile_strength=0200e3 \
    147 		--seed=$$seed \
    148 		"$@"; \
    149 	done
    150 	$(JULIA) $(PLOTSCRIPT) \
    151 		--nruns=$(NRUNS) \
    152 		"$@"
    153 
    154 $(PREFIX)_mu0.0_sigma_c200kPa_w8e3.pdf: $(SCRIPT)
    155 	for seed in `seq 1 $(NRUNS)`; do \
    156 		$(JULIA) $< \
    157 		--width=8e3 \
    158 		--r_min=1.0e3 --r_max=1.0e3 \
    159 		--mu_s=0.0 --mu_d=0.0 \
    160 		--mu_s_wall=0.0 --mu_d_wall=0.0 \
    161 		--tensile_strength=0200e3 \
    162 		--seed=$$seed \
    163 		"$@"; \
    164 	done
    165 	$(JULIA) $(PLOTSCRIPT) \
    166 		--nruns=$(NRUNS) \
    167 		"$@"
    168 
    169 $(PREFIX)_mu0.0_sigma_c200kPa_w9e3.pdf: $(SCRIPT)
    170 	for seed in `seq 1 $(NRUNS)`; do \
    171 		$(JULIA) $< \
    172 		--width=9e3 \
    173 		--r_min=1.0e3 --r_max=1.0e3 \
    174 		--mu_s=0.0 --mu_d=0.0 \
    175 		--mu_s_wall=0.0 --mu_d_wall=0.0 \
    176 		--tensile_strength=0200e3 \
    177 		--seed=$$seed \
    178 		"$@"; \
    179 	done
    180 	$(JULIA) $(PLOTSCRIPT) \
    181 		--nruns=$(NRUNS) \
    182 		"$@"
    183 
    184 $(PREFIX)_mu0.0_sigma_c200kPa_w1e4.pdf: $(SCRIPT)
    185 	for seed in `seq 1 $(NRUNS)`; do \
    186 		$(JULIA) $< \
    187 		--width=1e4 \
    188 		--r_min=1.0e3 --r_max=1.0e3 \
    189 		--mu_s=0.0 --mu_d=0.0 \
    190 		--mu_s_wall=0.0 --mu_d_wall=0.0 \
    191 		--tensile_strength=0200e3 \
    192 		--seed=$$seed \
    193 		"$@"; \
    194 	done
    195 	$(JULIA) $(PLOTSCRIPT) \
    196 		--nruns=$(NRUNS) \
    197 		"$@"
    198 
    199 compress:
    200 	tar cvfz seed1-vtp.tar.gz $(PREFIX)*seed1/*.vtp
    201 	tar cvfz seed1-vtu.tar.gz $(PREFIX)*seed1/*.vtu
    202 	tar cvfz plots.tar.gz *.png *.pdf
    203 
    204 .PHONY: clean
    205 clean:
    206 	$(RM) $(PREFIX)_*.run
    207 	$(RM) -rf $(PREFIX)_*pdf*/
    208 	$(RM) -rf test*pdf*/
    209 	for f in $(PREFIX)_*.vts; do $(RM) $$f; done
    210 	for f in $(PREFIX)_*.vtp; do $(RM) $$f; done
    211 	for f in $(PREFIX)_*.vtu; do $(RM) $$f; done
    212 	$(RM) $(PREFIX)_*.pdf
    213 	$(RM) $(PREFIX)_*.png
    214 	$(RM) $(PREFIX)_*.txt
    215 	$(RM) test*.pdf
    216 	$(RM) test*.txt