open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello everyone,
I'd like to know if someone has tried to programm a new bed load formula(other than those offered by SISYPHE)and what tips you can give me. Im trying to program a new bed load formula using qsform.f. If you have already used it,please send me also a copy to have a look.
Regards
Najib

Dear Najib,

please find attached a simple example of qsform.f (implementation of Grass formula).

Do not forget to include the following in your steering file (sisyphe):

BED LOAD = YES
BED-LOAD TRANSPORT FORMULA = 0 / user defined qb = Au^3

I hope it helps,

Pablo

i am triying also to program a simple bed load formula
In the first step i am trying to use your simple example of qsform.f (but no result). the same error "INITIALISATION EN CAS DE COUPLAGE : PART= forrtl: severe (157):Program Exception - access violation
or DEBIMP : PROBLEME SUR LA FRONTIERE 2
DONNER UN PROFIL DE VITESSES
DANS LE FICHIER DES CONDITIONS AUX LIMITES
OU VERIFIER LES HAUTEURS D'EAU

Hello,
the test case KD09/ includes the implementation of the Grass formula. You can run this test case to see if the code has been correctly installed and get inspiration for your own implementation.

Before that, please run your case with one of the already implemented sediment transport formulae and also be sure that the hydrodynamics has been correctly solved.

Cheers,

Pablo

Hello
First I would like to thank you for your answer.
As you told me I tried to run the test case KD09, with Telemac V6P1 but I get a lot of errors, I think I have to install V6P3.
For my simple case (study of bed load transport in a small rectangular channel) when I choose one of the already implemented sediment transport formula ( Van Rijn, Einstein, etc..) I have no problem. But when I try to use my formula or (qsform.f implementation of Grass formula), I have an error message.
Just as an indication, i tried to run qsform.f (implementation of Grass formula with the training case, also i have the same problem. DONNER UN PROFIL DE VITESSES
DANS LE FICHIER DES CONDITIONS AUX LIMITES
OU VERIFIER LES HAUTEURS D'EAU!!
I don’t know where is the problem ?

Hello

I try to follow what you have told me. My project works well with the formulas already established (van Rijn, Meyer, etc ...). But when I try to use my own formula I still have this error message.
INITIALISATION EN CAS DE COUPLAGE : PART= 0
forrtl: severe (157):Program Exception - access violation
Image PC Routine line Source
Out4216_win.exe 004011a1 QSFORM 74 te2fort.f
May you just tell me it’s due to what exactly??
Thank you for your understanding

Hello
First of all thank you for the sample,
I tried to program the formula for bed load transport in steep channels by Tetsuro Tsujimoto, but i face some problems while compiling it. It encounter some errors,Please have a look at the Qsform and the problem i get and give me some hints to correct it.
Thanking you in advance!
Regards
Najib

Hi
You cannot use bief objects directly in formulas like you're doing it.
You should use OV or OS functions for computations or use loops to compute each values
Take care on the length of your lines you shouldn't have more than 72 characters per lines...

Hope this helps

Hi,

Sorry for posting in this old thread. I'm a little bit confused.

It looks like the way the Grass formula is implemented is only one-dimensional (see also test case KD09)

I'm not sure, but it looks like all bedload fluxes are defined as 1d fluxes in Sisyphe. Later the x,y-dependency is computed via the velocity.


This means that the Grass formula should look like


Best regards,

Leo