open TELEMAC-MASCARET The mathematically superior suite of solvers

Dear all ;
I use telemac 2D version v6p2, internal coupling with sisyphe. In my steering file (for telemac 2d) I have used BOUNDARY CONDITIONS FILE =QH.qsl. In which I have used the key word SL (m) for elevation, Q (m3/s) for discharge, TR for tracer (written at top of text file). Telemac 2D work good.
My question are:
- What is the key word for bed-load discharge ?
- I can write this data of this bed-load discharge in the same above file QS.qsl?
Thank you in advance for your help and advice.
Best regards
Thong NGUYEN
HCMUT, Viet Nam

Dear Thong,
thanks for your post. At this moment in Sisyphe version 6.2, bedload discharge can be specified via: (i) the keyword PRESCRIBED SOLID DISCHARGES (in which you can prescribe a constant solid discharge expressed in m3/s); (ii) modifying the boundary conditions file for Sisyphe (in which you can prescribed a nodal solid discharge expressed in m2/s, where the values can vary from one node to another); and (iii) by modifying the subroutine conlit.f.
Providing a time-varying bedload discharge from a file has not been yet implemented.

To clarify all these points, we are preparing a documentation that will be available soon. In the meantime, I would suggest you to take a look to some of the previous forum post, where this topic has been aborded extensively.

See, e.g. www.opentelemac.org/index.php/kunena/17-...nflows-and-rigid-bed

www.opentelemac.org/index.php/kunena/17-...tration-imposed#6899

All the best,

Pablo

Hi,

It seems imposing time-varying bedload input can only be done through modification of conlit.f

There is some discussion here:
www.opentelemac.org/index.php/kunena/17-...n-parallel-mode#3607

Does anyone have an example of a modified conlit.f, ideally where QBOR is loaded in from a specified text file.

Should be reasonably straightforward (?), but don't want to re-invent the wheel if someone already has a successful implementation.

Thanks,
Dougal

Dear Pablo and Dougal ;
Thank you very much for your suggestion and advice. My program telec2D internal coupling with sisyphe work with option where I have modified in subroutine conlit.f in the loop for all boundary node as:
- choice of number of boundary by used variable : NUMLIQ%I(K) (already know when I prepare my network with Blue Kenue).
- Set up the value of variables: EBOR%ADR(1)%P%R(K)=0.1D-4 (for example)
QBOR%ADR(1)%P%R(K)=1.0D-4 (for example)
At boundary where I want described the values above, I have verified anh found value of LIEBOR%I(K)=KSORT (=4) and LIQBOR%I(K)=KENT (5) .
It seem le program work and done something. Unfortunately, I'm not ready to understand all informations from result file and compare with imposed condition at frontier.
Plus, I have a problem when I use key word DEBITS SOLIDES IMPOSES = 0.0;1.0 in my steering file for sisphe (one of the options as your suggestion above), le program always stop and information is :
DISIPM PROBLEME SUR LA FRONTIERE 2
DONNER UN PROFIL DE DEBIT DE SOLIDE
DANS LE FICHIER DES CONDITIONS AUX LIMITES
I don't know what should I do !
Again, thank you very much for your help.
Best regards
Thong NGUYEN
HCMUT, Viet Nam