code: def generate_gaussian_profile(seeing_fwhm):
      
      FWHM_G = 2 * math.sqrt(2 * math.log(2))
      sigma = seeing_fwhm / FWHM_G
      amplitude = 1.0 / (2 * math.pi * sigma * sigma)
      seeing_model = Gaussian2D(amplitude=amplitude,
                                x_mean=0.0,
                                y_mean=0.0,
                                x_stddev=sigma,
                                y_stddev=sigma)
      return seeing_model
docs: Generate a normalized Gaussian profile from its FWHM