added CalcDispersion routine and FDTDDispersion2D routine for q3 scalar case (q1 and q2 are vectors)