  In this work, we consider the general family of the so called ADER PNPM
schemes for the numerical solution of hyperbolic partial differential equations
with \textit{arbitrary} high order of accuracy in space and time.
  The family of one-step PNPM schemes was introduced in [Dumbser et al., JCP,
2008] and represents a unified framework for classical high order Finite Volume
(FV) schemes (N=0), the usual Discontinuous Galerkin (DG) methods (N=M), as
well as a new class of intermediate hybrid schemes for which a reconstruction
operator of degree M is applied over piecewise polynomial data of degree N with
M>N. In all cases with M >= N > 0 the PNPM schemes are linear in the sense of
Godunov, thus when considering phenomena characterized by discontinuities,
spurious oscillations may appear and even destroy the simulation. Therefore, in
this paper we present a new simple, robust and accurate a posteriori subcell
finite volume limiting strategy that is valid for the entire class of PNPM
schemes. The subcell FV limiter is activated only where it is needed, i.e. in
the neighborhood of shocks or other discontinuities, and is able to maintain
the resolution of the underlying high order PNPM schemes, due to the use of a
rather fine subgrid of 2N+1 subcells per space dimension.
  The paper contains a wide set of test cases for different hyperbolic PDE
systems, solved on adaptive Cartesian meshes (AMR) that show the capabilities
of the proposed method both on smooth and discontinuous problems, as well as
the broad range of its applicability. The tests range from compressible
gasdynamics over classical MHD to relativistic magnetohydrodynamics.
