open TELEMAC-MASCARET The mathematically superior suite of solvers

Other point:
apparently your listing printout says that the variable SEDMIX is unknown and will not be kept.
Suggestion (and this is maybe also simpler than adding a variable to an existing result file):
- use a conventional previous computation file, without the SEDMIX variable
- add the SEDMIX variable in your geometry file
- read in the geometry file in your Fortran file using Telemac-2D

USE DECLARATIONS_TELEMAC2D, ONLY: T2D_FILES, T2DGEO
USE DECLARATIONS_SISYPHE , ONLY : IELMT, MESH
...
Clemens

Dear Clemens,
thank you for your suggestion. I also think that the SEDMIX is not seen or read.
One question to your suggestion, what you mean by reading the geometry file in the fortran file using Telemac2d? could you please a bit explain it.
I did something: telemac2d.py (my geometry) and i got this result:
Kind regards
Najib

Hi
SEDMIX is not read because it's not a native variable of Telemac nor Sisyphe.
In the next version, a new development will be available to allow user to declare it's own variable in the steering file, those variable will be given in the geometry file. Have a look at the topic http://www.opentelemac.org/index.php/kunena/16-telemac-2d/7803-infiltration#16971
But before this you have to manually adapt the code ...

Regards

dear Costas,
the problem still remains. I think the SEDMIX is not read or seen from the geometry file. I have to find out why.

Najib

try to add the SEDMIX variable into your geometry file where you have the BOTTOM variable also.
Then the procedure is similar as you already did beside that you read in the T2DGEO file.
in init_compo.f you have to add
USE DECLARATIONS_TELEMAC2D, ONLY: T2D_FILES, T2DGEO
...

CALL BIEF_ALLVEC(1,SEDMIX,'SEDMIX ',IELMT,1,2,MESH)
ALLOCATE(W(NPOIN),STAT=ERR)
NAME='SEDMIX '
CALL FIND_IN_SEL(SEDMIX,NAME,T2D_FILES(T2DGEO)%LU,
& 'SERAFIN ',W,OK)

...

Hope this works,
Clemens

Before all, use output redirection in order to capture the listing into a file, e.g.:This way you can post the whole output and not the final screen.

Set DEBUGGER=1 and post the full output. The error that Clemens pointed regarding SEDMIX not read when reading the previous computation is expected. SEDMIX is read directly from the init_conpo routine and not suite_serafin.

Costas

Dear All,
I implemented your suggestions and here are all files.Please tell me, if I'm mistaken somewhere.
Now the output file says that there is a (Segmentation fault-invalid memory reference).
thanking you in advance.
najib

Hi
in your code, you allocate W(npoin) but npoin is not define!
This probably explain the segmentation fault

Try to replace npoin by MESH%NPOIN

Regards

thank you,your suggestion solved the Segmentation fault problem and the old problem appears again.
here are the new output and init_compo.f files.
Najib

Hi
Let the keyword DEBUGGER = 1 (you qhould probably also declare it in the sisyphe steering file)

The error message probably comes from the subroutine init_avai.f
you could add it to your fortran file and try to add some print to understand which variable as got a wrong value and then leads to this error

Regards