thinwall

Questions and remarks about code_saturne usage
Forum rules
Please read the forum usage recommendations before posting.
Post Reply
leguichet
Posts: 34
Joined: Wed Sep 03, 2014 1:00 pm

thinwall

Post by leguichet »

Version 3.0.6

I am using the thinwall to add internal boundary condition.
In cs_user_mesh.c, I select the name of the internal faces, for example, 'wall0'
In cs_user_boundary_condition.f90, I imposed the boundary condition on the thinwall as
call getfbr('wall0', nlelt, lstelt)

do ilelt = 1, nlelt
ifac = lstelt(ilelt)
izfppp(ifac) = izone
itypfb(ifac) = iparug
rcodcl(ifac,iu,3) = 0.1d0
rcodcl(ifac,iv,3) = 0.1d0

write(*,*) 'ifac=', ifac
write(nfecra,*) 'ifac=', ifac
enddo

Here, it appears a strange thing:
The 'ifac=', ifac is printed on the screen correctly, but not printed in "listing" at all.
Is that a bug?
Jacques Fontaine
Posts: 118
Joined: Mon Feb 20, 2012 2:07 pm

Re: thinwall

Post by Jacques Fontaine »

Hello,

It is not really a thinwall problem.

Is it multi-cores calculation ? If yes, maybe there is no face 'wall0' on rank 0 (rank 0 print in listing).
Regards,
Post Reply