TELEMAC-MASCARET Forum: flow throught a culvert (1/1)

Welcome, Guest

TOPIC: flow throught a culvert

flow throught a culvert 7 years 7 months ago #26171

adepoi OFFLINE Fresh Boarder Posts: 2	Hi, I need to extract flow evolution throught a culvert. I try control sections but it's not the answer because the culvert is incharge and under downstream control of sea-level. Thanks Regards, AP
	The administrator has disabled public write access.

flow throught a culvert 7 years 7 months ago #26187

c.coulet OFFLINE Moderator Posts: 3722 Thank you received: 1031	Hi No direct solution as the culvert is not in the mesh... You have to print the discharge in the listing or adapt the program to send the value to a file Regards PS: please, update your profile
	Christophe The administrator has disabled public write access.

flow throught a culvert 7 years 7 months ago #26199

Svensmolders
OFFLINE
Senior Boarder
Posts: 105
Thank you received: 20

Hi AP,

you can use the private variables foreseen in the code to be used by the enduser.

You could add the following keywords to your steering file:

NUMBER OF PRIVATE ARRAYS =1
NAMES OF PRIVATE VARIABLES ='CULVERTS M3/S
! make sure you ask output for the private variable
VARIABLES FOR GRAPHIC PRINTOUTS =U,V,S,H,TA1,PRIVE1 '

and then add the following lines of code to the telemac2d.f fortran routine. You have to add the lines after the call for the BUSE.f subroutine. (actually the code is for telemac3d so you might need to change the names a little bit; I think PRIVE2D will just be PRIVE)
! this writes a zero to all nodes of the mesh
CALL OS('X=0 ',X=PRIVE2D)
! then write discharge data to the specific culvert point in your mesh.
DO I=1,NBUSE
PRIVE2D%ADR(1)%P%R(ENTBUS(I))=-DBUS%R(I)
PRIVE2D%ADR(1)%P%R(SORBUS(I))= DBUS%R(I)
ENDDO
! the variable CULVERTS will then be present in your results file and you can extract the time series from there.

Hope this helps

Sven

The administrator has disabled public write access.