anchor.diff:  
       end subroutine
   
  +
  +    subroutine _prepare_mesh()
  +        implicit none
  +        
  +        integer :: ind,node_max
  +        
  +        node_max = max(nnode,nnoded)
  +        
  +        allocate(xyz(3,node_max),dxyz(6,node_max))
  +        
  +        do ind = 1,nnode
  +            xyz(1:3,ind) = xyztp(1:3,ind)
  +        end do
  +       ! didn't assing damptp 
  +        do ind = 1,nnoded
  +            dxyz(1:3,ind) = dxyztp(1:3,ind)
  +        end do
  +
  +    end subroutine 
   !       MESHFS4 + MESHBD 
   !
   !C *******************************************************************
  
positive.code: 
         NNODE=M
       NNB=NNODE-NNF
       write(11,*)  ' NNB=',NNB,' NNF=',NNF,' NNODE=',NNODE
       write(11,*) 
  !
  ! ===========================================================
  !
  
negative.code: 200     CONTINUE
  
         NNF=M     
         END IF    
  !
  ! -----------------------------------------
  