An efficient parallel code for the approximate solution of initial boundary value problems for hyperbolic balance laws is introduced. The method combines three modern numerical techniques: locally-adaptive upwind finite-volume methods on unstructured grids, parallelization based on non-overlapping domain decomposition, and dynamic load balancing. Key ingredient is a hierarchical mesh in three space dimensions.The proposed method is applied to the equations of compressible magnetohydrodynamics (MHD). Results for several testproblems with computable exact solution and for a realistic astrophysical simulation are shown.