open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello,

I get a bug while reading the SHIELDS PARAMETER (AC) in lecdon_sisyphe.f (infos read from the steering file).

Indeed it reads only my first value (over 8 indicated as : 0.0542;0.0304;0.0339;0.0382;0.0481;0.0698;0.1060;0.1851).

What is rather frustrating is that it reads okay others parameters like SEDIMENT DIAMETERS, SETTLING VELOCITIES,etc. which are treated the same way in lecdon_sisyphe.f.

By coincidence (?) Fudaa prepro only accepts a single value in the SHIELDS PARAMETER box (and PARAMETER is a singular noun = no 'S'!!).

I don't understand why I couldn't work with a Shields parameter for each grain classes.

Thank very much in advance.

Best regards.

Alexis

Hello,

There is no obvious reason, I checked that we read a number of values equal to the maximum of :

1) your number of Shields parameters
2) the number of sediment classes

This could trigger a bug if you give 8 numbers and have 9 classes, but if you give a correct number, i.e. the number of classes, it should be OK. You could take subroutine lecdon_sisyphe in your Fortran and print the value of DIMENS(2,6) just before the line filling AC, to understand how many values are seen by the library Damocles. It could be that you have a problem on the 9th number or a problem of line length.

With best regards,

Jean-Michel Hervouet

That is the thing: DIMENS(2,6) = 1 whereas it should be 8 !?

Cheers.


Alexis

Do you know where I can modify the numbers of values seen by Damocles (to change from 1 to 8)?

Thank you in advance.

Alexis

Hi
The number of values is linked to the number of parameters entered in your steering file for the keyword.
You probably have a problem with your steering file
Hope this helps

Hello,

I have ask for 8 classes and set 8 values for SEDIMENT DIAMETERS, SHIELDS PARAMETER and SETTLING VELOCITIES but for an obscure reason lecdon_sisyphe only reads the first value indicated in the steering file for SHIELDS PARAMETER.

Anyway, I will set this Shields parameter in my fortran file instead...

Thank you

Alexis

Hello,

I got it. The keyword Shields parameters has no size declared in the dictionary sisyphev6p1.dico (directory lib). We must add after the line 1276 (i.e. line INDEX = 6) the line :

TAILLE = 10

The value 10 is not very important, if you give seven numbers, Damocles will return 7 numbers, the important thing is that it is not 1 (if I remember well how Damocles is programmed). I'll do the correction through svn but it may take some time before it reaches the download section, you can do it yourself right now and tell me if it works.

With best regards,

Jean-Michel Hervouet