The HLLD Approximate Riemann Solver for Magnetospheric Simulation