For self-consistent ICRF heating simulations, four codes have been coupled: VMEC[1] provides a fully shaped, anisotropic 3D MHD equilibrium (bi- Maxwellian distribution function), transferred to Boozer coordinates by TERPSICHORE[2]. The full-wave code LEMan[3,4] provides the IC wave field, power deposition and wave numbers (new anisotropic dielectric tensor and upshifted k). Finally, VENUS[5,6] computes the evolution of the distribution function due to ICRF heating and Coulomb collisions on the background thermal plasma (using Monte Carlo operators for Coulomb collisions and ICRH).