add density to fdtd