problem when using "call getfbr"

[sylvain_serra]

Hello,

I have a strange problem with my boundary conditions.

I use combustion module and I want to impose the boundary conditions in usd3pc.f90.
I use:


CALL GETFBR('100',NLELT,LSTELT)
!==========

write(nfecra,*) 'juste apres nlelt = ', NLELT


In the listing, it gives :
juste apres nlelt = 0


I use CS 2.0.4 and gmsh 2.6.1, and create physicals surface and volume
Physical Surface(100) = {68};

If I try to use the GUI, the checkmesh is ok and I can read the output file for my boundary conditions without problems.

You can see my listpre.checkmesh in attachments

Thank you for any help you might bring me !

Sylvain

[Yvan Fournier]

Hello,

This seems strange, and I suspect a memory overwrite.

The best would be to run the case using the Valgrind memory debugger (setting VALGRIND in the runcase, if you have installed Valgrind onyour machine), so as to see if it detects an overwrite, and where.

You could also re-check or post your user subroutines.

Regards,

Yvan

[sylvain_serra]

Hello Yvan and thank you for your answer.

You could find in attachment my subroutine: usd3pc.f90

I don't find the problem...

[Yvan Fournier]

Hello,

An overwrite might be in an other file. Do you use other user subroutines, and did you test under Valgrind ?

Regards,

Yvan