The three‐dimensional (3D) turbulence code TOKAM3X‐EIRENE, coupling the 3D non‐isothermal version of TOKAM3X to the EIRENE Monte Carlo solver has been developed with the ability to simulate self‐consistently the interactions between large‐scale flows and turbulence both in limited and diverted plasmas, including recycling. This is especially important for diverted plasmas, where neutrals play a key role and where the recycling source is strongly dominant. The code package relies on the same interface as the Soledge2D‐EIRENE code, which retains state‐of‐the‐art plasma–wall interaction, as well as atomic and molecular physics. In this paper, we present the first results obtained in WEST divertor geometry, in laminar mode, with the aim of verifying the new code package. The divertor density regimes are recovered, and the code results are shown to be consistent with the results of the two‐point model, thus opening the way for turbulent simulations.