Resolution of Gibbs phenomenon using a modified pseudo-spectral operator in harmonic balance CFD solvers*