42       dimension BCtmp(nshg,ndof+7)
51       BCtmp = zero
55             where (nBC(:) .ne. 0) BCtmp(:,i) = BCinp(nBC(:),i)
61       if(any(BCtmp(:,12).ne.0)) then
65             acs(:,1) = cos(BCtmp(:,12)) 
66             acs(:,2) = sin(BCtmp(:,12)) 
75          call genwnm (iBC,  BCtmp,   x,   ilwork, iper, nsurf)
79          call genotwn (x,BCtmp, iBC, nsurf)
89          call genBC1 (BCtmp,  iBC,  BC)
95          BCtmpg=BCtmp
97          call genBC1 (BCtmp,  iBC,  BC)
161       subroutine genwnm (iBC, BCtmp,    x,   ilwork, iper, nsurf)  argument
186       dimension  BCtmp(nshg,ndof+7),      fBC(nshg,ndofBC),
215       BCtmpSAV=BCtmp
262                         if(nsurf(nn).eq.1) BCtmp(nn,4:6)=zero
263                         if(nsurf(nn).eq.2) BCtmp(nn,8:10)=zero
280                   fBC(:,1) = cos(BCtmp(:,12)) 
281                   fBC(:,2) = sin(BCtmp(:,12))
320      &              BCtmp(nn,4:6)=BCtmp(nn,4:6)+wnrm(nn,:)
322      &              BCtmp(nn,8:10)=BCtmp(nn,8:10)+wnrm(nn,:)
324      &              BCtmp(nn,4:6)=BCtmp(nn,4:6)+wnrm(nn,:)
362                   wnrm(nn,:)=BCtmp(nn,4:6)
373                   BCtmp(nn,7)=zero
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)
452                   BCtmp(nn,7)=zero
453                   BCtmp(nn,11)=zero
454                   wnrm(nn,:)=BCtmp(nn,4:6)+BCtmp(nn,8:10)
460                   wnrm(nn,:)=BCtmp(nn,4:6)+BCtmp(nn,8:10)
462                   BCtmp(nn,7)=zero
468                if(itwmod.eq.-2) BCtmp(nn,7)=zero ! no slip
470      &              wnrm(nn,:)=BCtmp(nn,4:6)
472      &              wnrm(nn,:)=BCtmp(nn,4:6)+BCtmp(nn,8:10)
485                BCtmp(nn,:)=BCtmpSAV(nn,:) ! ditto
712       subroutine genotwn(x, BCtmp, iBC, nsurf)  argument
746       real*8 x(numnp,nsd),BCtmp(nshg,ndof+7)  local
780      &                    BCtmp(nn,7).eq.zero) cycle
811      &               BCtmp(can,7).eq.zero) then  ! no slip nodes not allowed