Hello Xinjie,
I strongly recommend you to carefully read the TELEMAC-3D user manual before posting on this forum.
As David say, the feature already exists in TELEMAC-3D, you only need to apply what is written in the manuals. please see the dedicated section for keyword BINARY BOUNDARY DATA FILE (including the cautions, but first the name of the keyword!).
If you use BINARY DATA FILE 2, as written in the TELEMAC-3D user manual, "the user has to handle their reading within the FORTRAN FILE" and you have not done anything...
TELEMAC can only do what you write in the steering file and Fortran file...
Anyway, please also read the recommendations written in the TELEMAC-3D user manual and TELEMAC-3D reference manual for every keyword you fill in your steering file, it may prevent you from bad uses and strange results, e.g. I noticed:
- COEFFICIENT FOR HORIZONTAL DIFFUSION OF VELOCITIES = 0.0001 whereas you use k-epsilon,
- you use a Strickler law to model bottom friction. With this law, friction coefficient is quite low! It is the default value if you use Nikuradse law (which is the recommended law in TELEMAC-3D, as told in the recommendations appendix).
Hope this helps,
Chi-Tuan