/TITLE, Axisymmetric, radial numbering, node2dr.a
  ! create nodes and elements in 2D section
  ! 8 axial, 4 radial sets of elements
  ! radial 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,nrad,dia/2,0    ! lower-right node
FILL,1,nrad       ! bottom radial nodes 
NGEN,nax,10,1,nrad,1,0,(len/elax) ! axial nodes
/PNUM,node,1      ! plot the nodes
NPLOT
E,1,2,12,11       ! lower-left element
EGEN,elra,1,-1    ! bottom radial elements
EGEN,elax,10,-elra ! remaining elements
EPLOT
FINISH            ! preprocessor section