open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello, I'm having a problem using NOEROD.f

If I edit the rigid bed command line:
CALL OV( 'X=C ',ZR,ZF,ZF,-100.D0,NPOIN)

...to this do-loop:
DO I=1,NPOIN
IF(ZF(I).LT.2.D0) THEN
CALL OV( 'X=C ',ZR,ZF,ZF,-1.D0,I)
ELSE
CALL OV( 'X=C ',ZR,ZF,ZF,0.D0,I)
ENDIF
END DO

... I get 'divide by zero' or 'floating overflow' errors.

Any suggestions much appreciated.
Thanks,
Peter

Dear Peter,
this topic has been extensively discussed some weeks ago in the Sisyphe's post "Definition of non-erodable points".
You can also take a look to the validation test "005_rigid-bed"

All the best,

Pablo

Hello,

Please have a look at the guide for programming in the Telemac system, to understand what is subroutine OV and what it does. The loop from 1 to NPOIN is done by OV, NPOIN being the last argument. If you put I as last argument in a loop from 1 to NPOIN, imagine what it can do...

Regards,

JMH

Great,
Thanks both, working now. I modified ZR rather than subroutine OV
Peter