Over the past decades, multiple gyrokinetic codes have shown to be able to simulate turbulence and associated transport in the core of Tokamak devices. However, their application to the edge and scrape-off layer (SOL) region presents significant challenges. To date, only few codes and models have been adapted to SOL/edge conditions. To further study the SOL region in particular, with its steep temperature and density gradients as well as large fluctuation amplitudes, the full-f particle-in-cell code PICLS has been developed. PICLS is based on a full-f gyrokinetic model with linearized field equations, considers kinetic electrons and uses logical sheath boundary conditions. In the past, PICLS was verified by applying it to a well-studied 1D parallel transport problem during an edge-localized mode in the SOL under both collisionless and collisional conditions, for which a Lenard–Bernstein collision operator was implemented. PICLS recently was extended towards three spatial dimensions to study turbulence in open-field-line regions in slab and closed-field-line toroidal geometries. In this work, we will focus on the models and methods we used for extending the code towards three spatial dimensions, including validation efforts and comparisons with other existing codes in closed-field-line geometry.