if cfg.par.MANUAL_ORIENTATION == True:
image.set_viewing_angles(np.array(cfg.model.THETA),
np.array(cfg.model.PHI))
else:
image.set_viewing_angles(
np.linspace(0, 90, par.NTHETA).tolist() * par.NPHI,
np.repeat(np.linspace(0, 90, par.NPHI), par.NPHI))
image.set_track_origin('basic')
image.set_image_size(cfg.par.npix_x, cfg.par.npix_y)
image.set_image_limits(-dx / 2., dx / 2., -dy / 2., dy / 2.)
if cfg.par.SKIP_RT == False:
m_imaging.write(model.inputfile + '.image', overwrite=True)
m_imaging.run(model.outputfile + '.image',
mpi=True,
n_processes=par.n_MPI_processes,
overwrite=True)
convolve(model.outputfile + '.image', par.filterfiles, filter_data)
# Print a message in case that skip_rt debugging flag is set:
print('++++++++++++++++++++++++++++++++++++')
print(
'WARNING: SKIP RT is set in the parameters_master file - this is why your code didnt run'
)
print('++++++++++++++++++++++++++++++++++++')
if ds_type in ['gadget_hdf5', 'tipsy', 'arepo_hdf5']:
dump_data(reg, model)
if cfg.par.MANUAL_ORIENTATION == True:
image.set_viewing_angles(np.array(cfg.model.THETA),
np.array(cfg.model.PHI))
else:
image.set_viewing_angles(
np.linspace(0, 90, par.NTHETA).tolist() * par.NPHI,
np.repeat(np.linspace(0, 90, par.NPHI), par.NPHI))
image.set_track_origin('basic')
image.set_image_size(cfg.par.npix_x, cfg.par.npix_y)
image.set_image_limits(-dx / 2., dx / 2., -dy / 2., dy / 2.)
if cfg.par.SKIP_RT == False:
m_imaging.write(model.inputfile + '.image', overwrite=True)
m_imaging.run(model.outputfile + '.image',
mpi=True,
n_processes=par.n_MPI_processes,
overwrite=True)
convolve(model.outputfile + '.image', par.filterfiles, filter_data)
# Print a message in case that skip_rt debugging flag is set:
print('++++++++++++++++++++++++++++++++++++')
print(
'WARNING: SKIP RT is set in the parameters_master file - this is why your code didnt run'
)
print('++++++++++++++++++++++++++++++++++++')
dump_data(pf, model)