An operator splitting Legendre-tau spectral method for Maxwell's equations with nonlinear conductivity in two dimensions