3D simulations of surface harmonic generation with few-cycle laser pulses