Hello,
Tomawac is a spectral wave model, it resolves the propagation/transformation of a wave spectrum on each grid node for each time step. Hs, Tp, Dir are only integrated parameters partially representing a spectrum.
So, you (or the model by calling speini.f) must build a spectrum from these parameters on each boundary node for the model to work.
For the difference between the 2 versions, you need to investigate what you've modified in limwac subroutine to understand what differs and why, by putting "print" statements for instance.
Regards,
Laurent