open TELEMAC-MASCARET The mathematically superior suite of solvers

Thank you for reviewing my case file,
We found that the Fortran files just read salinity (Tracer), not H, U and V.
I am not familiar with Fortran; I wish I could fix the .f files and post them here.
Thank you again,
Cheers

Hi DHverson!
I've been having a bit of a hard time lately, but luckily I came across your post. The binary file was also generated by generate_ban python script. I kept encountering errors when I tried to read it using user_bord3d.f subroutine. I am very happy to know that you have successfully run the message. I wonder if you could share the user_bord3d.f subroutine and read_bin_frliq.f subroutine that you have successfully run in telemac3d before with me for learning. I will be looking forward to your reply!
Best regard
Xinjie

Hi Xinjie

If you're running in 3D then you should not need to modify any code. The standard read_bin_frliq.f and bord3d.f are already set to read binary boundary files. You should not need to use the user_bord3d.f subroutine.

For example lines 258-270 in bord3d.f:

! READ ELEVATIONS GIVEN IN THE BINARY BOUNDARY DATA FILE IF PRESENT
!

---

IF(T3D_FILES(T3DBND)%NAME(1:1).NE.' ') THEN
! GET H FROM THE BINARY BOUNDARY DATA FILE
! IT IS NECESSARILY IN THE FILE OTHERWISE
! READ_BIN_FRLIQ CALLS PLANTE
CALL READ_BIN_FRLIQ
& (T3_01%R,'ELEVATION Z ',
& AT,T3D_FILES(T3DBND)%LU,T3D_FILES(T3DBND)%FMT,
& READ_BIN_Z)
ELSE
READ_BIN_Z = .FALSE.
ENDIF

One thing to double check though is the timestamp of your binary boundary file. If you've used converter.py generate_bnd script then there's a bug that doesn't reproduce the correct time stamp.

Hope this points you in the right direction.

Cheers

David

Hi DHaverson
Thank you very much for your reply!
I have tried without any subroutine modification and still failed (if I remember correctly), but I will try your suggestion again now and hope I have a good result.
Best regard!
Xinjie

Hi DHaverson
I ran the case again (using the default subroutine) but I didn't get a good result.
And the error was something I've never seen before, and I don't know if you've seen anything like it.
Thank you for your busy reply.
Xinjie

Hi all
Telemac was able to read the binary boundary condition after I regenerated the new binary boundary file, but worse was the error in reading the data. I will put the errors in the picture below, hoping that someone can give me some help, and I will continue to find the cause of these errors.
Xinjie

No error in the screen capture joined.
Just an indication that the RHS of the system was 0 then X (the unknown) is 0...

Hi coulet
Yes, I think my statement is ambiguous, although the program can run, but the final result is still not correct. Including INITIAL MASS OF WATER IN THE DOMAIN: 0 This is very unusual in the calculation, such as X=0 I suspect that the value in the binary boundary file may not be correctly read (if not, please correct). Finally, I attach some parameters of my binary boundary file and its vector plot in tecplot.

You could redo the scan with --data option to check the values...
Initial mass = 0 means no water so maybe a problem in the bathymetry of geo file

Hi coulet
After scaning the file with - -data, I found that the data in the file did not seem to be abnormal.
Xinjie