PointingVectorDirectionFun(t, d, s, rho)
    t - direction of Poyting
    d - electric displacement direction
    s - unit vector along wavevector (k) - propagation direction
    rho - walkoff angle
    F - should be zero