Infragravity (IG) waves play an important role in the coastal zone affecting wave runup, overtopping, dune erosion and shipping. They originate from the sea-swell wave groups that force them either as bound IG waves or as free IG waves due to the breakpoint variation. The bound IG waves that are released within the surfzone (partly) reflect and are either refractively trapped or leaky where the latter can travel large distances and arrive at remote shorelines (Ardhuin et al., 2014, Rijnsdorp et al., 2021). As such the incident IG wave field at a particular coast consists of wave- group-forced (WGF) IG waves combining the bound and released IG waves (Herbers and Burton, 1997), and free infragravity (FIG) waves consisting of locally trapped and remote IG waves. To assess the potential impact of overtopping as well as dune erosion during extreme storm conditions sophisticated process models like XBeach (Roelvink et al., 2007) and SWASH (Zijlema et al., 2011) are often used. In addition to the sea-swell conditions and tide and surge levels these models require a boundary condition for the incident IG waves. Typically, this boundary condition is defined by the bound IG waves only using the equilibrium theory of Hasselmann (1962). This approach may lead to a significant over estimation in the presence of a sloping beach, but also an underestimation in the presence of incident FIG waves. Here we use an extended version of SWAN (Booij et al., 1999) to predict the total incident IG wave field during storm conditions that subsequently can be used as a boundary condition for these more detailed overtopping and dune erosion models. To verify the modeling approach, detailed comparisons with recently acquired observations at the Sand Engine of the directional free and forced IG wave field (Rutten et al., 2023 ) are used.