Generating gamma random field with given covariance matrix

I have to generate multivariate gamma distributions with given positive-definite covariance matrix. Anyone can suggest me a method?

Thanks,

library(clusterGeneration)