Hello Fay,
I had a quick look to your files. You can use a SUPG option of 0;0 because the FREE SURFACE GRADIENT COMPATIBILITY will do the job of preventing inf-sup conditions, also the time step seems very small. The most important remark is that you use the incident wave feature. I was just thinking that this was obsolete and about to propose to remove it in the next versions... As a matter of fact it should be superseded by the Thompson boundary conditions, in which case you would have to give in subroutine bord.f the values of depths and velocities at the entrance, according to your incident wave, as we do for tidal conditions. The theory in both cases is about the same and based on characteristics, but Thompson happened to be better on cases like the one called "gouttedo".
Otherwise I have no experience with seiches caused by waves, I knew only seiches caused by wind or atmospheric pressure.
With best regards,
Jean-Michel Hervouet