A well-known, benchmarked code, where a Bohm–gyro–Bohm transport model is complemented with an empirical scaling for the dynamics of Internal Transport Barriers (ITBs), is used to model the ITB oscillations often seen in advanced tokamak scenarios with a dominant fraction of bootstrap current.