from dust import setup_dust if not os.path.exists("models"): os.mkdir("models") if not os.path.exists("models/indiv"): os.mkdir("models/indiv") seed = -69102 for model_name in ["sky_updarm_gaussian_hole_morepah"]: for spectral_type in spectral_types: print "Processing %s..." % spectral_type m = Model("models/indiv/energy_%s_%s" % (model_name, spectral_type.replace(" ", "_"))) seed += 1 m.set_seed(seed) m = setup_sources(m, model_name, sky_model[model_name][0], spectral_type=spectral_type) m = setup_dust(m, **sky_model[model_name][1]) m.set_n_initial_iterations(3) m.set_n_photons(initial=100000000, imaging=0) m.set_enforce_energy_range(False) m.write()
if not os.path.exists("models"): os.mkdir("models") if not os.path.exists("models/basic"): os.mkdir("models/basic") seed = -43912 for model_name in sky_model: m = Model("models/basic/basic_%s" % model_name) seed += 1 m.set_seed(seed) m = setup_sources(m, model_name, sky_model[model_name][0]) m = setup_dust(m, **sky_model[model_name][1]) image = m.add_peeled_images(image=True, sed=False) image.set_image_size(1, 1) image.set_viewing_angles([90.0], [180.0]) image.set_image_limits(65.0, -65.0, -1.0, 1.0) image.set_inside_observer((8.5 * kpc, 0.0 * kpc, 0.015 * kpc)) image.set_wavelength_range(160, 3.0, 140.0) image.set_track_origin("detailed") image = m.add_peeled_images(image=True, sed=False) image.set_image_size(130, 1) image.set_viewing_angles([90.0], [180.0]) image.set_image_limits(65.0, -65.0, -1.0, 1.0) image.set_inside_observer((8.5 * kpc, 0.0 * kpc, 0.015 * kpc))