The inverse distribution function method for drawing randomly from normal andtruncated normal distributions is used to set up a Gibbs' sampler for the posterior densityfunction of the parameters of a truncated multivariate normal distribution. The sampler isapplied to shire level rainfall for five shires in Western Australia.