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