Adaptive reconstruction of multiple source intensities for parabolic PDEs