This investigation focuses on hyper reduction of CFD data that depends on a series of coordinates (space, time) and/or parameters (flow speed, flow direction, etc.) into a series of basis vectors and generalized coordinates occupying much less storage. The data, first mapped on a multidimensional rectangular grid, can be “shuffled” into two-dimensional arrays on which Proper Orthogonal Decomposition (POD) is applied. This shuffling can be done in various ways but also repeated leading to a succession of POD-based data reduction, termed Progressive Proper Orthogonal Decomposition (p-POD). The possibility of downsampling is also explored using autoregressive (AR) modeling. These techniques are demonstrated on ship airwake data from two ships which depends on the three spatial coordinates, time, the Wind-Over-Deck angles, and wind speed. The p-POD approach leads to a data reduction by a factor 265 while maintaining 99% of the energy. Alone, the AR based approach allows a data reduction by a factor of 12.5 on its own and 297 with the p-POD. The stated reduction factors are computed with respect to the data on the Cartesian extraction grid. If the mapping on this grid is also included, the data reduction achieved from the CFD mesh maybe even much larger.