C  This routine saves the deformed nodal data of model in unit 49
C  In .dat file use udump card to call this routine
C Once you have the data written you may run rezone.awk to rezone.
      subroutine impd(n,dd,td,xord,ff,v,a,ndeg,ncrd)
c* * * * * *
c
c     user subroutine for output of displacements.
c
c     n            node number
c     dd           array of incremental displacements of this node
c     td (ndeg,numnp)        array of total displacements of this node
c     xord (ncrd,numnp)      array of coordinates of this node
c     ndeg         number of degrees of freedom per node
c
c* * * * * *
      implicit real*8 (a-h,o-z)                                               dp
      dimension dd(ndeg),td(ndeg),xord(ncrd)
     1 ,ff(ndeg),v(1),a(1)
	dimension txord(3)
	include '../common/blnk'

C TXORD will be the coords + total displacements
	if(ndeg.eq.2) td(3)=0.
	do 1 i=1,ndeg
1	 txord(i) = xord(i) + td(i)
	write(49,2) n,txord
2	format(1x,i5,3e13.5)
      return
      end