Unsteady supersonic flow calculations for wing-body combinations using harmonic gradient method