Lines Matching refs:nn
105 nn=idxslpw(i)
106 call slpwBC(mpslpw(nn,1),mpslpw(nn,2),iBCg(nn),
107 & BCg(nn,:), BCtmpg(nn,:),
108 & iBC(nn), BC(nn,:),
109 & wlnorm(nn,:,:) )
129 nn = 0
132 nn = nn + 1
133 if(mod(nn,50).eq.1)
258 nn=ienb(j) ! global node number
259 if(firstvisit(nn)) then
260 firstvisit(nn)=.false.
261 nsurf(nn)=nsurf(nn)+1
262 if(nsurf(nn).eq.1) BCtmp(nn,4:6)=zero
263 if(nsurf(nn).eq.2) BCtmp(nn,8:10)=zero
265 wnrm(nn,:)=wnrm(nn,:)+elnrm
318 do nn = 1, nshg
319 if(nsurf(nn).eq.1)
320 & BCtmp(nn,4:6)=BCtmp(nn,4:6)+wnrm(nn,:)
321 if(nsurf(nn).eq.2)
322 & BCtmp(nn,8:10)=BCtmp(nn,8:10)+wnrm(nn,:)
323 if(nsurf(nn).gt.2)
324 & BCtmp(nn,4:6)=BCtmp(nn,4:6)+wnrm(nn,:)
334 do nn = 1, numnp ! loop over all nodes
336 if(nsurf(nn).eq.0) cycle
341 if(ibits(iBC(nn),3,3).eq.7 ) ikp=1
344 ixset=ibits(iBC(nn),3,1)
345 iyset=ibits(iBC(nn),4,1)
346 izset=ibits(iBC(nn),5,1)
347 iBC(nn)=iBC(nn)-ixset*8-iyset*16-izset*32
351 iBCSAV=iBC(nn)
357 if(nsurf(nn).eq.1) then ! 1 modeled wall
361 iBC(nn)=iBC(nn)+8 ! assume normal is x-dominated first
362 wnrm(nn,:)=BCtmp(nn,4:6)
363 if(abs(wnrm(nn,3)).gt.abs(wnrm(nn,2)))then ! z beats y
364 if(abs(wnrm(nn,3)).gt.abs(wnrm(nn,1)))then ! z beats x
365 iBC(nn)=iBC(nn)+24
368 if(abs(wnrm(nn,2)).ge.abs(wnrm(nn,3)))then ! y beats z
369 if(abs(wnrm(nn,2)).gt.abs(wnrm(nn,1)))then ! y beats x
370 iBC(nn)=iBC(nn)+8
373 BCtmp(nn,7)=zero
374 else if(nsurf(nn).eq.2) then
386 e1=BCtmp(nn,4:6)
391 e2=BCtmp(nn,8:10)
442 BCtmp(nn,4:6)=BCvecs(ivec(1),1:3)
443 BCtmp(nn,8:10)=BCvecs(ivec(2),1:3)
445 BCtmp(nn,4:6)=BCvecs(ivec(1),1:3)
446 BCtmp(nn,8:10)=BCvecs(ivec(3),1:3)
448 BCtmp(nn,4:6)=BCvecs(ivec(2),1:3)
449 BCtmp(nn,8:10)=BCvecs(ivec(3),1:3)
451 iBC(nn)=iBC(nn)+ixset*8+iyset*16+izset*32
452 BCtmp(nn,7)=zero
453 BCtmp(nn,11)=zero
454 wnrm(nn,:)=BCtmp(nn,4:6)+BCtmp(nn,8:10)
455 else if(nsurf(nn).gt.2) then
460 wnrm(nn,:)=BCtmp(nn,4:6)+BCtmp(nn,8:10)
461 iBC(nn)=iBC(nn)+56
462 BCtmp(nn,7)=zero
467 iBC(nn)=iBC(nn)+56 ! set 3-comp
468 if(itwmod.eq.-2) BCtmp(nn,7)=zero ! no slip
469 if(nsurf(nn).eq.1)
470 & wnrm(nn,:)=BCtmp(nn,4:6)
471 if(nsurf(nn).ge.2)
472 & wnrm(nn,:)=BCtmp(nn,4:6)+BCtmp(nn,8:10)
476 wmag=sqrt(wnrm(nn,1)*wnrm(nn,1)
477 & +wnrm(nn,2)*wnrm(nn,2)+wnrm(nn,3)*wnrm(nn,3))
478 wnrm(nn,:)=wnrm(nn,:)/wmag
484 iBC(nn)=iBCSAV+56 ! put it back to a comp3
485 BCtmp(nn,:)=BCtmpSAV(nn,:) ! ditto
775 nn=ienb(can)
777 if(nsurf(nn).ne.0) cycle
779 if(ibits(iBC(nn),3,3).eq.7 .and.
780 & BCtmp(nn,7).eq.zero) cycle
801 do nn = 1, nshg
802 if((otwn(nn).eq.nn).and.(nsurf(nn).ne.0)) then ! if a node on a
813 vec(:)=x(nn,:)-x(can,:)
817 otwn(nn)=can