Astropy now enables FITS file compression for data exceeding three dimensions. ndfits need to be updated accordingly.