em_fname = fname
for i in range(len(sys.argv)):
if sys.argv[i] == "-grid":
zscale = float(sys.argv[i + 1])
nz = int(sys.argv[i + 2])
nogrid_read = False
em_fname = fname + "_scale1e" + str(np.log10(zscale))
else:
fname = "test"
fname_to_open = fname + ".pf"
specfname = fname + ".spec"
print "ZEPHYR: models SV93 wind for parameter file", fname_to_open
sv.input_params(fname_to_open)
#################################################
# print out some simple luminosities
ledd = 1.3e38 * geo.mstar / MSOL
print "The eddington luminosity for this black hole is ", ledd
medd = ledd / (0.1 * C * C)
print "Which implies an accretion rate of ", medd / (MSOL / YR)
edd_frac = (geo.disk_mdot) / medd
print "This geometry is accreting at an eddington fraction of ", edd_frac
omega = np.cos(geo.sv_thetamin) - np.cos(geo.sv_thetamax)
if tau < 1e-6: escape = 1. elif tau < 10.0: escape = (1. - np.exp (-tau)) / tau else: escape = 1. / tau return escape ''' fname_to_open = sys.argv[1] sv.input_params(fname = fname_to_open) #define grid using input file parameters and user specified zscale, otherwise set to rstar if nogrid_read: zscale=geo.rstar/10.0 nz=30 xscale=geo.sv_rmin xcen,zcen = geo.get_grid(zscale, xscale, nz, nx, geo.wind_rmax, 1) xup, zup = geo.get_grid(zscale, xscale, nz, nx, geo.wind_rmax, 0)
em_fname = fname
for i in range(len(sys.argv)):
if sys.argv[i] == "-grid":
zscale = float(sys.argv[i + 1])
nz = int(sys.argv[i + 2])
nogrid_read = False
em_fname = fname + '_scale1e' + str(np.log10(zscale))
else:
fname = "test"
fname_to_open = fname + ".pf"
specfname = fname + ".spec"
print 'ZEPHYR: models SV93 wind for parameter file', fname_to_open
sv.input_params(fname_to_open)
#################################################
# print out some simple luminosities
ledd = 1.3e38 * geo.mstar / MSOL
print "The eddington luminosity for this black hole is ", ledd
medd = ledd / (0.1 * C * C)
print "Which implies an accretion rate of ", medd / (MSOL / YR)
edd_frac = (geo.disk_mdot) / medd
print "This geometry is accreting at an eddington fraction of ", edd_frac
omega = np.cos(geo.sv_thetamin) - np.cos(geo.sv_thetamax)
print "With these angles, the wind will subtend a fraction of ", omega, " possible viewing angles"
def P_escape(tau):
if tau < 1e-6:
escape = 1.
elif tau < 10.0:
escape = (1. - np.exp(-tau)) / tau
else:
escape = 1. / tau
return escape
'''
fname_to_open = sys.argv[1]
sv.input_params(fname = fname_to_open)
#define grid using input file parameters and user specified zscale, otherwise set to rstar
if nogrid_read:
zscale=geo.rstar/10.0
nz=30
xscale=geo.sv_rmin
xcen,zcen = geo.get_grid(zscale, xscale, nz, nx, geo.wind_rmax, 1)
xup, zup = geo.get_grid(zscale, xscale, nz, nx, geo.wind_rmax, 0)
