Thanks for your suggestion..

I am doing the same thing that you had stated. My question is that when i am calculating the effect of array by introducing floquet boundary condition. Then in that case what would be my integration surface, is it be the surface of nanoparticle or any surface which lie in the the far field domain??

]]>1. Calculate scattered energy by integration of the Poynting vector over an imaginary surface S around the array:

W_sca=intop1_S(-(emw.relPoavx*nx+ emw.relPoavy*ny+ emw.relPoavz*nz)) [W]

2. Calculate incident irradiance, defined as energy flux of the incident (background) wave:

P0=c0*epsilon*E0^2/2, [W/m^2]

where c0 is speed of light, epsilon is permittivity, E0 is background incident plane wave amplitude

3. Calculate scattering cross section:

sig_sca=W_sca/P0 [m^2]

]]>