We here developed a third-order accurate numerical method for scattering of 3D electromagnetic waves by doubly periodic structures. The method is an intuitively simple numerical scheme based on a boundary integral formulation. ...