SED_Modeling/tracy_benchmark.py at master · mjrfringes/SED_Modeling · GitHub

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
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
import models as m
import pickle
import time
from hyperion.util.convenience import OptThinRadius
'''
name = 'WL16-Ulrich-Env_Mass=0.1-Rc=2000-Lsun=5.6'
folder = 'WL16/'
WL16 = m.YSOModelSim(name=name,folder=folder,L_sun=5.6,env_mass=0.1,env_type='ulrich',rc=2000)
WL16.runModel()
WL16.plotModel(extinction=9.5,show=True)
WL16.modelDump()
'''
T=3000
env_rmin_list=[0.1]
env_rmax=14000
rc_list=[139]
l_sun_list=[0.8]#[25]
m_sun=0.5#2
env_mass_list=[0.05]
mdisk=0.004
disk_rmax=100
h0 =0.01
cav_theta=15.
disk_rmin=0.1
envtype='ulrich'
env_power=0.0
envmdot=3.5785507e-05
amb_rmin=1
amb_rmax=env_rmax
dustfilelist=['OH5.hdf5']
amb_dens = 1e-22
angles=[0.,10.,20.,30.,40.,50.,60.,70.,80.,90.]
angles2=[0.,0.,0.,0.,0.,0.,0.,0.,0.,0.]
foldername='WL6'
env=True
cav=True
ext=0
inc=1
cav_rho_0 = 1.6e-19
for env_rmin in env_rmin_list:
	for l_sun in l_sun_list:
		for env_mass in env_mass_list:
			for rc in rc_list:
				for dustfile in dustfilelist:
					start = time.time()
					name = foldername+'-'+envtype+'-Env_Mass='+str(env_mass)+'-Rmin='+str(env_rmin)+'-Rmax='+str(env_rmax)+\
					'-Lsun='+str(l_sun)+'-Msun='+str(m_sun)+'-pow='+str(env_power)+'-h0='+str(h0)+'-mdisk='+str(mdisk)+ \
					'-rc='+str(rc)+'-cav='+str(cav_theta)+'-disk_rmin='+str(disk_rmin)+'-envmdot='+str(envmdot)+ \
					'-disk_rmax='+str(disk_rmax)+'-T='+str(T)+'-'+dustfile.split('.')[0]+'-angle='+str(angles[0])+ \
					'-ext='+str(ext)+'-amb_rmax='+str(amb_rmax)+'-kmh'+'-OH5' #+'-nodisk'
					###
					folder = foldername+'/'
					#WL16=m.modelLoad(folder,name)

					WL6 = m.YSOModelSim(name=name,folder=folder,L_sun=l_sun,M_sun=m_sun,env_mass=env_mass,env_type=envtype,disk_rmax=disk_rmax,\
							env_rmin=env_rmin,rc=rc,env_power=env_power,disk_h_0=h0,disk_mass=mdisk,\
							cav=cav,cav_theta=cav_theta,disk_rmin=disk_rmin,mdot=envmdot,amb_rmin=amb_rmin,amb_rmax=amb_rmax,\
							env_rmax=env_rmax,T=T,dustfile=dustfile,amb_dens=amb_dens,angles=angles,angles2=angles2,env=env,cav_rho_0=cav_rho_0)
					WL6.runModel()
					print 'Elapsed time: ',time.time()-start,'s'
					WL6.plotModel(extinction=ext,show=True,inc=inc,sourcename='None')
					#time.sleep(2)
					#m.modelDump(WL16)


print 'Script is finished!'
'''
env_rmin_list=[80]
env_rmax=8000
rc_list=[76]
l_sun_list=[0.5]
env_mass_list=[0.1]
mdisk=0.0002
disk_rmax=100
h0 = 20
cav_theta=0.
disk_rmin=1
envtype='power'
env_power=-0.5
envmdot=2e-5
amb_rmin=1
dustfilelist=['OH5.hdf5']
amb_dens = 2e-22
angles=[0.,10.,20.,30.,40.,50.,60.,70.,80.,90.]
angles2=[0.,0.,0.,0.,0.,0.,0.,0.,0.,0.]
foldername='WL6'
env=True
cav=False
ext=0
'''