1 subroutine iBCupdate(iBCpart, ienb, iBCB) 2c 3c--------------------------------------------------------------------- 4c This is the subroutine update iBC for deformable wall case 5c where iBCB(:,1) = 16: Turbulence wall 6c make iBC(:) = iBC(:) + 8192 7c--------------------------------------------------------------------- 8c 9 include "common.h" 10 11 dimension iBCpart(nshg), iBCB(npro,ndiBCB), 12 & ienb(npro,nshl) 13 14 do iel = 1, npro 15 if(btest(iBCB(iel,1),4)) then ! turbulence wall 16 do inode = 1, nenbl 17 iglobal = ienb(iel,inode) 18 iBCpart(iglobal) = 8192 19 enddo 20 endif 21 enddo 22 23 return 24 end 25