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)