1*10167291SKenneth E. Jansen subroutine findTurbWall(iTurbWall) 2*10167291SKenneth E. Jansen 3*10167291SKenneth E. Jansenccc 4*10167291SKenneth E. Jansenc 5*10167291SKenneth E. Jansenccc 6*10167291SKenneth E. Jansen use pointer_data 7*10167291SKenneth E. Jansen include "common.h" 8*10167291SKenneth E. Jansen include "mpif.h" 9*10167291SKenneth E. Jansen integer, allocatable :: ienb(:) 10*10167291SKenneth E. Jansen integer iTurbWall(nshg) 11*10167291SKenneth E. Jansen 12*10167291SKenneth E. Jansenc--- Exclude all the nodes on turbulence wall----c 13*10167291SKenneth E. Jansen iTurbWall(:)=0 14*10167291SKenneth E. Jansen do iblk=1,nelblb ! loop element block 15*10167291SKenneth E. Jansen npro=lcblkb(1,iblk+1)-lcblkb(1,iblk) 16*10167291SKenneth E. Jansen nenbl=lcblkb(6,iblk) 17*10167291SKenneth E. Jansen nshl=lcblkb(9,iblk) 18*10167291SKenneth E. Jansen allocate(ienb(nshl)) 19*10167291SKenneth E. Jansen do i=1,npro ! loop element 20*10167291SKenneth E. Jansen iBCB1=miBCB(iblk)%p(i,1) 21*10167291SKenneth E. Jansen if(.not.btest(iBCB1,4))cycle 22*10167291SKenneth E. Jansen ienb(1:nshl)=mienb(iblk)%p(i,1:nshl) 23*10167291SKenneth E. Jansen do j=1,nenbl ! loop elemental vertex 24*10167291SKenneth E. Jansen nn=ienb(j) 25*10167291SKenneth E. Jansen iTurbWall(nn)=1 26*10167291SKenneth E. Jansen enddo ! end loop elemental vertex 27*10167291SKenneth E. Jansen enddo ! end loop element 28*10167291SKenneth E. Jansen deallocate(ienb) 29*10167291SKenneth E. Jansen enddo ! end loop element block 30*10167291SKenneth E. Jansen 31*10167291SKenneth E. Jansen return 32*10167291SKenneth E. Jansen end 33