open TELEMAC-MASCARET The mathematically superior suite of solvers

Dear All,

I am running a t2d-sisyphe internally coupled simulation in telemac v6p3.
The continued computation in t2d runs without problems, but if I want to run a continued computation coupled with sisyphe a couple of problems arise.

I make my restart file by exporting the last time-step of the previous run using FUDAA-PREPRO1.3-RC3.

My printout variables of the previous sisyphe computation are:
B,E,CS*,QSSUSP,TOB,*CONC,R,*ES,1A*,H

I checked previous posts to correct reported bugs:
#10874 add line in INIT_MIXTE.f
#7144 I added the updated INIT_TRANSPORT.f

Oddly, are the parameter names of the previous run wrong, in the sense that the new run cannot recognize them.
E.g. "LAYER CONCKNES" which should be "LAYER1 CONC KG/L "

I followed this post #3145 and used the nomvar_sisyphe.f to change the parameter names in the *.slf file(each variable= 32 characters, 16 for the name, the other 16 for the unit) unfortunately without success, see below error messages.

Please find below the ORIGINAL error message:
ITERATION 0 TIME: 0.0000 S
CALLING PRERES_TELEMAC2D
BACK FROM PRERES_TELEMAC2D
CALLING DESIMP
BACK FROM DESIMP
CALLING HPROPA
BACK FROM HPROPA
CALLING DE CHPCON
BACK FROM CHPCON
CALLING BILAN
BACK FROM BILAN
TELEMAC2D COUPLED WITH: SISYPHE
TELEMAC-2D: INTERNAL COUPLING WITH SISYPHE
PREMIER APPEL DE SISYPHE
PART 0 : INITIALISING SISYPHE
INBIEF (BIEF): NOT A VECTOR MACHINE (ACCORDING TO YOUR DATA)

LIQUID BOUNDARIES: 2
FONSTR : FRICTION COEFFICIENTS READ IN THE
GEOMETRY FILE
STRCHE (BIEF): NO MODIFICATION OF FRICTION

TITLE OF PREVIOUS COMPUTATION: grid

NAME: BED SHEAR STRESS UNIT:
NAME: BOTTOM UNIT: M
NAME: CONC MAS CL1 UNIT:
NAME: CONC MAS CL2 UNIT:
NAME: EVOLUTION UNIT:
NAME: FRACLAY 1 CL 1 UNIT:
NAME: FRACLAY 1 CL 2 UNIT:
NAME: FREE SURFACE UNIT: M
NAME: LAYER1 CONCKNES UNIT:
NAME: LAYER1 THICKNESS UNIT: M
NAME: LAYER2 CONC UNIT:
NAME: LAYER2 THICKNESS UNIT: M
NAME: QS SUSPENSION UNIT:
NAME: RIGID BED UNIT:
NAME: WATER DEPTH UNIT: M

SUITE_SERAFIN : READ OF RECORD 1

TIME OF RECORD: 435900.0 S

VARIABLE : BED SHEAR STRESS
UNKNOWN, IT WILL NOT BE KEPT

VARIABLE : CONC MAS CL1
UNKNOWN, IT WILL NOT BE KEPT

VARIABLE : CONC MAS CL2
UNKNOWN, IT WILL NOT BE KEPT

VARIABLE : EVOLUTION
UNKNOWN, IT WILL NOT BE KEPT

VARIABLE : LAYER1 CONCKNES
UNKNOWN, IT WILL NOT BE KEPT

VARIABLE : LAYER1 THICKNESSM
UNKNOWN, IT WILL NOT BE KEPT

VARIABLE : LAYER2 CONC
UNKNOWN, IT WILL NOT BE KEPT

VARIABLE : LAYER2 THICKNESSM
UNKNOWN, IT WILL NOT BE KEPT

VARIABLE : QS SUSPENSION
UNKNOWN, IT WILL NOT BE KEPT

VARIABLE : RIGID BED
UNKNOWN, IT WILL NOT BE KEPT

---

RESCUE : HYDRODYNAMIC FILE
PREVIOUS CALCULATION WITHOUT NON ERODABLE
BOTTOM

---

INITIALISATION EN CAS DE COUPLAGE : PART= 0
ERROR IN INIT-MIXTE:
THE SUM OF THICKNESS OF BED LAYERS IS DIFFERENT FROM ERODIBLE BED THICKNESS

PLANTE: PROGRAM STOPPED AFTER AN ERROR

AND the error message after replacing some of the variables with the right ones from nomvar_sisyphe.f, were interestingly some are recognized(e.g.:LAYER2 CONC KG/L) and some are not(e.g.LAYER1 CONC KG/L)

TITLE OF PREVIOUS COMPUTATION: grid

NAME: BED SHEAR STRESS UNIT: N/M2
NAME: BOTTOM UNIT: M
NAME: CONC MAS CL1 UNIT: G/L
NAME: CONC MAS CL2 UNIT: G/L
NAME: EVOLUTION UNIT: M
NAME: FRACLAY 1 CL 1 UNIT:
NAME: FRACLAY 1 CL 2 UNIT:
NAME: FREE SURFACE UNIT: M
NAME: LAYER1 CONC UNIT: KG/L
NAME: LAYER1 THICKNESS UNIT: M
NAME: LAYER2 CONC UNIT: KG/L
NAME: LAYER2 THICKNESS UNIT: M
NAME: QS SUSPENSION UNIT:
NAME: RIGID BED UNIT:
NAME: WATER DEPTH UNIT: M

SUITE_SERAFIN : READ OF RECORD 1

TIME OF RECORD: 435900.0 S

VARIABLE : BED SHEAR STRESSN/M2
IS IN THE FILE BUT WILL NOT BE READ

VARIABLE : LAYER1 CONC KG/L
UNKNOWN, IT WILL NOT BE KEPT

VARIABLE : LAYER1 THICKNESSM
UNKNOWN, IT WILL NOT BE KEPT

VARIABLE : LAYER2 CONC KG/L
IS IN THE FILE BUT WILL NOT BE READ

VARIABLE : LAYER2 THICKNESSM
UNKNOWN, IT WILL NOT BE KEPT

VARIABLE : QS SUSPENSION
UNKNOWN, IT WILL NOT BE KEPT

VARIABLE : RIGID BED
UNKNOWN, IT WILL NOT BE KEPT

---

RESCUE : HYDRODYNAMIC FILE
PREVIOUS CALCULATION WITHOUT NON ERODABLE
BOTTOM

---

INITIALISATION EN CAS DE COUPLAGE : PART= 0


Did anyone encounter such a problem already?
Is there a minimum set of print out variables of the previous run that computation continued with sisyphe needs?
Is there a better way to make a restart-file than using FUDAA?

Regards,
Christian

Dear Christian,
it seems to be a problem related with Fudaa. Did you try to launch a simulation "outside" the Fudaa environment?

I hope it helps,

Pablo

Dear Pablo,

Thank you for your fast reply.

I am not sure exactly what you mean by "launch a simulation "outside" the Fudaa environment"?

But, I only use FUDAA to make the PREVIOUS SEDIMENTOLOGICAL COMPUTATION FILE for the continued computation (*.slf or *.ser) otherwise my telemac runs(on a cluster) are not linked to FUDAA.

Another interesting thing is that also when looking at my results file in BLUE KENUE they have the "wrong names".(like LAYER1 CONCKNES)

Is there another option to make a PREVIOUS SEDIMENTOLOGICAL COMPUTATION FILE?

Best Regards,
Christian