; return the first 4 moments of data
function get_moments, data

  mean = mean(data,/double)
  stddev = stddev(data,/double)
  skewness = skewness(data,/double)
  kurtosis = kurtosis(data,/double)

  ret = [mean, stddev, skewness, kurtosis]

  return, ret

end