We use the Gutzwiller Monte Carlo approach to simulate the dissipative XYZ-model in the vicinity of a dissipative phase transition. This approach captures classical spatial correlations together with the full on-site quantum behavior, while neglecting non-local quantum effects. By considering finite two-dimensional lattices of various sizes, we identify a ferromagnetic and two paramagnetic phases, in agreement with earlier studies. The greatly reduced numerical complexity the Gutzwiller Monte Carlo approach facilitates efficient simulation of relatively large lattice sizes. The inclusion of the spatial correlations allows to describe critical behavior which is completely missed by the widely applied Gutzwiller decoupling of the density matrix.