Understanding the behaviour of solid solutions over wide ranges of temperature and pressure remains a major challenge to both theory and experiment. Here we report a detailed exchange Monte Carlo study using a classical ionic model of the model perovskite parascandolaite-neighborite (K,Na)MgF3 solid solution and its end-members for temperatures in the range 300-1000 K and pressures from 0-8 GPa. Full account is taken of the local environment of the individual cations, clustering and thermal effects. Properties considered include the crystal structure, phase transitions, the thermodynamics of mixing and the non-ideality of the solid solution. Clustering of the potassium ions is examined via a short-range order parameter. Where experimental data are available for comparison, agreement is very good.