Wave extreme events can be understood as the combination of storm-intensity, directionality and intra-time distribution. However, the dependence structure among these factors is still unclear. A methodology has been developed to model wave-storms whose components are linked together. The model is composed by three parts: an intensity module, a wave directionality module, and a intra-time distribution module. In the Storm-intensity sub-model, generalized Pareto distributions and hierarchical Archimedean copulas have been used to characterize the storm energy, unitary energy, peak wave-period and duration. In the Directionality and the Intra-time sub-models, the wave direction (at the peak of the storm) and the storm growth-decay rates are linked to the variables from the intensity model, respectively. The model is applied to the Catalan coast (NW Mediterranean). The outcomes denote spatial patterns that coincide with the state of knowledge. The proposed methodology is able to provide boundary conditions for wave and near-shore studies, saving computational time and establishing the dependence of the proposed variables. Such synthetic storms reproduce the inter-variable co-dependence of the original data.