Model reduction by proper orthogonal decomposition for lambda-omega systems