Describir: Estimation and use of a multivariate parametric model for simulating heteroscedastic, correlated, non-normal random variables: the case of corn-belt corn, soybeans and wheat yields