1!$Id:$ 2 subroutine pnumna(ix,nen1,nen,numel, ip) 3 4! * * F E A P * * A Finite Element Analysis Program 5 6!.... Copyright (c) 1984-2017: Regents of the University of California 7! All rights reserved 8 9 implicit none 10 11 include 'pbody.h' 12 13 integer nen1,nen, numel, i,nn 14 integer ix(nen1,*), ip(*) 15 16 save 17 18! Tag active nodes 19 20 do nn = 1,numel 21 if(ix(nen1-1,nn).ge.0 .and. 22 & (maplt.eq.0 .or. ix(nen1,nn).eq.maplt)) then 23 do i = 1,nen 24 if(ix(i,nn).gt.0) ip(ix(i,nn)) = 1 25 end do 26 endif 27 end do 28 29 end 30