A robust numerical method for a coupled system of singularly perturbed parabolic delay problems