<s> import numpy as np <EOL> from hyperion . model import AnalyticalYSOModel <EOL> from hyperion . util . constants import rsun , lsun , au , msun , yr , c <EOL> m = AnalyticalYSOModel ( ) <EOL> wav , fnu = np . loadtxt ( '<STR_LIT>' , unpack = True ) <EOL> nu = c / ( wav * <NUM_LIT> ) <EOL> m . star . radius = <NUM_LIT> * rsun <EOL> m . star . spectrum = ( nu , fnu ) <EOL> m . star . luminosity = lsun <EOL> m . star . mass = <NUM_LIT:0.5> * msun <EOL> disk = m . add_flared_disk ( ) <EOL> disk . mass = <NUM_LIT> * msun <EOL> disk . rmin = <NUM_LIT:7> * m . star . radius <EOL> disk . rmax = <NUM_LIT:200> * au <EOL> disk . r_0 = m . star . radius <EOL> disk . h_0 = <NUM_LIT> * disk . r_0 <EOL> disk . p = - <NUM_LIT:1.0> <EOL> disk . beta = <NUM_LIT> <EOL> disk . dust = '<STR_LIT>' <EOL> envelope = m . add_ulrich_envelope ( ) <EOL> envelope . rc = disk . rmax <EOL> envelope . mdot = <NUM_LIT> * msun / yr <EOL> envelope . rmin = <NUM_LIT:7> * m . star . radius <EOL> envelope . rmax = <NUM_LIT> * au <EOL> envelope . dust = '<STR_LIT>' <EOL> cavity = envelope . add_bipolar_cavity ( ) <EOL> cavity . power = <NUM_LIT> <EOL> cavity . theta_0 = <NUM_LIT:20> <EOL> cavity . r_0 = envelope . rmax <EOL> cavity . rho_0 = <NUM_LIT> * <NUM_LIT> <EOL> cavity . rho_exp = <NUM_LIT:0.> <EOL> cavity . dust = '<STR_LIT>' <EOL> m . set_raytracing ( True ) <EOL> m . set_mrw ( True , gamma = <NUM_LIT> ) <EOL> m . set_spherical_polar_grid_auto ( <NUM_LIT> , <NUM_LIT> , <NUM_LIT:1> ) <EOL> sed = m . add_peeled_images ( sed = True , image = False ) <EOL> sed . set_viewing_angles ( np . linspace ( <NUM_LIT:0.> , <NUM_LIT> , <NUM_LIT:10> ) , np . repeat ( <NUM_LIT> , <NUM_LIT:10> ) ) <EOL> sed . set_wavelength_range ( <NUM_LIT> , <NUM_LIT> , <NUM_LIT> ) <EOL> m . set_n_photons ( initial = <NUM_LIT> , imaging = <NUM_LIT> , <EOL> raytracing_sources = <NUM_LIT> , raytracing_dust = <NUM_LIT> ) <EOL> m . set_n_initial_iterations ( <NUM_LIT:10> ) <EOL> m . set_convergence ( True , percentile = <NUM_LIT> , absolute = <NUM_LIT> , relative = <NUM_LIT> ) <EOL> m . write ( '<STR_LIT>' ) <EOL> m . run ( '<STR_LIT>' , mpi = True ) </s>
