Hi,
If you are using the coadded images
then those images are already the sum of many individual images with different exposure times and noise.
The images have a pixel-by-pixel uncertainty map, stored as inverse-variances, in the "invvar" files.
For example, each pixel in the file
legacysurvey-0001p000-image-r.fits.fz
has an uncertainty sigma = 1/sqrt(invvar), where the invvars are in this file,
legacysurvey-0001p000-invvar-r.fits.fz
So if you want to stack images, you can do inverse-variance weighted stacks.
cheers,
dustin