Hi Cath,
Thank you very much for the tip, I have no experience in programming and I am struggling to find my way around telemac routines.
I managed to sort this issue and the next one regarding the FFORMAT addition in FIND_IN_SEL for v6p3. Now the file compiles successfully but the program fails to run.
Debugger shows that it fails when the program calls the NOEROD subroutine. The error is:
LIT: READ ERROR
ONE INTENED TO READ A RECORD OF 1 VALUES OF TYPE: R8
ON LOGICAL UNIT: 1
This is the main routine for reading the non-erodible bed from the geometry file:
CALL BIEF_ALLVEC(1,NOER,'NOER ',IELMT,1,2,MESH)
ALLOCATE(W(NPOIN),STAT=ERR)
NAME='NOER '
CALL FIND_IN_SEL(NOER,NAME,SIS_FILES(SISGEO)%LU,'SERAFIND',W,OK)
CALL OV( 'X=Y-Z ',ZR,ZF,NOER%R,0.D0,NPOIN)
DEALLOCATE(W)
What am I doing wrong?
Many thanks in advance!
Costas