Interpersonal relationships impact our health and well-being, whereby interpersonal similarity in for example personality or experienced emotions is an important facilitator. Until now, this interpersonal similarity has mostly been examined using variable-centered (i.e., one variable at a time) and cross-sectional approaches. However, interpersonal similarity often concerns multiple variables simultaneously (e.g., multiple personality related behaviors or emotions) and can be expected to change over time and contexts. We therefore propose a workflow that combines existing methods to optimally capture dyadic similarity in a multivariate, timepoint-specific way, by computing and modeling momentary profile similarity values. We illustrate the different steps of this workflow by applying it to an existing longitudinal data on the multivariate emotional experience of romantic couples. We will show how to compute momentary profile similarity values with two different indices for these data and how to explore and model the differences in these values across time by relating them to time-invariant and time-varying predictors.