Parallel implementation of a 3D subgridding FDTD algorithm for large simulations