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