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