Superconducting (SC) tokamak JT-60SA plays an essential role in fusion research and development by supporting and complementing the ITER project, providing directions to the DEMO design activity and fostering next generation scientists and engineers. Since the short circuit incident at the terminal joints of equilibrium field coil #1 during the integrated commissioning (IC) in March 2021, both EU and JA implementing agencies (IAs) have examined how to ensure safe operation of JT-60SA by mitigating the risk of possible discharge occurrence inside the cryostat. Based on the experience of the global Paschen tests, the IAs have established a strategy of risk mitigation measures, which is a combination of (i) reinforcement of insulation, (ii) avoiding unnecessary voltage application to the coil systems and (iii) immediate de-energization of the coils when deteriorated vacuum conditions are detected. Thanks to the considerable efforts of the Integrated Project Team members, the IC restarted in May 2023. After confirmation of the SC state of the coil systems (TF, EF and CS), the coil energization test and the plasma operation phase 1 (OP-1) started. The first plasma was successfully achieved on 23 October 2023 with a limited value of voltage and current applied to the coils. The plasma configuration control was also confirmed with low plasma current and low auxiliary heating power conditions. Based on the IO-F4E-QST collaboration, activities of JT-60SA have been shared with the IO and provided an important lesson for ITER assembly and commissioning, and will provide an outstanding contribution to fusion research at large. After OP-1, maintenance & enhancement phase 1 (M/E-1) starts from January 2024, in which in-vessel components are installed, and heating and diagnostic systems are extensively upgraded to allow a high power heating experiment planned in OP-2. In order to make the best use of JT-60SA, a newly organized JT-60SA experiment team will refine the research plan for the future high heating power operation phase.

The mission of WEST (tungsten-W Environment in Steady-state Tokamak) is to explore long pulse operation in a full tungsten (W) environment for preparing next-step fusion devices (ITER and DEMO) with a focus on testing the ITER actively cooled W divertor in tokamak conditions. Following the successful completion of phase 1 (2016-2021), phase 2 started in December 2022 with the lower divertor made entirely of actively cooled ITER-grade tungsten mono-blocks. A boronization prior the first plasma attempt allowed for a smooth startup with the new divertor. Despite the reduced operating window due to tungsten, rapid progress has been made in long pulse operation, resulting in discharges with a pulse length of 100 s and an injected energy of around 300 MJ per discharge. Plasma startup studies were carried out with equatorial boron nitride limiters to compare them with tungsten limiters, while Ion Cyclotron Resonance Heating assisted startup was attempted. High fluence operation in attached regime, which was the main thrust of the first campaigns, already showed the progressive build up of deposits and appearance of dust, impacting the plasma operation as the plasma fluence increased. In total, the cumulated injected energy during the first campaigns reached 43 GJ and the cumulated plasma time exceeded 5 h. Demonstration of controlled X-Point Radiator regime is also reported, opening a promising route for investigating plasma exhaust and plasma-wall interaction issues in more detached regime. This paper summarises the lessons learned from the manufacturing and the first operation of the ITER-grade divertor, describing the progress achieved in optimising operation in a full W environment with a focus on long pulse operation and plasma wall interaction.

In 2021 JET exploited its unique capabilities to operate with T and D-T fuel with an ITER-like Be/W wall (JET-ILW). This second major JET D-T campaign (DTE2), after DTE1 in 1997, represented the culmination of a series of JET enhancements-new fusion diagnostics, new T injection capabilities, refurbishment of the T plant, increased auxiliary heating, in-vessel calibration of 14 MeV neutron yield monitors-as well as significant advances in plasma theory and modelling in the fusion community. DTE2 was complemented by a sequence of isotope physics campaigns encompassing operation in pure tritium at high T-NBI power. Carefully conducted for safe operation with tritium, the new T and D-T experiments used 1 kg of T (vs 100 g in DTE1), yielding the most fusion reactor relevant D-T plasmas to date and expanding our understanding of isotopes and D-T mixture physics. Furthermore, since the JET T and DTE2 campaigns occurred almost 25 years after the last major D-T tokamak experiment, it was also a strategic goal of the European fusion programme to refresh operational experience of a nuclear tokamak to prepare staff for ITER operation. The key physics results of the JET T and DTE2 experiments, carried out within the EUROfusion JET1 work package, are reported in this paper. Progress in the technological exploitation of JET D-T operations, development and validation of nuclear codes, neutronic tools and techniques for ITER operations carried out by EUROfusion (started within the Horizon 2020 Framework Programme and continuing under the Horizon Europe FP) are reported in (Litaudon et al Nucl. Fusion accepted), while JET experience on T and D-T operations is presented in (King et al Nucl. Fusion submitted).

The fusion reaction between deuterium and tritium, D(T,n)4 He is the main source of energy in future thermonuclear reactors. Alpha-particles (4 He-ions) born with an average energy of 3.5 MeV transferring energy to the thermal plasma during their slowing down, should provide the self-sustained D-T plasma burn. The adequate confinement of alpha-particles is essential to provide efficient heating of the bulk plasma and steady burning of a reactor plasma. That is why the fusion-born alpha-particle studies have been a priority task in the second D-T experiments (DTE2) on the Joint European Torus (JET) to understand the main mechanisms of their slowing down, redistribution and losses and to develop optimal plasma scenarios. JET with Be-wall and W-divertor, enhanced auxiliary heating systems and improved energetic-particle diagnostic capabilities, producing significant population of alpha-particles, provided the possibility for comprehensive studying of the alpha-particle behaviour. Selected results of the confined and lost alpha-particle measurements, evidence of alpha-particle self-heating and assessments of the fusion performance are presented in this paper giving an opportunity for further modelling and extrapolation to the International Thermonuclear Experimental Reactor and burning plasma reactors.

Within the 9th European Framework programme, since 2021 EUROfusion is operating five tokamaks under the auspices of a single Task Force called ‘Tokamak Exploitation’. The goal is to benefit from the complementary capabilities of each machine in a coordinated way and help in developing a scientific output scalable to future largre machines. The programme of this Task Force ensures that ASDEX Upgrade, MAST-U, TCV, WEST and JET (since 2022) work together to achieve the objectives of Missions 1 and 2 of the EUROfusion Roadmap: i) demonstrate plasma scenarios that increase the success margin of ITER and satisfy the requirements of DEMO and, ii) demonstrate an integrated approach that can handle the large power leaving ITER and DEMO plasmas. The Tokamak Exploitation task force has therefore organized experiments on these two missions with the goal to strengthen the physics and operational basis for the ITER baseline scenario and for exploiting the recent plasma exhaust enhancements in all four devices (PEX: Plasma EXhaust) for exploring the solution for handling heat and particle exhaust in ITER and develop the conceptual solutions for DEMO. The ITER Baseline scenario has been developed in a similar way in ASDEX Upgrade, TCV and JET. Key risks for ITER such as disruptions and run-aways have been also investigated in TCV, ASDEX Upgrade and JET. Experiments have explored successfully different divertor configurations (standard, super-X, snowflakes) in MAST-U and TCV and studied tungsten melting in WEST and ASDEX Upgrade. The input from the smaller devices to JET has also been proven successful to set-up novel control schemes on disruption avoidance and detachment.

The presented data refer to the Shattered Pellet Injector (SPI) experiments carried out at JET in 2019-2020. This paper is a full journal version of the data originally presented as posters at TMPDM_2020 and EPS_2021. This paper presents various aspects of the interaction of pellets with plasma and associated disruptions. The experiment was performed with I p = (1.1-3.1) MA plasmas and mainly with Ne + D2 pellet composition, but also with Ar pellets. The Current Quench (CQ) time, tau 80-20, is the key characteristic of mitigation effectiveness. A pellet with a high content of Ne or Ar can reduce the CQ duration below the upper required JET threshold. Plasmas with high (thermal + internal poloidal magnetic) pre-disruptive plasma energy require a high content of Ne pellets to obtain a short CQ duration. Pellets with a small amount of Ne (and accordingly large amount of D), instead of causing a mitigated CQ, create the conditions for a 'cold' Vertical Displacement Events (VDE). The SPI was applied to plasma with different status: mainly to normal ('healthy') plasma, i.e. not prone to disruption, post-disruptive and VDE plasma. This study shows that SPI effectiveness in terms of CQ duration and, accordingly, EM loads does not depend on the state of the plasma, whether it is 'healthy' or post-disruptive plasma. SPI has been shown to reduce the axisymmetric vertical vessel reaction forces by about (30-40) % compared to unmitigated disruptions. On JET, the VDE, whether 'hot' or 'cold', always creates the conditions for a toroidal asymmetry in the plasma, so the VDE on the JET is referred to as Asymmetric VDE (AVDE). The interrupting of VDE and prevention of AVDE with SPI has been demonstrated. Thus, the effectiveness of disruption mitigation using SPI has been confirmed.

This work reports on error field identification and control studies carried out at the Joint European Torus (JET). To identify the n=1 error field source, the novel compass scan method, designed and tested in the DIII-D device (Paz-Soldan C. et al, Nuclear Fusion 62 (2022) 126007), has been applied at JET. This method also yields the correction currents for error field compensation, which have been included in the empirical EF controller. The main outcomes of error field correction includes the spin-up of born locked modes and the ability to explore a lower density regime than otherwise achievable in the 1.8 MA plasma current, 2.1 T toroidal magnetic field Ohmic scenario.

The paper reports the first demonstration of in situ laser-induced desorption - quadrupole mass spectrometry (LID-QMS) application on a large scale fusion device performed in summer 2023. LID-QMS allows direct measurements of the fuel inventory of plasma facing components without retrieving them from the fusion device. The diagnostic desorbs the retained gases by heating a 3 mm diameter spot on the wall using a 1 ms long laser pulse and detects them by QMS. Thus, it can measure the gas content at any wall position accessible to the laser. The successful LID-QMS application in laboratory scale and on medium size fusion devices has now been demonstrated on the larger scale and it is already foreseen as tritium monitor diagnostic in ITER. This in situ diagnostic gives direct access to retention physics on a short timescale instead of campaign-integrated measurements and can assess the space-resolvedefficacy of detritation methods. LID-QMS can be applied on many materials: on Be deposits like in JET, B deposits like in TEXTOR, C based materials or on bulk-W.

For its initial operational phase, ITER has until recently considered using non-nuclear hydrogen (H) or helium (He) plasmas to keep nuclear activation at low levels. To this end, the Tokamak Exploitation Task Force of the EUROfusion Consortium carried out dedicated experimental campaigns in He on the ASDEX Upgrade (AUG) and JET tokamaks in 2022, with particular emphasis put on the ELMy H-mode operation and plasma-wall interaction processes as well as comparison to H or deuterium (D) plasmas. Both in pure He and mixed He + H plasmas, H-mode operation could be reached but more effort was needed to obtain a stable plasma scenario than in H or D. Even if the power threshold for the LH transition was lower in He, entering the type-I ELMy regime appeared to require equally much or even more heating power than in H. Suppression of ELMs by resonant magnetic perturbations was studied on AUG but was only possible in plasmas with a He content below 19%; the reason for this unexpected behaviour remains still unclear and various theoretical approaches are being pursued to properly understand the physics behind ELM suppression. The erosion rates of tungsten (W) plasma-facing components were an order of magnitude larger than what has been reported in hydrogenic plasmas, which can be attributed to the prominent role of He2+ ions in the plasma. For the first time, the formation of nanoscale structures (W fuzz) was unambiguously demonstrated in H-mode He plasmas on AUG. However, no direct evidence of fuzz creation on JET was obtained despite the main conditions for its occurrence being met. The reason could be a delicate balance between W erosion by ELMs, competition between the growth and annealing of the fuzz, and coverage of the surface with co-deposits.

Experiments on ASDEX Upgrade (AUG) in 2021 and 2022 have addressed a number of critical issues for ITER and EU DEMO. A major objective of the AUG programme is to shed light on the underlying physics of confinement, stability, and plasma exhaust in order to allow reliable extrapolation of results obtained on present day machines to these reactor-grade devices. Concerning pedestal physics, the mitigation of edge localised modes (ELMs) using resonant magnetic perturbations (RMPs) was found to be consistent with a reduction of the linear peeling-ballooning stability threshold due to the helical deformation of the plasma. Conversely, ELM suppression by RMPs is ascribed to an increased pedestal transport that keeps the plasma away from this boundary. Candidates for this increased transport are locally enhanced turbulence and a locked magnetic island in the pedestal. The enhanced D-alpha (EDA) and quasi-continuous exhaust (QCE) regimes have been established as promising ELM-free scenarios. Here, the pressure gradient at the foot of the H-mode pedestal is reduced by a quasi-coherent mode, consistent with violation of the high-n ballooning mode stability limit there. This is suggestive that the EDA and QCE regimes have a common underlying physics origin. In the area of transport physics, full radius models for both L- and H-modes have been developed. These models predict energy confinement in AUG better than the commonly used global scaling laws, representing a large step towards the goal of predictive capability. A new momentum transport analysis framework has been developed that provides access to the intrinsic torque in the plasma core. In the field of exhaust, the X-Point Radiator (XPR), a cold and dense plasma region on closed flux surfaces close to the X-point, was described by an analytical model that provides an understanding of its formation as well as its stability, i.e., the conditions under which it transitions into a deleterious MARFE with the potential to result in a disruptive termination. With the XPR close to the divertor target, a new detached divertor concept, the compact radiative divertor, was developed. Here, the exhaust power is radiated before reaching the target, allowing close proximity of the X-point to the target. No limitations by the shallow field line angle due to the large flux expansion were observed, and sufficient compression of neutral density was demonstrated. With respect to the pumping of non-recycling impurities, the divertor enrichment was found to mainly depend on the ionisation energy of the impurity under consideration. In the area of MHD physics, analysis of the hot plasma core motion in sawtooth crashes showed good agreement with nonlinear 2-fluid simulations. This indicates that the fast reconnection observed in these events is adequately described including the pressure gradient and the electron inertia in the parallel Ohm’s law. Concerning disruption physics, a shattered pellet injection system was installed in collaboration with the ITER International Organisation. Thanks to the ability to vary the shard size distribution independently of the injection velocity, as well as its impurity admixture, it was possible to tailor the current quench rate, which is an important requirement for future large devices such as ITER. Progress was also made modelling the force reduction of VDEs induced by massive gas injection on AUG. The H-mode density limit was characterised in terms of safe operational space with a newly developed active feedback control method that allowed the stability boundary to be probed several times within a single discharge without inducing a disruptive termination. Regarding integrated operation scenarios, the role of density peaking in the confinement of the ITER baseline scenario (high plasma current) was clarified. The usual energy confinement scaling ITER98( p,y ) does not capture this effect, but the more recent H20 scaling does, highlighting again the importance of developing adequate physics based models. Advanced tokamak scenarios, aiming at large non-inductive current fraction due to non-standard profiles of the safety factor in combination with high normalised plasma pressure were studied with a focus on their access conditions. A method to guide the approach of the targeted safety factor profiles was developed, and the conditions for achieving good confinement were clarified. Based on this, two types of advanced scenarios (‘hybrid’ and ‘elevated’ q -profile) were established on AUG and characterised concerning their plasma performance.

This work reports on error field identification and control studies carried out at JET and insights on the development of the empirical EF controller for ITER. The empirical EF controller has been included in the JET real-time central controller following the execution of the non-disruptive compass scan method (Paz-Soldan C. et al., Nuclear Fusion 62 (2022) 126007), which allowed the identification of the EF source and the currents for error field compensation. When testing the empirical EF controller, born locked n = 1 modes have been observed to spin-up and a lower density regime has been explored in the 1.8 MA plasma current, 2.1 T toroidal magnetic field scenario than otherwise achievable. These experimental results demonstrate the benefits of EF correction. In preparation of EF correction studies in ITER, the empirical EF controller for ITER has been developed and integrated in the plasma control system database.

In magnetic fusion devices, undesired non-axisymmetric magnetic field perturbations, typically called error fields (EF), have been observed to have a detrimental effect on plasma stability and confinement and can lead to brute plasma terminations, i.e. plasma disruption events. The main strategies that can be adopted to minimize the effect of EFs on the plasma dynamics consist in a careful alignment of the coils, when assembling the fusion device, and, most commonly, in the use of EF correction coils which counteract the non-axisymmetric fields by prescribing properly designed correction currents. In this work, an assessment of the n=1 EF source in MAST-U is presented. When constructing the MAST-U device, an optimization of poloidal and divertor coil positions has been adopted to reduce the n=1 EF source. This optimization consisted in the application of coil shifts and tilts, of the order of mms and mrads, respectively. To investigate the presence of a residual n=1 EF dedicated studies have been performed during the first MAST-U campaign. The compass scan method was employed to identify the n=1 EF, relying on the detection of the locked mode (LM) onset, which proved to be a challenging task. Therefore, a methodology based on Transfer Functions (TFs) among each coil and the n=1 radial magnetic field has been developed which allows the detection of LM formation. Such method is described here, complemented with the experimental results achieved, which suggest that the intrinsic n=1 EF source on MAST-U is relatively small with respect to MAST. Indeed, the empirical correction currents for n=1 EF minimization are smaller, about 0.2 kA, than the ones used in MAST, 1 kA range. This proves that the optimal coil alignment for n=1 EF minimization has been a successful strategy in MAST-U.

Precise values for radiated energy in tokamak disruption experiments are needed to validate disruption mitigation techniques for burning plasma tokamaks like ITER and SPARC. Control room analysis of radiated power (P rad) on JET assumes axisymmetry, since fitting 3D radiation structures with limited bolometry coverage is an under-determined problem. In mitigated disruptions, radiation is toroidally asymmetric and 3D, due to fast-growing 3D MHD modes and localized impurity sources. To address this problem, Emis3D adopts a physics motivated forward modeling ('guess and check') approach, comparing experimental bolometry data to synthetic data from user-defined radiation structures. Synthetic structures are observed with the Cherab modeling framework and a best fit chosen using a reduced chi 2 statistic. 2D tomographic inversion models are tested, as well as helical flux tubes and 3D MHD simulated structures from JOREK. Two nominally identical pure neon shattered pellet injection (SPI) mitigated discharges in JET are analyzed. 2D tomographic inversions with added toroidal freedom are the best fits in the thermal quench (TQ) and current quench (CQ). In the pre-TQ, 2D reconstructions are statistically the best fits, but are likely over-optimized and do not capture the 3D radiation structure seen in fast camera images. The next-best pre-TQ fits are helical structures that extend towards the high-field side, consistent with an impurity flow under the magnetic nozzle effect also observed in JOREK simulations. Whole-disruption radiated fractions of 0.98+0.03/ -0.29 and 1.01+0.02/-0.17 are found, suggesting that the stored energy may have been fully mitigated by each SPI, although mitigation efficiencies well below ITER and SPARC requirements for high energy pulses are still within the large uncertainties. Emis3D is also used to validate JOREK SPI simulations, and confirms improvements in matching experiment from changes to impurity modeling. Time-dependent toroidal peaking factors are calculated and discussed.

JET, the world’s largest operating tokamak with unique Be/W wall and tritium handling capability, completed a Deuterium-Tritium (D-T) campaign in 2021 (Maggi et al 29th Fusion Energy Conf. ) following a decade of preparatory experiments, dedicated enhancements, technical rehearsals and training (Horton et al 2016 Fusion Eng. Des. 109–111 925). Operation with tritium raises significant technical, safety and scientific challenges not encountered in standard protium or deuterium operation. This contribution describes the tritium operational requirements, pulses and technical preparations, new operating procedures, lessons learned and details on the achieved operational availability and performance. The preparation and execution of the recent JET tritium experiments benefitted from the previous experience in 1991 (Preliminary Tritium Experiment), 1997 (DTE1 campaign) and 2003 (Trace Tritium Campaigns) and consisted of the following five phases: technical rehearsals and scenario preparation, tritium commissioning, 100% tritium campaign, D-T campaign (DTE2), tritium clean-up. Following the clean-up JET resumed normal operation and is currently undertaking a further D-T campaign (DTE3).

At the beginning of a machine operation, an assessment of the intrinsic error fields, spurious magnetic field perturbations which can affect plasma dynamics, is often carried out by executing the compass scan method [Scoville J.T. et al. Nucl. Fusion 43 250 (2003)]. This method relies on the application of 3D magnetic fields with various phases, induced by EF correction coils, to trigger a locked mode. The instant of locked mode onset allows the identification of the amplitude and phase of the intrinsic error field, from which the empirical correction currents for its minimization can be deduced. The presence of a locked mode needs to be carefully monitored during this study because of the potential disruptive mode behavior, especially in devices which can tolerate a maximum number of disruptions, as in SPARC and in ITER. A novel method, the so-called non-disruptive compass scan method [Paz-Soldan C. et al., Nuclear Fusion 54 (2014) 073013], avoids the disruption risk, as the name recalls, via magnetic island healing, i.e. stabilizing the locked mode. The magnetic island healing is achieved by switching off the error field correction coil current during the execution of the compass scan and asynchronously by increasing the plasma density. The crucial point of this new method is the detection of the locked mode to initiate the EFCC-density control actions. In this work, the locked mode detector adopted during non-disruptive compass scan test at JET is presented, together with brand-new locked mode metrics, which take into account the actual poloidal deformation due to a locked mode and a class of MHD instabilities, named Beta Alfve & PRIME;n Eigenmodes, that appear in the Mirnov signal in concomitance to the locked mode. The use of multiple metrics for locked mode detection during the execution of the non-disruptive compass scan increases the fidelity of the real-time control system to pinpoint the event, compensating possible magnetic probe failure, and initiate the control sequences to heal the magnetic island.

The recent deuterium–tritium campaign in JET-ILW (DTE2) has provided a unique opportunity to study the isotope dependence of the L-H power threshold in an ITER-like wall environment (Be wall and W divertor). Here we present results from dedicated L-H transition experiments at JET-ILW, documenting the power threshold in tritium and deuterium–tritium plasmas, comparing them with the matching deuterium and hydrogen datasets. From earlier experiments in JET-ILW it is known that as plasma isotopic composition changes from deuterium, through varying deuterium/hydrogen concentrations, to pure hydrogen, the value of the line averaged density at which the threshold is minimum, nˉe,min , increases, leading us to expect that nˉe,min (T) < nˉe,min (DT) < nˉe,min (D) < nˉe,min (H). The new power threshold data confirms these expectations in most cases, with the corresponding ordering of the minimum power thresholds. We present a comparison of this data to power threshold scalings, used for extrapolation to future devices such as ITER and DEMO.

This work describes the usage of Error Field Correction Coil (EFCC) system [Barlow I. et al. 2001 Fusion Eng. Des. 58-59 189], which is a set of 4 coils located external to the vessel of the JET device, with the aim of introducing non-axisymmetric n =1 magnetic field perturbations in various targeted plasma experiments. Besides being used to characterize and correct the intrinsic error field, the EFCCs system has been exploited to i) probe plasma stability in high-ON regimes through the MHD spectroscopy technique, allowing the identification of resistivewall mode stable operational scenarios, ii) study radiation asymmetry when testing the shattered pellet injector as the prime candidate for ITER disruption mitigation system, and iii) control edge localized modes, permitting a reduction in heat fluxes and carbon erosion of the divertor target plates.

The isotope effect on intrinsic rotation was studied at the Joint European Torus (JET) tokamak. With the unique capability of JET to operate with tritium (T), for the first time, experiments in hydrogen (H), deuterium (D) and T in Ohmic plasmas were compared. Two rotation reversals per isotope type are observed in plasma density scans spanning the linear and the saturated Ohmic confinement regimes. A clear isotope mass dependence is observed at the higher densities. The magnitude of the core rotation was found to depend on isotope mass, with stronger co-current rotation observed in H. Change on intrinsic rotation characteristics coexist with a stronger thermal energy confinement in T.

The required heating power, P L H , to access the high confinement regime (H-mode) in tritium containing plasmas is investigated in JET with ITER-like wall at a toroidal magnetic field of B t = 1.8 T and a plasma current of I p = 1.7 MA. P L H , also referred to as the L-H power threshold, is determined in plasmas of pure tritium as well as mixtures of hydrogen with tritium (H-T) and mixtures of deuterium with tritium (D-T), and is compared to the L-H power threshold in plasmas of pure hydrogen and pure deuterium. It is found that, for otherwise constant parameters, P L H is not the same in plasmas with the same effective isotope mass, A e f f , when they differ in their isotope composition. Thus, A e f f is not sufficient to describe the isotope effect of P L H in a consistent manner for all considered isotopes and isotope mixtures. The electron temperature profiles measured at the L-H transition in the outer half of the radius are very similar for all isotopes and isotope mixtures, despite the fact that the L-H power threshold varies by a factor of about six. This finding, together with the observation of an offset linear relation between the L-H power threshold, P L H , and an effective heat diffusivity, χ e f f , indicates that the composition-dependent heat transport in the low confinement mode (L-mode) determines, how much power is needed to reach the necessary electron temperatures at the edge, and hence P L H .

Plasma free boundary simulation codes including light transport and current diffusion are an important mean to develop plasma scenarios and control algorithms in view of the operation of large tokamaks with exploratory scientific and technological purposes like ITER and JT-60SA. In particular, the intrinsic flexibility of the latter will offer the possibility to study a variety of scenarios including advanced scenarios which could be extrapolated to a steady-state demonstration of fusion power plants. The development of a numerical simulator is the object of this paper. The ultimate objective of producing a reasonably fast "flight simulator"is twofold: provide to scientists a tool for integrated scenario design including control laws and supply to tokamak operators an instrument for pulse discharge training. The option explored is obtained by coupling METIS with CREATE-NL+.