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