forked from nx1/ModellingAccretionDisks
-
Notifications
You must be signed in to change notification settings - Fork 0
/
batch.py
45 lines (37 loc) · 1.59 KB
/
batch.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
# -*- coding: utf-8 -*-
import itertools
N_l = [1,5,10,30]
const_l = [1.1,1.5]
Q_factor_l = [0.005,0.025,0.1]
tMax_factor_l = [1.1,1.4]
H_R_l= [0.1,1.0]
M_0_start_l = [1.0,10.0]
biglist= list(itertools.product(N_l,const_l,Q_factor_l,tMax_factor_l,H_R_l,M_0_start_l ))
for i in range(182,len(biglist)):
print i
print biglist [i]
N = biglist[i][0]
const = biglist[i][1]
Q_factor = biglist[i][2]
tMax_factor = biglist[i][3]
H_R_= biglist[i][4]
M_0_start = biglist[i][5]
execfile("discfunc.py")
name0 = N_l.index(N)
name1 = const_l.index(const)
name2 = Q_factor_l.index(Q_factor)
name3 = tMax_factor_l.index(tMax_factor)
name4 = H_R_l.index(H_R_)
name5 = M_0_start_l.index(M_0_start)
NAMEpng = '%s%s%s%s%s%s.png' % (name0,name1,name2,name3,name4,name5)
NAMEeps = '%s%s%s%s%s%s.eps' % (name0,name1,name2,name3,name4,name5)
#LC
fig1.savefig('variable_graphs/LC/' + NAMEeps, format='eps', dpi=300, bbox_inches='tight', pad_inches = 0.0)
fig1.savefig('variable_graphs/LC/' + NAMEpng, format='png', dpi=300, bbox_inches='tight', pad_inches = 0.0)
#RMS
fig2.savefig('variable_graphs/RMS/' + NAMEeps, format='eps', dpi=300, bbox_inches='tight', pad_inches = 0.0)
fig2.savefig('variable_graphs/RMS/' + NAMEpng, format='png', dpi=300, bbox_inches='tight', pad_inches = 0.0)
#PSD
fig3.savefig('variable_graphs/PSD/' + NAMEeps, format='eps', dpi=300, bbox_inches='tight', pad_inches = 0.0)
fig3.savefig('variable_graphs/PSD/' + NAMEpng, format='png', dpi=300, bbox_inches='tight', pad_inches = 0.0)
print NAMEpng