Non-linear MHD simulations play an essential role in active research and understanding of tokamakplasmas for the realisation of a fusion power plant. The development of MHD codes like JOREK isa key aspect of this research effort. In this paper, we present a fully operational full-MHD modelin JOREK, a significant advancement from the reduced-MHD model used for previous studies.The model is presented in detail, and benchmarks are performed using both linear and non-linearsimulations, including comparisons between the new full-MHD model of JOREK and the previ-ously extensively studied reduced-MHD model, as well as results from the linear full-MHD codeCASTOR3D. It is shown that this new JOREK full-MHD model is numerically and physicallyreliable, even without the use of numerical stabilisation methods. Non-linear modelling results oftypical tokamak instabilities are presented, including disruption and ELM physics, most relevantto current open issues concerning future tokamaks like ITER and DEMO.