Hello,
if you know the coordinates and the wind speed values of your wind data, you can use the FASP subroutine to interpolate your wind data over the mesh: for example you can adapt and use the subroutine VENUTI.f.
If your wind data are unsteady, please refer to the subroutines LECDOI.f and NOUDON.f, for the temporal interpolation of the wind data. If the wind data are stationary, refer to the subroutine LECDON.f.
Another method is to convert the wind data you have in one of the formats that are automtically read by TOMAWAC (see Appendix 8 of the TOMAWAC manual). In this case you don't have to write or modify any TOMAWAC subroutine.
Hope it helps.
Kind regards,
Giovanni