Welcome, Guest
Username: Password: Remember me

TOPIC: Modifying noerod.f

Modifying noerod.f 9 years 10 months ago #15453

  • amanj2013
  • amanj2013's Avatar
  • OFFLINE
  • Expert Boarder
  • Posts: 211
  • Thank you received: 24
Hello everyone,

I may need your help with solving this issue. I am trying to simulate dam breaching in Telemac2D and sisyphe. I have a Case as shown in the attached image. I have zero elevation which is flume level and its non erodible and I have an embankment which is erodible. How to modified the noerod.f in order to set the zero elevation is rigid and above zero is erodible.

Thanks
AMANJ
Attachments:
The administrator has disabled public write access.

Modifying noerod.f 9 years 10 months ago #15454

  • cyamin
  • cyamin's Avatar
  • OFFLINE
  • openTELEMAC Guru
  • Posts: 997
  • Thank you received: 234
Hello Amanj,

Your case is quite easy to implement. noerod.f by default sets the non-erodible bed elevation as 100m below bottom level:
!     DEFAULT VALUE: ZR=ZF-100.D0
!
      CALL OV('X=Y+C   ',ZR,ZF,ZF,-100.D0,NPOIN)
You want to set the non-erodible bed at absolute level 0m, so you could make the following modification:
      CALL OV('X=C     ',ZR,ZF,ZF,0.D0,NPOIN)
You can find more details about subroutine OV in the "telemac guide for programming".

Edit:
I noticed that your leeward side is lower than 0, so you need to add an IF statement to set a different non-erodible bed level after the dam (e.g. when X>500).

Best Regards,
Costas
The administrator has disabled public write access.
The following user(s) said Thank You: amanj2013

Modifying noerod.f 9 years 10 months ago #15469

  • amanj2013
  • amanj2013's Avatar
  • OFFLINE
  • Expert Boarder
  • Posts: 211
  • Thank you received: 24
Thank you dear Costas,

I fixed and now is working well, but I have another problem, when i have rigid bed My outflow discharge which is computing through control section is much lower than that case of erodable bed even 10 times and the peak value is too small in non erodable while its unexpected because peak reach before reaching to zero bed.
The administrator has disabled public write access.
Moderators: Pablo, pavans

The open TELEMAC-MASCARET template for Joomla!2.5, the HTML 4 version.