combine output questions: how to avoid trimming and how to normalize by # of overlapping frames

astropy: 5.2.2
ccdproc: 2.4.0
numpy: 1.24.3es

I'm combining a dozen dithers of a 2x2 mosaic of 2x2k NIR detectors. It seemed quickest to use ccdproc machinery because each detector image in the mosaic has updated WCS using astrometry.net Astrometry.net confirms that the detectors have a slight skew to each other and optical changes etc.

My questions:
1. Unfortunately, after searching this forum etc, I cannot find how to prevent ccdproc.combine from trimming the registered stack down to the coverage of the first image. Instead I want the output to encompass the entire footprint of the input stack. Possible?
2. I also see no way written up on how to scale the result at each pixel by the number of images that overlapped there. 'average' in combine doesn't seem to do this and I don't know the values a priori so can't use combine's weights.

Everything else in ccdproc is working fine for me, thank you for this package.