1 subroutine npSucPatch(nsfnp,nmap) 2 3ccc 4c output: nsfnp, number of the nodal points which lie on the surface, local for one process 5c nmap, a map from nsfnp to nshg, map the node on the surface to the global nodal numbering for one process 6c 7ccc 8 use pointer_data 9 use wallData 10 include "common.h" 11 include "mpif.h" 12 integer nsfnp 13 integer nmap(nshg) 14 15 16c--- Start to find surface nodes by looking for nonzero wall normals --- c 17 nsfnp=0 18 do ind = 1, nshg 19 wnrmsq = wnorm(ind,1)**2 + wnorm(ind,2)**2 + wnorm(ind,3)**2 20 if(wnrmsq.gt.pt5) then 21 nsfnp=nsfnp+1 22 nmap(nsfnp)=ind 23 endif 24 enddo 25 26 return 27 end 28