Electromagnetic field components as integral expressions

I am trying to evaluate the electromagnetic field near the focus of high NA objective. The most common way to do it is to represent the electromagnetic fields in integral form , for instance -using Richards-Wolf diffraction theory.

Is it possible to use integrals in Ex, Ey, Ez input quantities of Port or Scattering Boundary Condition?
Is it possible to import 3D numerical data matrix for each of the Ex, Ey, Ez input quantities and to do some interpolation instead of integral ?