High-precision pulsar timing is highly dependent on the precise and accurate modelling of any effects that can potentially impact the data. In particular, effects that contain stochastic elements contribute to some level of corruption and complexity in the analysis of pulsar-timing data. It has been shown that commonly used solar wind models do not accurately account for variability in the amplitude of the solar wind on both short and long timescales. In this study, we test and validate a new, cutting-edge solar wind modelling method included in the enterprise software suite (widely used for pulsar noise analysis) through extended simulations. We use it to investigate temporal variability in LOFAR data. Our model testing scheme in itself provides an invaluable asset for pulsar timing array (PTA) experiments. Since, improperly accounting for the solar wind signature in pulsar data can induce false-positive signals, it is of fundamental importance to include in any such investigations. We employed a Bayesian approach utilising a continuously varying Gaussian process to model the solar wind. It uses a spherical approximation that modulates the electron density. This method, which we refer to as a solar wind Gaussian process (SWGP), has been integrated into existing noise analysis software, specifically enterprise . Our Validation of this model was performed through simulations. We then conduct noise analysis on eight pulsars from the LOFAR dataset, with most pulsars having a time span of $ 11$ years encompassing one full solar activity cycle. Furthermore, we derived the electron densities from the dispersion measure values obtained by the SWGP model. Our analysis reveals a strong correlation between the electron density at 1 AU and the ecliptic latitude (ELAT) of the pulsar. Pulsars with $|ELAT|< 3^ circ $ exhibit significantly higher average electron densities. Furthermore, we observed distinct temporal patterns in the electron densities in different pulsars. In particular, pulsars within $|ELAT|< 3^ circ $ exhibit similar temporal variations, while the electron densities of those outside this range correlate with the solar activity cycle. Notably, some pulsars exhibit sensitivity to the solar wind up to $45^ circ $ away from the Sun in LOFAR data. The continuous variability in electron density offered in this model represents a substantial improvement over previous models, that assume a single value for piece-wise bins of time. This advancement holds promise for solar wind modelling in future International Pulsar Timing Array (IPTA) data combinations.