An integrated model capable of self-consistent Ion Cyclotron Resonant Heating (ICRH) simulations has been developed. This model includes both full shaping and pressure effects, warm contributions to the dielectric tensor, pressure anisotropy and finite orbit width. It evolves the equilibrium, wave field and full hot particle distribution function until a self-consistent solution is found. This article describes the workings of the three codes VMEC, LEMan and VENUS and how they are linked for iterated computations in a code package we have named SCENIC. The package is thoroughly tested and it is demonstrated that a number of iterations have to be performed in order to find a consistent solution. Since the formulation of the problem can treat general 3D systems, we show a quasi-axisymmetric stellarator low power test case, and then concentrate on experimentally relevant Joint European Torus (JET) 2D configurations. (C) 2010 Elsevier B.V. All rights reserved.