How to calculate the green function of n number of dipoles placed randomly in photonic structures

HI,
We are trying to calculate the green function of n number of dipole source embedded randomly in fcc photonic crystal. We are attaching one file have look at it and please tell us how to proceed further. We have gone through one example that is already given in Knowledge exchange pcfcc.fsp (338.7 KB)