/TITLE, Axisymmetric, axial numbering, node2da.a
  ! create nodes and elements in 2D section
  ! 8 axial, 4 radial sets of elements
  ! axial node numbering

/COLOR,PBAK,off   ! black background
/PREP7            ! preprocessor section
ET,1,PLANE42,,,1  ! axisymmetric option
MP,EX,1,30e6      ! Young's modulus, steel
MP,GXY,1,11.2E6   ! shear modulus

len = 8           ! rod length
dia = 4           ! rod diameter
elax = 8          ! axial number of elements
nax = elax+1      ! axial number of nodes
elra = 4          ! radial number of elements
nrad = elra+1     ! radial number of nodes

N,1,0,0           ! lower-left node
N,nax,0,len       ! upper-left node
FILL,1,nax        ! axial nodes
NGEN,nrad,10,1,nax,1,(dia/2/elra),0 ! axial nodes
/PNUM,node,1      ! plot the nodes
NPLOT
E,1,11,12,2       ! lower-left element
EGEN,elax,1,-1    ! axial elements
EGEN,elra,10,-elax ! remaining elements
EPLOT
FINISH            ! preprocessor section