!$Id:$

subroutine plot2d(ie,ix,ip,x,xl,nie,ndm,nen,nen1,nume,n1,n2)

! * * F E A P * * A Finite Element Analysis Program

!.... Copyright (c) 1984-2024: Regents of the University of California

! All rights reserved

!-----[--.----+----.----+----.-----------------------------------------]

! Purpose: Two dimensional mesh plot routine

! Inputs:

! ie(nie,*) - Assembly data for material sets

! ix(nen1,*)- Element nodal connections

! ip(*) - Sorted element order

! x(ndm,*) - Nodal coordinates

! ndm - Dimension of x array

! nen - Number of nodes on element

! nen1 - Dimension of ix array

! nume - Number of elements/faces

! n1 - Color for plots

! n2 - Outline indicator

! Scratch:

! xl(ndm,*) - Element nodal coordinates

! Outputs:

! none - Plot output to screen/file

!-----[--.----+----.----+----.-----------------------------------------]

implicit none

include 'eldata.h'

include 'pbody.h'

include 'pdata2.h'

integer :: nie, ndm, nen, nen1, n1,n2, nume, isp

integer :: i, j, ii, n,nn

integer :: ie(nie,*),ix(nen1,*), ip(*)

real (kind=8) :: x(ndm,*),xl(ndm,*)

save

! Loop over elements to draw mesh

do nn = 1,nume

n = ip(nn)

if(n.gt.0) then

if(ix(nen1-1,n).ge.0) then

ma = ix(nen1,n)

! Plot correct material number: ma > 0 active material

if( maplt.eq.0 .or. ma.eq.maplt) then

if(n1.eq.0) then

icolr = ma + 1

elseif(n1.lt.0) then

icolr = mod(ix(nen1-1,n),7) + 1

else

icolr = n1

endif

call pppcol(icolr,0)

do i = 1,nen

ii = abs(ix(i,n))

if(ii.gt.0) then

nel = i

do j = 1,ndm

xl(j,i) = x(j,ii)

end do

else

do j = 1,ndm

xl(j,i) = 0.0d0

end do

endif

end do

! Check for a line element

if(ix(1,n).eq.ix(4,n) .and. ix(2,n).eq.ix(3,n)) nel = 2

if(n2.gt.0) then

isp = -1

else

isp = 1

endif

call plot9(ie(nie-1,ma),ix(1,n),xl,ndm,nel,isp)

endif

endif

endif

end do

end subroutine plot2d