parser.add_argument("-t", "--tsize", type=int, default=360)
parser.add_argument("-p", "--pad", type=int, default=60)
parser.add_argument("-a", "--apod-edge", type=int, default=30)
parser.add_argument("-S", "--nsigma", type=float, default=4)
parser.add_argument("--model-nsigma", type=float, default=4)
parser.add_argument("-l", "--lmin", type=float, default=1000)
parser.add_argument("--debug-tile", type=str, default=None)
parser.add_argument("--nmax", type=int, default=None)
parser.add_argument("-v", "--verbose", action="count", default=3)
parser.add_argument("-q", "--quiet", action="count", default=0)
parser.add_argument("-c", "--cont", action="store_true")
parser.add_argument("-P", "--npass", type=int, default=3)
parser.add_argument("--output-full-model", action="store_true")
args = parser.parse_args()
config = jointmap.read_config(args.config)
mapinfo = jointmap.Mapset(config, args.sel)
tsize = args.tsize # pixels
pad = args.pad # pixels
dtype = np.float64
ncomp = 1
comm = mpi.COMM_WORLD
signals = args.signals.split(",")
verbosity = args.verbose - args.quiet
highpass_alpha = 3
utils.mkdir(args.odir)
debug_tile = None if args.debug_tile is None else [int(w) for w in args.debug_tile.split(",")]
# Get the set of bounding boxes, after normalizing them
boxes = np.sort(np.array([d.box for d in mapinfo.datasets]),-2)
parser.add_argument("area")
parser.add_argument("odir")
parser.add_argument("-t", "--tsize", type=int, default=480)
parser.add_argument("-p", "--pad", type=int, default=240)
parser.add_argument("-C", "--ncomp", type=int, default=3)
parser.add_argument("-B", "--obeam", type=str, default=None)
parser.add_argument("-c", "--cont", action="store_true")
parser.add_argument("-v", "--verbose", action="store_true")
parser.add_argument("-W", "--wiener", action="store_true")
parser.add_argument("-m", "--mask", type=str, default=None)
parser.add_argument("--filter-mode", type=str, default="weight")
parser.add_argument("--cg-tol", type=float, default=1e-4)
parser.add_argument("--detrend", type=int, default=1)
args = parser.parse_args()
config = jointmap.read_config(args.config)
mapinfo = jointmap.Mapset(config, args.sel)
tsize = args.tsize # pixels
pad = args.pad # pixels
dtype = np.float64
ncomp = args.ncomp
comm = mpi.COMM_WORLD
utils.mkdir(args.odir)
# Get the set of bounding boxes, after normalizing them
boxes = np.sort(np.array([d.box for d in mapinfo.datasets]), -2)
# Read the cmb power spectrum, which is an effective noise
# component. T-only
#cl_path = os.path.join(os.path.dirname(args.config),config.cl_background)
#cl_bg = powspec.read_spectrum(cl_path)[0,0]