open TELEMAC-MASCARET The mathematically superior suite of solvers

Hi, I want to modify subroutine BUSE.f to operate a culvert using a pre-defined start head and end head (i.e. hydro-turbine operation).

Does anyone have any experience with this kind of setup?


cheers
Ed

Hi Ed,

You may find the code here useful: github.com/NHanousek/T2D-TRS it's for simulating tidal range schemes in Telemac, and should have most of the features you need.

Let me know if you need any help with it

Nick

Hi Nick,

I am extremely thankful for the code. It is beyond what I could have written with the time and Fortran knowledge that I have.

cheers,Ed

Hi Ed,

Thanks, I hope it is useful!

Add all three files, to your 'user_fortran'. You need all three for it to work. LECBUS is the initial loader of the culvert tools, and T2D-TRS contains various functions that are called by BUSE (to try and keep the implementation abstracted where possible).

All the best,
Nick

Hi Nick,

When I include the subroutines in the user fortran file, I am faced with the following error:
(Segmentation fault: address not mapped to object at address 0xffffffffd4a3b57c)

then in the specific node LOG file...
PROPIN: BOUNDARY POINT 13 UNKNOWN PARAMETER FOR V

Without the custom fortran file (same setup) the model runs without error. I have checked my boundary conditions and tidal database, all seems ok.

Hi Ed,

Without seeing your setup this is hard to diagnose. Are you using your own mesh etc?

Can you test run with one of the sample models (I suggest the simple one first) then you will know if it is an issue with the Telemac version (I ran on v8p3) and other updates have caused something to change

Update: I now need to save a new variable inside the BUSE.f routine for each time step (LT) and each culvert(N).
Can I do something like this...
PRIVE%ADR(1)%N(LT) = 2

Hi
PRIVE as the dimension of the mesh (NPOIN)...
This will crash if LT =0 or if LT became greater than NPOIN

Regards