; applies Gaussian smoothing to a 2D array
pro beam_smoothing, INPUT=input, OUTPUT=output, SIGMA=sigma

  ; sigma is the standard deviation of the Gaussian kernel
  ; EDGE_WRAP treats the data as a periodic map
  ; (which is the case for our input data)
  output = gauss_smooth(input, sigma, /EDGE_WRAP)

  return

end