This study presents a comprehensive methodology for the reanalysis, characterization and propagation of wave bound infragravity waves (1/300Hz to 1/30Hz) in the nearshore. A deep water short wave reanalysis has been modified analytically to include the energy input associated with infragravity waves (full spectrum), through the decomposition and energy balance of the radiation stress of waves in deep water. A hybrid clustering-numerical wave propagation downscaling technique was used, forced by the modified full wave spectra previously obtained, resulting in hourly series of high spatial resolution nearshore ocean waves, including short and bound long waves. The methodology has been applied to the northern coast of Spain and validated by comparing buoy and pressure gauge records in the coast, for short wave (agitation) and infragravity wave (resonance) responses in open harbours, with very satisfactory results.