Double Proper Orthogonal decomposition for representing and simulating turbulence fields