A novel model and simulation method for multivariate Gaussian fields involving nonlinear probabilistic dependencies and different variable-wise spatial variabilities