On the basis of the Neyman-Scott process a model is proposed to generate rainfall at two sites simultaneously. This model takes into account the correlation, which is produced by generating bivariate correlated variates. The cross correlation between the sites has been computed in order to estimate the parameters. Finally, the model is validated for three pairs of stations situated in the Swiss Plateau region.