Python code to create diffuse gas that follows a 2-D beta (Lorentz) profile and write to a FITS file.