Page 1 of 1

thinwall

Posted: Tue Mar 17, 2015 2:57 pm
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?

Re: thinwall

Posted: Wed Mar 18, 2015 2:23 pm
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,