Add command 'save array level FITS'

Generate FITS images of the rate matrix inverted to obtain level populations for the specified species, iteration, and zone. This command is a debugging aid.