open TELEMAC-MASCARET The mathematically superior suite of solvers

Hi, all. I am using parameter estimation keyword to identify the friction value of the case. The case is a classic L-shaped flume with observations. One of the observations was add to the user_mesures.f file. The basic modelling (without parameter estimation keyword) was run successfully. However, when the parameter estimation keyword was actived, i got the error messages as list below.
In addition, I am not sure the code of the observation information (including time and value) that added in the user_mesures.f has a right style.
Thanks.







***********************************************************************
DIFFERENT NUMBER OF PARALLEL PROCESSORS:
DECLARED BEFORE (CASE OF COUPLING ?): 1
TELEMAC-2D : 0
VALUE 1 IS KEPT

********************************************
* LECDON: *
* AFTER CALLING DAMOCLES *
* CHECKING OF DATA READ *
* IN THE STEERING FILE *
********************************************

BINARY RESULTS FILE NECESSARY
FOR PARAMETER ESTIMATION

PLANTE: PROGRAM STOPPED AFTER AN ERROR
RETURNING EXIT CODE: 2

---

MPI_ABORT was invoked on rank 0 in communicator MPI_COMM_WORLD
with errorcode 2.

NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
You may or may not see output from other processes, depending on
exactly when Open MPI kills them.

---

Traceback (most recent call last):
File "../../scripts/python3/telemac2d.py", line 7, in <module>
main('telemac2d')
File "/home/**/opentelemac/v8p1r1/scripts/python3/runcode.py", line 272, in main
run_study(cas_file, code_name, options)
File "/home/**/opentelemac/v8p1r1/scripts/python3/execution/run_cas.py", line 157, in run_study
run_local_cas(my_study, options)
File "/home/**/opentelemac/v8p1r1/scripts/python3/execution/run_cas.py", line 65, in run_local_cas
my_study.run(options)
File "/home/**/opentelemac/v8p1r1/scripts/python3/execution/study.py", line 610, in run
self.run_local()
File "/home/**/opentelemac/v8p1r1/scripts/python3/execution/study.py", line 445, in run_local
run_code(self.run_cmd, self.sortie_file)
File "/home/**/opentelemac/v8p1r1/scripts/python3/execution/run.py", line 182, in run_code
raise TelemacException('Fail to run\n'+exe)
utils.exceptions.TelemacException: Fail to run
/usr/local/bin/mpiexec -wdir /home/**/opentelemac/v8p1r1/mycase/20191127_Ushape/test.cas_2020-12-24-10h39min16s -n 1 /home/**/opentelemac/v8p1r1/mycase/20191127_Ushape/test.cas_2020-12-24-10h39min16s/out_telemac2d
***************************************************************************

Hi,

As mentioned in the error message, the keyword "BINARY RESULTS FILE" is missing. To help you, you can inspire yourself with the test case "$HOMETEL/example/telemac2d/estimation"

Best regards,

Cedric

Thanks for your kind reply.
Inspired with the test case (telemac2d/estimation), several keywords were added and the parameter was estimated by using reference files. The reference file, using the friction of 60, was generated using the telemac2d system.
***********************************************************
/REFERENCE FILE : f2d.slf
BINARY RESULTS FILE : rbi_estimation.slf
BINARY RESULTS FILE FORMAT : 'SERAFIN'
FORMATTED RESULTS FILE : rfo_estimation.txt
************************************************************




However, observational information (time and value of depth) added to user_mesures.f file was a more practical way to estimation the friction parameter instead of using reference file.
For the reasons, I deactived the keyword of REFERENCE FILE and observed depth was added to the user_mesures.f file.
************************************************************
IF (ABS(TT-720.D0) .LT. 1.D-4) THEN
HD%R(2086) = 0.39D0
UD%R(2086) = 0.18D0
VD%R(2086) = 0.320D0
ALPHA1%R(2086) = T1%R(2086)
ALPHA2%R(2086) = T1%R(2086)
ALPHA3%R(2086) = T1%R(2086)
WRITE(LU,*) 'this is a test'
END IF
************************************************************


However, the way of using user_mesures.f was failed (as list below) at the first interation of the parameter estimation .
*************************************************************
END OF TIME LOOP

NO BND ELEMENTS IN A 2D MESH

READ_MESH_INFO: TITLE=
NUMBER OF ELEMENTS: 8175
NUMBER OF POINTS: 4392

SINGLE PRECISION FORMAT (R4)

TITLE OF PREVIOUS COMPUTATION:

NAME: VELOCITY U UNIT: M/S
NAME: VELOCITY V UNIT: M/S
NAME: WATER DEPTH UNIT: M
NAME: FREE SURFACE UNIT: M
NAME: BOTTOM UNIT: M

READ_DATASET : READ OF RECORD 2

TIME OF RECORD: 60.00000 S
MEASUREMENTS TO IMPLEMENT IN USER_MESURES

PLANTE: PROGRAM STOPPED AFTER AN ERROR
RETURNING EXIT CODE: 2

---

MPI_ABORT was invoked on rank 0 in communicator MPI_COMM_WORLD
with errorcode 2.

NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
You may or may not see output from other processes, depending on
exactly when Open MPI kills them.

---

**************************************************************************

Hello lairuixun,

If you delete or comment the 2 linesin user_mesures subroutine, your computation will go further.
The call to subroutine PLANTE will write the line "PLANTE: PROGRAM STOPPED AFTER AN ERROR" in the listing and then the line "RETURNING EXIT CODE: " with a number of error.
Then STOP will stop the program with error messages depending on your configuration.

Anyway, as recommended in the TELEMAC-2D user manual, "it is strongly recommended to work only in permanent mode" i.e. PARAMETER ESTIMATION : 'FRICTION, STEADY' not only 'FRICTION'. The estimation TELEMAC-2D example is built like this.

Hope this helps,

Chi-Tuan

Thanks for your kind reply. It works now.