open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello,

I am using the trunk, revision 7307.

I want to write the spectrum of a point to continue a computation with ARTEMIS. I specify the abscissa/ordinate of ONLY ONE POINT point in the TOMAWAC steering file and get the following error (run is in parallel):

When I add a second point, TOMAWAC works, but then I cannot use the PUNCTUAL RESULTS FILE with ARTEMIS because it requires a file with only one point. At firsts it looks that there is a problem with the code that makes TOMAWAC look for a second point, even if it is not specified in the cas file.

Having a look at lecdon_tomawac.f there is the following comment:
So, now the question is whether it is an error with the code, OR the code has changed and the points have to be specified in another manner.

Best Regards,
Costas

Hello Costas,

I did this change when looking for another bug, when I realised that I had NO spectrum point and I was also stopped because a spectrum point was outside the domain.

In lecdon_tomawac.f I changed:

NPLEO=DIMEN(2,4) into:

IF(TROUVE(2,4).EQ.2) THEN
NPLEO=DIMEN(2,4)
ELSE
NPLEO=0
ENDIF

to secure the fact that if the keyword ABSCISSAE OF SPECTRUM PRINTOUT POINTS is not in the user steering file, no spectrum will be considered (I was afraid that the default size given in the dictionary would pop out instead of 0). Normally the original version NPLEO=DIMEN(2,4) should work and return 0 in this case.

To sum up: the number of spectrum points is deduced from the number of abscissae that you give.

Now we have to understand why you get NPLEO = 2 instead of 1, but maybe you will see after these explanations, otherwise I would need the case to investigate further.

With best regards,

Jean-Michel Hervouet

Hello Jean-Michel,

Thank you for your explanation. I will have another look at it. In the meantime, I attach my cas file:



Regards,
Costas

Hi Costas

did you find solution of your problem?

i think i have something similar in a way.
i ask for 19 points and get only 2 back.
the same cas gave 19 outputs with v7p0 but gives 2 with v7p1.
my search brought me to lecdon_tomawac.f > damocle.f > damoc.f but not further yet.

bests
Homayoon

Hello, Homayoon,
I have modified lecdon_tomawac.f and since then it works fine for me. See if it works for you as well:

Regards,
Costas

thank you Costas for sharing your lecdon_tomawac.f

but it didnt solve it. i will post this as a separate post. someone might have faced this.

bests
Homayoon

Having revisited this issue after testing the latest trunk revision, I noticed that the problem still persists. My modification for lecdon_tomawac was simply to hardcode 1 for NPLEO and problem solved for me. Obviously this is not right. The original lecdon_tomawac still returns NPLEO=2 even when I have only one abscissa defined.

Still, I cannot find an explanation. Can someone have a another look?

Best Regards,
Costas

Hello Costas,

For a reason i don't know, the size of Abcisse and ordinate is set to 2 by default, I still don't know in detail how the dictionnary is managed but putting tail : 1 instead of 2 make your case work.

Hope it will still work for other number ....

Regards

Thierry

Hello Thierry,

Thank you for looking into the problem. I have noticed the correction you did in the trunk, but then you undid it. I will wait for a solid solution before I give it a try in my current simulations.

Best Regards,
Costas

Hello Costas

yesterday i found the problem. There was not the word &DYN in the dictionnary, which tells to the reader to adapt dynamically the dimension of arrays read in the steering file

Best regard

Thierry