open TELEMAC-MASCARET The mathematically superior suite of solvers

Dear JMH,
Thank you very much. It seems that the error occures in the
DO KK=1,NPOIN
UDEL%r(KK)=u%r(KK)
write(lu,*) udel%r(KK)
VDEL%r(KK)=v%r(KK)
enddo

When I add a 'write(lu,*) udel%r(KK)' like that, it shows that udel%r(KK)=0.0000

Hello,

I got it, NPOIN here is a local variable that is not initialised, not the number of points in the mesh transmited in the arguments. So add:

NPOIN=H%DIM1 before using it, and don't forget to add the forgotten DM1%R(KK)=0.D0 into a loop.

The CALL OS to initialise H,U and V are still not useful if H, U and V are given by your files...

With best regards,

JMH

Dear JMH,
Thank you very much.It shows the same error do as you said. But it can work when I modified the code as this:
DO KK=1,NPOIN
UCONV%R(KK)=u%r(KK)
VCONV%R(KK)=v%r(KK)
DM1%R(KK)=0.D0
enddo
! PRINT*,'EXITING PROPAG'
! CALL OS('X=0 ',X=DM1)
! DO I=1,3*MESH%NELMAX
! ZCONV%R(I)=0.D0
! ENDDO
Can I do like that?
But what is the meaning of DM1%R(KK) and ZCONV%R?

Hello,

The lines :

CALL OS('X=0 ',X=DM1)
DO I=1,3*MESH%NELMAX
ZCONV%R(I)=0.D0
ENDDO

should work as well so I am not sure that there is no other mistake. In some options however DM1 and ZCONV will not be used and this will trigger no mistake.

With best regards,

Jean-Michel Hervouet

Dear JMH,
I see. Thank you very much. YOu helped me a lot. Thank you. I will try to find the problem.

Hello
I am trying the same situation, Did you success to do it?

Thx a lot