crop_maps=2**12,
epsilon=0.,
verbose=False)
elif params.projection == 'healpix':
## Save data on disk into fits file for later use in xpure
name_out = 'sim{:03d}_{}_{}_{}'.format(args.sim_number, params.tag,
params.name_instrument,
params.name_strategy)
write_maps_a_la_xpure(sky_out_tot,
name_out=name_out,
output_path='xpure/maps')
write_weights_a_la_xpure(sky_out_tot,
name_out=name_out,
output_path='xpure/masks',
epsilon=0.08,
HWP=False)
if args.inifile_xpure is not None:
## Import parameters from the user parameter file
params_xpure = import_string_as_module(args.inifile_xpure)
batch_file = 'sim{:03d}_{}_{}_{}.batch'.format(
args.sim_number, params.tag, params.name_instrument,
params.name_strategy)
create_batch(batch_file, name_out, params, params_xpure)
qsub = commands.getoutput('sbatch ' + batch_file)
print(qsub)
if rank == 0:
if params.projection == 'flat':
name_out = '{}_{}_{}'.format(params.tag, params.name_instrument,
params.name_strategy)
sky_out_tot.pickle_me('{}/sim{:03d}_{}.pkl'.format(
args.folder_out, args.sim_number, name_out),
shrink_maps=False,
crop_maps=2**12,
epsilon=0.,
verbose=False)
elif params.projection == 'healpix':
## Save data on disk into fits file for later use in xpure
name_out = 'sim{:03d}_{}_{}_{}'.format(args.sim_number, params.tag,
params.name_instrument,
params.name_strategy)
write_maps_a_la_xpure(sky_out_tot,
name_out=name_out,
output_path='{}/maps'.format(
args.folder_out))
write_weights_a_la_xpure(sky_out_tot,
name_out=name_out,
output_path='{}/masks'.format(
args.folder_out),
epsilon=0.08,
HWP=False)
MPI.COMM_WORLD.barrier()