add harmonic source: the suburtine for generate the harmonic source needs modification and a subroutine to sort also needed