  A new version of GBS (Ricci et al. Plasma Phys. Control. Fusion 54, 124047,
2012; Halpern et al. J. Comput. Phys. 315, 388-408, 2016; Paruta et al. Phys.
Plasmas 25, 112301, 2018) is described. GBS is a three-dimensional,
flux-driven, global, two-fluid turbulence code developed for the
self-consistent simulation of plasma turbulence and kinetic neutral dynamics in
the tokamak boundary. In the new version presented here, the simulation domain
is extended to encompass the whole plasma volume, avoiding an artificial
boundary with the core, hence retaining the core-edge-SOL interplay. A new
toroidal coordinate system is introduced to increase the code flexibility,
allowing for the simulation of arbitrary magnetic configurations (e.g.
single-null, double-null and snowflake configurations), which can also be the
result of the equilibrium reconstruction of an experimental discharge. The
implementation of a new iterative solver for the Poisson and Amp\`ere equations
is presented, leading to a remarkable speed-up of the code with respect to the
use of direct solvers, therefore allowing for efficient electromagnetic
simulations that avoid the use of the Boussinesq approximation. The
self-consistent kinetic neutral model, initially developed for limited
configurations, is ported to the magnetic configurations considered by the
present version of GBS and carefully optimized. A new MPI parallelisation is
implemented to evolve the plasma and neutral models in parallel, thus improving
the code scalability. The numerical implementation of the plasma and neutral
models is verified by means of the method of manufactured solutions. As an
example of the simulation capabilities of the new version of GBS, a simulation
of a TCV tokamak discharge is presented.
