Thanks Jean-Michel
I implemented your changes and found the segmentation error was caused by a silly mistake: I incorrectly defined time in the SELAFIN file as double precision instead of real. Please find the updated subroutine attached.
Unfortunately it still doesn't work: I have attached an image of the TELEMAC printout. condim_sisyphe.f correctly assigns the wave height, period and direction to the variables HRW, TWR and THETAWR as shown by the printout but then I get a MISSING WAVE DATA error message.
THETAWR, HWR and TWR are variable names in the subroutine call, however I did see that Pat used
USE DECLARATIONS_SISYPHE, ONLY: HW,TW,THETAW
and assigns values to these variables not THETAWR, HWR, TWR is this what I need to do? I don't understand, what is the difference between HWR and HW etc?
Thanks
Fay