High order parallel discontinuous Galerkin time domain method with curvilinear element and continuously varying material properties for Maxwell's equations