Analysis and compression of six-dimensional gyrokinetic datasets using higher order singular value decomposition