!       TITLE:  t_rplot_nodal.mac (VERIFICATION)
!
!---------------------------------------------------------------------------
!               Rev -   01/14/2001      D. Edwards
!---------------------------------------------------------------------------
!
!       Create large 3-D stress components to verify
!       formatting of prnsol output.
!
!       ct1cs:        Cartesian 1 Coordinate System
!       cyl1cs:       Cylindrical 1 Coordinate System
!       sph1cs:       Spherical 1 Coordinate System
!
!
!******************************* TABLE OF CONTENTS **************************
!
!           Model Assembly Preparation
!
!       1.  Parameter List
!       2.  Model Coordinate System
!       3.  Material and Element Type
!       4.  Solid Model
!
!************************** Model Assembly Preparation **********************
!
mod_assemb_prep

/title,Box for testing rplot_nodal

/units,bin
!
!<<<<<<<<<<<<<<<<<<<<<<<< 1.  Parameter List >>>>>>>>>>>>>>>>>>>>>>>>>>
!
t_rplot_nodal.parm
!
!<<<<<<<<<<<<<<<<<<<< 2.  Model Coordinate Systems >>>>>>>>>>>>>>>>>>>>
!
mod_coord_system,0,1,1,1,0,0,0
!
!---------------------- Material Properties --------------------
!
cres=1
UIMP,cres,EX, , ,10.3E6,          ! Material:  cres
UIMP,cres,DENS, , ,.097/386.4,
UIMP,cres,NUXY, , ,0.33,
!
!------------------------- Element Type ------------------------
!
ET,1,plane42
ET,2,solid45
!
!<<<<<<<<<<<<<<<<<<<<<<<<<<< 4.  Solid Model >>>>>>>>>>>>>>>>>>>>>>>>>>
!

csys,ct1cs

k,,0,0
k,,elen,0
k,,elen,elen
k,,0,elen

create_lines
isolate_lines

al,all

init_entity_num

CREATE_VOL-AREA,elen

csys,ct1cs
asel,s,loc,z,0                          ! Bottom surface area


!                  (1)      (2)     (3)     (4)      (5)       (6)     (7)   (8)
! MSH_ALGORITHM, perflsz, perfmsh, mstyp, mshdens, elmszovr, atargov, etpa, etpv,
!
!        elmnl, elmxl, mtrl, rl, ecs, volsel, sprat, radrot, varcang, ell
!          (9)    (10)  (11) (12) (13)  (14)   (15)    (16)     (17)  (18)
!             1 2 3 4 5 6 7 8 9      10
msh_algorithm,1,3, , , , ,1,2,esizel,5*esizel,cres,1

!
! Fix the bottom surface
!

nsel,s,loc,z,0
d,all,all

!
! Displace the top surface
!
nsel,s,loc,z,elen
d,all,uz,-.005

alls

fini

!
! Solve
!

/sol

solv

fini

/post1

! Select node along the z-axis
lsel,s,loc,z,elen/2
lsel,r,loc,x,0
lsel,r,loc,y,0
nsll,s,1

mkasc=1
fcn1='abs'
rplot_nodal,ct1cs,'z','u','x'