A Three Dimensional Finite Difference Time Domain Sub-Gridding Method