Stochastic weather simulation, or weather generators (WGs), have gained a wide acceptance and been used for a variety of purposes, including climate change studies and the evaluation of climate variability and uncertainty effects. The two major challenges in WGs are improving the estimation of interannual variability and reducing overdispersion in the synthetic series of simulated weather. The objective of this work is to develop a WG model of daily rainfall, incorporating a covariable that accounts for interannual variability, and apply it in three climate regions (arid, Mediterranean, and temperate) of Chile. Precipitation occurrence was modeled using a two-stage, first-order Markov chain, whose parameters are fitted with a generalized lineal model (GLM) using a logistic function. This function considers monthly values of the observed Sea Surface Temperature Anomalies of the Region 3.4 of El Nino-Southern Oscillation (ENSO index) as a covariable. Precipitation intensity was simulated with a mixed exponential distribution, fitted using a maximum likelihood approach. The stochastic simulation shows that the application of the approach to Mediterranean and arid climates largely eliminates the overdispersion problem, resulting in a much improved interannual variability in the simulated values.